Identification of nonsense variants in the genomes of 15 Murciano-Granadina bucks and analysis of their segregation in parent-offspring trios

Nonsense variants can inactivate gene function by causing the synthesis of truncated proteins or by inducing nonsense mediated decay of messenger RNAs. The occurrence of such variants in the genomes of livestock species is modulated by multiple demographic and selective factors. Even though nonsense variants can have causal effects on embryo lethality, abortions, and disease, their genomic distribution and segregation in domestic goats have not been characterized in depth yet. In this work, we have sequenced the genomes of 15 Murciano-Granadina bucks with an average coverage of 32.92× ± 1.45×. Bioinformatic analysis revealed 947 nonsense variants consistently detected with SnpEff and Ensembl-VEP. These variants were especially abundant in the 3' end of the protein-coding regions. Genes related to olfactory perception, ATPase activity coupled to transmembrane movement of substances, defense to virus, hormonal response, and sensory perception of taste were particularly enriched in nonsense variants. Seventeen nonsense variants expected to have harmful effects on fitness were genotyped in parent-offspring trios. We observed that several nonsense variants predicted to be lethal based on mouse knockout data did not have such effect, a finding that could be explained by the existence of multiple mechanisms counteracting lethality. These findings demonstrate that predicting the effects of putative nonsense variants on fitness is extremely challenging. As a matter of fact, such a goal could only be achieved by generating a high-quality telomere-to-telomere goat reference genome combined with carefully curated annotation and functional testing of promising candidate variants.