A complete BAC-based physical map of the Arabidopsis thaliana genome

The art of painting chromosome loops

How to get a metre of DNA into a tiny space while preserving its functional characteristics? This question seems easy to pose, but the answer is far from being trivial. Facing this riddle, salvation came from technical improvements in microscopy and in situ hybridisation techniques applied to cytogenetics. Here, we would like to look into the past at one of these pure cytogenetics articles that makes a breakthrough in addressing this question in plant science. Our choice fell on the work published two decades ago by Fransz et al. (2002). Besides the elegant manner in which DNA probes were organised to bring into light the out-looping arrangement of interphase chromosomes in Arabidopsis thaliana nuclei, this article perfectly illustrates that painting is not reserved to the fine art. As for whether emotional expression prioritised by artists can sometimes hide behind scientific empirical evidence, there is only a small step to make to the general case.

Sequencing of the Arabidopsis NOR2 reveals its distinct organization and tissue-specific rRNA ribosomal variants

Despite vast differences between organisms, some characteristics of their genomes are conserved, such as the nucleolus organizing region (NOR). The NOR is constituted of multiple, highly repetitive rDNA genes, encoding the catalytic ribosomal core RNAs which are transcribed from 45S rDNA units. Their precise sequence information and organization remain uncharacterized. Here, using a combination of long- and short-read sequencing technologies we assemble contigs of the Arabidopsis NOR2 rDNA domain. We identify several expressed rRNA gene variants which are integrated into translating ribosomes in a tissue-specific manner. These findings support the concept of tissue specific ribosome subpopulations that differ in their rRNA composition and provide insights into the higher order organization of NOR2.

The nucleolus organizing region (NOR) consists of multiple, highly repetitive rDNA genes. Here Sims et al. use both long- and short-read sequencing to determine the organization and sequence of Arabidopsis NOR2 rDNA and show that different rRNA gene variants are integrated into translating ribosomes in a tissue-specific manner.

Working with Bacterial Artificial Chromosomes (BACs) and Other High-Capacity Vectors

Genetic targeting of specific cell types is fundamentally important for modern molecular-genetic studies. The development of simple methods to engineer high-capacity vectors-in particular, bacterial artificial chromosomes (BACs)-for the preparation of transgenic lines that accurately express a gene of interest has resulted in commonplace usage of transgenic techniques in a wide variety of experimental systems. Here we provide a brief description of each of the four major types of large-capacity vectors, with a focus on the use of BAC vectors.

De Novo Sequencing and Hybrid Assembly of the Biofuel Crop Jatropha curcas L.: Identification of Quantitative Trait Loci for Geminivirus Resistance

Jatropha curcas is an important perennial, drought tolerant plant that has been identified as a potential biodiesel crop. We report here the hybrid de novo genome assembly of J. curcas generated using Illumina and PacBio sequencing technologies, and identification of quantitative loci for Jatropha Mosaic Virus (JMV) resistance. In this study, we generated scaffolds of 265.7 Mbp in length, which correspond to 84.8% of the gene space, using Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Additionally, 96.4% of predicted protein-coding genes were captured in RNA sequencing data, which reconfirms the accuracy of the assembled genome. The genome was utilized to identify 12,103 dinucleotide simple sequence repeat (SSR) markers, which were exploited in genetic diversity analysis to identify genetically distinct lines. A total of 207 polymorphic SSR markers were employed to construct a genetic linkage map for JMV resistance, using an interspecific F₂ mapping population involving susceptible J. curcas and resistant Jatropha integerrima as parents. Quantitative trait locus (QTL) analysis led to the identification of three minor QTLs for JMV resistance, and the same has been validated in an alternate F₂ mapping population. These validated QTLs were utilized in marker-assisted breeding for JMV resistance. Comparative genomics of oil-producing genes across selected oil producing species revealed 27 conserved genes and 2986 orthologous protein clusters in Jatropha. This reference genome assembly gives an insight into the understanding of the complex genetic structure of Jatropha, and serves as source for the development of agronomically improved virus-resistant and oil-producing lines.

Bioinformatics resources for deciphering the biogenesis and action pathways of plant small RNAs

The next-generation sequencing (NGS) technology has revolutionized our previous understanding of the plant genomes, relying on its innate advantages, such as high throughput and deep sequencing depth. In addition to the protein-coding gene loci, massive transcription signals have been detected within intergenic or intragenic regions. Most of these signals belong to non-coding ones, considering their weak protein-coding potential. Generally, these transcripts could be divided into long non-coding RNAs and small non-coding RNAs (sRNAs) based on their sequence length. In addition to the well-known microRNAs (miRNAs), many plant endogenous sRNAs were collectively referred to as small interfering RNAs. However, an increasing number of unclassified sRNA species are being discovered by NGS. The high heterogeneity of these novel sRNAs greatly hampered the mechanistic studies, especially on the clear description of their biogenesis and action pathways. Fortunately, public databases, bioinformatics softwares and NGS datasets are increasingly available for plant sRNA research. Here, by summarizing these valuable resources, we proposed a general workflow to decipher the RDR (RNA-dependent RNA polymerase)- and DCL (Dicer-like)-dependent biogenesis pathways, and the Argonaute-mediated action modes (such as target cleavages and chromatin modifications) for specific sRNA species in plants. Taken together, we hope that by summarizing a list of the public resources, this work will facilitate the plant biologists to perform classification and functional characterization of the interesting sRNA species.

Genome Mapping in Plant Comparative Genomics

Genome mapping produces fingerprints of DNA sequences to construct a physical map of the whole genome. It provides contiguous, long-range information that complements and, in some cases, replaces sequencing data. Recent advances in genome-mapping technology will better allow researchers to detect large (>1kbp) structural variations between plant genomes. Some molecular and informatics complications need to be overcome for this novel technology to achieve its full utility. This technology will be useful for understanding phenotype responses due to DNA rearrangements and will yield insights into genome evolution, particularly in polyploids. In this review, we outline recent advances in genome-mapping technology, including the processes required for data collection and analysis, and applications in plant comparative genomics.

De novo SNP discovery and genetic linkage mapping in poplar using restriction site associated DNA and whole-genome sequencing technologies

BACKGROUND

Restriction site associated DNA sequencing (RAD-seq), a next-generation sequencing technology, has greatly facilitated genetic linkage mapping studies in outbred species. RAD-seq is capable of discovering thousands of genetic markers for linkage mapping across many individuals, and can be applied in species with or without a reference genome. Although several analytical tools are available for RAD-seq data, alternative strategies are necessary for improving the marker quality and hence the genetic mapping accuracy.

RESULTS

We demonstrate a strategy for constructing dense genetic linkage maps in hybrid forest trees by combining RAD-seq and whole-genome sequencing technologies. We performed RAD-seq of 150 progeny and whole-genome sequencing of the two parents in an F1 hybrid population of Populus deltoides × P. simonii. Two rough references were assembled from the whole-genome sequencing reads of the two parents separately. Based on the parental reference sequences, 3442 high-quality single nucleotide polymorphisms (SNPs) were identified that segregate in the ratio of 1:1. The maternal linkage map of P. deltoides was constructed with 2012 SNPs, containing 19 linkage groups and spanning 4067.16 cM of the genome with an average distance of 2.04 cM between adjacent markers, while the male map of P. simonii consisted of 1430 SNPs and the same number of linkage groups with a total length of 4356.04 cM and an average interval distance of 3.09 cM. Collinearity between the parental linkage maps and the reference genome of P. trichocarpa was also investigated. Compared with the result on the basis of the existing reference genome, our strategy identified more high-quality SNPs and generated parental linkage groups that nicely match the karyotype of Populus.

CONCLUSIONS

The strategy of simultaneously using RAD and whole-genome sequencing technologies can be applied to constructing high-density genetic maps in forest trees regardless of whether a reference genome exists. The two parental linkage maps constructed here provide more accurate genetic resources for unraveling quantitative trait loci and accelerating molecular breeding programs, as well as for comparative genomics in Populus.

The Arabidopsis TAC Position Viewer: a high-resolution map of transformation-competent artificial chromosome (TAC) clones aligned with the Arabidopsis thaliana Columbia-0 genome

We present a high-resolution map of genomic transformation-competent artificial chromosome (TAC) clones extending over all Arabidopsis thaliana (Arabidopsis) chromosomes. The Arabidopsis genomic TAC clones have been valuable genetic tools. Previously, we constructed an Arabidopsis genomic TAC library consisting of more than 10,000 TAC clones harboring large genomic DNA fragments extending over the whole Arabidopsis genome. Here, we determined 13,577 end sequences from 6987 Arabidopsis TAC clones and mapped 5937 TAC clones to precise locations, covering approximately 90% of the Arabidopsis chromosomes. We present the large-scale data set of TAC clones with high-resolution mapping information as a Java application tool, the Arabidopsis TAC Position Viewer, which provides ready-to-go transformable genomic DNA clones corresponding to certain loci on Arabidopsis chromosomes. The TAC clone resources will accelerate genomic DNA cloning, positional walking, complementation of mutants and DNA transformation for heterologous gene expression.

Generation of knockout alleles by RFLP based BAC targeting of polymorphic embryonic stem cells

The isolation of germ line competent mouse Embryonic Stem (ES) cells and the ability to modify the genome by homologous recombination has revolutionized life science research. Since its initial discovery, several approaches have been introduced to increase the efficiency of homologous recombination, including the use of isogenic DNA for the generation of targeting constructs, and the use of Bacterial Artificial Chromosomes (BACs). BACs have the advantage of combining long stretches of homologous DNA, thereby increasing targeting efficiencies, with the possibilities delivered by BAC recombineering approaches, which provide the researcher with almost unlimited possibilities to efficiently edit the genome in a controlled fashion. Despite these advantages of BAC targeting approaches, a widespread use has been hampered, mainly because of the difficulties in identifying BAC-targeted knockout alleles by conventional methods like Southern Blotting. Recently, we introduced a novel BAC targeting strategy, in which Restriction Fragment Length Polymorphisms (RFLPs) are targeted in polymorphic mouse ES cells, enabling an efficient and easy PCR-based readout to identify properly targeted alleles. Here we provide a detailed protocol for the generation of targeting constructs, targeting of ES cells, and convenient PCR-based analysis of targeted clones, which enable the user to generate knockout ES cells of almost every gene in the mouse genome within a 2-month period.

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.

Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.)

Three linkage maps of flax (Linum usitatissimum L.) were constructed from populations CDC Bethune/Macbeth, E1747/Viking and SP2047/UGG5-5 containing between 385 and 469 mapped markers each. The first consensus map of flax was constructed incorporating 770 markers based on 371 shared markers including 114 that were shared by all three populations and 257 shared between any two populations. The 15 linkage group map corresponds to the haploid number of chromosomes of this species. The marker order of the consensus map was largely collinear in all three individual maps but a few local inversions and marker rearrangements spanning short intervals were observed. Segregation distortion was present in all linkage groups which contained 1-52 markers displaying non-Mendelian segregation. The total length of the consensus genetic map is 1,551 cM with a mean marker density of 2.0 cM. A total of 670 markers were anchored to 204 of the 416 fingerprinted contigs of the physical map corresponding to ~274 Mb or 74 % of the estimated flax genome size of 370 Mb. This high resolution consensus map will be a resource for comparative genomics, genome organization, evolution studies and anchoring of the whole genome shotgun sequence.

[Progress on whole genome sequencing in woody plants]

In recent years, the number of sequencing data of plant whole genome have been increasing rapidly and the whole genome sequencing has been also performed widely in woody plants. However, there are a set of obstacles in investigating the whole genome sequencing in woody plants, which include larger genome, complex genome structure, limitations of assembly, annotation, functional analysis, and restriction of the funds for scientific research. Therefore, to promote the efficiency of the whole genome sequencing in woody plants, the development and defect of this field should be analyzed. The three-generation sequencing technologies (i.e., Sanger sequencing, synthesis sequencing, and single molecule sequencing) were compared in our studies. The progress mainly focused on the whole genome sequencing in four woody plants (Populus, Grapevine, Papaya, and Apple), and the application of sequencing results also was analyzed. The future of whole genome sequencing research in woody plants, consisting of material selection, establishment of genetic map and physical map, selection of sequencing technology, bioinformatic analysis, and application of sequencing results, was discussed.

Genome physical mapping of polyploids: a BIBAC physical map of cultivated tetraploid cotton, Gossypium hirsutum L

Polyploids account for approximately 70% of flowering plants, including many field, horticulture and forage crops. Cottons are a world-leading fiber and important oilseed crop, and a model species for study of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. This study has addressed the concerns of physical mapping of polyploids with BACs and/or BIBACs by constructing a physical map of the tetraploid cotton, Gossypium hirsutum L. The physical map consists of 3,450 BIBAC contigs with an N50 contig size of 863 kb, collectively spanning 2,244 Mb. We sorted the map contigs according to their origin of subgenome, showing that we assembled physical maps for the A- and D-subgenomes of the tetraploid cotton, separately. We also identified the BIBACs in the map minimal tilling path, which consists of 15,277 clones. Moreover, we have marked the physical map with nearly 10,000 BIBAC ends (BESs), making one BES in approximately 250 kb. This physical map provides a line of evidence and a strategy for physical mapping of polyploids, and a platform for advanced research of the tetraploid cotton genome, particularly fine mapping and cloning the cotton agronomic genes and QTLs, and sequencing and assembling the cotton genome using the modern next-generation sequencing technology.

Integration of the Draft Sequence and Physical Map as a Framework for Genomic Research in Soybean (Glycine max (L.) Merr.) and Wild Soybean (Glycine soja Sieb. and Zucc.)

Soybean is a model for the legume research community because of its importance as a crop, densely populated genetic maps, and the availability of a genome sequence. Even though a whole-genome shotgun sequence and bacterial artificial chromosome (BAC) libraries are available, a high-resolution, chromosome-based physical map linked to the sequence assemblies is still needed for whole-genome alignments and to facilitate map-based gene cloning. Three independent G. max BAC libraries combined with genetic and gene-based markers were used to construct a minimum tiling path (MTP) of BAC clones. A total of 107,214 clones were assembled into 1355 FPC (FingerPrinted Contigs) contigs, incorporating 4628 markers and aligned to the G. max reference genome sequence using BAC end-sequence information. Four different MTPs were made for G. max that covered from 92.6% to 95.0% of the soybean draft genome sequence (gmax1.01). Because our purpose was to pick the most reliable and complete MTP, and not the MTP with the minimal number of clones, the FPC map and draft sequence were integrated and clones with unpaired BES were added to build a high-quality physical map with the fewest gaps possible (http://soybase.org). A physical map was also constructed for the undomesticated ancestor (G. soja) of soybean to explore genome variation between G. max and G. soja. 66,028 G. soja clones were assembled into 1053 FPC contigs covering approximately 547 Mbp of the G. max genome sequence. These physical maps for G. max and its undomesticated ancestor, G. soja, will serve as a framework for ordering sequence fragments, comparative genomics, cloning genes, and evolutionary analyses of legume genomes.

A comparative physical map reveals the pattern of chromosomal evolution between the turkey (Meleagris gallopavo) and chicken (Gallus gallus) genomes

Background

A robust bacterial artificial chromosome (BAC)-based physical map is essential for many aspects of genomics research, including an understanding of chromosome evolution, high-resolution genome mapping, marker-assisted breeding, positional cloning of genes, and quantitative trait analysis. To facilitate turkey genetics research and better understand avian genome evolution, a BAC-based integrated physical, genetic, and comparative map was developed for this important agricultural species.

Results

The turkey genome physical map was constructed based on 74,013 BAC fingerprints (11.9 × coverage) from two independent libraries, and it was integrated with the turkey genetic map and chicken genome sequence using over 41,400 BAC assignments identified by 3,499 overgo hybridization probes along with > 43,000 BAC end sequences. The physical-comparative map consists of 74 BAC contigs, with an average contig size of 13.6 Mb. All but four of the turkey chromosomes were spanned on this map by three or fewer contigs, with 14 chromosomes spanned by a single contig and nine chromosomes spanned by two contigs. This map predicts 20 to 27 major rearrangements distinguishing turkey and chicken chromosomes, despite up to 40 million years of separate evolution between the two species. These data elucidate the chromosomal evolutionary pattern within the Phasianidae that led to the modern turkey and chicken karyotypes. The predominant rearrangement mode involves intra-chromosomal inversions, and there is a clear bias for these to result in centromere locations at or near telomeres in turkey chromosomes, in comparison to interstitial centromeres in the orthologous chicken chromosomes.

Conclusion

The BAC-based turkey-chicken comparative map provides novel insights into the evolution of avian genomes, a framework for assembly of turkey whole genome shotgun sequencing data, and tools for enhanced genetic improvement of these important agricultural and model species.

Sequencing crop genomes: approaches and applications

Many challenges face plant scientists, in particular those working on crop production, such as a projected increase in population, decrease in water and arable land, changes in weather patterns and predictability. Advances in genome sequencing and resequencing can and should play a role in our response to meeting these challenges. However, several barriers prevent rapid and effective deployment of these tools to a wide variety of crops. Because of the complexity of crop genomes, de novo sequencing with next-generation sequencing technologies is a process fraught with difficulties that then create roadblocks to the utilization of these genome sequences for crop improvement. Collecting rapid and accurate phenotypes in crop plants is a hindrance to integrating genomics with crop improvement, and advances in informatics are needed to put these tools in the hands of the scientists on the ground.

Construction of papaya male and female BAC libraries and application in physical mapping of the sex chromosomes

Papaya is a major fruit crop in the tropics and has recently evolved sex chromosomes. Towards sequencing the papaya sex chromosomes, two bacterial artificial chromosome (BAC) libraries were constructed from papaya male and female genomic DNA. The female BAC library was constructed using restriction enzyme BstY I and consists of 36,864 clones with an average insert size of 104 kb, providing 10.3x genome equivalents. The male BAC library was constructed using restriction enzyme EcoR I and consists of 55,296 clones with an average insert size of 101 kb, providing 15.0x genome equivalents. The male BAC library was used in constructing the physical map of the male-specific region of the male Y chromosome (MSY) and in filling gaps and extending the physical map of the hermaphrodite-specific region of the Y^h chromosome (HSY) and the X chromosome physical map. The female BAC library was used to extend the X physical map gap. The MSY, HSY, and X physical maps offer a unique opportunity to study chromosomal rearrangements, Y chromosome degeneration, and dosage compensation of the papaya nascent sex chromosomes.

Physical mapping and BAC-end sequence analysis provide initial insights into the flax (Linum usitatissimum L.) genome

BACKGROUND

Flax (Linum usitatissimum L.) is an important source of oil rich in omega-3 fatty acids, which have proven health benefits and utility as an industrial raw material. Flax seeds also contain lignans which are associated with reducing the risk of certain types of cancer. Its bast fibres have broad industrial applications. However, genomic tools needed for molecular breeding were non existent. Hence a project, Total Utilization Flax GENomics (TUFGEN) was initiated. We report here the first genome-wide physical map of flax and the generation and analysis of BAC-end sequences (BES) from 43,776 clones, providing initial insights into the genome.

RESULTS

The physical map consists of 416 contigs spanning ~368 Mb, assembled from 32,025 fingerprints, representing roughly 54.5% to 99.4% of the estimated haploid genome (370-675 Mb). The N50 size of the contigs was estimated to be ~1,494 kb. The longest contig was ~5,562 kb comprising 437 clones. There were 96 contigs containing more than 100 clones. Approximately 54.6 Mb representing 8-14.8% of the genome was obtained from 80,337 BES. Annotation revealed that a large part of the genome consists of ribosomal DNA (~13.8%), followed by known transposable elements at 6.1%. Furthermore, ~7.4% of sequence was identified to harbour novel repeat elements. Homology searches against flax-ESTs and NCBI-ESTs suggested that ~5.6% of the transcriptome is unique to flax. A total of 4064 putative genomic SSRs were identified and are being developed as novel markers for their use in molecular breeding.

CONCLUSION

The first genome-wide physical map of flax constructed with BAC clones provides a framework for accessing target loci with economic importance for marker development and positional cloning. Analysis of the BES has provided insights into the uniqueness of the flax genome. Compared to other plant genomes, the proportion of rDNA was found to be very high whereas the proportion of known transposable elements was low. The SSRs identified from BES will be valuable in saturating existing linkage maps and for anchoring physical and genetic maps. The physical map and paired-end reads from BAC clones will also serve as scaffolds to build and validate the whole genome shotgun assembly.

An integrated physical, genetic and cytogenetic map of Brachypodium distachyon, a model system for grass research

The pooid subfamily of grasses includes some of the most important crop, forage and turf species, such as wheat, barley and Lolium. Developing genomic resources, such as whole-genome physical maps, for analysing the large and complex genomes of these crops and for facilitating biological research in grasses is an important goal in plant biology. We describe a bacterial artificial chromosome (BAC)-based physical map of the wild pooid grass Brachypodium distachyon and integrate this with whole genome shotgun sequence (WGS) assemblies using BAC end sequences (BES). The resulting physical map contains 26 contigs spanning the 272 Mb genome. BES from the physical map were also used to integrate a genetic map. This provides an independent vaildation and confirmation of the published WGS assembly. Mapped BACs were used in Fluorescence In Situ Hybridisation (FISH) experiments to align the integrated physical map and sequence assemblies to chromosomes with high resolution. The physical, genetic and cytogenetic maps, integrated with whole genome shotgun sequence assemblies, enhance the accuracy and durability of this important genome sequence and will directly facilitate gene isolation.

Rapid bacterial artificial chromosome modification for large-scale mouse transgenesis

We report here a high-throughput method for the modification of bacterial artificial chromosomes (BACs) that uses a novel two-plasmid approach. In this protocol, a vector modified in our laboratory to hold an R6Kγ origin of replication and a marker recombination cassette is inserted into a BAC in a single recombination step. Temporal control of recombination is achieved through the use of a second plasmid, pSV1.RecA, which possesses a recombinase gene and a temperature-sensitive origin of replication. This highly efficient protocol has allowed us to successfully modify more than 2,000 BACs, from which over 1,000 BAC transgenic mice have been generated. A complete cycle from BAC choice to embryo implantation takes about 5 weeks. Marker genes introduced into the mice include EGFP and EGFP-L10a. All vectors used in this project can be obtained from us by request, and the EGFP reporter mice are available through the Mutant Mouse Regional Resource Center (NINDS/GENSAT collection). CNS anatomical expression maps of the mice are available to the public at http://www.gensat.org/.

Generation of a BAC-based physical map of the melon genome

BACKGROUND

Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has high intra-specific genetic variation, morphologic diversity and a small genome size (450 Mb), which make this species suitable for a great variety of molecular and genetic studies that can lead to the development of tools for breeding varieties of the species. A number of genetic and genomic resources have already been developed, such as several genetic maps and BAC genomic libraries. These tools are essential for the construction of a physical map, a valuable resource for map-based cloning, comparative genomics and assembly of whole genome sequencing data. However, no physical map of any Cucurbitaceae has yet been developed. A project has recently been started to sequence the complete melon genome following a whole-genome shotgun strategy, which makes use of massive sequencing data. A BAC-based melon physical map will be a useful tool to help assemble and refine the draft genome data that is being produced.

RESULTS

A melon physical map was constructed using a 5.7 x BAC library and a genetic map previously developed in our laboratories. High-information-content fingerprinting (HICF) was carried out on 23,040 BAC clones, digesting with five restriction enzymes and SNaPshot labeling, followed by contig assembly with FPC software. The physical map has 1,355 contigs and 441 singletons, with an estimated physical length of 407 Mb (0.9 x coverage of the genome) and the longest contig being 3.2 Mb. The anchoring of 845 BAC clones to 178 genetic markers (100 RFLPs, 76 SNPs and 2 SSRs) also allowed the genetic positioning of 183 physical map contigs/singletons, representing 55 Mb (12%) of the melon genome, to individual chromosomal loci. The melon FPC database is available for download at http://melonomics.upv.es/static/files/public/physical_map/.

CONCLUSIONS

Here we report the construction of the first physical map of a Cucurbitaceae species described so far. The physical map was integrated with the genetic map so that a number of physical contigs, representing 12% of the melon genome, could be anchored to known genetic positions. The data presented is already helping to improve the quality of the melon genomic sequence available as a result of a project currently being carried out in Spain, adopting a whole genome shotgun approach based on 454 sequencing data.

Physical mapping in highly heterozygous genomes: a physical contig map of the Pinot Noir grapevine cultivar

BACKGROUND

Most of the grapevine (Vitis vinifera L.) cultivars grown today are those selected centuries ago, even though grapevine is one of the most important fruit crops in the world. Grapevine has therefore not benefited from the advances in modern plant breeding nor more recently from those in molecular genetics and genomics: genes controlling important agronomic traits are practically unknown. A physical map is essential to positionally clone such genes and instrumental in a genome sequencing project.

RESULTS

We report on the first whole genome physical map of grapevine built using high information content fingerprinting of 49,104 BAC clones from the cultivar Pinot Noir. Pinot Noir, as most grape varieties, is highly heterozygous at the sequence level. This resulted in the two allelic haplotypes sometimes assembling into separate contigs that had to be accommodated in the map framework or in local expansions of contig maps. We performed computer simulations to assess the effects of increasing levels of sequence heterozygosity on BAC fingerprint assembly and showed that the experimental assembly results are in full agreement with the theoretical expectations, given the heterozygosity levels reported for grape. The map is anchored to a dense linkage map consisting of 994 markers. 436 contigs are anchored to the genetic map, covering 342 of the 475 Mb that make up the grape haploid genome.

CONCLUSIONS

We have developed a resource that makes it possible to access the grapevine genome, opening the way to a new era both in grape genetics and breeding and in wine making. The effects of heterozygosity on the assembly have been analyzed and characterized by using several complementary approaches which could be easily transferred to the study of other genomes which present the same features.

Generation of transgenic rats using YAC and BAC DNA constructs

Transgenic rats with a simple plasmid vector smaller than 20 Kb show insufficient expression and tissue specificity of the introduced transgene. Vectors derived from yeast artificial chromosome (YAC) and bacterial artificial chromosome (BAC), consisting of DNA fragments up to approximately 1 Mb (YAC) and approximately 200 Kb (BAC), respectively, and containing various endogenous regulatory sequences, were expected to work well and showed expression profiles comparable to their endogenous counterparts in transgenic animals. While attempting to make transgenic rats using YAC and BAC vectors, we faced two problems: how to prepare sufficiently concentrated intact DNA and how to reliably microinject a large DNA fragment into the fragile pronuclear ova of the rat. After solving these problems, we were able to make transgenic rats by introducing YAC/BAC gene constructs (YACs/BACs) into the pronuclear ova. And then we examined the relative transcription rates of these genes in the transgenic rats. In this chapter, we focus on this injection process.

Arabidopsis FtsZ2-1 and FtsZ2-2 are functionally redundant, but FtsZ-based plastid division is not essential for chloroplast partitioning or plant growth and development

FtsZ1 and FtsZ2 are phylogenetically distinct families of FtsZ in plants that co-localize to mid-plastid rings and facilitate division of chloroplasts. In plants, altered levels of either FtsZ1 or FtsZ2 cause dose-dependent defects in chloroplast division; thus, studies on the functional relationship between FtsZ genes require careful manipulation of FtsZ levels in vivo. To define the functional relationship between the two FtsZ2 genes in Arabidopsis thaliana, FtsZ2-1 and FtsZ2-2, we expressed FtsZ2-1 in an ftsZ2-2 null mutant, and vice versa, and determined whether the chloroplast division defects were rescued in plants expressing different total levels of FtsZ2. Full rescue was observed when either the FtsZ2-1 or FtsZ2-2 level approximated total FtsZ2 levels in wild-type (WT). Additionally, FtsZ2-2 interacts with ARC6, as shown previously for FtsZ2-1. These data indicate that FtsZ2-1 and FtsZ2-2 are functionally redundant for chloroplast division in Arabidopsis. To rigorously validate the requirement of each FtsZ family for chloroplast division, we replaced FtsZ1 with FtsZ2 in vivo, and vice versa, while maintaining the FtsZ level in the transgenic plants equal to that of the total level in WT. Chloroplast division defects were not rescued, demonstrating conclusively that FtsZ1 and FtsZ2 are non-redundant for maintenance of WT chloroplast numbers. Finally, we generated ftsZ triple null mutants and show that plants completely devoid of FtsZ protein are viable and fertile. As plastids are presumably essential organelles, these findings suggest that an FtsZ-independent mode of plastid partitioning may occur in higher plants.

A physical map of the papaya genome with integrated genetic map and genome sequence

BACKGROUND

Papaya is a major fruit crop in tropical and subtropical regions worldwide and has primitive sex chromosomes controlling sex determination in this trioecious species. The papaya genome was recently sequenced because of its agricultural importance, unique biological features, and successful application of transgenic papaya for resistance to papaya ringspot virus. As a part of the genome sequencing project, we constructed a BAC-based physical map using a high information-content fingerprinting approach to assist whole genome shotgun sequence assembly.

RESULTS

The physical map consists of 963 contigs, representing 9.4x genome equivalents, and was integrated with the genetic map and genome sequence using BAC end sequences and a sequence-tagged high-density genetic map. The estimated genome coverage of the physical map is about 95.8%, while 72.4% of the genome was aligned to the genetic map. A total of 1,181 high quality overgo (overlapping oligonucleotide) probes representing conserved sequences in Arabidopsis and genetically mapped loci in Brassica were anchored on the physical map, which provides a foundation for comparative genomics in the Brassicales. The integrated genetic and physical map aligned with the genome sequence revealed recombination hotspots as well as regions suppressed for recombination across the genome, particularly on the recently evolved sex chromosomes. Suppression of recombination spread to the adjacent region of the male specific region of the Y chromosome (MSY), and recombination rates were recovered gradually and then exceeded the genome average. Recombination hotspots were observed at about 10 Mb away on both sides of the MSY, showing 7-fold increase compared with the genome wide average, demonstrating the dynamics of recombination of the sex chromosomes.

CONCLUSION

A BAC-based physical map of papaya was constructed and integrated with the genetic map and genome sequence. The integrated map facilitated the draft genome assembly, and is a valuable resource for comparative genomics and map-based cloning of agronomically and economically important genes and for sex chromosome research.

Markers and mapping revisited: finding your gene

This paper is an update of our earlier review (Jones et al., 1997, Markers and mapping: we are all geneticists now. New Phytologist 137: 165-177), which dealt with the genetics of mapping, in terms of recombination as the basis of the procedure, and covered some of the first generation of markers, including restriction fragment length polymorphisms (RFLPs), random amplified polymorphic DNA (RAPDs), simple sequence repeats (SSRs) and quantitative trait loci (QTLs). In the intervening decade there have been numerous developments in marker science with many new systems becoming available, which are herein described: cleavage amplification polymorphism (CAP), sequence-specific amplification polymorphism (S-SAP), inter-simple sequence repeat (ISSR), sequence tagged site (STS), sequence characterized amplification region (SCAR), selective amplification of microsatellite polymorphic loci (SAMPL), single nucleotide polymorphism (SNP), expressed sequence tag (EST), sequence-related amplified polymorphism (SRAP), target region amplification polymorphism (TRAP), microarrays, diversity arrays technology (DArT), single-strand conformation polymorphism (SSCP), denaturing gradient gel electrophoresis (DGGE), temperature gradient gel electrophoresis (TGGE) and methylation-sensitive PCR. In addition there has been an explosion of knowledge and databases in the area of genomics and bioinformatics. The number of flowering plant ESTs is c. 19 million and counting, with all the opportunity that this provides for gene-hunting, while the survey of bioinformatics and computer resources points to a rapid growth point for future activities in unravelling and applying the burst of new information on plant genomes. A case study is presented on tracking down a specific gene (stay-green (SGR), a post-transcriptional senescence regulator) using the full suite of mapping tools and comparative mapping resources. We end with a brief speculation on how genome analysis may progress into the future of this highly dynamic arena of plant science.

Gene expression trends and protein features effectively complement each other in gene function prediction

MOTIVATION

Genome-scale 'omics' data constitute a potentially rich source of information about biological systems and their function. There is a plethora of tools and methods available to mine omics data. However, the diversity and complexity of different omics data types is a stumbling block for multi-data integration, hence there is a dire need for additional methods to exploit potential synergy from integrated orthogonal data. Rough Sets provide an efficient means to use complex information in classification approaches. Here, we set out to explore the possibilities of Rough Sets to incorporate diverse information sources in a functional classification of unknown genes.

RESULTS

We explored the use of Rough Sets for a novel data integration strategy where gene expression data, protein features and Gene Ontology (GO) annotations were combined to describe general and biologically relevant patterns represented by If-Then rules. The descriptive rules were used to predict the function of unknown genes in Arabidopsis thaliana and Schizosaccharomyces pombe. The If-Then rule models showed success rates of up to 0.89 (discriminative and predictive power for both modeled organisms); whereas, models built solely of one data type (protein features or gene expression data) yielded success rates varying from 0.68 to 0.78. Our models were applied to generate classifications for many unknown genes, of which a sizeable number were confirmed either by PubMed literature reports or electronically interfered annotations. Finally, we studied cell cycle protein-protein interactions derived from both tandem affinity purification experiments and in silico experiments in the BioGRID interactome database and found strong experimental evidence for the predictions generated by our models. The results show that our approach can be used to build very robust models that create synergy from integrating gene expression data and protein features.

AVAILABILITY

The Rough Set-based method is implemented in the Rosetta toolkit kernel version 1.0.1 available at: http://rosetta.lcb.uu.se/

Simple sequence repeat marker development and genetic mapping in quinoa (Chenopodium quinoa Willd.)

Quinoa is a regionally important grain crop in the Andean region of South America. Recently quinoa has gained international attention for its high nutritional value and tolerances of extreme abiotic stresses. DNA markers and linkage maps are important tools for germplasm conservation and crop improvement programmes. Here we report the development of 216 new polymorphic SSR (simple sequence repeats) markers from libraries enriched for GA, CAA and AAT repeats, as well as 6 SSR markers developed from bacterial artificial chromosome-end sequences (BES-SSRs). Heterozygosity (H) values of the SSR markers ranges from 0.12 to 0.90, with an average value of 0.57. A linkage map was constructed for a newly developed recombinant inbred lines (RIL) population using these SSR markers. Additional markers, including amplified fragment length polymorphisms (AFLPs), two 11S seed storage protein loci, and the nucleolar organizing region (NOR), were also placed on the linkage map. The linkage map presented here is the first SSR-based map in quinoa and contains 275 markers, including 200 SSR. The map consists of 38 linkage groups (LGs) covering 913 cM. Segregation distortion was observed in the mapping population for several marker loci, indicating possible chromosomal regions associated with selection or gametophytic lethality. As this map is based primarily on simple and easily-transferable SSR markers, it will be particularly valuable for research in laboratories in Andean regions of South America.

An overview of the apple genome through BAC end sequence analysis

The apple, Malus x domestica Borkh., is one of the most important fruit trees grown worldwide. A bacterial artificial chromosome (BAC)-based physical map of the apple genome has been recently constructed. Based on this physical map, a total of approximately 2,100 clones from different contigs (overlapping BAC clones) have been selected and sequenced at both ends, generating 3,744 high-quality BAC end sequences (BESs) including 1,717 BAC end pairs. Approximately 8.5% of BESs contain simple sequence repeats (SSRs), most of which are AT/TA dimer repeats. Potential transposable elements are identified in approximately 21% of BESs, and most of these elements are retrotransposons. About 11% of BESs have homology to the Arabidopsis protein database. The matched proteins cover a broad range of categories. The average GC content of the predicted coding regions of BESs is 42.4%; while, that of the whole BESs is 39%. A small number of BES pairs were mapped to neighboring chromosome regions of A. thaliana and Populus trichocarpa; whereas, no pairs are mapped to the Oryza sativa genome. The apple has a higher degree of synteny with the closely related Populus than with the distantly related Arabidopsis. BAC end sequencing can be used to anchor a small proportion of the apple genome to the Populus and possibly to the Arabidopsis genomes.

Integration of cytogenetic and genetic linkage maps unveils the physical architecture of tomato chromosome 2

We report the integration of the linkage map of tomato chromosome 2 with a high-density bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH)-based cytogenetic map. The euchromatic block of chromosome 2 resides between 13 and 142 cM and has a physical length of 48.12 microm, with 1 microm equivalent to 540 kb. BAC-FISH resolved a pair of loci that were 3.7-3.9 Mb apart and were not resolved on the linkage map. Most of the regions had crossover densities close to the mean of approximately 200 kb/cM. Relatively hot and cold spots of recombination were unevenly distributed along the chromosome. The distribution of centimorgan/micrometer values was similar to the previously reported recombination nodule distribution along the pachytene chromosome. FISH-based physical maps will play an important role in advanced genomics research for tomato, including map-based cloning of agronomically important traits and whole-genome sequencing.

Simple sequence repeat marker development and genetic mapping in quinoa (Chenopodium quinoa Willd.)

Quinoa is a regionally important grain crop in the Andean region of South America. Recently quinoa has gained international attention for its high nutritional value and tolerances of extreme abiotic stresses. DNA markers and linkage maps are important tools for germplasm conservation and crop improvement programmes. Here we report the development of 216 new polymorphic SSR (simple sequence repeats) markers from libraries enriched for GA, CAA and AAT repeats, as well as 6 SSR markers developed from bacterial artificial chromosome-end sequences (BES-SSRs). Heterozygosity (H) values of the SSR markers ranges from 0.12 to 0.90, with an average value of 0.57. A linkage map was constructed for a newly developed recombinant inbred lines (RIL) population using these SSR markers. Additional markers, including amplified fragment length polymorphisms (AFLPs), two 11S seed storage protein loci, and the nucleolar organizing region (NOR), were also placed on the linkage map. The linkage map presented here is the first SSR-based map in quinoa and contains 275 markers, including 200 SSR. The map consists of 38 linkage groups (LGs) covering 913 cM. Segregation distortion was observed in the mapping population for several marker loci, indicating possible chromosomal regions associated with selection or gametophytic lethality. As this map is based primarily on simple and easily-transferable SSR markers, it will be particularly valuable for research in laboratories in Andean regions of South America.

A physical map of the heterozygous grapevine 'Cabernet Sauvignon' allows mapping candidate genes for disease resistance

BACKGROUND

Whole-genome physical maps facilitate genome sequencing, sequence assembly, mapping of candidate genes, and the design of targeted genetic markers. An automated protocol was used to construct a Vitis vinifera 'Cabernet Sauvignon' physical map. The quality of the result was addressed with regard to the effect of high heterozygosity on the accuracy of contig assembly. Its usefulness for the genome-wide mapping of genes for disease resistance, which is an important trait for grapevine, was then assessed.

RESULTS

The physical map included 29,727 BAC clones assembled into 1,770 contigs, spanning 715,684 kbp, and corresponding to 1.5-fold the genome size. Map inflation was due to high heterozygosity, which caused either the separation of allelic BACs in two different contigs, or local mis-assembly in contigs containing BACs from the two haplotypes. Genetic markers anchored 395 contigs or 255,476 kbp to chromosomes. The fully automated assembly and anchorage procedures were validated by BAC-by-BAC blast of the end sequences against the grape genome sequence, unveiling 7.3% of chimerical contigs. The distribution across the physical map of candidate genes for non-host and host resistance, and for defence signalling pathways was then studied. NBS-LRR and RLK genes for host resistance were found in 424 contigs, 133 of them (32%) were assigned to chromosomes, on which they are mostly organised in clusters. Non-host and defence signalling genes were found in 99 contigs dispersed without a discernable pattern across the genome.

CONCLUSION

Despite some limitations that interfere with the correct assembly of heterozygous clones into contigs, the 'Cabernet Sauvignon' physical map is a useful and reliable intermediary step between a genetic map and the genome sequence. This tool was successfully exploited for a quick mapping of complex families of genes, and it strengthened previous clues of co-localisation of major NBS-LRR clusters and disease resistance loci in grapevine.

CAUT lines: a novel resource for studies of cell autonomy in Arabidopsis

Plant development is critically dependent on the interactions between clonally unrelated cell layers. The cross-talk between layers can be addressed by studies of cell autonomy. Cell autonomy is a property of genetic mosaics composed of cells of differing genotypes. Broadly, if the phenotype of a mutant tissue reflects only its genotype and is unaffected by the presence of wild-type tissue, the trait is cell-autonomous. Conversely, if the phenotype of a mutant tissue reflects that of wild-type tissue in the mosaic, the trait is non-autonomous. Here we report a novel, versatile and robust method for studies of cell autonomy in Arabidopsis. Cell autonomy (CAUT) lines consist of a collection of homozygous stocks, each containing one of 76 mapped T-DNA inserts, each of which corrects the yellow ch-42 mutant to green (CH-42) by complementation. This has the effect of translocating the colour marker to 76 new locations around the genome. X-irradiation of heterozygous CAUT line seeds results in yellow sectors, with loss of the CH-42 transgene and adjacent wild-type genes. This property can be used to remove the wild-type copy of developmental genes in appropriate heterozygotes, resulting in yellow (ch-42) sectors that are hemizygous for the trait of interest. Such sectors can provide insight into cell autonomy. Experiments using the ap1, ap3, ag and clv1 mutants show that CAUT lines are useful in the study of cell autonomy.

BAC libraries construction from the ancestral diploid genomes of the allotetraploid cultivated peanut

BACKGROUND

Cultivated peanut, Arachis hypogaea is an allotetraploid of recent origin, with an AABB genome. In common with many other polyploids, it seems that a severe genetic bottle-neck was imposed at the species origin, via hybridisation of two wild species and spontaneous chromosome duplication. Therefore, the study of the genome of peanut is hampered both by the crop's low genetic diversity and its polyploidy. In contrast to cultivated peanut, most wild Arachis species are diploid with high genetic diversity. The study of diploid Arachis genomes is therefore attractive, both to simplify the construction of genetic and physical maps, and for the isolation and characterization of wild alleles. The most probable wild ancestors of cultivated peanut are A. duranensis and A. ipaënsis with genome types AA and BB respectively.

RESULTS

We constructed and characterized two large-insert libraries in Bacterial Artificial Chromosome (BAC) vector, one for each of the diploid ancestral species. The libraries (AA and BB) are respectively c. 7.4 and c. 5.3 genome equivalents with low organelle contamination and average insert sizes of 110 and 100 kb. Both libraries were used for the isolation of clones containing genetically mapped legume anchor markers (single copy genes), and resistance gene analogues.

CONCLUSION

These diploid BAC libraries are important tools for the isolation of wild alleles conferring resistances to biotic stresses, comparisons of orthologous regions of the AA and BB genomes with each other and with other legume species, and will facilitate the construction of a physical map.

Novel bacterial artificial chromosome vector pUvBBAC for use in studies of the functional genomics of Listeria spp

Bacterial artificial chromosome (BAC) vectors are important tools for microbial genome research. We constructed a novel BAC vector, pUvBBAC, for replication in both gram-negative and gram-positive bacterial hosts. The pUvBBAC vector was used to generate a BAC library for the facultative intracellular pathogen Listeria monocytogenes EGD-e. The library had insert sizes ranging from 68 to 178 kb. We identified two recombinant BACs from the L. monocytogenes pUvBBAC library that each contained the entire virulence gene cluster (vgc) of L. monocytogenes and transferred them to a nonpathogenic Listeria innocua strain. Recombinant L. innocua strains harboring pUvBBAC+vgc1 and pUvBBAC+vgc2 produced the vgc-specific listeriolysin (LLO) and actin assembly protein ActA and represent the first reported cloning of the vgc locus in its entirety. The use of the novel broad-host-range BAC vector pUvBBAC extends the versatility of this technology and provides a powerful platform for detailed functional genomics of gram-positive bacteria as well as its use in explorative functional metagenomics.

End-sequencing and characterization of silkworm (Bombyx mori) bacterial artificial chromosome libraries

Background

We performed large-scale bacterial artificial chromosome (BAC) end-sequencing of two BAC libraries (an EcoRI- and a BamHI-digested library) and conducted an in silico analysis to characterize the obtained sequence data, to make them a useful resource for genomic research on the silkworm (Bombyx mori).

Results

More than 94000 BAC end sequences (BESs), comprising more than 55 Mbp and covering about 10.4% of the silkworm genome, were sequenced. Repeat-sequence analysis with known repeat sequences indicated that the long interspersed nuclear elements (LINEs) were abundant in BamHI BESs, whereas DNA-type elements were abundant in EcoRI BESs. Repeat-sequence analysis revealed that the abundance of LINEs might be due to a GC bias of the restriction sites and that the GC content of silkworm LINEs was higher than that of mammalian LINEs. In a BLAST-based sequence analysis of the BESs against two available whole-genome shotgun sequence data sets, more than 70% of the BESs had a BLAST hit with an identity of ≥ 99%. About 14% of EcoRI BESs and about 8% of BamHI BESs were paired-end clones with unique sequences at both ends. Cluster analysis of the BESs clarified the proportion of BESs containing protein-coding regions.

Conclusion

As a result of this characterization, the identified BESs will be a valuable resource for genomic research on Bombyx mori, for example, as a base for construction of a BAC-based physical map. The use of multiple complementary BAC libraries constructed with different restriction enzymes also makes the BESs a more valuable genomic resource. The GenBank accession numbers of the obtained end sequences are DE283657–DE378560.

The bacterial artificial chromosome (BAC) library of the narrow-leafed lupin (Lupinus angustifolius L.)

The narrow-leafed lupin possesses valuable traits for environment-friendly agriculture and for the production of unconventional agricultural products. Despite various genetic and environmental studies, the breeding of improved cultivars has been slow due to the limited knowledge of its genomic structure. Further advances in genomics require, among other things, the availability of a genomic DNA library with large inserts. We report here on the construction of the first DNA library cloned in a BAC (bacterial artificial chromosome) vector from diploid Lupinus angustifolius L. cv. Sonet. The high molecular weight DNA used for its preparation was isolated from interphase nuclei that were purified by flow cytometry. The library comprises 55,296 clones and is ordered in 144x384-well microtitre plates. With an average insert size of 100 kb, the library represents six haploid genome equivalents. Thanks to the purification of the nuclei by flow cytometry, contamination with chloroplast DNA and mitochondrial DNA was negligible. The availability of a BAC library opens avenues for the development of a physical contig map and positional gene cloning, as well as for the analysis of the plant's genome structure and evolution.

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 +/- 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa, version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.

Comparative physical mapping between Oryza sativa (AA genome type) and O. punctata (BB genome type)

A comparative physical map of the AA genome (Oryza sativa) and the BB genome (O. punctata) was constructed by aligning a physical map of O. punctata, deduced from 63,942 BAC end sequences (BESs) and 34,224 fingerprints, onto the O. sativa genome sequence. The level of conservation of each chromosome between the two species was determined by calculating a ratio of BES alignments. The alignment result suggests more divergence of intergenic and repeat regions in comparison to gene-rich regions. Further, this characteristic enabled localization of heterochromatic and euchromatic regions for each chromosome of both species. The alignment identified 16 locations containing expansions, contractions, inversions, and transpositions. By aligning 40% of the punctata BES on the map, 87% of the punctata FPC map covered 98% of the O. sativa genome sequence. The genome size of O. punctata was estimated to be 8% larger than that of O. sativa with individual chromosome differences of 1.5-16.5%. The sum of expansions and contractions observed in regions >500 kb were similar, suggesting that most of the contractions/expansions contributing to the genome size difference between the two species are small, thus preserving the macro-collinearity between these species, which diverged approximately 2 million years ago.

A BAC-based physical map of the apple genome

Genome-wide physical mapping is an essential step toward investigating the genetic basis of complex traits as well as pursuing genomics research of virtually all plant and animal species. We have constructed a physical map of the apple genome from a total of 74,281 BAC clones representing approximately 10.5x haploid genome equivalents. The physical map consists of 2702 contigs, and it is estimated to span approximately 927 Mb in physical length. The reliability of contig assembly was evaluated by several methods, including assembling contigs using variable stringencies, assembling contigs using fingerprints from individual libraries, checking consensus maps of contigs, and using DNA markers. Altogether, the results demonstrated that the contigs were properly assembled. The apple genome-wide BAC-based physical map represents the first draft genome sequence not only for any member of the large Rosaceae family, but also for all tree species. This map will play a critical role in advanced genomics research for apple and other tree species, including marker development in targeted chromosome regions, fine-mapping and isolation of genes/QTL, conducting comparative genomics analyses of plant chromosomes, and large-scale genomics sequencing.

An integrated BAC and genome sequence physical map of Phytophthora sojae

Phytophthora spp. are serious pathogens that threaten numerous cultivated crops, trees, and natural vegetation worldwide. The soybean pathogen P. sojae has been developed as a model oomycete. Here, we report a bacterial artificial chromosome (BAC)-based, integrated physical map of the P. sojae genome. We constructed two BAC libraries, digested 8,681 BACs with seven restriction enzymes, end labeled the digested fragments with four dyes, and analyzed them with capillary electrophoresis. Fifteen data sets were constructed from the fingerprints, using individual dyes and all possible combinations, and were evaluated for contig assembly. In all, 257 contigs were assembled from the XhoI data set, collectively spanning approximately 132 Mb in physical length. The BAC contigs were integrated with the draft genome sequence of P. sojae by end sequencing a total of 1,440 BACs that formed a minimal tiling path. This enabled the 257 contigs of the BAC map to be merged with 207 sequence scaffolds to form an integrated map consisting of 79 superscaffolds. The map represents the first genome-wide physical map of a Phytophthora sp. and provides a valuable resource for genomics and molecular biology research in P. sojae and other Phytophthora spp. In one illustration of this value, we have placed the 350 members of a superfamily of putative pathogenicity effector genes onto the map, revealing extensive clustering of these genes.

Construction of a primary RH panel of Italian ryegrass genome via UV-induced protoplast fusion

Symmetric and asymmetric somatic hybrids were produced via protoplast fusion between common wheat ( TRITICUM AESTIVUM L.) cv. "Jinan 177" and Italian ryegrass ( LOLIUM MULTIFLORUM Lam.). The ryegrass without or with UV irradiation was used as a donor, providing a small amount of chromatin. In these somatic hybrids, most ryegrass chromosomes have been confirmed preferential elimination and the somatic hybrid calli and plants showed wheat-like morphology. Some of the hybrid lines were used for the analysis of distribution and heredity of donor DNA in the hybrid genome and the possibility of establishing a radiation hybrid (RH) panel of the ryegrass in the present experiment. These hybrids, subcultured for two and three years, retained the ryegrass DNA examined by RFLP and GISH analysis, respectively. Distribution of the ryegrass DNA in the wheat genomes of 20 single-cell individuals, randomly selected from hybrid cell lines produced, were analyzed by 21 ryegrass genome specific SSR markers. The average frequencies of molecular marker retention in symmetric hybrid lines (UV 0), as well as asymmetric hybrid lines from UV 30 s and 1 min were 10.88, 15.48 and 33.86, respectively. It was suggested that the UV dose increased the introgression of donor DNA into wheat genome. The ryegrass SSR fragments in most asymmetric hybrid cell lines remained stable over a period of 2 approximately 3 years. This revealed that those asymmetric somatic hybrids are suitable for the introgression of ryegrass DNA into wheat, and for RH panel and RH mapping.

Physical organization of the major duplication onBrassica oleraceachromosome O6 revealed through fluorescence in situ hybridization withArabidopsisandBrassicaBAC probes

The close relationship between Brassica oleracea and Arabidopsis thaliana has been used to explore the genetic and physical collinearity of the two species, focusing on an inverted segmental chromosome duplication within linkage group O6 of B. oleracea. Genetic evidence suggests that these segments share a common origin with a region of Arabidopsis chromosome 1. Brassica oleracea and Arabidopsis bacterial artificial chromosome probes have been used for fluorescence in situ hybridization analysis of B. oleracea pachytene chromosomes to further characterize the inverted duplication. This has been highly effective in increasing the local resolution of the cytogenetic map. We have shown that the physical order of corresponding genetic markers is highly conserved between the duplicated regions in B. oleracea and the physical lengths of the regions at pachytene are similar, while the genetic distances are considerably different. The physical marker order is also well conserved between Arabidopsis and B. oleracea, with only one short inversion identified. Furthermore, the relative physical distances between the markers in one segment of B. oleracea and Arabidopsis have stayed approximately the same. The efficacy of using fluorescence in situ hybridization, together with other forms of physical and genetic mapping, for elucidating such issues relating to synteny is discussed.Key words: collinearity, cytogenetic map, pachytene chromosomes, Brassica, Arabidopsis.

A 'Chinese Spring' wheat (Triticum aestivum L.) bacterial artificial chromosome library and its use in the isolation of SSR markers for targeted genome regions

A bacterial artificial chromosome (BAC) library was constructed from the bread wheat (Triticum aestivum L.) genotype 'Chinese Spring' ('CS'). The library consists of 395,136 clones with an estimated average insert size of 157 kb. This library provides an estimated 3.4-fold genome coverage for this hexaploid species. The genome coverage was confirmed by RFLP analysis of single-copy RFLP clones. The CS BAC library was used to develop simple sequence repeat (SSR) markers for targeted genome regions using five sequence-tagged-site (STS) markers designed from the chromosome arm of 3BS. The SSR markers for the targeted genome region were successfully obtained. However, similar numbers of new SSR markers were also generated for the other two homologous group 3 chromosomes. This data suggests that BAC clones belonging to all three chromosomes of homologous group 3 were isolated using the five STS primers. The potential impacts of these results on marker isolation in wheat and on library screening in general are discussed.

A fine physical map of the rice chromosome 5

A fine physical map of the rice (Oryza sativa spp. Japonica var. Nipponbare) chromosome 5 with bacterial artificial chromosome (BAC) and PI-derived artificial chromosome (PAC) clones was constructed through integration of 280 sequenced BAC/PAC clones and 232 sequence tagged site/expressed sequence tag markers with the use of fingerprinted contig data of the Nipponbare genome. This map consists of five contigs covering 99% of the estimated chromosome size (30.08 Mb). The four physical gaps were estimated at 30 and 20 kb for gaps 1-3 and gap 4, respectively. We have submitted 42.2-Mb sequences with 29.8 Mb of nonoverlapping sequences to public databases. BAC clones corresponding to telomere and centromere regions were confirmed by BAC-fluorescence in situ hybridization (FISH) on a pachytene chromosome. The genetically centromeric region at 54.6 cM was covered by a minimum tiling path spanning 2.1 Mb with no physical gaps. The precise position of the centromere was revealed by using three overlapping BAC/PACs for approximately 150 kb. In addition, FISH results revealed uneven chromatin condensation around the centromeric region at the pachytene stage. This map is of use for positional cloning and further characterization of the rice functional genomics.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

Whole-genome fractionation rapidly purifies DNA from centromeric regions

The condensed centromeric regions of higher eukaryotic chromosomes contain satellite sequences, transposons and retroelements, as well as transcribed genes that perform a variety of functions. These chromosomal domains nucleate kinetochores, mediate sister chromatid cohesion and inhibit recombination, yet their characterization has often lagged behind that of chromosome arms. Here, we describe a whole-genome fractionation technique that rapidly identifies bacterial artificial chromosome (BAC) clones derived from plant centromeric regions. This approach, which relies on hybridization of methylated genomic DNA, revealed BACs that correspond to the genetically mapped and sequenced Arabidopsis thaliana centromeric regions. Extending this method to other species in the Brassicaceae family identified centromere-linked clones and provided genome-wide estimates of methylated DNA abundance. Sequencing these clones will elucidate the changes that occur during plant centromere evolution. This genomic fractionation technique could identify centromeric DNA in genomes with similar methylation and repetitive DNA content, including those from crops and mammals.

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.