• Asian-Australasian Journal of Animal Sciences
• /
• v.29 no.3
• /
• pp.343-351
• /
• 2016

Although the chemical, physical, and nutritional properties of bovine milk have been extensively studied, only a few studies have attempted to characterize milk-synthesizing genes using RNA-seq data. RNA-seq data was collected from 21 Holstein samples, along with group information about milk production ability; milk yield; and protein, fat, and solid contents. Meta-analysis was employed in order to generally characterize genes related to milk production. In addition, we attempted to investigate the relationship between milk related traits, parity, and lactation period. We observed that milk fat is highly correlated with lactation period; this result indicates that this effect should be considered in the model in order to accurately detect milk production related genes. By employing our developed model, 271 genes were significantly (false discovery rate [FDR] adjusted p-value<0.1) detected as milk production related differentially expressed genes. Of these genes, five (albumin, nitric oxide synthase 3, RNA-binding region (RNP1, RRM) containing 3, secreted and transmembrane 1, and serine palmitoyltransferase, small subunit B) were technically validated using quantitative real-time polymerase chain reaction (qRT-PCR) in order to check the accuracy of RNA-seq analysis. Finally, 83 gene ontology biological processes including several blood vessel and mammary gland development related terms, were significantly detected using DAVID gene-set enrichment analysis. From these results, we observed that detected milk production related genes are highly enriched in the circulation system process and mammary gland related biological functions. In addition, we observed that detected genes including caveolin 1, mammary serum amyloid A3.2, lingual antimicrobial peptide, cathelicidin 4 (CATHL4), cathelicidin 6 (CATHL6) have been reported in other species as milk production related gene. For this reason, we concluded that our detected 271 genes would be strong candidates for determining milk production.

• Journal of Acupuncture Research
• /
• v.37 no.2
• /
• pp.128-135
• /
• 2020

Background: Classical acupuncture is being used in the treatment of rheumatoid arthritis (RA). To explore the biological response to acupuncture, a network-based analysis was performed on gene expression data collected from an animal model of RA treated with acupuncture. Methods: Gene expression data were obtained from published microarray studies on blood samples from rats with collagen induced arthritis (CIA) and non-CIA rats, both treated with manual acupuncture. The weighted gene co-expression network analysis was performed to identify gene clusters expressed in association with acupuncture treatment time and RA status. Gene ontology and pathway analyses were applied for functional annotation and network visualization. Results: A cluster of 347 genes were identified that differentially downregulated expression in association with acupuncture treatment over time; specifically in rats with CIA with module-RA correlation at 1 hour after acupuncture (-0.27; p < 0.001) and at 34 days after acupuncture (-0.33; p < 0.001). Functional annotation showed highly significant enrichment of porphyrin-containing compound biosynthetic processes (p < 0.001). The network-based analysis also identified a module of 140 genes differentially expressed between CIA and non-CIA in rats (p < 0.001). This cluster of genes was enriched for antigen processing and presentation of exogenous peptide antigen (p < 0.001). Other functional gene clusters previously reported in earlier studies were also observed. Conclusion: The identified gene expression networks and their hub-genes could help with the understanding of mechanisms involved in the pathogenesis of RA, as well understanding the effects of acupuncture treatment of RA.

• Asian-Australasian Journal of Animal Sciences
• /
• v.31 no.4
• /
• pp.595-606
• /
• 2018

Objective: The Fusarium mycotoxins of deoxynivalenol (DON) and zerolenone (ZEN) cause health hazards for both humans and farm animals. Therefore, the main intention of this study was to reveal DON and ZEN effects on the mRNA expression of pro-inflammatory cytokines and other immune related genes in the liver of piglets. Methods: In the present study, 15 six-week-old piglets were randomly assigned to the following three different dietary treatments for 4 weeks: control diet, diet containing 8 mg DON/kg feed, and diet containing 0.8 mg ZEN/kg feed. After 4 weeks, liver samples were collected and sequenced using RNA-Seq to investigate the effects of the mycotoxins on genes and gene networks associated with the immune systems of the piglets. Results: Our analysis identified a total of 249 differentially expressed genes (DEGs), which included 99 upregulated and 150 downregulated genes in both the DON and ZEN dietary treatment groups. After biological pathway analysis, the DEGs were determined to be significantly enriched in gene ontology terms associated with many biological pathways, including immune response and cellular and metabolic processes. Consistent with inflammatory stimulation due to the mycotoxin-contaminated diet, the following Kyoto encyclopedia of genes and genomes pathways, which were related to disease and immune responses, were found to be enriched in the DEGs: allograft rejection pathway, cell adhesion molecules, graft-versus-host disease, autoimmune thyroid disease (AITD), type I diabetes mellitus, human T-cell leukemia lymphoma virus infection, and viral carcinogenesis. Genome-wide expression analysis revealed that DON and ZEN treatments downregulated the expression of the majority of the DEGs that were associated with inflammatory cytokines (interleukin 10 receptor, beta, chemokine [C-X-C motif] ligand 9), proliferation (insulin-like growth factor 1, major facilitator superfamily domain containing 2A, insulin-like growth factor binding protein 2, lipase G, and salt inducible kinase 1), and other immune response networks (paired immunoglobulin-like type 2 receptor beta, Src-like-adaptor-1 [SLA1], SLA3, SLA5, SLA7, claudin 4, nicotinamide N-methyltransferase, thyrotropin-releasing hormone degrading enzyme, ubiquitin D, histone $H_2B$ type 1, and serum amyloid A). Conclusion: In summary, our results demonstrated that high concentrations DON and ZEN disrupt immune-related processes in the liver.

• Genes and Genomics
• /
• v.40 no.12
• /
• pp.1351-1361
• /
• 2018

A new cold tolerant germplasm resource named glutinous rice 89-1 (Gr89-1, Oryza sativa L.) can overwinter using axillary buds, with these buds being ratooned the following year. The overwintering seedling rate (OSR) is an important factor for evaluating cold tolerance. Many quantitative trait loci (QTLs) controlling cold tolerance at different growth stages in rice have been identified, with some of these QTLs being successfully cloned. However, no QTLs conferring to the OSR trait have been located in the perennial O. sativa L. To identify QTLs associated with OSR and to evaluate cold tolerance. 286 $F_{12}$ recombinant inbred lines (RILs) derived from a cross between the cold tolerant variety Gr89-1 and cold sensitive variety Shuhui527 (SH527) were used. A total of 198 polymorphic simple sequence repeat (SSR) markers that were distributed uniformly on 12 chromosomes were used to construct the linkage map. The gene ontology (GO) annotation of the major QTL was performed through the rice genome annotation project system. Three main-effect QTLs (qOSR2, qOSR3, and qOSR8) were detected and mapped on chromosomes 2, 3, and 8, respectively. These QTLs were located in the interval of RM14208 (35,160,202 base pairs (bp))-RM208 (35,520,147 bp), RM218 (8,375,236 bp)-RM232 (9,755,778 bp), and RM5891 (24,626,930 bp)-RM23608 (25,355,519 bp), and explained 19.6%, 9.3%, and 11.8% of the phenotypic variations, respectively. The qOSR2 QTL displayed the largest effect, with a logarithm of odds score (LOD) of 5.5. A total of 47 candidate genes on the qOSR2 locus were associated with 219 GO terms. Among these candidate genes, 11 were related to cell membrane, 7 were associated with cold stress, and 3 were involved in response to stress and biotic stimulus. OsPIP1;3 was the only one candidate gene related to stress, biotic stimulus, cold stress, and encoding a cell membrane protein. After QTL mapping, a total of three main-effect QTLs-qOSR2, qOSR3, and qOSR8-were detected on chromosomes 2, 3, and 8, respectively. Among these, qOSR2 explained the highest phenotypic variance. All the QTLs elite traits come from the cold resistance parent Gr89-1. OsPIP1;3 might be a candidate gene of qOSR2.

• Animal Bioscience
• /
• v.35 no.1
• /
• pp.115-125
• /
• 2022

Objective: Intramuscular fat (IMF) is a critical economic indicator of pork quality. Studies on IMF among different pig breeds have been performed via high-throughput sequencing, but comparisons within the same pig breed remain unreported. Methods: This study was performed to explore the gene profile and identify candidate long noncoding RNA (lncRNAs) and mRNAs associated with IMF deposition among Laiwu pigs with different IMF contents. Based on the longissimus dorsi muscle IMF content, eight pigs from the same breed and management were selected and divided into two groups: a high IMF (>12%, H) and low IMF group (<5%, L). Whole-transcriptome sequencing was performed to explore the differentially expressed (DE) genes between these two groups. Results: The IMF content varied greatly among Laiwu pig individuals (2.17% to 13.93%). Seventeen DE lncRNAs (11 upregulated and 6 downregulated) and 180 mRNAs (112 upregulated and 68 downregulated) were found. Gene Ontology analysis indicated that the following biological processes played an important role in IMF deposition: fatty acid and lipid biosynthetic processes; the extracellular signal-regulated kinase cascade; and white fat cell differentiation. In addition, the peroxisome proliferator-activated receptor, phosphatidylinositol-3-kinase-protein kinase B, and mammalian target of rapamycin pathways were enriched in the pathway analysis. Intersection analysis of the target genes of DE lncRNAs and mRNAs revealed seven candidate genes associated with IMF accumulation. Five DE lncRNAs and 20 DE mRNAs based on the pig quantitative trait locus database were identified and shown to be related to fat deposition. The expression of five DE lncRNAs and mRNAs was verified by quantitative real time polymerase chain reaction (qRT-PCR). The results of qRT-PCR and RNA-sequencing were consistent. Conclusion: These results demonstrated that the different IMF contents among pig individuals may be due to the DE lncRNAs and mRNAs associated with lipid droplets and fat deposition.

• Mycobiology
• /
• v.50 no.1
• /
• pp.66-78
• /
• 2022

The identification of oleaginous yeast species capable of simultaneously utilizing xylose and glucose as substrates to generate value-added biological products is an area of key economic interest. We have previously demonstrated that the Cutaneotrichosporon dermatis NICC30027 yeast strain is capable of simultaneously assimilating both xylose and glucose, resulting in considerable lipid accumulation. However, as no high-quality genome sequencing data or associated annotations for this strain are available at present, it remains challenging to study the metabolic mechanisms underlying this phenotype. Herein, we report a 39,305,439 bp draft genome assembly for C. dermatis NICC30027 comprised of 37 scaffolds, with 60.15% GC content. Within this genome, we identified 524 tRNAs, 142 sRNAs, 53 miRNAs, 28 snRNAs, and eight rRNA clusters. Moreover, repeat sequences totaling 1,032,129 bp in length were identified (2.63% of the genome), as were 14,238 unigenes that were 1,789.35 bp in length on average (64.82% of the genome). The NCBI non-redundant protein sequences (NR) database was employed to successfully annotate 11,795 of these unigenes, while 3,621 and 11,902 were annotated with the Swiss-Prot and TrEMBL databases, respectively. Unigenes were additionally subjected to pathway enrichment analyses using the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster of Orthologous Groups of proteins (COG), Clusters of orthologous groups for eukaryotic complete genomes (KOG), and Non-supervised Orthologous Groups (eggNOG) databases. Together, these results provide a foundation for future studies aimed at clarifying the mechanistic basis for the ability of C. dermatis NICC30027 to simultaneously utilize glucose and xylose to synthesize lipids.

• Journal of Veterinary Science
• /
• v.23 no.6
• /
• pp.90.01-90.13
• /
• 2022

Background: Insulin regulates glucose homeostasis and has important effects on metabolism, cell growth, and differentiation. Depending on the cell type and physiological context, insulin signal has specific pathways and biological outcomes in different tissues and cells. For studying the signal pathway of insulin on glycolipid metabolism in porcine embryonic fibroblast (PEF), we used high-throughput sequencing to monitor gene expression patterns regulated by insulin. Objectives: The goal of our research was to see how insulin affected glucose and lipid metabolism in PEFs. Methods: We cultured the PEFs with the addition of insulin and sampled them at 0, 48, and 72 h for RNA-Seq analysis in triplicate for each time point. Results: At 48 and 72 h, 801 and 1,176 genes were differentially expressed, respectively. Of these, 272 up-regulated genes and 264 down-regulated genes were common to both time points. Gene Ontology analysis was used to annotate the functions of the differentially expressed genes (DEGs), the biological processes related to lipid metabolism and cell cycle were dominant. And the DEGs were significantly enriched in interleukin-17 signaling pathway, phosphatidylinositol-3-kinase-protein kinase B signaling pathway, pyruvate metabolism, and others pathways related to lipid metabolism by Kyoto Encyclopedia of Genes and Genomes enrichment analysis. Conclusions: These results elucidate the transcriptomic response to insulin in PEF. The genes and pathways involved in the transcriptome mechanisms provide useful information for further research into the complicated molecular processes of insulin in PEF.

• Animal Bioscience
• /
• v.36 no.1
• /
• pp.144-155
• /
• 2023

Objective: Adipocyte differentiation is regulated by a variety of functional genes and noncoding RNAs. However, the role of miRNAs in lipid deposition of goat white adipose tissue is still unclear. Therefore, this study revealed the miRNA expression profile in goat subcutaneous adipocytes by sRNA-seq. Methods: The miRNA expressed in goat subcutaneous preadipocytes and the mature adipocytes were sequenced by sRNA-seq. The differentially expressed miRNAs (DEm) were screened and gene ontology (GO) and Kyoto encyclopedia for genes and genomes (KEGG) analyses were performed. Gain-of-function and loss-of-function combined with oil red O staining, Bodipy staining, and quantitative reverse-transcription polymerase chain reaction (qPCR) were utilized to determine the effect of miR-133a-3p on adipocyte differentiation. Results: A total of 218 DEm were screened out. The target genes of these DEm were significantly enriched in GO items such as biological regulation and in KEGG terms such as FAK signaling pathway and MAPK signaling pathway. qPCR verified that the expression trend of miRNA was consistent with miRNA-seq. The gain-of-function or loss-of-function of miR-133a-3p showed that it promoted or inhibited the accumulation of lipid droplets, and CCAAT enhancer binding protein α (C/EBPα) and C/EBPβ were extremely significantly up-regulated or down-regulated respectively (p<0.01), the loss-of-function also led to a significant down-regulation of peroxisome proliferator activated receptor gamma (PPARγ) (p<0.01). Conclusion: This study successfully identified miRNAs expression patterns in goat subcutaneous adipocytes, and functional identification indicates that miR-133a-3p is a positive regulator of the differentiation process of goat subcutaneous adipocytes. Our results lay the foundation for the molecular mechanism of lipid deposition in meat-source goats from the perspective of miRNA.

• Animal Bioscience
• /
• v.35 no.3
• /
• pp.367-376
• /
• 2022

Objective: The highly pathogenic avian influenza virus (HPAIV) is a threat to the poultry industry as well as the economy and remains a potential source of pandemic infection in humans. Antiviral genes are considered a potential factor for HPAIV resistance. Therefore, in this study, we investigated gene expression related to cytokine-cytokine receptor interactions by comparing resistant and susceptible Ri chicken lines for avian influenza virus infection. Methods: Ri chickens of resistant (Mx/A; BF2/B21) and susceptible (Mx/G; BF2/B13) lines were selected by genotyping the Mx dynamin like GTPase (Mx) and major histocompatibility complex class I antigen BF2 genes. These chickens were then infected with influenza A virus subtype H5N1, and their lung tissues were collected for RNA sequencing. Results: In total, 972 differentially expressed genes (DEGs) were observed between resistant and susceptible Ri chickens, according to the gene ontology and Kyoto encyclopedia of genes and genomes pathways. In particular, DEGs associated with cytokine-cytokine receptor interactions were most abundant. The expression levels of cytokines (interleukin-1β [IL-1β], IL-6, IL-8, and IL-18), chemokines (C-C Motif chemokine ligand 4 [CCL4] and CCL17), interferons (IFN-γ), and IFN-stimulated genes (Mx1, CCL19, 2'-5'-oligoadenylate synthase-like, and protein kinase R) were higher in H5N1-resistant chickens than in H5N1-susceptible chickens. Conclusion: Resistant chickens show stronger immune responses and antiviral activity (cytokines, chemokines, and IFN-stimulated genes) than those of susceptible chickens against HPAIV infection.

• Animal Bioscience
• /
• v.36 no.3
• /
• pp.404-416
• /
• 2023

Objective: Daweizi (DWZ) is a famous indigenous pig breed in China and characterized by tender meat and high fat percentage. However, the expression profiles and functions of transcripts in DWZ pigs is still in infancy. The object of this study was to depict the transcript profiles in DWZ pigs and screen the potential pathway influence adipogenesis and fat deposition, Methods: Histological analysis of backfat tissue was firstly performed between DWZ and lean-type Yorkshire pigs, and then RNA sequencing technology was utilized to explore miRNAs, lncRNAs and mRNAs profiles in backfat tissue. 18 differentially expressed (DE) transcripts were randomly selected for quantitative real-time polymerase chain reaction (QPCR) to validate the reliability of the sequencing results. Finally, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis were conducted to investigate the potential pathways influence adipocyte differentiation, adipogenesis and lipid metabolism, and a schematic model was further proposed. Results: A total of 1,625 differentially expressed transcripts were identified in DWZ pigs, including 27 upregulated and 45 downregulated miRNAs, 64 upregulated and 119 down-regulated lncRNA, 814 upregulated and 556 downregulated mRNAs. QPCR analysis exhibited strong consistency with the sequencing data. GO and KEGG analysis elucidated that the differentially expressed transcripts were mainly associated with cell growth and death, signal transduction, peroxisome proliferator-activated receptors (PPAR), AMP-activated protein kinase (AMPK), PI3K-Akt, adipocytokine and foxo signaling pathways, all of which are strongly involved in cell development, lipid metabolism and adipogenesis. Further analysis indicated that the BGIR9823_87926/miR-194a-5p/AQP7 network may be effective in the process of adipocyte differentiation or adipogenesis. Conclusion: Our study provides comprehensive insights into the regulatory network of backfat deposition and lipid metabolism in pigs from the point of view of miRNAs, lncRNAs and mRNAs.