Sequencing and analysis of the Mediterranean amphioxus (Branchiostoma lanceolatum) transcriptome

Discovering the deep evolutionary roots of serum amyloid A protein family

Deep evolutionary origin of the conserved animal serum amyloid A (SAA) apolipoprotein family leading to yet unknown highly similar SAA-like sequences occurring in certain bacterial genomes is demonstrated in this contribution. Horizontal gene transfer event of corresponding genes between gut bacteria and non-vertebrate animals was discovered in the reconstructed phylogenetic tree obtained with maximum likelihood and neighbor-joining methods, respectively. This detailed phylogeny based on totally 128 complete sequences comprised diverse serum amyloid A isoforms from various animal vertebrate and non-vertebrate phyla and also corresponding genes coding for highly similar proteins from animal gut bacteria. Typical largely conserved sequence motifs and a peculiar structural fold consisting mainly of four α-helices in a bundle within all reconstructed clades of the SAA protein family are discussed with respect to their supposed biological functions in various organisms that contain corresponding genes.

TransIntegrator: capture nearly full protein-coding transcript variants via integrating Illumina and PacBio transcriptomes

Genes have the ability to produce transcript variants that perform specific cellular functions. However, accurately detecting all transcript variants remains a long-standing challenge, especially when working with poorly annotated genomes or without a known genome. To address this issue, we have developed a new computational method, TransIntegrator, which enables transcriptome-wide detection of novel transcript variants. For this, we determined 10 Illumina sequencing transcriptomes and a PacBio full-length transcriptome for consecutive embryo development stages of amphioxus, a species of great evolutionary importance. Based on the transcriptomes, we employed TransIntegrator to create a comprehensive transcript variant library, namely iTranscriptome. The resulting iTrancriptome contained 91 915 distinct transcript variants, with an average of 2.4 variants per gene. This substantially improved current amphioxus genome annotation by expanding the number of genes from 21 954 to 38 777. Further analysis manifested that the gene expansion was largely ascribed to integration of multiple Illumina datasets instead of involving the PacBio data. Moreover, we demonstrated an example application of TransIntegrator, via generating iTrancriptome, in aiding accurate transcriptome assembly, which significantly outperformed other hybrid methods such as IDP-denovo and Trinity. For user convenience, we have deposited the source codes of TransIntegrator on GitHub as well as a conda package in Anaconda. In summary, this study proposes an affordable but efficient method for reliable transcriptomic research in most species.

Highly distinct genetic programs for peripheral nervous system formation in chordates

Background

Vertebrates develop their peripheral nervous system (PNS) from transient unique embryonic structures, the neural crest, and the ectodermal placodes that are located at the border of the forming central nervous system. By contrast, in the invertebrate chordates, amphioxus and ascidians, a large part of the PNS originates at the opposite of the embryo, in the ventral ectoderm. In both groups, a biphasic mechanism regulates ventral PNS formation: high BMP levels specify a neurogenic territory within which glutamatergic epidermal sensory neuron formation is controlled by the Notch pathway. Given these similarities and the phylogenetic relationships within chordates, it is likely that ventral PNS is an ancestral feature in chordates and that it has been lost in vertebrates.

Results

In order to get insights into the molecular control of ventral PNS formation and to test the hypothesis of their homology and potential contribution to the emergence of vertebrate PNS, we undertook a close comparison of ventral PNS formation in the ascidian Phallusia mammillata and the amphioxus Branchiostoma lanceolatum. Using timed RNA-seq series, we identified novel markers of the ventral PNS during different phases of its development in both species. By extensively determining the expression of paralogous and orthologous genes, we observed that only a minority of genes have a shared expression in the ventral PNS. However, a large fraction of ventral PNS orthologous genes are expressed in the dorsally forming PNS of vertebrates.

Conclusions

Our work has significantly increased the molecular characterization of ventral PNS formation in invertebrate chordates. The low observed conservation of gene expression in the ventral PNS suggests that the amphioxus and ascidian ventral PNS are either not homologous, or alternatively extensive drift has occurred in their regulatory mechanisms following a long period (600 My) of separate evolution and accelerated evolution in the ascidian lineage. The homology to genes expressed in the dorsally forming PNS of vertebrates suggests that ancestral sensory neurons gene networks have been redeployed in vertebrates.

Crosstalk between nitric oxide and retinoic acid pathways is essential for amphioxus pharynx development

During animal ontogenesis, body axis patterning is finely regulated by complex interactions among several signaling pathways. Nitric oxide (NO) and retinoic acid (RA) are potent morphogens that play a pivotal role in vertebrate development. Their involvement in axial patterning of the head and pharynx shows conserved features in the chordate phylum. Indeed, in the cephalochordate amphioxus, NO and RA are crucial for the correct development of pharyngeal structures. Here, we demonstrate the functional cooperation between NO and RA that occurs during amphioxus embryogenesis. During neurulation, NO modulates RA production through the transcriptional regulation of Aldh1a.2 that irreversibly converts retinaldehyde into RA. On the other hand, RA directly or indirectly regulates the transcription of Nos genes. This reciprocal regulation of NO and RA pathways is essential for the normal pharyngeal development in amphioxus and it could be conserved in vertebrates.

The Spindle Assembly Checkpoint Functions during Early Development in Non-Chordate Embryos

In eukaryotic cells, a spindle assembly checkpoint (SAC) ensures accurate chromosome segregation, by monitoring proper attachment of chromosomes to spindle microtubules and delaying mitotic progression if connections are erroneous or absent. The SAC is thought to be relaxed during early embryonic development. Here, we evaluate the checkpoint response to lack of kinetochore-spindle microtubule interactions in early embryos of diverse animal species. Our analysis shows that there are two classes of embryos, either proficient or deficient for SAC activation during cleavage. Sea urchins, mussels, and jellyfish embryos show a prolonged delay in mitotic progression in the absence of spindle microtubules from the first cleavage division, while ascidian and amphioxus embryos, like those of Xenopus and zebrafish, continue mitotic cycling without delay. SAC competence during early development shows no correlation with cell size, chromosome number, or kinetochore to cell volume ratio. We show that SAC proteins Mad1, Mad2, and Mps1 lack the ability to recognize unattached kinetochores in ascidian embryos, indicating that SAC signaling is not diluted but rather actively silenced during early chordate development.

Genetic regulation of amphioxus somitogenesis informs the evolution of the vertebrate head mesoderm

The evolution of vertebrates from an ancestral chordate was accompanied by the acquisition of a predatory lifestyle closely associated to the origin of a novel anterior structure, the highly specialized head. While the vertebrate head mesoderm is unsegmented, the paraxial mesoderm of the earliest divergent chordate clade, the cephalochordates (amphioxus), is fully segmented in somites. We have previously shown that fibroblast growth factor signalling controls the formation of the most anterior somites in amphioxus; therefore, unravelling the fibroblast growth factor signalling downstream effectors is of crucial importance to shed light on the evolutionary origin of vertebrate head muscles. By using a comparative RNA sequencing approach and genetic functional analyses, we show that several transcription factors, such as Six1/2, Pax3/7 and Zic, act in combination to ensure the formation of three different somite populations. Interestingly, these proteins are orthologous to key regulators of trunk, and not head, muscle formation in vertebrates. Contrary to prevailing thinking, our results suggest that the vertebrate head mesoderm is of visceral and not paraxial origin and support a multistep evolutionary scenario for the appearance of the unsegmented mesoderm of the vertebrates new 'head'.

Nature-inspired engineering of an F-type lectin for increased binding strength

Individual lectin-carbohydrate interactions are usually of low affinity. However, high avidity is frequently attained by the multivalent presentation of glycans on biological surfaces coupled with the occurrence of high order lectin oligomers or tandem repeats of lectin domains in the polypeptide. F-type lectins are l-fucose binding lectins with a typical sequence motif, HX(26)RXDX(4)R/K, whose residues participate in l-fucose binding. We previously reported the presence of a few eukaryotic F-type lectin domains with partial sequence duplication that results in the presence of two l-fucose-binding sequence motifs. We hypothesized that such partial sequence duplication would result in greater avidity of lectin-ligand interactions. Inspired by this example from Nature, we attempted to engineer a bacterial F-type lectin domain from Streptosporangium roseum to attain avid binding by mimicking partial duplication. The engineered lectin demonstrated 12-fold greater binding strength than the wild-type lectin to multivalent fucosylated glycoconjugates. However, the affinity to the monosaccharide l-fucose in solution was similar and partial sequence duplication did not result in an additional functional l-fucose binding site. We also cloned, expressed and purified a Branchiostoma floridae F-type lectin domain with naturally occurring partial sequence duplication and confirmed that the duplicated region with the F-type lectin sequence motif did not participate in l-fucose binding. We found that the greater binding strength of the engineered lectin from S. roseum was instead due to increased oligomerization. We believe that this Nature-inspired strategy might be useful for engineering lectins to improve binding strength in various applications.

A Phylogenomic Framework and Divergence History of Cephalochordata Amphioxus

Amphioxus, or cephalochordates, are often used as the living invertebrate proxy of vertebrate ancestors and are widely used as evolutionary biology models of chordates. However, their phylogeny, divergence history, and speciation characteristics remain poorly understood, and phylogenomic studies to explore these problems lacking entirely from the literature. Here, we determined a new transcriptome of Branchiostoma japonicum. Combined with mass sequences of all other 18 species, a 19-way phylogeny was constructed via multiple methods (ML, BI, PhyloBayes, and ASTRAL), consistently supporting a phylogeny of [(B. belcheri + B. japonicum) + (B. lanceolatum + B. floridae) + Asymmetron lucayanum] in amphioxus. Congruent phylogenetic signals were found across mitochondrial genes, 12S RNA, and complete mitochondrial genomes according to previous reports, indicating that 12S RNA may have potential as a molecular marker for phylogenetic analysis in amphioxus. Molecular dating analysis indicated a radiation of the cephalochordates during the Cretaceous (∼104-61 million years ago), supporting an association between the diversification and speciation of cephalochordates with continental drift and associated changes in their respective habitats during this time. The identified functional enrichment analysis for species-specific domains indicated that their function mainly involves immune response, apoptosis, and lipid metabolism and utilization, signaling that pathogens and changes of energy requirements are an important driving force for amphioxus speciation. This study represents the first large-scale phylogenomic analysis of most major amphioxus genera based on phylogenomic data, providing a new perspective on both phylogeny and divergence speciation of cephalochordates.

Morphological Stasis and Proteome Innovation in Cephalochordates

Lancelets, extant representatives of basal chordates, are prototypic examples of evolutionary stasis; they preserved a morphology and body-plan most similar to the fossil chordates from the early Cambrian. Such a low level of morphological evolution is in harmony with a low rate of amino acid substitution; cephalochordate proteins were shown to evolve slower than those of the slowest evolving vertebrate, the elephant shark. Surprisingly, a study comparing the predicted proteomes of Chinese amphioxus, Branchiostoma belcheri and the Florida amphioxus, Branchiostoma floridae has led to the conclusion that the rate of creation of novel domain combinations is orders of magnitude greater in lancelets than in any other Metazoa, a finding that contradicts the notion that high rates of protein innovation are usually associated with major evolutionary innovations. Our earlier studies on a representative sample of proteins have provided evidence suggesting that the differences in the domain architectures of predicted proteins of these two lancelet species reflect annotation errors, rather than true innovations. In the present work, we have extended these studies to include a larger sample of genes and two additional lancelet species, Asymmetron lucayanum and Branchiostoma lanceolatum. These analyses have confirmed that the domain architecture differences of orthologous proteins of the four lancelet species are because of errors of gene prediction, the error rate in the given species being inversely related to the quality of the transcriptome dataset that was used to aid gene prediction.

Retinoic acid signaling and neurogenic niche regulation in the developing peripheral nervous system of the cephalochordate amphioxus

The retinoic acid (RA) signaling pathway regulates axial patterning and neurogenesis in the developing central nervous system (CNS) of chordates, but little is known about its roles during peripheral nervous system (PNS) formation and about how these roles might have evolved. This study assesses the requirement of RA signaling for establishing a functional PNS in the cephalochordate amphioxus, the best available stand-in for the ancestral chordate condition. Pharmacological manipulation of RA signaling levels during embryogenesis reduces the ability of amphioxus larvae to respond to sensory stimulation and alters the number and distribution of ectodermal sensory neurons (ESNs) in a stage- and context-dependent manner. Using gene expression assays combined with immunohistochemistry, we show that this is because RA signaling specifically acts on a small population of soxb1c-expressing ESN progenitors, which form a neurogenic niche in the trunk ectoderm, to modulate ESN production during elongation of the larval body. Our findings reveal an important role for RA signaling in regulating neurogenic niche activity in the larval amphioxus PNS. Although only few studies have addressed this issue so far, comparable RA signaling functions have been reported for neurogenic niches in the CNS and in certain neurogenic placode derivatives of vertebrates. Accordingly, the here-described mechanism is likely a conserved feature of chordate embryonic and adult neural development.

Pax3/7 duplicated and diverged independently in amphioxus, the basal chordate lineage

The Pax3/7 transcription factor family is integral to developmental gene networks contributing to important innovations in vertebrate evolution, including the neural crest. The basal chordate lineage of amphioxus is ideally placed to understand the dynamics of the gene regulatory network evolution that produced these novelties. We report here the discovery that the cephalochordate lineage possesses two Pax3/7 genes, Pax3/7a and Pax3/7b. The tandem duplication is ancestral to all extant amphioxus, occurring in both Asymmetron and Branchiostoma, but originated after the split from the lineage leading to vertebrates. The two paralogues are differentially expressed during embryonic development, particularly in neural and somitic tissues, suggesting distinct regulation. Our results have implications for the study of amphioxus regeneration, neural plate and crest evolution, and differential tandem paralogue evolution.

Roles of Retinoic Acid Signaling in Shaping the Neuronal Architecture of the Developing Amphioxus Nervous System

The morphogen retinoic acid (RA) patterns vertebrate nervous systems and drives neurogenesis, but how these functions evolved remains elusive. Here, we show that RA signaling plays stage- and tissue-specific roles during the formation of neural cell populations with serotonin, dopamine, and GABA neurotransmitter phenotypes in amphioxus, a proxy for the ancestral chordate. Our data suggest that RA signaling restricts the specification of dopamine-containing cells in the ectoderm and of GABA neurons in the neural tube, probably by regulating Hox1 and Hox3 gene expression, respectively. The two Hox genes thus appear to serve distinct functions rather than to participate in a combinatorial Hox code. We were further able to correlate the RA signaling-dependent mispatterning of hindbrain GABA neurons with concomitant motor impairments. Taken together, these data provide new insights into how RA signaling and Hox genes contribute to nervous system as well as to motor control development in amphioxus and hence shed light on the evolution of these functions within vertebrates.

Nodal/Activin Pathway is a Conserved Neural Induction Signal in Chordates

Neural induction is the process through which pluripotent cells are committed to a neural fate. This first step of Central Nervous System formation is triggered by the "Spemann organizer" in amphibians and by homologous embryonic regions in other vertebrates. Studies in classical vertebrate models have produced contrasting views about the molecular nature of neural inducers and no unifying scheme could be drawn. Moreover, how this process evolved in the chordate lineage remains an unresolved issue. In this work, by using graft and micromanipulation experiments, we definitively establish that the dorsal blastopore lip of the cephalochordate amphioxus is homologous to the vertebrate organizer and is able to trigger the formation of neural tissues in a host embryo. In addition, we demonstrate that Nodal/Activin is the main signal eliciting neural induction in amphioxus, and that it also functions as a bona fide neural inducer in the classical vertebrate model Xenopus. Altogether, our results allow us to propose that Nodal/Activin was a major player of neural induction in the ancestor of chordates. This study further reveals the diversity of neural inducers deployed during chordate evolution and advocates against a universally conserved molecular explanation for this process.

Evolutionary history of the extant amphioxus lineage with shallow-branching diversification

Amphioxus or lancelets have been regarded as a key animal in understanding the origin of vertebrates. However, the evolutionary history within this lineage remains unexplored. As the amphioxus lineage has likely been separated from other chordates for a very long time and displays a marked left-right asymmetry, its evolutionary history is potentially helpful in better understanding chordate and vertebrate origins. We studied the phylogenetic relationships within the extant amphioxus lineage based on mitochondrial genomes incorporating new Asymmetron and Epigonichthys populations, and based on previously reported nuclear transcriptomes. The resulting tree patterns are consistent, showing the Asymmetron clade diverging first, followed by the Epigonichthys and Branchiostoma clades splitting. Divergence time estimates based on nuclear transcriptomes with vertebrate calibrations support a shallow diversification of the extant amphioxus lineage in the Tertiary. These estimates fit well with the closure of seaways between oceans by continental drift, ocean currents, and present geographical distributions, and suggest a long cryptic history from the origin of amphioxus to its most recent diversification. Deduced character polarities based on phylogenetic analyses suggest that the common ancestor of the extant amphioxus existed in a tiny epibenthic state with larva-like appearance of extant amphioxus, likely with ciliate epidermis.

Amphioxus Sp5 is a member of a conserved Specificity Protein complement and is modulated by Wnt/β-catenin signalling

A cluster of three Specificity Protein (Sp) genes (Sp1-4, Sp5 and Sp6-9) is thought to be ancestral in both chordates and the wider Eumetazoa. Sp5 and Sp6-9 gene groups are associated with embryonic growth zones, such as tailbuds, and are both Wnt/β-catenin signalling pathway members and targets. Currently, there are conflicting reports as to the number and identity of Sp genes in the cephalochordates, the sister group to the vertebrates and urochordates. We confirm the SP complement of Branchiostoma belcheri and Branchiostoma lanceolatum, as well as their genomic arrangement, protein domain structure and residue frequency. We assay Sp5 expression in B. lanceolatum embryos, and determine its response to pharmacologically increased β-catenin signalling. Branchiostoma possesses three Sp genes, located on the same genomic scaffold. Phylogenetic and domain structure analyses are consistent with their identification as SP1-4, SP5 and SP6-9, although SP1-4 contains a novel glutamine-rich N-terminal region. SP5 is expressed in axial mesoderm and neurectoderm, and marks the cerebral vesicle and presumptive pharynx. Early exposure to increased β-catenin caused ubiquitous SP5 expression in late gastrula, while later treatment at gastrula stages reduced SP5 expression in the posterior growth zone during axis elongation. Amphioxus possess a typical invertebrate eumetazoan SP complement, and SP5 expression in embryos is well conserved with vertebrate homologues. Its expression in the tailbud, a posterior growth zone, is consistent with expression seen in other bilaterians. Branchiostoma SP5 shows a dynamic response to Wnt/β-catenin signalling.

Prdm12 specifies V1 interneurons through cross-repressive interactions with Dbx1 and Nkx6 genes in Xenopus

V1 interneurons are inhibitory neurons that play an essential role in vertebrate locomotion. The molecular mechanisms underlying their genesis remain, however, largely undefined. Here, we show that the transcription factor Prdm12 is selectively expressed in p1 progenitors of the hindbrain and spinal cord in the frog embryo, and that a similar restricted expression profile is observed in the nerve cord of other vertebrates as well as of the cephalochordate amphioxus. Using frog, chick and mice, we analyzed the regulation of Prdm12 and found that its expression in the caudal neural tube is dependent on retinoic acid and Pax6, and that it is restricted to p1 progenitors, due to the repressive action of Dbx1 and Nkx6-1/2 expressed in the adjacent p0 and p2 domains. Functional studies in the frog, including genome-wide identification of its targets by RNA-seq and ChIP-Seq, reveal that vertebrate Prdm12 proteins act as a general determinant of V1 cell fate, at least in part, by directly repressing Dbx1 and Nkx6 genes. This probably occurs by recruiting the methyltransferase G9a, an activity that is not displayed by the amphioxus Prdm12 protein. Together, these findings indicate that Prdm12 promotes V1 interneurons through cross-repressive interactions with Dbx1 and Nkx6 genes, and suggest that this function might have only been acquired after the split of the vertebrate and cephalochordate lineages.

Homeodomain proteins: an update

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.

Identification, evolution and expression of an insulin-like peptide in the cephalochordate Branchiostoma lanceolatum

Insulin is one of the most studied proteins since it is central to the regulation of carbohydrate and fat metabolism in vertebrates and its expression and release are disturbed in diabetes, the most frequent human metabolic disease worldwide. However, the evolution of the function of the insulin protein family is still unclear. In this study, we present a phylogenetic and developmental analysis of the Insulin Like Peptide (ILP) in the cephalochordate amphioxus. We identified an ILP in the European amphioxus Branchiostoma lanceolatum that displays structural characteristics of both vertebrate insulin and Insulin-like Growth Factors (IGFs). Our phylogenetic analysis revealed that amphioxus ILP represents the sister group of both vertebrate insulin and IGF proteins. We also characterized both temporal and spatial expression of ILP in amphioxus. We show that ilp is highly expressed in endoderm and paraxial mesoderm during development, and mainly expressed in the gut of both the developing embryo and adult. We hypothesize that ILP has critical implications in both developmental processes and metabolism and could display IGF- and insulin-like functions in amphioxus supporting the idea of a common ancestral protein.

The Nodal signaling pathway controls left-right asymmetric development in amphioxus

Background

Nodal is an important determinant of the left-right (LR) body axis in bilaterians, specifying the right side in protostomes and non-chordate deuterostomes as opposed to the left side in chordates. Amphioxus represents an early-branching chordate group, rendering it especially useful for studying the character states that predate the origin of vertebrates. However, its anatomy, involving offset arrangement of axial structures, marked asymmetry of the oropharyngeal region, and, most notably, a mouth positioned on the left side, contrasts with the symmetric arrangement of the corresponding regions in other chordates.

Results

We show that the Nodal signaling pathway acts to specify the LR axis in the cephalochordate amphioxus in a similar way as in vertebrates. At early neurula stages, Nodal switches from initial bilateral to the left-sided expression and subsequently specifies the left embryonic side. Perturbation of Nodal signaling with small chemical inhibitors (SB505124 and SB431542) alters expression of other members of the pathway and of left/right-sided, organ-specific genes. Upon inhibition, larvae display loss of the innate alternation of both somites and axons of peripheral nerves and loss of left-sided pharyngeal structures, such as the mouth, the preoral pit, and the duct of the club-shaped gland. Concomitantly, the left side displays ectopic expression of otherwise right-sided genes, and the larvae exhibit bilaterally symmetrical morphology, with duplicated endostyle and club-shaped gland structures.

Conclusions

We demonstrate that Nodal signaling is necessary for establishing the LR embryonic axis and for developing profound asymmetry in amphioxus. Our data suggest that initial symmetry breaking in amphioxus and propagation of the pathway on the left side correspond with the situation in vertebrates. However, the organs that become targets of the pathway differ between amphioxus and vertebrates, which may explain the pronounced asymmetry of its oropharyngeal and axial structures and the left-sided position of the mouth.

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-6-5) contains supplementary material, which is available to authorized users.

Genome-wide DNA methylation profiles and their relationships with mRNA and the microRNA transcriptome in bovine muscle tissue (Bos taurine)

DNA methylation is a key epigenetic modification in mammals and plays important roles in muscle development. We sampled longissimus dorsi muscle (LDM) from a well-known elite native breed of Chinese Qinchuan cattle living within the same environment but displaying distinct skeletal muscle at the fetal and adult stages. We generated and provided a genome-wide landscape of DNA methylomes and their relationship with mRNA and miRNA for fetal and adult muscle studies. Integration analysis revealed a total of 77 and 1,054 negatively correlated genes with methylation in the promoter and gene body regions, respectively, in both the fetal and adult bovine libraries. Furthermore, we identified expression patterns of high-read genes that exhibit a negative correlation between methylation and expression from nine different tissues at multiple developmental stages of bovine muscle-related tissue or organs. In addition, we validated the MeDIP-Seq results by bisulfite sequencing PCR (BSP) in some of the differentially methylated promoters. Together, these results provide valuable data for future biomedical research and genomic and epigenomic studies of bovine skeletal muscle that may help uncover the molecular basis underlying economically valuable traits in cattle. This comprehensive map also provides a solid basis for exploring the epigenetic mechanisms of muscle growth and development.

SeaBase: a multispecies transcriptomic resource and platform for gene network inference

Marine and aquatic animals are extraordinarily useful as models for identifying mechanisms of development and evolution, regeneration, resistance to cancer, longevity and symbiosis, among many other areas of research. This is due to the great diversity of these organisms and their wide-ranging capabilities. Genomics tools are essential for taking advantage of these "free lessons" of nature. However, genomics and transcriptomics are challenging in emerging model systems. Here, we present SeaBase, a tool for helping to meet these needs. Specifically, SeaBase provides a platform for sharing and searching transcriptome data. More importantly, SeaBase will support a growing number of tools for inferring gene network mechanisms. The first dataset available on SeaBase is a developmental transcriptomic profile of the sea anemone Nematostella vectensis (Anthozoa, Cnidaria). Additional datasets are currently being prepared and we are aiming to expand SeaBase to include user-supplied data for any number of marine and aquatic organisms, thereby supporting many potentially new models for gene network studies. SeaBase can be accessed online at: http://seabase.core.cli.mbl.edu.

Is the distribution of the lancelet Branchiostoma caribaeum affected by sewage discharges? An analysis at multiple scales of variability

Spatial variation in the density and biomass of Branchiostoma caribaeum was analyzed along a sewage contamination gradient identified by fecal steroids in a subtropical estuary, southern Brazil. Sampling, repeated in the austral winter and summer, followed a hierarchical design nested at four spatial scales (sector>1 km; area>100 m; site>10 m; replicate<1 m). Density and biomass were significantly lower at sites characterized by high concentrations of fecal steroids. The best combinations of variables that explained the biological similarities among sites involved contamination indicators. Most of the variation of biological data was found at the smallest scales and could be related with the sediment texture. Our study highlighted the usefulness of a multi-scale perspective to evaluate distribution patterns of benthic invertebrates as a biological indication of environmental pollution. Gradient analyses at larger spatial scales may be invalidated by the patchy distribution of benthic fauna if they do not account for such small scale variability.

Broken colinearity of the amphioxus Hox cluster

Background

In most eumetazoans studied so far, Hox genes determine the identity of structures along the main body axis. They are usually linked in genomic clusters and, in the case of the vertebrate embryo, are expressed with spatial and temporal colinearity. Outside vertebrates, temporal colinearity has been reported in the cephalochordate amphioxus (the least derived living relative of the chordate ancestor) but only for anterior and central genes, namely Hox1 to Hox4 and Hox6. However, most of the Hox gene expression patterns in amphioxus have not been reported. To gain global insights into the evolution of Hox clusters in chordates, we investigated a more extended expression profile of amphioxus Hox genes.

Results

Here we report an extended expression profile of the European amphioxus Branchiostoma lanceolatum Hox genes and describe that all Hox genes, except Hox13, are expressed during development. Interestingly, we report the breaking of both spatial and temporal colinearity for at least Hox6 and Hox14, which thus have escaped from the classical Hox code concept. We show a previously unidentified Hox6 expression pattern and a faint expression for posterior Hox genes in structures such as the posterior mesoderm, notochord, and hindgut. Unexpectedly, we found that amphioxus Hox14 had the most divergent expression pattern. This gene is expressed in the anterior cerebral vesicle and pharyngeal endoderm. Amphioxus Hox14 expression represents the first report of Hox gene expression in the most anterior part of the central nervous system. Nevertheless, despite these divergent expression patterns, amphioxus Hox6 and Hox14 seem to be still regulated by retinoic acid.

Conclusions

Escape from colinearity by Hox genes is not unusual in either vertebrates or amphioxus and we suggest that those genes escaping from it are probably associated with the patterning of lineage-specific morphological traits, requiring the loss of those developmental constraints that kept them colinear.

Metal dealing at the origin of the Chordata phylum: the metallothionein system and metal overload response in amphioxus

Non-vertebrate chordates, specifically amphioxus, are considered of the utmost interest for gaining insight into the evolutionary trends, i.e. differentiation and specialization, of gene/protein systems. In this work, MTs (metallothioneins), the most important metal binding proteins, are characterized for the first time in the cephalochordate subphylum at both gene and protein level, together with the main features defining the amphioxus response to cadmium and copper overload. Two MT genes (BfMT1 and BfMT2) have been identified in a contiguous region of the genome, as well as several ARE (antioxidant response element) and MRE (metal response element) located upstream the transcribed region. Their corresponding cDNAs exhibit identical sequence in the two lancelet species (B. floridae and B. lanceolatum), BfMT2 cDNA resulting from an alternative splicing event. BfMT1 is a polyvalent metal binding peptide that coordinates any of the studied metal ions (Zn, Cd or Cu) rendering complexes stable enough to last in physiological environments, which is fully concordant with the constitutive expression of its gene, and therefore, with a metal homeostasis housekeeping role. On the contrary, BfMT2 exhibits a clear ability to coordinate Cd(II) ions, while it is absolutely unable to fold into stable Cu (I) complexes, even as mixed species. This identifies it as an essential detoxification agent, which is consequently only induced in emergency situations. The cephalochordate MTs are not directly related to vertebrate MTs, neither by gene structure, protein similarity nor metal-binding behavior of the encoded peptides. The closest relative is the echinoderm MT, which confirm proposed phylogenetic relationships between these two groups. The current findings support the existence in most organisms of two types of MTs as for their metal binding preferences, devoted to different biological functions: multivalent MTs for housekeeping roles, and specialized MTs that evolve either as Cd-thioneins or Cu-thioneins, according to the ecophysiological needs of each kind of organisms.