Origin and Dispersal History of Two Colonial Ascidian Clades in the Botryllus schlosseri Species Complex

Historical navigation routes in European waters leave their footprint on the contemporary seascape genetics of a colonial urochordate

Humans have intensively sailed the Mediterranean and European Atlantic waters throughout history, from the upper Paleolithic until today and centuries of human seafaring have established complex coastal and cross-seas navigation networks. Historical literature revealed three major long-lasting maritime routes (eastern, western, northern) with four commencing locations (Alexandria, Venice, Genoa, Gibraltar) and a fourth route (circum-Italian) that connected between them. Due to oceangoing and technological constraints, most voyages were coastal, lasted weeks to months, with extended resting periods, allowing the development of fouling organisms on ship hulls. One of the abiding travellers in maritime routes is the colonial ascidian Botryllus schlosseri already known since the eighteenth century in European and Mediterranean ports. This species, was almost certainly one of the common hull fouling travellers in all trade routes for centuries. Employing COI haplotypes (1008 samples) and microsatellite alleles (995 samples) on colonies sampled from 64 pan-European sites, present-day Botryllus populations in the Mediterranean Sea/European Atlantic revealed significant segregation between all four maritime routes with a conspicuous partition of the northern route. These results reveal that past anthropogenic transports of sedentary marine species throughout millennia long seafaring have left their footprint on contemporary seascape genetics of marine organisms.

Phylogeographic and Morphological Analysis of Botrylloides niger Herdman, 1886 from the Northeastern Mediterranean Sea

Botrylloides niger (class Ascidiacea) is an invasive marine filter-feeding invertebrate that is believed to originate from the West Atlantic region. This species of colonial tunicate has been observed in several locations along the coasts of Israel and around the Suez Canal, but it has not yet been reported on the coasts of the Northeastern Mediterranean Sea (NEMS), suggesting an ongoing Lessepsian migration. However, the extent of this invasion might be concealed by reports of other potentially misidentified species of Botrylloides, given that the strong morphological similarities within this genus renders taxonomical identification particularly challenging. In this study, we performed a phylogeographic and morphological analysis of B. niger in the NEMS. We collected 238 samples from 8 sampling stations covering 824 km of the coastlines of NEMS. We reported 14 different morphotypes, of which the orange-brown, orange, and brown-striped morphs were the most abundant. Using the mitochondrial cytochrome C oxidase I (COI) as a DNA barcode marker, we identified 4 haplotypes. The COI haplotypes clustered with the reference B. niger sequences from GenBank and differed significantly from the sister Botrylloides species. We confirmed our identification using three additional barcoding markers (Histone 3, 18S rRNA, and 28S rRNA), which all matched with over 99% similarity to reference sequences. In addition, we monitored a station for a year and conducted a temporal analysis of the collected colonies. The colonies were absent during the winter and spring, while new colonies were established in the summer and expanded during autumn. We performed demographic population analysis on our spatial data that identified a possible population subdivision at a sampling site, which might have been caused by local freshwater input. Herein, we present the first report on the presence of Botrylloides niger in the NEMS. This study represents a key step toward understanding the diversity and the propagation of this highly invasive species of colonial ascidians, both within the Mediterranean basin as well as globally.

Complementary genomic and epigenomic adaptation to environmental heterogeneity

While adaptation is commonly thought to result from selection on DNA sequence-based variation, recent studies have highlighted an analogous epigenetic component as well. However, the relative roles of these mechanisms in facilitating population persistence under environmental heterogeneity remain unclear. To address the underlying genetic and epigenetic mechanisms and their relationship during environmental adaptation, we screened the genomes and epigenomes of nine global populations of a predominately sessile marine invasive tunicate, Botryllus schlosseri, using reduced-representation methods. We detected clear population differentiation at the genetic and epigenetic levels. Patterns of genetic and epigenetic structure were significantly influenced by local environmental variables. Among these variables, minimum annual sea surface temperature was identified as the top explanatory variable for both genetic and epigenetic variation. However, patterns of population structure driven by genetic and epigenetic variation were somewhat distinct, suggesting possible autonomy of epigenetic variation. We found both shared and specific genes and biological pathways among genetic and epigenetic loci associated with environmental factors, consistent with complementary and independent contributions of genetic and epigenetic variation to environmental adaptation in this system. Collectively, these mechanisms may facilitate population persistence under environmental change and sustain successful invasions across novel environments.

First Record of Colonial Ascidian, Botrylloides diegensis Ritter and Forsyth, 1917 (Ascidiacea, Stolidobranchia, Styelidae), in South Korea

Botrylloides species are important members of the fouling community colonizing artificial substrates in harbors and marinas. During monitoring in 2017–2020 of non-indigenous species in Korea, one colonial ascidian species was distinctly different from other native colonial ascidians, such as B. violaceus and Botryllus schlosseri, in South Korea. This species was identified as B. diegensis. DNA barcodes with mitochondrial COI were used to identify one-toned and two-toned colonies of B. diegensis. Intraspecific variations between Korean and other regions of B. diegensis from the NCBI ranged from 0.0% to 1.3%. The Korean B. diegensis was clearly distinct from other species of Botrylloides at 15.8–24.2%. In phylogenetic analysis results, Korean B. diegensis was established as a single clade with other regions of B. diegensis and was clearly distinct from Korean B. violaceus. After reviewing previous monitoring data, it was found that two-toned B. diegensis was already found in six harbors by July 2017. It has now spread into 14 harbors along the coastal line of South Korea. This means that B. diegensis might have been introduced to South Korea between 1999 and 2016.

Phylogenomic and morphological relationships among the botryllid ascidians (Subphylum Tunicata, Class Ascidiacea, Family Styelidae)

Ascidians (Phylum Chordata, Class Ascidiacea) are a large group of invertebrates which occupy a central role in the ecology of marine benthic communities. Many ascidian species have become successfully introduced around the world via anthropogenic vectors. The botryllid ascidians (Order Stolidobranchia, Family Styelidae) are a group of 53 colonial species, several of which are widespread throughout temperate or tropical and subtropical waters. However, the systematics and biology of this group of ascidians is not well-understood. To provide a systematic framework for this group, we have constructed a well-resolved phylogenomic tree using 200 novel loci and 55 specimens. A Principal Components Analysis of all species described in the literature using 31 taxonomic characteristics revealed that some species occupy a unique morphological space and can be easily identified using characteristics of adult colonies. For other species, additional information such as larval or life history characteristics may be required for taxonomic discrimination. Molecular barcodes are critical for guiding the delineation of morphologically similar species in this group.

An elongated COI fragment to discriminate botryllid species and as an improved ascidian DNA barcode

Botryllids are colonial ascidians widely studied for their potential invasiveness and as model organisms, however the morphological description and discrimination of these species is very problematic, leading to frequent specimen misidentifications. To facilitate species discrimination and detection of cryptic/new species, we developed new barcoding primers for the amplification of a COI fragment of about 860 bp (860-COI), which is an extension of the common Folmer's barcode region. Our 860-COI was successfully amplified in 177 worldwide-sampled botryllid colonies. Combined with morphological analyses, 860-COI allowed not only discriminating known species, but also identifying undescribed and cryptic species, resurrecting old species currently in synonymy, and proposing the assignment of clade D of the model organism Botryllus schlosseri to Botryllus renierii. Importantly, within clade A of B. schlosseri, 860-COI recognized at least two candidate species against only one recognized by the Folmer's fragment, underlining the need of further genetic investigations on this clade. This result also suggests that the 860-COI could have a greater ability to diagnose cryptic/new species than the Folmer's fragment at very short evolutionary distances, such as those observed within clade A. Finally, our new primers simplify the amplification of 860-COI even in non-botryllid ascidians, suggesting their wider usefulness in ascidians.

Genetic analyses reveal cryptic diversity in the widely distributed Styela canopus (Ascidiacea:Styelidae)

The routine use of DNA sequencing techniques and phylogenetic analysis has resulted in the discovery of many cryptic species, especially in the oceans. The common, globally introduced species Styela canopus is suspected to be a complex of cryptic species because of its widespread distribution and variable external morphology. We tested this possibility using COI and ANT marker sequences to uncover the phylogenetic relationship among 19 populations, and to examine genetic variability as well as gene flow. We obtained 271 COI and 67 ANT sequences and found surprising diversity among the 19 populations (COI: π = 0.18, hd = 0.99; ANT: π = 0.13, hd = 0.95). Corresponding topologies were found using Bayesian inference and maximum likelihood for both simple locus (COI) and multilocus (COI + ANT) analyses and so the clades received strong support. We used simple (ABGD, bPTP, GMYC) and multiple (BSD) locus methods to delimit species. The simple locus methods indicated that the current Styela canopus comprises at least 15 species. The BSD method for concatenated data supported 7 of the 15 species. We suggest that S. canopus should be treated as the Styela canopus complex. The large number of cryptic species found, often with more than one clade found in sympatry, creates opportunities for better understanding reproductive isolation, hybridisation or speciation. As several lineages have already been introduced widely around the world, we must quickly understand their diversity and invasive abilities.

Hidden diversity in Antarctica: Molecular and morphological evidence of two different species within one of the most conspicuous ascidian species

The Southern Ocean is one of the most isolated marine ecosystems, characterized by high levels of endemism, diversity, and biomass. Ascidians are among the dominant groups in Antarctic benthic assemblages; thus, recording the evolutionary patterns of this group is crucial to improve our current understanding of the assembly of this polar ocean. We studied the genetic variation within Cnemidocarpa verrucosa sensu lato, one of the most widely distributed abundant and studied ascidian species in Antarctica. Using a mitochondrial and a nuclear gene (COI and 18S), the phylogeography of fifteen populations distributed along the West Antarctic Peninsula and Burdwood Bank/MPA Namuncurá (South American shelf) was characterized, where the distribution of the genetic distance suggested the existence of, at least, two species within nominal C. verrucosa. When reevaluating morphological traits to distinguish between genetically defined species, the presence of a basal disk in one of the genotypes could be a diagnostic morphological trait to differentiate the species. These results are surprising due to the large research that has been carried out with the conspicuous C. verrucosa with no differentiation between species. Furthermore, it provides important tools to distinguish species in the field and laboratory. But also, these results give new insights into patterns of differentiation between closely related species that are distributed in sympatry, where the permeability of species boundaries still needs to be well understood.

Impacts of climate change on geographical distributions of invasive ascidians

Ocean warming associated with global climate change renders marine ecosystems susceptible to biological invasions. Here, we used species distribution models to project habitat suitability for eight invasive ascidians under present-day and future climate scenarios. Distance to shore and maximum sea surface temperature were identified as the most important variables affecting species distributions. Results showed that eight ascidians might respond differently to future climate change. Alarmingly, currently colonized areas are much smaller than predicted, suggesting ascidians may expand their invasive ranges. Areas such as Americas, Europe and Western Pacific have high risks of receiving new invasions. In contrast, African coasts, excluding the Mediterranean side, are not prone to new invasions, likely due to the high sea surface temperature there. Our results highlight the importance of climate change impacts on future invasions and the need for accurate modelling of invasion risks, which can be used as guides to develop management strategies.

Toward a resolution of the cosmopolitan Botryllus schlosseri species complex (Ascidiacea, Styelidae): mitogenomics and morphology of clade E (Botryllus gaiae)

Botryllus schlosseri is a model colonial ascidian and a marine invader. It is currently recognized as a species complex comprising five genetically divergent clades, with clade A globally distributed and clade E found only in Europe. This taxon has also been recently redescribed by designation of a clade A specimen as the neotype. To clarify the taxonomic status of clade E and its relationship to clade A, we examine the entire mitochondrial genome and study the morphology of clade E. The mitogenome of clade E has an identical gene order to clade A, but substantially differs in the size of several non-coding regions. Remarkably, the nucleotide divergence of clade A-clade E is incompatible with the intraspecies ascidian divergence, but similar to the congeneric one and almost identical to the divergence between species once considered morphologically indistinguishable (e.g. the pair Ciona intestinalis (Linnaeus, 1767)-Ciona robusta Hoshino & Tokioka, 1967, and the pair Botrylloides niger Herdman, 1886-Botrylloides leachii (Savigny, 1816)). Clade E differs morphologically from the Botryllusschlosseri neotype mainly in the number and appearance of the stomach folds, and the shape of the anal opening, the first intestinal loop and the typhlosole. Our integrative taxonomical approach clearly distinguishes clade E as a species separate from Botryllusschlosseri, with unique morphological and molecular characters. Therefore, we here describe clade E as the new species Botryllus gaiae sp. nov.

Comparative feeding rates of native and invasive ascidians

Ascidians have a recent history of species introductions globally, often with strong ecological impacts. Comparisons of per capita effects of invaders and comparable natives are useful to assess such impacts. Here, we explore ingestion rates (IR) and clearance rates (CR) of Ciona intestinalis and Ciona robusta, co-occurring native and non-native ascidians, respectively, from Brittany, France. IR was positively related to food concentration, with the invader responding more strongly to increasing food concentration. CR also differed by species, with the invader demonstrating higher values. C. robusta exhibited a higher functional response (Type I) than did C. intestinalis (Type II). Relative impact measured using seasonal abundance and IR revealed that C. robusta has a much greater impact than C. intestinalis at all food concentrations tested, though the former has a constrained distribution which limits its regional impact. Nevertheless, when abundant, we expect C. robusta to exert a greater impact on algal foods.

De novo draft assembly of the Botrylloides leachii genome provides further insight into tunicate evolution

Tunicates are marine invertebrates that compose the closest phylogenetic group to the vertebrates. These chordates present a particularly diverse range of regenerative abilities and life-history strategies. Consequently, tunicates provide an extraordinary perspective into the emergence and diversity of these traits. Here we describe the genome sequencing, annotation and analysis of the Stolidobranchian Botrylloides leachii. We have produced a high-quality 159 Mb assembly, 82% of the predicted 194 Mb genome. Analysing genome size, gene number, repetitive elements, orthologs clustering and gene ontology terms show that B. leachii has a genomic architecture similar to that of most solitary tunicates, while other recently sequenced colonial ascidians have undergone genome expansion. In addition, ortholog clustering has identified groups of candidate genes for the study of colonialism and whole-body regeneration. By analysing the structure and composition of conserved gene linkages, we observed examples of cluster breaks and gene dispersions, suggesting that several lineage-specific genome rearrangements occurred during tunicate evolution. We also found lineage-specific gene gain and loss within conserved cell-signalling pathways. Such examples of genetic changes within conserved cell-signalling pathways commonly associated with regeneration and development that may underlie some of the diverse regenerative abilities observed in tunicates. Overall, these results provide a novel resource for the study of tunicates and of colonial ascidians.

A massive update of non-indigenous species records in Mediterranean marinas

The Mediterranean Sea is home to over 2/3 of the world's charter boat traffic and hosts an estimated 1.5 million recreational boats. Studies elsewhere have demonstrated marinas as important hubs for the stepping-stone transfer of non-indigenous species (NIS), but these unique anthropogenic, and typically artificial habitats have largely gone overlooked in the Mediterranean as sources of NIS hot-spots. From April 2015 to November 2016, 34 marinas were sampled across the following Mediterranean countries: Spain, France, Italy, Malta, Greece, Turkey and Cyprus to investigate the NIS presence and richness in the specialized hard substrate material of these marina habitats. All macroinvertebrate taxa were collected and identified. Additionally, fouling samples were collected from approximately 600 boat-hulls from 25 of these marinas to determine if boats host diverse NIS not present in the marina. Here, we present data revealing that Mediterranean marinas indeed act as major hubs for the transfer of marine NIS, and we also provide evidence that recreational boats act as effective vectors of spread. From this wide-ranging geographical study, we report here numerous new NIS records at the basin, subregional, country and locality level. At the basin level, we report three NIS new to the Mediterranean Sea (Achelia sawayai sensu lato, Aorides longimerus, Cymodoce aff. fuscina), and the re-appearance of two NIS previously known but currently considered extinct in the Mediterranean (Bemlos leptocheirus, Saccostrea glomerata). We also compellingly update the distributions of many NIS in the Mediterranean Sea showing some recent spreading; we provide details for 11 new subregional records for NIS (Watersipora arcuata, Hydroides brachyacantha sensu lato and Saccostrea glomerata now present in the Western Mediterranean; Symplegma brakenhielmi, Stenothoe georgiana, Spirobranchus tertaceros sensu lato, Dendostrea folium sensu lato and Parasmittina egyptiaca now present in the Central Mediterranean, and W. arcuata, Bemlos leptocheirus and Dyspanopeus sayi in the Eastern Mediterranean). We also report 51 new NIS country records from recreational marinas: 12 for Malta, 10 for Cyprus, nine for Greece, six for Spain and France, five for Turkey and three for Italy, representing 32 species. Finally, we report 20 new NIS records (representing 17 species) found on recreational boat-hulls (mobile habitats), not yet found in the same marina, or in most cases, even the country. For each new NIS record, their native origin and global and Mediterranean distributions are provided, along with details of the new record. Additionally, taxonomic characters used for identification and photos of the specimens are also provided. These new NIS records should now be added to the relevant NIS databases compiled by several entities. Records of uncertain identity are also discussed, to assess the probability of valid non-indigenous status.