High-resolution phylogeny for Helianthus (Asteraceae) using the 18S-26S ribosomal DNA external transcribed spacer

The sunflower genus, Helianthus, is recognized widely for the cultivated sunflower H. annuus and scientifically as a model organism for studying diploid and polyploid hybrid speciation, introgression, and genetic architecture. A resolved phylogeny for the genus is essential for the advancement of these scientific areas. In the past, phylogenetic relationships of the perennial species and polyploid hybrids have been particularly difficult to resolve. Using the external transcribed spacer region of the nuclear 18S-26S rDNA region, we reveal for the first time a highly resolved gene tree for Helianthus. Phylogenetic analysis allowed the determination of a monophyletic annual H. sect. Helianthus, a two-lineage polyphyletic H. sect. Ciliares, and the monotypic H. sect. Agrestis, all of which were nested within a large perennial and polyphyletic H. sect. Divaricati. The distribution of perennial polyploids and known annual diploid hybrids on this phylogeny suggested multiple independent hybrid speciation events that gave rise to at least four polyploids and three diploid hybrids. Also provided by this phylogeny was evidence for homoploid hybrid speciation outside H. sect. Helianthus. Finally, previous hypotheses about the secondary chemistry in the genus were tested in a phylogenetic framework to obtain a better understanding of the evolution of these compounds in Helianthus.     

Low intraspecific variation for genomic isolation between hybridizing sunflower species

Barriers to gene flow between species result from selection against foreign linkage blocks in hybrids. When the geographic ranges of taxa meet at multiple locations, the opportunity exists for variation in the genetic architecture of isolating barriers. Hybrid zones between two sunflower species (Helianthus annuus and H. petiolaris) in Nebraska and California exhibited remarkably similar patterns of introgression of mapped molecular markers. Congruence among hybrid zones may result from limited intraspecific variation at loci contributing to isolation and from similar selective effects of alleles in the heterospecific genetic background. The observed consistency of introgression patterns across distantly separated hybrid zones suggests that intrinsic forces predominate in determining hybrid zone dynamics and boundaries between these sunflower species.     

Adaptive introgression of herbivore resistance traits in the weedy sunflower Helianthus annuus

The role of hybridization in adaptive evolution is contentious. While many cases of adaptive trait introgression have been proposed, the relevant traits have rarely been identified, resulting in a lack of clear examples of this process. Here, we examine a purported case of adaptive introgression in which the annual sunflower Helianthus annuus annuus has captured alleles from a congener (Helianthus debilis) to form a stabilized hybrid, Helianthus annuus texanus. We tested the hypotheses that herbivore resistance traits have introgressed from H. debilis to H. annuus and have increased adaptation in the latter. In two common gardens, fitness (estimated by seed production) was on average 55% higher in H. a. texanus than in H. a. annuus. For H. a. texanus, three damage traits (of seven tested) differed significantly from the H. a. annuus parent in one or both sites and were shifted in the direction of the more resistant H. debilis. Natural selection favored H. a. annuusxH. debilis BC(1) hybrids (synthesized to mimic the ancestors of H. a. texanus) with H. debilis-like resistance to seed midges Neolasioptera helianthis and to receptacle/seed feeding Lepidoptera at one or both sites. Assuming similar herbivore pressures in the past, these results suggest that introgression of biotic resistance traits was important in the adaptation of H. annuus to central and southern Texas.     

Complex evolution in Arundinarieae (Poaceae: Bambusoideae): incongruence between plastid and nuclear GBSSI gene phylogenies

The monophyly of tribe Arundinarieae (the temperate woody bamboos) has been unequivocally recovered in previous molecular phylogenetic studies. In a recent phylogenetic study, 10 major lineages in Arundinarieae were resolved based on eight non-coding plastid regions, which conflicted significantly with morphological classifications both at the subtribal and generic levels. Nevertheless, relationships among and within the 10 lineages remain unclear. In order to further unravel the evolutionary history of Arundinarieae, we used the nuclear GBSSI gene sequences along with those of eight plastid regions for phylogenetic reconstruction, with an emphasis on Chinese species. The results of the plastid analyses agreed with previous studies, whereas 13 primary clades revealed in the GBSSI phylogeny were better resolved at the generic level than the plastid phylogeny. Our analyses also revealed many inconsistencies between the plastid DNA and the nuclear GBSSI trees. These results implied that the nuclear genome and the plastid genome had different evolutionary trajectories. The patterns of incongruence suggested that lack of informative characters, incomplete lineage sorting, and/or hybridization (introgression) could be the causes. Seven putative hybrid species were hypothesized, four of which are discussed in detail on the basis of topological incongruence, chromosome numbers, morphology, and distribution patterns, and those taxa probably resulted from homoploid hybrid speciation. Overall, our study indicates that the tribe Arundinarieae has undergone a complex evolution.     

The genomic organization of Ty3/gypsy-like retrotransposons in Helianthus (Asteraceae) homoploid hybrid species

The origin of new diploid, or homoploid, hybrid species is associated with rapid genomic restructuring in the hybrid neospecies. This mode of speciation has been best characterized in wild sunflower species in the genus Helianthus, where three homoploid hybrid species (H. anomalus, H. deserticola, and H. paradoxus) have independently arisen via ancient hybridization events between the same two parental species (H. annuus and H. petiolaris). Most previous work examining genomic restructuring in these sunflower hybrid species has focused on chromosomal rearrangements. However, the origin of all three homoploid hybrid sunflower species also is associated with massive proliferation events of Ty3/gypsy-like retrotransposons in the hybrid species' genomes. We compared the genomic organization of these elements in the parent species and two of the homoploid hybrid species using fluorescence in situ hybridization (FISH). We found a significant expansion of Ty3/gypsy-like retrotransposons confined to the pericentromeric regions of two hybrid sunflower species, H. deserticola and H. paradoxus. In contrast, we detected no significant increase in the frequency or extent of dispersed retrotransposon populations in the hybrid species within the resolution limits of our assay. We discuss the potential role that transposable element proliferation and localization plays in the evolution of homoploid hybrid species.     

THE INTERFACE BETWEEN PHYLOGENETICS AND POPULATION GENETICS: INVESTIGATING GENE TREES,SPECIES TREES,AND POPULATION DYNAMICS IN THE PHYLLOPHAGA FRATERNA SPECIES GROUP

This study uses traditional and contemporary phylogenetic and population genetic analyses to assess the causes of discordance (i.e., lineage sorting and introgression) among mitochondrial and nuclear gene trees for a clade of eastern North American scarab beetles (fraterna species group, genus Phyllophaga). I estimated gene trees using individual and combined analysis of one mitochondrial and two nuclear loci in MrBayes , and inferred a species tree using a hierarchical coalescent approach based on all loci in the program Best . Because hybridization violates the assumptions of Best , I tested for introgression by comparing species monophyly between the mitochondrial and nuclear gene trees based on the prediction that cytoplasmic genomes introgress more readily than nuclear genomes. Haplotype exclusivity was identified using Bayesian tests of monophyly and the genealogical sorting index. I used the results of the phylogenetic analyses and monophyly tests to develop an explicit hypothesis of introgression that could be tested in the program IMa. Results from these analyses provided evidence for introgression across clades within the fraterna group. The tiered analytical approach used in this study demonstrated how the use of multiple methods can identify when assumptions are violated and methods are prone to yield misleading results.     

Chromosomal and Genic Barriers to Introgression in Helianthus

The sexual transfer of genes between taxa possessing different structural karyotypes must involve the passage of genes through a chromosomal sterility barrier. Yet little is known about the effects of structural differences on gene introgression within or adjacent to the rearranged chromosomal fragments or about the patterns of introgression in collinear regions. Here, we employ 197 mapped molecular markers to study the effects of chromosomal structural differences on introgression in backcrossed progeny of the domesticated sunflower, Helianthus annuus, and its karyotypically divergent wild relative, H. petiolaris. Forty percent of the genome from the seven collinear linkages introgressed, whereas only 2.4% of the genome from the 10 rearranged linkages was transferred. Thus, chromosomal rearrangements appear to provide an effective mechanism for reducing or eliminating introgression in rearranged chromosomal segments. On the other hand, observations that 60% of the markers from within the collinear portion of the genome did not introgress suggests that genic factors also resist introgression in Helianthus. That is, selection against H. petiolaris genes in concert with linkage may have reduced or eliminated parts of the genome not protected by structural changes. Thus, barriers to introgression in Helianthus appear to include both chromosomal structural and genic factors.     

Reconstructing the origin of Helianthus deserticola: survival and selection on the desert floor

The diploid hybrid species Helianthus deserticola inhabits the desert floor, an extreme environment relative to its parental species Helianthus annuus and Helianthus petiolaris. Adaptation to the desert floor may have occurred via selection acting on transgressive, or extreme, traits in early hybrids between the parental species. We explored this possibility through a field experiment in the hybrid species' native habitat using H. deserticola, H. annuus, H. petiolaris, and two populations of early-generation (BC(2)) hybrids between the parental species, which served as proxies for the ancestral genotype of the ancient hybrid species. Character expression was evaluated for each genotypic class. Helianthus deserticola was negatively transgressive for stem diameter, leaf area, and flowering date, and the latter two traits are likely to be advantageous in a desert environment. The BC(2) hybrids contained a range of variation that overlapped these transgressive trait means, and an analysis of phenotypic selection revealed that some of the selective pressures on leaf size and flowering date, but not stem diameter, would move the BC(2) population toward the H. deserticola phenotype. Thus, H. deserticola may have originated from habitat-mediated directional selection acting on hybrids between H. annuus and H. petiolaris in a desert environment.     

Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Genome scan of hybridizing sunflowers from Texas (Helianthus annuus and H. debilis) reveals asymmetric patterns of introgression and small islands of genomic differentiation

Although the sexual transfer of genetic material between species (i.e. introgression) has been documented in many groups of plants and animals, genome‐wide patterns of introgression are poorly understood. Is most of the genome permeable to interspecific gene flow, or is introgression typically restricted to a handful of genomic regions? Here, we assess the genomic extent and direction of introgression between three sunflowers from the south‐central USA: the common sunflower, Helianthus annuus ssp. annuus; a near‐endemic to Texas, Helianthus debilis ssp. cucumerifolius; and their putative hybrid derivative, thought to have recently colonized Texas, H. annuus ssp. texanus. Analyses of variation at 88 genetically mapped microsatellite loci revealed that long‐term migration rates were high, genome‐wide and asymmetric, with higher migration rates from H. annuus texanus into the two parental taxa than vice versa. These results imply a longer history of intermittent contact between H. debilis and H. annuus than previously believed, and that H. annuus texanus may serve as a bridge for the transfer of alleles between its parental taxa. They also contradict recent theory suggesting that introgression should predominantly be in the direction of the colonizing species. As in previous studies of hybridizing sunflower species, regions of genetic differentiation appear small, whether estimated in terms of F_ST or unidirectional migration rates. Estimates of recent immigration and admixture were inconsistent, depending on the type of analysis. At the individual locus level, one marker showed striking asymmetry in migration rates, a pattern consistent with tight linkage to a Bateson–Dobzhansky–Muller incompatibility.     

Genealogical portraits of speciation in montane grasshoppers (genus Melanoplus) from the sky islands of the Rocky Mountains

Grasshoppers in the genus Melanoplus have undergone a radiation in the 'sky islands' of western North America, with many species originating during the Pleistocene. Despite their recent origins, phylogenetic analyses indicate that all the species exhibit monophyletic or paraphyletic gene trees. The objectives of this study were to determine whether the monophyletic genealogies are the result of a bottleneck at speciation and to investigate the extent to which the different phylogenetic states of eight species (i.e. monophyletic versus paraphyletic gene trees) can be ascribed to the effects of speciation. A coalescent simulation was used to test for a bottleneck at speciation in each species. The effective population sizes and demographic histories of species were compared across taxa to evaluate the possibility that the paraphyly versus monophyly of the species reflects differential rates of lineage loss rather than speciation mode. While coalescent analyses indicate that the monophyly of Melanoplus species might not be indicative of bottlenecks at speciation, the results suggest that the paraphyletic gene trees may reflect the demography of speciation, involving localized divergences in the ancestral species. With respect to different models of Pleistocene divergence, the data do not support a model of founder-effect speciation but are compatible with divergence in allopatric refugia.     

GENEALOGICAL DISCORDANCE AND PATTERNS OF INTROGRESSION AND SELECTION ACROSS A CRICKET HYBRID ZONE

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus , offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.     

Analysis of a dehydrin encoding gene and its phylogenetic utility in Helianthus

Dehydrins are ubiquitous plant proteins, synthesized in late stages of plant embryo development and following any environmental stress involving dehydration. With the aim to study the evolution of such a stress-responsive gene within Helianthus and to test the possibility of using this gene for phylogenetic studies, fragments of the same dehydrin gene were isolated by PCR and sequenced in 16 wild Helianthus species or subspecies. All isolated sequences included the typical dehydrin domains (Y, S and K), a portion of 3'-UTR and an intron, inserted in the same position within the S domain-encoding region. The number of nucleotide substitutions (both synonymous and nonsynonymous) was calculated keeping separate the different gene regions, and differences occur even among coding domains, indicating that evolutionary constraints act differently on each region. The occurrence of indels and/or insertions was also observed. At the deduced protein level, the calculation of isoelectric point, molecular weight and the percentage of alpha-helix showed a diversification of biochemical properties of this protein between annual and perennial Helianthus species. Phylogenetic trees were built by the maximum-likelihood, maximum-parsimony, and neighbor-joining methods. In all cases the same topology was observed; perennial and annual species form a supported clade, and H. annuus was separated from the other annuals and from perennials. These data support the use of this stress-responsive gene to study the phylogeny of Helianthus.     

The origin of ecological divergence in Helianthus paradoxus (Asteraceae): selection on transgressive characters in a novel hybrid habitat

Diploid hybrid speciation in plants is often accompanied by rapid ecological divergence between incipient neospecies and their parental taxa. One plausible means by which novel adaptation in hybrid lineages may arise is transgressive segregation, that is, the generation of extreme phenotypes that exceed those of the parental lines. Early generation (BC2) hybrids between two wild, annual sunflowers, Helianthus annuus and Helianthus petiolaris, were used to study directional selection on transgressive characters associated with the origin of Helianthus paradoxus, a diploid hybrid species adapted to extremely saline marshes. The BC2 plants descended from a single F1 hybrid backcrossed toward H. petiolaris. The strength of selection on candidate adaptive traits in the interspecific BC2 was measured in natural H. paradoxus salt marsh habitat. Positive directional selection was detected for leaf succulence and Ca uptake, two traits that are known to be important in salt stress response in plants. Strong negative directional selection operated on uptake of Na and correlated elements. A significant decrease in trait correlations over time was observed in the BC2 population for Na and Ca content, suggesting an adaptive role for increased Ca uptake coupled with increased net exclusion of Na from leaves. Patterns of directional selection in BC2 hybrids were concordant with character expression in the natural hybrid species, H. paradoxus, transplanted into the wild. Moreover, the necessary variation for generating the H. paradoxus phenotype existed only in the BC2 population, but not in samples of the two parental species, H. annuus and H. petiolaris. These results are consistent with the hypothesis that transgressive segregation of elemental uptake and leaf succulence contributed to the origin of salt adaptation in the diploid hybrid species H. paradoxus.     

Rampant gene exchange across a strong reproductive barrier between the annual sunflowers, Helianthus annuus and H. petiolaris

Plant species may remain morphologically distinct despite gene exchange with congeners, yet little is known about the genomewide pattern of introgression among species. Here we analyze the effects of persistent gene flow on genomic differentiation between the sympatric sunflower species Helianthus annuus and H. petiolaris. While the species are strongly isolated in testcrosses, genetic distances at 108 microsatellite loci and 14 sequenced genes are highly variable and much lower (on average) than for more closely related but historically allopatric congeners. Our analyses failed to detect a positive association between levels of genetic differentiation and chromosomal rearrangements (as reported in a prior publication) or proximity to QTL for morphological differences or hybrid sterility. However, a significant increase in differentiation was observed for markers within 5 cM of chromosomal breakpoints. Together, these results suggest that islands of differentiation between these two species are small, except in areas of low recombination. Furthermore, only microsatellites associated with ESTs were identified as outlier loci in tests for selection, which might indicate that the ESTs themselves are the targets of selection rather than linked genes (or that coding regions are not randomly distributed). In general, these results indicate that even strong and genetically complex reproductive barriers cannot prevent widespread introgression.     

Hybridization and genome size evolution: timing and magnitude of nuclear DNA content increases in Helianthus homoploid hybrid species

Hybridization and polyploidy can induce rapid genomic changes, including the gain or loss of DNA, but the magnitude and timing of such changes are not well understood. The homoploid hybrid system in Helianthus (three hybrid-derived species and their two parents) provides an opportunity to examine the link between hybridization and genome size changes in a replicated fashion. Flow cytometry was used to estimate the nuclear DNA content in multiple populations of three homoploid hybrid Helianthus species (Helianthus anomalus, Helianthus deserticola, and Helianthus paradoxus), the parental species (Helianthus annuus and Helianthus petiolaris), synthetic hybrids, and natural hybrid-zone populations. Results confirm that hybrid-derived species have 50% more nuclear DNA than the parental species. Despite multiple origins, hybrid species were largely consistent in their DNA content across populations, although H. deserticola showed significant interpopulation differences. First- and sixth-generation synthetic hybrids and hybrid-zone plants did not show an increase from parental DNA content. First-generation hybrids differed in DNA content according to the maternal parent. In summary, hybridization by itself does not lead to increased nuclear DNA content in Helianthus, and the evolutionary forces responsible for the repeated increases in DNA content seen in the hybrid-derived species remain mysterious.     

Examination of hemiplasy, homoplasy and phylogenetic discordance in chromosomal evolution of the Bovidae

Robertsonian chromosomal fusions predominate in shaping the genomes of many species of Bovidae. These and other cytogenetic data (from 52 taxa representing 51 species and 9 tribes of Bovidae) were (i) examined for usefulness in defining phylogenetic relationships and (ii) subsequently mapped to a consensus tree based on mitochondrial and nuclear DNA gene sequences with divergence dates of the corresponding species calculated from cytochrome b sequences. This permitted persistence time estimates for the various rearrangements. The chromosomal data resulted in an unsupported higher-level topology, but with recognition of the monophyly of some genera and tribes within Bovidae. The distribution and temporal spread of character states on the species tree is suggestive of a restricted role for hemiplasy (the retention of an ancestral chromosomal polymorphism through multiple speciation events) and for introgression (resulting from secondary contact among taxa), processes that can potentially lead to phylogenetic discordance. We conclude that the most probable interpretation for these data is that genuine karyotypic homoplasy predominates, but that hemiplasy (and/or introgression) is a realistic hypothesis for the observed patterns of several shared characters in Bovidae.     

Distribution of Ty3-gypsy- and Ty1-copia-like DNA sequences in the genus Helianthus and other Asteraceae.

Two repeated DNA sequences, pHaS13 and pHaS211, which revealed similarity to the int gene of Ty3-gypsy retrotransposons and the RNAse-H gene of Ty1-copia retroelements, respectively, were surveyed in Asteraceae species and within the genus Helianthus. Southern analysis of the genome of selected Asteraceae that belong to different tribes showed that pHaS13- and pHaS211-related subfamilies of gypsy- and copia-like retroelements are highly redundant only in Helianthus and, to a lesser extent, in Tithonia, a Helianthus strict relative. However, under low stringency posthybridization washes, bands were observed in almost all the other Asteraceae tested when pHaS13 was used as a probe, and in several species when pHaS211 was hybridized. FISH analysis of pHaS13 or pHaS211 probes was performed in species in which labelling was observed in Southern hybridizations carried out under high stringency conditions (Helianthus annuus, Tithonia rotundifolia, Ageratum spp., Leontopodium spp., Senecio vulgaris for pHaS13, and H. annuus, Tithonia rotundifolia, and S. vulgaris for pHaS211). Scattered labelling was observed over all metaphase chromosomes, indicating a large dispersal of both Ty3-gypsy- and Ty1-copia-like retroelements. However, preferential localization of Ty3-gypsy-like sequences at centromeric chromosome regions was observed in all of the species studies but one, even in species in which pHaS13-related elements are poorly represented. Ty1-copia-like sequences showed preferential localization at the chromosome ends only in H. annuus. To study the evolution of gypsy- and copia-like retrotransposons in Helianthus, cladograms were built based on the Southern blot hybridization patterns of pHaS13 or pHaS211 sequences to DNA digests of several species of this genus. Both cladograms agree in splitting the genomes studied into annuals and perennials. Differences that occurred within the clades of perennial and annual species between gypsy- and copia-like retroelements indicated that these retrotransposons were differentially active during Helianthus speciation, suggesting that the evolution of the 2 retroelement families was, within limits, independent.     

Genomic markers reveal introgressive hybridization in the Indo-West Pacific mangroves: a case study

Biodiversity of mangrove ecosystems is difficult to assess, at least partly due to lack of genetic verification of morphology-based documentation of species. Natural hybridization, on the one hand, plays an important role in evolution as a source of novel gene combinations and a mechanism of speciation. However, on the other hand, recurrent introgression allows gene flow between species and could reverse the process of genetic differentiation among populations required for speciation. To understand the dynamic evolutionary consequences of hybridization, this study examines genomic structure of hybrids and parental species at the population level. In the Indo-West Pacific, Bruguiera is one of the dominant mangrove genera and species ranges overlap extensively with one another. Morphological intermediates between sympatric Bruguiera gymnorrhiza and Bruguiera sexangula have been reported as a variety of B. sexangula or a new hybrid species, B. × rhynchopetala. However, the direction of hybridization and extent of introgression are unclear. A large number of species-specific inter-simple sequence repeat (ISSR) markers were found in B. gymnorrhiza and B. sexangula, and the additive ISSR profiling of B. × rhynchopetala ascertained its hybrid status and identified its parental origin. The varying degree of scatterness among hybrid individuals in Principal Coordinate Analysis and results from NewHybrids analysis indicate that B. × rhynchopetala comprises different generations of introgressants in addition to F(1)s. High genetic relatedness between B. × rhynchopetala and B. gymnorrhiza based on nuclear and chloroplast sequences suggests preferential hybrid backcrosses to B. gymnorrhiza. We conclude that B. × rhynchopetala has not evolved into an incipient hybrid species, and its persistence can be explained by recurrent hybridization and introgression. Genomic data provide insights into the hybridization dynamics of mangrove plants. Such information can assist in biodiversity assessment by helping detect novel taxa and/or define species boundaries.