Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains.

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Phylogenomics of endemic Australian Ulopinae (Hemiptera: Cicadomorpha: Cicadellidae).

Ulopinae is a distinctive subfamily of leafhoppers that is widely distributed across the Afrotropical, Palearctic, Indomalayan and Australasian regions. The ulopine fauna of Australia is entirely endemic and includes two tribes of striking appearance, the Ulopini and Cephalelini. Knowledge of these groups is fragmentary and in many instances, no information is available beyond original descriptions. We assess the monophyly, phylogenetic placement and species-level diversity of the Ulopini genus Austrolopa . Phylogenetic analyses based on sequence data from target nuclear loci (18S , 28S , H2A and H3 ) and mitochondrial genomes (15 genes) for 23 membracoid taxa yielded congruent topologies. Our results provide strong evidence for the monophyly of Ulopinae and a clade consisting of Ulopini + Cephalelini. However, a non-monophyletic Cephalelini arises from within a polyphyletic Ulopini. Austrolopa was strongly recovered as monophyletic in all analyses, a result also supported by morphological features. The genus currently includes six species, three of which are described based on morphological and molecular data: Austrolopa botanica , sp. nov. , Austrolopa rotunda , sp. nov. and Austrolopa sublima , sp. nov. A lectotype designation is provided for Austrolopa kingensis Evans, 1937, sp. reval. Our findings illustrate that the Australian Ulopinae is far more diverse than currently circumscribed and several species of Austrolopa are yet to be recognised. ZooBank: urn:lsid:zoobank.org:pub:1480285B-8F61-4659-A929-2B1EF3168868.

300 million years apart: the extreme case of macromorphological skeletal convergence between deltocyathids and a turbinoliid coral (Anthozoa, Scleractinia).

The integration of morphological and molecular lines of evidence has enabled the family Deltocyathidae to be erected to accommodate Deltocyathus species that were previously ascribed to the family Caryophylliidae. However, although displaying the same morphological characteristics as other species of Deltocyathus , molecular data suggested that D. magnificus was phylogenetically distant from Deltocyathidae, falling within the family Turbinoliidae instead. To elucidate the enigmatic evolutionary history of this species and skeletal microstructural features, the phylogenetic relationships of Deltocyathidae and Turbinoliidae were investigated using nuclear ultraconserved and exon loci and complete mitochondrial genomes. Both nuclear and mitochondrial phylogenomic reconstructions confirmed the position of D. magnificus within turbinolids. Furthermore, a novel mitochondrial gene order was uncovered for Deltocyathidae species. This gene order was not present in Turbinoliidae or in D. magnificus that both have the scleractinian canonical gene order, further indicating the taxonomic utility of mitochondrial gene order. D. magnificus is therefore formally moved to the family Turbinoliidae and accommodated in a new genus (Dennantotrochus Kitahara, Vaga & Stolarski, gen. nov.). Surprisingly, turbinolids and deltocyathids do not differ in microstructural organisation of the skeleton that consists of densely packed, individualised rapid accretion deposits and thickening deposits composed of fibres perpendicular to the skeleton surface. Therefore, although both families are clearly evolutionarily divergent, macromorphological features indicate a case of skeletal convergence while these may still share conservative biomineralisation mechanisms. ZooBank: urn:lsid:zoobank.org:pub:5F1C0E25-3CC6-4D1F-B1F0-CD9D0014678E.

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

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

Nuclear sequences of mitochondrial origin in domestic yak.

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the domestic yak genome. We obtained 499 alignment matches covering 340.2 kbp of the yak nuclear genome. After a merging step, we identified 167 NUMT regions with a total length of ~ 503 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic and mostly untranscribed. 98 different NUMT regions from domestic yak showed high homology with cow and/or wild yak genomes, suggesting selection or hybridization between domestic/wild yak and cow. To rule out the possibility that the identified NUMTs could be artifacts of the domestic yak genome assembly, we validated experimentally five NUMT regions by PCR amplification. As NUMT regions show high similarity to the mitochondrial genome can potentially pose a risk to domestic yak DNA mitochondrial studies, special care is therefore needed to select primers for PCR amplification of mitochondrial DNA sequences.

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes.

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing.

BACKGROUND: Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date. RESULTS: PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum. CONCLUSIONS: We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

A new long-read mitochondrial-genome protocol (PacBio HiFi) for haemosporidian parasites: a tool for population and biodiversity studies.

BACKGROUND: Studies on haemosporidian diversity, including origin of human malaria parasites, malaria's zoonotic dynamic, and regional biodiversity patterns, have used target gene approaches. However, current methods have a trade-off between scalability and data quality. Here, a long-read Next-Generation Sequencing protocol using PacBio HiFi is presented. The data processing is supported by a pipeline that uses machine-learning for analysing the reads. METHODS: A set of primers was designed to target approximately 6 kb, almost the entire length of the haemosporidian mitochondrial genome. Amplicons from different samples were multiplexed in an SMRTbell® library preparation. A pipeline (HmtG-PacBio Pipeline) to process the reads is also provided; it integrates multiple sequence alignments, a machine-learning algorithm that uses modified variational autoencoders, and a clustering method to identify the mitochondrial haplotypes/species in a sample. Although 192 specimens could be studied simultaneously, a pilot experiment with 15 specimens is presented, including in silico experiments where multiple data combinations were tested. RESULTS: The primers amplified various haemosporidian parasite genomes and yielded high-quality mt genome sequences. This new protocol allowed the detection and characterization of mixed infections and co-infections in the samples. The machine-learning approach converged into reproducible haplotypes with a low error rate, averaging 0.2% per read (minimum of 0.03% and maximum of 0.46%). The minimum recommended coverage per haplotype is 30X based on the detected error rates. The pipeline facilitates inspecting the data, including a local blast against a file of provided mitochondrial sequences that the researcher can customize. CONCLUSIONS: This is not a diagnostic approach but a high-throughput method to study haemosporidian sequence assemblages and perform genotyping by targeting the mitochondrial genome. Accordingly, the methodology allowed for examining specimens with multiple infections and co-infections of different haemosporidian parasites. The pipeline enables data quality assessment and comparison of the haplotypes obtained to those from previous studies. Although a single locus approach, whole mitochondrial data provide high-quality information to characterize species pools of haemosporidian parasites.

Complete mitochondrial genome of critically endangered catfish Hemibagrus punctatus (Jerdon, 1849) and comparative analysis for insights into the phylogeny of hemibagrids through mitogenomic approach.

BACKGROUND: Hemibagrus punctatus (Jerdon, 1849) is a critically endangered bagrid catfish endemic to the Western Ghats of India, whose population is declining due to anthropogenic activities. The current study aims to compare the mitogenome of H. punctatus with that of other Bagrid catfishes and provide insights into their evolutionary relationships. METHODS AND RESULTS: Samples were collected from Hemmige Karnataka, India. In the present study, the mitogenome of H. punctatus was successfully assembled, and its phylogenetic relationships with other Bagridae species were studied. The total genomic DNA of samples was extracted following the phenol-chloroform isoamyl alcohol method. Samples were sequenced, and the Illumina paired-end reads were assembled to a contig length of 16,517 bp. The mitochondrial genome was annotated using MitoFish and MitoAnnotator (Iwasaki et al., 2013). A robust phylogenetic analysis employing NJ (Maximum composite likelihood) and ASAP methods supports the classification of H. punctatus within the Bagridae family, which validates the taxonomic status of this species. In conclusion, this research enriches our understanding of H. punctatus mitogenome, shedding light on its evolutionary dynamics within the Bagridae family and contributing to the broader knowledge of mitochondrial genes in the context of evolutionary biology. CONCLUSIONS: The study's findings contribute to a better understanding of the mitogenome of H. punctatus and provide insights into the evolutionary relationships within other Hemibagrids.

Mitogenome architecture supports the non-monophyly of the cosmopolitan parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) recovered by nuclear and mitochondrial phylogenomics.

Mitochondrial DNA gene organisation is an important source of phylogenetic information for various metazoan taxa at different evolutionary timescales, though this has not been broadly tested for all insect groups nor within a phylogenetic context. The cosmopolitan subfamily Doryctinae is a highly diverse group of braconid wasps mainly represented by ectoparasitoids of xylophagous beetle larvae. Previous molecular studies based on Sanger and genome-wide (ultraconserved elements, UCE; and mitochondrial genomes) sequence data have recovered a non-monophyletic Doryctinae, though the relationships involved have always been weakly supported. We characterised doryctine mitogenomes and conducted separate phylogenetic analyses based on mitogenome and UCE sequence data of ~100 representative doryctine genera to assess the monophyly and higher-level classification of the subfamily. We identified rearrangements of mitochondrial transfer RNAs (tRNAs) that support a non-monophyletic Doryctinae consisting of two separate non-related clades with strong geographic structure ('New World' and 'Old World' clades). This geographic structure was also consistently supported by the phylogenetic analyses preformed with mitogenome and UCE sequence data. These results highlight the utility of the mitogenome gene rearrangements as a potential source of phylogenetic information at different evolutionary timescales.

A G-quadruplex-binding platinum complex induces cancer mitochondrial dysfunction through dual-targeting mitochondrial and nuclear G4 enriched genome.

BACKGROUND: G-quadruplex DNA (G4) is a non-canonical structure forming in guanine-rich regions, which play a vital role in cancer biology and are now being acknowledged in both nuclear and mitochondrial (mt) genome. However, the impact of G4-based targeted therapy on both nuclear and mt genome, affecting mt function and its underlying mechanisms remain largely unexplored. METHODS: The mechanisms of action and therapeutic effects of a G4-binding platinum(II) complex, Pt-ttpy, on mitochondria were conducted through a comprehensive approaches with in vitro and in vivo models, including ICP-MS for platinum measurement, PCR-based genetic analysis, western blotting (WB), confocal microscope for mt morphology study, extracellular flux analyzer, JC1 and Annexin V apoptosis assay, flow cytometry and high content microscope screening with single-cell quantification of both ROS and mt specific ROS, as well as click-chemistry for IF study of mt translation. Decipher Pt-ttpy effects on nuclear-encoded mt related genes expression were undertaken via RNA-seq, Chip-seq and CUT-RUN assays. RESULTS: Pt-ttpy, shows a highest accumulation in the mitochondria of A2780 cancer cells as compared with two other platinum(II) complexes with no/weak G4-binding properties, Pt-tpy and cisplatin. Pt-ttpy induces mtDNA deletion, copy reduction and transcription inhibition, hindering mt protein translation. Functional analysis reveals potent mt dysfunction without reactive oxygen species (ROS) induction. Mechanistic study provided first evidence that most of mt ribosome genes are highly enriched in G4 structures in their promoter regions, notably, Pt-ttpy impairs most nuclear-encoded mt ribosome genes' transcription through dampening the recruiting of transcription initiation and elongation factors of NELFB and TAF1 to their promoter with G4-enriched sequences. In vivo studies show Pt-ttpy's efficient anti-tumor effects, disrupting mt genome function with fewer side effects than cisplatin. CONCLUSION: This study underscores Pt-ttpy as a G4-binding platinum(II) complex, effectively targeting cancer mitochondria through dual action on mt and nuclear G4-enriched genomes without inducing ROS, offering promise for safer and effective platinum-based G4-targeted cancer therapy.

Real-Time Monitoring on the Chinese Giant Salamander Using RPA-LFD.

The Chinese giant salamander (Andrias davidianus), listed as an endangered species under "secondary protection" in China, faces significant threats due to ecological deterioration and the expansion of human activity. Extensive field investigations are crucial to ascertain the current status in the wild and to implement effective habitat protection measures to safeguard this species and support its population development. Traditional survey methods often fall short due to the elusive nature of the A. davidianus, presenting challenges that are time-consuming and generally ineffective. To overcome these obstacles, this study developed a real-time monitoring method that uses environmental DNA (eDNA) coupled with recombinase polymerase amplification and lateral flow strip (RPA-LFD). We designed five sets of species-specific primers and probes based on mitochondrial genome sequence alignments of A. davidianus and its close relatives. Our results indicated that four of these primer/probe sets accurately identified A. davidianus, distinguishing it from other tested caudata species using both extracted DNA samples and water samples from a tank housing an individual. This method enables the specific detection of A. davidianus genomic DNA at concentrations as low as 0.1 ng/mL within 50 min, without requiring extensive laboratory equipment. Applied in a field survey across four sites in Huangshan City, Anhui Province, where A. davidianus is known to be distributed, the method successfully detected the species at three of the four sites. The development of these primer/probe sets offers a practical tool for field surveying and monitoring, facilitating efforts in population recovery and resource conservation for A. davidianus.

The genetic identity of the Vedda: A language isolate of South Asia.

Linguistic data from South Asia identified several language isolates in the subcontinent. The Vedda, an indigenous population of Sri Lanka, are the least studied amongst them. Therefore, to understand the initial peopling of Sri Lanka and the genetic affinity of the Vedda with other populations in Eurasia, we extensively studied the high-resolution autosomal and mitogenomes from the Vedda population of Sri Lanka. Our autosomal analyses suggest a close genetic link of Vedda with the tribal populations of India despite no evidence of close linguistic affinity, thus suggesting a deep genetic link of the Vedda with these populations. The mitogenomic analysis supports this association by pointing to an ancient link with Indian populations. We suggest that the Vedda population is a genetically drifted group with limited gene flow from neighbouring Sinhalese and Sri Lankan Tamil populations. Interestingly, the genetic ancestry sharing of Vedda neglects the isolation-by-distance model. Collectively, the demography of Sri Lanka is unique, where Sinhalese and Sri Lankan Tamil populations excessively admixed, whilst Vedda largely preserved their isolation and deep genetic association with India.

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries).

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Asikli Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Asikli Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Asikli Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals.

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae).

BACKGROUND: Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes. RESULTS: We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit ß (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera. CONCLUSIONS: The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Evolution of myxozoan mitochondrial genomes: insights from myxobolids.

BACKGROUND: Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms. RESULTS: Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene. CONCLUSIONS: Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

Characterisation of the mitochondrial genome and phylogenetic analysis of Toxocara apodemi (Nematoda: Ascarididae).

We first sequenced and characterised the complete mitochondrial genome of Toxocara apodeme, then studied the evolutionary relationship of the species within Toxocaridae. The complete mitochondrial genome was amplified using PCR with 14 specific primers. The mitogenome length was 14303 bp in size, including 12 PCGs (encoding 3,423 amino acids), 22 tRNAs, 2 rRNAs, and 2 NCRs, with 68.38% A+T contents. The mt genomes of T. apodemi had relatively compact structures with 11 intergenic spacers and 5 overlaps. Comparative analyses of the nucleotide sequences of complete mt genomes showed that T. apodemi had higher identities with T. canis than other congeners. A sliding window analysis of 12 PCGs among 5 Toxocara species indicated that nad4 had the highest sequence divergence, and cox1 was the least variable gene. Relative synonymous codon usage showed that UUG, ACU, CCU, CGU, and UCU most frequently occurred in the complete genomes of T. apodemi. The Ka/Ks ratio showed that all Toxocara mt genes were subject to purification selection. The largest genetic distance between T. apodemi and the other 4 congeneric species was found in nad2, and the smallest was found in cox2. Phylogenetic analyses based on the concatenated amino acid sequences of 12 PCGs demonstrated that T. apodemi formed a distinct branch and was always a sister taxon to other congeneric species. The present study determined the complete mt genome sequences of T. apodemi, which provide novel genetic markers for further studies of the taxonomy, population genetics, and systematics of the Toxocaridae nematodes.