Complete mitochondrial genome of Skeletonema marinoi (Mediophyceae, Bacillariophyta), a clonal chain forming diatom in the west coast of Korea

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.

Identification of Medicinal Bidens Plants for Quality Control Based on Organelle Genomes

Bidens plants are annuals or perennials of Asteraceae and usually used as medicinal materials in China. They are difficult to identify by using traditional identification methods because they have similar morphologies and chemical components. Universal DNA barcodes also cannot identify Bidens species effectively. This situation seriously hinders the development of medicinal Bidens plants. Therefore, developing an accurate and effective method for identifying medicinal Bidens plants is urgently needed. The present study aims to use phylogenomic approaches based on organelle genomes to address the confusing relationships of medicinal Bidens plants. Illumina sequencing was used to sequence 12 chloroplast and eight mitochondrial genomes of five species and one variety of Bidens. The complete organelle genomes were assembled, annotated and analysed. Phylogenetic trees were constructed on the basis of the organelle genomes and highly variable regions. The organelle genomes of these Bidens species had a conserved gene content and codon usage. The 12 chloroplast genomes of the Bidens species were 150,489 bp to 151,635 bp in length. The lengths of the eight mitochondrial genomes varied from each other. Bioinformatics analysis revealed the presence of 50–71 simple sequence repeats and 46–181 long repeats in the organelle genomes. By combining the results of mVISTA and nucleotide diversity analyses, seven candidate highly variable regions in the chloroplast genomes were screened for species identification and relationship studies. Comparison with the complete mitochondrial genomes and common protein-coding genes shared by each organelle genome revealed that the complete chloroplast genomes had the highest discriminatory power for Bidens species and thus could be used as a super barcode to authenticate Bidens species accurately. In addition, the screened highly variable region trnS-GGA-rps4 could be also used as a potential specific barcode to identify Bidens species.

Comparative analysis of full-length mitochondrial genomes of five Skeletonema species reveals conserved genome organization and recent speciation

Background

Skeletonema species are prominent primary producers, some of which can also cause massive harmful algal blooms (HABs) in coastal waters under specific environmental conditions. Nevertheless, genomic information of Skeletonema species is currently limited, hindering advanced research on their role as primary producers and as HAB species. Mitochondrial genome (mtDNA) has been extensively used as “super barcode” in the phylogenetic analyses and comparative genomic analyses. However, of the 21 accepted Skeletonema species, full-length mtDNAs are currently available only for a single species, S. marinoi.

Results

In this study, we constructed full-length mtDNAs for six strains of five Skeletonema species, including S. marinoi, S. tropicum, S. grevillei, S. pseudocostatum and S. costatum (with two strains), which were isolated from coastal waters in China. The mtDNAs of all of these Skeletonema species were compact with short intergenic regions, no introns, and no repeat regions. Comparative analyses of these Skeletonema mtDNAs revealed high conservation, with a few discrete regions of high variations, some of which could be used as molecular markers for distinguishing Skeletonema species and for tracking the biogeographic distribution of these species with high resolution and specificity. We estimated divergence times among these Skeletonema species using 34 mtDNAs genes with fossil data as calibration point in PAML, which revealed that the Skeletonema species formed the independent clade diverging from Thalassiosira species approximately 48.30 Mya.

Conclusions

The availability of mtDNAs of five Skeletonema species provided valuable reference sequences for further evolutionary studies including speciation time estimation and comparative genomic analysis among diatom species. Divergent regions could be used as molecular markers for tracking different Skeletonema species in the fields of coastal regions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07999-z.

Complete mitochondrial genome of Rhizosolenia setigera (Coscinodiscophyceae, Bacillariophyta)

Rhizosolenia is a species-rich genus with 144 described species, many of which are harmful algal species (HABs) with significant negative ecological impact. Despite their significance in primary production and their potential to induce HABs, genome data of these species remain extremely limited. In this study, the complete mitochondrial genome (mtDNA) of Rhizosolenia setigera Brightwell 1858 was determined for the first time, which also represented the first mtDNA of the order Rhizosoleniales. The circular mtDNA was 34,792 bp in length with GC content of 23.28%. It encoded 63 genes including 35 protein-coding genes (PCGs), 24 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 2 conserved open reading frames (orfs). Phylogenetic analysis using concatenated PCGs revealed that R. setigera and Melosira undulate, which also belongs to the class Coscinodiscophyceae, clustered together as expected. However, comparison of these two mtDNAs revealed extensive genome rearrangement events, suggesting large evolutionary distance. The complete mtDNA of R. setigera will facilitate research on the phylogenetic relationship among Rhizosolenia species, which will in turn facilitate exploration of the evolutionary relationships in the class of Coscinodiscophyceae.

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

The cosmopolitan phytoplankton species Eucampia zodiacus is a common harmful algal bloom (HAB) species that have been found to cause HABs in essentially all coastal regions except the Polar regions. However, molecular information for this HAB species is limited with only a few molecular markers. In this project, we constructed the mitochondrial genome (mtDNA) of E. zodiacus, which was also the first mtDNA constructed for any species in the order Hemiaulales that includes 145 reported species (including two additional HAB species Cerataulina bicornis and Cerataulina pelagica). Comparative analysis of eight E. zodiacus strains revealed that they could not be distinguished using common molecular markers, suggesting that common molecular markers do not have adequate resolution for distinguishing E. zodiacus strains. However, these E. zodiacus strains could be distinguished using whole mtDNAs, suggesting the presence of different genotypes due to evolutionary divergence. Through comparative analysis of the mtDNAs of multiple E. zodiacus strains, we identified a new molecular marker ezmt1 that could adequately distinguish different E. zodiacus strains isolated in various coastal regions in China. This molecular marker ezmt1, which was ∼400 bp in size, could be applied to identify causative genotypes during E. zodiacus HABs through tracking the dynamic changes of genetic diversity of E. zodiacus in HABs.

The complete chloroplast and mitochondrial genomes of the diatom Nitzschia palea (Bacillariophyceae) demonstrate high sequence similarity to the endosymbiont organelles of the dinotom Durinskia baltica

Nitzschia palea is a common freshwater diatom used as a bioindicator because of its tolerance of polluted waterways. There is also evidence it may be the tertiary endosymbiont within the "dinotom" dinoflagellate Durinskia baltica. A putative strain of N. palea was collected from a pond on the University of Virginia's College at Wise campus and cultured. For initial identification, three markers were sequenced-nuclear 18S rDNA, the chloroplast 23S rDNA, and rbcL. Morphological characteristics were determined using light and scanning electron microscopy; based on these observations the cells were identified as N. palea and named strain "Wise." DNA from N. palea was deep sequenced and the chloroplast and mitochondrial genomes assembled. Single gene phylogenies grouped N. palea-Wise within a clearly defined N. palea clade and showed it was most closely related to the strain "SpainA3." The chloroplast genome of N. palea is 119,447 bp with a quadripartite structure, 135 protein-coding, 28 tRNA, and 3 rRNA genes. The mitochondrial genome is 37,754 bp with a single repeat region as found in other diatom chondriomes, 37 protein-coding, 23 tRNA, and 2 rRNA genes. The chloroplast genomes of N. palea and D. baltica have identical gene content, synteny, and a 92.7% pair-wise sequence similarity with most differences occurring in intergenic regions. The N. palea mitochondrial genome and D. baltica's endosymbiont mitochondrial genome also have identical gene content and order with a sequence similarity of 90.7%. Genome-based phylogenies demonstrated that D. baltica is more similar to N. palea than any other diatom sequence currently available. These data provide the genome sequences of two organelles for a widespread diatom and show they are very similar to those of Durinskia baltica's endosymbiont.

Zebrafish-based identification of the antiseizure nucleoside inosine from the marine diatom Skeletonema marinoi

With the goal of identifying neuroactive secondary metabolites from microalgae, a microscale in vivo zebrafish bioassay for antiseizure activity was used to evaluate bioactivities of the diatom Skeletonema marinoi, which was recently revealed as being a promising source of drug-like small molecules. A freeze-dried culture of S. marinoi was extracted by solvents with increasing polarities (hexane, dichloromethane, methanol and water) and these extracts were screened for anticonvulsant activity using a larval zebrafish epilepsy model with seizures induced by the GABAA antagonist pentylenetetrazole. The methanolic extract of S. marinoi exhibited significant anticonvulsant activity and was chosen for bioassay-guided fractionation, which associated the bioactivity with minor constituents. The key anticonvulsant constituent was identified as the nucleoside inosine, a well-known adenosine receptor agonist with previously reported antiseizure activities in mice and rat epilepsy models, but not reported to date as a bioactive constituent of microalgae. In addition, a UHPLC-HRMS metabolite profiling was used for dereplication of the other constituents of S. marinoi. Structures of the isolated compounds were elucidated by nuclear magnetic resonance and high-resolution spectrometry. These results highlight the potential of zebrafish-based screening and bioassay-guided fractionation to identify neuroactive marine natural products.

Comparative analysis of the mitochondrial genomes of six newly sequenced diatoms reveals group II introns in the barcoding region of cox1

Diatoms are the most diverse lineage of algae and at the base of most aquatic food webs, but only 11 of their mitochondrial genomes have been described. Herein, we present the mitochondrial genomes of six diatom species, including: Melosira undulata, Nitzschia alba, Surirella sp., Entomoneis sp., Halamphora coffeaeformis, and Halamphora calidilacuna. Comparison of these six genomes to the 11 currently published diatom mitochondrial genomes revealed a novel ubiquitous feature block consisting of tatC-orf157-rps11. The presence of intronic retrotransposable elements in the barcoding region of cox1 in the Halamphora genomes may explain historic difficulty (especially PCR) with cox1 as a universal barcode for diatoms. Our analysis suggests that high rates of variability in number and position of introns, in many commonly used coding sequences, prevent these from being universally viable as barcodes for diatoms. Therefore, we suggest researchers examine the chloroplast and/or nuclear genomes for universal barcoding markers.

Paleolimnology and resurrection ecology: The future of reconstructing the past

Paleolimnologists have utilized lake sediment records to understand historical lake and landscape development, timing and magnitude of environmental change at lake, watershed, regional and global scales, and as historical datasets to target watershed and lake management. Resurrection ecologists have long recognized lake sediments as sources of viable propagules ("seed or egg banks") with which to explore questions of community ecology, ecological response, and evolutionary ecology. Most researchers consider Daphnia as the primary model organism in these efforts, but many other aquatic biota, from viruses to macrophytes, similarly produce viable propagules that are incorporated in the sediment record but have been underutilized in resurrection ecology. The common goals shared by these two disciplines have led to mutualistic and synergistic collaborations-a development that must be encouraged to expand. We give an overview of the achievements of paleolimnology and the reconstruction of environmental history of lakes, review the untapped diversity of aquatic organisms that produce dormant propagules, compare Daphnia as a model of resurrection ecology with other organisms amenable to resurrection studies, especially diatoms, and consider new research directions that represent the nexus of these two fields.

Complete mitochondrial genome sequence of the freshwater diatom Asterionella formosa

We report the complete mitochondrial genome sequence of the freshwater diatom Asterionella formosa. The large 61.9 kb circular sequence encodes 34 proteins and 25 tRNAs that are universally conserved in other sequenced diatoms. We fully resolved a unique 24 kb region containing highly conserved repeated sequence units, possibly collocating with an origin of replication.

Complete mitochondrial genome of biraphid benthic diatom, Navicula ramosissima (Naviculales, Bacillariophyceae)

The complete mitochondrial DNA of biraphid benthic diatom, Navicula ramosissima TA439 was sequenced and characterized. The circular mitogenome contains 67 genes in 48,652 bp (31.1% GC), including 41 protein-coding, 24 transfer RNA (tRNA) and 2 rRNA genes. Twenty-four protein-coding sequences (CDS, 59%) have start with ATG codon and 17 CDS start with alternatives such as ATA (5), ATT (6), TTA (5) and TTG (1). The GC content of tRNA genes (42.1%) is relatively higher than those of the rRNA (35.2%) and CDS (30.5%). Three genes are consisted of multiple exons and introns, i.e. cox1 (three exons, two introns), rps11 (two exons, one intron), rrl (four exons, three introns). Phylogeny of diatoms based on mitogenome data (34 CDS, 8530 amino acids combined) supports the monophyly of Naviculales, including N. ramosissima (Naviculaceae), Berkeleya fennica (Berkeleyaceae), Fistulifera solaris (Stauroneidaceae) and Phaeodactylum tricornutum (Phaeodactylaceae). Mitogenome data may be useful for phylogenetic study of the diatoms and stramenopiles.