Characterization of the global transcriptome for Pyropia haitanensis (Bangiales, Rhodophyta) and development of cSSR markers

Recent biotechnological applications of value-added bioactive compounds from microalgae and seaweeds

Microalgae and seaweed have been consumed as food for several decades to combat starvation and food shortages worldwide. The most famous edible microalgae species are Nostoc, Spirulina, and Aphanizomenon, in addition to seaweeds, which are used in traditional medicine and food, such as Nori, which is one of the most popular foods containing Pyropia alga as a major ingredient. Recently, many applications use algae-derived polysaccharides such as agar, alginate, carrageenan, cellulose, fucoidan, mannan, laminarin, ulvan, and xylan as gelling agents in food, pharmaceuticals, and cosmetics industries. Moreover, pigments (carotenoids particularly astaxanthins, chlorophylls, and phycobilins), minerals, vitamins, polyunsaturated fatty acids, peptides, proteins, polyphenols, and diterpenes compounds are accumulated under specific cultivation and stress conditions in the algal cells to be harvested and their biomass used as a feedstock for the relevant industries and applications. No less critical is the use of algae in bioremediation, thus contributing significantly to environmental sustainability.This review will explore and discuss the various applications of microalgae and seaweeds, emphasising their role in bioremediation, recent products with algal added-value compounds that are now on the market, and novel under-developing applications such as bioplastics and nanoparticle production. Nonetheless, special attention is also drawn towards the limitations of these applications and the technologies applied, and how they may be overcome.

Identification of different proteins binding to Na, K-ATPase α1 in LPS-induced ARDS cell model by proteomic analysis

Background

Acute respiratory distress syndrome (ARDS) is characterized by refractory hypoxemia caused by accumulation of pulmonary fluid, which is related to inflammatory cell infiltration, impaired tight junction of pulmonary epithelium and impaired Na, K-ATPase function, especially Na, K-ATPase α1 subunit. Up until now, the pathogenic mechanism at the level of protein during lipopolysaccharide- (LPS-) induced ARDS remains unclear.

Methods

Using an unbiased, discovery and quantitative proteomic approach, we discovered the differentially expressed proteins binding to Na, K-ATPase α1 between LPS-A549 cells and Control-A549 cells. These Na, K-ATPase α1 interacting proteins were screened by co-immunoprecipitation (Co-IP) technology. Among them, some of the differentially expressed proteins with significant performance were identified and quantified by liquid chromatography-tandem mass spectrometry (LC–MS/MS). Data are available via ProteomeXchange with identifier PXD032209. The protein interaction network was constructed by the related Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Several differentially expressed proteins were validated by Western blot.

Results

Of identified 1598 proteins, 89 were differentially expressed proteins between LPS-A549 cells and Control-A549 cells. Intriguingly, protein–protein interaction network showed that there were 244 significantly enriched co-expression among 60 proteins in the group control-A549. while the group LPS-A549 showed 43 significant enriched interactions among 29 proteins. The related GO and KEGG analysis found evident phenomena of ubiquitination and deubiquitination, as well as the pathways related to autophagy. Among proteins with rich abundance, there were several intriguing ones, including the deubiquitinase (OTUB1), the tight junction protein zonula occludens-1 (ZO-1), the scaffold protein in CUL4B-RING ubiquitin ligase (CRL4B) complexes (CUL4B) and the autophagy-related protein sequestosome-1 (SQSTM1).

Conclusions

In conclusion, our proteomic approach revealed targets related to the occurrence and development of ARDS, being the first study to investigate significant differences in Na, K-ATPase α1 interacting proteins between LPS-induced ARDS cell model and control-A549 cell. These proteins may help the clinical diagnosis and facilitate the personalized treatment of ARDS.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12953-022-00193-3.

A draft genome of Drung cattle reveals clues to its chromosomal fusion and environmental adaptation

Drung cattle (Bos frontalis) have 58 chromosomes, differing from the Bos taurus 2n = 60 karyotype. To date, its origin and evolution history have not been proven conclusively, and the mechanisms of chromosome fusion and environmental adaptation have not been clearly elucidated. Here, we assembled a high integrity and good contiguity genome of Drung cattle with 13.7-fold contig N50 and 4.1-fold scaffold N50 improvements over the recently published Indian mithun assembly, respectively. Speciation time estimation and phylogenetic analysis showed that Drung cattle diverged from Bos taurus into an independent evolutionary clade. Sequence evidence of centromere regions provides clues to the breakpoints in BTA2 and BTA28 centromere satellites. We furthermore integrated a circulation and contraction-related biological process involving 43 evolutionary genes that participated in pathways associated with the evolution of the cardiovascular system. These findings may have important implications for understanding the molecular mechanisms of chromosome fusion, alpine valleys adaptability and cardiovascular function.

Origin, evolution, breeding, and omics of Apiaceae: a family of vegetables and medicinal plants

Many of the world's most important vegetables and medicinal crops, including carrot, celery, coriander, fennel, and cumin, belong to the Apiaceae family. In this review, we summarize the complex origins of Apiaceae and the current state of research on the family, including traditional and molecular breeding practices, bioactive compounds, medicinal applications, nanotechnology, and omics research. Numerous molecular markers, regulatory factors, and functional genes have been discovered, studied, and applied to improve vegetable and medicinal crops in Apiaceae. In addition, current trends in Apiaceae application and research are also briefly described, including mining new functional genes and metabolites using omics research, identifying new genetic variants associated with important agronomic traits by population genetics analysis and GWAS, applying genetic transformation, the CRISPR-Cas9 gene editing system, and nanotechnology. This review provides a reference for basic and applied research on Apiaceae vegetable and medicinal plants.

Influence of phosphorus on the uptake and biotransformation of arsenic in Porphyra haitanensis at environmental relevant concentrations

Edible seaweeds are rich in essential vitamins and minerals, which made them a popular food worldwide. Porphyra haitanensis is one of the most commonly consumed seaweeds with the known ability to accumulate a high level of total arsenic (As). A large number of articles have shown arsenic and phosphorus (P) interactions in microalgae due to the plant's inability to differentiate arsenate from phosphate. However, very limited information is available for edible seaweed at environmentally relevant concentrations. In this study, P. haitanensis was treated with arsenic as As^V (As1: 0.06 μM, As2: 0.4 μM, As3: 1.2 μM) and phosphorous (P1: 3.2 μM, P2: 13 μM) in a filtered seawater matrix under laboratory condition for six days. A better growth rate was found in seaweeds grown in P2 treatments. Moreover, superoxide dismutase (SOD) activity and malondialdehyde (MDA) content measurements revealed that a higher P concentration prevent seaweeds from lipid peroxidation and oxidative stress. Transcriptome studies indicated the As replacement to P has the ability to target seaweed cell membrane composition, transmembrane transport, DNA and ATP binding. The inorganic As (iAs) had a concentration of 0.54 to 4.45 mg/kg in P. haitanensis on Day 6 with As1, As2, and As3 treatments under low P regime (P1), which exceeds the limits of iAs concentration (0.1-0.5 mg/kg) in National Food Safety Standard-Limits of Pollutants in Food (GB 2762-2017). High P regime (P2) not only reduced the total As but also iAs effectively, even in the highest As treatment (As3), the iAs concentration was less than 0.5 mg/kg on Day 6. These findings provide a good insight for seafood safety guarantees and are important for the management of coastal artificial seaweed farming.

Transcriptome Co-expression Network Analysis Identifies Key Genes Regulating Conchosporangia Maturation of Pyropia haitanensis

Conchosporangia maturation is crucial for the yield of Pyropia/Porphyra. However, the molecular mechanisms underlying this process are poorly understood. In this study, we selected two strains of Pyropia haitanensis that show significant differences in conchosporangia maturation as materials to produce RNA-Seq libraries. Then, we identified key molecular pathways and genes involved in conchosporangia maturation by conducting a weighted gene co-expression network analysis. Two specific modules were identified, and included functions such as phosphorus metabolism, lipid metabolism, and the phosphatidylinositol signaling system. The hub genes that responded positively during conchosporangia maturation encoded diacylglycerol kinase (DGK) and phosphatidylinositol-3-phosphate-5-kinase, which are involved in the synthesis of phosphatidic acid, a key component of lipid metabolism. A full-length DGK sequence of P. haitanensis, designated as PhDGK1, was obtained by rapid-amplification of cDNA ends. Conserved motif and phylogenetic tree analyses showed that PhDGK1 belongs to DGK Cluster II. The transcript level of PhDGK1 increased during conchosporangia maturation in both strains, but increased earlier, and to higher levels, in the early-maturing strain than in the late-maturing strain. This pattern of gene expression was consistent with the patterns of maturity and changes in pigment contents. These results indicate that lipid metabolism plays a key role in regulating conchosporangia maturation in Pyropia spp., and that PhDGK1 might be a useful molecular marker for breeding new early-maturing strains.

Application of omics research in seaweeds with a focus on red seaweeds

Targeted 'omics' research for seaweeds, utilizing various computational and informatics frameworks, has the potential to rapidly develop our understanding of biological processes at the molecular level and contribute to solutions for the most pressing environmental and social issues of our time. Here, a systematic review into the current status of seaweed omics research was undertaken to evaluate the biological diversity of seaweed species investigated (red, green and brown phyla), the levels to which the work was undertaken (from full genome to transcripts, proteins or metabolites) and the field of research to which it has contributed. We report that from 1994 to 2021 the majority of seaweed omics research has been performed on the red seaweeds (45% of total studies), with more than half of these studies based upon two genera Pyropia and Gracilaria. A smaller number of studies examined brown seaweed (key genera Saccharina and Sargassum) and green seaweed (primarily Ulva). Overall, seaweed omics research is most highly associated with the field of evolution (46% of total studies), followed by the fields of ecology, natural products and their biosynthesis, omics methodology and seaweed-microbe interactions. Synthesis and specific outcomes derived from omics studies in the red seaweeds are provided. Together, these studies have provided a broad-scale interrogation of seaweeds, facilitating our ability to answer fundamental queries and develop applied outcomes. Crucial to the next steps will be establishing analytical tools and databases that can be more broadly utilized by practitioners and researchers across the globe because of their shared interest in the key seaweed genera.

The carbonate concentration mechanism of Pyropia yezoensis (Rhodophyta): evidence from transcriptomics and biochemical data

BACKGROUND

Pyropia yezoensis (Rhodophyta) is widely cultivated in East Asia and plays important economic, ecological and research roles. Although inorganic carbon utilization of P. yezoensis has been investigated from a physiological aspect, the carbon concentration mechanism (CCM) of P. yezoensis remains unclear. To explore the CCM of P. yezoensis, especially during its different life stages, we tracked changes in the transcriptome, photosynthetic efficiency and in key enzyme activities under different inorganic carbon concentrations.

RESULTS

Photosynthetic efficiency demonstrated that sporophytes were more sensitive to low carbon (LC) than gametophytes, with increased photosynthesis rate during both life stages under high carbon (HC) compared to normal carbon (NC) conditions. The amount of starch and number of plastoglobuli in cells corresponded with the growth reaction to different inorganic carbon (Ci) concentrations. We constructed 18 cDNA libraries from 18 samples (three biological replicates per Ci treatment at two life cycles stages) and sequenced these using the Illumina platform. De novo assembly generated 182,564 unigenes, including approximately 275 unigenes related to CCM. Most genes encoding internal carbonic anhydrase (CA) and bicarbonate transporters involved in the biophysical CCM pathway were induced under LC in comparison with NC, with transcript abundance of some PyCAs in gametophytes typically higher than that in sporophytes. We identified all key genes participating in the C4 pathway and showed that their RNA abundances changed with varying Ci conditions. High decarboxylating activity of PEPCKase and low PEPCase activity were observed in P. yezoensis. Activities of other key enzymes involved in the C4-like pathway were higher under HC than under the other two conditions. Pyruvate carboxylase (PYC) showed higher carboxylation activity than PEPC under these Ci conditions. Isocitrate lyase (ICL) showed high activity, but the activity of malate synthase (MS) was very low.

CONCLUSION

We elucidated the CCM of P. yezoensis from transcriptome and enzyme activity levels. All results indicated at least two types of CCM in P. yezoensis, one involving CA and an anion exchanger (transporter), and a second, C4-like pathway belonging to the PEPCK subtype. PYC may play the main carboxylation role in this C4-like pathway, which functions in both the sporophyte and gametophyte life cycles.

Comparison of the transcriptomes of different life history stages of the freshwater Rhodophyte Thorea hispida

Thorea hispida exclusively inhabits freshwater environments and is characterized by a triphasic life history. In this study, the organelle genomes and transcriptomes of different life history stages of T. hispida were examined using next generation sequencing. The chloroplast and mitochondrial genomes of the chantransia stage were 175,747 and 25,411 bp in length, respectively. The chantransia stage was highly similar to the gametophyte stage based on comparisons of organelle genomes and phylogenetic reconstruction. Transcriptomic comparisons of two stages found that ribosome-related genes were the most up-regulated in the gametophyte stage of T. hispida. Seven meiosis-specific genes, including SPO11 initiator of meiotic double-stranded breaks(spo11), meiotic nuclear divisions 1(mnd1), RAD51 recombinase(rad51), mutS homolog 4(msh4), mutS homolog 5(msh5), REC8 meiotic recombination protein(rec8), and DNA helicase Mer3(mer3), were differentially regulated between the two life history stages. The organelle genomes and transcriptomes from T. hispida provided in this study will be valuable for future studies of freshwater red algae.

Health Functionality and Quality Control of Laver (Porphyra, Pyropia): Current Issues and Future Perspectives as an Edible Seaweed

The growing interest in laver as a food product and as a source of substances beneficial to health has led to global consumer demand for laver produced in a limited area of northeastern Asia. Here we review research into the benefits of laver consumption and discuss future perspectives on the improvement of laver product quality. Variation in nutritional/functional values among product types (raw and processed (dried, roasted, or seasoned) laver) makes product-specific nutritional analysis a prerequisite for accurate prediction of health benefits. The effects of drying, roasting, and seasoning on the contents of both beneficial and harmful substances highlight the importance of managing laver processing conditions. Most research into health benefits has focused on substances present at high concentrations in laver (porphyran, Vitamin B12, taurine), with assessment of the expected effects of laver consumption. Mitigation of chemical/microbiological risks and the adoption of novel technologies to exploit under-reported biochemical characteristics of lavers are suggested as key strategies for the further improvement of laver product quality. Comprehensive analysis of the literature regarding laver as a food product and as a source of biomedical compounds highlights the possibilities and challenges for application of laver products.

Characterization of Nme5-Like Gene/Protein from the Red Alga Chondrus Crispus

The Nme gene/protein family of nucleoside diphosphate kinases (NDPK) was originally named after its member Nm23-H1/Nme1, the first identified metastasis suppressor. Human Nme proteins are divided in two groups. They all possess nucleoside diphosphate kinase domain (NDK). Group I (Nme1-Nme4) display a single type NDK domain, whereas Group II (Nme5-Nme9) display a single or several different NDK domains, associated or not associated with extra-domains. Data strongly suggest that, unlike Group I, none of the members of Group II display measurable NDPK activity, although some of them autophosphorylate. The multimeric form is required for the NDPK activity. Group I proteins are known to multimerize, while there are no data on the multimerization of Group II proteins. The Group II ancestral type protein was shown to be conserved in several species from three eukaryotic supergroups. Here, we analysed the Nme protein from an early branching eukaryotic lineage, the red alga Chondrus crispus. We show that the ancestral type protein, unlike its human homologue, was fully functional multimeric NDPK with high affinity to various types of DNA and dispersed localization throughout the eukaryotic cell. Its overexpression inhibits both cell proliferation and the anchorage-independent growth of cells in soft agar but fails to deregulate cell apoptosis. We conclude that the ancestral gene has changed during eukaryotic evolution, possibly in correlation with the protein function.

Insight into transketolase of Pyropia haitanensis under desiccation stress based on integrative analysis of omics and transformation

BACKGROUND

Pyropia haitanensis, distributes in the intertidal zone, can tolerate water losses exceeding 90%. However, the mechanisms enabling P. haitanensis to survive harsh conditions remain uncharacterized. To elucidate the mechanism underlying P. haitanensis desiccation tolerance, we completed an integrated analysis of its transcriptome and proteome as well as transgenic Chlamydomonas reinhardtii carrying a P. haitanensis gene.

RESULTS

P. haitanensis rapidly adjusted its physiological activities to compensate for water losses up to 60%, after which, photosynthesis, antioxidant systems, chaperones, and cytoskeleton were activated to response to severe desiccation stress. The integrative analysis suggested that transketolase (TKL) was affected by all desiccation treatments. Transgenic C. reinhardtii cells overexpressed PhTKL grew better than the wild-type cells in response to osmotic stress.

CONCLUSION

P. haitanensis quickly establishes acclimatory homeostasis regarding its transcriptome and proteome to ensure its thalli can recover after being rehydrated. Additionally, PhTKL is vital for P. haitanensis desiccation tolerance. The present data may provide new insights for the breeding of algae and plants exhibiting enhanced desiccation tolerance.

Regulatory mechanisms underlying the maintenance of homeostasis in Pyropia haitanensis under hypersaline stress conditions

Intertidal macroalgae are highly resistant to hypersaline stress conditions. However, the underlying mechanism remains unknown. In the present study, the mechanism behind Pyropia haitanensis responses to two hypersaline stress conditions [100‰ (HSS_100) and 110‰ (HSS_110)] was investigated via analyses of physiological and transcriptomic changes. We observed that the differences between the responses of Py. haitanensis to HSS_100 and HSS_110 conditions involved the following three aspects: osmotic regulation, ionic homeostasis, and adjustment to secondary stresses. First, the water retention of Py. haitanensis was maintained through increased expansin production under HSS_100 conditions, while cell wall pectin needed to be protected from hydrolysis via the increased abundance of a pectin methylesterase inhibitor under HSS_110 conditions. Meanwhile, Py. haitanensis achieved stable and rapid osmotic adjustments because of the coordinated accumulation of inorganic ions (K⁺, Na⁺, and Cl^-) and organic osmolytes (glycine betaine and trehalose) under HSS_100 conditions, but not under HSS_110 conditions. Second, Py. haitanensis maintained a higher K⁺/Na⁺ ratio under HSS_100 conditions than under HSS_110 conditions, mainly via the export of Na⁺ into the apoplast rather than compartmentalizing it into the vacuoles, and the enhanced uptake and retention of K⁺. However, K⁺/Na⁺ homeostasis was not completely disrupted during a short-term exposure to HSS_110 conditions. Finally, the Py. haitanensis antioxidant system scavenged more ROS and synthesized more heat shock proteins under HSS_100 conditions than under HSS_110 conditions, although thalli may have been able to maintain a certain redox balance during a short-term exposure to HSS_110 conditions. These differences may explain why Py. haitanensis can adapt to HSS_100 conditions rather than HSS_110 conditions, and also why the thalli exposed to HSS_110 conditions can recover after being transferred to normal seawater. Thus, the data presented herein may elucidate the mechanisms enabling Pyropia species to tolerate the sudden and periodic changes in salinity typical of intertidal systems.

Gene Expression Profiling in Ischemic Postconditioning to Alleviate Mouse Liver Ischemia/Reperfusion Injury

Ischemic postconditioning (IPO) attenuates hepatic ischemia/reperfusion (I/R) injury. However, little is known about the underlying biological pathophysiology, which could be, at least in part, informed by exploring the transcriptomic changes using next-generation RNA sequencing (RNA-Seq). In this study, 18 mice (C57BL/6) were involved and randomly assigned to three groups: normal (n=6), I/R (n=6, subjected to 70% hepatic I/R), and IR+IPO (n=6, applying IPO to mice with I/R injury). We randomly selected 3 mice per group and extracted their liver tissues for next-generation RNA-Seq. We performed a bioinformatics analysis for two comparisons: normal vs. I/R and I/R vs. IR+IPO. From the analysis, 2416 differentially expressed genes (DEGs) were identified (p < 0.05 and fold change ≥ 1.5). Gene ontology (GO) analysis revealed that these genes were mainly related to cellular metabolic processes, nucleic acids and protein binding processes. The enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for the DEGs were the mitogen-activated protein kinase (MAPK), IL-17 signalling pathway, regulating pluripotency of stem cells, and insulin resistance pathway. Validation of 12 selected DEGs by qRT-PCR showed that Cyr61, Atf3, Nr4a1, Gdf15, Osgin1, Egr1, Epha2, Dusp1, Dusp6, Gadd45a and Gadd45b were significantly amplified. Finally, a protein-protein interaction (PPI) network constructed to determine interactions of these 11 DEGs. In summary, by exploring gene expression profiling in regard to hepatic I/R and IPO using next-generation RNA-Seq, we suggested a few progression-related genes and pathways, providing some clues for future experimental research.

Transcriptome analysis of the typical freshwater rhodophytes Sheathia arcuata grown under different light intensities

The Rhodophyta Sheathia arcuata is exclusively distributed in freshwater, constituting an important component in freshwater flora. This study presents the first transcriptome profiling of freshwater Rhodophyta taxa. A total of 161,483 assembled transcripts were identified, annotated and classified into different biological categories and pathways based on BLAST against diverse databases. Different gene expression patterns were caused principally by different irradiances considering the similar water conditions of the sampling site when the specimens were collected. Comparison results of gene expression levels under different irradiances revealed that photosynthesis-related pathways significantly up-regulated under the weak light. Molecular responses for improved photosynthetic activity include the transcripts corresponding to antenna proteins (LHCA1 and LHCA4), photosynthetic apparatus proteins (PSBU, PETB, PETC, PETH and beta and gamma subunits of ATPase) and metabolic enzymes in the carbon fixation. Along with photosynthesis, other metabolic activities were also regulated to optimize the growing and development of S. arcuata under appropriate sunlight. Protein-protein interactive networks revealed the most responsive up-expressed transcripts were ribosomal proteins. The de-novo transcriptome assembly of S. arcuata provides a foundation for further investigation on the molecular mechanism of photosynthesis and environmental adaption for freshwater Rhodophyta.

Transcriptomic study to understand thermal adaptation in a high temperature-tolerant strain of Pyropia haitanensis

Pyropia haitanensis, a high-yield commercial seaweed in China, is currently undergoing increasing levels of high-temperature stress due to gradual global warming. The mechanisms of plant responses to high temperature stress vary with not only plant type but also the degree and duration of high temperature. To understand the mechanism underlying thermal tolerance in P. haitanensis, gene expression and regulation in response to short- and long-term temperature stresses (SHS and LHS) was investigated by performing genome-wide high-throughput transcriptomic sequencing for a high temperature tolerant strain (HTT). A total of 14,164 differential expression genes were identified to be high temperature-responsive in at least one time point by high-temperature treatment, representing 41.10% of the total number of unigenes. The present data indicated a decrease in the photosynthetic and energy metabolic rates in HTT to reduce unnecessary energy consumption, which in turn facilitated in the rapid establishment of acclimatory homeostasis in its transcriptome during SHS. On the other hand, an increase in energy consumption and antioxidant substance activity was observed with LHS, which apparently facilitates in the development of resistance against severe oxidative stress. Meanwhile, ubiquitin-mediated proteolysis, brassinosteroids, and heat shock proteins also play a vital role in HTT. The effects of SHS and LHS on the mechanism of HTT to resist heat stress were relatively different. The findings may facilitate further studies on gene discovery and the molecular mechanisms underlying high-temperature tolerance in P. haitanensis, as well as allow improvement of breeding schemes for high temperature-tolerant macroalgae that can resist global warming.

Transcriptome-wide identification of optimal reference genes for expression analysis of Pyropia yezoensis responses to abiotic stress

Background

Pyropia yezoensis, a marine red alga, is an ideal research model for studying the mechanisms of abiotic stress tolerance in intertidal seaweed. Real-time quantitative polymerase chain reaction (RT-qPCR) is the most commonly used method to analyze gene expression levels. To accurately quantify gene expression, selection and validation of stable reference genes is required.

Results

We used transcriptome profiling data from different abiotic stress treatments to identify six genes with relatively stable expression levels: MAP, ATPase, CGS1, PPK, DPE2, and FHP. These six genes and three conventional reference genes, UBC, EF1-α, and eif4A, were chosen as candidates for optimal reference gene selection. Five common statistical approaches (geNorm, ΔCt method, NormFinder, BestKeeper, and ReFinder) were used to identify the stability of each reference gene. Our results show that: MAP, UBC, and FHP are stably expressed in all analyzed conditions; CGS1 and UBC are stably expressed under conditions of dehydration stress; and MAP, UBC, and CGS1 are stably expressed under conditions of temperature stress.

Conclusion

We have identified appropriate reference genes for RT-qPCR in P. yezoensis under different abiotic stress conditions which will facilitate studies of gene expression under these conditions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4643-8) contains supplementary material, which is available to authorized users.

Candidate Genes as Biomarkers in Lipopolysaccharide-Induced Acute Respiratory Distress Syndrome Based on mRNA Expression Profile by Next-Generation RNA-Seq Analysis

Up until now, the regulation mechanism at the level of gene during lipopolysaccharide- (LPS-) induced acute respiratory distress syndrome (ARDS) remains unclear. The discovery of differentially expressed genes (DEGs) between LPS-induced ARDS rats and normal rats by next-generation RNA sequencing analysis is of particular interest for the current study. These DEGs may help clinical diagnosis of ARDS and facilitate the selection of the optimal treatment strategy. Randomly, 20 rats were equally divided into 2 groups, the control group and the LPS group. Three rats from each group were selected at random for RNA sequencing analysis. Sequence reads were obtained from Illumina HiSeq4000 and mapped onto the rat reference genome RN6 using Hisat2. We identified 5244 DEGs (Fold_Change > 1.5, and P < 0.05) in the lung tissues from LPS-treated rats compared with normal rats, including 1413 upregulated and 3831 downregulated expressed genes. Lots of chemokine family members were among the most upregulated genes in LPS group. Gene ontology (GO) analysis revealed that almost all of the most enriched and meaningful biological process terms were mainly involved in the functions like immune-inflammation response and the pathways like cytokine-cytokine receptor interaction. We also found that, as for GO molecular function terms, the enriched terms were mainly related to chemokines and cytokines. DEGs with fold change over 100 were verified by quantitative real-time polymerase chain reaction and reanalyzed by gene-gene coexpression network, and the results elucidated central roles of chemokines in LPS-induced ARDS. Our results revealed some new biomarkers for uncovering mechanisms and processes of ARDS.

An Efficient PCR-RFLP Method for the Rapid Identification of Korean Pyropia Species

The present study utilizes polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis using partial plastid rbcL and mitochondrial trnC–trnP gene sequences to distinguish the six representative Pyropia species produced via mariculture in Korea. The rbcL, trnC, and trnP sequences of 15 Pyropia species from the NCBI database were aligned to determine specific restriction enzyme sites of the six Pyropia species. To confirm the presence of restriction sites of eight enzymes, PCR amplicons were digested as follows: a 556 bp fragment within the rbcL region of chloroplast DNA was confirmed in P. yezoensis using BglI, whereas Tth111I, AvaII, BsrI, and BsaAI enzymes produced fragments of 664, 271, 600, and 510 bp, respectively, from the rps11–trnG region of mitochondrial DNA in P. seriata, P. dentata, P. suborbiculata, and P. haitanensis. In the case of P. pseudolinearis, HindIII, SacII, and SphI enzymes each had two cleavage sites, at positions 174 and 825, 788 and 211, and 397 and 602 bp, respectively. All six species were successfully distinguished using these eight restriction enzymes. Therefore, we propose that PCR-RFLP analysis is an efficient tool for the potential use of distinguishing between the six Pyropia species cultivated via mariculture in Korea.

Investigating right ovary degeneration in chick embryos by transcriptome sequencing

In asymmetric chick gonads, the left and right female gonads undergo distinct programs during development, generating a functional ovary on the left side only. Despite some progress being made in recent years, the mechanisms of molecular regulation remain incompletely understood, and little genomic information is available regarding the degeneration of the right ovary in the chick embryo testis. In this study, we performed transcriptome sequencing to investigate differentially expressed genes in the left and right ovaries and gene functions at two critical time points; embryonic days 6 (E6) and 10 (E10). Using high-throughput RNA-sequencing technologies, 539 and 1046 genes were identified as being significantly differentially expressed between 6R-VS-6L and 10R-VS-10L. Gene ontology analysis of the differentially expressed genes revealed enrichment in functional pathways. Among these, candidate genes associated with degeneration of the right ovary in the chick embryo were identified. Identification of a pathway involved in ovarian degeneration provides an important resource for the further study of its molecular mechanisms and functions.

Distribution, function and evolution characterization of microsatellite in Sargassum thunbergii (Fucales, Phaeophyta) transcriptome and their application in marker development

Using transcriptome data to mine microsatellite and develop markers has growingly become prevalent. However, characterizing the possible function of microsatellite is relatively rare. In this study, we explored microsatellites in the transcriptome of the brown alga Sargassum thunbergii and characterized the frequencies, distribution, function and evolution, and developed primers to validate these microsatellites. Our results showed that Tri-nucleotide is the most abundant, followed by di- and mono-nucleotide. The length of microsatellite was significantly affected by the repeat motif size. The density of microsatellite in the CDS region is significantly lower than that in the UTR region. The annotation of the transcripts containing microsatellite showed that 573 transcripts have GO terms and can be categorized into 42 groups. Pathways enrichment showed that microsatellites were significantly overrepresented in the genes involved in pathways such as Ubiquitin mediated proteolysis, RNA degradation, Spliceosome, etc. Primers flanking 961 microsatellite loci were designed, and among the 30 pairs of primer selected randomly for availability test, 23 were proved to be efficient. These findings provided new insight into the function and evolution of microsatellite in transcriptome, and the identified microsatellite loci within the annotated gene will be useful for developing functional markers in S. thunbergii.

Genome-wide expression profiles of Pyropia haitanensis in response to osmotic stress by using deep sequencing technology

BACKGROUND

Pyropia haitanensis is an economically important marine crop grown in harsh intertidal habitats of southern China; it is also an excellent model system for studying mechanisms of stress tolerance. To understand the molecular mechanisms underlying osmotic tolerance and adaptation to intertidal environments, a comprehensive analysis of genome-wide gene expression profiles in response to dehydration and rehydration in Py. haitanensis was undertaken using digital gene expression profile (DGE) approaches combined with de novo transcriptome sequencing.

RESULTS

RNA-sequencing of the pooled RNA samples from different developmental phases and stress treatments was performed, which generated a total of 47.7 million clean reads. These reads were de novo assembled into 28,536 unigenes (≥ 200 bp), of which 18,217 unigenes (63.83 %) were annotated in at least one reference database. DGE analysis was performed on four treatments (two biological replicates per treatment), which included moderate dehydration, severe dehydration, rehydration, and normal conditions. The number of raw reads per sample ranged from 12.47 to 15.79 million, with an average of 14.69 million reads per sample. After quality filtering, the number of clean reads per sample ranged from 11.83 to 15.04 million. All distinct sequencing reads were annotated using the transcriptome of Py. haitanensis as reference. A total of 1,681 unigenes showed significant differential expression between moderate dehydration and normal conditions, in which 977 genes were upregulated, and 704 genes were downregulated. Between severe dehydration and normal conditions, 1,993 unigenes showed significantly altered expression, which included both upregulated (1,219) and downregulated genes (774). In addition, 1,086 differentially expressed genes were detected between rehydration and normal conditions, of which 720 genes were upregulated and 366 unigenes were downregulated. Most gene expression patterns in response to dehydration differed from that of rehydration, except for the synthesis of unsaturated fatty acids, several transcription factor families, and molecular chaperones that have been collectively implicated in the processes of dehydration and rehydration in Py. haitanensis.

CONCLUSIONS

Taken together, these data provide a global high-resolution analysis of gene expression changes during osmotic stress that could potentially serve as a key resource for understanding the biology of osmotic acclimation in intertidal red seaweed.

Recent Major Advances of Biotechnology and Sustainable Aquaculture in China

Background:

Global aquaculture production has increased continuously over the last five decades, and particularly in China. Its aquaculture has become the fastest growing and most efficient agri-sector, with production accounting for more than 70% of the world’s aquaculture output. In the new century, with serious challenges regarding population, resources and the environment, China has been working to develop high-quality, effective, healthy, and sustainable blue agriculture through the application of modern biotechnology. Sound knowledge related to the biology and ecology of aquatic organisms has laid a solid foundation and provided the innovation and technology for rapid development of the aquaculture industry. Marine biotechnology, which is enabling solutions for ocean productivity and sustainability, has been promoted since the last decades of the 20th Century in China.

Objective:

In this article, priority areas of research, mainly genetic breeding, omics studies, novel production systems, biosecurity, bioprocesses and biorefinery, as well as the major progress of marine biotechnology R&D in China are reviewed.

Conclusion:

Current innovative achievements in China are not enough and the level and frequency of academic advancements must be improved. International cooperation and assistance remain crucial for the success of marine biotechnology.

Digital gene expression analysis of transcriptomes in lipopolysaccharide-induced acute respiratory distress syndrome

BACKGROUND

The mortality from acute respiratory distress syndrome (ARDS) is high, and its exact pathogenesis remains unclear, which forms a major obstacle for prevention and treatment of this disease. In the present study, we used digital gene expression (DGE) to detect the differentially expressed genes of the lung at 4h after lipopolysaccharide (LPS) exposure in a mouse model.

METHODS

Mice were treated with LPS or control saline by intratracheal instillation for 4h, and their lung tissues were collected for DGE analysis. We used a false discovery rate ≤0.001 and an absolute value of the log2 ratio≥1 as the thresholds for judging the significance of any difference in gene expression between the two members of each pair of mice.

RESULTS

We obtained 3,387,842 clean tags (i.e., after filtering to remove potentially erroneous tags) and about 84,513 corresponding distinct clean tags (i.e., types of tag). Approximately 91.20% of the clean tags could be mapped, and 82.71% could be uniquely mapped, to the reference tags, and 3.82% were unknown tags. At least 2200 differentially expressed genes were identified and analyzed for enrichment of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway. Twenty genes with the greatest difference in expression levels between the two members of every pair of mice were chosen. The majority of these genes are involved in signaling transduction, molecular adhesion, and metabolic pathways.

CONCLUSIONS

Using the powerful technology of DGE, we present, to our knowledge, the first in-depth transcriptomic analysis of mouse lungs after LPS exposure. We found some differentially expressed genes that might play important roles in the pathogenesis of ARDS.

Glycerol-3-phosphate metabolism plays a role in stress response in the red alga Pyropia haitanensis

Glycerol-3-phosphate (G3P) has been suggested as a novel regulator of plant defense signaling, however, its role in algal resistance remains largely unknown. The glycerol kinase (also designated as NHO1) and NAD-dependent G3P dehydrogenase (GPDH) are two key enzymes involved in the G3P biosynthesis. In our study, we cloned the full-length cDNA of NHO1 (NHO1Ph ) and GPDH (GPDHP h ) from the red alga Pyropia haitanensis (denoted as NHO1Ph and GPDHP h ) and examined their expression level under flagellin peptide 22 (flg22) stimulation or heat stress. We also measured the level of G3P and floridoside (a downstream product of G3P in P. haitanensis) under flg22 stimulation or heat stress. Both NHO1Ph and GPDHP h shared high sequence identity and structural conservation with their orthologs from different species, especially from red algae. Phylogenetic analysis showed that NHO1s and GPDHs from red algae were closely related to those from animals. Under flg22 stimulation or heat stress, the expression levels of NHO1Ph and GPDHP h were up-regulated, G3P levels increased, and the contents of floridoside decreased. But the floridoside level increased in the recovery period after heat stress. Taken together, we found that G3P metabolism was associated with the flg22-induced defense response and heat stress response in P. haitanensis, indicating the general conservation of defense response in angiosperms and algae. Furthermore, floridoside might also participate in the stress resistance of P. haitanensis.

Characterization of the global transcriptome for cotton (Gossypium hirsutum L.) anther and development of SSR marker

Cotton is an important fiber plant, and it's attractive to elucidate the molecular mechanism of anther development due to the close relationship between the anther fertility and boll-setting, and also fiber yield. In the present paper, 47.2 million paired-end reads with average length of 82.87 bp from the anthers of TM-1 (Gossypium hirsutum L.), a genetic standard line, were generated through transcriptome sequencing, and 210,965 unigenes of more than 100 bp were obtained. BLAST, KEGG, COG, and GO analyses showed that the genes were enriched in the processes of transcription, translation, and post-translation as well as hormone signal transduction, the transcription factor families, and cell wall-related genes mainly participating in cell expansion and carbohydrate metabolism. Further analysis identified 11,153 potential SSRs. A suit of 5122 primer pair sequences were designed, and 82 of 300 randomly selected primer pairs produced reproducible amplicons that were polymorphic among 22 cotton accessions from G. hirsutum, Gossypium barbadense and Gossypium arboreum. The UPGMA clustering analysis further confirmed high quality and effectiveness of these novel SSR markers. The present study provided insights into the transcriptome profile of the cotton and established a public information platform for functional genomics and molecular breeding.

Different responses to heat shock stress revealed heteromorphic adaptation strategy of Pyropia haitanensis (Bangiales, Rhodophyta)

Pyropia has a unique heteromorphic life cycle with alternation stages between thallus and conchocelis, which lives at different water temperatures in different seasons. To better understand the different adaptation strategies for temperature stress, we tried to observe comparative biochemical changes of Pyropia haitanensis based on a short term heat shock model. The results showed that: (1) At normal temperature, free-living conchocelis contains significantly higher levels of H₂O₂, fatty acid-derived volatiles, the copy number of Phrboh and Phhsp70 genes,the activities of NADPH oxidase and floridoside than those in thallus. The released H₂O₂ and NADPH oxidase activity of conchocelis were more than 7 times higher than those of thallus. The copy number of Phrboh in conchocelis was 32 times that in thallus. (2) After experiencing heat shock at 35°C for 30 min, the H₂O₂ contents, the mRNA levels of Phrboh and Phhsp70, NADPH oxidase activity and the floridoside content in thallus were all significantly increased. The mRNA levels of Phrboh increased 5.78 times in 5 min, NADPH oxidase activity increased 8.45 times in 20 min. (3) Whereas, in conchocelis, the changes in fatty acids and their down-stream volatiles predominated, significantly increasing levels of saturated fatty acids and decreasing levels of polyunsaturated fatty acids occurred, and the 8-carbon volatiles were accumulated. However, the changes in H₂O₂ content and expression of oxidant-related genes and enzymatic activity were not obvious. Overall, these results indicate that conchocelis maintains a high level of active protective apparatus to endure its survival at high temperature, while thallus exhibit typical stress responses to heat shock. It is concluded that Pyropia haitanensis has evolved a delicate strategy for temperature adaptation for its heteromorphic life cycle.

Analysis of horizontal genetic transfer in red algae in the post-genomics age

The recently published genome of the unicellular red alga Porphyridium purpureum revealed a gene-rich, intron-poor species, which is surprising for a free-living mesophile. Of the 8,355 predicted protein-coding regions, up to 773 (9.3%) were implicated in horizontal genetic transfer (HGT) events involving other prokaryote and eukaryote lineages. A much smaller number, up to 174 (2.1%) showed unambiguous evidence of vertical inheritance. Together with other red algal genomes, nearly all published in 2013, these data provide an excellent platform for studying diverse aspects of algal biology and evolution. This novel information will help investigators test existing hypotheses about the impact of endosymbiosis and HGT on algal evolution and enable comparative analysis within a more-refined, hypothesis-driven framework that extends beyond HGT. Here we explore the impacts of this infusion of red algal genome data on addressing questions regarding the complex nature of algal evolution and highlight the need for scalable phylogenomic approaches to handle the forthcoming deluge of sequence information.

De novo transcriptome analysis of Liriodendron chinense petals and leaves by Illumina sequencing

Liriodendron chinense (Hemsl.) Sarg is an endangered species and occupies a pivotal position in phylogenetic studies of flowering plants, while its genomic resources are limited. In this study, we performed transcriptome sequencing for L. chinense petals and leaves using the Illumina paired-end sequencing technique. Approximately 17.02-Gb clean reads were obtained, and de novo assembly generated 87,841 unigenes, with an average length of 778 bp. Of these, there were 65,535 (74.61%) unigenes with significant similarity to publically available plant protein sequences. There were 3386 genes identified as significant differentially expressed between petals and leaves, among them 2969 (87.68%) were up-regulated and 417 (12.31%) down-regulated in petals. Metabolic pathway analysis revealed that 25 unigenes were predicted to be responsible for the biosynthesis of carotenoids, with 7 genes differentially expressed between these two tissues. This report is the first to identify genes associated with carotenoid biosynthesis in Liriodendron and represents a valuable resource for future genomic studies on the endangered species L. chinense.

De novo transcriptome assembly for the tropical grass Panicum maximum Jacq

Guinea grass (Panicum maximum Jacq.) is a tropical African grass often used to feed beef cattle, which is an important economic activity in Brazil. Brazil is the leader in global meat exportation because of its exclusively pasture-raised bovine herds. Guinea grass also has potential uses in bioenergy production due to its elevated biomass generation through the C4 photosynthesis pathway. We generated approximately 13 Gb of data from Illumina sequencing of P. maximum leaves. Four different genotypes were sequenced, and the combined reads were assembled de novo into 38,192 unigenes and annotated; approximately 63% of the unigenes had homology to other proteins in the NCBI non-redundant protein database. Functional classification through COG (Clusters of Orthologous Groups), GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses showed that the unigenes from Guinea grass leaves are involved in a wide range of biological processes and metabolic pathways, including C4 photosynthesis and lignocellulose generation, which are important for cattle grazing and bioenergy production. The most abundant transcripts were involved in carbon fixation, photosynthesis, RNA translation and heavy metal cellular homeostasis. Finally, we identified a number of potential molecular markers, including 5,035 microsatellites (SSRs) and 346,456 single nucleotide polymorphisms (SNPs). To the best of our knowledge, this is the first study to characterize the complete leaf transcriptome of P. maximum using high-throughput sequencing. The biological information provided here will aid in gene expression studies and marker-assisted selection-based breeding research in tropical grasses.