BACKGROUND
the aanat family is part of the large and diverse superfamily of gcn5-like acetyltransferases which use accoa as the acetyl donor and share a common accoa binding fold  <cit> .
members of the aanat family share limited sequence identity and are divided into two groups: vertebrate aanats; and, the non-vertebrate aanats.
the latter are found in fungi, protists, and bacteria and lack defining physical characteristics of vertebrate aanat  <cit> .
the biological role of vertebrate aanat is to acetylate serotonin in the synthesis of melatonin   <cit> .
vertebrate aanat is associated with biological timing: daily changes in the activity of this enzyme regulate the daily rhythm in melatonin synthesis, which is essential for optimal temporal coordination of biological functions with night/day and seasonal changes and for photic entrainment  <cit> .
the central role of vertebrate aanat in biological timing has earned it the moniker 'the timezyme'  <cit> .
non-vertebrate aanats are thought to play a detoxifying role by neutralizing arylalkylamines  <cit>  and a role in dna biology by acetylating polyamines  <cit> .
genomes of vertebrates contain a single copy of the aanat gene, except for teleost fish, some of which possess up to three paralogs  <cit> , and cows which possess two paralogs .
members of the aanat family are not in the available genomes of hemichordates and urochordates, which leaves open the question of when vertebrate aanat first appeared in chordates.
a striking characteristic of vertebrate aanat is that it is consistently expressed at significant levels in only two tissues, both of which are photosensitive organs, the pineal gland and retina.
this fits with the evidence that pinealocytes and retinal photoreceptors evolved from a common ancestral photodetector  <cit> .
the vertebrate aanat has a neutral ph optimum and exhibits high selectivity for arylalkylamines  <cit> .
vertebrate aanats encode proteins that have several highly conserved structural characteristics  <cit> , which facilitate arylalkylamine acetylation and regulation.
these features include flanking regulatory regions, which mediate rapid changes in enzyme activity; a pair of histidines which facilitate catalysis  <cit> ; and, a proline-containing tripeptide in a floppy loop, which confers a high catalytic rate through an effect on substrate binding  <cit> .
vertebrate aanats also have high selectivity for arylalkylamines conferred by the binding pocket.
non-vertebrate type aanats are found in the genomes of most fungi, many unicellular eukaryotes, and a variety of bacteria  <cit> .
the proteins encoded by these genes do not contain the characteristic structural features of vertebrate aanats.
genes similar to non-vertebrate aanats are not found in vertebrates.
amphioxus aanats are of interest in understanding the evolution of vertebrate aanat because cephalochordates and vertebrates form a clade  <cit> .
accordingly, the study of the aanat from this organism might provide new insights into the transition from the non-vertebrate members of the aanat family to the vertebrate forms.
sequence in the trace archives has indicated that a non-vertebrate type member of the aanat family is present in amphioxus  <cit> .
the subsequent availability of the near complete amphioxus genome  <cit>  has allowed us to comprehensively analyze the aanat family in amphioxus, as presented here.
the findings of this study indicate that rapid expansion of non-vertebrate forms occurred in amphioxus and that a vertebrate form is absent.
the substrate preference profile of the encoded enzymes suggests to us that they are likely to play a role in detoxification or in biogenic amine inactivation, or both.
methods
animals
branchiostoma floridae  and branchiostoma lanceolatum  were used; the former were obtained from gulf specimen marine lab  and the latter from the bay of banyuls sur mer   <cit> .
b. floridae heads were used to generate the cdna used to clone full length bfaanatα and bfaanatδ' that were used for enzyme characterization.
b. lanceolatum were used for in vivo studies of gene expression and enzyme activity.
they were cultured in sand-filtered seawater under a natural lighting cycle and temperature, as indicated elsewhere  <cit> .
intensity of the light was  <dig> lux.
for the experiment with constant conditions, animals were housed for  <dig> days in constant light and/or darkness and then were used for the activity assay.
identification of amphioxus aanats
bfaanat homologs were identified through a tblastn search of genomic sequences available at: http://genome.jgi-psf.org <cit> .
seven genes for aanat homologs  for cloning of blaanat cdnas.
qpcr primers for specific gene products were based on the b. lanceolatum clones.
sequence similarity search and phylogenetic analysis
the aanat homologs were identified by searching the non-redundant database of protein sequences  using the psi blast program  <cit>  and by searching the genomes and proteomes available at jgi http://genome.jgi-psf.org/euk_cur <dig> html using the tblastn and blastp programs, respectively.
additionally, aanat homologs were detected in the genomes of the red algae gracilaria changii, segmented worm alvinella pompejana, mediterranean mussel mytilus galloprovincialis  and the basal fungus rhizopus oryzae http://www.broadinstitute.org/annotation/genome/rhizopus_oryzae/multihome.html.
representatives from different groups including bacteria, fungi, protists, and animals were selected for phylogenetic analysis, and a multiple alignment of the respective protein sequences was constructed using the muscle program  <cit> .
a maximum likelihood phylogenetic tree was constructed using the treefinder program  <cit>  by optimizing a default starting tree constructed using the neighbor-joining method with the whelan and goldman  empirical model of substitutions  <cit>  to select the best substitution model, we compared the likelihoods of trees constructed using  <dig> models: wag , jtt  <cit> , vt  <cit> , blosum  <cit> , dayhoff  <cit> , cprev  <cit> , rtrev  <cit> , pmb  <cit>  by running the treefinder program  <cit> ; this test showed the best fit with the data for the wag model.
the reliability of the internal tree branches was estimated using the lr-elw bootstrap method  <cit> .
a constrained large scale tree topology was created using the treeview program  <cit> .
alternative tree topologies were compared with the approximately unbiased test  <cit>  using the treefinder program.
evolutionary distances were estimated using the distance matrix created by the protdist program of phylip package  <cit>  based on a truncated alignment of aanats and using the jones-taylor-thornton  substitution model  <cit> .
the mean differences  and identities between the major taxonomic groups of aanats were calculated using the mega program  <cit> .
prokaryotic expression of recombinant proteins
pgex 4t- <dig>  was used for the recombinant expression of full-length bfaanatα and bfaanatδ'.
the constructed plasmid was sequenced to confirm identity.
the construct, which produces a gst fusion protein, was transformed into e. coli strain bl21gold  plyss .
the cells were grown at 37°c; when the od <dig> =  <dig> , the cultures were cooled to 25°c and isopropyl β-d-1-thiogalactopyranoside was added .
the cells were harvested by centrifugation after  <dig> hours of culturing  and resuspended in 2× pbs, ph  <dig> , containing  <dig> mm dtt  and a mixture of protease inhibitors .
the cells were then lysed by sonication and the resulting lysate was centrifuged .
the supernatant was mixed with glutathione-sepharose 4b affinity matrix  pre-equilibrated with buffer a. the suspension was agitated for  <dig> h and then packed into a glass column.
the column was washed with  <dig> column volumes of buffer a, followed by  <dig> column volumes of buffer containing  <dig> mm tris-hcl,  <dig>  m sodium citrate,  <dig> mm dtt and 10% glycerol, ph  <dig>  .
gst fusion protein was then eluted with  <dig> column volumes of buffer b containing  <dig> mm glutathione.
protein was concentrated and dialyzed against buffer containing  <dig>  m ammonium acetate,  <dig> mm nacl and  <dig> mm tcep; the resulting preparation was stored at -80°c.
where indicated, purified gst-free bfaanatα was used; otherwise, the gst-bfaanatα or gst-aanatδ' fusion product was used, consistent with previous studies  <cit> .
sds-page
proteins were resolved on preformed 14% tris-glycine  gels, using the manufacturer's protocol   <cit> .
the molecular mass of the proteins was estimated using rainbow™ markers .
fluorescence-based protein studies
fluorescence measurements were done on an iss pc <dig> photon-counting spectrofluorimeter  at 22°c.
to study binding properties of bfaanatα protein, the interaction between bfaanatα and coa-s-n-acetyl-7-hydroxynaphthylethylamine , a fluorescent bisubstrate inhibitor of vertebrate aanat, was analyzed  <cit> .
this was done by measuring coa-hne -dependent quenching of tryptophan fluorescence in  <dig>  m ammonium acetate,  <dig> mm nacl and  <dig> mm tcep, ph  <dig>  .
off-rates of coa-hne were monitored as described  <cit>  using protein preparations that were 90% saturated with coa-hne.
following a control incubation period, the probe was displaced with a 100-fold excess of another bisubstrate inhibitor, coa-s-n-acetyltryptamine , that is optically inactive in this analysis  <cit> .
the limiting value of koff was calculated by non-linear fitting  using the equation:
where a is the observed fluorescence anisotropy at time 't' after the initial addition of coa-t and koff is the rate constant.
determination of n-acetyltransferase activity
recombinant bfaanatα activity was measured radiochemically by measuring acetyl product formed from acetyl-coenzyme a  and an amine substrate as described  <cit> ; or, colorimetrically by measuring coash formed by incubation with accoa and an amine substrate.
coash is detected as the colored product formed by reaction with  <dig> '-dithio-bis   <cit> .
in the colorimetric assay, standard reactions were performed in a total volume of  <dig> μl containing: the substrate at indicated concentrations,  <dig>  mm accoa,  <dig>  mg/ml bovine serum albumin , and  <dig> mm ethylene diamine tetraacetic acid .
for analysis of gst-bfaanatα,  <dig>  μg of protein was included in  <dig>  m tris buffer ph  <dig> ; the activity of ovine gst-aanat, which was used as a reference standard  <cit> , was determined by incubating  <dig> ng of protein in  <dig>  m phosphate buffer, ph  <dig> .
incubation conditions were either  <dig> min at 37°c for amphioxus or  <dig> min at 37°c for ovine aanat.
enzymes were then inactivated by the addition of  <dig> μl stop solution ; absorbance  was measured after  <dig> min incubation at room temperature.
the vmax and km were calculated by non-linear fitting  using the equation:
where v is the observed enzyme velocity , vmax is the maximum enzyme velocity, s is the substrate concentration, and km is the michaelis-menten constant.
for the measurement of enzyme activity in amphioxus head tissue, the colorimetric assay  was used with  <dig> mm β-phenylethylamine  as a substrate and  <dig> μg of protein from tissue homogenized in assay buffer.
for each measurement a corresponding assay was performed in the presence of  <dig> mm coa-t  <cit> , and another in the absence of pea.
the controls used included an assay without tissue homogenate and another with a heat-inactivated  homogenate.
samples were collected at noon and midnight in the experiments studying the light/dark variations, and at noon in all other experiments.
quantitative measurement of b. lanceolatum aanat transcripts
total rna was extracted from head and from body using a trizol extraction kit .
<dig> μg was incubated with  <dig> unit of dnase i  for  <dig> min at 37°c.
following dnase inactivation , the sequence was reverse-transcribed using powerscript reverse transcriptase .
qpcr of the resulting cdnas was performed with the primer sets qaanat  using the 'light cycler fast start dna master sybr green i kit'  in a total volume of  <dig> μl as follows: 95°c , followed by  <dig> cycles of denaturation at 95°c , annealing at 60°c  and extension at 72°c .
qpcr was done using a light cycler  <dig>  .
amplification efficiency was measured using serial dilutions of plasmids containing a fragment of the corresponding blaanat gene.
data were normalized to the average value obtained with three housekeeping genes from b. lanceolatum: l <dig> , l <dig>  and actin .
the pcr conditions for the housekeeping genes were as indicated above except that the annealing temperature was 66°c for actin.
primers for the housekeeping genes are given in additional file  <dig>
tissue distribution of amphioxus aanatα studied by rt-pcr
the tissue distribution of amphioxus aanatα was determined by rt-pcr.
total rna from different parts of b. lanceolatum was extracted using the trizol method ;  <dig> μg was incubated with  <dig> unit of dnase i  for  <dig> min at 37°c.
following dnase inactivation , the rna was reverse-transcribed using powerscript reverse transcriptase  pcr amplifications of the resulting cdnas using advantage cdna polymerase  were performed with the primers for blaanatα , and the following conditions: 95°c  followed by  <dig> cycles of denaturation at 94°c , annealing at 66°c  and extension at 68°c .
in the controls, the template cdna was replaced by either water or rna that was not reverse-transcribed.
the pcr products were loaded in a 1% agarose gel, in the presence of dna/hinf i marker .
fragments of the expected size were extracted, subcloned in pgem-t easy  and sequenced.
tissue distribution of amphioxus aanatα studied by in toto hybridization
b. lanceolatum were collected at argelès-sur-mer  in july 2009; gametes were obtained by heat stimulation  <cit> .
eggs and embryos were fixed in 4% paraformaldehyde  in mops-egta buffer and processed for in toto hybridization as detailed elsewhere  <cit> , except that the chromogenic reaction was performed using bm purple   <cit> .
two digoxigenin labeled probes corresponding to the α or δζ subtypes were designed.
the former was obtained using the cloning primers caanatαf <dig> and caanatαr <dig>  and yielded a probe of  <dig> bp.
the latter was obtained using the cloning primers caanatδζf <dig> and caanatδζr <dig>  and yielded a probe of  <dig> bp.
hybridization treatments with the different sense and antisense probes were performed simultaneously.
RESULTS
sequence analysis, phyletic distribution, and phylogeny of the aanats
seven intronless aanat genes were identified in the b. floridae genome .
a second allele was found for  <dig> of the  <dig> genes in version v <dig>  of the b. floridae genome assembly.
the amino acid sequence identity among the amphioxus aanats ranged from  <dig> to 86% excluding allele/allele comparisons .
the structural features which characterize vertebrate aanats, including the flanking regulatory regions, a pair of neighboring histidines in the active site, and a tripeptide in the floppy loop, do not appear in the proteins encoded by the amphioxus aanats.
the motifs a and b, which characterize members of the gcn5-like superfamily of acetyltransferases are conserved in both amphioxus and vertebrate aanats .
the relative locations of the bfaanats on  <dig> assembly scaffolds are shown schematically in additional file  <dig>  the  <dig> genes comprise  <dig> groups based on sequence similarity - group 1: αη group 2: β group 3: γε group 4: δζ.
primers were prepared to amplify members of each of the four groups in b. lanceolatum; low sequence conservation prevented amplification of all seven.
the resulting pcr products were sequenced and designated blaanatα, blaanatα' blaanatγ, and blaanatγδ , according to the homology with individual bfaanat genes; representatives of the other genes were not obtained in this screen.
the clone blaanatδζ was designated this way because it could not be reliably ascribed to either the δ or ζ gene based on similarity to the b. floridae counterparts.
the amino acid sequence identity between the cloned pcr products for the b. lanceolatum α, α', γ and δζ isoforms and the corresponding b. floridae genes was 94%, 86%, 76%, and 73%, respectively.
comparisons of the amino acid sequences of the cloned blaanatα and blaanatδζ fragments to their b. floridae counterparts are shown in additional file  <dig>
examination of the phyletic distribution of the detected aanat family members shows that the family is represented not only in chordates and cephalochordates but also in some invertebrates such as the primitive metazoan trichoplax , the segmented worms capitella and alvinella, and the mollusk mytilus.
among non-metazoan eukaryotes, aanat homologs were detected in different fungal lineages including basidiomycetes, ascomycetes, and the basal fungi rhizopus oryzae and phycomyces blakesleeanus, and in many unicellular eukaryotes .
the protists in which aanat homologs were detected include representatives from three supergroups of eukaryotes, namely, plantae, chromalveolates and excavates  <cit> .
in plantae, aanat homologs were found in green and red algae; among chromalveolates, members of the aanat family were present in diatoms, oomycetes, and a haptophyte emiliania huxleyi; and, among excavates, an aanat homolog was identified in the heterolobosean naegleria gruberi.
aanat homologs are also present in a number of bacterial lineages including firmicutes, actinobacteria, proteobacteria and the bacteroidetes/chlorobi group .
aanats are relatively short proteins which complicates reliable phylogenetic analysis.
phylogenetic trees were constructed based on an alignment of  <dig> amino acid sequences .
eight treefinder runs were performed for  <dig> substitution models .
the ml tree constructed with the wag matrix was associated with the highest likelihood value.
the topology of this tree is incompatible with the species phylogeny .
the monophyly of all  <dig> amphioxus aanats was moderately supported and these sequences formed an apparent clade with mollusk and segmented worm that was even more strongly supported .
however, beyond this invertebrate branch, there were many anomalies in the tree.
in particular, vertebrates clustered with ascomycetes and bacteria, whereas the position of trichoplax remained uncertain.
the monophyly of fungi was not recovered, and ostreococii were separated from another green alga represented in the selected set of protists .
we tested the statistical compatibility of the aanat tree with the large-scale organismal phylogeny of eukaryotes using the constrained tree approach .
a constrained tree was created  in which all species were combined into four major clades in accord with the overall organismal phylogeny, namely, bacteria, protists, fungi and animals.
inside the animal branch, with trichoplax at the root, the amphioxus branch was placed as the sister group of vertebrates, and the segmented worms and mussel were positioned between chordates and trichoplax.
the fungal branch was placed as the sister group of animals.
the approximately unbiased test result  for the tree generated with this approach  shows that the "classical" constrained tree topology could not be rejected with confidence although the likelihood of that topology is much lower than that of the best tree shown in figure  <dig>
the results of the phylogenetic analysis of the aanats emphasize the major distinction between the vertebrate forms and those from other animals, presumably, owing to a dramatic acceleration of evolution at the base of the vertebrate branch.
beyond that conclusion, however, the results were quite uncertain due to the low information content of the alignment, between-branch differences in the evolutionary rates, and uneven sampling of different groups.
expression and characterization of bfaanatα protein
a comparison of the predicted amino acid sequences derived from the published genomic sequence and from our cloned bfaanatα  revealed differences in five amino acid residues scattered throughout .
these differences, most likely, reflect the high degree of polymorphism in b. floridae  <cit> .
the effects of these differences are unknown.
to determine the molecular mass of expressed bfaanatα, the gst-bfaanatα fusion protein was cleaved and bfaanatα was purified as described in the methods.
the final preparation contained a single major band of protein which migrated at ~ <dig> kda on sds-page .
maldi-tof analysis revealed the mass to be  <dig>  corresponding to the predicted mass based on amino acid composition.
using the same procedure we also expressed gst-bfaanatδ'.
to establish whether the expressed protein was correctly folded, we measured binding of a fluorescent bisubstrate inhibitor of the expressed protein, coa-hne; binding of the inhibitor to oaanat quenches fluorescence  <cit> .
here we found that bfaanatα quenched inhibitor fluorescence and that the expressed protein bound the inhibitor stoichiometrically .
moreover, binding was totally disrupted by treatment with  <dig> m guanidine.
accordingly, the purified preparation of bfaanatα is not denatured and possesses characteristics generally similar to those of oaanat, as regards binding of coa-hne.
to further characterize bfaanatα, off-rate analysis was done using the displacement strategy described previously  <cit> , in which bound coa-hne is displaced by coa-t, a non-fluorescent aanat inhibitor  <cit>  that does not interfere with detection of coa-hne .
this result shows that the off-rate observed for the coa-hne is >  <dig>  s- <dig>  similar to that of yeast aanat; these rates are much faster than that for oaanat  <cit> , consistent with the view that the proline-containing tripeptide in the floppy loop of vertebrate aanat enhances binding of substrate.
biochemical characterization of bfaanatα and bfaanatδ'
activity as a function of ph: gst-bfaanatα and gst-bfaanatδ' both exhibited a ph activity profile that differed from that of ovine aanat and yeast aanat as examined using  <dig> mm tryptamine as a substrate .
the ph optimum for gst-bfaanatα was determined to be  <dig> , which differs substantially from published data on fish, ovine and yeast aanat   <cit> .
the ph optimum for gst-bfaanatδ' seems to be even higher.
similar results were obtained using butylamine as a substrate .
effect of aanat bisubstrate inhibitor: enzyme activity was inhibited by the bisubstrate inhibitor coa-t  <cit>  .
the ic <dig> value  was found to be ~10-fold greater than that reported for ovine aanat  <cit> .
temperature dependence: the temperature optimum for bfaanatα was found to be 36°c .
this is similar to that of ovine aanat , yeast aanat  or pike aanat- <dig>   <cit> .
substrate preference: the relative capacities of oaanat, bfaanatα and bfaanatδ' to acetylate a diamine , alkylamines  and arylalkylamines  were determined at ph  <dig>  for oaanat, and  <dig>  for bfaanatα and bfaanatδ'  using gst-fusion proteins.
the activity of bfaanats towards the substrates tested  was relatively similar and fell within a 5-fold range of activity; the highest activity was observed with pea, however, bfaanatα and bfaanatδ' seemed to show a slight preference for hydrophobic or polar substrates, respectively.
the km and vmax values were determined for five biogenic arylalkylamines: octopamine, pea, tyramine, tryptamine and serotonin .
activity was measured at ph  <dig> .
values are calculated from data in figure  <dig>  vmax values are nmol/h/pmol enzyme and are given as the mean ± s.e.m.
.
for further details see methods.
in contrast to the relatively broad substrate preference profiles of bfaanats observed with the substrates tested, oaanat exhibited markedly higher activity towards the two indolamines, tryptamine and serotonin, and lower activity towards putrescine, propylamine and butylamine.
this points to a marked difference in substrate preference: bfaanatα and bfaanatδ' exhibit a broad capacity to acetylate amines, including arylalkylamines and a diamine, whereas oaanat exhibits a distinct preference for arylalkylamines with an aromatic group, including serotonin, pea and tryptamine.
control of the abundance of blaanat transcripts
in some vertebrate species, the nocturnal increase in aanat activity is accompanied by a ~100- fold increase in the abundance of transcripts  <cit> .
here we used qpcr to determine whether the abundance of blaanat mrnas varied on a night/day basis .
however statistically significant differences were only observed for blaanatα.
accordingly, although it does not appear that there are marked night/day differences in expression of blaanats, it is clear that small differences do occur.
analysis of tissue expression by rt-pcr
investigation of tissue expression of blaanatα by rt-pcr revealed the presence of mrna along the whole body .
the identity of the product was confirmed by sequencing .
daily changes in blaanat activity
homogenates of the heads of b. lanceolatum were used for measurement of aanat activity, using  <dig> mm pea as a substrate; coa-t was added to determine the contribution from aanat to total acetyltransferase activity.
in all conditions tested, activity was significantly higher in heads sampled at midday than in those sampled at midnight .
activity measured in the absence of added substrate was more than 50% as high as that measured in the presence of exogenous substrate ; in both cases, activity was nearly completely inhibited by coa-t, indicating that blaanat acetylated endogenous substrates during the assay.
in other studies this rhythm was confirmed, however a similar daily change was not observed in animals housed in constant darkness or constant lighting .
a detailed graph describing the inhibition of acetylation activity in homogenates by coa-t is shown in additional file  <dig>
substrates
comparison of samples collected during day and night.
the exogenous substrate is  <dig> mm pea;  <dig> mm coa-t has been used for inhibition.
data are given as the mean ± s.e.m.
.
acetylation activity under normal light/darkness changes  compared with constant light  and constant darkness  conditions.
the ph of the assay was  <dig> ; data are given as the mean ± s.e.m.
.
in toto hybridization
expression of one or both blaanats was detected by in toto hybridization at all stages examined, in varying patterns.
although it is difficult to quantify the colorimetric reactions, it seemed that the α form was more highly expressed than the δζ form at early stages, whereas the opposite seemed true at later stages.
the first specific hybridization signal was detected at the 8-cell embryo and was stronger in  <dig> out of the  <dig> cells.
at the blastula stage most of the cells exhibited a diffuse labeling with both probes; however, a group of cells exhibited a more intense dark signal .
at the early neurula stage, the labeling appeared diffuse and weak with the δζ probe, whereas some regionalization appeared with the α probe.
this tendency was even more apparent at the late neurula stage where the α probe gave a specific labeling in the antero-dorsal  and posterio-dorsal parts of the embryo, as well as antero-ventrally in the general area of the intestine.
at the premouth stage, no labeling with the α probe was detected, whereas the δζ probe hybridized from head to tail in a dorsal band, corresponding to the neural tube.
discussion
the publication of the near complete amphioxus genome has marked a new era in understanding chordate evolution  <cit> .
the studies presented in this report, which were made possible by this advance, focus on one family of genes and extends our understanding of this gene family in amphioxus and, potentially, in animals in general.
aanats in the amphioxus genome and evolution of the aanat family
analysis of the amphioxus genome indicates that it lacks a vertebrate-like aanat.
in contrast, our analysis shows that the amphioxus genome contains a set of closely related aanat homologs that are similar in their gene and protein architectures to homologs from invertebrates, protists, fungi, and bacteria.
none of these non-vertebrate homologs possesses the defining structural features of vertebrate aanats.
the discovery of the family of aanats in amphioxus sheds light on the timing of the emergence of the typical vertebrate aanats that, most likely, occurred subsequent to the radiation of the vertebrate lineage from the common ancestor with cephalochordates.
it appears that the change in the aanat gene architecture at the base of the vertebrate tree was accompanied by a substantial acceleration of the evolutionary rate which is the likely reason why the amphioxus aanat homologs are more similar to those from fungi and protists than to those from vertebrates.
the results of phylogenetic analysis of the aanat family presented here were not entirely conclusive owing, in part, to the relatively small number of phylogenetically informative positions in the multiple sequence alignment of the aanat family, and in part, to the non-uniformity of evolutionary rates, and perhaps, to a rapid divergence associated with evolutionary transitions.
the best ml  tree showed multiple deviations from the species tree but the topology corresponding to the species tree, including the cephalochordate-vertebrate clade, could not be rigorously rejected.
thus, the phylogeny of aanats might follow the organismal phylogeny of animals, taking into account the apparent multiple losses, e.g., in arthropods, nematodes, and echinodermata, as well as accelerations of evolution, particularly, in vertebrates.
previous analyses of the phyletic patterns and phylogenies of the aanats led to the suggestion that horizontal gene transfer, perhaps, even into the animal lineage or the animal-fungal ancestral lineage was a major contribution to the evolution of this family  <cit> .
that possibility still remains considering the patchy phyletic distribution of the aanats and the irregularities in the topologies of the phylogenetic trees, such as the polyphyly of fungi and the unusual assortment of protists.
however, given the lack of strong resolution in the trees, the major changes in evolutionary rates like the one in vertebrates, and the extreme propensity of the aanat family for gene loss, scenarios that involve unusual routes of horizontal transfer should be addressed with much caution.
biochemical features of bfaanatα and bfaanatδ'
the biochemical characteristics of bfaanatα and bfaanatδ' were analyzed in this study.
both appear to have significantly lower intrinsic activity at ph  <dig>  as compared to oaanat.
the substrate preference spectrum of the bfaanats is broader than that of oaanat.
as with other non-vertebrate aanats, these differences are likely to reflect differences in the primary structure of the enzyme, including the presence in vertebrate aanat of two histidines near the site of acetyl transfer.
these are thought to be part of a proton wire which accelerates catalysis by facilitating the expulsion of protons generated during acetyl transfer, thereby preventing their accumulation and the resulting acidification of the active site .
another important difference in the primary structure between vertebrate and amphioxus aanats occurs in a floppy loop which is critical for catalysis, apparently by increasing the dwell time of amine stubstrates in the binding pocket  <cit> .
the floppy loop of the vertebrate aanats is three residues longer than that of known aanat homologs, including amphioxus.
this addition appears to enhance catalysis by directing the movement of the floppy loop and preventing intramolecular interactions which inhibit binding.
as a result, the off-rate of substrate is substantially increased when these three residues are experimentally deleted  <cit> .
during the course of evolution these three residues were added and, as shown here, this was associated with a decrease in the off-rate relative to that of amphioxus and yeast aanats  <cit> .
regional expression of blaanat
our studies revealed that blaanatα is expressed throughout the body , in contrast to the pattern of expression of vertebrate aanat, which is primarily limited to structures located in the head: pineal gland, retina and in some cases the pituitary gland  <cit> .
the development of these structures is strongly dependent on expression of the transcription factor pax <dig>  pax <dig> expression in amphioxus is limited to the frontal organ and lamellar organ; the latter has been identified as a possible pineal precursor, based on its medial location and presence of photoreceptors in the tissue  <cit> .
in toto hybridization
the in toto hybridization analyses revealed expression of both α and δζ genes early in development in b. lanceolatum.
it seems that the relative strength of expression changes during development.
blaanatδζ was found to be expressed in a region close to the neural tube.
this observation and the finding that these bfaanats acetylate the biogenic amines serotonin and octopamine , suggest to us that it is reasonable to consider that these enzymes play a role in the metabolism of these neurotransmitters.
the possible roles of amphioxus aanats
the unexpected finding that several of the seven aanats in the amphioxus genome are expressed is in sharp contrast to the situation in vertebrates, most of which contain only a single aanat.
the existence of multiple bfaanat genes seems to be the product of several duplications of this gene because there is no evidence of genome duplication events within the amphioxus line.
thus, the aanats of amphioxus represents a typical lineage-specific expansion of a paralogous gene family, a phenomenon that is widely observed in eukaryotes and is generally thought to reflect specific physiological adaptations  <cit> .
the lineage-specific expansion of aanats in amphioxus is the most extensive so far detected for this family but smaller expansions are seen in trichoplax, as demonstrated here, and in teleost fish.
the finding that the amino acid sequences of some of the paralogous amphioxus aanats are less than 40% identical suggests distinctly different substrate preferences, stability and kinetic characteristics.
based on these findings, it seems reasonable to hypothesize that each of the bfaanat paralogs plays a unique biological role.
this function is obviously not the time keeping by melatonin synthesis as it is in vertebrates, since there is no circadian rhythm in the blaanat activity; moreover, no melatonin has been found in amphioxus , and attempts to identify melatonin receptors in amphioxus were unsuccessful  <cit> .
one likely function of bfaanats is detoxification.
as discussed above, this possibility is supported by the observation that bfaanatα and bfaanatδ' acetylate diamines and alkylamines, in addition to arylalkylamines.
amphioxus lives buried in sediment, and, under these conditions, a wide range of aanat substrates might be generated as a result of biodegradation.
these molecules could be nonspecifically ingested or otherwise absorbed by the organism and as a result might represent a challenge to amphioxus physiology.
under these circumstances, the presence of multiple copies of aanat would provide a critical survival advantage.
considering the similarity between the life styles of amphioxus and trichoplax, it is tempting to speculate that these animals experience similar pressures from toxic environmental molecules, driving the fixation of gene duplications in the aanat family.
it is also possible that the functions of one or more bfaanats evolved from detoxification and ph regulation to a more specialized role in neurobiology.
this suggestion is based on the evidence that demonstrated substrates of bfaanatα and bfaanatδ', serotonin and octopamine, are thought to play a role in neurotransmission in amphioxus  <cit> .
the acetylation of these compounds by bfaanats is likely to inactivate them, thereby contributing to the process of neurotransmission by termination of neural stimulation.
transmitter inactivation is essential for efficient neural signaling.
in closing, it should be emphasized that the study presented here characterized only two of the seven bfaanat proteins.
accordingly, it is possible that other family members have distinctly different substrate preferences and unique biological roles.
the implications of day/night changes in aanat activity
we found a daily rhythm in total aanat activity, with high values during the day.
the finding that a daily rhythm in head aanat activity did not persist under constant light or constant darkness conditions is consistent with the interpretation that the aanat rhythm is not driven by a circadian clock, such as that found to operate in the amphioxus head  <cit> .
CONCLUSIONS
the unexpected finding of multiple bfaanat homolog genes points to the importance of this gene family in amphioxus.
the function of these genes appears to be detoxification of exogenous and endogenous amines, the latter including neurotransmitters.
it will be of interest to determine the precise distribution of each gene, how their expression is regulated and how each gene contributes to amphioxus biology.
abbreviations
aanat: arylalkylamine n-acetyltransferase; oaanat: ovine aanat; scaanat: saccharomyces cerevisiae aanat; bfaanat: b. floridae aanat; blaanat: b. lanceolatum aanat; accoa: acetyl coenzyme a; au: arbitrary unit; coa-hne: hydroxynaphthylethylamine coenzyme a; coa-t: tryptamine coenzyme a; dtt: dithiothreitol; gst: glutathione s-transferase; pbs: phosphate buffer saline; pea: phenylethylamine; tcep: trisphosphine.
authors' contributions
jp performed biochemical, physico-chemical, molecular biological studies and drafted the manuscript, ss participated in biochemical studies, lb participated in biological studies, slc participated in molecular genetic studies and helped to draft the manuscript, jlw participated in biochemical studies, gb, pg and jf participated in biochemical and biological studies, mvo and evk performed sequence and phylogenetic analysis and dck conceived of the study, participated in its coordination and helped to draft the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1
primers used for cloning pcr, rt-pcr localization of blaanat expression and qpcr experiments.
cloning primers were based on b. floridae genomic sequence and used to amplify the indicated blaanats.
click here for file
additional file 2
pairwise comparison scores of protein sequences.
pairwise comparison scores based on clustalw multiple sequence alignment of deduced peptide sequences of all  <dig> amphioxus aanat homologs and their known alleles.
shaded boxes designate comparisons between alleles.
the coordinates of each gene and allele  in genomic assembly v <dig>  are : bfaanatα, 54: <dig> ,444- <dig> ,187; bfaanatα', 56: <dig> - <dig> ; bfaanatβ, 54: <dig> ,188- <dig> ,679; bfaanatβ', 56: <dig> - <dig> ; bfaanatγ, 54: <dig> ,683- <dig> ,174; bfaanatγ', 56: <dig> - <dig> ; bfaanatδ, 54: <dig> ,399- <dig> ,929; bfaanatδ', 54: <dig> - <dig> ; bfaanatε, 200: <dig> ,854- <dig> , <dig> ; bfaanatζ, 200: <dig> ,115- <dig> ,588; bfaanatη, 234: <dig> - <dig>  ; bfaanatη', 81: <dig> ,130- <dig> , <dig> .
click here for file
additional file 3
positioning of the  <dig> genes for bfaanat on the assembly scaffolds as given in b. floridae genome assembly v <dig> .
click here for file
additional file 4
comparison of deduced amino acid sequences of characterized aanats to reference genomic sequences.
a. aanatα.
b. aanatδ'.
denotes derived from the published b. floridae genomic sequence.
denotes cloned from b. floridae head cdna, and are the full length recombinant proteins used for enzyme characterization.
denotes a partial fragment cloned from b. lanceolatum head cdna and used as a probe for in toto hybridization.
highlighted residues are those that differ from the published b. floridae genomic sequence. "
--", gap inserted to preserve alignment; "...", missing sequence.
click here for file
additional file 5
list of species used for construction of phylogenic tree.
list of species containing aanat homologs used for the construction of the phylogenetic tree shown in figure  <dig>  "gi" numbers precede the species name, where available; the source of the sequence is given in parentheses.
sequences for which "gi" numbers are not available are given in additional file  <dig>
click here for file
additional file 6
fasta sequences of aanat homologs.
fasta sequences of those of aanat homologs used for the construction of the phylogenetic tree shown in figure  <dig>  but not having "gi" numbers, .
click here for file
additional file 7
multiple alignment used for construction of the phylogenic tree.
multiple sequence alignment of the conserved regions of selected aanat sequences from various species used for the construction of the phylogenetic tree in figure  <dig>  species included in this alignment are given in additional file  <dig>
click here for file
additional file 8
average evolutionary distance  between major taxonomic groups for aanat proteins based on the truncated alignment shown in additional file <dig>
click here for file
additional file 9
average differences  between major taxonomic groups for aanat proteins, based on the truncated alignment shown in additional file <dig>  calculated using the mega program <cit> .
click here for file
additional file 10
average percent identity between major taxonomic groups for aanat proteins based on the truncated alignment shown in additional file <dig>  calculated from additional file <dig>
click here for file
additional file 11
model of constrained large-scale taxonomic topology tree of selected aanat sequences.
click here for file
additional file 12
topology tree - normal view.
reconstructed constrained large-scale taxonomic topology tree of selected aanat sequences - normal view based on model shown on additional file  <dig>  the scale bar represents the number of substitutions per position; numbers in parenthesis show number of species.
the number at internal branches shows the bootstrap support .
click here for file
additional file 13
topology tree - radiation view.
reconstructed constrained large-scale taxonomic topology tree of selected aanat sequences - radiation view based on model shown on additional file  <dig>  the scale bar represents the number of substitutions per position; numbers in parenthesis show number of species.
the number at internal branches shows the bootstrap support .
click here for file
additional file 14
comparison of day/night expression of three blaanats detected by qpcr in the head and body of amphioxus.
primers used are given in additional file  <dig>  asterisk indicates statistical significance of day/night differences.
data are presented as the mean ± s.e.m.
.
for further details see methods.
click here for file
additional file 15
inhibition of acetylation activity.
inhibition of acetylation activity in homogenate from amphioxus heads at ph  <dig>  using  <dig> mm pea and  <dig>  mm accoa as a substrate and various concentrations of coa-t as an inhibitor.
data are presented as the mean ± s.e.m.
.
click here for file
acknowledgements
this work was supported by the nih intramural research program through the nichd, national institutes of health.
the costs of publication of this article were defrayed in part by the payment of page charges.
this article must therefore be hereby marked "advertisement" in accordance with  <dig> u.s.c.
section  <dig> solely to indicate this fact.
work was also supported by the cnrs, upmc and the french national research agency .
we are grateful to dr. dan l. sackett for help with fluorescence measurement and drs.
howard jaffe and sam j. clokie for help with the maldi-tof analysis of protein samples.

 the aanat family is part of the large and diverse superfamily of gcn5-like acetyltransferases which use accoa as the acetyl donor and share a common accoa binding fold. members of the aanat family share limited sequence identity and are divided into two groups: vertebrate aanats; and, the non-vertebrate aanats. the latter are found in fungi, protists, and bacteria and lack defining physical characteristics of vertebrate aanat. the biological role of vertebrate aanat is to acetylate serotonin in the synthesis of melatonin. vertebrate aanat is associated with biological timing: daily changes in the activity of this enzyme regulate the daily rhythm in melatonin synthesis, which is essential for optimal temporal coordination of biological functions with night/day and seasonal changes and for photic entrainment. the central role of vertebrate aanat in biological timing has earned it the moniker'the timezyme'. non-vertebrate aanats are thought to play a detoxifying role by neutralizing arylalkylamines and a role in dna biology by acetylating polyamines. genomes of vertebrates contain a single copy of the aanat gene, except for teleost fish, some of which possess up to three paralogs, and cows which possess two paralogs. members of the aanat family are not in the available genomes of hemichordates and urochordates, which leaves open the question of when vertebrate aanat first appeared in chordates. a striking characteristic of vertebrate aanat is that it is consistently expressed at significant levels in only two tissues, both of which are photosensitive organs, the pineal gland and retina. this fits with the evidence that pinealocytes and retinal photoreceptors evolved from a common ancestral photodetector. the vertebrate aanat has a neutral ph optimum and exhibits high selectivity for arylalkylamines. vertebrate aanats encode proteins that have several highly conserved structural characteristics, which facilitate arylalkylamine acetylation and regulation. these features include flanking regulatory regions, which mediate rapid changes in enzyme activity; a pair of histidines which facilitate catalysis; and, a proline-containing tripeptide in a floppy loop, which confers a high catalytic rate through an effect on substrate binding. vertebrate aanats also have high selectivity for arylalkylamines conferred by the binding pocket. non-vertebrate type aanats are found in the genomes of most fungi, many unicellular eukaryotes, and a variety of bacteria. the proteins encoded by these genes do not contain the characteristic structural features of vertebrate aanats. genes similar to non-vertebrate aanats are not found in vertebrates. amphioxus aanats are of interest in understanding the evolution of vertebrate aanat because cephalochordates and vertebrates form a clade. accordingly, the study of the aanat from this organism might provide new insights into the transition from the non-vertebrate members of the aanat family to the vertebrate forms. sequence in the trace archives has indicated that a non-vertebrate type member of the aanat family is present in amphioxus. the subsequent availability of the near complete amphioxus genome has allowed us to comprehensively analyze the aanat family in amphioxus, as presented here. the findings of this study indicate that rapid expansion of non-vertebrate forms occurred in amphioxus and that a vertebrate form is absent. the substrate preference profile of the encoded enzymes suggests to us that they are likely to play a role in detoxification or in biogenic amine inactivation, or both. methods animals branchiostoma floridae and branchiostoma lanceolatum were used; the former were obtained from gulf specimen marine lab and the latter from the bay of banyuls sur mer. b. floridae heads were used to generate the cdna used to clone full length bfaanatα and bfaanatδ' that were used for enzyme characterization. b. lanceolatum were used for in vivo studies of gene expression and enzyme activity. they were cultured in sand-filtered seawater under a natural lighting cycle and temperature, as indicated elsewhere. intensity of the light was lux. for the experiment with constant conditions, animals were housed for days in constant light and/or darkness and then were used for the activity assay. identification of amphioxus aanats bfaanat homologs were identified through a tblastn search of genomic sequences available at: http://genome.jgi-psf.org. seven genes for aanat homologs for cloning of blaanat cdnas. qpcr primers for specific gene products were based on the b. lanceolatum clones. sequence similarity search and phylogenetic analysis the aanat homologs were identified by searching the non-redundant database of protein sequences using the psi blast program and by searching the genomes and proteomes available at jgi http://genome.jgi-psf.org/euk_cur html using the tblastn and blastp programs, respectively. additionally, aanat homologs were detected in the genomes of the red algae gracilaria changii, segmented worm alvinella pompejana, mediterranean mussel mytilus galloprovincialis and the basal fungus rhizopus oryzae http://www.broadinstitute.org/annotation/genome/rhizopus_oryzae/multihome.html. representatives from different groups including bacteria, fungi, protists, and animals were selected for phylogenetic analysis, and a multiple alignment of the respective protein sequences was constructed using the muscle program. a maximum likelihood phylogenetic tree was constructed using the treefinder program by optimizing a default starting tree constructed using the neighbor-joining method with the whelan and goldman empirical model of substitutions to select the best substitution model, we compared the likelihoods of trees constructed using models: wag, jtt, vt, blosum, dayhoff, cprev, rtrev, pmb by running the treefinder program; this test showed the best fit with the data for the wag model. the reliability of the internal tree branches was estimated using the lr-elw bootstrap method. a constrained large scale tree topology was created using the treeview program. alternative tree topologies were compared with the approximately unbiased test using the treefinder program. evolutionary distances were estimated using the distance matrix created by the protdist program of phylip package based on a truncated alignment of aanats and using the jones-taylor-thornton substitution model. the mean differences and identities between the major taxonomic groups of aanats were calculated using the mega program. prokaryotic expression of recombinant proteins pgex 4t- was used for the recombinant expression of full-length bfaanatα and bfaanatδ'. the constructed plasmid was sequenced to confirm identity. the construct, which produces a gst fusion protein, was transformed into e. coli strain bl21gold plyss. the cells were grown at 37°c; when the od =, the cultures were cooled to 25°c and isopropyl β-d-1-thiogalactopyranoside was added. the cells were harvested by centrifugation after hours of culturing and resuspended in 2× pbs, ph, containing mm dtt and a mixture of protease inhibitors. the cells were then lysed by sonication and the resulting lysate was centrifuged. the supernatant was mixed with glutathione-sepharose 4b affinity matrix pre-equilibrated with buffer a. the suspension was agitated for h and then packed into a glass column. the column was washed with column volumes of buffer a, followed by column volumes of buffer containing mm tris-hcl, m sodium citrate, mm dtt and 10% glycerol, ph. gst fusion protein was then eluted with column volumes of buffer b containing mm glutathione. protein was concentrated and dialyzed against buffer containing m ammonium acetate, mm nacl and mm tcep; the resulting preparation was stored at -80°c. where indicated, purified gst-free bfaanatα was used; otherwise, the gst-bfaanatα or gst-aanatδ' fusion product was used, consistent with previous studies. sds-page proteins were resolved on preformed 14% tris-glycine gels, using the manufacturer's protocol. the molecular mass of the proteins was estimated using rainbow™ markers. fluorescence-based protein studies fluorescence measurements were done on an iss pc photon-counting spectrofluorimeter at 22°c. to study binding properties of bfaanatα protein, the interaction between bfaanatα and coa-s-n-acetyl-7-hydroxynaphthylethylamine, a fluorescent bisubstrate inhibitor of vertebrate aanat, was analyzed. this was done by measuring coa-hne -dependent quenching of tryptophan fluorescence in m ammonium acetate, mm nacl and mm tcep, ph. off-rates of coa-hne were monitored as described using protein preparations that were 90% saturated with coa-hne. following a control incubation period, the probe was displaced with a 100-fold excess of another bisubstrate inhibitor, coa-s-n-acetyltryptamine, that is optically inactive in this analysis. the limiting value of koff was calculated by non-linear fitting using the equation: where a is the observed fluorescence anisotropy at time't' after the initial addition of coa-t and koff is the rate constant. determination of n-acetyltransferase activity recombinant bfaanatα activity was measured radiochemically by measuring acetyl product formed from acetyl-coenzyme a and an amine substrate as described; or, colorimetrically by measuring coash formed by incubation with accoa and an amine substrate. coash is detected as the colored product formed by reaction with'-dithio-bis. in the colorimetric assay, standard reactions were performed in a total volume of μl containing: the substrate at indicated concentrations, mm accoa, mg/ml bovine serum albumin, and mm ethylene diamine tetraacetic acid. for analysis of gst-bfaanatα, μg of protein was included in m tris buffer ph; the activity of ovine gst-aanat, which was used as a reference standard, was determined by incubating ng of protein in m phosphate buffer, ph. incubation conditions were either min at 37°c for amphioxus or min at 37°c for ovine aanat. enzymes were then inactivated by the addition of μl stop solution; absorbance was measured after min incubation at room temperature. the vmax and km were calculated by non-linear fitting using the equation: where v is the observed enzyme velocity, vmax is the maximum enzyme velocity, s is the substrate concentration, and km is the michaelis-menten constant. for the measurement of enzyme activity in amphioxus head tissue, the colorimetric assay was used with mm β-phenylethylamine as a substrate and μg of protein from tissue homogenized in assay buffer. for each measurement a corresponding assay was performed in the presence of mm coa-t, and another in the absence of pea. the controls used included an assay without tissue homogenate and another with a heat-inactivated homogenate. samples were collected at noon and midnight in the experiments studying the light/dark variations, and at noon in all other experiments. quantitative measurement of b. lanceolatum aanat transcripts total rna was extracted from head and from body using a trizol extraction kit. μg was incubated with unit of dnase i for min at 37°c. following dnase inactivation, the sequence was reverse-transcribed using powerscript reverse transcriptase. qpcr of the resulting cdnas was performed with the primer sets qaanat using the'light cycler fast start dna master sybr green i kit' in a total volume of μl as follows: 95°c, followed by cycles of denaturation at 95°c, annealing at 60°c and extension at 72°c. qpcr was done using a light cycler. amplification efficiency was measured using serial dilutions of plasmids containing a fragment of the corresponding blaanat gene. data were normalized to the average value obtained with three housekeeping genes from b. lanceolatum: l, l and actin. the pcr conditions for the housekeeping genes were as indicated above except that the annealing temperature was 66°c for actin. primers for the housekeeping genes are given in additional file tissue distribution of amphioxus aanatα studied by rt-pcr the tissue distribution of amphioxus aanatα was determined by rt-pcr. total rna from different parts of b. lanceolatum was extracted using the trizol method; μg was incubated with unit of dnase i for min at 37°c. following dnase inactivation, the rna was reverse-transcribed using powerscript reverse transcriptase pcr amplifications of the resulting cdnas using advantage cdna polymerase were performed with the primers for blaanatα, and the following conditions: 95°c followed by cycles of denaturation at 94°c, annealing at 66°c and extension at 68°c. in the controls, the template cdna was replaced by either water or rna that was not reverse-transcribed. the pcr products were loaded in a 1% agarose gel, in the presence of dna/hinf i marker. fragments of the expected size were extracted, subcloned in pgem-t easy and sequenced. tissue distribution of amphioxus aanatα studied by in toto hybridization b. lanceolatum were collected at argelès-sur-mer in july 2009; gametes were obtained by heat stimulation. eggs and embryos were fixed in 4% paraformaldehyde in mops-egta buffer and processed for in toto hybridization as detailed elsewhere, except that the chromogenic reaction was performed using bm purple. two digoxigenin labeled probes corresponding to the α or δζ subtypes were designed. the former was obtained using the cloning primers caanatαf and caanatαr and yielded a probe of bp. the latter was obtained using the cloning primers caanatδζf and caanatδζr and yielded a probe of bp. hybridization treatments with the different sense and antisense probes were performed simultaneously.  sequence analysis, phyletic distribution, and phylogeny of the aanats seven intronless aanat genes were identified in the b. floridae genome. a second allele was found for of the genes in version v of the b. floridae genome assembly. the amino acid sequence identity among the amphioxus aanats ranged from to 86% excluding allele/allele comparisons. the structural features which characterize vertebrate aanats, including the flanking regulatory regions, a pair of neighboring histidines in the active site, and a tripeptide in the floppy loop, do not appear in the proteins encoded by the amphioxus aanats. the motifs a and b, which characterize members of the gcn5-like superfamily of acetyltransferases are conserved in both amphioxus and vertebrate aanats. the relative locations of the bfaanats on assembly scaffolds are shown schematically in additional file the genes comprise groups based on sequence similarity - group 1: αη group 2: β group 3: γε group 4: δζ. primers were prepared to amplify members of each of the four groups in b. lanceolatum; low sequence conservation prevented amplification of all seven. the resulting pcr products were sequenced and designated blaanatα, blaanatα' blaanatγ, and blaanatγδ, according to the homology with individual bfaanat genes; representatives of the other genes were not obtained in this screen. the clone blaanatδζ was designated this way because it could not be reliably ascribed to either the δ or ζ gene based on similarity to the b. floridae counterparts. the amino acid sequence identity between the cloned pcr products for the b. lanceolatum α, α', γ and δζ isoforms and the corresponding b. floridae genes was 94%, 86%, 76%, and 73%, respectively. comparisons of the amino acid sequences of the cloned blaanatα and blaanatδζ fragments to their b. floridae counterparts are shown in additional file examination of the phyletic distribution of the detected aanat family members shows that the family is represented not only in chordates and cephalochordates but also in some invertebrates such as the primitive metazoan trichoplax, the segmented worms capitella and alvinella, and the mollusk mytilus. among non-metazoan eukaryotes, aanat homologs were detected in different fungal lineages including basidiomycetes, ascomycetes, and the basal fungi rhizopus oryzae and phycomyces blakesleeanus, and in many unicellular eukaryotes. the protists in which aanat homologs were detected include representatives from three supergroups of eukaryotes, namely, plantae, chromalveolates and excavates. in plantae, aanat homologs were found in green and red algae; among chromalveolates, members of the aanat family were present in diatoms, oomycetes, and a haptophyte emiliania huxleyi; and, among excavates, an aanat homolog was identified in the heterolobosean naegleria gruberi. aanat homologs are also present in a number of bacterial lineages including firmicutes, actinobacteria, proteobacteria and the bacteroidetes/chlorobi group. aanats are relatively short proteins which complicates reliable phylogenetic analysis. phylogenetic trees were constructed based on an alignment of amino acid sequences. eight treefinder runs were performed for substitution models. the ml tree constructed with the wag matrix was associated with the highest likelihood value. the topology of this tree is incompatible with the species phylogeny. the monophyly of all amphioxus aanats was moderately supported and these sequences formed an apparent clade with mollusk and segmented worm that was even more strongly supported. however, beyond this invertebrate branch, there were many anomalies in the tree. in particular, vertebrates clustered with ascomycetes and bacteria, whereas the position of trichoplax remained uncertain. the monophyly of fungi was not recovered, and ostreococii were separated from another green alga represented in the selected set of protists. we tested the statistical compatibility of the aanat tree with the large-scale organismal phylogeny of eukaryotes using the constrained tree approach. a constrained tree was created in which all species were combined into four major clades in accord with the overall organismal phylogeny, namely, bacteria, protists, fungi and animals. inside the animal branch, with trichoplax at the root, the amphioxus branch was placed as the sister group of vertebrates, and the segmented worms and mussel were positioned between chordates and trichoplax. the fungal branch was placed as the sister group of animals. the approximately unbiased test result for the tree generated with this approach shows that the "classical" constrained tree topology could not be rejected with confidence although the likelihood of that topology is much lower than that of the best tree shown in figure the results of the phylogenetic analysis of the aanats emphasize the major distinction between the vertebrate forms and those from other animals, presumably, owing to a dramatic acceleration of evolution at the base of the vertebrate branch. beyond that conclusion, however, the results were quite uncertain due to the low information content of the alignment, between-branch differences in the evolutionary rates, and uneven sampling of different groups. expression and characterization of bfaanatα protein a comparison of the predicted amino acid sequences derived from the published genomic sequence and from our cloned bfaanatα revealed differences in five amino acid residues scattered throughout. these differences, most likely, reflect the high degree of polymorphism in b. floridae. the effects of these differences are unknown. to determine the molecular mass of expressed bfaanatα, the gst-bfaanatα fusion protein was cleaved and bfaanatα was purified as described in the methods. the final preparation contained a single major band of protein which migrated at ~ kda on sds-page. maldi-tof analysis revealed the mass to be corresponding to the predicted mass based on amino acid composition. using the same procedure we also expressed gst-bfaanatδ'. to establish whether the expressed protein was correctly folded, we measured binding of a fluorescent bisubstrate inhibitor of the expressed protein, coa-hne; binding of the inhibitor to oaanat quenches fluorescence. here we found that bfaanatα quenched inhibitor fluorescence and that the expressed protein bound the inhibitor stoichiometrically. moreover, binding was totally disrupted by treatment with m guanidine. accordingly, the purified preparation of bfaanatα is not denatured and possesses characteristics generally similar to those of oaanat, as regards binding of coa-hne. to further characterize bfaanatα, off-rate analysis was done using the displacement strategy described previously, in which bound coa-hne is displaced by coa-t, a non-fluorescent aanat inhibitor that does not interfere with detection of coa-hne. this result shows that the off-rate observed for the coa-hne is > s- similar to that of yeast aanat; these rates are much faster than that for oaanat, consistent with the view that the proline-containing tripeptide in the floppy loop of vertebrate aanat enhances binding of substrate. biochemical characterization of bfaanatα and bfaanatδ' activity as a function of ph: gst-bfaanatα and gst-bfaanatδ' both exhibited a ph activity profile that differed from that of ovine aanat and yeast aanat as examined using mm tryptamine as a substrate. the ph optimum for gst-bfaanatα was determined to be, which differs substantially from published data on fish, ovine and yeast aanat. the ph optimum for gst-bfaanatδ' seems to be even higher. similar results were obtained using butylamine as a substrate. effect of aanat bisubstrate inhibitor: enzyme activity was inhibited by the bisubstrate inhibitor coa-t. the ic value was found to be ~10-fold greater than that reported for ovine aanat. temperature dependence: the temperature optimum for bfaanatα was found to be 36°c. this is similar to that of ovine aanat, yeast aanat or pike aanat-. substrate preference: the relative capacities of oaanat, bfaanatα and bfaanatδ' to acetylate a diamine, alkylamines and arylalkylamines were determined at ph for oaanat, and for bfaanatα and bfaanatδ' using gst-fusion proteins. the activity of bfaanats towards the substrates tested was relatively similar and fell within a 5-fold range of activity; the highest activity was observed with pea, however, bfaanatα and bfaanatδ' seemed to show a slight preference for hydrophobic or polar substrates, respectively. the km and vmax values were determined for five biogenic arylalkylamines: octopamine, pea, tyramine, tryptamine and serotonin. activity was measured at ph. values are calculated from data in figure vmax values are nmol/h/pmol enzyme and are given as the mean ± s.e.m.. for further details see methods. in contrast to the relatively broad substrate preference profiles of bfaanats observed with the substrates tested, oaanat exhibited markedly higher activity towards the two indolamines, tryptamine and serotonin, and lower activity towards putrescine, propylamine and butylamine. this points to a marked difference in substrate preference: bfaanatα and bfaanatδ' exhibit a broad capacity to acetylate amines, including arylalkylamines and a diamine, whereas oaanat exhibits a distinct preference for arylalkylamines with an aromatic group, including serotonin, pea and tryptamine. control of the abundance of blaanat transcripts in some vertebrate species, the nocturnal increase in aanat activity is accompanied by a ~100- fold increase in the abundance of transcripts. here we used qpcr to determine whether the abundance of blaanat mrnas varied on a night/day basis. however statistically significant differences were only observed for blaanatα. accordingly, although it does not appear that there are marked night/day differences in expression of blaanats, it is clear that small differences do occur. analysis of tissue expression by rt-pcr investigation of tissue expression of blaanatα by rt-pcr revealed the presence of mrna along the whole body. the identity of the product was confirmed by sequencing. daily changes in blaanat activity homogenates of the heads of b. lanceolatum were used for measurement of aanat activity, using mm pea as a substrate; coa-t was added to determine the contribution from aanat to total acetyltransferase activity. in all conditions tested, activity was significantly higher in heads sampled at midday than in those sampled at midnight. activity measured in the absence of added substrate was more than 50% as high as that measured in the presence of exogenous substrate; in both cases, activity was nearly completely inhibited by coa-t, indicating that blaanat acetylated endogenous substrates during the assay. in other studies this rhythm was confirmed, however a similar daily change was not observed in animals housed in constant darkness or constant lighting. a detailed graph describing the inhibition of acetylation activity in homogenates by coa-t is shown in additional file substrates comparison of samples collected during day and night. the exogenous substrate is mm pea; mm coa-t has been used for inhibition. data are given as the mean ± s.e.m.. acetylation activity under normal light/darkness changes compared with constant light and constant darkness conditions. the ph of the assay was; data are given as the mean ± s.e.m.. in toto hybridization expression of one or both blaanats was detected by in toto hybridization at all stages examined, in varying patterns. although it is difficult to quantify the colorimetric reactions, it seemed that the α form was more highly expressed than the δζ form at early stages, whereas the opposite seemed true at later stages. the first specific hybridization signal was detected at the 8-cell embryo and was stronger in out of the cells. at the blastula stage most of the cells exhibited a diffuse labeling with both probes; however, a group of cells exhibited a more intense dark signal. at the early neurula stage, the labeling appeared diffuse and weak with the δζ probe, whereas some regionalization appeared with the α probe. this tendency was even more apparent at the late neurula stage where the α probe gave a specific labeling in the antero-dorsal and posterio-dorsal parts of the embryo, as well as antero-ventrally in the general area of the intestine. at the premouth stage, no labeling with the α probe was detected, whereas the δζ probe hybridized from head to tail in a dorsal band, corresponding to the neural tube. discussion the publication of the near complete amphioxus genome has marked a new era in understanding chordate evolution. the studies presented in this report, which were made possible by this advance, focus on one family of genes and extends our understanding of this gene family in amphioxus and, potentially, in animals in general. aanats in the amphioxus genome and evolution of the aanat family analysis of the amphioxus genome indicates that it lacks a vertebrate-like aanat. in contrast, our analysis shows that the amphioxus genome contains a set of closely related aanat homologs that are similar in their gene and protein architectures to homologs from invertebrates, protists, fungi, and bacteria. none of these non-vertebrate homologs possesses the defining structural features of vertebrate aanats. the discovery of the family of aanats in amphioxus sheds light on the timing of the emergence of the typical vertebrate aanats that, most likely, occurred subsequent to the radiation of the vertebrate lineage from the common ancestor with cephalochordates. it appears that the change in the aanat gene architecture at the base of the vertebrate tree was accompanied by a substantial acceleration of the evolutionary rate which is the likely reason why the amphioxus aanat homologs are more similar to those from fungi and protists than to those from vertebrates. the results of phylogenetic analysis of the aanat family presented here were not entirely conclusive owing, in part, to the relatively small number of phylogenetically informative positions in the multiple sequence alignment of the aanat family, and in part, to the non-uniformity of evolutionary rates, and perhaps, to a rapid divergence associated with evolutionary transitions. the best ml tree showed multiple deviations from the species tree but the topology corresponding to the species tree, including the cephalochordate-vertebrate clade, could not be rigorously rejected. thus, the phylogeny of aanats might follow the organismal phylogeny of animals, taking into account the apparent multiple losses, e.g., in arthropods, nematodes, and echinodermata, as well as accelerations of evolution, particularly, in vertebrates. previous analyses of the phyletic patterns and phylogenies of the aanats led to the suggestion that horizontal gene transfer, perhaps, even into the animal lineage or the animal-fungal ancestral lineage was a major contribution to the evolution of this family. that possibility still remains considering the patchy phyletic distribution of the aanats and the irregularities in the topologies of the phylogenetic trees, such as the polyphyly of fungi and the unusual assortment of protists. however, given the lack of strong resolution in the trees, the major changes in evolutionary rates like the one in vertebrates, and the extreme propensity of the aanat family for gene loss, scenarios that involve unusual routes of horizontal transfer should be addressed with much caution. biochemical features of bfaanatα and bfaanatδ' the biochemical characteristics of bfaanatα and bfaanatδ' were analyzed in this study. both appear to have significantly lower intrinsic activity at ph as compared to oaanat. the substrate preference spectrum of the bfaanats is broader than that of oaanat. as with other non-vertebrate aanats, these differences are likely to reflect differences in the primary structure of the enzyme, including the presence in vertebrate aanat of two histidines near the site of acetyl transfer. these are thought to be part of a proton wire which accelerates catalysis by facilitating the expulsion of protons generated during acetyl transfer, thereby preventing their accumulation and the resulting acidification of the active site. another important difference in the primary structure between vertebrate and amphioxus aanats occurs in a floppy loop which is critical for catalysis, apparently by increasing the dwell time of amine stubstrates in the binding pocket. the floppy loop of the vertebrate aanats is three residues longer than that of known aanat homologs, including amphioxus. this addition appears to enhance catalysis by directing the movement of the floppy loop and preventing intramolecular interactions which inhibit binding. as a result, the off-rate of substrate is substantially increased when these three residues are experimentally deleted. during the course of evolution these three residues were added and, as shown here, this was associated with a decrease in the off-rate relative to that of amphioxus and yeast aanats. regional expression of blaanat our studies revealed that blaanatα is expressed throughout the body, in contrast to the pattern of expression of vertebrate aanat, which is primarily limited to structures located in the head: pineal gland, retina and in some cases the pituitary gland. the development of these structures is strongly dependent on expression of the transcription factor pax pax expression in amphioxus is limited to the frontal organ and lamellar organ; the latter has been identified as a possible pineal precursor, based on its medial location and presence of photoreceptors in the tissue. in toto hybridization the in toto hybridization analyses revealed expression of both α and δζ genes early in development in b. lanceolatum. it seems that the relative strength of expression changes during development. blaanatδζ was found to be expressed in a region close to the neural tube. this observation and the finding that these bfaanats acetylate the biogenic amines serotonin and octopamine, suggest to us that it is reasonable to consider that these enzymes play a role in the metabolism of these neurotransmitters. the possible roles of amphioxus aanats the unexpected finding that several of the seven aanats in the amphioxus genome are expressed is in sharp contrast to the situation in vertebrates, most of which contain only a single aanat. the existence of multiple bfaanat genes seems to be the product of several duplications of this gene because there is no evidence of genome duplication events within the amphioxus line. thus, the aanats of amphioxus represents a typical lineage-specific expansion of a paralogous gene family, a phenomenon that is widely observed in eukaryotes and is generally thought to reflect specific physiological adaptations. the lineage-specific expansion of aanats in amphioxus is the most extensive so far detected for this family but smaller expansions are seen in trichoplax, as demonstrated here, and in teleost fish. the finding that the amino acid sequences of some of the paralogous amphioxus aanats are less than 40% identical suggests distinctly different substrate preferences, stability and kinetic characteristics. based on these findings, it seems reasonable to hypothesize that each of the bfaanat paralogs plays a unique biological role. this function is obviously not the time keeping by melatonin synthesis as it is in vertebrates, since there is no circadian rhythm in the blaanat activity; moreover, no melatonin has been found in amphioxus, and attempts to identify melatonin receptors in amphioxus were unsuccessful. one likely function of bfaanats is detoxification. as discussed above, this possibility is supported by the observation that bfaanatα and bfaanatδ' acetylate diamines and alkylamines, in addition to arylalkylamines. amphioxus lives buried in sediment, and, under these conditions, a wide range of aanat substrates might be generated as a result of biodegradation. these molecules could be nonspecifically ingested or otherwise absorbed by the organism and as a result might represent a challenge to amphioxus physiology. under these circumstances, the presence of multiple copies of aanat would provide a critical survival advantage. considering the similarity between the life styles of amphioxus and trichoplax, it is tempting to speculate that these animals experience similar pressures from toxic environmental molecules, driving the fixation of gene duplications in the aanat family. it is also possible that the functions of one or more bfaanats evolved from detoxification and ph regulation to a more specialized role in neurobiology. this suggestion is based on the evidence that demonstrated substrates of bfaanatα and bfaanatδ', serotonin and octopamine, are thought to play a role in neurotransmission in amphioxus. the acetylation of these compounds by bfaanats is likely to inactivate them, thereby contributing to the process of neurotransmission by termination of neural stimulation. transmitter inactivation is essential for efficient neural signaling. in closing, it should be emphasized that the study presented here characterized only two of the seven bfaanat proteins. accordingly, it is possible that other family members have distinctly different substrate preferences and unique biological roles. the implications of day/night changes in aanat activity we found a daily rhythm in total aanat activity, with high values during the day. the finding that a daily rhythm in head aanat activity did not persist under constant light or constant darkness conditions is consistent with the interpretation that the aanat rhythm is not driven by a circadian clock, such as that found to operate in the amphioxus head.  the unexpected finding of multiple bfaanat homolog genes points to the importance of this gene family in amphioxus. the function of these genes appears to be detoxification of exogenous and endogenous amines, the latter including neurotransmitters. it will be of interest to determine the precise distribution of each gene, how their expression is regulated and how each gene contributes to amphioxus biology. abbreviations aanat: arylalkylamine n-acetyltransferase; oaanat: ovine aanat; scaanat: saccharomyces cerevisiae aanat; bfaanat: b. floridae aanat; blaanat: b. lanceolatum aanat; accoa: acetyl coenzyme a; au: arbitrary unit; coa-hne: hydroxynaphthylethylamine coenzyme a; coa-t: tryptamine coenzyme a; dtt: dithiothreitol; gst: glutathione s-transferase; pbs: phosphate buffer saline; pea: phenylethylamine; tcep: trisphosphine. authors' contributions jp performed biochemical, physico-chemical, molecular biological studies and drafted the manuscript, ss participated in biochemical studies, lb participated in biological studies, slc participated in molecular genetic studies and helped to draft the manuscript, jlw participated in biochemical studies, gb, pg and jf participated in biochemical and biological studies, mvo and evk performed sequence and phylogenetic analysis and dck conceived of the study, participated in its coordination and helped to draft the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1 primers used for cloning pcr, rt-pcr localization of blaanat expression and qpcr experiments. cloning primers were based on b. floridae genomic sequence and used to amplify the indicated blaanats. click here for file additional file 2 pairwise comparison scores of protein sequences. pairwise comparison scores based on clustalw multiple sequence alignment of deduced peptide sequences of all amphioxus aanat homologs and their known alleles. shaded boxes designate comparisons between alleles. the coordinates of each gene and allele in genomic assembly v are: bfaanatα, 54:,444-,187; bfaanatα', 56: -; bfaanatβ, 54:,188-,679; bfaanatβ', 56: -; bfaanatγ, 54:,683-,174; bfaanatγ', 56: -; bfaanatδ, 54:,399-,929; bfaanatδ', 54: -; bfaanatε, 200:,854-,; bfaanatζ, 200:,115-,588; bfaanatη, 234: -; bfaanatη', 81:,130-,. click here for file additional file 3 positioning of the genes for bfaanat on the assembly scaffolds as given in b. floridae genome assembly v. click here for file additional file 4 comparison of deduced amino acid sequences of characterized aanats to reference genomic sequences. a. aanatα. b. aanatδ'. denotes derived from the published b. floridae genomic sequence. denotes cloned from b. floridae head cdna, and are the full length recombinant proteins used for enzyme characterization. denotes a partial fragment cloned from b. lanceolatum head cdna and used as a probe for in toto hybridization. highlighted residues are those that differ from the published b. floridae genomic sequence. "--", gap inserted to preserve alignment; "...", missing sequence. click here for file additional file 5 list of species used for construction of phylogenic tree. list of species containing aanat homologs used for the construction of the phylogenetic tree shown in figure "gi" numbers precede the species name, where available; the source of the sequence is given in parentheses. sequences for which "gi" numbers are not available are given in additional file click here for file additional file 6 fasta sequences of aanat homologs. fasta sequences of those of aanat homologs used for the construction of the phylogenetic tree shown in figure but not having "gi" numbers,. click here for file additional file 7 multiple alignment used for construction of the phylogenic tree. multiple sequence alignment of the conserved regions of selected aanat sequences from various species used for the construction of the phylogenetic tree in figure species included in this alignment are given in additional file click here for file additional file 8 average evolutionary distance between major taxonomic groups for aanat proteins based on the truncated alignment shown in additional file click here for file additional file 9 average differences between major taxonomic groups for aanat proteins, based on the truncated alignment shown in additional file calculated using the mega program. click here for file additional file 10 average percent identity between major taxonomic groups for aanat proteins based on the truncated alignment shown in additional file calculated from additional file click here for file additional file 11 model of constrained large-scale taxonomic topology tree of selected aanat sequences. click here for file additional file 12 topology tree - normal view. reconstructed constrained large-scale taxonomic topology tree of selected aanat sequences - normal view based on model shown on additional file the scale bar represents the number of substitutions per position; numbers in parenthesis show number of species. the number at internal branches shows the bootstrap support. click here for file additional file 13 topology tree - radiation view. reconstructed constrained large-scale taxonomic topology tree of selected aanat sequences - radiation view based on model shown on additional file the scale bar represents the number of substitutions per position; numbers in parenthesis show number of species. the number at internal branches shows the bootstrap support. click here for file additional file 14 comparison of day/night expression of three blaanats detected by qpcr in the head and body of amphioxus. primers used are given in additional file asterisk indicates statistical significance of day/night differences. data are presented as the mean ± s.e.m.. for further details see methods. click here for file additional file 15 inhibition of acetylation activity. inhibition of acetylation activity in homogenate from amphioxus heads at ph using mm pea and mm accoa as a substrate and various concentrations of coa-t as an inhibitor. data are presented as the mean ± s.e.m.. click here for file acknowledgements this work was supported by the nih intramural research program through the nichd, national institutes of health. the costs of publication of this article were defrayed in part by the payment of page charges. this article must therefore be hereby marked "advertisement" in accordance with u.s.c. section solely to indicate this fact. work was also supported by the cnrs, upmc and the french national research agency. we are grateful to dr. dan l. sackett for help with fluorescence measurement and drs. howard jaffe and sam j. clokie for help with the maldi-tof analysis of protein samples.
expression of vertebrate aanats is limited primarily to the pineal gland and retina, where it plays a role in controlling the circadian rhythm in melatonin synthesis. non-vertebrate aanats, which occur in fungi and protists, are thought to play a role in detoxification and are not known to be associated with a specific tissue. bfaanatα and bfaanatδ' have unusually alkaline optimal ph, more than two ph units higher than that of vertebrate aanats. analysis of the expression pattern of the two most structurally divergent aanats in branchiostoma lanceolatum revealed that they are expressed early in development and also in the adult at low levels throughout the body, possibly associated with the neural tube. the substrate selectivity profiles of bfaanatα and δ' are relatively broad, including alkylamines, arylalkylamines and diamines, in contrast to vertebrate forms, which selectively acetylate serotonin and other arylalkylamines. based on the role melatonin plays in biological timing, aanat has been given the moniker "the timezyme". based on these features, it appears that amphioxus aanats could play several roles, including detoxification and biogenic amine inactivation. blaanat activity is influenced by environmental lighting, but light/dark differences do not persist under constant light or constant dark conditions, indicating they are not circadian in nature. we have found that the amphioxus genome contains seven aanats, all having non-vertebrate type features. expression is clearly not exclusively associated with the proposed analogs of the pineal gland and retina.
 the arylalkylamine n-acetyltransferase family is divided into structurally distinct vertebrate and non-vertebrate groups. expression of vertebrate aanats is limited primarily to the pineal gland and retina, where it plays a role in controlling the circadian rhythm in melatonin synthesis. based on the role melatonin plays in biological timing, aanat has been given the moniker "the timezyme". non-vertebrate aanats, which occur in fungi and protists, are thought to play a role in detoxification and are not known to be associated with a specific tissue.  we have found that the amphioxus genome contains seven aanats, all having non-vertebrate type features. this and the absence of aanats from the genomes of hemichordates and urochordates support the view that a major transition in the evolution of the aanats may have occurred at the onset of vertebrate evolution. analysis of the expression pattern of the two most structurally divergent aanats in branchiostoma lanceolatum revealed that they are expressed early in development and also in the adult at low levels throughout the body, possibly associated with the neural tube. expression is clearly not exclusively associated with the proposed analogs of the pineal gland and retina. blaanat activity is influenced by environmental lighting, but light/dark differences do not persist under constant light or constant dark conditions, indicating they are not circadian in nature. bfaanatα and bfaanatδ' have unusually alkaline optimal ph, more than two ph units higher than that of vertebrate aanats.  the substrate selectivity profiles of bfaanatα and δ' are relatively broad, including alkylamines, arylalkylamines and diamines, in contrast to vertebrate forms, which selectively acetylate serotonin and other arylalkylamines. based on these features, it appears that amphioxus aanats could play several roles, including detoxification and biogenic amine inactivation. the presence of seven aanats in amphioxus genome supports the view that arylalkylamine and polyamine acetylation is important to the biology of this organism and that these genes evolved in response to specific pressures related to requirements for amine acetylation. 