Nucleotide sequences encoding glycinin B4 polypeptide of the cultivated soybean Glycine max and its presumed wild-type ancestor Glycine soja in connection with the origin of cultivated soybean

Nucleotide sequences of cDNAs encoding soybean glycinin B4 polypeptide were compared in three soybean cultivars and two plant introductions of wild soybean Glycine soja. Only two nucleotide substitutions were found in three cultivars G. max, as compared with G. max and G. soja having nucleotide sequences which contain four nucleotide substitutions. These data serve as additional evidence, at molecular level, indicating the origin of G. max from G. soja. On the other hand, the time period required for four nucleotide substitutions' accumulation, as calculated from parameters of molecular evolution of 11S globulins, is much longer than the term having passed after soybean domestication.     

Constituents of soybean cultivars differing in insect resistance

Leaf samples of two insect-resistant plant introductions (PI 227687 and PI 229358) and two susceptible (Ransom and Coker Hampton 266A) cultivars of soybeans, Glycine max, were analyzed at different growth stages for their contents of total nitrogen, carbohydrates, organic acids, and sterols. The two susceptible cultivars accumulated more total nitrogen and at a faster rate than did the two resistant plant introductions. At pod-filling, the two resistant cultivars had equivalent soluble carbohydrates and 33% more than the susceptible cultivars. The quantity of organic acids was essentially the same for the two susceptible cultivars. The resistant cultivars had distinctly different quantities from each other as well as from the susceptible cultivars. The quantity of total sterol of these soybean cultivars varied during the growth of the plant. The resistant cultivars accumulated sterol faster and by pod-filling contained from 20 to 50% more sterol than did the susceptible cultivars.     

Cloning and comparative analysis of the gene encoding diacylglycerol acyltransferase from wild type and cultivated soybean

Diacylglycerol acyltransferase (DGAT), as an important enzyme in triacylglycerol synthesis, catalyzes the final acylation of the Kennedy pathway. In the present study, the GmDGAT gene was cloned from Glycine max by using AtDGAT as a query to search against the soybean EST database and the rapid amplification of cDNA ends (RACE) method. Allelic genes were also isolated from 13 soybean accessions and the divergence of the deduced amino acid sequences were compared. The comparison reveals that although GmDGAT is a highly conserved protein, several differences of insertion/deletion were identified in the N-terminal region of the GmDGATs from various soybean accessions. In the C-terminal regions, a single amino acid mutation specific to both G. max and G. soja was also found. The GmDGAT genomic sequences were further cloned and the number and size of exons in the DGAT genomic sequence were very similar among different plant species, whereas the introns were more diverged. These results may have significance in elucidating the genetic diversity of the GmDGAT among the soybean subgenus.     

CO and O2 complexes of soybean leghemoglobins: pH effects upon infrared and visible spectra. Comparisons with CO and O2 complexes of myoglobin and hemoglobin.

The effects of pH upon infrared spectra [CO stretching frequency (vco) region] and visible spectra of the CO complexes of soybean leghemoglobins a, c1, and c2, sperm whale myoglobin, and human hemoglobin A are reported. The vco for leghemoglobin--CO complexes was 1947.5 cm-1 at neutral pH. At acid pH myoglobin-- and hemoglobin--CO complexes developed vco bands at 1966--1968 cm-1, whereas leghemoglobin--CO complexes developed vco bands at approximately 1957 cm-1. All pKapp co values determined by pH-dependent variation of vco fell in the range 4.0--4.6. The pKapp co values determined from visible spectra were consistent with vco-determined values except for that of myoglobin--CO (visible pKapp co = 5.8). The pKapp co values in the 4.0--4.6 range appear to be pK values of the distal histidines, while the visible pKapp co of myoglobin--CO appears to be the pK of a group other than the distal and proximal histidines. The data are consistent with a model in which protonation of the distal histidine permits protein-free heme FeCO geometry in leghemoglobin--CO complexes but not in myoglobin-- or hemoglobin--CO complexes. Thus the heme pockets of leghemoglobins appear to be more flexible than the heme pockets of myoglobin and hemoglobin. The effects of pH upon visible spectra of the O2 complexes of soybean leghemoglobins a, c1, and c2, sperm whale myoglobin, and human hemoglobin A also are reported. pKapp o2 values of approximately 5.5 (leghemoglobins) and 4.4 (hemoglobin) are probably the pK values of the distal histidines. Comparisons of pKapp o2 values with pKapp co values indicate a more flexible heme pocket in leghemoglobins than in hemoglobin. The O2 complex of leghemoglobin c2 differed significantly from the O2 complexes of leghemoglobins a and c1 in visible spectra and titration behavior. These differences might be associated with the small structural differences in the region between the E and F helixes of leghemoglobins.     

Biochemical responses of glyphosate resistant and susceptible soybean plants exposed to glyphosate

Glyphosate is a wide spectrum, non-selective, post-emergence herbicide. It acts on the shikimic acid pathway inhibiting 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), thus obstructing the synthesis of tryptophan, phenylalanine, tyrosine and other secondary products, leading to plant death. Transgenic glyphosate-resistant (GR) soybean [Glycine max (L.)] expressing an glyphosate-insensitive EPSPS enzyme has provided new opportunities for weed control in soybean production. The effect of glyphosate application on chlorophyll level, lipid peroxidation, catalase (CAT), ascorbate peroxidase (APX), guaiacol peroxidase (GOPX) and superoxide dismutase (SOD) activities, soluble amino acid levels and protein profile, in leaves and roots, was examined in two conventional (non-GR) and two transgenic (GR) soybean. Glyphosate treatment had no significant impact on lipid peroxidation, whilst the chlorophyll content decreased in only one non-GR cultivar. However, there was a significant increase in the levels of soluble amino acid in roots and leaves, more so in non-GR than in GR soybean cultivars. Root CAT activity increased in non-GR cultivars and was not altered in GR cultivars. In leaves, CAT activity was inhibited in one non-GR and one GR cultivar. GOPX activity increased in one GR cultivar and in both non-GR cultivars. Root APX activity increased in one GR cultivar. The soluble protein profiles as assessed by 1-D gel electrophoresis of selected non-GR and GR soybean lines were unaffected by glyphosate treatment. Neither was formation of new isoenzymes of SOD and CAT observed when these lines were treated by glyphosate. The slight oxidative stress generated by glyphosate has no relevance to plant mortality. The potential antioxidant action of soluble amino acids may be responsible for the lack of lipid peroxidation observed. CAT activity in the roots and soluble amino acids in the leaves can be used as indicators of glyphosate resistance.     

Expression of the narrow leaflet gene for yield and agronomic traits in soybean

Genes that affect plant form and function may be used to enhance the yield of soybean [Glycine max (L.) Merr.]. Most soybean cultivars have broad (ovate) leaflets. A single gene, ln, controls inheritance for the narrow leaflet characteristic. Narrow leaflet cultivars (ln/ln) also tend to have a higher percentage of four-seeded pods than ovate (Ln/Ln) leaflet cultivars. Heterozygous (Ln/ln) plants have a leaflet shape intermediate between narrow and ovate. Determining the agronomic effects of the narrow leaflet allele (ln) in the heterozygous (Ln/ln) condition in soybean may have applications in practical plant breeding. We studied an ovate leaflet and a narrow leaflet cultivar, crosses between them in the F(1) and F(2), and backcrosses to both cultivars. The ratio of leaflet width to leaflet length accurately distinguished among narrow, ovate, and intermediate leaflet plants in the F(2) and backcross generations. In the F(2) generation, differences occurred among plants with different leaflet morphology. Narrow leaflet plants produced more seeds per pod and lower seed weight than ovate leaflet plants. Narrow and ovate leaflet plants produced comparable numbers of pods per plant and plant yield. Compared to ovate leaflet plants, intermediate leaflet plants produced similar numbers of seeds per pod and seed weight. Intermediate leaflet plants produced significantly more pods per plant and plant yield than plants with either ovate or narrow leaflets. The heterozygous condition at the locus for leaflet morphology resulted in heterosis for plant yield and may be of benefit in association with commercialization and development of hybrid soybean.     

NMR studies of the conformations of leghemoglobins from soybean and lupin

Phase-sensitive two-dimensional NMR methods have been used to obtain extensive proton resonance assignments for the carbon monoxide complexes of lupin leghemoglobins I and II and soybean leghemoglobin a. The assigned resonances provide information on the solution conformations of the proteins, particularly in the vicinity of the heme. The structure of the CO complex of lupin leghemoglobin II in solution is compared with the X-ray crystal structure of the cyanide complex by comparison of observed and calculated ring current shifts. The structures are generally very similar but significant differences are observed for the ligand contact residues, Phe30, His63 and Val67, and for the proximal His97 ligand. Certain residues are disordered and adopt two interconverting conformations in lupin leghemoglobin II in solution. The proximal heme pocket structure is closely conserved in the lupin leghemoglobins I and II but small differences in conformation in the distal heme pocket are apparent. Larger conformational differences are observed when comparisons are made with the CO complex of soybean leghemoglobin. Altered protein-heme packing is indicated on the proximal side of the heme and some conformational differences are evident in the distal heme pocket. The small conformational differences between the three leghemoglobins probably contribute to the known differences in their O2 and CO association and dissociation kinetics. The heme pocket conformations of the three leghemoglobins are more closely related to each other than to sperm whale myoglobin. The most notable differences between the leghemoglobins and myoglobin are: (a) reduced steric crowding of the ligand binding site in the leghemoglobins, (b) different orientations of the distal histidine, and (c) small but significant differences in proximal histidine coordination geometry. These changes probably contribute to the large differences in ligand binding kinetics between the leghemoglobins and myoglobin.     

Cloning of a higher-plant plastid omega-6 fatty acid desaturase cDNA and its expression in a cyanobacterium.

Oligomers based on amino acids conserved between known plant omega-3 and cyanobacterium omega-6 fatty acid desaturases were used to screen an Arabidopsis cDNA library for related sequences. An identified clone encoding a novel desaturase-like polypeptide was used to isolate its homologs from Glycine max and Brassica napus. The plant deduced amino acid sequences showed less than 27% similarity to known plant omega-6 and omega-3 desaturases but more than 48% similarity to cyanobacterial omega-6 desaturase, and they contained putative plastid transit sequences. Thus, we deduce that the plant cDNAs encode the plastid omega-6 desaturase. The identity was supported by expression of the B. napus cDNA in cyanobacterium. Synechococcus transformed with a chimeric gene that contains a prokaryotic promoter fused to the rapeseed cDNA encoding all but the first 73 amino acids partially converted its oleic acid fatty acid to linoleic acid, and the 16:1(9c) fatty acid was converted primarily to 16:2(9c, 12) in vivo. Thus, the plant omega-6 desaturase, which utilizes 16:1(7c) in plants, can utilize 16:1(9c) in the cyanobacterium. The plastid and cytosolic homologs of plant omega-6 desaturases are much more distantly related than those of omega-3 desaturases.     

Indole-3-acetic acid oxidase and syringaldazine oxidase activities of peroxidase isozymes in soybean root nodules

Many isoperoxidases with indole-3-acetic acid oxidase (IAA) and syringaldazine oxidase activities were detected by polyacrylamide gel electrophoresis in soybean root nodules [ Glycine max (L.) Merrill, cv. Asgrow], detached at the onset of flowering. The kinetics of the two activities were studied with some of the isoperoxidases partially purified by ion exchange chromatography. IAA oxidase activity of the cationic isoforms showed a sigmoidal kinetic behaviour and a higher substrate affinity than the anionic ones, whereas typical saturation kinetics were found with an anionic fraction that contained leghemoglobins. So, nodule IAA oxidase activity may mainly be displayed by the cationic isoforms. These cationic isoperoxidases had high affinity towards syringaldazine and they also may be associated with cell wall rigidification.     

Highly variable patterns of linkage disequilibrium in multiple soybean populations

Prospects for utilizing whole-genome association analysis in autogamous plant populations appear promising due to the reported high levels of linkage disequilibrium (LD). To determine the optimal strategies for implementing association analysis in soybean (Glycine max L. Merr.), we analyzed the structure of LD in three regions of the genome varying in length from 336 to 574 kb. This analysis was conducted in four distinct groups of soybean germplasm: 26 accessions of the wild ancestor of soybean (Glycine soja Seib. et Zucc.); 52 Asian G. max Landraces, the immediate results of domestication from G. soja; 17 Asian Landrace introductions that became the ancestors of North American (N. Am.) cultivars, and 25 Elite Cultivars from N. Am. In G. soja, LD did not extend past 100 kb; however, in the three cultivated G. max groups, LD extended from 90 to 574 kb, likely due to the impacts of domestication and increased self-fertilization. The three genomic regions were highly variable relative to the extent of LD within the three cultivated soybean populations. G. soja appears to be ideal for fine mapping of genes, but due to the highly variable levels of LD in the Landraces and the Elite Cultivars, whole-genome association analysis in soybean may be more difficult than first anticipated.     

Tolerance to Heterodera glycines in Soybean

Fifty-four susceptible soybean, Glycine max, cultivars or plant introductions were evaluated for tolerance to H. glycines, the soybean cyst nematode (SCN). Seed yields of genotypes were compared in nematicide-treated (1,2-dibromo-3-chloropropane, 58 kg a.i./ha) and nontreated plots at two SCN-infested locations over 3 years. Distinct and consistent levels of tolerance to SCN were observed among soybean genotypes. PI 97100, an introduction from Korea, exhibited the highest level of tolerance with an average tolerance index ([yield in nontreated plot ÷ yield in nematicide-treated plot] × 100) of 96 over 2 years. Coker 156 and Wright had moderate levels of tolerance (range in index values 68 to 95) compared to the intolerant cuhivars Bragg and Coker 237 (range in index values 33 to 68). Most of the soybean genotypes evaluated were intolerant to SCN. The rankings of five genotypes for tolerance to SCN and Hoplolaimus columbus were similar. Tolerance for seed yield was more consistently correlated with tolerance for plant height (r = 0.55 to 0.64) than for seed weight (r = 0.23 to 0.65) among genotypes.     

Development of SNP markers and haplotype analysis of the candidate gene for rhg1, which confers resistance to soybean cyst nematode in soybean

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) is one of the most destructive pests in the cultivation of soybean (Glycine max (L.) Merr.) worldwide. Markers based on the SCN resistance gene will enable efficient marker-assisted selection (MAS). We sequenced the candidate gene rhg1 in six resistant and two susceptible soybean genotypes and identified 37 SNPs (single nucleotide polymorphisms) among the sequences, of which 11 were in the coding region. Seven of these 11 SNPs led to changes in the amino acid sequence of the gene. The amino acid sequence we obtained differs from the previously published one by a stretch of 26–27 amino acids. Six codominant allele-specific SNP markers based on agarose gel detection were developed and tested in 70 genotypes, among which occurred only nine different haplotypes. Two neutrality tests (Tajima’s D and Fu and Li’s F) were significant for the six SNP loci in the 70 genotypes, which is consistent with intensive directional selection. A strong LD pattern was detected among five SNPs except 2868T > C. Two SNPs (689C > A and 757C > T) formed one haplotype (689C-757C) that was perfectly associated with SCN resistance. The new allele-specific PCR markers located in the alleged sequence of the rhg1 candidate gene, combined with the microsatellite marker BACR-Satt309, will significantly improve the efficiency of MAS during the development of SCN-resistant cultivars.     

Molecular cloning and functional expression of soybean allene oxide synthases

A plant allene oxide synthase (AOS) reacting with 13S-hydroperoxy-9Z,11E,15Z-octadecatrienoic acid (13-HPOT), a lipoxygenase product of alpha-linolenic acid, provides an allene oxide which functions as an intermediate for jasmonic acid (JA) synthesis, making AOS a key enzyme regulating the JA level in plants. Although AOSs in various plants have been investigated, there is only limited information about AOSs in soybean (Glycine max). In this study, we cloned and characterized two soybean AOSs, GmAOS1 and GmAOS2, sharing 95% homology in the predicted amino acid sequences. GmAOS1 and GmAOS2 were composed of 564 and 559 amino acids respectively, with predicted N-terminal chloroplast-targeting signal peptides. Both AOSs expressed in Escherichia coli were selective for 13S-hydroperoxides of alpha-linolenic and linoleic acids, suggesting the potential of GmAOS1 and GmAOS2 to contribute to JA synthesis. GmAOS1 and GmAOS2 were expressed in leaves, stems, and roots, suggesting broad distribution in a soybean plant.     

Genetic divergence between North American ancestral soybean lines and introductions with resistance to soybean cyst nematode revealed by chloroplast haplotype

Domesticated soybean [Glycine max (L.) Merr.] is a major crop with an established ancestral relationship to wild soybean (Glycine soja Sieb. & Zucc.) native to Asia. Soybean genetic diversity can be assessed at different levels by identification of polymorphic alleles at genetic loci, in either the plastid or nuclear genomes. The objective of this study was to evaluate genetic diversity based on chloroplast haplotypes for soybean genotypes present in the USDA germplasm resource collection. Shared chloroplast haplotypes represent broad groups of genetic relatedness. Previous work categorized three-quarters of the cultivated soybeans from Asia into a single haplotype group. Our results confirmed the close relationship of North American soybean ancestors and G. max plant introductions previously identified as representing potential sources of soybean genetic variation with the finding that these genotypes belonged to a single chloroplast haplotype group. Genetic diversity was identified in soybean genotypes determined to have a high density of single nucleotide polymorphisms and in a screen of accessions with resistance to soybean cyst nematode. Characterization of soybean plant introduction lines into chloroplast haplotype group may be an important initial step in evaluating the appropriate use of particular soybean genotypes.     

Comparative study on amino acid sequences of Kunitz-type soybean trypsin inhibitors, Tia, Tib, and Tic

The amino acid sequences of three variants of the Kunitz-type trypsin inhibitors, Tia, Tib, and Tic, obtained from some cultivars of soybean were determined by conventional methods. All three inhibitors consisted of 181 amino acid residues. The differences in the amino acid sequences are as follows: Tia E12 G55 Y62 H71 S74 M114 L120 P137 L176; Tib S F N R V I T V; Tic E. The amino acid sequences of Pro(60)-Ser(61) and Asp(154)-Asp(155)-Gly(156)-His(157) of Tia reported previously (Koide & Ikenaka (1973) Eur. J. Biochem. 32, 417-431) were amended to Ser(60)-Pro(61) and His(154)-Asp-Asp-Gly(157), respectively.     

Peptide mapping and assessment of cryoprotective activity of 26/27-kDa dehydrin from soybean seeds

To characterize the molecular weight diversity of seed dehydrin among soybean cultivars, 26/27-kDa soybean dehydrins were purified and compared in peptide mapping patterns, partial amino acid sequences, and cryoprotective activity on enzyme. In reverse phase chromatograms of their trypsin digests, we detected several distinctive peaks, one of which was attributed to a part of the internal glycine-rich region. Partial amino acid sequences of peptide fragments from trypsin and S. aureus V8 protease cleavage were found to be identical to the Mat9 translation. The CP50 of purified 26/27-kDa dehydrins were estimated to be 0.30 and 0.11 microM, respectively.     

Development of soybean aphid genomic SSR markers using next generation sequencing

Simple sequence repeats (SSRs) or microsatellites are very useful molecular markers, owing to their locus-specific codominant and multiallelic nature, high abundance in the genome, and high rates of transferability across species. The soybean aphid (Aphis glycines Matsumura) has become the most damaging insect pest of soybean (Glycine max (L.) Merr.) in North America, since it was first found in the Midwest of the United States in 2000. Biotypes of the soybean aphid capable of colonizing newly developed aphid-resistant soybean cultivars have been recently discovered. Genetic resources, including molecular markers, to study soybean aphids are severely lacking. Recently developed next generation sequencing platforms offer opportunities for high-throughput and inexpensive genome sequencing and rapid marker development. The objectives of this study were (i) to develop and characterize genomic SSR markers from soybean aphid genomic sequences generated by next generation sequencing technology and (ii) to evaluate the utility of the SSRs for genetic diversity or relationship analyses. In total 128 SSR primer pairs were designed from sequences generated by Illumina GAII from a reduced representation library of A. glycines. Nearly 94% (120) of the primer pairs amplified SSR alleles of expected size and 24 SSR loci were polymorphic among three aphid samples from three populations. The polymorphic SSRs were successfully used to differentiate among 24 soybean aphids from Ohio and South Dakota. Sequencing of PCR products of two SSR markers from four aphid samples revealed that the allelic polymorphism was due to variation in the SSR repeats among the aphids. These markers should be particularly useful for genetic differentiation among aphids collected from soybean fields at different localities and regions. These SSR markers provide the soybean aphid research community with the first set of PCR-based codominant markers developed from the genomic sequences of A. glycines.     

Development of microsatellite markers in potato and their use in phylogenetic and fingerprinting analyses.

Three genomic libraries were constructed using a mixture of DNA from Solanum phureja Juz. & Buk., and S. chacoense Bitt. Two of the libraries were enriched for ATT and GT repeats (a 27-fold enrichment was achieved). In total, 3500 clones of the conventional library, 1,000 of the library enriched for ATT, and 12,000 of the one enriched for GT were screened with five different repeat motifs, and a total of 18 primer pairs was obtained. Another group of 12 primer pairs was obtained from the SSR-containing sequences in the public databases (18 SSR-containing sequences were utilized). From among 30 newly developed primer pairs, 12 previously published ones, and 12 pairs developed for tomato, 7 were used to identify 12 different potato cultivars and introductions, and 12 were used to study phylogenetic distance among seven wild and cultivated potato species. Two SSR markers were sufficient to discriminate the 12 cultivars. The mean number of alleles per polymorphic locus was 5 for the 12 cultivars and 4.5 for the seven species. The results obtained in this study confirm those achieved in similar studies in other plant species regarding the abundance and use of SSR markers in identifying species and cultivars.     

Divergence and allelomorphic relationship of a soybean virus resistance gene based on tightly linked DNA microsatellite and RFLP markers

The use of genetically diverse resistance sources is important in breeding for durable disease resistance. Detection and evaluation of resistance genes by conventional inheritance experiments, however, often require laborious screening and genetic testing. In the present study, a marker-assisted screening for resistance sources was initiated in soybean [Glycine max (L.) Merr] using one DNA microsatellite and two RFLP markers tightly linked to a soybean mosaic virus (SMV) resistance gene (Rsv1). The three marker loci were used to screen 67 diverse soybean cultivars, breeding lines, and plant introductions. Five variants were found at the microsatellite locus (HSP176L), and the two RFLP loci (pA186 and pK644a) near Rsv1 show a remarkably higher level of restriction polymorphism than Rsv1-independent RFLP loci. Several specific variants at the three marker loci were found to be correlated with virus resistance, among which HSP176L-2 can be detected by PCR, thus may be useful for germplasm screening. The grouping of the 67 accessions according to their multilocus marker variants agrees with the available pedigree information. When all, or most, of the cultivars within a given group with the same Rsv1-linked marker variant are resistant, their SMV resistance is most likely conferred by Rsv1. These putatively Rsv1-carrying groups contain a total of 38 SMV-resistant lines including six differential cultivars that are known to carry Rsv1. The remaining seven resistant accessions (Columbia, Holladay, Peking, Virginia, FFR-471, PI 507403, and PI 556949) do not carry resistance marker variants, and at least some of them could be sources of resistance genes independent of Rsv1.