Lamas-Toranzo 2019 Reduce Mosaicism Bovine Embryos
Lamas-Toranzo 2019 Reduce Mosaicism Bovine Embryos
Lamas-Toranzo 2019 Reduce Mosaicism Bovine Embryos
com/scientificreports
Genome modification at specific loci allows the ablation (knock-out, KO) or insertion (knock-in, KI) of specific
DNA sequences to unequivocally assess the role of a specific gene on a particular physiological process or to alter
the phenotype of an animal for diverse purposes. Unfortunately, technical limitations have largely restricted targeted
genome modification in mammals to the mouse model. These limitations derived from the extremely low efficiency
of the only available technique for targeted mutagenesis, homologous recombination (HR)1, which impedes its direct
application on embryos. The low efficiencies of HR could only be bypassed by the use of embryonic stem cells (ESC)
as an intermediary to generate targeted genetically modified mice2. Later, the development of somatic cell nuclear
transfer (SCNT) allowed the use of somatic cells as intermediaries for HR to generate genetically modified animals
in other species where truly pluripotent ESC were not available3, but the complexity and inefficiency of SCNT has
impeded its widespread use. The advent of site-specific endonucleases has situated genome modification at an attain-
able distance in these species, empowering scientist to perform loss-of-function experiments or to generate geneti-
cally modified animals for different applications4. Among the different site-specific endonucleases, CRISPR, the last
technology to be developed has become the method of choice given its ease of use and flexibility5.
CRISPR technology allows KO generation in a single-step by microinjection at the zygote stage6. For this purpose,
CRISPR is directed to the beginning of the coding region of the target gene, where it will induce a double-strand
break (DSB). The DSB can be repaired by either HR, which reconstitutes the target site allowing CRISPR recogni-
tion and thereby the generation of another DSB, or by non-homologous end joining repair (NHEJ), which often
generates random insertions or deletions (indels) at the target site. Those indels constitute a stable mutation, as they
impede target recognition by CRISPR, and can produce KO alleles, as those indels not multiple of three disrupt the
open reading frame (ORF) of the target gene, leading to a truncated protein. However, KO generation in one step
(i.e. one pregnancy) requires all alleles harboured by a given individual to be KO (i.e., not multiple of three) and,
therefore, a reduction in the number of alleles generated greatly increases the odds for direct KO generation.
Animal Reproduction Department, INIA, Madrid, Spain. Correspondence and requests for materials should be
addressed to P.B.-Á. (email: [email protected])
No. of embryos genotyped No. of embryos No. of embryos genotyped No. of mosaic
Group by PCR sequencing edited (%) by clonal sequencing embryos (%)
RNP-injected 0 hpi 23 20 (87.0) 20 6 (30.0)a
mRNA-injected 0 hpi 25 22 (88.0) 20 6 (30.0)a
mRNA-injected 10 hpi 24 20 (83.3) 20 7 (35.0)a
RNP-injected 20 hpi 6 5 (83.3) 5 5 (100)b
mRNA-injected 20 hpi 25 21 (84.0) 10 10 (100)b
Table 1. Genome edition and mosaicism rates following the alternative protocols tested. Different superscript
letters indicate significant differences based on Chi-square test (p < 0.05).
The possible generation of more than two alleles per individual following CRISPR edition was initially
overlooked, as seminal work in mice did not report this phenomenon6, and because in this species breeding of
founders to obtain an heterozygous F1 generation is the routine protocol for experimental purposes. Under ideal
conditions, zygote genome edition should occur at the 2n2c stage, resulting in two indels (alleles). However,
DNA replication occurs soon after fertilization, before pronuclei fusion, transitioning to the 2n4c stage where
genome edition can result into more than two alleles, a phenomenon called genetic mosaicism. DNA edition may
also occur at later developmental stages (after cleavage) if both unedited alleles and active ribonucleoprotein are
present. Mosaic individuals are composed by more than one genetic cell type and have been consistently found
in most publications that have performed allele screening following CRISPR direct injection in zygotes of diverse
species such as mice7,8, pigs9–19, goats20, sheep21–23, cattle25 and rabbits26–34. In this article, we report strategies
to reduce genetic mosaicism following CRISPR edition of bovine embryos based on early delivery of CRISPR
components.
Results
Development of a shortened IVF protocol. Conventional protocols of bovine in vitro fertilization (IVF)
entail the co-incubation of cumulus-oocyte complexes and spermatozoa for ~20 h35. Following co-incubation,
cumulus cells can be removed and then, the ooplasm of the presumptive zygotes can be clearly visualized and
microinjected. These IVF and micromanipulation conditions have been optimized to attain high developmental
rates, but the delivery of CRISPR components by microinjection at 20 h after the onset of IVF resulted in all edited
embryos being mosaic (Table 1). Aiming to explore the different possibilities of an earlier delivery of CRISPR
components, we first established the minimum gamete co-incubation time to achieve normal developmental rates
in a preliminary experiment. For this aim we reduced IVF time to 8, 9 or 10 h, observing that, in the case of the
semen used in these experiments, 10 h was the minimum gamete co-incubation time to achieve similar develop-
mental rates to the conventional 20 h IVF protocol (Fig. 1A,B and Table S1).
Characterization of S-phase in bovine zygotes. Once the minimum co-incubation time was estab-
lished, we analysed the timing of DNA replication in bovine IVF derived zygotes from the minimum IVF time
(10 hours post-insemination –hpi-) to 20 hpi in 2 h intervals by 5-Ethynyl-2′-deoxyuridine (EdU) incorporation
(Fig. 1C). The data obtained show that DNA replication occurs short after fertilization, with ~40% of the zygotes
already replicating its DNA at 10 hpi (Fig. 1D). Four hours later (14 hpi) most zygotes were on S-phase and at the
time point when gamete co-incubation ceases in conventional IVF procotols (20 hpi), ~40% of the zygotes are on
S-phase. These kinetics suggest that S-phase spans for about 10 h and that by the end of the conventional gamete
co-incubation (20 hpi), most zygotes have replicated its DNA or are very close to end their S-phase. In this con-
text, an earlier CRISPR deliver is required to prevent the appearance of more than two alleles in a given embryo.
Genome edition and mosaicism rates following different microinjection times. Genome edi-
tion rates were assessed in all four microinjected groups (3 microinjected with mRNA + sgRNA at 0, 10 and 20
hpi, and 1 microinjected with ribonucleoprotein at 0 hpi) and in 5 embryos obtained from a replicate microin-
jected with ribonucleoprotein at 20 hpi (Table 1). For this aim, a PCR product including CRISPR target site was
Figure 1. Developmental rates following shortened gamete co-incubation times and kinetics of S-phase in
bovine IVF zygotes. (A) Cleavage and (B) blastocyst rates obtained following gamete co-incubation for 8, 9, 10
or 20 h. (C) Representative pictures of S-phase analysis in bovine zygotes following EdU protocol; upper images
show a zygote that is not replicating its DNA (EdU negative), lower images depict a zygote on S-phase (EdU
positive). (D) Percentages of zygotes replicating its DNA from 10 to 20 hpi in 2 h intervals.
amplified and sequenced in at least 20 blastocysts per group. All groups displayed similar genome edition rates,
with more than 80% of the blastocyst produced from microinjected oocytes or zygotes being edited at the target
site. This sequencing strategy allowed to distinguish between edited and not-edited embryos, but does not pro-
vide the number of alleles harboured by a given embryo, as mixed sequencing peaks preclude allele identification
(Fig. 3). Allele identification was achieved by clonal sequencing, analysing 10 colonies per embryo in 20 embryos
on each of the three early microinjection groups (600 sequences) and 10 and 5 in 20 hpi group microinjected with
mRNA or ribonucleoprotein, respectively (150 sequences) (Fig. S1). Mosaicism rates were significantly reduced
in all early delivery groups compared to 20 hpi (70% reduction from 100% to ~30%, Table 1), but no differences
were noted between 0 or 10 hpi or between CRISPR injection formats. The number of embryos carrying WT (not
edited) alleles and carrying only alleles formed by indels not multiple of three (there would be KO embryos if the
target was allocated on a coding region) are shown in Table 2. As expected the percentage of embryos carrying
WT alleles was significantly higher in the group where all embryos were mosaic (20 hpi).
Genome edition and mosaicism rates targeting multiple coding loci. To further test the most
user-friendly of the two alternative microinjection protocols developed (0 hpi), we analysed genome edition
and mosaicism rates in blastocysts derived from mRNA (n = 26) or ribonucleoprotein microinjections (n = 23)
targeting two coding loci simultaneously. These loci were localized at the beginning of the coding region of two
genes encoding for the proteins CSN2 and PAEP, both present in cow milk and responsible for allergic reactions
in humans. In this experiment genome edition rates were around or above 90%, and all embryos edited on one
locus were also edited in the other except for one in the RNP group (Table 3). Mosaicism rates were reduced to
12–24%, and again no statistically significant differences were found on edition, mosaicism or KO generation
rates between the groups injected with mRNA or ribonucleoprotein (Tables 3 and 4). Embryos were deemed as
KO for a gene when all alleles detected were formed by indels not multiple of three.
Discussion
Mosaicism constitutes a major problem for direct genome modification following the microinjection of CRISPR
microinjection in zygotes. The presence of more than two alleles in a given individual reduces greatly the chances
of KO generation in one step, without the need of subsequent breeding of founder individuals, as all alleles gener-
ated must disrupt the ORF of the target gene. The generation of a KO individual in a single pregnancy is especially
relevant in livestock species, where, in contrast to mice, generation intervals can be counted by years. In this per-
spective, the requirement of breeding mosaic founders for the generation of bi-allelic animals impede the culmi-
nation of genome modification projects within a reasonable time frame in these species. Besides, KO generation
Figure 2. Developmental rates following the alternative microinjection protocols tested. (A) Schematic
representation of the three alternative protocols tested and the conventional 20 hpi protocol. (B) Cleavage and
(C) blastocyst rates obtained following the alternative protocols tested with or without microinjection. Different
letters indicate significant differences based on ANOVA (p < 0.05).
in one step would allow to elucidate gene function during embryo development without the need of a colony of
heterozygous founders. A myriad of publications that have performed allele screening following CRISPR direct
injection in zygotes of diverse species such as mice7,8, pigs9–19, goats20, sheep21–23, cattle25 and rabbits26–34 have
encountered this problem and yet, probably due to the time and resource consuming allele screening, no study
has systemically analysed mosaicism rates following different protocols.
The analysis of S-phase on bovine zygotes show that DNA replicates short after fertilization, providing a phys-
iological explanation for the high mosaicism rates obtained following conventional microinjection times, as four
DNA copies of the target sequence will be present when genome edition occurs. The timing and duration of
S-phase in bovine in vitro produced zygotes overtly agrees with previous studies36,37, and although there are differ-
ences between bulls in both the onset of S-phase37 and cleavage division38, the protocols tested in this study could
be applied to the use of different sires with minor (10 hpi) or no (0 hpi) modifications. Likewise, early delivery
of CRISPR components is expected to reduce mosaicism rates in other species. In this context, earlier CRISPR
delivery has been observed to diminish mosaicism rates following CRISPR electroporation in mice39, and small
data sets from sheep24 and human40 studies also suggest a reduction.
The early microinjection protocols tested achieved a ~70–90% reduction in mosaicism rates. This reduction
constitutes a remarkable advance for direct KO generation, but mosaicism was not completely abolished. The
use of ribonucleoprotein instead of RNAs did not reduced further mosaicism rates, suggesting that the time
required for RNA translation and ribonucleoprotein assembly was not a crucial factor for mosaicism appearance.
Figure 3. PCR and clonal sequencing of an edited non-mosaic embryo. Upper image show the sequencing
reaction of a PCR product, mixed peaks indicates edition but do not allow allele discrimination. Lower images
show the sequencing reaction of individual alleles following clonal sequencing.
No. of edited No. of edited non- No. of embryos containing No. of embryos containing only
Group embryos mosaic embryos (%) an unedited WT allele (%) indels non-multiple of 3 (%)
RNP-injected 0 hpi 20 14 (70.0)a 8 (40.0)a,b 3 (15.0)
mRNA-injected 0 hpi 20 14 (70.0)a 6 (30.0)a 6 (30.0)
mRNA-injected 10 hpi 20 13 (65.0)a 6 (30)a 6 (30.0)
RNP-injected 20 hpi 5 0 (0)n.d. 4 (80.0)n.d. 0 (0)
mRNA-injected 20 hpi 10 0 (0)b 8 (80.0)b 1 (10.0)
Table 2. Percentage of embryos containing unedited WT alleles or only indels non-multiple of 3. Different
superscript letters indicate significant differences based on Chi-square test (p < 0.05) between the four
microinjection groups.
No. of embryos genotyped by No. of embryos No. of mosaic embryos No. of embryos genotyped by No. of embryos No. of mosaic embryos
Group PCR sequencing for PAEP edited for PAEP (%) for PAEP (%) PCR sequencing for CSN2 edited for CSN2 (%) for CSN2 (%)
RNP-injected
23 21 (91.3) 3 (14.3) 23 20 (87.0) 3 (15.0)
0 hpi
mRNA-injected
26 25 (96.2) 3 (12.0) 26 25 (96.2) 6 (24.0)
0 hpi
Table 3. Genome edition and mosaicism rates following dual targeting (PAEP and CSN2) by mRNA or RNP
microinjections at 0 hpi. No significant differences were observed based on Chi-square test (p > 0.05).
Somehow surprisingly, CRISPR delivery to bovine oocytes (i.e., before fertilization) did not decrease mosaicism
rates further compared to the delivery to presumptive zygotes right after a reduced IVF protocol (10 hpi), sug-
gesting that CRISPR activity is very low before pronuclei formation. A possible explanation for this result is that
CRISPR may not be able to recognize its target locus before some degree of chromatin decondensation has been
achieved in the pronuclei, as chromatin accessibility exerts a great impact on Cas9 binding in vivo41,42. In this
sense, the high degree of DNA condensation in matured oocytes and especially on spermatozoa43 may impede an
earlier CRISPR-mediated genome edition.
In conclusion, early delivery of CRISPR components to bovine oocytes prior to IVF or zygotes following a
shortened IVF reduces mosaicism rates from 100% to ~10–30% while achieving similar genome edition and
developmental rates. Oocyte microinjection is more convenient schedule-wise than microinjection following a 10
h IVF protocol and both RNA or ribonucleoprotein delivery formats can be used, as both achieved similar results.
Methods
Bovine in vitro production. Ovaries were collected at local slaughterhouse and transported to the lab-
oratory within 2 h. Cumulus-oocyte complexes (COCs) were collected from 2 to 8 mm diameter follicles and
selected based on conventional morphological criteria44. In vitro maturation (IVM) was performed in TCM-199
supplemented with 10% (v/v) fetal calf serum (FCS) and 10 ng/ml epidermal growth factor at 39 °C under an
atmosphere of 5% CO2 in air with maximum humidity for 24 h. IVF was performed with frozen-thawed sperma-
®
tozoa from a single sire selected by Bovipure (Nidacon). Spermatozoa were diluted to a final concentration of
106 spermatozoa/ml and were co-incubated with mature oocytes in TALP medium at 39 °C under an atmosphere
No. of embryos No. of KO embryos No. of embryos No. of KO embryos No. of embryos edited No. of KO embryos for
Group edited for PAEP for PAEP (%) edited for CSN2 for CSN2 (%) for both PAEP and CSN2 both PAEP and CSN2 (%)
RNP-injected
21 11 (52.3) 20 8 (40.0) 20 5 (25.0)
0 hpi
mRNA-injected 0 hpi 25 8 (32.0) 25 14 (56.0) 25 5 (20.0)
Table 4. Percentage of KO embryos (i.e., harbouring only frame-disrupting alleles) following dual targeting
(PAEP and CSN2) by mRNA or RNP microinjections at 0 hpi.
of 5% CO2 in air with maximum humidity. Different co-incubation times were initially tested (8, 9, 10 and 20 h)
to determine the minimum co-incubation time (>80 embryos/group in two independent replicates). For micro-
injection experiments, gametes were co-incubated for 10 (10 hpi group) or 20 h (0 and 20 hpi groups) (Fig. 2A).
Microinjected and non-microinjected oocytes of 0 hpi groups were denuded by vortexing in PBS supplemented
with 300 µg/ml hyaluronidase for 3 min prior to fertilization, in order to visualize the ooplasm membrane for
microinjection. In 10 and 20 hpi groups intact COCs were used for IVF and cumulus cells were removed at the
end of the co-incubation time by vortexing in PBS for 3 min. Cumulus-free presumptive zygotes were cultured in
vitro in SOF media supplemented with 5% FCS under an atmosphere of 5% CO2 and 5% O2 in air with maximum
humidity. Cleavage rates were assessed at 48 hpi and blastocyst yield was recorded 9 days post-insemination (dpi)
in 3–4 independent replicates per group. Statistical differences in developmental rates were assessed by One Way
ANOVA using SigmaStat software.
Assessment of DNA replication in bovine zygotes. DNA synthesis was detected by the incorporation
of 5-Ethynyl-2′-deoxyuridine (EdU) (Click-iT EdU HCS assay kit, Invitrogen). Briefly, following 10, 12, 14, 18
and 20 hpi, 243 presumptive zygotes from 3 independent replicates were vortexed 3 min to remove cumulus cells
and incubated in 50 µM EdU for 30 min45. EdU labelling was detected following manufacturer instructions and
pronuclei were counterstained with 10 µg/ml Hoechst. Zygotes were observed on an epifluorescence inverted
microscope (Nikon Eclipse TE300) and considered to be on S-phase stage when they displayed at least one pro-
nucleus labelled with EdU.
CRISPR/Cas9 microinjection. Details for sgRNAs are provided on Table S3. sgRNAs were synthesized and
®
purified using Guide-it sgRNA In Vitro Transcription Kit (Takara). Capped polyadenylated Cas9 mRNA was
®
produced by in vitro transcription (mMESSAGE mMACHINE T7 ULTRA kit , Life Technologies) using as tem-
plate the plasmid pMJ920 (Addgene 42234) linearized with BbsI and treated with Antarctic phosphatase (NEB).
mRNA was purified using MEGAClear kit (Life Technologies). A solution of 300 ng/µl of mRNA and 100 ng/µl
of each sgRNA was used for RNA microinjection. For ribonucleoprotein injection, Guide-it Recombinant Cas9
(Takara) and sgRNA/s were mixed to a final concentration of 300 ng/µl and 60 ng/µl, respectively, and incubated
at 37 °C for 5 min to achieve ribonucleoprotein assembly prior to microinjection. Previous experiments were con-
ducted to determine that the concentrations used did not reduce developmental rates compared to sham (buffer)
injections. Microinjection was performed under a Nikon Diaphot TMD inverted microscope delivering 3–5 pl
into the ooplasm using a filament needle.
Embryo genotyping. CRISPR-injected embryos were kept in culture until day 9 after insemination. Zona
pellucida was removed from unhatched blastocysts by incubation in acid PBS (pH 2) to avoid any residual sper-
matozoa and to facilitate subsequent enzymatic digestion. Zona-free blastocysts were individually placed at the
bottom of a 0.2 ml PCR tube and stored at −80 °C until analysis. Each blastocyst was digested in 8 µl of Picopure ®
(ThermoFisher Scientific) for 1 h at 65 °C followed by inactivation at 95 °C for 10 min.
Genotyping of the intronic region was performed by clonal sequencing. PCR was performed on a 50 µl
reaction containing 4 µl of the inactivated digestion product under the following conditions: 95 °C for 2 min;
39 × (95 °C for 20 s, 60 °C for 30 s, 72 °C for 40 s); 72 °C for 5 min; hold at 8 °C. PCR was performed using primers
spanning the target sequence (Table S4). The PCR product from each blastocyst was purified using FavorPrep ™
PCR Purification Kit (Favorgen). The purified product was Sanger sequenced and analyzed for the presence of
indels. Consecutively, the purified PCR products from edited embryos, identified by a mixed sequence reaction
around the target site (Fig. 3), were analyzed by clonal sequencing. For that aim, the purified PCR product was
ligated into pMD20 T-vector (Takara) by Blunt/TA Ligase (NEB) and transformed into Escherichia coli DH5-α
competent cells. For each embryo analysed, ten plasmids containing the insert were Sanger sequenced.
Genotyping of PAEP and CSN2 target sequences was performed by deep sequencing. A first 30 cycle amplicon
PCR was performed on 25 µl reaction containing 2 µl of the inactivated digestion product using the conditions
detailed above and primers including Illumina overhangs (marked in italics in Table S4). Amplicons were purified
by AMPPure XP beads (Beckman Coulter) and libraries were prepared by an index PCR which added Illumina
adaptors and indexes identifying each embryo using Nextera XT (Illumina). Libraries were purified, pooled to
2 nM and sequenced on Illumina miSeq platform providing 250 bp paired-end sequencing reads. Individual
alleles were identified following QC filtering, mapping to reference and variant calling.
Embryos showing more than two alleles by clonal or deep sequencing were considered mosaic. Statistical
differences in genome edition and mosaicism rates were assessed by Chi-square test using SigmaStat software.
References
1. Brinster, R. L. et al. Targeted correction of a major histocompatibility class II E alpha gene by DNA microinjected into mouse eggs.
Proc Natl Acad Sci USA 86, 7087–7091 (1989).
2. Evans, M. J., Bradley, A., Kuehn, M. R. & Robertson, E. J. The ability of EK cells to form chimeras after selection of clones in G418
and some observations on the integration of retroviral vector proviral DNA into EK cells. Cold Spring Harb Symp Quant Biol 50,
685–689 (1985).
3. Schnieke, A. E. et al. Human factor IX transgenic sheep produced by transfer of nuclei from transfected fetal fibroblasts. Science 278,
2130–2133 (1997).
4. Lamas-Toranzo, I. et al. CRISPR is knocking on barn door. Reprod Domest Anim 52(Suppl 4), 39–47 (2017).
5. Jinek, M. et al. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science 337, 816–821 (2012).
6. Wang, H. et al. One-step generation of mice carrying mutations in multiple genes by CRISPR/Cas-mediated genome engineering.
Cell 153, 910–918 (2013).
7. Zhong, H., Chen, Y., Li, Y., Chen, R. & Mardon, G. CRISPR-engineered mosaicism rapidly reveals that loss of Kcnj13 function in
mice mimics human disease phenotypes. Sci Rep 5, 8366 (2015).
8. Noda, T., Oji, A. & Ikawa, M. Genome editing in mouse zygotes and embryonic stem cells by introducing sgRNA/Cas9 expressing
plasmids. Methods Mol Biol 1630, 67–80 (2017).
9. Hai, T., Teng, F., Guo, R., Li, W. & Zhou, Q. One-step generation of knockout pigs by zygote injection of CRISPR/Cas system. Cell
research 24, 372–375 (2014).
10. Sato, M. et al. Direct Injection of CRISPR/Cas9-Related mRNA into Cytoplasm of Parthenogenetically Activated Porcine Oocytes
Causes Frequent Mosaicism for Indel Mutations. Int J Mol Sci 16, 17838–17856 (2015).
11. Wang, Y. et al. Efficient generation of gene-modified pigs via injection of zygote with Cas9/sgRNA. Sci Rep 5, 8256 (2015).
12. Yu, H. H. et al. Porcine Zygote Injection with Cas9/sgRNA Results in DMD-Modified Pig with Muscle Dystrophy. Int J Mol Sci 17,
e1668 (2016).
13. Chuang, C. K. et al. Generation of GGTA1 Mutant Pigs by Direct Pronuclear Microinjection of CRISPR/Cas9 Plasmid Vectors. Anim
Biotechnol 28, 174–181 (2016).
14. Zhou, X. et al. Efficient Generation of Gene-Modified Pigs Harboring Precise Orthologous Human Mutation via CRISPR/Cas9-
Induced Homology-Directed Repair in Zygotes. Hum Mutat 37, 110–118 (2016).
15. Kang, J. T. et al. Biallelic modification of IL2RG leads to severe combined immunodeficiency in pigs. Reprod Biol Endocrinol 14, 74
(2016).
16. Petersen, B. et al. Efficient production of biallelic GGTA1 knockout pigs by cytoplasmic microinjection of CRISPR/Cas9 into
zygotes. Xenotransplantation 23, 338–346 (2016).
17. Park, K. E. et al. Targeted gene knock-in by CRISPR/Cas ribonucleoproteins in porcine zygotes. Sci Rep 7, 42458 (2017).
18. Burkard, C. et al. Precision engineering for PRRSV resistance in pigs: Macrophages from genome edited pigs lacking CD163 SRCR5
domain are fully resistant to both PRRSV genotypes while maintaining biological function. PLoS Pathog 13, e1006206 (2017).
19. Whitworth, K. M. et al. Zygote injection of CRISPR/Cas9 RNA successfully modifies the target gene without delaying blastocyst
development or altering the sex ratio in pigs. Transgenic Res 26, 97–107 (2017).
20. Wang, X. et al. Disruption of FGF5 in Cashmere Goats Using CRISPR/Cas9 Results in More Secondary Hair Follicles and Longer
Fibers. PLoS One 11, e0164640 (2016).
21. Crispo, M. et al. Efficient Generation of Myostatin Knock-Out Sheep Using CRISPR/Cas9 Technology and Microinjection into
Zygotes. PLoS One 10, e0136690 (2015).
22. Wang, X. et al. Multiplex gene editing via CRISPR/Cas9 exhibits desirable muscle hypertrophy without detectable off-target effects
in sheep. Sci Rep 6, 32271 (2016).
23. Zhang, X. et al. Disruption of the sheep BMPR-IB gene by CRISPR/Cas9 in in vitro-produced embryos. Theriogenology 91, 163–172
e162 (2017).
24. Vilarino, M. et al. CRISPR/Cas9 microinjection in oocytes disables pancreas development in sheep. Sci Rep 7, 17472 (2017).
25. Bevacqua, R. J. et al. Efficient edition of the bovine PRNP prion gene in somatic cells and IVF embryos using the CRISPR/Cas9
system. Theriogenology 86, 1886–1896 (2016).
26. Yan, Q. et al. Generation of multi-gene knockout rabbits using the Cas9/gRNA system. Cell Regen 3, 12 (2014).
27. Honda, A. et al. Single-step generation of rabbits carrying a targeted allele of the tyrosinase gene using CRISPR/Cas9. Exp Anim 64,
31–37 (2015).
28. Yuan, L. et al. CRISPR/Cas9-mediated GJA8 knockout in rabbits recapitulates human congenital cataracts. Sci Rep 6, 22024 (2016).
29. Song, J. et al. RS-1 enhances CRISPR/Cas9- and TALEN-mediated knock-in efficiency. Nat Commun 7, 10548 (2016).
30. Guo, R. et al. Generation and evaluation of Myostatin knock-out rabbits and goats using CRISPR/Cas9 system. Sci Rep 6, 29855
(2016).
31. Song, Y. et al. Efficient dual sgRNA-directed large gene deletion in rabbit with CRISPR/Cas9 system. Cell Mol Life Sci 73, 2959–2968
(2016).
32. Sui, T. et al. CRISPR/Cas9-mediated mutation of PHEX in rabbit recapitulates human X-linked hypophosphatemia (XLH). Hum
Mol Genet 25, 2661–2671 (2016).
33. Lv, Q. et al. Efficient Generation of Myostatin Gene Mutated Rabbit by CRISPR/Cas9. Sci Rep 6, 25029 (2016).
34. Yang, D. et al. Identification and characterization of rabbit ROSA26 for gene knock-in and stable reporter gene expression. Sci Rep
6, 25161 (2016).
35. Parrish, J. J. et al. Bovine in vitro fertilization with frozen-thawed semen. Theriogenology 25, 591–600 (1986).
36. Eid, L. N., Lorton, S. P. & Parrish, J. J. Paternal influence on S-phase in the first cell cycle of the bovine embryo. Biol Reprod 51,
1232–1237 (1994).
37. Comizzoli, P., Marquant-Le Guienne, B., Heyman, Y. & Renard, J. P. Onset of the first S-phase is determined by a paternal effect
during the G1-phase in bovine zygotes. Biol Reprod 62, 1677–1684 (2000).
38. Bermejo-Alvarez, P., Lonergan, P., Rath, D., Gutierrez-Adan, A. & Rizos, D. Developmental kinetics and gene expression in male and
female bovine embryos produced in vitro with sex-sorted spermatozoa. Reprod Fertil Dev 22, 426–436 (2010).
39. Hashimoto, M., Yamashita, Y. & Takemoto, T. Electroporation of Cas9 protein/sgRNA into early pronuclear zygotes generates non-
mosaic mutants in the mouse. Dev Biol 418, 1–9 (2016).
40. Ma, H. et al. Correction of a pathogenic gene mutation in human embryos. Nature 548, 413–419 (2017).
41. Wu, X. et al. Genome-wide binding of the CRISPR endonuclease Cas9 in mammalian cells. Nature Biotechnol 32, 670–676 (2014).
42. Kuscu, C., Arslan, S., Singh, R., Thorpe, J. & Adli, M. Genome-wide analysis reveals characteristics of off-target sites bound by the
Cas9 endonuclease. Nature Biotechnol 32, 677–683 (2014).
43. Kimmins, S. & Sassone-Corsi, P. Chromatin remodelling and epigenetic features of germ cells. Nature 434, 583–589 (2005).
44. Hawk, H. W. & Wall, R. J. Improved yields of bovine blastocysts from in vitro-produced oocytes. I. Selection of oocytes and zygotes.
Theriogenology 41, 1571–1583 (1994).
45. Lin, C. J., Koh, F. M., Wong, P., Conti, M. & Ramalho-Santos, M. Hira-mediated H3.3 incorporation is required for DNA replication
and ribosomal RNA transcription in the mouse zygote. Dev Cell 30, 268–279 (2014).
Acknowledgements
We thank the slaughterhouse “Transformación Ganadera de Leganés SA” for gently providing bovine ovaries to
conduct the experiments. This work has been funded by the projects RYC-2012-10193, AGL2014-58739-R and
AGL2017-84908-R from the Spanish Ministry of Economy and Competitiveness (MINECO) and the project
StG-757886-ELONGAN from the European Research Council. I.L.T. is funded by a FPI fellowship by MINECO.
Author Contributions
I.L.T., B.G.C., F.C.A., J.C.F., O.O. and P.B.A. participated on in vitro embryo production and embryo genotyping.
I.L.T. performed oocyte/zygote microinjections. P.B.A. and I.L.T. performed S-phase analysis, statistical analysis
and wrote the manuscript.
Additional Information
Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-51366-8.
Competing Interests: The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and
institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International
License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.