Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle.

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype

QTL-traits that overlap with polymorphic CNVRs in the Nelore population.

<p>Depicted are the most frequently overlapped traits. Red bars correspond to the observed number of polymorphic CNVRs overlapping with the respective trait, and cyan bars indicate the mean overlap when sampling 1000x random regions matching the polymorphic CNVRs in size and chromosomal location. The error bars indicate the standard deviation and permutation <i>p-</i>values are listed on the right.</p

Functional groups associated with polymorphic CNVRs (>1% of the population).

<p>The statistical significance level was set to 0.01. The enrichment <i>p</i>-values listed in the third column were computed using DAVID [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0157711#pone.0157711.ref062" target="_blank">62</a>].</p

Overlap of the detected Nelore CNVRs with known cattle CNVRs.

<p><b>A)</b> The Venn diagram shows the genomic size in Mb that is overlapped. Known CNVRs listed in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0157711#pone.0157711.s009" target="_blank">S3 Table</a> were reduced to unique and non-overlapping CNVRs. The overlapped fraction of 79 Mb corresponds to 1,387 Nelore CNVRs. <b>(B)</b> For these 1,387 CNVRs, the histogram shows the number of CNVRs (<i>y</i>-axis) overlapping with known CNVRs by the percentage shown on the <i>x-</i>axis. For example, genomic locations of >700 Nelore CNVRs overlap individually >90% with genomic locations of known CNVRs.</p

Shear force estimated breeding value (SF-EBV) distribution in each CN state for six significant CNV segments.

<p>Each dot represents an animal in the corresponding CN state (0-3n) on the <i>x</i>-axis and the assigned SF-EBV on the <i>y</i>-axis. The legend on the top right displays the color code for the CN state. See the main text for a detailed description of each segment and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0157711#pone.0157711.s017" target="_blank">S11 Table</a> for additional details.</p

Chromosomal distribution of 521 polymorphic CNVRs (>1% of the population).

<p>The regions were categorized into loss or gain of copy number and states for which both events were observed.</p