Abstract
During the evolutionary transition from land to water, cetaceans have undergone numerous critical challenges, with osmoregulation being the major one. Two subspecies of the narrow-ridged finless porpoise (Neophocaena asiaeorientalis), the freshwater Yangtze finless porpoise (N. a. asiaeorientalis, NAA) and the marine East Asian finless porpoise (N. a. sunameri, NAS), provide excellent subjects to understand the genetic basis of osmoregulatory divergence between freshwater and marine mammals. The kidney plays an important and well-established role in osmoregulation in marine mammals and thus, herein, we utilized RNA-seq to characterize the renal transcriptome and preliminarily analyze the divergence between the NAA and the NAS. Approximately 48.98 million clean reads from NAS and 49.40 million clean reads from NAA were obtained by RNA-Seq. And 73,449 (NAS) and 68,073 (NAA) unigenes were assembled. Among these annotations, 22,231 (NAS) and 21,849 (NAA) unigenes were annotated against the NCBI nr protein database. The ion channel complex GO term and four pathways were detected as relevant to osmoregulation by GO and KEGG pathway classification of these annotated unigenes. Although the endangered status of the study species prevented analysis of biological replicates, we identified nine differentially expressed genes (DEGs) that may be vital in the osmoregulation of the narrow-ridged finless porpoise and worthwhile for future studies. Of these DEGs, the differential expression and distribution of the aquaporin-2 (AQP2) in the collecting duct were verified using immunohistochemical experiments. Together, this work is the first report of renal transcriptome sequencing in cetaceans, and it will provide a valuable resource for future molecular genetics studies on cetacean osmoregulation.