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

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.

Impact of pod and seed photosynthesis on seed filling and canopy carbon gain in soybean

There is a limited understanding of the carbon assimilation capacity of nonfoliar green tissues and its impact on yield and seed quality since most photosynthesis research focuses on leaf photosynthesis. In this study, we investigate the photosynthetic efficiency of soybean (Glycine max) pods and seeds in a field setting and evaluate its effect on mature seed weight and composition. We demonstrate that soybean pod and seed photosynthesis contributes 13% to 14% of the mature seed weight. Carbon assimilation by soybean pod and seed photosynthesis can compensate for 81% of carbon loss through the respiration of the same tissues, and our model predicts that soybean pod and seed photosynthesis contributes up to 9% of the total daily carbon gain of the canopy. Chlorophyll fluorescence (CF) shows that the operating efficiency of photosystem II in immature soybean seeds peaks at the 10 to 100 mg seed weight stage, while that of immature pods peaks at the 75 to 100 mg stage. This study provides quantitative information about the efficiency of soybean pod and seed photosynthesis during tissue development and its impact on yield.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

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.

Natural and artificial sources of genetic variation used in crop breeding: A baseline comparator for genome editing

Traditional breeding has successfully selected beneficial traits for food, feed, and fibre crops over the last several thousand years. The last century has seen significant technological advancements particularly in marker assisted selection and the generation of induced genetic variation, including over the last few decades, through mutation breeding, genetic modification, and genome editing. While regulatory frameworks for traditional varietal development and for genetic modification with transgenes are broadly established, those for genome editing are lacking or are still evolving in many regions. In particular, the lack of “foreign” recombinant DNA in genome edited plants and that the resulting SNPs or INDELs are indistinguishable from those seen in traditional breeding has challenged development of new legislation. Where products of genome editing and other novel breeding technologies possess no transgenes and could have been generated via traditional methods, we argue that it is logical and proportionate to apply equivalent legislative oversight that already exists for traditional breeding and novel foods. This review analyses the types and the scale of spontaneous and induced genetic variation that can be selected during traditional plant breeding activities. It provides a base line from which to judge whether genetic changes brought about by techniques of genome editing or other reverse genetic methods are indeed comparable to those routinely found using traditional methods of plant breeding.

Identification of genetic loci conferring seed coat color based on a high-density map in soybean

Seed coat color is a typical evolutionary trait. Identification of the genetic loci that control seed coat color during the domestication of wild soybean could clarify the genetic variations between cultivated and wild soybean. We used 276 F₁₀ recombinant inbred lines (RILs) from the cross between a cultivated soybean (JY47) and a wild soybean (ZYD00321) as the materials to identify the quantitative trait loci (QTLs) for seed coat color. We constructed a high-density genetic map using re-sequencing technology. The average distance between adjacent markers was 0.31 cM on this map, comprising 9,083 bin markers. We identified two stable QTLs (qSC08 and qSC11) for seed coat color using this map, which, respectively, explained 21.933 and 26.934% of the phenotypic variation. Two candidate genes (CHS3C and CHS4A) in qSC08 were identified according to the parental re-sequencing data and gene function annotations. Five genes (LOC100786658, LOC100801691, LOC100806824, LOC100795475, and LOC100787559) were predicted in the novel QTL qSC11, which, according to gene function annotations, might control seed coat color. This result could facilitate the identification of beneficial genes from wild soybean and provide useful information to clarify the genetic variations for seed coat color in cultivated and wild soybean.

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.

A major gene for tolerance to cold-induced seed coat discoloration relieves viral seed mottling in soybean

In yellow soybeans, inhibition of seed coat pigmentation by RNA silencing of CHS genes is suppressed by low temperature and a viral suppressor, resulting in 'cold-induced seed coat discoloration' and 'seed mottling', respectively. Differences exist in the degree of cold-induced seed coat discoloration among Japanese yellow soybean cultivars; for example, Toyomusume is sensitive, Toyohomare has some tolerance, and Toyoharuka is highly tolerant. In this study, we compared the degree of seed mottling severity due to soybean mosaic virus (SMV) among these three soybean cultivars. Obvious differences were found, with the order of severity as follows: Toyohomare > Toyomusume > Toyoharuka. RNA gel blot analysis indicated that CHS transcript abundance in the seed coat, which was increased by SMV infection, was responsible for the severity of seed mottling. Quantitative reverse transcription PCR analysis revealed why mottling was most severe in SMV-infected Toyohomare: the SMV titer in its seed coat was higher than in the other two infected cultivars. We further suggest that a major gene (Ic) for tolerance to cold-induced seed coat discoloration can relieve the severity of seed mottling in SMV-infected Toyoharuka.

Gynoecy instability in cucumber (Cucumis sativus L.) is due to unequal crossover at the copy number variation-dependent Femaleness (F) locus

Cucumber, Cucumis sativus is an important vegetable crop, and gynoecy has played a critical role in yield increase of hybrid cucumber production. Cucumber has a unique genetic system for gynoecious sex expression, which is determined by the copy number variation (CNV)-based, dominant, and dosage-dependent femaleness (F) locus. However, this gynoecy expression system seems unstable since monecious plants could often be found in F-dependent gynoecious cucumber inbreds. We hypothesized that gynoecy instability (gynoecy loss) may be due to unequal crossing over (UCO) during meiosis among repeat units of the CNV. In this study, using high throughput genome resequencing, fiber-FISH and genomic qPCR analyses, we first confirmed and refined the structure of the F locus, which was a CNV of a 30.2-kb tandem repeat. Gynoecious plants contained three genes: CsACS1, CsACS1G, and CsMYB, of which CsACS1G is a duplication of CsACS1 but with a recombinant distal promoter that may contribute to gynoecy sex expression. In two large populations from self-pollinated gynoecious inbred lines, 'gynoecy loss' mutants were identified with similar mutation rates (~0.12%). We show that these monecious mutants have lost CsACS1G. In addition, we identified gynoecious lines in natural populations that carry two copies of CSACS1G. We proposed a model to explain gynoecy instability in F-dependent cucumbers, which is caused by UCO among CSACS1/G units during meiosis. The findings present a convincing case that the phenotypic variation of an economically important trait is associated with the dynamic changes of copy numbers at the F locus. This work also has important implications in cucumber breeding.