Loss-of-function mutations affecting a specific Glycine max R2R3 MYB transcription factor result in brown hilum and brown seed coats

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

Natural and artificial selection of multiple alleles revealed through genomic analyses

Genome-to-phenome research in agriculture aims to improve crops through in silico predictions. Genome-wide association study (GWAS) is potent in identifying genomic loci that underlie important traits. As a statistical method, increasing the sample quantity, data quality, or diversity of the GWAS dataset positively impacts GWAS power. For more precise breeding, concrete candidate genes with exact functional variants must be discovered. Many post-GWAS methods have been developed to narrow down the associated genomic regions and, ideally, to predict candidate genes and causative mutations (CMs). Historical natural selection and breeding-related artificial selection both act to change the frequencies of different alleles of genes that control phenotypes. With higher diversity and more extensive GWAS datasets, there is an increased chance of multiple alleles with independent CMs in a single causal gene. This can be caused by the presence of samples from geographically isolated regions that arose during natural or artificial selection. This simple fact is a complicating factor in GWAS-driven discoveries. Currently, none of the existing association methods address this issue and need to identify multiple alleles and, more specifically, the actual CMs. Therefore, we developed a tool that computes a score for a combination of variant positions in a single candidate gene and, based on the highest score, identifies the best number and combination of CMs. The tool is publicly available as a Python package on GitHub, and we further created a web-based Multiple Alleles discovery (MADis) tool that supports soybean and is hosted in SoyKB (https://soykb.org/SoybeanMADisTool/). We tested and validated the algorithm and presented the utilization of MADis in a pod pigmentation L1 gene case study with multiple CMs from natural or artificial selection. Finally, we identified a candidate gene for the pod color L2 locus and predicted the existence of multiple alleles that potentially cause loss of pod pigmentation. In this work, we show how a genomic analysis can be employed to explore the natural and artificial selection of multiple alleles and, thus, improve and accelerate crop breeding in agriculture.

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.

The domestication-associated L1 gene encodes a eucomic acid synthase pleiotropically modulating pod pigmentation and shattering in soybean

Pod coloration is a domestication-related trait in soybean, with modern cultivars typically displaying brown or tan pods, while their wild relative, Glycine soja, possesses black pods. However, the factors regulating this color variation remain unknown. In this study, we cloned and characterized L1, the classical locus responsible for black pods in soybean. By using map-based cloning and genetic analyses, we identified the causal gene of L1 and revealed that it encodes a hydroxymethylglutaryl-coenzyme A (CoA) lyase-like (HMGL-like) domain protein. Biochemical assays showed that L1 functions as a eucomic acid synthase and facilitates the synthesis of eucomic acid and piscidic acid, both of which contribute to coloration of pods and seed coats in soybean. Interestingly, we found that L1 plants are more prone to pod shattering under light exposure than l1 null mutants because dark pigmentation increases photothermal efficiency. Hence, pleiotropic effects of L1 on pod color and shattering, as well as seed pigmentation, likely contributed to the preference for l1 alleles during soybean domestication and improvement. Collectively, our study provides new insights into the mechanism of pod coloration and identifies a new target for future de novo domestication of legume crops.

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.

The Allele Catalog Tool: a web-based interactive tool for allele discovery and analysis

BACKGROUND

The advancement of sequencing technologies today has made a plethora of whole-genome re-sequenced (WGRS) data publicly available. However, research utilizing the WGRS data without further configuration is nearly impossible. To solve this problem, our research group has developed an interactive Allele Catalog Tool to enable researchers to explore the coding region allelic variation present in over 1,000 re-sequenced accessions each for soybean, Arabidopsis, and maize.

RESULTS

The Allele Catalog Tool was designed originally with soybean genomic data and resources. The Allele Catalog datasets were generated using our variant calling pipeline (SnakyVC) and the Allele Catalog pipeline (AlleleCatalog). The variant calling pipeline is developed to parallelly process raw sequencing reads to generate the Variant Call Format (VCF) files, and the Allele Catalog pipeline takes VCF files to perform imputations, functional effect predictions, and assemble alleles for each gene to generate curated Allele Catalog datasets. Both pipelines were utilized to generate the data panels (VCF files and Allele Catalog files) in which the accessions of the WGRS datasets were collected from various sources, currently representing over 1,000 diverse accessions for soybean, Arabidopsis, and maize individually. The main features of the Allele Catalog Tool include data query, visualization of results, categorical filtering, and download functions. Queries are performed from user input, and results are a tabular format of summary results by categorical description and genotype results of the alleles for each gene. The categorical information is specific to each species; additionally, available detailed meta-information is provided in modal popups. The genotypic information contains the variant positions, reference or alternate genotypes, the functional effect classes, and the amino-acid changes of each accession. Besides that, the results can also be downloaded for other research purposes.

CONCLUSIONS

The Allele Catalog Tool is a web-based tool that currently supports three species: soybean, Arabidopsis, and maize. The Soybean Allele Catalog Tool is hosted on the SoyKB website ( https://soykb.org/SoybeanAlleleCatalogTool/ ), while the Allele Catalog Tool for Arabidopsis and maize is hosted on the KBCommons website ( https://kbcommons.org/system/tools/AlleleCatalogTool/Zmays and https://kbcommons.org/system/tools/AlleleCatalogTool/Athaliana ). Researchers can use this tool to connect variant alleles of genes with meta-information of species.

Genetic mapping of QTLs controlling brown seed coat traits by genome resequencing in sesame (Sesamum indicum L.)

INTRODUCTION

Sesame seeds have become an irreplaceable source of edible oils and food products with rich nutrients and a unique flavor, and their metabolite contents and physiological functions vary widely across different seed coat colors. Although the quantitative trait loci (QTLs) for genetic variation in seed coat color have been extensively investigated, the identification of unique genetic loci for intermediate colors such as brown has not been reported due to their complexity.

METHODS

Here, we crossed the white sesame 'Yuzhi No. 8' (YZ8) and the brown sesame 'Yanzhou Erhongpi' (YZEHP) to construct a recombinant inbred line (RIL) population with consecutive self-fertilization for ten generations.

RESULTS

The selfed F1 seeds were brown which was controlled by a dominant gene. Based on the genotyping by whole-genome resequencing of the RILs, a major-effect QTL for brown coat color was identified through both bulk segregant analysis (BSA) and genetic linkage mapping in sesame, which was located within a 1.19 Mb interval on chromosome 6 (qBSCchr6). Moreover, we found that the YZEHP seed coat initially became pigmented at 20 days post-anthesis (DPA) and was substantially colored at 30 DPA. We screened 13 possible candidate genes based on the effects of genetic variants on protein coding and predicted gene functions. Furthermore, qRT‒PCR was used to verify the expression patterns of these genes in different post-anthesis developmental periods. We noted that in comparison to YZ8 seeds, YZEHP seeds had expression of SIN_1023239 that was significantly up-regulated 2.5-, 9.41-, 6.0-, and 5.9-fold at 15, 20, 25, and 30 DPA, respectively, which was consistent with the pattern of brown seed coat pigment accumulation.

DISCUSSION

This study identified the first major-effect QTL for the control of the brown seed coat trait in sesame. This finding lays the foundation for further fine mapping and cloning as well as investigating the regulatory mechanism of seed coat color in sesame.

NM, Jebakani K. S, Selvaraj S, Pramitha J. L, Selvaraj R, Petchiammal K. I, Kather Sheriff S, Thinakaran J, Rathinamoorthy S, Kumar P. R. Genetic manipulation of anti-nutritional factors in major crops for a sustainable diet in future

The consumption of healthy food, in order to strengthen the immune system, is now a major focus of people worldwide and is essential to tackle the emerging pandemic concerns. Moreover, research in this area paves the way for diversification of human diets by incorporating underutilized crops which are highly nutritious and climate-resilient in nature. However, although the consumption of healthy foods increases nutritional uptake, the bioavailability of nutrients and their absorption from foods also play an essential role in curbing malnutrition in developing countries. This has led to a focus on anti-nutrients that interfere with the digestion and absorption of nutrients and proteins from foods. Anti-nutritional factors in crops, such as phytic acid, gossypol, goitrogens, glucosinolates, lectins, oxalic acid, saponins, raffinose, tannins, enzyme inhibitors, alkaloids, β-N-oxalyl amino alanine (BOAA), and hydrogen cyanide (HCN), are synthesized in crop metabolic pathways and are interconnected with other essential growth regulation factors. Hence, breeding with the aim of completely eliminating anti-nutrition factors tends to compromise desirable features such as yield and seed size. However, advanced techniques, such as integrated multi-omics, RNAi, gene editing, and genomics-assisted breeding, aim to breed crops in which negative traits are minimized and to provide new strategies to handle these traits in crop improvement programs. There is also a need to emphasize individual crop-based approaches in upcoming research programs to achieve smart foods with minimum constraints in future. This review focuses on progress in molecular breeding and prospects for additional approaches to improve nutrient bioavailability in major crops.

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.

The mungbean VrP locus encoding MYB90, an R2R3-type MYB protein, regulates anthocyanin biosynthesis

Anthocyanins are water-soluble pigments present in several tissues/parts of plants. The pigments provide color and are wildly known for health benefits for human, insect attraction for plant pollination, and stress resistance in plants. Anthocyanin content variations in mungbean [Vigna radiata (L.) Wilczek] were first noticed a long time ago, but the genetic mechanism controlling the anthocyanins in mungbean remains unknown. An F₂ population derived from the cross between purple-hypocotyl (V2709) and green-hypocotyl (Sulv1) mungbeans was used to map the VrP locus controlling purple hypocotyl. The VrP locus was mapped to a 78.9-kb region on chromosome 4. Sequence comparison and gene expression analysis identified an R2R3-MYB gene VrMYB90 as the candidate gene for the VrP locus. Haplotype analysis using 124 mungbean accessions suggested that 10 single nucleotide polymorphisms (SNPs) in exon 3 may lead to an abolished expression of VrMYB90 and an absence of anthocyanin accumulation in the hypocotyl of Sulv1 and KPS2. The overexpression of VrMYB90 in mungbean hairy root, tobacco leaf, and Arabidopsis resulted in anthocyanin accumulation (purple color). Gene expression analysis demonstrated that VrMYB90 regulated anthocyanin accumulation in the hypocotyl, stem, petiole, and flowers, and the expression was sensitive to light. VrMYB90 protein may upregulate VrDFR encoding dihydroflavonol 4-reductase at the late biosynthesis step of anthocyanins in mungbeans. These results suggest that VrMYB90 is the dominator in the spatiotemporal regulation of anthocyanin biosynthesis. Our results provide insight into the biosynthesis mechanism of anthocyanin and a theoretical basis for breeding mungbeans.

A novel Synthetic phenotype association study approach reveals the landscape of association for genomic variants and phenotypes

INTRODUCTION

Genome-Wide Association Studies (GWAS) identify tagging variants in the genome that are statistically associated with the phenotype because of their linkage disequilibrium (LD) relationship with the causative mutation (CM). When both low-density genotyped accession panels with phenotypes and resequenced data accession panels are available, tagging variants can assist with post-GWAS challenges in CM discovery.

OBJECTIVES

Our objective was to identify additional GWAS evaluation criteria to assess correspondence between genomic variants and phenotypes, as well as enable deeper analysis of the localized landscape of association.

METHODS

We used genomic variant positions as Synthetic phenotypes in GWAS that we named "Synthetic phenotype association study" (SPAS). The extreme case of SPAS is what we call an "Inverse GWAS" where we used CM positions of cloned soybean genes. We developed and validated the Accuracy concept as a measure of the correspondence between variant positions and phenotypes.

RESULTS

The SPAS approach demonstrated that the genotype status of an associated variant used as a Synthetic phenotype enabled us to explore the relationships between tagging variants and CMs, and further, that utilizing CMs as Synthetic phenotypes in Inverse GWAS illuminated the landscape of association. We implemented the Accuracy calculation for a curated accession panel to an online Accuracy calculation tool (AccuTool) as a resource for gene identification in soybean. We demonstrated our concepts on three examples of soybean cloned genes. As a result of our findings, we devised an enhanced "GWAS to Genes" analysis (Synthetic phenotype to CM strategy, SP2CM). Using SP2CM, we identified a CM for a novel gene.

CONCLUSION

The SP2CM strategy utilizing Synthetic phenotypes and the Accuracy calculation of correspondence provides crucial information to assist researchers in CM discovery. The impact of this work is a more effective evaluation of landscapes of GWAS associations.

Fine-mapping and identification of a candidate gene controlling seed coat color in melon (Cucumis melo L. var. chinensis Pangalo)

MELO3C019554 encoding a homeobox protein (PHD transcription factor) is a candidate gene that involved in the formation of seed coat color in melon. Seed coat color is related to flavonoid content which is closely related to seed dormancy. According to the genetic analysis of a six-generation population derived from two parents (IC2508 with a yellow seed coat and IC2518 with a brown seed coat), we discovered that the yellow seed coat trait in melon is controlled by a single dominant gene, named CmBS-1. Bulked segregant analysis sequencing (BSA-Seq) revealed that the gene is located at 11,860,000-15,890,000 bp (4.03 Mb) on Chr 6. The F₂ population was genotyped using insertion-deletions (InDels), from which cleaved amplified polymorphic sequence (dCAPS) markers were derived to construct a genetic map. The gene was then fine-mapped to a 233.98 kb region containing 12 genes. Based on gene sequence analysis with two parents, we found that the MELO3C019554 gene encoding a homeobox protein (PHD transcription factor) had a nonsynonymous single nucleotide polymorphism (SNP) mutation in the coding sequence (CDS), and the SNP mutation resulted in the conversion of an amino acid (A → T) at residue 534. In addition, MELO3C019554 exhibited lower relative expression levels in the yellow seed coat than in the brown seed coat. Furthermore, we found that MELO3C019554 is related to 12 flavonoid metabolites. Thus, we predicted that MELO3C019554 is a candidate gene controlling seed coat color in melon. The study lays a foundation for further cloning projects and functional analysis of this gene, as well as marker-assisted selection breeding.

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.

A pubescence color gene enhances tolerance to cold-induced seed cracking in yellow soybean

In yellow soybean, severe cold weather causes seed cracking on the dorsal side. Yellow soybeans carry the I or iⁱ allele of the I locus and have a yellow (I) or pigmented (iⁱ ) hilum. We previously isolated an additional allele, designated as Ic, of the I locus, and reported that yellow soybeans with the IcIc genotype may be tolerant to cold-induced seed cracking. The Ic allele by itself, however, does not confer high tolerance. The association of a pubescence color gene (T) with suppression of low-temperature-induced seed coat deterioration has been previously reported. In the present study, we tested whether T is effective for the suppression of cold-induced seed cracking using two pairs of near-isogenic lines for the T locus in the iⁱiⁱ or IcIc background. In both backgrounds, the cracked seed rate of the near-isogenic line with the TT genotype was significantly lower than that with the tt genotype, which indicates that T has an inhibitory effect on cold-induced seed cracking. Furthermore, we also showed that gene pyramiding of Ic and T can improve tolerance to cold-induced seed cracking. Our findings should aid the development of highly SC-tolerant cultivars in soybean breeding programs.

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.

Domestication of Crop Metabolomes: Desired and Unintended Consequences

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

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.

Overexpression of GmGAMYB Accelerates the Transition to Flowering and Increases Plant Height in Soybean

The flowering time and plant height of soybean are important agronomic characters, which control the adaptability and yield of soybean. R2R3 MYB transcription factor plays an important regulatory role in plant growth and development. In this study, soybean GmGAMYB gene of R2R3-MYB type was induced by long-days (LDs). GmGAMYB showed higher transcriptional levels in the flowers, leaves and pods of soybean. Overexpression of GmGAMYB in transgenic soybean showed earlier flowering time and maturity in LDs and short-days (SDs). GmGAMYB interacted with GmGBP1 and might promote flowering time by up-regulating the expression of GmFULc gene in soybean. Moreover, the expression level of GmGAMYB was also induced by gibberellins (GAs) and the plant height of GmGAMYB-ox plants was significantly increased, which was caused by the enlargement of internode cell in stem. Furthermore, GmGAMYB overexpression led to increased GA sensitivity in the hypocotyl of soybean seedlings compared with WT. GmGAMYB may be a positive regulator of GA response of promoting plant height by up-regulating the expression of GmGA20ox gene in soybean. Together, our studies preliminarily showed that the partial functions of GmGAMYB in regulating flowering time and GA pathway.

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.

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.

Characterization of the soybean R2R3-MYB transcription factor GmMYB81 and its functional roles under abiotic stresses

R2R3-type MYBs are a key group of regulatory factors that control diverse developmental processes and stress tolerance in plants. Soybean is a major legume crop with the richness of seed protein and edible vegetable oil, and 244 R2R3-type MYBs have been identified in soybean. However, the knowledge regarding their functional roles has been greatly limited as yet. In this study, a novel R2R3-type MYB (GmMYB81) was functionally characterized in soybean, and it is closely related to two abiotic stress-associated regulators (AtMYB44 and AtMYB77). GmMYB81 transcripts not only differentially accumulated in soybean tissues and during embryo development, but also were significantly enhanced by drought, salt and cold stress. Histochemical GUS assay in Arabidopsis indicated that GmMYB81 promoter showed high activity in seedlings, rosette leaves, inflorescences, silique wall, mature anthers, roots, and germinating seeds. Further investigation indicated that over-expression of GmMYB81 in Arabidopsis caused auxin-associated phenotypes, including small flower and silique, more branch, and weakened apical dominance. Moreover, over-expression of GmMYB81 significantly elevated the rates of seed germination and green seedling under salt and drought stress, indicating that GmMYB81 might confer plant tolerance to salt and drought stress during seed germination. Additionally, protein interaction analysis showed that GmMYB81 interacts with the abiotic stress regulator GmSGF14l. Further observation indicated that they displayed similar expression patterns under drought and salt stress, suggesting GmMYB81 and GmSGF14l might cooperatively affect stress tolerance. These findings will facilitate future investigations of the regulatory mechanisms of GmMYB81 in response to plant stress tolerance, especially seed germination under abiotic stresses.

Integrated Metabolomics and Transcriptomics Analyses Reveal the Molecular Mechanisms Underlying the Accumulation of Anthocyanins and Other Flavonoids in Cowpea Pod (Vigna unguiculata L.)

As an important vegetable crop of the legume family, cowpea (Vigna unguiculata L.) is grown widely for its tender pod with good taste and nutrition. The purple cowpea pods attract more attention mainly for the eye-catching color and health-promoting ingredients. Initially, large quantities of two major anthocyanins (delphinidin 3-O-glucoside and cyanidin 3-O-glucoside) and nine kinds of flavonoids (most are quercetin-based flavonol glycosides) were separated and identified from purple cowpea pod by ultra-high performance liquid chromatography coupled with quadrupole Orbitrap high-resolution mass spectrometry. To study them systematically, two representative cowpea cultivars with a drastic difference in anthocyanin accumulation were further analyzed by the integration of metabolomics and transcriptomics. A total of 56 differentially accumulated metabolites and 4142 differentially expressed genes were identified, respectively. On the basis of the comprehensive analysis of multiomic data, it was shown that VuMYB90-1, VuMYB90-2, VuMYB90-3, VuCPC, VuMYB4, and endogenous bHLH and WD40 proteins coordinately control anthocyanin and flavonoid accumulation via transcriptional regulation of structural genes in purple cowpea pod.

Development of a minigenome cassette for Lettuce necrotic yellows virus: A first step in rescuing a plant cytorhabdovirus

Rhabdoviruses are enveloped negative-sense RNA viruses that have numerous biotechnological applications. However, recovering plant rhabdoviruses from cDNA remains difficult due to technical difficulties such as the need for concurrent in planta expression of the viral genome together with the viral nucleoprotein (N), phosphoprotein (P) and RNA-dependent RNA polymerase (L) and viral genome instability in E. coli. Here, we developed a negative-sense minigenome cassette for Lettuce necrotic yellows virus (LNYV). We introduced introns into the unstable viral ORF and employed Agrobacterium tumefaciens to co-infiltrate Nicotiana with the genes for the N, P, and L proteins together with the minigenome cassette. The minigenome cassette included the Discosoma sp. red fluorescent protein gene (DsRed) cloned in the negative-sense between the viral trailer and leader sequences which were placed between hammerhead and hepatitis delta ribozymes. In planta DsRed expression was demonstrated by western blotting while the appropriate splicing of introduced introns was confirmed by sequencing of RT-PCR product.

Loss-of-Function Mutation of Soybean R2R3 MYB Transcription Factor Dilutes Tawny Pubescence Color

Pubescence color of soybean is controlled by two genes, T and Td. In the presence of a dominant T allele, dominant and recessive alleles of the Td locus generate tawny and light tawny (or near-gray) pubescence, respectively. Flavones, responsible for pubescence color, are synthesized via two copies of flavone synthase II genes (FNS II-1 and FNS II-2). This study was conducted to map and clone the Td gene. Genetic and linkage analysis using an F2 population and F3 families derived from a cross between a Clark near-isogenic line with light tawny pubescence (genotype: TT tdtd) and a Harosoy near-isogenic line with tawny pubescence (TT TdTd) revealed a single gene for pubescence color around the end of chromosome 3. Genome sequence alignment of plant introductions revealed an association between premature stop codons in Glyma.03G258700 (R2R3 MYB transcription factor) and recessive td allele. Cultivars and lines having near-gray or light tawny pubescence and a gray pubescence cultivar with td allele had premature stop codons in the gene. These results suggest that Glyma.03G258700 corresponds to the Td gene. It was predominantly expressed in pubescence. Compared to a tawny pubescence line, a near-isogenic line with td allele produced extremely small amounts of transcripts of Glyma.03G258700, FNS II-1, and FNS II-2 in pubescence. The promoter of FNS II-1 and FNS II-2 shared cis-acting regulatory elements for binding of MYB proteins. These results suggest that the wild type of Glyma.03G258700 protein may bind to the promoter of FNS II genes and upregulate their expression, resulting in increased flavone content and deeper pubescence color. In contrast, mutated Glyma.03G258700 protein may fail to upregulate the expression of FNS II genes, resulting in decreased flavone content and dilute pubescence color.

Parallelism and convergence in post-domestication adaptation in cereal grasses

The selection of desirable traits in crops during domestication has been well studied. Many crops share a suite of modified phenotypic characteristics collectively known as the domestication syndrome. In this sense, crops have convergently evolved. Previous work has demonstrated that, at least in some instances, convergence for domestication traits has been achieved through parallel molecular means. However, both demography and selection during domestication may have placed limits on evolutionary potential and reduced opportunities for convergent adaptation during post-domestication migration to new environments. Here we review current knowledge regarding trait convergence in the cereal grasses and consider whether the complexity and dynamism of cereal genomes (e.g., transposable elements, polyploidy, genome size) helped these species overcome potential limitations owing to domestication and achieve broad subsequent adaptation, in many cases through parallel means. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Genetic diversity patterns and domestication origin of soybean

KEY MESSAGE

Genotyping data of a comprehensive Korean soybean collection obtained using a large SNP array were used to clarify global distribution patterns of soybean and address the evolutionary history of soybean. Understanding diversity and evolution of a crop is an essential step to implement a strategy to expand its germplasm base for crop improvement research. Accessions intensively collected from Korea, which is a small but central region in the distribution geography of soybean, were genotyped to provide sufficient data to underpin population genetic questions. After removing natural hybrids and duplicated or redundant accessions, we obtained a non-redundant set comprising 1957 domesticated and 1079 wild accessions to perform population structure analyses. Our analysis demonstrates that while wild soybean germplasm will require additional sampling from diverse indigenous areas to expand the germplasm base, the current domesticated soybean germplasm is saturated in terms of genetic diversity. We then showed that our genome-wide polymorphism map enabled us to detect genetic loci underlying flower color, seed-coat color, and domestication syndrome. A representative soybean set consisting of 194 accessions was divided into one domesticated subpopulation and four wild subpopulations that could be traced back to their geographic collection areas. Population genomics analyses suggested that the monophyletic group of domesticated soybeans was likely originated at a Japanese region. The results were further substantiated by a phylogenetic tree constructed from domestication-associated single nucleotide polymorphisms identified in this study.

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.

QTL sequencing strategy to map genomic regions associated with resistance to ascochyta blight in chickpea

Whole-genome sequencing-based bulked segregant analysis (BSA) for mapping quantitative trait loci (QTL) provides an efficient alternative approach to conventional QTL analysis as it significantly reduces the scale and cost of analysis with comparable power to QTL detection using full mapping population. We tested the application of next-generation sequencing (NGS)-based BSA approach for mapping QTLs for ascochyta blight resistance in chickpea using two recombinant inbred line populations CPR-01 and CPR-02. Eleven QTLs in CPR-01 and six QTLs in CPR-02 populations were mapped on chromosomes Ca1, Ca2, Ca4, Ca6 and Ca7. The QTLs identified in CPR-01 using conventional biparental mapping approach were used to compare the efficiency of NGS-based BSA in detecting QTLs for ascochyta blight resistance. The QTLs on chromosomes Ca1, Ca4, Ca6 and Ca7 overlapped with the QTLs previously detected in CPR-01 using conventional QTL mapping method. The QTLs on chromosome Ca4 were detected in both populations and overlapped with the previously reported QTLs indicating conserved region for ascochyta blight resistance across different chickpea genotypes. Six candidate genes in the QTL regions identified using NGS-based BSA on chromosomes Ca2 and Ca4 were validated for their association with ascochyta blight resistance in the CPR-02 population. This study demonstrated the efficiency of NGS-based BSA as a rapid and cost-effective method to identify QTLs associated with ascochyta blight in chickpea.

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.

Transcriptomic variation in proanthocyanidin biosynthesis pathway genes in soybean (Glycine spp.)

BACKGROUND

Proanthocyanidins are oligomeric or polymeric end products of flavonoid metabolic pathways starting with the central phenylpropanoid pathway. Although soybean (Glycine spp.) seeds represent a major source of nutrients for the human diet, as well as components for the cosmetics industry as a result of their high levels of flavonoid metabolites, including isoflavonoids, anthocyanins and proanthocyanidins, the genetic regulatory mechanisms underlying proanthocyanidin biosynthesis in soybean remain unclear.

RESULTS

We evaluated interspecific and intraspecific variability in flavonoid components in soybean using 43 cultivars, landraces and wild soybean accessions. We performed transcriptomic profiling of genes encoding enzymes involved in flavonoid biosynthesis using three soybean genotypes, Hwangkeum (elite cultivar), IT109098 (landrace) and IT182932 (wild accession), in seeds. We identified a Glycine max landrace, IT109098, with a proanthocyanidin content as high as that of wild soybean. Different homologous genes for anthocyanidin reductase, which is involved in proanthocyanidin biosynthesis, were detected as differentially expressed genes between IT109098 and IT182932 compared to Hwangkeum.

CONCLUSION

We detected major differences in the transcriptional levels of genes involved in the biosynthesis of proanthocyanidin and anthocyanin among genotypes beginning at the early stage of seed development. The results of the present study provide insights into the underlying genetic variation in proanthocyanidin biosynthesis among soybean genotypes. © 2017 Society of Chemical Industry.

GmMYB181, a Soybean R2R3-MYB Protein, Increases Branch Number in Transgenic Arabidopsis

Soybean (Glycine max) is an important economic crop that provides abundant oil and high quality protein for human beings. As the process of reproductive growth directly determines the crop seed yield and quality, we initiated studies to identify genes that regulate soybean floral organ development. One R2R3-MYB transcription factor gene, designated as GmMYB181, was found to be enriched in flowers based on microarray analysis and was further functionally investigated in transgenic Arabidopsis. GmMYB181 protein contains two MYB domains, which localized to the nucleus and displayed transcriptional activation in yeast hybrid system. Real-time quantitative PCR (qRT-PCR) results suggested GmMYB181 exclusively expressed in flower tissue. In Arabidopsis, overexpression of GmMYB181 altered the morphology of floral organs, fruit size and plant architecture, including outward curly sepals, smaller siliques, increased lateral branches and reduced plant height, indicating that GmMYB181 is involved in the development of reproductive organs and plays an important role in controlling plant architecture. Further, microarray analysis revealed that overexpressing GmMYB181 in Arabidopsis affected the expression of 3450 genes in mature flowers, including those involved in floral organ, seed/fruit development, and responded to different hormone signals.

Whole-Genome Resequencing Identifies the Molecular Genetic Cause for the Absence of a Gy5 Glycinin Protein in Soybean PI 603408

During ongoing proteomic analysis of the soybean (Glycine max (L.) Merr) germplasm collection, PI 603408 was identified as a landrace whose seeds lack accumulation of one of the major seed storage glycinin protein subunits. Whole genomic resequencing was used to identify a two-base deletion affecting glycinin 5 The newly discovered deletion was confirmed to be causative through immunological, genetic, and proteomic analysis, and no significant differences in total seed protein content were found to be due to the glycinin 5 loss-of-function mutation per se In addition to focused studies on this one specific glycinin subunit-encoding gene, a total of 1,858,185 nucleotide variants were identified, of which 39,344 were predicted to affect protein coding regions. In order to semiautomate analysis of a large number of soybean gene variants, a new SIFT 4G (Sorting Intolerant From Tolerated 4 Genomes) database was designed to predict the impact of nonsynonymous single nucleotide soybean gene variants, potentially enabling more rapid analysis of soybean resequencing data in the future.

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.

Accumulation of proanthocyanidins and/or lignin deposition in buff-pigmented soybean seed coats may lead to frequent defective cracking

MAIN CONCLUSION

Defective cracking frequently occurs in buff-pigmented soybean seed coats, where proanthocyanidins accumulate and lignin is deposited, suggesting that proanthocyanidins and/or lignin may change physical properties and lead to defective cracking. In the seed production of many yellow soybean (Glycine max) cultivars, very low percentages of self-pigmented seeds are commonly found. This phenomenon is derived from a recessive mutation of the I gene inhibiting seed coat pigmentation. In Japan, most of these self-pigmented seeds are buff-colored, and frequently show multiple defective cracks in the seed coat. However, it is not known why cracking occurs specifically in buff seed coats. In this study, quantitative analysis was performed between yellow and buff soybean seed coats. Compared with yellow soybeans, in which defective cracking rarely occurs, contents of proanthocyanidins (PAs) and lignin were significantly higher in buff seed coats. Histochemical data of PAs and lignin in the seed coats strongly supported this result. Measurements of the physical properties of seed coats using a texture analyzer showed that a hardness value was significantly decreased in the buff seed coats. These results suggest that PA accumulation and/or lignin deposition may affect the physical properties of buff seed coats and lead to the defective cracking. This work contributes to understanding of the mechanism of defective cracking, which decreases the seed quality of soybean and related legumes.

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.

Anthocyanins and flavonols are responsible for purple color of Lablab purpureus (L.) sweet pods

Lablab pods, as dietary vegetable, have high nutritional values similar to most of edible legumes. Moreover, our studies confirmed that purple lablab pods contain the natural pigments of anthocyanins and flavonols. Compared to green pods, five kinds of anthocyanins (malvidin, delphinidin and petunidin derivatives) were found in purple pods by HPLC-ESI-MS/MS and the major contents were delphinidin derivatives. Besides, nine kinds of polyphenol derivatives (quercetin, myricetin, kaempferol and apigenin derivatives) were detected by UPLC-ESI-MS/MS and the major components were quercetin and myricetin derivatives. In order to discover their molecular mechanism, expression patterns of biosynthesis and regulatory gens of anthocyanins and flavonols were investigated. Experimental results showed that LpPAL, LpF3H, LpF3'H, LpDFR, LpANS and LpPAP1 expressions were significantly induced in purple pods compared to green ones. Meanwhile, transcripts of LpFLS were more abundant in purple pods than green or yellow ones, suggestind that co-pigments of anthocyanins and flavonols are accumulated in purple pods. Under continuously dark condition, no anthocyanin accumulation was detected in purple pods and transcripts of LpCHS, LpANS, LpFLS and LpPAP1 were remarkably repressed, indicating that anthocyanins and flavonols biosynthesis in purple pods was regulated in light-dependent manner. These results indicate that co-pigments of anthocyanins and flavonols contribute to purple pigmentations of pods.

Domestication footprints anchor genomic regions of agronomic importance in soybeans

Present-day soybeans consist of elite cultivars and landraces (Glycine max, fully domesticated (FD)), annual wild type (Glycine soja, nondomesticated (ND)), and semi-wild type (semi-domesticated (SD)). FD soybean originated in China, although the details of its domestication history remain obscure. More than 500 diverse soybean accessions were sequenced using specific-locus amplified fragment sequencing (SLAF-seq) to address fundamental questions regarding soybean domestication. In total, 64,141 single nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) > 0.05 were found among the 512 tested accessions. The results indicated that the SD group is not a hybrid between the FD and ND groups. The initial domestication region was pinpointed to central China (demarcated by the Great Wall to the north and the Qinling Mountains to the south). A total of 800 highly differentiated genetic regions and > 140 selective sweeps were identified, and these were three- and twofold more likely, respectively, to encompass a known quantitative trait locus (QTL) than the rest of the soybean genome. Forty-three potential quantitative trait nucleotides (QTNs; including 15 distinct traits) were identified by genome-wide association mapping. The results of the present study should be beneficial for soybean improvement and provide insight into the genetic architecture of traits of agronomic importance.

CACTA-superfamily transposable element is inserted in MYB transcription factor gene of soybean line producing variegated seeds

The R gene of soybean, presumably encoding a MYB transcription factor, controls seed coat color. The gene consists of multiple alleles, R (black), r-m (black spots and (or) concentric streaks on brown seed), and r (brown seed). This study was conducted to determine the structure of the MYB transcription factor gene in a near-isogenic line (NIL) having r-m allele. PCR amplification of a fragment of the candidate gene Glyma.09G235100 generated a fragment of about 1 kb in the soybean cultivar Clark, whereas a fragment of about 14 kb in addition to fragments of 1 and 1.4 kb were produced in L72-2040, a Clark 63 NIL with the r-m allele. Clark 63 is a NIL of Clark with the rxp and Rps1 alleles. A DNA fragment of 13 060 bp was inserted in the intron of Glyma.09G235100 in L72-2040. The fragment had the CACTA motif at both ends, imperfect terminal inverted repeats (TIR), inverse repetition of short sequence motifs close to the 5' and 3' ends, and a duplication of three nucleotides at the site of integration, indicating that it belongs to a CACTA-superfamily transposable element. We designated the element as Tgm11. Overall nucleotide sequence, motifs of TIR, and subterminal repeats were similar to those of Tgm1 and Tgs1, suggesting that these elements comprise a family.

Identification of loci governing eight agronomic traits using a GBS-GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping-by-sequencing (GBS) approach was used to provide dense genome-wide marker coverage (>47,000 SNPs) for a panel of 304 short-season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.

Methylation affects transposition and splicing of a large CACTA transposon from a MYB transcription factor regulating anthocyanin synthase genes in soybean seed coats

We determined the molecular basis of three soybean lines that vary in seed coat color at the R locus which is thought to encode a MYB transcription factor. RM55-r(m) is homozygous for a mutable allele (r(m)) that specifies black and brown striped seeds; RM30-R* is a stable black revertant isoline derived from the mutable line; and RM38-r has brown seed coats due to a recessive r allele shown to translate a truncated MYB protein. Using long range PCR, 454 sequencing of amplicons, and whole genome re-sequencing, we determined that the variegated RM55-r(m) line had a 13 kb CACTA subfamily transposon insertion (designated TgmR*) at a position 110 bp from the beginning of Intron2 of the R locus, Glyma09g36983. Although the MYB encoded by R was expressed at only very low levels in older seed coats of the black revertant RM30-R* line, it upregulated expression of anthocyanidin synthase genes (ANS2, ANS3) to promote the synthesis of anthocyanins. Surprisingly, the RM30-R* revertant also carried the 13 kb TgmR* insertion in Intron2. Using RNA-Seq, we showed that intron splicing was accurate, albeit at lower levels, despite the presence of the 13 kb TgmR* element. As determined by whole genome methylation sequencing, we demonstrate that the TgmR* sequence was relatively more methylated in RM30-R* than in the mutable RM55-r(m) progenitor line. The stabilized and more methylated RM30-R* revertant line apparently lacks effective binding of a transposae to its subterminal repeats, thus allowing intron splicing to proceed resulting in sufficient MYB protein to stimulate anthocyanin production and thus black seed coats. In this regard, the TgmR* element in soybean resembles McClintock's Spm-suppressible and change-of-state alleles of maize. This comparison explains the opposite effects of the TgmR* element on intron splicing of the MYB gene in which it resides depending on the methylation state of the element.

Methylation affects transposition and splicing of a large CACTA transposon from a MYB transcription factor regulating anthocyanin synthase genes in soybean seed coats

We determined the molecular basis of three soybean lines that vary in seed coat color at the R locus which is thought to encode a MYB transcription factor. RM55-r(m) is homozygous for a mutable allele (r(m)) that specifies black and brown striped seeds; RM30-R* is a stable black revertant isoline derived from the mutable line; and RM38-r has brown seed coats due to a recessive r allele shown to translate a truncated MYB protein. Using long range PCR, 454 sequencing of amplicons, and whole genome re-sequencing, we determined that the variegated RM55-r(m) line had a 13 kb CACTA subfamily transposon insertion (designated TgmR*) at a position 110 bp from the beginning of Intron2 of the R locus, Glyma09g36983. Although the MYB encoded by R was expressed at only very low levels in older seed coats of the black revertant RM30-R* line, it upregulated expression of anthocyanidin synthase genes (ANS2, ANS3) to promote the synthesis of anthocyanins. Surprisingly, the RM30-R* revertant also carried the 13 kb TgmR* insertion in Intron2. Using RNA-Seq, we showed that intron splicing was accurate, albeit at lower levels, despite the presence of the 13 kb TgmR* element. As determined by whole genome methylation sequencing, we demonstrate that the TgmR* sequence was relatively more methylated in RM30-R* than in the mutable RM55-r(m) progenitor line. The stabilized and more methylated RM30-R* revertant line apparently lacks effective binding of a transposae to its subterminal repeats, thus allowing intron splicing to proceed resulting in sufficient MYB protein to stimulate anthocyanin production and thus black seed coats. In this regard, the TgmR* element in soybean resembles McClintock's Spm-suppressible and change-of-state alleles of maize. This comparison explains the opposite effects of the TgmR* element on intron splicing of the MYB gene in which it resides depending on the methylation state of the element.

Variation block-based genomics method for crop plants

BACKGROUND

In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits.

RESULTS

We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color.

CONCLUSIONS

We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

RNA-Seq profiling of a defective seed coat mutation in Glycine max reveals differential expression of proline-rich and other cell wall protein transcripts

The plant cell wall performs a number of essential functions including providing shape to many different cell types and serving as a defense against potential pathogens. The net pattern mutation creates breaks in the seed coat of soybean (Glycine max) because of ruptured cell walls. Using RNA-Seq, we examined the seed coat transcriptome from three stages of immature seed development in two pairs of isolines with normal or defective seed coat phenotypes due to the net pattern. The genome-wide comparative study of the transcript profiles of these isolines revealed 364 differentially expressed genes in common between the two varieties that were further divided into different broad functional categories. Genes related to cell wall processes accounted for 19% of the differentially expressed genes in the middle developmental stage of 100-200 mg seed weight. Within this class, the cell wall proline-rich and glycine-rich protein genes were highly differentially expressed in both genetic backgrounds. Other genes that showed significant expression changes in each of the isoline pairs at the 100-200 mg seed weight stage were xylem serine proteinase, fasciclin-related genes, auxin and stress response related genes, TRANSPARENT TESTA 1 (TT1) and other transcription factors. The mutant appears to shift the timing of either the increase or decrease in the levels of some of the transcripts. The analysis of these data sets reveals the physiological changes that the seed coat undergoes during the formation of the breaks in the cell wall.

The role of the testa during development and in establishment of dormancy of the legume seed

Timing of seed germination is one of the key steps in plant life cycles. It determines the beginning of plant growth in natural or agricultural ecosystems. In the wild, many seeds exhibit dormancy and will only germinate after exposure to certain environmental conditions. In contrast, crop seeds germinate as soon as they are imbibed usually at planting time. These domestication-triggered changes represent adaptations to cultivation and human harvesting. Germination is one of the common sets of traits recorded in different crops and termed the "domestication syndrome." Moreover, legume seed imbibition has a crucial role in cooking properties. Different seed dormancy classes exist among plant species. Physical dormancy (often called hardseededness), as found in legumes, involves the development of a water-impermeable seed coat, caused by the presence of phenolics- and suberin-impregnated layers of palisade cells. The dormancy release mechanism primarily involves seed responses to temperature changes in the habitat, resulting in testa permeability to water. The underlying genetic controls in legumes have not been identified yet. However, positive correlation was shown between phenolics content (e.g., pigmentation), the requirement for oxidation and the activity of catechol oxidase in relation to pea seed dormancy, while epicatechin levels showed a significant positive correlation with soybean hardseededness. myeloblastosis family of transcription factors, WD40 proteins and enzymes of the anthocyanin biosynthesis pathway were involved in seed testa color in soybean, pea and Medicago, but were not tested directly in relation to seed dormancy. These phenolic compounds play important roles in defense against pathogens, as well as affecting the nutritional quality of products, and because of their health benefits, they are of industrial and medicinal interest. In this review, we discuss the role of the testa in mediating legume seed germination, with a focus on structural and chemical aspects.

A R2R3-MYB transcription factor, GmMYB12B2, affects the expression levels of flavonoid biosynthesis genes encoding key enzymes in transgenic Arabidopsis plants

Isoflavones play diverse roles in plant-microbe interactions and are potentially important for human nutrition and health. To study the regulation of isoflavonoid synthesis in soybean, the R2R3-MYB transcription factor GmMYB12B2 was isolated and characterized. Yeast expression experiments demonstrated that GmMYB12B2 showed transcriptional activity. GmMYB12B2 was localized in the nucleus when it was transiently expressed in onion epidermal cells. Real-time quantitative PCR analysis revealed that GmMYB12B2 transcription was increased in roots and mature seeds compared with other organs. The gene expression level in immature embryos was consistent with the accumulation of isoflavones. CHS8 is a key enzyme in plant flavonoid biosynthesis. Transient expression experiments in soybean calli demonstrated that CHS8 was regulated by GmMYB12B2 and produced more fluorescence. The expression levels of some key enzymes in flavonoid biosynthesis were examined in transgenic Arabidopsis lines. The results showed that the expression levels of PAL1, CHS and FLS in transgenic plants were significantly higher than those in wild type plants. However, the expression level of DFR was lower, and the expression levels of CHI, F3H and F3'H were the same in all lines. GmMYB12B2 expression caused a constitutive increase in the accumulation of flavonoids in transgenic Arabidopsis lines compared with wild type plants.

Molecular mechanisms involved in convergent crop domestication

Domestication has helped to understand evolution. We argue that, vice versa, novel insights into evolutionary principles could provide deeper insights into domestication. Molecular analyses have demonstrated that convergent phenotypic evolution is often based on molecular changes in orthologous genes or pathways. Recent studies have revealed that during plant domestication the causal mutations for convergent changes in key traits are likely to be located in particular genes. These insights may contribute to defining candidate genes for genetic improvement during the domestication of new plant species. Such efforts may help to increase the range of arable crops available, thus increasing crop biodiversity and food security to help meet the predicted demands of the continually growing global population under rapidly changing environmental conditions.

In silico comparison of genomic regions containing genes coding for enzymes and transcription factors for the phenylpropanoid pathway in Phaseolus vulgaris L. and Glycine max L. Merr

Legumes contain a variety of phytochemicals derived from the phenylpropanoid pathway that have important effects on human health as well as seed coat color, plant disease resistance and nodulation. However, the information about the genes involved in this important pathway is fragmentary in common bean (Phaseolus vulgaris L.). The objectives of this research were to isolate genes that function in and control the phenylpropanoid pathway in common bean, determine their genomic locations in silico in common bean and soybean, and analyze sequences of the 4CL gene family in two common bean genotypes. Sequences of phenylpropanoid pathway genes available for common bean or other plant species were aligned, and the conserved regions were used to design sequence-specific primers. The PCR products were cloned and sequenced and the gene sequences along with common bean gene-based (g) markers were BLASTed against the Glycine max v.1.0 genome and the P. vulgaris v.1.0 (Andean) early release genome. In addition, gene sequences were BLASTed against the OAC Rex (Mesoamerican) genome sequence assembly. In total, fragments of 46 structural and regulatory phenylpropanoid pathway genes were characterized in this way and placed in silico on common bean and soybean sequence maps. The maps contain over 250 common bean g and SSR (simple sequence repeat) markers and identify the positions of more than 60 additional phenylpropanoid pathway gene sequences, plus the putative locations of seed coat color genes. The majority of cloned phenylpropanoid pathway gene sequences were mapped to one location in the common bean genome but had two positions in soybean. The comparison of the genomic maps confirmed previous studies, which show that common bean and soybean share genomic regions, including those containing phenylpropanoid pathway gene sequences, with conserved synteny. Indels identified in the comparison of Andean and Mesoamerican common bean 4CL gene sequences might be used to develop inter-pool phenylpropanoid pathway gene-based markers. We anticipate that the information obtained by this study will simplify and accelerate selections of common bean with specific phenylpropanoid pathway alleles to increase the contents of beneficial phenylpropanoids in common bean and other legumes.