- Research article
- Open Access
Key genes differential expressions and pathway involved in salt and water-deprivation stresses for renal cortex in camel
BMC Molecular Biology volume 20, Article number: 11 (2019)
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 , 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 . 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 . In kidney, the formation and excretion of urine are involved in renal cortex and then renal medulla . 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 . 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).
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.
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.
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 , 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.
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) . 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 , 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 , 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 , novel miRNA prediction by miREvo  and mirdeep2 , differential expression analysis by TPM  and DEGseq , small RNA annotation by Bos Taurus database, target gene prediction by miRanda, GO and KEGG pathway enrichment analyses of target genes by GOseq  and KOBAS (2.0) . The analyses procedures of lncRNA data were as follows: quality control, reads mapping to C. Bactrianus genome by TopHat (v2.0.9) , transcript assembly by Cufflinks , lncRNA and mRNA sorting by CNCI (v2) , CPC (0.9-r2) , Pfam Scan (v1.3) [41, 42], PhyloCSF (v20121028) , mRNA and lncRNA annotation by C. Bactrianus and Bos Taurus databases, differential expression analysis by Cuffdiff (http://cole-trapnell-lab.github.io/cufflinks/cuffdiff/index.html), target gene prediction by lncRNA gene upstream/downstream 100 kb and Pearson Correlation Coefficient , GO and KEGG pathway enrichment analyses of target genes and differentially expressed mRNAs by GOseq  and KOBAS (2.0) .
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
Kyoto Encyclopedia of Genes and Genomes
long non-coding RNA
neutral amino acids
- 3′ UTR:
3′ untranslated region
Roberts MBV. Biology: a functional approach. Nelson Thornes; 1986. p. 234–235, 241.
Schmidt-Nielsen K. The physiology of the camel. Sci Am. 1959;201(6):140–51.
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.
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.
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.
Geerling JC, Loewy AD. Central regulation of sodium appetite. Exp Physiol. 2008;93(2):177–209.
Hill RW, Wyse GA, Anderson M. Animal Physiology. Sunderland: Sinauer Associates, Inc; 2016. p. 792.
Alpern RJ, Moe OW, Caplan M, et al. Seldin and Giebisch’s the kidney: physiology and pathophysiology. Cambridge: Academic Press; 2013. p. 630–3.
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.
Wu H, Guang X, Al-Fageeh MB, et al. Camelid genomes reveal evolution and adaptation to desert environments. Nat Commun. 2014;5:5188.
Jeffares DC, Poole AM, Penny D. Relics from the RNA world. J Mol Evol. 1998;46(1):18–36.
Poole AM, Jeffares DC, Penny D. The path from the RNA world. J Mol Evol. 1998;46(1):1–17.
Poole A, Jeffares D, Penny D. Early evolution: prokaryotes, the new kids on the block. BioEssays. 1999;21(10):880–9.
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.
Hansen TB, Jensen TI, Clausen BH, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.
Vasudevan S, Steitz JA. AU-rich-element-mediated upregulation of translation by FXR1 and Argonaute 2. Cell. 2007;128(6):1105–18.
Vasudevan S, Tong Y, Steitz JA. Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007;318(5858):1931–4.
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.
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.
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.
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.
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.
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.
Dayton WR. Protein separation techniques. Reciprocal Meat Conf Proc. 1983;36:98–102.
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.
Leitch I. The feeding of camels. In: Imperial bureau of animal nutrition technical communication No. 13; 1940. p. 10.
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.
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.
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.
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.
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.
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.
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.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.
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.
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.
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.
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.
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.
Punta M, Coggill PC, Eberhardt RY, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40(D1):D290–301.
Bateman A, Birney E, Cerruti L, et al. The Pfam protein families database. Nucleic Acids Res. 2002;30(1):276–80.
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.
Ma L, Bajic VB, Zhang Z. On the classification of long non-coding RNAs. RNA Biol. 2013;10(6):924–33.
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.
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
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).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.