The Hill-Robertson effect: evolutionary consequences of weak selection and linkage in finite populations.

The 'Hill-Robertson (HR) effect' describes that linkage between sites under selection will reduce the overall effectiveness of selection in finite populations. Here we discuss the major concepts associated with the HR effect and present results of computer simulations focusing on the linkage effects generated by multiple sites under weak selection. Most models of linkage and selection forecast differences in effectiveness of selection between chromosomes or chromosomal regions involving a number of genes. The abundance and physical clustering of weakly selected mutations across genomes, however, justify the investigation of HR effects at a very local level and we pay particular attention to linkage effects among selected sites of the same gene. Overall, HR effects caused by weakly selected mutations predict differences in effectiveness of selection between genes that differ in exon-intron structures and across genes. Under this scenario, introns might play an advantageous role reducing intragenic HR effects. Finally, we summarize observations that are consistent with local HR effects in Drosophila, discuss potential consequences on population genetic studies and suggest future lines of research.

Differences in (G+C) content between species: a commentary on Forsdyke's "chromosomal viewpoint" of speciation.

Forsdyke (1999) has recently argued that differences in (G+C)%, or G+C content, may trigger new species formation. He further argues that the genic model has shortcomings that can be overcome by his "chromosomal" (hereafter, "G+C") model. We disagree on several counts. First, we do not accept that the genic model has the shortcomings suggested by Forsdyke. There is an abundance of empirical support for the contribution of individual genes, as well as of mapped chromosomal regions, to post-zygotic reproductive isolation (and Haldane's rule). Further, we argue that the G+C model suffers from the same theoretical difficulties as other speciation models based on underdominance. We also question the evidence Forsdyke uses to support his model. Finally, we describe analyses of G+C content in a well-studied model system of speciation (the Drosophila melanogaster species complex), the results of which are incompatible with the G+C model. Thus, while Forsdyke's G+C model cannot be explicitly ruled out, it is not directly supported by empirical data. In contrast, the genic model is well supported by empirical data, holds up on theoretical grounds, and does not require any assistance from the G+C model.

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 population genetics of the origin and divergence of the Drosophila simulans complex species.

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

Recent selection on synonymous codon usage in Drosophila.

Evidence from a variety of sources indicates that selection has influenced synonymous codon usage in Drosophila. It has generally been difficult, however, to distinguish selection that acted in the distant past from ongoing selection. However, under a neutral model, polymorphisms usually reflect more recent mutations than fixed differences between species and may, therefore, be useful for inferring recent selection. If the ancestral state is preferred, selection should shift the frequency distribution of derived states/site toward lower values; if the ancestral is unpreferred, selection should increase the number of derived states/site. Polymorphisms were classified as ancestrally preferred or unpreferred for several genes of D. simulans and D. melanogaster. A computer simulation of coalescence was employed to derive the expected frequency distributions of derived states/site under various modifications of the Wright-Fisher neutral model, and distributions of test statistics (t and Mann-Whitney U) were derived by appropriate sampling. One-tailed tests were applied to transformed frequency data to assess whether the two frequency distributions deviated from neutral expectations in the direction predicted by selection on codon usage. Several genes from D. simulans appear to be subject to recent selection on synonymous codons, including one gene with low codon bias, esterase-6. Selection may also be acting in D. melanogaster.

Mutation pressure, natural selection, and the evolution of base composition in Drosophila.

Genome sequencing in a number of taxa has revealed variation in nucleotide composition both among regions of the genome and among functional classes of sites in DNA. Mutational biases, biased gene conversion, and natural selection have been proposed as causes of this variation. Here, we review patterns of base composition in Drosophila DNA. Nucleotide composition in Drosophila melanogaster varys regionally, and base composition is correlated between introns and exons. Drosophila species also show striking patterns of non-random codon usage. Patterns of synonymous codon usage and the biochemistry of translation suggest that natural selection may act at 'silent' sites. A relationship between recombination rates and codon usage and comparisons of the evolutionary dynamics of silent mutations within and between species support natural selection discriminating among synonymous codons. The causes of regional base composition variation are less clear. Progress in functional studies of non-coding DNA, further investigations of genome patterns, and statistical tests based on evolutionary theory will lead to a greater understanding of the contributions of mutational processes and natural selection in patterning genome-wide nucleotide composition.

Patterns of base composition within the genes of Drosophila melanogaster.

Base composition is not uniform across the genome of Drosophila melanogaster. Earlier analyses have suggested that there is variation in composition in D. melanogaster on both a large scale and a much smaller, within-gene, scale. Here we present analyses on 117 genes which have reliable intron/exon boundaries and no known alternative splicing. We detect significant heterogeneity in G+C content among intron segments from the same gene, as well as a significant positive correlation between the intron and the third codon position G+C content within genes. Both of these observations appear to be due, in part, to an overall decline in intron and third codon position G+C content along Drosophila genes with introns. However, there is also evidence of an increase in third codon position G+C content at the start of genes; this is particularly evident in genes without introns. This is consistent with selection acting against preferred codons at the start of genes.

The effects of mutation and natural selection on codon bias in the genes of Drosophila.

Codon bias varies widely among the loci of Drosophila melanogaster, and some of this diversity has been explained by variation in the strength of natural selection. A study of correlations between intron and coding region base composition shows that variation in mutation pattern also contributes to codon bias variation. This finding is corroborated by an analysis of variance (ANOVA), which shows a tendency for introns from the same gene to be similar in base composition. The strength of base composition correlations between introns and codon third positions is greater for genes with low codon bias than for genes with high codon bias. This pattern can be explained by an overwhelming effect of natural selection, relative to mutation, in highly biased loci. In particular, this correlation is absent when examining fourfold degenerate sites of highly biased genes. In general, it appears that selection acts more strongly in choosing among fourfold degenerate codons than among twofold degenerate codons. Although the results indicate regional variation in mutational bias, no evidence is found for large scale regions of compositional homogeneity.

Reduced natural selection associated with low recombination in Drosophila melanogaster.

Synonymous codons are not used equally in many organisms, and the extent of codon bias varies among loci. Earlier studies have suggested that more highly expressed loci in Drosophila melanogaster are more biased, consistent with findings from several prokaryotes and unicellular eukaryotes that codon bias is partly due to natural selection for translational efficiency. We link this model of varying selection intensity to the population-genetics prediction that the effectiveness of natural selection is decreased under reduced recombination. In analyses of 385 D. melanogaster loci, we find that codon bias is reduced in regions of low recombination (i.e., near centromeres and telomeres and on the fourth chromosome). The effect does not appear to be a linear function of recombination rate; rather, it seems limited to regions with the very lowest levels of recombination. The large majority of the genome apparently experiences recombination at a sufficiently high rate for effective natural selection against suboptimal codons. These findings support models of the Hill-Robertson effect and genetic hitchhiking and are largely consistent with multiple reports of low levels of DNA sequence variation in regions of low recombination.

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex.

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1-kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58-0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.

DNA sequence variation at the period locus within and among species of the Drosophila melanogaster complex.

A 1.9-kilobase region of the period locus was sequenced in six individuals of Drosophila melanogaster and from six individuals of each of three sibling species: Drosophila simulans, Drosophila sechellia and Drosophila mauritiana. Extensive genealogical analysis of 174 polymorphic sites reveals a complex history. It appears that D. simulans, as a large population still segregating very old lineages, gave rise to the island species D. mauritiana and D. sechellia. Rather than considering these speciation events as having produced "sister" taxa, it seems more appropriate to consider D. simulans a parent species to D. sechellia and D. mauritiana. The order, in time, of these two phylogenetic events remains unclear. D. mauritiana supports a large number of polymorphisms, many of which are shared with D. simulans, and so appears to have begun and persisted as a large population. In contrast, D. sechellia has very little variation and seems to have experienced a severe population bottleneck. Alternatively, the low variation in D. sechellia could be due to recent directional selection and genetic hitchhiking at or near the per locus.

Evidence for genetic variation in the occurrence of the photoresponse of the Djungarian hamster, Phodopus sungorus.

The Djungarian hamster generally responds to a short-day photoperiod with a complex syndrome of physiological and behavioral changes; however, not all hamsters are photoresponsive. The phenotypic difference is, in part, genetically determined. Parent-offspring regression on a number of continuous and discontinuous measures indicated significant heritability for photoresponsiveness. Four generations of replicated bidirectional selection on a photoresponse index (PI) resulted in significant shifts in the percentage of responsive hamsters, although the average PI of responsive individuals was not significantly changed. Eight estimates of heritability ranged from 0.20 to 0.52. We hypothesize that the circadian system is responsible for the occurrence of the photoresponse, but that the extent of photoresponse is controlled by a separate functional system.

Evidence for independence of circadian characters and extent of photoresponsiveness in the Djungarian hamster, Phodopus sungorus.

Djungarian hamsters (Phodopus sungorus) exposed to a short-day photoperiod generally respond with a syndrome of physiological and behavioral changes, such as body weight loss and molt to a white pelage. The extent of the short-day-induced responses differs among individuals. Furthermore, some hamsters show no photoresponse. In this study, we sought to determine whether variation in the photoresponse would be associated with circadian function: whether phase angle or free-running period (tau) would differ between responsive and nonresponsive hamsters; and whether changes in these circadian characters would correlate with the extent of weight loss and molt (and the timing of molt onset) in photoresponsive hamsters. Adult hamsters were kept in a short-day photoperiod (9 hr light, 15 hr dark) for 14 weeks, during which time body weight and molt were measured biweekly. Hamsters were then transferred to cages equipped with running wheels; we measured the phase angle of activity onset under a short-day photoperiod and tau in constant dark. Hamsters exhibiting a short-day-induced molt had a significantly shorter tau and a less negative phase angle than nonmolting animals. Hamsters that exhibited weight loss also had a significantly less negative phase angle, but no difference in tau. No significant Pearson's or Spearman's correlation coefficients were found between extent (or timing) of the photoresponse and the circadian characters in responsive hamsters. Although these results indicate that threshold for photoresponsiveness is related to circadian function, the extent (and timing) of the photoresponse may not be.

Genetic analyses of photoresponsiveness in the Djungarian hamster, Phodopus sungorus.

Endotherms living at temperate and arctic latitudes must adjust their physiology and behavior in order to survive seasonal change. The Djungarian hamster uses photoperiod to cue annual cycles of reproduction and thermoregulation, and its responses to short photoperiod include loss of body weight and change in pelage color. Some individuals do not exhibit these responses when exposed to short days. In this study individual variation in photoresponsiveness is quantified, and four lines of evidence for a genetic component to that variation are provided. First, two separate breeding stocks differed in both the percent of animals responding to a short-day lighting regimen (SD) and in the degree and timing of their response. Second, analysis of variance within and between families of full sibs for a photoresponsive index, PI (body weight loss +2 (molt index -1] following 12 weeks in SD demonstrated a significant family resemblance (intraclass correlation of 0.36 +/- 0.03). Third, heritability estimates from regression of offspring scores on parent scores for body weight loss, molt index and PI after 12 weeks in SD were 0.34 +/- 0.13, 0.36 +/- 0.10 and 0.37 +/- 0.12, respectively, indicating a strong additive genetic component for the three characters. Finally, a significant response occurred after one generation of artificial selection for and against photoresponsiveness.