PubMed Health⌕ Search

Biomedical subjects

M Kreitman

Publications and source records attributed to M Kreitman.

At least 19 recordsLinked to original sources

Evidence for stabilizing selection in a eukaryotic enhancer element.

Eukaryotic gene expression is mediated by compact cis-regulatory modules, or enhancers, which are bound by specific sets of transcription factors. The combinatorial interaction of these bound transcription factors determines time- and tissue-specific gene activation or repression. The even-skipped stripe 2 element controls the expression of the second transverse stripe of even-skipped messenger RNA in Drosophila melanogaster embryos, and is one of the best characterized eukaryotic enhancers. Although even-skipped stripe 2 expression is strongly conserved in Drosophila, the stripe 2 element itself has undergone considerable evolutionary change in its binding-site sequences and the spacing between them. We have investigated this apparent contradiction, and here we show that two chimaeric enhancers, constructed by swapping the 5' and 3' halves of the native stripe 2 elements of two species, no longer drive expression of a reporter gene in the wildtype pattern. Sequence differences between species have functional consequences, therefore, but they are masked by other co-evolved differences. On the basis of these results, we present a model for the evolution of eukaryotic regulatory sequences.

Animals↗

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans.

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.

Animals↗

The correlation between intron length and recombination in drosophila. Dynamic equilibrium between mutational and selective forces.

Intron length is negatively correlated with recombination in both Drosophila melanogaster and humans. This correlation is not likely to be the result of mutational processes alone: evolutionary analysis of intron length polymorphism in D. melanogaster reveals equivalent ratios of deletion to insertion in regions of high and low recombination. The polymorphism data do reveal, however, an excess of deletions relative to insertions (i.e., a deletion bias), with an overall deletion-to-insertion events ratio of 1.35. We propose two types of selection favoring longer intron lengths. First, the natural mutational bias toward deletion must be opposed by strong selection in very short introns to maintain the minimum intron length needed for the intron splicing reaction. Second, selection will favor insertions in introns that increase recombination between mutations under the influence of selection in adjacent exons. Mutations that increase recombination, even slightly, will be selectively favored because they reduce interference among selected mutations. Interference selection acting on intron length mutations must be very weak, as indicated by frequency spectrum analysis of Drosophila intron length polymorphism, making the equilibrium for intron length sensitive to changes in the recombinational environment and population size. One consequence of this sensitivity is that the advantage of longer introns is expected to decrease inversely with the rate of recombination, thus leading to a negative correlation between intron length and recombination rate. Also in accord with this model, intron length differs between closely related Drosophila species, with the longest variant present more often in D. melanogaster than in D. simulans. We suggest that the study of the proposed dynamic model, taking into account interference among selected sites, might shed light on many aspects of the comparative biology of genome sizes including the C value paradox.

Animals↗

The population genetics of the origin and divergence of the Drosophila simulans complex species.

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.

Animals↗

Dynamics of disease resistance polymorphism at the Rpm1 locus of Arabidopsis.

The co-evolutionary 'arms race' is a widely accepted model for the evolution of host-pathogen interactions. This model predicts that variation for disease resistance will be transient, and that host populations generally will be monomorphic at disease-resistance (R-gene) loci. However, plant populations show considerable polymorphism at R-gene loci involved in pathogen recognition. Here we have tested the arms-race model in Arabidopsis thaliana by analysing sequences flanking Rpm1, a gene conferring the ability to recognize Pseudomonas pathogens carrying AvrRpm1 or AvrB. We reject the arms-race hypothesis: resistance and susceptibility alleles at this locus have co-existed for millions of years. To account for the age of alleles and the relative levels of polymorphism within allelic classes, we use coalescence theory to model the long-term accumulation of nucleotide polymorphism in the context of the short-term ecological dynamics of disease resistance. This analysis supports a 'trench warfare' hypothesis, in which advances and retreats of resistance-allele frequency maintain variation for disease resistance as a dynamic polymorphism.

Arabidopsis↗

Multiple origins of cyclodiene insecticide resistance in Tribolium castaneum (Coleoptera: Tenebrionidae).

The number of origins of pesticide resistance-associated mutations is important not only to our understanding of the evolution of resistance but also in modeling its spread. Previous studies of amplified esterase genes in a highly dispersive Culex mosquito have suggested that insecticide resistance-associated mutations (specifically a single-gene duplication event) can occur a single time and then spread throughout global populations. In order to provide data for resistance-associated point mutations, which are more typical of pesticide mechanisms as a whole, we studied the number of independent origins of cyclodiene insecticide resistance in the red flour beetle Tribolium castaneum. Target-site insensitivity to cyclodienes is conferred by single point mutations in the gene Resistance to dieldrin (Rdl), which codes for a subunit of a gamma-aminobutyric acid (GABA) receptor. These point mutations are associated with replacements of alanine 302 which render the receptor insensitive to block by the insecticide. We collected 141 strains of Tribolium worldwide and screened them for resistance. Twenty-four strains contained resistant individuals. After homozygosing 23 of these resistance alleles we derived a nucleotide sequence phylogeny of the resistant strains from a 694-bp section of Rdl, encompassing exon 7 (which contains the resistance-associated mutation) and part of a flanking intron. The phylogeny also included six susceptible alleles chosen at random from a range of geographical locations. Resistance alleles fell into six clades and three clades contained both resistant and susceptible alleles. Although statistical analysis provided support at only the 5-6% level, the pattern of variation in resistance alleles is more readily explained by multiple independent origins of resistance than by spread of a single resistance-associated mutation. For example, two resistance alleles differed from two susceptible alleles only by the resistance-associated mutation itself, suggesting that they form the susceptible ancestors and that resistance arose independently in several susceptible backgrounds. This suggests that in Tribolium Rdl, de novo mutations for resistance have arisen independently in several populations. Identical alleles were found in geographically distant regions as well, also implying that some Rdl alleles have been exported in stored grain. These differences from the Culex study may stem both from differences in the population genetics of Tribolium versus that of mosquitoes and differences in mutation rates associated with point mutations versus gene duplication events. The Tribolium data therefore suggest that multiple origins of insecticide resistance (associated with specific point mutations) may be more common than the spread of single events. These findings have implications for the way in which we model the evolution and spread of insecticide resistance genes and also suggest that parallel adaptive substitutions may not be uncommon in phyletic evolution.

Animals↗

The nonrandom location of synonymous codons suggests that reading frame-independent forces have patterned codon preferences.

Biased codon usage is common in eukaryotic and prokaryotic genes. Evidence from Escherichia, Saccharomyces, and Drosophila indicates that it favors translational efficiency and accuracy. However, to date no functional advantages have been identified in the codon-anticodon interactions involving the most frequently used (preferred) codons. Here we present evidence that forces not related to the individual codon-anticodon interaction may be involved in determining which synonymous codons are preferred or avoided. We show that the "off-frame" trinucleotide motif preferences inferrable from Drosophila coding regions are often in the same direction as Drosophila's "in-frame" codon preferences, i.e., its codon usage. The off-frame preferences were inferred from the nonrandomness of the location of confamilial synonymous codons along coding regions-a pattern often described as a context dependence of nucleotide choice at synonymous positions or as codon-pair bias. We relied on randomizations of the location of confamilial codons that do not alter, and cannot be influenced by, the encoded amino acid sequences, codon usage, or base composition of the genes examined. The statistically significant congruency of in-frame and off-frame trinucleotide preferences suggests that the same kind of reading-frame-independent force(s) may also influence synonymous codon choice. These forces may have produced biases in codon usage that then led to the evolution of the translational advantages of these motifs as preferred codons. Under this scenario, tRNA pool size differences between preferred and nonpreferred codons initially were evolved to track the default overrepresentation of codons with preferred motifs. The motif preference hypothesis can explain the structuring of codon preferences and the similarities in the codon usages of distantly related organisms.

Amino Acid Sequence↗

Coding sequence evolution.

Dramatic progress has been made in the past ten years in the development of statistical and experimental techniques for investigating features of molecular evolution. Applied to coding regions, these techniques have produced remarkable advances in our understanding of selection for codon usage but, ironically, have had little impact on our understanding of protein evolution. That may be about to change.

Animals↗

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila.

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.

Animals↗

The effect of mating system differences on nucleotide diversity at the phosphoglucose isomerase locus in the plant genus Leavenworthia.

To test the theoretical prediction that highly inbreeding populations should have low neutral genetic diversity relative to closely related outcrossing populations, we sequenced portions of the cytosolic phosphoglucose isomerase (PgiC) gene in the plant genus Leavenworthia, which includes both self-incompatible and inbreeding taxa. On the basis of sequences of intron 12 of this gene, the expected low diversity was seen in both populations of the selfers Leavenworthia uniflora and L. torulosa and in three highly inbreeding populations of L. crassa, while high diversity was found in self-incompatible L. stylosa, and moderate diversity in L. crassa populations with partial or complete self-incompatibility. In L. stylosa, the nucleotide diversity was strongly structured into three haplotypic classes, differing by several insertion/deletion sequences, with linkage disequilibrium between sequences of the three types in intron 12, but not in the adjacent regions. Differences between the three kinds of haplotypes are larger than between sequences of this gene region from different species. The haplotype divergence suggests the presence of a balanced polymorphism at this locus, possibly predating the split between L. stylosa and its two inbreeding sister taxa, L. uniflora and L. torulosa. It is therefore difficult to distinguish between different potential causes of the much lower sequence diversity at this locus in inbreeding than outcrossing populations. Selective sweeps during the evolution of these populations are possible, or background selection, or merely loss of a balanced polymorphism maintained by overdominance in the populations that evolved high selfing rates.

Base Sequence↗

Unusual haplotype structure at the proximal breakpoint of In(2L)t in a natural population of Drosophila melanogaster.

The existence of temporally stable frequency clines for In(2L)t in natural populations of Drosophila melanogaster suggests a role for selection in the maintenance of this polymorphism. We have collected nucleotide polymorphism data from the proximal breakpoint junction regions of In(2L)t to infer its evolutionary history. The finding of a novel LINE-like element near the In(2L)t breakpoint junction in sampled inverted chromosomes supports a transposable element-mediated origin for this inversion. An analysis of nucleotide variation in a Costa Rican population sample of standard and inverted chromosomes indicates a unique and relatively recent origin for In(2L)t. Additional In(2L)t alleles from three geographically diverse populations reveal no detectable geographic differentiation. Low levels of In(2L)t nucleotide polymorphism suggest a recent increase in the inversion's frequency in tropical populations. An unusual feature of our sample of standard alleles is a marked heterogeneity in levels of linkage disequilibrium among polymorphic sites across the breakpoint region. We introduce a test of neutral equilibrium haplotype structure that corrects both for multiple tests and for an arbitrarily chosen window size. It reveals that an approximately 1.4-kb region immediately spanning the breakpoint has fewer haplotypes than expected under the neutral model, given the expected level of recombination in this genomic region. Certain features of our data suggest that the unusual pattern in standard chromosomes is the product of selection rather than demography.

Animals↗

Genetic variation within and among populations of Arabidopsis thaliana.

We investigated levels of nucleotide polymorphism within and among populations of the highly self-fertilizing Brassicaceous species, Arabidopsis thaliana. Four-cutter RFLP data were collected at one mitochondrial and three nuclear loci from 115 isolines representing 11 worldwide population collections, as well as from seven commonly used ecotypes. The collections include multiple populations from North America and Eurasia, as well as two pairs of collections from locally proximate sites, and thus allow a hierarchical geographic analysis of polymorphism. We found no variation at the mitochondrial locus Nad5 and very low levels of intrapopulation nucleotide diversity at Adh, Dhs1, and Gpa1. Interpopulation nucleotide diversity was also consistently low among the loci, averaging 0.0014. gst, a measure of population differentiation, was estimated to be 0.643. Interestingly, we found no association between geographical distance between populations and genetic distance. Most haplotypes have a worldwide distribution, suggesting a recent expansion of the species or long-distance gene flow. The low level of polymorphism found in this study is consistent with theoretical models of neutral mutations and background selection in highly self-fertilizing species.

Arabidopsis↗

The correlation between synonymous and nonsynonymous substitutions in Drosophila: mutation, selection or relaxed constraints?

Codon usage bias, the preferential use of particular codons within each codon family, is characteristic of synonymous base composition in many species, including Drosophila, yeast, and many bacteria. Preferential usage of particular codons in these species is maintained by natural selection acting largely at the level of translation. In Drosophila, as in bacteria, the rate of synonymous substitution per site is negatively correlated with the degree of codon usage bias, indicating stronger selection on codon usage in genes with high codon bias than in genes with low codon bias. Surprisingly, in these organisms, as well as in mammals, the rate of synonymous substitution is also positively correlated with the rate of nonsynonymous substitution. To investigate this correlation, we carried out a phylogenetic analysis of substitutions in 22 genes between two species of Drosophila, Drosophila pseudoobscura and D. subobscura, in codons that differ by one replacement and one synonymous change. We provide evidence for a relative excess of double substitutions in the same species lineage that cannot be explained by the simultaneous mutation of two adjacent bases. The synonymous changes in these codons also cannot be explained by a shift to a more preferred codon following a replacement substitution. We, therefore, interpret the excess of double codon substitutions within a lineage as being the result of relaxed constraints on both kinds of substitutions in particular codons.

Animals↗

Nucleotide variation in the triosephosphate isomerase (Tpi) locus of Drosophila melanogaster and Drosophila simulans.

DNA sequence variation in a 1.1-kb region including the coding portion of the Tpi locus was examined in 25 homozygous third-chromosome lines of Drosophila melanogaster, nine lines of Drosophila simulans, and one line of Drosophila yakuba. Our data show that the widespread allozyme polymorphism observed in cosmopolitan D. melanogaster is due to a glutamic acid substitution occurring in a phylogenetically conserved lysine that has been identified as part of the "hinged-lid" active site of the enzyme. This observation suggests that the replacement polymorphism may have important functional consequences. One replacement polymorphism was also observed in D. simulans, although its functional relevance is more difficult to assess, since it affects a site that is not strongly conserved. This amino acid change in D. simulans is associated with a single lineage possessing seven unique silent substitutions, which may be indicative of balancing selection or population subdivision. The absence of fixed amino acid differences between D. melanogaster and D. simulans and only a single difference with D. yakuba suggests that triose phosphate isomerase is under strong functional constraint. Silent variation is slightly higher for D. melanogaster than for D. simulans. Finally, we outline the general lack of evidence for old balanced polymorphisms at allozyme loci in D. melanogaster.

Alleles↗

Functional analysis of eve stripe 2 enhancer evolution in Drosophila: rules governing conservation and change.

Experimental investigations of eukaryotic enhancers suggest that multiple binding sites and trans-acting regulatory factors are often required for wild-type enhancer function. Genetic analysis of the stripe 2 enhancer of even-skipped (eve), an important developmental gene in Drosophila, provides support for this view. Given the importance of even-skipped expression in early Drosophila development, it might be predicted that many structural features of the stripe 2 enhancer will be evolutionarily conserved, including the DNA sequences of protein binding sites and the spacing between them. To test this hypothesis, we compared sequences of the stripe 2 enhancer between four species of Drosophila: D. melanogaster, D. yakuba, D. erecta and D. pseudoobscura. Our analysis revealed a large number of nucleotide substitutions in regulatory protein binding sites for bicoid, hunchback, Kruppel and giant, as well as a systematic change in the size of the enhancer. Some of the binding sites in D. melanogaster are either absent or modified in other species. One functionally important bicoid-binding site in D. melanogaster appears to be recently evolved. We, therefore, investigated possible functional consequences of sequence differences among these stripe 2 enhancers by P-element-mediated transformation. This analysis revealed that the eve stripe 2 enhancer from each of the four species drove reporter gene expression at the identical time and location in D. melanogaster embryos. Double staining of native eve protein and transgene mRNA in early embryos showed that the reporter gene mimicked native eve expression and, in every case, produced sharply defined stripes at the blastoderm stage that were coincident with eve stripe 2 protein. We argue that stripe 2 eve expression in Drosophila evolution can be viewed as being under constant stabilizing selection with respect to the location of the anterior and posterior borders of the stripe. We further hypothesize that the stripe 2 enhancer is functionally robust, so that its evolution may be governed by the fixation of both slightly deleterious and adaptive mutations in regulatory protein binding sites as well as in the spacing between binding sites. This view allows for a slow but continual turnover of functionally important changes in the stripe 2 enhancer.

Animals↗

The molecular clock revisited: the rate of synonymous vs. replacement change in Drosophila.

Rates of synonymous and nonsynonymous substitution were investigated for 24 genes in three Drosophila species, D. pseudoobscura, D. subobscura, and D. melanogaster. D. pseudoobscura and D. subobscura, two distantly related members of the obscura clade, differ on average by 0.29 synonymous nucleotide substitutions per site. D. melanogaster differs from the two obscura species by an average of 0.81 synonymous substitutions per site. Using a method developed by Gillespie, we investigated the variance to mean ratio, or Index of Dispersion, R, of substitutions along the three species' branches to test the fundamental prediction of the neutral theory of molecular evolution, E(R) = 1. For nonsynonymous substitutions, the average R, Ra is 1.6, which is not significantly different from the neutral theory prediction. Only 5 of the 24 genes had significantly large Ra valves, and 12 of the genes had Ra estimates of less than one. In contrast, the Index of Dispersion for synonymous substitutions was significantly large for 12 of the 24 genes, with an average of R(s) = 4.4, also statistically significant. These findings contrast with results for mammals, which showed overdispersion of nonsynonymous substitutions, but not of synonymous substitutions. Weak selection acting to maintain codon bias in Drosophila, but not in mammals, may be important in explaining the high variance in the rate of synonymous substitutions in this group of organisms.

Animals↗