Definition of a High-Resolution Molecular Marker for Tracking the Genetic Diversity of the Harmful Algal Species Eucampia zodiacus Through Comparative Analysis of Mitochondrial Genomes

Genetic Diversity and Geographical Distribution of the Red Tide Species Coscinodiscus granii Revealed Using a High-Resolution Molecular Marker

Diatoms are responsible for approximately 40% of the global primary photosynthetic production and account for up to 20% of global carbon fixation. Coscinodiscus granii is a red tide forming species of the phylum Bacillariophyta that has been detected in a wide range of coastal regions, suggesting the possibility of the existence of high genetic diversity with differential adaptation. Common molecular markers including 18S rDNA, 16S rDNA, ITS, cox1, and rbcL do not provide sufficient resolution for distinguishing intra-species genetic diversity, hindering in-depth research on intra-species genetic diversity and their spatial and temporal dynamics. In this project, we aimed to develop molecular markers with high resolution and specificity for C. granii, attempting to identify different taxa of this species, which will set up a stage for subsequent functional assays. Comparative genomics analysis of the mtDNAs of C. granii strains identified a genomic region with high genomic variations, which was used to guide the development of a molecular marker with high resolution and high specificity. This new molecular marker, which was named cgmt1 (C. granii mitochondrial 1), was 376 bp in size and differentiated C. granii samples collected in coastal regions of China into three different clades. Preliminary analysis of field samples collected in various coastal regions in China revealed that C. granii clades were almost exclusively found in the Bohai Sea and the north Yellow Sea. This newly developed molecular marker cgmt1 could be used for tracking intra-species genetic diversity and biogeographic distribution of C. granii in different ecosystems.

The complete mitochondrial genome of Isochrysis galbana harbors a unique repeat structure and a specific trans-spliced cox1 gene

The haptophyte Isochrysis galbana is considered as a promising source for food supplements due to its rich fucoxanthin and polyunsaturated fatty acids content. Here, the I. galbana mitochondrial genome (mitogenome) was sequenced using a combination of Illumina and PacBio sequencing platforms. This 39,258 bp circular mitogenome has a total of 46 genes, including 20 protein-coding genes, 24 tRNA genes and two rRNA genes. A large block of repeats (~12.7 kb) was segregated in one region of the mitogenome, accounting for almost one third of the total size. A trans-spliced gene cox1 was first identified in I. galbana mitogenome and was verified by RNA-seq and DNA-seq data. The massive expansion of tandem repeat size and cis- to trans-splicing shift could be explained by the high mitogenome rearrangement rates in haptophytes. Strict SNP calling based on deep transcriptome sequencing data suggested the lack of RNA editing in both organelles in this species, consistent with previous studies in other algal lineages. To gain insight into haptophyte mitogenome evolution, a comparative analysis of mitogenomes within haptophytes and among eight main algal lineages was performed. A core gene set of 15 energy and metabolism genes is present in haptophyte mitogenomes, consisting of 1 cob, 3 cox, 7 nad, 2 atp and 2 ribosomal genes. Gene content and order was poorly conserved in this lineage. Haptophyte mitogenomes have lost many functional genes found in many other eukaryotes including rps/rpl, sdh, tat, secY genes, which make it contain the smallest gene set among all algal taxa. All these implied the rapid-evolving and more recently evolved mitogenomes of haptophytes compared to other algal lineages. The phylogenetic tree constructed by cox1 genes of 204 algal mitogenomes yielded well-resolved internal relationships, providing new evidence for red-lineages that contained plastids of red algal secondary endosymbiotic origin. This newly assembled mitogenome will add to our knowledge of general trends in algal mitogenome evolution within haptophytes and among different algal taxa.

Diatom Biodiversity and Speciation Revealed by Comparative Analysis of Mitochondrial Genomes

Diatoms (Bacillariophyta) constitute one of the most diverse and ecologically significant groups of phytoplankton, comprising 100,000-200,000 species in three classes Bacillariophyceae, Mediophyceae, and Coscinodiscophyceae. However, due to the limited resolution of common molecular markers including 18S rDNA, 28S rDNA, ITS, rbcL, and cox1, diatom biodiversity has not been adequately ascertained. Organelle genomes including mitochondrial genomes (mtDNAs) have been proposed to be "super barcodes" for distinguishing diatom species because of their rich genomic content, and the rapid progress of DNA sequencing technologies that has made it possible to construct mtDNAs with increasing throughout and decreasing cost. Here, we constructed complete mtDNAs of 15 diatom species including five Coscinodiscophyceae species (Guinardia delicatula, Guinardia striata, Stephanopyxis turris, Paralia sulcata, and Actinocyclus sp.), four Mediophyceae species (Hemiaulus sinensis, Odontella aurita var. minima, Lithodesmioides sp., and Helicotheca tamesis), and six Bacillariophyceae species (Nitzschia ovalis, Nitzschia sp., Nitzschia traheaformis, Cylindrotheca closterium, Haslea tsukamotoi, and Pleurosigma sp.) to test the practicality of using mtDNAs as super barcodes. We found that mtDNAs have much higher resolution compared to common molecular markers as expected. Comparative analysis of mtDNAs also suggested that mtDNAs are valuable in evolutionary studies by revealing extensive genome rearrangement events with gene duplications, gene losses, and gains and losses of introns. Synteny analyses of mtDNAs uncovered high conservation among species within an order, but extensive rearrangements including translocations and/or inversions between species of different orders within a class. Duplication of cox1 was discovered for the first time in diatoms in Nitzschia traheaformis and Haslea tsukamotoi. Molecular dating analysis revealed that the three diatom classes split 100 Mya and many diatom species appeared since 50 Mya. In conclusion, more diatom mtDNAs representing different orders will play great dividends to explore biodiversity and speciation of diatoms in different ecological regions.

Complete chloroplast genome of Eucampia zodiacus (Mediophyceae, Bacillariophyta)

The cosmopolitan phytoplankton species Eucampia zodiacus Ehrenberg 1839 is a common harmful algal bloom (HAB) species with significant negative ecological impact. However, molecular information for this HAB species is limited. In this study, the complete chloroplast genome (cpDNA) of E. zodiacus was constructed for the first time. The circular genome was 118,107 bp in length, containing a pair of inverted repeats (IR) (6991 bp each). The overall GC content of E. zodiacus cpDNA was 31.64%. It encoded 169 genes, including 131 protein-coding genes (PCGs), 30 tRNA genes, one tmRNA gene, one ncRNA gene, and six rRNA genes in IR regions. Phylogenetic analysis using concatenated PCGs of 56 diatom cpDNAs strongly supported that E. zodiacus was closely related to Cerataulina daemon, which belongs to the same family Hemiaulaceae according to AlgaeBase. Syntenic analysis revealed nearly identical gene order between the two cpDNAs, except for an inversion in the small single-copy (SSC) region.

The complete mitochondrial genome and phylogenetic analysis of Coscinodiscus wailesii (Coscinodiscophyceae, Bacillariophyta)

Coscinodiscus is a species-rich genus with about 400 species described, many of which are harmful algal bloom species with significant negative ecological impact. Despite of their importance in primary production and as harmful algal bloom species, genome data for species in this genus is limited. No mitochondrial genome (mtDNA) of any species in this genus has been reported. Here, the complete mtDNA sequence of the Coscinodiscus wailesii Gran & Angst 1931 was constructed and analyzed. The circular mtDNA was 36,071 bp in length, encoding 64 genes, including 34 protein coding genes (PCGs), 24 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and 4 conserved open reading frames (orfs). The overall AT content of C. wailesii mtDNA was 75.00%, which was slightly lower than that of Melosira undulate (78.40%). Maximum likelihood (ML) phylogenetic analysis using 29 shared protein-coding genes revealed that C. wailesii clustered well with M. undulata, which was the only species of class Coscinodiscophyceae whose mtDNA has been fully constructed. The complete mtDNAs of more Coscinodiscus species will be valuable for studying the evolutionary relationships among species in the genus Coscinodiscus and in the Class of Coscinodiscophyceae.