Long-term population genetic structure of an invasive urochordate: the ascidian Botryllus schlosseri

Geography and developmental plasticity shape post-larval thermal tolerance in the golden star tunicate, Botryllus schlosseri

Local adaptation and phenotypic plasticity play key roles in mediating organisms' ability to respond to spatiotemporal variation in temperature. These two processes often act together to generate latitudinal or elevational clines in acute temperature tolerance. Phenotypic plasticity is also subject to local adaptation, with the expectation that populations inhabiting more variable environments should exhibit greater phenotypic plasticity of thermal tolerance. Here we examine the potential for local adaptation and developmental plasticity of thermal tolerance in the widespread invasive tunicate Botryllus schlosseri. By comparing five populations across a thermal gradient spanning 4.4° of latitude in the northwest Atlantic, we demonstrate that warmer populations south of the Gulf of Maine exhibit significantly increased (∼0.2 °C) post-larval temperature tolerance relative to the colder populations within it. We also show that B. schlosseri post-larvae possess a high degree of developmental plasticity for this trait, shifting their median temperature of survival (LT50) upwards by as much as 0.18 °C per 1 °C increase in environmental temperature. Lastly, we found that populations vary in their degrees of developmental plasticity, with populations that experience more pronounced short-term temperature variability exhibiting greater developmental plasticity, suggesting the local adaptation of developmental plasticity. By comparing the thermal tolerance of populations across space and through time, we demonstrate how geography and developmental plasticity have shaped thermal tolerance in B. schlosseri. These results help inform our understanding of how species are able to adjust their thermal physiology in new environments, including those encountered during invasion and under increasingly novel climate conditions.

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.

Chimerism and population dieback alter genetic inference related to invasion pathways and connectivity of biofouling populations on artificial substrata

Disentangling pathways by which nonindigenous species expand and spread regionally remains challenging. Molecular ecology tools are often employed to determine the origins and spread of introduced species, but the complexities of some organisms may be reducing the efficacy of these tools. Some colonial species exhibit complexities by way of chimerism and winter colony regression, which may alter the genetic diversity of populations and mask the connectivity occurring among them. This study uses nuclear microsatellite data and simple GIS-based modeling to investigate the influence of chimerism and winter regression on the genetic diversity and patterns of genetic population connectivity among colonies of Didemnum vexillum on artificial substrates. Colonies sampled in summer were shown to form a metapopulation, with high levels of admixture, extreme outcrossing, and some substructure. These patterns were consistent within the subsampled winter colonies and with the inclusion of chimeric data. However, allelic richness and diversity were significantly different between winter and summer samples, altering interpretations relating to population connectivity and pelagic larval duration. This study demonstrates the importance of including seasonal sampling and imperative life history traits in genetic studies for clear interpretations and the successful management of introduced species.

Genetic homogeneity of the invasive lionfish across the Northwestern Atlantic and the Gulf of Mexico based on Single Nucleotide Polymorphisms

Despite the devastating impact of the lionfish (Pterois volitans) invasion on NW Atlantic ecosystems, little genetic information about the invasion process is available. We applied Genotyping by Sequencing techniques to identify 1,220 single nucleotide polymorphic sites (SNPs) from 162 lionfish samples collected between 2013 and 2015 from two areas chronologically identified as the first and last invaded areas in US waters: the east coast of Florida and the Gulf of Mexico. We used population genomic analyses, including phylogenetic reconstruction, Bayesian clustering, genetic distances, Discriminant Analyses of Principal Components, and coalescence simulations for detection of outlier SNPs, to understand genetic trends relevant to the lionfish’s long-term persistence. We found no significant differences in genetic structure or diversity between the two areas (F_STp-values > 0.01, and t-test p-values > 0.05). In fact, our genomic analyses showed genetic homogeneity, with enough gene flow between the east coast of Florida and Gulf of Mexico to erase previous signals of genetic divergence detected between these areas, secondary spreading, and bottlenecks in the Gulf of Mexico. These findings suggest rapid genetic changes over space and time during the invasion, resulting in one panmictic population with no signs of divergence between areas due to local adaptation.

Invasiveness, chimerism and genetic diversity

Adaptation for invasiveness should comprise the capability to exploit and prosper in a wide range of ecological conditions and is therefore expected to be associated with a certain level of genetic diversity. Paradoxically, however, invasive populations are established by only a few founders, resulting in low genetic diversity. As a conceivable way of attaining high genetic diversity and high variance of gene expression even when a small number of founders is involved in invasiveness, I suggest here chimerism, a fusion between different individuals-a common phenomenon found in numerous phyla. The composite entity offers the chimeric organism genetic flexibility and higher inclusive fitness that depends on the joint genomic fitness of the original partners. The ability to form a chimeric entity is also applied to subsequent generations, and consequently, the level of genetic diversity does not decline over generations of population establishment following invasion.

Phylogenetics, biogeography and population genetics of the ascidian Botryllus schlosseri in the Mediterranean Sea and beyond

The wide distribution of the ascidian Botryllus schlosseri along the Mediterranean coasts has been documented since the eighteenth century. However, despite copious documentation, analyses of dispersal modes and genetic profiles were limited to local populations or restricted regions. In order to get a pan-Mediterranean overview, 288 specimens from 11 populations of B. schlosseri from the western and eastern Mediterranean basins were sampled and analyzed using five microsatellite loci and COI sequences. Both molecular markers revealed high polymorphisms, with 182 microsatellites alleles and 54 COI haplotypes. Overall, F_st, D_est, and COI Ф_pt values were 0.146, 0.635 and 0.322, respectively, reflecting a high genetic diversity and a significant genetic structure as compared to other B. schlosseri populations worldwide, reflected by substantially higher values for effective number of alleles (N_e) in the Mediterranean. A phylogenetic analysis of the COI sequences resulted in four distinct clades and two molecular operational taxonomic units (OTUs). We recorded a stronger genetic structure among the populations of the eastern basin compared to the western basin (microsatellites F_st=0.217 versus 0.082; COI Ф_pt=0.416 versus 0.171), suggesting either a restricted connectivity between the basins or a stronger genetic drift in each basin. The occurrence of two OTUs and different ecological conditions may also contribute to this finding. Mean Nei's genetic distance in the eastern Mediterranean populations was more than three times higher compared to the western basin. No correlation was observed between geographic and genetic distances (Mantel test), suggesting that maritime transport is the main dispersal vector of B. schlosseri colonies. The possibility that the Mediterranean is a center of diversity for B. schlosseri, and probably its site of origin, is further discussed.

Stochasticity in space, persistence in time: genetic heterogeneity in harbour populations of the introduced ascidian Styela plicata

Spatio-temporal changes in genetic structure among populations provide crucial information on the dynamics of secondary spread for introduced marine species. However, temporal components have rarely been taken into consideration when studying the population genetics of non-indigenous species. This study analysed the genetic structure of Styela plicata, a solitary ascidian introduced in harbours and marinas of tropical and temperate waters, across spatial and temporal scales. A fragment of the mitochondrial gene Cytochrome Oxidase subunit I (COI) was sequenced from 395 individuals collected at 9 harbours along the NW Mediterranean coast and adjacent Atlantic waters (> 1,200 km range) at two time points 5 years apart (2009 and 2014). The levels of gene diversity were relatively low for all 9 locations in both years. Analyses of genetic differentiation and distribution of molecular variance revealed strong genetic structure, with significant differences among many populations, but no significant differences among years. A weak and marginally significant correlation between geographic distance and gene differentiation was found. Our results revealed spatial structure and temporal genetic homogeneity in S. plicata, suggesting a limited role of recurrent, vessel-mediated transport of organisms among small to medium-size harbours. Our study area is representative of many highly urbanized coasts with dense harbours. In these environments, the episodic chance arrival of colonisers appears to determine the genetic structure of harbour populations and the genetic composition of these early colonising individuals persists in the respective harbours, at least over moderate time frames (five years) that encompass ca. 20 generations of S. plicata.

Population genetics features for persistent, but transient, Botryllus schlosseri (Urochordata) congregations in a central Californian marina

The colonial tunicate Botryllus schlosseri is a globally distributed, invasive ascidian that has colonized the Californian coasts of the USA during the mid-late 1940s and has, since the late 1980s, spread north to Washington. This study analyzes the population genetic characteristics of transient populations residing at the Elkhorn Yacht-Club (EYC), in central California (seven sessions, 1996-2008), which suffered periodic catastrophes caused by episodic fresh-water floods and a single sampling session (in the year 2001) of five West-Coast populations using the mtDNA COI gene and five microsatellite markers. EYC microsatellite results were further compared with the closely situated but persistent population of the Santa Cruz Harbor (SCH) to understand the impact on EYC population regeneration processes after the 2005-flood catastrophe. All microsatellites were highly polymorphic, revealing a large number of unique alleles at different sampling dates. Whereas pairwise θ did not reveal significant differences between the EYC time-series samplings, the overall θ was significant, as it was between all the 2001 West Coast populations. The most likely cluster number was 3 for the EYC samples whereas two K values were obtained (2 and 5) for the 2001 samples. Tajima's D and Fu's/Fs tests did not reject the null hypothesis for COI neutral evolution, except for in the EYC-2000, 2007 and two 2001 samplings. The wide geographical range of the analyses has indicated that following the EYC 2005-flood catastrophe, newcomers could have originated from neighboring populations, from deep-water colonies that may have escaped the 2005 low salinity event, or less expectedly, from far away West-Coast populations, while revealing that the SCH population is the most probable source for the EYC population.

Population structure of the greenhouse whitefly, Trialeurodes vaporariorum (Westwood), an invasive species from the Americas, 60 years after invading China

Though the greenhouse whitefly, Trialeurodes vaporariorum (Westwood) (Hemiptera: Aleyrodidae) was introduced into China more than 60 years ago, the genetic diversity and structure of this exotic insect pest and virus vector have not been studied. To investigate the population genetic characteristics of this invasive species and to identify potential invasion routes, the genetic diversity and population structure of 17 collections of T. vaporariorum from nine provinces in China were analyzed using seven microsatellite loci. The results of the analyses indicated that the genetic diversity for the populations examined from the four provinces: Jilin, Ningxia, Guizhou and Qinghai, was lower than the genetic diversity of populations from the five provinces: Yunnan, Shandong, Shanxi, Liaoning, and Gansu. The T. vaporariorum populations analyzed in this study grouped as two distinct genetic clusters based on the analysis using STRUCTURE, whereas, 8 clusters were identified based on the BAPS analysis. Of the 136 genetic distance (Fst) values, 128 (94%) were associated with a significant exact test. However, there was no significant relationship between Fst and geographical distance. These results demonstrate that populations of T. vaporariorum in China exhibit significant genetic differentiation, indicating the likelihood that multiple introductions of T. vaporariorum into China have occurred. Also, the populations collected from the provinces of Jilin, Ningxia, Guizhou and Qinghai appear to represent secondary introductions originating from other Chinese provinces.

Can novel genetic analyses help to identify low-dispersal marine invasive species?

Genetic methods can be a powerful tool to resolve the native versus introduced status of populations whose taxonomy and biogeography are poorly understood. The genetic study of introduced species is presently dominated by analyses that identify signatures of recent colonization by means of summary statistics. Unfortunately, such approaches cannot be used in low-dispersal species, in which recently established populations originating from elsewhere in the species' native range also experience periods of low population size because they are founded by few individuals. We tested whether coalescent-based molecular analyses that provide detailed information about demographic history supported the hypothesis that a sea squirt whose distribution is centered on Tasmania was recently introduced to mainland Australia and New Zealand through human activities. Methods comparing trends in population size (Bayesian Skyline Plots and Approximate Bayesian Computation) were no more informative than summary statistics, likely because of recent intra-Tasmanian dispersal. However, IMa2 estimates of divergence between putatively native and introduced populations provided information at a temporal scale suitable to differentiate between recent (potentially anthropogenic) introductions and ancient divergence, and indicated that all three non-Tasmanian populations were founded during the period of European settlement. While this approach can be affected by inaccurate molecular dating, it has considerable (albeit largely unexplored) potential to corroborate nongenetic information in species with limited dispersal capabilities.