Publications

The spv operon is common to all Salmonella virulence plasmids. DNA hybridization analysis indicates that the spv region is limited in distribution to serovars of Salmonella enterica subspecies I, II, IIIa, IV, and VII and is absent from Salmonella bongori isolates. Among strains of subspecies II, IIIa, and VII, all isolates examined contained sequences that hybridized with the spv region. However, among isolates of subspecies I, DNA sequences capable of hybridizing with the spv region were found in some isolates of certain serovars. Furthermore, in isolates of subspecies I, the virulence plasmid was found in the same set of isolates as an F-related plasmid, as determined by the presence of the spv region of the virulence plasmid and the finO, traD, and repA sequences of the F-plasmid. The concordance of the virulence plasmid and all three F-plasmid sequences in subspecies I serovar Choleraesuis, Paratyphi, and Typhimurium is most easily explained if the spv region is carried in an F-related plasmid in these isolates. In contrast, among S. enterica subspecies II, IIIa, IV, and VII, the isolates that contain spv sequences did not hybridize with an F-related plasmid or any other identifiable plasmid. With the use of pulse-field gel electrophoresis, the spv region in subspecies II, IIIa, and VII was found to be encoded on the chromosome. Analysis of the phylogenetic distribution of spv among Salmonella isolates and comparative nucleotide sequence analysis of spvA and spvC suggests that the spv region was acquired very recently, after speciation of the salmonellae.

The pattern of genetic variation across the genome of Drosophila melanogaster is consistent with the occurrence of frequent 'selective sweeps', in which new favourable mutations become incorporated into the species so quickly that linked alleles can 'hitchhike' and also become fixed. Because of the hitchhiking of linked genes, it is generally difficult to identify the target of any putative selective sweep. Here, however, we identify a new gene in D. melanogaster that codes for a sperm-specific axonemal dynein subunit. The gene has a new testes-specific promoter derived from a protein-coding region in a gene encoding the cell-adhesion protein annexin X (AnnX), and it contains a new protein-coding exon derived from an intron in a gene encoding a cytoplasmic dynein intermediate chain (Cdic). The new transcription unit, designated Sdic (for sperm-specific dynein intermediate chain), has been duplicated about tenfold in a tandem array. Consistent with the selective sweep of this gene, the level of genetic polymorphism near Sdic is unusually low. The discovery of this gene supports other results that point to the rapid molecular evolution of male reproductive functions.

Most theoretical models in population genetics fail to deal in a realistic manner with the process of mutation. They are consequently not informative about the central evolutionary problem of the origin, progression, and limit of adaptation. Here we present an explicit distribution of phenotypes expected in an ensemble of populations under a mutation-selection-drift model that allows mutations with a distribution of adaptive values to occur randomly in time. The model of mutation is a geometrical model in which the effect of a new mutation is determined by a random angle in n dimensional space and in which the adaptive value (fitness) of an organism decreases as the square of the deviation of its phenotype from an optimum. Each new mutation is subjected to random genetic drift and fixed or lost according to its selective value and the effective population number. Time is measured in number of fixation events, so that, at any point in time, each population is regarded as genetically homogeneous. In this mutation-selection-drift model, among an ensemble of populations, the equilibrium average phenotype coincides with the optimum because the distribution of positive and negative deviations from the optimum is symmetrical. However, at equilibrium the mean of the absolute value of the deviation from the optimum equals square root of n-/8Ns), where n is the dimensionality of the trait space, N is the effective population size, and s is the selection coefficient against a mutation whose phenotype deviates by one unit from the optimum. Furthermore, at equilibrium, the average fitness across the ensemble of populations equals 1 - (n + 1)/8N. When n is sufficiently large, there is a strong mutation pressure toward the fixation of slightly deleterious mutations. This feature relates our model to the nearly neutral theory of molecular evolution.

With the increased popularity of zebrafish (Danio rerio) for mutagenesis studies, efficient methods for manipulation of its genome are needed. One approach is the use of a transposable element as a vector for gene transfer in this species. We report here the transformation of zebrafish and germ-line transmission of the mariner element from Drosophila mauritiana. The mariner element was selected because its transposition is independent of host-specific factors. One- to two-cell-stage zebrafish embryos were coinjected with a supercoiled plasmid carrying the nonautonomous mariner element peach and mRNA encoding the transposase. Surviving larvae were reared to adulthood, and the transmission of peach to the F1 generation was tested by PCR. Four of the 12 founders, following plasmid injections on 2 different days, transmitted the element to their progeny. Inheritance of the transgene from the F1 to the F2 generation showed a Mendelian pattern. No plasmid sequences were detected by PCR or Southern blot analysis, indicating transposition of peach rather than random integration of the plasmid DNA. These data provide evidence of transformation of a vertebrate with a transposable element and support the host-independent mechanism for transposition of the mariner element. We suggest this system could be used for insertional mutagenesis or for identifying active regions of the genome in the zebrafish.

In Drosophila, the availability of polytene chromosome maps and of sets of probes covering most regions of the chromosomes allows a direct comparison of the organization of the genome in different species. In this work, we report the localization, in Drosophila virilis, D. montana, and D. novamexicana, of > 100 bacteriophage PI clones containing approximately 65 kilobase inserts of genomic DNA from D. virilis. Each clone hybridizes with a single euchromatic site in either chromosome 1 or chromosome 3 in D. virilis. From these data, it is possible to estimate the minimum number of inversions required to transform the map positions of the probes in one species into the map positions of the same probes in a related species. The data indicate that, in the D. virilis species group, the X chromosome has up to four times the number of inversions as are observed in chromosome 3. The first photographic polytene chromosome maps for D. montana and D. novamexicana are also presented.

We analyzed a 5,770-bp genomic region of Drosophila virilis that contains a cluster of two maltase genes showing sequence similarity with genes in a cluster of three maltase genes previously identified in Drosophila melanogaster. The D. virilis maltase genes are designated Mav1 and Mav2. In addition to being different in gene number, the cluster of genes in D. virilis differs dramatically in intron-exon structure from the maltase genes in D. melanogaster, the transcriptional orientation of the genes in the cluster also differs between the species. Our findings support a model in which the maltase gene cluster in D. virilis and D. melanogaster evolved independently. Furthermore, while in D. melanogaster the maltase gene cluster lies only 10 kb distant from the larval cuticle gene cluster, the maltase and larval cuticle gene clusters in D. virilis are located very far apart and on a different chromosome than that expected from the known chromosome arm homologies between D. virilis and D. melanogaster. A region of the genome containing the maltase and larval cuticle gene clusters appears to have been relocated between nonhomologous chromosomes.

The analysis of patterns of genome evolution may help to evaluate the evolutionary forces that shape the composition and organization of the genome. Comparisons between the physical maps of divergent species can be used to identify conserved blocks of closely linked genes whose synteny is possibly under selective constraint. We have used in situ hybridization to determine the genomic position of 732 randomly selected clones from a bacteriophage P1 library of Drosophila virilis. The resulting map includes at least one clone in each of 69% of the subdivisions into which the D. virilis polytene chromosomes are divided. A subset of these clones was used to carry out a comparative physical analysis of chromosome 2 from D. virilis and from Drosophila montana. A number of discrepancies with the classical scenario of chromosome evolution were noted. The D. virilis P1 clones were also used to determine the physical relations between ten genes that are located in the X chromosome of Drosophila melanogaster between the markers crn (2F1) and omb (4C5-6). In this region, which is approximately 2 Mb in length, there have been at least six breakpoints since the divergence of the species, and six of the genes are found at widely scattered locations in the D. virilis X chromosome. However, a block of four functionally unrelated genes, including white, roughest, Notch, and dunce, seems to be conserved between the two species.

To determine whether nuclear rDNA sequences provide a useful means for assessing the structure of populations of Ixodes ticks, we compared variability among copies of an internal transcribed spacer (ITS-2) sequence within individual ticks to the variability between ticks. At least 4% of the nucleotides comprising this sequence vary among the copies present within individual ticks. ITS-2 diversity in each of two ticks is nearly half as great as that reported between ticks from geographically disparate populations. Because individual ticks retain ancestral polymorphism, ITS-2 variation does not accurately reflect descent relationships among these ticks. Sequencing single copies of PCR-amplified ITS-2 therefore does not permit assessment of the phylogenetic relationships among the I. ricinus-like ticks in eastern North America. We recommend caution in future analyses, and emphasize the importance of procedures designed to ensure that the many paralogous copies of the rDNA cistron have been sufficiently homogenized by concerted evolutionary processes. Such precautionary measures will make certain that phylogenetic trees based on these gene sequences reflect the phyletic relatedness of the biological species.

The mariner/Tc1 superfamily of transposable elements is one of the most diverse and widespread Class II transposable elements. Within the larger assemblage, the mariner-like elements (MLEs) and the Tc1-like elements (TLEs) are distinct families differing characteristically in the composition of the "D,D(35)E" cation-binding domain. Based on levels of sequence similarity, the elements in each family can be subdivided further into several smaller subfamilies. MLEs and TLEs both have an extraordinarily wide host range. They are abundant in insect genomes and other invertebrates and are found even in some vertebrate species including, in the case of mariner, humans, in which one element on chromosome 17p has been implicated as a hotspot of recombination. In spite of the extraordinary evolutionary success of the elements, virtually nothing is known about their mode of regulation within genomes. There is abundant evidence that the elements are disseminated to naive host genomes by horizontal transmission, and there is a substantial base of evidence for inference about the subsequent population dynamics. Studies of engineered mariner elements and induced mutations in the transposase have identified two mechanisms that may be operative in mariner regulation. One mechanism is overproduction inhibition, in which excessive wild-type transposase reduces the rate of excision of a target element. A second mechanism is dominant-negative complementation, in which certain mutant transposase proteins antagonize the activity of the wild-type transposase. The latter process may help explain why the vast majority of MLEs in nature undergo "vertical inactivation" by multiple mutations and, eventually, stochastic loss. There is also evidence that mariner/Tc1 elements can be mobilized in hybrid dysgenesis; in particular, certain dysgenic crosses in Drosophila virilis result in mobilization of a TLE designated Paris as well as the mobilization of other unrelated transposable elements.

Genetic analysis of eukaryote transposases and comparison with their prokaryote counterparts have been greatly hindered by difficulty in isolating mutations. We describe a simple eye-color screen that facilitates isolation and analysis of mutations in the mariner transposase in Drosophila melanogaster. Use of ethyl methanesulfonate and site-directed mutagenesis has identified 18 residues that are critical for in vivo excision of a target mariner element. When the mutations were examined in heterozygous mutant/nonmutant genotypes, more than half of the mutant transposase proteins were found to reduce the activity of the wild-type transposase, as assayed by the frequency of germline excision of a target element. Remarkably, transposase function is obliterated when the D,D(34)D acidic, ion-binding domain is replaced with the consensus sequence D,D(34)E found in the nematode Tc1 transposase and in many other transposases in the superfamily. A number of mutations strongly complement wild-type transposase in a dominant-negative manner, suggestive of subunit interactions in the excision reaction; these mutations are located in a small region that includes part of the D,D(34)D motif. Transposase function also is eliminated by a mutation in the inferred initiation codon and by a mutation in a putative nuclear localization signal.

We have studied the spatial distribution of IS1 elements in the genomes of natural isolates comprising the ECOR reference collection of Escherichia coli. We find evidence for nonrandomness at three levels. Many pairs of IS1 elements are in much closer proximity (< 10 kb) than can be accounted for by chance. IS1 elements in close proximity were identified by long-range PCR amplification of the genomic sequence between them. Each amplified region was sequenced and its map location determined by database screening of DNA hybridization. Among the ECOR strains with at least two IS1 elements, 54% had one or more pairs of elements separated by < 10 kb. We propose that this type of clustering is a result of "local hopping," in which we assume that a significant proportion of tranposition events leads to the insertion of a daughter IS element in the vicinity of the parental element. A second level of nonrandomness is found in strains with a modest number of IS1 elements that are mapped through the use of inverse PCR to amplify flanking genomic sequences: in these strains, the insertion sites tend to be clustered over a smaller region of chromosome than would be expected by chance. A third level of nonrandomness is observed in the composite distribution of IS elements across strains: among 20 mapped IS1 elements, none were found in the region of 48-77 minutes, a significant gap. One region of the E. coli chromosome, at 98 min, had a cluster of IS1 elements in seven ECOR strains of diverse phylogenetic origin. We deduce from sequence analysis that this pattern of distribution is a result of initial insertion in the most recent common ancestor of these strains and therefore not a hot spot of insertion. Analysis using long-range PCR with primers for IS2 and IS3 also yielded pairs of elements in close proximity, suggesting that these elements may also occasionally transpose by local hopping.

Seventy-one natural isolates obtained from a Salmonella reference collection were examined for the presence of plasmids closely related to the Escherichia coli F plasmid. The collection consists of several serovars of the S. enterica Typhimurium complex, subspecies I, to which 99% of pathogenic salmonellae belong. Molecular genetic techniques of DNA hybridization, along with PCR and DNA sequencing, were used to examine the occurrence, distribution, and genetic diversity of F-like plasmids among Salmonella strains. The F plasmid genes examined were finO, traD, traY, and repA, which map at dispersed positions on the F plasmid of E. coli. Comparative sequence analysis of each of the four genes in Salmonella plasmids showed them to be homologous (in some cases, virtually identical) to those found in F plasmids of E. coli natural isolates. Furthermore, the frequency of F-like plasmids in Salmonella strains was approximately the same as that observed in the E. coli Reference Collection. However, in Salmonella, the distribution was confined predominately to the serovars Typhimurium and Muenchen. The unexpected finding of a shared pool of F-like plasmids between S. enterica and E. coli demonstrates the significant role of conjugation in the histories of these important bacterial species.

The mariner/Tcl superfamily of transposable elements is widely distributed in animal genomes and is especially prevalent in insects. Their wide distribution results from their ability to be disseminated among hosts by horizontal transmission and also by their ability to persist in genomes through multiple speciation events. Although a great deal is known about the molecular mechanisms of transposition and excision, very little is known about the mechanisms by which transposition is controlled within genomes. The issue of mariner/Tcl regulation is critical in view of the great interest in these elements as vectors for germline transformation of insect pests and vectors of human disease. Several potentially important regulatory mechanisms have been identified in studies of genetically engineered mariner elements. One mechanism is overproduction inhibition, in which excessive wild-type transposase reduces the rate of excision of a target element. A second mechanism is mediated by certain mutant transposase proteins, which antagonize the activity of the wild-type transposase. The latter process may help explain why the vast majority of MLEs in nature undergo 'vertical inactivation' by multiple mutations and, eventually, stochastic loss. Another potential mechanism of regulation may result from transposase titration by defective elements that retain their DNA binding sites and ability to transpose. There is also evidence that some mariner/Tcl elements can be mobilized in a type of hybrid dysgenesis.

We have recently described a novel method of estimating neutral rates and patterns of spontaneous mutation (Petrov et al., 1996). This method takes advantage of the propensity of non-LTR retrotransposable elements to create non-functional, 'dead-on-arrival' copies as a product of transposition. Maximum parsimony analysis is used to separate the evolution of actively transposing lineages of a non-LTR element from the fate of individual inactive insertions, and thereby allows one to assess directly the relative rates of different types of mutation, including point substitutions, deletions and insertions. Because non-LTR elements enjoy wide phylogenetic distribution, this method can be used in taxa that do not harbor a significant number of bona fide pseudogenes, as is the case in Drosophila (Jeffs and Ashburner, 1991; Weiner et al., 1986). We used this method with Helena, a non-LTR retrotransposable element present in the Drosophila virilis species group. A striking finding was the virtual absence of insertions and remarkably high incidence of large deletions, which combine to produce a high overall rate of DNA loss. On average, the rate of DNA loss in D. virilis is approximately 75 times faster than that estimated for mammalian pseudogenes (Petrov et al., 1996). The high rate of DNA loss should lead to rapid elimination of non-essential DNA and thus may explain the seemingly paradoxical dearth of pseudogenes in Drosophila. Varying rates of DNA loss may also contribute to differences in genome size (Graur et al., 1989; Petrov et al., 1996), thus explaining the celebrated 'C-value' paradox (John and Miklos, 1988). In this paper we outline the theoretical basis of our method, examine the data from this perspective, and discuss potential problems that may bias our estimates.

A number of mechanisms have recently been described that might be important in restricting the level of activity of mariner-like transposable elements (MLEs) in natural populations. These mechanisms include overproduction inhibition, in which increasing the dose of transposase decreases net activity. Another mechanism is mediated by certain missense mutations, in which a mutant transposase protein impairs the activity of the wild-type transposase in heterozygous mutant/nonmutant genotypes. A further mechanism is the potential for transposase titration by defective elements that retain transposase binding activity. The issue of regulation is not only of theoretical importance in understanding the molecular and evolutionary genetics of MLEs, but also of practical significance in learning how best to use MLEs in the germline transformation of insect pests and disease vectors.

Genetic studies of the mariner transposable element Mos1 have revealed two novel types of regulatory mechanisms. In one mechanism, overproduction of the wild-type transposase reduces the overall level of transposase activity as assayed by the excision of a nonautonomous mariner target element. This mechanism is termed overproduction inhibition (OPI). Another mechanism is observed in a class of hypomorphic missense mutations in the transposase. In the presence of wild-type Mos1 transposase, these mutations exhibit dominant-negative complementation (DNC) that antagonizes the activity of the wild-type transposase. We propose that these regulatory mechanisms act at the level of the transposase protein subunits by promoting the assembly of oligomeric forms, or of mixed-subunit oligomers, that have reduced activity. We suggest that these regulatory mechanisms may apply generally to mariner-like elements (MLEs). Overproduction inhibition may help explain why the MLE copy number reaches very different levels in different species. Dominant-negative complementation may help explain why most naturally occurring copies of MLEs have been mutationally inactivated.

One implication of Kacser's analysis of complex metabolic systems is that mutations with small effects exist as a consequence of the typically small flux control coefficient relating enzyme activity to the rate of a metabolic process. Although a slightly detrimental mutation is somewhat less likely to become fixed by chance than a slightly favorable mutation, mutations that are slightly detrimental might be expected to be more numerous than favorable mutations owing to the previous incorporation of favorable mutations by a long history of natural selection. The result is that, as Ohta has pointed out, a significant fraction of mutations that are fixed in evolution are slightly detrimental. In the long run, the fixation of detrimental mutations in a gene increases the opportunity for the occurrence of a compensatory favorable mutation, either in the same gene or in an interacting gene. On a suitably long timescale, therefore, every gene incorporates favorable mutations that compensate for detrimental mutations previously fixed. This form of evolution is driven primarily by natural selection, but it results in no change or permanent improvement in enzymatic function.

A recently published study has identified a set of candidate genes for human diseases based on findings from Drosophila. Each human expressed sequence tag (EST) in a large database was compared with all known Drosophila genes. After eliminating matches between genes of already known function, the remaining sequences were mapped in the human genome. In each region, the phenotypes of all known human diseases were compared with the phenotypes of known Drosophila mutations in order to identify candidate genes for the human diseases. Are the correspondences real or coincidental?

A laboratory strain of Drosophila virilis was genetically transformed with a hobo vector carrying the miniwhite cassette using a helper plasmid with an hsp70-driven hobo transposase-coding sequence. The rate of transformation was 0.5% per fertile GO animal. Three transgenic insertions were cloned and characterized and found to be authentic hobo insertions. These results, together with the known widespread distribution of hobo in diverse insect species, suggest that hobo and related transposable elements may be of considerable utility in the germline transformation of insects other than D. melanogaster.