An anchored framework BAC map of mouse chromosome 11 assembled using multiplex oligonucleotide hybridization

A pair of gametologous genes provides further insights into avian comparative cytogenomics

Exploration of avian gametologous genes, i.e., homologous genes located on both the Z and W chromosomes, provides a crucial information about the underlying mechanism pertaining to the evolution of these chromosomes. The domestic chicken (Gallus gallus (Linnaeus 1758); GGA) traditionally serves as the primary reference subject of these comparative cytogenomic studies. Using bioinformatic, molecular (overgo BAC library scanning), and cytogenetic (BAC-based FISH) techniques, we have investigated in detail a pair of UBE2R2/UBE2R2L gametologs. By screening a gridded genomic jungle fowl BAC library, CHORI-261, with a short labeled UBE2R2L gene fragment called overgo probe, we detected seven specific clones. For three of them, CH261-019I23, CH261-105E16, and CH261-114G22, we identified their precise cytogenetic location on the Gallus gallus W chromosome (GGAW). They also co-localized with the UBAP2L2 gene on the, as was shown previously, along with the CH261-053P09 BAC clone also containing the GGAW-specific UBE2R2L DNA sequence. The fine mapping of the UBE2R2/UBE2R2L homologs in the chicken genome also shed the light on comparative cytogenetic aspects in birds. Our findings provided further evidence that bird genomes moderately changed only during evolution and are suitable for successful use of interspecies hybridization using both overgo-based BAC library screen and BAC-based FISH.

Preparation of high molecular weight gDNA and bacterial artificial chromosome (BAC) libraries in plants

Bacterial artificial chromosome (BAC) libraries are extremely valuable large-insert DNA libraries for physical mapping, positional cloning, comparative genomic analysis, complete genome sequencing, and evolutionary studies. Due to their stability and relative simplicity BAC libraries are most preferred over other approaches for cloning large genomic DNA fragments for large-insert libraries. Isolation of intact high molecular weight (HMW) DNA is a critical step underlying the success of large-insert genomic DNA library construction. It requires the isolation of purified nuclei, embedding them into LMP agarose plugs, restriction digestion of the plugs, and quite often size selection using PFGE and electro-elution of insert DNA. The construction of BAC libraries is complex and challenging for most molecular laboratories. To facilitate the construction of BAC libraries, we present a step-by-step protocol for isolation of HMW DNA and construction of plant BAC libraries.

A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor

The current limitations in genome sequencing technology require the construction of physical maps for high-quality draft sequences of large plant genomes, such as that of Aegilops tauschii, the wheat D-genome progenitor. To construct a physical map of the Ae. tauschii genome, we fingerprinted 461,706 bacterial artificial chromosome clones, assembled contigs, designed a 10K Ae. tauschii Infinium SNP array, constructed a 7,185-marker genetic map, and anchored on the map contigs totaling 4.03 Gb. Using whole genome shotgun reads, we extended the SNP marker sequences and found 17,093 genes and gene fragments. We showed that collinearity of the Ae. tauschii genes with Brachypodium distachyon, rice, and sorghum decreased with phylogenetic distance and that structural genome evolution rates have been high across all investigated lineages in subfamily Pooideae, including that of Brachypodieae. We obtained additional information about the evolution of the seven Triticeae chromosomes from 12 ancestral chromosomes and uncovered a pattern of centromere inactivation accompanying nested chromosome insertions in grasses. We showed that the density of noncollinear genes along the Ae. tauschii chromosomes positively correlates with recombination rates, suggested a cause, and showed that new genes, exemplified by disease resistance genes, are preferentially located in high-recombination chromosome regions.

Physical mapping of a pollen modifier locus controlling self-incompatibility in apricot and synteny analysis within the Rosaceae

S-locus products (S-RNase and F-box proteins) are essential for the gametophytic self-incompatibility (GSI) specific recognition in Prunus. However, accumulated genetic evidence suggests that other S-locus unlinked factors are also required for GSI. For instance, GSI breakdown was associated with a pollen-part mutation unlinked to the S-locus in the apricot (Prunus armeniaca L.) cv. 'Canino'. Fine-mapping of this mutated modifier gene (M-locus) and the synteny analysis of the M-locus within the Rosaceae are here reported. A segregation distortion loci mapping strategy, based on a selectively genotyped population, was used to map the M-locus. In addition, a bacterial artificial chromosome (BAC) contig was constructed for this region using overlapping oligonucleotides probes, and BAC-end sequences (BES) were blasted against Rosaceae genomes to perform micro-synteny analysis. The M-locus was mapped to the distal part of chr.3 flanked by two SSR markers within an interval of 1.8 cM corresponding to ~364 Kb in the peach (Prunus persica L. Batsch) genome. In the integrated genetic-physical map of this region, BES were mapped against the peach scaffold_3 and BACs were anchored to the apricot map. Micro-syntenic blocks were detected in apple (Malus × domestica Borkh.) LG17/9 and strawberry (Fragaria vesca L.) FG6 chromosomes. The M-locus fine-scale mapping provides a solid basis for self-compatibility marker-assisted selection and for positional cloning of the underlying gene, a necessary goal to elucidate the pollen rejection mechanism in Prunus. In a wider context, the syntenic regions identified in peach, apple and strawberry might be useful to interpret GSI evolution in Rosaceae.

Palaeohexaploid ancestry for Caryophyllales inferred from extensive gene-based physical and genetic mapping of the sugar beet genome (Beta vulgaris)

Sugar beet (Beta vulgaris) is an important crop plant that accounts for 30% of the world's sugar production annually. The genus Beta is a distant relative of currently sequenced taxa within the core eudicotyledons; the genomic characterization of sugar beet is essential to make its genome accessible to molecular dissection. Here, we present comprehensive genomic information in genetic and physical maps that cover all nine chromosomes. Based on this information we identified the proposed ancestral linkage groups of rosids and asterids within the sugar beet genome. We generated an extended genetic map that comprises 1127 single nucleotide polymorphism markers prepared from expressed sequence tags and bacterial artificial chromosome (BAC) end sequences. To construct a genome-wide physical map, we hybridized gene-derived oligomer probes against two BAC libraries with 9.5-fold cumulative coverage of the 758 Mbp genome. More than 2500 probes and clones were integrated both in genetic maps and the physical data. The final physical map encompasses 535 chromosomally anchored contigs that contains 8361 probes and 22 815 BAC clones. By using the gene order established with the physical map, we detected regions of synteny between sugar beet (order Caryophyllales) and rosid species that involves 1400-2700 genes in the sequenced genomes of Arabidopsis, poplar, grapevine, and cacao. The data suggest that Caryophyllales share the palaeohexaploid ancestor proposed for rosids and asterids. Taken together, we here provide extensive molecular resources for sugar beet and enable future high-resolution trait mapping, gene identification, and cross-referencing to regions sequenced in other plant species.

Characterization of the genome of bald cypress

BACKGROUND

Bald cypress (Taxodium distichum var. distichum) is a coniferous tree of tremendous ecological and economic importance. It is a member of the family Cupressaceae which also includes cypresses, redwoods, sequoias, thujas, and junipers. While the bald cypress genome is more than three times the size of the human genome, its 1C DNA content is amongst the smallest of any conifer. To learn more about the genome of bald cypress and gain insight into the evolution of Cupressaceae genomes, we performed a Cot analysis and used Cot filtration to study Taxodium DNA. Additionally, we constructed a 6.7 genome-equivalent BAC library that we screened with known Taxodium genes and select repeats.

RESULTS

The bald cypress genome is composed of 90% repetitive DNA with most sequences being found in low to mid copy numbers. The most abundant repeats are found in fewer than 25,000 copies per genome. Approximately 7.4% of the genome is single/low-copy DNA (i.e., sequences found in 1 to 5 copies). Sequencing of highly repetitive Cot clones indicates that most Taxodium repeats are highly diverged from previously characterized plant repeat sequences. The bald cypress BAC library consists of 606,336 clones (average insert size of 113 kb) and collectively provides 6.7-fold genome equivalent coverage of the bald cypress genome. Macroarray screening with known genes produced, on average, about 1.5 positive clones per probe per genome-equivalent. Library screening with Cot-1 DNA revealed that approximately 83% of BAC clones contain repetitive sequences iterated 103 to 104 times per genome.

CONCLUSIONS

The BAC library for bald cypress is the first to be generated for a conifer species outside of the family Pinaceae. The Taxodium BAC library was shown to be useful in gene isolation and genome characterization and should be an important tool in gymnosperm comparative genomics, physical mapping, genome sequencing, and gene/polymorphism discovery. The single/low-copy (SL) component of bald cypress is 4.6 times the size of the Arabidopsis genome. As suggested for other gymnosperms, the large amount of SL DNA in Taxodium is likely the result of divergence among ancient repeat copies and gene/pseudogene duplication.

Narrowing down the apricot Plum pox virus resistance locus and comparative analysis with the peach genome syntenic region

Sharka disease, caused by the Plum pox virus (PPV), is one of the main limiting factors for stone fruit crops worldwide. Only a few resistance sources have been found in apricot (Prunus armeniaca L.), and most studies have located a major PPV resistance locus (PPVres) on linkage group 1 (LG1). However, the mapping accuracy was not sufficiently reliable and PPVres was predicted within a low confidence interval. In this study, we have constructed two high-density simple sequence repeat (SSR) improved maps with 0.70 and 0.68 markers/cm, corresponding to LG1 of 'Lito' and 'Goldrich' PPV-resistant cultivars, respectively. Using these maps, and excluding genotype-phenotype incongruent individuals, a new binary trait locus (BTL) analysis for PPV resistance was performed, narrowing down the PPVres support intervals to 7.3 and 5.9 cm in 'Lito' and 'Goldrich', respectively. Subsequently, 71 overlapping oligonucleotides (overgo) probes were hybridized against an apricot bacterial artificial chromosome (BAC) library, identifying 870 single BACs from which 340 were anchored onto a map region of approximately 30-40 cm encompassing PPVres. Partial BAC contigs assigned to the two allelic haplotypes (resistant/susceptible) of the PPVres locus were built by high-information content fingerprinting (HICF). In addition, a total of 300 BAC-derived sequences were obtained, and 257 showed significant homology with the peach genome scaffold_1 corresponding to LG1. According to the peach syntenic genome sequence, PPVres was predicted within a region of 2.16 Mb in which a few candidate resistance genes were identified.

Polyploidization altered gene functions in cotton (Gossypium spp.)

Cotton (Gossypium spp.) is an important crop plant that is widely grown to produce both natural textile fibers and cottonseed oil. Cotton fibers, the economically more important product of the cotton plant, are seed trichomes derived from individual cells of the epidermal layer of the seed coat. It has been known for a long time that large numbers of genes determine the development of cotton fiber, and more recently it has been determined that these genes are distributed across At and Dt subgenomes of tetraploid AD cottons. In the present study, the organization and evolution of the fiber development genes were investigated through the construction of an integrated genetic and physical map of fiber development genes whose functions have been verified and confirmed. A total of 535 cotton fiber development genes, including 103 fiber transcription factors, 259 fiber development genes, and 173 SSR-contained fiber ESTs, were analyzed at the subgenome level. A total of 499 fiber related contigs were selected and assembled. Together these contigs covered about 151 Mb in physical length, or about 6.7% of the tetraploid cotton genome. Among the 499 contigs, 397 were anchored onto individual chromosomes. Results from our studies on the distribution patterns of the fiber development genes and transcription factors between the At and Dt subgenomes showed that more transcription factors were from Dt subgenome than At, whereas more fiber development genes were from At subgenome than Dt. Combining our mapping results with previous reports that more fiber QTLs were mapped in Dt subgenome than At subgenome, the results suggested a new functional hypothesis for tetraploid cotton. After the merging of the two diploid Gossypium genomes, the At subgenome has provided most of the genes for fiber development, because it continues to function similar to its fiber producing diploid A genome ancestor. On the other hand, the Dt subgenome, with its non-fiber producing D genome ancestor, provides more transcription factors that regulate the expression of the fiber genes in the At subgenome. This hypothesis would explain previously published mapping results. At the same time, this integrated map of fiber development genes would provide a framework to clone individual full-length fiber genes, to elucidate the physiological mechanisms of the fiber differentiation, elongation, and maturation, and to systematically study the functional network of these genes that interact during the process of fiber development in the tetraploid cottons.

The Tomato Sequencing Project, the first cornerstone of the International Solanaceae Project (SOL)

The genome of tomato (Solanum lycopersicum) is being sequenced by an international consortium of 10 countries (Korea, China, the United Kingdom, India, The Netherlands, France, Japan, Spain, Italy and the United States) as part of a larger initiative called the ‘International Solanaceae Genome Project (SOL): Systems Approach to Diversity and Adaptation’. The goal of this grassroots initiative, launched in November 2003, is to establish a network of information, resources and scientists to ultimately tackle two of the most significant questions in plant biology and agriculture: (1) How can a common set of genes/proteins give rise to a wide range of morphologically and ecologically distinct organisms that occupy our planet? (2) How can a deeper understanding of the genetic basis of plant diversity be harnessed to better meet the needs of society in an environmentally friendly and sustainable manner? The Solanaceae and closely related species such as coffee, which are included in the scope of the SOL project, are ideally suited to address both of these questions. The first step of the SOL project is to use an ordered BAC approach to generate a high quality sequence for the euchromatic portions of the tomato as a reference for the Solanaceae. Due to the high level of macro and micro-synteny in the Solanaceae the BAC-by-BAC tomato sequence will form the framework for shotgun sequencing of other species. The starting point for sequencing the genome is BACs anchored to the genetic map by overgo hybridization and AFLP technology. The overgos are derived from approximately 1500 markers from the tomato high density F2-2000 genetic map (http://sgn.cornell.edu/). These seed BACs will be used as anchors from which to radiate the tiling path using BAC end sequence data. Annotation will be performed according to SOL project guidelines. All the information generated under the SOL umbrella will be made available in a comprehensive website. The information will be interlinked with the ultimate goal that the comparative biology of the Solanaceae—and beyond—achieves a context that will facilitate a systems biology approach.

A draft physical map of a D-genome cotton species (Gossypium raimondii)

BACKGROUND

Genetically anchored physical maps of large eukaryotic genomes have proven useful both for their intrinsic merit and as an adjunct to genome sequencing. Cultivated tetraploid cottons, Gossypium hirsutum and G. barbadense, share a common ancestor formed by a merger of the A and D genomes about 1-2 million years ago. Toward the long-term goal of characterizing the spectrum of diversity among cotton genomes, the worldwide cotton community has prioritized the D genome progenitor Gossypium raimondii for complete sequencing.

RESULTS

A whole genome physical map of G. raimondii, the putative D genome ancestral species of tetraploid cottons was assembled, integrating genetically-anchored overgo hybridization probes, agarose based fingerprints and 'high information content fingerprinting' (HICF). A total of 13,662 BAC-end sequences and 2,828 DNA probes were used in genetically anchoring 1585 contigs to a cotton consensus genetic map, and 370 and 438 contigs, respectively to Arabidopsis thaliana (AT) and Vitis vinifera (VV) whole genome sequences.

CONCLUSION

Several lines of evidence suggest that the G. raimondii genome is comprised of two qualitatively different components. Much of the gene rich component is aligned to the Arabidopsis and Vitis vinifera genomes and shows promise for utilizing translational genomic approaches in understanding this important genome and its resident genes. The integrated genetic-physical map is of value both in assembling and validating a planned reference sequence.

A high-throughput strategy for screening of bacterial artificial chromosome libraries and anchoring of clones on a genetic map constructed with single nucleotide polymorphisms

BACKGROUND

Current techniques of screening bacterial artificial chromosome (BAC) libraries for molecular markers during the construction of physical maps are slow, laborious and often assign multiple BAC contigs to a single locus on a genetic map. These limitations are the principal impediment in the construction of physical maps of large eukaryotic genomes. It is hypothesized that this impediment can be overcome by screening multidimensional pools of BAC clones using the highly parallel Illumina GoldenGate assay.

RESULTS

To test the efficacy of the Golden Gate assay in BAC library screening, multidimensional pools involving 302976 Aegilops tauschii BAC clones were genotyped for the presence/absence of specific gene sequences with multiplexed Illumina GoldenGate oligonucleotide assays previously used to place single nucleotide polymorphisms on an Ae. tauschii genetic map. Of 1384 allele-informative oligonucleotide assays, 87.6% successfully clustered BAC pools into those positive for a BAC clone harboring a specific gene locus and those negative for it. The location of the positive BAC clones within contigs assembled from 199190 fingerprinted Ae. tauschii BAC clones was used to evaluate the precision of anchoring of BAC clones and contigs on the Ae. tauschii genetic map. For 41 (95%) assays, positive BAC clones were neighbors in single contigs. Those contigs could be unequivocally assigned to loci on the genetic map. For two (5%) assays, positive clones were in two different contigs and the relationships of these contigs to loci on the Ae. tauschii genetic map were equivocal. Screening of BAC libraries with a simple five-dimensional BAC pooling strategy was evaluated and shown to allow direct detection of positive BAC clones without the need for manual deconvolution of BAC clone pools.

CONCLUSION

The highly parallel Illumina oligonucleotide assay is shown here to be an efficient tool for screening BAC libraries and a strategy for high-throughput anchoring of BAC contigs on genetic maps during the construction of physical maps of eukaryotic genomes. In most cases, screening of BAC libraries with Illumina oligonucleotide assays results in the unequivocal relationship of BAC clones with loci on the genetic map.

The mouse hairy ears mutation exhibits an extended growth (anagen) phase in hair follicles and altered Hoxc gene expression in the ears

The mouse In(15)2Rl (hairy ears, Eh) mutation is a paracentric inversion of the distal half of chromosome 15 (Chr 15). Heterozygous Eh/+ mice display misshaped and hairy ears that have more and longer hair than the ears of their wild-type littermates. We mapped, cloned and sequenced both inversion breakpoints. No protein-coding transcript was disrupted by either breakpoint. The proximal breakpoint is located between syntrophin basic 1 (Sntb1) and hyaluronan synthase 2 (Has2), and the distal breakpoint maps between homeobox C4 (Hoxc4) and single-strand selective monofunctional uracil DNA glycosylase (Smug1), near the middle and the telomere ends of Chr 15, respectively. The inversion spans ~47 megabases. Our genetic analysis suggests that the hairy-ear phenotype is caused by the proximal breakpoint of the inversion-bearing Chr 15. Quantitative RNA analysis by real-time polymerase chain reaction for the genes flanking the breakpoint indicated no changes in expression levels except for some homeobox C (Hoxc) genes whose expression was elevated in developing and mature skin of the ears but not of other body regions. The increased hair length on the ears of Eh/+ mice was due to an extension of the anagen stage in the hair cycle, as determined by histological analysis. Our data indicate that the Eh phenotype arises from mis-expression of Hoxc genes.

Gene-rich islands for fiber development in the cotton genome

Cotton fiber is an economically important seed trichome and the world's leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.

An integrated genetic and physical map of homoeologous chromosomes 12 and 26 in Upland cotton (G. hirsutum L.)

Background

Upland cotton (G. hirsutum L.) is the leading fiber crop worldwide. Genetic improvement of fiber quality and yield is facilitated by a variety of genomics tools. An integrated genetic and physical map is needed to better characterize quantitative trait loci and to allow for the positional cloning of valuable genes. However, developing integrated genomic tools for complex allotetraploid genomes, like that of cotton, is highly experimental. In this report, we describe an effective approach for developing an integrated physical framework that allows for the distinguishing between subgenomes in cotton.

Results

A physical map has been developed with 220 and 115 BAC contigs for homeologous chromosomes 12 and 26, respectively, covering 73.49 Mb and 34.23 Mb in physical length. Approximately one half of the 220 contigs were anchored to the At subgenome only, while 48 of the 115 contigs were allocated to the Dt subgenome only. Between the two chromosomes, 67 contigs were shared with an estimated overall physical similarity between the two chromosomal homeologs at 40.0 %. A total of 401 fiber unigenes plus 214 non-fiber unigenes were located to chromosome 12 while 207 fiber unigenes plus 183 non-fiber unigenes were allocated to chromosome 26. Anchoring was done through an overgo hybridization approach and all anchored ESTs were functionally annotated via blast analysis.

Conclusion

This integrated genomic map describes the first pair of homoeologous chromosomes of an allotetraploid genome in which BAC contigs were identified and partially separated through the use of chromosome-specific probes and locus-specific genetic markers. The approach used in this study should prove useful in the construction of genome-wide physical maps for polyploid plant genomes including Upland cotton. The identification of Gene-rich islands in the integrated map provides a platform for positional cloning of important genes and the targeted sequencing of specific genomic regions.

Development of a pooled probe method for locating small gene families in a physical map of soybean using stress related paralogues and a BAC minimum tile path

BACKGROUND

Genome analysis of soybean (Glycine max L.) has been complicated by its paleo-autopolyploid nature and conserved homeologous regions. Landmarks of expressed sequence tags (ESTs) located within a minimum tile path (MTP) of contiguous (contig) bacterial artificial chromosome (BAC) clones or radiation hybrid set can identify stress and defense related gene rich regions in the genome. A physical map of about 2,800 contigs and MTPs of 8,064 BAC clones encompass the soybean genome. That genome is being sequenced by whole genome shotgun methods so that reliable estimates of gene family size and gene locations will provide a useful tool for finishing. The aims here were to develop methods to anchor plant defense- and stress-related gene paralogues on the MTP derived from the soybean physical map, to identify gene rich regions and to correlate those with QTL for disease resistance.

RESULTS

The probes included 143 ESTs from a root library selected by subtractive hybridization from a multiply disease resistant soybean cultivar 'Forrest' 14 days after inoculation with Fusarium solani f. sp. glycines (F. virguliforme). Another 166 probes were chosen from a root EST library (Gm-r1021) prepared from a non-inoculated soybean cultivar 'Williams 82' based on their homology to the known defense and stress related genes. Twelve and thirteen pooled EST probes were hybridized to high-density colony arrays of MTP BAC clones from the cv. 'Forrest' genome. The EST pools located 613 paralogues for 201 of the 309 probes used (range 1-13 per functional probe). One hundred BAC clones contained more than one kind of paralogue. Many more BACs (246) contained a single paralogue of one of the 201 probes detectable gene families. ESTs were anchored on soybean linkage groups A1, B1, C2, E, D1a+Q, G, I, M, H, and O.

CONCLUSION

Estimates of gene family sizes were more similar to those made by Southern hybridization than by bioinformatics inferences from EST collections. When compared to Arabidopsis thaliana there were more 2 and 4 member paralogue families reflecting the diploidized-tetraploid nature of the soybean genome. However there were fewer families with 5 or more genes and the same number of single genes. Therefore the method can identify evolutionary patterns such as massively extensive selective gene loss or rapid divergence to regenerate the unique genes in some families.

An improved method to identify BAC clones using pooled overgos

Hybridization using overgo probes is an established approach for screening arrayed bacterial artificial chromosome (BAC) libraries. We have improved the use of overgos by increasing the yield of positive clones using reduced levels of radioisotopes and enzyme. The strategy involves labeling with all four radiolabeled nucleotides in a hot pulse followed by a cold nucleotide chase and then extending the exposure time to compensate for reduced specific activity of the probes. The resulting cost savings and reduced human exposure to radiation make the use of highly pooled overgo probes a more attractive approach for screening of BAC libraries from organisms with large genomes.

Integration of hybridization-based markers (overgos) into physical maps for comparative and evolutionary explorations in the genus Oryza and in Sorghum

BACKGROUND

With the completion of the genome sequence for rice (Oryza sativa L.), the focus of rice genomics research has shifted to the comparison of the rice genome with genomes of other species for gene cloning, breeding, and evolutionary studies. The genus Oryza includes 23 species that shared a common ancestor 8-10 million years ago making this an ideal model for investigations into the processes underlying domestication, as many of the Oryza species are still undergoing domestication. This study integrates high-throughput, hybridization-based markers with BAC end sequence and fingerprint data to construct physical maps of rice chromosome 1 orthologues in two wild Oryza species. Similar studies were undertaken in Sorghum bicolor, a species which diverged from cultivated rice 40-50 million years ago.

RESULTS

Overgo markers, in conjunction with fingerprint and BAC end sequence data, were used to build sequence-ready BAC contigs for two wild Oryza species. The markers drove contig merges to construct physical maps syntenic to rice chromosome 1 in the wild species and provided evidence for at least one rearrangement on chromosome 1 of the O. sativa versus Oryza officinalis comparative map. When rice overgos were aligned to available S. bicolor sequence, 29% of the overgos aligned with three or fewer mismatches; of these, 41% gave positive hybridization signals. Overgo hybridization patterns supported colinearity of loci in regions of sorghum chromosome 3 and rice chromosome 1 and suggested that a possible genomic inversion occurred in this syntenic region in one of the two genomes after the divergence of S. bicolor and O. sativa.

CONCLUSION

The results of this study emphasize the importance of identifying conserved sequences in the reference sequence when designing overgo probes in order for those probes to hybridize successfully in distantly related species. As interspecific markers, overgos can be used successfully to construct physical maps in species which diverged less than 8 million years ago, and can be used in a more limited fashion to examine colinearity among species which diverged as much as 40 million years ago. Additionally, overgos are able to provide evidence of genomic rearrangements in comparative physical mapping studies.

Construction of BAC and BIBAC libraries from sunflower and identification of linkage group-specific clones by overgo hybridization

Complementary BAC and BIBAC libraries were constructed from nuclear DNA of sunflower cultivar HA 89. The BAC library, constructed with BamHI in the pECBAC1 vector, contains 107,136 clones and has an average insert size of 140 kb. The BIBAC library was constructed with HindIII in the plant-transformation-competent binary vector pCLD04541 and contains 84,864 clones, with an average insert size of 137 kb. The two libraries combined contain 192,000 clones and are equivalent to approximately 8.9 haploid genomes of sunflower (3,000 Mb/1C), and provide a greater than 99% probability of obtaining a clone of interest. The frequencies of BAC and BIBAC clones carrying chloroplast or mitochondrial DNA sequences were estimated to be 2.35 and 0.04%, respectively, and insert-empty clones were less than 0.5%. To facilitate chromosome engineering and anchor the sunflower genetic map to its chromosomes, one to three single- or low-copy RFLP markers from each linkage group of sunflower were used to design pairs of overlapping oligonucleotides (overgos). Thirty-six overgos were designed and pooled as probes to screen a subset (5.1x) of the BAC and BIBAC libraries. Of the 36 overgos, 33 (92%) gave at least one positive clone and 3 (8%) failed to hit any clone. As a result, 195 BAC and BIBAC clones representing 19 linkage groups were identified, including 76 BAC clones and 119 BIBAC clones, further verifying the genome coverage and utility of the libraries. These BAC and BIBAC libraries and linkage group-specific clones provide resources essential for comprehensive research of the sunflower genome.

Construction of a medium-density horse gene map

A medium-density map of the horse genome (Equus caballus) was constructed using genes evenly distributed over the human genome. Three hundred and twenty-three exonic primer pairs were used to screen the INRA and the CHORI-241 equine BAC libraries by polymerase chain reaction and by filter hybridization respectively. Two hundred and thirty-seven BACs containing equine gene orthologues, confirmed by sequencing, were isolated. The BACs were localized to horse chromosomes by fluorescent in situ hybridization (FISH). Overall, 165 genes were assigned to the equine genomic map by radiation hybrid (RH) (using an equine RH(5000) panel) and/or by FISH mapping. A comparison of localizations of 713 genes mapped on the horse genome and on the human genome revealed 59 homologous segments and 131 conserved segments. Two of these homologies (ECA27/HSA8 and ECA12p/HSA11p) had not been previously identified. An enhanced resolution of conserved and rearranged chromosomal segments presented in this study provides clarification of chromosome evolution history.

Characterization of 23 CC chemokine genes and analysis of their expression in channel catfish (Ictalurus punctatus)

Chemokines are a large family of chemotactic cytokines playing crucial roles in the innate immune response. CC chemokines constitute the largest subfamily of chemokines, with 28 CC chemokines identified from mammalian species. However, the status of CC chemokines in teleosts is yet to be determined. We previously identified 26 catfish CC chemokine cDNAs from catfish. In this study, we isolated and sequenced 23 channel catfish CC chemokine genes amounting to a total of over 56 kb of genomic sequences. Genomic organization of the 23 CC chemokine genes was determined by comparing the generated genomic sequences with the previously identified cDNA sequences. Microsatellites were identified from 16 catfish CC chemokine genes allowing them to be utilized for genome mapping. Structural analysis indicated conservation of genomic organization of CC chemokine genes, which may facilitate the establishment of orthologies. Expression of all known catfish CC chemokine transcripts was assessed in nine important tissues. Of the 26 catfish CC chemokine genes, 14 were universally expressed, six were widely expressed in many tissues, while six were highly tissue-specific.

Catfish CC chemokines: genomic clustering, duplications, and expression after bacterial infection with Edwardsiella ictaluri

Chemokines are a family of structurally related chemotactic cytokines that regulate the migration of leukocytes, under both physiological and inflammatory conditions. CC chemokines represent the largest subfamily of chemokines with 28 genes in mammals. Sequence conservation of chemokines between teleost fish and higher vertebrates is low and duplication and divergence may have occurred at a significantly faster rate than in other genes. One feature of CC chemokine genes known to be conserved is genomic clustering. CC chemokines are highly clustered within the genomes of human, mouse, and chicken. To exploit knowledge from comparative genome analysis between catfish and higher vertebrates, here we mapped to bacterial artificial chromosome (BAC) clones 26 previously identified catfish (Ictalurus sp.) chemokine cDNAs. Through a combination of hybridization and fluorescent fingerprinting, 18 fingerprinted contigs were assembled from BACs containing catfish CC chemokine genes. The catfish CC chemokine genes were found to be not only highly clustered in the catfish genome, but also extensively duplicated at various levels. Comparisons of the syntenic relationships of CC chemokines may help to explain the modes of duplication and divergence that resulted in the present repertoire of vertebrate CC chemokines. Here we have also analyzed the expression of the transcripts of the 26 catfish CC chemokines in head kidney and spleen in response to bacterial infection of Edwardsiella ictaluri, an economically devastating catfish pathogen. Such information should pinpoint research efforts on the CC chemokines most likely involved in inflammatory responses.

Comparative physical mapping links conservation of microsynteny to chromosome structure and recombination in grasses

Nearly finished sequences for model organisms provide a foundation from which to explore genomic diversity among other taxonomic groups. We explore genome-wide microsynteny patterns between the rice sequence and two sorghum physical maps that integrate genetic markers, bacterial artificial chromosome (BAC) fingerprints, and BAC hybridization data. The sorghum maps largely tile a genomic component containing 41% of BACs but 80% of single-copy genes that shows conserved microsynteny with rice and partially tile a nonsyntenic component containing 46% of BACs but only 13% of single-copy genes. The remaining BACs are centromeric (4%) or unassigned (8%). The two genomic components correspond to cytologically discernible "euchromatin" and "heterochromatin." Gene and repetitive DNA distributions support this classification. Greater microcolinearity in recombinogenic (euchromatic) than nonrecombinogenic (heterochromatic) regions is consistent with the hypothesis that genomic rearrangements are usually deleterious, thus more likely to persist in nonrecombinogenic regions by virtue of Muller's ratchet. Interchromosomal centromeric rearrangements may have fostered diploidization of a polyploid cereal progenitor. Model plant sequences better guide studies of related genomes in recombinogenic than nonrecombinogenic regions. Bridging of 35 physical gaps in the rice sequence by sorghum BAC contigs illustrates reciprocal benefits of comparative approaches that extend at least across the cereals and perhaps beyond.

Exploratory integration of peanut genetic and physical maps and possible contributions from Arabidopsis

Arachis hypogaea is a widely cultivated crop both as an oilseed and protein source. The genomic analysis of Arachis species hitherto has been limited to the construction of genetic maps; the most comprehensive one contains 370 loci over 2,210 cM in length. However, no attempt has been made to analyze the physical structure of the peanut genome. To investigate the practicality of physical mapping in peanut, we applied a total of 117 oligonucleotide-based probes ("overgos") derived from genetically mapped RFLP probes onto peanut BAC filters containing 182,784 peanut large-insert DNA clones in a multiplex experimental design; 91.5% of the overgos identified at least one BAC clone. In order to gain insights into the potential value of Arabidopsis genome sequence for studies in divergent species with complex genomes such as peanut, we employed 576 Arabidopsis-derived overgos selected on the basis of maximum homology to orthologous sequences in other plant taxa to screen the peanut BAC library. A total of 353 (61.3%) overgos detected at least one peanut BAC clone. This experiment represents the first steps toward the creation of a physical map in peanut and illustrates the potential value of leveraging information from distantly related species such as Arabidopsis for both practical applications such as comparative map-based cloning and shedding light on evolutionary relationships. We also evaluated the possible correlation between functional categories of Arabidopsis overgos and their success rates in hybridization to the peanut BAC library.

Mapping segmental and sequence variations among laboratory mice using BAC array CGH

We used arrays of 2069 BACs (1303 nonredundant autosomal clones) to map sequence variation among Mus spretus (SPRET/Ei and SPRET/Glasgow) and Mus musculus (C3H/HeJ, BALB/cJ, 129/J, DBA/2J, NIH, FVB/N, and C57BL/6) strains. We identified 80 clones representing 74 autosomal loci of copy number variation (|log(2)ratio| >/= 0.4). These variant loci distinguish laboratory strains. By FISH mapping, we determined that 63 BACs mapped to a single site on C57BL/6J chromosomes, while 17 clones mapped to multiple chromosomes (n = 16) or multiple sites on one chromosome (n = 1). We also show that small ratio changes (Delta log(2)ratio approximately 0.1) distinguish homozygous and heterozygous regions of the genome in interspecific backcross mice, providing an efficient method for genotyping progeny of backcrosses.

Serial segmental duplications during primate evolution result in complex human genome architecture

The human genome is particularly rich in low-copy repeats (LCRs) or segmental duplications (5%-10%), and this characteristic likely distinguishes us from lower mammals such as rodents. How and why the complex human genome architecture consisting of multiple LCRs has evolved remains an open question. Using molecular and computational analyses of human and primate genomic regions, we analyzed the structure and evolution of LCRs that resulted in complex architectural features of the human genome in proximal 17p. We found that multiple LCRs of different origins are situated adjacent to one another, whereas each LCR changed at different time points between >25 to 3-7 million years ago (Mya) during primate evolution. Evolutionary studies in primates suggested communication between the LCRs by gene conversion. The DNA transposable element MER1-Charlie3 and retroviral ERVL elements were identified at the breakpoint of the t(4;19) chromosome translocation in Gorilla gorilla, suggesting a potential role for transpositions in evolution of the primate genome. Thus, a series of consecutive segmental duplication events during primate evolution resulted in complex genome architecture in proximal 17p. Some of the more recent events led to the formation of novel genes that in human are expressed primarily in the brain. Our observations support the contention that serial segmental duplication events might have orchestrated primate evolution by the generation of novel fusion/fission genes as well as potentially by genomic inversions associated with decreased recombination rates facilitating gene divergence.

Comparative analysis of early embryonic sunflower cDNA libraries

To gain information concerning cell functions and activities during sunflower embryogenesis, an expressed sequence tag (EST) approach was used to analyse gene expression in the early stages of sunflower embryos development. Confocal microscopy observations of whole-mounted embryos allowed us to identify precisely the major steps of the zygotic embryonic development. A time-course analysis was then employed to collect the embryonic material. Three cDNA libraries were constructed from microdissected embryos, and three other cDNA libraries were created using a classical day after pollination schedule. A total of 7106 ESTs were produced and assembled. The total number of putative different genes represents about 43.1 (3064 tentative contigs and singlets) of the analysed sequences. The unigenes that showed similarity to proteins with known or predicted functions (50.3) were classified into 15 different functional categories. The functional profiles were found to be quite similar for all studied embryo stages but statistical analysis revealed that successive and coordinate sets of genes are expressed at each embryonic stage. The analysis allowed us to identify abundant and differentially expressed genes at the early stages of embryos development as well as some putatively interesting genes, showing strong similarities with genes playing key roles in plant and animal embryogenesis. The data presented in this study not only provide a first global overview of the genes expression profile during sunflower embryogenesis but also represent an original and valuable tool for developmental genomics studies on exalbuminous dicots.

Catfish hepcidin gene is expressed in a wide range of tissues and exhibits tissue-specific upregulation after bacterial infection

Antimicrobial peptides (AMPs) are important components of the host innate immune response against microbial invasion. The cysteine-rich AMPs such as defensin and hepcidin have been extensively studied from various organisms, but their role in disease defense in catfish is unknown. As a first step, we sequenced a hepcidin cDNA from both channel catfish and blue catfish, and characterized the channel catfish hepcidin gene. The channel catfish hepcidin gene consists of two introns and three exons that encode a peptide of 96 amino acids. The amino acid sequences and gene organization were conserved between catfish and other organisms. In contrast to its almost exclusive expression in the liver in humans, the channel catfish hepcidin gene was expressed in a wide range of tissues except brain. Its expression was detected early during embryonic and larval development, and induced after bacterial infection with Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC) in a tissue-specific manner. The upregulation was observed in the spleen and head kidney, but not in the liver. The expression of hepcidin was upregulated 1--3 days after challenge, but returned to normal levels at 7 days after challenge. The expression profile of the catfish hepcidin gene during the course of bacterial infection mirrors those of inflammatory proteins such as chemokines, suggesting an important role for hepcidin during inflammatory responses.

Comparative genome analysis of monocots and dicots, toward characterization of angiosperm diversity

The importance of angiosperms to sustaining humanity by providing a wide range of 'ecosystem services' warrants increased exploration of their genomic diversity. The nearly completed sequences for two species representing the major angiosperm subclasses, specifically the dicot Arabidopsis thaliana and the monocot Oryza sativa, provide a foundation for comparative analysis across the angiosperms. The angiosperms also exemplify some challenges to be faced as genomics makes new inroads into describing biotic diversity, in particular polyploidy (genome-wide chromatin duplication), and much larger genome sizes than have been studied to date.

The neurological mutant quaking(viable) is Parkin deficient

The mouse mutant quaking(viable) ( qk(v)) has been studied for almost four decades as a model for dysmyelination of the central nervous system (CNS). The genetic lesion associated with the qk(v) phenotype is a large deletion of approximately 1 Megabase on mouse Chromosome (Chr) 17. This deficiency alters the expression of transcripts from the qkI locus in oligodendrocytes, resulting in improper myelination of the CNS in animals homozygous for the deletion. To determine whether other genes within the deletion contribute to the quaking(viable) phenotype, we physically mapped and sequenced the deleted interval. We determined that the mouse Parkin gene, as well as the Parkin co-regulated gene ( Pacrg), lies within the qk(v) deletion. We determined that qk(v) mutants completely lack the expression of the Parkin gene product. Loss-of-function mutations in the human PARKIN gene cause autosomal juvenile Parkinson's disease (AR-JP). Our studies show that the deletion of Parkin in qk(v) brains does not result in the loss of dopaminergic neurons typical of AR-JP patients. Also, alpha-synuclein, a target of Parkin-dependent ubiquitination, does not accumulate in qk(v) mutant brains. Despite the lack of AR-JP-like neuropathology in qk(v) mice, this mutant may constitute a readily available model for the study of the cellular function of Parkin. This is the first report of a gene distinct from qkI affected by the qk(v) deletion. The discovery of the multigenic nature of this classical mouse mutation calls for the re-evaluation of its phenotypic characterization.

Anchoring 9,371 maize expressed sequence tagged unigenes to the bacterial artificial chromosome contig map by two-dimensional overgo hybridization

Our goal is to construct a robust physical map for maize (Zea mays) comprehensively integrated with the genetic map. We have used a two-dimensional 24 x 24 overgo pooling strategy to anchor maize expressed sequence tagged (EST) unigenes to 165,888 bacterial artificial chromosomes (BACs) on high-density filters. A set of 70,716 public maize ESTs seeded derivation of 10,723 EST unigene assemblies. From these assemblies, 10,642 overgo sequences of 40 bp were applied as hybridization probes. BAC addresses were obtained for 9,371 overgo probes, representing an 88% success rate. More than 96% of the successful overgo probes identified two or more BACs, while 5% identified more than 50 BACs. The majority of BACs identified (79%) were hybridized with one or two overgos. A small number of BACs hybridized with eight or more overgos, suggesting that these BACs must be gene rich. Approximately 5,670 overgos identified BACs assembled within one contig, indicating that these probes are highly locus specific. A total of 1,795 megabases (Mb; 87%) of the total 2,050 Mb in BAC contigs were associated with one or more overgos, which are serving as sequence-tagged sites for single nucleotide polymorphism development. Overgo density ranged from less than one overgo per megabase to greater than 20 overgos per megabase. The majority of contigs (52%) hit by overgos contained three to nine overgos per megabase. Analysis of approximately 1,022 Mb of genetically anchored BAC contigs indicates that 9,003 of the total 13,900 overgo-contig sites are genetically anchored. Our results indicate overgos are a powerful approach for generating gene-specific hybridization probes that are facilitating the assembly of an integrated genetic and physical map for maize.

A comparative map of bovine chromosome 19 based on a combination of mapping on a bacterial artificial chromosome scaffold map, a whole genome radiation hybrid panel and the human draft sequence

We have constructed a medium density physical map of bovine chromosome 19 using a combination of mapping loci on both a bovine bacterial artificial chromosome (BAC) scaffold map and a whole genome radiation hybrid (WGRH) panel. The resulting map contains 70 loci spanning the length of bovine chromosome 19. Three contiguous groups of BACs were identified on the basis of multiple loci mapping to individual BAC clones. Bovine chromosome 19 was found in this study to be comprised almost entirely from regions of human chromosome 17, with a small region putatively assigned to human chromosome 10. Fourteen breakpoints between the bovine and human chromosomes were detected, with a possibility of five more based on ordering of the WGRH map.

An EST survey of the sugarcane transcriptome

Its large genome and high polyploidy makes sugarcane (Saccharum spp.) a singularly challenging crop to study and improve using genetic approaches. To provide large numbers of functionally characterized candidate genes that might be tested for direct association (rather than distant linkage) with economically important traits, we sequenced the 5' ends of 9,216 clones from three cDNA libraries (apex, leaf and mature internode), representing 3,401 non-redundant sequences. About 57% of these sequences could be assigned a tentative function based on statistically significant similarity to previously characterized proteins or DNA sequences. Another 28% corresponded to previously identified, but uncharacterized, sequences. Some of the remaining unidentified sequences were predicted to be genes which could potentially be new to plants or unique to sugarcane. Comparisons of the sugarcane ESTs to a large sorghum EST database revealed similar compositions of expressed genes between some different tissues. Comparison to a detailed Arabidopsis protein database showed some highly conserved sequences, which might be useful DNA markers for pan-angiosperm comparative mapping. These EST sequences provide a foundation for many new studies to accelerate isolation of agronomically important genes from the cumbersome sugarcane genome.

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.

Methods for high throughput validation of amplified fragment pools of BAC DNA for constructing high resolution CGH arrays

BACKGROUND

The recent development of array based comparative genomic hybridization (CGH) technology provides improved resolution for detection of genomic DNA copy number alterations. In array CGH, generating spotting solution is a multi-step process where bacterial artificial chromosome (BAC) clones are converted to replenishable PCR amplified fragments pools (AFP) for use as spotting solution in a microarray format on glass substrate. With completion of the human and mouse genome sequencing, large BAC clone sets providing complete genome coverage are available for construction of whole genome BAC arrays. Currently, Southern hybridization, fluorescent in-situ hybridization (FISH), and BAC end sequencing methods are commonly used to identify the initial BAC clone but not the end product used for spotting arrays. The AFP sequencing technique described in this study is a novel method designed to verify the identity of array spotting solution in a high throughput manner.

RESULTS

We show here that Southern hybridization, FISH, and AFP sequencing can be used to verify the identity of final spotting solutions using less than 10% of the AFP product. Single pass AFP sequencing identified over half of the 960 AFPs analyzed. Moreover, using two vector primers approximately 90% of the AFP spotting solutions can be identified.

CONCLUSIONS

In this feasibility study we demonstrate that current methods for identifying initial BAC clones can be adapted to verify the identity of AFP spotting solutions used in printing arrays. Of these methods, AFP sequencing proves to be the most efficient for large scale identification of spotting solution in a high throughput manner.

A High-Density Genetic Recombination Map of Sequence-Tagged Sites for Sorghum, as a Framework for Comparative Structural and Evolutionary Genomics of Tropical Grains and Grasses

We report a genetic recombination map for Sorghum of 2512 loci spaced at average 0.4 cM (∼300 kb) intervals based on 2050 RFLP probes, including 865 heterologous probes that foster comparative genomics of Saccharum (sugarcane), Zea (maize), Oryza (rice), Pennisetum (millet, buffelgrass), the Triticeae (wheat, barley, oat, rye), and Arabidopsis. Mapped loci identify 61.5% of the recombination events in this progeny set and reveal strong positive crossover interference acting across intervals of ≤50 cM. Significant variations in DNA marker density are related to possible centromeric regions and to probable chromosome structural rearrangements between Sorghum bicolor and S. propinquum, but not to variation in levels of intraspecific allelic richness. While cDNA and genomic clones are similarly distributed across the genome, SSR-containing clones show different abundance patterns. Rapidly evolving hypomethylated DNA may contribute to intraspecific genomic differentiation. Nonrandom distribution patterns of multiple loci detected by 357 probes suggest ancient chromosomal duplication followed by extensive rearrangement and gene loss. Exemplifying the value of these data for comparative genomics, we support and extend prior findings regarding maize-sorghum synteny—in particular, 45% of comparative loci fall outside the inferred colinear/syntenic regions, suggesting that many small rearrangements have occurred since maize-sorghum divergence. These genetically anchored sequence-tagged sites will foster many structural, functional and evolutionary genomic studies in major food, feed, and biomass crops.

Characterization of the region containing the jcpk PKD gene on mouse Chromosome 10

The jcpk gene on mouse Chromosome 10 causes a severe, early onset form of polycystic kidney disease (PKD) when inherited in an autosomal recessive manner. In order to positionally clone this gene, high resolution genetic and radiation hybrid maps were generated along with a detailed physical map of the approximately 500-kb region containing the jcpk gene. Additionally, sixty-nine kidney-specific ESTs were evaluated as candidates for jcpk and subsequently localized throughout the mouse genome by radiation hybrid mapping analysis. Previous studies indicating non-complementation of the jcpk mutation and 67Gso, a new PKD translocation mutant had suggested that 67Gso represents a new allele of jcpk. Fluorescence in situ hybridization (FISH) analysis using key bacterial artificial chromosome clones from the jcpk critical region, refined the 67Gso breakpoint and provided support for the allelism of jcpk and 67Gso.

The use of DNA polymorphisms in genetic mapping

The introduction of molecular markers has revolutionized genetics. The range of polymorphisms that are available is increasing and the advent of large-scale cDNA and genomic sequencing is a source of an ever-increasing set of available markers. The ease with which any particular marker type can be applied to an experimental system depends, to some extent, on the amount of genomic information available for that system. However, comparative genomics is enabling a wider range of marker technology to be applied to relatively information-poor systems. The types of markers that are available include restriction fragment length polymorphisms, amplified fragment length polymorphisms, ransom amplified polymorphic DNAs, simple sequence repeats, single nucleotide polymorphisms and small insertions/deletions. The types of questions that can be addressed with these molecular markers include the generation of genetic and physical maps for the identification of interesting loci, the development of marker-based gene tags, map-based cloning of agronomically important genes, synteny mapping, marker-assisted selection and quantitative trait analysis. The continued development of technology including new high throughput methods, for example those being applied to single nucleotide polymorphisms, will change the ease with which current questions can be answered as well as enable new analyses that are presently impossible to undertake.

Exclusion of the lim homeodomain gene LHX4 as a candidate gene for pituitary dwarfism in German shepherd dogs

Pituitary dwarfism in the German shepherd dog is an autosomal recessive inherited abnormality. We tested the hypothesis that a variant of the LIM homeodomain gene LHX4 is responsible for the dwarfism phenotype. To this end, we isolated Bacterial Artificial Chromosome clones for the canine LHX4 gene. Southern blotting experiments showed that the LHX4 gene is a single copy gene in the canine genome. A complex CA-repeat was isolated from the BAC clones and was found to be polymorphic in German shepherd dogs. Genotyping 5 litters in which the dwarfism was segregating showed disconcordance between the inheritance of the dwarfism phenotype and the DNA marker. It is concluded that the LHX4 gene does not play a primary role in the pituitary dwarfism in the German shepherd dogs.

Structure and evolution of the Smith-Magenis syndrome repeat gene clusters, SMS-REPs

An approximately 4-Mb genomic segment on chromosome 17p11.2, commonly deleted in patients with the Smith-Magenis syndrome (SMS) and duplicated in patients with dup(17)(p11.2p11.2) syndrome, is flanked by large, complex low-copy repeats (LCRs), termed proximal and distal SMS-REP. A third copy, the middle SMS-REP, is located between them. SMS-REPs are believed to mediate nonallelic homologous recombination, resulting in both SMS deletions and reciprocal duplications. To delineate the genomic structure and evolutionary origin of SMS-REPs, we constructed a bacterial artificial chromosome/P1 artificial chromosome contig spanning the entire SMS region, including the SMS-REPs, determined its genomic sequence, and used fluorescence in situ hybridization to study the evolution of SMS-REP in several primate species. Our analysis shows that both the proximal SMS-REP (approximately 256 kb) and the distal copy (approximately 176 kb) are located in the same orientation and derived from a progenitor copy, whereas the middle SMS-REP (approximately 241 kb) is inverted and appears to have been derived from the proximal copy. The SMS-REP LCRs are highly homologous (>98%) and contain at least 14 genes/pseudogenes each. SMS-REPs are not present in mice and were duplicated after the divergence of New World monkeys from pre-monkeys approximately 40-65 million years ago. Our findings potentially explain why the vast majority of SMS deletions and dup(17)(p11.2p11.2) occur at proximal and distal SMS-REPs and further support previous observations that higher-order genomic architecture involving LCRs arose recently during primate speciation and may predispose the human genome to both meiotic and mitotic rearrangements.

Genes in a refined Smith-Magenis syndrome critical deletion interval on chromosome 17p11.2 and the syntenic region of the mouse

Smith-Magenis syndrome (SMS) is a multiple congenital anomaly/mental retardation syndrome associated with behavioral abnormalities and sleep disturbance. Most patients have the same approximately 4 Mb interstitial genomic deletion within chromosome 17p11.2. To investigate the molecular bases of the SMS phenotype, we constructed BAC/PAC contigs covering the SMS common deletion interval and its syntenic region on mouse chromosome 11. Comparative genome analysis reveals the absence of all three approximately 200-kb SMS-REP low-copy repeats in the mouse and indicates that the evolution of SMS-REPs was accompanied by transposition of adjacent genes. Physical and genetic map comparisons in humans reveal reduced recombination in both sexes. Moreover, by examining the deleted regions in SMS patients with unusual-sized deletions, we refined the minimal Smith-Magenis critical region (SMCR) to an approximately 1.1-Mb genomic interval that is syntenic to an approxiamtely 1.0-Mb region in the mouse. Genes within the SMCR and its mouse syntenic region were identified by homology searches and by gene prediction programs, and their gene structures and expression profiles were characterized. In addition to 12 genes previously mapped, we identified 8 new genes and 10 predicted genes in the SMCR. In the mouse syntenic region of the human SMCR, 16 genes and 6 predicted genes were identified. The SMCR is highly conserved between humans and mice, including 19 genes with the same gene order and orientation. Our findings will facilitate both the identification of gene(s) responsible for the SMS phenotype and the engineering of an SMS mouse model.

Genomics and plant breeding

Much of our most basic understanding of genetics has its roots in plant genetics and crop breeding. The study of plants has led to important insights into highly conserved biological process and a wealth of knowledge about development. Agriculture is now well positioned to take its share benefit from genomics. The primary sequences of most plant genes will be determined over the next few years. Informatics and functional genomics will help identify those genes that can be best utilized to crop production and quality through genetic engineering and plant breeding. Recent developments in plant genomics are reviewed.

Genome-wide detection of chromosomal imbalances in tumors using BAC microarrays

Chromosomal imbalances such as deletions and amplifications are common rearrangements in most tumors. Specific rearrangements are consistently associated with specific tumor types or stages, implicating the role of the genes in a region of chromosomal imbalance in tumor initiation and progression. The development of comparative genomic hybridization (CGH) has obviated the need to obtain metaphase spreads from tumors, so that the chromosomal imbalances in many solid tumors may be revealed using an extracted genomic DNA sample. However, the resolution of the cytogenetic method remains and the extreme technical difficulty of CGH has restricted its use. Conceptually, DNA microarray-based CGH is an obvious solution to all of the limitations of conventional CGH. Although arrays have been used for CGH studies, their success has been limited by poor specific signal-to-noise ratios. Here we demonstrate a microarray-based CGH method that allows reliable detection of chromosomal deletions and amplifications with high resolution. Our microarray system is fundamentally different from most current microarray technologies in that activated DNA is printed on natural glass surfaces while other systems almost exclusively focus on activating the surfaces, a strategy that invariably introduces hybridization backgrounds. The concept of using pre-modification may be generally applied for making arrays of other biological materials, as modifying the substrates will be more controllable in solution than on surfaces.

An SSLP marker-anchored BAC framework map of the mouse genome

We have constructed a BAC framework map of the mouse genome consisting of 2,808 PCR-confirmed BAC clusters, using a previously described method. Fingerprints of BACs from selected clusters confirm the accuracy of the map. Combined with BAC fingerprint data, the framework map covers 37% of the mouse genome.

Five-color-based high-information-content fingerprinting of bacterial artificial chromosome clones using type IIS restriction endonucleases

We have developed a high-information-content fingerprinting (HICF) system for bacterial artificial chromosome (BAC) clones using a Type IIS restriction endonuclease, HgaI, paired with a Type II restriction endonuclease, RsaI. In the method described, unknown five-base overhangs generated with HgaI are partially or fully sequenced by modified fluorescent dideoxy terminators. Using an in-lane size standard labeled with a fifth dye, fragments are characterized by both the size and the sequence of its terminal one to five bases. The enhanced information content associated with this approach significantly increases the accuracy and efficiency of detecting shared fragments among BAC clones. We have compared data obtained from this method to predicted HICF patterns of 10 fully sequenced BACs. We have further applied HICF to 555 BAC clones to assemble contigs spanning 16p11.2 to 16p13.1 of human chromosome 16.

Comparative feline genomics: a BAC/PAC contig map of the major histocompatibility complex class II region

The genome organization of the human major histocompatibility complex (MHC) will be best understood in a comparative evolutionary context. We describe here the construction of a physical map for the feline MHC. A large-insert domestic cat genomic DNA library was developed using a P1 artificial chromosome (PAC) with a genomic representation of 2.5x and an average insert size of 80 kb. A sequence-ready 660-kb bacterial artificial chromosome/PAC contig map of the domestic cat MHC class II region was constructed with a gene order similar to, but distinct from, that of human and mice: DPB/DPA, Ring3, DMB, TAP1, DOB, DRB2, DRA3, DRB1, DRA2, and DRA1. Fluorescence in situ hybridization analyses of selected class II PAC clones confirmed that the class II region lies in the pericentromeric region of cat chromosome B2. However, apparently unlike the human and mouse MHCs, the domestic cat DRA and DRB genes have undergone multiple duplications and the DQ region has been deleted.

BAC contig from a 3-cM region of mouse chromosome 11 surrounding Brca1

Even with the completion of a draft version of the human genome sequence only a fraction of the genes identified from this sequence have known functions. Chromosomal engineering in mouse cells, in concert with gene replacement assays to prove the functional significance of a given genomic region or gene, represents a rapid and productive means for understanding the role of a given set of genes. Both techniques rely heavily on detailed maps of chromosomal regions, initially to understand the scope of the regions being modified and finally to provide the cloned resources necessary to allow both finished sequencing and large insert complementation. This report describes the creation of a BAC clone contig on mouse chromosome 11 in a region showing conservation of synteny with sequences on human chromosome 17. We have created a detailed map of an approximately 3-cM region containing at least 33 genes through the use of multiple BAC mapping strategies, including chromosome walking and multiplex oligonucleotide hybridization and gap filling. The region described is one of the targets of a large effort to create a series of mice with regional deletions on mouse chromosome 11 (33-80 cM) that can subsequently be subjected to further mutagenesis.