1
|
Kafi Z, Masoudi AA, Torshizi RV, Ehsani A. Copy number variations affecting growth curve parameters in a crossbred chicken population. Gene 2024; 927:148710. [PMID: 38901536 DOI: 10.1016/j.gene.2024.148710] [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: 12/31/2023] [Revised: 06/01/2024] [Accepted: 06/17/2024] [Indexed: 06/22/2024]
Abstract
Copy number variations (CNVs) are key structural variations in the genome and may contribute to phenotypic differences. In this study, we used a F2 chicken population created from reciprocal crossing between fast-growing Arian broiler line and Urmia native chickens. The chickens were genotyped by 60 K SNP BeadChip, and PennCNV algorithm was used to detect genome-wide CNVs. The growth curve parameters of W0, k, L, Wf, Wi, ti and average GR were used as phenotypic data. The association between CNV and growth curve parameters was carried out using the CNVRanger R/Bioconductor package. Five CNV regions (CNVRs) were chosen for the validation experiment using qPCR. Gene enrichment analysis was done using WebGestalt. The STRING database was used to search for significant pathways. The results identified 966 CNVs and 600 CNVRs including 468 gains, 67 losses, and 65 both events on autosomal chromosomes. Validation of the CNVRs obtained from the qPCR assay were 79 % consistent with the prediction by PennCNV. A total of 43 significant CNVs were obtained for the seven growth curve parameters. The 416 genes annotated for significant CNVs. Six genes out of 416 genes were most related to growth curve parameters. These genes were LCP2, Dock2, CD80, CYFIP1, NIPA1 and NIPA2. Some of these genes in their biological process were associated with the growth, reproduction and development of cells or organs that ultimately lead to the growth of the body. The results of the study could pave the way for better understanding the molecular process of CNVs and growth curve parameters in birds.
Collapse
Affiliation(s)
- Zeinab Kafi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| | - Ali Akbar Masoudi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| | - Rasoul Vaez Torshizi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| | - Alireza Ehsani
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| |
Collapse
|
2
|
Lotfizadeh F, Masoudi AA, Vaez Torshizi R, Emrani H. Genome-wide association study of copy number variations with shank traits in a F 2 crossbred chicken population. Anim Genet 2024; 55:559-574. [PMID: 38764135 DOI: 10.1111/age.13447] [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: 10/02/2023] [Revised: 04/21/2024] [Accepted: 05/02/2024] [Indexed: 05/21/2024]
Abstract
Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.
Collapse
Affiliation(s)
- Fateme Lotfizadeh
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran
| | - Ali Akbar Masoudi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran
| | - Rasoul Vaez Torshizi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran
| | - Hossein Emrani
- Animal Science Research Institute of Iran, Agricultural Research, Education and Extension Organization, Karaj, Iran
| |
Collapse
|
3
|
Haqani MI, Nakano M, Nagano AJ, Nakamura Y, Tsudzuki M. Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing. Sci Rep 2023; 13:21307. [PMID: 38042890 PMCID: PMC10693557 DOI: 10.1038/s41598-023-48293-0] [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: 01/06/2023] [Revised: 10/10/2023] [Accepted: 11/24/2023] [Indexed: 12/04/2023] Open
Abstract
This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.
Collapse
Affiliation(s)
- Mohammad Ibrahim Haqani
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashi-Hiroshima, Hiroshima, 739-8525, Japan.
| | - Michiharu Nakano
- Faculty of Agriculture and Marine Sciences, Kochi University, Nankoku, Kochi, 783-8502, Japan
| | - Atsushi J Nagano
- Faculty of Agriculture, Ryukoku University, Otsu, Shiga, 520-2194, Japan
- Institute for Advanced Biosciences, Keio University, Yamagata, 997-0017, Japan
| | - Yoshiaki Nakamura
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashi-Hiroshima, Hiroshima, 739-8525, Japan
- Japanese Avian Bioresource Project Research Center, Hiroshima University, Higashi-Hiroshima, Hiroshima, 739-8525, Japan
| | - Masaoki Tsudzuki
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashi-Hiroshima, Hiroshima, 739-8525, Japan.
- Japanese Avian Bioresource Project Research Center, Hiroshima University, Higashi-Hiroshima, Hiroshima, 739-8525, Japan.
| |
Collapse
|
4
|
Phenotypic Analysis of Growth and Morphological Traits in Miniature Breeds of Japanese Indigenous Chickens. J Poult Sci 2022; 59:38-47. [PMID: 35125911 PMCID: PMC8791770 DOI: 10.2141/jpsa.0200110] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Accepted: 04/19/2021] [Indexed: 11/21/2022] Open
|
5
|
Genome-Wide Association Study Identifies 12 Loci Associated with Body Weight at Age 8 Weeks in Korean Native Chickens. Genes (Basel) 2021; 12:genes12081170. [PMID: 34440344 PMCID: PMC8394794 DOI: 10.3390/genes12081170] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Revised: 07/24/2021] [Accepted: 07/26/2021] [Indexed: 02/07/2023] Open
Abstract
Meat from Korean native chickens (KNCs) has high consumer demand; however, slow growth performance and high variation in body weight (BW) of KNCs remain an issue. Genome-wide association study (GWAS) is a powerful method to identify quantitative trait-associated genomic loci. A GWAS, based on a large-scale KNC population, is needed to identify underlying genetic mechanisms related to its growth traits. To identify BW-associated genomic regions, we performed a GWAS using the chicken 60K single nucleotide polymorphism (SNP) panel for 1328 KNCs. BW was measured at 8 weeks of age, from 2018 to 2020. Twelve SNPs were associated with BW at the suggestive significance level (p < 2.95 × 10−5) and located near or within 11 candidate genes, including WDR37, KCNIP4, SLIT2, PPARGC1A, MYOCD and ADGRA3. Gene set enrichment analysis based on the GWAS results at p < 0.05 (1680 SNPs) showed that 32 Gene Ontology terms and two Kyoto Encyclopedia of Genes and Genomes pathways, including regulation of transcription, motor activity, the mitogen-activated protein kinase signaling pathway, and tight junction, were significantly enriched (p < 0.05) for BW-associated genes. These pathways are involved in cell growth and development, related to BW gain. The identified SNPs are potential biomarkers in KNC breeding.
Collapse
|
6
|
Marchesi JAP, Ono RK, Cantão ME, Ibelli AMG, Peixoto JDO, Moreira GCM, Godoy TF, Coutinho LL, Munari DP, Ledur MC. Exploring the genetic architecture of feed efficiency traits in chickens. Sci Rep 2021; 11:4622. [PMID: 33633287 PMCID: PMC7907133 DOI: 10.1038/s41598-021-84125-9] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Accepted: 02/12/2021] [Indexed: 11/09/2022] Open
Abstract
Chicken feed efficiency (FE) traits are the most important economic traits in broiler production. Several studies evaluating genetic factors affecting food consumption in chickens are available. However, most of these studies identified genomic regions containing putative quantitative trait loci for each trait separately. It is still a challenge to find common gene networks related to these traits. Therefore, here, a genome-wide association study (GWAS) was conducted to explore candidate genomic regions responsible for Feed Intake (FI), Body Weight Gain (BWG) and Feed Conversion Ratio (FCR) traits and their gene networks. A total of 1430 broilers from an experimental population was genotyped with the high density Affymetrix 600K SNP array. A total of 119 associated SNPs located in 20 chromosomes were identified, where some of them were common in more than one FE trait. In addition, novel genomic regions were prospected considering the SNPs dominance effects and sex interaction, identifying putative candidate genes only when these effects were fit in the model. Relevant candidate genes such as ATRNL1, PIK3C2A, PTPRN2, SORCS3 and gga-mir-1759 were highlighted in this study helping to elucidate the genomic architecture of feed efficiency traits. These results provide new insights on the mechanisms underlying the consumption and utilization of food in chickens.
Collapse
Affiliation(s)
- Jorge Augusto Petroli Marchesi
- Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista "Júlio de Mesquita Filho", Jaboticabal, SP, 14884-900, Brazil.,Departamento de Genética, Universidade de São Paulo, Ribeirão Preto, SP, 14049-900, Brazil
| | - Rafael Keith Ono
- Embrapa Suínos e Aves, Concórdia, SC, 89715-899, Brazil.,Pamplona Alimentos S/A, Rio do Sul, SC, 89164-900, Brazil
| | | | | | | | - Gabriel Costa Monteiro Moreira
- Departamento de Zootecnia, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Av. Pádua Dias 11, Piracicaba, SP, 13419-900, Brazil
| | - Thaís Fernanda Godoy
- Departamento de Zootecnia, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Av. Pádua Dias 11, Piracicaba, SP, 13419-900, Brazil
| | - Luiz Lehmann Coutinho
- Departamento de Zootecnia, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Av. Pádua Dias 11, Piracicaba, SP, 13419-900, Brazil
| | - Danísio Prado Munari
- Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista "Júlio de Mesquita Filho", Jaboticabal, SP, 14884-900, Brazil
| | | |
Collapse
|
7
|
Seifi Moroudi R, Ansari Mahyari S, Vaez Torshizi R, Lanjanian H, Masoudi-Nejad A. Identification of new genes and quantitative trait locis associated with growth curve parameters in F2 chicken population using genome-wide association study. Anim Genet 2021; 52:171-184. [PMID: 33428266 DOI: 10.1111/age.13038] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/26/2020] [Indexed: 11/30/2022]
Abstract
The markers which are correlated with the growth curve parameters help in understanding the characteristics of individual growth during the rearing of livestock. This study aimed to identify a set of biomarkers through a GWAS for growth curve parameters in crossbred chickens using the Illumnia 60K chicken SNP Beadchip. Growth data were collected from a total of 301 birds from cross of a broiler line and native chickens. Using the Gompertz-Laird model, two growth curve parameters, the instantaneous growth rate per day (L) and the coefficient of relative growth or maturing index (k), were estimated. The L and k were used to estimate five derived parameters, namely asymptotic (mature) body weight, body weight at inflection point, age at the inflection point, average growth rate and maximum growth rate. These parameters were considered as phenotypic values in the GWAS based on generalized linear models. The results of the GWAS indicated 21 significant markers, which were located near or within 46 genes. A number of these genes, such as GH, RET, GRB14, FTSJ3 and CCK, are important for growth and meat quality in chickens, and some of them are growth related in other species such as sheep and cattle (GPI, XIRP2, GALNTL6, BMS1, THSD4, TRHDE, SHISA9, ACSL6 and DYNC1LI2). The other genes are associated with developmental biological pathways. These genes are particuarly related to body weight, average daily gain and growth QTL. The results of this study can shed light on the genetic mechanism of biological functions of growth factors in broiler chickens, which is useful for developing management practices and accelerating genetic progress in breeding programs.
Collapse
Affiliation(s)
- R Seifi Moroudi
- Department of Animal Science, College of Agriculture, Isfahan University of Technology, PO Box 841583111, Isfahan, Iran
| | - S Ansari Mahyari
- Department of Animal Science, College of Agriculture, Isfahan University of Technology, PO Box 841583111, Isfahan, Iran
| | - R Vaez Torshizi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, 14115-336, Iran
| | - H Lanjanian
- Laboratory of Systems Biology and Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, 1417614411, Iran
| | - A Masoudi-Nejad
- Laboratory of Systems Biology and Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, 1417614411, Iran
| |
Collapse
|
8
|
Asgari Z, Ehsani A, Masoudi AA, Vaez Torshizi R. Bayes factors revealed selection signature for time to market body weight in chicken: a genome-wide association study using BayesCpi methodology. ITALIAN JOURNAL OF ANIMAL SCIENCE 2021. [DOI: 10.1080/1828051x.2021.1965920] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Zeinab Asgari
- Department of Animal Science, Tarbiat Modares University, Tehran, Iran
| | - Alireza Ehsani
- Department of Animal Science, Tarbiat Modares University, Tehran, Iran
| | - Ali Akbar Masoudi
- Department of Animal Science, Tarbiat Modares University, Tehran, Iran
| | | |
Collapse
|
9
|
Emrani H, Masoudi AA, Vaez Torshizi R, Ehsani A. Genome-wide association study of shank length and diameter at different developmental stages in chicken F2 resource population. Anim Genet 2020; 51:722-730. [PMID: 32662094 DOI: 10.1111/age.12981] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/15/2020] [Indexed: 01/09/2023]
Abstract
In order to find SNPs and genes affecting shank traits, we performed a GWAS in a chicken F2 population of eight half-sib families from five hatches derived from reciprocal crosses between an Arian fast-growing line and an Urmia indigenous slow-growing chicken. A total of 308 birds were genotyped using a 60K chicken SNP chip. Shank traits including shank length and diameter were measured weekly from birth to 12 weeks of age. A generalized linear model and a compressed mixed linear model (CMLM) were applied to achieve the significant regions. The value of the average genomic inflation factor (λ statistic) of the CMLM model (0.99) indicated that the CMLM was more effective than the generalized linear model in controlling the population structure. The genes surrounding significant SNPs and their biological functions were identified from NCBI, Ensembl and UniProt databases. The results indicated that 12 SNPs at 12 different ages passed the LD-adjusted 5% Bonferroni significant threshold. Two SNPs were significant for shank length and nine SNPs were significant for shank diameter. The significant SNPs were located near to or inside 11 candidate genes. The results showed that a number of significant SNPs in the middle ages were higher than the rest. The MXRA8 gene was related to the significant SNP at week 1 that promotes proliferation of growth plate chondrocytes. A unique SNP of Gga_rs16689511 located on chicken Z chromosome within the LOC101747628 gene was related to shank length at three different ages of birds (weeks 8, 9 and 11). The significant SNPs for shank diameter were found at weeks 4 and 7 (four and five SNPs respectively). The identifications of SNPs and genes here could contribute to a better understanding of the genetic control of shank traits in chicken.
Collapse
Affiliation(s)
- H Emrani
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, PO Box 14115-336, Tehran, Iran
| | - A A Masoudi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, PO Box 14115-336, Tehran, Iran
| | - R Vaez Torshizi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, PO Box 14115-336, Tehran, Iran
| | - A Ehsani
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, PO Box 14115-336, Tehran, Iran
| |
Collapse
|
10
|
Hlongwane NL, Hadebe K, Soma P, Dzomba EF, Muchadeyi FC. Genome Wide Assessment of Genetic Variation and Population Distinctiveness of the Pig Family in South Africa. Front Genet 2020; 11:344. [PMID: 32457791 PMCID: PMC7221027 DOI: 10.3389/fgene.2020.00344] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2018] [Accepted: 03/23/2020] [Indexed: 12/12/2022] Open
Abstract
Genetic diversity is of great importance and a prerequisite for genetic improvement and conservation programs in pigs and other livestock populations. The present study provides a genome wide analysis of the genetic variability and population structure of pig populations from different production systems in South Africa relative to global populations. A total of 234 pigs sampled in South Africa and consisting of village (n = 91), commercial (n = 60), indigenous (n = 40), Asian (n = 5) and wild (n = 38) populations were genotyped using Porcine SNP60K BeadChip. In addition, 389 genotypes representing village and commercial pigs from America, Europe, and Asia were accessed from a previous study and used to compare population clustering and relationships of South African pigs with global populations. Moderate heterozygosity levels, ranging from 0.204 for Warthogs to 0.371 for village pigs sampled from Capricorn municipality in Eastern Cape province of South Africa were observed. Principal Component Analysis of the South African pigs resulted in four distinct clusters of (i) Duroc; (ii) Vietnamese; (iii) Bush pig and Warthog and (iv) a cluster with the rest of the commercial (SA Large White and Landrace), village, Wild Boar and indigenous breeds of Koelbroek and Windsnyer. The clustering demonstrated alignment with genetic similarities, geographic location and production systems. The PCA with the global populations also resulted in four clusters that where populated with (i) all the village populations, wild boars, SA indigenous and the large white and landraces; (ii) Durocs (iii) Chinese and Vietnamese pigs and (iv) Warthog and Bush pig. K = 10 (The number of population units) was the most probable ADMIXTURE based clustering, which grouped animals according to their populations with the exception of the village pigs that showed presence of admixture. AMOVA reported 19.92%-98.62% of the genetic variation to be within populations. Sub structuring was observed between South African commercial populations as well as between Indigenous and commercial breeds. Population pairwise F ST analysis showed genetic differentiation (P ≤ 0.05) between the village, commercial and wild populations. A per marker per population pairwise F ST analysis revealed SNPs associated with QTLs for traits such as meat quality, cytoskeletal and muscle development, glucose metabolism processes and growth factors between both domestic populations as well as between wild and domestic breeds. Overall, the study provided a baseline understanding of porcine diversity and an important foundation for porcine genomics of South African populations.
Collapse
Affiliation(s)
- Nompilo Lucia Hlongwane
- Biotechnology Platform, Agricultural Research Council, Onderstepoort, South Africa
- Discipline of Genetics, School of Life Sciences, University of KwaZulu-Natal, Pietermartizburg, South Africa
| | - Khanyisile Hadebe
- Biotechnology Platform, Agricultural Research Council, Onderstepoort, South Africa
| | - Pranisha Soma
- Animal Production Institute, Agricultural Research Council, Irene, South Africa
| | - Edgar Farai Dzomba
- Discipline of Genetics, School of Life Sciences, University of KwaZulu-Natal, Pietermartizburg, South Africa
| | | |
Collapse
|