A scaffolded and annotated reference genome of giant kelp (Macrocystis pyrifera).

Macrocystis pyrifera (giant kelp), is a brown macroalga of great ecological importance as a primary producer and structure-forming foundational species that provides habitat for hundreds of species. It has many commercial uses (e.g. source of alginate, fertilizer, cosmetics, feedstock). One of the limitations to exploiting giant kelp's economic potential and assisting in giant kelp conservation efforts is a lack of genomic tools like a high quality, contiguous reference genome with accurate gene annotations. Reference genomes attempt to capture the complete genomic sequence of an individual or species, and importantly provide a universal structure for comparison across a multitude of genetic experiments, both within and between species. We assembled the giant kelp genome of a haploid female gametophyte de novo using PacBio reads, then ordered contigs into chromosome level scaffolds using Hi-C. We found the giant kelp genome to be 537 MB, with a total of 35 scaffolds and 188 contigs. The assembly N50 is 13,669,674 with GC content of 50.37%. We assessed the genome completeness using BUSCO, and found giant kelp contained 94% of the BUSCO genes from the stramenopile clade. Annotation of the giant kelp genome revealed 25,919 genes. Additionally, we present genetic variation data based on 48 diploid giant kelp sporophytes from three different Southern California populations that confirms the population structure found in other studies of these populations. This work resulted in a high-quality giant kelp genome that greatly increases the genetic knowledge of this ecologically and economically vital species.

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.

Comparing and Evaluating Metagenome Assembly Tools from a Microbiologist's Perspective - Not Only Size Matters!

With the constant improvement in cost-efficiency and quality of Next Generation Sequencing technologies, shotgun-sequencing approaches -such as metagenomics- have nowadays become the methods of choice for studying and classifying microorganisms from various habitats. The production of data has dramatically increased over the past years and processing and analysis steps are becoming more and more of a bottleneck. Limiting factors are partly the availability of computational resources, but mainly the bioinformatics expertise in establishing and applying appropriate processing and analysis pipelines. Fortunately, a large diversity of specialized software tools is nowadays available. Nevertheless, choosing the most appropriate methods for answering specific biological questions can be rather challenging, especially for non-bioinformaticians. In order to provide a comprehensive overview and guide for the microbiological scientific community, we assessed the most common and freely available metagenome assembly tools with respect to their output statistics, their sensitivity for low abundant community members and variability in resulting community profiles as well as their ease-of-use. In contrast to the highly anticipated "Critical Assessment of Metagenomic Interpretation" (CAMI) challenge, which uses general mock community-based assembler comparison we here tested assemblers on real Illumina metagenome sequencing data from natural communities of varying complexity sampled from forest soil and algal biofilms. Our observations clearly demonstrate that different assembly tools can prove optimal, depending on the sample type, available computational resources and, most importantly, the specific research goal. In addition, we present detailed descriptions of the underlying principles and pitfalls of publically available assembly tools from a microbiologist's perspective, and provide guidance regarding the user-friendliness, sensitivity and reliability of the resulting phylogenetic profiles.

Nuclear DNA level and life cycle of kelps: Evidence for sex-specific polyteny in Macrocystis (Laminariales, Phaeophyceae).

Giant kelp, Macrocystis pyrifera (Linnaeus) C. Agardh, is the subject of intense breeding studies for marine biomass production and conservation of natural resources. In this context, six gametophyte pairs and a sporophyte offspring of Macrocystis from South America were analyzed by flow cytometry. Minimum relative DNA content per cell (1C) was found in five males. Unexpectedly, nuclei of all female gametophytes contained approximately double the DNA content (2C) of males; the male gametophyte from one locality also contained 2C, likely a spontaneous natural diploid variant. The results illustrate a sex-specific difference in nuclear DNA content among Macrocystis gametophytes, with the chromosomes of the females in a polytenic condition. This correlates with significantly larger cell sizes in female gametophytes compared to males and resource allocation in oogamous reproduction. The results provide key information for the interpretation of DNA measurements in kelp life cycle stages and prompt further research on the regulation of the cell cycle, metabolic activity, sex determination, and sporophyte development.

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.

Seascape drivers of Macrocystis pyrifera population genetic structure in the northeast Pacific.

At small spatial and temporal scales, genetic differentiation is largely controlled by constraints on gene flow, while genetic diversity across a species' distribution is shaped on longer temporal and spatial scales. We assess the hypothesis that oceanographic transport and other seascape features explain different scales of genetic structure of giant kelp, Macrocystis pyrifera. We followed a hierarchical approach to perform a microsatellite-based analysis of genetic differentiation in Macrocystis across its distribution in the northeast Pacific. We used seascape genetic approaches to identify large-scale biogeographic population clusters and investigate whether they could be explained by oceanographic transport and other environmental drivers. We then modelled population genetic differentiation within clusters as a function of oceanographic transport and other environmental factors. Five geographic clusters were identified: Alaska/Canada, central California, continental Santa Barbara, California Channel Islands and mainland southern California/Baja California peninsula. The strongest break occurred between central and southern California, with mainland Santa Barbara sites forming a transition zone between the two. Breaks between clusters corresponded approximately to previously identified biogeographic breaks, but were not solely explained by oceanographic transport. An isolation-by-environment (IBE) pattern was observed where the northern and southern Channel Islands clustered together, but not with closer mainland sites, despite the greater distance between them. The strongest environmental association with this IBE pattern was observed with light extinction coefficient, which extends suitable habitat to deeper areas. Within clusters, we found support for previous results showing that oceanographic connectivity plays an important role in the population genetic structure of Macrocystis in the Northern hemisphere.

Genetic and experimental evidence for a mixed-age, mixed-origin bank of kelp microscopic stages in southern California.

Laboratory studies have demonstrated that the microscopic stages of kelps can rapidly resume development from a delayed state. Like terrestrial seeds or aquatic resting eggs, banks of delayed kelp stages may supplement population recovery after periods of stress, playing an important role for kelp populations that experience adult sporophyte absences due to seasonal or interannual disturbances. We found that removing the microscopic stages from natural rock substratum could prevent the appearance of juvenile kelp sporophytes for three months and the establishment of a diverse kelp assemblage for over four months within a southern California kelp forest. Juveniles were observed within one month in plots where microscopic stages were left intact, which may confer an advantage for the resulting sporophytes as they attain larger sizes before later recruiting neighbors. Microsatellite diversity was high (expected heterozygosity HE approximately 0.9) for juveniles and adults within our sites. Using a microsatellite-based parentage analysis for the dominant kelp, Macrocystis pyrifera, we estimated that a portion of the new M. pyrifera sporophyte recruits had originated from their parents at least seven months after their parents had disappeared. Similar delay durations have been demonstrated in recent laboratory studies. Additionally, our results suggest that zoospore dispersal distances > 50 m may be supported by including additional microsatellite loci in the analysis. We propose a mixed-age and, potentially, a mixed-origin bank of M. pyrifera gametophytes promotes maximal genetic diversity in recovering populations and reduces population genetic subdivision and self-fertilization rates for intact populations by promoting the survival of zoospores dispersed > 10 m and during inhospitable environmental conditions.

Looking into the black box: simulating the role of self-fertilization and mortality in the genetic structure of Macrocystis pyrifera.

Patterns of spatial genetic structure (SGS), typically estimated by genotyping adults, integrate migration over multiple generations and measure the effective gene flow of populations. SGS results can be compared with direct ecological studies of dispersal or mating system to gain additional insights. When mismatches occur, simulations can be used to illuminate the causes of these mismatches. Here, we report a SGS and simulation-based study of self-fertilization in Macrocystis pyrifera, the giant kelp. We found that SGS is weaker than expected in M. pyrifera and used computer simulations to identify selfing and early mortality rates for which the individual heterozygosity distribution fits that of the observed data. Only one (of three) population showed both elevated kinship in the smallest distance class and a significant negative slope between kinship and geographical distance. All simulations had poor fit to the observed data unless mortality due to inbreeding depression was imposed. This mortality could only be imposed for selfing, as these were the only simulations to show an excess of homozygous individuals relative to the observed data. Thus, the expected data consistently achieved nonsignificant differences from the observed data only under models of selfing with mortality, with best fits between 32% and 42% selfing. Inbreeding depression ranged from 0.70 to 0.73. The results suggest that density-dependent mortality of early life stages is a significant force in structuring Macrocystis populations, with few highly homozygous individuals surviving. The success of these results should help to validate simulation approaches even in data-poor systems, as a means to estimate otherwise difficult-to-measure life cycle parameters.

Transcriptomic analysis of metabolic function in the giant kelp, Macrocystis pyrifera, across depth and season.

To increase knowledge of transcript diversity for the giant kelp, Macrocystis pyrifera, and assess gene expression across naturally occurring depth gradients in light, temperature and nutrients, we sequenced four cDNA libraries created from blades collected at the sea surface and at 18 m depth during the winter and summer. Comparative genomics cluster analyses revealed novel gene families (clusters) in existing brown alga expressed sequence tag data compared with other related algal groups, a pattern also seen with the addition of M. pyrifera sequences. Assembly of 228 Mbp of sequence generated c. 9000 isotigs and c. 12,000 open reading frames. Annotations were assigned using families of hidden Markov models for c. 11% of open reading frames; M. pyrifera had highest similarity to other members of the Phaeophyceae, namely Ectocarpus siliculosus and Laminaria digitata. Quantitative polymerase chain reaction of transcript targets verified depth-related differences in gene expression; stress response and light-harvesting transcripts, especially members of the LI818 (also known as LHCSR) family, showed high expression in the surface compared with 18 m depth, while some nitrogen acquisition transcripts (e.g. nitrite reductase) were upregulated at depth compared with the surface, supporting a conceptual biological model of depth-dependent physiology.

Isolation by oceanographic distance explains genetic structure for Macrocystis pyrifera in the Santa Barbara Channel.

Ocean currents are expected to be the predominant environmental factor influencing the dispersal of planktonic larvae or spores; yet, their characterization as predictors of marine connectivity has been hindered by a lack of understanding of how best to use oceanographic data. We used a high-resolution oceanographic model output and Lagrangian particle simulations to derive oceanographic distances (hereafter called transport times) between sites studied for Macrocystis pyrifera genetic differentiation. We build upon the classical isolation-by-distance regression model by asking how much additional variability in genetic differentiation is explained when adding transport time as predictor. We explored the extent to which gene flow is dependent upon seasonal changes in ocean circulation. Because oceanographic transport between two sites is inherently asymmetric, we also compare the explanatory power of models using the minimum or the mean transport times. Finally, we compare the direction of connectivity as estimated by the oceanographic model and genetic assignment tests. We show that the minimum transport time had higher explanatory power than the mean transport time, revealing the importance of considering asymmetry in ocean currents when modelling gene flow. Genetic assignment tests were much less effective in determining asymmetry in gene flow. Summer-derived transport times, in particular for the month of June, which had the strongest current speed, greatest asymmetry and highest spore production, resulted in the best-fit model explaining twice the variability in genetic differentiation relative to models that use geographic distance or habitat continuity. The best overall model also included habitat continuity and explained 65% of the variation in genetic differentiation among sites.

Habitat continuity and geographic distance predict population genetic differentiation in giant kelp.

Isolation by distance (IBD) models are widely used to predict levels of genetic connectivity as a function of Euclidean distance, and although recent studies have used GIS-landscape ecological approaches to improve the predictability of spatial genetic structure, few if any have addressed the effect of habitat continuity on gene flow. Landscape effects on genetic connectivity are even less understood in marine populations, where habitat mapping is particularly challenging. In this study, we model spatial genetic structure of a habitat-structuring species, the giant kelp Macrocystis pyrifera, using highly variable microsatellite markers. GIS mapping was used to characterize habitat continuity and distance between sampling sites along the mainland coast of the Santa Barbara Channel, and their roles as predictors of genetic differentiation were evaluated. Mean dispersal distance (sigma) and effective population size (Ne) were estimated by comparing our IBD slope with those from simulations incorporating habitat continuity and spore dispersal characteristics of the study area. We found an allelic richness of 7-50 alleles/locus, which to our knowledge is the highest reported for macroalgae. The best regression model relating genetic distance to habitat variables included both geographic distance and habitat continuity, which were respectively, positively and negatively related to genetic distance. Our results provide strong support for a dependence of gene flow on both distance and habitat continuity and elucidate the combination of Ne and a that explained genetic differentiation.