The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins

The origin and evolution of human glutaminases and their atypical C-terminal ankyrin repeats

On the basis of tissue-specific enzyme activity and inhibition by catalytic products, Hans Krebs first demonstrated the existence of multiple glutaminases in mammals. Currently, two human genes are known to encode at least four glutaminase isoforms. However, the phylogeny of these medically relevant enzymes remains unclear, prompting us to investigate their origin and evolution. Using prokaryotic and eukaryotic glutaminase sequences, we built a phylogenetic tree whose topology suggested that the multidomain architecture was inherited from bacterial ancestors, probably simultaneously with the hosting of the proto-mitochondrion endosymbiont. We propose an evolutionary model wherein the appearance of the most active enzyme isoform, glutaminase C (GAC), which is expressed in many cancers, was a late retrotransposition event that occurred in fishes from the Chondrichthyes class. The ankyrin (ANK) repeats in the glutaminases were acquired early in their evolution. To obtain information on ANK folding, we solved two high-resolution structures of the ANK repeat-containing C termini of both kidney-type glutaminase (KGA) and GLS2 isoforms (glutaminase B and liver-type glutaminase). We found that the glutaminase ANK repeats form unique intramolecular contacts through two highly conserved motifs; curiously, this arrangement occludes a region usually involved in ANK-mediated protein-protein interactions. We also solved the crystal structure of full-length KGA and present a small-angle X-ray scattering model for full-length GLS2. These structures explain these proteins' compromised ability to assemble into catalytically active supra-tetrameric filaments, as previously shown for GAC. Collectively, these results provide information about glutaminases that may aid in the design of isoform-specific glutaminase inhibitors.

Diversity and evolution of TIR-domain-containing proteins in bivalves and Metazoa: New insights from comparative genomics

The Toll/interleukin-1 receptor (TIR) domain has a fundamental role in the innate defence response of plants, vertebrate and invertebrate animals. Mostly found in the cytosolic side of membrane-bound receptor proteins, it mediates the intracellular signalling upon pathogen recognition via heterotypic interactions. Although a number of TIR-domain-containing (TIR-DC) proteins have been characterized in vertebrates, their evolutionary relationships and functional role in protostomes are still largely unknown. Due to the high abundance and diversity of TIR-DC proteins in bivalve molluscs, we investigated this class of marine invertebrates as a case study. The analysis of the available genomic and transcriptomic data allowed the identification of over 400 full-length sequences and their classification in protein families based on sequence homology and domain organization. In addition to TLRs and MyD88 adaptors, bivalves possess a surprisingly large repertoire of intracellular TIR-DC proteins, which are conserved across a broad range of metazoan taxa. Overall, we report the expansion and diversification of TIR-DC proteins in several invertebrate lineages and the identification of many novel protein families possibly involved in both immune-related signalling and embryonic development.

A microRNA-mRNA expression network during oral siphon regeneration in Ciona

Here we present a parallel study of mRNA and microRNA expression during oral siphon (OS) regeneration in Ciona robusta, and the derived network of their interactions. In the process of identifying 248 mRNAs and 15 microRNAs as differentially expressed, we also identified 57 novel microRNAs, several of which are among the most highly differentially expressed. Analysis of functional categories identified enriched transcripts related to stress responses and apoptosis at the wound healing stage, signaling pathways including Wnt and TGFβ during early regrowth, and negative regulation of extracellular proteases in late stage regeneration. Consistent with the expression results, we found that inhibition of TGFβ signaling blocked OS regeneration. A correlation network was subsequently inferred for all predicted microRNA-mRNA target pairs expressed during regeneration. Network-based clustering associated transcripts into 22 non-overlapping groups, the functional analysis of which showed enrichment of stress response, signaling pathway and extracellular protease categories that could be related to specific microRNAs. Predicted targets of the miR-9 cluster suggest a role in regulating differentiation and the proliferative state of neural progenitors through regulation of the cytoskeleton and cell cycle.

Genome-wide survey of miRNAs and their evolutionary history in the ascidian, Halocynthia roretzi

BACKGROUND

miRNAs play essential roles in the modulation of cellular functions via degradation and/or translation attenuation of target mRNAs. They have been surveyed in a single ascidian genus, Ciona. Recently, an annotated draft genome sequence for a distantly related ascidian, Halocynthia roretzi, has become available, but miRNAs in H. roretzi have not been previously studied.

RESULTS

We report the prediction of 319 candidate H. roretzi miRNAs, obtained through three complementary methods. Experimental validation suggests that more than half of these candidate miRNAs are expressed during embryogenesis. The majority of predicted H. roretzi miRNAs appear specific to ascidians or tunicates, and only 32 candidates, belonging to 25 families, are widely conserved across metazoans.

CONCLUSION

Our study presents a comprehensive identification of candidate H. roretzi miRNAs. This resource will facilitate the study of the mechanisms for miRNA-controlled gene regulatory networks during ascidian development. Further, our analysis suggests that the majority of Halocynthia miRNAs are specific to ascidian or tunicates, with only a small number of widely conserved miRNAs. This result is consistent with the general notion that animal miRNAs are less conserved between taxa than plant ones.

Molecular regionalization of the developing amphioxus neural tube challenges major partitions of the vertebrate brain

All vertebrate brains develop following a common Bauplan defined by anteroposterior (AP) and dorsoventral (DV) subdivisions, characterized by largely conserved differential expression of gene markers. However, it is still unclear how this Bauplan originated during evolution. We studied the relative expression of 48 genes with key roles in vertebrate neural patterning in a representative amphioxus embryonic stage. Unlike nonchordates, amphioxus develops its central nervous system (CNS) from a neural plate that is homologous to that of vertebrates, allowing direct topological comparisons. The resulting genoarchitectonic model revealed that the amphioxus incipient neural tube is unexpectedly complex, consisting of several AP and DV molecular partitions. Strikingly, comparison with vertebrates indicates that the vertebrate thalamus, pretectum, and midbrain domains jointly correspond to a single amphioxus region, which we termed Di-Mesencephalic primordium (DiMes). This suggests that these domains have a common developmental and evolutionary origin, as supported by functional experiments manipulating secondary organizers in zebrafish and mice.

The chordate ancestor possessed a single copy of the Brachyury gene for notochord acquisition

BACKGROUND

The T-box family transcription-factor gene, Brachyury, has two expression domains with discrete functions during animal embryogenesis. The primary domain, associated with the blastopore, is shared by most metazoans, while the secondary domain, involved in the notochord, is specific to chordates. In most animals, Brachyury is present in a single copy, but in cephalochordates, the most basal of the chordates, the gene is present in two copies, suggesting allotment of the two domains to each of the duplicates.

RESULTS

In order to clarify whether Brachyury duplication occurred in the common ancestor of chordates after which one of duplicates was lost in the urochordate and vertebrate lineages, we estimated phylogenetic relationships of Brachyury genes and examined the synteny of a Brachyury-containing genomic region of deuterostomes with decoded genomes. The monophyletic origin of tandemly arranged Brachyury genes of cephalochordates indicates that the tandem duplication occurred in the cephalochordate lineage, but not in the chordate ancestor.

CONCLUSIONS

Our results thus suggest that, in the common ancestor of chordates, a single copy of Brachyury acquired two expression domains and that the duplication was not involved in the acquisition of the notochord. However, in relation to regulatory mechanisms, both possibilities-namely a single copy with two domains and two copies with different domains-should be considered in future studies of Brachyury.

metaSPAdes: a new versatile metagenomic assembler

While metagenomics has emerged as a technology of choice for analyzing bacterial populations, the assembly of metagenomic data remains challenging, thus stifling biological discoveries. Moreover, recent studies revealed that complex bacterial populations may be composed from dozens of related strains, thus further amplifying the challenge of metagenomic assembly. metaSPAdes addresses various challenges of metagenomic assembly by capitalizing on computational ideas that proved to be useful in assemblies of single cells and highly polymorphic diploid genomes. We benchmark metaSPAdes against other state-of-the-art metagenome assemblers and demonstrate that it results in high-quality assemblies across diverse data sets.

Conservation of the Keap1-Nrf2 System: An Evolutionary Journey through Stressful Space and Time

The Keap1-Nrf2 system is an evolutionarily conserved defense mechanism against oxidative and xenobiotic stress. Its regulatory mechanisms, e.g., stress-sensing mechanism, proteasome-based regulation of Nrf2 activity and selection of target genes, have been elucidated mainly in mammals. In addition, emerging model animals, such as zebrafish, fruit fly and Caenorhabditis elegans, have been shown to have similar anti-stress systems to mammals, suggesting that analogous defense systems are widely conserved throughout the animal kingdom. Experimental evidence in lower animals provides important information beyond mere laboratory-confined utility, such as regarding how these systems transformed during evolution, which may help characterize the mammalian system in greater detail. Recent advances in genome projects of both model and non-model animals have provided a great deal of useful information toward this end. We herein review the research on Keap1-Nrf2 and its analogous systems in both mammals and lower model animals. In addition, by comparing the amino acid sequences of Nrf2 and Keap1 proteins from various species, we can deduce the evolutionary history of the anti-stress system. This combinatorial approach using both experimental and genetic data will suggest perspectives of approach for researchers studying the stress response.

Genetic signatures of natural selection in a model invasive ascidian

Invasive species represent promising models to study species’ responses to rapidly changing environments. Although local adaptation frequently occurs during contemporary range expansion, the associated genetic signatures at both population and genomic levels remain largely unknown. Here, we use genome-wide gene-associated microsatellites to investigate genetic signatures of natural selection in a model invasive ascidian, Ciona robusta. Population genetic analyses of 150 individuals sampled in Korea, New Zealand, South Africa and Spain showed significant genetic differentiation among populations. Based on outlier tests, we found high incidence of signatures of directional selection at 19 loci. Hitchhiking mapping analyses identified 12 directional selective sweep regions, and all selective sweep windows on chromosomes were narrow (~8.9 kb). Further analyses indentified 132 candidate genes under selection. When we compared our genetic data and six crucial environmental variables, 16 putatively selected loci showed significant correlation with these environmental variables. This suggests that the local environmental conditions have left significant signatures of selection at both population and genomic levels. Finally, we identified “plastic” genomic regions and genes that are promising regions to investigate evolutionary responses to rapid environmental change in C. robusta.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

How complexity increases in development: An analysis of the spatial-temporal dynamics of Gene expression in Ciona intestinalis

The increase in complexity in an embryo over developmental time is perhaps one of the most intuitive processes of animal development. It is also intuitive that the embryo becomes progressively compartmentalized over time and space. In spite of this intuitiveness, there are no systematic attempts to quantify how this occurs. Here, we present a quantitative analysis of the compartmentalization and spatial complexity of Ciona intestinalis over developmental time by analyzing thousands of gene expression spatial patterns from the ANISEED database. We measure compartmentalization in two ways: as the relative volume of expression of genes and as the disparity in gene expression between body parts. We also use a measure of the curvature of each gene expression pattern in 3D space. These measures show a similar increase over time, with the most dramatic change occurring from the 112-cell stage to the early tailbud stage. Combined, these measures point to a global pattern of increase in complexity in the Ciona embryo. Finally, we cluster the different regions of the embryo depending on their gene expression similarity, within and between stages. Results from this clustering analysis, which partially correspond to known fate maps, provide a global quantitative overview about differentiation and compartmentalization between body parts at each developmental stage.

Unsupervised detection of regulatory gene expression information in different genomic regions enables gene expression ranking

Background

The regulation of all gene expression steps (e.g., Transcription, RNA processing, Translation, and mRNA Degradation) is known to be primarily encoded in different parts of genes and in genomic regions in proximity to genes (e.g., promoters, untranslated regions, coding regions, introns, etc.). However, the entire gene expression codes and the genomic regions where they are encoded are still unknown.

Results

Here, we employ an unsupervised approach to estimate the concentration of gene expression codes in different non-coding parts of genes and transcripts, such as introns and untranslated regions, focusing on three model organisms (Escherichia coli, Saccharomyces cerevisiae, and Schizosaccharomyces pombe). Our analyses support the conjecture that regions adjacent to the beginning and end of ORFs and the beginning and end of introns tend to include higher concentration of gene expression information relatively to regions further away. In addition, we report the exact regions with elevated concentration of gene expression codes. Furthermore, we demonstrate that the concentration of these codes in different genetic regions is correlated with the expression levels of the corresponding genes, and with splicing efficiency measurements and meiotic stage gene expression measurements in S. cerevisiae.

Conclusion

We suggest that these discoveries improve our understanding of gene expression regulation and evolution; they can also be used for developing improved models of genome/gene evolution and for engineering gene expression in various biotechnological and synthetic biology applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1497-z) contains supplementary material, which is available to authorized users.

New insights into the evolutionary origins of the recombination-activating gene proteins and V(D)J recombination

The adaptive immune system of jawed vertebrates relies on V(D)J recombination as one of the main processes to generate the diverse array of receptors necessary for the recognition of a wide range of pathogens. The DNA cleavage reaction necessary for the assembly of the antigen receptor genes from an array of potential gene segments is mediated by the recombination-activating gene proteins RAG1 and RAG2. The RAG proteins have been proposed to originate from a transposable element (TE) as they share mechanistic and structural similarities with several families of transposases and are themselves capable of mediating transposition. A number of RAG-like proteins and TEs with sequence similarity to RAG1 and RAG2 have been identified, but only recently has their function begun to be characterized, revealing mechanistic links to the vertebrate RAGs. Of particular significance is the discovery of ProtoRAG, a transposon superfamily found in the genome of the basal chordate amphioxus. ProtoRAG has many of the sequence and mechanistic features predicted for the ancestral RAG transposon and is likely to be an evolutionary relative of RAG1 and RAG2. In addition, early observations suggesting that RAG1 is able to mediate V(D)J recombination in the absence of RAG2 have been confirmed, implying independent evolutionary origins for the two RAG genes. Here, recent progress in identifying and characterizing RAG-like proteins and the TEs that encode them is summarized and a refined model for the evolution of V(D)J recombination and the RAG proteins is presented.

Genome Editing of the Ascidian Ciona intestinalis with TALE Nuclease

The ascidian Ciona intestinalis is an important model animal for studying developmental mechanisms for constructing the chordate body. Although molecular and embryological techniques for manipulating Ciona genes were developed a long time ago, recent achievements of genome editing in this animal have innovated functional analyses of genes in Ciona. Particularly, knockout of genes in the G0 generation coupled with tissue-specific expression of TALENs enables us to rapidly address gene functions that were difficult using previous methods.

Transcriptional regulation of a horizontally transferred gene from bacterium to chordate

The horizontal transfer of genes between distantly related organisms is undoubtedly a major factor in the evolution of novel traits. Because genes are functionless without expression, horizontally transferred genes must acquire appropriate transcriptional regulations in their recipient organisms, although the evolutionary mechanism is not known well. The defining characteristic of tunicates is the presence of a cellulose containing tunic covering the adult and larval body surface. Cellulose synthase was acquired by horizontal gene transfer from Actinobacteria. We found that acquisition of the binding site of AP-2 transcription factor was essential for tunicate cellulose synthase to gain epidermal-specific expression. Actinobacteria have very GC-rich genomes, regions of which are capable of inducing specific expression in the tunicate epidermis as the AP-2 binds to a GC-rich region. Therefore, the actinobacterial cellulose synthase could have been potentiated to evolve its new function in the ancestor of tunicates with a higher probability than the evolution depending solely on a spontaneous event.

Evolutionary and Functional Diversification of the Vitamin D Receptor-Lithocholic Acid Partnership

The evolution, molecular behavior, and physiological function of nuclear receptors are of particular interest given their diverse roles in regulating essential biological processes. The vitamin D receptor (VDR) is well known for its canonical roles in calcium homeostasis and skeletal maintenance. Additionally, VDR has received an increased amount of attention due to the discovery of numerous non-calcemic functions, including the detoxification of lithocholic acid. Lithocholic acid is a toxic metabolite of chenodeoxycholic acid, a primary bile acid. The partnership between the VDR and lithocholic acid has been hypothesized to be a recent adaptation that evolved to mediate the detoxification and elimination of lithocholic acid from the gut. This partnership is speculated to be limited to higher vertebrates (birds and mammals), as lower vertebrates do not synthesize the parent compound of lithocholic acid. However, the molecular functions associated with the observed insensitivity of basal VDRs to lithocholic acid have not been explored. Here we characterize canonical nuclear receptor functions of VDRs from select species representing key nodes in vertebrate evolution and span a range of bile salt phenotypes. Competitive ligand binding assays revealed that the receptor's affinity for lithocholic acid is highly conserved across species, suggesting that lithocholic acid affinity is an ancient and non-adaptive trait. However, transient transactivation assays revealed that lithocholic acid-mediated VDR activation might have evolved more recently, as the non-mammalian receptors did not respond to lithocholic acid unless exogenous coactivator proteins were co-expressed. Subsequent functional assays indicated that differential lithocholic acid-mediated receptor activation is potentially driven by differential protein-protein interactions between VDR and nuclear receptor coregulator proteins. We hypothesize that the vitamin D receptor-lithocholic acid partnership evolved as a by-product of natural selection on the ligand-receptor partnership between the vitamin D receptor and the native VDR ligand: 1α,25-dihydroxyvitamin D3, the biologically active metabolite of vitamin D3.

The CNS connectome of a tadpole larva of Ciona intestinalis (L.) highlights sidedness in the brain of a chordate sibling

Left-right asymmetries in brains are usually minor or cryptic. We report brain asymmetries in the tiny, dorsal tubular nervous system of the ascidian tadpole larva, Ciona intestinalis. Chordate in body plan and development, the larva provides an outstanding example of brain asymmetry. Although early neural development is well studied, detailed cellular organization of the swimming larva's CNS remains unreported. Using serial-section EM we document the synaptic connectome of the larva's 177 CNS neurons. These formed 6618 synapses including 1772 neuromuscular junctions, augmented by 1206 gap junctions. Neurons are unipolar with at most a single dendrite, and few synapses. Some synapses are unpolarised, others form reciprocal or serial motifs; 922 were polyadic. Axo-axonal synapses predominate. Most neurons have ciliary organelles, and many features lack structural specialization. Despite equal cell numbers on both sides, neuron identities and pathways differ left/right. Brain vesicle asymmetries include a right ocellus and left coronet cells.

Localization of Neuropeptide Gene Expression in Larvae of an Echinoderm, the Starfish Asterias rubens

Neuropeptides are an ancient class of neuronal signaling molecules that regulate a variety of physiological and behavioral processes in animals. The life cycle of many animals includes a larval stage(s) that precedes metamorphic transition to a reproductively active adult stage but, with the exception of Drosophila melanogaster and other insects, research on neuropeptide signaling has hitherto largely focused on adult animals. However, recent advances in genome/transcriptome sequencing have facilitated investigation of neuropeptide expression/function in the larvae of protostomian (e.g., the annelid Platynereis dumerilii) and deuterostomian (e.g., the urochordate Ciona intestinalis) invertebrates. Accordingly, here we report the first multi-gene investigation of larval neuropeptide precursor expression in a species belonging to the phylum Echinodermata-the starfish Asterias rubens. Whole-mount mRNA in situ hybridization was used to visualize in bipinnaria and brachiolaria stage larvae the expression of eight neuropeptide precursors: L-type SALMFamide (S1), F-type SALMFamide (S2), vasopressin/oxytocin-type, NGFFYamide, thyrotropin-releasing hormone-type, gonadotropin-releasing hormone-type, calcitonin-type and corticotropin-releasing hormone-type. Expression of only three of the precursors (S1, S2, NGFFYamide) was observed in bipinnaria larvae but by the brachiolaria stage expression of all eight precursors was detected. An evolutionarily conserved feature of larval nervous systems is the apical organ and in starfish larvae this comprises the bilaterally symmetrical lateral ganglia, but only the S1 and S2 precursors were found to be expressed in these ganglia. A prominent feature of brachiolaria larvae is the attachment complex, comprising the brachia and adhesive disk, which mediates larval attachment to a substratum prior to metamorphosis. Interestingly, all of the neuropeptide precursors examined here are expressed in the attachment complex, with distinctive patterns of expression suggesting potential roles for neuropeptides in the attachment process. Lastly, expression of several neuropeptide precursors is associated with ciliary bands, suggesting potential roles for the neuropeptides derived from these precursors in control of larval locomotion and/or feeding. In conclusion, our findings provide novel perspectives on the evolution and development of neuropeptide signaling systems and neuroanatomical insights into neuropeptide function in echinoderm larvae.

Studying the genetic basis of speciation in high gene flow marine invertebrates

A growing number of genes responsible for reproductive incompatibilities between species (barrier loci) exhibit the signals of positive selection. However, the possibility that genes experiencing positive selection diverge early in speciation and commonly cause reproductive incompatibilities has not been systematically investigated on a genome-wide scale. Here, I outline a research program for studying the genetic basis of speciation in broadcast spawning marine invertebrates that uses a priori genome-wide information on a large, unbiased sample of genes tested for positive selection. A targeted sequence capture approach is proposed that scores single-nucleotide polymorphisms (SNPs) in widely separated species populations at an early stage of allopatric divergence. The targeted capture of both coding and non-coding sequences enables SNPs to be characterized at known locations across the genome and at genes with known selective or neutral histories. The neutral coding and non-coding SNPs provide robust background distributions for identifying FST-outliers within genes that can, in principle, identify specific mutations experiencing diversifying selection. If natural hybridization occurs between species, the neutral coding and non-coding SNPs can provide a neutral admixture model for genomic clines analyses aimed at finding genes exhibiting strong blocks to introgression. Strongylocentrotid sea urchins are used as a model system to outline the approach but it can be used for any group that has a complete reference genome available.

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea)

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni. Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

Embryo Microinjection and Electroporation in the Chordate Ciona intestinalis

Simple model organisms are instrumental for in vivo studies of developmental and cellular differentiation processes. Currently, the evolutionary distance to man of conventional invertebrate model systems and the complexity of genomes in vertebrates are critical challenges to modeling human normal and pathological conditions. The chordate Ciona intestinalis is an invertebrate chordate that emerged from a common ancestor with the vertebrates and may represent features at the interface between invertebrates and vertebrates. A common body plan with much simpler cellular and genomic composition should unveil gene regulatory network (GRN) links and functional genomics readouts explaining phenomena in the vertebrate condition. The compact genome of Ciona, a fixed embryonic lineage with few divisions and large cells, combined with versatile community tools foster efficient gene functional analyses in this organism. Here, we present several crucial methods for this promising model organism, which belongs to the closest sister group to vertebrates. We present protocols for transient transgenesis by electroporation, along with microinjection-mediated gene knockdown, which together provide the means to study gene function and genomic regulatory elements. We extend our protocols to provide information on how community databases are utilized for in silico design of gene regulatory or gene functional experiments. An example study demonstrates how novel information can be gained on the interplay, and its quantification, of selected neural factors conserved between Ciona and man. Furthermore, we show examples of differential subcellular localization in embryonic cells, following DNA electroporation in Ciona zygotes. Finally, we discuss the potential of these protocols to be adapted for tissue specific gene interference with emerging gene editing methods. The in vivo approaches in Ciona overcome major shortcomings of classical model organisms in the quest of unraveling conserved mechanisms in the chordate developmental program, relevant to stem cell research, drug discovery, and subsequent clinical application.

Distinguishing contemporary hybridization from past introgression with postgenomic ancestry-informative SNPs in strongly differentiated Ciona species

Biological introductions bring into contact species that can still hybridize. The evolutionary outcomes of such secondary contacts may be diverse (e.g. adaptive introgression from or into the introduced species) but are not yet well examined in the wild. The recent secondary contact between the non-native sea squirt Ciona robusta (formerly known as C. intestinalis type A) and its native congener C. intestinalis (formerly known as C. intestinalis type B), in the Western English Channel, provides an excellent case study to examine. To examine contemporary hybridization between the two species, we developed a panel of 310 ancestry-informative SNPs from a population transcriptomic study. Hybridization rates were examined on 449 individuals sampled in eight sites from the sympatric range and five sites from allopatric ranges. The results clearly showed an almost complete absence of contemporary hybridization between the two species in syntopic localities, with only one-first-generation hybrid and no other genotype compatible with recent backcrosses. Despite the almost lack of contemporary hybridization, shared polymorphisms were observed in sympatric and allopatric populations of both species. Furthermore, one allopatric population from SE Pacific exhibited a higher rate of shared polymorphisms compared to all other C. robusta populations. Altogether, these results indicate that the observed level of shared polymorphism is more probably the outcome of ancient gene flow spread afterwards at a worldwide scale. They also emphasize efficient reproductive barriers preventing hybridization between introduced and native species, which suggests hybridization should not impede too much the expansion and the establishment of the non-native species in its introduction range.

SOXE neofunctionalization and elaboration of the neural crest during chordate evolution

During chordate evolution, two genome-wide duplications facilitated acquisition of vertebrate traits, including emergence of neural crest cells (NCCs), in which neofunctionalization of the duplicated genes are thought to have facilitated development of craniofacial structures and the peripheral nervous system. How these duplicated genes evolve and acquire the ability to specify NC and their derivatives are largely unknown. Vertebrate SoxE paralogues, most notably Sox9/10, are essential for NC induction, delamination and lineage specification. In contrast, the basal chordate, amphioxus, has a single SoxE gene and lacks NC-like cells. Here, we test the hypothesis that duplication and divergence of an ancestral SoxE gene may have facilitated elaboration of NC lineages. By using an in vivo expression assay to compare effects of AmphiSoxE and vertebrate Sox9 on NC development, we demonstrate that all SOXE proteins possess similar DNA binding and homodimerization properties and can induce NCCs. However, AmphiSOXE is less efficient than SOX9 in transactivation activity and in the ability to preferentially promote glial over neuronal fate, a difference that lies within the combined properties of amino terminal and transactivation domains. We propose that acquisition of AmphiSoxE expression in the neural plate border led to NCC emergence while duplication and divergence produced advantageous mutations in vertebrate homologues, promoting elaboration of NC traits.

T-Box Genes and Developmental Gene Regulatory Networks in Ascidians

Ascidians are invertebrate chordates with a biphasic life cycle characterized by a dual body plan that displays simplified versions of chordate structures, such as a premetamorphic 40-cell notochord topped by a dorsal nerve cord and postmetamorphic pharyngeal slits. These relatively simple chordates are characterized by rapid development, compact genomes and ease of transgenesis, and thus provide the opportunity to rapidly characterize the genomic organization, developmental function, and transcriptional regulation of evolutionarily conserved gene families. This review summarizes the current knowledge on members of the T-box family of transcription factors in Ciona and other ascidians. In both chordate and nonchordate animals, these genes control a variety of morphogenetic processes, and their mutations are responsible for malformations and developmental defects in organisms ranging from flies to humans. In ascidians, T-box transcription factors are required for the formation and specialization of essential structures, including notochord, muscle, heart, and differentiated neurons. In recent years, the experimental advantages offered by ascidian embryos have allowed the rapid accumulation of a wealth of information on the molecular mechanisms that regulate the expression of T-box genes. These studies have also elucidated the strategies employed by these transcription factors to orchestrate the appropriate spatial and temporal deployment of the numerous target genes that they control.

Antagonism between β-catenin and Gata.a sequentially segregates the germ layers of ascidian embryos

Many animal embryos use nuclear β-catenin (nβ-catenin) during the segregation of endomesoderm (or endoderm) from ectoderm. This mechanism is thus likely to be evolutionarily ancient. In the ascidian embryo, nβ-catenin reiteratively drives binary fate decisions between ectoderm and endomesoderm at the 16-cell stage, and then between endoderm and margin (mesoderm and caudal neural) at the 32-cell stage. At the 16-cell stage, nβ-catenin activates endomesoderm genes in the vegetal hemisphere. At the same time, nβ-catenin suppresses the DNA-binding activity of a maternal transcription factor, Gata.a, through a physical interaction, and Gata.a thereby activates its target genes only in the ectodermal lineage. In the present study, we found that this antagonism between nβ-catenin and Gata.a also operates during the binary fate switch at the 32-cell stage. Namely, in marginal cells where nβ-catenin is absent, Gata.a directly activates its target, Zic-r.b (ZicL), to specify the marginal cell lineages. Thus, the antagonistic action between nβ-catenin and Gata.a is involved in two consecutive stages of germ layer segregation in ascidian embryos.

Fibronectin contributes to notochord intercalation in the invertebrate chordate, Ciona intestinalis

BACKGROUND

Genomic analysis has upended chordate phylogeny, placing the tunicates as the sister group to the vertebrates. This taxonomic rearrangement raises questions about the emergence of a tunicate/vertebrate ancestor.

RESULTS

Characterization of developmental genes uniquely shared by tunicates and vertebrates is one promising approach for deciphering developmental shifts underlying acquisition of novel, ancestral traits. The matrix glycoprotein Fibronectin (FN) has long been considered a vertebrate-specific gene, playing a major instructive role in vertebrate embryonic development. However, the recent computational prediction of an orthologous "vertebrate-like" Fn gene in the genome of a tunicate, Ciona savignyi, challenges this viewpoint suggesting that Fn may have arisen in the shared tunicate/vertebrate ancestor. Here we verify the presence of a tunicate Fn ortholog. Transgenic reporter analysis was used to characterize a Ciona Fn enhancer driving expression in the notochord. Targeted knockdown in the notochord lineage indicates that FN is required for proper convergent extension.

CONCLUSIONS

These findings suggest that acquisition of Fn was associated with altered notochord morphogenesis in the vertebrate/tunicate ancestor.

Automated detection of ncRNAs in the draft genome sequence of a colonial tunicate: the carpet sea squirt Didemnum vexillum

BACKGROUND

The colonial ascidian Didemnum vexillum, sea carpet squirt, is not only a key marine organism to study morphological ancestral patterns of chordates evolution but it is also of great ecological importance due to its status as a major invasive species. Non-coding RNAs, in particular microRNAs (miRNAs), are important regulatory genes that impact development and environmental adaptation. Beyond miRNAs, not much in known about tunicate ncRNAs.

RESULTS

We provide here a comprehensive homology-based annotation of non-coding RNAs in the recently sequenced genome of D. vexillum. To this end we employed a combination of several computational approaches, including blast searches with a wide range of parameters, and secondary structured centered survey with infernal. The resulting candidate set was curated extensively to produce a high-quality ncRNA annotation of the first draft of the D. vexillum genome. It comprises 57 miRNA families, 4 families of ribosomal RNAs, 22 isoacceptor classes of tRNAs (of which more than 72 % of loci are pseudogenes), 13 snRNAs, 12 snoRNAs, and 1 other RNA family. Additionally, 21 families of mitochondrial tRNAs and 2 of mitochondrial ribosomal RNAs and 1 long non-coding RNA.

CONCLUSIONS

The comprehensive annotation of the D. vexillum non-coding RNAs provides a starting point towards a better understanding of the restructuring of the small RNA system in ascidians. Furthermore it provides a valuable research for efforts to establish detailed non-coding RNA annotations for other recently published and recently sequences in tunicate genomes.

Ciona as a Simple Chordate Model for Heart Development and Regeneration

Cardiac cell specification and the genetic determinants that govern this process are highly conserved among Chordates. Recent studies have established the importance of evolutionarily-conserved mechanisms in the study of congenital heart defects and disease, as well as cardiac regeneration. As a basal Chordate, the Ciona model system presents a simple scaffold that recapitulates the basic blueprint of cardiac development in Chordates. Here we will focus on the development and cellular structure of the heart of the ascidian Ciona as compared to other Chordates, principally vertebrates. Comparison of the Ciona model system to heart development in other Chordates presents great potential for dissecting the genetic mechanisms that underlie congenital heart defects and disease at the cellular level and might provide additional insight into potential pathways for therapeutic cardiac regeneration.

The Simple Chordate Ciona intestinalis Has a Reduced Complement of Genes Associated with Fanconi Anemia

Fanconi anemia (FA) is a human genetic disease characterized by congenital defects, bone marrow failure, and increased cancer risk. FA is associated with mutation in one of 24 genes. The protein products of these genes function cooperatively in the FA pathway to orchestrate the repair of DNA interstrand cross-links. Few model organisms exist for the study of FA. Seeking a model organism with a simpler version of the FA pathway, we searched the genome of the simple chordate Ciona intestinalis for homologs of the human FA-associated proteins. BLAST searches, sequence alignments, hydropathy comparisons, maximum likelihood phylogenetic analysis, and structural modeling were used to infer the likelihood of homology between C. intestinalis and human FA proteins. Our analysis indicates that C. intestinalis indeed has a simpler and potentially functional FA pathway. The C. intestinalis genome was searched for candidates for homology to 24 human FA and FA-associated proteins. Support was found for the existence of homologs for 13 of these 24 human genes in C. intestinalis. Members of each of the three commonly recognized FA gene functional groups were found. In group I, we identified homologs of FANCE, FANCL, FANCM, and UBE2T/FANCT. Both members of group II, FANCD2 and FANCI, have homologs in C. intestinalis. In group III, we found evidence for homologs of FANCJ, FANCO, FANCQ/ERCC4, FANCR/RAD51, and FANCS/BRCA1, as well as the FA-associated proteins ERCC1 and FAN1. Evidence was very weak for the existence of homologs in C. intestinalis for any other recognized FA genes. This work supports the notion that C. intestinalis, as a close relative of vertebrates, but having a much reduced complement of FA genes, offers a means of studying the function of certain FA proteins in a simpler pathway than that of vertebrate cells.

Actin-Induced Structure in the Beta-Thymosin Family of Intrinsically Disordered Proteins

Thymosin β4 (Tβ4) is a 43-amino acid signature motif peptide that defines the beta-thymosin (βT) family of proteins. βTs are intrinsically unstructured in their free states and undergo disorder-to-order transitions in carrying out their biological functions. This property poses challenges in determining their 3D structures, mainly favoring structural studies on the complexes formed between βTs and their interaction partners. One of the βTs' primary binding partners is monomeric actin, a major component of the cytoskeleton in eukaryotic cells. Tβ4's role in this system is to maintain the highly concentrated pool of monomeric actin that can be accessed through profilin by actin filament nucleating machineries. Here, we give an account of the structures of βTs that have been illuminated by nuclear magnetic resonance (NMR) and X-ray crystallography. NMR has been the method of choice for probing regions that have intrinsic conformational preference within the largely disordered βTs in their native states in solution. X-ray crystallography has demonstrated at atomic detail how βTs interact with actin. Detailed analysis of these structures highlights the disorder-to-order transition of Tβ4 in binding to actin and its isoform specificity.

Evolution by gene loss

The recent increase in genomic data is revealing an unexpected perspective of gene loss as a pervasive source of genetic variation that can cause adaptive phenotypic diversity. This novel perspective of gene loss is raising new fundamental questions. How relevant has gene loss been in the divergence of phyla? How do genes change from being essential to dispensable and finally to being lost? Is gene loss mostly neutral, or can it be an effective way of adaptation? These questions are addressed, and insights are discussed from genomic studies of gene loss in populations and their relevance in evolutionary biology and biomedicine.

Functional diversification of sea urchin ABCC1 (MRP1) by alternative splicing

The multidrug resistance protein (MRP) family encodes a diverse repertoire of ATP-binding cassette (ABC) transporters with multiple roles in development, disease, and homeostasis. Understanding MRP evolution is central to unraveling their roles in these diverse processes. Sea urchins occupy an important phylogenetic position for understanding the evolution of vertebrate proteins and have been an important invertebrate model system for study of ABC transporters. We used phylogenetic analyses to examine the evolution of MRP transporters and functional approaches to identify functional forms of sea urchin MRP1 (also known as SpABCC1). SpABCC1, the only MRP homolog in sea urchins, is co-orthologous to human MRP1, MRP3, and MRP6 (ABCC1, ABCC3, and ABCC6) transporters. However, efflux assays revealed that alternative splicing of exon 22, a region critical for substrate interactions, could diversify functions of sea urchin MRP1. Phylogenetic comparisons also indicate that while MRP1, MRP3, and MRP6 transporters potentially arose from a single transporter in basal deuterostomes, alternative splicing appears to have been the major mode of functional diversification in invertebrates, while duplication may have served a more important role in vertebrates. These results provide a deeper understanding of the evolutionary origins of MRP transporters and the potential mechanisms used to diversify their functions in different groups of animals.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Tunicates: exploring the sea shores and roaming the open ocean. A tribute to Thomas Huxley

This review is a tribute to the remarkable contributions of Thomas Huxley to the biology of tunicates, the likely sister group of vertebrates. In 1851, the great biologist and philosopher published two landmark papers on pelagic tunicates in the Philosophical Transactions of the Royal Society. They were dedicated to the description of the adult anatomy and life cycle of thaliaceans and appendicularians, the pelagic relatives of ascidians. In the first part of this review, we discuss the novel anatomical observations and evolutionary hypotheses made by Huxley, which would have a lasting influence on tunicate biology. We also briefly comment on the more philosophical reflections of Huxley on individuality. In the second part, we stress the originality and relevance of past and future studies of tunicates in the resolution of major biological issues. In particular, we focus on the complex relationship between genotype and phenotype and the phenomenon of developmental system drift. We propose that more than 150 years after Huxley's papers, tunicate embryos are still worth studying in their own right, independently of their evolutionary proximity to vertebrates, as they provide original and crucial insights into the process of animal evolution. Tunicates are still at the forefront of biological research.

Genome-Wide Identification and Transcriptome-Based Expression Profiling of the Sox Gene Family in the Nile Tilapia (Oreochromis niloticus)

The Sox transcription factor family is characterized with the presence of a Sry-related high-mobility group (HMG) box and plays important roles in various biological processes in animals, including sex determination and differentiation, and the development of multiple organs. In this study, 27 Sox genes were identified in the genome of the Nile tilapia (Oreochromis niloticus), and were classified into seven groups. The members of each group of the tilapia Sox genes exhibited a relatively conserved exon-intron structure. Comparative analysis showed that the Sox gene family has undergone an expansion in tilapia and other teleost fishes following their whole genome duplication, and group K only exists in teleosts. Transcriptome-based analysis demonstrated that most of the tilapia Sox genes presented stage-specific and/or sex-dimorphic expressions during gonadal development, and six of the group B Sox genes were specifically expressed in the adult brain. Our results provide a better understanding of gene structure and spatio-temporal expression of the Sox gene family in tilapia, and will be useful for further deciphering the roles of the Sox genes during sex determination and gonadal development in teleosts.

Gut immunity in a protochordate involves a secreted immunoglobulin-type mediator binding host chitin and bacteria

Protochordate variable region-containing chitin-binding proteins (VCBPs) consist of immunoglobulin-type V domains and a chitin-binding domain (CBD). VCBP V domains facilitate phagocytosis of bacteria by granulocytic amoebocytes; the function of the CBD is not understood. Here we show that the gut mucosa of Ciona intestinalis contains an extensive matrix of chitin fibrils to which VCBPs bind early in gut development, before feeding. Later in development, VCBPs and bacteria colocalize to chitin-rich mucus along the intestinal wall. VCBP-C influences biofilm formation in vitro and, collectively, the findings of this study suggest that VCBP-C may influence the overall settlement and colonization of bacteria in the Ciona gut. Basic relationships between soluble immunoglobulin-type molecules, endogenous chitin and bacteria arose early in chordate evolution and are integral to the overall function of the gut barrier.

FastaHerder2: Four Ways to Research Protein Function and Evolution with Clustering and Clustered Databases

The accelerated growth of protein databases offers great possibilities for the study of protein function using sequence similarity and conservation. However, the huge number of sequences deposited in these databases requires new ways of analyzing and organizing the data. It is necessary to group the many very similar sequences, creating clusters with automated derived annotations useful to understand their function, evolution, and level of experimental evidence. We developed an algorithm called FastaHerder2, which can cluster any protein database, putting together very similar protein sequences based on near-full-length similarity and/or high threshold of sequence identity. We compressed 50 reference proteomes, along with the SwissProt database, which we could compress by 74.7%. The clustering algorithm was benchmarked using OrthoBench and compared with FASTA HERDER, a previous version of the algorithm, showing that FastaHerder2 can cluster a set of proteins yielding a high compression, with a lower error rate than its predecessor. We illustrate the use of FastaHerder2 to detect biologically relevant functional features in protein families. With our approach we seek to promote a modern view and usage of the protein sequence databases more appropriate to the postgenomic era.