Copy number variants in locally raised Chinese chicken genomes determined using array comparative genomic hybridization

Genome-wide investigation to assess copy number variants in the Italian local chicken population

BACKGROUND

Copy number variants (CNV) hold significant functional and evolutionary importance. Numerous ongoing CNV studies aim to elucidate the etiology of human diseases and gain insights into the population structure of livestock. High-density chips have enabled the detection of CNV with increased resolution, leading to the identification of even small CNV. This study aimed to identify CNV in local Italian chicken breeds and investigate their distribution across the genome.

RESULTS

Copy number variants were mainly distributed across the first six chromosomes and primarily associated with loss type CNV. The majority of CNV in the investigated breeds were of types 0 and 1, and the minimum length of CNV was significantly larger than that reported in previous studies. Interestingly, a high proportion of the length of chromosome 16 was covered by copy number variation regions (CNVR), with the major histocompatibility complex being the likely cause. Among the genes identified within CNVR, only those present in at least five animals across breeds (n = 95) were discussed to reduce the focus on redundant CNV. Some of these genes have been associated to functional traits in chickens. Notably, several CNVR on different chromosomes harbor genes related to muscle development, tissue-specific biological processes, heat stress resistance, and immune response. Quantitative trait loci (QTL) were also analyzed to investigate potential overlapping with the identified CNVR: 54 out of the 95 gene-containing regions overlapped with 428 QTL associated to body weight and size, carcass characteristics, egg production, egg components, fat deposition, and feed intake.

CONCLUSIONS

The genomic phenomena reported in this study that can cause changes in the distribution of CNV within the genome over time and the comparison of these differences in CNVR of the local chicken breeds could help in preserving these genetic resources.

Population Genomic Sequencing Delineates Global Landscape of Copy Number Variations that Drive Domestication and Breed Formation of in Chicken

Copy number variation (CNV) is an important genetic mechanism that drives evolution and generates new phenotypic variations. To explore the impact of CNV on chicken domestication and breed shaping, the whole-genome CNVs were detected via multiple methods. Using the whole-genome sequencing data from 51 individuals, corresponding to six domestic breeds and wild red jungle fowl (RJF), we determined 19,329 duplications and 98,736 deletions, which covered 11,123 copy number variation regions (CNVRs) and 2,636 protein-coding genes. The principal component analysis (PCA) showed that these individuals could be divided into four populations according to their domestication and selection purpose. Seventy-two highly duplicated CNVRs were detected across all individuals, revealing pivotal roles of nervous system (NRG3, NCAM2), sensory (OR), and follicle development (VTG2) in chicken genome. When contrasting the CNVs of domestic breeds to those of RJFs, 235 CNVRs harboring 255 protein-coding genes, which were predominantly involved in pathways of nervous, immunity, and reproductive system development, were discovered. In breed-specific CNVRs, some valuable genes were identified, including HOXB7 for beard trait in Beijing You chicken; EDN3, SLMO2, TUBB1, and GFPT1 for melanin deposition in Silkie chicken; and SORCS2 for aggressiveness in Luxi Game fowl. Moreover, CSMD1 and NTRK3 with high duplications found exclusively in White Leghorn chicken, and POLR3H, MCM9, DOCK3, and AKR1B1L found in Recessive White Rock chicken may contribute to high egg production and fast-growing traits, respectively. The candidate genes of breed characteristics are valuable resources for further studies on phenotypic variation and the artificial breeding of chickens.

Assessment of linkage disequilibrium patterns between structural variants and single nucleotide polymorphisms in three commercial chicken populations

BACKGROUND

Structural variants (SV) are causative for some prominent phenotypic traits of livestock as different comb types in chickens or color patterns in pigs. Their effects on production traits are also increasingly studied. Nevertheless, accurately calling SV remains challenging. It is therefore of interest, whether close-by single nucleotide polymorphisms (SNPs) are in strong linkage disequilibrium (LD) with SVs and can serve as markers. Literature comes to different conclusions on whether SVs are in LD to SNPs on the same level as SNPs to other SNPs. The present study aimed to generate a precise SV callset from whole-genome short-read sequencing (WGS) data for three commercial chicken populations and to evaluate LD patterns between the called SVs and surrounding SNPs. It is thereby the first study that assessed LD between SVs and SNPs in chickens.

RESULTS

The final callset consisted of 12,294,329 bivariate SNPs, 4,301 deletions (DEL), 224 duplications (DUP), 218 inversions (INV) and 117 translocation breakpoints (BND). While average LD between DELs and SNPs was at the same level as between SNPs and SNPs, LD between other SVs and SNPs was strongly reduced (DUP: 40%, INV: 27%, BND: 19% of between-SNP LD). A main factor for the reduced LD was the presence of local minor allele frequency differences, which accounted for 50% of the difference between SNP - SNP and DUP - SNP LD. This was potentially accompanied by lower genotyping accuracies for DUP, INV and BND compared with SNPs and DELs. An evaluation of the presence of tag SNPs (SNP in highest LD to the variant of interest) further revealed DELs to be slightly less tagged by WGS SNPs than WGS SNPs by other SNPs. This difference, however, was no longer present when reducing the pool of potential tag SNPs to SNPs located on four different chicken genotyping arrays.

CONCLUSIONS

The results implied that genomic variance due to DELs in the chicken populations studied can be captured by different SNP marker sets as good as variance from WGS SNPs, whereas separate SV calling might be advisable for DUP, INV, and BND effects.

Transcript expression profiling of fibromelanosis-related genes in black-bone chickens

1. The aim of the present study was to identify differentially expressed genes (DEGs) and metabolic pathways involved in this phenotype. Fibromelanosis is the most striking feature of black-bone chickens, such as the Silkie and Dongxiang indigenous breeds. Due to the accumulation of eumelanin in connective tissues, fibromelanosis manifests as black colouration of the skin, muscles, gut, and periosteum. Studies on fibromelanosis can provide useful information pertaining to human diseases and offer commercial value to the poultry industry. However, the genetic basis of fibromelanosis remains unclear.2. Digital gene expression analysis was performed on black and white skin samples collected from the HW1 black-bone chicken line to detect differences in genome-wide expression patterns. A total of >30 billion bp were sequenced, and 2,707,926,466 bp and 2,948,782,964 bp of clean data obtained for creation of libraries for black and white skin, respectively. In total, 252 DEGs from 15,508 mapped genes were identified with 83 up-regulated in white skin and 169 up-regulated in black skin.3. Gene ontology analysis highlighted that genes from the extracellular region and associated components were abundant among the DEGs. Pathway analysis revealed that many DEGs were linked to amino acid metabolism and the immune system. qRT-PCR validation using 14 genes showed good conformity with the sequence analysis of fibromelanosis-related genes.4. The results showed that L-dopachrometautomerase precursor (DCT), tyrosine aminotransferase (TAT), 4-hydroxyphenylpyruvate dioxygenase (HPD) from the tyrosine metabolism pathway, coagulation factor II (F2), fibrinogen beta chain (FGB), plasminogen (PLG) and complement component 7 (C7) from the complement and coagulation cascades were important genes in the fibromelanosis process in black-bone chickens. These candidate genes require further correlation analysis and functional verification.

Genetics and breeding of a black-bone and blue eggshell chicken line. 1. Body weight, skin color, and their combined selection

An experimental population of chickens was developed from the cross between 2 indigenous Chinese breeds, Dongxiang blue eggshell and Jiangshan black-bone. This breeding was aimed at eventually combining dark heavy black-bone body and blue eggshell, into a single dual-purpose breed. BW was recorded and skin L∗, a∗, and b∗ color parameters were measured by a Chroma Meter at several ages (56, 105, 150, 200, 250, and 300 d). At 250 d, 3 independent observers classified skin darkness using a 3-level visual scale (1 = light, 2 = intermediate, 3 = dark). The 7-level average visual skin darkness, calculated for each chicken, was highly correlated (−0.658 and −0.612 in females and males, respectively) with skin L∗ (lightness), indicating that the accurately measured L∗ is reliable and useful reverse expression of visual skin darkness of black-bone chickens. Mean BW and skin L∗ of both sexes increased with age, to 2,063 and 1,522 g in males and females, respectively, at 300 d, and to 63 and 55 L∗ units in males and females, respectively, at 250 d. The population's full-pedigree allowed estimating heritability and genetic correlations between traits. The heritability estimates of BW were similar in both sexes, increasing from around 0.25 at 56 d, to 0.53 to 0.60 at 150 d, and 0.57 to 0.62 at 300 d. Over the 5 ages, heritability estimates of skin L∗ were moderate to high, ranging from 0.45 to 0.58 in females, and from 0.31 to 0.65 in males, and the genetic correlations between BW and L∗ ranged mostly from 0.20 to 0.45. These low-to-moderate correlations between high BW and high L∗ (low darkness) are unfavorable; hence they were combined into an index, standardized BW minus standardized L∗, allowing future selection for high BW with low L∗. With high heritability of this index, 0.487 (females at 300 d) and 0.410 to 0.555 (males at 150 d or older), simultaneous improvements in BW and skin darkness appear to be feasible. The methodology used in this study can be useful in chicken populations experimentally bred for combination of high BW and other body characteristics.

Genome-wide characterization of copy number variations in the host genome in genetic resistance to Marek's disease using next generation sequencing

Background

Marek’s disease (MD) is a highly neoplastic disease primarily affecting chickens, and remains as a chronic infectious disease that threatens the poultry industry. Copy number variation (CNV) has been examined in many species and is recognized as a major source of genetic variation that directly contributes to phenotypic variation such as resistance to infectious diseases. Two highly inbred chicken lines, 6₃ (MD-resistant) and 7₂ (MD-susceptible), as well as their F₁ generation and six recombinant congenic strains (RCSs) with varied susceptibility to MD, are considered as ideal models to identify the complex mechanisms of genetic and molecular resistance to MD.

Results

In the present study, to unravel the potential genetic mechanisms underlying resistance to MD, we performed a genome-wide CNV detection using next generation sequencing on the inbred chicken lines with the assistance of CNVnator. As a result, a total of 1649 CNV regions (CNVRs) were successfully identified after merging all the nine datasets, of which 90 CNVRs were overlapped across all the chicken lines. Within these shared regions, 1360 harbored genes were identified. In addition, 55 and 44 CNVRs with 62 and 57 harbored genes were specifically identified in line 6₃ and 7₂, respectively. Bioinformatics analysis showed that the nearby genes were significantly enriched in 36 GO terms and 6 KEGG pathways including JAK/STAT signaling pathway. Ten CNVRs (nine deletions and one duplication) involved in 10 disease-related genes were selected for validation by using quantitative real-time PCR (qPCR), all of which were successfully confirmed. Finally, qPCR was also used to validate two deletion events in line 7₂ that were definitely normal in line 6₃. One high-confidence gene, IRF2 was identified as the most promising candidate gene underlying resistance and susceptibility to MD in view of its function and overlaps with data from previous study.

Conclusions

Our findings provide valuable insights for understanding the genetic mechanism of resistance to MD and the identified gene and pathway could be considered as the subject of further functional characterization.

First genome-wide CNV mapping in FELIS CATUS using next generation sequencing data

Background

Copy Number Variations (CNVs) have becoming very significant variants, representing a major source of genomic variation. CNVs involvement in phenotypic expression and different diseases has been widely demonstrated in humans as well as in many domestic animals. However, genome wide investigation on these structural variations is still missing in Felis catus. The present work is the first CNV mapping from a large data set of Next Generation Sequencing (NGS) data in the domestic cat, performed within the 99 Lives Consortium.

Results

Reads have been mapped on the reference assembly_6.2 by Maverix Biomics. CNV detection with cn.MOPS and CNVnator detected 592 CNVs. These CNVs were used to obtain 154 CNV Regions (CNVRs) with BedTools, including 62 singletons. CNVRs covered 0.26% of the total cat genome with 129 losses, 19 gains and 6 complexes. Cluster Analysis and Principal Component Analysis of the detected CNVRs showed that breeds tend to cluster together as well as cats sharing the same geographical origins. The 46 genes identified within the CNVRs were annotated.

Conclusion

This study has improved the genomic characterization of 14 cat breeds and has provided CNVs information that can be used for studies of traits in cats. It can be considered a sound starting point for genomic CNVs identification in this species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5297-2) contains supplementary material, which is available to authorized users.

A copy number variation generated by complicated organization of PCDHA gene cluster is associated with egg performance traits in Xinhua E-strain

In recent years, a mass of duplicated and deleted DNA sequences have been found in human and animal genomes following the prevalence of employing high-throughput sequencing and SNP array. However, few copy number variation (CNV) studies have been performed on egg performance traits of chicken. In this study, 17 loci reported in previous studies were selected for CNV detection in the Xinhua E-strain by using the CNVplex kit, and the detection results showed that locus14 exhibited CNV. Further association analysis indicated the copies of locus14 could be significantly associated with age at first egg (AFE; P < 0.0086) and egg number at 250 d (250EN; P < 0.036). DNA sequence amplification showed the loss of a 260-bp-long fragment in the upstream of locus14, which mainly occurred in normal or copy-gain individuals. The qPCR results showed that subjects with gain of copies could promote the total expression level of the PCDHA gene cluster in the pituitary gland of adult individuals. Additionally, PCR amplification with randomly combined primers revealed a larger number of chicken variable exons than that previously reported, indicating the complexity of the organization of the PCDHA gene cluster. Those variable exons are divergent in their distribution among the populations of Xinhua E-strain, Chahua, Tibetan, and Tulufan Game Chicken, and most individuals only possess part of variable exons. Overall, the copies of locus14 reflect the variable exon dosage effects on the total expression level of the PCDHA gene cluster, which may regulate the layer egg production by affecting the development of the neural system.

Detection of breed-specific copy number variations in domestic chicken genome

Copy number variations (CNVs) are important large-scale variants. They are widespread in the genome and may contribute to phenotypic variation. Detection and characterization of CNVs can provide new insights into the genetic basis of important traits. Here, we perform whole-genome short read sequence analysis to identify CNVs in two indigenous and commercial chicken breeds to evaluate the impact of the identified CNVs on breed-specific traits. After filtration, a total of 12 955 CNVs spanning (on average) about 9.42% of the chicken genome were found that made up 5467 CNV regions (CNVRs). Chicken quantitative trait loci (QTL) datasets and Ensembl gene annotations were used as resources for the estimation of potential phenotypic effects of our CNVRs on breed-specific traits. In total, 34% of our detected CNVRs were also detected in earlier CNV studies. These CNVRs partly overlap several previously reported QTL and gene ontology terms associated with some important traits, including shank length QTL in Creeper-specific CNVRs and body weight and egg production characteristics, as well as muscle and body organ growth, in the Arian commercial breed. Our findings provide new insights into the genomic structure of the chicken genome for an improved understanding of the potential roles of CNVRs in differentiating between breeds or lines.

Genomic variability in Mexican chicken population using copy number variants

BACKGROUND

Copy number variations are genome polymorphism that influence phenotypic variation and are an important source of genetic variation in populations. The aim of this study was to investigate genetic variability in the Mexican Creole chicken population using CNVs.

RESULTS

The Hidden Markov Model of the PennCNV software detected a total of 1924 CNVs in the genome of the 256 samples processed with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix). The mapped CNVs comprised 1538 gains and 386 losses, resulting at population level in 1216 CNV regions (CNVRs), of which 959 gains, 226 losses and 31 complex (i.e. containing both losses and gains). The CNVRs covered a total of 47 Mb of the whole genome sequence length, corresponding to 5.12% of the chicken galGal4 autosome assembly.

CONCLUSIONS

This study allowed a deep insight into the structural variation in the genome of unselected Mexican chicken population, which up to now has not been genetically characterized. The genomic study disclosed that the population, even if presenting extreme morphological variation, cannot be organized in differentiated genetic subpopulations. Finally this study provides a chicken CNV map based on the 600 K SNP chip array jointly with a genome-wide gene copy number estimates in a native unselected for more than 500 years chicken population.

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Genomic and genetic variability of six chicken populations using single nucleotide polymorphism and copy number variants as markers

Genomic and genetic variation among six Italian chicken native breeds (Livornese, Mericanel della Brianza, Milanino, Bionda Piemontese, Bianca di Saluzzo and Siciliana) were studied using single nucleotide polymorphism (SNP) and copy number variants (CNV) as markers. A total of 94 DNA samples genotyped with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix) were used in the analyses. The results showed the genetic and genomic variability occurring among the six Italian chicken breeds. The genetic relationship among animals was established with a principal component analysis. The genetic diversity within breeds was calculated using heterozygosity values (expected and observed) and with Wright's F-statistics. The individual-based CNV calling, based on log R ratio and B-allele frequency values, was done by the Hidden-Markov Model (HMM) of PennCNV software on autosomes. A hierarchical agglomerative clustering was applied in each population according to the absence or presence of definite CNV regions (CNV were grouped by overlapping of at least 1 bp). The CNV map was built on a total of 1003 CNV found in individual samples, after grouping by overlaps, resulting in 564 unique CNV regions (344 gains, 213 losses and 7 complex), for a total of 9.43 Mb of sequence and 1.03% of the chicken assembly autosome. All the approaches using SNP data showed that the Siciliana breed clearly differentiate from other populations, the Livornese breed separates into two distinct groups according to the feather colour (i.e. white and black) and the Bionda Piemontese and Bianca di Saluzzo breeds are closely related. The genetic variability found using SNP is comparable with that found by other authors in the same breeds using microsatellite markers. The CNV markers analysis clearly confirmed the SNP results.

A Complex Structural Variation on Chromosome 27 Leads to the Ectopic Expression of HOXB8 and the Muffs and Beard Phenotype in Chickens

Muffs and beard (Mb) is a phenotype in chickens where groups of elongated feathers gather from both sides of the face (muffs) and below the beak (beard). It is an autosomal, incomplete dominant phenotype encoded by the Muffs and beard (Mb) locus. Here we use genome-wide association (GWA) analysis, linkage analysis, Identity-by-Descent (IBD) mapping, array-CGH, genome re-sequencing and expression analysis to show that the Mb allele causing the Mb phenotype is a derived allele where a complex structural variation (SV) on GGA27 leads to an altered expression of the gene HOXB8. This Mb allele was shown to be completely associated with the Mb phenotype in nine other independent Mb chicken breeds. The Mb allele differs from the wild-type mb allele by three duplications, one in tandem and two that are translocated to that of the tandem repeat around 1.70 Mb on GGA27. The duplications contain total seven annotated genes and their expression was tested during distinct stages of Mb morphogenesis. A continuous high ectopic expression of HOXB8 was found in the facial skin of Mb chickens, strongly suggesting that HOXB8 directs this regional feather-development. In conclusion, our results provide an interesting example of how genomic structural rearrangements alter the regulation of genes leading to novel phenotypes. Further, it again illustrates the value of utilizing derived phenotypes in domestic animals to dissect the genetic basis of developmental traits, herein providing novel insights into the likely role of HOXB8 in feather development and differentiation.

Genetic variation is a key part for the study of evolution, development and differentiation. In domestic animals, many breeds display striking phenotypes that differentiate them from their wild ancestors. Several of these have been related to structural variations, including Fibromelanosis and Rose-comb in chickens, Double-muscled and Osteopetrosis in cattle, Cone degeneration in dogs, and White coat color in pigs. The feather is a type of skin appendages that exists in multiple variants on different body parts, and the derived feathering phenotypes in domestic birds are perfect resources to decipher the mechanisms regulating feather development and differentiation. Here we study the genetics of the Muffs and beard trait, a variant that alters the feather development in the facial area of chickens. We show that this phenotype is associated with a genomic structural variant that leads to an ectopic expression of HOXB8 in the facial skin during feather development. This is thus another example of how structural variants in the genome lead to novel, derived phenotypic changes in domestic animals and suggests an important role for HOXB8 in feather development.

Copy number variation identification and analysis of the chicken genome using a 60K SNP BeadChip

Copy number variation (CNV) is an important source of genetic variation in organisms and a main factor that affects phenotypic variation. A comprehensive study of chicken CNV can provide valuable information on genetic diversity and facilitate future analyses of associations between CNV and economically important traits in chickens. In the present study, an F2 full-sib chicken population (554 individuals), established from a cross between Xinghua and White Recessive Rock chickens, was used to explore CNV in the chicken genome. Genotyping was performed using a chicken 60K SNP BeadChip. A total of 1,875 CNV were detected with the PennCNV algorithm, and the average number of CNV was 3.42 per individual. The CNV were distributed across 383 independent CNV regions (CNVR) and covered 41 megabases (3.97%) of the chicken genome. Seven CNVR in 108 individuals were validated by quantitative real-time PCR, and 81 of these individuals (75%) also were detected with the PennCNV algorithm. In total, 274 CNVR (71.54%) identified in the current study were previously reported. Of these, 147 (38.38%) were reported in at least 2 studies. Additionally, 109 of the CNVR (28.46%) discovered here are novel. A total of 709 genes within or overlapping with the CNVR was retrieved. Out of the 2,742 quantitative trait loci (QTL) collected in the chicken QTL database, 43 QTL had confidence intervals overlapping with the CNVR, and 32 CNVR encompassed one or more functional genes. The functional genes located in the CNVR are likely to be the QTG that are associated with underlying economic traits. This study considerably expands our insight into the structural variation in the genome of chickens and provides an important resource for genomic variation, especially for genomic structural variation related to economic traits in chickens.

Identification of Copy Number Variations in Xiang and Kele Pigs

Xiang and Kele pigs are two well-known local Chinese pig breeds that possess rich genetic resources and have enormous economic and scientific value. We performed a comprehensive genomic analysis of the copy number variations (CNVs) in these breeds. CNVs are one of the most important forms of genomic variation and have profound effects on phenotypic variation. In this study, PorcineSNP60 genotyping data from 98 Xiang pigs and 22 Kele pigs were used to identify CNVs. In total, 172 candidate CNV regions (CNVRs) were identified, ranging from 3.19 kb to 8175.26 kb and covering 80.41 Mb of the pig genome. Approximately 56.40% (97/172) of the CNVRs overlapped with those identified in seven previous studies, and 43.60% (75/172) of the identified CNVRs were novel. Of the identified CNVRs, 82 (47 gain, 33 loss, and two gain-loss events that covered 4.58 Mb of the pig genome) were found only in a Xiang population with a large litter size. In contrast, 13 CNVRs (8 gain and 5 loss events) were unique to a Xiang population with small litter sizes, and 30 CNVRs (14 loss and 16 gain events) were unique to Kele pigs. The CNVRs span approximately 660 annotated Sus scrofa genes that are significantly enriched for specific biological functions, such as sensory perception, cognition, reproduction, ATP biosynthetic processes, and neurological processes. Many CNVR-associated genes, particularly the genes involved in reproductive traits, differed between the Xiang populations with large and small litter sizes, and these genes warrant further investigation due to their importance in determining the reproductive performance of Xiang pigs. Our results provide meaningful information about genomic variation, which may be useful in future assessments of the associations between CNVs and important phenotypes in Xiang and Kele pigs to ultimately help protect these rare breeds.

Improved Detection and Characterization of Copy Number Variations Among Diverse Pig Breeds by Array CGH

As a major component of genomic variation, copy number variations (CNVs) are considered as promising markers for some phenotypic and economically important traits in domestic animals. Using a custom-designed 1M array CGH (aCGH), we performed CNV discovery in 12 pig samples from one Asian wild boar population, six Chinese indigenous breeds, and two European commercial breeds. In total, we identified 758 CNV regions (CNVRs), covering 47.43 Mb of the pig genome sequence. Of the total porcine genes, 1295 genes were completely or partially overlapped with the identified CNVRs, which enriched in the terms related to sensory perception of the environment, neurodevelopmental processes, response to external stimuli, and immunity. Further probing the potential functions of these genes, we also found a suite of genes related important traits, which make them a promising resource for exploring the genetic basis of phenotype differences among diverse pig breeds. Compared with previous relevant studies, the current study highlights that different platforms can complement each other, and the combined implementation of different platforms is beneficial to achieve the most comprehensive CNV calls. CNVs detected in diverse populations herein are essentially complementary to the CNV map in the pig genome, which would be helpful for understanding the pig genome variants and investigating the associations between various phenotypes and CNVs.

Quantitative effect of a CNV on a morphological trait in chickens

Copy Number Variation has been associated with morphological traits, developmental defects or disease susceptibility. The autosomal dominant Pea-comb mutation in chickens is due to the massive amplification of a CNV in intron 1 of SOX5 and provides a unique opportunity to assess the effect of variation in the number of repeats on quantitative traits such as comb size and comb mass in Pea-comb chickens. The quantitative variation of comb size was estimated by 2D morphometry and the number of repeats (RQ) was estimated by qPCR, in a total of 178 chickens from 3 experimental lines, two of them showing segregation for the Pea-comb mutation. This study included only Pea-comb chickens. Analysis of variance showed highly significant effects of line and sex on comb measurements. Adult body weight (BW) and RQ were handled as covariates. BW significantly influenced comb mass but not comb size. RQ values significantly influenced comb size, and the linear regression coefficient was highest for heterozygous carriers: the higher the number of repeats, the smaller the comb size. A similar trend was observed for comb mass. The CNV contributed to 3.4% of the phenotypic variance of comb size in heterozygous carriers of the CNV, an order of magnitude frequently encountered for QTLs. Surprisingly, there was no such relationship between RQ values and comb size in the homozygous line. It may be concluded that heterozygosity for a CNV in a non-coding region may contribute to phenotypic plasticity.

Genome-wide copy number profiling using high-density SNP array in chickens

Phenotypic diversity is a direct consequence resulting mainly from the impact of underlying genetic variation, and recent studies have shown that copy number variation (CNV) is emerging as an important contributor to both phenotypic variability and disease susceptibility. Herein, we performed a genome-wide CNV scan in 96 chickens from 12 diversified breeds, benefiting from the high-density Affymetrix 600 K SNP arrays. We identified a total of 231 autosomal CNV regions (CNVRs) encompassing 5.41 Mb of the chicken genome and corresponding to 0.59% of the autosomal sequence. The length of these CNVRs ranged from 2.6 to 586.2 kb with an average of 23.4 kb, including 130 gain, 93 loss and eight both gain and loss events. These CNVRs, especially deletions, had lower GC content and were located particularly in gene deserts. In particular, 102 CNVRs harbored 128 chicken genes, most of which were enriched in immune responses. We obtained 221 autosomal CNVRs after converting probe coordinates to Galgal3, and comparative analysis with previous studies illustrated that 153 of these CNVRs were regarded as novel events. Furthermore, qPCR assays were designed for 11 novel CNVRs, and eight (72.73%) were validated successfully. In this study, we demonstrated that the high-density 600 K SNP array can capture CNVs with higher efficiency and accuracy and highlighted the necessity of integrating multiple technologies and algorithms. Our findings provide a pioneering exploration of chicken CNVs based on a high-density SNP array, which contributes to a more comprehensive understanding of genetic variation in the chicken genome and is beneficial to unearthing potential CNVs underlying important traits of chickens.

Genome-wide patterns of copy number variation in the diversified chicken genomes using next-generation sequencing

BACKGROUND

Copy number variation (CNV) is important and widespread in the genome, and is a major cause of disease and phenotypic diversity. Herein, we performed a genome-wide CNV analysis in 12 diversified chicken genomes based on whole genome sequencing.

RESULTS

A total of 8,840 CNV regions (CNVRs) covering 98.2 Mb and representing 9.4% of the chicken genome were identified, ranging in size from 1.1 to 268.8 kb with an average of 11.1 kb. Sequencing-based predictions were confirmed at a high validation rate by two independent approaches, including array comparative genomic hybridization (aCGH) and quantitative PCR (qPCR). The Pearson's correlation coefficients between sequencing and aCGH results ranged from 0.435 to 0.755, and qPCR experiments revealed a positive validation rate of 91.71% and a false negative rate of 22.43%. In total, 2,214 (25.0%) predicted CNVRs span 2,216 (36.4%) RefSeq genes associated with specific biological functions. Besides two previously reported copy number variable genes EDN3 and PRLR, we also found some promising genes with potential in phenotypic variation. Two genes, FZD6 and LIMS1, related to disease susceptibility/resistance are covered by CNVRs. The highly duplicated SOCS2 may lead to higher bone mineral density. Entire or partial duplication of some genes like POPDC3 may have great economic importance in poultry breeding.

CONCLUSIONS

Our results based on extensive genetic diversity provide a more refined chicken CNV map and genome-wide gene copy number estimates, and warrant future CNV association studies for important traits in chickens.

Identification and functional characterization of copy number variations in diverse chicken breeds

Background

The detection and functional characterization of genomic structural variations are important for understanding the landscape of genetic variation in the chicken. A recently recognized aspect of genomic structural variation, called copy number variation (CNV), is gaining interest in chicken genomic studies. The aim of the present study was to investigate the pattern and functional characterization of CNVs in five characteristic chicken breeds, which will be important for future studies associating phenotype with chicken genome architecture.

Results

Using a commercial 385 K array-based comparative genomic hybridization (aCGH) genome array, we performed CNV discovery using 10 chicken samples from four local Chinese breeds and the French breed Houdan chicken. The female Anka broiler was used as a reference. A total of 281 copy number variation regions (CNVR) were identified, covering 12.8 Mb of polymorphic sequences or 1.07% of the entire chicken genome. The functional annotation of CNVRs indicated that these regions completely or partially overlapped with 231 genes and 1032 quantitative traits loci, suggesting these CNVs have important functions and might be promising resources for exploring differences among various breeds. In addition, we employed quantitative PCR (qPCR) to further validate several copy number variable genes, such as prolactin receptor, endothelin 3 (EDN3), suppressor of cytokine signaling 2, CD8a molecule, with important functions, and the results suggested that EDN3 might be a molecular marker for the selection of dark skin color in poultry production. Moreover, we also identified a new CNVR (chr24: 3484617–3512275), encoding the sortilin-related receptor gene, with copy number changes in only black-bone chicken.

Conclusions

Here, we report a genome-wide analysis of the CNVs in five chicken breeds using aCGH. The association between EDN3 and melanoblast proliferation was further confirmed using qPCR. These results provide additional information for understanding genomic variation and related phenotypic characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-934) contains supplementary material, which is available to authorized users.

A genome-wide detection of copy number variation using SNP genotyping arrays in Beijing-You chickens

Copy number variation (CNV) has been recently examined in many species and is recognized as being a source of genetic variability, especially for disease-related phenotypes. In this study, the PennCNV software, a genome-wide CNV detection system based on the 60 K SNP BeadChip was used on a total sample size of 1,310 Beijing-You chickens (a Chinese local breed). After quality control, 137 high confidence CNVRs covering 27.31 Mb of the chicken genome and corresponding to 2.61 % of the whole chicken genome. Within these regions, 131 known genes or coding sequences were involved. Q-PCR was applied to verify some of the genes related to disease development. Results showed that copy number of genes such as, phosphatidylinositol-5-phosphate 4-kinase II alpha, PHD finger protein 14, RHACD8 (a CD8α- like messenger RNA), MHC B-G, zinc finger protein, sarcosine dehydrogenase and ficolin 2 varied between individual chickens, which also supports the reliability of chip-detection of the CNVs. As one source of genomic variation, CNVs may provide new insight into the relationship between the genome and phenotypic characteristics.

A cis-regulatory mutation of PDSS2 causes silky-feather in chickens

Silky-feather has been selected and fixed in some breeds due to its unique appearance. This phenotype is caused by a single recessive gene (hookless, h). Here we map the silky-feather locus to chromosome 3 by linkage analysis and subsequently fine-map it to an 18.9 kb interval using the identical by descent (IBD) method. Further analysis reveals that a C to G transversion located upstream of the prenyl (decaprenyl) diphosphate synthase, subunit 2 (PDSS2) gene is causing silky-feather. All silky-feather birds are homozygous for the G allele. The silky-feather mutation significantly decreases the expression of PDSS2 during feather development in vivo. Consistent with the regulatory effect, the C to G transversion is shown to remarkably reduce PDSS2 promoter activity in vitro. We report a new example of feather structure variation associated with a spontaneous mutation and provide new insight into the PDSS2 function.

The feather is an excellent model for evolution and development due to its complex structure and vast diversity. Some chickens have silky-feather because of a loss of hooklets in pennaceous feathers, while most chickens have the wild-type normal feather. Hooklets are formed in the last differentiation stage of the life cycle of a pennaceous feather. Chickens with silky-feather are homozygous for a recessive allele (hookless, h). Silkie chicken from China is one of the breeds showing the fascinating silky-feather phenotype and the breed has been known for hundreds of years. In this study, we mapped the silky-feather locus to an 18.9 kb interval and identified a single nucleotide polymorphism (SNP) completely associated with silky-feather. The causative mutation is located 103 base pairs upstream of the coding sequence of prenyl (decaprenyl) diphosphate synthase, subunit 2 (PDSS2). The expression of the PDSS2 gene is decreased in silky-feather skin during feather development in vivo. The silky-feather allele also reduces the PDSS2 promoter activity in vitro. This is the first report of feather structure variation associated with PDSS2 and provides new insight into molecular signaling in the late development stage of feather morphogenesis.