Department of Microbiology and Parasitology and Institute for Molecular Bioscience
Tropical Health Program, Austrailian Centre for International and Tropical Health and Nutrition, and Department of Zoology and Entomology, University of Queensland, Brisbane, Queensland, Australia
The mitochondrial (mt) genomes of most animals studied so far are circular, are about 16 kb in size, and have 13 protein-coding genes (atp6 and atp8 for ATP synthase subunits 6 and 8, cox1cox3 for cytochrome oxidase subunits 13, cob for cytochrome b, nad1nad6 and nad4L for NADH dehydrogenase subunits 16 and 4L), two rRNA genes (rrnL and rrnS for large and small rRNA subunits), 22 tRNA genes (one for each amino acid except for leucine and serine, which have two genes: trnL1 [anticodon sequence tag], trnL2 [taa], trnS1 [gct], and trnS2 [tga]), and one major noncoding region (Wolstenholme 1992
). The arrangement of these 37 genes, especially the arrangement of the large genes that encode proteins and rRNAs, is usually conserved in a phylum but varies substantially among phyla (Boore and Brown 1998
). However, the arrangement of protein-coding and rRNA genes has also been found to differ at lower taxonomic levels, e.g., in the family Ixodidae (Campbell and Barker 1998, 1999
) and in the genus Schistosoma (Le et al. 2000)
.
The arrangement of genes in mt genomes has been studied in more insects than any other group of invertebrates. So far, 15 species of insects have had their mt genomes sequenced completely: eight flies (order Diptera; Clary and Wolstenholme 1985
; Lewis, Farr, and Kaguni 1995
; Ballard 2000
; Lessinger et al. 2000
; Spanos et al. 2000
; Junqueira et al., unpublished but available in GenBank, accession number NC_002697); two mosquitoes (Diptera; Beard, Hamm, and Collins 1993
; Mitchell, Cockburn, and Seawright 1993
); the silkworm Bombyx mori (Lepidoptera; Lee et al., unpublished but available in GenBank, accession number AF149768); the honeybee Apis mellifera (Hymenoptera; Crozier and Crozier 1993
); the locust Locusta migratoria (Orthoptera; Flook, Rowell, and Gellissen 1995
); the bug Triatoma dimidiata (Hemiptera; Dotson and Beard 2001; accession number AF301594); and the wallaby louse Heterodoxus macropus (Phthiraptera; Shao, Campbell, and Barker 2001)
. Furthermore, part of the arrangement of mt genes is known for another 451 species of insects. The species of insects studied prior to this study were from 11 orders, but the vast majority were from four orders: Lepidoptera, Diptera, Hymenoptera, and Hemiptera. The first three orders are in the Endopterygota, whereas the Hemiptera is one of the four orders of the hemipteroid assemblage (Kristensen 1991
; fig. 1
). The arrangements of genes in the mt genomes of insects studied so far are very conserved since all species, except the wallaby louse, have the same arrangement of protein-coding and rRNA genes and most tRNA genes. Only the positions of a few tRNA genes differ, in particular, those in "hot spot" regions (Dowton and Austin 1999
), e.g., near the control region, and in the two clusters of tRNA genes, trnK-trnD and trnA-trnR-trnN-trnS1-trnE-trnF (the gene underlined here and those underlined elsewhere in this paper are transcribed from the minority strand; the two strands of the mt genomes of insects are differentiated as "majority strand" and "minority strand" according to the number of genes transcribed from them [Simon et al. 1994
]). The most common arrangement of the 37 genes in the mt genome, which is present in the fruit fly Drosophila yakuba, the bug T. dimidiata, and many other species, is inferred to be ancestral for insects (Boore, Lavrov, and Brown 1998
; Crease 1999
).
|
The wallaby louse H. macropus is the only species from the order Phthiraptera for which the gene arrangement in the mt genome has been studied. In stark contrast to all other insects studied, H. macropus has had many rearrangements of tRNA and protein-coding genes in the mt genome (Shao, Campbell, and Barker 2001)
. Prior to the present study, nothing was known about the arrangement of genes in the mt genomes of the other two orders of the hemipteroid assemblage, Thysanoptera and Psocoptera. We studied one species of Thysanoptera, seven species of Psocoptera, and another species of Hemiptera. We found that in addition to the numerous gene rearrangements in the mt genome of the wallaby louse, rearrangements of tRNA and protein-coding genes had also occurred in the Thysanoptera and the Psocoptera.
The species we studied were (1) Tectocoris diophthalmus (order Hemiptera, suborder Heteroptera, family Scutelleridae); (2) Thrips imaginis (Thysanoptera, Terebrantia, Thripidae); (3) Pteroxanium insularum (Psocoptera, Trogiomorpha, Lepidopsocidae); (4) an undescribed lepidopsocid species (Psocoptera, Trogiomorpha, Lepidopsocidae); (5) Caecilius quercus (Psocoptera, Psocomorpha, Caeciliusidae); (6) Caecilius concavistigma (Psocoptera, Psocomorpha, Caeciliusidae); (7) Clematostigma lunulata (Psocoptera, Psocomorpha, Psocidae); (8) Hemipsocus chloroticus (Psocoptera, Psocomorpha, Hemipsocidae); and (9) Pseudoscottiella tanei (Psocoptera, Psocomorpha, Pseudocaeciliidae). The cotton harlequin bug T. diophthalmus and the plague thrips T. imaginis were identified by Dianne Bell and Greg Daniels, respectively, Department of Zoology and Entomology, University of Queensland. Pteroxanium insularum was identified by Courtenay Smithers, Australian Museum, Sydney. One of us (E.R.S.) identified the six other psocopterans.
Live specimens were snap frozen in liquid nitrogen and stored at -70°C until DNA extraction. A CTAB method (Shahjahan et al. 1995
) or a DNeasy Tissue Kit (QIAGEN) was used to extract total genomic DNA. Fragments of mtDNA were amplified by PCR with the following primers: (1) C3-J-5476 and N4-N-8924 (3.5 kb), and N5-J-7567 and SR-N-14594 (7 kb) for T. diophthalmus; (2) C1-J-1718 and C3-N-5460 (3.2 kb) for T. imaginis; (3) C3-J-5476 and C2-N-3133 (5'-GATCATGAAAAAAAATTATTTGTTC-3') (616 bp) for P. insularum; (4) C1-J-2177 and C2-N-3661 (3.8 kb), and C2-J-3400 and N5-N-7736 (5'-TGTAAAATCTTATAATGCTGG-3') (2.4 kb) for the lepidopsocid species; and (5) C3-J-5357 (5'-TTTATAGCAACAGGTTTTCATGGA-3') and N5-N-7736 (1.8 kb), and N5-J-7567 and N4-N-8924 (2 kb) for C. quercus (only the sequences of primers reported here for the first time are shown; other primers and PCR conditions were reported in Shao, Campbell, and Barker 2001)
. C1-J-1718, C2-J-3400, C2-N-3661, C3-N-5460, and N4-N-8924 are conserved insect mt primers (Simon et al. 1994
); other primers were designed from the sequences we gathered and sequences in GenBank. The sizes of PCR products were estimated by electrophoresis in 1% agarose gel stained with ethidium bromide; a single, clear band was viewed from each PCR product. The PCR products from the bug T. diophthalmus were of the size expected for the ancestral arrangement of insects, whereas those from the thrips and the psocopterans were larger or smaller than expected. PCR products were purified with the QIAquick PCR Purification Kit and sequenced directly with the PCR primers and internal primers. An ABI Prism BigDye Terminator Kit was used for the sequencing reactions; sequencing products were precipitated and washed with sodium acetate and ethanol and resolved by an ABI 377 sequencer. The regions of the mt genomes reported here were sequenced from both strands for the thrips and three psocopterans, and mainly from one strand for the bug. Protein-coding and rRNA genes were identified by BLAST searches (Altschul et al. 1990
) of GenBank. tRNA genes were identified by the program tRNAscan-SE (Lowe and Eddy 1997
) or by eye. The nucleotide sequences are in GenBank under accession numbers AF335990AF335997.
Sequence analyses confirmed that the above PCR products were amplified from authentic mt targets, since (1) the base composition and the codon use were consistent with that of characterized insect mt genomes (A+T content ranged from 72% to 79% and codons with C+G at third positions range from 6% to 17%); (2) the lengths of fully sequenced cox2 and nad3 were as expected for mt genes; and (3) no frameshifts were found in the protein-coding regions when nucleotides were translated into amino acids.
Fifteen genes were identified for the bug T. diophthalmus (fig. 2A ). The arrangement of these genes is identical to the ancestral arrangement of insects. The thrips T. imaginis has the arrangement cox1-nad3-trnL2-cox2-trnG-trnK-cox3. If cox1, trnL2, cox2, trnK, and cox3 are in their ancestral positions, then three protein-coding genes (nad3, atp8, and atp6) and two tRNA genes (trnD and trnG) have changed their positions in this thrips relative to the ancestral arrangement of insects.
|
The gene rearrangements in the thrips, the psocopterans, and the wallaby louse suggest that the rate of gene rearrangement has accelerated in the mt genomes of three of the four orders of the hemipteroid insects: Thysanoptera, Psocoptera, and Phthiraptera. Of the 319 species of endopterygote insects that have been studied, only 20 species were found to have rearrangements of tRNA genes, and none had rearrangements of protein-coding or rRNA genes. Rearrangements have not been found in the 140 species of Hemiptera studied. However, we found rearrangements of tRNA and protein-coding genes in all of the Thysanoptera, Psocoptera, and Phthiraptera sequenced, and our PCR tests suggest that four other species of Psocoptera probably have gene rearrangements too. In the order Psocoptera, species from two different suborders have different arrangements of tRNA and protein-coding genes. Note that gene rearrangements have occurred not only in the "hot spot" regions of insects, but also in regions which are thought to be extremely conserved in animals, e.g., the rrnS-trnV-rrnL, nad4L-nad4-nad5, and cox1-cox2-atp8-atp6-cox3-nad3 regions.
Gene rearrangements in the mt genome provide great potential for the phylogenetic studies of the hemipteroid assemblage. Many phylogenetic relationships in the hemipteroid assemblage are contentious, e.g., whether or not the Psocoptera are monophyletic (Lyal 1985
) and which group is the sister group of the Thysanoptera (see fig. 1
). Some morphological features indicate a sister group relationship between the Thysanoptera and the Psocodea (Psocoptera plus Phthiraptera), whereas others indicate a sister group relationship between the Thysanoptera and the Hemiptera (Kristensen 1991
). The most recent analysis of the phylogeny of insect orders, from 18S, 28S, and morphological data, placed the Thysanoptera as the sister group of the Psocodea (Whiting et al. 1997
). It is intriguing that all of the Thysanoptera, Psocoptera, and Phthiraptera studied so far have had many gene rearrangements in their mt genomes, whereas the Hemiptera studied so far have the ancestral arrangement.
Why have the insects in the three hemipteroid orders had more gene rearrangements than other groups of insects? How have these rearrangements occurred? We cannot yet answer these questions. Dowton and Austin (1999)
noticed that in the Hymenoptera all mt gene rearrangements occurred after the parasitic lifestyle was adopted. The correlation between frequency of gene rearrangement and parasitism is supported by the numerous gene rearrangements in the mt genome of the wallaby louse, which is an obligate ectoparasite. In the well-studied order Diptera, gene rearrangements have been found only in three species of mosquitoes, which are also ectoparasites. Furthermore, inversions of mt genes in insects have been found only in mosquitoes, the wallaby louse, and species of Hymenoptera that evolved after the parasitic lifestyle was adopted. On the other hand, however, rearrangements of tRNA and protein-coding genes have been found in free-living insects like the thrips and the psocopterans studied by us.
Unlike the gene rearrangements of the wallaby louse, in which both translocations and inversions have occurred, only translocations were found in the thrips and the psocopterans studied here. Tandem duplications of a region followed by deletions of redundant copies of genes (Macey et al. 1997
) is a plausible mechanism for the rearrangements found in the thrips and the psocopterans. However, it is not clear how the translocations and inversions of genes in the mt genome of the wallaby louse occurred. The mechanism of tandem duplication and deletion can hardly account for the drastic translocations of genes and is not applicable to the inversions of genes. Lunt and Hyman (1997)
reported an intramitochondrial recombination in a nematode: there was double-strand breakage at both ends of a fragment in the major noncoding region, after which rejoining was indicated by the presence of a subgenomic minicircle about 250 bp long. This mechanism could account for small inversions, especially those that have occurred in the noncoding region. Whether or not fragments in a coding region, particularly large genes like nad1 (894 bp) and atp6 (657 bp) of the wallaby louse, could invert in this way awaits further studies. Comparisons between fully characterized mt genomes of a range of hemipterans, thrips, psocopterans, and lice should provide further insight into the mechanisms of gene rearrangement in animal mt genomes.
Acknowledgements
We thank Ellen Dotson and Charles Beard for sharing their data on T. dimidiata before it was published; Courtenay Smithers for collecting and identifying P. insularum; Steven Cameron for collecting and Greg Daniels for identifying T. imaginis; and Dianne Bell for collecting and identifying T. diophthalmus.
Footnotes
Ross Crozier, Reviewing Editor
1 Keywords: Paraneoptera
gene order
gene arrangement
mtDNA
evolution
phylogeny
2 Address for correspondence and reprints: Renfu Shao, Department of Microbiology and Parasitology, University of Queensland, Brisbane, Queensland, 4072, Australia. r.shao{at}mailbox.uq.edu.au
References
Altschul S. F., W. Gish, W. Miller, E. W. Myers, D. J. Lipman, 1990 Basic local alignment search tool J. Mol. Biol 215:403-410[ISI][Medline]
Ballard J. W. O., 2000 Comparative genomics of mitochondrial DNA in members of the Drosophila melanogaster subgroup J. Mol. Evol 51:48-63[ISI][Medline]
Beard C. B., D. M. Hamm, F. H. Collins, 1993 The mitochondrial genome of the mosquito Anopheles gambiae: DNA sequence, genome organization, and comparisons with mitochondrial sequences of other insects Insect Mol. Biol 2:103-104[Medline]
Boore J. L., W. M. Brown, 1998 Big trees from little genomes: mitochondrial gene order as a phylogenetic tool Curr. Opin. Genet. Dev 8:668-674[ISI][Medline]
Boore J. L., D. V. Lavrov, W. M. Brown, 1998 Gene translocation links insects and crustaceans Nature 392:667-668[ISI][Medline]
Campbell N. J. H., S. C. Barker, 1998 An unprecedented major rearrangement in an arthropod mitochondrial genome Mol. Biol. Evol 15:1786-1787
. 1999 The novel mitochondrial gene arrangement of the cattle tick, Boophilus microplus: fivefold tandem repetition of a coding region Mol. Biol. Evol 16:732-740[Abstract]
Clary D. O., D. R. Wolstenholme, 1985 The mitochondrial DNA molecule of Drosophila yakuba: nucleotide sequence, gene organization, and genetic code J. Mol. Evol 22:252-271[ISI][Medline]
Crease T. J., 1999 The complete sequence of the mitochondrial genome of Daphnia pulex (Cladocera: Crustacea) Gene 233:89-99[ISI][Medline]
Crozier R. H., Y. C. Crozier, 1993 The mitochondrial genome of the honeybee Apis mellifera: complete sequence and genome organization Genetics 133:97-117
Dotson E. M., C. B. Beard, 2001 Sequence and organization of the mitochondrial genome of the Chagas disease vector, Triatoma dimidiata Insect Mol. Biol. (in press).
Dowton M., A. D. Austin, 1999 Evolutionary dynamics of a mitochondrial rearrangement "hot spot" in the Hymenoptera Mol. Biol. Evol 16:298-309[Abstract]
Flook P. K., C. H. F. Rowell, G. Gellissen, 1995 The sequence, organisation and evolution of the Locusta migratoria mitochondrial genome J. Mol. Evol 41:928-941[ISI][Medline]
Kristensen N. P., 1991 Phylogeny of extant hexapods P. 135 in CSIRO, eds. The insects of Australia: a textbook for students and research workers. Melbourne University Press, Carlton, Victoria
Le T. H., D. Blair, T. Agarsuma, et al. (13 co-authors) 2000 Phylogenies inferred from mitochondrial gene ordersa cautionary tale from the parasitic flatworms Mol. Biol. Evol 17:1123-1125
Lessinger A. C., A. C. M. Junqueira, T. A. Lemos, E. L. Kemper, F. R. da Silva, A. L. Vettore, P. Arruda, A. M. L. Azeredo-Espin, 2000 The mitochondrial genome of the primary screwworm fly Cochliomyia hominivorax (Diptera: Calliphoridae) Insect Mol. Biol 9:521-529[ISI][Medline]
Lewis D. L., C. L. Farr, L. S. Kaguni, 1995 Drosophila melanogaster mitochondrial DNA: completion of the nucleotide sequence and evolutionary comparisons Insect Mol. Biol 4:263-278[ISI][Medline]
Lowe T. M., S. R. Eddy, 1997 tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence Nucleic Acids Res 25:955-964
Lunt D. H., B. C. Hyman, 1997 Animal mitochondrial DNA recombination Nature 387:247[ISI][Medline]
Lyal C. H. C., 1985 Phylogeny and classification of the Psocodea, with particular reference to the lice (Psocodea: Phthiraptera) Syst. Entomol 10:145-165[ISI]
Macey J. R., A. Larson, N. B. Anajeva, Z. Fang, T. J. Papenfuss, 1997 Two novel gene orders and the role of light-strand replication in the rearrangement of the vertebrate mitochondrial genome Mol. Biol. Evol 14:91-104[Abstract]
Mitchell S. E., A. F. Cockburn, J. A. Seawright, 1993 The mitochondrial genome of Anopheles quadrimaculatus species A: complete nucleotide sequence and gene organization Genome 36:1058-1073[ISI][Medline]
Shahjahan R. M, K. J. H. Roger, R. A. Leopold, J. D. Devault, 1995 Lower incubation-temperature increases yield of insect genomic DNA isolated by the CTAB method Biotechniques 19:332[ISI][Medline]
Shao R., N. J. H. Campbell, S. C. Barker, 2001 Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera) Mol. Biol. Evol 18:858-865
Simon C., F. Frati, A. Beckenbach, B. Crespi, H. Liu, P. Flook, 1994 Evolution, weighting and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers Ann. Entomol. Soc. Am 87:651-701[ISI]
Sorensen J. T., B. C. Campbell, R. J. Gill, J. D. Steffen-Campbell, 1995 Non-monophyly of Auchenorrhyncha ("Homoptera"), based upon 18S rDNA phylogeny: eco-evolutionary and cladistic implications within pre-Heteropterodea Hemiptera (S. L.) and a proposal for new monophyletic suborders Pan-Pacific Entomol 71:31-60
Spanos L., G. Koutroumbas, M. Kotsyfakis, C. Louis, 2000 The mitochondrial genome of the Mediterranean fruit fly, Ceratitis capitata Insect Mol. Biol 9:139-144[ISI][Medline]
Whiting M. F., J. C. Carpenter, Q. D. Wheeler, W. C. Wheeler, 1997 The Strepsiptera problem: phylogeny of the holometabolous insect orders inferred from 18S and 28S ribosomal DNA sequences and morphology Syst. Biol 46:1-68[ISI][Medline]
Wolstenholme D. R., 1992 Animal mitochondrial DNA: structure and evolution Int. Rev. Cytol 141:173-216[ISI][Medline]