CRUCIAN CARP (CARASSIUS AURATUS) STRAIN WITHOUT INTERMUSCULAR BONES AND BREEDING METHOD THEREOF

Abstract

The disclosure provides a method for breeding Carassius auratus strains without intermuscular bones, which comprises designing knockout target sites for the two copies of the bmp6, bmp6a and bmp6b, in the Carassius auratus genome, and then obtaining F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation through two rounds of gene knockout and screening, and then propagating by using the F2 generation individuals from the homozygous line with bmp6a and bmp6b double gene mutation to form a new Carassius auratus strain with intermuscular bone-deficient. Strains of Carassius auratus with less than 20 intermuscular bones and without intermuscular bones are obtained by the method of the present disclosure.

Claims

1. A method for breeding Carassius auratus strains without intermuscular bones, wherein the breeding of the Carassius auratus strains with intermuscular bone-deficient according to following steps: designing knockout target sites for two copies of bmp6 gene in the Carassius auratusn genome, bmp6a and bmp6b, respectively, and then obtaining F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation through two rounds of gene knockout and screening, and then using the F2 generation individuals from the homozygous line with bmp6a and bmp6b double gene mutation to propagate to form a new Carassius auratus strain without intermuscular bones; wherein, mixing a sgRNA corresponding to (bmp6a), (bmp6b) or (bmp6a and bmp6b) with Cas9 protein and microinjecting a mixture obtained into Carassius auratus embryos in the single cell stage, to perform a first round of gene knockout to construct a F0 generation population, culturing the F0 generation population for 3 to 5 months followed by PIT labeling and DNA extraction, then sequencing to determine alleles and mutation rates of somatic mutations of Carassius auratus, selecting F0 generation individuals with a somatic mutation rate of more than 95% as parents to prepare 0-generation fertilized eggs ; mixing the sgRNAs corresponding to (bmp6a), (bmp6b) or (bmp6a and bmp6b) with Cas9 protein and microinjecting a mixture obtained into the 0-generation fertilized eggs to perform a second round of gene knockout to construct a F1 generation population, culturing the F1 generation population for 3 to 5 months followed by PIT labeling and DNA extraction, then sequencing to determine alleles and mutation rates of somatic mutations of Carassius auratus, selecting F1 generation individuals from a somatic bmp6a and bmp6b double gene mutant line with a somatic mutation rate of more than 95% as parents for reproduction to construct F2 generation, then selecting F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation.

2. The method for breeding Carassius auratus strains without intermuscular bones according to claim 1, wherein the knockout target sites of bmp6 gene are shown in Table 1, TABLE-US-00003 TABLE1 sgRNAs Genome Gene Number coordinates Exon Targetsequence bmp6a CAA1 NC_039266.1: Exon1 GGTAGCAGTTCTG 5784476 CACTAGACGG (SEQIDNO:1) CAA2 NC_039266.1: Exon1 CCAGCGGAGGTTG 5784523 CGGACTCAGG (SEQIDNO:2) CAA3 NC_039266.1: Exon1 GAAAGAGATTCTG 5784562 TCCATACTGG (SEQIDNO:3) CAA4 NC_039266.1: Exon2 TCTTTATGGTGTCT 5815487 TCTCTGTGG (SEQIDNO:4) CAA5 NC_039266.1: Exon2 GGCCTCTCCCTCT 5815527 GGAATCTGGG (SEQIDNO:5) CAA6 NC_039266.1: Exon2 TGCAGAATTCAGG 5815558 ATCTACAAGG (SEQIDNO:6) bmp6b CAA7 NC_039291.1: Exon1 GGCAGGTAGTAGT 8127537 TCTGTACTGG (SEQIDNO:7) CAA8 NC_039291.1: Exon1 AGCCCAACTTCAT 8127477 TCATCGGAGG (SEQIDNO:8) CAA9 NC_03929.1: Exon1 GAAAGAGATCCTG 8127425 TCCATACTGG (SEQIDNO:9) CAA10 NC_03929.1: Exon2 GGCCTCTCCCTCT 8097264 GGAATCTGGG (SEQIDNO:10) CAA11 NC_039291.1: Exon2 AGCAGAATTCAGG 8097233 ATCTACAAGG (SEQIDNO:11)

3. The method for breeding Carassius auratus strains without intermuscular bones according to claim 2, wherein sgRNA forward primers designed for the knockout target sites of the bmp6 gene are shown in the table below, TABLE-US-00004 sgRNAs Gene Number sgRNAforwardprimer bmp6a CAA1 TTCTAATACGACTCACTATAGGTAGCAGTT CTGCACTAGAGTTTTAGAGCTAGA (SEQIDNO:12) CAA2 TTCTAATACGACTCACTATAGCCAGCGGAG GTTGCGGACTCGTTTTAGAGCTAGA (SEQIDNO:13) CAA3 TTCTAATACGACTCACTATAGGAAAGAGAT TCTGTCCATACGTTTTAGAGCTAGA (SEQIDNO:14) CAA4 TTCTAATACGACTCACTATAGTCTTTATGG TGTCTTCTCTGGTTTTAGAGCTAGA (SEQIDNO:15) CAA5 TTCTAATACGACTCACTATAGGCCTCTCCC TCTGGAATCTGTTTTAGAGCTAGA (SEQIDNO:16) CAA6 TTCTAATACGACTCACTATAGTGCAGAATT CAGGATCTACAGTTTTAGAGCTAGA (SEQIDNO:17) bmp6b CAA7 TTCTAATACGACTCACTATAGGCAGGTAGT AGTTCTGTACGTTTTAGAGCTAGA (SEQIDNO:18) CAA8 TTCTAATACGACTCACTATAGAGCCCAACT TCATTCATCGGGTTTTAGAGCTAGA (SEQIDNO:19) CAA9 TTCTAATACGACTCACTATAGGAAAGAGAT CCTGTCCATACGTTTTAGAGCTAGA (SEQIDNO:20) CAA10 TTCTAATACGACTCACTATAGGCCTCTCCC TCTGGAATCTGTTTTAGAGCTAGA (SEQIDNO:21) CAA11 TTCTAATACGACTCACTATAGAGCAGAATT CAGGATCTACAGTTTTAGAGCTAGA (SEQIDNO:22)

4. The method for breeding Carassius auratus strains without intermuscular bones according to claim 3, wherein using a sgRNA forward primer and a sgRNA reverse primer with the sequence of 5-GATCCGCACCGACTCGGTGCCACTTTTTCAAGTTGATAACGGACTAGCCTTATTTTA ACTTGCTATTTCTAGCTCTAAAAC-3 (SEQ ID NO:31) to synthesis sgRNA in vitro; a sgRNA synthesis and amplification volume is: 2.5 L of 10 M sgRNA forward and reverse primers, 25 L of 2Dream Taq Master mix, and supplementing the sgRNA synthesis and amplification volume to 50 L with enzyme-free water; sgRNA synthesis procedure in vitro is: denaturation at 95 C. for 3 min; then 30 cycles of 95 C. for 30 s, 58 C. for 30 s and 72 C. for 30 s, extension at 72 C. for 5 min; after synthesizing the sgRNA in vitro, purifying and recovering an obtained sgRNA PCR product, then using an RNA in vitro transcription kit to transcribe the sgRNA in vitro; establishing a 30 L reaction volume for sgRNA in vitro transcription for each target site: 1 g of sgRNA PCR recovered product, 10 L of NTP Buffer Mix, and 2 L of T7 RNA Polymerase Mix, supplementing the sgRNA in vitro transcription reaction volume to 30 L with enzyme-free water; transcribing at 37 C. for 4 h, adding 20 L of enzyme-free water after the reaction, mixing well, adding 2 L DNase I, digesting at 37 C. for 15 min to remove DNA.

5. The method for breeding Carassius auratus strains without intermuscular bones according to claim 4, wherein a method for the microinjection is: mixing the sgRNA synthesized in vitro and Cas9 protein in a molar concentration ratio of 3:1, then incubating at room temperature for 10 minutes, adding 25% phenol red and injecting the obtained product into Carassius auratus embryos in a single cell stage; wherein, mixing target site sgRNAs on each exon in equal amounts before injection, and a final concentration of each sgRNA is no less than 50 ng/L, an injection volume of each fertilized egg is 1 nL0.02 nL.

6. The method for breeding Carassius auratus strains without intermuscular bones according to claim 4, wherein a recovery concentration of the sgRNA PCR product is in a range of 100 to 160 ng/L.

7. The method for breeding Carassius auratus strains without intermuscular bones according to claim 4, wherein a recovery concentration of the transcribed sgRNA in vitro is a range of 800 to 3000 ng/L.

8. A Carassius auratus strain without intermuscular bones obtained by using the method for breeding Carassius auratus strains without intermuscular bones in any one of claims 1.

9. A reagent for knocking out the bmp6 gene in Carassius auratus genome, the reagent comprising sgRNA forward primers and sgRNA reverse primers designed for knockout target sites of bmp6 gene, the sgRNA forward primers and sgRNA reverse primers are shown in Table 1.

10. The reagent according to claim 9, wherein the sgRNA forward primers are shown in following table, TABLE-US-00005 sgRNAs Gene Number sgRNAforwardprimer bmp6a CAA1 TTCTAATACGACTCACTATAGGTAGCAGTTCT GCACTAGAGTTTTAGAGCTAGA (SEQIDNO:12) CAA2 TTCTAATACGACTCACTATAGCCAGCGGAGGT TGCGGACTCGTTTTAGAGCTAGA (SEQIDNO:13) CAA3 TTCTAATACGACTCACTATAGGAAAGAGATTC TGTCCATACGTTTTAGAGCTAGA (SEQIDNO:14) CAA4 TTCTAATACGACTCACTATAGTCTTTATGGTG TCTTCTCTGGTTTTAGAGCTAGA (SEQIDNO:15) CAA5 TTCTAATACGACTCACTATAGGCCTCTCCCTC TGGAATCTGTTTTAGAGCTAGA (SEQIDNO:16) CAA6 TTCTAATACGACTCACTATAGTGCAGAATTCA GGATCTACAGTTTTAGAGCTAGA (SEQIDNO:17) bmp6b CAA7 TTCTAATACGACTCACTATAGGCAGGTAGTAG TTCTGTACGTTTTAGAGCTAGA (SEQIDNO:18) CAA8 TTCTAATACGACTCACTATAGAGCCCAACTTC ATTCATCGGGTTTTAGAGCTAGA (SEQIDNO:19) CAA9 TTCTAATACGACTCACTATAGGAAAGAGATCC TGTCCATACGTTTTAGAGCTAGA (SEQIDNO:20) CAA10 TTCTAATACGACTCACTATAGGCCTCTCCCTC TGGAATCTGTTTTAGAGCTAGA (SEQIDNO:21) CAA11 TTCTAATACGACTCACTATAGAGCAGAATTCA GGATCTACAGTTTTAGAGCTAGA (SEQIDNO:22)

11. The reagent according to claim 9, wherein the nucleotide sequence of the sgRNA reverse primer is: 5-GATCCGCACCGACTCGGTGCCACTTTTTCAAGTTGATAACGGACTAGCCTTATTTTA ACTTGCTATTTCTAGCTCTAAAAC-3 (SEQ ID NO:31).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0029] FIG. 1 is a comparison picture of bone staining observation of a typical individual in Example 1, and arrows in the figure refer to intermuscular bones;

[0030] FIG. 2 is an X-ray diagram of wild-type Carassius auratus bone in Example 1, and the arrow in the figure indicates intermuscular bones;

[0031] FIG. 3 is an X-ray diagram of caudal bones of a mutant in Example 1, from which it can be seen that no intermuscular bone is present in the caudal muscle tissue;

[0032] FIG. 4 is an X-ray diagram of trunk bones of a mutant in Example 1, from which it can be seen that no intermuscular bone is present in back muscle tissue of trunk;

[0033] FIGS. 5A, 5B and 5C are diagrams showing the sequencing results of exon 1 of the new Carassius auratus strain bmp6a in Example 1; CAA-1, CAA-2, and CAA-3 are target sites of exon 1 of bmp6a; WT represents wild type, the number in sequence number is the PIT labeling number of mutant, and the letter followed by PIT labeling number is the label of mutant allele;

[0034] FIGS. 6A, 6B and 6C are diagrams showing the sequencing results of exon 1 of the new Carassius auratus strain bmp6b in Example 1; CAA-7, CAA-8, and CAA-9 are target sites of exon 1 of bmp6b; WT represents wild type, the number in sequence number is the PIT labeling number of mutants, and the letter followed by PIT labeling number is the label of mutant allele;

[0035] FIG. 7 is a comparison diagram of exon 1 protein sequences of the new Carassius auratus strain bmp6a in Example 1; WT represents wild type, the number in sequence number is the PIT labeling number of mutants, and the letter followed by PIT labeling number is the label of mutant allele;

[0036] FIG. 8 is a comparison diagram of exon 1 protein sequences of the new Carassius auratus strain bmp6b in Example 1; WT represents wild type, the number in sequence number is the PIT labeling number of mutants, and the letter followed by PIT labeling number is the label of mutant allele;

[0037] FIG. 9 is a micro-CT image of bones of wild-type and the new Carassius auratus strain in Example 1;

[0038] FIG. 10 is a viviperception diagram by X-ray of wild type, typical individuals with few intermuscular bones and typical individuals without intermuscular bones in Example 1, arrows indicate intermuscular bones.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0039] The present disclosure is further described in detail below in conjunction with examples. The following examples are intended to illustrate the present disclosure, but not to limit the scope herein.

[0040] The experimental methods used in the following examples are conventional methods unless otherwise specified. The materials, reagents, methods and instruments used, unless otherwise specified, are conventional materials, reagents, methods and instruments in the art, which are commercially available to those skilled in the art.

[0041] Embodiment 1: The breeding method of Carassius auratus strain without intermuscular bones of the present embodiment:

[0042] The target sites are designed to knock out two copies of bmp6 gene in the Carassius auratus genome, bmp6a and bmp6b, respectively, and then F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation are obtained through two rounds of gene knockout and screening, and then the F2 generation individuals from the homozygous line with bmp6a and bmp6b double gene mutation are propagated to form a new Carassius auratus strain without intermuscular bone;

[0043] Among them, the sgRNAs corresponding to (bmp6a), (bmp6b) or (bmp6a and bmp6b) are mixed with Cas9 protein, and the obtained mixture is microinjected into Carassius auratus embryos in the single cell stage, the first round of gene knockout is performed to construct the F0 generation population, and then the F0 generation population is cultured for 3 to 5 months followed by PIT labeling and DNA extraction then the alleles and mutation rates of somatic mutations of Carassius auratus are determined by sequencing, and the F0 generation individuals with the somatic mutation rate of more than 95% are selected as parents to prepare 0-generation fertilized eggs;

[0044] The sgRNAs corresponding to (bmp6a), (bmp6b) or (bmp6a and bmp6b) are mixed with Cas9 protein, and the obtained mixture is microinjected into 0-generation fertilized eggs, the second round of gene knockout is performed to construct the F1 generation population, and then the F1 generation population is cultured for 3 to 5 months followed by PIT labeling and DNA extraction, then the alleles and mutation rates of the somatic mutations of Carassius auratus are determined by sequencing, and the F1 generation individuals from a somatic bmp6a and bmp6b double gene mutant line with the somatic mutation rate of more than 95% are selected as parents for reproduction to construct an F2 generation, then F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation are selected.

[0045] Among them, the DNA extraction and sequencing method is as follows: 100 L of lysis solution is added (the lysis solution comprises 0.5 mg/mL proteinase K, 10 mM Tris (pH 8.0), 50 mM KCl, 0.3% Tween 20, 0.3% NP40) to 0.3-0.5cm.sup.2 fin clip samples; the lysis methods is: the PCR instrument is set at 55 C. for 6 h, 98 C. for 10 min for lysis; after lysis, the obtained product is mixed well, then centrifuged at 1000-2000 rpm for 2 min, and the supernatant is taken as a PCR amplification template.

[0046] PCR amplification: 2 pairs of primers are used for PCR amplification, of which, 1 pair of primers is the target site primers (the primer sequences used to amplify the target fragment are shown in Table 2 and concatenated by the M13 universal primer sequence at 5 end), and the other pair of primers is index primers (the M13 universal primer sequence and the index base sequence at the 5 end used to identify the sample). A two-step PCR amplification is performed. The PCR amplification volume in the first step is 10 L, including 1 L of genomic DNA supernatant, 0.5 L of 1 M forward and reverse target primers, and 5 L of 2Dream Taq Master Mix (Thermo Fisher, CA, USA), the volume is supplemented to 10 L with enzyme-free water; PCR program: 95 C. for 3 min; 10 cycles of 95 C. for 30 s, 60 C. for 30 s and 72 C. for 30 s; 72 C. for 2 min; PCR amplification in step 2 is that 0.5 L of 5 M forward and reverse index primers are added to the PCR amplification product obtained in step 1, the PCR program is set to 95 C. for 2 min; 6 cycles of 95 C. for 30 s, 58 C. for 30 s and 72 C. for 30 s; 15 cycles of 95 C. for 30 s, 72 C. for 30 s; 72 C. for 2 min.

[0047] After PCR amplification, the PCR products are mixed in equal amounts, and a DNA sequencing library is constructed with the TrueSeq Sample Preparation Kit, and is sequenced using Illumina MiSeq sequencing platform with 300 bp Pair-End mode. The data obtained by high-throughput sequencing are first used to identify amplification loci and samples using the method described by Tong et al..sup.[14]. The sequenced data of each sample in each PCR product is analyzed using the CRISPResso2.sup.[15] program to determine the alleles and mutation rates of somatic mutations.

[0048] Among them, Sanger sequencing is used for mutant detection: PCR amplification is performed on the samples using the primers shown in Table 2, and the volume is set to 25 L: 2 L of genomic DNA supernatant, 1 L of 10 M forward and reverse target primers, 12.5 L of 2Dream Taq Master Mix (Thermo Fisher, CA, USA), the volume is supplemented to 25 L with enzyme-free water; PCR program is 95 C. for 3 min; 35 cycles of 95 C. for 30 s, 60 C. for 30 s, 72 C. for 30 s; 72 C. for 5 min. PCR products are detected by 1.5% agarose gel electrophoresis. After detection, the products are purified and recovered, TA clone is performed, 2 L of PCR products are taken after colony PCR amplification, and detected with 8% polypropylene, and the colonies with different band sizes from the control are selected for Sanger sequencing.

[0049] Embodiment 2: The difference between this embodiment and Embodiment 1 is that the knockout target sites of the bmp6 gene are shown in Table 1 in this embodiment.

TABLE-US-00001 Knockouttargetsitesofbmp6gene sgRNAs Genome Gene Number coordinates Exon Targetsitessequence sgRNAforwardprimer bmp6a CAA1 NC_039266.1: Exon1 GGTAGCAGTTCTGCACTAGACGG TTCTAATACGACTCACTATAGGTAGCAG 5784476 (SEQIDNO:1) TTCTGCACTAGAGTTTTAGAGCTAGA (SEQIDNO:12) CAA2 NC_039266.1: Exon1 CCAGCGGAGGTTGCGGACTCAGG TTCTAATACGACTCACTATAGCCAGCGG 5784523 (SEQIDNO:2) AGGTTGCGGACTCGTTTTAGAGCTAGA (SEQIDNO:13) CAA3 NC_039266.1: Exon1 GAAAGAGATTCTGTCCATACTGG TTCTAATACGACTCACTATAGGAAAGAG 5784562 (SEQIDNO:3) ATTCTGTCCATACGTTTTAGAGCTAGA (SEQIDNO:14) CAA4 NC_039266.1: Exon2 TCTTTATGGTGTCTTCTCTGTGG TTCTAATACGACTCACTATAGTCTTTAT 5815487 (SEQIDNO:4) GGTGTCTTCTCTGGTTTTAGAGCTAGA (SEQIDNO:15) CAA5 NC_039266.1: Exon2 GGCCTCTCCCTCTGGAATCTGGG TTCTAATACGACTCACTATAGGCCTCTC 5815527 (SEQIDNO:5) CCTCTGGAATCTGTTTTAGAGCTAGA (SEQIDNO:16) CAA6 NC_039266.1: Exon2 TGCAGAATTCAGGATCTACAAGG TTCTAATACGACTCACTATAGTGCAGAA 5815558 (SEQIDNO:6) TTCAGGATCTACAGTTTTAGAGCTAGA (SEQIDNO:17) bmp6b CAA7 NC_039291.1: Exon1 GGCAGGTAGTAGTTCTGTACTGG TTCTAATACGACTCACTATAGGCAGGTA 8127537 (SEQIDNO:7) GTAGTTCTGTACGTTTTAGAGCTAGA (SEQIDNO:18) CAA8 NC_039291.1: Exon1 AGCCCAACTTCATTCATCGGAGG TTCTAATACGACTCACTATAGAGCCCAA 8127477 (SEQIDNO:8) CTTCATTCATCGGGTTTTAGAGCTAGA (SEQIDNO:19) CAA9 NC_03929.1: Exon1 GAAAGAGATCCTGTCCATACTGG TTCTAATACGACTCACTATAGGAAAGAG 8127425 (SEQIDNO:9) ATCCTGTCCATACGTTTTAGAGCTAGA (SEQIDNO:20) CAA10 NC_03929.1: Exon2 GGCCTCTCCCTCTGGAATCTGGG TTCTAATACGACTCACTATAGGCCTCTC 8097264 (SEQIDNO:10) CCTCTGGAATCTGTTTTAGAGCTAGA (SEQIDNO:21) CAA11 NC_039291.1: Exon2 AGCAGAATTCAGGATCTACAAGG TTCTAATACGACTCACTATAGAGCAGAA (SEQIDNO:11) TTCAGGATCTACAGTTTTAGAGCTAGA (SEQIDNO:22)

[0050] Others are the same as Embodiment 1.

TABLE-US-00002 PrimersforPCRamplification Sizeof Amplified amplified Primername region Forwardprimersequence Reverseprimersequence fragment CAA-imb1a-E1 bmp6a ATCACAGACTGGATGTAAAGATGA CATGAAGAGGGGTGCTGAA 250bp exon1 (SEQIDNO:23) (SEQIDNO:24) CAA-imb1a-E2 bmp6a GTTCGTCAAGCACACTAAAACT CTTTCTGGATGTTCCCCAA 250bp exon2 (SEQIDNO:25) (SEQIDNO:26) CAA-imb1b-E1 bmp6b TCCAGCATGAAGAGGGGTG CTGGATGTAAAGATGACAAGCA 247bp exon1 (SEQIDNO:27) (SEQIDNO:28) CAA-imb1b-E2 bmp6b ACAGGCAGCAGTTGACTTG GCTTTCACCACAGAGAAGACA 322bp exon2 (SEQIDNO:29) (SEQIDNO:30)

[0051] Embodiment 3: The difference between this embodiment and Embodiment 1 or 2 is that: the sgRNA forward primer shown in Table 1 and the sgRNA reverse primer with the sequence of 5-GATCCGCACCGACTCGGTGCCACTTTTTCAAGTTGATAACGGACTAGCCTT ATTTTAACTTGCTATTTCTAGCTCTAAAAC-3 (SEQ ID NO:31) are used to synthesize sgRNA in vitro. The sgRNA synthesis volume in vitro is: 2.5 L of 10 M sgRNA forward and reverse primers, 25 L of 2Dream Taq Master Mix (ThermoFisher, CA, USA), and sgRNA synthesis volume is supplemented to 50 L with enzyme-free water; a total of 100 L is amplified in 2 tubes. The sgRNA synthesis procedure in vitro is: denaturation at 95 C. for 3 min, followed by 30 cycles of 95 C. for 30 s, 58 C. for 30 s, 72 C. for 30 s; extension at 72 C. for 5 min; then the obtained amplified product is purified and recovered, followed by an in vitro transcription, and 30 L of reaction volume is established for each target site: 1 g of sgRNA PCR recovered product, 10 L of NTP Buffer Mix, 2 L of T7 RNA Polymerase Mix, the reaction volume is supplemented to 30 L with enzyme-free water; transcription was performed at 37 C. for 4 h, 20 L of enzyme-free water is added after the reaction, and 2L of DNase I was added after mixing. Digestion at 37 C. for 15 min to remove DNA. Others are the same as the embodiments 1 or 2.

[0052] The in vitro PCR products of sgRNA are detected by 1.5% agarose gel electrophoresis, and the length of products is 120 bp. After detection, the PCR products are purified and recovered with a PCR product purification and recovery kit (Exygen), and the concentration is determined with a Qubit 3 kit (Thermo Fisher, CA, USA) for use. The recovery concentration of sgRNA PCR products is in a range of 100 to 160 ng/L.

[0053] The in vitro transcribed sgRNAs are purified and recovered with an RNA purification kit (Qiagen). The concentration of the recovered product is determined by using a Qubit 3 kit (Thermo Fisher, CA, USA) and the recovered product is stored in a 80 C. refrigerator until use. The recovery concentration of in vitro transcribed sgRNA is in a range of 800 to 3000 ng/L.

[0054] Embodiment 4: The difference between this embodiment and Embodiment 1 or 2 or 3 is that the microinjection method in this embodiment is: the sgRNA synthesized in vitro and Cas9 protein (NEB M0646, MA, USA) are mixed in a molar ratio of 3:1, a mixture obtained is incubated at room temperature for 10 min, then 25% phenol red is added, and the obtained product is injected into Carassius auratus embryos in the single cell stage; among them, the target site sgRNAs on each exon are mixed in equal amounts and injected with a final concentration no less than 50 ng/L of each sgRNA, the control group is injected with 25% phenol red, and the injection volume of each fertilized egg is 1 nL0.02 nL. Others are the same as the Embodiment 1, 2 or 3.

Example 1

[0055] Breeding method of Carassius auratus strain without intermuscular bone: The target sites were designed to knock out the two copies of bmp6 gene in the Carassius auratus genome, bmp6a and bmp6b, respectively, and then F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation were obtained through two rounds of gene knockout and screening, and then the F2 generation individuals from the homozygous line with bmp6a and bmp6b double gene mutation were propagated to form a new Carassius auratus strain without intermuscular bones;

[0056] Among them, the sgRNAs corresponding to (bmp6a), (bmp6b) or (bmp6a and bmp6b) were mixed with Cas9 protein, and the obtained mixture was microinjected into Carassius auratus embryos in the single cell stage, the first round of gene knockout was performed to construct the F0 generation population, and then the F0 generation population was cultured for 3 to 5 months followed by PIT labeling and DNA extraction, then the alleles and mutation rates of somatic mutations of Carassius auratus was determined by sequencing, and the F0 generation individuals with the somatic mutation rate of more than 95% were selected as parents to prepare 0-generation fertilized eggs;

[0057] The sgRNAs corresponding to (bmp6a), (bmp6b) or (bmp6a and bmp6b) were mixed with Cas9 protein, and the obtained mixture was microinjected into the 0-generation fertilized eggs, the second round of gene knockout was performed to construct the F1 generation population, and then the F1 generation population was cultured for 3 to 5 months followed by PIT labeling and DNA extraction, then the alleles and mutation rates of the somatic mutations of Carassius auratus were determined by sequencing, and the F1 generation individuals from a somatic bmp6a and bmp6b double gene mutant line with the somatic mutation rate of more than 95% were selected as parents for reproduction to construct a F2 generation, then F2 generation individuals from a homozygous line with bmp6a and bmp6b double gene mutation were selected.

[0058] Among them, the DNA extraction and sequencing method was as follows: 100 L of lysis solution was added (the composition of the lysate is 0.5 mg/mL proteinase K, 10 mM Tris (pH 8.0), 50 mM KCl, 0.3% Tween 20, 0.3% NP40) to 0.3-0.5 cm.sup.2 fin clip samples; the lysis methods was: the PCR instrument was set at 55 C. for 6 h, 98 C. for 10 min for lysis; after lysis, the obtained product was mixed well, then centrifuged at 1000-2000 rpm for 2 min, and the supernatant was taken as the PCR amplification template. PCR amplification: 2 pairs of primers were used for PCR amplification, of which, 1 pair of primers was the target site primers (the primer sequences used to amplify the target fragment as shown in Table 2 and the 5-end M13 universal primer sequence), and the other pair of primers were index primers (the M13 universal primer sequence and the index base sequence at the 5 end used to identify the sample). PCR amplification was performed in 2 steps. The PCR amplification volume in the first step was 10 L, including 1 L of genomic DNA supernatant, 0.5 L of 1 M forward and reverse target primers, and 50 L of 2Dream Taq Master Mix (Thermo Fisher, CA, USA), the volume was supplemented to 10 L with enzyme-free water; PCR program: 95 C. for 3 min; 10 cycles of 95 C. for 30 s, 60 C. for 30 s and 72 C. for 30 s; 72 C. for 2 min; PCR amplification in step 2 is that 0.5 L of 5 M forward and reverse index primers were added to the PCR amplification product obtained in step 1, the PCR program was set to 95 C. for 2 min; 6 cycles of 95 C. for 30 s, 58 C. for 30 s and 72 C. for 30 s; 15 cycles of 95 C. for 30 s, 72 C. for 30 s; 72 C. for 2 min. After PCR amplification, the PCR products were mixed in equal amounts, and a DNA sequencing library was constructed with the TrueSeq Sample Preparation Kit, and was sequenced by using Illumina MiSeq sequencing platform with 300 bp Pair-End mode. The data obtained by high-throughput sequencing were first used to identify amplification loci and samples using the method described by Tong et al (Tong et al. 2018). The sequenced data of each sample in each PCR product was analyzed using the CRISPResso2 (Fiume et al. 2019) program to determine the alleles and mutation rates of somatic mutations.

[0059] Among them, Sanger sequencing was used for mutant detection: PCR amplification was performed on the samples using the primers shown in Table 2, and the system was set to 25 L: 2 L of genomic DNA supernatant, 1 L of 10 M forward and reverse target primers, 12.5 L of 2Dream Taq Master Mix (Thermo Fisher, CA, USA), the volume was supplemented to 25 L with enzyme-free water; PCR program was 95 C. for 3 min; 35 cycles of 95 C. for 30 s, 60 C. for 30 s, 72 C. for 30 s; 72 C. for 5 min. PCR products were detected by 1.5% agarose gel electrophoresis. After detection, the products were purified and recovered, TA clone was performed, 2 L of PCR products were taken after colony PCR amplification, and detected with 8% polypropylene, and the colonies with different size bands from the control were selected for Sanger sequencing.

[0060] Among them, the knockout target sites of the bmp6 gene were shown in Table 1.

[0061] The sgRNA forward primer was shown in Table 1 and the sgRNA reverse primer with the sequence of 5-GATCCGCACCGACTCGGTGCCACTTTTTCAAGTTGATAACGGACTAGCCTT ATTTTAACTTGCTATTTCTAGCTCTAAAAC-3 (SEQ ID NO:31) were used to synthesize sgRNA in vitro. The sgRNA synthesis volume in vitro was: 2.5 L of 10 M sgRNA forward and reverse primers, 25 L of 2Dream Taq Master mix (Thermo Fisher, CA, USA), the volume was supplemented to 50 L with enzyme-free water; a total of 100 L was amplified in 2 tubes. The sgRNA synthesis procedure in vitro was: denaturation at 95 C. for 3 min, followed by 30 cycles of 95 C. for 30 s, 58 C. for 30 s, 72 C. for 30 s; extension at 72 C. for 5 min; then the obtained amplified product was purified and recovered, followed by an in vitro transcription, and 30 L of reaction volume was established for each target site: 1 g of sgRNA PCR recovered product, 10 L of NTP Buffer Mix, and 2 L of T7 RNA Polymerase Mix, the volume was supplemented to 30 L with enzyme-free water; transcription was performed at 37 C. for 4 h, 20 L of enzyme-free water was added after the reaction, and 2 L of DNase I was added after mixing. Digestion at 37 C. for 15 min to remove DNA.

[0062] In vitro sgRNA PCR products were detected by 1.5% agarose gel electrophoresis, and the length of products was 120 bp. After detection, the sgRNA PCR products were purified and recovered with a PCR product purification and recovery kit (Exygen), and the concentration was determined with a Qubit 3 kit (Thermo Fisher, CA, USA) for use. The recovery concentration of the sgRNA PCR products was in a range of 100 to 160 ng/L.

[0063] The in vitro transcribed sgRNAs were purified and recovered with an RNA purification kit (Qiagen). The concentration of the recovered product was determined by using a Qubit 3 kit (Thermo Fisher, CA, USA) and the recovered product was stored in a 80 C. refrigerator until use. The recovery concentration of in vitro transcribed sgRNA was in a range of 800 to 3000 ng/L.

[0064] The microinjection method was as follows: sgRNA synthesized in vitro and Cas9 protein (NEB M0646, MA, USA) were mixed at a molar concentration ratio of 3:1, a mixture obtained was incubated at room temperature for 10 min, then 25% phenol red was added, and the obtained product was injected into Carassius auratus embryos in the single cell stage; among them, the target site sgRNAs on each exon were mixed in equal amounts and injected with a final concentration no less than 50 ng/L of each sgRNA, the control group was injected with 25% phenol red, and the injection volume of each fertilized egg is 1 nL0.02 nL.

[0065] The diploid Carassius auratus used in this example came from the Hulan Experiment Station of the Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences. The parents fish were cultured in a 700 m.sup.2 fish pond, and were injected with oxytocin drugs (4 g/kg luteinizing releasing hormone analogue, LRH-A2; 1 mg/k DOM, 100 units/kg chorionic gonadotropin HCG) for spawn. After artificial fertilization, the obtained fertilized eggs were incubated in a hatching tank of 24 cm20 cm24 cm, and the hatching water temperature was 22 C.-23 C. After hatching, larvae were fed with Artemia salina 4 times/d after the yolk was absorbed. The larvae were moved to an outdoor 500 m.sup.2 pond to culture when they were cultured in a hatching tank to the size of 2 cm-3 cm. Artificial diet was fed in the fish fingerling breeding stage, 2 times/d.

[0066] The bones of the new Carassius auratus strain with intermuscular bone-deficient in this example were observed by bone-staining, and it was found that the intermuscular bones were completely deficient in 5 samples, and the number of intermuscular bones in 4 samples was only 4 to 20 (as shown in FIG. 1 to 4, 9 and 10). The mutant alleles of above samples in exon 1 of bmp6a and exon 1 of bmp6b were detected, and it was found that mutations were observed in all three target sites of exon 1 (as shown in FIG. 5A to C and FIG. 6A to C). Pre-mature of stop codons in proteins was caused by gene knockout (as shown in FIGS. 7 and 8).

[0067] Although the present disclosure has been described in detail through the above example, it is only a part of embodiments of the present disclosure, rather than all embodiments. People can also obtain other embodiments according to the present embodiment without creativity. These embodiments all fall within the claimed scope of the present disclosure.