RIPCAL: a tool for alignment-based analysis of repeat-induced point mutations in fungal genomic sequences
© Hane and Oliver; licensee BioMed Central Ltd. 2008
Received: 28 April 2008
Accepted: 12 November 2008
Published: 12 November 2008
Repeat-induced point mutation (RIP) is a fungal-specific genome defence mechanism that alters the sequences of repetitive DNA, thereby inactivating coding genes. Repeated DNA sequences align between mating and meiosis and both sequences undergo C:G to T:A transitions. In most fungi these transitions preferentially affect CpA di-nucleotides thus altering the frequency of certain di-nucleotides in the affected sequences. The majority of previously published in silico analyses were limited to the comparison of ratios of pre- and post-RIP di-nucleotides in putatively RIP-affected sequences – so-called RIP indices. The analysis of RIP is significantly more informative when comparing sequence alignments of repeated sequences. There is, however, a dearth of bioinformatics tools available to the fungal research community for alignment-based RIP analysis of repeat families.
We present RIPCAL http://www.sourceforge.net/projects/ripcal, a software tool for the automated analysis of RIP in fungal genomic DNA repeats, which performs both RIP index and alignment-based analyses. We demonstrate the ability of RIPCAL to detect RIP within known RIP-affected sequences of Neurospora crassa and other fungi. We also predict and delineate the presence of RIP in the genome of Stagonospora nodorum – a Dothideomycete pathogen of wheat. We show that RIP has affected different members of the S. nodorum rDNA tandem repeat to different extents depending on their genomic contexts.
The RIPCAL alignment-based method has considerable advantages over RIP indices for the analysis of whole genomes. We demonstrate its application to the recently published genome assembly of S. nodorum.
Over 100 fungal genome sequences have been obtained or are in the pipeline  and next-generation sequencing technologies will further accelerate the accumulation of data over the next decade. This rapidly growing array of sequence information presents many new challenges for analysis. There is an urgent need to develop and implement efficient tools to describe features of new genomes. Repeat-induced point mutation (RIP) is one such area of fungal biology requiring efficient analytical tools. RIP is an irreversible genome defence mechanism first detected in Neurospora crassa [2, 3] and subsequently in Magnaporthe grisea [4, 5], Podospora anserina  and Leptosphaeria maculans . RIP is believed to be a defence against transposons, rendering them inactive and protecting sexual progeny from the expression of transposon genes.
Direct experimental observation of RIP requires both that the fungal species can be crossed under laboratory conditions and that the strain can be transformed with multiple copies of a transgene. Very few fungal species are amenable to such analysis and these procedures are slow in all cases. RIP-like processes can also be detected by in-silico analysis of repeated elements in whole or partial genomic sequences. Prior examples include Aspergillus fumigatus , Fusarium oxysporum [9–11], Aspergillus nidulans , Microbotryum violaceum , Magnaporthe oryzae , Aspergillus niger  and Penicillium chysogenum . We now have the opportunity to detect and measure RIP in silico from genomic sequences of diverse species.
The four possible CpN→TpN di-nucleotide RIP mutations and their reverse complements which form the basis for comparisons to determine the dominant form of RIP mutation in both alignment-based and statistical analyses.
RIP-indices are simple to calculate and do not require complete knowledge of the genome sequence or repeat families. They are also applicable to heavily mutated repeat families for which an alignment is not possible or questionable. However, RIP indices are insensitive tools which obscure many interesting features of RIP. These include the direction of RIP changes (i.e. which sequence is closer to the ancestral precursor of the RIP-affected sequence), the degree of RIP along the length of repeat alignments and differences in RIP profiles between members of the repeat class.
As RIP operates on aligned sequences, these questions are better answered using an alignment-based approach. Alignment-based analysis of RIP involves the multiple alignment of a repeat family and counting RIP mutations along the alignment for all sequences. This method has been previously used to identify RIP within the Ty1 transposon family of Microbotryum violaceum using the software tool Sequencher. Such manual calculation of RIP as was used by Hood et al  does not lend itself to whole genome RIP analysis. To enable a thorough, facile and automated analysis of RIP in the plethora of new fungal genomes, we have developed the free software tool RIPCAL (available at http://www.sourceforge.net/projects/ripcal. RIPCAL incorporates both RIP index and alignment-based methods. Its capabilities are demonstrated with examples taken from de novo-defined repeat families of the recently published Stagonospora nodorum genome, a major fungal pathogen of wheat [21, 22].
Validation of RIP detection by the alignment-based method
The RIPCAL alignment-based method was applied to both the 5S rDNA repeat family of Neurospora crassa, which is reportedly free from RIP mutation due to its short sequence length , and to the Tad1 transposons of N. crassa, which are reported to be heavily prone to CpA→TpA RIP mutation . The 5S rDNA and Tad1 repeat families served as negative and positive controls for RIP respectively. Analysis showed low levels of RIP mutation among 5S rDNAs, whereas high levels of RIP mutation were detected amongst Tad1 transposons as expected (Additional file 1). Interestingly, while CpA↔TpA changes were highly increased in the Tad1 family, these were overshadowed by a major increase in CpT↔TpT mutation, which has not been previously detected . This may be due to the fact that the former study compared Tad1 sequences between different strains of Neurospora crassa, whereas this comparison was restricted to all repeats within a single strain.
Identification of the dominant CpN to TpN di-nucleotide mutation in RIP-affected sequences
De novo RIP analysis of a fungal repeat unit first requires the identification of the most affected CpN di-nucleotides. The MATE transposon repeat family of Aspergillus nidulans and the Ty1 Copia-like transposon family of the Basidiomycete Microbotryum violaceum were analysed by RIPCAL. A. nidulans MATE repeats are reported to exhibit a dual preference for CpG→TpG and CpA→TpA RIP mutation in descending order of magnitude . The Ty1 repeats of M. violaceum were reported to exhibit a strong preference for CpG→TpG di-nucleotide RIP mutation . High levels of CpG→TpG and CpA→TpA RIP mutation were detected in the MATE transposons (Additional file 2). RIPCAL also detected the CpG→TpG bias in the Ty1 repeats of M. violaceum (Additional file 1). Hood et al have reported preferential mutation of the tri-nucleotide TpCpG to TpTpG in Ty1 , however RIPCAL is not currently designed to detect a tri-nucleotide RIP bias.
Di-nucleotide frequency and index analysis of RIP mutation in Stagonospora nodorum
Analysis of Stagonospora nodorum repeat families for evidence of RIP ranked by CpA↔TpA dominance.
1.70 ± 0.03
0.74 ± 0.02
Ubiquitin conjugating enzyme
1.75 ± 0.02
0.49 ± 0.02
Non LTR transposon
1.76 ± 0.06
0.56 ± 0.05
1.99 ± 0.03
0.48 ± 0.02
Non LTR transposon
1.35 ± 0.11
1.26 ± 0.13
2.68 ± 0.18
0.69 ± 0.04
Non-array rDNA repeats ≥ 1 kb
1.76 ± 0.07
0.58 ± 0.07
1.90 ± 0.06
0.40 ± 0.04
1.73 ± 0.08
0.27 ± 0.04
1.70 ± 0.03
0.52 ± 0.02
Sub-telomeric repeat/Transposon remnant
1.65 ± 0.04
0.44 ± 0.04
1.97 ± 0.18
0.44 ± 0.10
1.87 ± 0.19
0.56 ± 0.18
1.67 ± 0.04
0.46 ± 0.05
2.04 ± 0.03
0.38 ± 0.02
2.22 ± 0.13
0.39 ± 0.03
1.84 ± 0.07
0.36 ± 0.03
2.06 ± 0.07
0.33 ± 0.04
1.91 ± 0.03
0.74 ± 0.01
1.87 ± 0.04
0.33 ± 0.02
Sub-telomeric repeat/Gypsy-like transposon
1.92 ± 0.08
0.30 ± 0.03
2.08 ± 0.09
0.94 ± 0.02
rDNA repeats in tandem array
1.85 ± 0.03
0.28 ± 0.02
1.93 ± 0.07
0.27 ± 0.03
1.85 ± 0.09
0.31 ± 0.03
3.55 ± 0.39
0.25 ± 0.03
Non-array rDNA repeats < 1 kb
2.16 ± 0.15
0.31 ± 0.03
2.10 ± 0.18
0.24 ± 0.05
0.83 ± 0.01
1.25 ± 0.00
Genomic regions not corresponding to repeat matches
Alignment-based analysis of RIP mutation in Stagonospora nodorum
Repeat families of S. nodorum were aligned and scanned for RIP-like di-nucleotide changes using RIPCAL. RIP mutation statistics for all repeat families of S. nodorum are summarised in Additional file 2. Alignment-based analysis indicated that the dominant form of CpN-targeted RIP mutation in S. nodorum repeats was CpA to TpA as observed by index analysis. High levels of CpT to TpT mutation were also observed in some repeat classes (Additional file 2).
In this analysis we introduce a statistic called 'RIP dominance'. RIP dominance is the ratio of a particular CpN↔TpN RIP mutation over the sum of the other 3 alternative CpN↔TpN mutations within a multiple alignment (or sub-alignment). This was used to determine the relative strength of CpA to TpA type RIP mutations in S. nodorum (Table 2).
The alignment-based method employed by RIPCAL is an efficient, accurate and reliable method of RIP detection and characterisation. RIPCAL successfully detected the presence and absence of RIP in the positive and negative N. crassa control sequences. RIPCAL also accurately determined the preferential CpN mutation bias in RIP-affected sequences. The CpG bias in Ty1 repeats of M. violaceum and the dual CpG and CpA bias in MATE repeats of A. nidulans were also identified consistent with previously published results [13, 24].
Di-nucleotide frequency, RIP index and alignment-based analyses all indicated that CpA to TpA mutation was the dominant CpN-targeted mutation in the repeat families of S. nodorum. This preference is common to most known RIP-affected fungi. The high incidence of CpT to TpT mutation detected by alignment is less common, but has been observed in Magnaporthe grisea accompanying CpA-targeted mutation in RIP-affected sequences [4, 5]. However high levels of CpT to TpT mutation within S. nodorum short rDNA repeats, which are presumably unaffected by RIP, suggest that CpT-targeted mutation may not related to RIP in S. nodorum. Further experimental evidence is required to confirm to relevance of CpT to TpT mutation to RIP in S. nodorum and other Fungi.
RIPCAL alignment-based analysis displays the physical distribution of RIP along an alignment as shown in for the X0 repeat family in Figure 2 and the Y1/rDNA repeat family in Figure 4. This allows detection of individual repeats with anomalous changes, such as the single RIP-affected tandem rDNA repeat (Figure 4A). The lack of CpA to TpA mutation within the tandem rDNA repeats adds further supporting evidence for RIP-resistance within the rDNA nucleolus organiser region (NOR) [2, 25]. However, the RIP-affected tandem repeat, located at the terminus of the rDNA array suggests that protection from RIP within the NOR has a finite range.
The length of a S. nodorum repeat class and the degree of RIP mutation did not appear to be related (Table 2). This was highlighted by X48, a short sub-telomeric repeat, which had a high CpA↔TpA dominance score of 1.82. Its length of 275 bp was well below the 400 bp length considered the minimum for RIP in N. crassa  and the 280 bp length of the S. nodorum short rDNA repeats (which do not display CpA to TpA changes). Alignment-based analysis predicted that sub-telomeric repeats were among the most RIP-susceptible. This may explain the high CpA↔TpA dominance of X48 as chromosome ends may be physically more accessible to the molecular RIP machinery. Alternatively, the X48 repeat may be recognised in conjunction with adjacent repeats as a single unit. Unlike the NOR, fungal telomeres do not appear to be immune to RIP. RIP-like changes have also been reported in the sub-telomeric gene TLH of Magnaporthe oryzae .
We present RIPCAL as a versatile and efficient tool for the analysis of RIP which simplifies existing index-based analyses and adds alignment-based RIP analysis as a feasible alternative for whole genome analysis. These analyses highlight significant deficiencies in index-based methods of RIP detection. The alignment-based approach is biologically relevant and reveals novel features and predictions that can be tested experimentally in appropriate organisms. Sifting through the expected flood of fungal genome sequences for RIP-like phenomena may provide insights on fungal lifestyle, genomics and evolution.
RIPCAL has multiple modes of operation involving different combinations of RIP index and alignment-based methods. RIPCAL can be run in either command-line or graphical modes and is Perl-based. It is also compiled as a Windows executable. Dependent on the analysis method, RIPCAL accepts sequence input in Fasta format, pre-aligned sequence input in Fasta or ClustalW format and repeat coordinate input in either version 2 or 3 GFF format. If pre-aligned input is not provided, RIPCAL can interface with a local installation of ClustalW . Refer to Additional file 3 for more detailed information.
RIP index analysis
When using GFF input, RIP index data for repeat features was compared to a non-repetitive control family. If repeat family information is contained within the GFF input (via the target attribute) then this process was also separated by family. Fold changes between repeat families and the control were determined by ΔNpN = (repeat NpN count)/(control NpN count), where NpN represents any di-nucleotide combination.
RIP index sequence scan
Where two windows meeting the above criteria overlap, the predicted sub-region was extended (Additional file 3). Sub-regions were subject to a minimum size threshold (default 300 bp) reflecting the existence of an experimentally observed size threshold for RIP . Non-published indices were excluded by default, but can be employed as additional/replacement criteria using thresholds based on results obtained in this paper (Additional file 2). This method can be used to predict de novo ancient/non-repeated RIP-affected sequences. However, caution should be used with this method as the above threshold values are calibrated for RIP in N. crassa.
RIPCAL's alignment-based analysis indicates the presence, type and location of a putatively RIP-generated mutation within each copy of a repeat family. The input is accepted as Fasta or as both Fasta and GFF inputs. "Repeat_region" features in the GFF input were aligned by family via ClustalW (Additional file 4, Additional file 5). The prevalence of internal direct repeats within repeat families can result in poor alignment. Therefore the ClustalW default parameters have been adjusted for fast alignment, pairwise window length = 50 and k-tuple word-size = 2 to improve repeat family alignment. In some cases custom alignment parameters or manual alignment curation was used and is recommended. Sequence-only inputs are also accepted as pre-aligned Fasta files. It is assumed for sequence-only inputs that all sequences belong to the same family.
Aligned sequences are compared to a model sequence which can be either a sequence with highest total G:C content in the alignment, the alignment consensus or a user-defined sequence. The default model selection method is highest total G:C content. As RIP mutations deplete the G:C content, this default is assumed to select the least RIP-affected sequence as the model. RIPCAL also provides alternative methods of model selection, one of which is to define a majority consensus of the aligned sequences. The degenerate nucleotide code is used if two or more nucleotides are present in equal frequency (Additional file 3). The third option is for the model to be user-defined. This would be appropriate if the non-RIP-affected sequence was known, as in the case of experimentally transformed strains.
Following alignment and choice of model, the mutation frequencies are compared along the alignment for each sequence. Where the consensus sequence is degenerate, the probability of mutation at that location is added to the total count. The final output is a repeat family alignment and corresponding RIP frequency graph in GIF format. A summary of RIP mutation type versus total sequence divergence per sequence is also generated based on the alignment.
Validation of alignment-based RIP analysis
The alignment-based method was tested using the Tad1 transposon and 5S rDNA repeats from Neurospora crassa as positive and negative controls for detection of RIP mutation. These sequences [GenBank:L25662, GenBank:AF181821] were mapped to the N. crassa genome (release 7)  via RepeatMasker . The genomic matches were compared via RIPCAL for RIP mutation. Aspergillus nidulans MATE transposon sequences  [GenBank:.BK001592, GenBank:.BK001593, GenBank:.BK0015924, GenBank:.BK001595, GenBank:.BK001596, GenBank:X78051] were compared via RIPCAL using MATE-9 [GenBank:.BK001592] as the model for comparison to test for detection of non-classical (non Cpa→TpA) RIP mutation. RIP mutation of Ty1 Copia-like transposons of Mycrobotryum violaceum [PopSet:55418573] was also analysed using the degenerate consensus model to observe RIP detection in sequences with a known tri-nucleotide mutation bias .
RIP Analysis of S. nodorum de novo repeat families
Other CpN↔TpN dominance equations (Additional file 2) were of a similar format to the one above (8).
Time of Operation
All data was generated on a 2.99 GHz Dual-core ×64 Intel PC with 2 GB RAM. The combined run-time of the di-nucleotide and alignment-based analyses for the S. nodorum whole genome assembly was approximately 4 hours. Pre-aligned inputs with few sequences (i.e. < 20) can be expected to complete under a minute.
We thank Barbara Howlett and Thierry Rouxel for their expert advice and the Grains Research & Development Corporation (Barton, ACT, 2600, Australia) for providing the PhD scholarship to JKH.
- Soanes DM, Richards TA, Talbot NJ: Insights from sequencing fungal and oomycete genomes: what can we learn about plant disease and the evolution of pathogenicity? Plant Cell 2007, 19(11):3318–3326. 10.1105/tpc.107.056663PubMed CentralView ArticlePubMedGoogle Scholar
- Selker EU: Premeiotic instability of repeated sequences in Neurospora crassa . Annu Rev Genet 1990, 24: 579–613. 10.1146/annurev.ge.24.120190.003051View ArticlePubMedGoogle Scholar
- Selker E, Cambareri E, Jensen B, Haack K: Rearrangement of duplicated DNA in specialised cells of Neurospora . Cell 1987, 51: 741–752. 10.1016/0092-8674(87)90097-3View ArticlePubMedGoogle Scholar
- Ikeda K, Nakayashiki H, Kataoka T, Tamba H, Hashimoto Y, Tosa Y, Mayama S: Repeat-induced point mutation (RIP) in Magnaporthe grisea : implications for its sexual cycle in the natural field context. Mol Microbiol 2002, 45(5):1355–1364. 10.1046/j.1365-2958.2002.03101.xView ArticlePubMedGoogle Scholar
- Nakayashiki H, Nishimoto N, Ikeda K, Tosa Y, Mayama S: Degenerate MAGGY elements in a subgroup of Pyricularia grisea : a possible example of successful capture of a genetic invader by a fungal genome. Mol Gen Genet 1999, 261(6):958–966. 10.1007/s004380051044View ArticlePubMedGoogle Scholar
- Graia F, Lespinet O, Rimbault B, Dequard-Chablat M, Coppin E, Picard M: Genome quality control: RIP (repeat-induced point mutation) comes to Podospora . Mol Microbiol 2001, 40(3):586–595. 10.1046/j.1365-2958.2001.02367.xView ArticlePubMedGoogle Scholar
- Idnurm A, Howlett BJ: Analysis of loss of pathogenicity mutants reveals that repeat-induced point mutations can occur in the Dothideomycete Leptosphaeria maculans . Fungal Genet Biol 2003, 39(1):31–37. 10.1016/S1087-1845(02)00588-1View ArticlePubMedGoogle Scholar
- Neuveglise C, Sarfati J, Latge JP, Paris S: Afut1, a retrotransposon-like element from Aspergillus fumigatus . Nucleic Acids Res 1996, 24(8):1428–1434. 10.1093/nar/24.8.1428PubMed CentralView ArticlePubMedGoogle Scholar
- Hua-Van A, Hericourt F, Capy P, Daboussi MJ, Langin T: Three highly divergent subfamilies of the impala transposable element coexist in the genome of the fungus Fusarium oxysporum . Mol Gen Genet 1998, 259(4):354–362. 10.1007/s004380050822View ArticlePubMedGoogle Scholar
- Hua-Van A, Langin T, Daboussi MJ: Evolutionary history of the impala transposon in Fusarium oxysporum . Mol Biol Evol 2001, 18(10):1959–1969.View ArticlePubMedGoogle Scholar
- Julien J, Poirier-Hamon S, Brygoo Y: Foret1, a reverse transcriptase-like sequence in the filamentous fungus Fusarium oxysporum . Nucleic Acids Res 1992, 20(15):3933–3937. 10.1093/nar/20.15.3933PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen ML, Hermansen TD, Aleksenko A: A family of DNA repeats in Aspergillus nidulans has assimilated degenerated retrotransposons. Mol Genet Genomics 2001, 265(5):883–887. 10.1007/s004380100484View ArticlePubMedGoogle Scholar
- Hood ME, Katawczik M, Giraud T: Repeat-induced point mutation and the population structure of transposable elements in Microbotryum violaceum . Genetics 2005, 170(3):1081–1089. 10.1534/genetics.105.042564PubMed CentralView ArticlePubMedGoogle Scholar
- Farman ML: Telomeres in the rice blast fungus Magnaporthe oryzae : the world of the end as we know it. FEMS microbiology letters 2007, 273(2):125–132. 10.1111/j.1574-6968.2007.00812.xView ArticlePubMedGoogle Scholar
- Braumann I, Berg M, Kempken F: Repeat induced point mutation in two asexual fungi, Aspergillus niger and Penicillium chrysogenum . Curr Genet 2008.Google Scholar
- Cambareri E, Singer M, Selker E: Recurrence of repeat-induced point mutation (RIP) in Neurospora crassa . Genetics 1991., 1:Google Scholar
- Watters MK, Randall TA, Margolin BS, Selker EU, Stadler DR: Action of repeat-induced point mutation on both strands of a duplex and on tandem duplications of various sizes in Neurospora . Genetics 1999, 153(2):705–714.PubMed CentralPubMedGoogle Scholar
- Cambareri E, Jensen B, Schabtach E, Selker E: Repeat-induced G-C to A-T mutations in Neurospora . Science 1989, 244: 1571–1575. 10.1126/science.2544994View ArticlePubMedGoogle Scholar
- Margolin BS, Garrett-Engele PW, Stevens JN, Fritz DY, Garrett-Engele C, Metzenberg RL, Selker EU: A methylated Neurospora 5S rRNA pseudogene contains a transposable element inactivated by repeat-induced point mutation. Genetics 1998, 149(4):1787–1797.PubMed CentralPubMedGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, et al.: The genome sequence of the filamentous fungus Neurospora crassa . Nature 2003, 422(6934):859–868. 10.1038/nature01554View ArticlePubMedGoogle Scholar
- Hane JK, Lowe RGT, Solomon PS, Tan KC, Schoch CL, Spatafora JW, Crous PW, Kodira C, Birren BW, Galagan JE, et al.: Dothideomycete-plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum . Plant Cell 2007, 19(11):3347–3368. 10.1105/tpc.107.052829PubMed CentralView ArticlePubMedGoogle Scholar
- Solomon PS, Lowe RGT, Tan KC, Waters ODC, Oliver RP: Stagonospora nodorum : cause of Stagonospora nodorum blotch of wheat. Mol Plant Pathol 2006, 7(3):147–156. 10.1111/j.1364-3703.2006.00326.xView ArticlePubMedGoogle Scholar
- Cambareri EB, Helber J, Kinsey JA: Tad1–1, an active LINE-like element of Neurospora crassa . Mol Gen Genet 1994, 242(6):658–665. 10.1007/BF00283420View ArticlePubMedGoogle Scholar
- Clutterbuck AJ: MATE transposable elements in Aspergillus nidulans : evidence of repeat-induced point mutation. Fungal Genet Biol 2004, 41(3):308–316. 10.1016/j.fgb.2003.11.004View ArticlePubMedGoogle Scholar
- Perkins DD, Metzenberg RL, Raju NB, Selker EU, Barry EG: Reversal of a Neurospora translocation by crossing over involving displaced rDNA, and methylation of the rDNA segments that result from recombination. Genetics 1986, 114(3):791–817.PubMed CentralPubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al.: Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23(21):2947–2948. 10.1093/bioinformatics/btm404View ArticlePubMedGoogle Scholar
- Selker EU, Tountas NA, Cross SH, Margolin BS, Murphy JG, Bird AP, Freitag M: The methylated component of the Neurospora crassa genome. Nature 2003, 422(6934):893–897. 10.1038/nature01564View ArticlePubMedGoogle Scholar
- RepeatMasker Open-3.0[http://www.repeatmasker.org]
- Price AL, Jones NC, Pevzner PA: De novo identification of repeat families in large genomes. Bioinformatics 2005, 21(Suppl 1):i351–358. 10.1093/bioinformatics/bti1018View ArticlePubMedGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenetic and genome research 2005, 110(1–4):462–467. 10.1159/000084979View ArticlePubMedGoogle 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.