Endogenous, tissue-specific short interfering RNAs silence the chalcone synthase gene family in glycine max seed coats

Soybean2035: A decadal vision for soybean functional genomics and breeding

Soybean, the fourth most important crop in the world, uniquely serves as a source of both plant oil and plant protein for the world's food and animal feed. Although soybean production has increased approximately 13-fold over the past 60 years, the continually growing global population necessitates further increases in soybean production. In the past, especially in the last decade, significant progress has been made in both functional genomics and molecular breeding. However, many more challenges should be overcome to meet the anticipated future demand. Here, we summarize past achievements in the areas of soybean omics, functional genomics, and molecular breeding. Furthermore, we analyze trends in these areas, including shortages and challenges, and propose new directions, potential approaches, and possible outputs toward 2035. Our views and perspectives provide insight into accelerating the development of elite soybean varieties to meet the increasing demands of soybean production.

Dissecting the cellular architecture and genetic circuitry of the soybean seed

Seeds are complex structures composed of three regions, embryo, endosperm, and seed coat, with each further divided into subregions that consist of tissues, cell layers, and cell types. Although the seed is well characterized anatomically, much less is known about the genetic circuitry that dictates its spatial complexity. To address this issue, we profiled mRNAs from anatomically distinct seed subregions at several developmental stages. Analyses of these profiles showed that all subregions express similar diverse gene numbers and that the small gene numbers expressed subregion specifically provide information about the biological processes that occur in these seed compartments. In parallel, we profiled RNAs in individual nuclei and identified nuclei clusters representing distinct cell identities. Integrating single-nucleus RNA and subregion mRNA transcriptomes allowed most cell identities to be assigned to specific subregions and cell types and/or cell states. The number of cell identities exceeds the number of anatomically distinguishable cell types, emphasizing the spatial complexity of seeds. We defined gene coexpression networks that underlie distinct biological processes during seed development. We showed that network distribution among subregions and cell identities is highly variable. Some networks operate in single subregions and/or cell identities, and many coexpression networks operate in multiple subregions and/or cell identities. We also showed that single subregions and cell identities possess several networks. Together, our studies provide unique insights into the biological processes and genetic circuitry that underlie the spatial landscape of the seed.

Identifying Candidate Genes Related to Soybean (Glycine max) Seed Coat Color via RNA-Seq and Coexpression Network Analysis

BACKGROUND

The quality of soybeans is reflected in the seed coat color, which indicates soybean quality and commercial value. Researchers have identified genes related to seed coat color in various plants. However, research on the regulation of genes related to seed coat color in soybeans is rare.

METHODS

In this study, four lines of seed coats with different colors (medium yellow 14, black, green, and brown) were selected from the F2:5 population, with Beinong 108 as the female parent and green bean as the male parent, and the dynamic changes in the anthocyanins in the seed coat were stained with 4-dimethylaminocinnamaldehyde (DMACA) during the grain maturation process (20 days from grain drum to seed harvest). Through RNA-seq of soybean lines with four different colored seed coats at 30 and 50 days after seeding, we can further understand the key pathways and gene regulation modules between soybean seed coats of different colors.

RESULTS

DMACA revealed that black seed coat soybeans produce anthocyanins first and have the deepest staining. Clustering and principal component analysis (PCA) of the RNA-seq data divided the eight samples into two groups, resulting in 16,456 DEGs, including 5359 TFs. GO and KEGG enrichment analyses revealed that the flavonoid biosynthesis, starch and sucrose metabolism, carotenoid biosynthesis, and circadian rhythm pathways were significantly enriched. We also conducted statistical and expression pattern analyses on the differentially expressed transcription factors. Based on weighted gene coexpression network analysis (WGCNA), we identified seven specific modules that were significantly related to the four soybean lines with different seed coat colors. The connectivity and functional annotation of genes within the modules were calculated, and 21 candidate genes related to soybean seed coat color were identified, including six transcription factor (TF) genes and three flavonoid pathway genes.

CONCLUSIONS

These findings provide a theoretical basis for an in-depth understanding of the molecular mechanisms underlying differences in soybean seed coat color and provide new genetic resources.

ALTERED MERISTEM PROGRAM1 impairs RNA silencing by repressing the biogenesis of a subset of inverted repeat-derived siRNAs

RNA silencing negatively regulates gene expression at the transcriptional and posttranscriptional levels through DNA methylation, histone modification, mRNA cleavage, and translational inhibition. Small interfering RNAs (siRNAs) of 21 to 24 nucleotides are processed from double-stranded RNAs by Dicer-like (DCL) enzymes and play essential roles in RNA silencing in plants. Here, we demonstrated that ALTERED MERISTEM PROGRAM1 (AMP1) and its putative paralog LIKE AMP1 (LAMP1) impair RNA silencing by repressing the biogenesis of a subset of inverted repeat (IR)-derived siRNAs in Arabidopsis (Arabidopsis thaliana). AMP1 and LAMP1 inhibit Pol II-dependent IR gene transcription by suppressing ARGONAUTE 1 (AGO1) protein levels. Genetic analysis indicates that AMP1 acts upstream of RNA polymerase IV subunit 1 (NRPD1), RNA-dependent RNA polymerase 2 (RDR2), and DCL4, which are required for IR-induced RNA silencing. We also show that AMP1 and LAMP1 inhibit siRNA-mediated silencing in a different mechanism from that of AGO4 and DCL3. Together, these results reveal two previously unknown players in siRNA biogenesis from IRs-AGO1, which promotes IR transcription, and AMP1, which inhibits IR transcription indirectly through the repression of AGO1 expression.

Seed-Coat Pigmentation Plays a Crucial Role in Partner Selection and N2 Fixation in Legume-Root-Microbe Associations in African Soils

Legume-rhizobia symbiosis is the most important plant-microbe interaction in sustainable agriculture due to its ability to provide much needed N in cropping systems. This interaction is mediated by the mutual recognition of signaling molecules from the two partners, namely legumes and rhizobia. In legumes, these molecules are in the form of flavonoids and anthocyanins, which are responsible for the pigmentation of plant organs, such as seeds, flowers, fruits, and even leaves. Seed-coat pigmentation in legumes is a dominant factor influencing gene expression relating to N2 fixation and may be responsible for the different N2-fixing abilities observed among legume genotypes under field conditions in African soils. Common bean, cowpea, Kersting's groundnut, and Bambara groundnut landraces with black seed-coat color are reported to release higher concentrations of nod-gene-inducing flavonoids and anthocyanins compared with the Red and Cream landraces. Black seed-coat pigmentation is considered a biomarker for enhanced nodulation and N2 fixation in legumes. Cowpea, Bambara groundnut, and Kersting's bean with differing seed-coat colors are known to attract different soil rhizobia based on PCR-RFLP analysis of bacterial DNA. Even when seeds of the same legume with diverse seed-coat colors were planted together in one hole, the nodulating bradyrhizobia clustered differently in the PCR-RFLP dendrogram. Kersting's groundnut, Bambara groundnut, and cowpea with differing seed-coat colors were selectively nodulated by different bradyrhizobial species. The 16S rRNA amplicon sequencing also found significant selective influences of seed-coat pigmentation on microbial community structure in the rhizosphere of five Kersting's groundnut landraces. Seed-coat color therefore plays a dominant role in the selection of the bacterial partner in the legume-rhizobia symbiosis.

Pigmentation of soybean seed coats via a mutation that abolishes production of multiple-phased siRNAs of chalcone synthase genes

Lack of pigmentation in seed coats of soybean is caused by natural RNA silencing of chalcone synthase (CHS) genes. This phenomenon is an evolutionary consequence of structural changes in DNA that resulted in the production of double-stranded RNAs (dsRNAs) that trigger RNA degradation. Here we determined that a mutant with pigmented seed coats derived from a cultivar that lacked the pigmentation had a deletion between DNA regions ICHS1 and a cytochrome P450 gene; the deletion included GmIRCHS, a candidate gene that triggers CHS RNA silencing via production of CHS dsRNAs. We also characterized CHS short interfering RNAs (siRNAs) produced in the wild-type seed coats that had CHS RNA silencing. Phased 21-nt CHS siRNAs were detected in all 21 phases and were widely distributed in exon 2 of CHS7, which indicates commonality in the pattern of RNA degradation in natural CHS RNA silencing between distantly related species. These results with the similarities in the rearrangements found in spontaneous mutants suggest that the structural organization that generates dsRNAs that trigger phased siRNA production is vulnerable to further structural changes, which eventually abolish the induction of RNA silencing.

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

k-mer-based GWAS enhances the discovery of causal variants and candidate genes in soybean

Genome-wide association studies (GWAS) are powerful statistical methods that detect associations between genotype and phenotype at genome scale. Despite their power, GWAS frequently fail to pinpoint the causal variant or the gene controlling a given trait in crop species. Assessing genetic variants other than single-nucleotide polymorphisms (SNPs) could alleviate this problem. In this study, we tested the potential of structural variant (SV)- and k-mer-based GWAS in soybean by applying these methods as well as conventional SNP/indel-based GWAS to 13 traits. We assessed the performance of each GWAS approach based on loci for which the causal genes or variants were known from previous genetic studies. We found that k-mer-based GWAS was the most versatile approach and the best at pinpointing causal variants or candidate genes. Moreover, k-mer-based analyses identified promising candidate genes for loci related to pod color, pubescence form, and resistance to Phytophthora sojae. In our dataset, SV-based GWAS did not add value compared to k-mer-based GWAS and may not be worth the time and computational resources invested. Despite promising results, significant challenges remain regarding the downstream analysis of k-mer-based GWAS. Notably, better methods are needed to associate significant k-mers with sequence variation. Our results suggest that coupling k-mer- and SNP/indel-based GWAS is a powerful approach for discovering candidate genes in crop species.

An omics strategy increasingly improves the discovery of genetic loci and genes for seed-coat color formation in soybean

The phenotypic color of seeds is a complex agronomic trait and has economic and biological significance. The genetic control and molecular regulation mechanisms have been extensively studied. Here, we used a multi-omics strategy to explore the color formation in soybean seeds at a big data scale. We identified 13 large quantitative trait loci (QTL) for color with bulk segregating analysis in recombinant inbreeding lines. GWAS analysis of colors and decomposed attributes in 763 germplasms revealed associated SNP sites perfectly falling in five major QTL, suggesting inherited regulation on color during natural selection. Further transcriptomics analysis before and after color accumulation revealed 182 differentially expression genes (DEGs) in the five QTL, including known genes CHS, MYB, and F3'H involved in pigment accumulation. More DEGs with consistently upregulation or downregulation were identified as shared regulatory genes for two or more color formations while some DEGs were only for a specific color formation. For example, five upregulated DEGs in QTL qSC-3 were in flavonoid biosynthesis responsible for black and brown seed. The DEG (Glyma.08G085400) was identified in the purple seed only, which encodes gibberellin 2-beta-dioxygenase in the metabolism of colorful terpenoids. The candidate genes are involved in flavonoid biosynthesis, transcription factor regulation, gibberellin and terpenoid metabolism, photosynthesis, ascorbate and aldarate metabolism, and lipid metabolism. Seven differentially expressed transcription factors were also speculated that may regulate color formation, including a known MYB. The finds expand QTL and gene candidates for color formation, which could guide to breed better cultivars with designed colors.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01414-z.

Identification of candidate genes for soybean seed coat-related traits using QTL mapping and GWAS

Seed coat color is a typical morphological trait that can be used to reveal the evolution of soybean. The study of seed coat color-related traits in soybeans is of great significance for both evolutionary theory and breeding practices. In this study, 180 F₁₀ recombinant inbred lines (RILs) derived from the cross between the yellow-seed coat cultivar Jidou12 (ZDD23040, JD12) and the wild black-seed coat accession Y9 (ZYD02739) were used as materials. Three methods, single-marker analysis (SMA), interval mapping (IM), and inclusive composite interval mapping (ICIM), were used to identify quantitative trait loci (QTLs) controlling seed coat color and seed hilum color. Simultaneously, two genome-wide association study (GWAS) models, the generalized linear model (GLM) and mixed linear model (MLM), were used to jointly identify seed coat color and seed hilum color QTLs in 250 natural populations. By integrating the results from QTL mapping and GWAS analysis, we identified two stable QTLs (qSCC02 and qSCC08) associated with seed coat color and one stable QTL (qSHC08) related to seed hilum color. By combining the results of linkage analysis and association analysis, two stable QTLs (qSCC02, qSCC08) for seed coat color and one stable QTL (qSHC08) for seed hilum color were identified. Upon further investigation using Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, we validated the previous findings that two candidate genes (CHS3C and CHS4A) reside within the qSCC08 region and identified a new QTL, qSCC02. There were a total of 28 candidate genes in the interval, among which Glyma.02G024600, Glyma.02G024700, and Glyma.02G024800 were mapped to the glutathione metabolic pathway, which is related to the transport or accumulation of anthocyanin. We considered the three genes as potential candidate genes for soybean seed coat-related traits. The QTLs and candidate genes detected in this study provide a foundation for further understanding the genetic mechanisms underlying soybean seed coat color and seed hilum color and are of significant value in marker-assisted breeding.

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

QTL Mapping and Transcriptome Analysis Reveal Candidate Genes Regulating Seed Color in Brassica napus

Yellow seeds are desirable in rapeseed breeding because of their higher oil content and better nutritional quality than black seeds. However, the underlying genes and formation mechanism of yellow seeds remain unclear. Here, a novel yellow-seeded rapeseed line (Huangaizao, HAZ) was crossed with a black-seeded rapeseed line (Zhongshuang11, ZS11) to construct a mapping population of 196 F₂ individuals, based on which, a high-density genetic linkage map was constructed. This map, comprising 4174 bin markers, was 1618.33 cM in length and had an average distance of 0.39 cM between its adjacent markers. To assess the seed color of the F₂ population, three methods (imaging, spectrophotometry, and visual scoring) were used and a common major quantitative trait locus (QTL) on chromosome A09, explaining 10.91-21.83% of the phenotypic variance, was detected. Another minor QTL, accounting for 6.19-6.69% of the phenotypic variance, was detected on chromosome C03, only by means of imaging and spectrophotometry. Furthermore, a dynamic analysis of the differential expressions between the parental lines showed that flavonoid biosynthesis-related genes were down-regulated in the yellow seed coats at 25 and 35 days after flowering. A coexpression network between the differentially expressed genes identified 17 candidate genes for the QTL intervals, including a flavonoid structure gene, novel4557 (BnaC03.TT4), and two transcription factor genes, namely, BnaA09G0616800ZS (BnaA09.NFYA8) and BnaC03G0060200ZS (BnaC03.NAC083), that may regulate flavonoid biosynthesis. Our study lays a foundation for further identifying the genes responsible for and understanding the regulatory mechanism of yellow seed formation in Brassica napus.

Transcriptional Profile of Soybean Seeds with Contrasting Seed Coat Color

Soybean is the primary source of vegetable protein and is used for various purposes, mainly to feed animals. This crop can have diverse seed coat colors, varying from yellow, black, brown, and green to bicolor. Black seed coat cultivars have already been assigned as favorable for both seed and grain production. Thus, this work aimed to identify genes associated with soybean seed quality by comparing the transcriptomes of soybean seeds with contrasting seed coat colors. The results from RNA-seq analyses were validated with real-time PCR using the cultivar BRS 715A (black seed coat) and the cultivars BRS 413 RR and DM 6563 IPRO (yellow seed coat). We found 318 genes differentially expressed in all cultivars (freshly harvested seeds and seeds stored in cold chamber). From the in silico analysis of the transcriptomes, the following genes were selected and validated with RT-qPCR: ACS1, ACSF3, CYP90A1, CYP710A1, HCT, CBL, and SAHH. These genes are genes induced in the black seed coat cultivar and are part of pathways responsible for ethylene, lipid, brassinosteroid, lignin, and sulfur amino acid biosynthesis. The BRSMG 715A gene has almost 4times more lignin than the yellow seed coat cultivars. These attributes are related to the BRSMG 715A cultivar's higher seed quality, which translates to more longevity and resistance to moisture and mechanical damage. Future silencing studies may evaluate the knockout of these genes to better understand the biology of soybean seeds with black seed coat.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

Novel Seed Size: A Novel Seed-Developing Gene in Glycine max

Soybean-seed development is controlled in multiple ways, as in many known regulating genes. Here, we identify a novel gene, Novel Seed Size (NSS), involved in seed development, by analyzing a T-DNA mutant (S006). The S006 mutant is a random mutant of the GmFTL4pro:GUS transgenic line, with phenotypes with small and brown seed coats. An analysis of the metabolomics and transcriptome combined with RT-qPCR in the S006 seeds revealed that the brown coat may result from the increased expression of chalcone synthase 7/8 genes, while the down-regulated expression of NSS leads to small seed size. The seed phenotypes and a microscopic observation of the seed-coat integument cells in a CRISPR/Cas9-edited mutant nss1 confirmed that the NSS gene conferred small phenotypes of the S006 seeds. As mentioned in an annotation on the Phytozome website, NSS encodes a potential DNA helicase RuvA subunit, and no such genes were previously reported to be involved in seed development. Therefore, we identify a novel gene in a new pathway controlling seed development in soybeans.

Epigenetic management of self and non-self: lessons from 40 years of transgenic plants

Plant varieties exhibiting unstable or variegated phenotypes, or showing virus recovery have long remained a mystery. It is only with the development of transgenic plants 40 years ago that the epigenetic features underlying these phenomena were elucidated. Indeed, the study of transgenic plants that did not express the introduced sequences revealed that transgene loci sometimes undergo transcriptional gene silencing (TGS) or post-transcriptional gene silencing (PTGS) by activating epigenetic defenses that naturally control transposable elements, duplicated genes or viruses. Even when they do not trigger TGS or PTGS spontaneously, stably expressed transgenes driven by viral promoters set apart from endogenous genes in their epigenetic regulation. As a result, transgenes driven by viral promoters are capable of undergoing systemic PTGS throughout the plant, whereas endogenous genes can only undergo local PTGS in cells where RNA quality control is impaired. Together, these results indicate that the host genome distinguishes self from non-self at the epigenetic level, allowing PTGS to eliminate non-self, and preventing PTGS to become systemic and kill the plant when it is locally activated against deregulated self.

The dicing activity of DCL3 and DCL4 is negatively affected by flavonoids

The dicing activities of DCL3 and DCL4 are inhibited by accumulated metabolites in soybean leaves. Epicatechin and 7,4'-dihydroxyflavone inhibited Arabidopsis DCL3 and DCL4 in vitro. Flavonoids are major secondary metabolites in plants, and soybean (Glycine max L.) is a representative plant that accumulates flavonoids, including isoflavonoids, to high levels. Naturally-occurring RNA interference (RNAi) against the chalcone synthase (CHS) gene represses flavonoid (anthocyanin) biosynthesis in an organ-specific manner, resulting in a colorless (yellow) seed coat in many soybean cultivars. To better understand seed coat-specific naturally-occurring RNAi in soybean, we characterized soybean Dicer-like (DCL) 3 and 4, which play critical roles in RNAi. Using a previously established dicing assay, two dicing activities producing 24- and 21-nt siRNAs, corresponding to DCL3 and DCL4, respectively, were detected in soybean. Dicing activity was detected in colorless seed coats where RNAi against CHS genes was found, but no dicing activity was detected in leaves where CHS expression was prevalent. Biochemical analysis revealed that soybean leaves contained two types of inhibitors effective for Arabidopsis Dicers (AtDCL3 and AtDCL4), one of which was a heat-labile high molecular weight compound of 50 to 100 kD while another was a low molecular weight substance. We found that some flavonoids, such as epicatechin and 7,4'-dihydroxyflavone, inhibited both AtDCL3 and AtDCL4, but AtDCL4 was more sensitive to these flavonoids than AtDCL3. These results suggest that flavonoids inhibit the dicing activity of DCL4 and thereby attenuate RNAi in soybean leaves.

The Key to the Future Lies in the Past: Insights from Grain Legume Domestication and Improvement Should Inform Future Breeding Strategies

Crop domestication is a co-evolutionary process that has rendered plants and animals significantly dependent on human interventions for survival and propagation. Grain legumes have played an important role in the development of Neolithic agriculture some 12,000 years ago. Despite being early companions of cereals in the origin and evolution of agriculture, the understanding of grain legume domestication has lagged behind that of cereals. Adapting plants for human use has resulted in distinct morpho-physiological changes between the wild ancestors and domesticates, and this distinction has been the focus of several studies aimed at understanding the domestication process and the genetic diversity bottlenecks created. Growing evidence from research on archeological remains, combined with genetic analysis and the geographical distribution of wild forms, has improved the resolution of the process of domestication, diversification and crop improvement. In this review, we summarize the significance of legume wild relatives as reservoirs of novel genetic variation for crop breeding programs. We describe key legume features, which evolved in response to anthropogenic activities. Here, we highlight how whole genome sequencing and incorporation of omics-level data have expanded our capacity to monitor the genetic changes accompanying these processes. Finally, we present our perspective on alternative routes centered on de novo domestication and re-domestication to impart significant agronomic advances of novel crops over existing commodities. A finely resolved domestication history of grain legumes will uncover future breeding targets to develop modern cultivars enriched with alleles that improve yield, quality and stress tolerance.

Identification of ST1 reveals a selection involving hitchhiking of seed morphology and oil content during soybean domestication

Seed morphology and quality of cultivated soybean (Glycine max) have changed dramatically during domestication from their wild relatives, but their relationship to selection is poorly understood. Here, we describe a semi-dominant locus, ST1 (Seed Thickness 1), affecting seed thickness and encoding a UDP-D-glucuronate 4-epimerase, which catalyses UDP-galacturonic acid production and promotes pectin biosynthesis. Interestingly, this morphological change concurrently boosted seed oil content, which, along with up-regulation of glycolysis biosynthesis modulated by ST1, enabled soybean to become a staple oil crop. Strikingly, ST1 and an inversion controlling seed coat colour formed part of a single selective sweep. Structural variation analysis of the region surrounding ST1 shows that the critical mutation in ST1 existed in earlier wild relatives of soybean and the region containing ST1 subsequently underwent an inversion, which was followed by successive selection for both traits through hitchhiking during selection for seed coat colour. Together, these results provide direct evidence that simultaneously variation for seed morphology and quality occurred earlier than variation for seed coat colour during soybean domestication. The identification of ST1 thus sheds light on a crucial phase of human empirical selection in soybeans and provides evidence that our ancestors improved soybean based on taste.

Biosynthesis of Secondary Metabolites Based on the Regulation of MicroRNAs

MicroRNA (miRNA), a noncoding ribonucleic acid, is considered to be important for the progression of gene expression in plants and animals by rupture or translational repression of targeted mRNAs. Many types of miRNA regulate plant metabolism, growth, and response to biotic and abiotic factors. miRNA characterization helps to expose its function in regulating the process of post-transcriptional genetic regulation. There are a lot of factors associated with miRNA function, but the function of miRNA in the organic synthesis of by-products by natural products is not yet fully elucidated. The current review is aimed at observing and characterizing miRNAs and identifying those involved in the functioning of the biosynthesis of secondary metabolites in plants, with their use in controlled manipulation.

Post-transcriptional gene silencing of the chalcone synthase gene CHS causes corolla lobe-specific whiting of Japanese gentian

Post-transcriptional gene silencing of the chalcone synthase gene CHS specifically suppresses anthocyanin biosynthesis in corolla lobes and is responsible for the formation of a stripe type bicolor in Japanese gentian. The flower of Japanese gentian is a bell-shaped corolla composed of lobes and plicae, which is painted uniformly blue. However, the gentian cultivar 'Hakuju' shows bicolor phenotype (blue-white stripe corolla), in which anthocyanin accumulation is suppressed only in corolla lobes. Expression analysis indicated that steady-state levels of chalcone synthase (CHS) transcripts were remarkably reduced in corolla lobes compared with plicae during petal pigmentation initiation. However, no significant difference in expression levels of other flavonoid biosynthetic structural and regulatory genes was detected in its lobes and plicae. On feeding naringenin in white lobes, anthocyanin accumulation was recovered. Northern blotting probed with CHS confirmed the abundant accumulation of small RNAs in corolla lobes. Likewise, small RNA-seq analysis indicated that short reads from its lobes were predominantly mapped onto the 2nd exon region of the CHS gene, whereas those from the plicae were scarcely mapped. Subsequent infection with the gentian ovary ringspot virus (GORV), which had an RNA-silencing activity, showed the recovery of partial pigmentation in lobes. Hence, these results strongly suggested that suppressing anthocyanin accumulation in the lobes of bicolored 'Hakuju' was attributed to the specific degradation of CHS mRNA in corolla lobes, which was through post-transcriptional gene silencing (PTGS). Herein, we revealed the molecular mechanism of strip bicolor formation in Japanese gentian, and showed that PTGS of CHS was also responsible for flower color pattern in a floricultural plant other than petunia and dahlia.

Post-transcriptional gene silencing of CYP76AD controls betalain biosynthesis in bracts of bougainvillea

Betalain is one of four major plant pigments and shares some features with anthocyanin; however, no plant has been found to biosynthesize both pigments. Previous studies have reported that anthocyanin biosynthesis in some plants is regulated by post-transcriptional gene-silencing (PTGS), but the importance of PTGS in betalain biosynthesis remains unclear. In this study, we report the occurrence of PTGS in betalain biosynthesis in bougainvillea (Bougainvillea peruviana) 'Thimma', which produces bracts of three different color on the same plant, namely pink, white, and pink-white. This resembles the unstable anthocyanin pigmentation phenotype that is associated with PTGS, and hence we anticipated the presence of PTGS in the betalain biosynthetic pathway. To test this, we analysed pigments, gene expression, small RNAs, and transient overexpression. Our results demonstrated that PTGS of BpCYP76AD1, a gene encoding one of the betalain biosynthesis enzymes, is responsible for the loss of betalain biosynthesis in 'Thimma'. Neither the genetic background nor DNA methylation in the BpCYP76AD1 sequence could explain the induction of PTGS, implying that another locus controls the unstable pigmentation. Our results indicate that naturally occurring PTGS contributes to the diversification of color patterns not only in anthocyanin biosynthesis but also in betalain biosynthesis.

LIRBase: a comprehensive database of long inverted repeats in eukaryotic genomes

Small RNAs (sRNAs) constitute a large portion of functional elements in eukaryotic genomes. Long inverted repeats (LIRs) can be transcribed into long hairpin RNAs (hpRNAs), which can further be processed into small interfering RNAs (siRNAs) with vital biological roles. In this study, we systematically identified a total of 6 619 473 LIRs in 424 eukaryotic genomes and developed LIRBase (https://venyao.xyz/lirbase/), a specialized database of LIRs across different eukaryotic genomes aiming to facilitate the annotation and identification of LIRs encoding long hpRNAs and siRNAs. LIRBase houses a comprehensive collection of LIRs identified in a wide range of eukaryotic genomes. In addition, LIRBase not only allows users to browse and search the identified LIRs in any eukaryotic genome(s) of interest available in GenBank, but also provides friendly web functionalities to facilitate users to identify LIRs in user-uploaded sequences, align sRNA sequencing data to LIRs, perform differential expression analysis of LIRs, predict mRNA targets for LIR-derived siRNAs, and visualize the secondary structure of candidate long hpRNAs encoded by LIRs. As demonstrated by two case studies, collectively, LIRBase bears the great utility for systematic investigation and characterization of LIRs and functional exploration of potential roles of LIRs and their derived siRNAs in diverse species.

Increased copy number of gibberellin 2-oxidase 8 genes reduced trailing growth and shoot length during soybean domestication

Copy number variations (CNVs) play important roles in crop domestication. However, there is only very limited information on the involvement of CNVs in soybean domestication. Trailing growth and long shoots are soybean adaptations for natural habitats but cause lodging that hampers yield in cultivation. Previous studies have focused on Dt1/2 affecting the indeterminate/determinate growth habit, whereas the possible role of the gibberellin pathway remained unclear. In the present study, quantitative trait locus (QTL) mapping of a recombinant inbred population of 460 lines revealed a trailing-growth-and-shoot-length QTL. A CNV region within this QTL was identified, featuring the apical bud-expressed gibberellin 2-oxidase 8A/B, the copy numbers of which were positively correlated with expression levels and negatively with trailing growth and shoot length, and their effects were demonstrated by transgenic soybean and Arabidopsis thaliana. Based on the fixation index, this CNV region underwent intense selection during the initial domestication process.

Biosynthetic Pathway of Proanthocyanidins in Major Cash Crops

Proanthocyanidins (PAs) are a group of oligomers or polymers composed of monomeric flavanols. They offer many benefits for human fitness, such as antioxidant, anticancer, and anti-inflammatory activities. To date, three types of PA have been observed in nature: procyanidins, propelargonidins, and prodelphinidins. These are synthesized as some of the end-products of the flavonoid pathway by different consecutive enzymatic activities, from the same precursor—naringenin. Although the general biosynthetic pathways of PAs have been reported in a few model plant species, little is known about the species-specific pathways in major crops containing different types of PA. In the present study, we identified the species-specific pathways in 10 major crops, based on the presence/absence of flavanol-based intermediates in the metabolic pathway, and found 202 orthologous genes in the reference genomic database of each species, which may encode for key enzymes involved in the biosynthetic pathways of PAs. Parallel enzymatic reactions in the pathway are responsible for the ratio between PAs and anthocyanins, as well as among the three types of PAs. Our study suggests a promising strategy for molecular breeding, to regulate the content of PAs and anthocyanins and improve the nutritional quality of food sources globally.

Recent Advances in Heterologous Synthesis Paving Way for Future Green-Modular Bioindustries: A Review With Special Reference to Isoflavonoids

Isoflavonoids are well-known plant secondary metabolites that have gained importance in recent time due to their multiple nutraceutical and pharmaceutical applications. In plants, isoflavonoids play a role in plant defense and can confer the host plant a competitive advantage to survive and flourish under environmental challenges. In animals, isoflavonoids have been found to interact with multiple signaling pathways and have demonstrated estrogenic, antioxidant and anti-oncologic activities in vivo. The activity of isoflavonoids in the estrogen pathways is such that the class has also been collectively called phytoestrogens. Over 2,400 isoflavonoids, predominantly from legumes, have been identified so far. The biosynthetic pathways of several key isoflavonoids have been established, and the genes and regulatory components involved in the biosynthesis have been characterized. The biosynthesis and accumulation of isoflavonoids in plants are regulated by multiple complex environmental and genetic factors and interactions. Due to this complexity of secondary metabolism regulation, the export and engineering of isoflavonoid biosynthetic pathways into non-endogenous plants are difficult, and instead, the microorganisms Saccharomyces cerevisiae and Escherichia coli have been adapted and engineered for heterologous isoflavonoid synthesis. However, the current ex-planta production approaches have been limited due to slow enzyme kinetics and traditionally laborious genetic engineering methods and require further optimization and development to address the required titers, reaction rates and yield for commercial application. With recent progress in metabolic engineering and the availability of advanced synthetic biology tools, it is envisaged that highly efficient heterologous hosts will soon be engineered to fulfill the growing market demand.

Alleles disrupting LBD37-like gene by an 136 bp insertion show different distributions between green and purple cabbages (Brassica oleracea var. capitata)

BACKGROUND

In Arabidopsis thaliana (Arabidopsis), clade IIb lateral organ boundary domain (LBD) 37, 38, and 39 proteins negatively regulate anthocyanin biosynthesis and affect nitrogen responses.

OBJECTIVE

To investigate the possible role of LBD genes in anthocyanin accumulations among green and purple cabbages (Brassica oleracea var. capitata), we determined sequence variations and expression levels of cabbage homologs of Arabidopsis LBD37, 38, and 39.

METHODS

DNA and mRNA sequences of BoLBD37, BoLBD37L (BoLBD37-like), BoLBD38, BoLBD38L (BoLBD38-like), and BoLBD39 gene in cabbage were determined. Allelic variations of BoLBD37L alleles in cabbages, resulting from insertions, were validated by genomic DNA PCR. Gene expressions were examined by semi-quantitative reverse transcription (RT-PCR) or quantitative RT-PCR.

RESULTS

Based on the expression analyses, BoLBD37L with two alleles, BoLBD37L-G and BoLBD37L-P, was selected as a candidate gene important for differential anthocyanin accumulation. BoLBD37L-P contains an 136 base pair insertion in the 2nd exon, producing two splicing variants encoding truncated proteins. Most purple cabbage lines were found to have BoLBD37L-P allele as homozygotes or heterozygotes, and only two out of sixty-four purple cabbages were identified as BoLBD37L-G homozygotes. Expression analyses of anthocyanin biosynthesis genes and their upstream regulators, including BoLBD37L, suggest that truncated proteins encoded by splicing variants of BoLBD37L-P may disrupt the BoLBD37L function as repressor.

CONCLUSION

Difference in the C-terminal region of BoLBD37L-G and BolBD37L-P may affect the expression of downstream genes, BoMYB114L and BoTT8, resulting in differential anthocyanin accumulation.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Contrasting epigenetic control of transgenes and endogenous genes promotes post-transcriptional transgene silencing in Arabidopsis

Transgenes that are stably expressed in plant genomes over many generations could be assumed to behave epigenetically the same as endogenous genes. Here, we report that whereas the histone H3K9me2 demethylase IBM1, but not the histone H3K4me3 demethylase JMJ14, counteracts DNA methylation of Arabidopsis endogenous genes, JMJ14, but not IBM1, counteracts DNA methylation of expressed transgenes. Additionally, JMJ14-mediated specific attenuation of transgene DNA methylation enhances the production of aberrant RNAs that readily induce systemic post-transcriptional transgene silencing (PTGS). Thus, the JMJ14 chromatin modifying complex maintains expressed transgenes in a probationary state of susceptibility to PTGS, suggesting that the host plant genome does not immediately accept expressed transgenes as being epigenetically the same as endogenous genes.

Accumulating evidences point to a discrepancy in the epigenetic behaviour of transgenes and endogenous genes. Here, via characterization of mutants impaired in histone demethylases JMJ14 and IBM1, the authors show that transgenes and endogenous genes are regulated by different epigenetic mechanisms in Arabidopsis.

Soybean DICER-LIKE2 Regulates Seed Coat Color via Production of Primary 22-Nucleotide Small Interfering RNAs from Long Inverted Repeats

In plants, 22-nucleotide small RNAs trigger the production of secondary small interfering RNAs (siRNAs) and enhance silencing. DICER-LIKE2 (DCL2)-dependent 22-nucleotide siRNAs are rare in Arabidopsis (Arabidopsis thaliana) and are thought to function mainly during viral infection; by contrast, these siRNAs are abundant in many crops such as soybean (Glycine max) and maize (Zea mays). Here, we studied soybean 22-nucleotide siRNAs by applying CRISPR-Cas9 to simultaneously knock out the two copies of soybean DCL2, GmDCL2a and GmDCL2b, in the Tianlong1 cultivar. Small RNA sequencing revealed that most 22-nucleotide siRNAs are derived from long inverted repeats (LIRs) and disappeared in the Gmdcl2a/2b double mutant. De novo assembly of a Tianlong1 reference genome and transcriptome profiling identified an intronic LIR formed by the chalcone synthase (CHS) genes CHS1 and CHS3 This LIR is the source of primary 22-nucleotide siRNAs that target other CHS genes and trigger the production of secondary 21-nucleotide siRNAs. Disruption of this process in Gmdcl2a/2b mutants substantially increased CHS mRNA levels in the seed coat, thus changing the coat color from yellow to brown. Our results demonstrated that endogenous LIR-derived transcripts in soybean are predominantly processed by GmDCL2 into 22-nucleotide siRNAs and uncovered a role for DCL2 in regulating natural traits.

Dissecting seed pigmentation-associated genomic loci and genes by employing dual approaches of reference-based and k-mer-based GWAS with 438 Glycine accessions

The soybean is agro-economically the most important among all cultivated legume crops, and its seed color is considered one of the most attractive factors in the selection-by-breeders. Thus, genome-wide identification of genes and loci associated with seed colors is critical for the precision breeding of crop soybeans. To dissect seed pigmentation-associated genomic loci and genes, we employed dual approaches by combining reference-based genome-wide association study (rbGWAS) and k-mer-based reference-free GWAS (rfGWAS) with 438 Glycine accessions. The dual analytical strategy allowed us to identify four major genomic loci (designated as SP1-SP4 in this study) associated with the seed colors of soybeans. The k-mer analysis enabled us to find an important recombination event that occurred between subtilisin and I-cluster B in the soybean genome, which could describe a special structural feature of ii allele within the I locus (SP3). Importantly, mapping analyses of both mRNAs and small RNAs allowed us to reveal that the subtilisin-CHS1/CHS3 chimeric transcripts generate and act as an initiator towards 'mirtron (i.e., intron-harboring miRNA precursor)'-triggered silencing of chalcone synthase (CHS) genes. Consequently, the results led us to propose a working model of 'mirtron-triggered gene silencing (MTGS)' to elucidate a long-standing puzzle in the genome-wide CHS gene silencing mechanism. In summary, our study reports four major genomic loci, lists of key genes and genome-wide variations that are associated with seed pigmentation in soybeans. In addition, we propose that the MTGS mechanism plays a crucial role in the genome-wide silencing of CHS genes, thereby suggesting a clue to currently predominant soybean cultivars with the yellow seed coat. Finally, this study will provide a broad insight into the interactions and correlations among seed color-associated genes and loci within the context of anthocyanin biosynthetic pathways.

Overexpression of an ethylene-forming ACC oxidase (ACO) gene precedes the Minute Hilum seed coat phenotype in Glycine max

BACKGROUND

To elucidate features of seed development, we investigated the transcriptome of a soybean isoline from the germplasm collection that contained an introgressed allele known as minute hilum (mi) which confers a smaller hilum region where the seed attaches to the pod and also results in seed coat cracking surrounding the hilum region.

RESULTS

RNAs were extracted from immature seed from an extended hilum region (i.e., the hilum and a small ring of tissue surrounding the hilum in which the cracks form) at three different developmental stages:10-25, 25-50 and 50-100 mg seed fresh weight in two independent replicates for each stage. The transcriptomes of these samples from both the Clark isoline containing the mi allele (PI 547628, UC413, ii R t mi G), and its recurrent Clark 63 parent isoline (PI 548532, UC7, ii R T Mi g), which was used for six generations of backcrossing, were compared for differential expression of 88,648 Glyma models of the soybean genome Wm82.a2. The RNA sequence data obtained from the 12 cDNA libraries were subjected to padj value < 0.05 and at least two-fold expression differences to select with confidence genes differentially expressed in the hilum-containing tissue of the seed coat between the two lines. Glyma.09G008400 annotated as encoding an ethylene forming enzyme, ACC oxidase (ACO), was found to be highly overexpressed in the mi hilum region at 165 RPKMs (reads per kilobase per million mapped reads) compared to the standard line at just 0.03 RPKMs. Evidence of changes in expression of genes downstream of the ethylene pathway included those involved in auxin and gibberellin hormone action and extensive differences in expression of cell wall protein genes. These changes are postulated to determine the restricted hilum size and cracking phenotypes.

CONCLUSIONS

We present transcriptome and phenotypic evidence that substantially higher expression of an ethylene-forming ACO gene likely shifts hormone balance and sets in motion downstream changes resulting in a smaller hilum phenotype and the cracks observed in the minute hilum (mi) isoline as compared to its recurrent parent.

The Effects of Domestication on Secondary Metabolite Composition in Legumes

Legumes are rich in secondary metabolites, such as polyphenols, alkaloids, and saponins, which are important defense compounds to protect the plant against herbivores and pathogens, and act as signaling molecules between the plant and its biotic environment. Legume-sourced secondary metabolites are well known for their potential benefits to human health as pharmaceuticals and nutraceuticals. During domestication, the color, smell, and taste of crop plants have been the focus of artificial selection by breeders. Since these agronomic traits are regulated by secondary metabolites, the basis behind the genomic evolution was the selection of the secondary metabolite composition. In this review, we will discuss the classification, occurrence, and health benefits of secondary metabolites in legumes. The differences in their profiles between wild legumes and their cultivated counterparts will be investigated to trace the possible effects of domestication on secondary metabolite compositions, and the advantages and drawbacks of such modifications. The changes in secondary metabolite contents will also be discussed at the genetic level to examine the genes responsible for determining the secondary metabolite composition that might have been lost due to domestication. Understanding these genes would enable breeding programs and metabolic engineering to produce legume varieties with favorable secondary metabolite profiles for facilitating adaptations to a changing climate, promoting beneficial interactions with biotic factors, and enhancing health-beneficial secondary metabolite contents for human consumption.

Synthesis of Full-Length cDNA Infectious Clones of Soybean Mosaic Virus and Functional Identification of a Key Amino Acid in the Silencing Suppressor Hc-Pro

Soybean mosaic virus (SMV), which belongs to the Potyviridae, causes significant reductions in soybean yield and seed quality. In this study, both tag-free and reporter gene green fluorescent protein (GFP)-containing infectious clones for the SMV N1 strain were constructed by Gibson assembly and with the yeast homologous recombination system, respectively. Both infectious clones are suitable for agroinfiltration on the model host N. benthamiana and show strong infectivity for the natural host soybean and several other legume species. Both infectious clones were seed transmitted and caused typical virus symptoms on seeds and progeny plants. We used the SMV-GFP infectious clone to further investigate the role of key amino acids in the silencing suppressor helper component-proteinase (Hc-Pro). Among twelve amino acid substitution mutants, the co-expression of mutant 2-with an Asparagine→Leucine substitution at position 182 of the FRNK (Phe-Arg-Asn-Lys) motif-attenuated viral symptoms and alleviated the host growth retardation caused by SMV. Moreover, the Hc-Prom2 mutant showed stronger oligomerization than wild-type Hc-Pro. Taken together, the SMV infectious clones will be useful for studies of host-SMV interactions and functional gene characterization in soybeans and related legume species, especially in terms of seed transmission properties. Furthermore, the SMV-GFP infectious clone will also facilitate functional studies of both virus and host genes in an N. benthamiana transient expression system.

Pan-Genome of Wild and Cultivated Soybeans

Soybean is one of the most important vegetable oil and protein feed crops. To capture the entire genomic diversity, it is needed to construct a complete high-quality pan-genome from diverse soybean accessions. In this study, we performed individual de novo genome assemblies for 26 representative soybeans that were selected from 2,898 deeply sequenced accessions. Using these assembled genomes together with three previously reported genomes, we constructed a graph-based genome and performed pan-genome analysis, which identified numerous genetic variations that cannot be detected by direct mapping of short sequence reads onto a single reference genome. The structural variations from the 2,898 accessions that were genotyped based on the graph-based genome and the RNA sequencing (RNA-seq) data from the representative 26 accessions helped to link genetic variations to candidate genes that are responsible for important traits. This pan-genome resource will promote evolutionary and functional genomics studies in soybean.

Single-base deletion in GmCHR5 increases the genistein-to-daidzein ratio in soybean seed

Novel mutant alleles related to isoflavone content are useful for breeding programs to improve the disease resistance and nutritional content of soybean. However, identification of mutant alleles from high-density mutant libraries is expensive and time-consuming because soybean has a large, complicated genome. Here, we identified the gene responsible for increased genistein-to-daidzein ratio in seed of the mutant line F333ES017D9. For this purpose, we used a time- and cost-effective approach based on selective genotyping of a small number of F₂ plants showing the mutant phenotype with nearest-neighboring-nucleotide substitution-high-resolution melting analysis markers, followed by alignment of short reads obtained by next-generation sequencing analysis with the identified locus. In the mutant line, GmCHR5 harbored a single-base deletion that caused a change in the substrate flow in the isoflavone biosynthetic pathway towards genistein. Mutated GmCHR5 was expressed at a lower level during seed development than wild-type GmCHR5. Ectopic overexpression of GmCHR5 increased the production of daidzein derivatives in both the wild-type and mutant plants. The present strategy will be useful for accelerating identification of mutant alleles responsible for traits of interest in agronomically important crops.

Identification of iron-chelating phenolics contributing to seed coat coloration in soybeans (Glycine max (L.) Merr.) expressing aryloxyalkanoate dioxygenase-12

Soybeans (Glycine max (L.) Merr.) genetically modified to express aryloxyalkanoate dioxygenase-12 (AAD-12), an enzyme that confers resistance to the herbicide 2,4-D, can sometimes exhibit a darker seed coat coloration than equivalent unmodified soybeans. The biochemical basis for this coloration was investigated in a non-commercial transgenic event, DAS-411Ø4-7 that exhibited more pronounced AAD-12-associated seed coat coloration than the commercial event, DAS-444Ø6-6. Analysis of color-enriched seed coat fractions from DAS-411Ø4-7 showed that the color was due to localized accumulation of iron-chelating phenolics, particularly the isoflavone genistin, that are associated with seed coat pectic polysaccharide and produce a brown chromophore. The association between genistin, iron, and pectic polysaccharide was characterized using a variety of analytical methods. Darker seeds from commercial soybean event DAS-444Ø6-6 also show higher genistin content localized to the darker colored portions of the seed coat (with no increase in whole seed genistin levels).

Construction and comparison of three reference-quality genome assemblies for soybean

We report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

Host-derived gene silencing of parasite fitness genes improves resistance to soybean cyst nematodes in stable transgenic soybean

KEY MESSAGE

Soybean expressing small interfering RNA of SCN improved plant resistance to SCN consistently, and small RNA-seq analysis revealed a threshold of siRNA expression required for resistance ability. Soybean cyst nematode (SCN), Heterodera glycines, is one of the most destructive pests limiting soybean production worldwide, with estimated losses of $1 billion dollars annually in the USA alone. RNA interference (RNAi) has become a powerful tool for silencing gene expression. We report here that the expression of hairpin RNAi constructs, derived from two SCN genes related to reproduction and fitness, HgY25 and HgPrp17, enhances resistance to SCN in stably transformed soybean plants. The analyses of T3 to T5 generations of stable transgenic soybeans by molecular strategies and next-generation sequencing confirmed the presence of specific short interfering RNAs complementary to the target SCN genes. Bioassays performed on transgenic soybean lines targeting SCN HgY25 and HgPrp17 fitness genes showed significant reductions (up to 73%) for eggs/g root in the T3 and T4 homozygous transgenic lines. Targeted mRNAs of SCN eggs collected from the transgenic soybean lines were efficiently down-regulated, as confirmed by quantitative RT-PCR. Based on the small RNA-seq data and bioassays, it is our hypothesis that a threshold of small interfering RNA molecules is required to significantly reduce SCN populations feeding on the host plants. Our results demonstrated that host-derived gene silencing of essential SCN fitness genes could be an effective strategy for enhancing resistance in crop plants.

Alteration of flavonoid pigmentation patterns during domestication of food crops

Flavonoids are plant pigments that provide health benefits for human and animal consumers. Understanding why domesticated crops have altered pigmentation patterns and unraveling the molecular/genetic mechanisms that underlie this will facilitate the breeding of new (healthier) varieties. We present an overview of changes in flavonoid pigmentation patterns that have occurred during crop domestication and, where possible, link them to the molecular changes that brought about the new phenotypes. We consider species that lost flavonoid pigmentation in the edible part of the plant at some point during domestication (like cereals). We also consider the converse situation, for example eggplant (aubergine), which instead gained strong anthocyanin accumulation in the skin of the fruit during domestication, and some varieties of citrus and apple that acquired anthocyanins in the fruit flesh. Interestingly, the genes responsible for such changes are sometimes closely linked to, or have pleiotropic effects on, important domestication genes, suggesting accidental and perhaps inevitable changes of anthocyanin patterning during domestication. In other cases, flavonoid pigmentation patterns in domesticated crops are the result of cultural preferences, with examples being found in varieties of citrus, barley, wheat, and maize. Finally, and more recently, in some species, anthocyanins seem to have been the direct target of selection in a second wave of domestication that followed the introduction of industrial food processing.

Nonallelic homologous recombination events responsible for copy number variation within an RNA silencing locus

The structure of chalcone synthase (CHS) gene repeats in different alleles of the I (inhibitor) locus in soybean spawns endogenous RNA interference (RNAi) that leads to phenotypic change in seed coat color of this major agronomic crop. Here, we examined CHS gene copy number by digital PCR and single nucleotide polymorphisms (SNPs) through whole genome resequencing of 15 cultivars that varied in alleles of the I locus (I, ii , ik , and i) that control the pattern distribution of pigments in the seed coats. Lines homozygous for the ii allele had the highest copy number followed by the I and ik cultivars which were more related to each other than to the lines with ii alleles. Some of the recessive i alleles were spontaneous mutations, and each revealed a loss of copy number by digital PCR relative to the parent varieties. Amplicon sequencing and whole genome resequencing determined that the breakpoints of several ii to i mutations resulted from nonallelic homologous recombination (NAHR) events between CHS genes located in segmental duplications leading to large 138-kilobase deletions that erase the structure generating the CHS siRNAs along with eight other non-CHS genes. Functional hybrid CHS genes (designated CHS5:1) were formed in the process and represent rare examples of NAHR in higher plants that have been captured by examining spontaneous mutational events in isogenic mutant lines.

Integrated physiological and genomic analysis reveals structural variations and expression patterns of candidate genes for colored- and green-leaf poplar

Colored-leaf plants are increasingly popular and have been attracting more and more attentions. However, the molecular mechanism of leaf coloration in plants has not been fully understood. In this study, a colored-leaf cultivar of Populus deltoides (Caihong poplar, CHP) and green-leaf cultivar of Populus deltoides L2025 were used to explore the mechanism of leaf coloration through physiological and the whole genome resequencing analysis. The content of anthocyanins, total Chl, and carotenoids in the leaves of CHP and L2025 were evaluated. The ratio of anthocyanins to total Chl in CHP was 25.0 times higher than that in L2025; this could be attributed to the red leaf color of CHP. Based on the whole genome resequencing analysis, 951,421 polymorphic SNPs and 221,907 indels were screened between CHP and L2025. Using qRT-PCR analysis, three structural genes (flavonol synthase 1 family protein, UDP-glucose flavonoid 3-O-glucosyltransferase 3′ and flavonoid 3-O-galactosyl transferase family protein) and six transcription factors (MYB-related protein Myb4, transcription factor GAMYB, PtrMYB179, transcription factor bHLH53, transcription factor bHLH3, VARICOSE family protein) may be involved in the anthocyanin synthesis pathway, which could be used as candidate genes to explore the molecular regulation mechanism of leaf coloration in Populus deltoids, and could be used in molecular breeding in the future.

Isolation and Characterization of Three Chalcone Synthase Genes in Pecan (Carya illinoinensis)

Phenolics are a group of important plant secondary metabolites that have been proven to possess remarkable antioxidant activity and to be beneficial for human health. Pecan nuts are an excellent source of dietary phenolics. In recent years, many studies have focused on the separation and biochemical analysis of pecan phenolics, but the molecular mechanisms of phenolic metabolism in pecans have not been fully elucidated, which significantly hinders quality breeding research for this plant. Chalcone synthase (CHS) plays crucial roles in phenolic biosynthesis. In this study, three Carya illinoinensisCHSs (CiCHS1, CiCHS2, and CiCHS3), were isolated and analyzed. CiCHS2 and CiCHS3 present high expression levels in different tissues, and they are also highly expressed at the initial developmental stages of kernels in three pecan genotypes. A correlation analysis was performed between the phenolic content and CHSs expression values during kernel development. The results indicated that the expression variations of CiCHS2 and CiCHS3 are significantly related to changes in total phenolic content. Therefore, CiCHSs play crucial roles in phenolic components synthesis in pecan. We believe that the isolation of CiCHSs is helpful for understanding phenolic metabolism in C. illinoinensis, which will improve quality breeding and resistance breeding studies in this plant.

A reference-grade wild soybean genome

Efficient crop improvement depends on the application of accurate genetic information contained in diverse germplasm resources. Here we report a reference-grade genome of wild soybean accession W05, with a final assembled genome size of 1013.2 Mb and a contig N50 of 3.3 Mb. The analytical power of the W05 genome is demonstrated by several examples. First, we identify an inversion at the locus determining seed coat color during domestication. Second, a translocation event between chromosomes 11 and 13 of some genotypes is shown to interfere with the assignment of QTLs. Third, we find a region containing copy number variations of the Kunitz trypsin inhibitor (KTI) genes. Such findings illustrate the power of this assembly in the analysis of large structural variations in soybean germplasm collections. The wild soybean genome assembly has wide applications in comparative genomic and evolutionary studies, as well as in crop breeding and improvement programs.

Comprehensive RNA sequencing and co-expression network analysis to complete the biosynthetic pathway of coumestrol, a phytoestrogen

Coumestrol (CMS), a coumestan isoflavone, plays key roles in nodulation through communication with rhizobia, and has been used as phytoestrogens for hormone replacement therapy in humans. Because CMS content is controlled by multiple genetic factors, the genetic basis of CMS biosynthesis has remained unclear. We identified soybean genotypes with consistently high (Daewonkong) or low (SS0903-2B-21-1-2) CMS content over 2 years. We performed RNA sequencing of leaf samples from both genotypes at developmental stage R7, when CMS levels are highest. Within the phenylpropanoid biosynthetic pathway, 41 genes were tightly connected in a functional co-expression gene network; seven of these genes were differentially expressed between two genotypes. We identified 14 candidate genes involved in CMS biosynthesis. Among them, seven were annotated as encoding oxidoreductases that may catalyze the transfer of electrons from daidzein, a precursor of CMS. Two of the other genes, annotated as encoding a MYB domain protein and a MLP-like protein, may increase CMS accumulation in response to stress conditions. Our results will help to complete our understanding of the CMS biosynthetic pathway, and should facilitate development of soybean cultivars with high CMS content that could be used to promote the fitness of plants and human beings.

Genome-wide identification and localization of chalcone synthase family in soybean (Glycine max [L]Merr)

BACKGROUND

Soybean is a paleopolyploid that has undergone two whole genome duplication events. Gene duplication is a type of genomic change that can lead to novel functions of pre-existing genes. Chalcone synthase (CHS) is the plant-specific type III polyketide synthase that catalyzes the first committed step in (iso)flavonoid biosynthesis in plants.

RESULTS

Here we performed a genome-wide search of CHS genes in soybean, and identified 21 GmCHS loci containing 14 unique GmCHS (GmCHS1-GmCHS14) that included 5 newly identified GmCHSs (GmCHS10-GmCHS14). Furthermore, 3 copies of GmCHS3 and 2 copies of GmCHS4 were found in soybean. Analysis of gene structure of GmCHSs revealed the presence of a single intron in protein-coding regions except for GmCHS12 that contained 3 introns. Even though GmCHS genes are located on 8 different chromosomes, a large number of these genes are present on chromosome 8 where they form 3 distinct clusters. Expression analysis of GmCHS genes revealed tissue-specific expression pattern, and that some GmCHS isoforms localize in the cytoplasm and the nucleus while other isoforms are restricted to cytoplasm only.

CONCLUSION

Overall, we have identified 21 GmCHS loci with 14 unique GmCHS genes in the soybean genome. Their gene structures and genomic organization together with the spatio-temporal expression and protein localization suggest their importance in the production of downstream metabolites such as (iso)flavonoids and their derived phytoalexins.

Dissecting genomic hotspots underlying seed protein, oil, and sucrose content in an interspecific mapping population of soybean using high-density linkage mapping

The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high-resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K-SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K-SNP markers. Importantly, the bin map derived from next-generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication-related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome-wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL-based haplotype using whole-genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

Two amino acid changes in the R3 repeat cause functional divergence of two clustered MYB10 genes in peach

R2R3-MYB genes play a pivotal role in regulating anthocyanin accumulation. Here, we report two tandemly duplicated R2R3-MYB genes in peach, PpMYB10.1 and PpMYB10.2, with the latter showing lower ability to induce anthocyanin accumulation than the former. Site-directed mutation assay revealed two amino acid changes in the R3 repeat, Arg/Lys⁶⁶ and Gly/Arg⁹³, responsible for functional divergence between these two PpMYB10 genes. Anthocyanin-promoting activity of PpMYB10.2 was significantly increased by a single amino acid replacement of Arg⁹³ with Gly⁹³. However, either the Gly⁹³ → Arg⁹³ or Arg⁶⁶ → Lys⁶⁶ substitutions alone showed little impact on anthocyanin-promoting activity of PpMYB10.1, but simultaneous substitutions caused a significant decrease. Reciprocal substitution of Arg/Gly⁹³ could significantly alter binding affinity to PpbHLH3, while the Arg⁶⁶ → Lys⁶⁶ substitution is predicted to affect the folding of the MYB DNA-binding domain, instead of PpbHLH3-binding affinity. Overall, the change of anthocyanin-promoting activity was accompanied with that of bHLH-binding affinity, suggesting that DNA-binding affinity of R2R3-MYBs depends on their bHLH partners. Our study is helpful for understanding of functional evolution of R2R3-MYBs and their interaction with DNA targets.

Purple Brassica oleracea var. capitata F. rubra is due to the loss of BoMYBL2-1 expression

BACKGROUND

Water-soluble anthocyanin pigments are important ingredients in health-improving supplements and valuable for the food industry. Although great attention has been paid to the breeding and production of crops containing high levels of anthocyanin, genetic variation in red or purple cabbages (Brassica oleracea var. capitata F. rubra) has not yet been characterized at the molecular level. In this study, we identified the mechanism responsible for the establishment of purple color in cabbages.

RESULTS

BoMYBL2-1 is one of the regulatory genes in the anthocyanin biosynthesis pathway in cabbages. It is a repressor whose expression is inversely correlated to anthocyanin synthesis and is not detectable in purple cabbages. Sequence analysis of purple cabbages revealed that most lacked BoMYBL2-1 coding sequences, although a few had a substitution in the region of the promoter 347 bp upstream of the gene that was associated with an absence of BoMYBL2-1 expression. Lack of transcriptional activity of the substitution-containing promoter was confirmed using transgenic Arabidopsis plants transformed with promoter::GUS fusion constructs. The finding that the defect in BoMYBL2-1 expression was solely responsible for purple coloration in cabbages was further demonstrated using genomic PCR and RT-PCR analyses of many other structural and regulatory genes in anthocyanin biosynthesis. Molecular markers for purple cabbages were developed and validated using 69 cabbage lines.

CONCLUSION

Expression of BoMYBL2-1 was inversely correlated to anthocyanin content, and purple color in cabbages resulted from a loss of BoMYBL2-1 expression, caused by either the promoter substitution or deletion of the gene. This is the first report of molecular markers that distinguish purple cabbages. Such markers will be useful for the production of intraspecific and interspecific hybrids for functional foods, and for industrial purposes requiring high anthocyanin content.