The complete mitochondrial genomes of two rice planthoppers, Nilaparvata lugens and Laodelphax striatellus: conserved genome rearrangement in Delphacidae and discovery of new characteristics of atp8 and tRNA genes
© Zhang et al.; licensee BioMed Central Ltd. 2013
Received: 3 December 2012
Accepted: 19 June 2013
Published: 22 June 2013
Nilaparvata lugens (the brown planthopper, BPH) and Laodelphax striatellus (the small brown planthopper, SBPH) are two of the most important pests of rice. Up to now, there was only one mitochondrial genome of rice planthopper has been sequenced and very few dependable information of mitochondria could be used for research on population genetics, phylogeographics and phylogenetic evolution of these pests. To get more valuable information from the mitochondria, we sequenced the complete mitochondrial genomes of BPH and SBPH. These two planthoppers were infected with two different functional Wolbachia (intracellular endosymbiont) strains (wLug and wStri). Since both mitochondria and Wolbachia are transmitted by cytoplasmic inheritance and it was difficult to separate them when purified the Wolbachia particles, concomitantly sequencing the genome of Wolbachia using next generation sequencing method, we also got nearly complete mitochondrial genome sequences of these two rice planthoppers. After gap closing, we present high quality and reliable complete mitochondrial genomes of these two planthoppers.
The mitogenomes of N. lugens (BPH) and L. striatellus (SBPH) are 17, 619 bp and 16, 431 bp long with A + T contents of 76.95% and 77.17%, respectively. Both species have typical circular mitochondrial genomes that encode the complete set of 37 genes which are usually found in metazoans. However, the BPH mitogenome also possesses two additional copies of the trnC gene. In both mitochondrial genomes, the lengths of the atp8 gene were conspicuously shorter than that of all other known insect mitochondrial genomes (99 bp for BPH, 102 bp for SBPH). That two rearrangement regions (trnC-trnW and nad6-trnP-trnT) of mitochondrial genomes differing from other known insect were found in these two distantly related planthoppers revealed that the gene order of mitochondria might be conservative in Delphacidae. The large non-coding fragment (the A+T-rich region) putatively corresponding responsible for the control of replication and transcription of mitochondria contained a variable number of tandem repeats (VNTRs) block in different natural individuals of these two planthoppers. Comparison with a previously sequenced individual of SBPH revealed that the mitochondrial genetic variation within a species exists not only in the sequence and secondary structure of genes, but also in the gene order (the different location of trnH gene).
The mitochondrial genome arrangement pattern found in planthoppers was involved in rearrangements of both tRNA genes and protein-coding genes (PCGs). Different species from different genera of Delphacidae possessing the same mitochondrial gene rearrangement suggests that gene rearrangements of mitochondrial genome probably occurred before the differentiation of this family. After comparatively analyzing the gene order of different species of Hemiptera, we propose that except for some specific taxonomical group (e.g. the whiteflies) the gene order might have diversified in family level of this order. The VNTRs detected in the control region might provide additional genetic markers for studying population genetics, individual difference and phylogeographics of planthoppers.
The mitochondrial genome is the most commonly used molecular marker for phylogenetic studies, population genetics and dynamics, phylogeography  and even for mitochondrial related insecticide/acaricide resistance . It also can provide very favourable information at the genome level, such as the relative positions of different genes [3–7], strand asymmetry in nucleotide composition  and evolutionary patterns of the control region . Furthermore, some special structural features such as minicircles could exist in some specified groups of insects [10, 11].
The order Hemiptera is the largest group of the hemimetabolous insects  and includes three phylogenetically controversial suborders: Auchenorrhyncha, Sternorrhyncha, and Heteroptera . The complete mitochondrial genome sequences have been determined for more than 300 species of insects so far. Within Hemiptera, although 47 complete or nearly complete mitochondrial genomes are available at NCBI (as of October 1, 2012), still very few species in some groups, especially in the suborder Auchenorrhyncha (7 species) and Sternorrhyncha (9 species) had been sequenced. In the Delphacidae, only one species has had its mitochondrial genome sequenced so far . The phylogenetic relationships of this order at the superfamily and suborder levels are still controversial [14–18]. To solve these controversial phylogenetic problems in this order requires more information from mitochondrial genomes and even from nuclear genes .
The brown planthopper (BPH), Nilaparvata lugens and the small brown planthopper (SBPH), Laodelphax striatellus (Hemiptera: Delphacidae) are two of the most important pests in rice fields of many Asian countries. By transmitting virus diseases  and causing directly feeding damage, these two planthoppers cause serious yield losses of rice and economic loss for taking control measures each year . These two planthoppers were also infected with the most widespread endosymbiont Wolbachia which can manipulate the reproduction and affect the fitness of its host [23, 24]. While sequencing the genomes of two different phenotype-inducing Wolbachia strains, we also obtained nearly the complete mitochondrial genomes of BPH and SBPH. The Wolbachia strain wStri infecting SBPH can induce complete Cytoplasmic Incompatibility (CI), a reproductive modification that typically results in embryonic lethality between crosses of infected males and uninfected females. On the other hand, the Wolbachia strain wLug infecting BPH cannot induce CI at all.
After gap-closing, the complete mitochondrial genomes of these two most important rice piercing and sucking type pests, BPH and SBPH were determined. We analyzed the nucleotide composition, codon usage and genomic rearrangement of these two mitochondrial genomes and compared them with mitochondrial genomes of other hemipteran species. Some special characteristics of mitochondria of these two planthoppers (for example, the variable number of tandem repeats (VNTRs) in ‘A + T’ rich regions; location of trnH gene) were verified using natural populations. In addition, comparing mitochondrial genome of SBPH sequenced in this study with a previously sequenced individual revealed a remarkable variability between these two individuals.
Result and discussion
Genome organization, structure and composition
All of the tRNA genes have the typical cloverleaf secondary structure, except for trnS1 (AGN) of BPH, trnS2 (UCN) of SBPH which harbor a simple loop in the dihydrouridine (DHU) arm and trnH of SBPH could have a loop-stem structure like trnT of BPH (Additional file 2: Figure S1, Additional file 3: Figure S2). The dihydrouridine (DHU) arm of trnS1 (AGN) that forms a loop is found in many Hemiptera [17, 18, 27–30] and other metazoans . For SBPH, the secondary structure in the dihydrouridine (DHU) arm in the trnS1 (AGN) has a completed loop-stem structure while that of trnS2 (UCN) has a loop. In contrast, the trnS1 (AGN) has a loop and trnS2 (UCN) has a completed loop-stem structure in previously reported SBPH . Furthermore, it is similar that the trnH gene of SBPH which formed a weak loop-stem structure of the T-arm (TΨC) in this study possessed a stable loop-stem structural of the T-arm (TΨC) in previous study .
The length of all tRNA genes ranged from 56 bp for trnS2 (UCN) (SBPH) and 57 bp for trnS1 (AGN) (BPH) to 71 bp for trnK (BPH and SBPH) and trnV (BPH) (Additional file 1: Table S1 and Additional file 4: Table S2). As in other hemipteran species, BPH and SBPH harbor two rRNA genes, rrnS (748 bp for BPH and 747 bp for SBPH) and rrnL (1,219 bp for both). They are located between trnL1 (CUN) and an A+T-rich region, separated by trnV.
Mitochondrial gene order in Hemiptera
However, the arrangement of these 37 genes, especially the arrangement of large genes that encode proteins and rRNAs, is usually conserved in the suborders Auchenorrhyncha and Heteroptera, except for the unique-headed bug Stenopirates sp. . In this study, two long intergenic regions in upstream and downstream of nad6 gene which contain a poly “T” and poly “A” stretch might provide some clues of rearrangement trace. In addition to rearrangement events involving large genes, rearrangement of trnI and trnQ in Neuroctenus parus and the rearrangement of trnT and trnP in Physopelta gutta and Dysdercus cingulatus have also been reported . Except for these 4 species of Heteroptera, 2 planthoppers of Delphacidae in Auchenorrhyncha and 6 species of Aleyrodidae in the suborder Sternorrhyncha , the gene order of mitochondrial genomes of all the other species (27 species (22 families) in Heteroptera, 3 species (2 families) of Sternorrhyncha, 6 species (6 families) of Auchenorrhyncha) were found to be consistent with the mitochondrial gene order of the putative ancestral insect (Figure 3, Additional file 5: Table S3). Considering the fact that some families in the same superfamily have different gene arrangements of mitochondrial genome (Three families, Issidae, Fulgoridae and Flatidae of superfamily Fulgoroidea have the same gene order of mitochondria as the proposed ancestral insect, but members of the family Delphacidae have a different one), we proposed that excluding the whiteflies, the arrangement of genes in the mitochondria could be considered as a factor in inferring phylogenetic relationships of different species at the family level of hemipteran insects. In any case, further genome sequencing is necessary to establish whether this feature is a mitochondrial signature of the whole order.
Protein-coding genes (PCGs) and codon usage patterns
The mtDNAs of BPH and SBPH contain the full set of PCGs usually presented in animal mtDNA. Except for the nad6 gene, PCGs are arranged along the genome according to the proposed standard order of Insects  (Figure 2). All PCGs have typical start codons. Both of the putative start codons ATG and ATT are used in BPH and SBPH (Additional file 1: Table S1, Additional file 4: Table S2). Compared with SBPH, BPH has one more kind of initiation codon ATA in nad3 and atp6, while ATA was not observed in SBPH. The cox1, cox3, atp6, and nad5 genes of these two mtDNAs have incomplete stop codons. For the cox2 gene, these two planthoppers have different stop codons, TAA for BPH and TAG for SBPH (Additional file 1: Table S1, Additional file 4: Table S2). The atp8 and atp6 of BPH and SBPH are the only PCGs having different overlapping nucleotides (1 bp for BPH, 7bp for SBPH). This feature is common in hemipteran insects , and is similar to all lepidopteran mitochondrial genomes known  and many other animal mtDNAs .
It is worth mentioning that the lengths of atp8 gene of planthoppers (99 bp for BPH, 102 bp for SBPH) were the shortest in insects known so far. In fact, the longest atp8 gene (228 bp) was found in Bemisia tabaci (Hemiptera: Aleyrodidae) . In almost all other insects, the length ranged from 138 bp (Hemiptera: Psyllidae: Pachypsylla venusta)  to 183 bp (Lepidoptera: Lymantriidae: Lymantria dispar) . In most insects, atp8 gene had a length of about 160 bp. Taking into account the length of atp8 gene of planthoppers discovered in this study, it seems that the length distribution of atp8 gene in Hemiptera was wider (From 99 bp to 228 bp) than that of all other orders of insects reported so far. The sequence alignment of atp8 gene with other insects revealed that nearly 60 bp length of sequence in the 3′end of this gene was lost in these two planthoppers. Does the short length of this gene change the structure of its encoded protein? Does length change of this gene alter the functions or efficiency of this gene? These aspects are worth further studying.
Control region of mitochondrial genome of BPH and SBPH
In addition, both of these two mitochondrial genomes possess a long poly “T” stretch (23 bp) in the A+T-rich region. This poly “T” structure might be conservative in Delphacidae and be responsible for the control of replication and transcription of mitochondria.
Base composition and AT/GC-skew of mtDNA of Hemiptera
The composition of the majority strand (J-strand) of SBPH mtDNA is A = 7099 (43.2%), T = 5581 (33.97%), G = 1527 (9.29%) and C = 2224 (13.54%). The composition of the J strand of BPH mtDNA is A = 7393 (41.96%), T = 6165 (34.99%), G = 1661 (9.43%) and C = 2400 (13.62%). Similarly to most hemipteran insects, the nucleotide compositions of the two planthopper mitogenomes are significantly biased toward A and T (BPH: 76.95%; SBPH: 77.17%). Comparing with 48 mitochondrial genomes of hemipteran species, we found that there is considerable variation in base composition among different hemipteran species: ranging from 68.9% (Neuroctenus parus) to 86.3% (Aleurodicus dugesii). These two extreme A+T content species were inferring gene arrangement [32, 33]. But in different taxonomical groups of Hemiptera, the A+T% values show different ranges: Heteroptera (31 species, including 5 near complete mitochondrial genomes): 68.9-82.5%, average: 74.7±2.95%; Sternorrhyncha (9 species): 72-86.3%, average: 78.79±5.46%; Auchenorrhyncha (8 species): 75.1-78.4%, average: 76.86±0.95%. The average A+T% value for all analyzed mtDNAs set is 75.83 ± 3.67%.
The hemipteran AT-skews vary from −0.1812 (Tria-leurodes vaporariorum) to 0.2765 (Lycorma delicatula) with the BPH and SBPH mtDNA exhibiting a slight different AT-skew (0.0909 for BPH and 0.119171 for SBPH). The GC-skews of hemipteran species vary from −0.2827 (Acyrthosiphon pisum) to 0.2086 (Neomaskellia andropogonis) with the BPH and SBPH mtDNAs exhibiting very similar GC-skew (−0.1826 for BPH and −0.1842 for SBPH). Similarly, the AT/GC-skew values exhibit different ranges among three suborders: AT-skew: Heteroptera, 0.0642 - 0.2181, Sternorrhyncha, -0.1812 - 0.0956, Auchenorrhyncha, 0.0623 -T 0.2765; GC-skew: Heteroptera, -0.2656 -(−0.0965), Sternorrhyncha, -0.2810 - 0.2086, Auchenorrhyncha, -0.2827 - (−0.0783).
The mitochondrial genome is an effective data source for resolving deep-level phylogenetic problems , especially in intraordinal relationships, such as in Diptera , Hymenoptera , Orthoptera  and Hemiptera . Phylogenic analysis based on all three codon positions or first and second codon positions of PCGs can get similar results of phylogenetic relationships of Hemiptera [14, 16]. But based on the rRNA genes  or on nuclear gene combination with mitochondrial data  generated a different phylogenetic relationship. However, when comparing A+T content of all PCGs, rrnL gene and rrnS gene with the entire mitochondrial genome, PCGs show more significant positive correlation (R2 = 0.9674) than rrnL gene (R2 = 0.8873) and rrnS gene (R2 = 0.7189) (Additional file 6: Figure S3). The PCGs seem to better reflect the evolution of the entire mitochondrial genomes than the rRNA genes. Comparing with the nuclear data set, why the mitochondrial data set result in different phylogenies in Hemiptera? Actually, many factors can influence final tree : the quality of the sequences and the alignment, the amount of phylogenetic information, compositional heterogeneity, unreasonable marker choices, the accuracy of the evolutionary model and the efficiency of the tree search algorithm used. In this study, the data set of hemipteran insects show apparent base compositional heterogeneity, especially in the suborder Sternorrhyncha and Auchenorrhyncha (Figure 4). Nonstationary sequence evolution leading to unequal nucleotide composition can cause inference methods to fail and phylogenies to be inaccurate in beetles . Similarly, when we chose mitochondrial data to recover the phylogenies of Hemiptera, compositional heterogeneity of mitochondrial genomes among hemipteran species should also be considered.
Comparative analysis of two mitochondrial genomes of L. striatellus
Another significant difference of these two mitochondrial genomes of SBPH was the distinct lengths of four NADH dehydrogenase genes (nad2, nad4, nad5, nad6) and one ATP synthase gene (atp8). Because of the different locations of trnH gene, positions of the start codon of nad5 and nad6 genes and the end codon of nad4 and nad6 genes had been changed. Comparing the region (from nad5 to nad6) of mitochondrial genomes of SBPH with that of previous study, we found the differences were 57 bp deletion in the front of nad5 gene, 21 bp deletion at the end of nad4 gene and 84 bp insert in the front and 55 bp deletion at the end of nad6 gene. Furthermore, there also existed 21 bp absent in the nad2 gene and 60 bp absent at the end of atp8 gene (Additional file 9: Table S6).
Excluding five length distinct genes (nad2, nad4, nad5, nad6 and atp8) of mitochondrial genome of these two SBPH, the other eight PCGs and two ribosomal RNA genes represented significant difference in the evolution rates (Additional file 10: Figure S4). The nad4L gene which was located in the minority strand (N-strand) showed the fastest mutation rate among these genes. The mutation rates of the other PCGs ranged from 0.0022 to 0.0076. The evolution rate of rrnS gene (0.8%) was five times faster than rrnL gene (0.16%). In general, except the trnH and two trnS genes, the other tRNA genes were conserved (only few of them have single variant site).
Indeed, many insecticide resistances in arthropod pests correlated with genotype of mitochondria [2, 49–52] and membrane proteins of mitochondria  had been reported. It is not clear whether variations found in mitochondria of planthoppers are directly implicated in resistance. Furthermore, the phenomenon that insecticide selection on heteroplasmic individuals favors those carrying resistant haplotypes in natural population had been revealed in Tetranychus urticae. Whether similar mechanism referring to very rapid evolution and mutations being fixed occurred in planthoppers was also unknown. These aspects are well worth further study.
The rearrangement of mitochondrial genome of planthoppers reported in this study adds more examples of the rearranged pattern of genes in the mitochondrial genomes of hemipteran insects studied so far. The fact that the rearrangement was conserved in Delphacidae suggests that gene rearrangements of mitochondrial genome probably occurred before the differentiation of this family. Compared with other hemipterans and other insects, the extreme short lengths of atp8 genes found in BPH and SBPH and the additional repeats of trnC found in BPH raise the possibility of functional changes of these genes. The variable number of tandem repeats (VNTRs) detected in the control region might provide additional genetic markers for studying population genetics, individual difference and phylogeographics of planthoppers. The comparative analyses of two sequenced mitochondrial genomes of SBPH revealed that rearrangements of genes in the mitochondrial genome could have occurred in the same species. On account of the different locations of trnH gene, the start and end codon of its neighboring genes had been changed. Even the same tRNA gene (e.g. trnS1 (AGN) and trnS2 (UCN) genes) could form different secondary structures in these two different individuals of SBPH. All these characteristics reveal that great divergence of mitochondrial genome could occur in different individuals of a species.
Sample origin and DNA sequence assembly
Lab-reared populations of Laodelphax striatellus (Nanjing, Jiangsu province, June 2010) and Nilaparvata lugens (Sanya, Hainan province, August 2010) for mitochondrial genome sequence were collected in China. Other natural populations which were used to validating location of trnH of SBPH and detecting the variable number of tandem repeats (VNTRs) in control regions of BPH and SBPH were collected in different locations of China (Additional file 11: Table S7). Following sequencing the genome of Wolbachia infected in these two planthopper species, partial mtDNAs of the two planthoppers were also generated using Illumina HiSeq™ 2000. From the pair-end reads (read length: 90 bp) generated from the PCR-free random shotgun-sequencing libraries of average insert size 500 bp, one scaffold composed of three contigs of mitochondrial genome of SBPH (93,375 reads with approximately 520-fold coverage) and one scaffold composed of four contigs of mitochondrial genome of BPH (71,280 reads with approximately 400-fold coverage) were generated (Figure 1).
Gap closing-PCR amplification, cloning and sequencing
According to the flanking sequences of the gaps assembled from the filtered data, we designed 11 pairs of perfectly matched primers to connect the adjacent contigs (Additional file 7: Table S4). In order to reduce the error probability, we chose the high fidelity DNA polymerase Phusion® High-Fidelity DNA Polymerase (New England Biolabs) to perform the PCR amplification. The cycling conditions: 98°C 30 seconds, followed by 35 cycles of 10 s at 98°C, 30 s at different annealing temperatures, 1–2.5 min at 72°C depending on the size of amplicons, and the subsequent final elongation step at 72°C for 10 mins. The quality of PCR products was examined by 1% agarose gel electrophoresis. After purification, PCR products which were shorter than 800 bp were directly sequenced from both directions using the ABI 3730XL Genetic Analyzer (PE Applied Biosystems). Otherwise, the PCR fragments (>800 bp) were ligated into pEASY-Blunt cloning vector (Beijing TransGen Biotech) and resulting plasmid DNAs were isolated and used for sequencing. For each larger PCR product, at least two independent clones were sequenced to ensure that we obtained the consistent sequence. The flanking sequences of each gap were exactly matched with the sequences generated based on PCR method. For closing the gaps which contain the AT-rich region and detecting the variable number of tandem repeats (VNTRs) in natural population individuals, we chose the PrimeSTAR GXL DNA polymerase (TAKARA) to perform the PCR amplification. The complete sequences of mitochondrial genome of SBPH and BPH are deposited in the GenBank database under the Accession No. JX880068, JX880069 respectively.
Annotation and analysis
After finished the gap closing work, we annotated the mitochondrial genomes of BPH and SBPH. All PCGs were identified by ORF Finder implemented at the NCBI website with the invertebrate mitochondrial genetic codes. Finally pair-wise comparisons with orthologous proteins from other hemipteran insects were performed to better define the limits of PCGs. The transfer RNA (tRNA) genes were identified using both of the tRNAscan-SE  and ARWEN  programs and recognized manually as sequences having the appropriate anticodon and capable of folding into the typical cloverleaf secondary structure. The locations of the two rRNA genes (rrnL and rrnS) were determined based on alignments and secondary structures of rRNA sequences of other hemipteran species. The boundaries of the ribosomal rrnL gene were assumed to be delimited by the ends of the trnL1 (CUN) and trnV. The 3′ end of rrnS gene was assumed to be delimited by the start of trnV while the 5′ end was determined through comparison with orthologous genes of other hemipteran species sequenced so far. The base composition, the Relative Synonymous Codon Usage (RSCU) values of BPH and SBPH were calculated with MEGA 4 program .
Comparative genomics analysis
After downloaded other mitochondrioal genome sequence (Heteroptera, 31 species; Sternorrhyncha, 9 species; Auchenorrhyncha, 7 species, Additional file 5: Table S3), we compared the two planthopper species with other hemipteran species, including the gene order, the base composition and the GC/AT-skew. The GC-skew = (G-C)/(G+C) and AT-skew = (A-T)/(A+T) were used  to measure the base compositional difference between different taxa in Hemiptera. In addition, in order to reveal whether two trnH rearrangement types really exist in natural populations, we sequenced a region of mitochondrial genome encompassing the trnH gene from 309 individuals of SBPH which were collected from 15 natural populations and also checked the Wolbachia infection status of each individual (Additional file 11: Table S7). After got the sequences, we used Network 4.610 (http://www.fluxus-engineering.com) to calculate the frequency of each haplotype of the sequence of this region. The sequences of each haplotype are deposited in the GenBank database under the Accession No. KC006945 - KC006960.
We thank Dr. Ming-Long Yuan (Lanzhou University, Lanzhou, China) for helpful comments on the manuscript. We thank Dr. Hu Li (Department of Entomology, China Agricultural University, Beijing, China) for valuable communication and advice in statistical analysis. This study was supported in part by grants from the Special Fund for Agro-Scientific Research in the Public Interest (Nos. 200903051; 200803003) from the Ministry of Agriculture of China, and a grant-in-aid from the Specialized Research Fund for the Doctoral Program of Higher Education (SRFDP) from the Ministry of Education of China (Priority Development Area, No. 20110097130005).
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27 (8): 1767-1780. 10.1093/nar/27.8.1767.PubMed CentralView ArticlePubMed
- Van Leeuwen T, Vanholme B, Van Pottelberge S, Van Nieuwenhuyse P, Nauen R, Tirry L, Denholm I: Mitochondrial heteroplasmy and the evolution of insecticide resistance: non-Mendelian inheritance in action. Proc Natl Acad Sci USA. 2008, 105 (16): 5980-5985. 10.1073/pnas.0802224105.PubMed CentralView ArticlePubMed
- Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 392 (6677): 667-668. 10.1038/33577.View ArticlePubMed
- Boore JL, Brown WM: Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998, 8 (6): 668-674. 10.1016/S0959-437X(98)80035-X.View ArticlePubMed
- Boore JL: The use of genome-level characters for phylogenetic reconstruction. Trends Ecol Evol. 2006, 21 (8): 439-446. 10.1016/j.tree.2006.05.009.View ArticlePubMed
- Dowton M, Castro LR, Austin AD: Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: the examination of genome ‘morphology’. Invertebr Syst. 2002, 16 (3): 345-356. 10.1071/IS02003.View Article
- Dowton M, Austin AD: Evolutionary dynamics of a mitochondrial rearrangement “hot spot” in the hymenoptera. Mol Biol Evol. 1999, 16 (2): 298-309. 10.1093/oxfordjournals.molbev.a026111.View ArticlePubMed
- Wei SJ, Shi M, Chen XX, Sharkey MJ, van Achterberg C, Ye GY, He JH: New views on strand asymmetry in insect mitochondrial genomes. PLoS One. 2010, 5 (9): e12708-10.1371/journal.pone.0012708.PubMed CentralView ArticlePubMed
- Zhang DX, Szymura JM, Hewitt GM: Evolution and structural conservation of the control region of insect mitochondrial DNA. J Mol Evol. 1995, 40 (4): 382-391. 10.1007/BF00164024.View ArticlePubMed
- Shao RF, Kirkness EF, Barker SC: The single mitochondrial chromosome typical of animals has evolved into 18 minichromosomes in the human body louse, Pediculus humanus. Genome Res. 2009, 19 (5): 904-912. 10.1101/gr.083188.108.PubMed CentralView ArticlePubMed
- Cameron SL, Yoshizawa K, Mizukoshi A, Whiting MF, Johnson KP: Mitochondrial genome deletions and minicircles are common in lice (Insecta: Phthiraptera). BMC Genomics. 2011, 12: 394-10.1186/1471-2164-12-394.PubMed CentralView ArticlePubMed
- Schuh RT, Slater JA: True bugs of the world (Hemiptera: Heteroptera): classification and natural history. 1995, New York, the United States of America: Cornell University Press
- Carver M, Gross G, Woodward T: The Insects of Australia: a Textbook for Students and Research Workers. vol. 1. Hemiptera (bugs, leafhoppers, cicadas, aphids, scale insects etc.). 1991, Melbourne, Australia: Melbourne University Press, 429-509. 2
- Song N, Liang AP: Complete mitochondrial genome of the small brown planthopper, Laodelphax striatellus (Delphacidae: Hemiptera), with a novel gene order. Zool Sci. 2009, 26 (12): 851-860. 10.2108/zsj.26.851.View ArticlePubMed
- Cryan JR, Urban JM: Higher‒level phylogeny of the insect order Hemiptera: is Auchenorrhyncha really paraphyletic?. Syst Entomol. 2012, 37: 7-21. 10.1111/j.1365-3113.2011.00611.x.View Article
- Song N, Liang AP, Ma C: The complete mitochondrial genome sequence of the planthopper, Sivaloka damnosus. J Insect Sci. 2010, 10: 76-PubMed CentralView ArticlePubMed
- Song N, Liang A: The complete mitochondrial genome sequence of Geisha distinctissima (Hemiptera: Flatidae) and comparison with other hemipteran insects. Acta Biochim Biophys Sin (Shanghai). 2009, 41 (3): 206-216. 10.1093/abbs/gmp003.View Article
- Lee W, Kang J, Jung C, Hoelmer K, Lee SH, Lee S: Complete mitochondrial genome of brown marmorated stink bug Halyomorpha halys (Hemiptera: Pentatomidae), and phylogenetic relationships of hemipteran suborders. Mol Cells. 2009, 28 (3): 155-165. 10.1007/s10059-009-0125-9.View ArticlePubMed
- Talavera G, Vila R: What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny. BMC Evol Biol. 2011, 11 (1): 315-10.1186/1471-2148-11-315.PubMed CentralView ArticlePubMed
- Hibino H: Biology and epidemiology of rice viruses. Annu Rev Phytopathol. 1996, 34: 249-274. 10.1146/annurev.phyto.34.1.249.View ArticlePubMed
- Cheng J: Rice planthopper problems and relevant causes in China. Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia. Edited by: Heong KL, Hardy B. 2009, Los Baños, Philippines: International Rice Research Institute, 157-178.
- Hilgenboecker K, Hammerstein P, Schlattmann P, Telschow A, Werren JH: How many species are infected with Wolbachia?–A statistical analysis of current data. FEMS Microbiol Lett. 2008, 281 (2): 215-220. 10.1111/j.1574-6968.2008.01110.x.PubMed CentralView ArticlePubMed
- Zhang H, Zhang KJ, Hong XY: Population dynamics of noncytoplasmic incompatibility-inducing Wolbachia in Nilaparvata lugens and its effects on host adult life span and female fitness. Environ Entomol. 2010, 39 (6): 1801-1809. 10.1603/EN10051.View ArticlePubMed
- Noda H, Koizumi Y, Zhang Q, Deng K: Infection density of Wolbachia and incompatibility level in two planthopper species, Laodelphax striatellus and Sogatella furcifera. Insect Biochem Mol Biol. 2001, 31 (6–7): 727-737.View ArticlePubMed
- Shao R, Mitani H, Barker SC, Takahashi M, Fukunaga M: Novel mitochondrial gene content and gene arrangement indicate illegitimate inter-mtDNA recombination in the chigger mite, Leptotrombidium pallidum. J Mol Evol. 2005, 60 (6): 764-773. 10.1007/s00239-004-0226-1.View ArticlePubMed
- Jeyaprakash A, Hoy MA: The mitochondrial genome of the predatory mite Metaseiulus occidentalis (Arthropoda: Chelicerata: Acari: Phytoseiidae) is unexpectedly large and contains several novel features. Gene. 2007, 391 (1–2): 264-274.View ArticlePubMed
- Hua JM, Li M, Dong PZ, Cui Y, Xie Q, Bu WJ: Phylogenetic analysis of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha): evidence from mitochondrial genomes. BMC Evol Biol. 2009, 9: 134-10.1186/1471-2148-9-134.PubMed CentralView ArticlePubMed
- Li H, Liu H, Shi A, Stys P, Zhou X, Cai W: The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae). PLoS One. 2012, 7 (1): e29419-10.1371/journal.pone.0029419.PubMed CentralView ArticlePubMed
- Li H, Gao JY, Liu HY, Liu H, Liang AP, Zhou XG, Cai WZ: The architecture and complete sequence of mitochondrial genome of an assassin bug Agriosphodrus Dohrni (Hemiptera: Reduviidae). Int J Biol Sci. 2011, 7 (6): 792-804.PubMed CentralView ArticlePubMed
- Li H, Liu HY, Cao LM, Shi AM, Yang HL, Cai WZ: The complete mitochondrial genome of the damsel bug Alloeorhynchus bakeri (Hemiptera: Nabidae). Int J Biol Sci. 2012, 8 (1): 93-107.PubMed CentralView ArticlePubMed
- Lavrov DV, Brown WM, Boore JL: A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc Natl Acad Sci USA. 2000, 97 (25): 13738-13742. 10.1073/pnas.250402997.PubMed CentralView ArticlePubMed
- Hua J, Li M, Dong P, Cui Y, Xie Q, Bu W: Comparative and phylogenomic studies on the mitochondrial genomes of Pentatomomorpha (Insecta: Hemiptera: Heteroptera). BMC Genomics. 2008, 9: 610-10.1186/1471-2164-9-610.PubMed CentralView ArticlePubMed
- Thao ML, Baumann L, Baumann P: Organization of the mitochondrial genomes of whiteflies, aphids, and psyllids (Hemiptera, Sternorrhyncha). BMC Evol Biol. 2004, 4: 25-10.1186/1471-2148-4-25.PubMed CentralView ArticlePubMed
- Salvato P, Simonato M, Battisti A, Negrisolo E: The complete mitochondrial genome of the bag-shelter moth Ochrogaster lunifer (Lepidoptera, Notodontidae). BMC Genomics. 2008, 9: 331-10.1186/1471-2164-9-331.PubMed CentralView ArticlePubMed
- Zhu YJ, Zhou GL, Fang R, Ye J, Yi JP: The complete sequence determination and analysis of Lymantria dispar (Lepidoptera: Lymantriidae) mitochondrial genome. Plant Quarantine. 2010, 24 (4): 6-11.
- Sharp PM, Tuohy TM, Mosurski KR: Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986, 14 (13): 5125-5143. 10.1093/nar/14.13.5125.PubMed CentralView ArticlePubMed
- Dotson E, Beard CB: Sequence and organization of the mitochondrial genome of the Chagas disease vector, Triatoma dimidiata. Insect Mol Biol. 2001, 10 (3): 205-215. 10.1046/j.1365-2583.2001.00258.x.View ArticlePubMed
- Zhang DX, Hewitt GM: Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997, 25 (2): 99-120. 10.1016/S0305-1978(96)00042-7.View Article
- Cameron SL, Beckenbach AT, Dowton M, Whiting M: Evidence from mitochondrial genomics on interordinal relationships in insects. Arthropod Syst Phylogeny. 2006, 64 (1): 27-34.
- Cameron SL, Lambkin CL, Barker SC, Whiting MF: A mitochondrial genome phylogeny of Diptera: whole genome sequence data accurately resolve relationships over broad timescales with high precision. Syst Entomol. 2006, 32 (1): 40-59.View Article
- Cameron SL, Dowton M, Castro LR, Ruberu K, Whiting MF, Austin AD, Diement K, Stevens J: Mitochondrial genome organization and phylogeny of two vespid wasps. Genome. 2008, 51 (10): 800-808. 10.1139/G08-066.View ArticlePubMed
- Fenn JD, Song H, Cameron SL, Whiting MF: A preliminary mitochondrial genome phylogeny of Orthoptera (Insecta) and approaches to maximizing phylogenetic signal found within mitochondrial genome data. Mol Phylogenet Evol. 2008, 49 (1): 59-68. 10.1016/j.ympev.2008.07.004.View ArticlePubMed
- Sheffield NC, Song H, Cameron SL, Whiting MF: Nonstationary evolution and compositional heterogeneity in beetle mitochondrial phylogenomics. Syst Biol. 2009, 58 (4): 381-394. 10.1093/sysbio/syp037.View ArticlePubMed
- Baudry E, Bartos J, Emerson K, Whitworth T, Werren J: Wolbachia and genetic variability in the birdnest blowfly Protocalliphora sialia. Mol Ecol. 2003, 12 (7): 1843-1854. 10.1046/j.1365-294X.2003.01855.x.View ArticlePubMed
- Dean MD, Ballard KJ, Glass A, Ballard JW: Influence of two Wolbachia strains on population structure of East African Drosophila simulans. Genetics. 2003, 165 (4): 1959-1969.PubMed CentralPubMed
- Jiggins FM: Male-killing Wolbachia and mitochondrial DNA: selective sweeps, hybrid introgression and parasite population dynamics. Genetics. 2003, 164 (1): 5-12.PubMed CentralPubMed
- Shoemaker DD, Keller G, Ross KG: Effects of Wolbachia on mtDNA variation in two fire ant species. Mol Ecol. 2003, 12 (7): 1757-1771. 10.1046/j.1365-294X.2003.01864.x.View Article
- Shoemaker DD, Dyer KA, Ahrens M, McAbee K, Jaenike J: Decreased diversity but increased substitution rate in host mtDNA as a consequence of Wolbachia endosymbiont infection. Genetics. 2004, 168 (4): 2049-2058. 10.1534/genetics.104.030890.PubMed CentralView ArticlePubMed
- Kadir HA, Knowles CO: Inhibition of ATP dephosphorylation by acaricides with emphasis on the anti-ATPase activity of the carbodiimide metabolite of diafenthiuron. J Econ Entomol. 1991, 84 (3): 801-805.View ArticlePubMed
- Dekeyser MA: Acaricide mode of action. Pest Manag Sci. 2004, 61 (2): 103-110.View Article
- Song C, Scharf ME: Mitochondrial impacts of insecticidal formate esters in insecticide-resistant and insecticide-susceptible Drosophila melanogaster. Pest Manag Sci. 2009, 65 (6): 697-703. 10.1002/ps.1747.View ArticlePubMed
- Van Pottelberge S, Van Leeuwen T, Nauen R, Tirry L: Resistance mechanisms to mitochondrial electron transport inhibitors in a field-collected strain of Tetranychus urticae Koch (Acari: Tetranychidae). Bull Entomol Res. 2009, 99 (1): 23-31. 10.1017/S0007485308006081.View ArticlePubMed
- Guzov VM, Unnithan GC, Chernogolov AA, Feyereisen R: CYP12A1, a mitochondrial cytochrome P450 from the house fly. Arch Biochem Biophys. 1998, 359 (2): 231-240. 10.1006/abbi.1998.0901.View ArticlePubMed
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25 (5): 955-964.PubMed CentralView ArticlePubMed
- Laslett D, Canback B: ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bioinformatics. 2008, 24 (2): 172-175. 10.1093/bioinformatics/btm573.View ArticlePubMed
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMed
- Perna NT, Kocher TD: Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995, 41 (3): 353-358. 10.1007/BF01215182.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.