Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis.

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Cryptic diversity in southern African kelp.

The southern coast of Africa is one of the few places in the world where water temperatures are predicted to cool in the future. This endemism-rich coastline is home to two sister species of kelps of the genus Ecklonia maxima and Ecklonia radiata, each associated with specific thermal niches, and occuring primarily on opposite sides of the southern tip of Africa. Historical distribution records indicate that E. maxima has recently shifted its distribution ~ 70 km eastward, to sites where only E. radiata was previously reported. The contact of sister species with contrasting thermal affinities and the occurrence of mixed morphologies raised the hypothesis that hybridization might be occurring in this contact zone. Here we describe the genetic structure of the genus Ecklonia along the southern coast of Africa and investigate potential hybridization and cryptic diversity using a combination of nuclear microsatellites and mitochondrial markers. We found that both species have geographically discrete genetic clusters, consistent with expected phylogeographic breaks along this coastline. In addition, depth-isolated populations were found to harbor unique genetic diversity, including a third Ecklonia lineage. Mito-nuclear discordance and high genetic divergence in the contact zones suggest multiple hybridization events between Ecklonia species. Discordance between morphological and molecular identification suggests the potential influence of abiotic factors leading to convergent phenotypes in the contact zones. Our results highlight an example of cryptic diversity and hybridization driven by contact between two closely related keystone species with contrasting thermal affinities.

Chromosome stability of synthetic Triticum turgidum-Aegilops umbellulata hybrids.

BACKGROUND: Unreduced gamete formation during meiosis plays a critical role in natural polyploidization. However, the unreduced gamete formation mechanisms in Triticum turgidum-Aegilops umbellulata triploid F1 hybrid crosses and the chromsome numbers and compostions in T. turgidum-Ae. umbellulata F2 still not known. RESULTS: In this study, 11 T.turgidum-Ae. umbellulata triploid F1 hybrid crosses were produced by distant hybridization. All of the triploid F1 hybrids had 21 chromosomes and two basic pathways of meiotic restitution, namely first-division restitution (FDR) and single-division meiosis (SDM). Only FDR was found in six of the 11 crosses, while both FDR and SDM occurred in the remaining five crosses. The chromosome numbers in the 127 selfed F2 seeds from the triploid F1 hybrid plants of 10 crosses (no F2 seeds for STU 16) varied from 35 to 43, and the proportions of euploid and aneuploid F2 plants were 49.61% and 50.39%, respectively. In the aneuploid F2 plants, the frequency of chromosome loss/gain varied among genomes. The chromosome loss of the U genome was the highest (26.77%) among the three genomes, followed by that of the B (22.83%) and A (11.81%) genomes, and the chromosome gain for the A, B, and U genomes was 3.94%, 3.94%, and 1.57%, respectively. Of the 21 chromosomes, 7U (16.54%), 5 A (3.94%), and 1B (9.45%) had the highest loss frequency among the U, A, and B genomes. In addition to chromosome loss, seven chromosomes, namely 1 A, 3 A, 5 A, 6 A, 1B, 1U, and 6U, were gained in the aneuploids. CONCLUSION: In the aneuploid F2 plants, the frequency of chromosome loss/gain varied among genomes, chromsomes, and crosses. In addition to variations in chromosome numbers, three types of chromosome translocations including 3UL·2AS, 6UL·1AL, and 4US·6AL were identified in the F2 plants. Furthermore, polymorphic fluorescence in situ hybridization karyotypes for all the U chromosomes were also identified in the F2 plants when compared with the Ae. umbellulata parents. These results provide useful information for our understanding the naturally occurred T. turgidum-Ae. umbellulata amphidiploids.

Lifetime milk production of Holstein cattle in the humid tropics compared to Holstein-Gyr and Holstein-Brahman crosses.

Crossbred cattle are commonly used for milk production in the tropics, combining the potential benefits of pure breeds with the heterosis effects of the offspring. However, no comprehensive assessment of lifetime productivity for crossbred versus purebred cattle in low-altitude tropical environments has been carried out. The present study compares the lifetime productivity of purebred Holstein (HO, n = 17,269), Gyr (GY4, n = 435), and Brahman (BR4, n = 622) with crossbreds Gyr × Holstein (GY × HO, n = 5521) and Brahman×Holstein (BR × HO, n = 5429) cows from dairy farms located in low and medium altitude tropical regions in Costa Rica. The production traits of interest were age at first calving (AFC), days open (DO), milk production per lactation (TMP), lactation length (LLEN), age at culling (ACUL), and number of lactations (NLAC). Estimates of heterosis were also calculated. The AFC for GY × HO crosses (33-34 months) was not significantly different (p > .05) from HO (33.8 months). For BR × HO crosses, a significant (p < .05) decrease in AFC (BR3HO1 35.6 months, BR2HO2 34.5 months, and BR1H03 33.3 months) was observed as the fraction of HO breed increased. Estimates of heterosis for AFC were favourable for both crosses, of a magnitude close to 3%. The DO for F1 crosses (GY2HO2 94 days; BR2HO2 96 days) was significantly (p < .05) lower than HO (123 days). Estimates of heterosis for DO were also favourable and above 15% for both crosses. The TMP and LLEN were higher for HO (TMP = 5003 kg; LLEN = 324 days) compared with GY × HO (TMP = 4428 to 4773 kg; LLEN = 298 to 312 days) and BR × HO (TMP = 3950 to 4761 kg; LLEN = 273 to 313 days) crosses. Heterosis for TMP was favourable but low for both crosses, with a magnitude below 3.0%. The NLAC for HO (4.6 lactations) was significantly (p < .05) lower than F1 (GY2HO2, 5.8 lactations; BR2HO2, 5.4 lactations). Heterosis for NLAC was above 6.0% for both crosses. Overall, estimates of lifetime income over feed costs per cow on average were USD 2637 (30.3%) and USD 734 (8.4%) higher in F1 GY × HO and BR × HO, respectively, compared to HO. In conclusion, crossbred animals, specifically those with Gyr and Brahman genetics, extend the productive lifespan, increasing economic returns.

Meiotic drive against chromosome fusions in butterfly hybrids.

Species frequently differ in the number and structure of chromosomes they harbor, but individuals that are heterozygous for chromosomal rearrangements may suffer from reduced fitness. Chromosomal rearrangements like fissions and fusions can hence serve as a mechanism for speciation between incipient lineages, but their evolution poses a paradox. How can rearrangements get fixed between populations if heterozygotes have reduced fitness? One solution is that this process predominantly occurs in small and isolated populations, where genetic drift can override natural selection. However, fixation is also more likely if a novel rearrangement is favored by a transmission bias, such as meiotic drive. Here, we investigate chromosomal transmission distortion in hybrids between two wood white (Leptidea sinapis) butterfly populations with extensive karyotype differences. Using data from two different crossing experiments, we uncover that there is a transmission bias favoring the ancestral chromosomal state for derived fusions, a result that shows that chromosome fusions actually can fix in populations despite being counteracted by meiotic drive. This means that meiotic drive not only can promote runaway chromosome number evolution and speciation, but also that it can be a conservative force acting against karyotypic change and the evolution of reproductive isolation. Based on our results, we suggest a mechanistic model for why chromosome fusion mutations may be opposed by meiotic drive and discuss factors contributing to karyotype evolution in Lepidoptera.

The phased Solanum okadae genome and Petota pangenome analysis of 23 other potato wild relatives and hybrids.

Potato is an important crop in the genus Solanum section Petota. Potatoes are susceptible to multiple abiotic and biotic stresses and have undergone constant improvement through breeding programs worldwide. Introgression of wild relatives from section Petota with potato is used as a strategy to enhance the diversity of potato germplasm. The current dataset contributes a phased genome assembly for diploid S. okadae, and short read sequences and de novo assemblies for the genomes of 16 additional wild diploid species in section Petota that were noted for stress resistance and were of interest to potato breeders. Genome sequence data for three additional genomes representing polyploid hybrids with cultivated potato, and an additional genome from non-tuberizing S. etuberosum, which is outside of section Petota, were also included. High quality short reads assemblies were achieved with genome sizes ranging from 575 to 795 Mbp and annotations were performed utilizing transcriptome sequence data. Genomes were compared for presence/absence of genes and phylogenetic analyses were carried out using plastome and nuclear sequences.

Heterotic growth of hybrids of Arabidopsis thaliana is enhanced by elevated atmospheric CO2.

PREMISE: With the global atmospheric CO2 concentration on the rise, developing crops that can thrive in elevated CO2 has become paramount. We investigated the potential of hybridization as a strategy for creating crops with improved growth in predicted elevated atmospheric CO2. METHODS: We grew parent accessions and their F1 hybrids of Arabidopsis thaliana in ambient and elevated atmospheric CO2 and analyzed numerous growth traits to assess their productivity and underlying mechanisms. RESULTS: The heterotic increase in total dry mass, relative growth rate and leaf net assimilation rate was significantly greater in elevated CO2 than in ambient CO2. The CO2 response of net assimilation rate was positively correlated with the CO2 response of leaf nitrogen productivity and with that of leaf traits such as leaf size and thickness, suggesting that hybridization-induced changes in leaf traits greatly affected the improved performance in elevated CO2. CONCLUSIONS: Vegetative growth of hybrids seems to be enhanced in elevated CO2 due to improved photosynthetic nitrogen-use efficiency compared with parents. The results suggest that hybrid crops should be well-suited for future conditions, but hybrid weeds may also be more competitive.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages.

BACKGROUND: Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored. METHODS: In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages. RESULTS: Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles. CONCLUSION: These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

Some plasma biomarkers of residual feed intake in beef cattle remain consistent regardless of intake level.

This study investigated whether plasma biomarkers of residual feed intake (RFI), identified under ad libitum feeding conditions in beef cattle, remained consistent during feed restriction. Sixty Charolais crossbred young bulls were divided into two groups for a crossover study. Group A was initially fed ad libitum (first test) and then restricted (second test) on the same diet, while Group B experienced the opposite sequence. Blood samples were collected from the 12 most divergent RFI animals in each group at the end of the first test and again after the second test. 12 plasma variables consistently increased, while three consistently decreased during feed restriction (FDR < 0.05). Only two metabolites, α-aminoadipic acid for Group A and 5-aminovaleric acid for Group B, were associated with RFI independent of feed intake level (FDR < 0.05), demonstrating moderate-to-high repeatability across feeding levels (intraclass correlation coefficient ≥ 0.59). Notably, both metabolites belong to the same metabolic pathway: lysine degradation. These metabolites consistently correlated with RFI, irrespective of fluctuations in feed intake, indicating a connection to individual metabolic processes influencing feed efficiency. These findings suggest that a portion of RFI phenotypic variance is inherent to an individual's metabolic efficiency beyond variations in feed intake.

Is self-incompatibility a reproductive barrier for hybridization in a sympatric species?

PREMISE: Barriers at different reproductive stages contribute to reproductive isolation. Self-incompatibility (SI) systems that prevent self-pollination could also act to control interspecific pollination and contribute to reproductive isolation, preventing hybridization. Here we evaluated whether SI contributes to reproductive isolation among four co-occurring Opuntia species that flower at similar times and may hybridize with each other. METHODS: We assessed whether Opuntia cantabrigiensis, O. robusta, O. streptacantha, and O. tomentosa, were self-compatible and formed hybrid seeds in five manipulation treatments to achieve self-pollination, intraspecific cross-pollination, open pollination (control), interspecific crosses or apomixis, then recorded flowering phenology and synchrony. RESULTS: All species flowered in the spring with a degree of synchrony, so that two pairs of species were predisposed to interspecific pollination (O. cantabrigiensis with O. robusta, O. streptacantha with O. tomentosa). All species had distinct reproductive systems: Opuntia cantabrigiensis is self-incompatible and did not produce hybrid seeds as an interspecific pollen recipient; O. robusta is a dioecious species, which formed a low proportion of hybrid seeds; O. streptacantha and O. tomentosa are self-compatible and produced hybrid seeds. CONCLUSIONS: Opuntia cantabrigiensis had a strong pollen-pistil barrier, likely due to its self-incompatibility. Opuntia robusta, the dioecious species, is an obligate outcrosser and probably partially lost its ability to prevent interspecific pollen germination. Given that the self-compatible species can set hybrid seeds, we conclude that pollen-pistil interaction and high flowering synchrony represent weak barriers; whether reproductive isolation occurs later in their life cycle (e.g., germination or seedling survival) needs to be determined.

Phenotypic variation seems not to be associated with the genetic profile in Zygopetalum (Orchidaceae): a case study of a high-elevation rocky complex.

BACKGROUND: Hybridization associated with polyploidy studies is rare in the tropics. The genus Zygopetalum (Orchidaceae) was investigated here as a case study of Neotropical plants. In the rocky highlands of the Ibitipoca State Park (ISP), southeast Brazil, individuals with intermediate colors and forms between the species Z. maculatum and Z. triste were commonly identified. METHODS AND RESULTS: Chromosomal analysis and DNA quantity showed a uniform population. Regardless of the aspects related to the color and shape of floral structures, all individuals showed 2n = 96 chromosomes and an average of 14.05 pg of DNA. Irregularities in meiosis associated with chromosome number and C value suggest the occurrence of polyploidy. The genetic distance estimated using ISSR molecular markers revealed the existence of genetic variability not related to morphological clusters. Morphometric measurements of the flower pieces revealed that Z. maculatum shows higher variation than Z. triste although lacking a defined circumscription. CONCLUSION: The observed variation can be explained by the polyploid and phenotypic plasticity resulting from the interaction of the genotypes with the heterogeneous environments observed in this habitat.

Transgressive gene expression and expression plasticity under thermal stress in a stable hybrid zone.

Interspecific hybridization can lead to myriad outcomes, including transgressive phenotypes in which the hybrids are more fit than either parent species. Such hybrids may display important traits in the context of climate change, able to respond to novel environmental conditions not previously experienced by the parent populations. While this has been evaluated in an agricultural context, the role of transgressive hybrids under changing conditions in the wild remains largely unexplored; this is especially true regarding transgressive gene expression. Using the blue mussel species complex (genus Mytilus) as a model system, we investigated the effects of hybridization on temperature induced gene expression plasticity by comparing expression profiles in parental species and their hybrids following a 2-week thermal challenge. Hybrid expression plasticity was most often like one parent or the other (50%). However, a large fraction of genes (26%) showed transgressive expression plasticity (i.e. the change in gene expression was either greater or lesser than that of both parent species), while only 2% were intermediately plastic in hybrids. Despite their close phylogenetic relationship, there was limited overlap in the differentially expressed genes responding to temperature, indicating interspecific differences in the responses to high temperature in which responses from hybrids are distinct from both parent species. We also identified differentially expressed long non-coding RNAs (lncRNAs), which we suggest may contribute to species-specific differences in thermal tolerance. Our findings provide important insight into the impact of hybridization on gene expression under warming. We propose transgressive hybrids may play an important role in population persistence under future warming conditions.

Genetic introgression from commercial European pigs to the indigenous Chinese Lijiang breed and associated changes in phenotypes.

BACKGROUND: Gene flow is crucial for enhancing economic traits of livestock. In China, breeders have used hybridization strategies for decades to improve livestock performance. Here, we performed whole-genome sequencing of a native Chinese Lijiang pig (LJP) breed. By integrating previously published data, we explored the genetic structure and introgression of genetic components from commercial European pigs (EP) into the LJP, and examined the impact of this introgression on phenotypic traits. RESULTS: Our analysis revealed significant introgression of EP breeds into the LJP and other domestic pig breeds in China. Using a haplotype-based approach, we quantified introgression levels and compared EP to LJP and other Chinese domestic pigs. The results show that EP introgression is widely prevalent in Chinese domestic pigs, although there are significant differences between breeds. We propose that LJP could potentially act as a mediator for the transmission of EP haplotypes. We also examined the correlation between EP introgression and the number of thoracic vertebrae in LJP and identified VRTN and STUM as candidate genes for this trait. CONCLUSIONS: Our study provides evidence of introgressed European haplotypes in the LJP breed and describes the potential role of EP introgression on phenotypic changes of this indigenous breed.

Extant interspecific hybridization among trematodes within the Schistosoma haematobium species complex in Nigeria.

BACKGROUND: Natural interspecific hybridization between the human parasite (Schistosoma haematobium [Sh]) and bovine parasites (Schistosoma bovis [Sb], Schistosoma curassoni [Sc]) is increasingly reported in Africa. We developed a multi-locus PCR DNA-Seq strategy that amplifies two unlinked nuclear (transITS, BF) and two linked organellar genome markers (CO1, ND5) to genotype S. haematobium eggs collected from infected people in Ile Oluji/Oke Igbo, Ondo State (an agrarian community) and Kachi, Jigawa State (a pastoral community) in Southwestern and Northern Nigeria, respectively. PRINCIPAL FINDINGS: Out of a total of 219 urine samples collected, 57 were positive for schistosomes. All patients from Jigawa state possessed an Sh mitochondrial genome and were infected with a genetic profile consistent with an Sh x Sb hybrid based on sequences obtained at CO1, ND5, transITS and BF nuclear markers. Whereas samples collected from Ondo state were more varied. Mitonuclear discordance was observed in all 17 patients, worms possessed an Sb mitochondrial genome but one of four different genetic profiles at the nuclear markers, either admixed (heterozygous between Sh x Sc or Sh x Sb) at both markers (n = 10), Sh at BF and admixed at transITS (Sh x Sc) (n = 5), admixed (Sh x Sc) at BF and homozygous Sc at transITS (n = 1) or homozygous Sh at BF and homozygous Sc at transITS (n = 1). SIGNIFICANCE: Previous work suggested that zoonotic transmission of S. bovis in pastoral communities, where humans and animals share a common water source, is a driving factor facilitating interspecific hybridization. However, our data showed that all samples were hybrids, with greater diversity identified in Southwestern Nigeria, a non-pastoral site. Further, one patient possessed an S. bovis mitochondrial genome but was homozygous for S. haematobium at BF and homozygous for S. curassoni at transITS supporting at least two separate backcrosses in its origin, suggesting that interspecific hybridization may be an ongoing process.

Mapping seasonal migration in a songbird hybrid zone -- heritability, genetic correlations, and genomic patterns linked to speciation.

Seasonal migration is a widespread behavior relevant for adaptation and speciation, yet knowledge of its genetic basis is limited. We leveraged advances in tracking and sequencing technologies to bridge this gap in a well-characterized hybrid zone between songbirds that differ in migratory behavior. Migration requires the coordinated action of many traits, including orientation, timing, and wing morphology. We used genetic mapping to show these traits are highly heritable and genetically correlated, explaining how migration has evolved so rapidly in the past and suggesting future responses to climate change may be possible. Many of these traits mapped to the same genomic regions and small structural variants indicating the same, or tightly linked, genes underlie them. Analyses integrating transcriptomic data indicate cholinergic receptors could control multiple traits. Furthermore, analyses integrating genomic differentiation further suggested genes underlying migratory traits help maintain reproductive isolation in this hybrid zone.

Italian peninsula as a hybridization zone of Ixodes inopinatus and I. ricinus and the prevalence of tick-borne pathogens in I. inopinatus, I. ricinus, and their hybrids.

BACKGROUND: Ixodes inopinatus was described from Spain on the basis of morphology and partial sequencing of 16S ribosomal DNA. However, several studies suggested that morphological differences between I. inopinatus and Ixodes ricinus are minimal and that 16S rDNA lacks the power to distinguish the two species. Furthermore, nuclear and mitochondrial markers indicated evidence of hybridization between I. inopinatus and I. ricinus. In this study, we tested our hypothesis on tick dispersal from North Africa to Southern Europe and determined the prevalence of selected tick-borne pathogens (TBPs) in I. inopinatus, I. ricinus, and their hybrids. METHODS: Ticks were collected in Italy and Algeria by flagging, identified by sequencing of partial TROSPA and COI genes, and screened for Borrelia burgdorferi s.l., B. miyamotoi, Rickettsia spp., and Anaplasma phagocytophilum by polymerase chain reaction and sequencing of specific markers. RESULTS: Out of the 380 ticks, in Italy, 92 were I. ricinus, 3 were I. inopinatus, and 136 were hybrids of the two species. All 149 ticks from Algeria were I. inopinatus. Overall, 60% of ticks were positive for at least one TBP. Borrelia burgdorferi s.l. was detected in 19.5% of ticks, and it was significantly more prevalent in Ixodes ticks from Algeria than in ticks from Italy. Prevalence of Rickettsia spotted fever group (SFG) was 51.1%, with significantly greater prevalence in ticks from Algeria than in ticks from Italy. Borrelia miyamotoi and A. phagocytophilum were detected in low prevalence (0.9% and 5.2%, respectively) and only in ticks from Italy. CONCLUSIONS: This study indicates that I. inopinatus is a dominant species in Algeria, while I. ricinus and hybrids were common in Italy. The higher prevalence of B. burgdorferi s.l. and Rickettsia SFG in I. inopinatus compared with that in I. ricinus might be due to geographical and ecological differences between these two tick species. The role of I. inopinatus in the epidemiology of TBPs needs further investigation in the Mediterranean Basin.

Integrative taxonomic analyses reveal that rapid genetic divergence drives Acropora speciation.

Reef-building corals provide the structural basis for one of Earth's most spectacular and diverse but increasingly threatened ecosystems. The reef-building coral genus Acropora may have undergone substantial speciation during the Pleistocene climate and sea-level changes. Here, we aimed to evaluate the speciation history of four morphologically similar tabular Acropora species (Acropora aff. hyacinthus, A. cf. bifurcata, A. cf. cytherea, and A. cf. subulata) using an integrative approach with morphology, genetic, and reproduction methodology. Extensive morphological analyses showed that these four species are distinct and exhibited high gamete incompatibility, preventing hybridization. Furthermore, population structure and principal component analyses with SNPs (>60,000) indicated that these species were genetically distinct, and the ABBA-BABA test did not support introgression among these species. Many of their coding and noncoding RNA sequences showed high genetic variance at loci with high Fst values along the genome. Comparison of these orthologs with those of other Acropora species suggested that many of these genes are under positive selection, which could be associated with spawning time, gamete, and morphological divergence. Our findings show that the speciation of tabular Acropora occurred without hybridization, and the divergence accompanying the rapid evolution of genes in species-rich Acropora could be associated with speciation.