Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders.

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

The impact of inversions across 33,924 families with rare disease from a national genome sequencing project.

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

The barley pan-genome reveals the hidden legacy of mutation breeding.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Complex evolutionary processes maintain an ancient chromosomal inversion.

Genome re-arrangements such as chromosomal inversions are often involved in adaptation. As such, they experience natural selection, which can erode genetic variation. Thus, whether and how inversions can remain polymorphic for extended periods of time remains debated. Here we combine genomics, experiments, and evolutionary modeling to elucidate the processes maintaining an inversion polymorphism associated with the use of a challenging host plant (Redwood trees) in Timema stick insects. We show that the inversion is maintained by a combination of processes, finding roles for life-history trade-offs, heterozygote advantage, local adaptation to different hosts, and gene flow. We use models to show how such multi-layered regimes of balancing selection and gene flow provide resilience to help buffer populations against the loss of genetic variation, maintaining the potential for future evolution. We further show that the inversion polymorphism has persisted for millions of years and is not a result of recent introgression. We thus find that rather than being a nuisance, the complex interplay of evolutionary processes provides a mechanism for the long-term maintenance of genetic variation.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes.

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

A high-resolution map of small-scale inversions in the gibbon genome.

Gibbons are the most speciose family of living apes, characterized by a diverse chromosome number and rapid rate of large-scale rearrangements. Here we performed single-cell template strand sequencing (Strand-seq), molecular cytogenetics, and deep in silico analysis of a southern white-cheeked gibbon genome, providing the first comprehensive map of 238 previously hidden small-scale inversions. We determined that more than half are gibbon specific, at least fivefold higher than shown for other primate lineage-specific inversions, with a significantly high number of small heterozygous inversions, suggesting that accelerated evolution of inversions may have played a role in the high sympatric diversity of gibbons. Although the precise mechanisms underlying these inversions are not yet understood, it is clear that segmental duplication-mediated NAHR only accounts for a small fraction of events. Several genomic features, including gene density and repeat (e.g., LINE-1) content, might render these regions more break-prone and susceptible to inversion formation. In the attempt to characterize interspecific variation between southern and northern white-cheeked gibbons, we identify several large assembly errors in the current GGSC Nleu3.0/nomLeu3 reference genome comprising more than 49 megabases of DNA. Finally, we provide a list of 182 candidate genes potentially involved in gibbon diversification and speciation.

Systematic identification of gene-altering programmed inversions across the bacterial domain.

Programmed chromosomal inversions allow bacteria to generate intra-population genotypic and functional heterogeneity, a bet-hedging strategy important in changing environments. Some programmed inversions modify coding sequences, producing different alleles in several gene families, most notably in specificity-determining genes such as Type I restriction-modification systems, where systematic searches revealed cross phylum abundance. Yet, a broad, gene-independent, systematic search for gene-altering programmed inversions has been absent, and little is known about their genomic sequence attributes and prevalence across gene families. Here, identifying intra-species variation in genomes of over 35 000 species, we develop a predictive model of gene-altering inversions, revealing key attributes of their genomic sequence attributes, including gene-pseudogene size asymmetry and orientation bias. The model predicted over 11,000 gene-altering loci covering known targeted gene families, as well as novel targeted families including Type II restriction-modification systems, a protein of unknown function, and a fusion-protein containing conjugative-pilus and phage tail domains. Publicly available long-read sequencing datasets validated representatives of these newly predicted inversion-targeted gene families, confirming intra-population genetic heterogeneity. Together, these results reveal gene-altering programmed inversions as a key strategy adopted across the bacterial domain, and highlight programmed inversions that modify Type II restriction-modification systems as a possible new mechanism for maintaining intra-population heterogeneity.

Conservation of a Chromosome 8 Inversion and Exon Mutations Confirm Common Gulonolactone Oxidase Gene Evolution Among Primates, Including H. Neanderthalensis.

Ascorbic acid functions as an antioxidant and facilitates other biochemical processes such as collagen triple helix formation, and iron uptake by cells. Animals which endogenously produce ascorbic acid have a functional gulonolactone oxidase gene (GULO); however, humans have a GULO pseudogene (GULOP) and depend on dietary ascorbic acid. In this study, the conservation of GULOP sequences in the primate haplorhini suborder were investigated and compared to the GULO sequences belonging to the primates strepsirrhini suborder. Phylogenetic analysis suggested that the conserved GULOP exons in the haplorhini primates experienced a high rate of mutations following the haplorhini/strepsirrhini divergence. This high mutation rate has decreased during the evolution of the haplorhini primates. Additionally, indels of the haplorhini GULOP sequences were conserved across the suborder. A separate analysis for GULO sequences and well-conserved GULOP sequences focusing on placental mammals identified an in-frame GULO sequence in the Brazilian guinea pig, and a potential GULOP sequence in the pika. Similar to haplorhini primates, the guinea pig and lagomorph species have experienced a high substitution rate when compared to the mammals used in this study. A shared synteny to examine the conservation of local genes near GULO/GULOP identified a conserved inversion around the GULO/GULOP locus between the haplorhini and strepsirrhini primates. Fischer's exact test did not support an association between GULOP and the chromosomal inversion. Mauve alignment showed that the inversion of the length of the syntenic block that the GULO/GULOP genes belonged to was variable. However, there were frequent rearrangements around ~ 2 million base pairs adjacent to GULOP involving the KIF13B and MSRA genes. These data may suggest that genes acquiring deleterious mutations in the coding sequence may respond to these deleterious mutations with rapid substitution rates.

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding.

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

Nanopore long-read sequencing analysis reveals ZIC1 dysregulation caused by a de novo 3q inversion with a breakpoint located 7 kb downstream of ZIC1.

Zic family member 1 (ZIC1), a gene located on chromosome 3q24, encodes a transcription factor with zinc finger domains that is essential for the normal development of the cerebellum. Heterozygous loss-of-function of ZIC1 causes Dandy-Walker malformation, while heterozygous gain-of-function leads to a multiple congenital anomaly syndrome characterized by craniosynostosis, brain abnormalities, facial features, and learning disability. In this study, we present the results of genetic analysis of a male patient with clinically suspected Gomez-Lopez-Hernandez syndrome. The patient displayed multiple congenital abnormalities, including bicoronal craniosynostosis, characteristic facial features, cerebellar malformation with rhombencephalosynapsis, and temporal alopecia, and a de novo inversion of chromosome 3q. Breakpoint analysis using a Nanopore long-read sequencer revealed a breakpoint in the distal centromere of 3q24 located 7 kb downstream of the 3' untranslated region of ZIC1. On the basis of the clinical similarities, we concluded that the abnormalities in this patient were caused by the transcriptional dysregulation of ZIC1. We hypothesize the underlying molecular mechanisms of transcriptional dysregulation of ZIC1 such as the abnormalities in topologically associated domains encompassing ZIC1. This study highlights the usefulness of long-read sequencing in the analysis of de novo balanced chromosomal abnormalities.

Genome sequencing detects a balanced pericentric inversion with breakpoints that impact the DMD and upstream region of POU3F4 genes.

We describe a family with two maternal half-brothers both of whom presented with muscular dystrophy, autism spectrum disorder, developmental delay, and sensorineural hearing loss. The elder brother had onset of features at ~3 months of age, followed by clinical confirmation of muscular dystrophy at 3 years. Skeletal biopsy staining at 4.7 years showed an absence of dystrophin protein which prompted extensive molecular testing over 4 years that included gene panels, targeted single-gene assays, arrays, and karyotyping, all of which failed to identify a clinically significant variant in the DMD gene. At 10 years of age, clinical whole-genome sequencing (cWGS) was performed, which revealed a novel hemizygous ~50.7 Mb balanced pericentric inversion on chromosome X that disrupts the DMD gene in both siblings, consistent with the muscular dystrophy phenotype. This inversion also impacts the upstream regulatory region of POU3F4, structural rearrangements which are known to cause hearing loss. The unaffected mother is a heterozygous carrier for the pericentric inversion. This finding illustrates the ability of cWGS to detect a wide breadth of disease-causing genomic variations including large genomic rearrangements.

Conclusion of diagnostic odysseys due to inversions disrupting GLI3 and FBN1.

Many genetic testing methodologies are biased towards picking up structural variants (SVs) that alter copy number. Copy-neutral rearrangements such as inversions are therefore likely to suffer from underascertainment. In this study, manual review prompted by a virtual multidisciplinary team meeting and subsequent bioinformatic prioritisation of data from the 100K Genomes Project was performed across 43 genes linked to well-characterised skeletal disorders. Ten individuals from three independent families were found to harbour diagnostic inversions. In two families, inverted segments of 1.2/14.8 Mb unequivocally disrupted GLI3 and segregated with skeletal features consistent with Greig cephalopolysyndactyly syndrome. For one family, phenotypic blending was due to the opposing breakpoint lying ~45 kb from HOXA13 In the third family, long suspected to have Marfan syndrome, a 2.0 Mb inversion disrupting FBN1 was identified. These findings resolved lengthy diagnostic odysseys of 9-20 years and highlight the importance of direct interaction between clinicians and data-analysts. These exemplars of a rare mutational class inform future SV prioritisation strategies within the NHS Genomic Medicine Service and similar genome sequencing initiatives. In over 30 years since these two disease-gene associations were identified, large inversions have yet to be described and so our results extend the mutational spectra linked to these conditions.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species.

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present "haplotagging," a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline.

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

Tempo and mode in karyotype evolution revealed by a probabilistic model incorporating both chromosome number and morphology.

Karyotype, including the chromosome and arm numbers, is a fundamental genetic characteristic of all organisms and has long been used as a species-diagnostic character. Additionally, karyotype evolution plays an important role in divergent adaptation and speciation. Centric fusion and fission change chromosome numbers, whereas the intra-chromosomal movement of the centromere, such as pericentric inversion, changes arm numbers. A probabilistic model simultaneously incorporating both chromosome and arm numbers has not been established. Here, we built a probabilistic model of karyotype evolution based on the "karyograph", which treats karyotype evolution as a walk on the two-dimensional space representing the chromosome and arm numbers. This model enables analysis of the stationary distribution with a stable karyotype for any given parameter. After evaluating their performance using simulated data, we applied our model to two large taxonomic groups of fish, Eurypterygii and series Otophysi, to perform maximum likelihood estimation of the transition rates and reconstruct the evolutionary history of karyotypes. The two taxa significantly differed in the evolution of arm number. The inclusion of speciation and extinction rates demonstrated possibly high extinction rates in species with karyotypes other than the most typical karyotype in both groups. Finally, we made a model including polyploidization rates and applied it to a small plant group. Thus, the use of this probabilistic model can contribute to a better understanding of tempo and mode in karyotype evolution and its possible role in speciation and extinction.

Polymorphic Inversions Underlie the Shared Genetic Susceptibility of Obesity-Related Diseases.

The burden of several common diseases including obesity, diabetes, hypertension, asthma, and depression is increasing in most world populations. However, the mechanisms underlying the numerous epidemiological and genetic correlations among these disorders remain largely unknown. We investigated whether common polymorphic inversions underlie the shared genetic influence of these disorders. We performed an inversion association analysis including 21 inversions and 25 obesity-related traits on a total of 408,898 Europeans and validated the results in 67,299 independent individuals. Seven inversions were associated with multiple diseases while inversions at 8p23.1, 16p11.2, and 11q13.2 were strongly associated with the co-occurrence of obesity with other common diseases. Transcriptome analysis across numerous tissues revealed strong candidate genes for obesity-related traits. Analyses in human pancreatic islets indicated the potential mechanism of inversions in the susceptibility of diabetes by disrupting the cis-regulatory effect of SNPs from their target genes. Our data underscore the role of inversions as major genetic contributors to the joint susceptibility to common complex diseases.

The molecular mechanisms of recombinant chromosome 18 with parental pericentric inversions and a review of the literature.

Chromosomal rearrangements mostly result from non-allelic homologous recombination mediated by low-copy repeats (LCRs) or segmental duplications (SDs). Recent studies on recombinant chromosome 18 (rec (18)) have focused on diagnoses and clinical phenotypes. We diagnosed two cases of prenatal rec (18) and identified precise breakpoint intervals using karyotype and chromosomal microarray analyses. We analyzed the distribution characteristics of breakpoint repetitive elements to infer rearrangement mechanisms and reviewed relevant literature to identify genetic trends. Among the 12 families with 25 pregnancies analyzed, 68% rec (18), 24% spontaneous abortions, and 8% normal births were reported. In the 17 rec (18) cases, 65% presented maternal origin and 35% were paternal. Short-arm breakpoints at p11.31 were reported in 10 cases, whereas the long-arm breakpoints were located at q21.3 (6 cases) and q12 (4 cases). Breakpoints of pericentric inversions on chromosome 18 are concentrated in p11.31, q21.3, and q12 regions. Rearrangements at 18p11.31 are non-recurrent events. ALUs, LINE1s, and MIRs were enriched at the breakpoint regions (1.85 to 3.42-fold enrichment over the entire chromosome 18), while SDs and LCRs were absent. ALU subfamilies had sequence identities of 85.94% and 83.01% between two pair breakpoints. Small repetitive elements may mediate recombination-coupled DNA repair processes, facilitating rearrangements on chromosome 18. Maternal inversion carriers are more prone to abnormal recombination in prenatal families with rec (18). Recombinant chromosomes may present preferential segregation during gamete formation.

Experimental introgression in Drosophila: Asymmetric postzygotic isolation associated with chromosomal inversions and an incompatibility locus on the X chromosome.

Interspecific gene flow (introgression) is an important source of new genetic variation, but selection against it can reinforce reproductive barriers between interbreeding species. We used an experimental approach to trace the role of chromosomal inversions and incompatibility genes in preventing introgression between two partly sympatric Drosophila virilis group species, D. flavomontana and D. montana. We backcrossed F1 hybrid females from a cross between D. flavomontana female and D. montana male with the males of the parental species for two generations and sequenced pools of parental strains and their reciprocal second generation backcross (BC2 mon and BC2 fla) females. Contrasting the observed amount of introgression (mean hybrid index, HI) in BC2 female pools along the genome to simulations under different scenarios allowed us to identify chromosomal regions of restricted and increased introgression. We found no deviation from the HI expected under a neutral null model for any chromosome for the BC2 mon pool, suggesting no evidence for genetic incompatibilities in backcrosses towards D. montana. In contrast, the BC2 fla pool showed high variation in the observed HI between different chromosomes, and massive reduction of introgression on the X chromosome (large X-effect). This observation is compatible with reduced recombination combined with at least one dominant incompatibility locus residing within the X inversion(s). Overall, our study suggests that genetic incompatibilities arising within chromosomal inversions can play an important role in speciation.

Weak antagonistic fitness effects can maintain an inversion polymorphism.

The study of chromosomal inversion polymorphisms has received much recent attention, particularly in cases where inversions have drastic effects on phenotypes and fitness (e.g. lethality of homozygotes). Less attention has been paid to the question of the maintenance of inversion polymorphisms that show only weak effects. Here, we study the maintenance of such an inversion polymorphism that links 250 genes on chromosome Tgu11 in the zebra finch (Taeniopygia guttata). Based on data from over 6000 captive birds, we estimated the effects of this inversion on a wide range of fitness-related traits. We found that, compared with the ancestral allele A, the inverted allele D had small additive beneficial effects on male siring success and on female fecundity. These fitness-enhancing effects may explain the initial spread of the derived D allele (allele frequency 53%). However, individuals that were homozygous for D had a slightly lower survival rate, which may explain why the D allele has not spread to fixation. We used individual-based simulations to examine how an inversion polymorphism with such antagonistic fitness effects behaves over time. Our results indicate that polymorphisms become stabilized at an intermediate allele frequency if the inversion links an additively beneficial allele of small effect size to a recessive weakly deleterious mutation, overall resulting in weak net heterosis. Importantly, this conclusion remains valid over a wide range of selection coefficients against the homozygous DD (up to lethality), suggesting that the conditions needed to maintain the polymorphism may frequently be met. However, the simulations also suggest that in our zebra finch populations, the estimated recessive deleterious effect of the D allele (on survival in captivity) is not quite large enough to prevent fixation of the D allele in the long run. Estimates of fitness effects from free-living populations are needed to validate these results.