Intracellular zinc concentration and localization are strictly regulated by two main protein components, metallothioneins and membrane transporters. In mammalian cells, two membrane transporters family are involved in intracellular zinc homeostasis: the uptake transporters calledSLC39or Zip family and the efflux transporters calledSLC30or ZnT family. ZnT proteins are members of the cation diffusion facilitator (CDF) family of metal ion transporters.
From genomic databanks analysis, we identified the full-length sequences of two novelSLC30genes,SLC30A8andSLC30A10, extending theSLC30family to ten members. We used an expressed sequence tag (EST) data mining strategy to determine the pattern of ZnT genes expression in tissues.In silicoresults obtained for already studied ZnT sequences were compared to experimental data, previously published. We determined an overall good correlation with expression pattern obtained by RT-PCR or immunomethods, particularly for highly tissue specific genes.
The method presented herein provides a useful tool to complete gene families from sequencing programs and to produce preliminary expression data to select the proper biological samples for laboratory experimentation.
Zinc is involved in many cellular processes as a cofactor of numerous enzymes, nuclear factors and hormones and as an intra- and intercellular signal ion [1,2], and hence, is a very important component of cell viability. However, since both zinc excess and deficiency could be toxic, local intracellular zinc concentrations must be strictly regulated. The two main protein components involved in zinc homeostasis are metallothioneins, zinc transporters , and specific, gated, zinc permeable membrane spanning channels [4,5]. Metallothioneins play an important role in zinc transport, storage and distribution . Zinc transporters are transmembrane proteins, which ensure zinc ions carriage across biological membranes. Some transporters allow intracellular uptake of zinc, while others permit cellular efflux of zinc. Proteins involved in cellular uptake of zinc have been characterized in plants, yeast and mammals . In mammalian cells, seven homologous zinc export proteins, named ZnT-1 to -7 have been discovered (for review see ). These proteins are members of theSLC30solute carrier subfamily of the CDF family (Cation Diffusion Facilitator), and share the same predicted structure, with six membrane-spanning domains and a histidine-rich intracellular loop between helixes IV and V, excepted for ZnT-6 which retains a serine-rich loop . It is still controversial whether mammalian ZnT proteins are truly transporters or proteins controlling zinc transportation through other channels . However, recent works demonstrated that bacterial ZitB and CzcD proteins, two members of the CDF family are antiporters catalyzing the obligatory exchange of Zn2+or Cd2+for K+and H+with a 1:1 stoichiometry [10,11].
ZnT-1 is an ubiquitous zinc transporter located in the plasma membrane and ensures zinc efflux from the cell . ZnT-2 equally confers zinc resistance, although it is located in acidic endosomal/lysosomal vesicles and allows vesicular zinc accumulation inside the cell . ZnT-3 and ZnT-4 are more closely related to ZnT-2 than ZnT-1. ZnT-3 is tissue specific and mainly located in brain, in the membranes of zinc-rich synaptic vesicles within mossy fiber boutons of hippocampus  and in testis . Conversely, ZnT-4 is expressed ubiquitously , but higher levels of ZnT-4 are found in brain, mammary glands and epithelial cells . This transporter has been shown to be essential in mammary epithelia for regulating milk zinc content in mice . ZnT-5 is an ubiquitous zinc transporter localized in intracellular non-acidotropic vesicles and found to be abundantly expressed in pancreatic beta cells . A sixth member of the ZnT family, ZnT-6 has been described and is responsible for the relocation of cytoplasmic zinc into the trans Golgi network and the vesicular compartment . Recently, ZnT-7 was also described as a Golgi apparatus protein involved in accumulation of zinc .
From genomic databanks analysis, we identified two novelSLC30genes,SLC30A8andSLC30A10. During the preparation of this article, anotherSLC30gene,SLC30A9, appeared in Genbank [21,22] under the accession number BC016949, extending the family to 10 genes. However, the homology for this latter gene to the otherSLC30sequences is very low. To further characterize these new genes and prove the validity of this method, we took advantage of the ever-increasing wealth of information available through the human expressed sequence tag database (dbEST). Assuming that cDNA libraries used for EST sequencing are representative of all mRNA transcripts in a given tissue , we determinedSLC30family mRNA transcript levels in different tissues by EST database analysis for all the already known ZnTs (except for ZnT-9) and compared theirin silicoexpression profiles with experimental data on human tissues. For most cases, the experimental data correlate within silicoanalysis. Hence, this strategy provides valuable informations and the method presented herein is a useful tool to complete gene families from sequencing programs and to produce preliminary expression data before selecting the proper biological samples for laboratory experimentation.
Results and Discussion
An approach for discovering new genes is to search the whole human genome sequence for homologous sequences of known genes or known gene families byin silicomethods. Recent publications demonstrate the efficiency of this technique to find new genes [24,25]. Using the different already known ZnT cDNA and protein sequences in human, mouse or rat as a bait for a BLASTN or a TBLASTN search of the human genome databanks, we discovered two DNA sequences encoding new putative zinc transporters belonging to the ZnT family. These new genes were namedSLC30A8andSLC30A10, encoding the proteins designated ZnT-8 and ZnT-10 respectively. HumanSLC30A8cDNA was found in the contig AC027419, which allowed us to localize theSLC30A8gene to human chromosome 8 at the position q24.11 (Table1). The gene contained 8 exons, spanned 37 kb and is predicted to code for a 40.8 KDa protein (Fig.1). The sequence data reported for human ZnT-8 mRNA was submitted to Genbank under the accession number AY117411. HumanSLC30A10cDNA was found in the contig AC093562, which allowed us to localize theSLC30A10gene to human chromosome 1 at the position q41. The gene contained 4 exons, spanned 15 kb and was predicted to code for a 52.7 KDa protein (Fig.1). The sequence for ZnT-10 mRNA was submitted to Genbank under the accession number BK004163. We also localized ZnT-2 gene (SLC30A2) in human genome to chromosome 1 at the position p36.11 (Table1) by homology with the rat ZnT-2 sequence . The predicted cDNA and protein sequences are identical to the NM_032513 nucleotide and NP_115902 protein entries of Entrez database.
HumanSLC30family genes. Chromosomal localisations, number of exons, number of histidine residues between the fourth and the fifth predicted transmembrane domains (TMIV and TMV) and Genbank accession number ofSLC30genes.
Histidine residues between TMIV and TMV
Genbank Accession Number
a) this study nd: not determined
We then aligned the predicted sequence of the nine human ZnT proteins1with the ClustalW program. When compared, all the proteins of the family are predicted to have a conserved structure, with a common pattern composed of 6 transmembrane helices and a histidine-rich domain between helices IV and V. Both N- and C termini are predicted to be located on the cytoplasmic side of the plasma membrane. Alignments of amino acids composing the fifth and sixth transmembrane domains illustrate this homology (Fig.2A). However, despite very well conserved residues, the homology between amino acid sequences can differ from one protein to another. For example it was known that ZnT-5 exhibited 15 transmembrane domains, but the region homologous to the members of the family is located in the carboxyl-terminal portion and is predicted to adopt six membrane-spanning domains. The histidine-rich loop is replaced by a serine-rich loop for ZnT-6. We report the presence of a loop rich in basic residues for ZnT-10, while ZnT-8 keeps the characteristic histidine-rich loop (Fig.2B). The histidine content is also very different, from no histidine residue for ZnT-10 to 20 histidine residues for ZnT-7 (Table1).
Using the amino acids alignment, a phylogenetic tree for the 10 ZnT sequences was calculated by the neighbour-joining method (Fig.3). Zip-2, a zinc membrane transporter belonging to the SLC39 family was used as an outgroup. From this analysis, we can delineate three subfamilies: ZnT-1 and ZnT-10; ZnT-5 and ZnT-7; ZnT-2, ZnT-3, ZnT-4 and ZnT-8. This result was confirmed by similarity analysis of the amino acid protein sequences. The subfamily ZnT-2, -3, -4, -8 exhibited the highest homologies, with the highest score of 53.5 % between ZnT-2 and ZnT-8. The homology between ZnT-1 and ZnT-10 is high with a score of 48.3. But, ZnT-5, -6 and -7 are less homologous, with a highest score of 27.8 % between ZnT-5 and ZnT-7. ZnT-9 has the lowest homology with the other ZnTs. Despite an overall shared topological structure, the similarity between the subfamilies is relatively low.
Thein silicocharacterization ofSLC30tissue expression pattern was performed by an expressed sequence tag (EST) data mining strategy. The predictedSLC30transcripts (ORF, 5' and 3' UTRs) were queried against the human EST database using BLASTN. We obtained a total of 426 significant hits with a bit score >150 and an E-value <0.001. The ESTs sequenced from libraries prepared from pooled cDNA tissues or derived from other libraries were rejected for the analysis. 306 EST sequences were considered for tissue origin. The numbers of EST clustered for each Zinc transporter sequence are very different, ranging from 8 ESTs forSLC30A10to 100 forSLC30A5. The average number of ESTs per library ranges from 1.1 forSLC30A1to 5.6 forSLC30A8(Table2).
Statistics forin silicoanalysis ofSLC30genes tissue expression. The total EST are the significant ESTs (a bit score >150 and an E-value <0.001) resulting from BLAST analysis withSLC30sequences. The EST sequenced from libraries prepared from pooled cDNA tissues or derived from other libraries were rejected. The number of libraries was determined from information regarding each cDNA library retrieved from EST data.
Number of libraries
We foundSLC30A1expression in 18 tissues out of 36, indicating a very wide expression pattern (Fig.4). ZnT-1 was demonstrated to display a broad tissue distribution. It is particularly abundant in intestine, liver  and in the brain . On the other hand, no ZnT-1 transcripts were expressed at detectable levels inlamina propriaintestinal cells or in many kidney cells. ZnT-1 gene is controlled at the transcriptional level by zinc status. The elevation of extracellular zinc concentration results in a rapidly and dramatically increase ofSLC30A1mRNA levels, mediated by the transcription factor MTF-1, a sensor of zinc level . So, the basal evaluation ofSLC30A1expression by EST analysis may not reflect the real expression which depends on extracellular conditions.
ForSLC30A2expression, we calculated a highest level in placenta and high levels in eye, kidney and ovary. Experimental results indicate an expression of ZnT-2 in intestine, kidney, seminal vesicles and testis . In rats,SLC30A2mRNA expression is limited to small intestine, kidney, placenta and liver, whileSLC30A2mRNA levels were increased several fold only in small intestine, liver and kidney upon a single oral dose of zinc . The very high level ofSL30A2expression in placenta presumably indicated an important role of ZnT-2 transporter in zinc exchange between maternal tissues and foetus.
The results forSLC30A3display a good correlation betweenin silicoanalysis and experimental data. We determined a restricted expression with very high levels of expression in brain and testis, an expression pattern previously identified by northern blot and reverse transcriptase-PCR analysis . In brain,SLC30A3mRNA is most abundant in the cerebral cortex and in synaptic vesicle membranes within mossy fiber boutons in the hippocampus . Zinc is secreted from these vesicles in response to high frequency stimulations [29,30].
ESTs forSLC30A4were founded only in few tissues (B-cells, muscle, ovary, parathyroid gland, stomach and testis) and the correlation with experimental data is very poor. The highest level was calculated for testis tissue sample. Znt-4 was first thought to play an important role in milk secretion. A nonsense mutation, leading to a truncated form of ZnT-4, is responsible for the inherited zinc deficiency in the lethal milk (lm) mouse [17,31]. In thelmmouse, the maternal milk does not contain enough zinc for the newborn mice to live. ZnT-4 is constitutively expressed in human breast epithelial cells . However, in human no difference in ZnT-4 expression levels was observed between lactating and resting breasts. In rats, ZnT-4 is expressed ubiquitously and was refractory to changes in zinc intake . ZnT-4 is also expressed in polarized enterocytes, in which it is localized in the membrane of intracellular vesicles, the majority of which concentrates in the basal cytoplasmic region. The protein was not founded in proliferating cells of the crypt, but was detected in differenciated enterocytes of thevilli, the apparition corresponding to the junction crypt/villi.
From EST analysis results,SLC30A5is ubiquitously expressed, with high levels in kidney, liver, pancreas, brain, skin, bone marrow and T-cells. We determined the presence and the level ofSLC30A5mRNA by PCR amplification of cDNA libraries prepared from different human tissues [seeadditional file 1]. As expected from calculated data, aSLC30A5specific product was detected at a high level in nearly all kind of tissues, thus confirming the previously published results .
Fromin silicoanalysis,SLC30A6displayed the highest levels in germinal B-cells and colon and high levels in eye and lung.In vivo,SLC30A6mRNA has been detected in liver, brain, small intestine and kidney. Western blot analysis indicated that ZnT-6 is present in mouse brain, small intestine, kidney and lung .
Low levels of expression were calculated forSLC30A7, excepted in the colon and the eye.SLC30A7mRNA (Northern-blot) or PCR products were detected in the heart, liver, spleen, plasma blood leukocytes, small intestine, kidney, brain, lung, ovary, prostate and testis at a very low level ( and seeadditional file 1). Recently, we demonstrated an induction ofSLC30A7expression by extracellular zinc deficiency .
SLC30A8had a very high expression restricted to the pancreas. We detected a faint signal for four other tissues. We then analyzedSLC30A8gene expression by PCR using a panel of 24 cDNAs prepared from different tissues. A specific PCR product was only detected in pancreatic tissue sample (seeadditional file 2). This last result is highly correlated within silicoanalysis.
From EST analysis results,SLC30A10had a restricted expression to fetal liver and fetal brain. It is the first zinc transporter predicted to have a fetal restricted expression.SLC30A10andSLC30A1have a high homology. At birth, ZnT-1 protein is nearly undetectable and ZnT-1 expression increases at the end of the first postnatal week . So, we speculate that ZnT-10 could play a role comparable to that of ZnT-1 during fetal development.
From genomic databanks analysis, we identified two novelSLC30genes,SLC30A8andSLC30A10, extending theSLC30family to ten members. We determined an overall good correlation of ZnTin silicogene expression with expression patterns obtained by RT-PCR or immunomethods, particularly for highly tissue-specific genes. As the average number of ESTs recovered per library was relatively low (few copies of ZnT sequence per library), we can not definitively conclude that tissues without ESTs for a given ZnT do not express this gene at all. We have also to keep in mind that the zinc status of the cells and, hence, the adaptative mechanisms to extracellular zinc concentrations were usually unknown for sample tissues used for RT-PCR experiments or EST sequencing programs. In conclusion, this method provides a useful tool to complete gene families from sequencing programs and to produce preliminary expression data to select the proper biological samples for laboratory experimentation.
Genomic DNA and translated sequence analysis
BLASTN or TBLASTN  analysis were performed on the NCBI web server  and the ENSEMBL web server  with human ZnT-1 (AF323590 and Q9Y6M5), mouse ZnT-1 (U17132 and Q60738), rat ZnT-1 (U17133 and Q62720), rat ZnT-2 (U50927 and Q62941), human ZnT-3 (NM_003459 and Q99726), mouse ZnT-3 (U76007 and P97441), human ZnT-4 (AF025409 and O14863), mouse ZnT-4 (O35149), rat ZnT-4 (O55174), human ZnT-5 (AY089991 and AAM09099), human ZnT-6 (NM_017964 and NP_060434) and human ZnT-7 (AY094606 and AAM21969) DNA or protein sequences. The TMpred computer program  was used to analyze and predict the transmembrane potential of the predicted protein sequences .
EST analysis forin silicodetermination of SLC30 genes tissue expression
TheSLC30sequences (ORF, 5' and 3' UTRs) were used for a BLASTN search of the human EST database through NCBI BLAST web service . The significant ESTs (bit score >150 and E-value <0.001) were sorted and information regarding each cDNA library was retrieved from either the human Unigene databank  or from the respective company catalogue. Libraries prepared from pooled tissue samples or derived from other libraries were rejected from the analysis. The frequency of each mRNA transcript for a given tissue was calculated. The ESTs were also analyzed by the Gene2EST program  to precisely locate the 5' and 3' starting ends of the transcript and the spliced variants .
Sequences alignment and phylogeny
PredictedHomo sapiensZnT protein sequences were aligned using the clustalW program . For phylogenetic analysis, bootstrapping (2000 replicate sets) and calculation of the consensus tree were performed with the DAMBE program by the neighbour-joining method . Bootstrap analysis is based on multiple re-sampling of the original data and is the commonest method of estimating the degree of confidence in the topology of phylogenetic trees. Zip-2 protein sequence (NP_055394) was used as an outgroup.
Expression in human tissues
The presence and the level of SLC30A5, SLC30A7 and SLC30A8 expression were determined by PCR amplification of cDNA libraries (Origene Technologies, Boston, MA) prepared from different human tissues with Expand high fidelity DNA polymerase (Roche, Meylan, France). The specific primers were: 5'-CTG CTT TAG TCA TGG GAC TTT TTG C and 5'-TAG AAC CTG GCA GGC ATC TTT AAT C for ZnT-5, 5'-GAT GTC CAC CAT GAG AGA CTG CTT C and 5'-CAT AAA GTC CAG AAG TGC TGT TCC TG for ZnT-7, 5'-GAT GCT GCC CAC CTC TTA ATT GAC and 5'-CCA AGA CCA GGA TGG AAA AGA TGA for ZnT-8 and 5'-CCA AGG CCA ACC GCG AGA AGA TGA C and 5'-AGG GTA CAT GGT GGT GCC GCC AGA C for β-actin. The products were analyzed by agarose gel electrophoresis, stained with ethidium bromide and photographed under UV light with a CCD camera.
expressed sequence tag
Cation Diffusion Facilitator
This work has been supported by a grant from the "Centre Evian Pour l'Eau" to MS.
Hosie AM, Dunne EL, Harvey RJ, Smart TG:Zinc-mediated inhibition of GABA(A) receptors: discrete binding sites underlie subtype specificity.Nat Neurosci2003,6(4):362–369.View ArticlePubMed
Chimienti F, Aouffen M, Favier A, Seve M:Zinc homeostasis-regulating proteins: new drug targets for triggering cell fate.Curr Drug Targets2003,4(4):323–338.View ArticlePubMed
Jia Y, Jeng JM, Sensi SL, Weiss JH:Zn2+ currents are mediated by calcium-permeable AMPA/kainate channels in cultured murine hippocampal neurones.J Physiol2002,543:35–48.View ArticlePubMed
Kerchner GA, Canzoniero LM, Yu SP, Ling C, Choi DW:Zn2+ current is mediated by voltage-gated Ca2+ channels and enhanced by extracellular acidity in mouse cortical neurones.J Physiol2000,528:39–52.View ArticlePubMed
Cousins RJ, McMahon RJ:Integrative aspects of zinc transporters.J Nutr2000,130(5S Suppl):1384S–1387S.PubMed
Eide DJ:The SLC39 family of metal ion transporters.Pflugers Arch2003,14:769–800.
Palmiter RD, Huang L:Efflux and compartmentalization of zinc by members of the SLC30 family of solute carriers.Pflugers Arch2004,447(5):744–751.View ArticlePubMed
Ohana E, Segal D, Palty R, Ton-That D, Moran A, Sensi SL, Weiss JH, Hershfinkel M, Sekler I:A sodium zinc exchange mechanism is mediating extrusion of zinc in mammalian cells.J Biol Chem2004,279(6):4278–4284.View ArticlePubMed
Chao Y, Fu D:Kinetic study of the antiport mechanism of an Escherichia coli zinc transporter, ZitB.J Biol Chem2004,279(13):12043–12050.View ArticlePubMed
Guffanti AA, Wei Y, Rood SV, Krulwich TA:An antiport mechanism for a member of the cation diffusion facilitator family: divalent cations efflux in exchange for K+ and H+.Mol Microbiol2002,45(1):145–153.View ArticlePubMed
Palmiter RD, Findley SD:Cloning and functional characterization of a mammalian zinc transporter that confers resistance to zinc.Embo J1995,14(4):639–649.PubMed
Palmiter RD, Cole TB, Findley SD:ZnT–2, a mammalian protein that confers resistance to zinc by facilitating vesicular sequestration.Embo J1996,15(8):1784–1791.PubMed
Wenzel HJ, Cole TB, Born DE, Schwartzkroin PA, Palmiter RD:Ultrastructural localization of zinc transporter–3 (ZnT–3) to synaptic vesicle membranes within mossy fiber boutons in the hippocampus of mouse and monkey.Proc Natl Acad Sci U S A1997,94(23):12676–12681.View ArticlePubMed
Palmiter RD, Cole TB, Quaife CJ, Findley SD:ZnT–3, a putative transporter of zinc into synaptic vesicles.Proc Natl Acad Sci U S A1996,93(25):14934–14939.View ArticlePubMed
Liuzzi JP, Blanchard RK, Cousins RJ:Differential regulation of zinc transporter 1, 2, and 4 mRNA expression by dietary zinc in rats.J Nutr2001,131(1):46–52.PubMed
Huang L, Gitschier J:A novel gene involved in zinc transport is deficient in the lethal milk mouse.Nat Genet1997,17(3):292–297.View ArticlePubMed
Kambe T, Narita H, Yamaguchi-Iwai Y, Hirose J, Amano T, Sugiura N, Sasaki R, Mori K, Iwanaga T, Nagao M:Cloning and characterization of a novel mammalian zinc transporter, ZnT–5, abundantly expressed in pancreatic beta cells.J Biol Chem2002,271(21):19094–19055.
Huang L, Kirschke CP, Gitschier J:Functional characterization of a novel mammalian zinc transporter, ZnT6.J Biol Chem2002,277(29):26389–26385.View ArticlePubMed
Kirschke CP, Huang L:ZnT7, a novel mammalian zinc transporter, accumulates zinc in the Golgi apparatus.J Biol Chem2003,278(6):4096–4102.View ArticlePubMed
Sim DL, Chow VT:The novel human HUEL (C4orf1) gene maps to chromosome 4p12-p13 and encodes a nuclear protein containing the nuclear receptor interaction motif.Genomics1999,59(2):224–233.View ArticlePubMed
Pires Martins R, Leach RE, Krawetz SA:Whole-body gene expression by data mining.Genomics2001,72(1):34–42.View ArticlePubMed
Wittenberger T, Schaller HC, Hellebrand S:An expressed sequence tag (EST) data mining strategy succeeding in the discovery of new G-protein coupled receptors.J Mol Biol2001,307(3):799–813.View ArticlePubMed
Bathgate RA, Samuel CS, Burazin TC, Layfield S, Claasz AA, Reytomas IG, Dawson NF, Zhao C, Bond C, Summers RJ, Parry LJ, Wade JD, Tregear GW:Human relaxin gene 3 (H3) and the equivalent mouse relaxin (M3) gene. Novel members of the relaxin peptide family.J Biol Chem2002,277(2):1148–1157.View ArticlePubMed
McMahon RJ, Cousins RJ:Regulation of the zinc transporter ZnT–1 by dietary zinc.Proc Natl Acad Sci U S A1998,95(9):4841–4846.View ArticlePubMed
Sekler I, Moran A, Hershfinkel M, Dori A, Margulis A, Birenzweig N, Nitzan Y, Silverman WF:Distribution of the zinc transporter ZnT–1 in comparison with chelatable zinc in the mouse brain.J Comp Neurol2002,447(3):201–209.View ArticlePubMed
Langmade SJ, Ravindra R, Daniels PJ, Andrews GK:The transcription factor MTF–1 mediates metal regulation of the mouse ZnT1 gene.J Biol Chem2000,275(44):34803–34809.View ArticlePubMed
Li Y, Hough CJ, Suh SW, Sarvey JM, Frederickson CJ:Rapid translocation of Zn(2+) from presynaptic terminals into postsynaptic hippocampal neurons after physiological stimulation.J Neurophysiol2001,86(5):2597–2604.PubMed
Li Y, Hough CJ, Frederickson CJ, Sarvey JM:Induction of mossy fiber ––> Ca3 long-term potentiation requires translocation of synaptically released Zn2+.J Neurosci2001,21(20):8015–8025.PubMed
Piletz JE, Ganschow RE:Zinc deficiency in murine milk underlies expression of the lethal milk (lm) mutation.Science1978,199(4325):181–183.View ArticlePubMed
Michalczyk AA, Allen J, Blomeley RC, Ackland ML:Constitutive expression of hZnT4 zinc transporter in human breast epithelial cells.Biochem J2002,364:105–113.PubMed
Devergnas S, Chimienti F, Naud N, Pennequin A, Coquerel Y, Chantegrel J, Favier A, Seve M:Differential regulation of zinc efflux transporters ZnT–1, ZnT–5 and ZnT–7 gene expression by zinc levels: a real-time RT-PCR study.Biochem Pharmacol2004,in press.
Nitzan YB, Sekler I, Hershfinkel M, Moran A, Silverman WF:Postnatal regulation of ZnT–1 expression in the mouse brain.Brain Res Dev Brain Res2002,137(2):149–157.View ArticlePubMed
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ:Basic local alignment search tool.J Mol Biol1990,215(3):403–410.PubMed
Gemund C, Ramu C, Altenberg-Greulich B, Gibson TJ:Gene2EST: a BLAST2 server for searching expressed sequence tag (EST) databases with eukaryotic gene-sized queries.Nucleic Acids Res2001,29(6):1272–1277.View ArticlePubMed
Thompson JD, Higgins DG, Gibson TJ:CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.Nucleic Acids Res1994,22(22):4673–4680.View ArticlePubMed
Xia X, Xie Z:DAMBE: software package for data analysis in molecular biology and evolution.J Hered2001,92(4):371–373.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.