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

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

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.

Argonaute 5-mediated antiviral defense and viral counter-defense in Nicotiana benthamiana

The argonaute (AGO) family proteins play a crucial role in preventing viral invasions through the plant antiviral RNA silencing pathway, with distinct AGO proteins recruited for specific antiviral mechanisms. Our previous study revealed that Nicotiana benthamiana AGO5 (NbAGO5) expression was significantly upregulated in response to bamboo mosaic virus (BaMV) infection. However, the roles of NbAGO5 in antiviral mechanisms remained to be explored. In this research, we examined the antiviral functions of NbAGO5 in the infections of different viruses. It was found that the accumulation of NbAGO5 was induced not only at the RNA but also at the protein level following the infections of BaMV, potato virus X (PVX), tobacco mosaic virus (TMV), and cucumber mosaic virus (CMV) in N. benthamiana. To explore the antiviral mechanism and regulatory function of NbAGO5, we generated NbAGO5 overexpression (OE-NbAGO5) and knockout (nbago5) transgenic N. benthamiana lines. Our findings reveal that NbAGO5 provides defense against BaMV, PVX, TMV, and a mutant CMV deficient in 2b gene, but not against the wild-type CMV and turnip mosaic virus (TuMV). Through affinity purification and small RNA northern blotting, we demonstrated that NbAGO5 exerts its antiviral function by binding to viral small interfering RNAs (vsiRNAs). Moreover, we observed that CMV 2b and TuMV HC-Pro interact with NbAGO5, triggering its degradation via the 26S proteasome and autophagy pathways, thereby allowing these viruses to overcome NbAGO5-mediated defense. In addition, TuMV HC-Pro provides another line of counter-defense by interfering with vsiRNA binding by NbAGO5. Our study provides further insights into the antiviral RNA interference mechanism and the complex interplay between NbAGO5 and plant viruses.

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.

Genome-Wide Identification and Characterization of Argonaute, Dicer-like and RNA-Dependent RNA Polymerase Gene Families and Their Expression Analyses in Fragaria spp

In the growth and development of plants, some non-coding small RNAs (sRNAs) not only mediate RNA interference at the post-transcriptional level, but also play an important regulatory role in chromatin modification at the transcriptional level. In these processes, the protein factors Argonaute (AGO), Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) play very important roles in the synthesis of sRNAs respectively. Though they have been identified in many plants, the information about these gene families in strawberry was poorly understood. In this study, using a genome-wide analysis and a phylogenetic approach, 13 AGO, six DCL, and nine RDR genes were identified in diploid strawberry Fragaria vesca. We also identified 33 AGO, 18 DCL, and 28 RDR genes in octoploid strawberry Fragaria × ananassa, studied the expression patterns of these genes in various tissues and developmental stages of strawberry, and researched the response of these genes to some hormones, finding that almost all genes respond to the five hormone stresses. This study is the first report of a genome-wide analysis of AGO, DCL, and RDR gene families in Fragaria spp., in which we provide basic genomic information and expression patterns for these genes. Additionally, this study provides a basis for further research on the functions of these genes and some evidence for the evolution between diploid and octoploid strawberries.

Argonaute5 and its associated small RNAs modulate the transcriptional response during the rhizobia-Phaseolus vulgaris symbiosis

Both plant- and rhizobia-derived small RNAs play an essential role in regulating the root nodule symbiosis in legumes. Small RNAs, in association with Argonaute proteins, tune the expression of genes participating in nodule development and rhizobial infection. However, the role of Argonaute proteins in this symbiosis has been overlooked. In this study, we provide transcriptional evidence showing that Argonaute5 (AGO5) is a determinant genetic component in the root nodule symbiosis in Phaseolus vulgaris. A spatio-temporal transcriptional analysis revealed that the promoter of PvAGO5 is active in lateral root primordia, root hairs from rhizobia-inoculated roots, nodule primordia, and mature nodules. Transcriptional analysis by RNA sequencing revealed that gene silencing of PvAGO5 affected the expression of genes involved in the biosynthesis of the cell wall and phytohormones participating in the rhizobial infection process and nodule development. PvAGO5 immunoprecipitation coupled to small RNA sequencing revealed the small RNAs bound to PvAGO5 during the root nodule symbiosis. Identification of small RNAs associated to PvAGO5 revealed miRNAs previously known to participate in this symbiotic process, further supporting a role for AGO5 in this process. Overall, the data presented shed light on the roles that PvAGO5 plays during the root nodule symbiosis in P. vulgaris.

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 F10 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.

Analyses on the pigment composition of different seed coat colors in adzuki bean

Seed coat color is an important quality and domestication trait. The adzuki bean has more than a dozen seed coat colors closely associated with the anthocyanin and flavonoid metabolism pathways. In this study, we explored the pigment composition of 10 different seed coat color adzuki beans including red, black mottle on red, black mottle on gray, golden, green, black, ivory, brown, and light brown. The results showed that anthocyanins are the main pigment in adzuki bean seed coat. There were no carotenoid or pelargonidin derivatives in the seed coats of any accessions. Different colors of adzuki bean seed coat have different pigment compositions and the combination of procyanidins and anthocyanins affected seed coat color. The ivory seed coat had an extremely low proanthocyanidin and anthocyanin content. Only the green adzuki bean seed coats contained chlorophyll. Our results explain the pigment composition of the different seed coat colors and the combination of proanthocyanidins and anthocyanins affected seed coat color in adzuki bean. These results can provide a theoretical basis for the study of adzuki bean coloring mechanism.

Functional Divergence of G and Its Homologous Genes for Green Pigmentation in Soybean Seeds

The degradation of chlorophyll in mature soybean seeds is closely related to the development of their yellow color. In this study, we examined G, its homologue G-like (GL), and their mutant alleles and investigated the relationship between these genes and chlorophyll accumulation in the seed coats of mature seeds. Transient expression of G and GL proteins fused with green fluorescent protein revealed that both were localized in plastids. Overexpression of G resulted in the accumulation of chlorophyll in the seed coats and cotyledons of mature seeds, indicating that high expression levels of G result in chlorophyll accumulation that exceeds its metabolism in the seeds of yellow soybean. Analysis of near isogenic lines at the G locus demonstrated a significant difference in the chlorophyll content of the seed coats and cotyledons of mature seeds when G and mutant g alleles were expressed in the d1d2 stay-green genetic background, indicating that the G protein might repress the SGR-independent degradation of chlorophyll. We examined the distribution of mutant alleles at the G and GL loci among cultivated and wild soybean germplasm. The g allele was widely distributed in cultivated soybean germplasm, except for green seed coat soybean lines, all of which contained the G allele. The gl alleles were much fewer in number than the g alleles and were mainly distributed in the genetic resources of cultivated soybean from Japan. None of the landraces and breeding lines investigated in this study were observed to contain both the g and gl alleles. Therefore, in conclusion, the mutation of the G locus alone is essential for establishing yellow soybeans, which are major current soybean breeding lines.

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.

MYB transcription factors GmMYBA2 and GmMYBR function in a feedback loop to control pigmentation of seed coat in soybean

Soybean has undergone extensive selection pressures for seed nutrient composition and seed color during domestication, but the major genetic loci controlling seed coat color have not been completely understood, and the transcriptional regulation relationship among the loci remains elusive. Here, two major regulators, GmMYBA2 and GmMYBR, were functionally characterized as an anthocyanin activator and repressor, respectively. Ectopic expression of GmMYBA2 in soybean hairy roots conferred the enhanced accumulation of delphinidin and cyanidin types of anthocyanins in W1t and w1T backgrounds, respectively, through activating anthocyanin biosynthetic genes in the reported loci. The seed coat pigmentation of GmMYBA2-overexpressing transgenic plants in the W1 background mimicked the imperfect black phenotype (W1/w1, i, R, t), suggesting that GmMYBA2 was responsible for the R locus. Molecular and biochemical analysis showed that GmMYBA2 interacted with GmTT8a to directly activate anthocyanin biosynthetic genes. GmMYBA2 and GmMYBR might form a feedback loop to fine-tune seed coat coloration, which was confirmed in transgenic soybeans. Both GmTT8a and GmMYBR that were activated by GmMYBA2 in turn enhanced and obstructed the formation of the GmMYBA2-GmTT8a module, respectively. The results revealed the sophisticated regulatory network underlying the soybean seed coat pigmentation loci and shed light on the understanding of the seed coat coloration and other seed inclusions.

Characterization of the Common Genetic Basis Underlying Seed Hilum Size, Yield, and Quality Traits in Soybean

Developing high yielding cultivars with outstanding quality traits are perpetual objectives throughout crop breeding operations. Confoundingly, both of these breeding objectives typically involve working with complex quantitative traits that can be affected by genetic and environmental factors. Establishing correlations of these complex traits with more easily identifiable and highly heritable traits can simplify breeding processes. In this study, two parental soybean genotypes contrasting in seed hilum size, yield, and seed quality, as well as 175 F₉ recombinant inbred lines (RILs) derived from these parents, were grown in 3 years. The h² _b of four hilum size, two quality and two yield traits, ranged from 0.72 to 0.87. The four observed hilum size traits exhibited significant correlation (P < 0.05) with most of seed yield and quality traits, as indicated by correlation coefficients varying from -0.35 to 0.42, which suggests that hilum size could be considered as a proxy trait for soybean yield and quality. Interestingly, among 53 significant quantitative trait loci (QTLs) with logarithm of odds (LOD) values ranging from 2.51 to 6.69 and accounting for 6.40-16.10% of genetic variation, three loci encoding hilum size, qSH6.2, qSH8, and qSH10, colocated with QTLs for seed yield and quality traits, demonstrating that genes impacting seed hilum size colocalize in part with genes acting on soybean yield and quality. As a result of the breeding efforts and field observations described in this work, it is reasonable to conclude that optimizing hilum size through selection focused on a few QTLs may be useful for breeding new high yielding soybean varieties with favorable quality characteristics.

Mining Grapevine Downy Mildew Susceptibility Genes: A Resource for Genomics-Based Breeding and Tailored Gene Editing

Several pathogens continuously threaten viticulture worldwide. Until now, the investigation on resistance loci has been the main trend to understand the interaction between grapevine and the mildew causal agents. Dominantly inherited gene-based resistance has shown to be race-specific in some cases, to confer partial immunity, and to be potentially overcome within a few years since its introgression. Recently, on the footprint of research conducted in Arabidopsis, putative genes associated with downy mildew susceptibility have been discovered also in the grapevine genome. In this work, we deep-sequenced four putative susceptibility genes-namely VvDMR6.1, VvDMR6.2, VvDLO1, VvDLO2-in 190 genetically diverse grapevine genotypes to discover new sources of broad-spectrum and recessively inherited resistance. Identified Single Nucleotide Polymorphisms were screened in a bottleneck analysis from the genetic sequence to their impact on protein structure. Fifty-five genotypes showed at least one impacting mutation in one or more of the scouted genes. Haplotypes were inferred for each gene and two of them at the VvDMR6.2 gene were found significantly more represented in downy mildew resistant genotypes. The current results provide a resource for grapevine and plant genetics and could corroborate genomic-assisted breeding programs as well as tailored gene editing approaches for resistance to biotic stresses.

Whole-genome sequence diversity and association analysis of 198 soybean accessions in mini-core collections

We performed whole-genome Illumina resequencing of 198 accessions to examine the genetic diversity and facilitate the use of soybean genetic resources and identified 10 million single nucleotide polymorphisms and 2.8 million small indels. Furthermore, PacBio resequencing of 10 accessions was performed, and a total of 2,033 structure variants were identified. Genetic diversity and structure analysis congregated the 198 accessions into three subgroups (Primitive, World, and Japan) and showed the possibility of a long and relatively isolated history of cultivated soybean in Japan. Additionally, the skewed regional distribution of variants in the genome, such as higher structural variations on the R gene clusters in the Japan group, suggested the possibility of selective sweeps during domestication or breeding. A genome-wide association study identified both known and novel causal variants on the genes controlling the flowering period. Novel candidate causal variants were also found on genes related to the seed coat colour by aligning together with Illumina and PacBio reads. The genomic sequences and variants obtained in this study have immense potential to provide information for soybean breeding and genetic studies that may uncover novel alleles or genes involved in agronomically important traits.

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.

What does it take to be antiviral? An Argonaute-centered perspective on plant antiviral defense

RNA silencing is a major mechanism of constitutive antiviral defense in plants, mediated by a number of proteins, including the Dicer-like (DCL) and Argonaute (AGO) endoribonucleases. Both DCL and AGO protein families comprise multiple members. In particular, the AGO protein family has expanded considerably in different plant lineages, with different family members having specialized functions. Although the general mode of action of AGO proteins is well established, the properties that make different AGO proteins more or less efficient at targeting viruses are less well understood. In this report, we review methodologies used to study AGO antiviral activity and current knowledge about which AGO family members are involved in antiviral defense. In addition, we discuss what is known about the different properties of AGO proteins thought to be associated with this function.

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.

Mapping Quantitative Trait Loci for Soybean Seedling Shoot and Root Architecture Traits in an Inter-Specific Genetic Population

Wild soybean species (Glycine soja Siebold & Zucc.) comprise a unique resource to widen the genetic base of cultivated soybean [Glycine max (L.) Merr.] for various agronomic traits. An inter-specific mapping population derived from a cross of cultivar Williams 82 and PI 483460B, a wild soybean accession, was utilized for genetic characterization of root architecture traits. The objectives of this study were to identify and characterize quantitative trait loci (QTL) for seedling shoot and root architecture traits, as well as to determine additive/epistatic interaction effects of identified QTLs. A total of 16,469 single nucleotide polymorphisms (SNPs) developed for the Illumina beadchip genotyping platform were used to construct a high resolution genetic linkage map. Among the 11 putative QTLs identified, two significant QTLs on chromosome 7 were determined to be associated with total root length (RL) and root surface area (RSA) with favorable alleles from the wild soybean parent. These seedling root traits, RL (BARC_020495_04641 ~ BARC_023101_03769) and RSA (SNP02285 ~ SNP18129_Magellan), could be potential targets for introgression into cultivated soybean background to improve both tap and lateral roots. The RL QTL region harbors four candidate genes with higher expression in root tissues: Phosphofructokinase (Glyma.07g126400), Snf7 protein (Glyma.07g127300), unknown functional gene (Glyma.07g127900), and Leucine Rich-Repeat protein (Glyma.07g127100). The novel alleles inherited from the wild soybean accession could be used as molecular markers to improve root system architecture and productivity in elite soybean lines.

ARGONAUTE5 Represses Age-Dependent Induction of Flowering through Physical and Functional Interaction with miR156 in Arabidopsis

Flowering time is a finely tuned process in plants, in part controlled by the age-regulated microRNA156 (miR156), which functions by suppressing the transcripts of SQUAMOSA-PROMOTER BINDING LIKE (SPL) transcription factors. ARGONAUTE (AGO) proteins are essential effectors of miRNA-mediated gene regulation. However, which AGO(s) mediate(s) the control of flowering time remains unclear. Here, we demonstrate a role of AGO5 in controlling flowering time by modulating the expression of SPL transcription factors. We show that AGO5 interacts physically and functionally with miR156 and that ago5 mutants present an early flowering phenotype in Arabidopsis. Furthermore, in ago5 mutants, the repression of flowering caused by miR156 overexpression is largely reversed, whereas leaf morphology remains unaffected. Our results thus indicate a specific role for AGO5 in mediating miR156 activity in meristematic, but not vegetative, tissue. As such, our data suggest a spatiotemporal regulation of the miR156 aging pathway mediated through different AGO proteins in different tissues.

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.

Extensive profiling in Arabidopsis reveals abundant polysome-associated 24-nt small RNAs including AGO5-dependent pseudogene-derived siRNAs

In a reductionist perspective, plant silencing small (s)RNAs are often classified as mediating nuclear transcriptional gene silencing (TGS) or cytosolic posttranscriptional gene silencing (PTGS). Among the PTGS diagnostics is the association of AGOs and their sRNA cargos with the translation apparatus. In Arabidopsis, this is observed for AGO1 loaded with micro(mi)RNAs and, accordingly, translational-repression (TR) is one layer of plant miRNA action. Using AGO1:miRNA-mediated TR as a paradigm, we explored, with two unrelated polysome-isolation methods, which, among the ten Arabidopsis AGOs and numerous sRNA classes, interact with translation. We found that representatives of all three AGO-clades associate with polysomes, including the TGS-effector AGO4 and stereotypical 24-nt sRNAs that normally mediate TGS of transposons/repeats. Strikingly, approximately half of these annotated 24-nt siRNAs displayed unique matches in coding regions/introns of genes, and in pseudogenes, but not in transposons/repeats commonly found in their vicinity. Protein-coding gene-derived 24-nt sRNAs correlate with gene-body methylation. Those derived from pseudogenes belong to two main clusters defined by their parental-gene expression patterns, and are vastly enriched in AGO5, itself found on polysomes. Based on their tight expression pattern in developing and mature siliques, their biogenesis, and genomic/epigenomic features of their loci-of-origin, we discuss potential roles for these hitherto unknown polysome-enriched, pseudogene-derived siRNAs.

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.

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.

Genotyping of Soybean Cultivars With Medium-Density Array Reveals the Population Structure and QTNs Underlying Maturity and Seed Traits

Soybean was domesticated about 5,000 to 6,000 years ago in China. Although genotyping technologies such as genotyping by sequencing (GBS) and high-density array are available, it is convenient and economical to genotype cultivars or populations using medium-density SNP array in genetic study as well as in molecular breeding. In this study, 235 cultivars, collected from China, Japan, USA, Canada and some other countries, were genotyped using SoySNP8k iSelect BeadChip with 7,189 single nucleotide polymorphisms (SNPs). In total, 4,471 polymorphic SNP markers were used to analyze population structure and perform genome-wide association study (GWAS). The most likely K value was 7, indicating this population can be divided into 7 subpopulations, which is well in accordance with the geographic origins of cultivars or accession studied. The LD decay rate was estimated at 184 kb, where r2 dropped to half of its maximum value (0.205). GWAS using FarmCPU detected a stable quantitative trait nucleotide (QTN) for hilum color and seed color, which is consistent with the known loci or genes. Although no universal QTNs for flowering time and maturity were identified across all environments, a total of 30 consistent QTNs were detected for flowering time (R1) or maturity (R7 and R8) on 16 chromosomes, most of them were corresponding to known E1 to E4 genes or QTL region reported in SoyBase (soybase.org). Of 16 consistent QTNs for protein and oil contents, 11 QTNs were detected having antagonistic effects on protein and oil content, while 4 QTNs soly for oil content, and one QTN soly for protein content. The information gained in this study demonstrated that the usefulness of the medium-density SNP array in genotyping for genetic study and molecular breeding.

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.

Recent advances in flower color variation and patterning of Japanese morning glory and petunia

The Japanese morning glory (Ipomoea nil) and petunia (Petunia hybrida), locally called "Asagao" and "Tsukubane-asagao", respectively, are popular garden plants. They have been utilized as model plants for studying the genetic basis of floricultural traits, especially anthocyanin pigmentation in flower petals. In their long history of genetic studies, many mutations affecting flower pigmentation have been characterized, and both structural and regulatory genes for the anthocyanin biosynthesis pathway have been identified. In this review, we will summarize recent advances in the understanding of flower pigmentation in the two species with respect to flower hue and color patterning. Regarding flower hue, we will describe a novel enhancer of flavonoid production that controls the intensity of flower pigmentation, new aspects related to a flavonoid glucosyltransferase that has been known for a long time, and the regulatory mechanisms of vacuolar pH being a key determinant of red and blue coloration. On color patterning, we describe particular flower patterns regulated by epigenetic and RNA-silencing mechanisms. As high-quality whole genome sequences of the Japanese morning glory and petunia wild parents (P. axillaris and P. inflata, respectively) were published in 2016, further study on flower pigmentation will be accelerated.

Gene Silencing of Argonaute5 Negatively Affects the Establishment of the Legume-Rhizobia Symbiosis

The establishment of the symbiosis between legumes and nitrogen-fixing rhizobia is finely regulated at the transcriptional, posttranscriptional and posttranslational levels. Argonaute5 (AGO5), a protein involved in RNA silencing, can bind both viral RNAs and microRNAs to control plant-microbe interactions and plant physiology. For instance, AGO5 regulates the systemic resistance of Arabidopsis against Potato Virus X as well as the pigmentation of soybean (Glycine max) seeds. Here, we show that AGO5 is also playing a central role in legume nodulation based on its preferential expression in common bean (Phaseolus vulgaris) and soybean roots and nodules. We also report that the expression of AGO5 is induced after 1 h of inoculation with rhizobia. Down-regulation of AGO5 gene in P. vulgaris and G. max causes diminished root hair curling, reduces nodule formation and interferes with the induction of three critical symbiotic genes: Nuclear Factor Y-B (NF-YB), Nodule Inception (NIN) and Flotillin2 (FLOT2). Our findings provide evidence that the common bean and soybean AGO5 genes play an essential role in the establishment of the symbiosis with rhizobia.