Sequencing and analysis of the complete organellar genomes of Parmales, a closely related group to Bacillariophyta (diatoms)

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Evolutionary dynamics of plastomes in coscinodiscophycean diatoms revealed by comparative genomics

To understand the evolution of coscinodiscophycean diatoms, plastome sequences of six coscinodiscophycean diatom species were constructed and analyzed in this study, doubling the number of constructed plastome sequences in Coscinodiscophyceae (radial centrics). The platome sizes varied substantially in Coscinodiscophyceae, ranging from 119.1 kb of Actinocyclus subtilis to 135.8 kb of Stephanopyxis turris. Plastomes in Paraliales and Stephanopyxales tended to be larger than those in Rhizosoleniales and Coscinodiacales, which were due to the expansion of the inverted repeats (IRs) and to the marked increase of the large single copy (LSC). Phylogenomic analysis indicated that Paralia and Stephanopyxis clustered tightly to form the Paraliales-Stephanopyxales complex, which was sister to the Rhizosoleniales-Coscinodiscales complex. The divergence time between Paraliales and Stephanopyxales was estimated at 85 MYA in the middle Upper Cretaceous, indicating that Paraliales and Stephanopyxales appeared later than Coscinodiacales and Rhizosoleniales according to their phylogenetic relationships. Frequent losses of housekeeping protein-coding genes (PCGs) were observed in these coscinodiscophycean plastomes, indicating that diatom plastomes showed an ongoing reduction in gene content during evolution. Two acpP genes (acpP1 and acpP2) detected in diatom plastomes were found to be originated from an early gene duplication event occurred in the common progenitor after diatom emergence, rather than multiple independent gene duplications occurring in different lineages of diatoms. The IRs in Stephanopyxis turris and Rhizosolenia fallax-imbricata exhibited a similar trend of large expansion to the small single copy (SSC) and slightly small contraction from the LSC, which eventually led to the conspicuous increase in IR size. Gene order was highly conserved in Coscinodiacales, while multiple rearrangements were observed in Rhizosoleniales and between Paraliales and Stephanopyxales. Our results greatly expanded the phylogenetic breadth in Coscinodiscophyceae and gained novel insights into the evolution of plastomes in diatoms.

Comparative analysis of Thalassionema chloroplast genomes revealed hidden biodiversity

The cosmopolitan Thalassionema species are often dominant components of the plankton diatom flora and sediment diatom assemblages in all but the Polar regions, making important ecological contribution to primary productivity. Historical studies concentrated on their indicative function for the marine environment based primarily on morphological features and essentially ignored their genomic information, hindering in-depth investigation on Thalassionema biodiversity. In this project, we constructed the complete chloroplast genomes (cpDNAs) of seven Thalassionema strains representing three different species, which were also the first cpDNAs constructed for any species in the order Thalassionematales that includes 35 reported species and varieties. The sizes of these Thalassionema cpDNAs, which showed typical quadripartite structures, varied from 124,127 bp to 140,121 bp. Comparative analysis revealed that Thalassionema cpDNAs possess conserved gene content inter-species and intra-species, along with several gene losses and transfers. Besides, their cpDNAs also have expanded inverted repeat regions (IRs) and preserve large intergenic spacers compared to other diatom cpDNAs. In addition, substantial genome rearrangements were discovered not only among different Thalassionema species but also among strains of a same species T. frauenfeldii, suggesting much higher diversity than previous reports. In addition to confirming the phylogenetic position of Thalassionema species, this study also estimated their emergence time at approximately 38 Mya. The availability of the Thalassionema species cpDNAs not only helps understand the Thalassionema species, but also facilitates phylogenetic analysis of diatoms.

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.

Spatiotemporal Variations in Antarctic Protistan Communities Highlight Phytoplankton Diversity and Seasonal Dominance by a Novel Cryptophyte Lineage

The Andvord fjord in the West Antarctic Peninsula (WAP) is known for its productivity and abundant megafauna. Nevertheless, seasonal patterns of the molecular diversity and abundance of protistan community members underpinning WAP productivity remain poorly resolved. We performed spring and fall expeditions pursuing protistan diversity, abundance of photosynthetic taxa, and the connection to changing conditions. 18S rRNA amplicon sequence variant (ASV) profiles revealed diverse predatory protists spanning multiple eukaryotic supergroups, alongside enigmatic heterotrophs like the Picozoa. Among photosynthetic protists, cryptophyte contributions were notable. Analysis of plastid-derived 16S rRNA ASVs supported 18S ASV results, including a dichotomy between cryptophytes and diatom contributions previously reported in other Antarctic regions. We demonstrate that stramenopile and cryptophyte community structures have distinct attributes. Photosynthetic stramenopiles exhibit high diversity, with the polar diatom Fragilariopsis cylindrus, unidentified Chaetoceros species, and others being prominent. Conversely, ASV analyses followed by environmental full-length rRNA gene sequencing, electron microscopy, and flow cytometry revealed that a novel alga dominates the cryptophytes. Phylogenetic analyses established that TPG clade VII, as named here, is evolutionarily distinct from cultivated cryptophyte lineages. Additionally, cryptophyte cell abundance correlated with increased water temperature. Analyses of global data sets showed that clade VII dominates cryptophyte ASVs at Southern Ocean sites and appears to be endemic, whereas in the Arctic and elsewhere, Teleaulax amphioxeia and Plagioselmis prolonga dominate, although both were undetected in Antarctic waters. Collectively, our studies provide baseline data against which future change can be assessed, identify different diversification patterns between stramenopiles and cryptophytes, and highlight an evolutionarily distinct cryptophyte clade that thrives under conditions enhanced by warming.

Olisthodiscus represents a new class of Ochrophyta

The phylogenetic diversity of Ochrophyta, a diverse and ecologically important radiation of algae, is still incompletely understood even at the level of the principal lineages. One taxon that has eluded simple classification is the marine flagellate genus Olisthodiscus. We investigated Olisthodiscus luteus K-0444 and documented its morphological and genetic differences from the NIES-15 strain, which we described as Olisthodiscus tomasii sp. nov. Phylogenetic analyses of combined 18S and 28S rRNA sequences confirmed that Olisthodiscus constitutes a separate, deep, ochrophyte lineage, but its position could not be resolved. To overcome this problem, we sequenced the plastid genome of O. luteus K-0444 and used the new data in multigene phylogenetic analyses, which suggested that Olisthodiscus is a sister lineage of the class Pinguiophyceae within a broader clade additionally including Chrysophyceae, Synchromophyceae, and Eustigmatophyceae. Surprisingly, the Olisthodiscus plastid genome contained three genes, ycf80, cysT, and cysW, inherited from the rhodophyte ancestor of the ochrophyte plastid yet lost from all other ochrophyte groups studied so far. Combined with nuclear genes for CysA and Sbp proteins, Olisthodiscus is the only known ochrophyte possessing a plastidial sulfate transporter SulT. In addition, the finding of a cemA gene in the Olisthodiscus plastid genome and an updated phylogenetic analysis ruled out the previously proposed hypothesis invoking horizontal cemA transfer from a green algal plastid into Synurales. Altogether, Olisthodiscus clearly represents a novel phylogenetically distinct ochrophyte lineage, which we have proposed as a new class, Olisthodiscophyceae.

Complete mitochondrial genome of Rhizosolenia setigera (Coscinodiscophyceae, Bacillariophyta)

Rhizosolenia is a species-rich genus with 144 described species, many of which are harmful algal species (HABs) with significant negative ecological impact. Despite their significance in primary production and their potential to induce HABs, genome data of these species remain extremely limited. In this study, the complete mitochondrial genome (mtDNA) of Rhizosolenia setigera Brightwell 1858 was determined for the first time, which also represented the first mtDNA of the order Rhizosoleniales. The circular mtDNA was 34,792 bp in length with GC content of 23.28%. It encoded 63 genes including 35 protein-coding genes (PCGs), 24 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 2 conserved open reading frames (orfs). Phylogenetic analysis using concatenated PCGs revealed that R. setigera and Melosira undulate, which also belongs to the class Coscinodiscophyceae, clustered together as expected. However, comparison of these two mtDNAs revealed extensive genome rearrangement events, suggesting large evolutionary distance. The complete mtDNA of R. setigera will facilitate research on the phylogenetic relationship among Rhizosolenia species, which will in turn facilitate exploration of the evolutionary relationships in the class of Coscinodiscophyceae.

Draft Whole-Genome Sequence of Triparma laevis f. inornata (Parmales, Bolidophyceae), Isolated from the Oyashio Region, Western North Pacific Ocean

We present the first draft whole-genome sequence for the Parmales (Bolidophyceae, Heterokonta), a picoplanktonic sister group of diatoms, using a Triparma laevis f. inornata strain that was isolated from the Oyashio region in the western North Pacific Ocean.

We present the first draft whole-genome sequence for the Parmales (Bolidophyceae, Heterokonta), a picoplanktonic sister group of diatoms, using a Triparma laevis f. inornata strain that was isolated from the Oyashio region in the western North Pacific Ocean.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

In this study, the mitogenomes of two Paxillus species were assembled, annotated and compared. The two mitogenomes of Paxillus involutus and P. rubicundulus comprised circular DNA molecules, with the size of 39,109 bp and 41,061 bp, respectively. Evolutionary analysis revealed that the nad4L gene had undergone strong positive selection in the two Paxillus species. In addition, 10.64 and 36.50% of the repetitive sequences were detected in the mitogenomes of P. involutus and P. rubicundulus, respectively, which might transfer between mitochondrial and nuclear genomes. Large-scale gene rearrangements and frequent intron gain/loss events were detected in 61 basidiomycete species, which revealed large variations in mitochondrial organization and size in Basidiomycota. In addition, the insertion sites of the basidiomycete introns were found to have a base preference. Phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies, indicating that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships of Basidiomycota. This study is the first report on the mitogenomes of Paxillus, which will promote a better understanding of their contrasted ecological strategies, molecular evolution and phylogeny of these important ectomycorrhizal fungi and related basidiomycete species.

Complex origins of chloroplast membranes with photosynthetic machineries: multiple transfers of genes from divergent organisms at different times or a single endosymbiotic event?

The paradigm "cyanobacterial origin of chloroplasts" is currently viewed as an established fact. However, we may have to re-consider the origin of chloroplast membranes, because membranes are not replicated by their own. It is the genes for lipid biosynthetic enzymes that are inherited. In the current understandings, these enzymes became encoded by the nuclear genome as a result of endosymbiotic gene transfer from the endosymbiont. However, we previously showed that many enzymes involved in the synthesis of chloroplast peptidoglycan and glycolipids did not originate from cyanobacteria. Here I present results of comprehensive phylogenetic analysis of chloroplast enzymes involved in fatty acid and lipid biosynthesis, as well as additional chloroplast components related to photosynthesis and gene expression. Four types of phylogenetic relationship between chloroplast enzymes (encoded by the chloroplast and nuclear genomes) and cyanobacterial counterparts were found: type 1, chloroplast enzymes diverged from inside of cyanobacterial clade; type 2, chloroplast and cyanobacterial enzymes are sister groups; type 3, chloroplast enzymes originated from homologs of bacteria other than cyanobacteria; type 4, chloroplast enzymes diverged from eukaryotic homologs. Estimation of evolutionary distances suggested that the acquisition times of chloroplast enzymes were diverse, indicating that multiple gene transfers accounted for the chloroplast enzymes analyzed. Based on the results, I try to relax the tight logic of the endosymbiotic origin of chloroplasts involving a single endosymbiotic event by proposing alternative hypotheses. The hypothesis of host-directed chloroplast formation proposes that glycolipid synthesis ability had been acquired by the eukaryotic host before the acquisition of chloroplast ribosomes. Chloroplast membrane system could have been provided by the host, whereas cyanobacteria contributed to the genes for the genetic and photosynthesis systems, at various times, either before or after the formation of chloroplast membranes. The origin(s) of chloroplasts seems to be more complicated than the single event of primary endosymbiosis.

Dictyochophyceae Plastid Genomes Reveal Unusual Variability in Their Organization

Dictyochophyceae (silicoflagellates) are unicellular freshwater and marine algae (Heterokontophyta, stramenopiles). Despite their abundance in global oceans and potential ecological significance, discovered in recent years, neither nuclear nor organellar genomes of representatives of this group were sequenced until now. Here, we present the first complete plastid genome sequences of Dictyochophyceae, obtained from four species: Dictyocha speculum, Rhizochromulina marina, Florenciella parvula and Pseudopedinella elastica. Despite their comparable size and genetic content, these four plastid genomes exhibit variability in their organization: plastid genomes of F. parvula and P. elastica possess conventional quadripartite structure with a pair of inverted repeats, R. marina instead possesses two direct repeats with the same orientation and D. speculum possesses no repeats at all. We also observed a number of unusual traits in the plastid genome of D. speculum, including expansion of the intergenic regions, presence of an intron in the otherwise non-intron-bearing psaA gene, and an additional copy of the large subunit of RuBisCO gene (rbcL), the last of which has never been observed in any plastid genome. We conclude that despite noticeable gene content similarities between the plastid genomes of Dictyochophyceae and their relatives (pelagophytes, diatoms), the number of distinctive features observed in this lineage strongly suggests that additional taxa require further investigation.

Plastid Genomes and Proteins Illuminate the Evolution of Eustigmatophyte Algae and Their Bacterial Endosymbionts

Eustigmatophytes, a class of stramenopile algae (ochrophytes), include not only the extensively studied biotechnologically important genus Nannochloropsis but also a rapidly expanding diversity of lineages with much less well characterized biology. Recent discoveries have led to exciting additions to our knowledge about eustigmatophytes. Some proved to harbor bacterial endosymbionts representing a novel genus, Candidatus Phycorickettsia, and an operon of unclear function (ebo) obtained by horizontal gene transfer from the endosymbiont lineage was found in the plastid genomes of still other eustigmatophytes. To shed more light on the latter event, as well as to generally improve our understanding of the eustigmatophyte evolutionary history, we sequenced plastid genomes of seven phylogenetically diverse representatives (including new isolates representing undescribed taxa). A phylogenomic analysis of plastid genome-encoded proteins resolved the phylogenetic relationships among the main eustigmatophyte lineages and provided a framework for the interpretation of plastid gene gains and losses in the group. The ebo operon gain was inferred to have probably occurred within the order Eustigmatales, after the divergence of the two basalmost lineages (a newly discovered hitherto undescribed strain and the Pseudellipsoidion group). When looking for nuclear genes potentially compensating for plastid gene losses, we noticed a gene for a plastid-targeted acyl carrier protein that was apparently acquired by horizontal gene transfer from Phycorickettsia. The presence of this gene in all eustigmatophytes studied, including representatives of both principal clades (Eustigmatales and Goniochloridales), is a genetic footprint indicating that the eustigmatophyte-Phycorickettsia partnership started no later than in the last eustigmatophyte common ancestor.

Ontogenetic analysis of siliceous cell wall formation in Triparma laevis f. inornata (Parmales, Stramenopiles)

Triparma laevis f. inornata is a unicellular alga belonging to the Bolidophyceae, which is most closely related to diatoms. Like diatoms, T. laevis f. inornata has a siliceous cell wall. The cell wall of T. laevis f. inornata consists of four round plates (three shields and one ventral plate) and one dorsal and three girdle plates. But, unlike diatoms, T. laevis f. inornata cells can grow when concentrations of silica are depleted. We took advantage of this ability, using TEM to study the ontogeny of the siliceous plate, pattern center formation, and development. Two types of pattern centers (annulus and sternum) were observed in the early and middle stage of plate formation. During their formation, the annuli were initially crescent-shaped but eventually their ends fused to make a ring. Only outward silica deposition of the branching ribs occurred on the growing annulus until it became a ring, resulting in an unfilled circle inside the annulus. The pattern center of the shield plate was always an annulus, but in ventral plates both annulus and sternum were observed. The annuli and sterna in T. laevis f. inornata round plates were very similar to the annuli and sterna in diatom valves. These results suggested that the round plates of Parmales are homologous to diatom valves. This information on the plate ontogeny of T. laevis f. inornata provides new insights into the evolution of the siliceous cell wall in the Parmales and diatoms.

Comparative plastid genomics of Synurophyceae: inverted repeat dynamics and gene content variation

BACKGROUND

The Synurophyceae is one of most important photosynthetic stramenopile algal lineages in freshwater ecosystems. They are characterized by siliceous scales covering the cell or colony surface and possess plastids of red-algal secondary or tertiary endosymbiotic origin. Despite their ecological and evolutionary significance, the relationships amongst extant Synurophyceae are unclear, as is their relationship to most other stramenopiles.

RESULTS

Here we report a comparative analysis of plastid genomes sequenced from five representative synurophycean algae. Most of these plastid genomes are highly conserved with respect to genome structure and coding capacity, with the exception of gene re-arrangements and partial duplications at the boundary of the inverted repeat and single-copy regions. Several lineage-specific gene loss/gain events and intron insertions were detected (e.g., cemA, dnaB, syfB, and trnL).

CONCLUSIONS

Unexpectedly, the cemA gene of Synurophyceae shows a strong relationship with sequences from members of the green-algal lineage, suggesting the occurrence of a lateral gene transfer event. Using a molecular clock approach based on silica fossil record data, we infer the timing of genome re-arrangement and gene gain/loss events in the plastid genomes of Synurophyceae.

Plastid Genome of Dictyopteris divaricata (Dictyotales, Phaeophyceae): Understanding the Evolution of Plastid Genomes in Brown Algae

Dictyotophycidae is a subclass of brown algae containing 395 species that are distributed worldwide. A complete plastid (chloroplast) genome (ptDNA or cpDNA) had not previously been sequenced from this group. In this study, the complete plastid genome of Dictyopteris divaricata (Okamura) Okamura (Dictyotales, Phaeophyceae) was characterized and compared to other brown algal ptDNAs. This plastid genome was 126,099 bp in size with two inverted repeats (IRs) of 6026 bp. The D. divaricata IRs contained rpl21, making its IRs larger than representatives from the orders Fucales and Laminariales, but was smaller than that from Ectocarpales. The G + C content of D. divaricata (31.19%) was the highest of the known ptDNAs of brown algae (28.94-31.05%). Two protein-coding genes, rbcR and rpl32, were present in ptDNAs of Laminariales, Ectocarpales (Ectocarpus siliculosus), and Fucales (LEF) but were absent in D. divaricata. Reduced intergenic space (13.11%) and eight pairs of overlapping genes in D. divaricata ptDNA made it the most compact plastid genome in brown algae so far. The architecture of D. divaricata ptDNA showed higher similarity to that of Laminariales compared with Fucales and Ectocarpales. The difference in general features, gene content, and architecture among the ptDNAs of D. divaricata and LEF clade revealed the diversity and evolutionary trends of plastid genomes in brown algae.

The plastid genome of some eustigmatophyte algae harbours a bacteria-derived six-gene cluster for biosynthesis of a novel secondary metabolite

Acquisition of genes by plastid genomes (plastomes) via horizontal gene transfer (HGT) seems to be a rare phenomenon. Here, we report an interesting case of HGT revealed by sequencing the plastomes of the eustigmatophyte algae Monodopsis sp. MarTras21 and Vischeria sp. CAUP Q 202. These plastomes proved to harbour a unique cluster of six genes, most probably acquired from a bacterium of the phylum Bacteroidetes, with homologues in various bacteria, typically organized in a conserved uncharacterized putative operon. Sequence analyses of the six proteins encoded by the operon yielded the following annotation for them: (i) a novel family without discernible homologues; (ii) a new family within the superfamily of metallo-dependent hydrolases; (iii) a novel subgroup of the UbiA superfamily of prenyl transferases; (iv) a new clade within the sugar phosphate cyclase superfamily; (v) a new family within the xylose isomerase-like superfamily; and (vi) a hydrolase for a phosphate moiety-containing substrate. We suggest that the operon encodes enzymes of a pathway synthesizing an isoprenoid–cyclitol-derived compound, possibly an antimicrobial or other protective substance. To the best of our knowledge, this is the first report of an expansion of the metabolic capacity of a plastid mediated by HGT into the plastid genome.

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.