96 Validation of the candidate mutation responsible for embryonic lethality in Holstein haplotype 2 carriers

Abstract

Holstein haplotype 2 (HH2) is embryonic lethal and carried by 1.21% of the US Holstein population. Using next-generation sequencing, we identified a high-impact frameshift mutation in intraflagellar protein 80 (IFT80) as the putative causal mutation. In bovine embryos, IFT80 expression begins at the 8-cell stage and decreases by the blastocyst stage. We hypothesised that the loss of function of IFT80 early in development causes the lethal phenotype. The aim of this study was to mimic the mutation observed in vivo using a CRISPR-Cas9 approach to determine its effect on embryo development. Two guide RNAs (gRNAs) were designed to disrupt exon 11 (Ex11), one before and one after the known IFT80 mutation site, creating a 317-nucleotide (nt) cut to facilitate genotyping. Then, gRNAs annealed to a tracr-Cas9mRNA complex were delivered to 1-cell embryos by microinjection. Each replicate contained control embryos injected with only Cas9mRNA and treated embryos injected with gRNAs targeting IFT80. Embryos from each group were collected at the 8-cell stage for genotyping and gene expression analysis (n = 47), or on Day 8 to validate genotypes of embryos left to develop (n = 50). DNA sequences containing gRNA target sequences were amplified and visualised on an agarose gel. IFT80 expression was determined in biallelic embryos (n = 13) using quantitative PCR and normalized to GAPDH. Primers were designed for the transcript regions before and after gRNAs target sequences, exons 9 and 12, respectively. Expression data were analysed using SAS software (v. 9.4; SAS Institute Inc.) using PROC GLM and LSMEANS to determine expression differences. Biallelic samples (n = 9) were Sanger-sequenced (SS) and aligned with the reference sequence to determine exact cut sites. Protein amino acid (AA) sequences were predicted using SS data. Protein models were constructed using the I-Tasser platform, and then aligned and visualised using PyMol 2.4. Biallelic edits showed a significant decrease in exon 12 expression (P < 0.05), and no difference in exon 9 compared with controls (P > 0.05), indicating that the transcript was severely affected downstream of the edited sites. The reference protein model contained 777 AA, whereas the biallelic sample with the most accurate cut sites yielded a 385-AA protein, indicating that the mutation severely altered protein conformation and possible function. Embryos injected with CRISPR-Cas9 targeting Ex11 arrested at the 8-cell stage and failed to form blastocysts. Day 8 embryos were genotyped (n = 24) and 58% were biallelic, 21% were monoallelic, and 21% appeared wild-type. Given the high rate of edits, the observed embryonic arrest is likely due to disruption of IFT80, and wild-type embryos may contain small edits not visible by gel. In conclusion, generation of CRISPR-Cas9 IFT80 knockouts demonstrated that the frameshift mutation in Ex11 results in a seemingly nonfunctional protein that is responsible for the embryonic lethality seen in HH2 carriers. Future research is needed to determine how IFT80 regulates embryonic development.

This research was supported by USDA-NIFA National Needs Fellowship, USDA-NIFA AFRI Grant No. 2019-67015-28998.