Literature DB >> 34674375

Zona pellucida treatment before CRISPR/Cas9-mediated genome editing of porcine zygotes.

Zhao Namula1,2, Maki Hirata2, Quynh Anh Le2, Qingyi Lin2, Koki Takebayashi2, Naoaki Yoshimura2, Fuminori Tanihara2, Chommanart Thongkittidilok3, Takeshige Otoi2.   

Abstract

BACKGROUND: Increasing the permeability of the zona pellucida (ZP) of oocytes before CRISPR/Cas9 electroporation may improve the efficiency of gene editing; however, the effects of this approach on subsequent developmental processes are unclear. In this study, the effects of ZP treatment before electroporation on embryonic development and gene editing in porcine embryos were evaluated.
METHODS: The ZP of zygotes was weakened or removed by exposure to 0.5% actinase E, followed by electroporation of the Cas9 protein with guide RNA targeting GGTA1.
RESULTS: The blastocyst formation rate of ZP-free zygotes after electroporation was significantly lower (p < 0.05) than that of ZP-intact zygotes. The mutation rate in blastocysts from ZP-weakened zygotes was similar to that in ZP-intact zygotes, whereas ZP removal increased the mutation rate. The mutation efficiency in blastocysts from electroporated zygotes did not differ among ZP treatment groups.
CONCLUSIONS: Our results indicate that weakening the ZP does not affect the developmental competence, mutation rate, or mutation efficiency of electroporated zygotes, whereas ZP removal has a detrimental effect on embryonic development but may increase the mutation rate.
© 2021 The Authors. Veterinary Medicine and Science published by John Wiley & Sons Ltd.

Entities:  

Keywords:  CRISPR; electroporation; genome editing; porcine zygote; zona pellucida

Mesh:

Substances:

Year:  2021        PMID: 34674375      PMCID: PMC8788957          DOI: 10.1002/vms3.659

Source DB:  PubMed          Journal:  Vet Med Sci        ISSN: 2053-1095


INTRODUCTION

Pigs are considered one of the best animals for generating models of human diseases and for providing organs for xenotransplantation because their physiological and genetic backgrounds are similar to those of humans (Abkowitz et al., 1995; Meurens et al., 2012). Recently, genome editing technologies, such as the clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9) system, have allowed us to efficiently generate genetically engineered pigs in a short period of time (Tanihara et al., 2016; Wang et al., 2015). We established the gene editing by electroporation of Cas9 protein (GEEP) method for the highly efficient disruption of a targeted gene by electroporation (Tanihara et al., 2016). This is a simple application of the CRISPR/Cas9 system in which the Cas9 protein and guide RNA (gRNA) are introduced into in vitro‐fertilised zygotes by electroporation. The GEEP method does not require complex techniques associated with micromanipulation (such as microinjection) for the introduction of the CRISPR/Cas9 system into zygotes. Electroporation‐based protocols must confer high levels of gene editing while minimising damage to embryos. Therefore, optimizing electroporation conditions with respect to embryonic development and gene editing are critical. We have previously demonstrated that the mutation frequency in targeted genes and embryonic development are influenced by various electroporation parameters, such as the duration, number, polarity, and voltage of pulses (Nishio et al., 2018; Tanihara et al., 2016). Moreover, the concentration of Cas9 affects the mutation rate and gene editing efficiency in porcine embryos (Le et al., 2020). The efficiency of electroporation‐based protocols for a variety of simple editing schemes is well‐established; however, the application of this approach to more complex genome engineering, such as the insertion of fluorescent tags or conditional alleles, remains unclear. Therefore, complex editing schemes may require alternative approaches, such as the microinjection of macromolecules, and the zona pellucida (ZP) surrounding the plasma membrane of an oocyte may therefore present a barrier preventing the access of CRISPR/Cas9 components to the zygote by electroporation. In mice, the pre‐treatment of zygotes with Tyrode's acidic solution has been used to increase the permeability of the ZP for large RNA molecules during electroporation (Peng et al., 2012; Qin et al., 2015). However, the ZP is important for embryonic development, wherein the developmental potential of ZP‐free embryos is often lower than that of ZP‐intact embryos. Information on suitable conditions of the ZP for the introduction of the CRISPR/Cas9 system into porcine embryos by electroporation is limited. In the present study, we investigated whether ZP treatment before electroporation influences embryonic development and the targeted mutation efficiency in porcine embryos from zygotes that were edited via the GEEP method, using gRNA targeting the GGTA1 gene.

MATERIALS AND METHODS

No live animals were used in this study; accordingly, ethical approval was not required. All procedures were approved by the Animal Research Committee of Tokushima University.

OOCYTE COLLECTION, IN VITRO MATURATION, AND FERTILISATION

Oocyte collection, in vitro maturation, fertilisation, and embryo culture were carried out as described previously (Nishio et al., 2018). Pig ovaries were obtained from prepubertal, crossbreed gilts (Landrace × Large White × Duroc breeds) at a local slaughterhouse. Cumulus–oocyte complexes (COCs) were collected from 3 to 6 mm diameter follicles on the ovarian surface by the slicing method. Approximately 50 COCs were cultured in 500 μl of maturation medium consisting of tissue culture medium 199 (TCM 199) with Earle's salts (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 10% (v/v) porcine follicular fluid, 0.6 mM cysteine (Sigma‐Aldrich, St. Louis, MO, USA), 50 μM β‐mercaptoethanol (Wako Pure Chemical Industries Ltd., Osaka, Japan), 50 μM sodium pyruvate (Sigma‐Aldrich), 2 mg/ml d‐sorbitol (Wako Pure Chemical Industries Ltd.), 10 IU/ml equine chorionic gonadotropin (Kyoritu Seiyaku, Tokyo, Japan), 10 IU/ml human chorionic gonadotropin (Kyoritu Seiyaku), and 50 μg/ml gentamicin (Sigma‐Aldrich), covered with mineral oil (Sigma‐Aldrich) for 22 h in 4‐well dishes (Nunc A/S, Roskilde, Denmark). The COCs were transferred into a maturation medium without hormones and then cultured for an additional 22 h. The COCs were incubated at 39°C in a humidified incubator containing 5% CO2. For in vitro fertilisation (IVF), frozen‐thawed boar spermatozoa (Landrace breed) were transferred into 6 ml of porcine fertilisation medium (PFM; Research Institute for Functional Peptides Co., Yamagata, Japan) and washed by centrifugation at 500 × g for 5 min. The pelleted spermatozoa were resuspended in PFM and adjusted to 5 × 106 cells/ml. The matured oocytes were then transferred to the sperm‐containing PFM and co‐incubated in a humidified incubator containing 5% CO2, 5% O2, and 90% N2 for 5 h at 39°C. After co‐incubation, the putative zygotes were denuded from the cumulus cells and the attached spermatozoa by mechanical pipetting, transferred to porcine zygote medium (PZM‐5; Research Institute for the Functional Peptides Co.), and then cultured for 7 h until electroporation treatment. After 7 h of culture, zygotes were randomly assigned to each experimental group.

Zona pellucida treatment

Before electroporation, the ZP of zygotes (1‐cell stage) was weakened or removed by actinase E treatment (Figure 1a–c). Briefly, denuded zygotes were exposed to 0.5% (w/v) actinase E (Kaken‐Seiyaku Co., Tokyo, Japan) in Dulbecco's PBS (Thermo Fisher Scientific) for 3–15 s and then transferred to PZM‐5 to stop the activity of actinase E. Some zygotes were freed completely from their ZP by gentle pipetting. The ZP‐weakened or ZP‐free zygotes were incubated at 39°C in a humidified incubator containing 5% CO2, 5% O2, and 90% N2 for 1 h before electroporation.
FIGURE 1

Representative images of porcine zygotes with intact (a), weakened (b), and free (c) zona pellucida, and blastocysts derived from zygotes with intact (d), weakened (e), and free (f) zona pellucida. Blastocysts derived from zygotes with weakened zona pellucida lost the zona pellucida during embryonic development. Scale bar = 50 μm

Representative images of porcine zygotes with intact (a), weakened (b), and free (c) zona pellucida, and blastocysts derived from zygotes with intact (d), weakened (e), and free (f) zona pellucida. Blastocysts derived from zygotes with weakened zona pellucida lost the zona pellucida during embryonic development. Scale bar = 50 μm

Electroporation

Electroporation was performed as described previously (Tanihara et al., 2016). Briefly, the electrode (LF501PT1‐20; BEX Co. Ltd., Tokyo, Japan) was connected to a CUY21EDIT II electroporator (BEX Co. Ltd.) and placed under a stereoscopic microscope. The zygotes with or without ZP treatment were washed with Opti‐MEM I solution (Thermo Fisher Scientific) and placed in a line in the electrode gap on the chamber slide filled with 10 μl of nuclease‐free duplex buffer (Integrated DNA Technologies (IDT), Coralville, IA, USA) with 100 ng/μl gRNA (Alt‐R CRISPR crRNAs and tracrRNA, chemically modified and length‐optimised variants of the native guide RNA purchased from IDT) targeting GGTA1 and 100 ng/μl Cas9 protein (Guide‐it Recombinant Cas9; Takara Bio, Shiga, Japan). A gRNA targeting GGTA1 (5′‐GGATTAAACCAGTCCACTAG‐3′) resulting in a high biallelic mutation rate in our previous study was used for electroporation treatment (Tanihara et al., 2020). After electroporation with five 1‐msec pulses at 25 V, the zygotes with or without ZP treatment were washed with PZM‐5 and cultured in 4‐well dishes (∼35 embryos per well) for 3 days. The embryos were subsequently incubated in porcine blastocyst medium (PBM; Research Institute for the Functional Peptides Co.) for 4 days for evaluations of blastocyst formation and for blastocyst genotyping (Figure 1d–f). As a control, some zygotes were cultured in PZM‐5 and PBM for 7 days without electroporation. Zygotes and embryos were incubated at 39°C in a humidified incubator containing 5% CO2, 5% O2, and 90% N2.

Analysis of targeted gene sequences in embryos

Targeted genes in individual blastocysts were analysed as described previously (Tanihara et al., 2020). Briefly, genomic DNA was isolated from individual blastocysts by boiling in a 50 mM NaOH solution at 98°C for 10 min and then neutralised with 100 mM Tris HCl. After neutralisation, the genomic regions flanking the gRNA target sequences were amplified by polymerase chain reaction (PCR) using the following primer pair: GGTA1, 5′‐AAAAGGGGAGCACTGAACCT‐3′ (forward) and 5′‐ATCCGGACCCTGTTTTAAGG‐3′ (reverse). The PCR products were separated by agarose gel electrophoresis and purified using a Fast Gene Gel/PCR Extraction Kit (Nippon Genetics, Tokyo, Japan). The PCR products were directly sequenced by Sanger sequencing using the BigDye Terminator Cycle Sequencing Kit version 3.1 (Thermo Fisher Scientific K.K., Tokyo, Japan) and an ABI 3500 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). The tracking of indels by decomposition (TIDE) bioinformatics package was used to determine the genotype of each blastocyst (http://shinyapps.datacurators.nl/tide/) (Brinkman et al., 2018). Blastocysts that did not carry wild‐type (WT) sequences were classified as having biallelic mutations. Blastocysts that carried more than one type of mutation in addition to the WT sequence were classified as mosaics. Blastocysts that carried only the WT sequence were classified as WT (Figure 2). The mutation rate was defined as the ratio of the number of blastocysts carrying a mutation to the total number of sequenced blastocysts. The mutation efficiency was defined as the proportion of indel mutation events in blastocysts carrying biallelic or mosaic mutations.
FIGURE 2

Representative results of Sanger sequencing of blastocysts formed after electroporation with Cas9 protein and gRNAs. Dotted lines indicate the predicted Cas9 cleavage sites. Blue nucleotides indicate the target sequences, and red nucleotides indicate the protospacer adjacent motif (PAM) sequences. Biallelic, bi‐allelic mutation; mosaic: mosaic mutation; WT, wild‐type

Representative results of Sanger sequencing of blastocysts formed after electroporation with Cas9 protein and gRNAs. Dotted lines indicate the predicted Cas9 cleavage sites. Blue nucleotides indicate the target sequences, and red nucleotides indicate the protospacer adjacent motif (PAM) sequences. Biallelic, bi‐allelic mutation; mosaic: mosaic mutation; WT, wild‐type

Statistical analysis

Data for embryonic development were evaluated by an analysis of variance (ANOVA) followed by Fisher's protected least significant difference tests using STATVIEW (Abacus Concepts, Inc., Berkeley, CA, USA). All percentages were subjected to arcsine transformation before ANOVA. The percentages of mutations within all sequenced blastocysts were analysed by a chi‐squared test with Yates’ correction. Differences with a p‐value of ≤0.05 were regarded as significant.

RESULTS

The blastocyst formation rates of the electroporated zygotes were significantly lower (p < 0.05) than those of the control zygotes, irrespective of ZP treatment (Table 1). Moreover, the blastocyst formation rate of the ZP‐free zygotes after electroporation was significantly lower (p < 0.05) than that of ZP‐intact zygotes. In contrast, the mutation rate obtained by sequencing target sites in GGTA1 in blastocysts from ZP‐free zygotes was significantly higher (p < 0.05) than that in blastocysts from ZP‐intact zygotes. The biallelic mutation rates in the examined blastocysts and the mutation efficiency in the mutant blastocysts did not differ among the ZP treatment groups.
TABLE 1

Effects of zona pellucida (ZP) treatment before electroporation with the clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9) system on the development of zygotes and the mutation rate in blastocysts

No. (mean ± SEM) of embryosNo. (mean) of blastocysts
Zygote †† No. of zygotes examinedCleavedDeveloped to blastocystsNo. of blastocysts examinedMutant Biallelic mutant § Total mutation efficiency (mean ± SEM)
Control199170 (85.4 ± 2.3)34 (17.1 ± 1.9) a
ZP‐intact208171 (82.2 ± 3.0)23 (11.1 ± 2.4) b 207 (35.0) a 2 (10.0)69.6 ± 8.4
ZP‐weakened199163 (81.4 ± 5.0)16 (8.2 ± 1.7) b,c 127 (58.3) a,b 2 (16.7)71.9 ± 8.0
ZP‐free203175 (85.9 ± 4.1)6 (3.0 ± 0.9) c 66 (100) b 3 (50.0)88.9 ± 6.7

Four replicate trials were carried out.

Electroporation was performed by five 1‐msec pulses at 25 V. As a control, zygotes with an intact ZP were cultured without electroporation.

Proportions were calculated by dividing the number of mutant blastocysts by the total number of sequenced blastocysts.

Proportions were calculated by dividing the number of biallelic mutant blastocysts by the total number of sequenced blastocysts.

Mean proportions represent the proportion of indel mutation events in mutant blastocysts determined by the TIDE bioinformatics package.

a‐cValues with different superscript letters in the same column are significantly different (p < 0.05).

Effects of zona pellucida (ZP) treatment before electroporation with the clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9) system on the development of zygotes and the mutation rate in blastocysts Four replicate trials were carried out. Electroporation was performed by five 1‐msec pulses at 25 V. As a control, zygotes with an intact ZP were cultured without electroporation. Proportions were calculated by dividing the number of mutant blastocysts by the total number of sequenced blastocysts. Proportions were calculated by dividing the number of biallelic mutant blastocysts by the total number of sequenced blastocysts. Mean proportions represent the proportion of indel mutation events in mutant blastocysts determined by the TIDE bioinformatics package. a‐cValues with different superscript letters in the same column are significantly different (p < 0.05).

DISCUSSION

Electroporation of zygotes is an effective alternative to microinjection for CRISPR/Cas9‐mediated genome editing in pigs. Electroporation is a rapid procedure and can be applied to many zygotes simultaneously, enabling high‐throughput production. However, the successful CRISPR‐mediated generation of large genomic integrations by electroporation into zygotes has been limited by the inability of large molecules to reach the nucleus (Remy et al., 2017). Insufficient permeability of the nuclear membrane may also be responsible for low transgene integration upon electroporation. Moreover, in 1‐cell‐stage embryos, the ZP makes the introduction of macromolecules difficult. To overcome this physical barrier, the weakening of the ZP before electroporation has been performed to increase permeability for efficient introduction (Peng et al., 2012; Qin et al., 2015). To date, enzymatic digestion with pronase for ZP removal has been applied to oocytes and embryos in domestic and laboratory species, including cattle, pigs, and sheep (Li et al., 2013; Peura & Vajta, 2003). It has been suggested that ZP removal by pronase improves the quality of porcine parthenogenetic embryos by accelerating the speed of embryonic development and decreasing the number of apoptotic cells in blastocysts (Li et al., 2013). In the present study, therefore, we used actinase E (a synonym of pronase) for the ZP treatment. We found that ZP removal had a detrimental effect on the blastocyst formation rates of electroporated zygotes; however, the blastocyst formation rate was similar in ZP‐weakened and ZP‐intact zygotes. When the quality of blastocysts from the electroporated zygotes with and without the ZP treatment was examined (number of examined embryos in each group, 9–10), there were no differences in the total cell numbers of embryos among the groups (83.4–121.8 cells). However, the apoptotic rate of blastocysts from ZP‐weakened zygotes (7.0%) was significantly higher (p < 0.05) than that of blastocysts from ZP‐free zygotes (2.9%) but was similar to that from ZP‐free zygotes (6.5%) (data not shown). These results indicate that the ZP treatment may affect embryonic development or quality. It has been reported that, in mice, the weakening of the ZP by acidic Tyrode's solution does not compromise the survival of the electroporated embryos and improves access to embryos by CRISPR/Cas9 reagents, resulting in a high targeting efficiency (Peng et al., 2012; Qin et al., 2015). Accordingly, under identical electric pulses, reagents can enter the weakened ZP more easily. However, our results showed that the mutation rate in blastocysts from ZP‐weakened zygotes was similar to that in blastocysts from ZP‐intact zygotes, whereas the removal of the ZP increased the mutation rate in blastocysts. There were no differences in the biallelic mutation rates in blastocysts and the mutation efficiency in mutant blastocysts among the electroporation groups; however, half of the embryos from ZP‐free zygotes had biallelic mutations. Previous studies have suggested that the weakening of the ZP before electroporation is necessary for gene editing using Cas9 mRNA of about 4,500 nucleotides length but does not substantially improve the uptake of small molecules, like morpholinos of only 25 nucleotides in length (Peng et al., 2012; Qin et al., 2015). These studies are consistent with our results for mutation rates and efficiencies, demonstrating that the weakening of the ZP does not increase the mutation rate in blastocysts from zygotes electroporated with the Cas9 protein (160 kDa), which is much smaller than Cas9 mRNA (∼1500 kDa). Although electroporation variation, such as a higher voltage, may enhance the mutation efficiency in blastocysts from ZP‐weakened zygotes, our results indicated that ZP removal increased the mutation rate in embryos and biallelic mutation rates. In conclusion, the weakening of the ZP did not affect the developmental competence of electroporated zygotes and did not improve the mutation rate or mutation efficiency in blastocysts. In contrast, ZP removal had a detrimental effect on embryonic development but increased the mutation rate of embryos and biallelic mutations.

CONFLICT OF INTEREST

The authors declare no conflict of interest.

ETHICS STATEMENT

The authors confirm that the ethical policies of the journal, as noted on the journal's author guidelines page, have been adhered to and the appropriate ethical review committee approval has been received. All the animals involved in this study received humane care in compliance with the Guide for the Care and Use of Laboratory Animals prepared by the Institute of Laboratory Animal Resources, National Research Council. All procedures were approved by the Animal Research Committee of Tokushima University.

AUTHOR CONTRIBUTIONS

Writing – original draft preparation and investigation: Zhao Namula. Conceptualisation and writing – review & editing: Maki Hirata. Investigation: Quynh Anh Le, Qingyi Lin, and Koki Takebayashi. Conceptualisation and writing – review & editing: Naoaki Yoshimura. Investigation, formal analysis, and writing – review & editing: Fuminori Tanihara. Methodology and writing – review & editing: Chommanart Thongkittidilok. Conceptualisation, project administration, and writing – review & editing: Takeshige Otoi.

PEER REVIEW

The peer review history for this article is available at https://publons.com/publon/10.1002/vms3.659
  14 in total

1.  Effects of voltage strength during electroporation on the development and quality of in vitro-produced porcine embryos.

Authors:  K Nishio; F Tanihara; T-V Nguyen; T Kunihara; M Nii; M Hirata; T Takemoto; T Otoi
Journal:  Reprod Domest Anim       Date:  2017-11-14       Impact factor: 2.005

2.  Behavior of hematopoietic stem cells in a large animal.

Authors:  J L Abkowitz; M T Persik; G H Shelton; R L Ott; J V Kiklevich; S N Catlin; P Guttorp
Journal:  Proc Natl Acad Sci U S A       Date:  1995-03-14       Impact factor: 11.205

3.  Development and quality of porcine parthenogenetically activated embryos after removal of zona pellucida.

Authors:  Rong Li; Ying Liu; Hanne Skovsgaard Pedersen; Peter M Kragh; Henrik Callesen
Journal:  Theriogenology       Date:  2013-04-17       Impact factor: 2.740

4.  Effects of electroporation treatment using different concentrations of Cas9 protein with gRNA targeting Myostatin (MSTN) genes on the development and gene editing of porcine zygotes.

Authors:  Quynh A Le; Maki Hirata; Nhien T Nguyen; Koki Takebayashi; Manita Wittayarat; Yoko Sato; Zhao Namula; Masahiro Nii; Fuminori Tanihara; Takeshige Otoi
Journal:  Anim Sci J       Date:  2020-01       Impact factor: 1.749

5.  Efficient CRISPR/Cas9-Mediated Genome Editing in Mice by Zygote Electroporation of Nuclease.

Authors:  Wenning Qin; Stephanie L Dion; Peter M Kutny; Yingfan Zhang; Albert W Cheng; Nathaniel L Jillette; Ankit Malhotra; Aron M Geurts; Yi-Guang Chen; Haoyi Wang
Journal:  Genetics       Date:  2015-03-27       Impact factor: 4.562

6.  Generation of gene-edited rats by delivery of CRISPR/Cas9 protein and donor DNA into intact zygotes using electroporation.

Authors:  Séverine Remy; Vanessa Chenouard; Laurent Tesson; Claire Usal; Séverine Ménoret; Lucas Brusselle; Jean-Marie Heslan; Tuan Huan Nguyen; Jeremy Bellien; Jean Merot; Anne De Cian; Carine Giovannangeli; Jean-Paul Concordet; Ignacio Anegon
Journal:  Sci Rep       Date:  2017-11-29       Impact factor: 4.379

7.  Zona pellucida treatment before CRISPR/Cas9-mediated genome editing of porcine zygotes.

Authors:  Zhao Namula; Maki Hirata; Quynh Anh Le; Qingyi Lin; Koki Takebayashi; Naoaki Yoshimura; Fuminori Tanihara; Chommanart Thongkittidilok; Takeshige Otoi
Journal:  Vet Med Sci       Date:  2021-10-21

8.  Efficient delivery of DNA and morpholinos into mouse preimplantation embryos by electroporation.

Authors:  Hui Peng; Yongyan Wu; Yong Zhang
Journal:  PLoS One       Date:  2012-08-21       Impact factor: 3.240

9.  Somatic cell reprogramming-free generation of genetically modified pigs.

Authors:  Fuminori Tanihara; Tatsuya Takemoto; Eri Kitagawa; Shengbin Rao; Lanh Thi Kim Do; Akira Onishi; Yukiko Yamashita; Chisato Kosugi; Hitomi Suzuki; Shoichiro Sembon; Shunichi Suzuki; Michiko Nakai; Masakazu Hashimoto; Akihiro Yasue; Munehide Matsuhisa; Sumihare Noji; Tatsuya Fujimura; Dai-Ichiro Fuchimoto; Takeshige Otoi
Journal:  Sci Adv       Date:  2016-09-14       Impact factor: 14.136

10.  Efficient generation of GGTA1-deficient pigs by electroporation of the CRISPR/Cas9 system into in vitro-fertilized zygotes.

Authors:  Fuminori Tanihara; Maki Hirata; Nhien Thi Nguyen; Osamu Sawamoto; Takeshi Kikuchi; Masako Doi; Takeshige Otoi
Journal:  BMC Biotechnol       Date:  2020-08-18       Impact factor: 2.563

View more
  1 in total

1.  Zona pellucida treatment before CRISPR/Cas9-mediated genome editing of porcine zygotes.

Authors:  Zhao Namula; Maki Hirata; Quynh Anh Le; Qingyi Lin; Koki Takebayashi; Naoaki Yoshimura; Fuminori Tanihara; Chommanart Thongkittidilok; Takeshige Otoi
Journal:  Vet Med Sci       Date:  2021-10-21
  1 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.