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.     

Detection of copy number variants in the horse genome and examination of their association with recurrent laryngeal neuropathy

We used the data from a recently performed genome‐wide association study using the Illumina Equine SNP50 beadchip for the detection of copy number variants (CNVs) and examined their association with recurrent laryngeal neuropathy (RLN), an important equine upper airway disease compromising performance. A total of 2797 CNVs were detected for 477 horses, covering 229 kb and seven SNPs on average. Overlapping CNVs were merged to define 478 CNV regions (CNVRs). CNVRs, particularly deletions, were shown to be significantly depleted in genes. Fifty‐two of the 67 common CNVRs (frequency ≥ 1%) were validated by association mapping, Mendelian inheritance, and/or Mendelian inconsistencies. None of the 67 common CNVRs were significantly associated with RLN when accounting for multiple testing. However, a duplication on chromosome 10 was detected in 10 cases (representing three breeds) and two unphenotyped parents but in none of the controls. The duplication was embedded in an 8‐Mb haplotype shared across breeds.     

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.     

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.     

Occurrence and Genetic Characteristics of Cryptosporidium hominis and Cryptosporidium andersoni in Horses from Southwestern China

A total of 333 fecal specimens from horses in southwestern China were genotyped based on analysis of the small subunit rRNA (SSU rRNA) gene. Cryptosporidium hominis and Cryptosporidium andersoni were identified in 2 and 4 stool specimens, respectively. The identification of C. hominis was confirmed by sequence analysis of the 70‐kDa heat shock protein (HSP70) and oocyst wall protein (COWP) genes. Subtyping analysis of the 60‐kDa glycoprotein (GP60) gene sequence of C. hominis revealed a new rare subtype Id, named IdA15; only three Id isolates have been reported in humans to date. Multilocus sequence typing (MLST) analysis indicated that the C. andersoni subtype was A6, A5, A2, and A1 at the four minisatellite loci (MS1, MS2, MS3, and MS16, respectively). This is the first report to identify the presence of C. andersoni and C. hominis in horses in southwestern China and the first to identify a rare zoonotic subtype Id of C. hominis in horses. These findings suggest that infected horses may act as potential reservoirs of Cryptosporidium to transmit infections to humans.     

A genome‐wide association study for equine recurrent airway obstruction in European Warmblood horses reveals a suggestive new quantitative trait locus on chromosome 13

Recurrent airway obstruction (RAO), also known as heaves, is an asthma‐like respiratory disease. Its development is strongly influenced by environmental risk factors such as sensitization and exposure to moldy hay, straw bedding and stabling indoors. A hereditary component has been documented in previous studies; however, so far no causative genetic variant that influences the risk of developing RAO has been identified. In this study, we revised an existing dataset and selected 384 horses for genotyping on the Affymetrix high‐density equine SNP array. We performed an allelic case–control genome‐wide association study, which revealed a suggestively significant association on equine chromosome 13 at 32 843 309 bp. This SNP is located in the protein‐coding gene TXNDC11, which is possibly involved in the folding process of the multiprotein complexes DUOX1 and DUOX2. In humans, these proteins are known to take part in regulating the production of H₂O₂ in the respiratory tract epithelium as well as in MUC5AC mucin expression. Therefore, TXNDC11 may be considered a functional candidate gene, and further research is needed to explore its potential role in RAO‐affected horses.     

Genome‐wide association study of osteochondrosis in the tarsocrural joint of Dutch Warmblood horses identifies susceptibility loci on chromosomes 3 and 10

Equine osteochondrosis is a developmental joint disease that is a significant source of morbidity affecting multiple breeds of horse. The genetic variants underlying osteochondrosis susceptibility have not been established. Here, we describe the results of a genome‐wide association study of osteochondrosis using 90 cases and 111 controls from a population of Dutch Warmblood horses. We report putative associations between osteochondrosis and loci on chromosome 3 (BIEC2‐808543; P = 5.03 × 10^?7) and chromosome 10 (BIEC2‐121323; P = 2.62 × 10^?7).     

Copy number variation and variant discovery in Bullmastiff dogs

Identification of genomic variants within dogs is important for understanding genetic factors contributing to breed diversity and phenotypic traits. This study aimed to identify sources of variation in the Bullmastiff using high‐density signal intensity and whole‐genome sequence data. Close to 3000 copy number variants (CNVs) were identified in Bullmastiff dogs using Canine HD BeadChip data. When CNVs were collated, 82 CNV regions (CNVRs) were detected, 50% in transcribed regions encompassing 432 genes. Fifty of the CNVRs detected have not been reported in other breeds and represent potential breed‐specific variants. A proportion of the CNVR variants with predicted modifying effects on gene pathways may contribute to breed traits. Approximately 5 million putative variants per dog, inclusive of single nucleotide polymorphisms (SNPs), multi‐nucleotide polymorphisms (MNPs) and insertion and deletions (INDELs), were identified from DNA sequence data on a small number of animals. Identification of genetic variants in the Bullmastiff highlights sources of variation in the breed and molecular markers that will assist in future trait and disease investigations in dogs.     

Identification of Genome-Wide Copy Number Variations among Diverse Pig Breeds Using SNP Genotyping Arrays

Copy number variations (CNVs) are important forms of genetic variation complementary to SNPs, and can be considered as promising markers for some phenotypic and economically important traits or diseases susceptibility in domestic animals. In the present study, we performed a genome-wide CNV identification in 14 individuals selected from diverse populations, including six types of Chinese indigenous breeds, one Asian wild boar population, as well as three modern commercial foreign breeds. We identified 63 CNVRs in total, which covered 9.98 Mb of polymorphic sequence and corresponded to 0.36% of the genome sequence. The length of these CNVRs ranged from 3.20 to 827.21 kb, with an average of 158.37 kb and a median of 97.85 kb. Functional annotation revealed these identified CNVR have important molecular function, and may play an important role in exploring the genetic basis of phenotypic variability and disease susceptibility among pigs. Additionally, to confirm these potential CNVRs, we performed qPCR for 12 randomly selected CNVRs and 8 of them (66.67%) were confirmed successfully. CNVs detected in diverse populations herein are essential complementary to the CNV map in the pig genome, which provide an important resource for studies of genomic variation and the association between various economically important traits and CNVs.     

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.     

Genome‐wide association study identifies variants in the CAPN9 gene associated with umbilical hernia in pigs

Pig umbilical hernia (UH) affects pig welfare and brings considerable economic loss to the pig industry. To date, the molecular mechanisms underlying pig UH are still poorly understood. To identify potential loci for susceptibility to this disease, we performed a genome‐wide association study in an Erhualian × Shaziling F₂ intercross population. A total of 45 animals were genotyped using Illumina Porcine SNP60 BeadChips. We observed a SNP (rs80993347) located in the calpain‐9 (CAPN9) gene on Sus scrofa chromosome 14 that was significantly associated with UH (P = 1.97 × 10^?10). Then, we identified a synonymous mutation rs321865883 (g.20164T>C) in exon 10 of the CAPN9 gene that distinguished two affected individuals (CC) from their normal full‐sibs (TC). Finally, quantitative polymerase chain reaction was explored to investigate the mRNA expression profile of the CAPN9 gene in 12 tissues in Yorkshire pigs at different developmental stages (3, 90 and 180 days). CAPN9 showed high expression levels in the gastrointestinal tract at these three growth stages. The results of this study indicate that the CAPN9 gene might be implicated in UH. Further studies are required to establish a role of CAPN9 in pig UH.     

Genetic dissection of pre‐anthesis sub‐phase durations during the reproductive spike development of wheat

Flowering time is an important factor affecting grain yield in wheat. In this study, we divided reproductive spike development into eight sub‐phases. These sub‐phases have the potential to be delicately manipulated to increase grain yield. We measured 36 traits with regard to sub‐phase durations, determined three grain yield‐related traits in eight field environments and mapped 15 696 single nucleotide polymorphism (SNP, based on 90k Infinium chip and 35k Affymetrix chip) markers in 210 wheat genotypes. Phenotypic and genetic associations between grain yield traits and sub‐phase durations showed significant consistency (Mantel test; r = 0.5377, P < 0.001). The shared quantitative trait loci (QTLs) revealed by the genome‐wide association study suggested a close association between grain yield and sub‐phase duration, which may be attributed to effects on spikelet initiation/spikelet number (double ridge to terminal spikelet stage, DR‐TS) and assimilate accumulation (green anther to anthesis stage, GA‐AN). Moreover, we observed that the photoperiod‐sensitivity allele at the Ppd‐D1 locus on chromosome 2D markedly extended all sub‐phase durations, which may contribute to its positive effects on grain yield traits. The dwarfing allele at the Rht‐D1 (chromosome 4D) locus altered the sub‐phase duration and displayed positive effects on grain yield traits. Data for 30 selected genotypes (from among the original 210 genotypes) in the field displayed a close association with that from the greenhouse. Most importantly, this study demonstrated specific connections to grain yield in narrower time windows (i.e. the eight sub‐phases), rather than the entire stem elongation phase as a whole.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

ICAM‐1 and MKL‐1 polymorphisms impose considerable impacts on coronary heart disease occurrence

This study was aimed to explore the correlation of intercellular adhesion molecule‐1 (ICAM‐1) K469E and megakaryoblastic leukaemia factor‐1 (MKL‐1) ?184C/T polymorphisms with the susceptibility to coronary heart disease (CHD) in the Chinese Han population. 100 CHD patients and 91 healthy people that had no blood connection with each other were enrolled in this case‐control study. ICAM‐1 and MKL‐1 polymorphisms were genotyped by polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) approach. Multiple logistic regression was used to analyse the correlation between polymorphisms of ICAM‐1 and MKL‐1 and CHD susceptibility. Differences of genotype and allele frequencies of the two SNPs between case and control groups were analysed by chi‐square test. Odds ratios (ORs) and 95% confidence intervals (CIs) were indicated relative susceptibility of CHD. The distributions of ICAM‐1 and MKL‐1 polymorphisms in each group conformed to Hardy‐Weinberg equilibrium (HWE). After adjusting for traditional risk factors, the TT genotype frequency of MKL‐1 ?184C/T polymorphism was found significantly higher in case group than in control group (P < .05). Meanwhile, T allele frequency increased in case group compared with control group, and the differences had statistical significance (P = .04, OR = 2.34, 95% CI = 1.34‐5.26). Logistic regression analysis in this study proved that smoking, hypertension, diabetes and triglyceride (TG) were all risk factors for CHD ICAM‐1 K469E polymorphism has no association with the onset of CHD. But MKL‐1 ?184C/T polymorphism is associated with the risk of CHD and T allele might be a susceptibility factor for CHD.     

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.     

Validation of genetic polymorphisms on BTA14 associated with carcass trait in a commercial Hanwoo population

The objective of this study was to validate the association of significant SNPs identified from a previous genome‐wide association study with carcass weight (CWT) in a commercial Hanwoo population. We genotyped 13 SNPs located on BTA14 in 867 steers from Korea Hanwoo feedlot bulls. Of these 13 SNPs, five SNPs, namely rs29021868, rs110061498, rs109546980, rs42404006 and rs42303720, were found to be significantly associated (P < 0.001) with CWT. These five significant markers spanned the 24.3 to 29.4 Mb region of BTA14. The most significant marker (rs29021868) for CWT in this study had a 13.07 kg allele substitution effect and accounted for 2.4% of the additive genetic variance in the commercial Hanwoo population. The SNP marker rs109546980 was found to be significantly associated with both CWT (P < 0.001) and eye muscle area (P < 0.001) and could potentially be exploited for marker‐assisted selection in Hanwoo cattle. We also genotyped the ss319607402 variation, which maps to intron2 of PLAG1 gene and which is already reported to be associated with height, to identify any significant association with carcass weight; however, no such association was observed in this Hanwoo commercial population.