Distribution of new satellites and simple sequence repeats in annual and perennial Glycine species

The repeat sequences occupied more than 50 % of soybean genome. In order to understand where these repeat sequences distributed in soybean genome and its related Glycine species, we examined three new repeat sequences—soybean repeat sequence (SBRS1, SBRS2 and SBRS3), some nonspecific repeat sequences and 45S rDNA on several Glycine species, including annual and perennial accessions in this study. In the annual species, G. soja, signals for SBRS1 and ATT repeat can be found on each chromosome in GG genome, but those for SBRS2 and SBRS3 were located at three specific loci. In perennial Glycine species, these three SBR repeat frequently co-localized with 45S rDNA, two major 45S rDNA loci were found in all tetraploid species. However, an extra minor locus was found in one accession of the G. pescadrensis (Tab074), but not in another accession (Tab004). We demonstrate that some repetitive sequences are present in all Glycine species used in the study, but the abundancy is different in annual or perennial species. We suggest this study may provide additional information in investigations of the phylogeny in the Glycine species. Electronic supplementary material The online version of this article (doi:10.1186/s40529-015-0103-9) contains supplementary material, which is available to authorized users.


Background
Soybean (Glycine max L. Merrill) is the most important legume crop because of its economic value and biological features for research. Since soybean seed consists of high protein and oil contents, it has become one of the most important food crops for human and livestock, as well as a valuable resource for biofuel energy. Improving soybean varieties has been slow and limited because of the lack of genetic diversity (Keim et al. 1990; Krishnan et al. 2001). Wild germplasm may be a significant breeding resource for current varieties that have undergone several genetic bottlenecks and have shown limited genetic variability (Hyten et al. 2006;Keim et al. 1990). These wild germplasm harbor many traits of interest to soybean breeders, such as biotic and abiotic disease resistance (Newell and Hymowitz 1975), high seed protein content (Sebolt et al. 2000) and increased yield (Concibido et al. 2003). Thus, investigating Glycine species is important for soybean genetics and breeding.
There are around 30 species in Genus Glycine which had been classified into two subgenera, Soja and Glycine. The subgenus Soja (2n = 40) includes only G. max and its annual close related G. soja Sieb. and Zucc., which may intercross freely with soybean (Palmer et al. 1987;Singh and Hymowitz 1988). The subgenus Glycine contains all perennial species with various genome types and ploidies (2n = 2x = 40, 2n = 4x = 80 and some aneuploids) (Doyle et al. 2004;Grant et al. 1984;Newell and Hymowitz 1983;Singh and Hymowitz 1985). Early investigations of phylogenetic relationships in Glycine species based on their ability to produce fertile hybrids and the degree of meiotic chromosomes pairing (reviewed by Ratnaparkhe et al. 2011). This led to the designation of seven genome groups, i.e. reproductively compatible within species, while reproductively isolated from other genome groups (groups A to G, Singh and Hymowitz 1985;Hymowitz et al. 1998). Other data, including isozymes (Doyle et al. Open Access *Correspondence: bohsing@gate.sinica.edu.tw 1 Institute of Plant and Microbial Biology, Academia Sinica, Taipei 115, Taiwan Full list of author information is available at the end of the article 1986) and DNA phylogenetic analyses (Doyle et al. 1990a, b, c;Kollipara et al. 1997), suggested hybrid fertility as the primary criterion to identify species to genome groups and resulted in the designation of two additional genome groups, H and I. Besides, the situation in some species even more complicated, such as the G. tabacina and G. tomentella species complex, multiple separate 'races' are recognized that have not yet taxonomically arisen to species status, but that display enough molecular phylogenetic divergence from one another to be grouped into separate taxa (Doyle et al. 2004;Sherman-Broyles et al. 2014).
The Glycine species germplasm in Taiwan is abundant and important. Four species have been collected, including G. soja Siebold and Zucc. (GG genome, 2n = 40), G. tomentella Hayata (DDD 1 D 1 genome, 2n = 80), G. dolichocarpa Tateishi and Ohashi (A 6 A 6 DD genome, 2n = 80) and G. pescadrensis Hayata (A 6 A 6 B 3 B 3 genome, 2n = 80) (Hsing et al. 1995Tateishi and Ohashi 1992;Thseng et al. 1999;Tsai et al. 2001;Tsai 2006). G. soja is widely distributed in Russia, Korea, Japan, China and Taiwan (Singh and Hymowitz 1987), with Taiwan being the southernmost among these areas. The three wild perennial Glycine species, G. tomentella, G. dolichocarpa and G. pescadrensis, are found in Taiwan, the Ryukyu Islands, the Philippines, the South Pacific islands and Australia, with Taiwan being the northernmost of these areas (Hymowitz et al. 1998). All the remaining perennial Glycine species are found only in Australia (Tindale and Craven 1988). Therefore, the Ryukyu Islands, Taiwan, and its adjacent islands are unusual because they have representatives of both subgenera of Glycine. The collection and study of Glycine species in Taiwan and its adjacent islands thus provide information for better understanding the evolution and relationships between these two subgenera .
The soybean genome has been well documented as paleopolyploidy and manifested as allopolyploid diploidy (Gill et al. 2009;Mudge et al. 2005;Schlueter et al. 2007;Shoemaker et al. 1996;Tek et al. 2010). In addition, whole-genome sequencing (Williams 82) and resequencing projects on Glycine species were available recently Katayose et al. 2012;Kim et al. 2010;Lam et al. 2010;Schmutz et al. 2010;Weidner et al. 2012). However, little information is available on the genome composition of wild relatives; that is, cytological research of wild soybean is still limited.
Fluorescent in situ hybridization (FISH) is a powerful tool for cytological studies. FISH used with specific markers can be used to label specific chromosomes and target given genomes or particular regions on chromosomes (Cheng et al. 2007;Jiang and Gill 2006;Kopecky et al. 2008;Pinkel et al. 1986). Furthermore, nonspecific repeats such as simple repeat sequences (SSRs) may also be used as genomic and chromosomal markers among many species (Cuadrado et al. 2008;Cuadrado and Jouve 2007).
FISH has been used in cytological research on soybean. Several repeat sequences such as ribosomal RNA genes (rDNAs) (Shi et al. 1996), SB92 (Kolchinsky and Gresshoff 1995;Vahedian et al. 1995) and STR120 (Morgante et al. 1997) were labeled on soybean chromosomes and used in chromosome "painting" (Shi et al. 1996). FISH with bacterial artificial chromosome (BAC) clones used as probes suggested heterochromatic blocks on the chromosome pericentromeric region (Lin et al. 2005), and some research has provided evidence of chromosomelevel homeology in the paleopolyploid soybean genome (Pagel et al. 2004;Walling et al. 2006). In the study of soybean chromosome structure and genomic evolution, SB92 was previously suggested as a major repeat sequence in the G. max genome (Kolchinsky and Gresshoff 1995;Vahedian et al. 1995). Later, SB92 and its similar sequence SB91 were both defined as the centromere sequences of the GG genome by cytological features (Gill et al. 2009) and thus designated CentGm-1 (SB92) and CentGm-2 (SB91) (Tek et al. 2010). In addition, FISH revealed CentGm-1 or CentGm-2 present in about a half of all G. max chromosomes, which suggests allopolyploidy in the paleopolyploid genome (Gill et al. 2009). The G. max karyotyping system in mitotic metaphase was recently studied with a probe cocktail of BAC clones and partial centromere oligonucleotides (Findley et al. 2010), and this system has enabled the characterization of most of the chromosomal translocation lines (Findley et al. 2011).
The rDNAs have been widely used as reliable landmarks in chromosome karyotyping of several higher plants. In higher eukaryotes, the 45S and 5S rDNA loci are transcribed by different RNA polymerases and usually located in different positions on chromosomes (Srivastava and Schlessinger 1991). The number and location on chromosomes for both 45S and 5S may be diverse within related species (Chang et al. 2009;Chung et al. 2008;Pedrosa-Harand et al. 2006). In other words, karyotyping with rDNA loci may reflect the relationships of related species. Most of the recent cytological research on Glycine species has focused on the number of rDNA loci. The genome from diploid Glycine species, known as the paleopolyploidy, was predicted to consist of two 45S rDNA loci because it was ancient allopolyploidy; however, only one locus was observed on diploid Glycine species (Krishnan et al. 2001;Singh et al. 2001). Therefore, one 45S rDNA locus disappeared during soybean genome diploidization, whereby the ancient tetraploid became a new diploid paleopolyploidy.
Some studies revealed two 5S rDNA and two 45S rDNA loci in most allotetraploid genomes (Krishnan et al. 2001;Singh et al. 2001); however, an extra minor 45S rDNA locus was found in only two allotetraploid species (Krishnan et al. 2001). Therefore, the number of 45S loci among Glycine species or accessions was not as simple as the model of "diploid-doubled to tetraploid" of one locus for diploid and two loci for tetraploid species (Krishnan et al. 2001;Singh et al. 2001).
In the current study, we investigated the difference in repeat sequences between soybean and its wild related species. Here, we report three new soybean repeat sequences-SBRS1, SBRS2 and SBRS3-isolated from the soybean WSG data. In light of their distribution on chromosomes of GG genome, the three repeat sequences may be considered as FISH karyotyping markers for GG genome. We also labeled the three sequences in some wild Glycine species and found that they widely co-localized with 45S rDNA.

Plant materials
Glycine species used and locations where they were collected are in Table 1 along with their genome types and known USDA permanent plant introduction (PI) numbers or IL number designated by Dr. T. Hymowitz. All plant materials were previously described in detail . Cultivated soybean Shi-shi was kindly provided by Kaohsiung District Agricultural Research and Extension Station. Glycine species distributed in Taiwan and nearby islands were collected by Drs. J. S. Hsieh and Y. C. Huang, Department of Agronomy, National Taiwan University. Three diploid Glycine tomentella accessions originally from Papua New Guinea and Australia were kindly provided by Dr. T. Hymowitz, Department of Crop Science, University of Illinois, USA. These plants were grown in the greenhouse at Academia Sinica, Taiwan. Root tips of seedling were used for metaphase chromosome preparation.

Flow cytometry for estimating ploidy
Flow cytometric estimation of nuclear DNA content involved the MoFlo XDP Cell Sorter (Beckman Coulter, Fullerton, CA, USA), with propidium iodide (PI) as the fluorescent stain. A total of 20 mg of young leaves of each plant and standard nuclear fractions were prepared and mixed. Chicken erythrocyte nuclei (CEN; BioSure Inc., Grass Valley, CA, USA) and rainbow trout erythrocyte nuclei (TEN; NPE Systems Inc, Pembroke Pines, FL, USA) were used as internal standards and the DNA index (Vindelov et al. 1983). Leaf samples were chopped by using of a razor in a Petri dish containing 0.5 ml extraction buffer (Partec CyStain PI Absolute P Nuclei Extraction Buffer; Partec GMBH, Munster, Germany). The resulting extract was passed through two 30-μm filters before the addition of 0.5 ml Partec CyStain PI Absolute P Staining Buffer. Samples were kept in the dark for up to 30 min before flow cytometry. Six thousand cells were counted in each sample, with three replicates. The predicted genome size was calculated by assuming that 1 pg DNA = 965 Mbp (Bennett et al. 2000).

Predictions of repeat sequences in the soybean genome
Candidate repeat sequences were identified on the basis of the soybean WGS data . "eTandem" in EMBOSS 6.0 was used to scan all scaffold sequences, and the outputs were filtered by shell scripts and AWK script in Linux with repeat unit length >70 bp, homologous identity >70 % and copies >30. The longest Table 1 The accessions of the Glycine species used in this study, with their nuclear DNA content and collection locations a Accession number in our laboratory b 1 pg DNA = 965 Mbp (Bennett et al. 2000 repeats were used for primer design (Additional file 1: Table S1).

Syntheses of probes for FISH analysis
To prepare probes for FISH analysis, the pTA71 plasmid that contains a 9.1-kb fragment of 18S-5.8S-26S rDNA from common wheat (Gerlach and Bedbrook 1979) was labeled with biotin-16-dUTP by nick translation (Roche Diagnostics, Penzberg, Germany). Other repeat sequences were subcloned into pGEM-T Easy vector (Promega) by using Shi-shi genomic DNA and were labeled with digoxigenin-11-dUTP by using the Roche PCR DIG Probe synthesis kit (Roche).

Southern blot analysis
Genomic DNA (30 μg) of each Glycine accession was completely digested with restriction enzymes (EcoRI or NdeI, New England Biolabs) before separation by electrophoresis overnight on a 0.8 % agarose gel in 0.5X TBE buffer, then blotting onto Hybond-N+ nylon membranes (Amersham Pharmacia Biotech). Probe preparation, membrane hybridization, and signal detection followed the instructions of the ECL direct nucleic acid labeling and detection system (Amersham Pharmacia Biotech).

Chromosome preparation and FISH analysis
The terminal 0.5 cm of young roots was collected from individual seedlings, pretreated with 2 mM 8-hydroxyquinoline for 2 h at room temperature and fixed overnight in Farmer's fixative [95 % ethanol + glacial acetic acid (3:1 v/v)]. The fixed root tips were washed with distilled water before incubation in a cell wall digestive enzyme mixture of 6 % pectinase (Sigma Chemical, St. Louis, MO, USA) and 6 % cellulase (Onozuka R-10; Yakult Honsha, Tokyo, Japan) in 75 mM KCl (pH = 4.0) buffer at 37 °C for 1.5-3 h depending on the characteristics of each accession. After a wash with distilled water, the terminal 1 mm of softened tissues was squashed on a slide (Chang et al. 2009). Samples were pretreated with pepsin (1.25 mg/ ml for G. soja, 5 mg/ml for other accessions, in 10 mM HCl) at 37 °C for 3 to 8 h depending on the characteristics of each accession. FISH analysis was performed as described previously (Chung et al. 2008).

Repeat sequences from WGS data and the analysis
To investigate the repeat sequence component in the Glycine genome, we used the WGS data to find new repeat sequences. We identified 17 candidate repeat sequences by bioinformatics analysis. Several of these candidates belong to the LTR120 family (Morgante et al. 1997) and SB92 family (Kolchinsky and Gresshoff 1995). Finally, we chose 3 new repeat sequences, designated soybean repeat sequence1 (SBRS1), SBRS2 and SBRS3; their sequences and designed primers are in Additional files 1 and 2.

SBRS2
A 124 bp minisatellite sequence, represents approximately 6000 copies in the WGS data. About 95 % of the SBRS2 repeats were aligned into pseudomolecules. The SBRS2 is a head-to-tail tandem repeat with usually less than 10 copies in each continguous region and is dispersed on all of the 20 pseudomolecules in the WSG data. The most significant region rich in SBRS2 was located at Gm9, where 1095 copies of the SBRS2 clustered in a 270-kb region.

SBRS3
A 201 bp minisatellite sequence, represents 8206 copies in the WGS data; 8016 copies were aligned into pseudomolecules. It is also a tandem repeat, with usually less than 10 copies in each region, dispersed on all 20 pseudomolecules. Three regions especially contain a high density of SBRS3 repeats, including 2 on Gm14 (1512 copies were separated into 2 groups 26,000 kb apart) and the other one on Gm20 (400 copies clustered together).
To understand the existence and the distribution of these three repeat sequences in the Glycine genome, we used them as probes in Southern blot analyses with genomic DNA of the cultivated soybean and the collected Glycine species. All three repeats, SBRS1, SBRS2 and SBRS3, showed heavy ladder signals in the GG genome species, G. max and G. soja (Lanes 1-4 of Fig. 1; Additional file 3: Figure S1, Additional file 4: Figure S2). These heavy ladder signals are similar between G. max and G. soja, but the patterns among the three repeat sequences are different. However, the 3 repeat sequences showed single band and the similar size in other genome species, including G. tomentella, G. dolichocarpa and G. pescadrensis, a 3.8-and a 9-kbp band were found in EcoRI-and NdeI-digested DNA samples, respectively, of all perennial Glycine species tested (Lanes 5-12 of Fig. 1. Additional file 3: Figure S1, Additional file 4: Figure S2). Therefore, the three repeat sequences exist in each Glycine species tested, with high copies in GG genome, but with very low copies in other species.

FISH analysis with SBRS1, SBRS2, SBRS3 and 45S rDNA
FISH analysis was used in further evaluated the distribution of the three repeat sequences in the genome of Glycine species. We used the root tips of G. soja to represent the GG genome on FISH analysis because the cultivated soybean and G. soja both represent the GG genome, and band patterns on Southern blot analysis were similar (Fig. 1, Additional file 3: Figure S1, Additional file 4: Figure S2). The 45S rDNA probe was used as a control, and it was located at the distal ends of one pair of chromosomes of G. soja (Fig. 2a, green). SBRS1 signals were dispersed on every chromosome of G. soja as shown by metaphase (Panel A) and interphase (Panel B) chromosomes (Fig. 2a, red). Throughout the genome, the distribution of SBRS1 signals varied among chromosomes. SBRS2 signals were detected as clusters on two pairs of chromosomes of G. soja (Fig. 3). FISH signals for one of the SBRS2 clusters overlapped with 45S rDNA at the distal ends of one pair of chromosomes (Fig. 3a, yellowish green), and the second SBRS2 cluster was detected in the proximal regions of another pair of chromosomes (Fig. 3a, b, red). SBRS3 signals were detected on three pairs of chromosomes of G. soja (Fig. 4,  red), and none overlappped with 45S rDNA (Fig. 4,  green). One of the SBRS3 signals was detected at the subdistal ends of one pair of chromosomes, and the other 2 were detected in the proximal regions and the distal ends of another pair of chromosomes. We labeled probes of the three repeat sequences and 45S rDNA in allotetraploid Glycine species (i.e., Tom034, Tom039 and Tab 004; Fig. 5), and their diploid relatives (i.e., Tom052, Tom062 and Tom051; Additional file 5: Figure S3; Additional file 6: Figure S4, Additional file 7: Figure S5, Additional file 8: Figure S6). Most of the three tandem repeat sequences signals co-localized with 45S rDNA signals among the tetraploid and diploid relatives tested (Figs. 3, 5; Additional file 5: Figure S3, Additional file 6: Figure S4, Additional file 8: Figure S6). SBRS1 and SBRS3 signals did not co-localize with 45S rDNA signals in the 2 diploid genomes, GG genome (Figs. 2, 4) or DD genome (Additional file 7: Figure S5, Additional file 8: Figure S6). In the genome of Tom062 (DD genome), very weak SBRS1 signals located on many chromosomes (Additional file 7: Figure S5), SBRS2 signals co-localized with 45S rDNA (Additional file 8: Figure S6), and no substantial SBRS3 signal was detected. In addition, all tetraploid accessions showed two 45S rDNA loci, with a minor signal detected in the proximal regions of another pair of chromosomes of Tab 074 (Fig. 6). FISH analysis results are summarized in Table 2.

FISH analysis with SSR probe in Glycine species
We additionally used 5 types of SSR probes (i.e., ATT, AT, CAA, CT and CTT) in FISH analysis. However, only one SSR probe, with 33 copies of the repeat ATT, on G. soja, gave interesting results, and, similar to SBRS1, ATT signals were detected as dispersed signals on every chromosomes of the GG genome (Fig. 7, red). There was no ATT signal on the mitotic chromosomes of most other Glycine species except Tom034. Few very weak ATT probe signals were detected on this accession, but the signals were unstable. In addition, the signal distribution for these SSRs varied among chromosomes.

Discussion
The predicted repeat sequences were confirmed by cytological study The soybean WGS sequence project has set significantly new milestones for understanding the soybean genome and its evolutionary history Schmutz et al. 2010). In the current study, we identified three new repeat sequences from soybean WGS data. Using bioinformatics, we propose their characterization in the soybean genome. We used WGS as a resource in chromosome marker design for the cultivated soybean and other Glycine species, and found that use of the WGS data was powerful for searching for candidates of repeat sequences and designing probes for cytological research.
For the GG genome, the WGS data could be used to predict general copy numbers of the repeat sequences in each cluster and the location of repeat-sequence condensed regions. FISH analyses confirmed these bioinformatics findings (Figs. 2, 3, 4; Table 2). WGS data may be also used in cytological research for the wild related species. Combining bioinformatics and cytological strategies may help bridge the gap between cultivated species and related species. In the current study, we analysis the three new repeat sequences on many Glycine species by FISH. Similar studies had shown that FISH analysis would reveal heterochromatin  (Ananiev et al. 1998 for maize;Badaeva et al. 1996 for Aegilops species; Kulikova et al. 2004 for Medicago species). Because Glycine species consist of high chromosome numbers (38, 40, 78, or 80;Hymowitz et al. 1998), as well, these chromosomes are compact and small, the utilization of repetitive sequences as probe for FISH analysis may help on the cytogenesis and breeding analysis in the post-genomics era.

Possible evolutionary process for the three new repeat sequences among Glycine species
Differences in the distribution of repetitive DNA sequences were shown between Glycine species in the current study (summarized in Table 2). The chromosomal localization of these DNA probes was different.
The possible association of the these repeat sequences and rDNA loci may be due to (1) the three repeat sequences having existed in the rDNA of the genome of their common earlier ancestor, (2) all having expanded and dispersed in GG genome, and finally (3) SBRS1 and SBRS3 being reduced or disappeared in 45S rDNA loci during the GG (Soja001) and DD (Tom062) genome evolution and SBRS2 still existing in the 45S rDNA of all tested Glycine species (Table 2). Dolye and his colleague used global, gene-specific, and relaxed clock methods to date the polyploid history of soybean and suggested the date for the G. max/G. tomentella split averaged between 3.8 and 6.86 Mya across methods (Egan and Doyle 2010). It coincides well with the fact that the shared repeat sequences, but with different amount of such fragments, among the Glycine species. Many satellite sequences are homologs to the IGS subrepeats of rDNA from other plants genome, such as legumes (Falquet et al. 1997), potato (Stupar et al. 2002), tobacco (Lim et al. 2004;Volkov et al. 1999) andtomato (Jo et al. 2009). In tomato study, the insertion of a retrotransposon into the rDNA gene enhanced the amplification of the nearby IGS subrepeats. Those amplified subrepeats then dispersed into other loci and became repeat sequences (Jo et al. 2009). Other studies also proposed that some transposon elements, such as En/Spmlike transposons, could be involved in rDNA movement and as one of the key steps when the rDNA subrepeat turned into repeat sequences (Raskina et al. 2004a(Raskina et al. , 2008Schubert and Wobus 1985).

Diversity of the genome in Glycine species
We analyzed the DNA content of the nine Glycine accessions, five of them are diploid and another four are allotetraploid, by flow cytometry. The two G. pescadrensis accessions Tab 004 and Tab 074 collected from Penghu and Kingmen, differed in DNA content (Table 1). In addition, the DNA content of Soja001 and the cultivated soybean (GG genome) differs up to 22 % (Table 1). The DNA content of G. dolichocarpa (Tom039) and G. tomentella (Tom034) are more than the sum of their two parental diploid species (Table 1).
Research using Helianthus (Ungerer et al. 2006) suggested that retrotransposons may play an important role in the genome expansion process of allotetraploid evolution. The dynamics of transposon element driven genome evolution have been reported for several plant genomes (El Baidouri and Panaud 2013). In fact, retrotransposons occupy the largest part of the cultivated soybean genome ) and may play a role in Glycine genome size variations.
The number of 45S rDNA loci among wild soybeans is dynamic. In previous studies, the number of 45S rDNA loci among the tetraploid Glycine species was suggested to be the same (Krishnan et al. 2001;Singh et al. 2001). However, two major and one minor loci of 45S rDNA signals were detected in Tab 074 (G. pescadrensis, A 6 A 6 B 3 B 3 genome). The minor signals are found on the middle area of a pair of the chromosomes. Only two major loci of 45S rDNA were found for the other tetraploids, including Tab 004. These two G. pescadrensis accessions were collected from Penghu and Kingmen, and the distance between these two islands is about 140 km, the geographical isolation may explain the diversity of 45S rDNA loci. Thus,  the evolution of the 45S rDNA in the Glycine genome is more complex than the model of "diploid-doubled to tetraploid". Variations in the 45S rDNA loci among natural hybrids and related species were reported for Lycoris (Chang et al. 2009), Oryza (Chung et al. 2008), Triticeae (Kim et al. 1993), and Arabidopsis (Maluszynska and Heslopharrison 1993). As well, the ability to alter the number and location of rDNA loci among closely related species or even accessions was noted (Chung et al. 2008;Raskina et al. 2004b). In our study, the variation in number of rDNA loci in Glycine was more dramatic than in previous studies (Krishnan et al. 2001;Singh et al. 2001). Such differences in number or location of 45S rDNA may be used as chromosome markers.

Genome and chromosome markers of Glycine species
Specific repeat sequences, such as SBRS1-3, and nonspecific repeat sequences, such as SSR and rDNA, may be potential genome or chromosome markers. The nonspecific repeat sequences can also be used as specific markers, if the distributions of repeats are different enough among genomes. For instance, the SSRs, which are widely identified among species (Varshney et al. 2005), have been used to label the chromosomes of different genomes in hexaploidy Triticum aestivum (Cuadrado et al. 2000(Cuadrado et al. , 2008. We use nonspecific repeat sequences as probes for labeling chromosomes of several Glycine species. Signals of ATT were exclusively detected on the GG genome, with no difference between G. soja and G. max (data not shown). Although 2-3, weak ATT signals could be detected on the chromosome of Tom034 (data not shown), they were easy to be separated from the strong ones of GG genome.
The study of the Glycine genome is important for soybean improvement; combing bioinformatics and cytology techniques can advance investigations. Our study is just a beginning to show that the soybean WSG data and cytological methods can be used to find new repeat sequences, reveal the genome evolution process, and provide addition information in investigation of the phylogeny in the Glycine species.