Skip to main content

Widespread duplications in the genomes of laboratory stocks of Dictyostelium discoideum

Abstract

Background

Duplications of stretches of the genome are an important source of individual genetic variation, but their unrecognized presence in laboratory organisms would be a confounding variable for genetic analysis.

Results

We report here that duplications of 15 kb or more are common in the genome of the social amoeba Dictyostelium discoideum. Most stocks of the axenic 'workhorse' strains Ax2 and Ax3/4 obtained from different laboratories can be expected to carry different duplications. The auxotrophic strains DH1 and JH10 also bear previously unreported duplications. Strain Ax3/4 is known to carry a large duplication on chromosome 2 and this structure shows evidence of continuing instability; we find a further variable duplication on chromosome 5. These duplications are lacking in Ax2, which has instead a small duplication on chromosome 1. Stocks of the type isolate NC4 are similarly variable, though we have identified some approximating the assumed ancestral genotype. More recent wild-type isolates are almost without large duplications, but we can identify small deletions or regions of high divergence, possibly reflecting responses to local selective pressures. Duplications are scattered through most of the genome, and can be stable enough to reconstruct genealogies spanning decades of the history of the NC4 lineage. The expression level of many duplicated genes is increased with dosage, but for others it appears that some form of dosage compensation occurs.

Conclusion

The genetic variation described here must underlie some of the phenotypic variation observed between strains from different laboratories. We suggest courses of action to alleviate the problem.

Background

Genetic variation within a given species can extend from simple polymorphisms at single nucleotides to translocations, inversions and duplications affecting many genes. Recent work shows that such large-scale structural variation may be much more important than previously thought: for instance, the genomes of healthy human individuals may differ in copy number at hundreds of loci, that is, have distinct amplifications and deletions detectable by DNA microarray hybridizations [1–3]. These structural variations can have marked effects on phenotype as demonstrated by their association with pathologies of various kinds [4]. For instance, amplifications of alpha-synuclein cause a rare class of familial Parkinson's disease [5], and triplication of the trypsinogen locus can cause hereditary pancreatitis [6]. All sequence variation can, in principle, affect the function and regulation of genes and it is now possible to estimate the relative contribution of different kinds of mutation to changes in gene expression [7].

Similar variability can occur in laboratory organisms: inbred mouse strains show widespread copy number variation [8, 9], which can be associated with complex phenotypes [10]. Budding yeast grown for generations in particular culture conditions displayed experimentally induced variations, reproducibly accumulating copy number mutations on certain chromosomes [11]; strains selected to suppress a loss of function mutation develop particular segmental duplications [12]. Spontaneous translocations have also been observed genetically in Aspergillus nidulans [13] and Neurospora crassa [14, 15].

Copy number can influence phenotype through a proportional effect on mRNA abundance: aneuploidy, associated with direct increases in gene expression, is implicated in the antifungal drug resistance of certain Candida albicans strains [16]. These effects can also be of pathological significance: for instance, DNA copy number alteration is associated for many genes with altered gene expression in breast tumors [17, 18] and progression of colorectal cancer coincides with large scale changes on copy number that are broadly mirrored by similar changes in mRNA level of affected genes [19].

Dictyostelium discoideum is a widely used laboratory organism, particularly useful for examining problems in cell biology, developmental signaling, the evolution of altruism and the function of conserved genes [20, 21]. The organism grows as singled-celled amoebae, feeding on bacteria, and enters a multi-cellular stage when starved, to eventually produce a stalked fruiting body with a head of viable spores. Virtually all laboratory strains derive from the original type isolate from North Carolina, NC4, dating from 1933. Around 1970 two independent axenic strains - Ax2 and Ax3 - able to grow in complex media, were selected from NC4 [22, 23]. These and their descendents now form the great majority of strains in current use.

Dictyostelium cells can be maintained as vegetatively growing amoebae or stored over long periods either frozen, or as spores. Although a sexual cycle via macrocyst formation exists, it has not been used as a laboratory tool [24, 25]. Genetic exchanges are possible by a parasexual cycle, but are largely limited to chromosomal re-assortments with only a low frequency of recombination [26]. Today this cycle is not widely exploited. Most laboratory stocks therefore represent individual lineages that have become isolated from each other at various times in the past, and which may potentially have diverged from each other over time.

The published genome sequence of the Ax4 strain contains a large inverted segmental duplication on one chromosome [21, 27], which is absent in other lines, notably the type strain NC4, from which Ax4 is ultimately derived. Other genetically marked strains have also been reported to contain duplicated chromosomes, or chromosome segments [28–30] and there are cases where duplicated genes are reported in particular stocks [31, 32], but are only present as single copies in the sequenced genome. Pulsed field gel electrophoresis has also evidenced differences in chromosome size and number between certain strains [33].

These variations are of major practical importance to investigators, especially when they remain unknown, causing phenotypic differences between strains, and difficulties in genetic manipulation. We have surveyed a range of Dictyostelium laboratory strains and wild isolates using array comparative genomic hybridization and find that duplications are unfortunately widespread, such that the same strains, sourced from different laboratories, often differ substantially.

Results

Virtually all laboratory strains of D. discoideum derive from the original type isolate, NC4 [34], with only limited use being made of other wild isolates, such as V12. The axenic strains Ax2 and Ax3 are the most widely used and a particular lineage of Ax3, termed Ax4, has been fully sequenced [21]. A simplified family tree of this lineage is shown in Figure 1a. Axenic strains differ substantially from their parental NC4 stock: they grow more slowly on bacteria and produce smaller fruiting bodies, as is readily apparent from their plaque morphologies (Figure 1b,c). Amplifications and deletions (copy number variation) could be one source of this between-strain variability, in addition to small-scale mutation of individual genes and promoters.

Figure 1
figure 1

Relationships between the most commonly used Dictyostelium strains. (a) Simplified genealogical tree showing the relationships between common laboratory strains derived from NC4. The branch marked 'Ax3' is more complex than shown here: sub-lineages have been given the names KAx3 and Ax4. The auxotrophic strain DH1 was engineered in an 'Ax3' background, and JH10 from 'Ax4.' (b) Plaque morphologies. Cells were plated clonally in association with Klebsiella aerogenes on SM agar. Plaques were photographed after 4 days. Small DH1 plaques are indicated with arrowheads. Variation in diameter is a function of the rate of feeding and of the motility of the amoebae. Where the bacteria are cleared the amoebae aggregate in streams; this process had not yet begun in the slow-growing DH1 plaques. (c) Fruiting bodies. Wild type cells - in this instance NC4(Dee) - form larger, more robust fruiting bodies than axenic mutants.

To assess this potential source of variation, we used a custom-built DNA microarray to perform array comparative genome hybridization. In this procedure, DNA from a strain of interest and a reference strain is labeled with different dyes and the mixture hybridized to the array; after background subtraction the ratio of fluorescent signals gives the relative abundance of the DNA, which we normalize to 1 over the whole genome (log2 ratio of zero). Duplications should give a log2 ratio of 1 and deletions a large negative log2 ratio. In practice, cross-hybridization produces smaller than theoretically expected log2 ratios. Duplications can only be mapped to the nearest array marker, which average roughly 4 kb apart, and the procedure gives no information on chromosomal location of the duplication; their size is given as that of the region duplicated (thus the known duplication on chromosome 2 of Ax3 is referred to as 750 kb, not 1.5 Mb). The reference strain throughout was our version of Ax2 - called Ax2(Ka) - and other stocks were from the Dictyostelium Stock Center [35] or had been received into our laboratory in the past (Table 1).

Table 1 Strains used in this work

Duplications are frequent in laboratory stocks

We examined 11 examples of the Ax2, Ax3, and Ax4 axenic strains. As expected, all Ax3/4 strains share the known chromosome 2 duplication (Figure S1 in Additional data file 1) and we also identified a small duplication/amplification on chromosome 1, common to all Ax2 strains, as described below. Apart from this, 9 of the 11 strains possessed additional duplications, some of which are shared between several lines, indicating clear patterns of relationship. Selected duplications are shown in Figure 2; the sizes and locations of all are given in Table 2, and chromosomal locations are displayed schematically in Figure 3.

Table 2 Chromosomal locations of duplications and their distribution among strains
Figure 2
figure 2

Duplications are frequent in 'wild type' axenic strains. (a-e) Log2 ratios (each strain compared to the Ax2(Ka) reference) are indicated by vertical lines; array probes are ordered according to their chromosomal location given by dictyBase assembly version 2.5. The previously known Ax3 duplication involves the region of chromosome 2 between approximately 2.25 and 3 Mb, which is wholly contained within the region duplicated in Ax2(Wee).

Figure 3
figure 3

The distribution of amplifications across the genome. For each chromosome (depicted as arrows, with scale indicating Mb of sequence), different colored bars represent the segments duplicated, approximately to scale. Each feature is named according to the first column of Table 2, in which more precise data concerning size and location are given, along with the strains involved.

Four of the eleven strains carry unique duplications. Ax2(I) and KAx3(U) have duplications of parts of chromosome 1, of 274 and 62 kb, respectively (Figure 2a,b). Ax2(Wee) and Ax2-206 (a rarely used Ax2 clone from the Gerisch laboratory) bear larger 1,179 and 621 kb non-overlapping duplications from chromosome 2 (Figure 2c,d). The Ax2(Wee) duplication encompasses the Ax3 common duplication, plus around 400 kb to one side of it. This region is probably a hotspot, as three further, independent, duplications have been observed from expression profiling experiments comparing mutant with other strains (unpublished results). Ax2-206 also carries another large duplication of part of chromosome 6 (Figure 2e), within a larger region of log2 ratios greater than zero, but averaging less than we typically observe for regions present in two copies per genome. Ax2-214 (the standard Gerisch stock) and Ax2(M), ultimately deriving from the same laboratory, share a feature in the same region duplicated in Ax2(Wee) and Ax3 (Table 2). Log2 ratios in this feature are clearly shifted away from zero, but average less than 0.2. The basis of these 'sub-duplication' features is not known.

The auxotrophic mutant strains JH10 [36] and DH1 [37] - used as parental strains in molecular genetic studies - also show novel duplications. JH10 carries a unique 129 kb duplication of a segment of chromosome 2 (Table 2), while DH1 has two duplications, both shared with its parent Ax3(Dev) (Table 2, and below).

Ax2 has a small duplication/amplification

A segment of 11 genes on chromosome 1 is under-represented in all strains tested compared to Ax2(Ka) (log2 ratio between -0.5 and -1) except for other examples of Ax2 (Figure 4a). This is most easily explained by an approximately 26 kb amplification common to all Ax2 lines, which presumably occurred when the original strain was selected, analogous to the much larger Ax3 duplication. Ax2(Ka) appears to have two copies of this sequence, but all the other Ax2s tested show an increase compared to Ax2(Ka), indicating three or more copies. The approximate breakpoints of this feature were confirmed by quantitative real-time PCR (Figure 4b). The genes amplified in Ax2 strains are listed in Table S1 in Additional data file 4; notably, there are three protein kinases, as well as a formin and a potential transcription factor.

Figure 4
figure 4

A duplication common to Ax2 lines. (a) All Ax2 strains in our study plus selected other strains of NC4 and non-NC4 backgrounds are displayed. Each block is colored according to the log2 ratio for the comparisons of each strain with reference Ax2(Ka). Since log2 ratios are consistently greater than zero for the duplicated genes in examples of Ax2 other than the reference, we suggest that this region is amplified further in these strains. The genes plotted are: a, DDB0190411; b, DDB0190412; c, DDB0190413; d, DDB0201787 (probe 1); e, DDB0201787 (probe 2); f, DDB0190415; g, DDB0190416; h, DDB0201789; i, DDB0190418; j, DDB0216669; k, DDB0190421; l, DDB0190422; m, DDB0190424; n, DDB0190426; and o, DDB0190427. (b) The breakpoints of the duplication in Ax2(Ka) were confirmed by real-time quantitative PCR, in comparison with Ax4(Ku). Mean log2 ratios ± standard error are shown, summarizing, per gene, four pairwise comparisons of threshold cycles.

A segment of chromosome 5 is often duplicated in the Ax3 lineage

Seven strains descending from Ax3 share a small duplication of chromosome 5 sequence, including Ax3(Dev) and its offspring DH1, as mentioned above (Figure 5). The duplicated genes are listed in Table S2 in Additional data file 4. Also among this group are the parasexually derived strains XP99, NP81, and the latter's offspring HU32, which all derive some, but not all, of their chromosomes from Ax3 (Figure S2 in Additional data file 2). Curiously, this feature is present in Ax4(F) but absent in that strain's presumed offspring JH10, though these two strains are clearly related because they share a three gene deletion not observed in any other strain (see below). The chromosome 5 duplication is also absent in several other examples of the Ax3 lineage, notably Ax4(Ku). It seems that this duplication must have arisen in the Ax3 lineage, but is relatively unstable, and has been independently lost at least twice (this seems a more likely explanation than the possibility of separate duplication events in, and only in, the Ax3 lineage).

Figure 5
figure 5

A novel duplication present in a subset of the Ax3 lineage. Nine strains in our study are lineal descendants of Ax3, and one other carries one or more chromosomes from it. NC4A2, based on our evidence, also descends from Ax3. Of these 12 lines, 7 carry a near identical duplication of chromosome 5 sequence. The breakpoints are not entirely clear because of noise in the data, and it is possible that there is some difference between strains. The genes plotted here are: a, DDB0188657; b, DDB0219507; c, DDB0188659; d, DDB0188660; e, DDB0188661; f, DDB0188665; g, DDB0188667; h, DDB0216146; i, DDB0188671; j, DDB0188673; k, DDB0188674; l, DDB0188677; m, DDB0188678; n, DDB0188686; o, DDB0188687; and p, DDB0188688.

Strains used in parasexual genetics

Haploid Dictyostelium cells occasionally fuse to make fairly stable diploids, which can break down by random chromosome loss to reform haploids with a re-assorted chromosome complement. By selecting for diploid formation and breakdown, a workable parasexual system was developed for complementation testing and assigning markers to linkage groups [38]. However, this system sometimes produced anomalous results, to which unrecognized duplications might have contributed [39]. We therefore examined a number of strains dating from this parasexual era.

The most complicated pattern we have seen is given by NP81 and its offspring HU32. As well as multiple duplications, they also possess many contiguous regions of apparent gene loss (an example chromosome of each strain is shown in Figure S3 in Additional data file 3; all chromosomes show some stretches of gene loss). The log2 ratios in these regions are not extreme enough to suggest complete absence of the sequences, and in any case this is unlikely, given the likely presence of essential genes in these regions. They cannot represent duplications in the reference genome because the same DNA sample was used as reference in all hybridizations. We tentatively propose that these strains are degenerate diploids, hemizygous at the regions of apparent gene loss.

NP81 was selected for growth in the presence of the DNA damaging agent ethidium bromide [40] so it is not entirely surprising for its genome to show multiple abnormalities. In contrast, none of X22 [41], XP55 [42] and XP99 [43], which are derived from heavily mutagenized strains but not selected on ethidium bromide, show aberrations similar to NP81. There are no duplications discernible using our arrays in XP55 and X22, although XP99 has a unique one involving chromosome 6, as well as the smaller chromosome 5 feature it inherited from Ax3. The data for XP55 and X22 suggest that the once-standard methods of chemical mutagenesis and parasexual manipulation do not necessarily induce duplications at high frequency.

NC4A2 carries a duplication indistinguishable from the chromosome 2 duplication common to all Ax3 strains

NC4A2 is an axenic strain claimed to be directly selected from NC4, and in consequence, to have superior properties to the standard axenic strains [44]. However, two examples of this strain, obtained from different sources, both carry what appears to be the same chromosome 2 duplication as seen in Ax3 (Figure 6). Although regions of chromosome 2 have been duplicated independently several times, the breakpoints in this case are very similar (or indeed, the same) to those in Ax3; NC4A2 also lacks two other distinct duplications present in its presumed parent, NC4(Kn), as listed in Table 2. Thus, we believe that the strain currently designated as NC4A2 arose from inadvertent contamination by Ax3/Ax4 cells. There have been reports that its properties differ significantly from Ax4 (R Insall, personal communication), but in our hands its growth on bacteria and fruiting body morphology are much more similar to Ax4 than NC4 (not shown).

Figure 6
figure 6

NC4A2 lines contain a duplication of the same segment of chromosome 2 that is duplicated in Ax3. The duplication appears for the most part identical in all strains derived from Ax3. We show here (a) Kax3(U), (b) NC4A2(Kn), and (c) NC4A2(SC) because they display points of similarity not observed in the other examples of this lineage in our study. The point of inversion of this tandem inverse duplication is to the right of the plot, where some genes (log2 ratios negative) appear to have been deleted in both copies in NC4A2 and KAx3(U). At least one of these genes appears to have been lost in both copies in several other of the Ax3-lineage strains in our study, but unfortunately some of the probes for these genes were not printed well and so our data do not permit us to assess exactly how frequent these deletions are. A segment within the duplication towards the left-hand side appears to be present as a single copy in both NC4A2 lines and in KAx3(U); this runs from DDB0233427 to DDB0191242, and appears to be present in the expected two copies in all other Ax3 derived strains we have studied.

NC4A2 appears to be most closely related to KAx3(U), since both these strains have lost a segment of about 29 kb from one half of the inverted duplication on chromosome 2, which is now present as a single copy, and lack the other, novel, Ax3 duplication of chromosome 5 sequence. These and several other strains of the Ax3 lineage appear to have completely lost sequence near the point of inversion of the chromosome 2 duplication. The open reading frame designated DDB0217158 [45] is especially unstable. This mirrored region could be a target for recombination, leading to excision of segments. It is possible that the sequence of this region in Ax4(Ku), although apparently more complete than in some of its relatives, has also degenerated in the same way, resulting in the complete loss of some of the ancestral sequence.

Duplications are also frequent in different stocks of NC4

To test whether duplications are a peculiarity of axenically maintained stocks, we examined a number of stocks of NC4, their non-axenic parent. We particularly sought lines of known history: for instance, NC4(S) came from a vial of spores lyophilized in the Raper laboratory in 1969, which was finally opened in the Schaap laboratory in 2005 (P Schaap, personal communication) and NC4(L) came directly from Raper, but was received in the Loomis laboratory after the generation of Ax3 (W Loomis, personal communication). We were surprised to find that most of the NC4 lines also contain duplications, which predominate in the sub-telomeric regions of the chromosomes (Figure 7 and Table 2). Again, most duplications differed in location in different lines, the exception being NC4(Kn), a stock of NC4(L) taken by D Knecht when he left the Loomis laboratory. This retains the same duplication as NC4(L), without gaining any further duplications, showing both that this duplication arose early and that duplications are not necessarily common. This duplication had been previously detected by 'mapping using haploid amounts of DNA and the polymerase chain reaction' (HAPPY mapping) - the strain is just referred to as NC4 in the paper [21] - but our estimate of its length at 495 kb is larger than the earlier rough estimate of 300 kb.

Figure 7
figure 7

Duplications are also frequent in non-axenic wild types in the NC4 lineage. (a) NC4(L) chromosome (chr) 2, (b) NC4(B) chromosome 4, (c) NC4(S) chromosome 6. As well as the three strains shown, NC4(Wi) has a duplication overlapping that observed in NC4(S). It is not clear why duplications near the chromosome ends are especially frequent in these non-axenic lines, or indeed whether this is merely a sampling artifact. The parasexual segregant XP99 is not axenic and yet carries two duplications away from termini; however this strain inherited some of its chromosomes ultimately from Ax3 and we can be confident that one of these duplications occurred in that axenic ancestor.

Since these duplications differ from stock to stock, we assumed that the original NC4 isolate lacked all of them, and therefore attempted to recover an NC4 strain of this genomic structure. Finally, we found three lines, DdB(SC), DdB(Wel), and NC4(Dee), which are without any discernible duplication, though they do lack a small duplication believed to be present in the founding NC4 stock (see below). DdB is a clone of NC4 that was selected in the laboratory of M Sussman, and NC4(Dee) was obtained by R Deering in the late 1960s from Sussman, then maintained in his laboratory, before transfer to D Welker in 1977 (D Welker, personal communication).

Duplications in other wild isolates

The unexpected prevalence of duplications even among different stocks of NC4 might imply that the Dictyostelium genome is inherently unstable, or alternatively that instability is a consequence of laboratory culture. To examine this question we tested a number of other wild, little-cultured lines, including recent isolates made by D Francis at the site of the original type isolate at Little Butts Gap, North Carolina [46]. Only one of these seven strains shows evidence of a large duplication similar to those observed in laboratory strains (Table 2).

Two proximal derivatives of V12, another isolate from the wild that has been used as a standard non-axenic strain, were tested and found to be without such amplification: V12M2 is a clone of V12 chosen by G Gerisch and used for stalk cell inductions [47] and NP73 is an axenic derivative of V12 selected by K Williams (not shown). Two other wild strains, NYA64 and WS205, and a cycloheximide-resistant mutant derived from another wild isolate (A2cycR) also lack detectable duplications.

Most wild isolates have a two-gene duplication that has been lost in all axenic strains

Small duplications are difficult to distinguish from experimental noise at the level of replication used in this study. However, when present in a large enough portion of the sample they can still be reliably discerned. The notable example we found concerns two genes on chromosome 3. Remarkably, this duplication is found in all non-NC4 wild isolates tested (and A2cycR, a mutant derived from a Wisconsin wild isolate) and a subset of NC4 lines, including the mutant X22 (Figure 8). It is absent in NC4(Dee), the two DdB lines, NC4(B), two of the genetically marked non-axenic strains (XP55 and XP99), and all axenic lines tested. Note that the duplication is absent in NC4A2 but not its supposed parent NC4(Kn). Since it is extraordinarily unlikely that this clear division is the result of independent duplications in many different wild strains from different locations, we infer that this duplication was present in the original NC4 isolate, but a copy was lost in the DdB lineage of NC4 in the Sussman laboratory, and hence in all its descendants. The two duplicated genes encode a putative GATA transcription factor and a protein serine/threonine kinase; in between them is a small open reading frame that would encode a 46 amino acid polypeptide, if expressed.

Figure 8
figure 8

Two adjacent genes are duplicated in wild isolates but not in a subset of laboratory stocks, including all axenic strains. Shown here are all wild isolates in our study, plus selected laboratory strains descended from NC4. Note that some NC4 derived strains group with the wild strains, having two copies of these genes, but a subset (including all axenic strains) possess them in single copy.

Deleted or diverged genes

Since essential genes are likely to be interspersed at regular intervals along all chromosomes, deletions are likely to encompass only a few genes and, therefore, will be more difficult to detect than large duplications. Potential deleted (or highly diverged) genes, compared to reference Ax2(Ka), were taken as those with a log2 ratio of less than -3. Their occurrence, in all wild isolates and in those laboratory strains in which they occur, is plotted in Figure 9. Of the very few such genes among NC4-derived strains, one is the engineered deletion of the pyr5-6 locus in the uridine auxotrophic strain DH1. Also notable is an apparent deletion on chromosome 6 shared by Ax4(F) and JH10 (which are from the same laboratory, with JH10 deriving from Ax4(F)). The three genes affected all show homology to known genes in other organisms.

Figure 9
figure 9

Potentially deleted sequences. Probes with extreme log2 ratios (below -3 in any strain) are considered to reflect deletions in that strain. The data have been clustered both by gene and by strain, to display relationships better (but information about groupings of genes along chromosomes is scrambled somewhat). The dendrogram gives an indication of relationships between strains, but note that this is based on the extremely small set of 36 genes shown here, and so should not be over-interpreted. The candidate genes are: 1, DDB0205403; 2, DDB0188003; 3, DDB0188004; 4, DDB0168894; 5, DDB0218478; 6, DDB0188007; 7, DDB0188002; 8, DDB0215073; 9, DDB0191949 (probe 1); 10, DDB0191949 (probe 2); 11, DDB0206115; 12, DDB0206106; 13, DDB0206110; 14, DDB0206111; 15, DDB0191930; 16, DDB0206108; 17, DDB0206112; 18, DDB0206109; 19, DDB0167672; 20, DDB0219404; 21, DDB0187848; 22, DDB0185937; 23, DDB0188514; 24, DDB0203140; 25, DDB0186442; 26, DDB0206404; 27, DDB0218143; 28, DDB0219338; 29, DDB0202734; 30, DDB0188991; 31, DDB0217456; 32, DDB0184376; 33, DDB0184375; 34, DDB219746; 35, DDB0206525; 36, DDB0217158.

By far the majority of these putative deletions are from wild isolates other than NC4. These strains tend also to have many measurements between log2 ratios -1 and -3, which likely represent a mixture of deletions and polymorphic sequences; the strain containing the most of these overall is NYA64, a Japanese isolate, and the only one not from the USA in our panel. Although these potentially deleted loci have not been confirmed by other means, they are listed in Table 3 as a resource for defining non-essential genes or (in the wild isolates) particularly divergent loci. None, other than pyr5-6, has been previously characterized well enough to appear in the published literature.

Table 3 Putative deleted genes

Effect of gene dosage on mRNA abundance

The most likely phenotypic consequence of gene duplication is through effects on the relative abundance of mRNAs expressed from genes within the duplication. To assess this we compared gene expression in growing Ax2(Wee) and Ax2(Ka) cells using the same microarray, and examined the relative expression of genes within the large duplication possessed by Ax2(Wee) on chromosome 2 (Figure 10a). In this case, gene dosage clearly affects the abundance of a large subset of vegetative mRNAs within the duplication, with some mRNAs being roughly twice as abundant as in the unduplicated parent, while many others are little altered (Figure 10b). The distribution of log2 ratios within the duplication is non-normal (Shapiro-Wilk test, p = 3.6 × 10-6). Since the amplified region in Ax2(Wee) includes the entire region duplicated in the Ax3 lineage, we could also compare these data with an expression experiment comparing DH1 with Ax2 (unpublished work with T Soldati). The log2 ratios of the two comparisons correlated well (Figure 10c), and a high proportion of the little-altered genes were shared between the two cases (Fisher's exact test, p = 2.7 × 10-12 using a simple cut off of absolute log2 ratio less than 0.25; the p-value is similar if a cut off of 0.5 is used). Having excluded low intensity spots, and there being no obvious reason why these probes should have grossly underestimated differences in expression - for instance, there is very little correlation between mean log2 ratio and mean log intensity - we propose that, unlike budding yeast [12], Dictyostelium must possess modes of regulation that counteract increased gene dosage.

Figure 10
figure 10

Duplications are apparent also at the mRNA level. (a) Data from an experiment comparing mRNA from growing Ax2(Ka) and Ax2(Wee) ordered by chromosomal location, for chromosome 2 only, where Ax2(Wee) has a duplication. This is the only region displaying a striking shift away from log2 ratio = 0. (b) Histogram of log2 ratios within the region duplicated in Ax2(Wee). The distribution tends towards bimodality, with one clear peak near log2 ratio = 0 and a less distinct one towards log2 ratio = 1. Thus, a portion of genes is dosage-insensitive. (c) Correlation between log2 ratios from independent mRNA comparisons. The duplication carried by Ax2(Wee) entirely encompasses the segment duplicated in Ax3 and its derivatives, allowing us to compare the dosage-sensitivity of genes in this overlap in different strains. The data are clearly correlated, with genes sensitive to dosage in Ax2(Wee) tending also to be sensitive in DH1 (unpublished work with T Soldati; the Soldati laboratory strains are indicated by So). A similar set of genes appears to undergo dosage compensation in both strains.

Discussion

Our results reveal that duplications are common among all Dictyostelium laboratory stocks, except for those that have recently been isolated from the wild. This discovery brings an unwelcome element to phenotypic and molecular genetic experiments, but one that can be ameliorated by suitable precautions. We have also been able to glean useful information about Dictyostelium genealogy, identify some apparently non-essential genes, and loci that vary between wild-type stocks. Unexpectedly, our experiments also suggest that Dictyostelium has a form of gene-dosage compensation that acts on individual genes.

Prevalence and practical consequences of duplications

The frequency of duplications is such that any two laboratory stocks of Ax2 or Ax3/4 or NC4 are more likely than not to differ by at least one duplication or related feature (with a conservative lower limit for detection of 30 kb). These differences range from unique duplications, estimated to occur in 30-40% of laboratory strains of the NC4 lineage, both axenic (5 out of 15 strains used as wild type, including DH1, JH10 and NC4A2) and non-axenic strains (3 of 8 strains of NC4 or DdB) to features common to several stocks of a strain, but differing in others, such as that on chromosome 1 of Ax2. We have not studied transformed strains systematically but in comparisons of mRNA levels between mutants and wild-types, novel duplications can be discerned in 5 out of 30 independent mutant strains (unpublished observations). This is likely to be a minimal estimate since the greater noise in RNA experiments means that smaller duplications will be missed. We have not examined any strains generated by restriction mediated integration (REMI) insertional mutagenesis [48].

The two most commonly used strains - Ax2 and Ax3/4 - were already known to differ by the large, 750 kb duplication on chromosome 2 present in all Ax3 derivatives, but we find two additional copy number differences between the lineages. Ax2 has a duplication of chromosome 1 sequence of about 26 kb and Ax3/4 has a second, less stable duplication on chromosome 5 of about 55 kb. This latter feature is present only in a subset of Ax3 derivatives, notably being absent from the sequenced strain Ax4(Ku), and JH10 despite being present in the latter's parent, Ax4(F). It was almost certainly present in the early history of Ax3, since it is found in strain XP99, which inherited chromosome 5 from NP2, which in turn was selected from Ax3 in the early 1970s [49].

The presence of duplications in the genome brings several experimental problems. The increased gene dosage causes increased RNA levels for at least some of the genes in the duplication, with unpredictable phenotypic consequences. Depending on the mechanism of duplication, there may also be damage at the event boundaries and the duplications themselves may be unstable, giving long-term genetic instability, as seen to some extent with the two duplications in the Ax3/4 lineage.

In addition to this, genetic manipulation of duplicated genes is much more problematic: two hits are required to knock them both out, which is achievable [50], but more difficult than when there is only a single locus to target. Duplicated regions would be expected to be relatively resistant to mutagenesis, since for recessive mutations, independent hits in both copies of the duplicated gene would be required to see the full phenotype. Effectively, this makes duplicated genes almost completely recalcitrant to genetic screens based on REMI insertional mutagenesis [48], which works at a frequency of roughly one hit per genome.

Amelioration

It is clear that with good practice - strains stored as spores or frozen cells and working stocks renewed frequently - Dictyostelium stocks can be kept for many years without any detectable change. For instance, Ax2-214 stocks from the Gerisch and Malchow laboratories do not seem to have diverged despite a separation of more than 20 years, and nor does NC4(L) from the Loomis and Knecht laboratories. Although making a knock-out mutant may carry some risk of producing a duplication, the standard practice of comparing several independent clones, and rejecting any that are different, should greatly reduce the possibility of studying a phenotype that is due to secondary changes. Even so, it seems from the occurrence of duplications among knock-out stocks that it would be desirable to screen any new strain for copy number changes.

What is clearer than before is that it is risky to compare a mutant with anything other than its direct parent, or another mutant from the same parent. It was already known that the parental stocks Ax2 and Ax3/4 differ in fine detail when high-resolution phenotypic assays are applied, such as those for chemotaxis or patterns of developmental gene expression [51, 52] (see the discussion of strain history at dictyBase for examples [53]), but it is now apparent that differences will exist even between stocks of the same strain. It would be advantageous in future to focus on just a few parental stocks, preferably those with a minimum of duplications, thus allowing wider phenotypic comparisons of mutants.

Genealogy

Because the duplications and other features that we detect are quite stable over time, we can reconstruct some genealogies. The family tree of strains descended from the type isolate is given in Figure S2 of Additional data file 2. Our results broadly agree with the published histories of these strains, the one anomaly being NC4-A2 [44], which may originally have been an independent axenic isolate selected from NC4, but at some stage appears to have been contaminated and replaced by a line indistinguishable from Ax3.

We can also shed a little more light on the prehistory of the axenic strains. Most examples of NC4 and all the wild isolates we examined have a duplication of two genes on chromosome 3, which is missing in Ax2 and Ax3/4 and all their axenic derivatives - they have single copies of these two genes. Since Ax2 and Ax3 are of independent origin, they most likely share a common ancestor that also lacks the duplication. This is the case: strain DdB, their reported ancestor (see the discussion of strain history at dictyBase [45]; W Loomis, personal communication), has only single copies of these genes. DdB (also known as D. discoideum strain B, or NC4 strain B, not to be confused with the NC4(B) of this study) was a clone selected from NC4 in the laboratory of Maurice Sussman by 1967 at the latest [54, 55]. Assuming that other strains lacking this small duplication also descend from the same clone, then members of the DdB lineage in our study include NC4 (Dee), which was obtained from the Sussman laboratory in 1968 or earlier [56]; the axenic strains NP81 and HU32, which inherited chromosome 3 ultimately from Ax3 [57]; XP55 and XP99, both of which inherited chromosome 3 (bearing the bsgA mutation) from NP194. NP194 is the offspring of NP20, a mutant that arose spontaneously from the DdB stock in the laboratory of P Newell [55].

We can only infer the original genomic structure of the NC4 type isolate. It should have lacked the sporadic duplications found in several existing stocks, such as NC4(S) and NC4(L), but had the small duplication on chromosome 3 discussed above. Unfortunately, we have been unable to recover an NC4 stock with exactly these characteristics. DdB, originating from the Sussman laboratory, may be the closest approximation, lacking only the small duplication.

Mechanisms of duplication

Our method gives no information about duplication structure, and so we cannot even be sure that duplicates are on the same chromosome; experimental noise and the spacing of microarray probes also makes the precise determination of the breakpoints difficult. Both these factors limit our ability to identify possible mechanisms by which duplications arise.

However, the terminal position of several duplications suggests that these ones could have arisen from translocations. Some duplications appear to share one end to within a few kilobases, suggesting that certain regions are particularly prone to these events, presumably arising through the same mechanism. On chromosome two, duplications in JH10 and Ax2-206 both end near nucleotide 6,752,000, close to a complex repeat region associated with TRE3 repeats and two tRNA genes. The other breakpoint in Ax2-206 on this chromosome has a similar repeat. Therefore, a second possible mechanism is that the complex repeats found in scattered clusters throughout the genome [21, 58] might promote unequal crossing-over events, thus duplicating segments of DNA between them [59]. Consistent with this, the Ax2(Wee) duplication is flanked by DIRS-1 sequences, and a further seven duplications have repeat sequences within 10 kb of one breakpoint (duplications 2G, 2I, 4B, 6A, 6D, 6F and 6H in Table 2). A subset of retrotransposons insert preferentially next to tRNA genes [58] and the chromosome six duplications in Ax2-206 and NP81 both end near nucleotide 3,055,000, which is near a tRNA gene; and another duplication in NP81 on the same chromosome (6A) is also flanked by tRNAs. It is possible that these strains have repeats inserted at the tRNAs that are not present in the Ax4(Ku) genome. For the other duplications, no complex repeats or tRNA genes could be identified near either breakpoint, so other mechanisms must be considered for these events. Detailed investigation of the breakpoint sequences will be required to address these issues.

We presume that most duplications become fixed when populations are established from single cells, although increased copy number of certain loci may also confer a selective advantage in laboratory conditions. Their rarity in the wild, as judged from recent isolates, suggests there may be negative selection pressures against duplications, confirming indirectly that they are likely to have phenotypic consequences. A genetically-inferred duplication on chromosome 3 was suggested to be detrimental to growth [29], so it is notable that the two duplications on this chromosome in the present study are very small, of approximately 4 and 17 kb in length.

Deletions

Because essential genes are interspersed along all chromosomes, large deletions should be lethal, and indeed none are observed. A small number of NC4-derived strains have apparent deletions of at most three adjacent genes. This includes the engineered deletion of pyr5-6 in DH1 [37] but not of thyA in JH10, which is disrupted by an insertion rather than a deletion [36].

In wild isolates, apparent gene loss reflects a spectrum from sequence divergence to actual gene loss. More genes in wild isolates (other than NC4) pass our cutoff for putative deletions, but nevertheless the overall difference in gene repertoire is small, being at most 0.2% missing genes compared to NC4 (using a cutoff of log2 ratio < -3). As the array is effectively based on NC4, we cannot detect genes present in other isolates but not in NC4. Hints of geographical relationship are apparent: isolates from the type location cluster together most closely, and NYA64 from Japan (the only non-American isolate) is the most diverged when the log2 ratio cutoff is relaxed to allow detection of diverged sequence, rather than just those that are absent entirely.

Gene dosage affects mRNA abundance, but not simply

In metazoans there is clear evidence that gene dosage effects feed through to the phenotype of the organism, so much so that compensatory mechanisms exist to counteract them. However, this may not be the case in yeast, as strains of budding yeast either aneuploid for entire chromosomes or carrying a smaller segmental duplication show little or no detectable dosage compensation [12]. We find that mRNA abundance is increased over a duplication, but varies from gene to gene, with some showing little increase. This variation does not appear to be due just to experimental noise, especially as there is a very good correlation between expression levels over two, independent, duplications spanning a common region. This might reflect incidental effects of differing modes of regulation of these genes - for instance, an activating factor could be limiting, or a repressor superabundant - but conceivably compensatory mechanisms to counteract dosage increase might exist for some genes, especially since duplication events seem to be so frequent.

Conclusion

Our results show that previously unrecognized duplications are common in laboratory stocks of D. discoideum. These duplications are likely to alter the phenotype of the cells carrying them, and to cause problems, both when detailed phenotypes are examined, and in genetic manipulation. Nevertheless, stocks are relatively stable, and have been maintained over many years without apparent change, showing that with appropriate care, the problems caused by gratuitous duplications can be minimized. We also discerned possible genetic differences between wild isolates, and provide a list of potentially divergent loci.

Materials and methods

D. discoideum strains were grown on SM agar plates in association with Klebsiella aerogenes, or else in axenic growth medium as previously described [60]. Strains were either obtained from the Dictyostelium Stock Centre or received into this laboratory as noted in Table 1. Different examples of the same strain are distinguished by letter codes indicating laboratories of origin. All strains were stored, without cloning, immediately upon receipt, either as spores (4°C or -20°C on silica gel with renewal times of 5-15 years) or frozen in liquid nitrogen (indefinitely viable).

To extract genomic DNA, cells were starved over-night and resuspended in lysis buffer (20 mM Tris-HCl, 5 mM MgCl2, 0.32 M sucrose, 0.02% sodium azide, 1% Triton X-100, pH 7.4) at 4°C, vortexed and incubated at 4°C for 15 minutes. Nuclei were pelleted at 3,000 g for 10 minutes, resuspended in lysis buffer and pelleted again, before freezing the pellets on dry ice. Proteinase K (100 μl (20 mg/ml)) was added, followed immediately by 10 ml digestion buffer (10 mM Tris-HCl, 5 mM EDTA, 0.7% EDTA, pH 7.5), and the pellet resuspended by gentle trituration. The lysate was incubated for 1 h at 60°C and the DNA phenol-chloroform extracted. DNA (5 μg) was labeled with Cy3 or Cy5 by incorporation of dye-conjugated dCTP (GE Healthcare, Little Chalfont, UK) from random hexamer primed strand synthesis catalyzed by the Klenow fragment (Invitrogen, Paisley, UK). Labeled samples were paired with an Ax2 (Ka) sample labeled with the complementary dye and hybridized to custom DNA microarrays printed on Codelink slides (GE Healthcare). Between two and four replicate hybridizations, always including both dye orientations, were carried out for each strain.

RNA was extracted using RNeasy kits (Qiagen, Crawley, UK). Samples were labeled and hybridized essentially similarly except that 25 μg of total RNA was labeled using anchored oligo(dT) primers and Superscript III reverse transcriptase (Invitrogen). Five independent replicates were carried out, two in one dye orientation, and three in the other.

The DNA microarray consisted of 9,247 PCR products non-redundantly representing 8,579 genes, with little bias regarding chromosomal location. The design was based on the sequence of strain Ax4(Ku), and genomic DNA of this strain was used as the template for PCR amplifications to generate the probes. The PCR products are predominantly between 200 and 400 bp in length (with a small fraction between 150 and 200 bp) and located towards the 3' end of predicted genes. Although overall coverage is good (8,579 genes non-redundantly represented of 10,500-12,500 predicted), non-coding DNA, including ncRNA genes, is not covered at all (unless accidentally), and some regions are more sparsely covered than others. Probes are ordered according to dictyBase assembly version 2.5 [45]).

Microarrays were hybridized for 16-18 h at 42°C, washed at room temperature, and scanned using a Genepix 4000B scanner (Molecular Devices, Wokingham, UK). Images were quantified using Genepix 3 software, then background-subtracted using the Kooperberg function, normalized with the print-tip loess algorithm, and model-fitted using Limma [61–63]. Data from low intensity probes (defined as mean log2 intensity in the hybridizations involving a particular strain less than 6) were omitted; in such a case the median was taken of the log2 ratios of the ten probes either side of it as they are arranged chromosomally (except those of mean log2 intensity less than 6). For probes within 10 of the end of a chromosome a median of the 20 terminal probes was taken (except those of mean log2 intensity less than 6).

The known 750 kb duplication of the sequenced strain Ax4(Ku) was readily apparent (Figure S1 in Additional data file 1), standing out on average 5.4 standard deviations beyond the remainder of chromosome 2 probes. The mean log2 ratio of probes within the duplication is 0.60. A similar underestimation relative to the expected value of one for a duplication in a haploid strain was observed for many different duplications, and presumably reflects a hybridization or scanning artifact, perhaps caused by cross-hybridization of equal copy-number probes. Noise in the data was approximately the same in hybridizations of the same DNA labeled in both channels as in comparisons of different strains (as assessed by overall standard deviation) and so can be explained mostly by experimental error, and not, for example, by sequence variation.

Novel duplications were identified by direct inspection of the data, and on occasion rough breakpoint positions were estimated using DNAcopy [64]. Complex repeat sequences were identified by performing BLAST alignments [65] of the 20 kb sequence surrounding our estimate of the breakpoint position against a database of the sequences described in [66].

Quantitative real-time PCR was carried out using Superscript III and Platinum Taq (Invitrogen) and an Mx3000P thermal cycler (Stratagene, Amsterdam, The Netherlands). Standard curves were made for each primer pair using Ax4(Ku) genomic DNA, and four comparisons made between that template and genomic DNA of Ax2(Ka). The mean of these four log2 ratios thus obtained was then taken.

Array data are deposited at Array Express [67] as array design A-SGRP-3 and the data as submission E-TABM-394.

Additional data files

The following additional data are available with the online version of this paper. Additional data file 1 is Figure S1 showing the Ax3/4 duplication on chromosome 2. Additional data file 2 is Figure S2 showing the genealogy of NC4-derived strains. Additional data file 3 is Figure S3, showing multiple copy number variants on chromosome 5 of strains NP81 and HU32. Additional data file 4 contains Table S1, listing genes present on the Ax2 duplication, Table S2, listing genes on the chromosome 5 duplication present in many Ax3 lines, and the legends for Figures S1, S2 and S3.

Abbreviations

REMI mutagenesis:

REstriction Mediated Integration mutagenesis.

References

  1. Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Maner S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science. 2004, 305: 525-528. 10.1126/science.1098918.

    Article  PubMed  CAS  Google Scholar 

  2. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36: 949-951. 10.1038/ng1416.

    Article  PubMed  CAS  Google Scholar 

  3. Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, Cho EK, Dallaire S, Freeman JL, Gonzalez , Gratacas , Huang J, Kalaitzopoulos D, Komura D, MacDonald JR, Marshall CR, Mei R, Montgomery L, Nishimura K, Okamura K, Shen F, Somerville MJ, Tchinda J, Valsesia A, Woodwark C, Yang F, et al: Global variation in copy number in the human genome. Nature. 2006, 444: 444-454. 10.1038/nature05329.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Estivill X, Armengol L: Copy number variants and common disorders: filling the gaps and exploring complexity in genome-wide association studies. PLoS Genet. 2007, 3: 1787-1799. 10.1371/journal.pgen.0030190.

    Article  PubMed  CAS  Google Scholar 

  5. Singleton AB, Farrer M, Johnson J, Singleton A, Hague S, Kachergus J, Hulihan M, Peuralinna T, Dutra A, Nussbaum R, Lincoln S, Crawley A, Hanson M, Maraganore D, Adler C, Cookson MR, Muenter M, Baptista M, Miller D, Blancato J, Hardy J, Gwinn-Hardy K: alpha-Synuclein locus triplication causes Parkinson's disease. Science. 2003, 302: 841-10.1126/science.1090278.

    Article  PubMed  CAS  Google Scholar 

  6. Le Marechal C, Masson E, Chen JM, Morel F, Ruszniewski P, Levy P, Ferec C: Hereditary pancreatitis caused by triplication of the trypsinogen locus. Nat Genet. 2006, 38: 1372-1374. 10.1038/ng1904.

    Article  PubMed  CAS  Google Scholar 

  7. Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, de Grassi A, Lee C, Tyler-Smith C, Carter N, Scherer SW, Tavare S, Deloukas P, Hurles ME, Dermitzakis ET: Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007, 315: 848-853. 10.1126/science.1136678.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  8. Li J, Jiang T, Mao JH, Balmain A, Peterson L, Harris C, Rao PH, Havlak P, Gibbs R, Cai WW: Genomic segmental polymorphisms in inbred mouse strains. Nat Genet. 2004, 36: 952-954. 10.1038/ng1417.

    Article  PubMed  CAS  Google Scholar 

  9. Adams DJ, Dermitzakis ET, Cox T, Smith J, Davies R, Banerjee R, Bonfield J, Mullikin JC, Chung YJ, Rogers J, Bradley A: Complex haplotypes, copy number polymorphisms and coding variation in two recently divergent mouse strains. Nat Genet. 2005, 37: 532-536. 10.1038/ng1551.

    Article  PubMed  CAS  Google Scholar 

  10. Cutler G, Marshall LA, Chin N, Baribault H, Kassner PD: Significant gene content variation characterizes the genomes of inbred mouse strains. Genome Res. 2007, 17: 1743-1754. 10.1101/gr.6754607.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Dunham MJ, Badrane H, Ferea T, Adams J, Brown PO, Rosenzweig F, Botstein D: Characteristic genome rearrangements in experimental evolution of Saccharomyces cerevisiae. Proc Natl Acad Sci USA. 2002, 99: 16144-16149. 10.1073/pnas.242624799.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  12. Hughes TR, Roberts CJ, Dai H, Jones AR, Meyer MR, Slade D, Burchard J, Dow S, Ward TR, Kidd MJ, Friend SH, Marton MJ: Widespread aneuploidy revealed by DNA microarray expression profiling. Nat Genet. 2000, 25: 333-337. 10.1038/77116.

    Article  PubMed  CAS  Google Scholar 

  13. Kafer E: Translocations in stock strains of Aspergillus nidulans. Genetica. 1962, 33: 59-68. 10.1007/BF01725748.

    Article  PubMed  CAS  Google Scholar 

  14. Perkins DD: Chromosome rearrangements in Neurospora and other filamentous fungi. Adv Genet. 1997, 36: 239-398.

    Article  PubMed  CAS  Google Scholar 

  15. Perkins DD, Kinsey JA, Asch DK, Frederick GD: Chromosome rearrangements recovered following transformation of Neurospora crassa. Genetics. 1993, 134: 729-736.

    PubMed  CAS  PubMed Central  Google Scholar 

  16. Selmecki A, Forche A, Berman J: Aneuploidy and isochromosome formation in drug-resistant Candida albicans. Science. 2006, 313: 367-370. 10.1126/science.1128242.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Pollack JR, Sorlie T, Perou CM, Rees CA, Jeffrey SS, Lonning PE, Tibshirani R, Botstein D, Borresen-Dale AL, Brown PO: Microarray analysis reveals a major direct role of DNA copy number alteration in the transcriptional program of human breast tumors. Proc Natl Acad Sci USA. 2002, 99: 12963-12968. 10.1073/pnas.162471999.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Chin K, DeVries S, Fridlyand J, Spellman PT, Roydasgupta R, Kuo WL, Lapuk A, Neve RM, Qian Z, Ryder T, Chen F, Feiler H, Tokuyasu T, Kingsley C, Dairkee S, Meng Z, Chew K, Pinkel D, Jain A, Ljung BM, Esserman L, Albertson DG, Waldman FM, Gray JW: Genomic and transcriptional aberrations linked to breast cancer pathophysiologies. Cancer Cell. 2006, 10: 529-541. 10.1016/j.ccr.2006.10.009.

    Article  PubMed  CAS  Google Scholar 

  19. Tsafrir D, Bacolod M, Selvanayagam Z, Tsafrir I, Shia J, Zeng Z, Liu H, Krier C, Stengel RF, Barany F, Gerald WL, Paty PB, Domany E, Notterman DA: Relationship of gene expression and chromosomal abnormalities in colorectal cancer. Cancer Res. 2006, 66: 2129-2137. 10.1158/0008-5472.CAN-05-2569.

    Article  PubMed  CAS  Google Scholar 

  20. Kessin RH: Dictyostelium. 2001, Cambridge: Cambridge University Press

    Book  Google Scholar 

  21. Eichinger L, Pachebat JA, Glockner G, Rajandream MA, Sucgang R, Berriman M, Song J, Olsen R, Szafranski K, Xu Q, Tunggal B, Kummerfeld S, Madera M, Konfortov BA, Rivero F, Bankier AT, Lehmann R, Hamlin N, Davies R, Gaudet P, Fey P, Pilcher K, Chen G, Saunders D, Sodergren E, Davis P, Kerhornou A, Nie X, Hall N, Anjard C, et al: The genome of the social amoeba Dictyostelium discoideum. Nature. 2005, 435: 43-57. 10.1038/nature03481.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Watts DJ, Ashworth JM: Growth of myxamoebae of the cellular slime mould Dictyostelium discoideum in axenic culture. Biochem J. 1970, 119: 171-174.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Loomis WF: Sensitivity of Dictyostelium discoideum to nucleic acid analogues. Exp Cell Res. 1971, 64: 484-486. 10.1016/0014-4827(71)90107-8.

    Article  PubMed  CAS  Google Scholar 

  24. Wallace MA, Raper KB: Genetic exchanges in the macrocysts of Dictyostelium discoideum. J Gen Microbiol. 1979, 113: 327-337.

    Article  PubMed  CAS  Google Scholar 

  25. Francis D: High frequency recombination during the sexual cycle of Dictyostelium discoideum. Genetics. 1998, 148: 1829-1832.

    PubMed  CAS  PubMed Central  Google Scholar 

  26. Katz ER, Sussman M: Parasexual recombination in Dictyostelium discoideum: selection of stable diploid heterozygotes and stable haploid segregants. Proc Natl Acad Sci USA. 1972, 69: 495-498. 10.1073/pnas.69.2.495.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  27. Kuspa A, Maghakian D, Bergesch P, Loomis WF: Physical mapping of genes to specific chromosomes in Dictyostelium discoideum. Genomics. 1992, 13: 49-61. 10.1016/0888-7543(92)90201-3.

    Article  PubMed  CAS  Google Scholar 

  28. Williams KL, Robson GE, Welker DL: Chromosome fragments in Dictyostelium discoideum obtained from parasexual crosses between strains of different genetic background. Genetics. 1980, 95: 289-304.

    PubMed  CAS  PubMed Central  Google Scholar 

  29. Welker DL, Williams KL: Chromosome rearrangements in Dictyostelium discoideum. Genetics. 1982, 102: 711-723.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. Welker DL, Metz BA, Williams KL: Recessive lethal mutations and the maintenance of duplication bearing strains of Dictyostelium discoideum. Genetics. 1987, 115: 101-106.

    PubMed  CAS  PubMed Central  Google Scholar 

  31. Gamper M, Howard PK, Hunter T, Firtel RA: Multiple roles of the novel protein tyrosine phosphatase PTP3 during Dictyostelium growth and development. Mol Cell Biol. 1996, 16: 2431-2444.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Myers SA, Han JW, Lee Y, Firtel RA, Chung CY: A Dictyostelium homologue of WASP is required for polarized F-actin assembly during chemotaxis. Mol Biol Cell. 2005, 16: 2191-2206. 10.1091/mbc.E04-09-0844.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Cox EC, Vocke CD, Walter S, Gregg KY, Bain ES: Electrophoretic karyotype for Dictyostelium discoideum. Proc Natl Acad Sci USA. 1990, 87: 8247-8251. 10.1073/pnas.87.21.8247.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Raper KB: Dictyostelium discoideum, a new species of slime mold from decaying forest leaves. J Agr Res. 1935, 50: 135-147.

    Google Scholar 

  35. StockCenter. [http://www.dictybase.org/StockCenter/StockCenter.html]

  36. Hadwiger JA, Firtel RA: Analysis of Galpha4, a G-protein subunit required for multicellular development in Dictyostelium. Genes Devel. 1992, 6: 38-49. 10.1101/gad.6.1.38.

    Article  PubMed  CAS  Google Scholar 

  37. Caterina MJ, Milne JLS, Devreotes PN: Mutation of the third intracellular loop of the cAMP receptor, cAR1, of Dictyostelium yields mutants impaired in multiple signaling pathways. J Biol Chem. 1994, 269: 1523-1532.

    PubMed  CAS  Google Scholar 

  38. Newell PC: Genetics of the cellular slime molds. Annu Rev Genet. 1978, 12: 69-93. 10.1146/annurev.ge.12.120178.000441.

    Article  PubMed  CAS  Google Scholar 

  39. Welker DL, Williams KL: Translocations in Dictyostelium discoideum. Genetics. 1985, 109: 341-364.

    PubMed  CAS  PubMed Central  Google Scholar 

  40. Wright MD, Williams KL, Newell PC: Ethidium bromide resistance: a selective marker located on linkage group IV of Dictyostelium discoideum. J Gen Microbiol. 1977, 102: 423-426.

    Article  CAS  Google Scholar 

  41. Williams KL, Newell PC: A genetic study of aggregation in the cellular slime mould Dictyostelium discoideum using complementation analysis. Genetics. 1976, 82: 287-307.

    PubMed  CAS  PubMed Central  Google Scholar 

  42. Newell PC, Henderson RF, Mosses D, Ratner DI: Sensitivity to Bacillus subtilis: a novel system for selection of heterozygous diploids of Dictyostelium discoideum. J Gen Microbiol. 1977, 100: 207-211.

    Article  Google Scholar 

  43. Ratner DI, Newell PC: Linkage analysis in Dictyostelium discoideum using multiply marked tester strains: establishment of linkage group VII and the reassessment of earlier linkage data. J Gen Microbiol. 1978, 109: 225-236.

    Article  PubMed  CAS  Google Scholar 

  44. Shelden E, Knecht DA: Mutants lacking myosin II cannot resist forces generated during multicellular morphogenesis. J Cell Sci. 1995, 108: 1105-1115.

    PubMed  CAS  Google Scholar 

  45. GeneDB. [http://www.genedb.org/genedb/dicty/]

  46. Francis D, Eisenberg R: Genetic structure of a natural population of Dictyostelium discoideum, a cellular slime mold. Mol Ecol. 1993, 2: 385-392. 10.1111/j.1365-294X.1993.tb00031.x.

    Article  PubMed  CAS  Google Scholar 

  47. Town CD, Gross JD, Kay RR: Cell differentiation without morphogenesis in Dictyostelium discoideum. Nature. 1976, 262: 717-719. 10.1038/262717a0.

    Article  PubMed  CAS  Google Scholar 

  48. Kuspa A, Loomis WF: Tagging developmental genes in Dictyostelium by restriction enzyme-mediated integration of plasmid DNA. Proc Natl Acad Sci USA. 1992, 89: 8803-8807. 10.1073/pnas.89.18.8803.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Kessin RH, Williams KL, Newell PC: Linkage analysis in Dictyostelium discoideum using temperature-sensitive growth mutants selected with bromodeoxyuridine. J Bacteriol. 1974, 119: 776-783.

    PubMed  CAS  PubMed Central  Google Scholar 

  50. Sun TJ, Devreotes PN: Gene targeting of the aggregation stage cAMP receptor cAR1 in Dictyostelium. Genes Dev. 1991, 5: 572-582. 10.1101/gad.5.4.572.

    Article  PubMed  CAS  Google Scholar 

  51. Schilde C, Araki T, Williams H, Harwood A, Williams JG: GSK3 is a multifunctional regulator of Dictyostelium development. Development. 2004, 131: 4555-4565. 10.1242/dev.01330.

    Article  PubMed  CAS  Google Scholar 

  52. Pollitt AY, Blagg SL, Ibarra N, Insall RH: Cell motility and SCAR localisation in axenically growing Dictyostelium cells. Eur J Cell Biol. 2006, 85: 1091-1098. 10.1016/j.ejcb.2006.05.014.

    Article  PubMed  CAS  Google Scholar 

  53. dictyBase. [http://www.dictybase.org/]

  54. Brackenbury RW, Schindler J, Alexander S, Sussman M: A choice of morphogenetic pathways in Dictyostelium discoideum induced by the adenosine analog formycin-B. J Mol Biol. 1974, 90: 529-539. 10.1016/0022-2836(74)90233-2.

    Article  PubMed  CAS  Google Scholar 

  55. Mosses D, Williams KL, Newell PC: The use of mitotic crossing-over for genetic analysis in Dictyostelium discoideum: mapping of linkage group II. J Gen Microbiol. 1975, 90: 247-259.

    Article  Google Scholar 

  56. Deering RA: Dictyostelium discoideum: a gamma-ray resistant organism. Science. 1968, 162: 1289-1290. 10.1126/science.162.3859.1289.

    Article  PubMed  CAS  Google Scholar 

  57. Williams KL, Kessin RH, Newell PC: Parasexual genetics in Dictyostelium discoideum: Mitotic analysis of acriflavin resistance and growth in axenic medium. J Gen Microbiol. 1974, 84: 59-69.

    Article  PubMed  CAS  Google Scholar 

  58. Glockner G, Szafranski K, Winckler T, Dingermann T, Quail MA, Cox E, Eichinger L, Noegel AA, Rosenthal A: The complex repeats of Dictyostelium discoideum. Genome Res. 2001, 11: 585-594. 10.1101/gr.GR-1622RR.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  59. Goldberg ML, Sheen JY, Gehring WJ, Green MM: Unequal crossing-over associated with asymmetrical synapsis between nomadic elements in the Drosophila melanogaster genome. Proc Natl Acad Sci USA. 1983, 80: 5017-5021. 10.1073/pnas.80.16.5017.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  60. Kay RR: Cell differentiation in monolayers and the investigation of slime mold morphogens. Methods Cell Biol. 1987, 28: 433-448.

    Article  PubMed  CAS  Google Scholar 

  61. Kooperberg C, Fazzio TG, Delrow JJ, Tsukiyama T: Improved background correction for spotted DNA microarrays. J Comput Biol. 2002, 9: 55-66. 10.1089/10665270252833190.

    Article  PubMed  CAS  Google Scholar 

  62. Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-

    PubMed  Google Scholar 

  63. Smyth GK: Limma: Linear Models for Microarray Data. 2005, New York: Springer

    Google Scholar 

  64. Venkatraman ES, Olshen AB: A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics. 2007, 23: 657-663. 10.1093/bioinformatics/btl646.

    Article  PubMed  CAS  Google Scholar 

  65. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.

    Article  PubMed  CAS  Google Scholar 

  66. Glockner G, Eichinger L, Szafranski K, Pachebat JA, Bankier AT, Dear PH, Lehmann R, Baumgart C, Parra G, Abril JF, Guigo R, Kumpf K, Tunggal B, Cox E, Quail MA, Platzer M, Rosenthal A, Noegel AA, Barrell BG, Rajandream MA, Williams JG, Kay RR, Kuspa A, Gibbs R, Sucgang R, Muzny D, Desany B, Zeng K, Zhu B, De Jong P, et al: Sequence and analysis of chromosome 2 of Dictyostelium discoideum. Nature. 2002, 418: 79-85. 10.1038/nature00847.

    Article  PubMed  Google Scholar 

  67. ArrayExpress. [http://www.ebi.ac.uk/microarray-as/aer/]

  68. Williams KL: Characterization of dominant resistance to cobalt chloride in Dictyostelium discoideum and its use in parasexual genetic analysis. Genetics. 1978, 90: 37-47.

    PubMed  CAS  PubMed Central  Google Scholar 

  69. Williams KL: Isolation of strains of the cellular slime mold Dictyostelium discoideum capable of growing after a single passage in axenic medium. Appl Environ Microbiol. 1976, 32: 635-637.

    PubMed  CAS  PubMed Central  Google Scholar 

  70. Sussman M: Cellular differentiation in the slime mold. Growth in Living Systems. Edited by: Zarrow MX. 1960, New York: Basic Books, 221-239.

    Google Scholar 

  71. Sussman M, Sussman RR: Ploidal inheritance in Dictyostelium discoideum: stable haploid, stable diploid and metastable strains. J Gen Microbiol. 1962, 28: 417-429.

    Article  Google Scholar 

Download references

Acknowledgements

We should like to thank the Dictyostelium Stock Centre and Dennis Welker especially for strains, as well as the donors listed in Table 1. Régis Dieckmann and Navin Gopaldass for preparing DH1 and Ax2 RNA; and Nefeli Katsaridou, Kay Jagels and Theresa Feltwell for technical assistance. This work was supported by a Wellcome Trust Grant 06724, and core funding from the MRC.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Gareth Bloomfield.

Additional information

Authors' contributions

GB, JS and AI designed and built the Dictyostelium microarray. GB, AI and RRK conceived of the study. GB, YT and RRK prepared biological material. GB, JS and YT performed microarray experiments. GB and JS analyzed the data. GB and RRK prepared the paper.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Bloomfield, G., Tanaka, Y., Skelton, J. et al. Widespread duplications in the genomes of laboratory stocks of Dictyostelium discoideum. Genome Biol 9, R75 (2008). https://doi.org/10.1186/gb-2008-9-4-r75

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2008-9-4-r75

Keywords