First Evidence of Thalassochory in the Ficus Genus: Seed Dispersal Using the Kuroshio Oceanic Current

AIM

Plants distributed between southern Taiwan and the north of the Philippines are spread among numerous small islands in an area crossed by the powerful Kuroshio current. Oceanic currents can be effective seed-dispersal agents for coastal plant species. Moreover, the Luzon Strait is an area prone to tropical cyclones. The aim of this study is to look at the dispersal capability of an endangered coastal plant species, the Mearns fig (Ficus pedunculosa var. mearnsii), using both experimental and population genetics methods.

LOCATION

Southern Taiwan, the Philippines, and the islands between Luzon and Taiwan Island.

METHODS

This study combined two types of analysis, i.e., buoyancy experiments on syconia and double digest restriction-associated DNA sequencing (ddRAD), to analyze the population genetics of the Mearns fig.

RESULTS

We first discovered that mature Mearns fig syconia could float in seawater. They have a mean float duration of 10 days to a maximum of 21 days. Germination rates varied significantly between Mearns fig seeds that had undergone different durations of flotation treatment. Population genetic analysis shows a high degree of inbreeding among various Mearns fig populations. Moreover, no isolation by distance was found between the populations and individuals.

MAIN CONCLUSIONS

From our analysis of the genetic structure of the Mearns fig populations, we can clearly highlight the effect of the Kuroshio oceanic current on the seed dispersal of this fig tree. Comprehensive analysis has shown that Mearns fig seeds are still viable before the mature syconium sinks into the seawater, and so they could use the Kuroshio Current to float to the current population locations in Taiwan.

Genomic diversity and structure of a Neotropical microendemic fig tree

Genetic diversity is a key component of evolution, and unraveling factors that promote genetic differentiation in space and time is a central question in evolutionary biology. One of the most diverse and ecologically important tree genera in tropical forests worldwide is Ficus (Moraceae). It has been suggested that, given the great dispersal capacity of pollinating fig wasps (Chalcidoidea; Agaonidae), the spatial genetic structure, particularly in monoecious fig species, should be weak. However, no studies have addressed the factors that determine the genetic structure of Ficus species in regions of high geological, geographic, and climatic complexity, such as the Mexican Transition Zone. Using nuclear single nucleotide polymorphisms (5311 SNPs) derived from low-coverage whole genomes and 17 populations, we analyzed the population genomics of Ficus pringlei to characterize neutral and adaptive genetic variation and structure and its association with geographic barriers such as the Trans-Mexican Volcanic Belt, environmental heterogeneity, and wind connectivity. From genomic data of 71 individuals, high genetic diversity, and the identification of three genomic lineages were recorded (North, South, and Churumuco). The results suggest that genetic variation is primarily determined by climatic heterogeneity. Ficus pringlei populations from the north and south of the Trans-Mexican Volcanic Belt also exhibited minimal genetic differentiation (F ST = 0.021), indicating that this mountain range may not act as an insurmountable barrier to gene flow. Wind connectivity is also highlighted in structuring putative adaptive genetic variation, underscoring the intricate complexity of the various factors influencing genetic variation in the species. This study provides information on the possible mechanisms underlying the genetic variation of endemic species of the tropical dry forest of Western Mexico, such as F. pringlei.

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.

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.

Development of 19 universal microsatellite loci for three closely related Ficus species (Moraceae) by high-throughput sequencing

The genus Ficus is striking for its species diversity, ecological significance, and its often species-specific relationship with coevolved pollinating fig wasps, which has long fascinated biologists. The three closely related and generally co-distributed dioecious species Ficus hispida, F. heterostyla and F. squamosa provide an ideal system for the study of speciation, hybridization (caused by pollinator sharing) and comparative phylogeography to infer historical biogeography. We aimed to develop microsatellite markers for these allied species to facilitate the outlined study investigations. A DNA library was constructed from one F. heterostyla sample, and 19 microsatellite loci were developed based on high-throughput sequencing. These markers showed relatively high polymorphism in all three fig species. The mean number of alleles per locus was 3.594-5.286, and the mean observed and expected heterozygosity were 0.469-0.546 and 0.467-0.528, respectively. Principal coordinate, STRUCTURE and AMOVA analyses revealed different degrees of genetic differentiation within species, and, despite some observed genetic admixture, indicated the presence of clear boundaries between different species. In summary, we successfully developed universal microsatellite markers for three closely related Ficus species. These markers will be of great value for investigating patterns of biodiversity among the species in this model system for coevolutionary studies.

Janzen-Connell effects shape gene flow patterns and realized fitness in the tropical dioecious tree Spondias purpurea (ANACARDIACEAE)

Pollination and seed dispersal patterns determine gene flow within plant populations. In tropical forests, a high proportion of trees are dioecious, insect pollinated and dispersed by vertebrates. Dispersal vectors and density dependent factors may modulate realized gene flow and influence the magnitude of Fine Scale Genetic Structure (FSGS), affecting individual fitness. Spondias purpurea is a vertebrate-dispersed, insect-pollinated dioecious tropical tree. We assessed the influence of sex ratio, effective and realized gene flow on genetic diversity, FSGS and individual fitness within a 30 ha plot in the tropical dry forest reserve of Chamela-Cuixmala, Mexico. All individuals within the plot were tagged, geo-referenced and sampled for genetic analysis. We measured dbh and monitored sex expression during two reproductive seasons for all individuals. We collected seeds directly from maternal trees for effective pollen dispersal analysis, and analyzed established seedlings to assess realized pollen and seed dispersal. Nine microsatellite loci were used to describe genetic diversity parameters, FSGS and gene flow patterns among different size classes. A total of 354 individuals were located and classified into three size classes based on their dbh (<10, 10–20, and >20 cm). Population sex ratios were male biased and diametric size distributions differed among sexes, these differences may be the result of precocious male reproduction at early stages. Autocorrelation analyses indicate low FSGS (Fj <0.07) across all size classes. Long realized pollen and seed dispersal and differences among effective and realized gene flow were detected. In our study site low FSGS is associated with high gene flow levels. Effective and realized gene flow indicate a population recruitment curve indicating Janzen-Connell effects and suggesting fitness advantages for long-distance pollen and seed dispersal events.

Assessment of genetic diversity, population structure and sex identification in dioecious crop, Trichosanthes dioica employing ISSR, SCoT and SRAP markers

Twenty inter-simple sequence repeat (ISSR) and twenty two start codon targeted (SCoT) primers were employed to analyze genetic diversity and population structure among 52 Trichosanthes dioica Roxb. accessions collected from nine different eco-geographical regions of India. ISSR markers proved to be more informative in genetic diversity assessment and produced higher mean number of polymorphic bands (15.25 with 95.96% polymorphism) and polymorphic information content (PIC) value (0.47) compared to SCoT markers (12.55 polymorphic bands with 92.20% polymorphism and PIC: 0.45). Total genetic diversity (Ht) and genetic diversity within populations (Hs) in T. dioica accessions was found to be very high (0.45 and 0.43, respectively). AMOVA analysis also revealed higher genetic variation within populations (81%) than among them (19%). Among different T. dioica populations, very low genetic differentiation (Gst: 0.05) and high gene flow (Nm: 9.32) were observed. T. dioica populations of Bihar state were found to be highly diverse and Kolkata and Cuttack populations were least diverse. T. dioica male plants were more variable than females. UPGMA, Neighbor-Joining and population structure analyses divided T. dioica populations into three main clusters. First cluster comprised of Meerut population, second cluster included of Cuttack and Kolkata populations and populations of Bihar, Delhi and Kanpur occurred in third cluster. Genetic diversity was found to be strongly positively correlated with the latitude and strongly negatively correlated with annual mean rainfall of different T. dioica cultivated regions. For sex identification, one SRAP primer combination, 'Em-6/Me-4' amplified two molecular markers of around 230 and 290 bp specific to male T. dioica plants of Bihar, Kanpur, North Delhi and Meerut populations and were completely absent from female plants.

Resource dispersion influences dispersal evolution of highly insulated insect communities

Communities in which species are obligately associated with a single host are ideal to test adaptive responses of community traits to host-imposed selection because such communities are often highly insulated. Fig species provide oviposition resources to co-evolved fig-wasp communities. Dispersing fig-wasp communities move from one host plant to another for oviposition. We compared the spatial dispersion of two fig species and the dispersal capacities of their multitrophic wasp communities. Dispersal capacities were assessed by measuring vital dispersal correlates, namely tethered flight durations, somatic lipid contents and resting metabolic rates. We suggest that dispersal-trait distributions of congeneric wasp species across the communities are an adaptive response to host plant dispersion. Larger dispersal capacities of the entire multitrophic community are related to more widely dispersed resources. Our results provide evidence and a novel perspective for understanding the potential role of adaptation in whole-community dispersal-trait distributions.

Effects of zoochory on the spatial genetic structure of plant populations

Spatial genetic structure (SGS) of plants results from the nonrandom distribution of related individuals. SGS provides information on gene flow and spatial patterns of genetic diversity within populations. Seed dispersal creates the spatial template for plant distribution. Thus, in zoochorous plants, dispersal mode and disperser behaviour might have a strong impact on SGS. However, many studies only report the taxonomic group of seed dispersers, without further details. The recent increase in studies on SGS provides the opportunity to review findings and test for the influence of dispersal mode, taxonomic affiliation of dispersers and their behaviour. We compared the proportions of studies with SGS among groups and tested for differences in strength of SGS using Sp statistics. The presence of SGS differed among taxonomic groups, with reduced presence in plants dispersed by birds. Strength of SGS was instead significantly influenced by the behaviour of seed dispersal vectors, with higher SGS in plant species dispersed by animals with behavioural traits that result in short seed dispersal distances. We observed high variance in the strength of SGS in plants dispersed by animals that actively or passively accumulate seeds. Additionally, we found SGS was also affected by pollination and marker type used. Our study highlights the importance of vector behaviour on SGS even in the presence of variance created by other factors. Thus, more detailed information on the behaviour of seed dispersers would contribute to better understand which factors shape the spatial scale of gene flow in animal-dispersed plant species.

Spatial Scales of Genetic Structure in Free-Standing and Strangler Figs (Ficus, Moraceae) Inhabiting Neotropical Forests

Wind-borne pollinating wasps (Agaonidae) can transport fig (Ficus sp., Moraceae) pollen over enormous distances (> 100 km). Because of their extensive breeding areas, Neotropical figs are expected to exhibit weak patterns of genetic structure at local and regional scales. We evaluated genetic structure at the regional to continental scale (Panama, Costa Rica, and Peru) for the free-standing fig species Ficus insipida. Genetic differentiation was detected only at distances > 300 km (Jost´s D_est = 0.68 ± 0.07 & F_ST = 0.30 ± 0.03 between Mesoamerican and Amazonian sites) and evidence for phylogeographic structure (R_ST>>permuted R_ST) was only significant in comparisons between Central and South America. Further, we assessed local scale spatial genetic structure (SGS, d ≤ 8 km) in Panama and developed an agent-based model parameterized with data from F. insipida to estimate minimum pollination distances, which determine the contribution of pollen dispersal on SGS. The local scale data for F. insipida was compared to SGS data collected for an additional free-standing fig, F. yoponensis (subgenus Pharmacosycea), and two species of strangler figs, F. citrifolia and F. obtusifolia (subgenus Urostigma) sampled in Panama. All four species displayed significant SGS (mean Sp = 0.014 ± 0.012). Model simulations indicated that most pollination events likely occur at distances > > 1 km, largely ruling out spatially limited pollen dispersal as the determinant of SGS in F. insipida and, by extension, the other fig species. Our results are consistent with the view that Ficus develops fine-scale SGS primarily as a result of localized seed dispersal and/or clumped seedling establishment despite extensive long-distance pollen dispersal. We discuss several ecological and life history factors that could have species- or subgenus-specific impacts on the genetic structure of Neotropical figs.

The incidence and pattern of copollinator diversification in dioecious and monoecious figs

Differences in breeding system are associated with correlated ecological and morphological changes in plants. In Ficus, dioecy and monoecy are strongly associated with different suites of traits (tree height, population density, fruiting frequency, pollinator dispersal ecology). Although approximately 30% of fig species are pollinated by multiple species of fig-pollinating wasps, it has been suggested that copollinators are rare in dioecious figs. Here, we test whether there is a connection between the fig breeding system and copollinator incidence and diversification by conducting a meta-analysis of molecular data from pollinators of 119 fig species that includes new data from 15 Asian fig species. We find that the incidence of copollinators is not significantly different between monoecious and dioecious Ficus. Surprisingly, while all copollinators in dioecious figs are sister taxa, only 32.1% in monoecious figs are sister taxa. We present hypotheses to explain those patterns and discuss their consequences on the evolution of this mutualism.

Dioecy, more than monoecy, affects plant spatial genetic structure: the case study of Ficus

In this analysis, we attempt to understand how monoecy and dioecy drive spatial genetic structure (SGS) in plant populations. For this purpose, plants of the genus Ficus were used as a comparative model due to their particular characteristics, including high species diversity, variation in life histories, and sexual systems. One of the main issues we assessed is whether dioecious fig tree populations are more spatially genetically structured than monoecious populations. Using the Sp statistic, which allows for quantitative comparisons among different studies, we compared the extent of SGS between monoecious and dioecious Ficus species. To broaden our conclusions we used published data on an additional 27 monoecious and dioecious plant species. Furthermore, genetic diversity analyses were performed for two monoecious Ficus species using 12 microsatellite markers in order to strengthen our conclusions about SGS. Our results show that dioecy, more than monoecy, significantly contributes to SGS in plant populations. On average, the estimate of Sp was six times higher for dioecious Ficus species than monoecious Ficus species and it was two times higher in dioecious than monoecious plant species. Considering these results, we emphasize that the long-distance pollen dispersal mechanism in monoecious Ficus species seems to be the dominant factor in determining weak spatial genetic structure, high levels of genetic diversity, and lack of inbreeding. Although Ficus constitute a model species to study SGS, a more general comparison encompassing a wider range of plants is required in order to better understand how sexual systems affect genetic structure.

Transferability and characterization of simple sequence repeat markers from Anacardium occidentale to A. humile (Anacardiaceae)

Use of molecular markers can be limited by the high cost and extensive time required for their development. Transfer of simple sequence repeat (SSR) markers reduces the cost and time limitations and has allowed the use of these markers in a larger number of species. We tested 11 SSR markers previously developed for Anacardium occidentale on A. humile. The 11 loci were successfully amplified in A. humile. All loci were polymorphic and generated a mean of 5.4 alleles per locus. The observed heterozygosity was lower than the expected heterozygosity under Hardy-Weinberg equilibrium for most loci, with mean values of 0.463 and 0.696, respectively. The endogamy coefficients were positive and significant for seven loci. However, the combined probability of paternity exclusion was high, and the combined probability of genetic identity was low. None of the pairs of loci were in linkage disequilibrium. The informative power of these loci demonstrates that they are suitable for studies of diversity and genetic structure of natural populations of A. humile. In addition, the loci are suitable for estimating gene flow between populations, assessing species crossing preferences, and performing interspecific comparisons.

Geographic variation of floral scent in a highly specialized pollination mutualism

Floral scents are important signals for communication between plants and pollinators. Several studies have focused on interspecific variation of these signals, but little is known about intraspecific variation in flower scent, particularly for species with wide geographic distributions. In the highly specific mutualism between Ficus species and their pollinating wasps, chemical mediation is crucial for partner encounter. Several studies show that scents, i.e. blends of volatiles, are species-specific, but no studies address interpopulation variation of scents in fig pollination mutualisms, which often have broad geographic distributions. In this study, using absorption/desorption headspace techniques, we analyzed variation in floral scent composition among three populations of each of two widely distributed Asian Ficus species. We identified more than 100 different volatile organic compounds, predominantly terpenes. In both species, significant differences were found between scent bouquets of East Asian and Indian populations. These differences are discussed in relation to geographical barriers that could disrupt gene exchange between these two areas, thereby isolating Indian populations from those of Eastern Asia.