A protocol for isolating insect mitochondrial genomes: a case study of NUMT in Melipona flavolineata (Hymenoptera: Apidae)

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

Insights from Melipona bicolor hybrid genome assembly: a stingless bee genome with chromosome-level scaffold

BACKGROUND

The highly eusocial stingless bees are crucial pollinators of native and agricultural ecosystems. Nevertheless, genomic studies within this bee tribe remain scarce. We present the genome assembly of the stingless bee Melipona bicolor. This bee is a remarkable exception to the typical single-queen colony structure, since in this species, multiple queens may coexist and share reproductive duties, resulting in genetically diverse colonies with weak kinship connections. As the only known genuinely polygynous bee, M. bicolor's genome provides a valuable resource for investigating sociality beyond kin selection.

RESULTS

The genome was assembled employing a hybrid approach combining short and long reads, resulting in 241 contigs spanning 259 Mb (N50 of 6.2 Mb and 97.5% complete BUSCOs). Comparative analyses shed light on some evolutionary aspects of stingless bee genomics, including multiple chromosomal rearrangements in Melipona. Additionally, we explored the evolution of venom genes in M. bicolor and other stingless bees, revealing that, apart from two genes, the conserved repertoire of venom components remains under purifying selection in this clade.

CONCLUSION

This study advances our understanding of stingless bee genomics, contributing to the conservation efforts of these vital pollinators and offering insights into the evolutionary mechanisms driving their unique adaptations.

Evidence for morphological and genetic structuring of Plebeia flavocincta (Apidae: Meliponini) populations in Northeast Brazil

Tropical dry forests are increasingly threatened by human activities. In Northeast Brazil, the Caatinga, an area of tropical dry forests surrounded by wetlands, is considered globally unique among these habitats. In this region, the stingless bee Plebeia flavocincta is found in a variety of environmental, ecological, and demographic conditions. We aimed to characterize P. flavocincta populations within its natural range through wing geometric morphometrics and mitochondrial DNA (mtDNA) analyses. The characterization of population variability can clarify whether the species is morphologically and genetically diverse and whether populations are morphologically and genetically structured. We analyzed 673 samples by wing morphometry and 75 by cytochrome-b assays. Our results revealed P. flavocincta is genetically and morphologically diverse and populations are morphologically and genetically structured. Despite the differentiation between the two most morphologically distant populations, we verified a large overlap of morphological variation between all populations. The genetic analysis showed that the haplotypes were geographically structured into six clusters, four of which were located in coastal areas, and the remaining two in the drier inland region. The characterization of P. flavocincta populations is an important step toward decision-making in programs for the protection, management, and sustainable use of this stingless bee in local breeding efforts.

Five Mitochondrial Genomes of the Genus Eysarcoris Hahn, 1834 with Phylogenetic Implications for the Pentatominae (Hemiptera: Pentatomidae)

Simple Summary

Pentatominae is the largest subfamily in the Pentatomidae, and most of its species are considered important agricultural pests. The phylogenetic relationships of tribes within Pentatominae remain controversial despite the fact that many studies have been performed using various molecular markers. In this study, five mitogenomes of the genus Eysarcoris were sequenced and analyzed, and the phylogenetic relationships of tribes within Pentatominae were reconstructed. The gene arrangement of the five mitochondrial genomes were found to be conserved and identical to other heteropteran mitogenomes. Differences in start codon usage and tandem repeats within control regions were found between E. gibbosus and the other four Eysarcoris species. In addition, the phylogenetic analyses showed that E. gibbosus is the first diverging clade within Eysarcorini. The results support the proposal to transfer E. gibbosus to the Stagonomus, another genus of Eysarcorini. Our results clarified the phylogenetic relationships among tribes of Pentatominae and laid a foundation for the further studies of Pentatominae.

Abstract

Four complete mitogenomes of Eysarcoris rosaceus, E. montivagus, E. gibbosus, E. annamita and one near-complete mitochondrial genome of E. ventralis were sequenced and used to explore the phylogenetic relationships of tribes within the subfamily Pentatominae. The mitogenomes range from 15,422 to 16,043 base pairs (bp) in length and encode 37 genes, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes (21 in E. ventralis), and a control region. Similar to other heteropteran species, the AT contents of the sequenced species were higher than their GC contents. The most frequently used start/stop codon was ATN/TAA. GTG was only found in atp6 and atp8 of E. gibbosus. All transfer RNA genes (tRNAs) exhibit the typical cloverleaf secondary structure, except for the trnS1 and trnV, which lacks the stem of the DHU arm. The length and copy number of repeat units were conserved within Eysarcoris, with the exception of E. gibbosus. Phylogenetic analyses based on mitogenomes using both maximum likelihood (ML) and Bayesian inference (BI) methods strongly supported the relationship among tribes within Pentatominae and confirmed that Graphosoma should be an intermediate lineage of Pentatominae. The relationship between Eysarcoris and Carbula was strongly supported and combined with our previous geometric morphometrics and chromosomal studies, suggest the Eysarcoris should belong to the tribe Eyasrcorini. This work will help to enhance our understanding of mitochondrial genomic evolution and phylogenetic relationships in Pentatominae.

Fidelity of DNA polymerases in the detection of intraindividual variation of mitochondrial DNA

Here we investigated the consequences of PCR amplification errors in the identification of intraindividual mtDNA variation. The bumblebee Bombus morio was chosen as model for the COI gene amplification tests with two DNA polymerases (Taq and Q5) presenting different error rates. The amplifications using Taq resulted in a significant increase of singleton haplotypes per individual in comparison to Q5. The sequence characteristics indicated that Taq resulted haplotypes are mostly due to amplification errors. Studies focusing on intraindividual variability should address special attention to the DNA polymerase fidelity to avoid overestimation of heteroplasmic haplotypes.

Conserved numts mask a highly divergent mitochondrial-COI gene in a species complex of Australian stingless bees Tetragonula (Hymenoptera: Apidae)

Tetragonula carbonaria, Tetragonula davenporti, Tetragonula hockingsi and Tetragonula mellipes comprise a species complex of Australian stingless bee species known as the 'Carbonaria' group. The species are difficult to distinguish morphologically and the major species-defining characters relate to comb architecture and nest entrance ornamentation. The taxonomy of the group is further complicated by likely nuclear mitochondrial pseudogenes (numts) and inter-specific hybrids. Here we demonstrate the existence of COI numts and isolate and characterize the 'true' mt-COI gene in T. carbonaria and T. hockingsi. Numts were isolated from enriched-nuclear DNA extraction followed by PCR amplification and Sanger sequencing, and were recognized by the presence of deletions and/or premature stop codons in the translated sequences. The mt-COI sequences were obtained from NGS sequencing using purified mtDNA. In T. carbonaria, two numts (numt1 and numt2) were identified and a third (numt3) was identified in T. hockingsi. Numt2 and numt3 are similar (1.2% sequence divergence), indicating a recent common origin. The genetic distance between the mt-COI of the two Tetragonula species was higher than might be expected for closely related species, 16.5%, corroborating previous studies in which T. carbonaria and T. hockingsi were regarded as separate species. The three numts are more similar to the COI of other stingless bee species, including Australian Austroplebia australis and South American Melipona bicolor (81.7-83.9%) than to the mt-COI of their own species (70-71.4%). This is because the mt-COI of T. carbonaria and T. hockingsi differ greatly from other Meliponinae. Our findings explain some formerly puzzling aspects of Carbonaria biogeography, and misinterpreted amplifications.

Haplotype identification and detection of mitochondrial DNA heteroplasmy in Varroa destructor mites using ARMS and PCR-RFLP methods

In the present study, amplification refractory mutation system (ARMS) and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) methods were used for identification of recently described Serbia 1 (S1) and Peshter 1 (P1) mitochondrial haplotypes of Varroa destructor. Based on single nucleotide polymorphisms (SNPs) within cytochrome oxidase 1 (cox1) and cytochrome b (cytb) gene sequences, a total of 64 adult V. destructor females were analyzed from locations where the S1 and P1 haplotypes had been detected previously. Results of haplotype identification obtained by ARMS and PCR-RFLP methods were completely consistent with the sequencing data. Furthermore, in some analyzed samples the occurrence of site heteroplasmy at haplotype-defining sites was detected, as it was confirmed by double peaks in the sequence chromatograms. Neither mites with simultaneous nucleotide variability, nor those with combined SNP and heteroplasmy in cox1 and cytb were found. Given that this is the first occurrence of site heteroplasmy in V. destructor, the origin of this phenomenon and possible specific traits of heteroplasmic mites have yet to be determined.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness