Parental Lineages Research Articles

Genetic and karyotype divergence between parents affect clonality and sterility in hybrids.

Asexual reproduction can be triggered by interspecific hybridization, but its emergence is supposedly rare, relying on exceptional combinations of suitable genomes. To examine how genomic and karyotype divergence between parental lineages affect the incidence of asexual gametogenesis, we experimentally hybridized fishes (Cobitidae) across a broad phylogenetic spectrum, assessed by whole exome data. Gametogenic pathways generally followed a continuum from sexual reproduction in hybrids between closely related evolutionary lineages to sterile or inviable crosses between distant lineages. However, most crosses resulted in a combination of sterile males and asexually reproducing females. Their gametes usually experienced problems in chromosome pairing, but females also produced a certain proportion of oocytes with premeiotically duplicated genomes, enabling their development into clonal eggs. Interspecific hybridization may thus commonly affect cell cycles in a specific way, allowing the formation of unreduced oocytes. The emergence of asexual gametogenesis appears tightly linked to hybrid sterility and constitutes an inherent part of the extended speciation continuum.

Genetic and karyotype divergence between parents affect clonality and sterility in hybrids

Asexual reproduction can be triggered by interspecific hybridization, but its emergence is supposedly rare, relying on exceptional combinations of suitable genomes. To examine how genomic and karyotype divergence between parental lineages affect the incidence of asexual gametogenesis, we experimentally hybridized fishes (Cobitidae) across a broad phylogenetic spectrum, assessed by whole exome data. Gametogenic pathways generally followed a continuum from sexual reproduction in hybrids between closely related evolutionary lineages to sterile or inviable crosses between distant lineages. However, most crosses resulted in a combination of sterile males and asexually reproducing females. Their gametes usually experienced problems in chromosome pairing, but females also produced a certain proportion of oocytes with premeiotically duplicated genomes, enabling their development into clonal eggs. Interspecific hybridization may thus commonly affect cell cycles in a specific way, allowing the formation of unreduced oocytes. The emergence of asexual gametogenesis appears tightly linked to hybrid sterility and constitutes an inherent part of the extended speciation continuum.

Integrating morphological, molecular and cytogenetic data for F2 sea turtle hybrids diagnosis revealed balanced chromosomal sets.

Hybridization could be considered part of the evolutionary history of many species. The hybridization among sea turtle species on the Brazilian coast is atypical and occurs where nesting areas and reproductive seasons overlap. Integrated analysis of morphology and genetics is still scarce, and there is no evidence of the parental chromosome set distribution in sea turtle interspecific hybrids. In this study, chromosome markers previously established for pure sea turtle species were combined with morphological and molecular analyses aiming to recognize genetic composition and chromosome sets in possible interspecific hybrids initially identified by mixed morphology. The data showed that one hybrid could be an F2 individual among Caretta caretta × Eretmochelys imbricata × Chelonia mydas, and another is resulting from backcross between C. caretta × Lepidochelys olivacea. Native alleles of different parental lineages were reported in the hybrids, and, despite this, it was verified that the hybrid chromosome sets were still balanced. Thus, how sea turtle hybridism can affect genetic features in the long term is a concern, as the implications of the crossing-over in hybrid chromosomal sets and the effects on genetic function are still unpredictable.

Broadly Reactive SARS-CoV-2-Specific T-Cell Response and Participation of Memory B and T Cells in Patients with Omicron COVID-19 Infection.

January 2022 onward, India witnessed a sudden increase in Omicron COVID-19 infections, having a mild course that prompted us to identify the key host factors/immune molecules modulating disease course/outcomes. The current study evaluated the percentages of lymphocyte subsets by flowcytometry, SARS-CoV-2 specific T-cell immune response by ELISPOT, estimation of plasma cytokine/chemokine levels on a Bio-plex Multiplex Immunoassay System and anti-SARS-CoV-2 IgG levels by enzyme-linked immunosorbent assay in 19 mild Omicron infected patients, 45 mild SARS-CoV-2 (2020) patients and 36 uninfected controls from India. Natural killer cells, B and memory B cells were high in vaccinated and total Omicron-infected patients groups compared to the mild SARS-CoV-2 (2020) patient group, while CD8+ T cells were high in total Omicron-infected patients group compared to the uninfected control group (p < 0.05 each). Omicron-infected patients had T-cell response against SARS-CoV-2 whole virus, S1 proteins (wild type and delta variant) in 10 out of 17 (59%), 10 out of 17 (59%), and 8 out of 17 (47%), respectively. The current study of Omicron-infected patients elucidates broadly reactive antibody, T-cell response, and participation of memory B and T cells induced by vaccination/natural infection. The limited effect of Omicron's mutations on T-cell response is suggestive of protection from severity. Pro-inflammatory IL-6, IFN-γ, chemokines CCL-2, CCL-3, CCL-4, CCL-5, and IL-8 as potential biomarkers of Omicron infection may have future diagnostic importance. The cellular immune response data in Omicron-infected patients with parental Omicron lineage could serve as a starting point to define the readouts of protective immunity against circulating Omicron subvariants.

An effective artificial feeder for breeding Aedes aegypti (Diptera: Culicidae)

Mosquitoes, like other hematophagous insects, need to feed on blood to produce eggs. For laboratory colony maintenance, live animals are generally used, while various artificial feeders with different characteristics have been developed for experimental assays. In this study, a previously developed and pathogen-tested glass artificial feeder was selected and evaluated for rearing Aedes aegypti (Linnaeus). The device utilizes bovine gut as a membrane and human blood as a food source, maintaining a constant temperature of 37 °C through a thermostatic bath that enables water circulation around the feeder. Egg hatching, feeding rates, and fertility were analyzed for three cohorts of a parental lineage (P) and their respective filial lineages (F1) over a period of 4-5 weeks. The results demonstrated the efficacy of this methodology for rearing an Ae. aegypti colony up to the second generation, with average egg hatching rates of P= 0.92 and F1= 0.80, average feeding rates of P= 0.54 and F1= 0.73, and egg quantities obtained of P= 1809 and F1 = 8579.

Clonal diversity and zoonotic potential of MDR Escherichia coli isolated from poultry at different age intervals

ABSTRACT 1. A pool of 480 E. coli isolates of poultry (broilers and ducks) representing different time intervals (0, 10, 20 and 30 days) was selected for ribotyping and used to determine polymorphism of 16-23S ribosomal RNA intergenic space. All the isolates were multidrug-resistant (MDR). 2. Out of these, 10 isolates were tested for MultiLocus Sequence Typing (MLST) among which novel allelic combinations and therefore new sequence types were identified in seven isolates. 3. This work showed the changes in E. coli strains structure at farm level and individual bird level in host species raised on organised farms with similar parental lineage and environmental housing. The statistical results showed that the structure of variation is very different by farm, supporting a strong effect of location, which confirms the temporal clustering. 4. There were significant differences between E. coli strains in chickens and ducks, indicating host specificity of the E. coli strains. 5. Some of the pathogenic E. coli strains found using MLST belonged to ST735, ST2796 and a pandemic clone ST752 of ST10 clonal complex. The results strongly suggested the clonal expansion and establishment of specific MDR clones that have zoonotic relevance.

Integrative Genome-Based Survey of the SARS-CoV-2 Omicron XBB.1.16 Variant.

The XBB.1.16 SARS-CoV-2 variant, also known as Arcturus, is a recent descendant lineage of the recombinant XBB (nicknamed Gryphon). Compared to its direct progenitor, XBB.1, XBB.1.16 carries additional spike mutations in key antigenic sites, potentially conferring an ability to evade the immune response compared to other circulating lineages. In this context, we conducted a comprehensive genome-based survey to gain a detailed understanding of the evolution and potential dangers of the XBB.1.16 variant, which became dominant in late June. Genetic data indicates that the XBB.1.16 variant exhibits an evolutionary background with limited diversification, unlike dangerous lineages known for rapid changes. The evolutionary rate of XBB.1.16, which amounts to 3.95 × 10-4 subs/site/year, is slightly slower than that of its direct progenitors, XBB and XBB.1.5, which have been circulating for several months. A Bayesian Skyline Plot reconstruction suggests that the peak of genetic variability was reached in early May 2023, and currently, it is in a plateau phase with a viral population size similar to the levels observed in early March. Structural analyses indicate that, overall, the XBB.1.16 variant does not possess structural characteristics markedly different from those of the parent lineages, and the theoretical affinity for ACE2 does not seem to change among the compared variants. In conclusion, the genetic and structural analyses of SARS-CoV-2 XBB.1.16 do not provide evidence of its exceptional danger or high expansion capability. Detected differences with previous lineages are probably due to genetic drift, which allows the virus constant adaptability to the host, but they are not necessarily connected to a greater danger. Nevertheless, continuous genome-based monitoring is essential for a better understanding of its descendants and other lineages.

Spike protein mutations and structural insights of pangolin lineage B.1.1.25 with implications for viral pathogenicity and ACE2 binding affinity

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of COVID -19, is constantly evolving, requiring continuous genomic surveillance. In this study, we used whole-genome sequencing to investigate the genetic epidemiology of SARS-CoV-2 in Bangladesh, with particular emphasis on identifying dominant variants and associated mutations. We used high-throughput next-generation sequencing (NGS) to obtain DNA sequences from COVID-19 patient samples and compared these sequences to the Wuhan SARS-CoV-2 reference genome using the Global Initiative for Sharing All Influenza Data (GISAID). Our phylogenetic and mutational analyzes revealed that the majority (88%) of the samples belonged to the pangolin lineage B.1.1.25, whereas the remaining 11% were assigned to the parental lineage B.1.1. Two main mutations, D614G and P681R, were identified in the spike protein sequences of the samples. The D614G mutation, which is the most common, decreases S1 domain flexibility, whereas the P681R mutation may increase the severity of viral infections by increasing the binding affinity between the spike protein and the ACE2 receptor. We employed molecular modeling techniques, including protein modeling, molecular docking, and quantum mechanics/molecular mechanics (QM/MM) geometry optimization, to build and validate three-dimensional models of the S_D614G-ACE2 and S_P681R-ACE2 complexes from the predominant strains. The description of the binding mode and intermolecular contacts of the referenced systems suggests that the P681R mutation may be associated with increased viral pathogenicity in Bangladeshi patients due to enhanced electrostatic interactions between the mutant spike protein and the human ACE2 receptor, underscoring the importance of continuous genomic surveillance in the fight against COVID -19. Finally, the binding profile of the S_D614G-ACE2 and S_P681R-ACE2 complexes offer valuable insights to deeply understand the binding site characteristics that could help to develop antiviral therapeutics that inhibit protein–protein interactions between SARS-CoV-2 spike protein and human ACE2 receptor.

Phylogenomics insights into gene evolution, rapid species diversification, and morphological innovation of the apple tribe (Maleae, Rosaceae).

Maleae is one of the most widespread tribes of Rosaceae and includes several important fruit crops and ornamental plants. We used nuclear genes from 62 transcriptomes/genomes, including 26 newly generated transcriptomes, to reconstruct a well-supported phylogeny and study the evolution of fruit and leaf morphology and the possible effect of whole genome duplication (WGD). Our phylogeny recovered 11 well-supported clades and supported the monophyly of most genera (except Malus, Sorbus, and Pourthiaea) with at least two sampled species. A WGD was located to the most recent common ancestor (MRCA) of Maleae and dated to c. 54 million years ago (Ma) near the Early Eocene Climatic Optimum, supporting Gillenieae (x = 9) being a parental lineage of Maleae (x = 17) and including duplicate regulatory genes related to the origin of the fleshy pome fruit. Whole genome duplication-derived paralogs that are retained in specific lineages but lost in others are predicted to function in development, metabolism, and other processes. An upshift of diversification and innovations of fruit and leaf morphologies occurred at the MRCA of the Malinae subtribe, coinciding with the Eocene-Oligocene transition (c. 34 Ma), following a lag from the time of the WGD event. Our results provide new insights into the Maleae phylogeny, its rapid diversification, and morphological and molecular evolution.

Ancient hybridization in Curcuma (Zingiberaceae)-Accelerator or brake in lineage diversifications?

Hybridization is a widespread phenomenon in the evolution of plants and exploring its role is crucial to understanding diversification processes of many taxonomic groups. Recently, more attention is focused on the role of ancient hybridization that has repeatedly been shown as triggers of evolutionary radiation, although in some cases, it can prevent further diversification. The causes, frequency, and consequences of ancient hybridization remain to be explored. Here, we present an account of several events of ancient hybridization in turmeric, the economically important plant genus Curcuma (Zingiberaceae), which harbors about 130 known species. We analyzed 1094 targeted low-copy genes and plastomes obtained by next-generation sequencing of 37 species of Curcuma, representing the known genetic diversity and spanning the geographical distribution of the genus. Using phylogenetic network analysis, we show that the entire genus Curcuma as well as its most speciose lineage arose via introgression from the genus Pyrgophyllum and one of the extinct lineages, respectively. We also document a single event of ancient hybridization, with C. vamana as a product, that represents an evolutionary dead end. We further discuss distinct circumstances of those hybridization events that deal mainly with (in)congruence in chromosome counts of the parental lineages.

Environmentally independent selection for hybrids between divergent freshwater stickleback lineages in semi-natural ponds.

Hybridization following secondary contact of genetically divergent populations can influence the range expansion of invasive species, though specific outcomes depend on the environmental dependence of hybrid fitness. Here, using two genetically and ecologically divergent threespine stickleback lineages that differ in their history of freshwater colonization, we estimate fitness variation of parental lineages and hybrids in semi-natural freshwater ponds with contrasting histories of nutrient loading. In our experiment, we found that fish from the older freshwater lineage (Lake Geneva) and hybrids outperformed fish from the younger freshwater lineage (Lake Constance) in terms of both growth and survival, regardless of the environmental context of our ponds. Across all ponds, hybrids exhibited the highest survival. Although wild-caught adult populations differed in their functional and defence morphology, it is unclear which of these traits underlie the fitness differences observed among juveniles in our experiment. Overall, our work suggests that when hybrid fitness is insensitive to environmental conditions, as observed here, introgression may promote population expansion into unoccupied habitats and accelerate invasion success.

Exploring the Possible Role of Hybridization in the Evolution of Photosynthetic Pathways in Flaveria (Asteraceae), the Prime Model of C4 Photosynthesis Evolution

Flaveria (Asteraceae) is the prime model for the study of C4 photosynthesis evolution and seems to support a stepwise acquisition of the pathway through C3-C4 intermediate phenotypes, still existing in Flaveria today. Molecular phylogenies of Flaveria based on concatenated data matrices are currently used to reconstruct the complex sequence of trait shifts during C4 evolution. To assess the possible role of hybridization in C4 evolution in Flaveria, we re-analyzed transcriptome data of 17 Flaveria species to infer the extent of gene tree discordance and possible reticulation events. We found massive gene tree discordance as well as reticulation along the backbone and within clades containing C3-C4 intermediate and C4-like species. An early hybridization event between two C3 species might have triggered C4 evolution in the genus. The clade containing all C4 species plus the C4-like species F. vaginata and F. palmeri is highly supported in our phylogenetic analyses, but it might be of hybrid origin involving F. angustifolia and F. sonorensis (both C3-C4 intermediate) as parental lineages. Hybridization seems to be a driver of C4 evolution in Flaveria and likely promoted the fast acquisition of C4 traits. This new insight can be used in further exploring C4 evolution and can inform C4 bioengineering efforts.

SARS-CoV-2 Recombinants: Genomic Comparison between XBF and Its Parental Lineages.

Recombination events are very common and represent one of the primary drivers of RNA virus evolution. The XBF SARS-CoV-2 lineage is one of the most recently generated recombinants during the COVID-19 pandemic. It is a recombinant of BA.5.2.3 and BA.2.75.3, both descendants of lineages that caused many concerns (BA.5 and BA.2.75, respectively). Here, we performed a genomic survey focused on comparing the recombinant XBF with its parental lineages to provide a comprehensive assessment of the evolutionary potential, epidemiological trajectory, and potential risks. Genetic analyses indicated that although XBF initially showed the typical expansion depicted by a steep curve, causing several concerns, currently there is no indication of significant expansion potential or a contagion rate surpassing that of other currently active or previously prevalent lineages. BSP indicated that the peak has been reached around 19 October 2022 and then the genetic variability suffered slight oscillations until early 5 March 2023 when the population size reduced for the last time starting its last plateau that is still lasting. Structural analyses confirmed its reduced potential, also indicating that properties of NTDs and RBDs of XBF and its parental lineages present no significant difference. Of course, cautionary measures must still be taken and genome-based monitoring remains the best tool for detecting any important changes in viral genome composition.

Breakthrough SARS-CoV-2 Omicron Variant in Individuals Primed with Heterologous Vaccines Enhances Inhibition Performance of Neutralizing Antibody to BA.2 Parental Lineage.

This study aims to analyze the neutralization ability against Omicron parental variants in five clusters of individuals with different Coronavirus disease (COVID-19) immunity backgrounds, including individuals receiving a homologous or heterologous vaccine without prior infection, recovered patients with homologous or heterologous vaccination, and recovery patients without vaccination. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) surrogate virus neutralization assay was performed on serum samples. Spearman correlation analysis showed that the percent inhibition against Omicron B.1.1.529 and BA.2 was significantly related to the period of serum collection (r = 0.730 and 0.787, p < 0.001, respectively). Very strong correlation between percent inhibition of neutralizing antibody against Omicron B.1.1.529 and BA.2 variants (rs = 0.973, p < 0.001) was also observed. The neutralizing activity of the sera from recovery patients receiving homologous and heterologous vaccine against the wild-type, B.1.1.529, and BA.2 Omicron variants was significantly higher (p < 0.001) than that of recovery patients without vaccination. This study robustly showed that the breakthrough SARS-CoV-2 Omicron variant in individuals who received homologous and heterologous vaccines had a high level of neutralizing activity against B.1.1.529 and BA.2 parental lineage of XBB subvariants. Therefore, the next-generation COVID-19 vaccine against emerging variants is needed to improve resilience against ongoing variants, particularly for persons who have never been infected.

Survival rates and mortality risks of Plecturocebus cupreus at the California National Primate Research Center.

This article describes survivorship and explores factors affecting mortality risks in a captive colony of coppery titi monkeys (Plecturocebus cupreus) housed at the California National Primate Research Center (CNPRC), at UC Davis, in Davis, CA. We analyzed data collected on individuals since the colony's creation in the 1960s, with a sample of 600 animals with partially complete information (date of birth, age at death, body mass, parental lineage). We used three methods: (1) Kaplan-Meier regressions followed by a log-rank test to compare survival in male and female titi monkeys, (2) a breakpoint analysis to identify shifts in the survival curves, and (3) Cox regressions to test the effect of body mass change, parental pair tenure, and parental age on mortality risk. We found that males tend to have a longer median lifespan than females (14.9 and 11.4 years; p = 0.094) and that survival decreases earlier in males than in females during adulthood (9.8 and 16.2 years). A body mass loss of 10% from adulthood to the time of death led to a 26% higher risk of dying (p < 0.001) as compared to an individual with stable body mass. We found no evidence of sociobiological factors on mortality risks (parental age, parental pair tenure), but an exploratory analysis suggested that a higher rate of offspring conceptions increases mortality risks. This description of factors influencing survival and mortality in titi monkeys is a first step toward understanding aging in this species to consider titi monkeys as a primate model for socioemotional aging.

Demographic responses of hybridizing cinquefoils to changing climate in the Colorado Rocky Mountains.

Hybridization between taxa generates new pools of genetic variation that can lead to different environmental responses and demographic trajectories over time than seen in parental lineages. The potential for hybrids to have novel environmental tolerances may be increasingly important in mountainous regions, which are rapidly warming and drying due to climate change. Demographic analysis makes it possible to quantify within- and among-species responses to variation in climate and to predict population growth rates as those conditions change. We estimated vital rates and population growth in 13 natural populations of two cinquefoil taxa (Potentilla hippiana and P. pulcherrima) and their hybrid across elevation gradients in the Southern Rockies. Using three consecutive years of environmental and demographic data, we compared the demographic responses of hybrid and parental taxa to environmental variation across space and time. All three taxa had lower predicted population growth rates under warm, dry conditions. However, the magnitude of these responses varied among taxa and populations. Hybrids had consistently lower predicted population growth rates than P. hippiana. In contrast, hybrid performance relative to P. pulcherrima varied with population and climate, with the hybrid maintaining relatively stable growth rates while populations of P. pulcherrima shrank under warm, dry conditions. Our findings demonstrate that hybrids in this system are neither intrinsically unfit nor universally more vigorous than parents, suggesting that the demographic consequences of hybridization are context-dependent. Our results also imply that shifts to warmer and drier conditions could have particularly negative repercussions for P. pulcherrima, which is currently the most abundant taxon in the study area, possibly as a legacy of more favorable historical climates. More broadly, the distributions of these long-lived taxa are lagging behind their demographic trajectories, such that the currently less common P. hippiana could become the most abundant of the Potentilla taxa as this region continues to warm and dry.

Genetic features of Sri Lankan elephant, Elephas maximus maximus Linnaeus revealed by high throughput sequencing of mitogenome and ddRAD-seq.

Elephas maximus maximus Linnaeus, the Sri Lankan subspecies is the largest and the darkest among Asian elephants. Patches of depigmented areas with no skin color on the ears, face, trunk, and belly morphologically differentiate it from the others. The elephant population in Sri Lanka is now limited to smaller areas and protected under Sri Lankan law. Despite its ecological and evolutionary importance, the relationship between Sri Lankan elephants and their phylogenetic position among Asian elephants remains controversial. While identifying genetic diversity is the key to any conservation and management strategies, limited data is currently available. To address such issues, we analyzed 24 elephants with known parental lineages with high throughput ddRAD-seq. The mitogenome suggested the coalescence time of the Sri Lankan elephant at ~0.2 million years, and sister to Myanmar elephants supporting the hypothesis of the movement of elephants in Eurasia. The ddRAD-seq approach identified 50,490 genome-wide SNPs among Sri Lankan elephants. The genetic diversity within Sri Lankan elephants assessed with identified SNPs suggests a geographical differentiation resulting in three main clusters; north-eastern, mid-latitude, and southern regions. Interestingly, though it was believed that elephants from the Sinharaja rainforest are of an isolated population, the ddRAD-based genetic analysis clustered it with the north-eastern elephants. The effect of habitat fragmentation on genetic diversity could be further assessed with more samples with specific SNPs identified in the current study.

Hybrid origin of a primate, the gray snub-nosed monkey.

Hybridization is widely recognized as promoting both species and phenotypic diversity. However, its role in mammalian evolution is rarely examined. We report historical hybridization among a group of snub-nosed monkeys (Rhinopithecus) that resulted in the origin of a hybrid species. The geographically isolated gray snub-nosed monkey Rhinopithecus brelichi shows a stable mixed genomic ancestry derived from the golden snub-nosed monkey (Rhinopithecus roxellana) and the ancestor of black-white (Rhinopithecus bieti) and black snub-nosed monkeys (Rhinopithecus strykeri). We further identified key genes derived from the parental lineages, respectively, that may have contributed to the mosaic coat coloration of R. brelichi, which likely promoted premating reproductive isolation of the hybrid from parental lineages. Our study highlights the underappreciated role of hybridization in generating species and phenotypic diversity in mammals.

Comparative genomics reveals the hybrid origin of a macaque group.

Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.

Transgressive segregation in mating traits drives hybrid speciation.

Hybridization can rapidly generate novel genetic variation, which can promote ecological speciation by creating novel adaptive phenotypes. However, it remains unclear how hybridization, creating novel mating phenotypes (e.g., mating season, genitalia shapes, sexual displays, mate preferences), affects speciation especially when the phenotypes do not confer adaptive advantages. Here, based on individual-based evolutionary simulations, we propose that transgressive segregation of mating traits can drive incipient hybrid speciation. Simulations demonstrated that incipient hybrid speciation occurred most frequently when the hybrid population received moderate continued immigration from parental lineages causing recurrent episodes of hybridization. Recurrent hybridization constantly generated genetic variation, which promoted the rapid stochastic evolution of mating phenotypes in a hybrid population. The stochastic evolution continued until a novel mating phenotype came to dominate the hybrid population, which reproductively isolates the hybrid population from parental lineages. However, too frequent hybridization rather hindered the evolution of reproductive isolation by inflating the variation of mating phenotypes to produce phenotypes allowing mating with parental lineages. Simulations also revealed conditions for the long-term persistence of hybrid species after their incipient emergence. Our results suggest that recurrent transgressive segregation of mating phenotypes can offer a plausible explanation for hybrid speciation and radiations that involved little adaptive ecological divergence.