<title>Background</title> <p>The recent shift of genomic biology from conventional model organisms to evolutionarily relevant species has led to the questioning of numerous ideas about metazoan evolution. For instance, the recently released genome of the starlet anemone has revealed a striking conservation with its vertebrate counterparts despite an apparent morphological gap between these organisms [##REF##17615350##1##]. On the contrary, whereas the Hox gene clusters have been considered for a long time as structures strictly required for the development of the common bilaterian body plan, they were found to be disorganized or even dislocated in animals such as nematodes or urochordates [##REF##17553908##2##,##REF##15343333##3##]. These cases illustrate the interest of genomic insights from organisms that display either peculiar morphological characteristics or have key phylogenetic positions.</p> <p>Interestingly, chætognaths, also known as arrow worms, fulfill both of these criteria: they have one of the most intriguing sets of morphological and developmental characteristics among animals and their phylogenetic position was recently reevaluated as a pivotal one for the understanding of animal evolution [##REF##16890517##4##]. These free-living marine creatures represent one of the major predators of the zooplancton food-chain but the phylum is mainly known for its original mosaic of morphological characteristics that have puzzled zoologists for years [##UREF##0##5##]. Their nervous system exhibits typical protostome features, such as ventral nervous mid-body ganglions and circum-esophageal fibers [##REF##17511857##6##], whereas the enterocoelous formation of their body cavity and the secondary emergence of their mouth are embryological features traditionally related to deuterostomes [##UREF##1##7##]. Strikingly, this original body plan has been conserved since the lowermost Cambrian period as shown by convincing fossil evidence [##REF##12364798##8##,##REF##17254986##9##]. First attempts to position chætognaths using molecular phylogeny were difficult because small subunits (SSUs) and large subunits (LSUs) of ribosomal RNA genes display very fast evolutionary rates that hinder accurate tree reconstruction [##REF##16434216##10##, ####REF##8336548##11##, ##REF##9069174##12####9069174##12##]. Subsequent analysis of their mitochondrial genome prompted classification of chætognaths among protostomes, but their exact branching in this clade remains elusive [##REF##15249679##13##,##REF##15306659##14##]. The Hox genes of chætognaths are distinct from those typical of other protostomes: their original <italic>MedPost </italic>gene shares similarity with both median and posterior classes [##REF##12690453##15##] and the posterior Hox genes that were recently identified in these animals are neither related to the AbdB nor Post1/2 classes, which are specific for ecdysozoans and lophotrochozoans, respectively [##UREF##2##16##].</p> <p>Recently, the phylogenomic approach has provided the opportunity to sum up the phylogenetic signal from hundreds of genes and thereby to increase the resolution of the phylogenies [##REF##15861208##17##]. Two different phylogenomic studies involving different chætognath species and based on different samples of nuclear genes have assessed the phylogenetic position of chætognaths. They have both provided strong support for the inclusion of chætognaths within protostomes [##REF##15861208##17##, ####REF##16890510##18##, ##REF##16890509##19####16890509##19##]. Matus <italic>et al. </italic>[##REF##16890509##19##] suggested the branching of chætognaths at the base of lophotrochozoans on the basis of 72 nuclear genes described as valuable phylogenetic markers by Philippe <italic>et al</italic>. [##REF##15703236##20##]. Conversely, using a slightly larger taxonomic sampling and 78 ribosomal protein (RP) genes, Marlétaz <italic>et al. </italic>[##REF##16890510##18##] proposed that chætognaths are the sister group of all other protostomes. This last hypothesis has deep implications for the evolution of developmental patterns among bilaterians since it promotes the view that deuterostome-like developmental features such as enterocoely or a secondary mouth opening may be ancestral among bilaterians. Interestingly, recent insights into the structure of the nervous system of chætognaths suggest that these organisms have an intra-epidermal non-centralized nerve plexus, such as those observed in hemichordates or cnidarians [##REF##17511857##6##]. This is another example of a putative ancestral characteristic in this phylum. Then, both the phylogenetic position of chætognaths and their peculiar morphology and development indicate that these organisms are pivotal for the understanding of animal evolution.</p> <p>The expressed sequence tag (EST) approach provides an interesting opportunity to survey genomes and to perform comparisons between organisms. For instance, whole transcriptome comparisons based on ESTs initially suggested that the gene repertory shared by all metazoans is larger than expected [##REF##14680636##21##]. Moreover, in regard to the unexpected genetic complexity of cnidarians, the evolutionary extent of gene losses observed in nematodes and <italic>Drosophila </italic>remains to be defined [##REF##14680636##21##]. Through their original phylogenetic position, chætognaths offer the opportunity to check whether the ancestral protostome transcriptome has already undergone such gene losses or remains close to the ancestral bilaterian gene set conserved between vertebrates and cnidarians. Furthermore, the identification of a core set of metazoan conserved genes from a large range of organisms provides marker genes for phylogenomic analyses and signature genes as rare genomic changes, which could lead to a reevaluation of animal phylogeny [##REF##16919363##22##,##REF##11050348##23##].</p> <p>Here, we describe an overview of <italic>Spadella cephaloptera </italic>genomics through fine-scale mining of consistent transcriptomic data. Although the morphology of chætognaths has been extensively described, only a few molecular studies have focused on these strange organisms. The transcriptome of chætognaths reveals a strong similarity with that of other bilaterians. This comparative framework allowed detection of molecular signatures and stressed the usefulness of RPs as marker genes for phylogenomic reconstruction. Along with the structural RNAs, RPs are major components of the ribosome translation complex [##REF##12007402##24##]. They constitute a set of remarkably conserved genes among eukaryotes, which have not been significantly affected by lineage-specific duplication [##REF##15800040##25##]. We took advantage of their high levels of expression, which allowed the assembly of a large dataset with extensive taxon sampling using ESTs. We then investigated the origin of the polymorphisms observed within the EST collection in the light of genome duplication or cryptic speciation as alternative explanatory hypotheses. Lastly, we found evidence for trans-splicing mRNA maturation in chætognaths from this EST data. This original mRNA processing mechanism involves the addition of a spliced-leader sequence at the 5' extremity of transcripts. This mechanism has been discovered in several animal phyla by analyzing other EST collections [##REF##15797620##26##]. Interestingly, the occurrence of trans-splicing in chætognaths has deep implications for the evolutionary origin and functional significance of this mechanism.</p>
<title>Materials and methods</title> <title>cDNA library and EST sequencing</title> <p>The juvenile-staged cDNA library of <italic>S. cephaloptera</italic>, Bush 1851 was previously described in [##REF##16890510##18##]. Briefly, polyA+ RNAs were isolated from hatchling to one-day-old juveniles collected in the coastal area near Marseille (Le Brusc, La Ciotat). After reverse transcription and selection of longer transcripts by size fractionation, cDNAs were cloned into Lambda Triplex 2 vector (Clontech, Palo Alto, CA, USA). There were 12,324 clones sequenced from a 5' primer at Génoscope (CNS, France). A total of 390 sequences were discarded due to vector contamination or cloning artifacts, yielding a total 11,934 sequences.</p> <title>BAC library</title> <p>The BAC library was constructed by BioS&amp;T Inc. (Montreal, Québec, Canada) in the pIndigoBAC-5 vector (Epicentre, Madison, WI, USA) from adult <italic>S. cephalotera </italic>genomic DNA. Average insert size is 135 kb. The library was arrayed and selected clones were sequenced using a shotgun approach at Génoscope.</p> <title>Analysis of the EST collection</title> <p>The sequences were annotated through Blast searches [##REF##9254694##68##] against SwissProt and TrEMBL databases. The transcripts were searched for reliable coding sequences using ESTScan and subsequently sorted into coding and non-coding sequences [##REF##14534179##69##]. Transcripts from the same putative genes were grouped into clusters using CLOBB [##REF##12398795##70##]. The functional Gene Ontology classification was retrieved from SwissProt matches using Fatigo [##REF##14990455##71##]. Transcripts that bear splice-leader at their 5' extremity were searched with Blastn using the SL sequences obtained from a preliminary screen. Data parsing processes were conducted using Perl scripting language and the Bioperl resource [##REF##12368254##72##].</p> <title>Databases and searches for marker genes</title> <p>The set of databases compared with the chætognath transcriptome includes the transcriptomes of <italic>D. melanogaster </italic>and <italic>H. sapiens </italic>downloaded from Ensembl [##UREF##7##73##]. The partial transcriptome of <italic>Lumbricus rubellus </italic>as well as clade-level EST databases were obtained from the NCBI EST database using appropriate keywords [##UREF##8##74##]. The recently released genomes of several lophotrochozoans (<italic>Capitella capitata, Helobdella robusta, Lottia gigantea, Aplysia californica</italic>, <italic>Schmidtea mediterranea </italic>and <italic>Schistosoma mansoni</italic>) were accessed directly on the JGI web site [##UREF##9##75##] or retrieved from the NCBI trace archive [##UREF##10##76##]. The respective scores were visualized using Simitri [##REF##12584125##77##].</p> <title>Duplicated gene detection</title> <p>All the transcript sequences encoding similar protein coding sequences were clustered on the basis of a similar SwissProt hit or a reciprocal blast hit after conceptual translation by ESTScan [##REF##14534179##69##]. Gene clusters were analyzed using Clustalw alignment, phylogenetic analyses and visual checking for the inclusion of potential diverging sequences (Table S1 in Additional data file 2). The minimal EST number (N) per cluster to detect the occurrence of two duplicated paralogous sequences at a 5% probability is 6 as estimated by 1 - 2(1/2)<sup>N</sup>. Ka/Ks were computed using PAML [##REF##17483113##78##].</p> <title>Phylogenetic analyses</title> <p>The RP dataset was built using a composite database approach so as to select the least divergent sequences for each taxon [##REF##16890510##18##]. The species set involved in each taxon as well as the data on missing sequences are available in Table S1 in Additional data file 2. This dataset was analyzed using the CAT model as implemented in Phylobayes [##REF##15014145##79##]. Maximum-likelihood and Bayesian inferences were performed using PhyML [##REF##14530136##80##], Treefinder [##REF##15222900##81##] and MrBayes [##REF##12912839##82##], respectively, and both assumed a WAG+Γ4+I model. For bayesian analyses (MrBayes and Phylobayes), at least two chains were launched, their convergence was verified and the corresponding burn-in period discarded.</p> <title>PCR amplifications of duplicated markers</title> <p>The individuals of <italic>S. cephaloptera </italic>employed for polymorphism analysis were collected across the <italic>Posidonia </italic>seagrass, in the calanque of Sormiou, in the coastal area near Marseille, France. Copies of nuclear RP S36 and RP L40 and mitochondrial Cytb genes were amplified by PCR using allele-specific primers (Table S3 in Additional data file 2). Genomic DNA extractions were carried out on heads and bodies of <italic>S. cephaloptera </italic>individuals from the calanque of Sormiou using the Wizard SV Kit (Promega, Madison, WI, USA) and PCR reactions were carried out using GoTaq polymerase (Promega). PCR conditions were defined according to the melting temperature determined for each pair of primers. Fragments were subsequently gel-purified, cloned in pGEM-T easy (Promega) and sequenced. PCR artifacts were excluded by two approaches: first, amplifications were performed on head and body of the same individual, which allows checking of the congruence between the sequences; then, the accuracy of DNA polymerase was verified by amplifying and sequencing cloned DNA. These validations were confirmed by the recovery of the amplified sequences in the EST library and by the high levels of substitutions between sequences, which are too large to be interpreted as PCR artifacts.</p> <title>Population genetic analyses</title> <p>Pairwise sequence comparisons and molecular evolution analyses were performed with Mega2 software assuming the kimura 2 parameters model [##REF##11751241##83##]. The demographic history of the populations was evaluated by mismatch analyses using Tajima's D and Fu and Li's D neutrality test implemented in DNAsp software [##REF##14668244##84##].</p> <title>Accession numbers</title> <p>The set of sequences analyzed in this paper has been deposited in GenBank with the following accession numbers: <ext-link ext-link-type="gen" xlink:href="CU555081">CU555081</ext-link>-<ext-link ext-link-type="gen" xlink:href="CU563075">CU563075</ext-link> and <ext-link ext-link-type="gen" xlink:href="CU693608">CU693608</ext-link>-<ext-link ext-link-type="gen" xlink:href="CU693592">CU693592</ext-link> for the 11,934 ESTs, <ext-link ext-link-type="gen" xlink:href="CU672232">CU672232</ext-link> for the BAC 35YA21, and <ext-link ext-link-type="gen" xlink:href="EU529890">EU529890</ext-link>-<ext-link ext-link-type="gen" xlink:href="EU529968">EU529968</ext-link> for the sequences of RP L36 and L40 paralogs in nine individuals (Table ##TAB##1##2##).</p>
<title>Results and discussion</title> <title>Partial transcriptome of the chætognath <italic>S. cephaloptera</italic></title> <p>The sequencing of an EST collection of the juvenile-staged chætognath <italic>S. cephaloptera </italic>offered the opportunity to explore the transcriptome of this evolutionarily significant organism. The survey of sequence length and quality supported the accuracy of these data (Figure S1 in Additional data file 1). During these steps, we noticed that 16% of sequences match mitochondrial rRNA sequences (12S and 16S rRNAs, Figure ##FIG##0##1##) probably because the long polyadenine stretches of these rRNA molecules were isolated by the oligos-dT employed for mRNA isolation (see Materials and methods). We attempted to build clusters that gathered all transcripts from a unique gene so as to deal with a non-redundant partial transcriptome. However, the low complexity regions of some ESTs, which did not include an accurate open reading frame, hindered this process. Thus, ESTs were sorted into predicted coding and non-coding sequences using conceptual translation, and the coding transcripts were retained for comparative analyses. The overall content of the EST collection was evaluated using these steps (Figure ##FIG##0##1##). We noticed that up to 54% of the ESTs could be non-coding polyadenylated RNA, a striking figure that is, however, similar to that obtained for the human genome [##REF##17486121##27##]. The removal of non-coding sequences greatly improved clustering efficiency, yielding 1,447 clusters, of which 459 include more than one sequence (Figure S1 in Additional data file 1). A total of 694 of these clusters have significant matches within a protein database (TrEMBL, score &gt;50) and 250 have clear homologs in this database with an average of 72% identity (score &gt;150). Among the transcripts that match nuclear coding genes, the RP genes are largely represented compared to other genes similar to SwissProt entries (Figure ##FIG##0##1##).</p> <p>The average gene content of the library was checked regarding functional annotation as implemented in Gene Ontology [##REF##11997345##28##]. The <italic>S. cephaloptera </italic>library exhibited a broad diversity of functional classes with a majority of transcripts involved in metabolism or cellular activities and a non-negligible amount of transcripts involved in development (Figure S2 in Additional data file 1), which is consistent with the juvenile stage of the animals used. Hence, this EST collection contains representative, high quality sequences, providing suitable material for comparative analyses.</p> <title>Gene core conservation</title> <p>The set of non-redundant chætognath transcripts was compared with several databases using the Blast program. These databases first included sets of transcripts of representative species belonging to the most important clades of bilaterians: <italic>Drosophila melanogaster </italic>as an ecdysozoan, <italic>Lumbricus terrestris </italic>as a lophotrochozoan and <italic>Homo sapiens </italic>as a deuterostome. These comparisons were depicted through the plotting of respective similarity scores for all transcripts that have a significant match to at least one of these species (score &gt;150, Figure ##FIG##1##2##). This comparison demonstrated that a pool of 141 transcripts is strongly conserved between these distantly related species (Figure ##FIG##1##2a##). Conversely, 169 transcripts did not have significant matches in one or two of the species despite their strong similarity between chætognath and the remaining species. This lack of homologs is generally imputed to extensive gene loss [##REF##14680636##21##]. Therefore, further comparisons were performed to identify genes whose homology assignment and gene loss in a peculiar lineage were unambiguous. Interestingly, the number of transcripts that did not match to one or more databases decreased from 169 to 74 when the complete set of sequences available for each bilaterian clade was employed as the database, instead of only one representative species (Figure ##FIG##1##2b##). The lack of homologous matches in some species could then be explained by an increase in evolutionary rates, which could have weakened the sequence similarity signal [##REF##15543149##29##]. Additionally, the similarity level of matches increased when composite databases were employed (Figure ##FIG##1##2##), which supports the interest in this approach for phylogenomic reconstruction [##REF##16890510##18##].</p> <p>Two classes of genes provide reliable information for phylogeny inference (Figure ##FIG##1##2b##). Those that are highly shared between distantly related taxa constitute a set of conserved genes that are valuable markers for constructing phylogenomic datasets. In parallel, the genes that lack a homologous copy in one of the considered clades represent meaningful signature genes whose loss is attributable to a discrete event [##REF##11050348##23##].</p> <p>The candidates for signature genes are the genes inferred to be lost in one of the investigated clades (Figure ##FIG##1##2b##). Those candidates were carefully examined and their presence checked in the largest sets of available ESTs and full genome sequences. These data include the newly sequenced full genomes of lophotrochozoans and is assumed to include an exhaustive gene set in these species. Numerous candidate genes were invalidated because their homology relationships are disputable or because a homolog was retrieved from the full genome sequences surveyed. Among these candidates, the guanidinoacetate N-methyltransferase (GAMT) enzyme was recovered in the chætognath <italic>S. cephaloptera</italic>, in all studied deuterostomes, cnidarians and sister groups of metazoans (Figure S3 in Additional data file 1) but was not retrieved in any of the protostomes surveyed. Notably, this GAMT enzyme was also recovered in the acoel <italic>Convoluta pulchra</italic>, which was recently excluded from the protostomes [##REF##17684563##30##]. This enzyme catalyzes the key step of creatine synthesis, an activity that was previously checked biochemically in a variety of organisms but was not found in selected protostomes [##UREF##3##31##]. GAMT was later noticed as missing in <italic>D. melanogaster</italic>, <italic>Anopheles gambiae </italic>and <italic>Caeorhabditis elegans </italic>genomes [##REF##12364792##32##]. The presence of this ancient gene provides strong evidence for an early divergence of chætognaths from other protostomes. Indeed, the most parsimonious scenario states that this gene was lost in the protostome lineage after its split with chætognaths [##REF##16890510##18##].</p> <title>Selection of marker genes for metazoan phylogeny</title> <p>We attempted to evaluate the phylogenetic properties of the conserved genes that share equal levels of similarity with the main animal clades with respect to the convenience of their orthology assignment, their abundance in EST data and their molecular evolution properties. The main concerns when constructing phylogenomic-class datasets, especially from EST data, are the discarding of paralogous sequences, the removal of contaminants and the limitation of missing data. According to these criteria, we argue here that the set of RP genes is one of the best for setting up phylogenomic analysis in a large sample of taxa.</p> <p>Among the 694 chætognath genes similar to a database entry, only 267 genes have homologs in the three main clades of bilaterians (score &gt;150, Figure ##FIG##1##2b##). Copies of each selected marker were retrieved for all phyla studied for which EST data are available (Figure ##FIG##2##3##). In this way, the missing data were estimated through the occurrence of each gene in EST collections and preliminary phylogenetic analyses were carried out for all these independent alignments. Such controls unexpectedly highlighted putative paralogy problems for many candidate markers. If the orthologous transcript of a surveyed gene is missing in a non-exhaustive EST collection, a paralogous relative of this gene could be retrieved instead, with little chance of detection. Among candidate marker genes, RPs exhibit no ancient duplicates or out-paralogs and constitute a class of markers free from potential paralogy assignment problems [##REF##15800040##25##,##REF##12446146##33##]. Moreover, the gene-specific trees allowed detection of some contaminants in the EST collections, through the verification of unexpected clusterings in the tree (for example, several EST collections of parasitic organisms being contaminated by transcripts from their hosts).</p> <p>Next, the amount of missing data was estimated using these raw alignments and compared with the number of ESTs in each available collection (Figure ##FIG##2##3##). The positive correlation observed between the number of ESTs and the completeness of the dataset is stronger when dealing with a dataset composed of RPs. For instance, the 5,235 EST collection of tardigrades yielded a dataset that is 77% complete for RPs, but only 35% complete for non-ribosomal markers. Thus, their large representation in EST collections strengthens the usefulness of RPs as phylogenetic markers.</p> <title>Chætognaths within renewed metazoan phylogeny</title> <p>In order to assess the branching of chætognaths and to stress the usefulness of RP genes for phylogenomics, a RP dataset was assembled using the composite dataset approach [##REF##16890510##18##]. This method depends on the selection of the least diverging copy of each marker gene in each taxon, such as a phylum, and thus allows reduction of the branch lengths of composite taxa (Table S1 in Additional data file 2). To overcome previous problems, both taxon sampling and inference methods were improved. Several new phyla were included in this analysis and, in particular, numerous protostome groups: priapulids, platyhelminthes, nermerteans, ectoprocts, entoprocts and rotifers [##REF##17051155##34##, ####REF##17921486##35##, ##REF##18222945##36####18222945##36##]. Most rotifer sequences were retrieved from <italic>Oryza sativa </italic>(rice) ESTs, where they exist as contaminants, using their very specific splice-leader sequence as an anchor (see below and [##REF##15788744##37##]). Rotifers constitute a key phylum with respect to chætognaths because they were sometimes grouped together in the gnathifera clade on the basis of morphological criteria [##UREF##4##38##]. Alternatively, a splitting of lophotrochozoans into two main lineages, the platyzoans (uniting platyhelminthes and rotifers) and the trochozoans (mainly annelids, molluscs, lophophorates and nermertes) has been proposed [##REF##12116426##39##,##UREF##5##40##]. Otherwise, in addition to the traditional site-homogenous WAG model, we have assessed the phylogeny of bilaterians using the site-heterogeneous CAT model, which recently improved the limitation of the long-branch attraction artifact, a common pitfall in phylogenetic reconstruction [##REF##17012374##41##,##REF##17288577##42##]. The inclusion of the most recently released EST data for this large set of phyla led to a dataset including 11,730 amino acid positions and 25 taxa (Additional data file 4).</p> <p>The analysis of this dataset confirmed the branching of chætognaths at the base of the protostomes with significant support values for both the site-homogeneous WAG model and the site-heterogeneous CAT model (bootstrap proportion (BP) of 76 and posterior probability (PP) of 1; Figure ##FIG##3##4a,b##). The inclusion of chætognaths within protostomes is still firmly supported (BP 95, PP 1; Figure ##FIG##3##4##). The inclusion of new taxa strengthens support for both the ecdyozoa and lophotrochozoa clades but the exact relationships within these two clades remain elusive [##REF##17921486##35##,##REF##18222945##36##,##REF##18192187##43##]. Chætognaths and rotifers do not exhibit any peculiar affinities, prompting us to reject the gnathifera hypothesis [##UREF##4##38##]. Conversely, the branching of rotifers is problematic since this phylum is alternatively included in ecdysozoans and lophotrochozoans depending on the use of, respectively, site heterogeneous or homogeneous models (Figure ##FIG##3##4##). Thus, the clustering of platyhelminthes and rotifers in a platyzoa clade is supported by the WAG model but rejected by the CAT model, suggesting that this grouping may be somehow related to long-branch attraction (Figure ##FIG##3##4##). Alternatively, previous studies based on morphology and SSU genes have not argued for the ecdysozoan affinities of rotifers [##UREF##4##38##,##REF##12116426##39##]. Surprisingly, CAT model analysis no longer succeeds in recovering the monophyly of the deuterostomes (Figure ##FIG##3##4b##). Instead, it provides limited support for the successive divergence of chordates and ambulacrarians (echinoderms and hemichordates; PP 0.9; Figure ##FIG##3##4b##). This striking topology was recovered by an independent study using the same heterogeneous CAT model [##REF##18192187##43##] but was neither confirmed by WAG analyses (BP 89 for the monophyly of deuterostomes; Figure ##FIG##3##4a##) nor supported on morphological bases [##REF##17051155##34##,##UREF##4##38##]. One can consider that the two unexpected branchings of rotifers and deuterostomes may be related to some artifact affecting the CAT model, such as sensitivity toward compositional biases [##REF##18234708##44##]. Finally, the placozoan <italic>Trichoplax adherens </italic>surprisingly clustered within the poriferans, as a sister group of the homoscleromorphs (BP 91, PP 0.94; Figure ##FIG##3##4##), although this poriferan status has never been suggested before [##REF##16731622##45##,##REF##16299758##46##]. These challenging hypotheses will be investigated in further studies because they have deep implications for the evolution of metazoans (F Marlétaz <italic>et al</italic>., in progress).</p> <p>Through extended taxon sampling and improved substitution models, these analyses strongly confirm our previous statements about basal-protostome branching of chætognaths [##REF##16890510##18##] and exclude the basal-lophotrochozan hypothesis [##REF##16890509##19##]. Although some areas of bilaterian trees are sometimes incongruent depending on models and inference methods, the position of chætognaths remains remarkably stable throughout our analyses. Furthermore, this branching is not only supported by the presence of GAMT, an unambiguous molecular signature, but also by the posterior Hox genes of chætognaths that are not related to the classes specific to ecdysozoans (Abd-B) or lophotrochozoans (Post1/2) [##UREF##2##16##]. Finally, this topology was also recovered by independent studies involving alternative gene and taxon sampling [##REF##17684563##30##,##REF##17921486##35##,##REF##18192187##43##]. In a broader perspective, the strengthening of their phylogenetic position makes chætognaths a key model for comparative genomics among bilaterians.</p> <title>Genome duplication in the chætognath phylum</title> <p>The clustering of similar sequences indicated that alternative nucleotide forms are present among the transcripts encoding the same protein. Two distinct forms are observed in most cases, although three forms encode some proteins. These forms are separated by a large amount of molecular divergence and can also be distinguished by their different 5' and 3' untranslated regions (UTRs), suggesting that they correspond to different genes (Figure ##FIG##4##5## and Additional data files 5 and 6).</p> <p>Ka/Ks ratios were calculated for all pairs of diverging forms to consider the impact of the nucleotide divergence on the protein sequences. The values of Ka/Ks range from 0.001-0.154 with a median value of 0.004, which confirmed the strong conservation of amino acid sequence despite the large synonymous substitutions observed in some cases (Ks values range from 0.8-75; Table S2 in Additional data file 2). These distinct forms were mainly retrieved for the most highly expressed genes, among which RP genes are prominent (Table ##TAB##0##1##). We verified that the observed molecular divergence could not be explained by the clustering of distant paralogous sequences. For the genes that have clear homologs among metazoans, the sequences of alternative forms always cluster together in phylogenetic analyses and are thus strongly separated from homologous genes of other animals. For instance, the RUX genes have undergone an ancient duplication resulting in the RUX-E and RUX-G paralogs in all metazoans. Interestingly, chætognaths display up to three forms of RUX-E and two forms of RUX-G, all these forms being closely related (Figure S4 in Additional data file 1; Additional data file 7).</p> <p>Such a pattern could be explained by either the duplication of a large set of genes in the genome of chætognaths or, alternatively, it could be explained by the presence of cryptic species within the sampled population. In the first case, the observed differences would be attributed to the divergence between paralogous genes originating through the duplication, where the genome of one individual is expected to contain the two alternative nucleotide forms. In the second case, the observed genetic differences would be caused by the genetic divergence between the orthologous genes of several cryptic species spread among the population, where one individual is thus expected to contain only one of the alternative forms.</p> <p>This cryptic speciation hypothesis may be supported by the strong polymorphism also observed for all genes of the mitochondrial genome, which constitutes an independent lineage from the nuclear genome. For example, cytochrome b (Cytb) transcripts but also cytochrome oxydase I and III are split into distinct forms separated by large molecular distances (Figure ##FIG##4##5c##; Figure S4 in Additional data file 1; Additional data files 8-10), thus testifying to the presence of distinct mitochondrial lineages within the sampled population.</p> <p>To decide between these hypotheses, we designed a PCR screen to survey the alternative forms of selected markers in independent individuals. The genes for RPs L36 and L40 were targeted because they are nuclear genes displaying two alternative forms with the highest number of transcripts in the library (Table S2 in Additional data file 2). The mitochondrial Cytb gene served as an independent reference for the interpretation of results from nuclear genes. The three distinct forms of this strongly diverging mitochondrial gene were surveyed in all the individuals tested (Figure ##FIG##4##5c##). Chætognaths are hermaphroditic and, after fertilization, they store exogenous sperm in their sperm receptacles (Figure ##FIG##4##5a##), which makes it possible to amplify the DNA from another individual. Hence, in order to detect such contamination, we performed independent amplifications on heads, which are considered free from sperm contamination, as well as on the rest of the body, which contained sperm receptacles (Figure ##FIG##4##5a##). The experimental design made it possible to detect alternative forms through the amplification of specific DNA fragments of distinct sizes (Figure ##FIG##4##5b##; Table S3 in Additional data file 2). The PCR products were characterized by sequencing and nucleotide polymorphism was subsequently carefully examined. In addition to their nucleotide divergence in coding sequences, the distinct forms of nuclear genes for RPs L36 and L40 have alternative intron positions and lengths as well as differences in their 5' and 3' UTR regions (Figure ##FIG##4##5b##).</p> <p>Performed on nine individuals, the amplifications revealed the presence of the two forms of the nuclear genes for RPs L36 and L40 in each individual (Table ##TAB##1##2##). Conversely, only one form of the mitochondrial Cytb gene was amplified in each individual with the exception of the body of individual 1, which includes two forms, thus suggesting contamination by exogenous sperm (Table ##TAB##1##2##). The amplification of the divergent nucleotide forms within one individual indicates that the alternative nucleotide forms correspond to paralogous nuclear copies originating through past gene duplication events (Table ##TAB##0##1##). Conversely, the alternative forms of the mitochondrial gene correspond to variation within the population. Because some genes, such as that encoding Translationally controlled tumor protein (TCTP), do not present paralogous copies despite their high expression levels (112 TCTP transcripts in the EST collection; Table S2 in Additional data file 2), we addressed the extent of these duplications in evaluating the quantity of duplicated genes. If the clusters of transcripts encoding the same protein include all the transcripts from alternative paralogous genes and if those paralogous genes have similar levels of expression, the probability that transcripts from these paralogous genes are represented in a given cluster is related to the size of this cluster (see Materials and methods). Hence, all the clusters that include more than six transcripts have at least a 95% chance of including transcripts from the two copies if they exist. Such clusters of transcripts were all checked for paralogous copies through sequence alignments and trees. Paralogs were detected within 35 of the 66 clusters investigated, which suggests that up to 69% of chætognath genes are the products of duplications. These paralogs could have arisen through either a whole genome duplication (WGD) event followed by an extensive gene loss, or several segmental duplication events.</p> <p>The hypothesis of a WGD event is reinforced by the high occurrence of RPs among duplicated genes (Table ##TAB##0##1##). The trend to retain RP genes was previously observed after WGD for <italic>Paramecium tetraurelia</italic>, yeast and plants [##REF##17086204##47##, ####REF##15208398##48##, ##REF##10508730##49####10508730##49##] but is not a common occurrence in small-scale duplications. Conversely, it is difficult to understand why the paralogous genes have been retained after their duplication and maintained under purifying selection as emphasized by Ka/Ks values. This conclusion is in contradiction with the current view of gene destiny after genome duplication, which alternatively predicts that one of the gene duplicates is lost or undergoes the accumulation of substitutions [##REF##10101175##50##]. Using a genome-level dataset, similar findings were made about the strongly duplicated genome of <italic>Paramecium </italic>where the retention of duplicated genes was accounted in part by dosage compensation constraints [##REF##17086204##47##].</p> <p>The most plausible dating is that this duplication occurred before the diversification of the major chætognath lineages. Two copies of SSU and LSU were retrieved in members of the phylum dispersed all over the tree of chætognaths [##REF##16434216##10##, ####REF##8336548##11##, ##REF##9069174##12####9069174##12##]. Moreover, the survey of 226 ESTs available for <italic>Flaccisagitta enflata </italic>also revealed the presence of alternative nucleotide forms for some genes (data not shown), which would confirm that the duplication is not limited to SSU/LSU genes at this taxonomic scale. Further genome data would be required to date the duplication, for instance, in considering the Ks distribution of the set of paralogs [##REF##14757817##51##], and also to definitively state the nature of the duplication through the analysis of synteny in duplicated blocks of the genome. Nevertheless, this preliminary transcriptomic survey stresses the usefulness of the chætognaths to study phylum-level genome duplication events and the destiny of paralogous genes.</p> <title>Population genomics</title> <p>Beyond the molecular divergence between the coding sequences of duplicated paralogous genes, a subsequent survey of the genomic sequences of selected genes revealed that the level of polymorphism is strong within each paralogous gene (Table S4 in Additional data file 2). Multiple nucleotide substitutions as well as insertion/deletion events (indels) occurred within the introns of the four selected nuclear genes (paralogous copies of the genes for both RPs L36 and L40; Additional data files 11-14). Similarly, a large number of substitutions have accumulated in the various mitochondrial genes, thus revealing distinct mitochondrial lineages within the sampled population (Figure ##FIG##4##5c##; Figure S4 in Additional data file 1). However, these strong levels of divergence remain consistent with a population genetic structure because of the regular AT composition and the limited degree of saturation revealed by Ts/Tv ratios, singleton positions being essentially transition substitutions (Table S4 in Additional data file 2; Figure S6 in Additional data file 1).</p> <p>We attempted to determine the origin of this population genetic heterogeneity, which could, for instance, be due to a cryptic speciation or to a past hybridization. For this, the sequences of each individual were compared using phylogenetic trees and indels as discrete informative characteristics (Figure ##FIG##5##6##). For each marker gene, individual sequences split into several major clades supported by strong bootstrap and discrete indel events, which allows unambiguous identification of heterozygous individuals (Figure ##FIG##5##6##). For example, individual 4 is heterozygous for all markers and individuals 6, 9 and 3 are heterozygous for at least one marker. Moreover, the occurrence of several cases of putative recombinations between alleles highlights the heterozygous status of some individuals (individuals 3 and 4, Figure ##FIG##5##6b,d##). Notably, our PCR-based experimental design provided positive evidence only for heterozygosis because two amplifications (head and body) were carried out per individual, yielding 0.5 probability to detect heterozygosity. Heterozygous individuals could thus be even more abundant than observed. These heterozygous cases convincingly demonstrate that a shuffling occurs between the most divergent alleles of each gene, which constitutes strong evidence for interbreeding within the sampled population. This finding definitely excludes the possibility of cryptic speciation within this <italic>S. cephaloptera </italic>population. Alternatively, the panmixy hypothesis was confirmed by the unimodal distribution of pairwise divergences in mismatch analysis, which is consistent with constant population size and excludes a past hybridization event (Figure S6 in Additional data file 1). Finally, the distinct mitochondrial lineages are spread within the population but they are not correlated with any haplotype differentiation at the nuclear level, which is a strong argument against the cryptic speciation hypothesis. This type of mitochondrial diversity was previously discovered for the planktonic species <italic>Sagitta setosa </italic>but was also interpreted with difficulty [##REF##16968273##52##].</p> <p>Strikingly, these comparisons also highlighted molecular divergence between the head and the body of some individuals for each of the five markers investigated (Figure ##FIG##5##6## and Additional data file 4). Such substitutions cannot be explained by a heterozygous status of those individuals because sequences from head and body were firmly clustered in the tree (Figure ##FIG##5##6##). For example, individual 4 exhibits well-separated alleles present in both head and body but intra-individual substitution took place between head and body for both of these alleles (Figure ##FIG##5##6c##). This pattern of substitutions may be explained by the occurrence of somatic mutations during the life of individuals. This interpretation is corroborated by the large extent of intra-individual substitutions in all marker genes and all individuals. Somatic mutations are considered as rare conditions, mainly known from related disorders in humans [##REF##12360233##53##]. Less clear are the evolutionary implications and putative benefits of this phenomenon [##REF##15525396##54##]. They are sometimes suspected to play a prominent role in apoptosis and possibly in the regulation of cell division [##REF##15525396##54##]. Moreover, somatic mutations have been demonstrated to be more widespread in <italic>Drosophila </italic>than in mammals [##REF##17435764##55##], and are sometimes correlated with extensive chromosome rearrangement in the <italic>Drosophila </italic>lineage [##REF##17550304##56##]. However, little is known about the extent and importance of this process in the non-model organisms. In the case of the chætognath, somatic mutation could be due to the high mutation rates that seem to affect both germline and soma and could explain the divergence at the population and individual levels. The possible relationship of these accelerated mutation rates with structural reshaping of the genome after duplication deserves further evaluation.</p> <p>Notably, this level of somatic mutation generates a strong background noise that hinders the accurate interpretation of point mutations related to the diversity of haplotypes. Moreover, traditional hypotheses of population genetics are challenged by our findings: the genetic distances observed between individuals of a single population reach species-level without any evidence for cryptic speciation or past hybridization. In parallel, multiple mitochondrial lineages diverge and are spread and maintained within a single population [##REF##16968273##52##]. If such features are revealed as more widespread than expected, the data collected from numerous population genetic or phylogeographic studies require more cautious interpretation. This observation pleads for an increase of the sampling depth for population genetics, especially through genomic approaches [##REF##14631358##57##].</p> <title>Trans-splicing transcript maturation in chætognaths</title> <p>The survey of the chætognath EST collection detected a common 36 nucleotide motif shared at their 5' end between transcripts from unrelated genes. Selected cases showed that these short motifs are absent from upstream genomic regions of the genes (Figure ##FIG##6##7c##; Additional data files 15-17). These findings strongly suggest that mRNA of <italic>S. cephaloptera </italic>undergoes trans-splicing maturation. This mRNA processing occurs through the spliceosomal transfer of a small RNA molecule at the 5' end of the mRNA from an independent spliced-leader (SL) gene and has been described in numerous species [##REF##15797620##26##]. To exclude the possibility of an artifact, the 226 available EST sequences of <italic>F. enflata</italic>, another chætognath species, were screened using the identified SL sequences from <italic>S. cephaloptera</italic>. This survey led to the recovery of similar SL sequences at the 5' end of several cDNA sequences. As <italic>F. enflata </italic>and <italic>S. cephaloptera </italic>belong to the two main orders of chætognaths [##REF##16434216##10##], the trans-splicing mechanism is likely to be present in the whole chætognath phylum.</p> <p>An extensive study of trans-splicing in the nematode phylum has previously revealed the presence of different forms of SL sequences [##REF##17121468##58##]. Similarly, several different SL sequences were retrieved from <italic>S. cephalotera </italic>(Table ##TAB##2##3##). This number is consistent with the strong level of polymorphism previously observed for coding sequences but two distinct SL forms alone (SL1 and SL2) represent 87% of the SL sequences (Table ##TAB##2##3##). These forms do not exhibit any specificity for the different paralogs since SLs are added randomly to transcripts from these paralogs (Figure ##FIG##6##7##; Additional data files 5 and 6). In an attempt to understand the evolutionary history of trans-splicing within the chætognath phylum, the SL forms of <italic>F. enflata </italic>were compared with those of <italic>S. cephaloptera</italic>. This set of chætognath SL sequences splits into four main forms: two of them, SL1 and SL3, are present in both the species investigated whereas the two other forms, SL2 and SL4, are specific for <italic>S. cephaloptera </italic>and <italic>F. enflata</italic>, respectively (Table ##TAB##2##3##). Otherwise, neither of these chætognath SL forms display similarity with the SL of another phyla. This finding suggests that just as in nematodes, the evolution of alternative forms could have occurred at a relatively reduced taxonomic scale [##REF##17121468##58##].</p> <p>Within the EST collection, 2,914 sequences exhibit SL addition, which represents 30% of nuclear transcripts (Figure ##FIG##7##8a##). Among the SL population, 72% are coding transcripts, of which 46% have a homolog in SwissProt. The clustering of similar coding trans-spliced transcripts indicated that 41% of putative genes undergo SL addition. Furthermore, the relationship between trans-splicing and expression level was tested through the comparison of the number of ESTs per cluster of trans-spliced or non-trans-spliced transcripts. If we posit that this number can be considered as an estimate of the expression level, trans-spliced genes are significantly more expressed than others (Wilcoxon rank test, <italic>p </italic>&lt; 2.2e<sup>-16</sup>). For instance, among the 50 more expressed genes (that is, biggest EST clusters), only two are not trans-spliced. These values suggest that trans-splicing is involved in the regulation of a set of strongly expressed genes responsible for key cellular functions, for example, the RP set.</p> <p>Hence, the set of trans-spliced genes of <italic>S. cephaloptera </italic>was compared with those of other animals that carry out this RNA maturation process (Table ##TAB##3##4##). Particularly, trans-splicing was characterized in the nematode <italic>C. elegans </italic>model, for which a large EST set is available [##REF##17121468##58##,##REF##12560808##59##]. These data were especially useful for extensive comparisons with <italic>S. cephaloptera </italic>(Figure ##FIG##7##8b##). The amount of ESTs with a SL is slightly smaller in <italic>C. elegans </italic>(22%) than in <italic>S. cephaloptera </italic>(29%) and trans-spliced transcripts are missing from the very reduced non-coding ESTs, only 5% of the total EST collection in <italic>C. elegans</italic>. This comparison suggests that trans-splicing is at least as widespread in <italic>S. cephaloptera </italic>as in <italic>C. elegans</italic>. Furthermore, the set of genes affected by trans-splicing in <italic>S. cephaloptera </italic>was comparable with those of <italic>C. elegans</italic>. Among the 119 different genes that are both trans-spliced and annotated using SwissProt with high confidence, 79 have trans-spliced homologous genes in <italic>C. elegans </italic>(Additional data file 3). These values do not include the 78 RP set, which are conserved and trans-spliced in both species. While the molecular actors involved in trans-splicing do not exhibit evolutionary conservation [##REF##15797620##26##], the genes undergoing this kind of processing could otherwise be similar in distantly related species.</p> <p>Trans-splicing has been discovered in a set of organisms spread all over the phylogenetic tree of eukaryotes, but the molecular actors involved in this RNA maturation process do not exhibit any evolutionary conservation: the SL sequences cannot be aligned and the ribonucleoproteic machineries involved in this process have no homologs in other species [##REF##15797620##26##,##REF##15788744##37##,##REF##9297698##60##, ####REF##16822859##61##, ##REF##15972844##62####15972844##62##]. However, the discovery of SL trans-splicing in new animal phyla strongly argues for an ancient origin of this process, especially if these phyla have significant phylogenetic positions, such as chætognaths or acoels. Noticeably, we also found evidence for SL addition in recently released ESTs from the acoels (unpublished observations), worm-like animals that were recently excluded from the protostomes [##REF##17684563##30##]. Moreover, an extensive comparison of available genomic data revealed that the set of genes undergoing trans-splicing is evolutionarily conserved between those species, which is another indication of a putative common evolutionary origin of trans-splicing (Additional data file 3 and Table ##TAB##3##4##).</p> <title>Operonic transcription in <italic>S. cephaloptera</italic></title> <p>In the nematode chætognath <italic>C. elegans </italic>and the tunicate <italic>Ciona intestinalis</italic>, polycistronic mRNA molecules that contain two or more genes are transcribed from operon structures and subsequently resolved by trans-splicing [##REF##17121468##58##,##REF##12560808##59##]. We mined for similar eukaryotic operons in our pool of sequenced bacterial artificial chromosomes (BACs) in looking for clusters of genes with the same transcription orientation and very short intergenic regions. The trans-splicing maturation of these putative operons was confirmed by mapping the tran-spliced ESTs onto these genomic sequences. One operon including RP S14 and PCNA genes was successfully identified in the BAC 35A21 that otherwise contains a large number of genes undergoing trans-splicing (Figure ##FIG##6##7a##; Table S5 in Additional data file 2; Additional data files 15-17). EST mapping and <italic>ab initio </italic>gene prediction suggested that only three bases separate the polyadenylation site of the upstream RP S14 gene from the acceptor splice site at the PCNA downstream gene. This distance is even shorter than that in <italic>C. intestinalis </italic>operons [##REF##16822859##61##] and clearly excludes the possibility for a re-initiation of transcription between the genes (Figure ##FIG##6##7c##). This BAC sequence contains several other predicted genes that are closely clustered, share the same orientation and could thus belong to operonic structures (Figure ##FIG##6##7##; Table S5 in Additional data file 2). The transcripts of gene RP S14 include the two major SL 1 or 2 forms (Figure ##FIG##6##7c##). Similarly, transcripts from both paralogous copies of duplicated genes bear any of the alternative SL forms. Hence, contrary to <italic>C. elegans </italic>for which a SL form is specific for the upstream and the downstream genes in operons [##REF##17121468##58##], no specificity of the SL for gene position in the operon was recovered in <italic>S. cephaloptera </italic>(Figure ##FIG##6##7##).</p> <p>Hence, we have provided arguments for polycistronic transcription in another animal phylum, the chætognaths. Animal operons have previously been described in the nematodes and in the tunicate <italic>C. intestinalis </italic>and were often interpreted as an adaptation to the genome size reduction observed in both these organisms [##REF##15797620##26##]. However, the genome size of <italic>S. cephaloptera </italic>is rather large at 1.05 Gb [##UREF##6##63##] and the transcriptome of chætognaths is conserved among bilaterians, in terms of gene similarity and content. Furthermore, extensive study of genes co-transcribed in <italic>C. elegans </italic>led to the conclusion that the operons constitute functional units in grouping genes involved in similar pathways [##REF##12560808##59##]. Hence, the presence of operonic transcription in the chætognath, a basal protostome, strongly suggests that this mechanism should be more likely considered as an original mode of gene regulation, whose evolutionary importance may have been underestimated. Thanks to the pivotal comparative genomic abilities of chætognaths, the comparison of operon structure and content between nematodes and chætognaths should allow the extension of the functional results obtained in <italic>C. elegans </italic>to a larger set of organisms, including vertebrates. Such an approach would require further genomic data in the chætognath but will likely lead to the characterization of new gene regulatory networks, potentially leading to a better understanding of some genetic disorders [##REF##12560808##59##].</p>
<title>Results and discussion</title> <title>Partial transcriptome of the chætognath <italic>S. cephaloptera</italic></title> <p>The sequencing of an EST collection of the juvenile-staged chætognath <italic>S. cephaloptera </italic>offered the opportunity to explore the transcriptome of this evolutionarily significant organism. The survey of sequence length and quality supported the accuracy of these data (Figure S1 in Additional data file 1). During these steps, we noticed that 16% of sequences match mitochondrial rRNA sequences (12S and 16S rRNAs, Figure ##FIG##0##1##) probably because the long polyadenine stretches of these rRNA molecules were isolated by the oligos-dT employed for mRNA isolation (see Materials and methods). We attempted to build clusters that gathered all transcripts from a unique gene so as to deal with a non-redundant partial transcriptome. However, the low complexity regions of some ESTs, which did not include an accurate open reading frame, hindered this process. Thus, ESTs were sorted into predicted coding and non-coding sequences using conceptual translation, and the coding transcripts were retained for comparative analyses. The overall content of the EST collection was evaluated using these steps (Figure ##FIG##0##1##). We noticed that up to 54% of the ESTs could be non-coding polyadenylated RNA, a striking figure that is, however, similar to that obtained for the human genome [##REF##17486121##27##]. The removal of non-coding sequences greatly improved clustering efficiency, yielding 1,447 clusters, of which 459 include more than one sequence (Figure S1 in Additional data file 1). A total of 694 of these clusters have significant matches within a protein database (TrEMBL, score &gt;50) and 250 have clear homologs in this database with an average of 72% identity (score &gt;150). Among the transcripts that match nuclear coding genes, the RP genes are largely represented compared to other genes similar to SwissProt entries (Figure ##FIG##0##1##).</p> <p>The average gene content of the library was checked regarding functional annotation as implemented in Gene Ontology [##REF##11997345##28##]. The <italic>S. cephaloptera </italic>library exhibited a broad diversity of functional classes with a majority of transcripts involved in metabolism or cellular activities and a non-negligible amount of transcripts involved in development (Figure S2 in Additional data file 1), which is consistent with the juvenile stage of the animals used. Hence, this EST collection contains representative, high quality sequences, providing suitable material for comparative analyses.</p> <title>Gene core conservation</title> <p>The set of non-redundant chætognath transcripts was compared with several databases using the Blast program. These databases first included sets of transcripts of representative species belonging to the most important clades of bilaterians: <italic>Drosophila melanogaster </italic>as an ecdysozoan, <italic>Lumbricus terrestris </italic>as a lophotrochozoan and <italic>Homo sapiens </italic>as a deuterostome. These comparisons were depicted through the plotting of respective similarity scores for all transcripts that have a significant match to at least one of these species (score &gt;150, Figure ##FIG##1##2##). This comparison demonstrated that a pool of 141 transcripts is strongly conserved between these distantly related species (Figure ##FIG##1##2a##). Conversely, 169 transcripts did not have significant matches in one or two of the species despite their strong similarity between chætognath and the remaining species. This lack of homologs is generally imputed to extensive gene loss [##REF##14680636##21##]. Therefore, further comparisons were performed to identify genes whose homology assignment and gene loss in a peculiar lineage were unambiguous. Interestingly, the number of transcripts that did not match to one or more databases decreased from 169 to 74 when the complete set of sequences available for each bilaterian clade was employed as the database, instead of only one representative species (Figure ##FIG##1##2b##). The lack of homologous matches in some species could then be explained by an increase in evolutionary rates, which could have weakened the sequence similarity signal [##REF##15543149##29##]. Additionally, the similarity level of matches increased when composite databases were employed (Figure ##FIG##1##2##), which supports the interest in this approach for phylogenomic reconstruction [##REF##16890510##18##].</p> <p>Two classes of genes provide reliable information for phylogeny inference (Figure ##FIG##1##2b##). Those that are highly shared between distantly related taxa constitute a set of conserved genes that are valuable markers for constructing phylogenomic datasets. In parallel, the genes that lack a homologous copy in one of the considered clades represent meaningful signature genes whose loss is attributable to a discrete event [##REF##11050348##23##].</p> <p>The candidates for signature genes are the genes inferred to be lost in one of the investigated clades (Figure ##FIG##1##2b##). Those candidates were carefully examined and their presence checked in the largest sets of available ESTs and full genome sequences. These data include the newly sequenced full genomes of lophotrochozoans and is assumed to include an exhaustive gene set in these species. Numerous candidate genes were invalidated because their homology relationships are disputable or because a homolog was retrieved from the full genome sequences surveyed. Among these candidates, the guanidinoacetate N-methyltransferase (GAMT) enzyme was recovered in the chætognath <italic>S. cephaloptera</italic>, in all studied deuterostomes, cnidarians and sister groups of metazoans (Figure S3 in Additional data file 1) but was not retrieved in any of the protostomes surveyed. Notably, this GAMT enzyme was also recovered in the acoel <italic>Convoluta pulchra</italic>, which was recently excluded from the protostomes [##REF##17684563##30##]. This enzyme catalyzes the key step of creatine synthesis, an activity that was previously checked biochemically in a variety of organisms but was not found in selected protostomes [##UREF##3##31##]. GAMT was later noticed as missing in <italic>D. melanogaster</italic>, <italic>Anopheles gambiae </italic>and <italic>Caeorhabditis elegans </italic>genomes [##REF##12364792##32##]. The presence of this ancient gene provides strong evidence for an early divergence of chætognaths from other protostomes. Indeed, the most parsimonious scenario states that this gene was lost in the protostome lineage after its split with chætognaths [##REF##16890510##18##].</p> <title>Selection of marker genes for metazoan phylogeny</title> <p>We attempted to evaluate the phylogenetic properties of the conserved genes that share equal levels of similarity with the main animal clades with respect to the convenience of their orthology assignment, their abundance in EST data and their molecular evolution properties. The main concerns when constructing phylogenomic-class datasets, especially from EST data, are the discarding of paralogous sequences, the removal of contaminants and the limitation of missing data. According to these criteria, we argue here that the set of RP genes is one of the best for setting up phylogenomic analysis in a large sample of taxa.</p> <p>Among the 694 chætognath genes similar to a database entry, only 267 genes have homologs in the three main clades of bilaterians (score &gt;150, Figure ##FIG##1##2b##). Copies of each selected marker were retrieved for all phyla studied for which EST data are available (Figure ##FIG##2##3##). In this way, the missing data were estimated through the occurrence of each gene in EST collections and preliminary phylogenetic analyses were carried out for all these independent alignments. Such controls unexpectedly highlighted putative paralogy problems for many candidate markers. If the orthologous transcript of a surveyed gene is missing in a non-exhaustive EST collection, a paralogous relative of this gene could be retrieved instead, with little chance of detection. Among candidate marker genes, RPs exhibit no ancient duplicates or out-paralogs and constitute a class of markers free from potential paralogy assignment problems [##REF##15800040##25##,##REF##12446146##33##]. Moreover, the gene-specific trees allowed detection of some contaminants in the EST collections, through the verification of unexpected clusterings in the tree (for example, several EST collections of parasitic organisms being contaminated by transcripts from their hosts).</p> <p>Next, the amount of missing data was estimated using these raw alignments and compared with the number of ESTs in each available collection (Figure ##FIG##2##3##). The positive correlation observed between the number of ESTs and the completeness of the dataset is stronger when dealing with a dataset composed of RPs. For instance, the 5,235 EST collection of tardigrades yielded a dataset that is 77% complete for RPs, but only 35% complete for non-ribosomal markers. Thus, their large representation in EST collections strengthens the usefulness of RPs as phylogenetic markers.</p> <title>Chætognaths within renewed metazoan phylogeny</title> <p>In order to assess the branching of chætognaths and to stress the usefulness of RP genes for phylogenomics, a RP dataset was assembled using the composite dataset approach [##REF##16890510##18##]. This method depends on the selection of the least diverging copy of each marker gene in each taxon, such as a phylum, and thus allows reduction of the branch lengths of composite taxa (Table S1 in Additional data file 2). To overcome previous problems, both taxon sampling and inference methods were improved. Several new phyla were included in this analysis and, in particular, numerous protostome groups: priapulids, platyhelminthes, nermerteans, ectoprocts, entoprocts and rotifers [##REF##17051155##34##, ####REF##17921486##35##, ##REF##18222945##36####18222945##36##]. Most rotifer sequences were retrieved from <italic>Oryza sativa </italic>(rice) ESTs, where they exist as contaminants, using their very specific splice-leader sequence as an anchor (see below and [##REF##15788744##37##]). Rotifers constitute a key phylum with respect to chætognaths because they were sometimes grouped together in the gnathifera clade on the basis of morphological criteria [##UREF##4##38##]. Alternatively, a splitting of lophotrochozoans into two main lineages, the platyzoans (uniting platyhelminthes and rotifers) and the trochozoans (mainly annelids, molluscs, lophophorates and nermertes) has been proposed [##REF##12116426##39##,##UREF##5##40##]. Otherwise, in addition to the traditional site-homogenous WAG model, we have assessed the phylogeny of bilaterians using the site-heterogeneous CAT model, which recently improved the limitation of the long-branch attraction artifact, a common pitfall in phylogenetic reconstruction [##REF##17012374##41##,##REF##17288577##42##]. The inclusion of the most recently released EST data for this large set of phyla led to a dataset including 11,730 amino acid positions and 25 taxa (Additional data file 4).</p> <p>The analysis of this dataset confirmed the branching of chætognaths at the base of the protostomes with significant support values for both the site-homogeneous WAG model and the site-heterogeneous CAT model (bootstrap proportion (BP) of 76 and posterior probability (PP) of 1; Figure ##FIG##3##4a,b##). The inclusion of chætognaths within protostomes is still firmly supported (BP 95, PP 1; Figure ##FIG##3##4##). The inclusion of new taxa strengthens support for both the ecdyozoa and lophotrochozoa clades but the exact relationships within these two clades remain elusive [##REF##17921486##35##,##REF##18222945##36##,##REF##18192187##43##]. Chætognaths and rotifers do not exhibit any peculiar affinities, prompting us to reject the gnathifera hypothesis [##UREF##4##38##]. Conversely, the branching of rotifers is problematic since this phylum is alternatively included in ecdysozoans and lophotrochozoans depending on the use of, respectively, site heterogeneous or homogeneous models (Figure ##FIG##3##4##). Thus, the clustering of platyhelminthes and rotifers in a platyzoa clade is supported by the WAG model but rejected by the CAT model, suggesting that this grouping may be somehow related to long-branch attraction (Figure ##FIG##3##4##). Alternatively, previous studies based on morphology and SSU genes have not argued for the ecdysozoan affinities of rotifers [##UREF##4##38##,##REF##12116426##39##]. Surprisingly, CAT model analysis no longer succeeds in recovering the monophyly of the deuterostomes (Figure ##FIG##3##4b##). Instead, it provides limited support for the successive divergence of chordates and ambulacrarians (echinoderms and hemichordates; PP 0.9; Figure ##FIG##3##4b##). This striking topology was recovered by an independent study using the same heterogeneous CAT model [##REF##18192187##43##] but was neither confirmed by WAG analyses (BP 89 for the monophyly of deuterostomes; Figure ##FIG##3##4a##) nor supported on morphological bases [##REF##17051155##34##,##UREF##4##38##]. One can consider that the two unexpected branchings of rotifers and deuterostomes may be related to some artifact affecting the CAT model, such as sensitivity toward compositional biases [##REF##18234708##44##]. Finally, the placozoan <italic>Trichoplax adherens </italic>surprisingly clustered within the poriferans, as a sister group of the homoscleromorphs (BP 91, PP 0.94; Figure ##FIG##3##4##), although this poriferan status has never been suggested before [##REF##16731622##45##,##REF##16299758##46##]. These challenging hypotheses will be investigated in further studies because they have deep implications for the evolution of metazoans (F Marlétaz <italic>et al</italic>., in progress).</p> <p>Through extended taxon sampling and improved substitution models, these analyses strongly confirm our previous statements about basal-protostome branching of chætognaths [##REF##16890510##18##] and exclude the basal-lophotrochozan hypothesis [##REF##16890509##19##]. Although some areas of bilaterian trees are sometimes incongruent depending on models and inference methods, the position of chætognaths remains remarkably stable throughout our analyses. Furthermore, this branching is not only supported by the presence of GAMT, an unambiguous molecular signature, but also by the posterior Hox genes of chætognaths that are not related to the classes specific to ecdysozoans (Abd-B) or lophotrochozoans (Post1/2) [##UREF##2##16##]. Finally, this topology was also recovered by independent studies involving alternative gene and taxon sampling [##REF##17684563##30##,##REF##17921486##35##,##REF##18192187##43##]. In a broader perspective, the strengthening of their phylogenetic position makes chætognaths a key model for comparative genomics among bilaterians.</p> <title>Genome duplication in the chætognath phylum</title> <p>The clustering of similar sequences indicated that alternative nucleotide forms are present among the transcripts encoding the same protein. Two distinct forms are observed in most cases, although three forms encode some proteins. These forms are separated by a large amount of molecular divergence and can also be distinguished by their different 5' and 3' untranslated regions (UTRs), suggesting that they correspond to different genes (Figure ##FIG##4##5## and Additional data files 5 and 6).</p> <p>Ka/Ks ratios were calculated for all pairs of diverging forms to consider the impact of the nucleotide divergence on the protein sequences. The values of Ka/Ks range from 0.001-0.154 with a median value of 0.004, which confirmed the strong conservation of amino acid sequence despite the large synonymous substitutions observed in some cases (Ks values range from 0.8-75; Table S2 in Additional data file 2). These distinct forms were mainly retrieved for the most highly expressed genes, among which RP genes are prominent (Table ##TAB##0##1##). We verified that the observed molecular divergence could not be explained by the clustering of distant paralogous sequences. For the genes that have clear homologs among metazoans, the sequences of alternative forms always cluster together in phylogenetic analyses and are thus strongly separated from homologous genes of other animals. For instance, the RUX genes have undergone an ancient duplication resulting in the RUX-E and RUX-G paralogs in all metazoans. Interestingly, chætognaths display up to three forms of RUX-E and two forms of RUX-G, all these forms being closely related (Figure S4 in Additional data file 1; Additional data file 7).</p> <p>Such a pattern could be explained by either the duplication of a large set of genes in the genome of chætognaths or, alternatively, it could be explained by the presence of cryptic species within the sampled population. In the first case, the observed differences would be attributed to the divergence between paralogous genes originating through the duplication, where the genome of one individual is expected to contain the two alternative nucleotide forms. In the second case, the observed genetic differences would be caused by the genetic divergence between the orthologous genes of several cryptic species spread among the population, where one individual is thus expected to contain only one of the alternative forms.</p> <p>This cryptic speciation hypothesis may be supported by the strong polymorphism also observed for all genes of the mitochondrial genome, which constitutes an independent lineage from the nuclear genome. For example, cytochrome b (Cytb) transcripts but also cytochrome oxydase I and III are split into distinct forms separated by large molecular distances (Figure ##FIG##4##5c##; Figure S4 in Additional data file 1; Additional data files 8-10), thus testifying to the presence of distinct mitochondrial lineages within the sampled population.</p> <p>To decide between these hypotheses, we designed a PCR screen to survey the alternative forms of selected markers in independent individuals. The genes for RPs L36 and L40 were targeted because they are nuclear genes displaying two alternative forms with the highest number of transcripts in the library (Table S2 in Additional data file 2). The mitochondrial Cytb gene served as an independent reference for the interpretation of results from nuclear genes. The three distinct forms of this strongly diverging mitochondrial gene were surveyed in all the individuals tested (Figure ##FIG##4##5c##). Chætognaths are hermaphroditic and, after fertilization, they store exogenous sperm in their sperm receptacles (Figure ##FIG##4##5a##), which makes it possible to amplify the DNA from another individual. Hence, in order to detect such contamination, we performed independent amplifications on heads, which are considered free from sperm contamination, as well as on the rest of the body, which contained sperm receptacles (Figure ##FIG##4##5a##). The experimental design made it possible to detect alternative forms through the amplification of specific DNA fragments of distinct sizes (Figure ##FIG##4##5b##; Table S3 in Additional data file 2). The PCR products were characterized by sequencing and nucleotide polymorphism was subsequently carefully examined. In addition to their nucleotide divergence in coding sequences, the distinct forms of nuclear genes for RPs L36 and L40 have alternative intron positions and lengths as well as differences in their 5' and 3' UTR regions (Figure ##FIG##4##5b##).</p> <p>Performed on nine individuals, the amplifications revealed the presence of the two forms of the nuclear genes for RPs L36 and L40 in each individual (Table ##TAB##1##2##). Conversely, only one form of the mitochondrial Cytb gene was amplified in each individual with the exception of the body of individual 1, which includes two forms, thus suggesting contamination by exogenous sperm (Table ##TAB##1##2##). The amplification of the divergent nucleotide forms within one individual indicates that the alternative nucleotide forms correspond to paralogous nuclear copies originating through past gene duplication events (Table ##TAB##0##1##). Conversely, the alternative forms of the mitochondrial gene correspond to variation within the population. Because some genes, such as that encoding Translationally controlled tumor protein (TCTP), do not present paralogous copies despite their high expression levels (112 TCTP transcripts in the EST collection; Table S2 in Additional data file 2), we addressed the extent of these duplications in evaluating the quantity of duplicated genes. If the clusters of transcripts encoding the same protein include all the transcripts from alternative paralogous genes and if those paralogous genes have similar levels of expression, the probability that transcripts from these paralogous genes are represented in a given cluster is related to the size of this cluster (see Materials and methods). Hence, all the clusters that include more than six transcripts have at least a 95% chance of including transcripts from the two copies if they exist. Such clusters of transcripts were all checked for paralogous copies through sequence alignments and trees. Paralogs were detected within 35 of the 66 clusters investigated, which suggests that up to 69% of chætognath genes are the products of duplications. These paralogs could have arisen through either a whole genome duplication (WGD) event followed by an extensive gene loss, or several segmental duplication events.</p> <p>The hypothesis of a WGD event is reinforced by the high occurrence of RPs among duplicated genes (Table ##TAB##0##1##). The trend to retain RP genes was previously observed after WGD for <italic>Paramecium tetraurelia</italic>, yeast and plants [##REF##17086204##47##, ####REF##15208398##48##, ##REF##10508730##49####10508730##49##] but is not a common occurrence in small-scale duplications. Conversely, it is difficult to understand why the paralogous genes have been retained after their duplication and maintained under purifying selection as emphasized by Ka/Ks values. This conclusion is in contradiction with the current view of gene destiny after genome duplication, which alternatively predicts that one of the gene duplicates is lost or undergoes the accumulation of substitutions [##REF##10101175##50##]. Using a genome-level dataset, similar findings were made about the strongly duplicated genome of <italic>Paramecium </italic>where the retention of duplicated genes was accounted in part by dosage compensation constraints [##REF##17086204##47##].</p> <p>The most plausible dating is that this duplication occurred before the diversification of the major chætognath lineages. Two copies of SSU and LSU were retrieved in members of the phylum dispersed all over the tree of chætognaths [##REF##16434216##10##, ####REF##8336548##11##, ##REF##9069174##12####9069174##12##]. Moreover, the survey of 226 ESTs available for <italic>Flaccisagitta enflata </italic>also revealed the presence of alternative nucleotide forms for some genes (data not shown), which would confirm that the duplication is not limited to SSU/LSU genes at this taxonomic scale. Further genome data would be required to date the duplication, for instance, in considering the Ks distribution of the set of paralogs [##REF##14757817##51##], and also to definitively state the nature of the duplication through the analysis of synteny in duplicated blocks of the genome. Nevertheless, this preliminary transcriptomic survey stresses the usefulness of the chætognaths to study phylum-level genome duplication events and the destiny of paralogous genes.</p> <title>Population genomics</title> <p>Beyond the molecular divergence between the coding sequences of duplicated paralogous genes, a subsequent survey of the genomic sequences of selected genes revealed that the level of polymorphism is strong within each paralogous gene (Table S4 in Additional data file 2). Multiple nucleotide substitutions as well as insertion/deletion events (indels) occurred within the introns of the four selected nuclear genes (paralogous copies of the genes for both RPs L36 and L40; Additional data files 11-14). Similarly, a large number of substitutions have accumulated in the various mitochondrial genes, thus revealing distinct mitochondrial lineages within the sampled population (Figure ##FIG##4##5c##; Figure S4 in Additional data file 1). However, these strong levels of divergence remain consistent with a population genetic structure because of the regular AT composition and the limited degree of saturation revealed by Ts/Tv ratios, singleton positions being essentially transition substitutions (Table S4 in Additional data file 2; Figure S6 in Additional data file 1).</p> <p>We attempted to determine the origin of this population genetic heterogeneity, which could, for instance, be due to a cryptic speciation or to a past hybridization. For this, the sequences of each individual were compared using phylogenetic trees and indels as discrete informative characteristics (Figure ##FIG##5##6##). For each marker gene, individual sequences split into several major clades supported by strong bootstrap and discrete indel events, which allows unambiguous identification of heterozygous individuals (Figure ##FIG##5##6##). For example, individual 4 is heterozygous for all markers and individuals 6, 9 and 3 are heterozygous for at least one marker. Moreover, the occurrence of several cases of putative recombinations between alleles highlights the heterozygous status of some individuals (individuals 3 and 4, Figure ##FIG##5##6b,d##). Notably, our PCR-based experimental design provided positive evidence only for heterozygosis because two amplifications (head and body) were carried out per individual, yielding 0.5 probability to detect heterozygosity. Heterozygous individuals could thus be even more abundant than observed. These heterozygous cases convincingly demonstrate that a shuffling occurs between the most divergent alleles of each gene, which constitutes strong evidence for interbreeding within the sampled population. This finding definitely excludes the possibility of cryptic speciation within this <italic>S. cephaloptera </italic>population. Alternatively, the panmixy hypothesis was confirmed by the unimodal distribution of pairwise divergences in mismatch analysis, which is consistent with constant population size and excludes a past hybridization event (Figure S6 in Additional data file 1). Finally, the distinct mitochondrial lineages are spread within the population but they are not correlated with any haplotype differentiation at the nuclear level, which is a strong argument against the cryptic speciation hypothesis. This type of mitochondrial diversity was previously discovered for the planktonic species <italic>Sagitta setosa </italic>but was also interpreted with difficulty [##REF##16968273##52##].</p> <p>Strikingly, these comparisons also highlighted molecular divergence between the head and the body of some individuals for each of the five markers investigated (Figure ##FIG##5##6## and Additional data file 4). Such substitutions cannot be explained by a heterozygous status of those individuals because sequences from head and body were firmly clustered in the tree (Figure ##FIG##5##6##). For example, individual 4 exhibits well-separated alleles present in both head and body but intra-individual substitution took place between head and body for both of these alleles (Figure ##FIG##5##6c##). This pattern of substitutions may be explained by the occurrence of somatic mutations during the life of individuals. This interpretation is corroborated by the large extent of intra-individual substitutions in all marker genes and all individuals. Somatic mutations are considered as rare conditions, mainly known from related disorders in humans [##REF##12360233##53##]. Less clear are the evolutionary implications and putative benefits of this phenomenon [##REF##15525396##54##]. They are sometimes suspected to play a prominent role in apoptosis and possibly in the regulation of cell division [##REF##15525396##54##]. Moreover, somatic mutations have been demonstrated to be more widespread in <italic>Drosophila </italic>than in mammals [##REF##17435764##55##], and are sometimes correlated with extensive chromosome rearrangement in the <italic>Drosophila </italic>lineage [##REF##17550304##56##]. However, little is known about the extent and importance of this process in the non-model organisms. In the case of the chætognath, somatic mutation could be due to the high mutation rates that seem to affect both germline and soma and could explain the divergence at the population and individual levels. The possible relationship of these accelerated mutation rates with structural reshaping of the genome after duplication deserves further evaluation.</p> <p>Notably, this level of somatic mutation generates a strong background noise that hinders the accurate interpretation of point mutations related to the diversity of haplotypes. Moreover, traditional hypotheses of population genetics are challenged by our findings: the genetic distances observed between individuals of a single population reach species-level without any evidence for cryptic speciation or past hybridization. In parallel, multiple mitochondrial lineages diverge and are spread and maintained within a single population [##REF##16968273##52##]. If such features are revealed as more widespread than expected, the data collected from numerous population genetic or phylogeographic studies require more cautious interpretation. This observation pleads for an increase of the sampling depth for population genetics, especially through genomic approaches [##REF##14631358##57##].</p> <title>Trans-splicing transcript maturation in chætognaths</title> <p>The survey of the chætognath EST collection detected a common 36 nucleotide motif shared at their 5' end between transcripts from unrelated genes. Selected cases showed that these short motifs are absent from upstream genomic regions of the genes (Figure ##FIG##6##7c##; Additional data files 15-17). These findings strongly suggest that mRNA of <italic>S. cephaloptera </italic>undergoes trans-splicing maturation. This mRNA processing occurs through the spliceosomal transfer of a small RNA molecule at the 5' end of the mRNA from an independent spliced-leader (SL) gene and has been described in numerous species [##REF##15797620##26##]. To exclude the possibility of an artifact, the 226 available EST sequences of <italic>F. enflata</italic>, another chætognath species, were screened using the identified SL sequences from <italic>S. cephaloptera</italic>. This survey led to the recovery of similar SL sequences at the 5' end of several cDNA sequences. As <italic>F. enflata </italic>and <italic>S. cephaloptera </italic>belong to the two main orders of chætognaths [##REF##16434216##10##], the trans-splicing mechanism is likely to be present in the whole chætognath phylum.</p> <p>An extensive study of trans-splicing in the nematode phylum has previously revealed the presence of different forms of SL sequences [##REF##17121468##58##]. Similarly, several different SL sequences were retrieved from <italic>S. cephalotera </italic>(Table ##TAB##2##3##). This number is consistent with the strong level of polymorphism previously observed for coding sequences but two distinct SL forms alone (SL1 and SL2) represent 87% of the SL sequences (Table ##TAB##2##3##). These forms do not exhibit any specificity for the different paralogs since SLs are added randomly to transcripts from these paralogs (Figure ##FIG##6##7##; Additional data files 5 and 6). In an attempt to understand the evolutionary history of trans-splicing within the chætognath phylum, the SL forms of <italic>F. enflata </italic>were compared with those of <italic>S. cephaloptera</italic>. This set of chætognath SL sequences splits into four main forms: two of them, SL1 and SL3, are present in both the species investigated whereas the two other forms, SL2 and SL4, are specific for <italic>S. cephaloptera </italic>and <italic>F. enflata</italic>, respectively (Table ##TAB##2##3##). Otherwise, neither of these chætognath SL forms display similarity with the SL of another phyla. This finding suggests that just as in nematodes, the evolution of alternative forms could have occurred at a relatively reduced taxonomic scale [##REF##17121468##58##].</p> <p>Within the EST collection, 2,914 sequences exhibit SL addition, which represents 30% of nuclear transcripts (Figure ##FIG##7##8a##). Among the SL population, 72% are coding transcripts, of which 46% have a homolog in SwissProt. The clustering of similar coding trans-spliced transcripts indicated that 41% of putative genes undergo SL addition. Furthermore, the relationship between trans-splicing and expression level was tested through the comparison of the number of ESTs per cluster of trans-spliced or non-trans-spliced transcripts. If we posit that this number can be considered as an estimate of the expression level, trans-spliced genes are significantly more expressed than others (Wilcoxon rank test, <italic>p </italic>&lt; 2.2e<sup>-16</sup>). For instance, among the 50 more expressed genes (that is, biggest EST clusters), only two are not trans-spliced. These values suggest that trans-splicing is involved in the regulation of a set of strongly expressed genes responsible for key cellular functions, for example, the RP set.</p> <p>Hence, the set of trans-spliced genes of <italic>S. cephaloptera </italic>was compared with those of other animals that carry out this RNA maturation process (Table ##TAB##3##4##). Particularly, trans-splicing was characterized in the nematode <italic>C. elegans </italic>model, for which a large EST set is available [##REF##17121468##58##,##REF##12560808##59##]. These data were especially useful for extensive comparisons with <italic>S. cephaloptera </italic>(Figure ##FIG##7##8b##). The amount of ESTs with a SL is slightly smaller in <italic>C. elegans </italic>(22%) than in <italic>S. cephaloptera </italic>(29%) and trans-spliced transcripts are missing from the very reduced non-coding ESTs, only 5% of the total EST collection in <italic>C. elegans</italic>. This comparison suggests that trans-splicing is at least as widespread in <italic>S. cephaloptera </italic>as in <italic>C. elegans</italic>. Furthermore, the set of genes affected by trans-splicing in <italic>S. cephaloptera </italic>was comparable with those of <italic>C. elegans</italic>. Among the 119 different genes that are both trans-spliced and annotated using SwissProt with high confidence, 79 have trans-spliced homologous genes in <italic>C. elegans </italic>(Additional data file 3). These values do not include the 78 RP set, which are conserved and trans-spliced in both species. While the molecular actors involved in trans-splicing do not exhibit evolutionary conservation [##REF##15797620##26##], the genes undergoing this kind of processing could otherwise be similar in distantly related species.</p> <p>Trans-splicing has been discovered in a set of organisms spread all over the phylogenetic tree of eukaryotes, but the molecular actors involved in this RNA maturation process do not exhibit any evolutionary conservation: the SL sequences cannot be aligned and the ribonucleoproteic machineries involved in this process have no homologs in other species [##REF##15797620##26##,##REF##15788744##37##,##REF##9297698##60##, ####REF##16822859##61##, ##REF##15972844##62####15972844##62##]. However, the discovery of SL trans-splicing in new animal phyla strongly argues for an ancient origin of this process, especially if these phyla have significant phylogenetic positions, such as chætognaths or acoels. Noticeably, we also found evidence for SL addition in recently released ESTs from the acoels (unpublished observations), worm-like animals that were recently excluded from the protostomes [##REF##17684563##30##]. Moreover, an extensive comparison of available genomic data revealed that the set of genes undergoing trans-splicing is evolutionarily conserved between those species, which is another indication of a putative common evolutionary origin of trans-splicing (Additional data file 3 and Table ##TAB##3##4##).</p> <title>Operonic transcription in <italic>S. cephaloptera</italic></title> <p>In the nematode chætognath <italic>C. elegans </italic>and the tunicate <italic>Ciona intestinalis</italic>, polycistronic mRNA molecules that contain two or more genes are transcribed from operon structures and subsequently resolved by trans-splicing [##REF##17121468##58##,##REF##12560808##59##]. We mined for similar eukaryotic operons in our pool of sequenced bacterial artificial chromosomes (BACs) in looking for clusters of genes with the same transcription orientation and very short intergenic regions. The trans-splicing maturation of these putative operons was confirmed by mapping the tran-spliced ESTs onto these genomic sequences. One operon including RP S14 and PCNA genes was successfully identified in the BAC 35A21 that otherwise contains a large number of genes undergoing trans-splicing (Figure ##FIG##6##7a##; Table S5 in Additional data file 2; Additional data files 15-17). EST mapping and <italic>ab initio </italic>gene prediction suggested that only three bases separate the polyadenylation site of the upstream RP S14 gene from the acceptor splice site at the PCNA downstream gene. This distance is even shorter than that in <italic>C. intestinalis </italic>operons [##REF##16822859##61##] and clearly excludes the possibility for a re-initiation of transcription between the genes (Figure ##FIG##6##7c##). This BAC sequence contains several other predicted genes that are closely clustered, share the same orientation and could thus belong to operonic structures (Figure ##FIG##6##7##; Table S5 in Additional data file 2). The transcripts of gene RP S14 include the two major SL 1 or 2 forms (Figure ##FIG##6##7c##). Similarly, transcripts from both paralogous copies of duplicated genes bear any of the alternative SL forms. Hence, contrary to <italic>C. elegans </italic>for which a SL form is specific for the upstream and the downstream genes in operons [##REF##17121468##58##], no specificity of the SL for gene position in the operon was recovered in <italic>S. cephaloptera </italic>(Figure ##FIG##6##7##).</p> <p>Hence, we have provided arguments for polycistronic transcription in another animal phylum, the chætognaths. Animal operons have previously been described in the nematodes and in the tunicate <italic>C. intestinalis </italic>and were often interpreted as an adaptation to the genome size reduction observed in both these organisms [##REF##15797620##26##]. However, the genome size of <italic>S. cephaloptera </italic>is rather large at 1.05 Gb [##UREF##6##63##] and the transcriptome of chætognaths is conserved among bilaterians, in terms of gene similarity and content. Furthermore, extensive study of genes co-transcribed in <italic>C. elegans </italic>led to the conclusion that the operons constitute functional units in grouping genes involved in similar pathways [##REF##12560808##59##]. Hence, the presence of operonic transcription in the chætognath, a basal protostome, strongly suggests that this mechanism should be more likely considered as an original mode of gene regulation, whose evolutionary importance may have been underestimated. Thanks to the pivotal comparative genomic abilities of chætognaths, the comparison of operon structure and content between nematodes and chætognaths should allow the extension of the functional results obtained in <italic>C. elegans </italic>to a larger set of organisms, including vertebrates. Such an approach would require further genomic data in the chætognath but will likely lead to the characterization of new gene regulatory networks, potentially leading to a better understanding of some genetic disorders [##REF##12560808##59##].</p>
<title>Conclusion</title> <p>Chætognaths have been formerly known for their peculiar morphological characteristics. The fine-scale analysis of transcriptomic data reported here illuminates significant genomic features of this phylum, strengthening its very original status among bilaterians. These genome features may be shuffled into opposite categories: shared ancestral characteristics of bilaterians on one hand and lineage-specific rearrangements and divergences on the other. First, the chætognath transcriptome bears a core set of genes conserved within bilaterians. Together with recent evidence from lophotrochozoans such as annelids [##REF##16311335##64##], this result confirms that the extensive gene loss previously observed within insects or nematodes cannot be extended to the whole protostome lineage. The conservation of such a core gene set was exploited for phylogenomic inference through the definition of RP genes as high-grade marker genes for phylogeny. Subsequent analyses of these marker genes confirmed the basal position of chætognaths among protostomes, a position that has been previously supported [##REF##16890510##18##,##REF##17684563##30##,##REF##18192187##43##]. This position strongly impacts on the understanding of the evolution of embryological characteristics since it rejects some developmental characteristics at the stem of bilaterians, such as enterocoely or secondary mouth opening that were considered to be deuterostome-specific.</p> <p>This interpretation can be extended to the new genome-level characteristics that we have highlighted here. For example, the trans-splicing mechanism of mRNA processing as well as the operonic transcription depicted here were often interpreted as a secondary evolved condition related to peculiar adaptations in animals such as nematodes or tunicates [##REF##15797620##26##]. Contrarily, the occurrence of these mechanisms in chætognaths and acoels suggests that their evolutionary origins may be more ancient than originally expected. Together, these findings strongly argue for the relevance of chætognaths as a comparative genomics model because of their conserved gene set, their original phylogenetic position and their original features such as trans-splicing and operonic transcription. For example, the study of chætognath operons could be a promising approach to identify molecular pathways involved in human genetic diseases [##REF##12560808##59##].</p> <p>Some genetic and genomic features of chætognaths we have illustrated here are unique among the animals. We have discovered a genome duplication event followed by a high retention of duplicated genes that were maintained under strong purifying selection. The current view of the evolution of duplicated genes in animals or plants mainly predicts a loss of genes that have not evolved novel functions [##REF##10101175##50##,##REF##15568988##65##]. This view is contradicted by recent evidence from protists [##REF##17086204##47##], plants [##REF##15800040##25##] and the whole fungal kingdom [##REF##17805289##66##] that suggests new modalities of genome evolution and, in particular, strong retention and constraint of former duplicated genes. Hence, the chætognath genome duplication deserves further documentation and investigation as one possible special case among bilaterians. Another unique feature of chætognaths is the high level of divergence within the population. We demonstrated that this divergence is related to neither cryptic speciation nor past hybridization but conversely revealed very high mutation rates within the germ line of individuals but also within their soma. These observations have deep implications for the correct interpretation of numerous population genetics studies. This apparent genome plasticity with duplications and high mutation rates is in contrast to the strong morphological conservation of the phylum that has continued nearly unchanged since the early Cambrian period [##REF##12364798##8##,##REF##17254986##9##]. This uncoupling between morphological conservation and genome plasticity makes us wonder how the gene regulatory networks responsible for the establishment of the body organization of these animals could have remained stable despite such a genome dynamic?</p> <p>In summary, these findings promote the study of chætognath for orienting not only morphological but also genomic characteristics among bilaterians. With the recent rise of lophotrochozoans as promising model systems [##REF##17448990##67##], the chætognath phylum represents the next step of investigation within protostomes with its striking combination of distinct ancestral and derived features that could illuminate the evolution of bilaterians.</p>
<p>This is an open access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <p>The chætognath transcriptome reveals unusual genomic features in the evolution of this protostome and suggests that it could be used as a model organism for bilaterians.</p> <title>Background</title> <p>The chætognaths (arrow worms) have puzzled zoologists for years because of their astonishing morphological and developmental characteristics. Despite their deuterostome-like development, phylogenomic studies recently positioned the chætognath phylum in protostomes, most likely in an early branching. This key phylogenetic position and the peculiar characteristics of chætognaths prompted further investigation of their genomic features.</p> <title>Results</title> <p>Transcriptomic and genomic data were collected from the chætognath <italic>Spadella cephaloptera </italic>through the sequencing of expressed sequence tags and genomic bacterial artificial chromosome clones. Transcript comparisons at various taxonomic scales emphasized the conservation of a core gene set and phylogenomic analysis confirmed the basal position of chætognaths among protostomes. A detailed survey of transcript diversity and individual genotyping revealed a past genome duplication event in the chætognath lineage, which was, surprisingly, followed by a high retention rate of duplicated genes. Moreover, striking genetic heterogeneity was detected within the sampled population at the nuclear and mitochondrial levels but cannot be explained by cryptic speciation. Finally, we found evidence for trans-splicing maturation of transcripts through splice-leader addition in the chætognath phylum and we further report that this processing is associated with operonic transcription.</p> <title>Conclusion</title> <p>These findings reveal both shared ancestral and unique derived characteristics of the chætognath genome, which suggests that this genome is likely the product of a very original evolutionary history. These features promote chætognaths as a pivotal model for comparative genomics, which could provide new clues for the investigation of the evolution of animal genomes.</p>
<title>Abbreviations</title> <p>BAC, bacterial artificial chromosome; BP, bootstrap proportion; Cytb, cytochrome b; EST, expressed sequence tag; GAMT, guanidinoacetate N-methyltransferase; indels, insertion or deletion events; LSU, large subunit of ribosomal RNA; PP, posterior probability; RP, ribosomal protein; 
SL, splice-leader; SSU, small subunit of ribosomal RNA; UTR, untranslated region; WGD, whole genome duplication.</p> <title>Authors' contributions</title> <p>YLP and FM conceived the project and performed the analyses focused on phylogenomics, characterization of duplication and trans-splicing. YLP and XC were in charge of the construction of the cDNA and BAC libraries. YLP, FM, and AG carried out the experiments related to the study of polymorphism in the population. YP ensured the collection of animals required for construction of libraries and population genetics studies. PW supervised the sequencing project that SS, CD and GG carried out. FM wrote the paper and all authors agreed on the final version.</p> <title>Additional data files</title> <p>The following additional data are available. Additional data file ##SUPPL##0##1## includes supplementary Figures S1-S6 with legends: EST library statistics (Figure S1), Gene Ontology annotation (Figure S2), alignment of GAMT in selected taxa (Figure S3), trees illustrating the divergence between copies of nuclear and mitochondrial genes in the library (Figure S4), transition versus transversion ratios for the five targeted genes (Figure S5), and mismatch analyses and testing for the five targeted genes (Figure S6). Additional data file ##SUPPL##1##2## includes supplementary Tables S1-S5 with legends: composition of the composite taxa involved in phylogenomic analyses (Table S1), list of all clusters of transcripts corresponding to alternative forms or not with Ka/Ks ratios (Table S2), list of primers employed for PCR amplifications of alternative forms (Table S3), comparison of molecular evolution trends for four genes retrieved in nine individuals (Table S4), and annotation of the BAC 35A21 (Table S5). Additional data file ##SUPPL##2##3## lists the trans-spliced genes conserved between <italic>S. cephaloptera </italic>and <italic>C. elegans </italic>andhomologous to a SwissProt entry (Figure S6 in Additional data file ##SUPPL##0##1##). Additional data file ##SUPPL##3##4## is the concatenated gblocked alignment that clusters 77 RPs, yielding 24 taxa and 11,730 shuffled amino acid positions. This dataset was employed for phylogenomic reconstruction. Additional data files ##SUPPL##4##5## and ##SUPPL##5##6## are the transcript sequences of RP L36 (Additional data file ##SUPPL##4##5##) and RP L40 (Additional data file ##SUPPL##5##6##) exhibiting the alternative forms corresponding to paralogous copies and the position of primers designed to specifically amplify these forms. Additional data file ##SUPPL##6##7## is the alignment of amino acid sequences of RUX paralogs and variants from <italic>S. cephaloptera </italic>and some other eukaryotic taxa. Additional data file ##SUPPL##7##8## is the alignment of Cytb displaying selected variants retrieved from the EST library and the coding sequences of Cytb from nine individuals showing the differences between genetic lineages within the population (primers employed for the specific amplifications are included). Additional data files ##SUPPL##8##9## and ##SUPPL##9##10## are the alignments of cytochrome oxydase I and III transcripts retrieved from the EST library. Additional data files ##SUPPL##10##11##, ##SUPPL##11##12##, ##SUPPL##12##13## and ##SUPPL##13##14## are the alignments of genomic sequence of paralogs of RP L36 (paralogs 1 and 2, Additional data files ##SUPPL##10##11## and ##SUPPL##11##12##) and RP L40 (paralogs 1 and 2, Additional data files ##SUPPL##12##13## and ##SUPPL##13##14##) from 9 individuals that have been amplified from the body (B) or head (H) of individuals as indicated in the sequence name. Additional data files ##SUPPL##14##15##, ##SUPPL##15##16## and ##SUPPL##16##17## are EST/BAC alignments for four trans-spliced genes present in BAC 35A21, whose transcripts are represented in the EST library. Additional data file ##SUPPL##14##15## contains the RP S14 and PCNA genes gathered in an operonic structure and eight matching ESTs. Additional data file ##SUPPL##15##16## contains the COP9 gene and two matching ESTs. Additional data file ##SUPPL##16##17## contains an unidentified gene matching EST clones 6YJ23. All these alignments clearly exhibit lack of a SL motif in the genomic sequence.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Bastien Boussau and Cédric Finet for their insightful comments on the manuscript. We also thank Eliane Alivon for her help in the molecular biology experiments, Susan Cure and Andy Saurin for correcting the manuscript. This work was supported by the GIS Génomique Marine (CNRS), the Consortium National de Recherche en Génomique, BQR from Université de Provence and Université de la Méditerranée. FM is supported by a BDI CNRS Fellowship.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Overall composition of the EST collection. The annotation of transcripts is based on SwissProt (score &gt;150) and led to identification of mitochondrial genes. The conceptual translation of ESTs allowed detection of those that include coding sequences. The large portion of non-coding polyadenylate nuclear transcripts and RPs among nuclear transcripts is the most prominent aspect of this distribution as well as the unexpected presence of mitochondrial rRNAs (12 and 16S) related to their polyadenine stretches.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Visualization of relative similarity between the transcriptome of <italic>S. cephaloptera </italic>and <bold>(a) </bold>selected species or <bold>(b) </bold>corresponding clades: <italic>H. sapiens </italic>as a deuterostome, <italic>D. melanogaster </italic>as an ecdsyzoan and <italic>L. rubellus </italic>as a lophotrochozoan. The graphs are based on whole transcriptome Blast comparisons and the plotting of respective Blast scores was performed using Simitri [77] (cut-off score 150). Genes at the center of the plot are equally related to the three databases and hence represent valuable phylogenetic markers, whereas genes attracted by a node share a greater similarity with the corresponding database. Genes on the edge do not have a match in the database from the opposite vertex and those on the vertex only have a match in the corresponding database; these two types of genes constitute candidates for signature genes that have possibly been lost in a peculiar lineage. The color scale indicates the relevancy of scores.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>RP minimization of missing data in EST-based phylogenomic datasets. Dataset completeness was estimated for datasets composed of 78 RPs (red) or 115 other genes (green) retrieved from EST collections of a large range of sizes.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>The basal-protostome branching of chætognaths is confirmed through improved inference methods and expanded taxon sampling. A RP alignment of 11,730 positions (after GBlock filtration; see Additional data file 4) was analyzed using two classes of models. <bold>(a) </bold>Site-homogeneous model (WAG) implemented in a maximum-likelihood framework (PhyML [80] and Treefinder [81]). Similar topology and maximal posterior probabilities were obtained with Bayesian analyses using the same model (MrBayes). <bold>(b) </bold>Site-heterogeneous model (CAT) implemented in a bayesian framework (Phylobayes [79]). Plain colored circles denote nodes for which significant support values were obtained (likelihood ratio statistics based on expected-likelihood weights (LR-ELW) &gt;0.95 for site-homogenous and PP &gt;0.95 for site-heterogenous). Support values are indicated for selected nodes: LR-ELW statistics and bootstrap (bold type) for maximum likelihood (ML) using the WAG model and posterior probabilities for Bayesian inference using the CAT model.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>Alternative forms of selected markers amplified by PCR in order to assess the origin of polymorphism. <bold>(a) </bold>Localization of sperm within sperm receptacles (SR) and sperm ducts (SD) in the body of chætognath <italic>S. cephaloptera </italic>along with ovaries (Ov) and testis (Te). The double arrow indicates that head and body of individuals were split to perform independent PCR amplifications with the purpose of detecting possible contamination from the sperm genome. <bold>(b) </bold>Paralogous copies of nuclear genes RP L36 and L40 with their intron positions and average lengths, which are distinct in both cases (Additional data files 5 and 6). The names and positions of primers used for the amplification are also specified (Table S3 in Additional data file 2). <bold>(c) </bold>Relationships between alternative copies of Cytb retrieved within the ESTs with the three different forms detected by the designed primers (Additional data file 8). Boostrap proportions are indicated for selected nodes.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>Relationships between haplotypes of nine individuals, including distinct head and body sequences for four marker genes, including two pairs of paralogous sequences: <bold>(a) </bold>L36 form 1; <bold>(b) </bold>L36 form 2; <bold>(c) </bold>L40 form 1; <bold>(d) </bold>L40 form 2. The indels are plotted onto the branches (green lines). Noticeably, some individuals display mixed sequences from different haplotypes, which is explained by recombination events between alleles (purple star). The substitutions occurring between copies of the same allele in the head and body of individuals (blue branches) are assumed to be somatic mutations. These neighbor-joining trees were inferred assuming kimura 2 parameter distances from Additional data files 7-10. Boostrap proportions are indicated for selected nodes.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>Identified <italic>S. cephaloptera </italic>operon within the BAC 35YA21. <bold>(a) </bold>Structure of the 158 kb BAC 35A21, including the predicted genes and the mapping of ESTs that bear SLs (purple). Detailed EST/BAC alignments are provided as Additional data files 15-17. <bold>(b) </bold>Detailed structure of the identified chætognath operon with RP S14 and PCNA genes and the corresponding ESTs (purple) that exhibit SL sequences. UTR (genomic, light orange; EST, light purple) and coding sequences (genomic, orange; EST, purple). <bold>(c) </bold>Alignments of selected regions of the operon, beginning and end of genes showing genomic DNA (orange) and transcripts with their alternative SL forms. The very short distance encountered between the end and beginning of the two genes argues for polycistronic transcription.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p>Categories of trans-spliced transcripts for chætognath <bold>(a) </bold><italic>S. cephaloptera </italic>and <bold>(b) </bold>nematode <italic>C. elegans</italic>. The presence of a SL sequence is related to coding properties and homologous matches in SwissProt (score &gt;150) of the sequences. <italic>C. elegans </italic>exhibits less non-coding transcripts than <italic>S. cephaloptera</italic>.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Occurrence of paralogous gene copies for ribosomal and non-ribosomal genes</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center" colspan="3">Ribosomal protein genes</td><td align="center" colspan="3">Other genes</td></tr><tr><td/><td colspan="3"><hr/></td><td colspan="3"><hr/></td></tr><tr><td align="left">Inferred duplicates</td><td align="center">Gene number</td><td align="center">Percent selected genes</td><td align="center">Median EST number</td><td align="center">Gene number</td><td align="center">Percent selected genes</td><td align="center">Median EST number</td></tr></thead><tbody><tr><td align="left">1</td><td align="center">17</td><td align="center">37%</td><td align="center">12</td><td align="center">14</td><td align="center">70%</td><td align="center">16</td></tr><tr><td align="left">2</td><td align="center">28</td><td align="center">61%</td><td align="center">23</td><td align="center">3</td><td align="center">15%</td><td align="center">6</td></tr><tr><td align="left">3</td><td align="center">1</td><td align="center">2%</td><td align="center">23</td><td align="center">3</td><td align="center">15%</td><td align="center">8.5</td></tr><tr><td align="left">Total</td><td align="center">46</td><td align="center">-</td><td align="center">20.5</td><td align="center">20</td><td align="center">-</td><td align="center">8.5</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Distinct forms recovered from PCR amplification performed on heads and bodies of ten individuals for alternative marker genes</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center" colspan="2">L36</td><td align="center" colspan="2">L40</td><td align="center" colspan="3">Cytb</td></tr><tr><td/><td colspan="2"><hr/></td><td colspan="2"><hr/></td><td colspan="3"><hr/></td></tr><tr><td align="left">Individual</td><td align="center">1</td><td align="center">2</td><td align="center">1</td><td align="center">2</td><td align="center">1</td><td align="center">2</td><td align="center">3</td></tr></thead><tbody><tr><td align="left">1</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">/+</td><td align="center">+/+</td><td/></tr><tr><td align="left">2</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td/><td/></tr><tr><td align="left">3</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/2</td><td align="center">+</td><td/><td/></tr><tr><td align="left">4</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">2/2</td><td align="center">+/2</td><td align="center">+/+</td><td/><td/></tr><tr><td align="left">5</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td/><td/><td align="center">+/+</td></tr><tr><td align="left">6</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/2</td><td align="center">+/+</td><td/><td/></tr><tr><td align="left">7</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td/><td align="center">+/+</td><td/></tr><tr><td align="left">8</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td/><td/></tr><tr><td align="left">9</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td align="center">+/+</td><td/><td/></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Splice-leader isoforms isolated in two chætognath species: <italic>S. cephaloptera </italic>and <italic>F. enflata</italic></p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Form</td><td align="left">ID</td><td align="left">Sequence</td><td align="left">Species</td><td align="center">Abundance</td></tr></thead><tbody><tr><td align="left"><bold>Form 1</bold></td><td align="left"><bold>SL1</bold></td><td align="left">GGAAGCTAAAATTCTTTTA--TTTGCTT-AATTAAA</td><td align="left">Both</td><td align="center">&lt;</td></tr><tr><td/><td/><td/><td/><td/></tr><tr><td align="left"><bold>Form 2</bold></td><td align="left"><bold>SL2.0</bold></td><td align="left">GGAAGCTTATAATTGAGTAGTTTCAATTTGTTTAAA</td><td align="left">Both</td><td align="center">70.83%</td></tr><tr><td/><td align="left">SL2.1</td><td align="left">GGAAGCTTATAATTGAGTGGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.2</td><td align="left">GGAAGCTTATAATTGAGCAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.3</td><td align="left">GGAAGCTTATAATTGACTATTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.4</td><td align="left">GGAAGCTTATAATTGAGTATTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">1.44%</td></tr><tr><td/><td align="left">SL2.5</td><td align="left">GGAAGCTTATAATTGCGTATTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.6</td><td align="left">GGAAGCTTATAATTGCGTAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.7</td><td align="left">GGAAGCTTATAATTCATTAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.8</td><td align="left">GGAAGCTTATAATTCAGTAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.9</td><td align="left">GGAAGCTTATAATTGATTAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">2.17%</td></tr><tr><td/><td align="left">SL2.10</td><td align="left">GGAAGCTTATAACTGTTTAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.11</td><td align="left">GGAAGCTTATAATTGACTAGTTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL2.12</td><td align="left">GGAAGCTTATAATTGAGTAGCTTCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td/><td/><td/><td/></tr><tr><td align="left"><bold>Form 3</bold></td><td align="left">SL3.1</td><td align="left">GGAAGCCAAT-TTCTACTA-CTTCACTT-GTTTAAA</td><td align="left"><italic>F. enflata</italic></td><td align="center">-</td></tr><tr><td/><td align="left">SL3.2</td><td align="left">GGAAGCTAAT-TTCTACTA-CTTCACTT-GTTTAAA</td><td align="left"><italic>F. enflata</italic></td><td align="center">-</td></tr><tr><td/><td align="left">SL3.3</td><td align="left">GGAAGCTAAT-ATCTACTA-CTTCACTTTGTTTAAA</td><td align="left"><italic>F. enflata</italic></td><td align="center">-</td></tr><tr><td/><td/><td/><td/><td/></tr><tr><td align="left"><bold>Form 4</bold></td><td align="left">SL4.1</td><td align="left">GGAAGCTTATTATCAAGTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.2</td><td align="left">GGAAGCTTATTATTAAGTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">16.91%</td></tr><tr><td/><td align="left">SL4.3</td><td align="left">GGAAGCTTATTATTAAGTACTACCAAATTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.4</td><td align="left">GGAAGCTTATTATTAAGTACTACCAGTTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.5</td><td align="left">GGAAGCTTATTATTACGTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.6</td><td align="left">GGAAGCTTATTATTAATTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.7</td><td align="left">GGAAGCTTATTATTATTTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.8</td><td align="left">GGAAGCTTATTATTAACTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr><tr><td/><td align="left">SL4.9</td><td align="left">GGAAGCTTATTATTCAGTACTACCAATTTGTTTAAA</td><td align="left"><italic>S. cephaloptera</italic></td><td align="center">&lt;</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Distribution and modalities of trans-splicing among metazoans</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Phylum: species</td><td align="center">Trans-splicing</td><td align="center">Percent trans-spliced genes</td><td align="center">Number of SL forms</td><td align="center">Operons</td><td align="center">Percent genes in operons</td><td align="center">References</td></tr></thead><tbody><tr><td align="left">Cnidaria: <italic>Hydra vulgaris</italic></td><td align="center">+</td><td align="center">Very few*</td><td align="center">2</td><td align="center">?</td><td align="center">NA</td><td align="center">[85]</td></tr><tr><td align="left">Urochordata (tunicate): <italic>Ciona intestinalis</italic></td><td align="center">+</td><td align="center">~2%*</td><td align="center">1</td><td align="center">+</td><td align="center">~2%</td><td align="center">[61]</td></tr><tr><td align="left">Chætognatha: <italic>Spadella cephaloptera</italic>; <italic>Flaccisagitta enflata</italic></td><td align="center">+</td><td align="center">~55%</td><td align="center">2*</td><td align="center">+</td><td align="center">ND</td><td align="center">This study</td></tr><tr><td align="left">Acoela: <italic>Convoluta pulchra</italic></td><td align="center">+</td><td align="center">&gt;5%</td><td align="center">1</td><td align="center">?</td><td align="center">ND</td><td align="center">This study</td></tr><tr><td align="left">Rotifera: <italic>Adineta ricciae</italic>; <italic>Philodina </italic>sp.</td><td align="center">+</td><td align="center">ND</td><td align="center">1</td><td align="center">ND</td><td align="center">ND</td><td align="center">[37]</td></tr><tr><td align="left">Platyhelminth: <italic>Schistosoma mansoni</italic>; seven more species<sup>‡</sup></td><td align="center">+</td><td align="center">&gt;3%<sup>†</sup></td><td align="center">1-2</td><td align="center">+</td><td align="center">ND</td><td align="center">[60,62]</td></tr><tr><td align="left">Nematoda: <italic>Caenothabditis elegans</italic></td><td align="center">+</td><td align="center">~50%</td><td align="center">2</td><td align="center">+</td><td align="center">~20%</td><td align="center">[58,59]</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional data file 1</title><p>Figure S1: EST library statistics. Figure S2: Gene Ontology annotation. Figure S3: alignment of GAMT in selected taxa. Figure S4: trees illustrating the divergence between copies of nuclear and mitochondrial genes in the library. Figure S5: transition versus transversion ratios for the five targeted genes. Figure S6: mismatch analyses and testing for the five targeted genes.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional data file 2</title><p>Table S1: composition of the composite taxa involved in phylogenomic analyses. Table S2: list of all clusters of transcripts corresponding to alternative forms or not with Ka/Ks ratios. Table S3: list of primers employed for PCR amplifications of alternative forms. Table S4: comparison of molecular evolution trends for four genes retrieved in nine individuals. Table S5: annotation of the BAC 35A21.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional data file 3</title><p>Trans-spliced genes conserved between <italic>S. cephaloptera </italic>and <italic>C. elegans </italic>and homologous to a SwissProt entry (Figure S6 in Additional data file ##SUPPL##0##1##).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional data file 4</title><p>This dataset was employed for phylogenomic reconstruction.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional data file 5</title><p>Transcript sequences of RP L36 exhibiting the alternative forms corresponding to paralogous copies and the position of primers designed to specifically amplify these forms.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional data file 6</title><p>Transcript sequences of RP L40 exhibiting the alternative forms corresponding to paralogous copies and the position of primers designed to specifically amplify these forms.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional data file 7</title><p>Alignment of amino acid sequences of RUX paralogs and variants from <italic>S. cephaloptera </italic>and some other eukaryotic taxa.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional data file 8</title><p>Primers employed for the specific amplifications are included.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional data file 9</title><p>Alignment of cytochrome oxydase I transcripts retrieved from the EST library.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S10"><caption><title>Additional data file 10</title><p>Alignment of cytochrome oxydase III transcripts retrieved from the EST library.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S11"><caption><title>Additional data file 11</title><p>Sequences were amplified from the body (B) or head (H) of individuals as indicated in the sequence name.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S12"><caption><title>Additional data file 12</title><p>Sequences were amplified from the body (B) or head (H) of individuals as indicated in the sequence name.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S13"><caption><title>Additional data file 13</title><p>Sequences were amplified from the body (B) or head (H) of individuals as indicated in the sequence name.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S14"><caption><title>Additional data file 14</title><p>Sequences were amplified from the body (B) or head (H) of individuals as indicated in the sequence name.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S15"><caption><title>Additional data file 15</title><p>All the alignments in Additional data files ##SUPPL##14##15##, ##SUPPL##15##16##, ##SUPPL##16##17## clearly exhibit lack of a SL motif in the genomic sequence.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S16"><caption><title>Additional data file 16</title><p>All the alignments in Additional data files ##SUPPL##14##15##, ##SUPPL##15##16##, ##SUPPL##16##17## clearly exhibit lack of a SL motif in the genomic sequence.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S17"><caption><title>Additional data file 17</title><p>All the alignments in Additional data files ##SUPPL##14##15##, ##SUPPL##15##16##, ##SUPPL##16##17## clearly exhibit lack of a SL motif in the genomic sequence.</p></caption></supplementary-material>
<table-wrap-foot><p>A plus sign indicates that one copy was amplified and a numeral indicates the number of copies if more than one were amplified (size distinct alleles). The copies amplified in heads and bodies are separated by a slash (head/body).</p></table-wrap-foot> <table-wrap-foot><p>The abundance of the forms corresponds to the rate of the forms among all SLs and is calculated only for <italic>S. cephaloptera </italic>sequences. A 'less-than' sign indicates that the abundance is below 1%.</p></table-wrap-foot> <table-wrap-foot><p>*This figure was obtained by the analysis of a sample of several cDNA libraries sequenced from the 5' end. <sup>†</sup>This figure was obtained from 3' end sequenced cDNA libraries and may be greatly underestimated. <sup>‡</sup>For a precise list of these species see [62]. NA, not applicable; ND, not determined.</p></table-wrap-foot>
<graphic xlink:href="gb-2008-9-6-r94-1"/> <graphic xlink:href="gb-2008-9-6-r94-2"/> <graphic xlink:href="gb-2008-9-6-r94-3"/> <graphic xlink:href="gb-2008-9-6-r94-4"/> <graphic xlink:href="gb-2008-9-6-r94-5"/> <graphic xlink:href="gb-2008-9-6-r94-6"/> <graphic xlink:href="gb-2008-9-6-r94-7"/> <graphic xlink:href="gb-2008-9-6-r94-8"/>
<media xlink:href="gb-2008-9-6-r94-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S4.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S5.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S6.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S7.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S8.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S9.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S10.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S11.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S12.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S13.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S14.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S15.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S16.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-6-r94-S17.fast" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Hyman"], "given-names": ["LH"], "source": ["The Invertebrates, Vol 5, Smaller Coelomate Groups"], "year": ["1959"], "publisher-name": ["New York: McGraw-Hill"]}, {"surname": ["Kapp"], "given-names": ["H"], "article-title": ["The unique embryology of Chaetognatha."], "source": ["Zool Anz"], "year": ["2000"], "volume": ["239"], "fpage": ["263"], "lpage": ["266"]}, {"surname": ["Matus", "Halanych", "Martindale"], "given-names": ["DQ", "KM", "MQ"], "article-title": ["The "], "italic": ["Hox", "Flaccisagitta enflata"], "source": ["Integr Comp Biol"], "year": ["2007"], "volume": ["47"], "fpage": ["854"], "lpage": ["864"], "pub-id": ["10.1093/icb/icm077"]}, {"surname": ["Van Pilsum", "Stephens", "Taylor"], "given-names": ["JF", "GC", "D"], "article-title": ["Distribution of creatine, guanidinoacetate and the enzymes for their biosynthesis in the animal kingdom. Implications for phylogeny."], "source": ["Biochem J"], "year": ["1972"], "volume": ["126"], "fpage": ["325"], "lpage": ["345"]}, {"surname": ["Nielsen"], "given-names": ["C"], "source": ["Animal Evolution: Interelationships of the Living Phyla"], "year": ["2001"], "publisher-name": ["New York: Oxford University Press"]}, {"surname": ["Halanych"], "given-names": ["KM"], "article-title": ["The new view of animal phylogeny."], "source": ["Annu Rev Ecol Evol Syst"], "year": ["2004"], "volume": ["35"], "fpage": ["229"], "lpage": ["256"], "pub-id": ["10.1146/annurev.ecolsys.35.112202.130124"]}, {"article-title": ["Animal Genome Size Database"]}, {"article-title": ["Ensembl Genome Browser"]}, {"article-title": ["NCBI EST Database (dbEST)"]}, {"article-title": ["JGI Eukaryotic Genomics"]}, {"article-title": ["NCBI Trace Archive"]}]
Genome Biol. 2008 Jun 4; 9(6):R94