Study of the biosafety of genetically modified soybean in the center of its origin and diversity in the Far East of the Russian Federation

Wild soybean (Glycine soja Sieb. and Zucc.) is the nearest relative of cultivated soybean (G. max (L.) Merr.). Study of the population genetic structure of wild-growing relatives of genetically modified (GM) plants in the centers of their origin is one of the main procedures before introduction of GM crops in these areas. We studied the genetic variability of nine wild-growing soya populations of Primorskii krai using RAPD analysis. The level of genetic variability of G. soja was considerably higher than that of G. max. We analyzed phylogenetic relationships in the genus Glycine subgenus Soja using RAPD markers. Our data confirm the validity of allocation to G. gracilis of the rank of species.     

Evaluation of the geographical pattern of genetic diversity of Glycine soja and Glycine max based on four single copy nuclear gene loci: For conservation of soybean germplasm

Genetic diversity and its geographical patterns play a very important role in species conservation and exploitation. Here, nucleotide polymorphism patterns of four single copy nuclear gene loci in wild (Glycine soja) and cultivated soybean (Glycine max) populations from different geographical regions as well as their demographic history were analyzed. The results showed that: (1) Southern subpopulation has the highest, while central subpopulation revealed the lowest genetic diversity among three Chinese G. soja subpopulations. (2) Northern Chinese G. max subpopulation depicted higher genetic diversity than other two Chinese, Korean, Japanese and American G. max subpopulations. (3) Significant genetic differentiation (P < 0.001) was observed among Chinese G. soja subpopulations from three ecological zones. There was also a significant genetic differentiation（P < 0.01）between three Chinese and Japanese subpopulations of G. max. (4) The demographic dynamics revealed that effective population size of G. soja is expanding, while it was constant in G. max. G. soja is a useful germplasm resource to widen the genetic diversity of G. max. This study suggests that native populations of G. soja from different geo-ecological regions should be protected to conserve the genetic diversity.     

QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields

The objective of this study was to identify quantitative trait loci (QTL) affecting fitness of hybrids between wild soybean (Glycine soja) and cultivated soybean (Glycine max). Seed dormancy and seed number, both of which are important for fitness, were evaluated by testing artificial hybrids of G. soja × G. max in a multiple‐site field trial. Generally, the fitness of the F₁ hybrids and hybrid derivatives from self‐pollination was lower than that of G. soja due to loss of seed dormancy, whereas the fitness of hybrid derivatives with higher proportions of G. soja genetic background was comparable with that of G. soja. These differences were genetically dissected into QTL for each population. Three QTLs for seed dormancy and one QTL for total seed number were detected in the F₂ progenies of two diverse cross combinations. At those four QTLs, the G. max alleles reduced seed number and severely reduced seed survival during the winter, suggesting that major genes acquired during soybean adaptation to cultivation have a selective disadvantage in natural habitats. In progenies with a higher proportion of G. soja genetic background, the genetic effects of the G. max alleles were not expressed as phenotypes because the G. soja alleles were dominant over the G. max alleles. Considering the highly inbreeding nature of these species, most hybrid derivatives would disappear quickly in early self‐pollinating generations in natural habitats because of the low fitness of plants carrying G. max alleles.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

Uncovering signatures of selection in the soybean genome using SSR diversity near QTLs of agronomic importance

The cultivated soybean [Glycine max (L.) Merr.] is widely considered to descend from the wild soybean (G. soja Sieb. & Zucc.). This study was designed to evaluate the genetic variability and differentiation between G. soja and G. max, and to detect signatures of the selection that may have occurred during the domestication process from G. soja to G. max. A total of 192 G. soja accessions and 104 G. max accessions were genotyped using eight selected simple sequence repeat (SSR) markers assigned to three SSR groups. Four SSRs in group A were not located near any known QTL. Three SSRs in group B were associated with seed protein content, and an SSR in group C was associated with resistance to Sclerotinia stem rot. The number of alleles per locus and the level of genetic variability in G. soja were higher than those in G. max. A total of 122 out of 125 alleles were present in G. soja, but only 59 alleles were detected in G. max. The average gene diversity was 0.74 in G. soja and 0.64 in G. max. Four SSRs near QTLs of agronomic importance showed strong genetic differentiation and shift change in high frequency alleles in groups B and C between G. soja and G. max, revealing selection signatures that may reflect the domestication events and recent selective breeding. With reduced diversity in G. max, some undomesticated genes from G. soja should be prime candidates for introgression to increase the pool of diversity in G. max.     

A model to predict the frequency of integration of fitness-related QTLs from cultivated to wild soybean

With the proliferation of genetically modified (GM) products and the almost exponential growth of land use for GM crops, there is a growing need to develop quantitative approaches to estimating the risk of escape of transgenes into wild populations of crop relatives by natural hybridization. We assessed the risk of transgene escape by constructing a population genetic model based on information on fitness-related QTLs obtained from an F ₂ population of wild soybean G. soja × cultivated soybean Glycine max. Simulation started with ten F ₁ and 990 wild soybeans reproducing by selfing or outcrossing. Seed production was determined from the genetic effects of two QTLs for number of seeds (SN). Each seed survived winter according to the maternal genotype at three QTLs for winter survival (WS). We assumed that one neutral transgene was inserted at various sites and calculated its extinction rate. The presence of G. max alleles at SN and WS QTLs significantly decreased the probability of introgression of the neutral transgene at all insertion sites equally. The presence of G. max alleles at WS QTLs lowered the risk more than their presence at SN QTLs. Although most model studies have concentrated only on genotypic effects of transgenes, we show that the presence of fitness-related domestication genes has a large effect on the risk of transgene escape. Our model offers the advantage of considering the effects of both domestication genes and a transgene, and they can be widely applied to other wild × crop relative complexes.     

Amplified fragment length polymorphism (AFLP) in soybean: species diversity, inheritance, and near-isogenic line analysis

Amplified fragment length polymorphism (AFLP) analysis is a PCR-based technique capable of detecting more than 50 independent loci in a single PCR reaction. The objectives of the present study were to: (1) assess the extent of AFLP variation in cultivated (Gycine max L. Merr.) and wild soybean (G. soja Siebold & Zucc.), (2) determine genetic relationships among soybean accessions using AFLP data, and (3) evaluate the usefulness of AFLPs as genetic markers. Fifteen AFLP primer pairs detected a total of 759 AFLP fragments in a sample of 23 accessions of wild and cultivated soybean, with an average of 51 fragments produced per primer pair per accession. Two-hundred and seventy four fragments (36% of the total observed) were polymorphic, among which 127 (17%) were polymorphic in G. max and 237 (31%) were polymorphic in G. soja. F₂ segregation analysis of six AFLP fragments indicated that they segregate as stable Mendelian loci. The number of polymorphic loci detected per AFLP primer pair in a sample of 23 accessions ranged from 9 to 27. The AFLP phenotypic diversity values were greater in wild than in cultivated soybean. Cluster and principal component analyses using AFLP data clearly separated G. max and G. soja accessions. Within the G. max group, adapted soybean cultivars were tightly clustered, illustrating the relatively low genetic diversity present in cultivated soybean. AFLP analysis of four soybean near-isogenic lines (NILs) identified three AFLP markers putatively linked to a virus resistance gene from two sources. The capacity of AFLP analysis to detect thousands of independent genetic loci with minimal cost and time requirements makes them an ideal marker for a wide array of genetic investigations.     

Influence of age structure of plant populations on damping-off epidemics

Summary Germination synchrony may facilitate damping-off epidemics by creating a high density of uniformly susceptible individuals. We tested the hypothesis that synchronous germination causes increased seed and seedling mortality from damping-off in two legume species attacked by the fungal pathogen, Pythium aphanidermatum. Glycine max exhibited rapid, synchronous germination compared to its progenitor, G. soja, and suffered greater mortality from both pre-and postemergent damping-off in controlled environment experiments. However, when mixed-aged populations of G. max were created experimentally by staggering planting times, a significant increase in damping-off mortality occurred. In G. soja, which typically has mixed-aged populations due to asynchronous germination, experimental populations with an even-aged distribution also suffered increased damping-off mortality. Hence, the relationship between population age structure and damping-off mortality was species-specific. We propose that species differences in the duration of individual seedling susceptibility to disbase interact with population age structure to control the cutcome of damping-off epidemics.     

Polymorphism of the IPK1 gene among members of the genus Glycine

The structural polymorphism of the IPK1 gene involved in the biosynthesis of phytic acid was studied in 30 samples from three species of the genus Glycine (G. max, G. soja, and G. gracilis). Single nucleotide polymorphisms in the sequences of the first and the second introns of the IPK1 gene were found, differences in the length of microsatellite loci located in intron 2 were demonstrated, and a total of six polymorphic variants were identified. Data on the distribution of the identified alleles among the samples from the working set and the analysis and discussion of these data are given.     

Identification of putative QTL that underlie yield in interspecific soybean backcross populations

Glycine soja, the wild progenitor of soybean, is a potential source of useful genetic variation in soybean improvement. The objective of our study was to map quantitative trait loci (QTL) from G. soja that could improve the crop. Five populations of BC₂F₄-derived lines were developed using the Glycine max cultivar IA2008 as a recurrent parent and the G. soja plant introduction (PI) 468916 as a donor parent. There were between 57 and 112 BC₂F₄-derived lines in each population and a total of 468 lines for the five populations. The lines were evaluated with simple sequence repeat markers and in field tests for yield, maturity, plant height, and lodging. The field testing was done over 2 years and at two locations each year. Marker data were analyzed for linkage and combined with field data to identify QTL. Using an experimentwise significance threshold of P=0.05, four yield QTL were identified across environments on linkage groups C2, E, K, and M. For these yield QTL, the IA2008 marker allele was associated with significantly greater yield than the marker allele from G. soja. In addition, one lodging QTL, four maturity QTL, and five QTL for plant height were identified across environments. Of the 14 QTL identified, eight mapped to regions where QTL with similar effects were previously mapped. Many regions carrying the yield QTL were also significant for other traits, such as plant height and lodging. When the significance threshold was reduced and the data were analyzed with simple linear regression, four QTL with a positive allele for yield from G. soja were mapped. One epistatic interaction between two genetic regions was identified for yield using an experimentwise significance threshold of P=0.05. Additional research is needed to establish whether multiple trait associations are the result of pleiotropy or genetic linkage and to retest QTL with a positive effect from G. soja.Communicated by H.C. Becker     

Genetic diversity and linkage disequilibrium studies on a 3.1-Mb genomic region of chromosome 3B in European and Asian bread wheat (Triticum aestivum L.) populations

Genetic diversity and linkage disequilibrium (LD) were investigated in 376 Asian and European accessions of bread wheat (Triticum aestivum L.). After a first and rapid screening about diversity and genetic structure at the whole genome scale using 70 simple sequence repeats (SSRs), we focused on a sequenced contig (ctg954) of 3.1 Mb located on the short arm of chromosome 3B of cv. Chinese Spring, using 32 SSRs and 10 single nucleotide polymorphisms. This contig is part of a multiple fungal resistance region. Mean polymorphism information content value on the 32 SSRs was slightly higher in the Asian genepool (0.396) than that for the European (0.329) pool. Compared with results at the whole genome scale, data from this 3.1-Mb region indicated similar trends in genetic diversity indices between both genepools. Population structure and molecular variance analyses demonstrated significant genetic differentiation and geographical subdivision in both groups of accessions. Concerning LD at the contig level, the European population had a significantly higher mean r ² value (0.23) than the Asian population (0.18), indicating a stronger LD in the European material. With a mean of 1 marker every 74 kb, the resolution reached here allowed to perform a detailed comparative analysis of the LD and genetic diversity along the complete 3.1-Mb region in both genepools. A sliding-window approach revealed some interesting regions of the contig where LD is increasing when genetic diversity is decreasing. This study provides an in-depth understanding of molecular population genetics in European and Asian wheat gene pools, and prospects for association mapping of important sources of fungal disease resistance.     

Study of interactions between wild soybean subpopulations (<Emphasis Type="Italic">Glycine soja</Emphasis>) in the valley of the Tsukanovka river in the south of Russia’s Far East

Comparative analysis of the genetic structure of natural and anthropogenic populations of G. soja gives significant information about the formation of different populations and allows for the developing of measures for the preservation of the unique natural gene bank of wild soybean, which is a species closely related to cultivated soybean. In this study, ISSR-markers were used to carry out a comparative analysis of the genetic structure of natural and anthropogenic subpopulations of G. soja for studying the possible mutual influence of subpopulations of anthropogenic and natural phytocenosis on the formation of their genetic diversity and studying the genetic structure of natural subpopulations of wild soybean in the contact places between the two types of cenoses. As a result, the characteristics that describe the genetic diversity of the studied populations have been identified, and the important role of interaction between subpopulations of different phytocenoses in the formation of the spatial genetic structure of the population in the Tsukanovka river valley have been demonstrated.     

Population genetic structure,migration, and polyploidy origin of a medicinal species Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum, a member of family Cucurbitaceae, is a perennial creeping herb used as a traditional medicinal plant in China. In this study, six polymorphic nSSR and four EST‐SSR primers were used to genotype 1,020 individuals in 72 wild populations of G. pentaphyllum. The genetic diversity and population structure were investigated, and ecological niche modeling was performed to reveal the evolution and demographic history of its natural populations. The results show that G. pentaphyllum has a low level of genetic diversity and high level of variation among populations because of pervasive asexual propagation, genetic drift, and long‐term habitat isolation. Results of the Mantel test demonstrate that the genetic distance and geographic distance are significantly correlated among G. pentaphyllum natural populations. The populations can be divided into two clusters on the basis of genetic structure. Asymmetrical patterns of historical gene flow were observed among the clusters. For the contemporary, almost all the bidirectional gene flow of the related pairs was symmetrical with slight differences. Recent bottlenecks were experienced by 34.72% of the studied populations. The geographic range of G. pentaphyllum continues to expand northward and eastward from Hengduan Mountains. The present distribution of G. pentaphyllum is a consequence of its complex evolution. Polyploidy in G. pentaphyllum is inferred to be polygenetic. Finally, G. pentaphyllum is a species in need of protection, so in situ and ex situ measures should be considered in the future.     

Genetic divergence between the South Korean and Mongolian populations of the dung beetle,Gymnopleurus mopsus (Coleoptera: Scarabaeidae) based on mitochondrial cytochrome c oxidase subunit I (COI) gene sequences

The locally extinct dung beetle, Gymnopleurus mopsus Pallas, 1781 (Coleoptera: Scarabaeidae), has not been found in South Korea since the 1970s. This research was conducted to understand the genetic divergence between the South Korean and Mongolian populations of G. mopsus as a part of its reintroduction program in South Korea. The genetic distance and diversity were determined using the mitochondrial cytochrome c oxidase subunit I (COI) gene sequence (658 bp ) corresponding to the DNA barcode region. The mitochondrial COI gene sequences of 64 individuals of G. mopsus collected in South Korea (7 individuals) and Mongolia (57 individuals) showed a relatively high level of genetic diversity (nucleotide diversity, 0.0078 ± 0.0007; Haplotype diversity, 0.965 ± 0.017). The genetic distances between the South Korean and Mongolian populations lay within the intraspecific level and the phylogenetic reconstruction using the neighbor‐joining (NJ) method showed that all individuals belonged to a single clade. This result indicates that the current Mongolian population of G. mopsus is a good candidate source population to restore the locally extinct population of the species in South Korea.     

Investigation of the Gracilaria gracilis (Gracilariales,Rhodophyta) proteome response to nitrogen limitation

Inorganic nitrogen has been identified as the major growth‐limiting nutritional factor affecting Gracilaria gracilis populations in South Africa. Although the physiological mechanisms implemented by G. gracilis for adaption to low nitrogen environments have been investigated, little is known about the molecular mechanisms of these adaptions. This study provides the first investigation of G. gracilis proteome changes in response to nitrogen limitation and subsequent recovery. A differential proteomics approach employing two‐dimensional gel electrophoresis and liquid chromatography–tandem mass spectrometry was used to investigate G. gracilis proteome changes in response to nitrogen limitation and recovery. The putative identity of 22 proteins that changed significantly (P < 0.05) in abundance in response to nitrogen limitation and recovery was determined. The identified proteins function in a range of biological processes including glycolysis, photosynthesis, ATP synthesis, galactose metabolism, protein‐refolding and biosynthesis, nitrogen metabolism and cytoskeleton remodeling. The identity of fructose 1,6 biphosphate (FBP) aldolase was confirmed by western blot analysis and the decreased abundance of FBP aldolase observed with two‐dimensional gel electrophoresis was validated by enzyme assays and western blots. The identification of key proteins and pathways involved in the G. gracilis nitrogen stress response provide a better understanding of G. gracilis proteome responses to varying degrees of nitrogen limitation and is the first step in the identification of biomarkers for monitoring the nitrogen status of cultivated G. gracilis populations.     

Assessing the effects of rare alleles and linkage disequilibrium on estimates of genetic diversity in the chicken populations

Phenotypic diversity in poultry has been mainly driven by artificial selection and genetic drift. These led to the adaptation to the environment and the development of specific phenotypic traits of chickens in response to their economic use. This study evaluated genetic diversity within and between Russian breeds and populations using Illumina Chicken 60 K SNP iSelect BeadChip by analysing genetic differences between populations with Hudson's fixation index (F_ST statistic) and heterozygosity. We estimated the effect of rare alleles and linkage disequilibrium (LD) on these measurements. To assess the effect of LD on the genetic diversity population, we carried out the LD-based pruning (LD < 0.5 and LD < 0.1) for seven chicken populations combined (I) or separately (II). LD pruning was specific for different dataset groups. Because of the noticeably large sample size in the Russian White RG population, pruning was substantial for Dataset I, and F_ST values were only positive when LD < 0.1 pruning was applied. For Dataset II, the LD pruning results were confirmed by examining heterozygosity and alleles' frequency distribution. LD between single nucleotide polymorphisms was consistent across the seven chicken populations, except the Russian White RG population with the smallest r² values and the largest effective population size. Our findings suggest to study variability in each population LD pruning has to be carried separately not after merging to avoid bias in estimates.     

QTL associated with horizontal resistance to soybean cyst nematode in Glycine soja PI464925B

Soybean cyst nematode (Heterodera glycines Ichinohe; SCN) is the primary disease responsible for yield loss of soybean [Glycine max (L.) Merr.]. Resistant cultivars are an effective management tool; however, the sources currently available have common resistant genes. Glycine soja Sieb. and Zucc., the wild ancestor of domesticated soybean, represents a diverse germplasm pool with known SCN resistance. The objectives of this research were to: (1) determine the genetic variation and inheritance of SCN resistance in a G. max (‘S08-80’) × G. soja (PI464925B) F _4:5 recombinant inbred line (RIL) population; and (2) identify and evaluate quantitative trait loci (QTL) associated with SCN resistance. Transgressive segregation for resistance was observed, although neither parent was resistant to the Chatham and Ruthven SCN isolates. Broad sense heritability was 0.81 for the Ruthven and 0.91 for the Chatham isolate. Root dry weight was a significant covariate that influenced cyst counts. One RIL [female index (FI) = 5.2 ± 1.11] was identified as resistant to the Chatham isolate (FI < 10). Seventeen and three RILs infected with Chatham and Ruthven isolates, respectively, had mean adjusted cyst counts of zero. Unique and novel QTL, which derived resistance from G. soja, were identified on linkage groups I, K, and O, and individually explained 8, 7 and 5% (LOD = 2.1–2.7) of the total phenotypic variation, respectively. Significant epistatic interactions were found between pairs of SSR markers that individually may or may not have been associated with SCN resistance, which explained between 10 and 15% of the total phenotypic variation. Best-fit regression models explained 21 and 31% of the total phenotypic variation in the RIL population to the Chatham and Ruthven isolates, respectively. The results of this study help to improve the understanding of the genetic control of SCN resistance in soybean caused by minor genes resulting in horizontal resistance. The incorporation of the novel resistance QTL from G. soja could increase the durability of SCN-resistance in soybean cultivars, especially if major gene resistance breaks down.     

The salt resistance of wild soybean (Glycine soja Sieb. et Zucc. ZYD 03262) under NaCl stress is mainly determined by Na+ distribution in the plant

Understanding the salt resistance mechanism of wild soybean is important in improving salt tolerance of cultivated soybean. Therefore, we comparatively analyzed effects of NaCl on photosynthesis, antioxidant enzyme activity, and ion distribution in a cultivated (Glycine. max) and a wild (Glycine soja) soybean to study the salt resistance mechanism of the G. soja. The results showed that more Na⁺ was accumulated in the G. soja roots than in the G. max roots, but the Na⁺ in the G. soja leaves was much less than that observed in the G. max leaves. The Na⁺ concentrations in the G. soja leaves were not high enough to affect the photosynthetic apparatus, which was demonstrated by less inhibition of photosynthetic activity, stomatal conductance, carboxylation efficiency in the G. soja leaves than in the G. max leaves after treated with different concentrations of NaCl. Meanwhile, there were no significant changes in intercellular CO₂ concentration, maximum PSII quantum yield, and relative water content in the G. soja leaves after NaCl treatment, while they significantly decreased in the G. max leaves. The non-photochemical quenching and the activities of superoxide dismutase (EC 1.15.1.1) and ascorbate peroxidase (EC 1.11.1.11) in the G. soja leaves increased with the increasing of NaCl concentrations, whereas only the activity of superoxide dismutase increased in G. max leaves. Based on these results, we suggested that the G. soja is able to accumulate higher levels of Na⁺ in its roots, and prevent the transportation of Na⁺ to leaves to protect photosynthetic apparatus from salt damage.     

The genomic relationship between Glycine max (L.) Merr. and G. soja Sieb. and Zucc. as revealed by pachytene chromosome analysis

Summary This study was conducted with the objective of determining the genomic relationship between cultivated soybean (Glycine max) and wild soybean (G. soja) of the subgenus Soja, genus Glycine. Observations on cross-ability rate, hybrid viability, meiotic chromosome pairing, and pollen fertility in F ₁ hybrids of G. max × G. soja and reciprocals elucidated that both species hybridized readily and set mature putative hybrid pods, generated vigorous F₁ plants, had a majority of sporocytes that showed 18II + 1IV chromosome association at diakinesis and metaphase I, and had a pollen fertility that ranged from 49.2% to 53.3%. A quadrivalent was often associated with the nucleolus, suggesting that one of the chromosomes involved in the interchange is a satellited chromosome. Thus, G. max and G. soja genetic stocks used in this study have been differentiated by a reciprocal translocation. Pachytene analysis of F₁ hybrids helped construct chromosome maps based on chromosome length and euchromatin and heterochromatin distribution. Chromosomes were numbered in descending order of 1–20. Pachytene chromosomes in soybean showed heterochromatin distribution on either side of the centromeres. Pachytene analysis revealed small structural differences for chromosomes 6 and 11 which were not detected at diakinesis and metaphase I. This study suggests that G. max and G. soja carry similar genomes and validates the previously assigned genome symbol GG.Research supported in part by the Illinois Agricultural Experiment Station and U.S. Department of Agriculture Competitive Research Grant (85-CRCR-1-1616)