Visual pigments in a living fossil, the Australian lungfish Neoceratodus forsteri
© Bailes et al; licensee BioMed Central Ltd. 2007
Received: 02 April 2007
Accepted: 25 October 2007
Published: 25 October 2007
One of the greatest challenges facing the early land vertebrates was the need to effectively interpret a terrestrial environment. Interpretation was based on ocular adaptations evolved for an aquatic environment millions of years earlier. The Australian lungfish Neoceratodus forsteri is thought to be the closest living relative to the first terrestrial vertebrate, and yet nothing is known about the visual pigments present in lungfish or the early tetrapods.
Here we identify and characterise five visual pigments (rh1, rh2, lws, sws1 and sws2) expressed in the retina of N. forsteri. Phylogenetic analysis of the molecular evolution of lungfish and other vertebrate visual pigment genes indicates a closer relationship between lungfish and amphibian pigments than to pigments in teleost fishes. However, the relationship between lungfish, the coelacanth and tetrapods could not be absolutely determined from opsin phylogeny, supporting an unresolved trichotomy between the three groups.
The presence of four cone pigments in Australian lungfish suggests that the earliest tetrapods would have had a colorful view of their terrestrial environment.
Sarcopterygian fish gave rise to the first tetrapods and are represented today by the lungfishes (the Australian, Neoceratodus forsteri; the African, Protopterus spp.; and the South American, Lepidosiren paradoxa), and the coelacanth, Latimeria chalumnae. The relationship between all early Sarcopterygii remains controversial and highly debated despite the advent of phylogenetic analysis of nucleotide and amino acid sequences [8–10]. Some molecular analyses of sarcopterygian phylogeny reveal that the lungfish is more related to tetrapods than to the coelacanth, L. chalumnae , while others present an unresolved trichotomy between all three groups . Fossil forms of the lungfish family Ceratodontidae (genus Ceratodus) first appear in the fossil record in the Triassic period ( for review). The genus Neoceratodus (approx. 4 species, of which N. forsteri is the sole survivor) is found in the fossil record from the Lower Cretaceous period 135 million years ago (mya) and therefore N. forsteri lays claim to being the oldest surviving vertebrate genus . Consequently, the visual system of N. forsteri may represent an evolutionary design most closely reflecting that present just prior to the emergence of land vertebrates in the Devonian period.
The Australian lungfish Neoceratodus forsteri was thought to have poor eyesight due to its small eye size, low spatial resolving power [13, 14], sluggish behaviour in captivity [15–17] and ability to detect prey using electroreception . However, recent work on the retina of N. forsteri has revealed four morphologically-distinct photoreceptor types (one rod and at least three cones), some containing colored intracellular filters that are otherwise only found in terrestrial Orders [13, 19]. Although a partial sequence of the African lungfish Protopterus spp. rh1 opsin gene has been previously published [; Genbank: AF369054], nothing else is known about lungfish opsins. The other extant sarcopterygian fish, L. chalumnae, possesses only two functional opsin genes, rh1 and rh2 and lives in a photon-limited deep-sea environment . Conversely, N. forsteri inhabits a brightly lit, shallow freshwater habitat more similar to the environment from which terrestrial evolution occurred . This prompted us to investigate the complement of opsins expressed in N. forsteri in order to trace the evolution of photoreception in ancestral tetrapods.
Results and discussion
We have characterised the full length cDNA coding sequences of five visual opsin genes (rh1, rh2, lws, sws1 and sws2), one from each of the five image-forming vertebrate opsin groups, in N. forsteri. In addition, 3' RACE experiments reveal multiple transcripts of rh1 utilising different polyadenylation (polyA) signal sequences within the rh1 3' untranslated region (UTR). All five opsins are expressed in the retina of sub-adult fish and the deduced amino acid (aa) sequences yield polypeptides ranging in size between 351–356 aa's. Lungfish opsins share highly conserved residues known to be important in visual pigment function such as a glutamate counterion at aa site 113 (numbering follows that of bovine rhodopsin ) and a lysine residue forming the chromophore binding site at aa site 296.
More than ten independent clones (from PCR experiments using both degenerate primers and in 3' and 5' rapid amplification of cDNA ends) from at least two individual lungfish were sequenced for each opsin gene family and, while there was evidence of polymorphism in the gene pool (and/or possible sequencing error) with a variation between clones of up to 0.5%, there was no evidence for gene duplications in individuals. We therefore conclude that we did not sequence paralogous opsins within the five main vertebrate opsin groups in N. forsteri and that the lungfish genome encodes a single copy of each opsin gene. This is in contrast to opsin gene duplications in teleost fish. Ray-finned fish (Actinopterygii) underwent a whole genome duplication event around 350 mya, after the divergence of the Sarcopterygii, and many lineages of teleost fish are tetraploid ( and references therein), whereas N. forsteri is a diploid animal . In addition, some species of teleost fish appear to have undergone multiple opsin gene duplications independently . These have accumulated subsequent amino acid changes, resulting in differences in the maximum sensitivity of opsins. These opsins can then be preferentially expressed to fine-tune the animal's spectral sensitivity to environmental light, thus reflecting a degree of visual plasticity [3, 26]. For example, African cichlid fish determine their spectral sensitivity by means of preferential expression by up to seven available cone opsin genes [2, 26, 27].
Species and Genbank accession numbers  of opsin nucleotide sequences and deduced amino acids used in phylogenetic analyses.
Common name reference in Fig. 3
Genbank accession number
Italian wall lizard
African clawed frog
Phylogenetic analysis of rh opsin nucleotide sequences does not favour the coelacanth as a closer relative to tetrapods than lungfish and further supports an unresolved trichotomy between the coelacanth, lungfish and early tetrapods, which varies according to the gene family investigated and the method of analysis . For example, the phylogenetic tree produced using the NJ method with nucleotide sequences places lungfish rh2 together with tetrapod opsins rather than the coelacanth opsin gene. Within the rh1 rhodopsin group, however, lungfish and coelacanth sequences are placed together as a sister group to the tetrapod rod pigments. Bayesian inference favours the lungfish and coelacanth forming a sister group to the tetrapods from both rh1 and rh2 sequences (Fig. 4). Conversely, comparison of amino acid sequences using the NJ method [see Additional file 1] supports coelacanth opsins forming a sister group to both teleost and tetrapod opsins, while lungfish rh1 and rh2 genes are placed as both a sister group to teleost fish opsins (rh1) and tetrapod opsins (rh2).
The fish-tetrapod transition occurred in a space of <20 million years around 400 mya  and a large array of genes or whole genome sequences may therefore be needed to resolve the trichotomy between sarcopterygian fish . It is unfortunate that the coelacanth genome does not encode a functional sws1 pigment as sws1 genes have recently been proposed as a good marker for vertebrate phylogenies .
The lungfish visual system
Our results demonstrate that the full complement of the orthologues of the known vertebrate photoreceptive visual pigment genes are expressed in the retina of the Australian lungfish, the nearest, and most primitive, extant relative to the land vertebrates. Partial sequences of rh1 and lws have also been successfully amplified from the more derived African lungfish Protopterus sp. (unpublished data, H.J. Bailes, W.L. Davies, A.E.O. Trezise and S.P. Collin). An earlier molecular study of visual pigments in the retina of a sarcopterygian fish (the coelacanth Latimeria chalumnae) found only two opsins, rh1 and rh2 in addition to a pseudogene Ψsws1 . The coelacanth inhabits a photon-limited deep-sea environment and the move to this low-light, shortwave-shifted environment meant that the subsequent loss of three cone visual pigments (LWS, SWS1 & SWS2 used for color vision in bright-light environments in other vertebrates ) was not a selective disadvantage, revealing a correlation between loss of opsin genes and the spectral habitat/functional needs of a species . Conversely, the retention of all five vertebrate opsin families in N. forsteri suggests that the dipnoan lineage has lived in a brightly-lit colorful environment throughout its evolutionary history, making it a more appropriate model organism for the early tetrapod visual system than the coelacanth.
The characterization of four cone opsins reveals the potential for tetrachromatic vision in lungfish, although behavioral work is needed to verify if lungfish can discriminate objects based on differences in chromatic hue. Multiple opsins and a range of colored intraocular filters in the retina of N. forsteri [13, 19] suggest that it is adapted for diurnal vision, in contrast to earlier reports that adult lungfish are crepuscular [16, 35]. Lungfish have a mostly carnivorous diet  and may utilise color vision in prey capture or reproductive behaviour. These findings indicate that the first tetrapods probably possessed eyes adapted for chromatic diurnal vision, with all five opsins expressed in the retina. Colored oil droplets within the photoreceptors would have also filtered the incident light, enhancing color discrimination by reducing the spectral overlap of pigment absorbance curves [13, 36].
One adult lungfish was caught by hook and line in the Mary River near Tiaro, Queensland (Queensland Fisheries Management Authority Permit No. PRM01599G). Three subadult and two juvenile fish were bred in captivity and donated by Prof. Jean Joss from Macquarie University, Sydney, Australia. Animals were sacrificed using an overdose of benzocaine dissolved in acetone (according to animal ethics guidelines of the University of Queensland AEC No. ANAT/436/04/ARC). Australian lungfish are listed as 'Vulnerable' under the Australian Commonwealth Environment Protection and Biodiversity Conservation Act 1999 and as such only a limited number of animals were available for this work.
Dissected retinae were placed in RNAlater solution at 4°C. Total RNA was extracted using a Macherey-Nagel Nucleospin-RNA II kit (Machery-Nagel GmbH & Co. K.G.) for individual adult and subadult eyes, or from pooled left and right eyes for each juvenile fish used. Total RNA was converted to cDNA using Superscript II (Invitrogen Corp.) and random 9-mer or 16-mer oligo-dT primers. A series of degenerate primers were designed to conserved regions specific to each of the five vertebrate retinal opsin families (rh1, rh2, lws, sws1 and sws2). Primers were used in nested PCR on cDNA using standard methods . Amplified fragments were cloned into pBluescript vector (Stratagene Inc.) and sequenced at the Australian Genome Research Facility (AGRF Ltd.). Once sequenced fragments were obtained and assembled, rapid amplification of cDNA ends (RACE) at both the 5' and 3' ends was performed to produce the full-length sequence of visual opsin mRNA. Specific primers were then designed to the 5' and 3' ends of each opsin and used in conjunction with a proof-reading enzyme (Phusion; Finnzymes Oy) to verify the full length coding sequence of each opsin identified (rh1 forward TTA GGA GCT GCA ACC ATG AAC GGA ACA GAG, rh1 reverse [polyA transcript1] GCT TGT GGG TTT GTC TGC AGA TTG CAA TGG, rh1 reverse [polyA transcript2] CCG TTC TAT GCC TTC TCT ACC GGT TTC TTG, rh2 forward ACC AAC AGC AGT AGT GTA TTC GCA GCA AAG, rh2 reverse AGG GAT ACT TGG CTT GAG GAG ACT GAA GAG, lws forward ATA GAG ACA GAG AGG GAG AGA TGG CTG AAC, lws reverse CGC CGT ACA GTC ATT GCT TGT GAA ATA GTG, sws1 forward AGC AGA CAG AAG ATG TCA GGG GAA GAA GAG, sws1 reverse GCC ATA ACA CAA CTA AGG GGC CAT CAC TTC, sws2 forward CCG GGT TAC ACA CCA CTA CAA GTC AAC TAC, sws2 reverse AAT GGC TGG AGG AGA CCG AAG AGA CCT GAG). The verification of each opsin with these primer sets was carried out using cDNA from additional individuals from those used in original opsin identification experiments. In this way, each sequence was confirmed as transcribed in at least two individuals.
Full-length coding sequences were compared against known sequences from other fish and tetrapods in the Genbank database, and against an outgroup of Drosophila melanogaster rh4 (Table 1). A codon-matched alignment was carried out using ClustalW (European Bioinformatics Institute) and by manual inspection. Phylogenetic analysis of nucleotides and corresponding amino acids was performed using MEGA3  software and the Neighbour-joining (NJ) method  with 1000 bootstrap replications. The Tamura-Nei  model of DNA evolution was used, with complete deletion and a homogeneous pattern of nucleotide substitution among lineages and uniform rates of nucleotide substitution across all sites. The NJ method with 1000 bootstrap replications, a homogeneous pattern of substitution among lineages and uniform rates of evolution was also used to infer phylogeny from deduced amino acids using a Poisson correction substitution model. Phylogeny was also deduced using Bayesian inference via a Metropolis-coupled Markov chain Monte Carlo (MCMC) simulation using MrBayes 3.1 software [41, 42]. A general time-reversible model (GTR ) of DNA evolution was used, with a gamma-shaped rate variation with a proportion of invariable sites. Two simultaneous runs were performed for 300,000 generations with chains sampled taken every 1000 generations. The first 75 trees sampled (25%) were discarded as burnin. Consensus trees were extracted in Treeview 1.6.6 .
The authors wish to thank Robert Lucas of The University of Manchester, UK, for helpful comments on the manuscript. We are also grateful to Jean Joss of Macquarie University, Australia and Mike Bennett of the University of Queensland, Australia for the donation of juvenile and subadult lungfish. Mal Jepson and Helen Nicoll kindly helped in catching wild adult lungfish. This work was funded by an Australian Research Council grant (DP0209452) to SPC and AEOT, and an International Postgraduate Research Scholarship and Travel Scholarship from The University of Queensland to HJB.
- Yokoyama S, Yokoyama R: Adaptive evolution of photoreceptors and visual pigments in vertebrates. Annu Rev Ecol Syst. 1996, 27: 543-567. 10.1146/annurev.ecolsys.27.1.543.View ArticleGoogle Scholar
- Spady TC, Seehausen O, Loew ER, Jordan RC, Kocher TD, Carleton KL: Adaptive molecular evolution in the opsin genes of rapidly speciating Cichlid species. Mol Biol Evol. 2005, 22: 1412-1422. 10.1093/molbev/msi137.View ArticlePubMedGoogle Scholar
- Trezise AEO, Collin SP: Opsins: Evolution in waiting. Curr Biol. 2005, 15: R794-R796. 10.1016/j.cub.2005.09.025.View ArticlePubMedGoogle Scholar
- Nathans J, Thomas D, Hogness DS: Molecular genetics of human color vision: the genes encoding blue, green, and red pigments. Science. 1986, 232: 193-202. 10.1126/science.2937147.View ArticlePubMedGoogle Scholar
- Johnson RL, Grant KB, Zankel TC, Boehm MF, Merbs SL, Nathans J, Nakanishi K: Cloning and Expression of Goldfish Opsin Sequences. Biochemistry. 1993, 32: 208-214. 10.1021/bi00052a027.View ArticlePubMedGoogle Scholar
- Chang WS, Harris WA: Sequential genesis and determination of cone and rod photoreceptors in Xenopus. J Neurobiol. 1998, 35: 227-244. 10.1002/(SICI)1097-4695(19980605)35:3<227::AID-NEU1>3.0.CO;2-0.View ArticlePubMedGoogle Scholar
- Collin SP, Knight MA, Davies WL, Potter IC, Hunt DM, Trezise AEO: Ancient colour vision: multiple opsin genes in the ancestral vertebrates. Curr Biol. 2003, 13: 864-865. 10.1016/j.cub.2003.10.044.View ArticleGoogle Scholar
- Carroll RL: Patterns and Processes of Vertebrate Evolution. 1997, Cambridge, England: Cambridge University PressGoogle Scholar
- Brinkmann H, Venkatesh B, Brenner S, Meyer A: Nuclear protein-coding genes support lungfish and not the coelacanth as the closest living relatives of land vertebrates. Proc Natl Acad Sci USA. 2004, 101: 4899-4905. 10.1073/pnas.0400609101.View ArticleGoogle Scholar
- Takezaki N, Figueroa F, Zaleska-Rutczynska Z, Takahata N: The phylogenetic relationship of tetrapod, coelacanth, and lungfish revealed by the sequences of 44 nuclear genes. Mol Biol Evol. 2004, 21: 1512-1524. 10.1093/molbev/msh150.View ArticlePubMedGoogle Scholar
- Marshall CR: A List of Fossil and Extant Dipnoans. J Morphol Supp. 1986, 1: 15-23. 10.1002/jmor.1051900405.View ArticleGoogle Scholar
- Kemp A, Molnar RE: Neoceratodus forsteri from the lower Cretaceous of New South Wales, Australia. J Paleontol. 1981, 55: 211-217.Google Scholar
- Bailes HJ, Robinson SR, Trezise AEO, Collin SP: Morphology, characterization and distribution of retinal photoreceptors in the Australian lungfish Neoceratodus forsteri (Krefft, 1870). J Comp Neurol. 2006, 494: 381-397. 10.1002/cne.20809.View ArticlePubMedGoogle Scholar
- Bailes HJ, Trezise AEO, Collin SP: The number, morphology and distribution of retinal ganglion cells and optic axons in the Australian lungfish Neoceratodus forsteri (Krefft 1870). Vis Neurosci. 2006, 23: 257-273. 10.1017/S0952523806232103.View ArticlePubMedGoogle Scholar
- Dean B: Notes on the living specimens of the Australian lungfish, Ceratodus forsteri, in the Zoological Society's collection. Proc Zool Soc Lond. 1906, 168-178.Google Scholar
- Kemp A: The Biology of the Australian Lungfish, Neoceratodus forsteri (Krefft 1870). J Morphol. 1986, 181-198. 10.1002/jmor.1051900413. Supp 1
- Simpson R, Kind P, Brooks S: Trials of the Queensland Lungfish. Nat Aust Winter. 2002, 36-43.Google Scholar
- Watt M, Evans CS, Joss JMP: Use of electroreception during foraging by the Australian lungfish. Anim Behav. 1999, 58: 1039-1045. 10.1006/anbe.1999.1216.View ArticlePubMedGoogle Scholar
- Robinson SR: Early vertebrate color-vision. Nature. 1994, 367: 121-10.1038/367121a0. [http://www.nature.com/nature/journal/v367/n6459/pdf/367121a0.pdf]View ArticleGoogle Scholar
- Venkatesh B, Mark EV, Brenner S: Molecular synapomorphies resolve evolutionary relationships of extant jawed vertebrates. Proc Natl Acad Sci USA. 2001, 98: 11382-11387. 10.1073/pnas.201415598.PubMed CentralView ArticlePubMedGoogle Scholar
- Yokoyama S, Zhang H, Radlwimmer FB, Blow NS: Adaptive evolution of color vision of the Comoran coelacanth (Latimeria chalumnae). Proc Natl Acad Sci USA. 1999, 96: 6279-6284. 10.1073/pnas.96.11.6279.PubMed CentralView ArticlePubMedGoogle Scholar
- Palczewski K, Kumasaka T, Hori T, Behnke CA, Motoshima H, Fox BA, Le Trong I, Teller DC, Okada T, Stenkamp RE, Yamamoto M, Miyano M: Crystal structure of rhodopsin: a G protein-coupled receptor. Science. 2000, 289: 739-745. 10.1126/science.289.5480.739.View ArticlePubMedGoogle Scholar
- Christoffels A, Koh EGL, Chia J-M, Brenner S, Aparicio S, Venkatesh B: Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes. Mol Biol Evol. 2004, 21: 1146-1151. 10.1093/molbev/msh114.View ArticlePubMedGoogle Scholar
- Rock J, Eldridge M, Champion A, Johnston P, Joss J: Karyotype and nuclear DNA content of the Australian lungfish, Neoceratodus forsteri (Ceratodontidae: Dipnoi). Cytogenet Cell Genet. 1996, 73: 187-189.View ArticlePubMedGoogle Scholar
- Matsumoto Y, Fukamachi S, Mitani H, Kawamura S: Functional characterization of visual opsin repertoire in Medaka (Oryzias latipes). Gene. 2006, 371: 268-278. 10.1016/j.gene.2005.12.005.View ArticlePubMedGoogle Scholar
- Parry JW, Carleton KL, Spady T, Carboo A, Hunt DM, Bowmaker JK: Mix and match color vision: tuning spectral sensitivity by differential opsin gene expression in Lake Malawi cichlids. Curr Biol. 2005, 15: 1734-1739. 10.1016/j.cub.2005.08.010.View ArticlePubMedGoogle Scholar
- Spady TC, Parry JW, Robinson PR, Hunt DM, Bowmaker JK, Carleton KL: Evolution of the cichlid visual palette through ontogenetic subfunctionalization of the opsin gene arrays. Mol Biol Evol. 2006, 23: 1538-1547. 10.1093/molbev/msl014.View ArticlePubMedGoogle Scholar
- Al-Ubaidi MR, Pittler SJ, Champagne MS, Triantafyllos JT, McGinnis JF, Baehr W: Mouse opsin: Gene structure and molecular basis of multiple transcripts. J Biol Chem. 1990, 265: 20563-20569.PubMedGoogle Scholar
- Hara-Nishimura I, Kondo M, Nishimura M, Hara R, Hara T: Cloning and nucleotide sequence of cDNA for rhodopsin of the squid Todarodes pacificus. FEBS Lett. 1993, 317: 5-11. 10.1016/0014-5793(93)81480-N.View ArticlePubMedGoogle Scholar
- Petersen-Jones SM, Sohal AK, Sargan DR: Nucleotide sequence of the canine rod-opsin-encoding gene. Gene. 1994, 143: 281-284. 10.1016/0378-1119(94)90111-2.View ArticlePubMedGoogle Scholar
- Lim J, Chang JL, Tsai HJ: A second type of rod opsin cDNA from the common carp (Cyprinus carpio). Biochim Biophys Acta. 1997, 1352: 8-12.View ArticlePubMedGoogle Scholar
- Edwalds-Gilbert G, Veraldi KL, Milcarek C: Alternative poly (A) site selection in complex transcription units: means to an end?. Nucl Acids Res. 1997, 25: 2547-2561. 10.1093/nar/25.13.2547.PubMed CentralView ArticlePubMedGoogle Scholar
- von Schantz M, Lucas RJ, Foster RG: Circadian oscillation of photopigment transcript levels in the mouse retina. Brain Res (Mol Brain Res). 1999, 72: 108-114. 10.1016/S0169-328X(99)00209-0.View ArticleGoogle Scholar
- van Hazel I, Santini F, Muller J, Chang BSW: Short-wavelength sensitive opsin (SWS1) as a new marker for vertebrate phylogenetics. BMC Evol Biol. 2006, 6: 97-10.1186/1471-2148-6-97.PubMed CentralView ArticlePubMedGoogle Scholar
- Grigg GC: Studies on the Queensland lungfish, Neoceratodus forsteri (Krefft). III. Aerial respiration in relation to habits. Aust J Zool. 1965, 13: 413-421. 10.1071/ZO9650413.View ArticleGoogle Scholar
- Vorobyev M: Coloured oil droplets enhance colour discrimination. Proc R Soc Lond Ser B: Biol Sci. 2003, 270: 1255-1261. 10.1098/rspb.2003.2381.View ArticleGoogle Scholar
- Sambrook J, Russell D: Molecular cloning: A laboratory manual. 2001, Cold Spring Harbor: Cold Spring Harbor Laboratory PressGoogle Scholar
- Kumar S, Tamura T, Nei M: MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10: 512-526.PubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Lanave C, Preperata GS, Saccone C, Serio G: A new method for calculating evolutionary substitution rates. J Mol Evol. 1984, 20: 86-93. 10.1007/BF02101990.View ArticlePubMedGoogle Scholar
- Page RDM: TREEVIEW: An application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
- Bowmaker JK, Hunt DM: Molecular Biology of Photoreceptor Spectral Sensitivity. Adaptive Mechanisms in the Ecology of Vision. Edited by: Archer MA, Djamgoz MBA, Loew ER, Partridge JC, Vallerga S. 1999, Dordrecht: Kluwer Academic Publishers, 439-462.View ArticleGoogle Scholar
- The National Institutes of Health Genetic Sequence Database. [http://www.ncbi.nlm.nih.gov/Genbank/index.html]
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.