Drosophila melanogaster mitochondrial DNA, a novel organization and genetic code

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

The first complete mitochondrial genome sequences for Ulidiidae and phylogenetic analysis of Diptera

BACKGROUND

Tetanops sintenisi is a pest that mainly damages the root of quinoa (Chenopodium quinoa) and it is first discovered in China in 2018.

METHODS AND RESULTS

Here, the complete mitochondrial genome (mitogenome) of T. sintenisi was sequenced and compared with the mitogenomes of other Diptera species. The results revealed that the mitogenome of T. sintenisi is 15,763 bp in length (GenBank accession number: MT795181) and is comprised of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA genes, and a non-coding A + T-rich region (959 bp). The highly conserved gene arrangement of the mitogenome of T. sintenisi was identical to that of other Diptera insects. Twelve PCGs contained the typical insect start codon ATN, while cox1 had CGA as the start codon. The genes cox2, nad4, and nad1 contained an incomplete termination codon T; nad3, nad5, and cob contained the complete termination codon TAG; and the remaining seven PCGs contained the termination codon TAA. All tRNA genes were predicted to fold into the typical cloverleaf secondary structure. Phylogenetic analysis of 48 species based on the mitogenome sequence revealed that T. sintenisi clustered with the Tephritidae family, indicating that T. sintenisi and Tephritidae have a close phylogenetic relationship.

CONCLUSIONS

The phylogenetic relationship of T. sintenisi based on the mitogenome is consistent with the traditional morphological taxonomy, according to which T. sintenisi belongs to the family Otitidae, which is closely related to the family Muscidae.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae)

BACKGROUND

Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species.

AIM

 Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species.

RESULTS

We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample.

CONCLUSIONS

Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

Systematic identification of tRNA genes in Drosophila melanogaster

Transfer RNAs (tRNAs) are ubiquitous adapter molecules that link specific codons in messenger RNA (mRNA) with their corresponding amino acids during protein synthesis. The tRNA genes of Drosophila have been investigated for over half a century but have lacked systematic identification and nomenclature. Here, we review and integrate data within FlyBase and the Genomic tRNA Database (GtRNAdb) to identify the full complement of tRNA genes in the D. melanogaster nuclear and mitochondrial genomes. We apply a logical and informative nomenclature to all tRNA genes, and provide an overview of their characteristics and genomic features.

First complete mitochondrial genome of Rhodinia species (Lepidoptera: Saturniidae): genome description and phylogenetic implication

To explore the characteristics of the mitochondrial genome (mitogenome) of the squeaking silkmoths Rhodinia, a genus of wild silkmoths in the family Saturniidae of Lepidoptera, and reveal phylogenetic relationships, the mitogenome of Rhodinia fugax Butler was determined. This wild silkmoth spins a green cocoon that has potential significance in sericulture, and exhibits a unique feature that its larvae can squeak loudly when touched. The mitogenome of R. fugax is a circular molecule of 15,334 bp long and comprises 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and an A + T-rich region, consistent with previous observations of Saturniidae species. The 370-bp A + T-rich region of R. fugax contains no tandem repeat elements and harbors several features common to the Bombycidea insects, but microsatellite AT repeat sequence preceded by the ATTTA motif is not present. Mitogenome-based phylogenetic analysis shows that R. fugax belongs to Attacini, instead of Saturniini. This study presents the first mitogenome for Rhodinia genus.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

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

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.

The mitochondrial genome of Muga silkworm (Antheraea assamensis) and its comparative analysis with other lepidopteran insects

Muga (Antheraea assamensis) is an economically important silkmoth endemic to the states of Assam and Meghalaya in India and is the producer of the strongest known commercial silk. However, there is a scarcity of genomic and proteomic data for understanding the organism at a molecular level. Our present study is on decoding the complete mitochondrial genome (mitogenome) of A. assamensis using next generation sequencing technology and comparing it with other available lepidopteran mitogenomes. Mitogenome of A. assamensis is an AT rich circular molecule of 15,272 bp (A+T content ~80.2%). It contains 37 genes comprising of 13 protein coding genes (PCGs), 22 tRNA and 2 rRNA genes along with a 328 bp long control region. Its typical tRNA^Met-tRNA^Ile-tRNA^Gln arrangement differed from ancestral insects (tRNA^Ile-tRNA^Gln-tRNA^Met). Two PCGs cox1 and cox2 were found to have CGA and GTG as start codons, respectively as reported in some lepidopterans. Interestingly, nad4l gene showed higher transversion mutations at intra-species than inter-species level. All PCGs evolved under strong purifying selection with highest evolutionary rates observed for atp8 gene while lowest for cox1 gene. We observed the typical clover-leaf shaped secondary structures of tRNAs with a few exceptions in case of tRNA^Ser1 and tRNA^Tyr where stable DHU and TΨC loop were absent. A significant number of mismatches (35) were found to spread over 19 tRNA structures. The control region of mitogenome contained a six bp (CTTAGA/G) deletion atypical of other Antheraea species and lacked tandem repeats. Phylogenetic position of A. assamensis was consistent with the traditional taxonomic classification of Saturniidae. The complete annotated mitogenome is available in GenBank (Accession No. KU379695). To the best of our knowledge, this is the first report on complete mitogenome of A. assamensis.

EVOLUTION OF ASEXUALITY IN THE COSMOPOLITAN MARINE CLAM LASAEA

The marine clam genus Lasaea is unique among marine bivalves in that it contains both sexual and asexual lineages. We employed molecular tools to infer intrageneric relationships of geographically restricted sexual versus cosmopolitan asexual forms. Polymerase chain reaction primers were used to amplify and sequence homologous 624 nucleotide fragments of COIII from polyploid, asexual, direct-developing individuals representing northeastern Pacific, northeastern Atlantic, Mediterranean, southern Indian Ocean, and Australian populations. DNA sequences also were obtained from the two known diploid congeners, the Australian sexual, indirect developer, Lasaea australis, and an undescribed meiotic Australian direct developer. Estimated tree topologies did not support monophyly for polyploid asexual Lasaea lineages. A robust dichotomy was evident in all phylogenetic trees and each of the two main branches included one of the diploid meitoic Australian congeners. Lasaea australis clustered with two of the direct-developing, polyploid asexual haplotypes, one from Australia, the other from the northeastern Atlantic. Monophyly is supported for the diploid Australian direct-developing lineage together with the remaining polyploid asexual lineages from the northeastern Pacific, northeastern Atlantic, Mediterranean, and southern Indian Ocean. These results indicate that asexual Lasaea lineages are polyphyletic and may have resulted from multiple hybridization events. The high degree of genetic divergence of asexual lineages from co-clustering meiotic congeners (16%-22%) and among geographically restricted monophyletic clones (9%-11%) suggests that asexual Lasaea lineages may be exceptionally long lived.

PHYLOGENETIC ANALYSIS OF BREEDING SITE USE AND α-AMANITIN TOLERANCE WITHIN THE DROSOPHILA QUINARIA SPECIES GROUP

The Drosophila quinaria group is unusual within the genus in that it comprises both mycophagous and nonmycophagous species. DNA sequence data from three regions of the mitochondrial genome were used to infer relationships among four mycophagous species and three that breed on decaying water plants. Phylogenetic analysis of these species show that breeding in mushrooms and tolerance of high levels of α-amanitin were the ancestral states within the group. Thus, breeding in decaying water plants and intolerance of α-amanitin are derived conditions. We also found that the D. quinaria species group does not comprise separate mycophagous and nonmycophagous clades, but rather that (1) the shift from mushrooms to decaying plants occurred on at least two occasions; or (2) mycophagy reevolved within a lineage that had previously shifted to breeding on plants. The correlation between mycophagy and α-amanitin tolerance is perfect across the species we have examined, indicating that there is no detectable time lag between an ecological shift to a new breeding site and correlated changes in biochemical adaptation. The genetic distance between the mycophagous D. recens and the nonmycophagous D. quinaria indicates that these species split only about 1 M.Y.B.P. In terms of α-amanitin tolerance, D. recens and D. quinaria are typical of other ecologically similar species within the group. Thus, evolutionary changes in α-amanitin tolerance can evidently occur on the order of about 1 million yr. Our data also indicate that, in comparison to other groups of Drosophila, the quinaria species group may be undergoing an adaptive radiation.

Duplication and Remolding of tRNA Genes in the Mitochondrial Genome of Reduvius tenebrosus (Hemiptera: Reduviidae)

Most assassin bugs are predators that act as important natural enemies of insect pests. Mitochondrial (mt) genomes of these insects are double-strand circular DNAs that encode 37 genes. In the present study, we explore the duplication and rearrangement of tRNA genes in the mt genome of Reduvius tenebrosus, the first mt genome from the subfamily Reduviinae. The gene order rearranges from CR (control region)-trnI-trnQ-trnM-ND2 to CR-trnQ-trnI2-trnI1-trnM-ND2. We identified 23 tRNA genes, including 22 tRNAs commonly found in insects and an additional trnI (trnI2), which has high sequence similarity to trnM. We found several pseudo genes, such as pseudo-trnI, pseudo-CR, and pseudo-ND2, in the hotspot region of gene rearrangement (between the control region and ND2). These features provided evidence that this novel gene order could be explained by the tandem duplication/random loss (TDRL) model. The tRNA duplication/anticodon mutation mechanism further explains the presence of trnI2, which is remolded from a duplicated trnM in the TDRL process (through an anticodon mutation of CAT to GAT). Our study also raises new questions as to whether the two events proceed simultaneously and if the remolded tRNA gene is fully functional. Significantly, the duplicated tRNA gene in the mitochondrial genome has evolved independently at least two times within assassin bugs.

Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera:Culicidae), members of the coronator group

BACKGROUND

The Coronator Group currently encompasses six morphologically similar species (Culex camposi Dyar, Culex coronator Dyar and Knab, Culex covagarciai Forattini, Culex usquatus Dyar, Culex usquatissimus Dyar, and Culex ousqua Dyar). Culex coronator has been incriminated as a potential vector of West Nile Virus (WNV), Saint Louis Encephalitis Virus (SLEV), and Venezuelan Equine Encephalitis Virus (VEEV). The complete mitochondrial genome of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi was sequenced, annotated, and analyzed to provide genetic information about these species.

RESULTS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi varied from 15,573 base pairs in Cx. usquatus to 15,576 in Cx. coronator. They contained 37 genes (13 protein-encoding genes, 2 rRNA genes, and 22 tRNA genes) and the AT-rich control region. Comparative analyses of the 37 genes demonstrated the mitochondrial genomes to be composed of variable and conserved genes. Despite the small size, the ATP8, ATP6 plus NADH5 protein-encoding genes were polymorphic, whereas tRNAs and rRNAs were conserved. The control region contained some poly-T stretch. The Bayesian phylogenetic tree corroborated that both the Coronator Group and the Culex pipens complex are monophyletic taxa.

CONCLUSIONS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx. usquatissimus and Cx. camposi share the same gene composition and arrangement features that match to those reported for most Culicidae species. They are composed of the same 37 genes and the AT-rich control region, which contains poly-T stretches that may be involved in the functional role of the mitochondrial genome. Taken together, results of the dN/dS ratios, the sliding window analyses and the Bayesian phylogenetic analyses suggest that ATP6, ATP8 and NADH5 are promising genes to be employed in phylogenetic studies involving species of the Coronator Group, and probably other species groups of the subgenus Culex. Bayesian topology corroborated the morphological hypothesis of the Coronator Group as monophyletic lineage within the subgenus Culex.

The Complete Mitochondrial Genome of the Beet Webworm, Spoladea recurvalis (Lepidoptera: Crambidae) and Its Phylogenetic Implications

The complete mitochondrial genome (mitogenome) of the beet webworm, Spoladea recurvalis has been sequenced. The circular genome is 15,273 bp in size, encoding 13 protein-coding genes (PCGs), two rRNA genes, and 22 tRNA genes and containing a control region with gene order and orientation identical to that of other ditrysian lepidopteran mitogenomes. The nucleotide composition of the mitogenome shows a high A+T content of 80.9%, and the AT skewness is slightly negative (-0.023). All PCGs start with the typical ATN codons, except for COX1, which may start with the CGA codon. Nine of 13 PCGs have the common stop codon TAA; however, COX1, COX2 and ND5 utilize the T nucleotide and ND4 utilizes TA nucleotides as incomplete termination codons. All tRNAs genes are folded into the typical cloverleaf structure of mitochondrial tRNAs, except for the tRNA^Ser(AGY) gene, in which the DHU arm fails to form a stable stem-loop structure. A total of 157 bp intergenic spacers are scattered in 17 regions. The overlapping sequences are 42 bp in total and found in eight different locations. The 329 bp AT-rich region is comprised of non-repetitive sequences, including the motif ATAG, which is followed by a 14 bp poly-T stretch, a (AT₁₁ microsatellite-like repeat, which is adjacent to the motif ATTTA, and a 9 bp poly-A, which is immediately upstream from the tRNA^Met gene. Phylogenetic analysis, based on 13 PCGs and 13 PCGs+2 rRNAs using Bayesian inference and Maximum likelihood methods, show that the classification position of Pyraloidea is inconsistent with the traditional classification. Hesperioidea is placed within the Papilionoidea rather than as a sister group to it. The Pyraloidea is placed within the Macrolepidoptera with other superfamilies instead of the Papilionoidea.

The complete mitochondrial genome of Antheraea pernyi strain Yu6 (Lepidoptera: Saturniidae)

In this study, the complete mitochondrial genome of Antheraea pernyi strain Yu6 (Lepidoptera: Saturniidae) has been reported for the first time. It is a circular molecule of 15,569 bp in length, containing 37 typical coding genes and 1 non-coding AT-rich region. The overall composition of the mitogenome is A (39.27%), G (7.71%), C (12.01%), and T (41.01%). Its gene order and content are identical to the common type found in most insect mitogenomes. All protein coding genes (PCGs) start with a typical ATN initiation codon, except for the cox1 gene, which begins with TTAG codon. Nine genes used standard complete termination codon TAA, whereas the cox1, cox2, nad3, and nad5 genes end with single T. All tRNAs display typical secondary cloverleaf structures as those of other insects. Additionally, the 552 bp long AT-rich region contained several structures common to the other lepidopterons, such as some structures of repeated motifs and microsatellite-like elements. Phylogenetic analysis showed that the Antheraea pernyi Yu6 was close to Saturniidae.

The complete mitochondrial genome of Antheraea pernyi strain 731 (Lepidoptera: Saturniidae)

In this study, the complete mitochondrial genome of Antheraea pernyi strain 731 (Lepidoptera: Saturniidae) is determined for the first time. It is a circular molecule of 15,570 bp in length, with 37 typical coding genes and one non-coding A T-rich region. Its gene order and content are identical to the common type found in most insect mitogenomes. All protein-coding genes (PCGs) start with a typical ATN initiation codon, except for the cox1 gene, which begins with TTAG codon. Nine genes used standard complete termination codon TAA, whereas the cox1, cox2, nad3, and nad5 genes end with single T. All tRNAs display typical secondary cloverleaf structures as those of other insects. Additionally, the non-coding AT-rich region is 553 bp long, located between rrnS and trnM genes. It contains some structures of repeated motifs and microsatellite-like elements characteristic of the other lepidopterons. Phylogenetic analysis showed that the Antheraea pernyi 731 was close to Saturniidae.

The complete mitochondrial genome of Bombyx mori strain Yu39 (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Bombyx mori strain Yu39 (Lepidoptera: Bombycidae) is a circular molecule of 15,652 bp in length, containing 37 typical mitochondrial genes: 13 protein-coding genes (PCGs), 22 transfer RNAs, 2 ribosomal RNAs and a non-coding AT-rich region. Its gene order and arrangement are identical to the common type found in most insect mitogenomes. All PCGs start with a typical ATN codon, except for the cox1 gene, which begins with uncertained codon. All PCGs terminate in the common stop codon TAA, except for the cox1 and cox2, which use single T as their stop codons. The non-coding AT-rich region is 494-bp long, located between rrnS and trnM genes. It contains some structures of repeated motifs and microsatellite-like elements characteristic of the other lepidopterons.

The complete mitochondrial genome of Bombyx mori strain Yu5 (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Bombyx mori strain Yu5 (Lepidoptera: Bombycidae) is a circular molecule of 15,644  bp in length, containing 37 typical coding genes and one non-coding AT-rich region. The overall composition of the mitogenome was A (43.07%), G (7.34%), C (11.33%), and T (38.26%). Its gene order and content were identical to the common type found in most insect mitogenomes. All protein coding genes (PCGs) started with a typical ATN initiation codon, except for the cox1 gene, which began with CGA codon. Eleven genes used standard complete termination codon TAA, whereas the cox1 and cox2 genes ended with single T. All tRNA genes displayed typical secondary cloverleaf structures as those of other insects. Additionally, the 494  bp long AT-rich region contained several structures common to the other lepidopterons, such as some structures of repeated motifs and microsatellite-like elements and a poly-A element upstream of transfer RNA M (trnM) gene.

The complete mitochondrial genome of Bombyx mori strain Chunyun (Lepidoptera: Bombycidae)

In this study, the complete mitochondrial genome sequence of Bombyx mori strain Chunyun (Lepidoptera: Bombycidae) has been reported for the first time. The total length of the mitogenome was 15,659 bp with 37 typical animal mitochondrial genes and 1 non-coding A + T-rich region. Its gene content and order were identical to those of other lepidopteran mitochondrial genomes. All protein-coding genes (PCGs) were initiated by ATN codons except for the COI gene, which began with CGA codon. Eleven PCGs stopped with termination codon TAA, whereas the COI and COII genes ended with single T. All tRNAs have typical structures of insect mitochondrial tRNAs. The 494 bp AT-rich region contains several features common to other lepidopterans, such as some structures of repeated motifs and microsatellite-like elements and a poly-A element upstream of transfer RNA M (trnM) gene.

Sequencing of the mitochondrial genome of the avocado lace bug Pseudacysta perseae (Heteroptera, Tingidae) using a genome skimming approach

Lace bugs (Tingidae) are a family of phytophagous heteropterans, some of which are important agricultural and forestry pests. They currently comprise around 2500 species distributed worldwide, for which only one mitochondrial genome has been described so far. We sequenced the complete mitochondrial genome and the nuclear ribosomal gene segment of the avocado lace bug Pseudacysta perseae using a genome skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including the one of the sycamore lace bug Corythucha ciliata, were retrieved to allow for comparisons and phylogenetic analyses. P. perseae mitochondrial genome was determined to be 15,850 bp long, and presented the typical organisation of insect mitogenomes. The phylogenetic analysis placed P. perseae as a sister to C. ciliata but did not confirm the monophyly of Miroidae including Tingidae. Our results contradicted widely accepted phylogenetic hypothesis, which highlights the limits of analyses based on mitochondrial data only. Shotgun sequencing approaches should provide substantial improvements in harmonizing mitochondrial and nuclear databases.

The complete mitochondrial genome of Bombyx mori strain Jin6 (Lepidoptera: Bombycidae)

In this study, the complete mitochondrial genome of Bombyx mori strain Jin6 (Lepidoptera: Bombycidae) has been reported for the first time. It was a circular molecule of 15,648 bp in length, containing 37 typical coding genes and one non-coding AT-rich region. The overall composition of the mitogenome was A (43.05%), G (7.30%), C (11.35%), and T (38.29%). Its gene order and content were identical to the common type found in most insect mitogenomes. All protein coding genes (PCGs) started with a typical ATN initiation codon, except for the cox1 gene, which began with CGA codon. Eleven genes used standard complete termination codon TAA, whereas the cox1 and cox2 genes ended with single T. All tRNA genes displayed typical secondary cloverleaf structures as those of other insects. Additionally, the 494 bp long AT-rich region contained several structures common to the other lepidopterons, such as some structures of repeated motifs and microsatellite-like elements.

The complete mitochondrial genome of Bombyx mori strain Huayu (Lepidoptera: Bombycidae)

The complete mitochondrial genome of Bombyx mori strain Huayu (Lepidoptera: Bombycidae) is determined in this study. The genome was 15,666 bp long, with 37 typical animal mitochondrial genes and 1 non-coding A + T-rich region. Its gene content and order were identical to those of other lepidopteran mitochondrial genomes. All protein-coding genes (PCGs) were initiated by ATN codons except for the COI gene, which began with uncertained codon. Eleven PCGs stopped with termination codon TAA, whereas the COI and COII genes ended with single T. All tRNAs have typical structures of insect mitochondrial tRNAs. The 494 bp AT-rich region contains several features common to other lepidopterans, such as the motif ATAGA followed by an 18 bp poly-T stretch and an 11 bp poly-A element upstream of transfer RNA M (trnM) gene.

Comparative analysis of the mitochondrial genomes of Callitettixini Spittlebugs (Hemiptera: Cercopidae) confirms the overall high evolutionary speed of the AT-rich region but reveals the presence of short conservative elements at the tribal level

The present study compares the mitochondrial genomes of five species of the spittlebug tribe Callitettixini (Hemiptera: Cercopoidea: Cercopidae) from eastern Asia. All genomes of the five species sequenced are circular double-stranded DNA molecules and range from 15,222 to 15,637 bp in length. They contain 22 tRNA genes, 13 protein coding genes (PCGs) and 2 rRNA genes and share the putative ancestral gene arrangement of insects. The PCGs show an extreme bias of nucleotide and amino acid composition. Significant differences of the substitution rates among the different genes as well as the different codon position of each PCG are revealed by the comparative evolutionary analyses. The substitution speeds of the first and second codon position of different PCGs are negatively correlated with their GC content. Among the five species, the AT-rich region features great differences in length and pattern and generally shows a 2–5 times higher substitution rate than the fastest PCG in the mitochondrial genome, atp8. Despite the significant variability in length, short conservative segments were identified in the AT-rich region within Callitettixini, although absent from the other groups of the spittlebug superfamily Cercopoidea.

Shotgun assembly of the assassin bug Brontostoma colossus mitochondrial genome (Heteroptera, Reduviidae)

The complete mitochondrial genome of the assassin bug Brontostoma colossus (Distant, 1902) (Heteroptera: Reduviidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including five assassin bug species, were retrieved to allow for comparisons and phylogenetic analyses. The mitochondrial genome of B. colossus was determined to be 16,625 bp long, and consists of 13 protein-coding genes (PCGs), 23 transfer-RNA genes (tRNAs), two ribosomal-RNA genes (rRNAs), and one control region. The nucleotide composition is biased toward adenine and thymine (A+T=73.4%). Overall, architecture, nucleotide composition and genome asymmetry are similar among all available assassin bug mitogenomes. All PCGs have usual start-codons (Met and Ile). Three T and two TA incomplete termination codons were identified adjacent to tRNAs, which was consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA. All tRNAs exhibit the classic clover-leaf secondary structure except for tRNASer(AGN) in which the DHU arm forms a simple loop. Two notable features are present in the B. colossus mitogenome: (i) a 131 bp duplicated unit including the complete tRNAArg gene, resulting in 23 potentially functional tRNAs in total, and (ii) a 857 bp duplicated region comprising 277 bp of the srRNA gene and 580 bp of the control region. A phylogenetic analysis based on 55 true bug mitogenomes confirmed that B. colossus belongs to Reduviidae, but contradicted a widely accepted hypothesis. This highlights the limits of phylogenetic analyses based on mitochondrial data only.

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A+T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA(Ser)(AGN). The A. atlas A+T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.

Apocrine secretion in Drosophila salivary glands: subcellular origin, dynamics, and identification of secretory proteins

In contrast to the well defined mechanism of merocrine exocytosis, the mechanism of apocrine secretion, which was first described over 180 years ago, remains relatively uncharacterized. We identified apocrine secretory activity in the late prepupal salivary glands of Drosophila melanogaster just prior to the execution of programmed cell death (PCD). The excellent genetic tools available in Drosophila provide an opportunity to dissect for the first time the molecular and mechanistic aspects of this process. A prerequisite for such an analysis is to have pivotal immunohistochemical, ultrastructural, biochemical and proteomic data that fully characterize the process. Here we present data showing that the Drosophila salivary glands release all kinds of cellular proteins by an apocrine mechanism including cytoskeletal, cytosolic, mitochondrial, nuclear and nucleolar components. Surprisingly, the apocrine release of these proteins displays a temporal pattern with the sequential release of some proteins (e.g. transcription factor BR-C, tumor suppressor p127, cytoskeletal β-tubulin, non-muscle myosin) earlier than others (e.g. filamentous actin, nuclear lamin, mitochondrial pyruvate dehydrogenase). Although the apocrine release of proteins takes place just prior to the execution of an apoptotic program, the nuclear DNA is never released. Western blotting indicates that the secreted proteins remain undegraded in the lumen. Following apocrine secretion, the salivary gland cells remain quite vital, as they retain highly active transcriptional and protein synthetic activity.

Description of new mitochondrial genomes (Spodoptera litura, Noctuoidea and Cnaphalocrocis medinalis, Pyraloidea) and phylogenetic reconstruction of Lepidoptera with the comment on optimization schemes

We newly sequenced mitochondrial genomes of Spodoptera litura and Cnaphalocrocis medinalis belonging to Lepidoptera to obtain further insight into mitochondrial genome evolution in this group and investigated the influence of optimal strategies on phylogenetic reconstruction of Lepidoptera. Estimation of p-distances of each mitochondrial gene for available taxonomic levels has shown the highest value in ND6, whereas the lowest values in COI and COII at the nucleotide level, suggesting different utility of each gene for different hierarchical group when individual genes are utilized for phylogenetic analysis. Phylogenetic analyses mainly yielded the relationships (((((Bombycoidea + Geometroidea) + Noctuoidea) + Pyraloidea) + Papilionoidea) + Tortricoidea), evidencing the polyphyly of Macrolepidoptera. The Noctuoidea concordantly recovered the familial relationships (((Arctiidae + Lymantriidae) + Noctuidae) + Notodontidae). The tests of optimality strategies, such as exclusion of third codon positions, inclusion of rRNA and tRNA genes, data partitioning, RY recoding approach, and recoding nucleotides into amino acids suggested that the majority of the strategies did not substantially alter phylogenetic topologies or nodal supports, except for the sister relationship between Lycaenidae and Pieridae only in the amino acid dataset, which was in contrast to the sister relationship between Lycaenidae and Nymphalidae in Papilionoidea in the remaining datasets.

The complete mitochondrial genome of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of the Ailanthus silkmoth, Samia cynthia cynthia (Lepidoptera: Saturniidae) was determined. The circular genome is 15,345 bp long, and presents a typical gene organization and order for sequenced mitogenomes of Bombycidea species. The nucleotide composition of the genome is highly A+T biased, accounting for 79.86%. The AT skew of the genome is slightly negative, indicating the occurrence of more Ts than As, as found in other Saturniidae species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI and COII, which are tentatively designated by CGA and GTG, respectively, as observed in other insects. Four of 13 PCGs, including COI, COII, ATP6, and ND3, harbor the incomplete termination codons, T or TA. With an exception for tRNASer(AGN), all other tRNAs can form a typical clover-leaf structure of mitochondrial tRNA. The 359 bp A+T-rich region of S. c. cynthia contains non-repetitive sequences, but harbors several features common to the Bombycidea insects, including the motif ATAGA followed by a poly-T stretch of 19 bp, a microsatellite-like (AT)7 element preceded by the ATTTA motif, and a poly-A element upstream tRNAMet. The phylogenetic analyses support the morphology-based current hypothesis that Bombycidae and Saturniidae are monophyletic. Our result confirms that Saturniini and Attacini form a reciprocal monophyletic group within Saturniidae.

Sequence classification with side effect machines evolved via ring optimization

The explosion of available sequence data necessitates the development of sophisticated machine learning tools with which to analyze them. This study introduces a sequence-learning technology called side effect machines. It also applies a model of evolution which simulates the evolution of a ring species to the training of the side effect machines. A comparison is done between side effect machines evolved in the ring structure and side effect machines evolved using a standard evolutionary algorithm based on tournament selection. At the core of the training of side effect machines is a nearest neighbor classifier. A parameter study was performed to investigate the impact of the division of training data into examples for nearest neighbor assessment and training cases. The parameter study demonstrates that parameter setting is important in the baseline runs but had little impact in the ring-optimization runs. The ring optimization technique was also found to exhibit improved and also more reliable training performance. Side effect machines are tested on two types of synthetic data, one based on GC-content and the other checking for the ability of side effect machines to recognize an embedded motif. Three types of biological data are used, a data set with different types of immune-system genes, a data set with normal and retro-virally derived human genomic sequence, and standard and nonstandard initiation regions from the cytochrome-oxidase subunit one in the mitochondrial genome.

The complete mitochondrial genome of spittlebug Paphnutius ruficeps (Insecta: Hemiptera: Cercopidae) with a fairly short putative control region

The mitochondrial genome of the spittlebug Paphnutius ruficeps is a double-strand DNA circular molecule of 14,841 bp with a total A and T content of 73.8%. It is one of the shortest genomes among published hemipteran mitogenomes and encodes 13 protein-coding genes, 2 ribosome RNA genes and 22 transfer RNA (tRNA) genes. The gene order is consistent with the hypothesized ancestral arthropod genome arrangement. Most of the protein-coding genes use ATG as start and TAA as stop codon. The codons show an evident bias toward the nucleotides T and A at the third codon position and the most commonly used codons contain more A and T than their synonymous ones. The anticodons of the 22 tRNA genes are identical to those of the mitogenome of Philaenus spumarius, another studied spittlebug. All the tRNAs could be folded into traditional clover leaf secondary structures. The putative control region (traditionally called A + T-rich region) is the main non-coding part of the mitogenome. The AT content of this region (74.5%) is not significantly higher than that of the total mitogenome (73.8%) and slightly lower than that of the N-chain protein-coding genes (75.3%). The absence of repeat sequences as well as its short length is the most obvious characteristics of the mitochondrial genome of Paphnutius ruficeps compared with those of other published hemipteran species.

Coevolution of Drosophila melanogaster mtDNA and Wolbachia genotypes

Maternally inherited microorganisms can influence the mtDNA pattern of variation in hosts. This influence is driven by selection among symbionts and can cause the frequency of mitochondrial variants in the population to eventually increase or decrease. Wolbachia infection is common and widespread in Drosophila melanogaster populations. We compared genetic variability of D. melanogaster mitotypes with Wolbachia genotypes among isofemale lines associated with different geographic locations and time intervals to study coevolution of the mtDNA and Wolbachia. Phylogenetic analysis of D. melanogaster mtDNA revealed two clades diverged in Africa, each associated with one of the two Wolbachia genotype groups. No evidence of horizontal transmission of Wolbachia between maternal lineages has been found. All the mtDNA variants that occur in infected isofemale lines are found in uninfected isofemale lines and vice versa, which is indicative of a recent loss of infection from some maternal fly lineages and confirms a significant role of Wolbachia in the D. melanogaster mtDNA pattern of variation. Finally, we present a comparative analysis of biogeographic distribution of D. melanogaster mitotypes all over the world.

Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster

Wolbachia are maternally inherited symbiotic bacteria, commonly found in arthropods, which are able to manipulate the reproduction of their host in order to maximise their transmission. The evolutionary history of endosymbionts like Wolbachia can be revealed by integrating information on infection status in natural populations with patterns of sequence variation in Wolbachia and host mitochondrial genomes. Here we use whole-genome resequencing data from 290 lines of Drosophila melanogaster from North America, Europe, and Africa to predict Wolbachia infection status, estimate relative cytoplasmic genome copy number, and reconstruct Wolbachia and mitochondrial genome sequences. Overall, 63% of Drosophila strains were predicted to be infected with Wolbachia by our in silico analysis pipeline, which shows 99% concordance with infection status determined by diagnostic PCR. Complete Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with strict vertical transmission through the maternal cytoplasm and imperfect transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most recent common ancestor of all Wolbachia and mitochondrial genomes in D. melanogaster dates to around 8,000 years ago. We find evidence for a recent global replacement of ancestral Wolbachia and mtDNA lineages, but our data suggest that the derived wMel lineage arose several thousand years ago, not in the 20th century as previously proposed. Our data also provide evidence that this global replacement event is incomplete and is likely to be one of several similar incomplete replacement events that have occurred since the out-of-Africa migration that allowed D. melanogaster to colonize worldwide habitats. This study provides a complete genomic analysis of the evolutionary mode and temporal dynamics of the D. melanogaster–Wolbachia symbiosis, as well as important resources for further analyses of the impact of Wolbachia on host biology.

Host–microbe interactions play important roles in the physiology, development, and ecology of many organisms. Studying how hosts and their microbial symbionts evolve together over time is crucial for understanding the impact that microbes have on host biology. With the advent of high-throughput sequencing technologies, it is now possible to obtain complete genomic information for hosts and their associated microbes. Here we use whole-genome sequences from ∼300 strains of the fruitfly Drosophila melanogaster to reveal the evolutionary history of this model species and its intracellular bacterial symbiont Wolbachia. The major findings of this study are that Wolbachia in D. melanogaster is inherited strictly through the egg with no evidence of horizontal transfer from other species, that the genealogies of Wolbachia and mitochondrial genomes are virtually the same, and that both Wolbachia and mitochondrial genomes show evidence for a recent incomplete global replacement event, which has left remnant lineages in North America, Europe, and Africa. We also use the fact that Wolbachia and mitochondrial genomes have the same genealogy to estimate the rate of molecular evolution for Wolbachia, which allows us to put dates on key events in the history of this important host–microbe model system.

The complete mitochondrial genome of the pink stem borer, Sesamia inferens, in comparison with four other Noctuid moths

The complete 15,413-bp mitochondrial genome (mitogenome) of Sesamia inferens (Walker) (Lepidoptera: Noctuidae) was sequenced and compared with those of four other noctuid moths. All of the mitogenomes analyzed displayed similar characteristics with respect to gene content, genome organization, nucleotide comparison, and codon usages. Twelve-one protein-coding genes (PCGs) utilized the standard ATN, but the cox1 gene used CGA as the initiation codon; cox1, cox2, and nad4 genes had the truncated termination codon T in the S. inferens mitogenome. All of the tRNA genes had typical cloverleaf secondary structures except for trnS1(AGN), in which the dihydrouridine (DHU) arm did not form a stable stem-loop structure. Both the secondary structures of rrnL and rrnS genes inferred from the S. inferens mitogenome closely resembled those of other noctuid moths. In the A+T-rich region, the conserved motif "ATAGA" followed by a long T-stretch was observed in all noctuid moths, but other specific tandem-repeat elements were more variable. Additionally, the S. inferens mitogenome contained a potential stem-loop structure, a duplicated 17-bp repeat element, a decuplicated segment, and a microsatellite "(AT)(7)", without a poly-A element upstream of the trnM in the A+T-rich region. Finally, the phylogenetic relationships were reconstructed based on amino acid sequences of mitochondrial 13 PCGs, which support the traditional morphologically based view of relationships within the Noctuidae.

Complete mitochondrial genome of the free-living earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and phylogeny of Polyneoptera

The insect order Dermaptera, belonging to Polyneoptera, includes ∼2,000 extant species, but no dermapteran mitochondrial genome has been sequenced. We sequenced the complete mitochondrial genome of the free-living earwig, Challia fletcheri, compared its genomic features to other available mitochondrial sequences from polyneopterous insects. In addition, the Dermaptera, together with the other known polyneopteran mitochondrial genome sequences (protein coding, ribosomal RNA, and transfer RNA genes), were employed to understand the phylogeny of Polyneoptera, one of the least resolved insect phylogenies, with emphasis on the placement of Dermaptera. The complete mitochondrial genome of C. fletcheri presents the following several unusual features: the longest size in insects is 20,456 bp; it harbors the largest tandem repeat units (TRU) among insects; it displays T- and G-skewness on the major strand and A- and C-skewness on the minor strand, which is a reversal of the general pattern found in most insect mitochondrial genomes, and it possesses a unique gene arrangement characterized by a series of gene translocations and/or inversions. The reversal pattern of skewness is explained in terms of inversion of replication origin. All phylogenetic analyses consistently placed Dermaptera as the sister to Plecoptera, leaving them as the most basal lineage of Polyneoptera or sister to Ephemeroptera, and placed Odonata consistently as the most basal lineage of the Pterygota.

The mitochondrial genomes of Euphausia pacifica and Thysanoessa raschii sequenced using 454 next-generation sequencing, with a phylogenetic analysis of their position in the Malacostracan family tree

Euphausiid krill play a critical role in coastal and oceanic food webs, linking primary producers to upper trophic levels. In addition, some species support commercial fisheries worldwide. Despite their ecological importance, the genetics of these important species remain poorly described. To improve our understanding of the genetics of these ecological links, we sequenced the mitochondrial genomes of two species of North Pacific krill, Euphausia pacifica and Thysanoessa raschii, using long-range PCR and 454 GS Junior next-generation sequencing technology. The E. pacifica mitogenome (14,692 + base pairs (bp)) encodes 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and at least 22 transfer RNA (tRNA) genes. The T. raschii mitogenome (14,240 + bp) encodes 13 PCGs, two rRNA genes, and at least 19 tRNA genes. The gene order in both species is similar to that of E. superba. Comparisons between Bering Sea and Yellow Sea E. pacifica revealed a total of 644 variable sites. The most variable protein-coding gene were atp8 (7.55 %, 12 of 159 sites variable), nad4 (6.35 %, 85 variable sites) and nad6 (6.32 %, 33 variable sites). Phylogenetic analyses to assess the phylogenetic position of the Euphausiacea, using the concatenated nucleic acid sequences of E. pacifica and T. raschii along with 46 previously published malacostracan mitogenomes, support the monophyly of the order Decapoda and indicate that the Euphausiacea share a common ancestor with the Decapoda. Future research should utilize this sequence data to explore the population genetics and molecular ecology of these species.

Recovery of the mitochondrial COI barcode region in diverse Hexapoda through tRNA-based primers

Background

DNA barcoding uses a 650 bp segment of the mitochondrial cytochrome c oxidase I (COI) gene as the basis for an identification system for members of the animal kingdom and some other groups of eukaryotes. PCR amplification of the barcode region is a key step in the analytical chain, but it sometimes fails because of a lack of homology between the standard primer sets and target DNA.

Results

Two forward PCR primers were developed following analysis of all known arthropod mitochondrial genome arrangements and sequence alignment of the tRNA-W gene which was usually located within 200 bp upstream of the COI gene. These two primers were combined with a standard reverse primer (LepR1) to produce a cocktail which generated a barcode amplicon from 125 of 141 species that included representatives of 121 different families of Hexapoda. High quality sequences were recovered from 79% of the species including groups, such as scale insects, that invariably fail to amplify with standard primers.

Conclusions

A cocktail of two tRNA-W forward primers coupled with a standard reverse primer amplifies COI for most hexapods, allowing characterization of the standard barcode primer binding region in COI 5' as well as the barcode segment. The current results show that primers designed to bind to highly conserved gene regions upstream of COI will aid the amplification of this gene region in species where standard primers fail and provide valuable information to design a primer for problem groups.

Extensive variation and heteroplasmy in size of mitochondrial DNA among geographic populations of Drosophila melanogaster

Size variation and heteroplasmy in mitochondrial DNA (mtDNA) are relatively common in natural populations of Drosophila melanogaster. Of 92 isofemale lines of flies obtained from various geographic regions throughout the world, 75 lines were homoplasmic and showed a total of 12 different mtDNA size classes. The remaining 17 lines were heteroplasmic, each line carrying two different mtDNAs, and, in all but one case, the mtDNAs in these heteroplasmic lines differed in size; a total of nine size classes was represented among them. In cases where one type was predominant within an individual, it was usually the smaller mtDNA. This finding parallels what was observed in homoplasmic lines, in that the smaller mtDNAs were much more common than the larger variants in most populations. The data suggest a high rate of mutational occurrence of mtDNA size variants and some natural selection against them.

The complete mitochondrial genome of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae)

The complete mitochondrial genome (mitogenome) of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae) was determined. The genome is a circular molecule 15 481 bp long. It presents a typical gene organization and order for completely sequenced lepidopteran mitogenomes, but differs from the insect ancestral type for the placement of tRNA(Met). The nucleotide composition of the genome is also highly A + T biased, accounting for 80.38%, with a slightly positive AT skewness (0.010), indicating the occurrence of more As than Ts, as found in the Noctuoidea species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI, which is tentatively designated by the CGA codon as observed in other lepidopterans. Four of 13 PCGs harbor the incomplete termination codon, T or TA. All tRNAs have a typical clover-leaf structure of mitochondrial tRNAs, except for tRNA(Ser)(AGN), the DHU arm of which could not form a stable stem-loop structure. The intergenic spacer sequence between tRNA(Ser)(AGN) and ND1 also contains the ATACTAA motif, which is conserved across the Lepidoptera order. The H. cunea A+T-rich region of 357 bp is comprised of non-repetitive sequences, but harbors several features common to the Lepidoptera insects, including the motif ATAGA followed by an 18 bp poly-T stretch, a microsatellite-like (AT)(8) element preceded by the ATTTA motif, an 11 bp poly-A present immediately upstream tRNA(Met). The phylogenetic analyses support the view that the H. cunea is closerly related to the Lymantria dispar than Ochrogaster lunifer, and support the hypothesis that Noctuoidea (H. cunea, L. dispar, and O. lunifer) and Geometroidea (Phthonandria atrilineata) are monophyletic. However, in the phylogenetic trees based on mitogenome sequences among the lepidopteran superfamilies, Papillonoidea (Artogeia melete, Acraea issoria, and Coreana raphaelis) joined basally within the monophyly of Lepidoptera, which is different to the traditional classification.

Unique features of animal mitochondrial translation systems. The non-universal genetic code, unusual features of the translational apparatus and their relevance to human mitochondrial diseases

In animal mitochondria, several codons are non-universal and their meanings differ depending on the species. In addition, the tRNA structures that decipher codons are sometimes unusually truncated. These features seem to be related to the shortening of mitochondrial (mt) genomes, which occurred during the evolution of mitochondria. These organelles probably originated from the endosymbiosis of an aerobic eubacterium into an ancestral eukaryote. It is plausible that these events brought about the various characteristic features of animal mt translation systems, such as genetic code variations, unusually truncated tRNA and rRNA structures, unilateral tRNA recognition mechanisms by aminoacyl-tRNA synthetases, elongation factors and ribosomes, and compensation for RNA deficits by enlarged proteins. In this article, we discuss molecular mechanisms for these phenomena. Finally, we describe human mt diseases that are caused by modification defects in mt tRNAs.

The mitochondrial trn-cox1 locus: rapid evolution in Pompilidae and evidence of bias in cox1 initiation and termination codon usage

BACKGROUND

The translation initiation codon of the mitochondrial (mt) cytochrome oxidase I gene, cox1, frequently differs from conventional AUN triplets, and the termination codon is often incomplete.

METHODS

To investigate whether initiation and termination codon usage may depend on the rate of mt genome evolution, a fast evolving locus comprising cox1 and adjacent tRNA-encoding trn genes was sequenced from representatives of the major clades of spider wasps, Pompilidae (Hymenoptera), and then included in insect-wide analyses of codon usage in the context of mtDNA sequence stability.

RESULTS

The rapid evolution of Pompilidae trn-cax1 locus was particularly evident in the rearrangement of the trnC and trnY genes, the presence of a trn pseudogene in Psorthaspis, and the relatively high frequency of cox1 nucleotide substitutions. Such extensive genetic diversification, underscored by variability in the arrangement of trn genes even at the subfamily level, contrasted with the conventional form of cox1 initiation and termination codons. Similar preference for codon usage was identified for other fast evolving insect lineages unrelated to Pompilidae. Instead, conventional codons were often absent in slower evolving loci.

CONCLUSION

Analysis of codon usage across a spectrum of evolutionary rates indicates the existence of positive correlation between the incidence of cox1 canonical initiation and termination codons and relative rate of the trn-cox1 locus evolution.

Comparative genomics of Drosophila mtDNA: Novel features of conservation and change across functional domains and lineages

To gain insight on mitochondrial DNA (mtDNA) evolution, we assembled and analyzed the mitochondrial genomes of Drosophila erecta, D. ananassae, D. persimilis, D. willistoni, D. mojavensis, D. virilis and D. grimshawi together with the sequenced mtDNAs of the melanogaster subgroup. Genomic comparisons across the well-defined Drosophila phylogeny impart power for detecting conserved mtDNA regions that maintain metabolic function and regions that evolve uniquely on lineages. Evolutionary rate varies across intergenic regions of the mtDNA. Rapidly evolving intergenic regions harbor the majority of mitochondrial indel divergence. In contrast, patterns of nearly perfect conservation within intergenic regions reveal a refined set of nucleotides underlying the binding of transcription termination factors. Sequencing of 5' cDNA ends indicates that cytochrome C oxidase I (CoI) has a novel (T/C)CG start codon and that perfectly conserved regions upstream of two NADH dehydrogenase (ND) genes are transcribed and likely extend these protein sequences. Substitutions at synonymous sites in the Drosophila mitochondrial proteomes reflect a mutation process that is biased toward A and T nucleotides and differs between mtDNA strands. Differences in codon usage bias across genes reveal that weak selection at silent sites may offset the mutation bias. The mutation-selection balance at synonymous sites has also diverged between the Drosophila and Sophophora lineages. Rates of evolution are highly heterogeneous across the mitochondrial proteome, with ND accumulating many more amino acid substitutions than CO. These oxidative phosphorylation complex-specific rates of evolution vary across lineages and may reflect physiological and ecological change across the Drosophila phylogeny.

Characterization of mature mitochondrial transcripts in Drosophila, and the implications for the tRNA punctuation model in arthropods

The mitochondrial genome of Drosophila melanogaster is thought to be transcribed in about five polycistronic primary transcripts, which are processed into 11 mRNA, 22 tRNA and two rRNA species required for the genome's function. The tRNA punctuation model has been proposed to predict the cleavage sites used in this process. In this model, tRNAs are removed from the primary transcripts and the fragments that remain become mRNA and rRNA transcripts. Thus the 5' and 3' ends of the major gene transcripts are defined by the endpoints of the intervening tRNA sequences. We used 5' and 3' RACE, and circularization and RT-PCR methods to determine the sequences of both ends of all major gene transcripts of the D. melanogaster mitochondrial genome. In general, the tRNA punctuation model accurately predicts the 3' ends of most mRNA and rRNA molecules, even where there are non-coding residues present. Non-coding residues at the 5' end are evidently removed during RNA processing. The mRNAs begin precisely at the start codon for each gene. In particular, the 5' end of the cox1 gene is the first in-frame sense codon, UCG, implying that this codon serves as the start. In-frame TAA stop codons immediately preceding both the cox1 and cox2 genes may serve to prevent in-frame translation of these genes prior to the completion of processing, and are removed from the mature transcripts. Where multiple tRNA genes are present, as between nad3 and nad5, they are removed sequentially in a 3' to 5' direction.

Insect mitochondrial genomics 3: the complete mitochondrial genome sequences of representatives from two neuropteroid orders: a dobsonfly (order Megaloptera) and a giant lacewing and an owlfly (order Neuroptera)

We describe the complete mitochondrial genomes from representatives of two orders of the Neuropterida: a dobsonfly, Corydalus cornutus (Megaloptera: Corydalidae, GenBank Accession No. FJ171323), a giant lacewing Polystoechotes punctatus (Neuroptera: Polystoechotidae, FJ171325), and an owlfly, Ascaloptynx appendiculatus (Neuroptera: Ascalaphidae, FJ171324). The dobsonfly sequence is 15,687 base pairs with a major noncoding (A+T rich) region of approximately 967 bp. The gene content and organization of the dobsonfly is identical to that of most insects. The giant lacewing sequence is 16 036 bp with a major noncoding region of about 1123 bp, while the owlfly sequence is 15,877 bp with a major noncoding region of about 1066 bp. The two Neuroptera sequences include a transposition of two tRNA genes, tRNATrp and tRNACys. These tRNA genes are coded on opposite strands and overlap by seven residues in the standard insect mitochondrial gene arrangement. Thus, the transposition required a duplication of at least the region of overlap. It is likely that the transposition occurred by a duplication of both genes followed by deletion of one copy of each gene. Examination of this region in two other neuropteroid species, a snakefly, Agulla sp. (Raphidioptera: Raphidiidae), and an antlion, Myrmeleon immaculatus (Neuroptera: Myrmeleontidae), shows that the rearrangement is widespread in the order Neuroptera but not present in either of the other two orders of Neuropterida.

Comparative and phylogenomic studies on the mitochondrial genomes of Pentatomomorpha (Insecta: Hemiptera: Heteroptera)

Background

Nucleotide sequences and the gene arrangements of mitochondrial genomes are effective tools for resolving phylogenetic problems. Hemipteroid insects are known to possess highly reorganized mitochondrial genomes, but in the suborder Heteroptera (Insecta: Hemiptera), there was only one complete mitochondrial genome sequenced without gene rearrangement and the phylogeny of infraorder Pentatomomorpha in Heteroptera was still uncertain.

Results

Fifteen mitochondrial genomes of the suborder Heteroptera were sequenced. Gene rearrangements were found as follows: 1) tRNA-I and tRNA-Q switched positions in Aradidae, 2) tRNA-T and tRNA-P switched positions in Largidae and Pyrrhocoridae. Two recombination events were found in Alydidae and Malcidae. The other mt-genomes were organized in the same way as observed in Drosophila yakuba. The phylogenetic analyses of infraorder Pentatomomorpha based on the nucleotide sequence raised the hypothesis of (Aradoidea + (Pentatomoidea + (Pyrrhocoroidea + (Lygaeoidea + Coreoidea)))). The rearrangement of tRNA-T and tRNA-P also linked Largidae and Pyrrhocoridae together. Furthermore, the conserved sequence block in the unusual intergenic spacers between tRNA-H and ND4 favored the monophyly of Lygaeoidea. Tetranucleotide ATCA was inferred to be the initiation codon of ND2 in Cydnidae. No correlation was found between the rates of nucleotide substitution and gene rearrangement. CG content was significantly correlated with the nucleotide substitution rate of each gene. For ND1, there was a positive correlation (P < 0.01) between amino acids variations and hydrophobicity, but a negative correlation (P < 0.01) for ND6. No conserved sequence was found among the control regions and these regions were not always the most AT-rich region of the mt-genome.

Conclusion

Heteropteran insects are extremely complex groups worthy of further study because of the unusual tetranucleotide initiation codon and their great mt-genomic diversity, including gene rearrangements and recombinations. The mt-genome is a powerful molecular marker for resolving phylogeny at the level of the superfamily and family. Gene rearrangements were not correlated with nucleotide substitution rates. CG content variation caused the different evolutionary patterns among genes. For ND1, in many polar or nonpolar regions the specific identity of the amino acid residues might be more important than maintaining the polarity of these regions, while the opposite is true for ND6. Most sequences of the control regions did not appear to be important for regulatory functions. Finally, we suggest that the term "AT-rich regions" should not be used.