Mitogenome phylogenetics: the impact of using single regions and partitioning schemes on topology, substitution rate and divergence time estimation

Mitogenomic Insights into the Evolution, Divergence Time, and Ancestral Ranges of Coturnix Quails

The Old-World quails, Coturnix coturnix (common quail) and Coturnix japonica (Japanese quail), are morphologically similar yet occupy distinct geographic ranges. This study aimed to elucidate their evolutionary trajectory and ancestral distribution patterns through a thorough analysis of their mitochondrial genomes. Mitogenomic analysis revealed high structural conservation, identical translational mechanisms, and similar evolutionary pressures in both species. Selection analysis revealed significant evidence of positive selection across the Coturnix lineage for the nad4 gene tree owing to environmental changes and acclimatization requirements during its evolutionary history. Divergence time estimations imply that diversification among Coturnix species occurred in the mid-Miocene (13.89 Ma), and their current distributions were primarily shaped by dispersal rather than global vicariance events. Phylogenetic analysis indicates a close relationship between C. coturnix and C. japonica, with divergence estimated at 2.25 Ma during the Pleistocene epoch. Ancestral range reconstructions indicate that the ancestors of the Coturnix clade were distributed over the Oriental region. C. coturnix subsequently dispersed to Eurasia and Africa, and C. japonica to eastern Asia. We hypothesize that the current geographic distributions of C. coturnix and C. japonica result from their unique dispersal strategies, developed to evade interspecific territoriality and influenced by the Tibetan Plateau's geographic constraints. This study advances our understanding of the biogeographic and evolutionary processes leading to the diversification of C. coturnix and C. japonica, laying important groundwork for further research on this genus.

An integrative approach highlights the discrepancy in the genetic, phenotypic, and presumptive taxonomic structure of Phoxinus (Actinopterygii, Leuciscidae, Phoxininae) in Bulgaria

The present drainage network of Bulgaria is the result of a complex Neogene and Quaternary evolution. Karst, which has developed on 23% of the territory, further complicates the hydrological pattern. Fresh waters of Bulgaria drain into the Black Sea and the Aegean Sea basins and can be roughly divided into the Danube (Middle and Lower Danube), non-Danube Black Sea, East Aegean, and West Aegean hydrological regions. Phoxinus, a small leuciscid fish, has a mosaic distribution in all four of these regions, inhabiting small mountainous and semi-mountainous streams. Based on morphology, it was identified as three species, Phoxinus phoxinus in the Danube, Phoxinus strandjae in the non-Danube, and Phoxinus strymonicus in West Aegean region. Later, molecular data revealed Phoxinus csikii and Phoxinus lumaireul in the Middle Danube and P. csikii in the Lower Danube. Phoxinus has been the focus of many studies, showing a high molecular and morphological diversity, which is not entirely consistent with previous morphology-only-based taxonomic concepts. In this study, molecular (a mitochondrial marker and a nuclear marker) and morphological data from both historical and recently sampled collections were analysed to assess the applicability of the integrative approach in Phoxinus. The results showed a significant influence of the complex paleo- and recent hydrology on the currently observed genetic structure of the considered populations and species. Furthermore, the study also demonstrated a strong influence of phenotypic plasticity on the morphological analysis of Phoxinus and the lack of a clear differentiation between P. csikii and P. strandjae. A barcoded specimen was designated as neotype to fix the species named P. strandjae in the current taxonomic concept. Finally, a significant discordance between genetically delimited clades and phenotypic groups did not allow a proper delineation of the species distributed in Bulgaria, demonstrating that more molecular markers are needed for further taxonomic study of the Phoxinus complex.

The efficacy of single mitochondrial genes at reconciling the complete mitogenome phylogeny-a case study on dwarf chameleons

Although genome-scale data generation is becoming more tractable for phylogenetics, there are large quantities of single gene fragment data in public repositories and such data are still being generated. We therefore investigated whether single mitochondrial genes are suitable proxies for phylogenetic reconstruction as compared to the application of full mitogenomes. With near complete taxon sampling for the southern African dwarf chameleons (Bradypodion), we estimated and compared phylogenies for the complete mitogenome with topologies generated from individual mitochondrial genes and various combinations of these genes. Our results show that the topologies produced by single genes (ND2, ND4, ND5, COI, and COIII) were analogous to the complete mitogenome, suggesting that these genes may be reliable markers for generating mitochondrial phylogenies in lieu of generating entire mitogenomes. In contrast, the short fragment of 16S commonly used in herpetological systematics, produced a topology quite dissimilar to the complete mitogenome and its concatenation with ND2 weakened the resolution of ND2. We therefore recommend the avoidance of this 16S fragment in future phylogenetic work.

Phylogenomics of Neogastropoda: the backbone hidden in the bush

The molluscan order Neogastropoda encompasses over 15,000 almost exclusively marine species playing important roles in benthic communities and in the economies of coastal countries. Neogastropoda underwent intensive cladogenesis in early stages of diversification, generating a 'bush' at the base of their evolutionary tree, that has been hard to resolve even with high throughput molecular data. In the present study to resolve the bush, we use a variety of phylogenetic inference methods and a comprehensive exon capture dataset of 1,817 loci (79.6% data occupancy) comprising 112 taxa of 48 out of 60 Neogastropoda families. Our results show consistent topologies and high support in all analyses at (super)family level, supporting monophyly of Muricoidea, Mitroidea, Conoidea, and, with some reservations, Olivoidea and Buccinoidea. Volutoidea and Turbinelloidea as currently circumscribed are clearly paraphyletic. Despite our analyses consistently resolving most backbone nodes, three prove problematic: First, uncertain placement of Cancellariidae, as the sister group to either a Ficoidea-Tonnoidea clade, or to the rest of Neogastropoda, leaves monophyly of Neogastropoda unresolved. Second, relationships are contradictory at the base of the major 'core Neogastropoda' grouping. Third, coalescence-based analyses reject monophyly of the Buccinoidea in relation to Vasidae. We analysed phylogenetic signal of targeted loci in relation to potential biases, and we propose most probable resolutions in the latter two recalcitrant nodes. The uncertain placement of Cancellariidae may be explained by orthology violations due to differential paralog loss shortly after the whole genome duplication, which should be resolved with a curated set of longer loci.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Guidelines and quantitative standards for improved cetacean taxonomy using full mitochondrial genomes

In many organisms, especially those of conservation concern, traditional lines of evidence for taxonomic delineation, such as morphological data, are often difficult to obtain. In these cases, genetic data are often the only source of information available for taxonomic studies. In particular, population surveys of mitochondrial genomes offer increased resolution and precision in support of taxonomic decisions relative to conventional use of the control region or other gene fragments of the mitochondrial genome. To improve quantitative guidelines for taxonomic decisions in cetaceans, we build on a previous effort targeting the control region and evaluate, for whole mitogenome sequences, a suite of divergence and diagnosability estimates for pairs of recognized cetacean populations, subspecies, and species. From this overview, we recommend new guidelines based on complete mitogenomes, combined with other types of evidence for isolation and divergence, which will improve resolution for taxonomic decisions, especially in the face of small sample sizes or low levels of genetic diversity. We further use simulated data to assist interpretations of divergence in the context of varying forms of historical demography, culture, and ecology.

Time-calibrated phylogeny and full mitogenome sequence of the Galapagos sea lion (Zalophus wollebaeki) from scat DNA

Background

The Galapagos sea lion, Zalophus wollebaeki, is an endemic and endangered otariid, which is considered as a sentinel species of ecosystem dynamics in the Galapagos archipelago. Mitochondrial DNA is an important tool in phylogenetic and population genetic inference. In this work we use Illumina sequencing to complement the mitogenomic resources for Zalophus genus-the other two species employed Sanger sequencing-by a complete mitochondrial genome and a molecular clock of this species, which is not present in any case.

Materials and Methods

We used DNA obtained from a fresh scat sample of a Galapagos sea lion and shotgun-sequenced it on the Illumina NextSeq platform. The obtained raw reads were processed using the GetOrganelle software to filter the mitochondrial Zalophus DNA reads (∼16% survive the filtration), assemble them, and set up a molecular clock.

Results

From the obtained 3,511,116 raw reads, we were able to assemble a full mitogenome of a length of 16,676 bp, consisting of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNA), and two ribosomal RNAs (rRNA). A time-calibrated phylogeny confirmed the phylogenetic position of Z. wollebaeki in a clade with Z. californianus, and Z. japonicus, and sister to Z. californianus; as well as establishing the divergence time for Z. wollebaeki 0.65 million years ago. Our study illustrates the possibility of seamlessly sequencing full mitochondrial genomes from fresh scat samples of marine mammals.

Whole mitochondrial genomes reveal the relatedness of the browsing ant incursions in Australia

Global trade and human movements outspread animal species, for example ants, from their native habitats to new areas. This causes biosecurity concerns because an exotic ant might have adverse impacts on agriculture, the environment, or health; thus, incurring economic losses. The browsing ant, Lepisiota frauenfeldi, was first detected in 2013 at the Perth Airport. Since then, more discrete browsing ant infestations have been found in Perth and at the Ports of Darwin and Brisbane. This exotic ant has been deemed a significant pest in Australia and eradication efforts are underway. However, tackling this invasion requires an understanding of how these infestations are related. Are they same or separate or a combination of both? Here, we carried out a phylogenetic analysis using high-throughput sequencing data to determine their relatedness. Our results showed that each interstate incursion was separate. Furthermore, the Western Australian incursions might have two introductions. These findings are critical in devising effective biosecurity measures. However, we discovered that this information could only be revealed by analysing the whole mitochondrial genome; not by a single mitochondrial gene as typically done for species identification. Here, we sequenced 51 whole mitogenomes including three of its congener L. incisa for the first time, for tracing future infestations.

Molecular phylogeny of selected dorid nudibranchs based on complete mitochondrial genome

Dorid nudibranchs are a large group of mollusks with approximately 2,000 recorded species. Although agreement exists on the monophyletic nature of the dorid nudibranch group, the interfamily relationships of the suborder are subject to debate. Despite efforts to elucidate this issue using short molecular markers, the conclusiveness of the findings has been hindered by branching polytomy. Mitogenomes are known to be effective markers for use in phylogenetic investigations. In this study, eight mitogenomes of dorid nudibranchs were decoded and analyzed. Gene content and structure showed little change among species, reflecting the conserved mitogenomes of dorid nudibranchs. For most genes, the direction was typical for nudibranchs; nevertheless, tRNA^Cys had an inverse direction in Cadlinidae species. Phylogenetic trees based on nucleotide and amino acid datasets revealed a relatively consistent pattern of interfamily relationships with little difference for positions of Phyllidiidae and Cadlinidae. Species of Cadlinidae were clustered together and did not form a clade with Chromododidae. Additionally, Goniodorididae was sister to Aegiridae, whereas Discodoridae was sister to Dorididae. This finding was supported by tree topology test based on mitogenome data. The results of the present study indicate that complete mitogenomes are promising markers for investigating interfamily relationships among dorid nudibranchs.

Complete mitochondrial genomes and updated divergence time of the two freshwater clupeids endemic to Lake Tanganyika (Africa) suggest intralacustrine speciation

Background

The hydrogeological history of Lake Tanganyika paints a complex image of several colonization and adaptive radiation events. The initial basin was formed around 9–12 million years ago (MYA) from the predecessor of the Malagarasi–Congo River and only 5–6 MYA, its sub-basins fused to produce the clear, deep waters of today. Next to the well-known radiations of cichlid fishes, the lake also harbours a modest clade of only two clupeid species, Stolothrissatanganicae and Limnothrissamiodon. They are members of Pellonulini, a tribe of clupeid fishes that mostly occur in freshwater and that colonized West and Central-Africa during a period of high sea levels during the Cenozoic. There is no consensus on the phylogenetic relationships between members of Pellonulini and the timing of the colonization of Lake Tanganyika by clupeids.

Results

We use short-read next generation sequencing of 10X Chromium libraries to sequence and assemble the full mitochondrial genomes of S.tanganicae and L.miodon. We then use Maximum likelihood and Bayesian inference to place them into the phylogeny of Pellonulini and other clupeiforms, taking advantage of all available full mitochondrial clupeiform genomes. We identify Potamothrissaobtusirostris as the closest living relative of the Tanganyika sardines and confirm paraphyly for Microthrissa. We estimate the divergence of the Tanganyika sardines around 3.64 MYA [95% CI: 0.99, 6.29], and from P.obtusirostris around 10.92 MYA [95% CI: 6.37–15.48].

Conclusions

These estimates imply that the ancestor of the Tanganyika sardines diverged from a riverine ancestor and entered the proto-lake Tanganyika around the time of its formation from the Malagarasi–Congo River, and diverged into the two extant species at the onset of deep clearwater conditions. Our results prompt a more thorough examination of the relationships within Pellonulini, and the new mitochondrial genomes provide an important resource for the future study of this tribe, e.g. as a reference for species identification, genetic diversity, and macroevolutionary studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02085-8.

Mitogenomes resolve the phylogeography and divergence times within the endemic New Zealand Callaeidae (Aves: Passerida)

The biogeographical origins of the endemic birds of New Zealand (Aotearoa) are of great interest, particularly Palaeogene lineages such as Callaeidae, a passerine family characterized by brightly coloured wattles behind the beak and, in some cases, extreme sexual dimorphism in bill size and shape. Ancestral representatives of Callaeidae are thought to have split from their closest relatives outside New Zealand in the Oligocene, but little is known about the timing of divergences within the family. We present a fully dated molecular phylogeny of Callaeidae mitogenomes and discuss the biogeographical implications. Our results suggest that formation of Pliocene marine seaways, such as the Manawatu Strait, are likely to have played a significant role in the differentiation of North Island and South Island kōkako (Callaeas spp.) and saddlebacks/tīeke (Philesturnus spp.).

First mitochondrial genome of the Amazonian marsh rat Holochilus sciureus Wagner 1842 (Rodentia, Cricetidae, Sigmodontinae)

The Amazonian marsh rat, Holochilus sciureus, is a member of the subfamily Sigmodontinae, the second-largest subfamily of muroid rodents, with 410 species and ca. 84 genera in 12 tribes. This semiaquatic rodent is distributed in South America and is of great economic and epidemiological importance. In this study, we obtained the first mitochondrial genome of the genus Holochilus obtained from a tissue sample associated with a museum voucher specimen. The generated mitogenome sequence of H. sciureus is 16,358 bp length. It comprises a control region and a conserved set of 37 genes encoding for 2 rRNA genes, 22 tRNA genes and 13 protein-coding genes. We conducted a phylogenetic analysis that included H. sciureus and the only five other published mitochondrial genomes of this poorly studied subfamily of rodents.

The Evolutionary Dynamics and Epidemiological History of Hepatitis C Virus Genotype 6, Including Unique Strains from the Li Community of Hainan Island, China

Hepatitis C virus (HCV) is a highly diverse pathogen that frequently establishes a chronic long-term infection, but the origins and drivers of HCV diversity in the human population remain unclear. Previously unidentified strains of HCV genotype 6 (gt6) were recently discovered in chronically infected individuals of the Li ethnic group living in Baisha County, Hainan Island, China. The Li community, who were early settlers on Hainan Island, have a distinct host genetic background and cultural identity compared to other ethnic groups on the island and mainland China. In this report, we generated 33 whole virus genome sequences to conduct a comprehensive molecular epidemiological analysis of these novel gt6 strains in the context of gt6 isolates present in Southeast Asia. With the exception of one gt6a isolate, the Li gt6 sequences formed 3 novel clades from two lineages which constituted 3 newly assigned gt6 subtypes and 30 unassigned strains. Using Bayesian inference methods, we dated the most recent common ancestor for all available gt6 whole virus genome sequences to approximately 2767 BCE (95% HPD intervals, 3670 to 1397 BCE), which is far earlier than previous estimates. The substitution rate was 1.20 x 10-4 substitutions/site/year (s/s/y) and this rate varied across the genome regions, from 1.02 x 10-5 s/s/y in the 5ʹUTR region to 3.07 x 10-4 s/s/y in E2. Thus, our study on an isolated ethnic minority group within a small geographical area of Hainan Island has substantially increased the known diversity of HCV gt6, already acknowledged as the most diverse HCV genotype. The extant HCV gt6 sequences from this study were probably transmitted to the Li through at least three independent events dating perhaps from around 4000 years ago. This analysis describes deeper insight into basic aspects of HCV gt6 molecular evolution including the extensive diversity of gt6 sequences in the isolated Li ethnic group.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Extension of Mitogenome Enrichment Based on Single Long-Range PCR: mtDNAs and Putative Mitochondrial-Derived Peptides of Five Rodent Hibernators

Enriching mitochondrial DNA (mtDNA) for sequencing entire mitochondrial genomes (mitogenomes) can be achieved by single long-range PCR. This avoids interference from the omnipresent nuclear mtDNA sequences (NUMTs). The approach is currently restricted to the use of samples collected from humans and ray-finned fishes. Here, we extended the use of single long-range PCR by introducing back-to-back oligonucleotides that target a sequence of extraordinary homology across vertebrates. The assay was applied to five hibernating rodents, namely alpine marmot, Arctic and European ground squirrels, and common and garden dormice, four of which have not been fully sequenced before. Analysis of the novel mitogenomes focussed on the prediction of mitochondrial-derived peptides (MDPs) providing another level of information encoded by mtDNA. The comparison of MOTS-c, SHLP4 and SHLP6 sequences across vertebrate species identified segments of high homology that argue for future experimentation. In addition, we evaluated four candidate polymorphisms replacing an amino acid in mitochondrially encoded subunits of the oxidative phosphorylation (OXPHOS) system that were reported in relation to cold-adaptation. No obvious pattern was found for the diverse sets of mammalian species that either apply daily or multiday torpor or otherwise cope with cold. In summary, our single long-range PCR assay applying a pair of back-to-back primers that target a consensus sequence motif of Vertebrata has potential to amplify (intact) mitochondrial rings present in templates from a taxonomically diverse range of vertebrates. It could be promising for studying novel mitogenomes, mitotypes of a population and mitochondrial heteroplasmy in a sensitive, straightforward and flexible manner.

Sequencing of the Complete Mitochondrial Genome of Pingus sinensis (Spirurina: Quimperiidae): Gene Arrangements and Phylogenetic Implications

Despite several decades of intensive research on spirurine nematodes, molecular data on some of the main lineages are still absent, which makes taxonomic classification insufficiently resolved. In the present study, we sequenced the first complete mitogenome for the family Quimperiidae, belonging to P. sinensis (Spirurina: Quimperiidae), a parasite living in the intestines of snakehead (Ophiocephalus argus). The circular mitogenome is 13,874 bp long, and it contains the standard nematode gene set: 22 transfer RNAs, 2 ribosomal RNAs and 12 protein-coding genes. There are also two long non-coding regions (NCR), in addition to only 8 other intergenic regions, ranging in size from 1 to 58 bp. To investigate its phylogenetic position and study the relationships among other available Spirurina, we performed the phylogenetic analysis using Bayesian inference and maximum likelihood approaches by concatenating the nucleotide sequences of all 36 genes on a dataset containing all available mitogenomes of the suborder Spirurina from NCBI and compared with gene order phylogenies using the MLGO program. Both supported the closer relationship of Ascaridoidea to Seuratoidea than to Spiruroidea. Pingus formed a sister-group with the Cucullanus genus. The results provide a new insights into the relationships within Spirurina.

Complete mitochondrial genomes reveal robust phylogenetic signals and evidence of positive selection in horseshoe bats

BACKGROUND

In genus Rhinolophus, species in the Rhinolophus philippinensis and R. macrotis groups are unique because the horseshoe bats in these group have relatively low echolocation frequencies and flight speeds compared with other horseshoe bats with similar body size. The different characteristics among bat species suggest particular evolutionary processes may have occurred in this genus. To study the adaptive evidence in the mitochondrial genomes (mitogenomes) of rhinolophids, especially the mitogenomes of the species with low echolocation frequencies, we sequenced eight mitogenomes and used them for comparative studies of molecular phylogeny and adaptive evolution.

RESULTS

Phylogenetic analysis using whole mitogenome sequences produced robust results and provided phylogenetic signals that were better than those obtained using single genes. The results supported the recent establishment of the separate macrotis group. The signals of adaptive evolution discovered in the Rhinolophus species were tested for some of the codons in two genes (ND2 and ND6) that encode NADH dehydrogenases in oxidative phosphorylation system complex I. These genes have a background of widespread purifying selection. Signals of relaxed purifying selection and positive selection were found in ND2 and ND6, respectively, based on codon models and physicochemical profiles of amino acid replacements. However, no pronounced overlap was found for non-synonymous sites in the mitogenomes of all the species with low echolocation frequencies. A signal of positive selection for ND5 was found in the branch-site model when R. philippinensis was set as the foreground branch.

CONCLUSIONS

The mitogenomes provided robust phylogenetic signals that were much more informative than the signals obtained using single mitochondrial genes. Two mitochondrial genes that encoding proteins in the oxidative phosphorylation system showed some evidence of adaptive evolution in genus Rhinolophus and the positive selection signals were tested for ND5 in R. philippinensis. These results indicate that mitochondrial protein-coding genes were targets of adaptive evolution during the evolution of Rhinolophus species, which might have contributed to a diverse range of acoustic adaptations in this genus.

A multi-locus approach to elucidating the evolutionary history of the clingfish Tomicodon petersii (Gobiesocidae) in the Tropical Eastern Pacific

Marine species that are widely distributed in the Tropical Eastern Pacific (TEP) has served as a model for studying biogeographic patterns resulting from the effects of intraregional habitat discontinuities and oceanographic processes on the diversification and evolution of cryptobenthic reef fishes. Tomicodon petersii, a clingfish (Gobiesocidae) endemic to the TEP, is found on very shallow rocky reefs from central Mexico to northern Peru, and in the Cocos and Galapagos islands. We evaluated the effect of likely biogeographic barriers in different parts of the TEP on the diversification process of this species. We used one mitochondrial and three nuclear DNA markers from 112 individuals collected across the distribution range of T. petersii. Our phylogenetic results showed the samples constituted a monophyletic group, with three well-supported, allopatric subgroups: in the Mexican province, the Panamic province (from El Salvador to Ecuador), and the Galapagos Islands. The split between the Mexican and more southerly clades was estimated to occur at the end of the Miocene ca. 5.74 Mya, and the subsequent cladogenetic event separating the Galapagos population from the Panamic population at the junction of the Pliocene and Pleistocene, ca. 2.85 Mya. The species tree, Bayesian species delimitation tests (BPP), STACEY, and substantial genetic distances separating these three populations indicate that these three independent evolutionary units likely include two unnamed species. The cladogenetic events that promoted the formation of those genetically differentiated groups are consistent with disruptive effects on gene flow of habitat discontinuities and oceanographic processes along the mainland shoreline in the TEP and of ocean-island isolation, in conjunction with the species intrinsic life-history characteristics.

Mitochondrial genome of bronze-winged jacana (Metopidius indicus, Latham 1790)

We reported the mitochondrial genome (mitogenome) of bronze-winged jacana (Metopidius indicus, Latham 1790). The circular mitogenome was 17,208 base pairs (bp) in length, containing 13 protein-coding genes, two rRNAs, 22 tRNAs, and a non-coding control region. A DNA spacer 109 bp long was also detected between ND5 and Cytb. Phylogenetic analysis indicated that M. indicus was more closely related with the genera Himantopus, Jacana and Hydrophasianus. This annotated mitogenome reference can be utilized as a data resource for comparative mitogenomics of waders or shorebirds, with possible use in ecological and evolutionary studies.

Evolutionary history of Sundaland shrews (Eulipotyphla: Soricidae: Crocidura) with a focus on Borneo

The hyperdiverse shrew genus Crocidura is one of few small mammal genera distributed across Sundaland and all of its boundaries. This represents a rare opportunity to study the geological history of this region through the evolutionary history of these shrews. We generate a phylogeny of all recognized species of Sundaland Crocidura and show that most speciation events took place during the Pleistocene, prior to the inundation of the Sunda Shelf around 400 000 years ago. We find east–west differentiation within two separate lineages on Borneo, and that the current taxonomy of its two endemic species does not reflect evolutionary history, but ecophenotypic variation of plastic traits related to elevation. Sulawesi shrews are monophyletic, with a single notable exception: the black-footed shrew (C. nigripes). We show that the black-footed shrew diverged from its relatives on Borneo recently, suggesting a human-assisted breach of Wallace’s line. Overall, the number of Crocidura species, especially on Borneo, probably remains an underestimate.

Phylogenomic systematics of the spotted skunks (Carnivora, Mephitidae, Spilogale): Additional species diversity and Pleistocene climate change as a major driver of diversification

Four species of spotted skunks (Carnivora, Mephitidae, Spilogale) are currently recognized: Spilogale angustifrons, S. gracilis, S. putorius, and S. pygmaea. Understanding species boundaries within this group is critical for effective conservation given that regional populations or subspecies (e.g., S. p. interrupta) have experienced significant population declines. Further, there may be currently unrecognized diversity within this genus as some taxa (e.g., S. angustifrons) and geographic regions (e.g., Central America) never have been assessed using DNA sequence data. We analyzed species limits and diversification patterns in spotted skunks using multilocus nuclear (ultraconserved elements) and mitochondrial (whole mitogenomes and single gene analysis) data sets from broad geographic sampling representing all currently recognized species and subspecies. We found a high degree of genetic divergence among Spilogale that reflects seven distinct species and eight unique mitochondrial lineages. Initial divergence between S. pygmaea and all other Spilogale occurred in the Early Pliocene (∼ 5.0 million years ago). Subsequent diversification of the remaining Spilogale into an "eastern" and a "western" lineage occurred during the Early Pleistocene (∼1.5 million years ago). These two lineages experienced temporally coincident patterns of diversification at ∼0.66 and ∼0.35 million years ago into two and ultimately three distinct evolutionary units, respectively. Diversification was confined almost entirely within the Pleistocene during a timeframe characterized by alternating glacial-interglacial cycles, with the origin of this diversity occurring in northeastern Mexico and the southwestern United States of America. Mitochondrial-nuclear discordance was recovered across three lineages in geographic regions consistent with secondary contact, including a distinct mitochondrial lineage confined to the Sonoran Desert. Our results have direct consequences for conservation of threatened populations, or species, as well as for our understanding of the evolution of delayed implantation in this enigmatic group of small carnivores.

Considering Pleistocene North American wolves and coyotes in the eastern Canis origin story

The evolutionary origins and hybridization patterns of Canis species in North America have been hotly debated for the past 30 years. Disentangling ancestry and timing of hybridization in Great Lakes wolves, eastern Canadian wolves, red wolves, and eastern coyotes are most often partitioned into a 2-species model that assigns all ancestry to gray wolves and/or coyotes, and a 3-species model that includes a third, North American evolved eastern wolf genome. The proposed models address recent or sometimes late Holocene hybridization events but have largely ignored potential Pleistocene era progenitors and opportunities for hybridization that may have impacted the current mixed genomes in eastern Canada and the United States. Here, we re-analyze contemporary and ancient mitochondrial DNA genomes with Bayesian phylogenetic analyses to more accurately estimate divergence dates among lineages. We combine that with a review of the literature on Late Pleistocene Canis distributions to: (a) identify potential Pleistocene progenitors to southern North American gray wolves and eastern wolves; and (b) illuminate opportunities for ancient hybridization events. Specifically, we propose that Beringian gray wolves (C. lupus) and extinct large wolf-like coyotes (C. latrans orcutti) are likely progenitors to Mexican and Plains gray wolves and eastern wolves, respectively, and may represent a potentially unrecognized source of introgressed genomic variation within contemporary Canis genomes. These events speak to the potential origins of contemporary genomes and provide a new perspective on Canis ancestry, but do not negate current conservation priorities of dwindling wolf populations with unique genomic signatures and key ecologically critical roles.

Improving the phylogenetic resolution of Malaysian and Javan mahseer (Cyprinidae), Tor tambroides and Tor tambra: Whole mitogenomes sequencing, phylogeny and potential mitogenome markers

The true mahseer (Tor spp.) is one of the highest valued fish in the world due to its high nutritional value and great unique taste. Nevertheless, its morphological characterization and single mitochondrial gene phylogeny in the past had yet to resolve the ambiguity in its taxonomical classification. In this study, we sequenced and assembled 11 complete mahseer mitogenomes collected from Java of Indonesia, Pahang and Terengganu of Peninsular Malaysia as well as Sarawak of East Malaysia. The mitogenome evolutionary relationships among closely related Tor spp. samples were investigated based on maximum likelihood phylogenetic tree construction. Compared to the commonly used COX1 gene fragment, the complete COX1, Cytb, ND2, ND4 and ND5 genes appear to be better phylogenetic markers for genetic differentiation at the population level. In addition, a total of six population-specific mitolineage haplotypes were identified among the mahseer samples analyzed, which this offers hints towards its taxonomical landscape.

Ancient Divergence Driven by Geographic Isolation and Ecological Adaptation in Forest Dependent Sundaland Tree Squirrels

A surprising amount of hidden phylogenetic diversity exists in the small to medium size, drab colored squirrels of the genus Sundasciurus. This genus is endemic to Sundaland and the Philippines, where it is widespread. An earlier revision of this genus found that the high elevation ‘populations’ of the widespread, lowland slender squirrel (S. tenuis) were different species. Previous phylogenies based on mitochondrial cytochrome b sequences also suggested that the widespread, lowland Low’s squirrel (S. lowii) and the narrow endemic Fraternal squirrel (S. fraterculus) are not reciprocally monophyletic. Additionally, deep divergences have been identified between lineages within Low’s squirrel that date to the early Pliocene. Here we focus on evaluating the relationships and differences within and between populations of these two nominal species using whole mitochondrial genome sequences, nuclear intron sequences, and morphology. We reassess the taxonomy of this group, revalidate the species status of Robinson’s squirrel (Sundasciurus robinsoniBonhote, 1903) support the species level recognition of the Natuna squirrel (Sundasciurus natunensisThomas, 1895) and identify three other lineages that require further study. We estimate times of divergence and integrate geologic history to find that most of the divergences are pre-Pleistocene, and thus predate the Pleistocene flooding of Sundaland. Biogeographic, and ecological factors may have played a more important role than climatic factors in generating these patterns. While divergence in allopatry seems to be the main process driving speciation in lowland Sundaland squirrels (Sundasciurus), ecomorphological and behavioral adaptations in this clade suggest an important role of niche divergence.

Intra-specific comparison of mitochondrial genomes reveals host gene fragment exchange via intron mobility in Tremella fuciformis

Background

Mitochondrial genomic sequences are known to be variable. Comparative analyses of mitochondrial genomes can reveal the nature and extent of their variation.

Results

Draft mitochondrial genomes of 16 Tremella fuciformis isolates (TF01-TF16) were assembled from Illumina and PacBio sequencing data. Mitochondrial DNA contigs were extracted and assembled into complete circular molecules, ranging from 35,104 bp to 49,044 bp in size. All mtDNAs contained the same set of 41 conserved genes with identical gene order. Comparative analyses revealed that introns and intergenic regions were variable, whereas genic regions (including coding sequences, tRNA, and rRNA genes) were conserved. Among 24 introns detected, 11 were in protein-coding genes, 3 in tRNA genes, and the other 10 in rRNA genes. In addition, two mobile fragments were found in intergenic regions. Interestingly, six introns containing N-terminal duplication of the host genes were found in five conserved protein-coding gene sequences. Comparison of genes with and without these introns gave rise to the following proposed model: gene fragment exchange with other species can occur via gain or loss of introns with N-terminal duplication of the host genes.

Conclusions

Our findings suggest a novel mechanism of fungal mitochondrial gene evolution: partial foreign gene replacement though intron mobility.

Can Mitogenomes of the Northern Wheatear (Oenanthe oenanthe) Reconstruct Its Phylogeography and Reveal the Origin of Migrant Birds?

The Northern Wheatear (Oenanthe oenanthe, including the nominate and the two subspecies O. o. leucorhoa and O. o. libanotica) and the Seebohm’s Wheatear (Oenanthe seebohmi) are today regarded as two distinct species. Before, all four taxa were regarded as four subspecies of the Northern Wheatear. Their classification has exclusively been based on ecological and morphological traits, while their molecular characterization is still missing. With this study, we used next-generation sequencing to assemble 117 complete mitochondrial genomes covering O. o. oenanthe, O. o. leucorhoa and O. seebohmi. We compared the resolution power of each individual mitochondrial marker and concatenated marker sets to reconstruct the phylogeny and estimate speciation times of three taxa. Moreover, we tried to identify the origin of migratory wheatears caught on Helgoland (Germany) and on Crete (Greece). Mitogenome analysis revealed two different ancient lineages that separated around 400,000 years ago. Both lineages consisted of a mix of subspecies and species. The phylogenetic trees, as well as haplotype networks are incongruent with the present morphology-based classification. Mitogenome could not distinguish these presumed species. The genetic panmixia among present populations and taxa might be the consequence of mitochondrial introgression between ancient wheatear populations.

Mitogenomic phylogenetic analyses of Leptogorgia virgulata and Leptogorgia hebes (Anthozoa: Octocorallia) from the Gulf of Mexico provides insight on Gorgoniidae divergence between Pacific and Atlantic lineages

The use of genetics in recent years has brought to light the need to reevaluate the classification of many gorgonian octocorals. This study focuses on two Leptogorgia species-Leptogorgia virgulata and Leptogorgia hebes-from the northwestern Gulf of Mexico (GOM). We target complete mitochondrial genomes and mtMutS sequences, and integrate this data with previous genetic research of gorgonian corals to resolve phylogenetic relationships and estimate divergence times. This study contributes the first complete mitochondrial genomes for L. ptogorgia virgulata and L. hebes. Our resulting phylogenies stress the need to redefine the taxonomy of the genus Leptogorgia in its entirety. The fossil-calibrated divergence times for Eastern Pacific and Western Atlantic Leptogorgia species based on complete mitochondrial genomes shows that the use of multiple genes results in estimates of more recent speciation events than previous research based on single genes. These more recent divergence times are in agreement with geologic data pertaining to the formation of the Isthmus of Panama.

What does mitogenomics tell us about the evolutionary history of the Drosophila buzzatii cluster (repleta group)?

The Drosophila repleta group is an array of more than 100 species endemic to the “New World”, many of which are cactophilic. The ability to utilize decaying cactus tissues as breeding and feeding sites is a key aspect that allowed the successful diversification of the repleta group in American deserts and arid lands. Within this group, the Drosophila buzzatii cluster is a South American clade of seven closely related species in different stages of divergence, making them a valuable model system for evolutionary research. Substantial effort has been devoted to elucidating the phylogenetic relationships among members of the D. buzzatii cluster, including molecular phylogenetic studies that have generated ambiguous results where different tree topologies have resulted dependent on the kinds of molecular marker used. Even though mitochondrial DNA regions have become useful markers in evolutionary biology and population genetics, none of the more than twenty Drosophila mitogenomes assembled so far includes this cluster. Here, we report the assembly of six complete mitogenomes of five species: D. antonietae, D. borborema, D. buzzatii, two strains of D. koepferae and D. seriema, with the aim of revisiting phylogenetic relationships and divergence times by means of mitogenomic analyses. Our recovered topology using complete mitogenomes supports the hypothesis of monophyly of the D. buzzatii cluster and shows two main clades, one including D. buzzatii and D. koepferae (both strains), and the other containing the remaining species. These results are in agreement with previous reports based on a few mitochondrial and/or nuclear genes, but conflict with the results of a recent large-scale nuclear phylogeny, indicating that nuclear and mitochondrial genomes depict different evolutionary histories.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

Preserved collagen reveals species identity in archaeological marine turtle bones from Caribbean and Florida sites

Advancements in molecular science are continually improving our knowledge of marine turtle biology and evolution. However, there are still considerable gaps in our understanding, such as past marine turtle distributions, which can benefit from advanced zooarchaeological analyses. Here, we apply collagen fingerprinting to 130 archaeological marine turtle bone samples up to approximately 2500 years old from the Caribbean and Florida's Gulf Coast for faunal identification, finding the vast majority of samples (88%) to contain preserved collagen despite deposition in the tropics. All samples can be identified to species-level with the exception of the Kemp's ridley (Lepidochelys kempii) and olive ridley (L. olivacea) turtles, which can be separated to genus level, having diverged from one another only approximately 5 Ma. Additionally, we identify a single homologous peptide that allows the separation of archaeological green turtle samples, Chelonia spp., into two distinct groups, which potentially signifies a difference in genetic stock. The majority of the archaeological samples are identified as green turtle (Chelonia spp.; 63%), with hawksbill (Eretmochelys imbricata; 17%) and ridley turtles (Lepidochelys spp.; 3%) making up smaller proportions of the assemblage. There were no molecular identifications of the loggerhead turtle (Caretta caretta) in the assemblage despite 9% of the samples being morphologically identified as such, highlighting the difficulties in relying on morphological identifications alone in archaeological remains. Finally, we present the first marine turtle molecular phylogeny using collagen (I) amino acid sequences and find our analyses match recent phylogenies based on nuclear and mitochondrial DNA. Our results highlight the advantage of using collagen fingerprinting to supplement morphological analyses of turtle bones and support the usefulness of this technique for assessing their past distributions across the Caribbean and Florida's Gulf Coast, especially in these tropical environments where DNA preservation may be poor.

Sequencing and analysis of the complete mitogenome of Doriprismatica atromarginata (Cuvier, 1804)

Doriprismatica atromarginata is a common nudibranch in the tropical and sub-tropical Indo-Pacific region. In this study, the complete mitogenome of D. atromarginata from Vietnam was recorded for the first time. The circular mitogenome had a size of 14,421 bp and consisted of 37 genes (13 protein-coding genes, two ribosomal RNA genes, and 22 tRNA genes). Phylogenetic tree based on amino acid sequences of coding genes demonstrated that D. atromarginata has sister group relationship with the genus Chromodoris.

Hybridization between two high Arctic cetaceans confirmed by genomic analysis

In 1990, a skull from a morphologically unusual Monodontid was found in West Greenland and collected for the Natural History Museum of Denmark, University of Copenhagen. From its intermediate morphology, the skull was hypothesized to be a beluga/narwhal hybrid. If confirmed, the specimen would, to our knowledge, represent the sole evidence of hybridization between the only two toothed whale species endemic to the Arctic. Here we present genome-wide DNA sequence data from the specimen and investigate its origin using a genomic reference panel of eight belugas and eight narwhals. Our analyses reveal that the specimen is a male, first-generation hybrid between a female narwhal and a male beluga. We use stable carbon and nitrogen isotope analysis to investigate the dietary niche of the hybrid and find a higher δ¹³C value than in both belugas and narwhals, suggesting a foraging strategy unlike either parental species. These results further our understanding of the interaction between belugas and narwhals, and underscore the importance of natural history collections in monitoring changes in biodiversity. In addition, our study exemplifies how recent major advances in population genomic analyses using genotype likelihoods can provide key biological and ecological insights from low-coverage data (down to 0.05x).

Systematics, Phylogeny, and Evolution of Braconid Wasps: 30 Years of Progress

The parasitoid wasp family Braconidae is likely the second-most species-rich family in the animal kingdom. Braconid wasps are widely distributed and often encountered. They constitute one of the principal groups of natural enemies of phytophagous insects, of which many are serious pest species. The enormous biological diversification of braconid wasps has led to many homoplasies, which contributed widely to instabilities in historical classifications. Recent studies using combinations of genetic markers or total mitochondrial genomes allow for better founded groupings and will ultimately lead to a stable classification. We present the current status of the phylogenetics of the Braconidae in a historical perspective and our understanding of the effects on higher classification.

Sequencing of organellar genomes of Gymnomitrion concinnatum (Jungermanniales) revealed the first exception in the structure and gene order of evolutionary stable liverworts mitogenomes

BACKGROUND

Comparative analyses of chloroplast and mitochondrial genomes have shown that organelle genomes in bryophytes evolve slowly. However, in contrast to seed plants, the organellar genomes are yet poorly explored in bryophytes, especially among liverworts. Discovering another organellar genomes of liverwort species by sequencing provides new conclusions on evolution of bryophytes.

RESULTS

In this work, the organellar genomes of Gymnomitrion concinnatum liverwort were sequenced, assembled and annotated for the first time. The chloroplast genome displays, typical for most plants, quadripartite structure containing large single copy region (81,701 bp), two inverted repeat regions (8704 bp each) and small single copy region (20,179 bp). The gene order and content of chloroplast are very similar to other liverworts with minor differences observed. A total number of 739 and 222 RNA editing sites were predicted in chloroplast and mitochondrial genes of G. concinnatum. The mitochondrial genome gene content is also in accordance with liverworts except few alterations such as: intron loss in cox1 and atp1 genes. Nonetheless the analysis revealed that G. concinnatum mitogenome structure and gene order are rearranged in comparison with other mitogenomes of liverworts. The causes underlying such mitogenomic rearrangement were investigated and the probable model of recombination was proposed.

CONCLUSIONS

This study provide the overview of mitochondrial and chloroplast genome structure and gene order diversity of Gymnomitrion concinnatum against the background of known organellar genomes of liverworts. The obtained results cast doubt on the idea that mitogenome structure of early land plants is highly conserved as previous studies suggested. In fact is the very first case of recombination within, evolutionary stable, mitogenomes of liverworts.

Phylogenetic analysis of sooty grunter and other major freshwater fishes in the suborder Percoidei based on mitochondrial DNA

Perciformes is the largest order of fishes and vertebrates. Sooty grunter (Hephaestus fuliginosus) is an economic fish species in the Terapontidae family of Percoidei, a suborder within Perciformes. To conduct molecular-level analysis of the phylogenetic relationships between sooty grunter and major freshwater fishes in Percoidei, we analysed the entire sooty grunter mitochondrial genome sequence and obtained the mitochondrial genome information of 19 fishes from Terapontidae, Serranidae, and Centrarchidae families in Percoidei from GenBank. The complete length of the sooty grunter mitochondrial genome was 16,770 bp; it encoded 13 proteins, 2 rRNAs, 22 tRNAs, and a displacement loop (D-loop). Other than ND6 and eight tRNA genes that are encoded by the light strand, the majority of genes are encoded by the heavy strand. The sequence and distribution of sooty grunter mitochondrial-encoded genes and non-coding segment were similar to those of most vertebrates. The results of neighbour joining, maximum parsimony, and Bayesian inference analyses of the complete mitochondrial genome and six genes, including cytochrome oxidase I, cytochrome B, 12S rRNA, ND2, ND4, and ND5, were consistent. In the phylogenetic trees, fishes in Terapontidae and Centrarchidae formed monophyletic clades, whereas those in Serranidae were divided into two clades, each containing Lateolabrax and Siniperca species. Among the three freshwater fish species in Terapontidae, the freshwater Terapontidae were more closely related to jade perch than with silver perch, suggesting that freshwater Terapontidae fishes originate from marine fishes. In addition, the phylogenetic results indicated that Micropterus salmoides salmoides and Micropterus salmoides floridanus in Centrarchidae should be designated as two independent species, and Siniperca in Serranidae should be considered an independent family. The sooty grunter mitochondrial genome sequence obtained in this study could be used to conduct population genetic diversity and germplasm resource studies. Furthermore, the phylogenetic analysis results of freshwater fishes in Percoidei could provide a molecular basis for cross-breeding.

Fungal mitochondrial genomes and genetic polymorphisms

Mitochondria are the powerhouses of eukaryotic cells, responsible for ATP generation and playing a role in a diversity of cellular and organismal functions. Different from the majority of other intracellular membrane structures, mitochondria contain their own genetic materials that are capable of independent replication and inheritance. In this mini-review, we provide brief summaries of fungal mitochondrial genome structure, size, gene content, inheritance, and genetic variation. We pay special attention to the relative genetic polymorphisms of the mitochondrial vs nuclear genomes at the population level within individual fungal species. Among the 20 species/groups of species reviewed here, there is a range of variation among genes and species in the relative nuclear and mitochondrial genetic polymorphisms. Interestingly, most (15/20) showed a greater genetic diversity for nuclear genes and genomes than for mitochondrial genes and genomes, with the remaining five showing similar or slower nuclear genome genetic variations. This fungal pattern is different from the dominant pattern in animals, but more similar to that in plants. At present, the mechanisms for the variations among fungal species and the overall low level of mitochondrial sequence polymorphisms are not known. The increasing availability of population genomic data should help us reveal the potential genetic and ecological factors responsible for the observed variations.

Evolutionary progression of mitochondrial gene rearrangements and phylogenetic relationships in Strigidae (Strigiformes)

The bird mitogenome is generally considered to have a conservative genome size, consistent gene content, and similar gene order. As more mitogenomes are sequenced, mitochondrial (mt) gene rearrangements have been frequently identified among diverse birds. Within two genera (Bubo and Strix) of typical owls (Strigidae, Strigiformes), the rearrangement of the mt gene has been a subject of debate. In the current study, we first sequenced the whole mitogenomes of S. uralensis and B. scandiaca and resequenced the entire mitogenome of B. bubo. By combining our data with previously sequenced mitogenomes in Strigidae, we examined the mt gene rearrangements in the family and attempted to reconstruct the evolutionary progression of these rearrangements. The mitogenomes were then used to review the phylogenies of Strigidae. Most mitogenomes exhibited the ancestral gene order (A) in Strigidae. The ancestral gene order in the previously published mitogenome of B. bubo was found to be incorrect. We determined the mt gene order (the duplicate tRNA^Thr-CR, B) and discovered two additional mt gene orders (the duplicate tRNA^Glu-L-CR and CR, C and D) in the Bubo and Strix genera. Gene order B was likely derived from A by a tandem duplication of the region spanning from tRNA^Thr to CR. The other two modified gene orders, C and D, were likely derived from B by further degenerations or deletions of one copy of specific duplicated genes. We also preliminarily reconstructed the evolutionary progression of mt gene rearrangements and discussed maintenance of the duplicated CR in the genera. Additionally, the phylogenetic trees based on the mitogenomes supported the division of Strigidae into three subfamilies: Ninoxinae + (Surniinae + Striginae). Within the Striginae clade, the four genera formed a phylogenetic relationship: Otus + (Asio + (Bubo + Strix)). This suggests that Otus firstly diverges in their evolutionary history, and Bubo and Strix show a close relationship. B. bubo, B. blakistoni and B. scandiaca form a clade should be considered members of the same genus. The well-supported topology obtained in our Bayesian inference (BI) and maximum likelihood (ML) analyses of Strigid mitogenomes suggests that these genomes are informative for constructing phylogenetic relationships.

The Microtus voles: Resolving the phylogeny of one of the most speciose mammalian genera using genomics

Sequential rapid radiations pose some of the greatest difficulties in phylogenetics, especially when analysing only a small number of genetic markers. Given that most of the speciation events occur in quick succession at various points in time, this creates particular challenges in determining phylogenetic relationships, i.e. branching order and divergence times. With the development of high throughput sequencing, thousands of markers can now readily be used to tackle these issues. Microtus is a speciose genus currently composed of 65 species that evolved over the last 2 million years. Although it is a well-studied group, there is still phylogenetic uncertainty at various divergence levels. Building upon previous studies that generally used small numbers of mitochondrial and/or nuclear loci, in this genomic-scale study we used both mitochondrial and nuclear data to study the rapid radiation within Microtus, using partial mitogenomes and genotyping-by-sequencing (GBS) on seven species representing five Microtus subgenera and the main biogeographic ranges where this group occurs. Both types of genome (mitochondrial and nuclear) generated similar tree topologies, with a basal split of the Nearctic (M. ochrogaster) and Holarctic (M. oeconomus) species, and then a subdivision of the five Palearctic species into two subgroups. These data support the occurrence of two European radiations, one North American radiation, and a later expansion of M. oeconomus from Asia to both Europe and North America. We further resolved the positioning of M. cabrerae as sister group of M. agrestis and refute the claim that M. cabrerae should be elevated to its own genus (Iberomys). Finally, the data support ongoing speciation events, especially within M. agrestis, with high levels of genetic divergence between the three Evolutionarily Significant Units (ESUs) previously identified. Similar high levels of divergence were also found among ESUs within M. oeconomus and M. arvalis.

Mitochondrial genome diversity and population structure of two western honey bee subspecies in the Republic of South Africa

Apis mellifera capensis Eschscholtz and A.m. scutellata Lepeletier are subspecies of western honey bees that are indigenous to the Republic of South Africa (RSA). Both subspecies have invasive potential and are organisms of concern for areas outside their native range, though they are important bees to beekeepers, agriculture, and the environment where they are native. The aim of the present study was to examine genetic differentiation among these subspecies and estimate their phylogenetic relationships using complete mitochondrial genomes sequences. We used 25 individuals that were either assigned to one of the subspecies or designated hybrids using morphometric analyses. Phylogenetic analyses of mitogenome sequences by maximum likelihood (ML) and Bayesian inference identified a monophyletic RSA clade, subdivided into two clades. A haplotype network was consistent with the phylogenetic trees. However, members of both subspecies occurred in both clades, indicating that A.m. capensis and A.m. scutellata are neither reciprocally monophyletic nor do they exhibit paraphyly with one subspecies nested within the other subspecies. Furthermore, no mitogenomic features were diagnostic to either subspecies. All bees analyzed from the RSA expressed a substantial level of haplotype diversity (most samples had unique haplotypes) but limited nucleotide diversity. The number of variable codons across protein-coding genes (PCGs) differed among loci, with CO3 exhibiting the most variation and ATP6 the least.

Species delimitation of the North American orchard-spider Leucauge venusta (Walckenaer, 1841) (Araneae, Tetragnathidae)

The orchard spider, Leucauge venusta (Walckenaer, 1841) is one of the most common and abundant orb-weavers in North America. This species has a broad geographic distribution extending across tropical and temperate regions of the Americas from Canada to Brazil. Guided by a preliminary observation of the barcode gap between sequences from specimens of L. venusta collected in Florida and other North American localities, we collected across a transect through the southeastern USA to investigate the observed genetic divide. The dataset, complemented with additional samples from Mexico, and Brazil was analyzed for species delimitation using STACEY and bGMYC based on sequences from one nuclear (ITS2) and one mitochondrial marker (COI). The analyses clearly separate USA samples into two deeply divergent and geographically structured groups (north-south) which we interpret as two different species. We generated ecological niche models for these two groups rejecting a niche equivalence hypothesis for these lineages. Taxonomic changes are proposed based on these findings, Leucauge venusta is restricted to denote the northern clade, and its known distribution restricted to the USA. Leucauge argyrobapta (White, 1841) is removed from synonymy to denote the populations in Florida, Mexico and Brazil. Although the delimitation analyses suggest each of these geographic clusters within the L. argyrobapta samples represent different species, more specimens from Central and South America are needed to properly test the cohesion of L. argyrobapta populations.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.