Molecular Phylogeny and Historical Biogeography of the Neotropical Swarm-Founding Social Wasp Genus Synoeca (Hymenoptera: Vespidae)

Abstract

The Neotropical Region harbors high biodiversity and many studies on mammals, reptiles, amphibians and avifauna have investigated the causes for this pattern. However, there is a paucity of such studies that focus on Neotropical insect groups. Synoeca de Saussure, 1852 is a Neotropical swarm-founding social wasp genus with five described species that is broadly and conspicuously distributed throughout the Neotropics. Here, we infer the phylogenetic relationships, diversification times, and historical biogeography of Synoeca species. We also investigate samples of the disjoint populations of S. septentrionalis that occur in both northwestern parts of South America through Central American and the Brazilian Atlantic rainforests. Our results showed that the interspecific relationships for Synoeca could be described as follows: (S. chalibea + S. virginea) + (S. cyanea + (S. septentrionalis/S. surinama)). Notably, samples of S. septentrionalis and S. surinama collected in the Atlantic Forest were interrelated and may be the result of incomplete lineage sorting and/or mitochondrial introgression among them. Our Bayesian divergence dating analysis revealed recent Plio-Pleistocene diversification in Synoeca. Moreover, our biogeographical analysis suggested an Amazonian origin of Synoeca, with three main dispersal events subsequently occurring during the Plio-Pleistocene.

About the Authors:

Rodolpho Santos Telles Menezes

* E-mail: [email protected]

Affiliations Departamento de Ciências Biológicas, Universidade Estadual de Santa Cruz, Ilhéus, Bahia, Brazil, Department of Entomology, National Museum of Natural History, Smithsonian Institution, Washington, DC, United States of America

Seán Gary Brady

Affiliation: Department of Entomology, National Museum of Natural History, Smithsonian Institution, Washington, DC, United States of America

Antônio Freire Carvalho

Affiliation: Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, São Paulo, Brazil

Marco Antonio Del Lama

Affiliation: Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, São Paulo, Brazil

Marco Antônio Costa

Affiliation: Departamento de Ciências Biológicas, Universidade Estadual de Santa Cruz, Ilhéus, Bahia, Brazil

Introduction

The Neotropical Region is one of the most important biodiversity hotspots in the world, harboring more than half of the Earth’s remaining rainforests [1]. Recent studies documenting the diversification history of the Neotropical biota have elucidated both temporal and spatial biogeographic patterns, resulting in proposals that diversification events in the South American biota have been driven by the uplift of the Andes, marine incursions, or Pleistocene climate changes [2–5]. In addition, studies involving mammals and birds have allowed the inference of historical connections between the morphoclimatic domains of the South America Amazonian rainforest (AM) and the Atlantic Forest (AF) (e.g., [6–7]). Currently these two forests are separated by a dry corridor of open vegetation [8], also called the ‘Dry Diagonal’, that encompasses the Argentinean and Paraguayan Chaco, the Caatinga in northeastern Brazil, and the central Brazilian Cerrado [6,8].

Synoeca de Saussure, 1852 (Hymenoptera: Vespidae: Epiponini) is a Neotropical swarm-founding social wasp genus with five described species: Synoeca chalibea de Saussure, 1852; Synoeca virginea (Fabricius, 1804); Synoeca surinama (Linnaeus, 1767); Synoeca septentrionalis Richards, 1978; and Synoeca cyanea (Fabricius, 1775). Species in this genus display vigorous defense of their nests accompanied by painful sting. Their nests are directly attached to the substrate and have a convex envelope that may present regular transverse corrugations resembling the back of an armadillo (Dasypodidae), and so are known in many parts of Latin America as “cachicamas”, “armadillas”, “conchajonas”, “carachupa”, “caba-tatu” or “marimbondo-tatu” [9–14]. Synoeca species have been the focus of many studies about their behavior [10, 12–15], caste differentiation [16–17] and nest architecture [11, 18]. They are broadly distributed in the Neotropical region, ranging from central Mexico down to northern Argentina [9]. Synoeca virginea is restricted to AM; S. chalibea also occurs in AM, but can be found in Costa Rica and Panama as well. Synoeca surinama is widely distributed in South American rainforests and, conspicuously, in the Brazilian savannah (i.e., the Cerrado) as well. Synoeca septentrionalis is distributed from northwestern parts of South America through Central America and into central Mexico; however, recently specimens have been found in the AF. Synoeca cyanea is the only species that is restricted to the eastern portion of South America (ESA) [9, 19].

Several previous studies have sought to understand the phylogenetic relationships among Synoeca species [20–22]. These analyses used only morphological characters and supported the monophyly of the genus, with the following relationships among the species: S. chalibea + (S. virginea + (S. septentrionalis + (S. surinama + S. cyanea))) [20, 22]. However, species of Synoeca exhibit wide morphological variation, as indicated by Cely and Sarmiento [21] and the presence of multiple synonyms for species [9]. Moreover, the recent record of colonies of S. septentrionalis in the AF [19] raises questions about the taxonomy and evolutionary history of the group.

Here we infer the first molecular phylogeny for Synoeca species using mitochondrial and nuclear DNA sequences. We also provide new observations of morphological structures that can help diagnose some species. Furthermore, we investigate the status of the disjoint populations of S. septentrionalis. Bayesian divergence dating analysis calibrated with an a priori mutation rate of insect mitochondrial DNA was used to estimate a timescale for diversification events in the genus. We also provide a Bayesian analysis of historical biogeography that suggests an Amazonian origin for Synoeca.

Material and Methods

Taxon sampling

Our sampling included specimens from multiple localities within known species distributions. Table 1 and Fig. 1 show the collection sites of specimens used in this study. The map was generated using the software Quantum-GIS v1.8.0 (Open Source Geospatial Foundation Project, Beaverton, OR, USA). All specimens were preserved in ethanol prior to the molecular analyses and vouchers are deposited in the entomological collections at the Universidade Estadual de Santa Cruz, Ilhéus, Bahia, Brazil and the National Museum of Natural History (USNM, Washington, DC). All necessary research permits for fieldwork and collection of samples by RSTM and AFC were issued by the Brazilian Institute for Biodiversity Conservation (ICMBio) recorded by SISBio (permit numbers 22180 and 5374523). Field studies did not involve endangered or protected species. The samples from localities below were donated from different researchers: Wanda, Misiones, Argentina: 25°58'55" S, 54°34'58" W; El Tuito, Jalisco, Mexico: 20°19'02" N, 105°19'42" W; Chiapas, Mexico: 16°48'05" N, 93°05'18" W; Pandi, Cundinamarca, Colombia: 4°;15'32" N, 74°36'35" W; Santa Rosa, Guanacaste, Costa Rica: 10°50'06" N, 85°42'12" W".

[Figure omitted. See PDF.]

Table 1. Descriptive information for all taxa used in this phylogenetic study.

https://doi.org/10.1371/journal.pone.0119151.t001

[Figure omitted. See PDF.]

Fig 1. Map showing collection sites.

The numbers refer to the sampling locality codes shown in Table 1.

https://doi.org/10.1371/journal.pone.0119151.g001

Species identification and morphological observations

Species identification was based on the keys provided by Richards [9], Andena et al. [20], and Cely and Sarmiento [21]. We observed additional external morphological structures which helped in identifying specimens, as follows: punctation on the head; erect setae on the scape; and erect setae on the scutum (Fig. 2). Images of specimens were generated using a JVC KY–F75U digital camera mounted on a Leica Z16 APO steromicroscope. All images were edited using Photoshop CS4 (Version 11.0) (Adobe Inc.).

[Figure omitted. See PDF.]

Fig 2. External morphological structures observed.

View frontal of head, (a) S. surinama and (b) S. septentrionalis. View frontal of scape, (c) S. cyanea and (d) S. septentrionalis. Scutum in lateral view, (e) S. surinama and (f) S. septentrionalis.

https://doi.org/10.1371/journal.pone.0119151.g002

DNA extraction, amplification and sequencing

The mesosoma and/or hind legs were removed from each specimen and DNA was extracted using the phenol-chloroform method following the Han and McPheron protocol [23]. We amplified three mitochondrial (16S, cytochrome b, and cytochrome c oxidase I) and one nuclear (wingless) gene fragments by PCR using specific primers and amplification conditions (S1 Table). PCR products were purified using exonuclease I and shrimp alkaline phosphatase and directly sequenced in an ABI Prism 3730 (Applied Biosystems) sequencer (Laboratório de Biotecnologia da FCAV—UNESP de Jaboticabal, SP). PCR products were sequenced in both directions and sequence contigs were assembled using Sequencher 5.1 (Gene Code Corp., Ann Arbor, MI, USA). DNA sequences were aligned using Muscle 3.7 [24] (with default parameters) in MEGA 5.10 [25], with each of the four genes aligned separately. All sequences are deposited in GenBank and accession numbers are listed in Table 1.

Phylogenetic analyses

We included DNA data from 26 Synoeca samples with each of the five Synoeca species represented by multiple collection localities. We used species from four other genera as outgroups. Most phylogenetic analyses were based on a concatenated data matrix (1829 base pairs) of the four gene fragments. The most appropriate model of nucleotide evolution and the best-fitting partitioning scheme were selected using PartitionFinder v1.1.1 [26] under the Bayesian information criterion (BIC) (Table 2). Phylogenetic inference was conducted by Bayesian inference (BI) performed using MrBayes v3.2.2 [27] and maximum likelihood (ML) performed using GARLI v2.0 [28]. BI was also performed on each of the three mitochondrial gene fragments (wingless contained little phylogenetic information) separately to examine potential conflicts in phylogenetic signal among genes. All BI analyses consisted of two independent runs of 50 million generations each with four chains (temp = 0.1) and sampled every 1000 generations. The burn-in, convergence, and stationarity were assessed using Tracer v1.5 [29]. We removed the first 20% of sampled generations and combined the remaining generations to produce the maximum credibility tree. We conducted 1000 ML bootstrap replicates in GARLI under the same partitions and nucleotide models as in BI. Trees from all analyses were visualized using FigTree v1.4.0 program [30].

[Figure omitted. See PDF.]

Table 2. Data partitions and nucleotide substitution models suggested by PartitionFinder v1.1.1.

https://doi.org/10.1371/journal.pone.0119151.t002

We conducted an additional analysis in which we combined the information from all gene trees into a single tree (species tree), since data from multiple genes and multiple individuals per species can be useful for resolving species trees [31–32]. We used *BEAST v1.8.0 [33] to infer species tree. The recognized five species and S. septentrionalis samples from AF were assigned as ‘‘species” in the analysis. The *BEAST run consisted of 50 million generations, a Yule process for the species tree prior, a piecewise linear and constant root model for population size, randomly generated starting trees for each gene, and a burn-in of 20%.

Divergence time estimation

We inferred divergence times under a Bayesian framework using BEAST v.1.8.0 [34]. We generated the input file in BEAUTi using the two mitochondrial protein-coding genes (COI + CytB) and the substitution model (GTR + Γ) as selected by PartitionFinder. Only nucleotide data from ten specimens were included in these analyses in order to avoid missing data. We employed an uncorrelated lognormal relaxed clock model [35]. Clock models were unlinked, and substitution and tree models were linked among partitions. A Yule speciation process with a random starting tree was used for the tree prior. Given the poorly known fossil record for social wasps in general, with none belonging to Synoeca [36], we applied the Brower [37] mutation rate of mitochondrial genes (under a normal distributed prior). This mutation rate estimated at 2.3% My^-1 was based on a set of seven studies that provided age estimates of lineage splits ranging from 300 to 3,250,000 years ago. Two independent Markov chain Monte Carlo (MCMC) searches were conducted with 100 million generations each, with parameters sampled every 10,000 steps and a burn-in of 20%. We checked for convergence between runs and analysis performance with Tracer v1.5 using effective sample size (ESS) scores. The resulting trees were combined using TreeAnnotator v1.8.0 and the consensus tree with the divergence times was visualized in FigTree v1.4.0.

Historical biogeographic analysis

We performed the Bayesian Binary MCMC (BBM) method of biogeographical and ancestral state reconstruction implemented in RASP (Reconstruct Ancestral State in Phylogenies) 2.1b [38]. We used the tree obtained from *BEAST (species tree) and published occurrence data for the analyzed species [9] as input files for RASP. Thus, we assigned species distribution areas to geographical regions as follows: (A) Amazonian forest, (B) Middle America, (C) Atlantic forest, (D) Dry Diagonal (Cerrado, Chaco and Caatinga). The BBM analysis was run applying the model F81 + Γ and no outgroup was defined. We ran the analysis for 5 millions generations, sampled every 1000 generations with the first 1000 samples being discarded as burn-in.

Results

Morphological observations

We observed the presence of much punctation on the head of Synoeca septentrionalis and S. chalibea, in contrast to S. virginea, S. cyanea, and S. surinama, which showed little or no punctation on the head (Fig. 2a, b). Furthermore, S. septentrionalis have many erect setae on the scape (Fig. 2d), as opposed to S. virginea, S. chalibea, and S. cyanea, which do not present erect setae (Fig. 2c); S. surinama have few erect hairs on the scape. Another noteworthy morphological character is the presence of erect setae on the scutum of S. septentrionalis, S. cyanea, and S. chalibea that are absent from S. surinama and S. virginea (Fig. 2e, f). Moreover, the two S. chalibea samples studied here showed differences in body color, one form with a dark body (“dark form”) and the other form with the body entirely yellowish (“yellow form”).

Molecular phylogeny

Our concatenated data matrix of 1829 aligned nucleotide sites contained 416 variable sites, including 521 aligned base pairs (bp) of 16S (28.5% of data set), 472 bp of COI (25.8%), 433 bp of Cytb (23.7%), and 403 bp of Wg (22%). The BI and ML results from this concatenated data matrix are summarized in Fig. 3. The two types of analyses inferred the same topology with differences only in the degree of support for some nodes. Monophyly of the genus Synoeca is supported with high support [posterior probability (PP): 1 and bootstrap (BS): 0.99; Fig. 3]. The genus is divided into two major clades with high node support (PP: 1 and BS: 0.99; Fig. 3) containing: (1) S. chalibea and S. virginea; and (2) S. cyanea, S. septentrionalis, and S. surinama. Within the latter clade, S. cyanea is sister to S. septentrionalis and S. surinama with weak support (PP: 0.87 and BS: 0.52; Fig. 3). Synoeca septentrionalis samples from Middle America [henceforth called MA denoting samples from the Central America Rainforest region (Costa Rica), Southern North America (Mexico) and Northwestern South America (Colombia)] do not cluster with S. septentrionalis samples from AF. S. surinama samples from Acre, Amapá, and Distrito Federal (this group henceforth called AM/C denoting samples from both Amazon forest and Cerrado) do not cluster with S. surinama from AF. Moreover, S. surinama samples from the AF formed distinct clusters with S. septentrionalis from AF as following: the S. surinama sample from Bahia clusters with S. septentrionalis samples from Bahia, Espírito Santo, and São Paulo (this group henceforth called CAF denoting samples from Central Atlantic forest); the S. surinama sample from Sergipe clusters with S. septentrionalis sample from Pernambuco (this group henceforth called NAF denoting samples from the North Atlantic forest) (Fig. 3).

[Figure omitted. See PDF.]

Fig 3. Phylogeny of Synoeca reconstructed by Bayesian analysis.

Shown is the consensus tree resulting from analyses based on 1829 bp from 16S, CytB, COI and wingless. Maximum likelihood (ML) analysis resulted in the same topology. Numbers above or below branches indicate Bayesian posterior probabilities (first number) and ML bootstrap values (second number). Numbers in parentheses indicate the collection sites as shown in Fig. 1. Abbreviations for biogeographical units are as follows: AM, Amazonian forest; AM/C, Amazonian forest and Cerrado; NAF, North Atlantic forest; CAF, Central Atlantic forest; MA, Middle America; ESA, Eastern portion of South America.

https://doi.org/10.1371/journal.pone.0119151.g003

BI performed on each of the three mitochondrial gene fragments all show divergence into the two major clades in Synoeca and consistently place S. surinama into the three groups discussed above, but there are some potential conflicts regarding the placement of some S. septentrionalis specimens among the three gene fragments (S1 Fig.). The species tree generated by *BEAST was congruent with the concatenated BI and ML results with higher node support (PP: 1) for all ingroup clades (Fig. 4).

[Figure omitted. See PDF.]

Fig 4. Species tree generated by *BEAST.

Bayesian posterior probabilities indicated above or below branches. Branch colors follow Fig. 3.

https://doi.org/10.1371/journal.pone.0119151.g004

Divergence times and historical biogeography

The divergence dating analysis revealed a middle/late Miocene origin for Synoeca with subsequent diversification of extant species occurring in the Plio-Pleistocene (Fig. 5). The oldest divergence event was the split between the two major clades: (S. chalibea + S. virginea) and (S. cyanea + (S. septentrionalis MA + S. surinama/ S. septentrionalis AF)) at 7.87 million of years ago (mya) [95% of high posterior density (HPD): 11.16–5.44 mya]. The youngest diversification event was the separation of specimens of S. surinama and S. septentrionalis in the northern clade at 0.14 mya (95% of HPD: 0.02–0.31 mya). The biogeographic ancestral area analysis supports an Amazonian origin for Synoeca, with three independent dispersal-vicariance events from Amazonia to the Dry Diagonal and AF, one secondarily from Amazonia to Central America, and one from Amazonia to the AF (Fig. 6).

[Figure omitted. See PDF.]

Fig 5. Chronogram with divergence times generated by BEAST.

Divergence times between clades are annotated on each node and standard deviations [95% of high posterior density (HPD)] are represented by the blue bars.

https://doi.org/10.1371/journal.pone.0119151.g005

[Figure omitted. See PDF.]

Fig 6. Biogeographical analysis using Bayesian Binary MCMC, showing only the most likely states at each node.

https://doi.org/10.1371/journal.pone.0119151.g006

Discussion

Synoeca systematics

Despite the fact that Synoeca has only five species described, often some specimens are misidentified in collections. This problem was pointed out by Richards [9], for example concerning the misidentification of specimens of S. chalibea and S. virginea, which requires careful comparison. Indeed, the “yellow form” of S. chalibea is quite similar to S. virginea, but can be differentiated by the presence of punctation on the scutum and pronotum [20], as well as on the head. Also, specimens of S. septentrionalis from the AF show variation in the dark triangular area in the clypeus as verified by Menezes et al. [19], a major character used in the diagnosis of this species according to Cely and Sarmiento [21]. In some specimens, this dark area is totally absent, which may cause this species to be misidentified as S. cyanea, which also has a reddish clypeus. Considering this, specimens from AF identified as S. cyanea in collections need to be carefully verified. Despite the morphological similarities between S. septentrionalis AF and S. cyanea, the new morphological characters described here can be useful for correct species identification.

The division of Synoeca into two major clades is in disagreement with previous morphological phylogenetic analyses [20, 22]. However, morphological characters, such as punctation on the propodeum, clypeal-eye contact, malar space and wing color, might separate Synoeca into two groups (see [9, 20]) as shown here: (S. chalibea + S. virginea) + (S. cyanea + (S. surinama/S. septentrionalis)). Moreover, the previous placement of S. septentrionalis as sister to the clade (S. surinama + S. cyanea) is not supported by our molecular results. This clade in the morphological analyses is supported only by the absence of numerous erect outstanding setae on the first metasomal tergum and sternum (see [20]). However, Andena et al. [20] commented that this condition is homoplastic because it is found in several outgroups. Thus, this character may not be sufficient to establish relationships between the species within the clade S. cyanea + S. surinama/S. septentrionalis.

Historical biogeography

Our results also shed light on biogeographic patterns within Synoeca. Three Synoeca species occur in the AM (S. chalibea, S. virginea and S. surinama), S. septentrionalis occurs in the MA and AF, and only one species is restricted to the ESA, namely S. cyanea. Also, most of the species of the tribe to which Synoeca belongs (Epiponini) occurs in or is restricted to the AM [9]. Moreover, our analysis of ancestral area reconstruction supports an Amazonian origin of Synoeca (Fig. 6). Thus, the genus may have experienced three main colonization events from Amazonia during the Plio-Pleistocene.

The oldest inferred route probably occurred in southern South America between the AM and the AF. There is much evidence to support contact in the past between the AF and the AM [6, 7, 39–41]. Batalha-Filho et al. [7] combined phylogenetic and distributional data of avifauna and suggested old connections (middle to late Miocene) between AM and AF through the current southern Cerrado and Pantanal and the transition towards the Chaco and palm savannas of Bolivia and Paraguay. Also, it has been suggested that taxa representing divergences through this connection in the AF are spatially restricted to the southern AF and upland forests in southern Bahia, Minas Gerais, Espírito Santo, Rio de Janeiro and São Paulo [7]. Thus, the current geographical distribution of S. cyanea and its time of diversification at 3.46 mya (95% of HPD: 4.79–2.4 mya; Fig. 5) are consistent with this species colonizing by this route in southern South America.

A second route probably occurred relatively recently in direction to the Central America Rainforest region and southern North America via the Isthmus of Panama. The formation of Isthmus of Panama during the Pliocene at ~3.5 mya led to the Great American Biotic Interchange [2, 42]. The current distribution and diversification time of S. septentrionalis MA at 2.29 mya (95% of HPD: 3.15–1.56 mya; Fig. 5) agrees with this colonization route. Moreover, the lower species richness of Epiponini in the southern North America and Central America Rainforest regions compared to South America may be explained by recent dispersal via the Isthmus of Panama and, therefore, less time for species diversification.

The third route appears to have recently occurred in northeastern Brazil between AM and AF. Batalha-Filho et al. [7] suggest two connection pathways between AM and AF in northeastern Brazil, one through the coastal zones of Maranhão, Piauí, Ceará, and Rio Grande do Norte (Brazil), and another through Tocantins and Bahia (Brazil). Populations of S. surinama may have reached the AF by one of these routes and, after the breaking of this connection, experiencing consequent divergence between populations of AM/C and AF. Genetic differences among populations in both AM and AF are also found in the swarm-founding social wasp Angiopolybia pallens [43].

Atlantic Forest and S. septentrionalis

Despite the morphological similarity seen in S. septentrionalis specimens living in different regions, our molecular data suggest that S. septentrionalis and S. surinama specimens from the Atlantic Forest (AF) are interrelated and belong to a distinct lineage of Synoeca. Moreover, there appears to be a very recent division [at 0.38 mya (95% of HPD: 0.67–0.16 mya; Fig. 5)] between northern and central Atlantic forest groups (NAF and CAF). This division may be explained by the occurrence of refuges during the Pleistocene in the northern and central AF seen by Carnaval et al. [44] for amphibians. However, our molecular results may also be due to incomplete lineage sorting and/or mitochondrial introgression between S. septentrionalis AF and S. surinama AF. Nevertheless, there is insufficient morphological information to separate S. septentrionalis AF and S. septentrionalis MA and we suggest that they may represent two species with a very similar morphology (i.e., potential cryptic species). Further studies at the population level will be useful in characterizing the diversification processes occurring between members of S. surinama and S. septentrionalis.

Our molecular phylogenetic findings suggest that Synoeca species richness in the Neotropical Region may be underestimated due to morphological similarity and lack of broad geographical sampling. Further studies combining morphology, genetics, and population-level sampling should be seen as the main challenge in the future of phylogenetic research in social wasps as a whole. This work may result in the future recognition of additional species of social wasps in the Neotropics.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Trees constructed by Bayesian inference for each mitochondrial gene.

(a) 16S, (b) COI and (c) CytB.

https://doi.org/10.1371/journal.pone.0119151.s001

(TIF)

S1 Table. Target gene fragments, primer sequences, and origin of the primers used in this study.

https://doi.org/10.1371/journal.pone.0119151.s002

(DOCX)

Acknowledgments

We thank Alexandre Somavilla for samples from Amazonas and Maranhão, Elder F. Morato and Márcio Oliveira for logistic support in Acre and Manaus, respectively, Danilo Baia for samples from Amapá, Jaqueline Coimbra for samples from Distrito Federal, James M. Carpenter for samples from Argentina and Costa Rica, Lúcio Zardo for samples from Rio Grande do Sul, Gabriel A. R. Melo for samples from Mexico, Gisele G. Azevedo for samples from Maranhão, Floresta Nacional do Ibura (Nossa Senhora do Socorro, SE), Reserva Ecológica Carnijó (Moreno, PE) and Reserva Florestal Adolpho Ducke (Manaus, AM) for sampling permits in their areas. RSTM and SGB also thank Eugenia Okonski (NMNH) for helping with specimen imaging and the Laboratories of Analytical Biology (NMNH) for computational resources. We would like to thank Anthony Raw for his encouragement and suggestions for this work. We thank James M. Carpenter (AMNH, New York) for his suggestions on a previous version of the manuscript.

Author Contributions

Conceived and designed the experiments: RSTM SGB AFC MADL MAC. Performed the experiments: RSTM. Analyzed the data: RSTM SGB. Contributed reagents/materials/analysis tools: RSTM SGB AFC MADL MAC. Wrote the paper: RSTM SGB AFC MADL MAC. Collection of samples: RSTM AFC.

Citation: Menezes RST, Brady SG, Carvalho AF, Del Lama MA, Costa MA (2015) Molecular Phylogeny and Historical Biogeography of the Neotropical Swarm-Founding Social Wasp Genus Synoeca (Hymenoptera: Vespidae). PLoS ONE 10(3): e0119151. https://doi.org/10.1371/journal.pone.0119151

References

1. Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000; 403: 853–858. pmid:10706275

2. Hoorn C, Wesselingh FP, ter Steege H, Bermudez MA, Mora A, Sevink J, et al. Amazonia through time: andean uplift, climate change, landscape evolution, and biodiversity. Science. 2010; 330: 927–931. pmid:21071659

3. Hoorn C, Wesselingh F. Amazonia: landscape and species evolution: a look into the past. London: John Wiley & Sons; 2010.

4. Rull V. Speciation timing and neotropical biodiversity: the Tertiary—Quaternary debate in the light of molecular phylogenetic evidence. Mol Ecol. 2008; 17: 2722–2729. pmid:18494610

5. Rull V. Neotropical biodiversity: timing and potential drivers. Trends Ecol Evol. 2011; 26: 508–513. pmid:21703715

6. Costa LP. The historical bridge between the Amazon and the Atlantic forest of Brazil: a study of molecular phylogeography with small mammals. J Biogeogr. 2003; 30: 71–86.

7. Batalha-Filho H, Fjeldså J, Fabre PH, Miyaki CY. Connections between the Atlantic and the Amazonian forest avifaunas represent distinct historical events. J Ornithol. 2013; 154: 41–50.

8. Vanzolini PE. Ecological and geographical distribution of lizards in Pernambuco, north-eastern Brazil (Sauria). Pap Avulsos Zool. 1974; 28: 61–90.

9. Richards OW. The social wasps of the Americas (excluding the Vespinae). London: British Museum of Natural History; 1978.

10. Bermúdes EGC. Reprodução e dinâmica de populacão em Synoeca surinama (Hymenoptera: Vespidae). Acta Amazon. 1980a; 10: 679–690.

11. Bermúdes EGC. Orientação, arquitetura e construção dos ninhos de Synoeca surinama (L.) (Hymenoptera; Vespidae). Acta Amazon. 1980b; 10: 883–896.

12. Bermúdes EGC. Comportamento e causas da migração de Synoeca surinama L. (Hymenoptera, Vespidae), na amazonia brasileira. Rev Bras Entomol. 1981a; 26: 71–74.

13. Bermúdez EGC. Alarme e defesa no ninho de Synoeca surinama (Hymenoptera: Vespidae). Acta Amazon. 1981b; 11: 377–382.

14. Bermúdes EGC. Evidências sobre o comportamento de cópula dos machos de Synoeca surinama L. (Hymenoptera: Vespidae) num ninho em construção. Acta Amazon. 1982; 12: 665–666.

15. Elisei T, Ribeiro JC, Guimaraes DL, Prezoto F. Foraging activity and nesting of swarm-founding wasp Synoeca cyanea (Hymenoptera: Vespidae, Epiponini). Sociobiology. 2005; 46: 317–327.

16. Noda SCM, Shima SN, Noll FB. Morphological and physiological caste differences in Synoeca cyanea (Hymenoptera, Vespidae, Epiponini) according to the ontogenetic development of the colonies. Sociobiology. 2003; 41: 547–570.

17. Noll FB, Wenzel JW, Zucchi R. Evolution of Caste in Neotropical Swarm-Founding Wasps (Hymenoptera: Vespidae; Epiponini). Am Mus Novit. 2004; 3467: 1–24.

18. Wenzel JW. A generic key to the nests of hornets, yellowjackets, and paper wasps worldwide (Vespidae: Vespinae, Polistinae). Am Mus Novit. 1998; 3224: 1–39.

19. Menezes RST, Andena SR, Carvalho AF, Costa MA. First records of Synoeca septentrionalis Richards, 1978 (Hymenoptera, Vespidae, Epiponini) in the Brazilian Atlantic Rain Forest. Zookeys. 2011; 151: 75–78. pmid:22368453

20. Andena SR, Carpenter JM, Noll FB. A phylogenetic analysis of Synoeca de saussure, 1852, a neotropical genus of social wasps (Hymenoptera: Vespidae: Epiponini). Entomol Am. 2009; 115: 81–89.

21. Cely CC, Sarmiento CE. What about intraspecific variation? Reassessment of taxonomic and phylogenetic characters in the genus Synoeca de Saussure (Hymenoptera: Vespidae: Polistinae). Zootaxa. 2011; 2899: 43–59.

22. Carpenter JM, Andena SR, Noll FB, Wenzel JW. Well, what about intraspecific variation? Taxonomic and phylogenetic characters in the genus Synoeca de Saussure (Hymenoptera, Vespidae). Zootaxa. 2013; 3682: 421–431. pmid:25243297

23. Han H, McPheron BA. Molecular phylogenetic study of Tephritidae (Insecta: Diptera) using partial sequences of the Mitochondrial 16S Ribosomal DNA. Mol Phylogenet Evol. 1997; 7: 17–32. pmid:9007017

24. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32: 1792–97. pmid:15034147

25. Tamura K, Peterson D, Peterson D, Stecher G, Nei M, Kumar S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony method. Mol Biol Evol. 2011; 28: 2731–2739. pmid:21546353

26. Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012; 29: 1695–1701. pmid:22319168

27. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012; 61: 539–542. pmid:22357727

28. Zwickl DJ. Genetic Algorithm Approaches for the Phylogenetic Analysis of Large Biological Sequence Datasets under the Maximum Likelihood Criterion. Ph.D. Thesis, University of Texas at Austin, TX. 2006. Available: http://repositories.lib.utexas.edu/bitstream/handle/2152/2666/zwickld81846.pdf?sequence=2

29. Rambaut A, Drummond AJ. Tracer v.1.5; 2009. Available: http://beast.bio.ed.ac.uk/Tracer.

30. Rambaut A. FigTree: Tree Figure Drawing Tool v1.4.0; 2009. Available: http://tree.bio.ed.ac.uk/software/figtree.

31. Maddison WP, Knowles LL. Inferring phylogeny despite incomplete lineage sorting. Syst Biol. 2006; 55: 21–30. pmid:16507521

32. Knowles LL. Sampling strategies for species tree estimation. In: Knowles LL, Kubatko LS, editors. Estimating Species Trees: Practical and Theoretical Aspects. Wiley-Blackwell; 2010. pp. 163–173.

33. Heled J, Drummond AJ. Bayesian inference of species trees from multilocus data. Mol Biol Evol. 2010; 27: 570–580. pmid:19906793

34. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012; 29: 1969–1973. pmid:22367748

35. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. Plos Biol. 2006; 4, e88. pmid:16683862

36. Carpenter JM, Grimaldi DA. Social wasps in amber. Am Mus Novit. 1997; 3203.

37. Brower AVZ. Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from patterns of mitochondrial DNA Evolution. Proc Natl Acad Sci USA. 1994; 91: 6491–6495. pmid:8022810

38. Yu Y, Harris AJ, He XJ. RASP (Reconstruct Ancestral State in Phylogenies) 2.1b; 2014. Available: http://mnh.scu.edu.cn/soft/blog/RASP.

39. Oliveira PE, Barreto AMF, Suguio K. Late pleistocene/holocene climatic and vegetational history of the Brazilian caatinga: the fossil dunes of the middle São Francisco River. Palaeogeogr Palaeoclimatol Palaeoecol. 1999; 152: 319–337.

40. Behling H, Arz HW, Patzold J, Wefer G. Late Quaternary vegetational and climate dynamics in northeastern Brazil, inferences from marine core GeoB3104–1. Quat Sci Rev. 2000; 19: 981–994.

41. Auler AS, Wang A, Edwards RL, Cheng H, Cristalli PS, Smart ML, et al. Quaternary ecological and geomorphic changes associated with rainfall events in presently semi-arid northeastern Brazil. J Quat Sci. 2004; 19: 693–701.

42. Coates AG, Jackson JBC, Collins LS, Cronin TM, Dowsett HJ, Bybell LM, et al. Closure of the Isthmus of Panama the near-shore marine record of Costa Rica and western Panama. Geol Soc Am Bull. 1992; 104: 814–828.

43. Carvalho AF, Santos GMM, Menezes RST, Costa MA. Genetic diversity of Angiopolybia pallens (Lepeletier) (Hymenoptera, Vespidae, Polistinae) explained by the disjunction of South American rainforests. Genet Mol Res. 2014; 13: 89–94. pmid:24446290

44. Carnaval AC, Hickerson MJ, Haddad CFB, Rodrigues MT, Moritz C. Stability predicts genetic diversity in the Brazilian Atlantic Forest Hotspot. Science. 2009; 323: 785–789. pmid:19197066