Genetic changes in flowering and morphology in response to adaptation to a high-latitude environment in Arabidopsis lyrata

Local climate and vernalization sensitivity predict the latitudinal patterns of flowering onset in the crop wild relative Linum bienne Mill

BACKGROUND AND AIMS

The timing of flowering onset is often correlated with latitude, indicative of climatic gradients. Flowering onset in temperate species commonly requires exposure to cold temperatures, known as vernalization. Hence, population differentiation of flowering onset with latitude might reflect adaptation to the local climatic conditions experienced by populations.

METHODS

Within its western range, seeds from Linum bienne populations (the wild relative of cultivated Linum usitatissimum) were used to describe the latitudinal differentiation of flowering onset to determine its association with the local climate of the population. A vernalization experiment including different crop cultivars was used to determine how vernalization accelerates flowering onset, in addition to the vernalization sensitivity response among populations and cultivars. Additionally, genetic differentiation of L. bienne populations along the latitudinal range was scrutinized using microsatellite markers.

KEY RESULTS

Flowering onset varied with latitude of origin, with southern populations flowering earlier than their northern counterparts. Vernalization reduced the number of days to flowering onset, but vernalization sensitivity was greater in northern populations compared with southern ones. Conversely, vernalization delayed flowering onset in the crop, exhibiting less variation in sensitivity. In L. bienne, both flowering onset and vernalization sensitivity were better predicted by the local climate of the population than by latitude itself. Microsatellite data unveiled genetic differentiation of populations, forming two groups geographically partitioned along latitude.

CONCLUSIONS

The consistent finding of latitudinal variation across experiments suggests that both flowering onset and vernalization sensitivity in L. bienne populations are under genetic regulation and might depend on climatic cues at the place of origin. The association with climatic gradients along latitude suggests that the climate experienced locally drives population differentiation of the flowering onset and vernalization sensitivity patterns. The genetic population structure suggests that past population history could have influenced the flowering initiation patterns detected, which deserves further work.

The role of introgressive hybridization in shaping the geographically isolated gene pools of wax palm populations (genus Ceroxylon)

The speciation continuum is the process by which genetic groups diverge until they reach reproductive isolation. It has become common in the literature to show that this process is gradual and flickering, with possibly many instances of secondary contact and introgression after divergence has started. The level of divergence might vary among genomic regions due to, among others, the different forces and roles of selection played by the shared regions. Through hybrid capture, we sequenced ca. 4,000 nuclear regions in populations of six species of wax palms, five of which form a monophyletic group (genus Ceroxylon, Arecaceae: Ceroxyloideae). We show that in this group, the different populations show varying degrees of introgressive hybridization, and two of them are backcrosses of the other three 'pure' species. This is particularly interesting because these three species are dioecious, have a shared main pollinator, and have slightly overlapping reproductive seasons but highly divergent morphologies. Our work supports shows wax palms diverge under positive and background selection in allopatry, and hybridize due to secondary contact and inefficient reproductive barriers, which sustain genetic diversity. Introgressed regions are generally not under positive selection. Peripheral populations are backcrosses of other species; thus, introgressive hybridization is likely modulated by demographic effects rather than selective pressures. In general, these species might function as an 'evolutionary syngameon' where expanding, peripheral, small, and isolated populations maintain diversity by crossing with available individuals of other wax palms. In the Andean context, species can benefit from gained variation from a second taxon or the enhancement of population sizes by recreating a common genetic pool.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

Maintenance of Adaptive Dynamics and No Detectable Load in a Range-Edge Outcrossing Plant Population

During range expansion, edge populations are expected to face increased genetic drift, which in turn can alter and potentially compromise adaptive dynamics, preventing the removal of deleterious mutations and slowing down adaptation. Here, we contrast populations of the European subspecies Arabidopsis lyrata ssp. petraea, which expanded its Northern range after the last glaciation. We document a sharp decline in effective population size in the range-edge population and observe that nonsynonymous variants segregate at higher frequencies. We detect a 4.9% excess of derived nonsynonymous variants per individual in the range-edge population, suggesting an increase of the genomic burden of deleterious mutations. Inference of the fitness effects of mutations and modeling of allele frequencies under the explicit demographic history of each population predicts a depletion of rare deleterious variants in the range-edge population, but an enrichment for fixed ones, consistent with the bottleneck effect. However, the demographic history of the range-edge population predicts a small net decrease in per-individual fitness. Consistent with this prediction, the range-edge population is not impaired in its growth and survival measured in a common garden experiment. We further observe that the allelic diversity at the self-incompatibility locus, which ensures strict outcrossing and evolves under negative frequency-dependent selection, has remained unchanged. Genomic footprints indicative of selective sweeps are broader in the Northern population but not less frequent. We conclude that the outcrossing species A. lyrata ssp. petraea shows a strong resilience to the effect of range expansion.

Two reproductive traits show contrasting genetic architectures in Plantago lanceolata

In many species, temperature-sensitive phenotypic plasticity (i.e., an individual's phenotypic response to temperature) displays a positive correlation with latitude, a pattern presumed to reflect local adaptation. This geographical pattern raises two general questions: (a) Do a few large-effect genes contribute to latitudinal variation in a trait? (b) Is the thermal plasticity of different traits regulated pleiotropically? To address the questions, we crossed individuals of Plantago lanceolata derived from northern and southern European populations. Individuals naturally exhibited high and low thermal plasticity in floral reflectance and flowering time. We grew parents and offspring in controlled cool- and warm-temperature environments, mimicking what plants would encounter in nature. We obtained genetic markers via genotype-by-sequencing, produced the first recombination map for this ecologically important nonmodel species, and performed quantitative trait locus (QTL) mapping of thermal plasticity and single-environment values for both traits. We identified a large-effect QTL that largely explained the reflectance plasticity differences between northern and southern populations. We identified multiple smaller-effect QTLs affecting aspects of flowering time, one of which affected flowering time plasticity. The results indicate that the genetic architecture of thermal plasticity in flowering is more complex than for reflectance. One flowering time QTL showed strong cytonuclear interactions under cool temperatures. Reflectance and flowering plasticity QTLs did not colocalize, suggesting little pleiotropic genetic control and freedom for independent trait evolution. Such genetic information about the architecture of plasticity is environmentally important because it informs us about the potential for plasticity to offset negative effects of climate change.

Inflorescence shoot elongation, but not flower primordia formation, is photoperiodically regulated in Arabidopsis lyrata

BACKGROUND AND AIMS

Photoperiod contains information about the progress of seasons. Plants use the changing photoperiod as a cue for the correct timing of important life history events, including flowering. Here the effect of photoperiod on flowering in four Arabidopsis lyrata populations originating from different latitudes was studied, as well as expression levels of candidate genes for governing the between-population differences.

METHODS

Flowering of plants from four A. lyrata populations was studied in three different photoperiods after vernalization. Flowering development was separated into three steps: flower primordia formation, inflorescence shoot elongation and opening of the first flower. Circadian expression rhythms of the A. lyrata homologues of GIGANTEA (GI), FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (FKF1), CONSTANS (CO) and FLOWERING LOCUS T (FT) were studied in three of the populations in the intermediate (14 h) photoperiod treatment.

KEY RESULTS

Most plants in all populations formed visible flower primordia during vernalization. Further inflorescence development after vernalization was strongly inhibited by short days in the northern European population (latitude 61°N), only slightly in the central European population (49°N) and not at all in the North American populations (36°N and 42°N). In the 14 h daylength, where all plants from the three southernmost populations but only 60 % of the northernmost population flowered, the circadian expression rhythm of the A. lyrata FT was only detected in the southern populations, suggesting differentiation in the critical daylength for activation of the long-day pathway. However, circadian expression rhythms of A. lyrata GI, FKF1 and CO were similar between populations.

CONCLUSIONS

The results indicate that in A. lyrata, transition to flowering can occur through pathways independent of long days, but elongation of inflorescences is photoperiodically regulated.

Local adaptation and ecological differentiation under selection, migration, and drift in Arabidopsis lyrata

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole-genome-based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4Ne m < 11) and have small effective population sizes (Ne <  6000), suggesting that both migration and drift can counteract local selection. However, the three-year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.

Low genetic diversity contrasts with high phenotypic variability in heptaploid Spartina densiflora populations invading the Pacific coast of North America

Species can respond to environmental pressures through genetic and epigenetic changes and through phenotypic plasticity, but few studies have evaluated the relationships between genetic differentiation and phenotypic plasticity of plant species along changing environmental conditions throughout wide latitudinal ranges. We studied inter‐ and intrapopulation genetic diversity (using simple sequence repeats and chloroplast DNA sequencing) and inter‐ and intrapopulation phenotypic variability of 33 plant traits (using field and common‐garden measurements) for five populations of the invasive cordgrass Spartina densiflora Brongn. along the Pacific coast of North America from San Francisco Bay to Vancouver Island. Studied populations showed very low genetic diversity, high levels of phenotypic variability when growing in contrasted environments and high intrapopulation phenotypic variability for many plant traits. This intrapopulation phenotypic variability was especially high, irrespective of environmental conditions, for those traits showing also high phenotypic plasticity. Within‐population variation represented 84% of the total genetic variation coinciding with certain individual plants keeping consistent responses for three plant traits (chlorophyll b and carotenoid contents, and dead shoot biomass) in the field and in common‐garden conditions. These populations have most likely undergone genetic bottleneck since their introduction from South America; multiple introductions are unknown but possible as the population from Vancouver Island was the most recent and one of the most genetically diverse. S. densiflora appears as a species that would not be very affected itself by climate change and sea‐level rise as it can disperse, establish, and acclimate to contrasted environments along wide latitudinal ranges.

Trait variation in response to varying winter temperatures, diversity patterns and signatures of selection along the latitudinal distribution of the widespread grassland plant Arrhenatherum elatius

Across Europe, genetic diversity can be expected to decline toward the North because of stochastic and selective effects which may imply diminished phenotypic variation and less potential for future genetic adaptations to environmental change. Understanding such latitudinal patterns can aid provenance selection for breeding or assisted migration approaches. In an experiment simulating different winter temperatures, we assessed quantitative trait variation, genetic diversity, and differentiation for natural populations of the grass Arrhenatherum elatius originating from a large latitudinal gradient. In general, populations from the North grew smaller and had a lower flowering probability. Toward the North, the absolute plastic response to the different winter conditions as well as heritability for biomass production significantly declined. Genetic differentiation in plant height and probability of flowering were very strong and significantly higher than under neutral expectations derived from SNP data, suggesting adaptive differentiation. Differentiation in biomass production did not exceed but mirrored patterns for neutral genetic differentiation, suggesting that migration‐related processes caused the observed clinal trait variation. Our results demonstrate that genetic diversity and trait differentiation patterns for A. elatius along a latitudinal gradient are likely shaped by both local selection and genetic drift.

Role of seed germination in adaptation and reproductive isolation in Arabidopsis lyrata

Seed germination is an important developmental and life history stage. Yet, the evolutionary impact of germination has mainly been studied in the context of dormancy, or for its role in reproductive isolation between species. Here, we aim to examine multiple consequences of genetic divergence on germination traits between two Arabidopsis lyrata subspecies: ssp. petraea (Eurasia) and ssp. lyrata (North America). Postdormancy germination time, a potentially adaptive trait, showed differentiation between the populations, and quantitative trait loci (QTL) mapping revealed that the trait variation is mainly controlled by two antagonistic loci. These QTL areas contain several candidate genes with known function in postdormancy germination in A. thaliana. The sequence variation of three genes was consistent with differential selection, and they also included fixed nonsynonymous substitutions with potential to account for the phenotypic differentiation. We further show that the divergence between the subspecies has led to a slight but significant reduction in hybrid germination proportions, indicating incipient reproductive isolation. Comparison of reciprocal F₁ and F₂ progenies suggests that Bateson-Dobzhansky-Muller incompatibilities likely act through uniparentally inherited factors. Examination of genomewide transmission ratio distortion further revealed that cytonuclear interactions cause substantial pregermination inviability in the hybrids. These results confirm that seed germination has adaptive potential beyond the dormancy stage and that hybrid seed inviability can be one of the first reproductive barriers to arise during divergence.

Timing of shoot development transitions affects degree of perenniality in Arabidopsis lyrata (Brassicaceae)

BACKGROUND

Perenniality is best understood in quantitative terms, involving the relationship between production vs. turnover of meristems, biomass, or energy reserves. Previous quantitative trait locus (QTL) studies using divergent populations of the perennial rock cress Arabidopsis lyrata have shown that trade-offs in vegetative growth vs. reproduction are due to cascading effects of differences in early vegetative development, which contribute to local adaptation. However, details of the developmental differences and how they affect perenniality remained unclear. In this study, we investigated in detail the developmental differences in perenniality between populations. A. lyrata from Norway and North Carolina populations, representing contrasting environments and degrees of perenniality, were grown under controlled conditions, and data were collected on plant phenology and shoot-level development. We tested hypotheses that differences in perenniality involve strict allocation of lateral meristems to vegetative vs. reproductive fates, or alternatively quantitative effects of pre-reproductive vegetative development.

RESULTS

The two populations showed large differences in the degree of vegetative development on individual shoots prior to reproductive transitions. The number of leaves produced on shoots prior to bolting, and not strict meristem allocation or variation in apical dominance, was able to explain variation in the number of inflorescences on individual plants. These results suggested that allocation of time to shoot vegetative vs. reproductive development could be a major factor in resource allocation differences between the populations.

CONCLUSIONS

Based on these results and those of previous QTL studies, we propose a model in which the degree of shoot vegetative development shapes the developmental context for reproduction and subsequent vegetative growth in different environments. Climate-specific effects of shoot development patterns on reproductive output and survival may result in divergent evolutionary trajectories along a perenniality continuum, which may have broader relevance for plant life history evolution.

Genetic variation in the flowering and yield formation of timothy (Phleum pratense L.) accessions after different photoperiod and vernalization treatments

Timothy is a perennial forage grass grown commonly in Boreal regions. This study explored the effect of vernalization and photoperiod (PP) on flowering and growth characteristics and how this related to changes in expression of three flowering related genes in accessions from different geographic origin. Large variation was found in accessions in their vernalization and PP responses. In southern accessions vernalization response or requirement was not observed, the heading date remained unchanged, and plants flowered without vernalization. On the contrary, northern types had obligatory requirement for vernalization and long PP, but the tiller elongation did not require vernalization at 16-h PP. Longer vernalization or PP treatments reduced the genotypical differences in flowering. Moreover, the vernalization saturation progressed stepwise from main tiller to lateral tillers, and this process was more synchronized in southern accessions. The expression of PpVRN1 was associated with vernalization while PpVRN3 accumulated at long PP. A crucial role for PpVRN3 in the transition to flowering was supported as in southern accession the transcript accumulated in non-vernalized plants after transfer to 16-h PP, and the apices transformed to generative stage. Differences in vernalization requirements were associated with variation in expression levels of PpVRN1 and PpVRN3, with higher expression levels in southern type. Most divergent transcript accumulation of PpMADS10 was found under different vernalization conditions. These differences between accessions can be translated into agronomic traits, such as the tiller composition of canopy, which affects the forage yield. The southern types, with minimal vernalization response, have fast re-growth ability and rapidly decreasing nutritive value, whereas northern types grow slowly and have better quality. This information can be utilized in breeding for new cultivars for longer growing seasons at high latitudes.

Temperature-stress resistance and tolerance along a latitudinal cline in North American Arabidopsis lyrata

The study of latitudinal gradients can yield important insights into adaptation to temperature stress. Two strategies are available: resistance by limiting damage, or tolerance by reducing the fitness consequences of damage. Here we studied latitudinal variation in resistance and tolerance to frost and heat and tested the prediction of a trade-off between the two strategies and their costliness. We raised plants of replicate maternal seed families from eight populations of North American Arabidopsis lyrata collected along a latitudinal gradient in climate chambers and exposed them repeatedly to either frost or heat stress, while a set of control plants grew under standard conditions. When control plants reached maximum rosette size, leaf samples were exposed to frost and heat stress, and electrolyte leakage (PEL) was measured and treated as an estimate of resistance. Difference in maximum rosette size between stressed and control plants was used as an estimate of tolerance. Northern populations were more frost resistant, and less heat resistant and less heat tolerant, but-unexpectedly-they were also less frost tolerant. Negative genetic correlations between resistance and tolerance to the same and different thermal stress were generally not significant, indicating only weak trade-offs. However, tolerance to frost was consistently accompanied by small size under control conditions, which may explain the non-adaptive latitudinal pattern for frost tolerance. Our results suggest that adaptation to frost and heat is not constrained by trade-offs between them. But the cost of frost tolerance in terms of plant size reduction may be important for the limits of species distributions and climate niches.

Large-scale adaptive differentiation in the alpine perennial herb Arabis alpina

Information about the incidence and magnitude of local adaptation can help to predict the response of natural populations to a changing environment, and should be of particular interest in arctic and alpine environments where the effects of climate change are expected to be severe. To quantify adaptive differentiation in the arctic-alpine perennial herb Arabis alpina, we conducted reciprocal transplant experiments for 3 yr between Spanish and Scandinavian populations. At the sites of one Spanish and one Scandinavian population, we planted seedlings representing two Spanish and four Scandinavian populations, and recorded survival, flowering propensity and fecundity. The experiment was replicated in two subsequent years. The results demonstrate strong adaptive differentiation between A. alpina populations from the two regions. At the field site in Spain, survival and fruit production of Spanish populations were higher than those of Scandinavian populations, while the opposite was true at the site in Scandinavia, and these differences were consistent across years. By comparison, fitness varied little among populations from the same region. The results suggest that the magnitude and geographical scale of local adaptation need to be considered in predictions of the effects of global change on the dynamics of arctic and alpine plant populations.

On selection for flowering time plasticity in response to density

Different genotypes often exhibit opposite plastic responses in the timing of the onset of flowering with increasing plant density. In experimental studies, selection for accelerated flowering is generally found. By contrast, game theoretical studies predict that there should be selection for delayed flowering when competition increases. Combining different optimality criteria, the conditions under which accelerated or delayed flowering in response to density would be selected for are analysed with a logistic growth simulation model. To maximize seed production at the whole-stand level (simple optimization), selection should lead to accelerated flowering at high plant density, unless very short growing seasons select for similar onset of flowering at all densities. By contrast, selection of relative individual fitness will lead to delayed flowering when season length is long and/or growth rates are high. These different results give a potential explanation for the observed differences in direction of the plastic responses within and between species, including homeostasis, as a result of the effect of the variation in season length on the benefits of delayed flowering. This suggests that limited plasticity can evolve without the costs and limits that are currently thought to constrain the evolution of plasticity.

Variable response of three Trifolium repens ecotypes to soil flooding by seawater

BACKGROUND AND AIMS

Despite concerns about the impact of rising sea levels and storm surge events on coastal ecosystems, there is remarkably little information on the response of terrestrial coastal plant species to seawater inundation. The aim of this study was to elucidate responses of a glycophyte (white clover, Trifolium repens) to short-duration soil flooding by seawater and recovery following leaching of salts.

METHODS

Using plants cultivated from parent ecotypes collected from a natural soil salinity gradient, the impact of short-duration seawater soil flooding (8 or 24 h) on short-term changes in leaf salt ion and organic solute concentrations was examined, together with longer term impacts on plant growth (stolon elongation) and flowering.

KEY RESULTS

There was substantial Cl(-) and Na(+) accumulation in leaves, especially for plants subjected to 24 h soil flooding with seawater, but no consistent variation linked to parent plant provenance. Proline and sucrose concentrations also increased in plants following seawater flooding of the soil. Plant growth and flowering were reduced by longer soil immersion times (seawater flooding followed by drainage and freshwater inputs), but plants originating from more saline soil responded less negatively than those from lower salinity soil.

CONCLUSIONS

The accumulation of proline and sucrose indicates a potential for solute accumulation as a response to the osmotic imbalance caused by salt ions, while variation in growth and flowering responses between ecotypes points to a natural adaptive capacity for tolerance of short-duration seawater soil flooding in T. repens. Consequently, it is suggested that selection for tolerant ecotypes is possible should the predicted increase in frequency of storm surge flooding events occur.

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57-kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4-kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.

Investigating incipient speciation in Arabidopsis lyrata from patterns of transmission ratio distortion

Our understanding of the development of intrinsic reproductive isolation is still largely based on theoretical models and thorough empirical studies on a small number of species. Theory suggests that reproductive isolation develops through accumulation of epistatic genic incompatibilities, also known as Bateson-Dobzhansky-Muller (BDM) incompatibilities. We can detect these from marker transmission ratio distortion (TRD) in hybrid progenies of crosses between species or populations, where TRD is expected to result from selection against heterospecific allele combinations in hybrids. TRD may also manifest itself because of intragenomic conflicts or competition between gametes or zygotes. We studied early stage speciation in Arabidopsis lyrata by investigating patterns of TRD across the genome in F2 progenies of three reciprocal crosses between four natural populations. We found that the degree of TRD increases with genetic distance between crossed populations, but also that reciprocal progenies may differ substantially in their degree of TRD. Chromosomes AL6 and especially AL1 appear to be involved in many single- and two-locus distortions, but the location and source of TRD vary between crosses and between reciprocal progenies. We also found that the majority of single- and two-locus TRD appears to have a gametic, as opposed to zygotic, origin. Thus, while theory on BDM incompatibilities is typically illustrated with derived nuclear alleles proving incompatible in hybrid zygotes, our results suggest a prominent role for distortions emerging before zygote formation.