The transition from primary siRNAs to amplified secondary siRNAs that regulate chalcone synthase during development of Glycine max seed coats

Reducing chlorophyll levels in seed-filling stages results in higher seed nitrogen without impacting canopy carbon assimilation

Chlorophyll is the major light-absorbing pigment for plant photosynthesis. While evolution has been selected for high chlorophyll content in leaves, previous work suggests that domesticated crops grown in modern high-density agricultural environments overinvest in chlorophyll production, thereby lowering light use and nitrogen use efficiency. To investigate the potential benefits of reducing chlorophyll levels, we created ethanol-inducible RNAi tobacco mutants that suppress Mg-chelatase subunit I (CHLI) with small RNA within 3 h of induction and reduce chlorophyll within 5 days in field conditions. We initiated chlorophyll reduction later in plant development to avoid the highly sensitive seedling stage and to allow young plants to have full green leaves to maximise light interception before canopy formation. This study demonstrated that leaf chlorophyll reduction >60% during seed-filling stages increased tobacco seed nitrogen concentration by as much as 17% while canopy photosynthesis, biomass and seed yields were maintained. These results indicate that time-specific reduction of chlorophyll could be a novel strategy that decouples the inverse relationship between yield and seed nitrogen by utilising saved nitrogen from the reduction of chlorophyll while maintaining full carbon assimilation capacity.

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.

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.

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.

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.

Dissecting the transcriptional regulation of proanthocyanidin and anthocyanin biosynthesis in soybean (Glycine max)

Proanthocyanidins (PAs), also known as condensed tannins, are plant natural products that are beneficial for human and livestock health. As one of the largest grown crops in the world, soybean (Glycine max) is widely used as human food and animal feed. Many cultivated soybeans with yellow seed coats lack PAs or anthocyanins, although some soybean cultivars have coloured seed coats that contain these compounds. Here, we analyse the transcriptional control of PA and anthocyanin biosynthesis in soybean. Ectopic expression of the transcription factors (TFs) GmTT2A, GmTT2B, GmMYB5A or R in soybean hairy roots induced the accumulation of PAs (primarily in phloem tissues) or anthocyanins and led to up-regulation of 1775, 856, 1411 and 1766 genes, respectively, several of which encode enzymes involved in PA biosynthesis. The genes regulated by GmTT2A and GmTT2B partially overlapped, suggesting conserved but potentially divergent roles for these two TFs in regulating PA accumulation in soybean. The two key enzymes anthocyanidin reductase and leucoanthocyanidin reductase were differentially upregulated, by GmTT2A/GmTT2B and GmMYB5A, respectively. Transgenic soybean plants overexpressing GmTT2B or MtLAP1 (a proven up-regulator of the upstream reactions for production of precursors for PA biosynthesis in legumes) showed increased accumulation of PAs and anthocyanins, respectively, associated with transcriptional reprogramming paralleling the RNA-seq data collected in soybean hairy roots. Collectively, our results show that engineered PA biosynthesis in soybean exhibits qualitative and spatial differences from the better-studied model systems Arabidopsis thaliana and Medicago truncatula, and suggest targets for engineering PAs in soybean plants.

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.

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.

Post-transcriptional silencing of chalcone synthase is involved in phenotypic lability in petals and leaves of bicolor dahlia (Dahlia variabilis) 'Yuino'

Post-transcriptional gene silencing (PTGS) of a chalcone synthase ( DvCHS2 ) occurred in the white part of bicolor petals and flavonoid-poor leaves; however, it did not in red petals and flavonoid-rich leaves. Petal color lability is a prominent feature of bicolor dahlia cultivars, and causes plants to produce not only original bicolor petals with colored bases and pure white tips, but also frequently single-colored petals without white tips. In this study, we analysed the molecular mechanisms that are associated with petal color lability using the red-white bicolor cultivar 'Yuino'. Red single-colored petals lose their white tips as a result of recover of flavonoid biosynthesis. Among flavonoid biosynthetic genes including four chalcone synthase (CHS)-like genes (DvCHS1, DvCHS2, DvCHS3, and DvCHS4), DvCHS1 and DvCHS2 had significantly lower expression levels in the white part of bicolor petals than in red petals, while DvCHS3, DvCHS4, and other flavonoid biosynthetic genes had almost the same expression levels. Small RNAs from the white part of a bicolor petal were mapped onto DvCHS1 and DvCHS2, while small RNAs from a red single-colored petal were not mapped onto any of the four CHS genes. A relationship between petal color and leaf flavonoid accumulation has previously been demonstrated, whereby red petal-producing plants accumulate flavonoids in their leaves, while bicolor petal-producing plants tend not to. The expression level of DvCHS2 was down-regulated in flavonoid-poor leaves and small RNAs from flavonoid-poor leaves were mapped onto DvCHS2, suggesting that the down-regulation of DvCHS2 in flavonoid-poor leaves occurs post-transcriptionally. Genomic analysis also suggested that DvCHS2 is the key gene involved in bicolor formation. Together, these results suggest that post-transcriptional gene silencing of DvCHS2 plays a key role in phenotypic lability in this bicolor dahlia.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Contemporary Understanding of miRNA-Based Regulation of Secondary Metabolites Biosynthesis in Plants

Plant's secondary metabolites such as flavonoids, terpenoids, and alkaloids etc. are known for their role in the defense against various insects-pests of plants and for medicinal benefits in human. Due to the immense biological importance of these phytochemicals, understanding the regulation of their biosynthetic pathway is crucial. In the recent past, advancement in the molecular technologies has enabled us to better understand the proteins, enzymes, genes, etc. involved in the biosynthetic pathway of the secondary metabolites. miRNAs are magical, tiny, non-coding ribonucleotides that function as critical regulators of gene expression in eukaryotes. Despite the accumulated knowledge of the miRNA-mediated regulation of several processes, the involvement of miRNAs in regulating secondary plant product biosynthesis is still poorly understood. Here, we summarize the recent progress made in the area of identification and characterizations of miRNAs involved in regulating the biosynthesis of secondary metabolites in plants and discuss the future perspectives for designing the viable strategies for their targeted manipulation.

Identification and Validation of Loci Governing Seed Coat Color by Combining Association Mapping and Bulk Segregation Analysis in Soybean

Soybean seed coat exists in a range of colors from yellow, green, brown, black, to bicolor. Classical genetic analysis suggested that soybean seed color was a moderately complex trait controlled by multi-loci. However, only a couple of loci could be detected using a single biparental segregating population. In this study, a combination of association mapping and bulk segregation analysis was employed to identify genes/loci governing this trait in soybean. A total of 14 loci, including nine novel and five previously reported ones, were identified using 176,065 coding SNPs selected from entire SNP dataset among 56 soybean accessions. Four of these loci were confirmed and further mapped using a biparental population developed from the cross between ZP95-5383 (yellow seed color) and NY279 (brown seed color), in which different seed coat colors were further dissected into simple trait pairs (green/yellow, green/black, green/brown, yellow/black, yellow/brown, and black/brown) by continuously developing residual heterozygous lines. By genotyping entire F2 population using flanking markers located in fine-mapping regions, the genetic basis of seed coat color was fully dissected and these four loci could explain all variations of seed colors in this population. These findings will be useful for map-based cloning of genes as well as marker-assisted breeding in soybean. This work also provides an alternative strategy for systematically isolating genes controlling relative complex trait by association analysis followed by biparental mapping.

Investigation of a miRNA-Induced Gene Silencing Technique in Petunia Reveals Alterations in miR173 Precursor Processing and the Accumulation of Secondary siRNAs from Endogenous Genes

MIGS (miRNA-induced gene silencing) is a straightforward and efficient gene silencing technique in Arabidopsis. It works by exploiting miR173 to trigger the production of phasiRNAs (phased small interfering RNAs). MIGS can be used in plant species other than Arabidopsis by co-expression of miR173 and target gene fragments fused to an upstream miR173 target site. However, the efficiency and technical mechanisms have not been thoroughly investigated in other plants. In this work, two vectors, pMIGS-chs and pMIGS-pds, were constructed and transformed into petunia plants. The transgenic plants showed CHS (chalcone synthase) and PDS (phytoene desaturase) gene-silencing phenotypes respectively, indicating that MIGS functions in petunia. MIGS-chs plants were used to investigate the mechanisms of this technique in petunia. Results of 5′- RACE showed that the miR173 target site was cleaved at the expected position and that endogenous CHS genes were cut at multiple positions. Small RNA deep sequencing analysis showed that the processing of Arabidopsis miR173 precursors in MIGS-chs transgenic petunia plants did not occur in exactly the same way as in Arabidopsis, suggesting differences in the machinery of miRNA processing between plant species. Small RNAs in-phase with the miR173 cleavage register were produced immediately downstream from the cleavage site and out-of-phase small RNAs were accumulated at relatively high levels from processing cycle 5 onwards. Secondary siRNAs were generated from multiple sites of endogenous CHS-A and CHS-J genes, indicating that miR173 cleavage induced siRNAs have the same ability to initiate siRNA transitivity as the siRNAs functioning in co-suppression and hpRNA silencing. On account of the simplicity of vector construction and the transitive amplification of signals from endogenous transcripts, MIGS is a good alternative gene silencing method for plants, especially for silencing a cluster of homologous genes with redundant functions.

The Glycine max cv. Enrei Genome for Improvement of Japanese Soybean Cultivars

We elucidated the genome sequence of Glycine max cv. Enrei to provide a reference for characterization of Japanese domestic soybean cultivars. The whole genome sequence obtained using a next-generation sequencer was used for reference mapping into the current genome assembly of G. max cv. Williams 82 obtained by the Soybean Genome Sequencing Consortium in the USA. After sequencing and assembling the whole genome shotgun reads, we obtained a data set with about 928 Mbs total bases and 60,838 gene models. Phylogenetic analysis provided glimpses into the ancestral relationships of both cultivars and their divergence from the complex that include the wild relatives of soybean. The gene models were analyzed in relation to traits associated with anthocyanin and flavonoid biosynthesis and an overall profile of the proteome. The sequence data are made available in DAIZUbase in order to provide a comprehensive informatics resource for comparative genomics of a wide range of soybean cultivars in Japan and a reference tool for improvement of soybean cultivars worldwide.