Species specific gene expression dynamics during harmful algal blooms

Sexual reproduction during diatom bloom

Phytoplankton supports food webs in all aquatic ecosystems. Ecological studies highlighted the links between environmental variables and species successions in situ. However, the role of life cycle characteristics on phytoplankton community dynamics remains poorly characterized. In diatoms, sexual reproduction creates new genetic combinations and prevents excessive cell size miniaturization. It has been extensively studied in vitro but seldom in the natural environment. Here, analyzing metatranscriptomic data in the light of the expression patterns previously characterized in vitro, we identified a synchronized and transient sexual reproduction event during a bloom of the toxic diatom species Pseudo-nitzschia australis. Despite the complexity of environmental conditions encountered in situ, sexual reproduction appeared to be the strongest differential gene expression signal that occurred during the bloom. The potential link between environmental conditions and the initiation of sexual reproduction remain to be determined, but sexual reproduction probably had a major impact on the bloom dynamic.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

Marine heterotrophic flagellates (HF) are dominant bacterivores in the ocean, where they represent the trophic link between bacteria and higher trophic levels and participate in the recycling of inorganic nutrients for regenerated primary production. Studying their activity and function in the ecosystem is challenging since most of the HFs in the ocean are still uncultured. In the present work, we investigated gene expression of natural HF communities during bacterivory in four unamended seawater incubations.

RESULTS

The most abundant species growing in our incubations belonged to the taxonomic groups MAST-4, MAST-7, Chrysophyceae, and Telonemia. Gene expression dynamics were similar between incubations and could be divided into three states based on microbial counts, each state displaying distinct expression patterns. The analysis of samples where HF growth was highest revealed some highly expressed genes that could be related to bacterivory. Using available genomic and transcriptomic references, we identified 25 species growing in our incubations and used those to compare the expression levels of these specific genes. Video Abstract CONCLUSIONS: Our results indicate that several peptidases, together with some glycoside hydrolases and glycosyltransferases, are more expressed in phagotrophic than in phototrophic species, and thus could be used to infer the process of bacterivory in natural assemblages.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Aureococcus anophagefferens (Pelagophyceae) genomes improve evaluation of nutrient acquisition strategies involved in brown tide dynamics

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment.

De novo transcriptome assembly of the green alga Ankistrodesmus falcatus

Ankistrodesmus falcatus is a globally distributed freshwater chlorophyte that is a candidate for biofuel production, is used to study the effects of toxins on aquatic communities, and is used as food in zooplankton research. Each of these research fields is transitioning to genomic tools. We created a reference transcriptome for of A. falcatus using NextGen sequencing and de novo assembly methods including Trinity, Velvet-Oases, and EvidentialGene. The assembled transcriptome has a total of 17,997 contigs, an N50 value of 2,462, and a GC content of 64.8%. BUSCO analysis recovered 83.3% of total chlorophyte BUSCOs and 82.5% of the eukaryotic BUSCOs. A portion (7.9%) of these supposedly single-copy genes were found to have transcriptionally active, distinct duplicates. We annotated the assembly using the dammit annotation pipeline, resulting in putative functional annotation for 68.89% of the assembly. Using available rbcL sequences from 16 strains (10 species) of Ankistrodesmus, we constructed a neighbor-joining phylogeny to illustrate genetic distances of our A. falcatus strain to other members of the genus. This assembly will be valuable for researchers seeking to identify Ankistrodesmus sequences in metatranscriptomic and metagenomic field studies and in experiments where separating expression responses of zooplankton and their algal food sources through bioinformatics is important.