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.

Phylogeographic analyses of the 30 degrees S south-east Pacific biogeographic transition zone establish the occurrence of a sharp genetic discontinuity in the kelp Lessonia nigrescens: vicariance or parapatry?

Phylogeographic studies are lacking in the Southern Hemisphere, and in particular in the south-eastern Pacific. To infer the possible scenario for the debated biogeographic transition zone located at 30-33 degrees S along the Chilean coast, we investigated whether there is a concordance between the phylogeographic pattern and the biogeographic transition in the intertidal kelp Lessonia nigrescens whose distribution is continuous across this transition zone. Using a combination of four makers located in the three genomic compartments (chloroplast, mitochondria and nucleus), we showed the presence of two main divergent lineages, possibly cryptic species. There was an exact match of the phylogeographic break with the 30 degrees S biogeographic transition zone, suggesting a common origin. The combined information given by the multilocus approach and by the population analysis suggested the occurrence of a budding speciation, with a northward range expansion.

Phospholipases and Galactolipases Trigger Oxylipin-Mediated Wound-Activated Defence in the Red Alga Gracilaria chilensis against Epiphytes

We investigated the wound response of the commercially important red alga, Gracilaria chilensis, in order to obtain insight into its interaction with epiphytic pests. After wounding, the host releases free fatty acids as well as the hydroxylated eicosanoids, 8R-hydroxy eicosatetraenoic acid (8-HETE) and 7S,8R-dihydroxy eicosatetraenoic acid (7,8-di-HETE). While the release of free arachidonic acid and subsequent formation of 8-HETE is controlled by phospholipase A, 7,8-di-HETE production is independent of this lipase. This dihydroxylated fatty acid might be directly released from galactolipids. Physiologically relevant concentrations of oxylipins reduced spore settlement of Acrochaetium sp. (Rhodophyta, Acrochaetiaceae) and suppressed the development of hapteria in Ceramium rubrum (Rhodophyta, Ceramiaceae) when these model epiphytes were exposed to artificial surfaces that contained 8-HETE or 7,8-di-HETE. Thus, the immediate release of oxylipins can be seen as G. chilensis defence against epiphytes.

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

Background

The Chilean shoreline, a nearly strait line of coast expanding across 35 latitudinal degrees, represents an interesting region to assess historical processes using phylogeographic analyses. Stretching along the temperate section of the East Pacific margin, the region is characterized by intense geologic activity and has experienced drastic geomorphological transformations linked to eustatic and isostatic changes during the Quaternary. In this study, we used two molecular markers to evaluate the existence of phylogeographic discontinuities and detect the genetic footprints of Pleistocene glaciations among Patagonian populations of Mazzaella laminarioides, a low-dispersal benthic intertidal red seaweed that inhabits along ~3,700 km of the Chilean coastal rocky shore.

Results

Three main genetic lineages were found within M. laminarioides. They are distributed along the Chilean coast in strict parapatry. The deep divergence among lineages suggests that they could be considered putative genetic sibling species. Unexpectedly, genetic breaks were not strictly concordant with the biogeographic breaks described in the region. A Northern lineage was restricted to a broad transition zone located between 30°S and 33°S and showed signals of a recent bottleneck. The reduction of population size could be related to warm events linked to El Niño Southern Oscillation, which is known to cause massive seaweed mortality in this region. To the south, we propose that transient habitat discontinuities driven by episodic tectonic uplifting of the shoreline around the Arauco region (37°S-38°S); one of the most active forearc-basins in the South East Pacific; could be at the origin of the Central/South genetic break. The large beaches, located around 38°S, are likely to contribute to the lineages’ integrity by limiting present gene flow. Finally, the Southern lineage, occupies an area affected by ice-cover during the last glaciations. Phylogeny suggested it is a derived clade and demographic analyses showed the lineage has a typical signature of postglacial recolonization from a northern glacial refugium area.

Conclusions

Even if environmental adaptation could have strengthened divergence among lineages in M. laminarioides, low dispersal capacity and small population size are sufficient to generate phylogeographic discontinuities determined by genetic drift alone. Interestingly, our results confirm that seaweed population connectivity over large geographic scales does not rely only on dispersal capacity but also seem to depend highly on substratum availability and population density of the receiving locality.

Self-fertilization is the main sexual reproduction mechanism in native wine yeast populations

Saccharomyces cerevisiae is a model eukaryotic organism for classical genetics and genomics, and yet its ecology is still largely unknown. In this work, a population genetic analysis was performed on five yeast populations isolated from wine-making areas with different enological practices using simple sequence repeats and restriction fragment length polymorphism of mitochondrial DNA as molecular markers on 292 strains. In accordance with other studies, genome size estimation suggests that native S. cerevisiae strains are mainly homothallic and diploids. Analysis of mtDNA data showed that yeast populations from nonindustrial areas have 40% higher genetic diversity than populations isolated from industrial areas, demonstrating that industrial enological practices are likely to affect native yeast populations negatively by reducing its biodiversity. On the other hand, genetic differentiation analysis based on their microsatellite showed no correlation between genetic and geographic distance and a nonsignificant value when a Mantel test was applied. Finally, in the five populations studied, positive inbreeding (F(is)) values from 0.4 to 0.75, a low but significant level of linkage disequilibrium and a high number of multilocus genotypes were detected. These results strongly advocate that sexual reproduction is frequent enough to erase clonal signature in natural populations and that self-fertilization is the main mating system.

A fine-tuned global distribution dataset of marine forests

Species distribution records are a prerequisite to follow climate-induced range shifts across space and time. However, synthesizing information from various sources such as peer-reviewed literature, herbaria, digital repositories and citizen science initiatives is not only costly and time consuming, but also challenging, as data may contain thematic and taxonomic errors and generally lack standardized formats. We address this gap for important marine ecosystem-structuring species of large brown algae and seagrasses. We gathered distribution records from various sources and provide a fine-tuned dataset with ~2.8 million dereplicated records, taxonomically standardized for 682 species, and considering important physiological and biogeographical traits. Specifically, a flagging system was implemented to signal potentially incorrect records reported on land, in regions with limiting light conditions for photosynthesis, and outside the known distribution of species, as inferred from the most recent published literature. We document the procedure and provide a dataset in tabular format based on Darwin Core Standard (DwC), alongside with a set of functions in R language for data management and visualization.

Phylogeography of the Marine Otter (Lontra felina): historical and contemporary factors determining its distribution

The evolutionary history of a species can be revealed by phylogeographical analysis; nevertheless, not only historical but also contemporary processes can imprint on the distribution of genetic diversity. We report on the phylogeny of Lontra ssp. in South America, and the role of spatial heterogeneity in shaping the distribution and population structure of the endangered marine otter, Lontra felina. Analyzing a total of 2261 bp of mitochondrial DNA (mtDNA) revealed the recent divergence of L. felina from L. provocax. A strong population structure (Φ(st) = 0.83, P < 0.0001) and a significant pattern of isolation by distance were described for L. felina (n = 168) across a wide geographical distribution (13°53'S to 43°36'S). Lontra felina mtDNA phylogeny is composed of 2 main clades: a clade from Peru and another composed of Chilean haplotypes. Northern populations show different divergent lineages and higher genetic diversity when compared with more recently colonized southern populations. Furthermore, long sandy beaches seem to act as barriers to dispersal, creating 2 evolutionary significant units in agreement with subspecies previous description, and at least 5 different management units (MUs). At a fine spatial scale, the size of rocky seashore patches, the distance between patches and anthropogenic factors also play important roles in species gene flow.

High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga Ectocarpus

Deciphering the genetic architecture of adaptation of brown algae to environmental stresses such as temperature and salinity is of evolutionary as well as of practical interest. The filamentous brown alga Ectocarpus sp. is a model for the brown algae and its genome has been sequenced. As sessile organisms, brown algae need to be capable of resisting the various abiotic stressors that act in the intertidal zone (e.g. osmotic pressure, temperature, salinity, UV radiation) and previous studies have shown that an important proportion of the expressed genes is regulated in response to hyposaline, hypersaline or oxidative stress conditions. Using the double digest RAD sequencing method, we constructed a dense genetic map with 3,588 SNP markers and identified 39 QTLs for growth-related traits and their plasticity under different temperature and salinity conditions (tolerance to high temperature and low salinity). GO enrichment tests within QTL intervals highlighted membrane transport processes such as ion transporters. Our study represents a significant step towards deciphering the genetic basis of adaptation of Ectocarpus sp. to stress conditions and provides a substantial resource to the increasing list of tools generated for the species.

THE LESSONIA NIGRESCENS SPECIES COMPLEX (LAMINARIALES, PHAEOPHYCEAE) SHOWS STRICT PARAPATRY AND COMPLETE REPRODUCTIVE ISOLATION IN A SECONDARY CONTACT ZONE(1)

During secondary contact between phylogenetically closely related species (sibling species) having diverged in allopatry, the maintenance of species integrity depends on intrinsic and extrinsic reproductive barriers. In kelps (Phaeophyceae), the observations of hybrids in laboratory conditions suggest that reproductive isolation is incomplete. However, not all interspecific crosses are successful, and very few hybrids have been observed in nature, despite the co-occurrence of many kelp species in sympatry. This suggests that there are reproductive barriers that maintain species integrity. In this study, we characterized the fine genetic structure of a secondary contact zone to clarify the extent of reproductive isolation between two sister species. In Lessonia nigrescens Bory (Laminariales, Phaeophyta) species complex, two cryptic species have been recently found out from gene phylogenies, and-waiting for a formal taxonomic description-we used their geographic distribution to name them (northern and southern species). We studied 12 populations, distributed along 50 km of coastline, and employed two molecular approaches, assigning individuals to phylogenetic species according to a diagnostic mitochondrial marker (351 individuals analyzed) and quantifying interspecific gene flow with four microsatellite markers (248 individuals analyzed). No hybridization or introgression was revealed, indicating complete reproductive isolation in natural conditions. Unexpectedly, our study demonstrated that the two species were strictly segregated in space. This absence of co-occurrence along the contact zone can partially explain the lack of hybridization, raising new interesting questions as to the mechanisms that limit sympatry at small spatial scales.

SEX RATIO VARIATION IN THE LESSONIA NIGRESCENS COMPLEX (LAMINARIALES, PHAEOPHYCEAE): EFFECT OF LATITUDE, TEMPERATURE, AND MARGINALITY(1)

Little is known about variation of sex ratio, the proportion of males to females, in natural populations of seaweed, though it is a major determinant of the mating system. The observation of sexual chromosomes in kelps suggested that sex is partly genetically determined. However, it is probably not purely genetic since the sex ratio can be modified by environmental factors such as salinity or temperature. In this paper, sex ratio variation was studied in the kelp Lessonia nigrescens Bory complex, recently identified as two cryptic species occurring along the Chilean coast: one located north and the other south of the biogeographic boundary at latitude 29°-30° S. The life cycle of L. nigrescens is characterized by an alternation of microscopic haploid gametophytic individuals and large macroscopic fronds of diploid sporophytes. The sex ratio was recorded in progenies from 241 sporophytic individuals collected from 13 populations distributed along the Chilean coast in order (i) to examine the effect of an environmental gradient coupled with latitude, and (ii) to compare marginal populations to central populations of the two species. In addition, we tested the hypothesis that the sex ratios of the two cryptic species would be affected differently by temperature. First, our results demonstrate that sex ratio seems to be mainly genetically determined and temperature can significantly modify it. Populations of the northern species showed a lower frequency of males at 14°C than at 10°C, whereas populations of the southern species showed the opposite pattern. Second, both species displayed an increased variation in sex ratio at the range limits. This greater variation at the margins could be due either to differential mortality between sexes or to geographic parthenogenesis (asexual reproduction).

Development of PCR-Based Markers to Determine the Sex of Kelps

Sex discriminating genetic markers are commonly used to facilitate breeding programs in economically and ecologically important animal and plant species. However, despite their considerable economic and ecological value, the development of sex markers for kelp species has been very limited. In this study, we used the recently described sequence of the sex determining region (SDR) of the brown algal model Ectocarpus to develop novel DNA-based sex-markers for three commercially relevant kelps: Laminaria digitata, Undaria pinnatifida and Macrocystis pyrifera. Markers were designed within nine protein coding genes of Ectocarpus male and female (U/V) sex chromosomes and tested on gametophytes of the three kelp species. Seven primer pairs corresponding to three loci in the Ectocarpus SDR amplified sex-specific bands in the three kelp species, yielding at least one male and one female marker for each species. Our work has generated the first male sex-specific markers for L. digitata and U. pinnatifida, as well as the first sex markers developed for the genus Macrocystis. The markers and methodology presented here will not only facilitate seaweed breeding programs but also represent useful tools for population and demography studies and provide a means to investigate the evolution of sex determination across this largely understudied eukaryotic group.

Association of the estrogen receptor alpha gene polymorphisms with osteoporosis in the Mexican population

The estrogen receptor gene (ER alpha) has been implicated in the development of osteoporosis. In this study, the association of two ER alpha gene polymorphic markers (a TA dinucleotide repeat and a single nucleotide polymorphism, G2014A) with osteoporosis was tested in 70 osteoporotic women, 70 non-osteoporotic women and 500 subjects from the Mexican population. According to the genetic analysis of the Mexican population using eight unlinked polymorphic markers, we found that our population is structured into three subpopulations; therefore, the allele-phenotype relationship was analyzed with a statistical method that considered population stratification. We found that the G2014A polymorphism is associated with the presence of osteoporosis while the TA dinucleotide repeat is not. The G allele and the GG genotype frequencies of the G2014A marker were significantly higher in osteoporotic than in non-osteoporotic women. Likewise, subjects bearing the G allele in heterozygous or homozygous displayed lower values for lumbar bone mineral density and T score than those who did not present any G allele. The effect of confounders for osteoporosis on the association of G allele-osteoporosis was ruled out. In summary, we conclude that the G2014 polymorphism may become a useful marker for genetic studies of osteoporosis in the Mexican population.

Temperature effects on gametophyte life-history traits and geographic distribution of two cryptic kelp species

A major determinant of the geographic distribution of a species is expected to be its physiological response to changing abiotic variables over its range. The range of a species often corresponds to the geographic extent of temperature regimes the organism can physiologically tolerate. Many species have very distinct life history stages that may exhibit different responses to environmental factors. In this study we emphasized the critical role of the haploid microscopic stage (gametophyte) of the life cycle to explain the difference of edge distribution of two related kelp species. Lessonia nigrescens was recently identified as two cryptic species occurring in parapatry along the Chilean coast: one located north and the other south of a biogeographic boundary at latitude 29-30°S. Six life history traits from microscopic stages were identified and estimated under five treatments of temperature in eight locations distributed along the Chilean coast in order to (1) estimate the role of temperature in the present distribution of the two cryptic L. nigrescens species, (2) compare marginal populations to central populations of the two cryptic species. In addition, we created a periodic matrix model to estimate the population growth rate (λ) at the five temperature treatments. Differential tolerance to temperature was demonstrated between the two species, with the gametophytes of the Northern species being more tolerant to higher temperatures than gametophytes from the south. Second, the two species exhibited different life history strategies with a shorter haploid phase in the Northern species contrasted with considerable vegetative growth in the Southern species haploid stage. These results provide strong ecological evidence for the differentiation process of the two cryptic species and show local adaptation of the life cycle at the range limits of the distribution. Ecological and evolutionary implications of these findings are discussed.

Herbivore-induced chemical and molecular responses of the kelps Laminaria digitata and Lessonia spicata

Kelps are founding species of temperate marine ecosystems, living in intertidal coastal areas where they are often challenged by generalist and specialist herbivores. As most sessile organisms, kelps develop defensive strategies to restrain grazing damage and preserve their own fitness during interactions with herbivores. To decipher some inducible defense and signaling mechanisms, we carried out metabolome and transcriptome analyses in two emblematic kelp species, Lessonia spicata from South Pacific coasts and Laminaria digitata from North Atlantic, when challenged with their main specialist herbivores. Mass spectrometry based metabolomics revealed large metabolic changes induced in these two brown algae following challenges with their own specialist herbivores. Targeted metabolic profiling of L. spicata further showed that free fatty acid (FFA) and amino acid (AA) metabolisms were particularly regulated under grazing. An early stress response was illustrated by the accumulation of Sulphur containing amino acids in the first twelve hours of herbivory pressure. At latter time periods (after 24 hours), we observed FFA liberation and eicosanoid oxylipins synthesis likely representing metabolites related to stress. Global transcriptomic analysis identified sets of candidate genes specifically induced by grazing in both kelps. qPCR analysis of the top candidate genes during a 48-hours time course validated the results. Most of these genes were particularly activated by herbivore challenge after 24 hours, suggesting that transcriptional reprogramming could be operated at this time period. We demonstrated the potential utility of these genes as molecular markers for herbivory by measuring their inductions in grazed individuals of field harvested L. digitata and L. spicata. By unravelling the regulation of some metabolites and genes following grazing pressure in two kelps representative of the two hemispheres, this work contributes to provide a set of herbivore-induced chemical and molecular responses in kelp species, showing similar inducible responses upon specialist herbivores in their respective ecosystems.

Biodiversity of rocky intertidal benthic communities associated with copper mine tailing discharges in northern Chile

Copper mine tailings have been discharged around the city of Chanaral, in northern Chile, for more than 60 years. This report summarizes a 17-month long monitoring study of species richness and biodiversity at five intertidal sites around the point of the tailing discharge. Total dissolved copper in sites close to the point of discharge varied between 8.72 microg/l and 34.15 microg/l, showing that there has not been a significant reduction since 1994. However, species richness has increased, suggesting a possible recovery of the system. While diversity of sessile organisms correlates negatively with dissolved copper, diversity of mobile invertebrates did not correlate with the metal concentration. To explain the observed results we discuss the role of algal turf interference on the distribution of mobile invertebrates at reference sites, a top-down effect caused by the absence of carnivores at impacted sites, and an avoidance strategy by some species to reduce their contact with contaminated seawater.

Mates Matter: Gametophyte Kinship Recognition and Inbreeding in the Giant Kelp, Macrocystis pyrifera (Laminariales, Phaeophyceae)

Inbreeding, the mating between genetically related individuals, often results in reduced survival and fecundity of offspring, relative to outcrossing. Yet, high inbreeding rates are commonly observed in seaweeds, suggesting compensatory reproductive traits may affect the costs and benefits of the mating system. We experimentally manipulated inbreeding levels in controlled crossing experiments, using gametophytes from 19 populations of Macrocystis pyrifera along its Eastern Pacific coastal distribution (EPC). The objective was to investigate the effects of male-female kinship on female fecundity and fertility, to estimate inbreeding depression in the F1 progeny, and to assess the variability of these effects among different regions and habitats of the EPC. Results revealed that the presence and kinship of males had a significant effect on fecundity and fertility of female gametophytes. Females left alone or in the presence of sibling males express the highest gametophyte size, number, and size of oogonia, suggesting they were able to sense the presence and the identity of their mates before gamete contact. The opposite trend was observed for the production of embryos per female gametes, indicating higher costs of selfing and parthenogenesis than outcrossing on fertility. However, the increased fecundity compensated for the reduced fertility, leading to a stable overall reproductive output. Inbreeding also affected morphological traits of juvenile sporophytes, but not their heatwave tolerance. The male-female kinship effect was stronger in high-latitude populations, suggesting that females from low-latitude marginal populations might have evolved to mate with any male gamete to guarantee reproductive success.

Past climate-driven range shifts structuring intraspecific biodiversity levels of the giant kelp (Macrocystis pyrifera) at global scales

The paradigm of past climate-driven range shifts structuring the distribution of marine intraspecific biodiversity lacks replication in biological models exposed to comparable limiting conditions in independent regions. This may lead to confounding effects unlinked to climate drivers. We aim to fill in this gap by asking whether the global distribution of intraspecific biodiversity of giant kelp (Macrocystis pyrifera) is explained by past climate changes occurring across the two hemispheres. We compared the species' population genetic diversity and structure inferred with microsatellite markers, with range shifts and long-term refugial regions predicted with species distribution modelling (SDM) from the last glacial maximum (LGM) to the present. The broad antitropical distribution of Macrocystis pyrifera is composed by six significantly differentiated genetic groups, for which current genetic diversity levels match the expectations of past climate changes. Range shifts from the LGM to the present structured low latitude refugial regions where genetic relics with higher and unique diversity were found (particularly in the Channel Islands of California and in Peru), while post-glacial expansions following ~ 40% range contraction explained extensive regions with homogenous reduced diversity. The estimated effect of past climate-driven range shifts was comparable between hemispheres, largely demonstrating that the distribution of intraspecific marine biodiversity can be structured by comparable evolutionary forces across the global ocean. Additionally, the differentiation and endemicity of regional genetic groups, confers high conservation value to these localized intraspecific biodiversity hotspots of giant kelp forests.