First plastid phylogenomic study reveals potential cyto-nuclear discordance in the evolutionary history of Ficus L. (Moraceae)

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

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

Phylogeny of Leontopodium (Asteraceae) in China-with a reference to plastid genome and nuclear ribosomal DNA

The infrageneric taxonomy system, species delimitation, and interspecies systematic relationships of Leontopodium remain controversial and complex. However, only a few studies have focused on the molecular phylogeny of this genus. In this study, the characteristics of 43 chloroplast genomes of Leontopodium and its closely related genera were analyzed. Phylogenetic relationships were inferred based on chloroplast genomes and nuclear ribosomal DNA (nrDNA). Finally, together with the morphological characteristics, the relationships within Leontopodium were identified and discussed. The results showed that the chloroplast genomes of Filago, Gamochaeta, and Leontopodium were well-conserved in terms of gene number, gene order, and GC content. The most remarkable differences among the three genera were the length of the complete chloroplast genome, large single-copy region, small single-copy region, and inverted repeat region. In addition, the chloroplast genome structure of Leontopodium exhibited high consistency and was obviously different from that of Filago and Gamochaeta in some regions, such as matk, trnK (UUU)-rps16, petN-psbM, and trnE (UUC)-rpoB. All the phylogenetic trees indicated that Leontopodium was monophyletic. Except for the subgeneric level, our molecular phylogenetic results were inconsistent with the previous taxonomic system, which was based on morphological characteristics. Nevertheless, we found that the characteristics of the leaf base, stem types, and carpopodium base were phylogenetically correlated and may have potential value in the taxonomic study of Leontopodium. In the phylogenetic trees inferred using complete chloroplast genomes, the subgen. Leontopodium was divided into two clades (Clades 1 and 2), with most species in Clade 1 having herbaceous stems, amplexicaul, or sheathed leaves, and constricted carpopodium; most species in Clade 2 had woody stems, not amplexicaul and sheathed leaves, and not constricted carpopodium.

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

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

Origin and diversification of a Himalayan orchid genus Pleione

Pleione is an orchid endemically distributed in high mountain areas across the Hengduan Mountains (HDM), Himalayas, Southeast Asia and South of China. The unique flower shapes, rich colors and immense medicinal importance of Pleione are valuable ornamental and economic resources. However, the phylogenetic relationships and evolutionary history of the genus have not yet been comprehensively resolved. Here, the evolutionary history of Pleione was investigated using single-copy gene single nucleotide polymorphisms and chloroplast genome datasets. The data revealed that Pleione could be divided into five clades. Discordance in topology between the two phylogenetic trees and network and D-statistic analyses indicated the occurrence of reticulate evolution in the genus. The evolution could be attributed to introgression and incomplete lineage sorting. Ancestral area reconstruction suggested that Pleione was originated from the HDM. Uplifting of the HDM drove rapid diversification by creating conditions favoring rapid speciation. This coincided with two periods of consolidation of the Asian monsoon climate, which caused the first rapid diversification of Pleione from 8.87 to 7.83 Mya, and a second rapid diversification started at around 4.05 Mya to Pleistocene. The interaction between Pleione and climate changes, especially the monsoons, led to the current distribution pattern and shaped the dormancy characteristic of the different clades. In addition to revealing the evolutionary relationship of Pleione with orogeny and climate changes, the findings of this study provide insights into the speciation and diversification mechanisms of plants in the East Asian flora.

Phylogenomic evolutionary insights in the fern family Gleicheniaceae

The pantropical fern family Gleicheniaceae comprises approximately 157 species. Seven genera are currently recognized in the family, although their monophyly is still uncertain due to low sampling in phylogenetic studies. We examined the monophyly of the genera through extended sampling, using the first phylogenomic inference of the family including data from both nuclear and plastid genomes. Seventy-six samples were sequenced (70 Gleicheniaceae species and six outgroups) using high throughput sequencing, including all seven currently recognized genera. Plastid and nuclear data were recovered and assembled; the nuclear data was phased to reduce paralogy as well as hybrid noise in the final recovered topology. Maximum likelihood trees were built for each locus, and a concatenated dataset was built for both datasets. A species tree based on a multispecies coalescent model was generated, and divergence time analyses performed. We here present the first genomic phylogenetic inferences concerning Gleicheniaceae, confirming the monophyly of most genera except Sticherus, which we recovered as paraphyletic. Although most of the extant genera of Gleicheniaceae originated during the Mesozoic, several genera show Neogene and even Quaternary diversifications, and our results suggest that reticulation and polyploidy may have played significant roles during this diversification. However, some genera, such as Rouxopteris and Stromatopteris, appear to represent evolutionary relicts.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

The complete chloroplast genome sequence of Berchemia racemosa Siebold & Zucc. (Rhamnaceae), a rare plant species in Korea

Berchemia racemosa Siebold & Zucc., 1845 is a rare species distributed in restricted areas in the western Korean peninsula. In this study, the complete chloroplast genome (plastome) of B. racemosa was sequenced and assembled by Illumina paired-end sequencing. The plastome of B. racemosa was 161,187 bp in length and was quadripartite in structure, including a large single-copy (LSC) region of 89,503 bp, a small single-copy (SSC) region of 18,214 bp, and two inverted repeats of 26,735 bp. The GC content was 37.2%. The plastome of B. racemosa contains 130 genes, including eight ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and 85 protein-coding genes. Phylogenetic analysis using complete genome sequences showed that B. racemosa is most closely related to Berchemia flavescens.

Pollinator and host sharing lead to hybridization and introgression in Panamanian free-standing figs, but not in their pollinator wasps

Obligate pollination mutualisms, in which plant and pollinator lineages depend on each other for reproduction, often exhibit high levels of species specificity. However, cases in which two or more pollinator species share a single host species (host sharing), or two or more host species share a single pollinator species (pollinator sharing), are known to occur in current ecological time. Further, evidence for host switching in evolutionary time is increasingly being recognized in these systems. The degree to which departures from strict specificity differentially affect the potential for hybridization and introgression in the associated host or pollinator is unclear. We addressed this question using genome-wide sequence data from five sympatric Panamanian free-standing fig species (Ficus subgenus Pharmacosycea, section Pharmacosycea) and their six associated fig-pollinator wasp species (Tetrapus). Two of the five fig species, F. glabrata and F. maxima, were found to regularly share pollinators. In these species, ongoing hybridization was demonstrated by the detection of several first-generation (F1) hybrid individuals, and historical introgression was indicated by phylogenetic network analysis. By contrast, although two of the pollinator species regularly share hosts, all six species were genetically distinct and deeply divergent, with no evidence for either hybridization or introgression. This pattern is consistent with results from other obligate pollination mutualisms, suggesting that, in contrast to their host plants, pollinators appear to be reproductively isolated, even when different species of pollinators mate in shared hosts.

Comparison of chloroplast genomes and phylogenomics in the Ficus sarmentosa complex (Moraceae)

Due to maternal inheritance and minimal rearrangement, the chloroplast genome is an important genetic resource for evolutionary studies. However, the evolutionary dynamics and phylogenetic performance of chloroplast genomes in closely related species are poorly characterized, particularly in taxonomically complex and species-rich groups. The taxonomically unresolved Ficus sarmentosa species complex (Moraceae) comprises approximately 20 taxa with unclear genetic background. In this study, we explored the evolutionary dynamics, hotspot loci, and phylogenetic performance of thirteen chloroplast genomes (including eleven newly obtained and two downloaded from NCBI) representing the F. sarmentosa complex. Their sequence lengths, IR boundaries, repeat sequences, and codon usage were compared. Both sequence length and IR boundaries were found to be highly conserved. All four categories of long repeat sequences were found across all 13 chloroplast genomes, with palindromic and forward sequences being the most common. The number of simple sequence repeat (SSR) loci varied from 175 (F. dinganensis and F. howii) to 190 (F. polynervis), with the dinucleotide motif appearing the most frequently. Relative synonymous codon usage (RSCU) analysis indicated that codons ending with A/T were prior to those ending with C/T. The majority of coding sequence regions were found to have undergone negative selection with the exception of ten genes (accD, clpP, ndhK, rbcL, rpl20, rpl22, rpl23, rpoC1, rps15, and rps4) which exhibited potential positive selective signatures. Five hypervariable genic regions (rps15, ycf1, rpoA, ndhF, and rpl22) and five hypervariable intergenic regions (trnH-GUG-psbA, rpl32-trnL-UAG, psbZ-trnG-GCC, trnK-UUU-rps16 and ndhF-rpl32) were identified. Overall, phylogenomic analysis based on 123 Ficus chloroplast genomes showed promise for studying the evolutionary relationships in Ficus, despite cyto-nuclear discordance. Furthermore, based on the phylogenetic performance of the F. sarmentosa complex and F. auriculata complex, the chloroplast genome also exhibited a promising phylogenetic resolution in closely related species.

Plastid Phylogenomics of Paeonia and the Evolution of Ten Flower Types in Tree Peony

Paeonia suffruticosa Andr., a member of Paeoniaceae, is native to China. In its 1600 years' cultivation, more than 2000 cultivars for different purposes (ornamental, medicinal and oil use) have been inbred. However, there are still some controversies regarding the provenance of tree peony cultivars and the phylogenetic relationships between and within different cultivar groups. In this study, plastid genome sequencing was performed on 10 representative tree peony cultivars corresponding to 10 different flower types. Structure and comparative analyses of the plastid genomes showed that the total lengths of the chloroplast genome of the 10 cultivars ranged from 152,153 to 152,385 bp and encoded 84-88 protein-coding genes, 8 rRNAs and 31-40 tRNAs. The number of simple sequence repeats and interspersed repeat sequences of the 10 cultivars ranged from 65-68 and 40-42, respectively. Plastid phylogenetic relationships of Paeonia species/cultivars were reconstructed incorporating data from our newly sequenced plastid genomes and 15 published species, and results showed that subsect. Vaginatae was the closest relative to the central plains cultivar group with robust support, and that it may be involved in the formation of the group. Paeonia ostii was recovered as a successive sister group to this lineage. Additionally, eleven morphological characteristics of flowers were mapped to the phylogenetic skeleton to reconstruct the evolutionary trajectory of flower architecture in Paeoniaceae.

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

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

Compensatory Genetic and Transcriptional Cytonuclear Coordination in Allopolyploid Lager Yeast (Saccharomyces pastorianus)

Cytonuclear coordination between biparental-nuclear genomes and uniparental-cytoplasmic organellar genomes in plants is often resolved by genetic and transcriptional cytonuclear responses. Whether this mechanism also acts in allopolyploid members of other kingdoms is not clear. Additionally, cytonuclear coordination of interleaved allopolyploid cells/individuals within the same population is underexplored. The yeast Saccharomyces pastorianus provides the opportunity to explore cytonuclear coevolution during different growth stages and from novel dimensions. Using S. pastorianus cells from multiple growth stages in the same environment, we show that nuclear mitochondria-targeted genes have undergone both asymmetric gene conversion and growth stage-specific biased expression favoring genes from the mitochondrial genome donor (Saccharomyces eubayanus). Our results suggest that cytonuclear coordination in allopolyploid lager yeast species entails an orchestrated and compensatory genetic and transcriptional evolutionary regulatory shift. The common as well as unique properties of cytonuclear coordination underlying allopolyploidy between unicellular yeasts and higher plants offers novel insights into mechanisms of cytonuclear evolution associated with allopolyploid speciation.

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.

Phylogenomic analysis of 997 nuclear genes reveals the need for extensive generic re-delimitation in Caesalpinioideae (Leguminosae)

Subfamily Caesalpinioideae with ca. 4,600 species in 152 genera is the second-largest subfamily of legumes (Leguminosae) and forms an ecologically and economically important group of trees, shrubs and lianas with a pantropical distribution. Despite major advances in the last few decades towards aligning genera with clades across Caesalpinioideae, generic delimitation remains in a state of considerable flux, especially across the mimosoid clade. We test the monophyly of genera across Caesalpinioideae via phylogenomic analysis of 997 nuclear genes sequenced via targeted enrichment (Hybseq) for 420 species and 147 of the 152 genera currently recognised in the subfamily. We show that 22 genera are non-monophyletic or nested in other genera and that non-monophyly is concentrated in the mimosoid clade where ca. 25% of the 90 genera are found to be non-monophyletic. We suggest two main reasons for this pervasive generic non-monophyly: (i) extensive morphological homoplasy that we document here for a handful of important traits and, particularly, the repeated evolution of distinctive fruit types that were historically emphasised in delimiting genera and (ii) this is an artefact of the lack of pantropical taxonomic syntheses and sampling in previous phylogenies and the consequent failure to identify clades that span the Old World and New World or conversely amphi-Atlantic genera that are non-monophyletic, both of which are critical for delimiting genera across this large pantropical clade. Finally, we discuss taxon delimitation in the phylogenomic era and especially how assessing patterns of gene tree conflict can provide additional insights into generic delimitation. This new phylogenomic framework provides the foundations for a series of papers reclassifying genera that are presented here in Advances in Legume Systematics (ALS) 14 Part 1, for establishing a new higher-level phylogenetic tribal and clade-based classification of Caesalpinioideae that is the focus of ALS14 Part 2 and for downstream analyses of evolutionary diversification and biogeography of this important group of legumes which are presented elsewhere.

Comparative Analysis of the Complete Chloroplast Genomes of Eight Ficus Species and Insights into the Phylogenetic Relationships of Ficus

The genus Ficus is an evergreen plant, the most numerous species in the family Moraceae, and is often used as a food and pharmacy source. The phylogenetic relationships of the genus Ficus have been debated for many years due to the overlapping phenotypic characters and morphological similarities between the genera. In this study, the eight Ficus species (Ficus altissima, Ficus auriculata, Ficus benjamina, Ficus curtipes, Ficus heteromorpha, Ficus lyrata, Ficus microcarpa, and Ficus virens) complete chloroplast (cp) genomes were successfully sequenced and phylogenetic analyses were made with other Ficus species. The result showed that the eight Ficus cp genomes ranged from 160,333 bp (F. heteromorpha) to 160,772 bp (F. curtipes), with a typical quadripartite structure. It was found that the eight Ficus cp genomes had similar genome structures, containing 127 unique genes. The cp genomes of the eight Ficus species contained 89−104 SSR loci, which were dominated by mono-nucleotides repeats. Moreover, we identified eight hypervariable regions (trnS-GCU_trnG-UCC, trnT-GGU_psbD, trnV-UAC_trnM-CAU, clpP_psbB, ndhF_trnL-UAG, trnL-UAG_ccsA, ndhD_psaC, and ycf1). Phylogenetic analyses have shown that the subgenus Ficus and subgenus Synoecia exhibit close affinities and based on the results, we prefer to merge the subgenus Synoecia into the subgenus Ficus. At the same time, new insights into the subgeneric classification of the Ficus macrophylla were provided. Overall, these results provide useful data for further studies on the molecular identification, phylogeny, species identification and population genetics of speciation in the Ficus genus.

Comprehensive analysis of complete chloroplast genome and phylogenetic aspects of ten Ficus species

BACKGROUND

The large genus Ficus comprises approximately 800 species, most of which possess high ornamental and ecological values. However, its evolutionary history remains largely unknown. Plastome (chloroplast genome) analysis had become an essential tool for species identification and for unveiling evolutionary relationships between species, genus and other rank groups. In this work we present the plastomes of ten Ficus species.

RESULTS

The complete chloroplast (CP) genomes of eleven Ficus specimens belonging to ten species were determined and analysed. The full length of the Ficus plastome was nearly 160 kbp with a similar overall GC content, ranging from 35.88 to 36.02%. A total of 114 unique genes, distributed in 80 protein-coding genes, 30 tRNAs, and 4 rRNAs, were annotated in each of the Ficus CP genome. In addition, these CP genomes showed variation in their inverted repeat regions (IR). Tandem repeats and mononucleotide simple sequence repeat (SSR) are widely distributed across the Ficus CP genome. Comparative genome analysis showed low sequence variability. In addition, eight variable regions to be used as potential molecular markers were proposed for future Ficus species identification. According to the phylogenetic analysis, these ten Ficus species were clustered together and further divided into three clades based on different subgenera. Simultaneously, it also showed the relatedness between Ficus and Morus.

CONCLUSION

The chloroplast genome structure of 10 Ficus species was similar to that of other angiosperms, with a typical four-part structure. Chloroplast genome sizes vary slightly due to expansion and contraction of the IR region. And the variation of noncoding regions of the chloroplast genome is larger than that of coding regions. Phylogenetic analysis showed that these eleven sampled CP genomes were divided into three clades, clustered with species from subgenus Urostigma, Sycomorus, and Ficus, respectively. These results support the Berg classification system, in which the subgenus Ficus was further decomposed into the subgenus Sycomorus. In general, the sequencing and analysis of Ficus plastomes, especially the ones of species with no or limited sequences available yet, contribute to the study of genetic diversity and species evolution of Ficus, while providing useful information for taxonomic and phylogenetic studies of Ficus.

Dark-colored Oncocyclus irises in Israel analyzed by AFLP, whole chloroplast genome sequencing and species distribution modeling

The Haynei is one of seven species aggregates (clusters of species having similar flower morphology) recognized in section Oncocyclus of genus Iris. This aggregate, characterized by dark-colored flowers, is represented by six species in Israel and adjacent Jordan. There is, however, no knowledge of the genetic relationship of these species making verification of their taxonomic status impossible. We investigated genetic variation in this group using analysis of whole chloroplast genomes and amplified fragment length polymorphism (AFLP). We also used species distribution modeling (SDM) to predict species ranges under current climatic conditions. We found some population groups within the currently recognized species of section Oncocyclus to represent dramatically different genetic entities which devaluates a general trend of merging many previously recognized species of section Oncocyclus based on their flower morphology. Despite the importance of homoploid hybridization in this group’s evolution and some apparently sporadically happening inter-specific gene flow, the main evolutionary forces in Oncocyclus appear to be vicariance and spatial isolation. Our findings suggest that some of the currently recognized species in section Oncocyclus need revision. A revision must be based on genetic analyses allowing the reconstruction of ancestry and recognition of the importance of vicariance and spatial isolation in the evolution of this group. The implications of the present findings for conservation are discussed.

Genome-wide sequence data show no evidence of hybridization and introgression among pollinator wasps associated with a community of Panamanian strangler figs

The specificity of pollinator host choice influences opportunities for reproductive isolation in their host plants. Similarly, host plants can influence opportunities for reproductive isolation in their pollinators. For example, in the fig and fig wasp mutualism, offspring of fig pollinator wasps mate inside the inflorescence that the mothers pollinate. Although often host specific, multiple fig pollinator species are sometimes associated with the same fig species, potentially enabling hybridization between wasp species. Here, we study the 19 pollinator species (Pegoscapus spp.) associated with an entire community of 16 Panamanian strangler fig species (Ficus subgenus Urostigma, section Americanae) to determine whether the previously documented history of pollinator host switching and current host sharing predicts genetic admixture among the pollinator species, as has been observed in their host figs. Specifically, we use genome‐wide ultraconserved element (UCE) loci to estimate phylogenetic relationships and test for hybridization and introgression among the pollinator species. In all cases, we recover well‐delimited pollinator species that contain high interspecific divergence. Even among pairs of pollinator species that currently reproduce within syconia of shared host fig species, we found no evidence of hybridization or introgression. This is in contrast to their host figs, where hybridization and introgression have been detected within this community, and more generally, within figs worldwide. Consistent with general patterns recovered among other obligate pollination mutualisms (e.g. yucca moths and yuccas), our results suggest that while hybridization and introgression are processes operating within the host plants, these processes are relatively unimportant within their associated insect pollinators.

Synthesis of Nuclear and Chloroplast Data Combined With Network Analyses Supports the Polyploid Origin of the Apple Tribe and the Hybrid Origin of the Maleae-Gillenieae Clade

Plant biologists have debated the evolutionary origin of the apple tribe (Maleae; Rosaceae) for over a century. The "wide-hybridization hypothesis" posits that the pome-bearing members of Maleae (base chromosome number x = 17) resulted from a hybridization and/or allopolyploid event between progenitors of other tribes in the subfamily Amygdaloideae with x = 8 and x = 9, respectively. An alternative "spiraeoid hypothesis" proposed that the x = 17 of Maleae arose via the genome doubling of x = 9 ancestors to x = 18, and subsequent aneuploidy resulting in x = 17. We use publicly available genomic data-448 nuclear genes and complete plastomes-from 27 species representing all major tribes within the Amygdaloideae to investigate evolutionary relationships within the subfamily containing the apple tribe. Specifically, we use network analyses and multi-labeled trees to test the competing wide-hybridization and spiraeoid hypotheses. Hybridization occurred between an ancestor of the tribe Spiraeeae (x = 9) and an ancestor of the clade Sorbarieae (x = 9) + Exochordeae (x = 8) + Kerrieae (x = 9), giving rise to the clade Gillenieae (x = 9) + Maleae (x = 17). The ancestor of the Maleae + Gillenieae arose via hybridization between distantly related tribes in the Amygdaloideae (i.e., supporting the wide hybridization hypothesis). However, some evidence supports an aspect of the spiraeoid hypothesis-the ancestors involved in the hybridization event were likely both x = 9, so genome doubling was followed by aneuploidy to result in x = 17 observed in Maleae. By synthesizing existing genomic data with novel analyses, we resolve the nearly century-old mystery regarding the origin of the apple tribe. Our results also indicate that nuclear gene tree-species tree conflict and/or cytonuclear conflict are pervasive at several other nodes in subfamily Amygdaloideae of Rosaceae.

Exploring systematic biases, rooting methods and morphological evidence to unravel the evolutionary history of the genus Ficus (Moraceae)

Despite many attempts in the Sanger sequencing era, the phylogeny of fig trees remains unresolved, which limits our ability to analyze the evolution of key traits that may have contributed to their evolutionary and ecological success. We used restriction-site-associated DNA sequencing (c. 420 kb) and 102 morphological characters to elucidate the relationships between 70 species of Ficus. To increase phylogenetic information for higher-level relationships, we targeted conserved regions and assembled paired reads into long loci to enable the retrieval of homologous loci in outgroup genomes. We compared morphological and molecular results to highlight discrepancies and reveal possible inference bias. For the first time, we recovered a monophyletic subgenus Urostigma (stranglers) and a clade with all gynodioecious Ficus. However, we show, with a new approach based on iterative principal component analysis, that it is not (and will probably never be) possible to homogenize evolutionary rates and GC content for all taxa before phylogenetic inference. Four competing positions for the root of the molecular tree are possible. The placement of section Pharmacosycea as sister to other fig trees is not supported by morphological data and considered a result of a long-branch attraction artefact to the outgroups. Regarding morphological features and indirect evidence from the pollinator tree of life, the topology that divides Ficus into monoecious versus gynodioecious species appears most plausible. It seems most likely that the ancestor of fig trees was a freestanding tree and active pollination is inferred as the ancestral state, contrary to previous hypotheses. However, ambiguity remains on the ancestral breeding system. Despite morphological plasticity, we advocate restoring a central role to morphology in our understanding of the evolution of Ficus, as it can help detect systematic errors that appear more pronounced with larger molecular datasets.

Phylogeny, origin, and dispersal of Dubyaea (Asteraceae) based on Hyb-Seq data

Dubyaea DC. is a small genus of Asteraceae that is almost exclusively endemic to the Pan-Himalayan region. Within Dubyaea, phylogenetic relationships remain poorly understood. Here, our well-supported phylogeny based on Hyb-Seq data shows that all samples of Dubyaea in this study belong to a monophyletic group, which is sister to the clade of Soroseris, Syncalathium, and Nabalus. Dubyaea (s. str.) can be divided into three major clades, which are supported by flower color as well as morphological features of the stems and basal leaves. Based on our phylogenetic results, we performed biogeographic analyses and inferred that Dubyaea arose in the late Miocene in Hengduan Mountains and its eastern areas. Following its evolutionary origin, Dubyaea underwent diversification in situ as well as spread to the Himalayas.

Phylogenetics and comparative plastome genomics of two of the largest genera of angiosperms, Piper and Peperomia (Piperaceae)

Biological radiations provide unique opportunities to understand the evolution of biodiversity. One such radiation is the pepper plant family Piperaceae, an early-diverging and mega-diverse lineage that could serve as a model to study the diversification of angiosperms. However, traditional genetic markers lack sufficient variation for such studies, and testing hypotheses on poorly resolved phylogenetic frameworks becomes challenging. Limited genomic data is available for Piperaceae, which contains two of the largest genera of angiosperms, Piper (>2100 species) and Peperomia (>1300 species). To address this gap, we used genome skimming to assemble and annotate whole plastomes (152-161kbp) and >5kbp nuclear ribosomal DNA region from representatives of Piper and Peperomia. We conducted phylogenetic and comparative genomic analyses to study plastome evolution and investigate the role of hybridization in this group. Plastome phylogenetic trees were well resolved and highly supported, with a hard incongruence observed between plastome and nuclear phylogenetic trees suggesting hybridization in Piper. While all plastomes of Piper and Peperomia had the same gene content and order, there were informative structural differences between them. First, ycf1 was more variable and longer in Piper than Peperomia, extending well into the small single copy region by thousands of base pairs. We also discovered previously unknown structural variation in 14 out of 25 Piper taxa, tandem duplication of the trnH-GUG gene resulting in an expanded large single copy region. Other early-diverging angiosperms have a duplicated trnH-GUG, but the specific rearrangement we found is unique to Piper and serves to refine knowledge of relationships among early-diverging angiosperms. Our study demonstrates that genome skimming is an efficient approach to produce plastome assemblies for comparative genomics and robust phylogenies of species-rich plant genera.

The plastid genome of a narrowly distributed species Artocarpus petelotii (Moraceae)

The complete chloroplast genome sequence of Artocarpus petelotii was determined, a narrowly distributed species at high altitudes in the family Moraceae. To better determine its phylogenetic location with respect to the other Moraceae species, the complete plastid genome of A. petelotii was sequenced. The whole chloroplast genome is 161,009 bp in length, consisting of a pair of inverted repeat (IR) regions of 25,682 bp, one large single-copy (LSC) region of 89,552 bp, and one small single-copy (SSC) region of 20,093 bp. The overall GC content of the whole chloroplast genome is 35.8%. Further, maximum likelihood phylogenetic analysis was conducted using 26 complete plastomes of the Moraceae, which support close relationships among A. petelotii, A. nanchuanensis and A. heterophyllus.

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

Ficus (figs) and their agaonid wasp pollinators present an ecologically important mutualism that also provides a rich comparative system for studying functional co-diversification throughout its coevolutionary history (~75 million years). We obtained entire nuclear, mitochondrial, and chloroplast genomes for 15 species representing all major clades of Ficus. Multiple analyses of these genomic data suggest that hybridization events have occurred throughout Ficus evolutionary history. Furthermore, cophylogenetic reconciliation analyses detect significant incongruence among all nuclear, chloroplast, and mitochondrial-based phylogenies, none of which correspond with any published phylogenies of the associated pollinator wasps. These findings are most consistent with frequent host-switching by the pollinators, leading to fig hybridization, even between distantly related clades. Here, we suggest that these pollinator host-switches and fig hybridization events are a dominant feature of fig/wasp coevolutionary history, and by generating novel genomic combinations in the figs have likely contributed to the remarkable diversity exhibited by this mutualism.

Horizontal genome transfer by cell-to-cell travel of whole organelles

Recent work has revealed that both plants and animals transfer genomes between cells. In plants, horizontal transfer of entire plastid, mitochondrial, or nuclear genomes between species generates new combinations of nuclear and organellar genomes, or produces novel species that are allopolyploid. The mechanisms of genome transfer between cells are unknown. Here, we used grafting to identify the mechanisms involved in plastid genome transfer from plant to plant. We show that during proliferation of wound-induced callus, plastids dedifferentiate into small, highly motile, amoeboid organelles. Simultaneously, new intercellular connections emerge by localized cell wall disintegration, forming connective pores through which amoeboid plastids move into neighboring cells. Our work uncovers a pathway of organelle movement from cell to cell and provides a mechanistic framework for horizontal genome transfer.

The Complete Plastid Genome of Artocarpus camansi: A High Degree of Conservation of the Plastome Structure in the Family Moraceae

Understanding the plastid genome is extremely important for the interpretation of the genetic mechanisms associated with essential physiological and metabolic functions, the identification of possible marker regions for phylogenetic or phylogeographic analyses, and the elucidation of the modes through which natural selection operates in different regions of this genome. In the present study, we assembled the plastid genome of Artocarpus camansi, compared its repetitive structures with Artocarpus heterophyllus, and searched for evidence of synteny within the family Moraceae. We also constructed a phylogeny based on 56 chloroplast genes to assess the relationships among three families of the order Rosales, that is, the Moraceae, Rhamnaceae, and Cannabaceae. The plastid genome of A. camansi has 160,096 bp, and presents the typical circular quadripartite structure of the Angiosperms, comprising a large single copy (LSC) of 88,745 bp and a small single copy (SSC) of 19,883 bp, separated by a pair of inverted repeat (IR) regions each with a length of 25,734 bp. The total GC content was 36.0%, which is very similar to Artocarpus heterophyllus (36.1%) and other moraceous species. A total of 23,068 codons and 80 SSRs were identified in the A. camansi plastid genome, with the majority of the SSRs being mononucleotide (70.0%). A total of 50 repeat structures were observed in the A. camansi plastid genome, in contrast with 61 repeats in A. heterophyllus. A purifying selection signal was found in 70 of the 79 protein-coding genes, indicating that they have all been highly conserved throughout the evolutionary history of the genus. The comparative analysis of the structural characteristics of the chloroplast among different moraceous species found a high degree of similarity in the sequences, which indicates a highly conserved evolutionary model in these plastid genomes. The phylogenetic analysis also recovered a high degree of similarity between the chloroplast genes of A. camansi and A. heterophyllus, and reconfirmed the hypothesis of the intense conservation of the plastome in the family Moraceae.

Phylogenomics of the Andean Tetraploid Clade of the American Amaryllidaceae (Subfamily Amaryllidoideae): Unlocking a Polyploid Generic Radiation Abetted by Continental Geodynamics

One of the two major clades of the endemic American Amaryllidaceae subfam. Amaryllidoideae constitutes the tetraploid-derived (n = 23) Andean-centered tribes, most of which have 46 chromosomes. Despite progress in resolving phylogenetic relationships of the group with plastid and nrDNA, certain subclades were poorly resolved or weakly supported in those previous studies. Sequence capture using anchored hybrid enrichment was employed across 95 species of the clade along with five outgroups and generated sequences of 524 nuclear genes and a partial plastome. Maximum likelihood phylogenetic analyses were conducted on concatenated supermatrices, and coalescent-based species tree analyses were run on the gene trees, followed by hybridization network, age diversification and biogeographic analyses. The four tribes Clinantheae, Eucharideae, Eustephieae, and Hymenocallideae (sister to Clinantheae) are resolved in all analyses with > 90 and mostly 100% support, as are almost all genera within them. Nuclear gene supermatrix and species tree results were largely in concordance; however, some instances of cytonuclear discordance were evident. Hybridization network analysis identified significant reticulation in Clinanthus, Hymenocallis, Stenomesson and the subclade of Eucharideae comprising Eucharis, Caliphruria, and Urceolina. Our data support a previous treatment of the latter as a single genus, Urceolina, with the addition of Eucrosia dodsonii. Biogeographic analysis and penalized likelihood age estimation suggests an origin in the Cauca, Desert and Puna Neotropical bioprovinces for the complex in the mid-Oligocene, with more dispersals than vicariances in its history, but no extinctions. Hymenocallis represents the only instance of long-distance vicariance from the tropical Andean origin of its tribe Hymenocallideae. The absence of extinctions correlates with the lack of diversification rate shifts within the clade. The Eucharideae experienced a sudden lineage radiation ca. 10 Mya. We tie much of the divergences in the Andean-centered lineages to the rise of the Andes, and suggest that the Amotape-Huancabamba Zone functioned as both a corridor (dispersal) and a barrier to migration (vicariance). Several taxonomic changes are made. This is the largest DNA sequence data set to be applied within Amaryllidaceae to date.

The plastid genome of winter cropping plants Ficus tinctoria (Moraceae)

Ficus tinctoria subsp. gibbosa (Blume) Corner is a hemiepiphytic dioecious fig species of the genus Ficus in the family Moraceae. To better determine its phylogenetic location with respect to the other Ficus species, the complete plastid genome of F. tinctoria was sequenced. The whole plastome is 160,342 bp in length, consisting of a pair of inverted repeat (IR) regions of 25,859 bp, one large single-copy (LSC) region of 88,526 bp, and one small single-copy (SSC) region of 20,098 bp. The overall GC content of the whole plastome is 35.9%. Further, maximum likelihood (ML) phylogenetic analyze was conducted using 23 complete fig plastomes, which support close relationships among F. tinctoria, F. heteropleura, F. obscura, and F. deltoidea.

The plastid genome of a large hemiepiphytic plants Ficus altissima (Moraceae)

Ficus altissima Blume is a hemiepiphytic monoecious fig species of the genus Ficus in the family Moraceae. To better determine its phylogenetic location with respect to the other Ficus species, the complete plastid genome of F. altissima was sequenced. The whole plastome is 160,251 bp in length, consisting of a pair of inverted repeat (IR) regions of 25,886 bp, one large single-copy (LSC) region of 88,470 bp, and one small single-copy (SSC) region of 20,009 bp. The overall GC content of the whole plastome is 35.9%. Further, maximum likelihood phylogenetic analyses was conducted using 29 complete fig plastomes, which support close relationships among F. altissima, F. benjamina, F. microcarpa, and F. consociata.

Complete plastome sequencing resolves taxonomic relationships among species of Calligonum L. (Polygonaceae) in China

BACKGROUND

Calligonum (Polygonaceae) is distributed from southern Europe through northern Africa to central Asia, and is typically found in arid, desert regions. Previous studies have revealed that standard DNA barcodes fail to discriminate Calligonum species. In this study, the complete plastid genomes (plastome) for 32 accessions of 21 Calligonum species is sequenced to not only generate the first complete plastome sequence for the genus Calligonum but to also 1) Assess the ability of the complete plastome sequence to discern species within the group, and 2) screen the plastome sequence for a cost-effective DNA barcode that can be used in future studies to resolve taxonomic relationships within the group.

RESULTS

The whole plastomes of Calligonum species possess a typical quadripartite structure. The size of the Calligonum plastome is approximately 161 kilobase pairs (kbp), and encodes 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Based on ML phylogenetic tree analyses, the complete plastome has higher species identification (78%) than combinations of standard DNA barcodes (rbcL + matK + nrITS, 56%). Five newly screened gene regions (ndhF, trnS-G, trnC-petN, ndhF-rpl32, rpl32-trnL) had high species resolution, where ndhF and trnS-G were able to distinguish the highest proportion of Calligonum species (56%).

CONCLUSIONS

The entire plastid genome was the most effective barcode for the genus Calligonum, although other gene regions showed great potential as taxon-specific barcodes for species identification in Calligonum.

Pharmaceutical resource discovery from traditional medicinal plants: Pharmacophylogeny and pharmacophylogenomics

The worldwide botanical and medicinal culture diversity are astonishing and constitute a Pierian spring for innovative drug R&D. Here, the latest awareness and the perspectives of pharmacophylogeny and pharmacophylogenomics, as well as their expanding utility in botanical drug R&D, are systematically summarized and highlighted. Chemotaxonomy is based on the fact that closely related plants contain the same or similar chemical profiles. Correspondingly, it is better to combine morphological characters, DNA markers and chemical markers in the inference of medicinal plant phylogeny. Medicinal plants within the same phylogenetic groups may have the same or similar therapeutic effects, thus forming the core of pharmacophylogeny. Here we systematically review and comment on the versatile applications of pharmacophylogeny in (1) looking for domestic resources of imported drugs, (2) expanding medicinal plant resources, (3) quality control, identification and expansion of herbal medicines, (4) predicting the chemical constituents or active ingredients of herbal medicine and assisting in the identification and determination of chemical constituents, (5) the search for new drugs sorting out, and (6) summarizing and improving herbal medicine experiences, etc. Such studies should be enhanced within the context of deeper investigations of molecular biology and genomics of traditional medicinal plants, phytometabolites and metabolomics, and ethnomedicine-based pharmacological activity, thus enabling the sustainable conservation and utilization of traditional medicinal resources.

Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae)

BACKGROUND

Paris (Melanthiaceae) is an economically important but taxonomically difficult genus, which is unique in angiosperms because some species have extremely large nuclear genomes. Phylogenetic relationships within Paris have long been controversial. Based on complete plastomes and nuclear ribosomal DNA (nrDNA) sequences, this study aims to reconstruct a robust phylogenetic tree and explore historical biogeography and clade diversification in the genus.

RESULTS

All 29 species currently recognized in Paris were sampled. Whole plastomes and nrDNA sequences were generated by the genome skimming approach. Phylogenetic relationships were reconstructed using the maximum likelihood and Bayesian inference methods. Based on the phylogenetic framework and molecular dating, biogeographic scenarios and historical diversification of Paris were explored. Significant conflicts between plastid and nuclear datasets were identified, and the plastome tree is highly congruent with past interpretations of the morphology. Ancestral area reconstruction indicated that Paris may have originated in northeastern Asia and northern China, and has experienced multiple dispersal and vicariance events during its diversification. The rate of clade diversification has sharply accelerated since the Miocene/Pliocene boundary.

CONCLUSIONS

Our results provide important insights for clarifying some of the long-standing taxonomic debates in Paris. Cytonuclear discordance may have been caused by ancient and recent hybridizations in the genus. The climatic and geological changes since the late Miocene, such as the intensification of Asian monsoon and the rapid uplift of Qinghai-Tibet Plateau, as well as the climatic fluctuations during the Pleistocene, played essential roles in driving range expansion and radiative diversification in Paris. Our findings challenge the theoretical prediction that large genome sizes may limit speciation.

Estimating divergence times and ancestral breeding systems in Ficus and Moraceae

Background and Aims

Although dioecy, which characterizes only 6 % of angiosperm species, has been considered an evolutionary dead end, recent studies have demonstrated that this is not necessarily the case. Moraceae (40 genera, 1100 spp., including Ficus, 750 spp.) are particularly diverse in breeding systems (including monoecy, gynodioecy, androdioecy and dioecy) and thus represent a model clade to study macroevolution of dioecy.

Methods

Ancestral breeding systems of Ficus and Moraceae were inferred. To do so, a new dated phylogenetic tree of Ficus and Moraceae was first reconstructed by combining a revised 12-fossil calibration set and a densely sampled molecular data set of eight markers and 320 species. Breeding system evolution was then reconstructed using both parsimony and model-based (maximum likelihood and Bayesian) approaches with this new time scale.

Key Results

The crown group ages of Ficus and Moraceae were estimated in the Eocene (40.6-55.9 Ma) and Late Cretaceous (73.2-84.7 Ma), respectively. Strong support was found for ancestral dioecy in Moraceae. Although the ancestral state of Ficus remained particularly sensitive to model selection, the results show that monoecy and gynodioecy evolved from dioecy in Moraceae, and suggest that gynodioecy probably evolved from monoecy in Ficus.

Conclusions

Dioecy was found not to be an evolutionary dead end in Moraceae. This study provides a new time scale for the phylogeny and a new framework of breeding system evolution in Ficus and Moraceae.

Evolution of pollination by frugivorous birds in Neotropical Myrtaceae

Bird pollination is relatively common in the tropics, and especially in the Americas. In the predominantly Neotropical tribe Myrteae (Myrtaceae), species of two genera, Acca and Myrrhinium, offer fleshy, sugary petals to the consumption of birds that otherwise eat fruits, thus pollinating the plants in an unusual plant-animal interaction. The phylogenetic position of these genera has been problematic, and therefore, so was the understanding of the evolution of this interaction. Here we include new sequences of Myrrhinium atropurpureum in a comprehensive molecular phylogeny based on a balanced sample of two plastid and two nuclear markers, with the aim of providing the historical framework of pollination by frugivorous birds in Myrteae. We developed 13 flower and inflorescence characters that comprehensively depict the macroscopic morphological components of this interaction. Bayesian and parsimony phylogenies concur in placing both Acca and Myrrhinium in a clade with Psidium species; with Myrrhinium sister to Psidium. Mapping of morphological characters indicated some degree of convergence (e.g., fleshy petals, purplish display) but also considerable divergence in key characters that point to rather opposing pollination strategies and also different degrees of specialization in Acca versus Myrrhinium. Pollination by frugivorous birds represents a special case of mutualism that highlights the evolutionary complexities of plant-animal interactions.

Atlantic forests to the all Americas: Biogeographical history and divergence times of Neotropical Ficus (Moraceae)

Ficus (Moraceae) is well diversified in the Neotropics with two lineages inhabiting the wet forests of this region. The hemiepiphytes of section Americanae are the most diversified with c. 120 species, whereas section Pharmacosycea includes about 20 species mostly with a terrestrial habit. To reconstruct the biogeographical history and diversification of Ficus in the Americas, we produced a dated Bayesian phylogenetic hypothesis of Neotropical Ficus including two thirds of the species sequenced for five nuclear regions (At103, ETS, G3pdh, ITS/5.8S and Tpi). Ancestral range was estimated using all models available in Biogeobears and Binary State Speciation and Extinction analysis was used to evaluate the role of the initial habit and propagule size in diversification. The phylogenetic analyses resolved both Neotropical sections as monophyletic but the internal relationships between species in section Americanae remain unclear. Ficus started their diversification in the Neotropics between the Oligocene and Miocene. The genus experienced two bursts of diversification: in the middle Miocene and the Pliocene. Colonization events from the Amazon to adjacent areas coincide with the end of the Pebas system (10 Mya) and the connection of landmasses. Divergence of endemic species in the Atlantic forest is inferred to have happened after its isolation and the opening and consolidation of the Cerrado. Our results suggest a complex diversification in the Atlantic forest differing between postulated refuges and more instable areas in the South distribution of the forest. Finally the selection for initial hemiepiphytic habit and small to medium propagule size influenced the diversification and current distribution of the species at Neotropical forests marked by the historical instability and long-distance dispersal.