- 1Acceligen Inc., Eagan, MN, United States
- 2Department of Preventive Veterinary Medicine and Animal Reproduction, School of Agricultural and Veterinarian Sciences, São Paulo State University (Unesp), Jaboticabal, Brazil
- 3Recombinetics Inc., Eagan, MN, United States
- 4USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, United States
- 5Hypor, Hendrix Genetics, Boxmeer, Netherlands
Introduction: Most male pigs are surgically castrated to avoid puberty-derived boar taint and aggressiveness. However, this surgical intervention represents a welfare concern in swine production. Disrupting porcine KISS1 is hypothesized to delay or abolish puberty by inducing variable hypogonadotropism and thus preventing the need for castration.
Methods: To test this hypothesis, we generated the first KISS1-edited large animal using CRISPR/Cas9-ribonucleoproteins and single-stranded donor oligonucleotides. The targeted region preceded the sequence encoding a conserved core motif of kisspeptin. Genome editors were intracytoplasmically injected into 684 swine zygotes and transferred to 19 hormonally synchronized surrogate sows. In nine litters, 49 American Yorkshire and 20 Duroc liveborn piglets were naturally farrowed.
Results: Thirty-five of these pigs bore KISS1-disruptive alleles ranging in frequency from 5% to 97% and did not phenotypically differ from their wild-type counterparts. In contrast, four KISS1-edited pigs (two boars and two gilts) with disruptive allele frequencies of 96% and 100% demonstrated full hypogonadotropism, infantile reproductive tracts, and failed to reach sexual maturity. Change in body weight during development was unaffected by editing KISS1. Founder pigs partially carrying KISS1-disruptive alleles were bred resulting in a total of 53 KISS1+/+, 60 KISS1+/−, and 34 KISS1−/− F1 liveborn piglets, confirming germline transmission.
Discussion: Results demonstrate that a high proportion of KISS1 alleles in pigs must be disrupted before variation in gonadotropin secretion is observed, suggesting that even a small amount of kisspeptin ligand is sufficient to confer proper sexual development and puberty in pigs. Follow-on studies will evaluate fertility restoration in KISS1 KO breeding stock to fully realize the potential of KISS1 gene edits to eliminate the need for surgical castration.
1 Introduction
Surgical castration is routinely implemented in male piglets destined for pork production to prevent the development of boar taint and reduce androgen-driven behaviors (e.g., aggression and mounting) that increase the risk of injuries (Rault et al., 2011). Boar taint is a an off-odor and strong flavor found in meat from intact male pigs that consumers find unacceptable. Surgical castration is a welfare concern because it is considered painful and methods for effective analgesia are largely unavailable (Rault et al., 2011). Some alternatives to this management procedure are commercially available, such as immunization against gonadotropin-releasing hormone (GnRH) (Dunshea et al., 2001), but there are a number of constraints that limit their use (Bonneau and Weiler, 2019; Squires et al., 2020); and the goal of eliminating castration has not proven feasible (Backus et al., 2018).
We sought a genetic mechanism to avoid the need for surgical castration. Editing a gene to block sexual maturation was a promising approach as male pigs that remain prepubertal are not expected to develop boar taint and aggressive behavior. We selected the kisspeptin system based on its conserved role in initiating mammalian puberty (Lents, 2019; Uenoyama et al., 2019; Sobrino et al., 2022). Kisspeptin is a peptide encoded by the highly conserved KISS1 gene, that stimulates the release of GnRH and secretion of gonadotropins (Lents, 2019). Mutations in the kisspeptin receptor gene (KISS1R) result in hypogonadotropic hypogonadism (HH) and insufficient sexual maturity in humans (de Roux et al., 2003; Seminara et al., 2003; Semple et al., 2005). Similarly, knocking out either Kiss1 (d’Anglemont et al., 2007; Lapatto et al., 2007; Uenoyama et al., 2015; Ikegami et al., 2020) or Kiss1r (Funes et al., 2003; Seminara et al., 2003; Lapatto et al., 2007) genes in laboratory rodents results in pubertal failure and infertility owing to HH. Hypogonadotropic hypogonadism was also induced in pigs by Sonstegard et al. (2016), Sonstegard et al. (2017), who developed the first large animal model of impaired kisspeptin system when they knocked out KISS1R in pigs using TALENs (Tan et al., 2013), demonstrating that kisspeptin signaling is vital for sexual maturation of boars.
Humans, mice, and pigs with impaired KISS1/KISS1R genes have responded to exogenous GnRH, gonadotropins (Seminara et al., 2003; Sonstegard et al., 2017), or kisspeptin analogs (d’Anglemont et al., 2007; Lapatto et al., 2007), although some of these approaches only partially reversed the KISS1R KO phenotype in boars (Sonstegard et al., 2017). We hypothesized that KISS1 KO pigs would be phenocopies of KISS1R KO pigs, as observed in some Kiss1 KO mice (Lapatto et al., 2007). However, editing KISS1 should have the advantage of making fertility rescue less challenging because the ligand can be given exogenously. Therefore, KISS1 KO pigs were generated to test the hypothesis that the kisspeptin ligand is essential for sexual maturity in swine. This was accomplished by microinjecting CRISPR/Cas9-ribonucleoproteins (RNPs) and homology-directed repair (HDR) templates into in vivo fertilized porcine zygotes, to create mosaic founder (F0) pigs harboring KISS1-disruptive alleles. Pigs with partial disruption of this gene were used to produce piglets carrying KISS1-edited alleles without mosaicism.
2 Results
2.1 Generation of KISS1-edited pigs via zygote microinjection
To introduce loss-of-function mutations in swine KISS1, a single guide RNA (sgRNA) was designed to target the second exon of this gene, upstream of the sequence encoding the highly conserved KISS1 10-amino acid (aa) core (Figure 1A). The designed sgRNA and Cas9 proteins were co-injected into porcine zygotes with two HDR templates, which separately contained a stop codon (HDRSC) and a silent blocking mutation (HDRBM). These single-stranded donor oligonucleotides (ssODNs) were paired to increase the odds of generating pigs with edits capable of disrupting KISS1 (HDRSC) and/or maintaining KISS1 function (HDRBM; Figures 1B, C). A total of 684 zygotes were injected and transferred into 19 recipient sows. After full-term pregnancies, 69 liveborn (35 males and 34 females) and five stillborn piglets were farrowed in nine litters (Supplementary Table S1).
FIGURE 1. Generation and genotyping of KISS1-edited pigs via zygote injection of CRISPR/Cas9 and single-stranded donor oligonucleotides (ssODNs). (A) Graphics depicting the genomic structure of S. scrofa KISS1, the gene sites bound by the designed sgRNA as well as the homology-directed repair (HDR) templates, and encoding the conserved amino acid core of this protein. At the lower part, it is shown the kisspeptin sequence alignment among different mammalian species in the region flanking the highly conserved core of this protein. (B) Schematic representation of CRISPR/Cas9 RNPs targeting KISS1 besides both ssODNs, and sequences of these HDR templates used to introduce a stop codon and the HindIII restriction site (HDRSC) as well as synonymous SNPs, to change the PAM sequence, to reduce recutting and generate the AcuI restriction site (HDRBM). (C) Schematic of the intracytoplasmic microinjection of CRISPR components and ssODNs into porcine zygotes, transfer of injected embryos, obtention of pigs, and their genotyping for the targeted locus. (D) Bar graph showing gene-edited pigs grouped according to their KISS1-disruptive editing percentages and depicting the allele names and percentages in each individual. In the X-axis are the pig IDs, in which F stands for female, M for male, Y for Yorkshire, D for Duroc, and SB for stillborn.
Analyses of DNA sequences from tail tissues revealed that 53 (72%) piglets were edited, and 18 of these (24%) bore HDR-mediated edits (Table 1). At least 22 pigs bore multiple (>2) alleles. In pigs with two different alleles, 24 pigs had one wild type (WT) allele while four animals had two different editing-derived alleles. AmpSeq indicated that most of the pigs bearing two alleles were mosaic instead of heterozygous (monoallelic edit) or compound heterozygous (different biallelic edits). Three pigs (4MD, 71MY, and 73SB) presented identical biallelic edits, of which there was one pig for each HDR-intended mutation (Figure 1D).
TABLE 1. Genotype summary of pigs produced via microinjection into zygotes of CRISPR/Cas9 RNPs and ssODNs intending to disrupt KISS1.
Detected indels ranged in size from 1 to 68 bp for insertions and 1 to 211 bp for deletions (Supplementary Figure S1), whereas larger indel alleles would not be detected by our AmpSeq assay. Predicted translations indicate that five and 35 alleles contained non-sense and frameshift mutations, respectively. These alleles were predicted to disrupt the KISS1 aa core (Supplementary Figure S2) and were carried by 43 pigs at different frequencies. The percentage of reads attributed to gene-editing-derived alleles in a pig’s genotype was defined as editing percent. In all, 27 animals were generated with KISS1-disruptive alleles between 30% and 90%, with seven pigs having a mosaicism level greater than 90% (Figure 1D).
2.2 Phenotypic analysis
Four pigs (two boars, two gilts) displayed phenotypic characteristics consistent with hypogonadism. Testicular weight and size were dramatically reduced in two boars (51MY and 54MY) at 8 ± 0.6 months of age compared to control animals (Figures 2A–C; Supplementary Table S2). KISS1 WT alleles were absent in both of these two boars. In all other boars, typical variation in testicular volume throughout development was observed, and the variation was not related to the extent to which pig’s KISS1 was disrupted (Figure 2A). The KISS1 KO gilts had KISS1-disruptive allele frequencies of 96% and 100% (49FY and 32FY, respectively) with considerably smaller ovaries and visual absence of surface follicles. The uteri of these KISS1 KO gilts were underdeveloped when compared to age-matched controls (Figures 2D–F; Supplementary Table S3). These findings alongside behavioral observations and hormonal data (Supplementary Figure S3) suggested that the KISS1 KO pigs failed to become pubertal. The other KISS1-edited and WT pigs were considered sexually mature as they displayed genital sizes and external characteristics typical of boars and gilts that advanced through puberty, which is consistent with the hormonal profiles of these animals. Additionally, mosaic F0 pigs that were mated exhibited sexual behaviors, gamete production, and fertility.
FIGURE 2. Reproductive organ phenotyping and growth-related traits of KISS1-edited pigs. (A) Representation of individual testicular volumes collected across different time points, where each panel corresponds to a KISS1-disruption group. “Cry” stands for cryptorchid. (B) In vivo and ex vivo comparison of representative WT and KISS1 KO testicles. (C) Comparison of the testicular weight of WT and KISS1-edited boars. (D,E) Comparison of representative WT and KISS1 KO gilt reproductive tracts (D) and ovaries (E). (F) Follicle count and measurements of ovaries. In (C) and (F), values are individually represented for pigs with KISS1-disruptive editing percentages above 90% and as means for the other groups. (G,H) Body weights (LSmeans ± SE; age, p < .001) (G) and average daily gain of weights (LS means ± SE; birth weight, p < .001) (H) across the evaluated time points (±4 days) and periods, respectively. (I) Scatter plot of total average daily gain of weight (i.e., from birth to 160 days of age) with lines indicating the KISS1 KO pigs. In (G–I), results are grouped by the KISS1-disruptive editing percent.
Piglet birth weight was unaffected by sex, KISS1-disruption group, and their interaction (p ≥ .55). The main effect of age was the only fixed effect to influence body weight (p < .001), but average daily gain was influenced by age and piglet birth weight (p < .001). There was no significant interactive effect of age × KISS1-disruption group for body weight or average daily gain (p ≥ .33; Figures 2G, H). Body weight was lower (p = .04) in WT than 5%–90% KISS1-disruptive pigs at 160 days, but no KISS1-disruption groups differed in body weight or average daily gain at any age category evaluated (p ≥ .07). In this study, the pigs considered as WT were those that carried no KISS1-disruptive alleles, see descriptive statistics in Supplementary Table S4. Boar 15MD was classified in the group >90% KISS1-disruption based on initial sequencing, but his phenotype was remarkably similar to an unedited or partially mosaic boar (<90%). Resequencing revealed him to have 88% of KISS1-disruptive allele frequency instead of 95% as originally determined. Reclassification of this boar did not alter statistical significances or interpretation of the results. It cannot be ruled out that similar cases may have occurred. Nevertheless, no further inconsistencies were observed between genotype and expected phenotype aside from pig 67FY. This gilt had a KISS1-disruptive allele frequency of 97% but had no phenotypic abnormalities.
2.3 Testicular immunohistology
KISS1 KO and WT boars presented an architecture of seminiferous tubules that looked similar. The main difference between these two groups of 8-month-old pigs was related to the size of the observed structures, which were smaller in kisspeptin deficient boars. In the KISS1 KO group, there were no visible lumina and Sertoli cells were arrayed along the parietal edge of these small seminiferous tubules. Germ cells were located both in the adluminal area and toward the edge of the tubules while no spermatids were observed. In the WT boar, the germ cells had all homed to the basement membrane of the germinal epithelium, with none in the seminiferous tubule lumen; germ cells and Sertoli cells were mixed lining the tubule. DAPI staining revealed the presence of cells across the germinal epithelium and adluminal area only in the WT boar, and some sperm tails and elongated spermatids were observed. This indicates that, in this animal, germ cells were undergoing differentiation to eventually become spermatozoa (Figure 3).
FIGURE 3. Immunohistochemical analysis of germ and Sertoli cell markers in seminiferous tubules from WT and KISS1 KO pigs. Boars were 8 months old. DAPI was applied as nuclear staining and is represented in blue fluorescence. Combined images of structure, targeted proteins, and nuclei are shown in the last column. Scale bar = 100 µm.
2.4 Hormonal profiling
Hormone measurements were conducted every 4 weeks beginning when pigs were 40 days old and halved to a fortnight frequency from 130 to 280 days of age. Only follicle-stimulating hormone (FSH) of gilts showed an interaction of age with KISS1-disruptive editing percent (p = .05), but least-squares means of this interactive effect are displayed for all hormones in Figure 4 to describe patterns over time. Gilt FSH was lower for the >90% KISS1-disruption group than all others at 40 and 70 days of age (p ≤ .05; Figure 4B), but no differences between groups were observed thereafter (p ≥ .20). Serum concentrations of boar FSH were not affected by KISS1-disruptive editing percent (p = .39; Figure 4A). All hormones, except for boar FSH (p = .07), were influenced by the main effect of age (p < .001). Overall, boars with KISS1-disruptive allele frequency >90% had less (p < .001) serum luteinizing hormone (LH) and testosterone during all ages of development compared with other groups, which did not differ from one another (p ≥ .91; Figures 4C, E). Concentrations of LH in gilts with KISS1-disruption >90% were numerically less than in other groups, but this did not reach significance (p = .07; Figure 4D). Individual concentrations of serum hormones for some pigs belonging to the group with KISS1-disruptive allele frequency >90% are shown in Supplementary Figure S3.
FIGURE 4. Hormone profiles of KISS1-edited pigs. (A–D) Serum concentration (LSmeans ± SE) in boars and gilts of FSH [(A) and (B), respectively] and LH [(C) and (D), respectively] during development by KISS1-disruption group. (E) Developmental change in serum concentrations (LS means ± SE) of testosterone in boars according to the KISS1-disruptive editing percent.
2.5 F0-to-F1 transmission of KISS1-edited alleles
Intergenerational transmission of the KISS1-mediated castration free trait is required to provide a more effective alternative to end surgical castration. Semen from five F0 boars (13MY, 36MY, 61MY, 53MY, and 15MD) with KISS1-disruptive allele frequencies between 35% and 88%, was collected for analysis (Supplementary Table S5). KISS1 gene edits in sperm DNA were evaluated by NGS and compared to results obtained from the tail samples. Each mutant allele detected in the tail samples could be found in the sperm DNA, confirming germline transmission. However, apart from 13MY and 61MY, there was discordance ranging from 16% to 43% between the allele frequencies observed in the tail versus those in sperm DNA (Figure 5).
FIGURE 5. Confirmation of KISS1-edited allele transmission. (A) KISS1 allele frequencies of mosaic F0 boar candidates for breeding, calculated using DNA from tail tissue and spermatozoa. (B) Inherited KISS1 alleles by F1 piglets, and the number of genotyped pigs in each litter according to their sex-breed (first bar) or genotype category (second bar). “DEP” stands for disruptive editing percent and “Δ” means deletion.
Three F0 boars and 10 F0 gilts were selected for mating based highest levels of KISS1 disruption under 90 percent and having the lowest inbreeding coefficient. From the first matings, seven litters were farrowed (Figure 5) with an average of 12.9 ± 1.6 liveborn piglets per litter, for a total of 90 liveborn F1 piglets. The genotype distribution was, 13 KISS1−/−, 29 KISS1+/− and 40 KISS1+/+; 44 males and 38 females. For second matings, five and one of the previously chosen F0 gilts and boars, respectively, were bred again. Eighty liveborn F1 piglets were farrowed through five litters (35 males and 30 females; average of 16 ± 1.7 piglets per litter), from which 21 piglets were KISS1−/−, 31 KISS1+/− and 13 KISS1+/+ (Table 2). In nine piglets farrowed in four litters (from 47FY and 65FY), only the paternal alleles were detected (Figures 1D, 5), suggesting that these two mothers harbored alleles with large deletions that could not be detected with the implemented PCR test. Inherited alleles of genotyped F1 piglets are shown in Figure 5 according to the litter.
3 Discussion
In laboratory and livestock species, function of the kisspeptin system is needed for the proper regulation of gonadotropin secretion and reproduction (Seminara et al., 2003; d’Anglemont et al., 2007; Sonstegard et al., 2017). However, empirical evidence about the role of KISS1 in pig reproduction is required to postulate that disruption of the gene encoding the kisspeptin ligand can be useful as an alternative to surgical castration. Here, we have generated four pigs harboring KISS1 KO alleles in high frequencies that were phenocopies of humans and animals with HH. They were produced alongside KISS1-edited mosaic pigs that upon mating gave origin to 34 KISS1−/− F1 piglets, which confirmed the germline transmission of edited alleles. The generation of the KISS1 KO F0 individuals shows that the kisspeptin ligand plays a determinant role in initiating the puberty of pigs. We observed that some animals with >80% KISS1-disruption still achieved sexual maturity, suggesting that even low concentrations of this peptide can trigger puberty.
Mosaicism in gene-edited zygotes is affected by variation in concentrations of gene editing reagents (Zhou et al., 2015; Tanihara et al., 2019; Menchaca et al., 2020). In this study, we sought to produce heterozygous-like KO founders to enable their breeding and generation of mosaic-free offspring for eventual deeper phenotyping. To accomplish this, we optimized sgRNA/Cas9 and ssODN levels in vitro (data not shown) and also co-injected two ssODNs, one designed to knock out KISS1, and the second designed to repair KISS1 with silent alterations in the sgRNA biding site to prevent further cleavage. This combined approach resulted in intermediate editing with variable mosaicism in the F0 population, allowing puberty and breeding of subsequent generations. While mosaicism in this case was the intent, it is a common outcome in embryo editing and occurs when the genome editor persists past the one-cell stage. This results in unequal allele frequencies and/or more than two different alleles within the same individual (Mehravar et al., 2019; Hennig et al., 2020).
Although mosaicism can confound phenotyping in founder animals, when the latter is analyzed together with the genotypes, some clues about the function of the gene may be obtained. Four mosaic pigs with elevated KISS1-disruptive allele frequencies (>90%), displayed typical phenotypes of individuals with non-functional KISS1; hypogonadotropism, hypogonadism, failure to achieve puberty, and lack of gamete production (d’Anglemont et al., 2007; Lapatto et al., 2007; Uenoyama et al., 2015; Ikegami et al., 2020). In agreement with the heterozygosity of KISS1 mutations (Lapatto et al., 2007), pigs of intermediate KISS1-disruptive editing (30%–70%) appeared to be reproductively normal. The absence of HH was also observed in pigs with above-intermediate frequency (>70%–90%) and even in a gilt (67FY) with considerably high frequency (97%) of KISS1-disruptive alleles. Further, ovulation was confirmed by postmortem evaluation of ovaries and serum concentrations of progesterone (data not shown) in some gilts of >80% KISS1-disruption. This suggests that a low amount of unaltered KISS1 is sufficient to stimulate the hypothalamic release of GnRH to support gonadotropin secretion for gonadal maturation and puberty, as evidenced by the pattern of LH secretion between highly edited individuals. Studies with rodents suggest that activation of only 10%–20% of GnRH neurons is necessary for an ovulatory surge of LH to be generated (Greig and Weisz, 1973; Gosden and Everett, 1976; Herbison et al., 2008). It is unknown how many kisspeptin neurons are required to activate this proportion of GnRH neurons, but with the data herein it can be hypothesized it is few.
In this regard, it is important to consider that KISS1 allele frequencies clearly differed between tail tissue and sperm samples in three out of five mosaic F0 boars tested. If the same happened between hypothalamic and tail tissues, disparity could contribute to the absence of HH in some pigs classified as having high KISS1-disruptive mosaicism (>85%; e.g., 15MD and 67FY). This can occur when the tail tissue used for genotyping provides an allele frequency similar to the one observed in KO animals, but the hypothalamic region where KISS1 is produced has a lower KISS1-disruptive allele frequency to the point there are functional neurons to produce enough kisspeptin. Differences in allele frequencies among organs may explain why gilt 67FY was not KO with a KISS1-disruptive editing percent of 97% but that 49FY was KO with a KISS1-disruptive allele frequency of 96%. The discrepancy of KISS1 allele frequencies between tail tissue and semen samples ranged between 2% and 43%. Similarly, Tanihara et al. (2020) reported that their greatest difference in allele frequency was 38% when comparing different organs from a single mosaic pig derived from editing embryos with CRISPR/Cas9. Another consideration is that WT cells in mosaic embryos could have preferentially developed into neurons by the process of blastocyst complementation (Kobayashi et al., 2010). However, that outcome is unlikely since KISS1 KO does not prevent the formation of neurons, eliminating the selective pressure for WT neurons to out-compete KISS1 KO neurons.
Overall, patterns of serum hormones during development for both gilts and boars were consistent with expected responses for pigs (Allrich et al., 1982; Berardinelli et al., 1984; Lutz et al., 1984; Ford et al., 2001; Barb et al., 2010) with the exception of FSH in boars (Schinckel et al., 1984), but this latter observation was not due to disruption of KISS1 alleles. Taken together, there was no effect of gene editing on LH in gilts with >90% of KISS1-disruption, but FSH secretion early in development was suppressed in these animals. Serum concentrations of LH in boars with a frequency of KISS1-disruptive alleles >90% was severely suppressed throughout development. This in turn resulted in dramatic reduction in testosterone production and lack of age-related increases in circulating concentrations of testosterone expected for boars over a similar age range (Ford et al., 2001). Lower serum concentrations of testosterone in boars with high levels of KISS1 disruption portends that they would likely have a reduction in other androgen compounds (androstenone and skatole) that would cause boar taint (Zamaratskaia and Squires, 2009).
Phenotypical variability for hypogonadism, uterine weight, and other similar reproductive-related observations have been described in Kiss1 KO and Kiss1r KO rodents (d’Anglemont et al., 2007; Lapatto et al., 2007; Uenoyama et al., 2015). We are currently assessing these phenotypes in KISS1-disrupted F1 pigs lacking mosaicism. Although bilateral cryptorchidism was reported in a human and rats with impaired kisspeptin systems (Semple et al., 2005; Uenoyama et al., 2015), there is no clear link between these conditions. In this study, one pig presented unilateral cryptorchidism, but its frequency of KISS1-disruptive alleles was considerably lower (60%) than that of KISS1 KO boars (100%). Cryptorchidism is one of the most common congenital defects in pigs (Mattsson, 2011) and it that impaired testicular descent of boar 53MY is likely unrelated to KISS1 disruption.
Spermatozoa can typically be found in histological sections of pig testes by 120 days (Malmgren et al., 1996; Lee et al., 2014). Testicular histology of KISS1 KO pigs, at 255–262 days of age, showed a lack of spermatozoa, with atrophied seminiferous tubules, but otherwise normal in structure. Similar findings in KISS1R KO pigs were reported by Sonstegard et al. (2017), who observed hypogonadism and lack of gamete production in KO boars. At 8 months of age, spermatogenesis was observed in the seminiferous tubules of WT boars, while in KISS1 KO boars the germ cells were still randomly arranged within the tubules and failed to undergo spermatogenesis. This distribution pattern of germ cells is typical of pigs whose ages range between five and 60 days. Boar germ cells have been reported to complete migration to the basement membrane of the seminiferous tubule by 90 days of age (Lee et al., 2014). It is believed that the placement of germ cells at the basement membrane is an important environmental stimulus for further maturation (Nagano et al., 2000). Impaired germ cell relocation may be one of the events preventing spermatogenesis in KISS1 KO boars.
Breeding of mosaic F0 boars and gilts resulted in non-mosaic F1 piglets, some of them harboring KISS1-disruptive alleles in either the heterozygous or homozygous state. Two F0 mothers transmitted alleles to some of their piglets that were not amplified due to the deletion of at least one of the primer-annealing sites, which are required for genotyping of the targeted locus. This indicates the presence of on-target deletions larger than detected by our PCR assay. The occurrence of these large deletion events at a target site has been documented using CRISPR/Cas9-mediated editing of embryos in different species (Owens et al., 2019; Koppes et al., 2020; Korablev et al., 2020; Alanis-Lobato et al., 2021; Höijer et al., 2022). A number of aspects could affect the presence and abundance of larger indels at on-target sites such as the delivery method of the CRISPR/Cas9 system, its concentration, and the use of HDR donor templates (Wen et al., 2021; Höijer et al., 2022). It is thought that the differential kinetics of non-homologous end joining, HDR, or microhomology-mediated end joining repair machinery influence the frequency of large deletions (Wen and Zhang, 2022).
Generation of offspring from F0 pigs confirmed allele transmission. Mosaic founders with KISS1 disruption between 25% and 70% had normal pregnancy rates and litter sizes in the first and second parity. KISS1 is attributed a role in pregnancy as this peptide is believed to inhibit trophoblast migration and invasion in the placenta (Cao et al., 2019). In the case of boars, Zou et al. (2019) reported a positive relationship between KISS1 concentration in seminal plasma and total sperm count as well as total motile sperm count. Although minimal numbers of ejaculates were assessed in the current study, they were similar to what one would expect for commercial boars. Volume and quality of ejaculate increased as boars became fully trained to semen collection.
Variation in growth-related traits was not affected by the percent of KISS1 disruption (Figures 2G–I), and F0 KISS1 KO pigs were within the average for the swine industry. Similar results were found by Lapatto et al. (2007) but not by d'Anglemont et al. (2007) for the body weight of Kiss1 KO mice compared to WT mice. Likewise, Kiss1 KO rats were reported to have normal growth (Uenoyama et al., 2015). No growth impairment due to KISS1 disruption in swine was expected because, unlike in sheep and cattle (Kadokawa et al., 2008; Whitlock et al., 2008; Foradori et al., 2017), there is no link between kisspeptin and the somatotropic axis in pigs (Lents et al., 2008). This is relevant because adequate growth performance of KISS1-deficient pigs is necessary for wide acceptance in the swine industry.
Demonstrating that KISS1 KO in pigs prevents puberty supports conducting additional performance studies at the appropriate statistical power. This will enable comparison not only of growth, but also of other phenotypes of interest like feed efficiency and boar taint compound levels from KO and WT animals. While these studies are essential for producer and consumer acceptance, they are not warranted until a breeding strategy can be established. The results in this manuscript provide valuable histological and hormonal results for use as a baseline in follow-on fertility restoration studies.
Producers are not likely to adopt approaches in which perpetuating the trait involves gene editing each litter or insertion of foreign DNA sequences into the genome of the breeding stock. For example, knocking out the SRY gene to induce male-to-female sex reversal would require direct editing of embryos for every litter, and integration into the parental genomes of a CRISPR/Cas vector would be needed for the production of only-female litters (Kurtz et al., 2021). Likewise, editing genes to reduce the biosynthesis of boar taint precursors would be difficult to implement due to the large number of enzymes involved, and steroid biosynthetic precursors can flow to synthetic endpoints through multiple pathways (Rydhmer et al., 2006; Squires et al., 2019). Alternatives in which both boar taint and aggressiveness can be prevented and spread through selective breeding are more likely to be embraced by producers.
The current results demonstrate that kisspeptin ligand activation of the hypothalamic-pituitary-gonadal axis is required to induce puberty in pigs. Other mechanisms in this species cannot autonomously compensate for the loss of KISS1 function. These findings also confirm the inheritance of KISS1-disruptive alleles and suggest that a low amount of functional KISS1 may be required to trigger the onset of puberty in pigs. The goal of this and ongoing research with the F1 generation is to improve the welfare of millions of pigs globally, and the safety of the personnel who handle them, by preventing surgical castration of male piglets, while ensuring good pork quality. Achieving this will depend on further investigation to provide critical information about how to manage the reproductive phenotypes of these KISS1-edited animals.
4 Materials and methods
4.1 Ethics approval
All animal procedures followed institutional guidelines for the care and use of animals and were approved by the Recombinetics Institutional Animal Care and Use Committee (IACUC; approval number: RCI-1903-03A-A3).
4.2 Animals
The F0 pigs were of Duroc or American Yorkshire breeding, derived from commercial herds, and farrowed from first and second parity surrogate sows. The F0 boars and gilts chosen for breeding resulted in both purebred and crossbred animals. An industry standard fortified corn-soybean meal diet formulated to meet and exceed nutrient requirements was provided with ad libitum access to water.
4.3 Generation of KISS1-edited pigs using sgRNA/Cas9 RNPs and ssODNs
Sus scrofa KISS1 possesses a genomic structure of two exons that encode 138 aa. Amino acids 111–120 form a 10 aa core conserved across vertebrates that is essential for binding and activation of KISS1R. A sgRNA was designed to target a KISS1 exonic region preceding the sequence encoding this completely conserved aa core motif. The online tools Gene Sculpt Suite MENTHU (Ata et al., 2018) and microhomology predictor (Bae et al., 2014) were used to design the sgRNA with the sequence of 5′-GGTCCCCCGAGGGTTCGCCTCGG-3′ (PAM is underlined). crRNA and tracrRNA (IDT, IA, United States) were resuspended in injection buffer (5 mM Tris-HCl, 0.1 mM EDTA, pH 7.5) and incubated with HiFi Cas9 protein (IDT, IA, United States) to form RNP complexes. The assembled RNPs were co-injected in fertilized porcine zygotes with two ssODNs of 90 bp (IDT, IA, United States) that served as templates for HDR-mediated editing (Figure 1B). The first ssODN, HDRSC, introduces an 8-bp insertion to generate a premature termination codon followed by the HindIII restriction site. The second ssODN, HDRBM, produces two silent mutations at and close to the PAM site that reduce Cas9 recutting after HDRBM editing while also introducing the restriction site AcuI.
Single-cell embryos for microinjection were produced in vivo by synchronizing the ovulation of 71 sows with 200 μg of the GnRH agonist triptorelin acetate (OvuGel; JBS United Animal Health, IN, United States) 96 h (h) post-weaning. At 22 ± 2 h after triptorelin acetate treatment, sows were artificially inseminated with extended boar semen. The following day, sows were transferred to the slaughterhouse at the Meat Science Laboratory at the University of Minnesota, where the reproductive tracts were collected. The zygotes were recovered by flushing the oviduct with phosphate-buffered saline (1X PBS; Corning, NY, United States). A total of 836 in vivo fertilized zygotes were obtained. Single-cell embryos were immediately transferred into TCM-199 medium (Gibco, MD, United States) and 3–5 h after zygote recovery, 684 one-cell stage embryos were injected. Zygote microinjection was done into the cytoplasm with sgRNA (25 ng/μl)/Cas9 (50 ng/μl) RNPs and ssODNs (HDRSC: 33.3 ng/μl; HDRBM: 66.7 ng/μl) using an inverted microscope (Nikon, Tokyo, Japan), which was equipped with micromanipulators (Narishige, Tokyo, Japan) and an electronic microinjection system (Femtojet; Eppendorf, Hamburg, Germany). Injected zygotes were cultured for 12–16 h in PZM-3 medium (Yoshioka et al., 2002) covered with mineral oil, at 38°C in a humidity-controlled atmosphere containing 5% CO2 and 5% O2.
Within approximately 20–24 h from zygote collection, an average of 36 microinjected zygotes were surgically transferred into oviducts of 19 hormonally synchronized surrogate sows. To synchronize the estrous of recipient sows, they were treated with 18 mg/day of altrenogest (MATRIX; Merck & Co., NJ, United States) for 18 days, plus 1,000 IU of human chorionic gonadotropin (hCG; Chorulon, Merck & Co., NJ, United States) on day 19, and embryo transfers were carried out 48 h later. Pregnancies were confirmed and monitored on days 30 and 60 after the embryo transfers using an Aloka 500 Ultrasound Scanner (Aloka Co. Ltd., CT, United States). After 115–118 days of gestation, 10 recipient females farrowed naturally, and piglets were genotyped.
4.4 Genotyping
4.4.1 F0 and F1 pigs
Genomic DNA was extracted from tail docking biopsies of all (alive or dead) F0 and alive F1 piglets born using the DNeasy Blood and Tissue Kit (QIAGEN, MD, United States). DNA was subjected to PCR amplification of the CRISPR/Cas9 targeted locus with specific primers (F: 5′-GGATGAGCAAACGGTCCAGA-3′and R: 5′-CTCCCGGGTTTGAAGGTCTC-3′; WT amplicon size: 409 bp) using AccuStart II GelTrack PCR SuperMix (QuantaBio, MA, United States). Amplicons were purified with the QIAquick PCR Purification Kit (QIAGEN, MD, United States) and sent to Sanger sequencing (ACGT Inc., IL, United States) for both F0 and F1 pigs. The resulting trace files were analyzed with the bioinformatics package ICE (Conant et al., 2022). For a better quantification of the mutations and more accurate deciphering of genotype in pigs generated by embryo editing, amplicons of 45 F0 pigs were sequenced using Amplicon-EZ NGS (GENEWIZ, NJ, United States). These 45 piglets were chosen because the Sanger sequencing results indicated that they could have the potential to be used as parents to generate the KISS1-edited F1 generation, or that they could be KISS1 KO animals.
NGS data was demultiplexed with Geneious 7.1.91. Reads were analyzed using the CRISPR RGEN Tools Cas-Analyzer software (Park et al., 2017), which provided a detailed breakdown of the allele frequency for the targeted KISS1 locus. Predicted translations of the detected alleles, as well as DNA and protein sequence alignments, were generated utilizing Geneious 7.1.9; this information was used to estimate the potential of the generated alleles to disrupt KISS1. Two alleles were classified as of unknown translation since their corresponding Sanger trace files did not allow proper allele deciphering (Figure 1D).
4.4.2 Semen of KISS1-edited boars
Ejaculated semen samples from F0 boars were subjected to genomic DNA isolation using the QIAamp DNA Mini Kit and an adapted QIAamp Tissue Protocol (QIAGEN, MD, United States). DNA obtained from semen was PCR-amplified, amplicons were purified as described above, and PCR products were NGS-sequenced. Amplicons were submitted for NGS to the US Meat Animal Research Center (NE, United States) and downstream analyzed as described in the previous paragraph.
4.5 Phenotypic measures
Body weights of all F0 pigs were individually recorded within 48 h after birth (birth weight) and at 40, 70, 100, 130, and 160 days of age ±4 days. The growth rate was calculated as the average daily gain of weight for the time periods birth-40, 40–70, 70–100, 100–130, and 130–160 days of age. The total average daily gain of weight was estimated from birth to 160 days of age. Beginning at 40 days of age and continuing, at each time body weight was collected, the length and width of both testicles were individually measured with a vernier caliper (35-OD8; iGaging, CA, United States).
Testicular volume was estimated using the equation for a prolate spheroidal shape (Young et al., 1986):
Blood samples were collected from F0 pigs monthly starting at 40 days of age and biweekly from 130 to 280 days of age by jugular venipuncture into serum-separating tubes. Serum was obtained by centrifugation at 1,500 × g for 15 min (min) at 4°C, aliquoted, and stored at −80°C until analysis. Serum concentrations of LH (Kesner et al., 1987), FSH (Trout et al., 1992), testosterone (Ford et al., 2001), and progesterone (Calderón et al., 2017) were quantified with validated radioimmunoassay. Testosterone and progesterone assays were obtained from MP Biomedicals (Irvine, CA, United States). The reference standards for LH (AFP-10506A) and FSH (AFP-10640B) were provided by Dr. A. F. Parlow (Scientific Director of the NIH, NIDDK, National Hormone and Peptide Program, Torrance, CA, United States). Pools of porcine serum were included in each assay. For LH, concentrations of pools measured .81, 1.54, and 8.41 ng/ml with an average intra- and inter-assay coefficient of variation (CV) of 6.7% and 10.8%, respectively (n = 8 assays). For FSH, serum concentrations of pools ranged from 2.0 to 6.8 ng/ml with average intra- and inter-assay CV of 6.3% and 7.3%, respectively (n = 6 assays). For progesterone, a serum pool (midluteal phase) measuring 24.2 ng/ml had an intra- and inter-assay CV of 5.7% and 8.6%, respectively (n = 4 assays). Intra- and inter-assay CV were 10.7% and 13.8% (n = 6 assays), respectively, for a pool of boar serum with testosterone concentration of 5.3 pg/ml.
4.6 Gross morphology
Pigs were humanely euthanized according to established guidelines (American Veterinary Medical Association, 2020) for tissue collection. The length and width of ovaries and testes were individually measured with a vernier caliper, and the number of ovarian follicles was counted. Ovaries and testes were individually weighed. External and internal genitalia were observed for visible abnormalities. The number of pigs that were evaluated postmortem as well as their minimum and maximum ages at this point, for each KISS1-disruption group, are presented in Supplementary Tables S2, S3.
4.7 Histology and immunohistochemistry
Testes from F0 boars at 255–262 days of age were collected, bivalved, cut through the longest dimension, and 0.5 mm × 0.5 mm pieces containing sections of mediastinum and periphery were collected. These tissues were placed into 15 ml conical tubes with 10% Neutral Buffered Formalin (89370-94; VWR, PA, United States), and maintained at room temperature for 18 h. Tissues were then washed three times in PBS with rocking for 10 min. Samples were placed in 70% ethanol (V1001TP; Decon Laboratories Inc., PA, United States) and stored at room temperature until submission to Scientific Solutions (MN, United States) for paraffin embedding. Tissue sections of 20 microns were mounted on slides (Globe Scientific Inc., NJ, United States) and deparaffinized with Xylene (89370-008; VWR, PA, United States) followed by rehydration with increasing concentrations of 70%–100% ethanol. Antigens were retrieved by boiling the tissue in a microwave for 15 min with a citrate-based unmasking solution (1:100; H-3300; Vector Laboratories Inc., CA, United States). Membrane permeabilization was accomplished by incubating tissue sections in PBS containing 0.1% Triton X-100 (Sigma-Aldrich, MO, United States) for 10 min at room temperature. To block the binding of non-specific proteins, tissue sections were incubated with 10% normal goat serum (Sigma-Aldrich, MO, United States) in PBS containing 0.025% tween-20 (PBST; Bio-Rad, CA, United States) overnight at 4°C. Washes with PBS were done twice for 10 min after blocking, primary antibody stain, secondary stain, and DAPI staining.
Tissue sections were incubated overnight at 4°C with primary antibodies diluted to manufacturer recommendations with 5% normal goat serum-PBST. Vimentin (1:500 dilution; Abcam, Cambridge, UK) and GATA binding protein 4 (GATA4; 1:500 dilution; Cell Signaling Technology, MA, United States) were used as markers for Sertoli cells. Deleted in azoospermia-like protein (DAZL; 1:500 dilution; Abcam, Cambridge, UK) and protein gene product 9.5 (PGP9.5; 1:1,000 dilution; Abcam, Cambridge, UK) were used as germ cells markers. Tissue sections were incubated for 1 h at room temperature in secondary antibodies diluted 1:400 in 5% normal goat serum-PBST. These antibodies were Alexa Fluor 488 (anti-rabbit) or Alexa Fluor 594 (anti-mouse; Invitrogen, MA, United States). Details of antibodies used are in Supplementary Table S6. Nuclei were stained with Hoechst solution (Thermo Fisher Scientific Inc., MA, United States) diluted 1:10,000 in PBS for 15 min. Tissue sections were washed with PBS (5 min, room temperature) and coverslipped using ImmuMount (Fisher Scientific, MA, United States), followed by incubation at room temperature for 30–120 min. A Leica DM6000B microscope equipped with a Leica DFC7000 T camera (Leica Microsystems, Wetzlar, Germany) was used to view and image slides.
4.8 Statistical analysis
Descriptive statistics of all anatomical-related phenotypes and statistical analyses of body weight, average daily gain, and hormone concentration were estimated using SAS 9.4 (SAS Institute Inc., NC, United States). The WT group always consisted of all pigs that harbored 0% of KISS1-disruptive alleles, including individuals carrying KISS1-non-disruptive alleles at different frequencies as these did not differ phenotypically from unedited pigs.
4.8.1 Body weight and growth rate
The statistical model for birth weight (n = 58) included fixed effects of sex, KISS1-disruption group (WT, 5%–90%, or >90%), and their interaction along with a random litter effect (nine levels). Body weight and average daily gain (n = 287 each) were analyzed as repeated measures with fixed effects of sex, KISS1-disruption group, age (five levels), and their two- and three-way interactions. Birth weight was included as a linear covariate. A type I autoregressive (co)variance structure with heterogenous variance across age was chosen to model the random pig effect.
4.8.2 Hormone analysis
Box plots were initially constructed using all available observations of boar testosterone (n = 234), boar LH (n = 233), boar FSH (n = 235), gilt LH (n = 246), and gilt FSH (n = 247) and outliers were identified as those having values greater or less than three times the interquartile range within each sampling age. After removing outliers, 226, 231, 240, and 244 observations remained for analysis of boar testosterone, boar LH, gilt LH, and gilt FSH, respectively. Little variability was observed in boar FSH after 160 days of age and these data were removed, leaving 169 observations for analysis. A natural logarithm transformation was applied to all hormones, which were analyzed as repeated measures within sex with fixed effects of KISS1-disruption group (WT, 5%–90%, or >90%), age (up to nine levels), and their interaction. A type I autoregressive (co)variance structure with heterogenous variance across age was chosen to model the random pig effect. Least-squares means and their standard errors were back-transformed to the scale of measurement using the delta method.
4.9 Semen collection of KISS1-edited boars
F0 boars that were candidates to be bred were trained on a semen-collection dummy. When mosaic boars were between 272 and 329 days of age, semen was collected via the hand-gloved method. Ejaculates were quantified, samples were sent for NGS-mediated genotyping as outlined above, and replicate semen was extended for evaluation of sperm motility, morphology, and concentration with a light microscope (CxL; Labomed, CA, United States).
4.10 Artificial insemination using KISS1-edited boars and gilts
In mosaic gilts that were candidates to be mated due to their genotypes, estrus was synchronized using altrenogest for at least 14 days. Gilts that displayed classical signs of estrus were artificially inseminated with semen of F0 boars bearing KISS1-disruptive alleles. First matings by artificial insemination occurred when gilts were between 275 and 355 days of age and boars were 317 and 348 days old. A second mating was performed 9 months later. Semen collections and pregnancy diagnoses were done as previously mentioned.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was reviewed and approved by the Recombinetics Institutional Animal Care and Use Committee.
Author contributions
TS conceived the study. TS, CL, AH, and JF designed the study. SS, UG, JB, and AC did the gene editing experiments under the supervision of DC. PR-V, FO, SS, and UG coordinated and did the reproductive biotechnology techniques using microinjection. JF, KM, and JB genotyped the pigs. KM led and together with JF and AC participated in the collection of phenotypes. KM did the histology and immunohistochemistry. CL ran the hormonal assays and analyzed them. TM did the statistics for the phenotypic data. E-SK conducted genomics analyses. SL, TS, AH, KM, CL, and DK supervised, coordinated, and/or contributed to the logistics of the project. JF prepared Figures, Tables, and wrote the manuscript with revision and text editing contributions from CL. All authors contributed to the article and approved the final version.
Funding
The funding sources for this research were mainly awarded by The Foundation for Food and Agriculture Research (FFAR; grant number: 552176). Additional supporting funds were provided through a collaborative agreement between Hendrix Genetics and Acceligen Inc.
Acknowledgments
The authors sincerely thank Jun Chen, Chris Fahrenkrug, and Mandy Taisto for their help in the collection of phenotypic measurements and gross morphology data. We would also like to thank Shane Hodgson for the technical support in histology, Tim Smith for sequencing boar ejaculates, and the farm staff for their good animal care. JF was funded by the scholarship program Pasaporte a la Ciencia in the initial stages of this research.
Conflict of interest
The authors JF, KM, JB, PR-V, FO, SL, E-SK, and TS are employed by Acceligen Inc., a wholly owned subsidiary of Recombinetics Inc.; SS, UG, AC, and DC are employed by Recombinetics Inc. AH is an employee of Hendrix Genetics. These companies contributed supplemental funding for this study.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author Disclaimer
Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the USDA. The USDA prohibits discrimination in all its programs and activities on the basis of race, color, national origin, age, disability, and where applicable, sex, marital status, familial status, parental status, religion, sexual orientation, genetic information, political beliefs, reprisal, or because all or part of an individual’s income is derived from any public assistance program (not all prohibited bases apply to all programs). Persons with disabilities who require alternative means for communication of program information (Braille, large print, audiotape, etc.,) should contact USDA’s TARGET Center at (202) 720-2600 (voice and TDD). To file a complaint of discrimination, write to USDA, Director, Office of Civil Rights, 1400 Independence Avenue, S.W., Washington, DC 20250-9410, or call (800) 795-3272 (voice) or (202) 720-6382( TDD). USDA is an equal opportunity provider and employer.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1078991/full#supplementary-material
Footnotes
References
Alanis-Lobato, G., Zohren, J., McCarthy, A., Fogarty, N. M. E., Kubikova, N., Hardman, E., et al. (2021). Frequent loss of heterozygosity in CRISPR-Cas9–edited early human embryos. Proc. Natl. Acad. Sci. U. S. A. 118, e2004832117. doi:10.1073/pnas.2004832117
Allrich, R. D., Christenson, R. K., Ford, J. J., and Zimmerman, D. R. (1982). Pubertal development of the boar: testosterone, estradiol-17β, cortisol and LH concentrations before and after castration at various ages. J. Anim. Sci. 55, 1139–1146. doi:10.2527/jas1982.5551139x
American Veterinary Medical Association (2020). AVMA guidelines for the euthanasia of animals: 2020 edition. Available at: https://www.avma.org/sites/default/files/2020-02/Guidelines-on-Euthanasia-2020.pdf (Accessed February 24, 2022).
Ata, H., Ekstrom, T. L., Martínez-Gálvez, G., Mann, C. M., Dvornikov, A. V., Schaefbauer, K. J., et al. (2018). Robust activation of microhomology-mediated end joining for precision gene editing applications. PLoS Genet. 14, e1007652. doi:10.1371/journal.pgen.1007652
Backus, G., Higuera, M., Juul, N., Nalon, E., and de Briyne, N. (2018). Second progress report 2015 – 2017 on the European declaration on alternatives to surgical castration of pigs. Brussels Available at: https://www.boarsontheway.com/wp-content/uploads/2018/08/Second-progress-report-2015-2017-final-1.pdf (Accessed December 12, 2022).
Bae, S., Kweon, J., Kim, H. S., and Kim, J.-S. (2014). Microhomology-based choice of Cas9 nuclease target sites. Nat. Methods 11, 705–706. doi:10.1038/nmeth.3015
Barb, C. R., Hausman, G. J., and Kraeling, R. R. (2010). Luteinizing hormone secretion as influenced by age and estradiol in the prepubertal gilt. Anim. Reprod. Sci. 122, 324–327. doi:10.1016/j.anireprosci.2010.10.006
Berardinelli, J. G., Ford, J. J., Christenson, R. K., and Anderson, L. L. (1984). Luteinizing hormone secretion in ovariectomized gilts: effects of age, reproductive state and estrogen replacement. J. Anim. Sci. 58, 165–173. doi:10.2527/jas1984.581165x
Bonneau, M., and Weiler, U. (2019). Pros and cons of alternatives to piglet castration: welfare, boar taint, and other meat quality traits. Animals. 9, 884. doi:10.3390/ani9110884
Calderón, J. A., Vallet, J. L., Boyd, R. D., Lents, C. A., Prince, T. J., DeDecker, A. E., et al. (2017). Effect of feeding three lysine to energy diets on growth, body composition and age at puberty in replacement gilts. Anim. Reprod. Sci. 184, 1–10. doi:10.1016/j.anireprosci.2017.06.007
Cao, Y., Li, Z., Jiang, W., Ling, Y., and Kuang, H. (2019). Reproductive functions of Kisspeptin/KISS1R Systems in the periphery. Reprod. Biol. Endocrinol. 17, 65–69. doi:10.1186/s12958-019-0511-x
Conant, D., Hsiau, T., Rossi, N., Oki, J., Maures, T., Waite, K., et al. (2022). Inference of CRISPR edits from sanger trace data. CRISPR J. 5, 123–130. doi:10.1089/crispr.2021.0113
d’Anglemont, X., Fagg, L. A., Dixon, J. P. C., Day, K., Leitch, H. G., Hendrick, A. G., et al. (2007). Hypogonadotropic hypogonadism in mice lacking a functional Kiss1 gene. Proc. Natl. Acad. Sci. U. S. A. 104, 10714–10719. doi:10.1073/pnas.0704114104
de Roux, N., Genin, E., Carel, J.-C., Matsuda, F., Chaussain, J.-L., and Milgrom, E. (2003). Hypogonadotropic hypogonadism due to loss of function of the KiSS1-derived peptide receptor GPR54. Proc. Natl. Acad. Sci. U. S. A. 100, 10972–10976. doi:10.1073/pnas.1834399100
Dunshea, F. R., Colantoni, C., Howard, K., McCauley, I., Jackson, P., Long, K. A., et al. (2001). Vaccination of boars with a GnRH vaccine (Improvac) eliminates boar taint and increases growth performance. J. Anim. Sci. 79, 2524–2535. doi:10.2527/2001.79102524x
Foradori, C. D., Whitlock, B. K., Daniel, J. A., Zimmerman, A. D., Jones, M. A., Read, C. C., et al. (2017). Kisspeptin stimulates growth hormone release by utilizing neuropeptide Y pathways and is dependent on the presence of ghrelin in the ewe. Endocrinology 158, 3526–3539. doi:10.1210/en.2017-00303
Ford, J. J., Wise, T. H., Lunstra, D. D., and Rohrer, G. A. (2001). Interrelationships of Porcine X and Y chromosomes with pituitary gonadotropins and testicular size. Biol. Reprod. 65, 906–912. doi:10.1095/biolreprod65.3.906
Funes, S., Hedrick, J. A., Vassileva, G., Markowitz, L., Abbondanzo, S., Golovko, A., et al. (2003). The KiSS-1 receptor GPR54 is essential for the development of the murine reproductive system. Biochem. Biophys. Res. Commun. 312, 1357–1363. doi:10.1016/j.bbrc.2003.11.066
Gosden, R. G., Everett, J. W., and Tyrey, L. (1976). Luteinizing hormone requirements for ovulation in the pentobarbital-treated proestrous rat. Endocrinology 99, 1046–1053. doi:10.1210/endo-99-4-1046
Greig, F., and Weisz, J. (1973). Preovulatory levels of luteinizing hormone, the critical period and ovulation in rats. J. Endocrinol. 57, 235–245. doi:10.1677/joe.0.0570235
Han, X., Zhou, Y., Zeng, Y., Sui, F., Liu, Y., Tan, Y., et al. (2017). Effects of active immunization against GnRH versus surgical castration on hypothalamic-pituitary function in boars. Theriogenology 97, 89–97. doi:10.1016/j.theriogenology.2017.04.038
Hennig, S. L., Owen, J. R., Lin, J. C., Young, A. E., Ross, P. J., Van Eenennaam, A. L., et al. (2020). Evaluation of mutation rates, mosaicism and off target mutations when injecting Cas9 mRNA or protein for genome editing of bovine embryos. Sci. Rep. 10, 22309. doi:10.1038/s41598-020-78264-8
Herbison, A. E., Porteous, R., Pape, J.-R., Mora, J. M., and Hurst, P. R. (2008). Gonadotropin-releasing hormone neuron requirements for puberty, ovulation, and fertility. Endocrinology 149, 597–604. doi:10.1210/en.2007-1139
Höijer, I., Emmanouilidou, A., Östlund, R., van Schendel, R., Bozorgpana, S., Tijsterman, M., et al. (2022). CRISPR-Cas9 induces large structural variants at on-target and off-target sites in vivo that segregate across generations. Nat. Commun. 13, 627. doi:10.1038/s41467-022-28244-5
Ikegami, K., Goto, T., Nakamura, S., Watanabe, Y., Sugimoto, A., Majarune, S., et al. (2020). Conditional kisspeptin neuron-specific Kiss1 knockout with newly generated Kiss1-floxed and Kiss1-Cre mice replicates a hypogonadal phenotype of global Kiss1 knockout mice. J. Reprod. Dev. 66, 359–367. doi:10.1262/jrd.2020-026
Jacyno, E., Kawęcka, M., Pietruszka, A., and Sosnowska, A. (2015). Phenotypic correlations of testes size with semen traits and the productive traits of young boars. Reprod. Dom. Anim. 50, 926–930. doi:10.1111/rda.12610
Kadokawa, H., Matsui, M., Hayashi, K., Matsunaga, N., Kawashima, C., Shimizu, T., et al. (2008). Peripheral administration of kisspeptin-10 increases plasma concentrations of GH as well as LH in prepubertal Holstein heifers. J. Endocrinol. 196, 331–334. doi:10.1677/JOE-07-0504
Kesner, J. S., Kraeling, R. R., Rampacek, G. B., and Johnson, B. (1987). Absence of an estradiol-induced surge of luteinizing hormone in pigs receiving unvarying pulsatile gonadotropin-releasing hormone stimulation. Endocrinology 121, 1862–1869. doi:10.1210/endo-121-5-1862
Kobayashi, T., Yamaguchi, T., Hamanaka, S., Kato-Itoh, M., Yamazaki, Y., Ibata, M., et al. (2010). Generation of rat pancreas in mouse by interspecific blastocyst injection of pluripotent stem cells. Cell 142, 787–799. doi:10.1016/j.cell.2010.07.039
Koppes, E. A., Redel, B. K., Johnson, M. A., Skvorak, K. J., Ghaloul-Gonzalez, L., Yates, M. E., et al. (2020). A porcine model of phenylketonuria generated by CRISPR/Cas9 genome editing. JCI Insight 5, e141523. doi:10.1172/jci.insight.141523
Korablev, A., Lukyanchikova, V., Serova, I., and Battulin, N. (2020). On-target CRISPR/Cas9 activity can cause undesigned large deletion in mouse zygotes. Int. J. Mol. Sci. 21, 3604. doi:10.3390/ijms21103604
Kurtz, S., Lucas-Hahn, A., Schlegelberger, B., Göhring, G., Niemann, H., Mettenleiter, T. C., et al. (2021). Knockout of the HMG domain of the porcine SRY gene causes sex reversal in gene-edited pigs. Proc. Natl. Acad. Sci. U. S. A. 118, e2008743118. doi:10.1073/pnas.2008743118
Lapatto, R., Pallais, J. C., Zhang, D., Chan, Y.-M., Mahan, A., Cerrato, F., et al. (2007). Kiss1−/− mice exhibit more variable hypogonadism than Gpr54−/− mice. Endocrinology 148, 4927–4936. doi:10.1210/en.2007-0078
Lee, W.-Y., Lee, K.-H., Heo, Y.-T., Kim, N.-H., Kim, J.-H., Kim, J.-H., et al. (2014). Transcriptional coactivator undifferentiated embryonic cell transcription factor 1 expressed in spermatogonial stem cells: A putative marker of boar spermatogonia. Anim. Reprod. Sci. 150, 115–124. doi:10.1016/j.anireprosci.2014.09.010
Lents, C. A., Heidorn, N. L., Barb, C. R., and Ford, J. J. (2008). Central and peripheral administration of kisspeptin activates gonadotropin but not somatotropin secretion in prepubertal gilts. Reproduction 135, 879–887. doi:10.1530/REP-07-0502
Lents, C. A. (2019). Review: kisspeptin and reproduction in the pig. Animal 13, 2986–2999. doi:10.1017/S1751731119001666
Lutz, J. B., Rampacek, G. B., Kraeling, R. R., and Pinkert, C. A. (1984). Serum luteinizing hormone and estrogen profiles before puberty in the gilt. J. Anim. Sci. 58, 686–691. doi:10.2527/jas1984.583686x
Malmgren, L., Rodriguez-Martinez, H., and Einarsson, S. (1996). Attainment of spermatogenesis in Swedish cross-bred boars. Zentralbl. Veterinarmed. A 43, 169–179. doi:10.1111/j.1439-0442.1996.tb00442.x
Mattsson, P. (2011). Prevalence of congenital defects in Swedish Hampshire, Landrace and Yorkshire pig breeds and opinions on their prevalence in Swedish commercial herds. Uppsala: SLU. Dept. of Animal Breeding and Genetics. Second cycle, A2E.
Mehravar, M., Shirazi, A., Nazari, M., and Banan, M. (2019). Mosaicism in CRISPR/Cas9-mediated genome editing. Dev. Biol. 445, 156–162. doi:10.1016/j.ydbio.2018.10.008
Menchaca, A., dos Santos-Neto, P. C., Souza-Neves, M., Cuadro, F., Mulet, A. P., Tesson, L., et al. (2020). Otoferlin gene editing in sheep via CRISPR-assisted ssODN-mediated Homology Directed Repair. Sci. Rep. 10, 5995. doi:10.1038/s41598-020-62879-y
Nagano, R., Tabata, S., Nakanishi, Y., Ohsako, S., Kurohmaru, M., and Hayashi, Y. (2000). Reproliferation and relocation of mouse male germ cells (gonocytes) during prespermatogenesis. Anat. Rec. 258, 210–220. doi:10.1002/(SICI)1097-0185(20000201)258:2<210::AID-AR10>3.0.CO;2-X
Owens, D. D. G., Caulder, A., Frontera, V., Harman, J. R., Allan, A. J., Bucakci, A., et al. (2019). Microhomologies are prevalent at Cas9-induced larger deletions. Nucleic Acids Res. 47, 7402–7417. doi:10.1093/nar/gkz459
Park, J., Lim, K., Kim, J.-S., and Bae, S. (2017). Cas-analyzer: an online tool for assessing genome editing results using NGS data. Bioinformatics 33, 286–288. doi:10.1093/bioinformatics/btw561
Rault, J.-L., Lay, D. C., and Marchant-Forde, J. N. (2011). Castration induced pain in pigs and other livestock. Appl. Anim. Behav. Sci. 135, 214–225. doi:10.1016/j.applanim.2011.10.017
Rydhmer, L., Zamaratskaia, G., Andersson, H. K., Algers, B., Guillemet, R., and Lundström, K. (2006). Aggressive and sexual behaviour of growing and finishing pigs reared in groups, without castration. Acta Agric. Scand. Sect. A - Animal Sci. 56, 109–119. doi:10.1080/09064700601079527
Schinckel, A. P., Johnson, R. K., and Kittok, R. J. (1984). Relationships among measures of testicular development and endocrine function in boars. J. Anim. Sci. 58, 1255–1261. doi:10.2527/jas1984.5851255x
Seminara, S. B., Chatzidaki, E. E., Thresher, R. R., Acierno, J. S., Kuohung, W., Zahn, D., et al. (2003). The GPR54 gene as a regulator of puberty. N. Engl. J. Med. 349, 1614–1627. doi:10.1056/NEJMoa035322
Semple, R. K., Achermann, J. C., Ellery, J., Farooqi, I. S., Karet, F. E., Stanhope, R. G., et al. (2005). Two novel missense mutations in G protein-coupled receptor 54 in a patient with hypogonadotropic hypogonadism. J. Clin. Endocrinol. Metab. 90, 1849–1855. doi:10.1210/jc.2004-1418
Sinclair, P. A., Squires, E. J., Raeside, J. I., Britt, J. H., and Hedgpeth, V. G. (2001). The effect of early postnatal treatment with a gonadotropin-releasing hormone agonist on the developmental profiles of testicular steroid hormones in the intact male pig. J. Anim. Sci. 79, 1003–1010. doi:10.2527/2001.7941003x
Sobrino, V., Avendaño, M. S., Perdices-López, C., Jimenez-Puyer, M., and Tena-Sempere, M. (2022). Kisspeptins and the neuroendocrine control of reproduction: Recent progress and new frontiers in kisspeptin research. Front. Neuroendocrinol. 65, 100977. doi:10.1016/j.yfrne.2021.100977
Sonstegard, T. S., Carlson, D., Lancto, C., and Fahrenkrug, S. C. (2016). Precision animal breeding as a sustainable, non-GMO solution for improving animal production and welfare. 31, 316–317. Corpus ID: 37749692.
Sonstegard, T. S., Fahrenkrug, S. C., and Carlson, D. (2017). 307 Precision animal breeding to make genetically castrated animals for improved animal welfare and alternative breeding applications. J. Animal Sci. 95, 149–150. doi:10.2527/asasmw.2017.307
Squires, E. J., Bone, C., and Cameron, J. (2020). Pork production with entire males: directions for control of boar taint. Animals. 10, 1665. doi:10.3390/ani10091665
Squires, E. J., Gray, M. A., and Lou, Y. (2019). Effect of mutations in porcine CYB5A and CYP17A1 on the metabolism of pregnenolone. J. Steroid Biochem. Mol. Biol. 195, 105469. doi:10.1016/j.jsbmb.2019.105469
Tan, W., Carlson, D. F., Lancto, C. A., Garbe, J. R., Webster, D. A., Hackett, P. B., et al. (2013). Efficient nonmeiotic allele introgression in livestock using custom endonucleases. Proc. Natl. Acad. Sci. U. S. A. 110, 16526–16531. doi:10.1073/pnas.1310478110
Tanihara, F., Hirata, M., Nguyen, N. T., Le, Q. A., Hirano, T., and Otoi, T. (2019). Effects of concentration of CRISPR/Cas9 components on genetic mosaicism in cytoplasmic microinjected porcine embryos. J. Reprod. Dev. 65, 209–214. doi:10.1262/jrd.2018-116
Tanihara, F., Hirata, M., Nguyen, N. T., Sawamoto, O., Kikuchi, T., Doi, M., et al. (2020). Efficient generation of GGTA1-deficient pigs by electroporation of the CRISPR/Cas9 system into in vitro-fertilized zygotes. BMC Biotechnol. 20, 40–11. doi:10.1186/s12896-020-00638-7
Trout, W. E., Killen, J. H., Christenson, R. K., Schanbacher, B. D., and Ford, J. J. (1992). Effects of weaning on concentrations of inhibin in follicular fluid and plasma of sows. J. Reprod. Fertil. 94, 107–114. doi:10.1530/jrf.0.0940107
Uenoyama, Y., Inoue, N., Nakamura, S., and Tsukamura, H. (2019). Central mechanism controlling pubertal onset in mammals: A triggering role of Kisspeptin. Front. Endocrinol. 10, 312. doi:10.3389/fendo.2019.00312
Uenoyama, Y., Nakamura, S., Hayakawa, Y., Ikegami, K., Watanabe, Y., Deura, C., et al. (2015). Lack of pulse and surge modes and glutamatergic stimulation of luteinising hormone release in Kiss1 knockout rats. J. Neuroendocrinol. 27, 187–197. doi:10.1111/jne.12257
Wen, W., Quan, Z.-J., Li, S.-A., Yang, Z.-X., Fu, Y.-W., Zhang, F., et al. (2021). Effective control of large deletions after double-strand breaks by homology-directed repair and dsODN insertion. Genome Biol. 22, 236–322. doi:10.1186/s13059-021-02462-4
Wen, W., and Zhang, X.-B. (2022). CRISPR–Cas9 gene editing induced complex on-target outcomes in human cells. Exp. Hematol. 110, 13–19. doi:10.1016/j.exphem.2022.03.002
Whitlock, B. K., Daniel, J. A., Wilborn, R. R., Rodning, S. P., Maxwell, H. S., Steele, B. P., et al. (2008). Interaction of estrogen and progesterone on kisspeptin-10-stimulated luteinizing hormone and growth hormone in ovariectomized cows. Neuroendocrinology 88, 212–215. doi:10.1159/000146242
Yoshioka, K., Suzuki, C., Tanaka, A., Anas, I. M.-K., and Iwamura, S. (2002). Birth of piglets derived from porcine zygotes cultured in a chemically defined medium. Biol. Reprod. 66, 112–119. doi:10.1095/biolreprod66.1.112
Young, L. D., Leymaster, K. A., and Lunstra, D. D. (1986). Genetic variation in testicular development and its relationship to female reproductive traits in swine. J. Anim. Sci. 63, 17–26. doi:10.2527/jas1986.63117x
Zamaratskaia, G., and Squires, E. J. (2009). Biochemical, nutritional and genetic effects on boar taint in entire male pigs. Animal 3, 1508–1521. doi:10.1017/S1751731108003674
Zhou, X., Wang, L., Du, Y., Xie, F., Li, L., Liu, Y., et al. (2015). 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. doi:10.1002/humu.22913
Keywords: kisspeptin, pig puberty, animal welfare, boar taint, embryo editing, homology-directed repair, knockout
Citation: Flórez JM, Martins K, Solin S, Bostrom JR, Rodríguez-Villamil P, Ongaratto F, Larson SA, Ganbaatar U, Coutts AW, Kern D, Murphy TW, Kim E-S, Carlson DF, Huisman A, Sonstegard TS and Lents CA (2023) CRISPR/Cas9-editing of KISS1 to generate pigs with hypogonadotropic hypogonadism as a castration free trait. Front. Genet. 13:1078991. doi: 10.3389/fgene.2022.1078991
Received: 24 October 2022; Accepted: 05 December 2022;
Published: 04 January 2023.
Edited by:
Bjoern Petersen, Friedrich Loeffler Institute (FLI), GermanyReviewed by:
Martina Crispo, Institut Pasteur de Montevideo, UruguayJoaquin Gadea, University of Murcia, Spain
Copyright © 2023 Flórez, Martins, Solin, Bostrom, Rodríguez-Villamil, Ongaratto, Larson, Ganbaatar, Coutts, Kern, Murphy, Kim, Carlson, Huisman, Sonstegard and Lents. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tad S. Sonstegard, tad@acceligen.com