The mitochondrial DNA control region of Muscidae flies: evolution and structural conservation in a dipteran context

Features and evolution of control regions in leafroller moths (Lepidoptera: Tortricidae) inferred from mitochondrial genomes and phylogeny

The control region (CR) of the mitochondrial genome (mitogenome) represents a major noncoding fragment with several special structural features that are thought to be responsible for the initiation of mitogenome transcription and replication. However, few studies have revealed the evolutionary patterns of CR in the phylogenetic context. Here, we explain the characteristics and evolution of CR in Tortricidae, inferred from a mitogenome-based phylogeny. The first complete mitogenomes of the genera Meiligma and Matsumuraeses were sequenced. Both mitogenomes are double-stranded circular DNA molecules with lengths of 15,675 bp and 15,330 bp, respectively. Phylogenetic analyses derived from 13 protein-coding genes and two ribosomal RNAs showed that most tribes, including subfamilies Olethreutinae and Tortricinae, were recovered as monophyletic clades, similar to previous studies based on morphological or nuclear data. Moreover, comprehensive comparative analyses of the structural organization and role of tandem replications on the length variation and high AT content of CR sequences were conducted. The results reveal a significant positive correlation between the total length and AT content of tandem repeats and whole CR sequences in Tortricidae. The structural organization in CR sequences is diverse, even between closely related tribes, which demonstrates the plasticity of the mitochondrial DNA molecule in Tortricidae.

Characterization of the complete mitochondrial genome of Orthaga olivacea Warre (Lepidoptera Pyralidae) and comparison with other Lepidopteran insects

Orthaga olivacea Warre (Lepidoptera: Pyralidae) is an important agricultural pest of camphor trees (Cinnamomum camphora). To further supplement the known genome-level features of related species, the complete mitochondrial genome of Orthaga olivacea is amplified, sequenced, annotated, analyzed, and compared with 58 other species of Lepidopteran. The complete sequence is 15,174 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a putative control region. Base composition is biased toward adenine and thymine (79.02% A+T) and A+T skew are slightly negative. Twelve of the 13 PCGs use typical ATN start codons. The exception is cytochrome oxidase 1 (cox1) that utilizes a CGA initiation codon. Nine PCGs have standard termination codon (TAA); others have incomplete stop codons, a single T or TA nucleotide. All the tRNA genes have the typical clover-leaf secondary structure, except for trnS(AGN), in which dihydrouridine (DHU) arm fails to form a stable stem-loop structure. The A+T-rich region (293 bp) contains a typical Lepidopter motifs 'ATAGA' followed by a 17 bp poly-T stretch, and a microsatellite-like (AT)13 repeat. Codon usage analysis revealed that Asn, Ile, Leu2, Lys, Tyr and Phe were the most frequently used amino acids, while Cys was the least utilized. Phylogenetic analysis suggested that among sequenced lepidopteran mitochondrial genomes, Orthaga olivacea Warre was most closely related to Hypsopygia regina, and confirmed that Orthaga olivacea Warre belongs to the Pyralidae family.

Characterization of the complete mitogenome of Anopheles aquasalis, and phylogenetic divergences among Anopheles from diverse geographic zones

Whole mitogenome sequences (mtDNA) have been exploited for insect ecology studies, using them as molecular markers to reconstruct phylogenies, or to infer phylogeographic relationships and gene flow. Recent Anopheles phylogenomic studies have provided information regarding the time of deep lineage divergences within the genus. Here we report the complete 15,393 bp mtDNA sequences of Anopheles aquasalis, a Neotropical human malaria vector. When comparing its structure and base composition with other relevant and available anopheline mitogenomes, high similarity and conserved genomic features were observed. Furthermore, 22 mtDNA sequences comprising anopheline and Dipteran sibling species were analyzed to reconstruct phylogenies and estimate dates of divergence between taxa. Phylogenetic analysis using complete mtDNA sequences suggests that A. aquasalis diverged from the Anopheles albitarsis complex ~28 million years ago (MYA), and ~38 MYA from Anopheles darlingi. Bayesian analysis suggests that the most recent ancestor of Nyssorhynchus and Anopheles + Cellia was extant ~83 MYA, corroborating current estimates of ~79–100 MYA. Additional sampling and publication of African, Asian, and North American anopheline mitogenomes would improve the resolution of the Anopheles phylogeny and clarify early continental dispersal routes.

Comparative analysis of mitochondrial genomes among four species of muscid flies (Diptera: Muscidae) and its phylogenetic implications

Muscidae, commonly known as house flies and their close relatives, is one of the dipteran insects of recognized medical, veterinary, and ecological importance. Mitochondrial genomes (Mitogenomes) have been widely used for exploring phylogenetic analysis and taxonomic diagnosis due to the difficulty in distinguishing them morphologically. In this study, our complete mitogenomes of muscid flies were sequenced and aligned, which ranged from 15,117 bp (Synthesiomyia nudiseta) to 16,089 bp (Musca sorbens) in length, and contained a typical circular molecule comprising 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and a non-coding control region. The order and orientation of genes were identical with that from the ancestral insects. The phylogenetic analysis based on the mitochondrial genes indicated that the subfamily relationships within Muscidae were reconstructed as (Mydaeinae (Muscinae (Reinwardtiinae + Azeliinae))). Similar tree topologies were recovered from both Maximum Likelihood (ML) and Bayesian Inference (BI) analysis. Furthermore, we compared the phylogenetic analyses that were constructed using internal transcribed spacer 2 (ITS2), elongation factor-1α (EF-1α), 13 PCGs and 13 PCGs + ITS2 + EF-1α, respectively. Combined analysis of nuclear gene partitions improved support and resolution for resulting topologies but the positions of branches were obviously inconsistent due to limited species. More mitogenomes should be sequenced representing various taxonomic levels, especially close related species, which will enhance our understanding of phylogenetic relationships among muscids.

Hemiptera Mitochondrial Control Region: New Sights into the Structural Organization, Phylogenetic Utility, and Roles of Tandem Repetitions of the Noncoding Segment

As a major noncoding fragment, the control region (CR) of mtDNA is responsible for the initiation of mitogenome transcription and replication. Several structural features of CR sequences have been reported in many insects. However, comprehensive analyses on the structural organization and phylogenetic utility, as well as the role of tandem replications (TRs) on length variation, high A+T content, and shift of base skew of CR sequences are poorly investigated in hemipteran insects. In this study, we conducted a series of comparative analyses, using 116 samples covering all 11 infraorders of the five currently recognized monophyletic groups in the Hemiptera. Several structural elements (mononucleotide stretches containing conserved sequence blocks (CSBs), TRs, and GA-rich region) were identified in the mitochondrial control region in hemipteran insects, without showing a consistent location. The presence and absence of certain specific structural elements in CR sequences show the various structural organizations of that segment among the five monophyletic groups, which indicates the diversification of the control region’s structural organization in Hemiptera. Among the many groups within Hemiptera, eight monophyletic groups and three consistent phylogenetic trees were recovered, using CSBs datasets by maximum likelihood and Bayesian methods, which suggests the possible utility of CR sequences for phylogenetic reconstruction in certain groups of Hemiptera. Statistical analyses showed that TRs may contribute to the length variation, high AT content, and the shift of base skewing of CR sequences toward high AT content in the Hemiptera. Our findings enrich the knowledge of structural organization, phylogenetic utility, and roles of tandem replication of hemipteran CR, and provide a possible framework for mitochondrial control region analyses in hemimetabolous insects.

Characterization of the complete mitochondrial genome of Biston marginata (Lepidoptera: Geometridae) and phylogenetic analysis among lepidopteran insects

The complete mitochondrial genome (mitogenome) of Biston marginata (Lepidoptera: Geometridae) was determined and annotated. The circular genome is 15,470bp long and it contains the entire set of 37 genes usually present in lepidopteran mitogenomes. The nucleotide composition of the genome is highly A+T biased, accounting for 81.20%, with a slightly positive AT skewness (0.028), indicating the occurrence of more As than Ts, as found in other Geometridae species. Except for cox1 gene starts with non-canonical initial codon CGA, all protein-coding genes start with ATN codon. Three of the 13 PCGs (protein coding gene) had an incomplete termination codon, T or TA, while the others terminated with TAA. All tRNA genes are predicted to fold into typical clover-leaf secondary structure, except for the trnS1 (AGN), in which the dihydrouridine (DHU) arm could not form a stable stem-loop structure. The A+T-rich region of 343bp is comprised of non-repetitive sequences, but have several distinctive features, including the motif "ATAGA" followed by a 19bp poly-T stretch, a microsatellite-like (TA)₇ element next to the ATTTA motif. The phylogenetic analyses support the view that the B. marginata is closely related to the Biston pantrinaria, and confirm that Biston marginata belongs to the family Geometridae.

Revisiting Coleoptera a + T-rich region: structural conservation, phylogenetic and phylogeographic approaches in mitochondrial control region of bioluminescent Elateridae species (Coleoptera)

The control region (CR) or A + T-rich region in Coleoptera mt genome is poorly characterized, including the Elateroidea bioluminescent species. Here, we provided the first attempt to characterize and compare the structure and organization of the CR of different species within Elateridae. We also revisited some sequenced Coleoptera CR and observed consensus T-stretches, non-conserved sequences near the stem-loop and unusual inner tRNAs-like sequences. All these features are probably involved in the replication start of the mt genome. The phylogenetic relationships in Elateridae bioluminescent groups using partial sequence of CR showed the monophyly of Pyrearinus pumilus group and Pyrearinus as a polyphyletic genus, corroborating our previous results. The wider genetic variation obtained by CR analysis could separate two different lineages that occur within P. termitilluminans populations. In Elateridae, the CR exhibited high polymorphism within and between populations, which was also observed in other Coleoptera species, suggesting that the CR could be described as a suitable molecular marker to be applied in phylogenetic and phylogeographic studies.

The Phylogeny and Evolutionary Timescale of Muscoidea (Diptera: Brachycera: Calyptratae) Inferred from Mitochondrial Genomes

Muscoidea is a significant dipteran clade that includes house flies (Family Muscidae), latrine flies (F. Fannidae), dung flies (F. Scathophagidae) and root maggot flies (F. Anthomyiidae). It is comprised of approximately 7000 described species. The monophyly of the Muscoidea and the precise relationships of muscoids to the closest superfamily the Oestroidea (blow flies, flesh flies etc) are both unresolved. Until now mitochondrial (mt) genomes were available for only two of the four muscoid families precluding a thorough test of phylogenetic relationships using this data source. Here we present the first two mt genomes for the families Fanniidae (Euryomma sp.) (family Fanniidae) and Anthomyiidae (Delia platura (Meigen, 1826)). We also conducted phylogenetic analyses containing of these newly sequenced mt genomes plus 15 other species representative of dipteran diversity to address the internal relationship of Muscoidea and its systematic position. Both maximum-likelihood and Bayesian analyses suggested that Muscoidea was not a monophyletic group with the relationship: (Fanniidae + Muscidae) + ((Anthomyiidae + Scathophagidae) + (Calliphoridae + Sarcophagidae)), supported by the majority of analysed datasets. This also infers that Oestroidea was paraphyletic in the majority of analyses. Divergence time estimation suggested that the earliest split within the Calyptratae, separating (Tachinidae + Oestridae) from the remaining families, occurred in the Early Eocene. The main divergence within the paraphyletic muscoidea grade was between Fanniidae + Muscidae and the lineage ((Anthomyiidae + Scathophagidae) + (Calliphoridae + Sarcophagidae)) which occurred in the Late Eocene.

The first mitochondrial genome of the sepsid fly Nemopoda mamaevi Ozerov, 1997 (Diptera: Sciomyzoidea: Sepsidae), with mitochondrial genome phylogeny of cyclorrhapha

Sepsid flies (Diptera: Sepsidae) are important model insects for sexual selection research. In order to develop mitochondrial (mt) genome data for this significant group, we sequenced the first complete mt genome of the sepsid fly Nemopoda mamaevi Ozerov, 1997. The circular 15,878 bp mt genome is typical of Diptera, containing all 37 genes usually present in bilaterian animals. We discovered inaccurate annotations of fly mt genomes previously deposited on GenBank and thus re-annotated all published mt genomes of Cyclorrhapha. These re-annotations were based on comparative analysis of homologous genes, and provide a statistical analysis of start and stop codon positions. We further detected two 18 bp of conserved intergenic sequences from tRNAGlu-tRNAPhe and ND1-tRNASer(UCN) across Cyclorrhapha, which are the mtTERM binding site motifs. Additionally, we compared automated annotation software MITOS with hand annotation method. Phylogenetic trees based on the mt genome data from Cyclorrhapha were inferred by Maximum-likelihood and Bayesian methods, strongly supported a close relationship between Sepsidae and the Tephritoidea.

The mitochondrial genome of Dastarcus helophoroides (Coleoptera: Bothrideridae) and related phylogenetic analyses

The complete mitochondrial genome of Dastarcus helophoroides (Coleoptera: Bothrideridae) which consists of 13 PCGs, 22 tRNA genes, two rRNA genes and a non-coding region (D-loop), is sequenced for its nucleotide sequence of 15,878 bp (GenBank: KF811054.1). The genome has a typical gene order which is identical to other Coleoptera species. Except for COI gene generally starts with non-canonical initial codon, all protein-coding genes start with ATN codon and terminate with the stop codon TA(A) or TAG. The secondary structure of rrnL and rrnS consists of 48 helices (contains four newly proposed helices) and 35 helices (contains two newly proposed helices) respectively. All 22 tRNAs in D. helophoroides are predicted to fold into typical cloverleaf secondary structure, except trnS1 (AGN), in which the dihydrouracil arm (DHU arm) could not form stable stem-loop structure. Thirteen protein-coding genes (nucleotide dataset and nucleic acid dataset) of the available species (29 taxa) have been used to infer the phylogenetic relationships among these orders. Tenebrionoidea and Cucujoidea form a sister group, and D. helophoroides is classified into Cucujoidea (Bothrideridae). The study first research on the phylogenetic analyses involving to the D. helophoroides mitogenome, and the results strongly bolster the current morphology-based hypothesis.

The mitochondrial genome of Elodia flavipalpis Aldrich (Diptera: Tachinidae) and the evolutionary timescale of Tachinid flies

Tachinid flies are natural enemies of many lepidopteran and coleopteran pests of forests, crops, and fruit trees. In order to address the lack of genetic data in this economically important group, we sequenced the complete mitochondrial genome of the Palaearctic tachinid fly Elodia flavipalpis Aldrich, 1933. Usually found in Northern China and Japan, this species is one of the primary natural enemies of the leaf-roller moths (Tortricidae), which are major pests of various fruit trees. The 14,932-bp mitochondrial genome was typical of Diptera, with 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. However, its control region is only 105 bp in length, which is the shortest found so far in flies. In order to estimate dipteran evolutionary relationships, we conducted a phylogenetic analysis of 58 mitochondrial genomes from 23 families. Maximum-likelihood and Bayesian methods supported the monophyly of both Tachinidae and superfamily Oestroidea. Within the subsection Calyptratae, Muscidae was inferred as the sister group to Oestroidea. Within Oestroidea, Calliphoridae and Sarcophagidae formed a sister clade to Oestridae and Tachinidae. Using a Bayesian relaxed clock calibrated with fossil data, we estimated that Tachinidae originated in the middle Eocene.

The mitochondrial genome of the venomous cone snail Conus consors

Cone snails are venomous predatory marine neogastropods that belong to the species-rich superfamily of the Conoidea. So far, the mitochondrial genomes of two cone snail species (Conus textile and Conus borgesi) have been described, and these feed on snails and worms, respectively. Here, we report the mitochondrial genome sequence of the fish-hunting cone snail Conus consors and describe a novel putative control region (CR) which seems to be absent in the mitochondrial DNA (mtDNA) of other cone snail species. This possible CR spans about 700 base pairs (bp) and is located between the genes encoding the transfer RNA for phenylalanine (tRNA-Phe, trnF) and cytochrome c oxidase subunit III (cox3). The novel putative CR contains several sequence motifs that suggest a role in mitochondrial replication and transcription.

The complete mitochondrial genome of the citrus red mite Panonychus citri (Acari: Tetranychidae): high genome rearrangement and extremely truncated tRNAs

Background

The family Tetranychidae (Chelicerata: Acari) includes ~1200 species, many of which are of agronomic importance. To date, mitochondrial genomes of only two Tetranychidae species have been sequenced, and it has been found that these two mitochondrial genomes are characterized by many unusual features in genome organization and structure such as gene order and nucleotide frequency. The scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). Information on Tetranychidae mitochondrial genomes is quite important for phylogenetic evaluation and population genetics, as well as the molecular evolution of functional genes such as acaricide-resistance genes. In this study, we sequenced the complete mitochondrial genome of Panonychus citri (Family Tetranychidae), a worldwide citrus pest, and provide a comparison to other Acari.

Results

The mitochondrial genome of P. citri is a typical circular molecule of 13,077 bp, and contains the complete set of 37 genes that are usually found in metazoans. This is the smallest mitochondrial genome within all sequenced Acari and other Chelicerata, primarily due to the significant size reduction of protein coding genes (PCGs), a large rRNA gene, and the A + T-rich region. The mitochondrial gene order for P. citri is the same as those for P. ulmi and Tetranychus urticae, but distinctly different from other Acari by a series of gene translocations and/or inversions. The majority of the P. citri mitochondrial genome has a high A + T content (85.28%), which is also reflected by AT-rich codons being used more frequently, but exhibits a positive GC-skew (0.03). The Acari mitochondrial nad1 exhibits a faster amino acid substitution rate than other genes, and the variation of nucleotide substitution patterns of PCGs is significantly correlated with the G + C content. Most tRNA genes of P. citri are extremely truncated and atypical (44-65, 54.1 ± 4.1 bp), lacking either the T- or D-arm, as found in P. ulmi, T. urticae, and other Acariform mites.

Conclusions

The P. citri mitochondrial gene order is markedly different from those of other chelicerates, but is conserved within the family Tetranychidae indicating that high rearrangements have occurred after Tetranychidae diverged from other Acari. Comparative analyses suggest that the genome size, gene order, gene content, codon usage, and base composition are strongly variable among Acari mitochondrial genomes. While extremely small and unusual tRNA genes seem to be common for Acariform mites, further experimental evidence is needed.

Complete mtDNA genomes of Anopheles darlingi and an approach to anopheline divergence time

BACKGROUND

The complete sequences of the mitochondrial genomes (mtDNA) of members of the northern and southern genotypes of Anopheles (Nyssorhynchus) darlingi were used for comparative studies to estimate the time to the most recent common ancestor for modern anophelines, to evaluate differentiation within this taxon, and to seek evidence of incipient speciation.

METHODS

The mtDNAs were sequenced from mosquitoes from Belize and Brazil and comparative analyses of structure and base composition, among others, were performed. A maximum likelihood approach linked with phylogenetic information was employed to detect evidence of selection and a Bayesian approach was used to date the split between the subgenus Nyssorhynchus and other Anopheles subgenera.

RESULTS

The comparison of mtDNA sequences within the Anopheles darlingi taxon does not provide sufficient resolution to establish different units of speciation within the species. In addition, no evidence of positive selection in any protein-coding gene of the mtDNA was detected, and purifying selection likely is the basis for this lack of diversity. Bayesian analysis supports the conclusion that the most recent ancestor of Nyssorhynchus and Anopheles+Cellia was extant ~94 million years ago.

CONCLUSION

Analyses of mtDNA genomes of Anopheles darlingi do not provide support for speciation in the taxon. The dates estimated for divergence among the anopheline groups tested is in agreement with the geological split of western Gondwana (95 mya), and provides additional support for explaining the absence of Cellia in the New World, and Nyssorhynchus in the Afro-Eurasian continents.

Comparative mitochondrial genomics of freshwater mussels (Bivalvia: Unionoida) with doubly uniparental inheritance of mtDNA: gender-specific open reading frames and putative origins of replication

Doubly uniparental inheritance (DUI) of mitochondrial DNA in marine mussels (Mytiloida), freshwater mussels (Unionoida), and marine clams (Veneroida) is the only known exception to the general rule of strict maternal transmission of mtDNA in animals. DUI is characterized by the presence of gender-associated mitochondrial DNA lineages that are inherited through males (male-transmitted or M types) or females (female-transmitted or F types), respectively. This unusual system constitutes an excellent model for studying basic aspects of mitochondrial DNA inheritance and the evolution of mtDNA genomes in general. Here we compare published mitochondrial genomes of unionoid bivalve species with DUI, with an emphasis on characterizing unassigned regions, to identify regions of the F and M mtDNA genomes that could (i) play a role in replication or transcription of the mtDNA molecule and/or (ii) determine whether a genome will be transmitted via the female or the male gamete. Our results reveal the presence of one F-specific and one M-specific open reading frames (ORFs), and we hypothesize that they play a role in the transmission and/or gender-specific adaptive functions of the M and F mtDNA genomes in unionoid bivalves. Three major unassigned regions shared among all F and M unionoid genomes have also been identified, and our results indicate that (i) two of them are potential heavy-strand control regions (O(H)) for regulating replication and/or transcription and that (ii) multiple and potentially bidirectional light-strand origins of replication (O(L)) are present in unionoid F and M mitochondrial genomes. We propose that unassigned regions are the most promising candidate sequences in which to find regulatory and/or gender-specific sequences that could determine whether a mitochondrial genome will be maternally or paternally transmitted.

High divergence across the whole mitochondrial genome in the "pan-Antarctic" springtail Friesea grisea: evidence for cryptic species?

Collembola are one of the few hexapod groups adapted to live in the harsh environmental conditions of Antarctic terrestrial ecosystems. Diversity is limited to a few species that can be very abundant in coastal deglaciated sites. A remarkable lack of overlap in Collembola species composition is evident between Western and Eastern Antarctica, and Friesea grisea is currently the only species whose distribution is thought to span these two main regions of the continent. However, our analysis of the complete sequences of the mitochondrial genomes from specimens obtained from each of the two regions showed unexpected genetic divergence, well above the average levels observed between populations belonging to the same species, and so indicating that these are actually separate species, despite their lack of distinguishing morphology. Detailed analysis of the two genomes showed the presence of a non-coding region observed between trnS(uga) and nad1. Other features of these mitochondrial genomes, such as base compositional bias, secondary structure features of tRNAs and the presence of regulatory elements in the control region, are described and discussed from an evolutionary standpoint.

Complex patterns of genetic variability in populations of the New World screwworm fly revealed by mitochondrial DNA markers

Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae), the New World screwworm fly, is an important agent of traumatic myiasis, which is endemic in the Neotropical region and which has great economic impact on the development of the livestock industry. International efforts have been aimed at designing programmes to control and eradicate this species from endemic areas. Thorough knowledge of the population genetics of an insect pest is a fundamental component to ensuring the success of a pest management strategy because it enables the determination of an appropriate geographic scale for carrying out effective treatments. This study undertook an analysis of mtDNA polymerase chain reaction-restricted fragment length polymorphism (PCR-RFLP) in 34 populations of C. hominivorax from 10 countries, encompassing almost all the current distribution of the species. Results showed high levels of mitochondrial DNA variability (pi= 2.9%) and a complex pattern of population genetic structure for this species. Significant population structure (Phi st= 0.5234) and low variability were found in Caribbean populations, suggesting that, in general, islands constitute independent evolutionary entities connected by restricted gene flow. By contrast, high variability and low, but significant, differentiation was found among mainland populations (Phi st= 0.0483), which could not be attributed to geographic distance. Several processes may be acting to maintain the observed patterns, with different implications for establishing control programmes.

The complete mitochondrial genome of the Antarctic springtail Cryptopygus antarcticus (Hexapoda: Collembola)

BACKGROUND

Mitogenomics data, i.e. complete mitochondrial genome sequences, are popular molecular markers used for phylogenetic, phylogeographic and ecological studies in different animal lineages. Their comparative analysis has been used to shed light on the evolutionary history of given taxa and on the molecular processes that regulate the evolution of the mitochondrial genome. A considerable literature is available in the fields of invertebrate biochemical and ecophysiological adaptation to extreme environmental conditions, exemplified by those of the Antarctic. Nevertheless, limited molecular data are available from terrestrial Antarctic species, and this study represents the first attempt towards the description of a mitochondrial genome from one of the most widespread and common collembolan species of Antarctica.

RESULTS

In this study we describe the mitochondrial genome of the Antarctic collembolan Cryptopygus antarcticus Willem, 1901. The genome contains the standard set of 37 genes usually present in animal mtDNAs and a large non-coding fragment putatively corresponding to the region (A+T-rich) responsible for the control of replication and transcription. All genes are arranged in the gene order typical of Pancrustacea. Three additional short non-coding regions are present at gene junctions. Two of these are located in positions of abrupt shift of the coding polarity of genes oriented on opposite strands suggesting a role in the attenuation of the polycistronic mRNA transcription(s). In addition, remnants of an additional copy of trnL(uag) are present between trnS(uga) and nad1. Nucleotide composition is biased towards a high A% and T% (A+T = 70.9%), as typically found in hexapod mtDNAs. There is also a significant strand asymmetry, with the J-strand being more abundant in A and C. Within the A+T-rich region, some short sequence fragments appear to be similar (in position and primary sequence) to those involved in the origin of the N-strand replication of the Drosophila mtDNA.

CONCLUSION

The mitochondrial genome of C. antarcticus shares several features with other pancrustacean genomes, although the presence of unusual non-coding regions is also suggestive of molecular rearrangements that probably occurred before the differentiation of major collembolan families. Closer examination of gene boundaries also confirms previous observations on the presence of unusual start and stop codons, and suggests a role for tRNA secondary structures as potential cleavage signals involved in the maturation of the primary transcript. Sequences potentially involved in the regulation of replication/transcription are present both in the A+T-rich region and in other areas of the genome. Their position is similar to that observed in a limited number of insect species, suggesting unique replication/transcription mechanisms for basal and derived hexapod lineages. This initial description and characterization of the mitochondrial genome of C. antarcticus will constitute the essential foundation prerequisite for investigations of the evolutionary history of one of the most speciose collembolan genera present in Antarctica and other localities of the Southern Hemisphere.