New approach for fish breeding by chemical mutagenesis: establishment of TILLING method in fugu (Takifugu rubripes) with ENU mutagenesis
- Miwa Kuroyanagi†1,
- Takashi Katayama2,
- Tadashi Imai2,
- Yoshihisa Yamamoto2,
- Shin-ichi Chisada3,
- Yasutoshi Yoshiura3,
- Tomokazu Ushijima4,
- Tomonao Matsushita4,
- Masashi Fujita1,
- Aoi Nozawa1,
- Yuzuru Suzuki1,
- Kiyoshi Kikuchi1Email author and
- Hiroyuki Okamoto†5Email author
© Kuroyanagi et al.; licensee BioMed Central Ltd. 2013
Received: 8 February 2013
Accepted: 4 November 2013
Published: 13 November 2013
In fish breeding, it is essential to discover and generate fish exhibiting an effective phenotype for the aquaculture industry, but screening for natural mutants by only depending on natural spontaneous mutations is limited. Presently, reverse genetics has become an important tool to generate mutants, which exhibit the phenotype caused by inactivation of a gene. TILLING (Targeting Induced Local Lesions IN Genomes) is a reverse genetics strategy that combines random chemical mutagenesis with high-throughput discovery technologies for screening the induced mutations in target genes. Although the chemical mutagenesis has been used widely in a variety of model species and also genetic breeding of microorganisms and crops, the application of the mutagenesis in fish breeding has been only rarely reported.
In this study, we developed the TILLING method in fugu with ENU mutagenesis and high-resolution melting (HRM) analysis to detect base pair changes in target sequences. Fugu males were treated 3 times at weekly intervals with various ENU concentrations, and then the collected sperm after the treatment was used to fertilize normal female for generating the mutagenized population (F1). The fertilization and the hatching ratios were similar to those of the control and did not reveal a dose dependency of ENU. Genomic DNA from the harvested F1 offspring was used for the HRM analysis. To obtain a fish exhibiting a useful phenotype (e.g. high meat production and rapid growth), fugu myostatin (Mstn) gene was examined as a target gene, because it has been clarified that the mstn deficient medaka exhibited double-muscle phenotype in common with MSTN knockout mice and bovine MSTN mutant. As a result, ten types of ENU-induced mutations were identified including a nonsense mutation in the investigated region with HRM analysis. In addition, the average mutation frequency in fugu Mstn gene was 1 mutant per 297 kb, which is similar to values calculated for zebrafish and medaka TILLING libraries.
These results demonstrate that the TILLING method in fugu was established. We anticipate that this TILLING approach can be used to generate a wide range of mutant alleles, and be applicable to many farmed fish that can be chemically mutagenized.
KeywordsTILLING Fugu ENU HRM Myostatin Mutagenesis Fish breeding
Mutations are the basis of genetic variation and mutant populations are indispensable genetic resources in all organisms. This variation can be either naturally occurring or, in plants, animals and lower organisms, induced by chemical or physical treatments. Mutation induction, for example, has played an important role in the genetic improvement of crop species that are of economic significance.
In breeding of aquacultural fish species, it is a final goal to discover and generate the fish with various useful phenotypes for supplying to consumers. Currently genetic breeding is carried out with the search for natural variation in traits such as growth, productivity and disease resistance, and then improved strains are produced by genetic selection or marker-assisted breeding approaches . Thus the success of fish breeding solely depends on the mutations pre-existing in the focal population. The rate of natural spontaneous mutations in fish species is generally lower than 1.0×10-6 at specific loci. On the other hand, mutation frequency caused by radiation and chemical mutagenesis range from 1.0×10-3 to 3.9×10-3 at specific loci when in vivo spermatogonial treatments were carried out in zebrafish and medaka [2–4]. Given the high frequency of the mutation ratio, it is reasonable to consider that the chemical mutagenesis will be an efficient way to produce new mutants for future genetic improvement in aquacultural species.
TILLING (Targeting Induced Local Lesions IN Genomes) is a reverse genetics strategy that identifies induced mutations in specific genes of interest in chemically mutagenized populations with high-throughput discovery technologies. TILLING, first described in 2000 for mutation detection in Arabidopsis, is now used in a wide range of plants including soybean , rice , barley  and maize  as well as for animal model systems, including Drosophila, Caenorhabditis elegans, rat , medaka , zebrafish  and for the discovery of naturally occurring polymorphisms in humans .
The first step in TILLING is chemical mutagenesis. In plants, mutations have been induced by using ethylmethane sulfonate (EMS) and methylnitrosourea (MNU) [5, 7]. In animals, EMS is used in many invertebrate species for mutagenesis, whereas N-ethyl-N-nitrosourea (ENU) is used for males of vertebrates, and is the most widely employed method . ENU acts as an alkylating agent and transfers its ethyl group to nucleophilic nitrogen or oxygen sites on deoxyribonucleotides, leading to base mismatch during DNA replication. Single-base substitutions that resemble natural spontaneous mutations [17, 18], are mainly induced by ENU treatment and all genes are mutated at random . Because this protocol induces mutations in germ-line stem cells, mutations are transmitted at the rate for at least several months after the treatment. Also, the mutations are induced with premeiotic ENU treatment, so that the initial mutation, usually an ethylation of a base on one DNA strand, is fixed by DNA replication prior to production of differentiated sperm, leading to non-mosaic offspring in the next generation . On the other hand, it has been recently reported that some point mutations were induced by postmeiotic mutagenesis in Ctenopharyngodon idellus (grass carp) with ENU for mature sperm . In zebrafish, the postmeiotic ENU mutagenesis resulted in mosaic progeny, but produced a 10 fold increase in the frequency of induced mutations . This protocol in zebrafish has been also shown to induce mutations such as deletions and translocations triggered by chromosomal rearrangements . These results indicate that postmeitotic ENU treatment of male germ cells can induce a variety of mutations from point mutations to deletions and translocations.
In zebrafish and medaka, which are model fish, mutations are introduced in spermatogonia by soaking founder fish in ENU solution [13, 22]. However, this protocol is difficult to apply to many aquaculture species including fugu, as the body size of these fish is markedly larger than that of medaka and zebrafish. Thus alternative protocols are required for fish with a large body size.
After the ENU treatment, a large F1 population is generated to choose the mutations for purpose from many random heterozygous mutations in their genomes. As the next step in TILLING, the genomic DNA from these animals containing induced mutations is screened with discovery technologies. Sequencing is considered the standard for DNA-based mutation detection because it reveals the exact location and the type of mutations. Taniguchi et al. employed a direct sequencing approach even for the prescreening step . CEL I nuclease-mediated screening of heteroduplexes are widely carried out to detect mutations in TILLING [23, 24]. Both methods are effective, but the costs are substantial and it takes a relatively long time for the screening. Recently, high resolution melting (HRM) analysis has succeeded to detect point mutations in plant populations , medaka TILLING library  and human disease . HRM analysis makes use of characteristics in DNA as follows; the thermal stability of a DNA fragment is determined by its base sequence. When the DNA fragment contains an altered sequence, the duplex stability is changed, leading to different melting behavior. Thus changes in the sequence within the DNA fragment such as point mutation can be detected by the melting analysis.
To apply TILLING technology for breeding practices, it is critical to select a proper target gene in which mutations are likely to cause useful phenotypes such as rapid growth and disease resistance. One of the growth-related factors, Myostatin (Mstn) is a member of the transforming growth factor-β (TGF-β) superfamily, and functions as a negative regulator of skeletal muscle mass . Naturally occurring mutations in the Mstn gene have been identified in double-muscling cattle breeds [29, 30], mice , dogs  and even humans . Recently, we have demonstrated that the antagonistic role of Mstn against muscle growth is conserved from mammals to fish species . Indeed mstn-deficient medaka (mstnC315Y), which were generated by the TILLING method, show a significant increase in body weight in comparison with wild-type medaka, due to over-growth of muscles. Therefore, Mstn is a promising target gene for mutation breeding in fishes.
Here, in order to initiate targeted mutation breeding in fish improvement, we searched for efficient methods for TILLING combining ENU mutagenesis with HRM analysis in Takifugu rubripes, tiger pufferfish (fugu). Fugu males were treated with ENU via intraperitoneal injection each week for three weeks, and the sperm collected from the founder fish was used for in vitro fertilization to obtain the progeny. We then isolated genomic DNA from the progeny, and carried out HRM analysis targeting the exon 3 of fugu Mstn gene. Our screening successfully identified ten mutations including one nonsense mutation in fugu Mstn gene. The mutation frequency in fugu Mstn gene was similar to that identified in the medaka TILLING library, demonstrating that we were able to establish an effective protocol of the TILLING method in fugu. This protocol is likely to be applicable to various species of aquacultured fish.
Development of the ENU mutagenesis protocol for fugu TILLING
Effects of various ENU concentrations on fertilization and hatching ratios of fugu
ENU conc. (mg/kg body weight)
Fertilization ratio (%)
Hatching ratio (%)
Effects of gonadotropin treatment on timing of sperm collection and the sperm volume in fugu
Period after the last ENU treatment
ENU conc. (mg/kg)
No. of fish from which sperm was collected/No. of fish examined (%)*
Sperm volume (ml)
No. of fish from which sperm was collected/No. of fish examined (%)*
Sperm volume (ml)
No. of fish from which sperm was collected/No. of fish examined (%)*
Sperm volume (ml)
Gonadotropin treatment (2 weeks after the last ENU treatment)
Gonadotropin treatment (3 weeks after the last ENU treatment)
Non-gonadotropin treatment (control)
Target gene for generation of a fugu mutant with TILLING
Detection of ENU-induced mutations in fugu with HRM analysis and DNA sequencing
Summary of fugu MSTN mutations detected by HRM analysis
ENU conc. (mg/kg body weight)
Amino acid substitution
CTGGGACTGG(A > G)CTGCGACGAG
TTGCAGAAAT(A > G)CCCGCACACC
gCAACCTTTC(A > G)TGGAGGTGAA
GAGACCTGGG(A > G)CTGGACTGCG
CACCTGGTGA(A > G)CAAGGCCAAC
AGTGTGAGTA(C > A)ATGCACTTGC
CGCACACCCA(C > T)CTGGTGAACA
AGGGACCGCG(G > A)GCCCTTGCTG
CATGCTCTAC(T > C)TCAACCAAGA
CTCTACTTCA(A > G)CCAAGAACAG
Identification of Mstn-deficient mutation in fugu
ENU-induced mutation ratio in fugu
The mutation ratio in the mutagenesis is a key point to produce the mutants effectively. Therefore, the ratio was calculated using data detected in fugu Mstn gene by HRM analysis. We screened about 3.0 Mb in total and identified 10 independent mutations (Table 3). Therefore, the ENU-induced mutation ratio for fugu Mstn gene was one substitution per 297 kb.
Establishment of fugu TILLING
Recently, in vitro ENU mutagenesis of postmeiotic sperm was performed in grass carp . In their study, the induced mutations were investigated using growth-related genes (Igf-2a, Igf-2b, Mstn-1, Mstn-2, Fst-1 and Fst-2) and the estimated mutation ratio was around 0.41%. However, the mutations are apparently mosaic in F1 individuals. Mosaicism of the induced mutation was also reported in zebrafish when postmeiotic sperm was mutagenized . This situation makes recovery of the mutant line inefficient in contrast to the case of non-mosaic mutation induced by ENU-treatment of spermatogonia . The mosaicism in F1 individuals and resulting non-Mendelian transmission from F1 to F2 generations hinder some screening protocols. In addition, large-scale rearrangements of the genomic structure such as translocation of genomic segments have been reported following the mutagenesis of postmeiotic sperm . Therefore, the suitability of the in vitro ENU mutagenesis of postmeiotic sperm for mutation breeding in fish improvement remains to be verified. Although the in vitro ENU treatment is more convenient than the in vivo methods used in fugu and other model animals, it will be required to examine whether the induced mutations will be inherited by the progeny in the mutants.
In TILLING, high efficiency for ENU-induced mutations is required to obtain various phenotypes in gene of interest. A previous report indicated that the induced mutation frequency increased with ENU dose dependently in mouse . In this study, there was no difference in the number of surviving fish, fertilization and hatching ratios in ENU treatment with 7, 70 mg/kg (Table 2). ENU treatment with 120 mg/kg affected slightly on the number of surviving treated fugu. However, the number of the induced mutations did not tend to increase in each ENU concentration between 7 mg/kg and 120 mg/kg (Table 3). These results suggested that each ENU concentration in the method may become saturation in premeiotic sperm (spermatogonia and spermatocyte cells) of fugu.
Identification of fugu Mstn mutant
It has been reported that loss of the Mstn function in medaka and mice causes double-muscling phenotype with both hyperplasia and hypertrophy [28, 33]. In the present study, ENU-induced mutations in fugu Mstn were confirmed and in particular, the nonsense mutation (Y316X) was identified at C-terminal region of the peptide. The mstn mutant medaka (mstnC315Y) made by using the TILLING method  corresponded to a naturally occurring mutation (C313Y) in double-muscled Piedmontese cattle  (Figure 2B). The cysteine residue at the 315-position exists in C-terminal region of the Mstn peptide (Figures 2B and 5D) and is known to be one of the important amino acid residues involved in forming the cysteine knot structure in members of the TGF-β superfamily . The structure functions as a homodimer to elicit a biological function . Therefore, we speculate that the structural formation composed of six essential cysteines in fugu Mstn will be non functional by the truncated form of the peptide resulting from a point mutation at 317 of Tyr. Eventually, the mutation will be conceivable to cause loss of the Mstn function in fugu Mstn mutant. The double-muscling phenotype in medaka and bovine could be similarly expected for an Mstn-deficient mutant fugu.
Development of HRM analysis
For screening, there are a variety of mutation discovery technologies such as direct DNA sequencing, CEL I assay and HRM analysis [13, 23, 26]. We introduced the HRM analysis for screening system in fugu because there is no albino mutant in fugu for specific locus test (SLT) and the analysis has recently succeeded to detect point mutations as a screening step for the TILLING approach [25, 26]. But, there was room for improvement with respect to the cost and variability of output data. In this study, we have developed a system of HRM analysis as follows: cost-cutting with fluorescent intercalator; SYTO 9, introduction of the nested PCR for precision enhancement and avoidance of variability, speed-up with 384 samples treatment in a single experiment, and high uniformity of dispensed volume and removal of human error by using the automatic dispense device. As a result of the improvement, some mutations could be detected efficiently in the analysis (Figure 4 and Table 3). It is difficult for HRM analysis to investigate mutations with a long amplicon in comparison with the other methods, but the acquisition of melting profiles in HRM is not dependent on the location and number of mutations present in the analyzed fragment. It has been compared about the mutations detected by the screening approaches. The results indicate that HRM analysis is useful as resequencing in medaka  and more efficient than CEL I assay in mammals . Taken together, these demonstrate that HRM analysis is a reliable high-throughput approach for screening ENU-mutagenized individuals in potentially various species.
Mutation ratios in fugu TILLING
The mutation ratio in ENU mutagenesis has been calculated in a SLT using the albino mutant in zebrafish and medaka [21, 41]. The characteristics of the albino mutant in fish reveal mainly red eyes and one of the responsible genes is tyrosinase. In some aquacultured fish with albino mutants, it may be easy for these ENU-treated fish to examine the mutation ratio with SLT. However, albino mutants are not found in nature in all animal species.
Instead of SLT, the mutation ratio was calculated using base-pairs screened and the number of mutations detected in fugu Mstn gene by HRM analysis (Table 3). The frequency in fugu (1 mutation per 297 kb) is similar to the mutation frequency in medaka (1 mutant per 345 kb) and in zebrafish (1 mutant per 235 kb), which were calculated by resequencing data [13, 42]. This suggests that ENU mutagenesis in fugu was efficiently performed in our way and the protocol detailed here will potentially be useful for constructing a TILLING library.
CG to AT (or reverse) changes are often detected in HRM analysis, but the endogenous SNPs are also detected as the same changes and give the most significant difference in this analysis. In our study, there were no endogenous SNPs in the C-terminal region of the Mstn gene in both the investigated family of fugu and their parents (data not shown). Although G > C, C > G, T > A, A > T mutations were not observed in the 10 mutations detected in the Mstn gene (Table 3), it was reasonable results since the appearance frequency of their mutations are also rare in human genome .
Generation of an Mstn-deficient mutant in fugu
Our next step will be to generate the Mstn mutant for the purpose of serving the aquaculture industry. One of the key components for a successful TILLING project is a near-perfect ability to recover valuable mutations once they have been found. For the recovery, TILLING libraries can either be ‘living’ or cryopreserved as suggested by Moens et al. . A living library is one in which mutagenized F1 fish are held in pools in tanks while their fin-tip genomic DNA is screened. A cryopreserved library is one in which each genomic DNA sample corresponds to vials of frozen sperm that can be used in an in vitro fertilization to generate F2 heterozygotes. As a large amount of sperm can be stored from a single adult male of aquacultured fish, in general, the cryopreserved library is more desirable than ‘living’ library. However, generation time of many aquacultured fish including fugu is more than one year, making the construction of the cryopreserved library time-consuming. We are currently generating an Mstn-deficient fugu with ‘living’ library as the next step.
In the present study, the fugu TILLING method was established from the results in both ENU mutagenesis and a mutation search technology, HRM analysis. This is a first report that TILLING is an applicable method for aquacultured fish. The approach for the fish species might play an important role in producing a fish, which exhibits a valuable phenotype (e.g., rapid growth, disease resistance) for the aquaculture industry. We hope that the TILLING technique will be an alternative to the classic selection approach for improving traits in aquacultured fish species.
Fisheries Research Agency and The University of Tokyo rule that ethic approval is not necessary for study of fish. Nevertheless, this project was conducted in accordance with the Regulation for Animal Experiments in the Fisheries Research Agency, which stipulated in the Fundamental Guidelines for Proper Conduct of Animal Experiments and Related Activities in Research Institutions (Ministry of Agriculture, Forestry and Fisheries of Japan), and under the supervision of the Animal Care and Use Committee for the National Research Institute of Aquaculture. In this study, ENU and hormone injections for fugu were carried out under 0.05% anesthetic agent (2-phenoxy-ethanol) and hatching fry of fugu were harvested after they were soaked in sea water containing ice. All efforts were made to minimize suffering.
Experimental fish and eggs for in vitrofertilization
50 fugu males (body weight: around 1 kg) were purchased from a commercial supplier (Shihou Suisan, Minami-ise, Mie). The fish were kept at the Nansei station of National Research Institute until obtaining sperm from the ENU-treated fish. Three-year old fugu females (body weight: around 2–3 kg) were cultivated in 30–50 ton of seawater tanks at the Yashima station of National Center. Two independent fugu females were used to obtain eggs for in vitro fertilization. Mainly, the fertilization and harvest of the hatched fry were performed at the Yashima station. A part of the fertilized-eggs were packed into seawater without air and sent to Fisheries Laboratory attached to the University of Tokyo. The fertilized eggs were incubated at 18°C during 8 to 10 days until harvest of the hatched fry.
Mutagenesis with ENU
All working solutions of ENU (Conc.: 11 mg/ml, Sigma-Ardrich, St. Louis, MO USA) were freshly prepared following the methods described by Taniguchi et al. . ENU was prepared at three concentrations (7, 70, 120 mg per kg, body weight of fugu). For ENU mutagenesis, fugu males (around 1 kg of body weight) were used and injected into the intraperitoneal with the ENU concentrations prepared for 10 mM sodium phosphate buffer (pH 6.3). The injections were repeated 3 times at weekly intervals to the founder fish. After the last ENU injection, the fish were injected intramuscularly with human chorionic gonadotropin (HCG, 500 IU/kg body weight) or salmon pituitary homogenate (SPH, 10 mg/kg body weight), consisting of acetone-dried chum salmon pituitary homogenized in Ringer, was added to HCG because of increasing the success rate of spermiation in ENU-treated males. After the hormone treatment, the sperm was collected from each fish.
Cryopreservation of sperm
The freshly-collected sperm from each fugu male was suspended at the rate of 1 to 4 in freezing diluents [10% dimethyl sulphoxide (Sigma-Aldrich Co., Tokyo, Japan) in 90% fetal bovine serum (Thermo Scientific Co., Tokyo, Japan)] and cryopreserved as described by Ohta et al. . The amount of sperm held in each capillary (0.5 ml vol.) was enough to fertilize more than twenty thousand eggs.
About one million unfertilized eggs were squeezed out from a fugu female, which was induced ovulation by administration of an LHRH analogue (des Gly10 [D-Ala6]-LHRH, 400 μg/kg body weight). Fresh sperm, obtained from each fugu family line, was used in this method, and the sperm activity was examined under a light microscope. 0.4 ml of the sperm was added to 12 g of the eggs (around twenty thousand eggs), and these were mixed well. After adding seawater to activate the sperm for fertilization, the eggs were immediately mixed and stood for 5 minutes. The eggs were washed 3 times with new seawater to remove the extra sperm after discarding the remaining seawater. The fertilized eggs were incubated for 4 to 6 hours at 18°C. The eggs were moved to a hatching tank with running water and kept at 18°C until they hatched. At 8 to 10 days after fertilization, around ten thousand hatched fry were obtained in each family. These fry were harvested in a 50 ml tube containing 99.5% ethanol (EtOH) and stored at −30°C.
The fertilization ratio was determined at the four-cell stage with Sera solution (EtOH 6: formalin 3: acetic acid 1), which makes the fertilized eggs transparent. In addition, one hundred eggs were counted at melanophore stage in the fertilized eggs of each cross, and then incubated at 18°C until hatching. The hatching ratio was calculated as the number of the hatched fry per 100 fertilized eggs.
Cloning of fugu Mstngene
RNA was extracted from wild-type fugu skeletal muscle by RNeasy® Mini Kit (QIAGEN Sciences, MD, USA) according to the manufacturer’s instructions. cDNA was synthesized using PowerScript (BD Biosciences Clontech, Palo Alto, CA, USA). To identify fugu Mstn gene orthologs, we used the basic local alignment search tool (BLAST) to search the fugu genome database (http://www.ensembl.org/Takifugu_rubripes/blastview). Fugu Mstn cDNA sequences was amplified using a reverse-transcription PCR by gene specific primer pairs flanking the ORF: 5′UTR-F1, 5′-TAGAACGCATTTCCACGTCTTC-3′, 3′UTR-R1: 5′-TGGTTTGAGGCCTCTACTCTGG-3′. Amplicons were cloned into a sequencing vector, and full-length ORF sequences were obtained.
Identification of mutations by HRM analysis
Extraction of genomic DNA
In this study, 1,056 hatched fry of fugu were prepared in each family obtained from in vitro fertilization. A single fry was placed into each well of 8-well tubes containing 45 μl of 50 mM NaOH, and incubated at 95°C for 10 min. After incubation, 5 μl of 1 M Tris–HCl (pH 8.0) was added to the lysate immediately for neutralization. These aliquots were centrifuged at 1,100 × g and room temperature for 20 min to remove tissue debris. Then, each 40 μl of supernatant was harvested into 96-deep well plates and stored at −30°C. Although the genomic DNA stocks were crude extract and freeze-thawed repeatedly, we could use the stocks for our assay during 10 months at least.
To find the loss-of-function mutations in the fugu Mstn gene, the third exon of the gene was screened using the primers designed by Primer3 programs on the web (source code available at http://bioinfo.ut.ee/primer3-0.4.0/primer3/). The primer sets and amplicon sizes were as follows, the forward and reverse primers #1: 5′-CAGCCATGATCTAACTTCTGCAAG-3′ and 5′-GTTTGAGGCCTCTACTCTGGTTCT-3′, yielding a fragment of 586 bp; the forward and reverse primers #2: 5′-CTCGCTCACATGCATTCTGTC -3′ and 5′-CTTGTTCACCAGGTGGGTGT-3′, yielding a fragment of 268 bp; the forward and reverse primers #3: 5′-AAGCGCTACAAGGCCAACTA -3′ and 5′-AAAGCGCTGGAACTGGTAGA-3′, yielding a fragment of 281 bp. These primers were used for the 1st-PCR (primer #1), and for the 2nd-PCR (primers #2 and #3) [see PCR].
PCR amplicons for HRM were obtained by using a nested-PCR approach. The genomic DNA stocks and all PCR reaction mixtures were dispensed by a laboratory automation system, Biomek®3000 (Beckman Coulter, CA, USA).
The 1st-PCR (PCR1) was performed in 96-well PCR plates in 15 μl volumes. Reactions included 2 μl of the genomic DNA solution in 1x KOD-FX PCR buffer (Toyobo, Tokyo, Japan), with 400 μM of deoxynucleotide triphosphate, 0.1 U of KOD-FX polymerase (Toyobo), and 300 nM of each primer. The polymerase was able to increase precise amplicons with crude samples. The plates were centrifuged at 200 × g for 1 min and then PCR was performed in T-Professional Thermocycler (Biometra, Germany). The cycling conditions were set as follows: one cycle of 94°C for 2 min, followed by 40 cycles of 98°C for 10 sec; 60°C for 30 sec; 68°C for 40 sec. The PCR1 products were firstly diluted 10 fold for DNA sequencing and then finally 140 fold with sterilized water by using Biomek® 3000. The 140 fold diluted samples were used as a template for second step of PCR reactions.
The 2nd-PCR (PCR2) was performed in 384-well PCR plates (Roche Diagnostics, Basel, Switzerland) in 9 μl volumes. Reactions included 2 μl of the diluted PCR1 products in 1x Hot-Start Gene Taq PCR buffer (Nippon Gene, Tokyo, Japan), with 200 μM of deoxynucleotide triphosphate, 1.25 μM of SYTO 9 (Molecular Probes, OR, USA) as a fluorescent substance, 0.16 U of Hot-Start Gene Taq polymerase (Nippon Gene), and 200 nM of each primer. The plates were centrifuged at 200 × g for 1 min and then PCR was performed in a 384-well Light Cycler 480 System (Roche Diagnostics). The cycling conditions were set as follows: one cycle of 95°C for 5 min, followed by 25 cycles of 95°C for 10 sec; 60°C for 30 sec; 72°C for 20 sec.
Several samples of the PCR1 and PCR2 products were examined with a microchip electrophoresis system, MultiNA (Shimadzu, Kyoto, Japan) whether there are the extra major fragments except for the proper fragment or not, because non-specific PCR products reduce the performance of HRM assay.
After PCR, the plates were imaged in the 384-well Light Cycler 480 System (Roche Diagnostics). The PCR products in the plates were heated at 0.1°C/second, and fluorescence measurements were collected from 65°C to 98°C. Briefly, background-subtracted fluorescence curves were normalized between 0% and 100%. The curves were viewed on the subtraction plots to magnify differences in the shape of the melting curves. The subtraction plots are useful to cluster the samples into groups. Clustering of the melting curves for genotype identification was carried out using the Light-Cycler software.
A single base-pair change within the probe region causes a significant shift in melting temperature (Tm) of the amplicon. In the TILLING method, all the induced mutations are heterozygous. PCR amplification of heterozygotes followed by heat denaturing and annealing result in the formation of 4 duplexes: two homoduplexes and two heteroduplexes. Each duplex has a characteristic melting temperature, and the sum of all transitions can be detected by melting curve analysis. The heterozygous samples are identified by differences in shape of the melting curve.
1,056 of the genomic DNA in each fugu family were subjected to PCR and two amplicons (268 bp and 281 bp) for screening were examined in HRM. The genomic DNA showing melting curves different from those of the wild-type allele were selected as mutant candidates. As previously reported in the medaka TILLING method , a substantial number of SNPs might exist in the fugu genomic DNA used. SNPs were discriminated from true mutations by defining the same type of base change detected in more than two sequencing reactions as an SNP. Moreover, the genomic DNA in the pair of fugu parents extracted from the cut fins was examined by DNA sequencing to confirm whether it was the true mutations or not.
The positive genomic DNA were screened as candidate mutations by HRM analysis. The PCR1 products of the genomic DNA were already diluted 10 fold, and stocked at −30°C [see PCR]. These samples were used as a template for direct DNA sequencing. Sequencing reactions were carried out using a BigDye® terminator ver. 3.1 cycle sequencing kit (Applied Biosystems, CA, USA) and each forward and reverse primers. PCR cycling condition was set as follows: one cycle of 96°C for 60 seconds, followed by 25 cycles of 96°C for 10 seconds; 50°C for 5 seconds; 60°C for 120 seconds. Sequencing products cleaned up by magnetic bead-based purification system (Agencourt CleanSEQ, Beckman Coulter) were run on automated ABI 3130xl DNA analyzers (Applied Biosystems). Sequences were analyzed for the presence of heterozygous mutations using PolyPhred  and manual inspection of the mutated positions. All candidate mutations were confirmed with new PCR1 products and resequencing read. The PCR1 products were treated with an ExoSAP-IT kit (Affymetrix, Ohio, USA) to remove unused primers and nucleotides.
CEL I assay
To confirm the presence of induced mutations in the target sequence, CEL I assay was performed with genomic DNA containing the mutations detected by direct sequencing. CEL I enzyme was purified from celery according to Anai et al. . The enzyme activity for each batch of CEL I was determined experimentally using control samples.
PCR was performed in 9 μl volumes, and the primer sets and amplicon size were as follows, the forward and reverse primers: 5′-CTCGCTCACATGCATTCTGTC-3′ and 5′-AAAGCGCTGGAACTGGTAGA-3′, yielding a fragment of 465 bp. Reactions included 2 μl of the diluted PCR1 products in 1x Hot-Start Gene Taq PCR buffer (Nippon Gene), with 200 μM of deoxynucleotide triphosphate, 0.16 U of Hot-Start Gene Taq polymerase (Nippon Gene), and 200 nM of each primer. The cycling condition was set as follows: one cycle of 95°C for 5 min, followed by 25 cycles of 95°C for 10 sec; 60°C for 30 sec; 72°C for 40 sec.
Directly following the PCR, heteroduplex was formed by incubating at 98°C for 8 min, 80°C for 20 sec, 70 cycles of 80°C for 7 sec with a decrement of 0.3 sec per cycle. Next, PCR reactions were incubated with 1 μl CEL I enzyme solution in a total volume of 10 μl at 37°C for 60 min. CEL I reactions were stopped by adding 10 mM EDTA. The cleavage products were detected with MultiNA (Shimadzu).
Targeting Induced Local Lesions In Genomes
We would like to thank Miwa Osa (NRIA) for the excellent technical assistance during the experiment. This work was supported by the Program for Promotion of Basic and Applied Researches for Innovations in Bio-oriented Industry (BRAIN) and Grant-Scientific Research from Japanese Society for Promotion of Science.
- Moen T, Baranshi M, Sonesson AK, Kjøglum S: Confirmation and fine–mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): population-level associations between markers and trait. BMC Genomics. 2009, 10: 368-381. 10.1186/1471-2164-10-368.PubMed CentralView ArticlePubMedGoogle Scholar
- Driever W, Solnica-Krezel L, Schier AF, Neuhauss SC, Malicki J, Stemple DL, Stainier DY, Zwartkruis F, Abdelilah S, Rangini Z, et al: A genetic screen for mutations affecting embryogenesis in zebrafish. Development. 1996, 123: 37-46.PubMedGoogle Scholar
- Shima A, Shimada A: Development of a possible nonmammalian test system for radiation-induced germ-cell mutagenesis using a fish, the Japanese medaka (Oryzias latipes). Proc Natl Acad Sci U S A . 1991, 88 (6): 2545-2549. 10.1073/pnas.88.6.2545.PubMed CentralView ArticlePubMedGoogle Scholar
- Loosli F, Koster RW, Carl M, Kuhnlein R, Henrich T, Mucke M, Krone A, Wittbrodt J: A genetic screen for mutations affecting embryonic development in medaka fish (Oryzias latipes). Mech Dev. 2000, 97 (1–2): 133-139.View ArticlePubMedGoogle Scholar
- McCallum CM, Comai L, Greene EA, Henikoff S: Targeted screening for induced mutations. Nat Biotechnol. 2000, 18 (4): 455-457. 10.1038/74542.View ArticlePubMedGoogle Scholar
- Cooper JL, Till BJ, Laport RG, Darlow MC, Kleffner JM, Jamai A, El-Mellouki T, Liu S, Ritchie R, Nielsen N, et al: TILLING to detect induced mutations in soybean. BMC Plant Biol. 2008, 8: 9-10.1186/1471-2229-8-9.PubMed CentralView ArticlePubMedGoogle Scholar
- Suzuki T, Eiguchi M, Kumamaru T, Satoh H, Matsusaka H, Moriguchi K, Nagato Y, Kurata N: MNU-induced mutant pools and high performance TILLING enable finding of any gene mutation in rice. Mol Genet Genomics. 2008, 279 (3): 213-223. 10.1007/s00438-007-0293-2.View ArticlePubMedGoogle Scholar
- Slade AJ, Fuerstenberg SI, Loeffler D, Steine MN, Facciotti D: A reverse genetic, nontransgenic approach to wheat crop improvement by TILLING. Nat Biotechnol. 2005, 23 (1): 75-81. 10.1038/nbt1043.View ArticlePubMedGoogle Scholar
- Till BJ, Reynolds SH, Weil C, Springer N, Burtner C, Young K, Bowers E, Codomo CA, Enns LC, Odden AR, et al: Discovery of induced point mutations in maize genes by TILLING. BMC Plant Biol. 2004, 4: 12-10.1186/1471-2229-4-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Bentley A, MacLennan B, Calvo J, Dearolf CR: Targeted recovery of mutations in Drosophila. Genetics. 2000, 156 (3): 1169-1173.PubMed CentralPubMedGoogle Scholar
- Gilchrist EJ, O’Neil NJ, Rose AM, Zetka MC, Haughn GW: TILLING is an effective reverse genetics technique for Caenorhabditis elegans. BMC Genomics. 2006, 7: 262-10.1186/1471-2164-7-262.PubMed CentralView ArticlePubMedGoogle Scholar
- Smits BM, Mudde J, Plasterk RH, Cuppen E: Target-selected mutagenesis of the rat. Genomics. 2004, 83 (2): 332-334. 10.1016/j.ygeno.2003.08.010.View ArticlePubMedGoogle Scholar
- Taniguchi Y, Takeda S, Furutani-Seiki M, Kamei Y, Todo T, Sasado T, Deguchi T, Kondoh H, Mudde J, Yamazoe M, et al: Generation of medaka gene knockout models by target-selected mutagenesis. Genome Biol. 2006, 7 (12): R116-10.1186/gb-2006-7-12-r116.PubMed CentralView ArticlePubMedGoogle Scholar
- Sood R, English MA, Jones M, Mullikin J, Wang DM, Anderson M, Wu D, Chandrasekharappa SC, Yu J, Zhang J, et al: Methods for reverse genetic screening in zebrafish by resequencing and TILLING. Methods. 2006, 39 (3): 220-227. 10.1016/j.ymeth.2006.04.012.View ArticlePubMedGoogle Scholar
- Till BJ, Zerr T, Bowers E, Greene EA, Comai L, Henikoff S: High-throughput discovery of rare human nucleotide polymorphisms by Ecotilling. Nucleic Acids Res. 2006, 34 (13): e99-10.1093/nar/gkl479.PubMed CentralView ArticlePubMedGoogle Scholar
- Kennedy CL, O’Bryan MK: N-ethyl-N-nitrosourea (ENU) mutagenesis and male fertility research. Hum Reprod Update. 2006, 12 (3): 293-301. 10.1093/humupd/dmk004.View ArticlePubMedGoogle Scholar
- Mullins MC, Hammerschmidt M, Haffter P, Nusslein-Volhard C: Large-scale mutagenesis in the zebrafish: in search of genes controlling development in a vertebrate. Curr Biol. 1994, 4 (3): 189-202. 10.1016/S0960-9822(00)00048-8.View ArticlePubMedGoogle Scholar
- Imai Y, Feldman B, Schier AF, Talbot WS: Analysis of chromosomal rearrangements induced by postmeiotic mutagenesis with ethylnitrosourea in zebrafish. Genetics. 2000, 155 (1): 261-272.PubMed CentralPubMedGoogle Scholar
- Knapik EW: ENU mutagenesis in zebrafish-from genes to complex diseases. Mamm Genome. 2000, 11 (7): 511-519. 10.1007/s003350010098.View ArticlePubMedGoogle Scholar
- Jiang XY, Sun CF, Zhang QG, Zou SM: ENU-induced mutagenesis in grass carp (Ctenopharyngodon idellus) by treating mature sperm. PLoS One. 2011, 6 (10): e26475-10.1371/journal.pone.0026475.PubMed CentralView ArticlePubMedGoogle Scholar
- Riley BB, Grunwald DJ: Efficient induction of point mutations allowing recovery of specific locus mutations in zebrafish. Proc Natl Acad Sci U S A . 1995, 92 (13): 5997-6001. 10.1073/pnas.92.13.5997.PubMed CentralView ArticlePubMedGoogle Scholar
- Solnica-Krezel L, Schier AF, Driever W: Efficient recovery of ENU-induced mutations from the zebrafish germline. Genetics. 1994, 136 (4): 1401-1420.PubMed CentralPubMedGoogle Scholar
- Colbert T, Till BJ, Tompa R, Reynolds S, Steine MN, Yeung AT, McCallum CM, Comai L, Henikoff S: High-throughput screening for induced point mutations. Plant Physiol. 2001, 126 (2): 480-484. 10.1104/pp.126.2.480.PubMed CentralView ArticlePubMedGoogle Scholar
- Moens CB, Donn TM, Wolf-Saxon ER, Ma TP: Reverse genetics in zebrafish by TILLING. Brief Funct Genomic Proteomic. 2008, 7 (6): 454-459. 10.1093/bfgp/eln046.PubMed CentralView ArticlePubMedGoogle Scholar
- Gady AL, Hermans FW, Van de Wal MH, van Loo EN, Visser RG, Bachem CW: Implementation of two high through-put techniques in a novel application: detecting point mutations in large EMS mutated plant populations. Plant Methods. 2009, 5: 13-10.1186/1746-4811-5-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Ishikawa T, Kamei Y, Otozai S, Kim J, Sato A, Kuwahara Y, Tanaka M, Deguchi T, Inohara H, Tsujimura T, et al: High-resolution melting curve analysis for rapid detection of mutations in a Medaka TILLING library. BMC Mol Biol. 2010, 11: 70-10.1186/1471-2199-11-70.PubMed CentralView ArticlePubMedGoogle Scholar
- Li F, Niu B, Huang Y, Meng Z: Application of high-resolution DNA melting for genotyping in lepidopteran non-model species: Ostrinia furnacalis (Crambidae). PLoS One. 2012, 7 (1): e29664-10.1371/journal.pone.0029664.PubMed CentralView ArticlePubMedGoogle Scholar
- McPherron AC, Lawler AM, Lee SJ: Regulation of skeletal muscle mass in mice by a new TGF-beta superfamily member. Nature. 1997, 387 (6628): 83-90. 10.1038/387083a0.View ArticlePubMedGoogle Scholar
- Kambadur R, Sharma M, Smith TP, Bass JJ: Mutations in myostatin (GDF8) in double-muscled Belgian Blue and Piedmontese cattle. Genome Res. 1997, 7 (9): 910-916.PubMedGoogle Scholar
- McPherron AC, Lee SJ: Double muscling in cattle due to mutations in the myostatin gene. Proc Natl Acad Sci U S A . 1997, 94 (23): 12457-12461. 10.1073/pnas.94.23.12457.PubMed CentralView ArticlePubMedGoogle Scholar
- Mosher DS, Quignon P, Bustamante CD, Sutter NB, Mellersh CS, Parker HG, Ostrander EA: A mutation in the myostatin gene increases muscle mass and enhances racing performance in heterozygote dogs. PLoS Genet. 2007, 3 (5): e79-10.1371/journal.pgen.0030079.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuelke M, Wagner KR, Stolz LE, Hubner C, Riebel T, Komen W, Braun T, Tobin JF, Lee SJ: Myostatin mutation associated with gross muscle hypertrophy in a child. N Engl J Med. 2004, 350 (26): 2682-2688. 10.1056/NEJMoa040933.View ArticlePubMedGoogle Scholar
- Chisada S, Okamoto H, Taniguchi Y, Kimori Y, Toyoda A, Sakaki Y, Takeda S, Yoshiura Y: Myostatin-deficient medaka exhibit a double-muscling phenotype with hyperplasia and hypertrophy, which occur sequentially during post-hatch development. Dev Biol. 2011, 359 (1): 82-94. 10.1016/j.ydbio.2011.08.027.View ArticlePubMedGoogle Scholar
- Wansleeben C, van Gurp L, Feitsma H, Kroon C, Rieter E, Verberne M, Guryev V, Cuppen E, Meijlink F: An ENU-mutagenesis screen in the mouse: identification of novel developmental gene functions. PLoS One. 2011, 6 (4): e19357-10.1371/journal.pone.0019357.PubMed CentralView ArticlePubMedGoogle Scholar
- Russell WL, Kelly EM, Hunsicker PR, Bangham JW, Maddux SC, Phipps EL: Specific-locus test shows ethylnitrosourea to be the most potent mutagen in the mouse. Proc Natl Acad Sci U S A . 1979, 76 (11): 5818-5819. 10.1073/pnas.76.11.5818.PubMed CentralView ArticlePubMedGoogle Scholar
- Shibuya T, Horiya N, Matsuda H, Sakamoto K, Hara T: Dose-dependent induction of recessive mutations with N-ethyl-N-nitrosourea in primordial germ cells of male mice. Mutat Res. 1996, 357 (1–2): 219-224.View ArticlePubMedGoogle Scholar
- Berry C, Thomas M, Langley B, Sharma M, Kambadur R: Single cysteine to tyrosine transition inactivates the growth inhibitory function of piedmontese myostatin. Am J Physiol Cell Physiol. 2002, 283 (1): C135-141. 10.1152/ajpcell.00458.2001.View ArticlePubMedGoogle Scholar
- Henikoff S, Till BJ, Comai L: TILLING. Traditional mutagenesis meets functional genomics. Plant Physiol. 2004, 135 (2): 630-636. 10.1104/pp.104.041061.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee SJ, McPherron AC: Regulation of myostatin activity and muscle growth. Proc Natl Acad Sci U S A . 2001, 98 (16): 9306-9311. 10.1073/pnas.151270098.PubMed CentralView ArticlePubMedGoogle Scholar
- Bono C, Nuzzo D, Albeggiani G, Zizzo C, Francofonte D, Iemolo F, Sanzaro E, Duro G: Genetic screening of Fabry patients with EcoTILLING and HRM technology. BMC Res Notes. 2011, 4: 323-10.1186/1756-0500-4-323.PubMed CentralView ArticlePubMedGoogle Scholar
- Shima A, Shimada A: The Japanese medaka, Oryzias latipes, as a new model organism for studying environmental germ-cell mutagenesis. Environ Health Perspect. 1994, 102 (Suppl 12): 33-35. 10.1289/ehp.94102s1233.PubMed CentralView ArticlePubMedGoogle Scholar
- Wienholds E, Schulte-Merker S, Walderich B, Plasterk RH: Target-selected inactivation of the zebrafish rag1 gene. Science. 2002, 297 (5578): 99-102. 10.1126/science.1071762.View ArticlePubMedGoogle Scholar
- Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, et al: The sequence of the human genome. Science. 2001, 291: 1304-10.1126/science.1058040.View ArticlePubMedGoogle Scholar
- Ohta H: Cryopreservation of teleost spermatozoa. Cryobiology and Cryotechnology. 2008, 54 (1): 29-32.Google Scholar
- Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25 (14): 2745-2751. 10.1093/nar/25.14.2745.PubMed CentralView ArticlePubMedGoogle Scholar
- Anai T, Yamada T, Hideshima R, Kinoshita T, Rahman SM, Takagi Y: Two high-oleic-acid soybean mutants, M23 and KK21, have disrupted microsomal omega-6 fatty acid desaturase, encoded by GmFAD2-1a. Breed Sci. 2008, 58: 447-452. 10.1270/jsbbs.58.447.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.