Gene-rich plastid genomes of two parasitic red algal species,Laurencia australis and L. verruciformis (Rhodomelaceae,Ceramiales), and a taxonomic revision of Janczewskia

Parasitic red algae are an interesting system for investigating the genetic changes that occur in parasites. These parasites have evolved independently multiple times within the red algae. The functional loss of plastid genomes can be investigated in these multiple independent examples, and fine-scale patterns may be discerned. The only plastid genomes from red algal parasites known so far are highly reduced and missing almost all photosynthetic genes. Our study assembled and annotated plastid genomes from the parasites Janczewskia tasmanica and its two Laurencia host species (Laurencia elata and one unidentified Laurencia sp. A25) from Australia and Janczewskia verruciformis, its host species (Laurencia catarinensis), and the closest known free-living relative (Laurencia obtusa) from the Canary Islands (Spain). For the first time we show parasitic red algal plastid genomes that are similar in size and gene content to free-living host species without any gene loss or genome reduction. The only exception was two pseudogenes (moeB and ycf46) found in the plastid genome of both isolates of J. tasmanica, indicating potential for future loss of these genes. Further comparative analyses with the three highly reduced plastid genomes showed possible gene loss patterns, in which photosynthetic gene categories were lost followed by other gene categories. Phylogenetic analyses did not confirm monophyly of Janczewskia, and the genus was subsumed into Laurencia. Further investigations will determine if any convergent small-scale patterns of gene loss exist in parasitic red algae and how these are applicable to other parasitic systems.     

Reappraising plastid markers of the red algae for phylogenetic community ecology in the genomic era

Selection of appropriate genetic markers to quantify phylogenetic diversity is crucial for community ecology studies. Yet, systematic evaluation of marker genes for this purpose is scarcely done. Recently, the combined effort of phycologists has produced a rich plastid genome resource with taxonomic representation spanning all of the major lineages of the red algae (Rhodophyta). In this proof‐of‐concept study, we leveraged this resource by developing and applying a phylogenomic strategy to seek candidate plastid markers suitable for phylogenetic community analysis. We ranked the core genes of 107 published plastid genomes based on various sequence‐derived properties and their tree distance to plastid genome phylogenies. The resulting ranking revealed that the most widely used marker, rbcL, is not necessarily the optimal marker, while other promising markers might have been overlooked. We designed and tested PCR primers for several candidate marker genes, and successfully amplified one of them, rpoC1, in a taxonomically broad set of red algal specimens. We suggest that our general marker identification methodology and the rpoC1 primers will be useful to the phycological community for investigating the biodiversity and community ecology of the red algae.     

Comparative Analysis of the Complete Plastid Genome Sequence of the Red Alga Gracilaria tenuistipitata var. liui Provides Insights into the Evolution of Rhodoplasts and Their Relationship to Other Plastids

We sequenced to completion the circular plastid genome of the red alga Gracilaria tenuistipitata var. liui. This is the first plastid genome sequence from the subclass Florideophycidae (Rhodophyta). The genome is composed of 183,883 bp and contains 238 predicted genes, including a single copy of the ribosomal RNA operon. Comparisons with the plastid genome of Porphyra pupurea reveal strong conservation of gene content and order, but we found major genomic rearrangements and the presence of coding regions that are specific to Gracilaria. Phylogenetic analysis of a data set of 41 concatenated proteins from 23 plastid and two cyanobacterial genomes support red algal plastid monophyly and a specific evolutionary relationship between the Florideophycidae and the Bangiales. Gracilaria maintains a surprisingly ancient gene content in its plastid genome and, together with other Rhodophyta, contains the most complete repertoire of plastid genes known in photosynthetic eukaryotes.Supplementary material () is available for this article.[Reviewing Editor: Dr. W. Ford Doolittle]     

Plastid proteome prediction for diatoms and other algae with secondary plastids of the red lineage

The plastids of ecologically and economically important algae from phyla such as stramenopiles, dinoflagellates and cryptophytes were acquired via a secondary endosymbiosis and are surrounded by three or four membranes. Nuclear‐encoded plastid‐localized proteins contain N‐terminal bipartite targeting peptides with the conserved amino acid sequence motif ‘ASAFAP’. Here we identify the plastid proteomes of two diatoms, Thalassiosira pseudonana and Phaeodactylum tricornutum, using a customized prediction tool (ASAFind) that identifies nuclear‐encoded plastid proteins in algae with secondary plastids of the red lineage based on the output of SignalP and the identification of conserved ‘ASAFAP’ motifs and transit peptides. We tested ASAFind against a large reference dataset of diatom proteins with experimentally confirmed subcellular localization and found that the tool accurately identified plastid‐localized proteins with both high sensitivity and high specificity. To identify nucleus‐encoded plastid proteins of T. pseudonana and P. tricornutum we generated optimized sets of gene models for both whole genomes, to increase the percentage of full‐length proteins compared with previous assembly model sets. ASAFind applied to these optimized sets revealed that about 8% of the proteins encoded in their nuclear genomes were predicted to be plastid localized and therefore represent the putative plastid proteomes of these algae.     

Mutant swarms of a totivirus‐like entities are present in the red macroalga Chondrus crispus and have been partially transferred to the nuclear genome

Chondrus crispus Stackhouse (Gigartinales) is a red seaweed found on North Atlantic rocky shores. Electrophoresis of RNA extracts showed a prominent band with a size of around 6,000 bp. Sequencing of the band revealed several sequences with similarity to totiviruses, double‐stranded RNA viruses that normally infect fungi. This virus‐like entity was named C. crispus virus (CcV). It should probably be regarded as an extreme viral quasispecies or a mutant swarm since low identity (<65%) was found between sequences. Totiviruses typically code for two genes: one capsid gene (gag) and one RNA‐dependent RNA polymerase gene (pol) with a pseudoknot structure between the genes. Both the genes and the intergenic structures were found in the CcV sequences. A nonidentical gag gene was also found in the nuclear genome of C. crispus, with associated expressed sequence tags (EST) and upstream regulatory features. The gene was presumably horizontally transferred from the virus to the alga. Similar dsRNA bands were seen in extracts from different life cycle stages of C. crispus and from all geographic locations tested. In addition, similar bands were also observed in RNA extractions from other red algae; however, the significance of this apparently widespread phenomenon is unknown. Neither phenotype caused by the infection nor any virus particles or capsid proteins were identified; thus, the presence of viral particles has not been validated. These findings increase the known host range of totiviruses to include marine red algae.     

Molecular phylogenetics supports a clade of red algal parasites retaining native plastids: taxonomy and terminology revised

Parasitism is a life strategy that has repeatedly evolved within the Florideophyceae. Historically, the terms adelphoparasite and alloparasite have been used to distinguish parasites based on the relative phylogenetic relationship of host and parasite. However, analyses using molecular phylogenetics indicate that nearly all red algal parasites infect within their taxonomic family, and a range of relationships exist between host and parasite. To date, all investigated adelphoparasites have lost their plastid, and instead, incorporate a host‐derived plastid when packaging spores. In contrast, a highly reduced plastid lacking photosynthesis genes was sequenced from the alloparasite Choreocolax polysiphoniae. Here we present the complete Harveyella mirabilis plastid genome, which has also lost genes involved in photosynthesis, and a partial plastid genome from Leachiella pacifica. The H. mirabilis plastid shares more synteny with free‐living red algal plastids than that of C. polysiphoniae. Phylogenetic analysis demonstrates that C. polysiphoniae, H. mirabilis, and L. pacifica form a robustly supported clade of parasites, which retain their own plastid genomes, within the Rhodomelaceae. We therefore transfer all three genera from the exclusively parasitic family, Choreocolacaceae, to the Rhodomelaceae. Additionally, we recommend applying the terms archaeplastic parasites (formerly alloparasites), and neoplastic parasites (formerly adelphoparasites) to distinguish red algal parasites using a biological framework rather than taxonomic affiliation with their hosts.     

PHYLOGENY OF GRACILARIA LEMANEIFORMIS (RHODOPHYTA) BASED ON SEQUENCE ANALYSIS OF ITS SMALL SUBUNIT RIBOSOMAL RNA CODING REGION1

The small subunit ribosomal RNA (rRNA) sequence of Gracilaria lemaneiformis Bory Weber-van Bosse was inferred from analysis of rRNA coding regions that were amplified by the polymerase chain reaction method. Comparison of the G. lemaneiformis small subunit rRNA to homologous genes of diverse eukaryotes demonstrated that the red algal divergence was nearly simultaneous with the separation of plants, fungi, animals and many other protist lineages. This result conflicts with those of 5S rRNA sequence and plastid based phytogenies which suggest that red algae represent an early divergence in the eukaryotic line of descent. Further, algae appear to be of polyphyletic origin and red algae are unrelated to higher fungi.     

Phylogeny of Gracilariaceae (Rhodophyta): evidence from plastid and mitochondrial nucleotide sequences

Gracilariaceae are mostly pantropical red algae and include ~230 species in seven genera. Infrafamilial classification of the group has long been based on reproductive characters, but previous phylogenies have shown that traditionally circumscribed groups are not monophyletic. We performed phylogenetic analyses using two plastid (universal plastid amplicon and rbcL) and one mitochondrial (cox1) loci from a greatly expanded number of taxa to better assess generic relationships and understand patterns of character distributions. Our analyses produce the most well‐supported phylogeny of the family to date, and indicate that key characteristics of spermatangia and cystocarp type do not delineate genera as commonly suggested. Our results further indicate that Hydropuntia is not monophyletic. Given their morphological overlap with closely related members of Gracilaria, we propose that Hydropuntia be synonymized with the former. Our results additionally expand the known ranges of several Gracilariaceae species to include Brazil. Lastly, we demonstrate that the recently described Gracilaria yoneshigueana should be synonymized as G. domingensis based on morphological and molecular characters. These results demonstrate the utility of DNA barcoding for understanding poorly known and fragmentary materials of cryptic red algae.     

Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond

Advanced resources for genome‐assisted research in barley (Hordeum vulgare) including a whole‐genome shotgun assembly and an integrated physical map have recently become available. These have made possible studies that aim to assess genetic diversity or to isolate single genes by whole‐genome resequencing and in silico variant detection. However such an approach remains expensive given the 5 Gb size of the barley genome. Targeted sequencing of the mRNA‐coding exome reduces barley genomic complexity more than 50‐fold, thus dramatically reducing this heavy sequencing and analysis load. We have developed and employed an in‐solution hybridization‐based sequence capture platform to selectively enrich for a 61.6 megabase coding sequence target that includes predicted genes from the genome assembly of the cultivar Morex as well as publicly available full‐length cDNAs and de novo assembled RNA‐Seq consensus sequence contigs. The platform provides a highly specific capture with substantial and reproducible enrichment of targeted exons, both for cultivated barley and related species. We show that this exome capture platform provides a clear path towards a broader and deeper understanding of the natural variation residing in the mRNA‐coding part of the barley genome and will thus constitute a valuable resource for applications such as mapping‐by‐sequencing and genetic diversity analyzes.     

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C‐169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.     

Identification of a Marine Cyanophage in a Protist Single‐cell Metagenome Assembly

Analysis of microbial biodiversity is hampered by a lack of reference genomes from most bacteria, viruses, and algae. This necessitates either the cultivation of a restricted number of species for standard sequencing projects or the analysis of highly complex environmental DNA metagenome data. Single‐cell genomics (SCG) offers a solution to this problem by constraining the studied DNA sample to an individual cell and its associated symbionts, prey, and pathogens. We used SCG to study marine heterotrophic amoebae related to Paulinella ovalis (A. Wulff) P.W. Johnson, P.E. Hargraves & J.M. Sieburth (Rhizaria). The genus Paulinella is best known for its photosynthetic members such as P. chromatophora Lauterborn that is the only case of plastid primary endosymbiosis known outside of algae and plants. Here, we studied the phagotrophic sister taxa of P. chromatophora that are related to P. ovalis and found one SCG assembly to contain α‐cyanobacterial DNA. These cyanobacterial contigs are presumably derived from prey. We also uncovered an associated cyanophage lineage (provisionally named phage PoL_MC2). Phylogenomic analysis of the fragmented genome assembly suggested a minimum genome size of 200 Kbp for phage PoL_MC2 that encodes 179 proteins and is most closely related to Synechococcus phage S‐SM2. For this phage, gene network analysis demonstrates a highly modular genome structure typical of other cyanophages. Our work demonstrates that SCG is a powerful approach for discovering algal and protist biodiversity and for elucidating biotic interactions in natural samples.     

Comprehensive definition of genome features in Spirodela polyrhiza by high‐depth physical mapping and short‐read DNA sequencing strategies

Spirodela polyrhiza is a fast‐growing aquatic monocot with highly reduced morphology, genome size and number of protein‐coding genes. Considering these biological features of Spirodela and its basal position in the monocot lineage, understanding its genome architecture could shed light on plant adaptation and genome evolution. Like many draft genomes, however, the 158‐Mb Spirodela genome sequence has not been resolved to chromosomes, and important genome characteristics have not been defined. Here we deployed rapid genome‐wide physical maps combined with high‐coverage short‐read sequencing to resolve the 20 chromosomes of Spirodela and to empirically delineate its genome features. Our data revealed a dramatic reduction in the number of the rDNA repeat units in Spirodela to fewer than 100, which is even fewer than that reported for yeast. Consistent with its unique phylogenetic position, small RNA sequencing revealed 29 Spirodela‐specific microRNA, with only two being shared with Elaeis guineensis (oil palm) and Musa balbisiana (banana). Combining DNA methylation data and small RNA sequencing enabled the accurate prediction of 20.5% long terminal repeats (LTRs) that doubled the previous estimate, and revealed a high Solo:Intact LTR ratio of 8.2. Interestingly, we found that Spirodela has the lowest global DNA methylation levels (9%) of any plant species tested. Taken together our results reveal a genome that has undergone reduction, likely through eliminating non‐essential protein coding genes, rDNA and LTRs. In addition to delineating the genome features of this unique plant, the methodologies described and large‐scale genome resources from this work will enable future evolutionary and functional studies of this basal monocot family.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

Discrimination and Characterization of Two Mediterranean Species from the Laurencia Complex (Rhodomelacea) Using an NMR‐Based Metabolomic Approach

Generic and specific determination among the Laurencia complex is a challenging task. DNA barcoding combined with phenotypic investigations are mandatory for species differentiation. In this study, two morphologically different members of the Laurencia complex were investigated using untargeted ¹H‐NMR‐based metabolomics. Twenty‐one population samples were collected in order to evaluate both temporal and geographical homogeneity. Data obtained from ¹H‐NMR analysis followed by statistical analysis allowed a clear separation of all the samples into two groups. DNA mitochondrial tests confirmed this pattern and identified the two species as Laurenciella sp. and Laurencia obtusa. In addition, metabolites responsible of this discrimination were investigated directly in crude extracts by ¹³C‐NMR using an in‐house computer‐assisted method. The combination of both untargeted (¹H) and targeted (¹³C) NMR‐based metabolomic approaches proves to be a powerful and complementary approach to discriminate species from the Laurencia complex.