1
|
Pareek CS, Sachajko M, Kalra G, Sultana S, Szostak A, Chalaskiewicz K, Kepka-Borkowska K, Poławska E, Ogłuszka M, Pierzchała D, Starzyński R, Taniguchi H, Juszczuk-Kubiak E, Lepczyński A, Ślaska B, Kozera W, Czarnik U, Wysocki P, Kadarmideen HN, Te Pas MFW, Szyda J, Pierzchała M. Identification of trait-associated microRNA modules in liver transcriptome of pig fed with PUFAs-enriched supplementary diet. J Appl Genet 2025; 66:389-407. [PMID: 39546271 PMCID: PMC12000271 DOI: 10.1007/s13353-024-00912-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2024] [Revised: 08/16/2024] [Accepted: 10/11/2024] [Indexed: 11/17/2024]
Abstract
Dietary lipids provide energy, are cellular structural components, and are involved in physiological processes. Lipids are the dietary source in supplementary diet experiments in pigs. This study aims to investigate the dietary effects of PUFAs on the hepatic transcriptome and physiological pathways of two diets on two pig breeds. Polish Landrace (PL: n = 6) and six PLxDuroc (PLxD: n = 6) pigs were fed with a normal diet (n = 3) or PUFAs-enriched healthy diet (n = 3), and the hepatic miRNA profiles were studied for weighted gene co-expression network analysis biological interactions between gene networks and metabolic pathways of DE miRNA genes. The study identified trait-associated modules that were significantly associated with four phenotypic traits in the dietary groups of PL and PLxD: meat colour (a*), shoulder subcutaneous fat thickness, conductivity 24 h post-mortem (PE24), and ashes. Trait-wise, a large set of co-expressed miRNAs of porcine liver were identified in these trait-associated significant modules (9, 7, 2, and 8) in PL and PLxD. Each module is represented by a module eigengene (ME). Forty-four miRNAs out of 94 miRNAs interacted with 6719 statistically significant target genes with a target score > 90. The GO/pathway analysis showed association with pathways including regulation of metallopeptidase activity, sebaceous gland development, collagen fibril organization, WNT signalling, epithelial tube morphogenesis, etc. The study showed the differences in miRNA expression between the dietary groups of PL and PLxD breeds. Hub genes of discovered miRNA clusters can be considered predicted miRNA genes associated with PE24, meat colour, shoulder subcutaneous fat thickness, and ashes. Discovered target genes for miRNA clusters play significant roles in biological functions such as (i) muscle and body growth development, (ii) different cellular processes and developments, (iii) system development, and (iv) metabolic processes.
Collapse
Affiliation(s)
- C S Pareek
- Institute of Veterinary Medicine, Faculty of Biological and Veterinary Sciences, Nicolaus Copernicus University, 87-100, Toruń, Poland
- Division of Functional Genomics in Biological and Biomedical Research, Centre for Modern Interdisciplinary Technologies, Nicolaus Copernicus University, 87-100, Torun, Poland
| | - M Sachajko
- Institute of Veterinary Medicine, Faculty of Biological and Veterinary Sciences, Nicolaus Copernicus University, 87-100, Toruń, Poland
| | - G Kalra
- Institute of Veterinary Medicine, Faculty of Biological and Veterinary Sciences, Nicolaus Copernicus University, 87-100, Toruń, Poland
| | - S Sultana
- Institute of Veterinary Medicine, Faculty of Biological and Veterinary Sciences, Nicolaus Copernicus University, 87-100, Toruń, Poland
| | - A Szostak
- Department of Genomics and Biodiversity, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
| | - K Chalaskiewicz
- Department of Genomics and Biodiversity, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
| | - K Kepka-Borkowska
- Department of Genomics and Biodiversity, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
| | - E Poławska
- Department of Genomics and Biodiversity, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
| | - M Ogłuszka
- Department of Genomics and Biodiversity, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
| | - D Pierzchała
- Maria Sklodowska-Curie National Research Institute of Oncology, W.K. Roentgena 5 Str, 02-781, Warsaw, Poland
| | - R Starzyński
- Department of Molecular Biology, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
| | - H Taniguchi
- Department of Experimental Embryology, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland
- African Genome Center, Mohammed VI Polytechnic University, UM6P, Lot 660, Hay Moulay Rachid Ben Guerir, 43150, Morocco
| | - E Juszczuk-Kubiak
- Laboratory of Biotechnology and Molecular Engineering, Department of Microbiology Prof. Wacław, Dąbrowski Institute of Agriculture and Food Biotechnology - State Research Institute (IBPRS-PIB), Rakowiecka 36 Str, 02-532, Warsaw, Poland
| | - A Lepczyński
- Department of Physiology, Cytobiology and Proteomics, West Pomeranian University of Technology, K. Janickiego 32 Str, 71-270, Szczecin, Poland
| | - B Ślaska
- Faculty of Animal Sciences and Bioeconomy, University of Life Sciences in Lublin, Akademicka 13 Str, 20-950, Lublin, Poland
| | - W Kozera
- Department of Pig Breeding, Department of Animal Biochemistry and Biotechnology, Faculty of Animal Bio-Engineering, University of Warmia and Mazury in Olsztyn, Ul. M. Oczapowskiego 5 Str, 10-719, Olsztyn, Poland
| | - U Czarnik
- Department of Pig Breeding, Department of Animal Biochemistry and Biotechnology, Faculty of Animal Bio-Engineering, University of Warmia and Mazury in Olsztyn, Ul. M. Oczapowskiego 5 Str, 10-719, Olsztyn, Poland
| | - P Wysocki
- Department of Pig Breeding, Department of Animal Biochemistry and Biotechnology, Faculty of Animal Bio-Engineering, University of Warmia and Mazury in Olsztyn, Ul. M. Oczapowskiego 5 Str, 10-719, Olsztyn, Poland
| | - H N Kadarmideen
- Department of Animal and Veterinary Sciences, Aarhus University, Blichers Alle 20, 8830, Tjele, Denmark
| | - M F W Te Pas
- Wageningen Livestock Research, Wageningen University and Research, 6708 WD, Wageningen, The Netherlands
| | - J Szyda
- Biostatistics Group, Department of Genetics, Wrocław University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wrocław, Poland
| | - M Pierzchała
- Department of Genomics and Biodiversity, Institute of Genetics and Animal Biotechnology of the Polish Academy of Sciences, Ul. Postepu 36A Str, 05-552, Jastrzebiec, Magdalenka, Poland.
| |
Collapse
|
2
|
P. D. D, M. S, K.K. J, R. U, T.V. A, G. R, V.B. SK, Asaf M, Sebastian R. Role of microRNA, bta-miR-375 in Immune Sturdiness of Vechur: The Native Cattle Breed of Kerala, India. Heliyon 2023; 9:e22683. [PMID: 38213581 PMCID: PMC10782154 DOI: 10.1016/j.heliyon.2023.e22683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2023] [Revised: 10/04/2023] [Accepted: 11/16/2023] [Indexed: 01/13/2024] Open
Abstract
In the present study, next generation sequencing was employed to identify and explore the differential expression profiles of microRNAs (miRNAs) in peripheral blood mononuclear cells (PBMCs) of crossbred (B. taurus x B. indicus) and Vechur (B. indicus) cattle in response to the bacterial endotoxin-lipopolysaccharide (LPS). The PBMCs from adult apparently healthy female crossbred cows and Vechur cattle, a native cattle breed of Kerala, India were stimulated with 10 μg/mL of LPS for 6 h. Among the differentially expressed miRNAs, the expression of 13 miRNAs showed statistically significant up regulation while, significant decrease in the expression of 15 miRNAs was noticed in LPS treated PBMCs of Vechur cattle compared to crossbred cows. The expression profiling of miRNA, bta-miR-375, expression of which was found to be significantly down regulated in LPS treated PBMCs of Vechur cattle with respect to crossbred cattle by the NGS studies, is presented in the present manuscript. The decrease in expression of bta-miR-375 noticed by NGS was in accordance with the results of quantitative real time PCR assay. Functional gene enrichment analysis and pathway analysis revealed significant enrichment of predicted targets of bta-miR-375 in many immune related and cell signalling mechanisms. In addition, over representation of targets of bta-miR-375 was also noticed in pathogenesis of many of the bovine diseases. The study could also identify differences in the expression of cytokines, viz. Tumour Necrosis Factor Alpha (TNFα), Interleukin 4 (IL-4) and Interferon-γ (IFNγ) between LPS treated and untreated PBMCs of crossbred and Vechur cattle.
Collapse
Affiliation(s)
- Divya P. D.
- Department of Veterinary Biochemistry, CVAS, Pookode, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Shynu M.
- Department of Veterinary Biochemistry, CVAS, Mannuthy, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Jayavardhanan K.K.
- Department of Veterinary Biochemistry, CVAS, Mannuthy, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Uma R.
- Department of Veterinary Biochemistry, CVAS, Mannuthy, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Aravindakshan T.V.
- CASAGB, CVAS, Mannuthy, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Radhika G.
- CASAGB, CVAS, Mannuthy, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Sameer kumar V.B.
- Department of Genomic Science, Central University of Kerala, Kasaragod, India
| | - Muhasin Asaf
- Department of Animal Breeding and Genetics, CVAS, Pookode, Kerala Veterinary and Animal Sciences University, Kerala, India
| | - Renjith Sebastian
- Department of Veterinary Biochemistry, CVAS, Pookode, Kerala Veterinary and Animal Sciences University, Kerala, India
| |
Collapse
|
3
|
Aggio-Bruce R, Schumann U, Cioanca AV, Chen FK, McLenachan S, Heath Jeffery RC, Das S, Natoli R. Serum miRNA modulations indicate changes in retinal morphology. Front Mol Neurosci 2023; 16:1130249. [PMID: 36937046 PMCID: PMC10020626 DOI: 10.3389/fnmol.2023.1130249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 02/13/2023] [Indexed: 03/06/2023] Open
Abstract
Background Age-related macular degeneration (AMD) is the leading cause of vision loss in the developed world and the detection of its onset and progression are based on retinal morphological assessments. MicroRNA (miRNA) have been explored extensively as biomarkers for a range of neurological diseases including AMD, however differences in experimental design and the complexity of human biology have resulted in little overlap between studies. Using preclinical animal models and clinical samples, this study employs a novel approach to determine a serum signature of AMD progression. Methods Serum miRNAs were extracted from mice exposed to photo-oxidative damage (PD; 0, 1, 3 and 5 days), and clinical samples from patients diagnosed with reticular pseudodrusen or atrophic AMD. The expression of ~800 miRNAs was measured using OpenArray™, and differential abundance from controls was determined using the HTqPCR R package followed by pathway analysis with DAVID. MiRNA expression changes were compared against quantifiable retinal histological indicators. Finally, the overlap of miRNA changes observed in the mouse model and human patient samples was investigated. Results Differential miRNA abundance was identified at all PD time-points and in clinical samples. Importantly, these were associated with inflammatory pathways and histological changes in the retina. Further, we were able to align findings in the mouse serum to those of clinical patients. Conclusion In conclusion, serum miRNAs are a valid tool as diagnostics for the early detection of retinal degeneration, as they reflect key changes in retinal health. The combination of pre-clinical animal models and human patient samples led to the identification of a preliminary serum miRNA signature for AMD. This study is an important platform for the future development of a diagnostic serum miRNA panel for the early detection of retinal degeneration.
Collapse
Affiliation(s)
- Riemke Aggio-Bruce
- The John Curtin School of Medical Research, The Australian National University, Acton, ACT, Australia
- The School of Medicine and Psychology, Acton, ACT, Australia
| | - Ulrike Schumann
- The John Curtin School of Medical Research, The Australian National University, Acton, ACT, Australia
- The Save Sight Institute, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW, Australia
| | - Adrian V. Cioanca
- The John Curtin School of Medical Research, The Australian National University, Acton, ACT, Australia
| | - Fred K. Chen
- Centre of Ophthalmology and Visual Science, The University of Western Australia, Perth, WA, Australia
- Lions Eye Institute, Perth, WA, Australia
- Ophthalmology, Department of Surgery, University of Melbourne, East Melbourne, VIC, Australia
- Centre for Eye Research Australia, Royal Victorian Eye and Ear Hospital, East Melbourne, VIC, Australia
| | - Samuel McLenachan
- Centre of Ophthalmology and Visual Science, The University of Western Australia, Perth, WA, Australia
- Lions Eye Institute, Perth, WA, Australia
| | - Rachael C. Heath Jeffery
- Centre of Ophthalmology and Visual Science, The University of Western Australia, Perth, WA, Australia
- Lions Eye Institute, Perth, WA, Australia
- Centre for Eye Research Australia, Royal Victorian Eye and Ear Hospital, East Melbourne, VIC, Australia
| | - Shannon Das
- The John Curtin School of Medical Research, The Australian National University, Acton, ACT, Australia
| | - Riccardo Natoli
- The John Curtin School of Medical Research, The Australian National University, Acton, ACT, Australia
- The School of Medicine and Psychology, Acton, ACT, Australia
- *Correspondence: Riccardo Natoli,
| |
Collapse
|
4
|
Wang K, Huang XT, Miao YP, Bai XL, Jin F. MiR-148a-3p attenuates apoptosis and inflammation by targeting CNTN4 in atherosclerosis. ANNALS OF TRANSLATIONAL MEDICINE 2022; 10:1201. [PMID: 36544657 PMCID: PMC9761171 DOI: 10.21037/atm-22-3768] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Accepted: 10/18/2022] [Indexed: 11/18/2022]
Abstract
Background Atherosclerosis (AS) seriously affects human health. The role of microRNAs (miRNAs) in the pathogenesis and progression of AS has become a focus of research. Our goal was to identify the biological effect of differentially expressed miRNAs (DE-miRNAs) in AS. Methods To analyze differentially expressed genes (DEGs), including differentially expressed mRNAs (DE-mRNAs) and DE-miRNAs, in AS by using the Gene Expression Omnibus (GEO) database and limma package. DEGs protein-protein interaction (PPI) network and functional enrichment analysis were constructed by using the search tool for the retrieval of interacting genes/proteins (STRING) database, Cytoscape software and Cytoscape plugin "ClueGO2.5.6". We established a coexpression network of dysregulated miRNAs and mRNAs to predict the function of miRNAs by using miRWalk database and Pearson correlation coefficient (PCC) analysis. Cellular experiments were used to validate the results of bioinformatics. Results First, 69 common DEGs were obtained from datasets GSE43292 and GSE97210 using the limma package in R. Next, a DEG PPI network was constructed. Functional enrichment analysis of DEGs showed that 11 functional pathways were significantly enriched, such as positive regulation of monocyte chemotaxis. Seven common DE-miRNAs were obtained from the GSE99685 dataset and DE-mRNAs predicted miRNAs through the miRWalk database. The miRNA-mRNA network constructed using Cytoscape software suggested that miR-148a-3p targeted contactin 4 (CNTN4). Quantitative real-time polymerase chain reaction (qRT-PCR) assay results indicated that miR-148a-3p was downregulated and CNTN4 was upregulated in the THP-1 + phorbol 12-myristate 13-acetate (PMA) + oxidized low-density lipoprotein (oxLDL) group compared with the THP-1 + PMA group. qRT-PCR, flow cytometry, and enzyme-linked immunosorbent assay (ELISA) found that upregulated miR-148a-3p significantly inhibited the expression of CNTN4, cell apoptosis, and interleukin-6 (IL-6) and tumor necrosis factor-alpha (TNF-α) concentrations in oxLDL-induced THP-1 macrophages. In addition, a dual-luciferase reporter assay demonstrated that CNTN4 was a target gene of miR-148a-3p. Conclusions Overall, these findings suggested that miR-148a-3p inhibited oxLDL-induced cell apoptosis and inflammation via targeting CNTN4 in THP-1 macrophages.
Collapse
Affiliation(s)
- Kai Wang
- Department of Neurosurgery, The Affiliated Hospital of Inner Mongolia Medical University, Hohhot, China
| | - Xi-Tong Huang
- Department of Traditional Chinese Medicine, China Pharmaceutical University, Nanjing, China
| | - Yan-Ping Miao
- Department of Radiology, The Affiliated Hospital of Inner Mongolia Medical University, Hohhot, China
| | - Xiao-Long Bai
- Department of Radiology, The Affiliated Hospital of Inner Mongolia Medical University, Hohhot, China
| | - Feng Jin
- Department of Radiology, The Affiliated Hospital of Inner Mongolia Medical University, Hohhot, China
| |
Collapse
|
5
|
Circulating miRNAs as Epigenetic Mediators of Periodontitis and Preeclampsia Association. DISEASE MARKERS 2022; 2022:2771492. [PMID: 35860693 PMCID: PMC9293528 DOI: 10.1155/2022/2771492] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Revised: 06/20/2022] [Accepted: 06/30/2022] [Indexed: 11/18/2022]
Abstract
Objective Periodontal disease has been associated with pregnancy complications including preeclampsia. This bioinformatic study is aimed at investigating the possible role of circulating microRNAs (miRNAs) as mediators of the association between maternal periodontal disease and preeclampsia. Methods Peripheral blood miRNA profiles of periodontitis and controls were sought from Gene Expression Omnibus (GEO), and differential expression analysis was performed. Experimentally validated circulating miRNAs associated with preeclampsia were determined from the Human MicroRNA Disease Database (HMDD v3.0). Venn diagrams were drawn to identify shared circulating differential miRNAs (DEmiRNAs). Significantly enriched target genes, KEGG pathways, and Gene Ontology (GO) terms for the set of shared DEmiRNA were predicted using miRNA enrichment analysis and annotation tool (miEAA v 2.0). Additionally, the shared DEmiRNA-enriched target genes were analyzed for enriched WikiPathways, BioCarta metabolic pathways, and tissue proteins in the human proteome map. Results Among 183 circulating DEmiRNA in periodontitis and 60 experimentally validated miRNA in preeclampsia, 9 shared DEmiRNA were identified. The top among 32 overrepresented target genes included MAFB, PSAP, and CDK5RAP2, top among 14 enriched KEGG pathways were renin-angiotensin system and graft-versus-host disease, and that among enriched 44 GO profiles included “positive regulation of epidermal growth factor-activated receptor activity” and “sequestering of calcium ion.” In the overrepresented target gene set, among 10 enriched WikiPathways, the top included “NAD metabolism, sirtuins, and aging” and “regulation of Wnt/B-catenin signaling by small molecule compounds” and PPAR-related mechanisms was top among 13 enriched BioCarta metabolic pathways. Conclusion A circulating 9-DEmiRNA set was significantly linked to both periodontitis and preeclampsia. Enrichment analysis identified specific genes, pathways, and functional mechanisms, which may be epigenetically altered and thereby mediate the biological association of periodontitis and preeclampsia.
Collapse
|
6
|
Zhang J, Xu X, Chen H, Kang P, Zhu H, Ren H, Liu Y. Construction and analysis for dys-regulated lncRNAs and mRNAs in LPS-induced porcine PBMCs. Innate Immun 2021; 27:170-183. [PMID: 33504244 PMCID: PMC7882806 DOI: 10.1177/1753425920983869] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Long non-coding RNAs (lncRNAs) are emerging as key regulators in inflammation. However, their functions and profiles in LPS-induced inflammation in pigs are largely unknown. In this study, we profiled global lncRNA and mRNA expression changes in PBMCs treated with LPS using the lncRNA-seq technique. In total 43 differentially expressed (DE) lncRNAs and 1082 DE mRNAs were identified in porcine PBMCs after LPS stimulation. Functional enrichment analysis on DE mRNAs indicated these genes were involved in inflammation-related signaling pathways, including cytokine–cytokine receptor interaction, TNF-α, NF-κB, Jak-STAT and TLR signaling pathways. In addition, co-expression network and function analysis identified the potential lncRNAs related to inflammatory response and immune response. The expressions of eight lncRNAs (ENSSSCT00000045208, ENSSSCT00000051636, ENSSSCT00000049770, ENSSSCT00000050966, ENSSSCT00000047491, ENSSSCT00000049750, ENSSSCT00000054262 and ENSSSCT00000044651) were validated in the LPS-treated PBMCs by quantitative real-time PCR (qRT-PCR). In LPS-challenged piglets, we identified that expression of three lncRNAs (ENSSSCT00000051636, ENSSSCT00000049770, and ENSSSCT00000047491) was significantly up-regulated in liver, spleen and jejunum tissues after LPS challenge, which indicated that these lncRNAs might be important regulators for inflammation. This study provides the first lncRNA and mRNA transcriptomic landscape of LPS-mediated changes in porcine PBMCs, which might provide potential insights into lncRNAs involved in regulating inflammation in pigs.
Collapse
Affiliation(s)
- Jing Zhang
- Hubei Key Laboratory of Animal Nutrition and Feed Science, School of Animal Science and Nutrition Engineering, Wuhan Polytechnic University, Wuhan, China
| | - Xin Xu
- Hubei Key Laboratory of Animal Nutrition and Feed Science, School of Animal Science and Nutrition Engineering, Wuhan Polytechnic University, Wuhan, China
| | - Hongbo Chen
- Hubei Key Laboratory of Animal Nutrition and Feed Science, School of Animal Science and Nutrition Engineering, Wuhan Polytechnic University, Wuhan, China
| | - Ping Kang
- Hubei Key Laboratory of Animal Nutrition and Feed Science, School of Animal Science and Nutrition Engineering, Wuhan Polytechnic University, Wuhan, China
| | - Huiling Zhu
- Hubei Key Laboratory of Animal Nutrition and Feed Science, School of Animal Science and Nutrition Engineering, Wuhan Polytechnic University, Wuhan, China
| | - Hongyan Ren
- Hubei Key Laboratory of Animal Embryo Engineering and Molecular Breeding, Hubei Institute of Animal Science and Veterinary Medicine, Hubei Academy of Agricultural Sciences, Wuhan, China
| | - Yulan Liu
- Hubei Key Laboratory of Animal Nutrition and Feed Science, School of Animal Science and Nutrition Engineering, Wuhan Polytechnic University, Wuhan, China
| |
Collapse
|