SEARCH · Search PubMed
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Mechanisms responsible for the correlation between heterozygosity and recombination rate have been thought to be background selection and selective sweep. In addition to these two, diversity enhancing selection and/or gene conversion is suggested to be contributing to the correlation, by surveying some examples of highly polymorphic loci that tend to locate in high-recombination regions.
The Enhancer of split gene complex (E(spl)-C) is one of the more highly annotated gene regions in Drosophila, and the 12 genes within the complex help determine the spacing and patterning of adult bristles. Any E(spl)-C coding, transcribed, or cis-regulatory regions experiencing nonneutral evolution are strong candidates to harbor polymorphisms contributing to naturally occurring variation in bristle number. We confirm that the E(spl)-C is strongly conserved and show that 74% of regulatory elements previously identified in D. melanogaster are conserved in D. pseudoobscura. Regulatory elements in enhancer regions show lower nucleotide diversity and more rare polymorphisms compared with adjacent nonregulatory DNA, suggesting they are under purifying selection, and these effects are particularly pronounced when considering only conserved regulatory elements. The ratio of polymorphism to divergence was significantly different between binding sites and nonbinding sites for transcription factors within enhancer regions, suggesting the action of some form of selection. Too few polymorphisms in regions of the 3' UTR harboring regulatory motifs prevents adequate comparison of diversity and the polymorphism frequency spectrum between 3' UTR motif and nonmotif sequence. We identified at least two broad regions of the gene complex showing strong population subdivision among four populations, which is suggestive of local adaptation or background selection. Finally, two regions of the E(spl)-C exhibit low nucleotide diversity, a high level of rare polymorphisms, and an increase in linkage disequilibrium, which together suggest the action of positive selection. Notably, the gene m2 shows a significant deviation from neutrality by the McDonald-Kreitman test and resides in one of the two regions putatively experiencing a selective sweep. All sites in regions apparently visible to various selective forces are candidates for future work to determine their phenotypic effects.
Crop species experienced strong selective pressure directed at genes controlling traits of agronomic importance during their domestication and subsequent episodes of selective breeding. Consequently, these genes are expected to exhibit the signature of selection. We screened 501 maize genes for the signature of selection using microsatellites or simple sequence repeats (SSRs). We applied the Ewens-Watterson test, which can reveal deviations from a neutral-equilibrium model, as well as two nonequilibrium tests that incorporate the domestication bottleneck. We investigated two classes of SSRs: those known to be polymorphic in maize (Class I) and those previously classified as monomorphic in maize (Class II). Fifteen SSRs exhibited some evidence for selection in maize and 10 showed evidence under stringent criteria. The genes containing nonneutral SSRs are candidates for agronomically important genes. Because demographic factors can bias our tests, further independent tests of these candidates are necessary. We applied such an additional test to one candidate, which encodes a MADS box transcriptional regulator, and confirmed that this gene experienced a selective sweep during maize domestication. Genomic scans for the signature of selection offer a means of identifying new genes of agronomic importance even when gene function and the phenotype of interest are unknown.
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.
Two missense polymorphisms (E272K and L374F) of the AIM1 locus, encoding a melanocyte differentiation antigen, were shown to have a clear association with human ethnicities. These two nonpathogenic single nucleotide polymorphisms (SNPs) may be associated with human pigmentation variation. In this study, we investigated sequence variation in the coding region and exon-flanking sequence and found low genetic variation only in subjects of European descent. All four statistical tests applied to the 7.55-kb region surrounding the L374F polymorphism detected statistically significant deviations from selective neutrality in Europeans. In addition, haplotype analysis revealed that one haplotype carrying 374F was overrepresented in this population, and the low rate of variation, with some features of selective sweeps, was shown to be statistically significant. These results suggest that positive selection recently has been acting or has acted on at least this region of the melanogenic gene and that an advantageous haplotype spread rapidly in Europe.
Domesticated maize and its wild ancestor (teosinte) differ strikingly in morphology and afford an opportunity to examine the connection between strong selection and diversity in a major crop species. The tb1 gene largely controls the increase in apical dominance in maize relative to teosinte, and a region of the tb1 locus 5' to the transcript sequence was a target of selection during maize domestication. To better characterize the impact of selection at a major "domestication" locus, we have sequenced the upstream tb1 genomic region and systematically sampled nucleotide diversity for sites located as far as 163 kb upstream to tb1. Our analyses define a selective sweep of approximately 60-90 kb 5' to the tb1 transcribed sequence. The selected region harbors a mixture of unique sequences and large repetitive elements, but it contains no predicted genes. Diversity at the nearest 5' gene to tb1 is typical of that for neutral maize loci, indicating that selection at tb1 has had a minimal impact on the surrounding chromosomal region. Our data also show low intergenic linkage disequilibrium in the region and suggest that selection has had a minor role in shaping the pattern of linkage disequilibrium that is observed. Finally, our data raise the possibility that maize-like tb1 haplotypes are present in extant teosinte populations, and our findings also suggest a model of tb1 gene regulation that differs from traditional views of how plant gene expression is controlled.
Two patterns of plant defense gene evolution are emerging from molecular population genetic surveys. One is that specialist defenses experience stronger selection than generalist defenses. The second is that specialist defenses are more likely to be subject to balancing selection, i.e., evolve in a manner consistent with balanced-polymorphism or trench-warfare models of host-parasite coevolution. Because most of the data of specialist defenses come from Arabidopsis thaliana, we examined the genetic diversity and evolutionary history of three defense genes in two outcrossing species, the autotetraploid Zea perennis and its most closely related extant relative the diploid Z. diploperennis. Intraspecific diversity at two generalist defenses, the protease inhibitors wip1 and mpi, were consistent with a neutral model. Like previously studied genes in these taxa, wip1 and mpi harbored similar levels of diversity in Z. diploperennis and Z. perennis. In contrast, the specialist defense hm2 showed strong although distinctly different departures from a neutral model in the two species. Z. diploperennis appears to have experienced a strong and recent selective sweep. Using a rejection-sampling coalescent method, we estimate the strength of selection on Z. diploperennis hm2 to be approximately 3.0%, which is approximately equal to the strength of selection on tb1 during maize domestication. Z. perennis hm2 harbors three highly diverged alleles, two of which are found at high frequency. The distinctly different patterns of diversity may be due to differences in the phase of host-parasite coevolutionary cycles, although higher hm2 diversity in Z. perennis may also reflect reduced efficacy of selection in the autotetraploid relative to its diploid relative.
Patterns of DNA sequence diversity vary widely among genes encoding proteins that protect plants against pathogens and herbivores. Comparative studies may help determine whether these differences are due to the strength of selection acting on different types of defense, in different evolutionary lineages, or both. I analyzed sequence diversity at three chitinases, a well-studied component of defense, in two species of Zea and several Poaceae taxa. Although the Zea species are closely related and these genes code for proteins with similar biochemical function, patterns of diversity vary widely within and among species. Intraspecific diversity at chiB, chiI, and Z. mays ssp. parviglumis chiA are consistent with a neutral-equilibrium model whereas chiA had no segregating sites within Z. diploperennis--consistent with a recent and strong selective sweep. Codons identified as having diverged among Poaceae taxa in response to positive selection were significantly overrepresented among targets of selection in Arabis, suggesting common responses to selection in distantly related plant taxa. Divergence of the recent duplicates chiA and chiB is consistent with positive selection but relaxed constraint cannot be rejected. Weak evidence for adaptive divergence of these duplicated downstream components of defense contrasts with strong evidence for adaptive divergence of genes involved in pathogen recognition.
Genetic differentiation between the largely sympatric molecular forms M and S of Anopheles gambiae appears mostly limited to division 6 and part of division 5 of the X chromosome. This region is adjacent to the centromere and includes the rDNA that was used to define these forms. This localized differentiation between populations that experience gene flow strongly suggests that this region contains genes responsible for reproductive isolation. Regions adjacent to centromeres are known to experience less recombination in several species and it has recently been suggested that low recombination rates can facilitate the accumulation and maintenance of isolation genes in partially isolated populations. Therefore, we measured the recombination rate in division 5D/6 directly and estimate that it is at least 16-fold reduced across this region compared to the remainder of the X chromosome. Additionally, sequence data from four loci from field-collected mosquitoes from several West African countries show very strong differentiation between the molecular forms in division 5D/6, whereas none was observed in two loci elsewhere on the X chromosome. Furthermore, genetic variation was substantially lower in division 5D/6 compared to the two reference loci, and the inferred genealogies of the division 5D/6 genes show patterns consistent with selective sweeps. This suggests that the reduced recombination rate has increased the effect of selection on this region and that our data are consistent with the hypothesis that reduced recombination rates can play a role in the accumulation of isolation genes in the face of gene flow.
Mustached bats emit echolocation and communication calls containing both constant frequency (CF) and frequency-modulated (FM) components. Previously we found that 86% of neurons in the ventral division of the external nucleus of the inferior colliculus (ICXv) were directionally selective for linear FM sweeps and that selectivity was dependent on sweep rate. The ICXv projects to the suprageniculate nucleus (Sg) of the medial geniculate body. In this study, we isolated 37 single units in the Sg and measured their responses to best excitatory frequency (BEF) tones and linear 12-kHz upward and downward FM sweeps centered on the BEF. Sweeps were presented at durations of 30, 12, and 4 ms, yielding modulation rates of 400, 1,000, and 3,000 kHz/s. Spike count versus level functions were obtained at each modulation rate and compared with BEF controls. Sg units responded well to both tones and FM sweeps. BEFs clustered at 58 kHz, corresponding to the dominant CF component of the sonar signal. Spike count functions for both tones and sweeps were predominantly non-monotonic. FM directional selectivity was significant in 53-78% of the units, depending on modulation rate and level. Units were classified as up-selective (52%), down-selective (24%), or bi-directional (non-selective, 16%); a few units (8%) showed preferences that were either rate- or level-dependent. Most units showed consistent directional preferences at all SPLs and modulation rates tested, but typically showed stronger selectivity at lower sweep rates. Directional preferences were attributable to suppression of activity by sweeps in the non-preferred direction (~80% of units) and/or facilitation by sweeps in the preferred direction (~20-30%). Latencies for BEF tones ranged from 4.9 to 25.7 ms. Latencies for FM sweeps typically varied linearly with sweep duration. Most FM latency-duration functions had slopes ranging from 0.4 to 0.6, suggesting that the responses were triggered by the BEF. Latencies for BEF tones and FM sweeps were significantly correlated in most Sg units, i.e., the response to FM was temporally related to the occurrence of the BEF in the FM sweep. FM latency declined relative to BEF latency as modulation rate increased, suggesting that at higher rates response is triggered by frequencies in the sweep preceding the BEF. We conclude that Sg and ICXv units have similar, though not identical, response properties. Sg units are predominantly upsweep selective and could respond to either or both the CF and FM components in biosonar signals in a number of echolocation scenarios, as well as to a variety of communication sounds.
Sex chromosomes may provide a context for studying the local effects of mutation rate on molecular evolution, since the two types of sex chromosomes are generally exposed to different mutational environments in male and female germ lines. Importantly, recent studies of some vertebrates have provided evidence for a higher mutation rate among males than among females. Thus, in birds, the Z chromosome, which spends two thirds of its time in the male germ line, is exposed to more mutations than the female-specific W chromosome. We show here that levels of nucleotide diversity are drastically higher on the avian Z chromosome than in paralogous sequences on the W chromosome. In fact, no intraspecific polymorphism whatsoever was seen in about 3.4 kb of CHD1W intron sequence from a total of >150 W chromosome copies of seven different bird species. In contrast, the amount of genetic variability in paralogous sequences on the Z chromosome was significant, with an average pairwise nucleotide diversity (d) of 0.0020 between CHD1Z introns and with 37 segregating sites in a total of 3.8 kb of Z sequence. The contrasting levels of genetic variability on the avian sex chromosomes are thus in a direction predicted from a male-biased mutation rate. However, although a low gene number, as well as some other factors, argues against background selection and/or selective sweeps shaping the genetic variability of the avian W chromosome, we cannot completely exclude selection as a contributor to the low levels of variation on the W chromosome.
Neutral mutations may hitchhike to high frequency when they are situated close to sites under positive selection, generating local reductions in genetic diversity. This process is thought to be an important determinant of levels of genomic variation in natural populations. The size of genome regions affected by genetic hitchhiking is expected to be dependent on the strength of selection, but there is little empirical data supporting this prediction. Here, we compare microsatellite variation around two drug resistance genes (chloroquine resistance transporter (pfcrt), chromosome 7, and dihydrofolate reductase (dhfr), chromosome 4) in malaria parasite populations exposed to strong (Thailand) or weak selection (Laos) by anti-malarial drugs. In each population, we examined the point mutations underlying resistance and length variation at 22 (chromosome 4) or 25 (chromosome 7) microsatellite markers across these chromosomes. All parasites from Thailand carried the K76T mutation in pfcrt conferring resistance to chloroquine (CQ) and 2-4 mutations in dhfr conferring resistance to pyrimethamine. By contrast, we found both wild-type and resistant alleles at both genes in Laos. There were dramatic differences in the extent of hitchhiking in the two countries. The size of genome regions affected was smaller in Laos than in Thailand. We observed significant reduction in variation relative to sensitive parasites for 34-64 kb (2-4 cM) in Laos on chromosome 4, compared with 98-137 kb (6-8 cM) in Thailand. Similarly, on chromosome 7, we observed reduced variation for 34-69 kb (2-4 cM) around pfcrt in Laos, but for 195-268 kb (11-16 cM) in Thailand. Reduction in genetic variation was also less extreme in Laos than in Thailand. Most loci were monomorphic in a 12 kb region surrounding both genes on resistant chromosomes from Thailand, whereas in Laos, even loci immediately proximal to selective sites showed some variation on resistant chromosomes. Finally, linkage disequilibrium (LD) decayed more rapidly around resistant pfcrt and dhfr alleles from Laos than from Thailand. These results demonstrate that different realizations of the same selective sweeps may vary considerably in size and shape, in a manner broadly consistent with selection history. From a practical perspective, genomic regions containing resistance genes may be most effectively located by genome-wide association in populations exposed to strong drug selection. However, the lower levels of LD surrounding resistance alleles in populations under weak selection may simplify identification of functional mutations.
In finite populations, linkage disequilibria generated by the interaction of drift and directional selection (Hill-Robertson effect) can select for sex and recombination, even in the absence of epistasis. Previous models of this process predict very little advantage to recombination in large panmictic populations. In this article we demonstrate that substantial levels of linkage disequilibria can accumulate by drift in the presence of selection in populations of any size, provided that the population is subdivided. We quantify (i) the linkage disequilibrium produced by the interaction of drift and selection during the selective sweep of beneficial alleles at two loci in a subdivided population and (ii) the selection for recombination generated by these disequilibria. We show that, in a population subdivided into n demes of large size N, both the disequilibrium and the selection for recombination are equivalent to that expected in a single population of a size intermediate between the size of each deme (N) and the total size (nN), depending on the rate of migration among demes, m. We also show by simulations that, with small demes, the selection for recombination is stronger than both that expected in an unstructured population (m = 1 - 1/n) and that expected in a set of isolated demes (m = 0). Indeed, migration maintains polymorphisms that would otherwise be lost rapidly from small demes, while population structure maintains enough local stochasticity to generate linkage disequilibria. These effects are also strong enough to overcome the twofold cost of sex under strong selection when sex is initially rare. Overall, our results show that the stochastic theories of the evolution of sex apply to a much broader range of conditions than previously expected.
Natural selection, in the form of balancing selection or selective sweeps, can result in a decoupling of the amounts of molecular polymorphism and divergence. Thus natural selection can cause some areas of DNA sequence to have greater silent polymorphism, relative to divergence between species, than other areas. It would be useful to have a statistical test for heterogeneity in the polymorphism to divergence ratio across a region of DNA sequence, one that could identify heterogeneity greater than that expected from the neutral processes of mutation, drift, and recombination. The only currently available test requires that a region be arbitrarily divided into sections that are compared with each other, and the subjectivity of this division could be problematic. Here a test is proposed in which runs of polymorphic and fixed sites are counted, where a "run" is a set of one or more sites of one type preceded and followed by the other type. The number of runs is smaller than otherwise expected if polymorphisms are clumped together. By simulating neutral evolution and comparing the observed number of runs to the simulations, a statistical test is possible which does not require any a priori decisions about subdivision.
The maintenance of mobile DNA sequences in clonal organisms has been seen as a paradox. If selfish mobile sequences spread through genomes only by overreplication in transposition, then sexuality is necessary for their spread through populations. The persistence of bacterial transposable elements without obvious dominant selectable markers has previously been explained by horizontal transfer. However, advantageous insertions of mobile DNAs are known in bacteria. Here we model maintenance of an otherwise selfish mobile DNA element in a clonal species in which selection for null mutations occurs during one of two temporally alternating environments. Large areas of parameter space permit maintenance of mobile DNAs where, without selection, they would have gone extinct. Horizontal transfer diminishes, rather than enhances, mean copy number. In finite populations, effective population sizes are greatly reduced by selective sweeps, and mean copy number can be increased as the reduced variance in copy number results in reduced selection.
The developmental gene optomotor-blind (omb) encodes a T-box-containing transcription factor that has multiple roles in Drosophila development. Previous genetic analyses established that omb plays a key role in establishing the abdominal pigmentation pattern of Drosophila melanogaster. In this report we examine patterns of omb nucleotide variation in D. polymorpha, a species that is highly polymorphic for the phenotype of abdominal pigmentation. Haplotypes at this locus fall into two classes that are separated by six mutational steps; five of these mutational events result in amino acid changes. Two lines of evidence are consistent with a role for omb in the abdominal pigmentation polymorphism of D. polymorpha. First, we find that haplotype classes of omb are correlated with abdominal pigmentation phenotypes, as are microsatellite repeat numbers in the region. Second, tests of selection reveal that the two haplotype classes have been maintained by balancing selection. Within each class there is a significantly low amount of diversity, indicative of previous selective sweeps. An analysis including D. polymorpha's closest relatives (members of the cardini group) provides evidence for directional selection across species. Selection at this locus is expected if omb contributes to variation in abdominal pigmentation, since this trait is likely of ecological importance.
In geographic areas where pied and collared flycatchers (Ficedula hypoleuca and F. albicollis) breed in sympatry, hybridization occurs, leading to gene flow (introgression) between the two recently diverged species. Notably, while such introgression is observable at autosomal loci it is apparently absent at the Z chromosome, suggesting an important role for genes on the Z chromosome in creating reproductive isolation during speciation. To further understand the role of Z-linked loci in the formation of new species, we studied genetic variation of the two species from regions where they live in allopatry. We analyzed patterns of polymorphism and divergence in introns from 9 Z-linked and 23 autosomal genes in pied and collared flycatcher males. Average variation on the Z chromosome is greatly reduced compared to neutral expectations based on autosomal diversity in both species. We also observe significant heterogeneity between patterns of polymorphism and divergence at Z-linked loci and a relative absence of polymorphisms that are shared by the two species on the Z chromosome compared to the autosomes. We suggest that these observations may indicate the action of recurrent selective sweeps on the Z chromosome during the evolution of the two species, which may be caused by sexual selection acting on Z-linked genes. Alternatively, reduced variation on the Z chromosome could result from substantially higher levels of introgression at autosomal than at Z-linked loci or from a complex demographic history, such as a population bottleneck.
Malaria is a major human parasitic disease caused by four species of Plasmodium protozoa. Plasmodium vivax, the most widespread, affects millions of people across Africa, Asia, the Middle East, and Central and South America. We have studied the genetic variability of 13 microsatellite loci in 108 samples from 8 localities in Asia, Africa, South America, and New Guinea. Only one locus is polymorphic; nine are completely monomorphic, and the remaining three are monomorphic in all but one or two populations, which have a rare second allele. In contrast, Plasmodium falciparum displays extensive microsatellite polymorphism within and among populations. We further have analyzed, in 96 samples from the same 8 localities, 8 tandem repeats (TRs) located on a 100-kb contiguous chromosome segment described as highly polymorphic. Each locus exhibits 2-10 alleles in the whole sample but little intrapopulation polymorphism (1-5 alleles with a prevailing allele in most cases). Eight microsatellite loci monomorphic in P. vivax are polymorphic in three of five Plasmodium species related to P. vivax (two to seven individuals sampled). Plasmodium simium, a parasite of New World monkeys, is genetically indistinguishable from P. vivax. At 13 microsatellite loci and at 7 of the 8 TRs, both species share the same (or most common) allele. Scarce microsatellite polymorphism may reflect selective sweeps or population bottlenecks in recent evolutionary history of P. vivax; the differential variability of the TRs may reflect selective processes acting on particular regions of the genome. We infer that the world expansion of P. vivax as a human parasite occurred recently, perhaps <10,000 years ago.