Copy number variation-based genome wide association study reveals additional variants contributing to meat quality in Swine

Mapping and functional characterization of structural variation in 1060 pig genomes

BACKGROUND

Structural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence.

RESULTS

We present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies.

CONCLUSIONS

This SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

Finding Predictors of Leg Defects in Pigs Using CNV-GWAS

One of the most important areas of modern genome research is the search for meaningful relationships between genetic variants and phenotypes. In the livestock field, there has been research demonstrating the influence of copy number variants (CNVs) on phenotypic variation. Despite the wide range in the number and size of detected CNVs, a significant proportion differ between breeds and their functional effects are underestimated in the pig industry. In this work, we focused on the problem of leg defects in pigs (lumps/growths in the area of the hock joint on the hind legs) and focused on searching for molecular genetic predictors associated with this trait for the selection of breeding stock. The study was conducted on Large White pigs using three CNV calling tools (PennCNV, QuantiSNP and R-GADA) and the CNVRanger association analysis tool (CNV-GWAS). As a result, the analysis identified three candidate CNVRs associated with the formation of limb defects. Subsequent functional analysis suggested that all identified CNVs may act as potential predictors of the hock joint phenotype of pigs. It should be noted that the results obtained indicate that all significant regions are localized in genes (CTH, SRSF11, MAN1A1 and LPIN1) responsible for the metabolism of amino acids, fatty acids, glycerolipids and glycerophospholipids, thereby related to the immune response, liver functions, content intramuscular fat and animal fatness. These results are consistent with previously published studies, according to which a predisposition to the formation of leg defects can be realized through genetic variants associated with the functions of the liver, kidneys and hematological characteristics.

Detection distribution of CNVs of SNX29 in three goat breeds and their associations with growth traits

As a member of the SNX family, the goat sorting nexin 29 (SNX29) is initially identified as a myogenesis gene. Therefore, this study aimed to examine the polymorphism in the SNX29 gene and its association with growth traits. In this study, we used an online platform to predict the structures of the SNX29 protein and used quantitative real-time PCR to detect potential copy number variation (CNV) in Shaanbei white cashmere (SBWC) goats (n = 541), Guizhou black (GB) goats (n = 48), and Nubian (NB) goats (n = 39). The results showed that goat SNX29 protein belonged to non-secretory protein. Then, five CNVs were detected, and their association with growth traits was analyzed. In SBWC goats, CNV1, CNV3, CNV4, and CNV5 were associated with chest width and body length (P < 0.05). Among them, the CNV1 individuals with gain and loss genotypes were superior to those individuals with a median genotype, but CNV4 and CNV5 of individuals with the median genotype were superior to those with the loss and gain genotypes. In addition, individuals with the gain genotype had superior growth traits in CNV3. In brief, this study suggests that the CNV of SNX29 can be used as a molecular marker in goat breeding.

Copy Number Variation Regions Differing in Segregation Patterns Span Different Sets of Genes

Copy number variations regions (CNVRs) can be classified either as segregating, when found in both parents, and offspring, or non-segregating. A total of 65 segregating and 31 non-segregating CNVRs identified in at least 10 individuals within a dense pedigree of the Gochu Asturcelta pig breed was subjected to enrichment and functional annotation analyses to ascertain their functional independence and importance. Enrichment analyses allowed us to annotate 1018 and 351 candidate genes within the bounds of the segregating and non-segregating CNVRs, respectively. The information retrieved suggested that the candidate genes spanned by segregating and non-segregating CNVRs were functionally independent. Functional annotation analyses allowed us to identify nine different significantly enriched functional annotation clusters (ACs) in segregating CNVR candidate genes mainly involved in immunity and regulation of the cell cycle. Up to five significantly enriched ACs, mainly involved in reproduction and meat quality, were identified in non-segregating CNVRs. The current analysis fits with previous reports suggesting that segregating CNVRs would explain performance at the population level, whereas non-segregating CNVRs could explain between-individuals differences in performance.

Pangenome obtained by long-read sequencing of 11 genomes reveal hidden functional structural variants in pigs

Long-read sequencing (LRS) facilitates both the genome assembly and the discovery of structural variants (SVs). Here, we built a graph-based pig pangenome by incorporating 11 LRS genomes with an average of 94.01% BUSCO completeness score, revealing 206-Mb novel sequences. We discovered 183,352 nonredundant SVs (63% novel), representing 12.12% of the reference genome. By genotyping SVs in an additional 196 short-read sequencing samples, we identified thousands of population stratified SVs. Particularly, we detected 7,568 Tibetan specific SVs, some of which demonstrate significant population differentiation between Tibetan and low-altitude pigs, which might be associated with the high-altitude hypoxia adaptation in Tibetan pigs. Further integrating functional genomic data, the most promising candidate genes within the SVs that might contribute to the high-altitude hypoxia adaptation were discovered. Overall, our study generates a benchmark pangenome resource for illustrating the important roles of SVs in adaptive evolution, domestication, and genetic improvement of agronomic traits in pigs.

Genome-Wide Detection and Analysis of Copy Number Variation in Anhui Indigenous and Western Commercial Pig Breeds Using Porcine 80K SNP BeadChip

Copy number variation (CNV) is an important class of genetic variations widely associated with the porcine genome, but little is known about the characteristics of CNVs in foreign and indigenous pig breeds. We performed a genome-wide comparison of CNVs between Anhui indigenous pig (AHIP) and Western commercial pig (WECP) breeds based on data from the Porcine 80K SNP BeadChip. After analysis using the PennCNV software, we detected 3863 and 7546 CNVs in the AHIP and WECP populations, respectively. We obtained 225 (loss: 178, gain: 47) and 379 (loss: 293, gain: 86) copy number variation regions (CNVRs) randomly distributed across the autosomes of the AHIP and WECP populations, accounting for 10.90% and 22.57% of the porcine autosomal genome, respectively. Functional enrichment analysis of genes in the CNVRs identified genes related to immunity (FOXJ1, FOXK2, MBL2, TNFRSF4, SIRT1, NCF1) and meat quality (DGAT1, NT5E) in the WECP population; these genes were a loss event in the WECP population. This study provides important information on CNV differences between foreign and indigenous pig breeds, making it possible to provide a reference for future improvement of these breeds and their production performance.

Segregation patterns and inheritance rate of copy number variations regions assessed in a Gochu Asturcelta pig pedigree

Copy Number Variation Regions (CNVR) were subjected to pedigree analysis to contribute to the understanding of their segregation patterns. Up to 492 Gochu Asturcelta pig individuals forming 478 different parents-offspring trios (61 different families) were genotyped using the Axiom_PigHDv1 Array (658,692 SNPs). CNVR calling, performed using two different platforms (PennCNV and QuantiSNP), allowed to identify a total of 344 candidate CNVR on the 18 porcine autosomes covering about 106.8 Mb of the pig genome. Sixty-nine CNVR were identified, to some extent, in both the parents and the offspring and were classified as segregating CNVR. The other candidate CNVR were called in one or more progeny but in neither parent and classified either as singleton or recurrent de novo CNVR. Segregating CNVR were, on average, larger and more frequent than the recurrent de novo CNVR (444.8 kb vs 287.9 kb long and 34 vs 5 individuals, respectively). In any case, segregating CNVR did not conform to strict Mendelian inheritance patterns: estimates of average paternal and maternal transmission rates ranged from 11.0 % to 13.4 % and mean inheritance rate was below 21 %. This issue should be carefully considered when interpreting the results of CNV studies. Segregating CNVR, present across generations, are unlikely to be artifacts or false positives and can be hypothesized to be important at the population level.

Analysis of Copy Number Variation in the Whole Genome of Normal-Haired and Long-Haired Tianzhu White Yaks

Long-haired individuals in the Tianzhu white yak population are a unique genetic resource, and have important landscape value. Copy number variation (CNV) is an important source of phenotypic variation in mammals. In this study, we used resequencing technology to detect the whole genome of 10 long-haired Tianzhu white yaks (LTWY) and 10 normal-haired Tianzhu white yaks (NTWY), and analyzed the differences of CNV in the genome of LTWYs and NTWYs. A total of 110268 CNVs were identified, 2006 CNVRs were defined, and the distribution map of these CNVRs on chromosomes was constructed. The comparison of LTWYs and NTWYs identified 80 differential CNVR-harbored genes, which were enriched in lipid metabolism, cell migration and other functions. Notably, some differential genes were identified as associated with hair growth and hair-follicle development (e.g., ASTN2, ATM, COL22A1, GK5, SLIT3, PM20D1, and SGCZ). In general, we present the first genome-wide analysis of CNV in LTWYs and NTWYs. Our results can provide new insights into the phenotypic variation of different hair lengths in Tianzhu white yaks.

Copy Number Variations Contribute to Intramuscular Fat Content Differences by Affecting the Expression of PELP1 Alternative Splices in Pigs

Simple Summary

Copy number variation (CNV) is a type of variant that may influence meat quality of, for example intramuscular fat (IMF). In this study, a genome-wide association study (GWAS) was then performed between CNVs and IMF in a pig F2 resource population. A total of 19 CNVRs were found to be significantly associated with IMF. RNA-seq and qPCR validation results indicated that CNV150, which is located on the 3′UTR end of the proline, as well as glutamate and the leucine rich protein 1 (PELP1) gene may affect the expression of PELP1 alternative splices. We infer that the CNVR may influence IMF content by regulating the alternative splicing of the PELP1 gene and ultimately affects the structure of the PELP1 protein. These findings suggest a novel mechanistic approach for meat quality improvement in animals and the potential treatment of insulin resistance in human beings.

Abstract

Intramuscular fat (IMF) is a key meat quality trait. Research on the genetic mechanisms of IMF decomposition is valuable for both pork quality improvement and the treatment of obesity and type 2 diabetes. Copy number variations (CNVs) are a type of variant that may influence meat quality. In this study, a total of 1185 CNV regions (CNVRs) including 393 duplicated CNVRs, 432 deleted CNVRs, and 361 CNVRs with both duplicated and deleted status were identified in a pig F2 resource population using next-generation sequencing data. A genome-wide association study (GWAS) was then performed between CNVs and IMF, and a total of 19 CNVRs were found to be significantly associated with IMF. QTL colocation analysis indicated that 3 of the 19 CNVRs overlapped with known QTLs. RNA-seq and qPCR validation results indicated that CNV150, which is located on the 3′UTR end of the proline, as well as glutamate and the leucine rich protein 1 (PELP1) gene may affect the expression of PELP1 alternative splices. Sequence alignment and Alphafold2 structure prediction results indicated that the two alternative splices of PELP1 have a 23 AA sequence variation and a helix-fold structure variation. This region is located in the region of interaction between PELP1 and other proteins which have been reported to be significantly associated with fat deposition or insulin resistance. We infer that the CNVR may influence IMF content by regulating the alternative splicing of the PELP1 gene and ultimately affects the structure of the PELP1 protein. In conclusion, we found some CNVRs, especially CNV150, located in PELP1 that affect IMF. These findings suggest a novel mechanistic approach for meat quality improvement in animals and the potential treatment of insulin resistance in human beings.

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 association study reveals a quantitative trait locus and two candidate genes on Sus scrofa chromosome 5 affecting intramuscular fat content in Suhuai pigs

Intramuscular fat content (IFC) is an essential quantitative trait of meat, affecting multiple meat quality indicators. A certain amount of IFC could not only improve the sensory score of pork but also increase the flavour, tenderness, juiciness and shelf-life. To dissect the genetic determinants of IFC, two methods, including genome-wide efficient mixed-model analysis (GEMMA) and linkage disequilibrium adjusted kinships (LDAKs), were used to carry out genome-wide association studies for IFC in Suhuai pig population. A total of 14 and 18 significant single nucleotide polymorphisms (SNPs) were identified by GEMMA and LDAK, respectively. The results of these two methods were highly consistent and all 14 significant SNPs in GEMMA were detected by LDAK. Seven of the 18 SNPs reached the genome-wide significance level (P < 9.85E-07) while 11 cases reached the suggestive significance level (P < 1.77E-05). These significant SNPs were mainly distributed on Sus scrofa chromosome (SSC) 5, 3, and 7. Moreover, one locus resides in a 2.27 Mb (71.37-73.64 Mb) region on SSC5 harbouring 13 significant SNPs associated with IFC, and the lead SNP (rs81302978) also locates in this region. Linkage disequilibrium (LD) analysis showed that there were four pairs of complete LD (r² = 1) among these 13 SNPs, and the remaining 9 SNPs with incomplete LD (r² ≠ 1) were selected for subsequent analyses of IFC. Association analyses showed that 7 out of 9 SNPs were significantly associated with IFC (P < 0.05) in 330 Suhuai pigs, and the other 2 SNPs tended to reach a significant association level with IFC (P < 0.1). The phenotypic variance explained (PVE) range of these 9 SNPs was 0.92-3.55%. Meanwhile, the lead SNP was also significantly associated (rs81302978) with IFC (P < 0.05) in 378 commercial hybrid pigs (Pietrain × Duroc) × (Landrace × Yorkshire) (PDLY), and the PVE was 1.38%. Besides, two lipid metabolism-relevant candidate genes, the leucine rich repeat kinase 2 (LRRK2) and PDZ domain containing ring finger 4 (PDZRN4) were identified in the 2.27 Mb region on SSC5. In conclusion, our results may provide a set of markers useful for genetic improvement of IFC in pigs and will advance the genome selection process of IFC on pig breeding programmes.

Classical, Molecular, and Genomic Cytogenetics of the Pig, a Clinical Perspective

The chromosomes of the domestic pig (Sus scrofa domesticus) are known to be prone to reciprocal chromosome translocations and other balanced chromosome rearrangements with concomitant fertility impairment of carriers. In response to the remarkable prevalence of chromosome rearrangements in swine herds, clinical cytogenetics laboratories have been established in several countries in order to screen young boars for chromosome rearrangements prior to service. At present, clinical cytogenetics laboratories typically apply classical cytogenetics techniques such as giemsa-trypsin (GTG)-banding to produce high-quality karyotypes and reveal large-scale chromosome ectopic exchanges. Further refinements to clinical cytogenetics practices have led to the implementation of molecular cytogenetics techniques such as fluorescent in-situ hybridization (FISH), allowing for rearrangements to be visualized and breakpoints refined using fluorescently labelled painting probes. The next-generation of clinical cytogenetics include the implementation of DNA microarrays, and next-generation sequencing (NGS) technologies such as DNA sequencing to better explore tentative genome architecture changes. The implementation of these cytogenomics techniques allow the genomes of rearrangement carriers to be deciphered at the highest resolution, allowing rearrangements to be detected; breakpoints to be delineated; and, most importantly, potential gene implications of those chromosome rearrangements to be interrogated. Clinical cytogenetics has become an integral tool in the livestock industry, identifying rearrangements and allowing breeders to make informed breeding decisions.

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.

Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle

Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22-ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.

Detection and analysis of genome-wide copy number variation in the pig genome using an 80 K SNP Beadchip

Copy number variation (CNV) is an important source of genetic variability in human or animal genomes and play key roles in phenotypic diversity and disease susceptibility. In the present study, we performed a genome-wide analysis for CNV detection using SNP genotyping data of 857 Large White pigs. A total of 312 CNV regions (CNVRs) were detected with the PennCNV algorithm, which covered 57.76 Mb of the pig genome and correspond to 2.36% of the genome sequence. The length of the CNVRs on autosomes ranged from 1.77 Kb to 1.76 Mb with an average of 185.11 Kb. Of these, 220 completely or partially overlapped with 1,092 annotated genes, which enriched a wide variety of biological processes. Comparisons with previously reported pig CNVR revealed 92 (29.49%) novel CNVRs. Experimentally, 80% of CNVRs selected randomly were validated by quantitative PCR (qPCR). We also performed an association analysis between some of the CNVRs and reproductive traits, with results demonstrating the potential importance of CNVR61 and CNVR283 associated with litter sizes. Notably, the GPER1 gene located in CNVR61 plays a key role in reproduction. Our study is an important complement to the CNV map in the pig genome and provides valuable information for investigating the association between genomic variation and economic traits.

NTN1 Affects Porcine Intramuscular Fat Content by Affecting the Expression of Myogenic Regulatory Factors

Simple Summary

Intramuscular fat (IMF) is a key meat quality trait in the pork industry. In this study, we validated the effect of the copy number of Netrin-1 (NTN1-CNV) on Netrin-1 (NTN1) protein expression and explored the possible affective mechanism of NTN1 on IMF. The results indicated that NTN1-CNV may affect the expression of NTN1 protein by its gene dose, and the expression of NTN1 may affect the proliferation and differentiation of muscle cells by the AMP-activated protein kinase (AMPK) pathway and finally influence the IMF content.

Abstract

Intramuscular fat (IMF) content is an important economic trait for pork quality. Our previous results regarding the genome-wide association between IMF content and copy number variations (CNVs) indicated that the CNV within Netrin-1(NTN1-CNV) was significantly associated with IMF. In order to validate the effect of NTN1-CNV, we detected the Netrin-1 (NTN1) gene dose and protein expression content in the longissimus dorsi of different IMF content pigs using Western blotting and investigated the expression of NTN1 RNA in different tissues using real-time quantitative polymerase chain reaction (qPCR). The knock-down of the NTN1 gene in C2C12 and 3T3-L1 cells and over-expression in C2C12 cells during the proliferation and differentiation stage were also investigated to explore the possible pathway of action of NTN1. The results showed that in individuals with IMF content differences, the gene dose of NTN1 and the expression of NTN1 protein were also significantly different, which indicated that NTN1-CNV may directly affect IMF by its coding protein. NTN1 had the highest expression in pig longissimus dorsi and backfat tissues, which indicates that NTN1 may play an important role in muscle and fat tissues. The in vitro validation assay indicated that NTN1 silencing could promote the proliferation and inhibit the differentiation of C2C12 cells, with no effect on 3T3-L1 cells. Additionally, NTN1 over-expression could inhibit the proliferation and promote the differentiation of C2C12 cells. Combined with previous research, we conclude that NTN1-CNV may affect IMF by its gene dose, and the expression of NTN1 may affect the proliferation and differentiation of muscle cells by the AMP-activated protein kinase (AMPK) pathway and finally influence the IMF.

Integrative analysis of APOL3 gene CNV for adult cattle stature

Copy number variations (CNVs) have been identified as another important structural variation of genome. In recent years, a large amount of CNVRs have been identified in humans and animals. However, association and dosage effects studies of CNVs are very limited. Apolipoprotein L3 (APOL3) gene plays a central role in modulating gene transcription and is located within a CNVR that encompasses quantitative trait locis (QTLs) for economic traits like meat quality. Herein, we analyzed the CNV polymorphism of APOL3 in 421 individuals from five distinct cattle breeds, and then correlated their genotypes with growth traits. Association analysis revealed that the APOL3 CNV was significantly associated with hip height and cannon circumference of Xianan (XN) cattle (P < .01), and visibly associated with body slanting length and hucklebone width of Pinan (PN) cattle (P < .05). Overall, the data provide evidence for the functional role of APOL3 CNV and a basis for future applications in cattle breeding.

Genome-wide association study between CNVs and milk production traits in Valle del Belice sheep

Copy number variation (CNV) is a major source of genomic structural variation. The aim of this study was to detect genomic CNV regions (CNVR) in Valle del Belice dairy sheep population and to identify those affecting milk production traits. The GO analysis identified possible candidate genes and pathways related to the selected traits. We identified CNVs in 416 individuals genotyped using the Illumina OvineSNP50 BeadChip array. The CNV association using a correlation-trend test model was examined with the Golden Helix SVS 8.7.0 tool. Significant CNVs were detected when their adjusted p-value was <0.01 after false discovery rate (FDR) correction. We identified 7,208 CNVs, which gave 365 CNVRs after aggregating overlapping CNVs. Thirty-one CNVRs were significantly associated with one or more traits included in the analysis. All CNVRs, except those on OAR19, overlapped with quantitative trait loci (QTL), even if they were not directly related to the traits of interest. A total of 222 genes were annotated within the significantly associated CNVRs, most of which played important roles in biological processes related to milk production and health-related traits. Identification of the genes in the CNVRs associated with the studied traits will provide the basis for further investigation of their role in the metabolic pathways related to milk production and health traits.

Genome-wide patterns of copy number variations in Spodoptera litura

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.

Genome-wide copy number variant analysis reveals variants associated with 10 diverse production traits in Holstein cattle

BACKGROUND

Copy number variation (CNV) is an important type of genetic variation contributing to phenotypic differences among mammals and may serve as an alternative molecular marker to single nucleotide polymorphism (SNP) for genome-wide association study (GWAS). Recently, GWAS analysis using CNV has been applied in livestock, although few studies have focused on Holstein cattle.

RESULTS

We describe 191 CNV detected using intensity data from over 700,000 SNP genotypes generated with the BovineHD Genotyping BeadChip (Illumina, San Diego, CA) in 528 Holstein cows. The CNV were used for GWAS analysis of 10 important production traits of 473 cattle related to feed intake, milk quality, and female fertility, as well as 2 composite traits of net merit and productive life. In total, we detected 57 CNV associated (P < 0.05 after false discovery rate correction) with at least one of the 10 phenotypes. Focusing on feed efficiency and intake-related phenotypes of residual feed intake and dry matter intake, we detected a single CNV associated with both traits which overlaps a predicted olfactory receptor gene OR2A2 (LOC787786). Additionally, 2 CNV within the RXFP4 (relaxin/insulin like family peptide receptor 4) and 2 additional olfactory receptor gene regions, respectively, were associated with residual feed intake. The RXFP4 gene encodes a receptor for an orexigenic peptide, insulin-like peptide 5 produced by intestinal L cells, which is expressed by enteric neurons. Olfactory receptors are critical for transmitting the effects of odorants, contributing to the sense of smell, and have been implicated in participating in appetite regulation.

CONCLUSIONS

Our results identify CNV for genomic evaluation in Holstein cattle, and provide candidate genes, such as RXFP4, contributing to variation in feed efficiency and feed intake-related traits. These results indicate potential novel targets for manipulating feed intake-related traits of livestock.

A global analysis of CNVs in swine using whole genome sequence data and association analysis with fatty acid composition and growth traits

Copy number variations (CNVs) are important genetic variants complementary to SNPs, and can be considered as biomarkers for some economically important traits in domestic animals. In the present study, a genomic analysis of porcine CNVs based on next-generation sequencing data was carried out to identify CNVs segregating in an Iberian x Landrace backcross population and study their association with fatty acid composition and growth-related traits. A total of 1,279 CNVs, including duplications and deletions, were detected, ranging from 106 to 235 CNVs across samples, with an average of 183 CNVs per sample. Moreover, we detected 540 CNV regions (CNVRs) containing 245 genes. Functional annotation suggested that these genes possess a great variety of molecular functions and may play a role in production traits in commercial breeds. Some of the identified CNVRs contained relevant functional genes (e.g., CLCA4, CYP4X1, GPAT2, MOGAT2, PLA2G2A and PRKG1, among others). The variation in copy number of four of them (CLCA4, GPAT2, MOGAT2 and PRKG1) was validated in 150 BC1_LD (25% Iberian and 75% Landrace) animals by qPCR. Additionally, their contribution regarding backfat and intramuscular fatty acid composition and growth–related traits was analyzed. Statistically significant associations were obtained for CNVR112 (GPAT2) for the C18:2(n-6)/C18:3(n-3) ratio in backfat and carcass length, among others. Notably, GPATs are enzymes that catalyze the first step in the biosynthesis of both triglycerides and glycerophospholipids, suggesting that this CNVR may contribute to genetic variation in fatty acid composition and growth traits. These findings provide useful genomic information to facilitate the further identification of trait-related CNVRs affecting economically important traits in pigs.

CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation

Background

Copy number variations (CNVs) have been linked to different phenotypes in human, including many diseases. A genome-scale understanding of CNVs is available in a few plants but none are wild species, leaving a knowledge gap regarding their genome biology and evolutionary role. We developed a reliable CNV detection method for species lacking contiguous reference genome. We selected multiple probes within 14,078 gene sequences and developed comparative genome hybridization on arrays. Gene CNVs were assessed in three full-sib families from species with 20 Gb genomes, i.e., white and black spruce, and interior spruce - a natural hybrid.

Results

We discovered hundreds of gene CNVs in each species, 3612 in total, which were enriched in functions related to stress and defense responses and narrow expression profiles, indicating a potential role in adaptation. The number of shared CNVs was in accordance with the degree of relatedness between individuals and species. The genetically mapped subset of these genes showed a wide distribution across the genome, implying numerous structural variations. The hybrid family presented significantly fewer CNVs, suggesting that the admixture of two species within one genome reduces the occurrence of CNVs.

Conclusions

The approach we developed is of particular interest in non-model species lacking a reference genome. Our findings point to a role for CNVs in adaptation. Their reduced abundance in the hybrid may limit genetic variability and evolvability of hybrids. We propose that CNVs make a qualitatively distinct contribution to adaptation which could be important for short term change.

Electronic supplementary material

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

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

Genome-wide patterns of copy number variation in the Chinese yak genome

Background

Copy number variation (CNV) represents an important source of genetic divergence that can produce drastic phenotypic differences and may therefore be subject to selection during domestication and environmental adaptation. To investigate the evolutionary dynamics of CNV in the yak genome, we used a read depth approach to detect CNV based on genome resequencing data from 14 wild and 65 domestic yaks and determined CNV regions related to domestication and adaptations to high-altitude.

Results

We identified 2,634 CNV regions (CNVRs) comprising a total of 153 megabases (5.7 % of the yak genome) and 3,879 overlapping annotated genes. Comparison between domestic and wild yak populations identified 121 potentially selected CNVRs, harboring genes related to neuronal development, reproduction, nutrition and energy metabolism. In addition, we found 85 CNVRs that are significantly different between domestic yak living in high- and low-altitude areas, including three genes related to hypoxia response and six related to immune defense. This analysis shows that genic CNVs may play an important role in phenotypic changes during yak domestication and adaptation to life at high-altitude.

Conclusions

We present the first refined CNV map for yak along with comprehensive genomic analysis of yak CNV. Our results provide new insights into the genetic basis of yak domestication and adaptation to living in a high-altitude environment, as well as a valuable genetic resource that will facilitate future CNV association studies of important traits in yak and other bovid species.

Electronic supplementary material

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

SNP- and haplotype-based genome-wide association studies for growth, carcass, and meat quality traits in a Duroc multigenerational population

Background

The aim of the present study was to compare the power of single nucleotide polymorphism (SNP)-based genome-wide association study (GWAS) and haplotype-based GWAS for quantitative trait loci (QTL) detection, and to detect novel candidate genes affecting economically important traits in a purebred Duroc population comprising seven-generation pedigree. First, we performed a simulation analysis using real genotype data of this population to compare the power (based on the null hypothesis) of the two methods. We then performed GWAS using both methods and real phenotype data comprising 52 traits, which included growth, carcass, and meat quality traits.

Results

In total, 836 animals were genotyped using the Illumina PorcineSNP60 BeadChip and 14 customized SNPs from regions of known candidate genes related to the traits of interest. The power of SNP-based GWAS was greater than that of haplotype-based GWAS in a simulation analysis. In real data analysis, a larger number of significant regions was obtained by SNP-based GWAS than by haplotype-based GWAS. For SNP-based GWAS, 23 genome-wide significant SNP regions were detected for 17 traits, and 120 genome-wide suggestive SNP regions were detected for 27 traits. For haplotype-based GWAS, 6 genome-wide significant SNP regions were detected for four traits, and 11 genome-wide suggestive SNP regions were detected for eight traits. All genome-wide significant SNP regions detected by haplotype-based GWAS were located in regions also detected by SNP-based GWAS. Four regions detected by SNP-based GWAS were significantly associated with multiple traits: on Sus scrofa chromosome (SSC) 1 at 304 Mb; and on SSC7 at 35–39 Mb, 41–42 Mb, and 103 Mb. The vertnin gene (VRTN) in particular, was located on SSC7 at 103 Mb and was significantly associated with vertebrae number and carcass lengths. Mapped QTL regions contain some candidate genes involved in skeletal formation (FUBP3; far upstream element binding protein 3) and fat deposition (METTL3; methyltransferase like 3).

Conclusion

Our results show that a multigenerational pig population is useful for detecting QTL, which are typically segregated in a purebred population. In addition, a novel significant region could be detected by SNP-based GWAS as opposed to haplotype-based GWAS.

Electronic supplementary material

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