Molecular evolution of glycinin and β-conglycinin gene families in soybean (Glycine max L. Merr.)

Soybean genetics, genomics, and breeding for improving nutritional value and reducing antinutritional traits in food and feed

Soybean [Glycine max (L.) Merr.] is a globally important crop due to its valuable seed composition, versatile feed, food, and industrial end-uses, and consistent genetic gain. Successful genetic gain in soybean has led to widespread adaptation and increased value for producers, processors, and consumers. Specific focus on the nutritional quality of soybean seed composition for food and feed has further elucidated genetic knowledge and bolstered breeding progress. Seed components are historical and current targets for soybean breeders seeking to improve nutritional quality of soybean. This article reviews genetic and genomic foundations for improvement of nutritionally important traits, such as protein and amino acids, oil and fatty acids, carbohydrates, and specific food-grade considerations; discusses the application of advanced breeding technology such as CRISPR/Cas9 in creating seed composition variations; and provides future directions and breeding recommendations regarding soybean seed composition traits.

A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max)

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to β-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

Assessment of In Vitro Digestive Behavior of Lactic-Acid-Bacteria Fermented Soy Proteins: A Study Comparing Colloidal Solutions and Curds

This study investigated the effect of lactic-acid-bacteria fermentation on the microstructure and gastrointestinal digestibility of soy proteins using a digestomics approach. Fermented soy protein isolates (FSPIs) under varied fermentation-terminal pH demonstrated a colloidal solution (FSPI-7.0/6.0) or yogurt-like curd (FSPI-5.0/4.0) state. Cryo-electron microscopy figures demonstrated the loosely stacked layer of FSPI-7.0/6.0 samples, whereas a denser gel network was observed for FSPI-5.0/4.0 samples. Molecular interactions shifted from dominant ionic bonds to hydrophobic forces and disulfide bonds. The gastric/intestinal digestion demonstrated that the curd samples afforded a significantly low particle size and high-soluble protein and peptide contents in the medium and late digestive phases. A peptidomics study showed that the FSPI-6.0 digestate at early intestinal digestion had a high peptidome abundance, whereas FSPI curd digestates (FSPI-5.0/4.0) elicited a postponed but more extensive promotion during medium and late digestion. Glycinin G2/G4 and β-conglycinin α/α' subunits were the major subunits promoted by FSPI-curds. The spatial structures of glycinin G2 and β-conglycinin α subunits demonstrated variations located in seven regions. Glycinin G2 region 6 (A349-K356) and β-conglycinin α subunit region 7 (E556-E575), which were located at the interior of the 3D structure, were the key regions contributing to discrepancies at the late stage.

Confocal Fluorescence Microscopy Investigation for the Existence of Subdomains within Protein Storage Vacuoles in Soybean Cotyledons

In legumes, the seed storage proteins accumulate within specialized organelles called protein storage vacuoles (PSVs). In several plant species, PSVs are differentiated into subdomains that accumulate different kinds of proteins. Even though the existence of subdomains is common in cereals and legumes, it has not been reported in soybean PSVs. The two most abundant seed proteins of soybean, 7S and 11S globulins, have different temporal accumulation patterns and exhibit considerable solubility differences that could result in differential accretion of these proteins within the PSVs. Here, we employed confocal fluorescent microscopy to examine the presence or absence of subdomains within the soybean PSVs. Eosin-stained sections of FAA-fixed paraffin embedded soybean seeds, when viewed by confocal fluorescence microscopy, revealed the presence of intricate subdomains within the PSVs. However, fluorescence immunolabeling studies demonstrated that the 7S and 11S globulins were evenly distributed within the PSVs and failed to corroborate the existence of subdomains within the PSVs. Similarly, confocal scanning microscopy examination of free-hand, vibratome and cryostat sections also failed to demonstrate the existence of subdomains within PSVs. The subdomains, which were prominently seen in PSVs of FAA-fixed soybean seeds, were not observed when the seeds were fixed either in glutaraldehyde/paraformaldehyde or glutaraldehyde. Our studies demonstrate that the apparent subdomains observed in FAA-fixed seeds may be a fixation artifact.

Linkage and association study discovered loci and candidate genes for glycinin and β-conglycinin in soybean (Glycine max L. Merr.)

KEY MESSAGE

Linkage mapping and GWAS identified 67 QTLs related to soybean glycinin, β-conglycinin and relevant traits. Polymorphisms of the candidate gene Gy1 promoter were associated with the glycinin content in soybean. The major components of storage proteins in soybean seeds are glycinin and β-conglycinin, which play important roles in determining protein nutrition and soy food processing properties. Increasing the protein content while improving the ratio of glycinin to β-conglycinin is substantially important for soybean protein improvement. To investigate the genetic mechanism of storage protein subunits, 184 recombinant inbred lines (RILs) derived from a cross of Kefeng No. 1 and Nannong 1138-2 and 211 diverse soybean cultivars were used to detect loci related to glycinin (11S), β-conglycinin (7S), the sum of glycinin and β-conglycinin (SGC), and the ratio of glycinin to β-conglycinin (RGC). Sixty-seven QTLs and 11 hot genomic regions were identified as affecting the four traits. One genetic region (q10-1) on chromosome 10 was associated with multiple traits by both linkage and association analysis. Eight genes in 11 hot genomic regions might be related to soybean protein subunit. The candidate gene analysis showed that polymorphisms in Gy1 promoters were significantly correlated with the 11S content. The QTLs and candidate genes identified in the present study allow for further understanding the genetic basis of 11S and 7S regulation and provide useful information for marker-assisted selection (MAS) in soybean quality improvement.

Genome-wide association mapping for protein, oil and water-soluble protein contents in soybean

As a globally important legume crop, soybean provides excellent sources of protein and oil for human and livestock nutrition. Improving seed protein and oil contents has always been an important objective in soybean breeding. Water-soluble protein plays a significant role in the processing and efficacy of soybean protein. Here, a genome-wide association study (GWAS) of seed compositions (protein, oil, and water-soluble protein contents) was conducted using 211 diverse soybean accessions genotyped with a 355 K SoySNP array. Three, four, and five QTLs were identified related to the protein, oil, and water-soluble protein contents, respectively. Furthermore, five QTLs (qPC-15-1, qOC-8-1, qOC-12-1, qOC-20-1 and qWSPC-8-1) were detected in multiple environments. Analysis of the favorable alleles for oil and water-soluble protein contents showed that qOC-8-1 (qWSPC-8-1) exerted inverse effects on oil and water-soluble protein synthesis. Relative expression analysis suggested that Glyma.15G049200 in qPC-15-1 affects protein synthesis and Glyma.08G107800 in qOC-8-1 and qWSPC-8-1 might be involved in oil and water-soluble protein synthesis, producing opposite effects. The candidate genes and significant SNPs detected in the present study will allow a deeper understanding of the genetic basis for the regulation of protein, oil and water-soluble protein contents and provide important information that could be utilized in marker-assisted selection for soybean quality improvement.

Genome-wide association identifies several QTLs controlling cysteine and methionine content in soybean seed including some promising candidate genes

Soybean is an important source of protein, oil and carbohydrates, as well as other beneficial nutrients. A major function of proteins in nutrition is to supply adequate amounts of amino acids. Although they are essential for human nutrition, the sulfur-containing amino acids cysteine (Cys) and methionine (Met) are often limited and the genetic control of their content in soybean seeds is poorly characterized. This study aimed to characterize the phenotypic variation and identify quantitative trait loci (QTL) associated with Cys and Met content in a core set of 137 soybean lines, representative of the genetic diversity among Canadian short-season soybean, spanning maturity groups 000-II (MG000-II). Significant phenotypic differences were found among these lines for Cys, Met and Cys + Met content. Using both a mixed linear model and six multi-locus methods with a catalogue of 2.18 M SNPs, we report a total of nine QTLs and seventeen QTNs of which seven comprise promising candidate genes. This work allowed us to reproducibly detect multiple novel loci associated with sulfur-containing amino acid content. The markers and genes identified in this study may be useful for soybean genetic improvement aiming to increase Cys and Met content.

β-Conglycinin induces the formation of neutrophil extracellular traps dependent on NADPH oxidase-derived ROS, PAD4, ERK1/2 and p38 signaling pathways in mice

β-Conglycinin is one of the key thermostable anti-nutritional factors in soybean, which has strong immunogenicity that usually leads to weaning in some young animals such as piglets and calves and allergic reaction in rats. Neutrophils are involved in the pathogenesis of an allergy. However, the contribution of functional neutrophils to allergy needs to be clarified. The formation of neutrophil extracellular traps is a novel effector mechanism of neutrophils and has been extensively investigated in recent years. To the best of our knowledge, there is no information available on β-conglycinin-induced NETs. In this study, β-conglycinin-induced NET formation in mice was examined via immunofluorescence analysis and fluorescence microplate reader. The mechanism of β-conglycinin-induced NETs was investigated using inhibitors and fluorescent microplate methods. The results showed that β-conglycinin induced the classical characteristics of NETs, which mainly consist of DNA as the backbone and decorated with histones, myeloperoxidase (MPO) and neutrophil elastase (NE). Moreover, β-conglycinin significantly induced the formation of NETs in a dose-dependent way. NET degrading enzyme DNase I markedly reduced β-conglycinin-induced NETs, which suggests that β-conglycinin indeed triggered the release of NETs. Further investigation showed that the quantitation of NETs was markedly decreased by the inhibitors of reactive oxygen species (ROS)-derived-NADPH oxidase, ERK1/2, p38, Rac and PAD4 signaling pathways, indicating the crucial role of these signaling pathways in β-conglycinin-induced NETs. Furthermore, our findings revealed that β-conglycinin induced the formation of NETs, which is dependent on NADPH oxidase-derived ROS, ERK1/2, p38, Rac and PAD4 signaling pathways. This study is the first to demonstrate the underlying mechanisms of β-conglycinin-induced NET formation, and it could be helpful to understand diarrhea caused by β-conglycinin overexposure in young animals and provides the corresponding theoretical basis for clinical applications.

Protein analysis reveals differential accumulation of late embryogenesis abundant and storage proteins in seeds of wild and cultivated amaranth species

BACKGROUND

Amaranth is a plant naturally resistant to various types of stresses that produces seeds of excellent nutritional quality, so amaranth is a promising system for food production. Amaranth wild relatives have survived climate changes and grow under harsh conditions, however no studies about morphological and molecular characteristics of their seeds are known. Therefore, we carried out a detailed morphological and molecular characterization of wild species A. powellii and A. hybridus, and compared them with the cultivated amaranth species A. hypochondriacus (waxy and non-waxy seeds) and A. cruentus.

RESULTS

Seed proteins were fractionated according to their polarity properties and were analysed in one-dimensional gel electrophoresis (1-DE) followed by nano-liquid chromatography coupled to tandem mass spectrometry (nLC-MS/MS). A total of 34 differentially accumulated protein bands were detected and 105 proteins were successfully identified. Late embryogenesis abundant proteins were detected as species-specific. Oleosins and oil bodies associated proteins were observed preferentially in A. cruentus. Different isoforms of the granule-bound starch synthase I, and several paralogs of 7S and 11S globulins were also identified. The in silico structural analysis from different isoforms of 11S globulins was carried out, including new types of 11S globulin not reported so far.

CONCLUSIONS

The results provide novel information about 11S globulins and proteins related in seed protection, which could play important roles in the nutritional value and adaptive tolerance to stress in amaranth species.

Analysis of genes encoding seed storage proteins (SSPs) in chickpea (Cicer arietinum L.) reveals co-expressing transcription factors and a seed-specific promoter

Improvement of the quality and quantity of chickpea seed protein can be greatly facilitated by an understanding of the genic organization and the genetic architecture of the genes encoding seed storage proteins (SSPs). The aim of this study was to provide a comprehensive analysis of the chickpea SSP genes, putative co-expressing transcription factors (TFs), and to identify a seed-specific SSP gene promoter. A genome-wide identification of SSP genes in chickpea led to the identification of 21 non-redundant SSP encoding genes located on 6 chromosomes. Phylogenetic analysis grouped SSP genes into 3 subgroups where members within the same clade demonstrated similar motif composition and intron-exon organization. Tandem duplications were identified to be the major contributors to the expansion of the SSP gene family in chickpea. Co-expression analysis revealed 14 TFs having expression profiles similar to the SSP genes that included members of important TF families that are known to regulate seed development. Expression analysis of SSP genes and TFs revealed significantly higher expression in late stages of seed development as well as in high seed protein content (HPC) genotypes. In silico analysis of the promoter regions of the SSP encoding genes revealed several seed-specific cis-regulatory elements such as RY repeats, ACGT motifs, CAANTG, and GCN4. A candidate promoter was analyzed for seed specificity by generating stable transgenics in Arabidopsis. Overall, this study provides a useful resource to explore the regulatory networks involved in SSP synthesis and/or accumulation for utilization in developing nutritionally improved chickpea genotypes.

Genetic mapping and validation of the loci controlling 7S α' and 11S A-type storage protein subunits in soybean [Glycine max (L.) Merr.]

KEY MESSAGE

Four soybean storage protein subunit QTLs were mapped using bulked segregant analysis and an F2 population, which were validated with an F5 RIL population. The storage protein globulins β-conglycinin (7S subunit) and glycinin (11S subunits) can affect the quantity and quality of proteins found in soybean seeds and account for more than 70% of the total soybean protein. Manipulating the storage protein subunits to enhance soymeal nutrition and for desirable tofu manufacturing characteristics are two end-use quality goals in soybean breeding programs. To aid in developing soybean cultivars with desired seed composition, an F2 mapping population (n = 448) and an F5 RIL population (n = 180) were developed by crossing high protein cultivar 'Harovinton' with the breeding line SQ97-0263_3-1a, which lacks the 7S α', 11S A1, 11S A2, 11S A3 and 11S A4 subunits. The storage protein composition of each individual in the F2 and F5 populations were profiled using SDS-PAGE. Based on the presence/absence of the subunits, genomic DNA bulks were formed among the F2 plants to identify genomic regions controlling the 7S α' and 11S protein subunits. By utilizing polymorphic SNPs between the bulks characterized with Illumina SoySNP50K iSelect BeadChips at targeted genomic regions, KASP assays were designed and used to map QTLs causing the loss of the subunits. Soybean storage protein QTLs were identified on Chromosome 3 (11S A1), Chromosome 10 (7S α' and 11S A4), and Chromosome 13 (11S A3), which were also validated in the F5 RIL population. The results of this research could allow for the deployment of marker-assisted selection for desired storage protein subunits by screening breeding populations using the SNPs linked with the subunits of interest.

Development and Characterization of a Soybean Experimental Line Lacking the α' Subunit of β-Conglycinin and G1, G2, and G4 Glycinin

A soybean experimental line (BSH-3) devoid of a subset of seed storage proteins was developed by crossing a mutant donor line "HS99B" with a Chinese cultivar "Dongnong47" (DN47). One-dimensional and high-resolution 2-D gel electrophoresis revealed the absence of G1 (A1aB2), G2 (A2B1a), and G4 (A5A4B3) glycinin and the α' subunit of β-conglycinin in BSH-3 seeds. Despite the lack of these abundant seed proteins, BSH-3 seeds still accumulated 38% protein. BSH-3 seeds also accumulated high levels of free amino acids as compared with DN47 seeds, particularly arginine, and the amount of several essential amino acids were significantly elevated in BSH-3 seeds. Elevated accumulation of α and β-subunit of β-conglycinin, G5 glycinin, Kunitz trypsin inhibitor, and Bowman-Birk protease inhibitor indicates seed proteome rebalancing in BSH-3 seeds. Immunoblot analysis using sera from soybean allergic patients demonstrated the complete lack of a major allergen (α' subunit of β-conglycinin) in BSH-3 seeds. However, elevated levels of other allergens were found in BSH-3 seeds due to proteome rebalancing. Transmission electron microscopy observation of mature seeds of BSH-3 revealed striking differences in the appearance of the protein storage vacuoles when compared with DN47.

Plant Genome Duplication Database

Genome duplication, widespread in flowering plants, is a driving force in evolution. Genome alignments between/within genomes facilitate identification of homologous regions and individual genes to investigate evolutionary consequences of genome duplication. PGDD (the Plant Genome Duplication Database), a public web service database, provides intra- or interplant genome alignment information. At present, PGDD contains information for 47 plants whose genome sequences have been released. Here, we describe methods for identification and estimation of dates of genome duplication and speciation by functions of PGDD.The database is freely available at http://chibba.agtec.uga.edu/duplication/.

Quantitative Trait Loci (QTL) Mapping for Glycinin and β-Conglycinin Contents in Soybean (Glycine max L. Merr.)

Compared to β-conglycinin, glycinin contains 3-4 times the methionine and cysteine (sulfur-containing amino acids), accounting for approximately 40 and 30%, respectively, of the total storage protein in soybean. Increasing the soybean storage protein content while improving the ratio of glycinin to β-conglycinin is of great significance for soybean breeding and soy food products. The objective of this study is to analyze the genetic mechanism regulating the glycinin and β-conglycinin contents of soybean by using a recombinant inbred line (RIL) population derived from a cross between Kefeng No. 1 and Nannong 1138-2. Two hundred and twenty-one markers were used to map quantitative trait loci (QTLs) for glycinin (11S) and β-conglycinin (7S) contents, the ratio of glycinin to β-conglycinin (RGC), and the sum of glycinin and β-conglycinin (SGC). A total of 35 QTLs, 3 pairs of epistatic QTLs, and 5 major regions encompassing multiple QTLs were detected. Genes encoding the subunits of β-conglycinin were localized to marker intervals sat_418-satt650 and sat_196-sat_303, which are linked to RGC and SGC; marker sat_318, associated with 11S, 7S, and SGC, was located near Glyma10g04280 (Gy4), which encodes a subunit of glycinin. These results, which take epistatic interactions into account, will improve our understanding of the genetic basis of 11S and 7S contents and will lay a foundation for marker-assisted selection (MAS) breeding of soybean and improving the quality of soybean products.

Positive selection drives neofunctionalization of the UbiA prenyltransferase gene family

Gene duplication provides the key materials for new genes and novel functions. However, the mechanism underlying functional innovation remains unknown. In this study, we revealed the evolutionary pattern of the prenyltransferases of the UbiA gene family in 15 higher plants. Prenyltransferases of the UbiA gene family are involved in many important biological processes of both primary and secondary metabolism. Based on the phylogenetic relationships of the UbiA genes, seven subfamilies are classified. Confirming this classification, genes within each subfamily are characterized by similar exon numbers, exon lengths and patterns of motif combinations. Similar numbers of UbiA genes are found in different species within each subfamily except for Subfamily I, in which a Phaseoleae-specific expansion is detected in clade I-A. Homologous genes in clade I-A evolve rapidly, exchange sequences frequently and experience positive selection. Genes in clade I-A function as flavonoid prenyltransferase synthesis secondary compounds, while other genes from Subfamily I encode homogentisate phytyltransferase, which plays a role in primary metabolism. Thus, our results suggest that the secondary metabolism genes acquire new functions from those of primary metabolism through gene duplication and neofunctionalization driven by positive selection.

Structural and functional analysis of various globulin proteins from soy seed

Storage proteins of soybean mostly consist of globulins, which are classified according to their sedimentation coefficient. Among 4 major types: 2S, 7S, 11S, and 15S of globulins, 7S and 11S constitute major fraction. The 11S fraction consists only of glycinin and 7S fraction majorly consists of β-conglycinin, small amounts of γ-conglycinin and basic 7S globulin (Bg7S). Glycinin exist as a hexamer while β-conglycinin as a trimer and Bg7S as a tetramer. Glycinin subunits are coded by 5 genes of a family, whereas about 15 genes are present for β-conglycinin subunits. Bg7S gene is present in four copies in soybean genome. Synthesis of all proteins takes place as a single polypeptide chain, which is cleaved after folding to yield different chains or subunits. Glycinin and β-Conglycinin are made for storage purpose. However, Bg7S has potential xylanase inhibition activity and protein kinase activity. Primary structure of Bg7S reveals 12 conserved cysteine residues involved in forming 6 disulfide bonds, which provides appreciable stability to protein. Secondary structure is predominately rich in β-sheets with few alpha helices. Bg7S shares structural similarity with various aspartic-proteases. In this review, our aim is to discuss sequence, structure, and function of various globulins present in Glycine max.

Profiles of embryonic nuclear protein binding to the proximal promoter region of the soybean β-conglycinin α subunit gene

β-Conglycinin, a major component of seed storage protein in soybean, comprises three subunits: α, α' and β. The expression of genes for these subunits is strictly controlled during embryogenesis. The proximal promoter region up to 245 bp upstream of the transcription start site of the α subunit gene sufficiently confers spatial and temporal control of transcription in embryos. Here, the binding profile of nuclear proteins in the proximal promoter region of the α subunit gene was analysed. DNase I footprinting analysis indicated binding of proteins to the RY element and DNA regions including box I, a region conserved in cognate gene promoters. An electrophoretic mobility shift assay (EMSA) using different portions of box I as a probe revealed that multiple portions of box I bind to nuclear proteins. In addition, an EMSA using nuclear proteins extracted from embryos at different developmental stages indicated that the levels of major DNA-protein complexes on box I increased during embryo maturation. These results are consistent with the notion that box I is important for the transcriptional control of seed storage protein genes. Furthermore, the present data suggest that nuclear proteins bind to novel motifs in box I including 5'-TCAATT-3' rather than to predicted cis-regulatory elements.

Knockdown of the 7S globulin subunits shifts distribution of nitrogen sources to the residual protein fraction in transgenic soybean seeds

KEY MESSAGE

A platform of gene silencing by amiRNA had been established in fertile transgenic soybean. We demonstrated that knockdown of storage protein shifted the distribution of nitrogen sources in soybean seeds. Artificial microRNAs (amiRNAs) were designed using the precursor sequence of the endogenous soybean (Glycine max L. Merrill) miRNA gma-miR159a and expressed in transgenic soybean plants to suppress the biosynthesis of 7S globulin, which is one of the major storage proteins. Seed-specific expression of these amiRNAs (amiR-7S) resulted in a strong suppression of 7S globulin subunit genes and decreased accumulation of the 7S globulin subunits in seeds. Thus, the results demonstrate that a platform for gene silencing by amiRNA was first developed in fertile transgenic soybean plants. There was no difference in nitrogen, carbon, and lipid contents between amiR-7S and control seeds. Four protein fractions were collected from defatted mature seeds on the basis of solubility at different pH to examine the distribution of nitrogen sources and compensatory effects. In the whey and lipophilic fractions, nitrogen content was similar in amiR-7S and control seeds. Nitrogen content was significantly decreased in the major soluble protein fraction and increased in the residual fraction (okara) of the amiR-7S seeds. Amino acid analysis revealed that increased nitrogen compounds in okara were proteins or peptides rather than free amino acids. Our study indicates that the decrease in 7S globulin subunits shifts the distribution of nitrogen sources to okara in transgenic soybean seeds.

Comparative and evolutionary analysis of major peanut allergen gene families

Peanut (Arachis hypogaea L.) causes one of the most serious food allergies. Peanut seed proteins, Arah1, Arah2, and Arah3, are considered to be among the most important peanut allergens. To gain insights into genome organization and evolution of allergen-encoding genes, approximately 617 kb from the genome of cultivated peanut and 215 kb from a wild relative were sequenced including three Arah1, one Arah2, eight Arah3, and two Arah6 gene family members. To assign polarity to differences between homoeologous regions in peanut, we used as outgroups the single orthologous regions in Medicago, Lotus, common bean, chickpea, and pigeonpea, which diverged from peanut about 50 Ma and have not undergone subsequent polyploidy. These regions were also compared with orthologs in many additional dicot plant species to help clarify the timing of evolutionary events. The lack of conservation of allergenic epitopes between species, and the fact that many different proteins can be allergenic, makes the identification of allergens across species by comparative studies difficult. The peanut allergen genes are interspersed with low-copy genes and transposable elements. Phylogenetic analyses revealed lineage-specific expansion and loss of low-copy genes between species and homoeologs. Arah1 syntenic regions are conserved in soybean, pigeonpea, tomato, grape, Lotus, and Arabidopsis, whereas Arah3 syntenic regions show genome rearrangements. We infer that tandem and segmental duplications led to the establishment of the Arah3 gene family. Our analysis indicates differences in conserved motifs in allergen proteins and in the promoter regions of the allergen-encoding genes. Phylogenetic analysis and genomic organization studies provide new insights into the evolution of the major peanut allergen-encoding genes.

Identification and characterization of transcript polymorphisms in soybean lines varying in oil composition and content

BACKGROUND

Variation in seed oil composition and content among soybean varieties is largely attributed to differences in transcript sequences and/or transcript accumulation of oil production related genes in seeds. Discovery and analysis of sequence and expression variations in these genes will accelerate soybean oil quality improvement.

RESULTS

In an effort to identify these variations, we sequenced the transcriptomes of soybean seeds from nine lines varying in oil composition and/or total oil content. Our results showed that 69,338 distinct transcripts from 32,885 annotated genes were expressed in seeds. A total of 8,037 transcript expression polymorphisms and 50,485 transcript sequence polymorphisms (48,792 SNPs and 1,693 small Indels) were identified among the lines. Effects of the transcript polymorphisms on their encoded protein sequences and functions were predicted. The studies also provided independent evidence that the lack of FAD2-1A gene activity and a non-synonymous SNP in the coding sequence of FAB2C caused elevated oleic acid and stearic acid levels in soybean lines M23 and FAM94-41, respectively.

CONCLUSIONS

As a proof-of-concept, we developed an integrated RNA-seq and bioinformatics approach to identify and functionally annotate transcript polymorphisms, and demonstrated its high effectiveness for discovery of genetic and transcript variations that result in altered oil quality traits. The collection of transcript polymorphisms coupled with their predicted functional effects will be a valuable asset for further discovery of genes, gene variants, and functional markers to improve soybean oil quality.

PGDD: a database of gene and genome duplication in plants

Genome duplication (GD) has permanently shaped the architecture and function of many higher eukaryotic genomes. The angiosperms (flowering plants) are outstanding models in which to elucidate consequences of GD for higher eukaryotes, owing to their propensity for chromosomal duplication or even triplication in a few cases. Duplicated genome structures often require both intra- and inter-genome alignments to unravel their evolutionary history, also providing the means to deduce both obvious and otherwise-cryptic orthology, paralogy and other relationships among genes. The burgeoning sets of angiosperm genome sequences provide the foundation for a host of investigations into the functional and evolutionary consequences of gene and GD. To provide genome alignments from a single resource based on uniform standards that have been validated by empirical studies, we built the Plant Genome Duplication Database (PGDD; freely available at http://chibba.agtec.uga.edu/duplication/), a web service providing synteny information in terms of colinearity between chromosomes. At present, PGDD contains data for 26 plants including bryophytes and chlorophyta, as well as angiosperms with draft genome sequences. In addition to the inclusion of new genomes as they become available, we are preparing new functions to enhance PGDD.

The β-conglycinin deficiency in wild soybean is associated with the tail-to-tail inverted repeat of the α-subunit genes

β-conglycinin, a major seed protein in soybean, is composed of α, α', and β subunits sharing a high homology among them. Despite its many health benefits, β-conglycinin has a lower amino acid score and lower functional gelling properties compared to glycinin, another major soybean seed protein. In addition, the α, α', and β subunits also contain major allergens. A wild soybean (Glycine soja Sieb et Zucc.) line, 'QT2', lacks all of the β-conglycinin subunits, and the deficiency is controlled by a single dominant gene, Scg-1 (Suppressor of β-conglycinin). This gene was characterized using a soybean cultivar 'Fukuyutaka', 'QY7-25', (its near-isogenic line carrying the Scg-1 gene), and the F₂ population derived from them. The physical map of the Scg-1 region covered by lambda phage genomic clones revealed that the two α-subunit genes, a β-subunit gene, and a pseudo α-subunit gene were closely organized. The two α-subunit genes were arranged in a tail-to-tail orientation, and the genes were separated by 197 bp in Scg-1 compared to 3.3 kb in the normal allele (scg-1). In addition, small RNA was detected in immature seeds of the mutants by northern blot analysis using an RNA probe of the α subunit. These results strongly suggest that β-conglycinin deficiency in QT2 is controlled by post-transcriptional gene silencing through the inverted repeat of the α subunits.

Gene duplication and an accelerated evolutionary rate in 11S globulin genes are associated with higher protein synthesis in dicots as compared to monocots

BACKGROUND

Seed storage proteins are a major source of dietary protein, and the content of such proteins determines both the quantity and quality of crop yield. Significantly, examination of the protein content in the seeds of crop plants shows a distinct difference between monocots and dicots. Thus, it is expected that there are different evolutionary patterns in the genes underlying protein synthesis in the seeds of these two groups of plants.

RESULTS

Gene duplication, evolutionary rate and positive selection of a major gene family of seed storage proteins (the 11S globulin genes), were compared in dicots and monocots. The results, obtained from five species in each group, show more gene duplications, a higher evolutionary rate and positive selections of this gene family in dicots, which are rich in 11S globulins, but not in the monocots.

CONCLUSION

Our findings provide evidence to support the suggestion that gene duplication and an accelerated evolutionary rate may be associated with higher protein synthesis in dicots as compared to monocots.

MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity

MCScan is an algorithm able to scan multiple genomes or subgenomes in order to identify putative homologous chromosomal regions, and align these regions using genes as anchors. The MCScanX toolkit implements an adjusted MCScan algorithm for detection of synteny and collinearity that extends the original software by incorporating 14 utility programs for visualization of results and additional downstream analyses. Applications of MCScanX to several sequenced plant genomes and gene families are shown as examples. MCScanX can be used to effectively analyze chromosome structural changes, and reveal the history of gene family expansions that might contribute to the adaptation of lineages and taxa. An integrated view of various modes of gene duplication can supplement the traditional gene tree analysis in specific families. The source code and documentation of MCScanX are freely available at http://chibba.pgml.uga.edu/mcscan2/.