Recombination and selection in the maintenance of the adaptive value of inversions

A huge amount of data seem to confirm the adaptive value of inversions in Drosophila. The inhibition of recombination in heterokaryotypes mediated by inversions seems fundamental in maintaining their adaptive role. This study shows that recombination is highly suppressed in Drosophila subobscura because of chromosomal inversions, not only inside the inversions but also outside them. It seems that the region outside the inversion where recombination is inhibited is asymmetrical and independent of the inversion length. Despite the difficulty of crossovers taking place near inversion breakpoints, the only two recombination events detected inside inversions were located close to the breakpoint. Thus, selection could be largely responsible for the recombination reduction maintaining sets of adaptive alleles inside the inverted region. Heterokaryotype descendants were always in higher frequency than inbred or outbred homokaryotypes, regardless of the geographical origin of the chromosome, suggesting that chromosomes carrying the same arrangement, although with a different set of alleles for neutral markers, could be submitted to the same selection processes.     

Variation in Activities of Amylase Allozymes Associated with Chromosome Inversions in DROSOPHILA PSEUDOOBSCURA, D. PERSIMILIS and D. MIRANDA

Different electrophoretic alleles of amylase show associations with particular chromosome 3 inversions in D. pseudoobscura and D. persimilis. Relative adult amylase activities were compared in 37, 37 and 10 strains of D. pseudoobscura, D. persimilis and D. miranda, respectively. Strains carrying the same electrophoretic allele were compared by crossing these lines individually to a reference strain carrying a different electrophoretic mobility allele. This procedure allows comparisons among species, inversions, electromorphs and strains for genetic variation in amylase activity. F2 analysis established that the activity variation co-segregates with the structural amylase locus. This type of variation could be due to either structural gene differences or differences in closely linked, cis-acting regulatory regions. Variation has been detected among and within electrophoretic mobility classes. Moreover, this variation is clearly nonrandom and reveals more of the genetic structure associated with the chromosomal inversion phylogeny of D. pseudoobscura and D. persimilis. ----Some of the findings are: (1) Similar electromorphs in D. pseudoobscura and D. persimilis usually show different activities. These species show nearly complete differentiation of amylase alleles, based on activities. (2) D. persimilis has the broadest range of variation in amylase activity, about four-fold between the highest and lowest alleles. D. pseudoobscura and D. miranda are also polymorphic for activity, but have more constrained ranges of variation. D. miranda alleles show on the average about four times the activity of D. pseudoobscura alleles. (3) Some association of electrophoretic mobility and activity has been found. Alleles 1.09 of D. persimilis, as well as 1.43 and 1.55 of D. miranda, have relatively high activity. It may be that these high activity alleles are part of an adaptation to cooler habitats. (4) Within electrophoretic classes, associations of activities with inversions have been found. These are especially strong in D. persimilis. The 1.00 alleles in the ST, KL, MD and WT inversions, the 0.92 allele in the ST and MD inversions and the 1.09 allele in the WT and KL inversions have levels of activities that depend upon the arrangement in which they are located. These results demonstrate that suppression of recombination in inversion heterokaryotypes can result in extensive genic divergence between inversions.     

Allozyme divergence between different gene arrangements in Drosophila subobscura

Allozyme and chromosomal polymorphism were analysed simultaneously in three population samples of Drosophila subobscura from Scotland, Germany and Italy. The study revealed significant genic divergence between the various chromosomal gene arrangements for the loci Acph, Lap and Pep, which are located inside inversions. No difference was found with the loci Me and Odh, which are both on the same chromosome but proximal to the structurally polymorphic region. Samples of a given gene arrangement taken from the different localities are genetically uniform. The data indicate that the observed non-random associations between enzyme genes and chromosomal inversions are the result of inhibited recombination in heterokaryotypes and gradual differentiation of the separated gene pools, rather than selection and coadaptation.     

High Fitness of Heterokaryotypic Individuals Segregating Naturally within a Long-Standing Laboratory Population of Drosophila silvestris

Natural populations of Drosophila silvestris are polymorphic for inversions in one or more of four of the five major chromosome arms; laboratory stocks tend to retain this heterozygosity. A laboratory stock, U28T2, was started from a single naturally inseminated wild female caught at Kilauea Forest Reserve, Hawaii, in January 1977. Polytene analysis in 1980 showed the presence of three natural inversions in chromosome 4: k ² is distal, t is central and l² is proximal. The inversions are short but only short uncovered euchromatic sections exist at the distal and proximal ends. Periodic examinations through 1986 showed all three inversions to be persistent at moderately high frequencies. In 1984, a series of tests of mating performance of caged, mature males, taken at random as they eclosed from the stock, were followed by cytological testcrosses to females from a homokaryotypic stock. Only three of the eight possible haplotypes, k²/t/+ (A), +/+/l² ( a) and +/+/+ (a') were present. Tests of crossing over show none in males; in females, there is about 1% in each of the two regions between the inversions. Only one such apparent crossover haplotype was found among 1084 examined in samples from this stock. Thus, chromosome arrangements A, a and a' virtually behave as wholechromosome alleles in both sexes. Of 146 males marked and tested in cages, 61 produced progeny; the others failed to reproduce. Of 58 males and 80 females producing progeny and analyzed cytologically, there were, respectively, 49 and 59 heterokaryotypes. On the basis of frequencies calculated for fertilized eggs, 33.6 males and 46.3 females are expected. The facts suggest that individual males with the Aa karyotype are particularly successful in production of offspring. Adult females show an excess of Aa' as well as Aa. Such high fitness of heterokaryotypes in the effective breeding adults could be a major factor in the maintenance of stable chromosomal polymorphisms both in laboratory stocks and in nature. Although some of this heterosis is clearly ascribable to differential survival, the facts suggest that there is a substantial opportunity, indeed a likelihood, for a contribution from differential mating among surviving adults.     

Molecular population genetics of the rp49 gene region in different chromosomal inversions of Drosophila subobscura.

Nucleotide variation at the ribosomal protein 49 (rp49) gene region has been studied in 75 lines of Drosophila subobscura belonging to four chromosomal arrangements (Ost, O3+4, O3+4+8, and O3+4+23). The location of the rp49 gene region within the inversion loop differs among heterokaryotypes: it is very close to one of the breakpoints in heterozygotes involving Ost chromosomes, while it is in a more central position in all other heterokaryotypes. The distribution of nucleotide polymorphism in the different arrangements is consistent with a monophyletic origin of the inversions. The data also provide evidence that gene conversion and possibly double crossover are involved in shuffling nucleotide variation among gene arrangements. The analyses reveal that the level of genetic exchange is higher when the region is located in a more central position of the inverted fragment than when it is close to the breakpoints. The pairwise difference distributions as well as the negative values of Tajima''s and Fu and Li''s statistics further support the hypothesis that nucleotide variation within chromosomal arrangements still reflects expansion after the origin of the inversions. Under the expansion model, we have estimated the time of origin of the studied inversions.     

Genomic impacts of chromosomal inversions in parapatric Drosophila species

Chromosomal inversions impact genetic variation and facilitate speciation in part by reducing recombination in heterokaryotypes. We generated multiple whole-genome shotgun sequences of the parapatric species pair Drosophila pseudoobscura and Drosophila persimilis and their sympatric outgroup (Drosophila miranda) and compared the average pairwise differences for neutral sites within, just outside and far outside of the three large inversions. Divergence between D. pseudoobscura and D. persimilis is high inside the inversions and in the suppressed recombination regions extending 2.5 Mb outside of inversions, but significantly lower in collinear regions further from the inversions. We observe little evidence of decreased divergence predicted to exist in the centre of inversions, suggesting that gene flow through double crossovers or gene conversion is limited within the inversion, or selection is acting within the inversion to maintain divergence in the face of gene flow. In combination with past studies, we provide evidence that inversions in this system maintain areas of high divergence in the face of hybridization, and have done so for a substantial period of time. The left arm of the X chromosome and chromosome 2 inversions appear to have arisen in the lineage leading to D. persimilis approximately 2 Ma, near the time of the split of D. persimilis-D. pseudoobscura-D. miranda, but likely fixed within D. persimilis much more recently, as diversity within D. persimilis is substantially reduced inside and near these two inversions. We also hypothesize that the inversions in D. persimilis may provide an empirical example of the 'mixed geographical mode' theory of inversion origin and fixation, whereby allopatry and secondary contact both play a role.     

Genes in Neurospora That Suppress Recombination When They Are Heterozygous

Genes that suppress recombination when heterozygous have been found distributed as a polymorphism in wild and laboratory populations of Neurospora crassa. Three alleles, ss^E, ss^S and ss^C, are associated, respectively, with the three wild types Emerson, St. Lawrence 74A and Costa Rica A. It is proposed that ss (synaptic sequence) genes modulate recombination by determining the pairing closeness of DNA duplexes in the vicinity of the nit-2 locus. When heterozygous, ss suppresses recombination 2- to 20-fold within the nit-2 locus, which it adjoins, but crossing over in intervals flanking nit-2 is not affected. The magnitude of suppression depends upon the ss alleles involved, and ss acts multiplicatively with rec-1; together, these genes modulate recombination within the nit-2 locus over a range exceeding 100-fold. The ss effect is not attributable to gross chromosomal rearrangement, but could be due to small inversions or insertions, such as transposable elements.     

Gametic associations between inversion and allozyme polymorphisms in Drosophila buzzatii.

Gametic disequilibria between second chromosome polymorphic arrangements and seven linked allozyme loci were estimated in seven populations of Drosophila buzzatii from Argentina. Significant and consistent associations across populations were detected for Est-1, Est-2, Aldox, and XDH: Phenograms based on Nei's genetic distance showed that chromosomes carrying the 2ST arrangement were more similar to each other, irrespective of the population from which they were extracted, than to chromosomes carrying the derived 2J and 2JZ3.Restriction of recombination in heterokaryotypes seems to be the best explanation for the significant linkage disequilibria between inversions and the loci located inside the rearranged segments, for example, Est-1 and Aldox, or close to the break points, for example, Est-2. However, epistatic interactions between Xdh, which is outside the inversions and not near the break points, and loci tightly linked to the inversions, is the most likely explanation for the association between Xdh and chromosomal arrangements. Some of the associations detected in endemic Argentinean populations are coincident with data obtained in colonizing populations of the Old World and Australia. Thus historical processes that took place in the original area of the species' distribution can account for these linkage disequilibria in colonized populations of D. buzzatii.     

The evolutionary history of Drosophila buzzatii. XXXV. Inversion polymorphism and nucleotide variability in different regions of the second chromosome

Inversions are portions of a chromosome where the gene order is reversed relative to a standard reference orientation. Because of reduced levels of recombination in heterokaryotypes, inversions have a potentially important effect on patterns of nucleotide variability in those genomic regions close to, or included in, the inverted fragments. Here we report sequence variation at three anonymous regions (STSs) located at different positions in relation to second-chromosome inversion breakpoints in 29 isochromosomal lines derived from an Argentinean population of Drosophila buzzatii. In agreement with previous findings in Drosophila, gene flux (crossing over and/or gene conversion) between arrangements seems to appreciably increase as we approach the middle sections of inversion 2j, and patterns of nucleotide variability within, as well as genetic differentiation between chromosome arrangements, are comparable to those observed at the molecular marker outside the inverted fragments. On the other hand, nucleotide diversity near the proximal breakpoint of inversion 2j is reduced when contrasted with that found at the other regions, particularly in the case of derived inverted chromosomes. Using the data from the breakpoint, we estimate that the inversion polymorphism is approximately 1.63 N generations old, where N is the effective population size. An excess of low-frequency segregating polymorphisms is detected; mostly in the ancestral 2st arrangement and probably indicating a population expansion that predates the coalescent time of inversion 2j. Heterogeneity in mutation rates between the markers linked to the inversions may be sufficient to explain the different levels of nucleotide diversity observed. When considered in the context of other studies on patterns of variation relative to physical distance to inversion breakpoints, our data appear to be consistent with the conclusion that inversions are unlikely to be "long-lived" balanced polymorphisms.     

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.     

Mapping regions within cosmopolitan inversion In(3R)Payne associated with natural variation in body size in Drosophila melanogaster

Associations between genotypes for inversions and quantitative traits have been reported in several organisms, but little has been done to localize regions within inversions controlling variation in these traits. Here, we use an association mapping technique to identify genomic regions controlling variation in wing size within the cosmopolitan inversion In(3R)Payne in Drosophila melanogaster. Previous studies have shown that this inversion strongly influences variation in wing size, a trait highly correlated with body size. We found three alleles from two separate regions within In(3R)Payne with significant additive effects on wing size after the additional effect of the inversion itself had been taken into account. There were also several alleles with significant genotype-by-inversion interaction effects on wing size. None of the alleles tested had a significant additive effect on development time, suggesting different genes control these traits and that clinal patterns in them have therefore arisen independently. The presence of multiple regions within In(3R)Payne controlling size is consistent with the idea that inversions persist in populations because they contain multiple sets of locally adapted alleles, but more work needs to be done to test if they are indeed coadapted.     

Chromosome inversions, local adaptation and speciation

We study the evolution of inversions that capture locally adapted alleles when two populations are exchanging migrants or hybridizing. By suppressing recombination between the loci, a new inversion can spread. Neither drift nor coadaptation between the alleles (epistasis) is needed, so this local adaptation mechanism may apply to a broader range of genetic and demographic situations than alternative hypotheses that have been widely discussed. The mechanism can explain many features observed in inversion systems. It will drive an inversion to high frequency if there is no countervailing force, which could explain fixed differences observed between populations and species. An inversion can be stabilized at an intermediate frequency if it also happens to capture one or more deleterious recessive mutations, which could explain polymorphisms that are common in some species. This polymorphism can cycle in frequency with the changing selective advantage of the locally favored alleles. The mechanism can establish underdominant inversions that decrease heterokaryotype fitness by several percent if the cause of fitness loss is structural, while if the cause is genic there is no limit to the strength of underdominance that can result. The mechanism is expected to cause loci responsible for adaptive species-specific differences to map to inversions, as seen in recent QTL studies. We discuss data that support the hypothesis, review other mechanisms for inversion evolution, and suggest possible tests.     

Viabilities of originally natural O-chromosomal inversion homo- and heterokaryotypes in Drosophila subobscura

The influence of epistatic interactions of lethal and non-lethal genes upon viability of Drosophila inversion karyotypes is poorly known. In this paper we present comparative results for viabilities of 21 originally natural O-inversion homo- and 38 heterokaryotypes in a D. subobscura population. We observed strong heterotic effect in viability of O-lethal heterozygotes irrespective of different inversion backgrounds, which indicates a mechanism for protection of a great number of lethal genes. In O-non-lethal heterozygotes the heterotic effect in viability was absent irrespective of different inversion backgrounds. Our results showed a great extent of genetic load and high abundance of O-chromosomal arrangements in the population analyzed. It belongs to the set of central European populations. An epistasis of lethal genes present in O-inversion hetero- and not present in O-inversion homokaryotypes of moderate or low frequencies could be good example for co-adaptation of chromosomal inversions with regard to the genetic load. This represents a more efficient mechanism for limitation of genetic load than alternative mechanisms for protection of lethals. Except for lethal genes, possible epistatic interactions of mildly deleterious (subvital) genes, could also be responsible for limiting the extent of genetic load in natural D. subobscura populations. We demonstrated a non-random distribution of several combinations of viability classes among different O-inversion homo- and heterokaryotypes. As a consequence of that, the viabilities of the O-inversion homokaryotypes compared to heterokaryotypes were significantly higher at low frequencies than in moderate or high frequencies. This shows frequency-dependence as a mechanism of balancing selection for protection of O-chromosomal inversions in natural D. subobscura populations. In addition, the viabilities of the O-inversion homokaryotypes of lower frequency, compared to homokaryotypes of moderate or high frequency, were significantly higher. This again indicates the existence of supergene selection as another mechanism for protection of chromosomal inversions, as co-adapted complexes in natural D. subobscura populations.     

The Genetics of DROSOPHILA SUBOBSCURA Populations. V. a Study of Linkage Disequilibrium in Natural Populations between Genes and Inversions of the E Chromosome

The genetics of Hk and Est-9 complex gene have been studied in Drosophila subobscura. While Hk alleles mendelize normally, Est-9 is a complex locus consisting of several very closely linked genes with active and silent alleles. Both genes are located on chromosome E; a detailed genetic map was constructed with the help of visible markers and inversions. Both Hk and Est-9 are included in or are very near to inversions of the E chromosome. While Hk does not show linkage disequilibrium either with Est-9 or inversions, Est-9 does show disequilibrium in two ways: both with inversions and between different Est-9 genes. All natural populations studied show the same kind of association between Est-9 gene combinations and inversions. It is argued that these results are better explained by selection than by neutrality.     

Analysis of male recombination in third chromosomes of Drosophila melanogaster

Data on male recombination in twenty third-chromosomal lines of Drosophila melanogaster are presented. Frequencies of female and male recombination have been calculated in seven intervals along the third chromosome. The influence on male recombination (M.R.) exercised by different factors such as population origin (cellar, vineyard), the presence of heterozygous inversions and recessive lethal chromosomes, is analyzed. The results obtained lead to the main conclusion that M.R. is not affected by the presence of heterozygous inversions which reduce female recombination in the same lines. In the light of this effect, the possible mechanism operating on male recombination is discussed. Lethal chromosomes reduce significantly the number of male recombination events as compared with wild chromosomes. We have not obtained significant differences in male recombination frequencies between the cellar and the vineyard lines.     

Chromosomal inversion polymorphism of Anopheles funestus from forest villages of South Cameroon.

The polymorphism of paracentric inversions of Anopheles funestus polytene chromosomes was studied in three villages (Nkoteng, Obala, and Simbock) located in a forest area of South Cameroon in order to analyse the genetic structure of these populations. A total of 146-210 chromatids could be scored from specimens collected over about two years. A low degree of chromosomal polymorphism was observed with two floating inversions on chromosomal arm 2 (2h and 2d), and three fixed arrangements on arms 3 (3a and 3b), and 5 (5a). Such arrangement of inversions has never been recorded elsewhere so far. The chromosome analysis indicated that the population from Obala was in Hardy-Weinberg equilibrium, whereas the samples from Nkoteng and Simbock showed a significant excess and deficit of heterokaryotypes, respectively. Significant differences in inversion frequencies on chromosomal arm 2 among villages lying in contrasting eco-climatic settings suggested an adaptive role of these inversions.     

Cre-loxP Recombination System for Large Genome Rearrangements in Lactococcus lactis

We have used a new genetic strategy based on the Cre-loxP recombination system to generate large chromosomal rearrangements in Lactococcus lactis. Two loxP sites were sequentially integrated in inverse order into the chromosome either at random locations by transposition or at fixed points by homologous recombination. The recombination between the two chromosomal loxP sites was highly efficient (approximately 1 × 10⁻¹/cell) when the Cre recombinase was provided in trans, and parental- or inverted-type chromosomal structures were isolated after removal of the Cre recombinase. The usefulness of this approach was demonstrated by creating three large inversions of 500, 1,115, and 1,160 kb in size that modified the lactococcal genome organization to different extents. The Cre-loxP recombination system described can potentially be used for other gram-positive bacteria without further modification.     

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time‐consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome–environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described “islands of differentiation,” and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.     

Lack of correlation between crossing-over and chromosome distance between inversions in Drosophila ananassae

Two linked inversions, AL and ZE, located in the opposite limbs of the second chromosome of Drosophila ananassae are separated from each other by nearly 32% of the total length of the second chromosome. Crossing-over between these inversions when heterozygous was studied in females and males by the salivary-gland smear technique using karyotypically homozygous stocks. The results of recombination experiments show that there is a strong suppression of recombination between inversions when heterozygous, in spite of a large euchromatic distance available for crossing-over between them. Thus there is no correlation between chromosome distance and crossing-over between heterozygous inversions in the second chromosome of D. ananassae when studied cytologically.