Next Article in Journal
Designing Aedes (Diptera: Culicidae) Mosquito Traps: The Evolution of the Male Aedes Sound Trap by Iterative Evaluation
Previous Article in Journal
Nootkatone Is an Effective Repellent against Aedes aegypti and Aedes albopictus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variants in the Mitochondrial Genome Sequence of Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae)

by
Lindsey C. Perkin
1,
Timothy P. L. Smith
2 and
Brenda Oppert
3,*
1
USDA-ARS, Southern Plains Agricultural Research Center, College Station, TX 77845, USA
2
USDA-ARS, U.S. Meat Animal Research Center, Clay Center, NE 68933, USA
3
USDA-ARS, Center for Grain and Animal Health Research, Manhattan, KS 66502, USA
*
Author to whom correspondence should be addressed.
Insects 2021, 12(5), 387; https://doi.org/10.3390/insects12050387
Submission received: 5 April 2021 / Revised: 25 April 2021 / Accepted: 26 April 2021 / Published: 27 April 2021
(This article belongs to the Section Insect Molecular Biology and Genomics)

Abstract

:

Simple Summary

The lesser grain borer damages grains in storage worldwide. The major control method for this beetle is phosphine fumigation, but the increase in resistant populations has led to a loss in phosphine efficacy. Insect mitochondria are the major source of energy, and some phosphine-resistant insects have reduced energy production. Therefore, we want to understand whether changes in the mitochondrial genome may promote phosphine resistance in insects, but we need an accurate mitogenome sequence and annotation. We extracted and sequenced genomic DNA from a laboratory colony of the lesser grain borer and assembled and annotated the mitochondrial genome. The mitochondrial genome sequence was similar in structure to other insect mitochondria and encoded typical mitochondrial genes. We compared our predicted mitochondrial genome sequence to that of another lesser grain borer strain from Jingziguan (China). While there was mostly agreement among the two sequences, the data will be used to determine if key differences may suggest mutations in the two populations related to phosphine control pressure. However, differences also could be the result in different genome sequences and interpretations. The data will be useful as a research tool to examine the expression of mitochondrial genes in phosphine susceptible and -resistant insect populations.

Abstract

The lesser grain borer, Rhyzopertha dominica, is a coleopteran pest of stored grains and is mainly controlled by phosphine fumigation, but the increase in phosphine-resistant populations threatens efficacy. Some phosphine-resistant insects have reduced respiration, and thus studying the mitochondrial genome may provide additional information regarding resistance. Genomic DNA from an inbred laboratory strain of R. dominica was extracted and sequenced with both short (Illumina) and long (Pacific Biosciences) read technologies for whole genome sequence assembly and annotation. Short read sequences were assembled and annotated by open software to identify mitochondrial sequences, and the assembled sequence was manually annotated and verified by long read sequences. The mitochondrial genome sequence for R. dominica had a total length of 15,724 bp and encoded 22 trna genes, 2 rRNA genes, 13 protein coding genes (7 nad subunits, 3 cox, 2 atp, and 1 cytB), flanked by a long control region. We compared our predicted mitochondrial genome to that of another from a R. dominica strain from Jingziguan (China). While there was mostly agreement between the two assemblies, key differences will be further examined to determine if mutations in populations are related to insecticide control pressure, mainly that of phosphine. Differences in sequence data, assembly, and annotation also may result in different genome interpretations.

1. Introduction

The lesser grain borer, Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae), is a major pest of stored grains and grain products [1]. Phosphine fumigation is a major control method used in the management of infestations of stored product pests including R. dominica. Phosphine is a cost-effective fumigant, and the low cost has increased its usage, resulting in the selection and propagation of resistant populations [2].
Two nuclear genes have been identified in R. dominica, rph1 and rph2, that act synergistically to increase phosphine resistance when individuals carry both mutated genes [3,4]. Significant phosphine resistance occurs in insects with mutations in rph2, which encodes dihydrolipoamide dehydrogenase (DLD) [5]. The second gene, rph1, encodes cytochrome b5 fatty acid desaturase [6], and higher levels of resistance are found in insects with mutations in both genes.
Originally it was hypothesized that phosphine resistance was due to mutations in the mitochondrial genome [7]. The apparent association of altered respiration in phosphine-resistant insects, as has been reported in phosphine-resistant stored product insects including R. dominica [8], suggests that mitochondria also may be involved in resistance. A gene expression study indicated that genes related to mitochondrial functions were differentially expressed in phosphine-resistant Tribolium castaneum (Herbst) (Coleoptera: Tenebrionidae) [9].
Mitochondria are cellular organelles in eukaryotes that produce energy through oxidative phosphorylation (reviewed in [10]). Insect mitochondria contain a circular double-stranded DNA genome, mtDNA, that ranges in size from 15,000 to 18,000 bp [11]. Mitochondrial genomes are inherited through the female germline and their transcripts are translated using a modified genetic code [12]. Insect mitochondrial genomes are highly conserved, with 37 coding genes (13 protein-coding, two ribosomal RNA, and 22 transfer RNA genes) and a non-coding control region (CR) for initiating transcription and replication [11,13,14]. Non-protein coding RNA genes include trna (designated further with the amino acid that it transports), and rrn encoding a small mitochondrial ribosome subunit (rrnS) and large subunit (rrnL) [7]. Protein-coding genes encode polypeptides that integrate into the electron transfer chain, including complex I (6 subunits of NADH dehydrogenase, nad), complex III (cytochrome B subunit, cytB), complex IV (3 subunits of cytochrome c oxidase, cox), and complex V (2 ATP synthase subunits, atp).
During the preparation of this manuscript, the mitochondrial genome of a Chinese R. dominica strain was published [15]. The genome contained 15,862 bp with 74.36% AT content, including 13 protein-coding genes, 22 tRNA and two rRNA genes. In this study, we assembled and annotated another mitochondrial genome sequence for R. dominica from a different geographical origin. We identify differences between the mtDNA of U.S. and Chinese strains and discuss the potential relationship of the mitochondria and resistance in this stored product insect.

2. Materials and Methods

2.1. Insect Strain

The original R. dominica laboratory strain was collected from infested grain in Eastern Kansas in 1972 and has been reared at the Center for Grain and Animal Health Research in Manhattan, KS. An inbred strain from the laboratory colony of R. dominica (LGB Inbred D) was used as the source of genetic material. This strain was inbred from a single pair subculture for 20 generations, from December 2009 to March 2012.

2.2. Extraction of Nucleic Acids and Sequencing of the Mitochondrial Genome

For short read sequencing, genomic DNA was extracted from 20 mixed-sex R. dominica pupae of the LGB Inbred D strain using Quick-DNA™ Tissue/Insect Miniprep Kit (Zymo Research, Irvine, CA, USA). A portion of the DNA was sheared using a Covaris S220 for 400 base fragments as recommended by the manufacturer (Covaris Inc. Woburn, MA, USA). A library was prepared with the TruSeq® PCR-Free library preparation kit (Illumina Inc., San Diego CA, USA) and sequenced on a MiSeq instrument using a 2 × 300 base paired read v3 reagent kit. Size selection of a portion of the gDNA for long read sequencing was performed with a BluePippin instrument (Sage Science Inc., Beverly, MA, USA) using a 15kb lower cutoff value. Libraries for long read sequencing on the RSII platform were constructed using the SMRTbell™ Template Prep Kit 1.0 as recommended by the manufacturer (Pacific Biosciences, Menlo Park, CA, USA, USA). Four libraries were prepared from the same gDNA and sequenced on sixteen SMRT cells of the RSII using P5/C3 and P6/C4 chemistry (eight cells each).

2.3. Assembly and Annotation of the Mitochondrial Genome

A complete R. dominica mitochondrial genome was assembled from MiSeq short reads (accession #SRR12638565). The open-source software, NOVOPlasty version 2.6.7 [16], was used to separate the mitochondrial sequences from the nuclear sequences. This software requires only shallow sequence depth of genomic DNA and uses a genome skimming approach to extract the high-copy fraction of mitochondrial DNA from the low-copy fraction of nuclear DNA. A de novo assembly was initiated using the R. dominica cytochrome oxidase subunit 1 gene (COI; accession #KM450240). The genome was annotated using the web-based program MITOS version 1.0 [17] and visualized with GenomeVx version 1.0 [18]. Long read data were used for verification of assembled sequences (unpublished). Manual annotation of protein-coding gene sequences was done with Expasy Translate [19], BLAST [20], and Clustal Omega v. 1.2.4 [21]. The G+C content spanning the entire mitochondrial genome was assessed by an online tool (webgenetics.com, accessed on 20th September 2020).

3. Results

The short and long read sequencing produced 8.1 and 13.2 billion bases (Gb), respectively, representing approximately 17× and 58× coverage of the beetle genome. MiSeq reads corresponding to the mitochondrial genome were separated from nuclear genomic reads (see Methods), with approximately 0.39% of the 20,951,116 total reads (60,094 reads) representing putative mitochondrial reads. These reads produced a mitochondrial genome assembly (gMT) of the U.S. isolate of R. dominica with a total length of 15,724 bp and G+C content of 25.3% (accession #MW020612). Features of the gMT were annotated to identify 22 tRNA genes, 2 rRNA genes, 13 protein coding genes (7 nad subunits, 3 cox, 2 atp, and 1 cytB), and a non-coding AT-rich control region (1127 bp) (Figure 1). These data were compared to a previously reported gMT sequence of a Chinese isolate of R. dominica [15], which was slightly longer with 15,862 bp and a slightly higher G+C content of 25.6% but containing the same number of tRNA and rRNA genes (Supplemental Figure S1).

3.1. RNA Sequences

Genes encoding tRNA were the same in the two R. dominica gMT sequences, although some start codons were variable (Table 1). Single nucleotide variants (SNVs) were found only in trnaM. There was a difference in the directionality of trnaQ, which was on the reverse strand in this assembly but on the forward strand in the assembly from Ouyang et al. [15]. Predictions of tRNA secondary stem loop structures were similar to traditional tRNA (Supplemental Figure S2).
Predicted ribosomal sequences were more variable for rrnL (Table 1). The curated sequence from the Mitos prediction for rrnL had an earlier start (12,379) and stop codon (13,580), resulting in a coding sequence of 1202 vs. 1254 bp in the previous assembly. We used the MITOS prediction for rrnL, as we found that it gave higher BLAST maximum scores and lower e-values (data not shown). The sequences for rrnS only differed by a few nucleotides between the two gMT sequences. Predictions of rRNA secondary structures are provided in Supplemental Figure S3.

3.2. Protein Coding Sequences

Comparisons of protein-coding genes and gene products in the two R. dominica gMT sequences identified nonsynonymous substitutions due to SNVs, with the exception of nad4L (Supplemental Figure S4). SNVs in the other nad subunits were A-T and N-D (nad2), M-V (nad3), L-S (nad5), F-C and D-H (nad4), and L-M (nad1) (these reflect the protein predictions from our sequences versus the sequence prediction by Ouyang et al. [15]). SNVs also were found in other coding sequences: E-A (cox1), although we disagreed with the start codon and therefore there is also a difference in the N-terminus; M-T (atp8) and T-A (atp6); and V-I, M-I, M-V (cytB). The validation of these sequences was made by comparing with an alignment of the long-read sequences (Supplemental Figure S5). For every SNV, the long reads supported our predicted sequence by 13–47 long reads.
In some cases, we also found differences in N- and C-termini of predicted protein coding sequences of this assembly compared to that of Ouyang et al. [15] (Supplemental Figure S4). The coding sequences for cox1 were identical, except for a termination codon after lysine. An initiation codon of ATA has been reported for cox1 in Drosophila melanogaster [12]. For cox2, the stop codon was TCA with overlap in adjacent trnk. There also was ovelap with cox3 due to a potential stop codon within the coding sequence for trnG, but the protein coding sequence was highly similar (e−146) with cox3 from Dastarcus helophoroides (accession #YP_009049589). For atp8, we agreed with the initiation codon of ATT, but found little support in the NCBI database of other atp8 sequences for the earlier initiation codon, and many atp8 sequences have a MSPL N-terminus (data not shown). The coding sequence for nad5 was longer than the previous sequence, beginning at an earlier start codon but ending at the same stop codon. We retained the isolucene amino acid as the start codon, as well as additional amino acids YLNSL at the C-terminus of other coleopteran nad sequences (i.e., accession #ANJ70442).
The nad1 nucleotide sequences had high variability in a region where the amino acid prediction was YMNLF (Supplemental Figure S4). This region has multiple polyA and polyAT repeats, and therefore the predictions are less certain. However, we found 18 long reads to support our predicted nucleotide sequence for nad1 (Supplemental Figure S5). Our N-terminus sequence for nad1 also was slightly different than that of Ouyang et al. [15], also supported by long read data.

3.3. Control Region

The predicted control regions for the two R. dominica mitochondrial sequences were similar in the nucleotide start site (14,412 vs. 14,410) (Table 1). However, our control region length (1312 bp) was shorter than that of the comparative sequence (1449 bp) due to an additional 135 nucleotides after 14,714 bp (numbering from our sequence). There were two long reads that supported the previous sequence (Supplemental Figure S5), but the MITOS prediction did not include those 135 nucleotides, and a BLAST did not identify similar sequences in other organisms. In addition, the control regions in the two mitochondrial genome sequences had seventeen SNVs, with more variability than in other regions.
In evaluating the G+C content of the R. dominica gMT, we noticed a large increase in G+C content at the beginning of the predicted control region around 14,420 bp, followed by a sharp drop to an A+T-rich region (Supplementary Figure S6). The additional 135 nts in the control region of the previous R. dominica gMT sequence is at the transition of G+C- to A+T-rich sequences. The region with high A+T content was about 625 bp (14,850–15,475). A span of 21 polyTs starting at 14,930 nt in the A+T enriched sequence may suggest an origin of replication site at the end of the polyT run.

4. Discussion

We assembled and annotated mitochondrial gene sequences from R. dominica from short read data and used long read data to further validate sequences. The focus of our study was to understand the R. dominica mitochondrial genome in a laboratory strain from a control perspective, whereas the previous mitochondrial genome sequence of a strain from Jingziguan, China, was presented from an evolutionary biology interest [15]. We identified differences among the two assemblies and suggest that key differences from SNVs may result from nonsynonymous mutations in the two populations as a consequence of insecticide control pressure, mainly phosphine. Additionally, differences in sequence data, assembly, and annotation may result in different genome interpretations. We note that our mitochondrial genome sequence was from a highly inbred strain, and the low heterozygosity in our strain was supported by the long-read data. However, this strain has not been evaluated for resistance to phosphine.
In our assembly, trnaQ was in the reverse orientation but was forward in the assembly from Ouyang et al. [15]. The transcription orientation of mitochondrial trnas are highly conserved, and trnaQ generally has a positive orientation, but the orientation of trnaQ (-) in T. castaneum is an exception [22,23]. Thus, our annotation of the orientation of trnaQ (-) in R. dominica, as was predicted by MITOS, aligns with that exception.
Differences between the mitochondrial genome sequences of the two R. dominica strains occurred in protein-coding sequences, both in SNVs as well as interpretation of N- and C-termini. SNVs were found in all protein coding sequences from the two genome predictions, with the exception of nad4L and nad6. Differences were nonsynonymous and often resulted in amino acids with different charges or hydrophobicity that could alter the secondary structure and interactions of the subunits. SNVs may represent candidate variants that were selected by phosphine exposure, a hypothesis which could be tested in future studies by low-coverage sequencing of individuals (sufficient for genotyping of mitochondrial genomes) from phosphine-treated and phosphine-naïve locations. We also noted differences in N- and/or C-termini in cox2/cox3, atp8, and nad3/nad5 that were likely due to differences in annotation of the two gMT. High variability in the region of nad1 encoding YMNLF may suggest that the gene encoding that subunit is undergoing a higher rate of mutation. However, the poly A and AT sequences in that region also may indicate difficulty in assembly and accuracy in the sequence data.
The non-coding control region contains the longest and most variable sequence of the mitochondria [14]. Our annotated control region in the R. dominica gMT lacked a 135 nt span that was found in the sequence of Ouyang et al. [15]. The 135 nt sequence bordered the transition of a G+C enriched (>50%) sequence of approximately 250 bp to an A+T enriched sequence of approximately 625 bp. Since both gMT sequences contained this enriched G+C content, it is likely that the sequence is not an artifact, but the functional significance remains to be determined. The invertebrate control region is typified by an A+T-rich region that contains an origin of replication and transcription (Cameron, 2014). The origin of replication in T. castaneum was identified starting after a run of 17 polyT nt [24]. In our R. dominica gMT sequence, there was a span of 21 polyTs located in the middle of the A+T enriched region which also is likely an origin of replication site. However, there was no overlap in alignment of sequences from the gMT of R.dominica and T. castaneum due to high sequence variability.
A number of single point mutations in mitochondrial genes were found in phosphine-resistant Sitophilus oryzae [25]. In that paper, the authors reflected our earlier observations in T. castaneum [9], in that phosphine-resistant insects minimize energy production as a mechanism to avoid phosphine toxicity. In fact, in R. dominica adults, many mitochondrial genes are increased in expression when phosphine-susceptible insects are exposed to phosphine, whereas those genes decrease in expression in phosphine-resistant insects during phosphine exposure (unpublished data). It is not clear if the expression of mitochondrial genes is affected by mutations in previously described resistance loci [5,6]. Therefore, understanding mitochondrial genome sequence and structure, as well as gene expression, in phosphine-susceptible and -resistant insects could provide new insights into the control of resistant populations.

5. Conclusions

We sequenced, assembled, and annotated the mitochondrial genome of an inbred laboratory strain of a U.S. population of R. dominica, a problematic pest of stored grains. We compared this mitochondrial genome sequence to that of a recently published sequence from a strain from Jingziguan, China (Ouyang et al., 2019). Although the two assemblies were similar, key differences will be evaluated in the context of phosphine resistance. Additionally, differences in sequence data, assembly, and annotation may result in different genome interpretations.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12050387/s1, Figure S1: Alignment of the predicted sequence (accession #MW020612; top sequence with Rdo_MitoGenome) of the Rhyzopertha dominica mitochondrial genome and that predicted by Ouyang et al. [15] (accession #NC_042820), Figure S2: Predicted stem-loop structure for the trna genes from the mitochondrial genome of Rhyzopertha dominica, Figure S3: Predicted structures of rrna from the mitochondrial genome of Rhyzopertha dominica. (a), rrnL; (b), rrnS, Figure S4: Alignment of protein coding sequences (accession #MW020612; top sequence with Seq# in each alignment) of the Rhyzopertha dominica mitochondrial genome and that predicted by Ouyang et al. [15] (accession #NC_042820), Figure S5: Alignment of long reads (PacBio, ml#) to the predicted sequence (accession #MW020612; Rdo_MitoGenome) of the Rhyzopertha dominica mitochondrial genome and that predicted by Ouyang et al. [15] (accession #NC_042820), Figure S6: Plot of GC content (%) versus nt in the predicted sequence (accession #MW020612) of the Rhyzopertha dominica mitochondrial genome.

Author Contributions

L.C.P. and B.O. were involved in the conceptualization, methodology, and data visualization; B.O. applied a formal analysis, curation, and validation of the data; L.C.P., B.O., and T.P.L.S. contributed resources and were involved in aspects of writing, review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

LP was supported by ARS-USDA CRIS project 3091-22000-038-00-D; BO was supported by ARS-USDA CRIS project 5430-43000-025-00D; TPLS was supported by ARS-USDA CRIS Project 3040-31000-100-00D.

Institutional Review Board Statement

This research was performed in accordance with Kansas State University Research Compliance Office, Institutional Biosafety Committee registration number 1191, “Functional Genomics of Stored Product Insects”.

Data Availability Statement

The data presented in this study are openly available at the National Center for Biotechnology Information under accession numbers SRR12638565 and MW020612.

Acknowledgments

We thank Tom Morgan, Bob Lee, and Kelsey McClure for technical assistance. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. USDA is an equal opportunity provider and employer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Edde, P.A. A Review of the Biology and Control of Rhyzopertha dominica (F.) the Lesser Grain Borer. J. Stored Prod. Res. 2012, 48, 1–18. [Google Scholar] [CrossRef]
  2. Afful, E.; Elliott, B.; Nayak, M.K.; Phillips, T.W. Phosphine Resistance in North American Field Populations of the Lesser Grain Borer, Rhyzopertha dominica (Coleoptera: Bostrichidae). J. Econ. Entomol. 2018, 111, 463–469. [Google Scholar] [CrossRef]
  3. Schlipalius, D.I.; Chen, W.; Collins, P.J.; Nguyen, T.; Reilly, P.E.B.; Ebert, P.R. Gene Interactions Constrain the Course of Evolution of Phosphine Resistance in the Lesser Grain Borer, Rhyzopertha dominica. Heredity 2008, 100, 506–516. [Google Scholar] [CrossRef] [Green Version]
  4. Mau, Y.S.; Collins, P.J.; Daglish, G.J.; Nayak, M.K.; Pavic, H.; Ebert, P.R. The Rph1 Gene Is a Common Contributor to the Evolution of Phosphine Resistance in Independent Field Isolates of Rhyzopertha dominica. PLoS ONE 2012, 7, e31541. [Google Scholar] [CrossRef] [PubMed]
  5. Schlipalius, D.I.; Valmas, N.; Tuck, A.G.; Jagadeesan, R.; Ma, L.; Kaur, R.; Goldinger, A.; Anderson, C.; Kuang, J.; Zuryn, S.; et al. A Core Metabolic Enzyme Mediates Resistance to Phosphine Gas. Science 2012, 338, 807–810. [Google Scholar] [CrossRef] [PubMed]
  6. Schlipalius, D.I.; Tuck, A.G.; Jagadeesan, R.; Nguyen, T.; Kaur, R.; Subramanian, S.; Barrero, R.; Nayak, M.; Ebert, P.R. Variant Linkage Analysis Using de Novo Transcriptome Sequencing Identifies a Conserved Phosphine Resistance Gene in Insects. Genetics 2018, 209, 281–290. [Google Scholar] [CrossRef] [Green Version]
  7. Nakakita, H. The Mode of Action of Phosphine. J. Pestic. Sci. 1987, 12, 299–309. [Google Scholar] [CrossRef] [Green Version]
  8. Pimentel, M.A.G.; Faroni, L.R.D.; Tótola, M.R.; Guedes, R.N.C. Phosphine Resistance, Respiration Rate and Fitness Consequences in Stored-Product Insects. Pest Manag. Sci. 2007, 63, 876–881. [Google Scholar] [CrossRef] [Green Version]
  9. Oppert, B.; Guedes, R.N.C.; Aikins, M.J.; Perkin, L.; Chen, Z.; Phillips, T.W.; Zhu, K.Y.; Opit, G.P.; Hoon, K.; Sun, Y.; et al. Genes Related to Mitochondrial Functions Are Differentially Expressed in Phosphine-Resistant and -Susceptible Tribolium castaneum. BMC Genom. 2015, 16, 968. [Google Scholar] [CrossRef] [Green Version]
  10. Taanman, J.-W. The Mitochondrial Genome: Structure, Transcription, Translation and Replication. Biochim. Biophys. Acta BBA-Bioenerg. 1999, 1410, 103–123. [Google Scholar] [CrossRef] [Green Version]
  11. Cameron, S.L. Insect Mitochondrial Genomics: Implications for Evolution and Phylogeny. Annu. Rev. Entomol. 2014, 59, 95–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. de Bruijn, M.H.L. Drosophila melanogaster Mitochondrial DNA, a Novel Organization and Genetic Code. Nature 1983, 304, 234–240. [Google Scholar] [CrossRef]
  13. Boore, J.L. Animal Mitochondrial Genomes. Nucleic Acids Res. 1999, 27, 1767–1780. [Google Scholar] [CrossRef] [Green Version]
  14. Moritz, C.; Dowling, T.E.; Brown, W.M. Evolution of Animal Mitochondrial DNA: Relevance for Population Biology and Systematics. Annu. Rev. Ecol. Syst. 1987, 18, 269–292. [Google Scholar] [CrossRef]
  15. Ouyang, B.; Bai, Y.; Chen, J.; Wang, J.; Liang, S. Characterization of the Complete Mitochondrial Genome of Lesser Grain Borer Rhyzopertha dominica Fabricius (Insecta: Coleoptera: Bostrichidae) from Jingziguan. Mitochondrial DNA Part B 2019, 4, 3952–3953. [Google Scholar] [CrossRef] [Green Version]
  16. Dierckxsens, N.; Mardulyn, P.; Smits, G. NOVOPlasty: De Novo Assembly of Organelle Genomes from Whole Genome Data. Nucleic Acids Res. 2016, gkw955. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Bernt, M.; Donath, A.; Jühling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Pütz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de novo Metazoan Mitochondrial Genome Annotation. Mol. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef] [PubMed]
  18. Conant, G.C.; Wolfe, K.H. GenomeVx: Simple Web-Based Creation of Editable Circular Chromosome Maps. Bioinformatics 2008, 24, 861–862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Artimo, P.; Jonnalagedda, M.; Arnold, K.; Baratin, D.; Csardi, G.; de Castro, E.; Duvaud, S.; Flegel, V.; Fortier, A.; Gasteiger, E.; et al. ExPASy: SIB Bioinformatics Resource Portal. Nucleic Acids Res. 2012, 40, W597–W603. [Google Scholar] [CrossRef]
  20. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  21. Madeira, F.; Park, Y.m.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI Search and Sequence Analysis Tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Chandra, S.B.C.; Vlk, J.L.; Kapatral, V. Comparative Insect Mitochondrial Genomes: Differences despite Conserved Genome Synteny. Afr. J. Biotechnol. 2006, 5, 1308–1318. [Google Scholar]
  23. Liu, Q.-N.; Bian, D.-D.; Jiang, S.-H.; Li, Z.-X.; Ge, B.-M.; Xuan, F.-J.; Yang, L.; Li, C.-F.; Zhang, D.-Z.; Zhou, C.-L.; et al. The Complete Mitochondrial Genome of the Red Flour Beetle, Tribolium dastaneum (Coleoptera: Tenebrionidae). Mitochondrial DNA 2016, 27, 1525–1527. [Google Scholar] [CrossRef] [PubMed]
  24. Saito, S.; Tamura, K.; Aotsuka, T. Replication Origin of Mitochondrial DNA in Insects. Genetics 2005, 171, 1695–1705. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kim, K.; Yang, J.O.; Sung, J.-Y.; Lee, J.-Y.; Park, J.S.; Lee, H.-S.; Lee, B.-H.; Ren, Y.; Lee, D.-W.; Lee, S.-E. Minimization of Energy Transduction Confers Resistance to Phosphine in the Rice Weevil, Sitophilus oryzae. Sci. Rep. 2019, 9, 14605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Map of the mitochondrial genome of Rhyzopertha dominica, as predicted by Mitos. Photo used by permission (Lech Borowiec, https://www.coleoptera.org.uk/node/20008, accessed on 24 December 2020).
Figure 1. Map of the mitochondrial genome of Rhyzopertha dominica, as predicted by Mitos. Photo used by permission (Lech Borowiec, https://www.coleoptera.org.uk/node/20008, accessed on 24 December 2020).
Insects 12 00387 g001
Table 1. Comparison of Rhyzopertha dominica mitochondrial genome sequence MW020612 (this paper) and NC_042820.1. Start/stop are provided for protein-coding sequences. Bolded entries are differences among the two sequences.
Table 1. Comparison of Rhyzopertha dominica mitochondrial genome sequence MW020612 (this paper) and NC_042820.1. Start/stop are provided for protein-coding sequences. Bolded entries are differences among the two sequences.
GeneStrandFeatureStart, StopNC# (This Paper)NC_042820
StartStopLength (bp)StartStopLength (bp)
trnI+tRNA 1636306363
trnQtRNA 61129696012969
trnM+tRNA 1291976912819769
nad2+CDSATA, TAA19811909931971190993
trnW+tRNA 11891253651188125365
trnCtRNA 12461306611245130661
trnYtRNA 13061368631305136863
cox1+CDSATA, TAA137028981529136929001531
trnL2+tRNA 29012962622900296262
cox2+CDSATA, TCA2963364067829623638676
trnK+tRNA 36393709713638370971
trnD+tRNA 37093770623708377062
atp8+CDSATT, TAA3771392615637703926156
atp6+CDSATG, TAA3920458566639194585666
cox3+CDSATG, TTA4585537378945845368784
trnG+tRNA 53695430625368543062
nad3+CDSATA, TAG5431578435454305784354
trnA+tRNA 57835844625782584462
trnR+tRNA 58445908655843590865
trnN+tRNA 59085971645907597164
trnS+tRNA 59716037675970603767
trnE+tRNA 60386101646037610164
trnFtRNA 61006163646099616364
nad5CDSATT, TAA614678711726616378711708
trnHtRNA 78727934637871793463
nad4CDSATG, TTA793692521317793492521318
nad4LCDSATG, TAA9246951827392459518273
trnT+tRNA 95219582629520958262
trnPtRNA 95839645639582964563
nad6+CDSATT, TAA9647 10,129483964610,129483
cytb+CDSATG, TAG10,12911,268114010,12811,2681140
trnS2+tRNA 11,26711,3326611,26611,33266
nad1CDSATA, TAG11,35012,30395411,34912,300951
trnL1tRNA 12,30112,3626212,30012,36262
rrnLrRNA 12,37913,580120212,36213,6161254
trnVtRNA 13,61813,6796213,61613,67862
rrnSrRNA 13,68014,41173213,67814,410732
control region+ 14,41215,724131214,41015,8591449
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Perkin, L.C.; Smith, T.P.L.; Oppert, B. Variants in the Mitochondrial Genome Sequence of Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae). Insects 2021, 12, 387. https://doi.org/10.3390/insects12050387

AMA Style

Perkin LC, Smith TPL, Oppert B. Variants in the Mitochondrial Genome Sequence of Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae). Insects. 2021; 12(5):387. https://doi.org/10.3390/insects12050387

Chicago/Turabian Style

Perkin, Lindsey C., Timothy P. L. Smith, and Brenda Oppert. 2021. "Variants in the Mitochondrial Genome Sequence of Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae)" Insects 12, no. 5: 387. https://doi.org/10.3390/insects12050387

APA Style

Perkin, L. C., Smith, T. P. L., & Oppert, B. (2021). Variants in the Mitochondrial Genome Sequence of Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae). Insects, 12(5), 387. https://doi.org/10.3390/insects12050387

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop