Characterization and isolation of novel microsatellites from the Drosophila dunni subgroup

We have isolated and characterized 77 novel microsatellites from two species, Drosophila dunni and Drosophila nigrodunni, which are closely related Caribbean-island endemics from the Drosophila cardini species group. These species are very distantly related to all other Drosophila from which microsatellites have previously been characterized. We find that the average length of microsatellites isolated in these species is quite small, with an overall mean length of 9.8 repeat units for dinucleotide microsatellites in the two study species. The nucleotide composition of dinucleotides differs between the two species: D. nigrodunni has a predominance of (AC/GT)n repeats, whereas D. dunni has equal numbers of (AC/GT)n and (AG/CT)n repeats. Tri- and tetranucleotide repeats are not abundant in either species. We assayed the variability of eight microsatellites in a closely related third species, Drosophila arawakana, using wild-caught individuals from the island of Guadeloupe. We found the microsatellites to be extremely variable in this population, with observed heterozygosities ranging from 0.541 to 0.889. DNA amplification trials suggest that these eight microsatellites are widely conserved across the D. cardini group, with five of the eight producing amplification products in every species tested. However, the loci are very poorly conserved over greater phylogenetic distances. DNA amplification of the microsatellite loci was unreliable in members of the closely related Drosophila quinaria, Drosophila calloptera, Drosophila guarani and Drosophila tripunctata species groups. Furthermore, these microsatellites could not be detected in the genome of Drosophila melanogaster, despite the conservation of microsatellite flanking regions at some loci. These data indicate that Drosophila microsatellite loci are quite short lived over evolutionary timescales relative to many other taxa.     

Contrasting the efficacy of selection on the X and autosomes in Drosophila

To investigate the relative efficacy of both positive and purifying natural selection on the X chromosome and the autosomes in Drosophila, we compared rates and patterns of molecular evolution between these chromosome sets using the newly available alignments of orthologous genes from 12 species. Parameters that may influence the relative X versus autosomal substitution rates include the relative effective population sizes, the male and female germline mutation rates, the distribution of allelic effects on fitness, and the degree of dominance of novel mutations. Our analysis reveals that codon usage bias is consistently greater for X-linked genes, suggesting that purifying selection consistently has greater efficacy on the X chromosome than on the autosomes across the Drosophila phylogeny. However, our results are less consistent with respect to the efficacy of positive selection, with only some lineages showing a higher substitution rate on the X chromosome. This suggests that either the distribution of selective effects of mutations or other relevant parameters are sufficiently variable across species to tip the balance in different ways in individual lineages. These data suggest that rates of substitution are not solely governed by adaptive evolution. This genome-wide analysis provides a clear picture that the efficacy of selection varies intragenomically and that this effect is markedly more consistent across the phylogeny in the case of purifying selection. Our results also suggest that simple models that predict systematic differences in rates of evolution between the X and the autosomes can only be made to be compatible with these Drosophila data if the relevant population genetic parameters that drive substitution rates differ among species and chromosomal contexts.     

Microspatial structure of Drosophila melanogaster populations in Brazzaville: evidence of natural selection acting on morphometrical traits

Two genetically distinct habitat races of Drosophila melanogaster coexist in Brazzaville (Congo). One is the typical field type of Afrotropical populations, the other mainly breeds in beer residues in breweries. These two populations differ in their ethanol tolerance, in their allelic frequencies at several enzyme and microsatellite loci and in the composition of their cuticular hydrocarbons. The brewery population is quite similar to European temperate populations with regard to all these traits. Previous investigations of two morphological traits (ovariole number and sternopleural bristle number) failed to detect any difference between the two habitat races. Here we investigated other morphological traits (wing and thorax length, thorax pigmentation and female abdomen pigmentation). The reaction norms of these traits according to growth temperature were compared in the two Afrotropical habitat races and in a French temperate population. As expected, the French population was very different from the field African population: as a general rule, the brewery population (Kronenbourg) was intermediate in several aspects between the other two. We conclude that the strong selective forces that maintain the genetic divergence between the two habitat races also act on morphometrical traits, and the possible selective mechanisms are discussed.     

Molecular evolution of sex-biased genes in the Drosophila ananassae subgroup

Background  

Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura.     

Interactions between natural selection, recombination and gene density in the genes of Drosophila

In Drosophila, as in many organisms, natural selection leads to high levels of codon bias in genes that are highly expressed. Thus codon bias is an indicator of the intensity of one kind of selection that is experienced by genes and can be used to assess the impact of other genomic factors on natural selection. Among 13,000 genes in the Drosophila genome, codon bias has a slight positive, and strongly significant, association with recombination--as expected if recombination allows natural selection to act more efficiently when multiple linked sites segregate functional variation. The same reasoning leads to the expectation that the efficiency of selection, and thus average codon bias, should decline with gene density. However, this prediction is not confirmed. Levels of codon bias and gene expression are highest for those genes in an intermediate range of gene density, a pattern that may be the result of a tradeoff between the advantages for gene expression of close gene spacing and disadvantages arising from regulatory conflicts among tightly packed genes. These factors appear to overlay the more subtle effect of linkage among selected sites that gives rise to the association between recombination rate and codon bias.     

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.     

Contrasting mutual sexual selection on homologous signal traits in Drosophila serrata

The nature of male mating preferences, and how they differ from female mating preferences in species with conventional sex roles, has received little attention in sexual selection studies. We estimated the form and strength of sexual selection as a consequence of male and female mating preferences in a laboratory-based population of Drosophila serrata. The differences between sexual selection on male and female signal traits (cuticular hydrocarbons [CHCs]) were evaluated within a formal framework of linear and nonlinear selection gradients. Females tended to exert linear sexual selection on male CHCs, whereas males preferred intermediate female CHC phenotypes leading to convex (stabilizing) selection gradients. Possible mechanisms determining the nonlinear nature of sexual selection on female CHCs are proposed.     

Contrasting patterns of natural variation in global Drosophila melanogaster populations

Despite the popularity of Drosophila melanogaster in functional and evolutionary genetics, the global pattern of natural variation has not yet been comprehensively described in this species. For the first time, we report a combined survey using neutral microsatellites and mitochondrial sequence variation jointly. Thirty-five populations originating from five continents were compared. In agreement with previous microsatellite studies, sub-Saharan African populations were the most variable ones. Consistent with previous reports of a single 'out of Africa' habitat expansion, we found that non-African populations contained a subset of the African alleles. The pattern of variation detected for the mitochondrial sequences differed substantially. The most divergent haplotypes were detected in the Mediterranean region while Africa harboured most haplotypes, which were all closely related. In the light of the well-established African origin of D. melanogaster, our results cast severe doubts about the suitability of mtDNA for biogeographic inference in this model organism.     

Abdominal pigmentation in Drosophila melanogaster females from natural Indian populations

Females of Drosophila melanogaster collected from five geographically distant populations in India were analysed for the intensity of pigmentation in the 5th, 6th and 7th segments of the abdomen. In all three segments, this intensity was found to vary among individuals of any given population, and, furthermore, different populations differ with respect to this phenotypic trait. Statistical analysis revealed significant intra- and interpopulational variation. A clinical pattern was detected: females from populations closer to the equator tended to have lighter cuticle, in which case differences between the three segments could not be detected and all three segments responded both independently and jointly to latitudinal variation, as indicated by a statistically significant positive correlation between latitude and pigmentation score. This is the first report on abdominal pigmentation analysis in natural populations of D. melanogaster that provides evidence that phenotypic flexibility reflects temperature differences, as a result of which abdominal pigmentation shows geographic differentiation.     

Genetic basis of natural variation in body pigmentation in Drosophila melanogaster

Body pigmentation in insects and other organisms is typically variable within and between species and is often associated with fitness. Regulatory variants with large effects at bab1, t and e affect variation in abdominal pigmentation in several populations of Drosophila melanogaster. Recently, we performed a genome wide association (GWA) analysis of variation in abdominal pigmentation using the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP). We confirmed the large effects of regulatory variants in bab1, t and e; identified 81 additional candidate genes; and validated 17 candidate genes (out of 28 tested) using RNAi knockdown of gene expression and mutant alleles. However, these analyses are imperfect proxies for the effects of segregating variants. Here, we describe the results of an extreme quantitative trait locus (xQTL) GWA analysis of female body pigmentation in an outbred population derived from light and dark DGRP lines. We replicated the effects on pigmentation of 28 genes implicated by the DGRP GWA study, including bab1, t and e and 7 genes previously validated by RNAi and/or mutant analyses. We also identified many additional loci. The genetic architecture of Drosophila pigmentation is complex, with a few major genes and many other loci with smaller effects.     

Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes

Taste reception is fundamental to diet selection in many animals. The genetic basis underlying the evolution and diversity of taste reception, however, is not well understood. Recent discoveries of T1R sweet/umami receptor genes and T2R bitter receptor genes in humans and mice provided an opportunity to address this question. Here, we report the identification of 20 putatively functional T1R genes and 167 T2R genes from the genome sequences of nine vertebrates, including three fishes, one amphibian, one bird, and four mammals. Our comparative genomic analysis shows that orthologous T1R sequences are relatively conserved in evolution and that the T1R gene repertoire remains virtually constant in size across most vertebrates, except for the loss of the T1R2 sweet receptor gene in the sweet-insensitive chicken and the absence of all T1R genes in the tongueless western clawed frog. In contrast, orthologous T2R sequences are more variable, and the T2R repertoire diverges tremendously among species, from only three functional genes in the chicken to 49 in the frog. These evolutionary patterns suggest the relative constancy in the number and type of sweet and umami tastants encountered by various vertebrates or low binding specificities of T1Rs but a large variation in the number and type of bitter compounds detected by different species. Although the rate of gene duplication is much lower in T1Rs than in T2Rs, signals of positive selection are detected during the functional divergences of paralogous T1Rs, as was previously found among paralogous T2Rs. Thus, functional divergence and specialization of taste receptors generally occurred via adaptive evolution.     

Contrasting effects of natural selection on human and chimpanzee CC chemokine receptor 5

Human immunodeficiency virus type 1 (HIV-1) evolved via cross-species transmission of simian immunodeficiency virus (SIVcpz) from chimpanzees (Pan troglodytes). Chimpanzees, like humans, are susceptible to infection by HIV-1. However, unlike humans, infected chimpanzees seldom develop immunodeficiency when infected with SIVcpz or HIV-1. SIVcpz and most strains of HIV-1 require the cell-surface receptor CC chemokine receptor 5 (CCR5) to infect specific leukocyte subsets, and, subsequent to infection, the level of CCR5 expression influences the amount of HIV-1 entry and the rate of HIV-1 replication. Evidence that variants in the 5' cis-regulatory region of CCR5 (5'CCR5) affect disease progression in humans suggests that variation in CCR5 might also influence the response of chimpanzees to HIV-1/SIVcpz. To determine whether patterns of genetic variation at 5'CCR5 in chimpanzees are similar to those in humans, we analyzed patterns of DNA sequence variation in 37 wild-born chimpanzees (26 P. t. verus, 9 P. t. troglodytes, and 2 P. t. schweinfurthii), along with previously published 5'CCR5 data from 112 humans and 50 noncoding regions in the human and chimpanzee genomes. These analyses revealed that patterns of variation in 5'CCR5 differ dramatically between chimpanzees and humans. In chimpanzees, 5'CCR5 was less diverse than 80% of noncoding regions and was characterized by an excess of rare variants. In humans, 5'CCR5 was more diverse than 90% of noncoding regions and had an excess of common variants. Under a wide range of demographic histories, these patterns suggest that, whereas human 5'CCR5 has been subject to balancing selection, chimpanzee 5'CCR5 has been influenced by a selective sweep. This result suggests that chimpanzee 5'CCR5 might harbor or be linked to functional variants that influence chimpanzee resistance to disease caused by SIVcpz/HIV-1.     

Variation in Drosophila melanogaster central metabolic genes appears driven by natural selection both within and between populations

In this report, we examine the hypothesis that the drivers of latitudinal selection observed in the eastern US Drosophila melanogaster populations are reiterated within seasons in a temperate orchard population in Pennsylvania, USA. Specifically, we ask whether alleles that are apparently favoured in northern populations are also favoured early in the spring, and decrease in frequency from the spring to autumn with the population expansion. We use SNP data collected for 46 metabolic genes and 128 SNPs representing the central metabolic pathway and examine for the aggregate SNP allele frequencies whether the association of allele change with latitude and that with increasing days of spring–autumn season are reversed. Testing by random permutation, we observe a highly significant negative correlation between these associations that is consistent with this expectation. This correlation is stronger when we confine our analysis to only those alleles that show significant latitudinal changes. This pattern is not caused by association with chromosomal inversions. When data are resampled using SNPs for amino acid change the relationship is not significant but is supported when SNPs associated with cis-expression are only considered. Our results suggest that climate factors driving latitudinal molecular variation in a metabolic pathway are related to those operating on a seasonal level within populations.     

Molecular evolution of the 5'-flanking regions of the duplicated Amy genes in Drosophila melanogaster species subgroup

The nucleotide sequences of the 5'-flanking regions of the duplicated Amy genes in eight sibling species belonging to the melanogaster species subgroup are analyzed. In Drosophila melanogaster, a region of about 450 bp immediately upstream of the translation initiation site of the two paralogous genes (the proximal and distal genes) has sequence similarities. However, we could not detect any significant sequence similarity in the region more upstream than -450. This result indicates that the coding regions of the ancestral Amy gene were duplicated together with 450 bp of the 5'-flanking region as one unit. Multiple alignment of these 450-bp sequences in the proximal and distal genes of all eight species revealed a mosaic pattern of highly conserved and divergent regions. The conserved regions included almost all the putative regulatory elements identified in previous analyses of the sequences. A phylogenetic analysis of the aligned sequences shows that these 450-bp sequences are clustered into the proximal and the distal groups. As a whole, the divergence between groups in this region is very large in contrast to that in the coding regions. Based on the divergence between groups, the 450-bp region is divided into two subregions. We found that the ratios of the divergence between groups to that within groups differ in the two subregions. From these observations, we discuss a possibility of positive selection acting on the subregion immediately upstream of the Amy coding region to cause divergence of regulatory elements of the paralogous genes.