1
|
Davoudi P, Do DN, Rathgeber B, Colombo S, Sargolzaei M, Plastow G, Wang Z, Miar Y. Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data. J Anim Breed Genet 2024; 141:507-520. [PMID: 38389405 DOI: 10.1111/jbg.12859] [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: 11/30/2023] [Revised: 01/27/2024] [Accepted: 02/02/2024] [Indexed: 02/24/2024]
Abstract
The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.
Collapse
Affiliation(s)
- Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| | - Duy Ngoc Do
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| | - Bruce Rathgeber
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| | - Stefanie Colombo
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| | - Mehdi Sargolzaei
- Department of Pathobiology, University of Guelph, Guelph, Ontario, Canada
- Select Sires Inc., Plain City, Ohio, USA
| | - Graham Plastow
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
| | - Zhiquan Wang
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
| | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| |
Collapse
|
2
|
Christodoulides N, Urgiles VL, Guayasamin JM, Savage AE. Selection and Gene Duplication Associated With High-Elevation Diversification in Pristimantis, the Largest Terrestrial Vertebrate Genus. Genome Biol Evol 2024; 16:evae167. [PMID: 39109890 PMCID: PMC11342244 DOI: 10.1093/gbe/evae167] [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] [Accepted: 07/24/2024] [Indexed: 08/24/2024] Open
Abstract
The genus Pristimantis diversified in the tropical Andes mountains and is the most speciose genus of terrestrial vertebrates. Pristimantis are notable among frogs in that they thrive at high elevations (>2,000 m) and are direct developers without a tadpole stage. Despite their ecological significance, little is known about the genetic and physiological traits enabling their success. We conducted transcriptomic analysis on seven Pristimantis species sampled across elevations in the Ecuadorean Andes to explore three hypotheses for their success: (i) unique genes are under selection relative to all other frogs, (ii) common selection occurs across all direct developers, or (iii) common selection occurs across all high-elevation frog clades. Comparative analysis with 34 frog species revealed unique positive selection in Pristimantis genes related to aerobic respiration, hemostasis, signaling, cellular transportation of proteins and ions, and immunity. Additionally, we detected positive selection across all direct developers for genes associated with oxygenase activity and metal ion binding. While many genes under selection in Pristimantis were not positively selected in other high-elevation frog species, we identified some shared genes and pathways linked to lipid metabolism, innate immunity, and cellular redox processes. We observed more positive selection in duplicated- versus single-copy genes, while relaxed purifying selection was prevalent in single-copy genes. Notably, copy number of an innate immunity complement gene was positively correlated with Pristimantis species elevation. Our findings contribute novel insights into the genetic basis of adaptation in Pristimantis and provide a foundation for future studies on the evolutionary mechanisms leading to direct development and coping with high elevations.
Collapse
Affiliation(s)
| | - Veronica L Urgiles
- Department of Biology, University of Central Florida, Orlando, FL, USA
- Departamento de herpetologia, Instituto Nacional de Biodiversidad del Ecuador, Quito, Ecuador
| | - Juan M Guayasamin
- Colegio de Ciencias Biológicas y Ambientales COCIBA, Instituto Biósfera, Laboratorio de Biología Evolutiva, Universidad San Francisco de Quito USFQ, Quito, Ecuador
- Ingeniería en Biodiversidad y Recursos Genéticos, Centro de Biodiversidad y Cambio Climático BioCamb, Universidad Tecnológica Indoamérica, Quito, Ecuador
| | - Anna E Savage
- Department of Biology, University of Central Florida, Orlando, FL, USA
| |
Collapse
|
3
|
Davoudi P, Do DN, Rathgeber B, Colombo S, Sargolzaei M, Plastow G, Wang Z, Miar Y. Identification of consensus homozygous regions and their associations with growth and feed efficiency traits in American mink. BMC Genom Data 2024; 25:68. [PMID: 38982354 PMCID: PMC11234557 DOI: 10.1186/s12863-024-01252-8] [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: 03/15/2024] [Accepted: 07/01/2024] [Indexed: 07/11/2024] Open
Abstract
The recent chromosome-based genome assembly and the newly developed 70K single nucleotide polymorphism (SNP) array for American mink (Neogale vison) facilitate the identification of genetic variants underlying complex traits in this species. The objective of this study was to evaluate the association between consensus runs of homozygosity (ROH) with growth and feed efficiency traits in American mink. A subsample of two mink populations (n = 2,986) were genotyped using the Affymetrix Mink 70K SNP array. The identified ROH segments were included simultaneously, concatenated into consensus regions, and the ROH-based association studies were carried out with linear mixed models considering a genomic relationship matrix for 11 growth and feed efficiency traits implemented in ASReml-R version 4. In total, 298,313 ROH were identified across all individuals, with an average length and coverage of 4.16 Mb and 414.8 Mb, respectively. After merging ROH segments, 196 consensus ROH regions were detected and used for genome-wide ROH-based association analysis. Thirteen consensus ROH regions were significantly (P < 0.01) associated with growth and feed efficiency traits. Several candidate genes within the significant regions are known for their involvement in growth and body size development, including MEF2A, ADAMTS17, POU3F2, and TYRO3. In addition, we found ten consensus ROH regions, defined as ROH islands, with frequencies over 80% of the population. These islands harbored 12 annotated genes, some of which were related to immune system processes such as DTX3L, PARP9, PARP14, CD86, and HCLS1. This is the first study to explore the associations between homozygous regions with growth and feed efficiency traits in American mink. Our findings shed the light on the effects of homozygosity in the mink genome on growth and feed efficiency traits, that can be utilized in developing a sustainable breeding program for mink.
Collapse
Affiliation(s)
- Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Duy Ngoc Do
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Bruce Rathgeber
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Stefanie Colombo
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Mehdi Sargolzaei
- Department of Pathobiology, University of Guelph, Guelph, ON, Canada
- Select Sires Inc, Plain City, OH, USA
| | - Graham Plastow
- Department of Agricultural, Food and Nutritional Science, Livestock Gentec, University of Alberta, Edmonton, AB, Canada
| | - Zhiquan Wang
- Department of Agricultural, Food and Nutritional Science, Livestock Gentec, University of Alberta, Edmonton, AB, Canada
| | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada.
| |
Collapse
|
4
|
Han ZP, Yang RZ, Zhou W, Zhang LL, Wang JR, Liu CJ, Liu SD. Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert. BMC Genomics 2024; 25:681. [PMID: 38982349 PMCID: PMC11232224 DOI: 10.1186/s12864-024-10581-y] [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/29/2023] [Accepted: 06/28/2024] [Indexed: 07/11/2024] Open
Abstract
Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.
Collapse
Affiliation(s)
- Zhi-Peng Han
- College of Animal Science and Technology, Tarim University, Alar, 843300, China
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China
| | - Rui-Zhi Yang
- College of Animal Science and Technology, Tarim University, Alar, 843300, China
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China
| | - Wen Zhou
- College of Animal Science and Technology, Tarim University, Alar, 843300, China
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China
| | - Lu-Lu Zhang
- College of Animal Science and Technology, Tarim University, Alar, 843300, China
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China
| | - Jie-Ru Wang
- College of Animal Science and Technology, Tarim University, Alar, 843300, China
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China
| | - Chun-Jie Liu
- College of Animal Science and Technology, Tarim University, Alar, 843300, China
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China
| | - Shu-Dong Liu
- College of Animal Science and Technology, Tarim University, Alar, 843300, China.
- Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Alar, 843300, China.
| |
Collapse
|
5
|
Vostry L, Vostra-Vydrova H, Moravcikova N, Kasarda R, Margetin M, Rychtarova J, Drzaic I, Shihabi M, Cubric-Curik V, Sölkner J, Curik I. Genomic analysis of conservation status, population structure and admixture in local Czech and Slovak dairy goat breeds. J Dairy Sci 2024:S0022-0302(24)00937-8. [PMID: 38908686 DOI: 10.3168/jds.2023-24607] [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/2023] [Accepted: 05/30/2024] [Indexed: 06/24/2024]
Abstract
While dairy goat production, characterized by traditional production on small farms, is an important source of income in the Czech Republic and Slovakia, locally adapted breeds have not been fully consolidated over the last 100 years due to large fluctuations in population size and inconsistent breeding programs that allowed for different crossbreeding strategies. Our main objective in this study was therefore to assess the conservation status of 4 Czech (Alpine Goat, White Shorthair, Brown Shorthair and Czech Landrace) and one Slovak (Slovak White Shorthair) local goat breeds, to analyze their population structure and admixture, and to estimate their relatedness to several neighboring breeds. Our analyses included 142 goats belonging to 5 local breeds genotyped with the Illumina 50K BeadChip and 618 previously genotyped animals representing 15 goat breeds from Austria and Switzerland (all analyses based on 46,862 autosomal SNPs and 760 animals). In general, the conservation status of the Czech and Slovak local goat breeds was satisfactory, with the exception of the Brown Shorthair goat, as the analyzed parameters (heterozygosity, haplotype richness, ROH-based inbreeding and effective population size) were mostly above the median of 20 breeds. However, for all 5 Czech and Slovakian breeds, an examination of historical effective population size indicated a substantial decline about 8 to 22 generations ago. In addition, our study revealed that the Czech and Slovakian breeds are not fully consolidated; for instance, White Shorthair and Brown Shorthair were not clearly distinguishable. Considerable admixture, especially in Czech Landrace (effective number of parental clusters equal to 4.2), and low but numerous migration rates from other Austrian and Swiss breeds were found. These results provide valuable insights for future breeding programs and genetic diversity management of local Czech and Slovak goat breeds.
Collapse
Affiliation(s)
- Lubos Vostry
- Czech University of Life Science Prague, Kamycka 129, 16500 Prague, Czech Republic.
| | - Hana Vostra-Vydrova
- Czech University of Life Science Prague, Kamycka 129, 16500 Prague, Czech Republic
| | - Nina Moravcikova
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Radovan Kasarda
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Milan Margetin
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Jana Rychtarova
- Institute of Animal Science, Přátelství 815, 104 00 Prague, Czech Republic
| | - Ivana Drzaic
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Mario Shihabi
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Vlatka Cubric-Curik
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Johan Sölkner
- University of Natural Resources & Life Sciences Vienna, Gregor-Mendel-Strasse 33, 1180 Vienna, Austria
| | - Ino Curik
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia; Institute of Animal Sciences, Hungarian University of Agriculture and Life Sciences (MATE), Guba Sándor u. 40, 7400 Kaposvár, Hungary.
| |
Collapse
|
6
|
Zhao Q, Huang C, Chen Q, Su Y, Zhang Y, Wang R, Su R, Xu H, Liu S, Ma Y, Zhao Q, Ye S. Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat. Animals (Basel) 2024; 14:1246. [PMID: 38672394 PMCID: PMC11047310 DOI: 10.3390/ani14081246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Revised: 04/15/2024] [Accepted: 04/18/2024] [Indexed: 04/28/2024] Open
Abstract
Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.
Collapse
Affiliation(s)
- Qian Zhao
- Department of Animal Breeding and Reproduction, College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (Q.Z.); (C.H.)
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (Q.C.); (Y.S.); (Y.M.)
| | - Chang Huang
- Department of Animal Breeding and Reproduction, College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (Q.Z.); (C.H.)
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (Q.C.); (Y.S.); (Y.M.)
| | - Qian Chen
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (Q.C.); (Y.S.); (Y.M.)
| | - Yingxiao Su
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (Q.C.); (Y.S.); (Y.M.)
| | - Yanjun Zhang
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot 010018, China; (Y.Z.); (R.W.); (R.S.)
| | - Ruijun Wang
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot 010018, China; (Y.Z.); (R.W.); (R.S.)
| | - Rui Su
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot 010018, China; (Y.Z.); (R.W.); (R.S.)
| | - Huijuan Xu
- Chifeng Hanshan White Cashmere Goat Breeding Farm, Chifeng 024506, China; (H.X.); (S.L.)
| | - Shucai Liu
- Chifeng Hanshan White Cashmere Goat Breeding Farm, Chifeng 024506, China; (H.X.); (S.L.)
| | - Yuehui Ma
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (Q.C.); (Y.S.); (Y.M.)
| | - Qianjun Zhao
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (Q.C.); (Y.S.); (Y.M.)
| | - Shaohui Ye
- Department of Animal Breeding and Reproduction, College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (Q.Z.); (C.H.)
| |
Collapse
|
7
|
Manunza A, Ramírez-Díaz J, Rincón Flórez JC, Almeida de Oliveira T. Editorial: Environmental and genomic strategies for conservation and selection in small ruminants. Front Vet Sci 2024; 11:1396289. [PMID: 38628938 PMCID: PMC11018938 DOI: 10.3389/fvets.2024.1396289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 03/12/2024] [Indexed: 04/19/2024] Open
Affiliation(s)
- Arianna Manunza
- Institute of Agricultural Biology and Biotechnology, National Research Council (CNR), Milan, Italy
| | - Johanna Ramírez-Díaz
- Institute of Agricultural Biology and Biotechnology, National Research Council (CNR), Milan, Italy
| | | | | |
Collapse
|
8
|
Pérez S, Calvo JH, Calvete C, Joy M, Lobón S. Mitigation and animal response to water stress in small ruminants. Anim Front 2023; 13:81-88. [PMID: 37841759 PMCID: PMC10575300 DOI: 10.1093/af/vfad049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2023] Open
Affiliation(s)
- Sara Pérez
- Departamento de Ciencia Animal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), 50059 Zaragoza, España
| | - Jorge Hugo Calvo
- Departamento de Ciencia Animal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), 50059 Zaragoza, España
- Instituto Agroalimentario de Aragón – IA2 (CITA-Universidad deZaragoza)
- Fundación Agencia Aragonesa para la Investigación y el Desarrollo (ARAID), 50018 Zaragoza, España
| | - Carlos Calvete
- Departamento de Ciencia Animal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), 50059 Zaragoza, España
- Instituto Agroalimentario de Aragón – IA2 (CITA-Universidad deZaragoza)
| | - Margalida Joy
- Departamento de Ciencia Animal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), 50059 Zaragoza, España
- Instituto Agroalimentario de Aragón – IA2 (CITA-Universidad deZaragoza)
| | - Sandra Lobón
- Departamento de Ciencia Animal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), 50059 Zaragoza, España
- Instituto Agroalimentario de Aragón – IA2 (CITA-Universidad deZaragoza)
| |
Collapse
|