Pervasive gene conversion in chromosomal inversion heterozygotes

Chromosomal inversions shape recombination landscapes, and species differing by inversions may exhibit reduced gene flow in these regions of the genome. Though single crossovers within inversions are not usually recovered from inversion heterozygotes, the recombination barrier imposed by inversions is nuanced by noncrossover gene conversion. Here, we provide a genomewide empirical analysis of gene conversion rates both within species and in species hybrids. We estimate that gene conversion occurs at a rate of 1 × 10^–5 to 2.5 × 10^–5 converted sites per bp per generation in experimental crosses within Drosophila pseudoobscura and between D. pseudoobscura and its naturally hybridizing sister species D. persimilis. This analysis is the first direct empirical assessment of gene conversion rates within inversions of a species hybrid. Our data show that gene conversion rates in interspecies hybrids are at least as high as within‐species estimates of gene conversion rates, and gene conversion occurs regularly within and around inverted regions of species hybrids, even near inversion breakpoints. We also found that several gene conversion events appeared to be mitotic rather than meiotic in origin. Finally, we observed that gene conversion rates are higher in regions of lower local sequence divergence, yet our observed gene conversion rates in more divergent inverted regions were at least as high as in less divergent collinear regions. Given our observed high rates of gene conversion despite the sequence differentiation between species, especially in inverted regions, gene conversion has the potential to reduce the efficacy of inversions as barriers to recombination over evolutionary time.     

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.     

Genetic and Evolutionary Correlates of Fine-Scale Recombination Rate Variation in Drosophila persimilis

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed “recombinational neighborhoods,” where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5–7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.     

Mapping of within-species segregation distortion in Drosophila persimilis and hybrid sterility between D. persimilis and D. pseudoobscura

In contrast to the prevailing dogma in the 1990s, recent studies have suggested that an evolutionary history of segregation distortion within species may contribute to sterility in species hybrids. However, this recent work identified segregation distortion exclusively in species hybrids that may never have had an evolutionary history of segregation distortion in either parent species. We expand on previous work using a strain of Drosophila persimilis exhibiting segregation distortion within species to generate QTL maps for segregation distortion and hybrid sterility in crosses between D. persimilis and D. pseudoobscura. The maps localize regions along the XR contributing to both phenotypes, and they indicate one region of overlap between the two maps. This overlap could provide preliminary evidence for an association between segregation distortion within species and hybrid sterility, but the localizations are currently too broad to have confidence in this conclusion. This work is a first step towards possibly supporting a genetic conflict model of speciation in this system.     

Connecting recombination,nucleotide diversity,and species divergence in Drosophila

The association between recombination rate and nucleotide diversity provides compelling evidence for the action of natural selection across much of the Drosophila melanogaster genome. This conclusion is further supported by the lack of association between recombination rate and nucleotide divergence between species. However, studies of other species, including other Drosophila, have not always yielded the same results. Our recent study measured these parameters within the D. pseudoobscura species group using next-generation sequencing and high-throughput genotyping technologies. We documented fine-scale variation in crossover rate within D. pseudoobscura, and we observed that crossover variation was strongly associated with nucleotide diversity only when measured at a fine-scale. We also observed associations between crossover rate and sequence differences between D. pseudoobscura and its close relatives. These latter associations could have been driven in part by mutagenic effects associated with double-strand break repair, but we cannot exclude the possibility that it results primarily from shared ancestral polymorphisms. Overall, this work strongly underscores the importance of scale in testing for associations of recombination rate with other parameters, and it brings us one small step closer to understanding the role of natural selection and other evolutionary forces in shaping divergence among genomes.     

The effect of relative humidity on photoresponse in sympatric species of Drosophila: Short term exposure to desiccating environments

Photoresponses of isofemale strains of sympatric sibling species, Drosophila pseudoobscura and D. persimilis from MacDonald Ranch, California and D. pseudoobscura and D. miranda from Mather, California were compared at relative humidities of 95, 50 and 5% for a period of 30 min. The sympatric species could be distinguished behaviourally by their responses at the three humidities. Although the average photoresponses of all four populations were negative, D. pseudoobscura was significantly more photonegative than its sibling sympatric species in each habitat. There was no effect of humidity on light response over the short term duration of exposure to desiccating environments. A hypothesis is presented to account for this phenomenon.     

Gene Differences between the Sex Ratio and Standard Gene Arrangements of the X Chromosome in DROSOPHILA PERSIMILIS

The sibling species Drosophila pseudoobscura and D. persimilis each carry two gene arrangements in the right arm of the X chromosome; Standard (ST) and Sex Ratio (SR). The SR sequence of D. persimilis and the ST sequence of D. pseudoobscura have the same banding pattern. These cytologically identical arrangements carry different alleles at the esterase-5 (Est-5) and phosphoglucomutase-1 ( Pgm-1) loci. All the alleles on the SR arrangement of D. persimilis are also present on the ST arrangement of that species, and the same allele was most common in all cases. However, the allele frequencies at the Pgm-1 locus are significantly different between the two arrangements. Within the ST arrangement of D. persimilis, the alleles at the Est-5 and Pgm-1 loci are in linkage equilibrium.     

Principles of genome evolution in the Drosophila melanogaster species group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

A Comparative Study of the Esterase-5 Locus in DROSOPHILA PSEUDOOBSCURA,D. PERSIMILIS and D. MIRANDA

Electrophoretic phenotypes of the esterase-5 locus were examined in the sibling species D. pseudoobscura, D. persimilis and D. miranda. D. persimilis alleles were found to have uniformly higher charge on monomers than corresponding alleles of either D. pseudoobscura or D. miranda. Consequently, D. persimilis shares no alleles in common with either D. pseudoobscura or D. miranda, while the latter two species share a number of alleles. It was discovered that by increasing the concentration of acrylamide gel and increasing the length of migration, more allelic differences could be distinguished. Also more alleles were discovered by examining monomer mobility in addition to dimer mobility. In D. persimilis and D. miranda it was found that the previously known high frequency allelic classes broke down into several allelic classes. A test of goodness-of-fit to the infinite alleles model was done and a rough agreement with the model was found.     

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.     

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement‐specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.     

P1 clones from Drosophila melanogaster as markers to study the chromosomal evolution of Muller's A element in two species of the obscura group of Drosophila

Thirty P1 clones from the X chromosome (Muller's A element) of Drosophila melanogaster were cross-hybridized in situ to Drosophila subobscura and Drosophila pseudoobscura polytene chromosomes. An additional recombinant phage Dsuby was also used as a marker. Twenty-three (77%) of the P1 clones gave positive hybridization on D. pseudoobscura chromosomes bat only 16 (53%) did so with those of D. subobscura. Eight P1 clones gave more than one hybridization signal on D. pseudoobscura and/or D. subobscura chromosomes. All P1 clones and Dsuby hybridized on Muller's A element (X chromosome) of D. subobscura. In contrast, only 18 P1 clones and Dsuby hybridized on Muller's D element (XR chromosomal arm) of D. pseudoobscura; 4 additional P1 clones hybridized on Muller's D element (XR chromosomal arm) of this species and the remaining P1 clone gave on hybridization signal on each arm of the X chromosome. This latter clone may contain one breakpoint of a pericentric inversion that may account for the interchange of genetic material between Muller's A and D elements in D. pseudoobscura. In contrast to the rare interchange of genetic material between chromosomal elements, profound differences in the order and spacing of markers were detected between D. melanogaster, D. pseudoobscura and D. subobscura. In fact, the number of chromosomal segments delimited by identical markers and conserved between pairwise comparisons is small. Therefore, extensive reorganization within Muller's A element has been produced during the divergence of the three species. Rough estimates of the number of cytologically detectable inversions contributing to differentiation of Muller's A element were obtained. The most reliable of these estimates is that obtained from the D. pseudoobscura and D. melanogaster comparison since a greater number of markers have been mapped in both species. Tentatively, one inversion breakpoint about every 200 kb has been produced and fixed during the divergence of D. pseudoobscura and D. melanogaster.     

GENETICS OF STERILITY IN HYBRIDS BETWEEN TWO SUBSPECIES OF DROSOPHILA

Hybrids between D. pseudoobscura bogotana and D. pseudoobscura pseudoobscura are fertile except for males produced in one of the two reciprocal crosses. As there is no premating isolation between these subspecies, nonreciprocal male sterility represents the first step in speciation. Genetic analysis reveals two causes of hybrid F₁ sterility: a maternal effect and incompatibilities between chromosomes within males. The maternal effect appears to play the greatest role in hybrid sterility. The X chromosome has the largest effect on fertility of any chromosome, a ubiquitous result in analyses of hybrid sterility and inviability in Drosophila. This effect is entirely attributable to a region comprising less than 30% of the X chromosome. These results are compared to those from a similar study of D. pseudoobscura-D. persimilis hybrids, an older and more reproductively isolated species pair in the same lineage. Such comparisons may allow one to identify the genetic changes characterizing the early versus late stages of speciation.     

Multifaceted,Cross-Generational Costs of Hybridization in Sibling Drosophila species

Maladaptive hybridization, as determined by the pattern and intensity of selection against hybrid individuals, is an important factor contributing to the evolution of prezygotic reproductive isolation. To identify the consequences of hybridization between Drosophila pseudoobscura and D. persimilis, we estimated multiple fitness components for F1 hybrids and backcross progeny and used these to compare the relative fitness of parental species and their hybrids across two generations. We document many sources of intrinsic (developmental) and extrinsic (ecological) selection that dramatically increase the fitness costs of hybridization beyond the well-documented F1 male sterility in this model system. Our results indicate that the cost of hybridization accrues over multiple generations and reinforcement in this system is driven by selection against hybridization above and beyond the cost of hybrid male sterility; we estimate a fitness loss of >95% relative to the parental species across two generations of hybridization. Our findings demonstrate the importance of estimating hybridization costs using multiple fitness measures from multiple generations in an ecologically relevant context; so doing can reveal intense postzygotic selection against hybridization and thus, an enhanced role for reinforcement in the evolution of populations and diversification of species.     

A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

Principles of Genome Evolution in the Drosophila melanogaster Species Group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

Genic Heterogeneity at Two Alcohol Dehydrogenase Loci in DROSOPHILA PSEUDOOBSCURA and DROSOPHILA PERSIMILIS

A sequential electrophoretic survey of the second chromosome loci, alcohol dehydrogenase-6 (Adh-6) and octanol dehydrogenase ( Odh), was performed on 147 isochromosomal lines of Drosophila pseudoobscura and 60 lines of its sibling species, D. persimilis. Gels run with a variety of acrylamide concentrations and buffer pH's revealed the presence of 18 alleles of Adh-6 in the two species, where only eight had been previously detected by conventional electrophoretic methods. Only two alleles were added with our techniques to the previous total of nine in both species at the largely monomorphic Odh locus. Both enzymes show a predominance of one allele, with the other variants being fairly rare. There was no evidence of increased genetic divergence between the two species, but we found a striking increase in differentiation of Adh-6 alleles between the main body of D. pseudoobscura populations and the conspecific isolate from Bogotá, Colombia. These results are compared with our previous surveys of xanthine dehydrogenase in these species and discussed in reference to theories of genic polymorphism.     

Evidence for Extensive Genetic Differentiation between the Sex-Ratio and the Standard Arrangement of DROSOPHILA PSEUDOOBSCURA and D. PERSIMILIS and Identification of Hybrid Sterility Factors

This study deals with sex-ratio genes tightly linked within the Sex-Ratio inversion. By taking advantage of the fact that the Sex-Ratio chromosome of Drosophila persimilis [SR(B)] is homosequential to the Standard chromosome of D. pseudoobscura [ST(A)], we carried out two reciprocal introgression experiments. Individual segments of SR(B) or ST(A) were introgressed into the genome of D. pseudoobscura or D. persimilis, respectively. Males possessing a hybrid SR(B)-ST(A) X chromosome and a genetic background derived from either of the two species were tested for fertility and sex-ratio expression.—It was found that, in terms of the meiotic drive genes, the Sex-Ratio chromosome differs extensively from the Standard chromosome. Because recombinations of these genes result in a complete loss of sex-ratio expression, this finding lends strong support to the hypothesis of gene coadaptation. Coadaptation, in this context, is the advantage of being transmitted preferentially. In light of this finding, the evolution of the sex-ratio system in these two sibling species is discussed.—Introgression experiments also yielded information about hybrid sterility. With reciprocal introgression, sterility interactions were found to be "asymmetric." The asymmetry is fully expected from the viewpoint of evolution of postmating reproductive isolation.     

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

Background

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

Methodology/Principal Findings

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

Conclusions

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