Transcriptome analysis of long non-coding RNA and mRNA Profiles in VSV-infected BHK-21 Cells.

BACKGROUND: Vesicular stomatitis virus (VSV) is a typical non-segmented negative-sense RNA virus of the genus Vesiculovirus in the family Rhabdoviridae. VSV can infect a wide range of animals, including humans, with oral blister epithelial lesions. VSV is an excellent model virus with a wide range of applications as a molecular tool, a vaccine vector, and an oncolytic vector. To further understand the interaction between VSV and host cells and to provide a theoretical basis for the application prospects of VSV, we analyzed the expression of host differentially expressed genes (DEGs) during VSV infection using RNA-Seq. RESULTS: Our analyses found a total of 1015 differentially expressed mRNAs and 161 differentially expressed LncRNAs in BHK-21 cells infected with VSV for 24 h compared with controls. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment showed that the differentially expressed lncRNAs and their target genes were mainly concentrated in pathways related to apoptosis, cancer, disease, and immune system activation, including the TNF, P53, MAPK, and NF-kappaB signaling pathways. The differentially expressed lncRNA can modulate immune processes by regulating genes involved in these signaling transmissions. Ten randomly selected DEGs, namely, Il12rb2, F2, Masp2, Mcl1, FGF18, Ripk1, Fas, BMF, POLK, and JAG1, were validated using RT-qPCR. As predicted through RNA-Seq analysis, these DEGs underwent either up- or downregulation, suggesting that they may play key regulatory roles in the pathways mentioned previously. CONCLUSIONS: Our study showed that VSV infection alters the host metabolic network and activates immune-related pathways, such as MAPK and TNF. The above findings provide unique insights for further study of the mechanism of VSV-host interactions and, more importantly, provide a theoretical basis for VSV as an excellent vaccine carrier.

MiRNA-Seq reveals key MicroRNAs involved in fat metabolism of sheep liver.

There is a genetic difference between Hu sheep (short/fat-tailed sheep) and Tibetan sheep (short/thin-tailed sheep) in tail type, because of fat metabolism. Previous studies have mainly focused directly on sheep tail fat, which is not the main organ of fat metabolism. The function of miRNAs in sheep liver fat metabolism has not been thoroughly elucidated. In this study, miRNA-Seq was used to identify miRNAs in the liver tissue of three Hu sheep (short/fat-tailed sheep) and three Tibetan sheep (short/thin-tailed sheep) to characterize the differences in fat metabolism of sheep. In our study, Hu sheep was in a control group, we identified 11 differentially expressed miRNAs (DE miRNAs), including six up-regulated miRNAs and five down-regulated miRNAs. Miranda and RNAhybrid were used to predict the target genes of DE miRNAs, obtaining 3,404 target genes. A total of 115 and 67 GO terms as well as 54 and 5 KEGG pathways were significantly (padj < 0.05) enriched for predicted 3,109 target genes of up-regulated and 295 target genes of down-regulated miRNAs, respectively. oar-miR-432 was one of the most up-regulated miRNAs between Hu sheep and Tibetan sheep. And SIRT1 is one of the potential target genes of oar-miR-432. Furthermore, functional validation using the dual-luciferase reporter assay indicated that the up-regulated miRNA; oar-miR-432 potentially targeted sirtuin 1 (SIRT1) expression. Then, the oar-miR-432 mimic transfected into preadipocytes resulted in inhibited expression of SIRT1. This is the first time reported that the expression of SIRT1 gene was regulated by oar-miR-432 in fat metabolism of sheep liver. These results could provide a meaningful theoretical basis for studying the fat metabolism of sheep.

Transcriptome Comparison Reveals the Difference in Liver Fat Metabolism between Different Sheep Breeds.

Hu sheep and Tibetan sheep are two commonly raised local sheep breeds in China, and they have different morphological characteristics, such as tail type and adaptability to extreme environments. A fat tail in sheep is the main adipose depot in sheep, whereas the liver is an important organ for fat metabolism, with the uptake, esterification, oxidation, and secretion of fatty acids (FAs). Meanwhile, adaptations to high-altitude and arid environments also affect liver metabolism. Therefore, in this study, RNA-sequencing (RNA-seq) technology was used to characterize the difference in liver fat metabolism between Hu sheep and Tibetan sheep. We identified 1179 differentially expressed genes (DEGs) (Q-value < 0.05) between the two sheep breeds, including 25 fat-metabolism-related genes. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, 16 pathways were significantly enriched (Q-value < 0.05), such as the proteasome, glutamatergic synapse, and oxidative phosphorylation pathways. In particular, one of these pathways was enriched to be associated with fat metabolism, namely the thermogenesis pathway, to which fat-metabolism-related genes such as ACSL1, ACSL4, ACSL5, CPT1A, CPT1C, SLC25A20, and FGF21 were enriched. Then, the expression levels of ACSL1, CPT1A, and FGF21 were verified in mRNA and protein levels via qRT-PCR and Western blot analysis between the two sheep breeds. The results showed that the mRNA and protein expression levels of these three genes were higher in the livers of Tibetan sheep than those of Hu sheep. The above genes are mainly related to FAs oxidation, involved in regulating the oxidation of liver FAs. So, this study suggested that Tibetan sheep liver has a greater FAs oxidation level than Hu sheep liver. In addition, the significant enrichment of fat-metabolism-related genes in the thermogenesis pathway appears to be related to plateau-adaptive thermogenesis in Tibetan sheep, which may indicate that liver- and fat-metabolism-related genes have an impact on adaptive thermogenesis.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails.

BACKGROUND: Hu sheep and Tibetan sheep in China are characterized by fat tails and thin tails, respectively. Several transcriptomes have been conducted in different sheep breeds to identify the differentially expressed genes (DEGs) underlying this trait. However, these studies identified different DEGs in different sheep breeds. RESULTS: Hence, RNA sequencing was performed on Hu sheep and Tibetan sheep. We obtained a total of 45.57 and 43.82 million sequencing reads, respectively. Two libraries mapped reads from 36.93 and 38.55 million reads after alignment to the reference sequences. 2108 DEGs were identified, including 1247 downregulated and 861 upregulated DEGs. GO and KEGG analyses of all DEGs demonstrated that pathways were enriched in the regulation of lipolysis in adipocytes and terms related to the chemokine signalling pathway, lysosomes, and glycosaminoglycan degradation. Eight genes were selected for validation by RT-qPCR. In addition, the transfection of BMP2 overexpression into preadipocytes resulted in increased PPAR-γ expression and expression. BMP2 potentially induces adipogenesis through LOX in preadipocytes. The number of lipid drops in BMP2 overexpression detected by oil red O staining was also greater than that in the negative control. CONCLUSION: In summary, these results showed that significant genes (BMP2, HOXA11, PPP1CC and LPIN1) are involved in the regulation of adipogenesis metabolism and suggested novel insights into metabolic molecules in sheep fat tails.

Oar-miR-432 Regulates Fat Differentiation and Promotes the Expression of BMP2 in Ovine Preadipocytes.

The fat tail is a unique characteristic of sheep that represents energy reserves and is a complex adaptative mechanism of fat-tailed sheep to environmental stress. MicroRNA plays a significant role as regulators at the posttranscriptional level, but no studies have explained the molecular mechanisms of miRNA which regulate fat deposition in sheep tails. In this study, mRNA and miRNA analysis examined tail fat tissue from three Hu fat-tailed and three Tibetan thin-tailed sheep. After aligning to the reference sequences, 2,108 differentially expressed genes and 105 differential expression miRNAs were identified, including 1,247 up- and 861 downregulated genes and 43 up- and 62 downregulated miRNAs. Among these differentially expressed miRNAs, oar-miR-432 was one of the most downregulated miRNAs between Hu sheep and Tibetan sheep, and 712 genes were predicted to be targeted by oar-miR-432, 80 of which overlapped with DEGs. The Gene Ontology analysis on these genes showed that BMP2, LEP, GRK5, BMP7, and RORC were enriched in fat cell differentiation terms. The genes for BMP2 targeted by oar-miR-432 were examined using dual-luciferase assay. The oar-miR-432 mimic transfected into preadipocytes resulted in increased expression of BMP2. The marker gene PPAR-γ of fat differentiation had a lower expression than the negative control on days 0, 2, and 4 after induced differentiation. The decrease in the number of lipids in the oar-miR-432 mimic group detected by oil red O stain was also less than that in the negative control. This is the first study to reveal the fat mechanisms by which oar-miR-432 inhibits fat differentiation and promotes the expression of BMP2 in sheep tails.

Transcriptome reveals key microRNAs involved in fat deposition between different tail sheep breeds.

MicroRNA (miRNA) is a kind of noncoding RNA whose function involved in various biological processes in neuronal maturation and adipocyte cells, such as differentiation, proliferation, development, apoptosis, and metabolism. Herein, miRNA-Seq was used to identify miRNAs in the tail fat tissue of Hu sheep (short-fat-tailed) and Tibetan sheep (short-thin-tailed). In this study, 155 differentially expression miRNAs (DE miRNAs) were identified, including 78 up-regulated and 77 down-regulated. Among these DE miRNAs, 17 miRNAs were reported and related with lipid metabolism. MiRanda and RNAhybrid software were used to predict the target genes of DE miRNAs, obtaining the number of targeting relationships is 38553. Target genes were enriched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). 742 terms and 302 single pathways are enriched, including lipid metabolic process, response to lipid, cellular lipid catabolic process, lipid catabolic process, cellular lipid metabolic process, inositol lipid-mediated signaling, calcium channel activity, PI3K-Akt signaling pathway, MAPK signaling pathway, ECM-receptor interaction, AMPK signaling pathway, Wnt signaling pathway and TGF-beta signaling pathway. Notably, miR-379-5p was associated with tail fat deposition of sheep. Dual-Luciferase reporter assays showed miR-379-5p and HOXC9 had targeted relationship. The result of RT-qPCR showed that the expression trend of miR-379-5p and HOXC9 was opposite. miR-379-5p was down-regulated and highly expressed in tail adipose tissue of Tibetan sheep. HOXC9 was highly expressed in adipose tissue of Hu sheep. These results could provide a meaningful theoretical basis for studying the molecular mechanisms of sheep tail adipogenesis.

Genetic Signatures of Selection for Cashmere Traits in Chinese Goats.

Inner Mongolia and Liaoning cashmere goats in China are well-known for their cashmere quality and yield. Thus, they are great models for identifying genomic regions associated with cashmere traits. Herein, 53 Inner Mongolia cashmere goats, Liaoning cashmere goats and Huanghuai goats were genotyped, and 53,347 single-nucleotide polymorphisms (SNPs) were produced using the Illumina Caprine 50K SNP chip. Additionally, we identified some positively selected SNPs by analyzing Fst and XP-EHH. The top 5% of SNPs had selection signatures. After gene annotation, 222 and 173 candidate genes were identified in Inner Mongolia and Liaoning cashmere goats, respectively. Several genes were related to hair follicle development, such as TRPS1, WDR74, LRRC14, SPTLC3, IGF1R, PADI2, FOXP1, WNT10A and CSN3. Gene enrichment analysis of these cashmere trait-associated genes related 67 enriched signaling pathways that mainly participate in hair follicle development and stem cell pluripotency regulation. Furthermore, we identified 20 overlapping genes that were selected in both cashmere goat breeds. Among these overlapping genes, WNT10A and CSN3, which are associated with hair follicle development, are potentially involved in cashmere production. These findings may improve molecular breeding of cashmere goats in the future.

Selection Signatures Analysis Reveals Genes Associated with High-Altitude Adaptation in Tibetan Goats from Nagqu, Tibet.

Tibetan goat is an ancient breed, which inhabits the adverse conditions of the plateaus in China. To investigate the role of selection in shaping its genomes, we genotyped Tibetan goats (Nagqu Prefecture, above 4500 m) and three lowland populations (Xinjiang goats, Taihang goats and Huanghuai goats). The result of PCA, neighbor-joining (N-J) tree and model-based clustering showed that the genetic structure between the Tibetan goat and the three lowland populations has significant difference. As demonstrated by the di statistic, we found that some genes were related to the high-altitude adaptation of Tibetan goats. Functional analysis revealed that these genes were enriched in the VEGF (vascular endothelial growth factor) signaling pathway and melanoma, suggesting that nine genes (FGF2, EGFR, AKT1, PTEN, MITF, ENPEP, SIRT6, KDR, and CDC42) might have important roles in the high-altitude adaptation of Nagqu Tibetan goats. We also found that the LEPR gene was under the strongest selection (di value = 16.70), and it could induce upregulation of the hypoxic ventilatory response. In addition, five genes (LEPR, LDB1, EGFR, NOX4 and FGF2) with high di values were analyzed using q-PCR. Among them, we found that LEPR, LDB1 and FGF2 exhibited higher expression in the lungs of the Tibetan goats; LEPR, EGFR and LDB1 exhibited higher expression in the hearts of the Huanghuai goat. Our results suggest that LEPR, LDB1, EGFR and FGF2 genes may be related to the high-altitude adaptation of the goats. These findings improve our understanding of the selection of the high-altitude adaptability of the Nagqu Tibetan goats and provide new theoretical knowledge for the conservation and utilization of germplasm resources.

Verification and Analysis of Sheep Tail Type-Associated PDGF-D Gene Polymorphisms.

The aim of this study was to examine the correlation between the platelet-derived growth factor-D (PDGF-D) gene and sheep tail type character and explore the potential underlying mechanism. A total of 533 sheep were included in this study. Polymorphic sites were examined by Pool-seq, and individual genotype identification and correlation analysis between tail type data were conducted using the matrix-assisted laser desorption/ionization time-of-flight mass spectrometer (MALDI-TOF-MS) method. JASPART website was used to predict transcription factor binding sites in the promoter region with and without PDGF-D gene mutation. The effect of PDGF-D on adipogenic differentiation of sheep preadipocytes was investigated. Two single nucleotide polymorphism sites were identified: g.4122606 C > G site was significantly correlated with tail length, and g.3852134 C > T site was significantly correlated with tail width. g.3852134 C > T was located in the promoter region. Six transcription factor binding sites were eliminated after promoter mutation, and three new transcription factor binding sites appeared. Expression levels of peroxisome proliferator-activated receptor gamma (PPARγ) and lipoproteinlipase (LPL) were significantly up-regulated upon PDGF-D overexpression. Oil red O staining showed increased small and large oil drops in the PDGF-D overexpression group. Together these results indicate the PDGF-D gene is an important gene controlling sheep tail shape and regulating sheep tail fat deposition to a certain degree.

Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep.

With the intensified and large-scale development of sheep husbandry and global warming, sheep heat stress has become an increasingly important issue. However, little is known about the molecular mechanisms related to sheep responses to heat stress. In this study, transcriptomic analysis of liver tissues of sheep in the presence and absence of heat stress was conducted, with the goal of identifying genes and pathways related to regulation when under such stress. After a comparison with the sheep reference genome, 440,226,436 clean reads were obtained from eight libraries. A p-value ≤ 0.05 and fold change ≥ 2 were taken as thresholds for categorizing differentially expressed genes, of which 1137 were identified. The accuracy and reliability of the RNA-Seq results were confirmed by qRT-PCR. The identified differentially expressed genes were significantly associated with 419 GO terms and 51 KEGG pathways, which suggested their participation in biological processes such as response to stress, immunoreaction, and fat metabolism. This study's results provide a comprehensive overview of sheep heat stress-induced transcriptional expression patterns, laying a foundation for further analysis of the molecular mechanisms of sheep heat stress.