Phylogenetic analysis of the mitochondrial genomes in bees (Hymenoptera: Apoidea: Anthophila)

The Mitochondrial Genome of Melipona fasciculata (Apidae, Meliponini): Genome Organization and Comparative Analyses, Phylogenetic Implications and Divergence Time Estimations

The native stingless bee Melipona fasciculata is economically and ecologically important to the Brazilian Northeast, providing a sustainable source of income to family farmers and being considered an effective pollinator in most ecosystems and crops. This study describes, for the first time, the mitogenome of the species and its phylogenetic position. The mitochondrial genome was sequenced using a MiSeq Sequencer (Illumina Inc.) and compared with other GenBank bee mitogenomes. The length of the mitochondrial DNA, excluding most of the control region, is 14,753 bp, and contains 13 protein-coding genes (PCGs), 21 transfer RNAs, 2 ribosomal RNAs (16S and 12S), and 1 AT-rich region. The GC-content of the M. fasciculata mitogenome was 13.4%. Of the 36 coding regions, 12 tRNAs and 9 PCGs were encoded in the heavy strand, and 9 tRNAs, 4 PCGs and 2 rRNAs were encoded in the light strand. The relative orientation and gene order was the same as other stingless bee mitogenomes. Phylogenetic inference produced well-resolved relationships with high statistical support for concordant branch topologies, under different optimization schemes and model parameters, within and among Melipona, Bombus, Apis, and related clades of Hymenoptera. In general, our divergence time estimates, which were based on the concatenated gene sequences (PCGs + rRNAs) from various groups, overlapped estimations captured by Bayesian analysis from different studies. The divergence time among Melipona species was estimated to occur during the Oligocene, approximately 24 Mya (95% HPD 14-36 Mya). Our results represent a valuable addition to help understanding not only the taxonomy and evolution of Brazilian stingless bee species, but also to uncover historical dispersal and isolation patterns in Meliponinae.

Comparative analyses of mitogenomes in the social bees with insights into evolution of long inverted repeats in the Meliponini

The insect mitogenome is typically a compact circular molecule with highly conserved gene contents. Nonetheless, mitogenome structural variations have been reported in specific taxa, and gene rearrangements, usually the tRNAs, occur in different lineages. Because synapomorphies of mitogenome organizations can provide information for phylogenetic inferences, comparative analyses of mitogenomes have been given increasing attention. However, most studies use a very few species to represent the whole genus, tribe, family, or even order, overlooking potential variations at lower taxonomic levels, which might lead to some incorrect inferences. To provide new insights into mitogenome organizations and their implications for phylogenetic inference, this study conducted comparative analyses for mitogenomes of three social bee tribes (Meliponini, Bombini, and Apini) based on the phylogenetic framework with denser taxonomic sampling at the species and population levels. Comparative analyses revealed that mitogenomes of Apini and Bombini are the typical type, while those of Meliponini show diverse variations in mitogenome sizes and organizations. Large inverted repeats (IRs) cause significant gene rearrangements of protein coding genes (PCGs) and rRNAs in Indo-Malay/Australian stingless bee species. Molecular evolution analyses showed that the lineage with IRs have lower d N/ d S ratios for PCGs than lineages without IRs, indicating potential effects of IRs on the evolution of mitochondrial genes. The finding of IRs and different patterns of gene rearrangements suggested that Meliponini is a hotspot in mitogenome evolution. Unlike conserved PCGs and rRNAs whose rearrangements were found only in the mentioned lineages within Meliponini, tRNA rearrangements are common across all three tribes of social bees, and are significant even at the species level, indicating that comprehensive sampling is needed to fully understand the patterns of tRNA rearrangements, and their implications for phylogenetic inference.

Possible Epigenetic Origin of a Recurrent Gynandromorph Pattern in Megachile Wild Bees

Simple Summary

Gynandromorphs, i.e., individuals with a mix of male and female body parts, are known for many species of insects and other animals with separate sexes. This anomaly is generally regarded as the result of localized genetic mutations in sex-determining genes. We analyzed the specific mix of male and female characters in naturally occurring gynandromorphs of 21 species of the wild bee genus Megachile and found a recurrent pattern. Based on the regularity of this pattern, and the current knowledge on sex determination and sex differentiation in the relatively closely-related honey bee, we argue that the origin of these composite phenotypes is possibly epigenetic, rather than genetic, i.e., produced by some defects in the maintenance of the regulatory signals that control sex differentiation at the level of single cell lineages, rather than triggered by genetic mutations.

Abstract

Gynandromorphs, i.e., individuals with a mix of male and female traits, are common in the wild bees of the genus Megachile (Hymenoptera, Apoidea). We described new transverse gynandromorphs in Megachile pilidens Alfkeen, 1924 and analyze the spatial distribution of body parts with male vs. female phenotype hitherto recorded in the transverse gynandromorphs of the genus Megachile. We identified 10 different arrangements, nine of which are minor variants of a very general pattern, with a combination of male and female traits largely shared by the gynandromorphs recorded in 20 out of 21 Megachile species in our dataset. Based on the recurrence of the same gynandromorph pattern, the current knowledge on sex determination and sex differentiation in the honey bee, and the results of recent gene-knockdown experiments in these insects, we suggest that these composite phenotypes are possibly epigenetic, rather than genetic, mosaics, with individual body parts of either male or female phenotype according to the locally expressed product of the alternative splicing of sex-determining gene transcripts.

Molecular phylogeny, historical biogeography and revised classification of andrenine bees (Hymenoptera: Andrenidae)

The mining bee subfamily Andreninae (Hymenoptera: Andrenidae) is a widely distributed and diverse group of ground-nesting solitary bees, including numerous species known to be important pollinators. Most of the species diversity of Andreninae is concentrated in the mainly Holarctic genus Andrena, comprising ca. 1550 described species. The subfamily and especially the genus have remained relatively neglected by recent molecular phylogenetic studies, with current classifications relying largely on morphological characters. We sampled ultraconserved element (UCE) sequences from 235 taxa, including all andrenine genera and 98 out of 104 currently recognized Andrena subgenera. Using 419,858 aligned nucleotide sites from 1009 UCE loci, we present a comprehensive molecular phylogenetic analysis of the subfamily. Our analysis supports the recognition of seven distinct genera in the Andreninae: Alocandrena, Ancylandrena, Andrena, Cubiandrena, Euherbstia, Megandrena, and Orphana. Within the genus Andrena, present-day subgeneric concepts revealed high degrees of paraphyly and polyphyly, due to strong homoplasy of morphological characters, necessitating a thorough, extensive revision of the higher classification of the genus. Based on our findings, we place the subgenus Calcarandrena in synonymy with Andrena (Lepidandrena); Hyperandrena, Nemandrena, Scoliandrena, Tylandrena and Zonandrena with A. (Melandrena); Distandrena, Fumandrena and Proxiandrena with A. (Micrandrena); Carandrena with A. (Notandrena); Agandrena with A. (Plastandrena); Xiphandrena with A. (Scrapteropsis); and Platygalandrena and Poliandrena with A. (Ulandrena) (new synonymies). We additionally reestablish the groups known as Opandrena and Truncandrena as valid subgenera of Andrena. Our results also show that the MRCA of Andrena+Cubiandrena dispersed from the New World to the Palaearctic probably during the Eocene-early Oligocene, followed by 10-14 Neogene dispersal events from the Palaearctic to the Nearctic and 1-6 Neogene dispersals back into the Palaearctic, all within the genus Andrena.

Comparative Mitogenomic Analysis of Two Cuckoo Bees (Apoidea: Anthophila: Megachilidae) with Phylogenetic Implications

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.

Global systematic diversity, range distributions, conservation and taxonomic assessments of graylings (Teleostei: Salmonidae; Thymallus spp.)

Graylings (Thymallus) are among the less well-studied groups of salmonid fishes, especially across their Asian distribution range. Here we perform a comprehensive global review of their phylogeography, systematic diversity and range distributions, including biogeographic reconstruction and assessment of both conservation and taxonomic status of each species. Based on a mitogenomic phylogenetic analysis, three approaches to the delineation of molecular operational units, and evaluation of 15 a-priori defined species, we provide biological support for the recognition of 13 grayling species, plus two additional species tentatively. Several instances of paraphyly and its potential effect on systematic inferences are discussed. Overall, the genus displays increasing species diversity and decreasing range size from higher to lower latitudes and ancestral trait reconstruction supports an East Asian origin for extant diversity, most likely centred in the Amur River drainage. Europe’s colonization by Thymallus took place as early as the late Miocene, at least two colonisations of North America are supported, and multiple dispersal events likely took place into Western Siberia. The conservation status for the 15 taxa was estimated to be: 6 least concern, 1 near-threatened, 2 vulnerable, 3 endangered and 3 data deficient.

The complete mitogenome of Lasioglossum affine (Hymenoptera: Halictidae) and phylogenetic analysis

The complete mitogenome of Lasioglossum affine (Hymenoptera: Halictidae) was sequenced and analyzed. The whole mitogenome is 17,352 bp (AT%=84.1%) and encodes 37 typical eukaryotic mitochondrial genes, including 13 protein-coding genes (PCGs), 22 tRNAs, two rRNAs, and an AT-rich region. Further analysis found three gene rearrangements, where trn I-Q-M → trn M-I-Q, trn W-C-Y → trn C-W-Y, and trn K-D → trn D-K were shuffled. The phylogenetic relationships of 19 species of Hymenoptera were established using maximum-likelihood method based on 13 concatenated PCGs. The result showed that Lasioglossum affine is a sister of Lasioglossum sp. SJW-2017.

The complete mitogenome of Nomia chalybeata (Hymenoptera: Halictidae) and phylogenetic analysis

The complete mitogenome of Nomia chalybeata was sequenced. The mitochondrial length of N. chalybeata was 16,692 bp (AT content 85.4%), with 37 classic invertebrate mitochondrial genes (including 13 protein-coding genes, 22 transporter RNAs, and two ribosomal RNAs) and AT-rich region (AT content 91.7%). The maximum-likelihood phylogenetic relationship was constructed using 11 species from Hymenoptera. Through the phylogenetic relationship, our research team successfully used the molecular data of the mitochondrial genome to verify that N. chalybeata belongs to the family Halictidae, and also provides molecular data for the database of the family Halictidae.

The complete mitogenome of Habropoda rodoszkowskii (Hymenoptera: Apidae) and phylogenetic analysis

The mitogenome of Habropoda rodoszkowskii, the first complete mitogenome sequence of the genus Habropoda (hymenoptera: Apidae), was sequenced. The mitogenome is 18,497 bp (The proportion of A + T in 80.7%) long, with 37 classic eukaryotic mitochondrial genes (including 13 PCGs, 22 tRNAs, and 2 rRNAs) and an AT-rich region (The proportion of A + T in 78.2%). The Bayesian-inference and Maximum-likelihood phylogenetic relationship was constructed using 15 species from Hymenoptera. According to the phylogenetic tree, Habropoda rodoszkowskii converges with genus Nomada bees (Nomada flava and Nomada flavoguttata) to be supported. In addition, Habropoda rodoszkowskii is more closely related to Apidae than to Megahilidae and Colletidae.