Analysis of copy number variations in the sheep genome using 50K SNP BeadChip array

Copy number variations affecting growth curve parameters in a crossbred chicken population

Copy number variations (CNVs) are key structural variations in the genome and may contribute to phenotypic differences. In this study, we used a F2 chicken population created from reciprocal crossing between fast-growing Arian broiler line and Urmia native chickens. The chickens were genotyped by 60 K SNP BeadChip, and PennCNV algorithm was used to detect genome-wide CNVs. The growth curve parameters of W0, k, L, Wf, Wi, ti and average GR were used as phenotypic data. The association between CNV and growth curve parameters was carried out using the CNVRanger R/Bioconductor package. Five CNV regions (CNVRs) were chosen for the validation experiment using qPCR. Gene enrichment analysis was done using WebGestalt. The STRING database was used to search for significant pathways. The results identified 966 CNVs and 600 CNVRs including 468 gains, 67 losses, and 65 both events on autosomal chromosomes. Validation of the CNVRs obtained from the qPCR assay were 79 % consistent with the prediction by PennCNV. A total of 43 significant CNVs were obtained for the seven growth curve parameters. The 416 genes annotated for significant CNVs. Six genes out of 416 genes were most related to growth curve parameters. These genes were LCP2, Dock2, CD80, CYFIP1, NIPA1 and NIPA2. Some of these genes in their biological process were associated with the growth, reproduction and development of cells or organs that ultimately lead to the growth of the body. The results of the study could pave the way for better understanding the molecular process of CNVs and growth curve parameters in birds.

Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle

Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.

Genome-wide association study between copy number variation and feeding behavior, feed efficiency, and growth traits in Nellore cattle

BACKGROUND

Feeding costs represent the largest expenditures in beef production. Therefore, the animal efficiency in converting feed in high-quality protein for human consumption plays a major role in the environmental impact of the beef industry and in the beef producers' profitability. In this context, breeding animals for improved feed efficiency through genomic selection has been considered as a strategic practice in modern breeding programs around the world. Copy number variation (CNV) is a less-studied source of genetic variation that can contribute to phenotypic variability in complex traits. In this context, this study aimed to: (1) identify CNV and CNV regions (CNVRs) in the genome of Nellore cattle (Bos taurus indicus); (2) assess potential associations between the identified CNVR and weaning weight (W210), body weight measured at the time of selection (WSel), average daily gain (ADG), dry matter intake (DMI), residual feed intake (RFI), time spent at the feed bunk (TF), and frequency of visits to the feed bunk (FF); and, (3) perform functional enrichment analyses of the significant CNVR identified for each of the traits evaluated.

RESULTS

A total of 3,161 CNVs and 561 CNVRs ranging from 4,973 bp to 3,215,394 bp were identified. The CNVRs covered up to 99,221,894 bp (3.99%) of the Nellore autosomal genome. Seventeen CNVR were significantly associated with dry matter intake and feeding frequency (number of daily visits to the feed bunk). The functional annotation of the associated CNVRs revealed important candidate genes related to metabolism that may be associated with the phenotypic expression of the evaluated traits. Furthermore, Gene Ontology (GO) analyses revealed 19 enrichment processes associated with FF.

CONCLUSIONS

A total of 3,161 CNVs and 561 CNVRs were identified and characterized in a Nellore cattle population. Various CNVRs were significantly associated with DMI and FF, indicating that CNVs play an important role in key biological pathways and in the phenotypic expression of feeding behavior and growth traits in Nellore cattle.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing

Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.

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

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

Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds

Copy number variants (CNVs) are a type of genetic polymorphism which contribute to phenotypic variation in several species, including livestock. In this study, we used genomic data of 192 animals from 3 Iranian sheep breeds including 96 Baluchi sheep and 47 Lori-Bakhtiari sheep as fat-tailed breeds and 47 Zel sheep as thin-tailed sheep breed genotyped with Illumina OvineSNP50K Beadchip arrays. Also, for association test, 70 samples of Valle del Belice sheep were added to the association test as thin-tailed sheep breed. PennCNV and CNVRuler software were, respectively, used to study the copy number variation and genomic association analyses. We detected 573 and 242 CNVs in the fat and thin tailed breeds, respectively. In terms of CNV regions (CNVRs), these represented 328 and 187 CNVRs that were within or overlapping with 790 known Ovine genes. The CNVRs covered approximately 73.85 Mb of the sheep genome with average length 146.88 kb, and corresponded to 2.6% of the autosomal genome sequence. Five CNVRs were randomly chosen for validation, of which 4 were experimentally confirmed using Real time qPCR. Functional enrichment analysis showed that genes harbouring CNVs in thin-tailed sheep were involved in the adaptive immune response, regulation of reactive oxygen species biosynthetic process and response to starvation. In fat-tailed breeds these genes were involved in cellular protein modification process, regulation of heart rate, intestinal absorption, olfactory receptor activity and ATP binding. Association test identified one copy gained CNVR on chromosomes 6 harbouring two protein-coding genes HGFAC and LRPAP1. Our findings provide information about genomic structural changes and their association to the interested traits including fat deposition and environmental compatibility in sheep.

Genomic integration to identify molecular biomarkers associated with indicator traits of gastrointestinal nematode resistance in sheep

This study aimed to integrate GWAS and structural variants to propose possible molecular biomarkers related to gastrointestinal nematode resistance traits in Santa Inês sheep. The phenotypic records FAMACHA, haematocrit, white blood cell count, red blood cell count, haemoglobin, platelets and egg counts per gram of faeces were collected from 700 naturally infected animals, belonging to four Brazilian flocks. A total of 576 animals were genotyped using the Ovine SNP12k BeadChip and were imputed using a reference population with Ovine SNP50 BeadChip. The GWAS approaches were based on SNPs, haplotypes, CNVs and ROH. The overlapping between the significant genomic regions detected from all approaches was investigated, and the results were integrated using a network analysis. Genes related to the immune system were found, such as ABCB1, IL6, WNT5A and IRF5. Genomic regions containing candidate genes and metabolic pathways involved in immune responses, inflammatory processes and immune cells affecting parasite resistance traits were identified. The genomic regions, biological processes and candidate genes uncovered could lead to biomarkers for selecting more resilient sheep and improving herd welfare and productivity. The results obtained are the start point to identify molecular biomarkers related to indicator traits of gastrointestinal nematode resistance in Santa Inês sheep.

Breeding of African sheep reared under low-input/output smallholder production systems for trypanotolerance

Trypanosomiasis is a disease caused by unicellular protozoan parasites. Small ruminants succumb to trypanosomiasis in areas of high tsetse fly challenge, resulting in serious economic loss often to farmers in low-input smallholder systems. At present, trypanosomiasis is treated with trypanocidal drugs, but access to these can be limited, and increasing parasite resistance raises questions about their efficacy. The development of trypanotolerance in small ruminant flocks through targeted breeding strategies is considered a sustainable and economical option for controlling African trypanosomiasis. Recently, quantitative trait loci (QTLs) associated with trypanotolerance traits in sheep have been reported. The results of these studies form the basis for more studies to identify QTLs associated with trypanosomiasis resistance, particularly in African livestock species. For example, signatures of positive selection for trypanotolerance have been identified using genome-wide single-nucleotide polymorphism data. However, there are several challenges in performing genetic analyses using data from low-input smallholder systems, including a lack of recorded pedigree and production records and the need for large sample sizes when flock sizes are often fewer than 50 animals. Breeding strategies to improve trypanotolerance should also preserve existing genetic diversity as well as minimize excessive genetic introgression by trypanosusceptible breeds. This review discusses the possibilities of breeding for trypanosome tolerance/resistance in low-input/low-output small ruminant production systems. Potential challenges are outlined, and potential available genetic resources are described as a foundation for future work.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Identification of Copy Number Variations and Genetic Diversity in Italian Insular Sheep Breeds

Copy number variants (CNVs) are one of the major contributors to genetic diversity and phenotypic variation in livestock. The aim of this work is to identify CNVs and perform, for the first time, a CNV-based population genetics analysis with five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs with an average length of 1.81 Mb. The breeds showed similar mean numbers of CNVs, ranging from 20 (Sarda) to 27 (Comisana). A total of 365 CNV regions (CNVRs) were determined. The length of the CNVRs varied among breeds from 2.4 Mb to 124.1 Mb. The highest number of shared CNVRs was between Comisana and Pinzirita, and only one CNVR was shared among all breeds. Our results indicated that segregating CNVs expresses a certain degree of diversity across all breeds. Despite the low/moderate genetic differentiation among breeds, the different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups, similar to the previous studies based on single-nucleotide polymorphism markers. Gene enrichment was described for the 37 CNVRs selected, considering the top 10%. Out of 181 total genes, 67 were uncharacterized loci. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. Our results corroborated previous studies and showed that CNVs represent valuable molecular resources for providing useful information for separating the population and could be further used to explore the function and evolutionary aspect of sheep genome.

Two Different Copy Number Variations of the CLCN2 Gene in Chinese Cattle and Their Association with Growth Traits

Copy number variation (CNV) can affect gene function and even individual phenotypic traits by changing the transcription and translation level of related genes, and it also plays an important role in species evolution. Chloride voltage-gated channel 2 (CLCN2) encodes a voltage-gated chloride channel (CLC-2), which has a wide organ distribution and is ubiquitously expressed. Based on previous studies, we hypothesize that CLCN2 could be a candidate gene involved in cell volume regulation, transepithelial transport and cell proliferation. This study aimed to explore CNVs in the CLCN2 gene and investigate its association with growth traits in four Chinese cattle breeds (Yunling cattle, Xianan cattle, Qinchuan cattle and Pinan cattle). We identified there are two copy number variation regions (CNV1: 3600 bp, including exon 2-11; CNV2: 4800 bp, including exon 21-22) of the CLCN2 gene. The statistical analysis showed that the CNV1 mutation in the YL cattle population was significantly associated with cannon circumference (p < 0.01). The CNV2 mutation in the XN cattle population had a significant effect on body slanting length, chest girth and body weight (p < 0.05). In the YL cattle, the association analysis of CLCN2 gene CNV1 and CNV2 combination with cannon circumference was significant (p < 0.01). Our results provide evidence that CNV1 and CNV2 in CLCN2 are associated with growth traits in two different cattle populations and could be used as candidate markers for cattle molecular breeding.

Copy number variants in genomes of local sheep breeds from Russia

Copy number variants (CNVs) are genomic structural variations that contribute to many adaptive and economically important traits in livestock. In this study, we detected CNVs in 354 animals from 16 Russian indigenous sheep breeds and analysed their possible functional roles. Our analysis of the entire sample set resulted in 4527 CNVs forming 1450 CNV regions (CNVRs). When constructing CNVRs for individual breeds, a total of 2715 regions ranging from 88 in Groznensk to 337 in Osetin breeds were identified. To make interbreed CNVR frequency comparison possible, we also identified core CNVRs using CNVs with overlapping chromosomal locations found in different breeds. This resulted in 137 interbreed CNVRs with frequency >15% in at least one breed. Functional enrichment analysis of genes affected by CNVRs in individual breeds revealed 12 breeds with significant enrichments in olfactory perception, PRAME family proteins, and immune response. Function of genes affected by interbreed and breed-specific CNVRs revealed candidates related to domestication, adaptation to high altitudes and cold climates, reproduction, parasite resistance, milk and meat qualities, wool traits, fat storage, and fat metabolism. Our work is the first attempt to uncover and characterise the CNV makeup of Russian indigenous sheep breeds. Further experimental and functional validation of CNVRs would help in developing new and improving existing sheep breeds.

Copy number variations of the KAT6A gene are associated with body measurements of Chinese sheep breeds

Copy number variation (CNV) is one kind of genomic structure variations and presents as gains and losses of genomic fragments. More recently, we have made an atlas of CNV maps for livestock. In the future, it is a primary focus to determine the phenotypic effects of candidate CNVs. Lysine Acetyltransferase 6 A (KAT6A) is a protein coding gene and plays a critical role in many cellular processes. However, the effects of KAT6A CNVs on sheep body measurements remains unknown. In this study, we performed quantitative polymerase chain reaction (qPCR) to detect the presences and distributions of three CNV regions within KAT6A gene in 672 sheep from four Chinese breeds. Association analysis indicated that the three CNVs of KAT6A gene were significantly associated with body measurement(s) in Small-tailed Han sheep (STH) and Hu sheep (HU) (p < 0.05), while no effects on Large-tailed Han sheep (LTH) were observed (p > 0.05) were observed. Additionally, only one CNV was significantly associated with body measurement (body length) in Chaka sheep (CK) (p < 0.05). Our study provided evidence that the CNV(s) of KAT6A gene could be used as candidate marker(s) for molecular breedings of STH, HU, and CK breeds.

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 Copy Number Variations and Their Association With Distinct Phenotypes in the World's Sheep

Copy number variations (CNVs) are a major source of structural variation in mammalian genomes. Here, we characterized the genome-wide CNV in 2059 sheep from 67 populations all over the world using the Ovine Infinium HD (600K) SNP BeadChip. We tested their associations with distinct phenotypic traits by conducting multiple independent genome-wide tests. In total, we detected 7547 unique CNVs and 18,152 CNV events in 1217 non-redundant CNV regions (CNVRs), covering 245 Mb (∼10%) of the whole sheep genome. We identified seven CNVRs with frequencies correlating to geographical origins and 107 CNVRs overlapping 53 known quantitative trait loci (QTLs). Gene ontology and pathway enrichment analyses of CNV-overlapping genes revealed their common involvement in energy metabolism, endocrine regulation, nervous system development, cell proliferation, immune, and reproduction. For the phenotypic traits, we detected significantly associated (adjusted P < 0.05) CNVRs harboring functional candidate genes, such as SBNO2 for polycerate; PPP1R11 and GABBR1 for tail weight; AKT1 for supernumerary nipple; CSRP1, WNT7B, HMX1, and FGFR3 for ear size; and NOS3 and FILIP1 in Wadi sheep; SNRPD3, KHDRBS2, and SDCCAG3 in Hu sheep; NOS3, BMP1, and SLC19A1 in Icelandic; CDK2 in Finnsheep; MICA in Romanov; and REEP4 in Texel sheep for litter size. These CNVs and associated genes are important markers for molecular breeding of sheep and other livestock species.

Functional characterization of Copy Number Variations regions in Djallonké sheep

A total of 184 Djallonké (West African Dwarf) sheep of Burkina Faso were analysed for Copy Number Variations (CNV) using Ovine 50 K SNP BeadChip genotyping data and two different CNV calling platforms: PennCNV and QuantiSNP. Analyses allowed to identify a total of 63 candidate Copy Number Variations Regions (CNVR) on 11 different ovine chromosomes covering about 82.5 Mb of the sheep genome. Gene-annotation enrichment analysis allowed to identify a total of 751 potential candidate ovine genes located in the candidate CNVR bounds. Functional annotation allowed to identify five statistically significant Functional Clusters (FC; enrichment factor > 1.3) involving 61 candidate genes. All genes forming significantly enriched FC were located on ovine chromosome (OAR) 21. FC1 (22 genes including PAG4 and PAG6) and FC5 (three genes: CTSC, CTSW and CTSF), coding proteases (peptidases and cathepsins, respectively), were involved in reproductive performance and modulation of gestation. Both FC3 and FC4 were involved in inflammatory and immunologic response through coding serum amyloid A and B-box-type zinc finger proteins, respectively. Finally, FC2 consisted of 27 genes (including OR10G6 and OR8B8) involved in olfactory receptor activity, key for animals adapting to new food resources. CNVR identified on at least 15% of individuals were considered CNVR hotspots and further overlapped with previously reported quantitative trait loci (QTL). CNVR hotspots spanning genes putatively involved with lipid metabolism (SKP1, TCF7, JADE2, UBE2B and SAR1B) and differential expression in mammary gland (SEC24A and CDKN2AIPNL) on OAR5 and dairy traits (CCDC198 and SLC35F4) on OAR7 overlapped with QTL associated with lipid metabolism, milk protein yield and milk fat percentage. Information obtained from local sheep populations naturally adapted to harsh environments contributes to increase our understanding of the genomic importance of CNV.

A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing

Background

Copy number variation (CNV) is an important source of genetic variation that has a significant influence on phenotypic diversity, economically important traits and the evolution of livestock species. In this study, the genome-wide CNV distribution characteristics of 32 fine-wool sheep from three breeds were analyzed using resequencing.

Results

A total of 1,747,604 CNVs were detected in this study, and 7228 CNV regions (CNVR) were obtained after merging overlapping CNVs; these regions accounted for 2.17% of the sheep reference genome. The average length of the CNVRs was 4307.17 bp. “Deletion” events took place more frequently than “duplication” or “both” events. The CNVRs obtained overlapped with previously reported sheep CNVRs to variable extents (4.39–55.46%). Functional enrichment analysis showed that the CNVR-harboring genes were mainly involved in sensory perception systems, nutrient metabolism processes, and growth and development processes. Furthermore, 1855 of the CNVRs were associated with 166 quantitative trait loci (QTL), including milk QTLs, carcass QTLs, and health-related QTLs, among others. In addition, the 32 fine-wool sheep were divided into horned and polled groups to analyze for the selective sweep of CNVRs, and it was found that the relaxin family peptide receptor 2 (RXFP2) gene was strongly influenced by selection.

Conclusions

In summary, we constructed a genomic CNV map for Chinese indigenous fine-wool sheep using resequencing, thereby providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07387-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.

Copy number variation of the HPGDS gene in the Ashidan yak and its associations with growth traits

Copy number variation (CNV) is a structural variation at the submicroscopic level of the genome, which can affect gene-related phenotypes by changing genes dosage and transcript structure. Hematopoietic prostaglandin D synthase (HPGDS) is a member whose functions are closely related to weight gain and inflammatory diseases of the glutathione S-transferase (GSTs) family. In this study, the growth characteristics (body weight, withers height, body length, and chest girth) of 336 Ashidan yaks were monitored at four stages (6 months, 12 months, 18 months, and 30 months). In addition, CNV of the HPGDS gene was detected, discovered relationships of CNV with growth traits, and explored the level of gene expression. Based on the statistical analysis by IBM SPSS software, significant correlations were observed between HPGDS-CNV and body weight in 12-month-old yak (P < 0.01), 18-month-old yak (P < 0.001) and 30-month-old yak (P < 0.001) and body length in 18-month-old yak (P < 0.05) and 30-month-old yak (P < 0.05), respectively. Additionally, the individuals with gain copy number type performed better in body weight and body length than those with normal or loss copy number type. To our best of knowledge, this is the first time to make efforts to probe into the role of HPGDS-CNV and its interaction with livestock growth traits. Our results suggested that the CNV of the HPGDS gene may be an active candidate gene for the marker-assisted selection (MAS) of yaks.

Integrating Genome-Wide CNVs Into QTLs and High Confidence GWAScore Regions Identified Positional Candidates for Sheep Economic Traits

Copy number variations (CNVs) are important source of genetic variation, which can affect diverse economic traits through a variety of mechanisms. In addition, genome scan can identify many quantitative trait loci (QTLs) for the economic traits, while genome-wide association studies (GWAS) can localize genetic variants associated with the phenotypic variations. Here, we developed a method called GWAScore which collected GWAS summary data to identify potential candidates, and integrated CNVs into QTLs and high confidence GWAScore regions to detect crucial CNV markers for sheep growth traits. We got 197 candidate genes which were overlapping with the candidate CNVs. Some crucial genes (MYLK3, TTC29, HERC6, ABCG2, RUNX1, etc.) showed significantly elevated GWAScore peaks than other candidate genes. In this study, we developed the GWAScore method to excavate the potential value of candidate genes as markers for the sheep molecular breeding.

Genomic Structural Diversity in Local Goats: Analysis of Copy-Number Variations

Simple Summary

Copy-number variations (CNVs) are one of the widely dispersed forms of structural variations in mammalian genomes and are known to be present in genomic regions that regulate important physiological functions. In this study, CNV detection was performed starting from genotypic data of 120 individuals, belonging to four Sicilian dairy goat breeds, genotyped with the Illumina GoatSNP50 BeadChip array. Using PennCNV software, a total of 702 CNVs were identified in 107 individuals. These were merged in 75 CNV regions (CNVRs), i.e., regions containing CNVs overlapped by at least 1 base pair. Functional annotation of the CNVRs allowed the identification of 139 genes/loci within the most frequent CNVRs, which are involved in local adaptation, mild behaviour, immune response, reproduction, and olfactory receptors. This study provides insights into the genomic variations within these Italian goat breeds and should be of value for future studies to identify the relationships between this type of genetic variation and phenotypic traits.

Abstract

Copy-number variations (CNVs) are one of the widely dispersed forms of structural variations in mammalian genomes, and are present as deletions, insertions, or duplications. Only few studies have been conducted in goats on CNVs derived from SNP array data, and many local breeds still remain uncharacterized, e.g., the Sicilian goat dairy breeds. In this study, CNV detection was performed, starting from the genotypic data of 120 individuals, belonging to four local breeds (Argentata dell’Etna, Derivata di Siria, Girgentana, and Messinese), genotyped with the Illumina GoatSNP50 BeadChip array. Overall, 702 CNVs were identified in 107 individuals using PennCNV software based on the hidden Markov model algorithm. These were merged in 75 CNV regions (CNVRs), i.e., regions containing CNVs overlapped by at least 1 base pair, while 85 CNVs remained unique. The part of the genome covered by CNV events was 35.21 Mb (1.2% of the goat genome length). Functional annotation of the CNVRs allowed the identification of 139 genes/loci within the most frequent CNVRs that are involved in local adaptations, such as coat colour (ADAMTS20 and EDNRA), mild behaviour (NR3C2), immune response (EXOC3L4 and TNFAIP2), reproduction (GBP1 and GBP6), and olfactory receptors (OR7E24). This study provides insights into the genomic variations for these Sicilian dairy goat breeds and should be of value for future studies to identify the relationships between this type of genetic variation and phenotypic traits.

Genome-Wide Detection of CNVs and Association With Body Weight in Sheep Based on 600K SNP Arrays

Copy number variations (CNVs) are important genomic structural variations and can give rise to significant phenotypic diversity. Herein, we used high-density 600K SNP arrays to detect CNVs in two synthetic lines of sheep (DS and SHH) and in Hu sheep (a local Chinese breed). A total of 919 CNV regions (CNVRs) were detected with a total length of 48.17 Mb, accounting for 1.96% of the sheep genome. These CNVRs consisted of 730 gains, 102 losses, and 87 complex CNVRs. These CNVRs were significantly enriched in the segmental duplication (SD) region. A CNVR-based cluster analysis of the three breeds revealed that the DS and SHH breeds share a close genetic relationship. Functional analysis revealed that some genes in these CNVRs were also significantly enriched in the olfactory transduction pathway (oas04740), including members of the OR gene family such as OR6C76, OR4Q2, and OR4K14. Using association analyses and previous gene annotations, we determined that a subset of identified genes was likely to be associated with body weight, including FOXF2, MAPK12, MAP3K11, STRBP, and C14orf132. Together, these results offer valuable information that will guide future efforts to explore the genetic basis for body weight in sheep.

Copy Number Variation of the PIGY Gene in Sheep and Its Association Analysis with Growth Traits

Simple Summary

The PIGY (phosphatidylinositol glycan anchor biosynthesis class Y) gene is a member of the PIG gene family and encodes the glycosylphosphatidylinositol-N-acetylglucosaminyltransferase (GPI-GnT) complex. It initiates the biosynthesis of GPI and plays an important role in cell–cell interactions. Sequencing has revealed a 3600 bp copy number variation (CNV) in exon 2 of the PIGY gene in sheep, potentially altering a functional part of the protein. The CNV overlaps 28 quantitative trait loci that are relevant to some economic traits like muscle density and carcass weight. We screened for this CNV of the PIGY gene in 569 individuals, namely, 240 Chaka sheep (CKS), 168 Hu sheep (HS), and 161 small-tailed Han sheep (STHS), and analyzed the association between the presence of this CNV and sheep body size traits. The results showed that the loss-type CNV was more prevalent than other types in these three breeds, and there were significant effects of the PIGY gene CNV on body weight, chest circumference, and circumference of cannon bone of sheep. The results showed that sheep with gain-type CNV had better growth traits than those with other types. The findings reveal the relationship between the CNV of the PIGY gene and growth traits of sheep, suggesting that CNV could be utilized for improved molecular breeding of sheep.

Abstract

Copy number variation (CNV) is a type of genomic variation with an important effect on animal phenotype. We found that the PIGY gene contains a 3600 bp copy number variation (CNV) region located in chromosome 6 of sheep (Oar_v4.0 36,121,601–36,125,200 bp). This region overlaps with multiple quantitative trait loci related to phenotypes like muscle density and carcass weight. Therefore, in this study, the copy number variation of the PIGY gene was screened in three Chinese sheep breeds, namely, Chaka sheep (CKS, May of 2018, Wulan County, Qinghai Province, China), Hu sheep (HS, May of 2015, Mengjin County, Henan Province, China), and small-tailed Han sheep (STHS, May of 2016, Yongjing, Gansu Province, China). Association analyses were performed on the presence of CNV and sheep body size traits. We used real-time quantitative PCR (qPCR) to detect the CNV for association analysis. According to the results, the loss-type CNV was more common than other types in the three breeds (global average: loss = 61.5%, normal = 17.5%, and gain = 21.0%). The association analysis also showed significant effects of the PIGY gene CNV on body weight, chest circumference, and circumference of the cannon bone of sheep. Sheep with gain-type CNV had better growth traits than those with other types. The results indicate a clear relationship between the PIGY gene CNV and growth traits of sheep, suggesting the use of CNV as a new molecular breeding marker.

Novel copy number variation of the BAG4 gene is associated with growth traits in three Chinese sheep populations

Copy number variation (CNV) as an important source of genetic phenotypic and variation is related to complex phenotypic traits. The aim of this study was to investigate the potential associations of BAG4 (Bcl-2-associated athanogene 4) copy numbers variations with sheep growth traits in three Chinese sheep breeds (CKS, STHS, and HS). BAG4 is located within the stature and udder attachment quantitative trait loci (QTL) in sheep. Expression profiling revealed that the BAG4 gene was widely expressed in the tissues of sheep. The distribution of BAG4 gene copy number showed that the loss of copy number was more dominant in CKS and HS which was different from that in STHS. Statistical analysis revealed that the BAG4 CNV was significantly associated with body height in CKS (p < 0.05), with body slanting length in HS (p < 0.05), and with body height and hip cross height in STHS (p < 0.05). The χ² values showed significant differences in the BAG4 CNV distribution frequency between varieties. In conclusion, the results establish the association between BAG4 CNV and sheep traits and suggest that BAG4 CNV may be a promising marker for the molecular breeding of Chinese sheep.

Sustained Postnatal Skin Regeneration Upon Prenatal Application of Functionalized Collagen Scaffolds

Primary closure of fetal skin in spina bifida protects the spinal cord and improves clinical outcome, but is also associated with postnatal growth malformations and spinal cord tethering. In this study, we evaluated the postnatal effects of prenatally closed full-thickness skin defects in sheep applying collagen scaffolds with and without heparin/vascular endothelial growth factor/fibroblast growth factor 2, focusing on skin regeneration and growth. At 6 months, collagen scaffold functionalized with heparin, VEGF, and FGF2 (COL-HEP/GF) resulted in a 6.9-fold increase of the surface area of the regenerated skin opposed to 1.7 × for collagen only. Epidermal thickness increased 5.7-fold at 1 month, in line with high gene expression of S100 proteins, and decreased to 2.1 at 6 months. Increased adipose tissue and reduced scaffold degradation and number of myofibroblasts were observed for COL-HEP/GF. Gene ontology terms related to extracellular matrix (ECM) organization were enriched for both scaffold treatments. In COL-HEP/GF, ECM gene expression resembled native skin. Expression of hair follicle-related genes in COL-HEP/GF was comparable to native skin, and de novo hair follicle generation was indicated. In conclusion, in utero closure of skin defects using functionalized collagen scaffolds resulted in long-term skin regeneration and growth. Functionalized collagen scaffolds that grow with the child may be useful for prenatal treatment of closure defects like spina bifida. Impact statement Prenatal closure of fetal skin in case of spina bifida prevents damage to the spinal cord. Closure of the defect is challenging and may result in postnatal growth malformations. In this study, the postnatal effects of a prenatally applied collagen scaffold functionalized with heparin and vascular endothelial growth factor (VEGF)/fibroblast growth factor (FGF) were investigated. An increase of the surface area of regenerated skin ("growing with the child") and generation of hair follicles was observed. Gene expression levels resembled those of native skin with respect to the extracellular matrix and hair follicles. Overall, in utero closure of skin defects using heparin/VEGF/FGF functionalized collagen scaffolds results in long-term skin regeneration.

Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak

Simple Summary

Cell adhesion molecule 2 (CADM2), also known as synaptic cell adhesion molecule 2 (SYNCAM2), is the mediator of synaptic signals enriched in the brain. Overlaps between copy number variation (CNV) regions in CADM2 and quantitative trait loci (QTL) related to body weight have been clarified in a previous study. In this study, two loci were amplified in the CADM2 gene (CNV1: 235,915 bp, exon 1 and partial intron 1; CNV2: 60,430 bp, intron 9) to explore the relationship between CNV types in the CADM2 gene and growth traits in 350 Ashidan yaks. Association analysis illustrated that no significant effect was found on growth traits in CNV1. However, the CNV2 mutation had a significant effect on body weight at the sixth month (p < 0.05). Individuals with the gain-type copy number variation CNV2 were significantly superior to those with loss- or normal-type in terms of body weight (p < 0.05). In summary, this study confirmed that CADM2-CNVs affect growth traits in yaks, and may be candidate genes for successful yak breeding and genetics projects.

Abstract

Copy number variation (CNV) is currently accepted as a common source of genetic variation. It is reported that CNVs may influence the resistance to disease and complex economic traits, such as residual feed intake, muscle formation, and fat deposition in livestock. Cell adhesion molecule 2 (CADM2) is expressed widely in the brain and adipose tissue and can regulate body weight through the central nervous system. Growth traits are important economic traits for animal selection. In this study, we aimed to explore the effect of CADM2 gene copy number variants on yak growth traits. Here, two CNVs in the CADM2 gene were investigated using the quantitative polymerase chain reaction (qPCR), and the association of the CNVs with growth traits in yak was analyzed using statistical methods by SPSS software. Differences were considered significant if the p value was < 0.05. Statistical analysis indicated significant association of CADM2-CNV2 with the body weight of the Chinese Ashidan yak. A significant effect of CNV2 (p < 0.05) was found on body weight at 6 months. In CNV2, the gain-type copy number variation exhibited greater performance than the other variants, with greater body weight observed at 6 months (p < 0.05). To the best of our knowledge, this is the first attempt to investigate the function of CADM2-CNVs and their association with growth traits in animals. This may be a useful candidate marker in marker-assisted selection of yaks.

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.

Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip

Background

Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks.

Results

A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits.

Conclusions

Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.

A genome-wide single nucleotide polymorphism and copy number variation analysis for number of piglets born alive

Background

In this study we integrated the CNV (copy number variation) and WssGWAS (weighted single-step approach for genome-wide association) analyses to increase the knowledge about number of piglets born alive, an economically important reproductive trait with significant impact on production efficiency of pigs.

Results

A total of 3892 samples were genotyped with the Porcine SNP80 BeadChip. After quality control, a total of 57,962 high-quality SNPs from 3520 Duroc pigs were retained. The PennCNV algorithm identified 46,118 CNVs, which were aggregated by overlapping in 425 CNV regions (CNVRs) ranging from 2.5 Kb to 9718.4 Kb and covering 197 Mb (~ 7.01%) of the pig autosomal genome. The WssGWAS identified 16 genomic regions explaining more than 1% of the additive genetic variance for number of piglets born alive. The overlap between CNVR and WssGWAS analyses identified common regions on SSC2 (4.2–5.2 Mb), SSC3 (3.9–4.9 Mb), SSC12 (56.6–57.6 Mb), and SSC17 (17.3–18.3 Mb). Those regions are known for harboring important causative variants for pig reproductive traits based on their crucial functions in fertilization, development of gametes and embryos. Functional analysis by the Panther software identified 13 gene ontology biological processes significantly represented in this study such as reproduction, developmental process, cellular component organization or biogenesis, and immune system process, which plays relevant roles in swine reproductive traits.

Conclusion

Our research helps to improve the understanding of the genetic architecture of number of piglets born alive, given that the combination of GWAS and CNV analyses allows for a more efficient identification of the genomic regions and biological processes associated with this trait in Duroc pigs. Pig breeding programs could potentially benefit from a more accurate discovery of important genomic regions.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5687-0) contains supplementary material, which is available to authorized users.

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.

Connecting genome structural variation with complex traits in crop plants

Structural genome variation is a major determinant of useful trait diversity. We describe how genome analysis methods are enabling discovery of trait-associated structural variants and their potential impact on breeding. As our understanding of complex crop genomes continues to grow, there is growing evidence that structural genome variation plays a major role in determining traits important for breeding and agriculture. Identifying the extent and impact of structural variants in crop genomes is becoming increasingly feasible with ongoing advances in the sophistication of genome sequencing technologies, particularly as it becomes easier to generate accurate long sequence reads on a genome-wide scale. In this article, we discuss the origins of structural genome variation in crops from ancient and recent genome duplication and polyploidization events and review high-throughput methods to assay such variants in crop populations in order to find associations with phenotypic traits. There is increasing evidence from such studies that gene presence-absence and copy number variation resulting from segmental chromosome exchanges may be at the heart of adaptive variation of crops to counter abiotic and biotic stress factors. We present examples from major crops that demonstrate the potential of pangenomic diversity as a key resource for future plant breeding for resilience and sustainability.

Detection of copy number variations in brown and white layers based on genotyping panels with different densities

Background

Copy number variations (CNV) are an important source of genetic variation that has gained increasing attention over the last couple of years. In this study, we performed CNV detection and functional analysis for 18,719 individuals from four pure lines and one commercial cross of layer chickens. Samples were genotyped on four single nucleotide polymorphism (SNP) genotyping platforms, i.e. the Illumina 42K, Affymetrix 600K, and two different customized Affymetrix 50K chips. CNV recovered from the Affymetrix chips were identified by using the Axiom^® CNV Summary Tools and PennCNV software and those from the Illumina chip were identified by using the cnvPartition in the Genome Studio software.

Results

The mean number of CNV per individual varied from 0.50 to 4.87 according to line or cross and size of the SNP genotyping set. The length of the detected CNV across all datasets ranged from 1.2 kb to 3.2 Mb. The number of duplications exceeded the number of deletions for most lines. Between the lines, there were considerable differences in the number of detected CNV and their distribution. Most of the detected CNV had a low frequency, but 19 CNV were identified with a frequency higher than 5% in birds that were genotyped on the 600K panel, with the most common CNV being detected in 734 birds from three lines.

Conclusions

Commonly used SNP genotyping platforms can be used to detect segregating CNV in chicken layer lines. The sample sizes for this study enabled a detailed characterization of the CNV landscape within commercially relevant lines. The size of the SNP panel used affected detection efficiency, with more CNV detected per individual on the higher density 600K panel. In spite of the high level of inter-individual diversity and a large number of CNV observed within individuals, we were able to detect 19 frequent CNV, of which, 57.9% overlapped with annotated genes and 89% overlapped with known quantitative trait loci.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0428-4) contains supplementary material, which is available to authorized users.

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.

Copy number variation in livestock: A mini review

Copy number variation (CNV) is a phenomenon in which sections of the genome, ranging from one kilo base pair (Kb) to several million base pairs (Mb), are repeated and the number of repeats vary between the individuals in a population. It is an important source of genetic variation in an individual which is now being utilized rather than single nucleotide polymorphisms (SNPs), as it covers the more genomic region. CNVs alter the gene expression and change the phenotype of an individual due to deletion and duplication of genes in the copy number variation regions (CNVRs). Earlier, researchers extensively utilized SNPs as the main source of genetic variation. But now, the focus is on identification of CNVs associated with complex traits. With the recent advances and reduction in the cost of sequencing, arrays are developed for genotyping which cover the maximum number of SNPs at a time that can be used for detection of CNVRs and underlying quantitative trait loci (QTL) for the complex traits to accelerate genetic improvement. CNV studies are also being carried out to understand the evolutionary mechanism in the domestication of livestock and their adaptation to the different environmental conditions. The main aim of the study is to review the available data on CNV and its role in genetic variation among the livestock.

Mapping of genome-wide copy number variations in the Iranian indigenous cattle using a dense SNP data set

The objective of this study was to present the first map of the copy number variations (CNVs) in Iranian indigenous cattle based on a high-density single nucleotide polymorphism (SNP) dataset. A total of 90 individuals were genotyped using the Illumina BovineHD BeadChip containing 777 962 SNPs. The QuantiSNP algorithm was used to perform a genome-wide CNV detection across autosomal genome. After merging the overlapping CNV, a total of 221 CNV regions were identified encompassing 36.4 Mb or 1.44% of the bovine autosomal genome. The length of the CNV regions ranged from 3.5 to 2252.8 Kb with an average of 163.8 Kb. These regions included 147 loss (66.52%) and 74 gain (33.48%) events containing a total of 637 annotated Ensembl genes. Gene ontology analysis revealed that most of genes in the CNV regions were involved in environmental responses, disease susceptibility and immune system functions. Furthermore, 543 of these genes corresponded to the human orthologous genes, which involved in a wide range of biological functions. Altogether, 73% of the 221 CNV regions overlapped either completely or partially with those previously reported in other cattle studies. Moreover, novel CNV regions involved several quantitative trait loci (QTL)-related to adaptative traits of Iranian indigenous cattle. These results provided a basis to conduct future studies on association between CNV regions and phenotypic variations in the Iranian indigenous cattle.

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

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

Detection of breed-specific copy number variations in domestic chicken genome

Copy number variations (CNVs) are important large-scale variants. They are widespread in the genome and may contribute to phenotypic variation. Detection and characterization of CNVs can provide new insights into the genetic basis of important traits. Here, we perform whole-genome short read sequence analysis to identify CNVs in two indigenous and commercial chicken breeds to evaluate the impact of the identified CNVs on breed-specific traits. After filtration, a total of 12 955 CNVs spanning (on average) about 9.42% of the chicken genome were found that made up 5467 CNV regions (CNVRs). Chicken quantitative trait loci (QTL) datasets and Ensembl gene annotations were used as resources for the estimation of potential phenotypic effects of our CNVRs on breed-specific traits. In total, 34% of our detected CNVRs were also detected in earlier CNV studies. These CNVRs partly overlap several previously reported QTL and gene ontology terms associated with some important traits, including shank length QTL in Creeper-specific CNVRs and body weight and egg production characteristics, as well as muscle and body organ growth, in the Arian commercial breed. Our findings provide new insights into the genomic structure of the chicken genome for an improved understanding of the potential roles of CNVRs in differentiating between breeds or lines.

Genomic regions and pathways associated with gastrointestinal parasites resistance in Santa Inês breed adapted to tropical climate

Background

The aim of this study was to estimate variance components and to identify genomic regions and pathways associated with resistance to gastrointestinal parasites, particularly Haemonchus contortus, in a breed of sheep adapted to tropical climate. Phenotypes evaluations were performed to verify resistance to gastrointestinal parasites, and were divided into two categories: i) farm phenotypes, assessing body condition score (BCS), degree of anemia assessed by the famacha chart (FAM), fur score (FS) and feces consistency (FC); and ii) lab phenotypes, comprising blood analyses for hematocrit (HCT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin (HGB), platelets (PLT) and transformed (log₁₀) egg per gram of feces (EPG_log). A total of 576 animals were genotyped with the Ovine SNP12k BeadChip (Illumina, Inc.), that contains 12,785 bialleleic SNP markers. The variance components were estimated using a single trait model by single step genomic BLUP procedure.

Results

The overall linkage disequilibrium (LD) mean between pairs of markers measured by r² was 0.23. The overall LD mean between markers considering windows up to 10 Mb was 0.07. The mean LD between adjacent SNPs across autosomes ranged from 0.02 to 0.10. Heritability estimates were low for EPG_log (0.11), moderate for RBC (0.18), PLT (0.17) HCT (0.20), HGB (0.16) and WBC (0.22), and high for FAM (0.35). A total of 22, 21, 23, 20, 26, 25 and 23 windows for EPG_log for FAM, WBC, RBC, PLT, HCT and HGB traits were identified, respectively. Among the associated windows, 10 were shown to be common to HCT and HGB traits on OAR1, OAR2, OAR3, OAR5, OAR8 and OAR15.

Conclusion

The traits indicating gastrointestinal parasites resistance presented an adequate genetic variability to respond to selection in Santa Inês breed, and it is expected a higher genetic gain for FAM trait when compared to the others. The level of LD estimated for markers separated by less than 1 Mb indicated that the Ovine SNP12k BeadChip might be a suitable tool for identifying genomic regions associated with traits related to gastrointestinal parasite resistance. Several candidate genes related to immune system development and activation, inflammatory response, regulation of lymphocytes and leukocytes proliferation were found. These genes may help in the selection of animals with higher resistance to parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s40104-017-0190-4) contains supplementary material, which is available to authorized users.

Distribution and Functionality of Copy Number Variation across European Cattle Populations

Copy number variation (CNV), which is characterized by large-scale losses or gains of DNA fragments, contributes significantly to genetic and phenotypic variation. Assessing CNV across different European cattle populations might reveal genetic changes responsible for phenotypic differences, which have accumulated throughout the domestication history of cattle as consequences of evolutionary forces that act upon them. To explore pattern of CNVs across European cattle, we genotyped 149 individuals, that represent different European regions, using the Illumina Bovine HD Genotyping array. A total of 9,944 autosomal CNVs were identified in 149 samples using a Hidden Markov Model (HMM) as employed in PennCNV. Animals originating from several breeds of British Isles, and Balkan and Italian regions, on average, displayed higher abundance of CNV counts than Dutch or Alpine animals. A total of 923 CNV regions (CNVRs) were identified by aggregating CNVs overlapping in at least two animals. The hierarchical clustering of CNVRs indicated low differentiation and sharing of high-frequency CNVRs between European cattle populations. Various CNVRs identified in the present study overlapped with olfactory receptor genes and genes related to immune system. In addition, we also detected a CNV overlapping the Kit gene in English longhorn cattle which has previously been associated with color-sidedness. To conclude, we provide a comprehensive overview of CNV distribution in genome of European cattle. Our results indicate an important role of purifying selection and genomic drift in shaping CNV diversity that exists between different European cattle populations.

Haplotypes for Type, Degree, and Rate of Marbling in Cattle Are Syntenic with Human Muscular Dystrophy

Traditional analyses of a QTL on Bota 19 implicate a surfeit of candidates, but each is of marginal significance in explaining the deposition of healthy, low melting temperature fat within marbled muscle of Wagyu cattle. As an alternative approach, we have used genomic, multigenerational segregation to identify 14 conserved, ancestral 20 Mb haplotypes. These determine the degree and rate of marbling in Wagyu and other breeds of cattle. The melting temperature of intramuscular fat is highly heritable and traceable by haplotyping. Fortunately, for the production of healthy beef, some of these haplotypes are sufficiently penetrant to be expressed in heterozygous crossbreds, thereby allowing selection of sires which will improve the healthiness of beef produced under even harsh climatic conditions. The region of Bota 19 is syntenic to a region of Hosa 17 known to be important in muscle metabolism and in determining susceptibility to a form of human muscular dystrophy.

Genome-wide detection of copy number variation in Chinese indigenous sheep using an ovine high-density 600 K SNP array

Copy number variants (CNVs) represent a form of genomic structural variation underlying phenotypic diversity. In this study, we used the Illumina Ovine SNP 600 K BeadChip array for genome-wide detection of CNVs in 48 Chinese Tan sheep. A total of 1,296 CNV regions (CNVRs), ranging from 1.2 kb to 2.3 Mb in length, were detected, representing approximately 4.7% of the entire ovine genome (Oar_v3.1). We combined our findings with five existing CNVR reports to generate a composite genome-wide dataset of 4,321 CNVRs, which revealed 556 (43%) novel CNVRs. Subsequently, ten novel CNVRs were randomly chosen for further quantitative real-time PCR (qPCR) confirmation, and eight were successfully validated. Gene functional enrichment revealed that these CNVRs cluster into Gene Ontology (GO) categories of homeobox and embryonic skeletal system morphogenesis. One CNVR overlapping with the homeobox transcription factor DLX3 and previously shown to be associated with curly hair in sheep was identified as the candidate CNV for the special curly fleece phenotype in Tan sheep. We constructed a Chinese indigenous sheep genomic CNV map based on the Illumina Ovine SNP 600 K BeadChip array, providing an important addition to published sheep CNVs, which will be helpful for future investigations of the genomic structural variations underlying traits of interest in sheep.

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.