Soybean bacterial artificial chromosome contigs anchored with RFLPs: insights into genome duplication and gene clustering

Dynamic genetic features of chromosomes revealed by comparison of soybean genetic and sequence-based physical maps

Despite the intensive soybean [Glycine max (L.) Merrill] genome studies, the high chromosome number (20) of the soybean plant relative to many other major crops has hindered the development of a high-resolution genomewide genetic map derived from a single population. Here, we report such a map, which was constructed in an F15 population derived from a cross between G. max and G. soja lines using indel polymorphisms detected via a G. soja genome resequencing. By targeting novel indel markers to marker-poor regions, all marker intervals were reduced to under 6 cM on a genome scale. Comparison of the Williams 82 soybean reference genome sequence and our genetic map indicated that marker orders of 26 regions were discrepant with each other. In addition, our comparison showed seven misplaced and two absent markers in the current Williams 82 assembly and six markers placed on the scaffolds that were not incorporated into the pseudomolecules. Then, we showed that, by determining the missing sequences located at the presumed beginning points of the five major discordant segments, these observed discordant regions are mostly errors in the Williams 82 assembly. Distributions of the recombination rates along the chromosomes were similar to those of other organisms. Genotyping of indel markers and genome resequencing of the two parental lines suggested that some marker-poor chromosomal regions may represent introgression regions, which appear to be prevalent in soybean. Given the even and dense distribution of markers, our genetic map can serve as a bridge between genomics research and breeding programs.

Using microsatellites to understand the physical distribution of recombination on soybean chromosomes

Soybean is a major crop that is an important source of oil and proteins. A number of genetic linkage maps have been developed in soybean. Specifically, hundreds of simple sequence repeat (SSR) markers have been developed and mapped. Recent sequencing of the soybean genome resulted in the generation of vast amounts of genetic information. The objectives of this investigation were to use SSR markers in developing a connection between genetic and physical maps and to determine the physical distribution of recombination on soybean chromosomes. A total of 2,188 SSRs were used for sequence-based physical localization on soybean chromosomes. Linkage information was used from different maps to create an integrated genetic map. Comparison of the integrated genetic linkage maps and sequence based physical maps revealed that the distal 25% of each chromosome was the most marker-dense, containing an average of 47.4% of the SSR markers and 50.2% of the genes. The proximal 25% of each chromosome contained only 7.4% of the markers and 6.7% of the genes. At the whole genome level, the marker density and gene density showed a high correlation (R(2)) of 0.64 and 0.83, respectively with the physical distance from the centromere. Recombination followed a similar pattern with comparisons indicating that recombination is high in telomeric regions, though the correlation between crossover frequency and distance from the centromeres is low (R(2) = 0.21). Most of the centromeric regions were low in recombination. The crossover frequency for the entire soybean genome was 7.2%, with extremes much higher and lower than average. The number of recombination hotspots varied from 1 to 12 per chromosome. A high correlation of 0.83 between the distribution of SSR markers and genes suggested close association of SSRs with genes. The knowledge of distribution of recombination on chromosomes may be applied in characterizing and targeting genes.

A novel manual pooling system for preparing three-dimensional pools of a deep coverage soybean bacterial artificial chromosome library

Compared with hybridization-based techniques, polymerase chain reaction-based screening of large insert libraries has been used widely as it is fast, easy and sensitive. However, various pooling strategies are needed to ensure efficient screening. It is time-consuming and labourious to prepare three-dimensional pools for a deep coverage bacterial artificial chromosome (BAC) library of soybean (1.12 × 10(9)  bp) in the absence of robotic facility. In the present study, we describe a novel manual pooling system for preparing three-dimensional pools of a soybean BAC library. This simple technique enables a single researcher to construct three-dimensional pools for a deep-coverage (12 haploid genome equivalents) BAC library of soybean in less than 2 months without any robotic manipulation. When the prepared three-dimensional pools were screened with 29 polymerase chain reaction-based markers, an average of 9.2 clones per marker were identified. These identified clones will be useful either in quantitative trait loci gene isolation or in synteny study between soybean and other legumes including Lotus japonicus. This efficient pooling system could be applied to any other BAC libraries without the need for robotic manipulation.

Genome structure in soybean revealed by a genomewide genetic map constructed from a single population

A complete genetic linkage map of the soybean, in which sequence-based (SB) genetic markers are evenly distributed genomewide, was constructed from an F(12) population composed of 113 recombinant inbred lines derived from an interspecific cross involving Korean genotypes Hwangkeum and IT182932. Several approaches were employed for the development of 112 novel SB markers targeting both the gaps and the ends of the linkage groups (LGs). The resultant map harbored 20 well-resolved LGs presumed to correspond to the 20 pairs of soybean chromosomes. The map allowed us to identify the important chromosomal structures that were not observed in the integrated genetic maps, to identify the new potentially gene-rich regions, to detect segregation distortion regions within the whole genome, and to extend the ends of the LGs. The results will facilitate the further discovery of agronomically relevant genetic loci in the heretofore neglected chromosomal regions and should also provide some important links between the soybean genetic, physical, and genome sequence maps in the regions.

Genetic marker anchoring by six-dimensional pools for development of a soybean physical map

BACKGROUND

Integrated genetic and physical maps are extremely valuable for genomic studies and as important references for assembling whole genome shotgun sequences. Screening of a BAC library using molecular markers is an indispensable procedure for integration of both physical and genetic maps of a genome. Molecular markers provide anchor points for integration of genetic and physical maps and also validate BAC contigs assembled based solely on BAC fingerprints. We employed a six-dimensional BAC pooling strategy and an in silico approach to anchor molecular markers onto the soybean physical map.

RESULTS

A total of 1,470 markers (580 SSRs and 890 STSs) were anchored by PCR on a subset of a Williams 82 BstY I BAC library pooled into 208 pools in six dimensions. This resulted in 7,463 clones (approximately 1x genome equivalent) associated with 1470 markers, of which the majority of clones (6,157, 82.5%) were anchored by one marker and 1106 (17.5%) individual clones contained two or more markers. This contributed to 1184 contigs having anchor points through this 6-D pool screening effort. In parallel, the 21,700 soybean Unigene set from NCBI was used to perform in silico mapping on 80,700 Williams 82 BAC end sequences (BES). This in silico analysis yielded 9,835 positive results anchored by 4152 unigenes that contributed to 1305 contigs and 1624 singletons. Among the 1305 contigs, 305 have not been previously anchored by PCR. Therefore, 1489 (78.8%) of 1893 contigs are anchored with molecular markers. These results are being integrated with BAC fingerprints to assemble the BAC contigs. Ultimately, these efforts will lead to an integrated physical and genetic map resource.

CONCLUSION

We demonstrated that the six-dimensional soybean BAC pools can be efficiently used to anchor markers to soybean BACs despite the complexity of the soybean genome. In addition to anchoring markers, the 6-D pooling method was also effective for targeting BAC clones for investigating gene families and duplicated regions in the genome, as well as for extending physical map contigs.

A soybean transcript map: gene distribution, haplotype and single-nucleotide polymorphism analysis

The first genetic transcript map of the soybean genome was created by mapping one SNP in each of 1141 genes in one or more of three recombinant inbred line mapping populations, thus providing a picture of the distribution of genic sequences across the mapped portion of the genome. Single-nucleotide polymorphisms (SNPs) were discovered via the resequencing of sequence-tagged sites (STSs) developed from expressed sequence tag (EST) sequence. From an initial set of 9459 polymerase chain reaction primer sets designed to a diverse set of genes, 4240 STSs were amplified and sequenced in each of six diverse soybean genotypes. In the resulting 2.44 Mbp of aligned sequence, a total of 5551 SNPs were discovered, including 4712 single-base changes and 839 indels for an average nucleotide diversity of Theta= 0.000997. The analysis of the observed genetic distances between adjacent genes vs. the theoretical distribution based upon the assumption of a random distribution of genes across the 20 soybean linkage groups clearly indicated that genes were clustered. Of the 1141 genes, 291 mapped to 72 of the 112 gaps of 5-10 cM in the preexisting simple sequence repeat (SSR)-based map, while 111 genes mapped in 19 of the 26 gaps >10 cM. The addition of 1141 sequence-based genic markers to the soybean genome map will provide an important resource to soybean geneticists for quantitative trait locus discovery and map-based cloning, as well as to soybean breeders who increasingly depend upon marker-assisted selection in cultivar improvement.

Sequence conservation of homeologous bacterial artificial chromosomes and transcription of homeologous genes in soybean (Glycine max L. Merr.)

The paleopolyploid soybean genome was investigated by sequencing homeologous BAC clones anchored by duplicate N-hydroxycinnamoyl/benzoyltransferase (HCBT) genes. The homeologous BACs were genetically mapped to linkage groups C1 and C2. Annotation of the 173,747- and 98,760-bp BACs showed that gene conservation in both order and orientation is high between homeologous regions with only a single gene insertion/deletion and local tandem duplications differing between the regions. The nucleotide sequence conservation extends into intergenic regions as well, probably due to conserved regulatory sequences. Most of the homeologs appear to have a role in either transcription/DNA binding or cellular signaling, suggesting a potential preference for retention of duplicate genes with these functions. Reverse transcriptase-PCR analysis of homeologs showed that in the tissues sampled, most homeologs have not diverged greatly in their transcription profiles. However, four cases of changes in transcription were identified, primarily in the HCBT gene cluster. Because a mapped locus corresponds to a soybean cyst nematode (SCN) QTL, the potential role of HCBT genes in response to SCN is discussed. These results are the first sequenced-based analysis of homeologous BACs in soybean, a diploidized paleopolyploid.

Survey sequencing of soybean elucidates the genome structure, composition and identifies novel repeats

In order to expand our knowledge of the soybean genome and to create a useful DNA repeat sequence database, over 24 000 DNA fragments from a soybean [Glycine max (L.) Merr.] cv. Williams 82 genomic shotgun library were sequenced. Additional sequences came from over 29 000 bacterial artificial chromosome (BAC) end sequences derived from a BstI library of the cv. Williams 82 genome. Analysis of these sequences identified 348 different DNA repeats, many of which appear to be novel. To extend the utility of the work, a pilot study was also conducted using methylation filtration to estimate the hypomethylated, soybean gene space. A comparison between 8366 sequences obtained from a filtered library and 23 788 from an unfiltered library indicate a gene-enrichment of ~3.2-fold in the hypomethylated sequences. Given the 1.1-Gb soybean genome, our analysis predicts a ~343-Mb hypomethylated, gene-rich space.

Paleopolyploidy and gene duplication in soybean and other legumes

Two of the most important observations from whole-genome sequences have been the high rate of gene birth and death and the prevalence of large-scale duplication events, including polyploidy. There is also a growing appreciation that polyploidy is more than the sum of the gene duplications it creates, in part because polyploidy duplicates the members of entire regulatory networks. Thus, it may be important to distinguish paralogs that are produced by individual gene duplications from the homoeologous sequences produced by (allo)polyploidy. This is not a simple task, for several reasons, including the chromosomally cryptic nature of many duplications and the variable rates of gene evolution. Recent progress has been made in understanding patterns of gene and genome duplication in the legume family, specifically in soybean.

Chromosome-level homeology in paleopolyploid soybean (Glycine max) revealed through integration of genetic and chromosome maps

Soybean has 20 chromosome pairs that are derived from at least two rounds of genomewide duplication or polyploidy events although, cytogenetically, soybean behaves like a diploid and has disomic inheritance for most loci. Genetically anchored genomic clones were used as probes for fluorescence in situ hybridization (FISH) to determine the level of postpolyploid chromosomal rearrangements and to integrate the genetic and physical maps to (1) assign linkage groups to specific chromosomes, (2) assess chromosomal structure, and (3) determine the distribution of recombination along the length of a chromosome. FISH mapping of seven putatively gene-rich BACs from linkage group L (chromosome 19) revealed that most of the genetic map correlates to the highly euchromatic long arm and that there is extensive homeology with another chromosome pair, although colinearity of some loci does appear to be disrupted. Moreover, mapping of BACs containing high-copy sequences revealed sequestration of high-copy repeats to the pericentromeric regions of this chromosome. Taken together, these data present a model of chromosome structure in a highly duplicated but diploidized eukaryote, soybean.

Analysis and mapping of randomly chosen bacterial artificial chromosome clones from hexaploid bread wheat

The current view of wheat genome composition is that genes are compartmentalized into gene-rich and gene-poor regions. This model can be tested by analyzing randomly selected bacterial artificial chromosome (BAC) clones for gene content, followed by placement of these BACs onto physical and genetic maps. Map localization could be difficult for BACs that consist entirely of repeated elements. We therefore developed a technique where repeat junctions are used to generate unique markers. Four BAC clones from hexaploid wheat variety Chinese Spring were randomly selected and sequenced at 4- to 6-fold redundancy. About 50% of the BAC sequences corresponded to previously identified repeats, mainly LTR-retrotransposons, whereas most of the remaining DNA consisted of sequences with unknown origin or function. The average gene content was <1%, although each BAC contained one or two identified genes. Repeat boundaries were amplified and used to map each clone to a chromosome arm. Extrapolation from wheat-rice comparative knowledge suggests that three of the four BAC clones originate from "gene-rich" regions of the wheat genome. Nevertheless, because these BACs carry only a single gene (two BACs) or two genes (one BAC), the predicted gene density is approximately 1 gene per 75 kb, which is considerably lower than previously estimated gene densities (one gene per 5-20 kb) for gene-rich regions in wheat. This analysis of randomly selected wheat BAC clones suggests that genes are more evenly distributed in wheat than previously believed and substantiates the need for large-scale random BAC sequencing to determine wheat genome organization.