AT-rich Control Region Research Articles

Elucidating the Mitogenomic Blueprint of Pomadasys perotaei from the Eastern Atlantic: Characterization and Matrilineal Phylogenetic Insights into Haemulid Grunts (Teleostei: Lutjaniformes).

The parrot grunt fish, Pomadasys perotaei, has a limited distribution in the Eastern Atlantic Ocean and is an important species in marine capture fisheries across several West African countries. Despite its ecological and economic significance, the mitogenomic information for this species is lacking. This study utilized next-generation sequencing to generate the de novo mitogenome of P. perotaei from Eastern Atlantic. The resulting mitogenome is 16,691 base pairs and includes 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, and an AT-rich control region (CR). Most of the PCGs exhibit nonsynonymous (Ka) and synonymous (Ks) substitution rates of less than '1', indicating strong negative selection across haemulid fishes. The control region of Pomadasys species contains four conserved domains, as seen in other teleost's, with polymorphic nucleotides that can be used to study population structures through the amplification of short mitochondrial gene fragments. Additionally, Bayesian phylogenetic analysis based on PCGs revealed a non-monophyletic clustering pattern of Pomadasys within the haemulid matrilineal tree. Overall, the structural characterization and phylogenetic analysis enhance our understanding of the genetic composition and evolutionary history of Pomadasys species from the Indo-West Pacific and Eastern Atlantic Oceans.

Mitogenomic Architecture of Atlantic Emperor Lethrinus atlanticus (Actinopterygii: Spariformes): Insights into the Lineage Diversification in Atlantic Ocean.

The evolutionary history of emperors, particularly in the Atlantic and Indo-West Pacific Oceans, remains largely unmapped. This study explores the maternal lineage evolution of Lethrinids by examining the complete mitogenome of Lethrinus atlanticus, which is endemic to the Eastern Atlantic Ocean. Utilizing advanced next-generation sequencing, we found that the mitogenome spans 16,789 base pairs and encompasses 37 genes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and an AT-rich control region (CR). Our analysis indicates a preference for AT base pairs in the L. atlanticus mitogenome (53.10%). Most PCGs begin with the ATG codon, except for COI, which starts with GTG. Relative synonymous codon usage reveals high frequencies for alanine, leucine, proline, serine, and threonine. The ratio of nonsynonymous to synonymous substitutions suggests strong negative selection across all PCGs in Lethrinus species. Most transfer RNAs exhibit typical cloverleaf structures, with the exception of tRNA-serine (GCT), which lacks a dihydrouracil stem. Comparative analysis of conserved sequence blocks across the CRs of three Lethrinus species shows notable differences in length and nucleotide composition. Phylogenetic analysis using concatenated PCGs clearly distinguishes all Lethrinus species, including L. atlanticus, and sheds light on the evolutionary relationships among Spariformes species. The estimated divergence time of approximately 20.67 million years between L. atlanticus and its Indo-West Pacific relatives provides insights into their historical separation and colonization during the late Oligocene. The distribution of Lethrinids may be influenced by ocean currents and ecological factors, potentially leading to their speciation across the Eastern Atlantic and Indo-West Pacific. This study enhances our understanding of the genetic diversity and phylogenetic relationships within Lethrinus species. Further exploration of other emperor fish mitogenomes and comprehensive genomic data could provide vital insights into their genetic makeup, evolutionary history, and environmental adaptability in marine ecosystems globally.

The Complete Mitochondrial Genome and Phylogenetic Analysis of Rhagastis binoculata (Matsumura, 1909) (Lepidoptera: Sphingidae).

The mitochondrial genome (mitogenome) Rhagastis binoculata (Matsumura, 1909), an endemic moth species in Taiwan, was sequenced and analyzed. The complete circular mitogenome of R. binoculata is 15,303 bp and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region. The mitogenome has an overall nucleotide composition of 41.2% A, 11.9% C, 7.5% G, and 39.4% T, with an AT content of 80.6%. Of the protein-coding genes (PCGs), 12 start with ATG, ATT, and ATC, and COX1 starts with a "CGA" codon. All of the stop codons are "TAA, TAG, or T". Our phylogenetic analysis of 21 species of Sphingidae insects suggests that R. binoculata is clustered with Rhagastis mongoliana, which belongs to the subfamily Macroglossinae.

Phylogenetic and Comparative Genomics Study of Papilionidae Based on Mitochondrial Genomes.

Most species of Papilionidae are large and beautiful ornamental butterflies. They are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. Complete mitochondrial genomes (mitogenomes) have been widely used in phylogenetic studies of butterflies, but mitogenome knowledge within the family Papilionidae is limited, and its phylogeny is far from resolved. In this study, we first report the mitogenome of Byasa confusa from the subfamily Papilioninae of Papilionidae. The mitogenome of B. confusa is 15,135 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region (CR), closely mirroring the genomic structure observed in related butterfly species. Comparative analysis of 77 Papilionidae mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, Ka/Ks, and relative synonymous codon usage (RSCU) are all consistent with that of other reported butterfly mitogenomes. We conducted phylogenetic analyses using maximum-likelihood (ML) and Bayesian-inference (BI) methods, with 77 Papilionidae species as ingroups and two species of Nymphalidae and Lycaenidae as outgroups. The phylogenetic analysis indicated that B. confusa were clustered within Byasa. The phylogenetic trees show the monophyly of the subfamily Papilioninae and the tribes Leptocircini, Papilionini, and Troidini. The data supported the following relationships in tribe level on Papilioninae: (((Troidini + Papilionini) + Teinopalpini) + Leptocircini). The divergence time analysis suggests that Papilionidae originated in the late Creataceous. Overall, utilizing the largest number of Papilionidae mitogenomes sequenced to date, with the current first exploration in a phylogenetic analysis on Papilionidae (including four subfamilies), this study comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of the Papilionidae family.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae).

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

Matrilineal phylogeny and habitat suitability of the endangered spotted pond turtle (Geoclemys hamiltonii; Testudines: Geoemydidae): a two-dimensional approach to forecasting future conservation consequences.

The spotted pond turtle (Geoclemys hamiltonii) is a threatened and less explored species endemic to Bangladesh, India, Nepal, and Pakistan. To infer structural variation and matrilineal phylogenetic interpretation, the present research decoded the mitogenome of G. hamiltonii (16,509 bp) using next-generation sequencing technology. The mitogenome comprises 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one AT-rich control region (CR) with similar strand symmetry in vertebrates. The ATG was identified as a start codon in most of the PCGs except Cytochrome oxidase subunit 1 (cox1), which started with the GTG codon. The non-coding CR of G. hamiltonii was determined to have a unique structure and variation in different domains and stem-loop secondary structure as compared with other Batagurinae species. The PCGs-based Bayesian phylogeny inferred strong monophyletic support for all Batagurinae species and confirmed the sister relationship of G. hamiltonii with Pangshura and Batagur taxa. We recommend generating more mitogenomic data for other Batagurinae species to confirm their population structure and evolutionary relationships. In addition, the present study aims to infer the habitat suitability and habitat quality of G. hamiltonii in its global distribution, both in the present and future climatic scenarios. We identify that only 58,542 km2 (7.16%) of the total range extent (817,341 km2) is suitable for this species, along with the fragmented habitats in both the eastern and western ranges. Comparative habitat quality assessment suggests the level of patch shape in the western range is higher (71.3%) compared to the eastern range. Our results suggest a massive decline of approximately 65.73% to 70.31% and 70.53% to 75.30% under ssp245 and ssp585 future scenarios, respectively, for the years between 2021-2040 and 2061-2080 compared with the current distribution. The present study indicates that proper conservation management requires greater attention to the causes and solutions to the fragmented distribution and safeguarding of this endangered species in the Indus, Ganges, and Brahmaputra (IGB) river basins.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

First Mitogenome of Endangered Enteromius thysi (Actinopterygii: Cypriniformes: Cyprinidae) from Africa: Characterization and Phylogeny

The complete mitochondrial genome of endangered Enteromius thysi was determined from Cameroon in Western Africa. The genome was 16,688 bp in length, comprising 37 genes (13 PCGs, 2 rRNAs, 22 tRNAs, and an AT-rich control region). The heavy strand accommodates 28 genes (12 PCGs, 2 rRNAs, and 14 tRNAs), whereas the light strand holds 9 genes (NAD6 and 8 tRNAs). The E. thysi mitogenome is AT-biased (60.5%), as exhibited in other Enteromius species. Most of the PCGs start with the ATG initiation codon, except COI, with GTG, and seven PCGs end with the TAA termination codon, except some with an incomplete termination codon. Most of the tRNAs showed classical cloverleaf secondary structures, except tRNA-serine (trnS1). Bayesian phylogeny distinctly separated E. thysi from other congeners. The control regions of Enteromius species exhibited highly variable nucleotides, and parsimony-informative sites were found in the conserved sequence block-III (CSB-III) compared with other domains and a unique 9 bp tandem repeat (ATGCATGGT) in the variable-number tandem repeats (VNTRs) region of E. thysi. The present phylogeny with limited mitogenomes showed an uneven diversity and evolutionary patterns of Enteromius species distributed in the northwestern and southeastern riverine systems in Africa, which warrants further investigation. Based on the results of the present study, we recommend additional surveys, in-depth taxonomy, and the generation of more mitogenomes that could resolve the diversification pattern of Enteromius species in Africa.

Characterization of the complete mitogenome of Trachylophus sinensis (Coleoptera: Cerambycidae: Cerambycinae), the type species of Trachylophus and its phylogenetic implications

• The first complete mitogenome from the genus Trachylophus was determined using next-generation sequencing. • All available mitogenome of Cerambycinae were comparative and analysis to better understand of their mitogenomic evolution patterns. • The phylogenetic relationship of Cerambycinae was explored through all available mitogenome downloaded from GenBank. Complete mitochondrial genomes (mitogenomes) have long been proved as reliable markers for phylogenetic reconstruction among diverse animal groups, especially benefited from recent rapid development of sequencing techniques. However, the mitogenomes of many important clades remain poorly represented, which restricted the understanding of macroscale evolutionary history of these groups. Here, we sequenced and characterized the complete mitogenome of Trachylophus sinensis , a type species of the Trachylophus genus, which also represents the first sequenced mitogenome in this genus. The complete circular mitogenome was 15,746 bp in length, containing 37 typical genes and one noncoding AT-rich control region. The nucleotide composition of the mitogenome was highly A + T biased, accounting for 70.07% of the whole mitogenome with a slightly positive AT skewness (0.106). The 13 Protein coding genes (PCGs) used ATN as their start codons, except nad1 which used TTG. All tRNA genes were predicted with a characteristic cloverleaf secondary structure except trnS1 (AGN) , whose dihydrouridine (DHU) arm was replaced by a simple loop. Phylogenetic analyses recovered Cerambycinae as a monophyletic group with high node supports and the sister relationship between T. sinensis and Nadezhdiella cantori . However, we found that deeper nodes showed not strong support, which may be caused by limited taxa sampling in our study. More mitogenomes should be sequenced representing various taxonomic levels, especially closely related species, which will enhance our understanding of phylogenetic relationships among Cerambycinae.

The first complete mitochondrial genome of Matsucoccidae (Hemiptera, Coccoidea) and implications for its phylogenetic position.

The mitochondrial genome (mitogenome) has been extensively used to better understand the phylogenetic relationships within the hemipteran suborder Sternorrhyncha, but sequenced mitogenomes remain unavailable for the entire family Matsucoccidae to date. To address this, here we sequenced the complete mitogenome of Matsucoccusmatsumurae; the first for this family. The mitogenome is 15,360 bp in size and comprises the typical set of 37 mitochondrial genes and a large non-coding region (AT-rich region). Gene order, nucleotide composition and codon usage of protein-coding genes (PCGs) of M.matsumurae differ considerably from those of the other two sequenced Coccidae species. All PCGs were initiated by the ATN start codons and ended with the TAA/G or single T-- stop codons. Nine transfer RNA genes could be folded into typical clover-leaf secondary structures. The length and AT content of the ribosomal RNA genes are highly conserved in the Coccoidea mitogenomes. In contrast, the AT-rich control region is highly variable in size and in the number of tandem repeats. The sliding window analysis showed that the cox1 gene is the most conserved amongst the 13 PCGs, while the ratios of non-synonymous to synonymous substitution rates indicated that the evolution of this mitogenome has been dominated by positive selection. Phylogenetic analyses, based on nucleotide sequence data of 37 mitochondrial genes and amino acid sequence data of 13 PCGs using Bayesian Inference and Maximum Likelihood methods, showed that Matsucoccidae diverged before the Coccidae.

Complete mitochondrial genomes of Papilio nephelus chaon and Papilio epycides (Lepidoptera: Papilionidae: Papilioninae) and phylogenetic analysis

The complete mitochondrial genome (mitogenome) sequences of Papilio nephelus chaon and Papilio epycides were sequenced by Illumina and analyzed in this study. They are 15,287 bp and 15,012 bp in size, respectively, and contains 13 protein-coding genes (PCGs), 22 tRNA genes (tRNAs), 2 rRNA genes (rRNAs), and 1 AT-rich control region (CR). The phylogenetic relationships of 56 species in the Papilionidae were inferred based on concatenated nucleotide sequences by using Maximum Likelihood with the selected best-fit model GTR + F+R6. The phylogenetic analysis showed that P. nephelus chaon and P. epycides were located in the genus Papilio. This study provides a basis for further study on mitogenome and phylogenetics of the Papilionidae.

The Carniolan Honeybee from Slovenia-A Complete and Annotated Mitochondrial Genome with Comparisons to Closely Related Apis mellifera Subspecies.

Simple SummaryThe western honeybee, Apis mellifera, is a globally distributed bee species with many recognised subspecies, one of which is Apis mellifera carnica, the Carniolan honeybee. Apis m. carnica is native to southern Central Europe and parts of the Balkans, with the locus classicus in Slovenia. It is also widely popular with beekeepers in parts of Central and Northern Europe and other parts of the world, including the USA, Canada, and even New Zealand. In Slovenia, A. m. carnica is protected, with measures to conserve the subspecies’ autochthonous domestic population in place. Such efforts depend heavily upon genomic and phylogenetic information. In this study, we sequenced and annotated the mitochondrial genome of a specimen from Slovenia and compared the obtained data with a previously published sample of the A. m. carnica from Austria and the closely related Italian honeybee A. m. ligustica. We found several features unique to the new mitochondrial genome. We also phylogenetically analyzed the relationship between our sequence and the selected available A. mellifera mitochondrial sequences. The acquired position of the sequenced A. m. carnica from Slovenia on the phylogenetic tree brings new evidence for close relationships among C and O lineages and reflects their recent historical matrilinear ancestry.The complete mitochondrial genome of the Carniolan honeybee (Apis mellifera carnica) from Slovenia, a homeland of this subspecies, was acquired in two contigs from WGS data and annotated. The newly obtained mitochondrial genome is a circular closed loop of 16,447 bp. It comprises 37 genes (13 protein coding genes, 22 tRNA genes, and 2 rRNA genes) and an AT-rich control region. The order of the tRNA genes resembles the order characteristic of A. mellifera. The mitogenomic sequence of A. m. carnica from Slovenia contains 44 uniquely coded sites in comparison to the closely related subspecies A. m. ligustica and to A. m. carnica from Austria. Furthermore, 24 differences were recognised in comparison between A. m. carnica and A. m. ligustica subspecies. Among them, there are three SNPs that affect translation in the nd2, nd4, and cox2 genes, respectively. The phylogenetic placement of A. m. carnica from Slovenia within C lineage deviates from the expected position and changes the perspective on relationship between C and O lineages. The results of this study represent a valuable addition to the information available in the phylogenomic studies of A. mellifera—a pollinator species of worldwide importance. Such genomic information is essential for this local subspecies’ conservation and preservation as well as its breeding and selection.

The complete mitochondrial genome of major malaria vector Anopheles anthropophagus (Diptera: Culicidae) in China

Anopheles anthropophagus (Xu and Feng 1975) is the major vector of malaria in Eastern and Southern China. The species An. anthropophagus is considered a synonym of An. lesteri (Baisas & Hu, 1936), although they differ in several key biological characteristics. Here, we report the complete mitochondrial genome of An. anthropophagus for the first time. The mitogenome of An. anthropophagus is a typical circular, double-stranded molecule with a total length of 15,413 base pairs, and contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and an AT-rich control region. A phylogenetic analysis of the complete mitogenomes of 16 species of Anopheles (Culicidae) revealed that An. anthropophagus is closely related to An. sinensis (Wiedemann 1828), in the family Culicidae. The An. anthropophagus mitogenome provides new data for further taxonomic and phylogenetic studies of the genus Anopheles.

Next-generation sequencing of the mitochondrial genome of Bombus longipennis Friese, 1918 (Hymenoptera: Apidae)

Bombus longipennis is the species of Bumblebees (Hymenoptera: Apidae), which are important pollinators for wild plants and greenhouses crops. The complete mitochondrial genome (mitogenome) of B. longipennis was determined by next-generation sequencing. The mitogenome was 18,458 bp in size with 87.2% A + T content, containing 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an AT-rich control region (D-loop). Gene arrangement was found to be identical to those of other mitogenomes of bumblebees (e.g. Bombus terrestris and Bombus ignitus). All 13 PCGs initiated with typical ATN codons. Among them, 11 PCGs terminated with TAA or TAG; only cox2 and nad4 have incomplete stop codon T. All 22 tRNAs can be folded into typical cloverleaf structure. Phylogenetic analysis based on the concatenated nucleotide sequences of all 13 PCGs indicated that B. longipennis was more closely related to other species of subgenus Bombus, which clustered into a monophyletic group.

The complete mitochondrial genome of Dudusa sphingiformis (Lepidoptera: Notodontidae) and phylogenetic analysis

Dudusa sphingiformis is an important lepidopteran pest widely distributed in tropical and subtropical zones of Asia. In this paper, the complete mitochondrial genome (mitogenome) of D. sphingiformis was determined by next-generation sequencing. The mitogenome was 15,806 bp in length, comprising 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an AT-rich control region (D-loop). The gene arrangement of this mitogenome was identical to that of the previous studies of Notodontidae moths. Almost all the PCGs initiated with typical ATN codons, except for cox1 with CGA. Among them, nine PCGs terminated with TAA or TAG, while other four PCGs (cox1, cox2, nad5, and nad4) with incomplete stop codon T. All the 22 tRNAs had the typical cloverleaf structure, except for trnS1, whose dihydrouridine (DHU) arm forms a simple loop. Phylogenetic analysis based on the concatenated nucleotide sequences of 13 PCGs indicated that D. sphingiformis was more closely related to other species of family Notodontidae, forming a monophyletic group, with well-resolved relationships among five family of Noctuoidea.

The complete mitochondrial genome of Bombus filchnerae (Hymenoptera: Apidae) and phylogenetic analysis

Here, we report the complete mitochondrial genome (mitogenome) of Bombus filchnerae (Hymenoptera: Apidae). The genome size of B. filchnerae was 18,553 bp with 88.7% A + T content, containing 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an AT-rich control region (D-loop). tRNA rearrangement was observed in this mitochondrial genome when compared to those of other bumblebee species (e.g. Bombus breviceps and Bombus asiaticus). All the 13 PCGs initiated with typical ATN codons. Among them, 11 PCGs terminated with TAA, only nad4 and nad5 with incomplete stop codon TA and T, respectively. All the 22 tRNAs can be folded into typical cloverleaf structure, except for trnS1, whose dihydrouridine (DHU) arm forms a simple loop. The phylogenetic analysis based on the concatenated nucleotide sequences of all 13 PCGs indicated that B. filchnerae showed the closest relationship with Bombus pascuorum, forming a mono clade of the subgenus Thoracobombus, with well-resolved relationships among nine Bombus subgenera.

Characterizing the Mitogenome of the Endemic Bumblebee Subspecies from the Canary Islands for Conservation Purposes

The taxonomic status of Bombus terrestris subspecies is complex and has deep implications in the management of commercial bumblebees for crop pollination as well as in the establishment of appropriate conservation plans. Herein, the complete mitogenome of the endemic Canary Islands subspecies Bombus terrestris canariensis is newly sequenced and compared with available mitochondrial sequences in order to shed light into its taxonomic status. The mitochondrial genome was 17,300 bp in length and contained 37 genes, including 13 protein-coding genes (PCGs), two rRNAs, and 22 tRNAs and a partial sequence of the AT rich control region. The phylogenetic analysis of PCGs of the mitogenome was congruent with its subspecific status and a close relationship with the North African subspecies africanus as previously suggested. The sequencing of the mitogenome of B. t. canariensis provides useful genetic information to study the conservation genetics and genetic diversity of these island bumblebee populations.

Characterisation of the complete mitochondrial genome of the insect-parasitic nematode Heterorhabditis bacteriophora: an idiosyncratic gene order and the presence of multiple long non-coding regions

Summary We present here the complete mtDNA genome (mitogenome) of Heterorhabditis bacteriophora, an important biological control agent of soil-dwelling insect pests in agriculture and horticulture. This is the first description of a mitogenome for a member of the family Heterorhabditidae. The genome contains the typical chromadorean complement of 12 protein-coding genes, 22 tRNA genes and two rRNA genes. All genes are transcribed in the same direction and have a nucleotide composition high in A and T. For the entire genome, the nucleotide contents are 47.02% (T), 28.81% (A), 16.10% (G), 8.08% (C) and 75.83% (AT). Heterorhabditis bacteriophora has a unique, idiosyncratic gene arrangement. It differs from that of Caenorhabditis elegans in having a block of seven genes: trnQ-trnF-cytb-trnL1-cox3-trnT-nad4 translocated to a position between nad3 and nad5, as well as having a change in the position of the four tRNA block gene cluster, trnC-trnM-trnD-trnG, where trnC and trnM have switched places and trnD and trnG have translocated between nad4 and nad5 genes. The H. bacteriophora mitogenome is 18 128 bp long, and thus is ca 4 kb larger than the mitogenomes of most chromadoreans. This relatively large genome is due to the presence of five non-coding regions (NCR): NCR1 (114 bp), NCR2 (159 bp), NCR3 (498 bp), NCR4 (1917 bp) and NCR5 (2154 bp), which make up 26.7% of the genome. The NCR5 had the highest A + T content of 83.47% indicating that this region is the likely AT-rich control region. The complete 498 bp NCR3 sequence is duplicated in NCR4 and in NCR5 (the putative AT-rich control region). Such an organisation has not been reported previously in nematode mtDNA.

The complete mitochondrial genome of the bagworm from a tea plantation in China, Eumeta variegata (Lepidoptera: Psychidae)

The mitochondrial (mt) genome of Eumeta variegata Snellen (Psychidae) has been sequenced and annotated. The mt genome has a total length of 15,793 bp, consisting of 13 protein-coding genes, 22 tRNA, two rRNA genes, and an AT-rich control region (GenBank accession no. MN242985). The nucleotide composition was extremely AT-rich, and the AT content is 81.57%. The gene order is consistent with other sequenced mt genome of moths and butterflies from Ditrysia. The sequence similarity of E. variegate mt genomes between the specimen of China and South Korea is 98.38%, whereas the similarity between the specimen of China and Japan is 90.61%. The sequence of PCGs and rRNAs among different specimens are similar, and many differences are detected at the region of A + T-rich region and the tRNA block ‘ARNS1EF’.

First mtgenome sequences from three genera and phylogenetic relationships of the family Apidae based on mtgenome sequences (Hymenoptera: Apoidea)

In this study, we report the complete mitochondrial genomes (mtgenome) of Thyreus decorus, Ceratina okinawana and Amegilla calceifera, which are the first time of mtgenome report also for the genera Thyreus, Ceratina and Amegilla in the family Apidae. They contain 15,237, 15,207, and 17,728 bp, with AT content of 84.97%, 79.30%, and 84.63%, respectively. Each mtgenome includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA (16S and 12S rRNA) and an AT-rich control region (CR). The phylogenetic relationships of 45 species in the family were constructed using Bayesian Inference based on concatenated nucleotide sequences of 13 PCGs. Our study suggests that the subfamily Apinae is a paraphyletic group, with the genus Eucera claded into the subfamily Xylocopinae and the genera Amegilla and Thyreus into the subfamily Nomadinae. In Apinae, the genera Melipona and Bombus are significantly sister group, and the genus Apis is the sister group with Melipona + Bombus.