Rate variation of DNA sequence evolution in the Drosophila lineages

Genome-Wide Patterns of Sequence Divergence of Protein-Coding Genes Between Drosophila buzzatii and D. mojavensis

Evolutionary rates for protein-coding genes are determined not only by natural selection but also by multiple genomic factors including mutation rates, recombination, gene expression levels, and chromosomal location. To investigate the joint effects of different genomic determinants on protein evolution, we compared the coding sequences of 9017 single-copy orthologs between 2 cactophilic species from the Drosophila subgenus, Drosophila mojavensis and D. buzzatii, whose genomes have been previously sequenced. We assessed the impact of 7 genomic determinants, that is, chromosome type, recombination, chromosomal inversions, expression breadth, expression level, gene length, and the number of exons, on divergence rates of protein-coding genes to understand patterns of evolutionary variation. Integrative analysis of these factors revealed that 1) X-linked and autosomal genes evolve at significantly different rates in agreement with the faster-X hypothesis, 2) genes located on the dot chromosome and pericentromeric regions have higher divergence rates, 3) genes located at chromosomes with more fixed inversions have higher pairwise divergence than those located at nearly collinear chromosomes, and 4) gene expression patterns can be considered the strongest determinant of protein evolution. In addition, the number of exons and protein length had a significant effect on pairwise divergence at synonymous sites. All in all, our results show the relative importance of each genomic factor on the rates of protein evolution and functional constraint in these 2 cactophilic Drosophila species.

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.

Maximum-likelihood model averaging to profile clustering of site types across discrete linear sequences

A major analytical challenge in computational biology is the detection and description of clusters of specified site types, such as polymorphic or substituted sites within DNA or protein sequences. Progress has been stymied by a lack of suitable methods to detect clusters and to estimate the extent of clustering in discrete linear sequences, particularly when there is no a priori specification of cluster size or cluster count. Here we derive and demonstrate a maximum likelihood method of hierarchical clustering. Our method incorporates a tripartite divide-and-conquer strategy that models sequence heterogeneity, delineates clusters, and yields a profile of the level of clustering associated with each site. The clustering model may be evaluated via model selection using the Akaike Information Criterion, the corrected Akaike Information Criterion, and the Bayesian Information Criterion. Furthermore, model averaging using weighted model likelihoods may be applied to incorporate model uncertainty into the profile of heterogeneity across sites. We evaluated our method by examining its performance on a number of simulated datasets as well as on empirical polymorphism data from diverse natural alleles of the Drosophila alcohol dehydrogenase gene. Our method yielded greater power for the detection of clustered sites across a breadth of parameter ranges, and achieved better accuracy and precision of estimation of clusters, than did the existing empirical cumulative distribution function statistics.

Strong evidence for lineage and sequence specificity of substitution rates and patterns in Drosophila

Rates of single nucleotide substitution in Drosophila are highly variable within the genome, and several examples illustrate that evolutionary rates differ among Drosophila species as well. Here, we use a maximum likelihood method to quantify lineage-specific substitutional patterns and apply this method to 4-fold degenerate synonymous sites and introns from more than 8,000 genes aligned in the Drosophila melanogaster group. We find that within species, different classes of sequence evolve at different rates, with long introns evolving most slowly and short introns evolving most rapidly. Relative rates of individual single nucleotide substitutions vary approximately 3-fold among lineages, yielding patterns of substitution that are comparatively less GC-biased in the melanogaster species complex relative to Drosophila yakuba and Drosophila erecta. These results are consistent with a model coupling a mutational shift toward reduced GC content, or a shift in mutation-selection balance, in the D. melanogaster species complex, with variation in selective constraint among different classes of DNA sequence. Finally, base composition of coding and intronic sequences is not at equilibrium with respect to substitutional patterns, which primarily reflects the slow rate of the substitutional process. These results thus support the view that mutational and/or selective processes are labile on an evolutionary timescale and that if the process is indeed selection driven, then the distribution of selective constraint is variable across the genome.

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.

Strong regional heterogeneity in base composition evolution on the Drosophila X chromosome

Fluctuations in base composition appear to be prevalent in Drosophila and mammal genome evolution, but their timescale, genomic breadth, and causes remain obscure. Here, we study base composition evolution within the X chromosomes of Drosophila melanogaster and five of its close relatives. Substitutions were inferred on six extant and two ancestral lineages for 14 near-telomeric and 9 nontelomeric genes. GC content evolution is highly variable both within the genome and within the phylogenetic tree. In the lineages leading to D. yakuba and D. orena, GC content at silent sites has increased rapidly near telomeres, but has decreased in more proximal (nontelomeric) regions. D. orena shows a 17-fold excess of GC-increasing vs. AT-increasing synonymous changes within a small (approximately 130-kb) region close to the telomeric end. Base composition changes within introns are consistent with changes in mutation patterns, but stronger GC elevation at synonymous sites suggests contributions of natural selection or biased gene conversion. The Drosophila yakuba lineage shows a less extreme elevation of GC content distributed over a wider genetic region (approximately 1.2 Mb). A lack of change in GC content for most introns within this region suggests a role of natural selection in localized base composition fluctuations.

X chromosomes and autosomes evolve at similar rates in Drosophila: no evidence for faster-X protein evolution

Recent data from Drosophila suggest that a substantial fraction of amino acid substitutions observed between species are beneficial. If these beneficial mutations are on average partially recessive, then the rate of protein evolution is predicted to be faster for X-linked genes compared to autosomal genes (the "faster-X" hypothesis). We test this prediction by comparing rates of protein substitutions between orthologous genes, taking advantage of variations in chromosome fusions within the genus Drosophila. In members of the Drosophila melanogaster species group, the chromosomal arm 3L segregates as an ordinary autosome (i.e., two homologous copies in both males and females). However, in the Drosophila pseudoobscura species group, this chromosomal arm has become fused to the ancestral X chromosome and is hemizygous in males. The faster-X hypothesis predicts that protein evolution should be faster for genes on this chromosomal arm in the D. pseudoobscura lineage, relative to the D. melanogaster lineage. Here we combine new sequence data for 202 gene fragments in Drosophila miranda (in the pseudoobscura species group) with the completed genomes of D. melanogaster, D. pseudoobscura, and Drosophila yakuba to show that there are no detectable differences in rates of amino acid evolution for orthologous X-linked and autosomal genes. Our results imply that the contribution of the faster-X (if any) to the large-X effect on reproductive isolation in Drosophila is not due to a generally faster rate of protein evolution. The lack of a detectable faster-X effect in these species suggests either that beneficial amino acids are not partially recessive on average, or that adaptive evolution does not often use newly arising amino acid mutations.

Genomic tools and cDNA derived markers for butterflies

The Lepidoptera have long been used as examples in the study of evolution, but some questions remain difficult to resolve due to a lack of molecular genetic data. However, as technology improves, genomic tools are becoming increasingly available to tackle unanswered evolutionary questions. Here we have used expressed sequence tags (ESTs) to develop genetic markers for two Müllerian mimic species, Heliconius melpomene and Heliconius erato. In total 1363 ESTs were generated, representing 330 gene objects in H. melpomene and 431 in H. erato. User-friendly bioinformatic tools were used to construct a nonredundant database of these putative genes (available at http://www.heliconius.org), and annotate them with blast similarity searches, InterPro matches and Gene Ontology terms. This database will be continually updated with EST sequences for the Papilionideae as they become publicly available, providing a tool for gene finding in the butterflies. Alignments of the Heliconius sequences with putative homologues derived from Bombyx mori or other public data sets were used to identify conserved PCR priming sites, and develop 55 markers that can be amplified from genomic DNA in both H. erato and H. melpomene. These markers will be used for comparative linkage mapping in Heliconius and will have applications in other phylogenetic and genomic studies in the Lepidoptera.

Substitution patterns in alleles of immunoglobulin V genes in humans and mice

Immunoglobulins (Igs) constitute a subfamily of rapidly evolving proteins. It is postulated that this characteristic is due mainly to the participation of these proteins in highly diverse functions of recognition and defense. Although this vision of rapid evolution in Igs is widely accepted, various studies have demonstrated that diverse and contradictory forces not yet completely understood converge in the evolution of these receptors. In a recent study of the substitution patterns in the alleles that form the human IGHV locus, we found that the variation in genetic and structural information does not occur homogeneously among the different genes, nor among the regions and positions conforming said locus. In view of these results and of the importance of a better understanding of the basic evolutionary process in specific receptors (such as Igs) for both immunology and molecular evolution, it is important to explore the nature of the diversification process in these proteins in detail. In this work, therefore, we analyzed the substitution patterns in all the alleles reported for loci IGKV and IGLV in humans and mice, and we compared the results with those previously observed in the human IGHV locus. We found that the process of evolutionary variation of the Igs reflect the diversity of selective pressures operating on the different loci, genes, sub-regions and positions; for example, diversification through substitution is generally centered on CDRs, but only few positions inside the CDRs were frequently substituted. In spite of this general tendency, it is possible to observe differences in the degree of diversification among loci, families and genes. These tendencies to modify only certain attributes of IGV genes seem to be in agreement with differential strategies associated with the restrictions of the molecular immune recognition mechanism. The complexity of the evolutionary patterns observed in this study leads us to think that the predispositions observed herein may also be due in part to processes of DNA dynamics.

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.

Molecular clock-like evolution of human immunodeficiency virus type 1

The molecular clock hypothesis states that the rate of nucleotide substitution per generation is constant across lineages. If generation times were equal across lineages, samples obtained at the same calendar time would have experienced the same number of generations since their common ancestor. However, if sequences are not derived from contemporaneous samples, differences in the number of generations may be misinterpreted as variation in substitution rates and hence may lead to false rejection of the molecular clock hypothesis. A recent study has called into doubt the validity of clock-like evolution for HIV-1, using molecular sequences derived from noncontemporaneous samples. However, after separating their within-individual data according to sampling time, we found that what appeared to be nonclock-like behavior could be attributed, in most cases, to noncontemporaneous sampling, with contributions also likely to derive from recombination. Natural selection alone did not appear to obscure the clock-like evolution of HIV-1.

Patterns of Polymorphism and Divergence from Noncoding Sequences of Drosophila melanogaster and D. simulans: Evidence for Nonequilibrium Processes

Despite the fact that D. melanogaster and D. simulans have been the central model system for molecular population genetics, few data are available for noncoding regions. Here, we present an analysis of population genetic data from intergenic regions and comparisons of these data to previously collected data from introns and exons. Polymorphisms and fixations were categorized as A/T to G/C or G/C to A/T changes and were polarized by inferring the ancestral state using both parsimony and maximum likelihood. Noncoding fixations in both D. melanogaster and D. simulans were consistent with equilibrium base-composition evolution. However, polarized noncoding polymorphisms, revealed a different pattern. Although A/T to G/C and G/C to A/T polymorphisms in D. simulans were consistent with equilibrium, we observed a highly significant dearth of A/T to G/C polymorphisms in D. melanogaster introns but not in intergenic sequences. Such data could be explained by recent evolution of mutational biases associated with transcription or by lineage-specific selection on base composition. These data reveal the complexity of evolutionary processes acting even on noncoding DNA in Drosophila.

Complex evolution of orthologous and paralogous decarboxylase genes

The decarboxylases are involved in neurotransmitter synthesis in animals, and in pathways of secondary metabolism in plants. Different decarboxylase proteins are characterized for their different substrate specificities, but are encoded by homologous genes. We study, within a maximum-likelihood framework, the evolutionary relationships among dopa decarboxylase (Ddc), histidine decarboxylase (Hdc) and alpha-methyldopa hypersensitive (amd) in animals, and tryptophan decarboxylase (Wdc) and tyrosine decarboxylase (Ydc) in plants. The evolutionary rates are heterogeneous. There are differences between paralogous genes in the same lineages: 4.13 x 10(-10) nucleotide substitutions per site per year in mammalian Ddc vs. 1.95 in Hdc; between orthologous genes in different lineages, 7.62 in dipteran Ddc vs. 4.13 in mammalian Ddc; and very large temporal variations in some lineages, from 3.7 up to 54.9 in the Drosophila Ddc lineage. Our results are inconsistent with the molecular clock hypothesis.

Sequence divergence, functional constraint, and selection in protein evolution

The genome sequences of multiple species has enabled functional inferences from comparative genomics. A primary objective is to infer biological functions from the conservation of homologous DNA sequences between species. A second, more difficult, objective is to understand what functional DNA sequences have changed over time and are responsible for species' phenotypic differences. The neutral theory of molecular evolution provides a theoretical framework in which both objectives can be explicitly tested. Development of statistical tests within this framework has provided insight into the evolutionary forces that constrain and in some cases change DNA sequences and the resulting patterns that emerge. In this article, we review recent work on how functional constraint and changes in protein function are inferred from protein polymorphism and divergence data. We relate these studies to our understanding of the neutral theory and adaptive evolution.

Recruitment of the Proneural Gene scute to the Drosophila Sex-Determination Pathway

In flies, scute (sc) works with its paralogs in the achaete-scute-complex (ASC) to direct neuronal development. However, in the family Drosophilidae, sc also acquired a role in the primary event of sex determination, X chromosome counting, by becoming an X chromosome signal element (XSE)—an evolutionary step shown here to have occurred after sc diverged from its closest paralog, achaete (ac). Two temperature-sensitive alleles, scsisB2 and scsisB3, which disrupt only sex determination, were recovered in a powerful F1 genetic selection and used to investigate how sc was recruited to the sex-determination pathway. scsisB2 revealed 3′ nontranscribed regulatory sequences likely to be involved. The scsisB2 lesion abolished XSE activity when combined with mutations engineered in a sequence upstream of all XSEs. In contrast, changes in Sc protein sequence seem not to have been important for recruitment. The observation that the other new allele, scsisB3, eliminates the C-terminal half of Sc without affecting neurogenesis and that scsisB1, the most XSE-specific allele previously available, is a nonsense mutant, would seem to suggest the opposite, but we show that housefly Sc can substitute for fruit fly Sc in sex determination, despite lacking Drosophilidae-specific conserved residues in its C-terminal half. Lack of synergistic lethality among mutations in sc, twist, and dorsal argue against a proposed role for sc in mesoderm formation that had seemed potentially relevant to sex-pathway recruitment. The screen that yielded new sc alleles also generated autosomal duplications that argue against the textbook view that fruit fly sex signal evolution recruited a set of autosomal signal elements comparable to the XSEs.

Natural selection drives Drosophila immune system evolution

Evidence from disparate sources suggests that natural selection may often play a role in the evolution of host immune system proteins. However, there have been few attempts to make general population genetic inferences on the basis of analysis of several immune-system-related genes from a single species. Here we present DNA polymorphism and divergence data from 34 genes thought to function in the innate immune system of Drosophila simulans and compare these data to those from 28 nonimmunity genes sequenced from the same lines. Several statistics, including average K(A)/K(S) ratio, average silent heterozygosity, and average haplotype diversity, significantly differ between the immunity and nonimmunity genes, suggesting an important role for directional selection in immune system protein evolution. In contrast to data from mammalian immunoglobulins and other proteins, we find no strong evidence for the selective maintenance of protein diversity in Drosophila immune system proteins. This may be a consequence of Drosophila's generalized innate immune response.

Molecular population genetics of Xdh and the evolution of base composition in Drosophila

Few loci have been measured for DNA polymorphism and divergence in several species. Here we report such data from the protein-coding region of xanthine dehydrogenase (Xdh) in 22 species of Drosophila. Many of our samples were from closely related species, allowing us to confidently assign substitutions to individual lineages. Surprisingly, Xdh appears to be fixing more A/T mutations than G/C mutations in most lineages, leading to evolution of higher A/T content in the recent past. We found no compelling evidence for selection on protein variation, though some aspects of the data support the notion that a significant fraction of amino acid polymorphisms are slightly deleterious. Finally, we found no convincing evidence that levels of silent heterozygosity are associated with rates of protein evolution.

scuteexpression inCalliphora vicinareveals an ancestral pattern of longitudinal stripes on the thorax of higher Diptera

In Drosophila the stereotyped arrangement of sensory bristles on the notum is determined by the tightly regulated control of transcription of the achaete-scute (ac-sc) genes which are expressed in small proneural clusters of cells at the sites of each future bristle. Expression relies on a series of discrete cis-regulatory elements present in the ac-sc gene complex that are the target of the transcriptional activators pannier (pnr) and the genes of the iroquois complex. Stereotyped bristle patterns are common among species of acalyptrate Schizophora such as Drosophila, and are thought to have derived from an ancestral pattern of four longitudinal rows extending the length of the scutum, through secondary loss of bristles. To investigate evolutionary changes in bristle patterns and ac-sc regulation by pnr, we have isolated homologues of these genes from Calliphora vicina, a species of calyptrate Schizophora separated from Drosophila by at least 100 million years. Calliphora vicina displays a pattern of four rows of bristles on the scutum resembling the postulated ancestral one. We find that sc in Calliphora is expressed in two longitudinal stripes on the medial scutum that prefigure the development of the rows of acrostichal and dorsocentral bristles. This result suggests that a stripe-like expression pattern of sc may be an ancestral feature and may have preceded the evolution of proneural clusters. The implications for the evolution of the cis-regulatory elements responsible for sc expression in the proneural clusters of Drosophila, and function of Pnr are discussed.

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.

Neurogenesis in the spiderCupiennius salei

To uncover similarities and differences in neurogenesis in arthropod groups, we have studied the ventral neuroectoderm of the spider Cupiennius salei (Chelicerata, Aranea, Ctenidae). We found that invaginating cell groups arose sequentially, at stereotyped positions in each hemisegment and in separate waves, comparable with the generation of neuroblasts in Drosophila. However, we found no evidence for proliferating stem cells that would be comparable with the neuroblasts. Instead, the whole group of invaginating cells was directly recruited to the nervous system. The invagination process is comparable with Drosophila, with the cells attaining a bottle-shaped form with the nuclei moving inwards, while actin-rich cell processes remain initially connected to the surface of the epithelium. This general pattern is also found in another spider, Pholcus phalangioides, and appears thus to be conserved at least among the Araneae. We have identified two basic helix-loop-helix encoding genes – CsASH1 and CsASH2 – that share sequence similarities with proneural genes from other species. Functional analysis of the genes by double-stranded RNA interference revealed that CsASH1 was required for the formation of the invagination sites and the process of invagination itself, whereas CsASH2 seemed to be required for the differentiation of the cells into neurones. Our results suggest that the basic processes of neurogenesis, as well as proneural gene function is conserved among arthropods, apart of the lack of neuroblast-like stem cells in spiders.

Local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes

I present here evidence of remarkable local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes. The substitution pattern at 10 loci in the telomeric region of the X chromosome was studied for four species of the Drosophila melanogaster species subgroup. Drosophila orena and Drosophila erecta are clearly the most closely related species pair (the erecta complex) among the four species studied; however, the overall data at the 10 loci revealed a clear dichotomy in the silent substitution patterns between the AT-biased- substitution melanogaster and erecta lineages and the GC-biased-substitution yakuba and orena lineages, suggesting two or more independent changes in GC/AT substitution biases. More importantly, the results indicated a between- loci heterogeneity in GC/AT substitution bias in this small region independently in the yakuba and orena lineages. Indeed, silent substitutions in the orena lineage were significantly biased toward G and C at the consecutive yellow, lethal of scute, and asense loci, but they were significantly biased toward A and T at sta. The substitution bias toward G and C was centered in different areas in yakuba (significantly biased at EG:165H7.3, EG:171D11.2, and suppressor of sable). The similar silent substitution patterns in coding and noncoding regions, furthermore, suggested mutational biases as a cause of the substitution biases. On the other hand, previous study reveals that Drosophila yakuba has about 20-fold higher crossover frequencies in the telomeric region of the X chromosome than does D. melanogaster; this study revealed that the total genetic map length of the yakuba X chromosome was only about 1.5 times as large as that of melanogaster and that the map length of the X-telomeric y-sta region did not differ between Drosophila yakuba and D. erecta. Taken together, the data strongly suggested that an approximately 20- fold reduction in the X-telomeric crossover frequencies occurred in the ancestral population of D. melanogaster after the melanogaster-yakuba divergence but before the melanogaster-simulans divergence.

The correlation between intron length and recombination in drosophila. Dynamic equilibrium between mutational and selective forces

Intron length is negatively correlated with recombination in both Drosophila melanogaster and humans. This correlation is not likely to be the result of mutational processes alone: evolutionary analysis of intron length polymorphism in D. melanogaster reveals equivalent ratios of deletion to insertion in regions of high and low recombination. The polymorphism data do reveal, however, an excess of deletions relative to insertions (i.e., a deletion bias), with an overall deletion-to-insertion events ratio of 1.35. We propose two types of selection favoring longer intron lengths. First, the natural mutational bias toward deletion must be opposed by strong selection in very short introns to maintain the minimum intron length needed for the intron splicing reaction. Second, selection will favor insertions in introns that increase recombination between mutations under the influence of selection in adjacent exons. Mutations that increase recombination, even slightly, will be selectively favored because they reduce interference among selected mutations. Interference selection acting on intron length mutations must be very weak, as indicated by frequency spectrum analysis of Drosophila intron length polymorphism, making the equilibrium for intron length sensitive to changes in the recombinational environment and population size. One consequence of this sensitivity is that the advantage of longer introns is expected to decrease inversely with the rate of recombination, thus leading to a negative correlation between intron length and recombination rate. Also in accord with this model, intron length differs between closely related Drosophila species, with the longest variant present more often in D. melanogaster than in D. simulans. We suggest that the study of the proposed dynamic model, taking into account interference among selected sites, might shed light on many aspects of the comparative biology of genome sizes including the C value paradox.

Genetic analysis of bristle loss in hybrids between Drosophila melanogaster and D. simulans provides evidence for divergence of cis-regulatory sequences in the achaete-scute gene complex

The two closely related species of Drosophila, D. melanogaster and D. simulans, display an identical bristle pattern on the notum, but hybrids between the two are lacking a variable number of bristles. We show that the loss is temperature-dependent and provide evidence for two periods of temperature sensitivity. A first period of heat sensitivity occurs during larval development and corresponds to the time when the prepattern of expression of genes whose products activate achaete-scute in the proneural clusters preceding bristle precursor formation is established. A second period of cold sensitivity corresponds to the time of emergence of the bristle precursor cells and the maintenance of their neural fate, a process requiring high levels of Achaete-Scute. Expression of achaete-scute at these two critical periods depends on cis-regulatory elements of the achaete-scute complex (AS-C). The differences between males, which have only one copy of the X-linked AS-C from D. simulans, and females, which have copies from both parental species, are compared, together with the effects of crossing in different rearrangements of the D. melanogaster AS-C that delete regulatory and/or coding sequences. We provide evidence that bristle loss in the hybrids may result from a decrease in the level of transcription at the AS-C and argue that interaction between trans-acting factors and cis-regulatory elements within the AS-C has diverged between the two species.

Adaptive evolution of relish, a Drosophila NF-kappaB/IkappaB protein

NF-kappaB and IkappaB proteins have central roles in regulation of inflammation and innate immunity in mammals. Homologues of these proteins also play an important role in regulation of the Drosophila immune response. Here we present a molecular population genetic analysis of Relish, a Drosophila NF-kappaB/IkappaB protein, in Drosophila simulans and D. melanogaster. We find strong evidence for adaptive protein evolution in D. simulans, but not in D. melanogaster. The adaptive evolution appears to be restricted to the IkappaB domain. A possible explanation for these results is that Relish is a site of evolutionary conflict between flies and their microbial pathogens.

Local recombination and mutation effects on molecular evolution in Drosophila

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.

Unusual haplotype structure at the proximal breakpoint of In(2L)t in a natural population of Drosophila melanogaster

The existence of temporally stable frequency clines for In(2L)t in natural populations of Drosophila melanogaster suggests a role for selection in the maintenance of this polymorphism. We have collected nucleotide polymorphism data from the proximal breakpoint junction regions of In(2L)t to infer its evolutionary history. The finding of a novel LINE-like element near the In(2L)t breakpoint junction in sampled inverted chromosomes supports a transposable element-mediated origin for this inversion. An analysis of nucleotide variation in a Costa Rican population sample of standard and inverted chromosomes indicates a unique and relatively recent origin for In(2L)t. Additional In(2L)t alleles from three geographically diverse populations reveal no detectable geographic differentiation. Low levels of In(2L)t nucleotide polymorphism suggest a recent increase in the inversion's frequency in tropical populations. An unusual feature of our sample of standard alleles is a marked heterogeneity in levels of linkage disequilibrium among polymorphic sites across the breakpoint region. We introduce a test of neutral equilibrium haplotype structure that corrects both for multiple tests and for an arbitrarily chosen window size. It reveals that an approximately 1.4-kb region immediately spanning the breakpoint has fewer haplotypes than expected under the neutral model, given the expected level of recombination in this genomic region. Certain features of our data suggest that the unusual pattern in standard chromosomes is the product of selection rather than demography.

Molecular clock mirages

The hypothesis of the molecular clock proposes that molecular evolution occurs at rates that persist through time and across lineages, for a given gene. The neutral theory of molecular evolution predicts that the clock will be a Poisson process, with equal mean and variance. Experimental data have shown that the variance is typically larger than the mean. Hypotheses have been advanced to account for the hypervariance of molecular evolution. Four recent papers show that none of the predictive hypotheses that have been proposed can be generally maintained. The conclusion is that molecular evolution is dependent on the fickle process of natural selection. But it is a time-dependent process, so that accumulation of empirical data often yields an approximate clock, as a consequence of the expected convergence of large numbers.