1
|
Willoughby O, Karrow NA, Marques Freire Cunha S, Asselstine V, Mallard BA, Cánovas Á. Characterization of the Hepatic Transcriptome for Divergent Immune-Responding Sheep Following Natural Exposure to Gastrointestinal Nematodes. Genes (Basel) 2024; 15:713. [PMID: 38927648 PMCID: PMC11202434 DOI: 10.3390/genes15060713] [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/24/2024] [Revised: 05/18/2024] [Accepted: 05/22/2024] [Indexed: 06/28/2024] Open
Abstract
Infections with gastrointestinal nematodes (GINs) reduce the economic efficiency of sheep operations and compromise animal welfare. Understanding the host's response to GIN infection can help producers identify animals that are naturally resistant to infection. The objective of this study was to characterize the hepatic transcriptome of sheep that had been naturally exposed to GIN parasites. The hepatic transcriptome was studied using RNA-Sequencing technology in animals characterized as high (n = 5) or medium (n = 6) based on their innate immune acute-phase (AP) response phenotype compared with uninfected controls (n = 4), and with biased antibody-mediated (AbMR, n = 5) or cell-mediated (CMR, n = 5) adaptive immune responsiveness compared to uninfected controls (n = 3). Following the assessment of sheep selected for innate responses, 0, 136, and 167 genes were differentially expressed (DE) between high- and medium-responding animals, high-responding and uninfected control animals, and medium-responding and uninfected control animals, respectively (false discovery rate (FDR) < 0.05, and fold change |FC| > 2). When adaptive immune responses were assessed, 0, 53, and 57 genes were DE between antibody- and cell-biased animals, antibody-biased and uninfected control animals, and cell-biased and uninfected control animals, respectively (FDR < 0.05, |FC| > 2). Functional analyses identified enriched gene ontology (GO) terms and metabolic pathways related to the innate immune response and energy metabolism. Six functional candidate genes were identified for further functional and validation studies to better understand the underlying biological mechanisms of host responses to GINs. These, in turn, can potentially help improve decision making and management practices to increase the overall host immune response to GIN infection.
Collapse
Affiliation(s)
- Olivia Willoughby
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, Ontario Agriculture College, University of Guelph, 50 Stone Road E, Guelph, ON N1G 2W1, Canada; (O.W.); (N.A.K.); (S.M.F.C.); (V.A.)
| | - Niel A. Karrow
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, Ontario Agriculture College, University of Guelph, 50 Stone Road E, Guelph, ON N1G 2W1, Canada; (O.W.); (N.A.K.); (S.M.F.C.); (V.A.)
| | - Samla Marques Freire Cunha
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, Ontario Agriculture College, University of Guelph, 50 Stone Road E, Guelph, ON N1G 2W1, Canada; (O.W.); (N.A.K.); (S.M.F.C.); (V.A.)
| | - Victoria Asselstine
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, Ontario Agriculture College, University of Guelph, 50 Stone Road E, Guelph, ON N1G 2W1, Canada; (O.W.); (N.A.K.); (S.M.F.C.); (V.A.)
| | - Bonnie A. Mallard
- Department of Pathobiology, Ontario Veterinary College, University of Guelph, 50 Stone Road E, Guelph, ON N1G 2W1, Canada;
| | - Ángela Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, Ontario Agriculture College, University of Guelph, 50 Stone Road E, Guelph, ON N1G 2W1, Canada; (O.W.); (N.A.K.); (S.M.F.C.); (V.A.)
| |
Collapse
|
2
|
A V, Kumar A, Mahala S, Chandra Janga S, Chauhan A, Mehrotra A, Kumar De A, Ranjan Sahu A, Firdous Ahmad S, Vempadapu V, Dutt T. Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds. Gene 2024; 893:147950. [PMID: 37918549 DOI: 10.1016/j.gene.2023.147950] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 10/16/2023] [Accepted: 10/30/2023] [Indexed: 11/04/2023]
Abstract
In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.
Collapse
Affiliation(s)
- Vani A
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, UP, India
| | - Amit Kumar
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, UP, India.
| | - Sudarshan Mahala
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, UP, India
| | - Sarath Chandra Janga
- Luddy School of Informatics, Computing, and Engineering, Indiana University, IUPUI, Indianapolis, IN, USA
| | - Anuj Chauhan
- Livestock Production and Management, Indian Veterinary Research Institute, Bareilly, UP, India
| | | | - Arun Kumar De
- Central Island Agricultural Research Institute, Port Blair, Andaman and Nicobar Islands, India
| | - Amiya Ranjan Sahu
- Central Coastal Agricultural Research Institute, Old Goa, Goa, India
| | - Sheikh Firdous Ahmad
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, UP, India
| | - Varshini Vempadapu
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, UP, India
| | - Triveni Dutt
- Livestock Production and Management, Indian Veterinary Research Institute, Bareilly, UP, India
| |
Collapse
|
3
|
Perdomo-González DI, Id-Lahoucine S, Molina A, Cánovas A, Laseca N, Azor PJ, Valera M. Transmission ratio distortion detection by neutral genetic markers in the Pura Raza Española horse breed. Animal 2023; 17:101012. [PMID: 37950978 DOI: 10.1016/j.animal.2023.101012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 10/05/2023] [Accepted: 10/06/2023] [Indexed: 11/13/2023] Open
Abstract
Transmission Ratio Distortion (TRD) is a genetic phenomenon widely demonstrated in several livestock species, but barely in equine species. The TRD occurs when certain genotypes are over- or under-represented in the offspring of a particular mating and can be caused by a variety of factors during gamete formation or during embryonic development. For this study, 126 394 trios consisting of a stallion, mare, and offspring were genotyped using a panel of 17 neutral microsatellite markers recommended by the International Society for Animal Genetics for paternity tests and individual identification. The number of alleles available for each marker ranges from 13 to 18, been 268 the total number of alleles investigated. The TRDscan v.2.0 software was used with the biallelic procedure to identify regions with distorted segregation ratios. After completing the analysis, a total of 12 alleles (out of 11 microsatellites) were identified with decisive evidence for genotypic TRD; 3 and 9 with additive and heterosis patterns, respectively. In addition, 19 alleles (out of 10 microsatellites) were identified displaying allelic TRD. Among them, 14 and 5 were parent-unspecific and stallion-mare-specific TRD. Out of the TRD regions, 24 genes were identified and annotated, predominantly associated with cholesterol metabolism and homeostasis. These genes are often linked to non-specific symptoms like impaired fertility, stunted growth, and compromised overall health. The results suggest a significant impact on the inheritance of certain genetic traits in horses. Further analysis and validation are needed to better understand the TRD impact before the potential implementation in the horse breeding programme strategies.
Collapse
Affiliation(s)
| | - S Id-Lahoucine
- Department of Animal and Veterinary Science, Scotland's Rural College, Easter Bush, Edinburgh EH25 9RG, United Kingdom
| | - A Molina
- Departamento de Genética, Universidad de Córdoba, Córdoba 14014, Spain
| | - A Cánovas
- Center of Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - N Laseca
- Departamento de Genética, Universidad de Córdoba, Córdoba 14014, Spain
| | - P J Azor
- Real Asociación Nacional de Criadores de Caballos de Pura Raza Española (ANCCE), Sevilla 41014, Spain
| | - M Valera
- Departamento de Agronomía, ETSIA, Universidad de Sevilla, Sevilla 41005, Spain
| |
Collapse
|
4
|
Fonseca PAS, Lam S, Chen Y, Waters SM, Guan LL, Cánovas A. Multi-breed host rumen epithelium transcriptome and microbiome associations and their relationship with beef cattle feed efficiency. Sci Rep 2023; 13:16209. [PMID: 37758745 PMCID: PMC10533831 DOI: 10.1038/s41598-023-43097-8] [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: 12/15/2022] [Accepted: 09/19/2023] [Indexed: 09/29/2023] Open
Abstract
Understanding host-microbial interactions in the rumen and its influence on desirable production traits may lead to potential microbiota manipulation or genetic selection for improved cattle feed efficiency. This study investigated the host transcriptome and its correlation with the rumen archaea and bacteria differential abundance of two pure beef cattle breeds (Angus and Charolais) and one composite beef hybrid (Kinsella) divergent for residual feed intake (RFI; low-RFI vs. high-RFI). Using RNA-Sequencing of rumen tissue and 16S rRNA gene amplicon sequencing, differentially expressed genes (FDR ≤ 0.05, |log2(Fold-change) >|2) and differentially abundant (p-value < 0.05) archaea and bacteria amplicon sequence variants (ASV) were determined. Significant correlations between gene expression and ASVs (p-value < 0.05) were determine using Spearman correlation. Interesting associations with muscle contraction and the modulation of the immune system were observed for the genes correlated with bacterial ASVs. Potential functional candidate genes for feed efficiency status were identified for Angus (CCL17, CCR3, and CXCL10), Charolais (KCNK9, GGT1 and IL6), and Kinsella breed (ESR2). The results obtained here provide more insights regarding the applicability of target host and rumen microbial traits for the selection and breeding of more feed efficient beef cattle.
Collapse
Grants
- Beef Farmers of Ontario, Genome Canada and the Sustainable Beef and Forage Science Cluster funded by the Canadian Beef Cattle Check-Off, Beef Cattle Research Council (BCRC), Alberta Beef Producers, Alberta Cattle Feeders’ Association, Beef Farmers of Ontario, La Fédération des Productuers de bovins du Québec, and Agriculture and Agri-Food Canada’s Canadian Agricultural Partnership
- Ontario Ministry of Agriculture, Food, and Rural Affairs (OMAFRA), Ontario Ministry of Research and Innovation, and the Ontario Agri-Food Innovation Alliance
Collapse
Affiliation(s)
- P A S Fonseca
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - S Lam
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Y Chen
- Livestock Gentec, Department of Agriculture, Food & Nutritional Science, University of Alberta, Edmonton, AB, T6H 2P5, Canada
| | - S M Waters
- Teagasc, Animal and Bioscience Research Department, Animal and Grassland Research and Innovation Centre, Grange, Dunsany, C15 PW93, Co. Meath, Ireland
| | - L L Guan
- Livestock Gentec, Department of Agriculture, Food & Nutritional Science, University of Alberta, Edmonton, AB, T6H 2P5, Canada
| | - A Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
5
|
Nishihara K, van Niekerk J, Innes D, He Z, Cánovas A, Guan LL, Steele M. Transcriptome profiling revealed that key rumen epithelium functions change in relation to short-chain fatty acids and rumen epithelium-attached microbiota during the weaning transition. Genomics 2023; 115:110664. [PMID: 37286013 DOI: 10.1016/j.ygeno.2023.110664] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 05/12/2023] [Accepted: 05/31/2023] [Indexed: 06/09/2023]
Abstract
This study aims to characterize the functional changes of the rumen epithelium associated with ruminal short-chain fatty acid (SCFA) concentration and epithelium-attached microbes during the weaning transition in dairy calves. Ruminal SCFA concentrations were determined, and transcriptome and microbiota profiling in biopsied rumen papillae were obtained from Holstein calves before and after weaning using RNA- and amplicon sequencing. Metabolic pathway analysis showed that pathways related to SCFA metabolism and cell apoptosis were up- and down-regulated postweaning, respectively. Functional analysis showed that genes related to SCFA absorption, metabolism, and protective roles against oxidative stress were positively correlated with ruminal SCFA concentrations. The relative abundance of epithelium-attached Rikenellaceae RC9 gut group and Campylobacter was positively correlated with genes involved in SCFA absorption and metabolism, suggesting that these microbes can cooperatively affect host functions. Future research should examine the contribution of attenuated apoptosis on rumen epithelial functional shifts during the weaning transition.
Collapse
Affiliation(s)
- Koki Nishihara
- Department of Animal Biosciences, Animal Science and Nutrition, University of Guelph, Guelph, ON N1G 1Y2, Canada
| | - Jolet van Niekerk
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2P5, Canada
| | - David Innes
- Department of Animal Biosciences, Animal Science and Nutrition, University of Guelph, Guelph, ON N1G 1Y2, Canada
| | - Zhixiong He
- CAS Key Laboratory for Agro-Ecological Processes in Subtropical Region, National Engineering Laboratory for Pollution Control and Waste Utilization in Livestock and Poultry Production, South-Central Experimental Station of Animal Nutrition and Feed Science in Ministry of Agriculture, Hunan Provincial Engineering Research Center for Healthy Livestock and Poultry Production, Institute of Subtropical Agriculture, The Chinese Academy of Sciences, Changsha 410125, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - Le Luo Guan
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2P5, Canada
| | - Michael Steele
- Department of Animal Biosciences, Animal Science and Nutrition, University of Guelph, Guelph, ON N1G 1Y2, Canada.
| |
Collapse
|
6
|
Nishihara K, van Niekerk J, He Z, Innes D, Guan LL, Steele M. Reduction in mucosa thickness is associated with changes in immune function in the colon mucosa during the weaning transition in Holstein bull dairy calves. Genomics 2023; 115:110680. [PMID: 37454938 DOI: 10.1016/j.ygeno.2023.110680] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Revised: 07/06/2023] [Accepted: 07/12/2023] [Indexed: 07/18/2023]
Abstract
This study aims to characterize changes in the structure and the molecules related to immune function in the colon mucosa in dairy calves during the weaning transition (weaned at week 6 of age). Colon mucosa thickness, measured at week 5 to 8 and 12 of age, decreased for 2 weeks after weaning, but then recovered. Colon mucosa's transcriptome profiling at week 5, 7, and 12 of age was obtained using RNA-sequencing. Functional analysis showed that pathways related to immune function were up-regulated postweaning. A weighted gene co-expression network analysis identified 17 immune function related genes, expressed higher postweaning, which were negatively correlated with colon mucosa thickness, suggesting that these genes may be involved in colon mucosa inflammation and recovery from mucosa thickness decrement during the weaning transition. As such, it is important to determine the function of immune cells in the colon mucosa during the weaning transition in dairy calves.
Collapse
Affiliation(s)
- Koki Nishihara
- Department of Animal Biosciences, Animal Science and Nutrition, University of Guelph, Guelph, ON N1G 1Y2, Canada
| | - Jolet van Niekerk
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2P5, Canada
| | - Zhixiong He
- CAS Key Laboratory for Agro-Ecological Processes in Subtropical Region, Institute of Subtropical Agriculture, The Chinese Academy of Sciences, Changsha 410125, China
| | - David Innes
- Department of Animal Biosciences, Animal Science and Nutrition, University of Guelph, Guelph, ON N1G 1Y2, Canada
| | - Le Luo Guan
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2P5, Canada
| | - Michael Steele
- Department of Animal Biosciences, Animal Science and Nutrition, University of Guelph, Guelph, ON N1G 1Y2, Canada.
| |
Collapse
|
7
|
Laseca N, Cánovas Á, Valera M, Id-Lahoucine S, Perdomo-González DI, Fonseca PAS, Demyda-Peyrás S, Molina A. Genomic screening of allelic and genotypic transmission ratio distortion in horse. PLoS One 2023; 18:e0289066. [PMID: 37556504 PMCID: PMC10411798 DOI: 10.1371/journal.pone.0289066] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Accepted: 07/11/2023] [Indexed: 08/11/2023] Open
Abstract
The phenomenon in which the expected Mendelian inheritance is altered is known as transmission ratio distortion (TRD). The TRD analysis relies on the study of the transmission of one of the two alleles from a heterozygous parent to the offspring. These distortions are due to biological mechanisms affecting gametogenesis, embryo development and/or postnatal viability, among others. In this study, TRD phenomenon was characterized in horses using SNP-by-SNP model by TRDscan v.2.0 software. A total of 1,041 Pura Raza Español breed horses were genotyped with 554,634 SNPs. Among them, 277 horses genotyped in trios (stallion-mare-offspring) were used to perform the TRD analysis. Our results revealed 140 and 42 SNPs with allelic and genotypic patterns, respectively. Among them, 63 displayed stallion-TRD and 41 exhibited mare-TRD, while 36 SNPs showed overall TRD. In addition, 42 SNPs exhibited heterosis pattern. Functional analyses revealed that the annotated genes located within the TRD regions identified were associated with biological processes and molecular functions related to spermatogenesis, oocyte division, embryonic development, and hormonal activity. A total of 10 functional candidate genes related to fertility were found. To our knowledge, this is the most extensive study performed to evaluate the presence of alleles and functional candidate genes with transmission ratio distortion affecting reproductive performance in the domestic horse.
Collapse
Affiliation(s)
- Nora Laseca
- Department of Genetics, University of Cordoba, Córdoba, Spain
| | - Ángela Cánovas
- Center of Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, Ontario, Canada
| | - Mercedes Valera
- Department of Agronomy, School of Agronomy Engineering, University of Seville, Seville, Spain
| | - Samir Id-Lahoucine
- Department of Animal and Veterinary Science, Scotland’s Rural College, Aberdeen, Scotland, United Kingdom
| | | | | | | | - Antonio Molina
- Department of Genetics, University of Cordoba, Córdoba, Spain
| |
Collapse
|
8
|
Siurana A, Cánovas A, Casellas J, Calsamiglia S. Transcriptome Profile in Dairy Cows Resistant or Sensitive to Milk Fat Depression. Animals (Basel) 2023; 13:ani13071199. [PMID: 37048455 PMCID: PMC10093643 DOI: 10.3390/ani13071199] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 03/26/2023] [Accepted: 03/27/2023] [Indexed: 04/01/2023] Open
Abstract
Feeding linseed to dairy cows results in milk fat depression (MFD), but there is a wide range of sensitivity among cows. The objectives of this study were to identify target genes containing SNP that may play a key role in the regulation of milk fat synthesis in cows resistant or sensitive to MFD. Four cows were selected from a dairy farm after a switch from a control diet to a linseed-rich diet; two were resistant to MFD with a high milk fat content in the control (4.06%) and linseed-rich (3.90%) diets; and two were sensitive to MFD with the milk fat content decreasing after the change from the control (3.87%) to linseed-rich (2.52%) diets. Transcriptome and SNP discovery analyses were performed using RNA-sequencing technology. There was a large number of differentially expressed genes in the control (n = 1316) and linseed-rich (n = 1888) diets. Of these, 15 genes were detected as key gene regulators and harboring SNP in the linseed-rich diet. The selected genes MTOR, PDPK1, EREG, NOTCH1, ZNF217 and TGFB3 may form a network with a principal axis PI3K/Akt/MTOR/SREBP1 involved in milk fat synthesis and in the response to diets that induced MFD. These 15 genes are novel candidate genes to be involved in the resistance or sensitivity of dairy cows to milk fat depression.
Collapse
|
9
|
Identification of Novel mRNA Isoforms Associated with Acute Heat Stress Response Using RNA Sequencing Data in Sprague Dawley Rats. BIOLOGY 2022; 11:biology11121740. [PMID: 36552250 PMCID: PMC9774719 DOI: 10.3390/biology11121740] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 11/17/2022] [Accepted: 11/21/2022] [Indexed: 12/02/2022]
Abstract
The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This study identified the differentially expressed mRNA isoforms which narrowed down the most reliable DEG markers and molecular pathways that underlie the mechanisms of thermoregulation. This experiment was performed on Sprague Dawley rats housed at 22 °C (control group; CT), and three acute heat-stressed groups housed at 42 °C for 30 min (H30), 60 min (H60), and 120 min (H120). Earlier, we demonstrated that acute heat stress increased the rectal temperature of rats, caused abnormal changes in the blood biochemical parameters, as well as induced dramatic changes in the expression levels of genes through epigenetics and post-transcriptional regulation. Transcriptomic analysis using RNA-Sequencing (RNA-Seq) data obtained previously from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5764, and 4988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3549 and 120 in the liver, as well as 3078 and 220 in adrenal glands, respectively. About 35 genes were involved in the heat stress response, out of which, Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3, and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between the liver and adrenal glands. In conclusion, this study would enhance our understanding of the complex underlying mechanisms of acute heat stress, and the identified mRNA isoforms and genes can be used as potential candidates for thermotolerance selection in mammals.
Collapse
|
10
|
He J, Zhao B, Huang X, Fu X, Liu G, Tian Y, Wu C, Mao J, Liu J, Gun S, Tian K. Gene network analysis reveals candidate genes related with the hair follicle development in sheep. BMC Genomics 2022; 23:428. [PMID: 35672687 PMCID: PMC9175362 DOI: 10.1186/s12864-022-08552-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Accepted: 04/13/2022] [Indexed: 12/13/2022] Open
Abstract
Background Merino sheep are the most famous fine wool sheep in the world. They have high wool production and excellent wool quality and have attracted worldwide attention. The fleece of the Merino sheep is composed predominantly of wool fibers grown from secondary wool follicles. Therefore, it is necessary to study the development of hair follicles to understand the mechanism of wool production. The hair follicle is a complex biological system involved in a dynamic process governed by gene regulation. The hair follicle development process is very complex and poorly understood. The purpose of our research is to identify candidate genes related to hair follicle development, provide a theoretical molecular breeding basis for the cultivation of fine wool sheep, and provide a reference for the problems of hair loss and alopecia areata that affect human beings. Results We analyzed mRNAs data in skin tissues of 18 Merino sheep at four embryonic days (E65, E85, E105 and E135) and two postnatal days (P7 and P30). G1 to G6 represent hair follicles developmental at six stages (i.e. E65 to P30). We identified 7879 differentially expressed genes (DEGs) and 12623 novel DEGs, revealed different expression patterns of these DEGs at six stages of hair follicle development, and demonstrated their complex interactions. DEGs with stage-specific expression were significantly enriched in epidermal differentiation and development, hair follicle development and hair follicle morphogenesis and were enriched in many pathways related to hair follicle development. The key genes (LAMA5, WNT10A, KRT25, SOSTDC1, ZDHHC21, FZD1, BMP7, LRP4, TGFβ2, TMEM79, SOX10, ITGB4, KRT14, ITGA6, and GLI2) affecting hair follicle morphogenesis were identified by network analysis. Conclusion This study provides a new reference for the molecular basis of hair follicle development and lays a foundation for further improving sheep hair follicle breeding. Candidate genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine wool sheep. These results are a valuable resource for biological investigations of fleece evolution in animals. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08552-2.
Collapse
Affiliation(s)
- Junmin He
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Bingru Zhao
- College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Xixia Huang
- College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| | - Xuefeng Fu
- Key Laboratory of Genetics Breeding and Reproduction of the Fine Wool Sheep & Cashmere Goat in Xinjiang, Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi, China
| | - Guifen Liu
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Yuezhen Tian
- Key Laboratory of Genetics Breeding and Reproduction of the Fine Wool Sheep & Cashmere Goat in Xinjiang, Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi, China
| | - Cuiling Wu
- College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| | - Jingyi Mao
- College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| | - Jing Liu
- College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| | - Shuangbao Gun
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China.
| | - Kechuan Tian
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, Jinan, China.
| |
Collapse
|
11
|
Asselstine V, Medrano JF, Cánovas A. Identification of novel alternative splicing associated with mastitis disease in Holstein dairy cows using large gap read mapping. BMC Genomics 2022; 23:222. [PMID: 35305573 PMCID: PMC8934477 DOI: 10.1186/s12864-022-08430-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Accepted: 02/24/2022] [Indexed: 12/16/2022] Open
Abstract
Background Mastitis is a very common disease in the dairy industry that producers encounter daily. Transcriptomics, using RNA-Sequencing (RNA - Seq) technology, can be used to study the functional aspect of mastitis resistance to identify animals that have a better immune response to mastitis. When the cow has mastitis, not only genes but also specific mRNA isoforms generated via alternative splicing (AS) could be differentially expressed (DE), leading to the phenotypic variation observed. Therefore, the objective of this study was to use large gap read mapping to identify mRNA isoforms DE between healthy and mastitic milk somatic cell samples (N = 12). These mRNA isoforms were then categorized based on being 1) annotated mRNA isoforms for gene name and length, 2) annotated mRNA isoforms with different transcript length and 3) novel mRNA isoforms of non - annotated genes. Results Analysis identified 333 DE transcripts (with at least 2 mRNA isoforms annotated, with at least one being DE) between healthy and mastitic samples corresponding to 303 unique genes. Of these 333 DE transcripts between healthy and mastitic samples, 68 mRNA isoforms are annotated in the bovine genome reference (ARS.UCD.1.2), 249 mRNA isoforms had novel transcript lengths of known genes and 16 were novel transcript lengths of non - annotated genes in the bovine genome reference (ARS.UCD.1.2). Functional analysis including gene ontology, gene network and metabolic pathway analysis was performed on the list of 288 annotated and unique DE mRNA isoforms. In total, 67 significant metabolic pathways were identified including positive regulation of cytokine secretion and immune response. Additionally, numerous DE novel mRNA isoforms showed potential involvement with the immune system or mastitis. Lastly, QTL annotation analysis was performed on coding regions of the DE mRNA isoforms, identifying overlapping QTLs associated with clinical mastitis and somatic cell score. Conclusion This study identified novel mRNA isoforms generated via AS that could lead to differences in the immune response of Holstein dairy cows and be potentially implemented in future breeding programs. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08430-x.
Collapse
|
12
|
Fonseca PAS, Schenkel FS, Cánovas A. Genome-wide association study using haplotype libraries and repeated measures model to identify candidate genomic regions for stillbirth in Holstein cattle. J Dairy Sci 2022; 105:1314-1326. [PMID: 34998559 DOI: 10.3168/jds.2021-20936] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 09/24/2021] [Indexed: 11/19/2022]
Abstract
Reduced fertility is one of the main causes of economic losses on dairy farms, resulting in economic losses estimated at $938 per stillbirth case in Holstein herds. The identification of genomic regions associated with stillbirth could help to develop better management and breeding strategies aimed to reduce the frequency of undesirable gestation outcomes. Here, 10,570 cows and 50,541 birth records were used to perform a haplotype-based GWAS. A total of 41 significantly associated pseudo-SNPs (haplotypes within haplotype blocks converted to a binary classification) were identified after Bonferroni adjustment for multiple tests. A total of 117 positional candidate genes were annotated within or close (in a 200-kb interval) to significant pseudo-SNPs (haplotype blocks). The guilt-by-association functional prioritization identified 31 potential functional candidate genes for reproductive performance out of the 117 positional candidate genes annotated. These genes play crucial roles in biological processes associated with pregnancy persistence, fetus development, immune response, among others. These results helped us to better understand the genetic basis of stillbirth in dairy cattle and may be useful for the prediction of stillbirth in Holstein cattle, helping to reduce the related economic losses caused by this phenotype.
Collapse
Affiliation(s)
- P A S Fonseca
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - F S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - A Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada.
| |
Collapse
|
13
|
Dou J, Cánovas A, Brito LF, Yu Y, Schenkel FS, Wang Y. Comprehensive RNA-Seq Profiling Reveals Temporal and Tissue-Specific Changes in Gene Expression in Sprague-Dawley Rats as Response to Heat Stress Challenges. Front Genet 2021; 12:651979. [PMID: 33897767 PMCID: PMC8063118 DOI: 10.3389/fgene.2021.651979] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Accepted: 03/02/2021] [Indexed: 12/13/2022] Open
Abstract
Understanding heat stress physiology and identifying reliable biomarkers are paramount for developing effective management and mitigation strategies. However, little is known about the molecular mechanisms underlying thermal tolerance in animals. In an experimental model of Sprague–Dawley rats subjected to temperatures of 22 ± 1°C (control group; CT) and 42°C for 30 min (H30), 60 min (H60), and 120 min (H120), RNA-sequencing (RNA-Seq) assays were performed for blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120). A total of 53, 1,310, and 1,501 differentially expressed genes (DEGs) were significantly identified in the blood (P < 0.05 and |fold change (FC)| >2), liver (P < 0.01, false discovery rate (FDR)–adjusted P = 0.05 and |FC| >2) and adrenal glands (P < 0.01, FDR-adjusted P = 0.05 and |FC| >2), respectively. Of these, four DEGs, namely Junb, P4ha1, Chordc1, and RT1-Bb, were shared among the three tissues in CT vs. H120 comparison. Functional enrichment analyses of the DEGs identified in the blood (CT vs. H120) revealed 12 biological processes (BPs) and 25 metabolic pathways significantly enriched (FDR = 0.05). In the liver, 133 BPs and three metabolic pathways were significantly detected by comparing CT vs. H30, H60, and H120. Furthermore, 237 BPs were significantly (FDR = 0.05) enriched in the adrenal glands, and no shared metabolic pathways were detected among the different heat-stressed groups of rats. Five and four expression patterns (P < 0.05) were uncovered by 73 and 91 shared DEGs in the liver and adrenal glands, respectively, over the different comparisons. Among these, 69 and 73 genes, respectively, were proposed as candidates for regulating heat stress response in rats. Finally, together with genome-wide association study (GWAS) results in cattle and phenome-wide association studies (PheWAS) analysis in humans, five genes (Slco1b2, Clu, Arntl, Fads1, and Npas2) were considered as being associated with heat stress response across mammal species. The datasets and findings of this study will contribute to a better understanding of heat stress response in mammals and to the development of effective approaches to mitigate heat stress response in livestock through breeding.
Collapse
Affiliation(s)
- Jinhuan Dou
- Key Laboratory of Animal Genetics, Breeding and Reproduction, MARA, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China.,Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Angela Cánovas
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | - Ying Yu
- Key Laboratory of Animal Genetics, Breeding and Reproduction, MARA, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Flavio S Schenkel
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Yachun Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction, MARA, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| |
Collapse
|
14
|
Liu Y, Muniz MMM, Lam S, Song D, Zhang Y, Yin J, Cánovas A, Liu H. Gene expression profile of the taimen Hucho taimen in response to acute temperature changes. COMPARATIVE BIOCHEMISTRY AND PHYSIOLOGY D-GENOMICS & PROTEOMICS 2021; 38:100824. [PMID: 33743513 DOI: 10.1016/j.cbd.2021.100824] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2020] [Revised: 03/04/2021] [Accepted: 03/04/2021] [Indexed: 10/21/2022]
Abstract
The endangered cold-water fish species taimen (Hucho taimen) suffer acute temperature changes in culture and wild conditions. Understanding the effects of acute temperature changes on physiological processes of this species is essential for aquaculture practices and conservation. Liver transcriptomic profiles of the taimen (n = 24) exposed to acute temperature decrease (from 20 °C to 10 °C) and acute temperature increase (from 10 °C to 20 °C) was evaluated using high-throughput RNA-Sequencing. Samples were collected at day 0, 1, 7 and 35 in both treatments. Compared to day 0, the total numbers of differentially expressed genes (DEGs) in the taimen after acute temperature decrease were 173, 226 and 42 at day 1, 7 and 35, respectively, and the total numbers of DEGs following acute temperature increase were 260, 253 and 282 at day 1, 7 and 35, respectively. Particularly, 14 key regulatory genes were commonly found between both acute temperature treatments. Functional analysis based on the commonly identified DEGs revealed important metabolic pathways related to metabolism and immune function, suggesting a specific response mechanism of taimen against cold and heat shock. The results may assist in developing management strategies for stress mediation caused by acute temperature changes in the taimen and other cold water fish.
Collapse
Affiliation(s)
- Yang Liu
- Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin 150070, China; Key Laboratory of Aquatic Animal Diseases and Immune Technology of Heilongjiang Province, Harbin 150070, China
| | - Maria Malane Magalhães Muniz
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph N1G2W1, Canada
| | - Stephanie Lam
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph N1G2W1, Canada
| | - Dan Song
- Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin 150070, China
| | - Yongquan Zhang
- Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin 150070, China
| | - Jiasheng Yin
- Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin 150070, China
| | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph N1G2W1, Canada.
| | - Hongbai Liu
- Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin 150070, China; Key Laboratory of Aquatic Animal Diseases and Immune Technology of Heilongjiang Province, Harbin 150070, China.
| |
Collapse
|
15
|
Lam S, Miglior F, Fonseca PAS, Gómez-Redondo I, Zeidan J, Suárez-Vega A, Schenkel F, Guan LL, Waters S, Stothard P, Cánovas A. Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing. J Dairy Sci 2020; 104:1928-1950. [PMID: 33358171 DOI: 10.3168/jds.2020-18241] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Accepted: 08/29/2020] [Indexed: 12/15/2022]
Abstract
The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.
Collapse
Affiliation(s)
- S Lam
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - F Miglior
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - P A S Fonseca
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - I Gómez-Redondo
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - J Zeidan
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - A Suárez-Vega
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - F Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1
| | - L L Guan
- Department of Agriculture, Food & Nutritional Science, University of Alberta, Edmonton, Canada T6H 2P5
| | - S Waters
- Teagasc, Animal & Grassland Research and Innovation Centre, Grange, Dunsany, Ireland C15 PW93
| | - P Stothard
- Department of Agriculture, Food & Nutritional Science, University of Alberta, Edmonton, Canada T6H 2P5
| | - A Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada N1G 2W1.
| |
Collapse
|
16
|
Sweett H, Fonseca PAS, Suárez-Vega A, Livernois A, Miglior F, Cánovas A. Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle. Sci Rep 2020; 10:20102. [PMID: 33208801 PMCID: PMC7676258 DOI: 10.1038/s41598-020-75758-3] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Accepted: 10/16/2020] [Indexed: 12/20/2022] Open
Abstract
Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.
Collapse
Affiliation(s)
- H Sweett
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - P A S Fonseca
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - A Suárez-Vega
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - A Livernois
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada.,Department of Pathobiology, Ontario Veterinary College, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - F Miglior
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - A Cánovas
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
17
|
Cai C, Li M, Zhang Y, Meng S, Yang Y, Gao P, Guo X, Cao G, Li B. Comparative Transcriptome Analyses of Longissimus thoracis Between Pig Breeds Differing in Muscle Characteristics. Front Genet 2020; 11:526309. [PMID: 33329687 PMCID: PMC7717936 DOI: 10.3389/fgene.2020.526309] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2020] [Accepted: 10/15/2020] [Indexed: 01/23/2023] Open
Abstract
The two breeds, Mashen (MS; a northern China breed) and Large White (LW; a western lean breed) pigs, show important phenotypic differences in growth, adaptability, intramuscular fat (IMF) content, and energy metabolism since early developmental stage. The main aim of this study was the evaluation of longissimus thoracis muscle transcriptome profile of both genetic types to identify genes, pathways responsible for their differentiated phenotype. Longissimus thoracis of MS and LW pigs were sampled at 0 day (early stage), 90 days (middle stage), and 180 days (later stage) after birth. A total of 3,487 differentially expressed genes (DEGs) were identified at the three time points. Sample clustering analysis revealed the slower growth rate of MS than LW pigs. Gene expression pattern analysis revealed that multicellular organism growth genes (GHSR, EZR, FOXS1, DRD2, SH3PXD2B, CSF1, and TSHR) were involved in the fast growth rate of LW pigs at early stage. Furthermore, DEGs (ACACA, ACSF3, OXSM, CBR4, and HSD17B8) functionally related to fatty acid synthesis revealed that IMF accumulation occurred around 90 and up to 180 days. These DEGs provided valuable resource to study the phenotypic difference in longissimus thoracis muscle between MS and LW pigs.
Collapse
Affiliation(s)
- Chunbo Cai
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Meng Li
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Yanwei Zhang
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Shan Meng
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Yang Yang
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Pengfei Gao
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Xiaohong Guo
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Guoqing Cao
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| | - Bugao Li
- College of Animal Science, Shanxi Agricultural University, Shanxi, China
| |
Collapse
|
18
|
Lam S, Zeidan J, Miglior F, Suárez-Vega A, Gómez-Redondo I, Fonseca PAS, Guan LL, Waters S, Cánovas A. Development and comparison of RNA-sequencing pipelines for more accurate SNP identification: practical example of functional SNP detection associated with feed efficiency in Nellore beef cattle. BMC Genomics 2020; 21:703. [PMID: 33032519 PMCID: PMC7545862 DOI: 10.1186/s12864-020-07107-7] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Accepted: 09/28/2020] [Indexed: 12/14/2022] Open
Abstract
Background Optimization of an RNA-Sequencing (RNA-Seq) pipeline is critical to maximize power and accuracy to identify genetic variants, including SNPs, which may serve as genetic markers to select for feed efficiency, leading to economic benefits for beef production. This study used RNA-Seq data (GEO Accession ID: PRJEB7696 and PRJEB15314) from muscle and liver tissue, respectively, from 12 Nellore beef steers selected from 585 steers with residual feed intake measures (RFI; n = 6 low-RFI, n = 6 high-RFI). Three RNA-Seq pipelines were compared including multi-sample calling from i) non-merged samples; ii) merged samples by RFI group, iii) merged samples by RFI and tissue group. The RNA-Seq reads were aligned against the UMD3.1 bovine reference genome (release 94) assembly using STAR aligner. Variants were called using BCFtools and variant effect prediction (VeP) and functional annotation (ToppGene) analyses were performed. Results On average, total reads detected for Approach i) non-merged samples for liver and muscle, were 18,362,086.3 and 35,645,898.7, respectively. For Approach ii), merging samples by RFI group, total reads detected for each merged group was 162,030,705, and for Approach iii), merging samples by RFI group and tissues, was 324,061,410, revealing the highest read depth for Approach iii). Additionally, Approach iii) merging samples by RFI group and tissues, revealed the highest read depth per variant coverage (572.59 ± 3993.11) and encompassed the majority of localized positional genes detected by each approach. This suggests Approach iii) had optimized detection power, read depth, and accuracy of SNP calling, therefore increasing confidence of variant detection and reducing false positive detection. Approach iii) was then used to detect unique SNPs fixed within low- (12,145) and high-RFI (14,663) groups. Functional annotation of SNPs revealed positional candidate genes, for each RFI group (2886 for low-RFI, 3075 for high-RFI), which were significantly (P < 0.05) associated with immune and metabolic pathways. Conclusion The most optimized RNA-Seq pipeline allowed for more accurate identification of SNPs, associated positional candidate genes, and significantly associated metabolic pathways in muscle and liver tissues, providing insight on the underlying genetic architecture of feed efficiency in beef cattle.
Collapse
Affiliation(s)
- S Lam
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada
| | - J Zeidan
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada
| | - F Miglior
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada
| | - A Suárez-Vega
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada
| | - I Gómez-Redondo
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada.,Spanish National Institute for Agriculture and Food Research and Technology, Carretera de La Coruña, 28040, Madrid, Spain
| | - P A S Fonseca
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada
| | - L L Guan
- Department of Agriculture, Food & Nutritional Science, University of Alberta, Edmonton, Alberta, T6H 2P5, Canada
| | - S Waters
- Teagasc, Animal & Grassland Research and Innovation Centre, Grange, Dunsany, Co. Meath, C15 PW93, Ireland
| | - A Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road E, Guelph, Ontario, N1G2W1, Canada.
| |
Collapse
|
19
|
Carmelo VAO, Kadarmideen HN. Genetic variations (eQTLs) in muscle transcriptome and mitochondrial genes, and trans-eQTL molecular pathways in feed efficiency from Danish breeding pigs. PLoS One 2020; 15:e0239143. [PMID: 32941478 PMCID: PMC7498092 DOI: 10.1371/journal.pone.0239143] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Accepted: 08/31/2020] [Indexed: 01/08/2023] Open
Abstract
Feed efficiency (FE) is a key trait in pig production, as improvement in FE has positive economic and environmental impact. FE is a complex phenotype and testing animals for FE is costly. Therefore, there has been a desire to find functionally relevant single nucleotide polymorphisms (SNPs) as biomarkers, to improve our biological understanding of FE as well as accuracy of genomic prediction for FE. We have performed a cis- and trans- eQTL (expression quantitative trait loci) analysis, in a population of Danbred Durocs (N = 11) and Danbred Landrace (N = 27) using both a linear and ANOVA model based on muscle tissue RNA-seq. We analyzed a total of 1425x19179 or 2.7x107 Gene-SNP combinations in eQTL detection models for FE. The 1425 genes were from RNA-Seq based differential gene expression analyses using 25880 genes related to FE and additionally combined with mitochondrial genes. The 19179 SNPs were from applying stringent quality control and linkage disequilibrium filtering on genotype data using a GGP Porcine HD 70k SNP array. We applied 1000 fold bootstrapping and enrichment analysis to further validate and analyze our detected eQTLs. We identified 13 eQTLs with FDR < 0.1, affecting several genes found in previous studies of commercial pig breeds. Examples include MYO19, CPT1B, ACSL1, IER5L, CPT1A, SUCLA2, CSRNP1, PARK7 and MFF. The bootstrapping results showed statistically significant enrichment (p-value<2.2x10-16) of eQTLs with p-value < 0.01 in both cis and trans-eQTLs. Enrichment analysis of top trans-eQTLs revealed high enrichment for gene categories and gene ontologies associated with genomic context and expression regulation. This included transcription factors (p-value = 1.0x10-13), DNA-binding (GO:0003677, p-value = 8.9x10-14), DNA-binding transcription factor activity (GO:0003700,) nucleus gene (GO:0005634, p-value<2.2x10-16), negative regulation of expression (GO:0010629, p-value<2.2x10-16). These results would be useful for future genome assisted breeding of pigs to improve FE, and in the improved understanding of the functional mechanism of trans eQTLs.
Collapse
Affiliation(s)
- Victor A. O. Carmelo
- Quantitative Genomics, Bioinformatics and Computational Biology Group, Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Haja N. Kadarmideen
- Quantitative Genomics, Bioinformatics and Computational Biology Group, Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
- * E-mail:
| |
Collapse
|
20
|
González-Prendes R, Mármol-Sánchez E, Quintanilla R, Castelló A, Zidi A, Ramayo-Caldas Y, Cardoso TF, Manunza A, Cánovas Á, Amills M. About the existence of common determinants of gene expression in the porcine liver and skeletal muscle. BMC Genomics 2019; 20:518. [PMID: 31234802 PMCID: PMC6591854 DOI: 10.1186/s12864-019-5889-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2018] [Accepted: 06/07/2019] [Indexed: 12/16/2022] Open
Abstract
BACKGROUND The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population. RESULTS The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals. CONCLUSION Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.
Collapse
Affiliation(s)
- Rayner González-Prendes
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.,Departament de Producció Animal-Agrotecnio Center, Universitat de Lleida, 191 Rovira Roure, 25198, Lleida, Spain
| | - Emilio Mármol-Sánchez
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain
| | - Raquel Quintanilla
- Animal Breeding and Genetics Program, Institute for Research and Technology in Food and Agriculture (IRTA), Torre Marimon, 08140, Caldes de Montbui, Spain
| | - Anna Castelló
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Ali Zidi
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain
| | - Yuliaxis Ramayo-Caldas
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain
| | - Tainã Figueiredo Cardoso
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.,CAPES Foundation, Ministry of Education of Brazil, Brasilia D. F, Zip Code 70.040-020, Brazil
| | - Arianna Manunza
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain
| | - Ángela Cánovas
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road East, Guelph, Ontario, N1G 2W1, Canada
| | - Marcel Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain. .,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.
| |
Collapse
|
21
|
Mármol-Sánchez E, Quintanilla R, Cardoso TF, Jordana Vidal J, Amills M. Polymorphisms of the cryptochrome 2 and mitoguardin 2 genes are associated with the variation of lipid-related traits in Duroc pigs. Sci Rep 2019; 9:9025. [PMID: 31227735 PMCID: PMC6588565 DOI: 10.1038/s41598-019-45108-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2018] [Accepted: 05/24/2019] [Indexed: 01/06/2023] Open
Abstract
The genetic factors determining the phenotypic variation of porcine fatness phenotypes are still largely unknown. We investigated whether the polymorphism of eight genes (MIGA2, CRY2, NPAS2, CIART, ARNTL2, PER1, PER2 and PCK1), which display differential expression in the skeletal muscle of fasted and fed sows, is associated with the variation of lipid and mRNA expression phenotypes in Duroc pigs. The performance of an association analysis with the GEMMA software demonstrated that the rs330779504 SNP in the MIGA2 gene is associated with LDL concentration at 190 days (LDL2, corrected P-value = 0.057). Moreover, the rs320439526 SNP of the CRY2 gene displayed a significant association with stearic acid content in the longissimus dorsi muscle (LD C18:0, corrected P-value = 0.015). Both SNPs were also associated with the mRNA levels of the corresponding genes in the gluteus medius skeletal muscle. From a biological perspective these results are meaningful because MIGA2 protein plays an essential role in mitochondrial fusion, a process tightly connected with the energy status of the cell, while CRY2 is a fundamental component of the circadian clock. However, inclusion of these two SNPs in chromosome-wide association analyses demonstrated that they are not located at the peaks of significance for the two traits under study (LDL2 for rs330779504 and LD C18:0 for rs320439526), thus implying that these two SNPs do not have causal effects.
Collapse
Affiliation(s)
- Emilio Mármol-Sánchez
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Raquel Quintanilla
- Animal Breeding and Genetics Programme, Institute for Research and Technology in Food and Agriculture (IRTA), Torre Marimon, Caldes de Montbui, Spain
| | - Taina F Cardoso
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, Spain.,CAPES Foundation, Ministry of Education of Brazil, Brasilia, D. F., Brazil
| | - Jordi Jordana Vidal
- Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Marcel Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, Spain. .,Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, Bellaterra, Spain.
| |
Collapse
|
22
|
González-Prendes R, Quintanilla R, Mármol-Sánchez E, Pena RN, Ballester M, Cardoso TF, Manunza A, Casellas J, Cánovas Á, Díaz I, Noguera JL, Castelló A, Mercadé A, Amills M. Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles. BMC Genomics 2019; 20:170. [PMID: 30832586 PMCID: PMC6399881 DOI: 10.1186/s12864-019-5557-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2018] [Accepted: 02/22/2019] [Indexed: 12/23/2022] Open
Abstract
BACKGROUND Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles. RESULTS A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq. CONCLUSIONS We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.
Collapse
Affiliation(s)
- Rayner González-Prendes
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Raquel Quintanilla
- Animal Breeding and Genetics Program, Institute for Research and Technology in Food and Agriculture (IRTA), Rovira Roure 191, 25198 Lleida, Spain
| | - Emilio Mármol-Sánchez
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Ramona N. Pena
- Departament de Ciència Animal, Universitat de Lleida-Agrotecnio Centre, 25198 Lleida, Spain
| | - Maria Ballester
- Animal Breeding and Genetics Program, Institute for Research and Technology in Food and Agriculture (IRTA), Rovira Roure 191, 25198 Lleida, Spain
| | - Tainã Figueiredo Cardoso
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- CAPES Foundation, Ministry of Education of Brazil, Brasilia, DF 70.040-020 Brazil
| | - Arianna Manunza
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Joaquim Casellas
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Ángela Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, 50 Stone Road East, Guelph, Ontario N1G 2W1 Canada
| | - Isabel Díaz
- Institute for Research and Technology in Food and Agriculture (IRTA), Tecnologia dels Aliments, 17121 Monells, Spain
| | - José Luis Noguera
- Animal Breeding and Genetics Program, Institute for Research and Technology in Food and Agriculture (IRTA), Rovira Roure 191, 25198 Lleida, Spain
| | - Anna Castelló
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Anna Mercadé
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| | - Marcel Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
| |
Collapse
|
23
|
Velez-Irizarry D, Casiro S, Daza KR, Bates RO, Raney NE, Steibel JP, Ernst CW. Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs. BMC Genomics 2019; 20:3. [PMID: 30606113 PMCID: PMC6319002 DOI: 10.1186/s12864-018-5386-2] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 12/18/2018] [Indexed: 12/21/2022] Open
Abstract
Background Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs. Results Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat. Conclusion Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes. Electronic supplementary material The online version of this article (10.1186/s12864-018-5386-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
| | - Sebastian Casiro
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Kaitlyn R Daza
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Ronald O Bates
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Nancy E Raney
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Juan P Steibel
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA.,Department of Fisheries and Wildlife, Michigan State University, East Lansing, MI, 48824, USA
| | - Catherine W Ernst
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA.
| |
Collapse
|
24
|
Cardoso TF, Quintanilla R, Castelló A, González-Prendes R, Amills M, Cánovas Á. Differential expression of mRNA isoforms in the skeletal muscle of pigs with distinct growth and fatness profiles. BMC Genomics 2018; 19:145. [PMID: 29444639 PMCID: PMC5813380 DOI: 10.1186/s12864-018-4515-2] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2016] [Accepted: 01/31/2018] [Indexed: 01/03/2023] Open
Abstract
Background The identification of genes differentially expressed in the skeletal muscle of pigs displaying distinct growth and fatness profiles might contribute to identify the genetic factors that influence the phenotypic variation of such traits. So far, the majority of porcine transcriptomic studies have investigated differences in gene expression at a global scale rather than at the mRNA isoform level. In the current work, we have investigated the differential expression of mRNA isoforms in the gluteus medius (GM) muscle of 52 Duroc HIGH (increased backfat thickness, intramuscular fat and saturated and monounsaturated fatty acids contents) and LOW pigs (opposite phenotype, with an increased polyunsaturated fatty acids content). Results Our analysis revealed that 10.9% of genes expressed in the GM muscle generate alternative mRNA isoforms, with an average of 2.9 transcripts per gene. By using two different pipelines, one based on the CLC Genomics Workbench and another one on the STAR, RSEM and DESeq2 softwares, we have identified 10 mRNA isoforms that both pipelines categorize as differentially expressed in HIGH vs LOW pigs (P-value < 0.01 and ±0.6 log2fold-change). Only five mRNA isoforms, produced by the ITGA5, SEMA4D, LITAF, TIMP1 and ANXA2 genes, remain significant after correction for multiple testing (q-value < 0.05 and ±0.6 log2fold-change), being upregulated in HIGH pigs. Conclusions The increased levels of specific ITGA5, LITAF, TIMP1 and ANXA2 mRNA isoforms in HIGH pigs is consistent with reports indicating that the overexpression of these four genes is associated with obesity and metabolic disorders in humans. A broader knowledge about the functional attributes of these mRNA variants would be fundamental to elucidate the consequences of transcript diversity on the determinism of porcine phenotypes of economic interest. Electronic supplementary material The online version of this article (10.1186/s12864-018-4515-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Tainã Figueiredo Cardoso
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.,CAPES Foundation, Ministry of Education of Brazil, Brasilia D.F, 70.040-020, Brazil
| | - Raquel Quintanilla
- Animal Breeding and Genetics Programme, Institute for Research and Technology in Food and Agriculture (IRTA), Torre Marimon, 08140, Caldes de Montbui, Spain
| | - Anna Castelló
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.,Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain
| | - Rayner González-Prendes
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain
| | - Marcel Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus de la Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain. .,Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, Bellaterra, 08193, Barcelona, Spain.
| | - Ángela Cánovas
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada.
| |
Collapse
|
25
|
González-Prendes R, Quintanilla R, Amills M. Investigating the genetic regulation of the expression of 63 lipid metabolism genes in the pig skeletal muscle. Anim Genet 2017; 48:606-610. [DOI: 10.1111/age.12586] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/12/2017] [Indexed: 01/22/2023]
Affiliation(s)
- R. González-Prendes
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Campus de la Universitat Autònoma de Barcelona Bellaterra 08193 Spain
| | - R. Quintanilla
- Animal Breeding and Genetics Program; Institut de Recerca i Tecnologia Agroalimentàries (IRTA); Torre Marimon Caldes de Montbui 08140 Spain
| | - M. Amills
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Campus de la Universitat Autònoma de Barcelona Bellaterra 08193 Spain
- Departament de Ciència Animal i dels Aliments; Facultat de Veterinària; Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
| |
Collapse
|
26
|
Supakankul P, Kumchoo T, Mekchay S. Identification and characterization of novel single nucleotide polymorphism markers for fat deposition in muscle tissue of pigs using amplified fragment length polymorphism. ASIAN-AUSTRALASIAN JOURNAL OF ANIMAL SCIENCES 2016; 30:338-346. [PMID: 27608636 PMCID: PMC5337912 DOI: 10.5713/ajas.16.0200] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/09/2016] [Revised: 06/26/2016] [Accepted: 09/08/2016] [Indexed: 02/07/2023]
Abstract
Objective This study was conducted to identify and evaluate the effective single nucleotide polymorphism (SNP) markers for fat deposition in the longissimus dorsi muscles of pigs using the amplified fragment length polymorphism (AFLP) approach. Methods Sixty-four selective primer combinations were used to identify the AFLP markers in the 20 highest- and 20 lowest-intramuscular fat (IMF) content phenotypes. Five AFLP fragments were converted into simple codominant SNP markers. These SNP markers were tested in terms of their association with IMF content and fatty acid (FA) composition traits in 620 commercially crossbred pigs. Results The SSC7 g.4937240C>G marker showed an association with IMF content (p<0.05). The SSC9 g.5496647_5496662insdel marker showed a significant association with IMF content and arachidonic levels (p<0.05). The SSC10 g.71225134G>A marker revealed an association with palmitoleic and ω9 FA levels (p<0.05), while the SSC17 g.61976696G>T marker showed a significant association with IMF content and FA levels of palmitoleic, eicosenoic, arachidonic, monounsaturated fatty acids, and ω9 FA levels. However, no significant association of SSC8 g.47338181G>A was observed with any IMF and FA levels in this study. Conclusion Four SNP markers (SSC7 g.4937240C>G, SSC9 g.5496647_5496662insdel, SSC10 g.71225134G>A, and SSC17 g.61976696G>T) were found to be associated with IMF and/or FA content traits in commercially crossbred pigs. These findings provide evidence of the novel SNP markers as being potentially useful for selecting pigs with the desirable IMF content and FA composition.
Collapse
Affiliation(s)
- Pantaporn Supakankul
- Department of Animal and Aquatic Sciences, Faculty of Agriculture, Chiang Mai University, Chiang Mai 50200, Thailand.,Center of Excellence on Agricultural Biotechnology: (AG-BIO/PERDO-CHE), Bangkok 10900, Thailand.,Human and Animal Biotechnology Program, Graduate School, Chiang Mai University, Chiang Mai 50200, Thailand
| | - Tanavadee Kumchoo
- Department of Animal and Aquatic Sciences, Faculty of Agriculture, Chiang Mai University, Chiang Mai 50200, Thailand
| | - Supamit Mekchay
- Department of Animal and Aquatic Sciences, Faculty of Agriculture, Chiang Mai University, Chiang Mai 50200, Thailand
| |
Collapse
|
27
|
Identification of genes showing differential expression profile associated with growth rate in skeletal muscle tissue of Landrace weanling pig. J Genet 2016; 95:341-7. [DOI: 10.1007/s12041-016-0643-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
28
|
Pena RN, Quintanilla R, Manunza A, Gallardo D, Casellas J, Amills M. Application of the microarray technology to the transcriptional analysis of muscle phenotypes in pigs. Anim Genet 2014; 45:311-21. [DOI: 10.1111/age.12146] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/05/2014] [Indexed: 01/09/2023]
Affiliation(s)
- R. N. Pena
- Department of Animal Production; University of Lleida-Agrotecnio Center; 25198 Lleida Spain
| | | | - A. Manunza
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Universitat Autònoma de Barcelona; 08193 Bellaterra Spain
| | - D. Gallardo
- Departament de Ciència Animal i dels Aliments; Facultat de Veterinària; Universitat Autònoma de Barcelona; 08193 Bellaterra Spain
| | - J. Casellas
- Departament de Ciència Animal i dels Aliments; Facultat de Veterinària; Universitat Autònoma de Barcelona; 08193 Bellaterra Spain
| | - M. Amills
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Universitat Autònoma de Barcelona; 08193 Bellaterra Spain
- Departament de Ciència Animal i dels Aliments; Facultat de Veterinària; Universitat Autònoma de Barcelona; 08193 Bellaterra Spain
| |
Collapse
|
29
|
Molecular advances in QTL discovery and application in pig breeding. Trends Genet 2013; 29:215-24. [DOI: 10.1016/j.tig.2013.02.002] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2012] [Revised: 02/12/2013] [Accepted: 02/13/2013] [Indexed: 11/21/2022]
|