Genomic evidence of prevalent hybridization throughout the evolutionary history of the fig-wasp pollination mutualism

A chromosome-level genome assembly of Ficus benjamina, a fig tree with great ecological and ornamental value

Ficus benjamina, the weeping fig, is one of the most widely distributed and cultivated figs, with important ecological functions and landscape value. However, the lack of a reference genome has hindered molecular and functional research on this well-known fig-tree. Here we present a chromosome-scale genome assembly and annotation for F. benjamina, based on a combination of Illumina short-reads, PacBio subreads, and Hi-C sequencing data. The genome consists of 13 pseudochromosomes that contain 362.73 Mb of assembled sequences, with a contig N50 length of 25.76 Mb and a complete BUSCO score of 98.10%. In total, 28,840 protein-coding genes were identified, of which 96.22% were functionally annotated. Our study provides the first chromosome-level genome of F. benjamina, providing an important resource for exploring the genetic basis of its ecological and horticultural characters.

Cospeciation is not the dominant driver of plant-pollinator codiversification in specialized pollination systems

This study systematically rejects the long-standing notion of cospeciation as the dominant driver of codiversification between flowering plants and their specialist pollinators. Through cophylogenetic analysis of six classical specialized pollination systems, the research finds that cospeciation events are consistently outnumbered by non-cospeciation events, such as host-switch, duplication, and association losses. The findings support a more dynamic and diffuse codiversification paradigm, highlighting the importance of considering a broader range of evolutionary events in understanding plant-pollinator codiversification. This new understanding is robust across diverse pollination systems and has significant implications for conservation strategies in the face of environmental change.

Peering through the hedge: Multiple datasets yield insights into the phylogenetic relationships and incongruences in the tribe Lilieae (Liliaceae)

The increasing use of genome-scale data has significantly facilitated phylogenetic analyses, contributing to the dissection of the underlying evolutionary mechanisms that shape phylogenetic incongruences, such as incomplete lineage sorting (ILS) and hybridization. Lilieae, a prominent member of the Liliaceae family, comprises four genera and approximately 260 species, representing 43% of all species within Liliaceae. They possess high ornamental, medicinal and edible values. Yet, no study has explored the validity of various genome-scale data in phylogenetic analyses within this tribe, nor have potential evolutionary mechanisms underlying its phylogenetic incongruences been investigated. Here, transcriptome, Angiosperms353, plastid and mitochondrial data, were collected from 50 to 93 samples of Lilieae, covering all four recognized genera. Multiple datasets were created and used for phylogenetic analyses based on concatenated and coalescent-based methods. Evolutionary rates of different datasets were calculated, and divergence times were estimated. Various approaches, including coalescence simulation, Quartet Sampling (QS), calculation of concordance factors (gCF and sCF), as well as MSCquartets and reticulate network inference, were carried out to infer the phylogenetic discordances and analyze their underlying mechanisms using a reduced 33-taxon dataset. Despite extensive phylogenetic discordances among gene trees, robust phylogenies were inferred from nuclear and plastid data compared to mitochondrial data, with lower synonymous substitution detected in mitochondrial genes than in nuclear and plastid genes. Significant ILS was detected across the phylogeny of Lilieae, with clear evidence of reticulate evolution identified. Divergence time estimation indicated that most of lineages in Lilieae diverged during a narrow time frame (ranging from 5.0 Ma to 10.0 Ma), consistent with the notion of rapid radiation evolution. Our results suggest that integrating transcriptomic and plastid data can serve as cost-effective and efficient tools for phylogenetic inference and evolutionary analysis within Lilieae, and Angiosperms353 data is also a favorable choice. Mitochondrial data are more suitable for phylogenetic analyses at higher taxonomic levels due to their stronger conservation and lower synonymous substitution rates. Significant phylogenetic incongruences detected in Lilieae were caused by both incomplete lineage sorting (ILS) and reticulate evolution, with hybridization and "ghost introgression" likely prevalent in the evolution of Lilieae species. Our findings provide new insights into the phylogeny of Lilieae, enhancing our understanding of the evolution of species in this tribe.

Intraspecific Differentiation of Styrax japonicus (Styracaceae) as Revealed by Comparative Chloroplast and Evolutionary Analyses

Styrax japonicus is a medicinal and ornamental shrub belonging to the Styracaceae family. To explore the diversity and characteristics of the chloroplast genome of S. japonicus, we conducted sequencing and comparison of the chloroplast genomes of four naturally distributed S. japonicus. The results demonstrated that the four chloroplast genomes (157,914-157,962 bp) exhibited a typical quadripartite structure consisting of a large single copy (LSC) region, a small single copy (SSC) region, and a pair of reverse repeats (IRa and IRb), and the structure was highly conserved. DNA polymorphism analysis revealed that three coding genes (infA, psbK, and rpl33) and five intergene regions (petA-psbJ, trnC-petN, trnD-trnY, trnE-trnT, and trnY-trnE) were identified as mutation hotspots. These genetic fragments have the potential to be utilized as DNA barcodes for future identification purposes. When comparing the boundary genes, a small contraction was observed in the IR region of four S. japonicus. Selection pressure analysis indicated positive selection for ycf1 and ndhD. These findings collectively suggest the adaptive evolution of S. japonicus. The phylogenetic structure revealed conflicting relationships among several S. japonicus, indicating divergent evolutionary paths within this species. Our study concludes by uncovering the genetic traits of the chloroplast genome in the differentiation of S. japonicus variety, offering fresh perspectives on the evolutionary lineage of this species.

Hybridization, polyploidization, and morphological convergence make dozens of taxa into one chaotic genetic pool: a phylogenomic case of the Ficus erecta species complex (Moraceae)

The Ficus erecta complex, characterized by its morphological diversity and frequent interspecific overlap, shares pollinating fig wasps among several species. This attribute, coupled with its intricate phylogenetic relationships, establishes it as an exemplary model for studying speciation and evolutionary patterns. Extensive researches involving RADseq (Restriction-site associated DNA sequencing), complete chloroplast genome data, and flow cytometry methods were conducted, focusing on phylogenomic analysis, genetic structure, and ploidy detection within the complex. Significantly, the findings exposed a pronounced nuclear-cytoplasmic conflict. This evidence, together with genetic structure analysis, confirmed that hybridization within the complex is a frequent occurrence. The ploidy detection revealed widespread polyploidy, with certain species exhibiting multiple ploidy levels, including 2×, 3×, and 4×. Of particular note, only five species (F. abelii, F. erecta, F. formosana, F. tannoensis and F. vaccinioides) in the complex were proved to be monophyletic. Species such as F. gasparriniana, F. pandurata, and F. stenophylla were found to encompass multiple phylogenetically distinct lineages. This discovery, along with morphological comparisons, suggests a significant underestimation of species diversity within the complex. This study also identified F. tannoensis as an allopolyploid species originating from F. vaccinioide and F. erecta. Considering the integration of morphological, molecular systematics, and cytological evidences, it is proposed that the scope of the F. erecta complex should be expanded to the entire subsect. Frutescentiae. This would redefine the complex as a continuously evolving group comprising at least 33 taxa, characterized by blurred species boundaries, frequent hybridization and polyploidization, and ambiguous genetic differentiation.

Characterisation of the fig-fig wasp holobiont

Plants and animals have long been considered distinct kingdoms, yet here a 'plant-animal' is described. An extraordinary symbiosis in which neither organism can reproduce without the other, the fig tree (Ficus) provides the habitat for its exclusive pollinator: the fig wasp (Agaonidae). Characterising the 'fig-fig wasp holobiont' acknowledges, for the first time, 'plant-animal symbiogenesis'.

A telomere-to-telomere reference genome of ficus (Ficus hispida) provides new insights into sex determination

A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers' defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.

Echoes of ancient introgression punctuate stable genomic lineages in the evolution of figs

Studies investigating the evolution of flowering plants have long focused on isolating mechanisms such as pollinator specificity. Some recent studies have proposed a role for introgressive hybridization between species, recognizing that isolating processes such as pollinator specialization may not be complete barriers to hybridization. Occasional hybridization may therefore lead to distinct yet reproductively connected lineages. We investigate the balance between introgression and reproductive isolation in a diverse clade using a densely sampled phylogenomic study of fig trees (Ficus, Moraceae). Codiversification with specialized pollinating wasps (Agaonidae) is recognized as a major engine of fig diversity, leading to about 850 species. Nevertheless, some studies have focused on the importance of hybridization in Ficus, highlighting the consequences of pollinator sharing. Here, we employ dense taxon sampling (520 species) throughout Moraceae and 1,751 loci to investigate phylogenetic relationships and the prevalence of introgression among species throughout the history of Ficus. We present a well-resolved phylogenomic backbone for Ficus, providing a solid foundation for an updated classification. Our results paint a picture of phylogenetically stable evolution within lineages punctuated by occasional local introgression events likely mediated by local pollinator sharing, illustrated by clear cases of cytoplasmic introgression that have been nearly drowned out of the nuclear genome through subsequent lineage fidelity. The phylogenetic history of figs thus highlights that while hybridization is an important process in plant evolution, the mere ability of species to hybridize locally does not necessarily translate into ongoing introgression between distant lineages, particularly in the presence of obligate plant-pollinator relationships.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Comparative and phylogenetic analyses of the chloroplast genome reveal the taxonomy of the Morus genus

Mulberry (genus Morus) is an economically important woody plant with an altered ploidy level. The variable number of Morus species recognized by different studies indicates that the genus is in need of revision. In this study, the chloroplast (CP) genomes of 123 Morus varieties were de novo assembled and systematically analyzed. The 123 varieties represented six Morus species, namely, Morus alba, Morus nigra, Morus notabilis, Morus rubra, Morus celtidifolia, and Morus serrata. The Morus CP genome was found to be 158,969~159,548 bp in size with 125 genes, including 81 protein coding, 36 tRNA, and 8 rRNA genes. The 87 out of 123 mulberry accessions were assigned to 14 diverse groups with identical CP genome, which indicated that they are maternally inherited and share 14 common ancestors. Then 50 diverse CP genomes occurred in 123 mulberry accessions for further study. The CP genomes of the Morus genus with a quadripartite structure have two inverted repeat (IR) regions (25,654~25,702 bp) dividing the circular genome into a large single-copy (LSC) region (87,873~88,243 bp) and small single-copy (SSC) region (19,740~19,994 bp). Analysis of the phylogenetic tree constructed using the complete CP genome sequences of Morus revealed a monophyletic genus and that M. alba consisted of two clades, M. alba var. alba and M. alba var. multicaulis. The Japanese cultivated germplasms were derived from M. alba var. multicaulis. We propose that the Morus genus be classified into six species, M. nigra, M. notabilis, M. serrata, M. celtidifolia, M. rubra, and M. alba with two subspecies, M. alba var. alba and M. alba var. multicaulis. Our findings provide a valuable resource for the classification, domestication, and breeding improvement of mulberry.

Genome sequencing and comparative analysis of Ficus benghalensis and Ficus religiosa species reveal evolutionary mechanisms of longevity

Ficus benghalensis and Ficus religiosa are large woody trees well known for their long lifespan, ecological and traditional significance, and medicinal properties. To understand the genomic and evolutionary aspects of these characteristics, the whole genomes of these Ficus species were sequenced using 10x Genomics linked reads and Oxford Nanopore long reads. The draft genomes of F. benghalensis and F. religiosa comprised of 392.89 Mbp and 332.97 Mbp, respectively. We established the genome-wide phylogenetic positions of the two Ficus species with respect to 50 other Angiosperm species. Comparative evolutionary analyses with other phylogenetically closer Eudicot species revealed adaptive evolution in genes involved in key cellular mechanisms associated with prolonged survival including phytohormones signaling, senescence, disease resistance, and abiotic stress tolerance, which provide genomic insights into the mechanisms conferring longevity and suggest that longevity is a multifaceted phenomenon. This study also provides clues on the existence of CAM pathway in these Ficus species.

Comparative chloroplast genome analysis of Ficus (Moraceae): Insight into adaptive evolution and mutational hotspot regions

As the largest genus in Moraceae, Ficus is widely distributed across tropical and subtropical regions and exhibits a high degree of adaptability to different environments. At present, however, the phylogenetic relationships of this genus are not well resolved, and chloroplast evolution in Ficus remains poorly understood. Here, we sequenced, assembled, and annotated the chloroplast genomes of 10 species of Ficus, downloaded and assembled 13 additional species based on next-generation sequencing data, and compared them to 46 previously published chloroplast genomes. We found a highly conserved genomic structure across the genus, with plastid genome sizes ranging from 159,929 bp (Ficus langkokensis) to 160,657 bp (Ficus religiosa). Most chloroplasts encoded 113 unique genes, including a set of 78 protein-coding genes, 30 transfer RNA (tRNA) genes, four ribosomal RNA (rRNA) genes, and one pseudogene (infA). The number of simple sequence repeats (SSRs) ranged from 67 (Ficus sagittata) to 89 (Ficus microdictya) and generally increased linearly with plastid size. Among the plastomes, comparative analysis revealed eight intergenic spacers that were hotspot regions for divergence. Additionally, the clpP, rbcL, and ccsA genes showed evidence of positive selection. Phylogenetic analysis indicated that none of the six traditionally recognized subgenera of Ficus were monophyletic. Divergence time analysis based on the complete chloroplast genome sequences showed that Ficus species diverged rapidly during the early to middle Miocene. This research provides basic resources for further evolutionary studies of Ficus.

A chromosome-level genome assembly of the pollinating fig wasp Valisia javana

Fig wasp has always been thought the species-specific pollinator for their host fig (Moraceae, Ficus) and constitute a model system with its host to study co-evolution and co-speciation. The availability of a high-quality genome will help to further reveal the mechanisms underlying these characteristics. Here, we present a high-quality chromosome-level genome for Valisa javana developed by a combination of PacBio long-read and Illumina short-read. The assembled genome size is 296.34 Mb from 13 contigs with a contig N50 length of 26.76 kb. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid fig wasps living in syconium of its highly specific host. Protein-coding genes associated with chemosensory, detoxification and venom genes were identified. Several differentially expressed genes in transcriptome data of V. javana between odor-stimulated samples and the controls have been identified in some olfactory signal transduction pathways, e.g. olfactory transduction, cAMP, cGMP-PKG, Calcim, Ras and Rap1. This study provides a valuable genomic resource for a fig wasp, and sheds insight into further revealing the mechanisms underlying their adaptive traits to their hosts in different places and co-speciation with their host.

Phylogenomic approaches untangle early divergences and complex diversifications of the olive plant family

BACKGROUND

Deep-branching phylogenetic relationships are often difficult to resolve because phylogenetic signals are obscured by the long history and complexity of evolutionary processes, such as ancient introgression/hybridization, polyploidization, and incomplete lineage sorting (ILS). Phylogenomics has been effective in providing information for resolving both deep- and shallow-scale relationships across all branches of the tree of life. The olive family (Oleaceae) is composed of 25 genera classified into five tribes with tribe Oleeae consisting of four subtribes. Previous phylogenetic analyses showed that ILS and/or hybridization led to phylogenetic incongruence in the family. It was essential to distinguish phylogenetic signal conflicts, and explore mechanisms for the uncertainties concerning relationships of the olive family, especially at the deep-branching nodes.

RESULTS

We used the whole plastid genome and nuclear single nucleotide polymorphism (SNP) data to infer the phylogenetic relationships and to assess the variation and rates among the main clades of the olive family. We also used 2608 and 1865 orthologous nuclear genes to infer the deep-branching relationships among tribes of Oleaceae and subtribes of tribe Oleeae, respectively. Concatenated and coalescence trees based on the plastid genome, nuclear SNPs and multiple nuclear genes suggest events of ILS and/or ancient introgression during the diversification of Oleaceae. Additionally, there was extreme heterogeneity in the substitution rates across the tribes. Furthermore, our results supported that introgression/hybridization, rather than ILS, is the main factor for phylogenetic discordance among the five tribes of Oleaceae. The tribe Oleeae is supported to have originated via ancient hybridization and polyploidy, and its most likely parentages are the ancestral lineage of Jasmineae or its sister group, which is a "ghost lineage," and Forsythieae. However, ILS and ancient introgression are mainly responsible for the phylogenetic discordance among the four subtribes of tribe Oleeae.

CONCLUSIONS

This study showcases that using multiple sequence datasets (plastid genomes, nuclear SNPs and thousands of nuclear genes) and diverse phylogenomic methods such as data partition, heterogeneous models, quantifying introgression via branch lengths (QuIBL) analysis, and species network analysis can facilitate untangling long and complex evolutionary processes of ancient introgression, paleopolyploidization, and ILS.

Pollinator sharing, copollination, and speciation by host shifting among six closely related dioecious fig species

The obligate pollination mutualism between figs (Ficus, Moraceae) and pollinator wasps (Agaonidae, Hymenoptera) is a classic example of cospeciation. However, examples of phylogenetic incongruencies between figs and their pollinators suggest that pollinators may speciate by host shifting. To investigate the mechanism of speciation by host shifting, we examined the phylogenetic relationships and population genetic structures of six closely related fig species and their pollinators from southern China and Taiwan-Ryukyu islands using various molecular markers. The results revealed 1) an extraordinary case of pollinator sharing, in which five distinct fig species share a single pollinator species in southern China; 2) two types of copollination, namely, sympatric copollination by pollinator duplication or pollinator migration, and allopatric copollination by host migration and new pollinator acquisition; 3) fig species from southern China have colonized Taiwan repeatedly and one of these events has been followed by host shifting, reestablishment of host specificity, and pollinator speciation, in order. Based on our results, we propose a model for pollinator speciation by host shifting in which the reestablishment of host-specificity plays a central role in the speciation process. These findings provide important insights into understanding the mechanisms underlying pollinator speciation and host specificity in obligate pollination mutualism.

Phylogenetic and population genetic analyses reveal mechanisms of pollinator sharing, copollination and speciation by host-shift in fig-wasp mutualism.

Bursaphelenchus suri n. sp.: A second Bursaphelenchus syconial parasite of figs supports adaptive radiation among section Sycomorus figs

The nematode genus Bursaphelenchus is a highly divergent group. This genus mainly consists of mycophagous entomophilic species, but some species have specialized as obligate or facultative plant parasites, facultative insect parasites, or exhibit feeding dimorphism (phenotypic plasticity) leading to mycophagous and predatory forms. In the present study, a new Bursaphelenchus species, B. suri n. sp. was isolated from fresh syconia (figs) of Ficus sur and is described and illustrated based on its typological characters and molecular phylogenetic status. The new species is characterized by its highly derived feeding structures found in obligate plant parasites, lip possessing a labial disc and a long and thick stylet with a long conus and extremely well-developed basal swellings. In addition, slender body of both sexes is characteristic of the species. The new species is phylogenetically and typologically closely related to B. sycophilus, i.e., these two species share the characteristic feeding structures and form a well-supported clade within the B. fungivorus group in the genus. Biologically, these two species are both isolated from fresh figs of the section Sycomorus. However, the new species differs from B. sycophilus by the length of the female post-uterine sac and the shape of the male spicule, i.e., the new species has a long post-uterine sac and spicule condylus without dorsal recurvature. Thus, the new species is the second obligate fig parasite of the genus, and the evolutionary relationship between the B. suri n. sp. and B. sycophilus clade and section Sycomorus figs is hypothesized as an example of adaptive radiation with more species to be discovered.

The Syngameon Enigma

Despite their evolutionary relevance, multispecies networks or syngameons are rarely reported in the literature. Discovering how syngameons form and how they are maintained can give insight into processes such as adaptive radiations, island colonizations, and the creation of new hybrid lineages. Understanding these complex hybridization networks is even more pressing with anthropogenic climate change, as syngameons may have unique synergistic properties that will allow participating species to persist. The formation of a syngameon is not insurmountable, as several ways for a syngameon to form have been proposed, depending mostly on the magnitude and frequency of gene flow events, as well as the relatedness of its participants. Episodic hybridization with small amounts of introgression may keep syngameons stable and protect their participants from any detrimental effects of gene flow. As genomic sequencing becomes cheaper and more species are included in studies, the number of known syngameons is expected to increase. Syngameons must be considered in conservation efforts as the extinction of one participating species may have detrimental effects on the survival of all other species in the network.

Chloranthus genome provides insights into the early diversification of angiosperms

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

Intraspecies Genomic Divergence of a Fig Wasp Species Is Due to Geographical Barrier and Adaptation

Understanding how intraspecies divergence results in speciation has great importance for our knowledge of evolutionary biology. Here we applied population genomics approaches to a fig wasp species (Valisia javana complex sp 1) to reveal its intraspecies differentiation and the underlying evolutionary dynamics. With re-sequencing data, we prove the Hainan Island population (DA) of sp1 genetically differ from the continental ones, then reveal the differed divergence pattern. DA has reduced SNP diversity but a higher proportion of population-specific structural variations (SVs), implying a restricted gene exchange. Based on SNPs, 32 differentiated islands containing 204 genes were detected, along with 1,532 population-specific SVs of DA overlapping 4,141 genes. The gene ontology (GO) enrichment analysis performed on differentiated islands linked to three significant GO terms on a basic metabolism process, with most of the genes failing to enrich. In contrast, population-specific SVs contributed more to the adaptation than the SNPs by linking to 59 terms that are crucial for wasp speciation, such as host reorganization and development regulation. In addition, the generalized dissimilarity modeling confirms the importance of environment difference on the genetic divergence within sp1. Hence, we assume the genetic divergence between DA and the continent due to not only the strait as a geographic barrier, but also adaptation. We reconstruct the demographic history within sp1. DA shares a similar population history with the nearby continental population, suggesting an incomplete divergence. Summarily, our results reveal how geographic barriers and adaptation both influence the genetic divergence at population-level, thereby increasing our knowledge on the potential speciation of non-model organisms.

Adaptation of Fig Wasps (Agaodinae) to Their Host Revealed by Large-Scale Transcriptomic Data

Simple Summary

Research on fig wasps has made a considerable contribution to the understanding of insect–plant interactions. However, the molecular mechanisms underlying fig wasp host specificity are poorly understood. This study reports on a relatively large-scale transcriptomic dataset of 25 fig wasp species. We outline potential genetic mechanisms underlying the specific host adaptation by investigating changes in a gene family, in evolutionary rates, and in genes under positive selection. The transcriptome datasets reported here (1) provide new insights into the evolutionary diversification and host specificity of fig wasps and (2) contribute to a growing dataset on fig wasp genomics.

Abstract

Figs and fig wasps are highly species-specific and comprise a model system for studying co-evolution and co-speciation. The evolutionary relationships and molecular adaptations of fig wasps to their fig hosts are poorly understood, and this is in part due to limited sequence data. Here, we present large-scale transcriptomic datasets of 25 fig wasp species with the aim of uncovering the genetic basis for host specificity. Our phylogenetic results support the monophyly of all genera associated with dioecious figs, and two genera associated with monoecious figs, Eupristina and Platyscapa, were revealed to be close relatives. We identified gene loss and gain, potentially rapidly evolving genes, and genes under positive selection. Potentially functional changes were documented and we hypothesize as to how these may determine host specificity. Overall, our study provides new insights into the evolutionary diversification of fig wasps and contributes to our understanding of adaptation in this group.

Is speciation an unrelenting march to reproductive isolation?

Speciation is often portrayed as an "incomplete" or "incipient" process if two groups of organisms, technically distinguishable either by morphology or genetics, can exchange genes. The ultimate outcome of diversification, given this perspective, is complete reproductive isolation. But an increasing amount of evidence suggests that speciation is rarely complete and inter-fertility between different taxonomically accepted species is consistently maintained. In this issue of Molecular Ecology, Linan et al. (2021) provide results that bridge evolutionary processes from populations to phylogenies that indicate suites of closely related tree species in the Mascarene Islands actively exchange genes, evolving as a nested set of syngameons with a hierarchical pattern of interfertility. The deep insight into diversification provided by this study is particularly powerful because of the genomic scale of the data and the complete taxonomic sampling of an island clade evolving in situ. The prevalence of syngameon dynamics in a broad range of organisms indicates that we should adopt a fluid and comprehensive approach to defining evolutionary units for conservation and research. We should move beyond focusing on single endangered species in evolutionary and ecological isolation from other species but consider the entire network of potentially interfertile species and the potential for future adaptation and innovation, particularly in a human dominated world.

The evolution of parasitism from mutualism in wasps pollinating the fig, Ficus microcarpa, in Yunnan Province, China

Theory identifies factors that can undermine the evolutionary stability of mutualisms. However, theory's relevance to mutualism stability in nature is controversial. Detailed comparative studies of parasitic species that are embedded within otherwise mutualistic taxa (e.g., fig pollinator wasps) can identify factors that potentially promote or undermine mutualism stability. We describe results from behavioral, morphological, phylogenetic, and experimental studies of two functionally distinct, but closely related, Eupristina wasp species associated with the monoecious host fig, Ficus microcarpa, in Yunnan Province, China. One (Eupristina verticillata) is a competent pollinator exhibiting morphologies and behaviors consistent with observed seed production. The other (Eupristina sp.) lacks these traits, and dramatically reduces both female and male reproductive success of its host. Furthermore, observations and experiments indicate that individuals of this parasitic species exhibit greater relative fitness than the pollinators, in both indirect competition (individual wasps in separate fig inflorescences) and direct competition (wasps of both species within the same fig). Moreover, phylogenetic analyses suggest that these two Eupristina species are sister taxa. By the strictest definition, the nonpollinating species represents a "cheater" that has descended from a beneficial pollinating mutualist. In sharp contrast to all 15 existing studies of actively pollinated figs and their wasps, the local F. microcarpa exhibit no evidence for host sanctions that effectively reduce the relative fitness of wasps that do not pollinate. We suggest that the lack of sanctions in the local hosts promotes the loss of specialized morphologies and behaviors crucial for pollination and, thereby, the evolution of cheating.

Fitness costs for fig wasps that fail to pollinate their host Ficus perforata

Mutualisms are of fundamental ecological importance, but risk breaking down if one partner stops paying the costs yet still takes the benefits of the interaction. To prevent such cheating, many mutualisms have mechanisms that lower the fitness of uncooperative symbionts, often termed host sanctions. In mutualisms where the interacting partners are species-specific, we would expect to see coevolution of the levels of host sanctions and partner cooperation across species-pairs. In the mutualism between fig trees and their species-specific pollinating fig wasps, host sanctions vary greatly in strength, and wasp cooperation levels vary accordingly. Here I show experimentally that in Panamanian Ficus perforata (section Urostigma, Americana) there are fitness costs for wasps that do not pollinate. These fitness costs are caused by a combination of abortions of unpollinated figs and reduced proportion of wasp larvae that successfully develop to adults. The relative fitness of wasps that do not pollinate compared to wasps that pollinate is 0.59, leading to the intermediate sanction strength 0.41. Next, by screening pollinators of F. perforata I found that 1.9% of wasp individuals in natural populations failed to carry pollen. Across five actively pollinated Neotropical fig species and their pollinators, fig species with stronger host sanctions had fewer uncooperative wasps, as would be expected if sanctions promote cooperation.

Molecular mechanisms of mutualistic and antagonistic interactions in a plant-pollinator association

Many insects metamorphose from antagonistic larvae into mutualistic adult pollinators, with reciprocal adaptation leading to specialized insect-plant associations. It remains unknown how such interactions are established at molecular level. Here we assemble high-quality genomes of a fig species, Ficus pumila var. pumila, and its specific pollinating wasp, Wiebesia pumilae. We combine multi-omics with validation experiments to reveal molecular mechanisms underlying this specialized interaction. In the plant, we identify the specific compound attracting pollinators and validate the function of several key genes regulating its biosynthesis. In the pollinator, we find a highly reduced number of odorant-binding protein genes and an odorant-binding protein mainly binding the attractant. During antagonistic interaction, we find similar chemical profiles and turnovers throughout the development of galled ovules and seeds, and a significant contraction of detoxification-related gene families in the pollinator. Our study identifies some key genes bridging coevolved mutualists, establishing expectations for more diffuse insect-pollinator systems.