Fine-Scale Position Effects Shape the Distribution of Inversion Breakpoints in Drosophila melanogaster

The potential of inversions to accumulate balanced sexual antagonism is supported by simulations and Drosophila experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with In(3 R)K favoring survival and In(3 L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, In(3 L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas In(3 R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Inversions Can Accumulate Balanced Sexual Antagonism: Evidence from Simulations and Drosophila Experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with ln(3R)K favoring survival and ln(3L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, ln(3L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas ln(3R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Genomics of natural populations: gene conversion events reveal selected genes within the inversions of Drosophila pseudoobscura

When adaptive phenotypic variation or quantitative trait loci map within an inverted segment of a chromosome, researchers often despair because the suppression of crossing over will prevent the discovery of selective target genes that established the rearrangement. If an inversion polymorphism is old enough, then the accumulation of gene conversion tracts offers the promise that quantitative trait loci or selected loci within inversions can be mapped. The inversion polymorphism of Drosophila pseudoobscura is a model system to show that gene conversion analysis is a useful tool for mapping selected loci within inversions. D. pseudoobscura has over 30 different chromosomal arrangements on the third chromosome (Muller C) in natural populations and their frequencies vary with changes in environmental habitats. Statistical tests of five D. pseudoobscura gene arrangements identified outlier genes within inverted regions that had potentially heritable variation, either fixed amino acid differences or differential expression patterns. We use genome sequences of the inverted third chromosome (Muller C) to infer 98,443 gene conversion tracts for a total coverage of 142 Mb or 7.2× coverage of the 19.7 Mb chromosome. We estimated gene conversion tract coverage in the 2,668 genes on Muller C and tested whether gene conversion coverage was similar among arrangements for outlier vs non-outlier loci. Outlier genes had lower gene conversion tract coverage among arrangements than the non-outlier genes suggesting that selection removes exchanged DNA in the outlier genes. These data support the hypothesis that the third chromosome in D. pseudoobscura captured locally adapted combinations of alleles prior to inversion mutation events.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

The Effect of an Irradiation-Induced Recombination Suppressing Inversion on the Genetic Stability and Biological Quality of a White Eye-Based Aedes aegypti Genetic Sexing Strain

Aedes aegypti is the primary vector of diseases such as dengue, chikungunya, Zika fever, and yellow fever. The sterile insect technique (SIT) has been proposed as a species-specific and environment-friendly tool for the suppression of mosquito vector populations as a major component of integrated vector management strategies. As female mosquitoes are blood-feeders and may transmit pathogenic microorganisms, mosquito SIT depends on the release of sterile males. Genetic sexing strains (GSS) can be used for the efficient and robust separation of males from females. Two Ae. aegypti GSS were recently developed by exploiting eye colour mutations, resulting in the Red-eye GSS (RGSS) and the White-eye GSS (WGSS). In this study, we compared two WGSS, with and without the chromosomal inversion 35 (Inv35), and evaluated their biological quality, including genetic stability. Our results suggest that the WGSS/Inv35 presents a low recombination rate and long-term genetic stability when recombinants are removed from the colony (filtering) and a slow accumulation of recombinants when they are not removed from the colony (non-filtering). The two strains were similar with respect to fecundity, pupal and adult recovery rates, pupation curve, and pupal weight. However, differences were detected in fertility, survival rate of females, and flight ability of males. The WGSS/Inv35 presented lower fertility, higher survival rate of females, and better flight ability of males compared to the WGSS.

The relevance of chromatin architecture to genome rearrangements in Drosophila

DNA within chromosomes in the nucleus is non-randomly organized into chromosome territories, compartments and topologically associated domains (TADs). Chromosomal rearrangements have the potential to alter chromatin organization and modify gene expression leading to selection against these structural variants. Drosophila pseudoobscura has a wealth of naturally occurring gene arrangements that were generated by overlapping inversion mutations caused by two chromosomal breaks that rejoin the central region in reverse order. Unlike humans, Drosophila inversion heterozygotes do not have negative effects associated with crossing over during meiosis because males use achiasmate mechanisms for proper segregation, and aberrant recombinant meiotic products generated in females are lost in polar bodies. As a result, Drosophila populations are found to harbour extensive inversion polymorphisms. It is not clear, however, whether chromatin architecture constrains which inversions breakpoints persist in populations. We mapped the breakpoints of seven inversions in D. pseudoobscura to the TAD map to determine if persisting inversion breakpoints are more likely to occur at boundaries between TADs. Our results show that breakpoints occur at TAD boundaries more than expected by chance. Some breakpoints may alter gene expression within TADs supporting the hypothesis that position effects contribute to inversion establishment. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Comparative Genomics and Evolution of Avian Specialized Traits

Genomic data are important for understanding the origin and evolution of traits. Under the context of rapidly developing of sequencing technologies and more widely available genome sequences, researchers are able to study evolutionary mechanisms of traits via comparative genomic methods. Compared with other vertebrates, bird genomes are relatively small and exhibit conserved synteny with few repetitive elements, which makes them suitable for evolutionary studies. Increasing genomic progress has been reported on the evolution of powered flight, body size variation, beak morphology, plumage colouration, high-elevation colonization, migration, and vocalization. By summarizing previous studies, we demonstrate the genetic bases of trait evolution, highlighting the roles of small-scale sequence variation, genomic structural variation, and changes in gene interaction networks. We suggest that future studies should focus on improving the quality of reference genomes, exploring the evolution of regulatory elements and networks, and combining genomic data with morphological, ecological, behavioural, and developmental biology data.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

Inversion breakpoints and the evolution of supergenes

The coexistence of discrete morphs that differ in multiple traits is common within natural populations of many taxa. Such morphs are often associated with chromosomal inversions, presumably because the recombination suppressing effects of inversions help maintain alternate adaptive combinations of alleles across the multiple loci affecting these traits. However, inversions can also harbour selected mutations at their breakpoints, leading to their rise in frequency in addition to (or independent from) their role in recombination suppression. In this review, we first describe the different ways that breakpoints can create mutations. We then critically examine the evidence for the breakpoint-mutation and recombination suppression hypotheses for explaining the existence of discrete morphs associated with chromosomal inversions. We find that the evidence that inversions are favoured due to recombination suppression is often indirect. The evidence that breakpoints harbour mutations that are adaptive is also largely indirect, with the characterization of inversion breakpoints at the sequence level being incomplete in most systems. Direct tests of the role of suppressed recombination and breakpoint mutations in inversion evolution are thus needed. Finally, we emphasize how the two hypotheses of recombination suppression and breakpoint mutation can act in conjunction, with implications for understanding the emergence of supergenes and their evolutionary dynamics. We conclude by discussing how breakpoint characterization could improve our understanding of complex, discrete phenotypic forms in nature.

The Cyclically Seasonal Drosophila subobscura Inversion O7 Originated From Fragile Genomic Sites and Relocated Immunity and Metabolic Genes

Chromosome inversions are important contributors to standing genetic variation in Drosophila subobscura. Presently, the species is experiencing a rapid replacement of high-latitude by low-latitude inversions associated with global warming. Yet not all low-latitude inversions are correlated with the ongoing warming trend. This is particularly unexpected in the case of O7 because it shows a regular seasonal cycle that peaks in summer and rose with a heatwave. The inconsistent behavior of O7 across components of the ambient temperature suggests that is causally more complex than simply due to temperature alone. In order to understand the dynamics of O7, high-quality genomic data are needed to determine both the breakpoints and the genetic content. To fill this gap, here we generated a PacBio long read-based chromosome-scale genome assembly, from a highly homozygous line made isogenic for an O3 + 4 + 7 chromosome. Then we isolated the complete continuous sequence of O7 by conserved synteny analysis with the available reference genome. Main findings include the following: (i) the assembled O7 inversion stretches 9.936 Mb, containing > 1,000 annotated genes; (ii) O7 had a complex origin, involving multiple breaks associated with non-B DNA-forming motifs, formation of a microinversion, and ectopic repair in trans with the two homologous chromosomes; (iii) the O7 breakpoints carry a pre-inversion record of fragility, including a sequence insertion, and transposition with later inverted duplication of an Attacin immunity gene; and (iv) the O7 inversion relocated the major insulin signaling forkhead box subgroup O (foxo) gene in tight linkage with its antagonistic regulatory partner serine/threonine-protein kinase B (Akt1) and disrupted concerted evolution of the two inverted Attacin duplicates, reattaching them to dFOXO metabolic enhancers. Our findings suggest that O7 exerts antagonistic pleiotropic effects on reproduction and immunity, setting a framework to understand its relationship with climate change. Furthermore, they are relevant for fragility in genome rearrangement evolution and for current views on the contribution of breakage versus repair in shaping inversion-breakpoint junctions.