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

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.

Revision of Ardissoneaceae (Bacillariophyta, Mediophyceae) from Micronesian populations, with descriptions of two new genera, Ardissoneopsis and Grunowago, and new species in Ardissonea, Synedrosphenia and Climacosphenia

Ardissonea was resurrected from Synedra in 1986 and was included as a genus by Round, Crawford and Mann ("The Diatoms") in its own Family and Order. They commented that there might be several genera involved since the type species of the genus possesses a double-walled structure and other taxa placed in Ardissonea have only a single-walled structure. Two other genera of "big sticks," Toxarium and Climacosphenia, were placed in their own Families and Orders but share many characters with Ardissoneaceae, especially growth from a bifacial annulus. Eighteen taxa (11 new species) from Micronesia were compared with the literature and remnant material from Grunow's Honduras Sargassum sample to address the concepts of Ardissonea and Ardissoneaceae. Phylogenetic and morphological analyses showed three clades within Ardissonea sensu lato: Ardissonea emend. for the double-walled taxa, Synedrosphenia emend. and Ardissoneopsis gen. nov. for single-walled taxa. New species include Ardissoneadensistriata sp. nov.; Synedrospheniabikarensis sp. nov., S.licmophoropsis sp. nov., S.parva sp. nov., and S.recta sp. nov.; Ardissoneopsisfulgicans sp. nov., A.appressata sp. nov., and A.gracilis sp. nov. Transfers include Synedrospheniacrystallina comb. nov. and S.fulgens comb. nov. Synedraundosa, seen for the first time in SEM in Grunow's material, is transferred to Ardissoneopsisundosa comb. nov. Three more genera have similar structure: Toxarium, Climacosphenia and Grunowago gen. nov., erected for Synedrabacillaris and a lanceolate species, G.pacifica sp. nov. Morphological characters of Toxarium in our region support separation of Toxariumhennedyanum and T.undulatum and suggest additional species here and elsewhere. Climacospheniamoniligera was not found but we clarify its characters based on the literature and distinguish C.soulonalis sp. nov. from it. Climacospheniaelongata and a very long, slender C.elegantissima sp. nov., previously identified as C.elongata, were present along with C.scimiter. Morphological and molecular phylogenetics strongly suggested that all these genera belong in one family and we propose to include them in the Ardissoneacae and to reinstate the Order Ardissoneales Round.

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.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

Diploid genomic architecture of Nitzschia inconspicua, an elite biomass production diatom

A near-complete diploid nuclear genome and accompanying circular mitochondrial and chloroplast genomes have been assembled from the elite commercial diatom species Nitzschia inconspicua. The 50 Mbp haploid size of the nuclear genome is nearly double that of model diatom Phaeodactylum tricornutum, but 30% smaller than closer relative Fragilariopsis cylindrus. Diploid assembly, which was facilitated by low levels of allelic heterozygosity (2.7%), included 14 candidate chromosome pairs composed of long, syntenic contigs, covering 93% of the total assembly. Telomeric ends were capped with an unusual 12-mer, G-rich, degenerate repeat sequence. Predicted proteins were highly enriched in strain-specific marker domains associated with cell-surface adhesion, biofilm formation, and raphe system gliding motility. Expanded species-specific families of carbonic anhydrases suggest potential enhancement of carbon concentration efficiency, and duplicated glycolysis and fatty acid synthesis pathways across cytosolic and organellar compartments may enhance peak metabolic output, contributing to competitive success over other organisms in mixed cultures. The N. inconspicua genome delivers a robust new reference for future functional and transcriptomic studies to illuminate the physiology of benthic pennate diatoms and harness their unique adaptations to support commercial algae biomass and bioproduct production.

Extreme Enlargement of the Inverted Repeat Region in the Plastid Genomes of Diatoms from the Genus Climaconeis

We sequenced the plastid genomes of three diatoms from the genus Climaconeis, including two strains formerly designated as Climaconeis scalaris. At 208,097 and 216,580 bp, the plastid genomes of the latter strains are the largest ever sequenced among diatoms and their increased size is explained by the massive expansion of the inverted repeat region. Important rearrangements of gene order were identified among the two populations of Climaconeis cf. scalaris. The other sequenced Climaconeis chloroplast genome is 1.5 times smaller compared with those of the Climaconeis cf. scalaris strains and it features an usual quadripartite structure. The extensive structural changes reported here for the genus Climaconeis are compared with those previously observed for other algae and plants displaying large plastid genomes.

Evidence of horizontal gene transfer between land plant plastids has surprising conservation implications

BACKGROUND AND AIMS

Horizontal gene transfer (HGT) is an important evolutionary mechanism because it transfers genetic material that may code for traits or functions between species or genomes. It is frequent in mitochondrial and nuclear genomes but has not been demonstrated between plastid genomes of different green land plant species.

METHODS

We Sanger-sequenced the nuclear internal transcribed spacers (ITS1 and 2) and the plastid rpl16 G2 intron (rpl16). In five individuals with foreign rpl16 we also sequenced atpB-rbcL and trnLUAA-trnFGAA.

KEY RESULTS

We discovered 14 individuals of a moss species with typical nuclear ITSs but foreign plastid rpl16 from a species of a distant lineage. None of the individuals with three plastid markers sequenced contained all foreign markers, demonstrating the transfer of plastid fragments rather than the entire plastid genome, i.e. entire plastids were not transferred. The two lineages diverged 165-185 Myr BP. The extended time interval since lineage divergence suggests that the foreign rpl16 is more likely explained by HGT than by hybridization or incomplete lineage sorting.

CONCLUSIONS

We provide the first conclusive evidence of interspecific plastid-to-plastid HGT among land plants. Two aspects are critical: it occurred at several localities during the massive colonization of recently disturbed open habitats that were created by large-scale liming as a freshwater biodiversity conservation measure; and it involved mosses whose unique life cycle includes spores that first develop a filamentous protonema phase. We hypothesize that gene transfer is facilitated when protonema filaments of different species intermix intimately when colonizing disturbed early succession habitats.

Haslea silbo, A Novel Cosmopolitan Species of Blue Diatoms

Specimens of a new species of blue diatoms from the genus Haslea Simonsen were discovered in geographically distant sampling sites, first in the Canary Archipelago, then North Carolina, Gulf of Naples, the Croatian South Adriatic Sea, and Turkish coast of the Eastern Mediterranean Sea. An exhaustive characterization of these specimens, using a combined morphological and genomic approach led to the conclusion that they belong to a single new to science cosmopolitan species, Haslea silbo sp. nov. A preliminary characterization of its blue pigment shows similarities to marennine produced by Haslea ostrearia, as evidenced by UV-visible spectrophotometry and Raman spectrometry. Life cycle stages including auxosporulation were also observed, providing data on the cardinal points of this species. For the two most geographically distant populations (North Carolina and East Mediterranean), complete mitochondrial and plastid genomes were sequenced. The mitogenomes of both strains share a rare atp6 pseudogene, but the number, nature, and positions of the group II introns inside its cox1 gene differ between the two populations. There are also two pairs of genes fused in single ORFs. The plastid genomes are characterized by large regions of recombination with plasmid DNA, which are in both cases located between the ycf35 and psbA genes, but whose content differs between the strains. The two sequenced strains hosts three plasmids coding for putative serine recombinase protein whose sequences are compared, and four out of six of these plasmids were highly conserved.

Seasonal and Geographical Transitions in Eukaryotic Phytoplankton Community Structure in the Atlantic and Pacific Oceans

Much is known about how broad eukaryotic phytoplankton groups vary according to nutrient availability in marine ecosystems. However, genus- and species-level dynamics are generally unknown, although important given that adaptation and acclimation processes differentiate at these levels. We examined phytoplankton communities across seasonal cycles in the North Atlantic (BATS) and under different trophic conditions in the eastern North Pacific (ENP), using phylogenetic classification of plastid-encoded 16S rRNA amplicon sequence variants (ASVs) and other methodologies, including flow cytometric cell sorting. Prasinophytes dominated eukaryotic phytoplankton amplicons during the nutrient-rich deep-mixing winter period at BATS. During stratification (‘summer’) uncultured dictyochophytes formed ∼35 ± 10% of all surface plastid amplicons and dominated those from stramenopile algae, whereas diatoms showed only minor, ephemeral contributions over the entire year. Uncultured dictyochophytes also comprised a major fraction of plastid amplicons in the oligotrophic ENP. Phylogenetic reconstructions of near-full length 16S rRNA sequences established 11 uncultured Dictyochophyte Environmental Clades (DEC). DEC-I and DEC-VI dominated surface dictyochophytes under stratification at BATS and in the ENP, and DEC-IV was also important in the latter. Additionally, although less common at BATS, Florenciella-related clades (FC) were prominent at depth in the ENP. In both ecosystems, pelagophytes contributed notably at depth, with PEC-VIII (Pelagophyte Environmental Clade) and (cultured) Pelagomonas calceolata being most important. Q-PCR confirmed the near absence of P. calceolata at the surface of the same oligotrophic sites where it reached ∼1,500 18S rRNA gene copies ml^–1 at the DCM. To further characterize phytoplankton present in our samples, we performed staining and at-sea single-cell sorting experiments. Sequencing results from these indicated several uncultured dictyochophyte clades are comprised of predatory mixotrophs. From an evolutionary perspective, these cells showed both conserved and unique features in the chloroplast genome. In ENP metatranscriptomes we observed high expression of multiple chloroplast genes as well as expression of a selfish element (group II intron) in the psaA gene. Comparative analyses across the Pacific and Atlantic sites support the conclusion that predatory dictyochophytes thrive under low nutrient conditions. The observations that several uncultured dictyochophyte lineages are seemingly capable of photosynthesis and predation, raises questions about potential shifts in phytoplankton trophic roles associated with seasonality and long-term ocean change.

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.

Metabolic Innovations Underpinning the Origin and Diversification of the Diatom Chloroplast

Of all the eukaryotic algal groups, diatoms make the most substantial contributions to photosynthesis in the contemporary ocean. Understanding the biological innovations that have occurred in the diatom chloroplast may provide us with explanations to the ecological success of this lineage and clues as to how best to exploit the biology of these organisms for biotechnology. In this paper, we use multi-species transcriptome datasets to compare chloroplast metabolism pathways in diatoms to other algal lineages. We identify possible diatom-specific innovations in chloroplast metabolism, including the completion of tocopherol synthesis via a chloroplast-targeted tocopherol cyclase, a complete chloroplast ornithine cycle, and chloroplast-targeted proteins involved in iron acquisition and CO₂ concentration not shared between diatoms and their closest relatives in the stramenopiles. We additionally present a detailed investigation of the chloroplast metabolism of the oil-producing diatom Fistulifera solaris, which is of industrial interest for biofuel production. These include modified amino acid and pyruvate hub metabolism that might enhance acetyl-coA production for chloroplast lipid biosynthesis and the presence of a chloroplast-localised squalene synthesis pathway unknown in other diatoms. Our data provides valuable insights into the biological adaptations underpinning an ecologically critical lineage, and how chloroplast metabolism can change even at a species level in extant algae.

Evolution of light-independent protochlorophyllide oxidoreductase

The nonhomologous enzymes, the light-independent protochlorophyllide reductase (DPOR) and the light-dependent protochlorophyllide oxidoreductase (LPOR), catalyze the reduction of protochlorophyllide (Pchlide) to chlorophyllide (Chlide) in the penultimate step of biosynthesis of chlorophyll (Chl) required for photosynthetic light absorption and energy conversion. The two enzymes differ with respect to the requirement of light for catalysis and oxygen sensitivity. DPOR and LPOR initially evolved in the ancestral prokaryotic genome perhaps at different times. DPOR originated in the anoxygenic environment of the Earth from nitrogenase-like enzyme of methanogenic archaea. Due to the transition from anoxygenic to oxygenic photosynthesis in the prokaryote, the DPOR was mostly inactivated in the daytime by photosynthetic O₂ leading to the evolution of oxygen-insensitive LPOR that could function in the light. The primary endosymbiotic event transferred the DPOR and LPOR genes to the eukaryotic phototroph; the DPOR remained in the genome of the ancestor that turned into the plastid, whereas LPOR was transferred to the host nuclear genome. From an evolutionary point of view, several compelling theories that explain the disappearance of DPOR from several species cutting across different phyla are as follows: (i) pressure of the oxygenic environment; (ii) change in the light conditions and temperature; and (iii) lineage-specific gene losses, RNA editing, and nonsynonymous substitution. Certain primary amino acid sequence and the physiochemical properties of the ChlL subunit of DPOR have similarity with that of LPOR suggesting a convergence of these two enzymes in certain evolutionary event. The newly obtained sequence data from different phototrophs will further enhance the width of the phylogenetic information on DPOR.

Plastome-Wide Rearrangements and Gene Losses in Carnivorous Droseraceae

The plastid genomes of four related carnivorous plants (Drosera regia, Drosera erythrorhiza, Aldrovanda vesiculosa, and Dionaea muscipula) were sequenced to examine changes potentially induced by the transition to carnivory. The plastid genomes of the Droseraceae show multiple rearrangements, gene losses, and large expansions or contractions of the inverted repeat. All the ndh genes are lost or nonfunctional, as well as in some of the species, clpP1, ycf1, ycf2 and some tRNA genes. Uniquely, among land plants, the trnK gene has no intron. Carnivory in the Droseraceae coincides with changes in plastid gene content similar to those induced by parasitism and mycoheterotrophy, suggesting parallel changes in chloroplast function due to the similar switch from autotrophy to (mixo-) heterotrophy. A molecular phylogeny of the taxa based on all shared plastid genes indicates that the "snap-traps" of Aldrovanda and Dionaea have a common origin.

Sequencing and Analysis of Chrysanthemum carinatum Schousb and Kalimeris indica. The Complete Chloroplast Genomes Reveal Two Inversions and rbcL as Barcoding of the Vegetable

Chrysanthemum carinatum Schousb and Kalimeris indica are widely distributed edible vegetables and the sources of the Chinese medicine Asteraceae. The complete chloroplast (cp) genome of Asteraceae usually occurs in the inversions of two regions. Hence, the cp genome sequences and structures of Asteraceae species are crucial for the cp genome genetic diversity and evolutionary studies. Hence, in this paper, we have sequenced and analyzed for the first time the cp genome size of C. carinatum Schousb and K. indica, which are 149,752 bp and 152,885 bp, with a pair of inverted repeats (IRs) (24,523 bp and 25,003) separated by a large single copy (LSC) region (82,290 bp and 84,610) and a small single copy (SSC) region (18,416 bp and 18,269), respectively. In total, 79 protein-coding genes, 30 distinct transfer RNA (tRNA) genes, four distinct rRNA genes and two pseudogenes were found not only in C. carinatum Schousb but also in the K. indica cp genome. Fifty-two (52) and fifty-nine (59) repeats, and seventy (70) and ninety (90) simple sequence repeats (SSRs) were found in the C. carinatum Schousb and K. indica cp genomes, respectively. Codon usage analysis showed that leucine, isoleucine, and serine are the most frequent amino acids and that the UAA stop codon was the significantly favorite stop codon in both cp genomes. The two inversions, the LSC region ranging from trnC-GCA to trnG-UCC and the whole SSC region were found in both of them. The complete cp genome comparison with other Asteraceae species showed that the coding area is more conservative than the non-coding area. The phylogenetic analysis revealed that the rbcL gene is a good barcoding marker for identifying different vegetables. These results give an insight into the identification, the barcoding, and the understanding of the evolutionary model of the Asteraceae cp genome.

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.