Genome sequence of the palaeopolyploid soybean

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Glycine max polygalacturonase inhibiting protein 11 (GmPGIP11) functions in the root to suppress Heterodera glycines parasitism

Pathogen-secreted polygalacturonases (PGs) alter plant cell wall structure by cleaving the α-(1 → 4) linkages between D-galacturonic acid residues in homogalacturonan (HG), macerating the cell wall, facilitating infection. Plant PG inhibiting proteins (PGIPs) disengage pathogen PGs, impairing infection. The soybean cyst nematode, Heterodera glycines, obligate root parasite produces secretions, generating a multinucleate nurse cell called a syncytium, a byproduct of the merged cytoplasm of 200-250 root cells, occurring through cell wall maceration. The common cytoplasmic pool, surrounded by an intact plasma membrane, provides a source from which H. glycines derives nourishment but without killing the parasitized cell during a susceptible reaction. The syncytium is also the site of a naturally-occurring defense response that happens in specific G. max genotypes. Transcriptomic analyses of RNA isolated from the syncytium undergoing the process of defense have identified that one of the 11 G. max PGIPs, GmPGIP11, is expressed during defense. Functional transgenic analyses show roots undergoing GmPGIP11 overexpression (OE) experience an increase in its relative transcript abundance (RTA) as compared to the ribosomal protein 21 (GmRPS21) control, leading to a decrease in H. glycines parasitism as compared to the overexpression control. The GmPGIP11 undergoing RNAi experiences a decrease in its RTA as compared to the GmRPS21 control with transgenic roots experiencing an increase in H. glycines parasitism as compared to the RNAi control. Pathogen associated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) components are shown to influence GmPGIP11 expression while numerous agricultural crops are shown to have homologs.

Genomic insights into cytokinin oxidase/dehydrogenase (CKX) gene family, identification, phylogeny and synteny analysis for its possible role in regulating seed number in Pigeonpea (Cajanus cajan)

Cytokinin oxidase/dehydrogenase (CKX) regulates cytokinin levels in plants which are vital for plant growth and development. However, there is a paucity of evidence regarding their role in controlling embryo/seed development in pigeonpea. This study is perhaps the first to provide information on the identification and characterization of CKX genes in pigeonpea. A genome-wide analysis identified 18 CKX genes, each with distinct structure, expression patterns, and possible diverse functions. Domain analysis revealed the presence of the sequences including FAD and CK-Binding domain, and subcellular localization analysis showed that almost 50 % of them reside within the nucleus. They were observed to be located unevenly on chromosome numbers 2, 4, 6, 7, and 11 with a majority of them present on the scaffold. The 8 homologous pairs and various orthologous gene pairs provided further insights into their evolution pattern. Further, SNP/Indels variation in CKX genes and haplotype groups among contrasting genotypes for SNPP (seed number per pod) were analyzed. Spatiotemporal expression analysis revealed the significant expression pattern of CcCKX15, CcCKX17, and CcCKX2 in genotypes carrying low SNPP reiterating their role as negative regulators. These genes can be potential targets to undertake crop improvement programs in pigeonpea.

MutL homolog 1 participates in interference-sensitive meiotic crossover formation in soybean

MutL homolog 1 (MLH1), a member of the MutL homolog family, is required for normal recombination in most organisms. However, its role in soybean (Glycine max) remains unclear to date. Here, we characterized the Glycine max female and male sterility 1 (Gmfms1) mutation that reduces pollen grain viability and increases embryo sac abortion in soybean. Map-based cloning revealed that the causal gene of Gmfms1 is Glycine max MutL homolog 1 (GmMLH1), and CRISPR/Cas9 knockout approach further validated that disruption of GmMLH1 confers the female-male sterility phenotype in soybean. Loss of GmMLH1 function disrupted bivalent formation, leading to univalent mis-segregation during meiosis and ultimately to female-male sterility. The Gmmlh1 mutant showed about a 78.16% decrease in meiotic crossover frequency compared to the wild type. The residual chiasmata followed a Poisson distribution, suggesting that interference-sensitive crossover formation was affected in the Gmmlh1 mutant. Furthermore, GmMLH1 could interact with GmMLH3A and GmMLH3B both in vivo and in vitro. Overall, our work demonstrates that GmMLH1 participates in interference-sensitive crossover formation in soybean, and provides additional information about the conserved functions of MLH1 across plant species.

Water wisteria genome reveals environmental adaptation and heterophylly regulation in amphibious plants

Heterophylly is a phenomenon whereby an individual plant dramatically changes leaf shape in response to the surroundings. Hygrophila difformis (Acanthaceae; water wisteria), has recently emerged as a model plant to study heterophylly because of its striking leaf shape variation in response to various environmental factors. When submerged, H. difformis often develops complex leaves, but on land it develops simple leaves. Leaf complexity is also influenced by other factors, such as light density, humidity, and temperature. Here, we sequenced and assembled the H. difformis chromosome-level genome (scaffold N50: 60.43 Mb, genome size: 871.92 Mb), which revealed 36 099 predicted protein-coding genes distributed over 15 pseudochromosomes. H. difformis diverged from its relatives during the Oligocene climate-change period and expanded gene families related to its amphibious habit. Genes related to environmental stimuli, leaf development, and other pathways were differentially expressed in submerged and terrestrial conditions, possibly modulating morphological and physiological acclimation to changing environments. We also found that auxin plays a role in H. difformis heterophylly. Finally, we discovered candidate genes that respond to different environmental conditions and elucidated the role of LATE MERISTEM IDENTITY 1 (LMI1) in heterophylly. We established H. difformis as a model for studying interconnections between environmental adaptation and morphogenesis.

Genome-wide identification and expression profiling of growth‑regulating factor (GRF) and GRF‑interacting factor (GIF) gene families in chickpea and pigeonpea

The growth-regulating factor (GRF) and GRF-interacting factor (GIF) families encode plant-specific transcription factors and play vital roles in plant development and stress response processes. Although GRF and GIF genes have been identified in various plant species, there have been no reports of the analysis and identification of the GRF and GIF transcription factor families in chickpea (Cicer arietinum) and pigeonpea (Cajanus cajan). The present study identified seven CaGRFs, eleven CcGRFs, four CaGIFs, and four CcGIFs. The identified proteins were grouped into eight and three clades for GRFs and GIFs, respectively based on their phylogenetic relationships. A comprehensive in-silico analysis was performed to determine chromosomal location, sub-cellular localization, and types of regulatory elements present in the putative promoter region. Synteny analysis revealed that GRF and GIF genes showed diploid-polyploid topology in pigeonpea, but not in chickpea. Tissue-specific expression data at the vegetative and reproductive stages of the plant showed that GRFs and GIFs were strongly expressed in tissues like embryos, pods, and seeds, indicating that GRFs and GIFs play vital roles in plant growth and development. This research characterized GRF and GIF families and hints at their primary roles in the chickpea and pigeonpea growth and developmental process. Our findings provide potential gene resources and vital information on GRF and GIF gene families in chickpea and pigeonpea, which will help further understand the regulatory role of these gene families in plant growth and development.

Subfunctionalisation and self-repression of duplicated E1 homologues finetunes soybean flowering and adaptation

Soybean is a photoperiod-sensitive staple crop. Its photoperiodic flowering has major consequences for latitudinal adaptation and grain yield. Here, we identify and characterise a flowering locus named Time of flower 4b (Tof4b), which encodes E1-Like b (E1Lb), a homologue of the key soybean floral repressor E1. Tof4b protein physically associates with the promoters of two FLOWERING LOCUS T (FT) genes to repress their transcription and delay flowering to impart soybean adaptation to high latitudes. Three E1 homologues undergo subfunctionalisation and show differential subcellular localisation. Moreover, they all possess self-repression capability and each suppresses the two homologous counterparts. Subfunctionalisation and the transcriptional regulation of E1 genes collectively finetune flowering time and high-latitude adaptation in soybean. We propose a model for the functional fate of the three E1 genes after the soybean whole-genome duplication events, refine the molecular mechanisms underlying high-latitude adaption, and provide a potential molecular-breeding resource.

Development and implementation of nested single-nucleotide polymorphism (SNP) assays for breeding and genetic research applications

SoySNP50K and SoySNP6K are commonly used for soybean (Glycine max) genotyping. The SoySNP50K assay has been used to genetically analyze the entire USDA Soybean Germplasm Collection, while the SoySNP6K assay, containing a subset of 6000 single-nucleotide polymorphisms (SNPs) from SoySNP50K, has been used for quantitative trait loci mapping of different traits. To meet the needs for genomic selection, selection of parents for crosses, and characterization of breeding populations, especially early selection of ideal offspring from thousands of lines, we developed two assays, SoySNP3K and SoySNP1K, containing 3072 and 1252 SNPs, respectively, based on SoySNP50K and SoySNP6K mark sets. These two assays also contained the trait markers reported or contributed by soybean breeders. The SNPs in the SoySNP3K are a subset from SoySNP6K, while the SNPs in the SoySNP1K are a subset from SoySNP3K. These SNPs were chosen to reduce the SNP number in the large linkage blocks while capturing as much of the haplotype diversity as possible. They are highly polymorphic and of high quality. The mean minor allele frequencies of the SNPs in the southern and northern US elites were 0.25 and 0.27 for SoySNP3K, respectively, and 0.29 and 0.33 for SoySNP1K. The selected SNPs are a valuable source for developing targeted amplicon sequencing assay or beadchip assay in soybean. SoySNP3K and SoySNP1K assays are commercialized by Illumina Inc. and AgriPlex Genomics, respectively. Together with SoySNP50K and SoySNP6K, a series of nested assays with different marker densities will serve as additional low-cost genomic tools for genetic, genomic, and breeding research.

A chromosome-level genome assembly and evolutionary analysis of Coregonus ussuriensis Berg

Coregonus ussuriensis Berg, distributed widely in cold waters above 45° N latitude, is a savored freshwater whitefish that has been included in the list of endangered animals as a consequence of overfishing. Lack of genomic information seriously hampers evolutionary and genetic research on C. ussuriensis warranting the need to assemble a high-quality reference genome to promote its genetic breeding. We assembled and constructed a reference chromosome-level C. ussuriensis genome (sequence length, 2.51 Gb; contig N50 length, 4.27 Mb) using PacBio sequencing and Hi-C assembly technology, 3,109 contigs were assembled into scaffolds, resulting in a genome assembly with 40 chromosomes and a scaffold N50 length of 62.20 Mb. In addition, 43,320 protein-coding genes were annotated. The peak Ks position in the species comparison reflects the whole-genome replication event of C. ussuriensis. This chromosome-level genome provides reference data for further studies on the molecular breeding of C. ussuriensis.

Soybean symbiotic-nodule zonation and cell differentiation are defined by NIN2 signaling and GH3-dependent auxin homeostasis

Symbiotic nodules comprise two classes, indeterminate and determinate, defined by the presence/absence of apical meristem and developmental zonation. Why meristem and zonation are absent from determinate nodules remains unclear. Here, we define cell types in developing soybean nodules, highlighting the undifferentiated infection zones and differentiated nitrogen-fixation zones. Auxin governs infection zone maintenance. GRETCHEN HAGEN 3 (GH3) enzymes deactivate auxin by conjugation and promote cell differentiation. gh3 mutants increased undifferentiated cells and enlarged infection zones. The central symbiosis-transcription factor NIN2a activates GH3.1 to reduce auxin levels and facilitates cell differentiation. High auxin promotes NIN2a protein accumulation and enhances signaling, further deactivating auxin and depleting infection zones. Our findings shed light on the NIN2a-GH3-auxin module that drives soybean nodule cell differentiation. This study challenges our understanding of determinate nodule development and proposes that the regulation of nodule zonation offers valuable insights into broader mechanisms of cell differentiation across plant species.

Genome-Wide Investigation of the CRF Gene Family in Maize and Functional Analysis of ZmCRF9 in Response to Multiple Abiotic Stresses

The cytokinin response factors (CRFs) are pivotal players in regulating plant growth, development, and responses to diverse stresses. Despite their significance, comprehensive information on CRF genes in the primary food crop, maize, remains scarce. In this study, a genome-wide analysis of CRF genes in maize was conducted, resulting in the identification of 12 members. Subsequently, we assessed the chromosomal locations, gene duplication events, evolutionary relationships, conserved motifs, and gene structures of all ZmCRF members. Analysis of ZmCRF promoter regions indicated the presence of cis-regulatory elements associated with plant growth regulation, hormone response, and various abiotic stress responses. The expression patterns of maize CRF genes, presented in heatmaps, exhibited distinctive patterns of tissue specificity and responsiveness to multiple abiotic stresses. qRT-PCR experiments were conducted on six selected genes and confirmed the involvement of ZmCRF genes in the plant's adaptive responses to diverse environmental challenges. In addition, ZmCRF9 was demonstrated to positively regulate cold and salt tolerance. Ultimately, we explored the putative interaction partners of ZmCRF proteins. In summary, this systematic overview and deep investigation of ZmCRF9 provides a solid foundation for further exploration into how these genes contribute to the complex interplay of plant growth, development, and responses to stress.

Genetic Control of Tolerance to Drought Stress in Wild Soybean (Glycine soja) at the Vegetative and the Germination Stages

Drought stress, which is becoming more prevalent due to climate change, is a significant abiotic factor that adversely impacts crop production and yield stability. Cultivated soybean (Glycine max), a versatile crop for humans and animals, exhibits sensitivity to drought, resulting in reduced growth and development under drought conditions. However, few genetic studies have assessed wild soybean's (Glycine soja) response to drought stress. In this work, we conducted a genome-wide association study (GWAS) and analysis of wild soybean accessions to identify loci responsible for drought tolerance at the vegetative (n = 187) and the germination stages (n = 135) using the available resequencing data. The GWAS analysis of the leaf wilting score (LWS) identified eight single-nucleotide polymorphisms (SNPs) on chromosomes 10, 11, and 19. Of these, wild soybeans with both SNPs on chromosomes 10 (adenine) and 11 (thymine) produced lower LWS, indicating that these SNPs have an important role in the genetic effect on LWS for drought tolerance at the vegetative stage. At the germination stage, nine SNPs associated with five phenotypic measurements were identified on chromosomes 6, 9, 10, 13, 16, and 17, and the genomic regions identified at the germination stage were different from those identified for the LWS, supporting our previous finding that there may not be a robust correlation between the genes influencing phenotypes at the germination and vegetative stages. This research will benefit marker-assisted breeding programs aimed at enhancing drought tolerance in soybeans.

Mechanisms of vacuolar phosphate efflux supporting soybean root hair growth in response to phosphate deficiency

Phosphorus is an essential macronutrient for plant growth and development. In response to phosphate (Pi) deficiency, plants rapidly produce a substitutive amount of root hairs; however, the mechanisms underlying Pi supply for root hair growth remain unclear. Here, we observed that soybean (Glycine max) plants maintain a consistent level of Pi within root hairs even under external Pi deficiency. We therefore investigated the role of vacuole-stored Pi, a major Pi reservoir in plant cells, in supporting root hair growth under Pi-deficient conditions. Our findings indicated that two vacuolar Pi efflux (VPE) transporters, GmVPE1 and GmVPE2, remobilize vacuolar stored Pi to sustain cytosolic Pi content in root hair cells. Genetic analysis showed that double mutants of GmVPE1 and GmVPE2 exhibited reduced root hair growth under low Pi conditions. Moreover, GmVPE1 and GmVPE2 were highly expressed in root hairs, with their expression levels significantly upregulated by low Pi treatment. Further analysis revealed that GmRSL2 (ROOT HAIR DEFECTIVE 6-like 2), a transcription factor involved in root hair morphogenesis, directly binds to the promoter regions of GmVPE1 and GmVPE2, and promotes their expressions under low Pi conditions. Additionally, mutants lacking both GmRSL2 and its homolog GmRSL3 exhibited impaired root hair growth under low Pi stress, which was rescued by overexpressing either GmVPE1 or GmVPE2. Taken together, our study has identified a module comprising vacuolar Pi exporters and transcription factors responsible for remobilizing vacuolar Pi to support root hair growth in response to Pi deficiency in soybean.

Divergent evolution of NLR genes in the genus Glycine: impacts of annuals and perennials' life history strategies

Within the family Fabaceae, the genus Glycine is composed of two subgenera annuals (2n=40) and perennials. This life strategy transition may have differentially affected the evolution of various gene families. Its cultivated species G. max has high level of susceptibility to major pathogens including viruses, bacteria and fungi. Understanding nucleotide-binding domain leucine-rich repeat (NLR) genes evolution in soybean is of paramount importance due to their central role in plant immunity and their potential in improving disease resistance in soybean cultivars. In this study, we investigated the significance of this annual-perennial transition on the macroevolution of NLR genes in the genus Glycine. Our results reveal a remarkable distinction between annual species such as Glycine max and Glycine soja, which exhibit an expanded NLRome compared to perennial species (G. cyrtoloba, G. stenophita, G. dolichocarpa, G. falcata, G. syndetika, G. latifolia and G. tomentella). Our evolutionary timescale analysis pinpoints recent accelerated gene duplication events for this expansion, which occurred between 0.1 and 0.5 million years ago, driven predominantly by lineage-specific and terminal duplications. In contrast, perennials initially experienced significant contraction during the diploidisation phase following the Glycine-specific whole-genome duplication event (~10 million years ago). Despite the reduction in the NLRome, perennial lineages exhibit a unique and highly diversified repertoire of NLR genes with limited interspecies synteny. The investigation of gene gain and loss ratios revealed that this diversification resulted from the birth of novel genes following individual speciation events. Among perennials, G. latifolia, a well-known resistance resource, has the highest ratio of these novel genes in the tertiary gene pool. Our study suggests evolutionary mechanisms, including recombination and transposition, as potential drivers for the emergence of these novel genes. This study also provides evidence for the unbalanced expansion of the NLRome in the Dt subgenome compared with the At subgenome in the young allopolyploid G. dolichocarpa. To the best of our knowledge, this is the first study to investigate the effect of annuality and perenniality life transition on the evolution of NLR genes in the genus Glycine to identify its genomics resources for improving the resistance of soybean crop with global importance on the economy and food security.

A reference-grade genome of the xerophyte Ammopiptanthus mongolicus sheds light on its evolution history in legumes and drought-tolerance mechanisms

Plants that grow in extreme environments represent unique sources of stress-resistance genes and mechanisms. Ammopiptanthus mongolicus (Leguminosae) is a xerophytic evergreen broadleaf shrub native to semi-arid and desert regions; however, its drought-tolerance mechanisms remain poorly understood. Here, we report the assembly of a reference-grade genome for A. mongolicus, describe its evolutionary history within the legume family, and examine its drought-tolerance mechanisms. The assembled genome is 843.07 Mb in length, with 98.7% of the sequences successfully anchored to the nine chromosomes of A. mongolicus. The genome is predicted to contain 47 611 protein-coding genes, and 70.71% of the genome is composed of repetitive sequences; these are dominated by transposable elements, particularly long-terminal-repeat retrotransposons. Evolutionary analyses revealed two whole-genome duplication (WGD) events at 130 and 58 million years ago (mya) that are shared by the genus Ammopiptanthus and other legumes, but no species-specific WGDs were found within this genus. Ancestral genome reconstruction revealed that the A. mongolicus genome has undergone fewer rearrangements than other genomes in the legume family, confirming its status as a "relict plant". Transcriptomic analyses demonstrated that genes involved in cuticular wax biosynthesis and transport are highly expressed, both under normal conditions and in response to polyethylene glycol-induced dehydration. Significant induction of genes related to ethylene biosynthesis and signaling was also observed in leaves under dehydration stress, suggesting that enhanced ethylene response and formation of thick waxy cuticles are two major mechanisms of drought tolerance in A. mongolicus. Ectopic expression of AmERF2, an ethylene response factor unique to A. mongolicus, can markedly increase the drought tolerance of transgenic Arabidopsis thaliana plants, demonstrating the potential for application of A. mongolicus genes in crop improvement.

Oceanobacillus picturae alleviates cadmium stress and promotes growth in soybean seedlings

Cadmium (Cd) is a heavy metal that significantly impacts human health and the environment. Microorganisms play a crucial role in reducing heavy metal stress in plants; however, the mechanisms by which microorganisms enhance plant tolerance to Cd stress and the interplay between plants and microorganisms under such stress remain unclear. In this study, Oceanobacillus picturae (O. picturae) was isolated for interaction with soybean seedlings under Cd stress. Results indicated that Cd treatment alone markedly inhibited soybean seedling growth. Conversely, inoculation with O. picturae significantly improved growth indices such as plant height, root length, and fresh weight, while also promoting recovery in soil physiological indicators and pH. Metabolomic and transcriptomic analyses identified 157 genes related to aspartic acid, cysteine, and flavonoid biosynthesis pathways. Sixty-three microbial species were significantly associated with metabolites in these pathways, including pathogenic, adversity-resistant, and bioconductive bacteria. This research experimentally demonstrates, for the first time, the growth-promoting effect of the O. picturae strain on soybean seedlings under non-stress conditions. It also highlights its role in enhancing root growth and reducing Cd accumulation in the roots under Cd stress. Additionally, through the utilization of untargeted metabolomics, metagenomics, and transcriptomics for a multi-omics analysis, we investigated the impact of O. picturae on the soil microbiome and its correlation with differential gene expression in plants. This innovative approach unveils the molecular mechanisms underlying O. picturae's promotion of root growth and adaptation to Cd stress.

A spatially resolved multiomic single-cell atlas of soybean development

Cis-regulatory elements (CREs) precisely control spatiotemporal gene expression in cells. Using a spatially resolved single-cell atlas of gene expression with chromatin accessibility across ten soybean tissues, we identified 103 distinct cell types and 303,199 accessible chromatin regions (ACRs). Nearly 40% of the ACRs showed cell-type-specific patterns and were enriched for transcription factor (TF) motifs defining diverse cell identities. We identified de novo enriched TF motifs and explored conservation of gene regulatory networks underpinning legume symbiotic nitrogen fixation. With comprehensive developmental trajectories for endosperm and embryo, we uncovered the functional transition of the three sub-cell types of endosperm, identified 13 sucrose transporters sharing the DOF11 motif that were co-up-regulated in late peripheral endosperm and identified key embryo cell-type specification regulators during embryogenesis, including a homeobox TF that promotes cotyledon parenchyma identity. This resource provides a valuable foundation for analyzing gene regulatory programs in soybean cell types across tissues and life stages.

The Soybean Rpp3 Gene Encodes a TIR-NBS-LRR Protein that Confers Resistance to Phakopsora pachyrhizi

Soybean rust is an economically significant disease caused by the fungus Phakopsora pachyrhizi that negatively impacts soybean (Glycine max [L.] Merr.) production throughout the world. Susceptible plants infected by P. pachyrhizi develop tan-colored lesions on the leaf surface that give rise to funnel-shaped uredinia as the disease progresses. While most soybean germplasm is susceptible, seven genetic loci (Rpp1 to Rpp7) that provide race-specific resistance to P. pachyrhizi (Rpp) have been identified. Rpp3 was first discovered and characterized in the soybean accession PI 462312 (Ankur), and it was also determined to be one of two Rpp genes present in PI 506764 (Hyuuga). Genetic crosses with PI 506764 were later used to fine-map the Rpp3 locus to a 371-kb region on chromosome 6. The corresponding region in the susceptible Williams 82 (Wm82) reference genome contains several homologous nucleotide binding site-leucine rich repeat (NBS-LRR) genes. To identify Rpp3, we designed oligonucleotide primers to amplify Rpp3 candidate (Rpp3C) NBS-LRR genes at this locus from PI 462312, PI 506764, and Wm82 using polymerase chain reaction (PCR). Five Rpp3C genes were identified in both Rpp3-resistant soybean lines, and co-silencing these genes compromised resistance to P. pachyrhizi. Gene expression analysis and sequence comparisons of the Rpp3C genes in PI 462312 and PI 506764 suggest that a single candidate gene, Rpp3C3, is responsible for Rpp3-mediated resistance. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2024.

The Kirkhouse Trust: Successes and Challenges in Twenty Years of Supporting Independent, Contemporary Grain Legume Breeding Projects in India and African Countries

This manuscript reviews two decades of projects funded by the Kirkhouse Trust (KT), a charity registered in the UK. KT was established to improve the productivity of legume crops important in African countries and in India. KT's requirements for support are: (1) the research must be conducted by national scientists in their home institution, either a publicly funded agricultural research institute or a university; (2) the projects need to include a molecular biology component, which to date has mostly comprised the use of molecular markers for the selection of one or more target traits in a crop improvement programme; (3) the projects funded are included in consortia, to foster the creation of scientific communities and the sharing of knowledge and breeding resources. This account relates to the key achievements and challenges, reflects on the lessons learned and outlines future research priorities.

DEKR-SPrior: An Efficient Bottom-Up Keypoint Detection Model for Accurate Pod Phenotyping in Soybean

The pod and seed counts are important yield-related traits in soybean. High-precision soybean breeders face the major challenge of accurately phenotyping the number of pods and seeds in a high-throughput manner. Recent advances in artificial intelligence, especially deep learning (DL) models, have provided new avenues for high-throughput phenotyping of crop traits with increased precision. However, the available DL models are less effective for phenotyping pods that are densely packed and overlap in in situ soybean plants; thus, accurate phenotyping of the number of pods and seeds in soybean plant is an important challenge. To address this challenge, the present study proposed a bottom-up model, DEKR-SPrior (disentangled keypoint regression with structural prior), for in situ soybean pod phenotyping, which considers soybean pods and seeds analogous to human people and joints, respectively. In particular, we designed a novel structural prior (SPrior) module that utilizes cosine similarity to improve feature discrimination, which is important for differentiating closely located seeds from highly similar seeds. To further enhance the accuracy of pod location, we cropped full-sized images into smaller and high-resolution subimages for analysis. The results on our image datasets revealed that DEKR-SPrior outperformed multiple bottom-up models, viz., Lightweight-OpenPose, OpenPose, HigherHRNet, and DEKR, reducing the mean absolute error from 25.81 (in the original DEKR) to 21.11 (in the DEKR-SPrior) in pod phenotyping. This paper demonstrated the great potential of DEKR-SPrior for plant phenotyping, and we hope that DEKR-SPrior will help future plant phenotyping.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Comparative analysis of lipid and flavonoid biosynthesis between Pongamia and soybean seeds: genomic, transcriptional, and metabolic perspectives

BACKGROUND

Soybean (Glycine max) is a vital oil-producing crop. Augmenting oleic acid (OA) levels in soybean oil enhances its oxidative stability and health benefits, representing a key objective in soybean breeding. Pongamia (Pongamia pinnata), known for its abundant oil, OA, and flavonoid in the seeds, holds promise as a biofuel and medicinal plant. A comparative analysis of the lipid and flavonoid biosynthesis pathways in Pongamia and soybean seeds would facilitate the assessment of the potential value of Pongamia seeds and advance the genetic improvements of seed traits in both species.

RESULTS

The study employed multi-omics analysis to systematically compare differences in metabolite accumulation and associated biosynthetic genes between Pongamia seeds and soybean seeds at the transcriptional, metabolic, and genomic levels. The results revealed that OA is the predominant free fatty acid in Pongamia seeds, being 8.3 times more abundant than in soybean seeds. Lipidomics unveiled a notably higher accumulation of triacylglycerols (TAGs) in Pongamia seeds compared to soybean seeds, with 23 TAG species containing OA. Subsequently, we identified orthologous groups (OGs) involved in lipid biosynthesis across 25 gene families in the genomes of Pongamia and soybean, and compared the expression levels of these OGs in the seeds of the two species. Among the OGs with expression levels in Pongamia seeds more than twice as high as in soybean seeds, we identified one fatty acyl-ACP thioesterase A (FATA) and two stearoyl-ACP desaturases (SADs), responsible for OA biosynthesis, along with two phospholipid:diacylglycerol acyltransferases (PDATs) and three acyl-CoA:diacylglycerol acyltransferases (DGATs), responsible for TAG biosynthesis. Furthermore, we observed a significantly higher content of the flavonoid formononetin in Pongamia seeds compared to soybean seeds, by over 2000-fold. This difference may be attributed to the tandem duplication expansions of 2,7,4'-trihydroxyisoflavanone 4'-O-methyltransferases (HI4'OMTs) in the Pongamia genome, which are responsible for the final step of formononetin biosynthesis, combined with their high expression levels in Pongamia seeds.

CONCLUSIONS

This study extends beyond observations made in single-species research by offering novel insights into the molecular basis of differences in lipid and flavonoid biosynthetic pathways between Pongamia and soybean, from a cross-species comparative perspective.

WRKY Transcription Factor Responses and Tolerance to Abiotic Stresses in Plants

Plants are subjected to abiotic stresses throughout their developmental period. Abiotic stresses include drought, salt, heat, cold, heavy metals, nutritional elements, and oxidative stresses. Improving plant responses to various environmental stresses is critical for plant survival and perpetuation. WRKY transcription factors have special structures (WRKY structural domains), which enable the WRKY transcription factors to have different transcriptional regulatory functions. WRKY transcription factors can not only regulate abiotic stress responses and plant growth and development by regulating phytohormone signalling pathways but also promote or suppress the expression of downstream genes by binding to the W-box [TGACCA/TGACCT] in the promoters of their target genes. In addition, WRKY transcription factors not only interact with other families of transcription factors to regulate plant defence responses to abiotic stresses but also self-regulate by recognising and binding to W-boxes in their own target genes to regulate their defence responses to abiotic stresses. However, in recent years, research reviews on the regulatory roles of WRKY transcription factors in higher plants have been scarce and shallow. In this review, we focus on the structure and classification of WRKY transcription factors, as well as the identification of their downstream target genes and molecular mechanisms involved in the response to abiotic stresses, which can improve the tolerance ability of plants under abiotic stress, and we also look forward to their future research directions, with a view of providing theoretical support for the genetic improvement of crop abiotic stress tolerance.

Genome-Wide Analysis of Serine Carboxypeptidase-like Genes in Soybean and Their Roles in Stress Resistance

The serine carboxypeptidase-like (SCPL) gene family plays a crucial role in the regulation of plant growth, development, and stress response through activities such as acyltransferases in plant secondary metabolism pathways. Although SCPL genes have been identified in various plant species, their specific functions and characteristics in soybean (Glycine max) have not yet been studied. We identified and characterized 73 SCPL genes, grouped into three subgroups based on gene structure and phylogenetic relationships. These genes are distributed unevenly across 20 soybean chromosomes and show varied codon usage patterns influenced by both mutation and selection pressures. Gene ontology (GO) enrichment suggests these genes are involved in plant cell wall regulation and stress responses. Expression analysis in various tissues and under stress conditions, including the presence of numerous stress-related cis-acting elements, indicated that these genes have varied expression patterns. This suggests that they play specialized roles such as modulating plant defense mechanisms against nematode infections, enhancing tolerance to drought and high salinity, and responding to cold stress, thereby helping soybean adapt to environmental stresses. Moreover, the expression of specific GmSCPLs was significantly affected following exposure to nematode infection, drought, high salt (NaCl), and cold stresses. Our findings underscore the potential of SCPL genes in enhancing stress resistance in soybean, providing a valuable resource for future genetic improvement and breeding strategies.

The MADS-box transcription factor GmFULc promotes GmZTL4 gene transcription to modulate maturity in soybean

Flowering time and maturity are crucial agronomic traits that affect the regional adaptability of soybean plants. The development of soybean cultivars with early maturity adapted to longer days and colder climates of high latitudes is very important for ensuring normal ripening before frost begins. FUL belongs to the MADS-box transcription factor family and has several duplicated members in soybeans. In this study, we observed that overexpression of GmFULc in the Dongnong 50 cultivar promoted soybean maturity, while GmFULc knockout mutants exhibited late maturity. Chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) revealed that GmFULc could bind to the CArG, bHLH and homeobox motifs. Further investigation revealed that GmFULc could directly bind to the CArG motif in the promoters of the GmZTL3 and GmZTL4 genes. Overexpression of GmZTL4 promoted soybean maturity, whereas the ztl4 mutants exhibited delayed maturity. Moreover, we found that the cis element box 4 motif of the GmZTL4 promoter, a motif of light response elements, played an important role in controlling the growth period. Deletion of this motif shortened the growth period by increasing the expression levels of GmZTL4. Functional investigations revealed that short-day treatment promoted the binding of GmFULc to the promoter of GmZTL4 and inhibited the expression of E1 and E1Lb, ultimately resulting in the promotion of flowering and early maturation. Taken together, these findings suggest a novel photoperiod regulatory pathway in which GmFULc directly activates GmZTL4 to promote earlier maturity in soybean.

Soybean hypocotyl elongation is regulated by a MYB33-SWEET11/21-GA2ox8c module involving long-distance sucrose transport

The length of hypocotyl affects the height of soybean and lodging resistance, thus determining the final grain yield. However, research on soybean hypocotyl length is scarce, and the regulatory mechanisms are not fully understood. Here, we identified a module controlling the transport of sucrose, where sucrose acts as a messenger moved from cotyledon to hypocotyl, regulating hypocotyl elongation. This module comprises four key genes, namely MYB33, SWEET11, SWEET21 and GA2ox8c in soybean. In cotyledon, MYB33 is responsive to sucrose and promotes the expression of SWEET11 and SWEET21, thereby facilitating sucrose transport from the cotyledon to the hypocotyl. Subsequently, sucrose transported from the cotyledon up-regulates the expression of GA2ox8c in the hypocotyl, which ultimately affects the length of the hypocotyl. During the domestication and improvement of soybean, an allele of MYB33 with enhanced abilities to promote SWEET11 and SWEET21 has gradually become enriched in landraces and cultivated varieties, SWEET11 and SWEET21 exhibit high conservation and have undergone a strong purified selection and GA2ox8c is under a strong artificial selection. Our findings identify a new molecular pathway in controlling soybean hypocotyl elongation and provide new insights into the molecular mechanism of sugar transport in soybean.

GmDFB1, an ARM-repeat superfamily protein, regulates floral organ identity through repressing siRNA- and miRNA-mediated gene silencing in soybean

The development of flowers in soybean (Glycine max) is essential for determining the yield potential of the plant. Gene silencing pathways are involved in modulating flower development, but their full elucidation is still incomplete. Here, we conducted a forward genetic screen and identified an abnormal flower mutant, deformed floral bud1-1 (Gmdfb1-1), in soybean. We mapped and identified the causal gene, which encodes a member of the armadillo (ARM)-repeat superfamily. Using small RNA sequencing (sRNA-seq), we found an abnormal accumulation of small interfering RNAs (siRNAs) and microRNA (miRNAs) in the Gmdfb1 mutants. We further demonstrated that GmDFB1 interacts with the RNA exosome cofactor SUPER KILLER7 (GmSKI7). Additionally, GmDFB1 interacts with the PIWI domain of ARGONAUTE 1 (GmAGO1) to inhibit the cleavage efficiency on the target genes of sRNAs. The enhanced gene silencing mediated by siRNA and miRNA in the Gmdfb1 mutants leads to the downregulation of their target genes associated with flower development. This study revealed the crucial role of GmDFB1 in regulating floral organ identity in soybean probably by participating in two distinct gene silencing pathways.

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

A comparative plastome approach enhances the assessment of genetic variation in the Melilotus genus

BACKGROUND

Melilotus, a member of the Fabaceae family, is a pivotal forage crop that is extensively cultivated in livestock regions globally due to its notable productivity and ability to withstand abiotic stress. However, the genetic attributes of the chloroplast genome and the evolutionary connections among different Melilotus species remain unresolved.

RESULTS

In this study, we compiled the chloroplast genomes of 18 Melilotus species and performed a comprehensive comparative analysis. Through the examination of protein-coding genes, we successfully established a robust phylogenetic tree for these species. This conclusion is further supported by the phylogeny derived from single-nucleotide polymorphisms (SNPs) across the entire chloroplast genome. Notably, our findings revealed that M. infestus, M. siculus, M. sulcatus, and M. speciosus formed a distinct subgroup within the phylogenetic tree. Additionally, the chloroplast genomes of these four species exhibit two shared inversions. Moreover, inverted repeats were observed to have reemerged in six species within the IRLC. The distribution patterns of single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within protein-coding genes indicated that ycf1 and ycf2 accumulated nonconservative alterations during evolutionary development. Furthermore, an examination of the evolutionary rate of protein-coding genes revealed that rps18, rps7, and rpl16 underwent positive selection specifically in Melilotus.

CONCLUSIONS

We present a comparative analysis of the complete chloroplast genomes of Melilotus species. This study represents the most thorough and detailed exploration of the evolution and variability within the genus Melilotus to date. Our study provides valuable chloroplast genomic information for improving phylogenetic reconstructions and making biogeographic inferences about Melilotus and other Papilionoideae species.

High-integrity Pueraria montana var. lobata genome and population analysis revealed the genetic diversity of Pueraria genus

Pueraria montana var. lobata (P. lobata) is a traditional medicinal plant belonging to the Pueraria genus of Fabaceae family. Pueraria montana var. thomsonii (P. thomsonii) and Pueraria montana var. montana (P. montana) are its related species. However, evolutionary history of the Pueraria genus is still largely unknown. Here, a high-integrity, chromosome-level genome of P. lobata and an improved genome of P. thomsonii were reported. It found evidence for an ancient whole-genome triplication and a recent whole-genome duplication shared with Fabaceae in three Pueraria species. Population genomics of 121 Pueraria accessions demonstrated that P. lobata populations had substantially higher genetic diversity, and P. thomsonii was probably derived from P. lobata by domestication as a subspecies. Selection sweep analysis identified candidate genes in P. thomsonii populations associated with the synthesis of auxin and gibberellin, which potentially play a role in the expansion and starch accumulation of tubers in P. thomsonii. Overall, the findings provide new insights into the evolutionary and domestication history of the Pueraria genome and offer a valuable genomic resource for the genetic improvement of these species.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Deciphering PDH1's role in mung bean domestication: a genomic perspective on pod dehiscence

Mung bean (Vigna radiata) stands as a crucial legume crop in Asia, contributing to food security. However, our understanding of the underlying genetic foundation governing domesticated agronomic traits, especially those linked to pod architecture, remains largely unexplored. In this study, we delved into the genomic divergence between wild and domesticated mung bean varieties, leveraging germplasm obtained from diverse sources. Our findings unveiled pronounced variation in promoter regions (35%) between the two mung bean subpopulations, suggesting substantial changes in gene expression patterns during domestication. Leveraging transcriptome analysis using distinct reproductive stage pods and subpopulations, we identified candidate genes responsible for pod and seed architecture development, along with Genome-Wide Association Studies (GWAS) and Quantitative Trait Locus (QTL) analysis. Notably, our research conclusively confirmed PDH1 as a parallel domesticated gene governing pod dehiscence in legumes. This study imparts valuable insights into the genetic underpinnings of domesticated agronomic traits in mung bean, and simultaneously highlighting the parallel domestication of pivotal traits within the realm of legume crops.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Long noncoding RNAs underlie multiple domestication traits and leafhopper resistance in soybean

The origin and functionality of long noncoding RNA (lncRNA) remain poorly understood. Here, we show that multiple quantitative trait loci modulating distinct domestication traits in soybeans are pleiotropic effects of a locus composed of two tandem lncRNA genes. These lncRNA genes, each containing two inverted repeats, originating from coding sequences of the MYB genes, function in wild soybeans by generating clusters of small RNA (sRNA) species that inhibit the expression of their MYB gene relatives through post-transcriptional regulation. By contrast, the expression of lncRNA genes in cultivated soybeans is severely repressed, and, consequently, the corresponding MYB genes are highly expressed, shaping multiple distinct domestication traits as well as leafhopper resistance. The inverted repeats were formed before the divergence of the Glycine genus from the Phaseolus-Vigna lineage and exhibit strong structure-function constraints. This study exemplifies a type of target for selection during plant domestication and identifies mechanisms of lncRNA formation and action.

CircRNA: An emerging star in plant research: A review

CircRNAs are a class of covalently closed non-coding RNA formed by linking the 5' terminus and the 3' terminus after reverse splicing. CircRNAs are widely found in eukaryotes, and they are highly conserved, with spatio-temporal expression specificity and stability. CircRNAs can act as miRNA sponges to regulate the expression of downstream target genes, regulating the transcription of parental genes and some can even be translated into peptides or proteins. Research on circRNAs in plants is still in its infancy compared to that in animals. With the deepening of research, the results of a variety of plant circRNAs suggest that they play an important role in growth and development, and tolerance towards abiotic stresses such as salt, drought, low temperature, high temperature and other adverse environments. In this review paper, we elaborated the molecular characteristics, mechanism of action, function and bioinformatics databases of plant circRNAs, combined with the progress of circRNA research in animals, discussed the potential mechanism of action of plant circRNAs, and proposed the unsolved problems and prospects for future application of plant circRNAs.

Genome resources for three modern cotton lines guide future breeding efforts

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars ('UGA230', 'UA48' and 'CSX8308') and updating the 'TM-1' cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from 'Pima' cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

Physical Seed Dormancy in Legumes: Molecular Advances and Perspectives

Physical dormancy of seeds is a form of dormancy due to the presence of an impermeable seed coat layer, and it represents a feature for plants to adapt to environmental changes over an extended period of phylogenetic evolution. However, in agricultural practice, physical dormancy is problematic. because it prevents timely and uniform seed germination. Therefore, physical dormancy is an important agronomical trait to target in breeding and domestication, especially for many leguminous crops. Compared to the well-characterized physiological dormancy, research progress on physical dormancy at the molecular level has been limited until recent years, due to the lack of suitable research materials. This review focuses on the structure of seed coat, factors affecting physical dormancy, genes controlling physical dormancy, and plants suitable for studying physical dormancy at the molecular level. Our goal is to provide a plethora of information for further molecular research on physical dormancy.

Comparative Genome-Wide Identification of the Fatty Acid Desaturase Gene Family in Tea and Oil Tea

Camellia oil is valuable as an edible oil and serves as a base material for a range of high-value products. Camellia plants of significant economic importance, such as Camellia sinensis and Camellia oleifera, have been classified into sect. Thea and sect. Oleifera, respectively. Fatty acid desaturases play a crucial role in catalyzing the formation of double bonds at specific positions of fatty acid chains, leading to the production of unsaturated fatty acids and contributing to lipid synthesis. Comparative genomics results have revealed that expanded gene families in oil tea are enriched in functions related to lipid, fatty acid, and seed processes. To explore the function of the FAD gene family, a total of 82 FAD genes were identified in tea and oil tea. Transcriptome data showed the differential expression of the FAD gene family in mature seeds of tea tree and oil tea tree. Furthermore, the structural analysis and clustering of FAD proteins provided insights for the further exploration of the function of the FAD gene family and its role in lipid synthesis. Overall, these findings shed light on the role of the FAD gene family in Camellia plants and their involvement in lipid metabolism, as well as provide a reference for understanding their function in oil synthesis.

Bowman-Birk Inhibitor Mutants of Soybean Generated by CRISPR-Cas9 Reveal Drastic Reductions in Trypsin and Chymotrypsin Inhibitor Activities

Despite the high quality of soybean protein, raw soybeans and soybean meal cannot be directly included in animal feed mixtures due to the presence of Kunitz (KTi) and Bowman-Birk protease inhibitors (BBis), which reduces animal productivity. Heat treatment can substantially inactivate trypsin and chymotrypsin inhibitors (BBis), but such treatment is energy-intensive, adds expense, and negatively impacts the quality of seed proteins. As an alternative approach, we have employed CRISPR/Cas9 gene editing to create mutations in BBi genes to drastically lower the protease inhibitor content in soybean seed. Agrobacterium-mediated transformation was used to generate several stable transgenic soybean events. These independent CRISPR/Cas9 events were examined in comparison to wild-type plants using Sanger sequencing, proteomic analysis, trypsin/chymotrypsin inhibitor activity assays, and qRT-PCR. Collectively, our results demonstrate the creation of an allelic series of loss-of-function mutations affecting the major BBi gene in soybean. Mutations in two of the highly expressed seed-specific BBi genes lead to substantial reductions in both trypsin and chymotrypsin inhibitor activities.

A DNA demethylase reduces seed size by decreasing the DNA methylation of AT-rich transposable elements in soybean

Understanding how to increase soybean yield is crucial for global food security. The genetic and epigenetic factors influencing seed size, a major crop yield determinant, are not fully understood. We explore the role of DNA demethylase GmDMEa in soybean seed size. Our research indicates that GmDMEa negatively correlates with soybean seed size. Using CRISPR-Cas9, we edited GmDMEa in the Dongnong soybean cultivar, known for small seeds. Modified plants had larger seeds and greater yields without altering plant architecture or seed nutrition. GmDMEa preferentially demethylates AT-rich transposable elements, thus activating genes and transcription factors associated with the abscisic acid pathway, which typically decreases seed size. Chromosomal substitution lines confirm that these modifications are inheritable, suggesting a stable epigenetic method to boost seed size in future breeding. Our findings provide insights into epigenetic seed size control and suggest a strategy for improving crop yields through the epigenetic regulation of crucial genes. This work implies that targeted epigenetic modification has practical agricultural applications, potentially enhancing food production without compromising crop quality.

SOS1 tonoplast neo-localization and the RGG protein SALTY are important in the extreme salinity tolerance of Salicornia bigelovii

The identification of genes involved in salinity tolerance has primarily focused on model plants and crops. However, plants naturally adapted to highly saline environments offer valuable insights into tolerance to extreme salinity. Salicornia plants grow in coastal salt marshes, stimulated by NaCl. To understand this tolerance, we generated genome sequences of two Salicornia species and analyzed the transcriptomic and proteomic responses of Salicornia bigelovii to NaCl. Subcellular membrane proteomes reveal that SbiSOS1, a homolog of the well-known SALT-OVERLY-SENSITIVE 1 (SOS1) protein, appears to localize to the tonoplast, consistent with subcellular localization assays in tobacco. This neo-localized protein can pump Na+ into the vacuole, preventing toxicity in the cytosol. We further identify 11 proteins of interest, of which SbiSALTY, substantially improves yeast growth on saline media. Structural characterization using NMR identified it as an intrinsically disordered protein, localizing to the endoplasmic reticulum in planta, where it can interact with ribosomes and RNA, stabilizing or protecting them during salt stress.

Light control of three-dimensional chromatin organization in soybean

Higher-order chromatin structure is critical for regulation of gene expression. In plants, light profoundly affects the morphogenesis of emerging seedlings as well as global gene expression to ensure optimal adaptation to environmental conditions. However, the changes and functional significance of chromatin organization in response to light during seedling development are not well documented. We constructed Hi-C contact maps for the cotyledon, apical hook and hypocotyl of soybean subjected to dark and light conditions. The resulting high-resolution Hi-C contact maps identified chromosome territories, A/B compartments, A/B sub-compartments, TADs (Topologically Associated Domains) and chromatin loops in each organ. We observed increased chromatin compaction under light and we found that domains that switched from B sub-compartments in darkness to A sub-compartments under light contained genes that were activated during photomorphogenesis. At the local scale, we identified a group of TADs constructed by gene clusters consisting of different numbers of Small Auxin-Upregulated RNAs (SAURs), which exhibited strict co-expression in the hook and hypocotyl in response to light stimulation. In the hypocotyl, RNA polymerase II (RNAPII) regulated the transcription of a SAURs cluster under light via TAD condensation. Our results suggest that the 3D genome is involved in the regulation of light-related gene expression in a tissue-specific manner.

Circadian and photoperiodic regulation of the vegetative to reproductive transition in plants

As sessile organisms, plants must respond constantly to ever-changing environments to complete their life cycle; this includes the transition from vegetative growth to reproductive development. This process is mediated by photoperiodic response to sensing the length of night or day through circadian regulation of light-signaling molecules, such as phytochromes, to measure the length of night to initiate flowering. Flowering time is the most important trait to optimize crop performance in adaptive regions. In this review, we focus on interplays between circadian and light signaling pathways that allow plants to optimize timing for flowering and seed production in Arabidopsis, rice, soybean, and cotton. Many crops are polyploids and domesticated under natural selection and breeding. In response to adaptation and polyploidization, circadian and flowering pathway genes are epigenetically reprogrammed. Understanding the genetic and epigenetic bases for photoperiodic flowering will help improve crop yield and resilience in response to climate change.

Integrating Computational and Experimental Methods to Identify Novel Sweet Peptides from Egg and Soy Proteins

Sweetness in food delivers a delightful sensory experience, underscoring the crucial role of sweeteners in the food industry. However, the widespread use of sweeteners has sparked health concerns. This underscores the importance of developing and screening natural, health-conscious sweeteners. Our study represents a groundbreaking venture into the discovery of such sweeteners derived from egg and soy proteins. Employing virtual hydrolysis as a novel technique, our research entailed a comprehensive screening process that evaluated biological activity, solubility, and toxicity of the derived compounds. We harnessed cutting-edge machine learning methodologies, specifically the latest graph neural network models, for predicting the sweetness of molecules. Subsequent refinements were made through molecular docking screenings and molecular dynamics simulations. This meticulous research approach culminated in the identification of three promising sweet peptides: DCY(Asp-Cys-Tyr), GGR(Gly-Gly-Arg), and IGR(Ile-Gly-Arg). Their binding affinity with T1R2/T1R3 was lower than -15 kcal/mol. Using an electronic tongue, we verified the taste profiles of these peptides, with IGR emerging as the most favorable in terms of taste with a sweetness value of 19.29 and bitterness value of 1.71. This study not only reveals the potential of these natural peptides as healthier alternatives to traditional sweeteners in food applications but also demonstrates the successful synergy of computational predictions and experimental validations in the realm of flavor science.

The genomes of seven economic Caesalpinioideae trees provide insights into polyploidization history and secondary metabolite biosynthesis

The Caesalpinioideae subfamily contains many well-known trees that are important for economic sustainability and human health, but a lack of genomic resources has hindered their breeding and utilization. Here, we present chromosome-level reference genomes for the two food and industrial trees Gleditsia sinensis (921 Mb) and Biancaea sappan (872 Mb), the three shade and ornamental trees Albizia julibrissin (705 Mb), Delonix regia (580 Mb), and Acacia confusa (566 Mb), and the two pioneer and hedgerow trees Leucaena leucocephala (1338 Mb) and Mimosa bimucronata (641 Mb). Phylogenetic inference shows that the mimosoid clade has a much higher evolutionary rate than the other clades of Caesalpinioideae. Macrosynteny comparison suggests that the fusion and breakage of an unstable chromosome are responsible for the difference in basic chromosome number (13 or 14) for Caesalpinioideae. After an ancient whole-genome duplication (WGD) shared by all Caesalpinioideae species (CWGD, ∼72.0 million years ago [MYA]), there were two recent successive WGD events, LWGD-1 (16.2-19.5 MYA) and LWGD-2 (7.1-9.5 MYA), in L. leucocephala. Thereafter, ∼40% gene loss and genome-size contraction have occurred during the diploidization process in L. leucocephala. To investigate secondary metabolites, we identified all gene copies involved in mimosine metabolism in these species and found that the abundance of mimosine biosynthesis genes in L. leucocephala largely explains its high mimosine production. We also identified the set of all potential genes involved in triterpenoid saponin biosynthesis in G. sinensis, which is more complete than that based on previous transcriptome-derived unigenes. Our results and genomic resources will facilitate biological studies of Caesalpinioideae and promote the utilization of valuable secondary metabolites.

Auxin efflux carrier PsPIN4 identified through genome-wide analysis as vital factor of petal abscission

PIN-FORMED (PIN) proteins, which function as efflux transporters, play many crucial roles in the polar transportation of auxin within plants. In this study, the exogenous applications of auxin IAA and TIBA were found to significantly prolong and shorten the florescence of tree peony (Paeonia suffruticosa Andr.) flowers. This finding suggests that auxin has some regulatory influence in petal senescence and abscission. Further analysis revealed a total of 8 PsPINs distributed across three chromosomes, which could be categorized into two classes based on phylogenetic and structural analysis. PsPIN1, PsPIN2a-b, and PsPIN4 were separated into the "long" PIN category, while PsPIN5, PsPIN6a-b, and PsPIN8 belonged to the "short" one. Additionally, the cis-regulatory elements of PsPIN promoters were associated with plant development, phytohormones, and environmental stress. These genes displayed tissue-specific expression, and phosphorylation sites were abundant throughout the protein family. Notably, PsPIN4 displayed distinct and elevated expression levels in roots, leaves, and flower organs. Expression patterns among the abscission zone (AZ) and adjacent areas during various flowering stages and IAA treatment indicate that PsPIN4 likely influences the initiation of peony petal abscission. The PsPIN4 protein was observed to be co-localized on both the plasma membrane and the cell nucleus. The ectopic expression of PsPIN4 reversed the premature flower organs abscission in the Atpin4 and significantly protracted florescence when introduced to Col Arabidopsis. Our findings established a strong basis for further investigation of PIN gene biological functions, particularly concerning intrinsic relationship between PIN-mediated auxin polar.

Genome-Wide Identification of Members of the Soybean CBL Gene Family and Characterization of the Functional Role of GmCBL1 in Responses to Saline and Alkaline Stress

Calcium ions function as key messengers in the context of intracellular signal transduction. The ability of plants to respond to biotic and abiotic stressors is highly dependent on the calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) signaling network. Here, a comprehensive effort was made to identify all members of the soybean CBL gene family, leading to the identification of 15 total genes distributed randomly across nine chromosomes, including 13 segmental duplicates. All the GmCBL gene subfamilies presented with similar gene structures and conserved motifs. Analyses of the expression of these genes in different tissues revealed that the majority of these GmCBLs were predominantly expressed in the roots. Significant GmCBL expression and activity increases were also observed in response to a range of stress-related treatments, including salt stress, alkaline stress, osmotic stress, or exposure to salicylic acid, brassinosteroids, or abscisic acid. Striking increases in GmCBL1 expression were observed in response to alkaline and salt stress. Subsequent analyses revealed that GmCBL1 was capable of enhancing soybean salt and alkali tolerance through the regulation of redox reactions. These results offer new insight into the complex mechanisms through which the soybean CBL gene family regulates the responses of these plants to environmental stressors, highlighting promising targets for efforts aimed at enhancing soybean stress tolerance.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Activation of the native PHYTOENE SYNTHASE 1 promoter by modifying near-miss cis-acting elements induces carotenoid biosynthesis in embryogenic rice callus

KEY MESSAGE

Modification of silent latent endosperm-enabled promoters (SLEEPERs) allows the ectopic activation of non-expressed metabolic genes in rice callus Metabolic engineering in plants typically involves transgene expression or the mutation of endogenous genes. An alternative is promoter modification, where small changes in the promoter sequence allow genes to be switched on or off in particular tissues. To activate silent genes in rice endosperm, we screened native promoters for near-miss cis-acting elements that can be converted to endosperm-active regulatory motifs. We chose rice PHYTOENE SYNTHASE 1 (PSY1), encoding the enzyme responsible for the first committed step in the carotenoid biosynthesis pathway, because it is not expressed in rice endosperm. We identified six motifs within a 120-bp region, upstream of the transcriptional start site, which differed from endosperm-active elements by up to four nucleotides. We mutated four motifs to match functional elements in the endosperm-active BCH2 promoter, and this promoter was able to drive GFP expression in callus and in seeds of regenerated plants. The 4 M promoter was not sufficient to drive PSY1 expression, so we mutated the remaining two elements and used the resulting 6 M promoter to drive PSY1 expression in combination with a PDS transgene. This resulted in deep orange callus tissue indicating the accumulation of carotenoids, which was subsequently confirmed by targeted metabolomics analysis. PSY1 expression driven by the uncorrected or 4 M variants of the promoter plus a PDS transgene produced callus that lacked carotenoids. These results confirm that the adjustment of promoter elements can facilitate the ectopic activation of endogenous plant promoters in rice callus and endosperm and most likely in other tissues and plant species.

C-terminally encoded peptide-like genes are associated with the development of primary root at qRL16.1 in soybean

Root architecture traits are belowground traits that harness moisture and nutrients from the soil and are equally important to above-ground traits in crop improvement. In soybean, the root length locus qRL16.1 was previously mapped on chromosome 16. The qRL16.1 has been characterized by transcriptome analysis of roots in near-isogenic lines (NILs), gene expression analysis in a pair of lines contrasting with alleles of qRL16.1, and differential gene expression analysis in germplasm accessions contrasting with root length. Two candidate genes, Glyma.16g108500 and Glyma.16g108700, have shown relatively higher expression in longer root accessions than in shorter rooting accessions. The C-terminal domain of Glyma.16g108500 and Glyma.16g108700 is similar to the conserved domain of C-terminally encoded peptides (CEPs) that regulate root length and nutrient response in Arabidopsis. Two polymorphisms upstream of Glyma.16g108500 showed a significant association with primary root length and total root length traits in a germplasm set. Synthetic peptide assay with predicted CEP variants of Glyma.16g108500 and Glyma.16g108700 demonstrated their positive effect on primary root length. The two genes are root-specific in the early stage of soybean growth and showed differential expression only in the primary root. These genes will be useful for improving soybean to develop a deep and robust root system to withstand low moisture and nutrient regimes.