Gene Duplication, Shifting Selection, and Dosage Balance of Silicon Transporter Proteins in Marine and Freshwater Diatoms.

Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine-freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1-2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90-100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance.

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

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.

The genome of a nonphotosynthetic diatom provides insights into the metabolic shift to heterotrophy and constraints on the loss of photosynthesis.

Although most of the tens of thousands of diatom species are photoautotrophs, a small number of heterotrophic species no longer photosynthesize. We sequenced the genome of a nonphotosynthetic diatom, Nitzschia Nitz4, to determine how carbon metabolism was altered in the wake of this trophic shift. Nitzschia Nitz4 has retained its plastid and plastid genome, but changes associated with the transition to heterotrophy were cellular-wide and included losses of photosynthesis-related genes from the nuclear and plastid genomes, elimination of isoprenoid biosynthesis in the plastid, and remodeling of mitochondrial glycolysis to maximize adenosine triphosphte (ATP) yield. The genome contains a ß-ketoadipate pathway that may allow Nitzschia Nitz4 to metabolize lignin-derived compounds. Diatom plastids lack an oxidative pentose phosphate pathway (oPPP), leaving photosynthesis as the primary source of NADPH to support essential biosynthetic pathways in the plastid and, by extension, limiting available sources of NADPH in nonphotosynthetic plastids. The genome revealed similarities between nonphotosynthetic diatoms and apicomplexan parasites for provisioning NADPH in their plastids and highlighted the ancestral absence of a plastid oPPP as a potentially important constraint on loss of photosynthesis, a hypothesis supported by the higher frequency of transitions to parasitism or heterotrophy in lineages that have a plastid oPPP.

Transcriptional Response of Osmolyte Synthetic Pathways and Membrane Transporters in a Euryhaline Diatom During Long-term Acclimation to a Salinity Gradient.

How diatoms respond to fluctuations in osmotic pressure is important from both ecological and applied perspectives. It is well known that osmotic stress affects photosynthesis and can result in the accumulation of compounds desirable in pharmaceutical and alternative fuel industries. Gene expression responses to osmotic stress have been studied in short-term trials, but it is unclear whether the same mechanisms are recruited during long-term acclimation. We used RNA-seq to study the genome-wide transcription patterns in the euryhaline diatom, Cyclotella cryptica, following long-term acclimation to salinity that spanned the natural range of fresh to oceanic water. Long-term acclimated C. cryptica exhibited induced synthesis or repressed degradation of the osmolytes glycine betaine, taurine and dimethylsulfoniopropionate (DMSP). Although changes in proline concentration is one of the main responses in short-term osmotic stress, we did not detect a transcriptional change in proline biosynthetic pathways in our long-term experiment. Expression of membrane transporters showed a general tendency for increased import of potassium and export of sodium, consistent with the electrochemical gradients and dependence on co-transported molecules. Our results show substantial between-genotype differences in growth and gene expression reaction norms and suggest that the regulation of proline synthesis important in short-term osmotic stress might not be maintained in long-term acclimation. Further examination using time-course gene expression experiments, metabolomics and genetic validation of gene functions would reinforce patterns inferred from RNA-seq data.

Improved Reference Genome for Cyclotella cryptica CCMP332, a Model for Cell Wall Morphogenesis, Salinity Adaptation, and Lipid Production in Diatoms (Bacillariophyta).

The diatom, Cyclotella cryptica, is a well-established model species for physiological studies and biotechnology applications of diatoms. To further facilitate its use as a model diatom, we report an improved reference genome assembly and annotation for C. cryptica strain CCMP332. We used a combination of long- and short-read sequencing to assemble a high-quality and contaminant-free genome. The genome is 171 Mb in size and consists of 662 scaffolds with a scaffold N50 of 494 kb. This represents a 176-fold decrease in scaffold number and 41-fold increase in scaffold N50 compared to the previous assembly. The genome contains 21,250 predicted genes, 75% of which were assigned putative functions. Repetitive DNA comprises 59% of the genome, and an improved classification of repetitive elements indicated that a historically steady accumulation of transposable elements has contributed to the relatively large size of the C. cryptica genome. The high-quality C. cryptica genome will serve as a valuable reference for ecological, genetic, and biotechnology studies of diatoms.

A single loss of photosynthesis in the diatom order Bacillariales (Bacillariophyta).

PREMISE OF THE STUDY: Loss of photosynthesis is a common and often repeated trajectory in nearly all major groups of photosynthetic eukaryotes. One small subset of "apochloritic" diatoms in the genus Nitzschia have lost their ability to photosynthesize and require extracellular carbon for growth. Similar to other secondarily nonphotosynthetic taxa, apochloritic diatoms maintain colorless plastids with highly reduced plastid genomes. Although the narrow taxonomic breadth of apochloritic Nitzschia suggests a single loss of photosynthesis in their common ancestor, previous phylogenetic analyses suggested that photosynthesis was lost multiple times. METHODS: We analyzed genes from the nuclear, plastid, and mitochondrial genomes for a broad set of taxa to test whether photosynthesis was lost one or multiple times in Bacillariales. We also sequenced and characterized the plastid genome of a nonphotosynthetic Nitzschia species. KEY RESULTS: Phylogenetic analyses showed that genes from all three genetic compartments either supported or failed to reject monophyly of apochloritic Nitzschia species, consistent with a single loss of photosynthesis in this group. The plastid genomes of two apochloritic Nitzschia are highly similar in all respects, indicating streamlining of the plastid genome before the split of these two species. CONCLUSIONS: A better understanding of the phylogeny and ecology of apochloritic Nitzschia, together with emerging genomic resources, will help identify the factors that have driven and maintained the loss of photosynthesis in this group of diatoms. Finally, some habitats host diverse communities of co-occurring nonphotosynthetic diatoms, reflecting resource abundance or resource partitioning in ecologically favorable habitats.

Recurrent Loss, Horizontal Transfer, and the Obscure Origins of Mitochondrial Introns in Diatoms (Bacillariophyta).

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

Phylogenomics reveals an extensive history of genome duplication in diatoms (Bacillariophyta).

PREMISE OF THE STUDY: Diatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and primary productivity motivate comparisons to angiosperms, whose genomes have been inordinately shaped by whole-genome duplication (WGD). WGDs have been linked to speciation, increased rates of lineage diversification, and identified as a principal driver of angiosperm evolution. We synthesized a large but scattered body of evidence that suggests polyploidy may be common in diatoms as well. METHODS: We used gene counts, gene trees, and distributions of synonymous divergence to carry out a phylogenomic analysis of WGD across a diverse set of 37 diatom species. KEY RESULTS: Several methods identified WGDs of varying age across diatoms. Determining the occurrence, exact number, and placement of events was greatly impacted by uncertainty in gene trees. WGDs inferred from synonymous divergence of paralogs varied depending on how redundancy in transcriptomes was assessed, gene families were assembled, and synonymous distances (Ks) were calculated. Our results highlighted a need for systematic evaluation of key methodological aspects of Ks-based approaches to WGD inference. Gene tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades. CONCLUSIONS: Our results suggest that WGD has played a major role in the evolution of diatom genomes. We outline challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the impact of genome duplication in a group that likely harbors substantial genomic diversity.

Morphological diversity and phylogeny of the diatom genus Entomoneis (Bacillariophyta) in marine plankton: six new species from the Adriatic Sea.

The diatom genus Entomoneis is known from the benthos and plankton of marine, brackish, and freshwaters. Entomoneis includes diatoms with a bilobate keel elevated above the valve surface, a sigmoid canal raphe, and numerous girdle bands. Owing mostly to the scarcity of molecular data for a diverse set of species, the phylogeny of Entomoneis has not been investigated in depth. The few previous studies that included Entomoneis were focused on broader questions and the available data were from a small number of either unidentified Entomoneis or well-known species (e.g., E. paludosa). Since the first description of new species combining both molecular and morphological characters (E. tenera), we have continued to cultivate and investigate Entomoneis in the plankton of the Adriatic Sea. Combined multigene phylogeny (SSU rDNA sequences, rbcL, and psbC genes) and morphological observations (LM, SEM and TEM) revealed six new Entomoneis species supported by phylogenetic and morphological data: E. pusilla, E. gracilis, E. vilicicii, E. infula, E. adriatica, and E. umbratica. The most important morphological features for species delineation were cell shape, the degree and mode of torsion, valve apices, the appearance and structure of the transition between keel and valve body, the ultrastructure and the shape of the girdle bands, and the arrangement and density of perforations along the valve and valvocopulae. Our results highlight the underappreciated diversity of Entomoneis and call for a more in-depth morphological and molecular investigation of this genus especially in planktonic habitats.

Hoarding and horizontal transfer led to an expanded gene and intron repertoire in the plastid genome of the diatom, Toxarium undulatum (Bacillariophyta).

Although the plastid genomes of diatoms maintain a conserved architecture and core gene set, considerable variation about this core theme exists and can be traced to several different processes. Gene duplication, pseudogenization, and loss, as well as intracellular transfer of genes to the nuclear genome, have all contributed to variation in gene content among diatom species. In addition, some noncoding sequences have highly restricted phylogenetic distributions that suggest a recent foreign origin. We sequenced the plastid genome of the marine diatom, Toxarium undulatum, and found that the genome contains three genes (chlB, chlL, and chlN) involved in light-independent chlorophyll a biosynthesis that were not previously known from diatoms. Phylogenetic and syntenic data suggest that these genes were differentially retained in this one lineage as they were repeatedly lost from most other diatoms. Unique among diatoms and other heterokont algae sequenced so far, the genome also contains a large group II intron within an otherwise intact psaA gene. Although the intron is most similar to one in the plastid-encoded psaA gene of some green algae, high sequence divergence between the diatom and green algal introns rules out recent shared ancestry. We conclude that the psaA intron was likely introduced into the plastid genome of T. undulatum, or some earlier ancestor, by horizontal transfer from an unknown donor. This genome further highlights the myriad processes driving variation in gene and intron content in the plastid genomes of diatoms, one of the world's foremost primary producers.

Phylogeny, ecology, morphological evolution, and reclassification of the diatom orders Surirellales and Rhopalodiales.

The Surirellales and Rhopalodiales are large, widespread, and morphologically diverse groups of raphid pennate diatoms (Bacillariphyta) whose raphe, a structure that facilitates active motility, opens internally into a siliceous canal. We collected 202 representatives of the lineage and sequenced genes from the nuclear, plastid, and mitochondrial genomes to infer phylogenetic relationships as a basis for comparative study of ecology and morphological evolution as well as reclassification. The lineage was ancestrally marine, and we report the first evidence for a 'stepping stone' model of marine-freshwater transitions in which freshwater invasions were preceded by adaptation to intermediate brackish habitats. Phylogenetic comparative analyses also showed that the shift from an apical (e.g., Entomoneis) to transapical major axis of development (e.g., Surirella) did not have to proceed through subcircular intermediate forms (i.e., Campylodiscus). Rather, subcircular forms evolved both within lineages with longer apical axis or longer transapical axis. We also used the inferred phylogeny as a basis for genus-level reclassification of the lineage. Campylodiscus now includes the fastuosoid members of Surirella and Campylodiscus, but excludes other marine Campylodiscus which are now classified as Coronia. Surirella includes the Surirella striatula clade, Surirella Pinnatae group, and species formerly classified as Cymatopleura. We resurrected the genus Iconella to accommodate Stenopterobia and the robustoid members of Surirella and Campylodiscus. We broadened Epithemia to include members of the paraphyletic genus Rhopalodia. Finally, we discuss the challenges of constructing a classification that best leverages available phylogenetic data, while minimizing disruption to the research community and recognizing practical considerations stemming from the slow rate of progress on systematic studies of understudied organisms.

Dissecting signal and noise in diatom chloroplast protein encoding genes with phylogenetic information profiling.

Previous analyses of single diatom chloroplast protein-encoded genes recovered results highly incongruent with both traditional phylogenies and phylogenies derived from the nuclear encoded small subunit (SSU) gene. Our analysis here of six individual chloroplast genes (atpB, psaA, psaB, psbA, psbC and rbcL) obtained similar anomalous results. However, phylogenetic noise in these genes did not appear to be correlated, and their concatenation appeared to effectively sum their collective signal. We empirically demonstrated the value of combining phylogenetic information profiling, partitioned Bremer support and entropy analysis in examining the utility of various partitions in phylogenetic analysis. Noise was low in the 1st and 2nd codon positions, but so was signal. Conversely, high noise levels in the 3rd codon position was accompanied by high signal. Perhaps counterintuitively, simple exclusion experiments demonstrated this was especially true at deeper nodes where the 3rd codon position contributed most to a result congruent with morphology and SSU (and the total evidence tree here). Correlated with our empirical findings, probability of correct signal (derived from information profiling) increased and the statistical significance of substitutional saturation decreased as data were aggregated. In this regard, the aggregated 3rd codon position performed as well or better than more slowly evolving sites. Simply put, direct methods of noise removal (elimination of fast-evolving sites) disproportionately removed signal. Information profiling and partitioned Bremer support suggest that addition of chloroplast data will rapidly improve our understanding of the diatom phylogeny, but conversely also illustrate that some parts of the diatom tree are likely to remain recalcitrant to addition of molecular data. The methods based on information profiling have been criticized for their numerous assumptions and parameter estimates and the fact that they are based on quartets of taxa. Our empirical results support theoretical arguments that the simplifying assumptions made in these methods are robust to "real-life" situations.

Serial gene losses and foreign DNA underlie size and sequence variation in the plastid genomes of diatoms.

Photosynthesis by diatoms accounts for roughly one-fifth of global primary production, but despite this, relatively little is known about their plastid genomes. We report the completely sequenced plastid genomes for eight phylogenetically diverse diatoms and show them to be variable in size, gene and foreign sequence content, and gene order. The genomes contain a core set of 122 protein-coding genes, with 15 additional genes exhibiting complex patterns of 1) gene losses at varying phylogenetic scales, 2) functional transfers to the nucleus, 3) gene duplication, divergence, and differential retention of paralogs, and 4) acquisitions of putatively functional recombinase genes from resident plasmids. The newly sequenced genomes also contain several previously unreported genes, highlighting how poorly characterized diatom plastid genomes are overall. Genome size variation reflects major expansions of the inverted repeat region in some cases but, more commonly, large-scale expansions of intergenic regions, many of which contain unique open reading frames of likely foreign origin. Although many gene clusters are conserved across species, rearrangements appear to be frequent in most lineages.

Transcriptomic Insights into the Life History of Bolidophytes, the Sister Lineage to Diatoms.

Diatoms are perhaps the most diverse lineage of eukaryotic algae, with their siliceous cell wall and diplontic life history often considered to have played important roles in their extraordinary diversification. The characteristic diminution of the diatom cell wall over the course of vegetative growth provides a reliable, intrinsic trigger for sexual reproduction, establishing a direct link between the evolution of their cell-wall and life-history features. It is unclear, however, whether the diplontic life cycle of diatoms represents an ancestral or derived trait. This uncertainty is based in part on our lack of understanding of the life cycle of the sister lineage to diatoms, which includes a mix of two free-living and separately classified forms: naked biflagellate unicells in the genus Bolidomonas and silicified forms in the order Parmales. These two forms might represent different life-history stages, although directly establishing such links can be difficult. We sequenced transcriptomes for Bolidomonas and two diatoms and found that ~0.1% of the coding regions in the two diploid diatoms are heterozygous, whereas Bolidomonas is virtually devoid of heterozygous alleles, consistent with expectations for a haploid genome. These results suggest that Bolidomonas is haploid and predict that parmaleans represent the diploid phase of a haplodiplontic life cycle. These data fill an important gap in our understanding of the origin of the diplontic life history of diatoms, which may represent an evolutionarily derived, adaptive feature.

Origin and evolution of the canal raphe system in diatoms.

One lineage of pennate diatoms has a slit through the siliceous cell wall, called a "raphe," that functions in motility. Raphid pennate diatoms number in the perhaps tens of thousands of species, with the diversity of raphe forms potentially matching this number. Three lineages-the Bacillariales, Rhopalodiales, and Surirellales-possess a complex and presumably highly derived raphe that is physically separated from the cell interior, most often by a set of siliceous braces. Because the relationship among these three lineages is unclear, the number of origins of the canal raphe system and the homology of it and its constitutive parts among these lineages, is equally unclear. We reconstructed the phylogeny of raphid pennate diatoms and included, for the first time, members of all three canal raphid diatom lineages, and used the phylogeny to test specific hypotheses about the origin of the canal raphe. The canal raphe appears to have evolved twice, once in the common ancestor of Bacillariales and once in the common ancestor of Rhopalodiales and Surirellales, which form a monophyletic group in our analyses. These results recommend careful follow-up morphogenesis studies of the canal raphe in these two lineages to determine the underlying developmental basis for this remarkable case of parallel evolution.

Extensive reorganization of the plastid genome of Trifolium subterraneum (Fabaceae) is associated with numerous repeated sequences and novel DNA insertions.

The plastid genome of Trifolium subterraneum is 144,763 bp, about 20 kb longer than those of closely related legumes, which also lost one copy of the large inverted repeat (IR). The genome has undergone extensive genomic reconfiguration, including the loss of six genes (accD, infA, rpl22, rps16, rps18, and ycf1) and two introns (clpP and rps12) and numerous gene order changes, attributable to 14-18 inversions. All endpoints of rearranged gene clusters are flanked by repeated sequences, tRNAs, or pseudogenes. One unusual feature of the Trifolium subterraneum genome is the large number of dispersed repeats, which comprise 19.5% (ca. 28 kb) of the genome (versus about 4% for other angiosperms) and account for part of the increase in genome size. Nine genes (psbT, rbcL, clpP, rps3, rpl23, atpB, psbN, trnI-cau, and ycf3) have also been duplicated either partially or completely. rpl23 is the most highly duplicated gene, with portions of this gene duplicated six times. Comparisons of the Trifolium plastid genome with the Plant Repeat Database and searches for flanking inverted repeats suggest that the high incidence of dispersed repeats and rearrangements is not likely the result of transposition. Trifolium has 19.5 kb of unique DNA distributed among 160 fragments ranging in size from 30 to 494 bp, greatly surpassing the other five sequenced legume plastid genomes in novel DNA content. At least some of this unique DNA may represent horizontal transfer from bacterial genomes. These unusual features provide direction for the development of more complex models of plastid genome evolution.