Invasion genetics of the Ciona intestinalis species complex: from regional endemism to global homogeneity

Determining the degree of population connectivity and investigating factors driving genetic exchange at various geographical scales are essential to understanding population dynamics and spread potential of invasive species. Here, we explore these issues in the highly invasive vase tunicate, Ciona intestinalis, a species whose invasion history has been obscured by its poorly understood taxonomy and population genetics. Recent phylogenetic and comparative genomic studies suggest that C. intestinalis is a cryptic species complex consisting of at least three species. We reconstructed phylogenies based on both mitochondrial (cytochrome c oxidase subunit 3—NADH dehydrogenase subunit 1 region and NADH dehydrogenase subunit 4 gene) and nuclear (internal transcribed spacer 1) sequences, results of which support four major phylogroups corresponding to the previously reported spA, spB and Ciona spp. (spC) as well as an undescribed cryptic species (spD). While spC and spD remain restricted to their native ranges in the Mediterranean Sea and Black Sea, respectively, the highly invasive species (spA and spB) have disjunct global distributions. Despite extensive interspecific divergences, we identified low phylogeographical structure within these two invasive species. Haplotype network analyses revealed comparatively limited mutation steps among haplotypes within each species. Population genetic analyses based on two mtDNA fragments and eight unlinked microsatellites illustrated relatively low population differentiation and high population connectivity at both regional and continental scales in the two invasive species. Human‐mediated dispersal coupled with a high potential for natural dispersal is probably responsible for the observed genetic homogeneity.     

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²⁺ 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.     

Predation contributes to invasion resistance of benthic communities against the non-indigenous tunicate <Emphasis Type="Italic">Ciona intestinalis</Emphasis>

Marine anthropogenic structures offer novel niches for introduced species but their role in the subsequent invasion to natural habitats remains unknown. Upon arrival in new environments, invaders must overcome biotic resistance from native competitors and predators if they are to establish successfully in natural habitats. We tested the hypotheses that (1) artificial structures (e.g., suspended aquaculture installations) present a niche opportunity for invasive species by providing a refuge from native benthic predators, and (2) native predators in natural benthic habitats suppress successful colonization by invaders. A recruitment experiment showed that the ascidians Pyura chilensis (native) and Ciona intestinalis (invasive) could recruit to both suspended artificial structures and natural benthic habitats. Ciona, however, was only able to establish adult populations on artificial structures. In natural benthic habitats Ciona only recruited and grew in predator-exclusion cages, because without this protection predation prevented its establishment. In predation experiments, native invertebrate and fish predators removed all invasive ascidians (recruits and adults) in benthic habitats, which contrasted with the high adult survival of the native ascidian P. chilensis. The refuge from a number of benthic predators facilitates the establishment of large populations of invasive species on suspended structures. We present a conceptual model of the invasion processes that includes the anthropogenic structures as a transitional stepping-stone that facilitates invasion by enhancing and prolonging propagule supply to surrounding natural communities. Those established invaders might then overcome biotic resistance during time periods when populations of consumers or competitors are weakened by natural or anthropogenic disturbances. Our results suggest that the conservation of natural habitats with a high diversity of native predators can be an effective means to prevent the spread of invasive species growing on suspended structures.     

Biofouling ascidians on aquaculture gear as potential vectors of harmful algal introductions

Biofouling ascidians are ubiquitous in coastal ecosystems and are among the main colonizers of aquaculture gear. Our study tested the hypothesis that the transport, removal, and transfer of fouling ascidian species by aquaculturists provide a mechanism for concentration and distribution of harmful-algal cells to new areas. Wild-caught specimens of common, biofouling ascidian species (Styela clava, Ciona intestinalis, Molgula manhattensis, Botrylloides violaceus, Didemnum vexillum, and Botryllus schlosseri) were exposed individually to cultured strains of co-occurring harmful algae (Prorocentrum minimum, Alexandrium fundyense, Alexandrium monilatum, Karenia brevis, Aureococcus anophagefferens, or Heterosigma akashiwo) at simulated bloom cell densities of each HAB species. After feeding, ascidians were transferred to ultrafiltered seawater. Immediately after exposure, and after 24 and 48 h in ultrafiltered seawater, biodeposits were collected and observed microscopically for the presence of intact, potentially viable cells. Subsamples of biodeposits were transferred into culture tubes with ultrafiltered seawater and monitored for algal growth during 8 weeks. Cells of all HAB species were found to pass intact through the ascidian digestive system, remained viable, and in many cases were capable of re-establishing populations at least 48 h post-ingestion. The results of our study will inform industry and managers of the potential threat and ecological impact of spreading biofouling ascidians, and practices to mitigate adverse impacts. Additionally, these management practices have been formally incorporated into a new cost-share program developed to help shellfish producers prevent the further spread of ascidians and associated HAB species.     

Characterization of highly polymorphic nuclear microsatellite loci from the ascidian Ciona intestinalis

Eight nuclear polymorphic microsatellite markers were characterized from the ascidian Ciona intestinalis whole genome sequence. The behaviour of these loci was investigated against two geographically distinct populations: one from Plymouth, UK the other from the Fusaro Lagoon, Italy, both belonging to the type A Ciona cryptic species. The markers exhibited six to 29 alleles and average observed heterozygosity ranging from 0.06 to 0.73. These new microsatellite loci demonstrated to be valuable tools for both population genetic analysis at different scales and genetic identification of mutant phenotypes frequently encountered in Mediterranean populations of C. intestinalis.     

污损性海鞘的生态特点研究展望

海鞘生长快,繁殖迅速,能产生大量在短时间内附着的幼虫,是海洋污损生物群落中的重要成员,对海上人工设施会产生严重危害。污损性海鞘主要由悉尼海鞘(Ascidia sydneiensis)、史氏菊海鞘(Botryllus schlosseri)、米氏小叶鞘(Diplosoma listerianum)、柄瘤海鞘(Styela clava)、红贺海鞘(Herdmania momus)等9科29属103种组成,其中在太平洋海域64种、印度洋23种、大西洋44种,而北冰洋海域仅3种;另外,其附着污损具有明显的地域性和季节特点,并与深度有关。今后应加强污损性海鞘的生态调查和分类研究,阐明深海及两极海域附着污损特点,揭示幼虫附着变态过程的分子调控机理,完善幼虫采集培养技术,以期更好地掌握海鞘生物学特性与生态特点,丰富和发展海洋生态学内容,并为海洋污损生物的防除奠定基础,促进海洋经济产业的发展。     

The proteomic response of sea squirts (genus Ciona) to acute heat stress: a global perspective on the thermal stability of proteins

Congeners belonging to the genus Ciona have disparate distributions limited by temperature. Ciona intestinalis is more widespread with a cosmopolitan distribution ranging from tropical to sub-arctic zones, while Ciona savignyi is limited to temperate-latitudes of the northern Pacific Ocean. To compare the heat stress response between congeners, we quantified changes in protein expression using proteomics. Animals were exposed to 22 °C, 25 °C, and 28 °C for 6 h, then recovered at a control temperature (13 °C) for 16 h (high heat stress experiment). In a second experiment we exposed animals to lower levels of heat stress at 18 °C, 20 °C, and 23 °C, with a 16 °C control. A quantitative analysis, using 2D gel electrophoresis and MALDI-TOF/TOF mass spectrometry (with a 69% and 93% identification rate for Ciona intestinalis and Ciona savignyi, respectively), showed changes in a number of protein functional groups, including molecular chaperones, extracellular matrix proteins, calcium-binding proteins, cytoskeletal proteins and proteins involved in energy metabolism. Our results indicate that C. intestinalis maintains higher constitutive levels of molecular chaperones than C. savignyi, suggesting that it is prepared to respond faster to thermal stress. Systematic discrepancies between estimated versus predicted molecular masses of identified proteins differed between protein families and were more pronounced under high heat conditions, suggesting that thermal sensitivities are lower for cytoskeletal proteins and ATP-synthase than for any other protein group represented on 2D gels.     

Mitochondrial Genome of <Emphasis Type="Italic">Ciona savignyi</Emphasis> (Urochordata,Ascidiacea, Enterogona): Comparison of Gene Arrangement and tRNA Genes with <Emphasis Type="Italic">Halocynthia roretzi</Emphasis> Mitochondrial Genome

The complete nucleotide sequence of the urochordate Ciona savignyi (Ascidiacea, Enterogona) mitochondrial (mt) genome (14,737 bp) was determined. The Ciona mt genome does not encode a gene for ATP synthetase subunit 8 but encodes an additional tRNA^Gly gene (anticodon UCU), as is the case in another urochordate, Halocynthia roretzi (Ascidiacea, Pleurogona), mt genome. In addition, the Ciona mt genome encodes two tRNA^Met genes; anticodon CAT and anticodon TAT. The tRNA^Cys gene is thought to lack base pairs at the D-stem. Thus, the Ciona mt genome encodes 12 protein, 2 rRNA, and 24 tRNA genes. The gene arrangement of the Ciona mt genome differs greatly from those of any other metazoan mt genomes reported to date. Only three gene boundaries are shared between the Halocynthia and the Ciona mt genomes. Molecular phylogenetic analyses based on amino acid sequences of mt protein genes failed to demonstrate the monophyly of the chordates.     

Regional connectivity and coastal expansion: differentiating pre‐border and post‐border vectors for the invasive tunicate Styela clava

The dramatic increase in marine bio‐invasions, particularly of non‐indigenous ascidians, has highlighted the vulnerability of marine ecosystems and the productive sectors that rely on them. A critical issue in managing invasive species is determining the relative roles of ongoing introductions, versus the local movement of propagules from established source populations. Styela clava (Herdman, 1882), the Asian clubbed tunicate, once restricted to the Pacific shores of Asia and Russia, is now abundant throughout the northern and southern hemispheres and has had significant economic impact in at least one site of incursion. In 2005 S. clava was identified in New Zealand. The recent introduction of this species, coupled with its restricted distribution, provided an ideal model to compare and contrast the introduction and expansion process. In this study, the mitochondrial DNA cytochrome oxidase subunit I gene (COI) gene and 11 microsatellite markers were used to test the regional genetic structure and diversity of 318 S. clava individuals from 10 populations within New Zealand. Both markers showed significant differentiation between the northern and southern populations, indicative of minimal pre‐ or post‐border connectivity. Additional statistics further support pre‐ and post‐border differentiation among Port and Harbour populations (i.e. marinas and aquaculture farms). We conclude that New Zealand receives multiple introductions, and that the primary vector for pre‐border incursions and post‐border spread is most likely the extensive influx of recreational vessels that enter northern marinas independent of the Port. This is a timely reminder of the potential for hull‐fouling organisms to expand their range as climates change and open new pathways.     

Pathways of spread of the introduced ascidian Styela clava (Tunicata) in Northern Europe, as revealed by microsatellite markers

Styela clava, a solitary ascidian native to the NW Pacific, has become a conspicuous member of fouling communities in NW European waters. As its natural dispersal appears to be limited, the wide distribution of S. clava along coasts within its introduced range may be attributed to secondary spread assisted by human activities. Here, we used six microsatellite loci to examine the genetic diversity and extent of gene flow among S. clava populations in its European introduced range. Samples were collected from 21 populations within Europe (N = 808), 4 populations within the USA and two populations within the native range (Japan). Large variation in genetic diversity was observed among the European populations but were not explained either by the geographic distance from the first introduction area (i.e. Plymouth, UK) nor by the time elapsed since the introduction. No founder effect was observed in the introduced populations, except possibly in Puget Sound (USA). At least two different introductions occurred in Europe, identified as distinct genetic clusters: northern Danish populations (resembling one Japanese population), and the rest of Europe; a sample from Shoreham (England) possibly represents a third introduction. In North America, the population from the Atlantic was genetically similar to the majority of European populations, suggesting a European origin for populations on this seaboard, while populations from the Pacific coast were genetically similar to the same Japanese population as the Danish populations.     

Inverse correlation of population similarity and introduction date for invasive ascidians

The genomes of many marine invertebrates, including the purple sea urchin and the solitary ascidians Ciona intestinalis and Ciona savignyi, show exceptionally high levels of heterozygosity, implying that these populations are highly polymorphic. Analysis of the C. savignyi genome found little evidence to support an elevated mutation rate, but rather points to a large population size contributing to the polymorphism level. In the present study, the relative genetic polymorphism levels in sampled populations of ten different ascidian species were determined using a similarity index generated by AFLP analysis. The goal was to determine the range of polymorphism within the populations of different species, and to uncover factors that may contribute to the high level of polymorphism. We observe that, surprisingly, the levels of polymorphism within these species show a negative correlation with the reported age of invasive populations, and that closely related species show substantially different levels of genetic polymorphism. These findings show exceptions to the assumptions that invasive species start with a low level of genetic polymorphism that increases over time and that closely related species have similar levels of genetic polymorphism.     

Evolutionary history of the reef fish Anisotremus interruptus (Perciformes: Haemulidae) throughout the Tropical Eastern Pacific

The Tropical Eastern Pacific (TEP) is a dynamic coastal environment characterized by a complex system of oceanic processes and discontinuous rocky habitats. These features, in conjunction with the ecological and physiological characteristics of Anisotremus interruptus, might limit gene flow and shape the evolutionary history of the species. In this study, we investigate the evolutionary history of the reef fish A. interruptus (and its Atlantic sister species A. surinamensis) throughout its range in the TEP, using two mitochondrial (cox1 and cytb) and two nuclear markers (S7 and RAG1). We found three genetic groups of A. interruptus with recent divergence times from the Galapagos Archipelago, Revillagigedo Archipelago, the continental TEP, and A. surinamensis the sister specie from the Atlantic. The haplotype mtDNA networks show A. surinamensis in a central position with respect to Pacific genetic haplogroups, whereas nDNA networks show mixed haplotypes between the four genetic groups. In the species tree, A. surinamensis appears as the sister species of all the Pacific samples and the Galapagos Archipelago population emerges as a genetically distinctive group. The samples from the Revillagigedo Archipelago also constitute a genetic distinctive group, closely related to the continental samples. Continental individuals do not show significant genetic structure and exhibit a population expansion during the Pleistocene. The sandy gaps of the TEP not appear to act as barriers isolating populations of A. interruptus, whereas the open sea gap between the oceanic islands and the continental coast do.     

Dinucleotide microsatellite markers in the genus <Emphasis Type="Italic">Caulerpa</Emphasis>

Caulerpa spp. are clonal green marine algae which often act as invasive species when growing outside their native biogeographical borders. Over the two past decades, Caulerpa taxifolia has spread along the Mediterranean coast, presently occurring at 70 sites and covering nearly 3,000 ha of subtidal area. New genetic markers (microsatellites) have been developed to assess clonal structure and genetic diversity of recently established populations of the invasive species C. taxifolia and Caulerpa racemosa in comparison with populations of the native Caulerpa prolifera in the Mediterranean. Our results show that nine polymorphic markers have been developed for C. prolifera, seven for C. taxifolia, and three for C. racemosa. Genetic diversity in Caulerpa was assessed in two geographical scales: one at a population scale where 40 thalli units were collected from C. prolifera in Cala d’Or, Mallorca, Spain, and another at a species scale, where 30 sample units were analyzed for C. prolifera, 24 for C. taxifolia, and 24 for C. racemosa from different sites in the Mediterranean, Atlantic, and Pacific Ocean. Number of alleles, expected heterozygosity, and marker amplification success are provided in each case.     

Understanding the genetic structure of <Emphasis Type="Italic">Symplocos laurina</Emphasis> Wall. Populations using nuclear gene markers

To characterize the genetic diversity of present populations of Symplocos laurina, which grow in the montane forests in India, we analyzed the DNA sequences of a nuclear gene. Using the 881 bp sequence of cytosolic Glyceraldehyde-3-phosphate dehydrogenase gene, we detected 24 haplotypes among 195 individuals sampled from 14 populations. Two dominant haplotypes were distributed over the entire range of this species in India and several private haplotypes were found. Low genetic diversity within population, high differentiation, number of population specific haplotypes and deviation from neutral evolution characterized the present populations of S. laurina. An analysis of molecular variance indicated the presence of geographic structure within the haplotype distribution. The occurrence of S. laurina preglaciation in India is the most parsimonious explanation for the current geographic structure observed. The populations are presumably ancient and might have spread across its extant distribution range in India through a recent range expansion event.     

Species‐level phylogeographical history of the endemic species Calligonum roborovskii and its close relatives in Calligonum section Medusa (Polygonaceae) in arid north‐western China

Quaternary climatic oscillations appear to have influenced the genetic diversity and evolutionary history of arid‐adapted plants. To understand the processes involved and reveal evolutionary relationships, haplotypes were examined from Calligonum roborovskii, an endemic species occurring in the arid zones across the desert regions of north‐western China, and seven other species also from Calligonum section Medusa, including C. gobicum, C. mongolicum and the narrow endemic species C. ebi‐nuricum, C. pumilum, C. taklimakanense, C. trifarium and C. yengisaricum. Forty‐three haplotypes were identified in 422 individuals from 51 natural populations, from variation of two plastid DNA intergenic spacers (rpl32‐trnL and ycf6‐psbM). A high level of total genetic diversity was found across species for which more than two populations were examined, including C. gobicum, C. mongolicum, C. pumilum and C. roborovskii. A distinct isolation‐by‐distance pattern in each of these species was suggested by the Mantel test, indicating that restricted gene flow caused high genetic differentiation among populations. Three haplotypes were shared by two or three species each, but the other 40 haplotypes were species‐specific. The 43 haplotypes split into three major clades, but not species‐specific lineages; most of the Calligonum species were not reciprocally monophyletic, probably due to incomplete lineage sorting or introgression. The identified haplotypes were dated to 1.97 Mya (95% highest posterior density: 2.95–0.99 Mya) and diverged until the late Pleistocene, possibly linked to aridification and enlargement of deserts caused by climate changes. Variation of desert habitats during the Pleistocene might play a key role in causing the divergence.     

Process of invasiveness among exotic tunicates in Prince Edward Island, Canada

Over the past decade, four exotic tunicates (Styela clava, Ciona intestinalis, Botrylloides violaceus and Botryllus schlosseri) have been reported in the Brudenell estuary in Prince Edward Island (PEI), Canada. Styela clava was the first exotic tunicate to arrive in 1997, rapidly establishing, spreading, invading, and eventually becoming a nuisance in several estuaries of PEI. In the Brudenell estuary, S. clava remained the only exotic nuisance tunicate until 2003. In the fall of 2004, the vase tunicate C. intestinalis, was reported in low abundance, followed by the two colonial species, B. schlosseri and B. violaceus, reported in the spring of 2005. The abundance of C. intestinalis rapidly increased post-introduction, eventually replacing S. clava as the foremost nuisance species on mussel farms in the estuary. To date, C. intestinalis continues to colonize this estuary at epidemic proportions, resulting in the continuing drop of S. clava abundance. The current abundance of C. intestinalis is estimated at 5 cm⁻², which is similar to S. clava abundance at its height in 2003. The 2006 abundance of S. clava is estimated to have fallen to near 0 cm⁻². The dominance of C. intestinalis as a fouling organism on mussel farms is considered a serious threat to this aquaculture industry, mainly due to its unmanageable weight. The process of the detection, establishment, invasiveness, and eventual rise to nuisance level of exotic tunicates in the Brudenell River is presented.     

Genetic evidence for recent range fragmentation and severely restricted dispersal in the critically endangered Sierra Madre Sparrow, Xenospiza baileyi

Assessing patterns of genetic structure and diversity of threatened species has become an essential tool for determining conservation status and designing management strategies. We examine the genetic structure of the Sierra Madre sparrow (Xenospiza baileyi), a species restricted to fragmented patches of subalpine bunchgrass in three small isolated areas of northwestern and central Mexico. Coding and non-coding regions of mtDNA (1,878 bp) from individuals of the only three known populations revealed the existence of a single major lineage, with closely related haplotypes being shared between populations across the range. The sharing of haplotypes between the distant northwest and central populations (~800 km) suggests a recent fragmentation of a formerly contiguous population. Despite a lack of large-scale phylogeographic structure, haplotype frequencies at local scales revealed significant genetic differentiation and high F _ST values between all three remaining populations, even between localities separated by less than 12 km. These results suggest restricted gene flow and limited dispersal, likely due to the species’ inability to cross areas of unsuitable habitat. On the basis of genetic interchange and ecological equivalence criteria, we recommend that the species be managed as a single unit, permitting the strengthening of the small population in the northwest with individuals from central Mexico, and/or the translocation of individuals to new areas of suitable habitat.     

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype‐by‐sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.     

Isolation and characterization of twelve polymorphic microsatellite markers for the invasive ascidian Styela clava (Tunicata)

The invasive solitary ascidian Styela clava has spread extensively along European coasts since its first occurrence in the early 1950s. In order to characterize spatial and temporal patterns of genetic change during its establishment and subsequent spread, we developed 12 species‐specific loci from an enriched microsatellite library. Polymorphism was explored in one native and two introduced populations (N = 31 ± 1). Number of alleles per locus varied from two to 13 (mean = 7). The average expected heterozygosity within populations ranged from 0.539 to 0.580.