Genome-wide scan for oil quality reveals a coregulation mechanism of tocopherols and fatty acids in soybean seeds.

Tocopherols (vitamin E) play essential roles in human health because of their antioxidant activity, and plant-derived oils are the richest sources of tocopherols in the human diet. Although soybean (Glycine max) is one of the main sources of plant-derived oil and tocopherol in the world, the relationship between tocopherol and oil in soybean seeds remains unclear. Here, we focus on dissecting tocopherol metabolism with the long-term goal of increasing α-tocopherol content and soybean oil quality. We first collected tocopherol and fatty acid profiles in a soybean population (>800 soybean accessions) and found that tocopherol content increased during soybean domestication. A strong positive correlation between tocopherol and oil content was also detected. Five tocopherol pathway-related loci were identified using a metabolite genome-wide association study strategy. Genetic variations in three tocopherol pathway genes were responsible for total tocopherol content and composition in the soybean population through effects on enzyme activity, mainly caused by non-conserved amino acid substitution or changes in gene transcription level. Moreover, the fatty acid regulatory transcription factor GmZF351 directly activated tocopherol pathway gene expression, increasing both fatty acid and tocopherol contents in soybean seeds. Our study reveals the functional differentiation of tocopherol pathway genes in soybean populations and provides a framework for development of new soybean varieties with high α-tocopherol content and oil quality in seeds.

A route to de novo domestication of wild allotetraploid rice.

Cultivated rice varieties are all diploid, and polyploidization of rice has long been desired because of its advantages in genome buffering, vigorousness, and environmental robustness. However, a workable route remains elusive. Here, we describe a practical strategy, namely de novo domestication of wild allotetraploid rice. By screening allotetraploid wild rice inventory, we identified one genotype of Oryza alta (CCDD), polyploid rice 1 (PPR1), and established two important resources for its de novo domestication: (1) an efficient tissue culture, transformation, and genome editing system and (2) a high-quality genome assembly discriminated into two subgenomes of 12 chromosomes apiece. With these resources, we show that six agronomically important traits could be rapidly improved by editing O. alta homologs of the genes controlling these traits in diploid rice. Our results demonstrate the possibility that de novo domesticated allotetraploid rice can be developed into a new staple cereal to strengthen world food security.

Quantitative Trait Locus Analysis of Protein and Oil Content in Response to Planting Density in Soybean (Glycine max [L.] Merri.) Seeds Based on SNP Linkage Mapping.

Soybean varieties suitable for high planting density allow greater yields. However, the seed protein and oil contents, which determine the value of this crop, can be influenced by planting density. Thus, it is important to understand the genetic basis of the responses of different soybean genotypes to planting density. In this study, we quantified the protein and oil contents in a four-way recombinant inbred line (FW-RIL) soybean population under two planting densities and the response to density. We performed quantitative trait locus (QTL) mapping using a single nucleotide polymorphism (SNP) linkage map generated by inclusive composite interval mapping. We identified 14 QTLs for protein content and 17 for oil content at a planting density of 2.15 × 105 plant/ha (D1) and 14 QTLs for protein content and 20 for oil content at a planting density 3.0 × 105 plant/ha (D2). Among the QTLs detected, two oil-content QTLs was detected at both plant densities. In addition, we identified 38 QTLs for the responses of protein and oil contents to planting density. Of the QTLs detected, 70 were identified in previous studies, while 33 were newly identified. Fourty-five QTLs accounted for over 10% of the phenotypic variation of the corresponding trait, based on 23 QTLs at a marker interval distance of ~600 kb detected under different densities and with the responses to density difference. Pathway analysis revealed four candidate genes involved in protein and oil biosynthesis/metabolism. These results improve our understanding of the genetic underpinnings of protein and oil biosynthesis in soybean, laying the foundation for enhancing protein and oil contents and increasing yields in soybean.

Pan-Genome of Wild and Cultivated Soybeans.

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

Simultaneous changes in seed size, oil content and protein content driven by selection of SWEET homologues during soybean domestication.

Soybean accounts for more than half of the global production of oilseed and more than a quarter of the protein used globally for human food and animal feed. Soybean domestication involved parallel increases in seed size and oil content, and a concomitant decrease in protein content. However, science has not yet discovered whether these effects were due to selective pressure on a single gene or multiple genes. Here, re-sequencing data from >800 genotypes revealed a strong selection during soybean domestication on GmSWEET10a. The selection of GmSWEET10a conferred simultaneous increases in soybean-seed size and oil content as well as a reduction in the protein content. The result was validated using both near-isogenic lines carrying substitution of haplotype chromosomal segments and transgenic soybeans. Moreover, GmSWEET10b was found to be functionally redundant with its homologue GmSWEET10a and to be undergoing selection in current breeding, leading the the elite allele GmSWEET10b, a potential target for present-day soybean breeding. Both GmSWEET10a and GmSWEET10b were shown to transport sucrose and hexose, contributing to sugar allocation from seed coat to embryo, which consequently determines oil and protein contents and seed size in soybean. We conclude that past selection of optimal GmSWEET10a alleles drove the initial domestication of multiple soybean-seed traits and that targeted selection of the elite allele GmSWEET10b may further improve the yield and seed quality of modern soybean cultivars.

Elevation of soybean seed oil content through selection for seed coat shininess.

Many leguminous species have adapted their seed coat with a layer of powdery bloom that contains hazardous allergens and makes the seeds less visible, offering duel protection against potential predators 1 . Nevertheless, a shiny seed surface without bloom is desirable for human consumption and health, and is targeted for selection under domestication. Here we show that seed coat bloom in wild soybeans is mainly controlled by Bloom1 (B1), which encodes a transmembrane transporter-like protein for biosynthesis of the bloom in pod endocarp. The transition from the 'bloom' to 'no-bloom' phenotypes is associated with artificial selection of a nucleotide mutation that naturally occurred in the coding region of B1 during soybean domestication. Interestingly, this mutation not only 'shined' the seed surface, but also elevated seed oil content in domesticated soybeans. Such an elevation of oil content in seeds appears to be achieved through b1-modulated upregulation of oil biosynthesis in pods. This study shows pleiotropy as a mechanism underlying the domestication syndrome 2 , and may pave new strategies for development of soybean varieties with increased seed oil content and reduced seed dust.

Functional conservation and divergence of GmCHLI genes in polyploid soybean.

Polyploidy is prevalent in nature. As the fate of duplicated genes becomes more complicated when the encoded proteins function as oligomers, functional investigations into duplicated oligomer-encoding genes in polyploid genomes will facilitate our understanding of how traits are expressed. In this study, we identified GmCHLI1, a gene encoding the I subunit of magnesium (Mg)-chelatase, which functions in hexamers as responsible for the semi-dominant etiolation phenotype in soybean. Four GmCHLI copies derived from two polyploidy events were identified in the soybean genome. Further investigation with regard to expression patterns indicated that these four copies have diverged into two pairs; mutation in the other copy of the pair that includes GmCHLI1 also resulted in a chlorophyll-deficient phenotype. Protein interaction assays showed that these four GmCHLIs can interact with each other, but stronger interactions were found with mutated subunits. The results indicate that, in polyploidy, deficiency in each copy of duplicated oligomer-encoding genes could result in a mutant phenotype due to hetero-oligomer formation, which is different from the model of allelic dosage or functional redundancy. In addition, we interestingly found an increase in isoflavonoids in the heterozygous etiolated plants, which might be useful for improving soybean seed quality.

Global dissection of alternative splicing in paleopolyploid soybean.

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.

Concerted evolution of D1 and D2 to regulate chlorophyll degradation in soybean.

Polyploidy is a common phenomenon, particularly in plants. The soybean (Glycine max [L.] Merr.) genome has undergone two whole genome duplication (WGD) events. The conservation and divergence of duplicated gene pairs are major contributors to genome evolution. D1 and D2 are two unlinked, paralogous nuclear genes, whose double-recessive mutant (d1d1d2d2) results in chlorophyll retention, called 'stay-green'. Through molecular cloning and functional analyses, we demonstrated that D1 and D2 are homologs of the STAY-GREEN (SGR) genes from other plant species and were duplicated as a result of the most recent WGD in soybean. Transcriptional analysis showed that both D1 and D2 were more highly expressed in older tissues, and chlorophyll degradation and programmed cell death-related genes were suppressed in a d1d2 double mutant, this situation indicated that these genes are probably involved in the early stages of tissue senescence. Investigation of genes that flank D1 and D2 revealed that evolution within collinear duplicated blocks may affect the conservation of individual gene pairs within the blocks. Moreover, we found that a long terminal repeat retrotransposon, GmD2IN, resulted in the d2 mutation. Further analysis of this retrotransposon family showed that insertion in or near the coding regions can affect gene expression or splicing patterns, and may be an important force to promote the divergence of duplicated gene pairs.

Comprehensive analyses of microRNA gene evolution in paleopolyploid soybean genome.

miRNA genes are thought to undergo quick birth and death processes in genomes and the emergence of a MIRNA-like hairpin provides the base for functional miRNA gene formation. However, the factors affecting the formation of an active miRNA gene from a MIRNA-like hairpin within a genome remain unclear. We performed a genome-wide investigation of MIRNA-like hairpin accumulation, expression, structural changes and relationships with annotated genomic features in the paleopolyploid soybean genome. Our results showed that adjacent gene and transposable element content, rates of genetic recombination at location of emergence, along with its own gene structure divergence greatly affected miRNA gene evolution. Further investigation suggested that miRNA genes from different duplication sources followed distinct evolutionary trajectories and that the accumulation of MIRNA-like hairpins might be an important factor in causing long terminal repeat retrotransposons to lose activity during genome evolution.

Pericentromeric effects shape the patterns of divergence, retention, and expression of duplicated genes in the paleopolyploid soybean.

The evolutionary forces that govern the divergence and retention of duplicated genes in polyploids are poorly understood. In this study, we first investigated the rates of nonsynonymous substitution (Ka) and the rates of synonymous substitution (Ks) for a nearly complete set of genes in the paleopolyploid soybean (Glycine max) by comparing the orthologs between soybean and its progenitor species Glycine soja and then compared the patterns of gene divergence and expression between pericentromeric regions and chromosomal arms in different gene categories. Our results reveal strong associations between duplication status and Ka and gene expression levels and overall low Ks and low levels of gene expression in pericentromeric regions. It is theorized that deleterious mutations can easily accumulate in recombination-suppressed regions, because of Hill-Robertson effects. Intriguingly, the genes in pericentromeric regions-the cold spots for meiotic recombination in soybean-showed significantly lower Ka and higher levels of expression than their homoeologs in chromosomal arms. This asymmetric evolution of two members of individual whole genome duplication (WGD)-derived gene pairs, echoing the biased accumulation of singletons in pericentromeric regions, suggests that distinct genomic features between the two distinct chromatin types are important determinants shaping the patterns of divergence and retention of WGD-derived genes.

Genetic analysis of starch paste viscosity parameters in glutinous rice (Oryza sativa L.).

Starch paste viscosity plays an important role in estimating the cooking, eating, and processing quality of rice. The inheritance of starch paste viscosity in glutinous rice remains undefined. In the present study, 118 glutinous rice accessions were collected, and the genotypes of 17 starch synthesis-related genes (SSRG) were analyzed by using 43 gene-specific molecular markers. Association analysis indicated that 10 of 17 SSRGs were involved in controlling the rapid visco analyzer (RVA) profile parameters. Among these, the PUL gene was identified to play an important role in control of peak viscosity (PKV), hot paste viscosity (HPV), cool paste viscosity (CPV), breakdown viscosity (BDV), peak time (PeT), and paste temperature (PaT) in glutinous rice. Other SSRGs involved only a few RVA profile parameters. Furthermore, interactions between SSRGs were found being responsible for PeT, PaT, and BDV. Some of the RVA parameters, including PKV, HPV, CPV, CSV, and PaT, were mainly governed by single SSRG, whereas other parameters, such as BDV, SBV, and PeT, were controlled by a few SSRGs, functioning cooperatively. Further, three near-isogenic lines (NIL) of a japonica glutinous cv. Suyunuo as genetic background, with PUL, SSIII-1, and SSIII-2 alleles replaced with those of indica cv. Guichao 2, were employed to verify the genetic effects of the various genes, and the results were consistent with those obtained from the association analysis. These findings indicated that starch paste viscosity in glutinous rice had a complex genetic system, and the PUL gene played an important role in determining the RVA profile parameters in glutinous rice. These results provide important information for potentially improving the quality of glutinous rice.

Genome sequence of the palaeopolyploid soybean.

Soybean (Glycine max) is one of the most important crop plants for seed protein and oil content, and for its capacity to fix atmospheric nitrogen through symbioses with soil-borne microorganisms. We sequenced the 1.1-gigabase genome by a whole-genome shotgun approach and integrated it with physical and high-density genetic maps to create a chromosome-scale draft sequence assembly. We predict 46,430 protein-coding genes, 70% more than Arabidopsis and similar to the poplar genome which, like soybean, is an ancient polyploid (palaeopolyploid). About 78% of the predicted genes occur in chromosome ends, which comprise less than one-half of the genome but account for nearly all of the genetic recombination. Genome duplications occurred at approximately 59 and 13 million years ago, resulting in a highly duplicated genome with nearly 75% of the genes present in multiple copies. The two duplication events were followed by gene diversification and loss, and numerous chromosome rearrangements. An accurate soybean genome sequence will facilitate the identification of the genetic basis of many soybean traits, and accelerate the creation of improved soybean varieties.