Mapel, X. M.; Leonard, A. S.; Pausch, H.

Sequencing mapping cohorts with long-read technology is crucial to understand the impact of structural variants (SVs) on complex traits. Here, we obtained 4.86 terabases of HiFi reads with an average read N50 of 16.3 Kb from 120 Bos taurus taurus bulls, yielding a mean coverage depth of 13.5-fold. We genotyped 23.8 M small variants and 79.3 k SVs to perform association testing with molecular phenotypes derived from a subset of 117 bulls with total RNA sequencing data from testis tissue. We identified 27.3 k molecular QTL including 316 for which SVs were the top variants. This corresponds to a 2.1- and 5.6-fold enrichment of SVs among expression and splicing QTL, respectively. When considering SVs in perfect LD with the lead small variant, the enrichment increased to 6.1- and 12-fold for expression and splicing QTL, respectively. Imperfect genotyping for large SVs and other variants limited our ability to detect all SV top variants, suggesting that the true enrichment of SVs among molecular QTL may be even higher. These results demonstrate that SVs have profound impacts on gene expression and splicing variation in cattle but highlight the necessity of improved SV genotyping to fully leverage long-read sequencing cohorts for dissecting complex traits.