Plastid Genome of Dictyopteris divaricata (Dictyotales, Phaeophyceae): Understanding the Evolution of Plastid Genomes in Brown Algae

Comparative structure and evolution of the organellar genomes of Padina usoehtunii (Dictyotales) with the brown algal crown radiation clade

BACKGROUND

Organellar genomes have become increasingly essential for studying genetic diversity, phylogenetics, and evolutionary histories of seaweeds. The order Dictyotales (Dictyotophycidae), a highly diverse lineage within the Phaeophyceae, is long-term characterized by a scarcity of organellar genome datasets compared to orders of the brown algal crown radiation (Fucophycidae).

RESULTS

We sequenced the organellar genomes of Padina usoehtunii, a representative of the order Dictyotales, to investigate the structural and evolutionary differences by comparing to five other major brown algal orders. Our results confirmed previously reported findings that the rate of structural rearrangements in chloroplast genomes is higher than that in mitochondria, whereas mitochondrial sequences exhibited a higher substitution rate compared to chloroplasts. Such evolutionary patterns contrast with land plants and green algae. The expansion and contraction of the inverted repeat (IR) region in the chloroplast correlated with the changes in the number of boundary genes. Specifically, the size of the IR region influenced the position of the boundary gene rpl21, with complete rpl21 genes found within the IR region in Dictyotales, Sphacelariales and Ectocarpales, while the rpl21 genes in Desmarestiales, Fucales, and Laminariales span both the IR and short single copy (SSC) regions. The absence of the rbcR gene in the Dictyotales may indicate an endosymbiotic transfer from the chloroplast to the nuclear genome. Inversion of the SSC region occurred at least twice in brown algae. Once in a lineage only represented by the Ectocarpales in the present study and once in a lineage only represented by the Fucales. Photosystem genes in the chloroplasts experienced the strongest signature of purifying selection, while ribosomal protein genes in both chloroplasts and mitochondria underwent a potential weak purifying selection.

CONCLUSIONS

Variations in chloroplast genome structure among different brown algal orders are evolutionarily linked to their phylogenetic positions in the Phaeophyceae tree. Chloroplast genomes harbor more structural rearrangements than the mitochondria, despite mitochondrial genes exhibiting faster mutation rates. The position and the change in the number of boundary genes likely shaped the IR regions in the chloroplast, and the produced structural variability is important mechanistically to create gene diversity in brown algal chloroplast.

Evolutionary dynamics of plastomes in coscinodiscophycean diatoms revealed by comparative genomics

To understand the evolution of coscinodiscophycean diatoms, plastome sequences of six coscinodiscophycean diatom species were constructed and analyzed in this study, doubling the number of constructed plastome sequences in Coscinodiscophyceae (radial centrics). The platome sizes varied substantially in Coscinodiscophyceae, ranging from 119.1 kb of Actinocyclus subtilis to 135.8 kb of Stephanopyxis turris. Plastomes in Paraliales and Stephanopyxales tended to be larger than those in Rhizosoleniales and Coscinodiacales, which were due to the expansion of the inverted repeats (IRs) and to the marked increase of the large single copy (LSC). Phylogenomic analysis indicated that Paralia and Stephanopyxis clustered tightly to form the Paraliales-Stephanopyxales complex, which was sister to the Rhizosoleniales-Coscinodiscales complex. The divergence time between Paraliales and Stephanopyxales was estimated at 85 MYA in the middle Upper Cretaceous, indicating that Paraliales and Stephanopyxales appeared later than Coscinodiacales and Rhizosoleniales according to their phylogenetic relationships. Frequent losses of housekeeping protein-coding genes (PCGs) were observed in these coscinodiscophycean plastomes, indicating that diatom plastomes showed an ongoing reduction in gene content during evolution. Two acpP genes (acpP1 and acpP2) detected in diatom plastomes were found to be originated from an early gene duplication event occurred in the common progenitor after diatom emergence, rather than multiple independent gene duplications occurring in different lineages of diatoms. The IRs in Stephanopyxis turris and Rhizosolenia fallax-imbricata exhibited a similar trend of large expansion to the small single copy (SSC) and slightly small contraction from the LSC, which eventually led to the conspicuous increase in IR size. Gene order was highly conserved in Coscinodiacales, while multiple rearrangements were observed in Rhizosoleniales and between Paraliales and Stephanopyxales. Our results greatly expanded the phylogenetic breadth in Coscinodiscophyceae and gained novel insights into the evolution of plastomes in diatoms.

Effect of crude polysaccharide from seaweed, Dictyopteris divaricata (CDDP) on gut microbiota restoration and anti-diabetic activity in streptozotocin (STZ)-induced T1DM mice

Type-1 Diabetes Mellitus (T1DM) is regarded as a multifunctional, immune-related disease which causes massive destruction of islet β-cells in pancreas resulting in hyperglycemic, hypoinsulinemia and hyperlipidimic conditions. The aim of the present study, was to investigate the hypothesis that streptozotocin (STZ)-induced T1DM in Balb/c mice when treated with crude polysaccharide from seaweed, Dictyopteris divaricata (CDDP) depicts improvement in diabetes-related symptoms. Treatment with CDDP resulted in decreased body weight loss, improved food consumption and water intake disbalances. The CDDP effectively improved fasting blood glucose, oral glucose tolerance (OGTT), serum insulin, insulin secretion, rejuvenation of β-cells mass, serum lipid profile and pro-inflammatory cytokines levels. Additionally, treatment with CDDP increased the population of beneficial bacteria such as Firmicutes, Bacteroidetes and Lactobacillus at phylum, family and genus levels by 16S rRNA sequencing. Furthermore, immunohistological examination confirmed that CDDP reduces the inflammation and restored the structural morphology of colon and upraised the levels of insulin receptor substrate-1 (IRS-1), Mucin-2 (MUC-2) and tight-junction proteins (TJs) whereby maintaining the gut structures and barrier permeability. Thus, the above presented data, highlights the safe and therapeutic effects of crude polysaccharide (CDDP) from D. divaricata in the treatment and restoration of T1DM disorders and can be used as a food supplement alternative to diabetes medicine.

The organellar genomes of Silvetia siliquosa (Fucales, Phaeophyceae) and comparative analyses of the brown algae

The brown alga Silvetia siliquosa (Tseng et Chang) Serrão, Cho, Boo & Brawly is endemic to the Yellow-Bohai Sea and southwestern Korea. It is increasingly endangered due to habitat loss and excessive collection. Here, we sequenced the mitochondrial (mt) and chloroplast (cp) genomes of S. siliquosa. De novo assembly showed that the mt-genome was 36,036 bp in length, including 38 protein-coding genes (PCGs), 26 tRNAs, and 3 rRNAs, and the cp-genome was 124,991 bp in length, containing 139 PCGs, 28 tRNAs, and 6 rRNAs. Gene composition, gene number, and gene order of the mt-genome and cp-genome were very similar to those of other species in Fucales. Phylogenetic analysis revealed a close genetic relationship between S. siliquosa and F. vesiculosus, which diverged approximately 8 Mya (5.7-11.0 Mya), corresponding to the Late Miocene (5.3-11.6 Ma). The synonymous substitution rate of mitochondrial genes of phaeophycean species was 1.4 times higher than that of chloroplast genes, but the cp-genomes were more structurally variable than the mt-genomes, with numerous gene losses and rearrangements among the different orders in Phaeophyceae. This study reports the mt- and cp-genomes of the endangered S. siliquosa and improves our understanding of its phylogenetic position in Phaeophyceae and of organellar genomic evolution in brown algae.

Genomic Rearrangements and Sequence Evolution across Brown Algal Organelles

Organellar genomes serve as useful models for genome evolution and contain some of the most widely used phylogenetic markers, but they are poorly characterized in many lineages. Here, we report 20 novel mitochondrial genomes and 16 novel plastid genomes from the brown algae. We focused our efforts on the orders Chordales and Laminariales but also provide the first plastid genomes (plastomes) from Desmarestiales and Sphacelariales, the first mitochondrial genome (mitome) from Ralfsiales and a nearly complete mitome from Sphacelariales. We then compared gene content, sequence evolution rates, shifts in genome structural arrangements, and intron distributions across lineages. We confirm that gene content is largely conserved in both organellar genomes across the brown algal tree of life, with few cases of gene gain or loss. We further show that substitution rates are generally lower in plastid than mitochondrial genes, but plastomes are more variable in gene arrangement, as mitomes tend to be colinear even among distantly related lineages (with exceptions). Patterns of intron distribution across organellar genomes are complex. In particular, the mitomes of several laminarialean species possess group II introns that have T7-like ORFs, found previously only in mitochondrial genomes of Pylaiella spp. (Ectocarpales). The distribution of these mitochondrial introns is inconsistent with vertical transmission and likely reflects invasion by horizontal gene transfer between lineages. In the most extreme case, the mitome of Hedophyllum nigripes is ∼40% larger than the mitomes of close relatives because of these introns. Our results provide substantial insight into organellar evolution across the brown algae.

Comparative Analysis of Sequence Polymorphism in Complete Organelle Genomes of the ‘Golden Tide’ Seaweed Sargassum horneri between Korean and Chinese Forms

Drifting and inundating brown seaweed Sargassum horneri biomass is called “golden tide”, as it resembles golden massive algal blooms like green tides. This phenomenon occurs globally and its serious ecological impacts on coastal ecosystems have recently begun to be paid attention to. In the present study, by sequencing whole organelle genomes of Korean indigenous S. horneri, we aimed to develop novel molecular markers that can be used for differentiating indigenous from nonindigenous individuals. To this end, we analyzed sequence polymorphisms in mitochondrial (mt) and chloroplast (cp) genomes of two Korean benthic samples in comparison to Chinese ones as a reference. We mapped mt genomes of 34,620~34,628 bp and cp genomes of 123,982~124,053 bp for the Korean samples. In comparative analyses, mtDNA cytochrome c oxidase subunit II (cox2) gene showed the highest number of single nucleotide polymorphisms (SNPs) between Korean and Chinese individuals. NADH dehydrogenase subunit 7 (Nad7)-proline tRNA (trnP) intergenic spacer (IGS) in the mt genome showed a 14 bp insertion or deletion (indel) mutation. For the cp genome, we found a total of 54 SNPs, but its overall evolution rate was approximately four-fold lower than the mt genome. Interestingly, analysis of Ka/Ks ratio in the cp genome revealed a signature of positive selection on several genes, although only negative selection prevalent in mt genome. The ‘candidate’ genetic markers that we found can be applied to discriminate between Korean indigenous and nonindigenous individuals. This study will assist in developing a molecular-based early detection method for effectively managing nonindigenous S. horneri in Korean waters.

Novel module and hub genes of distinctive breast cancer associated fibroblasts identified by weighted gene co-expression network analysis

BACKGROUND

As abundant and heterogeneous stromal cells in tumor microenvironment, carcinoma-associated fibroblasts (CAFs) are critically involved in cancer progression.

METHODS

To identify co-expression module and hub genes of distinctive breast CAFs, weighted gene co-expression network analysis (WGCNA) was conducted based on the expression array results of CAFs from seven chemo-sensitive breast cancer (BC) patients and seven chemo-resistant ones before neo-adjuvant chemotherapy.

RESULTS

A total of 4916 genes were included in WGCNA, and 12 modules were determined. Module-trait assay showed that the blue module (cor = 0.97, P < 0.001) was associated with CAF-related chemo-resistance, which was enriched mainly as "inflammatory response", "interferon-gamma-mediated signaling" and "NIK/NF-kappaB signaling" pathways. Moreover, CXCL8, CXCL10, CXCL11, PLSCR1, RIPK2 and USP18 were found to be potentially associated with chemo-resistance related to CAFs and prognosis of BC.

CONCLUSIONS

Our current data offered valuable insights into the molecular mechanisms of distinctive breast CAFs, which was beneficial for revealing how chemo-resistance of BC was initiated.

Characterization of the complete chloroplast genome of an annual halophyte, Chenopodium glaucum (Amaranthaceae)

The complete chloroplast genome (plastome) of Chenopodium glaucum, an annual halophytic herb, was determined. The plastome was 152,191 bp in size, containing a large single-copy region (83,675 bp), a small single-copy region (18,130 bp), and two inverted repeats regions (25,193 bp). The overall GC content of this plastome was 37.2%. In total, 113 unique genes were annotated including 79 protein-coding genes (PCGs), 30 tRNAs and 4 rRNAs. Phylogenomic analysis showed that C. glaucum was sister to C. album.

Identification of modules and hub genes associated with platinum-based chemotherapy resistance and treatment response in ovarian cancer by weighted gene co-expression network analysis

High-grade serous ovarian carcinoma (HGSOC) is the most prevalent and malignant ovarian tumor.To identify co-expression modules and hub genes correlated with platinum-based chemotherapy resistant and sensitive HGSOC, we performed weighted gene co-expression network analysis (WGCNA) on microarray data of HGSOC with 12 resistant samples and 16 sensitive samples of GSE51373 dataset.A total of 5122 genes were included in WGCNA, and 16 modules were identified. Module-trait analysis identified that the module salmon (cor = 0.50), magenta (cor = 0.49), and black (cor = 0.45) were discovered associated with chemotherapy resistant, and the significance for these platinum-resistant modules were validated in the GSE63885 dataset. Given that the black module was validated to be the most related one, hub genes of this module, alcohol dehydrogenase 1B, cadherin 11, and vestigial like family member 3were revealed to be expressional related with platinum resistance, and could serve as prognostic markers for ovarian cancer.Our analysis might provide insight for molecular mechanisms of platinum-based chemotherapy resistance and treatment response in ovarian cancer.

The complete chloroplast genome of an annual halophyte herb, Suaeda glauca (Amaranthaceae)

The complete chloroplast genome (plastome) of Suaeda glauca, an annual halophytic herb, was determined in this study. The plastome was 149,807 bp in size, containing a large single-copy region (82,162 bp), a small single-copy region (18,191 bp), and two inverted repeats regions (24,727 bp). The overall GC content of this plastome was 36.5%. In total, 113 unique genes, including 79 protein-coding genes (PCGs), 30 tRNAs and 4 rRNAs, were annotated. Phylogenomic analysis showed that S. glauca was sister to other Suaeda species.

Characterization of the complete plastome of Atriplex centralasiatica (Chenopodiaceae), an annual halophytic herb

Atriplex centralasiatica, an annual halophytic herb, is one of the most important Chinese herbal medicines, forages and indicator plants for saline-alkali soil. In this study, we report the complete plastome of A. centralasiatica. The plastome was 152,237 bp in length and comprises a large single-copy region (83,721 bp), a small single-copy region (18,096 bp), and a pair of inverted repeats (25,210 bp). It encodes 113 unique genes, including 79 protein-coding genes (PCGs), 30 tRNAs and 4 rRNAs. The overall GC content of this plastome was 37.3%. Phylogenomic analysis based on 21 plastomes revealed that A. centralasiatica was closely related to the genus Chenopodium.

Characterization of the complete chloroplast genome of Suaeda salsa (Amaranthaceae/Chenopodiaceae), an annual succulent halophyte

Suaeda salsa, an annual succulent halophytic herb, is one of the major halophyte widely distributed in both saline inland and the intertidal zone. In this study, we report the complete chloroplast genome (plastome) of S. salsa. The plastome was 151,642 bp in length and comprises a large single-copy region (83,502 bp), a small single-copy region (17,780 bp), and a pair of inverted repeats (25,180 bp). It encodes 113 unique genes, including 79 protein-coding genes (PCGs), 30 tRNAs, and four rRNAs. The overall GC content of this plastome was 36.4%. Phylogenomic analysis based on 20 plastomes revealed that S. salsa was closely related to S. malacosperma.

Neuroprotective effect of Dictyopteris divaricata extract-capped gold nanoparticles against oxygen and glucose deprivation/reoxygenation

Combination therapy remains a promising approach to ameliorate cerebral ischemia injury. Nevertheless, the primary mechanism of the neuroprotective properties of Dictyopteris divaricata extract-capped gold nanoparticles (DD-GNPs) is not completely understood. DD-GNPs displayed maximum absorption at 525 nm and a diameter of 62.6 ± 1.2 nm, with a zeta potential value of -26.1 ± 0.6 mV. High resolution-transmission electron microscopy confirmed the spherical shape and average diameter (28.01 ± 2.03 nm). Crystalline structure and gold nanoparticle synthesis of DD-GNPs were determined by X-ray powder diffraction, and the presence of elemental gold was confirmed by energy-dispersive X-ray spectroscopy and Fourier transform-infrared spectroscopy. We examined the neuroprotective properties of DD-GNPs and explored their potential mechanisms in human SH-SY5Y neuroblastoma cells treated with oxygen and glucose deprivation/reoxygenation (OGD/R). We found that DD-GNPs inhibited OGD/R-induced release of lactate dehydrogenase (LDH), loss of cell viability, and production of reactive oxygen species. This neuroprotection was accompanied by regulation of apoptosis-related proteins, as indicated by decreased levels of cleaved-caspase-3, cleaved-PARP, cleaved-caspase-9, p53, p21, and Bax, as well as an increased level of Bcl-2. Notably, the neuroprotective effects of DD-GNPs were partially abolished by HO-1, NQO1, Nrf2, and AMPK knockdown. Our results established that DD-GNPs effectively attenuated OGD/R-stimulated neuronal injury, as evidenced by reduced neuronal injury. Even though the accumulating evidence has indicated the low toxicity and minimal side effects of GNPs, experimental clinical trials of DD-GNPs are still limited because of the lack of knowledge regarding the effects of DD-GNPs as neuroprotective agents against neurodegenerative diseases.

Phylostratigraphic analysis of gene co-expression network reveals the evolution of functional modules for ovarian cancer

Ovarian cancer (OV) is an extremely lethal disease. However, the evolutionary machineries of OV are still largely unknown. Here, we used a method that combines phylostratigraphy information with gene co-expression networks to extensively study the evolutionary compositions of OV. The present co-expression network construction yielded 18,549 nodes and 114,985 edges based on 307 OV expression samples obtained from the Genome Data Analysis Centers database. A total of 20 modules were identified as OV related clusters. The human genome sequences were divided into 19 phylostrata (PS), the majority (67.45%) of OV genes was already present in the eukaryotic ancestor. There were two strong peaks of the emergence of OV genes screened by hypergeometric test: the evolution of the multicellular metazoan organisms (PS5 and PS6, P value = 0.002) and the emergence of bony fish (PS11 and PS12, P value = 0.009). Hence, the origin of OV is far earlier than its emergence. The integrated analysis of the topology of OV modules and the phylogenetic data revealed an evolutionary pattern of OV in human, namely, OV modules have arisen step by step during the evolution of the respective lineages. New genes have evolved and become locked into a pathway, where more and more biological pathways are fixed into OV modules by recruiting new genes during human evolution.

Understanding the Evolution of Mitochondrial Genomes in Phaeophyceae Inferred from Mitogenomes of Ishige okamurae (Ishigeales) and Dictyopteris divaricata (Dictyotales)

To gain further insight into the evolution of mitochondrial genomes (mtDNAs) in Phaeophyceae, the first recorded characterization of an Ishigeophycidae mtDNA from Ishige okamurae (Yendo), and only the second recorded characterization of a Dictyotophycidae mtDNA from Dictyopteris divaricata (Okamura) Okamura are presented in this study. The 35,485 bp I. okamurae mtDNA contained 36 protein-coding genes (PCGs), 22 tRNAs, three rRNAs, and four open reading frames (orfs), and the 32,021 bp D. divaricata mtDNA harbored 35 PCGs, 25 tRNAs, three rRNAs, and three orfs. The A + T content in D. divaricata (61.69%) was the lowest recorded in sequenced brown algal mtDNAs. The I. okamurae mtDNA displayed unique genome features including an elevated start-codon usage bias for GTG, while the organization of D. divaricata mtDNA was identical to that of Dictyota dichotoma. Phylogenetic analysis based on the amino acid sequence dataset of 35 PCGs indicated that I. okamurae (Ishigeophycidae) diverged early from the Fucophycidae-Dictyotophycidae complex, which was confirmed by the comparative analysis of the mitogenome structure. The novel mitogenome data made available by this study have improved our understanding of the evolution, phylogenetics, and genomics of brown algae.

Identification of Salt Stress Responding Genes Using Transcriptome Analysis in Green Alga Chlamydomonas reinhardtii

Salinity is one of the most important abiotic stresses threatening plant growth and agricultural productivity worldwide. In green alga Chlamydomonas reinhardtii, physiological evidence indicates that saline stress increases intracellular peroxide levels and inhibits photosynthetic-electron flow. However, understanding the genetic underpinnings of salt-responding traits in plantae remains a daunting challenge. In this study, the transcriptome analysis of short-term acclimation to salt stress (200 mM NaCl for 24 h) was performed in C. reinhardtii. A total of 10,635 unigenes were identified as being differently expressed by RNA-seq, including 5920 up- and 4715 down-regulated unigenes. A series of molecular cues were screened for salt stress response, including maintaining the lipid homeostasis by regulating phosphatidic acid, acetate being used as an alternative source of energy for solving impairment of photosynthesis, and enhancement of glycolysis metabolism to decrease the carbohydrate accumulation in cells. Our results may help understand the molecular and genetic underpinnings of salt stress responses in green alga C. reinhardtii.