Resting cells of Skeletonema marinoi assimilate organic compounds and respire by dissimilatory nitrate reduction to ammonium in dark, anoxic conditions.

Diatoms can survive long periods in dark, anoxic sediments by forming resting spores or resting cells. These have been considered dormant until recently when resting cells of Skeletonema marinoi were shown to assimilate nitrate and ammonium from the ambient environment in dark, anoxic conditions. Here, we show that resting cells of S. marinoi can also perform dissimilatory nitrate reduction to ammonium (DNRA), in dark, anoxic conditions. Transmission electron microscope analyses showed that chloroplasts were compacted, and few large mitochondria had visible cristae within resting cells. Using secondary ion mass spectrometry and isotope ratio mass spectrometry combined with stable isotopic tracers, we measured assimilatory and dissimilatory processes carried out by resting cells of S. marinoi under dark, anoxic conditions. Nitrate was both respired by DNRA and assimilated into biomass by resting cells. Cells assimilated nitrogen from urea and carbon from acetate, both of which are sources of dissolved organic matter produced in sediments. Carbon and nitrogen assimilation rates corresponded to turnover rates of cellular carbon and nitrogen content ranging between 469 and 10,000 years. Hence, diatom resting cells can sustain their cells in dark, anoxic sediments by slowly assimilating and respiring substrates from the ambient environment.

Strain-specific metabarcoding reveals rapid evolution of copper tolerance in populations of the coastal diatom Skeletonema marinoi.

Phytoplankton have short generation times, flexible reproduction strategies, large population sizes and high standing genetic diversity, traits that should facilitate rapid evolution under directional selection. We quantified local adaptation of copper tolerance in a population of the diatom Skeletonema marinoi from a mining-exposed inlet in the Baltic Sea and in a non-exposed population 100 km away. We hypothesized that mining pollution has driven evolution of elevated copper tolerance in the impacted population of S. marinoi. Assays of 58 strains originating from sediment resting stages revealed no difference in the average tolerance to copper between the two populations. However, variation within populations was greater at the mining site, with three strains displaying hyper-tolerant phenotypes. In an artificial evolution experiment, we used a novel intraspecific metabarcoding locus to track selection and quantify fitness of all 58 strains during co-cultivation in one control and one toxic copper treatment. As expected, the hyper-tolerant strains enabled rapid evolution of copper tolerance in the mining-exposed population through selection on available strain diversity. Within 42 days, in each experimental replicate a single strain dominated (30%-99% abundance) but different strains dominated the different treatments. The reference population developed tolerance beyond expectations primarily due to slowly developing plastic response in one strain, suggesting that different modes of copper tolerance are present in the two populations. Our findings provide novel empirical evidence that standing genetic diversity of phytoplankton resting stage allows populations to evolve rapidly (20-50 generations) and flexibly on timescales relevant for seasonal bloom progressions.

Chloroplast magnesium transporters play essential but differential roles in maintaining magnesium homeostasis.

Magnesium (Mg2+) is essential for photosynthesis in the chloroplasts of land plants and algae. Being the central ion of chlorophyll, cofactor and activator of many photosynthetic enzymes including RuBisCO, magnesium-deficient plants may suffer from leaf chlorosis symptoms and retarded growth. Therefore, the chloroplast Mg2+ concentration is tightly controlled by magnesium transport proteins. Recently, three different transporters from two distinct families have been identified in the chloroplast inner envelope of the model plant Arabidopsis thaliana: MGT10, MGR8, and MGR9. Here, we assess the individual roles of these three proteins in maintaining chloroplast Mg2+ homeostasis and regulating photosynthesis, and if their role is conserved in the model green alga Chlamydomonas reinhardtii. Phylogenetic analysis and heterologous expression revealed that the CorC-like MGR8 and MGR9 transport Mg2+ by a different mechanism than the CorA-like MGT10. MGR8 and MGT10 genes are highest expressed in leaves, indicating a function in chloroplast Mg2+ transport. MGR9 is important for chloroplast function and plant adaptation in conditions of deficiency or excess of Mg2+. Transmission electron microscopy indicated that MGT10 plays a differential role in thylakoid stacking than MGR8 and MGR9. Furthermore, we report that MGR8, MGR9, and MGT10 are involved in building up the pH gradient across the thylakoid membrane and activating photoprotection in conditions of excess light, however the mechanism has not been resolved yet. While there are no chloroplast MGR-like transporters in Chlamydomonas, we show that MRS4 is a homolog of MGT10, that is required for photosynthesis and cell growth. Taken together, our findings reveal that the studied Mg2+ transporters play essential but differential roles in maintaining chloroplast Mg2+ homeostasis.

Clones on the run: The genomics of a recently expanded partially clonal species.

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

A metabarcoding analysis of the wrackbed microbiome indicates a phylogeographic break along the North Sea-Baltic Sea transition zone.

Sandy beaches are biogeochemical hotspots that bridge marine and terrestrial ecosystems via the transfer of organic matter, such as seaweed (termed wrack). A keystone of this unique ecosystem is the microbial community, which helps to degrade wrack and re-mineralize nutrients. However, little is known about this community. Here, we characterize the wrackbed microbiome as well as the microbiome of a primary consumer, the seaweed fly Coelopa frigida, and examine how they change along one of the most studied ecological gradients in the world, the transition from the marine North Sea to the brackish Baltic Sea. We found that polysaccharide degraders dominated both microbiomes, but there were still consistent differences between wrackbed and fly samples. Furthermore, we observed a shift in both microbial communities and functionality between the North and Baltic Sea driven by changes in the frequency of different groups of known polysaccharide degraders. We hypothesize that microbes were selected for their abilities to degrade different polysaccharides corresponding to a shift in polysaccharide content in the different seaweed communities. Our results reveal the complexities of both the wrackbed microbial community, with different groups specialized to different roles, and the cascading trophic consequences of shifts in the near shore algal community.

Improving the genome and proteome annotations of the marine model diatom Thalassiosira pseudonana using a proteogenomics strategy.

Diatoms are unicellular eukaryotic phytoplankton that account for approximately 20% of global carbon fixation and 40% of marine primary productivity; thus, they are essential for global carbon biogeochemical cycling and climate. The availability of ten diatom genome sequences has facilitated evolutionary, biological and ecological research over the past decade; however, a complimentary map of the diatom proteome with direct measurements of proteins and peptides is still lacking. Here, we present a proteome map of the model marine diatom Thalassiosira pseudonana using high-resolution mass spectrometry combined with a proteogenomic strategy. In-depth proteomic profiling of three different growth phases and three nutrient-deficient samples identified 9526 proteins, accounting for ~ 81% of the predicted protein-coding genes. Proteogenomic analysis identified 1235 novel genes, 975 revised genes, 104 splice variants and 234 single amino acid variants. Furthermore, our quantitative proteomic analysis experimentally demonstrated that a considerable number of novel genes were differentially translated under different nutrient conditions. These findings substantially improve the genome annotation of T. pseudonana and provide insights into new biological functions of diatoms. This relatively comprehensive diatom proteome catalog will complement available diatom genome and transcriptome data to advance biological and ecological research of marine diatoms. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-022-00161-y.

MSALigMap-A Tool for Mapping Active-Site Amino Acids in PDB Structures onto Known and Novel Unannotated Homologous Sequences with Similar Function.

MSALigMap (Multiple Sequence Alignment Ligand Mapping) is a tool for mapping active-site amino-acid residues that bind selected ligands on to target protein sequences of interest. Users can also provide novel sequences (unavailable in public databases) for analysis. MSALigMap is written in Python. There are several tools and servers available for comparing and mapping active-site amino-acid residues among protein structures. However, there has not previously been a tool for mapping ligand binding amino-acid residues onto protein sequences of interest. Using MSALigMap, users can compare multiple protein sequences, such as those from different organisms or clinical strains, with sequences of proteins with crystal structures in PDB that are bound with the ligand/drug and DNA of interest. This allows users to easily map the binding residues and to predict the consequences of different mutations observed in the binding site. The MSALigMap server can be accessed at https://albiorix.bioenv.gu.se/MSALigMap/HomePage.py.

Strain-specific transcriptional responses overshadow salinity effects in a marine diatom sampled along the Baltic Sea salinity cline.

The salinity gradient separating marine and freshwater environments represents a major ecological divide for microbiota, yet the mechanisms by which marine microbes have adapted to and ultimately diversified in freshwater environments are poorly understood. Here, we take advantage of a natural evolutionary experiment: the colonization of the brackish Baltic Sea by the ancestrally marine diatom Skeletonema marinoi. To understand how diatoms respond to low salinity, we characterized transcriptomic responses of acclimated S. marinoi grown in a common garden. Our experiment included eight strains from source populations spanning the Baltic Sea salinity cline. Gene expression analysis revealed that low salinities induced changes in the cellular metabolism of S. marinoi, including upregulation of photosynthesis and storage compound biosynthesis, increased nutrient demand, and a complex response to oxidative stress. However, the strain effect overshadowed the salinity effect, as strains differed significantly in their response, both regarding the strength and the strategy (direction of gene expression) of their response. The high degree of intraspecific variation in gene expression observed here highlights an important but often overlooked source of biological variation associated with how diatoms respond to environmental change.

A highly contiguous, scaffold-level nuclear genome assembly for the fever tree (Cinchona pubescens Vahl) as a novel resource for Rubiaceae research.

The Andean fever tree (Cinchona L.; Rubiaceae) is a source of bioactive quinine alkaloids used to treat malaria. C. pubescens Vahl is a valuable cash crop within its native range in northwestern South America, however, genomic resources are lacking. Here we provide the first highly contiguous and annotated nuclear and plastid genome assemblies using Oxford Nanopore PromethION-derived long-read and Illumina short-read data. Our nuclear genome assembly comprises 603 scaffolds with a total length of 904 Mbp (∼82% of the full genome based on a genome size of 1.1 Gbp/1C). Using a combination of de novo and reference-based transcriptome assemblies we annotated 72,305 coding sequences comprising 83% of the BUSCO gene set and 4.6% fragmented sequences. Using additional plastid and nuclear datasets we place C. pubescens in the Gentianales order. This first genomic resource for C. pubescens opens new research avenues, including the analysis of alkaloid biosynthesis in the fever tree.

Ligand Binding Site Comparison - LiBiSCo - a web-based tool for analyzing interactions between proteins and ligands to explore amino acid specificity within active sites.

Interaction between protein and ligands are ubiquitous in a biological cell, and understanding these interactions at the atom level in protein-ligand complexes is crucial for structural bioinformatics and drug discovery. Here, we present a web-based protein-ligand interaction application named Ligand Binding Site Comparison (LiBiSCo) for comparing the amino acid residues interacting with atoms of a ligand molecule between different protein-ligand complexes available in the Protein Data Bank (PDB) database. The comparison is performed at the ligand atom level irrespectively of having binding site similarity or not between the protein structures of interest. The input used in LiBiSCo is one or several PDB IDs of protein-ligand complex(es) and the tool returns a list of identified interactions at ligand atom level including both bonded and non-bonded interactions. A sequence profile for the interaction for each ligand atoms is provided as a WebLogo. The LiBiSco is useful in understanding ligand binding specificity and structural promiscuity among families that are structurally unrelated. The LiBiSCo tool can be accessed through https://albiorix.bioenv.gu.se/LiBiSCo/HomePage.py.

Post-glacial establishment of locally adapted fish populations over a steep salinity gradient.

Studies of colonization of new habitats that appear from rapidly changing environments are interesting and highly relevant to our understanding of divergence and speciation. Here, we analyse phenotypic and genetic variation involved in the successful establishment of a marine fish (sand goby, Pomatoschistus minutus) over a steep salinity drop from 35 PSU in the North Sea (NE Atlantic) to two PSU in the inner parts of the post-glacial Baltic Sea. We first show that populations are adapted to local salinity in a key reproductive trait, the proportion of motile sperm. Thereafter, we show that genome variation at 22,190 single nucleotide polymorphisms (SNPs) shows strong differentiation among populations along the gradient. Sequences containing outlier SNPs and transcriptome sequences, mapped to a draft genome, reveal associations with genes with relevant functions for adaptation in this environment but without overall evidence of functional enrichment. The many contigs involved suggest polygenic differentiation. We trace the origin of this differentiation using demographic modelling and find the most likely scenario is that at least part of the genetic differentiation is older than the Baltic Sea and is a result of isolation of two lineages prior to the current contact over the North Sea-Baltic Sea transition zone.

Deep learning suggests that gene expression is encoded in all parts of a co-evolving interacting gene regulatory structure.

Understanding the genetic regulatory code governing gene expression is an important challenge in molecular biology. However, how individual coding and non-coding regions of the gene regulatory structure interact and contribute to mRNA expression levels remains unclear. Here we apply deep learning on over 20,000 mRNA datasets to examine the genetic regulatory code controlling mRNA abundance in 7 model organisms ranging from bacteria to Human. In all organisms, we can predict mRNA abundance directly from DNA sequence, with up to 82% of the variation of transcript levels encoded in the gene regulatory structure. By searching for DNA regulatory motifs across the gene regulatory structure, we discover that motif interactions could explain the whole dynamic range of mRNA levels. Co-evolution across coding and non-coding regions suggests that it is not single motifs or regions, but the entire gene regulatory structure and specific combination of regulatory elements that define gene expression levels.

The round goby genome provides insights into mechanisms that may facilitate biological invasions.

BACKGROUND: The invasive benthic round goby (Neogobius melanostomus) is the most successful temperate invasive fish and has spread in aquatic ecosystems on both sides of the Atlantic. Invasive species constitute powerful in situ experimental systems to study fast adaptation and directional selection on short ecological timescales and present promising case studies to understand factors involved the impressive ability of some species to colonize novel environments. We seize the unique opportunity presented by the round goby invasion to study genomic substrates potentially involved in colonization success. RESULTS: We report a highly contiguous long-read-based genome and analyze gene families that we hypothesize to relate to the ability of these fish to deal with novel environments. The analyses provide novel insights from the large evolutionary scale to the small species-specific scale. We describe expansions in specific cytochrome P450 enzymes, a remarkably diverse innate immune system, an ancient duplication in red light vision accompanied by red skin fluorescence, evolutionary patterns of epigenetic regulators, and the presence of osmoregulatory genes that may have contributed to the round goby's capacity to invade cold and salty waters. A recurring theme across all analyzed gene families is gene expansions. CONCLUSIONS: The expanded innate immune system of round goby may potentially contribute to its ability to colonize novel areas. Since other gene families also feature copy number expansions in the round goby, and since other Gobiidae also feature fascinating environmental adaptations and are excellent colonizers, further long-read genome approaches across the goby family may reveal whether gene copy number expansions are more generally related to the ability to conquer new habitats in Gobiidae or in fish.

Development and characterization of an EMS-mutagenized wheat population and identification of salt-tolerant wheat lines.

BACKGROUND: Triticum aestivum (wheat) is one of the world's oldest crops and has been used for >8000 years as a food crop in North Africa, West Asia and Europe. Today, wheat is one of the most important sources of grain for humans, and is cultivated on greater areas of land than any other crop. As the human population increases and soil salinity becomes more prevalent, there is increased pressure on wheat breeders to develop salt-tolerant varieties in order to meet growing demands for yield and grain quality. Here we developed a mutant wheat population using the moderately salt-tolerant Bangladeshi variety BARI Gom-25, with the primary goal of further increasing salt tolerance. RESULTS: After titrating the optimal ethyl methanesulfonate (EMS) concentration, ca 30,000 seeds were treated with 1% EMS, and 1676 lines, all originating from single seeds, survived through the first four generations. Most mutagenized lines showed a similar phenotype to BARI Gom-25, although visual differences such as dwarfing, giant plants, early and late flowering and altered leaf morphology were seen in some lines. By developing an assay for salt tolerance, and by screening the mutagenized population, we identified 70 lines exhibiting increased salt tolerance. The selected lines typically showed a 70% germination rate on filter paper soaked in 200 mM NaCl, compared to 0-30% for BARI Gom-25. From two of the salt-tolerant OlsAro lines (OA42 and OA70), genomic DNA was sequenced to 15x times coverage. A comparative analysis against the BARI Gom-25 genomic sequence identified a total of 683,201 (OA42), and 768,954 (OA70) SNPs distributed throughout the three sub-genomes (A, B and D). The mutation frequency was determined to be approximately one per 20,000 bp. All the 70 selected salt-tolerant lines were tested for root growth in the laboratory, and under saline field conditions in Bangladesh. The results showed that all the lines selected for tolerance showed a better salt tolerance phenotype than both BARI Gom-25 and other local wheat varieties tested. CONCLUSION: The mutant wheat population developed here will be a valuable resource in the development of novel salt-tolerant varieties for the benefit of saline farming.

Phenomics reveals a novel putative chloroplast fatty acid transporter in the marine diatom Skeletonema marinoi involved in temperature acclimation.

Diatoms are the dominant phytoplankton in temperate oceans and coastal regions and yet little is known about the genetic basis underpinning their global success. Here, we address this challenge by developing the first phenomic approach for a diatom, screening a collection of randomly mutagenized but identifiably tagged transformants. Based upon their tolerance to temperature extremes, several compromised mutants were identified revealing genes either stress related or encoding hypothetical proteins of unknown function. We reveal one of these hypothetical proteins is a novel putative chloroplast fatty acid transporter whose loss affects several fatty acids including the two omega-3, long-chain polyunsaturated fatty acids - eicosapentaenoic and docosahexaenoic acid, both of which have medical importance as dietary supplements and industrial significance in aquaculture and biofuels. This mutant phenotype not only provides new insights into the fatty acid biosynthetic pathways in diatoms but also highlights the future value of phenomics for revealing specific gene functions in these ecologically important phytoplankton.

Whole Genome Sequence of Marinobacter salarius Strain SMR5, Shown to Promote Growth in its Diatom Host.

Attempts to obtain axenic cultures of the marine diatom Skeletonema marinoi often result in poor growth, indicating the importance of the microbiome to the growth of its host. In order to identify the precise roles played by these associated bacteria, individual strains were isolated, cultured and sequenced. We report the genome of one such strain - SMR5, isolated from a culture of S. marinoi strain R05AC sampled from top layer sediments of the Swedish west coast. Its genome of 4,630,160 bp consists of a circular chromosome and one circular plasmid, and 4,263 CDSs were inferred in the annotation. Comparison of 16S rRNA sequences and other markers, along with phylotaxonomic analysis, leads us to place strain SMR5 in the taxon Marinobacter salarius. Pathway analysis and previous experimental work suggest that this strain may produce a growth factor, as well as improve iron availability for its host via siderophores.

Friends With Benefits: Exploring the Phycosphere of the Marine Diatom Skeletonema marinoi.

Marine diatoms are the dominant phytoplankton in the temperate oceans and coastal regions, contributing to global photosynthesis, biogeochemical cycling of key nutrients and minerals and aquatic food chains. Integral to the success of marine diatoms is a diverse array of bacterial species that closely interact within the diffusive boundary layer, or phycosphere, surrounding the diatom partner. Recently, we isolated seven distinct bacterial species from cultures of Skeletonema marinoi, a chain-forming, centric diatom that dominates the coastal regions of the temperate oceans. Genomes of all seven bacteria were sequenced revealing many unusual characteristics such as the existence of numerous plasmids of widely varying sizes. Here we have investigated the characteristics of the bacterial interactions with S. marinoi, demonstrating that several strains (Arenibacter algicola strain SMS7, Marinobacter salarius strain SMR5, Sphingorhabdus flavimaris strain SMR4y, Sulfitobacter pseudonitzschiae strain SMR1, Yoonia vestfoldensis strain SMR4r and Roseovarius mucosus strain SMR3) stimulate growth of the diatom partner. Testing of many different environmental factors including low iron concentration, high and low temperatures, and chemical signals showed variable effects on this growth enhancement by each bacterial species, with the most significant being light quality in which green and blue but not red light enhanced the stimulatory effect on S. marinoi growth by all bacteria. Several of the bacteria also inhibited growth of one or more of the other bacterial strains to different extents when mixed together. This study highlights the complex interactions between diatoms and their associated bacteria within the phycosphere, and that further studies are needed to resolve the underlying mechanisms for these relationships and how they might influence the global success of marine diatoms.

Complete Genome Sequence of the Diatom-Associated Bacterium Sphingorhabdus sp. Strain SMR4y.

The bacterial strain SMR4y belongs to the diverse microbiome of the marine diatom Skeletonema marinoi strain R05AC. After assembly of its genome, presented here, and subsequent analyses, we placed it in the genus Sphingorhabdus This strain has a 3,479,724-bp circular chromosome (with 3,340 coding sequences) and no known plasmids.

Exploring Molecular Signs of Sex in the Marine Diatom Skeletonema marinoi.

Sexual reproduction plays a fundamental role in diatom life cycles. It contributes to increasing genetic diversity through meiotic recombination and also represents the phase where large-sized cells are produced to counteract the cell size reduction process that characterizes these microalgae. With the aim to identify genes linked to the sexual phase of the centric planktonic diatom Skeletonemamarinoi, we carried out an RNA-seq experiment comparing the expression level of transcripts in sexualized cells with that of large cells not competent for sex. A set of genes involved in meiosis were found upregulated. Despite the fact that flagellate gametes were observed in the sample, we did not detect the expression of genes involved in the synthesis of flagella that were upregulated during sexual reproduction in another centric diatom. A comparison with the set of genes changing during the first phases of sexual reproduction of the pennate diatom Pseudo-nitzschiamultistriata revealed the existence of commonalities, including the strong upregulation of genes with an unknown function that we named Sex Induced Genes (SIG). Our results further broadened the panel of genes that can be used as a marker for sexual reproduction of diatoms, crucial for the interpretation of metatranscriptomic datasets.

Genome Sequence of Kordia sp. Strain SMS9 Identified in a Non-Axenic Culture of the Diatom Skeletonema marinoi.

Initial efforts to sequence the genome of the marine diatom Skeletonema marinoi were hampered by the presence of genetic material from bacteria, and there was sufficient material from some of these bacteria to enable the assembly of full chromosomes. Here, we report the genome of strain SMS9, one such bacterial species identified in a non-axenic culture of S. marinoi strain ST54. Its 5,482,391 bp circular chromosome contains 4,641 CDSs, and has a G+C content of 35.6%. Based on 16S rRNA comparison, phylotaxonomic analysis, and the genome similarity metrics dDDH and OrthoANI, we place this strain in the genus Kordia, and to the best of our knowledge, this is the first Kordia species to be initially described from European waters. As attempts to culture this strain have failed, however, the specifics of its relationship with S. marinoi are still uncertain.