Estimates of linkage disequilibrium and the recombination parameter determined from segregating nucleotide sites in the alcohol dehydrogenase region of Drosophila pseudoobscura

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

Divergence and Gene Flow Between Fusarium subglutinans and F. temperatum Isolated from Maize in Argentina

Fusarium subglutinans and F. temperatum are two important fungal pathogens of maize whose distinctness as separate species has been difficult to assess. We isolated strains of these species from commercial and native maize varieties in Argentina and sequenced >28,000 loci to estimate genetic variation in the sample. Our objectives were to measure genetic divergence between the species, infer demographic parameters related to their split, and describe the population structure of the sample. When analyzed together, over 30% of each species' polymorphic sites (>2,500 sites) segregate as polymorphisms in the other. Demographic modeling confirmed the species split predated maize domestication, but subsequent between-species gene flow has occurred, with gene flow from F. subglutinans into F. temperatum greater than gene flow in the reverse direction. In F. subglutinans, little evidence exists for substructure or recent selective sweeps, but there is evidence for limited sexual reproduction. In F. temperatum, there is clear evidence for population substructure and signals of abundant recent selective sweeps, with sexual reproduction probably less common than in F. subglutinans. Both genetic variation and the relative number of polymorphisms shared between species increase near the telomeres of all 12 chromosomes, where genes related to plant-pathogen interactions often are located. Our results suggest that species boundaries between closely related Fusarium species can be semipermeable and merit further study. Such semipermeability could facilitate unanticipated genetic exchange between species and enable quicker permanent responses to changes in the agro-ecosystem, e.g., pathogen-resistant host varieties, new chemical and biological control agents, and agronomic practices.

Genome-wide analyses of a plant-specific LIM-domain gene family implicate its evolutionary role in plant diversification

The Arabidopsis DA1 genes appear to have multiple functions in regulating organ size and abiotic stress response, but the biological roles of its closely related genes remain unknown. Evolutionary analyses might provide some clues to aid in an understanding of their functional diversification. In this work, we characterized the molecular evolution and expressional diversification of DA1-like genes. Surveying 354 sequenced genomes revealed 142 DA1-like genes only in plants, indicating plant-specificity of these genes. The DA1-like protein modular structure was composed of two UIMs (ubiquitin interaction motifs), one LIM-domain (from lin-11, isl-1, and mec-3), and a conserved C-terminal, and was distinguishable from the already defined three groups of LIM-domain proteins. We further found that the DA1-like genes diverged into Classes I and II at the ancestor of seed plants and acquired 13 clade-specific residues during their evolutionary history. Moreover, diverse intron size evolution was noted following the transition from size-expandable introns to minimal ones, accompanying the emergence and diversification of angiosperms. Functional diversification as it relates to gene expression was further investigated in soybean. Glycine max DA1 genes showed diverse tissues expression patterns during development and had substantially varied abiotic stress response expression. Thus, variations in the coding regions, intron size, and gene expression contributed to the functional diversification of this gene family. Our data suggest that the evolution of the DA1-like genes facilitated the development of diverse molecular and functional diversification patterns to accompany the successful radiation of plants into diverse environments during evolution.

Within-host evolution of Staphylococcus aureus during asymptomatic carriage

BACKGROUND

Staphylococcus aureus is a major cause of healthcare associated mortality, but like many important bacterial pathogens, it is a common constituent of the normal human body flora. Around a third of healthy adults are carriers. Recent evidence suggests that evolution of S. aureus during nasal carriage may be associated with progression to invasive disease. However, a more detailed understanding of within-host evolution under natural conditions is required to appreciate the evolutionary and mechanistic reasons why commensal bacteria such as S. aureus cause disease. Therefore we examined in detail the evolutionary dynamics of normal, asymptomatic carriage. Sequencing a total of 131 genomes across 13 singly colonized hosts using the Illumina platform, we investigated diversity, selection, population dynamics and transmission during the short-term evolution of S. aureus.

PRINCIPAL FINDINGS

We characterized the processes by which the raw material for evolution is generated: micro-mutation (point mutation and small insertions/deletions), macro-mutation (large insertions/deletions) and the loss or acquisition of mobile elements (plasmids and bacteriophages). Through an analysis of synonymous, non-synonymous and intergenic mutations we discovered a fitness landscape dominated by purifying selection, with rare examples of adaptive change in genes encoding surface-anchored proteins and an enterotoxin. We found evidence for dramatic, hundred-fold fluctuations in the size of the within-host population over time, which we related to the cycle of colonization and clearance. Using a newly-developed population genetics approach to detect recent transmission among hosts, we revealed evidence for recent transmission between some of our subjects, including a husband and wife both carrying populations of methicillin-resistant S. aureus (MRSA).

SIGNIFICANCE

This investigation begins to paint a picture of the within-host evolution of an important bacterial pathogen during its prevailing natural state, asymptomatic carriage. These results also have wider significance as a benchmark for future systematic studies of evolution during invasive S. aureus disease.

Ancestral polymorphisms in Drosophila pseudoobscura and Drosophila miranda

Ancestral polymorphisms are defined as variants that arose by mutation prior to the speciation event that generated the species in which they segregate. Their presence may complicate the interpretation of molecular data and lead to incorrect phylogenetic inferences. They may also be used to identify regions of the genome that are under balancing selection. It is thus important to take into account the contribution of ancestral polymorphisms to variability within species and divergence between species. Here, we extend and improve a method for estimation of the proportion of ancestral polymorphisms within a species, and apply it to a dataset of 33 X-linked and 34 autosomal protein-coding genes for which sequence polymorphism data are available in both Drosophila pseudoobscura and Drosophila miranda, using Drosophila affinis as an outgroup. We show that a substantial proportion of both X-linked and autosomal synonymous variants in these two species are ancestral, and that a small number of additional genes with unusually high sequence diversity seem to have an excess of ancestral polymorphisms, suggestive of balancing selection.

RPW8 and resistance to powdery mildew pathogens in natural populations of Arabidopsis lyrata

It is not clear to what extent the orthologues of genes that are adaptively important in one species also contribute to adaptive variation in others. Here, we examine Arabidopsis lyrata to assess the functional and evolutionary significance of natural variation in an orthologue of the gene RPW8 known to be a major determinant of powdery mildew resistance in Arabidopsis thaliana. We assessed the sequence variation at RPW8 and the associated resistance reaction in populations of A. lyrata ssp. petraea. Neutrality tests were performed to understand the importance of local adaptation in maintaining variation at the locus. Highly truncated RPW8 proteins were frequent in all populations and were associated with an increased risk of susceptibility. Haplotypes encoding full-length proteins were highly significantly associated with resistance. There were no signatures of selection at the species-wide level, but some evidence for positive selection in two populations. RPW8 in A. lyrata appears to have a role in powdery mildew resistance, similar to its orthologue in A. thaliana. Unlike A. thaliana, A. lyrata contains a genetic component that can act independently of RPW8 to confer resistance to powdery mildew pathogens. Infrequent local selective sweeps may favour different alleles in different populations, and thereby contribute to the maintenance of species-wide variation at the locus.

Whole genome linkage disequilibrium maps in cattle

Background

Bovine whole genome linkage disequilibrium maps were constructed for eight breeds of cattle. These data provide fundamental information concerning bovine genome organization which will allow the design of studies to associate genetic variation with economically important traits and also provides background information concerning the extent of long range linkage disequilibrium in cattle.

Results

Linkage disequilibrium was assessed using r² among all pairs of syntenic markers within eight breeds of cattle from the Bos taurus and Bos indicus subspecies. Bos taurus breeds included Angus, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black and Limousin while Bos indicus breeds included Brahman and Nelore. Approximately 2670 markers spanning the entire bovine autosomal genome were used to estimate pairwise r² values. We found that the extent of linkage disequilibrium is no more than 0.5 Mb in these eight breeds of cattle.

Conclusion

Linkage disequilibrium in cattle has previously been reported to extend several tens of centimorgans. Our results, based on a much larger sample of marker loci and across eight breeds of cattle indicate that in cattle linkage disequilibrium persists over much more limited distances. Our findings suggest that 30,000–50,000 loci will be needed to conduct whole genome association studies in cattle.

Cladistic approaches to identifying determinants of variability in multifactorial phenotypes and the evolutionary significance of variation in the human genome

Genetic surveys based on detailed restriction site mapping or DNA sequencing allow one to identify many different classes of mutational change at the molecular level and to estimate the evolutionary history of the genetic variation (a haplotype tree). These two sources of information can be combined in a powerful fashion to test hypotheses about the evolutionary significance of genetic variation and to identify mutations that are associated with diseases. Hypotheses about selection on various classes of genetic variation can be tested by examining the distribution patterns of different mutational classes upon the haplotype tree. The power of this procedure can be enhanced if it is coupled with comparative data from other, closely related species. With respect to disease associations, all mutations that affect phenotypic variation in a population occurred at some point in the evolutionary history of the region of the gene containing the mutations. Even if this evolutionary history is estimated from mutations other than those causing phenotypic effects, the phenotypically important mutations are imbedded in this same evolutionary history. Hence, whole branches (clades) of the haplotype tree should display homogeneous phenotypic effects and this fact is utilized to search for phenotypic associations of haplotypes by using nested clades in a haplotype tree. This procedure has more power than alternatives that do not use evolutionary history, and it avoids several statistical and interpretative problems associated with single-marker analyses. All of these methods could be used more extensively if more human genetic surveys concentrated on greater genetic resolution in small DNA regions and included non-human apes.

A simple and robust statistical test for detecting the presence of recombination

Recombination is a powerful evolutionary force that merges historically distinct genotypes. But the extent of recombination within many organisms is unknown, and even determining its presence within a set of homologous sequences is a difficult question. Here we develop a new statistic, phi(w), that can be used to test for recombination. We show through simulation that our test can discriminate effectively between the presence and absence of recombination, even in diverse situations such as exponential growth (star-like topologies) and patterns of substitution rate correlation. A number of other tests, Max chi2, NSS, a coalescent-based likelihood permutation test (from LDHat), and correlation of linkage disequilibrium (both r2 and /D'/) with distance, all tend to underestimate the presence of recombination under strong population growth. Moreover, both Max chi2 and NSS falsely infer the presence of recombination under a simple model of mutation rate correlation. Results on empirical data show that our test can be used to detect recombination between closely as well as distantly related samples, regardless of the suspected rate of recombination. The results suggest that phi(w) is one of the best approaches to distinguish recurrent mutation from recombination in a wide variety of circumstances.

Identification of one intron loss and phylogenetic evolution of Dfak gene in the Drosophila melanogaster species group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5-10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.

Estimating selection on nonsynonymous mutations

The distribution of mutational effects on fitness is of fundamental importance for many aspects of evolution. We develop two methods for characterizing the fitness effects of deleterious, nonsynonymous mutations, using polymorphism data from two related species. These methods also provide estimates of the proportion of amino acid substitutions that are selectively favorable, when combined with data on between-species sequence divergence. The methods are applicable to species with different effective population sizes, but that share the same distribution of mutational effects. The first, simpler, method assumes that diversity for all nonneutral mutations is given by the value under mutation-selection balance, while the second method allows for stronger effects of genetic drift and yields estimates of the parameters of the probability distribution of mutational effects. We apply these methods to data on populations of Drosophila miranda and D. pseudoobscura and find evidence for the presence of deleterious nonsynonymous mutations, mostly with small heterozygous selection coefficients (a mean of the order of 10(-5) for segregating variants). A leptokurtic gamma distribution of mutational effects with a shape parameter between 0.1 and 1 can explain observed diversities, in the absence of a separate class of completely neutral nonsynonymous mutations. We also describe a simple approximate method for estimating the harmonic mean selection coefficient from diversity data on a single species.

Mechanisms of genetic exchange within the chromosomal inversions of Drosophila pseudoobscura

We have used the inversion system of Drosophila pseudoobscura to investigate how genetic flux occurs among the gene arrangements. The patterns of nucleotide polymorphism at seven loci were used to infer gene conversion events between pairs of different gene arrangements. We estimate that the average gene conversion tract length is 205 bp and that the average conversion rate is 3.4 x 10(-6), which is 2 orders of magnitude greater than the mutation rate. We did not detect gene conversion events between all combinations of gene arrangements even though there was sufficient nucleotide variation for detection and sufficient opportunity for exchanges to occur. Genetic flux across the inverted chromosome resulted in higher levels of differentiation within 0.1 Mb of inversion breakpoints, but a slightly lower level of differentiation in central inverted regions. No gene conversion events were detected within 17 kb of an inversion breakpoint suggesting that the formation of double-strand breaks is reduced near rearrangement breakpoints in heterozygotes. At least one case where selection rather than proximity to an inversion breakpoint is responsible for reduction in polymorphism was identified.

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

BACKGROUND

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

RESULTS

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

CONCLUSION

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.

Nucleotide variation in the Egfr locus of Drosophila melanogaster

The Epidermal growth factor receptor is an essential gene with diverse pleiotropic roles in development throughout the animal kingdom. Analysis of sequence diversity in 10.9 kb covering the complete coding region and 6.4 kb of potential regulatory regions in a sample of 250 alleles from three populations of Drosophila melanogaster suggests that the intensity of different population genetic forces varies along the locus. A total of 238 independent common SNPs and 20 indel polymorphisms were detected, with just six common replacements affecting >1475 amino acids, four of which are in the short alternate first exon. Sequence diversity is lowest in a 2-kb portion of intron 2, which is also highly conserved in comparison with D. simulans and D. pseudoobscura. Linkage disequilibrium decays to background levels within 500 bp of most sites, so haplotypes are generally restricted to up to 5 polymorphisms. The two North American samples from North Carolina and California have diverged in allele frequency at a handful of individual SNPs, but a Kenyan sample is both more divergent and more polymorphic. The effect of sample size on inference of the roles of population structure, uneven recombination, and weak selection in patterning nucleotide variation in the locus is discussed.

Nucleotide variation of seven genes in Drosophila kikkawai

We examined levels and patterns of nucleotide variation in 21 strains of Drosophila kikkawai from Miyako island, Japan for the partial regions of the following seven nuclear genes: Adh, Ddc, esc, ksr, Pgi, su(f), and Tpi. The nucleotide variation at total sites (pi(t)) ranged from 0.0013 in the ksr, to 0.0173 in the Adh. The nucleotide divergence at total sites (K(t)) between D. kikkawai and D. lini ranged from 0.0286 in the Tpi to 0.0687 in the su(f). The levels of nucleotide polymorphism and divergence were heterogeneous among the investigated gene regions. The HKA test, which tests imbalance between the intra and interspecific nucleotide variation, showed that the intraspecific nucleotide variation in the Pgi region was much lower than the interspecific variation, while intraspecific variation in the Tpi region was only slightly lower than interspecific variation. The MK test showed an excess of low frequency replacement polymorphic changes in the Adh region, suggesting that most replacement mutations are deleterious. Fay and Wu's test detected an excess of newly arisen variants in the Ddc region. In total, four of the seven gene regions showed significant deviation from the neutrality.

Nucleotide variation in the tinman and bagpipe homeobox genes of Drosophila melanogaster

The tinman (tin) and bagpipe (bap) genes are members of the NK homeobox gene family of Drosophila, so that tin occupies a higher position than bap in the regulatory hierarchy. Little is known about the level and pattern of genetic polymorphism in homeobox genes. We have analyzed nucleotide polymorphism in 27 strains of Drosophila melanogaster and one each of D. simulans and D. sechellia, within two closely linked regions encompassing a partial sequence of tin and the complete sequence of bap. The two genes exhibit different levels and patterns of nucleotide diversity. Two sets of sharply divergent sequence types are detected for tin. The haplotype structure of bap is more complex: about half of the sequences are identical (or virtually so), while the rest are fairly heterogeneous. The level of silent nucleotide variability is 0.0063 for tin but significantly higher, 0.0141, for bap, a level of polymorphism comparable to the most polymorphic structural genes of D. melanogaster. Recombination rate and gene conversion are also higher for bap than for tin. There is strong linkage disequilibrium, with the highest values in the introns of both genes and exon II of bap. The patterns of polymorphism in tin and bap are not compatible with an equilibrium model of selective neutrality. We suggest that negative selection and demographic history are the major factors shaping the pattern of nucleotide polymorphism in the tin and bap genes; moreover, there are clear indications of positive selection in the bap gene.

Effects of population structure and sex on association between serotonin receptors and Drosophila heart rate

As a first step toward population and quantitative genetic analysis of neurotransmitter receptors in Drosophila melanogaster, we describe the parameters of nucleotide variation in three serotonin receptors and their association with pupal heart rate. Thirteen kilobases of DNA including the complete coding regions of 5-HT1A, 5-HT1B, and 5-HT2 were sequenced in 216 highly inbred lines extracted from two North American populations in California and North Carolina. Nucleotide and amino acid polymorphism is in the normal range for Drosophila genes and proteins, and linkage disequilibrium decays rapidly such that haplotype blocks are typically only a few SNPs long. However, intron 1 of 5-HT1A consists of two haplotypes that are at significantly different frequencies in the two populations. Neither this region of the gene nor any of the common amino acid polymorphisms in the three loci associate with either heart rate or heart rate variability. A cluster of SNPs in intron 2 of 5-HT1A, including a triallelic site, do show a highly significant interaction between genotype, sex, and population. While it is likely that a combination of weak, complex selection pressures and population structure has helped shape variation in the serotonin receptors of Drosophila, much larger sampling strategies than are currently adopted in evolutionary genetics will be required to disentangle these effects.

Pleiotropic effect of disrupting a conserved sequence involved in a long-range compensatory interaction in the Drosophila Adh gene

Recent advances in experimental analyses of the evolution of RNA secondary structures suggest a more complex scenario than that typically considered by Kimura's classical model of compensatory evolution. In this study, we examine one such case in more detail. Previous experimental analysis of long-range compensatory interactions between the two ends of Drosophila Adh mRNA failed to fit the classical model of compensatory evolution. To further investigate and verify long-range pairing in Drosophila Adh with respect to models of compensatory evolution and its potential functional role, we introduced site-directed mutations in the Drosophila melanogaster Adh gene. We explore two alternative hypotheses for why previous analysis of long-range compensatory interactions failed to fit the classical model. Specifically, we investigate whether the disruption of a conserved short-range pairing within Adh exon 2 has an effect on Adh expression or if there is a dual functional role of a conserved sequence in the 3'-UTR in both long-range pairing and the negative regulation of Adh expression. We find that a classical result was not observed due to the pleiotropic effect of changing a nucleotide involved in both long-range base pairing and the negative regulation of gene expression.

Compensatory evolution in the human malaria parasite Plasmodium ovale

The fixation of neutral compensatory mutations in a population depends on the effective population size of the species, which can fluctuate dramatically within a few generations, the mutation rate, and the selection intensity associated with the individual mutations. We observe compensatory mutations and intermediate states in populations of the malaria parasite Plasmodium ovale. The appearance of compensatory mutations and intermediate states in P. ovale raises interesting questions about population structure that could have considerable impact on the control of the associated disease.

Low linkage disequilibrium indicative of recombination in foot-and-mouth disease virus gene sequence alignments

We have applied tests for detecting recombination to genes of foot-and-mouth disease virus (FMDV). Our approach estimated summary statistics of linkage disequilibrium (LD), which are sensitive to recombination. Using the genealogical relationships, rate heterogeneity and mutation parameters estimated from individual sets of aligned gene sequences, we simulated matching RNA sequence datasets without recombination. These simulated datasets allowed for recurrent mutations at any site to mimic homoplasy in virus sequence data and allow construction of null distributions for LD parameters expected in the absence of recombination. We tested for recombination in two ways: by comparing LD in observed data with corresponding null distributions obtained from simulated data; and by testing for a negative relationship between observed LD between pairs of polymorphic nucleotide sites and inter-site distance. We applied these tests to six FMDV datasets from four serotypes and found some evidence for recombination in all of them.

Nucleotide Variation in the tinman and bagpipe Homeobox Genes of Drosophila melanogaster

The tinman (tin) and bagpipe (bap) genes are members of the NK homeobox gene family of Drosophila, so that tin occupies a higher position than bap in the regulatory hierarchy. Little is known about the level and pattern of genetic polymorphism in homeobox genes. We have analyzed nucleotide polymorphism in 27 strains of Drosophila melanogaster and one each of D. simulans and D. sechellia, within two closely linked regions encompassing a partial sequence of tin and the complete sequence of bap. The two genes exhibit different levels and patterns of nucleotide diversity. Two sets of sharply divergent sequence types are detected for tin. The haplotype structure of bap is more complex: about half of the sequences are identical (or virtually so), while the rest are fairly heterogeneous. The level of silent nucleotide variability is 0.0063 for tin but significantly higher, 0.0141, for bap, a level of polymorphism comparable to the most polymorphic structural genes of D. melanogaster. Recombination rate and gene conversion are also higher for bap than for tin. There is strong linkage disequilibrium, with the highest values in the introns of both genes and exon II of bap. The patterns of polymorphism in tin and bap are not compatible with an equilibrium model of selective neutrality. We suggest that negative selection and demographic history are the major factors shaping the pattern of nucleotide polymorphism in the tin and bap genes; moreover, there are clear indications of positive selection in the bap gene.

Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene

Evidence for the evolutionary maintenance of a hairpin structure possibly involved in intron processing had been found in intron 1 of the alcohol dehydrogenase gene (Adh) in diverse Drosophila species. In this study, the putative hairpin structure was evaluated systematically in Drosophila melanogaster by elimination of either side of the stem using site-directed mutagenesis. The effects of these mutations and the compensatory double mutant on intron splicing efficiency and ADH protein production were assayed in Drosophila melanogaster Schneider L2 cells and germ-line transformed adult flies. Mutations that disrupt the putative hairpin structure right upstream of the intron branch point were found to cause a significant reduction in both splicing efficiency and ADH protein production. In contrast, the compensatory double mutant that restores the putative hairpin structure was indistinguishable from the WT in both splicing efficiency and ADH level. It was also observed by mutational analysis that a more stable secondary structure (with a longer stem) in this intron decreases both splicing efficiency and ADH protein production. Implications for RNA secondary structure and intron evolution are discussed.

Evolutionary genomics of inversions in Drosophila pseudoobscura: evidence for epistasis

Drosophila pseudoobscura harbors a rich polymorphism for paracentric inversions on the third chromosome, and the clines in the inversion frequencies across the southwestern United States indicate that strong natural selection operates on them. Isogenic inversion strains were made from isofemale lines collected from four localities, and eight molecular markers were mapped on the third chromosome. Nucleotide diversity was measured for these loci and formed the basis of an evolutionary genomic analysis. The loci were differentiated among inversions. The inversions did not show significant differences among populations, however, likely the result of extensive gene flow among populations. Some loci had significant reductions in nucleotide diversity within inversions compared with interspecies divergence, suggesting that these loci are near inversion breakpoints or are near targets of directional selection. Linkage disequilibrium (LD) levels tended to decrease with distance between loci, indicating that some genetic exchange occurs among gene arrangements despite the presence of inversions. In some cases, however, adjacent genes had low levels of interlocus LD and loosely linked genes had high levels of interlocus LD, suggesting strong epistatic selection. Our results support the hypothesis that the inversions of D. pseudoobscura have emerged as suppressors of recombination to maintain positive epistatic relationships among loci within gene arrangements that developed as the species adapted to a heterogeneous environment.

Sequence variation of alcohol dehydrogenase (Adh) paralogs in cactophilic Drosophila

This study focuses on the population genetics of alcohol dehydrogenase (Adh) in cactophilic Drosophila. Drosophila mojavensis and D. arizonae utilize cactus hosts, and each host contains a characteristic mixture of alcohol compounds. In these Drosophila species there are two functional Adh loci, an adult form (Adh-2) and a larval and ovarian form (Adh-1). Overall, the greater level of variation segregating in D. arizonae than in D. mojavensis suggests a larger population size for D. arizonae. There are markedly different patterns of variation between the paralogs across both species. A 16-bp intron haplotype segregates in both species at Adh-2, apparently the product of an ancient gene conversion event between the paralogs, which suggests that there is selection for the maintenance of the intron structure possibly for the maintenance of pre-mRNA structure. We observe a pattern of variation consistent with adaptive protein evolution in the D. mojavensis lineage at Adh-1, suggesting that the cactus host shift that occurred in the divergence of D. mojavensis from D. arizonae had an effect on the evolution of the larval expressed paralog. Contrary to previous work we estimate a recent time for both the divergence of D. mojavensis and D. arizonae (2.4 +/- 0.7 MY) and the age of the gene duplication (3.95 +/- 0.45 MY).

Intron presence-absence polymorphism in Drosophila driven by positive Darwinian selection

Comparisons of intron-exon structures between homologous genes in different eukaryotic species have revealed substantial variation in the number of introns. These observations imply that, in each case, an intron presence-absence polymorphism must have existed in the past. Such a polymorphism, created by a recent intron-loss mutation, is reported here in a eukaryotic organism. This gene structure, detected in the jingwei (jgw) gene, segregates at high frequency (77%) in natural populations of Drosophila teissieri and is associated with a marked change in mRNA levels. Furthermore, the intron loss does not result from a mRNA-mediated mechanism as is usually proposed, but from a partial deletion at the DNA level that also results in the addition of four new amino acids to the JGW protein. Population genetic analyses of the pattern of nucleotide variation surrounding the intron polymorphism indicate the action of positive Darwinian selection on the intron-absent variant. Forward simulations suggest that the intensity of this selection is weak to moderate, roughly equal to the selection intensity on most replacement mutations in Drosophila.

Distinguishing between selection and population expansion in an experimental lineage of bacteriophage T7

Experimental evolution of short-lived organisms offers the opportunity to study the dynamics of polymorphism over time in a controlled environment. Here, we characterize DNA polymorphism data over time for four genes in bacteriophage T7. Our experiment ran for 2500 generations and populations were sampled after 500, 2000, and 2500 generations. We detect positive selection, purifying ("negative") selection, and population expansion in our experiment. We also present a statistical test that is able to distinguish demographic from selective events, processes that are hard to identify individually because both often produce an excess of rare mutations. Our "heterogeneity test" modifies common statistics measuring the frequency spectrum of polymorphism (e.g., Fu and Li's D) by looking for processes producing different patterns on nonsynonymous and synonymous mutations. Test results agree with the known conditions of the experiment, and we are therefore confident that this test offers a tool to evaluate natural populations. Our results suggest that instances of segregating deleterious mutations may be common, but as yet undetected, in nature.

Nucleotide variation of the duplicated amylase genes in Drosophila kikkawai

We examined levels and patterns of the nucleotide polymorphism of the Amylase genes with a head-to-head duplication in Drosophila kikkawai. The levels of variation in D. kikkawai were comparable to those in Drosophila melanogaster. Tajima's test, Fu and Li's test, HKA test, and MK test did not show significant departure from neutrality. We found an excess of replacement changes in the within-locus class, representing polymorphism in one of the duplicated genes, compared with the between-locus class, representing polymorphism shared between the duplicated genes. Most replacement changes in the within-locus class were singletons. These results suggest that most replacement changes are deleterious. A contrasting evolutionary pattern, involving concerted evolution in the coding regions but differential evolution in the 5'-flanking regions, was observed. However, unlike the duplicated Amy genes of D. melanogaster, the coding regions of the duplicated genes in D. kikkawai tended to diverge. Using Ohta's model of the small multigene family, we found that recombination (interchromosomal equal crossing-over) rate was one order higher than gene conversion (unequal crossing-over) rate, resulting in a considerable but incomplete homogenization of the duplicated coding regions. Linkage disequilibria were found in the intron as well as within and around the regulatory cis-element sequences of one of the duplicated genes (Amy1). The possible causes of these linkage disequilibria were discussed.

Patterns of linkage disequilibrium in the human genome

Particular alleles at neighbouring loci tend to be co-inherited. For tightly linked loci, this might lead to associations between alleles in the population a property known as linkage disequilibrium (LD). LD has recently become the focus of intense study in the hope that it might facilitate the mapping of complex disease loci through whole-genome association studies. This approach depends crucially on the patterns of LD in the human genome. In this review, we draw on empirical studies in humans and Drosophila, as well as simulation studies, to assess the current state of knowledge about patterns of LD, and consider the implications for the use of LD as a mapping tool.

Testing for recombination in a short nuclear DNA sequence of the European meadow grasshopper, Chorthippus parallelus

Single-copy nuclear DNA sequences have high potential as a source of genetic markers for population analyses. However, the difficulties that arise when haplotypes that are the product of recombinational rearrangements are present require additional consideration. Two statistical methods for identifying potential recombinants by detecting anomalies in the distribution of variable sites along sequences were used to screen sequences from a single-copy nuclear DNA fragment, cpnl-1, of the European meadow grasshopper (Chorthippus parallelus). Five of the 71 haplotypes in the cpnl-1 data set showed nonrandom distribution of polymorphic sites using both methods. The second method pinpointed an additional four haplotypes. Estimates of the rate of recombination in the entire data set were obtained using standard methods. It is concluded that cpnl-1 haplotypes have been involved in recombination or gene conversion events at a rate more than twice the mutation rate. This confirms that recombination and gene conversion are significant factors in the generation of haplotype variation in nuclear gene sequences. The cpnl-1 haplotypes identified by the tests were present only in populations that have had recent contact; the Balkan and Turkish refugial populations and their post-glacial colonies to the north. This is discussed in relation to the phylogenetic inferences drawn from the same data in a previous report.

Protein variation in Adh and Adh-related in Drosophila pseudoobscura. Linkage disequilibrium between single nucleotide polymorphisms and protein alleles

A 3.5-kb segment of the alcohol dehydrogenase (Adh) region that includes the Adh and Adh-related genes was sequenced in 139 Drosophila pseudoobscura strains collected from 13 populations. The Adh gene encodes four protein alleles and rejects a neutral model of protein evolution with the McDonald-Kreitman test, although the number of segregating synonymous sites is too high to conclude that adaptive selection has operated. The Adh-related gene encodes 18 protein haplotypes and fails to reject an equilibrium neutral model. The populations fail to show significant geographic differentiation of the Adh-related haplotypes. Eight of 404 single nucleotide polymorphisms (SNPs) in the Adh region were in significant linkage disequilibrium with three ADHR protein alleles. Coalescent simulations with and without recombination were used to derive the expected levels of significant linkage disequilibrium between SNPs and 18 protein haplotypes. Maximum levels of linkage disequilibrium are expected for protein alleles at moderate frequencies. In coalescent models without recombination, linkage disequilibrium decays between SNPs and high frequency haplotypes because common alleles mutate to haplotypes that are rare or that reach moderate frequency. The implication of this study is that linkage disequilibrium mapping has the highest probability of success with disease-causing alleles at frequencies of 10%.

Selection intensity against deleterious mutations in RNA secondary structures and rate of compensatory nucleotide substitutions

A two-locus model of reversible mutations with compensatory fitness interactions is presented; single mutations are assumed to be deleterious but neutral in appropriate combinations. The expectation of the time of compensatory nucleotide substitutions is calculated analytically for the case of tight linkage between sites. It is shown that selection increases the substitution time dramatically when selection intensity Ns > 1, where N is the diploid population size and s the selection coefficient. Computer simulations demonstrate that recombination increases the substitution time, but the effect of recombination is small when selection is weak. The amount of linkage disequilibrium generated in the process of compensatory substitution is also investigated. It is shown that significant linkage disequilibrium is expected to be rare in natural populations. The model is applied to the mRNA secondary structure of the bicoid 3' untranslated region of Drosophila. It is concluded that average selection intensity Ns against single deleterious mutations is not likely to be much larger than 1.

Genetic trends in a population evolving antibiotic resistance

The evolution of antibiotic resistance provides a well-documented, rapid, and recent example of a selection driven process that has occurred in many bacterial species. An exhaustive collection of Moraxella catarrhalis that spans a transition to chromosomally encoded penicillin resistance was used to analyze genetic changes accompanying the transition. The population was characterized by high haplotypic diversity with 148 distinct haplotypes among 372 isolates tested at three genomic regions. The power of a temporally stratified sample from a single population was highlighted by the finding of high genetic diversity throughout the transition to resistance, population numbers that remained high over time, and no evidence of departures from neutrality in the allele frequency spectra throughout the transition. The direct temporal analysis documented the persistence, antibiotic status, and haplotypic identity of strains undergoing apparent clonal expansions. Several haplotypes that were beta-lactamase nonproducers in early samples converted to producers in later years. Maintenance of genetic diversity and haplotype conversions from sensitive to resistant supported the hypothesis that penicillin resistance determinants spread to a diverse array of strains via horizontal exchange. Genetic differentiation between sample years, estimated by F(ST), was increasing at a rate that could cause complete haplotype turnover in less than 150 years. Widespread linkage disequilibrium among sites within one locus (copB) suggested recent mutation followed by clonal expansion. Nonrandom associations between haplotypes and resistance phenotypes provided further evidence of clonal expansion for some haplotypes. Nevertheless, the population structure was far from clonal as evidenced by a relatively low frequency of disequilibria both within sites at a second locus (M46) as well as between loci. The haplotype-antibiotic resistance association that was accompanied by gradual haplotype turnover is consistent with a hypothesis of genetic drift at marker loci with directional selection at the resistance locus.

Do mitochondria recombine in humans?

Until very recently, mitochondria were thought to be clonally inherited through the maternal line in most higher animals. However, three papers published in 2000 claimed population-genetic evidence of recombination in human mitochondrial DNA. Here I review the current state of the debate. I review the evidence for the two main pathways by which recombination might occur: through paternal leakage and via a mitochondrial DNA sequence in the nuclear genome. There is no strong evidence for either pathway, although paternal leakage seems a definite possibility. However, the population-genetic evidence, although not conclusive, is strongly suggestive of recombination in mitochondrial DNA. The implications of non-clonality for our understanding of human and mitochondrial evolution are discussed.

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.

A genome-wide departure from the standard neutral model in natural populations of Drosophila

We analyze nucleotide polymorphism data for a large number of loci in areas of normal to high recombination in Drosophila melanogaster and D. simulans (24 and 16 loci, respectively). We find a genome-wide, systematic departure from the neutral expectation for a panmictic population at equilibrium in natural populations of both species. The distribution of sequence-based estimates of 2Nc across loci is inconsistent with the assumptions of the standard neutral theory, given the observed levels of nucleotide diversity and accepted values for recombination and mutation rates. Under these assumptions, most estimates of 2Nc are severalfold too low; in other words, both species exhibit greater intralocus linkage disequilibrium than expected. Variation in recombination or mutation rates is not sufficient to account for the excess of linkage disequilibrium. While an equilibrium island model does not seem to account for the data, more complicated forms of population structure may. A proper test of alternative demographic models will require loci to be sampled in a more consistent fashion.

Molecular population genetics of X-linked genes in Drosophila pseudoobscura

This article presents a nucleotide sequence analysis of 500 bp determined in each of five X-linked genes, runt, sisterlessA, period, esterase 5, and Heat-shock protein 83, in 40 Drosophila pseudoobscura strains collected from two populations. Estimates of the neutral migration parameter for the five loci show that gene flow among D. pseudoobscura populations is sufficient to homogenize inversion frequencies across the range of the species. Nucleotide diversity at each locus fails to reject a neutral model of molecular evolution. The sample of 40 chromosomes included six Sex-ratio inversions, a series of three nonoverlapping inversions that are associated with a strong meiotic drive phenotype. The selection driven by the Sex-ratio meiotic drive element has not fixed variation across the X chromosome of D. pseudoobscura because, while significant linkage disequilibrium was observed within the sisterlessA, period, and esterase 5 genes, we did not find evidence for nonrandom association among loci. The Sex-ratio chromosome was estimated to be 25,000 years old based on the decomposition of linkage disequilibrium between esterase 5 and Heat-shock protein 83 or 1 million years old based on the net divergence of esterase 5 between Standard and Sex-ratio chromosomes. Genetic diversity was depressed within esterase 5 within Sex-ratio chromosomes, while the four other genes failed to show a reduction in heterozygosity in the Sex-ratio background. The reduced heterogeneity in esterase 5 is due either to its location near one of the Sex-ratio inversion breakpoints or that it is closely linked to a gene or genes responsible for the Sex-ratio meiotic drive system.

Molecular evolution near a two-locus balanced polymorphism

Balancing selection at one locus can increase the amount of selectively neutral variation within neighboring genomic regions. Discrete phenotypic polymorphisms studied in natural populations are frequently determined by sets of interacting genes instead of alternative alleles at single loci. We extend coalescent theory to investigate balancing selection on combinations of linked genes. We find that variation at neutral sites is increased across a much larger genomic region relative to the single-locus models: the entire region lying between the two loci in balanced combination is affected to some degree. Epistatic selection maintains these high levels of neutral variation because it directly opposes the homogenizing effect of recombination. The results of the theory are discussed in relation to published gene sequence data, primarily from Drosophila.

Contrasting patterns of molecular evolution of the genes on the new and old sex chromosomes of Drosophila miranda

In organisms with chromosomal sex determination, sex is determined by a set of dimorphic sex chromosomes that are thought to have evolved from a set of originally homologous chromosomes. The chromosome inherited only through the heterogametic sex (the Y chromosome in the case of male heterogamety) often exhibits loss of genetic activity for most of the genes carried on its homolog and is hence referred to as degenerate. The process by which the proto-Y chromosome loses its genetic activity has long been the subject of much speculation. We present a DNA sequence variation analysis of marker genes on the evolving sex chromosomes (neo-sex chromosomes) of Drosophila miranda. Due to its relatively recent origin, the neo-Y chromosome of this species is presumed to be still experiencing the forces responsible for the loss of its genetic activity. Indeed, several previous studies have confirmed the presence of some active loci on this chromosome. The genes on the neo-Y chromosome surveyed in the current study show generally lower levels of variation compared with their counterparts on the neo-X chromosome or an X-linked gene. This is in accord with a reduced effective population size of the neo-Y chromosome. Interestingly, the rate of replacement nucleotide substitutions for the neo-Y linked genes is significantly higher than that for the neo-X linked genes. This is not expected under a model where the faster evolution of the X chromosome is postulated to be the main force driving the degeneration of the Y chromosome.

Linkage disequilibrium and recombination in hominid mitochondrial DNA

The assumption that human mitochondrial DNA is inherited from one parent only and therefore does not recombine is questionable. Linkage disequilibrium in human and chimpanzee mitochondrial DNA declines as a function of the distance between sites. This pattern can be attributed to one mechanism only: recombination.

RNA secondary structure and compensatory evolution

The classic concept of epistatic fitness interactions between genes has been extended to study interactions within gene regions, especially between nucleotides that are important in maintaining pre-mRNA/mRNA secondary structures. It is shown that the majority of linkage disequilibria found within the Drosophila Adh gene are likely to be caused by epistatic selection operating on RNA secondary structures. A recently proposed method of RNA secondary structure prediction based on DNA sequence comparisons is reviewed and applied to several types of RNAs, including tRNA, rRNA, and mRNA. The patterns of covariation in these RNAs are analyzed based on Kimura's compensatory evolution model. The results suggest that this model describes the substitution process in the pairing regions (helices) of RNA secondary structures well when the helices are evolutionarily conserved and thermodynamically stable, but fails in some other cases. Epistatic selection maintaining pre-mRNA/mRNA secondary structures is compared to weak selective forces that determine features such as base composition and synonymous codon usage. The relationships among these forces and their relative strengths are addressed. Finally, our mutagenesis experiments using the Drosophila Adh locus are reviewed. These experiments analyze long-range compensatory interactions between the 5' and 3' ends of Adh mRNA, the different constraints on secondary structures in introns and exons, and the possible role of secondary structures in RNA splicing.

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.

Recombination and selection at Brassica self-incompatibility loci

In Brassica species, self-incompatibility is controlled genetically by haplotypes involving two known genes, SLG and SRK, and possibly an as yet unknown gene controlling pollen incompatibility types. Alleles at the incompatibility loci are maintained by frequency-dependent selection, and diversity at SLG and SRK appears to be very ancient, with high diversity at silent and replacement sites, particularly in certain "hypervariable" portions of the genes. It is important to test whether recombination occurs in these genes before inferences about function of different parts of the genes can be made from patterns of diversity within their sequences. In addition, it has been suggested that, to maintain the relationship between alleles within a given S-haplotype, recombination is suppressed in the S-locus region. The high diversity makes many population genetic measures of recombination inapplicable. We have analyzed linkage disequilibrium within the SLG gene of two Brassica species, using published coding sequences. The results suggest that intragenic recombination has occurred in the evolutionary history of these alleles. This is supported by patterns of synonymous nucleotide diversity within both the SLG and SRK genes, and between domains of the SRK gene. Finally, clusters of linkage disequilibrium within the SLG gene suggest that hypervariable regions are under balancing selection, and are not merely regions of relaxed selective constraint.

High recombination rate in natural populations of Plasmodium falciparum

Malaria parasites are sexually reproducing protozoa, although the extent of effective meiotic recombination in natural populations has been debated. If meiotic recombination occurs frequently, compared with point mutation and mitotic rearrangement, linkage disequilibrium between polymorphic sites is expected to decline with increasing distance along a chromosome. The rate of this decline should be proportional to the effective meiotic recombination rate in the population. Multiple polymorphic sites covering a 5-kb region of chromosome 9 (the msp1 gene) have been typed in 547 isolates from six populations in Africa to test for such a decline and estimate its rate in populations of Plasmodium falciparum. The magnitude of two-site linkage disequilibrium declines markedly with increasing molecular map distance between the sites, reaching nonsignificant levels within a map range of 0.3-1.0 kb in five of the populations and over a larger map distance in the population with lowest malaria endemicity. The rate of decline in linkage disequilibrium over molecular map distance is at least as rapid as that observed in most chromosomal regions of other sexually reproducing eukaryotes, such as humans and Drosophila. These results are consistent with the effective recombination rate expected in natural populations of P. falciparum, predicted on the basis of the underlying molecular rate of meiotic crossover and the coefficient of inbreeding caused by self-fertilization events. This is conclusive evidence to reject any hypothesis of clonality or low rate of meiotic recombination in P. falciparum populations. Moreover, the data have major implications for the design and interpretation of population genetic studies of selection on P. falciparum genes.

A highly conserved sequence in the 3'-untranslated region of the drosophila Adh gene plays a functional role in Adh expression

Phylogenetic analysis identified a highly conserved eight-base sequence (AAGGCTGA) within the 3'-untranslated region (UTR) of the Drosophila alcohol dehydrogenase gene, Adh. To examine the functional significance of this conserved motif, we performed in vitro deletion mutagenesis on the D. melanogaster Adh gene followed by P-element-mediated germline transformation. Deletion of all or part of the eight-base sequence leads to a twofold increase in in vivo ADH enzymatic activity. The increase in activity is temporally and spatially general and is the result of an underlying increase in Adh transcript. These results indicate that the conserved 3'-UTR motif plays a functional role in the negative regulation of Adh gene expression. The evolutionary significance of our results may be understood in the context of the amino acid change that produces the ADH-F allele and also leads to a twofold increase in ADH activity. While there is compelling evidence that the amino acid replacement has been a target of positive selection, the conservation of the 3'-UTR sequence suggests that it is under strong purifying selection. The selective difference between these two sequence changes, which have similar effects on ADH activity, may be explained by different metabolic costs associated with the increase in activity.

DNA variability and recombination rates at X-linked loci in humans

We sequenced 11,365 bp from introns of seven X-linked genes in 10 humans, one chimpanzee, and one orangutan to (i) provide an average estimate of nucleotide diversity (pi) in humans, (ii) investigate whether there is variation in pi among loci, (iii) compare ratios of polymorphism to divergence among loci, and (iv) provide a preliminary test of the hypothesis that heterozygosity is positively correlated with the local rate of recombination. The average value for pi was low 0.063%, SE = 0.036%, about one order of magnitude smaller than for Drosophila melanogaster, the species for which the best data are available. Among loci, pi varied by over one order of magnitude. Statistical tests of neutrality based on ratios of polymorphism to divergence or based on the frequency spectrum of variation within humans failed to reject a neutral, equilibrium model. However, there was a positive correlation between heterozygosity and rate of recombination, suggesting that the joint effects of selection and linkage are important in shaping patterns of nucleotide variation in humans.