A high-resolution HAPPY map of Dictyostelium discoideum chromosome 6

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

Assembly of complex genomes using short reads remains a major challenge, which usually yields highly fragmented assemblies. Generation of ultradense linkage maps is promising for anchoring such assemblies, but traditional linkage mapping methods are hindered by the infrequency and unevenness of meiotic recombination that limit attainable map resolution. Here we develop a sequencing-based "in vitro" linkage mapping approach (called RadMap), where chromosome breakage and segregation are realized by generating hundreds of "subhaploid" fosmid/bacterial-artificial-chromosome clone pools, and by restriction site-associated DNA sequencing of these clone pools to produce an ultradense whole-genome restriction map to facilitate genome scaffolding. A bootstrap-based minimum spanning tree algorithm is developed for grouping and ordering of genome-wide markers and is implemented in a user-friendly, integrated software package (AMMO). We perform extensive analyses to validate the power and accuracy of our approach in the model plant Arabidopsis thaliana and human. We also demonstrate the utility of RadMap for enhancing the contiguity of a variety of whole-genome shotgun assemblies generated using either short Illumina reads (300 bp) or long PacBio reads (6-14 kb), with up to 15-fold improvement of N50 (∼816 kb-3.7 Mb) and high scaffolding accuracy (98.1-98.5%). RadMap outperforms BioNano and Hi-C when input assembly is highly fragmented (contig N50 = 54 kb). RadMap can capture wide-range contiguity information and provide an efficient and flexible tool for high-resolution physical mapping and scaffolding of highly fragmented assemblies.

A deep coverage Dictyostelium discoideum genomic DNA library replicates stably in Escherichia coli

The natural history of the amoeba Dictyostelium discoideum has inspired scientific inquiry for seventy-five years. A genetically tractable haploid eukaryote, D. discoideum appeals as a laboratory model as well. However, certain rote molecular genetic tasks, such as PCR and cloning, are difficult due to the AT-richness and low complexity of its genome. Here we report on the construction of a ~20 fold coverage D. discoideum genomic library in Escherichia coli, cloning 4-10 kilobase partial restriction fragments into a linear vector. End-sequencing indicates that most clones map to the six chromosomes in an unbiased distribution. Over 70% of these clones contain at least one complete open reading frame. We demonstrate that individual clones and library composition are stable over multiple replication cycles. Our library will enable numerous molecular biological applications and the completion of additional species' genome sequences, and suggests a path towards the long-elusive goal of genetic complementation.

Insights into the genome structure and copy-number variation of Eimeria tenella

BACKGROUND

Eimeria is a genus of parasites in the same phylum (Apicomplexa) as human parasites such as Toxoplasma, Cryptosporidium and the malaria parasite Plasmodium. As an apicomplexan whose life-cycle involves a single host, Eimeria is a convenient model for understanding this group of organisms. Although the genomes of the Apicomplexa are diverse, that of Eimeria is unique in being composed of large alternating blocks of sequence with very different characteristics - an arrangement seen in no other organism. This arrangement has impeded efforts to fully sequence the genome of Eimeria, which remains the last of the major apicomplexans to be fully analyzed. In order to increase the value of the genome sequence data and aid in the effort to gain a better understanding of the Eimeria tenella genome, we constructed a whole genome map for the parasite.

RESULTS

A total of 1245 contigs representing 70.0% of the whole genome assembly sequences (Wellcome Trust Sanger Institute) were selected and subjected to marker selection. Subsequently, 2482 HAPPY markers were developed and typed. Of these, 795 were considered as usable markers, and utilized in the construction of a HAPPY map. Markers developed from chromosomally-assigned genes were then integrated into the HAPPY map and this aided the assignment of a number of linkage groups to their respective chromosomes. BAC-end sequences and contigs from whole genome sequencing were also integrated to improve and validate the HAPPY map. This resulted in an integrated HAPPY map consisting of 60 linkage groups that covers approximately half of the estimated 60 Mb genome. Further analysis suggests that the segmental organization first seen in Chromosome 1 is present throughout the genome, with repeat-poor (P) regions alternating with repeat-rich (R) regions. Evidence of copy-number variation between strains was also uncovered.

CONCLUSIONS

This paper describes the application of a whole genome mapping method to improve the assembly of the genome of E. tenella from shotgun data, and to help reveal its overall structure. A preliminary assessment of copy-number variation (extra or missing copies of genomic segments) between strains of E. tenella was also carried out. The emerging picture is of a very unusual genome architecture displaying inter-strain copy-number variation. We suggest that these features may be related to the known ability of this parasite to rapidly develop drug resistance.

Development and initial characterization of a HAPPY panel for mapping the X. tropicalis genome

HAPPY mapping was designed to pursue the analysis of approximately random HAPloid DNA breakage samples using the PolYmerase chain reaction for mapping genomes. In the present study, we improved the method and integrated two other molecular techniques into the process: whole genome amplification and the Sequenom SNP (single nucleotide polymorphism) genotyping assay in order to facilitate whole genome mapping of X. tropicalis. The former technique amplified enough DNA materials to genotype a large number of markers, while the latter allowed for relatively high throughput marker genotyping with multiplex assays on the HAPPY lines. A total of 58 X. tropicalis genes were genotyped on an initial panel of 383 HAPPY lines, which contributed to formation of a working panel of 146 lines. Further genotyping of 29 markers on the working panel led to construction of a HAPPY map for the X. tropicalis genome. We believe that our improved HAPPY method described in the present study has paved the way for the community to map different genomes with a simple, but powerful approach.

Enhancing radiation hybrid mapping through whole genome amplification

Radiation hybrid (RH) mapping is limited by the inherent genomic instability of RH clones entailing both, limited DNA sample amounts and genomic heterogeneity of the clones. Here the instability of RH clones is quantified and the suitability of the multiple strand displacement whole genome amplification method (WGA) for radiation hybrid mapping is assessed. To quantify the instability of RH clones, eleven clones of a 10,000(Rad) rhesus macaque radiation hybrid panel were passaged ten times and analyzed by interspersed repeat sequence specific quantitative PCR and by genotyping of 46 macaque chromosome 5 STS markers. The quantitative PCR data indicate an average loss of 55% of the donor DNA over 10 passages. Over the same period, a dropout of 46.2% of the STS markers was observed. These data indicate a genome wide half-life of the donor DNA of 8.7 passages and of 10.6 passages for the chromosome 5 markers. The genotyping data of the genomic RH DNA were compared to three sets of WGA experiments: 1) single wgaDNA amplifications, 2) six WGA replicates, and 3) re-amplification of wga DNA. The assays demonstrated concordance rates of 97.6%, 98% and 99.3%, respectively, and indicated the marker specificity of some repeated WGA dropouts. The study confirms that WGA is suitable for RH mapping studies should enable the accurate analysis of almost an infinite numbers of markers. WGA will allow the analysis of earliest RH clone passages, thus limiting their heterogeneity and RH mapping artifacts.

Old can be new again: HAPPY whole genome sequencing, mapping and assembly

During the last three decades, both genome mapping and sequencing methods have advanced significantly to provide a foundation for scientists to understand genome structures and functions in many species. Generally speaking, genome mapping relies on genome sequencing to provide basic materials, such as DNA probes and markers for their localizations, thus constructing the maps. On the other hand, genome sequencing often requires a high-resolution map as a skeleton for whole genome assembly. However, both genome mapping and sequencing have never come together in one pipeline. After reviewing mapping and next-generation sequencing methods, we would like to share our thoughts with the genome community on how to combine the HAPPY mapping technique with the new-generation sequencing, thus integrating two systems into one pipeline, called HAPPY pipeline. The pipeline starts with preparation of a HAPPY panel, followed by multiple displacement amplification for producing a relatively large quantity of DNA. Instead of conventional marker genotyping, the amplified panel DNA samples are subject to new-generation sequencing with barcode method, which allows us to determine the presence/absence of a sequence contig as a traditional marker in the HAPPY panel. Statistical analysis will then be performed to infer how close or how far away from each other these contigs are within a genome and order the whole genome sequence assembly as well. We believe that such a universal approach will play an important role in genome sequencing, mapping, and assembly of many species; thus advancing genome science and its applications in biomedicine and agriculture.

Sequencing and analysis of chromosome 1 of Eimeria tenella reveals a unique segmental organization

Eimeria tenella is an intracellular protozoan parasite that infects the intestinal tracts of domestic fowl and causes coccidiosis, a serious and sometimes lethal enteritis. Eimeria falls in the same phylum (Apicomplexa) as several human and animal parasites such as Cryptosporidium, Toxoplasma, and the malaria parasite, Plasmodium. Here we report the sequencing and analysis of the first chromosome of E. tenella, a chromosome believed to carry loci associated with drug resistance and known to differ between virulent and attenuated strains of the parasite. The chromosome--which appears to be representative of the genome--is gene-dense and rich in simple-sequence repeats, many of which appear to give rise to repetitive amino acid tracts in the predicted proteins. Most striking is the segmentation of the chromosome into repeat-rich regions peppered with transposon-like elements and telomere-like repeats, alternating with repeat-free regions. Predicted genes differ in character between the two types of segment, and the repeat-rich regions appear to be associated with strain-to-strain variation.

Understanding anthelmintic resistance: The need for genomics and genetics

Anthelmintic resistance is a major problem for the control of many parasitic nematode species and has become a major constraint to livestock production in many parts of the world. In spite of its increasing importance, there is still a poor understanding of the molecular and genetic basis of resistance. It is unclear which mutations contribute most to the resistance phenotype and how resistance alleles arise, are selected and spread in parasite populations. The main strategy used to identify mutations responsible for anthelmintic resistance has been to undertake experimental studies on candidate genes. These genes have been chosen predominantly on the basis of our knowledge of drug mode-of-action and the identification of mutations that can confer resistance in model organisms. The application of these approaches to the analysis of benzimidazole and ivermectin resistance is reviewed and the reasons for their relative success or failure are discussed. The inherent limitation of candidate gene studies is that they rely on very specific and narrow assumptions about the likely identity of resistance-associated genes. In contrast, forward genetic and functional genomic approaches do not make such assumptions, as illustrated by the successful application of these techniques in the study of insecticide resistance. Although there is an urgent need to apply these powerful approaches to anthelmintic resistance research, the basic methodologies and resources are still lacking. However, these are now being developed for the trichostrongylid nematode Haemonchus contortus and the current progress and research priorities in this area are discussed.

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.

Microarray-based comparative genomic analyses of the human malaria parasite Plasmodium falciparum using Affymetrix arrays

Microarray-based comparative genomic hybridization (CGH) provides a powerful tool for whole genome analyses and the rapid detection of genomic variation that underlies virulence and disease. In the field of Plasmodium research, many of the parasite genomes that one might wish to study in a high throughput manner are not laboratory clones, but clinical isolates. One of the key limitations to the use of clinical samples in CGH, however, is the miniscule amounts of genomic DNA available. Here we describe the successful application of multiple displacement amplification (MDA), a non-PCR-based amplification method that exhibits clear advantages over all other currently available methods. Using MDA, CGH was performed on a panel of NF54 and IT/FCR3 clones, identifying previously published deletions on chromosomes 2 and 9 as well as polymorphism in genes associated with disease pathology.

The genome of the social amoeba Dictyostelium discoideum

The social amoebae are exceptional in their ability to alternate between unicellular and multicellular forms. Here we describe the genome of the best-studied member of this group, Dictyostelium discoideum. The gene-dense chromosomes of this organism encode approximately 12,500 predicted proteins, a high proportion of which have long, repetitive amino acid tracts. There are many genes for polyketide synthases and ABC transporters, suggesting an extensive secondary metabolism for producing and exporting small molecules. The genome is rich in complex repeats, one class of which is clustered and may serve as centromeres. Partial copies of the extrachromosomal ribosomal DNA (rDNA) element are found at the ends of each chromosome, suggesting a novel telomere structure and the use of a common mechanism to maintain both the rDNA and chromosomal termini. A proteome-based phylogeny shows that the amoebozoa diverged from the animal-fungal lineage after the plant-animal split, but Dictyostelium seems to have retained more of the diversity of the ancestral genome than have plants, animals or fungi.

Revamp a model-status and prospects of the Dictyostelium genome project

International efforts are underway that aim at determining the complete genome sequence of the social amoeba Dictyostelium discoideum. As strategy, a whole chromosome shotgun (WCS) approach was chosen and each of the six Dictyostelium chromosomes was assigned to project partners. The project is well advanced, chromosome 2 was recently published, and it is expected that the sequences of chromosomes 1 and 6 and a gene catalogue for the complete genome will be available at the end of this year. The genome sequence, together with powerful molecular genetic tools, will undoubtedly further accelerate Dictyostelium research into a number of fundamental biological processes that are common to a wide range of eukaryotes. Furthermore, it will constitute the basis for genome-wide functional analyses. The integration of results from these studies should ultimately lead to a better understanding of the complex networks that govern cellular behavior and development.

Integrated mapping, chromosomal sequencing and sequence analysis of Cryptosporidium parvum

The apicomplexan Cryptosporidium parvum is one of the most prevalent protozoan parasites of humans. We report the physical mapping of the genome of the Iowa isolate, sequencing and analysis of chromosome 6, and approximately 0.9 Mbp of sequence sampled from the remainder of the genome. To construct a robust physical map, we devised a novel and general strategy, enabling accurate placement of clones regardless of clone artefacts. Analysis reveals a compact genome, unusually rich in membrane proteins. As in Plasmodium falciparum, the mean size of the predicted proteins is larger than that in other sequenced eukaryotes. We find several predicted proteins of interest as potential therapeutic targets, including one exhibiting similarity to the chloroquine resistance protein of Plasmodium. Coding sequence analysis argues against the conventional phylogenetic position of Cryptosporidium and supports an earlier suggestion that this genus arose from an early branching within the Apicomplexa. In agreement with this, we find no significant synteny and surprisingly little protein similarity with Plasmodium. Finally, we find two unusual and abundant repeats throughout the genome. Among sequenced genomes, one motif is abundant only in C. parvum, whereas the other is shared with (but has previously gone unnoticed in) all known genomes of the Coccidia and Haemosporida. These motifs appear to be unique in their structure, distribution and sequences.

Crawling into a new era-the Dictyostelium genome project

The social amoeba Dictyostelium discoideum is a well-established model organism for the study of basic aspects of differentiation, signal transduction, phagocytosis, cytokinesis and cell motility. Its genome is being sequenced by an international consortium using a whole chromosome shotgun approach. The pacemaker of the D.discoideum genome project has been chromosome 2, the largest chromosome, which at 8 Mb represents approximately 25% of the genome and whose sequence and analysis have been published recently. Chromosomes 1 and 6 are close to being finished. To accelerate completion of the genome sequence, the next step in the project will be a whole-genome assembly followed by the analysis of the complete gene content. The completed genome sequence and its analysis provide the basis for genome-wide functional studies. It will position Dictyostelium at the same level as other model organisms and further enhance its experimental attractiveness.

HAPPY mapping in a plant genome: reconstruction and analysis of a high-resolution physical map of a 1.9 Mbp region of Arabidopsis thaliana chromosome 4

HAPPY mapping is an in vitro approach for defining the order and spacing of DNA markers directly on native genomic DNA. This cloning-free technique is based on analysing the segregation of markers amplified from high molecular weight genomic DNA which has been broken randomly and 'segregated' by limiting dilution into subhaploid samples. It is a uniquely versatile tool, allowing for the construction of genome maps with flexible ranges and resolutions. Moreover, it is applicable to plant genomes, for which many of the techniques pioneered in animal genomes are inapplicable or inappropriate. We report here its demonstration in a plant genome by reconstructing the physical map of a 1.9 Mbp region around the FCA locus of Arabidopsis thaliana. The resulting map, spanning around 10% of chromosome 4, is in excellent agreement with the DNA sequence and has a mean marker spacing of 16 kbp. We argue that HAPPY maps of any required resolution can be made immediately and with relatively little effort for most plant species and, furthermore, that such maps can greatly aid the construction of regional or genome-wide physical maps.

Chemotaxis and cell differentiation in Dictyostelium

Dictyostelium genome sequencing predicts an unexpectedly large number of genes. Many are absent from yeast but present in animals and presumably support cellular abilities not found in yeast. Prominent amongst these abilities is chemotaxis, where great strides are being made in understanding how cells orient in a gradient and mobilise their cytoskeleton for movement. In multicellular development, a regulatory scheme for proportioning prespore and prestalk-O cells has emerged.

Sequence of Plasmodium falciparum chromosomes 1, 3-9 and 13

Since the sequencing of the first two chromosomes of the malaria parasite, Plasmodium falciparum, there has been a concerted effort to sequence and assemble the entire genome of this organism. Here we report the sequence of chromosomes 1, 3-9 and 13 of P. falciparum clone 3D7--these chromosomes account for approximately 55% of the total genome. We describe the methods used to map, sequence and annotate these chromosomes. By comparing our assemblies with the optical map, we indicate the completeness of the resulting sequence. During annotation, we assign Gene Ontology terms to the predicted gene products, and observe clustering of some malaria-specific terms to specific chromosomes. We identify a highly conserved sequence element found in the intergenic region of internal var genes that is not associated with their telomeric counterparts.

Synteny comparison between apes and human using fine-mapping of the genome

Comparing the genomes of the great apes and human should provide novel information concerning the origins of humankind. Relative to the great apes, the human karyotype has one fewer chromosome pair, as human chromosome 2 derived from the telomeric fusion of two ancestral primate chromosomes. To identify the genomic rearrangements that accompanied human speciation, we initiated a comparative study between human, chimpanzee, and gorilla. Using the HAPPY mapping method, an acellular adaptation of the radiation hybrid method, we mapped a few hundred markers on the human, chimpanzee, and gorilla genomes. This allowed us to identify several chromosome rearrangements, in particular a pericentric inversion and a translocation. We precisely localized the synteny breakpoint that led to the formation of human chromosome 2. This breakpoint was confirmed by FISH mapping.

Sequence and analysis of chromosome 2 of Dictyostelium discoideum

The genome of the lower eukaryote Dictyostelium discoideum comprises six chromosomes. Here we report the sequence of the largest, chromosome 2, which at 8 megabases (Mb) represents about 25% of the genome. Despite an A + T content of nearly 80%, the chromosome codes for 2,799 predicted protein coding genes and 73 transfer RNA genes. This gene density, about 1 gene per 2.6 kilobases (kb), is surpassed only by Saccharomyces cerevisiae (one per 2 kb) and is similar to that of Schizosaccharomyces pombe (one per 2.5 kb). If we assume that the other chromosomes have a similar gene density, we can expect around 11,000 genes in the D. discoideum genome. A significant number of the genes show higher similarities to genes of vertebrates than to those of other fully sequenced eukaryotes. This analysis strengthens the view that the evolutionary position of D. discoideum is located before the branching of metazoa and fungi but after the divergence of the plant kingdom, placing it close to the base of metazoan evolution.

Genomic database resources for Dictyostelium discoideum

Dictyostelium is an attractive model system for the study of mechanisms basic to cellular function or complex multicellular developmental processes. Recent advances in Dictyostelium genomics have generated a wide spectrum of resources. However, much of the current genomic sequence information is still not currently available through GenBank or related databases. Thus, many investigators are unaware that extensive sequence data from Dictyostelium has been compiled, or of its availability and access. Here, we discuss progress in Dictyostelium genomics and gene annotation, and highlight the primary portals for sequence access, manipulation and analysis (http://genome.imb-jena.de/dictyostelium/; http://dictygenome.bcm.tmc.edu/; http://www.sanger. ac.uk/Projects/D_discoideum/; http://www.csm.biol. tsukuba.ac.jp/cDNAproject.html).