Mammalian mitochondrial D-loop region structural analysis: identification of new conserved sequences and their functional and evolutionary implications

Mitochondrial and Nuclear DNA Variants in Amyotrophic Lateral Sclerosis: Enrichment in the Mitochondrial Control Region and Sirtuin Pathway Genes in Spinal Cord Tissue

Amyotrophic Lateral Sclerosis (ALS) is a progressive disease with prevalent mitochondrial dysfunctions affecting both upper and lower motor neurons in the motor cortex, brainstem, and spinal cord. Despite mitochondria having their own genome (mtDNA), in humans, most mitochondrial genes are encoded by the nuclear genome (nDNA). Our study aimed to simultaneously screen for nDNA and mtDNA genomes to assess for specific variant enrichment in ALS compared to control tissues. Here, we analysed whole exome (WES) and whole genome (WGS) sequencing data from spinal cord tissues, respectively, of 6 and 12 human donors. A total of 31,257 and 301,241 variants in nuclear-encoded mitochondrial genes were identified from WES and WGS, respectively, while mtDNA reads accounted for 73 and 332 variants. Despite technical differences, both datasets consistently revealed a specific enrichment of variants in the mitochondrial Control Region (CR) and in several of these genes directly associated with mitochondrial dynamics or with Sirtuin pathway genes within ALS tissues. Overall, our data support the hypothesis of a variant burden in specific genes, highlighting potential actionable targets for therapeutic interventions in ALS.

Comparative mitochondrial genome brings insights to slight variation in gene proportion and large intergenic spacer and phylogenetic relationship of mudskipper species

Fish mitochondrial genome have been largely studied worldwide for evolutionary and other genetic purposes and the structure and gene organization are commonly conservative. However, several studies have demonstrated that this scenario may present variations in some taxa, showing differentiation on the gene rearrangement. In this study, the complete mitogenome of terrestrial fish Boleophthalmus dussumieri was generated and compared with other species of the Exudercidae fishes. The newly complete mitogenome generated is circular and 16,685 bp of length, and it contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA genes (tRNAs), and one control region (CR), with high conservative structure, like other Mudskippers. Most of the PCG showed similar codon usage bias. The gene length was found to be different specially for the CR, 12S rRNA gene and ND5 gene in some taxon. All the Boleophthalmus species showed a gene duplication in the CR, except for B. dussumieri, and they presented a long intergenic spacer specially on the tRNA-Pro/ OH Tandem duplication/random loss (TDRL) and dimer-mitogenome and nonrandom loss (DMNL) are suitable to explain the mitogenome rearrangement observed in this study. The phylogenetic analysis well supported the monophyly of all mudskipper species and the analysis positioned the Periophthalmus clade as the most basal of the terrestrial fishes. This finding provides basis and brings insights for gene variation, gene rearrangements and replications showing evidence for variety of mitochondrial structure diversity within mudskippers.

Mitochondrial DNA mutations in HIV-exposed uninfected infants following the cessation of triple antiretroviral drugs

OBJECTIVES

Mitochondrial mutations in HIV-exposed uninfected (HEU) infants after cessation of ART are rarely studied. We analysed a group of HEU newborns born to mothers with late HIV diagnosis who received three doses of ART immediately after birth. We observed mitochondrial DNA (mtDNA) mutations at different times of withdrawal.

METHODS

The study was based on a clinical trial conducted from 2015 to 2020. Newborns of the intervention group who met the criteria for this study received triple antiretroviral drugs, zidovudine + lamivudine + nevirapine, within 2 h after the birth, as post-partum prophylaxis, and at 14 days were switched to zidovudine + lamivudine + lopinavir/ritonavir, which was continued until 6 weeks of age. From August to November 2019, blood samples from HEU infants were also collected after ceasing 12 months of ART, and analysed for mtDNA.

RESULTS

Our study found that mtDNA mutations remained prevalent in HEU infants a few years after three ARTs were stopped immediately after birth. Among them, D-loop, ND1 and CYTB are the first three mutated regions during different withdrawal periods. This pattern of mutations is similar to, but not exactly consistent with, HIV-infected children receiving standard ART.

CONCLUSIONS

Further studies are needed to determine the effects of these mutations on the development of HEU infants and whether stopping ART leads to the restoration of mitochondrial function.

Development of mitochondrial gene-editing strategies and their potential applications in mitochondrial hereditary diseases: a review

Mitochondrial DNA (mtDNA) is a critical genome contained within the mitochondria of eukaryotic cells, with many copies present in each mitochondrion. Mutations in mtDNA often are inherited and can lead to severe health problems, including various inherited diseases and premature aging. The lack of efficient repair mechanisms and the susceptibility of mtDNA to damage exacerbate the threat to human health. Heteroplasmy, the presence of different mtDNA genotypes within a single cell, increases the complexity of these diseases and requires an effective editing method for correction. Recently, gene-editing techniques, including programmable nucleases such as restriction endonuclease, zinc finger nuclease, transcription activator-like effector nuclease, clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeats-associated 9 and base editors, have provided new tools for editing mtDNA in mammalian cells. Base editors are particularly promising because of their high efficiency and precision in correcting mtDNA mutations. In this review, we discuss the application of these techniques in mitochondrial gene editing and their limitations. We also explore the potential of base editors for mtDNA modification and discuss the opportunities and challenges associated with their application in mitochondrial gene editing. In conclusion, this review highlights the advancements, limitations and opportunities in current mitochondrial gene-editing technologies and approaches. Our insights aim to stimulate the development of new editing strategies that can ultimately alleviate the adverse effects of mitochondrial hereditary diseases.

Comparative Mitogenomics Reveals Cryptic Species in Sillago ingenuua McKay, 1985 (Perciformes: Sillaginidae)

It is unreliable to identify marine fishes only by external morphological features. Species misidentification brings great challenges to fishery research, resource monitoring and ecomanagement. Sillago ingenuua is an important part of commercial marine fishes, and in which, the morphological differences between different groups are not obvious. Here, we compared different geographical groups of S. ingenuua which were collected from Xiamen, Dongshan, Keelung, Songkhla and Java. The results showed that all samples of S. ingenuua were similar in external morphological characteristics and the shape of the swim bladder, but there were two distinctive lineages which were flagged as cryptic species based on DNA barcoding. The comparative mitogenomic results showed that S. ingenuua A and S. ingenuua B were identical in structural organization and gene arrangement. Their nucleotide composition and codon usage were also similar. A phylogenetic analysis was performed based on 13 concatenated PCGs from eight Sillago species. The results showed that the genetic distance between S. ingenuua A and S. ingenuua B was large (D = 0.069), and this genetic distance was large enough to reveal that S. ingenuua A and S. ingenuua B might be different species.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

mtDNA CR Evidence Indicates High Genetic Diversity of Captive Forest Musk Deer in Shaanxi Province, China

Forest musk deer (Moschus berezovskii) are endangered ruminants whose adult males secrete musk. China has been breeding forest musk deer artificially since the 1950s in an effort to restore wild populations, with Shaanxi and Sichuan provinces as the two main sites for captive breeding. Genetic diversity is a significant indicator that determines the long-term viability and status of a population, particularly for species at risk of extinction. In this study, we analyzed the current genetic makeup of seven captive forest musk deer populations in the Shaanxi province, using the mitochondrial DNA (mtDNA) control region (CR) as the molecular marker. We sequenced 604 bp of mtDNA CR, with an average content of A+T higher than G+C. We observed 111 variable sites and 39 different haplotypes from 338 sequences. The nucleotide diversity (Pi) and haplotype diversity (Hd) were 0.02887 and 0.908, respectively. Genetic differentiation between these populations was not significant, and the populations might not have experienced rapid growth. By combining our sequences with previous ones, we identified 65 unique haplotypes with 26 rare haplotypes and estimated a total of 90 haplotypes in Shaanxi province captive populations. The Shaanxi province and Sichuan province obtained 88 haplotypes, the haplotypes from the two populations were mixed together, and the two populations showed moderate genetic differentiation. Our findings suggested that captive forest musk deer populations in the Shaanxi province had high genetic diversity, with a rich founder population of about 90 maternal lines. Additionally, managers could develop genetic management plans for forest musk deer based on the haplotype database. Overall, our study will provide insights and guidelines for the conservation of genetic diversity in captive forest musk deer populations in the Shaanxi province.

A new type of tandem repeats in Myotis petax (Chiroptera, Vespertilionidae) mitochondrial control region

BACKGROUND

Tandem repeats in mitochondrial DNA control region are known to different animal taxa, including bat species of the family Vespertilionidae. The long R1-repeats in the bat ETAS-domain are often presented in a variable copy number and may exhibit both inter-individual and intra-individual sequence diversity. The function of repeats in the control region is still unclear, but it has been shown that repetitive sequences in some animal groups (shrews, cats and sheep) may include parts of ETAS1 and ETAS2 conservative blocks of mitochondrial DNA.

METHODS AND RESULTS

Analysis of the control region sequences for 31 Myotis petax specimens allowed the identification of the inter-individual variability and clarification of the composition of the R1-repeats. The copy number of the R1-repeats varies from 4 to 7 in individuals. The specimens examined do not exhibit a size heteroplasmy previously described for Myotis species. The unusual short 30 bp R1-repeats have been detected in M. petax for the first time. The ten specimens from Amur Region and Primorsky Territory have one or two copies of these additional repeats.

CONCLUSIONS

It was determined that the R1-repeats in M. petax control region consist of parts of the ETAS1 and ETAS2 blocks. The origin of the additional repeats seems to be related to the 51 bp deletion in the central part of the R1-repeat unit and subsequent duplication. Comparison of repetitive sequences in the control region of closely-related Myotis species identified the occurrence of incomplete repeats also resulting from the short deletions, but distinct from additional repeats of M. petax.

The mitochondrial genomes of big-eared bats, Macrotus waterhousii and Macrotus californicus (Chiroptera: Phyllostomidae: Macrotinae)

In the species-rich family Phyllostomidae, the genus Macrotus ('big eared' bats) contains only two species; Macrotus waterhousii, distributed in western, central, and southern Mexico, Guatemala and some Caribbean Islands, and Macrotus californicus, distributed in the southwestern USA, and in the Baja California peninsula and the state of Sonora in Mexico. In this study, we sequenced and assembled the mitochondrial genome of Macrotus waterhousii and characterized in detail this genome and that of the congeneric M. californicus. Then, we examined the phylogenetic position of Macrotus in the family Phyllostomidae based on protein coding genes (PCGs). The AT-rich mitochondrial genomes of M. waterhousii and M. californicus are 16,792 and 16,691 bp long, respectively, and each encode 13 PCGs, 22 tRNA genes, 2 rRNA genes, and a putative non-coding control region 1,336 and 1,232 bp long, respectively. Mitochondrial synteny in Macrotus is identical to that reported before for all other cofamilial species. In the two studied species, all tRNAs exhibit a 'typical' cloverleaf secondary structure with the exception of trnS1, which lacks the D arm. A selective pressure analysis demonstrated that all PCGs are under purifying selection. The CR of the two species feature three domains previously reported in other mammals, including bats: extended terminal associated sequences (ETAS), central (CD), and conserved sequence block (CSB). A phylogenetic analysis based on the 13 mitochondrial PCGs demonstrated that Macrotus is monophyletic and the subfamily Macrotinae is a sister group of all remaining phyllostomids in our analysis, except Micronycterinae. The assembly and detailed analysis of these mitochondrial genomes represents a step further to continue improving the understanding of phylogenetic relationships within the species-rich family Phyllostomidae.

Comparative mitochondrial genomics and phylogenetics for species of the snakehead genus Channa Scopoli, 1777 (Perciformes: Channidae)

The family Channidae, members of which are commonly known as snakehead fish, includes 53 Channa species and three Parachanna species. In this study, we characterized mitochondrial genomes (mitogenomes) of five Channa species (C. andrao, C. bleheri, C. ornatipinnis, C. pulchra, and C. stewartii) for the first time. We compared the mitogenomes with the mitogenomes of 11 other Channidae fish. The newly sequenced mitogenomes were 16,714 - 16,895 bp in length and contained 37 typical genes [13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs) and 22 transfer RNA genes (tRNAs)]. Positive AT-skews and negative GC-skews were found in the mitogenomes of Channidae. Most PCGs started with the conventional start codon, ATN; however, the sequence of the stop codon was variable. There was no essential difference in relative synonymous codon usage (RSCU) among the Channidae mitogenomes. The fastest-evolving gene atp8 and slowest-evolving gene cox1 were identified using Ka/Ks and pairwise relative genetic distance (p-distance). The displacement loop (D-loop) regions showed great variability, which affected the size of Channa mitogenomes. One origin of replication on the light strand (O_L) region was found between trnN and trnC in the mitogenomes of Channidae. Phylogenetic analysis revealed three new sister pairs (C. andrao + C. bleheri, C. ornatipinnis + C. pulchra, and C. stewartii + C. gachua). Phylogenetic relationships established between the five Channa species based on mitogenomes were also supported by their morphological characteristics and geographical distribution. The novel information we obtained about these mitogenomes will contribute to elucidating the complex relationships among Channa species.

Comparative Mitogenome Analyses Uncover Mitogenome Features and Phylogenetic Implications of the Parrotfishes (Perciformes: Scaridae)

In order to investigate the molecular evolution of mitogenomes among the family Scaridae, the complete mitogenome sequences of twelve parrotfish species were determined and compared with those of seven other parrotfish species. The comparative analysis revealed that the general features and organization of the mitogenome were similar among the 19 parrotfish species. The base composition was similar among the parrotfishes, with the exception of the genus Calotomus, which exhibited an unusual negative AT skew in the whole mitogenome. The PCGs showed similar codon usage, and all of them underwent a strong purifying selection. The gene rearrangement typical of the parrotfishes was detected, with the tRNAMet inserted between the tRNAIle and tRNAGln, and the tRNAGln was followed by a putative tRNAMet pseudogene. The parrotfish mitogenomes displayed conserved gene overlaps and secondary structure in most tRNA genes, while the non-coding intergenic spacers varied among species. Phylogenetic analysis based on the thirteen PCGs and two rRNAs strongly supported the hypothesis that the parrotfishes could be subdivided into two clades with distinct ecological adaptations. The early divergence of the sea grass and coral reef clades occurred in the late Oligocene, probably related to the expansion of sea grass habitat. Later diversification within the coral reef clade could be dated back to the Miocene, likely associated with the geomorphology alternation since the closing of the Tethys Ocean. This work provided fundamental molecular data that will be useful for species identification, conservation, and further studies on the evolution of parrotfishes.

Features and evolution of control regions in leafroller moths (Lepidoptera: Tortricidae) inferred from mitochondrial genomes and phylogeny

The control region (CR) of the mitochondrial genome (mitogenome) represents a major noncoding fragment with several special structural features that are thought to be responsible for the initiation of mitogenome transcription and replication. However, few studies have revealed the evolutionary patterns of CR in the phylogenetic context. Here, we explain the characteristics and evolution of CR in Tortricidae, inferred from a mitogenome-based phylogeny. The first complete mitogenomes of the genera Meiligma and Matsumuraeses were sequenced. Both mitogenomes are double-stranded circular DNA molecules with lengths of 15,675 bp and 15,330 bp, respectively. Phylogenetic analyses derived from 13 protein-coding genes and two ribosomal RNAs showed that most tribes, including subfamilies Olethreutinae and Tortricinae, were recovered as monophyletic clades, similar to previous studies based on morphological or nuclear data. Moreover, comprehensive comparative analyses of the structural organization and role of tandem replications on the length variation and high AT content of CR sequences were conducted. The results reveal a significant positive correlation between the total length and AT content of tandem repeats and whole CR sequences in Tortricidae. The structural organization in CR sequences is diverse, even between closely related tribes, which demonstrates the plasticity of the mitochondrial DNA molecule in Tortricidae.

Mitochondrial Genomes Assembled from Non-Invasive eDNA Metagenomic Scat Samples in Critically Endangered Mammals

The abundance of many large-bodied vertebrates, both in marine and terrestrial environments, has declined substantially due to global and regional climate stressors that define the Anthropocene. The development of genetic tools that can serve to monitor population’s health non-intrusively and inform strategies for the recovery of these species is crucial. In this study, we formally evaluate whether whole mitochondrial genomes can be assembled from environmental DNA (eDNA) metagenomics scat samples. Mitogenomes of four different large vertebrates, the panda bear (Ailuropoda melanoleuca), the moon bear (Ursus thibetanus), the Java pangolin (Manis javanica), and the the North Atlantic right whale (Eubalaena glacialis) were assembled and circularized using the pipeline GetOrganelle with a coverage ranging from 12x to 480x in 14 out of 18 different eDNA samples. Partial mitochondrial genomes were retrieved from three other eDNA samples. The complete mitochondrial genomes of the studied species were AT-rich and comprised 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a putative D-loop/control region. Synteny observed in all assembled mitogenomes was identical to that reported for specimens of the same and other closely related species. This study demonstrates that it is possible to assemble accurate whole mitochondrial chromosomes from eDNA samples (scats) using forthright bench and bioinformatics workflows. The retrieval of mitochondrial genomes from eDNA samples represents a tool to support bioprospecting, bio-monitoring, and other non-intrusive conservation strategies in species considered ‘vulnerable’, ‘endangered’, and/or ‘critically endangered’ by the IUCN Red List of Threatened Species.

Changes in Mitochondrial Epigenome in Type 2 Diabetes Mellitus

Type 2 Diabetes Mellitus is a major chronic metabolic disorder in public health. Due to mitochondria's indispensable role in the body, its dysfunction has been implicated in the development and progression of multiple diseases, including Type 2 Diabetes mellitus. Thus, factors that can regulate mitochondrial function, like mtDNA methylation, are of significant interest in managing T2DM. In this paper, the overview of epigenetics and the mechanism of nuclear and mitochondrial DNA methylation were briefly discussed, followed by other mitochondrial epigenetics. Subsequently, the association between mtDNA methylation with T2DM and the challenges of mtDNA methylation studies were also reviewed. This review will aid in understanding the impact of mtDNA methylation on T2DM and future advancements in T2DM treatment.

Next-Generation Sequencing of the Complete Huaibei Grey Donkey Mitogenome and Mitogenomic Phylogeny of the Equidae Family

The Huaibei grey donkey (HGD) is an endangered species and a vital native breed in Anhui Province, China. However, its complete mitogenome, phylogeny, and maternal origin remain unclear. The objectives of this study were to detect the genetic diversity of the HGD and investigate its phylogenetic relationship with other breeds to inform conservation management. The complete mitogenome of the HGD was sequenced through next-generation sequencing, and the most variable region in the mitochondrial DNA displacement-loop (D-loop) was amplified via a polymerase chain reaction (PCR). Next, we used the median-joining network (MJN) to calculate the genetic relationships among populations and the neighbor-jointing method to build a phylogenetic tree and speculate as to its origin. The results showed that the mitogenome contains 22 tRNAs, 2 rRNAs, 13 PCGs, and 1 D-loop region. Analyzing the D-loop region of the HGDs, we identified 23 polymorphic sites and 11 haplotypes. The haplotype and nucleotide diversity were 0.87000 (Hd) and 0.02115 (Pi), respectively. The MJN analysis indicated that the HGD potentially has two maternal lineages, and phylogenetic analysis indicated that the Somali lineage could be the most probable domestication center for this breed. Therefore, our mitogenome analysis highlights the high genetic diversity of the HGD, which may have originated from the Somali wild ass, as opposed to the Asian wild ass. This study will provide a useful resource for HGD conservation and breeding.

Single-molecule Sequencing of an Animal Mitochondrial Genome Reveals Chloroplast-like Architecture and Repeat-mediated Recombination

Recent advances in long-read sequencing technology have allowed for single-molecule sequencing of entire mitochondrial genomes, opening the door for direct investigation of the mitochondrial genome architecture and recombination. We used PacBio sequencing to reassemble mitochondrial genomes from two species of New Zealand freshwater snails, Potamopyrgus antipodarum and Potamopyrgus estuarinus. These assemblies revealed a ∼1.7 kb structure within the mitochondrial genomes of both species that was previously undetected by an assembly of short reads and likely corresponding to a large noncoding region commonly present in the mitochondrial genomes. The overall architecture of these Potamopyrgus mitochondrial genomes is reminiscent of the chloroplast genomes of land plants, harboring a large single-copy (LSC) region and a small single-copy (SSC) region separated by a pair of inverted repeats (IRa and IRb). Individual sequencing reads that spanned across the Potamopyrgus IRa-SSC-IRb structure revealed the occurrence of a "flip-flop" recombination. We also detected evidence for two distinct IR haplotypes and recombination between them in wild-caught P. estuarinus, as well as extensive intermolecular recombination between single-nucleotide polymorphisms in the LSC region. The chloroplast-like architecture and repeat-mediated mitochondrial recombination we describe here raise fundamental questions regarding the origins and commonness of inverted repeats in cytoplasmic genomes and their role in mitochondrial genome evolution.

Comparison and Phylogenetic Analysis of Mitochondrial Genomes of Talpidae Animals

Talpidae is a model group for evolutionary studies due to their highly specialized morphologies and diverse lifestyles. Mitochondrial genomes are molecular markers commonly used in species evolution and phylogenetic studies. In this study, the complete mitochondrial genome sequence of Scaptochirus moschatus was obtained by Illumina NovaSeq sequencing. The complete mitochondrial genomes of 14 Talpidae species (including Scaptochirus moschatus obtained in the present study) and the cytochrome b (Cyt b) gene sequences of 48 Talpidae species were downloaded from the NCBI database for comparison and phylogenetic studies to analyze the phylogenetic relationships and to find the possible reasons of the niche differentiation and ecotype specialization of Talpidae animals. The results showed that the mitochondrial genome sequences of 14 species belonging to the family Talpidae were 16,528 to 16,962 bp, all containing 13 protein-coding genes, 22 tRNA, two rRNA, and a non-coding region (control region). The difference in the number of repetitive repeats in the control region is responsible for the difference in the length of Talpidae mitochondrial genome sequences. Combining the divergence time of Talpidae animals with the geological history, it is found that the niche differentiation and ecotype divergence of Talpidae is closely related to historically global climate changes. Semi-aquatic groups diverged in the early Oligocene (about 31.22 MYA), probably in response to the global climate transition from warm to cool. During the early Miocene (about 19.54 MYA), some species of Talpidae moved to underground habitats and formed fossorial groups that were adept at digging due to the effects of the glaciation. In the middle Miocene (about 16.23 MYA), some Talpidae animals returned to the ground and formed semi-fossorial shrew moles as global climate warming again.

Mitochondrial DNA Diversity of Mesocricetus auratus and Other Cricetinae Species among Cricetidae Family

Unique anatomical and physiological features have made hamster species desirable research models. Comparative genomics and phylogenetic analysis of the hamster family members to clarify their evolution and genetic relationship, can provide a genetic basis for the comprehension of the variable research results obtained using different hamster models. The Syrian golden hamster (Mesocricetus auratus) is the most widely used species. In this study, we sequenced the complete mitochondrial genome (mitogenome) of M. auratus, compared it with the mitogenome of other Cricetinae subfamily species, and defined its phylogenetic position in the Cricetidae family. Our results show that the mitogenome organization, gene arrangement, base composition, and genetic analysis of the protein coding genes (PCGs) of M. auratus are similar to those observed in previous reports on Cricetinae species. Nonetheless, our analysis clarifies some striking differences of M. auratus relative to other subfamily members, namely distinct codon usage frequency of TAT (Tyr), AAT (Asn), and GAA (Glu) and the presence of the conserved sequence block 3 (CSB-3) in the control region of M. auratus mitogenome and other hamsters (not found in Arvicolinae). These results suggest the particularity of amino acid codon usage bias of M. auratus and special regulatory signals for the heavy strand replication in Cricetinae. Additionally, Bayesian inference/maximum likelihood (BI/ML) tree shows that Cricetinae and Arvicolinae are sister taxa sharing a common ancestor, and Neotominae split prior to the split between Cricetinae and Arvicolinae. Our results support taxonomy revisions in Cricetulus kamensis and Cricetulus migratorius, and further revision is needed within the other two subfamilies. Among the hamster research models, Cricetulus griseus is the species with highest sequence similarity and closer genetic relationship with M. auratus. Our results show mitochondrial DNA diversity of M. auratus and other Cricetinae species and provide genetic basis for judgement of different hamster models, promoting the development and usage of hamsters with regional characteristics.

The complete mitochondrial genome of Talpa martinorum (Mammalia: Talpidae), a mole species endemic to Thrace: genome content and phylogenetic considerations

The complete mitogenome sequence of Talpa martinorum, a recently described Balkan endemic mole, was assembled from next generation sequence data. The mitogenome is similar to that of the three other Talpa species sequenced to date, being 16,835 bp in length, and containing 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, an origin of L-strand replication, and a control region or D-loop. Compared to other Talpa mitogenomes sequenced to date, that of T. martinorum differs in the length of D-loop and stop codon usage. TAG and T-- are the stop codons for the ND1 and ATP8 genes, respectively, in T. martinorum, whilst TAA acts as a stop codon for both ND1 and ATP8 in the other three Talpa species sequenced. Phylogeny reconstructions based on Maximum Likelihood and Bayesian inference analyses yielded phylogenies with similar topologies, demonstrating that T. martinorum nests within the western lineage of the genus, being closely related to T. aquitania and T. occidentalis.

Intraspecies variation of the mitochondrial genome: An evaluation for phylogenetic approaches based on the conventional choices of genes and segments on mitogenome

Intraspecies nucleotide sequence variation is a key to understanding the evolutionary history of a species, such as the geographic distribution and population structure. To date, numerous phylogenetic and population genetics studies have been conducted based on the sequences of a gene or an intergenic region on the mitochondrial genome (mtDNA), such as cytochrome c oxidase subunits or the D-loop. To evaluate the credibility of the usage of such ‘classic’ markers, we compared the phylogenetic inferences based on the analyses of the partial and entire mtDNA sequences. Importantly, the phylogenetic reconstruction based on the short marker sequences did not necessarily reproduce the tree topologies based on the analyses of the entire mtDNA. In addition, analyses on the datasets of various organisms revealed that the analyses based on the classic markers yielded phylogenetic trees with poor confidence in all tested cases compared to the results based on full-length mtDNA. These results demonstrated that phylogenetic analyses based on complete mtDNA sequences yield more insightful results compared to those based on mitochondrial genes and segments. To ameliorate the shortcomings of the classic markers, we identified a segment of mtDNA that may be used as an ‘approximate marker’ to closely reproduce the phylogenetic inference obtained from the entire mtDNA in the case of mammalian species, which can be utilized to design amplicon-seq-based studies. Our study demonstrates the importance of the choice of mitochondrial markers for phylogenetic analyses and proposes a novel approach to choosing appropriate markers for mammalian mtDNA that reproduces the phylogenetic inferences obtained from full-length mtDNA.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

Complete mitogenomes of Turkish tree squirrels, Sciurus anomalus and S. vulgaris, (Sciuridae: Rodentia: Mammalia) and their phylogenetic status within the tribe Sciurini

The genus Sciurus, a member of the family Sciuridae, is widely distributed in the Holarctic region. To better understand mitogenomic characteristics and to reveal internal phylogenetic relationships of the genus, 20 complete mitogenomes of Turkish tree squirrels were successfully sequenced for the first time, including 19 for S. anomalus (from 16,505 bp to 16,510 bp) and one for S. vulgaris (16,511 bp). The mitogenomes of two species were AT-biased. All tRNAs for two species displayed a typical clover-leaf structure, except for tRNA^Ser(AGY). The tRNA Serine1 (S1)-GCT structure lacked the dihydrouridine (DHU) loop and stem. Based on mitogenomic dataset for phylogeny of Sciurinae, phylogenetic analyses (Bayesian Inference and Maximum Likelihood) did not support monophyly of Sciurus and proposed that S. anomalus, the most basal taxa in the Sciurini tribe, had at least five mitogenome lineages, which were also supported by network analysis. The dissimilarities among the five linegaes of S. anomalus ranged from 0.0042 (0.42%) to 0.0062 (0.62%) using K2P sequence pairwise distances. In addition to this mitogenomic analysis result, phylogenetic analyses using the CYTB + D-loop dataset proposed the existence of at least nine lineages for S. anomalus, which was different than those of the previous studies. The current study proposed that the use of mitogenomic data for reconstructing the phylogeny of Turkey' Sciurus holds an important value for revealing evolutionary relationships.

Comparative mitogenomics of the genus Motacilla (Aves, Passeriformes) and its phylogenetic implications

The genus Motacilla belongs to Motacillidae (Passeriformes), where mitochondrial features are poorly understood and phylogeny is controversial. Whole mitochondrial genome (mitogenome) data and large taxon sampling are considered to be ideal strategies to obtain this information. We generated four complete mitogenomes of M.flava, M.cinerea, M.alba and Dendronanthusindicus, and made comparative analyses of Motacilla species combined with mitogenome data from GenBank, and then reconstructed phylogenetic trees based on 37 mitochondrial genes. The mitogenomes of four mitogenome sequences exhibited the same gene order observed in most Passeriformes species. Comparative analyses were performed among all six sampled Motacilla mitogenomes. The complete mitogenomes showed A-skew and C-skew. Most protein-coding genes (PCGs) start with an ATG codon and terminate with a TAA codon. The secondary structures of RNAs were similar among Motacilla and Dendronanthus. All tRNAs except for trnS(agy) could be folded into classic clover-leaf structures. Three domains and several conserved boxes were detected. Phylogenetic analysis of 90 mitogenomes of Passeriformes using maximum likelihood (ML) and Bayesian inference (BI) revealed that Motacilla was a monophyletic group. Among Motacilla species, M.flava and M.tschutschensis showed closer relationships, and M.cinerea was located in a basal position within Motacilla. These data provide important information for better understanding the mitogenomic characteristics and phylogeny of Motacilla.

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNA^Ile, ND1, tRNA^Leu2, 16S, tRNA^Val, 12S, tRNA^Phe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

First genomic resource for an endangered neotropical mega-herbivore: the complete mitochondrial genome of the forest-dweller (Baird's) tapir (Tapirus bairdii)

Baird's tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as 'endangered' on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG's). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 protein-coding genes (PCG's), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (O_L), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (O_H). A majority of the PCG's and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG's show values <1, suggesting that all these PCGs experience strong purifying selection. A maximum-likelihood phylogenetic analysis supports the monophyly of the genus Tapirus and the order Perissodactyla. The complete annotation and analysis of the mitochondrial genome of T. bairdii will contribute to a better understanding of the population genomic diversity and structure of this species, and it will assist in the conservation and protection of its dwindling populations.

Complete Mitochondrial Genomes of Five Racerunners (Lacertidae: Eremias) and Comparison with Other Lacertids: Insights into the Structure and Evolution of the Control Region

Comparative studies on mitochondrial genomes (mitogenomes) as well as the structure and evolution of the mitochondrial control region are few in the Lacertidae family. Here, the complete mitogenomes of five individuals of Eremias scripta (2 individuals), Eremias nikolskii, Eremias szczerbaki, and Eremias yarkandensis were determined using next-generation sequencing and were compared with other lacertids available in GenBank. The circular mitogenomes comprised the standard set of 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and a long non-coding control region (CR). The extent of purifying selection was less pronounced for the COIII and ND2 genes in comparison with the rest of the PCGs. The codons encoding Leucine (CUN), Threonine, and Isolecucine were the three most frequently present. The secondary structure of rRNA of Lacertidae (herein, E. scripta KZL15 as an example) comprised four domains and 28 helices for 12S rRNA, with six domains and 50 helices for 16S rRNA. Five types and twenty-one subtypes of CR in Lacertidae were described by following the criteria of the presence and position of tandem repeats (TR), termination-associated sequence 1 (TAS1), termination-associated sequence 2 (TAS2), conserved sequence block 1 (CBS1), conserved sequence block 2 (CSB2), and conserved sequence block 3 (CSB3). The compositions of conserved structural elements in four genera, Acanthodactylus, Darevskia, Eremias, and Takydromus, were further explored in detail. The base composition of TAS2 – TATACATTAT in Lacertidae was updated. In addition, the motif “TAGCGGCTTTTTTG” of tandem repeats in Eremias and the motif ”GCGGCTT” in Takydromus were presented. Nucleotide lengths between CSB2 and CSB3 remained 35 bp in Eremias and Darevskia. The phylogenetic analyses of Lacertidae recovered the higher-level relationships among the three subfamilies and corroborated a hard polytomy in the Lacertinae phylogeny. The phylogenetic position of E. nikolskii challenged the monophyly of the subgenus Pareremias within Eremias. Some mismatches between the types of CR and their phylogeny demonstrated the complicated evolutionary signals of CR such as convergent evolution. These findings will promote research on the structure and evolution of the CR and highlight the need for more mitogenomes in Lacertidae.

Starting the engine of the powerhouse: mitochondrial transcription and beyond

Mitochondria are central hubs for cellular metabolism, coordinating a variety of metabolic reactions crucial for human health. Mitochondria provide most of the cellular energy via their oxidative phosphorylation (OXPHOS) system, which requires the coordinated expression of genes encoded by both the nuclear (nDNA) and mitochondrial genomes (mtDNA). Transcription of mtDNA is not only essential for the biogenesis of the OXPHOS system, but also generates RNA primers necessary to initiate mtDNA replication. Like the prokaryotic system, mitochondria have no membrane-based compartmentalization to separate the different steps of mtDNA maintenance and expression and depend entirely on nDNA-encoded factors imported into the organelle. Our understanding of mitochondrial transcription in mammalian cells has largely progressed, but the mechanisms regulating mtDNA gene expression are still poorly understood despite their profound importance for human disease. Here, we review mechanisms of mitochondrial gene expression with a focus on the recent findings in the field of mammalian mtDNA transcription and disease phenotypes caused by defects in proteins involved in this process.

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.

Comparative Mitogenome Analysis of Gerbils and the Mitogenome Phylogeny of Gerbillinae (Rodentia: Muridae)

To enrich the mitogenomic database of Gerbillinae (Rodentia: Muridae), mitogenomes of three gerbils from different genera, Meriones tamariscinus (16,393 bp), Brachiones przewalskii (16,357 bp), and Rhombomys opimus (16,352 bp), were elaborated and compared with those of other gerbils in the present study. The three gerbil mitogenomes consisted of 2 ribosomal RNA genes, 13 protein-coding genes (PCGs), 22 transfer RNA genes, and one control region. Here, gerbil mitogenomes have shown unique characteristics in terms of base composition, codon usage, non-coding region, and the replication origin of the light strand. There was no significant correlation between the nucleotide percentage of G + C and the phylogenetic status in gerbils, and between the GC content of PCGs and the leucine count. Phylogenetic relationships of the subfamily Gerbillinae were reconstructed by 7 gerbils that represented four genera based on concatenated mitochondrial DNA data using both Bayesian Inference and Maximum Likelihood. The phylogenetic analysis indicated that M. tamariscinus was phylogenetically distant from the genus Meriones, but has a close relationship with R. opimus. B. przewalskii was closely related to the genus Meriones rather than that of R. opimus.

Complete mitogenome of the endangered and endemic Nicobar treeshrew (Tupaia nicobarica) and comparison with other Scandentians

The Nicobar treeshrew (Tupaia nicobarica) is an endangered small mammal endemic to the Nicobar Island of the Andaman Sea, India regarded as an alternative experimental animal model in biomedical research. The present study aimed to assemble the first mitochondrial genome of T. nicobarica to elucidate its phylogenetic position with respect to other Scandentians. The structure and variation of the novel mitochondrial genome were analyzed and compared with other Scandentians. The complete mitogenome (17,164 bp) encodes 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA (rRNAs), and one control region (CR). Most of the genes were encoded on majority strand, except nad6 and eight tRNAs. The nonsynonymous/synonymous ratio in all PCGs indicates strong negative selection among all Tupaiidae species. The comparative study of CRs revealed the occurrence of tandem repeats (CGTACA) found in T. nicobarica. The phylogenetic analyses (Maximum Likelihood and Bayesian Inference) showed distinct clustering of T. nicobarica with high branch supports and depict a substantial divergence time (12–19 MYA) from the ancestor lineage of Tupaiidae. The 16S rRNA dataset corroborates the taxonomic rank of two subspecies of T. nicobarica from the Great and Little Nicobar Islands. In the future, whole nuclear genome sequencing is necessary to further improve our understanding of evolutionary relationships among treeshrews, and will have implications for biomedical research.

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

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

Decreased mitochondrial D-loop region methylation mediates an increase in mitochondrial DNA copy number in CADASIL

Background

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a typical neurodegenerative disease associated with mitochondrial dysfunction. Methylation of the D-loop region and mitochondrial DNA copy number (mtDNAcn) play a critical role in the maintenance of mitochondrial function. However, the association between D-loop region methylation, mtDNAcn and CADASIL remains unclear.

Methods

Overall, 162 individuals were recruited, including 66 CADASIL patients and 96 age- and sex-matched controls. After extracting genomic DNA from the peripheral white blood cells, levels of D-loop methylation and mtDNAcn were assessed using MethylTarget sequencing and real-time PCR, respectively.

Results

We observed increased mtDNAcn and decreased D-loop methylation levels in CADASIL patients compared to the control group, regardless of gender stratification. Besides, we found a negative correlation between D-loop methylation levels and mtDNAcn. Mediation effect analysis shows that the proportion of the association between mtDNAcn and CADASIL that is mediated by D-loop methylation is 11.6% (95% CI 5.6, 22.6). After gender stratification, the proportions of such associations that are mediated by D-loop methylation in males and females were 7.2% (95% CI 2.4, 19.8) and 22.0% (95% CI 7.4, 50.1), respectively.

Conclusion

Decreased methylation of the D-loop region mediates increased mtDNAcn in CADASIL, which may be caused by a compensatory mechanism of mitochondrial dysfunction in patients with CADASIL.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01225-z.

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

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

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

The complete mitochondrial genome of the Honduran white bat Ectophylla alba (Allen 1982) (Chiroptera: Phyllostomidae)

The Honduran white bat, Ectophylla alba (Allen 1982), is one of eight species belonging to the family Phyllostomidae that exclusively roosts in tents. Due to its restricted distribution, habitat specificity, and diet requirements, E. alba has been strongly affected by habitat loss and fragmentation during the last decade. In this study, we developed the first genomic resource for this species; we assembled and analyzed in detail the complete mitochondrial genome of E. alba. The mitogenome of E. alba is 16,664 bp in length and is comprised of 13 protein coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a putative Control Region (CR) 1,232 bp in length. Gene arrangement in the mitochondrial chromosome of E. alba is identical to that reported before in other species of co-familiar bats. All PCGs are under purifying selection, with atp8 experiencing the least selective pressure. In all PCGs, codons ending with adenine are preferred over others ending in thymine and cytosine. Except tRNA-Serine 1, all tRNAs exhibit a cloverleaf secondary structure. The CR of E. alba exhibits three domains commonly described in other mammals, including bats; extended terminal associated sequences (ETAS), central, and conserved sequence block (CSB). A ML phylogenetic reconstruction of the family Phyllostomidae based on all 13 mitochondrial PCGs confirms the monophyletic status of the subfamily Sternodermatinae and indicates the close relationship between E. alba and the genus Artibeus. This is the first genomic resource developed for E. alba and represents the first step to improving our understanding of the genomic underpinnings involved in the evolution of specialization as well as acclimatization and adaptation to local and global change of specialist bats.

Genetic diversity and genetic differentiation of Megalobrama populations inferred by mitochondrial markers

BACKGROUND

Megalobrama is economically one of the most important freshwater fish genera in China. However, phylogenetic relationships among M. amblycephala, M. skolkovii, M. hoffmanni and M. pellegrini remain unresolved.

OBJECTIVE

To explore the genetic diversity and phylogenetic relationship of Megalobrama populations belonging to all four species.

METHODS

The concatenated sequences of mitochondrial cytochrome b (Cytb) and control region (CR) were used to analyze the genetic variation, genetic differentiation and population expansion of 15 Megalobrama populations.

RESULTS

The study showed that haplotype diversity and nucleotide diversity of M. hoffmanni and M. skolkovii were high, and that M. hoffmanni was the most genetically divergent of the four species. Haplotype network analysis revealed that M. hoffmanni and M. amblycephala formed a monophyletic group each, while M. skolkovii and M. pellegrini clustered together. There was a high genetic differentiation among the four Megalobrama species, and genetic distance among populations was not affected by geographical distance. Additionally, the results indicated that there was gene flow between the Liangzi Lake (LZL) population and Jinsha River Reservoir (JS) population. Also, Zhaoqing (ZQ) population of M. hoffmanni might have experienced a population expansion.

CONCLUSION

Our study verifies genetic diversity and differentiation of Megalobrama populations, and these findings will represent a significant contribution to the conservation and utilization of germplasm resources of Megalobrama.

The signals of selective constraints on the mitochondrial non-coding control region: insights from comparative mitogenomics of Clupeoid fishes

The vertebrate mitochondrial genome is characterized by an exceptional organization evolving towards a reduced size. However, the persistence of a non-coding and highly variable control region is against this evolutionary trend that is explained by the presence of conserved sequence motifs or binding sites for nuclear-organized proteins that regulate mtDNA maintenance and expression. We performed a comparative mitogenomic investigation of the non-coding control region to understand its evolutionary patterns in Clupeoid fishes which are widely distributed across oceans of the world, exhibiting exemplary evolutionary potential. We confirmed the ability of sequence flanking the conserved sequence motifs in the control region to form stable secondary structures. The existence of evolutionarily conserved secondary structures without primary structure conservation suggested the action of selective constraints towards maintaining the secondary structure. The functional secondary structure is maintained by retaining the frequency of discontinuous AT and TG repeats along with compensatory base substitutions in the stem forming regions which can be considered as a selective constraint. The nucleotide polymorphism along the flanking regions of conserved sequence motifs can be explained as errors during the enzymatic replication of secondary structure-forming repeat elements. The evidence for selective constraints on secondary structures emphasizes the role of the control region in mitogenome function. Maintenance of high frequency of discontinuous repeats can be proposed as a model of adaptive evolution against the mutations that break the secondary structure involved in the efficient regulation of mtDNA functions substantiating the efficient functioning of the control region even in a high nucleotide polymorphism environment.

Characterization and Comparative Analysis of Complete Mitogenomes of Three Cacatua Parrots (Psittaciformes: Cacatuidae)

Cacatua alba, Cacatua galerita, and Cacatua goffiniana are parrots of the family Cacatuidae. Wild populations of these species are declining with C. alba listed by the International Union for the Conservation of Nature and Natural Resources (IUCN) as Endangered. In this study, complete mitogenomes were sequenced for a comparative analysis among the Cacatua species, and a detailed analysis of the control region. Mitogenome lengths of C. alba,C. galerita, and C. goffiniana were 18,894, 18,900, and 19,084 bp, respectively. They included 13 protein coding genes, two ribosomal RNA genes, 24 transfer RNA genes, three degenerated genes, and two control regions. Ten conserved motifs were found in three domains within each of the two control regions. For an evolution of duplicated control regions of Cacatua, domain I and the 3′ end of domain III experienced an independent evolution, while domain II and most of the regions of domain III was subjected to a concerted evolution. Based on a phylogenetic analysis of 37 mitochondrial genes, the genus Cacatua formed a well-supported, monophyletic, crown group within the Cacatuidae. Molecular dating results showed that Cacatua diverged from other genera of Cacatuinae in the middle of Miocene.

Complete mitochondrial genome sequence of Lepus yarkandensis Günther, 1875 (Lagomorpha, Leporidae): characterization and phylogenetic analysis

Lepusyarkandensis is a national second-class protected animal endemic to China and distributed only in the hot and arid Tarim Basin in Xinjiang. We sequenced and described the complete mitogenome of L.yarkandensis to analyze its characteristics and phylogeny. The species’ DNA is a 17,047 bp circular molecule that includes 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. The overall base composition was as follows: A, 31.50%; T, 29.40%; G, 13.30% and C, 25.80%, with a high A+T bias of 60.9%. In the PCGs, ND6 had deviation ranges for AT skew (–0.303) and GC skew (0.636). The Ka/Ks values of ND1 (1.067) and ND6 (1.352) genes were >1, indicating positive selection, which might play an important role in the adaptation of L.yarkandensis to arid and hot environments. The conserved sequence block, the central conserved domain, and the extended termination-associated sequences of the control region and their features were identified and described. The phylogenetic tree based on the complete mitogenome showed that L.yarkandensis was closely related to the sympatric Lepustibetanuspamirensis. These novel datasets of L.yarkandensis can supply basic data for phylogenetic studies of Lepus spp., apart from providing essential and important resource for further genetic research and the protection of this species.

Modernizing the Toolkit for Arthropod Bloodmeal Identification

Simple Summary

The ability to identify the source of vertebrate blood in mosquitoes, ticks, and other blood-feeding arthropod vectors greatly enhances our knowledge of how vector-borne pathogens are spread. The source of the bloodmeal is identified by analyzing the remnants of blood remaining in the arthropod at the time of capture, though this is often fraught with challenges. This review provides a roadmap and guide for those considering modern techniques for arthropod bloodmeal identification with a focus on progress made in the field over the past decade. We highlight genome regions that can be used to identify the vertebrate source of arthropod bloodmeals as well as technological advances made in other fields that have introduced innovative new ways to identify vertebrate meal source based on unique properties of the DNA sequence, protein signatures, or residual molecules present in the blood. Additionally, engineering progress in miniaturization has led to a number of field-deployable technologies that bring the laboratory directly to the arthropods at the site of collection. Although many of these advancements have helped to address the technical challenges of the past, the challenge of successfully analyzing degraded DNA in bloodmeals remains to be solved.

Abstract

Understanding vertebrate–vector interactions is vitally important for understanding the transmission dynamics of arthropod-vectored pathogens and depends on the ability to accurately identify the vertebrate source of blood-engorged arthropods in field collections using molecular methods. A decade ago, molecular techniques being applied to arthropod blood meal identification were thoroughly reviewed, but there have been significant advancements in the techniques and technologies available since that time. This review highlights the available diagnostic markers in mitochondrial and nuclear DNA and discusses their benefits and shortcomings for use in molecular identification assays. Advances in real-time PCR, high resolution melting analysis, digital PCR, next generation sequencing, microsphere assays, mass spectrometry, and stable isotope analysis each offer novel approaches and advantages to bloodmeal analysis that have gained traction in the field. New, field-forward technologies and platforms have also come into use that offer promising solutions for point-of-care and remote field deployment for rapid bloodmeal source identification. Some of the lessons learned over the last decade, particularly in the fields of DNA barcoding and sequence analysis, are discussed. Though many advancements have been made, technical challenges remain concerning the prevention of sample degradation both by the arthropod before the sample has been obtained and during storage. This review provides a roadmap and guide for those considering modern techniques for arthropod bloodmeal identification and reviews how advances in molecular technology over the past decade have been applied in this unique biomedical context.

Complete Mitochondrial Genome of Three Species of the Genus Microtus (Arvicolinae, Rodentia)

The 65 species of the genus Microtus have unusual sex-related genetic features and a high rate of karyotype variation. However, only nine complete mitogenomes for these species are currently available. We describe the complete mitogenome sequences of three Microtus, which vary in length from 16,295 bp to 16,331 bp, contain 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes and a control region. The length of the 13 PCGs and the coded proteins is the same in all three species, and the start and stop codons are conserved. The non-coding regions include the L-strand origin of replication, with the same sequence of 35 bp, and the control region, which varies between 896 bp and 930 bp in length. The control region includes three domains (Domains I, II and III) with extended termination-associated sequences (ETAS-1 and ETAS-2) in Domain I. Domain II and Domain III include five (CSB-B, C, D, E and F) and three (CSB-1, CSB-2, and CSB-3) conserved sequence blocks, respectively. Phylogenetic reconstructions using the mitochondrial genomes of all the available Microtus species and one representative species from another genus of the Arvicolinae subfamily reproduced the established phylogenetic relationships for all the Arvicolinae genera that were analyzed.

Complete mitochondrial genome of a blue-tailed skink Plestiodon capito (Reptilia, Squamata, Scincidae) and comparison with other Scincidae lizards

Vertebrate mitochondrial genomes (mitogenomes) are valuable for studying phylogeny, evolutionary genetics and genomics. To date, however, compared to other vertebrate groups, our knowledge about the mitogenomes of skinks (the family Scincidae), even of reptile, has been relatively limited. In the present study, we determined the complete mitogenome of a blue-tailed skink Plestiodon capito for the first time, and compared it with other skinks available in GenBank. The circular genome is 17,344 bp long, showing a typical vertebrate pattern with 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one control region (CR). The gene organization, nucleotide composition, and codon usage are similar to those from skinks previously published. Twelve out of 13 PCGs initiates with canonical start codon (ATG), while COX1 starts with GTG. The codon usage analysis revealed a preferential use of the Leu^CUN (Leu1), Pro, and Thr codons with the A/U ending. All tRNAs in P. capito were predicted to fold into typical clover-leaf secondary structure, except tRNA-Ser ^AGY. The secondary structures of 12S rRNA and 16S rRNA comprises 34 helices and 56 helices, respectively. The alignment of the Plesitodon species CRs exhibited high genetic variability and rich A + T content. Besides, variable types and numbers of tandem repeat units were also identified in the CR of Plestiodon. Phylogenetic analyses recovered P. capito as the sister species to P. tunganus; monophyly of the Scincidae is well supported. Our results will help to better understand structure and evolution of the mitochondrial DNA control region in reptiles as well as the evolutionary status of P. capito, and to lay foundation for further phylogenetic study of skinks in a mitogenomic framework.

Characterization of three new mitochondrial genomes of Coraciiformes (Megaceryle lugubris, Alcedo atthis, Halcyon smyrnensis) and insights into their phylogenetics

Coraciiformes contains more than 200 species with great differences on external morphology and life-style. The evolutionary relationships within Coraciiformes and the phylogenetic placement of Coraciiformes in Aves are still questioned. Mitochondrial genome (mitogenome) sequences are popular markers in molecular phylogenetic studies of birds. This study presented the genome characteristics of three new mitogenomes in Coraciiformes and explored the phylogenetic relationships among Coraciiformes and other five related orders with mitogenome data of 30 species. The sizes of three mitogenomes were 17,383 bp (Alcedo atthis), 17,892 bp (Halcyon smyrnensis) and 17,223 bp (Megaceryle lugubris). Each mitogenome contained one control region and 37 genes that were common in vertebrate mitogenomes. The organization of three mitogenomes was identical to the putative ancestral gene order in Aves. Among 13 available Coraciiform mitogenomes, 12 protein coding genes showed indications of negative selection, while the MT-ND6 presented sign of positive selection or relaxed purifying selection. The phylogenetic results supported that Upupidae and Bucerotidae should be separated from Coraciiformes, and that Coraciiformes is more closely related to Piciformes than to Strigiformes, Trogoniformes and Cuculiformes. Our study provide valuable data for further phylogenetic investigation of Coraciiformes.

Wildlife forensics: A boon for species identification and conservation implications

Wildlife trade and fraudulence in food, artefacts and cosmetic industries had raised serious concern in protection of the wild faunal diversity. Lack of proper tools and molecular based techniques for identification of wild species are some of the major constrains faced by the judiciary and law enforcement agencies while framing charges against poachers and illicit agitator. The emergence of wildlife forensics serves as a boon in solving long pending cases of wildlife crimes. Wildlife forensics have proven to be fast, accurate and reliable criminal investigation processes with comprehensive coverage and easy accessibility. It has also helped resolving taxonomic disputes, determining spatiotemporal genetic divergence, evolutionary history, origins and even endemism. Collaboration among inter-disciplinary fields has even led to engineered signature markers and phylogenetics for several species. Development in fields of genetics, molecular and evolutionary biology and other omics techniques have further contributed in accurate identification of species. Wildlife forensics, with the support of proper international mega database units for population reference, will be fundamental in wildlife investigations through its unlimited information sharing ability. The efficient conservation of species will, however, require a collaborative approach consisting of national policy makers, local stakeholders and implementation agencies in addition to experts from the scientific communities.

Structural Variation of the Turtle Mitochondrial Control Region

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.