==== Front Mitochondrial DNA B Resour Mitochondrial DNA B Resour Mitochondrial DNA. Part B, Resources 2380-2359 Taylor & Francis 33366401 10.1080/23802359.2019.1693920 1693920 Version of Record Research Article Mitogenome Announcement The mitochondrial genome of the Spanish honey bee, Apis mellifera iberiensis (Insecta: Hymenoptera: Apidae), from Portugal L. Boardman et al. Mitochondrial DNA Part B https://orcid.org/0000-0002-7311-9913 Boardman Leigh a https://orcid.org/0000-0001-6448-5025 Eimanifar Amin a* https://orcid.org/0000-0001-5449-5481 Kimball Rebecca b https://orcid.org/0000-0003-1643-5212 Braun Edward b Fuchs Stefan c Grünewald Bernd c Ellis James D. a a Honey Bee Research and Extension Laboratory, Entomology and Nematology Department, University of Florida, Gainesville, FL, USA; b Department of Biology, University of Florida, Gainesville, FL, USA; c Institut für Bienenkunde, Polytechnische Gesellschaft, Goethe-Universität Frankfurt am Main, Oberursel, Germany * Currently an Independent Senior Research Scientist, Industrial District, Easton, MD, USA. CONTACT Leigh Boardman lboardman@ufl.eduHoney Bee Research and Extension Laboratory, Entomology and Nematology Department, University of Florida, Gainesville, FL, USA 9 12 2019 2020 5 1 1718 © 2019 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group. 2019 The Author(s) https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract The Spanish honey bee Apis mellifera iberiensis, had a mitochondrial genome of 16,560 bp. It consisted of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and an AT-rich control region. The sample was from Portugal and its mitogenome resembled those of the African (A)-lineage honey bee subspecies. It was most closely related to other North African honey bees, namely A. m. sahariensis and A. m. intermissa. Keywords Apis mellifera iberica A-lineage NGS mitogenome United States Department of Agriculture, Animal and Plant Health Inspection Service10.13039/100000199USDA-APHISThis work was supported by the United States Department of Agriculture, Animal and Plant Health Inspection Service (USDA-APHIS) under the cooperative agreements AP17PPQS&T00C172 and AP18PPQS&T00C116. ==== Body The Spanish honey bee, Apis mellifera iberiensis (Engel 1999), is found on the Iberian Peninsula (Spain and Portugal), as well as some nearby islands. Honey bees of this subspecies can have either African (A)-lineage or west Mediterranean and northwest European (M)-lineage mitochondrial haplotypes, with honey bees in the southwest of its range being A-lineage and those in the northeast being M-lineage (Smith et al. 1991; Cánovas et al. 2008; Chávez-Galarza et al. 2015). Here we sequenced an A. m. iberiensis worker bee sampled from Portugal (presumed to represent A-lineage A. m. iberiensis). The sample was obtained from the Ruttner Bee Collection at the Bee Research Institute in Oberursel (Voucher 1964, Dr. W.S. Sheppard, 1991, Evora, Portugal, 38°33 N 7°54 W). The GenBank accession number is MN585110. Genomic DNA was extracted and quantified following Eimanifar et al. (2017) before genomic library preparation and next-generation sequencing (PE-150bp, Illumina Hi-Seq 3000/4000, San Diego, CA) were completed. The quality of sequencing reads was assessed using FastQC (Andrews 2010) and data were trimmed with Trimmomatic (Bolger et al. 2014). Stringent mapping was performed in Geneious Prime 2019.0.4 (Kearse et al. 2012) following Boardman et al. (2019) with the R1 paired trimmomatic output mapped to the reference genome with the highest pairwise identity – in this case Apis mellifera sahariensis (MF351881). The assembled mitogenome was annotated in mitos2 (Bernt et al. 2013) and manually adjusted in Geneious. Phylogenetic relationships were determined using a Mesquite v3.5 (Maddison and Maddison 2018) alignment of the 13 protein-coding genes (PCGs) and two ribosomal RNA (rRNA) regions and RAxML 8.2.10 GTRGAMMA model (1000 bootstrap replicates, -f a option, Stamatakis 2014) on CIPRES Science Gateway v. 3.3 (Miller et al. 2010). The P-distances were generated with PAUP 4.0a (Swofford 2003). The mitogenome of A. m. iberiensis was 16,560 bp long and consisted of 43.2% A, 41.6% T, 9.6% C, and 5.6% G. The location of the 13 PCGs, two rRNAs, and 22 transfer RNAs (tRNAs) resembled arthropod mitogenomes, with nine PCGs (nad2, co1, co2, atp8, atp6, co3, nad3, nad6, and cytb) on the light strand, and four PCGs (nad1, nad4, nad4l, and nad5) and the two rRNAs on the heavy strand. The PCGs used four different start codons: ATT, ATG, ATA, or ATC, and all 13 ended with the TAA stop codon. Atp8 and atp6 overlapped, with 19 shared nucleotides. The 16S rRNA was 1,324 bp with 84% AT, while 12S rRNA was 785 bp with 81.1% AT. The 22 tRNAs were 63 to 78 bp long with the shortest tRNA being glutamate and serine, and the longest, threonine. Phylogenetically, this sample was confirmed as being part of the A-lineage (Figure 1). It was most closely related to North African subspecies A. m. sahariensis (P-distance: 0.00183) and Apis mellifera intermissa (P-distance: 0.00198). Southern African subspecies, Apis mellifera scutellata, Apis mellifera adansonii, Apis mellifera capensis, and Apis mellifera monticola, have P-distances greater than 0.004. The A. mellifera subspecies with the highest P-distance from A. m. iberiensis was Apis mellifera mellifera (P-distance: 0.01395). Figure 1. Phylogenetic tree showing the relationship between Apis mellifera iberiensis (GenBank: MN585110) and 23 other honey bee mitochondrial genomes. The sequences from 13 protein-coding genes and two rRNAs were used to build the tree. The tree is midpoint rooted. Node labels indicate bootstrap values and unlabeled lineages are 100%. The GenBank accession numbers are given after the species names. Subspecies of African (A)-lineage are denoted with (A). Disclosure statement The authors report no conflicts of interest. ==== Refs References Andrews S. 2010. FastQC: a quality control tool for high throughput sequence data. [Accessed 2018 Aug]. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ Bernt M, Pütz J, Florentz C, Donath A, Jühling F, Stadler PF, Externbrink F, Fritzsch G, Middendorf M. 2013. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 69 (2 ):313–319.22982435 Boardman L, Eimanifar A, Kimball RT, Braun EL, Fuchs S, Grunewald B, Ellis JD. 2019. The mitochondrial genome of Apis mellifera simensis (Hymenoptera: Apidae), an Ethiopian honey bee. Mitochondrial DNA Part B. doi:10.1080/23802359.2019.1693307 Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 30 (15 ):2114–2120.24695404 Cánovas F, De La Rúa P, Serrano J, Galián J. 2008. Geographical patterns of mitochondrial DNA variation in Apis mellifera iberiensis (Hymenoptera: Apidae). J Zool Syst Evol Res. 46 (1 ):24–30. Chávez-Galarza J, Henriques D, Johnston JS, Carneiro M, Rufino J, Patton JC, Pinto MA. 2015. Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome-wide nuclear variations provide support for secondary contact from historical refugia. Mol Ecol. 24 (12 ):2973–2992.25930679 Eimanifar A, T. Kimball R, L. Braun E, Fuchs S, Grünewald B, Ellis JD. 2017. The complete mitochondrial genome of Apis mellifera meda (Insecta: Hymenoptera: Apidae). Mitochondrial DNA Part B. 2 (1 ):268–269.33473795 Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. 2012. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 28 (12 ):1647–1649.22543367 Maddison WP, Maddison DR. 2018. Mesquite: a modular system for evolutionary analysis. Version 3.5. [Accessed 2018 Aug]. https://www.mesquiteproject.org/ Miller MA, Pfeiffer W, Schwartz T. 2010. Creating the CIPRES Science Gateway for Inference of Large Phylogenetic Trees. In: Proceedings of the Gateway Computing Environments Workshop (GCE). New Orleans, LA. New York (NY): IEEE; p. 1–8. Smith DR, Palopoli MF, Taylor BR, Garnery L, Cornuet J-M, Solignac M, Brown WM. 1991. Geographical overlap of two mitochondrial genomes in Spanish honeybees (Apis mellifera iberica). J Hered. 82 (2 ):96–100.2013694 Stamatakis A. 2014. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 30 (9 ):1312–1313.24451623 Swofford DL. 2003. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. Sunderland (MA): Sinauer Associates.