Rosane G. Collevatti
Jacqueline S. Lima
Liliana Ballesteros-Mejia- Laboratório de Genética and Biodiversidade, ICB, Universidade Federal de Goiás (UFG), Goiânia, Brazil
Restricted gene flow may lead to the loss of genetic diversity and higher genetic differentiation among populations, but the genetic consequences of megafauna extinction for plant populations still remain to be assessed. We performed a phylogenetic-independent meta-analysis across 102 Neotropical plants to test the hypothesis that plant species with megafaunal seed dispersal syndrome have a lower genetic diversity and a higher genetic differentiation than those without it. We classified as megafauna-dependent plant species those that potentially relied only on megafauna to seed dispersal, and as megafauna-independent those that relied on megafauna and other seed dispersers. Our data comprised 98 studies using microsatellite markers. We found no statistical difference in genetic diversity and differentiation between plants with megafauna and non-megafauna seed dispersal syndrome, although the statistical power to detect differences in genetic differentiation was low. Moreover, we found no statistical difference between megafauna-dependent and megafauna-independent plant species. We then used generalized linear mixed models and phylogenetic generalized least square models to investigate the effects of megafaunal seed dispersal syndromes and reproductive traits on variation in genetic diversity and genetic differentiation. We found no effect of megafaunal syndrome, rather, reproductive traits, such as pollination mode, mating, and breeding systems, showed significant effects. Our findings show that the genetic studies of Neotropical plants performed so far show no difference in genetic diversity and differentiation in plants with megafaunal compared to those with non-megafaunal seed dispersal syndromes. Our results also provide evidence pointing out that plant species with megafaunal seed dispersal syndromes may have used different strategies to counterbalance the extinction of their mutualistic megafauna dispersers, such as the dispersal by extant mammals that may promote long-distance seed dispersal. Our results also reinforce the importance of pollination to long-distance gene flow in Neotropical plants.
Introduction
Neotropical plants producing large fruits and seeds or seeds with thick and hard endocarps or coats embedded in indehiscent fruits are known as biological anachronisms (Janzen and Martin, 1982), adapted to a now extinct megafauna (mammals >1,000 kg) that once dominated Pleistocene Neotropical landscapes. Nowadays, seeds of several Neotropical species matching Janzen and Martin’s (1982) megafauna syndrome remain below their mother trees without being dispersed. It raises the hypothesis that these plants have lost their seed dispersers, especially the long-distance seed-dispersing species. Megafauna species may have promoted long-distance dispersal by massive seed dispersal (Pires et al., 2017), that is, defecating large amounts of seeds over large areas, or by carrying large individual seeds over long distances (Janzen and Martin, 1982; Guimarães et al., 2008; Hansen and Galetti, 2009; Pires et al., 2017).
In the Neotropics, fruits of megafauna syndrome species are consumed by many extant mammals, such as agouti, but only extant large mammals, such as tapirs, large bats, Atelinae monkeys, manned wolf, and deer, can act as long-distance seed dispersers by endozoochory (Tabarelli and Peres, 2002; Guimarães et al., 2008; Hansen and Galetti, 2009). In addition, germination rates of seeds swallowed and defecated by larger mammals may be higher than those of seeds consumed by smaller mammals because of greater scarification (Westoby et al., 1996; Moles and Westoby, 2004; Bodmer and Ward, 2006).
Even though extant frugivores may disperse seeds with megafaunal syndrome, they might not counterbalance the effects of the extinct megafauna species following the biological anachronism reasoning (Guimarães et al., 2008; Hansen and Galetti, 2009) because of both smaller body sizes and gut retention time (Van Soest, 1996; Clauss et al., 2007; Clauss et al., 2008). Predictions of the effects of megafauna extinction include a decrease in the number of seeds successfully dispersed away from the maternal plant, a decrease in recruitment because of high mortality of nondispersed seeds close to the mother plant, and restricted gene flow because of the loss of long-distance dispersal agents (Janzen and Martin, 1982; Collevatti et al., 2003; Guimarães et al., 2008; Hansen and Galetti, 2009). Nevertheless, plant species linked to megafauna syndromes may have counterbalanced the extinction of their long-distance seed dispersers by relying on scatter-hoarding rodents (Jansen et al., 2012).
Megafauna herbivores may have also affected community structure because of their foraging behavior, causing disturbance in vegetation community, creating gaps, and reducing vegetation density, thus increasing spatial heterogeneity (Johnson, 2009). Because of high deposition of urine and feces, they might also have increased nutrient recycling, affecting community structure, and reduced fire frequency because of grazing (e.g., Hester et al., 2006; Gill et al., 2009; Gill, 2014; Malhi et al., 2016). Although some of the consequences of the extinction of seed dispersers have been assessed for community structure (Hester et al., 2006; Gill et al., 2009), nutrient cycling (Malhi et al., 2016) and plant-frugivore interactions (Gill, 2014), the genetic consequences for plant populations still remain to be evaluated. Many authors have predicted that the loss of megafauna dispersal may lead to a lower gene flow among populations, leading in turn to the loss of genetic diversity and higher genetic differentiation among populations (e.g., Collevatti et al., 2003; Guimarães et al., 2008; Pires et al., 2014; Malhi et al., 2016); however, those predictions have not been tested so far. Alternatively, the effects of reproductive traits (RT), such as pollination, mating, or breeding mode, have not been tested either. In the Tropics, plant species with outcrossing mating systems show mixed results regarding genetic differentiation (Bawa, 1992). Pollination modes and mating systems have been related to patterns of genetic diversity and differentiation in Neotropical plants (Ballesteros-Mejia et al., 2016). For instance, Neotropical plants with outcrossing mating systems have a lower genetic differentiation (FST), and beetle-pollinated species have a higher allelic richness (AR) (Ballesteros-Mejia et al., 2016). Thus, RT may account for differences in genetic diversity and differentiation among Neotropical plants, despite seed dispersal syndrome (Ballesteros-Mejia et al., 2016).
Herein, we address whether Neotropical plant species linked to megafauna syndromes have indeed a lower genetic diversity and a higher genetic differentiation among populations using a phylogenetic meta-analysis framework. In addition, we test whether megafauna seed dispersal syndrome or other RT, such as breeding system, mating system, and pollination mode, are responsible for variations in genetic diversity and differentiation among species. We model the variation in genetic parameters using generalized linear mixed models (GLMM) and phylogenetic generalized least squares (pGLS) to account for phylogenetic non-independence. Our analyses are based on published data on the polymorphism at nuclear microsatellite molecular markers.
Materials and Methods
Data Selection
We compiled our database by performing an exhaustive search for published studies on population genetics of Neotropical plant species. We used the online databases Web of Science platform (ISI, www.webofknowledge.com), Scopus (http://www.elsevier.com/online-tools/scopus), and Portal de Periódicos Capes (http://www.periodicos.capes.gov.br). The online search engine Google (www.scholar.google.com) was also used to identify scientific reports, theses, and other gray literature. Our data compilation includes publications dating from 1945 (first register in ISI) to December 2016. We carried out the search with the most common keywords in the area: “phylogeography,” “population genetics,” “genetic diversity,” and “genetic structure,” combined with (and) “Neotropics,” “Neotropical tree*,” “Neotropical plant*.” We considered the Neotropics as the region comprising the Neotropical floristic region (Cox and Moore, 1985), including southern Florida, lowlands in Mexico, Central America, Caribe, and South America. We retained only the studies that estimated genetic parameters using nuclear microsatellite markers to avoid the effect of different evolution modes among molecular markers.
We first classified plant species into two categories: species that exhibit megafauna dispersal syndrome, that is, those species that had their seeds potentially dispersed by megafauna, and species that do not exhibit megafauna syndrome, that is, those species that have their seeds dispersed by other factors (wind, water, authochory) or animal species (zoochory) but not megafauna (Table S1). Moreover, we classified plant species as megafauna dependent, that is, species that potentially relied only on megafauna to disperse, and megafauna independent, that is, species that relied on megafauna and other seed dispersers (Table S1). To classify the plant species, we used the definitions proposed by Guimarães et al. (2008) and an expert opinion (Mauro Galetti, personal communication). Thus, we followed the literature on megafauna syndrome classification, avoiding new classifications that could lead to data misinterpretation. Guimarães et al. (2008) defined megafauna syndrome fruits as small or large fleshy fruits (4–10 cm in diameter) with up to five large seeds (generally >2.0 cm diameter), which they call “type I fruits,” and large fleshy fruits (> 10 cm diameter) with many small seeds (> 100), which they call “type II fruits.” This classification excludes large fruits without a fleshy pulp. In our data set, megafauna-dependent species includes “type I fruits” (Table S1).
To model the effects of megafauna seed dispersal syndromes and other RT on variations in genetic diversity and genetic differentiation among populations, we compiled information about RT, such as mating and breeding systems and pollination mode. Trait data were obtained from botanical reviews and articles about pollination system and plant reproduction. We then classified the species according to their pollination syndrome as anemophily (pollinated by wind), chiropterophily (pollinated by bats), entomophily (pollinated by insects), and ornithophily (pollinated by birds). For mating systems, we found species with outcrossing and mixed systems. For breeding systems, we found monoecious, dioecious, and hermaphrodite species.
Genetic Parameters
To measure genetic differentiations, we compiled Wright’s (1951) FST. To measure genetic diversity, we obtained Nei’s (1978) genetic diversity (He) and AR (Mousadik and Petit, 1996). As articles often do not report the same genetic parameters, sample sizes might vary among the parameters analyzed.
Effect of Megafauna Seed Dispersal Syndromes in Genetic Parameters
To analyze the effect of megafauna seed dispersal syndromes on genetic diversity and differentiation of Neotropical plant species, we divided the whole database into three different sets to account for differences in seed size and morphology caused by the dispersal syndrome. For example, hydrochoric species usually have large seeds, whereas autochoric species usually have very small seeds compared to mammal-dispersed seeds. Set 1 included only species dispersed by mammals, set 2 included only species with zoochoric dispersal syndrome (including mammals), and set 3 included all species in the data sets, that is, zoochoric, wind-dispersed, autochoric, and hydrochoric species.
Reproductive and seed traits related to seed dispersal mode (e.g., seed size and number) in Angiosperm have significant phylogenetic signal and evolutionary constraints imposed by reproductive structures (Jordano, 1995). Many traits in seed plants have evolutionary interdependencies that may constrain evolutionary and adaptive responses (Donoghue, 1989). Thus, for each of the three data sets, we performed phylogenetic non-independent (PH-NoIND) and phylogenetic-independent (PH-IND) meta-analyses (Lajeunesse, 2009) using random-effect models, fitting species identity as a random factor. We compared the overall pooled effect sizes (δ) for each data set using Hedges’s d (Hedges and Olkin, 1985), which measures the magnitude and direction of experimental outcomes in standardized units. We obtained their variances and 95% confidence intervals (CI) and performed the Z-test for non-zero effects (Hedges, 1992; Hedges and Olkin, 1985). The Z value is the weighted sum of square caused by the regression model (Hedges, 1992). Random-effect models assume that replicates come from different distributions and add an additional variance component (τ) to each replicate as an estimate of between-replicate variances. We evaluated the fit of PH-NoIND and PH-IND meta-analysis models using Akaike information criterion (AIC) based on Butler and King (2004). The model with the lowest AICc (AIC corrected by sample size and number of parameters; Burnham and Anderson, 2002) was considered as the more plausible to explain observed patterns. Analyses were performed using phyloMeta v.1.2 (Lajeunesse, 2011) implemented in R v. 3.3.2 (R Core Team, 2015).
Phylogenetic Analysis
To account for phylogenetic non-independence of megafauna syndromes and RT on meta-analysis and on regression-based models (see below), we built a phylogenetic hypothesis of all species included in the analysis using the internal master tree Phylomatic treeR20120829 from the platform Phylomatic (Webb and Donoghue, 2005). In the absence of information about branch length, it was set to a value of 1. Most of the phylogenetic relationships within plants at deeper levels are well–resolved; therefore, polytomies at terminal nodes (species level) would not affect the results (Purvis and Garland, 1993; Swenson, 2009).
Syndromes in Genetic Parameters
Because PH-NoIND meta-analysis models fitted better for all three data sets (see Results), yielding very similar results to phylogenetic independent meta-analyses, we only mentioned important data for PH-IND in the main text but presented details in the Supporting Information (Figures S1 and S2, Tables S3–S5). We tested whether species with megafauna seed dispersal syndromes have a higher genetic differentiation and a lower genetic diversity compared to species without megafauna syndrome with PH-NoIND meta-analysis. For each data set and genetic parameter, we assessed the mean difference and 95% CI (Schwarzer et al., 2015) between species with megafauna seed dispersal syndromes and no megafauna syndromes. We also compared megafauna-dependent versus megafauna-independent plant species and megafauna-dependent versus species with no megafauna syndromes. Analyses were performed using meta package (Schwarzer et al., 2007; Schwarzer et al., 2015) implemented in R v. 3.3.2 (R Core Team, 2015). Moreover, we ran a power analysis to determine the statistical power of our comparison given the restrictions in sample size (Valentine et al., 2010; Quintana, 2015). This analysis was performed with metapower R-script (Quintana and Tiebel, 2018).
Modeling the Effect of Megafauna Seed Dispersal Syndromes and Reproductive Traits on Genetic Parameters
To investigate the effects of megafauna seed dispersal syndromes and RT on genetic diversity and differentiation, we used generalized linear mixed models (GLMM). Pollination mode and breeding system were treated as multistate categorical variables, whereas mating system (i.e., mixed or outcrossing) and megafauna syndromes (yes or no) were treated as binary variables. RT and megafauna syndromes were fitted as fixed factors, whereas species identity was fitted as random factors. When a species was studied with the same marker more than once, we calculated the average values for each genetic parameter (10 cases in 102 species). We fitted separate models for each genetic parameter (FST, He, and AR) and for each one of the three data sets. Analyses were carried out using MCMCglmm package (Hadfield, 2010) implemented in R version 3.3.2 (R Core Team, 2015) that uses a Bayesian framework with Markov Chain and Monte Carlo approximation algorithms. It was run with a Gausian distribution and a total of 80,000 iterations (burn-in 20,000 chains). Results were summarized by the mean of the posterior distribution and 95% CI, indicating direction, strength, and significance of the effects. Modeling of megafauna-dependent and -independent species was not possible because of the small sample size.
We used the phylogenetic tree to test whether megafauna seed dispersal syndromes and RT have a phylogenetic signal for the species in our data sets, that is, phylogenetic-related species tend to be more similar than expected by chance (Blomberg and Garland, 2002) using the Abouheif’s (Abouheif, 1999; Pavoinea et al., 2008) proximity test implemented in the R package adephylo (Jombart et al., 2010). Then, we fitted pGLS models (Martins and Hansen, 1997) to verify whether results obtained by the GLMM were robust enough so the pattern persisted after accounting for phylogenetic relationships. pGLS analyses were carried out using the package carper (Orme et al., 2012) in R environment.
Results
Effect of Megafauna Seed Dispersal Syndromes in Genetic Parameters
We used a database with a total of 102 species (Table S1) to perform PH-NoIND and PH-IND meta-analyses, comprising 98 studies (Table S2). All effect sizes were significantly different from zero for both metrics of meta-analyses and for all data sets (p < 0.001 for all comparisons, Figures S1 and S2, Tables S3–S5). PH-NoIND meta-analyses’ metrics showed lower AICc (Akaike information criterion corrected for sample size and number of parameters) and fitted better to our data for all analyses (Figures S1 and S2, Tables S3–S5).
Contrary to expectation, metrics of PH-NoIND meta-analyses showed no significantly higher genetic differentiation for species with megafauna seed dispersal syndrome than for species without megafauna syndrome (Figure 1A, Table S6) for all data sets. However, power analyses reveal that statistical power given both effect size and sample size is low (Tables S3–S5). Furthermore, none of the genetic diversity parameters studied (He and AR) was significantly lower for species with megafauna syndrome than for no megafauna syndrome (Figure 1A, Table S6). Megafauna-dependent and megafauna-independent plant species showed no significant differences in any of the genetic parameters studied (Figure 1B, Table S6). Moreover, megafauna-dependent and plants with no megafauna syndrome also showed no significant differences in any of the genetic parameters analyzed (Figure 1C, Table S6). Power analysis indicates that results for genetic diversity parameters (He and AR) have enough statistical power (Tables S3–S5) to support results of the meta-analysis.
Figure 1 Effect of megafauna seed dispersal syndromes on genetic parameters across the three data sets for plant species with megafauna and no megafauna syndromes based on phylogenetic-dependent meta-analysis. (A) Effect size for species with megafauna versus no megafauna syndrome dispersal. (B) Effect size for megafauna-dependent versus megafauna-independent species. (C) Effect size for megafauna-dependent versus no megafauna syndrome species.
Modeling the Effect of Megafauna Seed Dispersal Syndromes and Reproductive Traits on Genetic Parameters
RT rather than megafauna seed dispersal syndromes showed significant effects on genetic parameters for Neotropical plant species in GLMMs.
Genetic differentiation (FST) was lower for plant species with outcrossing mating system when all data sets or only zoochoric species were analyzed (Tables S7–S9). Genetic diversity (He) was significantly lower for species pollinated by insects when considering only mammal-dispersed or zoochoric species (Tables S7–S9) and significantly higher for species with outcrossing mating systems. AR was not significantly related to any RT studied.
Megafauna seed dispersal syndromes and pollination mode showed significant phylogenetic signals in our data sets (Table S10). However, mating and breeding systems had no phylogenetic signal most likely caused by incomplete taxon sampling because few species have been studied. Analyses accounting for phylogenetic non-independence (pGLS) confirmed the results of GLMM for genetic differentiation (Tables S11–S12). FST values are lower for species with outcrossing mating systems, considering all data sets and zoochoric species data set (Tables S11–S12). Genetic diversity (He) was significantly higher for species with outcrossing system for all data sets (Tables S11–S12), but the relationship of He and pollination mode was not recovered.
Discussion
Our findings show no significant signal of megafauna extinction in genetic diversity and differentiation of Neotropical plants. Plants with megafaunal seed dispersal syndromes showed similar values of genetic differentiation and genetic diversity compared to species without megafaunal dispersal syndromes. Moreover, the lack of genetic signal was not caused by the dependency level on megafauna dispersal as megafauna-dependent species showed no difference in genetic parameters when compared to species with non-megafaunal syndromes or to megafauna-independent species. It is important to note that statistical power was low for genetic differentiation comparisons but very high for the other genetic parameters.
The lack of genetic signal may be caused by the dispersal rescue effect of other disperser species. Extant mammals, such as tapir and deer, may move seeds over long distances and may have long enough gut retention, potentially leading to long-distance dispersal (e.g., Vellend et al., 2003; Clauss et al., 2010). In addition, the role of scatter-hoarding rodents in long-distance dispersal has been underestimated until recently (Jansen et al., 2012), and these mammals may be important in long-distance dispersal to many Neotropical plant species, counterbalancing the past role of the extinct megafauna, maintaining genetic diversity and differentiation among populations at the same levels of species with non-megafaunal syndromes.
Furthermore, pollination may be more important than seed dispersal in shaping patterns in genetic differentiation and diversity in the Neotropics. In fact, many works show higher contribution of pollen dispersal than seed dispersal to gene flow in Angiosperms (e.g., Petit et al., 2005) and in Neotropical plants (e.g., Hamilton and Miller, 2002; Collevatti et al., 2003; Ballesteros-Mejia et al., 2016). Indeed, our findings show significant effects of RT, such as pollination mode, mating, and breeding systems, in genetic differentiation and diversity, strengthening the importance of RT in long-distance dispersal (Ballesteros-Mejia et al., 2016).
It is important to note that the lack of evidence of an effect of megafauna seed dispersal syndrome on genetic diversity and differentiation does not represent evidence of a lack of effect of megafauna because confounding effects, such as species demographic history, may hinder the detection of any genetic signal of loss of long-distance seed dispersal. The extinction of megafauna may be relatively recent compared to the generation time of the plant species and may not be long enough to generate a genetic effect. Most species studied so far are trees (89 in 102 species) with 15–20 years of generation time (i.e., time to first reproduction). The last megafauna extinct at c. 6 ka, corresponding to c. 400 generations, which may not be enough to reduce genetic diversity and increase genetic differentiation in long-lived trees. In addition, many Neotropical trees experienced retraction or expansion in both their geographical range and effective population size during glacial periods, such as the species with megafaunal seed dispersal syndromes Caryocarbrasiliense (Collevatti et al., 2012), Dipteryxalata (Collevatti et al., 2013), and Mauritia flexuosa (Lima et al., 2014) and the species with non-megafaunal syndromes Tabebuia impetiginosa (Collevatti et al., 2012), Dalbergiamiscolobium (Novaes et al., 2013), Tabebuia aurea (Collevatti et al., 2015), and Tabebuia roseoalba (Melo et al., 2016). Thus, demographic history may have also shaped the current geographical distribution of genetic diversity and differentiation, bewildering the effects of megafauna extinction. The data available so far have not enough statistical power to detect effects of megafauna syndrome in genetic differentiation. Hence, future studies may shed more light on the effects of megafauna extinction in genetic differentiation.
It is worth noticing however that the definition of megafauna seed dispersal syndromes is controversial because it includes fruits and seeds with highly contrasting characteristics (Howe, 1985). For instance, plants with fruits and large seeds, such as Mauritia flexuosa, Dipteryxalata, Caryocarbrasiliense, are considered species with megafauna seed dispersal syndromes but so are plants with large or medium fruits with small hard-coated seeds, such as Annona crassiflora, Annanascomosus, Enterolobiumcyclocarpum, and Solanum lycocarpum (Janzen and Martin, 1982; Guimarães et al., 2008). These megafauna syndrome species with very different fruit and seed characteristics may have used different dispersal agents for long-distance dispersal to counterbalance the extinction of their megafauna mutualistic seed dispersers or may rely on secondary dispersers that improve long-distance dispersal (Levey et al., 2002). In addition, megafauna seed dispersal syndromes may be adaptations to another environmental condition. For instance, large seed size may also be the outcome of adaptation to drought or nutrient-poor soils (Westoby et al., 1996; Moles and Westoby, 2004), which may explain large seeds in many savanna species. Species with megafauna syndromes may also have survived the extinction of megafauna via other traits that also promote long-distance dispersal and recruitment but are not related to fruit or seed characteristics (Johnson, 2009). For example, many Neotropical species have high capacity for vegetative propagation or sprouting or are long lived, such as most savanna tree species, or are able to establish beneath the parent plant (Donatti et al., 2007; Collevatti and Hay, 2011).
In conclusion, the genetic studies of Neotropical plants show no difference in genetic diversity and differentiation in species with megafaunal seed dispersal syndromes compared to non-megafaunal. Rather, a combination between RT, such as pollination mode and mating and breeding systems, along with demographic history may be more important in shaping the current genetic diversity and differentiation patterns in Neotropical plants.
Author Contributions
RC conceived and funded the work. RC and JL obtained the data. LB-M and JL carried out analyses. RC wrote the original draft, and all authors contributed to the manuscript and approved the final version.
Funding
This work was supported by the research network Rede Cerrado CNPq/PPBio (project no. 457406/2012-7), PROCAD/CAPES (project no 88881.068425/2014-01) and CAPES Ciências sem Fronteira (project CSF-PAJT/CAPES no. 88881.030318/2013-01). JSL receives a fellowship from PROCAD/CAPES and LBM receives a fellowship from CAPES Ciências sem Fronteira.
Conflict of Interest Statement
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.
Acknowledgments
RC has continuously been supported by CNPq and CAPES grants and scholarships whose assistance we gratefully acknowledge.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00788/full#supplementary-material
Supporting information
Supporting Information S1 | Tables S1–S12 with information of data set and meta-analysis and GLMM and pGLS results.
Supporting Information S2 | Figures S1 and S2 with information of meta-analysis results.
References
Abouheif, E. (1999). A method for testing the assumption of phylogenetic independence in comparative data. Evol. Ecol. Res. 1, 895–909.
Ballesteros-Mejia, L., Lima, N. E., Lima-Ribeiro, M. S., Collevatti, R. G. (2016). Pollination mode and mating system explain patterns in genetic differentiation in Neotropical plants. PLoS One 11, e0158660. doi: 10.1371/journal.pone.0158660
Bawa, K. S. (1992). Mating systems, genetic differentiation and speciation in tropical rain forest plants. Biotropica 24, 250–255. doi: 10.2307/2388519
Blomberg, S. P., Garland, T., Jr. (2002). Tempo and mode in evolution: phylogeneticinertia, adaptation and comparative methods. J. Evol. Biol. 15, 899–910. doi: 10.1046/j.1420-9101.2002.00472.x
Bodmer, R., Ward, D. (2006). Frugivory in large mammalian herbivores, in Large Herbivore Ecology, Ecosystem Dynamics and Conservation. Eds. Danell, K., Bergstrom, R., Duncan, P., Pastor, J. (Cambridge: Cambridge University Press), p. 232–260. doi: 10.1017/CBO9780511617461.010
Burnham, K. K. P., Anderson, D. R. D. (2002). Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer.
Butler, M. A., King, A. A. (2004). Phylogenetic comparative analysis: a modeling approach for adaptive evolution. Am. Nat. 164, 683–695. doi: 10.1086/426002
Clauss, M., Schwarm, A., Ortmann, S., Streich, W. J., Hummel, J. (2007). A case of non-scaling in mammalian physiology? Body size, digestive capacity, food intake, and ingesta passage in mammalian herbivores. Comp. Biochem. Physiol. A Physiol. 148, 249–265. doi: 10.1016/j.cbpa.2007.05.024
Clauss, M., Kaiser, T., Hummel, J. (2008). The morphophysiological adaptations of browsing and grazing mammals, in The Ecology of Browsing and Grazing. Eds. Gordon, I. J., Prins, H. H. T. (Heidelberg: Springer), 47–88. doi: 10.1007/978-3-540-72422-3_3
Clauss, M., Lang-Deuerling, S., Müller, D. W. H., Kienzle, E., Steuer, P., Hummel, J. (2010). Retention of fluid and particles in captive tapirs (Tapirus sp). Comp. Biochem. Physiol. A Physiol. 157, 95–101. doi: 10.1016/j.cbpa.2010.03.029
Collevatti, R. G., Grattapaglia, D., Hay, J. D. (2003). Evidences for multiple maternal lineages of Caryocarbrasiliense populations in the Brazilian Cerrado based on the analysis of chloroplast DNA sequences and microsatellite haplotype variation. Mol. Ecol. 12, 105–115. doi: 10.1046/j.1365-294X.2003.01701.x
Collevatti, R. G., Hay, J. D. (2011). Kin structure and genotype-dependent mortality: a study using the Neotropical tree Caryocarbrasiliense. J. Ecol. 99, 757–763. doi: 10.1111/j.1365-2745.2011.01796.x
Collevatti, R. G., Lima-Ribeiro, M. S., Souza-Neto, A. C., Franco, A. A., Oliveira, G., Terribile, L. C. (2012). Recovering the demographical history of a Brazilian cerrado tree species Caryocarbrasiliense: coupling ecological niche modeling and coalescent analyses. Nat. Conservacao 10, 169–176. doi: 10.4322/natcon.2012.024
Collevatti, R. G., Terrible, L. V., Lima-Ribeiro, M. S., Nabout, J. C., Oliveira, G., Rangel, T. F., et al. (2012). A coupled phylogeographic and species distribution modeling approach recovers the demographic history of a Neotropical seasonally dry forest tree species. Mol. Ecol. 21, 5845–5863. doi: 10.1111/mec.12071
Collevatti, R. G., Telles, M. P. C., Nabout, J. C., Chaves, J. L., Soares, T. N. (2013). Demographic history and the low genetic diversity in Dipteryxalata (Fabaceae) from Brazilian Neotropical savannas. Heredity 111, 97–105. doi: 10.1038/hdy.2013.23
Collevatti, R. G., Terribile, L. C., Rabelo, S. G., Lima-Ribeiro, M. S. (2015). Relaxed random walk model coupled with ecological niche modeling unravel the dispersal dynamics of a Neotropical savanna tree species in the deeper Quaternary. Front. Plant Sci. 6, 653. doi: 10.3389/fpls.2015.00653
Cox, C. B., Moore, P. D. (1985). Biogeography: an ecological and evolutionary approach. 4th ed. London: Blackwell Scientific Publications, p. 244.
Donatti, C. I., Galetti, M., Pizo, M. A., Guimarães, P. R., Jr., Jordano, P. (2007). Living in the land of ghosts: fruit traits and the importance of large mammals as seed dispersers in the Pantanal, Brazil, in Frugivory and seed dispersal: theory and applications in a changing world. Eds. Dennis, A., Green, R., Schupp, E. W., Wescott, D. A. (Wallingford: Commonwealth Agricultural Bureau International), 104–123. doi: 10.1079/9781845931650.0104
Donoghue, M. J. (1989). Phylogenies and the analysis of evolutionary sequences, with examples from seed plants. Evolution 43, 1137–1156. doi: 10.1111/j.1558-5646.1989.tb02565.x
Gill, J. L., Williams, J. W., Jackson, S. T., Lininger, K. B., Robinson, G. S. (2009). Pleistocene megafaunal collapse, novel plant communities, and enhanced fire regimes in North America. Science 326, 1100–1103. doi: 10.1126/science.1179504
Gill, J. L. (2014). Ecological impacts of the late Quaternary megaherbivore extinctions. New Phytol. 201, 1163–1169. doi: 10.1111/nph.12576
Guimarães, P. R., Jr., Galetti, M., Jordano, P. (2008). Seed dispersal anachronisms: rethinking the fruits extinct megafauna ate. PLoS One 3, e1745. doi: 10.1371/journal.pone.0001745
Hamilton, M. B., Miller, J. R. (2002). Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure. Genetics 162, 1897–1909.
Hansen, D. M., Galetti, M. (2009). The forgotten megafauna. Science 324, 42–43. doi: 10.1126/science.1172393
Hadfield, J. D. (2010). MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J. Stat. Softw. 33, 1–22. doi: 10.18637/jss.v033.i02
Hester, A. J., Bergman, M., Iason, G. R., Moen, J. (2006). Impacts of large herbivores on plant community structure and dynamics, in Large herbivore ecology, ecosystem dynamics and conservation. Eds. Danell, K., Bergstrom, R., Duncan, P., Pastor, J. (Cambridge: Cambridge University Press), p. 97–141. doi: 10.1017/CBO9780511617461.006
Janzen, D. H., Martin, P. S. (1982). Neotropical anachronisms: the fruits the Gomphoteres ate. Science 215, 19–27. doi: 10.1126/science.215.4528.19
Jansen, P. A., Hirsch, B. T., Emsens, W. J., Zamora-Gutierrez, V., Wikelski, M., Kays, R. (2012). Thieving rodents as substitute dispersers of megafaunal seeds. Proc. Natl. Acad. Sci. 109, 12610–12615. doi: 10.1073/pnas.1205184109
Johnson, C. N. (2009). Ecological consequences of Late Quaternary extinctions of megafauna. Proc. R. Soc. B 276, 2509–2519. doi: 10.1098/rspb.2008.1921
Jombart, T., Balloux, F., Dray, S. (2010). adephylo: new tools for investigating the phylogenetic signal in biological traits. Bioinformatics 26, 1907–1909. doi: 10.1093/bioinformatics/btq292
Jordano, P. (1995). Angiosperm fleshy fruits and seed dispersers: a comparative analysis of adaptation and constraints in plant–animal interactions. Am. Nat. 145, 163–191. doi: 10.1086/285735
Lajeunesse, M. J. (2009). Meta-analysis and the comparative phylogenetic method. Am. Nat. 174, 369–381. doi: 10.1086/603628
Lajeunesse, M. J. (2011). phyloMeta: a program for phylogenetic comparative analyses with meta-analysis. Bioinformatics 27, 2603–2604. doi: 10.1093/bioinformatics/btr438
Levey, D. J., Silva, W. R., Galetti, M. (2002). Seed Dispersal and Frugivory: ecology, evolution, and conservation. Wallingford, UK: CAB International, 500p.
Lima, E. N., Lima-Ribeiro, M. S., Tinoco, C. F., Terribile, L. C., Collevatti, R. G. (2014). Phylogeography and ecological niche modelling, coupled with the fossil pollen record, unravel the demographic history of a Neotropical swamp palm through the Quaternary. J. Biogeogr. 41, 673–686. doi: 10.1111/jbi.12269
Malhi, Y., Doughty, C. E., Galetti, M., Smith, F. A., Svenning, J. C., Terborgh, J. W. (2016). Megafauna and ecosystem function from the Pleistocene to the Anthropocene. Proc. Natl. Acad. Sci. 113, 838–846. doi: 10.1073/pnas.1502540113
Martins, E. P., Hansen, T. F. (1997). Phylogenies and the comparative method: a general approach to incorporating phylogenetic information into analysis of interspecific data. Am. Nat. 149, 646–667. doi: 10.1086/286013
Melo, W. A., Lima-Ribeiro, M. S., Terribile, L. C., Collevatti, R. G. (2016). Coalescent simulation and paleodistribution modeling for Tabebuia rosealba do not support South American dry forest refugia hypothesis. PLoS One 11, e0159314. doi: 10.1371/journal.pone.0159314
Moles, A. T., Westoby, M. (2004). Seedling survival and seed size: a synthesis of the literature. J. Ecol. 92, 372–383. doi: 10.1111/j.0022-0477.2004.00884.x
Mousadik, A., Petit, R. J. (1996). High level of genetic differentiation for allelic richness among populations of the argan tree [Arganiaspinosa (L). Skeels] endemic to Morocco. Theor. Appl. Genet. 92, 832–839. doi: 10.1007/BF00221895
Nei, M. (1978). Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89, 583–590.
Novaes, R. M. L., Ribeiro, R. A., Lemos-Filho, J. P., Lovato, M. B. (2013). Concordance between phylogeographical and biogeographical patterns in the Brazilian cerrado: diversification of the endemic tree Dalbergiamiscolobium (Fabaceae). PLoS One 8, e82198. doi: 10.1371/journal.pone.0082198
Orme, C. D. L., Freckleton, R. P., Thopas, G. H., Petzoldt, T., Fritz, S. A. (2012). The caper package: comparative analysis of phylogenetics and evolution in R. Available at ftP://cran.rproject.org/pub/R/web/packages/caper)vignettes/caper.pdf.
Pavoinea, S., Ollierb, S., Pontiera, D., Chessela, D. (2008). Testing for phylogenetic signal in phenotypic traits: new matrices of phylogenetic proximities. Theor. Popul. Biol. 73, 79–91. doi: 10.1016/j.tpb.2007.10.001
Petit, R. J., Duminil, J., Fineschi, S., Hampe, A., Salvini, D., Vendramin, G. G. (2005). Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Mol. Ecol. 14, 689–711. doi: 10.1111/j.1365-294X.2004.02410.x
Pires, M. M., Guimarães, P. R., Galetti, M., Jordano, P. (2017). Pleistocene megafaunal extinctions and the functional loss of long-distance seed-dispersal services. Ecography 41, 153–163. doi: 10.1111/ecog.03163
Pires, M. M., Galetti, M., Donatti, C. I., Pizo, M. A., Dirzo, R., Guimarães, P. R. (2014). Reconstructing past ecological networks: the reconfiguration of seed-dispersal interactions after megafaunal extinction. Oecologia 175, 1247–1256. doi: 10.1007/s00442-014-2971-1
Purvis, A., Garland, T., Jr. (1993). Polytomies in comparative analyses of continuous characters. Syst. Biol. 42, 569–575. doi: 10.1093/sysbio/42.4.569
Quintana, D. (2015). From pre-registration to publication: a non-technical primer for conducting a meta-analysis to synthesize correlational data. Front. Psychol. 6, 1549. doi: 10.3389/fpsyg.2015.01549
Quintana, D., Tiebel, J. (2018). How to calculate statistical power for your meta-analysis. Accessed April 2019. Free available at https://osf.io/w4xrs/.
R Core Team. (2015). R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing, ISBN: ISBN. 3-900051-07-0.
Schwarzer, G., Carpenter, J. R., Rücker, G. (2015). Meta-Analysis with R (Use-R)!. Switzerland: Springer International Publishing. Available at http://www.springer.com/gp/book/9783319214153. doi: 10.1007/978-3-319-21416-0
Swenson, N. G. (2009). Phylogenetic resolution and quantifying the phylogenetic diversity and dispersion of communities. PLoS One 4, e4390. doi: 10.1371/journal.pone.0004390
Tabarelli, M., Peres, C. A. (2002). Abiotic and vertebrate seed dispersal in the Brazilian Atlantic forest: implications for forest regeneration. Biol. Conserv. 106, 165–176. doi: 10.1016/S0006-3207(01)00243-9
Webb, C. O., Donoghue, M. J. (2005). Phylomatic: tree assembly for applied phylogenetics. Mol. Ecol. Notes 5, 181–183. doi: 10.1111/j.1471-8286.2004.00829.x
Westoby, M., Leishman, M., Lord, J. (1996). Comparative ecology of seed size and dispersal. Philos. Trans. R. Soc. B. 351, 1309–1318. doi: 10.1098/rstb.1996.0114
Wright, S. (1951). The genetic structure of populations. Annu. Eugen. 15, 323–354. doi: 10.1111/j.1469-1809.1949.tb02451.x
Valentine, J. C., Pigott, T. D., Rothstein, H. R. (2010). How many studies do you need?: a primer on statistical power for meta-analysis. J. Educ.Behav. Stat. 35, 215–247. doi: 10.3102/1076998609346961
Van Soest, P. J. (1996). Allometry and ecology of feeding behavior and digestive capacity in herbivores: a review. Zool. Biol. 15, 455–479. doi: 10.1002/(SICI)1098-2361(1996)15:5<455::AID-ZOO3>3.0.CO;2-A
Keywords: angiosperms, genetic diversity, general linear mixed models, meta-analysis, Neotropics, phylogenetic generalized least square models, phylogenetic meta-analysis, phylogenetic signal
Citation: Collevatti RG, Lima JS and Ballesteros-Mejia L (2019) Megafauna Seed Dispersal in the Neotropics: A Meta-Analysis Shows No Genetic Signal of Loss of Long-Distance Seed Dispersal. Front. Genet. 10:788. doi: 10.3389/fgene.2019.00788
Received: 11 December 2018; Accepted: 25 July 2019;
Published: 05 September 2019.
Edited by:
Margaret Byrne, Department of Biodiversity, Conservation and Attractions (DBCA), AustraliaReviewed by:
Rodrigo Medel, University of Chile, ChileAnna Traveset, Spanish National Research Council (CSIC), Spain
Copyright © 2019 Collevatti, Lima and Ballesteros-Mejia. 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: Rosane G. Collevatti, rosanegc68@hotmail.com