The complete chloroplast genome sequences of three Broussonetia species and comparative analysis within the Moraceae.

Background: Species of Broussonetia (family Moraceae) are commonly used to make textiles and high-grade paper. The distribution of Broussonetia papyrifera L. is considered to be related to the spread and location of humans. The complete chloroplast (cp) genomes of B. papyrifera, Broussonetia kazinoki Sieb., and Broussonetia kaempferi Sieb. were analyzed to better understand the status and evolutionary biology of the genus Broussonetia. Methods: The cp genomes were assembled and characterized using SOAPdenovo2 and DOGMA. Phylogenetic and molecular dating analysis were performed using the concatenated nucleotide sequences of 35 species in the Moraceae family and were based on 66 protein-coding genes (PCGs). An analysis of the sequence divergence (pi) of each PCG among the 35 cp genomes was conducted using DnaSP v6. Codon usage indices were calculated using the CodonW program. Results: All three cp genomes had the typical land plant quadripartite structure, ranging in size from 160,239 bp to 160,841 bp. The ribosomal protein L22 gene (RPL22) was either incomplete or missing in all three Broussonetia species. Phylogenetic analysis revealed two clades. Clade 1 included Morus and Artocarpus, whereas clade 2 included the other seven genera. Malaisia scandens Lour. was clustered within the genus Broussonetia. The differentiation of Broussonetia was estimated to have taken place 26 million years ago. The PCGs' pi values ranged from 0.0005 to 0.0419, indicating small differences within the Moraceae family. The distribution of most of the genes in the effective number of codons plot (ENc-plot) fell on or near the trend line; the slopes of the trend line of neutrality plots were within the range of 0.0363-0.171. These results will facilitate the identification, taxonomy, and utilization of the Broussonetia species and further the evolutionary studies of the Moraceae family.

A chromosome-level genome assembly of Artocarpus nanchuanensis (Moraceae), an extremely endangered fruit tree.

Artocarpus nanchuanensis (Moraceae), which is naturally distributed in China, is a representative and extremely endangered tree species. In this study, we obtained a high-quality chromosome-scale genome assembly and annotation information for A. nanchuanensis using integrated approaches, including Illumina, Nanopore sequencing platform, and Hi-C. A total of 128.71 Gb of raw Nanopore reads were generated from 20-kb libraries, and 123.38 Gb of clean reads were obtained after filtration with 160.34× coverage depth and a 17.48-kb average read length. The final assembled A. nanchuanensis genome was 769.44 Mb with a 2.09 Mb contig N50, and 99.62% (766.50 Mb) of the assembled data was assigned to 28 pseudochromosomes. In total, 39,596 genes (95.10%, 39,596/41,636) were successfully annotated, and 129 metabolic pathways were detected. Plants disease resistance/insect resistance genes, plant-pathogen interaction metabolic pathways, and abundant biosynthesis pathways of vitamins, flavonoid, and gingerol were detected. Unigene reveals the basis of species-specific functions, and gene family in contraction and expansion generally implies strong functional differences in the evolution. Compared with other related species, a total of 512 unigenes, 309 gene families in contraction, and 559 gene families in expansion were detected in A. nanchuanensis. This A. nanchuanensis genome information provides an important resource to expand our understanding of the unique biological processes, nutritional and medicinal benefits, and evolutionary relationship of this species. The study of gene function and metabolic pathway in A. nanchuanensis may reveal the theoretical basis of a special trait in A. nanchuanensis and promote the study and utilization of its rare medicinal value.

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.

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

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

Speciation slowing down in widespread and long-living tree taxa: insights from the tropical timber tree genus Milicia (Moraceae).

The long generation time and large effective size of widespread forest tree species can result in slow evolutionary rate and incomplete lineage sorting, complicating species delimitation. We addressed this issue with the African timber tree genus Milicia that comprises two morphologically similar and often confounded species: M. excelsa, widespread from West to East Africa, and M. regia, endemic to West Africa. We combined information from nuclear microsatellites (nSSRs), nuclear and plastid DNA sequences, and morphological systematics to identify significant evolutionary units and infer their evolutionary and biogeographical history. We detected five geographically coherent genetic clusters using nSSRs and three levels of genetic differentiation. First, one West African cluster matched perfectly with the morphospecies M. regia that formed a monophyletic clade at both DNA sequences. Second, a West African M. excelsa cluster formed a monophyletic group at plastid DNA and was more related to M. regia than to Central African M. excelsa, but shared many haplotypes with the latter at nuclear DNA. Third, three Central African clusters appeared little differentiated and shared most of their haplotypes. Although gene tree paraphyly could suggest a single species in Milicia following the phylogenetic species concept, the existence of mutual haplotypic exclusivity and nonadmixed genetic clusters in the contact area of the two taxa indicate strong reproductive isolation and, thus, two species following the biological species concept. Molecular dating of the first divergence events showed that speciation in Milicia is ancient (Tertiary), indicating that long-living tree taxa exhibiting genetic speciation may remain similar morphologically.

Differential effects of landscape-level environmental features on genetic structure in three codistributed tree species in Central America.

Landscape genetic studies use spatially explicit population genetic information to determine the physical and environmental causes of population genetic structure on regional scales. Comparative studies that identify common barriers to gene flow across multiple species within a community are important to both understand the evolutionary trajectories of populations and prioritize habitat conservation. Here, we use a comparative landscape genetic approach to ask whether gradients in temperature or precipitation seasonality structure genetic variation across three codistributed tree species in Central America, or whether a simpler (geographic distance) or more complex, species-specific environmental niche model is necessary to individually explain population genetic structure. Using descriptive statistics and causal modelling, we find that different factors best explain genetic distance in each of the three species: environmental niche distance in Bursera simaruba, geographic distance in Ficus insipida and historical barriers to gene flow or cryptic reproductive barriers for Brosimum alicastrum. This study confirms suggestions from previous studies of Central American tree species that imply that population genetic structure of trees in this region is determined by complex interactions of both historical and current barriers to gene flow.

Phenological patterns in a natural population of a tropical timber tree species, Milicia excelsa (Moraceae): Evidence of isolation by time and its interaction with feeding strategies of dispersers.

PREMISE OF THE STUDY: Population genetic structuring over limited timescales is commonly viewed as a consequence of spatial constraints. Indirect approaches have recently revealed reproductive isolation resulting from flowering time (so-called isolation by time, IBT). Since phenological processes can be subject to selection, the persistence of flowering asynchrony may be due to opposing selective pressures during mating, dispersal, and regeneration phases. Our study aimed to investigate phenology, fruit handling by animals, and their interaction in a timber tree species, Milicia excelsa. METHODS: We analyzed phenological data collected over 6 years on 69 genotyped trees in a Cameroonian natural rainforest complemented by data from germination trials and field observations of dispersers. KEY RESULTS: Initiation of flowering was correlated with variation in temperature and relative humidity, but was also affected by genetic factors: pairwise differences in flowering time between nearby individuals correlated with kinship coefficient, and earliness of flowering remained stable over time. A decrease in mean seed production per fruit with increasing flowering time suggests selection against late bloomers. However, germination rate was not affected by seed collection date, and the main seed disperser, the bat Eidolon helvum, seemed to increase in abundance at the end of the reproductive season and preferred trees in open habitats where early and late bloomers are expected. CONCLUSIONS: The pairwise approach performs well in detecting IBT. The persistence of different mating pools in such a case may result from a trade off between selective forces during the mating and seed dispersal processes.

Genetic variation in a tropical tree species influences the associated epiphytic plant and invertebrate communities in a complex forest ecosystem.

Genetic differences among tree species, their hybrids and within tree species are known to influence associated ecological communities and ecosystem processes in areas of limited species diversity. The extent to which this same phenomenon occurs based on genetic variation within a single tree species, in a diverse complex ecosystem such as a tropical forest, is unknown. The level of biodiversity and complexity of the ecosystem may reduce the impact of a single tree species on associated communities. We assessed the influence of within-species genetic variation in the tree Brosimum alicastrum (Moraceae) on associated epiphytic and invertebrate communities in a neotropical rainforest. We found a significant positive association between genetic distance of trees and community difference of the epiphytic plants growing on the tree, the invertebrates living among the leaf litter around the base of the tree, and the invertebrates found on the tree trunk. This means that the more genetically similar trees are host to more similar epiphyte and invertebrate communities. Our work has implications for whole ecosystem conservation management, since maintaining sufficient genetic diversity at the primary producer level will enhance species diversity of other plants and animals.

Forest refugia revisited: nSSRs and cpDNA sequences support historical isolation in a wide-spread African tree with high colonization capacity, Milicia excelsa (Moraceae).

The impact of the Pleistocene climate oscillations on the structure of biodiversity in tropical regions remains poorly understood. In this study, the forest refuge theory is examined at the molecular level in Milicia excelsa, a dioecious tree with a continuous range throughout tropical Africa. Eight nuclear microsatellites (nSSRs) and two sequences and one microsatellite from chloroplast DNA (cpDNA) showed a deep divide between samples from Benin and those from Lower Guinea. This suggests that these populations were isolated in separate geographical regions, probably for several glacial cycles of the Pleistocene, and that the nuclear gene pools were not homogenized despite M. excelsa's wind-pollination syndrome. The divide could also be related to seed dispersal patterns, which should be largely determined by the migration behaviour of M. excelsa's main seed disperser, the frugivorous bat Eidolon helvum. Within Lower Guinea, a north-south divide, observed with both marker types despite weak genetic structure (nSSRs: F(ST) = 0.035, cpDNA: G(ST) = 0.506), suggested the existence of separate Pleistocene refugia in Cameroon and the Gabon/Congo region. We inferred a pollen-to-seed dispersal distance ratio of c.1.8, consistent with wide-ranging gene dispersal by both wind and bats. Simulations in an Approximate Bayesian Computation framework suggested low nSSR and cpDNA mutation rates, but imprecise estimates of other demographic parameters, probably due to a substantial gene flow between the Lower Guinean gene pools. The decline of genetic diversity detected in some Gabonese populations could be a consequence of the relatively recent establishment of a closed canopy forest, which could negatively affect M. excelsa's reproductive system.

Secondary metabolites of Bagassa guianensis Aubl. wood: a study of the chemotaxonomy of the Moraceae family.

In order to explain the durability of the Moraceae plant family, phytochemistry of Bagassa guianensis was performed. Ethyl acetate extract was obtained from the heartwood and 18 secondary metabolites were isolated, including 6 moracins [6-O-methyl-moracin M, 6-O-methyl-moracin N and moracin Z; previously identified: moracin M, moracin N and moracin P], 8 stilbenoids [presently identified: (-)-epialboctalol and arachidin 4; previously identified: alboctalol, trans-resveratrol, arachidin 2, trans-oxyresveratrol and artogomezianol], 3 previously identified flavonoids, steppogenin, katuranin and dihydromorin, beta-sitosterol and resorcinol. Previous studies suggest that stilbenoids are responsible for the natural durability of wood. Our study has determined that B. guianensis is closely related to Morus sp. in phylogeny and should be included in the Moreae sensu stricto tribe of the Moraceae family.

Spatial genetic structure in Milicia excelsa (Moraceae) indicates extensive gene dispersal in a low-density wind-pollinated tropical tree.

In this study, we analysed spatial genetic structure (SGS) patterns and estimated dispersal distances in Milicia excelsa (Welw.) C.C. Berg (Moraceae), a threatened wind-pollinated dioecious African tree, with typically low density (approximately 10 adults/km(2)). Eight microsatellite markers were used to type 287 individuals in four Cameroonian populations characterized by different habitats and tree densities. Differentiation among populations was very low. Two populations in more open habitat did not display any correlation between genetic relatedness and spatial distance between individuals, whereas significant SGS was detected in two populations situated under continuous forest cover. SGS was weak with a maximum S(p)-statistic of 0.006, a value in the lower quartile of SGS estimates for trees in the literature. Using a stepwise approach with Bayesian clustering methods, we demonstrated that SGS resulted from isolation by distance and not colonization by different gene pools. Indirect estimates of gene dispersal distances ranged from sigma(g) = 1 to 7.1 km, one order of magnitude higher than most estimates found in the literature for tropical tree species. This result can largely be explained by life-history traits of the species. Milicia excelsa exhibits a potentially wide-ranging wind-mediated pollen dispersal mechanism as well as very efficient seed dispersal mediated by large frugivorous bats. Estimations of gene flow suggested no major risk of inbreeding because of reduction in population density by exploitation. Different strategy of seed collection may be required for reforestation programmes among populations with different extent of SGS.

Biogeography and divergence times in the mulberry family (Moraceae).

The biogeographical history of the mulberry family (Moraceae) was investigated using phylogenetic inferences from nuclear and chloroplast DNA, molecular dating with multiple fossil calibrations, and independent geological evidence. The Moraceae are centered in the tropics which has invited the hypothesis that the family has Gondwanan origins and extant distribution is the result of vicariance due to the break-up of Gondwana. However, the cosmopolitan distribution of Moraceae suggests a more complicated biogeographical history. The timing and location of Moraceae diversification also bears on the origin of the fig pollination mutualism, a model for the study of coevolution and specialization. Recent molecular dating of pollinating fig wasps suggested that an ancient Gondwanan origin coupled with vicariance and dispersal could account for the present day distribution of the mutualism. Here, we provide the first assessment of this hypothesis based on dating of figs and their relatives. Minimum age estimates suggest that the Moraceae had diversified by at least the mid-Cretaceous and major clades including the figs may have radiated during the Tertiary after the break-up of Gondwanaland. Molecular evidence together with Eurasian fossils suggest that the early diversification of Moraceae in Eurasia and subsequent migration into the southern hemisphere is at least as plausible as the Gondwanan hypothesis. These findings invite a reevaluation of the biogeography of fig pollination and highlight the need for incorporating multiple sources of evidence in biogeographical reconstructions.