Plasticity of animal genome architecture unmasked by rapid evolution of a pelagic tunicate

Sex-specific chromatin landscapes in an ultra-compact chordate genome

Background

In multicellular organisms, epigenome dynamics are associated with transitions in the cell cycle, development, germline specification, gametogenesis and inheritance. Evolutionarily, regulatory space has increased in complex metazoans to accommodate these functions. In tunicates, the sister lineage to vertebrates, we examine epigenome adaptations to strong secondary genome compaction, sex chromosome evolution and cell cycle modes.

Results

Across the 70 MB Oikopleura dioica genome, we profiled 19 histone modifications, and RNA polymerase II, CTCF and p300 occupancies, to define chromatin states within two homogeneous tissues with distinct cell cycle modes: ovarian endocycling nurse nuclei and mitotically proliferating germ nuclei in testes. Nurse nuclei had active chromatin states similar to other metazoan epigenomes, with large domains of operon-associated transcription, a general lack of heterochromatin, and a possible role of Polycomb PRC2 in dosage compensation. Testis chromatin states reflected transcriptional activity linked to spermatogenesis and epigenetic marks that have been associated with establishment of transgenerational inheritance in other organisms. We also uncovered an unusual chromatin state specific to the Y-chromosome, which combined active and heterochromatic histone modifications on specific transposable elements classes, perhaps involved in regulating their activity.

Conclusions

Compacted regulatory space in this tunicate genome is accompanied by reduced heterochromatin and chromatin state domain widths. Enhancers, promoters and protein-coding genes have conserved epigenomic features, with adaptations to the organization of a proportion of genes in operon units. We further identified features specific to sex chromosomes, cell cycle modes, germline identity and dosage compensation, and unusual combinations of histone PTMs with opposing consensus functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0110-4) contains supplementary material, which is available to authorized users.

Morphology, ecology, and molecular biology of a new species of giant larvacean in the eastern North Pacific: Bathochordaeus mcnutti sp. nov

Bathochordaeus mcnutti sp. nov. is described from the mesopelagic northeast Pacific Ocean (Monterey Bay, California, USA). Larvaceans in the genus Bathochordaeus are large, often abundant zooplankters found throughout much of the world ocean, but until recently it was unclear whether more than a single species of Bathochordaeus existed. Using remotely operated vehicles, we have made hundreds of in situ observations, compiled two decades of time-series data, and carefully collected enough specimens to determine that three species of Bathochordaeus occur in Monterey Bay: B. charon (Chun), B. stygius (Garstang), and B. mcnutti sp. nov. Bathochordaeus mcnutti is readily distinguished from its two congeners by the distinct blue outline visible around the periphery of its tail, and by other aspects of its morphology, ecology, and genetics. The abundance of larvaceans means they are ecologically important as particle processors. Species within the genus, Bathochordaeus, comprise the largest of described larvaceans.

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.

Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages

Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization. Data collected in the last decade have challenged this view by revealing considerable diversity in animal mitochondrial genome organization. Much of this diversity has been found in nonbilaterian animals (phyla Cnidaria, Ctenophora, Placozoa, and Porifera), which, from a phylogenetic perspective, form the main branches of the animal tree along with Bilateria. Within these groups, mt-genomes are characterized by varying numbers of both linear and circular chromosomes, extra genes (e.g. atp9, polB, tatC), large variation in the number of encoded mitochondrial transfer RNAs (tRNAs) (0-25), at least seven different genetic codes, presence/absence of introns, tRNA and mRNA editing, fragmented ribosomal RNA genes, translational frameshifting, highly variable substitution rates, and a large range of genome sizes. This newly discovered diversity allows a better understanding of the evolutionary plasticity and conservation of animal mtDNA and provides insights into the molecular and evolutionary mechanisms shaping mitochondrial genomes.

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.

Intron gain by tandem genomic duplication: a novel case in a potato gene encoding RNA-dependent RNA polymerase

The origin and subsequent accumulation of spliceosomal introns are prominent events in the evolution of eukaryotic gene structure. However, the mechanisms underlying intron gain remain unclear because there are few proven cases of recently gained introns. In an RNA-dependent RNA polymerase (RdRp) gene, we found that a tandem duplication occurred after the divergence of potato and its wild relatives among other Solanum plants. The duplicated sequence crosses the intron-exon boundary of the first intron and the second exon. A new intron was detected at this duplicated region, and it includes a small previously exonic segment of the upstream copy of the duplicated sequence and the intronic segment of the downstream copy of the duplicated sequence. The donor site of this new intron was directly obtained from the small previously exonic segment. Most of the splicing signals were inherited directly from the parental intron/exon structure, including a putative branch site, the polypyrimidine tract, the 3' splicing site, two putative exonic splicing enhancers, and the GC contents differed between the intron and exon. In the widely cited model of intron gain by tandem genomic duplication, the duplication of an AGGT-containing exonic segment provides the GT and AG splicing sites for the new intron. Our results illustrate that the tandem duplication model of intron gain should be diverse in terms of obtaining the proper splicing signals.

Ancient origins of metazoan gonadotropin-releasing hormone and their receptors revealed by phylogenomic analyses

The discovery of genes related to gonadotropin-releasing hormones (GnRH) and their receptors from diverse species has driven important advances in comparative endocrinology. However, our view of the evolutionary histories and nomenclature of these gene families has become inconsistent as several different iterations of GnRH and receptor relationships have been proposed. Whole genome sequence data are now available for most of the major lineages of animals, and an exhaustive view of the phylogenies of GnRH and their receptors is now possible. In this paper, we leverage data from publically available whole genome sequences to present a new phylogenomic analysis of GnRH and GnRH receptors and the distant relatives of each across metazoan phylogeny. Our approach utilizes a phylogenomics pipeline that searches data from 36 whole genome sequences and conducts phylogenetic analyses of gene trees. We provide a comprehensive analysis of the major groupings of GnRH peptides, related hormones and their receptors and provide some suggestions for a new nomenclature. Our study provides a framework for understanding the functional diversification of this family of neuromodulatory peptides and their receptors.

Coelimination and Survival in Gene Network Evolution: Dismantling the RA-Signaling in a Chordate

The bloom of genomics is revealing gene loss as a pervasive evolutionary force generating genetic diversity that shapes the evolution of species. Outside bacteria and yeast, however, the understanding of the process of gene loss remains elusive, especially in the evolution of animal species. Here, using the dismantling of the retinoic acid metabolic gene network (RA-MGN) in the chordate Oikopleura dioica as a case study, we combine approaches of comparative genomics, phylogenetics, biochemistry, and developmental biology to investigate the mutational robustness associated to biased patterns of gene loss. We demonstrate the absence of alternative pathways for RA-synthesis in O. dioica, which suggests that gene losses of RA-MGN were not compensated by mutational robustness, but occurred in a scenario of regressive evolution. In addition, the lack of drastic phenotypic changes associated to the loss of RA-signaling provides an example of the inverse paradox of Evo-Devo. This work illustrates how the identification of patterns of gene coelimination-in our case five losses (Rdh10, Rdh16, Bco1, Aldh1a, and Cyp26)-is a useful strategy to recognize gene network modules associated to distinct functions. Our work also illustrates how the identification of survival genes helps to recognize neofunctionalization events and ancestral functions. Thus, the survival and extensive duplication of Cco and RdhE2 in O. dioica correlated with the acquisition of complex compartmentalization of expression domains in the digestive system and a process of enzymatic neofunctionalization of the Cco, while the surviving Aldh8 could be related to its ancestral housekeeping role against toxic aldehydes.

"Wrecks of Ancient Life": Genetic Variants Vetted by Natural Selection

The Genetics Society of America’s George W. Beadle Award honors individuals who have made outstanding contributions to the community of genetics researchers and who exemplify the qualities of its namesake as a respected academic, administrator, and public servant. The 2015 recipient is John Postlethwait. He has made groundbreaking contributions in developing the zebrafish as a molecular genetic model and in understanding the evolution of new gene functions in vertebrates. He built the first zebrafish genetic map and showed that its genome, along with that of distantly related teleost fish, had been duplicated. Postlethwait played an integral role in the zebrafish genome-sequencing project and elucidated the genomic organization of several fish species. Postlethwait is also honored for his active involvement with the zebrafish community, advocacy for zebrafish as a model system, and commitment to driving the field forward.

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.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants

Background

Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification.

Results

We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages.

Conclusions

Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

Electronic supplementary material

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

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.

Flux Control in Glycolysis Varies Across the Tree of Life

Biochemical thought posits that rate-limiting steps (defined here as points of flux control) are strongly selected as points of pathway regulation and control and are thus expected to be evolutionarily conserved. Conversely, population genetic thought based upon the concepts of mutation-selection-drift balance at the pathway level might suggest variation in flux controlling steps over evolutionary time. Glycolysis, as one of the most conserved and best characterized pathways, was studied to evaluate its evolutionary conservation. The flux controlling step in glycolysis was found to vary over the tree of life. Further, phylogenetic analysis suggested at least 60 events of gene duplication and additional events of putative positive selection that might alter pathway kinetic properties. Together, these results suggest that even with presumed largely negative selection on pathway output on glycolysis, the co-evolutionary process under the hood is dynamic.

Models for gene duplication when dosage balance works as a transition state to subsequent neo-or sub-functionalization

Background

Dosage balance has been described as an important process for the retention of duplicate genes after whole genome duplication events. However, dosage balance is only a temporary mechanism for duplicate gene retention, as it ceases to function following the stochastic loss of interacting partners, as dosage balance itself is lost with this event. With the prolonged period of retention, on the other hand, there is the potential for the accumulation of substitutions which upon release from dosage balance constraints, can lead to either subsequent neo-functionalization or sub-functionalization. Mechanistic models developed to date for duplicate gene retention treat these processes independently, but do not describe dosage balance as a transition state to eventual functional change.

Results

Here a model for these processes (dosage plus neofunctionalization and dosage plus subfunctionalization) has been built within an existing framework. Because of the computational complexity of these models, a simpler modeling framework that captures the same information is also proposed. This model is integrated into a phylogenetic birth-death model, expanding the range of available models.

Conclusions

Including further levels of biological reality in methods for gene tree/species tree reconciliation should not only increase the accuracy of estimates of the timing and evolutionary history of genes but can also offer insight into how genes and genomes evolve. These new models add to the tool box for characterizing mechanisms of duplicate gene retention probabilistically.

Discovery of sea urchin NGFFFamide receptor unites a bilaterian neuropeptide family

Neuropeptides are ancient regulators of physiology and behaviour, but reconstruction of neuropeptide evolution is often difficult owing to lack of sequence conservation. Here, we report that the receptor for the neuropeptide NGFFFamide in the sea urchin Strongylocentrotus purpuratus (phylum Echinodermata) is an orthologue of vertebrate neuropeptide-S (NPS) receptors and crustacean cardioactive peptide (CCAP) receptors. Importantly, this has facilitated reconstruction of the evolution of two bilaterian neuropeptide signalling systems. Genes encoding the precursor of a vasopressin/oxytocin-type neuropeptide and its receptor duplicated in a common ancestor of the Bilateria. One copy of the precursor retained ancestral features, as seen in highly conserved vasopressin/oxytocin–neurophysin-type precursors. The other copy diverged, but this took different courses in protostomes and deuterostomes. In protostomes, the occurrence of a disulfide bridge in neuropeptide product(s) of the precursor was retained, as in CCAP, but with loss of the neurophysin domain. In deuterostomes, we see the opposite scenario—the neuropeptides lost the disulfide bridge, and neurophysin was retained (as in the NGFFFamide precursor) but was subsequently lost in vertebrate NPS precursors. Thus, the sea urchin NGFFFamide precursor and receptor are ‘missing links’ in the evolutionary history of neuropeptides that control ecdysis in arthropods (CCAP) and regulate anxiety in humans (NPS).

Extensive Mitochondrial mRNA Editing and Unusual Mitochondrial Genome Organization in Calcaronean Sponges

One of the unusual features of DNA-containing organelles in general and mitochondria in particular is the frequent occurrence of RNA editing [1]. The term "RNA editing" refers to a variety of mechanistically unrelated biochemical processes that alter RNA sequence during or after transcription [2]. The editing can be insertional, deletional, or substitutional and has been found in all major types of RNAs [3, 4]. Although mitochondrial mRNA editing is widespread in some eukaryotic lineages [5-7], it is rare in animals, with reported cases limited both in their scope and in phylogenetic distribution [8-11] (see also [12]). While analyzing genomic data from calcaronean sponges Sycon ciliatum and Leucosolenia complicata, we were perplexed by the lack of recognizable mitochondrial coding sequences. Comparison of genomic and transcriptomic data from these species revealed the presence of mitochondrial cryptogenes whose transcripts undergo extensive editing. This editing consisted of single or double uridylate (U) insertions in pre-existing short poly(U) tracts. Subsequent analysis revealed the presence of similar editing in Sycon coactum and the loss of editing in Petrobiona massiliana, a hypercalcified calcaronean sponge. In addition, mitochondrial genomes of at least some calcaronean sponges were found to have a highly unusual architecture, with nearly all genes located on individual and likely linear chromosomes. Phylogenetic analysis of mitochondrial coding sequences revealed accelerated rates of sequence evolution in this group. The latter observation presents a challenge for the mutational-hazard hypothesis [13], which posits that mRNA editing should not occur in lineages with an elevated mutation rate.

Evolutionary loss of peroxisomes--not limited to parasites

Background

Peroxisomes are ubiquitous eukaryotic organelles that compartmentalize a variety of metabolic pathways that are primarily related to the oxidative metabolism of lipids and the detoxification of reactive oxygen species. The importance of peroxisomes is underscored by serious human diseases, which are caused by disorders in peroxisomal functions. Some eukaryotic lineages, however, lost peroxisomes. These organisms are mainly anaerobic protists and some parasitic lineages including Plasmodium and parasitic platyhelminths. Here we performed a systematic in-silico analysis of peroxisomal markers among metazoans to assess presence of peroxisomes and peroxisomal enzymes.

Results

Our analyses reveal an obvious loss of peroxisomes in all tested flukes, tapeworms, and parasitic roundworms of the order Trichocephalida. Intriguingly, peroxisomal markers are absent from the genome of the free-living tunicate Oikopleura dioica, which inhabits oxygen-containing niches of sea waters. We further map the presence and predicted subcellular localization of putative peroxisomal enzymes, showing that in organisms without the peroxisomal markers the set of these enzymes is highly reduced and none of them contains a predicted peroxisomal targeting signal.

Conclusions

We have shown that several lineages of metazoans independently lost peroxisomes and that the loss of peroxisomes was not exclusively associated with adaptation to anaerobic habitats and a parasitic lifestyle. Although the reason for the loss of peroxisomes from O. dioica is unclear, organisms lacking peroxisomes, including the free-living O. dioica, share certain typical r-selected traits: high fecundity, limited ontogenesis and relatively low complexity of the gene content. We hypothesize that peroxisomes are generally the first compartment to be lost during evolutionary reductions of the eukaryotic cell.

Reviewers

This article was reviewed by Michael Gray and Nick Lane.

Electronic supplementary material

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

Extracting phylogenetic signal and accounting for bias in whole-genome data sets supports the Ctenophora as sister to remaining Metazoa

BACKGROUND

Understanding the phylogenetic relationships among major lineages of multicellular animals (the Metazoa) is a prerequisite for studying the evolution of complex traits such as nervous systems, muscle tissue, or sensory organs. Transcriptome-based phylogenies have dramatically improved our understanding of metazoan relationships in recent years, although several important questions remain. The branching order near the base of the tree, in particular the placement of the poriferan (sponges, phylum Porifera) and ctenophore (comb jellies, phylum Ctenophora) lineages is one outstanding issue. Recent analyses have suggested that the comb jellies are sister to all remaining metazoan phyla including sponges. This finding is surprising because it suggests that neurons and other complex traits, present in ctenophores and eumetazoans but absent in sponges or placozoans, either evolved twice in Metazoa or were independently, secondarily lost in the lineages leading to sponges and placozoans.

RESULTS

To address the question of basal metazoan relationships we assembled a novel dataset comprised of 1080 orthologous loci derived from 36 publicly available genomes representing major lineages of animals. From this large dataset we procured an optimized set of partitions with high phylogenetic signal for resolving metazoan relationships. This optimized data set is amenable to the most appropriate and computationally intensive analyses using site-heterogeneous models of sequence evolution. We also employed several strategies to examine the potential for long-branch attraction to bias our inferences. Our analyses strongly support the Ctenophora as the sister lineage to other Metazoa. We find no support for the traditional view uniting the ctenophores and Cnidaria. Our findings are supported by Bayesian comparisons of topological hypotheses and we find no evidence that they are biased by long-branch attraction.

CONCLUSIONS

Our study further clarifies relationships among early branching metazoan lineages. Our phylogeny supports the still-controversial position of ctenophores as sister group to all other metazoans. This study also provides a workflow and computational tools for minimizing systematic bias in genome-based phylogenetic analyses. Future studies of metazoan phylogeny will benefit from ongoing efforts to sequence the genomes of additional invertebrate taxa that will continue to inform our view of the relationships among the major lineages of animals.

Molecular signatures that are distinctive characteristics of the vertebrates and chordates and supporting a grouping of vertebrates with the tunicates

Members of the phylum Chordata and the subphylum Vertebrata are presently distinguished solely on the basis of morphological characteristics. The relationship of the vertebrates to the two non-vertebrate chordate subphyla is also a subject of debate. Analyses of protein sequences have identified multiple conserved signature indels (CSIs) that are specific for Chordata or for Vertebrata. Five CSIs in 4 important proteins are specific for the Vertebrata, whereas two other CSIs are uniquely found in all sequenced chordate species including Ciona intestinalis and Oikapleura dioica (Tunicates) as well as Branchiostoma floridae (Cephalochordates). The shared presence of these molecular signatures by all vertebrates/chordate species, but in no other animal taxa, strongly indicates that the genetic changes represented by the identified CSIs diagnose monophyletic groups. Two other discovered CSIs are uniquely shared by different vertebrate species and by either one (Ciona intestinalis) or both tunicate (Ciona and Oikapleura) species, but they are not found in Branchiostoma or other animal species. Specific presence of these CSIs in different vertebrates and either one or both tunicate species provides strong independent evidence that the vertebrate species are more closely related to the urochordates (tunicates) than to the cephalochordates.

Comparative analyses of developmental transcription factor repertoires in sponges reveal unexpected complexity of the earliest animals

Developmental transcription factors (DTFs) control development of animals by affecting expression of target genes, some of which are transcription factors themselves. In bilaterians and cnidarians, conserved DTFs are involved in homologous processes such as gastrulation or specification of neurons. The genome of Amphimedon queenslandica, the first sponge to be sequenced, revealed that only a fraction of these conserved DTF families are present in demosponges. This finding was in line with the view that morphological complexity in the animal lineage correlates with developmental toolkit complexity. However, as the phylum Porifera is very diverse, Amphimedon's genome may not be representative of all sponges. The recently sequenced genomes of calcareous sponges Sycon ciliatum and Leucosolenia complicata allowed investigations of DTFs in a sponge lineage evolutionarily distant from demosponges. Surprisingly, the phylogenetic analyses of identified DTFs revealed striking differences between the calcareous sponges and Amphimedon. As these differences appear to be a result of independent gene loss events in the two sponge lineages, the last common ancestor of sponges had to possess a much more diverse repertoire of DTFs than extant sponges. Developmental expression of sponge homologs of genes involved in specification of the Bilaterian endomesoderm and the neurosensory cells suggests that roles of many DTFs date back to the last common ancestor of all animals. Strikingly, even DTFs displaying apparent pan-metazoan conservation of sequence and function are not immune to being lost from individual species genomes. The quest for a comprehensive picture of the developmental toolkit in the last common metazoan ancestor is thus greatly benefitting from the increasing accessibility of sequencing, allowing comparisons of multiple genomes within each phylum.

Do echinoderm genomes measure up?

Echinoderm genome sequences are a corpus of useful information about a clade of animals that serve as research models in fields ranging from marine ecology to cell and developmental biology. Genomic information from echinoids has contributed to insights into the gene interactions that drive the developmental process at the molecular level. Such insights often rely heavily on genomic information and the kinds of questions that can be asked thus depend on the quality of the sequence information. Here we describe the history of echinoderm genomic sequence assembly and present details about the quality of the data obtained. All of the sequence information discussed here is posted on the echinoderm information web system, Echinobase.org.

Evolutionary volatile Cysteines and protein disorder in the fast evolving tunicate Oikopleura dioica

Cysteine (Cys) is regarded as the most conservative amino acid in nature, something that does not occur in the tunicate Oikopleura dioica, where this amino acid is one of the fastest evolving. In this work we analyze some of the causes of this intriguing absence of conservation. Considering the well-known stabilizing role of Cys, it was first investigated whether the lack of conservation was accompanied by an increase in intrinsic protein disorder. In contrast to expectations, it was found that O. dioica is the chordate that has the lowest levels of intrinsic disorder, while vertebrates (represented by Bos taurus) contain the most disordered proteins. Oikopleura proteins are shorter than their homologs in other Chordates (Ciona and B. taurus proteins are respectively 11% and 18% longer). This process of protein shortening was more intense in intrinsic disordered regions. As a result proteins became not only shorter but also more compact. It is also reported here that the conservation/divergence behavior of Cys depends on whether they are located in ordered or disordered regions. In the four species analyzed, disordered Cys are majorly (> 75%) not conserved at all. Ordered Cys instead, are much more free to diverge in Oikopleura than in the other chordates. We hypothesize that the preferential deletion of disordered regions resulted in a decreased protein disorder and a direct elimination (by deletion) of many ancestral Cys. Besides, the alterations (shortening or complete elimination) of some disordered regions (loops/random coils) probably promoted further Cys evolutionary volatility, because some ancestral Cys (and other amino acids which play a role in stability like Trp) located outside deleted regions became redundant due to the loss of their stabilizing partners.

Trans-splicing in metazoans: A link to translational control?

The trans-splicing of a spliced-leader RNA to a subset of mRNAs is a phenomenon that occurs in many species, including Caenorhabditis elegans, and yet the driving force for its evolution in disparate groups of animals remains unclear. Polycistronic mRNA resulting from the transcription of operons is resolved via trans-splicing, but operons comprise only a sub-set of trans-spliced genes. Using the marine chordate, Oikopleura dioica, we recently tested the hypothesis that metazoan operons accelerate recovery from growth arrest. We found no supporting evidence for this in O. dioica. Instead we found a striking relationship between trans-splicing and maternal mRNA in O. dioica, C. elegans and the ascidian, Ciona intestinalis. Furthermore, in O. dioica and C. elegans, we found evidence to suggest a role for mTOR signaling in the translational control of growth-related, trans-spliced maternal mRNAs. We propose that this may be a mechanism for adjusting egg number in response to nutrient levels in these species.

Novel Introner-Like Elements in fungi Are Involved in Parallel Gains of Spliceosomal Introns

Spliceosomal introns are key components of the eukaryotic gene structure. Although they contributed to the emergence of eukaryotes, their origin remains elusive. In fungi, they might originate from the multiplication of invasive introns named Introner-Like Elements (ILEs). However, so far ILEs have been observed in six fungal species only, including Fulvia fulva and Dothistroma septosporum (Dothideomycetes), arguing against ILE insertion as a general mechanism for intron gain. Here, we identified novel ILEs in eight additional fungal species that are phylogenetically related to F. fulva and D. septosporum using PCR amplification with primers derived from previously identified ILEs. The ILE content appeared unique to each species, suggesting independent multiplication events. Interestingly, we identified four genes each containing two gained ILEs. By analysing intron positions in orthologues of these four genes in Ascomycota, we found that three ILEs had inserted within a 15 bp window that contains regular spliceosomal introns in other fungal species. These three positions are not the result of intron sliding because ILEs are newly gained introns. Furthermore, the alternative hypothesis of an inferred ancestral gain followed by independent losses contradicts the observed degeneration of ILEs. These observations clearly indicate three parallel intron gains in four genes that were randomly identified. Our findings suggest that parallel intron gain is a phenomenon that has been highly underestimated in ILE-containing fungi, and likely in the whole fungal kingdom.

Maternal and zygotic transcriptomes in the appendicularian, Oikopleura dioica: novel protein-encoding genes, intra-species sequence variations, and trans-spliced RNA leader

RNA sequencing analysis was carried out to characterize egg and larval transcriptomes in the appendicularian, Oikopleura dioica, a planktonic chordate, which is characterized by rapid development and short life cycle of 5 days, using a Japanese population of the organism. De novo transcriptome assembly matched with 16,423 proteins corresponding to 95.4% of the protein-encoding genes deposited in the OikoBase, the genome database of the Norwegian population. Nucleotide and amino acid sequence identities between the Japanese and Norwegian O. dioica were estimated to be around 91.0 and 94.8%, respectively. We discovered 175 novel protein-encoding genes: 144 unigenes were common to both the Japanese and Norwegian populations, whereas 31 unigenes were not found in the OikoBase genome reference. Among the total 12,311 unigenes, approximately 63% were detected in egg-stage RNAs, whereas 99% were detected in larval stage RNAs; 3772 genes were up-regulated, and 1336 genes were down-regulated more than four-fold in the larvae. Gene ontology analyses characterized gene activities in these two developmental stages. We found a messenger RNA (mRNA) 5' trans-spliced leader, which was observed in 40.8% of the total unique transcripts. It showed preferential linkage to adenine at the 5' ends of the downstream exons. Trans-splicing was observed more frequently in egg mRNAs compared with larva-specific mRNAs.

Intron Invasions Trace Algal Speciation and Reveal Nearly Identical Arctic and Antarctic Micromonas Populations

Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence–absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture (Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica—a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga.

Embryonic expression of endogenous retroviral RNAs in somatic tissues adjacent to the Oikopleura germline

Selective pressure to maintain small genome size implies control of transposable elements, and most old classes of retrotransposons are indeed absent from the very compact genome of the tunicate Oikopleura dioica. Nonetheless, two families of retrotransposons are present, including the Tor elements. The gene organization within Tor elements is similar to that of LTR retrotransposons and retroviruses. In addition to gag and pol, many Tor elements carry a third gene encoding viral envelope-like proteins (Env) that may mediate infection. We show that the Tor family contains distinct classes of elements. In some classes, env mRNA is transcribed from the 5′LTR as in retroviruses. In others, env is transcribed from an additional promoter located downstream of the 5′LTR. Tor Env proteins are membrane-associated glycoproteins which exhibit some features of viral membrane fusion proteins. Whereas some elements are expressed in the adult testis, many others are specifically expressed in embryonic somatic cells adjacent to primordial germ cells. Such embryonic expression depends on determinants present in the Tor elements and not on their surrounding genomic environment. Our study shows that unusual modes of transcription and expression close to the germline may contribute to the proliferation of Tor elements.

Evolution of the metazoan mitochondrial replicase

The large number of complete mitochondrial DNA (mtDNA) sequences available for metazoan species makes it a good system for studying genome diversity, although little is known about the mechanisms that promote and/or are correlated with the evolution of this organellar genome. By investigating the molecular evolutionary history of the catalytic and accessory subunits of the mtDNA polymerase, pol γ, we sought to develop mechanistic insight into its function that might impact genome structure by exploring the relationships between DNA replication and animal mitochondrial genome diversity. We identified three evolutionary patterns among metazoan pol γs. First, a trend toward stabilization of both sequence and structure occurred in vertebrates, with both subunits evolving distinctly from those of other animal groups, and acquiring at least four novel structural elements, the most important of which is the HLH-3β (helix-loop-helix, 3 β-sheets) domain that allows the accessory subunit to homodimerize. Second, both subunits of arthropods and tunicates have become shorter and evolved approximately twice as rapidly as their vertebrate homologs. And third, nematodes have lost the gene for the accessory subunit, which was accompanied by the loss of its interacting domain in the catalytic subunit of pol γ, and they show the highest rate of molecular evolution among all animal taxa. These findings correlate well with the mtDNA genomic features of each group described above, and with their modes of DNA replication, although a substantive amount of biochemical work is needed to draw conclusive links regarding the latter. Describing the parallels between evolution of pol γ and metazoan mtDNA architecture may also help in understanding the processes that lead to mitochondrial dysfunction and to human disease-related phenotypes.

Diversity of animal immune receptors and the origins of recognition complexity in the deuterostomes

Invertebrate animals are characterized by extraordinary diversity in terms of body plan, life history and life span. The past impression that invertebrate immune responses are controlled by relatively simple innate systems is increasingly contradicted by genomic analyses that reveal significant evolutionary novelty and complexity. One accessible measure of this complexity is the multiplicity of genes encoding homologs of pattern recognition receptors. These multigene families vary significantly in size, and their sequence character suggests that they vary in function. At the same time, certain aspects of downstream signaling appear to be conserved. Here, we analyze five major classes of immune recognition receptors from newly available animal genome sequences. These include the Toll-like receptors (TLR), Nod-like receptors (NLR), SRCR domain scavenger receptors, peptidoglycan recognition proteins (PGRP), and Gram negative binding proteins (GNBP). We discuss innate immune complexity in the invertebrate deuterostomes, which was first recognized in sea urchins, within the wider context of emerging genomic information across animal phyla.

Trans-splicing and operons in metazoans: translational control in maternally regulated development and recovery from growth arrest

Polycistronic mRNAs transcribed from operons are resolved via the trans-splicing of a spliced-leader (SL) RNA. Trans-splicing also occurs at monocistronic transcripts. The phlyogenetically sporadic appearance of trans-splicing and operons has made the driving force(s) for their evolution in metazoans unclear. Previous work has proposed that germline expression drives operon organization in Caenorhabditis elegans, and a recent hypothesis proposes that operons provide an evolutionary advantage via the conservation of transcriptional machinery during recovery from growth arrested states. Using a modified cap analysis of gene expression protocol we mapped sites of SL trans-splicing genome-wide in the marine chordate Oikopleura dioica. Tiled microarrays revealed the expression dynamics of trans-spliced genes across development and during recovery from growth arrest. Operons did not facilitate recovery from growth arrest in O. dioica. Instead, we found that trans-spliced transcripts were predominantly maternal. We then analyzed data from C. elegans and Ciona intestinalis and found that an enrichment of trans-splicing and operon gene expression in maternal mRNA is shared between all three species, suggesting that this may be a driving force for operon evolution in metazoans. Furthermore, we found that the majority of known terminal oligopyrimidine (TOP) mRNAs are trans-spliced in O. dioica and that the SL contains a TOP-like motif. This suggests that the SL in O. dioica confers nutrient-dependent translational control to trans-spliced mRNAs via the TOR-signaling pathway. We hypothesize that SL-trans-splicing provides an evolutionary advantage in species that depend on translational control for regulating early embryogenesis, growth and oocyte production in response to nutrient levels.

Recurrent loss of specific introns during angiosperm evolution

Numerous instances of presence/absence variations for introns have been documented in eukaryotes, and some cases of recurrent loss of the same intron have been suggested. However, there has been no comprehensive or phylogenetically deep analysis of recurrent intron loss. Of 883 cases of intron presence/absence variation that we detected in five sequenced grass genomes, 93 were confirmed as recurrent losses and the rest could be explained by single losses (652) or single gains (118). No case of recurrent intron gain was observed. Deep phylogenetic analysis often indicated that apparent intron gains were actually numerous independent losses of the same intron. Recurrent loss exhibited extreme non-randomness, in that some introns were removed independently in many lineages. The two larger genomes, maize and sorghum, were found to have a higher rate of both recurrent loss and overall loss and/or gain than foxtail millet, rice or Brachypodium. Adjacent introns and small introns were found to be preferentially lost. Intron loss genes exhibited a high frequency of germ line or early embryogenesis expression. In addition, flanking exon A+T-richness and intron TG/CG ratios were higher in retained introns. This last result suggests that epigenetic status, as evidenced by a loss of methylated CG dinucleotides, may play a role in the process of intron loss. This study provides the first comprehensive analysis of recurrent intron loss, makes a series of novel findings on the patterns of recurrent intron loss during the evolution of the grass family, and provides insight into the molecular mechanism(s) underlying intron loss.

The spliceosomal introns are nucleotide sequences that interrupt coding regions of eukaryotic genes and are removed by RNA splicing after transcription. Recent studies have reported several examples of possible recurrent intron loss or gain, i.e., introns that are independently removed from or inserted into the identical sites more than once in an investigated phylogeny. However, the frequency, evolutionary patterns or other characteristics of recurrent intron turnover remain unknown. We provide results for the first comprehensive analysis of recurrent intron turnover within a plant family and show that recurrent intron loss represents a considerable portion of all intron losses identified and intron loss events far outnumber intron gain events. We also demonstrate that recurrent intron loss is non-random, affecting only a small number of introns that are repeatedly lost, and that different lineages show significantly different rates of intron loss. Our results suggest a possible role of DNA methylation in the process of intron loss. Moreover, this study provides strong support for the model of intron loss by reverse transcriptase mediated conversion of genes by their processed mRNA transcripts.

Retroduplication and loss of parental genes is a mechanism for the generation of intronless genes in Ciona intestinalis and Ciona savignyi

Tunicates, the sister clade of vertebrates, have miniature genomes and numerous intronless genes compared to other animals. It is still unclear how the tunicates acquired such a large number of intronless genes. Here, we analyzed sequences and intron-exon organizations of homologous genes from two closely related tunicates, Ciona intestinalis and Ciona savignyi. We found seven cases in which ancestral introns of a gene were completely lost in a species after their divergence. In four cases, both the intronless copy and the intron-containing copy were present in the genome, indicating that the intronless copy was generated by retroduplication. In the other three cases, the intron-containing copy was absent, implying it was lost after retroduplication. This result suggests that retroduplication and loss of parental genes is a major mechanism for the accumulation of intronless genes in tunicates.

Mutation studies in ascidians: a review

Historically, mutations have had a significant impact on the study of developmental processes and phenotypic evolution. Lesions in DNA are created by artificial methods or detected by natural genetic variation. Random mutations are then ascribed to genetic change by direct sequencing or positional cloning. Tunicate species of the ascidian genus Ciona represent nearly fully realized model systems in which gene function can be investigated in depth. Additionally, tunicates are valuable organisms for the study of naturally occurring mutations due to the capability to exploit genetic variation down to the molecular level. Here, we summarize the available information about how mutations are studied in ascidians with examples of insights that have resulted from these applications. We also describe notions and methodologies that might be useful for the implementation of easy and tight procedures for mutations studies in Ciona.

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

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

Arthropods are the most abundant animals on earth. Among them, insects clearly dominate on land, whereas crustaceans hold the title for the most diverse invertebrates in the oceans. Much is known about the biology of these groups, not least because of genomic studies of the fruit fly Drosophila, the water flea Daphnia, and other species used in research. Here we report the first genome sequence from a species belonging to a lineage that has previously received very little attention—the myriapods. Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes, one of which—Strigamia maritima—we have sequenced here. We find that the genome of this centipede retains more characteristics of the presumed arthropod ancestor than other sequenced insect genomes. The genome provides access to many aspects of myriapod biology that have not been studied before, suggesting, for example, that they have diversified receptors for smell that are quite different from those used by insects. In addition, it shows specific consequences of the largely subterranean life of this particular species, which seems to have lost the genes for all known light-sensing molecules, even though it still avoids light.

Marine invertebrate xenobiotic-activated nuclear receptors: their application as sensor elements in high-throughput bioassays for marine bioactive compounds

Developing high-throughput assays to screen marine extracts for bioactive compounds presents both conceptual and technical challenges. One major challenge is to develop assays that have well-grounded ecological and evolutionary rationales. In this review we propose that a specific group of ligand-activated transcription factors are particularly well-suited to act as sensors in such bioassays. More specifically, xenobiotic-activated nuclear receptors (XANRs) regulate transcription of genes involved in xenobiotic detoxification. XANR ligand-binding domains (LBDs) may adaptively evolve to bind those bioactive, and potentially toxic, compounds to which organisms are normally exposed to through their specific diets. A brief overview of the function and taxonomic distribution of both vertebrate and invertebrate XANRs is first provided. Proof-of-concept experiments are then described which confirm that a filter-feeding marine invertebrate XANR LBD is activated by marine bioactive compounds. We speculate that increasing access to marine invertebrate genome sequence data, in combination with the expression of functional recombinant marine invertebrate XANR LBDs, will facilitate the generation of high-throughput bioassays/biosensors of widely differing specificities, but all based on activation of XANR LBDs. Such assays may find application in screening marine extracts for bioactive compounds that could act as drug lead compounds.

Ascidian mitogenomics: comparison of evolutionary rates in closely related taxa provides evidence of ongoing speciation events

Ascidians are a fascinating group of filter-feeding marine chordates characterized by rapid evolution of both sequences and structure of their nuclear and mitochondrial genomes. Moreover, they include several model organisms used to investigate complex biological processes in chordates. To study the evolutionary dynamics of ascidians at short phylogenetic distances, we sequenced 13 new mitogenomes and analyzed them, together with 15 other available mitogenomes, using a novel approach involving detailed whole-mitogenome comparisons of conspecific and congeneric pairs. The evolutionary rate was quite homogeneous at both intraspecific and congeneric level, and the lowest congeneric rates were found in cryptic (morphologically undistinguishable) and in morphologically very similar species pairs. Moreover, congeneric nonsynonymous rates (dN) were up to two orders of magnitude higher than in intraspecies pairs. Overall, a clear-cut gap sets apart conspecific from congeneric pairs. These evolutionary peculiarities allowed easily identifying an extraordinary intraspecific variability in the model ascidian Botryllus schlosseri, where most pairs show a dN value between that observed at intraspecies and congeneric level, yet consistently lower than that of the Ciona intestinalis cryptic species pair. These data suggest ongoing speciation events producing genetically distinct B. schlosseri entities. Remarkably, these ongoing speciation events were undetectable by the cox1 barcode fragment, demonstrating that, at low phylogenetic distances, the whole mitogenome has a higher resolving power than cox1. Our study shows that whole-mitogenome comparative analyses, performed on a suitable sample of congeneric and intraspecies pairs, may allow detecting not only cryptic species but also ongoing speciation events.

Intron evolution in Neurospora: the role of mutational bias and selection

We used comparative and population genomics to study intron evolutionary dynamics in the fungal model genus Neurospora. For our investigation, we used well-annotated genomes of N. crassa, N. discreta, and N. tetrasperma, and 92 resequenced genomes of N. tetrasperma from natural populations. By analyzing the four well-annotated genomes, we identified 9495 intron sites in 7619 orthologous genes. Our data supports nonhomologous end joining (NHEJ) and tandem duplication as mechanisms for intron gains in the genus and the RT-mRNA process as a mechanism for intron loss. We found a moderate intron gain rate (5.78–6.89 × 10⁻¹³ intron gains per nucleotide site per year) and a high intron loss rate (7.53–13.76 × 10⁻¹⁰ intron losses per intron sites per year) as compared to other eukaryotes. The derived intron gains and losses are skewed to high frequencies, relative to neutral SNPs, in natural populations of N. tetrasperma, suggesting that selection is involved in maintaining a high intron turnover. Furthermore, our analyses of the association between intron population-level frequency and genomic features suggest that selection is involved in shaping a 5′ intron position bias and a low intron GC content. However, intron sequence analyses suggest that the gained introns were not exposed to recent selective sweeps. Taken together, this work contributes to our understanding of the importance of mutational bias and selection in shaping the intron distribution in eukaryotic genomes.

Guidelines for the nomenclature of genetic elements in tunicate genomes

Tunicates are invertebrate members of the chordate phylum, and are considered to be the sister group of vertebrates. Tunicates are composed of ascidians, thaliaceans, and appendicularians. With the advent of inexpensive high-throughput sequencing, the number of sequenced tunicate genomes is expected to rise sharply within the coming years. To facilitate comparative genomics within the tunicates, and between tunicates and vertebrates, standardized rules for the nomenclature of tunicate genetic elements need to be established. Here we propose a set of nomenclature rules, consensual within the community, for predicted genes, pseudogenes, transcripts, operons, transcriptional cis-regulatory regions, transposable elements, and transgenic constructs. In addition, the document proposes guidelines for naming transgenic and mutant lines.