Morphometric and genetic differentiation in isolated populations of the endangered Mesoamerican stingless bee Melipona yucatanica (Hymenoptera: Apoidea) suggest the existence of a two species complex

Phylogenetic and Morphological Characteristics Reveal Cryptic Speciation in Stingless Bee, Tetragonula laeviceps s.l. Smith 1857 (Hymenoptera; Meliponinae)

Tetragonula laeviceps sensu lato (s.l.) Smith 1857 has the most complicated nomenclatural history among the Tetragonula genera. The objective of this study was to investigate whether T. laeviceps s.l. individuals with worker bees are grouped in the same or nearly the same morphological characteristics and have similar COI haplotype cluster groups. A total of 147 worker bees of T. laeviceps s.l. were collected from six sampling sites in Sabah (RDC, Tuaran, Kota Marudu, Putatan, Kinarut and Faculty of Sustainable Agriculture (FSA)), but only 36 were selected for further studies. These specimens were first classified according to the most obvious morphological characteristics, i.e., hind tibia color, hind basitarsus color and body size. Group identification was based on morphological characteristics important for distinguishing the four groups within T. laeviceps s.l. The four groups of T. laeviceps s.l. had significantly different body trait measurements for the TL (total length), HW (head width), HL (head length), CEL (compound eye length), CEW (compound eye width), FWLT (forewing length, including tegula), FWW (forewing width), FWL (forewing length), ML (mesoscutum length), MW (mesoscutum width), SW (mesoscutellum width), SL (mesoscutellum length), HTL = (hind tibia length), HTW (hind tibia width), HBL (hind basitarsus length) and HBW (hind basitarsus width) (p < 0.001). Body color included HC (head color), CC (clypeus color), ASC (antennae scape color), CFPP (Clypeus and frons plumose pubescence), HTC (hind tibia color), BSC (basitarsus color), SP (leg setae pubescence), SP (Thorax mesoscutellum pubescence), SPL (thorax mesoscutellum pubescence length) and TC (thorax color) (p < 0.05). The most distinctive features of the morphological and morphometric characteristics measured by PCA and LDA biplot that distinguish Group 1 (TL6-1, TL6-2 and TL6-3) from the other groups were the yellowish-brown ASC and the dark brown TC. Group 2 (haplotypes TL2-1, TL2-2 and TL2-3 and TL4-1, TL4-2 and TL4-3) had a dark brown ASC and a black TC, while Group 3 (haplotypes TL11-1, TL11-2 and TL11-3) had a blackish-brown ASC, a black TC and the largest TL, FWW and FWL. As for phylogenetic relationships, 12 out of 36 haplotypes showed clear separation with good bootstrap values (97-100%). The rest of the haplotypes did not show clear differentiation between subclades that belonged together, regardless of their morphology and morphometric characteristics. This suggests that the combination of DNA barcoding for species identification and phylogenetic analysis, as well as traditional methods based on morphological grouping by body size and body color, can be reliably used to determine intraspecific variations within T. laeviceps s.l.

Wing geometry and genetic analyses reveal contrasting spatial structures between male and female Aedes aegypti (L.) (Diptera: Culicidae) populations in metropolitan Manila, Philippines

The population genetic structure of Aedes aegypti (Linnaeus, 1762) has been studied in order to understand its role as an efficient vector. Several studies utilized an integrative approach; to combine genetic and phenotypic data to determine its population structure but these studies have only focused on female populations. To address this particular gap, our study compared the population variability and structuring between its male and female populations using phenotypic and genetic data from a highly-urbanized and dengue-endemic region of the Philippines, Metropolitan Manila. Five mosquito populations comprised of female (n = 137) and male (n = 49) adult mosquitoes were used in this study. All mosquito individuals underwent geometric morphometric (26 landmarks), and genetic (11 microsatellite loci) analyses. Results revealed that F_ST estimates (genetic) were 0.055 and 0.009 while Q_ST estimates (phenotypic) were 0.318 and 0.309 in in male and female populations, respectively. Wing shape variation plots showed that male populations were distinctly separated from each other while female populations overlapped. Similarly, discriminant analysis of principal components using genetic data revealed that male populations were also distinctly separated from each other while female populations showed near-overlapping populations. Genetic and phenetic dendrograms showed the formation of two groups in male populations but no groups in female populations. Further analysis indicated a significant correlation (r = 0.68, p = 0.02) between the genetic and phenetic distances of male populations. Bayesian analysis using genetic data also detected multiple clusters in male (K = 3) and female (K = 2) populations, while no clusters were detected using the phenotypic data from both sexes. Our results revealed contrasting phenotypic and genetic patterns between male and female Ae. aegypti, indicating that male populations were more spatially structured than female populations.

Microsatellite Marker Discovery in the Stingless Bee Uruçu-Amarela (Melipona rufiventris Group, Hymenoptera, Meliponini) for Population Genetic Analysis

The species Melipona rufiventris Lepeletier, 1836 is a Brazilian native stingless bee that is part of a species complex known as the ‘rufiventris group’, making it difficult to distinguish between the different species. Populations in this group are facing a severe decline, leading to the risk of local extinction, and therefore, their conservation should be treated as a major concern. This study describes the first set of tri- and tetranucleotide microsatellite markers, using next-generation sequencing technology for use in the identification of genetic diversity and population structure in the ‘rufiventris group’. A total of 16 microsatellite loci displayed polymorphism. Analysis of the whole data set (n = 50) detected 63 alleles in all loci, ranging from 2 to 7 with a mean of 3.9 alleles/locus. A genetic diversity analysis revealed high values for population differentiation estimates (F_ST = 0.252, R_ST = 0.317, and D_EST = 0.284) between the Atlantic Forest, Cerrado, and Caatinga biomes. An additional evidence for genetic divergence among populations was also found in the ’rufiventris group’; these should be treated as separate conservation units or even as separate species. These microsatellite markers have demonstrated a strong potential for assessing population discrimination in this threatened stingless bee group.

Patterns of sexual size dimorphism in stingless bees: Testing Rensch's rule and potential causes in highly eusocial bees (Hymenoptera: Apidae, Meliponini)

Eusocial insects offer a unique opportunity to analyze the evolution of body size differences between sexes in relation to social environment. The workers, being sterile females, are not subject to selection for reproductive function providing a natural control for parsing the effects of selection on reproductive function (i.e., sexual and fecundity selection) from other kinds of natural selection. Patterns of sexual size dimorphism (SSD) and testing of Rensch's rule controlling for phylogenetic effects were analyzed in the Meliponini or stingless bees. Theory predicts that queens may exhibit higher selection for fecundity in eusocial taxa, but contrary to this, we found mixed patterns of SSD in Meliponini. Non-Melipona species generally have a female-biased SSD, while all analyzed species of Melipona showed a male-biased SSD, indicating that the direction and magnitude of the selective pressures do not operate in the same way for all members of this taxon. The phylogenetic regressions revealed that the rate of divergence has not differed between the two castes of females and the males, that is, stingless bees do not seem to follow Rensch's rule (a slope >1), adding this highly eusocial taxon to the various solitary insect taxa not conforming with it. Noteworthy, when Melipona was removed from the analysis, the phylogenetic regressions for the thorax width of males on queens had a slope significantly smaller than 1, suggesting that the evolutionary divergence has been larger in queens than males, and could be explained by stronger selection on female fecundity only in non-Melipona species. Our results in the stingless bees question the classical explanation of female-biased SSD via fecundity and provide a first evidence of a more complex determination of SSD in highly eusocial species. We suggest that in highly eusocial taxa, additional selection mechanisms, possibly related to individual and colonial interests, could influence the evolution of environmentally determined traits such as body size.

Molecular Characterization of the Indigenous Stingless Bees (Tetragonula spp. Complex) Using ISSR Marker from Southern Peninsular India

India is a country bestowed enormously with stingless bees, but genetic information about them is extremely minimal. This study focused to tap the geographic allocation, genetic variability, and differentiation among Tetragonula species complexes from natural and semi-urban habitats. Genetic analyses were assessed among 36 contrasting genotypes utilizing 20 ISSR primers. The dual combination exquisitely and productively amplified 245 DNA fragments at the loci, of which 240 bands were polymorphic (97.95%). Low to moderate level of genetic differentiation was detected from different estimators (Ht 0.29, G' _STest 0.16, D _est 0.072, F _ST 0.14, and Nm 2.68). Hierarchical clustering analysis aided to partition the individual genotypes into its respective five species group formed, aided by substantial bootstrap support values, but differing under morphological identification. It also provided valuable insight into the moderate eco-genetic diversity (H 0.39) prevailing from geographically scattered inhabitants. Potential exploitation of hyper-variable ISSR marker turned out fairly as a promising technique for finding valid polymorphisms and infers relevant variations. This baseline information enhances our understanding of the genetic status of the indigenous species from the country.

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population-specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.

Vulnerable Habitats Alter African Meliponine Bee's (Hymenoptera: Apidae) Assemblages in an Eastern Afromontane Biodiversity Hotspot

Habitat degradation has over time formed synergy with other factors to contribute to dwindling populations of both fauna and flora by altering their habitats. The disturbance of natural habitats affects the diversity of both vertebrates and invertebrates by altering both feeding and nesting sites for which organisms are known to depend on for survival. Little is known of the extent to which vulnerable habitats could shape the diversity of most indigent pollinators such as African meliponine bee species in tropical ecosystems. This study was conducted to determine how disturbance could shape the natural occurrence of African meliponine bee species in different ecological habitats of Taita Hills, leading to changes in their diversity. A total of four species depicted by the Renyi diversity profile was recorded in five of the six main habitat types surveyed, and a further extrapolation with Shannon index (EH) also predicted the highest species richness of 4.24 in a deciduous habitat type. These meliponine bee species (Hypotrigona gribodoi, Hypotrigona ruspolii, Meliponula ferruginea (black), and Plebeina hildebrandti) were observed to be unevenly distributed across all habitats, further indicating that mixed deciduous habitat was more diverse than acacia-dominated bush lands, grasslands, and exotic forest patches. Geometric morphometrics categorized all four meliponine bee species into two major clusters-cluster 1 (H gribodoi, H ruspolii, M ferruginea (black)) and cluster 2 (P hildebrandti)-and further discriminated populations against the 4 potential habitats they are likely to persist or survive in. Each habitat appeared to consist of a cluster of subpopulations and may possibly reveal ecotypes within the four meliponine populations. This has revealed that unprecedented conversions of natural habitats to agroecosystems are a key driving factor causing increased habitat isolation and vulnerability in this Afromontane region which may potentially distort local assemblages of native pollinators, such as meliponine bee species.

Distribution and predictors of wing shape and size variability in three sister species of solitary bees

Morphological traits can be highly variable over time in a particular geographical area. Different selective pressures shape those traits, which is crucial in evolutionary biology. Among these traits, insect wing morphometry has already been widely used to describe phenotypic variability at the inter-specific level. On the contrary, fewer studies have focused on intra-specific wing morphometric variability. Yet, such investigations are relevant to study potential convergences of variation that could highlight micro-evolutionary processes. The recent sampling and sequencing of three solitary bees of the genus Melitta across their entire species range provides an excellent opportunity to jointly analyse genetic and morphometric variability. In the present study, we first aim to analyse the spatial distribution of the wing shape and centroid size (used as a proxy for body size) variability. Secondly, we aim to test different potential predictors of this variability at both the intra- and inter-population levels, which includes genetic variability, but also geographic locations and distances, elevation, annual mean temperature and precipitation. The comparison of spatial distribution of intra-population morphometric diversity does not reveal any convergent pattern between species, thus undermining the assumption of a potential local and selective adaptation at the population level. Regarding intra-specific wing shape differentiation, our results reveal that some tested predictors, such as geographic and genetic distances, are associated with a significant correlation for some species. However, none of these predictors are systematically identified for the three species as an important factor that could explain the intra-specific morphometric variability. As a conclusion, for the three solitary bee species and at the scale of this study, our results clearly tend to discard the assumption of the existence of a common pattern of intra-specific signal/structure within the intra-specific wing shape and body size variability.

Beekeeping practices and geographic distance, not land use, drive gene flow across tropical bees

Across the globe, wild bees are threatened by ongoing natural habitat loss, risking the maintenance of plant biodiversity and agricultural production. Despite the ecological and economic importance of wild bees and the fact that several species are now managed for pollination services worldwide, little is known about how land use and beekeeping practices jointly influence gene flow. Using stingless bees as a model system, containing wild and managed species that are presumed to be particularly susceptible to habitat degradation, here we examine the main drivers of tropical bee gene flow. We employ a novel landscape genetic approach to analyse data from 135 populations of 17 stingless bee species distributed across diverse tropical biomes within the Americas. Our work has important methodological implications, as we illustrate how a maximum-likelihood approach can be applied in a meta-analysis framework to account for multiple factors, and weight estimates by sample size. In contrast to previously held beliefs, gene flow was not related to body size or deforestation, and isolation by geographic distance (IBD) was significantly affected by management, with managed species exhibiting a weaker IBD than wild ones. Our study thus reveals the critical importance of beekeeping practices in shaping the patterns of genetic differentiation across bee species. Additionally, our results show that many stingless bee species maintain high gene flow across heterogeneous landscapes. We suggest that future efforts to preserve wild tropical bees should focus on regulating beekeeping practices to maintain natural gene flow and enhancing pollinator-friendly habitats, prioritizing species showing a limited dispersal ability.

Genetic variability in captive populations of the stingless bee Tetragonisca angustula

Low genetic variability has normally been considered a consequence of animal husbandry and a major contributing factor to declining bee populations. Here, we performed a molecular analysis of captive and wild populations of the stingless bee Tetragonisca angustula, one of the most commonly kept species across South America. Microsatellite analyses showed similar genetic variability between wild and captive populations However, captive populations showed lower mitochondrial genetic variability. Male-mediated gene flow, transport and division of nests are suggested as the most probable explanations for the observed patterns of genetic structure. We conclude that increasing the number of colonies kept through nest divisions does not negatively affect nuclear genetic variability, which seems to be maintained by small-scale male dispersal and human-mediated nest transport. However, the transport of nests from distant localities should be practiced with caution given the high genetic differentiation observed between samples from western and eastern areas. The high genetic structure verified is the result of a long-term evolutionary process, and bees from distant localities may represent unique evolutionary lineages.

The integrative taxonomic approach reveals host specific species in an encyrtid parasitoid species complex

Integrated taxonomy uses evidence from a number of different character types to delimit species and other natural groupings. While this approach has been advocated recently, and should be of particular utility in the case of diminutive insect parasitoids, there are relatively few examples of its application in these taxa. Here, we use an integrated framework to delimit independent lineages in Encyrtus sasakii (Hymenoptera: Chalcidoidea: Encyrtidae), a parasitoid morphospecies previously considered a host generalist. Sequence variation at the DNA barcode (cytochrome c oxidase I, COI) and nuclear 28S rDNA loci were compared to morphometric recordings and mating compatibility tests, among samples of this species complex collected from its four scale insect hosts, covering a broad geographic range of northern and central China. Our results reveal that Encyrtus sasakii comprises three lineages that, while sharing a similar morphology, are highly divergent at the molecular level. At the barcode locus, the median K2P molecular distance between individuals from three primary populations was found to be 11.3%, well outside the divergence usually observed between Chalcidoidea conspecifics (0.5%). Corroborative evidence that the genetic lineages represent independent species was found from mating tests, where compatibility was observed only within populations, and morphometric analysis, which found that despite apparent morphological homogeneity, populations clustered according to forewing shape. The independent lineages defined by the integrated analysis correspond to the three scale insect hosts, suggesting the presence of host specific cryptic species. The finding of hidden host specificity in this species complex demonstrates the critical role that DNA barcoding will increasingly play in revealing hidden biodiversity in taxa that present difficulties for traditional taxonomic approaches.