Isolation with differentiation followed by expansion with admixture in the tunicate Pyura chilensis

Chemical compositions and nutritional profiles of two edible tunicate species (Halocynthia roretzi and Halocynthia aurantium)

As an abundant marine bioresource, tunicates could be exploited in the food industry. However, limited knowledge of their chemical composition and nutritional profiles prohibited further application. In this study, two common edible tunicate species, Halocynthia roretzi (HR) and Halocynthia aurantium (HA), were subjected to comprehensive composition analysis in terms of moisture, protein, lipids, cellulose, ash, amino acids, fatty acids, non-cellulose carbohydrates and minerals. Reddish HR was much bigger than purple HA with respect to body length and weight, and their moisture fell within 82.98 %-90.92 %. The non-edible outer shell part (OS) and edible internal organs part (IO) had a dry weight ratio of around 3:2 for both two species. Generally, for both HR and HA, IO was more abundant in protein and lipids. In contrast, OS had much higher cellulose contents, confirming the better suitability of IO as a nutritional seafood. IO was richer in essential amino acids and unsaturated fatty acids, while OS had more abundant saturated fatty acids. The detected non-cellulose monosugars ranged from 0.47 % to 1.18 % and indicated the presence of some sulfated glycans. IO of HR had higher contents of essential minerals, such as Cu, Zn, and Fe, while IO of HA showed a higher K content. To sum up, this study identified the chemical composition and nutritional profile variations among different tunicate species and various dissected parts, guiding the development of specific strategies to exploit tunicates for proper food applications.

Shedding light on variation in reproductive success through studies of population genetic structure in a Southeast Pacific Coast mussel

Phylogeography often focuses on the spatial dimension of genetic diversity, rarely including the temporal dynamics occurring interannually among local populations, which can provide insight into past variations in reproductive success. Currently, there is an intense aquaculture industry of Mytilus spp. on the Southeast Pacific Coast which depends entirely on the spat released by natural populations forming a relevant and sensitive social-ecological system. Temporal and spatial spat variability from natural mussel beds could be related to interannual reproductive dynamics with variable reproductive success and recruitment, which leave genetic signatures. Temporal and spatial genetic structure was evaluated in six natural beds in the Southeast Pacific (from 39°25'S to 43°07'S) on the most abundant and widespread Mytilus lineage detected, Mytilus cf. chilensis, in 4 consecutive years. Analyses included data from >180 individuals per year, with a total of 751 (mitochondrial COI) and 747 (nuclear H1) individuals, respectively. Overall, both markers showed high haplotype diversity and low spatial and temporal genetic differentiation. Likely, the high dispersal capacity of Mytilus cf. chilensis maintains population homogeneity and prevents diversity erosion. The slight differences in genetic variance of COI were better explained by differences among sites (space), and conversely, the H1 genetic variance was better explained by interannual (temporal) comparisons, which could explain temporal variability in spat availability. This study highlights the important insights achieved with the evaluation of both temporal and spatial population genetic structures in marine species with high reproductive output, which can condition the success and sustainability of the relevant social-ecological system.

The right tool for the right question: contrasting biogeographic patterns in the notothenioid fish Harpagifer spp. along the Magellan Province

Molecular-based analysis has become a fundamental tool to understand the role of Quaternary glacial episodes. In the Magellan Province in southern South America, ice covering during the last glacial maximum (20 ka) radically altered the landscape/seascape, speciation rates and distribution of species. For the notothenioid fishes of the genus Harpagifer, in the area are described two nominal species. Nevertheless, this genus recently colonized South America from Antarctica, providing a short time for speciation processes. Combining DNA sequences and genotyping-by-sequencing SNPs, we evaluated the role of Quaternary glaciations over the patterns of genetic structure in Harpagifer across its distribution in the Magellan Province. DNA sequences showed low phylogeographic structure, with shared and dominant haplotypes between nominal species, suggesting a single evolutionary unit. SNPs identified contrastingly two groups in Patagonia and a third well-differentiated group in the Falkland/Malvinas Islands with limited and asymmetric gene flow. Linking the information of different markers allowed us to infer the relevance of postglacial colonization mediated by the general oceanographic circulation patterns. Contrasting rough- and fine-scale genetic patterns highlights the relevance of combined methodologies for species delimitation, which, depending on the question to be addressed, allows discrimination among phylogeographic structure, discarding incipient speciation, and contemporary spatial differentiation processes.

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.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Secondary contacts and genetic admixture shape colonization by an amphiatlantic epibenthic invertebrate

Research on the genetics of invasive species often focuses on patterns of genetic diversity and population structure within the introduced range. However, a growing body of literature is demonstrating the need to study how native genotypes affect both ecological and evolutionary mechanisms within the introduced range. Here, we used genotyping-by-sequencing to study both native and introduced ranges of the amphiatlantic marine invertebrate Ciona intestinalis. A previous study using microsatellites analysed samples collected along the Swedish west coast and showed the presence of genetically distinct lineages in deep and shallow waters. Using 1,653 single nucleotide polymorphisms (SNPs) from newly collected samples (285 individuals), we first confirmed the presence of this depth-defined genomic divergence along the Swedish coast. We then used approximate Bayesian computation to infer the historical relationship among sites from the North Sea, the English Channel and the northwest Atlantic and found evidence of ancestral divergence between individuals from deep waters off Sweden and individuals from the English Channel. This divergence was followed by a secondary contact that led to a genetic admixture between the ancestral populations (i.e., deep Sweden and English Channel), which originated the genotypes found in shallow Sweden. We then revealed that the colonization of C. intestinalis in the northwest Atlantic was as a result of an admixture between shallow Sweden and the English Channel genotypes across the introduced range. Our results showed the presence of both past and recent genetic admixture events that together may have promoted the successful colonizations of C. intestinalis. Our study suggests that secondary contacts potentially reshape the evolutionary trajectories of invasive species through the promotion of intraspecific hybridization and by altering both colonization patterns and their ecological effects in the introduced range.

Fine-scale hierarchical genetic structure and kinship analysis of the ascidian Pyura chilensis in the southeastern Pacific

Studying population structure and genetic diversity at fine spatial scales is key for a better understanding of demographic processes that influence population connectivity. This is particularly important in marine benthic organisms that rely on larval dispersal to maintain connectivity among populations. Here, we report the results of a genetic survey of the ascidian Pyura chilensis from three localities along the southeastern Pacific. This study follows up on a previous report that described a genetic break in this region among localities only 20 km apart. By implementing a hierarchical sampling design at four spatial levels and using ten polymorphic microsatellite markers, we test whether differences in fine-scale population structure explain the previously reported genetic break. We compared genetic spatial autocorrelations, as well as kinship and relatedness distributions within and among localities adjacent to the genetic break. We found no evidence of significant autocorrelation at the scale up to 50 m despite the low dispersal potential of P. chilensis that has been reported in the literature. We also found that the proportion of related individuals in close proximity (<1 km) was higher than the proportion of related individuals further apart. These results were consistent in the three localities. Our results suggest that the spatial distribution of related individuals can be nonrandom at small spatial scales and suggests that dispersal might be occasionally limited in this species or that larval cohorts can disperse in the plankton as clustered groups. Overall, this study sheds light on new aspects of the life of this ascidian as well as confirms the presence of a genetic break at 39°S latitude. Also, our data indicate there is not enough evidence to confirm that this genetic break can be explained by differences in fine-scale genetic patterns among localities.

Genetic and morphological divergence at a biogeographic break in the beach-dwelling brooder Excirolana hirsuticauda Menzies (Crustacea, Peracarida)

BACKGROUND

There is a biogeographic break located at 30°S in the southeast Pacific, in a coastal area of strong environmental discontinuities. Several marine benthic taxa with restricted dispersal have a coincident phylogeographic break at 30°S, indicating that genetic structure is moulded by life history traits that limit gene flow and thereby promote divergence and speciation. In order to evaluate intraspecific divergence at this biogeographic break, we investigated the genetic and morphological variation of the directly developing beach isopod Excirolana hirsuticauda along 1900 km of the southeast Pacific coast, across 30°S.

RESULTS

The COI sequences and microsatellite data both identified a strong discontinuity between populations of E. hirsuticauda to the north and south of 30°S, and a second weaker phylogeographic break at approximately 35°S. The three genetic groups were evidenced by different past demographic and genetic diversity signatures, and were also clearly distinguished with microsatellite data clustering. The COI sequences established that the genetic divergence of E. hirsuticauda at 30°S started earlier than divergence at 35°. Additionally, the three groups have different past demographic signatures, with probable demographic expansion occurring earlier in the southern group (south of 35°S), associated with Pleistocene interglacial periods. Interestingly, body length, multivariate morphometric analyses, and the morphology of a fertilization-related morphological character in males, the appendix masculina, reinforced the three genetic groups detected with genetic data.

CONCLUSIONS

The degree of divergence of COI sequences, microsatellite data, and morphology was concordant and showed two geographic areas in which divergence was promoted at differing historical periods. Variation in the appendix masculina of males has probably promoted reproductive isolation. This variation together with gene flow restrictions promoted by life history traits, small body size, oceanographic discontinuities and sandy-beach habitat continuity, likely influenced species divergence at 30°S in the southeast Pacific coast. The degree of genetic and morphological differentiation of populations to the north and south of 30°S suggests that E. hirsuticauda harbours intraspecific divergence consistent with reproductive isolation and an advanced stage of speciation. The speciation process within E. hirsuticauda has been shaped by both restrictions to gene flow and a prezygotic reproductive barrier.

Lineage divergence, local adaptation across a biogeographic break, and artificial transport, shape the genetic structure in the ascidian Pyura chilensis

Marine benthic organisms inhabit a heterogeneous environment in which connectivity between populations occurs mainly through dispersive larval stages, while local selective pressures acting on early life history stages lead to non-random mortality, shaping adaptive genetic structure. In order to test the influence of local adaptation and neutral processes in a marine benthic species with low dispersal, in this study we used Genotyping by Sequencing technology to compare the neutral and putatively selected signals (neutral and outlier loci, respectively) in SNPs scattered throughout the genome in six local populations of the commercially exploited ascidian Pyura chilensis along the southeast Pacific coast (24°–42°S). This species is sessile as an adult, has a short-lived larval stage, and may also be dispersed by artificial transport as biofouling. We found that the main signal in neutral loci was a highly divergent lineage present at 39°S, and a subjacent signal that indicated a separation at 30°S (north/south), widely reported in the area. North/south separation was the main signal in outlier loci, and the linage divergence at 39°S was subjacent. We conclude that the geographic structure of the genetic diversity of outlier and neutral loci was established by different strengths of environmental, historical and anthropogenic factors.

Large-scale gene flow in the barnacle Jehlius cirratus and contrasts with other broadly-distributed taxa along the Chilean coast

We evaluate the population genetic structure of the intertidal barnacle Jehlius cirratus across a broad portion of its geographic distribution using data from the mitochondrial cytochrome oxidase I (COI) gene region. Despite sampling diversity from over 3,000 km of the linear range of this species, there is only slight regional structure indicated, with overall Φ CT of 0.036 (p < 0.001) yet no support for isolation by distance. While these results suggest greater structure than previous studies of J. cirratus had indicated, the pattern of diversity is still far more subtle than in other similarly-distributed species with similar larval and life history traits. We compare these data and results with recent findings in four other intertidal species that have planktotrophic larvae. There are no clear patterns among these taxa that can be associated with intertidal depth or other known life history traits.

Phylogeographic structure and deep lineage diversification of the red alga Chondrus ocellatus Holmes in the Northwest Pacific

A major goal of phylogeographic analysis using molecular markers is to understand the ecological and historical variables that influence genetic diversity within a species. Here, we used sequences of the mitochondrial Cox1 gene and nuclear internal transcribed spacer to reconstruct its phylogeography and demographic history of the intertidal red seaweed Chondrus ocellatus over most of its geographical range in the Northwest Pacific. We found three deeply separated lineages A, B and C, which diverged from one another in the early Pliocene-late Miocene (c. 4.5-7.7 Ma). The remarkably deep divergences, both within and between lineages, appear to have resulted from ancient isolations, accelerated by random drift and limited genetic exchange between regions. The disjunct distributions of lineages A and C along the coasts of Japan may reflect divergence during isolation in scattered refugia. The distribution of lineage B, from the South China Sea to the Korean Peninsula, appears to reflect postglacial recolonizations of coastal habitats. These three lineages do not coincide with the three documented morphological formae in C. ocellatus, suggesting that additional cryptic species may exist in this taxon. Our study illustrates the interaction of environmental variability and demographic processes in producing lineage diversification in an intertidal seaweed and highlights the importance of phylogeographic approaches for discovering cryptic marine biodiversity.

Transcriptome characterization of the ascidian Pyura chilensis using 454-pyrosequencing data from two distant localities on the southeast Pacific

This study describes the results from transcriptomes sequenced by 454-pyrosequencing from two populations separated by 10° of latitude of the endemic tunicate Pyura chilensis. Most transcripts were assembled in 43,972 contigs with an average length of 842 nucleotides. De novo assembly revealed that less than 30% of the contigs were annotated to Gene Ontology terms. A total of 71,662 single nucleotide polymorphisms (SNPs) were detected in 14,712 contigs. In silico differential expression of contigs annotated for SNPs revealed several genes differentially expressed in individuals collected from both populations. The present genomic resource will provide the basis to perform functional genomics on the species via the localization of genomic markers that can aid in determining levels of local adaptation, overall genetic structure and the genetic assessment of restocking programs for this species.

Phylogeographic structure in benthic marine invertebrates of the southeast Pacific coast of Chile with differing dispersal potential

The role of dispersal potential on phylogeographic structure, evidenced by the degree of genetic structure and the presence of coincident genetic and biogeographic breaks, was evaluated in a macrogeographic comparative approach along the north-central coast of Chile, across the biogeographic transition zone at 30°S. Using 2,217 partial sequences of the mitochondrial Cytochrome Oxidase I gene of eight benthic invertebrate species along ca. 2,600 km of coast, we contrasted dispersal potential with genetic structure and determined the concordance between genetic divergence between biogeographic regions and the biogeographic transition zone at 30°S. Genetic diversity and differentiation highly differed between species with high and low dispersal potential. Dispersal potential, sometimes together with biogeographic region, was the factor that best explained the genetic structure of the eight species. The three low dispersal species, and one species assigned to the high dispersal category, had a phylogeographic discontinuity coincident with the biogeographic transition zone at 30°S. Furthermore, coalescent analyses based on the isolation-with-migration model validate that the split between biogeographic regions north and south of 30°S has a historic origin. The signatures of the historic break in high dispersers is parsimoniously explained by the homogenizing effects of gene flow that have erased the genetic signatures, if ever existed, in high dispersers. Of the four species with structure across the break, only two had significant albeit very low levels of asymmetric migration across the transition zone. Historic processes have led to the current biogeographic and phylogeographic structure of marine species with limited dispersal along the north-central coast of Chile, with a strong lasting impact in their genetic structure.