BACKGROUND
synaptotagmin  proteins participate in regulated membrane fusion events in multicellular organisms.
syt research dates back to  <dig>  when the first syt was identified as an integral protein of synaptic vesicles  <cit> .
since then, much effort has gone into discovering the function and the detailed mechanism of action of this protein, syt <dig>  numerous experimental approaches have shown that syt <dig> is crucially involved in fast neurotransmitter release at synapses .
syt <dig> binds calcium, serving as the calcium sensor which triggers synaptic vesicle exocytosis.
syt <dig> also binds to the neuronal snare proteins which are required for membrane fusion, as well as to membranes directly.
the primary structure of syt <dig>  <cit>  revealed three important features: an n-terminal transmembrane  domain, which serves to anchor the protein in the synaptic vesicle, plus two, tandem, c-terminal, cytoplasmic c <dig> domains  which specify the calcium, snare and membrane binding properties.
crystallographic studies of some syt c <dig> domains  <cit>  have revealed the exact nature of their calcium binding abilities.
structural bioinformatics  <cit>  has shown that c <dig> domains are very common in the eukaryotic protein repertoire, but not all act to bind calcium.
c <dig> domains can occur singly or as multiple copies in a given protein  <cit> .
the particular domain organisation of syt <dig> is important for its role in synaptic vesicle exocytosis.
the tm domain tethers the protein to the membrane.
the c2a and c2b domains, endowed with individual capacities to bind calcium, snare proteins, other accessory proteins and phospholipids, act both independently and synergistically, at different stages during the life cycle of the synaptic vesicle, to promote or inhibit fusion  <cit> .
after the genetic code for syt <dig> was discovered  <cit>  genes similar to syt <dig> were sought and found.
biochemical methods which were employed initially, led to an increase in the size of the rodent syt family, from one member to thirteen members and beyond.
syt relatives were also identified and studied in d. melanogaster and c. elegans.
all of these homologous syt genes were found to encode proteins with a common domain architecture: an n-terminal tm domain connected by a variable length, poorly conserved linker sequence, to well conserved, tandem, c-terminal c2a and c2b domains and this domain architecture was taken to define the family.
functional studies with the proteins in vivo and in vitro, indicate that some respond to calcium but others do not.
whether or not they respond to calcium, syt proteins are able to regulate membrane fusion due to common properties of their c <dig> domains  <cit> .
studies of the anatomical expression patterns of syt genes show that syt <dig> is abundantly expressed in nervous tissue  <cit> .
in d. melanogaster, syt <dig> expression is neuron specific but other syt genes are expressed elsewhere  <cit> .
in rodents, most syt genes are expressed in the brain  <cit> .
with the publication of whole genome sequences from multicellular organisms, it became possible to identify and compare complete genome complements of syt genes using computational methods  <cit> .
the phylogeny of these genes could be examined by comparing sequences from different organisms.
seventeen syt homologues were found in each of the complete h. sapiens and m. musculus genome sequences.
these were named syt <dig>  to syt <dig>  according to the nomenclature conventions of the hugo gene nomenclature committee  <cit>  .
some of these genes however, encode syt proteins which lack the requisite domain architecture for inclusion in a syt family defined by domain structure.
the complete genome sequences of c. elegans, d. melanogaster and a. thaliana were analysed, together with draft genome sequences from other eukaryotes, to assess their syt coding potential  <cit> .
proteins which share the stereotypical syt domain architecture, are encoded by distinct  gene families in plants and in animals  <cit>  further confounding the notion of a homologous syt family defined in terms of domain architecture.
functional similarity between plant and animal syt proteins, at the level of calcium and membrane binding, has been demonstrated  <cit>  in keeping with the notion that protein form specifies function.
proteins have been classified as syt family members by the presence or absence of suitable characteristics in terms of protein form and function  without studious regard to gene ancestry.
legitimately but confusingly, family membership can be defined in different ways: in terms of descent from a common ancestor  or in terms of the common possession of distinguishing features.
since it is not straightforward to demonstrate that homology is the cause of the similarity between biological entities such as proteins, it is often just assumed, leading to the confusion of two different notions of a protein family.
please see  <cit>  for explanations of the terminology of homology.
examination of the genomic specification of eukaryotic proteins can help to reveal their phylogenetic relationships because the evolutionary conservation of gene structure provides additional information on which to base an inference of homology.
families may thereby be resolved as those related by homology and those otherwise related.
ambiguity still remains however, because depending on the exact biological context, members of a family, however defined, can sometimes be considered as functionally interchangeable and thus worthy of a common functional identity.
thus, genes and proteins can acquire multiple identities depending on the perspective from which they are considered.
universal agreement about nomenclature for genes and proteins which is suitably clear and meaningful is a challenge for the future, but norms for gene nomenclature according to phylogenetic relationships are established  <cit> .
it may be the case  that members of a gene family related by heredity, share a common functional role  <cit> .
however, due to the dynamics of eukaryotic genome evolution  <cit>  the functional attributes of gene relatives can diverge.
it is also possible, that separate start points, through advantageous gain of function, can lead through selection, to convergence upon a common form and function.
while it is practicable to assign family membership to genes or proteins based on their apparent hereditary relationships, the functional implications of family membership cannot simply be deduced.
more than a decade of research effort has been aimed at understanding how syt <dig> works.
considerable effort has also been directed at understanding the functions of other syt proteins.
current understanding however, is insufficient to predict the functions of a whole genome complement of syt proteins, or indeed to predict what a genome complement of syt proteins might consist of.
to help improve this lack of understanding, i have put together an open access resource intended to serve future research aimed at understanding the biological meaning of the hereditary and functional relationships among some of the metazoan syt-like proteins.
i have collected and compiled manually, information about homologous syt, esyt, rph3a, rph3al, doc <dig> and dblc <dig> genes from  <dig> metazoan genome sequences spanning a wide evolutionary range.
i have not attempted to collect all genes capable of encoding proteins with similarity on any level to syt proteins.
such a collection would ultimately include perhaps, all eukaryotic c <dig> domain proteins.
my collection is restricted to homologous syt genes plus a small group of genes already identified as similar but not homologous to syt genes  <cit> .
the collection does include a new group of homologous genes  which are present in the genomes of marine invertebrates and encode proteins similar to syt proteins.
the information in this collection may interest investigators in the field of evolutionary biology, with examples of dynamic genome evolution including whole gene duplication, partial duplication, gene fission, acquisition of novel coding exons, gene inactivation, extreme sequence divergence, intron mobilization and a variety of routes to the expression of altered gene products.
for experimentalists who want to understand how syt proteins act to effect the biology of animal nervous systems, the collection provides a utility with great evolutionary depth.
in order to make progress in neuroscience, appropriate model organisms must be chosen.
the choice of model organism is fundamental to the type of science which can be pursued.
the range of organisms and genes in this collection, offer novel possibilities for future descriptive, comparative and hypothesis driven research.
in order to make the large amount of information described by this paper available for public scrutiny, it is available as  <dig> additional files which can be accessed online.
RESULTS
origin of this gene collection
the origin of this manual gene curation project lay in an attempt to annotate some of the plant ntmc <dig> genes i had identified previously  <cit> .
i chose the unique ntmc <dig> nomenclature to emphasise the phylogenetic relationships among these plant genes and to distinguish them from animal syt genes.
plant ntmc2t <dig>  ntmc2t <dig> and ntmc2t <dig> genes resemble animal syt genes in having the same domain architecture.
that this shared domain architecture implies functions in common, has now been demonstrated experimentally  <cit> .
the a. thaliana ntmc2t <dig>  protein  possesses calcium and membrane binding activities which allow it to function in a plasma membrane repair pathway induced by stress.
while using the annotation facilities at doe-jgi  <cit>  i realised that the evolutionary range of genome sequences available for inspection at doe-jgi, could be used to search for earlier origins of the embryophyte ntmc <dig> genes and metazoan syt genes.
the genome sequence of m. brevicollis  <cit>  has demonstrated that, as theorised  <cit>  the evolutionary transition from unicellularity to multicellularity involved abundant domain shuffling in proteins.
m. brevicollis, a unicellular protist, closely related but basal to metazoans, does not possess syt genes.
it does have a gene with a degree of gene structure similarity, domain architecture similarity and amino acid sequence similarity to plant ntmc <dig> and animal esyt proteins xm_ <dig>  proteins with similarity to ntmc <dig> and esyt proteins  are present in other unicellular eukaryotes  <cit> .
it seems reasonable therefore, to propose that while ntmc <dig> and esyt genes could ultimately be homologous, having attained their current forms through descent from a unicellular ancestor, syt genes, distinguished by their conserved gene structure, more likely represent a metazoan novelty.
on finding uniquely metazoan syt homologues, with their characteristic gene structures, in the genomes of t. adhaerens and n. vectensis  <cit>  i decided to compile afresh, the syt repertoire of the currently available genomes of evolutionary significance.
in addition to syt genes, i set out to gather esyt, rph3a, rph3al and doc <dig> homologues.
these gene families, while clearly distinct from the syt family at the gene structure level, encode proteins with high level sequence similarity to the tandem c <dig> domain region of syt  or with similarity, in common with ntmc <dig> proteins, at a more crude domain architecture level .
because of these similarities in form, there will be common functions.
a collection like this can be used to further define these functions, by revealing the patterns of sequence conservation  and sequence divergence, apparent in:  <dig>  the individual domains of an orthologous group of proteins,  <dig>  the equivalent domain in different groups, not necessarily homologous.
such knowledge can provide clues about how these proteins function as molecular machines in cells and about how these functions have developed over evolutionary time.
i used the annotation facilities at doe-jgi to create and annotate gene models where possible.
full details of all genes are listed in additional files  <dig>   <dig>   <dig> and  <dig>  figures  <dig>   <dig> and  <dig> summarise the genes in additional files  <dig>   <dig> and  <dig> respectively.
gene identification
i searched for homologues among the metazoan genome sequences at doe-jgi  <cit>  ncbi  <cit>  and ucsc  <cit> .
i employed a variety of methods, first to identify, then to piece together, gene sequences for homologues of syt, esyt, rph3a,rph3al and doc <dig>  these are generally large, multi-exon genes which encompass regions of high  and low  levels of sequence conservation.
as a first step, i used blatp  <cit>  or tblastn  <cit>  with amino acid sequence probes representing already known homologues, to identify genomic loci with potential to encode a similar protein.
i then examined a 3-frame translation  <cit>  of the genomic region to look at the possible exon-intron structure.
intron positions and phases are generally well conserved among these genes and serve as a useful diagnostic tool.
in cases where the draft genome sequence was highly fragmented, homologues were pieced together exon by exon.
where sequence conservation between amino acid probe and novel genomic locus was high, it was straightforward to identify exons and introns.
where sequence conservation was low, the conservation of intron position and phase could serve as a guide to judge possible exon boundaries.
in all cases, when transcript sequences from the locus under examination, or from a homologous locus in a closely related organism, were available, these were used to help identify exons.
occasionally, trace archive reads were sought to confirm exon sequences, where there was a gap or where there appeared to be frameshift errors in a section of draft sequence for example.
details of such instances are noted in additional files  <dig>   <dig> and  <dig>
where transcript sequences indicated the expression of alternative variants, i collected the variants.
in the case of m. musculus and h. sapiens, where there is abundant transcript variation, particularly at certain genomic loci, i collected only those variants indicated by two or more transcript sequences.
where evidence of alternative splicing exists but is not represented by full length transcripts  the resultant gene predictions require validation, as do all gene predictions for which no transcript evidence exists.
for genomic loci with conserved variant exons, which could be employed to extend the length of the encoded protein but where transcript evidence was lacking, the gene prediction includes all of the conserved exons.
transcript evidence for the gene predictions is listed in additional files  <dig>   <dig>   <dig> and  <dig>
once i had gathered a set of crude gene predictions, i compared the sequences in order to refine the predictions.
exon-intron junctions, gaps, novel regions of sequence conservation useful as probes to fill gaps, were all carefully examined using multiple alignments  <cit> .
in this way, it was possible to extend and improve the gene predictions significantly.
many of the gene predictions in this collection are still incomplete, at the n-terminal region of syt proteins in particular.
when there was no transcript coverage and no detectable sequence conservation or very short coding exons, it was not possible to make a prediction.
all incomplete predictions are noted in additional files  <dig>   <dig>   <dig> and  <dig>
sequence comparison: global comparison and nomenclature
having collected  <dig> genes from  <dig> genomes, i used clustalw <dig> at ebi  <cit>  to compare subsets of amino acid sequences, in order to examine their relationships of similarity overall.
i compared sets which excluded expressed variants but which included one sequence per gene, whether complete or not.
comparisons of invertebrate sequences  vertebrate sequences  invertebrate  plus mammalian sequences  and invertebrate  plus a subset of vertebrate sequences  were used to examine relationships and assign a suitable name to each gene.
i assigned a common name to those genes which clustered together in conserved groups.
when a conserved group included a homologue from m. musculus/h.
sapiens, the assigned name was chosen to be consistent with that in m. musculus/h.
sapiens.
for those syt genes which did not consistently cluster together and are not conserved among different organisms, i assigned them numbers upward from the last conserved group.
this process led to the identification of  <dig> conserved groups of syt genes plus a further  <dig> unique syt genes named with the syt stem symbol plus a number identifier, in order to provide a recognisable name  for each.
this naming system allows further sampling of the tree of life by gene sequencing to reveal:  <dig>  more unique syt genes, from syt <dig> upwards, all members of a recognisably homologous family and  <dig>  previously unrecognisable patterns of conservation among syt genes, from syt <dig> upwards.
i contacted the relevant organism nomenclature authorities  <cit>  to discuss this nomenclature.
this wide consultation led to agreement to change fam <dig>  a temporary hgnc symbol for this poorly characterised gene family, to esyt  <cit> .
where current approved nomenclature is discordant with that proposed here, the approved nomenclature is indicated within brackets, in additional files  <dig>   <dig>   <dig> and figures  <dig>   <dig> and  <dig>  there is inconsistency in the application of my naming rationale in drosophila and caenorhabditis, where it was more practical to stick with the approved nomenclature than propose new names.
as yet, and in comparison to vertebrates, invertebrates are poorly represented by organised gene nomenclature authorities.
my proposed nomenclature for a large number of genes in this collection, therefore, awaits sanction.
note that due to differing nomenclature conventions, the gene symbols in different organisms, range from having no capital letters to all capital letters.
the generic form, which i have used to label the genes from organisms outside the purview of nomenclature authorities, is first letter in capitals.
this designation will be used here, to describe the gene in general, without reference to a specific organism.
deciding whether relationships of orthology or some other kind exist among a group of similar genes, can be complicated.
the processes of duplication, deletion, and rearrangement of genes, plus the action of selection or absence of selection on gene sequences, can obscure relationships.
whether absence from a genome assembly reflects true absence, or incompleteness of the assembly, presents a further problem of some seriousness  <cit> .
in m. musculus and h. sapiens there are several groups of duplicated syt genes: group  <dig>  group  <dig>  group  <dig>  and group  <dig> .
each group contains genes which share a distinct gene structure.
additional files  <dig>   <dig>   <dig>   <dig> and figure  <dig>  indicate that the parental gene of group  <dig> is syt <dig>  the parental gene of group  <dig> is syt <dig>  that of group  <dig> is syt <dig> and that of group  <dig> is syt <dig>  in each of these cases, the parental gene is more similar to a gene present in an evolutionarily more ancient genome than it is to another member of the mammalian group, which is to say, orthologous  <cit> .
orthologous genes in different organisms, are related by vertical descent from a common ancestor.
relevant pairwise comparison scores are listed in figure  <dig>  it is worthwhile trying to determine the phylogenetic relationships between the syt genes of mammals and those of other animals because it is the rodent syt proteins which are best characterised functionally.
however, because mammalian organisms and their brains in particular, represent biological systems of such complexity, modelling these systems by using simpler systems is essential.
the relationships and nomenclature proposed here, are in good general agreement with the syt data at treefam  <cit>  which does not include many of the invertebrate syt genes in this collection.
this gene collection includes several previously unrecognised syt groups which are absent from m. musculus/h.
sapiens.
members of the syt <dig> group appear in vertebrates and invertebrates, but the syt <dig> group has yet to be found beyond vertebrates.
neither of these groups are well conserved, both genes having diverged significantly over time.
the syt <dig> group ranges from members encoding: tandem c <dig> domains  a single c2b domain  or no c <dig> domains at all .
syt <dig> appears in only two hexapod species in this collection.
the other syt groups absent in vertebrates, are sytalpha and syt <dig>  of these, the sytalpha group has a higher level of sequence conservation.
the sytalpha designation applies to all members of the conserved group which contains the d. melanogaster sytalpha gene .
there are many invertebrate syt genes, related to some degree, at the gene structure level, which are not conserved among different phyla .
n. vectensis, for example, has a very large number ofunique and divergent syt genes .
most of these encode stereotypical syt proteins, but some have only one c <dig> domain and some lack a tm domain .
the caenorhabditis syt repertoire, in common with many others in this collection, demonstrates that variant syt proteins which lack the tm domain, are usually included.
they are encoded either by separate genes or by alternative splicing.
caenorhabditis snt- <dig> is very similar to snt- <dig> at the amino acid sequence level, but does not encode a tm domain.
in the ciona syt repertoire, it is syt <dig> which produces a variant lacking the tm domain, but by alternative splicing.
in mammals, syt <dig> does not encode a tm domain, but a number of other mammalian syt genes which do encode a tm domain, express splice variants which lack the tm domain .
transcript variants are listed in additional files  <dig>   <dig>  and  <dig>
sequence comparison: gene duplications
the duplicate genes in the fish genomes presented a nomenclature challenge and i am very grateful for help from the zfin nomenclature coordinator  <cit>  in choosing the most appropriate names for these genes based on gene structure similarity and synteny.
a number of genes present as a single copy in mammals  are present as duplicates in the fish genomes.
in addition, syt <dig> and syt <dig> are duplicated in d. rerio .
it is notable that the single rph3a of tetrapods appears to have undergone duplication and fission, giving rise to rph3al and doc2b genes .
rph3al and doc2b are linked as a tandem pair in tetrapod genomes .
invertebrates have a single, equivalent rph, but no doc <dig> genes.
it is possible that d. pulex rphl represents an n-terminal rph fission product.
in the fish genomes, there are duplicate rph3a genes, but only in the lineage leading to d. rerio, has a supposed rph3a duplication led to the retention of linked rph3al and doc2b fission products .
in the other fish genomes, doc2b genes are present, along with the paralogous doc2a and doc2d genes, as well as duplicate rph3aa and rph3ab genes, so doc2b genes are not always linked to a rph3al fission product.
apart from the gene duplications already noted, other paralogous gene expansions are evident in this collection .
in the case of esyt gene duplications, retroposition of processed transcripts has occasionally led to gene duplication .
sequence comparison: multiple alignments
in order to illustrate the conserved features of the proteins in this collection, i used multalin  <cit>  to construct multiple alignments of subsets of sequences, including expressed variants .
the conserved groups of vertebrates and invertebrates are compared separately but the syt <dig> group  contains members from vertebrates and invertebrates.
all of the syt sequences from n. vectensis  are compared.
these multiple alignments illustrate the patterns of sequence conservation  distinctive of different groups.
common patterns of alternative splicing which serve to alter particular domains, are also illustrated.
where calcium coordinating residues are absolutely conserved  these are indicated by arrows.
these absolutely conserved sequences are highlighted in red in figures  <dig>   <dig> and  <dig>  in the many more cases where calcium coordinating positions are not absolutely conserved, there is a range, down to zero conservation, of the ten sites present in the absolutely conserved, calcium responsive syt proteins.
a conserved region 20- <dig> residues upstream of the c2a domain, often recognisable as a variant of this motif: grikpely is indicated where present.
conservation of this motif was noticed previously  <cit>  but it is obvious here, that this region is much more widely conserved.
its function is yet to be investigated.
the position and phase of introns are also illustrated, in order to show the gene structures distinctive of different family groups, as well as instances of intron mobilization.
intron mobilization within a genus is specifically noted in additional file  <dig>  intron phase is indicated by colour: black for phase  <dig>  red for phase + <dig> and blue for phase + <dig>
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
the ten conserved calcium coordinating positions are indicated by arrows.
the c2ab region is highly conserved whereas the n-terminal portion has much more sequence variation.
expressed variants are included to demonstrate the regions altered by alternative splicing and rna editing.
the region between the tm domain and the c2a domain is commonly altered, for instance in s. purpuratus, b. floridae, a. gambiae and d. melanogaster.
the next region expressed as alternative forms, is the region between c2a and c2b, as indicated by d. pulex and a. gambiae variants.
the third region to be altered, by alternative splicing or by rna editing, is the central portion of the c2b domain.
in caenorhabditis, this recoding alters not only the identity of positions within this exon of  <dig> amino acids length in all other cases, but extends its length by  <dig> or  <dig> amino acids.
additional file  <dig> also includes the snt- <dig> sequences from caenorhabditis.
the snt- <dig> gene encodes a protein very similar to the snt- <dig> protein which is the syt <dig> orthologue.
the snt- <dig> proteins however, lack the n-terminal membrane anchoring portion.
the snt- <dig> proteins also differ in the region between c2a and c2b.
additional file  <dig> shows a multiple alignment of the vertebrate syt <dig> protein sequences in this collection.
the ten conserved calcium coordinating positions are indicated by arrows.
the c2ab region is highly conserved whereas the n-terminal portion has a little more sequence variation.
expressed variants are included to demonstrate the regions altered by alternative splicing.
this time, the alteration is restricted to the region between the tm and c2a domain.
conserved n-glycosylation consensus sites are indicated by blue boxes.
in the first coding exon, the n-glycosylation site and the conserved upstream o-glycosylated threonine, are known to be important functional sites in vivo  <cit> .
additional files  <dig>   <dig> and  <dig> show multiple alignments of the vertebrate syt <dig>  syt <dig> and syt <dig> protein sequences, respectively.
for syt <dig>  syt <dig> and syt <dig> genes, there is transcript evidence of alternatively expressed forms which lack the tm domain.
syt <dig>  syt <dig>  syt <dig> and syt <dig> are all related by duplication, with syt <dig> being the parental gene .
in the fish genomes, there are yet more duplicates.
in order to distinguish the relationships among these duplicates, pairwise comparison scores  and syntenic relationships  were examined.
the top scores highlighted in blue in figure  <dig>  indicate that d. rerio syt1a is orthologous to the syt <dig> gene of other vertebrates and that d. rerio syt <dig> is orthologous to the syt <dig> gene of other vertebrates.
the pairwise comparison scores do not however, indicate orthology for the other duplicates.
figure  <dig> shows that relationships of orthology, are indeed present.
the neighbouring genes of syt <dig>  syt <dig>  syt <dig> and syt <dig>  in the genomes of h. sapiens, x. tropicalis and d. rerio, reveal conserved syntenic blocks in each genome.
this synteny indicates that the four genomic segments share ancestry.
orthologous relationships exist between: h. sapiens syt <dig> and the duplicated syt1a/syt1b in d. rerio, between syt <dig> and syt2/syt2b and between syt <dig> and syt5a/syt5b, respectively.
since there is so much interest in the mechanism of syt <dig> function, close examination of the syt <dig> orthologues and paralogues in this collection, to see how they have been shaped by evolution, may prove productive.
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
the variants of d. pulex differ at the c-terminal end and those of a. gambiae differ at the n-terminal end.
additional file  <dig> shows a multiple alignment of the vertebrate syt <dig> protein sequences.
the variants of g. gallus differ in the region between c2a and c2b.
additional file  <dig> shows a multiple alignment of the protein sequences of syt <dig>  the vertebrate paralogue of syt <dig>  here, the variants in g. gallus, m. musculus and h. sapiens differ in the c2b domain.
recent functional studies of mouse syt <dig>  <cit>  have built on initial gene knockout results  <cit>  which suggested that syt <dig> dysfunction could be responsible for some human psychiatric disease.
syt <dig> dysfunction has also been implicated in human psychiatric disease  <cit> .
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
here, there is evidence from d. melanogaster, of variant forms lacking the tm domain.
additional file  <dig> shows the vertebrate syt <dig> sequences.
in comparison to the invertebrate syt <dig> genes, the vertebrate syt <dig> genes contain additional coding capacity in the region between the tm and c2a domains.
several additional exons which can be alternatively spliced, are available to extend the length of the region linking the tm and c2a domains.
a study of constitutive versus alternative splicing, has examined the syt <dig> alternative exons  <cit> .
recent gene knockout studies in mice  <cit>  appear to confirm a role for syt <dig> in glucose homeostasis.
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
the conserved motif upstream of the c2a domain is indicated.
additional file  <dig> shows the vertebrate syt <dig> sequences.
the variants expressed by m. musculus and h. sapiens differ at one or other end.
additional files  <dig>   <dig> and  <dig> show the protein sequences of the vertebrate syt <dig> paralogues, syt <dig>  syt <dig> and syt <dig>  respectively.
all of these proteins retain all ten calcium coordinating residues and are thus likely to have been selected to act as calcium sensors.
they also retain the motif of unknown function upstream of the c2a domain.
transcript evidence from x. tropicalis, m. musculus and h. sapiens indicates the expression of alternative forms of syt <dig> which lack the tm domain.
the c-terminal region of syt <dig> is also expressed as alternative forms.
functional studies have examined the conserved cysteine residues of the n-terminal region which are unique to this paralogous family in vertebrates  <cit> .
syt <dig> is thought to play a role in sperm acrosomal exocytosis  <cit> .
mouse syt <dig> has been identified as upregulated during meiosis in fetal ovaries  <cit> .
note that some investigators working with syt proteins, use different synonyms for the protein products of certain syt genes.
the officially approved gene nomenclature is used throughout this paper.
in the functional literature, the protein products of syt <dig> and syt <dig> have most often been referred to differently  <cit> .
syt <dig> and syt <dig>  so far poorly characterised, have also been prone to a similar lack of recognition  <cit> .
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
while there is some variation in length, of sequences from different organisms in the c2b region, the more pronounced difference is a trend towards increased length of the region between tm and c2a domains in insects.
additional file  <dig> shows the vertebrate syt <dig> sequences.
transcript evidence indicates the expression of forms varying in the n-terminal region in h. sapiens.
little functional work has been done so far with syt <dig>  <cit> .
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
in most of these genes, there are two alternative n-terminal coding exons, only one of which encodes a tm domain.
these sequences also retain the motif of unknown function just upstream of the c2a domain.
additional file  <dig> shows the vertebrate syt <dig> sequences.
transcript sequences from h. sapiens indicate the expression of a form lacking the tm domain.
a recent study suggests that human syt <dig> may function as a tumour supressor  <cit> .
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
in these sequences, it is the c2a domain which is more highly conserved than the c2b domain.
the motif upstream of the c2a domain is present.
there is a large amount of variation in amino acid sequence, intron position and intron phase in the n-terminal portion of these sequences.
additional file  <dig> shows the vertebrate syt <dig> sequences.
transcript evidence from m. musculus and h. sapiens indicates the expression of variants differing at the c-terminal end.
transcripts from h. sapiens indicate variant forms at the n-terminal end, which lack the tm domain.
very little is known about the function of syt <dig>  <cit> .
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
transcript evidence indicates the expression of variants in c. savignyi, c. intestinalis and s. purpuratus which differ at the n-terminal end and lack the tm domain.
additional file  <dig> shows the vertebrate syt <dig> sequences which contain a tm domain.
additional file  <dig> shows the vertebrate syt <dig> sequences which lack a tm domain.
additional file  <dig> shows vertebrate syt <dig> sequences.
transcript evidence in x. tropicalis, m. musculus and h. sapiens indicates the expression of variant forms of syt <dig> with altered c2b domains.
human syt <dig> transcripts also differ at the n-terminal region, resulting in forms which lack the tm domain.
loss of function of human syt <dig> has been associated with neurodevelopmental abnormalities  <cit> .
in the human genome, there is a repeat of a portion of the syt <dig> gene  on chromosome  <dig>  apparently the result of retroposition of a syt <dig> transcript.
the chromosome  <dig> locus has been named syt14l by hgnc  <cit> .
syt14l is classed as a pseudogene by ncbi  and is not included in this collection.
syt <dig> homologues
additional file  <dig> shows a multiple alignment of the invertebrate syt <dig> protein sequences in this collection.
a possible motif upstream of the c2a domain is indicated.
additional file  <dig> shows the vertebrate syt <dig> sequences.
transcript evidence in m. musculus and h. sapiens indicates the expression of variants altered at n-terminal and c-terminal ends.
a possible motif upstream of the c2a domain is indicated.
the protein product of syt <dig> has been referred to as b/k protein, but not much is known about its function  <cit> .
sytalpha homologues
additional file  <dig> shows the sytalpha sequences in this collection.
a possible motif upstream of the c2a domain is indicated.
the only functional study relevant to sytalpha is its anatomical localisation in d. melanogaster  <cit> .
syt <dig> homologues
additional file  <dig> shows the syt <dig> sequences in this collection.
these proteins are so far completely uncharacterised.
syt <dig> homologues
additional file  <dig> shows the vertebrate syt <dig> sequences in this collection.
these proteins are so far completely uncharacterised.
syt <dig> homologues
additional file  <dig> shows the invertebrate syt <dig> sequences in this collection.
these proteins are so far completely uncharacterised.
syt homologues in n. vectensis
additional file  <dig> shows all of the syt proteins of n. vectensis, demonstrating enormous sequence variation.
these proteins are so far completely uncharacterised.
dblc <dig> homologues
additional file  <dig> shows the dblc <dig> sequences in this collection.
transcript variants are indicated.
these proteins have not previously been reported and are so far completely uncharacterised.
dblc <dig> sequences are detectably similar to syt sequences at the amino acid sequence level and occur in the genomes of marine invertebrates in this collection.
while they lack an n-terminal tm domain, they have a tandem c <dig> domain architecture in common with syt proteins.
dblc <dig> genes lack similarity to syt genes at the gene structure level, so a new gene symbol was chosen to represent these genes.
the dblc <dig> designation represents a fuller description: double c <dig> domain protein.
dblc <dig> is thus similar but recognisably different from doc <dig>  the designation for the homologous family of double c <dig> domain proteins of vertebrates.
esyt homologues
additional file  <dig> shows the invertebrate esyt <dig> sequences in this collection.
transcript variants are indicated.
the arthropod esyt <dig> genes include mutually exclusive alternative exons which serve to alter the specification of the first c <dig> domain.
additional file  <dig> shows the sequence of the first of these alternatives only.
additional file  <dig> shows all versions of this alternatively coded region.
additional file  <dig> shows the vertebrate esyt <dig> sequences.
transcript variants are indicated.
in the fish genomes, esyt <dig> has duplicated, resulting in esyt1a and esyt1b.
the esyt1b duplicate has acquired an internal duplication which contributes an additional four c <dig> domains to the product of this gene  <cit> .
within this duplicated section, the pink dot marks the position of an apparent phase  <dig> intron loss in the tetraodontiform fish.
in each of the three other occurrences  of this portion of the esyt1b gene, the phase  <dig> intron is present.
this intron is also present in each of the two repeats of this sequence, in the esyt1a genes.
additional file  <dig> shows the vertebrate esyt <dig> sequences.
transcript variants are indicated.
in the fish genomes, esyt <dig> has duplicated, resulting in esyt2a and esyt2b.
additional file  <dig> shows the vertebrate esyt <dig> sequences.
transcript variants are indicated.
similarity between syt genes and what are now named esyt genes, was first evident from genome sequence comparisons  <cit> .
in trying to classify and annotate the genes within the human genome, hgnc  <cit>  noticed that these genes formed a distinct gene family, separate from syt genes.
as there was little functional data associated with these genes, a temporary nomenclature with no functional implication was provided: fam62a, fam62b and fam62c.
the functional data which did exist, was in the form of a cdna representing rat esyt <dig>  which had been cloned during an investigation of adipocyte proteins and found to encode a membrane bound c <dig> domain protein with similarity to proteins in plants  <cit> .
further genome sequence comparisons indicated that syt, fam <dig> and plant ntmc <dig> genes encode proteins with amino acid sequence similarity and domain architecture similarity, but belong to distinct gene families  <cit> .
apart from  <cit>  and until recently, functional data associated with plant ntmc <dig> and animal fam <dig> gene families had been lacking.
an initial investigation of proteins encoded by the human fam62a, fam62b and fam62c genes, has now been published and the authors named these proteins esyt <dig>  esyt <dig> and esyt <dig>  to stand for extended synaptotagmin-like proteins  <cit> .
consultation with the nomenclature committees during the preparation of this paper, led to the decision to implement the esyt nomenclature in place of fam <dig>  in an initial pair of studies on one member of the plant ntmc <dig> gene family  <cit>  the authors put forward the name syt <dig> for this particular gene and gene product.
as discussed in the background section of this paper, when different communities make naming decisions based on the criteria most relevant to them  it is inevitable that multiple names will be used to identify the same entities.
rph3a homologues
additional file  <dig> shows the invertebrate rabphilin sequences in this collection.
in caenorhabditis, alternative n-terminal exons can be used.
sequence conservation among rabphilins is high in the n-terminal, rab binding portion  <cit>  and the c-terminal tandem c <dig> domain portion, but a large middle portion is very poorly conserved.
this means that where transcript confirmation is absent, the gene predictions across this portion are unlikely to be accurate.
additional file  <dig> shows the vertebrate rph3a sequences.
transcript variants are indicated.
in vertebrates, all  <dig> calcium coordinating positions are absolutely conserved, but not all are conserved in invertebrates.
in additional file  <dig>  the positions of these  <dig> amino acids are indicated by pink dots.
the position of the fifth calcium coordinating residue in the c2a domain is not strictly conserved, but in most cases a suitable aspartate or glutamate residue is present one residue earlier.
i am not aware of biochemical evidence for calcium binding by invertebrate rabphilin proteins, but genetic evidence  <cit>  suggests that c. elegans rbf- <dig> at least, appears to function similarly to its mammalian homologues.
for this reason, in figures  <dig>  and  <dig>  i have counted the fifth position as present when a suitable aspartate or glutamate residue occurs at the usual spacing or one residue earlier.
rph3al homologues
additional file  <dig> shows vertebrate rph3al sequences.
transcript variants are indicated.
in the functional literature, products of rph3al genes have been referred to as noc <dig>  <cit> .
as outlined above, rph3al genes appear to be the result of duplication and fission of an ancestral rph3a gene, producing linked rph3al and doc2b genes in tetrapods and in d. rerio.
rph3al represents the n-terminal portion of the ancestral rph3a gene and doc2b represents the c-terminal portion, as seems clear from a comparison of the gene structures of rph3a, rph3al and the linked doc2b genes .
as also indicated above, doc2b genes do not always accompany rph3al genes, as they exist in other fish genomes where rph3al genes do not occur .
doc <dig> homologues
additional file  <dig> shows the vertebrate doc <dig> protein sequences.
transcript variants are indicated as well as the positions of the calcium coordinating residues, which are conserved in all of the doc <dig> proteins except doc2g.
i have included potential doc2g products from the human gene locus.
although spliced transcripts are produced from this locus, they cannot encode full length proteins and it is not known whether protein products are translated from these transcripts.
it is clear from the shared gene structure of doc <dig> genes and rph3a genes, that doc <dig> genes represent the c-terminal portion of a rph3a gene.
since rph genes exist in invertebrates, but doc <dig> genes do not, it seems reasonable to suggest that doc <dig> genes arose via duplication and fission of a rph3a gene.
evidence for this, in the form of the two linked fission products  is present in the genomes of d. rerio, x. tropicalis, a. carolinensis, g. gallus, t. guttata, o. anatinus, m. domestica, m. musculus and h. sapiens.
in this collection, the fish genomes reveal a novel lineage of doc <dig> genes, which in consultation with zfin staff  <cit>  were given the gene symbol doc2d, to stand for doc <dig> delta.
the first doc <dig> gene, cloned in  <dig>  <cit>  was soon found to represent the first member of a family of three genes in rodents.
the protein products of these three genes were named doc <dig> alpha, doc <dig> beta and doc <dig> gamma  <cit> .
the conserved region present in the n-terminal portion of doc <dig> sequences, specifies the munc <dig> binding capacity of doc <dig> proteins  <cit> .
CONCLUSIONS
this paper describes a manually curated collection of genes pertinent to intercellular trafficking in multicellular animals.
the collection includes  <dig> genes from the genome sequences of  <dig> marine invertebrates   <dig> genes from the genome sequences of  <dig> ecdysozoans, including the  <dig> drosophila and  <dig> caenorhabditis genomes  plus  <dig> genes from  <dig> vertebrate genomes .
when viewed overall  the collection implies that particular homologues which are always present in the genome repertoire, serve to specify functions required for animal life.
the basis for this view is that these genes have not been lost, but have been retained and conserved, over the long period of metazoan evolution.
such a view about the essentiality of gene function, differs from a synthetic biology view, which considers the genes necessary to specify an operable system which functions outside the natural world.
it differs too, from a reductionist molecular biology view where, if a gene can be deleted in the laboratory setting, without killing the organism, it is not considered essential.
such gene deletion experiments have been done with syt <dig> in different organisms  <cit>  and according this view, syt <dig> is not strictly essential to animal life because snt- <dig> null mutants are not lethal in caenorhabditis.
in drosophila and mus however, syt <dig> null mutants are indeed lethal as they cannot survive to reproduce as adults.
while my previous genomic comparison of syt genes  <cit>  could have been interpreted as indicating a correlation between increased numbers of syt genes and increased organism complexity  <cit>  the analysis here, which is much more comprehensive, definitively rules that conclusion out.
this collection of genes offers prospects for fruitful speculation about the functional roles of the syt repertoires of different animals and includes a great range of biological complexity.
the conclusions drawn from evolutionary genomics approaches are always provisional, because they depend on the quality of information available  and are thus subject to review when better information becomes available.
conclusions drawn from reassembled and reannotated genomes will inevitably be more accurate than those drawn from the first draft.
the d. rerio genome being repeat rich, is taking time to assemble into a complete version.
among the genomes in this collection, the draft d. rerio genome ranks first in terms of gene number.
d. rerio surely presents a useful model system with which to investigate the functional divergence of close paralogues  <cit> .
currently, it is the rodent syt proteins and, with the advent of transgenic mouse studies, their respective genes, which garner most interest and have been best characterised.
work with syt genes in other animal species has provided information supporting hypotheses about how syt genes fulfil their functional roles  <cit> .
to date, these studies indicate a primary role for syt <dig> in specifying a neuron specific synaptic vesicle protein which serves to trigger exocytosis.
the genomic specification of syt <dig> expression patterns, has been investigated and utilised in ascidian species  <cit> .
in c. elegans, the anatomical expression patterns of the mutually exclusive alternative exon variants of syt <dig> have been mapped  <cit> .
this gene collection shows that such mutually exclusive exon deployment has arisen on a number of occasions in invertebrates, altering c <dig> domains in syt <dig>  syt <dig> and esyt <dig> proteins .
the repeated evolution of this type of c <dig> domain alteration, suggests that careful comparison of the altered forms might prove a useful focus for investigations aimed at understanding how these c <dig> domains operate as molecular machines.
the feature most widely shared among the great diversity of syt proteins in this collection, is sequence conservation of the c-terminal c2b region, which suggests that a defining function of syt proteins is located here.
the regulatory machinery controlling variant expression can be illuminated by comparative analyses of genomic sequences, as was done to investigate the rna editing of syt <dig> which occurs in hexapods  <cit> .
another comparative analysis  <cit>  examining the genomic specification of alternative splicing of syt <dig>  concluded that splice variants expressed in d. melanogaster, c. pipiens, b. mori, t. castaneum and p. humanus depend on correctly coded intronic, cis-acting regulatory sequences acting on the requisite splice donor/acceptor sequences.
in this collection, there is evidence from a. gambiae, d. melanogaster, d. pulex, s. purpuratus, b. floridae, d. rerio, x. tropicalis, a. carolinensis, g. gallus, m. musculus and h. sapiens of rna editing and alternative splicing of syt <dig> .
these sequences may aid future efforts to identify the cis-sequences involved in transcript processing.
the biochemical understanding of transcript production and processing, acquired at great pains over many years  is not quite up to the task of extracting the full meaning from the flood of new transcript data from high throughput sequencing.
an enormous amount of experimental work will need to be done to understand the varied roles of transcripts in complex organisms.
the problem of identifying and annotating different kinds of transcript sequences is discussed in  <cit> .
work to assess the functional significance of alternatively expressed syt genes has barely begun  <cit>  so it is difficult at present, to judge the importance of the transcript variants listed in this collection.
evolutionary conservation of variant production, probably argues for functional significance.
new developments in mass spectrometry  <cit>  are bound to aid attempts to examine complex mixtures of proteins such as synaptic vesicle preparations  <cit> .
it is more likely however, that the functional impact of syt protein variants will eventually be revealed by research focussed on syt proteins of known importance.
the significance of the conservation of intron position and phase within syt genes could be investigated by gene targeting in transgenic mice, now a common route to investigate syt roles  <cit> .
the best conserved intron, a signature feature of syt genes, is the phase + <dig> intron at the start of the third beta strand of the c2b domain  <cit>  .
a single point mutation to inactivate this intron acceptor dinucleotide would be predicted to have devastating consequences for the expression of the host gene.
in contrast to the approach of disabling a gene by creating a large deletion, the production of a null mutant by such a small intervention would be unlikely to cause collateral damage.
the chances of cryptic splice acceptors being deployed to restore normal function to such a targeted syt splice site point mutant, seem remote, based on the importance of this region to syt protein function.
the presence of this intron in so many syt genes could be interpreted as evidence that evolutionary alteration of this site is hazardous.
an example of a gene inactivating mutation which is not rescued by altered splicing, is present in this gene collection.
the h. sapiens gene which is orthologous to m. musculus doc2g has acquired a frameshift mutation in the fifth coding exon which introduces a stop codon.
human transcript sequences indicate that an alternative intron acceptor dinucleotide capable of restoring the reading frame, which exists  <dig> nucleotides upstream of the exon  <dig> acceptor, is not used.
the regulatory sequences necessary to select this alternative acceptor are not present within this short intron .
it is not known what protein products, if any, are produced by this gene in h. sapiens and it appears likely that the human gene is non functional.
this observation underscores the importance of splicing signals which are poorly understood at present, as well as the importance of transcript analysis to verify gene predictions.
in sea urchin species, the syt <dig> protein has been shown to act in cortical granule exocytosis  <cit> .
syt <dig> products are used as pan neuronal markers in echinoderms, hemichordates and urochordates  <cit>  and the role of syt <dig> in synaptic vesicle exocytosis is assumed.
few comparative studies have attempted to map the expression patterns of the syt repertoire of a particular organism but all of the syt genes of d. melanogaster  <cit>  and a large fraction of those of r. norvegicus  <cit>  have been examined in this way.
t. adhaerens is thought to lack cells equivalent to neurons, but it has a well conserved syt <dig> as well as genes encoding the requisite neuronal snares and complexin for fast neurotransmission .
although never observed in the laboratory, t. adhaerens is expected to produce gametes in the wild  <cit>  so a role for the t. adhaerens syt <dig> similar to that of sea urchin syt <dig> in cortical granule exocytosis is possible, alongside any yet to be discovered 'neural' role.
n. vectensis, which has a bona fide nervous system, also has syt <dig> and the requisite neuronal snares and complexin .
n. vectensis is quite unusual among the metazoans collected here, in having so many highly divergent syt genes .
of all the n. vectensis syt genes, only syt <dig> and syt <dig> encode all ten calcium coordinating positions, indicating their likely role in calcium triggered membrane fusion.
future studies to determine the anatomical expression patterns of the n. vectensis syt genes in vivo, will be of great interest to further understand the evolutionary development of animal nervous systems.
initial study of the molecular anatomy of the n. vectensis nervous system has examined the expression patterns of neurotransmitters, neuropeptides, transcription factor and cell adhesion genes  <cit> .
synaptotagmin gene expression was used to label differentiated neurons in platynereis, in a study of the evolutionary development of the central nervous system of animals  <cit> .
capitella, like platynereis, is an annelid member of the lophotrochozoa, but has the advantage as a model organism, of a whole genome sequence  <cit> .
study of neurogenesis in capitella is now underway  <cit> .
in situ hybridization analysis of syt <dig> expression in capitella, shows a broad neural pattern .
it will be important in future, to describe and compare anatomical expression patterns and functions of syt genes, both conserved and unique, in vertebrates and invertebrates.
such comparisons will be essential to improve, or reduce confidence in functional annotation of genomes based on phylogenetic homology, or similarity otherwise defined.
good annotation of genomes is a necessary step in the process of biocuration  <cit>  a new field aimed at using twenty first century sequencing and information technology to make accessible to bench scientists, the vast amount of information they require, distilled, in order to take forward wisely informed hypothesis driven research.
community deposition of a variety of types of biological data into appropriate databases, database efforts to curate, update and integrate this data with relevant data in other databases, will if pursued successfully, lead to great progress in biology in the years ahead.
for the time being, the manual gene collection made available here, can be used immediately, by bench scientists engaged in hypothesis driven investigation of membrane trafficking by these proteins, to widen the bounds of a field already stretching from biophysical study in vitro at the atomic level, to phenotypic assessment of gene function in model organisms.
authors' contributions
the author had sole responsibility for the design and execution of this work.
supplementary material
additional file 1
comprehensive information about the marine invertebrate genes in this collection.
this spreadsheet file  contains full details of each gene identified in the marine invertebrate genomes in this collection.
click here for file
additional file 2
comprehensive information about the ecdysozoan genes in this collection.
this spreadsheet file  contains full details of each gene identified in the ecdysozoan genomes in this collection.
click here for file
additional file 3
comprehensive information about the vertebrate genes in this collection.
this spreadsheet file  contains full details of each gene identified in the vertebrate genomes in this collection.
click here for file
additional file 4
comprehensive information about the snare and complexin genes identified in t. adhaerens and n. vectensis.
this spreadsheet file  contains full details of the genes identified in t. adhaerens and n. vectensis which are similar to neuronal snares and complexin.
click here for file
additional file 5
dendrogram of relationships among the invertebrate sequences in this collection.
additional file  <dig> displays the guide tree of the clustalw <dig> comparison of the invertebrate sequences in this collection, excluding variants, totalling  <dig> sequences.
genes which encode mutually exclusive alternative exons are highlighted with a green box.
conserved groups of syt genes which have not previously been described, are highlighted with a blue box.
click here for file
additional file 6
dendrogram of relationships among the vertebrate sequences in this collection.
additional file  <dig> displays the guide tree of the clustalw <dig> comparison of the vertebrate sequences in this collection, excluding variants, totalling  <dig> sequences.
mouse genes are highlighted with a red box.
conserved groups of syt genes which have not previously been described, are highlighted with a blue box.
click here for file
additional file 7
dendrogram of relationships among the invertebrate and mammalian sequences in this collection.
additional file  <dig> displays the guide tree of the clustalw <dig> comparison of the invertebrate sequences  excluding variants, plus mammalian sequences, excluding variants, totalling  <dig> sequences.
mouse genes are highlighted with a red box.
conserved groups of syt genes which have not previously been described, are highlighted with a blue box.
click here for file
additional file 8
dendrogram of relationships among the invertebrate and vertebrate sequences in this collection.
additional file  <dig> displays the guide tree of the clustalw <dig> comparison of the invertebrate sequences  excluding variants, plus a subset of vertebrate sequences, excluding variants, totalling  <dig> sequences.
mouse genes are highlighted with a red box.
conserved groups of syt genes which have not previously been described, are highlighted with a blue box.
click here for file
additional file 9
alignment of the invertebrate syt <dig> sequences, plus caenorhabditis snt- <dig>  amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 10
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
conserved n-glycosylation consensus sites are indicated by blue boxes.
the conserved threonine, which can be o-glycosylated is also indicated by a blue box.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 11
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
the x. tropicalis splice variant is included and highlighted with a black dot where it differs.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 12
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 13
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 14
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 15
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
the g. gallus splice variant is included and highlighted with a black dot where it differs.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 16
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 17
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
the d. melanogaster splice variant is included and highlighted with a black dot where it differs.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 18
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 19
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 20
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 21
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 22
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 23
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
click here for file
additional file 24
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
because of their differing positions, tm domains are highlighted in blue.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 25
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
the h. sapiens splice variant is included and highlighted with a black dot where it differs.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 26
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants which specify an alternative n-terminus lacking a tm domain, are indicated on top, with a dotted line indicating where they join the common c-terminal portion.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 27
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
the h. sapiens splice variant is included and highlighted with a black dot where it differs.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 28
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
because of their differing positions, tm domains are highlighted in blue.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 29
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 30
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants which specify alternative n-termini lacking tm domains, are included and highlighted with black dots where their sequences differ.
tm domains are highlighted in blue.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 31
alignment of the vertebrate syt <dig> sequences which have tm domains.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 32
alignment of the vertebrate syt <dig> sequences which lack tm domains.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 33
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 34
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 35
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
a possible motif just upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 36
alignment of the invertebrate sytalpha sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 37
alignment of the syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 38
alignment of the vertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
the x residue indicates where a portion of sequence is missing.
click here for file
additional file 39
alignment of the invertebrate syt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
tm sequences are highlighted in blue.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 40
alignment of the n. vectensis syt sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
tm sequences are highlighted in blue.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 41
alignment of the invertebrate dblc <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
click here for file
additional file 42
alignment of the invertebrate esyt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 43
alignment of the mutually exclusive alternative esyt <dig> exons.
in addition to the alternatively coded exon which is highlighted, the exon upstream and the exon downstream are shown.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
blue bars indicate phase + <dig> introns.
click here for file
additional file 44
alignment of the vertebrate esyt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
the middle portion is only present in the fish esyt1b sequences.
a pink dot in this portion, marks an intron loss in t. rubripes and t. nigroviridis.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 45
alignment of the vertebrate esyt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 46
alignment of the vertebrate esyt <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
some intron positions are marked with dotted lines.
in these cases, transcript sequence covers a gap in the genomic sequence and the intron presence is assumed.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 47
alignment of the invertebrate rabphilin sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are indicated.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by pink dots at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 48
alignment of the vertebrate rabphilin  sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 49
alignment of the vertebrate rph3al sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
x residues indicate where a portion of sequence is missing.
click here for file
additional file 50
alignment of the vertebrate doc <dig> sequences.
amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment.
splice variants are included and highlighted with black dots where they differ.
intron position and phase is indicated with a coloured bar between amino acids.
black bars indicate phase  <dig> introns.
red bars indicate phase + <dig> introns.
blue bars indicate phase + <dig> introns.
some intron positions are marked with dotted lines.
in these cases, transcript sequence covers a gap in the genomic sequence and the intron presence is assumed.
the five conserved acidic amino acids in each c <dig> domain are indicated by black arrows at the top of the alignment.
x residues indicate where a portion of sequence is missing.
click here for file
acknowledgements
this work was supported by the uk medical research council.
i wish to thank hugh pelham for making it possible for me to complete this work.

 synaptotagmin proteins participate in regulated membrane fusion events in multicellular organisms. syt research dates back to when the first syt was identified as an integral protein of synaptic vesicles. since then, much effort has gone into discovering the function and the detailed mechanism of action of this protein, syt numerous experimental approaches have shown that syt is crucially involved in fast neurotransmitter release at synapses. syt binds calcium, serving as the calcium sensor which triggers synaptic vesicle exocytosis. syt also binds to the neuronal snare proteins which are required for membrane fusion, as well as to membranes directly. the primary structure of syt revealed three important features: an n-terminal transmembrane domain, which serves to anchor the protein in the synaptic vesicle, plus two, tandem, c-terminal, cytoplasmic c domains which specify the calcium, snare and membrane binding properties. crystallographic studies of some syt c domains have revealed the exact nature of their calcium binding abilities. structural bioinformatics has shown that c domains are very common in the eukaryotic protein repertoire, but not all act to bind calcium. c domains can occur singly or as multiple copies in a given protein. the particular domain organisation of syt is important for its role in synaptic vesicle exocytosis. the tm domain tethers the protein to the membrane. the c2a and c2b domains, endowed with individual capacities to bind calcium, snare proteins, other accessory proteins and phospholipids, act both independently and synergistically, at different stages during the life cycle of the synaptic vesicle, to promote or inhibit fusion. after the genetic code for syt was discovered genes similar to syt were sought and found. biochemical methods which were employed initially, led to an increase in the size of the rodent syt family, from one member to thirteen members and beyond. syt relatives were also identified and studied in d. melanogaster and c. elegans. all of these homologous syt genes were found to encode proteins with a common domain architecture: an n-terminal tm domain connected by a variable length, poorly conserved linker sequence, to well conserved, tandem, c-terminal c2a and c2b domains and this domain architecture was taken to define the family. functional studies with the proteins in vivo and in vitro, indicate that some respond to calcium but others do not. whether or not they respond to calcium, syt proteins are able to regulate membrane fusion due to common properties of their c domains. studies of the anatomical expression patterns of syt genes show that syt is abundantly expressed in nervous tissue. in d. melanogaster, syt expression is neuron specific but other syt genes are expressed elsewhere. in rodents, most syt genes are expressed in the brain. with the publication of whole genome sequences from multicellular organisms, it became possible to identify and compare complete genome complements of syt genes using computational methods. the phylogeny of these genes could be examined by comparing sequences from different organisms. seventeen syt homologues were found in each of the complete h. sapiens and m. musculus genome sequences. these were named syt to syt according to the nomenclature conventions of the hugo gene nomenclature committee. some of these genes however, encode syt proteins which lack the requisite domain architecture for inclusion in a syt family defined by domain structure. the complete genome sequences of c. elegans, d. melanogaster and a. thaliana were analysed, together with draft genome sequences from other eukaryotes, to assess their syt coding potential. proteins which share the stereotypical syt domain architecture, are encoded by distinct gene families in plants and in animals further confounding the notion of a homologous syt family defined in terms of domain architecture. functional similarity between plant and animal syt proteins, at the level of calcium and membrane binding, has been demonstrated in keeping with the notion that protein form specifies function. proteins have been classified as syt family members by the presence or absence of suitable characteristics in terms of protein form and function without studious regard to gene ancestry. legitimately but confusingly, family membership can be defined in different ways: in terms of descent from a common ancestor or in terms of the common possession of distinguishing features. since it is not straightforward to demonstrate that homology is the cause of the similarity between biological entities such as proteins, it is often just assumed, leading to the confusion of two different notions of a protein family. please see for explanations of the terminology of homology. examination of the genomic specification of eukaryotic proteins can help to reveal their phylogenetic relationships because the evolutionary conservation of gene structure provides additional information on which to base an inference of homology. families may thereby be resolved as those related by homology and those otherwise related. ambiguity still remains however, because depending on the exact biological context, members of a family, however defined, can sometimes be considered as functionally interchangeable and thus worthy of a common functional identity. thus, genes and proteins can acquire multiple identities depending on the perspective from which they are considered. universal agreement about nomenclature for genes and proteins which is suitably clear and meaningful is a challenge for the future, but norms for gene nomenclature according to phylogenetic relationships are established. it may be the case that members of a gene family related by heredity, share a common functional role. however, due to the dynamics of eukaryotic genome evolution the functional attributes of gene relatives can diverge. it is also possible, that separate start points, through advantageous gain of function, can lead through selection, to convergence upon a common form and function. while it is practicable to assign family membership to genes or proteins based on their apparent hereditary relationships, the functional implications of family membership cannot simply be deduced. more than a decade of research effort has been aimed at understanding how syt works. considerable effort has also been directed at understanding the functions of other syt proteins. current understanding however, is insufficient to predict the functions of a whole genome complement of syt proteins, or indeed to predict what a genome complement of syt proteins might consist of. to help improve this lack of understanding, i have put together an open access resource intended to serve future research aimed at understanding the biological meaning of the hereditary and functional relationships among some of the metazoan syt-like proteins. i have collected and compiled manually, information about homologous syt, esyt, rph3a, rph3al, doc and dblc genes from metazoan genome sequences spanning a wide evolutionary range. i have not attempted to collect all genes capable of encoding proteins with similarity on any level to syt proteins. such a collection would ultimately include perhaps, all eukaryotic c domain proteins. my collection is restricted to homologous syt genes plus a small group of genes already identified as similar but not homologous to syt genes. the collection does include a new group of homologous genes which are present in the genomes of marine invertebrates and encode proteins similar to syt proteins. the information in this collection may interest investigators in the field of evolutionary biology, with examples of dynamic genome evolution including whole gene duplication, partial duplication, gene fission, acquisition of novel coding exons, gene inactivation, extreme sequence divergence, intron mobilization and a variety of routes to the expression of altered gene products. for experimentalists who want to understand how syt proteins act to effect the biology of animal nervous systems, the collection provides a utility with great evolutionary depth. in order to make progress in neuroscience, appropriate model organisms must be chosen. the choice of model organism is fundamental to the type of science which can be pursued. the range of organisms and genes in this collection, offer novel possibilities for future descriptive, comparative and hypothesis driven research. in order to make the large amount of information described by this paper available for public scrutiny, it is available as additional files which can be accessed online.  origin of this gene collection the origin of this manual gene curation project lay in an attempt to annotate some of the plant ntmc genes i had identified previously. i chose the unique ntmc nomenclature to emphasise the phylogenetic relationships among these plant genes and to distinguish them from animal syt genes. plant ntmc2t ntmc2t and ntmc2t genes resemble animal syt genes in having the same domain architecture. that this shared domain architecture implies functions in common, has now been demonstrated experimentally. the a. thaliana ntmc2t protein possesses calcium and membrane binding activities which allow it to function in a plasma membrane repair pathway induced by stress. while using the annotation facilities at doe-jgi i realised that the evolutionary range of genome sequences available for inspection at doe-jgi, could be used to search for earlier origins of the embryophyte ntmc genes and metazoan syt genes. the genome sequence of m. brevicollis has demonstrated that, as theorised the evolutionary transition from unicellularity to multicellularity involved abundant domain shuffling in proteins. m. brevicollis, a unicellular protist, closely related but basal to metazoans, does not possess syt genes. it does have a gene with a degree of gene structure similarity, domain architecture similarity and amino acid sequence similarity to plant ntmc and animal esyt proteins xm_ proteins with similarity to ntmc and esyt proteins are present in other unicellular eukaryotes. it seems reasonable therefore, to propose that while ntmc and esyt genes could ultimately be homologous, having attained their current forms through descent from a unicellular ancestor, syt genes, distinguished by their conserved gene structure, more likely represent a metazoan novelty. on finding uniquely metazoan syt homologues, with their characteristic gene structures, in the genomes of t. adhaerens and n. vectensis i decided to compile afresh, the syt repertoire of the currently available genomes of evolutionary significance. in addition to syt genes, i set out to gather esyt, rph3a, rph3al and doc homologues. these gene families, while clearly distinct from the syt family at the gene structure level, encode proteins with high level sequence similarity to the tandem c domain region of syt or with similarity, in common with ntmc proteins, at a more crude domain architecture level. because of these similarities in form, there will be common functions. a collection like this can be used to further define these functions, by revealing the patterns of sequence conservation and sequence divergence, apparent in: the individual domains of an orthologous group of proteins, the equivalent domain in different groups, not necessarily homologous. such knowledge can provide clues about how these proteins function as molecular machines in cells and about how these functions have developed over evolutionary time. i used the annotation facilities at doe-jgi to create and annotate gene models where possible. full details of all genes are listed in additional files and figures and summarise the genes in additional files and respectively. gene identification i searched for homologues among the metazoan genome sequences at doe-jgi ncbi and ucsc. i employed a variety of methods, first to identify, then to piece together, gene sequences for homologues of syt, esyt, rph3a,rph3al and doc these are generally large, multi-exon genes which encompass regions of high and low levels of sequence conservation. as a first step, i used blatp or tblastn with amino acid sequence probes representing already known homologues, to identify genomic loci with potential to encode a similar protein. i then examined a 3-frame translation of the genomic region to look at the possible exon-intron structure. intron positions and phases are generally well conserved among these genes and serve as a useful diagnostic tool. in cases where the draft genome sequence was highly fragmented, homologues were pieced together exon by exon. where sequence conservation between amino acid probe and novel genomic locus was high, it was straightforward to identify exons and introns. where sequence conservation was low, the conservation of intron position and phase could serve as a guide to judge possible exon boundaries. in all cases, when transcript sequences from the locus under examination, or from a homologous locus in a closely related organism, were available, these were used to help identify exons. occasionally, trace archive reads were sought to confirm exon sequences, where there was a gap or where there appeared to be frameshift errors in a section of draft sequence for example. details of such instances are noted in additional files and where transcript sequences indicated the expression of alternative variants, i collected the variants. in the case of m. musculus and h. sapiens, where there is abundant transcript variation, particularly at certain genomic loci, i collected only those variants indicated by two or more transcript sequences. where evidence of alternative splicing exists but is not represented by full length transcripts the resultant gene predictions require validation, as do all gene predictions for which no transcript evidence exists. for genomic loci with conserved variant exons, which could be employed to extend the length of the encoded protein but where transcript evidence was lacking, the gene prediction includes all of the conserved exons. transcript evidence for the gene predictions is listed in additional files and once i had gathered a set of crude gene predictions, i compared the sequences in order to refine the predictions. exon-intron junctions, gaps, novel regions of sequence conservation useful as probes to fill gaps, were all carefully examined using multiple alignments. in this way, it was possible to extend and improve the gene predictions significantly. many of the gene predictions in this collection are still incomplete, at the n-terminal region of syt proteins in particular. when there was no transcript coverage and no detectable sequence conservation or very short coding exons, it was not possible to make a prediction. all incomplete predictions are noted in additional files and sequence comparison: global comparison and nomenclature having collected genes from genomes, i used clustalw at ebi to compare subsets of amino acid sequences, in order to examine their relationships of similarity overall. i compared sets which excluded expressed variants but which included one sequence per gene, whether complete or not. comparisons of invertebrate sequences vertebrate sequences invertebrate plus mammalian sequences and invertebrate plus a subset of vertebrate sequences were used to examine relationships and assign a suitable name to each gene. i assigned a common name to those genes which clustered together in conserved groups. when a conserved group included a homologue from m. musculus/h. sapiens, the assigned name was chosen to be consistent with that in m. musculus/h. sapiens. for those syt genes which did not consistently cluster together and are not conserved among different organisms, i assigned them numbers upward from the last conserved group. this process led to the identification of conserved groups of syt genes plus a further unique syt genes named with the syt stem symbol plus a number identifier, in order to provide a recognisable name for each. this naming system allows further sampling of the tree of life by gene sequencing to reveal: more unique syt genes, from syt upwards, all members of a recognisably homologous family and previously unrecognisable patterns of conservation among syt genes, from syt upwards. i contacted the relevant organism nomenclature authorities to discuss this nomenclature. this wide consultation led to agreement to change fam a temporary hgnc symbol for this poorly characterised gene family, to esyt. where current approved nomenclature is discordant with that proposed here, the approved nomenclature is indicated within brackets, in additional files and figures and there is inconsistency in the application of my naming rationale in drosophila and caenorhabditis, where it was more practical to stick with the approved nomenclature than propose new names. as yet, and in comparison to vertebrates, invertebrates are poorly represented by organised gene nomenclature authorities. my proposed nomenclature for a large number of genes in this collection, therefore, awaits sanction. note that due to differing nomenclature conventions, the gene symbols in different organisms, range from having no capital letters to all capital letters. the generic form, which i have used to label the genes from organisms outside the purview of nomenclature authorities, is first letter in capitals. this designation will be used here, to describe the gene in general, without reference to a specific organism. deciding whether relationships of orthology or some other kind exist among a group of similar genes, can be complicated. the processes of duplication, deletion, and rearrangement of genes, plus the action of selection or absence of selection on gene sequences, can obscure relationships. whether absence from a genome assembly reflects true absence, or incompleteness of the assembly, presents a further problem of some seriousness. in m. musculus and h. sapiens there are several groups of duplicated syt genes: group group group and group. each group contains genes which share a distinct gene structure. additional files and figure indicate that the parental gene of group is syt the parental gene of group is syt that of group is syt and that of group is syt in each of these cases, the parental gene is more similar to a gene present in an evolutionarily more ancient genome than it is to another member of the mammalian group, which is to say, orthologous. orthologous genes in different organisms, are related by vertical descent from a common ancestor. relevant pairwise comparison scores are listed in figure it is worthwhile trying to determine the phylogenetic relationships between the syt genes of mammals and those of other animals because it is the rodent syt proteins which are best characterised functionally. however, because mammalian organisms and their brains in particular, represent biological systems of such complexity, modelling these systems by using simpler systems is essential. the relationships and nomenclature proposed here, are in good general agreement with the syt data at treefam which does not include many of the invertebrate syt genes in this collection. this gene collection includes several previously unrecognised syt groups which are absent from m. musculus/h. sapiens. members of the syt group appear in vertebrates and invertebrates, but the syt group has yet to be found beyond vertebrates. neither of these groups are well conserved, both genes having diverged significantly over time. the syt group ranges from members encoding: tandem c domains a single c2b domain or no c domains at all. syt appears in only two hexapod species in this collection. the other syt groups absent in vertebrates, are sytalpha and syt of these, the sytalpha group has a higher level of sequence conservation. the sytalpha designation applies to all members of the conserved group which contains the d. melanogaster sytalpha gene. there are many invertebrate syt genes, related to some degree, at the gene structure level, which are not conserved among different phyla. n. vectensis, for example, has a very large number ofunique and divergent syt genes. most of these encode stereotypical syt proteins, but some have only one c domain and some lack a tm domain. the caenorhabditis syt repertoire, in common with many others in this collection, demonstrates that variant syt proteins which lack the tm domain, are usually included. they are encoded either by separate genes or by alternative splicing. caenorhabditis snt- is very similar to snt- at the amino acid sequence level, but does not encode a tm domain. in the ciona syt repertoire, it is syt which produces a variant lacking the tm domain, but by alternative splicing. in mammals, syt does not encode a tm domain, but a number of other mammalian syt genes which do encode a tm domain, express splice variants which lack the tm domain. transcript variants are listed in additional files and sequence comparison: gene duplications the duplicate genes in the fish genomes presented a nomenclature challenge and i am very grateful for help from the zfin nomenclature coordinator in choosing the most appropriate names for these genes based on gene structure similarity and synteny. a number of genes present as a single copy in mammals are present as duplicates in the fish genomes. in addition, syt and syt are duplicated in d. rerio. it is notable that the single rph3a of tetrapods appears to have undergone duplication and fission, giving rise to rph3al and doc2b genes. rph3al and doc2b are linked as a tandem pair in tetrapod genomes. invertebrates have a single, equivalent rph, but no doc genes. it is possible that d. pulex rphl represents an n-terminal rph fission product. in the fish genomes, there are duplicate rph3a genes, but only in the lineage leading to d. rerio, has a supposed rph3a duplication led to the retention of linked rph3al and doc2b fission products. in the other fish genomes, doc2b genes are present, along with the paralogous doc2a and doc2d genes, as well as duplicate rph3aa and rph3ab genes, so doc2b genes are not always linked to a rph3al fission product. apart from the gene duplications already noted, other paralogous gene expansions are evident in this collection. in the case of esyt gene duplications, retroposition of processed transcripts has occasionally led to gene duplication. sequence comparison: multiple alignments in order to illustrate the conserved features of the proteins in this collection, i used multalin to construct multiple alignments of subsets of sequences, including expressed variants. the conserved groups of vertebrates and invertebrates are compared separately but the syt group contains members from vertebrates and invertebrates. all of the syt sequences from n. vectensis are compared. these multiple alignments illustrate the patterns of sequence conservation distinctive of different groups. common patterns of alternative splicing which serve to alter particular domains, are also illustrated. where calcium coordinating residues are absolutely conserved these are indicated by arrows. these absolutely conserved sequences are highlighted in red in figures and in the many more cases where calcium coordinating positions are not absolutely conserved, there is a range, down to zero conservation, of the ten sites present in the absolutely conserved, calcium responsive syt proteins. a conserved region 20- residues upstream of the c2a domain, often recognisable as a variant of this motif: grikpely is indicated where present. conservation of this motif was noticed previously but it is obvious here, that this region is much more widely conserved. its function is yet to be investigated. the position and phase of introns are also illustrated, in order to show the gene structures distinctive of different family groups, as well as instances of intron mobilization. intron mobilization within a genus is specifically noted in additional file intron phase is indicated by colour: black for phase red for phase + and blue for phase + syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. the ten conserved calcium coordinating positions are indicated by arrows. the c2ab region is highly conserved whereas the n-terminal portion has much more sequence variation. expressed variants are included to demonstrate the regions altered by alternative splicing and rna editing. the region between the tm domain and the c2a domain is commonly altered, for instance in s. purpuratus, b. floridae, a. gambiae and d. melanogaster. the next region expressed as alternative forms, is the region between c2a and c2b, as indicated by d. pulex and a. gambiae variants. the third region to be altered, by alternative splicing or by rna editing, is the central portion of the c2b domain. in caenorhabditis, this recoding alters not only the identity of positions within this exon of amino acids length in all other cases, but extends its length by or amino acids. additional file also includes the snt- sequences from caenorhabditis. the snt- gene encodes a protein very similar to the snt- protein which is the syt orthologue. the snt- proteins however, lack the n-terminal membrane anchoring portion. the snt- proteins also differ in the region between c2a and c2b. additional file shows a multiple alignment of the vertebrate syt protein sequences in this collection. the ten conserved calcium coordinating positions are indicated by arrows. the c2ab region is highly conserved whereas the n-terminal portion has a little more sequence variation. expressed variants are included to demonstrate the regions altered by alternative splicing. this time, the alteration is restricted to the region between the tm and c2a domain. conserved n-glycosylation consensus sites are indicated by blue boxes. in the first coding exon, the n-glycosylation site and the conserved upstream o-glycosylated threonine, are known to be important functional sites in vivo. additional files and show multiple alignments of the vertebrate syt syt and syt protein sequences, respectively. for syt syt and syt genes, there is transcript evidence of alternatively expressed forms which lack the tm domain. syt syt syt and syt are all related by duplication, with syt being the parental gene. in the fish genomes, there are yet more duplicates. in order to distinguish the relationships among these duplicates, pairwise comparison scores and syntenic relationships were examined. the top scores highlighted in blue in figure indicate that d. rerio syt1a is orthologous to the syt gene of other vertebrates and that d. rerio syt is orthologous to the syt gene of other vertebrates. the pairwise comparison scores do not however, indicate orthology for the other duplicates. figure shows that relationships of orthology, are indeed present. the neighbouring genes of syt syt syt and syt in the genomes of h. sapiens, x. tropicalis and d. rerio, reveal conserved syntenic blocks in each genome. this synteny indicates that the four genomic segments share ancestry. orthologous relationships exist between: h. sapiens syt and the duplicated syt1a/syt1b in d. rerio, between syt and syt2/syt2b and between syt and syt5a/syt5b, respectively. since there is so much interest in the mechanism of syt function, close examination of the syt orthologues and paralogues in this collection, to see how they have been shaped by evolution, may prove productive. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. the variants of d. pulex differ at the c-terminal end and those of a. gambiae differ at the n-terminal end. additional file shows a multiple alignment of the vertebrate syt protein sequences. the variants of g. gallus differ in the region between c2a and c2b. additional file shows a multiple alignment of the protein sequences of syt the vertebrate paralogue of syt here, the variants in g. gallus, m. musculus and h. sapiens differ in the c2b domain. recent functional studies of mouse syt have built on initial gene knockout results which suggested that syt dysfunction could be responsible for some human psychiatric disease. syt dysfunction has also been implicated in human psychiatric disease. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. here, there is evidence from d. melanogaster, of variant forms lacking the tm domain. additional file shows the vertebrate syt sequences. in comparison to the invertebrate syt genes, the vertebrate syt genes contain additional coding capacity in the region between the tm and c2a domains. several additional exons which can be alternatively spliced, are available to extend the length of the region linking the tm and c2a domains. a study of constitutive versus alternative splicing, has examined the syt alternative exons. recent gene knockout studies in mice appear to confirm a role for syt in glucose homeostasis. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. the conserved motif upstream of the c2a domain is indicated. additional file shows the vertebrate syt sequences. the variants expressed by m. musculus and h. sapiens differ at one or other end. additional files and show the protein sequences of the vertebrate syt paralogues, syt syt and syt respectively. all of these proteins retain all ten calcium coordinating residues and are thus likely to have been selected to act as calcium sensors. they also retain the motif of unknown function upstream of the c2a domain. transcript evidence from x. tropicalis, m. musculus and h. sapiens indicates the expression of alternative forms of syt which lack the tm domain. the c-terminal region of syt is also expressed as alternative forms. functional studies have examined the conserved cysteine residues of the n-terminal region which are unique to this paralogous family in vertebrates. syt is thought to play a role in sperm acrosomal exocytosis. mouse syt has been identified as upregulated during meiosis in fetal ovaries. note that some investigators working with syt proteins, use different synonyms for the protein products of certain syt genes. the officially approved gene nomenclature is used throughout this paper. in the functional literature, the protein products of syt and syt have most often been referred to differently. syt and syt so far poorly characterised, have also been prone to a similar lack of recognition. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. while there is some variation in length, of sequences from different organisms in the c2b region, the more pronounced difference is a trend towards increased length of the region between tm and c2a domains in insects. additional file shows the vertebrate syt sequences. transcript evidence indicates the expression of forms varying in the n-terminal region in h. sapiens. little functional work has been done so far with syt. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. in most of these genes, there are two alternative n-terminal coding exons, only one of which encodes a tm domain. these sequences also retain the motif of unknown function just upstream of the c2a domain. additional file shows the vertebrate syt sequences. transcript sequences from h. sapiens indicate the expression of a form lacking the tm domain. a recent study suggests that human syt may function as a tumour supressor. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. in these sequences, it is the c2a domain which is more highly conserved than the c2b domain. the motif upstream of the c2a domain is present. there is a large amount of variation in amino acid sequence, intron position and intron phase in the n-terminal portion of these sequences. additional file shows the vertebrate syt sequences. transcript evidence from m. musculus and h. sapiens indicates the expression of variants differing at the c-terminal end. transcripts from h. sapiens indicate variant forms at the n-terminal end, which lack the tm domain. very little is known about the function of syt. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. transcript evidence indicates the expression of variants in c. savignyi, c. intestinalis and s. purpuratus which differ at the n-terminal end and lack the tm domain. additional file shows the vertebrate syt sequences which contain a tm domain. additional file shows the vertebrate syt sequences which lack a tm domain. additional file shows vertebrate syt sequences. transcript evidence in x. tropicalis, m. musculus and h. sapiens indicates the expression of variant forms of syt with altered c2b domains. human syt transcripts also differ at the n-terminal region, resulting in forms which lack the tm domain. loss of function of human syt has been associated with neurodevelopmental abnormalities. in the human genome, there is a repeat of a portion of the syt gene on chromosome apparently the result of retroposition of a syt transcript. the chromosome locus has been named syt14l by hgnc. syt14l is classed as a pseudogene by ncbi and is not included in this collection. syt homologues additional file shows a multiple alignment of the invertebrate syt protein sequences in this collection. a possible motif upstream of the c2a domain is indicated. additional file shows the vertebrate syt sequences. transcript evidence in m. musculus and h. sapiens indicates the expression of variants altered at n-terminal and c-terminal ends. a possible motif upstream of the c2a domain is indicated. the protein product of syt has been referred to as b/k protein, but not much is known about its function. sytalpha homologues additional file shows the sytalpha sequences in this collection. a possible motif upstream of the c2a domain is indicated. the only functional study relevant to sytalpha is its anatomical localisation in d. melanogaster. syt homologues additional file shows the syt sequences in this collection. these proteins are so far completely uncharacterised. syt homologues additional file shows the vertebrate syt sequences in this collection. these proteins are so far completely uncharacterised. syt homologues additional file shows the invertebrate syt sequences in this collection. these proteins are so far completely uncharacterised. syt homologues in n. vectensis additional file shows all of the syt proteins of n. vectensis, demonstrating enormous sequence variation. these proteins are so far completely uncharacterised. dblc homologues additional file shows the dblc sequences in this collection. transcript variants are indicated. these proteins have not previously been reported and are so far completely uncharacterised. dblc sequences are detectably similar to syt sequences at the amino acid sequence level and occur in the genomes of marine invertebrates in this collection. while they lack an n-terminal tm domain, they have a tandem c domain architecture in common with syt proteins. dblc genes lack similarity to syt genes at the gene structure level, so a new gene symbol was chosen to represent these genes. the dblc designation represents a fuller description: double c domain protein. dblc is thus similar but recognisably different from doc the designation for the homologous family of double c domain proteins of vertebrates. esyt homologues additional file shows the invertebrate esyt sequences in this collection. transcript variants are indicated. the arthropod esyt genes include mutually exclusive alternative exons which serve to alter the specification of the first c domain. additional file shows the sequence of the first of these alternatives only. additional file shows all versions of this alternatively coded region. additional file shows the vertebrate esyt sequences. transcript variants are indicated. in the fish genomes, esyt has duplicated, resulting in esyt1a and esyt1b. the esyt1b duplicate has acquired an internal duplication which contributes an additional four c domains to the product of this gene. within this duplicated section, the pink dot marks the position of an apparent phase intron loss in the tetraodontiform fish. in each of the three other occurrences of this portion of the esyt1b gene, the phase intron is present. this intron is also present in each of the two repeats of this sequence, in the esyt1a genes. additional file shows the vertebrate esyt sequences. transcript variants are indicated. in the fish genomes, esyt has duplicated, resulting in esyt2a and esyt2b. additional file shows the vertebrate esyt sequences. transcript variants are indicated. similarity between syt genes and what are now named esyt genes, was first evident from genome sequence comparisons. in trying to classify and annotate the genes within the human genome, hgnc noticed that these genes formed a distinct gene family, separate from syt genes. as there was little functional data associated with these genes, a temporary nomenclature with no functional implication was provided: fam62a, fam62b and fam62c. the functional data which did exist, was in the form of a cdna representing rat esyt which had been cloned during an investigation of adipocyte proteins and found to encode a membrane bound c domain protein with similarity to proteins in plants. further genome sequence comparisons indicated that syt, fam and plant ntmc genes encode proteins with amino acid sequence similarity and domain architecture similarity, but belong to distinct gene families. apart from and until recently, functional data associated with plant ntmc and animal fam gene families had been lacking. an initial investigation of proteins encoded by the human fam62a, fam62b and fam62c genes, has now been published and the authors named these proteins esyt esyt and esyt to stand for extended synaptotagmin-like proteins. consultation with the nomenclature committees during the preparation of this paper, led to the decision to implement the esyt nomenclature in place of fam in an initial pair of studies on one member of the plant ntmc gene family the authors put forward the name syt for this particular gene and gene product. as discussed in the background section of this paper, when different communities make naming decisions based on the criteria most relevant to them it is inevitable that multiple names will be used to identify the same entities. rph3a homologues additional file shows the invertebrate rabphilin sequences in this collection. in caenorhabditis, alternative n-terminal exons can be used. sequence conservation among rabphilins is high in the n-terminal, rab binding portion and the c-terminal tandem c domain portion, but a large middle portion is very poorly conserved. this means that where transcript confirmation is absent, the gene predictions across this portion are unlikely to be accurate. additional file shows the vertebrate rph3a sequences. transcript variants are indicated. in vertebrates, all calcium coordinating positions are absolutely conserved, but not all are conserved in invertebrates. in additional file the positions of these amino acids are indicated by pink dots. the position of the fifth calcium coordinating residue in the c2a domain is not strictly conserved, but in most cases a suitable aspartate or glutamate residue is present one residue earlier. i am not aware of biochemical evidence for calcium binding by invertebrate rabphilin proteins, but genetic evidence suggests that c. elegans rbf- at least, appears to function similarly to its mammalian homologues. for this reason, in figures and i have counted the fifth position as present when a suitable aspartate or glutamate residue occurs at the usual spacing or one residue earlier. rph3al homologues additional file shows vertebrate rph3al sequences. transcript variants are indicated. in the functional literature, products of rph3al genes have been referred to as noc. as outlined above, rph3al genes appear to be the result of duplication and fission of an ancestral rph3a gene, producing linked rph3al and doc2b genes in tetrapods and in d. rerio. rph3al represents the n-terminal portion of the ancestral rph3a gene and doc2b represents the c-terminal portion, as seems clear from a comparison of the gene structures of rph3a, rph3al and the linked doc2b genes. as also indicated above, doc2b genes do not always accompany rph3al genes, as they exist in other fish genomes where rph3al genes do not occur. doc homologues additional file shows the vertebrate doc protein sequences. transcript variants are indicated as well as the positions of the calcium coordinating residues, which are conserved in all of the doc proteins except doc2g. i have included potential doc2g products from the human gene locus. although spliced transcripts are produced from this locus, they cannot encode full length proteins and it is not known whether protein products are translated from these transcripts. it is clear from the shared gene structure of doc genes and rph3a genes, that doc genes represent the c-terminal portion of a rph3a gene. since rph genes exist in invertebrates, but doc genes do not, it seems reasonable to suggest that doc genes arose via duplication and fission of a rph3a gene. evidence for this, in the form of the two linked fission products is present in the genomes of d. rerio, x. tropicalis, a. carolinensis, g. gallus, t. guttata, o. anatinus, m. domestica, m. musculus and h. sapiens. in this collection, the fish genomes reveal a novel lineage of doc genes, which in consultation with zfin staff were given the gene symbol doc2d, to stand for doc delta. the first doc gene, cloned in was soon found to represent the first member of a family of three genes in rodents. the protein products of these three genes were named doc alpha, doc beta and doc gamma. the conserved region present in the n-terminal portion of doc sequences, specifies the munc binding capacity of doc proteins.  this paper describes a manually curated collection of genes pertinent to intercellular trafficking in multicellular animals. the collection includes genes from the genome sequences of marine invertebrates genes from the genome sequences of ecdysozoans, including the drosophila and caenorhabditis genomes plus genes from vertebrate genomes. when viewed overall the collection implies that particular homologues which are always present in the genome repertoire, serve to specify functions required for animal life. the basis for this view is that these genes have not been lost, but have been retained and conserved, over the long period of metazoan evolution. such a view about the essentiality of gene function, differs from a synthetic biology view, which considers the genes necessary to specify an operable system which functions outside the natural world. it differs too, from a reductionist molecular biology view where, if a gene can be deleted in the laboratory setting, without killing the organism, it is not considered essential. such gene deletion experiments have been done with syt in different organisms and according this view, syt is not strictly essential to animal life because snt- null mutants are not lethal in caenorhabditis. in drosophila and mus however, syt null mutants are indeed lethal as they cannot survive to reproduce as adults. while my previous genomic comparison of syt genes could have been interpreted as indicating a correlation between increased numbers of syt genes and increased organism complexity the analysis here, which is much more comprehensive, definitively rules that conclusion out. this collection of genes offers prospects for fruitful speculation about the functional roles of the syt repertoires of different animals and includes a great range of biological complexity. the conclusions drawn from evolutionary genomics approaches are always provisional, because they depend on the quality of information available and are thus subject to review when better information becomes available. conclusions drawn from reassembled and reannotated genomes will inevitably be more accurate than those drawn from the first draft. the d. rerio genome being repeat rich, is taking time to assemble into a complete version. among the genomes in this collection, the draft d. rerio genome ranks first in terms of gene number. d. rerio surely presents a useful model system with which to investigate the functional divergence of close paralogues. currently, it is the rodent syt proteins and, with the advent of transgenic mouse studies, their respective genes, which garner most interest and have been best characterised. work with syt genes in other animal species has provided information supporting hypotheses about how syt genes fulfil their functional roles. to date, these studies indicate a primary role for syt in specifying a neuron specific synaptic vesicle protein which serves to trigger exocytosis. the genomic specification of syt expression patterns, has been investigated and utilised in ascidian species. in c. elegans, the anatomical expression patterns of the mutually exclusive alternative exon variants of syt have been mapped. this gene collection shows that such mutually exclusive exon deployment has arisen on a number of occasions in invertebrates, altering c domains in syt syt and esyt proteins. the repeated evolution of this type of c domain alteration, suggests that careful comparison of the altered forms might prove a useful focus for investigations aimed at understanding how these c domains operate as molecular machines. the feature most widely shared among the great diversity of syt proteins in this collection, is sequence conservation of the c-terminal c2b region, which suggests that a defining function of syt proteins is located here. the regulatory machinery controlling variant expression can be illuminated by comparative analyses of genomic sequences, as was done to investigate the rna editing of syt which occurs in hexapods. another comparative analysis examining the genomic specification of alternative splicing of syt concluded that splice variants expressed in d. melanogaster, c. pipiens, b. mori, t. castaneum and p. humanus depend on correctly coded intronic, cis-acting regulatory sequences acting on the requisite splice donor/acceptor sequences. in this collection, there is evidence from a. gambiae, d. melanogaster, d. pulex, s. purpuratus, b. floridae, d. rerio, x. tropicalis, a. carolinensis, g. gallus, m. musculus and h. sapiens of rna editing and alternative splicing of syt. these sequences may aid future efforts to identify the cis-sequences involved in transcript processing. the biochemical understanding of transcript production and processing, acquired at great pains over many years is not quite up to the task of extracting the full meaning from the flood of new transcript data from high throughput sequencing. an enormous amount of experimental work will need to be done to understand the varied roles of transcripts in complex organisms. the problem of identifying and annotating different kinds of transcript sequences is discussed in. work to assess the functional significance of alternatively expressed syt genes has barely begun so it is difficult at present, to judge the importance of the transcript variants listed in this collection. evolutionary conservation of variant production, probably argues for functional significance. new developments in mass spectrometry are bound to aid attempts to examine complex mixtures of proteins such as synaptic vesicle preparations. it is more likely however, that the functional impact of syt protein variants will eventually be revealed by research focussed on syt proteins of known importance. the significance of the conservation of intron position and phase within syt genes could be investigated by gene targeting in transgenic mice, now a common route to investigate syt roles. the best conserved intron, a signature feature of syt genes, is the phase + intron at the start of the third beta strand of the c2b domain. a single point mutation to inactivate this intron acceptor dinucleotide would be predicted to have devastating consequences for the expression of the host gene. in contrast to the approach of disabling a gene by creating a large deletion, the production of a null mutant by such a small intervention would be unlikely to cause collateral damage. the chances of cryptic splice acceptors being deployed to restore normal function to such a targeted syt splice site point mutant, seem remote, based on the importance of this region to syt protein function. the presence of this intron in so many syt genes could be interpreted as evidence that evolutionary alteration of this site is hazardous. an example of a gene inactivating mutation which is not rescued by altered splicing, is present in this gene collection. the h. sapiens gene which is orthologous to m. musculus doc2g has acquired a frameshift mutation in the fifth coding exon which introduces a stop codon. human transcript sequences indicate that an alternative intron acceptor dinucleotide capable of restoring the reading frame, which exists nucleotides upstream of the exon acceptor, is not used. the regulatory sequences necessary to select this alternative acceptor are not present within this short intron. it is not known what protein products, if any, are produced by this gene in h. sapiens and it appears likely that the human gene is non functional. this observation underscores the importance of splicing signals which are poorly understood at present, as well as the importance of transcript analysis to verify gene predictions. in sea urchin species, the syt protein has been shown to act in cortical granule exocytosis. syt products are used as pan neuronal markers in echinoderms, hemichordates and urochordates and the role of syt in synaptic vesicle exocytosis is assumed. few comparative studies have attempted to map the expression patterns of the syt repertoire of a particular organism but all of the syt genes of d. melanogaster and a large fraction of those of r. norvegicus have been examined in this way. t. adhaerens is thought to lack cells equivalent to neurons, but it has a well conserved syt as well as genes encoding the requisite neuronal snares and complexin for fast neurotransmission. although never observed in the laboratory, t. adhaerens is expected to produce gametes in the wild so a role for the t. adhaerens syt similar to that of sea urchin syt in cortical granule exocytosis is possible, alongside any yet to be discovered'neural' role. n. vectensis, which has a bona fide nervous system, also has syt and the requisite neuronal snares and complexin. n. vectensis is quite unusual among the metazoans collected here, in having so many highly divergent syt genes. of all the n. vectensis syt genes, only syt and syt encode all ten calcium coordinating positions, indicating their likely role in calcium triggered membrane fusion. future studies to determine the anatomical expression patterns of the n. vectensis syt genes in vivo, will be of great interest to further understand the evolutionary development of animal nervous systems. initial study of the molecular anatomy of the n. vectensis nervous system has examined the expression patterns of neurotransmitters, neuropeptides, transcription factor and cell adhesion genes. synaptotagmin gene expression was used to label differentiated neurons in platynereis, in a study of the evolutionary development of the central nervous system of animals. capitella, like platynereis, is an annelid member of the lophotrochozoa, but has the advantage as a model organism, of a whole genome sequence. study of neurogenesis in capitella is now underway. in situ hybridization analysis of syt expression in capitella, shows a broad neural pattern. it will be important in future, to describe and compare anatomical expression patterns and functions of syt genes, both conserved and unique, in vertebrates and invertebrates. such comparisons will be essential to improve, or reduce confidence in functional annotation of genomes based on phylogenetic homology, or similarity otherwise defined. good annotation of genomes is a necessary step in the process of biocuration a new field aimed at using twenty first century sequencing and information technology to make accessible to bench scientists, the vast amount of information they require, distilled, in order to take forward wisely informed hypothesis driven research. community deposition of a variety of types of biological data into appropriate databases, database efforts to curate, update and integrate this data with relevant data in other databases, will if pursued successfully, lead to great progress in biology in the years ahead. for the time being, the manual gene collection made available here, can be used immediately, by bench scientists engaged in hypothesis driven investigation of membrane trafficking by these proteins, to widen the bounds of a field already stretching from biophysical study in vitro at the atomic level, to phenotypic assessment of gene function in model organisms. authors' contributions the author had sole responsibility for the design and execution of this work. supplementary material additional file 1 comprehensive information about the marine invertebrate genes in this collection. this spreadsheet file contains full details of each gene identified in the marine invertebrate genomes in this collection. click here for file additional file 2 comprehensive information about the ecdysozoan genes in this collection. this spreadsheet file contains full details of each gene identified in the ecdysozoan genomes in this collection. click here for file additional file 3 comprehensive information about the vertebrate genes in this collection. this spreadsheet file contains full details of each gene identified in the vertebrate genomes in this collection. click here for file additional file 4 comprehensive information about the snare and complexin genes identified in t. adhaerens and n. vectensis. this spreadsheet file contains full details of the genes identified in t. adhaerens and n. vectensis which are similar to neuronal snares and complexin. click here for file additional file 5 dendrogram of relationships among the invertebrate sequences in this collection. additional file displays the guide tree of the clustalw comparison of the invertebrate sequences in this collection, excluding variants, totalling sequences. genes which encode mutually exclusive alternative exons are highlighted with a green box. conserved groups of syt genes which have not previously been described, are highlighted with a blue box. click here for file additional file 6 dendrogram of relationships among the vertebrate sequences in this collection. additional file displays the guide tree of the clustalw comparison of the vertebrate sequences in this collection, excluding variants, totalling sequences. mouse genes are highlighted with a red box. conserved groups of syt genes which have not previously been described, are highlighted with a blue box. click here for file additional file 7 dendrogram of relationships among the invertebrate and mammalian sequences in this collection. additional file displays the guide tree of the clustalw comparison of the invertebrate sequences excluding variants, plus mammalian sequences, excluding variants, totalling sequences. mouse genes are highlighted with a red box. conserved groups of syt genes which have not previously been described, are highlighted with a blue box. click here for file additional file 8 dendrogram of relationships among the invertebrate and vertebrate sequences in this collection. additional file displays the guide tree of the clustalw comparison of the invertebrate sequences excluding variants, plus a subset of vertebrate sequences, excluding variants, totalling sequences. mouse genes are highlighted with a red box. conserved groups of syt genes which have not previously been described, are highlighted with a blue box. click here for file additional file 9 alignment of the invertebrate syt sequences, plus caenorhabditis snt- amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 10 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. conserved n-glycosylation consensus sites are indicated by blue boxes. the conserved threonine, which can be o-glycosylated is also indicated by a blue box. x residues indicate where a portion of sequence is missing. click here for file additional file 11 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. the x. tropicalis splice variant is included and highlighted with a black dot where it differs. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 12 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 13 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 14 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 15 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. the g. gallus splice variant is included and highlighted with a black dot where it differs. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 16 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 17 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. the d. melanogaster splice variant is included and highlighted with a black dot where it differs. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 18 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 19 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 20 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 21 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 22 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 23 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. click here for file additional file 24 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. because of their differing positions, tm domains are highlighted in blue. x residues indicate where a portion of sequence is missing. click here for file additional file 25 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. the h. sapiens splice variant is included and highlighted with a black dot where it differs. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 26 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants which specify an alternative n-terminus lacking a tm domain, are indicated on top, with a dotted line indicating where they join the common c-terminal portion. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 27 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. the h. sapiens splice variant is included and highlighted with a black dot where it differs. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 28 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. because of their differing positions, tm domains are highlighted in blue. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 29 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 30 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants which specify alternative n-termini lacking tm domains, are included and highlighted with black dots where their sequences differ. tm domains are highlighted in blue. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 31 alignment of the vertebrate syt sequences which have tm domains. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 32 alignment of the vertebrate syt sequences which lack tm domains. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 33 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 34 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 35 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. a possible motif just upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 36 alignment of the invertebrate sytalpha sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the widely conserved motif of unknown function, just upstream of the c2a domain, is indicated. x residues indicate where a portion of sequence is missing. click here for file additional file 37 alignment of the syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 38 alignment of the vertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. the x residue indicates where a portion of sequence is missing. click here for file additional file 39 alignment of the invertebrate syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. tm sequences are highlighted in blue. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 40 alignment of the n. vectensis syt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. tm sequences are highlighted in blue. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 41 alignment of the invertebrate dblc sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. click here for file additional file 42 alignment of the invertebrate esyt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 43 alignment of the mutually exclusive alternative esyt exons. in addition to the alternatively coded exon which is highlighted, the exon upstream and the exon downstream are shown. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. blue bars indicate phase + introns. click here for file additional file 44 alignment of the vertebrate esyt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. the middle portion is only present in the fish esyt1b sequences. a pink dot in this portion, marks an intron loss in t. rubripes and t. nigroviridis. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 45 alignment of the vertebrate esyt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 46 alignment of the vertebrate esyt sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. some intron positions are marked with dotted lines. in these cases, transcript sequence covers a gap in the genomic sequence and the intron presence is assumed. x residues indicate where a portion of sequence is missing. click here for file additional file 47 alignment of the invertebrate rabphilin sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are indicated. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by pink dots at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 48 alignment of the vertebrate rabphilin sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file additional file 49 alignment of the vertebrate rph3al sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. x residues indicate where a portion of sequence is missing. click here for file additional file 50 alignment of the vertebrate doc sequences. amino acid position is marked every hundred amino acids approximately, at the top of each page of the alignment. splice variants are included and highlighted with black dots where they differ. intron position and phase is indicated with a coloured bar between amino acids. black bars indicate phase introns. red bars indicate phase + introns. blue bars indicate phase + introns. some intron positions are marked with dotted lines. in these cases, transcript sequence covers a gap in the genomic sequence and the intron presence is assumed. the five conserved acidic amino acids in each c domain are indicated by black arrows at the top of the alignment. x residues indicate where a portion of sequence is missing. click here for file acknowledgements this work was supported by the uk medical research council. i wish to thank hugh pelham for making it possible for me to complete this work.
the evolutionary range of genome sequences now available, can shed more light on synaptotagmin gene phylogeny and clarify family relationships. synaptotagmin proteins were first identified in nervous tissue, residing in synaptic vesicles. only a subset of synaptotagmin proteins appear able to function as calcium triggers. esyt genes were formerly named fam the species in this collection are trichoplax, nematostella, capitella, helobdella, lottia, ciona, strongylocentrotus, branchiostoma, ixodes, daphnia, acyrthosiphon, tribolium, nasonia, apis, anopheles, drosophila, caenorhabditis, takifugu, tetraodon, gasterosteus, oryzias, danio, xenopus, anolis, gallus, taeniopygia,ornithorhynchus, monodelphis, mus and homo. this collection of genes offers prospects for fruitful speculation about the functional roles of the synaptotagmin repertoires of different animals and includes a great range of biological complexity. all of the data described in this paper is available as additional files. with reference to this gene collection, functional relationships among syt, esyt, rph3a, rph3al, doc and dblc genes, which encode similar proteins, can better be assessed in future. current understanding does not explain why this is so. further difficulty arises as families are defined quite legitimately in different ways: by common descent or by common possession of distinguishing features. these members have been investigated intensively, but other family members are not well understood, making it difficult to grasp the meaning of family membership in functional terms.
 synaptotagmin proteins were first identified in nervous tissue, residing in synaptic vesicles. synaptotagmins were subsequently found to form a large family, some members of which play important roles in calcium triggered exocytic events. these members have been investigated intensively, but other family members are not well understood, making it difficult to grasp the meaning of family membership in functional terms. further difficulty arises as families are defined quite legitimately in different ways: by common descent or by common possession of distinguishing features. one definition does not necessarily imply the other. the evolutionary range of genome sequences now available, can shed more light on synaptotagmin gene phylogeny and clarify family relationships. the aim of compiling this open access collection of synaptotagmin and synaptotagmin-like sequences, is that its use may lead to greater understanding of the biological function of these proteins in an evolutionary context.  metazoan genomes were examined and their complement of syt, esyt, rph3a, rph3al, doc and dblc genes identified. all of the sequences were compared, named, then examined in detail. esyt genes were formerly named fam the species in this collection are trichoplax, nematostella, capitella, helobdella, lottia, ciona, strongylocentrotus, branchiostoma, ixodes, daphnia, acyrthosiphon, tribolium, nasonia, apis, anopheles, drosophila, caenorhabditis, takifugu, tetraodon, gasterosteus, oryzias, danio, xenopus, anolis, gallus, taeniopygia,ornithorhynchus, monodelphis, mus and homo. all of the data described in this paper is available as additional files.  only a subset of synaptotagmin proteins appear able to function as calcium triggers. syt syt and syt are ancient conserved synaptotagmins of this type. some animals carry extensive repertoires of synaptotagmin genes. other animals of no less complexity, carry only a small repertoire. current understanding does not explain why this is so. the biological roles of many synaptotagmins remain to be understood. this collection of genes offers prospects for fruitful speculation about the functional roles of the synaptotagmin repertoires of different animals and includes a great range of biological complexity. with reference to this gene collection, functional relationships among syt, esyt, rph3a, rph3al, doc and dblc genes, which encode similar proteins, can better be assessed in future. 