Jenna Kropp Schmidt1*
Yun Hee Kim1
Nick Strelchenko1Sarah R. Gierczic1Derek Pavelec2
Thaddeus G. Golos1,3,4†
Igor I. Slukvin1,5,6†- 1Wisconsin National Primate Research Center, University of Wisconsin-Madison, Madison, WI, United States
- 2University of Wisconsin Biotechnology Center, University of Wisconsin-Madison, Madison, WI, United States
- 3Department of Comparative Biosciences, University of Wisconsin-Madison, Madison, WI, United States
- 4Department of Obstetrics and Gynecology, University of Wisconsin-Madison, Madison, WI, United States
- 5Department of Pathology and Laboratory Medicine, University of Wisconsin-Madison, Madison, WI, United States
- 6Department of Cell and Regenerative Biology, University of Wisconsin-Madison, Madison, WI, United States
Introduction: Genome editing by CRISPR-Cas9 approaches offers promise for introducing or correcting disease-associated mutations for research and clinical applications. Nonhuman primates are physiologically closer to humans than other laboratory animal models, providing ideal candidates for introducing human disease-associated mutations to develop models of human disease. The incidence of large chromosomal anomalies in CRISPR-Cas9-edited human embryos and cells warrants comprehensive genotypic investigation of editing outcomes in primate embryos. Our objective was to evaluate on- and off-target editing outcomes in CCR5 CRISPR-Cas9-targeted Mauritian cynomolgus macaque embryos.
Methods: DNA isolated from individual blastomeres of two embryos, along with paternal and maternal DNA, was subjected to whole genome sequencing (WGS) analysis.
Results: Large deletions were identified in macaque blastomeres at the on-target site that were not previously detected using PCR-based methods. De novo mutations were also identified at predicted CRISPR-Cas9 off-target sites.
Discussion: This is the first report of WGS analysis of CRISPR-Cas9-targeted nonhuman primate embryonic cells, in which a high editing efficiency was coupled with the incidence of editing errors in cells from two embryos. These data demonstrate that comprehensive sequencing-based methods are warranted for evaluating editing outcomes in primate embryos, as well as any resultant offspring to ensure that the observed phenotype is due to the targeted edit and not due to unidentified off-target mutations.
1 Introduction
Advances in genome editing, particularly using CRISPR-Cas9 technology, have facilitated the introduction and correction of disease-associated mutations in animal and cell culture models. Non-human primates (NHPs) are superior for modeling human diseases as they share similar aspects of immune, neuro-, and reproductive physiology and are ideal for transplant and neurodevelopmental disorder research. The interest in creating NHP models of human disease has been augmented by the need to better define the etiology of a disease and for the development of treatments and therapeutics (Dray et al., 2018; Abbott et al., 2019; Moshiri et al., 2019; Tapmeier et al., 2021; Ozirmak et al., 2022). For example, resistance to human immunodeficiency virus (HIV) has been observed in human patients with a 32 base pair deletion in the CCR5 gene (CCR5-Δ32). CCR5 serves as an HIV co-receptor (Dean et al., 1996; Liu et al., 1996; Samson et al., 1996). Transplantation of hematopoietic stem cells (HSCs) containing the CCR5-Δ32 mutation to human HIV patients led to the cure of HIV infection in some but not all cases (Hütter et al., 2009; Allers et al., 2011; Henrich et al., 2014; Cummins et al., 2017). Generating NHPs with CCR5-deletions would aid in determining the mechanisms of HIV elimination following transplantation of allogeneic HSCs with CCR5 mutations and the development of clinical protocols for reproducible HIV cure (Schmidt et al., 2022b).
Genome editing approaches to create gene disruption in NHP embryos have been successful, yet evidence of CRISPR-Cas9-induced chromosomal anomalies in mammalian cells and embryos warrants further investigation of embryonic editing outcomes in primate embryos. CRISPR-Cas9 editing has been shown to result in large scale deletions (up to 6 kb) and whole chromosome loss leading to genomic instability in mouse embryonic stem cells (Kosicki et al., 2018) and embryos (Adikusuma et al., 2018; Papathanasiou et al., 2021). CRISPR-Cas9 editing in human embryos has resulted in the loss of the targeted allele (Zuccaro et al., 2020) and also segmental chromosome losses (Alanis-Lobato et al., 2021). Moreover, large scale deletions in human embryos at an off-target site were also observed (Zuccaro et al., 2020). Loss of heterozygosity surrounding the on-target site is another consequence of CRISPR-Cas9 targeting observed in human embryos (Zuccaro et al., 2020; Alanis-Lobato et al., 2021). Collectively, these studies have revealed undesired on- and off-target mutations that arise when using wild-type Cas9 for gene correction in human and mouse embryos.
We previously demonstrated CRISPR-Cas9 editing of CCR5 in Mauritian cynomolgus macaque (MCM, Macaca fascicularis) embryos using PCR-based methods to confirm successful targeting of the locus (Schmidt et al., 2020). The objective of the present study was to comprehensively evaluate on- and off-target editing in CRISPR-Cas9-edited MCM embryos using whole genome sequencing (WGS) methods to survey individual blastomeres. Molecular analysis revealed large-scale deletions contributing to greater mosaicism within individual embryos than was previously identified using PCR-based methods. Given that large-scale on- and off-target mutations might hinder establishment of a viable pregnancy, further optimization of macaque embryo editing to avoid targeting errors would be essential to facilitate generation of novel NHP models for human diseases.
2 Materials and methods
2.1 Animals
Methods for deriving the MCM embryos analyzed in the present study were previously reported (Schmidt et al., 2020). Parental DNA was obtained from a female (12 years) and male (6 years) MCM used in that study. All procedures were performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals and under the approval of the University of Wisconsin-Madison College of Letters and Sciences and Vice Chancellor’s Office for Research and Graduate Education Institutional Animal Care and Use Committee.
2.2 Isolation of CRISPR-Cas9 injected embryo DNA
In vitro fertilized Mauritian cynomolgus macaque embryos were produced as previously described (Schmidt et al., 2020). Briefly, one-cell stage embryos were microinjected with Cas9 complexed with two sgRNAs targeting exon 2 of the CCR5 gene to form the ribonucleoprotein (RNP). Embryos were cultured individually in a microwell of a CultureCoin MIRI-TL dish (Esco Medical, Denmark) containing 25 μL of Global medium overlaid with mineral oil and the culture dish was placed in a MIRI TL Time-Lapse incubator (Esco Medical) to monitor embryo development. Individual blastomeres from two embryos arrested at the 6-cell and 9-cell stage were isolated. The zona pellucida was removed by treatment with 1 mg/ml of activated pronase E (Sigma-Aldrich, cat no: P2730) and the embryo was washed in calcium- and magnesium-free PBS with .2% EDTA and 1 mg/ml human albumin (MP Biomedicals, cat no: 823051). Individual blastomeres were then dissociated by gently pipetting the embryo. DNA was amplified from individual blastomeres using a REPLI-G single cell kit (Qiagen, cat no: 150343).
2.3 PCR analysis of CCR5 on-target editing in individual blastomeres
PCR amplifications were performed as previously described to assess the targeted CCR5 region (Schmidt et al., 2020) using primers that either amplified short (613 bp) or long (2,925 bp) amplicons surrounding the predicted cut sites. PCR products were run on 1.2%–1.5% agarose gels at 120 V. On-target CCR5 editing was determined by visualizing either the expected wild-type (unmodified) CCR5 PCR amplicon size of 613 bp or a biallelic mutation producing a product of 415 bp. To evaluate large-scale deletions near the on-target site, a long-range CCR5 PCR was performed. PCR primer sequences and expected amplicon sizes are listed in Supplementary Table S1. PCR reactions were performed using the Q5 Hot start High-Fidelity DNA polymerase kit following manufacturer recommendations. Gel electrophoresis was performed using standard methods to visualize the amplicons.
2.4 Isolation of parental DNA
Blood draws from both the oocyte and semen donors of the in vitro-produced embryos were performed to obtain parental DNA from peripheral blood mononuclear cells. Genomic DNA was isolated from blood cells using a Quick-DNA Miniprep kit (Zymo Research, cat no: D3024).
2.5 DNA quality assessment
DNA quality was assessed at the University of Wisconsin Biotechnology Center’s by the NexGen DNA Sequencing Core using an Agilent Femto Pulse system (Agilent, Santa Clara, CA) to confirm a uniform yield of DNA product with the average product length of greater than 9.4 kb.
2.6 Whole genome sequencing and analysis
Whole genome sequencing was performed by the University of Wisconsin-Madison Biotechnology Center using the Illumina short-read platform and a NovaSeq 6000 instrument. Reads were trimmed to remove sequencing adapters and low quality base calls using the software skewer (Jiang et al., 2014) and then mapped to the Macaca fascicularis reference genome, M_fascicularis_5.0, using an Illumina Dynamic Read Analysis for GENomics (DRAGEN) Bio-IT platform version 3.7. Small variant and calling was performed using DRAGEN. Variants from control (parental) samples were used to filter and identify de novo mutations. Variant annotation was performed using SNEPeff tool that will predict synonymous or non-synonymous amino acid changes, gains or losses of start/stop codons, and frame shifts due to insertions or deletions. Structural analysis was performed using Parliament2 (Zarate et al., 2020) and only those called by at least two callers were included. Variants with lower quality that were filtered out as well as non-filtered variants are included as potential candidate mutations. De novo structural variants were those identified in blastomeres that were not present in the parental sequences. Short read sequencing is not ideally suited for calling structural variants, hence low quality-filtered out variants as well non-filtered variants are both included as potential candidates. Integrated Genomics Viewer software (https://software.broadinstitute.org/software/igv/download) was used to view the WGS data.
2.7 Off-target analysis
Potential CRISPR-Cas9 off-target regions were identified using the Cas-OFFinder tool (Bae et al., 2014) (http://www.rgenome.net/cas-offinder/) and allowing for three mismatches. Regions of interest were then evaluated in the WGS dataset to see if de novo mutations were present in individual blastomeres compared to the parental DNA. The presence of de novo mutations in three predicted off-target genes were assessed by Sanger sequencing of PCR amplicons containing the region of interest for individual blastomeres. DNA obtained from a wild-type cynomolgus macaque iPSC line was sequenced in parallel. PCR reactions were performed using the Q5 Hot start High-Fidelity DNA polymerase kit following manufacturer recommendations and the reactions were cleaned up using a Gel extraction and PCR clean up kit (IBI, cat no: IB47010). PCR primer sequences are listed in Supplementary Table S1. Sanger sequencing reactions were carried out by Functional Biosciences Inc., Madison, Wisconsin and the sequencing data was analyzed using the 4Peaks (https://nucleobytes.com/4peaks/index.html) application.
3 Results
3.1 WGS of individual blastomeres produces variable sequence coverage
To functionally delete CCR5 in macaque embryos, we designed gRNAs that would encompass a 24-bp deletion that has previously been shown to be essential for expressing CCR5 in non-human primates (Chen et al., 1998). Our previous cell-based editing experiments confirmed successful on-target editing with functional deletion of the CCR5 gene in both human (Kang et al., 2015) and macaque (D'Souza et al., 2022) iPSCs. A schematic diagram of the targeting region is shown in Figure 1A. In our initial report describing targeting of this region in MCM embryos, PCR-based methods were used to evaluate CRISPR-Cas9 targeting of the CCR5 locus (Schmidt et al., 2020). Two embryos were dissociated into individual blastomeres and DNA was isolated for PCR evaluation and single-cell WGS. PCR and gel electrophoresis revealed editing mosaicism in each embryo (Figures 1B,C), although PCR signal was undetected in one and three blastomeres from embryo 4 and 5, respectively (Figure 1B).
FIGURE 1. Evaluation of CCR5 editing in MCM embryos using PCR-based methods (A). Schematic of the CCR5 gene including gRNA targeting and the forward (F1) and reverse (R1) primer. Wild-type (WT) product is 613 bp in length, whereas a biallelic deletion mutation (MT) produces a 415 bp product. Dashed line with the wild-type sequence indicates the targeting region (B). Gel electrophoresis images of PCR products from blastomeres of embryos 4 and 5. A positive control reaction with DNA from an unmanipulated control embryo and a no template negative control (NC) were included. The colored dots above each lane indicate the editing outcome as indicated in 1C. The PCR and gel electrophoresis results were provided in our initial report describing CCR5 editing in MCM embryos (Schmidt et al., 2020) (C). Diagram summarizing the PCR-based editing outcome and biallelic editing efficiency within each embryo.
Single-cell DNA amplification and WGS was performed on DNA from all six blastomeres of embryo 4 and 8 of 9 blastomeres from embryo 5. In addition, DNA isolated from peripheral blood mononuclear cells from the sire and dam of the embryos was sequenced in parallel. Chromosomal coverage varied across individual blastomeres and chromosomes ranging from .81–77.77-fold coverage, whereas the parental sequence coverage was at a depth of ∼30x, as expected for somatic cells. A sequencing coverage of 30x is interpreted as the genome was sequenced ∼30 times. Figure 2A shows the mean coverage and range of sequence depth across chromosomes for each sample and illustrates the variability in coverage in blastomeres compared to parental DNA isolated from peripheral blood cells. The CCR5 gene resides on chromosome 2, a chromosome that greatly varied in sequence depth coverage across blastomeres (Figure 2A). Sequence coverages for each sample by chromosome are listed in Supplementary Table S2. Blastomere 5-2 had an atypical distribution of GC content and was excluded from all analyses. For the remaining blastomeres, the sequence coverage at on- and off-target regions was taken into consideration and when limited sequences were observed at the region of interest the WGS result was deemed as not conclusive.
FIGURE 2. WGS identification of genotypes not identified by PCR-based methods (A). The mean and range in sequencing coverage depth (x-fold coverage) are plotted for each blastomere and parental DNA. The mean and range for all chromosomes is shown in the upper graph and the lower graph shows chromosome 2; CCR5 resides on chromosome 2 (B). Comparison of CCR5 editing outcomes by PCR versus WGS analysis. Blue text indicates cells where WGS identified deletions not previously identified by PCR. HET: heterozygous, HOM: homozygous (C). Sequence coverage at the CCR5 targeting region in blastomeres 5-4 and 5-8 where homozygous/biallelic edits were observed by WGS. The WGS viewer software indicates potential deletions with red bars and in the alignment tracks of 5-8, these are present around the target region and were minimal to absent in the parental coverage map. The vertical black box indicates the expected deleted region between the gRNA target sites.
3.2 WGS identified additional on-target deletions
WGS confirmed the genotypes identified using PCR-based methods for most of the blastomeres that had detectable PCR signal (7 out of 9) and determined the genotype for one blastomere in which the CCR5 region could not be amplified by PCR (Figure 2B). Sequences that spanned the target region are indicative of wild-type sequences, whereas a deletion was inferred if there was a break in the sequence coverage. Representative examples of wild-type (WT), heterozygous (HET) and homozygous deletion (HOM DEL) genotypes as determined by WGS are illustrated in Supplementary Figure S1. WGS sequence coverage at the CCR5 targeting region also identified deletions not previously observed using PCR-based methods (Figure 2B). If there was poor sequence coverage at the target site, the WGS genotype could not be determined and was deemed inconclusive. Blastomere 5-4 was identified to be HET by PCR-based methods, however, when looking at the sequencing alignment, no sequences spanned the region between the gRNA sites indicating that the blastomere contained a biallelic deletion (Figure 2C). Regardless of the unexpected gel band pattern of 5-4, the HOM DEL was confirmed via Sanger sequencing of the amplicons isolated from the three lowest bands in the agarose gel; each amplicon contained the expected ∼200 bp deletion and no WT sequences were detected (Supplementary Material S2). PCR using the standard primer pair previously failed to identify the genotype of blastomere 5-8, whereas WGS coverage indicated deletions spanning the gRNA sites and revealed large-scale deletions that encompassed the PCR primer sequences (Figure 2C). Moreover, the depth of coverage was reduced by approximately half at the 5’ end, confirming that one allele contained a large-scale deletion upstream of the first gRNA site.
3.3 Identification and validation of large-scale deletions at the CCR5 targeting region
De novo structural variants that were not present in the parental DNA but were within individual blastomeres of both embryos 4 and 5 were identified by WGS (Table 1). As short read sequencing platforms are not ideally suited for identifying structural variants, those that were called in at least two variant callers are listed. The on-target deletion between the gRNA sites was identified as a variant in 4-4. Large-scale deletions that span the target region were identified in blastomeres 4-6, 5-4, and a similar deletion of ∼5.2 kb was seen in both 5-5 and 5-8. Several inversions were also identified and within each embryo an inversion was unique to a pair of blastomeres.
TABLE 1. De novo structural variants at the CCR5 locus identified in individual blastomeres.
The PCR primer sequences used for initially genotyping blastomeres were within the deleted regions, hence these deletions could not be identified using PCR-based methods. The positions of the large-scale deletions identified by WGS in blastomeres 4-6, 5-5, and 5-8 are illustrated in Figure 3A. PCR-based methods using primers flanking each deleted region followed by gel electrophoresis confirmed the presence of the large-scale deletions in 4-6, 5-5, and 5-8 (Figure 3B). The large-scale deletions were confirmed by Sanger sequencing of PCR amplicons for the deletions detected in 4-6, 5-5, and 5-8 with one exception (Supplementary Material S1). Poor sequencing data did not allow for verification of the ∼5.2 kb deletion in the 5.8 blastomere whereas the deleted sequence was confirmed in blastomere 5.5 (Supplementary Material S1).
FIGURE 3. Identification and verification of a large-scale deletion identified by WGS in blastomere 4-6 (A). Schematic diagram of deletions identified by WGS in blastomeres 4-6, 5-5, and 5-8, and the positions of the primers to evaluate the ∼756 bp deletion in blastomere 4-6. Primers are denoted as F or R and their sequences are listed in Supplementary Table S1 (B). Gel electrophoresis image of the PCR amplicons generated by primers that flank each deletion. A no template negative control (NC) reaction was also included.
3.4 WGS detection of on-target INDELs in blastomere 5-4
The presence of insertions and/or deletions (INDELs) within the on-target region was evaluated. Single nucleotide variants (SNVs) were identified and considered to be de novo mutations if they were present in the blastomere but not detected in either parent. The number of de novo SNVs detected in each blastomere is provided in Supplementary Table S3. Homozygous insertions of 13 bp, 5 bp, 11 bp, and 33 bp were identified near the gRNA one cut site in blastomere 5-4 that were not identified in the cynomolgus macaque reference genome nor the parental DNA sequences (Figure 4). This was the only blastomere with INDEL formation near a cut site.
FIGURE 4. Insertions in blastomere 5-4 at the CCR5 target site. De novo insertions identified in blastomere 5-4 that were not identified in the parental DNA. The gRNA 1 sequence is indicated by the dark gray bar and the blue indicates the PAM sequence (5′-TTG). Homozygous (HOM) variants (VAR), specifically insertions, were detected near the predicted DNA cut site (indicated by scissors) that were not detected in the cynomolgus macaque reference genome (REF) nor the parental DNA sequence.
3.5 Off-target edits detected at loci with restricted sequence homology to the gRNAs
To assess the feasibility of detecting off-target mutations by WGS, we utilized the in silico Cas-OFFinder tool to identify off-target sites based on sequence homology to the gRNA sequences and the total number of off-target sites when allowing for up to 9 mismatches between the gRNA and off-target sequence is provided in Supplementary Table S4. While off-target editing can occur at sites with up to 4 mismatches (Pattanayak et al., 2013; Cromer et al., 2022), we focused on assessing the WGS data at 93 predicted off-target regions that contained 2 or 3 mismatches which occurs more commonly. De novo mutations within individual blastomeres were identified by WGS in 16 predicted off-targeted sites of which 7 were located within genes and nine were located in intergenic regions (Table 2). The specific mutations identified by WGS at each off-target site are provided in Supplementary Material S2.
TABLE 2. De novo mutations identified by WGS in predicted off-target regions.
Blastomeres with mutations at off-target regions identified by WGS were subjected to PCR and Sanger sequencing, with the exception that all blastomeres of embryo 5 were sequenced for the SFMTB2 region. Sanger sequences are shown in Table 3. The T to G point mutation at the predicted region in the NFASC gene was identified in blastomeres of both embryos 4 and 5. This mutation is likely a spontaneously occurring mutation as the parental genotype was T/T, embryo 1 blastomeres were G/G and embryo 5 blastomeres carried either a T/T, G/G or T/G genotype. De novo deletions near the predicted cut site of the off-target regions in the SFMBT2 and LIPC genes suggests off-target editing by the CRISPR-Cas9 RNP (Table 3). Furthermore, blastomeres of both embryos show deletions in SFMBT2 and pairs of blastomeres in embryo 5 have the same 7 or 8 bp deletion suggestive that the editing occurred in a previous cell division similar to mutations introduced by on-target CRISPR-Cas9 targeting.
TABLE 3. Sanger-sequencing of amplicons containing predicted off-target mutations identified by WGS.
3.6 Structural variants detected at off-target sites
The presence of structural variants at predicted off-target sites was investigated for the off-target regions in the NFASC, SFMBT2, and LIPC genes. Large-scale deletions, inversions and duplications were identified at these sites by WGS (Table 4). An inversion and deletion at the NFASC off-target region was observed in each blastomere of embryos 4 and 5, and an inversion identified at the LIPC site was shared between blastomeres 4-5 and 4-6. A smaller deletion of 210 bp was detected in blastomere 4-1. Fewer off-target structural variants were shared among blastomeres within an embryo compared to those that were identified at the on-target site. A summary of on- and off-target mutations identified by WGS in each blastomere is provided in Supplementary Table S5.
TABLE 4. De novo structural variants at predicted off-target sites identified in individual blastomeres by WGS.
4 Discussion
In the present study, comprehensive assessment of CRISPR-Cas9 targeting in MCM blastomeres by WGS confirmed editing mosaicism, revealed undesired on- and off-target editing events in NHP embryos, including large scale deletions, and resolved genotypes at the on-target sites that were previously undetected using PCR-based methods. INDELs were observed at on- and predicted off-target sites, where sequence disruption was confirmed by Sanger sequencing for two off-target regions. WGS analysis also provided insight into the timing of CRISPR-Cas9 targeting as identical structural variants and de novo mutations were shared in pairs of blastomeres but were not identified in the majority of blastomeres suggesting that editing was delayed and did not occur at the one-cell stage. While CRISPR-Cas9 can introduce mutations at disease-associated loci in NHP embryos, the occurrence of unexpected editing events requires rigorous assessment of not only embryos, but also the offspring to confirm that any resulting phenotype is not due to off-target effects.
Undesired editing events at the on-target site, including large-scale deletions, have been observed in human and mouse embryos targeted with wild-type Cas9 nuclease (Adikusuma et al., 2018; Zuccaro et al., 2020; Alanis-Lobato et al., 2021; Papathanasiou et al., 2021). In the present study, the expected 198 bp deletion as well as 756 bp, 925 bp and ∼5.2 kb deletions were detected at the CCR5 on-target site through WGS structural variant analysis of individual NHP embryonic cells. Deletions have been previously described in reports of CRISPR-Cas9 targeting in NHP embryos that were transferred to surrogates and produced edited offspring, including an ∼11.5 kb deletion in SHANK3 in one cynomolgus monkey (Zhao et al., 2017), a ∼7.2 kb deletion in PINK1 in two rhesus monkeys (Yang et al., 2019), and 920 bp at the OCT4 knock-in site in a cynomolgus monkey (Cui et al., 2018). The SHANK3 mutant died in utero at 125 days of gestation (term is 165 days) and the two PINK1 mutants were triplets that died days after birth. These studies demonstrated that implantation and pregnancies can be achieved despite the presence of large-scale editing errors. Furthermore, on-target deletions have now been observed across studies in both NHP embryos and tissues where different genes were targeted by wild-type CRISPR-Cas9, which necessitates refinement of genome editing tools for creating precise disease-associated mutations.
In the current study, we identified pairs of blastomeres of the same embryo which contained similar structural variants, however these variants were not identical suggesting that they arose from separate editing events (e.g., blastomeres 5-5 and 5-8 shared a ∼5.2 kb deletion and 5-4 and 5-5 shared a ∼35,341 kb inversion near the CCR5 targeting site). Mosaicism in blastomere genotypes suggests that editing was delayed and did not occur at the one-cell stage as the deletion was not detected in all cells of the embryo. Editing mosaicism has been observed in tissues of edited NHPs produced from embryo transfer of CRISPR-Cas9 targeted embryos (Niu et al., 2014; Chen et al., 2015; Tsukiyama et al., 2019; Zhou et al., 2019).
A goal of the present study was to assess the feasibility of using WGS to assess off-target effects in individual blastomeres. An in silico based approach guided the nomination of potential off-target regions based on sequence homology to the gRNAs allowing up to three mismatches. Following CCR5 targeting, mutations were observed in the genes SFMBT2 and LIPC by WGS and the introduction of INDELS was confirmed by Sanger sequencing. Blastomeres of both embryos displayed sequence disruption via INDEL formation with the presence of 4, 7 or 9 bp deletions in some of the cells at the predicted off-target site within the SFMBT2 gene. These results confirmed that off-target editing could be assessed by WGS in individual blastomeres, although we used a biased in silico method that relied on assessing targets with sequence homology and did not evaluate potential targets with greater than three mismatches. Additional in silico nominated targets should be evaluated to fully assess the impact of off-target editing as CRISPR-Cas9 cleavage can occur at off-target sites with up to four mismatches (Pattanayak et al., 2013; Cromer et al., 2022). Moreover, unbiased methods that survey the whole genome without prior knowledge or prediction of sequence homology would be more informative, yet there is not a current superior method or technique for this analysis (Chaudhari et al., 2020; Atkins et al., 2021). In vitro-based off-target analysis methods in the future could be adapted for single-cells, but with current use of whole genome amplification (WGA) and an incomplete reference genome, a large number of false positives might be called due to errors incurred during WGA or due to differences in the reference assembly.
Relatively few off-target mutations have been identified in studies that have generated edited NHP offspring by transfer of CRISPR-Cas9 targeted embryos to surrogate embryo recipients. A 2 bp deletion in one off-target region was reported in an edited cynomolgus monkey (Cui et al., 2018) and one intronic and two intergenic INDELS were identified in two edited rhesus monkeys (Wang et al., 2018). Luo et al. (2019) identified de novo mutations that the authors thought were not introduced by CRISPR-Cas9 and rather could be attributed to natural spontaneous generational mutations or that were due to technical noise. In these previous NHP studies, WGS analysis was performed on DNA obtained from cells or tissues of live offspring or miscarried fetuses whereas here we reported WGS on individual blastomeres of NHP embryos. The higher incidence of off-target editing in this study could be explained by the editing efficiency of the RNP or it is possible that such significant off-target errors could have been embryonic lethal and therefore not present in offspring that survived to or near term. Editing has shown to be more rapid and efficient when targeting with an RNP versus Cas9 mRNA in NHP embryos (Midic et al., 2017). To mitigate potential off-target editing events, the following strategies could be implemented: 1) microinjection of the RNP at the time of fertilization (Lamas-Toranzo et al., 2019), 2) use of a Cas9 nuclease modified to improve specificity (Huang et al., 2022), or 3) use of a base or prime editing approach that does not result in a double-stranded DNA break (Zeballos and Gaj, 2020).
Off-target INDELs and segmental chromosome errors introduced by CRISPR-Cas9 have been observed in human embryos (Zuccaro et al., 2020). The formation of INDELs has been reported at a predicted off-target site that had two mismatches to the gRNA sequence (Zuccaro et al., 2020). Moreover, segmental chromosome errors were detected near predicted off-target sites and were often restricted to one cell, hence the authors concluded that the events likely occurred during the second or third cell cycle. While the present study did not focus on whole or segmental chromosomal errors, structural variants were detected at off-target sites that could be investigated in the future.
Large-scale mutations introduced by CRISPR-Cas9 editing in human and mouse embryos have shown to contribute to genomic instability through segmental and whole chromosomal loss (Adikusuma et al., 2018; Zuccaro et al., 2020; Alanis-Lobato et al., 2021; Papathanasiou et al., 2021). Unrepaired double-stranded DNA breaks at the CRISPR-Cas9 cleavage site have shown to result in fragmented chromosomes leading to chromosome mis-segregation and micronuclei formation in human cell lines and mouse embryos (Leibowitz et al., 2021; Papathanasiou et al., 2021). In human cleavage stage embryos, failure to replicate the genome before entry into mitotic divisions contributes to poorer embryo quality due to aneuploidy associated with chromosomal fragmentation and the formation of a micronucleus (Palmerola et al., 2022). Human and NHP embryos naturally have a higher incidence of aneuploidy where partial or whole chromosomes encapsulated by micronuclei may be present in fragmented cells of the embryos (Daughtry et al., 2019; Palmerola et al., 2022). Cellular fragmentation was previously observed in >65% of in vitro fertilized cleavage stage rhesus macaque embryos, and when analyzing individual blastomeres from 50 embryos, 73.5% showed whole and/or partial chromosomal losses or gains (Daughtry et al., 2019). The incidence of chromothripsis was not assessed in the present study given the limitations in interpreting WGA artifacts in light of an incomplete reference genome, in addition to the difficulty in discerning whether chromosome loss was due to CRISPR-Cas9 targeting or naturally-occurring cellular fragmentation events. As CRISPR-Cas9 editing errors can result in chromosomal disruption and/or elimination (Leibowitz et al., 2021; Papathanasiou et al., 2021) it is plausible that the embryo development may be negatively impacted and could explain the poor embryo transfer rate observed in our previous study (Schmidt et al., 2020) and in general the low live birth rates of CRISPR-Cas9 targeted NHPs (Schmidt et al., 2022a; Schmidt et al., 2022b).
4.1 Limitations of the study
Single-cell WGA can introduce amplification bias potentially limiting the interpretation of the variants identified in this study. Several studies comparing commercially available sc-WGA kits revealed differences in reproducibility, error rates, target coverage, read depth distribution and allele drop out, however, the REPLI-G sc-WGA kit used in the current study was shown to have a high mapping rage (>90%), be reproducible and have a lower error rate (Borgström et al., 2017; Biezuner et al., 2021) compared to other kits. In the present study, sequence coverage was variable across chromosomes within and across individual blastomeres compared to the parental DNA that did not undergo WGA. A limitation to the present study is the lack of analysis of unmanipulated control blastomeres to assess the rate of errors or artifacts incurred due to WGA. Translocation events were not called for this reason and only structural variants identified by two callers were reported in this study. While it is uncertain whether the structural variants identified by WGS are due to CRISPR-Cas9 targeting, the variants were detected in pairs of blastomeres that were processed through independent WGA reactions suggesting that the mutation occurred during an early cleavage division or that similar sites of the genome are reproducibly prone to amplification errors during the WGA process. Recent advances in WGA technology such as primary template-directed amplification (Gonzalez-Pena et al., 2021) or linear amplification through transposon insertion (Chen et al., 2017) have shown to amplify single-cell genomes with more uniformity and accuracy. While these technologies are not incorporated into a commercial kit, they could be implemented in future single blastomere WGS studies for greater accuracy and reproducibility.
5 Conclusion
Overall, utilizing a WGS approach to determine CRISPR-Cas9 editing outcomes allows for the identification of edits not identified by PCR. In this study, WGS revealed the incidence of on-target large-scale deletions and INDEL formation at off-target sites. Imprecise editing could hinder the development of an NHP disease model that both genocopies and phenocopies the disease. The consequences of undesired editing events on gene expression of off-target and/or neighboring genes was not evaluated in this study, but should be considered in future studies. Based on evidence from human and mouse embryonic targeting by CRISPR-Cas9, it is likely that chromosomal damage incurred early in embryo development could negatively impact embryo viability. It remains unclear whether a reduced concentration or volume of the RNP delivered to the one-cell embryo would have a dose-dependent impact on on-target errors. Additional studies are needed to optimize embryonic editing by wild-type Cas9 and/or to use alternative next-generation Cas9 nucleases that do not create a double-stranded DNA break (Komor et al., 2016; Zeballos and Gaj, 2020). Regardless, WGS analysis should be implemented to thoroughly characterize editing genotypes in NHP models generated through this technology.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA880597.
Ethics statement
The animal study was reviewed and approved by the University of Wisconsin-Madison College of Letters and Sciences and Vice Chancellor’s Office for Research and Graduate Education Institutional Animal Care and Use Committee.
Author contributions
JS, NS, TG, and IS contributed to the conception and design of the study. NS, YK, SG, and JS collected and analyzed the data. DP performed the bioinformatic analysis. JS, TG, and IS drafted the manuscript. The publication’s contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIH.
Funding
This research was funded by NIH grants R24OD021322 awarded to IS and TG, P51OD011106 awarded to the Wisconsin National Primate Research Center, University of Wisconsin-Madison, and K99 HD099154-01 awarded to JS. This research was conducted at a facility constructed with support from Research Facilities Improvement Program grant numbers RR15459-01 and RR020141-01.
Acknowledgments
The authors would like to thank Dr. Jeffrey Rogers for his guidance in the analysis and for his insightful review of the manuscript. We would also like to acknowledge Logan Keding for his contribution to the graphic design of Figure 1.
Conflict of interest
The 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgeed.2022.1031275/full#supplementary-material
Abbreviations
WGS, Whole genome sequencing; NHP, Non-human primate; HIV, Human immunodeficiency virus; HSC, Hematopoietic stem cell; MCM, Mauritian cynomolgus macaque; RNP, Ribonucleoprotein; CNV, copy number variation; WT, Wild-type; HET, Heterozygous; HOM, Homozygous; DEL, Deletion; InDel, Insertions and/or deletions; SNVs, Single nucleotide variants.
References
Abbott, D., Rogers, J., Dumesic, D., and Levine, J. (2019). Naturally occurring and experimentally induced rhesus macaque models for polycystic ovary syndrome: Translational gateways to clinical application. Med. Sci. (Basel, Switz. 7, 107. doi:10.3390/medsci7120107
Adikusuma, F., Adikusuma, F., Piltz, S., and Corbett, M. A. (2018). Large deletions induced by Cas9 cleavage. Nature 560, E8–E9. doi:10.1038/s41586-018-0380-z
Alanis-Lobato, G., Zohren, J., Mccarthy, A., Fogarty, N., 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
Allers, K., Hütter, G., Hofmann, J., Loddenkemper, C., Rieger, K., Thiel, E., et al. (2011). Evidence for the cure of HIV infection by CCR5Δ32/Δ32 stem cell transplantation. Blood 117, 2791–2799. doi:10.1182/blood-2010-09-309591
Atkins, A., Chung, C., Allen, A., Dampier, W., Gurrola, T., Sariyer, I., et al. (2021). Off-target analysis in gene editing and applications for clinical translation of CRISPR/Cas9 in HIV-1 therapy. Front. genome Ed. 3, 673022. doi:10.3389/fgeed.2021.673022
Bae, S., Park, J., and Kim, J. (2014). Cas-OFFinder: A fast and versatile algorithm that searches for potential off-target sites of Cas9 RNA-guided endonucleases. Bioinforma. Oxf. Engl. 30, 1473–1475. doi:10.1093/bioinformatics/btu048
Biezuner, T., Raz, O., Amir, S., Milo, L., Ada, R., Fried, Y., et al. (2021). Comparison of seven single cell whole genome amplification commercial kits using targeted sequencing. Sci. Rep. 11, 17171. doi:10.1038/s41598-021-96045-9
Borgström, E., Paterlini, M., Mold, J., Frisen, J., and Lundeberg, J. (2017). Comparison of whole genome amplification techniques for human single cell exome sequencing. PloS one 12, e0171566. doi:10.1371/journal.pone.0171566
Chaudhari, H., Penterman, J., Whitton, H., Spencer, S., Flanagan, N., Lei Zhang, M., et al. (2020). Evaluation of homology-independent CRISPR-cas9 off-target assessment methods. CRISPR J. 3, 440–453. doi:10.1089/crispr.2020.0053
Chen, C., Xing, D., Tan, L., Li, H., Zhou, G., Huang, L., et al. (2017). Single-cell whole-genome analyses by linear amplification via transposon insertion (LIANTI). Sci. (New York, N.Y.) 356, 189–194. doi:10.1126/science.aak9787
Chen, Y., Zheng, Y., Kang, Y., Yang, W., Niu, Y., Guo, X., et al. (2015). Functional disruption of the dystrophin gene in rhesus monkey using CRISPR/Cas9. Hum. Mol. Genet. 24, 3764–3774. doi:10.1093/hmg/ddv120
Chen, Z., Kwon, D., Jin, Z., Monard, S., Telfer, P., Jones, M., et al. (1998). Natural infection of a homozygous delta24 CCR5 red-capped mangabey with an R2b-tropic simian immunodeficiency virus. J. Exp. Med. 188, 2057–2065. doi:10.1084/jem.188.11.2057
Cromer, M., Majeti, K., Rettig, G., Murugan, K., Kurgan, G., Hampton, J., et al. (2022). Comparative analysis of CRISPR off-target activity discovery tools following ex vivo editing of CD34+ hematopoietic stem and progenitor cells. bioRxiv.
Cui, Y., Niu, Y., Zhou, J., Chen, Y., Cheng, Y., Li, S., et al. (2018). Generation of a precise Oct4-hrGFP knockin cynomolgus monkey model via CRISPR/Cas9-assisted homologous recombination. Cell Res. 28, 383–386. doi:10.1038/cr.2018.10
Cummins, N., Rizza, S., Litzow, M., Hua, S., Lee, G., Einkauf, K., et al. (2017). Extensive virologic and immunologic characterization in an HIV-infected individual following allogeneic stem cell transplant and analytic cessation of antiretroviral therapy: A case study. PLoS Med. 14, e1002461. doi:10.1371/journal.pmed.1002461
D'Souza, S., Kumar, A., Weinfurter, J., Park, M., Maufort, J., Tao, L., et al. (2022). Generation of SIV-resistant T cells and macrophages from nonhuman primate induced pluripotent stem cells with edited CCR5 locus. Stem Cell Rep. 17, 953–963. doi:10.1016/j.stemcr.2022.03.003
Daughtry, B., Rosenkrantz, J., Lazar, N., Fei, S., Redmayne, N., Torkenczy, K., et al. (2019). Single-cell sequencing of primate preimplantation embryos reveals chromosome elimination via cellular fragmentation and blastomere exclusion. Genome Res. 29, 367–382. doi:10.1101/gr.239830.118
Dean, M., Carrington, M., Winkler, C., Huttley, G., Smith, M., Allikmets, R., et al. (1996). Genetic restriction of HIV-1 infection and progression to AIDS by a deletion allele of the CKR5 structural gene. Hemophilia growth and development study, multicenter AIDS cohort study, multicenter hemophilia cohort study, san francisco city cohort, ALIVE study. Sci. (New York, N.Y.) 273, 1856–1862. doi:10.1126/science.273.5283.1856
Dray, B., Raveendran, M., Harris, R., Benavides, F., Gray, S., Perez, C., et al. (2018). Mismatch repair gene mutations lead to lynch syndrome colorectal cancer in rhesus macaques. Genes & cancer 9, 142–152. doi:10.18632/genesandcancer.170
Gonzalez-Pena, V., Natarajan, S., Xia, Y., Klein, D., Carter, R., Pang, Y., et al. (2021). Accurate genomic variant detection in single cells with primary template-directed amplification. Proc. Natl. Acad. Sci. U. S. A. 118, e2024176118. doi:10.1073/pnas.2024176118
Henrich, T., Hanhauser, E., Marty, F., Sirignano, M., Keating, S., Lee, T., et al. (2014). Antiretroviral-free HIV-1 remission and viral rebound after allogeneic stem cell transplantation: Report of 2 cases. Ann. Intern. Med. 161, 319–327. doi:10.7326/M14-1027
Huang, X., Yang, D., Zhang, J., Xu, J., and Chen, Y. (2022). Recent advances in improving gene-editing specificity through CRISPR-cas9 nuclease engineering. Cells 11, 2186. doi:10.3390/cells11142186
Hütter, G., Nowak, D., Mossner, M., Ganepola, S., Müssig, A., Allers, K., et al. (2009). Long-term control of HIV by CCR5 Delta32/Delta32 stem-cell transplantation. N. Engl. J. Med. 360, 692–698. doi:10.1056/NEJMoa0802905
Jiang, H., Lei, R., Ding, S. W., and Zhu, S. (2014). Skewer: A fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinforma. 15, 182. doi:10.1186/1471-2105-15-182
Kang, H., Minder, P., Park, M., Mequitta, W., Torbett, B., and Slukvin, I. I. (2015). CCR5 disruption in induced pluripotent stem cells using CRISPR/Cas9 provides selective resistance of immune cells to CCR5-tropic HIV-1 virus. Mol. Ther. Nucleic Acids 4, e268. doi:10.1038/mtna.2015.42
Komor, A. C., Kim, Y. B., Packer, M. S., Zuris, J. A., and Liu, D. R. (2016). Programmable editing of a target base in genomic DNA without double-stranded DNA cleavage. Nature 533, 420–424. doi:10.1038/nature17946
Kosicki, M., Tomberg, K., and Bradley, A. (2018). Repair of double-strand breaks induced by CRISPR-Cas9 leads to large deletions and complex rearrangements. Nat. Biotechnol. 36, 765–771. doi:10.1038/nbt.4192
Lamas-Toranzo, I., Galiano-Cogolludo, B., Cornudella-Ardiaca, F., Cobos-Figueroa, J., Ousinde, O., and Bermejo-Álvarez, P. (2019). Strategies to reduce genetic mosaicism following CRISPR-mediated genome edition in bovine embryos. Sci. Rep. 9, 14900. doi:10.1038/s41598-019-51366-8
Leibowitz, M. L., Papathanasiou, S., Doerfler, P., Blaine, L., Sun, L., Yao, Y., et al. (2021). Chromothripsis as an on-target consequence of CRISPR-Cas9 genome editing. Nat. Genet. 53, 895–905. doi:10.1038/s41588-021-00838-7
Liu, R., Paxton, W., Choe, S., Ceradini, D., Martin, S., Horuk, R., et al. (1996). Homozygous defect in HIV-1 coreceptor accounts for resistance of some multiply-exposed individuals to HIV-1 infection. Cell 86, 367–377. doi:10.1016/s0092-8674(00)80110-5
Luo, X., He, Y., Zhang, C., He, X., Yan, L., Li, M., et al. (2019). Trio deep-sequencing does not reveal unexpected off-target and on-target mutations in Cas9-edited rhesus monkeys. Nat. Commun. 10, 5525. doi:10.1038/s41467-019-13481-y
Midic, U., Hung, P. H., Vincent, K. A., Goheen, B., Schupp, P. G., Chen, D. D., et al. (2017). Quantitative assessment of timing, efficiency, specificity and genetic mosaicism of CRISPR/Cas9-mediated gene editing of hemoglobin beta gene in rhesus monkey embryos. Hum. Mol. Genet. 26, 2678–2689. doi:10.1093/hmg/ddx154
Moshiri, A., Chen, R., Kim, S., Harris, R. A., Li, Y. M., Raveendran, M., et al. (2019). A nonhuman primate model of inherited retinal disease. J. Clin. Investigation 129, 863–874. doi:10.1172/JCI123980
Niu, Y., Shen, B., Cui, Y., Chen, Y., Wang, J., Wang, L., et al. (2014). Generation of gene-modified cynomolgus monkey via Cas9/RNA-mediated gene targeting in one-cell embryos. Cell 156, 836–843. doi:10.1016/j.cell.2014.01.027
Ozirmak, L. N., Gray, S., Bowen, C., Reyes-Uribe, L., Dray, B., Deng, N., et al. (2022). Comparative molecular genomic analyses of a spontaneous rhesus macaque model of mismatch repair-deficient colorectal cancer. PLoS Genet. 18, e1010163. doi:10.1371/journal.pgen.1010163
Palmerola, K., Amrane, S., De Los Angeles, A., Xu, S., Wang, N., De Pinho, J., et al. (2022). Replication stress impairs chromosome segregation and preimplantation development in human embryos. Cell 185, 2988–3007.e20. doi:10.1016/j.cell.2022.06.028
Papathanasiou, S., Markoulaki, S., Blaine, L., Leibowitz, M., Zhang, C., Jaenisch, R., et al. (2021). Whole chromosome loss and genomic instability in mouse embryos after CRISPR-Cas9 genome editing. Nat. Commun. 12, 5855. doi:10.1038/s41467-021-26097-y
Pattanayak, V., Lin, S., Guilinger, J., Ma, E., Doudna, J., and Liu, D. (2013). High-throughput profiling of off-target DNA cleavage reveals RNA-programmed Cas9 nuclease specificity. Nat. Biotechnol. 31, 839–843. doi:10.1038/nbt.2673
Samson, M., Libert, F., Doranz, B., Rucker, J., Liesnard, C., Farber, C., et al. (1996). Resistance to HIV-1 infection in caucasian individuals bearing mutant alleles of the CCR-5 chemokine receptor gene. Nature 382, 722–725. doi:10.1038/382722a0
Schmidt, J., Jones, K., Van Vleck, T., and Emborg, M. (2022a). Modeling genetic diseases in nonhuman primates through embryonic and germline modification: Considerations and challenges. Sci. Transl. Med. 14, eabf4879. doi:10.1126/scitranslmed.abf4879
Schmidt, J. K., Strelchenko, N., Park, M. A., Kim, Y. H., Mean, K. D., Schotzko, M. L., et al. (2020). Genome editing of CCR5 by CRISPR-Cas9 in Mauritian cynomolgus macaque embryos. Sci. Rep. 10, 18457. doi:10.1038/s41598-020-75295-z
Schmidt, J., Reynolds, M., Golos, T., and Slukvin, I. (2022b). CRISPR/Cas9 genome editing to create nonhuman primate models for studying stem cell therapies for HIV infection. Retrovirology 19, 17. doi:10.1186/s12977-022-00604-5
Tapmeier, T., Rahmioglu, N., Lin, J., De Leo, B., Obendorf, M., Raveendran, M., et al. (2021). Neuropeptide S receptor 1 is a nonhormonal treatment target in endometriosis. Sci. Transl. Med. 13, eabd6469. doi:10.1126/scitranslmed.abd6469
Tsukiyama, T., Kobayashi, K., Nakaya, M., Iwatani, C., Seita, Y., Tsuchiya, H., et al. (2019). Monkeys mutant for PKD1 recapitulate human autosomal dominant polycystic kidney disease. Nat. Commun. 10, 5517. doi:10.1038/s41467-019-13398-6
Wang, S., Ren, S., Bai, R., Xiao, P., Zhou, Q., Zhou, Y., et al. (2018). No off-target mutations in functional genome regions of a CRISPR/Cas9-generated monkey model of muscular dystrophy. J. Biol. Chem. 293, 11654–11658. doi:10.1074/jbc.AC118.004404
Yang, W., Liu, Y., Tu, Z., Xiao, C., Yan, S., Ma, X., et al. (2019). CRISPR/Cas9-mediated PINK1 deletion leads to neurodegeneration in rhesus monkeys. Cell Res. 29, 334–336. doi:10.1038/s41422-019-0142-y
Zarate, S., Carrol, A., Mahmoud, M., Krasheninina, O., Jun, G., Salerno, W. J., et al. (2020). Parliament2: Accurate structural variant calling at scale. GigaScience 9, giaa145. giaa145. doi:10.1093/gigascience/giaa145
Zeballos, M. A., and Gaj, T. (2020). Next-generation CRISPR technologies and their applications in gene and cell therapy. Trends Biotechnol. 39, 692–705.
Zhao, H., Tu, Z., Xu, H., Yan, S., Yan, H., Zheng, Y., et al. (2017). Altered neurogenesis and disrupted expression of synaptic proteins in prefrontal cortex of SHANK3-deficient non-human primate. Cell Res. 27, 1293–1297. doi:10.1038/cr.2017.95
Zhou, Y., Sharma, J., Ke, Q., Landman, R., Yuan, J., Chen, H., et al. (2019). Atypical behavior and connectivity in SHANK3-mutant macaques. Nature 570, 326–331. doi:10.1038/s41586-019-1278-0
Keywords: CRISPR- Cas9, whole genome sequencing (WGS), embryo, macaque, CCR5
Citation: Schmidt JK, Kim YH, Strelchenko N, Gierczic SR, Pavelec D, Golos TG and Slukvin II (2023) Whole genome sequencing of CCR5 CRISPR-Cas9-edited Mauritian cynomolgus macaque blastomeres reveals large-scale deletions and off-target edits. Front. Genome Ed. 4:1031275. doi: 10.3389/fgeed.2022.1031275
Received: 29 August 2022; Accepted: 15 December 2022;
Published: 12 January 2023.
Edited by:
Ayal Hendel, Bar-Ilan University, IsraelReviewed by:
Adam Ameur, Uppsala University, SwedenGavin Kurgan, Integrated DNA Technologies, United States
Roberto Nitsch, AstraZeneca, Sweden
Copyright © 2023 Schmidt, Kim, Strelchenko, Gierczic, Pavelec, Golos and Slukvin. 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: Jenna Kropp Schmidt, amtyb3BwQHdpc2MuZWR1
†These authors share senior authorship