The South America encompasses the highest levels of biodiversity found anywhere in the world and its rich biota is distributed among many different biogeographical regions. However, many regions of South America are still poorly studied, including its xeric environments, such as the threatened Caatinga and Cerrado phytogeographical domains. In particular, the effects of Quaternary climatic events on the demography of endemic species from xeric habitats are poorly understood. The present study uses an integrative approach to reconstruct the evolutionary history of Partamona rustica, an endemic stingless bee from dry forest diagonal in Brazil, in a spatial-temporal framework. In this sense, we sequenced four mitochondrial genes and genotyped eight microsatellite loci. Our results identified two population groups: one to the west and the other to the east of the São Francisco River Valley (SFRV). These groups split in the late Pleistocene, and the Approximate Bayesian Computation approach and phylogenetic reconstruction indicated that P. rustica originated in the west of the SFRV, subsequently colonising eastern region. Our tests of migration detected reduced gene flow between these groups. Finally, our results also indicated that the inferences both from the genetic data analyses and from the spatial distribution modelling are compatible with historical demographic stability.

Background Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste. Results The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication. Conclusions The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Gene variation and the differentiation of two populations of Partamona seridoensis (Hymenoptera: Apidae: Meliponini) from the Caatinga biome, a semiarid ecosystem unique to Brazil, were estimated through allozymic and microsatellite analyses. These populations exhibited similar low degrees of enzyme gene variation. Observed genotype frequencies at the allozyme and microsatellite loci were in accordance with Hardy-Weinberg equilibrium in the two populations. Both markers demonstrated that the two populations are not genetically homogeneous and must be considered distinct populations. The occurrence of private alleles at the allozyme and microsatellite loci corroborates this differentiation, sustaining the hypothesis of a low level of interpopulation gene flow. The phenotypic segregations clearly demonstrated that the progeny inside each nest were the result of mating between the queen of the colony and only one male.

Partamona seridoensis is an endemic stingless bee from the Caatinga, a Neotropical dry forest in northeastern Brazil. Like other stingless bees, this species plays an important ecological role as a pollinator. The aim of the present study was to investigate the genetic structure and evolutionary history of P. seridoensis across its current geographic range. Workers from 84 nests from 17 localities were analyzed for COI and Cytb genic regions. The population structure tests (Bayesian phylogenetic inference, AMOVA and haplotype network) consistently characterized two haplogroups (northwestern and eastern), with little gene flow between them, generating a high differentiation between them as well as among the populations within each haplogroup. The Mantel test revealed no isolation by distance. No evidence of a potential geographic barrier in the present that could explain the diversification between the P. seridoensis haplogroups was found. However, Pleistocene climatic changes may explain this differentiation, since the initial time for the P. seridoensis lineages diversification took place during the mid-Pleistocene, specifically the interglacial period, when the biota is presumed to have been more associated with dry conditions and had more restricted, fragmented geographical distribution. This event may have driven diversification by isolating the two haplogroups. Otherwise, the climatic changes in the late Pleistocene must not have drastically affected the population dynamics of P. seridoensis, since the Bayesian Skyline Plot did not reveal any substantial fluctuation in effective population size in either haplogroup. Considering its importance and the fact that it is an endemic bee from a very threatened Neotropical dry forest, the results herein could be useful to the development of conservation strategies for P. seridoensis.

No abstract