Comparative mitogenomic analysis of mirid bugs (Hemiptera: Miridae) and evaluation of potential DNA barcoding markers

Complete mitochondrial genome of Stethoconus japonicus (Hemiptera: Miridae): Insights into the evolutionary traits within the family Miridae

The mitochondrial (mt) genome sequence of insects possesses numerous evolutionary traits. To better understand the evolution of mt genomes within the family Miridae, the complete mt genome of the predatory Japanese plant bug Stethoconus japonicus Schumacher was sequenced before undertaking a comparative analysis of all reported plant bug mt genomes. The mt genome of S. japonicus is a closed-circular and double-stranded DNA molecule of 16,274 bp (GenBank: MK341530), which consists of 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and a putative control region (CR). Consistent with other plant bugs, the mt genome of S. japonicus is strongly AT-biased (73.49 %) with A-skewed (0.202) and C-skewed (-0.248). All 13 PCGs initiate translation using ATN codons and TAA served as complete stop codons for eight PCGs, which as incomplete stop codon "T-" for cox1, nad1, nad5-6 and "TA-" for cox2. Regarding other features, all 22 tRNAs could be folded into typical cloverleaf secondary structures. The control region is 1,717 bp and contains a long tandem repeat sequence of a 165 bp unit repeated six times. Similar sequence with variable number of tandemly repeated units from intra-genus CRs is a distinct characteristic of plant bug mt genomes. Phylogenetic relationships of 15 bugs were eventually analyzed based on Maximum likelihood (ML) and Bayesian inference (BI) methods using 17 mt genome sequences. In the phylogenetic trees, species from a same genus or subfamily are clustered into a branch with high supporting values.and the result suggest that Deraeocorinae is more closely related to Mirinae than Bryocorinae. Finally, this study revealed that mutation of tRNA anticodon is a useful phylogenetic marker that could be of significance for studies of evolutionary patterns.

Molecular identification and phylogenetic confirmation of Sarcocystis species in slaughtered camels in Al-Najif province, Iraq

Background

Sarcocystis is an intracellular parasite of particular importance as it infects many domestic animals as camels that play the role of intermediate host for the parasite.

Aim

This study aimed to identify Sarcocystis species in camels by molecular assay with confirmation of local isolates by phylogenetic analysis.

Methods

A total of 200 slaughtered camels (Camelus dromedarius) that were slaughtered in Al-Najaf province (Iraq) abattoirs from October (2021) to July (2022) were subjected to collect the fresh tissues from four organs (esophagus, diaphragm, skeletal muscle, and heart), to be tested later by the conventional polymerase chain reaction (PCR). Then, a total of 20 positive genomic DNA samples were sequenced, named, got specific access numbers (OP785703.1 to OP785722.1), and compared with the NCBI-GenBank isolates.

Results

Targeting Cox1 gene, 80% of collected tissues were found positive by the conventional PCR assay. Phylogenetic tree analysis revealed that the local Sarcocystis isolates were identical to Indian S. cameli isolates at 99.70%-99.90%. Significantly, an increase in Sarcocystis infection was seen in the esophagus compared to the diaphragm, skeletal muscle, and heart; older (>4 years) than younger (≤4 years) camels, and in females more than males.

Conclusion

To the best of our knowledge, this represents the first molecular study in Iraq that identifies Sarcocystis cameli in camels. However, additional epidemiological and molecular studies in camel populations as well as in other domestic and wild animals appeared to be necessary.

Comparative mitogenomics and phylogenetic analyses of the genus Menida (Hemiptera, Heteroptera, Pentatomidae)

In order to explore the genetic diversity and phylogenetic relationship of the genus Menida Motschulsky, 1861 and reveal the molecular evolution of the family Pentatomidae, subfamily Pentatominae, complete mitochondrial genomes of three species of Menida were sequenced, and the phylogenetic relationships of tribes within the subfamily Pentatominae were studied based on these results. The mitochondrial genomes of Menidamusiva (Jakovlev, 1876), M.lata Yang, 1934, and M.metallica Hsiao & Cheng, 1977 were 16,663 bp, 16,463 bp, and 16,418 bp, respectively, encoding 37 genes and including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and a control region. The mitochondrial genome characteristics of Menida were compared and analyzed, and the phylogenetic tree of the Pentatominae was constructed based on the mitochondrial genome datasets using Bayesian inference (BI) and maximum likelihood (MI) methods. The results showed that gene arrangements, nucleotide composition, codon preference, gene overlaps, and RNA secondary structures were highly conserved within the Menida and had more similar characteristics in Pentatominae. The phylogenetic analysis shows a highly consistent topological structure based on BI and ML methods, which supported that the genus Menida belongs to the Pentatominae and is closely related to Hoplistoderini. The examined East Asian species of Menida form a monophyletic group with the internal relationships: (M.musiva + (M.lata + (M.violacea + M.metallica))). In addition, these results support the monophyly of Eysarcorini and Strachiini. Placosternum and Cappaeini are stable sister groups in the evolutionary branch of Pentatominae. The results of this study enrich the mitochondrial genome databases of Pentatominae and have significance for further elucidation of the phylogenetic relationships within the Pentatominae.

Mitogenomic Features and Evolution of the Nile River Dominant Tilapiine Species (Perciformes: Cichlidae)

To better understand the diversity and evolution of cichlids, we sequenced, assembled, and annotated the complete mitochondrial genomes of three Nile tilapiine species (Coptodon zillii, Oreochromis niloticus, and Sarotherodon galilaeus) dominating the Nile River waters. Our results showed that the general mitogenomic features were conserved among the Nile tilapiine species. The genome length ranged from 16,436 to 16,631 bp and a total of 37 genes were identified (two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and 1 control region). The ND6 was the only CDS that presented a negative AT skew and a positive GC skew. The most extended repeat sequences were in the D-loop followed by the pseudogenes (trnS^GCU). The ND5 showed relatively high substitution rates whereas ATP8 had the lowest substitution rate. The codon usage bias displayed a greater quantity of NNA and NNC at the third position and anti-bias against NNG. The phylogenetic relationship based on the complete mitogenomes and CDS was able to differentiate the three species as previously reported. This study provides new insight into the evolutionary connections between various subfamilies within cichlids while providing new molecular data that can be applied to discriminate between Nile tilapiine species and their populations.

Molecular and Developmental Signatures of Genital Size Macro-Evolution in Bugs

Our understanding of the genetic architecture of phenotypic traits has experienced drastic growth over the last years. Nevertheless, the majority of studies associating genotypes and phenotypes have been conducted at the ontogenetic level. Thus, we still have an elusive knowledge of how these genetic-developmental architectures evolve themselves and how their evolution is mirrored in the phenotypic change across evolutionary time. We tackle this gap by reconstructing the evolution of male genital size, one of the most complex traits in insects, together with its underlying genetic architecture. Using the order Hemiptera as a model, spanning over 350 million years of evolution, we estimate the correlation between genitalia and three features: development rate, body size, and rates of DNA substitution in 68 genes associated with genital development. We demonstrate that genital size macro-evolution has been largely dependent on body size and weakly influenced by development rate and phylogenetic history. We further revealed significant correlations between mutation rates and genital size for 19 genes. Interestingly, these genes have diverse functions and participate in distinct signaling pathways, suggesting that genital size is a complex trait whose fast evolution has been enabled by molecular changes associated with diverse morphogenetic processes. Our data further demonstrate that the majority of DNA evolution correlated with the genitalia has been shaped by negative selection or neutral evolution. Thus, in terms of sequence evolution, changes in genital size are predominantly facilitated by relaxation of constraints rather than positive selection, possibly due to the high pleiotropic nature of the morphogenetic genes.

Comparative mitogenomic analyses provide evolutionary insights into the retrolateral tibial apophysis clade (Araneae: Entelegynae)

The retrolateral tibial apophysis (RTA) clade is the largest spider lineage within Araneae. To better understand the diversity and evolution, we newly determined mitogenomes of ten RTA species from six families and performed a comparative mitogenomics analysis by combining them with 40 sequenced RTA mitogenomes available on GenBank. The ten mitogenomes encoded 37 typical mitochondrial genes and included a large non-coding region (putative control region). Nucleotide composition and codon usage were well conserved within the RTA clade, whereas diversity in sequence length and structural features was observed in control region. A reversal of strand asymmetry in nucleotide composition, i.e., negative AT-skews and positive GC-skews, was observed in each RTA species, likely resulting from mitochondrial gene rearrangements. All protein-coding genes were evolving under purifying selection, except for atp8 whose Ka/Ks was larger than 1, possibly due to positive selection or selection relaxation. Both mutation pressure and natural selection might contribute to codon usage bias of 13 protein-coding genes in the RTA lineage. Phylogenetic analyses based on mitogenomic data recovered a family-level phylogeny within the RTA; {[(Oval calamistrum clade, Dionycha), Marronoid clade], Sparassidae}. This study characterized RTA mitogenomes and provided some new insights into the phylogeny and evolution of the RTA clade.

Intraspecific comparison of mitochondrial genomes reveals the evolution in medicinal fungus Ganoderma lingzhi

Several mitogenomes of the genus Ganoderma have been assembled, but intraspecific comparisons of mitogenomes in Ganoderma lingzhi have not been reported. In this study, 19 G. lingzhi mitogenomes were assembled and analyzed combined with three mitogenomes of G. lingzhi from GenBank in term of the characteristics, evolution, and phylogeny. The results showed that the mitogenomes of the G. lingzhi strains are closed circular ranging from 49.23 kb to 68.37 kb. The genetic distance, selective pressure, and base variation indicate that the 14 common protein coding genes were highly conserved. The differences in introns, open reading frames, and repetitive sequences in the mitogenome were the main factors leaded to the variations in mitogenome. The introns were horizontally transferred in mitogenomes, and the differences between introns in the same insertion, which were primarily caused by the repetitive sequence, showed that the introns may be under degeneration. Besides, the frequent insertion and deletion of introns showed an evolutionary rate faster than protein coding genes. Phylogenetic analysis showed that the G. lingzhi strains gathered with high support, and those with the same intron distribution law had closer clustering relationships.

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.

The mitochondrial genomes of Tortricidae: nucleotide composition, gene variation and phylogenetic performance

BACKGROUND

Mitochondrial genomes (mitogenomes) have greatly improved our understanding of the backbone phylogeny of Lepidoptera, but few studies on comparative mitogenomics below the family level have been conducted. Here, we generated 13 mitogenomes of eight tortricid species, reannotated 27 previously reported mitogenomes, and systematically performed a comparative analysis of nucleotide composition, gene variation and phylogenetic performance.

RESULTS

The lengths of completely sequenced mitogenomes ranged from 15,440 bp to 15,778 bp, and the gene content and organization were conserved in Tortricidae and typical for Lepidoptera. Analyses of AT-skew and GC-skew, the effective number of codons and the codon bias index all show a base bias in Tortricidae, with little heterogeneity among the major tortricid groups. Variations in the divergence rates among 13 protein-coding genes of the same tortricid subgroup and of the same PCG among tortricid subgroups were detected. The secondary structures of 22 transfer RNA genes and two ribosomal RNA genes were predicted and comparatively illustrated, showing evolutionary heterogeneity among different RNAs or different regions of the same RNA. The phylogenetic uncertainty of Enarmoniini in Tortricidae was confirmed. The synonymy of Bactrini and Olethreutini was confirmed for the first time, with the representative Bactrini consistently nesting in the Olethreutini clade. Nad6 exhibits the highest phylogenetic informativeness from the root to the tip of the resulting tree, and the combination of the third coding positions of 13 protein-coding genes shows extremely high phylogenetic informativeness.

CONCLUSIONS

This study presents 13 mitogenomes of eight tortricid species and represents the first detailed comparative mitogenomics study of Tortricidae. The results further our understanding of the evolutionary architectures of tortricid mitogenomes and provide a basis for future studies of population genetics and phylogenetic investigations in this group.

DNA barcoding and nutritional analysis as a tool for promoting the market of inland fish species

The increasing world market demand for seafood requires an expansion of product categories available to consumers. Inland fish are usually considered having unmarked taste and are less appreciated by consumers; thus, they have low commercial value. Therefore, the marketing of the lake's fresh and processed fish is limited to the local market and consumers are currently uninformed and mistrustful about these species. In this study, six different fish species were caught in the Fondi lake (Lazio, central Italy): Anguilla anguilla, Tinca tinca, Carassius gibelio, Cyprinus carpio, Micropterus salmoides, Chelon ramada. All the samples were subjected to nutritional and DNA barcoding analysis. Moisture, protein, fat, carbohydrates, ash, and sodium content were measured. As regards the fatty acids profile, the most abundant were MUFAs with the highest value in Anguilla anguilla (45.97%). Oleic acid (C18: 1 n9 cis) was particularly high in Cyprinus carpio (55.46%). The fraction of polyunsaturated fatty acids (PUFA) revealed a higher DHA content (C22: 6 n3) in Anguilla anguilla than the other species (>12 %) while Chelon ramada presented both higher EPA content (C 20: 5 n3) and total fraction of omega 3 PUFAs. Concerning molecular analysis, a 655 bp fragment of cytochrome C oxidase subunit I (COI) gene was successfully used for the identification at the species level using both BOLD and BLAST public databases. The present study gives the basis for improving the knowledge and promoting inland fish' market and traceability along the supply chain.

Five Mitochondrial Genomes of the Genus Eysarcoris Hahn, 1834 with Phylogenetic Implications for the Pentatominae (Hemiptera: Pentatomidae)

Simple Summary

Pentatominae is the largest subfamily in the Pentatomidae, and most of its species are considered important agricultural pests. The phylogenetic relationships of tribes within Pentatominae remain controversial despite the fact that many studies have been performed using various molecular markers. In this study, five mitogenomes of the genus Eysarcoris were sequenced and analyzed, and the phylogenetic relationships of tribes within Pentatominae were reconstructed. The gene arrangement of the five mitochondrial genomes were found to be conserved and identical to other heteropteran mitogenomes. Differences in start codon usage and tandem repeats within control regions were found between E. gibbosus and the other four Eysarcoris species. In addition, the phylogenetic analyses showed that E. gibbosus is the first diverging clade within Eysarcorini. The results support the proposal to transfer E. gibbosus to the Stagonomus, another genus of Eysarcorini. Our results clarified the phylogenetic relationships among tribes of Pentatominae and laid a foundation for the further studies of Pentatominae.

Abstract

Four complete mitogenomes of Eysarcoris rosaceus, E. montivagus, E. gibbosus, E. annamita and one near-complete mitochondrial genome of E. ventralis were sequenced and used to explore the phylogenetic relationships of tribes within the subfamily Pentatominae. The mitogenomes range from 15,422 to 16,043 base pairs (bp) in length and encode 37 genes, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes (21 in E. ventralis), and a control region. Similar to other heteropteran species, the AT contents of the sequenced species were higher than their GC contents. The most frequently used start/stop codon was ATN/TAA. GTG was only found in atp6 and atp8 of E. gibbosus. All transfer RNA genes (tRNAs) exhibit the typical cloverleaf secondary structure, except for the trnS1 and trnV, which lacks the stem of the DHU arm. The length and copy number of repeat units were conserved within Eysarcoris, with the exception of E. gibbosus. Phylogenetic analyses based on mitogenomes using both maximum likelihood (ML) and Bayesian inference (BI) methods strongly supported the relationship among tribes within Pentatominae and confirmed that Graphosoma should be an intermediate lineage of Pentatominae. The relationship between Eysarcoris and Carbula was strongly supported and combined with our previous geometric morphometrics and chromosomal studies, suggest the Eysarcoris should belong to the tribe Eyasrcorini. This work will help to enhance our understanding of mitochondrial genomic evolution and phylogenetic relationships in Pentatominae.

Characterization of the Complete Mitochondrial Genome of Basidiomycete Yeast Hannaella oryzae: Intron Evolution, Gene Rearrangement, and Its Phylogeny

In this study, the mitogenome of Hannaella oryzae was sequenced by next-generation sequencing (NGS) and successfully assembled. The H. oryzae mitogenome comprised circular DNA molecules with a total size of 26,444 bp. We found that the mitogenome of H. oryzae partially deleted the tRNA gene transferring cysteine. Comparative mitogenomic analyses showed that intronic regions were the main factors contributing to the size variations of mitogenomes in Tremellales. Introns of the cox1 gene in Tremellales species were found to have undergone intron loss/gain events, and introns of the H. oryzae cox1 gene may have different origins. Gene arrangement analysis revealed that H. oryzae contained a unique gene order different from other Tremellales species. Phylogenetic analysis based on a combined mitochondrial gene set resulted in identical and well-supported topologies, wherein H. oryzae was closely related to Tremella fuciformis. This study represents the first report of mitogenome for the Hannaella genus, which will allow further study of the population genetics, taxonomy, and evolutionary biology of this important phylloplane yeast and other related species.

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.

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 first eleven mitochondrial genomes from the ectomycorrhizal fungal genus (Boletus) reveal intron loss and gene rearrangement

In the present study, eleven novel complete mitogenomes of Boletus were assembled and compared. The eleven complete mitogenomes were all composed of circular DNA molecules, with sizes ranging from 32,883 bp to 48,298 bp. The mitochondrial gene arrangement of Boletus varied greatly from other Boletales mitogenomes, and gene position reversal were observed frequently in the evolution of Boletus. Across the 15 core protein-coding genes (PCGs) tested, atp9 had the least and rps3 had the largest genetic distances among the eleven Boletus species, indicating varied evolution rates of core PCGs. In addition, the Ka/Ks value for nad3 gene was >1, suggesting that this gene was subject to possible positive selection pressure. Comparative mitogenomic analysis indicated that the intronic region was significantly correlated with the size of mitogenomes in Boletales. Two large-scale intron loss events were detected in the evolution of Boletus. Phylogenetic analyses based on a combined mitochondrial gene dataset yielded a well-supported (BPP ≥ 0.99; BS =100) phylogenetic tree for 72 Agaricomycetes, and the Boletus species had a close relationship with Paxillus. This study served as the first report on complete mitogenomes in Boletus, which will further promote investigations of the genetics, evolution and phylogeny of the Boletus genus.

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.

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

The complete mitochondrial genome of medicinal fungus Taiwanofungus camphoratus reveals gene rearrangements and intron dynamics of Polyporales

Taiwanofungus camphoratus is a highly valued medicinal mushroom that is endemic to Taiwan, China. In the present study, the mitogenome of T. camphoratus was assembled and compared with other published Polyporales mitogenomes. The T. camphoratus mitogenome was composed of circular DNA molecules, with a total size of 114,922 bp. Genome collinearity analysis revealed large-scale gene rearrangements between the mitogenomes of Polyporales, and T. camphoratus contained a unique gene order. The number and classes of introns were highly variable in 12 Polyporales species we examined, which proved that numerous intron loss or gain events occurred in the evolution of Polyporales. The Ka/Ks values for most core protein coding genes in Polyporales species were less than 1, indicating that these genes were subject to purifying selection. However, the rps3 gene was found under positive or relaxed selection between some Polyporales species. Phylogenetic analysis based on the combined mitochondrial gene set obtained a well-supported topology, and T. camphoratus was identified as a sister species to Laetiporus sulphureus. This study served as the first report on the mitogenome in the Taiwanofungus genus, which will provide a basis for understanding the phylogeny and evolution of this important fungus.

A Rapid Method for the Identification of Fresh and Processed Pagellus erythrinus Species against Frauds

The commercialization of porgies or seabreams of the family Sparidae has greatly increased in the last decade, and some valuable species have become subject to seafood substitution. DNA regions currently used for fish species identification in fresh and processed products belong to the mitochondrial (mt) genes cytochrome b (Cytb), cytochrome c oxidase I (COI), 16S and 12S. However, these markers amplify for fragments with lower divergence within and between some species, failing to provide informative barcodes. We adopted comparative mitogenomics, through the analysis of complete mtDNA sequences, as a compatible approach toward studying new barcoding markers. The intent is to develop a specific and rapid assay for the identification of the common pandora Pagellus erythrinus, a sparid species frequently subject to fraudulent replacement. The genetic diversity analysis (Hamming distance, p-genetic distance, gene-by-gene sequence variability) between 16 sparid mtDNA genomes highlighted the discriminating potential of a 291 bp NAD2 gene fragment. A pair of species-specific primers were successfully designed and tested by end-point and real-time PCR, achieving amplification only in P. erythrinus among several fish species. The use of the NAD2 barcoding marker provides a rapid presence/absence method for the identification of P. erythrinus.

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.

Comparative Mitogenome Analysis Reveals Mitochondrial Genome Differentiation in Ectomycorrhizal and Asymbiotic Amanita Species

In this present study, we assembled and analyzed the mitogenomes of two asymbiotic and six ectomycorrhizal Amanita species based on next-generation sequencing data. The size of the eight Amanita mitogenomes ranged from 37,341 to 137,428 bp, and we considered introns to be one of the main factors contributing to the size variation of Amanita. The introns of the cox1 gene experienced frequent gain/loss events in Amanita; and the intron position class cox1P386 was lost in the six ectomycorrhizal Amanita species. In addition, ectomycorrhizal Amanita species had more repetitive sequences and fewer intergenic sequences than asymbiotic Amanita species in their mitogenomes. Large-scale gene rearrangements were detected in the Amanita species we tested, including gene displacements and inversions. On the basis of the combined mitochondrial gene set, we reconstructed the phylogenetic relationships of 66 Basidiomycetes. The six ectomycorrhizal Amanita species were of single origin, and the two saprophytic Amanita species formed two distinct clades. This study is the first to elucidate the functions of the mitogenome in the evolution and ecological adaptation of Amanita species.

MtOrt: an empirical mitochondrial amino acid substitution model for evolutionary studies of Orthoptera insects

BACKGROUND

Amino acid substitution models play an important role in inferring phylogenies from proteins. Although different amino acid substitution models have been proposed, only a few were estimated from mitochondrial protein sequences for specific taxa such as the mtArt model for Arthropoda. The increasing of mitochondrial genome data from broad Orthoptera taxa provides an opportunity to estimate the Orthoptera-specific mitochondrial amino acid empirical model.

RESULTS

We sequenced complete mitochondrial genomes of 54 Orthoptera species, and estimated an amino acid substitution model (named mtOrt) by maximum likelihood method based on the 283 complete mitochondrial genomes available currently. The results indicated that there are obvious differences between mtOrt and the existing models, and the new model can better fit the Orthoptera mitochondrial protein datasets. Moreover, topologies of trees constructed using mtOrt and existing models are frequently different. MtOrt does indeed have an impact on likelihood improvement as well as tree topologies. The comparisons between the topologies of trees constructed using mtOrt and existing models show that the new model outperforms the existing models in inferring phylogenies from Orthoptera mitochondrial protein data.

CONCLUSIONS

The new mitochondrial amino acid substitution model of Orthoptera shows obvious differences from the existing models, and outperforms the existing models in inferring phylogenies from Orthoptera mitochondrial protein sequences.

Evolutionary rates of and selective constraints on the mitochondrial genomes of Orthoptera insects with different wing types

Orthoptera is the most diverse order of polyneopterans, and the forewing and hindwing of its members exhibit extremely variability from full length to complete loss in many groups; thus, this order provides a good model for studying the effects of insect flight ability on the evolutionary constraints on and evolutionary rate of the mitochondrial genome. Based on a data set of mitochondrial genomes from 171 species, including 43 newly determined, we reconstructed Orthoptera phylogenetic relationships and estimated the divergence times of this group. The results supported Caelifera and Ensifera as two monophyletic groups, and revealed that Orthoptera originated in the Carboniferous (298.997 Mya). The date of divergence between the suborders Caelifera and Ensifera was 255.705 Mya, in the late Permian. The major lineages of Acrididae seemed to have radiated in the Cenozoic, and the six patterns of rearrangement of 171 Orthoptera mitogenomes mostly occurred in the Cretaceous and Cenozoic. Based on phylogenetic relationships and ancestral state reconstruction, we analysed the evolutionary selection pressure on and evolutionary rate of mitochondrial protein-coding genes (mPCGs). The results indicated that during approximately 300 Mya of evolution, these genes experienced purifying selection to maintain their function. Flightless orthopteran insects accumulated more non-synonymous mutations than flying species and experienced more relaxed evolutionary constraints. The different wing types had different evolutionary rates, and the mean evolutionary rate of Orthoptera mitochondrial mPCGs was 13.554 × 10^-9 subs/s/y. The differences in selection pressures and evolutionary rates observed between the mitochondrial genomes suggested that functional constraints due to locomotion play an important role in the evolution of mitochondrial DNA in orthopteran insects with different wing types.

Mitogenome evolution in ladybirds: Potential association with dietary adaptation

Dietary shifts can alter the relative availability of different nutrients and are therefore associated with metabolic adaptation in animals. The Coccinellidae (ladybirds) exhibits three major types of feeding habits and provides a useful model to study the effects of dietary changes on the evolution of mitogenomes, which encode proteins directly involved in energy metabolism. Here, mitogenomes of three coccinellid species were newly sequenced. These data were combined with other ten previously sequenced coccinellid mitogenomes to explore the relationship between mitogenome evolution and diets. Our results indicate that mitogenomic data can be effectively used to resolve phylogenetic relationships of Coccinellidae. Strong codon usage bias in coccinellid mitogenomes was predominantly determined by nucleotide composition. The 13 mitochondrial protein-coding genes (PCGs) globally evolved under negative constraints, with some PCGs showing a stronger purifying selection. Six PCGs (nad3, nad4L, and nad5 from Complex I; cox1 and cox3 from Complex IV; and atp6 from Complex V) displayed signs of positive selection. Of these, adaptive changes in cox3 were potentially associated with metabolic differences resulting from dietary shifts in Coccinellidae. Our results provide insights into the adaptive evolution of coccinellid mitogenomes in response to both dietary shifts and other life history traits.

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.

Characterization of the mitochondrial genome of the pathogenic fungus Scytalidium auriculariicola (Leotiomycetes) and insights into its phylogenetics

Scytalidium auriculariicola is the causative pathogen of slippery scar disease in the cultivated cloud ear fungus, Auricularia polytricha. In the present study, the mitogenome of S. auriculariicola was sequenced and assembled by next-generation sequencing technology. The circular mitogenome is 96,857 bp long and contains 56 protein-coding genes, 2 ribosomal RNA genes, and 30 transfer RNA genes (tRNAs). The high frequency of A and T used in codons contributed to the high AT content (73.70%) of the S. auriculariicola mitogenome. Comparative analysis indicated that the base composition and the number of introns and protein-coding genes in the S. auriculariicola mitogenome varied from that of other Leotiomycetes mitogenomes, including a uniquely positive AT skew. Five distinct groups were found in the gene arrangements of Leotiomycetes. Phylogenetic analyses based on combined gene datasets (15 protein-coding genes) yielded well-supported (BPP = 1) topologies. A single-gene phylogenetic tree indicated that the nad4 gene may be useful as a molecular marker to analyze the phylogenetic relationships of Leotiomycetes species. This study is the first report on the mitochondrial genome of the genus Scytalidium, and it will contribute to our understanding of the population genetics and evolution of S. auriculariicola and related species.

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.

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 first complete mitochondrial genome from the family Hygrophoraceae (Hygrophorus russula) by next-generation sequencing and phylogenetic implications

Hygrophorus russula (Schaeff.) Kauffman is an edible ectomycorrhizal fungus that is widely distributed in the world. In this study, the mitogenome of H. russula was sequenced and assembled. The mitogenome of H. russula is composed of circular DNA molecules, with a total size of 55,769 bp. Further analysis indicated that the frequent use of A and T in codons contributes to the high AT content (80.87%) in the H. russula mitogenome. Comparative analysis indicated that the length and base composition of the core protein-encoding genes, and the number of tRNA genes in the H. russula mitogenome varied from that of other Agaricales mitogenomes. Gene arrangement analysis revealed a novel gene order in the H. russula mitogenome. In addition, the expansion of the mitogenome in Agaricales was found to be closely related to the increase in the number of introns. Phylogenetic analysis of the combined mitochondrial gene set showed strong support for tree topologies, and H. russula was determined to be relatively distant from other Agaricales species. This study is the first report on the mitogenome of a member of genus Hygrophorus as well as family Hygrophoraceae, which improves our understanding of mitochondrial differentiation and evolution in the important ectomycorrhizal fungi Hygrophorus species.

Characterization and comparative mitogenomic analysis of six newly sequenced mitochondrial genomes from ectomycorrhizal fungi (Russula) and phylogenetic analysis of the Agaricomycetes

In this study, the mitochondrial genomes of six Russula species were sequenced using next generation sequencing. The six mitogenomes were all composed of circular DNA molecules, with lengths ranging from 40,961 bp to 69,423 bp. The length and number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among the six mitogenomes. The increased number and total size of introns likely contributed to the size expansion of mitogenomes in some Russula species. Gene synteny analysis revealed some gene rearrangements among the six mitochondrial genomes. The nad4L gene had the lowest K2P genetic distance of the 15 core PCGs among the six Russula species, indicating that this gene was highly conserved. The Ka/Ks values for all 15 core PCGs were <1, suggesting that they were all subject to purifying selection. Phylogenetic analyses based on two gene datasets (15 core PCGs, and 15 core PCGs + rnl + rns) recovered identical and well-supported trees. In addition, cox1 was identified as a potential single-gene molecular marker for the phylogenetic analysis of relationships among Agaricomycetes species. This study provides the first report of mitogenomes from the Russulaceae family and facilitates the investigation of population genetics and evolution of other ectomycorrhizal fungi.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Complete Mitochondrial Genome of Dinorhynchus dybowskyi (Hemiptera: Pentatomidae: Asopinae) and Phylogenetic Analysis of Pentatomomorpha Species

Dinorhynchus dybowskyi (Hemiptera: Pentatomidae: Asopinae) is used as a biological control agent against various insect pests for its predatory. In the present study, the complete mitochondrial genome (mitogenome) of the species was sequenced using the next-generation sequencing technology. The results showed that the mitogenome is 15,952 bp long, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and a control region. Furthermore, the gene order and orientation of this mitogenome are identical to those of most heteropterans. There are 21 intergenic spacers (of length 1-28 bp) and 13 overlapping regions (of length 1-23 bp) throughout the genome. The control region is 1,291 bp long. The start codon of the PCGs is ATN, except cox1 (TTG), and stop codon is TAA, except nad1 (TAG). The 22 tRNAs exhibit a typical cloverleaf secondary structure, except trnS1, which lacks a dihydrouridine (DHU) arm and trnV, where the DHU arm forms a simple loop. The analyses based on nucleotide sequences of the 13 PCGs by Bayesian Inference and maximum likelihood methods. The results support the monophyly of five superfamilies Aradoidea, Pentatomoidea, Pyrrhocoroidea, Lygaeoidea, and Coreoidea. Within Pentatomoidea, the relationship observed is as follows: (Plataspidae + Urostylididae) + (Pentatomidae + (Acanthosomatidae + (Cydnidae + (Scutelleridae + (Dinidoridae + Tessaratomidae))))), and D. dybowskyi was placed in Pentatomidae and close to Eurydema gebleri.