The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.This work was supported by the following grants: NHGRIU54HG003273 to R.A.G; EU Marie Curie ITN #215781 “Evonet” to M.A.; a Wellcome Trust Value in People (VIP) award to C.B. and Wellcome Trust graduate studentship WT089615MA to J.E.G; Marine rhythms of Life” of the University of Vienna, an FWF (http://www.fwf.ac.at/) START award (#AY0041321) and HFSP (http://www.hfsp.org/) research grant (#RGY0082/2010) to KT-­‐R; MFPL Vienna International PostDoctoral Program for Molecular Life Sciences (funded by Austrian Ministry of Science and Research and City of Vienna, Cultural Department -­‐Science and Research to T.K; Direct Grant (4053034) of the Chinese University of Hong Kong to J.H.L.H.; NHGRI HG004164 to G.M.; Danish Research Agency (FNU), Carlsberg Foundation, and Lundbeck Foundation to C.J.P.G.; U.S. National Institutes of Health R01AI55624 to J.H.W.; Royal Society University Research fellowship to F.M.J.; P.D.E. was supported by the BBSRC via the Babraham Institute;This is the final version of the article. It first appeared from PLOS via http://dx.doi.org/10.1371/journal.pbio.100200

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history

Ancestral protein kinases are extensively lost during arthropod evolution.

<p><i>S. maritima</i> is an exception and retains the largest number of ancestral kinases. Numbers of kinase subfamilies in selected species are shown in parentheses after species names. The gains, losses, and inferred content of common ancestors are listed on internal branches. Kinases found in at least two species from human, <i>C. elegans</i> and <i>Nematostella vectenesis</i> were used as an outgroup.</p

Presence and absence of immunity genes in different arthropods.

<p>Counts of immune genes are shown for <i>S. maritima</i>, <i>D. pulex</i><a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.1002005#pbio.1002005-McTaggart1" target="_blank">[131]</a>, <i>A. mellifera</i><a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.1002005#pbio.1002005-Evans1" target="_blank">[86]</a>, <i>T. castaneum</i>, <i>Anopheles gambiae</i>, and <i>D. melanogaster</i><a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.1002005#pbio.1002005-Dasmahapatra1" target="_blank">[132]</a>. ∼, identity of the gene is uncertain; -, not investigated.</p

Expansion of chemosensory receptor families.

<p>(A) Phylogenetic relationships among <i>S. maritima</i> (Smar), <i>I. scapularis</i> (Isca), <i>D. pulex</i> (Dpul), and a few insect GRs that encode for sugar, fructose, and carbon dioxide receptors (Dmel, <i>D. melanogaster</i>, and Amel, <i>A. mellifera</i>). (B) Phylogenetic relationships among <i>S. maritima</i>, <i>I. scapularis</i>, and a few <i>D. melanogaster</i> IRs and IgluR genes (the suffix at the end of the protein names indicates: i, incomplete and p, pseudogene).</p

Dscam diversity caused either by gene and/or exon duplication in different Metazoa.

<p>^aOnly canonical Dscam paralogues were considered. ^bIn <i>D. melanogaster</i> and <i>D. pulex</i> the paralogue Dscam-L2 has two Ig7 alternative coding exons. ^cPotential number of Dscam isoforms, circulating in one individual, produced by mutually exclusive alternative splicing of duplicated exons.</p

Frequency histogram of CpG_(o/e) observed in <i>S. maritima</i> gene bodies.

<p>The y-axis depicts the number of genes with the specific CpG_(o/e) values given on the x-axis. The distribution of CpG_(o/e) in <i>S. maritima</i> is a trimodal distribution, with a low-CpG_(o/e) peak consistent with the presence of historical DNA methylation in <i>S. maritima</i> and the presence of a high CpG_(o/e) peak. The data underlying this plot are available in <a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.1002005#pbio.1002005.s068" target="_blank">File S4</a>.</p

Instances of homeobox gene clustering and linkage.

<p>Instances of homeobox gene clustering and linkage.</p

Conserved macro synteny signal between <i>S. maritima</i> and the chordate lancelet <i>B. floridae</i> clustered into ancestral linkage groups.

<p>Each dot represents a pair of genes, one in <i>B. floridae</i>, one in <i>S. maritima</i>, assigned to the same gene family by our orthology analysis. The ancestral linkage group identifiers refer to groups of scaffolds from the <i>S. maritima</i> (SmALG) or <i>B. floridae</i> (BfALG) assemblies, as detailed in <a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.1002005#pbio.1002005.s066" target="_blank">File S2</a>. The identification of ALGs is described in the SI. Note that two <i>S. maritima</i> scaffolds were divided across ALGs, and so appear multiple times in <a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.1002005#pbio.1002005.s066" target="_blank">File S2</a>.</p

Arthropod phylogenetic tree (with nematode outgroup) showing selected events of gene loss, gene gain, and gene family expansions.

<p>Main taxa are listed on the tips, with representative species for which there is a fully sequenced genome listed below. Major nodes are also named. Data from the genome of <i>S. maritima</i> allow us to map when in arthropod evolution these events occurred, even when these events did not occur on the centipede lineage. A plausible node for the occurrence of each event is marked and colour-coded, with the possible range marked with a thin line of the same colour. The events, listed from left to right are: (1) Dscam alternative splicing as a strategy for increasing immune diversity is known from <i>D. melanogaster</i>, as well as the crustacean <i>D. pulex</i>, and thus probably evolved in the lineage leading to pancrustacea, after the split from centipedes. (2) Several wnt genes have been lost in holometabolous insects, leaving only seven of the 13 ancestral families. This loss occurred gradually over arthropod evolution, but reached its peak at the base of the Holometabola. (3) Selenoproteins are rare in insects. The presence of a large number of selenoproteins in <i>S. maritima</i> as well as in other non-insect arthropods suggests that the loss of many selenoproteins occurred at the base of the Insecta. (4) Expansion of chemosensory gene families occurred independently in different arthropod lineages as they underwent terrestrialisation. The OR family is expanded in insects only. (5) Chemosensory genes of the GR and IR genes have undergone a lineage specific expansion in the genome of <i>S. maritima</i>. As these are probably also linked with terrestrialisation we suggest that this expansion happened at the base of the Chilopoda, but it could have also occurred later in the lineage leading to <i>S. maritima</i>. (6) Cuticular proteins of the RR-1 family are numerous in the <i>S. maritima</i> genome. They are found in other arthropods, but not in chelicerates nor in any non-arthropod species. This suggests that the RR-1 subfamily evolved at the base of the Mandibulata. (7) The genome of <i>S. maritima</i> has a large complement of wnt genes, but is missing <i>wnt8</i>. Since this gene is found in the Diplopod <i>G. marginata</i> (a species without a fully sequenced genome), the loss most likely occurred at the base of the Chilopoda. (8) Unlike the situation in <i>D. melanogaster</i>, immune diversity in the <i>S. maritima</i> genome is achieved through multiple copies of the Dscam gene. This expansion of the family could have happened at any time after the split between Myriapoda and Pancrustacea. (9) Both circadian rhythm genes and many light receptors are missing in <i>S. maritima</i>. These losses are most likely due to the subterranean life style of geophilomorph centipedes and are probably specific to this group. However, we cannot rule out the possibility that they were lost somewhere in the lineage leading to myriapods. (10) The existence of JH signalling in <i>S. maritima</i> as well as in all other arthropods studied to date strengthens the idea that this signalling system evolved with the exoskeleton of arthropods, though its origins could be even more ancient and date back to the origin of moulting at the base of the Ecdysozoa.</p