Exogenous selection shapes germination behaviour and seedling traits of populations at different altitudes in a Senecio hybrid zone

Background and Aims

The Senecio hybrid zone on Mt Etna, Sicily, is characterized by steep altitudinal clines in quantitative traits and genetic variation. Such clines are thought to be maintained by a combination of ‘endogenous’ selection arising from genetic incompatibilities and environment-dependent ‘exogenous’ selection leading to local adaptation. Here, the hypothesis was tested that local adaptation to the altitudinal temperature gradient contributes to maintaining divergence between the parental species, S. chrysanthemifolius and S. aethnensis.

Methods

Intra- and inter-population crosses were performed between five populations from across the hybrid zone and the germination and early seedling growth of the progeny were assessed.

Key Results

Seedlings from higher-altitude populations germinated better under low temperatures (9–13 °C) than those from lower altitude populations. Seedlings from higher-altitude populations had lower survival rates under warm conditions (25/15 °C) than those from lower altitude populations, but also attained greater biomass. There was no altitudinal variation in growth or survival under cold conditions (15/5 °C). Population-level plasticity increased with altitude. Germination, growth and survival of natural hybrids and experimentally generated F₁s generally exceeded the worse-performing parent.

Conclusions

Limited evidence was found for endogenous selection against hybrids but relatively clear evidence was found for divergence in seed and seedling traits, which is probably adaptive. The combination of low-temperature germination and faster growth in warm conditions might enable high-altitude S. aethnensis to maximize its growth during a shorter growing season, while the slower growth of S. chrysanthemifolius may be an adaptation to drought stress at low altitudes. This study indicates that temperature gradients are likely to be an important environmental factor generating and maintaining adaptive divergence across the Senecio hybrid zone on Mt Etna.     

RECENT ECOLOGICAL SELECTION ON REGULATORY DIVERGENCE IS SHAPING CLINAL VARIATION IN SENECIO ON MOUNT ETNA

The hybrid zone on Mount Etna (Sicily) between Senecio aethnensis and Senecio chrysanthemifolius (two morphologically and physiologically distinct species) is a classic example of an altitudinal cline. Hybridization at intermediate altitudes and gradients in phenotypic and life‐history traits occur along altitudinal transects of the volcano. The cline is considered to be a good example of ecological selection with species differences arising by divergent selection opposing gene flow. However, the possibility that the cline formed from recent secondary contact following an allopatric phase is difficult to exclude. We demonstrate a recent split between S. aethnensis and S. chrysanthemifolius (as recent as ～32,000 years ago) and sufficient gene flow (2Nm > 1) to have prevented divergence (implicating a role for diversifying selection in the maintenance of the cline). Differentially expressed genes between S. aethnensis and S. chrysanthemifolius exhibit significantly higher genetic divergence relative to “expression invariant” controls, suggesting that species differences may in part be mediated by divergent selection on differentially expressed genes involved with altitude‐related adaptation. The recent split time and the absence of fixed differences between these two ecologically distinct species suggest the rapid evolution to an altitudinal cline involving selection on both sequence and expression variation.     

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

THE POPULATION GENETICS OF SPOROPHYTIC SELF‐INCOMPATIBILITY IN THREE HYBRIDIZING SENECIO (ASTERACEAE) SPECIES WITH CONTRASTING POPULATION HISTORIES

Hybridization generates evolutionary novelty and spreads adaptive variation. By promoting outcrossing, plant self‐incompatibility (SI) systems also favor interspecific hybridization because the S locus is under strong negative frequency‐dependent balancing selection. This study investigates the SI mating systems of three hybridizing Senecio species with contrasting population histories. Senecio aethnensis and S. chrysanthemifolius native to Sicily, form a hybrid zone at intermediate altitudes on Mount Etna, and their neo‐homoploid hybrid species, S. squalidus, has colonized disturbed urban habitats in the UK during the last 150 years. We show that all three species express sporophytic SI (SSI), where pollen incompatibility is controlled by the diploid parental genome, and that SSI is inherited and functions normally in hybrids. Large‐scale crossing studies of wild sampled populations allowed direct comparison of SSI between species and found that the main impacts of colonization in S. squalidus compared to Sicilian Senecio was a reduced number of S alleles, increased S allele frequencies, and increased interpopulation S allele sharing. In general, many S alleles were shared between species and the S locus showed reduced intra‐ and interspecific population genetic structure compared to molecular genetic markers, indicative of enhanced effective gene flow due to balancing selection.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Genome-wide associations between hybrid sterility QTL and marker transmission ratio distortion

Marker transmission ratio distortion (TRD) in genetic mapping populations is frequently ascribed to selection against allelic combinations that cause hybrid incompatibility. Accordingly, genomic regions of TRD should be nonrandomly associated (colocated) with loci that underlie hybrid incompatibility. To directly test this hypothesis, we evaluated the genome-wide qualitative and quantitative agreement between chromosomal regions exhibiting marker TRD and those known to contain hybrid incompatibility quantitative trait locus (QTL). Incompatibility data came from a near-isogenic line (NIL) analysis of pollen and seed sterility in a cross between two Solanum (formerly Lycopersicon) species. We assessed (1) whether these incompatibility loci are colocated with markers that show significant TRD in two earlier generations preceding these introgression lines and (2) whether the magnitude of marker distortion quantitatively matches the estimated strength of selection against each incompatibility locus. We found evidence that TRD regions are chromosomally colocated with hybrid incompatibility loci more frequently than is expected by chance: pollen sterility QTLs were most closely associated with distorted heterozygote frequencies in later-generation backcrosses. Nonetheless, there was no evidence for an association between TRD and seed sterility and little evidence of a quantitative association between the magnitude of marker TRD and the fitness effects of heterospecific alleles at each chromosomal location. We propose and test a model (the "dance partner" model) to explain several cases where regions of TRD are not associated with hybrid incompatibility loci. Under this model, some NILs containing greater than one heterospecific introgression may not express hybrid incompatibility phenotypes because they carry both appropriate genetic dance partners required for a fully functional interaction. Accordingly, negative interactions expressed in earlier backcross generations are masked in these double-introgression NILs. Based on this model, we identify the location of several new putative pairwise interactors underlying hybrid incompatibility in this species cross.     

The genomics of incompatibility factors and sex determination in hybridizing species of Cottus (Pisces)

Cottus rhenanus and Cottus perifretum have formed hybrid lineages and narrow hybrid zones that can be best explained through the action of natural selection. However, the underlying selective forces as well as their genomic targets are not well understood. This study identifies genomic regions in the parental species that cause hybrid incompatibilities and tests whether these manifest in a sex-specific manner to learn about processes that affect natural hybridization in Cottus. Interspecific F₂ crosses were analyzed for 255 markers for genetic mapping and to detect transmission distortion as a sign for genetic incompatibilities. The Cottus map consists of 24 linkage groups with a total length of 1575.4 cM. A male heterogametic (XY) sex determination region was found on different linkage groups in the two parental species. Genetic incompatibilities were incomplete, varied among individuals and populations and were not associated with the heterogametic sex. The variance between populations and individuals makes it unlikely that there are species-specific incompatibility loci that could affect the gene pool of natural hybrids in a simple and predictable way. Conserved synteny with sequenced fish genomes permits to genetically study the Cottus genome through the transfer of genomic information from the model fish species. Homology relationships of candidate genomic regions in Cottus indicate that sex determination is not based on the same genomic regions found in other fish species. This suggests a fast evolutionary turnover of the genetic basis of sex determination that, together with the small size of the heterogametic regions, may contribute to the absence of fitness effects related to the Haldane''s rule.     

Genomic collinearity and the genetic architecture of floral differences between the homoploid hybrid species Iris nelsonii and one of its progenitors,Iris hexagona

Hybrid speciation represents a relatively rapid form of diversification. Early models of homoploid hybrid speciation suggested that reproductive isolation between the hybrid species and progenitors primarily resulted from karyotypic differences between the species. However, genic incompatibilities and ecological divergence may also be responsible for isolation. Iris nelsonii is an example of a homoploid hybrid species that is likely isolated from its progenitors primarily by strong prezygotic isolation, including habitat divergence, floral isolation and post-pollination prezygotic barriers. Here, we used linkage mapping and quantitative trait locus (QTL) mapping approaches to investigate genomic collinearity and the genetic architecture of floral differences between I. nelsonii and one of its progenitor species I. hexagona. The linkage map produced from this cross is highly collinear with another linkage map produced between I. fulva and I. brevicaulis (the two other species shown to have contributed to the genomic makeup of I. nelsonii), suggesting that karyotypic differences do not contribute substantially to isolation in this homoploid hybrid species. Similar to other studies of the genetic architecture of floral characteristics, at least one QTL was found that explained >20% variance in each color trait, while minor QTLs were detected for each morphological trait. These QTLs will serve as hypotheses for regions under selection by pollinators.     

Isolation and characterization of microsatellite loci in Senecio

Interspecific hybridization has resulted in the recent origin of several hybrid Senecio taxa at diploid, tetraploid and hexaploid levels. As part of research aimed at constructing and comparing genomic maps of each of these taxa and their parents, we have isolated microsatellite loci from genomic DNA libraries of S. vulgaris and S. squalidus. Primers of 35 loci amplified microsatellites resolved in agarose gels from one or more of S. vulgaris, S. squalidus, S. aethnensis and S. chrysanthemifolius. Approximately 71% of primers amplified a product in all four species. A survey of microsatellite variation in S. chrysanthemifolius over a subset of 14 loci resolved 2–11 alleles per locus in polyacrylamide gels with expected heterozygosity (H_E) ranging from 0.26 to 0.87.     

Genetic differentiation and species cohesion in two widespread Central American Begonia species

Begonia is one of the ten largest plant genera, with over 1500 species. This high species richness may in part be explained by weak species cohesion, which has allowed speciation by divergence in allopatry. In this study, we investigate species cohesion in the widespread Central American Begonia heracleifolia and Begonia nelumbiifolia, by genotyping populations at microsatellite loci. We then test for post-zygotic reproductive barriers using experimental crosses, and assess whether sterility barriers are related to intraspecific changes in genome size, indicating major genome restructuring between isolated populations. Strong population substructure was found for B. heracleifolia (F_ST=0.364, F′_ST=0.506) and B. nelumbiifolia (F_ST=0.277, F′_ST=0.439), and Bayesian admixture analysis supports the division of most populations into discrete genetic clusters. Moderate levels of inferred selfing (B. heracleifolia s=0.40, B. nelumbiifolia s=0.62) and dispersal limitation are likely to have contributed to significant genetic differentiation (B. heracleifolia Jost''s D=0.274; B. nelumbiifolia D=0.294). Interpopulation crosses involving a divergent B. heracleifolia population with a genome size ∼10% larger than the species mean had a ∼20% reduction in pollen viability compared with other outcrosses, supporting reproductive isolation being polymorphic within the species. The population genetic data suggest that Begonia populations are only weakly connected by gene flow, allowing reproductive barriers to accumulate between the most isolated populations. This supports allopatric divergence in situ being the precursor of speciation in Begonia, and may also be a common speciation mechanism in other tropical herbaceous plant groups.     

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 F₂ 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.     

Spatial and temporal genetic structure in a hybrid cordgrass invasion

Invasive hybrids and their spread dynamics pose unique opportunities to study evolutionary processes. Invasive hybrids of native Spartina foliosa and introduced S. alterniflora have expanded throughout San Francisco Bay intertidal habitats within the past 35 years by deliberate plantation and seeds floating on the tide. Our goals were to assess spatial and temporal scales of genetic structure in Spartina hybrid populations within the context of colonization history. We genotyped adult and seedling Spartina using 17 microsatellite loci and mapped their locations in three populations. All sampled seedlings were hybrids. Bayesian ordination analysis distinguished hybrid populations from parent species, clearly separated the population that originated by plantation from populations that originated naturally by seed and aligned most seedlings within each population. Population genetic structure estimated by analysis of molecular variance was substantial (F_ST=0.21). Temporal genetic structure among age classes varied highly between populations. At one population, the divergence between adults and 2004 seedlings was low (F_ST=0.02) whereas at another population this divergence was high (F_ST=0.26). This latter result was consistent with local recruitment of self-fertilized seed produced by only a few parental plants. We found fine-scale spatial genetic structure at distances less than ∼200 m, further supporting local seed and/or pollen dispersal. We posit a few self-fertile plants dominating local recruitment created substantial spatial genetic structure despite initial long-distance, human dispersal of hybrid Spartina through San Francisco Bay. Fine-scale genetic structure may more strongly develop when local recruits are dominated by the offspring of a few self-fertile plants.     

Functional divergence caused by ancient positive selection of a Drosophila hybrid incompatibility locus

Interspecific hybrid lethality and sterility are a consequence of divergent evolution between species and serve to maintain the discrete identities of species. The evolution of hybrid incompatibilities has been described in widely accepted models by Dobzhansky and Muller where lineage-specific functional divergence is the essential characteristic of hybrid incompatibility genes. Experimentally tractable models are required to identify and test candidate hybrid incompatibility genes. Several Drosophila melanogaster genes involved in hybrid incompatibility have been identified but none has yet been shown to have functionally diverged in accordance with the Dobzhansky-Muller model. By introducing transgenic copies of the X-linked Hybrid male rescue (Hmr) gene into D. melanogaster from its sibling species D. simulans and D. mauritiana, we demonstrate that Hmr has functionally diverged to cause F1 hybrid incompatibility between these species. Consistent with the Dobzhansky-Muller model, we find that Hmr has diverged extensively in the D. melanogaster lineage, but we also find extensive divergence in the sibling-species lineage. Together, these findings implicate over 13% of the amino acids encoded by Hmr as candidates for causing hybrid incompatibility. The exceptional level of divergence at Hmr cannot be explained by neutral processes because we use phylogenetic methods and population genetic analyses to show that the elevated amino-acid divergence in both lineages is due to positive selection in the distant past—at least one million generations ago. Our findings suggest that multiple substitutions driven by natural selection may be a general phenomenon required to generate hybrid incompatibility alleles.     

GENETIC ARCHITECTURE AND POSTZYGOTIC REPRODUCTIVE ISOLATION: EVOLUTION OF BATESON–DOBZHANSKY–MULLER INCOMPATIBILITIES IN A POLYGENIC MODEL

The Bateson–Dobzhansky–Muller model predicts that postzygotic isolation evolves due to the accumulation of incompatible epistatic interactions, but few studies have quantified the relationship between genetic architecture and patterns of reproductive divergence. We examined how the direction and magnitude of epistatic interactions in a polygenic trait under stabilizing selection influenced the evolution of hybrid incompatibilities. We found that populations evolving independently under stabilizing selection experienced suites of compensatory allelic changes that resulted in genetic divergence between populations despite the maintenance of a stable, high‐fitness phenotype. A small number of loci were then incompatible with multiple alleles in the genetic background of the hybrid and the identity of these incompatibility loci changed over the evolution of the populations. For F₁ hybrids, reduced fitness evolved in a window of intermediate strengths of epistatic interactions, but F₂ and backcross hybrids evolved reduced fitness across weak and moderate strengths of epistasis due to segregation variance. Strong epistatic interactions constrained the allelic divergence of parental populations and prevented the development of reproductive isolation. Because many traits with varying genetic architectures must be under stabilizing selection, our results indicate that polygenetic drift is a plausible hypothesis for the evolution of postzygotic reproductive isolation.     

Reproductive and genetic roles of the maternal progenitor in the origin of common wheat (Triticum aestivum L.)

Common wheat (Triticum aestivum L., AABBDD genome) is thought to have emerged through natural hybridization between Triticum turgidum L. (AABB genome) and Aegilops tauschii Coss. (DD genome). Hybridization barriers and doubling of the trihaploid F₁ hybrids’ genome (ABD) via unreduced gamete fusion had key roles in the process. However, how T. turgidum, the maternal progenitor, was involved in these mechanisms remains unknown. An artificial cross‐experiment using 46 cultivated and 31 wild T. turgidum accessions and a single Ae. tauschii tester with a very short genetic distance to the common wheat D genome was conducted. Cytological and quantitative trait locus analyses of F₁ hybrid genome doubling were performed. The crossability and ability to cause hybrid inviability did not greatly differ between the cultivars and wild accessions. The ability to cause hybrid genome doubling was higher in the cultivars. Three novel T. turgidum loci for hybrid genome doubling, which influenced unreduced gamete production in F₁ hybrids, were identified. Cultivated T. turgidum might have increased the probability of the emergence of common wheat through its enhanced ability to cause genome doubling in F₁ hybrids with Ae. tauschii. The ability enhancement might have involved alterations at a relatively small number of loci.     

Altitudinal gradients,plant hybrid zones and evolutionary novelty

Altitudinal gradients are characterized by steep changes of the physical and biotic environment that present challenges to plant adaptation throughout large parts of the world. Hybrid zones may form where related species inhabit different neighbouring altitudes and can facilitate interspecific gene flow and potentially the breakdown of species barriers. Studies of such hybrid zones can reveal much about the genetic basis of adaptation to environmental differences stemming from changes in altitude and the maintenance of species divergence in the face of gene flow. Furthermore, owing to recombination and transgressive effects, such hybrid zones can be sources of evolutionary novelty. We document plant hybrid zones associated with altitudinal gradients and emphasize similarities and differences in their structure. We then focus on recent studies of a hybrid zone between two Senecio species that occur at high and low altitude on Mount Etna, Sicily, showing how adaptation to local environments and intrinsic selection against hybrids act to maintain it. Finally, we consider the potential of altitudinal hybrid zones for generating evolutionary novelty through adaptive introgression and hybrid speciation. Examples of homoploid hybrid species of Senecio and Pinus that originated from altitudinal hybrid zones are discussed.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.     

Transmission ratio distortion in an interspecific cross between Fusarium circinatum and Fusarium subglutinans

Previously, an interspecific cross between Fusarium circinatum and Fusarium subglutinans was used to generate a genetic linkage map. A ca. 55 % of genotyped markers displayed transmission ratio distortion (TRD) that demonstrated a genome-wide distribution. The working hypothesis for this study was that TRD would be non-randomly distributed throughout the genetic linkage map. This would indicate the presence of distorting loci. Using a genome-wide threshold of α = 0.01, 79 markers displaying TRD were distributed on all 12 linkage groups (LGs). Eleven putative transmission ratio distortion loci (TRDLs), spanning eight LGs, were identified in regions containing three or more adjacent markers displaying distortion. No epistatic interactions were observed between these TRDLs. Thus, it is uncertain whether the genome-wide TRD was due to linkage between markers and genomic regions causing distortion. The parental origins of markers followed a non-random distribution throughout the linkage map, with LGs containing stretches of markers originating from only one parent. Thus, due to the nature of the interspecific cross, the current hypothesis to explain these observations is that the observed genome-wide segregation was caused by the high level of genomic divergence between the parental isolates. Therefore, homologous chromosomes do not align properly during meiosis, resulting in aberrant transmission of markers. This also explains previous observations of the preferential transmission of F. subglutinans alleles to the F₁ progeny.     

Rapid genomic DNA changes in allotetraploid fish hybrids

Rapid genomic change has been demonstrated in several allopolyploid plant systems; however, few studies focused on animals. We addressed this issue using an allotetraploid lineage (4nAT) of freshwater fish originally derived from the interspecific hybridization of red crucian carp (Carassius auratus red var., ♀, 2n=100) × common carp (Cyprinus carpio L., ♂, 2n=100). We constructed a bacterial artificial chromosome (BAC) library from allotetraploid hybrids in the 20th generation (F₂₀) and sequenced 14 BAC clones representing a total of 592.126 kb, identified 11 functional genes and estimated the guanine–cytosine content (37.10%) and the proportion of repetitive elements (17.46%). The analysis of intron evolution using nine orthologous genes across a number of selected fish species detected a gain of 39 introns and a loss of 30 introns in the 4nAT lineage. A comparative study based on seven functional genes among 4nAT, diploid F₁ hybrids (2nF₁) (first generation of hybrids) and their original parents revealed that both hybrid types (2nF₁ and 4nAT) not only inherited genomic DNA from their parents, but also demonstrated rapid genomic DNA changes (homoeologous recombination, parental DNA fragments loss and formation of novel genes). However, 4nAT presented more genomic variations compared with their parents than 2nF₁. Interestingly, novel gene fragments were found for the iqca1 gene in both hybrid types. This study provided a preliminary genomic characterization of allotetraploid F₂₀ hybrids and revealed evolutionary and functional genomic significance of allopolyploid animals.     

Airborne‐pollen pool and mating pattern in a hybrid zone between Pinus pumila and P. parviflora var. pentaphylla

The reproductive isolation barriers and the mating patterns among Pinus pumila, P. parviflora var. pentaphylla and their hybrids were examined by flowering phenology and genetic assays of three life stages: airborne‐pollen grains, adults and seeds, in a hybrid zone on Mount Apoi, Hokkaido, Japan. Chloroplast DNA composition of the airborne‐pollen was determined by single‐pollen polymerase chain reaction. Mating patterns were analysed by estimating the molecular hybrid index of the seed parent, their seed embryos and pollen parents. The observation of flowering phenology showed that the flowering of P. pumila precedes that of P. parviflora var. pentaphylla by about 6 to 10 days within the same altitudinal ranges. Although this prezygotic isolation barrier is effective, the genetic assay of airborne‐pollen showed that the two pine species, particularly P. pumila, still have chances to form F₁ hybrid seeds. Both parental species showed a strong assortative mating pattern; F₁ seeds were found in only 1.4% of seeds from P. pumila mother trees and not at all in P. parviflora var. pentaphylla. The assortative mating was concluded as the combined result of flowering time differentiation and cross‐incompatibility. In contrast to the parental species, hybrids were fertilized evenly by the two parental species and themselves. The breakdown of prezygotic barriers (intermediate flowering phenology) and cross‐incompatibility may account for the unselective mating. It is suggested that introgression is ongoing on Mount Apoi through backcrossing between hybrids and parental species, despite strong isolation barriers between the parental species.