Sequence and phylogenetic analysis of the complete mitochondrial genome of Lasiopodomys mandarinus mandarinus (Arvicolinae, Rodentia)

Novel insights into adaptive evolution based on the unusual AT-skew in Acheilognathus gracilis mitogenome and phylogenetic relationships of bitterling

Acheilognathus gracilis, a bitterling species, distribute in lower reaches of Yangtze River. They are identified as the top-priority bitterling species for conservation as having high evolutionary distinctiveness and are at risk of extinction. In present study, we first sequenced the complete mitogenome of A. gracilis and analyzed its phylogenetic position using 13 PCGs. The A. gracilis mitogenome is 16,774 bp in length, including 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, a control region and the origin of the light strand replication. The overall base composition of A. gracilis in descending order is T 27.9 %, A 27.7 %, C 26.1 % and G 18.3 %, shows a unusual AT-skew with slightly negative. Further investigation revealed A. gracilis uses excess T over A in NADH dehydrogenase 5 (nd5), whereas the most of other bitterlings are biased toward to use A not T, implying there is likely to be unique strategy of adaptive evolution in A. gracilis. We also compared 13 PCGs of 30 bitterling mitogenomes and the results exhibit highly conservative. Phylogenetic trees constructed by 13 PCGs strongly support the monophyly of Acheilognathus and the paraphyly of Rhodeus and Tanakia. Current results will provide valuable information for follow-up research on conservation of species facing with serious population decline and can provide novel insights into the phylogenetic analysis and evolutionary biology research.

Phylogeny and evolution of Lasiopodomys in subfamily Arvivolinae based on mitochondrial genomics

The species of Lasiopodomys Lataste 1887 with their related genera remains undetermined owing to inconsistent morphological characteristics and molecular phylogeny. To investigate the phylogenetic relationship and speciation among species of the genus Lasiopodomys, we sequenced and annotated the whole mitochondrial genomes of three individual species, namely Lasiopodomys brandtii Radde 1861, L. mandarinus Milne-Edwards 1871, and Neodon (Lasiopodomys) fuscus Büchner 1889. The nucleotide sequences of the circular mitogenomes were identical for each individual species of L. brandtii, L. mandarinus, and N. fuscus. Each species contained 13 protein-coding genes (PCGs), 22 transfer RNAs, and 2 ribosomal RNAs, with mitochondrial genome lengths of 16,557 bp, 16,562 bp, and 16,324 bp, respectively. The mitogenomes and PCGs showed positive AT skew and negative GC skew. Mitogenomic phylogenetic analyses suggested that L. brandtii, L. mandarinus, and L. gregalis Pallas 1779 belong to the genus Lasiopodomys, whereas N. fuscus belongs to the genus Neodon grouped with N. irene. Lasiopodomys showed the closest relationship with Microtus fortis Büchner 1889 and M. kikuchii Kuroda 1920, which are considered as the paraphyletic species of genera Microtus. T_MRCA and niche model analysis revealed that Lasiopodomys may have first appeared during the early Pleistocene epoch. Further, L. gregalis separated from others over 1.53 million years ago (Ma) and then diverged into L. brandtii and L. mandarinus 0.76 Ma. The relative contribution of climatic fluctuations to speciation and selection in this group requires further research.

Whole mitochondrial genome of long-clawed mole vole (Prometheomys schaposchnikowi) from Turkey, with its phylogenetic relationships

The mitogenome of Prometheomys schaposchnikowi was characterized for the first time as a circular DNA molecule (16.284 bp), containing 37 coding and 2 non-coding regions. In the mitogenome, ND6 and 8 tRNA genes were encoded on the light chain, while 12 PCGs, 14 tRNAs, 2 rRNAs, D-loop and O_L were encoded on the heavy chain. The most common initiation codon in PCGs was ATG. As in many mammals, incomplete stop codons in P. schaposchnikowi were in the COX3, ND1 and ND4. Phylogenetic relationships were revealed using Bayesian method and the 13 PCGs. Seven genera (Arvicola, Dicrostonyx, Lasiopodomys, Myodes, Ondatra, Proedromys and Prometheomys) formed a monophyletic group, while Eothenomys, Microtus and Neodon were paraphyletic. P. schaposchnikowi constituted the most basal group within Arvicolinae. Divergence time estimation suggested that P. schaposchnikowi diversified during the Miocene (16.28 Mya). Further molecular studies are needed to test the distinctiveness and diversity of the genus Prometheomys.

Evolution of Hemoglobin Genes in a Subterranean Rodent Species (Lasiopodomys mandarinus)

The Mandarin vole (Lasiopodomys mandarinus), a typical subterranean rodent, has undergone hematological adaptations to tolerate the hypoxic/hypercapnic underground environment. Hemoglobin (Hb) genes encode respiratory proteins functioning principally in oxygen binding and transport to various tissues and organs. To investigate the evolution of α- and β-hemoglobin (Hb) in subterranean rodent species, we sequenced Hb genes of the Mandarin vole and the related aboveground Brandt's vole (L. brandtii). Sequencing showed that in both voles, α-globin was encoded by a cluster of five functional genes in the following linkage order: HBZ, HBA-T1, HBQ-T1, HBA-T2, and HBQ-T2; among these, HBQ-T2 is a pseudogene in both voles. The β-globin gene cluster in both voles also included five functional genes in the following linkage order: HBE, HBE/HBG, HBG, HBB-T1, and HBB-T2. Phylogenetic analysis revealed that the Mandarin vole underwent convergent evolution with its related aboveground species (Brandt's vole) but not with other subterranean rodent species. Selection pressure analyses revealed that α- and β-globin genes are under strong purifying selection (ω < 1), and branch-site analyses identified positive selection sites on HBAT-T1 and HBB-T1 in different subterranean rodent species. This suggests that the adaptive evolution of these genes enhanced the ability of Hb to store and transport oxygen in subterranean rodent species. Our findings highlight the critical roles of Hb genes in the evolution of hypoxia tolerance in subterranean rodent species.

Expression patterns of clock genes in the hypothalamus and eye of two Lasiopodomys species

To investigate the relationship between light sensing systems in the eye and circadian oscillators in the hypothalamus of subterranean rodents, we studied subterranean Mandarin voles (Lasiopodomys mandarinus) that spend their entire lives under dark conditions with degenerated eyes, and compared oscillatory expression patterns of clock genes in the hypothalamus and eye between Mandarin voles and their aboveground relatives, Brandt's voles (L. brandtii). Individuals of both vole species were kept under a 12-h light/12-h dark condition or continuous dark condition for 4 weeks. In both species, the expressions of most genes showed significant cosine rhythmicity in the hypothalamus but relatively weak rhythmicity in the eye. The number of rhythmic genes in the eye of Mandarin voles increased under the dark condition, but the opposite trend was observed in the eye of Brandt's voles. The expression levels of most clock genes in the hypothalamus of both vole species did not significantly differ between the two conditions, but unlike in Mandarin voles, these expression levels significantly decreased in the eye of Brandt's voles kept under the dark condition. In both vole species, the peak phase of most clock genes exhibited advanced or invariant change in the hypothalamus under the dark condition, and the peak phase of most clock genes showed consistent changes between the eye and hypothalamus of Mandarin voles. However, most clock genes in the eye showed a delayed phase in Brandt's voles kept under the dark condition. In conclusion, the hypothalamus plays an important role in both vole species irrespective of the light condition. However, the expression patterns of clock genes in the eye differed between the vole species, indicating that each species adapted differently to their environments.

Transcriptome analysis of the response provided by Lasiopodomys mandarinus to severe hypoxia includes enhancing DNA repair and damage prevention

Background

Severe hypoxia induces a series of stress responses in mammals; however, subterranean rodents have evolved several adaptation mechanisms of energy metabolisms and O₂ utilization for hypoxia. Mammalian brains show extreme aerobic metabolism. Following hypoxia exposure, mammals usually experience irreversible brain damage and can even develop serious diseases, such as hypoxic ischemic encephalopathy and brain edema. To investigate mechanisms underlying the responses of subterranean rodents to severe hypoxia, we performed a cross-species brain transcriptomic analysis using RNA sequencing and identified differentially expressed genes (DEGs) between the subterranean rodent Lasiopodomys mandarinus and its closely related aboveground species L. brandtii under severe hypoxia (5.0% O₂, 6 h) and normoxia (20.9% O₂, 6 h).

Results

We obtained 361 million clean reads, including 69,611 unigenes in L. mandarinus and 69,360 in L. brandtii. We identified 359 and 515 DEGs by comparing the hypoxic and normoxia groups of L. mandarinus and L. brandtii, respectively. Gene Ontology (GO) analysis showed that upregulated DEGs in both species displayed similar terms in response to severe hypoxia; the main difference is that GO terms of L. brandtii were enriched in the immune system. However, in the downregulated DEGs, GO terms of L. mandarinus were enriched in cell proliferation and protein transport and those of L. brandtii were enriched in nuclease and hydrolase activities, particularly in terms of developmental functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that upregulated DEGs in L. mandarinus were associated with DNA repair and damage prevention as well as angiogenesis and metastasis inhibition, whereas downregulated DEGs were associated with neuronal synaptic transmission and tumor-associated metabolic pathways. In L. brandtii, upregulated KEGG pathways were enriched in the immune, endocrine, and cardiovascular systems and particularly in cancer-related pathways, whereas downregulated DEGs were associated with environmental information processing and misregulation in cancers.

Conclusions

L. mandarinus has evolved hypoxia adaptation by enhancing DNA repair, damage prevention, and augmenting sensing, whereas L. brandtii showed a higher risk of tumorigenesis and promoted innate immunity toward severe hypoxia. These results reveal the hypoxic mechanisms of L. mandarinus to severe hypoxia, which may provide research clues for hypoxic diseases.

Mitochondrial genome of the garfish Hyporhamphus quoyi (Beloniformes: Hemiramphidae) and phylogenetic relationships within Beloniformes based on whole mitogenomes

Mitochondrial DNA (mtDNA) can provide genome-level information (e.g. mitochondrial genome structure, phylogenetic relationships and codon usage) for analyzing molecular phylogeny and evolution of teleostean species. The species in the order Beloniformes have commercial importance in recreational fisheries. In order to further clarify the phylogenetic relationship of these important species, we determined the complete mitochondrial genome (mitogenome) of garfish Hyporhamphus quoyi of Hemiramphidae within Beloniformes. The mitogenome was 16,524 bp long and was typical of other teleosts mitogenomes in size and content. Thirteen PCGs started with the typical ATG codon (with exception of the cytochrome coxidase subunit 1 (cox1) gene with GTG). All tRNA sequences could be folded into expected cloverleaf secondary structures except for tRNASer (AGN) which lost a dihydrouracil (DHU) stem. The control region was 866 bp in length, which contained some conserved sequence blocks (CSBs) common to Beloniformes. The phylogenetic relationship between 26 fish Beloniformes species was analyzed based on the complete nucleotide and amino acid sequences of 13 PCGs by two different inference methods (Maximum Likelihood and Bayesian Inference). Phylogenetic analyses revealed Hemiramphidae as the sister group to Exocoetidae and it is a paraphyletic grouping. Our results may provide useful information on mitogenome evolution of teleostean species.

The first two complete mitochondrial genomes for the family Triglidae and implications for the higher phylogeny of Scorpaeniformes

The mitochondrial genome (mitogenome) can provide useful information for analyzing phylogeny and molecular evolution. Scorpaeniformes is one of the most diverse teleostean orders and has great commercial importance. To develop mitogenome data for this important group, we determined the complete mitogenomes of two gurnards Chelidonichthys kumu and Lepidotrigla microptera of Triglidae within Scorpaeniformes for the first time. The mitogenomes are 16,495 bp long in C. kumu and 16,610 bp long in L. microptera. Both the mitogenomes contain 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and two non-coding regions. All PCGs are initiated by ATG codons, except for the cytochrome coxidase subunit 1 (cox1) gene. All of the tRNA genes could be folded into typical cloverleaf secondary structures, with the exception of tRNASer(AGN) lacks a dihydrouracil (DHU) stem. The control regions are both 838 bp and contain several features common to Scorpaeniformes. The phylogenetic relationships of 33 fish mitogenomes using Bayesian Inference (BI) and Maximum Likelihood (ML) based on nucleotide and amino acid sequences of 13 PCGs indicated that the mitogenome sequences could be useful in resolving higher-level relationship of Scorpaeniformes. The results may provide more insight into the mitogenome evolution of teleostean species.