Skip to main content

The generation of chromosomal deletions to provide extensive coverage and subdivision of the Drosophila melanogaster genome

Abstract

Background

Chromosomal deletions are used extensively in Drosophila melanogaster genetics research. Deletion mapping is the primary method used for fine-scale gene localization. Effective and efficient deletion mapping requires both extensive genomic coverage and a high density of molecularly defined breakpoints across the genome.

Results

A large-scale resource development project at the Bloomington Drosophila Stock Center has improved the choice of deletions beyond that provided by previous projects. FLP-mediated recombination between FRT-bearing transposon insertions was used to generate deletions, because it is efficient and provides single-nucleotide resolution in planning deletion screens. The 793 deletions generated pushed coverage of the euchromatic genome to 98.4%. Gaps in coverage contain haplolethal and haplosterile genes, but the sizes of these gaps were minimized by flanking these genes as closely as possible with deletions. In improving coverage, a complete inventory of haplolethal and haplosterile genes was generated and extensive information on other haploinsufficient genes was compiled. To aid mapping experiments, a subset of deletions was organized into a Deficiency Kit to provide maximal coverage efficiently. To improve the resolution of deletion mapping, screens were planned to distribute deletion breakpoints evenly across the genome. The median chromosomal interval between breakpoints now contains only nine genes and 377 intervals contain only single genes.

Conclusions

Drosophila melanogaster now has the most extensive genomic deletion coverage and breakpoint subdivision as well as the most comprehensive inventory of haploinsufficient genes of any multicellular organism. The improved selection of chromosomal deletion strains will be useful to nearly all Drosophila researchers.

Background

Chromosomal deletions are important to experimental genetic analysis in two fundamental ways. First, deletions fail to complement loss-of-function mutations in genes located in the chromosomal region of the deletion. This noncomplementation is the basis for using deletions to map mutations to specific chromosomal regions, to screen for new mutations in closely linked sets of genes and to assess the allelic strengths of new mutations. Second, heterozygous deletions can enhance or suppress mutant phenotypes. Reducing the copy number of one gene can modify the phenotype caused by abnormal expression of another gene involved in the same biological process. Modifier screens are a powerful way to identify suites of genes involved in related genetic pathways. Drosophila melanogaster tolerates heterozygous deletions of large numbers of genes quite well and its chromosomes are easy to manipulate in vivo; consequently, deletions are used more extensively in genetic experiments with D. melanogaster than other organisms. Indeed, fine-scale mapping approaches based on meiotic recombination such as SNP mapping are not used widely in fly research, because they are less convenient and precise. Consequently, chromosomal deletion stocks are the most heavily used class of stocks distributed by the Bloomington Drosophila Stock Center (BDSC). The demand for a better selection of deletions led us at the BDSC to undertake the large-scale research resource development project described here.

Many methods have been developed for isolating deletions in Drosophila. Early methods employed chemical mutagens or irradiation to produce random chromosomal breakpoints. Deletions in specific chromosomal regions were identified by their failure to complement a recessive mutation, by removal of a gain-of-function mutation, or from dominant effects of reduced gene copy number. Obtaining a deletion of a particular size in a desired location was difficult because the mutageneses were inefficient and the positions of breakpoints could not be controlled. Later methods relied on transposons to target deletions to particular chromosomal regions. Deletions flanking P and Minos elements result from defective transposition events when attempts are made to remobilize the transposons and deletions flanking Hobo elements result from recombination between Hobo copies following local transpositions [1–6]. While one to two orders of magnitude more efficient than chemical or irradiation screens, transposon remobilization screens still have the disadvantage of breakpoint unpredictability. The positions of breakpoints are typically confined to a small chromosomal region, but they cannot be directed with single-nucleotide certainty. Furthermore, mechanistic constraints make it difficult to obtain large deletions in most remobilization screens.

A more recent method for generating deletions employs FLP recombinase to catalyze recombination between FRT sequence elements carried on transposable elements [7]. It is one to two orders of magnitude more efficient than transposon remobilization, the positions of breakpoints can be predicted at the nucleotide level and there are no inherent limitations on deletion size other than aneuploidy effects. This method has been adopted enthusiastically by Drosophila geneticists because it has made the isolation of large, multigene deletions and small, single-gene knockouts relatively easy.

There have been three large-scale projects generating deletions using FLP-FRT technology. The DrosDel Project generated a large collection of insertions of the FRT-bearing transposons P{RS3} and P{RS5} [8] and used them in screens that provided 357 deletion stocks now in the BDSC collection [9]. Exelixis, Inc. generated a larger insertion collection of the FRT-bearing constructs P{XP}, PBac{RB} and PBac{WH} [10] and isolated 433 deletions now distributed by the BDSC [5]. These two deletion collections were generated independently and the deletions from these projects being distributed by the BDSC together provide 78% coverage of Drosophila euchromatin.

The third deletion project was undertaken at the BDSC and it is the subject of this paper. We began generating deletions before FLP-FRT technology for deletion screening was generally available. At the time, genomic coverage of approximately 71% was provided by chemical- or irradiation-induced deletions. We were able to provide novel coverage of 5 to 7% of the genome using a P element remobilization approach [5], but we retooled our project when the Exelixis FRT-bearing insertions were placed into the public domain [10, 11]. As we will show, our efforts using the improved technology have resulted in total euchromatic coverage of 98.4%.

To provide such extensive deletion coverage, it was necessary to identify all genes needed in two copies for normal viability and fertility. These haploinsufficient genes are the only barriers to complete deletion coverage in the absence of compensating chromosomal duplications and our work provides the first comprehensive inventory of these genes in Drosophila melanogaster. Most remaining gaps in deletion coverage contain a haploinsufficient gene and, because we flanked these genes as closely as possible with pairs of deletions, the gaps are all quite small.

To make it easier for geneticists to take advantage of this broad deletion coverage, the BDSC distributes the 'Deficiency Kit', a selected set of deletions providing maximal genomic coverage with a minimal number of stocks. Here we will describe a new, improved Deficiency Kit composed primarily of molecularly defined deletions from the three deletion projects. It replaces the original Deficiency Kit that consisted of older deletions characterized only at the level of chromosome banding.

While the fraction of the genome covered by deletions in aggregate determines their effectiveness in initially localizing genes, the density of deletion breakpoints determines the resolution of deletion mapping. The more finely the genome is subdivided by breakpoints, the easier it is to map a mutation or genetic modifier to a specific transcription unit. For this reason, we have also worked to improve the distribution of breakpoints across the genome. Because of our efforts, the median interval between breakpoints is now only nine genes and nearly 400 single-gene intervals have been defined.

In brief, we accomplished three goals. First, we improved deletion coverage of the D. melanogaster genome substantially and, in the process, provided a complete catalog of haplolethal and haplosterile genes. Second, we updated the Bloomington Deficiency Kit by replacing most chemical- or irradiation-induced deletions with molecularly defined deletions isolated from FLP-FRT screens. Finally, we improved the resolution of deletion mapping by increasing the density of molecularly defined breakpoints across the genome. The result is a collection of deletions at the BDSC for the rapid and efficient mapping of mutations and modifier loci that is unparalleled among multicellular organisms.

Results

Generating molecularly defined deletions by FLP-mediated recombination

The system for generating deletions devised at Exelixis involves the use of three FRT-bearing transposon constructs (Figure 1a) [5, 10]. Each construct carries the miniwhite marker and at least one FRT sequence. Deletions can be recovered by combining pairs of insertions in trans in the presence of an inducible genomic source of FLP recombinase. FLP-mediated recombination generates a deletion if the FRT sequences are arranged in the same relative orientation. If the miniwhite markers lie internal to the FRT sequences, deletion-bearing progeny will have white eyes from elimination of both copies of miniwhite upon recombination (Figure 1b). If they lie external to the FRTs, deletion-bearing progeny often have darker eyes than flies in either single-FRT progenitor strain from the presence of two copies of miniwhite (Figure 1c). We generated 87% of our deletions using construct combinations eliminating miniwhite markers; the remaining deletions were isolated in progeny bearing two copies of miniwhite. Pairs of insertions can also be used to recover deletions with no net gain or loss of miniwhite markers (Figure 1d), but such screens were not attempted in this project.

Figure 1
figure 1

Generating molecularly defined deletions using Exelixis FRT -bearing transposon insertions. (a) The structure of the P{XP}, PBac{RB}and PBac{WH} constructs. The miniwhite-marked P element or piggyBac (PBac) constructs carry one or two FRT sequences with the indicated orientations. (miniwhite is a version of the white gene engineered for compact size.) P{XP} and PBac{WH} constructs also carry UAS sequences oriented to allow GAL4-induced expression of genes near the genomic insertion sites of the constructs. One UAS sequence in P{XP} can be removed by FLP-mediated recombination and it is likely that this cassette is absent from most deletion chromosomes, though we did not assay for it in our deletion chromosomes. (b-d) Simplified diagrams of FLP-mediated recombination events generating deletion chromosomes with different miniwhite copy numbers. (b) Our most frequently used screening strategy, where deletions are identified based on loss of miniwhite and the resulting white eye color. (c) The alternative strategy of identifying deletions based on increased miniwhite copy number relative to progenitor chromosomes and the resulting darker eye color. (d) Deletions can be recovered without a decrease or increase in miniwhite copy number, though we did not undertake such screens.

Approximately 16,500 FRT-bearing transposon insertions on the X chromosome and autosomes were generated and placed into public distribution by Exelixis [10, 11]. The density of insertions afforded by such a large collection gave us flexibility in designing deletions. As we will describe below, convenient insertions were usually available so that we could flank haploinsufficient genes closely with pairs of deletions in our efforts to maximize deletion coverage and improve the distribution of breakpoints. In some chromosomal regions, however, the scarcity of insertions limited our efforts.

We isolated 793 deletions. They are described in detail in Additional file 1 and their breakpoint positions are shown relative to annotated genes in Additional file 2. Because we were primarily concerned with breakpoint placement, we did not standardize the size of the deletions. We did not generate deletions with heterochromatic breakpoints because most heterochromatic P{XP}, PBac{WH} and PBac{RB} insertions have not been localized unambiguously. In the screens that were successful, one deletion was recovered for every approximately 50 progeny screened with a range of roughly 1 in 10 to 1 in 1,000. We retained only one deletion from each screen for the BDSC collection.

We confirmed the existence of deletions by demonstrating that genes located between the progenitor insertions were missing. We usually designed deletions so they encompassed at least one gene with an existing loss-of-function mutation to use in complementation tests. Only two deletions designed early in the project were small enough to avoid removing vital genes. We also used preexisting deletions in complementation tests with new deletions if the progenitor FRT-bearing insertions were not themselves homozygous lethal to demonstrate that at least one vital gene was deleted. If an appropriate complementation test was not available, we verified the presence of a deletion by polytene cytology or by demonstrating that a particular internal DNA fragment could not be PCR amplified from the deletion chromosome. Eighteen new deletions showed dominant phenotypes expected for deleting known haploinsufficient genes (Additional file 1).

We also verified that the expected deletion event occurred using the 'hybrid PCR' method [5]. In brief, PCR primers within each FRT-bearing progenitor insertion can amplify a fragment only if they are juxtaposed by FLP-mediated recombination. We tested 93% of the deletions lacking miniwhite markers by this method and all gave the expected result (Additional file 1). We used a positive hybrid PCR test as the sole criterion for verifying only six deletions; all other deletions were also confirmed on the basis of deleted genes or sequences.

The Exelixis insertions were generated in a uniform genetic background [10] and we maintained this background in most deletion screens as described in Parks et al.[5]. Noninsertion chromosomes specific to our screens were substituted into the standard background prior to use. Minimizing genetic variability in the deletion stocks makes them more valuable in the analysis of background-sensitive phenotypes such as behavioral traits. It was necessary to establish 27 deletions in stock with nonstandard balancers due to dominant effects of the deletions on viability or fertility or due to noncomplementation with mutations on standard balancers (Additional file 1).

Deletion stocks were placed into public distribution as soon as the deletions were verified; consequently, many deletions have been in widespread use for years and have contributed to published research (for examples, [12–18]). Details of each deletion screen have been documented in FlyBase and references to the individual reports are given in Additional file 1. FlyBase developed a graphical interface for viewing molecularly defined deletions during our project [19] and the BDSC website lists available deletions [20].

We conducted 708 screens that produced no deletions and we were often able to determine the cause of the failure. First, the rate of recovering deletions was often low. We were able to generate most deletions with modest-sized screens, but some deletions required larger screens. Four percent of deletions were recovered only after additional rounds of screening. Five percent of the screens failed even though both progenitor insertions had been used in successful screens, there was no haploinsufficient gene in the region and the deletion would not have been particularly large, suggesting small screen size as the most likely explanation for failure. Second, the insertions in many stocks did not map as described. We saw evidence that insertions were mismapped in 12% of the unsuccessful screens. (All eight of the fourth chromosome deletion screens were in this category, suggesting that a high proportion of the fourth chromosome insertion stocks are not as represented. We abandoned plans for additional screens without validating the remaining putative fourth chromosome insertion stocks.) Third, in 6% of the failed screens, putative deletions could not be established in stock due to lethality and/or sterility, even though no haploinsufficient genes existed in the region. We suspect that the abnormal juxtaposition of genic regions upon deletion results in deleterious neomorphic phenotypes more often than is generally appreciated. Fourth, 5% of unsuccessful screens failed because flies in the intermediate crosses were weak, sterile or dead. We attribute these problems to noncomplementation and dominant interactions between mutations in the stocks. Fifth, 10% of unsuccessful screens failed because we intentionally deleted haploinsufficient genes. (We will discuss these experiments below.) The remaining screens probably failed from small screen sizes or mismapped insertions. They all involved at least one insertion stock that had not been used successfully in another screen. From these experiences, we recommend that anyone planning FLP-FRT screens verify the positions of progenitor insertions, plan screens larger than our routine screen sizes and, if possible, attempt to create deletions of a region using more than one combination of progenitor insertions.

False positives were recovered in many screens. For example, white-eyed progeny carrying no deletion were recovered in screens designed to eliminate the miniwhite markers from both progenitor insertions upon FLP-mediated recombination (Figure 1b). While false positives did not prevent the recovery of deletions, they increased the workload in the verification step. At least 9% of screens successfully producing deletions also produced false positives. Inexplicably, X chromosome screens were much more likely to give false positives than autosomal screens (49% versus 4%). At least 23% of the screens producing no deletions gave false positives, again, with a disproportionate number involving the X chromosome rather than an autosome (56% versus 15%). False positives were also seen in Exelixis and DrosDel screens (A Parks, personal communication) [9]. Consequently, we recommend testing multiple putative deletions in all screens. Though the false positives have not been characterized in detail, we surmise they originate from the repair of DNA damage associated with abnormal FLP-mediated recombination events.

Improving deletion coverage of the genome

The first goal of our project was to improve deletion coverage of the genome. As summarized in Table 1 and shown in detail in Additional file 2, 98.4% coverage of euchromatic genes is now available using deletions distributed by the BDSC regardless of their origin. The deletions we generated (the 'BSC' deletions) encompass 81.5% of euchromatic genes as compared to 65.2% for the DrosDel deletions at Bloomington and 47.4% for the Exelixis deletions. Together, the molecularly defined deletions from the three large-scale projects provide 94.9% coverage with 17.9% of euchromatic genes removed exclusively by BSC deletions. The remaining genomic coverage is provided by chemical- or irradiation-induced deletions (3.4%) and FRT-derived deletions from individual investigators (0.1%).

Table 1 Percent coverage of euchromatic genes by deletions

The single biological impediment to complete deletion coverage is the existence of genes that must be present in more than one copy for viability and fertility. Deletions of these haplolethal and haplosterile genes cannot be recovered in the absence of a compensating chromosomal duplication. (Obviously, haploinsufficiency pertains to X-linked genes only in females where two gene copies are normally present.) Duplications had been isolated for some haplolethal and haplosterile genes before we began our work, but generating new duplications in order to screen for deletions was beyond the scope of this project. Instead, we generated pairs of deletions to flank haplolethal and haplosterile genes as closely as possible as a way to maximize deletion coverage. This strategy allowed us to improve coverage substantially, but it precluded 100% coverage. The sizes of gaps in coverage were determined by the availability of FRT-bearing transgene insertions; consequently, some genes closely linked to haplolethal and haplosterile genes also remain undeleted. We chose not to use preexisting duplications to screen for deletions of haploinsufficient genes because we wanted to preserve a consistent genetic background in our deletion stocks. We were also aware of the fact that deletions of haploinsufficient genes have limited experimental value: the phenotypes of haploinsufficient genes often confound the interpretation of complementation tests and they enhance so many other mutant phenotypes that they are generally viewed as nuisances in modifier screens.

If we consider the distribution of deletion breakpoints for the DrosDel, Exelixis and BSC deletions, we see there is a total of 81 gaps in deletion coverage across the entire euchromatic genome (Additional file 2). Fifty-five gaps contain haploinsufficient genes, five gaps (including three with haploinsufficient genes) lie adjacent to centric heterochromatin where screening was not possible, three gaps (including one with a haploinsufficient and one adjacent centric heterochromatin) lie on the fourth chromosome where we recovered no deletions and 23 gaps are not covered despite the absence of haploinsufficient genes. The median size of the gaps is only 7 genes and the gaps range in size from 1 to 36 genes (1 to 27 if a large fourth chromosome gap is disregarded). We made multiple attempts to provide deletion coverage for most gaps lacking haplolethal or haplosterile genes, but we were unsuccessful for the reasons outlined above.

Fortunately, many gaps in coverage between BSC, Exelixis and DrosDel deletions are covered by older deletions or FRT-derived deletions from individual investigators. We have no clear explanation for the failure of FLP-FRT screens in regions where older deletions provided coverage. Although the haplosterility or haplolethality of some genes might depend on genetic background, it is unlikely this is the case for a large number of genes. Screens with different progenitor insertions to the ones we used might prove successful. When deletions from all sources are considered - including those maintained in stock with a duplication - only 35 gaps in coverage remain. Thirty-two gaps contain haploinsufficient genes, three gaps (including one with a haploinsufficient gene) lie adjacent to centric heterochromatin, and two gaps (including one adjacent to centric heterochromatin) lie on the fourth chromosome. The median gap size is 7 genes with a range of 1 to 18 genes. We verified that the older deletions span gaps with complementation tests and PCR amplifications. In many cases, these verification experiments also allowed us to map deletion breakpoints of the older deletions more precisely (Additional file 2). In addition, the breakpoints of three gap-covering deletions (Df(2L)C144, Df(2R)CX1 and Df(3L)ZN47) were mapped with comparative genomic hybridization microarrays at the same time breakpoints of X chromosome duplications from another BDSC resource project were localized [21].

Identifying haploinsufficient genes

Although it was necessary to identify haplolethal and haplosterile genes in order to improve deletion coverage, most were poorly mapped and few had been associated with annotated genes when we began our project. Consequently, the process of planning deletion screens was intertwined with the task of identifying haploinsufficient genes. When there were no hints about the molecular identity of a haploinsufficient gene, we faced the tedious process of iterative deletion screening to narrow its location.

Fortunately, most haploinsufficient genes in Drosophila encode protein components of cytoplasmic ribosomes (Rp genes) or translation initiation factors (eIF genes) and are associated with the Minute syndrome, a characteristic set of phenotypes including short, thin bristles and slower development [22]. All Minute genes are associated with some degree of haplolethality and haplosterility, but the phenotypes of many Minute genes are not severe enough to prevent the recovery of deletions in the absence of duplications. We assured that all Minute genes were flanked as closely as possible with deletions and we attempted to delete 34 Minute genes to assess their viability and fertility. Deletions for 7 Minute genes could be established in simple heterozygous stocks, but deletions for the remaining 27 could not be recovered or heterozygous stocks were too weak to maintain long term. We did not attempt to delete Minute genes where evidence already existed for strong lethality or sterility effects. Table 2 shows Minute genes with viability or fertility effects strong enough to prevent the recovery of deletions in the absence of a duplication; Table 3 shows Minute genes with weaker effects. Because Minute phenotypes lie on a continuum and vary by genetic background, genes with intermediate effects might be classified differently in other experiments. Additional files 2 and 3 provide details of deletion coverage near Minute genes. This inventory updates the number of Minute loci given in Marygold et al.[22] from 65 to 66 with the reclassification of RpL23.

Table 2 Haplolethal and haplosterile genes
Table 3 Minute genes not associated with strong viability and fertility effects

Nine haplolethal or haplosterile loci do not encode ribosome-associated proteins (Table 2; Additional files 2 and 3). Four of the non-Minute genes (Fs(1)10A, Hdl, the unnamed haplolethal in 67DE and Ms(3)88C) have not been mapped to single transcription units, but the DrosDel Project localized Ms(3)88C to a small region by successive rounds of screens [9] and we localized Fs(1)10A, Hdl and the 67DE haplolethal similarly. The Tpl locus corresponds to a repetitive gene cluster [23]. We defined a small region within the Tpl cluster that could not be further narrowed by deleting from the proximal and distal sides using the available FRT insertions, but it is possible that viability simply requires a minimal number of gene repeats. Thirty-five of our failed screens were devoted to mapping these non-Minute haploinsufficient loci.

Tables 2 and 3 present comprehensive inventories of haplolethal, haplosterile and Minute genes in D. melanogaster, but we also compiled a list of haploinsufficient genes associated with developmental or cellular phenotypes that do not affect viability or fertility enough to prevent the recovery of deletions (Table 4; Additional file 2). All were removed by molecularly defined deletions, but we made no attempt to flank these genes closely with deletions. Only the Hup locus has not been mapped to a transcription unit. We did not identify new haploinsufficient genes with obvious external morphological phenotypes and few are likely to exist. Additional haploinsufficient genes with cellular phenotypes may have been described in the literature and others will probably be identified in the future, so Table 4 is undoubtedly incomplete. We do not know how many haploinsufficient genes might exist with moderate viability or fertility effects similar to the Minute genes in Table 3 yet have no obvious external morphological effects.

Table 4 Haploinsufficient loci with developmental or cellular phenotypes not associated with strong haplolethality or haplosterility

The non-Minute haploinsufficient genes listed in Tables 2 and 4 implicate no single biological pathway or process, though there are a few discernable themes. The largest subset encodes the muscle components actin (Act88F), myosin (Mhc and Mlc2) and tropomyosin (Tm2) as well as a group of closely linked, muscle-related genes regulated by a haplolethal sequence within an intron of the Troponin I (wupA) gene [24]. This subset may also include Hdl, which may correspond to Troponin T (up) [21, 25]. Like the ribosomal protein genes, these genes may be particularly dosage sensitive because muscle assembly requires minimal levels or a particular stoichiometry of component proteins [22, 24, 26]. Other definable subsets encode homeodomain proteins (Abd-B, Dll, Scr and Ubx), Notch pathway components (Dl, H and N), Polycomb group repressor proteins (Pc and Pcl), apoptosis regulators (lok and p53) and melanin biosynthetic enzymes (b and e).

Assembling a new Deficiency Kit

For years, the BDSC has maintained and distributed a set of deletions providing maximal genomic coverage using the fewest deletions deemed practical. The 'Deficiency Kit' was established before FLP-FRT deletions existed and deletions were added to the original deletion set when they provided novel coverage. Consequently, the kit consisted mostly of deletions with breakpoints mapped only at the level of polytene cytology. With large collections of molecularly defined deletions available, the kit had become outdated.

As the second goal of our project, we created a new Deficiency Kit and began distributing it in July 2009 while we were still screening for new deletions. We chose a tiling path of FRT-derived deletions from the three deletion projects and from individual investigators to provide maximal coverage. The remaining gaps in coverage were filled with older deletions where possible. In many cases, the older deletions were maintained in stock with duplications, allowing us to delete many regions containing haploinsufficient loci. To create a kit better suited for genetic modifier screens, we included the deletion pairs most closely flanking haploinsufficient genes even when they were removed completely by older deletions. Older deletions also allowed us to cover heterochromatic regions where FLP-FRT deletions currently do not exist and the fourth chromosome where few molecularly defined deletions exist. As we generated new deletions, we updated and improved the kit. Table 5 shows a summary of genomic coverage provided by the current kit. Additional file 4 lists the contents of the Deficiency Kit and Additional file 2 shows the deletion breakpoints relative to annotated genes. Information on the Deficiency Kit is also provided on the BDSC website [27].

Table 5 The Bloomington Deficiency Kit

Because the average size of the DrosDel, Exelixis and BSC deletions is smaller than the average size of the deletions in the previous Deficiency Kit, the new kit is larger (468 versus 271 deletions). Nevertheless, the advantages of molecularly mapped breakpoints and consistent genetic backgrounds offset the inconvenience of handling more stocks in experiments. As we will show in the next section, once a phenotype is mapped to a chromosomal interval using the new Deficiency Kit, its position can be narrowed to a very small interval or even a single transcription unit using additional molecularly defined deletions.

Extensive subdivision of the genome with deletion breakpoints

When a gene is mapped with deletions, it is usually localized to a chromosomal interval defined by adjacent breakpoints of different deletions rather than the breakpoints of a single deletion. For example, a mutation would map between the distal breakpoints of two overlapping deletions if it complemented the proximal deletion but not the distal deletion. Consequently, the critical factor determining the resolution of deletion mapping is the density of deletion breakpoints rather than deletion size.

The third goal of our project was to improve the distribution of deletion breakpoints across the genome and our target was to generate deletions so that the intervals between adjacent breakpoints would contain no more than a dozen protein-coding genes. We examined the distribution of breakpoints provided by the DrosDel and Exelixis deletions and undertook screens to interpose new breakpoints between them (Additional file 2). Figure 2 shows the distribution of interval sizes provided by deletions from the three projects. To construct this distribution, we assigned every euchromatic gene a position between adjacent deletion breakpoints based on where it would map uniquely if a mutation in the gene were complementation tested with the deletions. These hypothetical complementation tests allowed us to assign a gene to a single interval even when a deletion breakpoint fell within the gene. The median interval size for the entire genome is nine genes. The median interval size for chromosome arms X, 2R, 3L and 3R is also nine genes, but the median size for 2L, where there are more FRT-derived deletions, is only seven genes.

Figure 2
figure 2

Frequency distribution of the number of genes between molecularly defined deletion breakpoints. The number of genes in a chromosomal interval between adjacent deletion breakpoints is shown on the x-axis and the number of intervals in the Drosophila genome with those sizes is shown on the y-axis. Because BSC, Exelixis and DrosDel deletions overlap extensively, the intervals were usually defined by breakpoints of different deletions. The median interval size is only nine genes. For simplicity, the Stellate gene cluster in chromosomal region 12E, the histone gene cluster in 39D and the 5S rRNA gene cluster in 56E were counted as single genes.

Many intervals exceeded our target size of 12 genes, though the largest interval contains only 36 genes (it lies on the fourth chromosome, where our screens were unsuccessful; the largest interval elsewhere contains 27 genes). The largest intervals could not be subdivided because no FRT-bearing transposon insertions existed in the interval or repeated attempts to subdivide the region with existing insertions failed. Other intervals exceeded our target size because gene annotations changed in the course of the project or because we counted genes for noncoding RNAs such as microRNAs and small nucleolar RNAs, which we did not consider in planning the deletion screens due to their poor annotation at the time.

As Figure 2 shows, it is now possible to map most mutations and genetic modifiers to very small groups of candidate transcription units with simple crosses. The deletion breakpoints even define 377 single-gene intervals. This degree of genomic subdivision makes deletion mapping the method of choice for most gene localization experiments in D. melanogaster.

Discussion

The enormous benefit of using FLP-mediated recombination between FRT-bearing transposon insertions to generate chromosomal deletions is absolute predictability. Deletion breakpoints can be positioned with single-nucleotide certainty. Homologous recombination methods [28] are the only other approaches in D. melanogaster allowing one to plan deletion screens with such precision. The advent of these molecularly defined deletions fundamentally changed the experience of deletion mapping. Previously, the breakpoints of most deletions had been mapped only at the resolution of polytene chromosome analysis, which, even at its best, positioned breakpoints in relatively large chromosomal intervals. The uncertainty of breakpoint positions in older deletions made high-resolution mapping difficult. In contrast, deletion mapping with FLP-FRT deletions is straightforward and unambiguous. The obvious advantages of these deletions led to their immediate and enthusiastic adoption by Drosophila geneticists. While the older deletions are still valuable, they are not used as frequently as they once were. Approximately 85% of deletion stocks currently ordered from the BDSC are FRT-derived deletions.

Our project at the BDSC was one of three large-scale projects generating FLP-FRT deletions. The Exelixis and DrosDel Projects also generated extensive genomic coverage, but our project is distinguished by its emphasis on improving deletion coverage by minimizing undeleted regions containing haploinsufficient genes and by its emphasis on genomic subdivision. Haplolethal and haplosterile loci were not a big consideration in the planning of the Exelixis deletions while the relatively low density of the DrosDel FRT-bearing transposons across the genome did not, on average, allow the DrosDel deletions to confine haploinsufficient genes to small undeleted intervals. Our project stressed flanking haploinsufficient genes as closely as possible with pairs of deletions as the most practical approach to improving deletion coverage even before we started using the FLP-FRT system. The deletions we generated by the hybrid element insertion method [5] were also designed to provide novel coverage. Our attention to haploinsufficient genes is largely responsible for improving overall genome coverage beyond that provided by either the older deletions or the DrosDel and Exelixis collections. It has also helped produce the most comprehensive inventory of haploinsufficient genes to date in Drosophila or any other multicellular eukaryote.

The first chromosomal deletion was discovered in D. melanogaster [29] and Drosophila now boasts the best genomic deletion coverage and subdivision of any animal. Total coverage is even roughly equivalent to that in the yeasts Saccharomyces cerevisiae and Schizosaccaromyces pombe, where high-throughput homologous recombination approaches are available [30, 31]. Despite this accomplishment, our approach of flanking haploinsufficient genes as closely as possible with deletions, rather than first producing new duplications and subsequently using them in deletion screens, precluded complete genomic coverage. The recent generation of X and second chromosome duplications [9, 21] and the development of ΦC31 transformation technology for the creation of new multigene duplications [25, 32] will make it much easier to isolate deletions for chromosomal regions containing haploinsufficient genes. Further improvements in deletion coverage will doubtlessly occur as these resources and techniques are exploited.

It is remarkable how few haploinsufficient genes are present in the D. melanogaster genome - particularly when one considers that the vast majority fall into a single class, the Minute genes, and produce a consistent set of phenotypes. There are very few haplolethal or haplosterile genes of any kind. In fact, most deletions in Drosophila appear to have no dominant phenotypic effect and, on average, one-twentieth of a chromosome arm must be deleted before significant aneuploidy-associated phenotypes are seen [33, 34]. This stands in stark contrast to humans, where haploinsufficiency causes hundreds of genetic diseases [35–38]. Even though the Drosophila genome has never been systematically surveyed for weak fitness effects arising from reduced gene dosage and one could argue that subtle mutant phenotypes are more easily recognized in humans, the tolerance Drosophila shows to hypoploidy must reflect fundamental differences in the abilities of organisms to compensate for changes in gene dosage.

The extensive genomic coverage and subdivision provided by molecularly defined deletions in Drosophila presents researchers with the ability to map mutations and genetic modifiers to very small chromosomal intervals with a few rounds of simple crosses. While it is sometimes necessary to generate additional deletions for further localization, it is often possible to identify an obvious candidate based on its sequence characteristics or previously characterized mutant phenotypes. Gene identification is much easier and more efficient than it was even a decade ago.

The mapping of X-linked mutations and modifiers benefits from the X chromosome deletions described here as well as duplications of X chromosomes that were generated at the BDSC in a concurrent project [21]. We have now isolated Y-linked duplications for 97% of the X chromosome (data not shown). These duplications can be paired with X chromosome deletions to allow inheritance of the deletions from fathers (indeed many stocks with deletion-duplication pairs are already distributed by the BDSC) and they can be used to map mutations by phenotypic rescue and identify hyperploid modifiers. Using the molecularly defined X chromosome deletions, the Y-linked X chromosome duplications and interchromosomal duplications of X segments generated by ΦC31 transformation [25], phenotypes can be mapped with near single-gene resolution.

Conclusions

We have presented the results of a large-scale resource development project at the BDSC. We substantially improved the overall genomic coverage provided by chromosomal deletions as well as the density of deletion breakpoints. Our efforts will simplify experiments utilizing deletions to map mutations or identify dosage-sensitive enhancers and suppressors of mutant phenotypes. To aid researchers in conducting such experiments efficiently, we organized deletions providing maximal genomic coverage into a Deficiency Kit. We also made a complete inventory of haplolethal and haplosterile genes and compiled information about other haploinsufficient genes in D. melanogaster in the process of improving deletion coverage. This information should contribute to future studies of gene dosage relationships and aneuploidy effects.

Chromosomal deletions are one of the most versatile and important genetic tools available to Drosophila researchers and the stocks generated in this project will contribute to the success of experimental studies for years to come. Resource development projects like ours and their targeted funding propel research by creating new experimental opportunities and they are vital to the continued vigor and advancement of scientific inquiry in genetic model organisms.

Materials and methods

Sources of stocks

Most P{XP}, PBac{WH} and PBac{RB} insertions used in deletion screens were obtained directly from Exelixis, Inc. The remainder came from the BDSC and the Exelixis collection at Harvard University. Other stocks came from Exelixis, the BDSC, the Harvard collection, the Szeged Stock Centre, the Drosophila Genetic Resource Center at the Kyoto Institute of Technology and the laboratory of Michael Ashburner at the University of Cambridge.

Chromosomal positions

All genomic coordinates and gene counts are based on Genome Release 5.16. Cytological band positions were predicted from Release 5 coordinates using FlyBase map conversion tables [39, 40]. We artificially defined the euchromatin/heterochromatin boundary as the most proximal extent of the contiguous genomic assembly for each chromosome arm in Genome Release 5.

The genomic coordinates of P{XP}, PBac{WH} and PBac{RB} insertions were initially determined at Exelixis and most were published in Thibault et al.[10]. Coordinates of unpublished insertions were obtained directly from Exelixis. Flanking sequences were reanalyzed by Roger Hoskins and Joe Carlson of Lawrence Berkeley National Laboratory and they submitted revised coordinates for many insertions to FlyBase. We used these reanalyzed coordinates where available. Otherwise, we used coordinates from a preliminary reanalysis (R Hoskins, personal communication) or the original Exelixis coordinates in that order of preference. The origins of coordinates are recorded in Additional file 1.

Deletion screens

Parks et al.[5] described methods for using P{XP}, PBac{WH} and PBac{RB} insertions to isolate chromosomal deletions in a standard genetic background. We followed the basic schemes developed at Exelixis for deletion screening, but we introduced some alternative stocks to improve the vigor of flies in the crosses and final deletion stocks. We built these stocks by carefully introgressing balancer and marker chromosomes into the standard background using chromosome substitution schemes patterned on Craymer [41]. The source of the Y and fourth chromosomes was not controlled. (Details of stock constructions will be provided upon request.) Here we present the crosses we eventually found to work best for the isolation of X, second and third chromosome deletions. We identified deletions based on loss of miniwhite or the presence of two copies of miniwhite; see Park et al.[5] and Results for the rationale. Tn{w+}1 and Tn{w+}2 are used generically to denote different miniwhite-marked transposon (Tn) insertions of P{XP}, PBac{WH} and PBac{RB}. Deletions were named with a BSC (Bloomington Stock Center) prefix followed by a number in the order they were placed into public distribution.

X chromosome deletions

G0: w1118 Tn{w+}1 ♀♀ × w1118/Y; MKRS, P{hsFLP}86E/TM6B, Tb1 ♂♂ (9 vials of 8 ♀♀ × 5 ♂♂, subcultured once).

G1: w1118 Tn{w+}2 ♀♀ × w1118 Tn{w+}1/Y; +/MKRS, P{hsFLP}86E ♂♂ (9 vials of 8 ♀♀ × 5 ♂♂, subcultured once).

G2: w1118 Tn{w+}2/w1118 Tn{w+}1; +/MKRS, P{hsFLP}86E ♀♀ × FM7h/Dp(2;Y)G, P{hs-hid}Y ♂♂ (9 vials of 8 ♀♀ × 5 ♂♂, subcultured once; these females were heat shocked as larvae at 37°C for 1 hour on days 4 to 8 after egg lay).

G3: Df(1)BSC, w1118/FM7h ♀ × FM7h/Dp(2;Y)G, P{hs-hid}Y ♂♂ (deletion-bearing females were recognized as white-eyed or darker-eyed flies, depending on the screen; we typically established stocks of five independent putative deletions for testing).

Although the temperature-sensitive lethal P{hs-hid} construct can be used to eliminate unwanted progeny classes, it was incidental to our choice of stocks in the preceding and succeeding crosses and not used. We include it here only for completeness.

Second chromosome deletions {3rd level heading}

G0: P{hsFLP}1, y1 w1118; P{hs-hid}2, wgSp-1/CyO ♀♀ × w1118; Tn{w+}1 ♂♂ (5 vials of 8 ♀♀ × 5 ♂♂, subcultured once).

G1: P{hsFLP}1, y1 w1118/w1118; CyO/Tn{w+}1 ♀♀ × w1118 /Y; Tn{w+}2 ♂♂ (5 vials of 8 ♀♀ × 5 ♂♂, subcultured once).

G2: w1118; wgSp-1/SM6a ♀♀ × P{hsFLP}1, y1 w1118/Y; Tn{w+}1/Tn{w+}2 ♂♂ (5 vials of 8 ♀♀ × 5 ♂♂, subcultured once; these males were heat shocked as larvae at 37°C for 1 hour on days 4 to 8 after egg lay).

G3: w1118; wgSp-1/SM6a ♀♀ × w1118/Y; Df(2)BSC/SM6a ♂ (deletion-bearing males were recognized as white-eyed or darker-eyed flies, depending on the screen; we typically established stocks of five independent putative deletions for testing.)

G4: w1118; Df(2)BSC/SM6a ♀♀ × w1118; Df(2)BSC/SM6a ♂♂.

Third chromosome deletions {3rd level heading}

G0: P{hsFLP}1, y1 w1118; Dr1/TM3, Sb1 ♀♀ × w1118; Tn{w+}1 ♂♂ (5 vials of 8 ♀♀ × 5 ♂♂, subcultured once).

G1: P{hsFLP}1, y1 w1118/w1118; TM3, Sb1/Tn{w+}1 ♀♀ × w1118/Y; Tn{w+}2 ♂♂ (5 vials of 8 ♀♀ × 5 ♂♂, subcultured once).

G2: w1118; P{hs-hid}3, Dr1/TM6C, cu1 Sb1 ♀♀ × P{hsFLP}1, y1 w1118/Y; Tn{w+}1/Tn{w+}2 ♂♂

(5 vials of 8 ♀♀ × 5 ♂♂, subcultured once; these males were heat shocked as larvae at 37°C for 1 hour on days 4 to 8 after egg lay).

G3: w1118; P{hs-hid}3, Dr1/TM6C, cu1 Sb1 ♀♀ × w1118/Y; Df(3)BSC/TM6C, cu1 Sb1 ♂ (deletion-bearing males were recognized as white-eyed or darker-eyed flies, depending on the screen; we typically established stocks of five independent putative deletions for testing).

G4: w1118; Df(3)BSC/TM6C, cu1 Sb1 ♀♀ × w1118; Df(3)BSC/TM6C, cu1 Sb1 ♂♂.

It was necessary to establish stocks of a few deletions in nonstandardized backgrounds to overcome problems arising from noncomplementation or haploinsufficiency, but these stocks are noted in Additional file 1 and BDSC stock records.

Verifying new deletions

Four approaches were used to verify that the planned deletion was recovered from a screen: complementation tests, polytene chromosome cytology, 'hybrid PCR' and PCR with primers flanking transposable element insertions. The complementation tests were routine and require no elaboration. Polytene chromosomes were analyzed in standard lacto-aceto-orcein preparations [42] using the maps of Lefevre [43] and Sorsa et al.[44, 45].

The 'hybrid PCR' approach was described in Parks et al.[5]. Briefly, primers were used to amplify a DNA fragment from the recombinant transposable element formed upon FLP-mediated recombination between two progenitor FRT-bearing insertions. We used the primers described in Parks et al.[5] with the substitution of a new primer (5'-GCTTCTAAACGCTTACGCATAAAC GATG-3') for the RB3' plus and RB3' minus primer. We were unable to establish Hybrid PCR primers and conditions to verify deletions identified from an increase in miniwhite copy number; consequently, we tested only deletions detected by loss of miniwhite markers.

To verify that a particular chromosomal region was removed by a deletion, we designed PCR primers flanking the insertion site of a transposon located within the deleted region. From crosses, we isolated progeny carrying both the transposon and the putative deficiency and isolated DNA. With short extension times, a PCR fragment is amplified only when there is no transposon between the primer sites; consequently, a fragment is recovered if the noninsertion chromosome lacks a deletion. Since this approach verifies a deletion by a negative result, we repeated the test with three to five independent DNA preparations and appropriate positive controls. Transposon insertions and primer sequences used in these tests are given in Additional file 1. Additional details are available from Cook et al.[21] where similar tests were used to confirm chromosomal duplications.

Stocks were placed into public distribution by the BDSC as soon as we verified deletions. We submitted a description of each deletion to FlyBase and Additional file 1 provides citations to these reports.

Stock availability

Stocks described in this paper may be obtained from the BDSC. Information on ordering, web pages devoted to deletion collections and the Deficiency Kit and lists of other stocks are available on the BDSC website [46]. FlyBase [39] provides information on individual genetic components with cross references to BDSC stocks.

Abbreviations

BDSC:

Bloomington Drosophila Stock Center

FLP:

FLP recombinase

FRT :

FLP recognition target

PCR:

polymerase chain reaction

SNP:

single-nucleotide polymorphism.

References

  1. Gray YH, Tanaka MM, Sved JA: P-element-induced recombination in Drosophila melanogaster: hybrid element insertion. Genetics. 1996, 144: 1601-1610.

    PubMed  CAS  PubMed Central  Google Scholar 

  2. Preston CR, Sved JA, Engels WR: Flanking duplications and deletions associated with P-induced male recombination in Drosophila. Genetics. 1996, 144: 1623-1638.

    PubMed  CAS  PubMed Central  Google Scholar 

  3. Huet F, Lu JT, Myrick KV, Baugh LR, Crosby MA, Gelbart WM: A deletion-generator compound element allows deletion saturation analysis for genomewide phenotypic annotation. Proc Natl Acad Sci USA. 2002, 99: 9948-9953. 10.1073/pnas.142310099.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Mohr SE, Gelbart WM: Using the P{wHy} hybrid transposable element to disrupt genes in region 54D-55B in Drosophila melanogaster. Genetics. 2002, 162: 165-176.

    PubMed  CAS  PubMed Central  Google Scholar 

  5. Parks AL, Cook KR, Belvin M, Dompe NA, Fawcett R, Huppert K, Tan LR, Winter CG, Bogart KP, Deal JE, Deal-Herr ME, Grant D, Marcinko M, Miyazaki WY, Robertson S, Shaw KJ, Tabios M, Vysotskaia V, Zhao L, Andrade RS, Edgar KA, Howie E, Killpack K, Milash B, Norton A, Thao D, Whittaker K, Winner MA, Friedman L, Margolis J, et al: Systematic generation of high-resolution deletion coverage of the Drosophila melanogaster genome. Nat Genet. 2004, 36: 288-292. 10.1038/ng1312.

    Article  PubMed  CAS  Google Scholar 

  6. Metaxakis A, Oehler S, Klinakis A, Savakis C: Minos as a genetic and genomic tool in Drosophila melanogaster. Genetics. 2005, 171: 571-581. 10.1534/genetics.105.041848.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Golic KG, Golic MM: Engineering the Drosophila genome: chromosome rearrangements by design. Genetics. 1996, 144: 1693-1711.

    PubMed  CAS  PubMed Central  Google Scholar 

  8. Ryder E, Blows F, Ashburner M, Bautista-Llacer R, Coulson D, Drummond J, Webster J, Gubb D, Gunton N, Johnson G, O'Kane CJ, Huen D, Sharma P, Asztalos Z, Baisch H, Schulze J, Kube M, Kittlaus K, Reuter G, Maroy P, Szidonya J, Rasmuson-Lestander A, Ekstrom K, Dickson B, Hugentobler C, Stocker H, Hafen E, Lepesant JA, Pflugfelder G, Heisenberg M, et al: The DrosDel collection: a set of P-element insertions for generating custom chromosomal aberrations in Drosophila melanogaster. Genetics. 2004, 167: 797-813. 10.1534/genetics.104.026658.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Ryder E, Ashburner M, Bautista-Llacer R, Drummond J, Webster J, Johnson G, Morley T, Chan YS, Blows F, Coulson D, Reuter G, Baisch H, Apelt C, Kauk A, Rudolph T, Kube M, Klimm M, Nickel C, Szidonya J, Maroy P, Pal M, Rasmuson-Lestander A, Ekstrom K, Stocker H, Hugentobler C, Hafen E, Gubb D, Pflugfelder G, Dorner C, Mechler B, et al: The DrosDel deletion collection: a Drosophila genomewide chromosomal deficiency resource. Genetics. 2007, 177: 615-629. 10.1534/genetics.107.076216.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Thibault ST, Singer MA, Miyazaki WY, Milash B, Dompe NA, Singh CM, Buchholz R, Demsky M, Fawcett R, Francis-Lang HL, Ryner L, Cheung LM, Chong A, Erickson C, Fisher WW, Greer K, Hartouni SR, Howie E, Jakkula L, Joo D, Killpack K, Laufer A, Mazzotta J, Smith RD, Stevens LM, Stuber C, Tan LR, Ventura R, Woo A, Zakrajsek I, et al: A complementary transposon tool kit for Drosophila melanogaster using P and piggyBac. Nat Genet. 2004, 36: 283-287. 10.1038/ng1314.

    Article  PubMed  CAS  Google Scholar 

  11. Anonymous: Release the flies. Nat Genet. 2004, 36: 205-205.

    Article  Google Scholar 

  12. Jang C, Lee G, Chung J: LKB1 induces apical trafficking of Silnoon, a monocarboxylate transporter, in Drosophila melanogaster. J Cell Biol. 2008, 183: 11-17. 10.1083/jcb.200807052.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Natzle JE, Kiger JA, Green MM: Bursicon signaling mutations separate the epithelial-mesenchymal transition from programmed cell death during Drosophila melanogaster wing maturation. Genetics. 2008, 180: 885-893. 10.1534/genetics.108.092908.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Xie X, Fischer JA: On the roles of the Drosophila KASH domain proteins Msp-300 and Klarsicht. Fly. 2008, 2: 74-81.

    Article  PubMed  Google Scholar 

  15. Dickman DK, Davis GW: The schizophrenia susceptibility gene dysbindin controls synaptic homeostasis. Science. 2009, 326: 1127-1130. 10.1126/science.1179685.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Cooper MT, Conant AW, Kennison JA: Molecular genetic analysis of Chd3 and polytene chromosome region 76B-D in Drosophila melanogaster. Genetics. 2010, 185: 811-822. 10.1534/genetics.110.115121.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Gehrke S, Imai Y, Sokol N, Lu B: Pathogenic LRRK2 negatively regulates microRNA-mediated translational repression. Nature. 2010, 466: 637-641. 10.1038/nature09191.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Wright AP, Fox AN, Johnson KG, Zinn K: Systematic screening of Drosophila deficiency mutatations for embryonic phenotypes and orphan receptor ligands. PLoS One. 2010, 5: e12288-10.1371/journal.pone.0012288.

    Article  PubMed  PubMed Central  Google Scholar 

  19. FlyBase: D.melanogaster Gene Disruptions, Stocks and Reagents. [http://flybase.org/cgi-bin/gbrowse/dmel/]

  20. Bloomington Drosophila Stock Center: Deletion List. [http://flystocks.bio.indiana.edu/Browse/df/dftop.htm]

  21. Cook RK, Deal ME, Deal JA, Garton RD, Brown CA, Ward ME, Andrade RS, Spana EP, Kaufman TC, Cook KR: A new resource for characterizing X-linked genes in Drosophila melanogaster: systematic coverage and subdivision of the X chromosome with nested, Y-linked duplications. Genetics. 2010, 186: 1095-1109. 10.1534/genetics.110.123265.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Marygold SJ, Roote J, Reuter G, Lambertsson A, Ashburner M, Millburn GH, Harrison PM, Yu Z, Kenmochi N, Kaufman TC, Leevers SJ, Cook KR: The ribosomal protein genes and Minute loci of Drosophila melanogaster. Genome Biol. 2007, 8: R216-10.1186/gb-2007-8-10-r216.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Dorer DR, Rudnick JA, Moriyama EN, Christensen AC: A family of genes clustered at the Triplo-lethal locus of Drosophila melanogaster has an unusual evolutionary history and significant synteny with Anopheles gambiae. Genetics. 2003, 165: 613-621.

    PubMed  CAS  PubMed Central  Google Scholar 

  24. Prado A, Canal I, Ferrus A: The haplolethal region at the 16F gene cluster of Drosophila melanogaster: structure and function. Genetics. 1999, 151: 163-175.

    PubMed  CAS  PubMed Central  Google Scholar 

  25. Venken KJT, Popodi E, Holtzman SL, Schulze KL, Park S, Carlson JW, Hoskins RA, Bellen HJ, Kaufman TC: A molecularly defined duplication set for the X chromosome of Drosophila melanogaster. Genetics. 2010, 186: 1111-1125. 10.1534/genetics.110.121285.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Homyk T, Emerson CP: Functional interactions between unlinked muscle genes within haploinsufficient regions of the Drosophila genome. Genetics. 1988, 119: 105-121.

    PubMed  Google Scholar 

  27. Bloomington Drosophila Stock Center: Deficiency Kit information. [http://flystocks.bio.indiana.edu/Browse/df/dfkit-info.htm]

  28. Huang J, Zhou W, Dong W, Hong Y: Targeted engineering of the Drosophila genome. Fly. 2009, 3: 274-277.

    Article  PubMed  CAS  Google Scholar 

  29. Bridges CB: Deficiency. Genetics. 1917, 2: 445-465.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. Kim D-U, Hayles J, Kim D, Wood V, Park H-O, Won M, Yoo H-S, Duhig T, Nam M, Palmer G, Han S, Jeffery L, Baek S-T, Lee H, Shim YS, Lee M, Kim L, Heo K-S, Noh EJ, Lee A-R, Jang Y-J, Chung K-S, Choi S-J, Park J-Y, Park Y, Kim HM, Park S-K, Park H-J, Kang E-J, Kim HB, et al: Analysis of a genome-wide set of gene deletions in the fission yeast Schizosaccharomyces pombe. Nat Biotechnol. 2010, 28: 1308-

    Article  CAS  Google Scholar 

  31. Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, Arkin AP, Astromoff A, El-Bakkoury M, Bangham R, Benito R, Brachat S, Campanaro S, Curtiss M, Davis K, Deutschbauer A, Entian K-D, Flaherty P, Foury F, Garfinkel DJ, Gerstein M, Gotte D, Guldener U, Hegemann JH, Hempel S, Herman Z, et al: Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002, 418: 387-391. 10.1038/nature00935.

    Article  PubMed  CAS  Google Scholar 

  32. Venken KJ, He Y, Hoskins RA, Bellen HJ: P[acman]: a BAC transgenic platform for targeted insertion of large DNA fragments in D. melanogaster. Science. 2006, 314: 1747-1751. 10.1126/science.1134426.

    Article  PubMed  CAS  Google Scholar 

  33. Lindsley DL, Sandler L, Baker BS, Carpenter AT, Denell RE, Hall JC, Jacobs PA, Miklos GL, Davis BK, Gethmann RC, Hardy RW, Steven AH, Miller M, Nozawa H, Parry DM, Gould-Somero M: Segmental aneuploidy and the genetic gross structure of the Drosophila genome. Genetics. 1972, 71: 157-184.

    PubMed  CAS  PubMed Central  Google Scholar 

  34. Ashburner M, Golic KG, Hawley RS: Drosophila: a Laboratory Handbook. 2005, Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press, 2

    Google Scholar 

  35. Seidman JG, Seidman C: Transcription factor haploinsufficiency: when half a loaf is not enough. J Clin Invest. 2002, 109: 451-455.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Dang VT, Kassahn KS, Marcos AE, Ragan MA: Identification of human haploinsufficient genes and their genomic proximity to segmental duplications. Eur J Hum Genet. 2008, 16: 1350-1357. 10.1038/ejhg.2008.111.

    Article  PubMed  CAS  Google Scholar 

  37. Huang N, Lee I, Marcotte EM, Hurles ME: Characterising and predicting haploinsufficiency in the human genome. PLoS Genet. 2010, 6: e1001154-10.1371/journal.pgen.1001154.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Kondrashov FA, Koonin EV: A common framework for understanding the origin of genetic dominance and evolutionary fates of gene duplications. Trends Genet. 2004, 20: 287-290. 10.1016/j.tig.2004.05.001.

    Article  PubMed  CAS  Google Scholar 

  39. FlyBase: A Database of Drosophila Genes & Genomes. [http://flybase.org/]

  40. Tweedie S, Ashburner M, Falls K, Leyland P, McQuilton P, Marygold S, Millburn G, Osumi-Sutherland D, Schroeder A, Seal R, Zhang H: FlyBase: enhancing Drosophila Gene Ontology annotations. Nucleic Acids Res. 2009, 37: D555-559. 10.1093/nar/gkn788.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Craymer L: A procedure for constructing isogenic stocks. Drosophila Information Service. 1984, 60: 78-79.

    Google Scholar 

  42. Carpenter AT: Salivary chromosome analysis of aberrations. Methods Mol Biol. 2004, 247: 257-277.

    PubMed  Google Scholar 

  43. Lefevre G: A photographic representation and interpretation of the polytene chromosomes of Drosophila melanogaster salivary glands. Genetics and Biology of Drosophila. Edited by: Ashburner M, Novitski E. 1976, New York: Academic Press, 1a: 31-66.

    Google Scholar 

  44. Sorsa V: Chromosome Maps of Drosophila. 1988, Boca Raton, FL: CRC Press

    Google Scholar 

  45. Electron micrograph maps of Drosophila melanogaster polytene chromosomes. [http://www.helsinki.fi/~saura/EM/index.html]

  46. Bloomington Drosophila Stock Center. [http://flystocks.bio.indiana.edu/]

  47. Mohr OL: Character changes caused by mutation of an entire region of a chromosome in Drosophila. Genetics. 1919, 4: 275-282.

    PubMed  CAS  PubMed Central  Google Scholar 

  48. Lindsley DL, Grell EH: Genetic Variations of Drosophila melanogaster. 1968, Washington, DC: Carnegie Institution of Washington

    Google Scholar 

  49. Gergen JP, Wieschaus E: Dosage requirements for runt in the segmentation of Drosophila embryos. Cell. 1986, 45: 289-299. 10.1016/0092-8674(86)90393-4.

    Article  PubMed  CAS  Google Scholar 

  50. Lewis EB: The relation of repeats to position effect in Drosophila melanogaster. Genetics. 1945, 30: 137-166.

    PubMed  CAS  PubMed Central  Google Scholar 

  51. Gao Z, Joseph E, Ruden DM, Lu X: Drosophila Pkd2 is haploid-insufficient for mediating optimal smooth muscle contractility. J Biol Chem. 2004, 279: 14225-14231. 10.1074/jbc.M312223200.

    Article  PubMed  CAS  Google Scholar 

  52. Bridges CB, Brehme KS: The Mutants of Drosophila melanogaster. 1944, Washington, DC: Carnegie Institution of Washington

    Google Scholar 

  53. Bernstein SI, Mogami K, Donady JJ, Emerson CP: Drosophila muscle myosin heavy chain encoded by a single gene in a cluster of muscle mutations. Nature. 1983, 302: 393-397. 10.1038/302393a0.

    Article  PubMed  CAS  Google Scholar 

  54. Kurzhals RL, Titen SWA, Xie HB, Golic KG: Chk2 and p53 are haploinsufficient with dependent and independent functions to eliminate cells after telomere loss. PLoS Genet. 2011, 7: e1002103-10.1371/journal.pgen.1002103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Green MM: A study in gene action using different dosages and alleles of vestigial in Drosophila melanogaster. Genetics. 1946, 31: 1-20.

    PubMed Central  Google Scholar 

  56. Bridges CB, Skoog EN, Li JC: Genetical and cytological studies of a deficiency (Notopleural) in the second chromosome of Drosophila melanogaster. Genetics. 1936, 21: 788-795.

    PubMed  CAS  PubMed Central  Google Scholar 

  57. Duncan IM: Polycomblike: a gene that appears to be required for the normal expression of the Bithorax and Antennapedia gene complexes of Drosophila melanogaster. Genetics. 1982, 102: 49-70.

    PubMed  CAS  PubMed Central  Google Scholar 

  58. Li JC: The effect of chromosome aberrations on development in Drosophila melanogaster. Genetics. 1927, 12: 1-58.

    PubMed  CAS  PubMed Central  Google Scholar 

  59. Sunkel CE, Whittle JRS: Brista: a gene involved in the specification and differentiation of distal cephalic and thoracic structures in Drosophila melanogaster. Roux's Arch Dev Biol. 1987, 196: 124-132. 10.1007/BF00402034.

    Article  Google Scholar 

  60. Wieschaus E, Nusslein-Volhard C, Kluding H: Krüppel, a gene whose activity is required early in the zygotic genome for normal embryonic segmentation. Dev Biol. 1984, 104: 172-186. 10.1016/0012-1606(84)90046-0.

    Article  PubMed  CAS  Google Scholar 

  61. Harris D, Orme C, Kramer J, Namba L, Champion M, Palladino MJ, Natzle J, Hawley RS: A deficiency screen of the major autosomes identifies a gene (matrimony) that is haplo-insufficient for achiasmate segregation in Drosophila oocytes. Genetics. 2003, 165: 637-652.

    PubMed  CAS  PubMed Central  Google Scholar 

  62. Denell RE: Homoeosis in Drosophila. II. a genetic analysis of Polycomb. Genetics. 1978, 90: 277-289.

    PubMed  CAS  PubMed Central  Google Scholar 

  63. Kaufman TC, Lewis R, Wakimoto B: Cytogenetic analysis of chromosome 3 in Drosophila melanogaster: the homoeotic gene complex in polytene chromosome interval 84A-B. Genetics. 1980, 94: 115-133.

    PubMed  CAS  PubMed Central  Google Scholar 

  64. Karlik CC, Fyrberg EA: An insertion within a variably spliced Drosophila tropomyosin gene blocks accumulation of only one encoded isoform. Cell. 1985, 41: 57-66. 10.1016/0092-8674(85)90061-3.

    Article  PubMed  CAS  Google Scholar 

  65. Karlik CC, Coutu MD, Fyrberg EA: A nonsense mutation within the act88F actin gene disrupts myofibril formation in Drosophila indirect flight muscles. Cell. 1984, 38: 711-719. 10.1016/0092-8674(84)90266-6.

    Article  PubMed  CAS  Google Scholar 

  66. Lewis EB: Genes and developmental pathways. Am Zool. 1963, 3: 33-56.

    Article  Google Scholar 

  67. Bridges CB, Morgan TH: The third-chromosome group of mutant characters of Drosophila melanogaster. Carnegie Inst Publ. 1923, 327: 1-251.

    Google Scholar 

  68. Jarecki J, Johnson E, Krasnow MA: Oxygen regulation of airway branching in Drosophila is mediated by branchless FGF. Cell. 1999, 99: 211-220. 10.1016/S0092-8674(00)81652-9.

    Article  PubMed  CAS  Google Scholar 

  69. Warmke JW, Kreuz AJ, Falkenthal S: Co-localization to chromosome bands 99E1-3 of the Drosophila melanogaster myosin light chain-2 gene and a haplo-insufficient locus that affects flight behavior. Genetics. 1989, 122: 139-151.

    PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgements

This project would have been impossible without the work of Exelixis scientists who designed the FLP-FRT deletion screening system and the generosity of Exelixis leadership in releasing strains into the public domain. We thank John Roote, Ed Ryder, Steve Russell, Michael Ashburner and other DrosDel Project members for their cooperation, Steven Marygold and Andrew Lambertsson for conversations about haploinsufficient genes, Roger Hoskins, Joe Carlson and Bob Levis for characterizing transgene insertions, Eric Spana for contributing microarray data and our FlyBase colleagues for designing the Gbrowse aberrations display and curating deletion information. We are grateful for the help and support of our Bloomington Stock Center coworkers Annette Parks, Kevin Bogart, Millie Winner, Nela Popadic, Megan Ward, Adam Brown, Russell Garton, Luke Jacobus and especially Kathy Matthews. This work was supported primarily by National Center for Research Resources grant RR014106 to KRC with additional funding from NSF Division of Biological Infrastructure Grant DBI-9816125 to KRC and TCK, Indiana Genomics Initiative funding to TCK and Bloomington Stock Center funds.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kevin R Cook.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RKC, SJC, JAD, RAC, MED and JMG planned and conducted genetic screens and characterized new deletions. KRC initiated and supervised the project, characterized deletions, analyzed data, prepared stocks for public distribution and wrote the manuscript. RKC analyzed data and revised the manuscript. TCK provided support and guidance. All authors read and approved the final manuscript.

Electronic supplementary material

13059_2012_2771_MOESM1_ESM.XLSX

Additional file 1: Breakpoints and characterization of BSC deletions. This file provides full information for deletions isolated in this project, including progenitor FRT insertions and their locations, selection criteria in screens, verification tests and references to FlyBase documentation. (XLSX 170 KB)

13059_2012_2771_MOESM2_ESM.XLSX

Additional file 2: Map of Xand fourth chromosome deletion breakpoints. This spreadsheet shows the breakpoints and progenitor FRT insertions of Exelixis, DrosDel and BSC deletions positioned relative to the proximal and distal ends of annotated genes. It also shows the breakpoints of other deletions in the BDSC Deficiency Kit. The filled-in cells in columns to the right depict deletions graphically. Similarly, the positions of haploinsufficient genes are highlighted in colored rows to show their relationships to the deletions. Notes explain how we estimated the positions of breakpoints for deletions in the Deficiency Kit that have not been molecularly characterized. (XLSX 8 MB)

13059_2012_2771_MOESM3_ESM.XLSX

Additional file 3: Haploinsufficient genes with viability and fertility effects. This table gives details of the mapping and identification of haplolethal and haplosterile genes in D. melanogaster, including cytological locations, hemizygous phenotypes, deletions removing the genes, deletions most closely flanking the genes and the number of closely linked genes sharing the same deletion breakpoint interval. (XLSX 42 KB)

13059_2012_2771_MOESM4_ESM.XLSX

Additional file 4: Bloomington DrosophilaStock Center Deficiency Kit stocks. This file lists the component deletions in the BDSC Deficiency Kit with their breakpoints and stock information. (XLSX 47 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and permissions

About this article

Cite this article

Cook, R.K., Christensen, S.J., Deal, J.A. et al. The generation of chromosomal deletions to provide extensive coverage and subdivision of the Drosophila melanogaster genome. Genome Biol 13, R21 (2012). https://doi.org/10.1186/gb-2012-13-3-r21

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2012-13-3-r21

Keywords