Aerobic photosynthetic organisms such as algae produce reactive oxygen species (ROS) as by-products of metabolism. ROS damage biomolecules such as proteins and lipids in cells, but also act as signaling molecules. The mechanisms that maintain the metabolic balance in aerobic photosynthetic organisms and how the cells specifically respond to different levels of ROS are unclear. Glutathione peroxidase (GPX) enzymes detoxify hydrogen peroxide or organic hydroperoxides, and thus are important components of the antioxidant system. In this study, we employed a Chlamydomonas reinhardtii glutathione peroxidase knockout (gpx5) mutant to identify the genetic response to singlet oxygen (¹O₂) generated by the photosensitizer rose bengal (RB). To this end, we compared the transcriptomes of the parental strain CC4348 and the gpx5 mutant sampled before, and 1 h after, the addition of RB. Functional annotation of differentially expressed genes showed that genes encoding proteins related to ROS detoxification, stress-response-related molecular chaperones, and ubiquitin–proteasome pathway genes were upregulated in CC4338. When GPX5 was mutated, higher oxidative stress specifically induced the TCA cycle and enhanced mitochondrial electron transport. Transcription of selenoproteins and flagellar-associated proteins was depressed in CC4348 and the gpx5 mutant. In addition, we found iron homeostasis played an important role in maintaining redox homeostasis, and we uncovered the relationship between ¹O₂ stress and iron assimilation, as well as selenoproteins. Based on the observed expression profiles in response to different levels of oxidative stress, we propose a model for dose-dependent responses to different ROS levels in Chlamydomonas.

High-throughput RNA sequencing unveiled the complexity of transcriptome and significantly increased the records of long noncoding RNAs (lncRNAs), which were reported to participate in a variety of biological processes. Identification of lncRNAs is a key step in lncRNA analysis, and a bunch of bioinformatics tools have been developed for this purpose in recent years. While these tools allow us to identify lncRNA more efficiently and accurately, they may produce inconsistent results, making selection a confusing issue. We compared the performance of 41 analysis models based on 14 software packages and different data sets, including high-quality data and low-quality data from 33 species. In addition, computational efficiency, robustness, and joint prediction of the models were explored. As a practical guidance, key points for lncRNA identification under different situations were summarized. In this investigation, no one of these models could be superior to others under all test conditions. The performance of a model relied to a great extent on the source of transcripts and the quality of assemblies. As general references, FEELnc_all_cl, CPC, and CPAT_mouse work well in most species while COME, CNCI, and lncScore are good choices for model organisms. Since these tools are sensitive to different factors such as the species involved and the quality of assembly, researchers must carefully select the appropriate tool based on the actual data. Alternatively, our test suggests that joint prediction could behave better than any single model if proper models were chosen. All scripts/data used in this research can be accessed at http://bioinfo.ihb.ac.cn/elit.

Background: Grass carp (Ctenopharyngodon idellus) are important species in Asian aquaculture. A draft genome for grass carp has already been published in 2015. However, there is still a requirement for a suitable genetic linkage map to arrange scaffolds on chromosomal frameworks. QTL analysis is a powerful tool to detect key locations for quantitative traits, especially in aquaculture. There no growth related QTLs of grass carp have been published yet. Even the growth trait is one of the focuses in grass carp culture.Results: In this study, a pair of distantly related parent grass carps and their 100 six-month-old full-sib offspring were used to construct a high-density genetic map with 6429 single nucleotide polymorphisms (SNPs) by 2b-RAD technology. The total length of the consensus map is 5553.43 cM with the average marker interval of 1.92 cM. The map has a good collinearity with both the grass carp draft genome and the zebrafish genome, and it assembled 89.91% of the draft genome to a chromosomal level. Additionally, according to the growth-related traits of progenies, 30 quantitative trait loci (QTLs), including 7 for body weight, 9 for body length, 5 for body height and 9 for total length, were identified in 16 locations on 5 linkage groups. The phenotypic variance explained for these QTLs varies from 13.4 to 21.6%. Finally, 17 genes located in these regions were considered to be growth-related because they either had functional mutations predicted from the resequencing data of the parents.Conclusion: A high density genetic linkage map of grass carp was built and it assembled the draft genome to a chromosomal level. Thirty growth related QTLs were detected. After the cross analysis of Parents resequencing data, 17 candidate genes were obtained for further researches.

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Intermuscular bones (IBs) are small spicule-like bones in the muscular septum of fish, which affect their edible and economic value. The molecular mechanism of IB development is still uncertain. Numerous studies have shown that the ceRNA network, which is composed of mRNA, lncRNA, and miRNA, plays an important regulatory role in bone development. In this study, we compared the mRNA, lncRNA, and miRNA expression profiles in different IB development segments of zebrafish. The development of IBs includes two main processes, which are formation and growth. A series of genes implicated in the formation and growth of IBs were identified through gene differential expression analysis and expression pattern analysis. Functional enrichment analysis showed that the functions of genes implicated in the regulation of the formation and growth of IBs were quite different. Ribosome and oxidative phosphorylation signaling pathways were significantly enriched during the formation of IBs, suggesting that many proteins are required to form IBs. Several pathways known to be associated with bone development have been shown to play an important role in the growth of IBs, including calcium, ECM-receptor interaction, Wnt, TGF-β, and hedgehog signaling pathways. According to the targeting relationship and expression correlation of mRNA, lncRNA, and miRNA, the ceRNA networks associated with the growth of IBs were constructed, which comprised 33 mRNAs, 9 lncRNAs, and 7 miRNAs. This study provides new insight into the molecular mechanism of the development of IBs.

Grass carp is an important commercial fish widely cultured in China. Large range of temperature, in particular extremely low temperature, has dramatic effects on the aquaculture of this teleost. However, there is relatively little research on the molecular responses in the fish exposed to cold. Given the limited vision of approaches targeting individual genes, we investigated the transcriptome profiles of brain in response to cold in order to comprehensively characterize molecular mechanisms behind it. This study indicated that the estrogen signaling pathway was inhibited in brain when grass carp acclimated to low temperature, while terpenoid backbone biosynthesis pathway and steroid biosynthesis pathway were significantly activated. Such a result implied the crucial role of cholesterol in cold acclimation. Moreover, plenty of differentially expressed genes associated with spliceosomes were enriched during cooling process, which suggested alternative splicing may be involved in the regulation of biological process in acclimation to temperature changes. In researches on extremely low-temperature tolerance, we identified four genes (DUSP1, HSPA6, NR4A1 and GADD45B) associated with MAPK signaling pathway. The four genes, extensively up-regulated at 4°C and remained relatively low expression at moderate temperature, were closely related with extremely cold condition. Further examination of the candidate genes can provide insights into the mechanisms of grass carp to endure extremely low temperature in the winter.