An initial map of chromosomal segmental copy number variations in the chicken

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

The prolactin receptor: A cross-species comparison of gene structure, transcriptional regulation, tissue-specificity, and genetic variation

The conserved and multifaceted functions of prolactin (PRL) are coordinated through varied distribution and expression of its cell-surface receptor (PRLR) across a range of tissues and physiological states. The resultant heterogeneous expression of PRLR mRNA and protein across different organs and cell types supports a wide range of PRL-regulated processes including reproduction, lactation, development, and homeostasis. Genetic variation within the PRLR gene also accounts for several phenotypes impacting agricultural production and human pathology. The goal of this review is to highlight the many elements that control differential expression of the PRLR across tissues, and the various phenotypes that exist across species due to variation in the PRLR gene.

Impact of STAT5A-CNVs on growth traits in Chinese beef cattle breeds

CNVs, which are a type of structural variation, make a substantial impact on diverse characteristics in multiple species. Q-PCR and data association analysis were used for STAT5A gene copy in this study. This study aimed to investigate the copy number variation (CNV) of the STAT5A gene in seven Chinese cattle breeds, namely Qinchuan cattle, Xianan cattle, Yunling cattle, Ji'an cattle, Jiaxian Red cattle, Qaidam cattle, and Guyuan yellow cattle. Blood samples were collected for CNV typing, and the correlation between CNV type and growth traits was analyzed using SPSS 23.0 software and ANOVA. The findings revealed variations in the distribution of different copy number types among the different cattle breeds. Furthermore, association analysis demonstrated a positive impact of CNV in the STAT5A gene on cattle growth: in the JX, individuals with duplication types exhibited superior performance in terms of rump length (P < 0.05). Conversely, normal GY cattle demonstrated better body height and abdomen circumference (P < 0.05), while QD cattle exhibited a significant correlation between weight and body length with normal individuals (P < 0.05). Moreover, QC bovine duplication individuals outperformed other types, with copy number variation significantly associated with chest depth, chest width, and body length (P < 0.05). The results validate the correlation between copy number variation (CNV) of the STAT5A gene and growth characteristics in five different cattle breeds, providing a reliable benchmark for the purpose of cattle breeding.

Whole-genome sequencing and genomic analysis of Norduz goat (Capra hircus)

Artificial and natural selective breeding of goats has resulted in many different goat breeds all around the world. Norduz goat is one of these breeds, and it is a local goat breed of Turkey. The goats are favorable due to pre-weaning viability and reproduction values compared to the regional breeds. Development in sequencing technologies has let to understand huge genomic structures and complex phenotypes. Until now, such a comprehensive study has not been carried out to understand the genomic structure of the Norduz goats, yet. In the study, the next-generation sequencing was carried out to understand the genomic structure of Norduz goat. Real-time PCR was used to evaluate prominent CNVs in the Norduz goat individuals. Whole genome of the goat was constructed with an average of 33.1X coverage level. In the stringent filtering condition, 9,757,980 SNPs, 1,536,715 InDels, and 290 CNVs were detected in the Norduz goat genome. Functional analysis of high-impact SNP variations showed that the classical complement activation biological process was affected significantly in the goat. CNVs in the goat genome were found in genes related to defense against viruses, immune response, and cell membrane transporters. It was shown that GBP2, GBP5, and mammalian ortholog GBP1, which are INF-stimulated GTPases, were found to be high copy numbers in the goats. To conclude, genetic variations mainly in immunological response processes suggest that Norduz goat is an immunologically improved goat breed and natural selection could take an important role in the genetical improvements of the goats.

Effects of Genetic Variation of the Sorting Nexin 29 (SNX29) Gene on Growth Traits of Xiangdong Black Goat

Previous studies have found that the copy number variation (CNV) and insertion/deletion (indels) located in the sorting nexin 29 (SNX29) gene, which is an important candidate gene related to meat production and quality, are associated with growth traits of African goats and Shaanbei white cashmere goats. However, the genetic effects of SNX29 genetic variation on growth traits of Xiangdong black (XDB) goat (a representative meat goat breed in China) are still unclear. The purpose of this study was to detect the mRNA expression level of SNX29 and to explore the genetic effects of CNV and indel within SNX29 on growth traits and gene expression in XDB goat. The SNX29 mRNA expression profile showed that the SNX29 was highly expressed in adipose tissues, indicating that the SNX29 gene could play a key role in subcutaneous adipose deposition of XDB goat. 17 bp indel (g.10559298-10559314), 21 bp indel (g.10918982-10919002) and CNV were detected in 516 individuals of XDB goat by PCR or qPCR. The association analysis of SNX29 CNV with growth traits in XDB goats showed that SNX29 CNV was significantly correlated with chest circumference and abdominal circumference (p < 0.01), and the normal type of SNX29 CNV goat individuals were more advantageous. For the mRNA expression of SNX29 gene, individuals with SNX29 copy number normal type had a higher trend than that of SNX29 gene with copy number gain type in longissimus dorsi muscle (p = 0.07), whereas individuals with SNX29 copy number gain type had a higher trend in abdominal adipose (p = 0.09). Overall, these results suggested that the SNX29 gene could play an important role in growth and development of XDB goats and could be used for marker-assisted selection (MAS) in XDB goats.

Genome-wide evaluation of copy gain and loss variations in three Afghan sheep breeds

Copy number variation (CNV) is one of the main sources of variation between different individuals that has recently attracted much researcher interest as a major source for heritable variation in complex traits. The aim of this study was to identify CNVs in Afghan indigenous sheep consisting of three Arab, Baluchi, and Gadik breeds using genomic arrays containing 53,862 single nucleotide polymorphism (SNP) markers. Data were analyzed using the Hidden Markov Model (HMM) of PennCNV software. In this study, out of 45 sheep studied, 97.8% (44 animals) have shown CNVs. In total, 411 CNVs were observed for autosomal chromosomes and the entire sequence length of around 144 Mb was identified across the genome. The average number of CNVs per each sheep was 9.13. The identified CNVs for Arab, Baluchi, and Gadik breeds were 306, 62, and 43, respectively. After merging overlapped regions, a total of 376 copy number variation regions (CNVR) were identified, which are 286, 50, and 40 for Arab, Baluchi, and Gadik breeds, respectively. Bioinformatics analysis was performed to identify the genes and QTLs reported in these regions and the biochemical pathways involved by these genes. The results showed that many of these CNVRs overlapped with the genes or QTLs that are associated with various pathways such as immune system development, growth, reproduction, and environmental adaptions. Furthermore, to determine a genome-wide pattern of selection signatures in Afghan sheep breeds, the unbiased estimates of FST was calculated and the results indicated that 37 of the 376 CNVRs (~ 10%) have been also under selection signature, most of those overlapped with the genes influencing production, reproduction and immune system. Finally, the statistical methods used in this study was applied in an external dataset including 96 individuals of the Iranian sheep breed. The results indicated that 20 of the 114 CNVRs (18%) identified in Iranian sheep breed were also identified in our study, most of those overlapped with the genes influencing production, reproduction and immune system. Overall, this is the first attempts to develop the genomic map of loss and gain variation in the genome of Afghan indigenous sheep breeds, and may be important to shed some light on the genomic regions associated with some economically important traits in these breeds.

Towards the detection of copy number variation from single sperm sequencing in cattle

BACKGROUND

Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms.

RESULTS

We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull's family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats.

CONCLUSION

Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.

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.

Genome-wide detection of CNVs and their association with performance traits in broilers

Background

Copy number variations (CNVs) are a major type of structural genomic variants that underlie genetic architecture and phenotypic variation of complex traits, not only in humans, but also in livestock animals. We identified CNVs along the chicken genome and analyzed their association with performance traits. Genome-wide CNVs were inferred from Affymetrix® high density SNP-chip data for a broiler population. CNVs were concatenated into segments and association analyses were performed with linear mixed models considering a genomic relationship matrix, for birth weight, body weight at 21, 35, 41 and 42 days, feed intake from 35 to 41 days, feed conversion ratio from 35 to 41 days and, body weight gain from 35 to 41 days of age.

Results

We identified 23,214 autosomal CNVs, merged into 5042 distinct CNV regions (CNVRs), covering 12.84% of the chicken autosomal genome. One significant CNV segment was associated with BWG on GGA3 (q-value = 0.00443); one significant CNV segment was associated with BW35 (q-value = 0.00571), BW41 (q-value = 0.00180) and BW42 (q-value = 0.00130) on GGA3, and one significant CNV segment was associated with BW on GGA5 (q-value = 0.00432). All significant CNV segments were verified by qPCR, and a validation rate of 92.59% was observed. These CNV segments are located nearby genes, such as KCNJ11, MyoD1 and SOX6, known to underlie growth and development. Moreover, gene-set analyses revealed terms linked with muscle physiology, cellular processes regulation and potassium channels.

Conclusions

Overall, this CNV-based GWAS study unravels potential candidate genes that may regulate performance traits in chickens. Our findings provide a foundation for future functional studies on the role of specific genes in regulating performance in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07676-1.

Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

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.

Copy number variation of CADM2 gene revealed its association with growth traits across Chinese Capra hircus (goat) populations

Copy number variations (CNVs) are the wide structural variations ranging from 50 bp to several Mb at genome which can affect gene expression and further impacting growth and development traits of livestock. Comparing with single nucleotide polymorphisms (SNPs), CNVs can better explain the genetic and phenotypic diversity, are increasingly important in biological research. As a member of immunoglobulin super-family, cell adhesion molecule 2 (CADM2) plays a vital role in cancer development and metabolic regulation. Here, we tested the CNV of CADM2 gene in 443 goats across five breeds (Guizhou white goat, GZW; Guizhou black goat, GZB; Africa Nubian goat, AN; Boer goat × Huai goat, BH; Boer goat, BG) and detected its association with phenotypic traits. Subsequently, we analyzed the CADM2 gene expression level in different tissues of NB goats (n = 3, Nubian × Black) and the transcriptional expression in lung is much higher than others. The results showed that the CNV of CADM2 has a significant association with withers height and body length in GZB goat (P < 0.01), in which individuals with type of deletion were superior to those with duplication or normal type in term of body hight and body length (P < 0.01). In summary, this study confirmed the association between CNV of CADM2 gene and growth traits, and our research data indicated the CADM2-CNV may considered as a prospective candidate for the molecular marker-assisted selection breeding of goat growth traits, which conducived to accelerating the genetic amelioration in Chinese goats.

The Association of the Copy Number Variation of the MLLT10 Gene with Growth Traits of Chinese Cattle

Simple Summary

Copy number variation is a common genetic polymorphism, mainly represented by submicroscopic levels of deletion and duplication, caused by rearrangement of the genome. It is well-known that the copy number variation of a gene is associated with growth traits of livestock. In this study, we detected the correlation between the copy number variation of the the MLLT10 gene and the growth traits of Chinese yellow cattle. We found that the copy number variation of the MLLT10 gene has a significant influence on hip width, rump length, hucklebone width, and cannon bone circumference of some Chinese yellow cattle breeds. The results provide preliminary suggestions for Chinese yellow cattle breeding and new insights about the future of copy number variation (CNV) as a new promising molecular marker in animal breeding.

Abstract

Copy number variation is a part of genomic structural variation and has caused widespread concern. According to the results of high-throughput screening of the MLLT10 gene, we found that the copy number variation region of the MLLT10 gene was correlated with bovine growth traits. We aimed to detect the MLLT10 gene copy number variation and provide materials for the Chinese yellow cattle breed. In this study, the SPSS software was used to analyze the correlation among the copy number type of six different cattle breeds (i.e., Qinchuan, Xianan, Jiaxian, Yanbian, Sinan, Yunling) and the corresponding growth traits. The results showed the following: In Qinchuan cattle, the copy number duplication type was greater than the deletion and normal types; in Xianan cattle, the copy number duplication and normal types were less as compared with the deletion type; and in Yunling cattle, the frequency of the duplication type was dominant among the three types of copy number variants. The correlation analysis result showed that there is a significant correlation between the copy number variation (CNV) of the MLLT10 gene and the growth traits of three cattle breeds. Furthermore, correlation analysis showed that MLLT10 CNV had positive effects on growth traits such as hip width, rump length, hucklebone width, and cannon bone circumference (p < 0.05). This study provides a basis for the molecular-assisted marker breeding of cattle and contributes to the breeding of cattle.

Linc-GALMD1 Regulates Viral Gene Expression in the Chicken

A rapidly increasing number of reports on dysregulated long intergenic non-coding RNA (lincRNA) expression across numerous types of cancers indicates that aberrant lincRNA expression may be a major contributor to tumorigenesis. Marek’s disease (MD) is a T cell lymphoma of chickens induced by Marek’s disease virus (MDV). Although we have investigated the roles of lincRNAs in bursa tissue of MDV-infected chickens in previous studies, the molecular mechanisms of lincRNA functions in T cells remain poorly understood. In the present study, Linc-GALMD1 was identified from CD4+ T cells and MSB1 cells, and its expression was significantly downregulated in MD-resistant line of birds in response to MDV challenge. Furthermore, loss-of-function experiments indicated that linc-GALMD1 significantly affected the expression of 290 genes in trans. Through integrated analysis of differentially expressed genes (DEGs) induced by MDV and linc-GALMD1, we found that IGLL1 gene expression levels had a positive correlation with the degree of MD infection and could potentially serve as an indicator for clinical diagnosis of MD. Moreover, an interaction between MDV and linc-GALMD1 was also observed. Accordingly, chicken embryonic fibroblast cells were inoculated with MDV with and without the linc-GALMD1 knockdown, and the data showed that linc-GALMD1 could repress MDV gene expression during the course of MDV infection. These findings uncovered a role of linc-GALMD1 as a viral gene regulator and suggested a function of linc-GALMD1 contributing to tumor suppression by coordinating expression of MDV genes and tumor-related genes and regulating immune responses to MDV infection.

Associations of ORMDL1 gene copy number variations with growth traits in four Chinese sheep breeds

Copy number variations (CNVs) are gains and losses of genomic sequence of more than 50 bp between two individuals of a species. Also, CNV is considered to be one of the main elements affecting the phenotypic diversity and evolutionary adaptation of animals. ORMDL sphingolipid biosynthesis regulator 1 (ORMDL1) is a protein-coding gene associated with diseases and development. In our study, the polymorphism of ORMDL1 gene copy numbers in four Chinese sheep breeds (abbreviated CK, HU, STH, and LTH) was detected. In addition, we analyzed the transcriptional expression level of ORMDL1 gene in different tissues of sheep and examined the association of ORMDL1 CNV with growth traits. The statistical analysis revealed that ORMDL1 CNV was remarkably correlated with body height, heart girth, and circumference of cannon bone in HU sheep (P<0.05), and there are significant effects on body weight, body height, body length, chest depth, and height of hip cross in STH sheep (P<0.05). In conclusion, our results provide a basis for the relationship between CNV of ORMDL1 gene and sheep growth traits, suggesting that ORMDL1 CNV may be considered a promising marker for the molecular breeding of Chinese sheep.

Identification of Copy Number Variation in Domestic Chicken Using Whole-Genome Sequencing Reveals Evidence of Selection in the Genome

Simple Summary

Chickens have been bred for meat and egg production as a source of animal protein. With the increase of productivity as the main purpose of domestication, factors such as metabolism and immunity were boosted, which are detectable signs of selection on the genome. This study focused on copy number variation (CNV) to find evidence of domestication on the genome. CNV was detected from whole-genome sequencing of 65 chickens including Red Jungle Fowl, broilers, and layers. After that, CNV region, the overlapping region of CNV between individuals, was made to identify which genomic regions showed copy number differentiation. The 663 domesticated-specific CNV regions were associated with various functions such as metabolism and organ development. Also, by performing population differentiation analyses such as clustering analysis and ANOVA test, we found that there are a lot of genomic regions with different copy number patterns between broilers and layers. This result indicates that different genetic variations can be found, depending on the purpose of artificial selection and provides considerations for future animal breeding.

Abstract

Copy number variation (CNV) has great significance both functionally and evolutionally. Various CNV studies are in progress to find the cause of human disease and to understand the population structure of livestock. Recent advances in next-generation sequencing (NGS) technology have made CNV detection more reliable and accurate at whole-genome level. However, there is a lack of CNV studies on chickens using NGS. Therefore, we obtained whole-genome sequencing data of 65 chickens including Red Jungle Fowl, Cornish (broiler), Rhode Island Red (hybrid), and White Leghorn (layer) from the public databases for CNV region (CNVR) detection. Using CNVnator, a read-depth based software, a total of 663 domesticated-specific CNVRs were identified across autosomes. Gene ontology analysis of genes annotated in CNVRs showed that mainly enriched terms involved in organ development, metabolism, and immune regulation. Population analysis revealed that CN and RIR are closer to each other than WL, and many genes (LOC772271, OR52R1, RD3, ADH6, TLR2B, PRSS2, TPK1, POPDC3, etc.) with different copy numbers between breeds found. In conclusion, this study has helped to understand the genetic characteristics of domestic chickens at CNV level, which may provide useful information for the development of breeding systems in chickens.

Copy Number Variation of the SHE Gene in Sheep and Its Association with Economic Traits

Simple Summary

Src Homology 2 Domain Containing E (SHE) is a protein coding gene, and member of the SH2 domain-containing family. Sequencing revealed a 2000 bp copy number variation in the SHE gene. There is overlap between SHE copy number variation (CNV) and quantitative trait loci related to milk fat percentage and bone density. A total of 750 sheep, including Chaka sheep (CKS), Hu sheep (HS), Small Tail Han sheep (STHS) and Large Tail Han sheep (LTHS) were available to detect the CNV of the SHE gene and correlate these gene variations with economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS. Association analysis showed that gain and normal copy number types performed better in body length (p < 0.05), circumference of cannon bone (p < 0.05), heart girth (p < 0.05), chest width (p < 0.05) and high at the cross (p < 0.05) in CKS, HS and STHS. Chi-square analyses found significant variation in the CNV of the SHE gene, so it varies greatly between varieties. These findings clarified the relationship between the CNV of the SHE gene and the economic traits in these four kinds of sheep, and provide a reference for sheep breeding.

Abstract

Copy number variation (CNV) caused by gene rearrangement is an important part of genomic structural variation. We found that the copy number variation region of the Src Homology 2 Domain Containing E (SHE) gene correlates with a quantitative trait locus of sheep related to milk fat percentage and bone density. The aim of our study was to detect the copy number variation of the SHE gene in four sheep breeds and to conduct a correlation analysis with economic traits, hoping to provide some reference for sheep breeding. In this study, we examined 750 sheep from four Chinese breeds: Chaka sheep (CKS), Hu sheep (HS), Large Tail Han sheep (LTHS) and Small Tail Han sheep (STHS). We used qPCR to evaluate the copy number of the SHE gene, and then used general linear models to analyze the associations between CNV and economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS individuals. Association analyses showed that gain and normal copy number types were correlated to body length, circumference of cannon bone, heart girth, chest width and high at the cross in CKS, HS and STHS (p < 0.05), but this association was not observed for LTHS. Chi-square values (χ²) found prominent differences in CNV distribution among the studied breeds. Overall, the CNV of the SHE gene may be an important consideration for sheep molecular breeding.

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.

Integrating CNVs into meta-QTL identified GBP4 as positional candidate for adult cattle stature

Copy number variation (CNV) of DNA sequences, functionally significant but yet fully ascertained, is believed to confer considerable increments in unexplained heritability of quantitative traits. Identification of phenotype-associated CNVs (paCNVs) therefore is a pressing need in CNV studies to speed up their exploitation in cattle breeding programs. Here, we provided a new avenue to achieve this goal that is to project the published CNV data onto meta-quantitative trait loci (meta-QTL) map which connects causal genes with phenotypes. Any CNVs overlapping meta-QTL therefore will be potential paCNVs. This study reported potential paCNVs in Bos taurus autosome 3 (BTA3). Notably, overview indexes and CNVs both highlighted a narrower region (BTA3 54,500,000-55,000,000 bp, named BTA3_INQTL_6) within one constructed meta-QTL. Then, we ascertained guanylate-binding protein 4 (GBP4) among the nine positional candidate genes was significantly associated with adult cattle stature, including body weight (BW, P < 0.05) and withers height (WHT, P < 0.05), fitting GBP4 CNV either with three levels or with six levels in the model. Although higher copy number downregulated the mRNA levels of GBP2 (P < 0.05) and GBP4 (P < 0.05) in 1-Mb window (54.0-55.0 Mb) in muscle and adipose, additional analyses will be needed to clarify the causality behind the ascertained association.

Genome-wide detection of CNVs associated with beak deformity in chickens using high-density 600K SNP arrays

Beak deformity (crossed beaks) is found in several indigenous chicken breeds including Beijing-You studied here. Birds with deformed beaks have reduced feed intake and poor production performance. Recently, copy number variation (CNV) has been examined in many species and is recognized as a source of genetic variation, especially for disease phenotypes. In this study, to unravel the genetic mechanisms underlying beak deformity, we performed genome-wide CNV detection using Affymetrix chicken high-density 600K data on 48 deformed-beak and 48 normal birds using penncnv. As a result, two and eight CNV regions (CNVRs) covering 0.32 and 2.45 Mb respectively on autosomes were identified in deformed-beak and normal birds respectively. Further RT-qPCR studies validated nine of the 10 CNVRs. The ratios of six CNVRs were significantly different between deformed-beak and normal birds (P < 0.01). Within these six regions, three and 21 known genes were identified in deformed-beak and normal birds respectively. Bioinformatics analysis showed that these genes were enriched in six GO terms and one KEGG pathway. Five candidate genes in the CNVRs were further validated using RT-qPCR. The expression of LRIG2 (leucine rich repeats and immunoglobulin like domains 2) was lower in birds with deformed beaks (P < 0.01). Therefore, the LRIG2 gene could be considered a key factor in view of its known functions and its potential roles in beak deformity. Overall, our results will be helpful for future investigations of the genomic structural variations underlying beak deformity in chickens.

Associations of GBP2 gene copy number variations with growth traits and transcriptional expression in Chinese cattle

Copy number variations (CNVs) recently have been recognized as another important genetic variability followed single nucleotide polymorphisms (SNPs). The guanylate binding protein 2 (GBP2) gene plays an important role in cell proliferation. This study was performed to determine the presence of GBP2 CNV (relative to Angus cattle) in 466 individuals representing six main cattle breeds from China, identify its relationship with growth, and explore the biological effects of gene expression. There were two CNV regions in the GBP2 gene, for three types, CNV1 loss type (relative to Angus cattle) was more frequent in XN than other breeds, and CNV2 loss type (relative to Angus cattle) was more frequent in XN and CDM than other breeds. Though the GBP2 gene copy number presented no correlation with the transcriptional expression of JX (P > .05), but the transcriptional expression in heart is higher than other tissues, and the copy number in muscles and fat of JX is higher than others breeds. Statistical analysis revealed that the GBP2 gene CNV1 and CNV2 were significantly associated with growth traits (P < .05). In conclusion, this research established the correlations between CNVs of GBP2 gene and growth traits in different cattle breeds, and our results suggested that the CNVs in GBP2 gene may be considered markers for the molecular breeding of Chinese beef cattle.

Genome-wide analysis reveals differential selection involved with copy number variation in diverse Chinese Cattle

Copy number variations (CNVs) are defined as deletions, insertions, and duplications between two individuals of a species. To investigate the diversity and population-genetic properties of CNVs and their diverse selection patterns, we performed a genome-wide CNV analysis using high density SNP array in Chinese native cattle. In this study, we detected a total of 13,225 CNV events and 3,356 CNV regions (CNVRs), overlapping with 1,522 annotated genes. Among them, approximately 71.43 Mb of novel CNVRs were detected in the Chinese cattle population for the first time, representing the unique genomic resources in cattle. A new V _i statistic was proposed to estimate the region-specific divergence in CNVR for each group based on unbiased estimates of pairwise V _ST . We obtained 12 and 62 candidate CNVRs at the top 1% and top 5% of genome-wide V _i value thresholds for each of four groups (North, Northwest, Southwest and South). Moreover, we identified many lineage-differentiated CNV genes across four groups, which were associated with several important molecular functions and biological processes, including metabolic process, response to stimulus, immune system, and others. Our findings provide some insights into understanding lineage-differentiated CNVs under divergent selection in the Chinese native cattle.

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species, white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.

Long intergenic non-coding RNA GALMD3 in chicken Marek's disease

Long intergenic non-coding RNAs (lincRNAs) are transcribed from non-coding DNA sequences. Studies have revealed that aberrant expressions of lincRNAs are associated with various types of cancers and neurological disorders. Marek's disease (MD) is a highly contagious T-cell lymphoid neoplasia of chicken induced by Marek's disease virus (MDV). In this study, we first identified and validated linc-GALMD3 highly expressed in MDV-infected CD4+ T cells by RNA-Seq and qRT-PCR. By RNA-Seq analysis in MDCC-MSB1 cells after loss of function of linc-GALMD3 by shRNA, we found that linc-GALMD3 could positively cis-regulate its downstream gga-miR-223 gene expression. In contrast, it could trans-regulate the 748 differentially expressed genes (FDR < 0.01) that were mainly enriched into mitochondrial structure and cell cycle processes using GO analysis. Of these, the most significantly expressed gene EPYC might cause iris lesion in MD. The other eight genes, NDUFA4, NDUFB6, NDUFV1, NDUFS8, SDHB, UQCRC1, UQCRC2, and COX7A2, actively participated in oxidative phosphorylation in mitochondrial dysfunction and cell death. Most importantly, we found that the MDV replication was repressed when linc-GALMD3 was knocked down in CEF cells. Our results suggested that linc-GALMD3 might be a critical regulator in chicken MD and could be used as a candidate-promising mark for MD prevention, diagnosis, and treatment.

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.

Association study and expression analysis of CYP4A11 gene copy number variation in Chinese cattle

The identification of copy number variations (CNVs) allow us to explore genomic polymorphisms. In recent years, significant progress in understanding CNVs has been made in studies of human and animals, however, association and expression studies of CNVs are still in the early stage. It was previously reported that the Cytochrome P-450 4A11 (CYP4A11) gene is located within a copy number variable region (CNVR) that encompasses quantitative trait loci (QTLs) for economic traits like meat quality and milk production. So, this study was performed to determine the presence of CYP4A11 CNV in six distinct cattle breeds, identify its relationship with growth, and explore the biological effects of gene expression. For three CYP4A11 CNV types, Normal was more frequent than Gain or Loss. Association analysis revealed a positive effect of CYP4A11 copy number on growth traits (P < 0.05). One-way analysis of variance (ANOVA) analysis revealed that more CYP4A11 copies increased the gene expression level. Moreover, overexpression of CYP4A11 in vitro revealed its effect on lipid deposit. The data provide evidence for the functional role of CYP4A11 CNV and provide the basis for future applications in cattle breeding.

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.

Detection and validation of structural variations in bovine whole-genome sequence data

BACKGROUND

Several examples of structural variation (SV) affecting phenotypic traits have been reported in cattle. Currently the identification of SV from whole-genome sequence data (WGS) suffers from a high false positive rate. Our aim was to construct a high quality set of SV calls in cattle using WGS data. First, we tested two SV detection programs, Breakdancer and Pindel, and the overlap of these methods, on simulated sequence data to determine their precision and sensitivity. We then identified population SV from WGS of 252 Holstein and 64 Jersey bulls based on the overlapping calls from the two programs. In addition, we validated an overlapped SV set in 28 twice-sequenced Holstein individuals, and in another two validated sets (one for each breed) that were transmitted from sire to son. We also tested whether highly conserved gene sets across eukaryotes and recently expanded gene families in bovine were depleted and enriched, respectively, for SV.

RESULTS

In empirical WGS data, 17,518 SV covering 27.36 Mb were found in the Holstein population and 4285 SV covering 8.74 Mb in the Jersey population, of which 4.62 Mb of SV overlapped between Holsteins and Jerseys. A total of 11,534 candidate SV covering 5.64 Mb were validated in the 28 twice-sequenced individuals, while 3.49 and 0.67 Mb of SV were validated from Holstein and Jersey sire-son transmission, respectively. Only eight of 237 core eukaryotic genes had at least a 50-bp overlap with an SV from our validated sets, suggesting that conserved genes are depleted for SV (p < 0.05). In addition, we observed that recently expanded gene families were significantly more associated with SV than other genes. Long interspersed nuclear elements-1 were enriched for deletions when compared to the rest of the genome (p = 0.0035).

CONCLUSIONS

We reported SV from 252 Holstein and 64 Jersey individuals. A considerable proportion of Jersey population SV (53.5%) were also found in Holstein. In contrast, about 76.90% sire-son transmission validated SV were present in Jerseys and Holsteins. The enrichment of SV in expanding gene families suggests that SV can be a source of genetic variation for evolution.

Duplication of chicken defensin7 gene generated by gene conversion and homologous recombination

Defensins constitute an evolutionary conserved family of cationic antimicrobial peptides that play a key role in host innate immune responses to infection. Defensin genes generally reside in complex genomic regions that are prone to structural variation, and defensin genes exhibit extensive copy number variation in humans and in other species. Copy number variation of defensin genes was examined in inbred lines of Leghorn and Fayoumi chickens, and a duplication of defensin7 was discovered in the Fayoumi breed. Analysis of junction sequences confirmed the occurrence of a simple tandem duplication of defensin7 with sequence identity at the junction, suggesting nonallelic homologous recombination between defensin7 and defensin6 The duplication event generated two chimeric promoters that are best explained by gene conversion followed by homologous recombination. Expression of defensin7 was not elevated in animals with two genes despite both genes being transcribed in the tissues examined. Computational prediction of promoter regions revealed the presence of several putative transcription factor binding sites generated by the duplication event. These data provide insight into the evolution and possible function of large gene families and specifically, the defensins.

Deletion of Indian hedgehog gene causes dominant semi-lethal Creeper trait in chicken

The Creeper trait, a classical monogenic phenotype of chicken, is controlled by a dominant semi-lethal gene. This trait has been widely cited in the genetics and molecular biology textbooks for illustrating autosomal dominant semi-lethal inheritance over decades. However, the genetic basis of the Creeper trait remains unknown. Here we have utilized ultra-deep sequencing and extensive analysis for targeting causative mutation controlling the Creeper trait. Our results indicated that the deletion of Indian hedgehog (IHH) gene was only found in the whole-genome sequencing data of lethal embryos and Creeper chickens. Large scale segregation analysis demonstrated that the deletion of IHH was fully linked with early embryonic death and the Creeper trait. Expression analysis showed a much lower expression of IHH in Creeper than wild-type chickens. We therefore suggest the deletion of IHH to be the causative mutation for the Creeper trait in chicken. Our findings unravel the genetic basis of the longstanding Creeper phenotype mystery in chicken as the same gene also underlies bone dysplasia in human and mouse, and thus highlight the significance of IHH in animal development and human haploinsufficiency disorders.

Comparative analyses across cattle genders and breeds reveal the pitfalls caused by false positive and lineage-differential copy number variations

We compared CNV region (CNVR) results derived from 1,682 Nellore cattle with equivalent results derived from our previous analysis of Bovine HapMap samples. By comparing CNV segment frequencies between different genders and groups, we identified 9 frequent, false positive CNVRs with a total length of 0.8 Mbp that were likely caused by assembly errors. Although there was a paucity of lineage specific events, we did find one 54 kb deletion on chr5 significantly enriched in Nellore cattle. A few highly frequent CNVRs present in both datasets were detected within genomic regions containing olfactory receptor, ATP-binding cassette, and major histocompatibility complex genes. We further evaluated their impacts on downstream bioinformatics and CNV association analyses. Our results revealed pitfalls caused by false positive and lineage-differential copy number variations and will increase the accuracy of future CNV studies in both taurine and indicine cattle.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits

Background

Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance.

Results

Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively.

Conclusions

Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2752-9) contains supplementary material, which is available to authorized users.

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.

Genome-wide identification of copy number variations between two chicken lines that differ in genetic resistance to Marek's disease

Background

Copy number variation (CNV) is a major source of genome polymorphism that directly contributes to phenotypic variation such as resistance to infectious diseases. Lines 6₃ and 7₂ are two highly inbred experimental chicken lines that differ greatly in susceptibility to Marek’s disease (MD), and have been used extensively in efforts to identify the genetic and molecular basis for genetic resistance to MD. Using next generation sequencing, we present a genome-wide assessment of CNVs that are potentially associated with genetic resistance to MD.

Methods

Three chickens randomly selected from each line were sequenced to an average depth of 20×. Two popular software, CNVnator and Pindel, were used to call genomic CNVs separately. The results were combined to obtain a union set of genomic CNVs in the two chicken lines.

Results

A total of 5,680 CNV regions (CNVRs) were identified after merging the two datasets, of which 1,546 and 1,866 were specific to the MD resistant or susceptible line, respectively. Over half of the line-specific CNVRs were shared by 2 or more chickens, reflecting the reduced diversity in both inbred lines. The CNVRs fixed in the susceptible lines were significantly enriched in genes involved in MAPK signaling pathway. We also found 67 CNVRs overlapping with 62 genes previously shown to be strong candidates of the underlying genes responsible for the susceptibility to MD.

Conclusions

Our findings provide new insights into the genetic architecture of the two chicken lines and additional evidence that MAPK signaling pathway may play an important role in host response to MD virus infection. The rich source of line-specific CNVs is valuable for future disease-related association studies in the two chicken lines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2080-5) contains supplementary material, which is available to authorized users.

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.

Uncovering genomic features and maternal origin of korean native chicken by whole genome sequencing

The Korean Native Chicken (KNC) is an important endemic biological resource in Korea. While numerous studies have been conducted exploring this breed, none have used next-generation sequencing to identify its specific genomic features. We sequenced five strains of KNC and identified 10.9 million SNVs and 1.3 million InDels. Through the analysis, we found that the highly variable region common to all 5 strains had genes like PCHD15, CISD1, PIK3C2A, and NUCB2 that might be related to the phenotypic traits of the chicken such as auditory sense, growth rate and egg traits. In addition, we assembled unaligned reads that could not be mapped to the reference genome. By assembling the unaligned reads, we were able to present genomic sequences characteristic to the KNC. Based on this, we also identified genes related to the olfactory receptors and antigen that are common to all 5 strains. Finally, through the reconstructed mitochondrial genome sequences, we performed phylogenomic analysis and elucidated the maternal origin of the artificially restored KNC. Our results revealed that the KNC has multiple maternal origins which are in agreement with Korea's history of chicken breed imports. The results presented here provide a valuable basis for future research on genomic features of KNC and further understanding of KNC's origin.

Global patterns of apparent copy number variation in birds revealed by cross-species comparative genomic hybridization

There is a growing interest in copy number variation (CNV) and the recognition of its importance in phenotype, disease, adaptation and speciation. CNV data is usually ascertained by array-CGH within-species, but similar inter-species comparisons have also been made in primates, mice and domestic mammals. Here, we conducted a broad appraisal of putative cross-species CNVs in birds, 16 species in all, using the standard array-CGH approach. Using a chicken oligonucleotide microarray, we detected 790 apparent CNVs within 135 unique regions and developed a bioinformatic tool 'CNV Analyser' for analysing and visualising cross-species data sets. We successfully addressed four hypotheses as follows: (a) Cross-species CNVs (compared to chicken) are, as suggested from preliminary evidence, smaller and fewer in number than in mammals; this 'dogma' was rejected in the light of the new evidence. (b) CNVs in birds are likely to have a functional effect through an association with genes; a large proportion of detected regions (70 %) were indeed associated with genes (suggesting functional significance), however, not necessarily more so than in mammals. (c) There are more CNVs in birds with more rearranged karyotypes; this hypothesis was rejected. Indeed, Falco species contained fewer than most with relatively standard (chicken-like) karyotypes. (d) There are more CNVs per megabase on micro-chromosomes than macrochromosomes; this hypothesis was accepted. Indeed, in species with rearranged karyotypes characterised by chromosomal fusions, the fused former microchromosomes still 'behaved' as though they were their microchromosomal ancestors. Gene ontology analysis of CNVRs revealed enrichment in immune response and antigen presentation genes and five CNVRs were perfectly correlated with the unique loss of sexual dichromatism in one Galliformes species.

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.

Copy number variation in the horse genome

We constructed a 400K WG tiling oligoarray for the horse and applied it for the discovery of copy number variations (CNVs) in 38 normal horses of 16 diverse breeds, and the Przewalski horse. Probes on the array represented 18,763 autosomal and X-linked genes, and intergenic, sub-telomeric and chrY sequences. We identified 258 CNV regions (CNVRs) across all autosomes, chrX and chrUn, but not in chrY. CNVs comprised 1.3% of the horse genome with chr12 being most enriched. American Miniature horses had the highest and American Quarter Horses the lowest number of CNVs in relation to Thoroughbred reference. The Przewalski horse was similar to native ponies and draft breeds. The majority of CNVRs involved genes, while 20% were located in intergenic regions. Similar to previous studies in horses and other mammals, molecular functions of CNV-associated genes were predominantly in sensory perception, immunity and reproduction. The findings were integrated with previous studies to generate a composite genome-wide dataset of 1476 CNVRs. Of these, 301 CNVRs were shared between studies, while 1174 were novel and require further validation. Integrated data revealed that to date, 41 out of over 400 breeds of the domestic horse have been analyzed for CNVs, of which 11 new breeds were added in this study. Finally, the composite CNV dataset was applied in a pilot study for the discovery of CNVs in 6 horses with XY disorders of sexual development. A homozygous deletion involving AKR1C gene cluster in chr29 in two affected horses was considered possibly causative because of the known role of AKR1C genes in testicular androgen synthesis and sexual development. While the findings improve and integrate the knowledge of CNVs in horses, they also show that for effective discovery of variants of biomedical importance, more breeds and individuals need to be analyzed using comparable methodological approaches.

Genomes of individuals in a species vary in many ways, one of which is DNA copy number variation (CNV). This includes deletions, duplications, and complex rearrangements typically larger than 50 base-pairs. CNVs are part of normal genetic variation contributing to phenotypic diversity but can also be pathogenic and associated with diseases and disorders. In order to distinguish between the two, detailed knowledge about CNVs in the species of interest is needed. Here we studied the genomes of 38 normal horses of 16 diverse breeds, and identified 258 CNV regions. We integrated our findings with previously published horse CNVs and generated a composite dataset of ∼1400 CNVRs. Despite this large number, our analysis shows that CNV research in horses needs further improvement because the current data are based on 10% of horse breeds and that most CNVRs are study-specific and require validation. Finally, we analyzed CNVs in horses with disorders of sexual development and found in two male pseudo-hermaphrodites a large deletion disrupting a group of genes involved in sex hormone metabolism and sexual differentiation. The findings underline the possible role of CNVs in complex disorders such as development and reproduction.

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.