Skip to main content
  • Research article
  • Open access
  • Published:

Key genes differential expressions and pathway involved in salt and water-deprivation stresses for renal cortex in camel



Camels possess the characteristics of salt- and drought-resistances, due to the long-time adaption to the living environment in desert. The camel resistance research on transcriptome is rare and deficient, especially reabsorption in renal cortex. Non-coding RNAs are normally considered as the RNA molecules that are not translated into proteins, their current roles remain mostly in regulation of information flux from DNA to protein, further on normal life activities and diseases. In order to reveal the mysterious veil of the post-transcriptional regulation of ncRNAs in renal cortex for the first time as far as we know, we designed and carried out the experiment of salt stress and water-deprivation stress in camel.


By means of RNA-seq in renal cortex of Alxa Bactrian Camel (Camelus bactrianus), we identified certain significantly differential RNAs, including 4 novel lncRNAs, 11 miRNAs and 13 mRNAs under salt stress, 0 lncRNAs, 18 miRNAs and 14 mRNAs under water-deprivation stress. By data analysis, the response pathway of post-transcriptional regulation concerning salt and water-deprivation stresses was put forward, involving preventing sodium from entering the cell, purifying of water and compensating neutral amino acids by miR-193b, miR-542-5p interaction with SLC6A19 mRNA.


Based on the resistance-related lncRNAs, miRNAs, and mRNAs, we proposed the post-transcriptional regulation pathway to explain how camels respond to salt and water-deprivation stresses in the ncRNAs regulation level of renal cortex for the first time, thus hoping to provide a theoretical basis for therapy of disease that is similar to high blood pressure in humans.


Camel, a kind of animal that living in the desert, has a series of anti-stress characteristics such as the resistance of heat, salt and drought [1,2,3]. The previous resistance researches of camel have reported that plasma osmolality and renal sodium excretion increased after the saline loading [4], while water-deprivation is correlated with rise in plasma sodium concentration, urine osmolality and sodium excretion, and with concomitant decrease in plasma volume and urine production [5]. Based on the studies in humans and other model species, salt metabolism is regulated by neuroregulation and humoral regulation, especially the two processes of intestinal absorption and renal reabsorption. Either in intestinal epithelial cells or in renal tubules, Na+ absorption requires proteins, such as solute carrier family 26 member 3, 6 and 9, EnAC complex and Na+/K+ ATP enzymes, and Na+ absorption is regulated by renin, angiotensin, aldosterone, vasopressin and other hormones in body fluids [6]. In kidney, the formation and excretion of urine are involved in renal cortex and then renal medulla [7]. The solute reabsorption containing Na+ reabsorption is involved in the function of renal cortex with high hydraulic conductivity of cortical peritubular capillaries [8, 9]. The water reservation mechanism of renal medulla in camel has also been revealed in our previous study [10]. But the molecular resistance regulation of these crucial features is rarely demonstrated under salt and water-deprivation stress in the renal cortex, especially in the level of non-coding RNA transcriptional regulation. Non-coding RNAs (ncRNAs) are normally considered as the RNA molecules that are not translated into proteins, their current roles remain mostly in regulation of information flux from DNA to protein [11,12,13]. In summary, we selected long non-coding RNA (lncRNA) and micro RNA (miRNA), the two kinds of ncRNA, for the investigation of regulatory role of mRNA, to explore the camel’s regulation in response to the environmental stress of salt and water-deprivation. Furthermore, concerning that the resistance-related genes and pathway may be interrelated between camel and human, we hope to provide inspirations of diagnostic and therapeutic strategies for human diseases such as high blood pressure from scientific experimental evidences.


Differently expressed protein-coding genes under salt and water-deprivation stresses

Salt resistance-related genes with significantly different expression were involved in 13 genes, including up-regulated CORO1C, LOC105076684, PICALM, HERC4, LOC105071055, HNRNPK and down-regulated WIPF2, ERBB3, ZC3H7A, TOM1L1, SLC6A19, TENM1, LOC105078699 through RNA-seq analysis of camel renal cortex (Fig. 1a–c; Additional file 1: Table S1). Under water-deprivation stress, there were 14 differently expressed genes which referred to up-regulated CORO1C, KTN1, CDH11, PICALM, HERC4, LOC105071055, HNRNPK and down-regulated SCMH1, SCAF4, HSPA6, SLC6A19, LOC105061856, LOC105078699, GSTT2 in renal cortex of camel (Fig. 1a, b, d; Additional file 1: Table S2). Any pathways associated with sodium/water-metabolism were not sorted out in Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis, except Mineral absorption with SLC6A19 (Fig. 1e, f). Moreover, SLC6A19 was also discovered in 7 overlapped genes under salt and water-deprivation stresses using Venn diagram (Fig. 1g).

Fig. 1
figure 1

Differentially expressed mRNAs, miRNAs, lncRNAs and the analysis of protein-coding genes within KEGG pathways. a Values of log2(fold change) of up-regulated mRNAs and miRNAs under salt and water-deprivation stresses. b Values of − log2(fold change) of down-regulated mRNAs, miRNAs and lncRNAs under salt and water-deprivation stresses. c Number of differentially expressed protein-coding genes under salt stress. d Number of differentially expressed protein-coding genes under water-deprivation stress. e Enriched KEGG pathways of protein-coding genes in the renal cortex under salt stress. f Enriched KEGG pathways of protein-coding genes in the renal cortex under water-deprivation stress. g Seven overlapped protein-coding genes of differential expression under salt and water-deprivation stresses

miRNAs identification and candidate interaction with miRNA–mRNA

The significantly differential expression of miRNAs under salt and water-deprivation stresses was obtained by RNA-seq and sequence alignment with Bos Taurus’s miRNA (Fig. 2a, b; Additional file 1: Tables S1, S2). Four significantly up-regulated miRNAs (miR-365-3p, miR-193b, miR-664b and miR-483) and seven significantly down-regulated miRNAs (miR-101, miR-143, miR-378c, miR-122, miR-452, miR-1246 and novel 1) were detected in the camel renal cortex under salt stress (Fig. 1a, b). Eight significantly up-regulated miRNAs (miR-125a, miR-92a, miR-204, miR-196a, miR-155, miR-424-3p, miR-224 and miR-30f) and ten significantly down-regulated miRNAs (miR-101, miR-378, miR-133a, miR-378c, miR-106b, miR-19b, miR-133b, miR-3601, miR-542-5p and novel 1) were discovered in the camel renal cortex under water-deprivation stress (Fig. 1a, b). Mineral absorption was enriched by KEGG pathway analysis of target mRNAs in terms of sodium/water-metabolism (Fig. 2c, d). Nevertheless, SLC6A19 mRNA was a corresponding target of miR-193b and miR-542-5p by miRNA target prediction using miRanda.

Fig. 2
figure 2

Number of differentially expressed miRNAs and KEGG pathways of target mRNAs. a Number of differentially expressed miRNAs under salt stress. b Number of differentially expressed miRNAs under water-deprivation stress. c Enriched KEGG pathways of target mRNAs in the renal cortex under salt stress. d Enriched KEGG pathways of target mRNAs in the renal cortex under water-deprivation stress

Novel lncRNAs under salt stress and water-deprivation stress

For camel renal cortex under salt stress and water-deprivation stress, the differentially expressed lncRNA was presented by RNA-seq (Fig. 3a, b; Additional file 1: Tables S1, S2). The data of RNA-seq indicated that four novel significantly down-regulated lncRNAs (LNC002600, LNC000062, LNC001899 and LNC000331) were detected under salt stress, utilizing sequence alignment with Bos Taurus’s lncRNA (Figs. 1a, 3c; Additional file 1: Table S3). Whereas, there was no significant difference in the expression of lncRNAs under water-deprivation stress (Fig. 3b). The co-located and co-expressed mRNAs were predicted, and KEGG pathway analysis of the mRNAs was carried out accordingly with the novel lncRNAs by bioinformatics analysis (Fig. 3d, e). There are two explanations for the RNA-seq results that differentially expressed lncRNAs were not detected under water-deprivation stress. Firstly, the renal cortex may not be a major regulation section of water-deprivation. Secondly, we did not sequence circular RNAs (circRNAs) which may be differentially expressed and execute essential bio-function, not lncRNAs.

Fig. 3
figure 3

Four novel lncRNAs and KEGG pathways of co-located and co-expressed mRNAs. a Number of differentially expressed lncRNAs under salt stress. b Number of differentially expressed lncRNAs under water-deprivation stress. c Detail genes information of differentially expressed lncRNAs. d, e Enriched KEGG pathways of co-located and co-expressed mRNAs in the renal cortex under salt stress


In order to explore the corporate regulation of camel to deal with the salt and water-deprivation stresses, we found a directly related mRNA with sodium transport in the common 7 protein-coding genes, namely SLC6A19. For the SLC6A19, we conducted a combined analysis of lncRNA–miRNA–mRNA in the post-transcriptional regulation level by the competing endogenous RNA model [14, 15]. However, we did not discover the co-located and co-expressed lncRNAs with SLC6A19. But there were two miRNAs, up-regulated miR-193b under salt stress and down-regulated miR-542-5p under water-deprivation stress, targeting SLC6A19 mRNA by target gene prediction. In miRNA–mRNA interaction model [16], miR-193b will bind to 3′ untranslated region (3′ UTR) of SLC6A19 mRNA to inhibit the translation of SLC6A19. According to the miRNA-mediated mechanism of extending mRNA half-life [17, 18], the down-regulated expression of miR-542-5p will fail to prevent SLC6A19 mRNA from degradation. Based on previous studies, neutral amino acids (NAAs) are transported by broad neutral (0) amino acid transporter 1 (B0AT1), encoded by SLC6A19, which is a kind of Na+-dependent transporters in the kidney [19, 20]. Glucose transporters (GLUTs) are a large group of membrane proteins. Their essential function is to accelerate the transport of glucose through a plasma membrane [21, 22]. By bioinformatics analysis of screened resistance-related genes, the post-transcriptional regulation pathway under salt and water-deprivation stresses in renal cortex of camel is illustrated (Fig. 4): the inhibition of SLC6A19 mRNA by up-regulated miR-193b and down-regulated miR-542-5p prevents excessive sodium ions from entering the renal cortical cells, and compensates for NAAs by tricarboxylic acid cycling of glucose which enters the cells according to the concentration gradient. PICALM which was also detected in the common 7 mRNAs, can encode phosphatidylinositol-binding clathrin assembly protein and participate in clathrin-mediated endocytosis [23, 24]. The endocytosis of NAAs-contained proteins also assists in compensating NAAs via proteinase.

Fig. 4
figure 4

The post-transcriptional regulation pathway to positively respond to salt and water-deprivation stresses in the renal cortex of camel

The above regulation pathway can be able to reduce excessive sodium ions into the cells of renal cortex, and reduce the amount of glucose which is well soluble in water to improve the amount of intracellular NAAs and purify water since the isoelectric point of NAAs is close to 7 (low solubility in water) [25]. The absorption of NAAs in camel renal cortex may not directly transport free NAAs into the cell through B0AT1, but count on direct endocytosis of the protein and tricarboxylic acid (TCA) cycle of glucose to achieve the purpose of compensating neutral amino acids.


Stress treatment of camel

Nine Alxa Bactrian Camels were randomly divided into three groups of salt stress (experimental group), water-deprivation stress (experimental group), and free diet (control group), and treated for a period of 24 days  (Additional file 2: Figure S1). Each group included three camels. Salt stress group was treated as follows: salt intake base was 200 g/day, and increased by 100 g every 3 days. Consequently, the formula of salt intake corresponding to salt stress days is as follows: an = 200 g + 100 g * int[(n − 1)/3], an refers to the amount of feeding salt, n means salt stress days [26, 27]. Water-deprivation group was free feed intake but fasting water.

Sample collection and preparation

We applied mercy-killing to the nine Alxa Bactrian Camels of three groups by carotid artery letting blood after intramuscular injection of 0.5 mg/kg xylazine [28], then collected renal cortex tissues separately. The renal cortex samples were stored in 1.5 mL frozen tubes and deposited in liquid nitrogen (− 196 °C) for preservation. The samples were sent to Novogene Corporation (Tianjin, China) by freezing in the dry ice for the next procedure.

RNA isolation and quantification

Total RNA was extracted from renal cortex tissues of the nine Alxa Bactrian Camels using RNeasy Mini Kit (QIAGEN, Germany). RNA degradation and contamination were monitored on 1% agarose gels. RNA purity was checked using the NanoPhotometer spectrophotometer (IMPLEN, CA, USA). RNA concentration was measured using Qubit RNA Assay Kit in Qubit 2.0 Fluorometer (Life Technologies, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit with the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA).

Library preparation for small RNA and lncRNA sequencing

A total amount of 3 μg total RNA per sample was used separately as input material for the small RNA and lncRNA library. Sequencing libraries of small RNA were generated using NEBNext Multiplex Small RNA Library Prep Set for Illumina (NEB, USA.). For lncRNA library preparation, firstly, ribosomal RNA was removed by Epicentre Ribo-zero rRNA Removal Kit (Epicentre, USA), and rRNA free residue was cleaned up by ethanol precipitation. Subsequently, sequencing libraries of lncRNA were generated using the rRNA-depleted RNA by NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, USA).

Clustering and sequencing for small RNA and lncRNA

The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq SR Cluster Kit v3-cBot-HS (Illumia) for small RNA and TruSeq PE Cluster Kit v3-cBot-HS (Illumia) for lncRNA. By pooled RNA-seq [29], the library constructions of small RNA and lncRNA were performed on an Illumina HiSeq 2500 platform and Illumina HiSeq 4000 platform, respectively.

Data analysis for miRNA and lncRNA

The analyses procedures of miRNA data were as follows: quality control, reads mapping to Camelus Bactrianus genome by Bowtie [30], novel miRNA prediction by miREvo [31] and mirdeep2 [32], differential expression analysis by TPM [33] and DEGseq [34], small RNA annotation by Bos Taurus database, target gene prediction by miRanda, GO and KEGG pathway enrichment analyses of target genes by GOseq [35] and KOBAS (2.0) [36]. The analyses procedures of lncRNA data were as follows: quality control, reads mapping to C. Bactrianus genome by TopHat (v2.0.9) [37], transcript assembly by Cufflinks [38], lncRNA and mRNA sorting by CNCI (v2) [39], CPC (0.9-r2) [40], Pfam Scan (v1.3) [41, 42], PhyloCSF (v20121028) [43], mRNA and lncRNA annotation by C. Bactrianus and Bos Taurus databases, differential expression analysis by Cuffdiff (, target gene prediction by lncRNA gene upstream/downstream 100 kb and Pearson Correlation Coefficient [44], GO and KEGG pathway enrichment analyses of target genes and differentially expressed mRNAs by GOseq [35] and KOBAS (2.0) [36].


Based on the lncRNA, miRNA, and mRNA data by RNA-seq and incorporated previous studies, we identified four novel salt-resistance-related lncRNAs in renal cortex and proposed the ncRNAs-related post-transcriptional regulation pathway to explain how camels respond to salt stress and water-deprivation stress for the first time, in light of differentially expressed SLC6A19, miR-193b and miR-542-5p. We are hoping to provide a theoretical basis for healing disease that is similar to high blood pressure for humans. Due to the defects in the camel genome function, we selected the database of the bovine which is close to the camel for the annotations of lncRNA and miRNA in this experiment, so we still need more evidences to optimize the results in the future .



broad neutral (0) amino acid transporter 1


competing endogenous RNAs


circular RNAs


glucose transporters


Kyoto Encyclopedia of Genes and Genomes


long non-coding RNA


micro RNA


neutral amino acids


non-coding RNAs


tricarboxylic acid

3′ UTR:

3′ untranslated region


  1. Roberts MBV. Biology: a functional approach. Nelson Thornes; 1986. p. 234–235, 241.

  2. Schmidt-Nielsen K. The physiology of the camel. Sci Am. 1959;201(6):140–51.

    Article  CAS  Google Scholar 

  3. Rehan S, Qureshi AS. Microscopic evaluation of the heart, kidneys and adrenal glands of one-humped camel calves (Camelus dromedarius) using semi automated image analysis system. J Camel Pract Res. 2006;13(2):123–7.

    Google Scholar 

  4. Benlamlih S, Dahlborn K, Filali RZ, Hossaini-Hilali J. Fluid retention after oral loading with water or saline in camels. Am J Physiol. 1992;262(5 Pt 2):R915–20.

    CAS  PubMed  Google Scholar 

  5. Ben Goumi M, Riad F, Giry J, de la Farge F, Safwate A, Davicco MJ, Barlet JP. Hormonal control of water and sodium in plasma and urine of camels during dehydration and rehydration. Gen Comp Endocrinol. 1993;89(3):378–86.

    Article  CAS  Google Scholar 

  6. Geerling JC, Loewy AD. Central regulation of sodium appetite. Exp Physiol. 2008;93(2):177–209.

    Article  CAS  Google Scholar 

  7. Hill RW, Wyse GA, Anderson M. Animal Physiology. Sunderland: Sinauer Associates, Inc; 2016. p. 792.

    Google Scholar 

  8. Alpern RJ, Moe OW, Caplan M, et al. Seldin and Giebisch’s the kidney: physiology and pathophysiology. Cambridge: Academic Press; 2013. p. 630–3.

    Google Scholar 

  9. Lee CJ, Gardiner BS, Ngo JP, et al. Accounting for oxygen in the renal cortex: a computational study of factors that predispose the cortex to hypoxia. Am J Physiol Renal Physiol. 2017;313(2):F218–36.

    Article  CAS  Google Scholar 

  10. Wu H, Guang X, Al-Fageeh MB, et al. Camelid genomes reveal evolution and adaptation to desert environments. Nat Commun. 2014;5:5188.

    Article  CAS  Google Scholar 

  11. Jeffares DC, Poole AM, Penny D. Relics from the RNA world. J Mol Evol. 1998;46(1):18–36.

    Article  CAS  Google Scholar 

  12. Poole AM, Jeffares DC, Penny D. The path from the RNA world. J Mol Evol. 1998;46(1):1–17.

    Article  CAS  Google Scholar 

  13. Poole A, Jeffares D, Penny D. Early evolution: prokaryotes, the new kids on the block. BioEssays. 1999;21(10):880–9.

    Article  CAS  Google Scholar 

  14. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Cell. 2011;146(3):353–8.

    Article  CAS  Google Scholar 

  15. Hansen TB, Jensen TI, Clausen BH, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.

    Article  CAS  Google Scholar 

  16. Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.

    Article  CAS  Google Scholar 

  17. Vasudevan S, Steitz JA. AU-rich-element-mediated upregulation of translation by FXR1 and Argonaute 2. Cell. 2007;128(6):1105–18.

    Article  CAS  Google Scholar 

  18. Vasudevan S, Tong Y, Steitz JA. Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007;318(5858):1931–4.

    Article  CAS  Google Scholar 

  19. Kleta R, Romeo E, Ristic Z, Ohura T, Stuart C, Arcos-Burgos M, Dave MH, Wagner CA, Camargo SR, Inoue S, Matsuura N, et al. Mutations in SLC6A19, encoding B0AT1, cause Hartnup disorder. Nat Genet. 2004;36(9):999–1002.

    Article  CAS  Google Scholar 

  20. Broer A, Klingel K, Kowalczuk S, Rasko JE, Cavanaugh J, Broer S. Molecular cloning of mouse amino acid transport system B0, a neutral amino acid transporter related to Hartnup disorder. J Biol Chem. 2004;279(23):24467–76.

    Article  CAS  Google Scholar 

  21. Fukumoto H, Seino S, Imura H, Seino Y, Bell GI. Characterization and expression of human HepG2/erythrocyte glucose-transporter gene. Diabetes. 1988;37(5):657–61.

    Article  CAS  Google Scholar 

  22. Kasahara M, Inui K, Takano M, Hori R. Distinction of three types of d-glucose transport systems in animal cells. Biochem Biophys Res Commun. 1985;132(2):490–6.

    Article  CAS  Google Scholar 

  23. Ando K, Tomimura K, Sazdovitch V, Suain V, Yilmaz Z, et al. Level of PICALM, a key component of clathrin-mediated endocytosis, is correlated with levels of phosphotau and autophagy-related proteins and is associated with tau inclusions in AD, PSP and Pick disease. Neurobiol Dis. 2016;94:32–43.

    Article  CAS  Google Scholar 

  24. Thomas RS, Henson A, Gerrish A, Jones L, et al. Decreasing the expression of PICALM reduces endocytosis and the activity of β-secretase: implications for Alzheimer’s disease. BMC Neurosci. 2016;17(1):50.

    Article  Google Scholar 

  25. Dayton WR. Protein separation techniques. Reciprocal Meat Conf Proc. 1983;36:98–102.

    Google Scholar 

  26. Assad F, El-Sherif MMA. Effect of drinking saline water and feed shortage on adaptive responses of sheep and camels. Small Ruminant Res. 2002;45(3):279–90.

    Article  Google Scholar 

  27. Leitch I. The feeding of camels. In: Imperial bureau of animal nutrition technical communication No. 13; 1940. p. 10.

  28. Custer R, Kramer L, Kennedy S, Bush M. Hematologic effects of xylazine when used for restraint of Bactrian camels. J Am Vet Med Assoc. 1977;171(9):899–901.

    CAS  PubMed  Google Scholar 

  29. Konczal M, Koteja P, Stuglik MT, et al. Accuracy of allele frequency estimation using pooled RNA-Seq. Mol Ecol Resour. 2014;14(2):381–92.

    Article  CAS  Google Scholar 

  30. Langmead B, Trapnell C, Pop M, et al. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.

    Article  Google Scholar 

  31. Wen M, Shen Y, Shi S, et al. miREvo: an integrative microRNA evolutionary analysis platform for next-generation sequencing experiments. BMC Bioinform. 2012;13:140.

    Article  CAS  Google Scholar 

  32. Friedländer MR, Mackowiak SD, Li N, et al. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.

    Article  Google Scholar 

  33. Zhou L, Chen J, Li Z, et al. Integrated profiling of MicroRNAs and mRNAs: MicroRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PLoS ONE. 2010;5(12):e15224.

    Article  CAS  Google Scholar 

  34. Wang L, Feng Z, Wang X, et al. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.

    Article  Google Scholar 

  35. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.

    Article  Google Scholar 

  36. Mao X, Cai T, Olyarchuk JG, et al. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787–93.

    Article  CAS  Google Scholar 

  37. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.

    Article  Google Scholar 

  38. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.

    Article  CAS  Google Scholar 

  39. Liang S, Luo H, Bu D, et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013;41(17):e166.

    Article  Google Scholar 

  40. Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007;35(suppl_2):W345–9.

    Article  Google Scholar 

  41. Punta M, Coggill PC, Eberhardt RY, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40(D1):D290–301.

    Article  CAS  Google Scholar 

  42. Bateman A, Birney E, Cerruti L, et al. The Pfam protein families database. Nucleic Acids Res. 2002;30(1):276–80.

    Article  CAS  Google Scholar 

  43. Lin MF, Jungreis I, Kellis M. PhyloCSF: a comparative genomics method to distinguish protein coding and non-coding regions. Bioinformatics. 2011;27(13):i275–82.

    Article  CAS  Google Scholar 

  44. Ma L, Bajic VB, Zhang Z. On the classification of long non-coding RNAs. RNA Biol. 2013;10(6):924–33.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

HZ and DZ conceived the study and designed the experiments. DZ and YC participated in sampling and performed the experiments. YC analyzed the data and wrote the manuscript. HZ and DZ provided valuable insights in the discussion and revision of the manuscript. All authors read and approved the final manuscript.


We thank Z. Dong, X. Liu, Y. Ling and J. Pan for camel’s feeding and experimental treatment. We are grateful for sampling assistance of F. Meng, Y. Liu, L. Li, T. Li, S. Wang, J. Pan et al.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analyzed during this study is included in this article.

Consent to publish

Not applicable.

Ethics approval and consent to participate

Procedures involving animals, their care and humane kill were conducted in conformity with Guidelines on the Humane Treatment of Laboratory Animals (HTLA Pub. Chapter 2-6, revised 2006 in China) and was approved by Animal Care and Use Committee of the Inner Mongolia Agricultural University. The informed consent of experimental research was achieved by Shenyuan Wang (as animal owner).


This work was supported by National Natural Science Foundation of China (31560313 and 31260271).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Dong Zhang or Huanmin Zhou.

Additional files

Additional file 1: Table S1.

The differentially expressed mRNAs, miRNAs and lncRNAs of renal cortex under salt stress. Table S2. The differentially expressed mRNAs and miRNAs of renal cortex under water-deprivation stress. Table S3. Sequence of significantly down-regulated four novel lncRNAs in the renal cortex of camel under salt stress.

Additional file 2: Figure S1.

Normalized mean plasma Na+ concentration of camel under salt stress, water-deprivation stress and free diet.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cao, Y., Zhang, D. & Zhou, H. Key genes differential expressions and pathway involved in salt and water-deprivation stresses for renal cortex in camel. BMC Molecular Biol 20, 11 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: