Chloroplast genomes as a tool to resolve red algal phylogenies: a case study in the Nemaliales

Description of Naviculavanseea sp. nov. (Naviculales, Naviculaceae), a new species of diatom from the highly alkaline Lake Van (Republic of Türkiye) with complete characterisation of its organellar genomes and multigene phylogeny

The current article describes Naviculavanseeasp. nov., a new species of diatom from Lake Van, a highly alkaline lake in Eastern Anatolia (Türkiye). The description is based on light and scanning electron microscopy performed on two monoclonal cultures. The complete nuclear rRNA clusters and plastid genomes have been sequenced for these two strains and the complete mitogenome for one of them. The plastome of both strains shows the probable loss of a functional ycf35 gene. They also exhibit two IB4 group I introns in their rrl, each encoding for a putative LAGLIDADG homing endonuclease, with the first L1917 IB4 intron reported amongst diatoms. The Maximum Likelihood phylogeny inferred from a concatenated alignment of 18S, rbcL and psbC distinguishes N.vanseea sp. nov. from the morphologically similar species Naviculacincta and Naviculamicrodigitoradiata.

Complete chloroplast genome of the marine red alga Rhodochorton tenue (Rhodochortonaceae, Rhodophyta) from San Juan Island, Washington

We present the complete chloroplast genome sequence of Rhodochorton tenue from San Juan Island, Washington. The chloroplast genome of R. tenue is 192,037 bp in length, contains 244 genes, and is similar in content to Acrochaetium secundatum. Rhodochorton tenue is genetically distinct from Rhodochorton purpureum from the North Atlantic Ocean.

A novel class of xylanases specifically degrade marine red algal β1,3/1,4-mixed-linkage xylan

Xylans are polysaccharides composed of xylose and include β1,4-xylan, β1,3-xylan, and β1,3/1,4-mixed-linkage xylan (MLX). MLX is widely present in marine red algae and constitutes a significant organic carbon in the ocean. Xylanases are hydrolase enzymes that play an important role in xylan degradation. While a variety of β1,4-xylanases and β1,3-xylanases involved in the degradation of β1,4-xylan and β1,3-xylan have been reported, no specific enzyme has yet been identified that degrades MLX. Herein, we report the characterization of a new MLX-specific xylanase from the marine bacterium Polaribacter sp. Q13 which utilizes MLX for growth. The bacterium secretes xylanases to degrade MLX, among which is Xyn26A, an MLX-specific xylanase that shows low sequence similarities (<27%) to β1,3-xylanases in the glycoside hydrolase family 26 (GH26). We show that Xyn26A attacks MLX precisely at β1,4-linkages, following a β1,3-linkage toward the reducing end. We confirm that Xyn26A and its homologs have the same specificity and mode of action on MLX, and thus represent a new xylanase group which we term as MLXases. We further solved the structure of a representative MLXase, AlXyn26A. Structural and biochemical analyses revealed that the specificity of MLXases depends critically on a precisely positioned β1,3-linkage at the -2/-1 subsite. Compared to the GH26 β1,3-xylanases, we found MLXases have evolved a tunnel-shaped cavity that is fine-tuned to specifically recognize and hydrolyze MLX. Overall, this study offers a foremost insight into MLXases, shedding light on the biochemical mechanism of bacterial degradation of MLX.

Phylogenomic analysis of pseudocryptic diversity reveals the new genus Deltalsia (Rhodomelaceae, Rhodophyta)

Molecular analyses, in combination with morphological studies, provide invaluable tools for delineating red algal taxa. However, molecular datasets are incomplete and taxonomic revisions are often required once additional species or populations are sequenced. The small red alga Conferva parasitica was described from the British Isles in 1762 and then reported from other parts of Europe. Conferva parasitica was traditionally included in the genus Pterosiphonia (type species P. cloiophylla in Schmitz and Falkenberg 1897), based on its morphological characters, and later transferred to Symphyocladia and finally to Symphyocladiella using molecular data from an Iberian specimen. However, although morphological differences have been observed between specimens of Symphyocladiella parasitica from northern and southern Europe they have yet to be investigated in a phylogenetic context. In this study, we collected specimens from both regions, studied their morphology and analyzed rbcL and cox1 DNA sequences. We determined the phylogenetic position of a British specimen using a phylogenomic approach based on mitochondrial and plastid genomes. Northern and southern European populations attributed to S. parasitica represent different species. Symphyocladiella arecina sp. nov. is proposed for specimens from southern Europe, but British specimens were resolved as a distant sister lineage to the morphologically distinctive Amplisiphonia, so we propose the new genus Deltalsia for this species. Our study highlights the relevance of using materials collected close to the type localities for taxonomic reassessments, and showcases the utility of genome-based phylogenies for resolving classification issues in the red algae.

Divergence times and plastid phylogenomics within the intron-rich order Erythropeltales (Compsopogonophyceae, Rhodophyta)

The advent of high-throughput sequencing (HTS) has allowed for the use of large numbers of coding regions to produce robust phylogenies. These phylogenies have been used to highlight relationships at ancient diversifications (subphyla, class) and highlight the evolution of plastid genome structure. The Erythropeltales are an order in the Compsopogonophyceae, a group with unusual plastid genomes but with low taxon sampling. We use HTS to produce near complete plastid genomes of all genera, and multiple species within some genera, to produce robust phylogenies to investigate character evolution, dating of divergence in the group, and plastid organization, including intron patterns. Our results produce a fully supported phylogeny of the genera in the Erythropeltales and suggest that morphologies (upright versus crustose) have evolved multiple times. Our dated phylogeny also indicates that the order is very old (~800 Ma), with diversification occurring after the ice ages of the Cryogenian period (750-635 Ma). Plastid gene order is congruent with phylogenetic relationships and suggests that genome architecture does not change often. Our data also highlight the abundance of introns in the plastid genomes of this order. We also produce a nearly complete plastid genome of Tsunamia transpacifica (Stylonematophyceae) to add to the taxon sampling of genomes of this class. The use of plastid genomes clearly produces robust phylogenetic relationships that can be used to infer evolutionary events, and increased taxon sampling, especially in less well-known red algal groups, will provide additional insights into their evolution.

Six Newly Sequenced Chloroplast Genomes From Trentepohliales: The Inflated Genomes, Alternative Genetic Code and Dynamic Evolution

Cephaleuros is often known as an algal pathogen with 19 taxonomically valid species, some of which are responsible for red rust and algal spot diseases in vascular plants. No chloroplast genomes have yet been reported in this genus, and the limited genetic information is an obstacle to understanding the evolution of this genus. In this study, we sequenced six new Trentepohliales chloroplast genomes, including four Cephaleuros and two Trentepohlia. The chloroplast genomes of Trentepohliales are large compared to most green algae, ranging from 216 to 408 kbp. They encode between 93 and 98 genes and have a GC content of 26-36%. All new chloroplast genomes were circular-mapping and lacked a quadripartite structure, in contrast to the previously sequenced Trentepohlia odorata, which does have an inverted repeat. The duplicated trnD -GTC, petD, and atpA genes in C. karstenii may be remnants of the IR region and shed light on its reduction. Chloroplast genes of Trentepohliales show elevated rates of evolution, strong rearrangement dynamics and several genes display an alternative genetic code with reassignment of the UGA/UAG codon presumably coding for arginine. Our results present the first whole chloroplast genome of the genus Cephaleuros and enrich the chloroplast genome resources of Trentepohliales.

The complete chloroplast genome and characteristics analysis of Callistemon rigidus R.Br

Callistemon rigidus R.Br. one of the traditional Chinese medicinal plants, is acrid-flavored and mild-natured, with the prominent effects reducing swelling, resolving phlegm, and dispelling rheumatism. Clinically, it has been commonly used to treat cold, cough and asthma, pain and swelling from impact injuries, eczema, rheumatic arthralgia. The chloroplast genome study on Callistemon rigidus R.Br. is a few seen. This study demonstrates the data collected from the assembly and annotation of the chloroplast (cp) genome of Callistemon rigidus R.Br., followed by furthers comparative analysis with the cp genomes of closely related species. C. rigidus R.Br. showed a cp genome in the size of 158, 961 bp long with 36.78% GC content, among which a pair of inverted repeats (IRs) of 26, 671 bp separated a large single-copy (LSC) region of 87, 162 bp and a small single-copy (SSC) region of 18, 457 bp. Altogether 131 genes were hosted, including 37 transfer RNAs, 8 ribosomal RNAs, and 86 protein-coding genes. 284 simple sequence repeats (SSRs) were also marked out. A comparative analysis of the genome structure and the sequence data of closely related species unveiled the conserved gene order in the IR and LSC/SSC regions, a quite constructive finding for future phylogenetic research. Overall, this study providing C. rigidus R.Br. genomic resources could positively contribute to the evolutionary study and the phylogenetic reconstruction of Myrtaceae.

In Silico Analysis of ACE Inhibitory Peptides from Chloroplast Proteins of Red Alga Grateloupia asiatica

Inhibition of angiotensin I-converting enzyme (ACE) is one of the key factors to repress high blood pressure. Although many studies have been reported that seaweed protein hydrolysates showed the ACE inhibitory activity, the comprehensive understanding of the relationship was still unclear. In this study, we employed chloroplast genome for in silico analysis and compared it with in vitro experiments. We first extracted water-soluble proteins (WSP) from red alga Grateloupia asiatica, which contained mainly PE, PC, APC, and Rbc, and prepared WSP hydrolysate by thermolysin, resulting that the hydrolysate showed ACE inhibitory activity. Then, we determined the complete chloroplast genome of G. asiatica (187,518 bp: 206 protein-coding genes, 29 tRNA, and 3 rRNA) and clarified the amino acid sequences of main WSP, i.e., phycobiliproteins and Rubisco, to perform in silico analysis. Consequently, 190 potential ACE inhibitory peptides existed in the main WSP sequences, and 21 peptides were obtained by in silico thermolysin digestion. By comparing in vitro and in silico analyses, in vitro ACE inhibitory activity was correlated to the IC₅₀ value from in silico digestion. Therefore, in silico approach provides insight into the comprehensive understanding of the potential bioactive peptides from seaweed proteins.

Relative expression analysis of light-harvesting genes in the freshwater alga Lympha mucosa (Batrachospermales, Rhodophyta)

Members of the freshwater red algal order Batrachospermales are often described as shade-adapted. Nevertheless, recent ecophysiological studies have demonstrated species-level differences in acclimation to a range of irradiances. Lympha mucosa occurs in open and shaded portions of temperate streams and is abundant during summer months, suggesting it tolerates high and low irradiances. Specimens of L. mucosa were collected from open (sun-acclimated) or shaded (shade-acclimated) sites and exposed to low (<20 μmol photons · m^-2 · s^-1 ) or high (220 μmol photon · m^-2 · s^-1 ) light for 72 h to examine mechanisms of photoacclimation at the transcriptional level. High-throughput sequence data were used to design specific primers for genes involved with light harvesting and these were quantified with qPCR. The greatest significant difference in transcript abundances was observed in the psaA gene (Photosystem I P700 apoprotein), and site-type had an effect on these responses. Shade-acclimated thalli were 22-fold down-regulated at high light, whereas sun-acclimated thalli were only 5-fold down-regulated. Another gene involved with Photosystem I (petF ferredoxin) was down-regulated at high light, but only individuals from the shaded site were significantly different (4-fold). In thalli from both sites, cpeA (Phycoerythrin alpha chain) was down-regulated at high light. Although not statistically significant, patterns consistent with previous physiological and transcriptomic studies were uncovered, namely the inverse response of transcriptional activity in genes that encode phycobiliproteins. In support of previous ecophysiological studies of freshwater red algae, these data indicate significant transcriptional changes involving Photosystem I and phycobiliprotein synthesis are required to tolerate and grow at various irradiances.

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.

Morphological evolution and classification of the red algal order Ceramiales inferred using plastid phylogenomics

The order Ceramiales contains about one third of red algal diversity and it was classically classified into four families according to morphology. The first phylogenies based on one or two molecular markers were poorly supported and failed to resolve these families as monophyletic. Nine families are currently recognized, but relationships within and among them are poorly understood. We produced a well-resolved phylogeny for the Ceramiales using plastid genomes for 80 (28 newly sequenced) representative species of the major lineages. Three of the previously recognized families were resolved as independent monophyletic lineages: Ceramiaceae, Wrangeliaceae and Rhodomelaceae. By contrast, our results indicated that the other six families require reclassification. We propose the new order Inkyuleeales, a new circumscription of the Callithamniaceae to include the Spyridiaceae, and a new concept of the Delesseriaceae that includes the Sarcomeniaceae and the Dasyaceae. We also investigated the evolution of the thallus structure, which has been important in the classical delineation of families. The ancestor of the Ceramiales was a monosiphonous filament that evolved into more complex morphologies several times independently during the evolutionary history of this hyperdiverse lineage.

The Complete Plastid Genome of Magnolia zenii and Genetic Comparison to Magnoliaceae species

Magnolia zenii is a critically endangered species known from only 18 trees that survive on Baohua Mountain in Jiangsu province, China. Little information is available regarding its molecular biology, with no genomic study performed on M. zenii until now. We determined the complete plastid genome of M. zenii and identified microsatellites. Whole sequence alignment and phylogenetic analysis using BI and ML methods were also conducted. The plastome of M. zenii was 160,048 bp long with 39.2% GC content and included a pair of inverted repeats (IRs) of 26,596 bp that separated a large single-copy (LSC) region of 88,098 bp and a small single-copy (SSC) region of 18,757 bp. One hundred thirty genes were identified, of which 79 were protein-coding genes, 37 were transfer RNAs, and eight were ribosomal RNAs. Thirty seven simple sequence repeats (SSRs) were also identified. Comparative analyses of genome structure and sequence data of closely-related species revealed five mutation hotspots, useful for future phylogenetic research. Magnolia zenii was placed as sister to M. biondii with strong support in all analyses. Overall, this study providing M. zenii genomic resources will be beneficial for the evolutionary study and phylogenetic reconstruction of Magnoliaceae.

Organellar genomics: a useful tool to study evolutionary relationships and molecular evolution in Gracilariaceae (Rhodophyta)

Gracilariaceae has a worldwide distribution including numerous economically important species. We applied high-throughput sequencing to obtain organellar genomes (mitochondria and chloroplast) from 10 species of Gracilariaceae and, combined with published genomes, to infer phylogenies and compare genome architecture among species representing main lineages. We obtained similar topologies between chloroplast and mitochondrial genomes phylogenies. However, the chloroplast phylogeny was better resolved with full support. In this phylogeny, Melanthalia intermedia is sister to a monophyletic clade including Gracilaria and Gracilariopsis, which were both resolved as monophyletic genera. Mitochondrial and chloroplast genomes were highly conserved in gene synteny, and variation mainly occurred in regions where insertions of plasmid-derived sequences (PDS) were found. In mitochondrial genomes, PDS insertions were observed in two regions where the transcription direction changes: between the genes cob and trnL, and trnA and trnN. In chloroplast genomes, PDS insertions were in different positions, but generally found between psdD and rrs genes. Gracilariaceae is a good model system to study the impact of PDS in genome evolution due to the frequent presence of these insertions in organellar genomes. Furthermore, the bacterial leuC/leuD operon was found in chloroplast genomes of Gracilaria tenuistipitata, G. chilensis, and M. intermedia, and in extrachromosomal plasmid of G. vermiculophylla. Phylogenetic trees show two different origins of leuC/leuD: genes found in chloroplast and plasmid were placed with proteobacteria, and genes encoded in the nucleus were close to Viridiplantae and cyanobacteria.

Reassessment of the classification of Bryopsidales (Chlorophyta) based on chloroplast phylogenomic analyses

The Bryopsidales is a morphologically diverse group of mainly marine green macroalgae characterized by a siphonous structure. The order is composed of three suborders - Ostreobineae, Bryopsidineae, and Halimedineae. While previous studies improved the higher-level classification of the order, the taxonomic placement of some genera in Bryopsidineae (Pseudobryopsis and Lambia) as well as the relationships between the families of Halimedineae remains uncertain. In this study, we re-assess the phylogeny of the order with datasets derived from chloroplast genomes, drastically increasing the taxon sampling by sequencing 32 new chloroplast genomes. The phylogenies presented here provided good support for the major lineages (suborders and most families) in Bryopsidales. In Bryopsidineae, Pseudobryopsis hainanensis was inferred as a distinct lineage from the three established families allowing us to establish the family Pseudobryopsidaceae. The Antarctic species Lambia antarctica was shown to be an early-branching lineage in the family Bryopsidaceae. In Halimedineae, we revealed several inconsistent phylogenetic positions of macroscopic taxa, and several entirely new lineages of microscopic species. A new classification scheme is proposed, which includes the merger of the families Pseudocodiaceae, Rhipiliaceae and Udoteaceae into a more broadly circumscribed Halimedaceae, and the establishment of tribes for the different lineages found therein. In addition, the deep-water genus Johnson-sea-linkia, currently placed in Rhipiliopsis, was reinstated based on our phylogeny.

Plastid genome analysis of three Nemaliophycidae red algal species suggests environmental adaptation for iron limited habitats

The red algal subclass Nemaliophycidae includes both marine and freshwater taxa that contribute to more than half of the freshwater species in Rhodophyta. Given that these taxa inhabit diverse habitats, the Nemaliophycidae is a suitable model for studying environmental adaptation. For this purpose, we characterized plastid genomes of two freshwater species, Kumanoa americana (Batrachospermales) and Thorea hispida (Thoreales), and one marine species Palmaria palmata (Palmariales). Comparative genome analysis identified seven genes (ycf34, ycf35, ycf37, ycf46, ycf91, grx, and pbsA) that were different among marine and freshwater species. Among currently available red algal plastid genomes (127), four genes (pbsA, ycf34, ycf35, ycf37) were retained in most of the marine species. Among these, the pbsA gene, known for encoding heme oxygenase, had two additional copies (HMOX1 and HMOX2) that were newly discovered and were reported from previously red algal nuclear genomes. Each type of heme oxygenase had a different evolutionary history and special modifications (e.g., plastid targeting signal peptide). Based on this observation, we suggest that the plastid-encoded pbsA contributes to the iron controlling system in iron-deprived conditions. Thus, we highlight that this functional requirement may have prevented gene loss during the long evolutionary history of red algal plastid genomes.

Organization of plastid genomes in the freshwater red algal order Batrachospermales (Rhodophyta)

Little is known about genome organization in members of the order Batrachospermales, and the infra-ordinal relationship remains unresolved. Plastid (cp) genomes of seven members of the freshwater red algal order Batrachospermales were sequenced, with the following aims: (i) to describe the characteristics of cp genomes and compare these with other red algal groups; (ii) to infer the phylogenetic relationships among these members to better understand the infra-ordinal classification. Cp genomes of Batrachospermales are large, with several cases of gene loss, they are gene-dense (high gene content for the genome size and short intergenic regions) and have highly conserved gene order. Phylogenetic analyses based on concatenated nucleotide genome data roughly supports the current taxonomic system for the order. Comparative analyses confirm data for members of the class Florideophyceae that cp genomes in Batrachospermales is highly conserved, with little variation in gene composition. However, relevant new features were revealed in our study: genome sizes in members of Batrachospermales are close to the lowest values reported for Florideophyceae; differences in cp genome size within the order are large in comparison with other orders (Ceramiales, Gelidiales, Gracilariales, Hildenbrandiales, and Nemaliales); and members of Batrachospermales have the lowest number of protein-coding genes among the Florideophyceae. In terms of gene loss, apcF, which encodes the allophycocyanin beta subunit, is absent in all sequenced taxa of Batrachospermales. We reinforce that the interordinal relationships between the freshwater orders Batrachospermales and Thoreales within the Nemaliophycidae is not well resolved due to limited taxon sampling.

Complete sequence and comparative analysis of the chloroplast genome of Plinia trunciflora

Plinia trunciflora is a Brazilian native fruit tree from the Myrtaceae family, also known as jaboticaba. This species has great potential by its fruit production. Due to the high content of essential oils in their leaves and of anthocyanins in the fruits, there is also an increasing interest by the pharmaceutical industry. Nevertheless, there are few studies focusing on its molecular biology and genetic characterization. We herein report the complete chloroplast (cp) genome of P. trunciflora using high-throughput sequencing and compare it to other previously sequenced Myrtaceae genomes. The cp genome of P. trunciflora is 159,512 bp in size, comprising inverted repeats of 26,414 bp and single-copy regions of 88,097 bp (LSC) and 18,587 bp (SSC). The genome contains 111 single-copy genes (77 protein-coding, 30 tRNA and four rRNA genes). Phylogenetic analysis using 57 cp protein-coding genes demonstrated that P. trunciflora, Eugenia uniflora and Acca sellowiana form a cluster with closer relationship to Syzygium cumini than with Eucalyptus. The complete cp sequence reported here can be used in evolutionary and population genetics studies, contributing to resolve the complex taxonomy of this species and fill the gap in genetic characterization.

Analysis of chloroplast genomes and a supermatrix inform reclassification of the Rhodomelaceae (Rhodophyta)

With over a thousand species, the Rhodomelaceae is the most species-rich family of red algae. While its genera have been assigned to 14 tribes, the high-level classification of the family has never been evaluated with a molecular phylogeny. Here, we reassess its classification by integrating genome-scale phylogenetic analysis with observations of the morphological characters of clades. In order to resolve relationships among the main lineages of the family we constructed a phylogeny with 55 chloroplast genomes (52 newly determined). The majority of branches were resolved with full bootstrap support. We then added 266 rbcL, 125 18S rRNA gene and 143 cox1 sequences to construct a comprehensive phylogeny containing nearly half of all known species in the family (407 species in 89 genera). These analyses suggest the same subdivision into higher-level lineages, but included many branches with moderate or poor support. The circumscription for nine of the 13 previously described tribes was supported, but the Lophothalieae, Polysiphonieae, Pterosiphonieae and Herposiphonieae required revision, and five new tribes and one resurrected tribe were segregated from them. Rhizoid anatomy is highlighted as a key diagnostic character for the morphological delineation of several lineages. This work provides the most extensive phylogenetic analysis of the Rhodomelaceae to date and successfully resolves the relationships among major clades of the family. Our data show that organellar genomes obtained through high-throughput sequencing produce well-resolved phylogenies of difficult groups, and their more general application in algal systematics will likely permit deciphering questions about classification at many taxonomic levels.

Phylogenetic position of the coral symbiont Ostreobium (Ulvophyceae) inferred from chloroplast genome data

The green algal genus Ostreobium is an important symbiont of corals, playing roles in reef decalcification and providing photosynthates to the coral during bleaching events. A chloroplast genome of a cultured strain of Ostreobium was available, but low taxon sampling and Ostreobium's early-branching nature left doubt about its phylogenetic position. Here, we generate and describe chloroplast genomes from four Ostreobium strains as well as Avrainvillea mazei and Neomeris sp., strategically sampled early-branching lineages in the Bryopsidales and Dasycladales respectively. At 80,584 bp, the chloroplast genome of Ostreobium sp. HV05042 is the most compact yet found in the Ulvophyceae. The Avrainvillea chloroplast genome is ~94 kbp and contains introns in infA and cysT that have nearly complete sequence identity except for an open reading frame (ORF) in infA that is not present in cysT. In line with other bryopsidalean species, it also contains regions with possibly bacteria-derived ORFs. The Neomeris data did not assemble into a canonical circular chloroplast genome but a large number of contigs containing fragments of chloroplast genes and showing evidence of long introns and intergenic regions, and the Neomeris chloroplast genome size was estimated to exceed 1.87 Mb. Chloroplast phylogenomics and 18S nrDNA data showed strong support for the Ostreobium lineage being sister to the remaining Bryopsidales. There were differences in branch support when outgroups were varied, but the overall support for the placement of Ostreobium was strong. These results permitted us to validate two suborders and introduce a third, the Ostreobineae.

Origin and evolutionary history of freshwater Rhodophyta: further insights based on phylogenomic evidence

Freshwater representatives of Rhodophyta were sampled and the complete chloroplast and mitochondrial genomes were determined. Characteristics of the chloroplast and mitochondrial genomes were analyzed and phylogenetic relationship of marine and freshwater Rhodophyta were reconstructed based on the organelle genomes. The freshwater member Compsopogon caeruleus was determined for the largest chloroplast genome among multicellular Rhodophyta up to now. Expansion and subsequent reduction of both the genome size and GC content were observed in the Rhodophyta except for the freshwater Compsopogon caeruleus. It was inferred that the freshwater members of Rhodophyta occurred through diverse origins based on evidence of genome size, GC-content, phylogenomic analysis and divergence time estimation. The freshwater species Compsopogon caeruleus and Hildenbrandia rivularis originated and evolved independently at the inland water, whereas the Bangia atropurpurea, Batrachospermum arcuatum and Thorea hispida are derived from the marine relatives. The typical freshwater representatives Thoreales and Batrachospermales are probably derived from the marine relative Palmaria palmata at approximately 415–484 MYA. The origin and evolutionary history of freshwater Rhodophyta needs to be testified with more organelle genome sequences and wider global sampling.