Intraspecific rearrangement of mitochondrial genome suggests the prevalence of the tandem duplication-random loss (TDLR) mechanism in Quasipaa boulengeri

Integrative taxonomy, phylogenetics and historical biogeography of subgenus Aeschyntelus Stål, 1872 (Hemiptera: Heteroptera: Rhopalidae)

The subgenus Aeschyntelus includes six species that show variations in body color and shape, thus making it difficult to identify them based on morphological identification alone. To date, no genetic study has evaluated species within this genus. Herein, we collected 171 individuals from 90 localities of Rhopalus and employed an integrative taxonomic approach that incorporated morphological data, mitochondrial genomic data (COI, whole mitochondrial data) and nuclear genomic data (18S + 28S rRNAs, nuclear genome-wide SNPs) to delineate species boundaries. Our analyses confirmed the status of nine described species of Rhopalus and proposed the recognition of one new species known as Rhopalus qinlinganus sp. nov., which is classified within the subgenus Aeschyntelus. Discrepancies arising from nuclear and mitochondrial data suggest the presence of mito-nuclear discordance. Specifically, mitochondrial data indicated admixture within Clade A, comprising R. kerzhneri and R. latus, whereas genome-wide SNPs unambiguously identified two separate species, aligning with morphological classification. Conversely, mitochondrial data clearly distinguished Clade B- consisting of R. sapporensis into two lineages, whereas genome-wide SNPs unequivocally identified a single species. Our study also provides insights into the evolutionary history of Aeschyntelus, thus indicating that it likely originated in East Asia during the middle Miocene. The development of Aeschyntelus biodiversity in the southwestern mountains of China occurred via an uplift-driven diversification process. Our findings highlight the necessity of integrating both morphological and multiple molecular datasets for precise species identification, particularly when delineating closely related species. Additionally, it reveals the important role of mountain orogenesis on speciation within the southwestern mountains of China.

Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera)

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny

In the present study, three mitogenomes from the Bipolaris genus (Bipolaris maydis, B. zeicola, and B. oryzae) were assembled and compared with the other two reported Bipolaris mitogenomes (B. oryzae and B. sorokiniana). The five mitogenomes were all circular DNA molecules, with lengths ranging from 106,403 bp to 135,790 bp. The mitogenomes of the five Bipolaris species mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). The PCG length, AT skew and GC skew showed large variability among the 13 PCGs in the five mitogenomes. Across the 13 core PCGs tested, nad6 had the least genetic distance among the 16 Pleosporales species we investigated, indicating that this gene was highly conserved. In addition, the Ka/Ks values for all 12 core PCGs (excluding rps3) were < 1, suggesting that these genes were subject to purifying selection. Comparative mitogenomic analyses indicate that introns were the main factor contributing to the size variation of Bipolaris mitogenomes. The introns of the cox1 gene experienced frequent gain/loss events in Pleosporales species. The gene arrangement and collinearity in the mitogenomes of the five Bipolaris species were almost highly conserved within the genus. Phylogenetic analysis based on combined mitochondrial gene datasets showed that the five Bipolaris species formed well-supported topologies. This study is the first report on the mitogenomes of B. maydis and B. zeicola, as well as the first comparison of mitogenomes among Bipolaris species. The findings of this study will further advance investigations into the population genetics, evolution, and genomics of Bipolaris species.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

The first two mitochondrial genomes from Apiotrichum reveal mitochondrial evolution and different taxonomic assignment of Trichosporonales

Apiotrichum is a diverse anamorphic basidiomycetous yeast genus, and its mitogenome characterization has not been revealed. In this study, we assembled two Apiotrichum mitogenomes and compared them with mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. The mitogenomes of Apiotrichum gracile and A. gamsii comprised circular DNA molecules, with sizes of 34,648 bp and 38,096 bp, respectively. Intronic regions were found contributed the most to the size expansion of A. gamsii mitogenome. Comparative mitogenomic analysis revealed that 6.85-38.89% of nucleotides varied between tRNAs shared by the two Apiotrichum mitogenomes. The GC content of all core PCGs in A. gamsii was lower than that of A. gracile, with an average low value of 4.97%. The rps3 gene differentiated the most among Agaricomycotina, Pucciniomycotina and Ustilaginomycotina species, while nad4L gene was the most conserved in evolution. The Ka/Ks values for cob and rps3 genes were > 1, indicating the two genes may be subjected to positive selection in Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. Frequent intron loss/gain events and potential intron transfer events have been detected in evolution of Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. We further detected large-scale gene rearrangements between the 19 mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina, and fifteen of the 17 mitochondrial genes shared by Apiotrichum varied in gene arrangements. Phylogenetic analyses based on maximum likelihood and Bayesian inference methods using a combined mitochondrial gene dataset revealed different taxonomic assignment of two Apiotrichum species, wherein A. gamsii had a more closely relationship with Trichosporon asahii. This study served as the first report on mitogenomes from the genus Apiotrichum, which promotes the understanding of evolution, genomics, and phylogeny of Apiotrichum.

The Complete Mitochondrial Genome of Box Tree Moth Cydalima perspectalis and Insights into Phylogenetics in Pyraloidea

To resolve and reconstruct phylogenetic relationships within Pyraloidea based on molecular data, the mitochondrial genome (mitogenome) was widely applied to understand phylogenetic relations at different taxonomic levels. In this research, a complete mitogenome of Cydalima perspectalis was recorded, and the phylogenetic position of C. perspectalis was inferred based on the sequence in combination with other available sequence data. According to the research, the circular mitochondrial genome is 15,180 bp in length. It contains 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), 13 typical protein-coding genes (PCGs), and a non-coding control region. The arrangement of a gene of the C. perspectalis mitogenome is not the same as the putative ancestral arthropod mitogenome. All of the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is undertaken by CGA. Five genes have incomplete stop codons that contain only ‘T’. All tRNA genes display a typical clover–leaf structure of mitochondrial tRNA, except for trnS1 (AGN). The control region contained an ‘ATAGG(A)’-like motif followed by a poly-T stretch. Based on the mitochondrial data, phylogenetic analysis within Pyraloidea was carried out using Bayesian inference (BI) and maximum likelihood (ML) analyses. Phylogenetic analysis showed that C. perspectalis is more closely related to Pygospila tyres within Spilomelinae than those of Crambidae and Pyraloidea.

First two mitochondrial genomes for the order Filobasidiales reveal novel gene rearrangements and intron dynamics of Tremellomycetes

In the present study, two mitogenomes from the Filobasidium genus were assembled and compared with other Tremellomycetes mitogenomes. The mitogenomes of F. wieringae and F. globisporum both comprised circular DNA molecules, with sizes of 27,861 bp and 71,783 bp, respectively. Comparative mitogenomic analysis revealed that the genetic contents, tRNAs, and codon usages of the two Filobasidium species differed greatly. The sizes of the two Filobasidium mitogenomes varied greatly with the introns being the main factor contributing to mitogenome expansion in F. globisporum. Positive selection was observed in several protein-coding genes (PCGs) in the Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina species, including cob, cox2, nad2, and rps3 genes. Frequent intron loss/gain events were detected to have occurred during the evolution of the Tremellomycetes mitogenomes, and the mitogenomes of 17 species from Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina have undergone large-scale gene rearrangements. Phylogenetic analyses based on Bayesian inference and the maximum likelihood methods using a combined mitochondrial gene set generated identical and well-supported phylogenetic trees, wherein Filobasidium species had close relationships with Trichosporonales species. This study, which is the first report on mitogenomes from the order Filobasidiales, provides a basis for understanding the genomics, evolution, and taxonomy of this important fungal group.

Comparative mitochondrial genome analyses reveal conserved gene arrangement but massive expansion/contraction in two closely related Exserohilum pathogens

Exserohilum turcicum and E. rostratum, two closely related fungal species, are both economically important pathogens but have quite different target hosts (specific to plants and cross-kingdom infection, respectively). In the present study, complete circular mitochondrial genomes of the two Exserohilum species were sequenced and de novo assembled, which mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). Comparative analyses indicated that these two fungi had significant mitogenomic collinearity and consistent mitochondrial gene arrangement, yet with vastly different mitogenome sizes, 264,948 bp and 64,620 bp, respectively. By contrast with the 17 introns containing 17 intronic ORFs (one-to-one) in the E. rostratum mitogenome, E. turcicum involved far more introns (70) and intronic ORFs (126), which was considered as the main contributing factors of their mitogenome expansion/contraction. Within the generally intron-rich gene cox1, a total of 18 and 10 intron position classes (Pcls) were identified separately in the two mitogenomes. Moreover, 16.16% and 10.85% ratios of intra-mitogenomic repetitive regions were detected in E. turcicum and E. rostratum, respectively. Based on the combined mitochondrial gene dataset, we established a well-supported topology of phylogeny tree of 98 ascomycetes, implying that mitogenomes may act as an effective molecular marker for fungal phylogenetic reconstruction. Our results served as the first report on mitogenomes in the genus Exserohilum, and would have significant implications in understanding the origin, evolution and pathogenic mechanisms of this fungal lineage.

Comparative analysis of the complete mitochondrial genomes of four cordyceps fungi

Cordyceps is a large group of entomogenous, medicinally important fungi. In this study, we sequenced, assembled, and annotated the entire mitochondrial genome of Ophiocordyceps xuefengensis, in addition to comparing it against other three complete cordyceps mitogenomes that were previously published. Comparative analysis indicated that the four complete mitogenomes are all composed of circular DNA molecules, although their sizes significantly differ due to high variability in intron and intergenic region sizes in the Ophiocordyceps sinensis and O. xuefengensis mitogenomes. All mitogenomes contain 14 conserved genes and two ribosomal RNA genes, but varying numbers of tRNA introns. The Ka/Ks ratios for all 14 PCGs and rps3 were all less than 1, indicating that these genes have been subject to purifying selection. Phylogenetic analysis was conducted using concatenated amino acid and nucleotide sequences of the 14 PCGs and rps3 using two different methods (Maximum Likelihood and Bayesian analysis), revealing highly supported relationships between O. xuefengensis and other Ophiocordyceps species, in addition to a close relationship with O. sinensis. Further, the analyses indicated that cox1 and rps3 play important roles in population differentiation. These mitogenomes will allow further study of the population genetics, taxonomy, and evolutionary biology of medicinally important cordyceps species.

Comparative genomic analysis of vertebrate mitochondrial reveals a differential of rearrangements rate between taxonomic class

Vertebrate mitochondrial genomes have been extensively studied for genetic and evolutionary purposes, these are normally believed to be extremely conserved, however, different cases of gene rearrangements have been reported. To verify the level of rearrangement and the mitogenome evolution, we performed a comparative genomic analysis of the 2831 vertebrate mitochondrial genomes representing 12 classes available in the NCBI database. Using a combination of bioinformatics methods, we determined there is a high number of errors in the annotation of mitochondrial genes, especially in tRNAs. We determined there is a large variation in the proportion of rearrangements per gene and per taxonomic class, with higher values observed in Actinopteri, Amphibia and Reptilia. We highlight that these are results for currently available vertebrate sequences, so an increase in sequence representativeness in some groups may alter the rearrangement rates, so in a few years it would be interesting to see if these rates are maintained or altered with the new mitogenome sequences. In addition, within each vertebrate class, different patterns in rearrangement proportion with distinct hotspots in the mitochondrial genome were found. We also determined that there are eleven convergence events in gene rearrangement, nine of which are new reports to the scientific community.

Mountain frog species losing out to climate change around the Sichuan Basin

Amphibians are particularly vulnerable to climate changes that are expected to cause habitat fragmentation and loss and, ultimately, local extirpations. However, little is known about how the interaction between climate change and fragmentation may impede the ability of amphibians to adapt to climate change. Here, we used the iconic mountain frog Quasipaa boulengeri as an indicator species to extrapolate climate-driven shifts in its habitat availability and connectivity in central and southern China according to the minimum and maximum representative concentration pathways. The models projected an average habitat loss of 36%-71% and the in situ and ex situ climate-change refugia to be 29%-64% and 5%-18% of the present-day suitable habitats, respectively. An increase in habitat fragmentation was reflected in a 51% decrease in core patch size, a 9% increase in the mean least-cost path (LCP) length, and a 19% increase in the cost-weighted distance. These climate-driven shifts varied spatially around the Sichuan Basin, with those in the southeast of the Basin being the most pronounced habitat and connectivity losses and those along the Basin being relatively optimistic. The effectiveness of refugia may only be maintained through a narrow passageway along the southern Sichuan Basin because of the presence of LCPs over time. Our results emphasize the need to understand how climate change and connectivity will jointly affect the distribution of mountain amphibians and to accordingly adopt conservation strategies. Further, our findings highlight the importance of identifying and preserving climate-change refugia and habitat connectivity for species persistence and conservation planning.

Comparative Analysis of Mitochondrial Genomes in Two Subspecies of the Sunwatcher Toad-Headed Agama (Phrynocephalus helioscopus): Prevalent Intraspecific Gene Rearrangements in Phrynocephalus

Intraspecific rearrangements of mitochondrial genomes are rarely reported in reptiles, even in vertebrates. The sunwatcher toad-headed agama, Phryncoephalus helioscopus, can serve as an excellent model for investigating the dynamic mitogenome structure at intraspecific level. To date, seven subspecies of P. helioscopus are well recognized, but little is known about the mitogenomic evolution among different subspecies. In this study, complete mitogenomes of subspecies P. helioscopus varius II and P. helioscopus cameranoi were determined by next-generation sequencing, and another P. helioscopus varius I retrieved from GenBank was compiled for comparative analysis. The nucleotide composition and the codon usage are similar to those previously published from toad-headed agamas. P. helioscopus varius II and P. helioscopus cameranoi have 23 tRNA genes, including standard 22 tRNA genes and one extra tRNA-Phe (tRNA-Phe duplication). Gene order and phylogenetic analyses in the genus Phrynocephalus support prevalent intraspecific gene rearrangement in P. helioscopus and other congener species including P. erythrurus, P. vlangalii, and P. forsythii. Six different mitochondrial gene arrangements are observed in Phrynocephalus. Overall, the occurrence of rearrangements may result from multiple independent structural dynamic events. The split of the two subspecies in P. helioscopus was dated at approximately 2.34 million years ago (Ma). Two types of gene rearrangements are found in the three mitogenomes of P. helioscopus, and this intraspecific rearrangement phenomenon can be explained by the tandem duplication/random loss (TDRL) model. Post duplication, the alternative loss types can occur in 0.23–0.72 Ma, suggesting that the duplication and fixation of these rearrangements can occur quite quickly. These findings highlight the need for more mitogenomes at the population level in order to better understand the potentially rampant intraspecific mitogenomic reorganization in Phrynocephalus.

Comparative Mitogenomic Analysis Reveals Dynamics of Intron Within and Between Tricholoma Species and Phylogeny of Basidiomycota

The genus of Tricholoma is a group of important ectomycorrhizal fungi. The overlapping of morphological characteristics often leads to the confusion of Tricholoma species classification. In this study, the mitogenomes of five Tricholoma species were sequenced based on the next-generation sequencing technology, including T. matsutake SCYJ1, T. bakamatsutake, T. terreum, T. flavovirens, and T. saponaceum. These five mitogenomes were all composed of circular DNA molecules, with sizes ranging from 49,480 to 103,090 bp. Intergenic sequences were considered to be the main factor contributing to size variations of Tricholoma mitogenomes. Comparative mitogenomic analysis showed that the introns of the Agaricales mitogenome experienced frequent loss/gain events. In addition, potential gene transfer was detected between the mitochondrial and nuclear genomes of the five species of Tricholoma. Evolutionary analysis showed that the rps3 gene of the Tricholoma species was under positive selection or relaxed selection in the evolutionary process. In addition, large-scale gene rearrangements were detected between some Tricholoma species. Phylogenetic analysis using the Bayesian inference and maximum likelihood methods based on a combined mitochondrial gene set yielded identical and well-supported tree topologies. This study promoted the understanding of the genetics, evolution, and phylogeny of the Tricholoma genus and related species.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

Quantification and evolution of mitochondrial genome rearrangement in Amphibians

BACKGROUND

Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution.

RESULTS

Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation.

CONCLUSION

It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus

In the present study, the complete mitogenome of Turbinellus floccosus was sequenced, assembled, and compared with other basidiomycete mitogenomes. The mitogenome of T. floccosus consists of a circular DNA molecule, with a size of 62,846 bp. Gene arrangement analysis indicated that large-scale gene rearrangements occurred in the levels of family and genus of basidiomycete species, and the mitogenome of T. floccosus contained a unique gene order. A significant correlation between the number of introns and the mitochondrial genome size of Basidiomycota were detected (P < 0.01). A total of 896 introns were detected in the core protein-coding genes (PCGs) of 74 basidiomycete species, and the cox1 gene was the largest host gene of basidiomycete introns. Intron position class (Pcls) P383 in the cox1 gene was the most common intron in Basidiomycota, which distributed in 40 of 74 basidiomycete species. In addition, frequent intron loss/gain events were detected in basidiomycete species. More than 50% of bases around insertion sites (- 15 bp to 15 bp) of Pcls from different species were conservative, indicating site preferences of intron insertions in Basidiomycota. Further analysis showed that 76.09% of introns tended to insert downstream to a T base in Basidiomycota. Phylogenetic analysis for 74 basidiomycetes indicated mitochondrial genes are effective molecular markers for phylogeny of basidiomycetes. The study served as the first report on the mitogenome from the family Gomphaceae, which will help to understand the intron origin and evolution in Basidiomycota. KEY POINTS: • The mitogenome of Turbinellus floccosus had a unique gene arrangement. • Intron loss/gain events were detected in the 74 basidiomycete species. • Introns tend to insert downstream of a T base in basidiomycete mitogenomes.

Comparative mitochondrial genome analysis reveals intron dynamics and gene rearrangements in two Trametes species

Trametes species are efficient wood decomposers that are widespread throughout the world. Mitogenomes have been widely used to understand the phylogeny and evolution of fungi. Up to now, two mitogenomes from the Trametes genus have been revealed. In the present study, the complete mitogenomes of two novel Trametes species, Trametes versicolor and T. coccinea, were assembled and compared with other Polyporales mitogenomes. Both species contained circular DNA molecules, with sizes of 67,318 bp and 99,976 bp, respectively. Comparative mitogenomic analysis indicated that the gene number, length and base composition varied between the four Trametes mitogenomes we tested. In addition, all of the core protein coding genes in Trametes species were identified and subjected to purifying selection. The mitogenome of T. coccinea contained the largest number of introns among the four Trametes species tested, and introns were considered the main factors contributing to size variations of Polyporales. Several novel introns were detected in the Trametes species we assembled, and introns identified in Polyporales were found to undergo frequent loss/gain events. Large-scale gene rearrangements were detected between closely related Trametes species, including gene inversions, insertions, and migrations. A well-supported phylogenetic tree for 77 Basidiomycetes was obtained based on the combined mitochondrial gene set using 2 phylogenetic inference methods. The results showed that mitochondrial genes are effective molecular markers for understanding the phylogeny of Basidiomycetes. This study is the first to report the mitogenome rearrangement and intron dynamics of Trametes species, which shed light on the evolution of Trametes and other related species.

The mitogenomes of two saprophytic Boletales species (Coniophora) reveals intron dynamics and accumulation of plasmid-derived and non-conserved genes

The order Boletales is a group of fungi with complex life styles, which include saprophytic and ectomycorrhizal mushroom-forming fungi. In the present study, the complete mitogenomes of two saprophytic Boletales species, Coniophora olivacea, and C. puteana, were assembled and compared with mitogenomes of ectomycorrhizal Boletales. Both mitogenomes comprised circular DNA molecules with sizes of 78,350 bp and 79,655 bp, respectively. Comparative mitogenomic analysis indicated that the two saprophytic Boletales species contained more plasmid-derived (7 on average) and unknown functional genes (12 on average) than the four ectomycorrhizal Boletales species previously reported. In addition, the core protein coding genes, nad2 and rps3, were found to be subjected to positive selection pressure between some Boletales species. Frequent intron gain/loss events were detected in Boletales and Basidiomycetes, and several novel intron classes were found in two Coniophora species. A total of 33 introns were detected in C. olivacea, and most were found to have undergone contraction in the C. olivacea mitogenome. Mitochondrial genes of Coniophora species were found to have undergone large-scale gene rearrangements, and the accumulation of intra-genomic repeats in the mitogenome was considered as one of the main contributing factors. Based on combined mitochondrial gene sets, we obtained a well-supported phylogenetic tree for 76 Basidiomycetes, demonstrating the utility of mitochondrial gene analysis for inferring Basidiomycetes phylogeny. The study served as the first report on the mitogenomes of the family Coniophorineae, which will help to understand the origin and evolution patterns of Boletales species with complex lifestyles.

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6-trnE-Cytb-D-loop-trnL2 (CUN)-ND5-D-loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.

The 256 kb mitochondrial genome of Clavaria fumosa is the largest among phylum Basidiomycota and is rich in introns and intronic ORFs

In the present study, the complete mitogenome of Clavaria fumosa, was sequenced, assembled, and compared. The complete mitogenome of C. fumosa is 256,807 bp in length and is the largest mitogenomes among all Basidiomycota mitogenomes reported. Comparative mitogenomic analysis indicated that the C. fumosa mitogenome contained the most introns and intronic ORFs among all fungal mitogenomes. Large intergenic regions, intronic regions, accumulation of repeat sequences and plasmid-derived genes together promoted the size expansion of the C. fumosa mitogenome. In addition, the rps3 gene was found subjected to positive selection between some Agaricales species. We found frequent intron gain/loss events in Agaricales mitogenomes, and four novel intron classes were detected in the C. fumosa mitogenome. Large-scale gene rearrangements were found occurred in Agaricales species and the C. fumosa mitogenome had a unique gene arrangement which differed from other Agaricales species. Phylogenetic analysis for 76 Basidiomycetes based on combined mitochondrial gene sets indicated that mitochondrial genes could be used as effective molecular markers for reconstructing evolution of Basidiomycota. The study served as the first report on the mitogenomes of the family Clavariaceae, which will promote the understanding of the genetics, evolution and taxonomy of C. fumosa and related species.

The First Mitochondrial Genome for Geastrales (Sphaerobolus stellatus) Reveals Intron Dynamics and Large-Scale Gene Rearrangements of Basidiomycota

In this study, the mitogenome of artillery fungus, Sphaerobolus stellatus, was assembled and compared with other Basidiomycota mitogenomes. The Sphaerobolus stellatus mitogenome was composed of circular DNA molecules, with a total size of 152,722 bp. Accumulation of intergenic and intronic sequences contributed to the Sphaerobolus stellatus mitogenome becoming the fourth largest mitogenome among Basidiomycota. We detected large-scale gene rearrangements in Basidiomycota mitogenomes, and the Sphaerobolus stellatus mitogenome contains a unique gene order. The quantity and position classes of intron varied between 75 Basidiomycota species we tested, indicating frequent intron loss/gain events occurred in the evolution of Basidiomycota. A novel intron position classes (P1281) was detected in the Sphaerobolus stellatus mitogenome, without any homologous introns from other Basidiomycota species. A pair of fragments with a total length of 9.12 kb in both the nuclear and mitochondrial genomes of Sphaerobolus stellatus was detected, indicating possible gene transferring events. Phylogenetic analysis based on the combined mitochondrial gene set obtained well-supported tree topologies (Bayesian posterior probabilities ≥ 0.99; bootstrap values ≥98). This study served as the first report on the mitogenome from the order Geastrales, which will promote the understanding of the phylogeny, population genetics, and evolution of the artillery fungus, Sphaerobolus stellatus.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

In this study, the mitogenomes of two Paxillus species were assembled, annotated and compared. The two mitogenomes of Paxillus involutus and P. rubicundulus comprised circular DNA molecules, with the size of 39,109 bp and 41,061 bp, respectively. Evolutionary analysis revealed that the nad4L gene had undergone strong positive selection in the two Paxillus species. In addition, 10.64 and 36.50% of the repetitive sequences were detected in the mitogenomes of P. involutus and P. rubicundulus, respectively, which might transfer between mitochondrial and nuclear genomes. Large-scale gene rearrangements and frequent intron gain/loss events were detected in 61 basidiomycete species, which revealed large variations in mitochondrial organization and size in Basidiomycota. In addition, the insertion sites of the basidiomycete introns were found to have a base preference. Phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies, indicating that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships of Basidiomycota. This study is the first report on the mitogenomes of Paxillus, which will promote a better understanding of their contrasted ecological strategies, molecular evolution and phylogeny of these important ectomycorrhizal fungi and related basidiomycete species.

The Origin and Evolution of Chromosomal Reciprocal Translocation in Quasipaa boulengeri (Anura, Dicroglossidae)

Chromosomal rearrangements have long fascinated evolutionary biologists for being widely implicated in causing genetic differentiation. Suppressed recombination has been demonstrated in various species with inversion; however, there is controversy over whether such recombination suppression would facilitate divergence in reciprocal translocation with reduced fitness. In this study, we used the spiny frog, Quasipaa boulengeri, whose western Sichuan Basin populations exhibit translocation polymorphisms, to test whether the genetic markers on translocated (rearranged) or normal chromosomes have driven this genetic differentiation. We also investigated its overall genetic structure and the possibility of chromosomal fixation. Whole-chromosome painting and genetic structure clustering suggested a single origin of the translocation polymorphisms, and high-throughput sequencing of rearranged chromosomes isolated many markers with known localizations on chromosomes. Using these markers, distinct patterns of gene flow were found between rearranged and normal chromosomes. Genetic differentiation was only found in the translocated chromosomes, not in normal chromosomes or the mitochondrial genome. Hybrid unfitness cannot explain the genetic differentiation, as then the differentiation would be observed throughout the whole genome. Our results suggest that suppressed recombination drives genetic differentiation into a balanced chromosomal polymorphism. Mapping to a reference genome, we found that the region of genetic differentiation covered a wide range of translocated chromosomes, not only in the vicinity of chromosomal breakpoints. Our results imply that the suppressed recombination region could be extended by accumulation of repetitive sequences or capture of alleles that are adapted to the local environment, following the spread and/or fixation of chromosomal rearrangement.

Exceptional Enlargement of the Mitochondrial Genome Results from Distinct Causes in Different Rain Frogs (Anura: Brevicipitidae: Breviceps)

The mitochondrial (mt) genome of the bushveld rain frog (Breviceps adspersus, Brevicipitidae, Afrobatrachia) is the largest (28.8 kbp) among the vertebrates investigated to date. The major cause of genome size enlargement in this species is the duplication of multiple genomic regions. To investigate the evolutionary lineage, timing, and process of mt genome enlargement, we sequenced the complete mtDNAs of two congeneric rain frogs, B. mossambicus and B. poweri. The mt genomic organization, gene content, and gene arrangements of these two rain frogs are very similar to each other but differ from those of B. adspersus. The B. mossambicus mt genome (22.5 kbp) does not differ significantly from that of most other afrobatrachians. In contrast, the B. poweri mtDNA (28.1 kbp) is considerably larger: currently the second largest among vertebrates, after B. adspersus. The main causes of genome enlargement differ among Breviceps species. Unusual elongation (12.5 kbp) of the control region (CR), a single major noncoding region of the vertebrate mt genome, is responsible for the extremely large mt genome in B. poweri. Based on the current Breviceps phylogeny and estimated divergence age, it can be concluded that the genome enlargements occurred independently in each species lineage within relatively short periods. Furthermore, a high nucleotide substitution rate and relaxation of selective pressures, which are considered to be involved in changes in genome size, were also detected in afrobatrachian lineages. Our results suggest that these factors were not direct causes but may have indirectly affected mt genome enlargements in Breviceps.

The complete mitochondrial genomes of two model ectomycorrhizal fungi (Laccaria): features, intron dynamics and phylogenetic implications

Laccaria amethystine and L. bicolor have served as model species for studying the life history and genetics of ectomycorrhizal fungi. However, the characterizations and variations of their mitogenomes are still unknown. In the present study, the mitogenomes of the two Laccaria species were assembled, annotated, and compared. The two mitogenomes of L. amethystine and L. bicolor comprised circular DNA molecules, with the sizes of 65,156 bp and 95,304 bp, respectively. Genome collinearity analysis revealed large-scale gene rearrangements between the two Laccaria species. Comparative mitogenome analysis indicated the introns of cox1 genes in Agaricales experienced frequent lost/gain eveants, which promoted the organization and size variations in Agaricales mitogenomes. Evolutionary analysis indicated the core protein-coding genes in the two mitogenomes were subject to strong pressure of purifying selection. Phylogenetic analysis using the Bayesian inference (BI) and Maximum likelihood (ML) methods based on a combined mitochondrial gene set resulted in identical and well-supported tree topologies, wherein the two Laccaria species were most closely related to Coprinopsis cinerea. This study severed as the first study on the mitogenomes of Laccaria species, which promoted a comprehensive understanding of the genetics and evolution of the model ectomycorrhizal fungi.

The Complete Mitogenome of Pyrrhocoris tibialis (Hemiptera: Pyrrhocoridae) and Phylogenetic Implications

We determined the complete mitogenome of Pyrrhocoris tibialis (Hemiptera: Heteroptera: Pyrrhocoridae) to better understand the diversity and phylogeny within Pentatomomorpha, which is the second largest infra-order of Heteroptera. Gene content, gene arrangement, nucleotide composition, codon usage, ribosomal RNA (rRNA) structures, and sequences of the mitochondrial transcription termination factor were well conserved in Pyrrhocoroidea. Different protein-coding genes have been subject to different evolutionary rates correlated with the G + C content. The size of control regions (CRs) was highly variable among mitogenomes of three sequenced Pyrrhocoroidea species, with the P. tibialis CR being the largest. All the transfer RNA genes found in Pyrrhocoroidea had the typical clover leaf secondary structure, except for trnS1 (AGN), which lacked the dihydrouridine arm and possessed an unusual anticodon stem (9 bp vs. the normal 5 bp). A total of three different phylogenetic relationships among the five super-families of Pentatomomorpha were obtained using three analytical methods (MrBayes and RAxML under site-homogeneous models and PhyloBayes under a site-heterogeneous CAT + GTR model) and two mitogenomic datasets (nucleotides and amino acids). The tree topology test using seven methods statistically supported a phylogeny of (Aradoidea + (Pentatomoidea + (Lygaeoidea + (Pyrrhocoroidea + Coreoidea)))) as the best topology, as recognized by both RAxML and MrBayes based on the two datasets.

Comparative Mitochondrial Genome Analysis of Two Ectomycorrhizal Fungi (Rhizopogon) Reveals Dynamic Changes of Intron and Phylogenetic Relationships of the Subphylum Agaricomycotina

In the present study, we assembled and compared two mitogenomes from the Rhizopogon genus. The two mitogenomes of R. salebrosus and R. vinicolor comprised circular DNA molecules, with the sizes of 66,704 bp and 77,109 bp, respectively. Comparative mitogenome analysis indicated that the length and base composition of protein coding genes (PCGs), rRNA genes and tRNA genes varied between the two species. Large fragments aligned between the mitochondrial and nuclear genomes of both R. salebrosus (43.41 kb) and R. vinicolor (12.83 kb) indicated that genetic transfer between mitochondrial and nuclear genomes has occurred over evolutionary time of Rhizopogon species. Intronic regions were found to be the main factors contributing to mitogenome expansion in R. vinicolor. Variations in the number and type of introns in the two mitogenomes indicated that frequent intron loss/gain events occurred during the evolution of Rhizopogon species. Phylogenetic analyses based on Bayesian inference (BI) and Maximum likelihood (ML) methods using a combined mitochondrial gene set yielded identical and well-supported tree topologies, wherein Rhizopogon species showed close relationships with Agaricales species. This is the first study of mitogenomes within the genus Rhizopogon, and it provides a basis for understanding the evolution and differentiation of mitogenomes from the ectomycorrhizal fungal genus.

The Complete Mitochondrial Genome and Expression Profile of Mitochondrial Protein-Coding Genes in the Bisexual and Parthenogenetic Haemaphysalis longicornis

The tick Haemaphysalis longicornis is widely distributed in eastern Asia, New Zealand and Australia, and is well-known as a vector of multiple zoonotic pathogens. This species exhibits two reproductive strategies, bisexual and obligate parthenogenetic reproduction. Hence, in the current study, the complete mitochondrial genomes of the bisexual and parthenogenetic populations were assembled and analyzed, and the expression of the mitochondrial protein-coding genes was evaluated and compared between the two reproductive populations. The results indicated that the length of the mitochondrial genomes of the two reproductive populations is 14,694 and 14,693 bp in the bisexual and parthenogenetic populations, respectively. The AT content in the mitochondrial genome of the bisexual and obligate parthenogenetic population reached 77.22 and 77.34%, respectively. The phylogenetic tree was constructed combining 13 protein-coding genes, which showed that the genetic distance between the bisexual and parthenogenetic populations was less than that between the subspecies. The expression of the mitochondrial protein-coding genes was quantitatively analyzed at different feeding status for the bisexual and parthenogenetic populations, and the results showed significant differences in the expression patterns of these genes, suggesting that they might trigger specific energy utilization mechanisms due to their different reproductive strategies and environmental pressures.

Suppressed Recombination of Sex Chromosomes Is Not Caused by Chromosomal Reciprocal Translocation in Spiny Frog (Quasipaa boulengeri)

Chromosome rearrangements (CRs) are perceived to be related to sex chromosome evolution, but it is a matter of controversy whether CRs are the initial causative mechanism of suppressed recombination for sex differentiation. The early stages of sex chromosome evolution in amphibians may represent intermediate states of differentiation, and if so, they potentially shed light on the ultimate cause of suppressed recombination and the role of CRs in sex chromosome differentiation. In this paper, we showed that sex determination differs among 16 populations of spiny frog (Quasipaa boulengeri), in which individuals have normal and rearranged chromosomes caused by reciprocal translocation. In eastern areas, without translocation, genetic differentiation between sexes was relatively low, suggesting unrestricted recombination. In comparison, in western populations that have both normal and translocated chromosomes, a male-heterogametic system and lack of X-Y recombination were identified by male-specific alleles and heterozygote excess. However, such genetic differentiation between sexes in western populations was not directly related to karyotypes, as it was found in individuals with both normal and translocated karyotypes. In the western Sichuan Basin, male-specific and translocation-specific allelic frequency distributions suggested that recombination of sex-differentiation ceased in all populations, but recombination suppression caused by translocation did not exist in some populations. Combined with phylogenetic inference, this indicated that the establishment of sex-linkage had taken place independently of reciprocal translocation, and translocation was not the ultimate cause of sex chromosome differentiation. Furthermore, comparison of the genetic diversity of alleles on Y chromosomes, X chromosomes, and autosomes in western populations showed a reduction of effective population size on sex chromosomes, which may be caused by reciprocal translocation. It indicates that, although it is not the ultimate cause of recombination suppression, reciprocal translocation may enhance sex chromosome differentiation.

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.

Complete mitochondrial genomes of Nanorana taihangnica and N. yunnanensis (Anura: Dicroglossidae) with novel gene arrangements and phylogenetic relationship of Dicroglossidae

BACKGROUND

Complete mitochondrial (mt) genomes have been used extensively to test hypotheses about microevolution and to study population structure, phylogeography, and phylogenetic relationships of Anura at various taxonomic levels. Large-scale mt genomic reorganizations have been observed among many fork-tongued frogs (family Dicroglossidae). The relationships among Dicroglossidae and validation of the genus Feirana are still problematic. Hence, we sequenced the complete mt genomes of Nanorana taihangnica (=F. taihangnica) and N. yunnanensis as well as partial mt genomes of six Quasipaa species (dicroglossid taxa), two Odorrana and two Amolops species (Ranidae), and one Rhacophorus species (Rhacophoridae) in order to identify unknown mt gene rearrangements, to investigate the validity of the genus Feirana, and to test the phylogenetic relationship of Dicroglossidae.

RESULTS

In the mt genome of N. taihangnica two trnM genes, two trnP genes and two control regions were found. In addition, the trnA, trnN, trnC, and trnQ genes were translocated from their typical positions. In the mt genome of N. yunnanensis, three control regions were found and eight genes (ND6, trnP, trnQ, trnA, trnN, trnC, trnY and trnS genes) in the L-stand were translocated from their typical position and grouped together. We also found intraspecific rearrangement of the mitochondrial genomes in N. taihangnica and Quasipaa boulengeri. In phylogenetic trees, the genus Feirana nested deeply within the clade of genus Nanorana, indicating that the genus Feirana may be a synonym to Nanorana. Ranidae as a sister clade to Dicroglossidae and the clade of (Ranidae + Dicroglossidae) as a sister clade to (Mantellidae + Rhacophoridae) were well supported in BI analysis but low bootstrap in ML analysis.

CONCLUSIONS

We found that the gene arrangements of N. taihangnica and N. yunnanensis differed from other published dicroglossid mt genomes. The gene arrangements in N. taihangnica and N. yunnanensis could be explained by the Tandem Duplication and Random Loss (TDRL) and the Dimer-Mitogenome and Non-Random Loss (DMNR) models, respectively. The invalidation of the genus Feirana is supported in this study.

Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders

Background

Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication.

Results

We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years.

Conclusions

The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4358-2) contains supplementary material, which is available to authorized users.

Comparative mitochondrial genome analysis of Spilarctia subcarnea and other noctuid insects

This study was performed to better understand the phylogenetic relationships within the lepidopteran superfamily Noctuoidea. The mitochondrial genome (mitogenome) has been extensively used for studying phylogenetic relationships at different taxonomic levels. In this study, the complete mitogenome of Spilarctia subcarnea (Noctuoidea: Erebidae) was sequenced and annotated. The mitogenome is 15,441bp in length, containing 13 typical protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a noncoding control region (CR). The order and orientation of genes of S. subcarnea mitogenome with the order trnM-trnI-trnQ-nad2 is different from the ancestral insects in which trnM is located between trnQ and nad2 (trnI-trnQ-trnM-nad2). The phylogenetic relationships based on mitochondrial sequences using Bayesian inference and Maximum likelihood methods showed that S. subcarnea was closely related to Lemyra melli, supporting that S. subcarnea belongs to Erebidae. These analyses confirm that Lymantriidae should be included as subfamilies within Erebidae. The Erebidae was sister to (Nolidae+(Euteliidae+Noctuidae)); Notodontidae is sister to the other families of Noctuoidea in our study.

Next-generation sequencing of mixed genomic DNA allows efficient assembly of rearranged mitochondrial genomes in Amolops chunganensis and Quasipaa boulengeri

Recent improvements in next-generation sequencing (NGS) technologies can facilitate the obtainment of mitochondrial genomes. However, it is not clear whether NGS could be effectively used to reconstruct the mitogenome with high gene rearrangement. These high rearrangements would cause amplification failure, and/or assembly and alignment errors. Here, we choose two frogs with rearranged gene order, Amolops chunganensis and Quasipaa boulengeri, to test whether gene rearrangements affect the mitogenome assembly and alignment by using NGS. The mitogenomes with gene rearrangements are sequenced through Illumina MiSeq genomic sequencing and assembled effectively by Trinity v2.1.0 and SOAPdenovo2. Gene order and contents in the mitogenome of A. chunganensis and Q. boulengeri are typical neobatrachian pattern except for rearrangements at the position of "WANCY" tRNA genes cluster. Further, the mitogenome of Q. boulengeri is characterized with a tandem duplication of trnM. Moreover, we utilize 13 protein-coding genes of A. chunganensis, Q. boulengeri and other neobatrachians to reconstruct the phylogenetic tree for evaluating mitochondrial sequence authenticity of A. chunganensis and Q. boulengeri. In this work, we provide nearly complete mitochondrial genomes of A. chunganensis and Q. boulengeri.