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

Selective sweep analysis of the adaptability of the Yarkand hare (Lepus yarkandensis) to hot arid environments using SLAF-seq

The Yarkand hare (Lepus yarkandensis) inhabits arid desert areas and is endemic to China. It has evolved various adaptations to survive in hot arid environments, including stress responses, the ability to maintain water homeostasis and heat tolerance. Here, we performed a selective sweep analysis to identify the candidate genes for adaptation to hot arid environments in the Yarkand hare. A total of 397 237 single-nucleotide polymorphisms were obtained from 80 Yarkand hares, which inhabit hot arid environments, and 36 Tolai hares (Lepus tolai), which inhabit environments with a mild climate, via specific-locus amplified fragment sequencing. We identified several candidate genes that were associated with the heat stress response (HSPE1), oxidative stress response (SLC23A and GLRX2), immune response (IL1R1 and IRG1), central nervous system development (FGF13, THOC2, FMR1 and MECP2) and regulation of water homeostasis (CDK1) according to fixation index values and θπ ratios in the selective sweep analysis, and six of these genes (GLRX2, IRG1, FGF13, FMR1, MECP2 and CDK1) are newly discovered genes. To the best of our knowledge, this is the first study to identify candidate genes for adaptation to hot arid environments in the Yarkand hare. The results of this study enhance our understanding of the adaptation of the Yarkand hare to hot arid environments and will aid future studies aiming to functionally verify these candidate genes.

Sequencing and description of the complete mitochondrial genome of Limatus durhamii (Diptera: Culicidae)

The genus Limatus (Diptera: Culicidae) are wild mosquitoes belonging to the Sabethini tribe that occurs in tropical countries and is related to transmission cycles of Orthobunyavirus (Bunyaviridae), particularly in the Amazon region. Given the unavailability of information related to evolutionary biology and molecular taxonomy aspects of this genus, we report here the first complete sequencing of the mitochondrial genome of Limatus durhamii Theobald, 1901. The NextSeq 500 platform was used for sample sequencing, and the mitochondrial sequence obtained was 14,875 bp long, comprising 37 functional subunits (13 PCGs, 22 tRNA and 02 rRNA). The phylogeny reconstructed by maximum likelihood based on the concatenation of all 13 PCGs corroborated the known taxonomic classification based most on aspects of the external morphology and few molecular studies. The data and information produced here may be useful in the future development of taxonomic and evolutionary studies for the genus, as well as the Culicidae family itself.

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

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

Description of mitochon genome and phylogenetic considerations of Sabethes bipartipes, Sabethes cyaneus, Sabethes quasicyaneus, and Sabethes tarsopus (Diptera: Culicidae)

The genus Sabethes (Diptera: Culicidae) comprises species of great epidemiological relevance, particularly involved in transmission cycles of the Yellow fever virus in South America. Given the unavailability of information related to aspects of evolutionary biology and molecular taxonomy of species of this genus of mosquitoes, we report here the first sequencing of the mitochondrial genomes of Sabethes bipartipes, Sabethes cyaneus, Sabethes tarsopus, and Sabethes quasicyaneus. The sequences obtained showed an average length of 14,920 bp, comprising 37 functional genes (13 PCGs, 22 tRNA, and 02 rRNA). The phylogenies reconstructed by Maximum likelihood and Bayesian inference methods, based on the concatenated sequences of all 13 PCGs, produced similar topologies and strongly supported the monophyletic relationship between the Sabethes subgenera, corroborating the known taxonomic classification based on aspects of the external morphology of the taxa assessed. The data and information produced from the Sabethes species evaluated here may be useful for future taxonomic and evolutionary studies of the genus, as well as the Culicidae family.

Expression Regulation of Water Reabsorption Genes and Transcription Factors in the Kidneys of Lepus yarkandensis

Lepus yarkandensis is a desert-dwelling animal that has various adaptations to cope with drought. The kidney maintains water and acid-base balance mainly through the vasopressin-regulated water reabsorption pathway and proximal tubular bicarbonate reabsorption pathway. In this study, we compared the differentially expressed genes (DEGs) and transcription factors in the kidneys of L. yarkandensis and Oryctolagus cuniculus to explore the relationship between the DEGs in kidneys and the animals’ adaptations. Transcriptome sequencing data were used to predict the differentially-expressed water reabsorption genes and their transcription factors. Quantitative real-time PCR, immunohistochemistry, and western blotting were used to detect and verify the expression of DEGs in the kidney at mRNA and protein levels. Transcriptome analysis of the kidney of L. yarkandensis and O. cuniculus showed that 6,610 genes were up-regulated and 5,727 genes down-regulated in data shared by both species. According to the data, 232 transcription factors and their corresponding target genes were predicted, from which genes and transcription factors related to renal water reabsorption were screened. Quantitative RT-PCR results showed AQP1, AQP2, ADCY3, HIF1A, CREB3, and NFATc1 had higher expression in the L. yarkandensis kidney; in comparison, FXYD2 mRNA expression levels were lower. In western blotting, transcription factors HIF1A, NFATc1, NF-κB1, and critical genes ADCY3, ATPA1, and SLC4A4, were highly expressed in the kidneys of L. yarkandensis. Immunohistochemical results showed that the ADCY3 protein was in the basolateral membrane of the collecting duct, the ATP1A1 protein was in the basolateral membrane and medulla of proximal tubules, and the SLC4A4 protein was in the basolateral membrane of proximal tubules. According to these results can be inferred that HIF1A, NFATc1, and NF-κB1 play a certain role in regulating the expression of genes related to water reabsorption in the kidney of L. yarkandensis, thus improving the water reclamation efficiency of L. yarkandensis, so as to adapt to the arid desert environment.

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.