Genomic resources of Ascidiella aspersa and comparative analysis across tunicates reveal class-level features and evolutionary diversification

This study establishes high-quality genomic and transcriptomic resources for the invasive tunicate *Ascidiella aspersa* and conducts a comprehensive comparative analysis across 35 other tunicate species to reveal class-level genomic features, evolutionary diversification patterns, and a new phylogenetic hypothesis, all accessible through the newly constructed TUNOME database.

The Gist

Imagine the ocean is a vast library, and inside it lives a tiny, transparent creature called the sea squirt (Ascidiella aspersa). For a long time, scientists knew this creature was special because it's a distant cousin to humans (we are both chordates), and its babies are so see-through you can watch their organs form in real-time. But there was a problem: we didn't have the "instruction manual" (the genome) for how to build them. It was like trying to fix a complex machine without the blueprint.

This paper is the story of a team of scientists who finally wrote that blueprint, and then used it to compare sea squirts with their entire extended family tree to see how they evolved.

Here is the breakdown of their adventure, explained simply:

1. Writing the Instruction Manual (The Genome Assembly)

Think of a genome as a massive cookbook containing every recipe needed to build a living thing. For A. aspersa, this cookbook was missing, or at least very messy.

• The Challenge: The scientists had to piece together millions of tiny DNA fragments (like a giant jigsaw puzzle) to create a complete, high-quality book.
• The Result: They successfully assembled a 306.5 million "letter" long cookbook. They didn't just stop at the text; they also added "highlighters" and "footnotes" (annotations) to explain what each gene does. They found about 24,500 genes.
• Why it matters: Now, instead of guessing how this transparent creature works, scientists have a clear map. This makes A. aspersa a new "super-model" for studying how animals develop, especially because its eggs are so clear.

2. The Family Reunion (Comparative Analysis)

Once they had the manual for A. aspersa, the team didn't stop there. They went to the library and pulled out the manuals for 37 other species of sea squirts and their relatives (tunicates).

• The Cleanup: Many of these other manuals were incomplete, written in a foreign language, or missing pages. The team used powerful computer tools to "translate" and "repair" these manuals, creating a standardized set of instructions for 38 different species.
• The Database (TUNOME): They built a giant, online "Google" for these sea creatures called TUNOME. Now, any scientist in the world can type in a gene name and instantly see how it looks in 38 different species, compare them, and download the data. It's like having a universal translator for the sea squirt family.

3. The Evolutionary Detective Work

With all these manuals in hand, the scientists started looking for patterns, like a detective solving a mystery.

• The Family Tree: They built a new family tree. It turns out the relationships between different groups of sea squirts are a bit more complicated than we thought. Some groups that looked like they were cousins might actually be more distant, and vice versa.
• The "Missing Tools" Mystery: They noticed something weird. Some groups of sea squirts (specifically the "salps" and "larvaceans") were missing entire chapters of their DNA repair manuals.
  • Analogy: Imagine a car factory that suddenly stops making the tools needed to fix a flat tire. Over time, the cars in that factory might start looking very different or driving in strange ways because they can't repair themselves the same way others can. The scientists think losing these "repair tools" might be why some sea squirt families evolved so rapidly and look so different from their ancestors.
• The "Accent" of DNA (GC Content): DNA is made of four letters: A, T, C, and G. Some species use a lot of G and C (like a heavy accent), while others use mostly A and T.
  • The scientists found that A. aspersa has a very "heavy" accent (high GC content) compared to its neighbors. This isn't just a random quirk; it changes how the genes are read and might be a signature of how this specific family evolved.

4. Why Should You Care?

You might wonder, "Why do we care about a transparent sea squirt?"

• The Human Connection: Because sea squirts are our distant cousins, studying them helps us understand how our own bodies (and our ancestors) developed.
• The Transparency Trick: Since A. aspersa embryos are see-through, they are perfect for watching development happen. Now that we have the "instruction manual," we can figure out exactly which genes make them transparent. This could lead to new discoveries in biology and medicine.
• The Toolkit: By fixing the manuals for 38 species, the scientists gave the whole scientific community a better toolbox. Instead of struggling with broken data, researchers can now focus on asking big questions about evolution, development, and how life adapts to the ocean.

In a nutshell: This paper is about taking a messy, missing instruction manual for a cool, transparent sea creature, fixing it, and then using it to clean up the manuals for its entire family. The result is a new, shared library that helps us understand how these animals evolved, why they look the way they do, and what secrets they hold about the history of life on Earth.

Technical Summary

1. Problem Statement

• Lack of Genomic Resources: Despite Ascidiella aspersa being a cosmopolitan invasive species with remarkably transparent embryos (making it a promising model for developmental biology), it lacks a high-quality, annotated genome assembly.
• Incomplete Comparative Landscape: While several tunicate genomes exist, many lack functional annotations or are of low quality. There has been no comprehensive comparative analysis across the three major tunicate classes (Ascidiacea, Thaliacea, and Appendicularia) to understand class-level genomic features, evolutionary diversification, and phylogenetic relationships.
• Need for Non-Model Resources: Current evo-devo studies are heavily reliant on a few model organisms (e.g., Ciona intestinalis, Halocynthia roretzi). Expanding genomic resources to non-model species is essential for broader evolutionary insights.

2. Methodology

The study employed a multi-faceted approach combining de novo assembly, comparative genomics, and database construction:

• Genome Assembly (A. aspersa):
  • Data: Used sperm (haploid) DNA to simplify assembly. Generated ~46.7 Gb of PacBio long reads and ~70 Gb of Illumina short reads.
  • Assembly: Employed hybrid assembly using tools like wtdbg2, Flye, Canu, and SMARTdenovo. Selected the wtdbg2 assembly for its superior contiguity (N50) and polished it with Illumina data.
  • Quality Control: Assessed completeness using BUSCO (Eukaryota dataset).
• Transcriptome Profiling:
  • Sequenced RNA from 9 adult organs and 6 developmental stages (embryos/larvae).
  • Assembled using Trinity and aligned to the genome using Hisat2.
• Gene Model Construction:
  • Pipeline: Used a combination of ab initio prediction (BRAKER2) and homology-based prediction (GeMoMa) using RNA-seq evidence and reference proteomes from related species.
  • Scope: Constructed or updated gene models for 38 tunicate species (including the new A. aspersa model) spanning 3 classes, 5 orders, and 12 families.
  • Filtering: Integrated models were filtered to retain the longest splice variants and remove genes with premature stop codons.
• Comparative & Phylogenomic Analysis:
  • Phylogeny: Used OrthoFinder to identify orthogroups. Constructed species trees using STAG (with paralogs) and IQ-TREE (maximum likelihood with single-copy orthologs) to resolve relationships between Phlebobranchia and Aplousobranchia.
  • Genomic Statistics: Analyzed genome size, GC content (whole genome vs. coding sequences), and codon usage bias (ENC-GC3 plots, PCA).
  • Gene Family Evolution: Used CAFÉ to identify rapidly expanded or contracted gene families and EnrichGO for functional enrichment.
• Database Construction: Developed TUNOME, a web-based database integrating genome sequences, gene models, functional annotations, and ortholog analysis for all 38 species.

3. Key Contributions

1. High-Quality A. aspersa Genome: Delivered a 306.5 Mb de novo genome assembly with 24,504 predicted genes and a BUSCO completeness score of 92.2%.
2. Pan-Tunicate Gene Models: Reconstructed high-quality gene models for 38 tunicate species, significantly improving BUSCO scores for 27 species (achieving >80% completeness).
3. TUNOME Database: Created a centralized, searchable resource allowing for BLAST searches, ortholog analysis, and expression visualization across diverse tunicate taxa.
4. Evolutionary Insights: Provided a new phylogenetic hypothesis regarding the relationships within Phlebobranchia and Aplousobranchia and identified class-specific genomic signatures.

4. Key Results

• Genome Assembly Statistics:

  • A. aspersa assembly: 306.5 Mb, 1,411 contigs, N50 of ~2.06 Mb.
  • Gene model (KAS25): 24,504 genes, 30,856 transcripts.
  • Functional annotation: 18,636 genes annotated (GO, KEGG, UniProt).

• Phylogenetic Relationships:

  • Confirmed the paraphyly of Ascidiacea (with Thaliacea nested within).
  • Proposed a new topology where Aplousobranchia is nested within Phlebobranchia.
  • Identified a sister-group relationship between Ascidiidae and Perophoridae, and between Corellidae and Cionidae.

• Genomic Diversity (GC Content & Size):

  • Genome Size: Varied 17-fold (64.3 Mb in Oikopleura dioica to 1.1 Gb in Trididemnum nubilum).
  • GC Content: Whole-genome GC ranged from 25.3% (Thalia democratica) to 42.7% (Polycarpa mamillaris).
  • Coding Sequence (CDS) GC: Generally 5–10% higher than whole-genome GC. A. aspersa exhibited the highest CDS GC content (47.1%) among all 38 species, despite having an intermediate whole-genome GC.
  • Codon Usage: A. aspersa showed distinct codon usage bias compared to other Ascidiidae (e.g., preference for TGC over TGT for Cysteine), likely driven by selection pressure rather than just mutational bias.

• Gene Family Evolution:

  • Thaliacea and Botrylloides: Showed rapid contraction of gene families, particularly those related to DNA repair (e.g., Base Excision Repair genes like SMUG1, MUTYH, MPG). This loss correlates with their high genomic variability and AT-richness (in some species).
  • A. aspersa: 51 gene families were rapidly expanded (enriched in DNA binding and transcriptional regulation), and 315 families were unique to the species.

• Functional Insights from TUNOME:

  • DNA Repair: Confirmed the loss of specific DNA repair pathways in Appendicularia and Thaliacea, suggesting a mechanism for their high genomic plasticity.
  • Egg Transparency: Identified 259 genes specific to Ascidiidae. Top highly expressed genes in A. aspersa eggs include those involved in cell cycle regulation, calcium signaling, and Nucleotide Excision Repair (NER), potentially protecting transparent eggs from UV damage.

5. Significance

• Model Organism Expansion: Establishes A. aspersa as a viable, genetically tractable model for developmental biology, offering an alternative to European-restricted models like Phallusia mammillata.
• Evolutionary Framework: The study provides a robust phylogenomic framework for Tunicata, resolving long-standing debates about the relationships between major clades and highlighting class-specific evolutionary trajectories (e.g., GC content shifts).
• Mechanistic Insights: Links genomic features (GC content, codon bias, DNA repair gene loss) to biological traits such as genome size variation, developmental speed, and environmental adaptation (e.g., UV protection in transparent eggs).
• Resource Availability: The TUNOME database democratizes access to tunicate genomics, enabling researchers to perform comparative analyses across non-model species, thereby accelerating evo-devo research beyond the traditional few model organisms.