Skip to main content

ORIGINAL RESEARCH article

Front. Vet. Sci., 14 April 2022
Sec. Veterinary Pharmacology and Toxicology
This article is part of the Research Topic Environmental Contaminants and Animal Health: Analysis, Toxicity, and Mitigation View all 8 articles

Whole Transcriptome Profiling of the Effects of Cadmium on the Liver of the Xiangxi Yellow Heifer

Updated
\nYameng Wei&#x;Yameng Wei1Kangle Yi&#x;Kangle Yi2Caomeihui ShenCaomeihui Shen1Xue ChenXue Chen1Tariq IqbalTariq Iqbal1Maosheng CaoMaosheng Cao1Tong ChenTong Chen1Yang LuoYang Luo2Jianbo LiJianbo Li2Xu ZhouXu Zhou1Chunjin LiChunjin Li1Lu Chen
Lu Chen1*
  • 1College of Animal Sciences, Jilin University, Changchun, China
  • 2Grassland and Herbivore Research Laboratory, Hunan Animal Husbandry and Veterinary Research Institute, Changsha, China

Cadmium (Cd) is a major heavy metal toxicant found in industrial zones. Humans and animals are exposed to it through their diet, which results in various physiological problems. In the current study, the toxic effects of Cd on the liver were investigated by whole-transcriptome sequencing (RNA-seq) of the livers of Xiangxi heifers fed a diet with excess Cd. We randomly divided six healthy heifers into two groups. The first group received a control diet, whereas the second group received Cd-exceeding diets for 100 days. After 100 days, the livers were collected. A total of 551 differentially expressed mRNAs, 24 differentially expressed miRNAs, and 169 differentially expressed lncRNAs were identified (p < 0.05, |log2FC| >1). Differentially expressed genes (DEGs) were analyzed by gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. We found that under Cd exposure, DEGs were enriched in the adenosine 5'-monophosphate–activated protein kinase pathway, which is involved in autophagy regulation, and the peroxisome proliferator–activated receptor pathway, which is involved in lipid metabolism. In addition, the apolipoprotein A4 gene, which has anti-inflammatory and antioxidant effects, the anti-apoptotic gene ATPase H+/K+ transporting the nongastric alpha2 subunit, and the cholesterol metabolism–associated gene endothelial lipase gene were significantly downregulated. C–X–C motif chemokine ligand 3, cholesterol 7α-hydroxylase, and stearoyl-CoA desaturase, which are involved in the development of fatty liver, were significantly upregulated. These genes revealed the main effects of Cd on the liver of Xiangxi yellow heifers. The current study provides insightful information regarding the DEGs involved in autophagy regulation, apoptosis, lipid metabolism, anti-inflammation, and antioxidant enzyme activity. These may serve as useful biomarkers for predicting and treating Cd-related diseases in the future.

Introduction

Heavy metals accumulate significantly in the soils of most parts of China, with cadmium (Cd) being the most prevalent heavy metal contaminant. Contaminated areas are mainly located in the central, southern, and southwestern parts of China, and Hunan is one of the most seriously contaminated provinces (1). There are several large mines and smelting areas in the Hunan Province, which is known for its non–ferrous metal mining activities and is a typical area of dense low-temperature deposits in China (2). Heavy metal contamination in soil is considered a potential cancer-causing health risk for Hunan residents. Moreover, heavy metal pollution in rivers used for agricultural irrigation in Hunan Province is considered to indirectly affect the health of humans and animals (3). Garden soil, rice soil, vegetables, and rice in the Hunan mining area are contaminated with heavy metals (4). Cd pollution is considered to be the most serious heavy metal pollution in Hunan. The Cd content in both rice straw and endosperm samples in Hunan is higher than that in the Yangtze River delta (5). The western part of Hunan is known as Xiangxi. The Xiangxi yellow heifer is the main farmed animal in the Xiangxi region, which mainly consumes local forage. This is highly significant for the development of the Xiangxi farming industry. The beneficial characteristics of Xiangxi yellow heifers include mild temperament, resistance to heat, medium size, well-developed and strong bones, tender meat, and rich nutrition (6). Beef is one of the main food sources for the local people. Studies have shown that metal accumulation in livestock can be transferred to humans through the food chain (7, 8).

Heavy metals greatly endanger the health of animals and humans through environmental exposure and the food chain. Cd contamination has caused several vicious incidents, including the Cd rice incident in Hunan, China, and the Itai-Itai disease public health incident in Japan. In China, Hunan, Henan, Anhui, and Jilin provinces have high levels of Cd in the soil. Cd in Hunan farmland exceeds the soil environmental background value by 30-fold (9). Cd is absorbed mainly by the respiratory and digestive tracts. Moreover, Cd adversely affects the liver, kidneys, lungs, immune system, reproductive system, and cardiovascular system (10). When Cd-containing foods and water are mistakenly consumed, Cd is transported from the intestine to the liver, causing functional disorders (11). Studies have shown that Cd interferes with collagen metabolism and extracellular signal transduction in the liver matrix, as well as cell adhesion, growth, and migration (1214). Cd also causes morphological and structural changes in the liver, such as blurred hepatic trabeculae, increased cell size, hepatocyte deformation, nuclear consolidation, lipid droplet aggregation, and mono-nuclear cell infiltration (15). Hepatic mitochondrial function has been reported to be impaired in a rat model of Cd exposure (1618). Cd exposure inhibits mitochondrial oxidative phosphorylation, resulting in a decrease in ATP. In addition, high concentrations of Cd inhibit basal respiration (19). Moreover, Cd can directly damage the mitochondrial structures in the liver, causing mitochondrial degeneration, vacuolation, and disruption of the mitochondrial membrane structure (20, 21). In the presence of severe mitochondrial damage, the process of mitochondrial autophagy is activated. Metallothionein (MT), which binds to Cd, protects organisms from Cd toxicity. Cd decreases the levels of MT and antioxidant glutathione (GSH) and inhibits the activity of antioxidant enzymes (22). Therefore, it can induce oxidative stress (23, 24). Cd can also induce systemic inflammatory responses involving innate, adaptive, and mucosal immunity (25). In addition, the increase in Ca2+ levels accelerates reactive oxygen species (ROS) formation, also resulting in oxidative stress (26, 27). Cd induces endoplasmic reticulum stress and DNA mutations by inhibiting the DNA repair system and inducing protein damage (2830). Cd damages the adenosine triphosphate system and the tricarboxylic acid cycle, which can lead to an increase in various metabolites related to lipid metabolism and triglycerides and disrupted energy metabolism (18). Cd exposure is a potential risk factor for the development of diabetes mellitus, non-alcoholic steatohepatitis (NASH), and non-alcoholic fatty liver disease (NAFLD). Cd exacerbates hepatocyte lipid storage through dysregulation of autophagy, which aggravates hepatic steatosis and leads to hyperlipidemia (31). Cd is also associated with an increased risk of prostate, breast, lung, and hepatocellular carcinoma (3234).

RNA-seq, characterized by low background noise and a large detection range, can provide a large amount of information on gene expression levels, facilitating the study of genetic mechanisms. Currently, it is the preferred method for studying gene expression (35, 36). Some studies have claimed that heifers are potential biological carriers of heavy metal contamination in soil (8). An association exists between metal levels in heifer tissues and the risk of human exposure. In addition, the liver is the most prominent target organ of Cd toxicity (37).

Therefore, in this study, the Hunan landmark product Xiangxi yellow heifers were used as experimental animals. Then, RNA-seq was used to study the effects of feeding Cd overload diets on the expression of mRNAs, miRNAs, and lncRNAs in the liver of Xiangxi yellow heifers to understand their transcriptional biology. At present, there are numerous studies on the hepatotoxicity of Cd, but there is a relative lack of information on the hepatic transcriptomics of Cd-exposed Chinese Xiangxi yellow heifers. Hunan Province is rich in mineral resources. However, Cd pollution is present. Cd has great potential to endanger the health of the local dominant breed of Xiangxi yellow heifers. Therefore, this study hypothesized that certain transcription factors and metabolic pathways would be affected by Cd exposure during foraging in Chinese Xiangxi yellow heifers. The current experiment aimed to find the pathological changes that occurred in the liver of Xiangxi yellow cattle after feeding Cd-containing forage. This experiment also aimed to analyze the gene expression at the RNA level and the pathways that are significantly enriched in Xiangxi yellow cattle exposed to cadmium using transcriptomics. Our findings provide new insights into the potential health risks of Cd, improve the current understanding of the toxic effects of Cd on the liver, and help to prevent Cd toxicity.

Materials and Methods

Materials, Animals, and Treatment

The experiments were conducted at the National Origin Farm of Xiangxi Cattle, which is located 10 km south of Huayuan County, Hunan, China. All experimental protocols and procedures used in this study followed Chinese standards for the use and care of research animals (SOP-017).

A total of six healthy female crossbred (Angus bull × Xiangxi local cattle) heifers weaned (BW = 78.50 ± 18.45 kg) at 6 months of age were selected from a paternal halfsib family. They were randomly allocated to the control group (control) or cadmium-treated group (cadmium), with three heifers in each group. Heifers were housed in an open-sided, naturally well-ventilated barn, and the groups were separated using railings. Thus, three heifers within each group were housed together in one pen. The feeding experiments were conducted for 109 days. The first 9 days were the pre-feeding period, during which the test heifers were dewormed and acclimatized to the diet. The following 100 days were the regular experimental feeding periods. During the experiments, the diets of both groups consisted of concentrate, paddy, and corn straw silage (Table 1). These diets were mixed daily and provided to the cattle a total mixed ration. In China, the standard for agricultural soil cadmium is 0.3 mg/kg (38). The corn silage used for the cadmium-treated group was made from a whole-plant corn species with strong cadmium enrichment ability (Yuqing 23). This corn was planted in soil with excessive cadmium content (3–5 mg/kg) in Yonghe Town, Liuyang City, Hunan Province. The corn silage used for the control group was made from normal whole-plant corn that was planted in soil with a Cd content of 0.2 mg/kg. The Cd concentrations of the diets in the control and cadmium groups were 0.72 and 6.74 mg/kg, respectively. The orts were removed one time daily before the morning feed, and a new feed was delivered two times daily at 7:00 a.m. and 17:00 a.m. Cattle within each group shared one feed trough. All the cattle had free access to water. A part of the liver was placed in 10% formalin and used to prepare histological sections for morphological analysis. The remaining liver tissue samples were frozen at a −80°C in a refrigerator for RNA-seq.

TABLE 1
www.frontiersin.org

Table 1. Feed formulation for the control group and cadmium-treated group.

Preparation of Liver Tissue Sections

Liver tissues were fixed in 10% formalin for 24–48 h, followed by dehydration in different grades of ethanol. Next, 6-μm-thick sections were cut from the prepared liver blocks, affixed to glass slides, rehydrated, and stained with hematoxylin and eosin. The prepared slides were observed under a microscope equipped with a microphotographic system.

Measuring Cd Concentration in Liver Tissues Using Atomic Absorption Spectrophotometry

Different concentrations (0, 0.1, 0.2, 0.4, 0.6, 1.0, 2.0, 3.0, and 4.0 μg/mL) of HNO3 solution were prepared by the National Center for Analysis and Testing of Non-ferrous Metals and Electronic Materials.

First, 0.5 g of the liver was dried in a desiccator at 75°C for 1 h. Then, 3 mL of the digestion solution (HNO3:HCLO = 7:3) was added and heated to near evaporation. The residue was dissolved in 1.0% HNO3. The samples to be measured and cadmium standard solutions with different standard concentrations were placed in the Z-5000 graphite furnace atomic absorption spectrophotometer. Finally, the determination was carried out separately according to a pre-set procedure (drying at 80–140°C for 40 s; ashing at 300°C for 20 s; and atomization at 1,500°C for 5 s). Finally, a concentration–absorbance curve was automatically generated by the machine.

Library Preparation and Sequencing

Liver tissue samples frozen at −80°C were used for RNA-seq. A 50–100 mg sample was collected from each group and added to 0.5 ml TRIzol (Invitrogen, Carlsbad, CA, United States) on ice. Then, the samples were mixed using a homogenizer (Cebo, Shanghai, China), and 0.5 mL TRIzol was added. The tissue and TRIzol were transferred together into a 1.5 mL EP tube and left at 37°C for 5 min. The sample was centrifuged at 12,000 × g for 5 min and the sediment was discarded. Then, 0.2 mL of chloroform was added. The samples were shaken for 15 s and incubated at room temperature for 3 s. The samples were then centrifuged at 4°C for 20 min at 12,000 × g. When the samples appeared to be stratified, the RNA was mainly in the upper aqueous phase; therefore, the aqueous phase was transferred to a new tube. Immediately after, 100% isopropanol (0.5 mL) was added. The samples were mixed by inversion and left at 37°C for 10 min. The supernatant was removed by centrifugation at 12,000 × g for 10 min at 4°C. Then, 1 mL of 80% ethanol was added and the sample was mixed. After centrifugation at 7,500 × g for 5 min at 4°C, the supernatant was discarded. The samples were then left at room temperature for approximately 5–10 min. Finally, 50 μL of DEPC water was added to solubilize the RNA, and the RNA solution was obtained by centrifugation under a light flow. Total RNA concentration and quality were determined using a NanoDrop spectrophotometer and an Agilent 2,100 Bioanalyzer (Thermo Fisher Scientific, MA, United States).

Two new libraries, lncRNA/mRNA and miRNA, were constructed. In the lncRNA/mRNA library, the Ribo-Zero Magnetic Kit (Epicenter) was used to remove 5.8S rRNA, 18S rRNA, 28S rRNA, 12S rRNA, and 16S rRNA from total RNA. Single-stranded cDNA was synthesized using random primers and reverse transcriptase. DNA polymerase I and RNase H (Invitrogen) were used to synthesize the second-strand cDNA. At this point, the RNA template was removed and dTTP was replaced with dUTP. The double-stranded cDNA product was then subjected to 3'-adenylation and splice ligation. The products were subjected to several rounds of PCR and thermal denaturation to a single strand and then cyclized by splint oligo to obtain single-stranded circular DNA. Paired-end sequencing (100 bp) was performed using a BGISEQ-500/MGISEQ-2000 system (BGI-Shenzhen, China). Each sample from the lncRNA/mRNA library yielded an average of 13.65 G of data.

For miRNA sequencing, 18–30 nt of small RNA (Ladder Marker, Takara) was isolated and purified from total RNA by polyacrylamide gel electrophoresis (PAGE) gel cutting. The 5-adenylated, 3-blocked single-stranded DNA junction was ligated to the 3' end of the purified RNA, and the 5' junction was ligated to the 5' end of the RNA. RNA with 3' and 5' junctions attached was synthesized into cDNA by reverse transcription extension using RT primers bearing UMI tags. The cDNA was subsequently amplified by PCR. The PCR products were separated by PAGE to obtain fragments in the range of 110–130 bp. Fragments were purified using a QIAquick Gel Extraction Kit (QIAGEN, CA, United States) and quality-controlled using an Agilent 2,100 Bioanalyzer (Thermo Fisher Scientific, United States). The purified library products were subsequently sequenced at the 50 bp single-end on the DNBseq platform (BGI-Shenzhen, China). The raw data obtained from sequencing were converted into raw sequence data (raw data or raw reads) by base calling and stored in the fastq file format, which contains the sequence of reads and sequencing quality information. The miRNA library yielded an average of 23.77 M data per sample. miRanda (3.3a) and TargetScan (7.1) were used to predict the target genes of the miRNAs, and the David database (6.8) was used for pathway enrichment analysis of the miRNA target genes.

Transcriptome Data Analysis

The sequencing data in the fastq format were filtered using SOAPnuke (v1.5.2) (39) for reads with unknown base N content >5%, reads with splice contamination, and reads without insert tags. The clean reads obtained were then mapped to the reference genome Bos taurus (vGCF_000003205.7_Btau_5.0.1) using Bowtie2 (v2.2.5) (40). mRNA expression was calculated using RSEM (v1. 2.12), which uses the normalization method fragments per kilobase of exon model per million mapped fragments (FPKM). miRNA expression was calculated using reads of exon model per million mapped reads values. Finally, the DESeq2 (41) method was used to determine the differential expression of two groups of genes (three biological replicates per group) based on a negative binomial distribution model with high sensitivity and accuracy. The obtained P was calculated using the likelihood-ratio test to represent the difference between the two groups of samples at the inclusion level. For DEG analysis, a P ≤ 0.05 and | log2FC |>1 were considered significant.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Enrichment Analyses

To annotate the biological functions of the genes, all differentially expressed mRNAs (DEmRNAs) were aligned against the KEGG and GO databases. The number of genes per GO term and the enriched KEGG pathways were calculated. The candidate genes were then compared to all background genes of the species using the phyper function (https://en.wikipedia.org/wiki/Hypergeometric_distribution) based on a hypergeometric test to obtain the p. The p-values were corrected using the Benjamini–Hochberg algorithm to obtain FDR values. FDR ≤ 0.05 was used as the threshold value. Finally, the GO terms and KEGG pathways that met this condition were recognized as significantly enriched compared with the background genes of this species.

In the histogram of GO enrichment, the X-axis represents the number of genes annotated to a particular term. The Y-axis represents the number of selected genes annotated to a particular term. In the KEGG pathway diagram, the rich factor on the X-axis represents the ratio of the number of DEGs in the metabolic pathway to the number of genes annotated in the pathway, with larger values indicating greater enrichment. The Y-axis represents the pathway to which the gene was annotated. Box plots and variable shear plots were constructed to determine the quality of the RNA-seq data. The box plot demonstrates the distribution of gene expression levels for each sample and displays the dispersion of the data distribution, which was generated using the OmicStudio tools (https://www.omicstudio.cn/tool). rMATS (v3.2.5) (42) was used to detect splicing events in the samples.

RNA Extraction

Liver tissues of Xiangxi yellow heifers were removed from the −80°C refrigerator. Total RNA was extracted according to the instructions. A small portion of the liver was homogenized in 1 mL of TRIzol. Then, 200 μL of chloroform was added and the sample was shaken vigorously. The sample was centrifuged at 1,200 × g for 3 min and the supernatant was transferred to an absorption column. Next, 400 μL of 96% ethanol was added and recentrifuged at 12,000 × g for 1 min. The supernatant was removed from the absorption column and 500 μL of 3 M sodium acetate trihydrate solution was added. It was then subjected to centrifugation at 12,000 × g for 1 min. Finally, 70% alcohol (400 μL) was added and the mixture was centrifuged at 12,000 × g in triplicate. The supernatant was removed from the absorption column, and 46–50 μL of DEPC was added to the RNA pellet. Total RNA was collected in Eppendorf tubes and centrifuged for 1 min at 10,000 × g. Total RNA content was measured using a NanoDrop spectrophotometer and an Agilent 2,100 Bioanalyzer (Thermo Fisher Scientific, MA, United States).

Quantitative RT–PCR

cDNA was synthesized by reverse transcription using the Superscript III reverse transcription kit (ABI Invitrogen, United States). qRT-PCR was performed for genes associated with Cd exposure using a fluorescent quantitative PCR instrument (Applied Biosystems, United States). Specific quantitative primers (KuMei, Changchun, China) for the 18 transcripts are shown in Supplementary Table S1. The total volume of the qRT-PCR mixture was 20 μL, which contained 10 μL of SYBR (Takara Bio), 1 μL of the forward primer, 1 μL of the reverse primer, 2 μL of the cDNA, and 6 μL of ddH2O. The reaction conditions were as follows: 1 cycle at 95°C for 5 min, 40 cycles at 95°C for 10 s, 58°C for 20 s, 72°C for 20 s, and 1 cycle at 95°C for 15 s. The magnitude of the fold differences in gene expression levels was calculated using the 2−ΔΔCT method. An independent-samples t-test was used to test the difference between the cadmium-treated and control group. Statistical significance was set at p < 0.05.

Statistical Analysis

t-Tests were performed on all experimental data in triplicate using IBM SPSS Statistics 22 and GraphPad Prism 5. All results are expressed as mean ± SEM, and the significance level was set at * p < 0.05, ** p < 0.01, and *** p < 0.001.

Results

Cadmium Accumulation and Pathological Changes in the Liver

The Cd content in the liver of heifers in the Cd-treated group was significantly higher than that in the control group. This finding indicates that Cd accumulated in the liver of Xiangxi yellow heifers fed a high-Cd diet for 100 days (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. The content of cadmium in liver tissues. Cadmium content in the liver tissues of the control group (Control) and cadmium-treated group (Cadmium), as obtained by graphite furnace atomic absorption spectrophotometry (μg/kg). Data were analyzed by independent samples t-test. The significance level was set at * p < 0.05.

Histopathological observations reveal that the liver of heifers in the cadmium-treated group had expanded sinusoids compared to the control group (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. H and E staining. Pathological changes in the liver of control heifers (Control) and cadmium-treated heifers (Cadmium) by H and E staining (40×, 100×). Arrows indicate expanded hepatic sinusoids.

Analysis of RNA-seq Libraries

To evaluate the toxic effects of Cd on the liver of Xiangxi yellow heifer, the liver was collected for whole transcriptome analysis. As shown in Table 2, a total of 132–140 million raw reads were obtained. After removing reads with an unknown base N content >5% and low-quality reads, 131–141 million clean reads remained, and 13–14 Gb of clean bases were acquired.

TABLE 2
www.frontiersin.org

Table 2. Data filtering for RNA-seq.

Clean reads were mapped to the reference genome. The localization rate was ~91.8% (Table 3). Approximately 88.6% of the reads could be uniquely mapped to the Bos taurus genome. The box plot represents the degree of dispersion and volatility of the data distribution (Figure 3A). Figure 3A shows that the boxes were symmetrical about the middle line, indicating that the data were normally distributed. The smaller height of the boxes was reflected by the smaller degree of fluctuation in the data. Good standardization of the data was indicated by the consistency of the data distribution of the samples in Figure 3A. As shown in Figure 3, mRNA was mainly generated by exon skipping to produce splice isoforms of different mRNAs and regulate gene expression.

TABLE 3
www.frontiersin.org

Table 3. Reference genome mapping.

FIGURE 3
www.frontiersin.org

Figure 3. Box plots and variable shear events. (A) Box plots of the gene expression levels: the X-axis is the sample name, the Y-axis is log10 (FPKM+1), and each region of the box plot corresponds to five statistics (upper limit, upper quartile, median, lower quartile, and lower limit, respectively, from top to bottom, where upper and lower limits are not considered outliers). (B) Variable Shear Events: alternative 5' splicing site (A5SS), alternative 3' splicing site (A3SS), mutually exclusive exons (MXE), retained intron (RI), skipped exon (SE).

DEGs Analysis

Cluster analysis of all DEGs was performed to observe gene expression patterns. Heat maps were generated, with red and blue representing upregulated and downregulated genes, respectively (Figure 4). The overall expression pattern of genes in the Cd-treated group was quite different from that of the control group, and the three samples within the same group showed similar expression patterns. This indicates that the samples from both groups were reproducible, demonstrating that the derived data were credible.

FIGURE 4
www.frontiersin.org

Figure 4. Heatmap of DEmRNAs, DElncRNAs, and DEmiRNAs clustering analysis. Red represents upregulated genes, and blue represents downregulated genes.

Volcano plots of the DEGs showed the differences in the distribution of gene expression levels between control and Cd-treated samples. A total of 551 DEmRNAs were identified in the control and cadmium-treated bovine livers. Of these, 318 genes were upregulated, and 233 genes were downregulated. Compared with the control group, 24 differentially expressed miRNAs (DEmiRNAs) were identified. Of these, 13 genes were upregulated and 11 downregulated in the cadmium-treated group. In total, 169 differentially expressed lncRNAs (DElncRNAs), were identified. Of these, 73 genes were upregulated and 96 genes downregulated (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. Volcano plot of DEmRNAs, DElncRNAs, and DEmiRNAs. Red represents upregulated genes, and blue represents downregulated genes.

Functional Analyses of DEGs

A total of 551 significant DEGs in the liver, were subjected to GO and KEGG pathway analyses to determine the enrichment of functions and action pathways by DEGs. Based on the results of the functional classification of DEGs in the categories of biological process (BP), cellular component (CC), and molecular function (MF), 92 GO terms were obtained. The terms of interest were included in the first 30 enriched GO terms (including oxidation–reduction process, cholesterol biosynthetic process, cell cycle, cyclin B1–CDK1 complex, cyclin-dependent protein kinase holoenzyme complex, and oxidoreductase activity) (Figure 6A). Table 4 lists the DEmRNAs involved in these functionally related events.

FIGURE 6
www.frontiersin.org

Figure 6. GO and KEGG pathway analyses. (A) Histograms of GO enrichment, KEGG functional classification of DEmRNAs. (B) Significantly enriched KEGG pathways (padj <0.05).

TABLE 4
www.frontiersin.org

Table 4. DEGs included in the GO terms.

The KEGG pathways that were significantly enriched (padj <0.05) were metabolic, peroxisome proliferator–activated receptor (PPAR) signaling, AMPK signaling, and steroid biosynthesis pathways (Figure 6B; Table 5). Table 6 shows the DEmRNAs involved in these signaling pathways in Cd-exposed livers of Xiangxi yellow heifers and the positions of these DEmRNAs in the pathways (Figure 7).

TABLE 5
www.frontiersin.org

Table 5. Results of the KEGG enrichment analysis of DEGs.

TABLE 6
www.frontiersin.org

Table 6. DEGs included in the PPAR signaling pathway, metabolic pathways, and AMPK signaling pathway.

FIGURE 7
www.frontiersin.org

Figure 7. The location of the DEGmRNAs in the significantly enriched pathways. Diagram of the PPAR signaling pathway (A) and AMPK signaling pathway (B). The red pentagram represents DEGmRNAs.

qPCR Analysis of Liver DEGs in Response to Cd

Real-time fluorescence quantitative PCR was used to validate the 11 DEmRNAs, four DEmiRNAs, and four DElncRNAs selected for this study. Among the DEmRNAs, stearoyl-CoA desaturase (SCD), voltage-gated channel auxiliary subunit gamma 4 (CACNG4), CYP7A1, insulin-like growth factor 1 (IGF1), and C–X–C motif chemokine ligand 3 (CXCL3) were significantly upregulated. The difference in CXCL3 expression was extremely significant. In addition, LIPG, ATPase H+/K+ transporting the non-gastric alpha2 subunit (ATP12A), apolipoprotein A (ApoA4), bone morphogenetic protein 8A (BMP8A), and bone morphogenetic protein 8B (BMP8B) were significantly downregulated DEmRNAs. Among the DEmiRNAs, bta-miR-12051, bta-miR-211, bta-miR-222, and bta-miR-11986c were significantly downregulated. Among the DElncRNAs, 107132706, 104972497, 790183, and 104973640 were significantly upregulated. It is worth mentioning that the difference in the expression levels of 104972497 and 790183 was extremely significant (Figure 8). The qRT-PCR results for these DEGs were consistent with the sequencing data.

FIGURE 8
www.frontiersin.org

Figure 8. The relative expression of DEGs. The relative expression of the DEGmRNAs (ApoA4, SCD, CYP7A1, CXCL3, LIPG, IGF1, ATP12A, CACNG4, BMP8A, and BMP8B), DEmiRNAs (bta-miR-12051, bta-miR-211, bta-miR-222, and bta-miR-11986c), and DElncRNAs (107132706, 104972497, 790183, and 104973640). Data (n = 3) are expressed as the mean ± SEM. For the cadmium-treated group compared with control group: *p < 0.05, **p < 0.01, ***p < 0.001.

Discussion

Hunan is rich in mineral resources, but the problem of Cd-contaminated soil is severe. Cd is a toxic heavy metal with high mobility in soil, and can greatly endanger the health of humans and animals when it enters their bodies through diet and water. The Xiangxi yellow heifer is a special breed in Hunan that feeds on local forage. Therefore, they are highly susceptible to Cd contamination. Chronic Cd exposure is mainly deposited in the liver and can induce liver inflammation (43), and even acute and chronic liver lesions (44). Moreover, redox imbalance, mitochondrial dysfunction, fatty acid oxidation inhibition, and apoptosis are considered the main adverse effects of Cd exposure (26, 45, 46). In this study, the whole transcriptome of liver tissues of Xiangxi yellow heifers fed high and low Cd concentrations was determined, and DEGs were analyzed to study the mechanisms of Cd toxicity in the heifer liver. The critical DEmRNAs involved in Cd detoxification between the control and cadmium-treated groups were analyzed using KEGG and GO analyses.

The predominant toxic effect of Cd exposure is ROS production, which promotes the production of mitochondrial phagosomes that engulf damaged mitochondria (47). Cd exposure has been reported to increase the expression of adenosine 5'-monophosphate (AMP)–activated protein kinase (AMPK) (48). AMPK is also involved in autophagy regulation (49). This is consistent with the KEGG enrichment results of the current study. The effect of Cd on mitochondrial autophagy in the liver of Xiangxi yellow heifers has been demonstrated. SCD is enriched in the AMPK signaling pathway and mediates fatty acid metabolism in the liver and adipogenesis (50). SCD-1, a common isoform of this enzyme (51), is a crucial gene in hepatic steatosis and plays a straightforward role in intracellular fat accumulation in NAFLD. In addition, excessive accumulation of liver fat leads to liver inflammation (52). Research has found that the mRNA levels of SCD-1 are also significantly higher in rats administered with CdCl2 in their drinking water (53). NAFLD is associated with elevated mRNA levels of SCD-1. High Cd exposure significantly upregulated SCD expression in our study. This demonstrates that the entry of Cd into the bodies of humans and animals can seriously disrupt the metabolic system and even cause metabolic diseases. In the current study, metabolic pathways were enriched with the most DEmRNAs. The endothelial lipase gene (LIPG) can hydrolyze phosphatidylcholine (PC) in HDL and release free fatty acids (FFAs) (54). In the present study, LIPG expression was significantly downregulated, suggesting that HDL hydrolysis was inhibited in the liver. Thus, cholesterol in the surrounding tissues could not be transported to bile acids in time, and fat accumulated in the liver. Cholesterol 7α-hydroxylase (CYP7A1) is in the metabolic pathway and belongs to the cytochrome P450 enzyme (CYP450) family. It is correlated with inflammatory cytokines and lysosomal function (55, 56). CYP7A1 overexpression strongly induces autophagy (57), and CYP7A1 is required for autophagy induction. In the present study, the expression levels of CYP7A1 mRNA were significantly altered, as shown by qRT-PCR. This indicates that Cd induces autophagy in the liver.

The PPAR signaling pathway was also significantly enriched in the liver. PPARs are a class of nuclear receptors that bind to fatty acids and their derivatives (58). PPAR-α has been reported to have transcriptional regulatory effects on genes involved in gluconeogenesis and lipid transport, and PPAR-γ controls glucose metabolism, adipogenesis, and adipocyte differentiation. Part of PPAR-γ inhibits inflammatory responses by competitively inhibiting inflammatory pathways and inflammatory mediator production (59). PPAR-β promotes hepatic fatty acid oxidation and limits inflammation (6062). Evidence suggests that PPAR-γ is closely related to the toxicity of heavy metals and that Cr inhibits fibroblast differentiation by reducing PPAR-γ (63). SCD and CYP7A1 genes enriched in the PPAR signaling pathway are closely related to cellular inflammation and adipogenesis in humans and animals. In this study, the mRNA expression levels of SCD and CYP7A1 were significantly upregulated in heifers under chronic Cd exposure. These findings show that this pathway is indeed involved in the toxic effects of heavy metals. In addition to the DEmRNAs in these three significantly different pathways, there were several DEGs related to liver physiological functions after Cd exposure (CACNG4, CXCL3, ATP12A, APOA4, BMP8A, BMP8B, and IGF-1). CACNG4 regulates voltage-gated calcium channels, and studies have shown that Cd exposure increases intracellular Ca2+ concentration (64). Elevated Ca2+ channel activity leads to increased CACNG4 expression (65, 66). The elevated levels of CACNG4 expression in this study suggest an imbalance in cellular Ca2+ concentrations in the liver owing to Cd exposure. CXCL3 is a pro-inflammatory cytokine that is induced alongside MMP12 in mice undergoing oxidative stress (67). CXCL3 exerts its anti-apoptotic effect mainly by counteracting the reversal of the cationic gradient during apoptosis (68). The significant upregulation of CXCL3 in this study demonstrates that Cd causes liver inflammation. ATP12A is a non-gastric, functionally active H+/K+-ATPase that is expressed mainly in the human and porcine respiratory tract, human prostate tissue, and bone marrow mono-nuclear cells. (6971). In the present study, ATP12A expression was measured in the livers of heifers. The marked downregulation of ATP12A expression levels suggests that Cd reversed the cationic gradient in the liver, leading to apoptosis. APOA4, a lipid-binding protein, is involved in the synthesis of cholesterol, steroids, and other lipids, and improves glucose homeostasis (72). The anti-inflammatory and antioxidant activities of APOA4 are crucial for maintaining hepatic homeostasis (73, 74). In our study, APOA4 was downregulated, indicating that Cd injury affected fat metabolism, inflammation, and oxidative stress in the liver. Furthermore, BMP8A regulates ROS homeostasis via Nrf2 (75), and BMP8B is involved in the regulation of growth traits in heifers (76). In addition, IGF-1 modulates cell proliferation and apoptosis through the PI3/AKT and MAPK pathways (77). Overall, these seven genes were not enriched in significantly different pathways in the transcriptome data. Nevertheless, they indicated that Cd exposure is involved in hepatotoxicity regarding several crucial metabolic pathways and life processes, including inflammation, apoptosis, oxidative stress, and lipid metabolism. These findings provide insight into the mechanisms underlying Cd-induced hepatotoxicity in mammals.

Studies have demonstrated that miRNAs and lncRNAs can regulate the cell cycle, apoptosis, DNA damage, and inflammation in the mechanisms underlying Cd toxicity (7881). Partial DElncRNAs and DEmiRNAs were selected for further analyses. They were verified by RT-PCR and matched with the data obtained by sequencing, which proved the accuracy and reliability of the sequencing data. The target mRNAs of the DElncRNAs and DEmiRNAs were enriched in the AMPK signaling pathway, glycerophospholipid metabolism, metabolic pathways, and peroxisomes. Therefore, the DElncRNAs, DEmiRNAs, and DEmRNAs identified in this study may be used as biomarkers of Cd toxicity in the future, which is of major research significance.

Conclusions

In the present study, we determined the gene expression patterns of Cd-exposed livers and DEGs in Xiangxi yellow heifers. The AMPK pathway involved in autophagy regulation, and the PPAR pathway involved in lipid metabolism, were significantly enriched in the livers of Xiangxi yellow cattle fed with high Cd diets. DEmRNAs included the anti-inflammatory and antioxidant APOA4, anti-apoptotic ATP12A, cholesterol metabolism–related LIPG, CXCL3 involved in fatty liver development, cholesterol-metabolizing enzyme CYP7A1, and lipid metabolism–related SCD. DEmiRNAs included bta-miR-12051, bta-miR-211, bta-miR-222, and bta-miR-11986c. DElncRNAs included 107132706, 104972497, 790183, and 104973640. In summary, these genes play a crucial role in Cd-induced liver damage and are critical for future studies of liver diseases associated with Cd exposure.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA801770; https://www.ncbi.nlm.nih.gov/, PRJNA801802.

Ethics Statement

The animal study was reviewed and approved by the Tab of Animal Experimental Ethical Inspection, JLU. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author Contributions

YW and KY performed the experiments and analyzed the data. All authors read and approved the final version of the manuscript.

Funding

This work was supported by the Key Research and Development Program of Jilin Province (20210202048NC and 20210202103NC), the Key Research and Development Project of Hunan Province, China (Nos. 2020NK2066 and 2022NK2025), the Construction of Modern Agricultural Industrial Technology System in Hunan Province (Hunan Financial Agriculture Guide 2019 [97]), and the Xiangxi Yellow Cattle Engineering Technology Center (2019TP2010).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.846662/full#supplementary-material

References

1. Chen H, Teng Y, Lu S, Wang Y, Wang J. Contamination features and health risk of soil heavy metals in China. Sci Total Environ. (2015) 512:143–53. doi: 10.1016/j.scitotenv.2015.01.025

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Yong hua L, Lin sheng Y, Yan fang J, Hong fei S. Lead uptake by plant in soil-plant system at lead-zinc deposit area of western Hunan province. Environ Sci Technol. (2008) 29:197–201. doi: 10.13227/j.hjkx.2008.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Jie S, Li M, Gan M, Zhu J, Yin H, Liu X. Microbial functional genes enriched in the Xiangjiang river sediments with heavy metal contamination. BMC Microbiol. (2016) 16:800. doi: 10.1186/s12866-016-0800-x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Liua H, Probst A, Liao B. Metal contamination of soils and crops affected by the chenzhou leadzinc mine spill (Hunan, China). Sci Total Environ. (2005) 339:153–66. doi: 10.1016/j.scitotenv.2004.07.030

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Paul N, Ming L, Yong-GZ. Occurrence and partitioning of cadmium, arsenic and lead in mine impacted paddy ricehunan, China. Environ Sci Technol. (2009) 43:637–42. doi: 10.1021/es802412r

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Hong Y, Liang Q. Study on the resource characteristics an breeding technology of xiangxi yellow cattle [in Chinese]. Mod Animal Husbandry Technol. (2021) 7:22–4. doi: 10.19369/j.cnki.2095-9737.2021.07.009

CrossRef Full Text | Google Scholar

7. Roggeman S, de Boeck G, De Cock H, Blust R, Bervoets L. Accumulation and detoxification of metals and arsenic in tissues of cattle (Bos Taurus), and the risks for human consumption. Sci Total Environ. (2014) 467:175–84. doi: 10.1016/j.scitotenv.2013.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Alonso ML Benedito JL Miranda M Castillo C Herncattle Herncattle cumulation and niil hebiomonitors of soil arsenic copper and zinc concentrations in galicia (Nw Spain). Arch Environ Contam Toxicol. (2002) 43:103–8. doi: 10.1007/s00244-002-1168-5

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Tianyue P, Jieling W. Study on the contamination status, risk and management strategies of heavy metal contaminated agricultural soil (in chinese). In: Proceedings of the 2021 Annual Science and Technology Conference of the Chinese Society of Environmental Science - Environmental Engineering Technology Innovation and Application Session. Tianjin: Peng Tianyue, Wang Jieling (2021).

Google Scholar

10. T021 and Application Sessionchnology Ince of the Cparameters as biomarkers of cadmium and lead exposure and effects in wild wood mice (Apodemus Sylvaticus) living along a pollution gradient. Chemosphere. (2015) 138:940–6. doi: 10.1016/j.chemosphere.2014.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ramírez-Elías MG, Kolosovas-Machuca ES, Kershenobich D, Guzmán C, Escobedo G, González FJ. Evaluation of liver fibrosis using Raman spectroscopy and infrared thermography: a pilot study. Photodiagnosis Photodyn Ther. (2017) 19: :278–83. doi: 10.1016/j.pdpdt.2017.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

12. de Lima EC, de Moura CFG, Silva MJD, Vilegas W, Santamarina AB, Pisani LP, et al. Therapeutical properties of mimosa caesalpiniifolia in rat liver intoxicated with cadmium. Environ Sci Poll Res. (2020) 27:10981–9. doi: 10.1007/s11356-019-07455-2

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Cupertino MC, Costa KLC, Santos DCM, Novaes RD, Condessa SS, Neves AC, et al. Long-lasting morphofunctional remodelling of liver parenchyma and stroma after a single exposure to low and moderate doses of cadmium in rats. Int J Exp Pathol. (2013) 94:343–51. doi: 10.1111/iep.12046

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Li Y, Shen R, Wu H, Yu L, Wang Z, Wang D. Liver changes induced by cadmium poisoning distinguished by confocal raman imaging. Spectrochimica Acta Part A: Mol Biomol Spectro. (2020) 225:117483. doi: 10.1016/j.saa.2019.117483

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Koz483 10. mica Acta Part A poisoning distinguished by confocal raman imaginrich extract from the berries of aronia Melanocarpa L. on collagen metabolism in the liver: a study in an in vivo model of human environmental exposure to cadmium. Nutrients. (2020) 12:2766. doi: 10.3390/nu12092766

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Elisabetta Casalino GC. Cesare sblano and clemente landriscina molecular inhibitory mechanisms of antioxidant enzymes in rat liver and kidney by cadmium. Toxicology. (2002) 179:37–50. doi: 10.1016/s0300-483x(02)00245-7

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Branca JJV, Fiorillo C, Carrino D, Paternostro F, Taddei N, Gulisano M, et al. Cadmium-induced oxidative stress: focus on the central nervous system. Antioxidants. (2020) 9:492. doi: 10.3390/antiox9060492

PubMed Abstract | CrossRef Full Text | Google Scholar

18. He X, Gao J, Hou H, Qi Z, Chen H, Zhang X-X. Inhibition of mitochondrial fatty acid oxidation contributes to development of non-alcoholic fatty liver disease induced by environmental cadmium exposure. Environ Sci Technol. (2019) 53:13992–4000. doi: 10.1021/acs.est.9b05131

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Belyaeva EA, Korotkov SM. Mechanism of primary Cd2+-induced rat liver mitochondria dysfunction: discrete modes of Cd2+ action on calcium and thiol-dependent domains. Toxicol Appl Pharmacol. (2003) 192:56–68. doi: 10.1016/s0041-008x(03)00255-2

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Oguz E. Transmission electron microscopy study of the effects of cadmium and copper on fetal rat liver tissue. Biol Trace Elem Res. (2007) 115:127–35. doi: 10.1007/BF02686024

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Eşrefoglu M, Gül M, Dogru MI, Dogru A, Yürekli M. Adrenomedullin fails to reduce cadmium-induced oxidative damage in rat liver. Exp Toxicol Pathol. (2007) 58:367–74. doi: 10.1016/j.etp.2006.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Samarghandian S, Azimi-Nezhad M, Shabestari MM, Azad FJ, Farkhondeh T, Bafandeh F. Effect of chronic exposure to cadmium on serum lipid, lipoprotein and oxidative stress indices in male rats. Interdiscip Toxicol. (2015) 8:151–4. doi: 10.1515/intox-2015-0023

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Nemmiche S. Oxidative signaling response to cadmium exposure. Toxicol Sci. (2016) 156:kfw222. doi: 10.1093/toxsci/kfw222

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Srinivasan P, Li YH, Hsu DZ, Su SB, Liu MY. Ostensibly ineffectual doses of cadmium and lipopolysaccharide causes liver damage in rats. Hum Exp Toxicol. (2010) 30:624–35. doi: 10.1177/0960327110376553

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hossein-Khannazer N, Azizi G, Eslami S, Alhassan Mohammed H, Fayyaz F, Hosseinzadeh R, et al. The effects of cadmium exposure in the induction of inflammation. Immunopharmacol Immunotoxicol. (2019) 42:1–8. doi: 10.1080/08923973.2019.1697284

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Ren L, Qi K, Zhang L, Bai Z, Ren C, Xu X, et al. Glutathione might attenuate cadmium-induced liver oxidative stress and hepatic stellate cell activation. Biol Trace Elem Res. (2019) 191:443–52. doi: 10.1007/s12011-019-1641-x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Th1641443e Elem Resmium-cellular signaling cascades: to be or not to be? Toxicol Appl Pharmacol. (2009) 238:221–39. doi: 10.1016/j.taap.2009.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Huo J, Dong A, Yan J, Dong A. Effects of cadmium on the activities of alt and ast as well as the content of tp in plasma of freshwater turtle mauremys reevesii. Environ Sci Pollution Res. (2020) 27:18025–8. doi: 10.1007/s11356-020-08338-7

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Huo J, Dong A, Yan J, Dong A. Effects of cadmium on the gene transcription of the liver in the freshwater turtle (Chinemys Reevesii). Environ Sci Pollution Res. (2020) 27:8431–8. doi: 10.1007/s11356-019-07432-9

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Zhu JY, Chan KM. Mechanism of cadmium-induced cytotoxicity on the zfl zebrafish liver cell line. Metallomics. (2012) 4:1064. doi: 10.1039/c2mt20134h

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Rosales-Cruz P, Domínguez-Pérez M, Reyes-Zárate E, Bello-Monroy O, Enríquez-Cortina C, Miranda-Labra R, et al. Cadmium exposure exacerbates hyperlipidemia in cholesterol-overloaded hepatocytes via autophagy dysregulation. Toxicology. (2018) 399:41–51. doi: 10.1016/j.tox.2018.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Satarug S. Long-term exposure to cadmium in food and cigarette smoke, liver effects and hepatocellular carcinoma. Curr Drug Metab. (2012) 13:257–71. doi: 10.2174/138920012799320446

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Julin B, Wolk A, Bergkvist L. Dietary cadmium exposure and risk of postmenopausal breast cancer a population-based prospective cohort study. Cancer Res. (2012) 72:1459–66. doi: 10.1158/0008-5472.CAN-11-0735

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Rapisarda V, Miozzi E, Loreto C, Matera S, Fenga C, Avola R, et al. Cadmium exposure and prostate cancer: insights, mechanisms and perspectives. Front Biosci. (2018) 23:1687–700. doi: 10.2741/4667

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Dijk EL, Jaszczyszyn Y, Thermes C. Library preparation methods for next-generation sequencing: tone down the bias. Exp Cell Res. (2014) 322:12–20. doi: 10.1016/j.yexcr.2014.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Li S, Tighe SW, Nicolet CM, Grove D, Levy S, Farmerie W, et al. Multi-platform assessment of transcriptome profiling using rna-seq in the abrf next-generation sequencing study. Nat Biotechnol. (2014) 32:915–25. doi: 10.1038/nbt.2972

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Matovi15echnolsment Matovi15echnolsment of transcriptome profiling using rna-seqoxidative stress induced by lead and/or cadmium in blood liver and kidneys. Food Chem Toxicol. (2015) 78:130–40. doi: 10.1016/j.fct.2015.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Soil Environmental Quality Risk Control Standard for Soil Contamination of Agricultural Land in China (2018). Available online at: http://www.mee.gov.cn/ywgz/fgbz/bz/bzwb/trhj/201807/t20180703_446029.shtml

Google Scholar

39. Li R, Li Y, Kristiansen K, Wang J. Soap: short oligonucleotide alignment program. Bioinformatics. (2008) 24:713–4. doi: 10.1093/bioinformatics/btn025

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. (2012) 9:357–9. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for rna-seq data with deseq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Shen S, Park JW, wLu Z-x, Lin L, Henry MD, Wu YN, et al. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-seq data. Proc Nat Acad Sci. (2014) 111:E55931E601. doi: 10.1073/pnas.1419161111

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Souza V, Escobar MdC, Gction of differential alternative splicing fro, et al. Acute cadmium exposure enhances Ap-1 DNA binding and induces cytokines expression and heat shock protein 70 in Hepg2 cells. Toxicology. (2004) 197:213–28. doi: 10.1016/j.tox.2004.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hyder O, Chung M, Cosgrove D, Herman JM Li Z, Firoozmand A, et al. Cadmium exposure and liver disease among us adults. J Gastrointestinal Surg. (2013) 17:1265–73. doi: 10.1007/s11605-013-2210-9

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Zou H, Wang T, Yuan J, Sun J, Yuan Y, Gu J, et al. Cadmium-induced cytotoxicity in mouse liver cells is associated with the disruption of autophagic flux via inhibiting the fusion of autophagosomes and lysosomes. Toxicol Lett. (2020) 321:32–43. doi: 10.1016/j.toxlet.2019.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Jin X, Jia T, Liu R, Xu S. The Antagonistic effect of selenium on cadmium-induced apoptosis via ppar-ct i3k/Akt pathway in chicken pancreas. J Hazard Mater. (2018) 357:355–62. doi: 10.1016/j.jhazmat.2018.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Zhang J. Teaching the basics of autophagy and mitophagy to redox biologists—mechanisms and experimental approaches. Redox Biol. (2015) 4:242–59. doi: 10.1016/j.redox.2015.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Chen M, Li X, Fan R, Yang J, Jin X, Hamid S, et al. Cadmium induces Bnip3-dependent autophagy in chicken spleen by modulating Mir-33-Ampk axis. Chemosphere. (2018) 194:396–402. doi: 10.1016/j.chemosphere.2017.12.026

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Villanueva-Paz M Cotn chicken ido-Maraver J Oropesa-mechanisms and experimental approaches. AMPK regulation of cell growth, apoptosis, autophagy, and bioenergetics. Experientia Supplementum. (2016) 107:45–71. doi: 10.1007/978-3-319-43589-3_3

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Zhu X, Bian H, Wang L, Sun X, Xu X, Yan H, et al. Berberine attenuates non-alcoholic hepatic steatosis through the Ampk-Srebp-1c-Scd1 pathway. Free Rad Biol Med. (2019) 141:192–204. doi: 10.1016/j.freeradbiomed.2019.06.019

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Powell DA. An overview of patented small molecule stearoyl coenzyme-a desaturase inhibitors (2009–2013). Expert Opin Ther Pat. (2014) 24:155–75. doi: 10.1517/13543776.2014.851669

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Liu X, Strable MS, Ntambi JM. Stearoyl coa desaturase 1: role in cellular inflammation and stress. Adv Nutr. (2011) 2:15–22. doi: 10.3945/an.110.000125

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Zhang J, Wang Y, Fu L, Feng YJ, Ji YL, Wang H, et al. Subchronic cadmium exposure upregulates the mrna level of genes associated to hepatic lipid metabolism in adult female Cd1 mice. J Applied Toxicol. (2018) 38:1026–35. doi: 10.1002/jat.3612

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Schaltenberg N, John C, Heine M, Haumann F, Rinninger F, Scheja L, et al. Endothelial lipase is involved in cold-induced high-density lipoprotein turnover and reverse cholesterol transport in mice. Front Cardiovasc Med. (2021) 8:8235. doi: 10.3389/fcvm.2021.628235

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Song K-H, Li T, Owsley E, Strom S, Chiang JYL. Bile acids activate fibroblast growth factor 19 signaling in human hepatocytes to inhibit cholesterol 7α-hydroxylase gene expression. Hepatology. (2009) 49:297–305. doi: 10.1002/hep.22627

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Song K-H, Chiang JYL. Glucagon and camp inhibit cholesterol 7a-hydroxylase (Cyp7a1) gene expression in human hepatocytes: discordant regulation of bile acid synthesis and gluconeogenesis. Hepatology. (2006) 43:117–25. doi: 10.1002/hep.20919

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Wang Y, Ding Y, Li J, Chavan H, Matye D, Ni H-M, et al. Targeting the enterohepatic bile acid signaling induces hepatic autophagy via a Cyp7a1–Akt–Mtor axis in mice. Cellular Mol Gastroenterol Hepatol. (2017) 3:245–60. doi: 10.1016/j.jcmgh.2016.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Caillaud M, Patel NH, White A, Wood M, Contreras KM, Toma W, et al. Targeting peroxisome proliferator-activated receptor-eroxisome prolireduce paclitaxel-induced peripheral neuropathy. Brain Behav Immun. (2021) 93:172–85. doi: 10.1016/j.bbi.2021.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Qi L, Chen Y, Shi K, Ma H, Wei S, Sha Z. Combining of transcriptomic and proteomic data to mine immune-related genes and proteins in the liver of cynoglossus semilaevis challenged with vibrio anguillarum. Comparative Biochem Physiol Part D: Genom Proteom. (2021) 39:100864. doi: 10.1016/j.cbd.2021.100864

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Fuchs C, Traussnigg S, Trauner M. Nuclear receptor modulation for the treatment of non-alcoholic fatty liver disease. Semin Liver Dis. (2016) 36:069–86. doi: 10.1055/s-0036-1571296

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Gross B, Pawlak M, Lefebvre P, Staels B. Ppars in Obesity-Induced T2dm, Dyslipidaemia and Nafld. Nature Rev Endocrinol. (2016) 13:36–49. doi: 10.1038/nrendo.2016.135

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Norikazu Maeda MT, Tohru F, Shinji K, Hitoshi N, Ken K, Hiroyuki N, et al. Pparγ ligands increase expression and plasma concentrations of adiponectin, an adipose-derived protein. Diabetes. (2001) 50:2094–9. doi: 10.2337/diabetes.50.9.2094

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Martini CN, Brandani JN, Gabrielli M, Vila MdC. Effect of hexavalent chromium on proliferation and differentiation to adipocytes of 3t3-l1 fibroblasts. Toxicol in vitro. (2014) 28:700–6. doi: 10.1016/j.tiv.2014.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Kania E Pajhromium Pajhromium on proliferation and homeostasis and ER stress in control of autophagy in cancer cells. Biomed Res Int. (2015) 2015:1–12. doi: 10.1155/2015/352794

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Kanwar N, Carmine-Simmen K, Nair R, Wang C, Moghadas-Jafari S, Blaser H, et al. Amplification of a calcium channel subunit cacng4 increases breast cancer metastasis. EBioMedicine. (2020) 52:102646. doi: 10.1016/j.ebiom.2020.102646

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Cao Y, Li R, Li Y, Zhang T, Wu N, Zhang J, et al. Identification of transcription factor-gene regulatory network in acute myocardial infarction. Heart, Lung and Circulation. (2017) 26:343–53. doi: 10.1016/j.hlc.2016.06.1209

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Yi-Jhu Lu, Yee-Jee Jan, Chih-Hui Chin. Expression of nik-related kinase in smooth muscle cells attenuates vascular inflammation and intimal hyperplasia. Aging. (2020) 12:7511–33. doi: 10.18632/aging.103104

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Jakab M, Hofer S, Ravasio A, Huber F, Schmidt S, Hitzl W, et al. The putative role of the non-gastric H+/K+-Atpase Atp12a (Atp1al1) as anti-apoptotic ion transporter: effect of the H+/K+Atpase inhibitor Sch28080 on butyrate-stimulated myelomonocytic Hl-60 cells. Cellular Physiol Biochem. (2014) 34:1507–26. doi: 10.1159/000366355

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Viral S. Shah, David K Meyerholz, Xiao Xiao. Tang airway acidification initiates host defense abnormalities in cystic fibrosis mice. Science. (2016) 351:503–7. doi: 10.1126/science.aad5589

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Gorrieri G, Scudieri P, Caci E, Schiavon M, Tomati V, Sirci F, et al. Goblet Cell hyperplasia requires high bicarbonate transport to support mucin release. Sci Reports. (2016) 6:36016. doi: 10.1038/srep36016

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Jakab M, Ketterl N, Fhigh bicarbonate transport to supSirci, et al. The H+/K+ Atpase inhibitor Sch-28080 inhibits insulin secretion and induces cell death in Ins-1e rat insulinoma cells. Cellular Physiol Biochem. (2017) 43:1037–51. doi: 10.1159/000481701

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Wang Z, Wang L, Zhang Z, Feng L, Song X, Wu J. Apolipoprotein a-Iv involves in glucose and lipid metabolism of rat. Nutr Metabol. (2019) 16:367. doi: 10.1186/s12986-019-0367-2

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Xu L, Shi L, Liu L, Liang R, Li Q, Li J, et al. Analysis of liver proteome and identification of critical proteins affecting milk fat, protein, and lactose metabolism in dariy cattle with Itraq. Proteomics. (2019) 19:1800387. doi: 10.1002/pmic.201800387

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Wang Y, Yang Z, Wei Y, Li X, Li S. Apolipoprotein A4 regulates the immune response in carbon tetrachloride-induced chronic liver injury in mice. Int Immunopharmacol. (2021) 90:107222. doi: 10.1016/j.intimp.2020.107222

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Yu Yp, Cai Lc, Wang Xy, Cheng Sy, Zhang Dm, Jian Wg, et al. Bmp8a promotes survival and drug resistance via Nrf2/Trim24 signaling pathway in clear cell renal cell carcinoma. Cancer Sci. (2020) 111:1555–66. doi: 10.1111/cas.14376

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Cao XK, Wang J, Lan XY, Lei CZ, Zhang CL, Qi XL, et al. Genetic variants in Bmp8b gene are associated with growth traits in chinese native cattle. Gene. (2013) 532:115–20. doi: 10.1016/j.gene.2013.09.059

PubMed Abstract | CrossRef Full Text | Google Scholar

77. J09. 05915associated with growth traits in chinese native cattlesignaling pathway in viral infections. Viruses. (2021) 13:1488. doi: 10.3390/v13081488

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Zhou Z, Huang Z, Chen B. Lncrnactrat00000446135 is a novel biomarker of cadmium toxicity in 16hbe cells, rats, and cd-exposed workers and regulates DNA damage and repair. Toxicol Res. (2020) 9:823–34. doi: 10.1093/toxres/tfaa088/6029993

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Huang Q, Lu Q, Chen B, Shen H, Liu Q, Zhou Z, et al. Lncrna-Malat1 as a novel biomarker of cadmium toxicity regulates cell proliferation and apoptosis. Toxicol Res. (2017) 6:361–71. doi: 10.1039/c6tx00433d

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Yuan W, Liu L, Liang L, Huang K, Deng Y, Dong M, et al. Mir-122-5p and Mir-326-3p: potential novel biomarkers for early detection of cadmium exposure. Gene. (2020) 724:144156. doi: 10.1016/j.gene.2019.144156

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Chen J, Zhang S, Tong J, Teng X, Zhang Z, Li S, et al. Whole transcriptome-based mirna-mrna network analysis revealed the mechanism of inflammation-immunosuppressive damage caused by cadmium in common carp spleens. Sci Total Environ. (2020) 717:137081. doi: 10.1016/j.scitotenv.2020.137081

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: RNA-seq, cadmium, liver, gene ontology, Kyoto Encyclopedia of Genes

Citation: Wei Y, Yi K, Shen C, Chen X, Iqbal T, Cao M, Chen T, Luo Y, Li J, Zhou X, Li C and Chen L (2022) Whole Transcriptome Profiling of the Effects of Cadmium on the Liver of the Xiangxi Yellow Heifer. Front. Vet. Sci. 9:846662. doi: 10.3389/fvets.2022.846662

Received: 31 December 2021; Accepted: 07 March 2022;
Published: 14 April 2022.

Edited by:

Ahmed Abdeen, Benha University, Egypt

Reviewed by:

Ameer Megahed, University of Illinois at Urbana–Champaign, United States
Lin Wang, Shandong Agricultural University, China

Copyright © 2022 Wei, Yi, Shen, Chen, Iqbal, Cao, Chen, Luo, Li, Zhou, Li and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lu Chen, luchen@jlu.edu.cn

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.