DNA sequence variation and the recombinational landscape in Drosophila pseudoobscura: a study of the second chromosome

Positive Selection and Functional Divergence at Meiosis Genes That Mediate Crossing Over Across the Drosophila Phylogeny

Meiotic crossing over ensures proper segregation of homologous chromosomes and generates genotypic diversity. Despite these functions, little is known about the genetic factors and population genetic forces involved in the evolution of recombination rate differences among species. The dicistronic meiosis gene, mei-217/mei-218, mediates most of the species differences in crossover rate and patterning during female meiosis between the closely related fruitfly species, Drosophila melanogaster and D. mauritiana The MEI-218 protein is one of several meiosis-specific mini-chromosome maintenance (mei-MCM) proteins that form a multi-protein complex essential to crossover formation, whereas the BLM helicase acts as an anti-crossover protein. Here we study the molecular evolution of five genes- mei-218, the other three known members of the mei-MCM complex, and Blm- over the phylogenies of three Drosophila species groups- melanogaster, obscura, and virilis We then use transgenic assays in D. melanogaster to test if molecular evolution at mei-218 has functional consequences for crossing over using alleles from the distantly related species D. pseudoobscura and D. virilis Our molecular evolutionary analyses reveal recurrent positive selection at two mei-MCM genes. Our transgenic assays show that sequence divergence among mei-218 alleles from D. melanogaster, D. pseudoobscura, and D. virilis has functional consequences for crossing over. In a D. melanogaster genetic background, the D. pseudoobscura mei-218 allele nearly rescues wildtype crossover rates but alters crossover patterning, whereas the D. virilis mei-218 allele conversely rescues wildtype crossover patterning but not crossover rates. These experiments demonstrate functional divergence at mei-218 and suggest that crossover rate and patterning are separable functions.

Recent selection on the Y-to-dot translocation in Drosophila pseudoobscura

The Drosophila pseudoobscura dot chromosome acquired genes from the ancestral Drosophila Y chromosome in a Y-to-dot translocation event that occurred between 12.7 and 20.8 Ma. The formerly Y-linked genes mostly retained their testis-specific expression but shrank drastically in size, mostly through intron reduction, since becoming part of the dot chromosome in this species. We investigated the impact of this translocation on the evolution of the both the Y-to-dot translocated region and the original segments of the dot chromosome in D. pseudoobscura. Our survey of polymorphism and divergence across the chromosome reveals a reduction in variation, a deletion polymorphism segregating at high frequency, and a shift in the frequency spectra, all consistent with a history of recent selective sweeps in the Y-to-dot translocated region but not on the rest of the dot chromosome. We do find evidence for recombination primarily as gene conversion on the dot chromosome; however, predicted recombination events are restricted to the part of the dot chromosome outside the translocation. It therefore appears that recombination has resulted in a degree of decoupling between the ancestral Y region and the conserved region of the dot chromosome.

Molecular population genetics of inversion breakpoint regions in Drosophila pseudoobscura

Paracentric inversions in populations can have a profound effect on the pattern and organization of nucleotide variability along a chromosome. Regions near inversion breakpoints are expected to have greater levels of differentiation because of reduced genetic exchange between different gene arrangements whereas central regions in the inverted segments are predicted to have lower levels of nucleotide differentiation due to greater levels of genetic flux among different karyotypes. We used the inversion polymorphism on the third chromosome of Drosophila pseudoobscura to test these predictions with an analysis of nucleotide diversity of 18 genetic markers near and away from inversion breakpoints. We tested hypotheses about how the presence of different chromosomal arrangements affects the pattern and organization of nucleotide variation. Overall, markers in the distal segment of the chromosome had greater levels of nucleotide heterozygosity than markers within the proximal segment of the chromosome. In addition, our results rejected the hypothesis that the breakpoints of derived inversions will have lower levels of nucleotide variability than breakpoints of ancestral inversions, even when strains with gene conversion events were removed. High levels of linkage disequilibrium were observed within all 11 breakpoint regions as well as between the ends of most proximal and distal breakpoints. The central region of the chromosome had the greatest levels of linkage disequilibrium compared with the proximal and distal regions because this is the region that experiences the highest level of recombination suppression. These data do not fully support the idea that genetic exchange is the sole force that influences genetic variation on inverted chromosomes.

Fine-scale heterogeneity in crossover rate in the garnet-scalloped region of the Drosophila melanogaster X chromosome

Homologous recombination affects myriad aspects of genome evolution, from standing levels of nucleotide diversity to the efficacy of natural selection. Rates of crossing over show marked variability at all scales surveyed, including species-, population-, and individual-level differences. Even within genomes, crossovers are nonrandomly distributed in a wide diversity of taxa. Although intra- and intergenomic heterogeneities in crossover distribution have been documented in Drosophila, the scale and degree of crossover rate heterogeneity remain unclear. In addition, the genetic features mediating this heterogeneity are unknown. Here we quantify fine-scale heterogeneity in crossover distribution in a 2.1-Mb region of the Drosophila melanogaster X chromosome by localizing crossover breakpoints in 2500 individuals, each containing a single crossover in this specific X chromosome region. We show 90-fold variation in rates of crossing over at a 5-kb scale, place this variation in the context of several aspects of genome evolution, and identify several genetic features associated with crossover rates. Our results shed new light on the scale and magnitude of crossover rate heterogeneity in D. melanogaster and highlight potential features mediating this heterogeneity.

The effects of deleterious mutations on evolution at linked sites

The process of evolution at a given site in the genome can be influenced by the action of selection at other sites, especially when these are closely linked to it. Such selection reduces the effective population size experienced by the site in question (the Hill-Robertson effect), reducing the level of variability and the efficacy of selection. In particular, deleterious variants are continually being produced by mutation and then eliminated by selection at sites throughout the genome. The resulting reduction in variability at linked neutral or nearly neutral sites can be predicted from the theory of background selection, which assumes that deleterious mutations have such large effects that their behavior in the population is effectively deterministic. More weakly selected mutations can accumulate by Muller's ratchet after a shutdown of recombination, as in an evolving Y chromosome. Many functionally significant sites are probably so weakly selected that Hill-Robertson interference undermines the effective strength of selection upon them, when recombination is rare or absent. This leads to large departures from deterministic equilibrium and smaller effects on linked neutral sites than under background selection or Muller's ratchet. Evidence is discussed that is consistent with the action of these processes in shaping genome-wide patterns of variation and evolution.

Odorant receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura lineages

BACKGROUND

In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ∼37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster--OR65c and OR67a--and D. pseudoobscura--OR65b1 and OR67c.

CONCLUSIONS

The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.

A recombination survey using microsatellites: the O chromosome of Drosophila subobscura

Recombination plays an important role in species adaptation since it acts as an evolutionary force that can influence genome pattern organization. However, recombination can be detrimental in some situations, causing the breakdown of some adaptive gene combinations such as coadapted gene complexes. Genetic and cytological chromosome maps allow recombination throughout the genome to be analyzed. In this study we compare the recombination rate of two types of homokaryotypic lines of D. subobscura (O(ST) and O( 3+4 )) using a set of at least 13 microsatellite loci. The genetic maps obtained present similar lengths: 184 and 196 cM for O(ST) and O( 3+4 ) chromosomes, respectively. For most pairs of markers analyzed, a sample size of about 150 individuals appeared sufficient to obtain appropriate recombination values, with the exception of markers located in the same cytological band. Recombination rates seemed to be fairly uniform along the O chromosome, but some regional differences were observed. Several recombination hot and coldspots were detected, and their numbers were different in the homokaryotypic line types (O(ST) and O( 3+4 )). This variability could be attributed to differences between the genetic content of the two arrangements or to differences between the lines.

Advances in maize genomics and their value for enhancing genetic gains from breeding

Maize is an important crop for food, feed, forage, and fuel across tropical and temperate areas of the world. Diversity studies at genetic, molecular, and functional levels have revealed that, tropical maize germplasm, landraces, and wild relatives harbor a significantly wider range of genetic variation. Among all types of markers, SNP markers are increasingly the marker-of-choice for all genomics applications in maize breeding. Genetic mapping has been developed through conventional linkage mapping and more recently through linkage disequilibrium-based association analyses. Maize genome sequencing, initially focused on gene-rich regions, now aims for the availability of complete genome sequence. Conventional insertion mutation-based cloning has been complemented recently by EST- and map-based cloning. Transgenics and nutritional genomics are rapidly advancing fields targeting important agronomic traits including pest resistance and grain quality. Substantial advances have been made in methodologies for genomics-assisted breeding, enhancing progress in yield as well as abiotic and biotic stress resistances. Various genomic databases and informatics tools have been developed, among which MaizeGDB is the most developed and widely used by the maize research community. In the future, more emphasis should be given to the development of tools and strategic germplasm resources for more effective molecular breeding of tropical maize products.

Estimation of fine-scale recombination intensity variation in the white-echinus interval of D. melanogaster

Accurate assessment of local recombination rate variation is crucial for understanding the recombination process and for determining the impact of natural selection on linked sites. In Drosophila, local recombination intensity has been estimated primarily by statistical approaches, by estimating the local slope of the relationship between the physical and genetic maps. However, these estimates are limited in resolution and, as a result, the physical scale at which recombination intensity varies in Drosophila is largely unknown. Although there is some evidence suggesting as much as a 40-fold variation in crossover rate at a local scale in D. pseudoobscura, little is known about the fine-scale structure of recombination rate variation in D. melanogaster. Here we experimentally examine the fine-scale distribution of crossover events in a 1.2-Mb region on the D. melanogaster X chromosome using a classic genetic mapping approach. Our results show that crossover frequency is significantly heterogeneous within this region, varying approximately 3.5-fold. Simulations suggest that this degree of heterogeneity is sufficient to affect levels of standing nucleotide diversity, although the magnitude of this effect is small. We recover no statistical association between empirical estimates of nucleotide diversity and recombination intensity, which is likely due to the limited number of loci sampled in our population genetic data set. However, codon bias is significantly negatively correlated with fine-scale recombination intensity estimates, as expected. Our results shed light on the relevant physical scale to consider in evolutionary analyses relating to recombination rate and highlight the motivations to increase the resolution of the recombination map in Drosophila.

Nucleotide polymorphism and within-gene recombination in Daphnia magna and D. pulex, two cyclical parthenogens

Theory predicts that partially asexual organisms may make the "best of both worlds": for the most part, they avoid the costs of sexual reproduction, while still benefiting from an enhanced efficiency of selection compared to obligately asexual organisms. There is, however, little empirical data on partially asexual organisms to test this prediction. Here we examine patterns of nucleotide diversity at eight nuclear loci in continentwide samples of two species of cyclically parthenogenetic Daphnia to assess the effect of partial asexual reproduction on effective population size and amount of recombination. Both species have high nucleotide diversities and show abundant evidence for recombination, yielding large estimates of effective population sizes (300,000-600,000). This suggests that selection will act efficiently even on mutations with small selection coefficients. Divergence between the two species is less than one-tenth of previous estimates, which were derived using a mitochondrial molecular clock. As the two species investigated are among the most distantly related species of the genus, this suggests that the genus Daphnia may be considerably younger than previously thought. Daphnia has recently received increased attention because it is being developed as a model organism for ecological and evolutionary genomics. Our results confirm the attractiveness of Daphnia as a model organism, because the high nucleotide diversity and low linkage disequilibrium suggest that fine-scale mapping of genes affecting phenotypes through association studies should be feasible.

Fine-scale mapping of recombination rate in Drosophila refines its correlation to diversity and divergence

Regional rates of recombination often correlate with levels of nucleotide diversity, and either selective or neutral hypotheses can explain this relationship. Regional recombination rates also correlate with nucleotide differences between human and chimpanzee, consistent with models where recombination is mutagenic; however, a lack of correlation is observed in the Drosophila melanogaster group, consistent with models invoking natural selection. Here, we revisit the relationship among recombination, diversity, and interspecies difference by generating empirical estimates of these parameters in Drosophila pseudoobscura. To measure recombination rate, we genotyped 1,294 backcross hybrids at 50 markers across the largest assembled linkage group in this species. Genome-wide diversity was estimated by sequencing a second isolate of D. pseudoobscura at shallow coverage. Alignment to the sequenced genome of the closely related species, Drosophila persimilis, provided nucleotide site orthology. Our findings demonstrate that scale is critical in determining correlates to recombination rate: fine-scale cross-over rate estimates are far stronger predictors of both diversity and interspecies difference than broad-scale estimates. The correlation of fine-scale recombination rate to diversity and interspecies difference appears to be genome-wide, evidenced by examination of an X-linked region in greater detail. Because we observe a strong correlation of cross-over rate with interspecies difference, even after correcting for segregating ancestral variation, we suggest that both mutagenic and selective forces generate these correlations, the latter in regions of low crossing over. We propose that it is not cross-overs per se that are mutagenic, but rather repair of DNA double-strand break precursors via crossing over and gene conversion.

Divergence between the Drosophila pseudoobscura and D. persimilis genome sequences in relation to chromosomal inversions

As whole-genome sequence assemblies accumulate, a challenge is to determine how these can be used to address fundamental evolutionary questions, such as inferring the process of speciation. Here, we use the sequence assemblies of Drosophila pseudoobscura and D. persimilis to test hypotheses regarding divergence with gene flow. We observe low differentiation between the two genome sequences in pericentromeric and peritelomeric regions. We interpret this result as primarily a remnant of the correlation between levels of variation and local recombination rate observed within populations. However, we also observe lower differentiation far from the fixed chromosomal inversions distinguishing these species and greater differentiation within and near these inversions. This finding is consistent with models suggesting that chromosomal inversions facilitate species divergence despite interspecies gene flow. We also document heterogeneity among the inverted regions in their degree of differentiation, suggesting temporal differences in the origin of each inverted region consistent with the inversions arising during a process of divergence with gene flow. While this study provides insights into the speciation process using two single-genome sequences, it was informed by lower throughput but more rigorous examinations of polymorphism and divergence. This reliance highlights the need for complementary genomic and population genetic approaches for tackling fundamental evolutionary questions such as speciation.

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.

Genome scanning for interspecific differentiation between two closely related oak species [Quercus robur L. and Q. petraea (Matt.) Liebl.]

Interspecific differentiation values (G(ST)) between two closely related oak species (Quercus petraea and Q. robur) were compiled across different studies with the aim to explore the distribution of differentiation at the genome level. The study was based on a total set of 389 markers (isozymes, AFLPs, SCARs, microsatellites, and SNPs) for which allelic frequencies were estimated in pairs of populations sampled throughout the sympatric distribution of the two species. The overall distribution of G(ST) values followed an L-shaped curve with most markers exhibiting low species differentiation (G(ST) < 0.01) and only a few loci reaching >10% levels. Twelve percent of the loci exhibited significant G(ST) deviations to neutral expectations, suggesting that selection contributed to species divergence. Coding regions expressed higher differentiation than noncoding regions. Among the 389 markers, 158 could be mapped on the 12 linkage groups of the existing Q. robur genetic map. Outlier loci with large G(ST) values were distributed over 9 linkage groups. One cluster of three outlier loci was found within 0.51 cM; but significant autocorrelation of G(ST) was observed at distances <2 cM. The size and distribution of genomic regions involved in species divergence are discussed in reference to hitchhiking effects and disruptive selection.

Evolutionarily stable infection by a male-killing endosymbiont in Drosophila innubila: molecular evidence from the host and parasite genomes

Maternally inherited microbes that spread via male-killing are common pathogens of insects, yet very little is known about the evolutionary duration of these associations. The few examples to date indicate very recent, and thus potentially transient, infections. A male-killing strain of Wolbachia has recently been discovered in natural populations of Drosophila innubila. The population-level effects of this infection are significant: approximately 35% of females are infected, infected females produce very strongly female-biased sex ratios, and the resulting population-level sex ratio is significantly female biased. Using data on infection prevalence and Wolbachia transmission rates, infected cytoplasmic lineages are estimated to experience a approximately 5% selective advantage relative to uninfected lineages. The evolutionary history of this infection was explored by surveying patterns of polymorphism in both the host and parasite genomes, comparing the Wolbachia wsp gene and the host mtDNA COI gene to five host nuclear genes. Molecular data suggest that this male-killing infection is evolutionarily old, a conclusion supported with a simple model of parasite and mtDNA transmission dynamics. Despite a large effective population size of the host species and strong selection to evolve resistance, the D. innubila-Wolbachia association is likely at a stable equilibrium that is maintained by imperfect maternal transmission of the bacteria rather than partial resistance in the host species.

Patterns of selection on synonymous and nonsynonymous variants in Drosophila miranda

We have investigated patterns of within-species polymorphism and between-species divergence for synonymous and nonsynonymous variants at a set of autosomal and X-linked loci of Drosophila miranda. D. pseudoobscura and D. affinis were used for the between-species comparisons. The results suggest the action of purifying selection on nonsynonymous, polymorphic variants. Among synonymous polymorphisms, there is a significant excess of synonymous mutations from preferred to unpreferred codons and of GC to AT mutations. There was no excess of GC to AT mutations among polymorphisms at noncoding sites. This suggests that selection is acting to maintain the use of preferred codons. Indirect evidence suggests that biased gene conversion in favor of GC base pairs may also be operating. The joint intensity of selection and biased gene conversion, in terms of the product of effective population size and the sum of the selection and conversion coefficients, was estimated to be approximately 0.65.

The incomplete natural history of mitochondria

Mitochondrial DNA (mtDNA) has been used to study molecular ecology and phylogeography for 25 years. Much important information has been gained in this way, but it is time to reflect on the biology of the mitochondrion itself and consider opportunities for evolutionary studies of the organelle itself and its ecology, biochemistry and physiology. This review has four sections. First, we review aspects of the natural history of mitochondria and their DNA to show that it is a unique molecule with specific characteristics that differ from nuclear DNA. We do not attempt to cover the plethora of differences between mitochondrial and nuclear DNA; rather we spotlight differences that can cause significant bias when inferring demographic properties of populations and/or the evolutionary history of species. We focus on recombination, effective population size and mutation rate. Second, we explore some of the difficulties in interpreting phylogeographical data from mtDNA data alone and suggest a broader use of multiple nuclear markers. We argue that mtDNA is not a sufficient marker for phylogeographical studies if the focus of the investigation is the species and not the organelle. We focus on the potential bias caused by introgression. Third, we show that it is not safe to assume a priori that mtDNA evolves as a strictly neutral marker because both direct and indirect selection influence mitochondria. We outline some of the statistical tests of neutrality that can, and should, be applied to mtDNA sequence data prior to making any global statements concerning the history of the organism. We conclude with a critical examination of the neglected biology of mitochondria and point out several surprising gaps in the state of our knowledge about this important organelle. Here we limelight mitochondrial ecology, sexually antagonistic selection, life-history evolution including ageing and disease, and the evolution of mitochondrial inheritance.

A High-Density Genetic Recombination Map of Sequence-Tagged Sites for Sorghum, as a Framework for Comparative Structural and Evolutionary Genomics of Tropical Grains and Grasses

We report a genetic recombination map for Sorghum of 2512 loci spaced at average 0.4 cM (∼300 kb) intervals based on 2050 RFLP probes, including 865 heterologous probes that foster comparative genomics of Saccharum (sugarcane), Zea (maize), Oryza (rice), Pennisetum (millet, buffelgrass), the Triticeae (wheat, barley, oat, rye), and Arabidopsis. Mapped loci identify 61.5% of the recombination events in this progeny set and reveal strong positive crossover interference acting across intervals of ≤50 cM. Significant variations in DNA marker density are related to possible centromeric regions and to probable chromosome structural rearrangements between Sorghum bicolor and S. propinquum, but not to variation in levels of intraspecific allelic richness. While cDNA and genomic clones are similarly distributed across the genome, SSR-containing clones show different abundance patterns. Rapidly evolving hypomethylated DNA may contribute to intraspecific genomic differentiation. Nonrandom distribution patterns of multiple loci detected by 357 probes suggest ancient chromosomal duplication followed by extensive rearrangement and gene loss. Exemplifying the value of these data for comparative genomics, we support and extend prior findings regarding maize-sorghum synteny—in particular, 45% of comparative loci fall outside the inferred colinear/syntenic regions, suggesting that many small rearrangements have occurred since maize-sorghum divergence. These genetically anchored sequence-tagged sites will foster many structural, functional and evolutionary genomic studies in major food, feed, and biomass crops.

A survey of chromosomal and nucleotide sequence variation in Drosophila miranda

There have recently been several studies of the evolution of Y chromosome degeneration and dosage compensation using the neo-sex chromosomes of Drosophila miranda as a model system. To understand these evolutionary processes more fully, it is necessary to document the general pattern of genetic variation in this species. Here we report a survey of chromosomal variation, as well as polymorphism and divergence data, for 12 nuclear genes of D. miranda. These genes exhibit varying levels of DNA sequence polymorphism. Compared to its well-studied sibling species D. pseudoobscura, D. miranda has much less nucleotide sequence variation, and the effective population size of this species is inferred to be several-fold lower. Nevertheless, it harbors a few inversion polymorphisms, one of which involves the neo-X chromosome. There is no convincing evidence for a recent population expansion in D. miranda, in contrast to D. pseudoobscura. The pattern of population subdivision previously observed for the X-linked gene period is not seen for the other loci, suggesting that there is no general population subdivision in D. miranda. However, data on an additional region of period confirm population subdivision for this gene, suggesting that local selection is operating at or near period to promote differentiation between populations.

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.

SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines

BACKGROUND

Recent studies of ancestral maize populations indicate that linkage disequilibrium tends to dissipate rapidly, sometimes within 100 bp. We set out to examine the linkage disequilibrium and diversity in maize elite inbred lines, which have been subject to population bottlenecks and intense selection by breeders. Such population events are expected to increase the amount of linkage disequilibrium, but reduce diversity. The results of this study will inform the design of genetic association studies.

RESULTS

We examined the frequency and distribution of DNA polymorphisms at 18 maize genes in 36 maize inbreds, chosen to represent most of the genetic diversity in U.S. elite maize breeding pool. The frequency of nucleotide changes is high, on average one polymorphism per 31 bp in non-coding regions and 1 polymorphism per 124 bp in coding regions. Insertions and deletions are frequent in non-coding regions (1 per 85 bp), but rare in coding regions. A small number (2-8) of distinct and highly diverse haplotypes can be distinguished at all loci examined. Within genes, SNP loci comprising the haplotypes are in linkage disequilibrium with each other.

CONCLUSIONS

No decline of linkage disequilibrium within a few hundred base pairs was found in the elite maize germplasm. This finding, as well as the small number of haplotypes, relative to neutral expectation, is consistent with the effects of breeding-induced bottlenecks and selection on the elite germplasm pool. The genetic distance between haplotypes is large, indicative of an ancient gene pool and of possible interspecific hybridization events in maize ancestry.

What has QTL mapping taught us about plant domestication?

The aim of this paper is to survey the general area of quantitative trait locus (QTL) mapping, and its specific impact on current understanding of plant domestication. Plant domestication is not only of historical interest, but is also of ongoing importance as changing human needs and availability of nonrenewable resources impel continuing (and perhaps even accelerated) investigation of prospective new crops. New genomic tools applied in conjunction with now-established approaches such as QTL mapping are opening new doors into searches for the 'footprints' of domestication, and promise to accelerate and streamline the identification of specific genes integral to domestication(s), building on early successes. Better understanding of plant domestication promises to enhance knowledge about the developmental basis of some of the more striking evolutionary events known, to guide efforts to catalog plant biodiversity, and to accelerate progress in improving existing and new crops to sustain humanity. Contents Summary 591 I. Introduction 592 II. A backdrop: QTL mapping basics 593 III. The tempo of domestication 596 IV. Domestication and polyploidy 601 V. New approaches to identifying the footprints of domestication 603 VI. Perspectives 605 Acknowledgements 606 References 606.

Population, evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.

Inferring the history of speciation from multilocus DNA sequence data: the case of Drosophila pseudoobscura and close relatives

The divergence of Drosophila pseudoobscura from its close relatives, D. persimilis and D. pseudoobscura bogotana, was examined using the pattern of DNA sequence variation in a common set of 50 inbred lines at 11 loci from diverse locations in the genome. Drosophila pseudoobscura and D. persimilis show a marked excess of low-frequency variation across loci, consistent with a model of recent population expansion in both species. The different loci vary considerably, both in polymorphism levels and in the levels of polymorphisms that are shared by different species pairs. A major question we address is whether these patterns of shared variation are best explained by gene flow or by persistence since common ancestry. A new test of gene flow, based on patterns of linkage disequilibrium, is developed. The results from these, and other tests, support a model in which D. pseudoobscura and D. persimilis have exchanged genes at some loci. However, the pattern of variation suggests that most gene flow, although occurring after speciation began, was not recent. There is less evidence of gene flow between D. pseudoobscura and D. p. bogotana. The results are compared with recent work on the genomic locations of genes that contribute to reproductive isolation between D. pseudoobscura and D. persimilis. We show that there is a good correspondence between the genomic regions associated with reproductive isolation and the regions that show little or no evidence of gene flow.

Patterns of genetic variation in Mendelian and complex traits

This review discusses the prospects for understanding the genetic basis of complex traits in humans. We take the view that work done on Drosophila melanogaster can serve as a model for understanding complex traits in humans, and the literature on this model system, as well as on humans, is reviewed. The prospects for success in understanding the genetic basis of complex traits depend, in part, on the nature of the forces acting on genetic variation. We suggest that different experimental approaches should be undertaken for traits caused by common genetic variants versus those arising from rare genetic variants.

Consequences of recombination rate variation on quantitative trait locus mapping studies. Simulations based on the Drosophila melanogaster genome

We examine the effect of variation in gene density per centimorgan on quantitative trait locus (QTL) mapping studies using data from the Drosophila melanogaster genome project and documented regional rates of recombination. There is tremendous variation in gene density per centimorgan across this genome, and we observe that this variation can cause systematic biases in QTL mapping studies. Specifically, in our simulated mapping experiments of 50 equal-effect QTL distributed randomly across the physical genome, very strong QTL are consistently detected near the centromeres of the two major autosomes, and few or no QTL are often detected on the X chromosome. This pattern persisted with varying heritability, marker density, QTL effect sizes, and transgressive segregation. Our results are consistent with empirical data collected from QTL mapping studies of this species and its close relatives, and they explain the "small X-effect" that has been documented in genetic studies of sexual isolation in the D. melanogaster group. Because of the biases resulting from recombination rate variation, results of QTL mapping studies should be taken as hypotheses to be tested by additional genetic methods, particularly in species for which detailed genetic and physical genome maps are not available.

Divergent structure and function of the bicoid gene in Muscoidea fly species

We have investigated the evolution of the bicoid (bcd) gene in fly species of the Muscoidea Superfamily. We obtained the complete bcd sequence from the housefly Musca domestica and found polymorphism in the coding region among Musca strains. In addition to Musca, we cloned most of the bcd coding sequences from two blowfly species Calliphora vicina and Lucilia sericata. The 5' and 3' regulatory regions flanking the Musca bcd gene are widely diverged in sequence from Drosophila; however, some important sequence motifs identified in Drosophila bcd are present. The predicted RNA secondary structures of the 3' UTRs are similar, despite sequence divergence. Comparison of Bicoid (Bcd) proteins shows a serine-rich domain of unknown function is present in the Muscoidea species, but is absent in other species. The in vivo function of bcd in Musca was tested by RNAi to mimic loss of function phenotype. We obtained a head defect phenotype similar to weak bcd alleles of Drosophila. Although our comparisons initially suggest functional conservation between species, closer inspection reveals significant differences. Divergence of structural motifs, such as regulatory elements in flanking regions and conservation of protein domains in some species but not in others, points to functional divergence between species. We suggest that the larger embryonic size in Muscoidea species restricts the morphogenetic activity of a weak Bcd activator, which has evolved a more specialized role in head determination and lost some functions in thoracic development.

Evolutionary history of microsatellites in the obscura group of Drosophila

The evolutionary origins of microsatellites are not well understood. Some investigators have suggested that point mutations that expand repeat arrays beyond a threshold size trigger microsatellites to become variable. However, little empirical data has been brought forth on this and related issues. In this study, we examine the evolutionary history of microsatellites in six species within the obscura group of Drosophila, tracing changes in microsatellite alleles using both PCR product size and sequence data. We found little evidence supporting a general role of point mutations triggering initial microsatellite expansion, and no consistent threshold size for expansion was observed. Flanking region length variation was extensive when alleles were sequenced in distantly related species, and some species possessed altogether different repeat arrays between the same primer binding sites. Our results suggest extreme caution in using microsatellite allele sizes for phylogenetic analyses or to infer divergences between populations.

The genetics of reproductive isolation and the potential for gene exchange between Drosophila pseudoobscura and D. persimilis via backcross hybrid males

Hybrid male sterility, hybrid inviability, sexual isolation, and a hybrid male courtship dysfunction reproductively isolate Drosophila pseudoobscura and D. persimilis. Previous studies of the genetic bases of these isolating mechanisms have yielded only limited information about how much and what areas of the genome are susceptible to interspecies introgression. We have examined the genetic basis of these barriers to gene exchange in several thousand backcross hybrid male progeny of these species using 14 codominant molecular genetic markers spanning the five chromosomes of these species, focusing particularly on the autosomes. Hybrid male sterility, hybrid inviability, and the hybrid male courtship dysfunction were all associated with X-autosome interactions involving primarily the inverted regions on the left arm of the X-chromosome and the center of the second chromosome. Sexual isolation from D. pseudoobscura females was primarily associated with the left arm of the X-chromosome, although both the right arm and the center of the second chromosome also contributed to it. Sexual isolation from D. persimilis females was primarily associated with the second chromosome. The absence of isolating mechanisms being associated with many autosomal regions, including some large inverted regions that separate the strains, suggests that these phenotypes may not be caused by genes spread throughout the genome. We suggest that gene flow between these species via hybrid males may be possible at loci spread across much of the autosomes.

Toward integration of comparative genetic, physical, diversity, and cytomolecular maps for grasses and grains, using the sorghum genome as a foundation

The small genome of sorghum (Sorghum bicolor L. Moench.) provides an important template for study of closely related large-genome crops such as maize (Zea mays) and sugarcane (Saccharum spp.), and is a logical complement to distantly related rice (Oryza sativa) as a "grass genome model." Using a high-density RFLP map as a framework, a robust physical map of sorghum is being assembled by integrating hybridization and fingerprint data with comparative data from related taxa such as rice and using new methods to resolve genomic duplications into locus-specific groups. By taking advantage of allelic variation revealed by heterologous probes, the positions of corresponding loci on the wheat (Triticum aestivum), rice, maize, sugarcane, and Arabidopsis genomes are being interpolated on the sorghum physical map. Bacterial artificial chromosomes for the small genome of rice are shown to close several gaps in the sorghum contigs; the emerging rice physical map and assembled sequence will further accelerate progress. An important motivation for developing genomic tools is to relate molecular level variation to phenotypic diversity. "Diversity maps," which depict the levels and patterns of variation in different gene pools, shed light on relationships of allelic diversity with chromosome organization, and suggest possible locations of genomic regions that are under selection due to major gene effects (some of which may be revealed by quantitative trait locus mapping). Both physical maps and diversity maps suggest interesting features that may be integrally related to the chromosomal context of DNA-progress in cytology promises to provide a means to elucidate such relationships. We seek to provide a detailed picture of the structure, function, and evolution of the genome of sorghum and its relatives, together with molecular tools such as locus-specific sequence-tagged site DNA markers and bacterial artificial chromosome contigs that will have enduring value for many aspects of genome analysis.

Population genomics: genome-wide sampling of insect populations

Modern population genetics underwent a major paradigm shift during the last decade of the 20th century with the discovery that thousands of genes of known function and position in a genome can be analyzed simultaneously in a single individual. The impact of this technology on insect population genetics is potentially profound. Sampling distributions of genetic statistics can now be derived from many individual loci or among many segregating sites within a gene. Inferences regarding random mating, gene flow, effective population sizes, disequilibrium, and relatedness among populations can now be based on patterns of variation at many loci. More importantly, genome-wide sampling enables population geneticists to distinguish effects that act on the whole genome from those that act on individual loci or nucleotides. We introduce the term "population genomics" to describe the process of simultaneous sampling of numerous variable loci within a genome and the inference of locus-specific effects from the sample distributions. The four critical assumptions implicit in the population genomics approach are explained in detail. Studies adopting this paradigm are reviewed, and the steps necessary to complete a population genomics study are outlined.

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.

Contrasting patterns of nonneutral evolution in proteins encoded in nuclear and mitochondrial genomes

We report that patterns of nonneutral DNA sequence evolution among published nuclear and mitochondrially encoded protein-coding loci differ significantly in animals. Whereas an apparent excess of amino acid polymorphism is seen in most (25/31) mitochondrial genes, this pattern is seen in fewer than half (15/36) of the nuclear data sets. This differentiation is even greater among data sets with significant departures from neutrality (14/15 vs. 1/6). Using forward simulations, we examined patterns of nonneutral evolution using parameters chosen to mimic the differences between mitochondrial and nuclear genetics (we varied recombination rate, population size, mutation rate, selective dominance, and intensity of germ line bottleneck). Patterns of evolution were correlated only with effective population size and strength of selection, and no single genetic factor explains the empirical contrast in patterns. We further report that in Arabidopsis thaliana, a highly self-fertilizing plant with effectively low recombination, five of six published nuclear data sets also exhibit an excess of amino acid polymorphism. We suggest that the contrast between nuclear and mitochondrial nonneutrality in animals stems from differences in rates of recombination in conjunction with a distribution of selective effects. If the majority of mutations segregating in populations are deleterious, high linkage may hinder the spread of the occasional beneficial mutation.

Maize as a model for the evolution of plant nuclear genomes

The maize genome is replete with chromosomal duplications and repetitive DNA. The duplications resulted from an ancient polyploid event that occurred over 11 million years ago. Based on DNA sequence data, the polyploid event occurred after the divergence between sorghum and maize, and hence the polyploid event explains some of the difference in DNA content between these two species. Genomic rearrangement and diploidization followed the polyploid event. Most of the repetitive DNA in the maize genome is retrotransposable elements, and they comprise 50% of the genome. Retrotransposon multiplication has been relatively recent-within the last 5-6 million years-suggesting that the proliferation of retrotransposons has also contributed to differences in DNA content between sorghum and maize. There are still unanswered questions about repetitive DNA, including the distribution of repetitive DNA throughout the genome, the relative impacts of retrotransposons and chromosomal duplication in plant genome evolution, and the hypothesized correlation of duplication events with transposition. Population genetic processes also affect the evolution of genomes. We discuss how centromeric genes should, in theory, contain less genetic diversity than noncentromeric genes. In addition, studies of diversity in the wild relatives of maize indicate that different genes have different histories and also show that domestication and intensive breeding have had heterogeneous effects on genetic diversity across genes.

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.