Transcript analysis of dinoflagellate plastid gene minicircles

Chloroplast mRNAs are 3' polyuridylylated in the Green Alga Pithophora roettleri (Cladophorales)

Species within the green algal order Cladophorales have an unconventional plastome structure where individual coding regions or small numbers of genes occur as linear single-stranded DNAs folded into hairpin structures. Another group of photosynthetic organisms with an equivalently reduced chloroplast genome are the peridinin dinoflagellates of the Alveolata eukaryotic lineage whose plastomes are mini-circles carrying one or a few genes required for photosynthesis. One unusual aspect of the Alveolata is the polyuridylylation of mRNA 3' ends among peridinin dinoflagellates and the chromerid algae. This study was conducted to understand if an unconventional highly reduced plastome structure co-occurs with unconventional RNA processing. To address this, the 5' and 3' mRNA termini of the known chloroplast genes of Pithophora roettleri (order Cladophorales) were analyzed for evidence of post-transcriptional processing. Circular Reverse Transcriptase PCR (cRT-PCR) followed by deep sequencing of the amplicons was used to analyze 5' and 3' mRNA termini. Evidence of several processing events were collected, most notably the 3' termini of six of the eight genes were polyuridylylated, which has not been reported for any lineage outside of the Alveolata. Other processing events include poly(A) and heteropolymeric 3' additions, 5' primary transcript start sites, as well as the presence of circularized RNAs. Five other species representing other green algal lineages were also tested and poly(U) additions appear to be limited to the order Cladophorales. These results demonstrate that chloroplast mRNA polyuridylylation is not the sole provenance of photosynthetic alveolates and may have convergently evolved in two distinct photosynthetic lineages.

PhyloMagnet: fast and accurate screening of short-read meta-omics data using gene-centric phylogenetics

Motivation

Metagenomic and metatranscriptomic sequencing have become increasingly popular tools for producing massive amounts of short-read data, often used for the reconstruction of draft genomes or the detection of (active) genes in microbial communities. Unfortunately, sequence assemblies of such datasets generally remain a computationally challenging task. Frequently, researchers are only interested in a specific group of organisms or genes; yet, the assembly of multiple datasets only to identify candidate sequences for a specific question is sometimes prohibitively slow, forcing researchers to select a subset of available datasets to address their question. Here, we present PhyloMagnet, a workflow to screen meta-omics datasets for taxa and genes of interest using gene-centric assembly and phylogenetic placement of sequences.

Results

Using PhyloMagnet, we could identify up to 87% of the genera in an in vitro mock community with variable abundances, while the false positive predictions per single gene tree ranged from 0 to 23%. When applied to a group of metagenomes for which a set of metagenome assembled genomes (MAGs) have been published, we could detect the majority of the taxonomic labels that the MAGs had been annotated with. In a metatranscriptomic setting, the phylogenetic placement of assembled contigs corresponds to that of transcripts obtained from transcriptome assembly.

Availability and implementation

PhyloMagnet is built using Nextflow, available at github.com/maxemil/PhyloMagnet and is developed and tested on Linux. It is released under the open source GNU GPL licence and documentation is available at phylomagnet.readthedocs.io. Version 0.5 of PhyloMagnet was used for all benchmarking experiments.

Supplementary information

Supplementary data are available at Bioinformatics online.

Integrated Genomic and Transcriptomic Analysis of the Peridinin Dinoflagellate Amphidinium carterae Plastid

The plastid genomes of peridinin-containing dinoflagellates are highly unusual, possessing very few genes, which are located on small chromosomal elements termed "minicircles". These minicircles may contain genes, or no recognisable coding information. Transcripts produced from minicircles may undergo unusual processing events, such as the addition of a 3' poly(U) tail. To date, little is known about the genetic or transcriptional diversity of non-coding sequences in peridinin dinoflagellate plastids. These sequences include empty minicircles, and regions of non-coding DNA in coding minicircles. Here, we present an integrated plastid genome and transcriptome for the model peridinin dinoflagellate Amphidinium carterae, identifying a previously undescribed minicircle. We also profile transcripts covering non-coding regions of the psbA and petB/atpA minicircles. We present evidence that antisense transcripts are produced within the A. carterae plastid, but show that these transcripts undergo different end cleavage events from sense transcripts, and do not receive 3' poly(U) tails. The difference in processing events between sense and antisense transcripts may enable the removal of non-coding transcripts from peridinin dinoflagellate plastid transcript pools.

Transcription of the apicoplast genome

Many members of the Apicomplexa contain a remnant chloroplast, known as an apicoplast. The apicoplast encodes numerous genes, and loss of the organelle is lethal. Here, we present a summary of what is known about apicoplast transcription. Unlike plant chloroplasts, there is a single RNA polymerase, and initial transcription is polycistronic. RNA is then cleaved into tRNA, mRNA and rRNA molecules. Significant levels of antisense transcription have been reported, together with a single case of RNA editing. Polycistronic transcription is also observed in the related algae Chromera and Vitrella, which retain a photosynthetic chloroplast. Surprisingly, a polyU tail is added to Chromera and Vitrella transcripts which encode proteins involved in photosynthesis. No such tail is added to Plasmodium transcripts. Transcription in the Apicomplexa is remarkably similar to that seen in the chloroplast of the related peridinin dinoflagellate algae, reflecting the common evolutionary origins of the organelle.

Transcriptome Analysis of Scrippsiella trochoidea CCMP 3099 Reveals Physiological Changes Related to Nitrate Depletion

Dinoflagellates are a major component of marine phytoplankton and many species are recognized for their ability to produce harmful algal blooms (HABs). Scrippsiella trochoidea is a non-toxic, marine dinoflagellate that can be found in both cold and tropic waters where it is known to produce “red tide” events. Little is known about the genomic makeup of S. trochoidea and a transcriptome study was conducted to shed light on the biochemical and physiological adaptations related to nutrient depletion. Cultures were grown under N and P limiting conditions and transcriptomes were generated via RNAseq technology. De novo assembly reconstructed 107,415 putative transcripts of which only 41% could be annotated. No significant transcriptomic response was observed in response to initial P depletion, however, a strong transcriptional response to N depletion was detected. Among the down-regulated pathways were those for glutamine/glutamate metabolism as well as urea and nitrate/nitrite transporters. Transcripts for ammonia transporters displayed both up- and down-regulation, perhaps related to a shift to higher affinity transporters. Genes for the utilization of DON compounds were up-regulated. These included transcripts for amino acids transporters, polyamine oxidase, and extracellular proteinase and peptidases. N depletion also triggered down regulation of transcripts related to the production of Photosystems I & II and related proteins. These data are consistent with a metabolic strategy that conserves N while maximizing sustained metabolism by emphasizing the relative contribution of organic N sources. Surprisingly, the transcriptome also contained transcripts potentially related to secondary metabolite production, including a homolog to the Short Isoform Saxitoxin gene (sxtA) from Alexandrium fundyense, which was significantly up-regulated under N-depletion. A total of 113 unique hits to Sxt genes, covering 17 of the 34 genes found in C. raciborskii were detected, indicating that S. trochoidea has previously unrecognized potential for the production of secondary metabolites with potential toxicity.

Diversity of transcripts and transcript processing forms in plastids of the dinoflagellate alga Karenia mikimotoi

Plastids produce a vast diversity of transcripts. These include mature transcripts containing coding sequences, and their processing precursors, as well as transcripts that lack direct coding functions, such as antisense transcripts. Although plastid transcriptomes have been characterised for many plant species, less is known about the transcripts produced in other plastid lineages. We characterised the transcripts produced in the fucoxanthin-containing plastids of the dinoflagellate alga Karenia mikimotoi. This plastid lineage, acquired through tertiary endosymbiosis, utilises transcript processing pathways that are very different from those found in plants and green algae, including 3' poly(U) tail addition, and extensive substitutional editing of transcript sequences. We have sequenced the plastid transcriptome of K. mikimotoi, and have detected evidence for divergent evolution of fucoxanthin plastid genomes. We have additionally characterised polycistronic and monocistronic transcripts from two plastid loci, psbD-tRNA (Met)-ycf4 and rpl36-rps13-rps11. We find evidence for a range of transcripts produced from each locus that differ in terms of editing state, 5' end cleavage position, and poly(U) tail addition. Finally, we identify antisense transcripts in K. mikimotoi, which appear to undergo different processing events from the corresponding sense transcripts. Overall, our study provides insights into the diversity of transcripts and processing intermediates found in plastid lineages across the eukaryotes.

Evolution of chloroplast transcript processing in Plasmodium and its chromerid algal relatives

It is well understood that apicomplexan parasites, such as the malaria pathogen Plasmodium, are descended from free-living algae, and maintain a vestigial chloroplast that has secondarily lost all genes of photosynthetic function. Recently, two fully photosynthetic relatives of parasitic apicomplexans have been identified, the ‘chromerid’ algae Chromera velia and Vitrella brassicaformis, which retain photosynthesis genes within their chloroplasts. Elucidating the processes governing gene expression in chromerid chloroplasts might provide valuable insights into the origins of parasitism in the apicomplexans. We have characterised chloroplast transcript processing pathways in C. velia, V. brassicaformis and P. falciparum with a focus on the addition of an unusual, 3′ poly(U) tail. We demonstrate that poly(U) tails in chromerids are preferentially added to transcripts that encode proteins that are directly involved in photosynthetic electron transfer, over transcripts for proteins that are not involved in photosynthesis. To our knowledge, this represents the first chloroplast transcript processing pathway to be associated with a particular functional category of genes. In contrast, Plasmodium chloroplast transcripts are not polyuridylylated. We additionally present evidence that poly(U) tail addition in chromerids is involved in the alternative processing of polycistronic precursors covering multiple photosynthesis genes, and appears to be associated with high levels of transcript abundance. We propose that changes to the chloroplast transcript processing machinery were an important step in the loss of photosynthesis in ancestors of parasitic apicomplexans.

Chloroplasts contain their own genomes, containing two broad functional types of gene: genes encoding proteins directly involved in photosynthesis, and genes with a non-photosynthesis function, such as cofactor biosynthesis, assembly of protein complexes, or expression of the chloroplast genome. Thus far, to our knowledge, no chloroplast gene expression pathways in any lineage have been found to target one functional category of gene specifically. Here, we show that a chloroplast RNA processing pathway – the addition of a 3′ poly(U) tail – is specifically associated with photosynthesis genes in two species of algae, the ‘chromerids’ Chromera and Vitrella. The addition of the poly(U) tail enables the precise processing of mature photosynthesis gene transcripts from precursor RNA, and is likely to be essential for expression of the chromerid photosynthesis machinery. The chromerid algae are the closest photosynthetic relatives of a parasitic group of eukaryotes, the apicomplexans, which include the malaria pathogen Plasmodium. Apicomplexans are descended from algae, and retain a reduced chloroplast, which contains genes only of non-photosynthesis function. We have confirmed that 3′ poly(U) tails are not added to Plasmodium chloroplast transcripts. The expression pathways associated with photosynthesis genes have therefore been lost in the evolution of the apicomplexan chloroplast, and this loss could potentially have driven the transition from photosynthesis to parasitism.

Transcript patterns of chloroplast-encoded genes in cultured Symbiodinium spp. (Dinophyceae): testing the influence of a light shift and diel periodicity

Microalgae possess numerous cellular mechanisms specifically employed for acclimating the photosynthetic pathways to changes in the physical environment. Despite the importance of coral-dinoflagellate symbioses, little focus has been given as to how the symbiotic algae (Symbiodinium spp.) regulate the expression of their photosynthetic genes. This study used real-time PCR to investigate the transcript abundance of the plastid-encoded genes, psbA (encoding the D1 protein of photosystem II) and psaA (encoding the P700 protein in photosystem I), within the cultured Symbiodinium ITS-2 (internal transcribed spacer region) types A20 and A13. Transcript abundance was monitored during a low to high-light shift, as well as over a full diel light cycle. In addition, psaA was characterized in three isolates (A20, A13, and D4-5) and noted as another example of a dinoflagellate plastid gene encoded on a minicircle. In general, the overall incongruence of transcript patterns for both psbA and psaA between the Symbiodinium isolates and other models of transcriptionally controlled chloroplast gene expression (e.g., Pisum sativum [pea], Sinapis alba [mustard seedling], and Synechocystis sp. PCC 6803 [cyanobacteria]) suggests that Symbiodinium is reliant on posttranscriptional mechanisms for homeostatic regulation of its photosynthetic proteins.

Polyuridylylation and processing of transcripts from multiple gene minicircles in chloroplasts of the dinoflagellate Amphidinium carterae

Although transcription and transcript processing in the chloroplasts of plants have been extensively characterised, the RNA metabolism of other chloroplast lineages across the eukaryotes remains poorly understood. In this paper, we use RT-PCR to study transcription and transcript processing in the chloroplasts of Amphidinium carterae, a model peridinin-containing dinoflagellate. These organisms have a highly unusual chloroplast genome, with genes located on multiple small 'minicircle' elements, and a number of idiosyncratic features of RNA metabolism including transcription via a rolling circle mechanism, and 3' terminal polyuridylylation of transcripts. We demonstrate that transcription occurs in A. carterae via a rolling circle mechanism, as previously shown in the dinoflagellate Heterocapsa, and present evidence for the production of both polycistronic and monocistronic transcripts from A. carterae minicircles, including several regions containing ORFs previously not known to be expressed. We demonstrate the presence of both polyuridylylated and non-polyuridylylated transcripts in A. carterae, and show that polycistronic transcripts can be terminally polyuridylylated. We present a model for RNA metabolism in dinoflagellate chloroplasts where long polycistronic precursors are processed to form mature transcripts. Terminal polyuridylylation may mark transcripts with the correct 3' end.

Organization and expression of organellar genomes

Protist mitochondrial genomes show a very wide range of gene content, ranging from three genes for respiratory chain components in Apicomplexa and dinoflagellates to nearly 100 genes in Reclinomonas americana. In many organisms the rRNA genes are fragmented, although still functional. Some protist mitochondria encode a full set of tRNAs, while others rely on imported molecules. There is similarly a wide variation in mitochondrial genome organization, even among closely related groups. Mitochondrial gene expression and control are generally poorly characterized. Transcription probably relies on a 'viral-type' RNA polymerase, although a 'bacterial-type' enzyme may be involved in some cases. Transcripts are heavily edited in many lineages. The chloroplast genome generally shows less variation in gene content and organization, although greatly reduced genomes are found in dinoflagellate algae and non-photosynthetic organisms. Genes in the former are located on small plasmids in contrast to the larger molecules found elsewhere. Control of gene expression in chloroplasts involves transcriptional and post-transcriptional regulation. Redox poise and the ATP/ADP ratio are likely to be important determinants. Some protists have an additional extranuclear genome, the nucleomorph, which is a remnant nucleus. Nucleomorphs of two separate lineages have a number of features in common.

Biosynthesis and molecular genetics of polyketides in marine dinoflagellates

Marine dinoflagellates are the single most important group of algae that produce toxins, which have a global impact on human activities. The toxins are chemically diverse, and include macrolides, cyclic polyethers, spirolides and purine alkaloids. Whereas there is a multitude of studies describing the pharmacology of these toxins, there is limited or no knowledge regarding the biochemistry and molecular genetics involved in their biosynthesis. Recently, however, exciting advances have been made. Expressed sequence tag sequencing studies have revealed important insights into the transcriptomes of dinoflagellates, whereas other studies have implicated polyketide synthase genes in the biosynthesis of cyclic polyether toxins, and the molecular genetic basis for the biosynthesis of paralytic shellfish toxins has been elucidated in cyanobacteria. This review summarises the recent progress that has been made regarding the unusual genomes of dinoflagellates, the biosynthesis and molecular genetics of dinoflagellate toxins. In addition, the evolution of these metabolic pathways will be discussed, and an outlook for future research and possible applications is provided.

Identification of transcribed and persistent variants of the psbA gene carried by plastid minicircles in a dinoflagellate

Plastid genomes of the peridinin-containing dinoflagellates are composed of a limited number of genes that are contained separately on small circular molecules (minicircles). It has been shown that occasionally aberrant minicircles are generated, but they are usually small and contain little coding information. In this study, we discovered multiple variants, a putative "gene family", of the gene psbA in the plastid minicircles of the dinoflagellate Alexandrium tamarense, which have persisted for almost 3 years in culture. Each variant, like the ordinary psbA, existed on distinct minicircles of similar size (5-6 kb). These psbA variants retained all, or almost all, the coding sequence of the ordinary gene, and all four were transcribed and edited after transcription, even though they could not encode the entire protein due to intervening or translocated sequences. Repeat elements were generally found in the relatively large non-coding region of these minicircles. Each psbA variant might have been generated by DNA recombination and/or replication slippage, as for previously reported aberrant minicircles. The fact that these minicircles are transcribed, individually edited and maintained in the genome suggests that they are functionally important, although their precise roles remain unclear.

Long transcripts from dinoflagellate chloroplast minicircles suggest "rolling circle" transcription

The chloroplast genome of a dinoflagellate consists of a group of small circular DNA molecules (minicircles), most of which carry a single gene. With RT-PCR, primer extension, and Northern analyses, we show that the entire minicircle is transcribed and that some minicircles can produce RNAs larger than themselves. Using an RNA ligase-mediated rapid amplification of cDNA ends method, we were able to detect large processed precursors that are generated by endonucleolytic cleavage of an even longer molecule. This cleavage produces the mature mRNA 3'-end and at the same time the 5'-end of the precursor. The tRNAs encoded by the petD and psbE minicircles appear to be processed in the same way. We propose a "rolling circle" model for chloroplast transcription in which transcription would proceed continuously around the minicircular DNA to produce transcripts larger than the minicircle itself. These transcripts would be further processed into discrete mature mRNAs and tRNAs.

Substitutional editing of Heterocapsa triquetra chloroplast transcripts and a folding model for its divergent chloroplast 16S rRNA

Substitutional editing increases genomic plasticity by changing or modifying bases at the RNA level. In this study we sequenced 10 mature chloroplast mRNAs, the chloroplast 16S rRNA and a partial chloroplast 23S rRNA from the dinoflagellate Heterocapsa triquetra, and found multiple types of substitutional editing, with A-to-G editing predominating. A-to-G editing of mRNAs converts two unusual AUA start codons into conventional AUG start codons, but three AUA start codons are not edited, showing that this dinoflagellate chloroplast has three possible start codons: AUG, AUA and UUG. To analyze the editing effects on rRNAs, we computationally predicted the secondary structure of the 16S rRNA based on the E. coli model. There are twenty editing sites in well-conserved regions of the secondary structure and eleven out of them restore conservation with other models. Moreover, A-to-G editing sites are frequently found in loop regions rather than double-stranded regions, suggesting that the A-to-G editing mechanism in dinoflagellate chloroplasts is different from that responsible for animal nuclear A-to-I(G) editing. The model of the edited 16S rRNA derived by the comparative method shares conserved secondary structural elements with other 16S rRNAs in spite of its very divergent primary sequence, supporting its role as a functional component of the chloroplast ribosome.

The replication of plastid minicircles involves rolling circle intermediates

Plastid genomes of peridinin-containing dinoflagellates are unique in that its genes are found on multiple circular DNA molecules known as ‘minicircles’ of ∼2–3 kb in size, carrying from one to three genes. The non-coding regions (NCRs) of these minicircles share a conserved core region (250–500 bp) that are AT-rich and have several inverted or direct repeats. Southern blot analysis using an NCR probe, after resolving a dinoflagellate whole DNA extract in pulsed-field gel electrophoresis (PFGE), revealed additional positive bands (APBs) of 6–8 kb in size. APBs preferentially diminished from cells treated with the DNA-replication inhibitor aphidicolin, when compared with 2–3 kb minicircles, implicating they are not large minicircles. The APBs are also exonuclease III-sensitive, implicating the presence of linear DNA. These properties and the migration pattern of the APBs in a 2D-gel electrophoresis were in agreement with a rolling circle type of replication, rather than the bubble-forming type. Atomic force microscopy of 6–8 kb DNA separated by PFGE revealed DNA intermediates with rolling circle shapes. Accumulating data thus supports the involvement of rolling circle intermediates in the replication of the minicircles.

The origin of plastids

It is generally accepted that plastids first arose by acquisition of photosynthetic prokaryotic endosymbionts by non-photosynthetic eukaryotic hosts. It is also accepted that photosynthetic eukaryotes were acquired on several occasions as endosymbionts by non-photosynthetic eukaryote hosts to form secondary plastids. In some lineages, secondary plastids were lost and new symbionts were acquired, to form tertiary plastids. Most recent work has been interpreted to indicate that primary plastids arose only once, referred to as a 'monophyletic' origin. We critically assess the evidence for this. We argue that the combination of Ockham's razor and poor taxon sampling will bias studies in favour of monophyly. We discuss possible concerns in phylogenetic reconstruction from sequence data. We argue that improved understanding of lineage-specific substitution processes is needed to assess the reliability of sequence-based trees. Improved understanding of the timing of the radiation of present-day cyanobacteria is also needed. We suggest that acquisition of plastids is better described as the result of a process rather than something occurring at a discrete time, and describe the 'shopping bag' model of plastid origin. We argue that dinoflagellates and other lineages provide evidence in support of this.

Dinoflagellates: a mitochondrial genome all at sea

Dinoflagellate algae are notorious for their highly unusual organization of nuclear and chloroplast genomes. Early studies on the dinoflagellate mitochondrial genome indicated that it encodes the same three protein-coding genes found in Plasmodium spp., but with a complex organization and transcript editing. Recent work has extended this view, showing that the dinoflagellate mitochondrial genome contains a wide array of gene fragments and genes interspersed with noncoding inverted repeats. The genome seems to require noncanonical start and stop codons, as well as high levels of editing, trans-splicing and the addition of oligonucleotide caps at the 5' and 3' ends of transcripts. Despite its small coding content, the dinoflagellate mitochondrial genome is one of the most complex known.