High-resolution melting curve analysis for rapid detection of mutations in a Medaka TILLING library
- Tomoko Ishikawa†1,
- Yasuhiro Kamei†1,
- Shinji Otozai2,
- Jinhyong Kim1,
- Ayuko Sato3,
- Yoshikazu Kuwahara4,
- Minoru Tanaka5,
- Tomonori Deguchi6,
- Hidenori Inohara2,
- Tohru Tsujimura3 and
- Takeshi Todo1Email author
© Ishikawa et al; licensee BioMed Central Ltd. 2010
Received: 16 June 2010
Accepted: 15 September 2010
Published: 15 September 2010
During the last two decades, DNA sequencing has led to the identification of numerous genes in key species; however, in most cases, their functions are still unknown. In this situation, reverse genetics is the most suitable method to assign function to a gene. TILLING (Targeting Induced Local Lesions IN Genomes) is a reverse-genetic strategy that combines random chemical mutagenesis with high-throughput discovery of the induced mutations in target genes. The method has been applied to a variety of plant and animal species. Screening of the induced mutations is the most important step in TILLING. Currently, direct sequencing or nuclease-mediated screening of heteroduplexes is widely used for detection of mutations in TILLING. Both methods are useful, but the costs are substantial and turnaround times are relatively long. Thus, there is a need for an alternative method that is of higher throughput and more cost effective.
In this study, we developed a high resolution melting (HRM) assay and evaluated its effectiveness for screening ENU-induced mutations in a medaka TILLING library. We had previously screened mutations in the p53 gene by direct sequencing. Therefore, we first tested the efficiency of the HRM assay by screening mutations in p53, which indicated that the HRM assay is as useful as direct sequencing. Next, we screened mutations in the atr and atm genes with the HRM assay. Nonsense mutations were identified in each gene, and the phenotypes of these nonsense mutants confirmed their loss-of-function nature.
These results demonstrate that the HRM assay is useful for screening mutations in TILLING. Furthermore, the phenotype of the obtained mutants indicates that medaka is an excellent animal model for investigating genome stability and gene function, especially when combined with TILLING.
Our understanding of the basic mechanisms underlying most biological processes has been transformed by the systematic application of mutational analysis. Traditionally, forward genetics, driven by the identification of mutant phenotypes, has been the most widely used approach. On the other hand, genome sequencing projects over the past few decades have identified numerous genes in key species, and the completion of these sequences produced a situation in which most of the genes are known, but most of their phenotypes are obscure. In this situation, reverse genetics, which provides targeted inactivation of genes identified by sequence analysis followed by phenotype analysis of the mutant, has become an important tool for many biologists. In mice, reverse genetics is usually carried out using homologous recombination in embryonic stem cells, which allow a precise mutation to be constructed in nearly any gene. However, embryonic stem cells are only available in a limited number of organisms. Thus, a general method that is applicable to many organisms would be in great demand, and several approaches have been tried.
One of these approaches is TILLING (Targeting Induced Local Lesions IN Genomes). TILLING is a reverse-genetic strategy that combines random chemical mutagenesis with high-throughput discovery of the induced mutations in target genes. The method is general and, following its original application to the model plant Arabidopsis thaliana[1, 2], has been applied to a variety of plant and animal species including maize, lotus, barley, wheat, Drosophila, zebrafish, and medaka [3–9]. The first step in TILLING is chemical mutagenesis. For mutagenesis in animals, males are mutagenized using N-ethyl-N-nitrosourea (ENU) and then used to generate a large population of F1 animals that consequently harbor many random heterozygous mutations in their genomes. Next, the DNA from these animals is analyzed for mutations in a specific gene of interest. Once a mutation is identified, homozygous mutant animals can be obtained by crossing progeny from heterozygous F1 matings.
Among the vertebrates, small laboratory fish are suitable for the study of gene function due to their ease of handling, large numbers of progeny per generation, and, in particular, their translucent embryos. In many species, embryos develop outside the mother's body, enabling easy visual inspection and manipulation of their tissues and cells. One such fish is the zebrafish, Danio rerio, which is the most widely used laboratory fish. The success of forward genetics in the past two decades has established the zebrafish as the premier vertebrate for the study of gene function. Medaka, Oryzias latipes, is another small laboratory fish that has been used as an experimental model animal since the 1920s. Medaka has a small genome size (one-half of that of zebrafish) and is phylogenetically distinct from zebrafish, having diverged about 110 million years ago , making it useful for comparisons of conserved and divergent gene function in teleost evolution. For these reasons, both the zebrafish and medaka are currently widely used in comparative mutagenesis, carcinogenesis, and genomic instability studies. Shimada and Shima  showed that the mutation response of the medaka male germ cell is comparable to that of the mouse, and therefore medaka could serve as a vertebrate model system of DNA damage response in germ cells.
All cells have elaborate mechanisms to protect their genomes. DNA can be damaged by reactive metabolic byproducts and by environmental mutagens, and the response to DNA damage and its repair are crucial for cell viability and disease prevention. The DNA-damage response is a signal-transduction pathway that coordinates cell-cycle transitions, DNA replication, DNA repair, and apoptosis [12, 13]. The major regulators of the DNA-damage response are two protein kinases, ataxia-telangiectasia mutated (ATM) and ATM and Rad3-related (ATR) [14–17]. These two large kinases have significant sequence homology and play important roles in the initiation of cell cycle checkpoints by targeting an overlapping set of substrates that promote cell-cycle arrest and DNA repair . Another key player that mediates cell cycle checkpoints is the tumor suppressor p53, which is activated and stabilized by ATM and serves as a transcription factor for the downstream cell cycle regulator . The p53 pathway depends on the transcription of its downstream genes and its activity is therefore delayed after DNA damage. Thus, its role lies in the maintenance of the cell cycle checkpoint.
The most widely used reverse genetics' tool in zebrafish and medaka is undoubtedly morpholinos, but they are not used as a substitute for mutations because they can only be used in transient methods and for early developmental stages. Recently, TILLING was applied to zebrafish, and several mutants have been described . The success of this reverse genetics approach in zebrafish allowed us to apply TILLING to medaka. Once males have been mutagenizaed with ENU, the TILLING procedure can be divided into two steps, construction of a library and screening of the library. The library consists of pair sets of genomic DNA and frozen sperm of F1 offspring from each ENU-mutagenized male. We established a medaka TILLING library from 5771 F1 male fish derived from ENU-mutagenized G0 [9, 20].
In the next screening step, two methods are widely used for the detection of mutations, direct sequencing  and Cel I nuclease-mediated screening for heteroduplex formation with wild-type and mutant alleles . Both methods are effective; however, costs are substantial and turnaround time is relatively long. These disadvantages have resulted in the development of alternative methods that are more cost effective, faster, and easier to perform. The high resolution melting (HRM) analysis is a recently developed, promising technology used for the detection of variations in DNA [21, 22]. It is an in-tube method that can be performed in a fast, cheap, and robust manner. 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, which can be identified with HRM analysis. During HRM analysis, melting curves are produced using intercalating DNA dyes that fluoresce in the presence of double-stranded DNA and a specialized instrument designed to monitor fluorescence during heating [23, 24]. When the temperature increases, the DNA-intercalating dye is released from the DNA and the fluorescence decreases. This process produces a characteristic melting profile that can be monitored with precision. Changes in the sequence within the DNA fragment, as in SNPs or mutations, alter the melting profile. The HRM assay has successfully been used to detect point mutations in ethyl methanesulfonate-mutated plant populations .
In the present study, we developed HRM assays to evaluate the effectiveness of this methodology for screening ENU-induced mutations in the TILLING library. We first used the HRM assay to screen mutations in exons 5 and 6 of the p53 gene because the same library was previously screened with the direct sequencing method , enabling comparison of the two methods. We detected mutations in exons 5 and 6 of p53 with efficiency equivalent to that of the direct sequencing method. Then we screened for mutations in the atr and atm genes. We identified 15 mutations in atr and 11 in atm, including one nonsense mutation in each gene. The phenotypes of these nonsense mutants confirmed their loss-of-function nature. The results demonstrate the utility of the HRM assay for mutant screening using TILLING.
Screening of mutations in p53 gene
Mutations identified in the present study.
Mutation rate (Base-pairs/One mutation)
p53 mutations detected by sequencing and HRM.
Type of mutation or amino acid change*
5'-TGGCCCAGTA(T > A)TTTGAAGACC-3'
5'-CCGCAGGTAA(A > T)GCTGCCCCAG-3'
5'-CTACATGTGT(A > G)ACAGCTCGTG-3'
5'-TACATGTGTA(A > G)CAGCTCGTGC-3'
5'-GTGTAACAGC(T > C)CGTGCATGGG-3'
5'-GTGTAACAGCT(C > A)GTGCATGGG-3'
5'-TCTGGAAACC(G > T)AGTAAGTTTA-3'
5'-GGAAACCGAG(T > C)AAGTTTAGTC-3'
Screening of mutations in atr and atm genes
Summary of ATM and ATR mutations detected by HRM.
Amino acid change
5'-TGAAGCAGGA(T > A)GAGTCTGGCT-3'
5'-TGCAGAGGGC(T > G)GGAACATGAT-3'
5'-GGCTGGAACA(T > A)GATAAAGCTA-3'
5'-AGCTACTGAG(A > G)GAAAGGTAAA-3'
5'-TGAAAACTTC(A > G)GACGGCTCCT-3'
5'-AACTTCAGAC(G > A)GCTCCTTCGA-3'
5'-ATGGGACAAT(G > A)TCTTCAGGTA-3'
5'-CCACTGCTGT(C > A)GGCGTTGCAC-3'
5'-CTGTACGTGC(T > A)GCGCTGCCTC-3'
5'-CGTGCTGCGC(T > A)GCCTCAGGGA-3'
5'-CCACACAGAG(T > C)CTCTCAGCCT-3'
5'-ACCATGCCCA(T > C)CGGGGAGTTT-3'
5'-ACCTCAGGAT(T > A)GGACCAGCCT-3'
5'-TTGCAAAAAC(T > C)TCAGACCGGT-3'
5'-TGCATGGAAC(G > A)ATTCCTAGAC-3'
5'-CATGTGAACC(T > A)GACCAGGGAA-3'
5'-CAGGCGGTTC(T > C)CCGTTCTCTG-3'
5'-GCGGTTCTCC(G > A)TTCTCTGACC-3'
5'-CCTGGTCGCG(C > A)GTCTGAAGGA-3'
5'-GGCTTTGACG(A > G)GTGTGATGGC-3'
5'-GTGATGGCTC(T > G)GATGCGCCTG-3'
5'-CGGCCTCCGA(T > C)ACAAAGAAGA-3'
5'-AAAGAAGACT(T > A)CCCTCTGCTG-3'
5'-AGCCTGGTGT(C > A)GGTGCTGCTG-3'
5'-CGGCTCCCTT(T > A)GAGGAGAAAC-3'
5'-CATCCGCCGG(T > C)CTTCCATGAG-3'
atr is required for cellular proliferation
These viable fry were used to verify the knockout phenotype. We could not confirm the knockout phenotype at the protein level in the atr S1339X/S1339X , due to the lack of an appropriate anti-ATR antibody that recognizes the medaka protein. However, we were able to show the knockout phenotype by quantitative PCR. Quantitative PCR was carried out using two sets of PCR primers; each was designed to amplify the 5' or 3' region of the atr transcript. Both primer sets yielded the same results. The atr transcript levels decreased to ~70% of the wild-type for heterozygous mutants and ~40% for homozygous mutants, presumably due to nonsense-mediated decay of the transcript containing the premature stop codon (Fig. 4C). Finally, we confirmed the nonsense mutation in the transcript by sequencing the cDNA from the homozygous mutant (Fig. 4B).
The atmS444X mutant is sensitive to radiation, but partially infertile
The critical step in the TILLING procedure is the screening for mutations, but to successfully identify mutations in the target genes, it is necessary to have an efficient mutation discovery system. The TILLING library is constructed from F1s derived from chemical mutagen-treated G0s, and thus the mutations induced in each F1 are heterozygous; therefore, methods that quickly distinguish heteroduplex DNA from the homoduplex counterpart can be applied to the discovery of these mutations. HRM was recently introduced as a screening method for mutation detection and reported to be useful for heteroduplex detection [21, 22]. It is an in-tube method that can be performed in a fast, cheap, and robust manner. To evaluate the HRM assay as an alternative mutation screening method for TILLING, in this report we used the HRM assay to screen a total of 17.3 Mb in a medaka TILLING library. We identified 55 independent mutations (Table 1) and found that the average mutation frequency of the library was 1 mutation per 314 kb. The value is slightly higher than that obtained by direct sequencing the same library; 1 mutation per 345 kb . The higher rate of mutation detection by the HRM assay relative to that by the direct sequencing method is unexpected because direct sequencing is the gold standard for identifying mutations and thus would be expected to detect mutations more efficiently than any other method. Misclassification of SNPs by the HRM assay might be the reason for the discrepancy. The genomic DNA library used contains a substantial level of SNPs. When the same base change was identified more than twice, it was classified not as a mutation but as a SNP (see Materials and Methods). In the HRM assay, screening was carried out using 96-well PCR plates. Melting curve analysis was performed on each plate, but comparison between plates could not be performed. If one type of SNP appeared only once in one plate and appeared more than twice in other plates, the former was misclassified as a mutation and confirmed by sequencing, whereas the latter were correctly classified as SNPs and not subjected to checking by sequencing. However, misclassification could not have been common because the mutation detection rate was only slightly higher in the HRM assay than in the direct sequencing method. Despite this issue of classification, our results demonstrate the utility of the HRM assay for the screening of mutations in the TILLING library.
In this report, we used the HRM assay to identify medaka mutants producing truncated forms of the ATM or ATR proteins, atmS444X and atrS1339X, respectively. Atm and Atr orchestrate overlapping DNA damage responses to different forms of DNA strand discontinuities. However, knockout mouse models have also suggested that Atr, in contrast to Atm, is essential for viability, and that Atr-null mice die at an early embryonic stage [27, 28]. Interestingly, atr S1339X/S1339X medaka can continue normal development during embryogenesis and hatch successfully, but the hatched fry all die before day 30 postfertilization. We tried to establish a culture cell line from atr S1339X/S1339X embryos, but the cells died within a few weeks and no cell line was obtained (data not shown), a finding similar to results with cells from Atr-null mice [27, 28]. It is not clear at present why atr-deficient embryos can survive for several weeks; however, the viable atr fry provide a good model system for analyzing atr function at the organismal level.
Meiosis is a specialized cell division program to generate hemizygous gametes. An important feature of meiosis is the exchange of genetic information between maternal and paternal chromatids, a process called meiotic recombination [33, 34]. Meiotic recombination is initiated by Spo11-catalyzed meiotic DNA double-strand breaks (DSBs), followed by their repair using an intact homologous non-sister chromatid . DSBs are highly hazardous for genome integrity, but meiotic cells deliberately introduce them into their genome in order to initiate homologous recombination. Therefore, DSBs left unrepaired are checked by a stringent surveillance mechanism, the meiotic recombination checkpoint, which delays meiosis I until DSB repair is completed. Another important genome surveillance mechanism is the DNA damage checkpoint, which senses DSBs that arise at unpredictable locations as a consequence of DNA damage during both mitosis and meiosis [13, 33, 36, 37]. Mechanistically, the meiotic recombination checkpoint is related to the DNA damage checkpoint; however, many aspects of the meiotic recombination checkpoint are considerably less well understood than those of mitotic DNA damage checkpoint. Given the well established role of ATM in DSB repair in mitosis, ATM might play an important role in the signaling and processing of the programmed meiotic DSBs. In this study, we found defective meiosis in the medaka atm S444X mutant. The striking difference between Atm mice and the mutant fish is the fertility of the latter; the Atm mice are completely infertile, whereas the mutant fish are partially infertile, and 25% of females and 57% of males remain fertile. In mice, Atm deficiency results in early arrest in prophase I followed by apoptotic degeneration, which results in the complete absence of mature gametes and spermatocytes in various stages of degeneration [29, 30, 38]. Similar premeiotic arrest was observed in tissue sections of the infertile mutant fish, but the spermatocytes did not show any stage of degeneration (Fig. 7B). Instead, a high percentage of abnormally large spermatids and/or sperm were noted in the testes of the fertile mutants. These observations suggest that, in medaka, many cells bearing unrepaired meiotic DSBs are removed by apoptosis, but a substantial fraction proceed to completion of meiosis, resulting in the production of abnormal sperm. In fact, a similar partially sterile phenotype has been reported in Arabidopsis ATM mutants , whose ATM meiocytes did not appear to undergo apoptosis before the end of meiosis despite the observed chromosomal fragmentation during anaphase I and II. The histology presented in the current report suggests that the infertility in atm-deficient medaka is possibly due to disorganized gonads in both males and females. At present, the precise role of the meiotic recombination checkpoint defect in atm medaka in causing the reproductive defects is not clear. More precise characterization of the abnormal sperm and oocytes should provide important information about the meiotic recombination checkpoint in atm mutant fish.
Finally, we should emphasize again the possible involvement of undetected ENU-induced mutations in the observed phenotype. The single major drawback of TILLING as a method for reverse genetics is that a mutation of interest that is identified in any given F1 fish is only one of many heterozygous mutations in its genome. By multiplying the estimated size of medaka genome by the mutation rate detected in the present TILLING genomic DNA library, an estimate of the number of ENU-induced mutations per F1 fish is ~2300. More than 80% of these mutations were removed after backcrossing to wild-type fish three times, however significant levels of mutations still remain. As a result, we cannot exclude possibility that the remaining mutations affect the phenotype. Most strong evidence for the direct relationship between the observed phenotype and the mutated gene is the phenotype rescue with the wild-type gene. Further studies including the phenotype rescue with BAC (Bacterial artificial chromosome) clone, will clarify a possible involvement of these unknown mutations.
ATM contains three domains that are important for its function, a FAT domain, a protein kinase domain, and a FATC domain, which are positioned in the C-terminal half of the protein. In addition to these domains, several known ATM substrates bind to the substrate binding site, a region near the N-terminus of the protein . The N-terminus is a crucial region of the ATM because its deletion inactivates the protein [40, 41]. In this study, we demonstrated the presence of a short atm transcript variant in medaka cells that encodes an isoform (isoform 2) lacking the N-terminal 767 residues of the full-length ATM protein. A similar short isoform has also been reported in human cells (GenBank Accession No. NM 138242 or Ensemble Transcript: ENST00000452508). These short isoforms lack the substrate binding site and thus cannot complement the radiosensitivity of atm-deficient mutants. Mutations in ataxia telangiectasia patients are well characterized: over 200 different ATM mutations have been described, and approximately 85% of the mutants reported in ataxia telangiectasia patients are predicted to truncate the protein, but the mutations are mapped uniformly over the gene, which includes the N-terminal region truncated in the short isoform [41, 42]. Therefore, at present we do not know the function of the short isoform, and there is no answer to why the short isoform is conserved in both human and fish. However, the medaka atm S444X/S444X mutant provides an excellent model for analyzing the function of this short atm isoform.
In the present study we demonstrate that the HRM assay is useful for screening mutations in TILLING. The present study also demonstrates that mutations in medaka atr affect viability and in atm affect fertility and radiosensitivity, similar to findings in mice and humans. However, the phenotypes of these mutants show milder effects in medaka than in mammals. On the basis of our results and those of others, we conclude that medaka is a useful model for the analysis of genome stability, and can be used to complement studies in zebrafish, mouse and human.
Medaka TILLING Library
Mutagenesis was carried out as described previously [9, 43]. Briefly, 102 male Kyoto-Cab, a substrain of Cab, were mutagenized by 3 consecutive treatments every week with 3 mM ENU (Sigma-Aldrich, St. Louis, MO). One month after the last ENU treatment, 97 surviving fertile fish were crossed with wild-type females to obtain the mutagenized F1 male library. The number of offspring produced from a single mutagenized male founder varied from 1 to 239. The sperm from each F1 medaka was cryopreserved, as described previously (Medaka Book, http://www.shigen.nig.ac.jp/medaka/medakabook/). After removal of the testes, the fish was kept at -80°C until DNA was extracted. The posterior half of the fish was incubated overnight at 55°C in lysis buffer containing 10 mM Tris-HCl (pH 7.5), 1 mM EDTA (pH 8.0), 120 mM sodium chloride, 12 mM sodium citrate (pH 7.0), 1% SDS, and 200 μg/mL proteinase K. The lysate was extracted with phenol and precipitated with isopropanol. The DNA pellet was dissolved in 1 mL TE (10 mM Tris pH 7.5, 1 mM EDTA, pH 8.0). The concentration was adjusted to 4 ng/μL, aliquoted into 96-deep-well plates, and stored at -20°C. Finally, sets of frozen sperm and genomic DNA were prepared from 5771 F1 males and used as the TILLING library.
Identification of mutations by HRM analysis
BLAST searches of medaka databases (http://www.ensembl.org/Oryzias_latipes/index.html) using the human ATR or ATM gene revealed one homologous gene each (Olatr or Olatm, respectively) in the medaka genome and no paralogous genes. The degrees of conservation with the human ATR or ATM are 59.4% or 51.9% identity, respectively. Each genomic DNA in the library was amplified by PCR with gene-specific primers. To design the gene-specific primers, parts of the cDNA sequences of atm and atr were determined by RT-PCR and used to retrieve the genomic sequences from two Medaka databases (http://www.ensembl.org/Oryzias_latipes/index.html and http://medaka.utgenome.org/). The following 3 criteria were used for selection of the target exon for each gene: 1) location in the N-terminal region, 2) correspondence to the region encoding conserved amino acid residues, and 3) sufficient length (more than 150 base pairs). The first two criteria were employed to obtain null-function mutants. The third criterion was selected for high cost performance. We designed the primer set using Light Scanner Primer Design software (Idaho Technology, Salt Lake City, UT, USA). The designed primers were further checked for simulated melting curves using MeltSim 1.0-beta for Windows (http://bioinformatics.org/meltsim/wiki/). The primers used are summarized in Additional file 2, Table S1.
PCR reactions were performed in 96-well microtiter plates in 7-μL volumes. Reactions included 7 ng of genomic DNA in 1 × KOD Plus PCR buffer (Toyobo, Tokyo), with 1 mM MgSO4, 200 μM of deoxynucleotide triphosphate, 0.6 × LC Green PLUS (Idaho Technology), 0.14 U of KOD Plus polymerase (Toyobo, Tokyo), and 300 nM of each primer. DMSO (10% vol/vol) was added when the amplified region was GC-rich. Reactions were overlaid with 10 μL of mineral oil (Nacalai, Kyoto, Japan), the plates were centrifuged (1500 × g for 3 min), and PCR was performed in an iCycler (Bio-Rad Laboratories, Hercules, CA) or PXE0.2 Thermal Cycler (Thermo Scientific, Rockford, IL). The cycling and melting conditions were set as follows: one cycle of 94°C for 2 minutes, followed by 45 cycles of 94°C for 15 seconds; a primer-specific annealing temperature for 25 seconds, and then 68°C for 25 seconds; and a final denaturing and re-annealing step (one cycle of 94°C for 30 seconds followed by rapid cooling to 28°C. The final step was included to redistribute DNA strands derived from mutant and wild-type alleles in order to maximize heteroduplex formation.
After PCR, the plates were centrifuged (1500 × g for 3 min) and imaged in a 96-well LightScanner (Idaho Technology). The plates were heated at 0.3°C/second, and fluorescence measurements were collected from 65°C to 98°C. Melting curves were analyzed as described previously  using the LightScanner software. Briefly, after exponential background subtraction, melting curves were normalized between 0% and 100%. Normalized and temperature-overlaid curves were viewed on the subtraction plots to magnify differences in the shapes of the melting curves. The subtraction plots were generated by subtracting each curve from the mean wild-type curve, defined as the most common genotype. The subtraction plot helps to cluster the samples into groups. Clustering of the melting curves for genotype identification was performed manually using the LightScanner software on high sensitivity.
The fluorescence data generated during DNA melting can be analyzed based on the melting temperature (Tm) or on the shape of the melting curve. Most homozygous sequence changes produce a Tm shift compared to the wild type . In contrast, heterozygous samples are identified not by product Tm, but by differences in melting curve configuration . In the TILLING library, all the induced mutations were heterozygous. PCR amplification of heterozygotes followed by heat denaturing and annealing resulted 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 observed by melting curve analysis. The subtraction plot was used to detect changes in the contour of the melting curve.
HRM screening of the TILLING library was carried out in two steps. First, each of the 5771 genomic DNAs was assayed by HRM, and the genomic DNAs showing melting curves different from the wild-type allele were selected as mutant candidates (first positive). The HRM assay was repeated for the first positive genomic DNAs, and the genomic DNAs showing melting curves different from the wild-type allele were picked up again as second positives. Proof of induced mutation was obtained by sequencing the PCR product of the second positive genomic DNAs. As reported previously , substantial numbers of SNPs exist in the genomic DNA of the library used. Therefore, the second positive DNAs included both SNPs and true mutations. SNPs were discriminated from true mutations by defining the same type of base change detected in more than two sequencing reactions as an SNP.
PCR products containing the second positive HRM amplicons were purified using a Sephadex G50 (Fine DNA Grade) column and then used as a template for the sequencing reaction. Sequencing reactions were carried out by using BigDye Terminator version 3.1 (Applied Biosystems, Foster City, CA) and the ABI 3730xl sequencing platform.
Generation of atm and atr mutant fish
F2 fish carrying the atmS444X or atrS1339X mutation were obtained by artificial insemination using the frozen sperm from the TILLING library. Artificial insemination was performed as described previously . Briefly, about 100 unfertilized eggs were obtained from the wild-type females. A single glass capillary containing 10 μL of sperm from the F1 fish identified in the screening was removed from the liquid nitrogen and thawed at ambient temperature. Immediately after thawing, the content was placed in balanced salt solution (BSS; 0.65% sodium chloride, 0.04% potassium chloride, 0.02% magnesium sulfate heptahydrate, 0.02% calcium chloride dihydrate, 0.00005% phenol red, 0.01% sodium hydrogen carbonate, pH 7.3) and incubated with eggs for 20 minutes with occasional pipetting. Eggs that were not fertilized were removed 3 h later, and BSS was replaced with 0.03% Red Sea salt water. The eggs were incubated at 28°C until hatching. Artificial insemination was conducted using atmS444X or atrS1339X sperm with a Cab or HdrR female, respectively. Each resultant fish was crossed out twice to wild-type Cab or HdrR fish, respectively. Three independent in-crosses of heterozygous progeny were then performed, from which homozygous mutant males were obtained and used for this study.
Fish were genotyped by PCR analysis or sequencing of fin-clip DNA. For PCR genotyping, PCR reactions were performed using two primer sets, wild-type and mutant, which have the wild-type or mutant allele sequence, respectively, at the 3' end. The primers used for genotyping are listed in Additional file 4, Table S2. For genotyping by sequencing, PCR reactions were performed with the same primers as used for the screening of the library.
Primary cell culture
Cell cultures were established from embryos as described previously (http://www.shigen.nig.ac.jp/medaka/medakabook/). Cells were cultured at 27°C. Briefly the embryos were sterilized for at least 40 s in Dakin's solution (58 mM NaOCl, 6 mM HCl, and 93 mM NaHCO3), washed in phosphate-buffered saline (PBS), resterilized for 10 s in 70% ethanol, and washed in PBS. The chorion and yolk sac were removed from each embryo, transferred to L-15 medium, and minced using a 22G needle connected to a 1-mL syringe to obtain fibroblast-like cells. Cells were cultured in L-15 medium supplemented with 10% fetal bovine serum (Gibco BRL, Grand Island, NY), 50 mg/mL streptomycin, 50 U/mL penicillin, and 10 mM 2-[4-(-hydroxyethyl)-1-piperazinyl] ethanesulfonic acid (HEPES; pH 7.5), and incubated at 27°C.
Sensitivity to X-ray irradiation
Cultured cells were exposed to γ-rays from a 137Cs source at a dose rate of 0.93 Gy/min (Gamma Cell 6000 Elan, MDS Nordion, Ottawa, Canada). OLCAB-e3 cells were irradiated with γ-rays at the dose of 100 Gy and used for feeder cells. After irradiation, atm +/+ or atm S444X/S444X cells were cultured in 6-well plates in the presence of feeder cells (7 × 105 cells/well) and incubated at 27°C. Two weeks after irradiation, the cells were fixed with ethanol and stained with 0.006% crystal violet, the formed colonies were counted, and the survival rate was calculated.
Survival of ATR-/- fish
atr +/S1339X fishes were intercrossed and eggs were cultured at 27°C in a 10-cm plastic dish. After culture for 7, 10, 15, 20, 25, or 30 days, the genotypes of surviving eggs or fry (80-100 for each fraction) were determined by PCR or by sequencing.
The gonads were dissected and fixed in Bouin's solution. Fixed tissues were embedded in paraffin blocks. Cross-sections were cut at a 5-μm thickness and stained with hematoxylin-eosin. Sections were examined and photographed under light microscopy.
Reverse-transcription PCR and quantitative PCR
RNA was isolated from cultured cells established from mutant or wild-type fish using Sepasol RNAI Super (Nacalai Tesque) and reverse transcribed by using the ReverTra AceR qPCR kit (Toyobo). For atr and atm quantitative PCR, cDNA-specific primers were designed to amplify either the 5' or 3' region of the atm and atr transcripts (5' regions were within exons 1-3, 3-4, 7-8 for atr and exons 2-3 for atm; 3' regions were within exons 47-48 for atr and exons 33-34 for atm). The sequences of all oligonucleotides that were used are listed in Additional file 4, Table S2.
Contribution of Authors
TI and YKa participated in the design and performance of all experiments described in this paper. In addition, TI wrote the first draft of the first manuscript. SO, JK, and HI participated in mutant screening with the HRM assay. AS and TTs performed histology, and TD participated in identification of the atr phenotype. YKu and MT provided expertise on histological examinations of fish gonads and supervised their interpretation. HI and TTo secured funding for the experiments. TTo conceived the project; designed the experimental plan; and supervised the design, execution, and interpretation of all experiments.
We thank Yasuko Morita, Kumi Nakamura, and Junko Ishikawa for the excellent technical assistance. This work was supported in part by Grants-in-Aid for Scientific Research (16201011 and 19101002 to TT, and 18681007 to TI).
- McCallum CM, Comai L, Greene EA, Henikoff S: Targeting induced local lesions IN genomes (TILLING) for plant functional genomics. Plant Physiol. 2000, 123 (2): 439-442. 10.1104/pp.123.2.439View ArticlePubMedPubMed CentralGoogle 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.480View ArticlePubMedPubMed CentralGoogle 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-12View ArticlePubMedPubMed CentralGoogle Scholar
- Perry JA, Wang TL, Welham TJ, Gardner S, Pike JM, Yoshida S, Parniske M: A TILLING reverse genetics tool and a web-accessible collection of mutants of the legume Lotus japonicus. Plant Physiol. 2003, 131 (3): 866-871. 10.1104/pp.102.017384View ArticlePubMedPubMed CentralGoogle Scholar
- Caldwell DG, McCallum N, Shaw P, Muehlbauer GJ, Marshall DF, Waugh R: A structured mutant population for forward and reverse genetics in Barley (Hordeum vulgare L.). Plant J. 2004, 40 (1): 143-150. 10.1111/j.1365-313X.2004.02190.xView 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/nbt1043View ArticlePubMedGoogle Scholar
- Winkler S, Schwabedissen A, Backasch D, Bokel C, Seidel C, Bonisch S, Furthauer M, Kuhrs A, Cobreros L, Brand M, et al: Target-selected mutant screen by TILLING in Drosophila. Genome Res. 2005, 15 (5): 718-723. 10.1101/gr.3721805View ArticlePubMedPubMed CentralGoogle Scholar
- Wienholds E, van Eeden F, Kosters M, Mudde J, Plasterk RH, Cuppen E: Efficient target-selected mutagenesis in zebrafish. Genome Res. 2003, 13 (12): 2700-2707. 10.1101/gr.1725103View ArticlePubMedPubMed CentralGoogle 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-r116View ArticlePubMedPubMed CentralGoogle Scholar
- Wittbrodt J, Shima A, Schartl M: Medaka--a model organism from the far East. Nat Rev Genet. 2002, 3 (1): 53-64. 10.1038/nrg704View ArticlePubMedGoogle 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 USA. 1991, 88 (6): 2545-2549. 10.1073/pnas.88.6.2545View ArticlePubMedPubMed CentralGoogle Scholar
- Hartwell LH, Weinert TA: Checkpoints: controls that ensure the order of cell cycle events. Science. 1989, 246 (4930): 629-634. 10.1126/science.2683079View ArticlePubMedGoogle Scholar
- Zhou BB, Elledge SJ: The DNA damage response: putting checkpoints in perspective. Nature. 2000, 408 (6811): 433-439. 10.1038/35044005View ArticlePubMedGoogle Scholar
- Lobrich M, Jeggo PA: The impact of a negligent G2/M checkpoint on genomic instability and cancer induction. Nat Rev Cancer. 2007, 7 (11): 861-869. 10.1038/nrc2248View ArticlePubMedGoogle Scholar
- Lavin MF: Ataxia-telangiectasia: from a rare disorder to a paradigm for cell signalling and cancer. Nat Rev Mol Cell Biol. 2008, 9 (10): 759-769. 10.1038/nrm2514View ArticlePubMedGoogle Scholar
- Cimprich KA, Cortez D: ATR: an essential regulator of genome integrity. Nat Rev Mol Cell Biol. 2008, 9 (8): 616-627. 10.1038/nrm2450View ArticlePubMedPubMed CentralGoogle Scholar
- Shiloh Y: The ATM-mediated DNA-damage response: taking shape. Trends Biochem Sci. 2006, 31 (7): 402-410. 10.1016/j.tibs.2006.05.004View ArticlePubMedGoogle Scholar
- Abraham RT: PI 3-kinase related kinases: 'big' players in stress-induced signaling pathways. DNA Repair (Amst). 2004, 3 (8-9): 883-887. 10.1016/j.dnarep.2004.04.002View ArticleGoogle Scholar
- Canman CE, Lim DS, Cimprich KA, Taya Y, Tamai K, Sakaguchi K, Appella E, Kastan MB, Siliciano JD: Activation of the ATM kinase by ionizing radiation and phosphorylation of p53. Science. 1998, 281 (5383): 1677-1679. 10.1126/science.281.5383.1677View ArticlePubMedGoogle Scholar
- Mitani H, Kamei Y, Fukamachi S, Oda S, Sasaki T, Asakawa S, Todo T, Shimizu N: The medaka genome: why we need multiple fish models in vertebrate functional genomics. Genome Dyn. 2006, 2: 165-182. full_textView ArticlePubMedGoogle Scholar
- Wittwer CT, Reed GH, Gundry CN, Vandersteen JG, Pryor RJ: High-resolution genotyping by amplicon melting analysis using LCGreen. Clin Chem. 2003, 49 (6 Pt 1): 853-860. 10.1373/49.6.853View ArticlePubMedGoogle Scholar
- Erali M, Voelkerding KV, Wittwer CT: High resolution melting applications for clinical laboratory medicine. Exp Mol Pathol. 2008, 85 (1): 50-58. 10.1016/j.yexmp.2008.03.012View ArticlePubMedPubMed CentralGoogle Scholar
- Herrmann MG, Durtschi JD, Bromley LK, Wittwer CT, Voelkerding KV: Amplicon DNA melting analysis for mutation scanning and genotyping: cross-platform comparison of instruments and dyes. Clin Chem. 2006, 52 (3): 494-503. 10.1373/clinchem.2005.063438View ArticlePubMedGoogle Scholar
- Herrmann MG, Durtschi JD, Wittwer CT, Voelkerding KV: Expanded instrument comparison of amplicon DNA melting analysis for mutation scanning and genotyping. Clin Chem. 2007, 53 (8): 1544-1548. 10.1373/clinchem.2007.088120View 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-13View ArticlePubMedPubMed CentralGoogle Scholar
- Reed GH, Wittwer CT: Sensitivity and specificity of single-nucleotide polymorphism scanning by high-resolution melting analysis. Clin Chem. 2004, 50 (10): 1748-1754. 10.1373/clinchem.2003.029751View ArticlePubMedGoogle Scholar
- Brown EJ, Baltimore D: ATR disruption leads to chromosomal fragmentation and early embryonic lethality. Genes Dev. 2000, 14 (4): 397-402.PubMedPubMed CentralGoogle Scholar
- de Klein A, Muijtjens M, van Os R, Verhoeven Y, Smit B, Carr AM, Lehmann AR, Hoeijmakers JH: Targeted disruption of the cell-cycle checkpoint gene ATR leads to early embryonic lethality in mice. Curr Biol. 2000, 10 (8): 479-482. 10.1016/S0960-9822(00)00447-4View ArticlePubMedGoogle Scholar
- Barlow C, Hirotsune S, Paylor R, Liyanage M, Eckhaus M, Collins F, Shiloh Y, Crawley JN, Ried T, Tagle D, et al: Atm-deficient mice: a paradigm of ataxia telangiectasia. Cell. 1996, 86 (1): 159-171. 10.1016/S0092-8674(00)80086-0View ArticlePubMedGoogle Scholar
- Xu Y, Ashley T, Brainerd EE, Bronson RT, Meyn MS, Baltimore D: Targeted disruption of ATM leads to growth retardation, chromosomal fragmentation during meiosis, immune defects, and thymic lymphoma. Genes Dev. 1996, 10 (19): 2411-2422. 10.1101/gad.10.19.2411View ArticlePubMedGoogle Scholar
- Elson A, Wang Y, Daugherty CJ, Morton CC, Zhou F, Campos-Torres J, Leder P: Pleiotropic defects in ataxia-telangiectasia protein-deficient mice. Proc Natl Acad Sci USA. 1996, 93 (23): 13084-13089. 10.1073/pnas.93.23.13084View ArticlePubMedPubMed CentralGoogle Scholar
- Saito D, Morinaga C, Aoki Y, Nakamura S, Mitani H, Furutani-Seiki M, Kondoh H, Tanaka M: Proliferation of germ cells during gonadal sex differentiation in medaka: Insights from germ cell-depleted mutant zenzai. Dev Biol. 2007, 310 (2): 280-290. 10.1016/j.ydbio.2007.07.039View ArticlePubMedGoogle Scholar
- Longhese MP, Bonetti D, Guerini I, Manfrini N, Clerici M: DNA double-strand breaks in meiosis: checking their formation, processing and repair. DNA Repair (Amst). 2009, 8 (9): 1127-1138. 10.1016/j.dnarep.2009.04.005View ArticleGoogle Scholar
- Petronczki M, Siomos MF, Nasmyth K: Un menage a quatre: the molecular biology of chromosome segregation in meiosis. Cell. 2003, 112 (4): 423-440. 10.1016/S0092-8674(03)00083-7View ArticlePubMedGoogle Scholar
- Keeney S, Giroux CN, Kleckner N: Meiosis-specific DNA double-strand breaks are catalyzed by Spo11, a member of a widely conserved protein family. Cell. 1997, 88 (3): 375-384. 10.1016/S0092-8674(00)81876-0View ArticlePubMedGoogle Scholar
- Terzoudi GI, Manola KN, Pantelias GE, Iliakis G: Checkpoint abrogation in G2 compromises repair of chromosomal breaks in ataxia telangiectasia cells. Cancer Res. 2005, 65 (24): 11292-11296. 10.1158/0008-5472.CAN-05-2148View ArticlePubMedGoogle Scholar
- Hartwell LH, Kastan MB: Cell cycle control and cancer. Science. 1994, 266 (5192): 1821-1828. 10.1126/science.7997877View ArticlePubMedGoogle Scholar
- Hamer G, Kal HB, Westphal CH, Ashley T, de Rooij DG: Ataxia telangiectasia mutated expression and activation in the testis. Biol Reprod. 2004, 70 (4): 1206-1212. 10.1095/biolreprod.103.024950View ArticlePubMedGoogle Scholar
- Garcia V, Bruchet H, Camescasse D, Granier F, Bouchez D, Tissier A: AtATM is essential for meiosis and the somatic response to DNA damage in plants. Plant Cell. 2003, 15 (1): 119-132. 10.1105/tpc.006577View ArticlePubMedPubMed CentralGoogle Scholar
- Young DB, Jonnalagadda J, Gatei M, Jans DA, Meyn S, Khanna KK: Identification of domains of ataxia-telangiectasia mutated required for nuclear localization and chromatin association. J Biol Chem. 2005, 280 (30): 27587-27594. 10.1074/jbc.M411689200View ArticlePubMedGoogle Scholar
- Lavin MF, Scott S, Gueven N, Kozlov S, Peng C, Chen P: Functional consequences of sequence alterations in the ATM gene. DNA Repair (Amst). 2004, 3 (8-9): 1197-1205. 10.1016/j.dnarep.2004.03.011View ArticleGoogle Scholar
- Meyn MS: Ataxia-telangiectasia, cancer and the pathobiology of the ATM gene. Clin Genet. 1999, 55 (5): 289-304. 10.1034/j.1399-0004.1999.550501.xView ArticlePubMedGoogle Scholar
- Furutani-Seiki M, Sasado T, Morinaga C, Suwa H, Niwa K, Yoda H, Deguchi T, Hirose Y, Yasuoka A, Henrich T, et al: A systematic genome-wide screen for mutations affecting organogenesis in Medaka, Oryzias latipes. Mech Dev. 2004, 121 (7-8): 647-658. 10.1016/j.mod.2004.04.016View ArticlePubMedGoogle Scholar
- Montgomery J, Wittwer CT, Palais R, Zhou L: Simultaneous mutation scanning and genotyping by high-resolution DNA melting analysis. Nat Protoc. 2007, 2 (1): 59-66. 10.1038/nprot.2007.10View ArticlePubMedGoogle Scholar
- Palais RA, Liew MA, Wittwer CT: Quantitative heteroduplex analysis for single nucleotide polymorphism genotyping. Anal Biochem. 2005, 346 (1): 167-175. 10.1016/j.ab.2005.08.010View ArticlePubMedGoogle Scholar
- Graham R, Liew M, Meadows C, Lyon E, Wittwer CT: Distinguishing different DNA heterozygotes by high-resolution melting. Clin Chem. 2005, 51 (7): 1295-1298. 10.1373/clinchem.2005.051516View ArticlePubMedGoogle Scholar
- Kamei Y, Itou J, Oda S, Masui M, Kim JH, Ishikawa T, Yuba S, Kinoshita M, Mitani H, Todo T: Development of a convenient in vitro fertilization method using interspecific hybrids between Oryzias latipes and Oryzias curvinotus. Dev Growth Differ. 2007, 49 (9): 721-730. 10.1111/j.1440-169X.2007.00966.xView ArticlePubMedGoogle Scholar