The C. savignyi genetic map and its integration with the reference sequence facilitates insights into chordate genome evolution

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

The Genome of the "Sea Vomit" Didemnum vexillum

Tunicates are the sister group of vertebrates and thus occupy a key position for investigations into vertebrate innovations as well as into the consequences of the vertebrate-specific genome duplications. Nevertheless, tunicate genomes have not been studied extensively in the past, and comparative studies of tunicate genomes have remained scarce. The carpet sea squirt Didemnum vexillum, commonly known as "sea vomit", is a colonial tunicate considered an invasive species with substantial ecological and economical risk. We report the assembly of the D. vexillum genome using a hybrid approach that combines 28.5 Gb Illumina and 12.35 Gb of PacBio data. The new hybrid scaffolded assembly has a total size of 517.55 Mb that increases contig length about eightfold compared to previous, Illumina-only assembly. As a consequence of an unusually high genetic diversity of the colonies and the moderate length of the PacBio reads, presumably caused by the unusually acidic milieu of the tunic, the assembly is highly fragmented (L50 = 25,284, N50 = 6539). It is sufficient, however, for comprehensive annotations of both protein-coding genes and non-coding RNAs. Despite its shortcomings, the draft assembly of the "sea vomit" genome provides a valuable resource for comparative tunicate genomics and for the study of the specific properties of colonial ascidians.

Chromosomal Inversion Polymorphisms in Two Sympatric Ascidian Lineages

Chromosomal rearrangements can reduce fitness of heterozygotes and can thereby prevent gene flow. Therefore, such rearrangements can play a role in local adaptation and speciation. In particular, inversions are considered to be a major potential cause for chromosomal speciation. There are two closely related, partially sympatric lineages of ascidians in the genus Ciona, which we call type-A and type-B animals in the present study. Although these invertebrate chordates are largely isolated reproductively, hybrids can be found in wild populations, suggesting incomplete prezygotic barriers. Although the genome of type-A animals has been decoded and widely used, the genome for type-B animals has not been decoded at the chromosomal level. In the present study, we sequenced the genomes of two type-B individuals from different sides of the English Channel (in the zone of sympatry with type-A individuals) and compared them at the chromosomal level with the type-A genome. Although the overall structures were well conserved between type A and type B, chromosomal alignments revealed many inversions differentiating these two types of Ciona; it is probable that the frequent inversions have contributed to separation between these two lineages. In addition, comparisons of the genomes between the two type-B individuals revealed that type B had high rates of inversion polymorphisms and nucleotide polymorphisms, and thus type B might be in the process of differentiation into multiple new types or species. Our results suggest an important role of inversions in chromosomal speciation of these broadcasting spawners.

A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

A draft genome of the striped catfish, Pangasianodon hypophthalmus, for comparative analysis of genes relevant to development and a resource for aquaculture improvement

Background

The striped catfish, Pangasianodon hypophthalmus, is a freshwater and benthopelagic fish common in the Mekong River delta. Catfish constitute a valuable source of dietary protein. Therefore, they are cultured worldwide, and P. hypophthalmus is a food staple in the Mekong area. However, genetic information about the culture stock, is unavailable for breeding improvement, although genetics of the channel catfish, Ictalurus punctatus, has been reported. To acquire genome sequence data as a useful resource for marker-assisted breeding, we decoded a draft genome of P. hypophthalmus and performed comparative analyses.

Results

Using the Illumina platform, we obtained both nuclear and mitochondrial DNA sequences. Molecular phylogeny using the mitochondrial genome confirmed that P. hypophthalmus is a member of the family Pangasiidae and is nested within a clade including the families Cranoglanididae and Ictaluridae. The nuclear genome was estimated at approximately 700 Mb, assembled into 568 scaffolds with an N50 of 14.29 Mbp, and was estimated to contain ~ 28,600 protein-coding genes, comparable to those of channel catfish and zebrafish. Interestingly, zebrafish produce gadusol, but genes for biosynthesis of this sunscreen compound have been lost from catfish genomes. The differences in gene contents between these two catfishes were found in genes for vitamin D-binding protein and cytosolic phospholipase A₂, which have lost only in channel catfish. The Hox cluster in catfish genomes comprised seven paralogous groups, similar to that of zebrafish, and comparative analysis clarified catfish lineage-specific losses of A5a, B10a, and A11a. Genes for insulin-like growth factor (IGF) signaling were conserved between the two catfish genomes. In addition to identification of MHC class I and sex determination-related gene loci, the hypothetical chromosomes by comparison with the channel catfish demonstrated the usefulness of the striped catfish genome as a marker resource.

Conclusions

We developed genomic resources for the striped catfish. Possible conservation of genes for development and marker candidates were confirmed by comparing the assembled genome to that of a model fish, Danio rerio, and to channel catfish. Since the catfish genomic constituent resembles that of zebrafish, it is likely that zebrafish data for gene functions is applicable to striped catfish as well.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5079-x) contains supplementary material, which is available to authorized users.

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.

Exploiting the extraordinary genetic polymorphism of ciona for developmental genetics with whole genome sequencing

Studies in tunicates such as Ciona have revealed new insights into the evolutionary origins of chordate development. Ciona populations are characterized by high levels of natural genetic variation, between 1 and 5%. This variation has provided abundant material for forward genetic studies. In the current study, we make use of deep sequencing and homozygosity mapping to map spontaneous mutations in outbred populations. With this method we have mapped two spontaneous developmental mutants. In Ciona intestinalis we mapped a short-tail mutation with strong phenotypic similarity to a previously identified mutant in the related species Ciona savignyi. Our bioinformatic approach mapped the mutation to a narrow interval containing a single mutated gene, α-laminin3,4,5, which is the gene previously implicated in C. savignyi. In addition, we mapped a novel genetic mutation disrupting neural tube closure in C. savignyi to a T-type Ca(2+) channel gene. The high efficiency and unprecedented mapping resolution of our study is a powerful advantage for developmental genetics in Ciona, and may find application in other outbred species.

Genetic and genomic toolbox of the chordate Ciona intestinalis

The experimental malleability and unique phylogenetic position of the sea squirt Ciona intestinalis as part of the sister group to the vertebrates have helped establish these marine chordates as model organisms for the study of developmental genetics and evolution. Here we summarize the tools, techniques, and resources available to the Ciona geneticist, citing examples of studies that employed such strategies in the elucidation of gene function in Ciona. Genetic screens, germline transgenesis, electroporation of plasmid DNA, and microinjection of morpholinos are all routinely employed, and in the near future we expect these to be complemented by targeted mutagenesis, homologous recombination, and RNAi. The genomic resources available will continue to support the design and interpretation of genetic experiments and allow for increasingly sophisticated approaches on a high-throughput, whole-genome scale.

Transcriptional enhancers in ascidian development

The study of cis-regulatory DNAs that control developmental gene expression is integral to the modeling of comprehensive genomic regulatory networks for embryogenesis. Ascidian embryos provide a unique opportunity for the analysis of cis-regulatory DNAs with cellular resolution in the context of a simple but typical chordate body plan. Here, we review landmark studies that have laid the foundations for the study of transcriptional enhancers, among other cis-regulatory DNAs, and their roles in ascidian development. The studies using ascidians of the Ciona genus have capitalized on a unique electroporation technique that permits the simultaneous transfection of hundreds of fertilized eggs, which develop rapidly and express transgenes with little mosaicism. Current studies using the ascidian embryo benefit from extensively annotated genomic resources to characterize transcript models in silico. The search for functional noncoding sequences can be guided by bioinformatic analyses combining evolutionary conservation, gene coexpression, and combinations of overrepresented short-sequence motifs. The power of the transient transfection assays has allowed thorough dissection of numerous cis-regulatory modules, which provided insights into the functional constraints that shape enhancer architecture and diversification. Future studies will benefit from pioneering stable transgenic lines and the analysis of chromatin states. Whole genome expression, functional and DNA binding data are being integrated into comprehensive genomic regulatory network models of early ascidian cell specification with a single-cell resolution that is unique among chordate model systems.

Evolutionary crossroads in developmental biology: the tunicates

The tunicates, or urochordates, constitute a large group of marine animals whose recent common ancestry with vertebrates is reflected in the tadpole-like larvae of most tunicates. Their diversity and key phylogenetic position are enhanced, from a research viewpoint, by anatomically simple and transparent embryos, compact rapidly evolving genomes, and the availability of powerful experimental and computational tools with which to study these organisms. Tunicates are thus a powerful system for exploring chordate evolution and how extreme variation in genome sequence and gene regulatory network architecture is compatible with the preservation of an ancestral chordate body plan.

Integration of the genetic map and genome assembly of fugu facilitates insights into distinct features of genome evolution in teleosts and mammals

The compact genome of fugu (Takifugu rubripes) has been used widely as a reference genome for understanding the evolution of vertebrate genomes. However, the fragmented nature of the fugu genome assembly has restricted its use for comparisons of genome architecture in vertebrates. To extend the contiguity of the assembly to the chromosomal level, we have generated a comprehensive genetic map of fugu and anchored the scaffolds of the assembly to the 22 chromosomes of fugu. The map consists of 1,220 microsatellite markers that provide anchor points to 697 scaffolds covering 86% of the genome assembly (http://www.fugu-sg.org/). The integrated genome map revealed a higher recombination rate in fugu compared with other vertebrates and a wide variation in the recombination rate between sexes and across chromosomes of fugu. We used the extended assembly to explore recent rearrangement events in the lineages of fugu, Tetraodon, and medaka and compared them with rearrangements in three mammalian (human, mouse, and opossum) lineages. Between the two pufferfishes, fugu has experienced fewer chromosomal rearrangements than Tetraodon. The gene order is more highly conserved in the three teleosts than in mammals largely due to a lower rate of interchromosomal rearrangements in the teleosts. These results provide new insights into the distinct patterns of genome evolution between teleosts and mammals. The consolidated genome map and the genetic map of fugu are valuable resources for comparative genomics of vertebrates and for elucidating the genetic basis of the phenotypic diversity of ~25 species of Takifugu that evolved within the last 5 My.

Embryological methods in ascidians: the Villefranche-sur-Mer protocols

Ascidians (marine invertebrates: urochordates) are thought to be the closest sister groups of vertebrates. They are particularly attractive models because of their non-duplicated genome and the fast and synchronous development of large populations of eggs into simple tadpoles made of about 3,000 cells. As a result of stereotyped asymmetric cleavage patterns all blastomeres become fate restricted between the 16- and 110 cell stage through inheritance of maternal determinants and/or cellular interactions. These advantageous features have allowed advances in our understanding of the nature and role of maternal determinants, inductive interactions, and gene networks that are involved in cell lineage specification and differentiation of embryonic tissues. Ascidians have also contributed to our understanding of fertilization, cell cycle control, self-recognition, metamorphosis, and regeneration. In this chapter we provide basic protocols routinely used at the marine station in Villefranche-sur-Mer using the cosmopolitan species of reference Ciona intestinalis and the European species Phallusia mammillata. These two models present complementary advantages with regard to molecular, functional, and imaging approaches. We describe techniques for basic culture of embryos, micro-injection, in vivo labelling, micro-manipulations, fixation, and immuno-labelling. These methods allow analysis of calcium signals, reorganizations of cytoplasmic and cortical domains, meiotic and mitotic cell cycle and cleavages as well as the roles of specific genes and cellular interactions. Ascidians eggs and embryos are also an ideal material to isolate cortical fragments and to isolate and re-associate individual blastomeres. We detail the experimental manipulations which we have used to understand the structure and role of the egg cortex and of specific blastomeres during development.

Ciona genetics

Ascidians, such as Ciona, are invertebrate chordates with simple embryonic body plans and small, relatively non-redundant genomes. Ciona genetics is in its infancy compared to many other model systems, but it provides a powerful method for studying this important vertebrate outgroup. Here we give basic methods for genetic analysis of Ciona, including protocols for controlled crosses both by natural spawning and by the surgical isolation of gametes; the identification and propagation of mutant lines; and strategies for positional cloning.

doublesex/mab3 related-1 (dmrt1) is essential for development of anterior neural plate derivatives in Ciona

Ascidian larvae have a hollow, dorsal central nervous system that shares many morphological features with vertebrate nervous systems yet is composed of very few cells. We show here that a null mutation in the gene dmrt1 in the ascidian Ciona savignyi results in profound abnormalities in the development of the sensory vesicle (brain), as well as other anterior ectodermal derivatives, including the palps and oral siphon primordium (OSP). Although the phenotype of the mutant embryos is variable, the majority have a complete loss of the most anterior structures (palps and OSP) and extensive disruption of sensory structures, such as the light-sensitive ocellus, in the sensory vesicle. dmrt1 is expressed early in the blastula embryo in a small group of presumptive ectodermal cells as they become restricted to anterior neural, OSP and palp fates. Despite the early and restricted expression of dmrt1, we were unable, using several independent criteria, to observe a defect in the mutant embryos until the early tailbud stage. We speculate that the variability and late onset in the phenotype may be due to partially overlapping activities of other gene products.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a “covarion-like” pattern of correlated changes in evolutionary rates, and switches in conserved residues (“conserved but different”). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than “constant but different,” but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Computational analysis of Ciona intestinalis operons

Operons are clusters of genes that are co-regulated from a common promoter. Operons are typically associated with prokaryotes, although a small number of eukaryotes have been shown to possess them. Among metazoans, operons have been extensively characterized in the nematode Caenorhabditis elegans in which ∼15% of the total genes are organized into operons. The most recent genome assembly for the ascidian Ciona intestinalis placed ∼20% of the genes (2909 total) into 1310 operons. The majority of these operons are composed of two genes, while the largest are composed of six. Here is reported a computational analysis of the genes that comprise the Ciona operons. Gene ontology (GO) terms were identified for about two-thirds of the operon-encoded genes. Using the extensive collection of public EST libraries, estimates of temporal patterns of gene expression were generated for the operon-encoded genes. Lastly, conservation of operons was analyzed by determining how many operon-encoded genes were present in the ascidian Ciona savignyi and whether these genes were organized in orthologous operons. Over 68% of the operon-encoded genes could be assigned one or more GO terms and 697 of the 1310 operons contained genes in which all genes had at least one GO term. Of these 697 operons, GO terms were shared by all of the genes within 146 individual operons, suggesting that most operons encode genes with unrelated functions. An analysis of operon gene expression from nine different EST libraries indicated that for 587 operons, all of the genes that comprise an individual operon were expressed together in at least one EST library, suggesting that these genes may be co-regulated. About 50% (74/146) of the operons with shared GO terms also showed evidence of gene co-regulation. Comparisons with the C. savignyi genome identified orthologs for 1907 of 2909 operon genes. About 38% (504/1310) of the operons are conserved between the two Ciona species. These results suggest that like C. elegans, operons in Ciona are comprised of a variety of genes that are not necessarily related in function. The genes in only 50% of the operons appear to be co-regulated, suggesting that more complex gene regulatory mechanisms are likely operating.

miRTRAP, a computational method for the systematic identification of miRNAs from high throughput sequencing data

A novel method for prediction of miRs from deep sequencing data. Its utility is demonstrated when applied to Ciona data.

MicroRNAs (miRs) have been broadly implicated in animal development and disease. We developed a novel computational strategy for the systematic, whole-genome identification of miRs from high throughput sequencing information. This method, miRTRAP, incorporates the mechanisms of miR biogenesis and includes additional criteria regarding the prevalence and quality of small RNAs arising from the antisense strand and neighboring loci. This program was applied to the simple chordate Ciona intestinalis and identified nearly 400 putative miR loci.

Characterization of a corticotropin-releasing factor (CRF)/diuretic hormone-like peptide from tunicates: insight into the origins of the vertebrate CRF family

The corticotropin-releasing factor (CRF) superfamily of peptides includes the four paralogous vertebrate peptide lineages of CRF, urotensin-1/urocortin/sauvagine, urocortin 2 and urocortin 3, as well as the diuretic hormones (DH) and peptides of the arthropods. However, there are considerable sequence differences between the group of vertebrate peptides and those of the arthropods, notably insects. Because of the likely incidence of the formation of paralogous forms within the insects and the great variation in primary structures among these peptides, establishing the structure of the ancestral vertebrate version has not been possible. We screened the genomes of the tunicates, Ciona intestinalis and Ciona savignii, in silico, using the various conserved motifs of both the vertebrate CRF paralogues and the insect diuretic hormone sequences to identify the structure of the Ciona CRF/DH-like peptide genes. A single peptide gene was found in both genomes that possessed motifs reflective of both groups of peptides. These structures suggest a single CRF-like peptide was inherited by vertebrates and possibly chordates. Moreover, the conserved structure of the CRF peptide may have become constrained once it became associated with the regulation of the hypothalamus-pituitary-adrenal/interrenal axis.

Construction of a high-resolution genetic linkage map and comparative genome analysis for the reef-building coral Acropora millepora

A high-resolution genetic linkage map for the coral Acropora millepora is constructed and compared with other metazoan genomes, revealing syntenic blocks.

Background

Worldwide, coral reefs are in decline due to a range of anthropogenic disturbances, and are now also under threat from global climate change. Virtually nothing is currently known about the genetic factors that might determine whether corals adapt to the changing climate or continue to decline. Quantitative genetics studies aiming to identify the adaptively important genomic loci will require a high-resolution genetic linkage map. The phylogenetic position of corals also suggests important applications for a coral genetic map in studies of ancestral metazoan genome architecture.

Results

We constructed a high-resolution genetic linkage map for the reef-building coral Acropora millepora, the first genetic map reported for any coral, or any non-Bilaterian animal. More than 500 single nucleotide polymorphism (SNP) markers were developed, most of which are transferable in populations from Orpheus Island and Great Keppel Island. The map contains 429 markers (393 gene-based SNPs and 36 microsatellites) distributed in 14 linkage groups, and spans 1,493 cM with an average marker interval of 3.4 cM. Sex differences in recombination were observed in a few linkage groups, which may be caused by haploid selection. Comparison of the coral map with other metazoan genomes (human, nematode, fly, anemone and placozoan) revealed synteny regions.

Conclusions

Our study develops a framework that will be essential for future studies of adaptation in coral and it also provides an important resource for future genome sequence assembly and for comparative genomics studies on the evolution of metazoan genome structure.