Sex-dependent and -independent transcriptional changes during haploid phase gametogenesis in the sugar kelp Saccharina latissima


Authors: Gareth A. Pearson aff001;  Neusa Martins aff001;  Pedro Madeira aff001;  Ester A. Serrão aff001;  Inka Bartsch aff002
Authors place of work: Centre for Marine Sciences (CCMAR)-CIMAR, University of Algarve, Portugal aff001;  Alfred-Wegener-Institute, Helmholtz Center for Polar and Marine Research, Am Handelshafen, Germany aff002
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: 10.1371/journal.pone.0219723

Summary

In haplodiplontic lineages, sexual reproduction occurs in haploid parents without meiosis. Although widespread in multicellular lineages such as brown algae (Phaeophyceae), haplodiplontic gametogenesis has been little studied at the molecular level. We addressed this by generating an annotated reference transcriptome for the gametophytic phase of the sugar kelp, Saccharina latissima. Transcriptional profiles of microscopic male and female gametophytes were analysed at four time points during the transition from vegetative growth to gametogenesis. Gametogenic signals resulting from a switch in culture irradiance from red to white light activated a core set of genes in a sex-independent manner, involving rapid activation of ribosome biogenesis, transcription and translation related pathways, with several acting at the post-transcriptional or post-translational level. Additional genes regulating nutrient acquisition and key carbohydrate-energy pathways were also identified. Candidate sex-biased genes under gametogenic conditions had potentially key roles in controlling female- and male-specific gametogenesis. Among these were several sex-biased or -specific E3 ubiquitin-protein ligases that may have important regulatory roles. Females specifically expressed several genes that coordinate gene expression and/or protein degradation, and the synthesis of inositol-containing compounds. Other female-biased genes supported parallels with oogenesis in divergent multicellular lineages, in particular reactive oxygen signalling via an NADPH-oxidase. Males specifically expressed the hypothesised brown algal sex-determining factor. Male-biased expression mainly involved upregulation of genes that control mitotic cell proliferation and spermatogenesis in other systems, as well as multiple flagella-related genes. Our data and results enhance genome-level understanding of gametogenesis in this ecologically and economically important multicellular lineage.

Keywords:

Biology and life sciences – Physiology – Reproductive physiology – Gametogenesis – Biochemistry – Proteins – DNA-binding proteins – Motor proteins – Protein domains – Genetics – Gene expression – Genomics – Genome analysis – Transcriptome analysis – Cell biology – Molecular motors – Motor proteins – Organisms – Eukaryota – Plants – Algae – Seaweed – Computational biology – Genome analysis – Transcriptome analysis – Medicine and health sciences – Physiology – Reproductive physiology – Gametogenesis

Introduction

Gametogenesis is a fundamental process for sexual reproduction. In diploid organisms, the sequential activation of two major developmental programs are coordinated; expression changes causing a switch from mitotic cell division (i.e., growth) to reductive meiotic division producing haploid cells, and subsequent activation of the morphogenetic program to produce differentiated haploid (male and female) gametes [1]. However, haplodiplontic lineages deviate from this pattern, because sex is expressed in the haploid phase of the life cycle, thereby uncoupling gametogenesis from meiosis (see [2,3]). Haploid phase gametogenesis can thereby be viewed as a relatively simplified process, limited to the integration of gametogenic signals to fix cell-fate decisions and produce gametes.

Brown algae in the family Laminariales (kelps, sensu stricto) maintain an extremely heteromorphic haplodiplontic life cycle alternating between haploid microscopic gametophytes and diploid macroscopic sporophytes. Kelp sporophytes are the most developmentally complex and largest members of the Phaeophyceae, an independently evolved eukaryotic multicellular lineage [4]. In contrast, gametophytes are cryptic, either free-living or endophytic [5] with a highly simplified filamentous morphology specialised for gamete production. In addition to their unique evolutionary position on the tree of multicellular life, kelp gametophytes are an attractive model system for the study of gametogenesis. Gametophytes can be isolated and maintained in long-term vegetative culture (i.e., for decades; see [6,7]). Gametogenesis can be induced in the laboratory by blue light [8] or Fe addition [9], and development to fertilization is rapid under optimal conditions (ca. 10 days). Sexual development proceeds independently in male and female gametophytes (dioicy), with external fertilization by biflagellate sperm that are discharged and attracted in response to the release by eggs of the pheromone lamoxirene [10].

Although evidence for genetic sex determination has long existed [11], the UV sex determination system in Ectocarpus [12] and related brown algae including Laminariales [13] was confirmed only recently, as genomic resources for brown algae have become available [4, 14]. Thus, haploid sex chromosome evolution across a broad evolutionary range of brown algae has recently advanced considerably [13], reviewed by [15], and some studies have looked at gene expression underlying sexual dimorphism [16, 17]. However, information about gene expression changes during gametogenesis in haplodiplontic and oogamous kelps remains limited (but see [18]).

There are both ecological and economic rationales for studying the genetic and genomic underpinnings of kelp sexual development. The fundamental role of kelps in structuring marine forest ecosystems across cold to temperate regions of both hemispheres is increasingly threatened by climate change-induced range shifts, over-exploitation, and habitat destruction [1923]. The sugar kelp Saccharina latissima is a highly successful amphi-boreal species [24], in which a diverse habitat range has promoted the development of distinct morpho- and ecotypes [24, 25] that may indicate a large reservoir of functional genetic variation with potential use in commercial breeding programs [26]. It also suggests that as range shifts gather pace, elements of this diversity could be lost.

To address this knowledge gap, we have undertaken the first transcriptional study of male and female gametogenesis in S. latissima. Our results provide evidence both for common elements expressed during the commitment to gametogenesis from the vegetative phase, as well as sex-specific developmental programs in males and females.

Materials and methods

Experimental culture conditions and sampling

Unialgal female and male gametophyte cultures of Saccharina latissima, each derived from the meiospores of a single sporophyte, were isolated from Oslo, Norway (AWI seaweed culture collection: ♀ 3301, ♂ 3300) and Spitsbergen, Svalbard (AWI seaweed culture collection: ♀ 3124, ♂ 3123). They were maintained in a vegetative stage at 15°C in 3 μmol photons m-2 s-1 of red light (RL; LED Mitras daylight 150 controlled by ProfiLux 3, GHL Advanced Technology, Kaiserslautern, Germany) under a 16:8 h light:dark cycle in sterile full strength Provasoli enriched seawater (PES; [27]). Gametophyte vegetative growth was enhanced to have sufficient material for the gametogenesis induction experiment and RNA extraction; female and male gametophyte material from both strains was transferred to seawater enriched with 8×10−4 M N and 1×10−4 M P, at 17°C, under continuous irradiance of 6 μmol photons m-2 s-1 of RL for 6 weeks. The seawater medium was renewed weekly. The irradiance conditions chosen for vegetative gametophyte growth were optimal based on initial experiments indicating improved culture health (mortality, qualitative assessment of pigmentation) in RL compared to WL (pers. obs.). Nutrient conditions were adapted from [28].

Induction of gametophyte reproduction

To induce gametogenesis, female and male gametophytes from both strains (biological replicates) were gently separated into small fragments (~1 mm in length) and cultured separately into Petri dishes (9.5 cm diameter, height 5.5 cm) with ½ strength PES at 10 ºC. The irradiance was 15 μmol photons m-2 s-1 of white light (WL; LED MITRAS lightbars 150 Daylight, GHL Advanced Technology, Kaiserslautern, Germany) in a 16:8 h light:dark cycle. Male and female gametophytes were sampled for RNA (flash-frozen in liquid nitrogen and stored at -80 ºC) after 1, 6 and 8 days in WL to follow the developmental stages of gametogenesis. Control gametophytes were sampled under RL growth conditions described above (0 days vegetative control). All samples were taken before gametophytes released eggs and sperm, based on observations of the cultures (Fig 1).

<h2>Cultured <i>Saccharina latissima</i> gametophytes used in this study.</h2>
Fig. 1.

Cultured Saccharina latissima gametophytes used in this study.


Vegetative filaments of A) males and B) females under RL vegetative growth conditions. The same male C) and female D) cultures after 8 d in WL gametogenic conditions. Solid and open arrows indicate sites of antheridial and oogonial development, respectively. Scale bar = 50 μm.

RNAseq analysis and de novo reference transcriptome assembly

Total RNA was extracted from lyophilized tissue equivalent to between 100–200 mg FW male and female gametophyte culture per sample following established protocols [29] and sequenced by a service provider (100 bp paired-end Illimina HiSeq 4000; BGI, China).

The raw sequence data were evaluated with standard quality control tools (FastQC v0.11.7; https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Prior to de novo transcriptome assembly, quality-filtered male and female reads were digitally normalized to reduce redundancy, equalize k-mer coverage and remove rare k-mers potentially arising from sequencing errors [30]. Assemblies (Velvet-Oases; [31]) were built over the k-mer range 21–61 (step size of 10), and these 5 assemblies were subsequently merged with transfuse (https://github.com/cboursnell/transfuse). The merged transcriptome was queried by Diamond [32] in BLASTX mode against Stramenopile proteins (subset from NCBI nr); contigs with top blast hits against Phaeophyceae were retained, and the remaining contigs containing potential contaminants were removed from the analysis.

A proportion of contigs from brown algae de novo transcriptomes may be “polycistronic”, likely resulting from brown algal genome structure, in which closely adjacent genes are transcribed from alternate strands [4]. Therefore, to produce the final reference transcriptome, putative open reading frames (ORFs) were identified with FragGeneScan [33], and clustered at 97% nucleotide identity with VSEARCH [34]. Sequences with length < 200 nt were discarded.

Transcriptome completeness

We assessed the completeness of the final merged, Phaeophyceae-screened reference transcriptome with BUSCO v2.0 [35], by querying predicted proteins against the eukaryotic reference database, using online resources provided by CyVerse https://www.cyverse.org/.

Differential expression analysis

High quality reads were mapped onto the reference transcriptome using the RSEM (v1.2.31) wrapper script and Bowtie2 [36, 37]. Expected count data were analysed in Bioconductor 3.8 using edgeR and limma [38, 39]. A total of 11,916 transcripts were retained for analysis after filtering for transcripts with > 4 counts per million (CPM, approximately 20 reads/transcript) in at least 6 samples. Samples from the two available strains of S. latissima (SLO—Oslofjord and SLS–Spitzbergen) were used as biological replicates to investigate transcriptome expression profiles in response to the factors “sex” (two levels; male and female [M and F]) and “time” (four levels; vegetative growth in RL [= day 0], and 1, 6 and 8 days following a transfer to WL to initiate gametogenesis). Differential expression (DE) was used to identify up- and down-regulated transcripts between groups defined by combinations of sex and time (false discovery rate [FDR] ≤ 0.05). The analysis was repeated on KEGG orthology (KO) terms, after summing transcript reads corresponding to unique KO entries.

Differentially expressed genes reported in this study were subject to additional confirmatory phylogenetic checks, since known pathogens and parasites of brown algae are themselves stramenopiles (oomycetes, labyrinthulids), while other phototrophic stramenopile lineages such as diatoms may be present in small numbers in culture. To do this, amino acid alignments were made (Muscle; [40]) from all blastx hits (Expect-value cutoff ≤ 1e-10), alignments were trimmed with GBlocks [41] before phylogenetic tree construction (PhyML, LG model and aLRT branch support; [42]). Some suspect contigs sister to labyrinthulids or in very few cases diatoms were screened and removed from further analysis.

Heatmaps (ComplexHeatmap R package; [43]) were prepared from a matrix of all DE transcripts or KO terms for all samples, without applying a fold change cutoff. The clustering distance and clustering methods were “spearman” and “ward.D2”, respectively. Clusters were identified with NbClust in R, focusing on the Gap statistic method with 500 bootstrap replicates. Clustering was used as a guide to identify major patterns and gene / annotation lists from the data. Venn diagrams and MA plots (log ratio versus mean expression) were generated in R using the limma and ggplot2 [44] packages.

Functional annotation and gene set enrichment

Transcripts were functionally annotated by Diamond (BLASTX mode) comparisons against Ectocarpus Ec32 strain proteins. Gene ontology (GO) terms were appended to transcripts based on annotation data for the Ectocarpus genome (https://bioinformatics.psb.ugent.be/gdb/ectocarpusV2/). Predicted protein sequences derived from S. latissima reference transcripts were annotated against the InterPro [45] and KEGG [46] databases. Gene set enrichment analysis (GSEA) of GO terms was performed using clusterProfiler [47] in Bioconductor 3.8.

Transcripts with homology to genes from the male and female sex-determining regions (sdr) of Ectocarpus and their expression level in S. latissima were collated from functional annotation (Diamond BLASTX) and mapping (RSEM) data. Uniquely-expressed male and female transcripts were defined as those with average mapping ≥ 4 transcripts per million (TPM) in the expressed sex and ≤ 1 TPM in the non-expressed sex, identified from RSEM expression results.

Results and discussion

Transcriptome sequencing and expression profiles

High throughput sequencing resulted in 52.5 Gb of high quality paired-end read data (266.2 million male and 259.2 million female), with an average of 3.28 Gb (32.8 ± 1.04 million reads) per sample (S1 Table). Digital normalization reduced the dataset to 50.7 M reads for assembly with Velvet-Oases. After merging individual k-mer assemblies with transfuse 323,917 contigs remained with an N50 of 1,509 bp. After extracting and screening predicted open reading frames for Phaeophyceae-specific top hits, the final reference transcriptome consisted of 34,002 contigs with an N50 of 1,281 bp (S1 Table). Transcriptome completeness indicated that 64.1% of Ectocarpus protein coding genes were represented. According to BUSCO v2.0 analysis, the transcriptome was 74.9% complete (with 25.9% duplication level; S1 Table) and is therefore reasonably complete for a de novo transcriptome lacking representation from the sporophyte stage. This value rose to > 90% when fragmented or partial matches were considered. For comparison, running the same analysis with Ectocarpus strain Ec32 V2 proteins resulted in an estimate of 82.5% completeness and 9.1% duplication.

Differential expression (DE; over- or under-representation between sexes or between timepoints within a sex) was detected for 1,122 of the 11,916 transcripts analysed in edgeR (9.4%). A total of 3,429 KEGG orthologues were identified in the dataset, of which 521 (15.2%) were differentially expressed (a full list of KEGG-annotated DE genes is provided in S2 Table). Heatmap clustering of DE KO terms revealed several distinct patterns of gene regulation (Fig 2). The clustering indicated a primary divergence in expression between RL and WL, overriding differences due to gender. Further clustering of samples was primarily sex-dependent both during vegetative growth (RL) and gametogenesis. Transcript clusters revealed the presence of both female and male “constitutive” genes (Fig 2, clusters 1 and 8), with others representing transcripts predominantly up- or down-regulated during WL exposure (clusters 6 and 2, respectively), or showing transient or sex-dependent regulation (clusters 3–5,7).

<h2>Cluster analysis of differential gene expression.</h2>
Fig. 2.

Cluster analysis of differential gene expression.


Heatmap of 512 KEGG-annotated genes showing differential expression between gametophyte transcriptomes, either between sexes (F = female; M = male gametophytes) or timepoints (time = 0, 1, 6, 8 d). Expression values for each KEGG gene (row) are normalized across all samples (columns) by Z-score. Both column and row clustering were applied, and distinct gene clusters identified by the Gap statistic method are shown to illustrate the major expression patterns observed in the data.

Transcriptional profiles of female gametophytes are shown in Fig 3, as comparisons between vegetative (RL) growth conditions and the subsequent timecourse of gametogenesis after 1–8 days in WL. The majority of up-regulated KEGG-annotated genes appeared early (by 1 day) following exposure to WL (105 of the total 125). Of this total, 67 KEGG genes were uniquely up-regulated after 1 day, with only 20 genes additionally up-regulated after 6 and/or 8 days (Fig 3, upper Venn diagram). In male gametophytes, “early responsive” KEGG genes also dominated (Fig 4; 117 of 217), but the contribution of “late-responsive” genes was considerably higher than in females (100 additional KEGG genes after 6 and/or 8 days).

<h2>Female gametophyte expression.</h2>
Fig. 3.

Female gametophyte expression.


MA plots (log expression ratio vs. mean average expression) comparing female gametophyte gene expression (as KEGG-annotated genes) under vegetative growth (d0) with gametogenic conditions after a) 1 d, b) 6 d, and c) 8 d in WL. Each point on the plots represents a unique KEGG gene, with differentially expressed genes shown as larger orange points (edgeR, FDR < 0.05). The plots show expression on the x-axes as average Log2 counts per million (CPM), and the ratio of RL/WL expression as Log2(RL/WL) is shown on the y-axes. Venn diagrams summarise the expression changes across timepoints for genes up-regulated (upper Venn) and down-regulated (lower Venn) in gametogenic (WL) compared with vegetative (RL) conditions. A full list of KEGG gene annotations can be found in S2 Table.

<h2>Male gametophyte expression.</h2>
Fig. 4.

Male gametophyte expression.


MA plots (log expression ratio vs. mean average expression) comparing male gametophyte gene expression (as KEGG-annotated genes) under vegetative growth (d0) with gametogenic conditions after a) 1 d, b) 6 d, and c) 8 d in WL. Each point on the plots represents a unique KEGG gene, with differentially expressed genes shown as larger orange points (edgeR, FDR < 0.05). The plots show expression on the x-axes as average Log2 counts per million (CPM), and the ratio of RL/WL expression as Log2(RL/WL) is shown on the y-axes. Venn diagrams summarise the expression changes across timepoints for genes up-regulated (upper Venn) and down-regulated (lower Venn) in gametogenic (WL) compared with vegetative (RL) conditions. A full list of KEGG gene annotations can be found in S2 Table.

Functional annotation

Gene set enrichment analysis (GSEA) of GO terms comparing the sexes showed that females were enriched for “ribosome”, “translation” and related terms throughout the timecourse, suggesting that transcripts encoding ribosomal proteins are relatively more abundant in females during both vegetative growth and gametogenesis (S3 Table). After 1 and 6 days in WL, translational activity in females was also enriched (“translation initiation factor activity”, “eukaryotic translation initiation factor 3 complex”; S3 Table). We also found GO:0055114 “oxidation-reduction process” overrepresented in female gametophytes (on day 6). This gene set contains several uniquely or highly over-expressed genes in females, important among which seem to be group of genes related to wounding/pathogenesis. These include transcripts for a lipoxygenase, a respiratory burst oxidase, and a manganese SOD, all of which were female specific (with no detectable expression in male gametophytes).

Under vegetative growth conditions, male gametophytes showed enrichment of several terms, including “galactosylceramide sulfotransferase activity” and “golgi”, containing several transcripts annotated as Galactose-3-O-sulfotransferases. This latter group of transcripts included 2 members uniquely expressed in males throughout gametogenesis (see “Gender-biased or -specific expression in males” below). Male vegetative gametophytes were also enriched in “ubiquitin-protein transferase activity” and “regulation of transcription, DNA-templated”. For the latter, a range of transcription factors were represented, including heat shock factors (HSF), TFIIB, TFIIH, CCR4-Not (nuclear transcription from RNA polymerase II), sigma-70 (plastid transcription), as well as genes involved in chromatin remodelling (e.g., TAZ-type Zn finger protein, histone acetyltransferase). TALE-like homeodomain transcription factors were also represented, with homology to the Ectocarpus ORO life-cycle regulator [48, 49] as well as Ec-04_000450. Homologues for SAM, a second TALE transcription factor key to the regulation of life-cycle phase in Ectocarpus, which heterodimerizes with ORO [49] were not detected in our dataset. Following transfer to WL to induce gametogenesis, no overrepresented GO terms were found in males (S3 Table).

Gametogenesis involves regulation of key gender-independent “early response” genes

A majority of “early responsive” genes (i.e., significant expression changes after 1 day in WL) were common to both sexes, with 69 KEGG genes (or approximately two thirds) up-regulated in both males and females (S2 Table). These genes provide insights into general cellular, developmental and metabolic processes triggered during early gametogenesis.

From the expression changes in response to WL it can be inferred that the gametogenesis developmental program involves rapid activation of ribosome, transcription and translation related pathways. Furthermore, it was notable that several genes were involved in post-transcriptional or epigenetic regulation. An interesting example is SETD6 (K05302), a member of the SET domain family of protein lysine methyltransferases (PKMTs). The only differentially expressed member of six SET domain proteins identified in our dataset, SETD6 was immediately up-regulated under WL in both male and female gametophytes (Table 1). Reported targets of SETD6 methylation suggest a key role in gametogenesis. By methylating a histone H2A variant, SETD6 is directly implicated in the control of cellular differentiation in mouse embryonic stem cells [50]. SETD6 also functions in cell cycle regulation, via nonhistone methylation of the positive regulator of mitosis PLK1 (polo-like kinase 1), thereby exerting control over the rate of cell division [51].

Tab. 1.

Core “early responsive” KEGG genes in female and male gametophytes.

<h2>Core “early responsive” KEGG genes in female and male gametophytes.</h2>

In total, 11 KEGG genes belonging to the eukaryotic ribosome biogenesis pathway were up-regulated, including genes involved in rRNA 2´-O-methyation (NOP56, NOP58) and pseudouridylation (DKC1, NHP2) in the nucleolus, the co-translational acetylation of proteins (NAT10), the nuclear chaperone of 60S rRNA (midasin), as well as members of the UTP-B and MPP10 complexes associated with 90S pre-ribosomal RNA. A gene essential for 16S rRNA maturation and assembly of 30S ribosome subunits (ERAL1), as well as several genes involved in post-transcriptional modification (pseudouridylation) of tRNA and 23S rRNA were also up-regulated (DUS1, PUS7, rluB). Among other early responsive genes were several DEAD box RNA helicases, required for mitochondrial intron splicing (MSS116), and nuclear rRNA synthesis and processing (DBP3, DDX21), a Pumilio domain homology protein implicated in translational control of mRNA, as well as several translation initiation factors (Table 1 and S2 Table). Ribosome biogenesis and protein synthesis gene networks play key roles in controlling germline stem cell differentiation in animals [52], and mutations affecting rRNA processing also affect plant gametogenesis and embryogenesis [53, 54]. Although our data are limited to observations of transcript levels, the evidence for involvement of similar pathways across independent multicellular lineages appears compelling. The central role for transcriptional and translational control in gametogenesis suggests that proteomic comparisons would complement RNAseq-based approaches. Another promising future direction might be mutant screens in the model brown algal system (such as Ectocarpus) to analyse phenotypes impaired in gamete formation.

The rapid regulation/reorganization of the transcriptional and translational machinery was accompanied by metabolic adjustment. Notably, several key genes in nitrogen metabolism were upregulated (Table 1). These included strong up-regulation (ca. 20-fold) of a nitrate/nitrite transporter (NRT), suggesting an increased requirement for inorganic N acquisition for gametogenesis, reminiscent of the response of diatoms under N-deprivation [55]. Although N and P levels in the medium were lower after transfer to WL gametogenic conditions, and we cannot exclude an effect on uptake rates, the medium was still nutrient replete. N-acquisition was accompanied by upregulation of ammonium assimilation genes contributing to cellular pools of L-glutamate including both chloroplastic Fd-GOGAT (GLT1) and cytosolic NADPH-GOGAT (gltD), and possibly GDH (but see [56, 57]). Biosynthesis of several amino-acids depends on available glutamate (e.g., alanine, aspartate and glutamate metabolism, arginine metabolism) and these pathways also showed evidence of upregulation during gametogenesis.

Pyruvate is both the product of glycolysis and starting metabolite for gluconeogenesis, and acts as a key metabolite in interacting energy and carbon metabolism pathways. In the penultimate and final steps of glycolysis to generate pyruvate, enolase (ENO) and pyruvate kinase (PK) were upregulated, as was pyruvate, phosphate dikinase (PPDK) which reverses the action of PK to generate phosphoenolpyruvate (S2 Table). Pyruvate and glutamate are also produced by transamination of alanine by alanine transaminase (ALT), which was strongly upregulated (S2 Table). Pyruvate is the precursor for acetyl CoA entering the TCA (citric acid cycle), which was moderately upregulated together with pyruvate and diydrolipoamide dehydrogenases (PDHA, DLAT and DLD).

We observed no clear pathway-level responses among down-regulated “early responsive” genes, although we noted a transitory down-regulation of light-harvesting (LHCA1 and 4) and photosynthesis proteins (psbQ, psbU). This was, however, followed by subsequent up-regulation at days 6–8, and can likely be attributed to increased excitation pressure on PS2 caused by lower temperature and higher light after the transition from RL to WL [58].

Gender-biased or -specific expression in females

The identification of “core” genes constitutively over- (or uniquely) expressed in the two genders can provide insight into sex-specific genetic networks and developmental programs. In female gametophytes, 15 KEGG genes were statistically assigned to the core set (Fig 5A, 5C–5F). These included several genes with no detectable expression in male gametophytes that thus may provide insight into the control of female gametogenesis (Table 2). They include a transcript with repeated RCC1 (regulator of chromosome condensation) motifs and annotated as KEGG orthologue HERC3 (K10614, an E3 ubiquitin-protein ligase). We also identified a female-specific Ran GTPase (SL_90811), for which RCC1 is the guanine nucleotide exchange factor. Together they are thought to play an important role in nucleocytoplasmic transport and mitotic regulation [59]. Another potential HERC3 interacting protein is KCTD9 (K21919), which in humans interacts with the E3 ubiquitin ligase complex to mediate the ubiquination of target proteins for degradation, and which was expressed (at low levels) in WL only in female gametophytes. The pre-mRNA splicing factor SYF1 (an LSM3 isoform and component of the spliceosome) was also uniquely expressed in female gametophytes (Fig 5A, 5C–5F). Taken together, these results suggest that several factors are involved in the control of female-specific gene expression and the female gametogenic developmental program in S. latissima.

<h2>Comparative expression in males and females during gametogenesis.</h2>
Fig. 5.

Comparative expression in males and females during gametogenesis.


Venn diagrams showing numbers of differentially expressed genes (KEGG-annotated genes) upregulated at each experimental timepoint in females a) and males b). Differentially expressed genes under vegetative conditions only are highlighted in red, “core” genes overexpressed in all conditions in one of the sexes are highlighted in bold, while WL-responsive genes are highlighted in purple (females) or blue (males). MA plots (log expression ratio vs. mean average expression) of female vs. male gametophyte expression of KEGG-annotated genes in vegetative growth conditions c) 0 days, and after d) 1 day, e) 6 days and f) 8 days of culture under WL gametogenic conditions. Each point represents a unique KEGG gene, with differentially expressed shown as larger green points (edgeR; FDR < 0.05) and gene labels are colour-coded according to the Venn diagrams to the left (vegetative = red; “core” = black; WL-responsive = purple or blue for females and males, respectively). The plots show expression on the x-axes as average Log2 counts per million (CPM), and the ratio of female/male expression as Log2(F/M) is shown on the y-axes. KEGG gene information: 1) Female vegetative (RL): lplA: lipoate-protein ligase; ATG101: autophagy-related protein 101; ATP23: mitochondrial inner membrane protease; NOP15: nucleolar protein 15; TOM40: mitochondrial import receptor subunit; SEC13: protein transport protein; fabI: enoyl-[acyl-carrier protein] reductase I; RRP41: exosome complex component; PSMB5: 20S proteasome subunit beta 5; EIF2S1: translation initiation factor 2 subunit 1. 2) Female “core” up-regulated: HERC3: E3 ubiquitin-protein ligase; phoB1: two-component system, OmpR family, alkaline phosphatase synthesis response regulator; pepM: phosphoenolpyruvate phosphomutase; CTSE: cathepsin E; NOX2: NADPH oxidase 2; NAGLT1: MFS transporter, FHS family, Na+ dependent glucose transporter 1; SLC24A5: solute carrier family 24 (sodium/potassium/calcium exchanger), member 5; SYF1: pre-mRNA-splicing factor; ECT2: protein ECT2; KCTD9: BTB/POZ domain-containing protein; LSM3: U6 snRNA-associated Sm-like protein; INO1: myo-inositol-1-phosphate synthase; ispE: 4-diphosphocytidyl-2-C-methyl-D-erythritol kinase; MSH2: DNA mismatch repair protein; HPGDS: prostaglandin-H2 D-isomerase / glutathione transferase. 3) Male “core” up-regulated: HMGB2: high mobility group protein B2; SHPRH: E3 ubiquitin-protein ligase; GAL3ST3: galactose-3-O-sulfotransferase 3. 4) Female WL-responsive: ABHD11: abhydrolase domain-containing protein 11; ABHD17: abhydrolase domain-containing protein 17; ACACA: acetyl-CoA carboxylase / biotin carboxylase 1; accC: acetyl-CoA carboxylase, biotin carboxylase subunit; AGPHD1: hydroxylysine kinase; APTX: aprataxin; ASPM: abnormal spindle-like microcephaly-associated protein; CEP41: centrosomal protein; CHMP4: charged multivesicular body protein 4; D2HGDH: D-2-hydroxyglutarate dehydrogenase; DECR2: peroxisomal 2,4-dienoyl-CoA reductase; DEPDC5: DEP domain-containing protein 5; DPM1: dolichol-phosphate mannosyltransferase; dxr: 1-deoxy-D-xylulose-5-phosphate reductoisomerase; EIF2B2: translation initiation factor eIF-2B subunit beta; gdhA: glutamate dehydrogenase (NADP+); groES: chaperonin GroES; GSTK1: glutathione S-transferase kappa 1; HDAC11: histone deacetylase 11; ilvH: acetolactate synthase I/III small subunit; ITGB1: integrin beta 1; LYRM4: LYR motif-containing protein 4; MCM7: DNA replication licensing factor; MINDY3_4: ubiquitin carboxyl-terminal hydrolase; msrB: peptide-methionine (R)-S-oxide reductase; mtnD: 1,2-dihydroxy-3-keto-5-methylthiopentene dioxygenase; NDUFAB1: NADH dehydrogenase (ubiquinone) 1 alpha/beta subcomplex 1; NIT1: deaminated glutathione amidase; NTO1: NuA3 HAT complex component; nuoB: NADH-quinone oxidoreductase subunit B; nusG: transcriptional antiterminator; PDCD5: programmed cell death protein 5; pepD: dipeptidase D; petF: ferredoxin; PEX4: peroxin-4; PGK: phosphoglycerate kinase; POLB: DNA polymerase beta; PPIH: peptidyl-prolyl isomerase H (cyclophilin H); RP-L15: large subunit ribosomal protein L15; SLC25A4S: solute carrier family 25 (mitochondrial adenine nucleotide translocator); thiG: thiazole synthase; trxA: thioredoxin 1; USP34: ubiquitin carboxyl-terminal hydrolase 34; XDH: xanthine dehydrogenase/oxidase; XRCC2: DNA-repair protein; yggS: PLP dependent protein; yhbH: putative sigma-54 modulation protein. 5) Male WL-responsive: ALG5: dolichyl-phosphate beta-glucosyltransferase; ANO7: anoctamin-7; APC7: anaphase-promoting complex subunit 7; argC: N-acetyl-gamma-glutamyl-phosphate reductase; ARL2BP: ADP-ribosylation factor-like protein 2-binding protein; ARL3: ADP-ribosylation factor-like protein 3; ARL8: ADP-ribosylation factor-like protein 8; AURKX: aurora kinase; AVIL: advillin; BBS9: Bardet-Biedl syndrome 9 protein; CCNT: cyclin T; CCT6: T-complex protein 1 subunit zeta; CETN3: centrin-3; CK: creatine kinase; CNOT1: CCR4-NOT transcription complex subunit 1; DEK: protein DEK; DNAAF2: dynein assembly factor 2, axonemal; DNAI1: dynein intermediate chain 1, axonemal; DNAJB1: DnaJ homolog subfamily B member 1; DNAL1: dynein light chain 1, axonemal; DNALI: dynein light intermediate chain, axonemal; DUR3: urea-proton symporter; DYNC2H: dynein heavy chain 2, cytosolic; DYNC2LI: dynein light intermediate chain 2, cytosolic; DYX1C1: dyslexia susceptibility 1 candidate gene 1 protein; ELF2C: eukaryotic translation initiation factor 2C; ERN1: serine/threonine-protein kinase/endoribonuclease IRE1; GATM: glycine amidinotransferase; GINS4: GINS complex subunit 4; HECTD3: E3 ubiquitin-protein ligase; HNRNPR: heterogeneous nuclear ribonucleoprotein R; HS3ST5: [heparan sulfate]-glucosamine 3-sulfotransferase 5; IFT22: intraflagellar transport protein 22; IFT43: intraflagellar transport protein 43; IFT46: intraflagellar transport protein 46; IFT80: intraflagellar transport protein 80; IFT81: intraflagellar transport protein 81; IFT172: intraflagellar transport protein 172; ILK: integrin-linked kinase; ITPR1: inositol 1,4,5-triphosphate receptor type 1; KIF4_21_27: kinesin family member; KIF13: kinesin family member 13; lysC: aspartate kinase; MBD4: methyl-CpG-binding domain protein 4; MCM3: DNA replication licensing factor; MYO7A: myosin VIIa; ORC1: origin recognition complex subunit 1; PELI: pellino; PIN1: peptidyl-prolyl cis-trans isomerase NIMA-interacting 1; POLA1: DNA polymerase alpha subunit A; PPP1R42: protein phosphatase 1 regulatory subunit 42; PRKX: protein kinase X; PRPF19: pre-mRNA-processing factor 19; RAB28: Ras-related protein; RDH12: retinol dehydrogenase 12; RENT2: regulator of nonsense transcripts 2; REV1: DNA repair protein; RNF115_126: E3 ubiquitin-protein ligase; RVB2: RuvB-like protein 2; SLC25A28_37: solute carrier family 25 (mitochondrial iron transporter); SMC4: structural maintenance of chromosome 4; TRAF3IP1: TRAF3-interacting protein 1; TUBB: tubulin beta; U2AF1: splicing factor U2AF 35 kDa subunit; UBE2D: ubiquitin-conjugating enzyme E2 D.

Tab. 2.

Over-expressed genes in female gametophytes.

<h2>Over-expressed genes in female gametophytes.</h2>

Two contigs annotated as myo-inositol-1-phosphate synthase (INO1; K01858) were present in our dataset, homologous to Ectocarpus Ec-26_005440 (with 95 and 74% amino acid identity, respectively). However, while the most homologous copy was equally expressed in male and female gametophytes, the other had higher overall expression but was undetectable in males (Table 2). Amino acid translations indicated that both encode functional full-length INO1 gene copies, based on proteinalignments and InterPro annotations. Furthermore, phylogenetic analysis of available stramenopile INO1 proteins placed both copies in the brown algal clade (S1 File). Given the level of sequence variation between S. latissima INO1 copies, it should in future be possible to verify these data from genomic DNA, and study their expression in more detail using qPCR. INO1 is a key enzyme catalysing the first (and rate limiting) step in the production of inositol-containing compounds from D-glucose 6-phosphate, thereby playing a central role in phospholipid biosynthesis and phosphatidylinositol signalling [60] and various other processes such as phosphate storage in plants, tolerance to abiotic stress, and morphogenesis [6164]. We were, however, unable to identify additional gene expression changes in support of specific downstream pathways involving inositol.

Reactive oxygen species (ROS) are emerging as important signalling molecules that regulate a wide range of physiological and developmental processes, in addition to their potential cytotoxic effects. An important role of ROS signalling during female gametogenesis has been implicated in both mammalian [65] and plant oogenesis [66]. In S. latissima we found two genes with female-specific expression that suggest ROS signalling may also be important in brown algal oogenesis (Table 2); a membrane localized superoxide generating NADPH oxidase (NOX; K21421) was constitutively expressed in female gametophytes, and an Fe-Mn SOD (SOD2; K04564), which showed a 7.5-fold increase in expression after transfer from RL to WL. SOD has been implicated in the regulation of ROS levels in plant megagametogenesis [67]. These observations suggest the working hypothesis that signalling via an oxidative burst may be involved in S. latissima oogenesis.

Perhaps related to putative ROS signalling discussed above, we identified two KEGG genes involved in the prostaglandin (PG) synthesis pathway, prostaglandin-E synthase 2 (PTGES2) and prostaglandin-H2 D-isomerase (HPGDS). Significantly, 9 out of 12 contigs for HPGDS were strongly and constitutively female-biased (Fig 5C–5F; Table 2). PGs are eicosanoids derived enzymatically from 20-carbon polyunsaturated fatty acids (PUFAs), universally present in animals where they have diverse hormonal effects. PG have been detected in brown algae, first identified from Laminariales in the context of stress responses to heavy metal exposure [68]. More recently, gene sequences for PTGES2 and HPGDS were confirmed from diatoms (unicellular heterokont algae) [69]. However, we believe this is the first report of sex-biased expression of HPGDS in algae, which clearly deserves further investigation.

Gender-biased or -specific expression in males

Only three KEGG genes were detected as constitutively overexpressed in male relative to female gametophytes (Fig 5B, 5C–5F). Prominent among these was a high mobility group (HMGB2) protein, represented by a single contig uniquely expressed in males. HMG proteins are transcription factors involved in sex-determination in both animals and fungi [70, 71]. The S. latissima sequence was homologous to Ectocarpus Ec-13_001750, a locus first identified from the male sex-determining region (SDR) in the genome of Ectocarpus [12] and implicated as the potential male-determining factor in these brown algae [13].

Unsurprisingly, an array of genes with mainly late-onset upregulation (days 6 and 8) confirmed the important role of cell proliferation in male gametogenesis (production of large numbers of sperm relative to eggs) at this stage in the time course, and revealed several parallels at the molecular level with plant and animal gametogenesis (Fig 5B, 5C–5F; S2 Table). A total of 4 E3 ubiquitin-protein ligases were over-expressed in males. One of these, SHPRH (K15710) was constitutively male-specific (i.e., TPM < 1 in females; S2 Table; Fig 5C–5F). SHPRH is thought to promote error-free replication in proliferating cells [72], including during mammalian spermatogenesis [73]. The constitutively male-biased expression of this gene may point to a key upstream role in male development. Two of the three other male-biased E3 ubiquitin-protein ligases were members of the RING-type (K110964, K11982), and all have putative roles in cell cycle progression or membrane trafficking. The sex-biased expression of various distinct E3 ubiquitin-protein ligases (with both male- and female-biased members identified) highlights their potential importance as specific and selective regulatory components of sexual development in brown algae, in common with other systems [74]. Several known markers for cell proliferation were male-biased; aurora kinase, a GINS complex subunit, Centrin 3 and MCM3, as well as protein DEK (variously involved in regulation of chromatin structure, epigenetic modification, and transcription) [75]. These data indicate that several common elements of the cell proliferation molecular machinery in multicellular lineages are invoked during spermatogenesis in S. latissima.

The sulfotransferase (ST) GAL3ST3 (galactose-3-O-sulfotransferase 3) was represented by 5 contigs with homology to known brown algal proteins. The most highly expressed member was male-specific in S. latissima, likely reflecting differences in the regulation of cell wall structure between male and female gametophytes. It was shown more than 30 years ago in S. latissima that the antheridial cap, which ruptures to allow the explosive release of sperm in response to female pheromone release, contains sulphated polysaccharides [76] and is therefore a potential site of action for these gene products in males. It was also shown by [76] that antheridial discharge requires Na+, with a rapidity suggesting the possible involvement of a membrane action potential. While not supported statistically, expression of a voltage-gated Na+ channel (K04834; NAV1.2) was detected only in males after 8 days, and is therefore potentially interesting in the context of sperm release.

Expression of flagella-associated genes (mainly) in males

Prominent among male-biased genes expressed after 8 days in WL (Fig 5F) were a number contributing to flagellar development in the antheridia. Several KEGG genes for intraflagellar transport (IFT) and dynein assembly proteins showed low but significantly up-regulated expression after day 8 in WL (Fig 5E and 5F; Table 3), while expression in females was generally very low or not detected, and invariant over the timecourse. The IFT system consists of protein complexes responsible for flagellar assembly and maintenance [77], whose upregulation would be expected during flagellar biogenesis. Several other flagella-related genes (axonemal dyneins, ciliogenesis-related proteins) were upregulated in male gametophytes after day 8 in WL (Table 3), suggesting that antheridial maturation and sperm development were ongoing processes at this time. We comprehensively analysed the expression of flagella-related genes in both female and male gametophyte transcriptomes, based on annotated proteins from the genome of Ectocarpus Ec32 [4] and previous proteomic analysis of brown algal flagella [78]. Of the 93 flagella-related genes identified in the Ectocarpus genome [4], 91 homologues were identified in our S. latissima transcriptome (S4 Table). Expression of flagellar-related genes was highly male-biased, with 45 genes expressed only in males (i.e., female read mapping ≤ 1 TPM, defined as the threshold for expression). However, a further 45 genes were expressed by both sexes, and female-specific contigs were found for 3 genes (S4 Table). In future, it will be important to obtain male and female gametophyte genomes for deep RNAseq mapping to explore sex-specific transcriptional variation and/or post-transcriptional variation (e.g., alternative splice variants).

Tab. 3.

Flagella-related genes over-expressed in male gametophytes.

<h2>Flagella-related genes over-expressed in male gametophytes.</h2>

While flagella are associated chiefly with motile reproductive cells of brown algae, such as sperm, and sexual reproduction in Laminariales is oogamous, it has been known for over 30 years that eggs of Laminariales retain the ability to produce flagella [79]. These lack a function in motility, however, while possessing distinct features (an absence of mastigonemes and divergent basal body structure). Rather than purely vestigial organelles remaining after the evolutionary transition from anisogamy to oogamy, their persistence post-fertilization suggests an important role in anchoring the developing zygote to the oogonium, potentially influencing zygote polarity in early sporophyte development [80]. It remains for future studies to determine the detailed timing and extent of possible female flagella-related gene expression, perhaps at later stages of egg development.

Conclusions

This study provides the first overview of expression differences during gametogenesis in kelps in the context of sex-bias and developmental trajectory. Although preliminary, the results provided an unexpectedly rich picture of extensive changes in transcriptional profiles triggered by gametogenic conditions in S. latissima. Transcriptome profiling revealed that a common immediate response of both sexes to gametogenic conditions overrides sex-specific transcriptional changes, with specific gene annotations suggesting an important role for post-transcriptional and epigenetic regulation of ribosome biogenesis, cell proliferation and differentiation. In this phase, key genes involved in nutrient (N) assimilation and control of intersecting energy and carbon metabolic pathways also showed evidence of regulatory changes.

Consistent with the idea that female development may be the default pathway in brown algal UV sexual systems [12], a number of female-biased or unique genes were identified with roles in gene expression regulation, while males uniquely expressed the putative sex-determining HMG factor protein, hypothesised to repress female development. Male-biased gene expression appears largely to function in coordinating cell proliferation during sperm production, and for sperm flagella biosynthesis. The results also highlight the potential roles of E3 ubiquitin-protein ligases in sex-specific gametogenic pathways, and of ROS signalling in female gametogenesis, both with either ancient or convergent parallels in evolutionarily divergent multicellular lineages.

These and other observations (e.g., on timing and sex-specificity of flagella development), whilst preliminary, suggest a fascinating array of genes, pathways and processes that may be targeted in future functional and comparative studies of brown algal gametogenesis.

Supporting information

S1 File [2017]
Alignment and phylogenetic analysis of Myo-inositol-1-phosphate synthase (INO1) predicted proteins from stramenopiles.

S1 Table [xlsx]
Read data, assembly and annotation statistics.

S2 Table [xlsx]
KEGG-annotated gene lists for differentially expressed transcripts in .

S3 Table [xlsx]
Results of gene set enrichment analysis (GSEA).

S4 Table [xlsx]
Flagella-related proteins in . .


Zdroje

1. van Werven FJ, Amon A. Regulation of entry into gametogenesis. Philos Trans R Soc B Biol Sci. 2011;366(1584):3521–31.

2. Richerd S, Couvet D, Valéro M. Evolution of the alternation of haploid and diploid phases in life cycles. II. Maintenance of the haplo‐diplontic cycle. J Evol Biol. 1993;6(2):263–80.

3. Mable BK, Otto SP. The evolution of life cycles with haploid and diploid phases. Vol. 20, BioEssays. John Wiley & Sons, Inc; 1998.

4. Cock JM, Sterck L, Rouzé P, Scornet D, Allen AE, Amoutzias G, et al. The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature. 2010 Jun;465(7298):617–21. doi: 10.1038/nature09016 20520714

5. Garbary DJ, Kim KY, Klinger T, Duggins D. Red algae as hosts for endophytic kelp gametophytes. Mar Biol. 1999 Oct;135(1):35–40.

6. Bartsch I, Wiencke C, Bischof K, Buchholz CM, Buck BH, Eggert A, et al. The genus Laminaria sensu lato: Recent insights and developments. Eur J Phycol. 2008;43(1):1–86.

7. Martins N, Pearson GA, Gouveia L, Tavares AI, Serrão EA, Bartsch I. Hybrid vigour for thermal tolerance in hybrids between the allopatric kelps Laminaria digitata and L. pallida with contrasting thermal affinities. European Journal of Phycology. 2019; in press.

8. Luning K, Dring MJ. Reproduction, growth and photosynthesis of gametophytes of Laminaria saccharina grown in blue and red light. Mar Biol. 1975;29:195–200.

9. Motomura T, Sakai Y. Ultrastructural studies of gametogenesis in Laminaria angustata (Laminariales, Phaeophyta) regulated by iron concentration in the medium. Phycologia. 1984 Sep 6;23(3):331–43.

10. Müller DG, Gassmann G, Lüning K. Isolation of a spermatozoid-releasing and -attracting substance from female gametophytes of Laminaria digitata. Nature. 1979;279(5712):430–1. doi: 10.1038/279430a0 16068180

11. Müller DG. Relative sexuality in Ectocarpus siliculosus. Arch Microbiol. 1976;109(1–2):89–94.

12. Ahmed S, Cock JM, Pessia E, Luthringer R, Cormier A, Robuchon M, et al. A Haploid System of Sex Determination in the Brown Alga Ectocarpus sp. Curr Biol. 2014 Aug;1–13.

13. Lipinska AP, Toda NRT, Heesch S, Peters AF, Cock JM, Coelho SM. Multiple gene movements into and out of haploid sex chromosomes. Genome Biol. 2017 Dec 8;18(1):104. doi: 10.1186/s13059-017-1201-7 28595587

14. Ye N, Zhang X, Miao M, Fan X, Zheng Y, Xu D, et al. Saccharina genomes provide novel insight into kelp biology. Nat Commun. 2015;6:6986. doi: 10.1038/ncomms7986 25908475

15. Coelho SM, Mignerot L, Cock JM. Origin and evolution of sex-determination systems in the brown algae. New Phytol. 2019;

16. Martins MJF, Mota CF, Pearson GA. Sex-biased gene expression in the brown alga Fucus vesiculosus. BMC Genomics. 2013 Jan;14:294. doi: 10.1186/1471-2164-14-294 23634783

17. Lipinska A, Cormier A, Luthringer R, Peters AF, Corre E, Gachon CMM, et al. Sexual dimorphism and the evolution of sex-biased gene expression in the brown alga Ectocarpus. Mol Biol Evol. 2015;32(6):1581–97. doi: 10.1093/molbev/msv049 25725430

18. Shan TF, Pang SJ, Li J, Li X. De novo transcriptome analysis of the gametophyte of Undaria pinnatifida (Phaeophyceae). J Appl Phycol. 2015;27(2):1011–9.

19. Bartsch I, Vogt J, Pehlke C, Hanelt D. Prevailing sea surface temperatures inhibit summer reproduction of the kelp Laminaria digitata at Helgoland (North Sea). J Phycol. 2013;49(6):1061–73. doi: 10.1111/jpy.12125 27007627

20. Wernberg T, Smale DA, Tuya F, Thomsen MS, Langlois TJ, de Bettignies T, et al. An extreme climatic event alters marine ecosystem structure in a global biodiversity hotspot. Nat Clim Chang. 2013 Jan 22;3(1):78–82.

21. Raybaud V, Beaugrand G, Goberville E, Delebecq G, Destombe C, Valero M, et al. Decline in Kelp in West Europe and Climate. Thrush S, editor. PLoS One. 2013 Jun 26;8(6):e66044. doi: 10.1371/journal.pone.0066044 23840397

22. Robuchon M, Le Gall L, Mauger S, Valero M. Contrasting genetic diversity patterns in two sister kelp species co-distributed along the coast of Brittany, France. Mol Ecol. 2014;23(11):2669–85. doi: 10.1111/mec.12774 24787679

23. Hargrave MS, Foggo A, Pessarrodona A, Smale DA. The effects of warming on the ecophysiology of two co-existing kelp species with contrasting distributions. Oecologia. 2017 Feb 23;183(2):531–43. doi: 10.1007/s00442-016-3776-1 27878385

24. Neiva J, Paulino C, Nielsen MM, Krause-Jensen D, Saunders GW, Assis J, et al. Glacial vicariance drives phylogeographic diversification in the amphi-boreal kelp Saccharina latissima. Sci Rep. 2018 Dec 18;8(1):1112. doi: 10.1038/s41598-018-19620-7 29348650

25. Gerard V a., Du Bois KR. Temperature ecotypes near the southern boundary of the kelp Laminaria saccharina. Mar Biol. 1988 Apr;97(4):575–80.

26. Breton TS, Nettleton JC, O’Connell B, Bertocci M. Fine-scale population genetic structure of sugar kelp, Saccharina latissima (Laminariales, Phaeophyceae), in eastern Maine, USA. Phycologia. 2018 Jan 8;57(1):32–40. doi: 10.2216/17-72.1 29170569

27. Provasoli L. Media and prospects for the cultivation of marine algae. In: Watanabe A, Hattori A, editors. Cultures and Collections of Algae. Proceedings of the US Japan Conference, Hakone. Japanese Society of Plant Physiology, Tokyo; 1966. p. 63–75.

28. Zhang QS, Qu SC, Cong YZ, Luo SJ, Tang XX. High throughput culture and gametogenesis induction of Laminaria japonica gametophyte clones. J Appl Phycol. 2008;20:205–11.

29. Pearson G, Lago-Leston A, Valente M, Serrão E. Simple and rapid RNA extraction from freeze-dried tissue of brown algae and seagrasses. Eur J Phycol. 2006;41:97–104.

30. Brown CT, Howe A, Zhang Q, Pyrkosz AB, Brom TH. A Reference-Free Algorithm for Computational Normalization of Shotgun Sequencing Data. arXiv. 2012 Mar 21;1–18.

31. Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012 Apr 15;28(8):1086–92. doi: 10.1093/bioinformatics/bts094 22368243

32. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015 Jan 17;12(1):59–60. doi: 10.1038/nmeth.3176 25402007

33. Rho M, Tang H, Ye Y. FragGeneScan: predicting genes in short and error-prone reads. Nucleic Acids Res. 2010 Nov;38(20):e191. doi: 10.1093/nar/gkq747 20805240

34. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016 Oct 18;4:e2584–e2584. doi: 10.7717/peerj.2584 27781170

35. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva E V, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015 Jun 9;31(19):3210–2. doi: 10.1093/bioinformatics/btv351 26059717

36. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011 Jan;12:323. doi: 10.1186/1471-2105-12-323 21816040

37. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012 Apr 4;9(4):357–9. doi: 10.1038/nmeth.1923 22388286

38. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010 Jan 1;26(1):139–40. doi: 10.1093/bioinformatics/btp616 19910308

39. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015/01/20. 2015 Apr 20;43(7):e47–e47. doi: 10.1093/nar/gkv007 25605792

40. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004 Mar 19;32(5):1792–7. doi: 10.1093/nar/gkh340 15034147

41. Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000 Apr;17(4):540–52. doi: 10.1093/oxfordjournals.molbev.a026334 10742046

42. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010 May;59(3):307–21. doi: 10.1093/sysbio/syq010 20525638

43. Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–9. doi: 10.1093/bioinformatics/btw313 27207943

44. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York; 2016.

45. Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40. doi: 10.1093/bioinformatics/btu031 24451626

46. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000 Jan 1;28(1):27–30. doi: 10.1093/nar/28.1.27 10592173

47. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. Omi A J Integr Biol. 2012;16(5):284–7.

48. Coelho SM, Godfroy O, Arun A, Le G, Peters AF, Cock JM. OUROBOROS is a master regulator of the gametophyte to sporophyte life cycle transition in the brown alga Ectocarpus. Proc Natl Acad Sci USA. 2011;108(28):11518–23. doi: 10.1073/pnas.1102274108 21709217

49. Arun A, Coelho SM, Peters AF, Bourdareau S, Pérès L, Scornet D, et al. Convergent recruitment of TALE homeodomain life cycle regulators to direct sporophyte development in land plants and brown algae. Elife. 2019;8:1–25.

50. Binda O, Sevilla A, LeRoy G, Lemischka IR, Garcia BA, Richard S. SETD6 monomethylates H2AZ on lysine 7 and is required for the maintenance of embryonic stem cell self-renewal. Epigenetics. 2013 Feb 27;8(2):177–83. doi: 10.4161/epi.23416 23324626

51. Feldman M, Vershinin Z, Goliand I, Elia N, Levy D. The methyltransferase SETD6 regulates Mitotic progression through PLK1 methylation. Proc Natl Acad Sci. 2019 Jan 22;116(4):1235–40. doi: 10.1073/pnas.1804407116 30622182

52. Sanchez CG, Teixeira FK, Czech B, Preall JB, Zamparini AL, Seifert JRK, et al. Regulation of Ribosome Biogenesis and Protein Synthesis Controls Germline Stem Cell Differentiation. Cell Stem Cell. 2016;18(2):276–90. doi: 10.1016/j.stem.2015.11.004 26669894

53. Li H-J, Liu N-Y, Shi D-Q, Liu J, Yang W-C. YAO is a nucleolar WD40-repeat protein critical for embryogenesis and gametogenesis in Arabidopsis. BMC Plant Biol. 2010 Aug;10(1):169.

54. Liu M, Shi D-Q, Yuan L, Liu J, Yang W-C. SLOW WALKER3, Encoding a Putative DEAD-box RNA Helicase, is Essential for Female Gametogenesis in Arabidopsis. J Integr Plant Biol. 2010 Sep 1;52(9):817–28. doi: 10.1111/j.1744-7909.2010.00972.x 20738726

55. Alipanah L, Rohloff J, Winge P, Bones AM, Brembu T. Whole-cell response to nitrogen deprivation in the diatom Phaeodactylum tricornutum. J Exp Bot. 2015;66(20):6281–96. doi: 10.1093/jxb/erv340 26163699

56. Zehr JP, Falkowski PG. Pathway of ammonium assimilation in a marine diatom determined with the radiotracer 13 N. J Phycol. 1988 Dec;24(4):588–91.

57. Guerra LT, Levitan O, Frada MJ, Sun JS, Falkowski PG, Dismukes GC. Regulatory branch points affecting protein and lipid biosynthesis in the diatom Phaeodactylum tricornutum. Biomass and Bioenergy. 2013 Dec;59:306–15.

58. Maxwell DP, Falk S, Trick CG, Huner NPA. Growth at Low Temperature Mimics High-Light Acclimation in Chlorella vulgaris. Plant Physiol. 1994 Jun;105(2):535–43. doi: 10.1104/pp.105.2.535 12232221

59. Cruz C, Ventura F, Bartrons R, Rosa JL. HERC3 binding to and regulation by ubiquitin. FEBS Lett. 2001 Jan 12;488(1–2):74–80. doi: 10.1016/s0014-5793(00)02371-1 11163799

60. Majumder AL, Johnson MD, Henry SA. 1l-myo-Inositol-1-phosphate synthase. Biochim Biophys Acta—Lipids Lipid Metab. 1997;1348(1):245–56.

61. Smart CC, Fleming AJ. A plant gene with homology to D-myo-inositol-3-phosphate synthase is rapidly and spatially up-regulated during an abscisic-acid-induced morphogenic response in Spirodela polyrrhiza. Plant J. 1993 Aug;4(2):279–93. doi: 10.1046/j.1365-313x.1993.04020279.x 8220483

62. Hegeman CE, Good LL, Grabau EA. Expression of d-myo-Inositol-3-Phosphate Synthase in Soybean. Implications for Phytic Acid Biosynthesis. Plant Physiol. 2001;125(4):1941–8. doi: 10.1104/pp.125.4.1941 11299373

63. Donahue JL, Alford SR, Torabinejad J, Kerwin RE, Nourbakhsh A, Ray WK, et al. The Arabidopsis thaliana Myo-Inositol 1-Phosphate Synthase1 Gene Is Required for Myo-inositol Synthesis and Suppression of Cell Death. Plant Cell. 2010;22(3):888–903. doi: 10.1105/tpc.109.071779 20215587

64. Kusuda H, Koga W, Kusano M, Oikawa A, Saito K, Hirai MY, et al. Ectopic expression of myo-inositol 3-phosphate synthase induces a wide range of metabolic changes and confers salt tolerance in rice. Plant Sci. 2015;232:49–56. doi: 10.1016/j.plantsci.2014.12.009 25617323

65. Covarrubias L, Hernández-García D, Schnabel D, Salas-Vidal E, Castro-Obregón S. Function of reactive oxygen species during animal development: Passive or active? Dev Biol. 2008;320(1):1–11. doi: 10.1016/j.ydbio.2008.04.041 18555213

66. Jiménez-Quesada MJ, Traverso JÁ, Alché J de D. NADPH Oxidase-Dependent Superoxide Production in Plant Reproductive Tissues. Front Plant Sci. 2016;7(March):1–13.

67. Martin MV, Fiol DF, Sundaresan V, Zabaleta EJ, Pagnussat GC. oiwa, a Female Gametophytic Mutant Impaired in a Mitochondrial Manganese-Superoxide Dismutase, Reveals Crucial Roles for Reactive Oxygen Species during Embryo Sac Development and Fertilization in Arabidopsis. Plant Cell. 2013;25(5):1573–91. doi: 10.1105/tpc.113.109306 23653473

68. Ritter A, Goulitquer S, Salaün J-P, Tonon T, Correa JA, Potin P. Copper stress induces biosynthesis of octadecanoid and eicosanoid oxygenated derivatives in the brown algal kelp Laminaria digitata. New Phytol. 2008 Dec;180(4):809–21. doi: 10.1111/j.1469-8137.2008.02626.x 18823315

69. Di Dato V, Orefice I, Amato A, Fontanarosa C, Amoresano A, Cutignano A, et al. Animal-like prostaglandins in marine microalgae. ISME J. 2017 Jul 28;11(7):1722–6. doi: 10.1038/ismej.2017.27 28350392

70. Herpin A, Schartl M. Plasticity of gene-regulatory networks controlling sex determination: of masters, slaves, usual suspects, newcomers, and usurpators. EMBO Rep. 2015 Oct 1;16(10):1260–74. doi: 10.15252/embr.201540667 26358957

71. Li W, Sullivan TD, Walton E, Averette AF, Sakthikumar S, Cuomo CA, et al. Identification of the Mating-Type (MAT) Locus That Controls Sexual Reproduction of Blastomyces dermatitidis. Eukaryot Cell. 2013 Jan;12(1):109–17. doi: 10.1128/EC.00249-12 23143684

72. Unk I, Hajdú I, Fátyol K, Szakál B, Blastyák A, Bermudez V, et al. Human SHPRH is a ubiquitin ligase for Mms2-Ubc13-dependent polyubiquitylation of proliferating cell nuclear antigen. Proc Natl Acad Sci. 2006;103(48):18107–12. doi: 10.1073/pnas.0608595103 17108083

73. Sun J, Yomogida K, Sakao S, Yamamoto H, Yoshida K, Watanabe K, et al. Rad18 is required for long-term maintenance of spermatogenesis in mouse testes. Mech Dev. 2009 Mar 1;126(3–4):173–83. doi: 10.1016/j.mod.2008.11.004 19068231

74. Richburg JH, Myers JL, Bratton SB. The role of E3 ligases in the ubiquitin-dependent regulation of spermatogenesis. Semin Cell Dev Biol. 2014 Jun 1;30:27–35. doi: 10.1016/j.semcdb.2014.03.001 24632385

75. Smith EA, Krumpelbeck EF, Jegga AG, Prell M, Matrka MM, Kappes F, et al. The nuclear DEK interactome supports multi-functionality. Proteins Struct Funct Bioinforma. 2018 Jan;86(1):88–97.

76. Maier I. New aspects of pheromone-triggered spermatozoid release in Laminaria digitata (Phaeophyta). Protoplasma. 1982;113(2):137–43.

77. Rosenbaum JL, Witman GB. Intraflagellar transport. Nat Rev Mol Cell Biol. 2002;3(11):813–25. doi: 10.1038/nrm952 12415299

78. Fu G, Nagasato C, Oka S, Cock JM, Motomura T. Proteomics Analysis of Heterogeneous Flagella in Brown Algae (Stramenopiles). Protist. 2014;165(5):662–75. doi: 10.1016/j.protis.2014.07.007 25150613

79. Motomura T, Sakai Y. The occurrence of flagellated eggs in Laminaria angustata (Phaeophyta, Laminariales). J Phycol. 1988;24:282–5.

80. Klochkova TA, Motomura T, Nagasato C, Klimova A V., Kim GH. The role of egg flagella in the settlement and development of zygotes in two Saccharina species. Phycologia. 2019 Mar 4;58(2):145–53.


Článek vyšel v časopise

PLOS One


2019 Číslo 9