Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.     

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.     

Age-specific fitness components in hybrid females of Drosophila pseudoobscura and D. persimilis

Most models of hybridization assume that hybrids are less fit than their parental taxa. In contrast, some researchers have explored the possibility that hybrid individuals may actually have higher fitness and so play an important role in the generation of new species or adaptations. By estimating age-specific fitness components, we can determine not only how hybrid fitness differs from parental taxa, but also whether the fitness of hybrids relative to parental taxa changes with age. Here we describe an analysis of age-specific fitness traits in two species of Drosophila, D. pseudoobscura and D. persimilis, and their F1 hybrids. At early ages, hybrid females lay as many eggs as parental individuals, on average, but produce far fewer offspring. By late ages, in contrast, parental taxa show a steep decline in production not seen in hybrids, such that hybrids produce more offspring, on average, than parental taxa. Furthermore, egg-adult survival in hybrids is negatively correlated with egg density, whereas these traits are only weakly correlated in parental taxa. The results are limited somewhat by the fact that we analyze only two strains, and that these may be partially inbred. Nonetheless, the results are certainly illustrative, pointing out not only that at least some hybrid individuals may be as fit or fitter than parental taxa, but also that the difference between hybrids and parental taxa varies with age.     

Selection in heterogeneous environments maintains the gene arrangement polymorphism of Drosophila pseudoobscura

Chromosomal rearrangements may play an important role in how populations adapt to a local environment. The gene arrangement polymorphism on the third chromosome of Drosophila pseudoobscura is a model system to help determine the role that inversions play in the evolution of this species. The gene arrangements are the likely target of strong selection because they form classical clines across diverse geographic habitats, they cycle in frequency over seasons, and they form stable equilibria in population cages. A numerical approach was developed to estimate the fitness sets for 15 gene arrangement karyotypes in six niches based on a model of selection-migration balance. Gene arrangement frequencies in the six different niches were able to reach a stable meta-population equilibrium that matched the observed gene arrangement frequencies when recursions used the estimated fitnesses with a variety of initial inversion frequencies. These analyses show that a complex pattern of selection is operating in the six niches to maintain the D. pseudoobscura gene arrangement polymorphism. Models of local adaptation predict that the new inversion mutations were able to invade populations because they held combinations of two to 13 local adaptation loci together.     

Fine scale structural variants distinguish the genomes of Drosophila melanogaster and D. pseudoobscura

Background  

A primary objective of comparative genomics is to identify genomic elements of functional significance that contribute to phenotypic diversity. Complex changes in genome structure (insertions, duplications, rearrangements, translocations) may be widespread, and have important effects on organismal diversity. Any survey of genomic variation is incomplete without an assessment of structural changes.     

Genetics of Male and Female Sterility in Hybrids of Drosophila pseudoobscura and D. persimilis

The genetic basis of male and female sterility in hybrids of Drosophila pseudoobscura-Drosophila persimilis was studied using backcross analysis. Previous studies indirectly assessed male fertility by measuring testis size; these studies concluded that male sterility results from an X chromosome-autosome imbalance. By directly scoring for the production of motile sperm, male sterility is shown to be largely due to an incompatibility between genes on the X and Y chromosomes of these two species. These species have diverged at a minimum of nine loci affecting hybrid male fertility. Semisterility of hybrid females appears to result from an X chromosome-cytoplasm interaction; the X chromosome thus has the largest effect on sterility in both male and female hybrids. This is apparently the first analysis of the genetic basis of female sterility, or of sterility/inviability affecting both sexes, in an animal hybridization.     

Positive Selection on Nucleotide Substitutions and Indels in Accessory Gland Proteins of the Drosophila pseudoobscura Subgroup

Genes encoding reproductive proteins often diverge rapidly due to positive selection on nucleotide substitutions. While this general pattern is well established, the extent to which specific reproductive genes experience similar selection in different clades has been little explored, nor have possible targets of positive selection other than nucleotide substitutions, such as indels, received much attention. Here, we inspect for the signature of positive selection in the genes encoding five accessory gland proteins (Acps) (Acp26Aa, Acp32CD, Acp53Ea, Acp62F, and Acp70A) originally described from Drosophila melanogaster but with recognizable orthologues in the D. pseudoobscura subgroup. We compare patterns of selection within the D. psuedoobscura subgroup to those in the D. melanogaster subgroup. Similar patterns of positive selection were found in Acp26Aa and Acp62F in the two subgroups, while Acp53Ea and Acp70A experienced purifying selection in both subgroups. These proteins have thus remained targets for similar types of selection over long (>21-MY) periods of time. We also found several indel substitutions and polymorphisms in Acp26Aa and Acp32CD. These indels occur in the same regions as positively selected nucleotide substitutions for Acp26Aa in the D. pseudoobscura subgroup but not in the D. melanogaster subgroup. Rates of indel substitution within Acp26Aa in the D. pseudoobscura subgroup were up to several times those in noncoding regions of the Drosophila genome. This suggests that indel substitutions may be under positive selection and may play a key role in the divergence of some Acps. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Willis Swanson]     

Significance of Population Size on the Fixation of Nonsynonymous Mutations in Genes Under Varying Levels of Selection Pressure

Previous studies observed a higher ratio of divergences at nonsynonymous and synonymous sites (ω = d_N/d_S) in species with a small population size compared to that estimated for those with a large population size. Here we examined the theoretical relationship between ω, effective population size (N_e), and selection coefficient (s). Our analysis revealed that when purifying selection is high, ω of species with small N_e is much higher than that of species with large N_e. However the difference between the two ω reduces with the decline in selection pressure (s → 0). We examined this relationship using primate and rodent genes and found that the ω estimated for highly constrained genes of primates was up to 2.9 times higher than that obtained for their orthologous rodent genes. Conversely, for genes under weak purifying selection the ω of primates was only 17% higher than that of rodents. When tissue specificity was used as a proxy for selection pressure we found that the ω of broadly expressed genes of primates was up to 2.1-fold higher than that of their rodent counterparts and this difference was only 27% for tissue specific genes. Since most of the nonsynonymous mutations in constrained or broadly expressed genes are deleterious, fixation of these mutations is influenced by N_e. This results in a higher ω of these genes in primates compared to those from rodents. Conversely, the majority of nonsynonymous mutations in less-constrained or tissue-specific genes are neutral or nearly neutral and therefore fixation of them is largely independent of N_e, which leads to the similarity of ω in primates and rodents.     

Chromosomal and behavioral studies of Mexican Drosophila. III. Inversion polymorphism of D. pseudoobscura

Four new gene arrangements of chromosome 3 of Drosophila pseudoobscura are reported, as well as an updated phylogenetic sequence of the inversion polymorphism in this species. Evidence that the Tree Line gene arrangement of chromosome 3 was the original form in D. pseudoobscura, and that Central Mexico is the center of distribution of this species is discussed.     

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.     

Fast Accumulation of Nonsynonymous Mutations on the Female-Specific W Chromosome in Birds

Following cessation of recombination during sex chromosome evolution, the nonrecombining sex chromosome is affected by a number of degenerative forces, possibly resulting in the fixation of deleterious mutations. This might take place because of weak selection against recessive or partly recessive deleterious mutations due to permanent heterozygosity of nonrecombining chromosomes. Furthermore, population genetic processes, such as selective sweeps, background selection, and Muller’s ratchet, result in a reduction in N_e, which increase the likelihood of fixation of deleterious mutations. Theory thus predicts that nonrecombining genes should show increased levels of nonsynonymous (d_N) to synonymous substitutions (d_S). We tested this in an avian system by estimating the ratio between d_N and d_S in six gametologous gene pairs located on the Z chromosome and the nonrecombining, female-specific W chromosome. In comparisons, we found a significantly higher d_N/d_S ratio for the W-linked than the Z-linked copy in three of the investigated genes. In a concatenated alignment of all six genes, the d_N/d_S ratio was six times higher for W-linked than Z-linked genes. By using human and mouse as outgroup in maximum likelihood analyses, W-linked genes were found to evolve differently compared with their Z-linked gametologues and outgroup sequences. This seems not to be a consequence of functional diversification because d_N/d_S ratios between gametologous gene copies were consistently low. We conclude that deleterious mutations are accumulating at a high rate on the avian W chromosome, probably as a result of the lack of recombination in this female-specific chromosome. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

Evolutionary changes in the organization of the major LCP gene cluster during sex chromosomal differentiation in the sibling speciesDrosophila persimilis,D. pseudoobscura andD. miranda

The major larval cuticle protein (LCP) genes I–IV ofDrosophila melanogaster are clustered on the right arm of the second chromosome. By cross-hybridization we cloned the corresponding genes from three different members of theobscura group:D. persimilis, D. pseudoobscura andD. miranda. InD. pseudoobscura andD. persimilis the gene cluster maps to autosome3. In contrast, inD. miranda it was found on theX2 andY sex chromosome. Hence, this exceptional karyotypic situation offers a unique opportunity to analyse the molecular processes underlying the phenomenon of chromosome degeneration. Comparison of LCP genes I–IV in theX2 andY chromosomal region inD. miranda revealed extensive DNA rearrangements at the latter. TheY chromosomal LCP cluster is characterized by DNA insertions which are absent in the correspondingX2 chromosomal DNA, suggesting that these DNA sequences must have invaded this area. In addition, part of the analysedY chromosomal region is duplicated.     

Cytogenetic mapping of marker genes on the chromosome elements C and E of Drosophila pseudoobscura and D. subobscura

Enzyme loci, visible marker genes and -cloned DNA-sequences from a D. miranda library were mapped cytologically on the chromosome elements C and E of D. pseudoobscura and D. subobscura. New data are incorporated into the linkage maps of the two species. Homologous segments can now be localized in the polytene chromosomes with these markers. A comparison of the chromosome elements E of D. melanogaster and D. subobscura shows 12 conserved subsections which have been rearranged by paracentric inversions in the evolution of the two lineages.     

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.     

Telomere repeats within the neo-Y-chromosome of Drosophila miranda

The clone Dmir1098, isolated from a genomic lambda library of Drosophila miranda labels exclusively the tips of the giant chromosomes in the highly polytenized nuclei of the female larval salivary glands. However, the in situ hybridizations to male metaphase plates, using the same probe, reveal a massive labeling block within the neo-Y-chromosome in addition to the labeling blocks at both chromosome ends. From the comparison with the Y chromosome labeling pattern of D. pseudoobscura, a sibling species to D. miranda, an end-to-end fusion mechanism involving the telomere repeats would be the most straightforward explanation for the karyotype change in D. miranda.     

Microsatellite variation in populations of Drosophila pseudoobscura and Drosophila persimilis

We have isolated, characterized and mapped 33 dinucleotide, three trinucleotide and one tetranucleotide repeat loci from the four major chromosomes of Drosophila pseudoobscura. Average inferred repeat unit length of the dinucleotide repeats is 12 repeat units, similar to D. melanogaster. Assays of D. pseudoobscura and populations of its sibling species, D. persimilis, using 10 of these loci show extremely high levels of variation compared with similar studies of dinucleotide repeat variation in D. melanogaster populations. The high levels of variation are consistent with an average mutation rate of approximately 10(-6) per locus per generation and an effective population size of D. pseudoobscura approximately four times larger than that of D. melanogaster. Consistent with allozymes and nucleotide sequence polymorphism, the dinucleotide repeat loci reveal minimal structure across four populations of D. pseudoobscura. Finally, our preliminary recombinational mapping of 24 of these microsatellites suggests that the total recombinational genome size may be larger than previously inferred using morphological mutant markers.