Xiaomei Wang
Shuai Liao2,3
Zhen Zhang
Hongqing Li- 1School of Life Sciences, East China Normal University, Shanghai, China
- 2Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- 3South China National Botanical Garden, Guangzhou, China
- 4College of Architecture and Urban Planning, Tongji University, Shanghai, China
- 5School of Life and Environmental Sciences, Shaoxing University, Shaoxing, China
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.
1 Introduction
Species complexes, as species groups characterized by overlapping morphologies and incomplete genetic differentiation, emerge through various speciation processes, including hybridization, polyploidization, introgression, and incomplete lineage sorting (Wang, 1998; Nosil et al., 2009; Stankowski and Ravinet, 2021). Given that speciation is a protracted and continuous phenomenon, the species within a complex often stem from diverse speciation events or stage of speciation processes. Therefore, linking microevolutionary mechanisms (such as polyploidization and hybridization) with macroevolutionary patterns (like speciation) is vital for a thorough understanding of their speciation and evolutionary histories (Nosil et al., 2009; Christie and Strauss, 2018; Pinheiro et al., 2018; Lu et al., 2022). However, phylogenetic and evolutionary interpretations based on morphology in species complexes frequently encounter uncertainties. This is due to intrinsic or extrinsic factors that may induce morphological variation resembling each other. Morphological traits are often continuous, but their diagnosis is limited by taxonomic thresholds. Additionally, the rapid, repeated evolution of adaptive phenotypes further complicates the situation (Fišer et al., 2018; Prata et al., 2018; Galtier, 2019). The paucity of reliable molecular markers at lower taxonomic levels further hinders the evolutionary inference within species complexes (Avise, 2012; Galtier, 2019). Moreover, challenges such as hybridization, incomplete lineage sorting, and polyploidization often result in that gene trees inadequately represent true species trees (Degnan and Rosenberg, 2006; Yang and Rannala, 2010; Villaverde et al., 2018). To address these issues, we intent to obtain whole-genome data, together with morphological, cytological, and ecological evidences. This multidisciplinary strategy may be pivotal in accurate species delimitation within a complex, elucidating the phylogenetic relationships of morphological variations or taxa.
Taxonomically complex groups (TCGs) often exhibit characteristics that complicate their classification into discrete species, such as uniparental reproduction, diverse ploidy levels, and hybridization. The genus Ficus (Moraceae) is a notable example of a TCG, primarily owing to its significant intraspecific variability (Ennos et al., 2005; Mahima et al., 2020). This genus is distinguished by extensive species diversity and a rich evolutionary history, especially influenced by specific symbiotic pollinators, leading multiple morphologically indistinct species complexes (Berg, 2007; Wang et al., 2016; Su et al., 2022). For example, the F. obliqua–F. rubiginosa complex consists of three species, all of which are pollinated by closely related fig wasps (Wiebes, 1994; Dixon et al., 2001). The F. sarmentosa complex, encompassing at least 21 climbing taxa adapted to limestone habitats, is characterized by widespread interspecific hybridization, resulting in highly convoluted phylogenetic relationships (Zhang et al., 2023). Similarly, the F. auriculata complex (Zhang et al., 2018) and the “hairy-fig” complex (Lu et al., 2016) both display extensive genetic sharing and morphological continua among their respective groups. The F. petiolaris complex, endemic to Mexico, is distinguished by its morphological continuity, leading to its classification as a single species, F. petiolaris (Piedra-Malagón et al., 2011). These instances underscore that classification within Ficus can be significantly influenced by frequent hybridization and convergent morphological trends.
The F. erecta complex, along with its allies, constitutes a complex species group within subsect. Frutescentiae of subg. Ficus, predominantly distributed across the subtropical regions of China. These taxa are characterized by their dioecious syconia, predominantly small shrub-like habit, the distinct Terminalia-branching pattern, and axillary figs (Berg et al., 2005). Their foliage is typically narrow and papery, and may exhibit variations in pilosity. The accurate taxonomic identification within this group encounters considerable hurdles due to the scarcity of comprehensive specimen sets (including diverse morphological forms, both genders, and syconia at different developmental stages) and pronounced morphological plasticity, which has notably hindered advancements in the taxonomic, systematic, and evolutionary ecological researches of this group. The most extensive study hitherto conducted by Lu et al. (Lu et al., 2017) employed nuclear genes (ITS+ETS) and simple sequence repeats (SSR) to analyze the phylogeny of approximately 30 taxa within this subsection. This seminal work delineated a core clade, with limited internal branch support and characterized by overlapping morphological variations, as the F. erecta complex, comprising 14 species and 3 cryptic species, culminating in 17 taxonomic entities (F. abelii, F. boninsimae, F. erecta, F. fengkaiensis, F. formosana, F. fusuiensis, F. gasparriniana var. laceratifolia, F. gasparriniana var. viridescens, F. iidaiana, F. nishimurae, F. pandurata, F. periptera, F. pyriformis, F. sinociliata, F. stenophylla, F. tannoensis and F. vaccinioides). However, extensive comparative morphological anatomical investigations have yet to discover distinctive unique features which are capable of reliably differentiating most species within the complex (Zhou and Gilbert, 2003; Gui, 2013; Gao, 2018). This issue is particularly evident in taxa exhibiting varietal morphological diversity, such as F. gasparriniana and F. pandurata. Besides, a series of studies have underscored the prevalence of hybridization and pollinator sharing events within the complex. Notably, early hybridization events have been identified between F. erecta and some species endemic to Ogasawara Islands in Japan, such as F. boninsimae, F. nishimurae, and F. iidaiana. Additionally, F. erecta is presumed to share pollinators with F. tannoensis and F. vaccinioides from Taiwan, as well as the more widely distributed F. formosana (Kusumi et al., 2012; Wachi et al., 2016; Su et al., 2022). Nevertheless, the phylogenetic relationships within the complex remain largely unresolved due to limitations in sampling adequacy or the paucity of informative genetic loci. This is particularly evident in the ambiguous relationships among the varieties of F. gasparriniana, the distinct genetic diversification of F. stenophylla, and the weak support for existing phylogenetic branches (Lu et al., 2017). Consequently, the delineation of species within the complex is still in need of a more robust and comprehensive approach.
In Ficus, typical diploid plants possess a chromosome count of 2n = 26, as established in prior studies (Ohri and Khoshoo, 1987; Guasmi et al., 2006). Within the F. erecta complex, the occurrence of polyploidy in multiple species has been documented, further complicating taxonomy and obscuring the understanding of interspecific phylogenetic relationships and evolutionary history. Notably, F. iidaiana has been identified as a tetraploid (Ono, 1990; Kusumi et al., 2012). Additionally, genome sequencing has indicated the existence of tetraploid individuals in the F. erecta population, with genome size approximating 590 Mbp (Shirasawa et al., 2020). The presence of such polyploidy events not only complicates taxonomy of the F. erecta complex but also obscures the understanding of the interspecific phylogenetic relationships and its evolutionary history.
Chloroplast genomes, functioning as super-barcodes, are pivotal for species identification and are instrumental in discovering new taxa and cryptic species (Ji et al., 2019; Liu et al., 2021). Additionally, reduced-representation genomic techniques, such as restriction-site associated DNA sequencing (RADseq), have been widely employed in phylogenetic studies, particularly for their efficacy in elucidating complex reticulate evolutionary events, including introgression and hybridization in closely related groups (Eaton and Ree, 2013; Eaton et al., 2015; Razkin et al., 2016; Maguilla et al., 2017; Paetzold et al., 2019). In this study, we employed skimming sequencing to acquire complete chloroplast genomes, complemented by single nucleotide polymorphisms (SNPs) derived from RADseq. This dual genomic approach is aimed at reconstructing a robust phylogenetic tree and then detect hybridization and introgression events for the F. erecta complex. Combining these phylogenomic analyses and extra chromosome ploidy detection, the study finally tried to figure out cryptic species, the extent of species differentiation, and interspecific phylogenetic relationships within the F. erecta complex. This comprehensive phylogenomic analysis of the F. erecta complex is pivotal not only for unraveling its complex evolutionary history but also for a broader understanding of plant speciation and evolutionary dynamics. Moreover, this study sets a precedent for the meticulous delineation of taxa within taxonomically challenging groups.
2 Materials and methods
2.1 Sampling strategy
The sampling strategy was meticulously devised to capture the extensive morphological diversity inherent in the complex. Particular attention was focused on species ambiguously classified as “the varieties of Ficus gasparriniana”, “F. stenophylla”, “F. pandurata” and specimens presenting morphological identification challenges within the complex. A representative set of 3–5 samples was collected for each species, with a minority represented by a single sample. Samples from species with restricted distributions were predominantly sourced from their respective type localities. This study successfully acquired samples from all taxa within the complex, with the exception of three cryptic species exclusively documented from type specimens in Southwest China and endemic species from Japan’s Ogasawara Islands. Informed by Lu et al.’s study (Lu et al., 2017), this study further incorporated representative samples from six related species, which are potential candidates for hybridization events within the complex. Evolutionary analyses were conducted on all samples, employing RADseq and genome skimming sequencing techniques. Individuals exhibiting markedly larger leaves and figs, indicative of potential polyploidy, were specifically selected for chromosome ploidy examination. The voucher specimens have been archived in the Herbarium of East China Normal University (HSNU). Detailed information on the samples is systematically cataloged in Figure 1 and Table S1. (The sample name marked with a question mark, indicates that it is a cryptic species. Cryptic species refers to taxa that cannot be readily distinguished morphologically, yet evidence indicates they are on different evolutionary trajectories.).
Figure 1 Geographical distribution map of the collected experimental materials.
2.2 DNA extraction, library construction, and sequencing
Total DNA was extracted utilizing a modified CTAB protocol (Doyle and Doyle, 1987), and its integrity was confirmed through agarose gel electrophoresis and detected by NanoDrop 2000 Spectrophotometer. Subsequently, DNA samples meeting quality standards were sent to BGI (Shenzhen, China) for library construction and sequencing. The chloroplast genome sequencing required a minimum DNA quantity of 1.5 μg with a concentration threshold of 20 ng/μL, while RADseq library sequencing stipulated a DNA minimum of 1 μg and a concentration of at least 12.5 ng/μL.
In the skimming sequencing approach, DNA purification was followed by fragmentation to generate short-insert (≤800 bp) paired-end reads. Sequencing tags facilitated the indexing of each sample’s genomic DNA for library construction. Whole-genome skimming sequencing was executed on the HiSeq™ 2000 platform (Illumina, San Diego, California, USA), yielding a minimum of 2 G of raw data per sample with 125 bp pair-end reads.
For RADseq library construction, EcoRI was selected based on pre-digestion experiments for enzymatic digestion, and sample-specific barcodes were employed to prepare RAD libraries. Sequencing was performed on the HiSeq™ 2000 platform (Illumina, San Diego, California, USA) utilizing PE150 paired-end sequencing, ensuring a data yield of no less than 1 G per sample. Sequencing data underwent quality control via FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) (Andrews, 2010). For subsequent assembly and analysis, read1 was chosen.
2.3 Data assembly
The initial trials with various assembly softwares revealed GetOrganelle as the superior tool in terms of assembly quality (Freudenthal et al., 2020). Thus, all samples were assembled using GetOrganelle (Jin et al., 2020), with F. religiosa (GenBank accession no. NC_033979) serving as the reference genome. The assembly parameters were meticulously optimized as followed: K-mer values were set to 55, 75, 95, 115, and 127 (Jin et al., 2020), with a word size of 85 and R at 15. For samples that did not automatically circularize, manual circularization was performed in Bandage v.0.8.0 (Wick et al., 2015), referencing other complete circular sequences from the same species to bridge gaps. Sequence alignments were executed using MAFFT v7.490 (Katoh and Standley, 2013) under default settings. Chloroplast genome annotation was conducted via Plastid Genome Annotator (PGA) software (https://github.com/quxiaojian/PGA) (Qu et al., 2019), using F. religiosa (NC_033979) (Bruun-Lund et al., 2017) and F. sarmentosa (NC_061976) (Zhang et al., 2022) as reference genomes. The tRNA genes’ boundaries were accurately determined using the tRNAscan-SE server (http://lowelab.ucsc.edu/tRNAscan-SE/) (Chan and Lowe, 2019). The whole plastid genome was aligned and used for phylogenetic analysis, with one of the two inverted regions and psbZ, and psbZ-trnG-GCC (The chloroplast genomes of three samples is incomplete, with one gap in this region.) removed.
Considering its robustness in processing indels and estimating heterozygosity (H) and sequencing error rates (E), ipyrad (Eaton and Overcast, 2020) was selected for RADseq data assembly using the single end of the paired-end sequences (R1) after control via FastQC. The data filtering was conducted in ipyrad through replacing bases with a Phred quality score (Q) below 33 with an ambiguous base (“N”), discarding reads exceeding 10% “N”. Clustering, based on sequence similarity, was conducted using UCLAST algorithm in USEARCH software tool (Edgar, 2010), followed by estimating individual hybrid rate (H) and sequencing error frequency (E) (Lynch, 2008), filtering out paralogous and highly repetitive sequences to attain consensus sequences. To mitigate the clustering values’ and sample coverage’s impacts on the phylogeny, we tested cluster_thread threshold of 0.85, 0.90, and 0.95, with 17 min_sample_loci gradients set from low to high under each threshold. This approach generated 51 sets of unlinked SNP datasets for the phylogenetic tree reconstruction. Consistency across phylogenetic tree topologies derived from all 51 datasets was observed, with only minor variations in shallow nodes and support values. The optimal data matrix for subsequent analyses was identified by calculating the Robinson-Foulds distance between different trees using RF.dist function in the phangorn package in R (Robinson and Foulds, 1981). The dataset with cluster_thread at 0.90 and min_sample_loci at 31 (i.e., c90min31), which exhibited the minimum Robinson-Foulds distance relative to others (Table S2), was selected for the phylogenetic tree construction and subsequent analyses (Table S3).
2.4 Phylogenetic analysis
Phylogenetic trees were reconstructed using both Maximum Likelihood (ML) and Bayesian Inference (BI) approaches with both the plastid and SNP datasets, as implemented in PhyloSuite (Zhang et al., 2020). For ML analysis, IQ-TREE 2 (Minh et al., 2020) was employed. Node support was assessed through 10,000 ultrafast bootstrap (UFboot) replications. ModelFinder (Kalyaanamoorthy et al., 2017), integrated within IQ-TREE 2, was utilized for optimal base substitution model selection based on the Bayesian Information Criterion (BIC) (Luo et al., 2010). Bayesian phylogenetic trees were generated using Mrbayes 3.2.6 (Ronquist et al., 2012). The Markov Chain Monte Carlo (MCMC) parameters were set as follows: total generations (mcmcp ngen) at 3,000,000 for plastid and 10,000,000 for SNPs datasets, sampling frequency (samplefreq) at every 1,000 generation, and four simultaneous runs (nruns). The initial 25% of samples were discarded as burn-in, and the remaining samples were used to calculate posterior probabilities (PP). The ModelFinder was again used to determine the optimal base substitution model, constrained to the 24 models compatible with MrBayes (MrBayes model parameters: -m TEST -mset mrbayes). F. pumila was designated as the outgroup.
Phylogenetic trees derived from the chloroplast whole-genome dataset and the RADseq dataset (c90min31 dataset, as shown in Table S3) were compared using a mirror-image approach. Visualization was conducted in R utilizing the phytools, ape, and maps packages. Node support was categorized into high (ML/UFboots ≥ 98 or BI/PP ≥ 0.98), moderate (98 > ML/UFboots ≥ 95 or 0.98 > BI/PP ≥ 0.90), and low/no support (ML/UFboots < 95 or BI/PP < 0.90).
Another phylogenetic tree was constructed using the Tetrad tool (https://ipyrad.readthedocs.io/en/latest/API-analysis/cookbook-tetrad.html) within ipyrad, which based on the SVDQuartets algorithm of Chifman and Kubatko (Chifman and Kubatko, 2014), following a coalescent-based approach. The analysis involved all quartets, with 100 bootstrap replications. F. pumila served as the outgroup. This method allowed the inference of both the optimal tree and the majority-rule consensus tree, along with the calculation of branch support.
2.5 Analysis of population genetic structure
In this study, the population genetic structure of the F. erecta complex was analyzed using ADMIXTURE v1.3.0 (Alexander et al., 2009). The analysis was tailored to precisely infer the genetic structure of the F. erecta complex, thus, the number of inferred ancestral populations (K) was methodically varied from 2 to 18, with 1000 bootstraps using default optimization method. The determination of the optimal K value was based on the minimization of the cross-validation error (CV error). The ADMIXTURE results, including the best and second-best K values, were visualized using R.
2.6 Ploidy level determination
Leaf samples, dried using silica gel and weighing 20–50 mg, were processed for ploidy analysis. Each sample was finely chopped in 1 mL of LB01 buffer. Subsequently, the supernatant, containing cellular material, was gently aspirated and filtered through a 400-mesh nylon net into a 2 mL centrifuge tube. Post-centrifugation, the supernatant was discarded, and the cell pellet was stained with 80 μg/mL propidium iodide (PI) and 80 μg/mL RNase for 30 min in an ice bath. Diploid F. pumila (CNA0019284), with a genome size of 315.7 Mbp, sample ID 20200034_1, was utilized as an external reference (Wang et al., 2021). The stained cells samples were analyzed using a MoFlo-XDP high-speed flow cytometer (Beckman Coulter Inc., USA), operating at a sheath pressure of 60 psi with a 70 μm ceramic nozzle. Nucleic acid PI fluorescence was detected employing a solid-state laser (488 nm) and a 530-/40-nm HQ bandpass filter. Each sample was analyzed for a minimum of 100 cells to ensure statistical relevance. Samples exhibiting atypical fluorescence values underwent a minimum of two analyses to ascertain result reliability.
3 Results
3.1 Phylogenetic tree reconstruction of the complex
For the phylogenetic analysis, RADseq and genome skimming data were obtained from 65 samples. The RADseq dataset, specifically the c90min31 (cluster_thread 0.90 and min_sample_loci 31) unlinked-SNPs matrix, spanned a length of 33,808 bp, encompassing 8,235 parsimony-informative sites. The optimal nucleotide substitution models of this dataset for ML (Maximum Likelihood) and BI (Bayesian Inference) phylogenetic analyses, determined by the BIC (Bayesian information cirtation), were TVM+R2+F+ASC and GTR+F+ASC, respectively. The chloroplast genome dataset presented a matrix length of 137,043 bp, featuring 5,708 variable sites and 980 parsimony-informative sites. The models deemed most suitable for ML/BI analyses of this dataset, following the BIC guidelines, were K3Pu+F+I+I+R3 and GTR+F +G4, respectively. The resulting phylogenetic trees, derived from RADseq and cpDNA datasets, are illustrated in Figure 2. Additionally, the tetrad tree based on the RADseq data is depicted in Figure 3.
Figure 2 Phylogenetic relationships within the F. erecta complex. The left panel portrays the tree derived from the RAD c90min31 dataset, while the right panel is based on the complete chloroplast genome dataset. Branch supports are denoted as ML tree/UFboots above and BI/PP below. Sample coloration is consistent with Figure 4, with diploids in green, triploids in blue, tetraploids in red, and uncertain ploidy levels in lighter hues.
Figure 3 The tetrad tree of the F. erecta complex. Clades c–q are in alignment with those delineated in Figure 2. Arrows indicate the samples whose phylogenetic positioning diverges from those in the ML/BI phylogenetic tree, as depicted in Figure 2.
Most “species” in the nature are on the way to the final speciation stage. The integrative species concept combining morphology and phylogeny is used in this article. Species have unique identifiable features. The RADseq-derived phylogenetic tree exhibited robust topology, with high support values across most nodes. Clade a encompassed samples from almost all taxa, barring F. ischnopoda and partial samples of F. periptera. Monophyly was observed in subsets of species, such as F. abelii (clade o), F. erecta (clade l), F. formosana (clade i), F. tannoensis (clade k), and F. vaccinioides (clade j), each with uniform ML/BI support of 100/1. F. pyriformis (clade d) forms a singular branch with a support value of 100/1. However, its variety, F. pyriformis var. hirtinervis, clustered with F. abelii (clade o) under high support, sharing traits like yellowish leaf undersides and densely haired young branches and petioles. F. pandurata (including var. angustifolia) diverged into two parts: clade h (support value100/1), featuring typical specimens with thinner leaves and fruits lacking white surface lumps, and clade g (support rate 100/1), characterized by thicker leaves and fruits with white lumps. F. stenophylla bifurcated into two highly supported branches (clades p and n), aligning with Lu et al., 2017. Clade p exhibited typical leaf morphology, sparse hair, and pronounced lateral veins on the leaf underside, contrasting with clade n (Tengchong, Yunnan) whose morphological traits were more akin to F. abelii (clade o). The samples of the three varieties of F. gasparriniana (var. esquirolii, var. laceratifolia and var. viridescens) were segregated into three distinct clades (e, f, m), with support values of 73/1, 100/1, and 100/1, in which only F. gasparriniana var. viridescens is monophyly (clade f). Clades e and f, exhibiting lumped and densely white-pubescent fruit surfaces, respectively, were closely related to F. heteromorpha but distant to clade m (robust plants) formed a sister relationship with one sample (20170011) of F. pandurata. Notably, one sample (20190041) of F. gasparriniana var. laceratifolia was sister to F. formosana (clade i).
In the phylogenetic tree derived from the complete chloroplast genomes, species within the F. erecta complex aggregated into two principal clades, clade A and clade B, both receiving robust support values of 100/1. Notably, only F. vaccinioides and F. tannoensis independently clustered within clade A, each achieving a unanimous support value of 100/1. F. stenophylla exhibited bifurcation, embedded within clade A (exhibiting leaf morphology akin to F. abelii, with a support value of 100/1) and clade B (representing the archetypal F. stenophylla, with a support value of 98/1). In contrast to the RAD-derived phylogenetic tree, the cpDNA tree positioned the samples corresponding to clade e (RAD tree) within clade A, where they formed a sister relationship with one sample (20170119) of F. heteromorpha (100/1). Meanwhile, the branches corresponding to clade f (RAD tree, 100/1) and clade m (RAD tree, 100/1) were situated within clade B. Furthermore, the samples of F. formosana and F. erecta were dispersed throughout clade A without clustering, whereas the samples of F.?pandurata (corresponding to clade g in RAD tree), F. periptera, and F. ischnopoda resided in clade B yet did not form distinct clusters. Additionally, several species, such as F. abelii, F. pandurata, F. pyriformis, and F. gasparriniana, had their samples distributed across both clade A and clade B.
The tetrad tree, generated from SVDQuartets algorithm (Chifman and Kubatko, 2014) via Tetrad tool in ipyrad (Eaton and Overcast, 2020), demonstrated remarkably high bootstrap support values across its structure (Figure 3), dividing into three primary limbs. When compared with the gene tree derived from RADseq data (Figure 2, left side), the tetrad tree revealed improved clustering for certain taxa. Notably, in the gene tree, F. pyriformis (clade d) and clade c were sister branches, together forming a sister group to F. stenophylla (clade p), which was, in turn, sister to F. periptera+F. ischnopoda, though the latter two did not form a distinct cluster. Conversely, in the tetrad tree, F. periptera and F. ischnopoda first clustered together and then formed a sister group with F. stenophylla (clade p). This collective cluster subsequently aligned as a sister group to F. pyriformis (clade d), ultimately constituting the sister group to clade c. Further distinctions in the tetrad tree included the sample (20170181) of F. pyriformis var. hirtinervis positioned distantly from F. abelii and the sample (20170119) of F. heteromorpha nested within clade e (Figure 3). Additionally, the sample (20170011) of F. pandurata did not cluster with clade m. Intriguingly, the sample (20190041) of F. gasparriniana var. laceratifolia clustered within clade m, comprising samples from the same species, rather than with clade i (Figure 3).
3.2 Genetic structure analysis of the F. erecta complex
CV error metrics (Table S4, Figure S1) identified 11 as the best-fit number of genetic clusters (K), with 12 as a secondary optimum (Figure 4). At K=11, distinct genetic backgrounds were evident for several species, including F. abelii, F. pyriformis, F. stenophylla, F. vaccinioides, and F. formosana, indicating well-defined genetic separations (Figure 4). Notably, F. erecta and F. pandurata samples displayed a shared genetic background (light purple), yet at K=12. The samples of F. gasparriniana (as depicted in clades e, f in Figures 2, 3), F.?pandurata, F. heteromorpha, and F. chapaensis shared a uniform genetic background (dark purple), hinting at close genetic relationships. Intriguingly, F.?stenophylla, morphologically akin to F. abelii, exhibited a genetic background significantly divergent from F. stenophylla. The genetic make-up of F. tannoensis samples was notably heterogeneous, reflecting an admixture of genetic contributions from F. vaccinioides and F. erecta. Furthermore, the genetic profiles of three F. periptera samples from Nonghua, Guangxi (20160380, 20190135), and Malipo, Yunnan (20190121) suggested hybrid origins. Additional samples exhibiting mixed genetic backgrounds included F. pyr. var. hirtinervis20170181, F. gas. var. laceratifolia20190041, F. pandurata20170011, and F. formosana20190220 (Figure 4).
Figure 4 ADMIXTURE results for K values of 11 and 12.
3.3 Analysis of chromosomal ploidy levels in the F. erecta complex
A comprehensive ploidy level analysis was conducted on 68 samples of the F. erecta complex, yielding 98 valid measurements (Table S5). To ensure data consistency, some samples were subjected to multiple tests. The number of cells counted per sample ranging from 95 to 4,746 in individual assays. Relative ploidy levels were calculated by comparing the average fluorescence mean values of each sample with those of the diploid reference species, F. pumila (CNA0019284). Figure 5 graphically represents these relative ploidy values. The majority of the F. erecta complex samples were diploid, such as F. pandurata (clade h, Figure 2), F.?pandurata (clade g), F. stenophylla (clade p), and F.?stenophylla (clade n). F. gasparriniana exhibited a complex ploidy pattern, with 2×, 3×, and 4× levels observed (Figures 2, 5). Samples from clade e and clade f (Figure 2) were consistently diploid, while clade m samples exhibited 2× and 3×, predominantly triploid. An unclustered sample (20190041) was identified as tetraploid. The chromosome ploidy of F. abelii and F. periptera was quite unique, with their relative ploidy values significantly lower than other species. The first one, F. abelii displayed two ploidy levels, with relative values of 1.26 and 2.46, suggestive of a 1:2 ratio. The second, F. periptera exhibited three distinct ploidy levels with relative values of 1.29, 1.84, and 2.51, approximating a 2:3:4 ratio. Whereas Ficus tannoensis and F. pedunculosa exhibited an approximate 4× ploidy level (Figure 5).
Figure 5 Relative ploidy values of the F. erecta complex. The y-axis represents the relative ploidy values, and the x-axis lists the sample names. The color coding (green for diploids, blue for triploids, and red for tetraploids) facilitates easy identification of the ploidy level, with lighter shades indicating uncertain ploidy.
4 Discussion and conclusion
4.1 Morphological variation and genetic differentiation in the F. erecta complex
The F. erecta complex is characterized by substantial intraspecific variability and indistinct interspecific boundary, leading to pervasive morphological overlap (Figure 6, which were integrated from type specimens, botanical descriptions, actual observations, and comparative taxonomic anatomy). The syconia (fig fruits) of the complex are generally uniform in morphology and size (1–2 cm), with inflorescences varying in shape from spherical to pear-shaped or spindle-shaped. Typically, the peduncles and ovary stalks exist or not. The lack of significant morphological variation among inflorescences and uniform phenological characteristics render leaf shape variations pivotal in species circumscription within the complex. However, reliance on foliar features alone often proves insufficient for accurate species identification.
Figure 6 Morphological relationships among various taxonomic units of the F. erecta complex and its close relatives.
Several widespread taxa within the complex, exemplified by F. erecta、F. formosana and F. pandurata, display considerable leaf shape polymorphism, mirroring similar variation patterns. Leaves predominantly range from elliptical to inversely elliptical, with an aspect ratio of less than 3, but gradually transition to elongated forms with aspect ratios exceeding 10 and vertical side veins (Figure 6, green and blue lines). Some leaves also exhibit split tips. Comparable variations are observed in closely related species like F. heteromorpha, F. pingtangensis, F. chapaensis, and F. daimingshanensis. Typical rheophytic species in the complex, including F. abelii、F. pyriformis、F. peritera、F. ischnopoda and F. stenophylla, share lanceolate leaf and sparse indumentum (exception: F. abelii), with elongated peduncle. However, these characteristics exhibited inconspicuous interspecific variation among them (Figure 6, green dotted line). F. gasparriniana, comprising multiple varieties, display the most pronounced variation, with F. gasparriniana var. esquirolii and var. laceratifolia primarily differentiated by leaf apex division or not. F. gasparriniana var. viridescens, with prominent indumentum, contrasts significantly from other varieties. Simultanously, these varieties of F. gasparriniana also shares morphological features with other groups in the complex. Additionally, differentiation in leaf texture and indumentum is evident within the complex, varying from papery to leathery or membranous and from dense to glabrous, yet these are insufficient for distinguishing the numerous taxonomic units. Therefore, careful consideration is necessary when identifying groups based on leaf morphology, demanding a comprehensive understanding of the group’s variation range (Keay, 1998).
Species within the F. erecta complex often exhibit sympatric or parapatric distributions, conducive to frequent gene flow (Whittemore and Schaal, 1991; Wang et al., 2021). Most species in the complex have similar ostiole sizes and overlapping fruiting phenology, facilitating the sharing of pollinating fig wasps. In extensively studied regions like Southeast China, pollinator sharing is commonplace. For instance, in Hong Kong, the pollinators of F. pyriformis and F. variolosa are almost indistinguishable (Hill, 1967), in Taiwan, F. formosana and F. vaccinioides share the pollinator Blastophaga nipponica (Wachi et al., 2016), and in Guangdong, five species including F. erecta, F. formosana, F. abelii, F. pyriformis and F. variolosa share the pollinator Blastophaga silvestriana (Su et al., 2022), leading to frequent hybridization within the F. erecta complex. Hybridization is also observed between closely related species like F. boninsimae and F. nishimurae (Yokoyama, 2003) (Figure 7, a summary of genetic relationships based on existing molecular evidence, including published results and the results of this study).
Figure 7 Genetic relationships among various taxonomic units within the F. erecta complex.
According to our study (Figures 2–4), high-frequency hybridization or gene flow events were also evident in taxa from other geographical areas, such as between F. abelii and F. pyriformis, F. periptera and F. ischnopoda, and F. formosana and F. gasparriniana in Southwest China. Both phylogenetictrees from SNPs datasets displayed several individuals in unusual positions (e.g., F. pyr. var. hirtinervis20170181, F. pandurata20170011), whose phylogenetic misplacement typically aligned with hybridization events (Figures 2, 3 (Bruun-Lund et al., 2017)) or genome polyploidization (Figure 5). Moreover, multiple ploidy levels were observed within several species (e.g., F. abelii, F. gasparriniana, F. pandurata and F. periptera). These findings collectively indicate blurred interspecific genetic boundaries in the F. erecta complex, leading to the formation of a tightly-knit species network (Figure 7). Nonetheless, there was a notable lack of alignment between the morphological relationship network (Figure 6) and the molecular systematic network (Figure 7). This divergence may be attributed to the limits of phylogenetic trees, which are impacted by gene incongruence arising from hybridization, horizontal gene transfer, and gene duplication and loss (Maddison, 1997; Massatti et al., 2016). The influence of phenotypic plasticity (West-Eberhard, 1989; Duminil and Di Michele, 2009; Fujita et al., 2012; Prata et al., 2018) and convergent evolution (Losos, 2011) may also contribute to these complex patterns.
4.2 Expansion of the F. erecta complex to encompass the entire subsect. Frutescentiae
The F. erecta complex was originally delineated based on a phylogenetically well-supported branch encompassing 11 taxonomic groups, augmented by selected cryptic species and island endemics, as per morphological characteristics (Lu et al., 2017). Lu et al., however, failed to establish a stable phylogenetic structure for taxa external to the complex, hinting at its potential species inclusions (Lu et al., 2017). Notably, F. gasparriniana is deemed to a member of the complex, while its variety, F. gasparriniana var. esquirolii, originally classified outside the complex in the study of Lu et al. In addressing this, our study incorporated additional samples from F. gasparriniana var. esquirolii, F. ischnopoda, F. heteromorpha, F. pingtangensis, F. chapaensis and F. pedunculosa, revealing their integration within the complex (Figures 2, 3). Specifically, F. gasparriniana var. esquirolii and F. heteromorpha clusters closely with other F. gasparriniana varieties (clade f), and F. pandurata (Figure 2, clade e, Figure 3, clade e). The substantial gene flow among these species confirmed their inclusion in the F. erecta complex. The affiliation of F. periptera, an established member of Ficus erecta complex, with the geographically and ecologically aligned F. ischnopoda, coupled with observed hybridization at the type locality (samples F. periptera20160380, F. periptera20190121 and F. periptera20190135), further supports expanding the complex’s scope (Figure 4). Moreover, F. ischnopoda’s alignment with the complex’s phylogenetic tree suggests its inclusion.
While the study did not encompass all related taxa in the subsect. Frutescentiae, such as F. oleifolia, F. variolosa, F. filicauda, and F. trivia, the phylogenetic proximity of the newly incorporated samples (e.g., sister group relationships) to these species had been evidenced in the previous study (Lu et al., 2017), notably the integration of F. pedunculosa (Figures 2, 3), previously deemed a sister group to the remaining subsection (Lu et al., 2017). Thus, it is worth believing that there is no clear genetic gap between the complex and allies out of the complex. Besides, there is no obvious morphological difference between the complex and allies out of the complex according to comparative taxonomic anatomy (Gao, 2018). Consequently, the F. erecta complex should be broadened to include the entire subsect. Frutescentiae. The previous studies have identified subsect. Frutescentiae as a robust monophyletic group within the genus Ficus, characterized by their dioecious syconia, predominantly small shrub-like habit, the distinct Terminalia-branching pattern, axillary figs, staminate flowers mostly scattered or near the ostiole,stamens 2 or 3, representing one of its nine primary branches (Berg, 2003; Gardner et al., 2023). Addressing taxonomic issues within any individual species of the complex necessitates a comprehensive understanding of the genetic backdrop of the entire subsection, highlighting the intertwined evolutionary narratives within this group.
4.3 Widespread polyploidy in the F. erecta complex
Our study uncovered the prevalence of polyploidization within the F. erecta complex. For instance, F. gasparriniana, along with its varieties, demonstrated a range of ploidy levels (2×, 3×, 4×). Tetraploidy has been detected in F. pandurata (F. pandurata 20170011), F. tannoensis and F. erecta. Relatively low nucleic acid PI fluorescence value together with a 1:2 ratio suggests possible diploidy and tetraploidy but with a smaller genome size in F. abelii (Figure 5). Similar low absolute value in F. periptera presumed to exhibit 2×, 3×, and 4× ploidy levels. Due to limited sampling in ploidy determination, other species within the complex may exhibit undetected polyploidy.
The peculiar phylogenetic positions of certain individuals, such as F. pandurata20170011 and F. gasparriniana var. laceratifolia20190041, which did not cluster with other samples of their species (Figures 2, 3), can be attributed to their allopolyploid nature. This phenomenon exemplifies the role of allopolyploidization in promoting rapid species differentiation (Yang et al., 2023) and complex evolution (Hörandl, 2022; Karbstein et al., 2022; Yang et al., 2023). The genome duplication of unreduced gametes from diploids plants through a triploid bridge is a key pathway of polyploids, as demonstrated by triploids in F. gasparriniana (clade m) and F. periptera (Figures 2, 5), suggesting ongoing polyploid formation in these taxa. Polyploidy occurring within multiple species resulting in species with multiple ploidy levels. Furthermore, considering polyploidization and/or hybridization often coexist with apomixis (Hojsgaard and Hörandl, 2019; Hörandl, 2022), the presence of apomixis in species like F. abelii, F. erecta, F. gasparriniana (clade m), F. pandurata, F. pedunculosa, F. periptera, and F. tannoensis is in need of further investigation. Our unpublished observations indicate the apparent occurrence of apomixis in F. gasparriniana (clade m).
4.4 Phylogenetic relationships among taxa and taxonomic resolution in the F. erecta complex
4.4.1 Ficus gasparriniana contains at least three independent lineages
Phylogenetic analyses revealed the subdivision of F. gasparriniana into three distinct clades, designated as clades e, f, and m, which do not correspond directly with the three recognized varieties of the species. Clades e and f exhibit closer phylogenetic affinity, contrasting with the relatively distant clade m. The ploidy levels within these clades vary, with clades e and f being diploid and clade m predominantly triploid with containing diploid individuals. Morphologically, clade e parallels F. heteromorpha, characterized by thin papery leaves, fruits with white lumps on the surface, and reddish-brown young branches and petioles. In contrast, clade f exhibits leafs densely covered with soft white hair, fruits similarly enveloped with soft white hair devoid of lumps, and green young branches and petioles, underscoring clear morphological distinctions between these two clades. The cpDNA-based phylogenetic tree segregates clade e within clade A, whereas clades f and m are nested within clade B, suggesting substantial genetic divergence between clades e and f. Clade m not only differs genetically and in ploidy levels but also manifests larger plant, leaf, and fruit sizes, coupled with evidence of apomixis, markedly distinguishing it from the former clades. Given the absent samples of F. gasparriniana var. gasparriniana, it can be inferred that F. gasparriniana encompasses at least three independent lineages, confirming its status as a polymorphic species. However, the redefinition of F. gasparriniana is contingent upon an extensive elucidation of the phylogenetic intricacies within the complex. This revelation about F. gasparriniana underscores the complex taxonomic intricacies within the F. erecta complex. The presence of distinct genetic lineages and morphological variabilities necessitates meticulous and comprehensive investigations for accurate species delineation and an in-depth understanding of their evolutionary histories.
4.4.2 Ficus pandurata contains two lineages
Ficus pandurata bifurcated into two distinct lineages: clade h and clade g, as revealed through RADseq-based analysis, both exhibiting diploid chromosomal structures (Figure 2). Genetic structuring analysis indicated that clade h shared a genetic affinity with F. erecta, while clade g aligned more closely with the genetic backgrounds of F. gasparriniana (clades e and f) and F. heteromorpha. A detailed morphological assessment revealed divergent characteristics between these clades. Clade h, encompassing F. pandurata and its varieties, is characterized by a thinner leaf texture and an absence of white lumps on the fruit surface. Conversely, clade g exhibits a thicker leaf morphology and typically bears fruits with tuberculate surfaces, delineating clear morphological distinctions from F. pandurata and its varieties, thus clade g represents a hitherto undescribed cryptic species within the F. erecta complex.
4.4.3 “Ficus stenophylla” samples from Tengchong, Yunnan (clade n), potentially represent a novel species
In this study, samples identified as F. stenophylla from Tengchong, Yunnan (clade n), exhibit a distinct phylogenetic lineage in the RADseq-derived tree, divergent from the canonical F. stenophylla (clade p) and showing closer affinity to F. abelii. This observation corroborates the findings of Lu et al (Lu et al., 2017), which utilized two nuclear markers (ITS+ETS). In the cpDNA phylogeny, F. abelii samples are distributed across Clade A and Clade B, whereas the F. stenophylla samples from clade n and clade p cluster distinctly in Clade A and Clade B, respectively (Figure 2). ADMIXTURE analysis at K=11 reveals divergent genetic compositions between two clades of F. stenophylla (n and p) and F. abelii (Figure 4). Morphological examination of clade n samples reveals an inverted ovate-lanceolate leaf shape, subdued lateral veins on the lower surface of blade, and elongated elliptical and often dehiscent syconia, bearing resemblance to F. abelii. Conversely, the prototypical F. stenophylla (clade p) displays linear-lanceolate leaves with pronounced lateral veins on the lower surface and slightly pubescent, elliptical to spherical syconia. These morphological distinctions, coupled with the phylogenetic and genetic data, strongly suggest that clade n represents a hitherto undescribed cryptic species within the F. erecta complex.
4.4.4 Discussion on the species status in other taxa of the F. erecta complex
The four species, F. tannoensis, F. vaccinioides, F. erecta, and F. formosana, each formed distinct monophyletic groups on both the RADseq-based phylogenetic trees (Figures 2, 3), constituting a bigger branch with strong support (100/1). This clustering suggests that they are closely related yet independent species. Notably, F. tannoensis exhibits hybrid origins, sharing genetic backgrounds with F. vaccinioide and F. erecta (Figure 4, K=11, K=12). Morphologically, F. tannoensis shares rare liana traits with F. vaccinioides, both native to southern Taiwan in China. However, its leaf shape variation is strikingly similar to F. erecta, albeit with noticeably smaller leaves. This is potentially influenced by its other parent, F. vaccinioides, which possesses the smallest leaves in the complex. Phylogenetically, F. tannoensis is positioned between F. vaccinioide and F. erecta. Furthermore, ploidy analysis confirmed tetraploidy in entire F. tannoensis samples in this study. Thus, it is plausible that F. tannoensis represents an allopolyploid species derived from F. vaccinioide and F. erecta, a finding that contributes novel insights into speciation processes within Ficus and highlights allopolyploidization as a mechanism driving its rich species diversity.
In the RADseq-derived phylogenetic tree, samples of F. ischnopoda and F. periptera failed to coalesce into monophyletic clusters. Intriguingly, some specimens from these taxa aligned as sister groups to other divergent clades within the complex (Figure 2). In stark contrast, the tetrad tree, constructed utilizing the SVDQuartets algorithm, presented a unified, well-supported branch comprising samples of both F. ischnopoda and F. periptera (Figure 3). Considered in light of the SVDQuartets algorithm is proficiency in resolving phylogenetic structures amid gene flow scenarios (Long and Kubatko, 2018), the tetrad tree topology should more accurately mirror the true evolutionary history of these taxa. This conjecture was further corroborated by the evident genetic background sharing between F. ischnopoda and F. periptera (Figure 4). Ploidy analysis underscores the occurrence of polyploidy phenomena (3× and 4×, Figure 5) exclusive to F. periptera. Coupled with the taxonomic challenge posed by their low morphological distinctiveness, the taxonomic statuses of these two species necessitates more rigorous scrutiny.
Previous molecular phylogenetic analyses suggested that F. iidaiana, F. nishimurae, and F. boninsimae, endemic to Japan’s Ogasawara Islands, exhibited close phylogenetic relationships with F. erecta (Azuma et al., 2010) F. iidaiana stands out with its larger stature, longer lamina and petioles, and more elongated fruits compared with F. nishimurae and F. boninsimae (Yokoyama and Iwatsuki, 1998; Yokoyama, 2006). It is also unique as the sole tetraploid among these species (Ono, 1990). F. nishimurae and F. boninsimae, while morphologically similar, exhibit divergences in some quantitative traits (such as leaf and fruit sizes) and habitat preferences. Notably, hybridization events between these two species have been documented (Kusumi et al., 2012). Despite their morphological similarities and instances of hybridization, each species maintains distinct characteristics, justifying their separate taxonomic statuses. However, this research did not include samples from these species, indicating a gap in comprehensive data that could further clarify their biogeographic histories within the complex.
Ficus fusuiensis resembles some variant individuals of F. pandurata, but F. fusuiensis can be distinguished by its cuneate leaf base and ovate-lanceolate leaf shape. Similarly, F. fengkaiensis and F. formosana demonstrate morphological parallels, yet F. fengkaiensis is characterized by broader leaves, elongated peduncle within the flowers. F. sinociliata closely resembles a variant form of F. erecta var. erecta. The primary distinguishing features of F. sinociliata include prominently marginal ciliate and reduced fruit size. It is noteworthy that only type specimens of these three species have been documented to date (Chang, 1983, Chang, 1984), with a lack of subsequent collection records. The observed high morphological variability, combined with extensive hybridization and the occurrence of polyploidization within the F. erecta complex, raises the possibility that these species may represent either hybrid progeny or occasional morphological variants of related species within the complex.
Chen et al. (Chen et al., 2020). recently described F. jingningensis, a novel species within the F. erecta complex. It features narrow lanceolate leaves with a cordate base based on both type specimens and our collectings from the type locality. However, a close comparison of leaf texture, indumentum patterns, and fruit morphology reveals a striking resemblance to F. erecta, suggesting that F. jingningensis may represent a specialized morphological variant within this species.
4.5 Ficus erecta complex is a rare and huge evolutionally joint species group with frequent reticulate evolution and polyploidization
Traditional molecular markers, typically sourced from singular organelles or limited nuclear DNA regions (Freeland et al., 2015), often fall short in reconstructing complex reticulate relationships due to constraints like uniparental inheritance, heterogeneous evolutionary history of the markers, and/or low variability. These challenges have historically hindered the precise delimitation of species within complexes (Rittmeyer and Austin, 2012; Kirschner et al., 2015; Rothfels, 2021).
This study, expanding upon Lu et al. (2017), leveraged vast cpDNA and RADseq data to amass a robust set of genetic markers. The resultant phylogenetic tree of the F. erecta complex exhibited pronounced nuclear-cytoplasmic conflicts (Figure 2) and notable discordances between SNPs trees which based on different tree construction methods (Figures 2, 3). This analysis not only identified multiple species in the complex with hybrid origins but also unraveled intricate patterns of reticulate evolution (Figure 7), thereby elucidating the phylogenetic relationships and taxonomic statuses of some taxa within the complex.
These phylogenomic instances evidence the faciliation of extensive RAD and cpDNA sequencing data in reconstructing phylogenetic relationships, identifying cryptic species, and resolving the complex evolutionary histories of species or taxonomic groups (Maddison, 1997; Zhou et al., 2018; Spriggs et al., 2019; Zhou et al., 2020). Additionally, this approach contributes to our understanding of the origins and evolutionary histories of polyploid species (He et al., 2023). Nevertheless, the study reveals challenges in achieving complete clustering of certain species on the phylogenetic trees, with several branches exhibiting low support values (Figure 2), and an overall lack of enhanced species resolution. This highlights that the greater the variability encompassed within a complex and the more comprehensive the sampling strategy, the more likely it is to reveal transitional morphological characteristics, hybridization events, and polyploidization phenomena. These factors collectively contribute to blurred species boundaries and amplified complexity in interspecific relationships (Duminil and Di Michele, 2009). Such scenarios underscore the conceptualization of a species complex as an assembly of interconnected and interdependent populations, each exhibiting unique morphological or genetic traits, yet demonstrating a continuum of variation across populations (Derkarabetian and Hedin, 2014; Hori et al., 2014). Extreme instances like the F. erecta complex, encompassing dozens of taxa, remain relatively rare. Future research, leveraging comprehensive sampling and population genomic data, is essential to elucidate the complex’s evolutionary history, speciation mechanisms, and species delineation. Such endeavors represent a contemporary approach to the longstanding biological challenge of species classification in the modern scientific era.
Data availability statement
The datasets presented in this study can be found in online repositories. The raw data for Ficus erecta complex generated for this study can be found at NCBI https://www.ncbi.nlm.nih.gov/sra/PRJNA1059135. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
XW: Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing. SL: Investigation, Resources, Supervision, Validation, Funding acquisition, Writing – review & editing. ZZ: Resources, Software, Writing – review & editing, Investigation. JZ: Resources, Validation, Writing – review & editing. LM: Data curation, Writing – review & editing. HL: Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Natural Science Foundation of China (Grant No. 31270242 and Grant No. 32300178).
Acknowledgments
We are grateful to Dr. Isaac Overcast at California Academy of Sciences for his support in the use of ipyrad, Dr Martin R. Smith at Durham University for his assistance in calculating the distance of phylogenetic trees, Dr. Jian-Jun Jin at Kunming Institute of Botany (Chinese Academy of Sciences) on chloroplast genome assembly, Mr. Su-Hua Yang at Institute of Botany (Chinese Academy of Sciences) for ploidy level determination, Dr. Zhi-Hui Su and Ms. Sasaki Ayako at Osaka University, Dr. Hsy-Yu Tzeng, PhD candidate Po-An Chou at National Chung Hsing University, Dr. Bo Pan and Dr. Gang Wang at Xishuangbanna Tropical Botanical Garden (Chinese Academy of Sciences), Ms. Hai-Yan Gao, Mei-Jiao Zhang, Yu-Ke Bi, Mr. Chen Fang and others for their support in wild collecting. The authors would like to thank TopEdit (www.topeditsci.com) for its linguistic assistance during the preparation of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1354812/full#supplementary-material
References
Alexander, D. H., Novembre, J., Lange, K. (2009). Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664. doi: 10.1101/gr.094052.109
Andrews, S. (2010). FastQC: a quality control tool for high throughput sequence data. (Cambridge, UK: Babraham Bioinformatics, Babraham Institute). Retrieved from https://www.bioinformatics. babraham.ac.uk/projects/fastqc/ [Accessed November 23, 2019].
Avise, J. C. (2012). Molecular markers, natural history and evolution (New York: Springer Science & Business Media).
Azuma, H., Harrison, R. D., Nakamura, K., Su, Z. H. (2010). Molecular phylogenies of figs and fig-pollinating wasps in the Ryukyu and Bonin (Ogasawara) islands, Japan. Genes Genet. Syst. 85, 175–192. doi: 10.1266/ggs.85.177
Berg, C. C. (2003). Flora Malesiana precursor for the treatment of Moraceae 3: Ficus subgenus Ficus. Blumea-Biodiversity Evol. Biogeography Plants 48, 529–550. doi: 10.3767/000651903X489537
Berg, C. C. (2007). Proposals for treating four species complexes in Ficus subgenus Urostigma section Americanae (Moraceae). Blumea-Biodiversity Evol. Biogeography Plants 52, 295–312. doi: 10.3767/000651907X609034
Berg, C. C., Corner, E. J. H., Nooteboom, H. P. (2005). “Flora malesiana. Series I, seed plants. Volume 17, part 2,” in Moraceae (Ficus) (Leiden: Nationaal Herbarium Nederland).
Bruun-Lund, S., Clement, W. L., Kjellberg, F., Rønsted, N. (2017). First plastid phylogenomic study reveals potential cyto-nuclear discordance in the evolutionary history of Ficus L. (Moraceae). Mol. Phylogenet. Evol. 109, 93–104. doi: 10.1016/j.ympev.2016.12.031
Chan, P. P., Lowe, T. M. (2019). tRNAscan-SE: searching for tRNA genes in genomic sequences. Methods Mol. Biol. 1962, 1–14. doi: 10.1007/978-1-4939-9173-0_1
Chang, S. S. (1984). New taxa of moraceae from China and Vietnam. Acta Phytotaxonomica Sin. 22, 64–76.
Chen, Z. H., Mei, X. D., Chen, F., Xie, W. Y., Li, G. Y. (2020). Ficus jingningensis (Moraceae), A new species from zhejiang. J. Zhejiang Forestry Sci. Technol. 40, 52–54. doi: 10.3969/j.issn.1001-3776.2020.05.009
Chifman, J., Kubatko, L. (2014). Quartet inference from SNP data under the coalescent model. Bioinformatics 30, 3317–3324. doi: 10.1093/bioinformatics/btu530
Christie, K., Strauss, S. Y. (2018). Along the speciation continuum: Quantifying intrinsic and extrinsic isolating barriers across five million years of evolutionary divergence in California jewelflowers. Evolution 72, 1063–1079. doi: 10.1111/evo.13477
Degnan, J. H., Rosenberg, N. A. (2006). Discordance of species trees with their most likely gene trees. PloS Genet. 2, e68. doi: 10.1371/journal.pgen.0020068
Derkarabetian, S., Hedin, M. (2014). Integrative taxonomy and species delimitation in harvestmen: a revision of the western North American genus Sclerobunus (Opiliones: Laniatores: Travunioidea). PloS One 9, e104982. doi: 10.1371/journal.pone.0104982
Dixon, D. J., Jackes, B. R., Bielig, L. M. (2001). Figuring out the figs: the Ficus obliqua-Ficus rubiginosa complex (Moraceae: Urostigma sect. Malvanthera). Aust. Syst. Bot. 14, 133–154. doi: 10.1071/SB99029
Doyle, J. J., Doyle, J. L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical bulletin. 19 (1), 11–15.
Duminil, J., Di Michele, M. (2009). Plant species delimitation: a comparison of morphological and molecular markers. Plant Biosyst. 143, 528–542. doi: 10.1080/11263500902722964
Eaton, D. A., Hipp, A. L., González-Rodríguez, A., Cavender-Bares, J. (2015). Historical introgression among the American live oaks and the comparative nature of tests for introgression. Evolution 69, 2587–2601. doi: 10.1111/evo.12758
Eaton, D. A., Overcast, I. (2020). Ipyrad: Interactive assembly and analysis of RADseq datasets. Bioinformatics 36, 2592–2594. doi: 10.1093/bioinformatics/btz966
Eaton, D. A., Ree, R. H. J. (2013). Inferring phylogeny and introgression using RADseq data: an example from flowering plants (Pedicularis: Orobanchaceae). Syst. Biol. 62, 689–706. doi: 10.1093/sysbio/syt032
Ennos, R. A., French, G. C., Hollingsworth, P. M. (2005). Conserving taxonomic complexity. Trends Ecol. Evol. 20, 164–168. doi: 10.1016/j.tree.2005.01.012
Fišer, C., Robinson, C. T., Malard, F. (2018). Cryptic species as a window into the paradigm shift of the species concept. Mol. Ecol. 27, 613–635. doi: 10.1111/mec.14486
Freeland, J. R., Kirk, H., Petersen, S. D. (2015). Molecular ecology (Hoboken: A John Wiley & Sons, Ltd), 35–75.
Freudenthal, J. A., Pfaff, S., Terhoeven, N., Korte, A., Ankenbrand, M. J., Förster, F. (2020). A systematic comparison of chloroplast genome assembly tools. Genome Biol. 21, 1–21. doi: 10.1186/s13059-020-02153-6
Fujita, M. K., Leaché, A. D., Burbrink, F. T., McGuire, J. A., Moritz, C. (2012). Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol. Evol. 27, 480–488. doi: 10.1016/j.tree.2012.04.012
Galtier, N. (2019). Delineating species in the speciation continuum: A proposal. Evol. Appl. 12, 657–663. doi: 10.1111/eva.12748
Gao, H. Y. (2018). Clustering Analyses of Ficus erecta Complex and its Allies Based on Comparative Taxonomic Anatomy (Shanghai (China: East China Normal University).
Gardner, E. M., Bruun-Lund, S., Niissalo, M., Chantarasuwan, B., Clement, W. L., Geri, C., et al. (2023). Echoes of ancient introgression punctuate stable genomic lineages in the evolution of figs. PNAS 120, e2222035120. doi: 10.1073/pnas.2222035120
Guasmi, F., Ferchichi, A., Fares, K., Touil, L. (2006). Identification and differentiation of Ficus carica L. cultivars using inter simple sequence repeat markers. Afr. J. Biotechnol. 5, 1370–1374.
Gui, P. (2013). Morphological and anatomical studies on ficus subg. Ficus (Moraceae) in China (Shanghai (China: East China Normal University).
He, L., Guo, F. Y., Cai, X. J., Chen, H. P., Lian, C. L., Wang, Y., et al. (2023). Evolutionary origin and establishment of a dioecious diploid-tetraploid complex. Mol. Ecol. 32, 2732–2749. doi: 10.1111/mec.16902
Hill, D. S. (1967). Figs (Ficus spp.) and fig-wasps (Chalcidoidea). J. Nat. Hist. 1, 413–434. doi: 10.1080/00222936700770401
Hojsgaard, D., Hörandl, E. (2019). The rise of apomixis in natural plant populations. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00358
Hörandl, E. (2022). Novel approaches for species concepts and delimitation in polyploids and hybrids. Plants 11, 204. doi: 10.3390/plants11020204
Hori, K., Tono, A., Fujimoto, K., Kato, J., Ebihara, A., Watano, Y., et al. (2014). Reticulate evolution in the apogamous Dryopteris varia complex (Dryopteridaceae, subg. Erythrovariae, sect. Variae) and its related sexual species in Japan. J. Plant Res. 127, 661–684. doi: 10.1007/s10265-014-0652-0
Ji, Y. H., Liu, C. K., Yang, Z. Y., Yang, L. F., He, Z. S., Wang, H. C., et al. (2019). Testing and using complete plastomes and ribosomal DNA sequences as the next generation DNA barcodes in Panax (Araliaceae). Mol. Ecol. 19, 1333–1345. doi: 10.1111/1755-0998.13050
Jin, J. J., Yu, W. B., Yang, J. B., Song, Y., DePamphilis, C. W., Yi, T. S., et al. (2020). GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 21, 241. doi: 10.1186/s13059-020-02154-5
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K., Von Haeseler, A., Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285
Karbstein, K., Tomasello, S., Hodač, L., Wagner, N., Marinček, P., Barke, B. H., et al. (2022). Untying Gordian knots: Unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex. New Phytol. 235, 2081–2098. doi: 10.1111/nph.18284
Katoh, K., Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Kirschner, J., Zaveska Drabkova, L., Štěpánek, J., Uhlemann, I. (2015). Towards a better understanding of the Taraxacum evolution (Compositae–Cichorieae) on the basis of nrDNA of sexually reproducing species. Plant Syst. Evol. 301, 1135–1156. doi: 10.1007/s00606-014-1139-0
Kusumi, J., Azuma, H., Tzeng, H. Y., Chou, L. S., Peng, Y. Q., Nakamura, K., et al. (2012). Phylogenetic analyses suggest a hybrid origin of the figs (Moraceae: Ficus) that are endemic to the Ogasawara (Bonin) Islands, Japan. Mol. Phylogenet. Evol. 63, 168–179. doi: 10.1016/j.ympev.2012.01.004
Liu, Z. F., Ma, H., Ci, X. Q., Li, L., Song, Y., Liu, B., et al. (2021). Can plastid genome sequencing be used for species identification in Lauraceae? Bot. J. Linn. Soc 197, 1–14. doi: 10.1093/botlinnean/boab018
Long, C., Kubatko, L. (2018). The effect of gene flow on coalescent-based species-tree inference. Syst. Biol. 67, 770–785. doi: 10.1093/sysbio/syy020
Losos, J. B. (2011). Convergence, adaptation, and constraint. Evolution 65, 1827–1840. doi: 10.1111/j.1558-5646.2011.01289.x
Lu, J., Gui, P., Lu, Z. L., Zhang, L. F., Tian, H. Z., Gilbert, M. G., et al. (2016). Phylogenetic analysis and taxonomic delimitation of the “hairy-fig” complex of Ficus sect. Eriosycea (Moraceae) in China. Phytotaxa 261, 121–136. doi: 10.11646/phytotaxa.261.2.2
Lu, W. X., Hu, X. Y., Wang, Z. Z., Rao, G. Y. (2022). Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China. Cladistics 38, 663–683. doi: 10.1111/cla.12514
Lu, Z. L., Zhang, Z., Zhou, Q. M., Lu, J., Tian, H. Z., Li, H. Q. (2017). Molecular analyses of Ficus erecta and its allies within the subsection Frutescentiae (Moraceae). Plant Syst. Evol. 303, 603–614. doi: 10.1007/s00606-017-1393-z
Luo, A. R., Qiao, H. J., Zhang, Y. Z., Shi, W. F., Ho, S. Y., Xu, W. J., et al. (2010). Performance of criteria for selecting evolutionary models in phylogenetics: a comprehensive study based on simulated datasets. BMC Evol. Biol. 10, 1–13. doi: 10.1186/1471-2148-10-242
Lynch, M. (2008). Estimation of nucleotide diversity, disequilibrium coefficients, and mutation rates from high-coverage genome-sequencing projects. Mol. Biol. Evol. 25, 2409–2419. doi: 10.1093/molbev/msn185
Maddison, W. P. (1997). Gene trees in species trees. Syst. Biol. 46, 523–536. doi: 10.1093/sysbio/46.3.523
Maguilla, E., Escudero, M., Hipp, A. L., Luceño, M. (2017). Allopatric speciation despite historical gene flow: Divergence and hybridization in Carex furva and C. lucennoiberica (Cyperaceae) inferred from plastid and nuclear RAD-seq data. Mol. Ecol. 26, 5646–5662. doi: 10.1111/mec.14253
Mahima, K., Sudhakar, J. V., Ramalingam, S. (2020). Taxonomic delimitation of endemic Ficus amplocarpa and Ficus dalhousiae Complexes (Moraceae) by DNA barcoding. Phytotaxa 436, 21–35. doi: 10.11646/phytotaxa.436.1.2
Massatti, R., Reznicek, A. A., Knowles, L. L. (2016). Utilizing RADseq data for phylogenetic analysis of challenging taxonomic groups: A case study in Carex sect. Racemosae. Am. J. Bot. 103, 337–347. doi: 10.3732/ajb.1500315
Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., Von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015
Nosil, P., Harmon, L. J., Seehausen, O. (2009). Ecological explanations for (incomplete) speciation. Trends Ecol. Evol. 24, 145–156. doi: 10.1016/j.tree.2008.10.011
Ohri, D., Khoshoo, T. N. (1987). Nuclear DNA contents in the genus Ficus (Moraceae). Plant Syst. Evol. 156, 1–4. doi: 10.1007/BF00937196
Ono, M. (1990). Chromosome number of some endemic plant species of the Bonin Islands, III. Ogasawara Res. 16, 1–27.
Paetzold, C., Wood, K. R., Eaton, D. A., Wagner, W. L., Appelhans, M. S. (2019). Phylogeny of Hawaiian Melicope (Rutaceae): RAD-seq resolves species relationships and reveals ancient introgression. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01074
Piedra-Malagón, E. M., Sosa, V., Ibarra-Manríquez, G. (2011). Clinal variation and species boundaries in the Ficus petiolaris complex (Moraceae). Syst. Bot. 36, 80–87. doi: 10.1600/036364411X553153
Pinheiro, F., Dantas-Queiroz, M. V., Palma-Silva, C. (2018). Plant species complexes as models to understand speciation and evolution: a review of South American studies. Crit. Rev. Plant Sci. 37, 54–80. doi: 10.1080/07352689.2018.1471565
Prata, E. M., Sass, C., Rodrigues, D. P., Domingos, F. M., Specht, C. D., Damasco, G., et al. (2018). Towards integrative taxonomy in Neotropical botany: disentangling the Pagamea guianensis species complex (Rubiaceae). Bot. J. Linn. Soc 188, 213–231. doi: 10.1093/botlinnean/boy051
Qu, X. J., Moore, M. J., Li, D. Z., Yi, T. S. (2019). PGA: a software package for rapid, accurate, and flexible batch annotation of plastomes. Plant Methods 15, 50. doi: 10.1186/s13007-019-0435-7
Razkin, O., Sonet, G., Breugelmans, K., Madeira, M. J., Gómez-Moliner, B. J., Backeljau, T. (2016). Species limits, interspecific hybridization and phylogeny in the cryptic land snail complex Pyramidula: the power of RADseq data. Mol. Phylogenet. Evol. 101, 267–278. doi: 10.1016/j.ympev.2016.05.002
Rittmeyer, E. N., Austin, C. C. (2012). The effects of sampling on delimiting species from multi-locus sequence data. Mol. Phylogenet. Evol. 65, 451–463. doi: 10.1016/j.ympev.2012.06.031
Robinson, D. F. F., Foulds, L. R. (1981). Comparison of phylogenetic trees. Math. Biosci. 53, 131–147. doi: 10.1016/0025-5564(81)90043-2
Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D. L., Darling, A., Höhna, S., et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542. doi: 10.1093/sysbio/sys029
Shirasawa, K., Yakushiji, H., Nishimura, R., Morita, T., Jikumaru, S., Ikegami, H., et al. (2020). The Ficus erecta genome aids Ceratocystis canker resistance breeding in common fig (F.carica). Plant J. 102, 1313–1322. doi: 10.1111/tpj.14703
Spriggs, E. L., Eaton, D. A., Sweeney, P. W., Schlutius, C., Edwards, E. J., Donoghue, M. J. (2019). Restriction-site-associated DNA sequencing reveals a cryptic Viburnum species on the North American coastal plain. Syst. Biol. 68, 187–203. doi: 10.1093/sysbio/syy084
Stankowski, S., Ravinet, M. (2021). Defining the speciation continuum. Evolution 75, 1256–1273. doi: 10.1111/evo.14215
Su, Z. H., Sasaki, A., Kusumi, J., Chou, P. A., Tzeng, H. Y., Li, H. Q., et al. (2022). Pollinator sharing, copollination, and speciation by host shifting among six closely related dioecious fig species. Commun. Biol. 5, 284. doi: 10.1038/s42003-022-03223-0
Villaverde, T., Pokorny, L., Olsson, S., Rincón-Barrado, M., Johnson, M. G., Gardner, E. M., et al. (2018). Bridging the micro-and macroevolutionary levels in phylogenomics: Hyb-Seq solves relationships from populations to species and above. New Phytol. 220, 636–650. doi: 10.1111/nph.15312
Wachi, N., Kusumi, J., Tzeng, H. Y., Su, Z. H. (2016). Genome-wide sequence data suggest the possibility of pollinator sharing by host shift in dioecious figs (Moraceae, Ficus). Mol. Ecol. 25, 5732–5746. doi: 10.1111/mec.13876
Wang, Z. R. (1998). Disscussion on Concept and Chinese translation of the terms' Biosystematics' and 'Complex'. Acta Phytotaxonomica Sin. 36, 569–571.
Wang, G., Cannon, C. H., Chen, J. (2016). Pollinator sharing and gene flow among closely related sympatric dioecious fig taxa. Proc. R. Soc Ser. B Biol. Sci. 283, 20152963. doi: 10.1098/rspb.2015.2963
Wang, R., Yang, Y., Jing, Y., Segar, S. T., Zhang, Y., Wang, G., et al. (2021). Molecular mechanisms of mutualistic and antagonistic interactions in a plant–pollinator association. Nat. Ecol. Evol. 5, 974–986. doi: 10.1038/s41559-021-01469-1
Wang, G., Zhang, X., Herre, E. A., McKey, D., MaChado, C. A., Yu, W. B., et al. (2021). Genomic evidence of prevalent hybridization throughout the evolutionary history of the fig-wasp pollination mutualism. Nat. Commun. 12, 718. doi: 10.1038/s41467-021-20957-3
West-Eberhard, M. J. (1989). Phenotypic plasticity and the origins of diversity. Annu. Rev. Ecol. Systematics 20, 249–278. doi: 10.1146/annurev.es.20.110189.001341
Whittemore, A. T., Schaal, B. A. (1991). Interspecific gene flow in sympatric oaks. PNAS 88, 2540–2544. doi: 10.1073/pnas.88.6.2540
Wick, R. R., Schultz, M. B., Zobel, J., Holt, K. E. (2015). Bandage: interactive visualization of de novo genome assemblies. Bioinformatics 31, 3350–3352. doi: 10.1093/bioinformatics/btv383
Wiebes, J. T. (1994). The Indo-Australian agaoninae (pollinators of figs). Koninklijke Nederlandse Akademie van Wetenschappen Verhandelingen Afdeling Natuurkunde Tweede Reeks 92, 1–208.
Yang, L. H., Harris, A. J., Wen, F., Li, Z., Feng, C., Kong, H. H., et al. (2023). Phylogenomic analyses reveal an allopolyploid origin of core didymocarpinae (Gesneriaceae) followed by rapid radiation. Syst. Biol. 72, 1064–1083. doi: 10.1093/sysbio/syad029
Yang, Z. H., Rannala, B. (2010). Bayesian species delimitation using multilocus sequence data. PNAS 107, 9264–9269. doi: 10.1073/pnas.0913022107
Yokoyama, J. (2003). Cospeciation of figs and fig-wasps: a case study of endemic species pairs in the Ogasawara Islands. Popul. Ecol. 45, 249–256. doi: 10.1007/s10144-003-0166-4
Yokoyama, J. (2006). “Ficus,” in Flora of Japan 2a. Eds. Iwatsuki, K., Boufford, D. E., Ohba, H. (KODANSHA, Tokyo), 68–74.
Yokoyama, J., Iwatsuki, K. (1998). A faunal survey of fig-wasps (Chalcidoidea: Hymenoptera) distributed in Japan and their associations with figs (Ficus: Moraceae). Entomol. Sci. 1, 37–46.
Zhang, Z., Chen, Y. J., Mei, L., Zou, L., Zhang, M. J., Li, H. Q., et al. (2023). Phylogeny of the Ficus sarmentosa complex (Moraceae) uncovered an omitted species and the rank elevation of F. pubigera var. Maliformis. Phytotaxa 579, 17–31. doi: 10.11646/phytotaxa.579.1.2
Zhang, D., Gao, F. L., Jakovlić, I., Zou, H., Zhang, J., Li, W. X., et al. (2020). PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 20, 348–355. doi: 10.1111/1755-0998.13096
Zhang, L. F., Zhang, Z., Wang, X. M., Gao, H., Tian, H. Z., Li, H. Q. (2018). Molecular phylogeny of the Ficus auriculata complex (Moraceae). Phytotaxa 362, 39–54. doi: 10.11646/phytotaxa.362.1.3
Zhang, Z., Zhang, D. S., Zou, L., Yao, C. Y. (2022). Comparison of chloroplast genomes and phylogenomics in the Ficus sarmentosa complex (Moraceae). PloS One 17, e0279849. doi: 10.1371/journal.pone.0279849
Zhou, Z. K., Gilbert, M. G. (2003). “Moraceae,” in Flora of China. Eds. Wu, C. Y., Raven, P. H. (Science Press, Beijing), 21–73.
Zhou, W. B., Ji, X., Obata, S., Pais, A., Dong, Y. B., Peet, R., et al. (2018). Resolving relationships and phylogeographic history of the Nyssa sylvatica complex using data from RAD-seq and species distribution modeling. Mol. Phylogenet. Evol. 126, 1–16. doi: 10.1016/j.ympev.2018.04.001
Keywords: Ficus erecta complex, RADseq, plastomes, phylogeny, species delimitation
Citation: Wang X, Liao S, Zhang Z, Zhang J, Mei L and Li H (2024) Hybridization, polyploidization, and morphological convergence make dozens of taxa into one chaotic genetic pool: a phylogenomic case of the Ficus erecta species complex (Moraceae). Front. Plant Sci. 15:1354812. doi: 10.3389/fpls.2024.1354812
Received: 13 December 2023; Accepted: 08 March 2024;
Published: 26 March 2024.
Edited by:
Tingshuang Yi, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Juliano Lino Ferreira, Embrapa Pecuária Sul, BrazilJames Cohen, Weber State University, United States
Copyright © 2024 Wang, Liao, Zhang, Zhang, Mei and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongqing Li, aHFsaUBiaW8uZWNudS5lZHUuY24=