Functional and population genetic features of copy number variations in two dairy cattle populations
Background Copy Number Variations (CNVs) are gain or loss of DNA segments that are known to play a role in shaping a wide range of phenotypes. In this study, we used two dairy cattle populations, Holstein Friesian and Jersey, to discover CNVs using the Illumina BovineHD Genotyping BeadChip aligned to the ARS-UCD1.2 assembly. The discovered CNVs were investigated for their functional impact and their population genetics features. Results We discovered 14,272 autosomal CNVs, which were aggregated into 1,755 CNV regions (CNVR) from 451 animals. These CNVRs together cover 2.8% of the bovine autosomes. The assessment of the functional impact of CNVRs showed that rare CNVRs (MAF < 0.01) are more likely to overlap with genes, than common CNVRs (MAF ≥ 0.05). The Population differentiation index (Fst) based on CNVRs revealed multiple highly diverged CNVRs between the two breeds. Some of these CNVRs overlapped with candidate genes such as MGAM and ADAMTS17 genes , which are related to starch digestion and body size, respectively. Lastly, linkage disequilibrium (LD) between CNVRs and BovineHD BeadChip SNPs was generally low, close to 0, although common deletions (MAF ≥ 0.05) showed slightly higher LD ( r 2 =~0.1 at 10kb distance) than the rest. Nevertheless, this LD is still lower than SNP-SNP LD ( r 2 =~0.5 at 10kb distance). Conclusions Our analyses showed that CNVRs detected using BovineHD BeadChip arrays are likely to be functional. This finding indicates that CNVs can potentially disrupt the function of genes and thus might alter phenotypes. Also, the population differentiation index revealed two candidate genes, MGAM and ADAMTS17 , which hint at adaptive evolution between the two populations. Lastly, low CNVR-SNP LD implies that genetic variation from CNVs might not be fully captured in routine animal genetic evaluation, which relies solely on SNP markers. Keywords : Copy number variations, Bos taurus , Linkage disequilibrium, population genetics
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
This is a list of supplementary files associated with this preprint. Click to download.
Posted 15 Jan, 2020
On 28 Jan, 2020
On 14 Jan, 2020
On 14 Jan, 2020
On 08 Jan, 2020
On 08 Jan, 2020
On 11 Dec, 2019
On 10 Dec, 2019
On 09 Dec, 2019
On 09 Dec, 2019
On 06 Dec, 2019
On 21 Nov, 2019
On 20 Nov, 2019
On 20 Nov, 2019
Received 22 Oct, 2019
On 22 Oct, 2019
Received 11 Oct, 2019
On 30 Sep, 2019
On 23 Sep, 2019
Invitations sent on 22 Sep, 2019
On 03 Sep, 2019
On 09 Aug, 2019
On 08 Aug, 2019
On 07 Aug, 2019
Functional and population genetic features of copy number variations in two dairy cattle populations
Posted 15 Jan, 2020
On 28 Jan, 2020
On 14 Jan, 2020
On 14 Jan, 2020
On 08 Jan, 2020
On 08 Jan, 2020
On 11 Dec, 2019
On 10 Dec, 2019
On 09 Dec, 2019
On 09 Dec, 2019
On 06 Dec, 2019
On 21 Nov, 2019
On 20 Nov, 2019
On 20 Nov, 2019
Received 22 Oct, 2019
On 22 Oct, 2019
Received 11 Oct, 2019
On 30 Sep, 2019
On 23 Sep, 2019
Invitations sent on 22 Sep, 2019
On 03 Sep, 2019
On 09 Aug, 2019
On 08 Aug, 2019
On 07 Aug, 2019
Background Copy Number Variations (CNVs) are gain or loss of DNA segments that are known to play a role in shaping a wide range of phenotypes. In this study, we used two dairy cattle populations, Holstein Friesian and Jersey, to discover CNVs using the Illumina BovineHD Genotyping BeadChip aligned to the ARS-UCD1.2 assembly. The discovered CNVs were investigated for their functional impact and their population genetics features. Results We discovered 14,272 autosomal CNVs, which were aggregated into 1,755 CNV regions (CNVR) from 451 animals. These CNVRs together cover 2.8% of the bovine autosomes. The assessment of the functional impact of CNVRs showed that rare CNVRs (MAF < 0.01) are more likely to overlap with genes, than common CNVRs (MAF ≥ 0.05). The Population differentiation index (Fst) based on CNVRs revealed multiple highly diverged CNVRs between the two breeds. Some of these CNVRs overlapped with candidate genes such as MGAM and ADAMTS17 genes , which are related to starch digestion and body size, respectively. Lastly, linkage disequilibrium (LD) between CNVRs and BovineHD BeadChip SNPs was generally low, close to 0, although common deletions (MAF ≥ 0.05) showed slightly higher LD ( r 2 =~0.1 at 10kb distance) than the rest. Nevertheless, this LD is still lower than SNP-SNP LD ( r 2 =~0.5 at 10kb distance). Conclusions Our analyses showed that CNVRs detected using BovineHD BeadChip arrays are likely to be functional. This finding indicates that CNVs can potentially disrupt the function of genes and thus might alter phenotypes. Also, the population differentiation index revealed two candidate genes, MGAM and ADAMTS17 , which hint at adaptive evolution between the two populations. Lastly, low CNVR-SNP LD implies that genetic variation from CNVs might not be fully captured in routine animal genetic evaluation, which relies solely on SNP markers. Keywords : Copy number variations, Bos taurus , Linkage disequilibrium, population genetics
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5