Copy number variants in Italian Large White pigs detected using high‐density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance‐tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with penncnv software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by penncnv were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.     

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.     

Identification of copy number variations in three Chinese horse breeds using 70K single nucleotide polymorphism BeadChip array

Copy number variation (CNV), an essential form of genetic variation, has been increasingly recognized as one promising genetic marker in the analysis of animal genomes. Here, we used the Equine 70K single nucleotide polymorphism genotyping array for the genome‐wide detection of CNVs in 96 horses from three diverse Chinese breeds: Debao pony (DB), Mongolian horse (MG) and Yili horse (YL). A total of 287 CNVs were determined and merged into 122 CNV regions (CNVRs) ranging from 199 bp to 2344 kb in size and distributed in a heterogeneous manner on chromosomes. These CNVRs were integrated with seven existing reports to generate a composite genome‐wide dataset of 1558 equine CNVRs, revealing 69 (56.6%) novel CNVRs. The majority (69.7%) of the 122 CNVRs overlapped with 438 genes, whereas 30.3% were located in intergenic regions. Most of these genes were associated with common CNVRs, which were shared by divergent horse breeds. As many as 60, 42 and 91 genes overlapping with the breed‐specific ss were identified in DB, MG and YL respectively. Among these genes, FGF11, SPEM1, PPARG, CIDEB, HIVEP1 and GALR may have potential relevance to breed‐specific traits. These findings provide valuable information for understanding the equine genome and facilitating association studies of economically important traits with equine CNVRs in the future.     

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves)

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves—an asthmalike inflammatory disease in horses. Analysis of 16 RAO‐susceptible (cases) and six RAO‐resistant (control) horses on a custom‐made whole‐genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO‐affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.     

Copy number variation detection using SNP genotyping arrays in three Chinese pig breeds

We performed genome‐wide CNV detection based on SNP genotyping data of 96 Chinese‐native Tibetan, Dahe and Wuzhishan pigs. These pigs are particularly interesting because of their excellent adaptation to hypoxia or small body size, which facilitates the use of them as models of different human diseases in addition to valuable agricultural animals. A total of 105 CNV regions (CNVRs) were identified, encompassing 16.71 Mb of the pig genome. Seven of 10 (70%) CNVRs selected randomly were validated by quantitative real‐time PCR. Comparison with previous studies revealed 25 (23.81%) novel CNVRs, indicating that CNV coverage of the pig genome is still incomplete and there exists large diversity between pig breeds. Functional analysis of genes located in these CNVRs confirmed the high representation of genes involved in sensory perception, neurological system processes and other basic metabolic processes. In addition, the majority of these CNVRs were detected to span reported pig QTL that affect various traits, which highlighted three biologically interesting genes with copy number changes (i.e., ANKRD34B, FAM110B and ABCG1). These genes may have economic importance in pig breeding and are worth being further investigated. We also obtained some CNVRs harboring genes that had human orthologs involved in human diseases such as cardiovascular disease and Alzheimer's disease. The findings of this study are a significant extension of the coverage of CNVRs in the pig genome and provide valuable resources for follow‐up‐associated studies of CNVs in pig complex traits as well as important implications of human diseases.     

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.

Electronic supplementary material

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

A genome‐wide scan for copy number variations using high‐density single nucleotide polymorphism array in Simmental cattle

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome‐wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high‐density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.     

Copy number variants in a highly inbred Iberian porcine strain

We carried out a comprehensive genomic analysis of porcine copy number variants (CNVs) based on whole‐genome SNP genotyping data and provided new measures of genomic diversity (number, length and distribution of CNV events) for a highly inbred strain (the Guadyerbas strain). This strain represents one of the most ancient surviving populations of the Iberian breed, and it is currently in serious danger of extinction. CNV detection was conducted on the complete Guadyerbas population, adjusted for genomic waves, and used strict quality criteria, pedigree information and the latest porcine genome annotation. The analysis led to the detection of 65 CNV regions (CNVRs). These regions cover 0.33% of the autosomal genome of this particular strain. Twenty‐nine of these CNVRs were identified here for the first time. The relatively low number of detected CNVRs is in line with the low variability and high inbreeding estimated previously for this Iberian strain using pedigree, microsatellite or SNP data. A comparison across different porcine studies has revealed that more than half of these regions overlap with previously identified CNVRs or multicopy regions. Also, a preliminary analysis of CNV detection using whole‐genome sequence data for four Guadyerbas pigs showed overlapping for 16 of the CNVRs, supporting their reliability. Some of the identified CNVRs contain relevant functional genes (e.g., the SCD and USP15 genes), which are worth being further investigated because of their importance in determining the quality of Iberian pig products. The CNVR data generated could be useful for improving the porcine genome annotation.     

A genome‐wide SNP scan reveals two loci associated with the chicken resistance to Marek's disease

Marek's disease (MD) is a neoplastic disease in chickens, caused by the Marek's disease virus (MDV). To investigate host genetic resistance to MD, we conducted a genome‐wide association study (GWAS) on 67 MDV‐infected chickens based on a case and control design, including 57 susceptible chickens in the case group and 10 resistant chickens as controls. After searching 38 655 valid genomic markers, two SNPs were found to be associated with host resistance to MD. One SNP, rs14527240, reaching chromosome‐wide significance level (P < 0.01) was located in the SPARC‐related modular calcium‐binding 1 (SMOC1) gene on GGA5. The other one, GGaluGA156129, reaching genome‐wide significance (P < 0.05), was located in the protein tyrosine phosphatase, non‐receptor type 3 (PTPN3) gene on GGA2. In addition, expression patterns of these two genes in spleens were detected by qPCR. The expression of SMOC1 was significantly up‐regulated (P < 0.05), whereas the expression of PTNP3 did not show significance when the case group was compared with the control group. Up‐regulation of SMOC1 in susceptible spleens suggests its important roles in MD tumorigenesis. This is the first study to investigate MD‐resistant loci, and it demonstrates the power of GWASs for mapping genes associated with MD resistance.     

Microanatomy of passerine hard‐cornified tissues: Beak and claw structure of the black‐capped chickadee (Poecile atricapillus)

The microanatomy of healthy beaks and claws in passerine birds has not been well described in the literature, despite the importance of these structures in avian life. Histological processing of hard‐cornified tissues is notoriously challenging and only a few reports on effective techniques have been published. An emerging epizootic of beak deformities among wild birds in Alaska and the Pacific Northwest region of North America recently highlighted the need for additional baseline information about avian hard‐cornified structures. In this study, we examine the beak and claw of the Black‐capped Chickadee (Poecile atricapillus), a common North American passerine that is affected by what has been described as “avian keratin disorder.” We use light and scanning electron microscopy and high‐magnification radiography to document the healthy microanatomy of these tissues and identify features of functional importance. We also describe detailed methods for histological processing of avian hard‐cornified structures and discuss the utility of special stains. Results from this study will assist in future research on the functional anatomy and pathology of hard‐cornified structures and will provide a necessary reference for ongoing investigations of avian keratin disorder in Black‐capped Chickadees and other wild passerine species. © 2011 Wiley Periodicals, Inc.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Copy number variation associates with mortality in long‐lived individuals: a genome‐wide assessment

Copy number variants (CNVs) represent a significant source of genetic variation in the human genome and have been implicated in numerous diseases and complex traits. To date, only a few studies have investigated the role of CNVs in human lifespan. To investigate the impact of CNVs on prospective mortality at the extreme end of life, where the genetic component of lifespan appears most profound, we analyzed genomewide CNV data in 603 Danish nonagenarians and centenarians (mean age 96.9 years, range 90.0–102.5 years). Replication was performed in 500 long‐lived individuals from the Leiden Longevity Study (mean age 93.2 years, range 88.9–103.4 years). First, we assessed the association between the CNV burden of each individual (the number of CNVs, the average CNV length, and the total CNV length) and mortality and found a significant increase in mortality per 10 kb increase in the average CNV length, both for all CNVs (hazard ratio (HR) = 1.024, P = 0.002) and for duplications (HR = 1.011, P = 0.005), as well as per 100 kb increase in the total length of deletions (HR = 1.009, P = 0.0005). Next, we assessed the relation between specific deletions and duplications and mortality. Although no genome–wide significant associations were discovered, we identified six deletions and one duplication that showed consistent association with mortality in both or either of the sexes across both study populations. These results indicate that the genome–wide CNV burden, specifically the average CNV length and the total CNV length, associates with higher mortality in long‐lived individuals.     

Genome‐wide analysis of TNF‐alpha response in pigs challenged with porcine circovirus 2b

Tumor necrosis factor alpha (TNF‐α) is a pro‐inflammatory cytokine with a role in activating adaptive immunity to viral infections. By inhibiting the capacity of plasmacytoid dendritic cells to produce interferon‐α and TNF‐α, porcine circovirus 2 (PCV2) limits the maturation of myeloid dendritic cells and impairs their ability to recognize viral and bacterial antigens. Previously, we reported QTL for viremia and immune response in PCV2‐infected pigs. In this study, we analyzed phenotypic and genetic relationships between TNF‐α protein levels, a potential indicator of predisposition to PCV2 co‐infection, and PCV2 susceptibility. Following experimental challenge with PCV2b, TNF‐α reached the peak at 21 days post‐infection (dpi), at which time a difference was observed between pigs that expressed extreme variation in viremia and growth (P < 0.10). A genome‐wide association study (n = 297) revealed that genotypes of 56 433 SNPs explained 73.9% of the variation in TNF‐α at 21 dpi. Major SNPs were identified on SSC8, SSC10 and SSC14. Haplotypes based on SNPs from a SSC8 (9 Mb) 1‐Mb window were associated with variation in TNF‐α (P < 0.02), IgG (P = 0.05) and IgM (P < 0.13) levels at 21 dpi. Potential overlap of regulatory mechanisms was supported by the correlations between genomic prediction values of TNF‐α and PCV2 antibodies (21 dpi, r > 0.22), viremia (14–21 dpi, P > 0.29) and viral load (r = 0.31, P < 0.0001). Characterization of the QTL regions uncovered genes that could influence variation in TNF‐α levels as well as T‐ and B‐cell development, which can affect disease susceptibility.     

Real‐time PCR detection of Enterococcus faecalis associated with amyloid arthropathy

Aims: To develop a RT‐PCR method for detection of the multilocus sequence type 82 of Enterococcus faecalis associated with amyloid arthropathy (AA) in layers. Methods and Results: Bacteria were selected from lesions including AA in layers. The primers were designed based on the phosphate ATP binding cassette transporter (pstS) and xanthine phosphoribosyltransferase (xpt) genes and first tested against three isolates with known base pairs at the specific sites. Subsequently, 12 isolates were selected from our collection by one researcher, and RT‐PCR was performed blinded. The sequence type (ST) was then confirmed by multilocus sequence analysis. Two single‐nucleotide polymorphisms in the pstS and xpt genes allowed an unambiguous identification of ST82. As an alternative to DNA extraction, a boiling method for release of DNA from cells was used. Conclusions: The real‐time PCR targeting ST82 enables rapid screening of Ent. faecalis cultured from suspect cases with results available after a few hours, much faster than multilocus sequence typing and pulse field gel electrophoresis. Significance and Impact of the Study: The new method allows a rapid screening of isolates with results available after only few hours. This RT‐PCR method could be a useful tool for molecular epidemiological studies on the spread of arthropathic and amyloidogenic Ent. faecalis within and between birds more efficiently.     

The non‐coding RNA OTUB1‐isoform2 promotes ovarian tumour progression and predicts poor prognosis

Ovarian cancer is the leading malignancy of the female reproductive system and is associated with inconspicuous early invasion and metastasis. We have previously reported that the oncogene OTUB1 plays a crucial role in ovarian cancer progression, but the role of its isoform, the non‐coding RNA OTUB1‐isoform2, in ovarian cancer is still elusive. Here, we reported that OTUB1‐isoform2 expression in ovarian cancer tissues was significantly higher than that in the paired paratumorous tissues (P < .01). The patients with high expression of OTUB1‐isoform2 had larger tumours than those with low expression (P < .05). The high expression of OTUB1‐isoform2 was correlated with the involvement of bilateral ovaries (P < .05), lymph node metastasis (P < .05), vascular invasion (P < .05), greater omentum involvement (P < .01), fallopian tube involvement (P < .05), advanced FIGO stages (P < .01) and recurrence (P < .01). Moreover, OTUB1‐isoform2 served as an independent negative prognostic predictor for disease‐free survival (DFS) and disease‐specific survival (DSS). Overexpression of OTUB1‐isoform2 in the ovarian cancer cells stimulated cell proliferation, migration and invasion both in vitro and in vivo. In summary, our study suggested that OTUB1‐isoform2 is a novel prognostic biomarker with independent oncogenic functions for ovarian cancer.     

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30–35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.