1
|
Gozdek M, Mucha S, Prostek A, Sadkowski T. Selected Monogenic Genetic Diseases in Holstein Cattle-A Review. Genes (Basel) 2024; 15:1052. [PMID: 39202412 PMCID: PMC11353376 DOI: 10.3390/genes15081052] [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: 07/18/2024] [Revised: 08/03/2024] [Accepted: 08/05/2024] [Indexed: 09/03/2024] Open
Abstract
Genetic disorders arise from alterations in the hereditary information encoded in DNA, leading to potential detrimental effects on the well-being and vitality of organisms. Within the bovine population, genetic conditions inherited in an autosomal recessive manner are frequently associated with particular breeds. In recent years, several recessive haplotypes and a few causative mutations have been discovered in Holstein cattle: CDH (Holstein cholesterol deficiency), haplotypes with a homozygous deficiency in Holstein (HH1, HH3, HH4, HH5, HH6 and HH7), BLAD (bovine leukocyte adhesion deficiency) and DUMPS (deficiency of uridine monophosphate synthase). All of these diseases are inherited in an autosomal recessive manner. From a breeding perspective, recessive mutations specifically exhibit considerable detrimental effects and are a significant problem for breeders, exposing them to economic losses. Individual mutations can cause embryo death at any stage of pregnancy. Only genetic research and conscious selection of animals for mating will lead to a reduction in the number of carriers and elimination of mutations from the population.
Collapse
Affiliation(s)
- Marta Gozdek
- Department of Physiological Sciences, Institute of Veterinary Medicine, Warsaw University of Life Sciences, 02-776 Warsaw, Poland; (M.G.); (A.P.)
- Polish Federation of Cattle Breeders and Dairy Farmers, 00-515 Warsaw, Poland;
| | - Sebastian Mucha
- Polish Federation of Cattle Breeders and Dairy Farmers, 00-515 Warsaw, Poland;
| | - Adam Prostek
- Department of Physiological Sciences, Institute of Veterinary Medicine, Warsaw University of Life Sciences, 02-776 Warsaw, Poland; (M.G.); (A.P.)
| | - Tomasz Sadkowski
- Department of Physiological Sciences, Institute of Veterinary Medicine, Warsaw University of Life Sciences, 02-776 Warsaw, Poland; (M.G.); (A.P.)
| |
Collapse
|
2
|
Panigrahi M, Rajawat D, Nayak SS, Jain K, Vaidhya A, Prakash R, Sharma A, Parida S, Bhushan B, Dutt T. Genomic insights into key genes and QTLs involved in cattle reproduction. Gene 2024; 917:148465. [PMID: 38621496 DOI: 10.1016/j.gene.2024.148465] [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: 02/25/2024] [Revised: 04/05/2024] [Accepted: 04/10/2024] [Indexed: 04/17/2024]
Abstract
From an economic standpoint, reproductive characteristics are fundamental for sustainable production, particularly for monotocous livestock like cattle. A longer inter-calving interval is indicative of low reproductive capacity. This issue changes the dynamics of current and future lactations since it necessitates more inseminations, veterinary care, and hormone interventions. Various reproductive phenotypes, including ovulation, mating, fertility, pregnancy, embryonic growth, and calving-related traits, are observed in dairy cattle, and these traits have been associated with several QTLs. Calving ease, age at puberty, scrotal circumference, and inseminations per conception have been associated with 4437, 10623, 10498, and 2476 Quantitative Trait Loci (QTLs), respectively. This data offers valuable insights into enhancing and comprehending reproductive traits in livestock breeding. Studying QTLs associated with reproductive traits has far-reaching implications across various fields, from agriculture and animal husbandry to human health, evolutionary biology, and conservation. It provides the foundation for informed breeding practices, advances in biotechnology, and a deeper understanding of the genetic underpinnings of reproduction.
Collapse
Affiliation(s)
- Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India.
| | - Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Sonali Sonejita Nayak
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Karan Jain
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Ayushi Vaidhya
- Division of Pharmacology & Toxicology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Ravi Prakash
- Division of Pharmacology & Toxicology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Anurodh Sharma
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Subhashree Parida
- Division of Pharmacology & Toxicology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Triveni Dutt
- Livestock Production and Management Section, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| |
Collapse
|
3
|
Saini T, Chauhan A, Ahmad SF, Kumar A, Vaishnav S, Singh S, Mehrotra A, Bhushan B, Gaur GK, Dutt T. Elucidation of population stratifying markers and selective sweeps in crossbred Landlly pig population using genome-wide SNP data. Mamm Genome 2024; 35:170-185. [PMID: 38485788 DOI: 10.1007/s00335-024-10029-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Accepted: 01/23/2024] [Indexed: 05/29/2024]
Abstract
The present study was aimed at the identification of population stratifying markers from the commercial porcine SNP 60K array and elucidate the genome-wide selective sweeps in the crossbred Landlly pig population. Original genotyping data, generated on Landlly pigs, was merged in various combinations with global suid breeds that were grouped as exotic (global pig breeds excluding Indian and Chinese), Chinese (Chinese pig breeds only), and outgroup pig populations. Post quality control, the genome-wide SNPs were ranked for their stratifying power within each dataset in TRES (using three different criteria) and FIFS programs and top-ranked SNPs (0.5K, 1K, 2K, 3K, and 4K densities) were selected. PCA plots were used to assess the stratification power of low-density panels. Selective sweeps were elucidated in the Landlly population using intra- and inter-population haplotype statistics. Additionally, Tajima's D-statistics were calculated to determine the status of balancing selection in the Landlly population. PCA plots showed 0.5K marker density to effectively stratify Landlly from other pig populations. The A-score in DAPC program revealed the Delta statistic of marker selection to outperform other methods (informativeness and FST methods) and that 3000-marker density was suitable for stratification of Landlly animals from exotic pig populations. The results from selective sweep analysis revealed the Landlly population to be under selection for mammary (NAV2), reproductive efficiency (JMY, SERGEF, and MAP3K20), body conformation (FHIT, WNT2, ASRB, DMGDH, and BHMT), feed efficiency (CSRNP1 and ADRA1A), and immunity (U6, MYO3B, RBMS3, and FAM78B) traits. More than two methods suggested sweeps for immunity and feed efficiency traits, thus giving a strong indication for selection in this direction. The study is the first of its kind in Indian pig breeds with a comparison against global breeds. In conclusion, 500 markers were able to effectively stratify the breeds. Different traits under selective sweeps (natural or artificial selection) can be exploited for further improvement.
Collapse
Affiliation(s)
- Tapendra Saini
- Animal Genetics Division, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
| | - Anuj Chauhan
- Swine Production Farm, LPM Section, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India.
| | - Sheikh Firdous Ahmad
- Animal Genetics Division, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
| | - Amit Kumar
- Animal Genetics Division, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
| | - Sakshi Vaishnav
- Animal Genetics Division, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
| | - Shivani Singh
- Swine Production Farm, LPM Section, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
| | | | - Bharat Bhushan
- Animal Genetics Division, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
| | - G K Gaur
- Swine Production Farm, LPM Section, ICAR-Indian Veterinary Research Institute, Izatnagar, 243122, India
- ADG Animal Production & Breeding, ICAR, New Delhi, 110001, India
| | - Triveni Dutt
- Indian Veterinary Research Institute, Izatnagar, 243122, India
| |
Collapse
|
4
|
Dias MS, Pedrosa VB, Rocha da Cruz VA, Silva MR, Batista Pinto LF. Genome-wide association and functional annotation analysis for the calving interval in Nellore cattle. Theriogenology 2024; 218:214-222. [PMID: 38350227 DOI: 10.1016/j.theriogenology.2024.01.034] [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: 11/05/2023] [Revised: 01/24/2024] [Accepted: 01/25/2024] [Indexed: 02/15/2024]
Abstract
Calving interval (CI) measures the number of days between two consecutive calves of the same cow, and previous studies based on phenotype and pedigree data reported low heritability for this trait. However, the genetic architecture of CI in the Nellore breed was not evaluated based on genomic data. Thus, this study aimed to estimate the heritability based on genomic data and carry out a genome-wide association study (GWAS) for CI in the Nellore breed, using 12,599 pedigree records, 5078 CI records, and 3818 animals genotyped with 50k SNPchip panel. Both quality control and GWAS were performed in BLUPF90 family packages, which use the single-step genomic best linear unbiased predictor (ssGBLUP) method. The average CI was 427.6 days, with a standard deviation of 106.9 and a total range of 270-730 days. The heritability estimate was 0.04 ± 0.04. The p-values of GWAS analysis resulted in a genomic inflation factor (lambda) of 1.08. The only significant SNP (rs136725686) at the genome-wide level (p-value = 1.53E-06) was located on BTA13. Other 19 SNPs were significant at the chromosome-wide level, distributed on BTA1, 2, 3, 6, 10, 13, 14, 17, 18, 22, and 26. Functional annotation analysis found thirty-six protein-coding genes, including genes related to cell cycle (RAD21, BCAR3), oocyte function (LHX8, CLPX, UTP23), immune system (TXK, TEC, NFATC2), endocrine function (LRRFIP2, GPR158), estrous cycle (SLC38A7), and female fertility (CCK, LYZL4, TRAK1, FOXP1, STAC). Therefore, CI is a complex trait with small heritability in Nellore cattle, and various biological processes may be involved with the genetic architecture of CI in Nellore cattle.
Collapse
Affiliation(s)
- Mayra Silva Dias
- Federal University of Bahia, Animal Science Department, Av. Milton Santos, 500, Ondina, Salvador, BA, 40170-110, Brazil.
| | | | | | - Marcio Ribeiro Silva
- Melhore Animal and Katayama Agropecuaria Lda, Guararapes, SP, 16700-000, Brazil.
| | - Luis Fernando Batista Pinto
- Federal University of Bahia, Animal Science Department, Av. Milton Santos, 500, Ondina, Salvador, BA, 40170-110, Brazil.
| |
Collapse
|
5
|
Turner KA, Achinger L, Kong D, Kluczynski DF, Fishman EL, Phillips A, Saltzman B, Loncarek J, Harstine BR, Avidor-Reiss T. Abnormal centriolar biomarker ratios correlate with unexplained bull artificial insemination subfertility: a pilot study. Sci Rep 2023; 13:18338. [PMID: 37884598 PMCID: PMC10603076 DOI: 10.1038/s41598-023-45162-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: 07/25/2023] [Accepted: 10/17/2023] [Indexed: 10/28/2023] Open
Abstract
The mechanisms underlying male infertility are poorly understood. Most mammalian spermatozoa have two centrioles: the typical barrel-shaped proximal centriole (PC) and the atypical fan-like distal centriole (DC) connected to the axoneme (Ax). These structures are essential for fertility. However, the relationship between centriole quality and subfertility (reduced fertility) is not well established. Here, we tested the hypothesis that assessing sperm centriole quality can identify cattle subfertility. By comparing sperm from 25 fertile and 6 subfertile bulls, all with normal semen analyses, we found that unexplained subfertility and lower sire conception rates (pregnancy rate from artificial insemination in cattle) correlate with abnormal centriolar biomarker distribution. Fluorescence-based Ratiometric Analysis of Sperm Centrioles (FRAC) found only four fertile bulls (4/25, 16%) had positive FRAC tests (having one or more mean FRAC ratios outside of the distribution range in a group's high-quality sperm population), whereas all of the subfertile bulls (6/6, 100%) had positive FRAC tests (P = 0.00008). The most sensitive biomarker was acetylated tubulin, which had a novel labeling pattern between the DC and Ax. These data suggest that FRAC and acetylated tubulin labeling can identify bull subfertility that remains undetected by current methods and may provide insight into a novel mechanism of subfertility.
Collapse
Affiliation(s)
- Katerina A Turner
- Department of Biological Sciences, College of Natural Sciences and Mathematics, University of Toledo, 3050 W. Towerview Blvd, Toledo, OH, 43606, USA
| | - Luke Achinger
- Department of Biological Sciences, College of Natural Sciences and Mathematics, University of Toledo, 3050 W. Towerview Blvd, Toledo, OH, 43606, USA
| | - Dong Kong
- Laboratory of Protein Dynamics and Signaling, Center for Cancer Research, National Institutes of Health, National Cancer Institute, Frederick, MD, USA
| | - Derek F Kluczynski
- Department of Biological Sciences, College of Natural Sciences and Mathematics, University of Toledo, 3050 W. Towerview Blvd, Toledo, OH, 43606, USA
| | - Emily Lillian Fishman
- Department of Biological Sciences, College of Natural Sciences and Mathematics, University of Toledo, 3050 W. Towerview Blvd, Toledo, OH, 43606, USA
| | - Audrey Phillips
- Department of Biological Sciences, College of Natural Sciences and Mathematics, University of Toledo, 3050 W. Towerview Blvd, Toledo, OH, 43606, USA
| | - Barbara Saltzman
- Department of Population Health, College of Health and Human Services, University of Toledo, Toledo, OH, USA
| | - Jadranka Loncarek
- Laboratory of Protein Dynamics and Signaling, Center for Cancer Research, National Institutes of Health, National Cancer Institute, Frederick, MD, USA
| | | | - Tomer Avidor-Reiss
- Department of Biological Sciences, College of Natural Sciences and Mathematics, University of Toledo, 3050 W. Towerview Blvd, Toledo, OH, 43606, USA.
- Department of Urology, College of Medicine and Life Sciences, University of Toledo, Toledo, OH, USA.
| |
Collapse
|
6
|
Marrella MA, Biase FH. A multi-omics analysis identifies molecular features associated with fertility in heifers (Bos taurus). Sci Rep 2023; 13:12664. [PMID: 37542054 PMCID: PMC10403585 DOI: 10.1038/s41598-023-39858-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Accepted: 08/01/2023] [Indexed: 08/06/2023] Open
Abstract
Infertility or subfertility is a critical barrier to sustainable cattle production, including in heifers. The development of heifers that do not produce a calf within an optimum window of time is a critical factor for the profitability and sustainability of the cattle industry. In parallel, heifers are an excellent biomedical model for understanding the underlying etiology of infertility because well-nourished heifers can still be infertile, mostly because of inherent physiological and genetic causes. Using a high-density single nucleotide polymorphism (SNP) chip, we collected genotypic data, which were analyzed using an association analysis in PLINK with Fisher's exact test. We also produced quantitative transcriptome data and proteome data. Transcriptome data were analyzed using the quasi-likelihood test followed by the Wald's test, and the likelihood test and proteome data were analyzed using a generalized mixed model and Student's t-test. We identified two SNPs significantly associated with heifer fertility (rs110918927, chr12: 85648422, P = 6.7 × 10-7; and rs109366560, chr11:37666527, P = 2.6 × 10-5). We identified two genes with differential transcript abundance (eFDR ≤ 0.002) between the two groups (Fertile and Sub-Fertile): Adipocyte Plasma Membrane Associated Protein (APMAP, 1.16 greater abundance in the Fertile group) and Dynein Axonemal Intermediate Chain 7 (DNAI7, 1.23 greater abundance in the Sub-Fertile group). Our analysis revealed that the protein Alpha-ketoglutarate-dependent dioxygenase FTO was more abundant in the plasma collected from Fertile heifers relative to their Sub-Fertile counterparts (FDR < 0.05). Lastly, an integrative analysis of the three datasets identified a series of molecular features (SNPs, gene transcripts, and proteins) that discriminated 21 out of 22 heifers correctly based on their fertility category. Our multi-omics analyses confirm the complex nature of female fertility. Very importantly, our results also highlight differences in the molecular profile of heifers associated with fertility that transcend the constraints of breed-specific genetic background.
Collapse
Affiliation(s)
- Mackenzie A Marrella
- School of Animal Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA
| | - Fernando H Biase
- School of Animal Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA.
| |
Collapse
|
7
|
Sölzer N, May K, Yin T, König S. Genomic analyses of claw disorders in Holstein cows: Genetic parameters, trait associations, and genome-wide associations considering interactions of SNP and heat stress. J Dairy Sci 2022; 105:8218-8236. [PMID: 36028345 DOI: 10.3168/jds.2022-22087] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Accepted: 05/31/2022] [Indexed: 11/19/2022]
Abstract
The aim of the present study was an in-depth genomic analysis to understand the genomic mechanisms of the 3 claw disorders dermatitis digitalis (DD), interdigital hyperplasia (HYP), and sole ulcer (SU). In this regard, we estimated genetic parameters based on genomic relationship matrices, performed genome-wide association studies, annotated potential candidate genes, and inferred genetic associations with breeding goal traits considering the most important chromosomal segments. As a further novelty of this study, we inferred possible SNP × heat stress interactions for claw disorders. The study consisted of 17,264 first-lactation Holstein Friesian cows kept in 50 large-scale contract herds. The disease prevalence was 15.96, 2.36, and 8.20% for DD, HYP, and SU, respectively. The remaining breeding goal traits consisted of type traits of the feet and leg composite, female fertility, health traits, and 305-d production traits. The final genotype data set included 44,474 SNPs from the 17,264 genotyped cows. Heritabilities for DD, HYP, and SU were estimated in linear and threshold models considering the genomic relationship matrix (G matrix). Genetic correlations with breeding goal traits based on G were estimated in a series of bivariate linear models, which were verified via SNP effect correlations for specific chromosome segments (i.e., segments harboring potential candidate genes for DD, HYP, and SU). Genome-wide association studies were performed for all traits in a case-control design by applying a single SNP linear mixed model. Furthermore, for DD, HYP, and SU, we modeled SNP × heat stress interactions in genome-wide association studies. Single nucleotide polymorphism-based heritabilities were 0.04 and 0.08 for DD, 0.03 and 0.10 for SU, and 0.03 and 0.23 for HYP from linear and threshold models, respectively. The genetic correlations between DD, HYP, and SU with conformation traits from the feet and leg composite were positive throughout, indicating the value of indirect selection on conformation traits to improve claw health. Genetic correlations between DD, SU, and HYP with other breeding goal traits indicated impaired female fertility, impaired udder health status, and productivity decline of diseased cows. Genetic correlations among DD, SU, and HYP were moderate to large, indicating that different claw disorders have similar genetic mechanisms. Nevertheless, we identified disease-specific potential candidate genes, and genetic associations based on the surrounding SNPs partly differed from the genetic correlations. Especially for candidate genes contributing to 2 traits simultaneously, correlations based on SNP effects from the respective chromosome segment were close to 1 or to -1. In this regard, we annotated the candidate genes KRT33A and KRT33B for HYP and DD, KIF27 for HYP and calving to first insemination, and MAN1A1 for SU and the production traits. For SNP × heat stress interactions, we identified significant SNPs on BTA 2, 4, 5, 7, 8, 9, 13, 22, 25, and 28, and we annotated the potential candidate genes FSIP2, CLCN1, ADGRV1, DOP1A, THBD, and RHOBTB1. Results indicate gene-specific mechanisms of the claw disorders only in specific environments.
Collapse
Affiliation(s)
- Niklas Sölzer
- Institute of Animal Breeding and Genetics, Justus-Liebig-University Giessen, 35390 Giessen, Germany
| | - Katharina May
- Institute of Animal Breeding and Genetics, Justus-Liebig-University Giessen, 35390 Giessen, Germany
| | - Tong Yin
- Institute of Animal Breeding and Genetics, Justus-Liebig-University Giessen, 35390 Giessen, Germany
| | - Sven König
- Institute of Animal Breeding and Genetics, Justus-Liebig-University Giessen, 35390 Giessen, Germany.
| |
Collapse
|
8
|
Hay E, Toghiani S, Roberts AJ, Paim T, Kuehn LA, Blackburn HD. Genetic architecture of a composite beef cattle population. J Anim Sci 2022; 100:6623572. [PMID: 35771897 DOI: 10.1093/jas/skac230] [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: 01/27/2022] [Accepted: 06/28/2022] [Indexed: 11/15/2022] Open
Abstract
Composite breeds are widely used in the beef industry. Composites allow producers to combine desirable traits from the progenitor breeds and simplify herd management, without repeated crossbreeding and maintenance of purebreds. In this study, genomic information was used to evaluate the genetic composition and characteristics of a three-breed beef cattle composite. This composite population referred to as Composite Gene Combination (CGC) consisted of 50% Red Angus, 25% Charolais, 25% Tarentaise. A total of 248 animals were used in this study CGC (n=79), Red Angus (n=61), Charolais (n=79) and Tarentaise (n=29). All animals were genotyped with 777k HD panel. Principal component and ADMIXTURE analyses were carried out to evaluate the genetic structure of CGC animals. The ADMIXTURE revealed the proportion of Tarentaise increased to approximately 57% while Charolais decreased to approximately 5%, and Red Angus decreased to 38% across generations. To evaluate these changes in the genomic composition across different breeds and in CGC across generations runs of homozygosity (ROH) were conducted. This analysis showed Red Angus to have the highest total length of ROH segments per animal with a mean of 349.92 Mb and lowest in CGC with a mean of 141.10 Mb. Furthermore, it showed the formation of new haplotypes in CGC around the sixth generation. Selection signatures were evaluated through Fst and HapFlk analyses. Several selection sweeps in CGC were identified especially in chromosomes 5 and 14 which have previously been reported to be associated with coat color and growth traits. The study supports our previous findings that progenitor combinations are not stable over generations and that either direct or natural selection plays a role in modifying the progenitor proportions. Furthermore, the results showed that Tarentaise contributed useful attributes to the composite in a cool semi-arid environment and suggests a re-exploration of this breed's role may be warranted.
Collapse
Affiliation(s)
- E Hay
- USDA Agricultural Research Service, Fort Keogh Livestock and Range Research Laboratory, Miles City, MT 59301, USA
| | - S Toghiani
- USDA Agricultural Research Service, Beltsville Agricultural Research Center, Beltsville, MD, 20705, USA
| | - A J Roberts
- USDA Agricultural Research Service, Fort Keogh Livestock and Range Research Laboratory, Miles City, MT 59301, USA
| | - T Paim
- Instituto Federal de Educação, Ciência e Tecnologia Goiano, Campus Rio Verde, Rio Verde, Goias, Brazil
| | - L A Kuehn
- USDA, Agricultural Research Service, US Meat Animal Research Center, Clay Center, 68933, USA
| | - H D Blackburn
- National Center for Genetic Resources Preservation, USDA, Fort Collins, CO, 80521, USA
| |
Collapse
|
9
|
Chen SY, Schenkel FS, Melo ALP, Oliveira HR, Pedrosa VB, Araujo AC, Melka MG, Brito LF. Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes. BMC Genomics 2022; 23:331. [PMID: 35484513 PMCID: PMC9052698 DOI: 10.1186/s12864-022-08555-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 04/12/2022] [Indexed: 02/06/2023] Open
Abstract
Background Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic. Results There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits. Conclusions These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08555-z.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, 611130, Sichuan, China
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Ana L P Melo
- Department of Reproduction and Animal Evaluation, Rural Federal University of Rio de Janeiro, Seropédica, RJ, 23897-000, Brazil
| | - Hinayah R Oliveira
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Victor B Pedrosa
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Department of Animal Sciences, State University of Ponta Grossa, Ponta Grossa, PR, 84030-900, Brazil
| | - Andre C Araujo
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA
| | - Melkaye G Melka
- Department of Animal and Food Science, University of Wisconsin River Falls, River Falls, WI, 54022, USA
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA. .,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
10
|
Crum TE, Schnabel RD, Decker JE, Taylor JF. Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds. Front Genet 2021; 12:758394. [PMID: 34733318 PMCID: PMC8558500 DOI: 10.3389/fgene.2021.758394] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Accepted: 09/27/2021] [Indexed: 11/14/2022] Open
Abstract
Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.
Collapse
Affiliation(s)
- Tamar E Crum
- Division of Animal Sciences, University of Missouri, Columbia, MO, United States
| | - Robert D Schnabel
- Division of Animal Sciences, University of Missouri, Columbia, MO, United States.,Informatics Institute, University of Missouri, Columbia, MO, United States
| | - Jared E Decker
- Division of Animal Sciences, University of Missouri, Columbia, MO, United States.,Informatics Institute, University of Missouri, Columbia, MO, United States
| | - Jeremy F Taylor
- Division of Animal Sciences, University of Missouri, Columbia, MO, United States
| |
Collapse
|
11
|
Gutierrez-Reinoso MA, Aponte PM, Garcia-Herreros M. Genomic Analysis, Progress and Future Perspectives in Dairy Cattle Selection: A Review. Animals (Basel) 2021; 11:599. [PMID: 33668747 PMCID: PMC7996307 DOI: 10.3390/ani11030599] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2020] [Revised: 11/23/2020] [Accepted: 11/24/2020] [Indexed: 12/16/2022] Open
Abstract
Genomics comprises a set of current and valuable technologies implemented as selection tools in dairy cattle commercial breeding programs. The intensive progeny testing for production and reproductive traits based on genomic breeding values (GEBVs) has been crucial to increasing dairy cattle productivity. The knowledge of key genes and haplotypes, including their regulation mechanisms, as markers for productivity traits, may improve the strategies on the present and future for dairy cattle selection. Genome-wide association studies (GWAS) such as quantitative trait loci (QTL), single nucleotide polymorphisms (SNPs), or single-step genomic best linear unbiased prediction (ssGBLUP) methods have already been included in global dairy programs for the estimation of marker-assisted selection-derived effects. The increase in genetic progress based on genomic predicting accuracy has also contributed to the understanding of genetic effects in dairy cattle offspring. However, the crossing within inbred-lines critically increased homozygosis with accumulated negative effects of inbreeding like a decline in reproductive performance. Thus, inaccurate-biased estimations based on empirical-conventional models of dairy production systems face an increased risk of providing suboptimal results derived from errors in the selection of candidates of high genetic merit-based just on low-heritability phenotypic traits. This extends the generation intervals and increases costs due to the significant reduction of genetic gains. The remarkable progress of genomic prediction increases the accurate selection of superior candidates. The scope of the present review is to summarize and discuss the advances and challenges of genomic tools for dairy cattle selection for optimizing breeding programs and controlling negative inbreeding depression effects on productivity and consequently, achieving economic-effective advances in food production efficiency. Particular attention is given to the potential genomic selection-derived results to facilitate precision management on modern dairy farms, including an overview of novel genome editing methodologies as perspectives toward the future.
Collapse
Affiliation(s)
- Miguel A. Gutierrez-Reinoso
- Facultad de Ciencias Agropecuarias y Recursos Naturales, Carrera de Medicina Veterinaria, Universidad Técnica de Cotopaxi (UTC), Latacunga 05-0150, Ecuador
- Laboratorio de Biotecnología Animal, Departamento de Ciencia Animal, Facultad de Ciencias Veterinarias, Universidad de Concepción (UdeC), Chillán 3780000, Chile
| | - Pedro M. Aponte
- Colegio de Ciencias Biológicas y Ambientales (COCIBA), Universidad San Francisco de Quito (USFQ), Quito 170157, Ecuador
- Campus Cumbayá, Instituto de Investigaciones en Biomedicina “One-health”, Universidad San Francisco de Quito (USFQ), Quito 170157, Ecuador
| | - Manuel Garcia-Herreros
- Instituto Nacional de Investigação Agrária e Veterinária (INIAV), 2005-048 Santarém, Portugal
| |
Collapse
|
12
|
Liu L, Zhou J, Chen CJ, Zhang J, Wen W, Tian J, Zhang Z, Gu Y. GWAS-Based Identification of New Loci for Milk Yield, Fat, and Protein in Holstein Cattle. Animals (Basel) 2020; 10:E2048. [PMID: 33167458 PMCID: PMC7694478 DOI: 10.3390/ani10112048] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2020] [Revised: 11/01/2020] [Accepted: 11/03/2020] [Indexed: 12/20/2022] Open
Abstract
High-yield and high-quality of milk are the primary goals of dairy production. Understanding the genetic architecture underlying these milk-related traits is beneficial so that genetic variants can be targeted toward the genetic improvement. In this study, we measured five milk production and quality traits in Holstein cattle population from China. These traits included milk yield, fat, and protein. We used the estimated breeding values as dependent variables to conduct the genome-wide association studies (GWAS). Breeding values were estimated through pedigree relationships by using a linear mixed model. Genotyping was carried out on the individuals with phenotypes by using the Illumina BovineSNP150 BeadChip. The association analyses were conducted by using the fixed and random model Circulating Probability Unification (FarmCPU) method. A total of ten single-nucleotide polymorphisms (SNPs) were detected above the genome-wide significant threshold (p < 4.0 × 10-7), including six located in previously reported quantitative traits locus (QTL) regions. We found eight candidate genes within distances of 120 kb upstream or downstream to the associated SNPs. The study not only identified the effect of DGAT1 gene on milk fat and protein, but also discovered novel genetic loci and candidate genes related to milk traits. These novel genetic loci would be an important basis for molecular breeding in dairy cattle.
Collapse
Affiliation(s)
- Liyuan Liu
- School of Agriculture, Ningxia University, Yinchuan 750021, Ningxia, China; (L.L.); (J.Z.); (J.Z.)
- Department of Crop and Soil Sciences, Washington State University, Pullman, Washington, DC 99164, USA;
| | - Jinghang Zhou
- School of Agriculture, Ningxia University, Yinchuan 750021, Ningxia, China; (L.L.); (J.Z.); (J.Z.)
- Department of Crop and Soil Sciences, Washington State University, Pullman, Washington, DC 99164, USA;
| | - Chunpeng James Chen
- Department of Crop and Soil Sciences, Washington State University, Pullman, Washington, DC 99164, USA;
| | - Juan Zhang
- School of Agriculture, Ningxia University, Yinchuan 750021, Ningxia, China; (L.L.); (J.Z.); (J.Z.)
| | - Wan Wen
- Animal Husbandry Workstation, Yinchuan 750001, Ningxia, China; (W.W.); (J.T.)
| | - Jia Tian
- Animal Husbandry Workstation, Yinchuan 750001, Ningxia, China; (W.W.); (J.T.)
| | - Zhiwu Zhang
- Department of Crop and Soil Sciences, Washington State University, Pullman, Washington, DC 99164, USA;
| | - Yaling Gu
- School of Agriculture, Ningxia University, Yinchuan 750021, Ningxia, China; (L.L.); (J.Z.); (J.Z.)
| |
Collapse
|
13
|
Galliou JM, Kiser JN, Oliver KF, Seabury CM, Moraes JGN, Burns GW, Spencer TE, Dalton J, Neibergs HL. Identification of Loci and Pathways Associated with Heifer Conception Rate in U.S. Holsteins. Genes (Basel) 2020; 11:genes11070767. [PMID: 32650431 PMCID: PMC7397161 DOI: 10.3390/genes11070767] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Revised: 06/27/2020] [Accepted: 07/06/2020] [Indexed: 12/15/2022] Open
Abstract
Heifer conception rate (HCR) is defined as the percentage of inseminated heifers that become pregnant at each service. The genome-wide association analyses in this study focused on identifying the loci associated with Holstein heifer (n = 2013) conception rate at first service (HCR1) and the number of times bred (TBRD) to achieve a pregnancy. There were 348 unique loci associated (p < 5 × 10−8) with HCR1 and 615 unique loci associated (p < 5 × 10−8) with TBRD. The two phenotypes shared 302 loci, and 56 loci were validated in independent cattle populations. There were 52 transcription factor binding sites (TFBS) and 552 positional candidate genes identified in the HCR1- and TBRD-associated loci. The positional candidate genes and the TFBS associated with HCR1 and TBRD were used in the ingenuity pathway analysis (IPA). In the IPA, 11 pathways, 207 master regulators and 11 upstream regulators were associated (p < 1.23 × 10−5) with HCR1 and TBRD. The validated loci associated with both HCR1 and TBRD make good candidates for genomic selection and further investigations to elucidate the mechanisms associated with subfertility and infertility.
Collapse
Affiliation(s)
- Justine M. Galliou
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA 99164, USA; (J.M.G.); (J.N.K.); (K.F.O.)
| | - Jennifer N. Kiser
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA 99164, USA; (J.M.G.); (J.N.K.); (K.F.O.)
| | - Kayleen F. Oliver
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA 99164, USA; (J.M.G.); (J.N.K.); (K.F.O.)
| | - Christopher M. Seabury
- Department of Veterinary Pathobiology, College of Veterinary Medicine, Texas A&M University, College Station, TX 77843, USA;
| | - Joao G. N. Moraes
- Division of Animal Sciences, University of Missouri, Columbia, MO 65211, USA; (J.G.N.M.); (G.W.B.); (T.E.S.)
| | - Gregory W. Burns
- Division of Animal Sciences, University of Missouri, Columbia, MO 65211, USA; (J.G.N.M.); (G.W.B.); (T.E.S.)
| | - Thomas E. Spencer
- Division of Animal Sciences, University of Missouri, Columbia, MO 65211, USA; (J.G.N.M.); (G.W.B.); (T.E.S.)
| | - Joseph Dalton
- Department of Animal and Veterinary Sciences, University of Idaho, Caldwell, ID 83844, USA;
| | - Holly L. Neibergs
- Department of Animal Sciences and Center for Reproductive Biology, Washington State University, Pullman, WA 99164, USA; (J.M.G.); (J.N.K.); (K.F.O.)
- Correspondence:
| |
Collapse
|
14
|
Sigdel A, Liu L, Abdollahi-Arpanahi R, Aguilar I, Peñagaricano F. Genetic dissection of reproductive performance of dairy cows under heat stress. Anim Genet 2020; 51:511-520. [PMID: 32363588 PMCID: PMC7383808 DOI: 10.1111/age.12943] [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] [Accepted: 03/31/2020] [Indexed: 02/06/2023]
Abstract
Heat stress negatively impacts the reproductive performance of dairy cows. The main objective of this study was to dissect the genetic basis underlying dairy cow fertility under heat stress conditions. Our first goal was to estimate genetic components of cow conception across lactations considering heat stress. Our second goal was to reveal individual genes and functional gene‐sets that explain a cow’s ability to conceive under thermal stress. Data consisted of 74 221 insemination records on 13 704 Holstein cows. Multitrait linear repeatability test‐day models with random regressions on a function of temperature–humidity index values were used for the analyses. Heritability estimates for cow conception under heat stress were around 2–3%, whereas genetic correlations between general and thermotolerance additive genetic effects were negative and ranged between −0.35 and −0.82, indicating an unfavorable relationship between cows’ ability to conceive under thermo‐neutral vs. thermo‐stress conditions. Whole‐genome scans identified at least six genomic regions on BTA1, BTA10, BTA11, BTA17, BTA21 and BTA23 associated with conception under thermal stress. These regions harbor candidate genes such as BRWD1, EXD2, ADAM20, EPAS1, TAOK3, and NOS1, which are directly implicated in reproductive functions and cellular response to heat stress. The gene‐set enrichment analysis revealed functional terms related to fertilization, developmental biology, heat shock proteins and oxidative stress, among others. Overall, our findings contribute to a better understanding of the genetics underlying the reproductive performance of dairy cattle under heat stress conditions and point out novel genomic strategies for improving thermotolerance and fertility via marker‐assisted breeding.
Collapse
Affiliation(s)
- A Sigdel
- Department of Animal Sciences, University of Florida, Gainesville, FL, 32611, USA
| | - L Liu
- Department of Animal Sciences, University of Florida, Gainesville, FL, 32611, USA
| | - R Abdollahi-Arpanahi
- Department of Animal Sciences, University of Florida, Gainesville, FL, 32611, USA
| | - I Aguilar
- Instituto Nacional de Investigación Agropecuaria, Montevideo, 11100, Uruguay
| | - F Peñagaricano
- Department of Animal Sciences, University of Florida, Gainesville, FL, 32611, USA.,University of Florida Genetics Institute, University of Florida, Gainesville, FL, 32611, USA
| |
Collapse
|
15
|
Mastrangelo S, Ben Jemaa S, Ciani E, Sottile G, Moscarelli A, Boussaha M, Montedoro M, Pilla F, Cassandro M. Genome-wide detection of signatures of selection in three Valdostana cattle populations. J Anim Breed Genet 2020; 137:609-621. [PMID: 32219904 DOI: 10.1111/jbg.12476] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Revised: 02/28/2020] [Accepted: 03/08/2020] [Indexed: 12/15/2022]
Abstract
The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (FST ) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock.
Collapse
Affiliation(s)
- Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Slim Ben Jemaa
- Laboratoire des Productions Animales et Fourragères, Institut National de la Recherche Agronomique de Tunisie, Université de Carthage, Ariana, Tunisia
| | - Elena Ciani
- Dipartimento di Bioscienze Biotecnologie e Biofarmaceutica, University of Bari, Bari, Italy
| | - Gianluca Sottile
- Dipartimento Scienze Economiche, Aziendali e Statistiche, University of Palermo, Palermo, Italy
| | - Angelo Moscarelli
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Mekki Boussaha
- INRAE, AgroParisTech, GABI, Université Paris Saclay, Jouy-en-Josas, France
| | - Marina Montedoro
- Istituto Sperimentale Italiano Lazzaro Spallanzani, Rivolta d'Adda, Italy
| | - Fabio Pilla
- Dipartimento di Agricoltura, Ambiente e Alimenti, University of Molise, Campobasso, Italy
| | - Martino Cassandro
- Dipartimento di Agronomia Animali Alimenti Risorse naturali e Ambiente, University of Padova, Legnaro, Italy
| |
Collapse
|