The first complete mitochondrial genome of Dacus longicornis (Diptera: Tephritidae) using next-generation sequencing and mitochondrial genome phylogeny of Dacini tribe

Revisiting the History and Biogeography of Bactrocera oleae and Other Olive-Feeding Fruit Flies in Africa and Asia

Bactrocera oleae Rossi causes important agricultural losses in olive growing regions. Despite its economic relevance, the expansion history and biogeography of B. oleae and other olive-feeding fruit flies remain unclear. We used mitogenomic data of Bactrocera species from a broad geographic range to explore possible historical biogeographic patterns in B. oleae and other olive-feeding fruit flies. Our data suggest that (1) the transition from stenophagy on Oleaceae to oligophagy on Olea most likely occurred in Africa more than 6 million years ago (MYA), and (2) the subsequent transition to monophagy on Olea europaea took place in Asia or in Africa in the Early Pleistocene. Our results support equally the hypotheses that the ancestors of modern B. oleae underwent two waves of migration from Asia to Africa or that they zigzagged between Asia and Africa.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genomes of Two Tiny Necrophagous Phorid Flies, Metopina sagittata and Puliciphora borinquenensis (Diptera: Phoridae)

The mitochondrial genome is frequently used for species identification and phylogenetic studies. In this study, we first sequenced and annotated the complete mitochondrial genomes of two phorid species that are forensically important in buried or enclosed environments: Metopina sagittata (Liu) and Puliciphora borinquenensis (Wheeler). The complete mitochondrial genome sequences of M. sagittata and P. borinquenensis were 15,640 bp with an A+T content of 75.97% and 15,429 bp with an A+T content of 75.38%, respectively. Their circular genomes both contained 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region located between rrnS and trnI which was 808 bp for M. sagittata and 746 bp for P. borinquenensis. All the PCGs of both species started with ATN codons except for cox1 which used TTG codon. In addition to the common stop codon TAA and TAG, the incomplete stop codon T was used in two PCGs (cox1 and nad4) of M. sagittata and five PCGs (cox1, cox2, cox3, nad5, and nad4) of P. borinquenensis. There were 3 and 10 mismatched base pairs in the tRNA secondary structures from M. sagittata and P. borinquenensis, respectively. Both maximum likelihood and Bayesian inference analyses indicated that Platypezidae and Phoridae are sister taxa. M. sagittata is closely related to P. borinquenensis within the subfamily Metopininae. This work enhances the databases of Phoridae genomes and contributes to the further study of species identification and phylogenetics of this family.

Complete mitochondrial genome of Dacus vijaysegarani and phylogenetic relationships with congeners and other tephritid fruit flies (Insecta: Diptera)

Background

Tephritid fruit flies of the genus Dacus are members of the tribe Dacini, subfamily Dacinae. There are some 274 species worldwide, distributed in Africa and the Asia-Pacific. To date, only five complete mitochondrial genomes (mitogenomes) of Dacus fruit flies have been published and are available in the GenBank.

Methods and results

In view of the lack of study on their mitogenome, we sequenced (by next generation sequencing) and annotated the complete mitogenome of D. vijaysegarani from Malaysia to determine its features and phylogenetic relationship. The whole mitogenome of D. vijaysegarani has identical gene order with the published mitogenomes of the genus Dacus, with 13 protein-coding genes, two rRNA genes, 22 tRNAs, a non-coding A + T rich control region, and intergenic spacer and overlap sequences. Phylogenetic analysis based on 15 mitochondrial genes (13 PCGs and two rRNA genes), reveals Dacus, Zeugodacus and Bactrocera forming a distinct clade. The genus Dacus forms a monophyletic group in the subclade containing also the Zeugodacus group; this Dacus-Zeugodacus subclade is distinct from the Bactrocera subclade. D. (Mellesis) vijaysegarani forms a lineage with D. (Mellesis) trimacula in the subcluster containing also the lineage of D. (Mellesis) conopsoides and D. (Callantra) longicornis. D. (Dacus) bivittatus and D. (Didacus) ciliatus form a distinct subcluster. Based on cox1 sequences, the Malaysia and Vietnam taxa of D. vijaysegarani may not be conspecific.

Conclusions

Overall, the mitochondrial genome of D. vijaysegarani provided essential molecular data that could be useful for further studies for species diagnosis, evolution and phylogeny research of other tephritid fruit flies in the future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06608-2.

Mitochondrial genome sequencing and phylogeny of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis (Diptera: Culicidae)

The genus Haemagogus (Diptera: Culicidae) comprises species of great epidemiological relevance, involved in transmission cycles of the Yellow fever virus and other arboviruses in South America. So far, only Haemagogus janthinomys has complete mitochondrial sequences available. Given the unavailability of information related to aspects of the evolutionary biology and molecular taxonomy of this genus, we report here, the first sequencing of the mitogenomes of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis. The mitogenomes showed an average length of 15,038 bp, average AT content of 79.3%, positive AT-skews, negative GC-skews, and comprised 37 functional subunits (13 PCGs, 22 tRNA, and 02 rRNA). The PCGs showed ATN as start codon, TAA as stop codon, and signs of purifying selection. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogenetic analyzes of Bayesian inference and Maximum Likelihood, based on concatenated sequences from all 13 PCGs, produced identical topologies and strongly supported the monophyletic relationship between the Haemagogus and Conopostegus subgenera, and corroborated with the known taxonomic classification of the evaluated taxa, based on external morphological aspects. The information produced on the mitogenomes of the Haemagogus species evaluated here may be useful in carrying out future taxonomic and evolutionary studies of the genus.

Development of internal COI primers to improve and extend barcoding of fruit flies (Diptera: Tephritidae: Dacini)

Accurate species-level identifications underpin many aspects of basic and applied biology; however, identifications can be hampered by a lack of discriminating morphological characters, taxonomic expertise or time. Molecular approaches, such as DNA "barcoding" of the cytochrome c oxidase (COI) gene, are argued to overcome these issues. However, nuclear encoding of mitochondrial genes (numts) and poor amplification success of suboptimally preserved specimens can lead to erroneous identifications. One insect group for which these molecular and morphological problems are significant are the dacine fruit flies (Diptera: Tephritidae: Dacini). We addressed these issues associated with COI barcoding in the dacines by first assessing several "universal" COI primers against public mitochondrial genome and numt sequences for dacine taxa. We then modified a set of four primers that more closely matched true dacine COI sequence and amplified two overlapping portions of the COI barcode region. Our new primers were tested alongside universal primers on a selection of dacine species, including both fresh preserved and decades-old dry specimens. Additionally, Bactrocera tryoni mitochondrial and nuclear genomes were compared to identify putative numts. Four numt clades were identified, three of which were amplified using existing universal primers. In contrast, our new primers preferentially amplified the "true" mitochondrial COI barcode in all dacine species tested. The new primers also successfully amplified partial barcodes from dry specimens for which full length barcodes were unobtainable. Thus we recommend these new primers be incorporated into the suites of primers used by diagnosticians and quarantine labs for the accurate identification of dacine species.

The First Complete Mitochondrial Genomes of Two Sibling Species from Nitidulid Beetles Pests

Carpophilus dimidiatus (Fabricius, 1792) and Carpophilus pilosellus Motschulsky are two sibling species and economically important storage pests worldwide. The first complete mitochondrial (mt) genomes of both were sequenced using next-generation sequencing. The mt genomes of C. dimidiatus and C. pilosellus are circular, with total lengths of 15,717 bp and 15,686 bp, respectively. Gene order and content for both species are similar to what has been observed in ancestral insects and consist of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region. Comparing the mt genome data of C. dimidiatus and C. pilosellus, they are similar in organization, arrangement patterns, GC contents, transfer RNA (tRNA) secondary structures, and codon usage patterns. Small differences were noted with regards to the nucleotide similarity of coding regions and the control region. This is the first publication of the complete mitochondrial genomes of two sibling species. The mt genome sequences can supplement the nuclear markers of the Carpophilus genus in research species identification, system evolution, and population genetic structure, and also will be valuable molecular marker for further genetic studies.

Τhe Complete Mitochondrial Genome of Bactrocera carambolae (Diptera: Tephritidae): Genome Description and Phylogenetic Implications

Bactrocera carambolae is one of the approximately 100 sibling species of the Bactrocera dorsalis complex and considered to be very closely related to B. dorsalis. Due to their high morphological similarity and overlapping distribution, as well as to their economic impact and quarantine status, the development of reliable markers for species delimitation between the two taxa is of great importance. Here we present the complete mitochondrial genome of B. carambolae sourced from its native range in Malaysia and its invaded territory in Suriname. The mitogenome of B. carambolae presents the typical organization of an insect mitochondrion. Comparisons of the analyzed B. carambolae sequences to all available complete mitochondrial sequences of B. dorsalis revealed several species-specific polymorphic sites. Phylogenetic analysis based on Bactrocera mitogenomes supports that B. carambolae is a differentiated taxon though closely related to B. dorsalis. The present complete mitochondrial sequences of B. carambolae could be used, in the frame of Integrative Taxonomy, for species discrimination and resolution of the phylogenetic relationships within this taxonomically challenging complex, which would facilitate the application of species-specific population suppression strategies, such as the sterile insect technique.

Complete mitochondrial genome of Dacus trimacula (Diptera: Tephritidae) using next-generation sequencing from China

The complete mitochondrial genome (mitogenome) of the fruit fly species Dacus trimacula (Diptera: Tephritidae: Dacinae) are sequenced and annotated. The mitochondrial genome is 15,851 bp (GenBank No. MK940811) has an A + T content of 72.8% (A 39.2%; C 17.0%; G 10.2%, and T 33.6%), which is the classical structure for insect mitogenome. All PCGs started with ATN and stopped with TAN. The phylogenetic tree confirms that D. trimacula clustered with D. longicornis as the sister group to Zeugodacus. This study enriches the mitogenomes of the fruit flies.

The first complete mitochondrial genome of Bactrocera tsuneonis (Miyake) (Diptera: Tephritidae) by next-generation sequencing and its phylogenetic implications

Bactrocera tsuneonis (Miyake), generally known as the Japanese orange fly, is considered to be a major pest of commercial citrus crops. It has a limited distribution in China, Japan and Vietnam, but it has the potential to invade areas outside of Asia. More genetic information of B. tsuneonis should be obtained in order to develop effective methodologies for rapid and accurate molecular identification due to the difficulty of distinguishing it from Bactrocera minax based on morphological features. We report here the whole mitochondrial genome of B. tsuneonis sequenced by next-generation sequencing. This mitogenome sequence had a total length of 15,865 bp, a typical circular molecule comprising 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and a non-coding region (A + T-rich control region). The structure and organization of the molecule were typical and similar compared with the published homologous sequences of other fruit flies in Tephritidae. The phylogenetic analyses based on the mitochondrial genome data presented a close genetic relationship between B. tsuneonis and B. minax. This is the first report of the complete mitochondrial genome of B. tsuneonis, and it can be used in further studies of species diagnosis, evolutionary biology, prevention and control.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

A conserved motif within cox 2 allows broad detection of economically important fruit flies (Diptera: Tephritidae)

The genera Anastrepha, Bactrocera, Ceratitis, Dacus and Rhagoletis in the family Tephritidae order Diptera are economically important, worldwide distributed and cause damage to a large number of commercially produced fruits and vegetables. China had regulated these five genera as quarantine pests, including the species Carpomya vesuviana. An accurate molecular method not depending on morphology able to detect all the quarantine fruit flies simultaneously is required for quarantine monitoring. This study contributes a comparative analysis of 146 mitochondrial genomes of Diptera species and found variable sites at the mt DNA cox2 gene only conserved in economically important fruit flies species. Degenerate primers (TephFdeg/TephR) were designed specific for the economically important fruit flies. A 603 bp fragment was amplified after testing each of the 40 selected representative species belonging to each economically important Tephritid genera, no diagnostic fragments were detected/amplified in any of the other Tephritidae and Diptera species examined. PCR sensitivity assays demonstrated the limit of detection of targeted DNA was 0.1 ng/μl. This work contributes an innovative approach for detecting all reported economically important fruit flies in a single-step PCR specific for reported fruit fly species of quarantine concern in China.

Complete mitochondrial genome of Zeugodacus tau (Insecta: Tephritidae) and differentiation of Z. tau species complex by mitochondrial cytochrome c oxidase subunit I gene

The tephritid fruit fly Zeugodacus tau (Walker) is a polyphagous fruit pest of economic importance in Asia. Studies based on genetic markers indicate that it forms a species complex. We report here (1) the complete mitogenome of Z. tau from Malaysia and comparison with that of China as well as the mitogenome of other congeners, and (2) the relationship of Z. tau taxa from different geographical regions based on sequences of cytochrome c oxidase subunit I gene. The complete mitogenome of Z. tau had a total length of 15631 bp for the Malaysian specimen (ZT3) and 15835 bp for the China specimen (ZT1), with similar gene order comprising 37 genes (13 protein-coding genes-PCGs, 2 rRNA genes, and 22 tRNA genes) and a non-coding A + T-rich control region (D-loop). Based on 13 PCGs and 15 mt-genes, Z. tau NC_027290 (China) and Z. tau ZT1 (China) formed a sister group in the lineage containing also Z. tau ZT3 (Malaysia). Phylogenetic analysis based on partial sequences of cox1 gene indicates that the taxa from China, Japan, Laos, Malaysia, Bangladesh, India, Sri Lanka, and Z. tau sp. A from Thailand belong to Z. tau sensu stricto. A complete cox1 gene (or 13 PCGs or 15 mt-genes) instead of partial sequence is more appropriate for determining phylogenetic relationship.

The chromosomes and the mitogenome of Ceratitis fasciventris (Diptera: Tephritidae): two genetic approaches towards the Ceratitis FAR species complex resolution

Ceratitis fasciventris is a serious agricultural pest of the Tephritidae family that belongs to the African Ceratitis FAR species complex. Species limits within the FAR complex are obscure and multidisciplinary approaches have attempted to resolve phylogenetic relationships among its members. These studies support the existence of at least three additional species in the complex, C. anonnae, C. rosa and C. quilicii, while they indicate the presence of two structured populations (F1 and F2) within the C. fasciventris species. In the present study we present the mitotic karyotype, polytene chromosome maps, in situ hybridization data and the complete mitochondrial genome sequence of an F2 population of C. fasciventris. This is the first polytene chromosome map and complete mitogenome of a member of the FAR complex and only the second reported for the Ceratitis genus. Both polytene chromosomes and mitochondrial sequence could provide valuable information and be used as reference for comparative analysis among the members of the complex towards the clarification of their phylogenetic relationships.