Distribution of restriction site polymorphism within the chloroplast genome of the genus Glycine, subgenus Soja

Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Pinaceae, the largest family of conifers, has a diversified organization of chloroplast (cp) genomes with two typical highly reduced inverted repeats (IRs). In the current study, we determined the complete sequence of the cp genome of an economically and ecologically important conifer tree, the loblolly pine (Pinus taeda L.), using Illumina paired-end sequencing and compared the sequence with those of other pine species. The results revealed a genome size of 121,531 base pairs (bp) containing a pair of 830-bp IR regions, distinguished by a small single copy (42,258 bp) and large single copy (77,614 bp) region. The chloroplast genome of P. taeda encodes 120 genes, comprising 81 protein-coding genes, four ribosomal RNA genes, and 35 tRNA genes, with 151 randomly distributed microsatellites. Approximately 6 palindromic, 34 forward, and 22 tandem repeats were found in the P. taeda cp genome. Whole cp genome comparison with those of other Pinus species exhibited an overall high degree of sequence similarity, with some divergence in intergenic spacers. Higher and lower numbers of indels and single-nucleotide polymorphism substitutions were observed relative to P. contorta and P. monophylla, respectively. Phylogenomic analyses based on the complete genome sequence revealed that 60 shared genes generated trees with the same topologies, and P. taeda was closely related to P. contorta in the subgenus Pinus. Thus, the complete P. taeda genome provided valuable resources for population and evolutionary studies of gymnosperms and can be used to identify related species.

Comparative analysis of complete plastid genomes from wild soybean (Glycine soja) and nine other Glycine species

The plastid genomes of different plant species exhibit significant variation, thereby providing valuable markers for exploring evolutionary relationships and population genetics. Glycine soja (wild soybean) is recognized as the wild ancestor of cultivated soybean (G. max), representing a valuable genetic resource for soybean breeding programmes. In the present study, the complete plastid genome of G. soja was sequenced using Illumina paired-end sequencing and then compared it for the first time with previously reported plastid genome sequences from nine other Glycine species. The G. soja plastid genome was 152,224 bp in length and possessed a typical quadripartite structure, consisting of a pair of inverted repeats (IRa/IRb; 25,574 bp) separated by small (178,963 bp) and large (83,181 bp) single-copy regions, with a 51-kb inversion in the large single-copy region. The genome encoded 134 genes, including 87 protein-coding genes, eight ribosomal RNA genes, and 39 transfer RNA genes, and possessed 204 randomly distributed microsatellites, including 15 forward, 25 tandem, and 34 palindromic repeats. Whole-plastid genome comparisons revealed an overall high degree of sequence similarity between G. max and G. gracilis and some divergence in the intergenic spacers of other species. Greater numbers of indels and SNP substitutions were observed compared with G. cyrtoloba. The sequence of the accD gene from G. soja was highly divergent from those of the other species except for G. max and G. gracilis. Phylogenomic analyses of the complete plastid genomes and 76 shared genes yielded an identical topology and indicated that G. soja is closely related to G. max and G. gracilis. The complete G. soja genome sequenced in the present study is a valuable resource for investigating the population and evolutionary genetics of Glycine species and can be used to identify related species.

Genome re-sequencing of semi-wild soybean reveals a complex Soja population structure and deep introgression

Semi-wild soybean is a unique type of soybean that retains both wild and domesticated characteristics, which provides an important intermediate type for understanding the evolution of the subgenus Soja population in the Glycine genus. In this study, a semi-wild soybean line (Maliaodou) and a wild line (Lanxi 1) collected from the lower Yangtze regions were deeply sequenced while nine other semi-wild lines were sequenced to a 3-fold genome coverage. Sequence analysis revealed that (1) no independent phylogenetic branch covering all 10 semi-wild lines was observed in the Soja phylogenetic tree; (2) besides two distinct subpopulations of wild and cultivated soybean in the Soja population structure, all semi-wild lines were mixed with some wild lines into a subpopulation rather than an independent one or an intermediate transition type of soybean domestication; (3) high heterozygous rates (0.19-0.49) were observed in several semi-wild lines; and (4) over 100 putative selective regions were identified by selective sweep analysis, including those related to the development of seed size. Our results suggested a hybridization origin for the semi-wild soybean, which makes a complex Soja population structure.

Fine-scale phylogenetic structure and major events in the history of the current wild soybean (Glycine soja) and taxonomic assignment of semi-wild type (Glycine gracilis Skvortz.) within the Chinese subgenus Soja

Wild and cultivated species of soybeans have coexisted for 5000 years in China. Despite this long history, there is very little information on the genetic relationship of Glycine soja and G. max. To gain insight into the major events in the history of the subgenus Soja, we examined 20 simple sequence repeat (SSR) markers of a large number of accessions (910). The results showed no significant differences between wild and semi-wild soybeans in genetic diversity but significant differences between G. soja and G. max. Ancestry and cluster analyses revealed that semi-wild soybeans should belong to the wild category and not to G. max. Our results also showed that differentiation had occurred not only among G. soja, G. gracilis, and G. max but also within G. soja and within G. gracilis. Glycine soja had 3 clear genetic categories: typical small-seeded (≤2.0 g 100-seed weight), dual-origin middle-seeded (2.0-2.5 g), and large-seeded plants (2.51-3.0 g). These last were genetically close to G. gracilis, their defining some traits having been acquired mainly by introgression from soybeans. Small-seeded G. gracilis (3.01-3.5 g) were genetically different from larger seeded ones (from 3.51 to 4.0 to over 10 g). Seed size predominated over seed coat color in evolutionary degree. Typical and large-seeded G. soja were found to have 0.7% and 12% introgressive cultivar genes, respectively. The genetic boundary of G. gracilis was at the range of 2.51-3.0 g of G. soja. In the great majority of wild accessions, traits such as white flowers, gray pubescences, no-seed bloom, and colored seed coats were likely introgressive from domesticated soybeans.

The origin and fate of morphological intermediates between wild and cultivated soybeans in their natural habitats in Japan

The spread of transgenes into the genome of wild soybean is a concern when transgenic and wild soybeans are planted sympatrically. The objectives of this study were to investigate the origin and fate of morphological intermediates between wild and cultivated soybeans in their natural habitats in Japan. Twenty nuclear microsatellite and two chloroplast dCAPS markers were used to evaluate genetic variation of 468 wild, 17 intermediate, and 12 cultivated soybean samples collected from six sites between 2003 and 2006. Allelic differentiation of microsatellite markers between wild and cultivated soybeans was sufficient to detect their hybrids. Based on levels of observed heterozygosity, intermediate soybean plants were from two generations: either F(1) or an early segregating generation. Genetic admixture analysis and parentage assignment analysis revealed that the parents of all intermediate soybean plants could be assigned to a particular wild soybean plant and late-maturing cultivar. The chloroplast DNA haplotypes revealed that all intermediate soybean plants originated from gene flow from cultivated to wild soybeans at all sites. Based on monitoring at both the phenotypic and molecular levels, hybrids quickly disappeared from natural habitats, and secondary gene flow from these plants to wild soybean was not detected. Thus, while gene flow from transgenic soybean into wild soybean can occur, gene introgression appears to be rare in natural habitats in Japan. This is the first report on the detection of gene flow from cultivated to wild soybean at the molecular level.

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.

Phylogenetic relationships of the mitochondrial genomes in the genus Glycine subgenus Soja

Restriction fragment length polymorphisms (RFLPs) of mitochondrial DNA (mtDNA) of wild and cultivated soybeans were analyzed to study their phylogenetic relationships. The observed number of differences in hybridization profiles greatly varied (1-17 patterns) with both the mtDNA probes and the restriction enzymes that were used. A cladistic analysis was conducted based on the RFLP data. In the parsimonious tree, four distinct groups appeared among 20 accessions of the subgenus Soja representing 20 mitochondrial genome types. Common features with regard to geographic distributions in natural populations in East Asia were observed among the mitochondrial genome types of wild soybean that belonged to the same group: one clade consisted of genome types IIg and VIIg that are detected with very low frequencies; another clade consisted of genome types Ic, Id, Ie, and Ik whose distributions are highly biased mainly in Japan. The genome types that are widely distributed in East Asia such as IVa, IVb, and Va were not grouped into the same clade. The mitochondrial genome types IIIb, IVb, and IVc, in which two different chloroplast genome types exist, belonged to the same clade. Possible changes in mitochondrial genomes during the expansion of the distribution of wild soybeans in East Asia were discussed.

The nuclear ribosomal DNA intergenic spacers of wild and cultivated soybean have low variation and cryptic subrepeats

The intergenic ribosomal DNA spacers (IGSs) from cultivated soybean (Glycine max) and wild soybean (Glycine soja) were sequenced and compared with six other legumes. These IGS sequences were 1821 bp in length in G. soja and G. max cultivars Arksoy, Ransom, and Tokyo, and 1823 bp long in the G. max cultivar Columbus. These represent the smallest published plant IGS sequences to date. Two clones from each of the above five cultivars were sequenced and only 22 sites (1.2%) were polymorphic, thereby supporting previous work that showed low genetic variation in cultivated soybean. The amount of variation observed between different clones derived from the same individual was equal to the amount seen between different cultivars. The soybean IGS sequence was aligned with six other published legume sequences and two homologous regions were identified. The first spans positions 706-1017 in the soybean IGS sequence and ends at a putative promoter site that appears conserved among all legumes. The second is located within the 5' external transcribed spacer, spans positions 1251-1823 in soybean, and includes sequences first identified as subrepeats IV-1 and IV-2 in Vicia angustifolia. Sequences homologous to these two subrepeats were identified among all legume species examined and are here designated "cryptic subrepeats" (CS-1, CS-2) given the range in similarity value (79-96% for CS-1 and 60-95% for CS-2). Comparisons of CS-1 and CS-2 sequences within individual species show that divergence (substitutional mutations, insertions, and deletions) is sufficiently high to obscure recognition of the repeat nature of these sequences by routine dot plot analytical methods. The lack of subrepeats in the 5' half of the soybean IGSs raises questions regarding the role they play in transcription termination or enhancement.Key words: Glycine, Fabaceae, ribosomal DNA, promoter.

Sequence variation of chloroplast DNA that involves EcoRI and ClaI restriction site polymorphisms in soybean

Restriction fragment length polymorphisms (RFLPs) of EcoRI-and ClaI-digested chloroplast DNA (cpDNA) within the genus Glycine subgenus Soja were characterized. Two mutations were found to be responsible for the EcoRI and ClaI restriction site polymorphisms, and both were located in a region in which many ribosomal protein genes are clustered. This region is within the large single copy region of cpDNA and is located close to an inverted repeat. The locations of restriction sites of EcoRI and ClaI in the cpDNA region were analyzed by DNA gel-blot analyses and PCR amplification, which were followed by sequencing analyses. The EcoRI site polymorphism was found to have occurred in the intergenic spacer between rps11 and rpl36, while the ClaI site polymorphism was located within the 3' part of the coding region of rps3. The mutations that cause EcoRI and ClaI polymorphisms were both found to be single base substitutions. In addition to these polymorphisms, novel sequence variations in soybean cpDNA were detected near the sites of these mutations. Previously, it was shown that cultivated soybeans could be classified into three groups (I, II, and III) based on their cpDNA RFLPs. A comparison of the cpDNA sequences of soybeans in the present study was consistent with the notion that the cpDNA of group II soybeans is an intermediate between the cpDNAs of groups I and II.

A single nuclear locus phylogeny of soybean based on DNA sequence

Soybean [Glycine max (L.) Merr.] evolution was examined by sequencing portions of the restriction fragment length polymorphism (RFLP) locus A-199a of 21 taxa from the Glycininae and 1 from the Phaseoleae. Four hundred nucleotides were determined in each, aligned, and then compared for these taxa. Within the annual soybean subgenus (Soja), the four accessions differed at as many as 2.2% of the nucleotides. Among 13 perennial soybean species (subgenus Glycine), nucleotide variation ranged from 1.7% to 8.4%. The nucleotide difference between the two soybean subgenera was 3.0-7.0%. Nucleotide variation between the genus Glycine and the related genera of Neonotonia, Amphicarpa, Teramnus, and Phaseolus ranged from 8.2% to 16.4%. In addition to nucleotide substitutions, insertions/deletions (indels) differences were also observed and were consistent with nucleotide-based analysis. Cladistic analysis of the A-199a sequences was performed using Wagner parsimony to construct a soybean phylogeny. Sixteen equally parsimonious trees were produced from these data. The trees were 246 steps in length with a consistency index of 0.78. Indels distribution upon the consensus topology revealed a pattern congruent with the nucleotide-based phylogeny. The current taxonomic status of the soybean subgenera and the related genera of Neonotonia, Amphicarpa, and Teramnus were well-supported and appear monophyletic in this analysis. Homoplasy within the subgenus Glycine led to a lack of resolved topology for many of these 13 taxa. However, the Glycine clade topology was consistent with phylogenies proposed using crossing experiments and cpDNA RFLPs. These genera were arranged from ancestral to derived as: Teramnus, Amphicarpa, Neonotonia, and Glycine when Phaseolus vulgaris was used as an outgroup.

Chloroplast DNA as an evolutionary marker in thePhaseolus vulgaris complex

We have analyzed the changes occurring in the chloroplast DNA (cpDNA) of taxa belonging to thePhaseolus vulgaris complex to help clarify relationships among species of this complex. Two restriction maps for 11 restriction enzymes comprising the whole chloroplast genome from a wildP. vulgaris and a wildP. coccineus accession were constructed. These maps allowed us to compare a total of 330 restriction sites between the two genomes in order to identify polymorphisms, assess the type of mutations detected, and identify regions of high variability. A region, located in the large single-copy region near the borders with the inverted repeats, accounted for a large portion of the variation. Most of the mutations detected were due to restriction sites gains or losses. Variable and conserved regions were then evaluated in 30 accessions belonging to taxa of theP. vulgaris complex. Phylogenetic analyses were made using parsimony methods. Conclusions obtained from such analyses were the following: (1) there was high cpDNA variability withinP. coccineus but not inP. vulgaris. (2)P. coccineus subsp.glabellus showed a very distinct cpDNA type that strongly suggests that it actually belongs to a different but as yet undetermined section of the genus. Our cpDNA observations are supported by distinctive morphological traits and reproductive biology of this taxon. (3) InP. coccineus subsp.darwinianus (also classified asP. polyanthus), the cpDNA lineage was in disagreement with data obtained from nuclear markers and suggested a reticulated origin by hybridization betweenP. coccineus as the male parent and an ancestralP. polyanthus type, closely allied toP. vulgaris, as the seed parent. This initial cross was presumably followed by repeated backcrossing toP. coccineus. Our cpDNA studies illustrate the importance of molecular markers in elucidating phylogenetic relationships. They also indicate that accurate phylogenies will require analyses of both nuclear and cytoplasmic genomes.

Chloroplast DNA evidence for non-random selection of females in an outcrossed population of soybeans [Glycine max (L.)]

Restriction fragment length polymorphisms (RFLPs) were used to assess chloroplast DNA (cpDNA) variation in a population of soybeans subjected to continuous cycles of forced outcrossing. This population was derived by crossing 39 female lines with four male-sterile (Ms2ms2) maintainer lines and advancing each generation by selecting only outcrossed seed borne on male-sterile (ms2ms2) plants. Analysis of the original 39 female lines revealed three groups based on cpDNA RFLPs. These three groups had been previously documented in soybeans, and the distribution of these groups among the female parents of this population was similar to that observed in germ plasm surveys of soybean. Thirty-four of the female parents had group I cpDNA, 3 had group II, and 2 had group III. Plants collected from this population after seven cycles of outcrossing were scored for four morphological traits (flower color, pubescence color, seed color, and pubescence type) known to be controlled by alleles at single nuclear loci. The frequencies of the phenotypes observed in this study indicated that the population underwent random mating with respect to flower and pubescence color, but deviated from random mating at the other two loci. Analysis of 158 of these same plants collected from the population after seven cycles of outcrossing revealed no individuals with group II or group III cpDNAs. The fixation of the group I cpDNA marker in this outcrossing population was judged to result primarily from selection against individuals in the population with the rare cpDNAs.

Chloroplast DNA restriction fragment length polymorphism in Sequoia sempervirens D. Don Endl., Pseudotsuga menziesii (Mirb.) Franco, Calocedrus decurrens (Torr.), and Pinus taeda L

The extent and type of chloroplast DNA restriction fragment length polymorphism was determined among individual tree samples of coast redwood, Douglas fir, incense-cedar, and loblolly pine. A total of 107 trees was surveyed for three restriction enzymes (BamHI, EcoRI, HindIII) and six chloroplast DNA probes from petunia (P3, P4, P6, P8, P10, S8). The probes comprise 64% of the petunia chloroplast genome. Polymorphisms were detected in all species but loblolly pine. Coast redwood and incense-cedar had a small number of rare variants, whereas Douglas fir had one highly polymorphic region of insertions/deletions in sequences revealed by the P6 probe from petunia. The mutation hotspot is currently being studied by DNA sequence analysis.