Biosystematics and evolutionary relationships of perennial Triticeae species revealed by genomic analyses

Understanding the classification and biosystematics of species in Triticeae Dumort., an economically important tribe in the grass family (Poaceae), is not an easy task, particularly for some perennial species. Does genomic analysis facilitate the understanding of evolutionary relationships of these Triticeae species? We reviewed literature published after 1984 to address questions concerning: (1) genome relationships among the monogenomic diploid species; (2) progenitors of the unknown Y genome in Elymus polyploids, X genome in Thinopyrum intermedium, and Xm genome in Leymus; and (3) genome constitutions of some perennial Triticeae species that were unknown or misidentified. A majority of publications have substantiated the close affinity of the E^b and E^e genomes in Th. bessarabicumand Th. elongatum, supporting the use of a common basic genome symbol. The E genome is close to the St genome of Pseudoroegneria and ABD genomes ofTriticum/Aegilops complex, providing an explanation for transferring genes from the E to ABD genomes with relative ease. Although the solid proof is still lacking, theW, P, and especially Xp genomes are possible origins for the Y genome of polyploid Elymus. The absence of the E genome and the allopolyploidy nature of tetraploidLeymus species have been unequivocally confirmed by both cytogenetic and molecular studies. However, the donor of the Xm genomes of Leymus was only speculated to be related to the P genome of Agropyron and F genome of Eremopyrum. Intermediate wheatgrass (Th. intermedium) has been extensively studied. The presence of the St (as the previously designated X) genome in Th. intermedium is now unequivocal. Its two more closely related E₁ and E₂ genomes are shown to be older versions of the E genome rather than the current E^b and E^e genomes. Speciation of Th. intermedium was similar to that of Triticum aestivum, in which the J^s/E^s(like B) genomes had the greatest differentiation from the current J (E^b) genome owning to repetitive sequences of the V genome, whereas its St (like D) had the least differentiation from the current St genome. Species with unknown or misidentified genomes have been correctly designated, including those with the ESt, StP, StPY,StWY, EStP, HW, StYHW, and NsXm genomes. Some of those species have been transferred to and renamed in appropriate genera.     

Phylogenetic analysis of North American Elymus and the monogenomic Triticeae (Poaceae) using three chloroplast DNA data sets.

Although the monogenomic genera of the Triticeae have been analyzed in numerous biosystematic studies, the allopolyploid genera have not been as extensively studied within a phylogenetic framework. We focus on North American species of Elymus, which, under the current genomic system of classification, are almost all allotetraploid, combining the St genome of Pseudoroegneria with the H genome of Hordeum. We analyze new and previously published chloroplast DNA data from Elymus and from most of the monogenomic genera of the Triticeae in an attempt to identify the maternal genome donor of Elymus. We also present a cpDNA phylogeny for the monogenomic genera that includes more data than, and thus builds on, those previously published. The chloroplast DNA data indicate that Pseudoroegneria is the maternal genome donor to all but one of the Elymus individuals. There is little divergence among the Elymus and Pseudoroegneria chloroplast genomes, and as a group, they show little divergence from the rest of the Triticeae. Within the monogenomic Triticeae, the problematic group Thinopyrum is resolved as monophyletic on the chloroplast DNA tree. At the intergeneric level, the data reveal several deeper-level relationships that were not resolved by previous cpDNA trees.     

Phylogenetic relationships in Elymus (Poaceae: Triticeae) based on the nuclear ribosomal internal transcribed spacer and chloroplast trnL-F sequences

To estimate the phylogenetic relationship of polyploid Elymus in Triticeae, nuclear ribosomal internal transcribed spacer (ITS) and chloroplast trnL-F sequences of 45 Elymus accessions containing various genomes were analysed with those of five Pseudoroegneria (St), two Hordeum (H), three Agropyron (P) and two Australopyrum (W) accessions. The ITS sequences revealed a close phylogenetic relationship between the polyploid Elymus and species from the other genera. The ITS and trnL-F trees indicated considerable differentiation of the StY genome species. The trnL-F sequences revealed an especially close relationship of Pseudoroegneria to all Elymus species included. Both the ITS and trnL-F trees suggested multiple origins and recurrent hybridization of Elymus species. The results suggested that: the St, H, P, and W genomes in polyploid Elymus were donated by Pseudoroegneria, Hordeum, Agropyron and Australopyrum, respectively, and the St and Y genomes may have originated from the same ancestor; Pseudoroegneria was the maternal donor of the polyploid Elymus; and some Elymus species showed multiple origin and experienced recurrent hybridization.     

Molecular diversity of the 5S nuclear ribosomal DNA in Campeiostachys with StHY haplome constitution

To detect the genomic constitutions and investigate the evolutionary relationships between Campeiostachys Drobov and Elymus L. species, we have cloned and analyzed 271 5S nuclear ribosomal DNA sequences from 27 accessions of these species, mostly of Chinese origin. We identified Long H1, Short S1, and Long Y1 unit classes in nine Campeiostachys or Elymus species. The identification of the three orthologous unit classes was confirmed by the neighbor‐joining tree of each unit class from PAUP and the phylogeny tree of three unit classes from MrBayes. The results suggested that these Elymus species comprise StYH haplomes and should be included in Campeiostachys. The phylogeny tree showed a clear separation between the S1 unit class and Y1 unit class. However, Y1 unit class sequences formed a sister clade to the S1 unit class, implying that although the St and Y haplomes might have some affinity, they are distinct from one another. The phylogeny tree also indicated that the five species in sect. Turczaninovia (C. dahurica var. cylindrica, C. dahurica var. dahurica, C. dahurica var. tangutorum, E. purpuraristatus, and E. dahuricus Turcz. ex Griseb. var. violeus C. P. Wang & H. L. Yang) might share a more recent common ancestor, whereas the four species in sect. Elymus (C. nutans, E. breviaristatus (Keng) Keng ex Keng f., E. sinosubmuticus (Keng) Keng f., and E. atratus (Nevski) Hand.‐Mazz.) share a close relationship. By identifying only one type of unit class for each haplome, we propose that the 5S nuclear ribosomal DNA sequences of species within Campeiostachys might have undergone haplome‐specific concerted evolution.     

Reticulate evolution, introgression, and intertribal gene capture in an allohexaploid grass

Recent molecular phylogenetic studies of polyploid plants have successfully clarified complex patterns of reticulate evolution. In this study of Elymus repens, an allohexaploid member of the wheat tribe Triticeae, chloroplast and nuclear DNA data reveal an extreme reticulate pattern, revealing at least five distinct gene lineages coexisting within the species, acquired through a possible combination of allohexaploidy and introgression from both within and beyond the Triticeae. Earlier cytogenetic studies of E. repens suggested that Hordeum (genome H) and Pseudoroegneria (St) were genome donors to E. repens. Chloroplast DNA data presented here (from the rpoA gene and from the region between trnT and trnF) identify three potential maternal genome donors (Pseudoroegneria, Thinopyrum, and Dasypyrum), and information from previous molecular work suggests that, of these, Pseudoroegneria is the most likely maternal donor. Nuclear starch synthase gene data indicate that both Hordeum and Pseudoroegneria have contributed to the nuclear genome of E. repens, in agreement with cytogenetic data. However, these data also show unexpected contributions from Taeniatherum, and from two additional donors of unknown identity. One of the sequences of unknown origin falls within the Triticeae, but is not closely associated with any of the sampled diploid genera. The second falls outside of the clade containing Triticeae and its outgroup Bromus, suggesting the acquisition of genetic material from a surprisingly divergent source. Bias toward the amplification of certain starch synthase variants has complicated attempts to thoroughly sample from within individuals, but the data clearly indicate a complex pattern of reticulate evolution, consistent not only with allohexaploidy, but also with introgression from unexpectedly divergent sources.     

Multiple origins of allopolyploid wheatgrass Elymus caninus revealed by RPB2, PepC and TrnD/T genes

We examined evolutionary mechanisms in the tetraploid Elymus caninus by comparing the phylogenetic relationships of 21 accessions suggested by sequence data from two single copy nuclear genes, the largest subunit of RNA polymerase II (RPB2) and phosphoenolpyruvate carboxylase (pepC), and one non-coding chloroplast region, TrnD/T. Elymus caninus is known combining two different genomes, an St genome and an H genome. Data from two single copy nuclear genes showed that there are two versions of the St genome in the species, St1 and St2. Most accessions combined one of these versions with an H genome version but two accessions had both versions of the St sequence for RPB2. This suggests that the RPB2gene may have been duplicated without chromosome doubling, possibly induced by transposable element. Our data also indicate that the H genome sequences in E. caninus have multiple origins, and a close phylogenetic relationship between Hordeum bogdanii and H sequences in some accessions of E. caninus. Thus, it is more likely that H. bogdanii is one of the major donors of the H copy in E. caninus. The maternal origin of E. caninus is the St genome species. There was no correlation between the geographic origin of the accessions and their sequence divergence.     

Phylogenetic relationships and the maternal donor of Roegneria (Triticeae: Poaceae) based on three nuclear DNA sequences (ITS,Acc1, and Pgk1) and one chloroplast region (trnL-F)

Roegneria is a polyploid perennial genus in the tribe Triticeae. Some species of Roegneria are morphologically similar to genus Elymus and have been classified in Elymus. To investigate the delimitation and phylogenetic relationships of Roegneria, nuclear (ITS, Acc1, and Pgk1) and chloroplast (trnL–trnF) DNA regions were sequenced for 38 allopolyploid species and 32 diploid species of Triticeae. Phylogenetic analyses of nuclear DNA revealed that all Roegneria species were included in the St and Y genome clades, and that the Y genome was closely related to the V and Xp genomes. The chloroplast DNA dataset showed that Roegneria species were grouped with Pseudoroegneria species. The Pseudoroegneria species from the Middle East (P. libanotica and P. tauri) and Central Asia (P. strigosa) were more closely related to Roegneria species. The results suggested that: (i) the species containing the St and Y genomes should be segregated from Elymus and treated as a distinct genus, Roegneria, based on the genomic constitution; (ii) P. libanotica, P. tauri, and/or P. strigosa potentially served as the maternal donor of the St genome in Roegneria; (iii) The Y genome of Roegneria originated from a diploid Y genome species, and the V and Xp genomes may have contributed to Y genome formation; (iv) among Roegneria species of previously uncertain genomic constitution, R. seriotina was tetraploid and possessed the StY genomes, E. calcicolus was hexaploid with the StYH genomic constitution and should be classified in Campeiostachys, R. glaucifolia possessed the StStY genomes, and R. tschimganica had the genomic constitution St₁St₂Y.     

Molecular phylogeny of RPB2 gene reveals multiple origin, geographic differentiation of H genome, and the relationship of the Y genome to other genomes in Elymus species

It has been hypothesized from isozymic and cytological studies of Elymus species that the Old and New World taxa may be of separate origin of the H genome in the StH genome species. To test this hypothesis, and estimate the phylogenetic relationships of polyploid Elymus species within the Triticeae, the second largest subunit of RNA polymerase II (RPB2) sequence of 36 Elymus accessions containing StH or StY genomes was analyzed with those of Pseudoroegneria (St), Hordeum (H), Agropyron (P), Australopyrum (W), Lophopyrum(Ee), Thinopyrum(Eb) and Dasypyrum (V). Our data indicated that the H genome in Elymus species is differentiated in accordance with geographical origin, and that the Eurasian and American StH genome species have independent alloploid origins with different H-genome donors. Phylogenetic analysis of Y genome sequences with other genome donors (St, H, P, W) of Elymus revealed that W and P genomes are sister to Y genome with a 87% bootstrap support, and that StY and StH species group might have acquired their RPB2 St sequences from distinct Pseudoroegneria gene pools. Our data did not support the suggestion that the St and Y genomes have the same origin as put forward in a previous study using ITS data. Our result provides some insight on the origin of Y genome and its relationship to other genomes in Elymus.     

Genome constitution of Elymus tangutorum (Poaceae: Triticeae) inferred from meiotic pairing behavior and genomic in situ hybridization

Elymus tangutorum (Nevski) Hand.-Mazz (2n = 6x = 42) is a perennial species in the tribe Triticeae, which distributes in Nepal and north and northwest China. However, the genome constitution of E. tangutorum is controversial and its taxonomic status is not clear. Hybridizations of E. tangutorum were carried out with E. wawawaiensis J. R. Carlson & Barkworth (StH), Roegneria grandis Keng (StY), and E. dahuricus Turcz. ex Griseb. (StYH). Meiotic pairing in the hybrids E. tangutorum × E. wawawaiensis (StH), E. tangutorum × R. grandis (StY), and E. tangutorum × E. dahuricus (StYH) averaged 10.48, 11.12, and 20.92 bivalents per cell, respectively. The results suggested that E. tangutorum is an allohexaploid and contains the StYH genomes. Results of genomic in situ hybridization analysis strongly supported the chromosome pairing data. Therefore, E. tangutorum should be treated as Campeiostachys dahurica var. tangutorum (Nevski) B. R. Baum, J. L. Yang & C. Yen. Intergenomic rearrangements of E. tangutorum may be affected by environmental factors.     

Origin of the H genome in StH-genomic Elymus species based on the single-copy nuclear gene DMC1

Previous studies have suggested that the H haplome in Elymus could originate from different diploid Hordeum species, however, which diploid species best represent the parental species remains unanswered. The focus of this study seeks to pinpoint the origin of the H genome in Elymus. Allopolyploid Elymus species that contain the StH genome were analyzed together with diploid Hordeum species and a broad sample of diploid genera in the tribe Triticeae using DMC1 sequences. Both parsimony and maximum likelihood analyses well separated the American Hordeum species, except Hordeum brachyantherum subsp. californicum, from the H genome of polyploid Elymus species. The Elymus H-genomic sequences were formed into different groups. Our data suggested that the American Horedeum species, except H. brachyantherum subsp. californicum, are not the H-genomic donor to the Elymus species. Hordeum brevisubulatum subsp. violaceum was the progenitor species to Elymus virescens, Elymus confusus, Elymus lanceolatus, Elymus wawawaiensis, and Elymus caninus. Furthermore, North American H. brachyantherum subsp. californicum was a progenitor of the H genome to Elymus hystrix and Elymus cordilleranus. The H genomes in Elymus canadensis, Elymus sibiricus, and Elymus multisetus were highly differentiated from the H genome in Hordeum and other Elymus species. The H genome in both North American and Eurasian Elymus species was contributed by different Hordeum species.     

Allohexaploidy, introgression, and the complex phylogenetic history of Elymus repens (Poaceae)

The phylogenetic position of hexaploid Elymus repens within the tribe Triticeae (Poaceae) was examined using cloned sequences from the low-copy nuclear genes encoding phosphoenolpyruvate carboxylase (pepC) and beta-amylase. A previous analysis of E. repens using data from the nuclear granule-bound starch synthase I (GBSSI) gene had yielded five phylogenetically distinct gene copies, two more than expected from hexaploidy alone. The three gene trees share three distinct E. repens clades, suggesting that E. repens contains three phylogenetically divergent genomes, contributed by Hordeum, Pseudoroegneria, and an unknown donor. The two additional GBSSI sequences, including one that was apparently derived from outside of the tribe, appear to reflect past introgression of GBSSI sequences into the E. repens genome. On all three trees, the Hordeum-like E. repens sequences are polyphyletic within Hordeum, and the trees are in conflict with regard to the placement of these sequences within Hordeum, highlighting multiple contributions from Hordeum to E. repens.     

Molecular cytogenetics of Leymus: Mapping the Ns genome-specific repetitive sequences

The objective of this paper is to summarize the work in my group on FISH (fluorescent in situ hybridization) mapping of Ns-specific repetitive DNA sequences fromLeymus and discuss the results in the context of classification based on the genome system currently used among Triticeae researchers. The key question here is whether the genome composition of a tetraploid Leymus species should be NsXm or NsNs (Ns₁Ns₂). Different types of Leymus-specific dispersed retroelement-like repeats have been isolated and characterized. Because the sequences occur in significantly high copy number in Leymus, based on strong hybridization signal in Southern blots, they are considered essentially specific to Leymus. They are also abundant in Psathyrostachys, the progenitor of Ns genome in Leymus. These dispersed repeats are found to distribute over the whole of all Leymus chromosomes, without any differentiation between chromosomes that have been suggested to be of different genomic origins, meaning that all genomes in Leymus are the same. GISH (genomic in situ hybridization) experiments on Leymus chromosomes using Psathyrostachys genomic DNA as probes further support the NsNs (Ns₁Ns₂) genome constitution for Leymus. The Xm genome of an unknown origin might have been there in the beginning of the allopolyploidization process, but the Ns genome-specific elements must have spread predominantly and rapidly across genomes, thus homogenizing the nuclear genomes of Leymus. I present here for the first time evidence that Ns-specific dispersed repeats can spread in a very short time, from Leymus over to wheat in Triticum × Leymus hybrids growing in artificial conditions.     

A study of 28 Elymus species using repetitive DNA sequences.

Four repetitive DNA sequences cloned from the barley (Hordeum vulgare) genome and common for different Triticeae species were used for a molecular study of phylogenetic relationships among 28 Elymus species. Two wild Hordeum species (H genome), two Pseudoroegneria species (S genome), Agropyron cristatum (P genome), and Australopyrum velutinum (W genome) were included as genomic representatives for the genomes that supposedly were involved in the evolution of the genus Elymus. Our results are essentially congruent with the genomic classification system. This study demonstrates that Elymus is not a monophyletic genus. Based on an analysis of Southern blot hybridization we could discriminate between SY and SH species owing to the strong specific hybridization pattern of the H genome. Hexaploid SYH species gave a hybridization pattern similar to SH species for the same reason. The results support the genomic composition of Elymus batalinii as SYP and also indicated the presence of at least one H genome in Elymus enysii with a hitherto unknown genomic constitution. Elymus erianthus had a hybridization pattern distinctly different from all other species in the investigation. Key words : Elymus, RFLP, phylogeny, repetitive DNA.     

Phylogeny and molecular evolution of the rbcL gene of St genome in Elymus sensu lato (Poaceae: Triticeae)

Analysis of the patterns and levels of diversity in duplicate gene not only traces evolutionary history of polyploids, but also provides insight into how the evolutionary process differs between lineages and between homoeologous loci within lineages. Elymus sensu lato is a group of allopolyploid species, which share a common St genome and with the different combinations of H, Y, P, and W genomes. To estimate the evolutionary process of the rbcL gene in species of Elymus s. l. and its putative dioploid relatives, 74 sequences were obtained from 21 species of Elymus s. l. together with 24 diploid taxa representing 19 basic genomes in Triticeae. Phylogeny and sequence diversity pattern analysis suggested that (1) species of Pseudoroegneria (Nevski) Á. Löve might serve as the maternal donor of the species of Elymus s. l; (2) differentiation of St genome were shown in the species of Elymus s. l. following polyploidy event; (3) divergences within the species might associate with geographic diversity and morphological variability; (4) differences in the levels and patterns of nucleotide diversity of the rbcL gene implied that the St genome lineages in the species of Elymus s. l. have differently evolutionary potentials.     

Comparison of Acetyl-CoA carboxylase 1 (Acc-1) gene diversity among different Triticeae genomes

It has widely been documented that life form and mating system have significant influences on genetic diversity. In the tribe Triticeae, several genera contain both annual and perennial species, whereas other genera comprise strictly annual or perennial species. It was suggested that Triticeae annuals have originated from Triticeae perennials. The present study aims to analyze nucleotide diversity of Acc-1 gene among different Triticeae genomes, and attempts to link effects of life history (annuals and perennials) and mating systems. The nucleotide diversity of 364 Acc-1 sequences in Triticeae species was characterized. The highest estimates of nucleotide diversity values (π = 0.01919, θ = 0.03515) were found for the Ns genome among the genomes analyzed. Nucleotide diversities in the D genome and Ns genome of polyploids are higher than those in respective genomes of diploids, while in the St genome of polyploids, it is lower than that in the St genome of diploids. The averaged π value (0.013705) in the genomes of perennials is more than twice of the value (0.00508) in the genomes of annuals. The averaged π value (0.01323) in the genomes of outcrossing species is two-fold of the value (0.005664) in the genomes of selfer. Our results suggested that the evolutionary history and mating system may play an important role in determining nucleotide diversity of Acc-1 gene in each genome.     

Phylogeny and molecular evolution of the Acc1 gene within the StH genome species in Triticeae (Poaceae)

To estimate the phylogeny and molecular evolution of a single-copy gene encoding plastid acetyl-CoA carboxylase (Acc1) within the StH genome species, two Acc1 homoeologous sequences were isolated from nearly all the sampled StH genome species and were analyzed with those from 35 diploid taxa representing 19 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) the StH genome species from the same areas or neighboring geographic regions are closely related to each other; (2) the Acc1 gene sequences of the StH genome species from North America and Eurasia are evolutionarily distinct; (3) Dasypyrum has contributed to the nuclear genome of Elymus repens and Elymus mutabilis; (4) the StH genome polyploids have higher levels of sequence diversity in the H genome homoeolog than the St genome homoeolog; and (5) the Acc1 sequence may evolve faster in the polyploid species than in the diploids. Our result provides some insight on evolutionary dynamics of duplicate Acc1 gene, the polyploidy speciation and phylogeny of the StH genome species.     

Intergeneric crosses of Psathyrostachys huashanica with Elymus spp. and cytogenetic studies of the hybrids with E. tsukushiensis (Poaceae, Triticeae)

Intergeneric crosses have been made between Psathyrostachys huashanica (2x= 14, NN) and three Elymus species, namely, E. pendulinus (2n = 4x = 28, SSYY), E. ciliaris (2n = 4x = 28, SSYY), and E. tsukushiensis (2n = 6x = 42, SSHHYY). Three accessions of E. tsukushiensis from different localities crossed with P. huashanica produced adult hybrid plants. Although completely sterile, the hybrid plants developed rather vigorously, and were morphologically intermediate between the two parents. Chromosome configuration of 24.62I + 1.51II + 0.03III per cell characterized meioses of the hybrids, but there existed slight variations among different combinations. These results indicate that there is little or no chromosome homoeology between "N" genome of P. huashanica and "S", "H" or "Y" genomes of E. tsukushiensis.     

Comparison of gene flow among species that occur within the same geographic locations versus gene flow among populations within species reveals introgression among several Elymus species

One of the challenges in evolutionary biology is to understand the evolution of speciation with incomplete reproductive isolation as many taxa have continued gene flow both during and after speciation. Comparison of population structure between sympatric and allopatric populations can reveal specific introgression and determine if introgression occurs in a unidirectional or bidirectional manner. Simple sequence repeat markers were used to characterize sympatric and allopatric population structure of plant species, Elymus alaskanus (Scribn. and Merr.) Löve, E. caninus L., E. fibrosus (Schrenk) Tzvel., and E. mutabilis (Drobov) Tzvelev. Our results showed that genetic diversity (H_E) at species level is E. caninus (0.5355) > E. alaskanus (0.4511) > E. fibrosus (0.3924) > E. mutabilis (0.3764), suggesting that E. caninus and E. alaskanus are more variable than E. fibrosus and E. mutabilis. Gene flow between species that occurs within the same geographic locations versus gene flow between populations within species was compared to provide evidence of introgression. Our results indicated that gene flow between species that occur within the same geographic location is higher than that between populations within species, suggesting that gene flow resulting from introgressive hybridization might have occurred among the sympatric populations of these species, and may play an important role in partitioning of genetic diversity among and within populations. The migration rate from E. fibrosus to E. mutabilis is highest (0.2631) among the four species studied. Asymmetrical rates of gene flow among four species were also observed. The findings highlight the complex evolution of these four Elymus species.     

Ubiquity of the St chloroplast genome in St-containing Triticeae polyploids.

Interspecific hybridization occurs between Tritceae species in the grass family (Poaceae) giving rise to allopolyploid species. To examine bias in cytoplasmic DNA inheritance in these hybridizations, the sequence of the 3' end of the chloroplast ndhF gene was compared among 29 allopolyploid Triticeae species containing the St nuclear genome in combination with the H, I, Ns, P, W, Y, and Xm nuclear genomes. These ndhF sequences were also compared with those from diploid or allotetraploid Triticeae species having the H, I, Ns, P, W, St, and Xm genomes. The cpDNA sequences were highly similar among diploid, allotetraploid, allohexaploid, and allooctoploid Triticeae accessions containing the St nuclear genome, with 0-6-nucleotide (nt) substitutions (0-0.8%) occurring between pairs of species. Neighbor-joining analysis of the sequences showed that the ndhF DNA sequences from species containing the St nuclear genome formed a strongly supported clade. The data indicated a strong preference for cpDNA inheritance from the St nuclear genome-containing parent in hybridizations between Triticeae species. This preference was independent of the presence of the H, I, Ns, P, W, and Xm nuclear genomes, the geographic distribution of the species, and the mode of reproduction. The data suggests that hybridizations having the St-containing parent as the female may be more successful.     

An HMM-Based Comparative Genomic Framework for Detecting Introgression in Eukaryotes

One outcome of interspecific hybridization and subsequent effects of evolutionary forces is introgression, which is the integration of genetic material from one species into the genome of an individual in another species. The evolution of several groups of eukaryotic species has involved hybridization, and cases of adaptation through introgression have been already established. In this work, we report on PhyloNet-HMM—a new comparative genomic framework for detecting introgression in genomes. PhyloNet-HMM combines phylogenetic networks with hidden Markov models (HMMs) to simultaneously capture the (potentially reticulate) evolutionary history of the genomes and dependencies within genomes. A novel aspect of our work is that it also accounts for incomplete lineage sorting and dependence across loci. Application of our model to variation data from chromosome 7 in the mouse (Mus musculus domesticus) genome detected a recently reported adaptive introgression event involving the rodent poison resistance gene Vkorc1, in addition to other newly detected introgressed genomic regions. Based on our analysis, it is estimated that about 9% of all sites within chromosome 7 are of introgressive origin (these cover about 13 Mbp of chromosome 7, and over 300 genes). Further, our model detected no introgression in a negative control data set. We also found that our model accurately detected introgression and other evolutionary processes from synthetic data sets simulated under the coalescent model with recombination, isolation, and migration. Our work provides a powerful framework for systematic analysis of introgression while simultaneously accounting for dependence across sites, point mutations, recombination, and ancestral polymorphism.