Ancient DNA Analysis of Mid-Holocene Individuals from the Northwest Coast of North America Reveals Different Evolutionary Paths for Mitogenomes

To gain a better understanding of North American population history, complete mitochondrial genomes (mitogenomes) were generated from four ancient and three living individuals of the northern Northwest Coast of North America, specifically the north coast of British Columbia, Canada, current home to the indigenous Tsimshian, Haida, and Nisga’a. The mitogenomes of all individuals were previously unknown and assigned to new sub-haplogroup designations D4h3a7, A2ag and A2ah. The analysis of mitogenomes allows for more detailed analyses of presumed ancestor–descendant relationships than sequencing only the HVSI region of the mitochondrial genome, a more traditional approach in local population studies. The results of this study provide contrasting examples of the evolution of Native American mitogenomes. Those belonging to sub-haplogroups A2ag and A2ah exhibit temporal continuity in this region for 5000 years up until the present day. Of possible associative significance is that archaeologically identified house structures in this region maintain similar characteristics for this same period of time, demonstrating cultural continuity in residence patterns. The individual dated to 6000 years before present (BP) exhibited a mitogenome belonging to sub-haplogroup D4h3a. This sub-haplogroup was earlier identified in the same general area at 10300 years BP on Prince of Wales Island, Alaska, and may have gone extinct, as it has not been observed in any living individuals of the Northwest Coast. The presented case studies demonstrate the different evolutionary paths of mitogenomes over time on the Northwest Coast.     

Next-generation sequencing,phylogenetic signal and comparative mitogenomic analyses in Metacrangonyctidae (Amphipoda: Crustacea)

Background

Comparative mitochondrial genomic analyses are rare among crustaceans below the family or genus level. The obliged subterranean crustacean amphipods of the family Metacrangonyctidae, found from the Hispaniola (Antilles) to the Middle East, including the Canary Islands and the peri-Mediterranean region, have an evolutionary history and peculiar biogeography that can respond to Tethyan vicariance. Indeed, recent phylogenetic analysis using all protein-coding mitochondrial sequences and one nuclear ribosomal gene have lent support to this hypothesis (Bauzà-Ribot et al. 2012).

Results

We present the analyses of mitochondrial genome sequences of 21 metacrangonyctids in the genera Metacrangonyx and Longipodacrangonyx, covering the entire geographical range of the family. Most mitogenomes were attained by next-generation sequencing techniques using long-PCR fragments sequenced by Roche FLX/454 or GS Junior pyro-sequencing, obtaining a coverage depth per nucleotide of up to 281×. All mitogenomes were AT-rich and included the usual 37 genes of the metazoan mitochondrial genome, but showed a unique derived gene order not matched in any other amphipod mitogenome. We compare and discuss features such as strand bias, phylogenetic informativeness, non-synonymous/synonymous substitution rates and other mitogenomic characteristics, including ribosomal and transfer RNAs annotation and structure.

Conclusions

Next-generation sequencing of pooled long-PCR amplicons can help to rapidly generate mitogenomic information of a high number of related species to be used in phylogenetic and genomic evolutionary studies. The mitogenomes of the Metacrangonyctidae have the usual characteristics of the metazoan mitogenomes (circular molecules of 15,000-16,000 bp, coding for 13 protein genes, 22 tRNAs and two ribosomal genes) and show a conserved gene order with several rearrangements with respect to the presumed Pancrustacean ground pattern. Strand nucleotide bias appears to be reversed with respect to the condition displayed in the majority of crustacean mitogenomes since metacrangonyctids show a GC-skew at the (+) and (-) strands; this feature has been reported also in the few mitogenomes of Isopoda (Peracarida) known thus far. The features of the rRNAs, tRNAs and sequence motifs of the control region of the Metacrangonyctidae are similar to those of the few crustaceans studied at present.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-566) contains supplementary material, which is available to authorized users.     

Characterization of the complete mitochondrial genomes of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae)

The complete mitochondrial genomes (mitogenomes) of Cnaphalocrocis medinalis and Chilo suppressalis (Lepidoptera: Pyralidae) were determined and analyzed. The circular genomes were 15,388 bp long for C. medinalis and 15,395 bp long for C. suppressalis. Both mitogenomes contained 37 genes, with gene order similar to that of other lepidopterans. Notably, 12 protein-coding genes (PCGs) utilized the standard ATN, but the cox1 gene used CGA as the initiation codon; the cox1, cox2, and nad4 genes in the two mitogenomes had the truncated termination codons T, T, and TA, respectively, but the nad5 gene was found to use T as the termination codon only in the C. medinalis mitogenome. Additionally, the codon distribution and Relative Synonymous Codon Usage of the 13 PCGs in the C. medinalis mitogenome were very different from those in other pyralid moth mitogenomes. Most of the tRNA genes had typical cloverleaf secondary structures. However, the dihydrouridine (DHU) arm of the trnS1(AGN) gene did not form a stable stem-loop structure. Forty-nine helices in six domains, and 33 helices in three domains were present in the secondary structures of the rrnL and rrnS genes of the two mitogenomes, respectively. There were four major intergenic spacers, except for the A+T-rich region, spanning at least 12 bp in the two mitogenomes. The A+T-rich region contained an 'ATAGT(A)'-like motif followed by a poly-T stretch in the two mitogenomes. In addition, there were a potential stem-loop structure, a duplicated 25-bp repeat element, and a microsatellite '(TA)(13)' observed in the A+T-rich region of the C. medinalis mitogenome. A poly-T motif, a duplicated 31-bp repeat element, and a 19-bp triplication were found in the C. suppressalis mitogenome. However, there are many differences in the A+T-rich regions between the C. suppressalis mitogenome sequence in the present study and previous reports. Finally, the phylogenetic relationships of these insects were reconstructed based on amino acid sequences of mitochondrial 13 PCGs using Bayesian inference and maximum likelihood methods. These molecular-based phylogenies support the traditional morphologically based view of relationships within the Pyralidae.     

Next‐generation sequencing workflow for assembly of nonmodel mitogenomes exemplified with North Pacific albatrosses (Phoebastria spp.)

Use of complete mitochondrial genomes (mitogenomes) can greatly increase the resolution achievable in phylogeographic and historical demographic studies. Using next‐generation sequencing methods, it is now feasible to efficiently sequence mitogenomes of large numbers of individuals once a reference mitogenome is available. However, assembling the initial mitogenomes of nonmodel organisms can present challenges, for example, in birds, where mtDNA is often subject to gene rearrangements and duplications. We developed a workflow based on Illumina paired‐end, whole‐genome shotgun sequencing, which we used to generate complete 19‐kilobase mitogenomes for each of three species of North Pacific albatross, a group of birds known to carry a tandem duplication. Although this duplication had been described previously, our procedure did not depend on this prior knowledge, nor did it require a closely related reference mitogenome (e.g. a mammalian mitogenome was sufficient). We employed an iterative process including de novo assembly, reference‐guided assembly and gap closing, which enabled us to detect duplications, determine gene order and identify sequence for primer positioning to resolve any mitogenome ambiguity (via minimal targeted Sanger sequencing). We present full mtDNA annotations, including 22 tRNAs, 2 rRNAs, 13 protein‐coding genes, a control region and a duplicated feature for all three species. Pairwise comparisons supported previous hypotheses regarding the phylogenetic relationships within this group and occurrence of a shared tandem duplication. The resulting mitogenome sequences will enable rapid, high‐throughput NGS mitogenome sequencing of North Pacific albatrosses via direct reference‐guided assembly. Moreover, our approach to assembling mitogenomes should be applicable to any taxon.     

Phylogeography of mitochondrial haplogroup D1: An early spread of subhaplogroup D1j from Central Argentina

We analyzed the patterns of variation of haplogroup D1 in central Argentina, including new data and published information from other populations of South America. Almost 28% (107/388) of the individuals sampled in the region belong to haplogroup D1, whereas more than 52% of them correspond to the recently described subhaplogroup D1j (Bodner et al.: Genome Res 22 (2012) 811–820), defined by the presence of additional transitions at np T152C–C16242T–T16311C to the nodal D1 motif. This lineage was found at high frequencies across a wide territory with marked geographical–ecological differences. Additionally, 12 individuals present the mutation C16187T that defines the recently named subhaplogroup D1g (Bodner et al.: Genome Res 22 (2012) 811–820), previously described in populations of Patagonia and Tierra del Fuego. Based on our results and additional data already published, we postulate that the most likely origin of subhaplogroup D1j is the region of Sierras Pampeanas, which occupies the center and part of the northwestern portion of Argentina. The extensive yet restricted geographical distribution, the relatively large internal diversity, and the absence or low incidence of D1j in other regions of South America suggest the existence of an ancient metapopulation covering the Sierras Pampeanas, being this lineage its genetic signature. Further support for a scenario of local origin for D1j in the Sierras Pampeanas stems from the fact that early derivatives from a putative ancestral lineage carrying the transitions T16311C–T152C have only been found in this region, supporting the hypothesis that it might represent an ancestral motif previous to the appearance of D1j‐specific change C16242T. © 2012 Wiley Periodicals, Inc.     

The complete mitogenome assemblies of 10 diploid potato clones reveal recombination and overlapping variants

The potato mitogenome is complex and to understand various biological functions and nuclear–cytoplasmic interactions, it is important to characterize its gene content and structure. In this study, the complete mitogenome sequences of nine diploid potato clones along with a diploid Solanum okadae clone were characterized. Each mitogenome was assembled and annotated from Pacific Biosciences (PacBio) long reads and 10X genomics short reads. The results show that each mitogenome consists of multiple circular molecules with similar structure and gene organization, though two groups (clones 07506-01, DW84-1457, 08675-21 and H412-1 in one group, and clones W5281-2, 12625-02, 12120-03 and 11379-03 in another group) could be distinguished, and two mitogenomes (clone 10908-06 and OKA15) were not consistent with those or with each other. Significant differences in the repeat structure of the 10 mitogenomes were found, as was recombination events leading to multiple sub-genomic circles. Comparison between individual molecules revealed a translocation of ∼774 bp region located between a short repeat of 40 bp in molecule 3 of each mitogenome, and an insertion of the same in molecule 2 of the 10908-06 mitogenome. Finally, phylogenetic analyses revealed a close relationship between the mitogenomes of these clones and previously published potato mitogenomes.     

Prospects on the evolutionary mitogenomics of plants: A case study on the olive family (Oleaceae)

The mitogenome is rarely used to reconstruct the evolutionary history of plants, contrary to nuclear and plastid markers. Here, we evaluate the usefulness of mitochondrial DNA for molecular evolutionary studies in Oleaceae, in which cases of cytoplasmic male sterility (CMS) and of potentially contrasted organelle inheritance are known. We compare the diversity and the evolution of mitochondrial and chloroplast genomes by focusing on the olive complex and related genera. Using high‐throughput techniques, we reconstructed complete mitogenomes (ca. 0.7 Mb) and plastomes (ca. 156 kb) for six olive accessions and one Chionanthus. A highly variable organization of mitogenomes was observed at the species level. In olive, two specific chimeric genes were identified in the mitogenome of lineage E3 and may be involved in CMS. Plastid‐derived regions (mtpt) were observed in all reconstructed mitogenomes. Through phylogenetic reconstruction, we demonstrate that multiple integrations of mtpt regions have occurred in Oleaceae, but mtpt regions shared by all members of the olive complex derive from a common ancestor. We then assembled 52 conserved mitochondrial gene regions and complete plastomes of ten additional accessions belonging to tribes Oleeae, Fontanesieae and Forsythieae. Phylogenetic congruence between topologies based on mitochondrial regions and plastomes suggests a strong disequilibrium linkage between both organellar genomes. Finally, while phylogenetic reconstruction based on plastomes fails to resolve the evolutionary history of maternal olive lineages in the Mediterranean area, their phylogenetic relationships were successfully resolved with complete mitogenomes. Overall, our study demonstrates the great potential of using mitochondrial DNA in plant phylogeographic and metagenomic studies.     

Tandem Duplication and Random Loss for mitogenome rearrangement in Symphurus (Teleost: Pleuronectiformes)

Background

The mitochondrial genomes (mitogenomes) of flatfishes (Pleuronectiformes) exhibit highly diversified types of large-scale gene rearrangements. We have reported that the mitogenomes of Crossorhombus azureus (Bothidae), Samariscus latus (Samaridae) and Cynoglossus fishes (Cynoglossidae) show different types of gene rearrangements.

Results

In the present study, the complete mitogenomes of two Symphurus species (Cynoglossidae), Symphurus plagiusa and Symphurus orientalis, were determined. The gene order in the S. plagiusa mitogenome is the same as that of a typical vertebrate (without any gene rearrangements). Surprisingly, large-scale gene rearrangements have occurred in S. orientalis. In the rearranged fragment from the control region (CR) to the WANCY tRNA cluster (tRNA cluster of tRNA-W, tRNA-A, tRNA-N, tRNA-C and tRNA-Y) in the S. orientalis mitogenome, tRNA-V and tRNA-M have been translocated to the 3’ end of the 16S rRNA gene, with six large intergenic spacers over 20 bp in length. In addition, an origin for light-strand replication (O_L) structure that is typically located in the WANCY region was absent in both the S. plagiusa and S. orientalis mitogenomes. It is generally recognized that a sequence in the WANCY region that encodes tRNAs forms a hairpin structure (O_L-like structure) and can act as the O_L when the typical locus is lost. Moreover, an additional O_L-like structure was identified near the control region in the S. plagiusa mitogenome.

Conclusions

The positions of the intergenic spacers and the rearranged genes of the S. orientalis mitogenome strongly indicate that the mechanism underlying the rearrangement of this mitogenome was Tandem Duplication and Random Loss. Additionally, two O_L-like regions substituting for the typical locus were found in the S. plagiusa mitogenome. We speculate that the ancestral mitogenomes of S. plagiusa and S. orientalis also had this characteristic, such that if both O_L-like structures functioned during mitochondrial replication, they could initiate duplicate replications of the light strand (L-strand), leading to duplication of the region between the two structures. We consider that this mechanism may account for the gene duplication that occurred during the gene rearrangement process in the evolution of the ancestral mitogenome to the S. orientalis mitogenome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1581-6) contains supplementary material, which is available to authorized users.     

Differentiation of the Mitochondrial Subhaplogroup U4 in the Populations of Eastern Europe,Ural, and Western Siberia: Implication to the Genetic History of the Uralic Populations

Phylogenetic relationships between the sequences of the mitochondrial DNA (mtDNA) hypervariable segment 1, belonging to subhaplogroup U4, were examined in the populations of Eastern Europe, Ural, and Northwest Siberia. It was shown that the frequency of subhaplogroup U4, as well as its proportion in the U-component of the gene pools, increased eastwards, reaching maximum values in the populations of Northwest Siberia. Phylogenetic analysis it was showed that the appearance of specific U4-lineage (16113C–16356–16362) in the ancestors of Mansi was most likely caused by its divergence from the East European cluster 16356–16362 in the Late Upper Paleolithic (18566 ± 12915 years before present). Other U4 mtDNA lineages (16189–16356 and 16311–16356), typical mostly of the indigenous populations of Northwest Siberia (Mansi, Nganasans, and Kets) may have formed during the Neolithic–early Bronze Age (6055 ± 3599 years before present, on average). It seems likely that the isolation of ancient populations inhabiting the region between the Ob' and Yenisei rivers was the key factor, providing the appearance of the unique Caucasoid mtDNA lineages in their gene pools. These results were consistent with the traditional point of view on the mixed origin of the Finno-Ugric populations of the Volga–Ural region and West Siberia, resulted from the genetic relationships between the populations of Europe and Asia.     

Mitochondrial DNA from El Mirador Cave (Atapuerca,Spain) Reveals the Heterogeneity of Chalcolithic Populations

Previous mitochondrial DNA analyses on ancient European remains have suggested that the current distribution of haplogroup H was modeled by the expansion of the Bell Beaker culture (ca 4,500–4,050 years BP) out of Iberia during the Chalcolithic period. However, little is known on the genetic composition of contemporaneous Iberian populations that do not carry the archaeological tool kit defining this culture. Here we have retrieved mitochondrial DNA (mtDNA) sequences from 19 individuals from a Chalcolithic sample from El Mirador cave in Spain, dated to 4,760–4,200 years BP and we have analyzed the haplogroup composition in the context of modern and ancient populations. Regarding extant African, Asian and European populations, El Mirador shows affinities with Near Eastern groups. In different analyses with other ancient samples, El Mirador clusters with Middle and Late Neolithic populations from Germany, belonging to the Rössen, the Salzmünde and the Baalberge archaeological cultures but not with contemporaneous Bell Beakers. Our analyses support the existence of a common genetic signal between Western and Central Europe during the Middle and Late Neolithic and points to a heterogeneous genetic landscape among Chalcolithic groups.     

The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

ABSTRACT: BACKGROUND: Lepidoptera encompasses more than 160,000 described species that have been classified into 45-48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the most derived lepidopteran lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to primitive lepidopteran lineages and conducted a comparative mitogenomic analysis across Lepidoptera. RESULTS: The mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,814 bp long with an A+T content of 81.28% and 82.33%, respectively. Different tandem repeats in the A+T-rich region mainly account for the size difference between the two mitogenomes. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. The 1,584-bp sequence from rrnS to nad2 was also determined for Thitarodes sp.QL, which has no repetitive sequence in the A+T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A+T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif 'ATAGA' and poly-T stretch in the A+T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes. The phylogenetic analysis supports that Hepialoidea, represented by T. renzhiensis and T. yunnanensis, occupies a basal position in the currently sampled seven superfamilies. The relationships of the other six superfamilies are (((((Bombycoidea + Geometroidea) + Noctuoidea) + Pyraloidea) + Papilionoidea) + Tortricoidea). CONCLUSION: The mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event occurred after Lepidoptera diverged from other holometabolous insect orders. Phylogenetic analysis based on mitogenome sequences is a power tool for addressing phylogenetic relationships among major Lepidoptera superfamilies. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.     

The complete mitochondrial genome sequence of the world's largest fish,the whale shark (Rhincodon typus), and its comparison with those of related shark species

The whale shark (Rhincodon typus) is the largest extant species of fish, belonging to the order Orectolobiformes. It is listed as a “vulnerable” species on the International Union for Conservation of Nature (IUCN)'s Red List of Threatened Species, which makes it an important species for conservation efforts. We report here the first complete sequence of the mitochondrial genome (mitogenome) of the whale shark obtained by next-generation sequencing methods. The assembled mitogenome is a 16,875 bp circle, comprising of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region. We also performed comparative analysis of the whale shark mitogenome to the available mitogenome sequences of 17 other shark species, four from the order Orectolobiformes, five from Lamniformes and eight from Carcharhiniformes. The nucleotide composition, number and arrangement of the genes in whale shark mitogenome are the same as found in the mitogenomes of the other members of the order Orectolobiformes and its closest orders Lamniformes and Carcharhiniformes, although the whale shark mitogenome had a slightly longer control region. The availability of mitogenome sequence of whale shark will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this species.     

The complete mitochondrial genome of Mahanta tanyae compared with other zygaenoid moths (Lepidoptera: Zygaenoidea)

The complete mitochondrial genome (mitogenome) of Mahanta tanyae was sequenced and extensively compared with all seven additionally reported zygaenoid mitogenomes. The M. tanyae mitogenome is circular, double-stranded, and 15,323 bp long. Gene content, gene order, and orientation are all typical of Lepidoptera, despite the existence of gene rearrangements for some other zygaenoid mitogenomes. Comparative analyses further showed that the incomplete termination codon T is consistently recognized in the mitochondrial cox1, cox2 and nad4 genes of all zygaenoid species, as well as in the nad5 gene in two limacodid species. Among 13 protein-coding genes, nad6 exhibits the highest evolutionary rate. The structure for each tRNA is highly conserved, including loss of the dihydorouidine (DHU) arm in trnS1 (AGN), but remarkable nucleotide variation exists, primarily in the pseudouridine (TψC) loops. Interestingly, in four species of Zygaenidae, the anticodons for trnS1 (AGN) are consistently UCU, instead of the routinely used codon GCU, in all three species of Limacodidae. In the intergenic region between trnS2 and nad1, a short sequence before the motif “ATACTAA” is present in the M. tanyae mitogenome that is unique among reported zygaenoid mitogenomes. In the A + T-rich region between the motif “ATTTA” and the microsatellite (AT)_n element, some nucleotides were present for most zygaenoid mitogenomes, which is, to our knowledge, rare even in reported lepidopteran mitogenomes. Phylogenetic analyses based on the combined 37 mitochondrial genes confirmed the position of M. tanyae in Limacodidae of the Zygaenoidea.     

Gene rearrangements in gekkonid mitochondrial genomes with shuffling,loss, and reassignment of tRNA genes

Background

Vertebrate mitochondrial genomes (mitogenomes) are 16–18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

Results

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNA^Arg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNA^Arg gene appears to have been converted to a tRNA^Gln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNA^Leu genes, as well as a translocation of the tRNA^Ala gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNA^Glu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNA^Glu gene adjacent to the 5’ end of the major noncoding region.

Conclusions

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-930) contains supplementary material, which is available to authorized users.     

Positioning the Red Deer (Cervus elaphus) Hunted by the Tyrolean Iceman into a Mitochondrial DNA Phylogeny

In the last years several phylogeographic studies of both extant and extinct red deer populations have been conducted. Three distinct mitochondrial lineages (western, eastern and North-African/Sardinian) have been identified reflecting different glacial refugia and postglacial recolonisation processes. However, little is known about the genetics of the Alpine populations and no mitochondrial DNA sequences from Alpine archaeological specimens are available. Here we provide the first mitochondrial sequences of an Alpine Copper Age Cervus elaphus. DNA was extracted from hair shafts which were part of the remains of the clothes of the glacier mummy known as the Tyrolean Iceman or Ötzi (5,350–5,100 years before present). A 2,297 base pairs long fragment was sequenced using a mixed sequencing procedure based on PCR amplifications and 454 sequencing of pooled amplification products. We analyzed the phylogenetic relationships of the Alpine Copper Age red deer''s haplotype with haplotypes of modern and ancient European red deer. The phylogenetic analyses showed that the haplotype of the Alpine Copper Age red deer falls within the western European mitochondrial lineage in contrast with the current populations from the Italian Alps belonging to the eastern lineage. We also discussed the phylogenetic relationships of the Alpine Copper Age red deer with the populations from Mesola Wood (northern Italy) and Sardinia.     

The mitogenomes of three beetles (Coleoptera: Polyphaga: Cucujiformia): New gene rearrangement and phylogeny

The mitochondrial genome (mitogenome) has been extensively used for studying phylogenetic relationships at different taxonomic levels. Several molecular analyses have been performed, but the phylogenetic relationships among infraorders in Polyphaga have not been well resolved. In this work, three nearly complete mitogenomes of Coleoptera, Sitophilus oryzae, Oryzaephilus surinamensis and Callosobruchus chinensis, were determined. The O. surinamensis and S. oryzae mitogenomes harbor gene content typical of other Polyphaga mitogenomes, while a gene rearrangement (trnQ) was found in the C. chinensis mitogenome. The mitogenomes of these three Coleoptera species each consist of approximately 13 protein-coding genes, 22 tRNA genes, two rRNA genes and one A + T-rich region. Phylogenetic analysis within Polyphaga was carried out based on mitochondrial data. The phylogenetic results within Polyphaga support the basal position of Cyphon sp., which belonged to Scirtoidea, Elateriformia. Within Cucujiformia, monophyletic Curculionoidea, Chrysomeloidea and Tenebrionoidea were confirmed.     

Mitogenomics of the suborder Cottoidei (Teleostei: Perciformes): Improved assemblies,mitogenome features,phylogeny, and ecological implications

We determined the mitogenome of Cyclopterus lumpus using a hybrid sequencing approach, and another four closely related species in the Liparidae based on available next-generation sequence data. We found that the mitogenome of C. lumpus was 17,266 bp in length, where the length and organisation were comparable to those reported for cottoids. However, we found a GC-homopolymer region in the intergenic space between tRNA^Leu2 and ND1 in liparids and cyclopterids. Phylogenetic reconstruction confirmed the monophyly of infraorders and firmly supported a sister-group relationship between Cyclopteridae and Liparidae. Purifying selection was the predominant force in the evolution of cottoid mitogenomes. There was significant evidence of relaxed selective pressures along the lineage of deep-sea fish, while selection was intensified in the freshwater lineage. Overall, our analysis provides a necessary expansion in the availability of mitogenomic sequences and sheds light on mitogenomic adaptation in Cottoidei fish inhabiting different aquatic environments.     

Mitochondrial genomes of five Hyphessobrycon tetras and their phylogenetic implications

To date, the taxonomic status and phylogenetic affinities within Hyphessobrycon, even among other genera in Characidae, remain unclear. Here, we determined five new mitochondrial genomes (mitogenomes) of Hyphessobrycon species (H. elachys, H. flammeus, H. pulchripinnis, H. roseus, and H. sweglesi). The mitogenomes were all classical circular structures, with lengths ranging from 16,008 to 17,224 bp. The type of constitutive genes and direction of the coding strand that appeared in the mitogenomes were identical to those of other species in Characidae. The highest value of the Ka/Ks ratio within 13 protein‐coding genes (PCGs) was found in ND2 with 0.83, suggesting that they were subject to purifying selection in the Hyphessobrycon genus. Comparison of the control region sequences among seven Hyphessobrycon fish revealed that repeat units differ in length and copy number across different species, which led to sharp differences in mitogenome sizes. Phylogenetic trees based on the 13 PCGs did not support taxonomic relationships, as the Hyphessobrycon fish mixed with those from other genera. These data were combined to explore higher level relationships within Characidae and could aid in the understanding of the evolution of this group.