Organization,expression and evolution of a disease resistance gene cluster in soybean

PCR amplification was previously used to identify a cluster of resistance gene analogues (RGAs) on soybean linkage group J. Resistance to powdery mildew (Rmd-c), Phytophthora stem and root rot (Rps2), and an ineffective nodulation gene (Rj2) map within this cluster. BAC fingerprinting and RGA-specific primers were used to develop a contig of BAC clones spanning this region in cultivar "Williams 82" [rps2, Rmd (adult onset), rj2]. Two cDNAs with homology to the TIR/NBD/LRR family of R-genes have also been mapped to opposite ends of a BAC in the contig Gm_Isb001_091F11 (BAC 91F11). Sequence analyses of BAC 91F11 identified 16 different resistance-like gene (RLG) sequences with homology to the TIR/NBD/LRR family of disease resistance genes. Four of these RLGs represent two potentially novel classes of disease resistance genes: TIR/NBD domains fused inframe to a putative defense-related protein (NtPRp27-like) and TIR domains fused inframe to soybean calmodulin Ca(2+)-binding domains. RT-PCR analyses using gene-specific primers allowed us to monitor the expression of individual genes in different tissues and developmental stages. Three genes appeared to be constitutively expressed, while three were differentially expressed. Analyses of the R-genes within this BAC suggest that R-gene evolution in soybean is a complex and dynamic process.     

BAC contig development by fingerprint analysis in soybean.

We constructed a soybean bacterial artificial chromosome (BAC) library suitable for map-based cloning and physical mapping in soybean. This library consists of approximately 40 000 clones (4-5 genome equivalents) stored individually in 384-well microtiter dishes. A random sampling of 224 clones yielded an average insert size of 150 kb, giving a 98% probability of recovering any specific sequence. We screened the library for seven single or very low copy genie or genomic sequences using the polymerase chain reaction (PCR) and found between one and seven BACs for each of the seven sequences. When testing the library with a portion of the soybean psbA chloroplast gene, we found less than 1% chloroplast DNA representation. We also screened the library for eight different classes of disease resistance gene analogs (RGAs) and identified BACs containing all RGAs except class 8. We arranged nine of the class 1 RGA BACs and six of the class 3 RGA BACs into individual contigs based on fingerprint patterns observed after Southern probing of restriction digests of the member BACs with a class-specific sequence. This resulted in the partial localization of the different multigene family sequences without precise definition of their exact positions. Using PCR-based end rescue techniques and RFLP mapping of BAC ends, we mapped individual BACs of each contig onto linkage group J of the soybean public map. The class 1 contig mapped to the region on linkage group J that contains several disease resistance genes. The class 1 contig extended approximately 400 kb. The arrangement of the BACs within this contig has been confirmed using PCR. One end of the class 1 contig core BAC mapped to two positions on linkage group J and cosegregated with two class 1 RGA loci, suggesting that this segment is within an area of regional duplication.     

A physical map of the 20q12-q13.1 region associated with type 2 diabetes

Several recent genetic studies have suggested linkage of Type 2 diabetes (non-insulin-dependent diabetes mellitus) susceptibility to a region of chromosome 20q12-q13.1. To facilitate the identification and cloning of a diabetes susceptibility gene(s) in this region, we have constructed correlated radiation hybrid and YAC/BAC contig physical maps of the region. A high-resolution radiation hybrid map encompassing 9.5 Mb between the PLC and the CEBPB genes was constructed using 68 markers: 25 polymorphic markers, 15 known genes, 21 ESTs, and 7 random genomic sequences. The physical order of the polymorphic markers within this radiation hybrid map is consistent with published genetic maps. A YAC/BAC contig that gives continuous coverage between PLC and CEBPB was also constructed. This contig was constructed from 24 YACs, 34 BACs, and 1 P1 phage clone onto which 71 markers were mapped: 23 polymorphic markers, 12 genes, 24 ESTs, and 12 random genomic sequences. The radiation hybrid map and YAC/BAC physical map enable precise mapping of newly identified transcribed sequences and polymorphic markers that will aid in linkage and linkage disequilibrium studies and facilitate identification and cloning of candidate Type 2 diabetes susceptibility genes residing in 20q12-q13.1.     

Resolving the aphid resistance locus Sd-1 on a BAC contig within a sub-telomeric region of Malus linkage group 7.

Aphids cause serious physical and economic damage to most major crops throughout the world, and there is a pressing requirement to isolate genes conferring aphid resistance. The Sd-1 locus in Malus spp. (apple) confers resistance against the rosy leaf-curling aphid (Dysaphis devecta Wlk.), and was recently positioned within a 1.3-cM region on linkage group 7, flanked by molecular markers. These markers were used as a basis for development of a BAC contig spanning the locus, together with adapter-mediated amplification of flanking sequences to obtain BAC insert-end sequences, and fingerprinting of BAC clones. Approximately 800 kb of the Sd-1 genomic region was covered by 19 overlapping BACs, with an average insert size of 75-150 kb. The physical-genetic distance ratio was estimated at 460 kb/cM, although the distribution of recombination events was irregular with respect to estimated physical distance. Recombinant analysis and development of new markers allowed Sd-1 to be positioned within an interval of approximately 180 kb located on either of two overlapping BACs. From one of these, an insert end sequence showed a significant degree of similarity to nucleotide binding site-leucine rich repeat (NBS-LRR) resistance genes. Fluorescent in situ hybridization (FISH) of BAC clones within the contig enabled positioning and orientation of the locus within a euchromatic region, very close to the telomere of linkage group 7.     

Expression and genome organization of resistance gene analogs in soybean.

Sequence analysis of cloned plant disease-resistance genes reveals a number of conserved domains. Researchers have used these domains to amplify analogous sequences, resistance gene analogs (RGAs), from soybean and other crops. Many of these RGAs map in close proximity to known resistance genes. While this technique is useful in identifying potential disease resistance loci, identifying the functional resistance gene from a cluster of homologs requires sequence information from outside of these conserved domains. To study RGA expression and to determine the extent of their similarity to other plant resistance genes, two soybean cDNA libraries (root and epicotyl) were screened by hybridization with RGA class-specific probes. cDNAs hybridizing to RGA probes were detected in each library. Two types of cDNAs were identified. One type was full-length and contained several disease-resistance gene (R-gene) signatures. The other type contained several deletions within these signatures. Sequence analyses of the cDNA clones placed them in the Toll-Interleukin-1 receptor, nucleotide binding domain, and leucine-rich repeat family of disease-resistance genes. Using clone-specific primers from within the 3' end of the LRRs, we were able to map two cDNA clones (LM6 and MG13) to a BAC contig that is known to span a cluster of disease-resistance genes.     

High-resolution human/goat comparative map of the goat polled/intersex syndrome (PIS): the human homologue is contained in a human YAC from HSA3q23

The genetic and cytogenetic map around the chromosome 1 region shown to be linked with polledness and intersexuality (PIS) in the domestic goat (Capra hircus) was refined. For this purpose, a goat BAC library was systematically screened with primers from human coding sequences, scraped chromosome 1 DNA, bovine microsatellites from the region, and BAC ends. All the BACs (n = 30) were mapped by fluorescence in situ hybridization (FISH) on goat chromosome 1q41-q45. The genetic mapping of 30 new goat polymorphic markers, isolated from these BACs, made it possible to reduce the PIS interval to a region of less than 1 cM on goat chromosome 1q43. The PIS locus is now located between the two genes ATP1B and COP, which both map to 3q23 in humans. Genetic, cytogenetic, and comparative data suggest that the PIS region is now probably circumscribed to an approximately 1-Mb DNA segment for which construction of a BAC contig is in progress. In addition, a human YAC contig encompassing the blepharophimosis-ptosis-epicanthus-inversus region was mapped by FISH to goat chromosome 1q43. This human disease, mapped to HSA 3q23 and affecting the development and maintenance of ovarian function, could be a potential candidate for goat PIS.     

Physical and cDNA mapping in the DBH region of human chromosome 9q34

Chromosome 9q34 has been extensively studied and mapped due to the presence of known disease genes, principally tuberous sclerosis 1 (TSC1), in this region. During the course of our mapping of this region we constructed a 555-kb contig beginning approximately 50 kb proximal to the dopamine-beta-hydroxylase (DBH) gene and extending, with one small deletion, distal to the D9S114 marker. The contig consists of 11 P1 clones, four PAC clones, one BAC clone and six cosmid clones and contains 27 new nonpolymorphic STSs. We have found the region to be unstable in P1, PAC and BAC cloning vehicles and have identified several deleted genomic clones. In addition, we have isolated and mapped the 3' portions of three putative genes located within or immediately distal to the DBH gene, including one large gene that runs on the opposite strand to DBH and utilizes portions of two DBH exons. The genomic clones of the contig, cDNAs and new STSs will be useful reagents for the further study and mapping of this region.     

Generation of a soybean BAC library, and identification of DNA sequences tightly linked to the Rps1-k disease resistance gene

 A soybean bacterial artificial chromosome (BAC) library, comprising approximately 45 000 clones, was constructed from high-molecular-weight nuclear DNA of cultivar Williams 82, which carries the Rps1-k gene for resistance against Phytophthora sojae. The library is stored in 130 pools with about 350 clones per pool. Completeness of the library was evaluated for 21 random sequences including four markers linked to the Rps1 locus and 16 cDNAs. We identified pools containing BACs for all sequences except for one cDNA. Additionally, when screened for possible contaminating BAC clones carrying chloroplast genes, no sequences homologous to two barley chloroplast genes were found. The estimated average insert size of the BAC clones was about 105 kb. The library comprises about four genome equivalents of soybean DNA. Therefore, this gives a probability of 0.98 of finding a specific sequence from this library. This library should be a useful resource for the positional cloning of Rps1-k, and other soybean genes. We have also evaluated the feasibility of an RFLP-based screening procedure for the isolation of BAC clones specific for markers that are members of repetitive sequence families, and are linked to the Rps1-k gene. We show that BAC clones isolated for two genetically linked marker loci, Tgmr and TC1-2, are physically linked. Application of this method in expediting the map-based cloning of a gene, especially from an organism, such as soybean, maize and wheat, with a complex genome is discussed. Received: 12 May 1998/Accepted: 24 August 1998     

Integration of animal linkage and BAC contig maps using overgo hybridization

The alignment of genome linkage maps, defined primarily by segregation of sequence-tagged site (STS) markers, with BAC contig physical maps and full genome sequences requires high throughput mechanisms to identify BAC clones that contain specific STS. A powerful technique for this purpose is multi-dimensional hybridization of "overgo" probes. The probes are chosen from available STS sequence data by selecting unique probe sequences that have a common melting temperature. We have hybridized sets of 216 overgo probes in subset pools of 36 overgos at a time to filter-spotted chicken BAC clone arrays. A four-dimensional pooling strategy, including one degree of redundancy, has been employed. This requires 24 hybridizations to completely assign BACs for all 216 probes. Results to date are consistent with about a 10% failure rate in overgo probe design and a 15-20% false negative detection rate within a group of 216 markers. Three complete rounds of overgo hybridization, each to sets of about 39,000 BACs (either BAMHI or ECORI partial digest inserts) generated a total of 1853 BAC alignments for 517 mapped chicken genome STS markers. These data are publicly available, and they have been used in the assembly of a first generation BAC contig map of the chicken genome.     

Soybean genomic survey: BAC-end sequences near RFLP and SSR markers.

We are building a framework physical infrastructure across the soybean genome by using SSR (simple sequence repeat) and RFLP (restriction fragment length polymorphism) markers to identify BACs (bacterial artificial chromosomes) from two soybean BAC libraries. The libraries were prepared from two genotypes, each digested with a different restriction enzyme. The BACs identified by each marker were grouped into contigs. We have obtained BAC- end sequence from BACs within each contig. The sequences were analyzed by the University of Minnesota Center for Computational Genomics and Bioinformatics using BLAST algorithms to search nucleotide and protein databases. The SSR-identified BACs had a higher percentage of significant BLAST hits than did the RFLP-identified BACs. This difference was due to a higher percentage of hits to repetitive-type sequences for the SSR-identified BACs that was offset in part, however, by a somewhat larger proportion of RFLP-identified significant hits with similarity to experimentally defined genes and soybean ESTs (expressed sequence tags). These genes represented a wide range of metabolic functions. In these analyses, only repetitive sequences from SSR-identified contigs appeared to be clustered. The BAC-end sequences also allowed us to identify microsynteny between soybean and the model plants Arabidopsis thaliana and Medicago truncatula. This map-based approach to genome sampling provides a means of assaying soybean genome structure and organization.     

Construction of a physical and transcript map for a 1-Mb genomic region containing the urofacial (Ochoa) syndrome gene on 10q23-q24 and localization of the disease gene within two overlapping BAC clones (<360 kb).

Urofacial (Ochoa) syndrome is an autosomal recessive disease characterized by distorted facial expression and urinary abnormalities. Previously, we mapped the UFS gene to chromosome 10q23-q24 and narrowed the interval to one YAC clone of 1410 kb. Here, we have constructed a BAC/PAC contig of the 1-Mb region using STS content mapping with 42 BAC/PAC-end sequences, 9 previously reported and 16 newly identified microsatellite markers, and 14 EST markers. A total of 26 polymorphic microsatellite markers were genotyped for 31 UFS patients from Colombia and 2 patients from the United States. Haplotype analyses suggest that the UFS gene is located within two overlapping BAC clones, a region of <360 kb of DNA sequence. We tested 42 EST markers previously mapped to the D10S1709-D10S603 interval against the BAC/PAC contig and identified 11 ESTs located in the 1-Mb region. Four of the 11 ESTs mapped to the 360-kb UFS critical region. Shotgun sequencing of the two BAC clones and BLASTN search of the EST databases revealed 3 other ESTs contained in the UFS critical region. These results will facilitate the cloning and identification of the UFS gene.     

An integrated map of human 6q22.3-q24 including a 3-Mb high-resolution BAC/PAC contig encompassing a QTL for fetal hemoglobin

Genetic studies have previously assigned a quantitative trait locus (QTL) for hemoglobin F and F cells to a region of approximately 4 Mb between the markers D6S408 and D6S292 on chromosome 6q23. An initial yeast artificial chromosome contig of 13 clones spanning this region was generated. Further linkage analysis of an extended kindred refined the candidate interval to 1-2 cM, and key recombination events now place the QTL within a region of <800 kb. We describe a high-resolution bacterial clone contig spanning 3 Mb covering this critical region. The map consists of 223 bacterial artificial chromosome (BAC) and 100 P1 artificial chromosome (PAC) clones ordered by sequence-tagged site (STS) content and restriction fragment fingerprinting with a minimum tiling path of 22 BACs and 1 PAC. A total of 194 STSs map to this interval of 3 Mb, giving an average marker resolution of approximately one per 15 kb. About half of the markers were novel and were isolated in the present study, including three CA repeats and 13 single nucleotide polymorphisms. Altogether 24 expressed sequence tags, 6 of which are unique genes, have been mapped to the contig.     

Estimates of conserved microsynteny among the genomes of <Emphasis Type="Italic">Glycine max</Emphasis>, <Emphasis Type="Italic">Medicago truncatula</Emphasis> and <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

A growing body of research indicates that microsynteny is common among dicot genomes. However, most studies focus on just one or a few genomic regions, so the extent of microsynteny across entire genomes remains poorly characterized. To estimate the level of microsynteny between Medicago truncatula (Mt) and Glycine max (soybean), and also among homoeologous segments of soybean, we used a hybridization strategy involving bacterial artificial chromosome (BAC) contigs. A Mt BAC library consisting of 30,720 clones was screened with a total of 187 soybean BAC subclones and restriction fragment length polymorphism (RFLP) probes. These probes came from 50 soybean contig groups, defined as one or more related BAC contigs anchored by the same low-copy probe. In addition, 92 whole soybean BAC clones were hybridized to filters of HindIII-digested Mt BAC DNA to identify additional cases of cross-hybridization after removal of those soybean BACs found to be repetitive in Mt. Microsynteny was inferred when at least two low-copy probes from a single soybean contig hybridized to the same Mt BAC or when a soybean BAC clone hybridized to three or more low-copy fragments from a single Mt BAC. Of the 50 soybean contig groups examined, 54% showed microsynteny to Mt. The degree of conservation among 37 groups of soybean contigs was also investigated. The results indicated substantial conservation among soybean contigs in the same group, with 86.5% of the groups showing at least some level of microsynteny. One contig group was examined in detail by a combination of physical mapping and comparative sequencing of homoeologous segments. A TBLASTX similarity search was performed between 1,085 soybean sequences on the 50 BAC contig groups and the entire Arabidopsis genome. Based on a criterion of sequence homologues <100 kb apart, each with an expected value of < or =1e-07, seven of the 50 soybean contig groups (14%) exhibited microsynteny with Arabidopsis.     

TOPAAS, a tomato and potato assembly assistance system for selection and finishing of bacterial artificial chromosomes

We have developed the software package Tomato and Potato Assembly Assistance System (TOPAAS), which automates the assembly and scaffolding of contig sequences for low-coverage sequencing projects. The order of contigs predicted by TOPAAS is based on read pair information; alignments between genomic, expressed sequence tags, and bacterial artificial chromosome (BAC) end sequences; and annotated genes. The contig scaffold is used by TOPAAS for automated design of nonredundant sequence gap-flanking PCR primers. We show that TOPAAS builds reliable scaffolds for tomato (Solanum lycopersicum) and potato (Solanum tuberosum) BAC contigs that were assembled from shotgun sequences covering the target at 6- to 8-fold coverage. More than 90% of the gaps are closed by sequence PCR, based on the predicted ordering information. TOPAAS also assists the selection of large genomic insert clones from BAC libraries for walking. For this, tomato BACs are screened by automated BLAST analysis and in parallel, high-density nonselective amplified fragment length polymorphism fingerprinting is used for constructing a high-resolution BAC physical map. BLAST and amplified fragment length polymorphism analysis are then used together to determine the precise overlap. Assembly onto the seed BAC consensus confirms the BACs are properly selected for having an extremely short overlap and largest extending insert. This method will be particularly applicable where related or syntenic genomes are sequenced, as shown here for the Solanaceae, and potentially useful for the monocots Brassicaceae and Leguminosea.     

Diversity and evolution of a non-TIR-NBS sequence family that clusters to a chromosomal ’’hotspot” for disease resistance genes in soybean

In soybean, genes controlling resistance to numerous diseases have been shown to cluster to regions on several chromosomes. One such vital chromosomal region is on the soybean molecular linkage group (MLG) F flanked by the RFLP markers K644 and B212. Here, genes controlling resistance to bacterial blight, Phytophthora root rot, and several viral diseases, as well as QTLs conditioning resistance to corn earworm, root knot nematode, and white mold have been mapped. We have previously identified two classes (b and j) of disease resistance-related nucleotide binding site (NBS) sequences that localize to this cluster of genes. Using both cDNA and genomic analyses, we have studied one multi-gene family of sequences representing the previously reported class j NBS of soybean. This class of NBS resembles the RPS2-like NBS sequences. RPS2 and similar resistance genes are referred to as non-TIR because they do not encode motifs homologous to the Toll-Interleukin-1 region (TIR). By designing PCR primers that specifically target these non-TIR-NBS encoding sequences, we have amplified at least six class j sequence members from soybean. In addition, we have conducted genomic and cDNA library screenings to identify additional class j members. In all, we have characterized 12 class j NBS sequence members. These members have been mapped within a 2-cM region of the soybean F linkage group. We have also identified homoeologous chromosomal regions on linkage groups A2 and E that contain class j NBS sequences. A BLAST search of the GenBank database has shown that non-TIR NBS sequences are present across the legume family. We have compared these non-TIR sequences from other legumes with the soybean clones to assess the level of diversity within this class of disease resistance-related sequences.