Molecular and phenotypic variation of the white locus region in Drosophila melanogaster

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.

Evolution of genome structure in the Drosophila simulans species complex

The rapid evolution of repetitive DNA sequences, including satellite DNA, tandem duplications, and transposable elements, underlies phenotypic evolution and contributes to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented and misassembled in most contemporary genome assemblies. We generated highly contiguous de novo reference genomes for the Drosophila simulans species complex (D. simulans, D. mauritiana, and D. sechellia), which speciated ∼250,000 yr ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive sequences between these four species. We find that at least 15% of the D. simulans complex species genomes fail to align uniquely to D. melanogaster owing to structural divergence-twice the number of single-nucleotide substitutions. We also find rapid turnover of satellite DNA and extensive structural divergence in heterochromatic regions, whereas the euchromatic gene content is mostly conserved. Despite the overall preservation of gene synteny, euchromatin in each species has been shaped by clade- and species-specific inversions, transposable elements, expansions and contractions of satellite and tRNA tandem arrays, and gene duplications. We also find rapid divergence among Y-linked genes, including copy number variation and recent gene duplications from autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.

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.

These aren't the loci you'e looking for: Principles of effective SNP filtering for molecular ecologists

Sequencing reduced-representation libraries of restriction site-associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

Massive habitat-specific genomic response in D. melanogaster populations during experimental evolution in hot and cold environments

Experimental evolution in combination with whole-genome sequencing (evolve and resequence [E&R]) is a promising approach to define the genotype-phenotype map and to understand adaptation in evolving populations. Many previous studies have identified a large number of putative selected sites (i.e., candidate loci), but it remains unclear to what extent these loci are genuine targets of selection or experimental noise. To address this question, we exposed the same founder population to two different selection regimes-a hot environment and a cold environment-and quantified the genomic response in each. We detected large numbers of putative selected loci in both environments, albeit with little overlap between the two sets of candidates, indicating that most resulted from habitat-specific selection. By quantifying changes across multiple independent biological replicates, we demonstrate that most of the candidate SNPs were false positives that were linked to selected sites over distances much larger than the typical linkage disequilibrium range of Drosophila melanogaster. We show that many of these mid- to long-range associations were attributable to large segregating inversions and confirm by computer simulations that such patterns could be readily replicated when strong selection acts on rare haplotypes. In light of our findings, we outline recommendations to improve the performance of future Drosophila E&R studies which include using species with negligible inversion loads, such as D. mauritiana and D. simulans, instead of D. melanogaster.

Long range linkage disequilibrium across the human genome

Long-range linkage disequilibria (LRLD) between sites that are widely separated on chromosomes may suggest that population admixture, epistatic selection, or other evolutionary forces are at work. We quantified patterns of LRLD on a chromosome-wide level in the YRI population of the HapMap dataset of single nucleotide polymorphisms (SNPs). We calculated the disequilibrium between all pairs of SNPs on each chromosome (a total of >2×10(11) values) and evaluated significance of overall disequilibrium using randomization. The results show an excess of associations between pairs of distant sites (separated by >0.25 cM) on all of the 22 autosomes. We discuss possible explanations for this observation.

A genome-wide, fine-scale map of natural pigmentation variation in Drosophila melanogaster

Various approaches can be applied to uncover the genetic basis of natural phenotypic variation, each with their specific strengths and limitations. Here, we use a replicated genome-wide association approach (Pool-GWAS) to fine-scale map genomic regions contributing to natural variation in female abdominal pigmentation in Drosophila melanogaster, a trait that is highly variable in natural populations and highly heritable in the laboratory. We examined abdominal pigmentation phenotypes in approximately 8000 female European D. melanogaster, isolating 1000 individuals with extreme phenotypes. We then used whole-genome Illumina sequencing to identify single nucleotide polymorphisms (SNPs) segregating in our sample, and tested these for associations with pigmentation by contrasting allele frequencies between replicate pools of light and dark individuals. We identify two small regions near the pigmentation genes tan and bric-à-brac 1, both corresponding to known cis-regulatory regions, which contain SNPs showing significant associations with pigmentation variation. While the Pool-GWAS approach suffers some limitations, its cost advantage facilitates replication and it can be applied to any non-model system with an available reference genome.

Phenotypic variation is abundant in natural populations, but its genetic basis is not always well-understood. Here, we examine the genetic basis of body pigmentation in Drosophila, a trait with a long history of study in Drosophila genetics and evolution. We conducted the first genome-wide scan for polymorphism associated with pigmentation variation in a large natural sample of D. melanogaster, and found SNPs near two genes, tan and bric-à-brac 1, affecting the trait. The SNPs associated with pigmentation variation in these genes appear to act by affecting the regulation of the pigmentation genes, rather than their protein coding sequence.

Adaptation of Drosophila to a novel laboratory environment reveals temporally heterogeneous trajectories of selected alleles

The genomic basis of adaptation to novel environments is a fundamental problem in evolutionary biology that has gained additional importance in the light of the recent global change discussion. Here, we combined laboratory natural selection (experimental evolution) in Drosophila melanogaster with genome-wide next generation sequencing of DNA pools (Pool-Seq) to identify alleles that are favourable in a novel laboratory environment and traced their trajectories during the adaptive process. Already after 15 generations, we identified a pronounced genomic response to selection, with almost 5000 single nucleotide polymorphisms (SNP; genome-wide false discovery rates < 0.005%) deviating from neutral expectation. Importantly, the evolutionary trajectories of the selected alleles were heterogeneous, with the alleles falling into two distinct classes: (i) alleles that continuously rise in frequency; and (ii) alleles that at first increase rapidly but whose frequencies then reach a plateau. Our data thus suggest that the genomic response to selection can involve a large number of selected SNPs that show unexpectedly complex evolutionary trajectories, possibly due to nonadditive effects.

Genomic variation in natural populations of Drosophila melanogaster

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5'- and 3'-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

Long-term and short-term evolutionary impacts of transposable elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host-TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host-pathogen coevolution.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Molecular population genomics: a short history

Population genomics is the study of the amount and causes of genome-wide variability in natural populations, a topic that has been under discussion since Darwin. This paper first briefly reviews the early development of molecular approaches to the subject: the pioneering unbiased surveys of genetic variability at multiple loci by means of gel electrophoresis and restriction enzyme mapping. The results of surveys of levels of genome-wide variability using DNA resequencing studies are then discussed. Studies of the extent to which variability for different classes of variants (non-synonymous, synonymous and non-coding) are affected by natural selection, or other directional forces such as biased gene conversion, are also described. Finally, the effects of deleterious mutations on population fitness and the possible role of Hill–Robertson interference in shaping patterns of sequence variability are discussed.

Genealogical lineage sorting leads to significant, but incorrect Bayesian multilocus inference of population structure

Over the past decades, the use of molecular markers has revolutionized biology and led to the foundation of a new research discipline—phylogeography. Of particular interest has been the inference of population structure and biogeography. While initial studies focused on mtDNA as a molecular marker, it has become apparent that selection and genealogical lineage sorting could lead to erroneous inferences. As it is not clear to what extent these forces affect a given marker, it has become common practice to use the combined evidence from a set of molecular markers as an attempt to recover the signals that approximate the true underlying demography. Typically, the number of markers used is determined by either budget constraints or by statistical power required to recognize significant population differentiation. Using microsatellite markers from Drosophila and humans, we show that even large numbers of loci (>50) can frequently result in statistically well-supported, but incorrect inference of population structure using the software baps. Most importantly, genomic features, such as chromosomal location, variability of the markers, or recombination rate, cannot explain this observation. Instead, it can be attributed to sampling variation among loci with different realizations of the stochastic lineage sorting. This phenomenon is particularly pronounced for low levels of population differentiation. Our results have important implications for ongoing studies of population differentiation, as we unambiguously demonstrate that statistical significance of population structure inferred from a random set of genetic markers cannot necessarily be taken as evidence for a reliable demographic inference.

Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster

Determining the genetic basis of environmental adaptation is a central problem of evolutionary biology. This issue has been fruitfully addressed by examining genetic differentiation between populations that are recently separated and/or experience high rates of gene flow. A good example of this approach is the decades-long investigation of selection acting along latitudinal clines in Drosophila melanogaster. Here we use next-generation genome sequencing to reexamine the well-studied Australian D. melanogaster cline. We find evidence for extensive differentiation between temperate and tropical populations, with regulatory regions and unannotated regions showing particularly high levels of differentiation. Although the physical genomic scale of geographic differentiation is small--on the order of gene sized--we observed several larger highly differentiated regions. The region spanned by the cosmopolitan inversion polymorphism In(3R)P shows higher levels of differentiation, consistent with the major difference in allele frequencies of Standard and In(3R)P karyotypes in temperate vs. tropical Australian populations. Our analysis reveals evidence for spatially varying selection on a number of key biological processes, suggesting fundamental biological differences between flies from these two geographic regions.

Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics

We investigated the historical demography of Anopheles albimanus using mosquitoes from five countries and three different DNA regions, the mitochondrial cytochrome oxidase subunit I gene (COI), the single copy nuclear white gene and the ribosomal internal transcribed spacer two (ITS2). All the molecular markers supported the taxonomic status of a single species of An. albimanus. Furthermore, agreement between the COI and the white genes suggested a scenario of Pleistocene geographic fragmentation (i.e., population contraction) and subsequent range expansion across southern Central America.

Transposable elements in natural populations of Drosophila melanogaster

Transposable elements (TEs) are families of small DNA sequences found in the genomes of virtually all organisms. The sequences typically encode essential components for the replicative transposition sequences of that TE family. Thus, TEs are simply genomic parasites that inflict detrimental mutations on the fitness of their hosts. Several models have been proposed for the containment of TE copy number in outbreeding host populations such as Drosophila. Surveys of the TEs in genomes from natural populations of Drosophila have played a central role in the investigation of TE dynamics. The early surveys indicated that a typical TE insertion is rare in a population, which has been interpreted as evidence that each TE is selected against. The proposed mechanisms of this natural selection are reviewed here. Subsequent and more targeted surveys identify heterogeneity among types of TEs and also highlight the large role of homologous and possibly ectopic crossing over in the dynamics of the Drosophila TEs. The recent discovery of germline-specific RNA interference via the piwi-interacting RNA pathway opens yet another interesting mechanism that may be critical in containing the copy number of TEs in natural populations of Drosophila. The expected flood of Drosophila population genomics is expected to rapidly advance understanding of the dynamics of TEs.

Pronounced inter- and intrachromosomal variation in linkage disequilibrium across the zebra finch genome

The extent of nonrandom association of alleles at two or more loci, termed linkage disequilibrium (LD), can reveal much about population demography, selection, and recombination rate, and is a key consideration when designing association mapping studies. Here, we describe a genome-wide analysis of LD in the zebra finch (Taeniopygia guttata) using 838 single nucleotide polymorphisms and present LD maps for all assembled chromosomes. We found that LD declined with physical distance approximately five times faster on the microchromosomes compared to macrochromosomes. The distribution of LD across individual macrochromosomes also varied in a distinct pattern. In the center of the macrochromosomes there were large blocks of markers, sometimes spanning tens of mega bases, in strong LD whereas on the ends of macrochromosomes LD declined more rapidly. Regions of high LD were not simply the result of suppressed recombination around the centromere and this pattern has not been observed previously in other taxa. We also found evidence that this pattern of LD has remained stable across many generations. The variability in LD between and within chromosomes has important implications for genome wide association studies in birds and for our understanding of the distribution of recombination events and the processes that govern them.

Seasonal changes in the long-distance linkage disequilibrium in Drosophila melanogaster

Seasonal environmental changes have the potential to influence the genetic structure of species with a short generation time, such as Drosophila. We previously found the seasonal change in linkage disequilibrium (LD) between the chemoreceptor (Cr) genes in a local Japanese population (Kyoto [KY]). This could be caused by fluctuation in the population size or selection in temporally heterogeneous environments or both. Here, we analyzed the scale of LD between 51 X-linked polymorphisms (10 Cr and 41 non-Cr gene markers) in the 2 seasonal samples from the KY population and an autumn sample from 106 localities in and around Japan (Ja03au). Many of the non-Cr genes have receptor function but fewer functional connections to each other. The magnitude of LD in Ja03au did not significantly differ from that in the KY autumn sample. The lack of local differentiation was confirmed in an autumn sample from another local Japanese population. On the other hand, the magnitude of LD was significantly larger in spring than in autumn in the 2 independent KY samples. This suggests that reduction in the population size during winter increased the magnitude of LD in spring in the mainland population in Japan. Long-distance LD could be a useful measure for assessing seasonal fluctuation in effective population size.

Population analysis using the nuclear white gene detects Pliocene/Pleistocene lineage divergence within Anopheles nuneztovari in South America

Anopheles (Nyssorhynchus) nuneztovari Gabaldón (Diptera: Culicidae), a locally important malaria vector in some regions of South America, has been hypothesized to consist of at least two cryptic incipient species. We investigated its phylogeographic structure in several South American localities to determine the number of lineages and levels of divergence using the nuclear white gene, a marker that detected two recently diverged genotypes in the primary Neotropical malaria vector Anopheles darlingi Root. In An. nuneztovari, five distinct lineages (1-5) were elucidated: (1) populations from northeastern and central Amazonia; (2) populations from Venezuela east and west of the Andes; (3) populations from Colombia and Venezuela west of the Andes; (4) southeastern and western Amazonian Brazil populations, and (5) southeastern and western Amazonian Brazil and Bolivian populations. There was a large amount of genetic differentiation among these lineages. The deepest and earliest divergence was found between lineage 3 and lineages 1, 2 and 4, which probably accounts for the detection of lineage 3 in some earlier studies. The multiple lineages within Amazonia are partially congruent with previous mtDNA and ITS2 data, but were undetected in many earlier studies, probably because of their recent (Pleistocene) divergence and the differential mutation rates of the markers. The estimates for the five lineages, interpreted as recently evolved or incipient species, date to the Pleistocene and Pliocene. We hypothesize that the diversification in An. nuneztovari is the result of an interaction between the Miocene/Pliocene marine incursion and Pleistocene climatic changes leading to refugial isolation. The identification of cryptic lineages in An. nuneztovari could have a significant impact on local vector control measures.

Linkage disequilibrium--understanding the evolutionary past and mapping the medical future

Linkage disequilibrium--the nonrandom association of alleles at different loci--is a sensitive indicator of the population genetic forces that structure a genome. Because of the explosive growth of methods for assessing genetic variation at a fine scale, evolutionary biologists and human geneticists are increasingly exploiting linkage disequilibrium in order to understand past evolutionary and demographic events, to map genes that are associated with quantitative characters and inherited diseases, and to understand the joint evolution of linked sets of genes. This article introduces linkage disequilibrium, reviews the population genetic processes that affect it and describes some of its uses. At present, linkage disequilibrium is used much more extensively in the study of humans than in non-humans, but that is changing as technological advances make extensive genomic studies feasible in other species.

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.

Drosophila bristles and the nature of quantitative genetic variation

Numbers of Drosophila sensory bristles present an ideal model system to elucidate the genetic basis of variation for quantitative traits. Here, we review recent evidence that the genetic architecture of variation for bristle numbers is surprisingly complex. A substantial fraction of the Drosophila genome affects bristle number, indicating pervasive pleiotropy of genes that affect quantitative traits. Further, a large number of loci, often with sex- and environment-specific effects that are also conditional on background genotype, affect natural variation in bristle number. Despite this complexity, an understanding of the molecular basis of natural variation in bristle number is emerging from linkage disequilibrium mapping studies of individual candidate genes that affect the development of sensory bristles. We show that there is naturally segregating genetic variance for environmental plasticity of abdominal and sternopleural bristle number. For abdominal bristle number this variance can be attributed in part to an abnormal abdomen-like phenotype that resembles the phenotype of mutants defective in catecholamine biosynthesis. Dopa decarboxylase (Ddc) encodes the enzyme that catalyses the final step in the synthesis of dopamine, a major Drosophila catecholamine and neurotransmitter. We found that molecular polymorphisms at Ddc are indeed associated with variation in environmental plasticity of abdominal bristle number.

DNA sequence polymorphism and divergence at the erect wing and suppressor of sable loci of Drosophila melanogaster and D. simulans

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (approximately 50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high FST values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.

African sequence variation accounts for most of the sequence polymorphism in non-African Drosophila melanogaster

We compared the sequence polymorphism of 12 genomic fragments in six geographically dispersed African populations to one European Drosophila melanogaster population. On the basis of one African and one European population half of these fragments have strongly reduced levels of variability outside of Africa. Despite this striking difference in European variation, we detected no significant difference in African variation between the two fragment classes. The joint analysis of all African populations indicated that all high-frequency European alleles are of African origin. We observed a negative Tajima's D in all African populations, with three populations deviating significantly from neutral equilibrium. Low, but statistically significant, population differentiation was observed among the African populations. Our results imply that the population structure and demographic past of African D. melanogaster populations need to be considered for the inference of footprints of selection in non-African populations.

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.

Interlocus nonrandom association of polymorphisms in Drosophila chemoreceptor genes

Some forms of multilocus selection with epistasis, such as truncation selection, can effectively reduce the mutation load [Kondrashov, A. S. (1988) Nature 336, 435-440]. Many quantitative characters, including complex genetic diseases, are likely to be subject to these types of selection. However, direct measurement of selection in natural populations is difficult and the effect of epistasis on within-species variations remains unclear. Epistatic interaction in the fitness effect can generate linkage disequilibrium (LD). Therefore, we may detect the action of natural selection from its amount and pattern. Here, we report a large number of interlocus nonrandom associations between polymorphisms in 98 Drosophila chemoreceptor genes. LD was examined in two fly samples collected at the same location, but in different seasons. The amount of LD was much larger in the spring sample than in the autumn one. The between-sample difference was much more striking for the replacement polymorphisms than for the silent polymorphisms. This difference between the replacement and silent polymorphisms could not be attributed to differences in the mean marker distances. We also found a significant excess of associations between one frequent and one less common allele for the replacement polymorphisms, but not for the silent polymorphisms. It is unlikely that a simple seasonal bottleneck could explain all these differences in the scale of LD between the samples and between the replacement and silent polymorphisms. Natural selection is suggested to play a significant role in shaping the pattern of LD observed in this study.

World-wide survey of an Accord insertion and its association with DDT resistance in Drosophila melanogaster

Previous work showed that insecticide resistance in Drosophila melanogaster is correlated with the insertion of an Accord-like element into the 5' region of the cytochrome P450 gene, Cyp6g1. Here, we study the distribution of the Accord-like element in 673 recently collected D. melanogaster lines from 34 world-wide populations. We also examine the extent of microsatellite variability along a 180-kilobase (kb) genomic region of chromosome II encompassing the resistance gene. We confirm a 100% correlation of the Accord insertion with insecticide resistance and a significant reduction in variability extending at least 20 kb downstream of the Cyp6g1 gene. The frequency of the Accord insertion differs significantly between East African (32-55%) and nonAfrican (85-100%) populations. This pattern is consistent with a selective sweep driving the Accord insertion close to fixation in nonAfrican populations as a result of the insecticide resistance phenotype it confers. This study confirms that hitchhiking mapping can be used to identify beneficial mutations in natural populations.

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.

No association between mitochondrial DNA haplotypes and a female-limited mimicry phenotype in Papilio glaucus

Alternative alleles at a locus on the W chromosome of Papilio glaucus (causing dark or yellow wing colors, respectively) underlie a female-limited mimicry polymorphism thought to be maintained by balancing selection. In species with heterogametic females (i.e., the ZZ-male/ZW-female sex chromosome system), the mitochondrial DNA and the W chromosome are genetically linked because they are both maternally transmitted. We investigate the association of COI and COII mitochondrial DNA haplotypes with alternative W-linked phenotypes. Surprisingly, we find no congruence between mitochondrial DNA genealogies and inferred W-linked color alleles in P. glaucus. Using a maximum-likelihood phylogenetic approach, we reject the hypothesis of monophyly for dark-morph mitochondrial DNA lineages, even in the presence of putative low-frequency mimicry suppressor alleles or alternative melanizing factors. The most likely genealogical tree topologies assume more than one exchange event between mitochondrial DNA cytotype and the W-linked color morph. These results suggest that there is either paternal leakage of mitochondrial DNA or that more than two W-linked alleles underlie the alternative color morphs. Using data from an additional mitochondrial DNA locus, ND5, we show that pairwise linkage disequilibrium decays with physical distance between polymorphic sites. This finding suggests that genetic exchanges between maternal and paternal mitochondrial DNAs may have contributed to the lack of association we observe between phenotype and genotype.

hairy: A Quantitative Trait Locus for Drosophila Sensory Bristle Number

Advances in medicine, agriculture, and an understanding of evolution depend on resolving the genetic architecture of quantitative traits, which is challenging since variation for complex traits is caused by multiple interacting quantitative trait loci (QTL) with small and conditional effects. Here, we show that the key developmental gene, hairy (h), is a QTL for Drosophila sternopleural bristle number, a model quantitative trait. Near-isoallelic lines (NIL) for the h gene region exhibited significant variation in bristle number and failed to complement a hairy mutation. Sequencing 10 h alleles from a single population revealed 330 polymorphic sites in ∼10 kb. Genotypes for 25 of these and 14 additional sites in the flanking regions were determined for the 57 NIL and associated with variation in bristle number in four genetic backgrounds. A highly significant association was found for a complicated insertion/deletion polymorphism upstream of the transcription start site. This polymorphism, present in 17.5% of the h alleles, was associated with an increase of 0.5 bristle and accounted for 31% of the genetic variance in bristle number in the NIL.

A coalescent-based method for detecting and estimating recombination from gene sequences

Determining the amount of recombination in the genealogical history of a sample of genes is important to both evolutionary biology and medical population genetics. However, recurrent mutation can produce patterns of genetic diversity similar to those generated by recombination and can bias estimates of the population recombination rate. Hudson 2001 has suggested an approximate-likelihood method based on coalescent theory to estimate the population recombination rate, 4N(e)r, under an infinite-sites model of sequence evolution. Here we extend the method to the estimation of the recombination rate in genomes, such as those of many viruses and bacteria, where the rate of recurrent mutation is high. In addition, we develop a powerful permutation-based method for detecting recombination that is both more powerful than other permutation-based methods and robust to misspecification of the model of sequence evolution. We apply the method to sequence data from viruses, bacteria, and human mitochondrial DNA. The extremely high level of recombination detected in both HIV1 and HIV2 sequences demonstrates that recombination cannot be ignored in the analysis of viral population genetic data.

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.

The genetic architecture of quantitative traits

Phenotypic variation for quantitative traits results from the segregation of alleles at multiple quantitative trait loci (QTL) with effects that are sensitive to the genetic, sexual, and external environments. Major challenges for biology in the post-genome era are to map the molecular polymorphisms responsible for variation in medically, agriculturally, and evolutionarily important complex traits; and to determine their gene frequencies and their homozygous, heterozygous, epistatic, and pleiotropic effects in multiple environments. The ease with which QTL can be mapped to genomic intervals bounded by molecular markers belies the difficulty in matching the QTL to a genetic locus. The latter requires high-resolution recombination or linkage disequilibrium mapping to nominate putative candidate genes, followed by genetic and/or functional complementation and gene expression analyses. Complete genome sequences and improved technologies for polymorphism detection will greatly advance the genetic dissection of quantitative traits in model organisms, which will open avenues for exploration of homologous QTL in related taxa.

The nature of quantitative genetic variation for Drosophila longevity

Longevity is a typical quantitative trait: the continuous variation in life span observed in natural populations is attributable to genetic variation at multiple quantitative trait loci (QTL), environmental sensitivity of QTL alleles, and truly continuous environmental variation. To begin to understand the genetic architecture of longevity at the level of individual QTL, we have mapped QTL for Drosophila life span that segregate between two inbred strains that were not selected for longevity. A mapping population of 98 recombinant inbred lines (RIL) was derived from these strains, and life span of virgin male and female flies measured under control culture conditions, chronic heat and cold stress, heat shock and starvation stress, and high and low density larval environments. The genotypes of the RIL were determined for polymorphic roo transposable element insertion sites, and life span QTL were mapped using composite interval mapping methods. A minimum of 19 life span QTL were detected by recombination mapping. The life span QTL exhibited strong genotype by sex, genotype by environment, and genotype by genotype (epistatic) interactions. These interactions complicate mapping efforts, but evolutionary theory predicts such properties of segregating QTL alleles. Quantitative deficiency mapping of four longevity QTL detected in the control environment by recombination mapping revealed a minimum of 11 QTL in these regions. Clearly, longevity is a complex quantitative trait. In the future, linkage disequilibrium mapping can be used to determine which candidate genes in a QTL region correspond to the genetic loci affecting variation in life span, and define the QTL alleles at the molecular level.

Quantitative trait loci in Drosophila

Phenotypic variation for quantitative traits results from the simultaneous segregation of alleles at multiple quantitative trait loci. Understanding the genetic architecture of quantitative traits begins with mapping quantitative trait loci to broad genomic regions and ends with the molecular definition of quantitative trait loci alleles. This has been accomplished for some quantitative trait loci in Drosophila. Drosophila quantitative trait loci have sex-, environment- and genotype-specific effects, and are often associated with molecular polymorphisms in non-coding regions of candidate genes. These observations offer valuable lessons to those seeking to understand quantitative traits in other organisms, including humans.

Linkage disequilibria and the site frequency spectra in the su(s) and su(w(a)) regions of the Drosophila melanogaster X chromosome

Over the last decade, surveys of DNA sequence variation in natural populations of several Drosophila species and other taxa have established that polymorphism is reduced in genomic regions characterized by low rates of crossing over per physical length. Parallel studies have also established that divergence between species is not reduced in these same genomic regions, thus eliminating explanations that rely on a correlation between the rates of mutation and crossing over. Several theoretical models (directional hitchhiking, background selection, and random environment) have been proposed as population genetic explanations. In this study samples from an African population (n = 50) and a European population (n = 51) were surveyed at the su(s) (1955 bp) and su(w(a)) (3213 bp) loci for DNA sequence polymorphism, utilizing a stratified SSCP/DNA sequencing protocol. These loci are located near the telomere of the X chromosome, in a region of reduced crossing over per physical length, and exhibit a significant reduction in DNA sequence polymorphism. Unlike most previously surveyed, these loci reveal substantial skews toward rare site frequencies, consistent with the predictions of directional hitchhiking and random environment models and inconsistent with the general predictions of the background selection model (or neutral theory). No evidence for excess geographic differentiation at these loci is observed. Although linkage disequilibrium is observed between closely linked sites within these loci, many recombination events in the genealogy of the sampled alleles can be inferred and the genomic scale of linkage disequilibrium, measured in base pairs between sites, is the same as that observed for loci in regions of normal crossing over. We conclude that gene conversion must be high in these regions of low crossing over.

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.

Consequences of recombination on traditional phylogenetic analysis

We investigate the shape of a phylogenetic tree reconstructed from sequences evolving under the coalescent with recombination. The motivation is that evolutionary inferences are often made from phylogenetic trees reconstructed from population data even though recombination may well occur (mtDNA or viral sequences) or does occur (nuclear sequences). We investigate the size and direction of biases when a single tree is reconstructed ignoring recombination. Standard software (PHYLIP) was used to construct the best phylogenetic tree from sequences simulated under the coalescent with recombination. With recombination present, the length of terminal branches and the total branch length are larger, and the time to the most recent common ancestor smaller, than for a tree reconstructed from sequences evolving with no recombination. The effects are pronounced even for small levels of recombination that may not be immediately detectable in a data set. The phylogenies when recombination is present superficially resemble phylogenies for sequences from an exponentially growing population. However, exponential growth has a different effect on statistics such as Tajima's D. Furthermore, ignoring recombination leads to a large overestimation of the substitution rate heterogeneity and the loss of the molecular clock. These results are discussed in relation to viral and mtDNA data sets.

Both naturally occurring insertions of transposable elements and intermediate frequency polymorphisms at the achaete-scute complex are associated with variation in bristle number in Drosophila melanogaster

A restriction enzyme survey of a 110-kb region including the achaete scute complex (ASC) examined 14 polymorphic molecular markers in a sample of 56 naturally occurring chromosomes. Large insertions as a class were associated with a reduction in both sternopleural and abdominal bristle number, supporting deleterious mutation-selection equilibrium models for the maintenance of quantitative genetic variation. Two polymorphic sites were independently associated with variation in bristle number measured in two genetic backgrounds as assessed by a permutation test. A 6-bp deletion near sc alpha is associated with sternopleural bristle number variation in both sexes and a 3.4-kb insertion between sc beta and sc gamma is associated with abdominal bristle number variation in females. Under an additive genetic model, the small deletion polymorphism near sc alpha accounts for 25% of the total X chromosome genetic variation in sternopleural bristle number, and the 3.4 kb insertion accounts for 22% of the total X chromosome variation in female abdominal bristle number. The observation of common polymorphisms associated with variation in bristle number is more parsimoniously explained by models that incorporate balancing selection or assume variants affecting bristle number are neutral, than mutation-selection equilibrium models.

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.

Genetic variation in rates of nondisjunction: association of two naturally occurring polymorphisms in the chromokinesin nod with increased rates of nondisjunction in Drosophila melanogaster

Genetic variation in nondisjunction frequency among X chromosomes from two Drosophila melanogaster natural populations is examined in a sensitized assay. A high level of genetic variation is observed (a range of 0.006-0.241). Two naturally occurring variants at the nod locus, a chromokinesin required for proper achiasmate chromosome segregation, are significantly associated with an increased frequency of nondisjunction. Both of these polymorphisms are found at intermediate frequency in widely distributed natural populations. To account for these observations, we propose a general model incorporating unique opportunities for meiotic drive during female meiosis. The oötid competition model can account for both high mean rates of female-specific nondisjunction in Drosophila and humans as well as the standing genetic variation in this critical fitness character in natural populations.

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.

Two sites in the Delta gene region contribute to naturally occurring variation in bristle number in Drosophila melanogaster

A restriction enzyme survey of a 57-kb region including the gene Delta uncovered 53 polymorphic molecular markers in a sample of 55 naturally occurring chromosomes. A permutation test, which assesses the significance of the molecular marker with the largest effect on bristle variation in four genetic backgrounds relative to permuted data-sets, found two sites that were independently associated with variation in bristle number. A common site in the second intron of Delta affected only sternopleural bristle number, and another common site in the fifth intron affected only abdominal bristle number in females. Under an additive genetic model, the polymorphism in the second intron may account for 12% of the total genetic variation in sternopleural bristle number due to third chromosomes, and the site in the fifth intron may account for 6% of the total variation in female abdominal bristle number due to the third chromosomes. These results suggest the following: (1) models that incorporate balancing selection are more consistent with observations than deleterious mutation-selection equilibrium models, (2) mapped quantitative trait loci of large effect may not represent a single variable site at a genetic locus, and (3) linkage disequilibrium can be used as a tool for understanding the molecular basis of quantitative variation.