TY - JOUR
T1 - A genome-wide analysis of copy number variation in Murciano-Granadina goats
AU - Guan, Dailu
AU - Martínez, Amparo
AU - Castelló, Anna
AU - Landi, Vincenzo
AU - Luigi-Sierra, María Gracia
AU - Fernández-Álvarez, Javier
AU - Cabrera, Betlem
AU - Delgado, Juan Vicente
AU - Such, Xavier
AU - Jordana, Jordi
AU - Amills, Marcel
PY - 2020/8/8
Y1 - 2020/8/8
N2 - Background: In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions.Results: Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 Bead-Chip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 x 10(-10)), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 x 10(-04)) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 x 10(-03)).Conclusions: A previous study reported that the average number of CNVR per goat breed was similar to 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.
AB - Background: In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions.Results: Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 Bead-Chip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 x 10(-10)), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 x 10(-04)) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 x 10(-03)).Conclusions: A previous study reported that the average number of CNVR per goat breed was similar to 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.
KW - ASIP GENE
KW - BLOOD-FLOW
KW - HIDDEN-MARKOV MODEL
KW - INSULIN
KW - MAP
KW - SUSCEPTIBILITY
KW - WHITE
UR - http://www.scopus.com/inward/record.url?scp=85089261828&partnerID=8YFLogxK
U2 - 10.1186/s12711-020-00564-4
DO - 10.1186/s12711-020-00564-4
M3 - Artículo
C2 - 32770942
AN - SCOPUS:85089261828
SN - 1297-9686
VL - 52
SP - 44
JO - Genetics Selection And Evolution
JF - Genetics Selection And Evolution
IS - 1
M1 - 44
ER -