Characterization of the complete mitochondrial genome of Corvus corone orientalis

In this study, the complete mitochondrial genome of Corvus corone orientalis was assembled through next-generation sequencing data. This circular mitochondrial genome of C. corone orientalis is 16,947 bp in length and has a base composition of A (30.8%), T (24.7%), C (29.9%), and G (14.5%), demonstrating a bias of higher AT content (55.5%) than GC content (44.5%). The mitochondrial genome contains a typically conserved structure among bird mitogenomes, encoding 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), two ribosomal RNA genes (12S rRNA and 16S rRNA), and a control region (D-loop region). Except ND6, all other PCGs were located on the H-strand. ATP8 gene and ATP6 gene were overlapped by 8 bp. The whole mt genome of C. corone orientalis and other Corvoidea mitogenomes (24 species, in total) were used for phylogenetic analysis. The result indicated C. corone orientalis has the closest relationship with Corvus cornix cornix (NC_024698) and clustered within clade of genus Corvus.

The complete mitochondrial genome of a leaf roller, Eudemis lucina (Lepidoptera: Tortricidae)

The leaf roller, Eudemis lucina, is a potential pest of Quercus in East Asia. In this study, we described the complete mitochondrial genome of this species by high-throughput sequencing. The mitochondrial genome is found to be a circular molecule of 16,056 bp in length, which consisted of 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and a non-coding control region (A + T-rich region). The A + T content is 80.5% for the whole mitogenome. All PCGs are initiated by ATN codons, except for COI which is initiated by the CGA codon. Eight PCGs use a typical stop codon of TAA, whereas the remaining PCGs use incomplete stop codon of T–– or TA–. The non-coding control region is 1013 bp and located between s-rRNA and Met-tRNA.

Mitochondrial genome of a brachypterous species in Meconematinae: Acosmetura nigrogeniculata and its phylogenetic implication

Acosmetura nigrogeniculata (Liu and Wang, 1998) is a brachypterous species in Meconematinae, which is only distributed in China. In this study, the complete mitochondrial genome of A. nigrogeniculata was determined and annotated. The 16,271 bp circular genome contained 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region. The overall base composition was 36.4% A, 34.8% T, 18.4% C, and 10.5% G, exhibiting obvious anti-G and AT bias (71.2%). The general genomic characters including nucleotides composition, gene arrangement, and codon usage were similar to those of other Meconematinae species. Phylogenetic analysis of all nine Meconematinae species indicated that the newly sequenced species were clustered closely with the brachypterous species Pseudosmetura snjiensis.

The complete mitochondrial genome of Aquila nipalensis and its phylogenetic position

Mitochondrial genome sequences are valuable resources for systematics and conservation biology studies. In this paper, we present the complete mitogenome of Aquila nipalensis which was 18,450 bp in length. The gene content and arrangement were typical for avian mtDNA. The overall A + T content of was 54.1%, and the AT skew was calculated as 0.12 for the complete mitogenome of A. nipalensis. The maximum-likelihood (ML) tree based on the concatenated 12 protein-coding genes (PCGs) revealed the basal phylogenetic position of A. nipalensis in Aquila.

The complete mitochondrial genome of Rhynchocypris oxycephalus (Teleostei: Cyprinidae) and its phylogenetic implications

Rhynchocypris oxycephalus (Teleostei: Cyprinidae) is a typical small cold water fish, which is distributed widely and mainly inhabits in East Asia. Here, we sequenced and determined the complete mitochondrial genome of R. oxycephalus and studied its phylogenetic implication. R. oxycephalus mitogenome is 16,609 bp in length (GenBank accession no.: MH885043), and it contains 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and two noncoding regions (the control region and the putative origin of light-strand replication). 12 PCGs started with ATG, while COI used GTG as the start codon. The secondary structure of tRNA-Ser (AGN) lacks the dihydrouracil (DHU) arm. The control region is 943bp in length, with a termination-associated sequence, six conserved sequence blocks (CSB-1, CSB-2, CSB-3, CSB-D, CSB-E, CSB-F), and a repetitive sequence. Phylogenetic analysis was performed with maximum likelihood and Bayesian methods based on the concatenated nucleotide sequence of 13 PCGs and the complete sequence without control region, and the result revealed that the relationship between R. oxycephalus and R. percnurus is closest, while the relationship with R. kumgangensis is farthest. The genus Rhynchocypris is revealed as a polyphyletic group, and R. kumgangensis had distant relationship with other Rhynchocypris species. In addition, COI and ND2 genes are considered as the fittest DNA barcoding gene in genus Rhynchocypris. This work provides additional molecular information for studying R. oxycephalus conservation genetics and evolutionary relationships.

Population substructure and signals of divergent adaptive selection despite admixture in the sponge Dendrilla antarctica from shallow waters surrounding the Antarctic Peninsula

Antarctic shallow-water invertebrates are exceptional candidates to study population genetics and evolution, because of their peculiar evolutionary history and adaptation to extreme habitats that expand and retreat with the ice sheets. Among them, sponges are one of the major components, yet population connectivity of none of their many Antarctic species has been studied. To investigate gene flow, local adaptation and resilience to near-future changes caused by global warming, we sequenced 62 individuals of the sponge Dendrilla antarctica along the Western Antarctic Peninsula (WAP) and the South Shetlands (spanning ~900 km). We obtained information from 577 double digest restriction site-associated DNA sequencing (ddRADseq)-derived single nucleotide polymorphism (SNP), using RADseq techniques for the first time with shallow-water sponges. In contrast to other studies in sponges, our 389 neutral SNPs data set showed high levels of gene flow, with a subtle substructure driven by the circulation system of the studied area. However, the 140 outlier SNPs under positive selection showed signals of population differentiation, separating the central-southern WAP from the Bransfield Strait area, indicating a divergent selection process in the study area despite panmixia. Fourteen of these outliers were annotated, being mostly involved in immune and stress responses. We suggest that the main selective pressure on D. antarctica might be the difference in the planktonic communities present in the central-southern WAP compared to the Bransfield Strait area, ultimately depending on sea-ice control of phytoplankton blooms. Our study unveils an unexpectedly long-distance larval dispersal exceptional in Porifera, broadening the use of genome-wide markers within nonmodel Antarctic organisms.

Potential of Transcript Editing Across Mitogenomes of Early Land Plants Shows Novel and Familiar Trends

RNA editing alters the identity of nucleotides in an RNA sequence so that the mature transcript differs from the template defined in the genome. This process has been observed in chloroplasts and mitochondria of both seed and early land plants. However, the frequency of RNA editing in plant mitochondria ranges from zero to thousands of editing sites. To date, analyses of RNA editing in mitochondria of early land plants have been conducted on a small number of genes or mitochondrial genomes of a single species. This study provides an overview of the mitogenomic RNA editing potential of the main lineages of these two groups of early land plants by predicting the RNA editing sites of 33 mitochondrial genes of 37 species of liverworts and mosses. For the purpose of the research, we newly assembled seven mitochondrial genomes of liverworts. The total number of liverwort genera with known complete mitogenome sequences has doubled and, as a result, the available complete mitogenome sequences now span almost all orders of liverworts. The RNA editing site predictions revealed that C-to-U RNA editing in liverworts and mosses is group-specific. This is especially evident in the case of liverwort lineages. The average level of C-to-U RNA editing appears to be over three times higher in liverworts than in mosses, while the C-to-U editing frequency of the majority of genes seems to be consistent for each gene across bryophytes.

Characterization of the Complete Mitochondrial Genome of a Whipworm Trichuris skrjabini (Nematoda: Trichuridae)

The complete mitochondrial (mt) genome of Trichuris skrjabini has been determined in the current study and subsequently compared with closely related species by phylogenetic analysis based on concatenated datasets of mt amino acid sequences. The whole mt genome of T. skrjabini is circular and 14,011 bp in length. It consists of a total of 37 genes including 13 protein coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNAs) genes, and two non-coding regions. The gene arrangement and contents were consistent with other members of the Trichuridae family including Trichuris suis, Trichuris trichiura, Trichuris ovis, and Trichuris discolor. Phylogenetic analysis based on concatenated datasets of amino acids of the 12 PCGs predicted the distinctiveness of Trichuris skrjabini as compared to other members of the Trichuridae family. Overall, our study supports the hypothesis that T. skrjabini is a distinct species. The provision of molecular data of whole mt genome of T. skrjabini delivers novel genetic markers for future studies of diagnostics, systematics, population genetics, and molecular epidemiology of T. skrjabini.

Biparental inheritance of mitochondrial DNA in humans is not a common phenomenon

PURPOSE

A recent report has raised the possibility of biparental mitochondrial DNA (mtDNA) inheritance, which could lead to concerns by health-care professionals and patients regarding investigations and genetic counseling of families with pathogenic mitochondrial DNA variants. Our aim was to examine the frequency of this phenomenon by investigating a cohort of patients with suspected mitochondrial disease.

METHODS

We studied genome sequencing (GS) data of DNA extracted from blood samples of 41 pediatric patients with suspected mitochondrial disease and their parents.

RESULTS

All of the mtDNA variants in the probands segregated with their mother or were apparently de novo. There were no variants that segregated only with the father and none of these families showed evidence of biparental inheritance of their mtDNA.

CONCLUSION

Paternal mitochondrial transmission is unlikely to be a common occurrence and therefore at this point we would not recommend changes in clinical practice.

Signals of positive selection in mitochondrial protein-coding genes of woolly mammoth: Adaptation to extreme environments?

The mammoths originated in warm and equatorial Africa and later colonized cold and high-latitude environments. Studies on nuclear genes suggest that woolly mammoth had evolved genetic variations involved in processes relevant to cold tolerance, including lipid metabolism and thermogenesis, and adaptation to extremely varied light and darkness cycles. The mitochondria is a major regulator of cellular energy metabolism, thus the mitogenome of mammoths may also exhibit adaptive evolution. However, little is yet known in this regard. In this study, we analyzed mitochondrial protein-coding genes (MPCGs) sequences of 75 broadly distributed woolly mammoths (Mammuthus primigenius) to test for signatures of positive selection. Results showed that a total of eleven amino acid sites in six genes, namely ND1, ND4, ND5, ND6, CYTB, and ATP6, displayed strong evidence of positive selection. Two sites were located in close proximity to proton-translocation channels in mitochondrial complex I. Biochemical and homology protein structure modeling analyses demonstrated that five amino acid substitutions in ND1, ND5, and ND6 might have influenced the performance of protein-protein interaction among subunits of complex I, and three substitutions in CYTB and ATP6 might have influenced the performance of metabolic regulatory chain. These findings suggest metabolic adaptations in the mitogenome of woolly mammoths in relation to extreme environments and provide a basis for further tests on the significance of the variations on other systems.

Diaphorina citri (Hemiptera: Psylloidea) in China: Two Invasion Routes and Three Transmission Paths

Diaphorina citri Kuwayama (Hemiptera: Liviidae) is one of the most common pests impacting citrus orchards in southern China. Samples of D. citri were collected in southern China in order to systematically explore the genetic architecture of the species. Mitochondrial cytochrome b (Cytb) and cytochrome coxidase subunit I (COI) were amplified by polymerase chain reaction (PCR) which allowed highlighting low haplotype and nucleotide diversities among the population. Two clades could be observed in the haplotype network. Moreover, Bayesian and maximum parsimony phylogenetic trees were constructed based on the sequences of Cytb and COI. Here, we report on the significant genetic variation of the species when comparing southwestern China with other regions of southern China (southern and southeastern). This analysis also suggested that the genetic structure of D. citri in China originates may from long-term climate fluctuations concomitant with recent disturbances resulting from human activity. Combined with previous data, the present work indicates that D. citri potentially entered China through two distinct invasion routes and spread within the country via three transmission paths.

Arms races with mitochondrial genome soft sweeps in a gynodioecious plant, Plantago lanceolata

Biological situations involving conflict can create arms race situations with repeated fixations of different functional variants, producing selective sweeps and lowering neutral diversity in genome regions linked to the functional locus. However, they can sometimes lead to balancing selection, potentially creating long coalescent times for sites with functionally different variants, and, if recombination occurs rarely, for extended haplotypes carrying such variants. We tested between these possibilities in a gynodioecious plant, Plantago lanceolata, in which cytoplasmic male-sterility factors conflict with nuclear restorers of male fertility. We find low mitochondrial diversity, which does not support very long-term coexistence of highly diverged mitochondrial haplotypes. Interestingly, however, we found a derived haplotype that is associated with male fertility in a restricted geographic region, and that has fixed differences from the ancestral sequence in several genes, suggesting that it did not arise very recently. Taken together, the results suggest arms race events that involved "soft" selective sweeps involving a moderately old-established haplotype, consistent with the frequency fluctuations predicted by theoretical models of gynodioecy.

Sequential horizontal gene transfers from different hosts in a widespread Eurasian parasitic plant, Cynomorium coccineum

PREMISE

Parasitic plants with large geographic ranges, and different hosts in parts of their range, may acquire horizontally transferred genes (HGTs), which might sometimes leave a footprint of gradual host and range expansion. Cynomorium coccineum, the only member of the Saxifragales family Cynomoriaceae, is a root holoparasite that occurs in water-stressed habitats from western China to the Canary Islands. It parasitizes at least 10 angiosperm families from different orders, some of them only in parts of its range. This parasite therefore offers an opportunity to trace HGTs as long as parasite-host pairs can be obtained and sequenced.

METHODS

By sequencing mitochondrial, plastid, and nuclear loci from parasite-host pairs from throughout the parasite's range and with prior information from completely assembled mitochondrial and plastid genomes, we detected 10 HGTs of five mitochondrial genes.

RESULTS

The 10 HGTs appear to have occurred sequentially as C. coccineum expanded from East to West. Molecular-clock models yield Cynomorium stem ages between 66 and 156 Myr, with relaxed clocks converging on 66-67 Myr. Chinese Sapindales, probably Nitraria, were the first source of transferred genes, followed by Iranian and Mediterranean Caryophyllales. The most recently acquired gene appears to come from a Tamarix host in the Iberian Peninsula.

CONCLUSIONS

Data on HGTs that have accumulated over the past 15 years, along with this discovery of multiple HGTs within a single widespread species, underline the need for more whole-genome data from parasite-host pairs to investigate whether and how transferred copies coexist with, or replace, native functional genes.

North-facing slopes and elevation shape asymmetric genetic structure in the range-restricted salamander Plethodon shenandoah

Species with narrow environmental tolerances are often distributed within fragmented patches of suitable habitat, and dispersal among these subpopulations can be difficult to directly observe. Genetic data can help quantify gene flow between localities, which is especially important for vulnerable species with a disjunct range. The Shenandoah salamander (Plethodon shenandoah) is a federally endangered species known only from three mountaintops in Virginia, USA. To reconstruct the evolutionary history and population connectivity of this species, we generated both mitochondrial and nuclear data using sequence capture from individuals collected across all three mountaintops. Applying population and landscape genetic methods, we found strong population structure that was independent of geographic distance. Both the nuclear markers and mitochondrial genomes indicated a deep split between the most southern population and the genetically similar central and northern populations. Although there was some mitochondrial haplotype-splitting between the central and northern populations, there was admixture in nuclear markers. This is indicative of either a recent split or current male-biased dispersal among mountain isolates. Models of landscape resistance found that dispersal across north-facing slopes at mid-elevation levels best explain the observed genetic structure among populations. These unexpected results highlight the importance of incorporating landscape features in understanding and predicting the movement and fragmentation of this range-restricted salamander species across space.

Mitochondrial Genome Fragmentation Unites the Parasitic Lice of Eutherian Mammals

Organelle genome fragmentation has been found in a wide range of eukaryotic lineages; however, its use in phylogenetic reconstruction has not been demonstrated. We explored the use of mitochondrial (mt) genome fragmentation in resolving the controversial suborder-level phylogeny of parasitic lice (order Phthiraptera). There are approximately 5000 species of parasitic lice in four suborders (Amblycera, Ischnocera, Rhynchophthirina, and Anoplura), which infest mammals and birds. The phylogenetic relationships among these suborders are unresolved despite decades of studies. We sequenced the mt genomes of eight species of parasitic lice and compared them with 17 other species of parasitic lice sequenced previously. We found that the typical single-chromosome mt genome is retained in the lice of birds but fragmented into many minichromosomes in the lice of eutherian mammals. The shared derived feature of mt genome fragmentation unites the eutherian mammal lice of Ischnocera (family Trichodectidae) with Anoplura and Rhynchophthirina to the exclusion of the bird lice of Ischnocera (family Philopteridae). The novel clade, namely Mitodivisia, is also supported by phylogenetic analysis of mt genome and cox1 gene sequences. Our results demonstrate, for the first time, that organelle genome fragmentation is informative for resolving controversial high-level phylogenies.

Phylogenetic position of the presumably extinct slender-billed curlew, Numenius tenuirostris

The high-capacity DNA analysis of museum samples opens new opportunities, associated with the investigation of extinct species evolution. Here, the complete mitochondrial genome of the presumably extinct bird species, the slender-billed curlew Numenius tenuirostris (Charadriiformes: Scolopacidae) is presented. Our results showed that mitochondrial DNA (mtDNA) is 16,705 base pairs (bp) in length and contain 13 protein-coding genes, two rRNA genes, and 22 tRNA genes. The overall base composition of the genome is 30.8% - A, 29.8% - C, 25.4% - T, 14.0% - G, and without a significant GC bias of 43.7%. Phylogenetic analyses based on the cytochrome B (cytB) gene and the whole mtDNA sequences revealed that N. tenuirostris had a close genetic relationship to Eurasian curlew (N. arquata), Far Eastern curlew (N. madagascariensis), and long-billed curlew - N. americanus. Besides, it reveals that Numenius genus is genetically distant from other Scolopacidae taxons. Together, these results provide a clear genetic perspective into the speciation process among the curlew genus members and points to a clear taxonomic position of N. tenuirostris.

Comparative Analysis of the Complete Mitochondrial Genomes for Development Application

This present research work reports the comparative analysis of the entire nucleotide sequence of mitochondrial genomes of Serranochromis robustus and Buccochromis nototaenia and phylogenetic analyses of their protein-coding genes in order to establish their phylogenetic relationship within Cichlids. The mitochondrial genomes of S. robustus and B. nototaenia are 16,583 and 16,580 base pairs long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes, and one control region (D-loop) which is 888 and 887 base pairs long, respectively, showing the same gene order and identical number of gene or regions with other well-elucidated mitogenomes of Cichlids. However, with exception of cytochrome-c oxidase subunit-1 (COX-1) gene, all the identified PCGs were initiated by ATG-codons. Structurally, 11 tRNA genes in B. nototaenia species and 9 tRNA genes in S. robustus species, folded into typical clover-leaf secondary structure created by the regions of self-complementarity within tRNA. All the 22 tRNA genes in both species lack variable loop. Moreover, 28 genes which include 12-protein-coding genes are encoded on the H-strand and the remaining 9 genes including one protein-coding gene are encoded on the L-strand. Thirteen sequences of concatenated mitochondrial protein-coding genes were aligned using MUSCLE, and the phylogenetic analyses performed using maximum likelihood and Bayesian inference showed that S. robustus and B. nototaenia had a broad phylogenetic relationship. These results may be a useful tool in resolving higher-level relationships in organisms and a useful dataset for studying the evolution of the Cichlidae mitochondrial genome, since Cichlids are well-known model species in the study of evolutionary biology, because of their extreme morphological, biogeographical, parental care behavior for eggs and larvae and phylogenetic diversities.

The complete mitogenome of Helixpomatia and the basal phylogeny of Helicinae (Gastropoda, Stylommatophora, Helicidae)

A complete mitochondrial genome of the Roman snail Helixpomatia Linnaeus, 1758 has been sequenced. The length and gene order correspond to that of other available helicid mitogenomes. We used the mitogenome sequence to reappraise the relationships among the four presumed principal groups of the helicid subfamily Helicinae. The results support the idea that the subfamily is divided between two western Palaearctic diversification centres: Iberian Peninsula and western Maghreb in the west, and Anatolia, the Aegean and Caucasus in the east. One group, the tribe Helicini, diversified in the east and the remaining three currently recognised tribes in the west. However, the exact relationships among lineages of the non-Helicini tribes could not be resolved.

The Complete Mitochondrial Genome Sequences of the Philomycus bilineatus (Stylommatophora: Philomycidae) and Phylogenetic Analysis

The mitochondrial genome (mitogenome) can provide information for phylogenetic analyses and evolutionary biology. We first sequenced, annotated, and characterized the mitogenome of Philomycus bilineatus in this study. The complete mitogenome was 14,347 bp in length, containing 13 protein-coding genes (PCGs), 23 transfer RNA genes, two ribosomal RNA genes, and two non-coding regions (A + T-rich region). There were 15 overlap locations and 18 intergenic spacer regions found throughout the mitogenome of P. bilineatus. The A + T content in the mitogenome was 72.11%. All PCGs used a standard ATN as a start codon, with the exception of cytochrome c oxidase 1 (cox1) and ATP synthase F0 subunit 8 (atp8) with TTG and GTG. Additionally, TAA or TAG was identified as the typical stop codon. All transfer RNA (tRNA) genes had a typical clover-leaf structure, except for trnS1 (AGC), trnS2 (TCA), and trnK (TTT). A phylogenetic analysis with another 37 species of gastropods was performed using Bayesian inference, based on the amino acid sequences of 13 mitochondrial PCGs. The results indicated that P. bilineatus shares a close ancestry with Meghimatium bilineatum. It seems more appropriate to reclassify it as Arionoidea rather than Limacoidea, as previously thought. Our research may provide a new meaningful insight into the evolution of P. bilineatus.

Complete Sequence of the Mitochondrial Genome of Spirometra ranarum: Comparison with S. erinaceieuropaei and S. decipiens

This study was undertaken to determine the complete mitochondrial DNA sequence and structure of the mitochondrial genome of Spirometra ranarum, and to compare it with those of S. erinaceieuropaei and S. decipiens. The aim of this study was to provide information of the species level taxonomy of Spirometra spp. using the mitochondrial genomes of 3 Spirometra tapeworms. The S. ranarum isolate originated from Myanmar. The mitochondrial genome sequence of S. ranarum was compared with that of S. erinaceieuropaei (GenBank no. KJ599680) and S. decipiens (Gen-Bank no. KJ599679). The complete mtDNA sequence of S. ranarum comprised 13,644 bp. The S. ranarum mt genome contained 36 genes comprising 12 protein-coding genes, 22 tRNAs and 2 rRNAs. The mt genome lacked the atp8 gene, as found for other cestodes. All genes in the S. ranarum mitochondrial genome are transcribed in the same direction and arranged in the same relative position with respect to gene loci as found for S. erinaceieuropaei and S. decipiens mt genomes. The overall nucleotide sequence divergence of 12 protein-coding genes between S. ranarum and S. decipiens differed by 1.5%, and 100% sequence similarity was found in the cox2 and nad6 genes, while the DNA sequence divergence of the cox1, nad1, and nad4 genes of S. ranarum and S. decipiens was 2.2%, 2.1%, and 2.6%, respectively.

Mode and Rate of Evolution of Haemosporidian Mitochondrial Genomes: Timing the Radiation of Avian Parasites

Haemosporidians are a diverse group of vector-borne parasitic protozoa that includes the agents of human malaria; however, most of the described species are found in birds and reptiles. Although our understanding of these parasites’ diversity has expanded by analyses of their mitochondrial genes, there is limited information on these genes’ evolutionary rates. Here, 114 mitochondrial genomes (mtDNA) were studied from species belonging to four genera: Leucocytozoon, Haemoproteus, Hepatocystis, and Plasmodium. Contrary to previous assertions, the mtDNA is phylogenetically informative. The inferred phylogeny showed that, like the genus Plasmodium, the Leucocytozoon and Haemoproteus genera are not monophyletic groups. Although sensitive to the assumptions of the molecular dating method used, the estimated times indicate that the diversification of the avian haemosporidian subgenera/genera took place after the Cretaceous–Paleogene boundary following the radiation of modern birds. Furthermore, parasite clade differences in mtDNA substitution rates and strength of negative selection were detected. These differences may affect the biological interpretation of mtDNA gene lineages used as a proxy to species in ecological and parasitological investigations. Given that the mitochondria are critically important in the parasite life cycle stages that take place in the vector and that the transmission of parasites belonging to particular clades has been linked to specific insect families/subfamilies, this study suggests that differences in vectors have affected the mode of evolution of haemosporidian mtDNA genes. The observed patterns also suggest that the radiation of haemosporidian parasites may be the result of community-level evolutionary processes between their vertebrate and invertebrate hosts.

Mitochondrial Genomes Provide Insights into the Phylogeny of Culicomorpha (Insecta: Diptera)

Culicomorpha is a monophyletic group containing most bloodsucking lower dipterans, including many important vectors of pathogens. However, the higher-level phylogenetic relationships within Culicomorpha are largely unresolved, with multiple competing hypotheses based on molecular sequence data. Here we sequenced four nearly complete mitochondrial (mt) genomes representing four culicomorph families, and combined these new data with published mt genomes to reconstruct the phylogenetic relationships of all eight extant culicomorph families. We estimated phylogenies using four datasets and three methods. We also used four-cluster likelihood mapping to study potential incongruent topologies supported by the different datasets and phylogenetic questions generated by the previous studies. The results showed that a clade containing Ceratopogonidae, Thaumaleidae and Simuliidae was the sister group to all other Culicomorpha; in another clade, the Dixidae was basal to the remaining four families; Chaoboridae, Corethrellidae and Culicidae formed a monophyletic group and the Chironomidae was the sister group to this clade; Culicidae and Corethrellidae were sister groups in all trees. Our study provides novel mt genome data in Culicomorpha for three new family representatives, and the resulting mt phylogenomic analysis helps to resolve the phylogeny and taxonomy of Culicomorpha.

The Mitochondrial Genomes of Neuropteridan Insects and Implications for the Phylogeny of Neuroptera

The higher-level phylogeny of Neuroptera is explored here based on the newly determined mitochondrial genomic data, with a special focus on the interfamilial relationships of this group. Despite considerable progress in our understanding of neuropteran relationships, several mutually exclusive hypotheses have come out according to morphology-based analyses and molecular sequence data. The evaluation of these hypotheses is hampered by the limited taxonomic coverage of previous studies. In this paper, we sequenced four mitochondrial genomes to improve the taxonomic sampling for families: Myrmeleontidae, Ascalaphidae and outgroup Corydalidae. Phylogenetic analyses were run using various inference methods to (1) confirm that Coniopterygidae is sister to all other Neuroptera; (2) place Hemerobiidae as sister to Chrysopidae; (3) support the monophyly of Myrmeleontiformia and define its interfamilial relationships; and (4) recover Myrmeleontidae as paraphyletic due to the nested Ascalaphidae.

Characterization of the Complete Mitochondrial Genome of Ostertagia trifurcata of Small Ruminants and its Phylogenetic Associations for the Trichostrongyloidea Superfamily

The complete mitochondrial (mt) genome of Ostertagia trifurcata, a parasitic nematode of small ruminants, has been sequenced and its phylogenetic relationship with selected members from the superfamily Trichostrongyloidea was investigated on the basis of deduced datasets of mt amino acid sequences. The entire mt genome of Ostertagia trifurcata is circular and 14,151 bp in length. It consists of a total of 36 genes comprising 12 genes coding for proteins (PCGs), 2 genes for ribosomal RNA (rRNA), 22 transfer RNA (tRNA) genes and 2 non-coding regions, since all genes are transcribed in the same direction. The phylogenetic analysis based on the concatenated datasets of predicted amino acid sequences of the 12 protein coding genes supported monophylies of the Haemonchidae, Dictyocaulidae and Molineidae families, but rejected monophylies of the Trichostrongylidae family. The complete characterization and provision of the mtDNA sequence of Ostertagia trifurcata provides novel genetic markers for molecular epidemiological investigations, systematics, diagnostics and population genetics of Ostertagia trifurcata and its correspondents.

Single-nucleotide polymorphism markers in mitochondrial genomes for identifying Varroa destructor-resistant and -susceptible strains of Apis mellifera (Hymenoptera: Apidae)

Mitogenome sequences have a high potential for possessing single-nucleotide polymorphisms (SNPs) that can be used to identify different strains of an organism bred based on maternal lines. The European honey bee, Apis mellifera ligustica (Hymenoptera: Apidae), with a high-hygienic behaviour (HHB) against the external parasitic mite Varroa destructor has been bred for several years in Korea. To distinguish this strain from low-hygienic behaviour (LHB) strains, the complete mitogenome of the two strains were sequenced using next-generation sequencing techniques to detect SNPs. The two mitogenomes with lengths of 16,449 and 16,426 base pairs (bp) in the HHB and LHB strains, respectively, contained a typical set of genes (13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes, plus one non-coding region), exhibited similar-nucleotide compositions, and had an identical gene arrangement compared to other available A. mellifera mitogenomes. The major differences between the HHB and LHB strains included the length of the intergenic spacer sequences located at the COIII and trnG junction (88 vs. 70 bp) and ND4 and ND4L junction (45 vs. 33 bp) and the presence or absence of a duplicated sequence block (CTTTTTTAAAAAAATAAAAA) in the A + T-rich region. Comparison of the mitogenome sequences from the two strains of A. m. ligustica revealed 23 SNPs in 11 protein-coding genes which were confirmed by sequencing of 10 randomly selected individuals from each strain, indicating the usefulness of these SNP markers for identifying the HHB strain of A. m. ligustica. Therefore, mitogenome sequences are a promising genome source for detecting SNP markers, particularly those in inbred female lines.