GWAS conducted directly on imputed whole genome sequence have led to the identification of numerous genetic variants associated with agronomic traits in cattle. However, such variants are often simply markers in linkage disequilibrium with the actual causal variants, which is a limiting factor for the development of accurate genomic predictions. It is possible to identify causal variants by integrating information on how variants impact gene expression into GWAS output. RNA splicing plays a major role in regulating gene expression. Thus, assessing the effect of variants on RNA splicing may explain their function. Here, we use a high-throughput strategy to functionally analyse putative splice-disrupting variants in the bovine genome. Using GWAS, massively parallel reporter assay and deep learning algorithms designed to predict splice-disrupting variants, we identify 38 splice-disrupting variants associated with complex traits in cattle, three of which could be classified as causal. Our results indicate that splice-disrupting variants are widely found in the quantitative trait loci related to these phenotypes. Using our combined approach, we also assess the validity of splicing predictors originally developed to analyse human variants in the context of the bovine genome.
Full author list & citation
Mathieu Charles, Nicolas Gaiani, Marie-Pierre Sanchez, Mekki Boussaha, Chris Hozé, Didier Boichard, Dominique Rocha, Arnaud Boulling. Functional impact of splicing variants in the elaboration of complex traits in cattle. 2025-04-24. https://doi.org/10.1038/s41467-025-58970-5
The authors transfected the allele-resolved Vex-seq plasmid library into HEK293T cells in three biological replicates. The library contained positive controls, bovine GWAS candidates (Var.GWAS), and prediction-enriched bovine variants (Var.P); exon inclusion was quantified as PSI and the alternate-minus-reference change as ΔPSI.
The authors transfected the allele-resolved Vex-seq plasmid library into immortalized bovine mammary epithelial MAC-T cells in three biological replicates. The library contained positive controls, bovine GWAS candidates (Var.GWAS), and prediction-enriched bovine variants (Var.P); exon inclusion was quantified as PSI and the alternate-minus-reference change as ΔPSI.
Original supplemental and deposited inputs retained for this study. Download files individually or together as a ZIP; nested folders are preserved. Source reuse terms apply, and sequencing reads may be omitted.