Chloroplast Genomes of Five Oedogonium Species: Genome Structure, Phylogenetic Analysis and Adaptive Evolution

doi:10.21203/rs.3.rs-449423/v1

Download PDF

Research Article

Chloroplast Genomes of Five Oedogonium Species: Genome Structure, Phylogenetic Analysis and Adaptive Evolution

https://doi.org/10.21203/rs.3.rs-449423/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background

The order Oedogoniales can be divided into three genera, Oedogonium, Oedocladium, and Bulbochaete based on traditional morphological criteria. While several molecular phylogenetic studies have suggested that both Oedogonium and Oedocladium may not be monophyletic, broader taxon sampling and large amounts of molecular data acquisition could help to resolve the phylogeny and evolutionary problems of this order. This study determined five chloroplast (cp) genomes of Oedogonium species and aimed to provide further information on cp genome for a better understanding of the phylogenetic and evolutionary relationships of the order Oedogoniales.

Results

The five Oedogonium cp genomes showed typical quadripartite and circular structures, and were relatively conserved in their structure, gene synteny, and inverted repeats boundaries in general, except for small variation in genome sizes, AT contents, introns, and repeats. Phylogenetic analyses based on 54 cp protein-coding genes examined by maximum likelihood and Bayesian analyses using amino acid and nucleotide datasets indicated that both Oedocladium and Oedogonium are polyphyletic groups. A positively selected gene (psbA) was identified in the two Oedocladium species and the terrestrial Oedogonium species, indicating that terrestrial Oedogoniales taxa may have undergone adaptive evolution to adjust to the difference in light intensity between aquatic and terrestrial habitats.

Conclusions

Our results enrich the data on cp genomes of the genus Oedogonium. The availability of these cp genomes can help in understanding the cp genome characteristics and resolve phylogenetic and evolutionary relationships of the order Oedogoniales.

Epigenetics & Genomics

Oedogoniales

chloroplast genome

synteny

polyphyletic group

adaptive evolution

The order Oedogoniales includes three genera: Oedogonium Link ex Hirn, Oedocladium Stahl, and Bulbochaete Agardh [1-4]. More than 600 species have been described in this order, most of which can be found in fresh waters throughout the world, although Oedocladium species are mainly found on soil surfaces, a few species of Oedogonium are found in moist soil surfaces. The presence of branches and hairs are the genus-level characteristics to distinguish this order; Oedogonium has simple and unbranched filaments, Bulbochaete has bulb-based hairs, and Oedocladium has branched filaments [4-13]. While some molecular phylogenetic studies on Oedogoniales have suggested that this order is monophyletic, both Oedogonium and Oedocladium do not appear to be monophyletic, and the morphological criteria of Oedogoniales do not define natural groups, making its evolutionary position unclear [14-19].

Chloroplast (cp) genomes have been found to be ideal for phylogenetic analysis and molecular evolution studies owing to advantages such as low evolution rate and maternal inheritance [20-24], and plastome has been increasingly used for phylogenetic and evolutionary studies of green algae. For example, Claude Lemieux et al. [25] conducted cp phylogenetic analysis based on the cp genes of 61 chlorophytes and revealed that Trebouxiophyceae is not monophyletic. Zhang et al. [26] demonstrated the adaptive mechanism of sea-ice environment by analyzing the molecular evolution of an Antarctic sea ice alga Chlamydomonas sp. based on cp protein-coding genes. However, only four cp genomes of Oedogoniales are currently available in public databases [27, 28], restricting the phylogenetic analysis and molecular evolution studies based on cp genomes of this group.

Nucleotide substitution rates are often used as the criterion to reflect selection pressure. While nonsynonymous substitution rates (dN) can cause amino acid change, synonymous substitution rates (dS) do not cause amino acid change. The dN/dS ratio is the measure of natural selection acting on the protein. According to Yang [29], dN/dS < 1 denotes negative purifying selection, dN/dS = 1 signifies neutral evolution, and dN/dS > 1 indicates positive selection [30]. As most of the plastid protein-coding genes undergo negative or purifying selection to maintain their function, they are conserved and have a low dN/dS ratio. However, some genes might undergo positive selection in response to environmental changes, consequently presenting relatively high dN/dS ratio [31-34].

In this study, the cp genomes of five Oedogonium species, were sequenced and an in-depth analysis of these genomes, including comparative analysis with previously reported Oedocladium and Oedogonium cp genomes, was performed. Furthermore, phylogenetic analysis and evolutionary study of the order Oedogoniales were conducted based on cp protein-coding genes and a positively selected gene was identified in Oedocladium species The results of this study could be useful to understand the phylogenetic and evolutionary relationships of Oedogoniales.

Species identification

Description

Oedogonium dentireticulatum

Nannandrous, gynandrosporous; vegetative cells cylindrical; oogonium single, subglobose, poriferous, pore median to inframedian, oospore of the same form as oogonium, nearly or completely filling the oogonium, outer layer reticulate and dentate, teeth spreading form reticulations; suffultory cells slightly inflated or inflated; androsporangia single or up to 2 seriate, scattered; dwarf males on suffultory cells, stipes slightly curved, antheridia exteriors, 1 or 2 continuous; top and base of the filaments often slender; apical and basal cells not observed; vegetative cells generally 5–9 times as long as their width (Figture 1A–D).

Vegetative cells: 14–19 × 35–91 μm

Oogonium: 33–48 × 33–51 μm

Oospore with dentate teeth: 33–49 × 33–51 μm

Oospore without dentate teeth: 27–37 × 29–36 μm

Androsporangia: 13–15 × 12–14 μm

Antheridia: 7 × 11 μm.

Oedogonium crispum

Monoecious; oogonium usually single, obovoid-globose, operculate, division superior; oospore globose or subglobose, nearly or completely filling the oogonium, spore wall smooth; antheridium single, subepigynous; sperms 2, division horizontal; terminal cells apically obtuse; basal cells elongate; vegetative cells generally 3–4.5 times as long as their width (Figture 1E–I).

Vegetative cells: 10–14 × 24–36 μm

Oogonium: 33–38 × 30–41 μm

Oospore: 28–34 × 27–32 μm

Antheridium: 5–8 × 10–12 μm.

Oedogonium capilliforme

Dioecious, macrandrous; oogonium single, slightly inflated, rarely 2 continuous, obovoid to subovoid, with superior pore; oospore variable, ovoid-globose, subglobose, or globose, nearly or completely filling the oogonium, spore wall smooth; antheridium 2–7 in a series, often alternating with the vegetative cell; sperms 2, division horizontal; terminal cells apically apiculate; basal cells elongate; vegetative cells generally 1–3 times as long as their width (Figture 2A–E).

Female vegetative cells: 20–29 × 54–95 μm; male vegetative cells: 19–24 × 42–62 μm

Oogonium: 45–53 × 48–69 μm

Oospore: 44–49 × 42–50 μm

Antheridium: 14–17 × 9–13 μm.

The characteristics of all the five Oedogonium species were examined, and strains FACHB-3309, FACHB-3310, and FACHB-3312 were identified as Oe. dentireticulatum, Oe. crispum, and Oe. Capilliforme, respectively. With regard to strains FACHB-3311 (Figture 2F) and FACHB-3313 (Figture 2G), the entire sexual features could not be observed; however, the filaments of both of these strains were unbranched, indicating that they obviously belonged to the genus Oedogonium. In particular, strain FACHB-3313 exhibited unbranched rhizoids that resembled those of Oedocladium.

General characteristics and comparison of Oedogoniales cp genomes

Table 1 summarizes the cp genomes characteristics of the five newly included Oedogonium species, three reported Oedocladium taxa and one Oedogonium species. The complete cp genomes of the nine species of Oedogoniales ranged from 146,367 bp (Oe. crispum) to 204,438 bp (O. carolinianum) in length. All of the five Oedogonium cp genomes displayed typical circular mapping with a large single copy (LSC) region (76,475–98,887 bp), a small single copy (SSC) region (43,305–58,055 bp), and two inverted repeats (IR) regions (12,808–35,492 bp) (Supplementary Figs S1–S5). The overall AT content in each cp genome was comparable and showed a little difference among the species, ranging from 69.98% (strain FACHB-3311) to 72.66% (O. prescottii); besides, difference was noted in coding proportion, which varied from 51.4% (O. carolinianum) to 69.5% (O. prescottii). The cp genomes of six Oedogonium species were moderately compact relative to those of the Oedocladium species. The number of genes located in the plus or minus strand showed some differences. All the cp genomes contained 68 protein-coding genes and three rRNA genes, except for the cp genome of Oe. cardiacum, which had two additional genes (dpoB and int) located in the IR region. With respect to tRNA, the cp genomes showed slight difference as follows: Oe. cardiacum exhibited two additional trnR(ccu) located in the IR regions and Oe. dentireticulatum (strain FACHB-3309) presented an additional trnR(ccu) in the LSC region; Oe. sp. (strain FACHB-3313) contained two additional trnR(ccg) in the IR regions and O. carolinianum had an additional trnR(ccg) in the LSC region; and O. carolinianum had an additional trnS(gga) in the LSC region. Sequence repeats of more than 30 bp were less frequent (3.9%–4.9%) in the cp genomes of the five Oedogonium species when compared with those in the two O. carolinianum cp genomes, but were more frequent, when compared with those in the Oe. cardiacum cp genome.

Introns content and insertion sites

The introns content and insertion sites of the nine Oedogoniales cp genomes are listed in Table 1 and Supplementary Tables S1 and S2. The nine cp genomes significantly differed with respect to the introns content. Oe. sp. (strain FACHB-3311) had the maximum introns content with 26 group I introns and 11 group II introns. When compared with the other six Oedogonium cp genomes, multiple intron losses were observed in the cp genome of Oe. crispum (strain FACHB-3310), with four group I introns in trnL(uaa), psbC, atpA, and psbD, respectively, and four group II introns in psbI, petD, psaC, and psaB, respectively. Besides, similar to O. prescottii, Oe. crispum also exhibited introns losses in psbA. Oe. sp. (strain FACHB-3311) presented two additional group II introns in chlB and chlL, introns were first observed in them. All the nine cp genomes included group I introns in trnL(uaa), which is common across all algal lineages and is considered to originate from the common ancestor of cp [35]. The nine cp genomes showed a certain variation in insertion sites. The common group I introns in trnL(uaa) and group II introns in petB, psaC, and psbI (only strain FACHB-3311 lost the intron in psbI) showed the same insertion sites. With regard to the other genes with introns, the insertion sites in different species showed similarities and variations. For instance, in psbA, the number of introns (introns in psbA are all group I) differed among the species, whereas the insertion sites of the first intron in Oe. dentireticulatum (strain FACHB-3309), Oe. sp. (strain FACHB-3311), and Oe. sp. (strain FACHB-3313) were identical. The two O. carolinianum were the same; however, the insertion site of the first intron in Oe. capilliforme was similar to that of the fourth intron in Oe. dentireticulatum and sp. (strain FACHB-3311).

Synteny analysis and average nucleotide identity analysis

ProgressiveMauve was used to analyze the Oedogoniales cp genomes synteny, with Oe. cardiacum used as a reference to compare gene order among the cp genomes (Fig. 3). More than 19 locally collinear blocks (LCBs) were identified in the cp genomes of the nine species of Oedogoniales, including six taxa from Oedogonium and three taxa from Oedocladium. The nine cp genomes showed high degree of syntenic conservation on the whole, with Oe. capilliforme exhibiting high similarity to Oe. cardiacum, and Oe. dentireticulatum resembling Oe. sp. (strain FACHB-3311). However, some rearrangements and inversions were still observed among certain short LCBs mainly owing to the inversion or loss of introns. The genes order and number were almost identical except for that an inversion between trnE(uuc) and petL with a length of less than 3 kb and including the genes petD and trnR(ucg) was detected in O. carolinianum (MT364369) and O. carolinianum (NC_031510).

The average nucleotide identity (ANI) of the nine species of Oedogoniales was calculated using FastANI (Supplementary Fig. S6). Oe. crispum showed high ANI with Oe. dentireticulatum and Oe. sp. (strain FACHB-3311) (90.64% and 90.56%, respectively), Oe. dentireticulatum was similar to Oe. sp. (strain FACHB-3311) with 92.57% ANI, and Oe. capilliforme was similar to Oe. cardiacum with 97.03% ANI.

IR expansion and contraction

The IR boundary regions of the nine species of Oedogoniales were compared as illustrated in Fig. 4. Oe. cardiacum and Oe. capilliforme (strain FACHB-3312) showed larger IRs reaching 35,000 bp, whereas Oe. crispum and O. prescottii exhibited smaller IRs reaching 13,284 and 12,808 bp, respectively. The IRs of all the nine cp genomes contained the same four protein-coding genes, three tRNAs, and three rRNAs. However, in Oe. capilliforme and Oe. sp. (strain FACHB-3313), an additional trnR(ccg) was observed between psbA and rbcL; Oe. cardiacum included two additional protein-coding genes (int and dpoB) and one tRNA (trnR(ccu)); and the IRa of four cp genomes included parts of the 5’-end of ccsA (390 bp in Oe. cardiacum, Oe. capilliforme, and O. prescottii and 383 bp in Oe. crispum).

The nine Oedogoniales cp genomes showed high conservation at four regional boundaries, with little variation. The LSC/IRb junctions (JLBs) in the cp genomes of Oe. cardiacum, Oe. capilliforme, O. prescottii, and Oe. sp. (strain FACHB-3311) were located in trnR(ucu); as a result, 2 bp of the 3’-end of this gene were a part of the IR region. In Oe. sp. (strain FACHB-3311), the IR region contained 6 bp of the 3’-end of trnR(ucu), and in the other five cp genomes, the LSC/IRb boundaries occurred between trnR(ccu) and psbA. The IRb/SSC boundaries in all the nine cp genomes occurred between trnL(caa) and psaA, and the SSC/IRa junctions were located in rpoA. The IRa/LSC junctions (JSAs) of the two O. carolinianum cp genomes occurred between psbA and ccsA, while those of the other seven genomes were located in ccsA, with 390 bp of the 5’-end of this gene being a part of the IR region in Oe. cardiacum, Oe. capilliforme, and O. prescottii, and 383, 388, 608, and 389 bp of the 5’-end of this gene being a part of the IR region in Oe. crispum, strain FACHB-3313, strain FACHB-3311, and Oe. dentireticulatum, respectively.

Phylogenetic analysis and adaptive evolution analysis

Phylogenetic assays based on 54 cp protein-coding genes were conducted using maximum likelihood (ML) and Bayesian analyses with amino acid and nucleotide datasets, which generated two kinds of phylogenetic trees showing the same results (Figs 5-6). Phylogenetic trees based on amino acid and nucleotide datasets both indicated that the nine species of Oedogoniales clustered into three clades Oe. sp. (MW250873) formed a separate clade with absolute high support value, the two O. carolinianum clustered together and formed another clade, and the other six Oedogoniales formed the third clade. With regard to the third clade, the two datasets showed a little difference in the location of O. prescottii. Based on nucleotide dataset, O. prescottii clustered with Oe. cardiacum and Oe. capilliforme, whereas according to the amino acid dataset, O. prescottii clustered with Oe. dentireticulatum, Oedogonium sp. (MW250875), and Oe. crispum. These results indicated that both Oedocladium and Oedogonium are polyphyletic, which is in accordance with that reported in a previous study [36].

Positive selection analysis was performed based on branch-site model, and the null and alternative models were compared. The null model considered that the foreground branch only has dN/dS = 1, and the alternative model assumed that sites on the foreground branch have dN/dS > 1 (positive selection). When the two Oedocladium species and MW250875 were labelled as the foreground branch, the FDR-adjusted P value of psbA was less than 0.05. When the foreground branch only included the three Oedocladium spp., the FDR-adjusted P value of psbA was less than 0.05. Based on Bayes empirical Bayes (BEB) assay, the two methods were noted to indicate that psbA may possibly contain sites under positive selection, with 291SER showing posterior probability higher than 99%. However, owing to the lack of related functional sites information on closely related species such as Chlamydomonas reinhardtii and Stigeoclonium helveticum in UniProt, the positively selected sites of psbA require further investigation.

In this study, we investigated five Oedogonium isolates from China, of which strains FACHB-3309, FACHB-3310, and FACHB-3312 were identified as Oe. dentireticulatum, Oe. crispum, and Oe. capilliforme, respectively. Strains FACHB-3311 and FACHB-3313 were considered to belong to the genus Oedogonium owing to their unbranched filaments; however, they could not be identified at species level owing to their lack of entire sexual characters.

Comparative analyses of the nine Oedogoniales cp genomes showed highly conserved structures and gene numbers. The cp genomes of the newly sequenced five Oedogonium species were found to share the same structure as the previously reported Oedogoniales cp genomes, and the structures of the tetrad were not altered, but were different from the other two orders in the OCC clade (the IR is obliterated in the reported cp genomes in Chaetophorales and Floydiella of Chaetopeltidales). It has been indicated IR loss may be a synapomorphy marking the common ancestry of Chaetophorales and Chaetopeltidales [37]. The total length of these cp genomes was observed to vary within a relatively large range, extending from 146,367 bp (Oe. crispum) to 204,438 bp (O. carolinianum), which may be the result of contraction and expansion of IR regions and proportion of non-coding sequences, such as the introns content. Furthermore, the nine cp genomes showed highly conserved protein-coding genes and rRNAs number; however, they presented a slight difference in the tRNAs content. With regard to the introns content, the nine cp genomes exhibited relative variation, and the number of group I introns significantly differed, mainly owing to the diversity in the introns in psbA. In particular, introns (group II) were observed for the first time in chlB and chlB in Oe. dentireticulatum. All the nine cp genomes retained the group I introns in trnL(uaa) and group II introns in petD and psaC, and shared the same insertion sites. With regard to the other genes with introns, the insertion sites of different species showed similarities and variations.

Synteny analyses revealed a relatively high degree of syntenic conservation among the nine cp genomes, and only one inversed segment was detected in O. carolinianum FACHB-2453 and O. carolinianum UTEX LB 1686. The other variations were mainly owing to the introns, and no structural variation was observed in the six Oedogonium species. The results of FastANI also supported the findings of synteny analyses, indicating that Oe. capilliforme had high similarity with Oe. cardiacum, and Oe. dentireticulatum resembled strain FACHB-3311.

IR regions are the most conserved regions in the cp genomes. Frequent expansions and contractions at the junctions of SSC and LSC with IRs illustrate the relationships among the taxa and have been recognized as evolutionary signals [38-42]. The nine species of Oedogoniales examined in the present study showed only a few variations at the junctions. When compared with the two O. carolinianum, O. prescottii showed higher similarities to the five Oedogonium species, and the five Oedogonium species were similar to each other. The IR regions of O. prescottii and Oe. crispum presented a contraction, when compared with those of the other Oedogoniales taxa, and the cp genomes of both O. prescottii and Oe. crispum exhibited the shortest length. Previous studies have indicated that IR expansion and contraction frequently result in variations in genome size, which can be applied to phylogenetics and genome evolution analyses [39, 40, 43], and gene conversion during speciation is considered to be responsible for small IR expansions or contractions [38, 40, 44-46].

Phylogenetic studies based on 54 cp protein-coding genes assayed using ML and Bayesian analyses with amino acid and nucleotide datasets showed that Oedocladium and Oedogonium are polyphyletic, which is in accordance with that reported previously. However, the location of O. prescottii remained uncertain, and the support value based on nucleotide dataset was not very high at the basal node, probably owing to the lack of sufficient representative taxa for this group as well as different evolutionary rates of the amino acid sequence and nucleotide data. Previous studies have proposed that larger sample sizes can substantially improve the phylogenetic results [47].

Positively selected genes are known to play a key role in adaptation to different environments and speciation [48-52], and it is necessary to understand the adaptive evolutionary history of Oedocladium species. The results of the present study showed that 291SER of psbA may be under positive selection with posterior probability higher than 99%. The genus Oedocladium (terrestrial) is presumed to have partly originated from Oedogonium species, which grow on moist soil surface and present underground filaments with slightly unbranched rhizoids [9]. The psbA encodes the photosystem II reaction center protein D1, which is one of the two reaction center proteins of photosystem II. Photosystem II is the first link in the chain of photosynthesis, and captures photons and uses the energy to extract electrons from water molecules [53]. It has been reported that the genes in the cp genome (including psbA) of Curcuma sp. show adaptive evolution to adapt to the changes in light conditions [54], and that the green alga Chlamydomonas sp. ICE-L underwent adaptive evolution to adapt to extreme polar environment [26]. We speculate that the Oedocladium species and terrestrial Oedogonium species could have partly originated from the aquatic Oedogonium species, and might have undergone adaptive evolution during this process to adapt to the difference in light intensity between aquatic and terrestrial habitats. Nevertheless, more genomic data, especially for terrestrial species, may help to verify these hypotheses and further understand the phylogenetic and evolutionary relationships of the order Oedogoniales.

The present study determined the cp genomes of five Oedogonium speciesand revealed that the overall structure and gene contents of the Oedogoniales cp genomes were relatively conserved, except for some variations in genome sizes, AT contents, introns, and repeats. Phylogenetic analysis based on 54 cp protein-coding genes indicated that both Oedogonium and Oedocladium are polyphyletic. The positively selected gene in the two Oedocladium species was identified, and the terrestrial Oedogonium species were speculated to have undergone adaptive evolution to adapt to the difference in light intensity between aquatic and terrestrial habitats. These findings not only strengthen our understanding of Oedogoniales cp genomes, but also help us to comprehend the phylogenetic and evolutionary relationships of the order Oedogoniales.

Sampling, culture conditions, DNA extraction, and morphological observation

Species identification

The strains described in this study were isolated from water or soil samples, and have been deposited to the Freshwater Algae Culture Collection at the Institute of Hydrobiology (FACHB collection), Wuhan, Hubei Province, China. Strain FACHB-3309 was collected from a paddy field in Hechuan (29°50’15” N, 106°12’46” E), Chongqing Province, China, in March 2019. Strain FACHB-3310 was collected from a pond in Lvliang (37°34’20” N, 112°12’29.25” E), Shanxi Province, China, in July 2018. Strain FACHB-3311 was collected from a pond in Wuhan (30°3’46” N, 114°23’56” E), Hubei Province, China, in June 2019. Strain FACHB-3312 was collected from a ditch in Wuhan (30°33’2” N, 114°25’48” E), Hubei Province, China, in April 2019. Strain FACHB-3313 was collected from damp soil in a park in Haikou (20°2’23” N, 110°21’1” E), Hainan Province, China, in December 2018. All the strains were grown at 25°C in liquid BG11 medium under a 12/12-h light/dark cycle. An Olympus BX53 (Tokyo, Japan) light microscope equipped with an Olympus DP80 digital camera and CellSens standard image analysis software (Tokyo, Japan) were used for morphological examination.

Library preparation, sequencing, genome assembly, and annotation

A NEB Next Ultra DNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA) was used for preparing sequencing libraries, which were sequenced on an Illumina NovaSeq 6000 platform by a commercial provider (Novogene, Beijing, China). The methods of genome assembly and annotation have been described elsewhere [55]. The data were trimmed using SOAPnuke software [56] and assembled using SPAdes [57]. The resulting assembly contigs were considered to have originated from the cp genome if the (1) BLAST searches in publicly available cp genomes returned Chlorophyta species with significant e-values (1e-5); (2) GC content of the contigs was less than 45% (the GC content of previously sequenced green algal cp genomes is typically less than 45%); and (3) sequencing depth was more than 100-fold coverage. Subsequently, trimmed reads were aligned to the resulting assembly contigs using BWA-MEM [58]. If the reads mapped to two contigs, the order of the contigs was determined and one sequence was produced, which was confirmed by Sanger dideoxy sequencing. The tRNA genes were identified using tRNAscan-SE [59]. BLAST was used to refine the annotation results. Intron boundaries were determined by comparing intron-containing genes with homologs without introns, and intron subgroup affiliation was determined by modelling intron secondary structures [60, 61] using RNAweasel tool [62]. Forward and palindromic repeats larger than 30 bp were searched using Vmatch software (http://www.vmatch.de/) with the options -f -p -l -h -allmax and masked in the genome sequence by RepeatMasker (http://repeatmasker.org) running under the NCBI/RMBLAST (2.9.0+) search engine (http://blast.advbiocomp.com). The annotated sequences have been deposited to the NCBI GenBank database under the accession numbers MW250871–MW250875 (corresponding to strains FACHB-3309–FACHB-3313, respectively). Genome maps were generated using OrganellarGenomeDRAW [63].

Phylogenetic analysis

Phylogenetic analysis of the algal strains was performed by examining the sequences of cp protein-coding genes based on amino acid and nucleotide datasets. The amino acid and nucleotide datasets of the cp genomes were assembled using the following 54 protein-coding genes: atpA, atpB, atpE, atpF, atpH, atpI, cemA, chlB, chlL, chlN, clpP, petB, petD, petG, petL, psaA, psaB, psaC, psaJ, psbA, psbB, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ, rpl14, rpl16, rpl2, rpl20, rpl23, rpl36, rpl5, rps11, rps12, rps14, rps18, rps19, rps2, rps3, rps7, rps8, rps9, tufA, ycf12, ycf3, ycf4. The genes used in the amino acid dataset were aligned using MAFFT 7.0 [64], and those employed in the nucleotide dataset were additionally aligned using the MUSCLE function of MEGA7 [65] with the option “align codons” [66]. Ambiguous regions were removed from each alignment using trimAl 1.2 [67] with the option gt = 1. Evolutionary models and partitions of the datasets were determined using PartitionFinder 2 [68], and the best partitions are shown in Table 2. ML and Bayesian analyses were used for inferring phylogenies. IQ-TREE web server [69] was employed to perform ML analysis with 1000 ultrafast bootstraps [70] and 1000 SH-aLRT tests [70, 71] to examine nodal support. Bayesian analysis was conducted using MrBayesv3.2.6 [72], and the dataset was partitioned as shown in Table 2. Markov chain Monte Carlo analyses were run with four Markov chains (three heated, one cold) for 1,000,000 generations, and trees were sampled every 1000 generations. In each round of calculation, a fixed number of samples (burn-in = 1 000) was discarded from the beginning of the chain. An alignment of the cp genome sequences of all the species of Oedogoniales was generated using Mauve ver. 2.3.1 with the progressive mode [73]. FastANI [74] was employed to determine the ANI of all the cp genomes.

Evolutionary analysis

Evolutionary analysis based on 54 cp protein-coding genes was conducted using the CODEML program of PAML v4.9 [29]. Branch-site model was utilized to find genes that possibly underwent positive selection. The improved branch-site model (model = 2, Nsites = 2) was used to detect signatures of positive selection on individual codons in a specific branch [75]. The three Oedocladium species and the terrestrial Oedogonium sp. (strain FACHB-3313) were set as the foreground branch. The null model assumed that no positive selection occurred on the foreground branch (fix_omega = 1, omega = 1), and the alternative model assumed that sites on the foreground branch were under positive selection (fix_omega = 0, omega = 2). LRT were used to test model fit and Chi-square test was applied for examining the P values. A correction was performed for multiple testing using an FDR criterion, and BEB method was employed to statistically identify sites under positive selection. Genes with an FDR-adjusted P < 0.05 were considered as putatively selected.

OCC clade: Oedogoniales, chaetophorales, and chaetopeltidales; IR: Inverted repeats; ML: Maximum likelihood; cpDNA: Chloroplast DNA; Oe. : Oedogonium; O. : Oedocladium.

Acknowledgements

Special thanks go to Benwen Liu and Huan Zhu f or their helpful advice. The research was supported by the Wuhan Branch, Supercomputing Center, Chinese Academy of Sciences, China.

Authors' Contributions

Q. Xiong: original concept, culture experiments, data analysis, writing and editing manuscript; Y. Hu, W. Lv, H. Wang: analyzed and interpreted the data. All authors revised, read, and approved the final version of the manuscript.

Funding

This work was supported by Key Research Program of Frontier Sciences, CAS (Grant No. QYZDY-SSWSMC029-6) and the National Natural Science Foundation of China (grant no. 31750001).

Availability of data and materials

All data generated and analysed during this study are included in this published article and its supplementary information files. Raw sequencing data of all species are available from the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/). Accession numbers of cpDNA: MW250871-MW250875

MW250871: https://www.ncbi.nlm.nih.gov/nuccore/MW250871

MW250872: https://www.ncbi.nlm.nih.gov/nuccore/MW250872

MW250873: https://www.ncbi.nlm.nih.gov/nuccore/MW250873

MW250874: https://www.ncbi.nlm.nih.gov/nuccore/MW250874

MW250875: https://www.ncbi.nlm.nih.gov/nuccore/MW250875

Ethics approval and consent to participate

Not applicable.

Consent to publish

Not applicable.

Competing interests

The authors declare that they have no competing interests

Author details

¹Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China.

²University of Chinese Academy of Sciences, Beijing 100039, China.

³Yangtze River Basin Ecological Environment Monitoring and Scientific Research Center, Yangtze River Basin Ecological Environment Supervision and Administration Bureau, Ministry of Ecological Environment, Wuhan, China.

⁴Department of Biology, Taiyuan Normal University, Jinzhong 030619, PR China

⁵State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China.

Agardh, CA. Synopsis Algarum Scandinaviae, adjecta dispositione universali Algarum. Lundae; 1817.
De Bary A. Ueber die Algengattungen Oedogonium und Bulbochaete. Frankfurt: Abhandlung Senckenberg. Naturf. Ges; 1854;1:29-105.
Stahl, E. Oedocladium protonema, eine neue Oedogoniaceen-Gattung. Wiss. Bot; 1891;23: 339-348.
Hirn, KE. Monographie und Iconographie der Oedogoniaceen. Acta Soc Scienti Fennicae; 1900;27: 1-395.
Tiffany, LH. North American Flora: Oedogoniales. New York: Botanical Garden; 1937.
Gemeinhardt, K. Oedogoniales. In: Rabenhorst’s Kryptogamenflora von Deutschland und der Schweiz, 12. Leipzig. 1939.
Islam, NAKM, Sarma, P. Two new species of terrestrial Oedogonium from east Pakistan. Trans Amer Micros Soc. 1963;82: 74-77.
Gauthier-Lièvre, L. Oedogoniacées Africaines. Verlag von J. Cramer, Stuttgart; 1964.
Jao, CC. Monographia Oedogoniales Sinicae. Beijing: Science Press; 1979.
Mrozińska, T. Oedogoniophyceae: Oedogoniales. In: Süßwasserflora von Mitteleuropa 14;Chlorophyta VI. Stuttgart: Gustav Fischer Verlag; 1985.
van Den Hoek, C, Mann, DG, Jahns, UM. Algae: An Introduction to Phycology. Cambridge: Cambridge University Press; 1995.
Graham, LE, Wilcox, LW. Algae. Prentice Hall, Upper Saddle River, NJ; 2000.
Liu GX, Hu ZY: Predominant occurrence of apical cell divisions in Oedogoniampakistanense and its phylogenetic significance. Phycologia. 2004;43(6):669-671.
Booton GC, Floyd GL, Fuerst PA. Origins and affinities of the filamentous green algal orders Chaetophorales and Oedogoniales based on 18S rRNA gene sequences. J. Phycol. 1998;34(2):312-318.
Buchheim MA, Michalopulos EA, Buchheim JA. Phylogeny of the Chlorophyceae with special reference to the Sphaeropleales: A study of 18S and 26S rDNA data. J. Phycol. 2001;37(5):819-835.
Krienitz, L, Hegewald, E, Hepperle, D, Wolf, M. The systematics of coccoid green algae. 18S rRNA gene sequence data versus morphology. Biologia. 2003;58: 437-446.
Shoup S, Lewis LA.Polyphyletic origin of parallel basal bodies in swimming cells of Chlorophycean green algae (Chlorophyta). J. Phycol. 2003;39(4):789-796.
Alberghina, JS, Vigna, MS, Confalonieri, VA. Phylogenetic position of the Oedogoniales within the green algae (Chlorophyta) and the evolution of the absolute orientation of the flagellar apparatus. Plant Syst. Evol. 2006;261: 151-163.
Mei H, Luo W, Liu GX, Hu ZY. Phylogeny of Oedogoniales (Chlorophyceae, Chlorophyta) inferred from 18S rDNA sequences with emphasis on the relationships in the genus Oedogonium based on ITS-2 sequences. Plant Syst. Evol. 2007;265(3-4):179-191.
Ravi V, Khurana JP, Tyagi AK, Khurana P. An update on chloroplast genomes. Plant Syst. Evol. 2008;271(1-2):101-122.
Burke SV, Grennan CP, Duvall MR. Plastome sequences of two New World bamboos—Arundinaria gigantea and Cryptochloa strictiflora (Poaceae)—extend phylogenomic understanding of Bambusoideae. Am. J. Bot. 2012;99(12):1951-1961.
Dong W, Xu C, Cheng T, Zhou S. Complete Chloroplast Genome of Sedum sarmentosum and Chloroplast Genome Evolution in Saxifragales. PLoS ONE. 2013;8(10).
Huang H, Shi C, Liu Y, Mao S-Y, Gao L-Z. Thirteen Camellia chloroplast genome sequences determined by high-throughput sequencing: genome structure and phylogenetic relationships. BMC Evol. Biol. 2014;14.
Yi D-K, Lee H-L, Sun B-Y, Chung MY, Kim K-J. The complete chloroplast DNA sequence of (Araliaceae); Comparative evolutionary analyses with other three asterids. Mol Cells 2012;33(5):497-508.
Lemieux C, Otis C, Turmel M. Chloroplast phylogenomic analysis resolves deep-level relationships within the green algal class Trebouxiophyceae. BMC Evol. Biol. 2014;14.
Zhang Z, An M, Miao J, Gu Z, Liu C, Zhong B. The Antarctic sea ice alga Chlamydomonas sp. ICE-L provides insights into adaptive patterns of chloroplast evolution. BMC Plant Biol. 2018;18.
Brouard J-S, Otis C, Lemieux C, Turmel M. Chloroplast DNA sequence of the green alga Oedogonium cardiacum (Chlorophyceae): Unique genome architecture, derived characters shared with the Chaetophorales and novel genes acquired through horizontal transfer. BMC Genomics. 2008;9.
Brouard J-S, Turmel M, Otis C, Lemieux C. Proliferation of group II introns in the chloroplast genome of the green alga Oedocladium carolinianum (Chlorophyceae). Peerj. 2016;4.
Yang Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 2007;24(8):1586-1591.
Hu Y, Xing W, Song H, Zhu H, Liu G, Hu Z. Evolutionary Analysis of Unicellular Species in Chlamydomonadales Through Chloroplast Genome Comparison With the Colonial Volvocine Algae. Front Microbiol. 2019;10.
Smith DR. Mutation Rates in Plastid Genomes. They Are Lower than You Might Think. Genome Biol Evol. 2015;7(5):1227-1234.
Iram S, Hayat MQ, Tahir M, Gul A, Abdullah, Ahmed I. Chloroplast Genome Sequence of Artemisia scoparia: Comparative Analyses and Screening of Mutational Hotspots. Plants. 2019;8(11).
Henriquez CL, Abdullah, Ahmed I, Carlsen MM, Zuluaga A, Croat TB, McKain MR. Molecular evolution of chloroplast genomes in Monsteroideae (Araceae). Planta. 2020;251(3).
Besendahl A, Qiu YL, Lee JH, Palmer JD, Bhattacharya D. The cyanobacterial origin and vertical transmission of the plastid tRNA(Leu) group-I intron. Curr. Genet. 2000;37(1):12-23.
Xiong Q, Hu Y, Liu B, Zhu H, Liu G, Hu Z. Chloroplast genomes and phylogenetic analysis of two species of Oedocladium (Oedogoniales, Chlorophyta). Eur. J. Phycol. 2021.
Watanabe S, Fucikova K, Lewis LA, Lewis PO. Hiding in plain sight: Koshicola spirodelophila gen. et sp nov (Chaetopeltidales, Chlorophyceae), a novel green alga associated with the aquatic angiosperm Spirodela polyrhiza. Am. J. Bot. 2016;103(5):865-875.
Khakhlova O, Bock R. Elimination of deleterious mutations in plastid genomes by gene conversion. Plant J. 2006;46(1):85-94.
Raubeson LA, Peery R, Chumley TW, Dziubek C, Fourcade HM, Boore JL, Jansen RK. Comparative chloroplast genomics: analyses including new sequences from the angiosperms Nuphar advena and Ranunculus macranthus. BMC Genomics. 2007;8.
Wang R-J, Cheng C-L, Chang C-C, Wu C-L, Su T-M, Chaw S-M. Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots. BMC Evol. Biol. 2008;8.
Park I, Yang S, Kim WJ, Noh P, Lee HO, Moon BC. The Complete Chloroplast Genomes of Six Ipomoea Species and Indel Marker Development for the Discrimination of Authentic Pharbitidis Semen (Seeds of I-nil or I-purpurea). Front Plant Sci. 2018;9.
Lu Q, Ye W, Lu R, Xu W, Qiu Y. Phylogenomic and Comparative Analyses of Complete Plastomes of Croomia and Stemona (Stemonaceae). Int J Mol Sci. 2018;19(8).
Kim KJ, Lee HL. Complete chloroplast genome sequences from Korean ginseng (Panax schinseng Nees) and comparative analysis of sequence evolution among 17 vascular plants. Dna Res. 2004;11:247-261
Goulding SE, Olmstead RG, Morden CW, Wolfe KH. Ebb and flow of the chloroplast inverted repeat. Mol Gen Genet. 1996;252(1-2):195-206.
Meng X-X, Xian Y-F, Xiang L, Zhang D, Shi Y-H, Wu M-L, Dong G-Q, Ip S-P, Lin Z-X, Wu L et al. Complete Chloroplast Genomes from Sanguisorba: Identity and Variation Among Four Species. Molecules. 2018;23(9).
Choi I-S, Jansen R, Ruhlman T. Lost and Found. Return of the Inverted Repeat in the Legume Clade Defined by Its Absence. Genome Biol Evol. 2019;11(4):1321-1333.
Pollock DD, Zwickl DJ, McGuire JA, Hillis DM. Increased taxon sampling is advantageous for phylogenetic inference. Syst. Biol. 2002;51(4):664-671.
Wang L, Wuyun T-n, Du H, Wang D, Cao D. Complete chloroplast genome sequences of Eucommia ulmoides: genome structure and evolution. Tree Genet. Genom. 2016;12(1).
Ma Q, Li S, Bi C, Hao Z, Sun C, Ye N. Complete chloroplast genome sequence of a major economic species, Ziziphus jujuba (Rhamnaceae). Curr. Genet. 2017;63(1):117-129.
Fan W-B, Wu Y, Yang J, Shahzad K, Li Z-H. Comparative Chloroplast Genomics of Dipsacales Species: Insights Into Sequence Variation, Adaptive Evolution, and Phylogenetic Relationships. Front Plant Sci. 2018;9.
Gao C, Deng Y, Wang J. The Complete Chloroplast Genomes of Echinacanthus Species (Acanthaceae): Phylogenetic Relationships, Adaptive Evolution, and Screening of Molecular Markers. Front Plant Sci. 2019;9.
Wu Y, Liu F, Yang D-G, Li W, Zhou X-J, Pei X-Y, Liu Y-G, He K-L, Zhang W-S, Ren Z-Y et al. Comparative Chloroplast Genomics of Gossypium Species: Insights Into Repeat Sequence Variations and Phylogeny. Front Plant Sci. 2018;9.
Ferreira KN, Iverson TM, Maghlaoui K, Barber J, Iwata S. Architecture of the photosynthetic oxygen-evolving center. Science. 2004;303(5665):1831-1838.
Gui L, Jiang S, Xie D, Yu L, Huang Y, Zhang Z, Liu Y. Analysis of complete chloroplast genomes of Curcuma and the contribution to phylogeny and adaptive evolution. Gene. 2020;732.
Hu Y, Xing W, Song H, Liu G, Hu Z. Analysis of mitochondrial and chloroplast genomes in two volvocine algae: Eudorina elegans and Eudorina cylindrica (Volvocaceae, Chlorophyta). Eur. J. Phycol. 2019;54(2):193-205.
Chen Y, Chen Y, Shi C, Huang Z, Zhang Y, Li S, Li Y, Ye J, Yu C, Li Z et al. SOAPnuke. a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience. 2017;7(1).
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Son P, Prjibelski AD et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. J. Comput. Biol. 2012;19(5):455-477.
Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv preprint, arXiv. 2013;1-3.
Lowe TM, Chan PP. tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016;44(W1):W54-W57.
Michel F, Umesono K, Ozeki H. Comparative and functional anatomy of group II catalytic introns--a review. Gene. 1989;82(1):5-30.
Michel F, Westhof E. Modelling of the three-dimensional architecture of group I catalytic introns based on comparative sequence analysis. J. Mol. Biol. 1990;216(3):585-610.
Lang BF, Laforest M-J, Burger G. Mitochondrial introns: a critical view. Trends Genet. 2007;23(3):119-125.
Lohse M, Drechsel O, Bock R. OrganellarGenomeDRAW (OGDRAW): a tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr. Genet. 2007;52(5-6):267-274.
Katoh K, Standley DM. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013;30(4):772-780.
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016;33(7):1870-1874.
Edgar, RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792-1797.
Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25(15):1972-1973.
Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: New Methods for Selecting Partitioned Models of Evolution for Molecular and Morphological Phylogenetic Analyses. Mol. Biol. Evol. 2017;34(3):772-773.
Trifinopoulos J, Lam-Tung N, von Haeseler A, Minh BQ. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016;44(W1):W232-W235.
Bui Quang M, Minh Anh Thi N, von Haeseler A. Ultrafast Approximation for Phylogenetic Bootstrap. Mol. Biol. Evol. 2013;30(5):1188-1195.
Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst. Biol. 2010;59(3):307-321.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space. Syst. Biol. 2012;61(3):539-542.
Darling AE, Mau B, Perna NT. progressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement. PLoS ONE. 2010;5(6).
Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun. 2018;9.
Zhang JZ, Nielsen R, Yang ZH. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol. Biol. Evol. 2005;22(12):2472-2479.
Lemieux C, Otis C, Turmel M. Ancestral chloroplast genome in Mesostigmaviride reveals an early branch of green plant evolution. Nature. 2000;Suppl 403:649-652

Table 1. General features of nine oedogonialean chloroplast genomes.

Genomic Feature	O. prescottii	O. carolinianum (FACHB-2453)	O. carolinianum	Oe. cardiacum	Oe. dentireticulatum	Oe. crispum	Oe.sp. (FACHB-3311)	Oe. capilliforme	Oe.sp. (FACHB-3313)
Size (bp)
Total	154,9 78	200,832	204,438	196,547	159,341	146,367	187,104	195,349	179,946
IR	12,808	22,275	23,748	35,492	19,159	13,284	25,841	35,138	26,999
LSC	80,821	98,462	98,887	80,363	77,718	76,475	86,403	80,003	80,347
SSC	48,542	57,820	58,055	45,200	43,305	43,324	49,019	45,070	45,601
A+T(%)	72.66	70.51	70.2	70.5	71.46	72.28	69.98	70.64	71.23
Coding proportion^a	67.0%	52.7%	51.4%	55.9%	64.1%	69.5%	54.8%	52.8%	55.9%
Gene (+/-)^bc	49/60	50/61	50/61	52/63	50/59	50/61	49/60	49/60	50/61
Protein-coding genes(number/proportion)	68 29/43	68 28/44	68 28/44	70 31/45	68 29/43	68 29/43	68 29/43	68 29/43	68 29/43
rRNA(number/proportion)	3	3	3	3	3/3	3/3	3/3	3/3	3/3
tRNA(number/proportion)	28 17/14	30 19/14	30 19/14	29 18/15	29 18/14	29 18/15	28 17/14	28 17/14	29 18/15
Introns	6	15	17	21	22	8	37	24	30
Group I (no.)	1	5	7	17	18	4	26	20	25
Group II (no.)	5	10	10	4	4	4	11	4	5
Repeats^d(%)	4.8	8.9	11.3	1.3	3.9	4.0	4.7	3.3	4.9
Accession number	MT364368	MT364369	NC_031510	EU677193	MW250871	MW250872	MW250873	MW250874	MW250875

^aThe coding proportion only includes all annotated protein-, rRNA-, and tRNA-coding regions; ^bGene and CDS numbers do not include ORF genes; ^cThe plus-minus sign means number of genes in plus strain (left side of slash) or minus strain (rightside of slash). ^dNon-overlapping repeat elements were mapped on each genome with RepeatMasker using as input sequences the repeats ≥30 bp identified with Vmatch.

Table 2. Partition scheme of 54 concatenated chloroplast protein-coding genes used in this study.

	Amino acid data sets		Nucleotide data sets
Subset	Best model	Partition scheme	Best model	Partition scheme

1	LG+I+G	rps12, chlL, atpB, atpA	GTR+I+G	atpA, rps19, tufA, atpF, atpI, ycf12, petL
2	LG+G	atpF, chlN, atpE	GTR+I+G	psaJ, psbZ, psbB, psaB, petD, psbD, psbK, petB, ycf3, psbH, rpl16, atpB, rpl14
3	MTZOA+G+F	psbM, petB, psbI, psbE, psbT, psbN, atpH, psbD, psbB, psaA, psaB	GTR+I+G	ycf4, rps12, rpl2, atpE
4	CPREV+G	rps19, atpI, chlB, rpl16, rpl2, rpl5	GTR+G	atpH, psbA
5	CPREV+G	rps14, rpl20, rps9, cemA, ycf4	GTR+G	psbM, psbI, psaC, psbL, psbF, psbE, cemA, rpl36, rps14
6	JTT+G+F	rps18, rps8, rpl23, rps3, clpP	GTR+I+G	psbN, chlB, chlL, rps11, rpl5, chlN
7	MTZOA+G	psaJ, psbF, psbJ, psbK, psbH, petL, psbZ, petD	GTR+G	clpP, rps7
8	LG	petG, ycf3	GTR+G	psbJ, psbT, petG
9	PMB+G	rpl36, psbA, psaC	GTR+I+G	psaA
10	LG	rps11, psbL, rpl14	GTR+G	rps2, rps8, rps9, rpl20, rpl23, rps18
11	LG+G	rps2	GTR+G	rps3
12	FLU+G	rps7
13	LG4M+G	tufA
14	MTZOA	ycf12

No competing interests reported.

SupplementaryfigureS1.GenemapofOedogoniumdentireticulatumchloroplastgenome.pdf
Gene map of Oe. dentireticulatum chloroplast genomes.
SupplementaryfigureS2.GenemapofOedogoniumcrispumchloroplastgenome.pdf
Gene map of Oe. crispum chloroplast genomes.
SupplementaryfigureS3.GenemapofOedogniumsp.FACHB3311chloroplastgenome.pdf
Gene map of Oe. sp. (FACHB-3311) chloroplast genomes.
SupplementaryfigureS4.GenemapofOe.capilliformechloroplastgenomes.pdf
Gene map of Oe. capilliforme chloroplast genomes.
SupplementaryfigureS5.GenemapofOe.sp.FACHB3313chloroplastgenomes.pdf
Gene map of Oe. sp. (FACHB-3313) chloroplast genomes. Arrows show the direction of transcription. The same colour block shows the functional gene group (legend at bottom left). Transfer RNAs are represented by their one-letter amino acid code. The grey circle on the inside shows a graph of the GC content.
SupplementaryfigureS6.HeatmapofANIvaluesofthenineOedogonialescpgenomes.pdf
Heat map of ANI values of the nine Oedogoniales cp genomes.
SupplementarytableS1.InsertionsitesofgroupIintronsofthenineOedogoniaescpgenomes..docx
SupplementarytableS2.InsertionsitesofgroupIIintronsofthenineOedogoniaescpgenomes..docx

Download PDF

Editorial decision: Major revision
20 May, 2021
Reviews received at journal
16 May, 2021
Reviewers agreed at journal
12 May, 2021
Reviewers agreed at journal
12 May, 2021
Reviewers invited by journal
12 May, 2021
Editor assigned by journal
12 May, 2021
Editor invited by journal
12 May, 2021
Submission checks completed at journal
12 May, 2021
First submitted to journal
22 Apr, 2021

You are reading this latest preprint version

Chloroplast Genomes of Five Oedogonium Species: Genome Structure, Phylogenetic Analysis and Adaptive Evolution

Status:

Version 1

Abstract

Figures

Background

Results

Discussion

Conclusion

Methods

Abbreviations

Declarations

References

Tables

Additional Declarations

Supplementary Files

Status:

Version 1