Genome-Wide Copy Number Variations Using SNP Genotyping in a Mixed Breed Swine Population

Genome wide copy number variations using Porcine 60K SNP Beadchip in Landlly pigs

In the present study, Porcine 60K SNP genotype data from 69 Landlly pigs were used to explore Copy Number Variations (CNVs) across the autosomes. A total of 386 CNVs were identified using Hidden Markov Model (HMM) in PennCNV software, which were subsequently aggregated to 115 CNV regions (CNVRs). Among the total detected CNVRs, 58 gain, 49 were loss type while remaining 8 events were both gain and loss types. Identified CNVRs covered 12.5 Mb (0.55%) of Sus scrofa reference 11.1 genome. Comparison of our results with previous investigations on pigs revealed that approximately 75% CNVRs were novel, which may be due to differences in genetic background, environment and implementation of artificial selection in Landlly pigs. Functional annotation and pathway analysis showed the significant enrichment of 267 well-annotated Sus scrofa genes in CNVRs. These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity. Additionally, KIT and FUT1 were two major genes detected on CNVR in our population. This investigation provided a comprehensive overview of CNV distribution in the Indian porcine genome for the first time, which may be useful for further investigating the association of important quantitative traits in Landlly pigs.Highlights115 CNVRs were identified in 69 Landlly pig population.Approximately 75% detected CNVRs were novel for Landlly population.Significant enrichment of 267 well-annotated Sus scrofa genes observed in these CNVRs.These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity.Comprehensive CNV map in the Indian porcine genome developed for the first time.

Genome-Wide Assessment Characteristics of Genes Overlapping Copy Number Variation Regions in Duroc Purebred Population

Copy number variations (CNVs) are important structural variations that can cause significant phenotypic diversity. Reliable CNVs mapping can be achieved by identification of CNVs from different genetic backgrounds. Investigations on the characteristics of overlapping between CNV regions (CNVRs) and protein-coding genes (CNV genes) or miRNAs (CNV-miRNAs) can reveal the potential mechanisms of their regulation. In this study, we used 50 K SNP arrays to detect CNVs in Duroc purebred pig. A total number of 211 CNVRs were detected with a total length of 118.48 Mb, accounting for 5.23% of the autosomal genome sequence. Of these CNVRs, 32 were gains, 175 losses, and four contained both types (loss and gain within the same region). The CNVRs we detected were non-randomly distributed in the swine genome and were significantly enriched in the segmental duplication and gene density region. Additionally, these CNVRs were overlapping with 1,096 protein-coding genes (CNV-genes), and 39 miRNAs (CNV-miRNAs), respectively. The CNV-genes were enriched in terms of dosage-sensitive gene list. The expression of the CNV genes was significantly higher than that of the non-CNV genes in the adult Duroc prostate. Of all detected CNV genes, 22.99% genes were tissue-specific (TSI > 0.9). Strong negative selection had been underway in the CNV-genes as the ones that were located entirely within the loss CNVRs appeared to be evolving rapidly as determined by the median dN plus dS values. Non-CNV genes tended to be miRNA target than CNV-genes. Furthermore, CNV-miRNAs tended to target more genes compared to non-CNV-miRNAs, and a combination of two CNV-miRNAs preferentially synergistically regulated the same target genes. We also focused our efforts on examining CNV genes and CNV-miRNAs functions, which were also involved in the lipid metabolism, including DGAT1, DGAT2, MOGAT2, miR143, miR335, and miRLET7. Further molecular experiments and independent large studies are needed to confirm our findings.

Genome-wide detection of CNV regions and their potential association with growth and fatness traits in Duroc pigs

Background

In the process of pig breeding, the average daily gain (ADG), days to 100 kg (AGE), and backfat thickness (BFT) are directly related to growth rate and fatness. However, the genetic mechanisms involved are not well understood. Copy number variation (CNV), an important source of genetic diversity, can affect a variety of complex traits and diseases and has gradually been thrust into the limelight. In this study, we reported the genome-wide CNVs of Duroc pigs using SNP genotyping data from 6627 animals. We also performed a copy number variation region (CNVR)-based genome-wide association studies (GWAS) for growth and fatness traits in two Duroc populations.

Results

Our study identified 953 nonredundant CNVRs in U.S. and Canadian Duroc pigs, covering 246.89 Mb (~ 10.90%) of the pig autosomal genome. Of these, 802 CNVRs were in U.S. Duroc pigs with 499 CNVRs were in Canadian Duroc pigs, indicating 348 CNVRs were shared by the two populations. Experimentally, 77.8% of nine randomly selected CNVRs were validated through quantitative PCR (qPCR). We also identified 35 CNVRs with significant association with growth and fatness traits using CNVR-based GWAS. Ten of these CNVRs were associated with both ADG and AGE traits in U.S. Duroc pigs. Notably, four CNVRs showed significant associations with ADG, AGE, and BFT, indicating that these CNVRs may play a pleiotropic role in regulating pig growth and fat deposition. In Canadian Duroc pigs, nine CNVRs were significantly associated with both ADG and AGE traits. Further bioinformatic analysis identified a subset of potential candidate genes, including PDGFA, GPER1, PNPLA2 and BSCL2.

Conclusions

The present study provides a necessary supplement to the CNV map of the Duroc genome through large-scale population genotyping. In addition, the CNVR-based GWAS results provide a meaningful way to elucidate the genetic mechanisms underlying complex traits. The identified CNVRs can be used as molecular markers for genetic improvement in the molecular-guided breeding of modern commercial pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07654-7.

Genome-wide association study between copy number variation regions and carcass- and meat-quality traits in Nellore cattle

Context Indicine breeds are the main source of beef products in tropical and subtropical regions. However, genetic improvement for carcass- and meat-quality traits in zebu cattle have been limited and genomics studies concerning structural variations that influence these traits are essential. Aim The aim of this study was to perform a genome-wide association study between copy number variation regions (CNVRs) and carcass- and meat quality-traits in Nellore cattle. Methods In total, 3794 animals, males and females included, were genotyped using a 777962 single-nucleotide polymorphism platform of BovineHD BeadChip (777k; Illumina Inc.). Of these, 1751 Nellore bulls were slaughtered at 24 months of age for further carcass beef analysis. The following traits were studied: beef tenderness, marbling, rib-eye area, backfat thickness and meat colour (lightness, redness and yellowness). The CNV detection was conducted through PennCNV software. The association analyses were performed using CNVRuler software. Key results Several identified genomic regions were linked to quantitative trait loci associated with fat deposition (FABP7) and lipid metabolism (PPARA; PLA2 family; BCHE), extracellular matrix (INS; COL10A1), contraction (SLC34A3; TRDN) and muscle development (CAPZP). The gene-enrichment analyses highlighted biological mechanisms directly related to the metabolism and synthesis of lipids and fatty acids. Conclusions The large number of potential candidate genes identified within the CNVRs, as well as the functions and pathways identified, should help better elucidate the genetic mechanisms involved in the expression of beef and carcass traits in Nellore cattle. Several CNVRs harboured genes that might have a functional impact to improve the beef and carcass traits. Implications The results obtained contribute to upgrade the sensorial and organoleptic attributes of Nellore cattle and make feasible the genetic improvement of carcass- and meat-quality traits.

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.

A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle

BACKGROUND

The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population.

RESULTS

The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals.

CONCLUSION

Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

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.

Genetic analysis of teat number in pigs reveals some developmental pathways independent of vertebra number and several loci which only affect a specific side

Background

Number of functional teats is an important trait in commercial swine production. As litter size increases, the number of teats must also increase to supply nutrition to all piglets. Therefore, a genome-wide association analysis was conducted to identify genomic regions that affect this trait in a commercial swine population. Genotypic data from the Illumina Porcine SNP60v1 BeadChip were available for 2951 animals with total teat number (TTN) records. A subset of these animals (n = 1828) had number of teats on each side recorded. From this information, the following traits were derived: number of teats on the left (LTN) and right side (RTN), maximum number of teats on a side (MAX), difference between LTN and RTN (L − R) and absolute value of L − R (DIF). Bayes C option of GENSEL (version 4.61) and 1-Mb windows were implemented. Identified regions that explained more than 1.5% of the genomic variation were tested in a larger group of animals (n = 5453) to estimate additive genetic effects.

Results

Marker heritabilities were highest for TTN (0.233), intermediate for individual side counts (0.088 to 0.115) and virtually nil for difference traits (0.002 for L − R and 0.006 for DIF). Each copy of the VRTN mutant allele increased teat count by 0.35 (TTN), 0.16 (LTN and RTN) and 0.19 (MAX). 15, 18, 13 and 18 one-Mb windows were detected that explained more than 1.0% of the genomic variation for TTN, LTN, RTN, and MAX, respectively. These regions cumulatively accounted for over 50% of the genomic variation of LTN, RTN and MAX, but only 30% of that of TTN. Sus scrofa chromosome SSC10:52 Mb was associated with all four count traits, while SSC10:60 and SSC14:54 Mb were associated with three count traits. Thirty-three SNPs accounted for nearly 39% of the additive genetic variation in the validation dataset. No effect of piglet sex or percentage of males in litter was detected, but birth weight was positively correlated with TTN.

Conclusions

Teat number is a heritable trait and use of genetic markers would expedite selection progress. Exploiting genetic variation associated with teat counts on each side would enhance selection focused on total teat counts. These results confirm QTL on SSC4, seven and ten and identify a novel QTL on SSC14.

Electronic supplementary material

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

A study of vertebra number in pigs confirms the association of vertnin and reveals additional QTL

Background

Formation of the vertebral column is a critical developmental stage in mammals. The strict control of this process has resulted in little variation in number of vertebrae across mammalian species and no variation within most mammalian species. The pig is quite unique as considerable variation exists in number of thoracic vertebrae as well as number of lumbar vertebrae. At least two genes have been identified that affect number of vertebrae in pigs yet considerable genetic variation still exists. Therefore, a genome-wide association (GWA) analysis was conducted to identify additional genomic regions that affect this trait.

Results

A total of 1883 animals were phenotyped for the number of ribs and thoracolumbar vertebrae as well as successfully genotyped with the Illumina Porcine SNP60 BeadChip. After data editing, 41,148 SNP markers were included in the GWA analysis. These animals were also phenotyped for kyphosis. Fifty-three 1 Mb windows each explained at least 1.0 % of the genomic variation for vertebrae counts while 16 regions were significant for kyphosis. Vertnin genotype significantly affected vertebral counts as well. The region with the largest effect for number of lumbar vertebrae and thoracolumbar vertebrae were located over the Hox B gene cluster and the largest association for thoracic vertebrae number was over the Hox A gene cluster. Genetic markers in significant regions accounted for approximately 50 % of the genomic variation. Less genomic variation for kyphosis was described by QTL regions and no region was associated with kyphosis and vertebra counts.

Conclusions

The importance of the Hox gene families in vertebral development was highlighted as significant associations were detected over the A, B and C families. Further evaluation of these regions and characterization of variants within these genes will expand our knowledge on vertebral development using natural genetic variants segregating in commercial swine.

Electronic supplementary material

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