Species replacement along a linear coastal habitat: phylogeography and speciation in the red alga Mazzaella laminarioides along the south east Pacific

Impact of persistent barrier to gene flow and catastrophic events on red algae evolutionary history along the Chilean coast

Historical vicariance events, linked to the existence of stable physical barriers to gene flow, generate concordant genetic breaks in co-distributed species while stochastic processes (e.g., costal uplift) could cause species-specific genetic breaks as a result of local strong demographic bottlenecks or extinction. In Chile, previous studies show that the area of the 30°S-33°S could correspond to a stable barrier to gene flow that have affected the genetic structure of various algae and marine invertebrates. Here we sequenced two organellar genes (COI and rbcL) in four taxonomically accepted co-distributed red seaweeds species characterized by a low dispersal potential: Mazzaella laminarioides, M. membranacea, Asterfilopsis disciplinalis, and Ahnfeltiopsis vermicularis. Our results revealed the existence of ten strongly differentiated linages in the taxa studied. Strong genetic breaks, concordant in both space and time (divergence estimated to have occurred some 2.9-12.4 million years ago), were observed between taxa distributed across the 33°S. Conversely, in the Central/South part of the Chilean coast, the localization of the genetic breaks/sub-structure observed varied widely (36°S, 38°S, 39°S, and 40°S). These results suggest that a major historical vicariance event has modeled the genetic structure of several Chilean marine organisms in the north of the Chilean coast during the mid-Miocene, while more recent stochastic events and genetic drift could be the driving forces of genetic divergence/structuration in the central-southern part of the coast.

Linking Acrosome Size and Genetic Divergence in an Inter-Oceanic Mussel from the Pacific and Atlantic Coasts: A Case of Incipient Speciation?

In recent years, advances in analyses of the sperm morphology and genetics of Perumytilus purpuratus have allowed to two evolutionary scenarios for this mussel to be suggested: (1) the scenario of cryptic species and (2) the scenario of incipient or in progress speciation. For a better understanding of the evolutionary history of P. purpuratus, we performed extensive sampling along a latitudinal gradient of ca. 7180 km of coastline-from the Southern Pacific Ocean to the Atlantic Ocean-and we delved deeper into the sperm morphology of P. purpuratus, exploring its association with the phylogeny and population genetics to determine whether the variability in sperm traits between the northern and southern regions was a signal of cryptic or incipient species. Overall, our results showed that sperm sizes were strongly correlated with the genetic structure in males of P. purpuratus. We identified at 37° S on the Pacific coast a coincident break of both sperm size and genetic disruption that can be explained by historical events and postglacial recolonization as causal phenomena for the observed divergences. Furthermore, evidence of genetic admixture between lineages was found at 38° S, suggesting the presence of an introgressive hybridization zone and incomplete reproductive isolation in an in fraganti or incipient speciation process.

Biotopes of the intertidal zone in Clarence Island (south of the Strait of Magellan)

Background

The characteristics of the Strait of Magellan promote the formation of unique environments, with diverse habitats and marine organisms. This fragmentation of the landscape generates diverse little-explored ecological associations, especially in the zone of sub-Antarctic islands of the Tierra del Fuego archipelago. One way to address this lack of knowledge is through the biotope characterization methodology, with ecological units composed of the habitat and the communities associated with these environments, obtaining data and information on the dominant and incidental taxonomic groups. This is a good research model to conduct baseline studies in coastal benthic marine environments.

New information

A data set in Darwin Core standard is presented of the species that make up the intertidal biotopes of Clarence Island (Tierra del Fuego Archipelago, south of the Strait of Magellan). This includes 50 identified species and the specific coordinates for each sampled location, with a total of 1400 georeferenced records. Mollusks were the most diverse taxon with 21 species, followed by algae (14 species). Sessile organisms such as the barnacles Elminiuskingii and Austromegabalanuspsittacus predominate in these ecosystems, followed by bivalve mollusks such as Choromytiluschorus and Mytiluschilensis, which together with Nacellamagellanica and the alga Hildenbrandia sp. make up more than 50% of the total records. The inclusion of biotope patterns in this study complements the information on benthic marine flora and fauna in the intertidal zone, including new records for the coast in the Clarence Island area, which is within the boundary of the Kawésqar National Park.

Genetic Diversity and Phylogeography of a Turf-Forming Cosmopolitan Marine Alga, Gelidium crinale (Gelidiales, Rhodo-Phyta)

Cosmopolitan species are rare in red algae, which have a low-dispersal capacity unless they are dispersed by human-mediated introductions. Gelidium crinale, a turf-forming red alga, has a widespread distribution in tropical and temperate waters. To decipher the genetic diversity and phylogeography of G. crinale, we analyzed mitochondrial COI-5P and plastid rbcL sequences from collections in the Atlantic, Indian, and Pacific Oceans. Phylogenies of both markers statistically supported the monophyly of G. crinale, with a close relationship to G. americanum and G. calidum from the Western Atlantic. Based on the molecular analysis from these materials, Pterocladia heteroplatos from India is here merged with G. crinale. Phylogeny and TCS networks of COI-5P haplotypes revealed a geographic structure of five groups: (i) Atlantic-Mediterranean, (ii) Ionian, (iii) Asian, (iv) Adriatic-Ionian, and (v) Australasia-India-Tanzania-Easter Island. The most common ancestor of G. crinale likely diverged during the Pleistocene. The Bayesian Skyline Plots suggested the pre-LGM population expansion. Based on geographical structure, lineage-specific private haplotypes, the absence of shared haplotypes between lineages, and AMOVA, we propose that the cosmopolitan distribution of G. crinale has been shaped by Pleistocene relicts. The survival of the turf species under environmental stresses is briefly discussed.

Assessing diversity of King Crab Lithodes spp. in the south‐eastern pacific using phylogeny and molecular species delimitation methods

The purpose of this study was to test the hypothesis that the genetic diversity of commercially significant species of King Crabs (Lithodes spp.) along the south‐eastern Pacific (SEP) comprises different independent evolutionary units (IEUs) with spatially isolated distribution. Nine localities from inner and open waters along the SEP Chilean coast (39°S‐55°S) were sampled. We analyzed sequences from 173 individuals for the mitochondrial gene Cytochrome oxidase I (COX‐I), 151 individuals for the Internal Transcribed Spacer 1 (ITS) and 135 for the structural ribosomal RNA (28S). Genetic delimitation was performed through three analytical methods: ABGD, GMYC, and its Bayesian implementation, bGMYC. Bayesian phylogenetic analyses and haplotype networks were also performed. Divergence time between clades was assessed for the COX‐I marker and estimated from known evolutionary rates for this marker in other crustacean species and fossil calibration from other Anomuran species. Delimitation analyses, phylogenetic analyses, and mitochondrial haplotype networks suggested the presence of two deeply divergent mitochondrial lineages of Lithodes in the SEP, referred to as Clade1 and Clade 2. Nuclear markers showed low phylogenetic resolution and therefore were unsuitable for molecular species delimitation. Divergence time analysis of the mitochondrial lineages suggests a separation between Clades of approximately 2.3 Mya. The divergence time obtained suggested that Pliocene glaciations and deglaciations cycles could be involved in hybridization events between Lithodes IEUs at southern tip of South American coasts. The different frequencies of Lithodes haplotypes in inner and open water environments along SEP coasts could be explained by events such as the last glacial maximum or by differences in the adaptation of each clade to different environments. These findings support the necessity of evaluating the taxonomic status of Lithodes individuals found along SEP coasts under an integrative taxonomy approach or through markers with other evolution rates than those already used.

Independent Effects of Species Removal and Asynchrony on Invariability of an Intertidal Rocky Shore Community

Ecological stability depends on interactions between different levels of biological organization. The insurance effects occur when increasing species diversity leads to more temporally invariable (i.e., more stable) community-level properties, due in part to asynchronous population-level fluctuations. While the study of insurance effects has received considerable attention, the role of dominant species that contribute with particular functional traits across different level of organizations is less understood. Using a field-based manipulative experiment, we investigated how species richness and different types of parameters at the population level, such as the invariability of dominants, population invariability, and population asynchrony, influence the community invariability. The experiment involved the repetitive removal of the canopy forming alga Mazzaella laminarioides (hereafter “Mazzaella”) during 32 months in two rocky intertidal sites of northern-central Chile. We predicted that the invariability of dominants enhances community invariability, that the effect of multispecies population-level parameters on community invariability are dependent on species richness, and that subdominant algae are unable to fully compensate the loss of canopies of the dominant species. Biomass of algae and mobile invertebrates was quantified over time. We observed independent effects of Mazzaella removal and community-wide asynchrony on community invariability. While canopy removal reduced community invariability, population asynchrony boosted community invariability regardless of the presence of canopies. In addition, filamentous and foliose algae were unable to compensate the loss of biomass triggered by the experimental removal of Mazzaella. Canopy removal led to a severe decrement in the biomass of macrograzers, while, at the same time, increased the biomass of mesograzers. Asynchrony stemmed from compensatory trophic responses of mesograzers to increased abundances of opportunistic algae. Thus, further work on consumer-resource interactions will improve our understanding of the links between population- and community-level aspects of stability.

Meridionella gen. nov., a New Genus of Cystocloniaceae (Gigartinales, Rhodophyta) from the Southern Hemisphere, Including M. obtusangula comb. nov. and M. antarctica sp. nov

The classification of Cystoclonium obtusangulum has been questioned since the species was first described by Hooker and Harvey as Gracilaria? obtusangula. The objective of this study was to provide the first comprehensive taxonomic analysis of Cystoclonium obtusangulum, based on DNA sequences coupled with morphological observations made on syntype specimens and new collections. Sequence divergences of rbcL, UPA, and COI-5P, and maximum-likelihood phylogenies for rbcL and 18S demonstrated that specimens identified as Cystoclonium obtusangulum represent a clade of two distinct species that are distantly related to the generitype Cystoclonium purpureum. A new genus, Meridionella gen. nov., is proposed for this clade. The two species placed in this new genus were morphologically indistinguishable cryptic species, but have disjunct distributions, with Meridionella obtusangula comb. nov. found from temperate to cold coasts of South America and the Falkland Islands and Meridionella antarctica sp. nov., occurring in Antarctic waters. Vegetative and reproductive characters of Meridionella gen. nov. are described, and implications of our results for the biogeography of the family Cystocloniaceae are discussed.

Comparative phylogeography of two Agarophyton species in the New Zealand archipelago

Molecular studies have reported the coexistence of two species of Agarophyton in New Zealand: the newly described A.transtasmanicum with an apparently restricted distribution to some sites in the North Island, and the more widespread A.chilense. Here, we compared the distribution, genetic diversity, and structure of both Agarophyton species throughout the archipelago using sequences of the nuclear Internal Transcribed Spacer 2 (ITS2) marker. Agarophyton chilense's distribution was continuous and extensive along the North and South Islands, Stewart Island, and Chatham Island, and the genetic clusters were mostly concordant with boundaries between biogeographic regions. In contrast, specimens of A.transtasmanicum were collected in four sites broadly distributed in both the North and South Islands, with no clear spatial structure of the genetic diversity. Populations, where the species co-occurred, tended to display similar levels in genetic diversity for the two species. Demographic inferences supported a postglacial demographic expansion for two A.chilense genetic clusters, one present in the South Island and the eastern coast of the North Island, and the other present in northern South Island. A third genetic cluster located on the western coast of the North Island had a signature of long-term demographic stability. For A.transtasmanicum, the skyline plot also suggested a postglacial demographic expansion. Last, we developed a new molecular tool to quickly and easily distinguish between the two Agarophyton species, which could be used to ease future fine-scale population studies, especially in areas where the two species coexist.

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.

Novel co-occurrence of functionally redundant consumers induced by range expansion alters community structure

Ongoing climate change is shifting the geographic distributions of some species, potentially imposing rapid changes in local community structure and ecosystem functioning. Besides changes in population-level interspecific interactions, such range shifts may also cause changes in functional structure within the host assemblages, which can result in losses or gains in ecosystem functions. Because consumer-resource dynamics are central to community regulation, functional reorganization driven by introduction of new consumer species can have large consequences on ecosystem functions. Here we experimentally examine the extent to which the recent poleward range expansion of the intertidal grazer limpet Scurria viridula along the coast of Chile has altered the role of the resident congeneric limpet S. zebrina, and whether the net collective impacts, and functional structure, of the entire herbivore guild have been modified by the introduction of this new member. We examined the functional role of Scurria species in controlling ephemeral algal cover, bare rock availability, and species richness and diversity, and compared the effects in the region of range overlap against their respective "native" abutted ranges. Experiments showed depression of per capita effects of the range-expanded species within the region of overlap, suggesting environmental conditions negatively affect individual performance. In contrast, effects of S. zebrina were commonly invariant at its range edge. When comparing single species versus polycultures, effects on bare rock cover were altered by the presence of the other Scurria species, suggesting competition between Scurria species. Importantly, although the magnitude of S. viridula effects at the range overlap was reduced, its addition to the herbivore guild seems to complement and intensify the role of the guild in reducing green algal cover, species richness and increasing bare space provision. Our study thus highlights that range expansion of an herbivore can modify the functional guild structure in the recipient community. It also highlights the complexity of predicting how functional structure may change in the face of natural or human-induced range expansions. There is a need for more field-based examination of regional functional compensation, complementarity, or inhibition before we can construct a conceptual framework to anticipate the consequences of species range expansions.

Contrasting patterns of genetic structure and phylogeography in the marine agarophytes Gelidiophycus divaricatus and G. freshwateri (Gelidiales, Rhodophyta) from East Asia

The evolutionary and population demographic history of marine red algae in East Asia is poorly understood. Here, we reconstructed the phylogeographies of two upper intertidal species endemic to East Asia, Gelidiophycus divaricatus and G. freshwateri. Phylogenetic and phylogeographic inferences of 393 mitochondrial cox1, 128 plastid rbcL, and 342 nuclear ITS2 sequences were complemented with ecological niche models. Gelidiophycus divaricatus, a southern species adapted to warm water, is characterized by a high genetic diversity and a strong geographical population structure, characteristic of stable population sizes and sudden reduction to recent expansion. In contrast, G. freshwateri, a northern species adapted to cold temperate conditions, is genetically relatively homogeneous with a shallow population structure resulting from steady population growth and recent equilibrium. The overlap zone of the two species roughly matches summer and winter isotherms, indicating that surface seawater temperature is a key feature influencing species range. Unidirectional genetic introgression was detected at two sites on Jeju Island where G. divaricatus was rare while G. freshwateri was common, suggesting the occurrence of asymmetric natural hybrids, a rarely reported event for rhodophytes. Our results illustrate that Quaternary climate oscillations have left strong imprints on the current day genetic structure and highlight the importance of seawater temperature and sea level change in driving speciation in upper intertidal seaweed species.

Dual influence of terrestrial and marine historical processes on the phylogeography of the Brazilian intertidal red alga Gracilaria caudata

In this study, we explored how past terrestrial and marine climate changes have interacted to shape the phylogeographic patterns of the intertidal red seaweed Gracilaria caudata, an economically important species exploited for agar production in the Brazilian north-east. Seven sites were sampled along the north-east tropical and south-east sub-tropical Brazilian coast. The genetic diversity and structure of G. caudata was inferred using a combination of mitochondrial (COI and cox2-3), chloroplast (rbcL) and 15 nuclear microsatellite markers. A remarkable congruence between nuclear, mitochondrial and chloroplast data revealed clear separation between the north-east (from 03° S to 08° S) and the south-east (from 20° S to 23° S) coast of Brazil. These two clades differ in their demographic histories, with signatures of recent demographic expansions in the north-east and divergent populations in the south-east, suggesting the maintenance of several refugia during the last glacial maximum due to sea-level rise and fall. The Bahia region (around 12° S) occupies an intermediate position between both clades. Microsatellites and mtDNA markers showed additional levels of genetic structure within each sampled site located south of Bahia. The separation between the two main groups in G. caudata is likely recent, probably occurring during the Quaternary glacial cycles. The genetic breaks are concordant with (i) those separating terrestrial refugia, (ii) major river outflows and (iii) frontiers between tropical and subtropical regions. Taken together with previously published eco-physiological studies that showed differences in the physiological performance of the strains from distinct locations, these results suggest that the divergent clades in G. caudata correspond to distinct ecotypes in the process of incipient speciation and thus should be considered for the management policy of this commercially important species.

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.

Genetic and morphological differentiation of Porphyra and Pyropia species (Bangiales, Rhodophyta) coexisting in a rocky intertidal in Central Chile

A recent molecular taxonomic study along the Chilean coast (18° S-53° S) described 18 candidate species of bladed Bangiales of which only two were formally described. Few studies focused on local genetic and morphological diversity of bladed Bangiales and attempted to determine their intertidal distribution in contrasting habitats, and none were performed in Chile. To delimit intertidal distributions of genetic species, 66 samples of bladed Bangiales were collected at Maitencillo (32° S) in four zones: a rocky platform, a rocky wall, and two boulders zones surrounded by sandy and rocky bottoms, respectively. These samples were identified based on sequences of the mitochondrial COI and chloroplast rbcL markers. We also collected 87 specimens for morphological characterization of the most common species, rapidly assessing their putative species identity using newly developed species-diagnostic (PCR-RFLP) markers. Eight microscopic and two macroscopic morphological traits were measured. We described and named three of four species that predominate in Maitencillo (including Pyropia orbicularis): Pyropia variabilis Zapata, Meynard, Ramírez, Contreras-Porcia, sp. nov., Porphyra luchea Meynard, Ramírez, Contreras-Porcia sp. nov., and Porphyra longissima Meynard, Ramírez, Contreras-Porcia, sp. nov. With the exception of Po. longissima restricted to boulders surrounded by sandy bottom, and a morphotype of Py. variabilis restricted to rocky walls, the other species/morphotypes have overlapping intertidal distributions. Except for Po. longissima, which is clearly differentiated morphologically (longest and thinnest blades), we conclude that morphology is not sufficient to differentiate bladed Bangiales. Our findings underscore the importance of refining our knowledge of intrinsic and environmental determinants on the distribution of bladed Bangiales.

Hidden diversity: Phylogeography of genus Ototyphlonemertes Diesing, 1863 (Ototyphlonemertidae: Hoplonemertea) reveals cryptic species and high diversity in Chilean populations

Ototyphlonemertes is a cosmopolitan genus of meiofaunal nemerteans. Their morphological characters are insufficient to reliably identify and delimit species. Consequently, some of the species are considered cosmopolitan despite anticipated low dispersion capability of the adults and a short planktonic larval phase. Indeed, recent studies show that some species actually comprise cryptic species, and populations are connected by stochastic events of long-distance dispersion. Based solely on morphological traits, a Lactea and a Pallida morph of Ototyphlonemertes are recognized here from collections at eight and five locations respectively along the Chilean coast. To assess the phylogeographic patterns of their populations, two mitochondrial markers (COI and COX3) of 162 specimens of Lactea and 25 of Pallida were sequenced. Final sequences are 605bp and 362bp for COI and COX3, respectively. Results from phylogenetic and haplotype network analyses suggest that the Lactea morph comprises up to three independent evolutionary units (one with only COX3 sequences). A COI gene tree including other previously published Ototyphlonemertes sequences groups the Chilean Lactea with other Lactea, while the Chilean Pallida is grouped with other Pallida. Different structuring and gene flow patterns found for the four groups support the hypothesis that these are four independent evolutionary entities with different ecological, dispersal and demographical characteristics.

The variable routes of rafting: stranding dynamics of floating bull kelp Durvillaea antarctica (Fucales, Phaeophyceae) on beaches in the SE Pacific

Dispersal on floating seaweeds depends on availability, viability, and trajectories of the rafts. In the southern hemisphere, the bull kelp Durvillaea antarctica is one of the most common floating seaweeds, but phylogeographic studies had shown low connectivity between populations from continental Chile, which could be due to limitations in local supply and dispersal of floating kelps. To test this hypothesis, the spatiotemporal dynamics of kelp strandings were examined in four biogeographic districts along the Chilean coast (28°-42°S). We determined the biomass and demography of stranded individuals on 33 beaches for three subsequent years (2013, 2014, 2015) to examine whether rafting is restricted to certain districts and seasons (winter or summer). Stranded kelps were found on all beaches. Most kelps had only one stipe (one individual), although we also frequently found coalesced holdfasts with mature males and females, which would facilitate successful rafting dispersal, gamete release, and reproduction upon arrival. High biomasses of stranded kelps occurred in the northern-central (30°S-33°S) and southernmost districts (37°S-42°S), and lower biomasses in the northernmost (28°S-30°S) and southern-central districts (33°S-37°S). The highest percentages and sizes of epibionts (Lepas spp.), indicative of prolonged floating periods, were found on stranded kelps in the northernmost and southernmost districts. Based on these results, we conclude that rafting dispersal can vary regionally, being more common in the northernmost and southernmost districts, depending on intrinsic (seaweed biology) and extrinsic factors (shore morphology and oceanography) that affect local supply of kelps and regional hydrodynamics.

Microsatellite markers and cytoplasmic sequences reveal contrasting pattern of spatial genetic structure in the red algae species complex Mazzaella laminarioides

Mazzaella laminarioides is a common haploid-diploid red alga that forms dense beds. This alga has a wide distributional range, covering 3,500 km of the Chilean coast, but is restricted to high rocky intertidal zones. Recently, the existence of three highly divergent genetic lineages was demonstrated for this taxon, and two cytoplasmic markers were used to determine that these lineages are distributed in strict parapatry. Here, using 454 next-generation sequencing, we developed polymorphic microsatellite loci that cross amplify in all three cytoplasmic lineages. Six sites (i.e., two sites within each lineage) were analyzed using nine microsatellite loci. Our work shows that, although substantial cytoplasmic differentiation occurs within M. laminarioides, the microsatellite loci did not retrieve three nuclear genetic clusters as expected. Indeed, while the northernmost and southernmost cytoplasmic lineages form two strongly divergent nuclear groups characterized by diagnostic alleles, the third cytoplasmic lineage did not form a third nuclear independent group. It is possible that inter-lineage gene exchange has occurred, particularly at sites along the contact zone between the different cytoplasmic lineages. This nuclear-cytoplasmic incongruence in M. laminarioides could be explained by incomplete lineage sorting of the nuclear genes or asymmetric introgressive hybridization between the lineages. Finally, highly significant heterozygote deficiencies (suggesting occurrence of intergametophytic selfing) were observed in the three small northernmost sites while the large southernmost sites generally approached panmixia.

The oceanic concordance of phylogeography and biogeography: a case study in Notochthamalus

Dispersal and adaptation are the two primary mechanisms that set the range distributions for a population or species. As such, understanding how these mechanisms interact in marine organisms in particular - with capacity for long-range dispersal and a poor understanding of what selective environments species are responding to - can provide useful insights for the exploration of biogeographic patterns. Previously, the barnacle Notochthamalus scabrosus has revealed two evolutionarily distinct lineages with a joint distribution that suggests an association with one of the two major biogeographic boundaries (~30°S) along the coast of Chile. However, spatial and genomic sampling of this system has been limited until now. We hypothesized that given the strong oceanographic and environmental shifts associated with the other major biogeographic boundary (~42°S) for Chilean coastal invertebrates, the southern mitochondrial lineage would dominate or go to fixation in locations further to the south. We also evaluated nuclear polymorphism data from 130 single nucleotide polymorphisms to evaluate the concordance of the signal from the nuclear genome with that of the mitochondrial sample. Through the application of standard population genetic approaches along with a Lagrangian ocean connectivity model, we describe the codistribution of these lineages through a simultaneous evaluation of coastal lineage frequencies, an approximation of larval behavior, and current-driven dispersal. Our results show that this pattern could not persist without the two lineages having distinct environmental optima. We suggest that a more thorough integration of larval dynamics, explicit dispersal models, and near-shore environmental analysis can explain much of the coastal biogeography of Chile.

Facilitative Effect of a Generalist Herbivore on the Recovery of a Perennial Alga: Consequences for Persistence at the Edge of Their Geographic Range

Understanding the impacts of consumers on the abundance, growth rate, recovery and persistence of their resources across their distributional range can shed light on the role of trophic interactions in determining species range shifts. Here, we examined if consumptive effects of the intertidal grazer Scurria viridula positively influences the abundance and recovery from disturbances of the alga Mazzaella laminarioides at the edge of its geographic distributions in northern-central Chilean rocky shores. Through field experiments conducted at a site in the region where M. laminarioides overlaps with the polar range edge of S. viridula, we estimated the effects of grazing on different life stages of M. laminarioides. We also used long-term abundance surveys conducted across ~700 km of the shore to evaluate co-occurrence patterns of the study species across their range overlap. We found that S. viridula had positive net effects on M. laminarioides by increasing its cover and re-growth from perennial basal crusts. Probability of occurrence of M. laminarioides increased significantly with increasing density of S. viridula across the range overlap. The negative effect of S. viridula on the percentage cover of opportunistic green algae—shown to compete for space with corticated algae—suggests that competitive release may be part of the mechanism driving the positive effect of the limpet on the abundance and recovery from disturbance of M. laminarioides. We suggest that grazer populations contribute to enhance the abundance of M. laminarioides, facilitating its recolonization and persistence at its distributional range edge. Our study highlights that indirect facilitation can determine the recovery and persistence of a resource at the limit of its distribution, and may well contribute to the ecological mechanisms governing species distributions and range shifts.

The bladed Bangiales (Rhodophyta) of the South Eastern Pacific: Molecular species delimitation reveals extensive diversity

A molecular taxonomic study of the bladed Bangiales of the South Eastern Pacific (coast of Chile) was undertaken based on sequence data of the mitochondrial COI and chloroplast rbcL for 193 specimens collected from Arica (18°S) in the north to South Patagonia (53°S) in the south. The results revealed for the first time that four genera, Porphyra, Pyropia, Fuscifolium and Wildemania were present in the region. Species delimitation was determined based on a combination of a General Mixed Yule Coalescence model (GMYC) and Automatic Barcode Gap Discovery (ABGD) coupled with detection of monophyly in tree reconstruction. The overall incongruence between the species delimitation methods within each gene was 29%. The GMYC method led to over-splitting groups, whereas the ABGD method had a tendency to lump groups. Taking a conservative approach to the number of putative species, at least 18 were recognized and, with the exception of the recently described Pyropia orbicularis, all were new to the Chilean flora. Porphyra and Pyropia were the most diverse genera with eight 'species' each, whereas only a 'single' species each was found for Fuscifolium and Wildemania. There was also evidence of recently diverging groups: Wildemania sp. was distinct but very closely related to W. amplissima from the Northern Hemisphere and raises questions in relation to such disjunct distributions. Pyropia orbicularis was very closely related to two other species, making species delimitation very difficult but provides evidence of an incipient speciation. The difference between the 'species' discovered and those previously reported for the region is discussed in relation to the difficulty of distinguishing species based on morphological identification.

Genes Left Behind: Climate Change Threatens Cryptic Genetic Diversity in the Canopy-Forming Seaweed Bifurcaria bifurcata

The global redistribution of biodiversity will intensify in the coming decades of climate change, making projections of species range shifts and of associated genetic losses important components of conservation planning. Highly-structured marine species, notably brown seaweeds, often harbor unique genetic variation at warmer low-latitude rear edges and thus are of particular concern. Here, a combination of Ecological Niche Models (ENMs) and molecular data is used to forecast the potential near-future impacts of climate change for a warm-temperate, canopy forming seaweed, Bifurcaria bifurcata. ENMs for B. bifurcata were developed using marine and terrestrial climatic variables, and its range projected for 2040-50 and 2090-2100 under two greenhouse emission scenarios. Geographical patterns of genetic diversity were assessed by screening 18 populations spawning the entire distribution for two organelle genes and 6 microsatellite markers. The southern limit of B. bifurcata was predicted to shift northwards to central Morocco by the mid-century. By 2090-2100, depending on the emission scenario, it could either retreat further north to western Iberia or be relocated back to Western Sahara. At the opposing margin, B. bifurcata was predicted to expand its range to Scotland or even Norway. Microsatellite diversity and endemism were highest in Morocco, where a unique and very restricted lineage was also identified. Our results imply that B. bifurcata will maintain a relatively broad latitudinal distribution. Although its persistence is not threatened, the predicted extirpation of a unique southern lineage or even the entire Moroccan diversity hotspot will erase a rich evolutionary legacy and shrink global diversity to current (low) European levels. NW Africa and similarly understudied southern regions should receive added attention if expected range changes and diversity loss of warm-temperate species is not to occur unnoticed.

Contrasting patterns of population structure and demographic history in cryptic species of Bostrychia intricata (Rhodomelaceae, Rhodophyta) from New Zealand

Spatial patterns of genetic diversity provide insight into the demography and history of species. Morphologically similar but genetically distinct "cryptic" species are increasingly being recognized in marine organisms through molecular analyses. Such species are, on closer inspection, often discovered to display contrasting life histories or occasionally minor morphological differences; molecular tools can thus be useful indicators of diversity. Bostrychia intricata, a marine red alga, is widely distributed throughout the Southern Hemisphere and comprises many cryptic species. We used mitochondrial cytochrome c oxidase I gene sequences to assess the genetic variation, population genetic structure, and demographic history of B. intricata in New Zealand. Our results supported the existence of three cryptic species of B. intricata (N2, N4, and N5) in New Zealand. Cryptic species N4, which was found throughout New Zealand, showed a higher genetic diversity and wider distribution than the other two species, which were only found in the North Island and northern South Island. Our analyses showed low to moderate genetic differentiation among eastern North Island populations for cryptic species N2, but high differentiation among North and South Island populations for N4, suggesting different population structure between these cryptic species. Data also indicated that N2 has recently undergone population expansion, probably since the Last Glacial Maximum (LGM), while the higher genetic diversity in N4 populations suggests persistence in situ through the LGM. The contrasting population structures and inferred demographic histories of these species highlight that life history can vary greatly even among morphologically indistinguishable taxa.

Latitudinal discontinuity in thermal conditions along the nearshore of central-northern Chile

Over the past decade, evidence of abrupt latitudinal changes in the dynamics, structure and genetic variability of intertidal and subtidal benthic communities along central-northern Chile has been found consistently at 30–32°S. Changes in the advective and thermal environment in nearshore waters have been inferred from ecological patterns, since analyses of in situ physical data have thus far been missing. Here we analyze a unique set of shoreline temperature data, gathered over 4–10 years at 15 sites between 28–35°S, and combine it with satellite-derived winds and sea surface temperatures to investigate the latitudinal transition in nearshore oceanographic conditions suggested by recent ecological studies. Our results show a marked transition in thermal conditions at 30–31°S, superimposed on a broad latitudinal trend, and small-scale structures associated with cape-and-bay topography. The seasonal cycle dominated temperature variability throughout the region, but its relative importance decreased abruptly south of 30–31°S, as variability at synoptic and intra-seasonal scales became more important. The response of shoreline temperatures to meridional wind stress also changed abruptly at the transition, leading to a sharp drop in the occurrence of low-temperature waters at northern sites, and a concurrent decrease in corticated algal biomass. Together, these results suggest a limitation of nitrate availability in nearshore waters north of the transition. The localized alongshore change results from the interaction of latitudinal trends (e.g., wind stress, surface warming, inertial period) with a major headland-bay system (Punta Lengua de Vaca at 30.25°S), which juxtaposes a southern stretch of coast characterized by upwelling with a northern stretch of coast characterized by warm surface waters and stratification. This transition likely generates a number of latitude-dependent controls on ecological processes in the nearshore that can explain species-specific effects, and add strength to the suggestion of an oceanography-driven, major spatial transition in coastal communities at 30–31°S.

Scorched mussels (BIVALVIA: MYTILIDAE: BRACHIDONTINAE) from the temperate coasts of South America: phylogenetic relationships, trans-Pacific connections and the footprints of Quaternary glaciations

This study addresses aspects of the phylogeny and phylogeography of scorched mussels (BIVALVIA: MYTILIDAE: BRACHIDONTINAE) from southern South America (Argentina and Chile), as well as their ecophylogenetic implications. Relationships were inferred from sequences of two nuclear (28S and 18S) and one mitochondrial (COI) genes, using Bayesian and maximum likelihood analyses. Our results indicate that the monophyletic BRACHIDONTINAE include three well supported clades: [i] Brachidontes Swainson (=Hormomya Mörch), [ii] Ischadium Jukes-Browne+Geukensia van de Poel, and [iii] Austromytilus Laseron+Mytilisepta Habe (usually considered a member of the SEPTIFERINAE)+Perumytilus Olsson. Species of clade [iii] are distributed along the temperate coasts of the Pacific Ocean. Available evidence supports divergence between Austromytilus (Australia) and Perumytilus (South American) following the breakup of Australian, Antarctic and South American shelves. Four brachidontins occur in southern South America: Brachidontes rodriguezii (d'Orbigny), B. granulatus (Hanley), and two genetically distinct clades of Perumytilus. The latter are confined to the Chile-Peru (North Clade) and Magellanic (South Clade) Biogeographic Provinces, respectively warm- and cold-temperate. The South Clade is the only brachidontin restricted to cold-temperate waters. Biogeographic considerations and the fossil record prompted the hypothesis that the South Clade originated from the North Clade by incipient peripatric differentiation, followed by isolation during the Quaternary glaciations, genetic differentiation in the non-glaciated coasts of eastern Patagonia, back-expansion over southern Chile following post-LGM de-glaciation, and development of a secondary contact zone between the two clades in south-central Chile. Evidence of upper Pleistocene expansion of the South Clade parallels similar results on other organisms that have colonized coastal ecosystems from eastern Patagonia since the LGM, apparently occupying free ecological space. We emphasize that the assembly of communities cannot be explained solely in terms of environmental drivers, as history also matters.

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.

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

Background

Pyura chilensis, a tunicate commercially exploited as food resource in Chile, is subject to management strategies, including restocking. The goal of this study was to examine the genetic structure of P. chilensis using information from a mitochondrial gene (Cytochrome Oxidase I, COI) and a nuclear gene (Elongation 1 alpha, EF1a), to characterize the geographic distribution of genetic diversity and differentiation, and to identify the main processes that have shaped it. We analyzed 268 and 208 sequences of COI and EF1a, respectively, from samples of eight local populations covering ca. 1800 km.

Results

For Pyura chilensis, partial sequences of the gene COI revealed three highly supported haplogroups that diverged 260000 to 470000 years ago. Two haplogroups currently are widely distributed and sympatric, while one is dominant only in Los Molinos (LM, 39°50′S). The two widespread COI haplogroups underwent a geographic expansion during an interglacial period of the Late Pleistocene ca. 100000 years ago. The nuclear gene was less divergent and did not resolve the COI haplogroups. Bayesian clustering of the nuclear gene’s SNPs revealed that individuals from the two widespread COI haplogroups were mostly assigned to two of the three detected clusters and had a marked degree of admixture. The third cluster predominated in LM and showed low admixture. Haplotypic diversity of both genes was very high, there was no isolation by distance, and most localities were genetically undifferentiated; only LM was consistently differentiated with both genes analyzed.

Conclusions

Pyura chilensis has less genetic structure than expected given its life history, which could be a consequence of dispersal on ship hulls. The only differentiated local population analyzed was LM. Coincidentally, it is the one furthest away from main maritime routes along the coast of Chile.

The use of mitochondrial and nuclear markers allowed detection of divergent mitochondrial haplogroups in P. chilensis, two of which revealed nuclear admixture. The genetic structure of P. chilensis has likely been shaped by Pleistocene’s climatic effect on sea level leading to population contraction with isolation, followed by geographic range expansions with concomitant secondary contact and admixture.