A BAC-based physical map of the major autosomes of Drosophila melanogaster

Drosophila melanogaster: a case study of a model genomic sequence and its consequences

The sequencing and annotation of the Drosophila melanogaster genome, first published in 2000 through collaboration between Celera Genomics and the Drosophila Genome Projects, has provided a number of important contributions to genome research. By demonstrating the utility of methods such as whole-genome shotgun sequencing and genome annotation by a community "jamboree," the Drosophila genome established the precedents for the current paradigm used by most genome projects. Subsequent releases of the initial genome sequence have been improved by the Berkeley Drosophila Genome Project and annotated by FlyBase, the Drosophila community database, providing one of the highest-quality genome sequences and annotations for any organism. We discuss the impact of the growing number of genome sequences now available in the genus on current Drosophila research, and some of the biological questions that these resources will enable to be solved in the future.

Analysis of pyrimidine catabolism in Drosophila melanogaster using epistatic interactions with mutations of pyrimidine biosynthesis and beta-alanine metabolism

The biochemical pathway for pyrimidine catabolism links the pathways for pyrimidine biosynthesis and salvage with beta-alanine metabolism, providing an array of epistatic interactions with which to analyze mutations of these pathways. Loss-of-function mutations have been identified and characterized for each of the enzymes for pyrimidine catabolism: dihydropyrimidine dehydrogenase (DPD), su(r) mutants; dihydropyrimidinase (DHP), CRMP mutants; beta-alanine synthase (betaAS), pyd3 mutants. For all three genes, mutants are viable and fertile and manifest no obvious phenotypes, aside from a variety of epistatic interactions. Mutations of all three genes disrupt suppression by the rudimentary gain-of-function mutation (r(Su(b))) of the dark cuticle phenotype of black mutants in which beta-alanine pools are diminished; these results confirm that pyrimidines are the major source of beta-alanine in cuticle pigmentation. The truncated wing phenotype of rudimentary mutants is suppressed completely by su(r) mutations and partially by CRMP mutations; however, no suppression is exhibited by pyd3 mutations. Similarly, su(r) mutants are hypersensitive to dietary 5-fluorouracil, CRMP mutants are less sensitive, and pyd3 mutants exhibit wild-type sensitivity. These results are discussed in the context of similar consequences of 5-fluoropyrimidine toxicity and pyrimidine catabolism mutations in humans.

Genetic interactions of separase regulatory subunits reveal the diverged Drosophila Cenp-C homolog

Faithful transmission of genetic information during mitotic divisions depends on bipolar attachment of sister kinetochores to the mitotic spindle and on complete resolution of sister-chromatid cohesion immediately before the metaphase-to-anaphase transition. Separase is thought to be responsible for sister-chromatid separation, but its regulation is not completely understood. Therefore, we have screened for genetic loci that modify the aberrant phenotypes caused by overexpression of the regulatory separase complex subunits Pimples/securin and Three rows in Drosophila. An interacting gene was found to encode a constitutive centromere protein. Characterization of its centromere localization domain revealed the presence of a diverged CENPC motif. While direct evidence for an involvement of this Drosophila Cenp-C homolog in separase activation at centromeres could not be obtained, in vivo imaging clearly demonstrated that it is required for normal attachment of kinetochores to the spindle.

The S6KII (rsk) gene of Drosophila melanogaster differentially affects an operant and a classical learning task

In an attempt to dissect classical and operant conditioning in Drosophila melanogaster, we have isolated the gene for ribosomal S6 kinase II (S6KII). This enzyme is part of a family of serine-threonine kinases that in mammals have been implicated in the MAPK (mitogen-activated protein kinase) signaling cascade controlling (among other processes) synaptic plasticity (long-term potentiation/long-term depression) and memory formation. The human homolog rsk2 has been linked to mental retardation (Coffin-Lowry syndrome). Mutant analysis in Drosophila shows that S6KII serves different functions in operant place learning and classical (pavlovian) olfactory conditioning. Whereas in the null mutant only pavlovian olfactory learning is affected, a P-element insertion mutant reducing the amount of S6KII only affects operant place learning. A mutant lacking part of the N-terminal kinase domain and performing poorly in both learning tasks is dominant in the operant paradigm and recessive in the pavlovian paradigm. The behavioral defects in the pavlovian task can be rescued by the genomic S6KII transgene. Overexpression of S6KII in wild type has a dominant-negative effect on the operant task that is rescued by the null mutant, whereas in the pavlovian task overexpression may even enhance learning performance.

A BAC-based physical map of the Drosophila buzzatii genome

Large-insert genomic libraries facilitate cloning of large genomic regions, allow the construction of clone-based physical maps, and provide useful resources for sequencing entire genomes. Drosophila buzzatii is a representative species of the repleta group in the Drosophila subgenus, which is being widely used as a model in studies of genome evolution, ecological adaptation, and speciation. We constructed a Bacterial Artificial Chromosome (BAC) genomic library of D. buzzatii using the shuttle vector pTARBAC2.1. The library comprises 18,353 clones with an average insert size of 152 kb and an approximately 18x expected representation of the D. buzzatii euchromatic genome. We screened the entire library with six euchromatic gene probes and estimated the actual genome representation to be approximately 23x. In addition, we fingerprinted by restriction digestion and agarose gel electrophoresis a sample of 9555 clones, and assembled them using FingerPrint Contigs (FPC) software and manual editing into 345 contigs (mean of 26 clones per contig) and 670 singletons. Finally, we anchored 181 large contigs (containing 7788 clones) to the D. buzzatii salivary gland polytene chromosomes by in situ hybridization of 427 representative clones. The BAC library and a database with all the information regarding the high coverage BAC-based physical map described in this paper are available to the research community.

Chriz, a chromodomain protein specific for the interbands of Drosophila melanogaster polytene chromosomes

Polytene interphase chromosomes are compacted into a series of bands and interbands reflecting their organization into independent chromosomal domains. In order to understand chromosomal organization, we set out to study the role of proteins that are selective for interbands. Here we describe the Drosophila melanogaster chromodomain protein Chriz that is coimmunoprecipitated with the zinc finger protein Z4. Both proteins colocalize exclusively to the interbands on Drosophila polytene chromosomes. Like Z4, Chriz is ubiquitously expressed throughout development and is associated with chromatin in all interphase nuclei. Following dissociation from chromatin, early in mitosis Chriz binds to the centrosomes and to the mitotic spindle. Newly induced amorphic Chriz alleles are early lethal, and ubiquitous overexpression of Chriz is lethal as well. Available Chriz hypomorphs which survive until pupal stage have a normal chromosomal phenotype. Reducing Z4 protein does not affect Chriz binding to polytene chromosomes and vice versa. Z4 is still chromosomally bound when Chriz protein is depleted by RNA interference.

Genome physical mapping with large-insert bacterial clones by fingerprint analysis: methodologies, source clone genome coverage, and contig map quality

Genome physical mapping with large-insert clones by fingerprint analysis is becoming an active area of genomics research. Here, we report two new capillary electrophoresis-based fingerprinting methods for genome physical mapping and the effects of different fingerprinting methods and source clone genome coverage on quality physical map construction revealed by computer simulations and laboratory experiments. It was shown that the manual sequencing gel-based two-enzyme fingerprinting method consistently generated larger and more accurate contigs, followed by the new capillary electrophoresis-based three-enzyme method, the new capillary electrophoresis-based five-enzyme (SNaPshot) method, the agarose gel-based one-enzyme method, and the automatic sequencing gel-based four-enzyme method, in descending order, when 1% or fewer questionable clones were allowed. Analysis of clones equivalent to 5x, 8x, 10x, and 15x genomes using the fingerprinting methods revealed that as the number of clones increased from 5x to 10x, the contig length rapidly increased for all methods. However, when the number of clones was increased from 10x to 15x coverage, the contig length at best increased at a lower rate or even decreased. The results will provide useful knowledge and strategies for effective construction of quality genome physical maps for advanced genomics research.

Genome physical mapping from large-insert clones by fingerprint analysis with capillary electrophoresis: a robust physical map of Penicillium chrysogenum

Physical mapping with large-insert clones is becoming an active area of genomics research, and capillary electrophoresis (CE) promises to revolutionize the physical mapping technology. Here, we demonstrate the utility of the CE technology for genome physical mapping with large-insert clones by constructing a robust, binary bacterial artificial chromosome (BIBAC)-based physical map of Penicillium chrysogenum. We fingerprinted 23.1x coverage BIBAC clones with five restriction enzymes and the SNaPshot kit containing four fluorescent-ddNTPs using the CE technology, and explored various strategies to construct quality physical maps. It was shown that the fingerprints labeled with one or two colors, resulting in 40-70 bands per clone, were assembled into much better quality maps than those labeled with three or four colors. The selection of fingerprinting enzymes was crucial to quality map construction. From the dataset labeled with ddTTP-dROX, we assembled a physical map for P.chrysogenum, with 2-3 contigs per chromosome and anchored the map to its chromosomes. This map represents the first physical map constructed using the CE technology, thus providing not only a platform for genomic studies of the penicillin-producing species, but also strategies for efficient use of the CE technology for genome physical mapping of plants, animals and microbes.

Meiotic recombination in Drosophila females depends on chromosome continuity between genetically defined boundaries

In the pairing-site model, specialized regions on each chromosome function to establish meiotic homolog pairing. Analysis of these sites could provide insights into the mechanism used by Drosophila females to form a synaptonemal complex (SC) in the absence of meiotic recombination. These specialized sites were first established on the X chromosome by noting that there were barriers to crossover suppression caused by translocation heterozygotes. These sites were genetically mapped and proposed to be pairing sites. By comparing the cytological breakpoints of third chromosome translocations to their patterns of crossover suppression, we have mapped two sites on chromosome 3R. We have performed experiments to determine if these sites have a role in meiotic homolog pairing and the initiation of recombination. Translocation heterozygotes exhibit reduced gene conversion within the crossover-suppressed region, consistent with an effect on the initiation of meiotic recombination. To determine if homolog pairing is disrupted in translocation heterozygotes, we used fluorescent in situ hybridization to measure the extent of homolog pairing. In wild-type oocytes, homologs are paired along their entire lengths prior to accumulation of the SC protein C(3)G. Surprisingly, translocation heterozygotes exhibited homolog pairing similar to wild type within the crossover-suppressed regions. This result contrasted with our observations of c(3)G mutant females, which were found to be defective in pairing. We propose that each Drosophila chromosome is divided into several domains by specialized sites. These sites are not required for homolog pairing. Instead, the initiation of meiotic recombination requires continuity of the meiotic chromosome structure within each of these domains.

Genetic and bioinformatic analysis of 41C and the 2R heterochromatin of Drosophila melanogaster: a window on the heterochromatin-euchromatin junction

Genomic sequences provide powerful new tools in genetic analysis, making it possible to combine classical genetics with genomics to characterize the genes in a particular chromosome region. These approaches have been applied successfully to the euchromatin, but analysis of the heterochromatin has lagged somewhat behind. We describe a combined genetic and bioinformatics approach to the base of the right arm of the Drosophila melanogaster second chromosome, at the boundary between pericentric heterochromatin and euchromatin. We used resources provided by the genome project to derive a physical map of the region, examine gene density, and estimate the number of potential genes. We also carried out a large-scale genetic screen for lethal mutations in the region. We identified new alleles of the known essential genes and also identified mutations in 21 novel loci. Fourteen complementation groups map proximal to the assembled sequence. We used PCR to map the endpoints of several deficiencies and used the same set of deficiencies to order the essential genes, correlating the genetic and physical map. This allowed us to assign two of the complementation groups to particular "computed/curated genes" (CGs), one of which is Nipped-A, which our evidence suggests encodes Drosophila Tra1/TRRAP.

Characterization of the grappa gene, the Drosophila histone H3 lysine 79 methyltransferase

We have identified a novel gene named grappa (gpp) that is the Drosophila ortholog of the Saccharomyces cerevisiae gene Dot1, a histone methyltransferase that modifies the lysine (K)79 residue of histone H3. gpp is an essential gene identified in a genetic screen for dominant suppressors of pairing-dependent silencing, a Polycomb-group (Pc-G)-mediated silencing mechanism necessary for the maintenance phase of Bithorax complex (BX-C) expression. Surprisingly, gpp mutants not only exhibit Pc-G phenotypes, but also display phenotypes characteristic of trithorax-group mutants. Mutations in gpp also disrupt telomeric silencing but do not affect centric heterochromatin. These apparent contradictory phenotypes may result from loss of gpp activity in mutants at sites of both active and inactive chromatin domains. Unlike the early histone H3 K4 and K9 methylation patterns, the appearance of methylated K79 during embryogenesis coincides with the maintenance phase of BX-C expression, suggesting that there is a unique role for this chromatin modification in development.

General method for the modification of different BAC types and the rapid generation of BAC transgenic mice

Most genome projects have relied on the sequencing of bacterial artificial chromosomes (BACs), which encompass 100-300 kb of genomic DNA. As a consequence, several thousand BAC clones are now mapped to the human and mouse genome. It is therefore possible to identify in silico a BAC clone that carries a particular gene and obtain it commercially. Given the large size of BACs, most if not all regulatory sequences of a gene are present and can be used to direct faithful and tissue-specific expression of heterologous genes in vitro in cell cultures and in vivo in BAC-transgenic mice. We describe here an optimized and comprehensive protocol to select, modify, and purify BACs in order to generate BAC-transgenic mice. Importantly, this protocol includes a method to generate, within 2 days, complex plasmid cassettes required to modify BACs, and to efficiently modify different types of BACs selected from the two major BAC libraries available. Altogether, using a combination of genomic database analysis, overlap PCR cloning, and BAC recombination in bacteria, our approach allows for the rapid and reliable generation of "pseudo knockin" mice. genesis 38:39-50, 2004.

GDR (Genome Database for Rosaceae): integrated web resources for Rosaceae genomics and genetics research

BACKGROUND

Peach is being developed as a model organism for Rosaceae, an economically important family that includes fruits and ornamental plants such as apple, pear, strawberry, cherry, almond and rose. The genomics and genetics data of peach can play a significant role in the gene discovery and the genetic understanding of related species. The effective utilization of these peach resources, however, requires the development of an integrated and centralized database with associated analysis tools.

DESCRIPTION

The Genome Database for Rosaceae (GDR) is a curated and integrated web-based relational database. GDR contains comprehensive data of the genetically anchored peach physical map, an annotated peach EST database, Rosaceae maps and markers and all publicly available Rosaceae sequences. Annotations of ESTs include contig assembly, putative function, simple sequence repeats, and anchored position to the peach physical map where applicable. Our integrated map viewer provides graphical interface to the genetic, transcriptome and physical mapping information. ESTs, BACs and markers can be queried by various categories and the search result sites are linked to the integrated map viewer or to the WebFPC physical map sites. In addition to browsing and querying the database, users can compare their sequences with the annotated GDR sequences via a dedicated sequence similarity server running either the BLAST or FASTA algorithm. To demonstrate the utility of the integrated and fully annotated database and analysis tools, we describe a case study where we anchored Rosaceae sequences to the peach physical and genetic map by sequence similarity.

CONCLUSIONS

The GDR has been initiated to meet the major deficiency in Rosaceae genomics and genetics research, namely a centralized web database and bioinformatics tools for data storage, analysis and exchange. GDR can be accessed at http://www.genome.clemson.edu/gdr/.

Creation of a minimal tiling path of genomic clones for Drosophila: provision of a common resource

On the basis of shotgun subclone libraries used in the sequencing of the Drosophila melanogaster genome, a minimal tiling path of subclones across much of the genome was determined. About 320,000 shotgun clones for chromosomes X(12–20), 2R, 2L, 3R, and 4 were available from the Berkeley Drosophila Genome Project. The clone inserts have an average length of 3.4 kb and are amenable to standard PCR amplification. The resulting tiling path covers 86.2% of chromosome X(12–20), 86.2% of chromosomal arm 2R, 79.0% of 2L, 89.6% of 3R, and 80.5% of chromosome 4. In total, the 25,135 clones represent 76.7 Mb—equivalent to about 67% of the genome —and would be suitable for producing a microarray on a single slide.

Integrated and sequence-ordered BAC- and YAC-based physical maps for the rat genome

As part of the effort to sequence the genome of Rattus norvegicus, we constructed a physical map comprised of fingerprinted bacterial artificial chromosome (BAC) clones from the CHORI-230 BAC library. These BAC clones provide approximately 13-fold redundant coverage of the genome and have been assembled into 376 fingerprint contigs. A yeast artificial chromosome (YAC) map was also constructed and aligned with the BAC map via fingerprinted BAC and P1 artificial chromosome clones (PACs) sharing interspersed repetitive sequence markers with the YAC-based physical map. We have annotated 95% of the fingerprint map clones in contigs with coordinates on the version 3.1 rat genome sequence assembly, using BAC-end sequences and in silico mapping methods. These coordinates have allowed anchoring 358 of the 376 fingerprint map contigs onto the sequence assembly. Of these, 324 contigs are anchored to rat genome sequences localized to chromosomes, and 34 contigs are anchored to unlocalized portions of the rat sequence assembly. The remaining 18 contigs, containing 54 clones, still require placement. The fingerprint map is a high-resolution integrative data resource that provides genome-ordered associations among BAC, YAC, and PAC clones and the assembled sequence of the rat genome.

Integration of animal linkage and BAC contig maps using overgo hybridization

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

Genetic and Bioinformatic Analysis of 41C and the 2R Heterochromatin of Drosophila melanogaster: A Window on the Heterochromatin-Euchromatin Junction

Genomic sequences provide powerful new tools in genetic analysis, making it possible to combine classical genetics with genomics to characterize the genes in a particular chromosome region. These approaches have been applied successfully to the euchromatin, but analysis of the heterochromatin has lagged somewhat behind. We describe a combined genetic and bioinformatics approach to the base of the right arm of the Drosophila melanogaster second chromosome, at the boundary between pericentric heterochromatin and euchromatin. We used resources provided by the genome project to derive a physical map of the region, examine gene density, and estimate the number of potential genes. We also carried out a large-scale genetic screen for lethal mutations in the region. We identified new alleles of the known essential genes and also identified mutations in 21 novel loci. Fourteen complementation groups map proximal to the assembled sequence. We used PCR to map the endpoints of several deficiencies and used the same set of deficiencies to order the essential genes, correlating the genetic and physical map. This allowed us to assign two of the complementation groups to particular “computed/curated genes” (CGs), one of which is Nipped-A, which our evidence suggests encodes Drosophila Tra1/TRRAP.

A BAC- and BIBAC-based physical map of the soybean genome

Genome-wide physical maps are crucial to many aspects of advanced genome research. We report a genome-wide, bacterial artificial chromosome (BAC) and plant-transformation-competent binary large-insert plasmid clone (hereafter BIBAC)-based physical map of the soybean genome. The map was constructed from 78001 clones from five soybean BAC and BIBAC libraries representing 9.6 haploid genomes and three cultivars, and consisted of 2905 BAC/BIBAC contigs, estimated to span 1408 Mb in physical length. We evaluated the reliability of the map contigs using different contig assembly strategies, independent contig building methods, DNA marker hybridization, and different fingerprinting methods, and the results showed that the contigs were assembled properly. Furthermore, we tested the feasibility of integrating the physical map with the existing soybean composite genetic map using 388 DNA markers. The results further confirmed the nature of the ancient tetraploid origin of soybean and indicated that it is feasible to integrate the physical map with the linkage map even though greater efforts are needed. This map represents the first genome-wide, BAC/BIBAC-based physical map of the soybean genome and would provide a platform for advanced genome research of soybean and other legume species. The inclusion of BIBACs in the map would streamline the utility of the map for positional cloning of genes and QTLs, and functional analysis of soybean genomic sequences.

Characterization ofDrosophila fruitless-gal4 transgenes reveals expression in male-specificfruitless neurons and innervation of male reproductive structures

The fruitless (fru) gene acts in the central nervous system (CNS) of Drosophila melanogaster to establish male sexual behavior. Genetic dissection of the locus has shown that one of the fru gene's promoter, P1, controls the spatial and temporal expression of male-specific FruM proteins critical to determining stereotypical male sexual behavior. By using the Gal4-expression system, we show that a 16-kb fragment of the fru P1 promoter's 5' regulatory region drives the expression of Gal4 in a subset of FruM-expressing neurons within both the pupal and adult CNS. Colocalization of FruM and a Gal4-responsive reporter shows that the fru(P1)-gal4 fusion construct generates expression in both previously characterized FruM-expressing neurons as well as within cells of both the CNS and the peripheral nervous system that have not been demonstrated as FruM-expressing. Gal4-expressing neurons are shown to innervate abdominal organs directly relevant to fru function; specifically, the muscle of Lawrence (MOL) and the male internal reproductive organs. Innervations of the latter are shown to originate from identified FruM-serotonergic neurons. Furthermore, we show that the MOL neuromuscular junction is sexually dimorphic. Finally, we describe Gal4 expression in neurites innervating male reproductive structures that are hypothesized to be targets of fru function. Isolation of the regulatory sequences controlling the expression of fru in the CNS, therefore, provides a potent tool for the manipulation of FruM-expressing neurons and for understanding the cellular basis of Drosophila reproductive behavior.

Locating sequence on FPC maps and selecting a minimal tiling path

This study discusses three software tools, the first two aid in integrating sequence with an FPC physical map and the third automatically selects a minimal tiling path given genomic draft sequence and BAC end sequences. The first tool, FSD (FPC Simulated Digest), takes a sequenced clone and adds it back to the map based on a fingerprint generated by an in silico digest of the clone. This allows verification of sequenced clone positions and the integration of sequenced clones that were not originally part of the FPC map. The second tool, BSS (Blast Some Sequence), takes a query sequence and positions it on the map based on sequence associated with the clones in the map. BSS has multiple uses as follows: (1) When the query is a file of marker sequences, they can be added as electronic markers. (2) When the query is draft sequence, the results of BSS can be used to close gaps in a sequenced clone or the physical map. (3) When the query is a sequenced clone and the target is BAC end sequences, one may select the next clone for sequencing using both sequence comparison results and map location. (4) When the query is whole-genome draft sequence and the target is BAC end sequences, the results can be used to select many clones for a minimal tiling path at once. The third tool, pickMTP, automates the majority of this last usage of BSS. Results are presented using the rice FPC map, BAC end sequences, and whole-genome shotgun from Syngenta.

Software for automated analysis of DNA fingerprinting gels

Here we describe software tools for the automated detection of DNA restriction fragments resolved on agarose fingerprinting gels. We present a mathematical model for the location and shape of the restriction fragments as a function of fragment size, with model parameters determined empirically from "marker" lanes containing molecular size standards. Automated identification of restriction fragments involves several steps, including: image preprocessing, to put the data in a form consistent with a linear model; marker lane analysis, for determination of the model parameters; and data lane analysis, a procedure for detecting restriction fragment multiplets while simultaneously determining the amplitude curve that describes restriction fragment amplitude as a function of mobility. In validation experiments conducted on fingerprinted and sequenced Bacterial Artificial Chromosome (BAC) clones, sensitivity and specificity of restriction fragment identification exceeded 96% on restriction fragments ranging in size from 600 base pairs (bp) to 30,000 bp. The integrated suite of software tools, written in MATLAB and collectively called BandLeader, is in use at the BC Cancer Agency Genome Sciences Centre (GSC) and the Washington University Genome Sequencing Center, and has been provided to the Wellcome Trust Sanger Institute and the Whitehead Institute. Employed in a production mode at the GSC, BandLeader has been used to perform automated restriction fragment identification for more than 850,000 BAC clones for mouse, rat, bovine, and poplar fingerprint mapping projects.

Construction of a BAC library and generation of BAC end sequence-tagged connectors for genome sequencing of the African malaria mosquito Anopheles gambiae

A Bacterial Artificial Chromosome (BAC) genomic DNA library of Anopheles gambiae, the major human malaria vector in sub-Saharan Africa, was constructed and characterized. This library (ND-TAM) is composed of 30,720 BAC clones in eighty 384-well plates. The estimated average insert size of the library is 133 kb, with an overall genome coverage of approximately 14-fold. The ends of approximately two-thirds of the clones in the library were sequenced, yielding 32,340 pair-mate ends. A statistical analysis (G-test) of the results of PCR screening of the library indicated a random distribution of BACs in the genome, although one gap encompassing the white locus on the X-chromosome was identified. Furthermore, combined with another previously constructed BAC library (ND-1), ~2,000 BACs have been physically mapped by polytene chromosomal in situ hybridization. These BAC end pair mates and physically mapped BACs have been useful for both the assembly of a fully sequenced A. gambiae genome and for linking the assembled sequence to the three polytene chromosomes. This ND-TAM library is now publicly available at both http://www.malaria.mr4.org/mr4pages/index.html/ and http://hbz.tamu.edu/, providing a valuable resource to the mosquito research community.

Highly efficient modification of bacterial artificial chromosomes (BACs) using novel shuttle vectors containing the R6Kgamma origin of replication

Bacterial artificial chromosome (BAC) mediated transgenesis has proven to be a highly reliable way to obtain accurate transgene expression for in vivo studies of gene expression and function. A rate-limiting step in use of this technology to characterize large numbers of genes has been the process with which BACs can be modified by homologous recombination in Escherichia coli. We report here a highly efficient method for modifying BACs by using a novel set of shuttle vectors that contain the R6Kgamma origin for DNA replication, the E. coli RecA gene for recombination, and the SacB gene for negative selection. These new vectors greatly increased the ease with which one can clone the shuttle vectors, as well as screen for co-integrated and resolved clones. Furthermore, we simplify the shuttle vector cloning to one step by incorporation of a "built-in" resolution cassette for rapid removal of the unwanted vector sequences. This new system has been used to modify a dozen BACs. It is well suited for efficient production of modified BACs for use in a variety of in vivo studies.

Genomics in pure and applied entomology

Genomics is the study of the structure and function of the genome: the set of genetic information encoded in the DNA of the nucleus and organelles of an organism. It is a dynamic field that combines traditional paths of inquiry with new approaches that would have been impossible without recent technological developments. Much of the recent focus has been on obtaining the sequence of entire genomes, determining the order and organization of the genes, and developing libraries that provide immediate physical access to any desired DNA fragment. This has enabled functional studies on a genome-wide level, including analysis of the genetic basis of complex traits, quantification of global patterns of gene expression, and systematic gene disruption projects. The successful contribution of genomics to problems in applied entomology requires the cooperation of the private and public sectors to build upon the knowledge derived from the Drosophila genome and effectively develop models for other insect Orders.

Programmed cell death takes flight: genetic and genomic approaches to gene discovery in Drosophila

Programmed cell death (PCD) is an essential and wide-spread physiological process that results in the elimination of cells. Genes required to carry out this process have been identified, and many of these remain the subjects of intense investigation. Here, we describe PCD, its functions, and some of the consequences when it goes awry. We review PCD in the model system, the fruit fly, Drosophila melanogaster, with a particular emphasis on cell death gene discovery resulting from both genetics and genomics-based approaches.

Finishing a whole-genome shotgun: release 3 of the Drosophila melanogaster euchromatic genome sequence

BACKGROUND

The Drosophila melanogaster genome was the first metazoan genome to have been sequenced by the whole-genome shotgun (WGS) method. Two issues relating to this achievement were widely debated in the genomics community: how correct is the sequence with respect to base-pair (bp) accuracy and frequency of assembly errors? And, how difficult is it to bring a WGS sequence to the accepted standard for finished sequence? We are now in a position to answer these questions.

RESULTS

Our finishing process was designed to close gaps, improve sequence quality and validate the assembly. Sequence traces derived from the WGS and draft sequencing of individual bacterial artificial chromosomes (BACs) were assembled into BAC-sized segments. These segments were brought to high quality, and then joined to constitute the sequence of each chromosome arm. Overall assembly was verified by comparison to a physical map of fingerprinted BAC clones. In the current version of the 116.9 Mb euchromatic genome, called Release 3, the six euchromatic chromosome arms are represented by 13 scaffolds with a total of 37 sequence gaps. We compared Release 3 to Release 2; in autosomal regions of unique sequence, the error rate of Release 2 was one in 20,000 bp.

CONCLUSIONS

The WGS strategy can efficiently produce a high-quality sequence of a metazoan genome while generating the reagents required for sequence finishing. However, the initial method of repeat assembly was flawed. The sequence we report here, Release 3, is a reliable resource for molecular genetic experimentation and computational analysis.

Heterochromatic sequences in a Drosophila whole-genome shotgun assembly

BACKGROUND

Most eukaryotic genomes include a substantial repeat-rich fraction termed heterochromatin, which is concentrated in centric and telomeric regions. The repetitive nature of heterochromatic sequence makes it difficult to assemble and analyze. To better understand the heterochromatic component of the Drosophila melanogaster genome, we characterized and annotated portions of a whole-genome shotgun sequence assembly.

RESULTS

WGS3, an improved whole-genome shotgun assembly, includes 20.7 Mb of draft-quality sequence not represented in the Release 3 sequence spanning the euchromatin. We annotated this sequence using the methods employed in the re-annotation of the Release 3 euchromatic sequence. This analysis predicted 297 protein-coding genes and six non-protein-coding genes, including known heterochromatic genes, and regions of similarity to known transposable elements. Bacterial artificial chromosome (BAC)-based fluorescence in situ hybridization analysis was used to correlate the genomic sequence with the cytogenetic map in order to refine the genomic definition of the centric heterochromatin; on the basis of our cytological definition, the annotated Release 3 euchromatic sequence extends into the centric heterochromatin on each chromosome arm.

CONCLUSIONS

Whole-genome shotgun assembly produced a reliable draft-quality sequence of a significant part of the Drosophila heterochromatin. Annotation of this sequence defined the intron-exon structures of 30 known protein-coding genes and 267 protein-coding gene models. The cytogenetic mapping suggests that an additional 150 predicted genes are located in heterochromatin at the base of the Release 3 euchromatic sequence. Our analysis suggests strategies for improving the sequence and annotation of the heterochromatic portions of the Drosophila and other complex genomes.

Fruits of human genome project and private venture, and their impact on life science

A small knowledge base was created by organizing the Human Genome Project (HGP) and its related issues in "Science" magazines between 1996 and 2000. This base revealed the stunning achievement of HGP and a private venture and its impact on today's biology and life science. In the mid-1990, they encouraged the development of advanced high throughput automated DNA sequencers and the technologies that can analyse all genes at once in a systematic fashion. Using these technologies, they completed the genome sequence of human and various other organisms. These fruits opened the door to comparative genomics, functional genomics, the interdisprinary field between computer and biology, and proteomics. They have caused a shift in biological investigation from studying single genes or proteins to studying all genes or proteins at once, and causing revolutional changes in traditional biology, drug discovery and therapy. They have expanded the range of potential drug targets and have facilitated a shift in drug discovery programs toward rational target-based strategies. They have spawned pharmacogenomics that could give rise to a new generation of highly effective drugs that treat causes, not just symptoms. They should also cause a migration from the traditional medications that are safe and effective for every members of the population to personalized medicine and personalized therapy.

Bacterial artificial chromosome-based physical map of the rice genome constructed by restriction fingerprint analysis

Genome-wide physical mapping with bacteria-based large-insert clones (e.g., BACs, PACs, and PBCs) promises to revolutionize genomics of large, complex genomes. To accelerate rice and other grass species genome research, we developed a genome-wide BAC-based map of the rice genome. The map consists of 298 BAC contigs and covers 419 Mb of the 430-Mb rice genome. Subsequent analysis indicated that the contigs constituting the map are accurate and reliable. Particularly important to proficiency were (1) a high-resolution, high-throughput DNA sequencing gel-based electrophoretic method for BAC fingerprinting, (2) the use of several complementary large-insert BAC libraries, and (3) computer-aided contig assembly. It has been demonstrated that the fingerprinting method is not significantly influenced by repeated sequences, genome size, and genome complexity. Use of several complementary libraries developed with different restriction enzymes minimized the "gaps" in the physical map. In contrast to previous estimates, a clonal coverage of 6.0-8.0 genome equivalents seems to be sufficient for development of a genome-wide physical map of approximately 95% genome coverage. This study indicates that genome-wide BAC-based physical maps can be developed quickly and economically for a variety of plant and animal species by restriction fingerprint analysis via DNA sequencing gel-based electrophoresis.

Strategies for the systematic sequencing of complex genomes

Recent spectacular advances in the technologies and strategies for DNA sequencing have profoundly accelerated the detailed analysis of genomes from myriad organisms. The past few years alone have seen the publication of near-complete or draft versions of the genome sequence of several well-studied, multicellular organisms - most notably, the human. As well as providing data of fundamental biological significance, these landmark accomplishments have yielded important strategic insights that are guiding current and future genome-sequencing projects.

Analysis of conserved noncoding DNA in Drosophila reveals similar constraints in intergenic and intronic sequences

Comparative genomic approaches to gene and cis-regulatory prediction are based on the principle that differential DNA sequence conservation reflects variation in functional constraint. Using this principle, we analyze noncoding sequence conservation in Drosophila for 40 loci with known or suspected cis-regulatory function encompassing >100 kb of DNA. We estimate the fraction of noncoding DNA conserved in both intergenic and intronic regions and describe the length distribution of ungapped conserved noncoding blocks. On average, 22%-26% of noncoding sequences surveyed are conserved in Drosophila, with median block length approximately 19 bp. We show that point substitution in conserved noncoding blocks exhibits transition bias as well as lineage effects in base composition, and occurs more than an order of magnitude more frequently than insertion/deletion (indel) substitution. Overall, patterns of noncoding DNA structure and evolution differ remarkably little between intergenic and intronic conserved blocks, suggesting that the effects of transcription per se contribute minimally to the constraints operating on these sequences. The results of this study have implications for the development of alignment and prediction algorithms specific to noncoding DNA, as well as for models of cis-regulatory DNA sequence evolution.

From first base: the sequence of the tip of the X chromosome of Drosophila melanogaster, a comparison of two sequencing strategies

We present the sequence of a contiguous 2.63 Mb of DNA extending from the tip of the X chromosome of Drosophila melanogaster. Within this sequence, we predict 277 protein coding genes, of which 94 had been sequenced already in the course of studying the biology of their gene products, and examples of 12 different transposable elements. We show that an interval between bands 3A2 and 3C2, believed in the 1970s to show a correlation between the number of bands on the polytene chromosomes and the 20 genes identified by conventional genetics, is predicted to contain 45 genes from its DNA sequence. We have determined the insertion sites of P-elements from 111 mutant lines, about half of which are in a position likely to affect the expression of novel predicted genes, thus representing a resource for subsequent functional genomic analysis. We compare the European Drosophila Genome Project sequence with the corresponding part of the independently assembled and annotated Joint Sequence determined through "shotgun" sequencing. Discounting differences in the distribution of known transposable elements between the strains sequenced in the two projects, we detected three major sequence differences, two of which are probably explained by errors in assembly; the origin of the third major difference is unclear. In addition there are eight sequence gaps within the Joint Sequence. At least six of these eight gaps are likely to be sites of transposable elements; the other two are complex. Of the 275 genes in common to both projects, 60% are identical within 1% of their predicted amino-acid sequence and 31% show minor differences such as in choice of translation initiation or termination codons; the remaining 9% show major differences in interpretation.

A bacterial artificial chromosome library for sequencing the complete human genome

A 30-fold redundant human bacterial artificial chromosome (BAC) library with a large average insert size (178 kb) has been constructed to provide the intermediate substrate for the international genome sequencing effort. The DNA was obtained from a single anonymous volunteer, whose identity was protected through a double-blind donor selection protocol. DNA fragments were generated by partial digestion with EcoRI (library segments 1--4: 24-fold) and MboI (segment 5: sixfold) and cloned into the pBACe3.6 and pTARBAC1 vectors, respectively. The quality of the library was assessed by extensive analysis of 169 clones for rearrangements and artifacts. Eighteen BACs (11%) revealed minor insert rearrangements, and none was chimeric. This BAC library, designated as "RPCI-11," has been used widely as the central resource for insert-end sequencing, clone fingerprinting, high-throughput sequence analysis and as a source of mapped clones for diagnostic and functional studies.

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

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

Identification of Waldo-A and Waldo-B, two closely related non-LTR retrotransposons in Drosophila

We have identified two novel, closely related subfamilies of non-long-terminal-repeat (non-LTR) retrotransposons in Drosophila melanogaster, the Waldo-A and Waldo-B subfamilies, that are in the same lineage as site-specific LTR retrotransposons of the R1 clade. Both contain potentially active copies with two large open reading frames, having coding capacities for a nucleoprotein as well as endonuclease and reverse transcriptase activities. Many copies are truncated at the 5' end, and most are surrounded by target site duplications of variable lengths. Elements of both subfamilies have a nonrandom distribution in the genome, often being inserted within or very close to (CA)(n) arrays. At the DNA level, the longest elements of Waldo-A and Waldo-B are 69% identical on their entire length, except for the 5' untranslated regions, which have a mosaic organization, suggesting that one arose from the other following new promoter acquisition. This event occurred before the speciation of the D. melanogaster subgroup of species, since both Waldo-A and Waldo-B coexist in other species of this subgroup.

Characterization of a MEN1 ortholog from Drosophila melanogaster

Multiple endocrine neoplasia type 1 (MEN1) is a familial cancer syndrome characterized by tumors of the parathyroid, entero-pancreatic neuroendocrine and pituitary tissues and caused by inactivating mutations in the MEN1 gene. Menin, the 610-amino acid nuclear protein encoded by MEN1, binds to the transcription factor JunD and can repress JunD-induced transcription. We report here the identification of a MEN1 ortholog in Drosophila melanogaster, Menin1, that encodes a 763 amino acid protein sharing 46% identity with human menin. Additionally, 69% of the missense mutations and in-frame deletions reported in MEN1 patients appear in amino acid residues that are identical in the Drosophila and human protein, suggesting the importance of the conserved regions. Drosophila Menin1 gene transcripts use alternative polyadenylation sites resulting in 4.3 and 5-kb messages. The 4.3-kb transcript appears to be largely maternal, while the 5-kb transcript appears mainly zygotic. The binding of Drosophila menin to human JunD or Drosophila Jun could not be demonstrated by the yeast two-hybrid analysis. The identification of the MEN1 ortholog from Drosophila melanogaster will provide an opportunity to utilize Drosophila genetics to enhance our understanding of the function of human menin.

The Drosophila genome

The past year has been a spectacular one for Drosophila research. The sequencing and annotation of the Drosophila melanogaster genome has allowed a comprehensive analysis of the first three eukaryotes to be sequenced-yeast, worm and fly-including an analysis of the fly's influences as a model for the study of human disease. This year has also seen the initiation of a full-length cDNA sequencing project and the first analysis of Drosophila development using high-density DNA microarrays containing several thousand Drosophila genes. For the first time homologous recombination has been demonstrated in flies and targeted gene disruptions may not be far off.

Contigs built with fingerprints, markers, and FPC V4.7

Contigs have been assembled, and over 2800 clones selected for sequencing for human chromosomes 9, 10 and 13. Using the FPC (FingerPrinted Contig) software, the contigs are assembled with markers and complete digest fingerprints, and the contigs are ordered and localised by a global framework. Publicly available resources have been used, such as, the 1998 International Gene Map for the framework and the GSC Human BAC fingerprint database for the majority of the fingerprints. Additional markers and fingerprints are generated in-house to supplement this data. To support the scale up of building maps, FPC V4.7 has been extended to use markers with the fingerprints for assembly of contigs, new clones and markers can be automatically added to existing contigs, and poorly assembled contigs are marked accordingly. To test the automatic assembly, a simulated complete digest of 110 Mb of concatenated human sequence was used to create datasets with varying coverage, length of clones, and types of error. When no error was introduced and a tolerance of 7 was used in assembly, the largest contig with no false positive overlaps has 9534 clones with 37 out-of-order clones, that is, the starting coordinates of adjacent clones are in the wrong order. This paper describes the new features in FPC, the scenario for building the maps of chromosomes 9, 10 and 13, and the results from the simulation.

Contigs built with fingerprints, markers, and FPC V4.7

Contigs have been assembled, and over 2800 clones selected for sequencing for human chromosomes 9, 10 and 13. Using the FPC (FingerPrinted Contig) software, the contigs are assembled with markers and complete digest fingerprints, and the contigs are ordered and localised by a global framework. Publicly available resources have been used, such as, the 1998 International Gene Map for the framework and the GSC Human BAC fingerprint database for the majority of the fingerprints. Additional markers and fingerprints are generated in-house to supplement this data. To support the scale up of building maps, FPC V4.7 has been extended to use markers with the fingerprints for assembly of contigs, new clones and markers can be automatically added to existing contigs, and poorly assembled contigs are marked accordingly. To test the automatic assembly, a simulated complete digest of 110 Mb of concatenated human sequence was used to create datasets with varying coverage, length of clones, and types of error. When no error was introduced and a tolerance of 7 was used in assembly, the largest contig with no false positive overlaps has 9534 clones with 37 out-of-order clones, that is, the starting coordinates of adjacent clones are in the wrong order. This paper describes the new features in FPC, the scenario for building the maps of chromosomes 9, 10 and 13, and the results from the simulation.

A physical map of the polytenized region (101EF-102F) of chromosome 4 in Drosophila melanogaster

Chromosome 4, the smallest autosome ( approximately 5 Mb in length) in Drosophila melanogaster contains two major regions. The centromeric domain ( approximately 4 Mb) is heterochromatic and consists primarily of short, satellite repeats. The remaining approximately 1.2 Mb, which constitutes the banded region (101E-102F) on salivary gland polytene chromosomes and contains the identified genes, is the region mapped in this study. Chromosome walking was hindered by the abundance of moderately repeated sequences dispersed along the chromosome, so we used many entry points to recover overlapping cosmid and BAC clones. In situ hybridization of probes from the two ends of the map to polytene chromosomes confirmed that the cloned region had spanned the 101E-102F interval. Our BAC clones comprised three contigs; one gap was positioned distally in 102EF and the other was located proximally at 102B. Twenty-three genes, representing about half of our revised estimate of the total number of genes on chromosome 4, were positioned on the BAC contigs. A minimal tiling set of the clones we have mapped will facilitate both the assembly of the DNA sequence of the chromosome and a functional analysis of its genes.

From sequence to chromosome: the tip of the X chromosome of D. melanogaster

One of the rewards of having a Drosophila melanogaster whole-genome sequence will be the potential to understand the molecular bases for structural features of chromosomes that have been a long-standing puzzle. Analysis of 2.6 megabases of sequence from the tip of the X chromosome of Drosophila identifies 273 genes. Cloned DNAs from the characteristic bulbous structure at the tip of the X chromosome in the region of the broad complex display an unusual pattern of in situ hybridization. Sequence analysis revealed that this region comprises 154 kilobases of DNA flanked by 1.2-kilobases of inverted repeats, each composed of a 350-base pair satellite related element. Thus, some aspects of chromosome structure appear to be revealed directly within the DNA sequence itself.

The genome sequence of Drosophila melanogaster

The fly Drosophila melanogaster is one of the most intensively studied organisms in biology and serves as a model system for the investigation of many developmental and cellular processes common to higher eukaryotes, including humans. We have determined the nucleotide sequence of nearly all of the approximately 120-megabase euchromatic portion of the Drosophila genome using a whole-genome shotgun sequencing strategy supported by extensive clone-based sequence and a high-quality bacterial artificial chromosome physical map. Efforts are under way to close the remaining gaps; however, the sequence is of sufficient accuracy and contiguity to be declared substantially complete and to support an initial analysis of genome structure and preliminary gene annotation and interpretation. The genome encodes approximately 13,600 genes, somewhat fewer than the smaller Caenorhabditis elegans genome, but with comparable functional diversity.

The Drosophila genome sequence: implications for biology and medicine

The 120-megabase euchromatic portion of the Drosophila melanogaster genome has been sequenced. Because the genome is compact and many genetic tools are available, and because fly cell biology and development have much in common with mammals, this sequence may be the Rosetta stone for deciphering the human genome.

A whole-genome assembly of Drosophila

We report on the quality of a whole-genome assembly of Drosophila melanogaster and the nature of the computer algorithms that accomplished it. Three independent external data sources essentially agree with and support the assembly's sequence and ordering of contigs across the euchromatic portion of the genome. In addition, there are isolated contigs that we believe represent nonrepetitive pockets within the heterochromatin of the centromeres. Comparison with a previously sequenced 2.9- megabase region indicates that sequencing accuracy within nonrepetitive segments is greater than 99. 99% without manual curation. As such, this initial reconstruction of the Drosophila sequence should be of substantial value to the scientific community.