Comprehensive analysis of full-length transcriptomes of Schizothorax prenanti by single-molecule long-read sequencing

Full-length transcriptome analysis of a bloom-forming dinoflagellate Prorocentrum shikokuense (Dinophyceae)

Prorocentrum shikokuense (formerly P. donghaiense) is a pivotal dinoflagellate species associating with the HABs in the East China Sea. The complexity of its large nuclear genome hindered us from understanding its genomic characteristics. Full-length transcriptome sequencing offers a practical solution to decipher the physiological mechanisms of a species without the reference genome. In this study, we employed single-molecule real-time (SMRT) sequencing technology to sequence the full-length transcriptome of Prorocentrum shikokuense. We successfully generated 41.73 Gb of clean SMRT sequencing reads and isolated 105,249 non-redundant full-length non-chimeric reads. Our trial has led to the identification of 11,917 long non-coding RNA transcripts, 514 alternative splicing events, 437 putative transcription factor genes from 17 TF gene families, and 34,723 simple sequence repeats. Additionally, a total of 78,265 open reading frames were identified, of them 15,501 were the protein coding sequences. This dataset is valuable for annotating P. shikokuense genome, and will contribute significantly to the in-depth studies on the molecular mechanisms underlining the dinoflagellate bloom formation.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing

BACKGROUND

Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates.

METHODS

In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species.

RESULTS

A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

PacBio Full-Length Transcriptome of a Tetraploid Sinocyclocheilus multipunctatus Provides Insights into the Evolution of Cavefish

Sinocyclocheilus multipunctatus is a second-class nationally protected wild animal in China. As one of the cavefish, S. multipunctatus has strong adaptability to harsh subterranean environments. In this study, we used PacBio SMRT sequencing technology to generate a first representative full-length transcriptome for S. multipunctatus. Sequence clustering analysis obtained 232,126 full-length transcripts. Among all transcripts, 40,487 were annotated in public databases, while 70,300 microsatellites, 2384 transcription factors, and 16,321 long non-coding RNAs were identified. The phylogenetic tree showed that S. multipunctatus shows a closer relationship to Carassius auratus and Cyprinus carpio, phylogenetically diverging from the common ancestor ~14.74 million years ago (Mya). We also found that between 15.6 and 17.5 Mya, S. multipunctatus also experienced an additional whole-genome duplication (WGD) event, which may have promoted the species evolution of S. multipunctatus. Meanwhile, the overall rates of evolutionary of polyploid S. multipunctatus were significantly higher than those of the other cyprinids, and 220 positively selected genes (PSGs) were identified in two sub-genomes of S. multipunctatus. These PSGs are likely to fulfill critical roles in the process of adapting to diverse cave environments. This study has the potential to facilitate future investigations into the genomic characteristics of S. multipunctatus and provide valuable insights into revealing the evolutionary history of polyploid S. multipunctatus.

Vascular Endothelial Growth Factor A (VEGFA) Regulates Hepatic Lipid and Glycogen Metabolism in Schizothorax prenanti

Vascular endothelial growth factor A (VEGFA) plays important roles in angiogenesis, inflammatory response as well as energy metabolism in mammals. However, its effect on glycolipid metabolism in fish has not been reported. In this study, we cloned and characterized the vegfa gene of Schizothorax prenanti (S. prenanti). vegfa expression was significantly higher in liver and muscle than that in other tissues. Then, the VEGFA recombinant protein was expressed in Escherichia coli and obtained after purification. VEGFA i.p. injection significantly increased the serum glucose and TG content compared with the control group. Moreover, VEGFA protein aggravated the glycogen and lipid deposition in the liver of S. prenanti. In addition, we found that VEGFA treatment increased hepatocyte glycogen and lipid droplet content and increased the levels of pAMPKα (T172). Furthermore, AMPKα inhibition attenuated the ability of VEGFA to induce TG and glycogen accumulation. These results demonstrate that VEGFA regulates hepatic lipid and glycogen metabolism through AMPKα in S. prenanti, which may contribute to a better understanding of VEGFA functions in the glycolipid metabolism of fish.

The sequencing of the key genes and end products in the TLR4 signaling pathway from the kidney of Rana dybowskii exposed to Aeromonas hydrophila

Infectious diseases caused by Aeromonas hydrophila (AH) have reduced the populations of Rana dybowskii). However, little is known about the immune response of R. dybowskii against AH infections. The toll-like receptor (TLR) signaling pathway has been identified as a critical component in innate immunity, responsible for identifying pathogen-associated molecular patterns in pathogens. Our study used the next-generation sequencing technique and single-molecule long-read sequencing to determine the structures of transcript isoforms and functions of genes in the kidneys of R. dybowskii, as well as identify and validate the related genes in the TLR4 signaling pathway. In total, 628,774 reads of inserts were identified, including 300,053 full-length non-chimeric reads and 233,592 non-full-length reads. Among the transcriptome sequences, 124 genes were identified as homologs of known genes in the TLR4 pathway especially inflammatory cytokines and receptors. Our findings shed light on the structures and functions of R. dybowskii genes exposed to AH and confirm the presence of both MyD88-dependent and independent pathways in R. dybowskii. Our work reveals how various functional proteins in amphibians at the initial stage of immune response are activated and complete their corresponding functions in a short time.

Comprehensive transcripts analysis based on single-molecule real-time sequencing and Illumina sequencing provides insights into the mining of Toll-like receptor family in Schizothorax lissolabiatus

Schizothorax lissolabiatus is an economically important cold-water fish species in southwestern China. Because of water pollution and habitat destruction, the number of wild populations has dramatically decreased. In this study, we used PacBio single-molecule real-time (SMRT) sequencing and Illumina sequencing to generate the first full-length transcriptome and transcriptome, respectively. A total of 19 310 polished consensus reads (PC) were obtained, with an average length of 1379 bp and an N50 length of 1485 bp. Meanwhile, 12 253 transcripts were successfully annotated as known homologous genes. The pathway annotation indicated that the enrichment and expression of most genes were mainly related to membrane, signal transduction and binding, and immune response. Furthermore, we identified 16 Toll-like receptors (TLRs) by mining the data from the transcripts. Phylogeny analysis showed that S. lissolabiatus TLR genes (slTLRs) supported the classification of TLRs into six families as in other vertebrates. Selection pressure analyses showed that 16 slTLRs revealed purification selection at the overall evolutionary selection. Further, positive selection signals were still detected in eight slTLRs, and most of the positive selection sites were located in the leucine-rich repeat region (LRR domain) associated with the recognition of pathogenic microorganisms, indicating that the function of these slTLR genes may be affected. Tissue specific expression analysis showed all slTLRs are present in kidney, spleen and liver but the relative expression varied among tissues. In conclusion, this study not only provided a valuable resource of transcripts for further research on S. lissolabiatus, but also contributed to improve the current understanding of the evolutionary history of immune-related genes and the TLR gene family in S. lissolabiatus.

Full-Length Transcriptome Analysis of the Ichthyotoxic Harmful Alga Heterosigma akashiwo (Raphidophyceae) Using Single-Molecule Real-Time Sequencing

The raphidophyte Heterosigma akashiwo is a harmful algal species. The bloom of this organism has been associated with the massive mortality of fish in many coastal waters. To investigate the molecular mechanism of H. akashiwo blooms, having a reliable reference transcriptome of this species is essential. Therefore, in this study, a full-length transcriptome of H. akashiwo was obtained by single-molecule real-time sequencing. In total, 45.44 Gb subread bases were generated, and 16,668 unigenes were obtained after the sequencing data processing. A total of 8666 (52.00%) unigenes were successfully annotated using seven public databases. Among them, mostly phosphorus and nitrogen metabolism genes were detected. Moreover, there were 300 putative transcription factors, 4392 putative long non-coding RNAs, and 7851 simple sequence repeats predicted. This study provides a valuable reference transcriptome for understanding how H. akashiwo blooms at a molecular level.

Full-length transcriptome analysis of the bloom-forming dinoflagellate Akashiwo sanguinea by single-molecule real-time sequencing

The dinoflagellate Akashiwo sanguinea is a harmful algal species and commonly observed in estuarine and coastal waters around the world. Harmful algal blooms (HABs) caused by this species lead to serious environmental impacts in the coastal waters of China since 1998 followed by huge economic losses. However, the full-length transcriptome information of A. sanguinea is still not fully explored, which hampers basic genetic and functional studies. Herein, single-molecule real-time (SMRT) sequencing technology was performed to characterize the full-length transcript in A. sanguinea. Totally, 83.03 Gb SMRT sequencing clean reads were generated, 983,960 circular consensus sequences (CCS) with average lengths of 3,061 bp were obtained, and 81.71% (804,016) of CCS were full-length non-chimeric reads (FLNC). Furthermore, 26,461 contigs were obtained after being corrected with Illumina library sequencing, with 20,037 (75.72%) successfully annotated in the five public databases. A total of 13,441 long non-coding RNA (lncRNA) transcripts, 3,137 alternative splicing (AS) events, 514 putative transcription factors (TFs) members from 23 TF families, and 4,397 simple sequence repeats (SSRs) were predicted, respectively. Our findings provided a sizable insights into gene sequence characteristics of A. sanguinea, which can be used as a reference sequence resource for A. sanguinea draft genome annotation, and will contribute to further molecular biology research on this harmful bloom algae.

Dynamic Transcriptomic Profiling During Liver Development in Schizothorax Prenanti

Liver is an important organ for glucose and lipid metabolism, immunity, and detoxification in fish. However, the gene regulatory network of postnatal liver development still remains unknown in teleost fish. In this study, we performed transcriptome analysis on the liver of S. prenanti at three stages. A total of 1692 differentially expressed genes (DGEs) were identified across three liver developmental stages. The oil red O staining and PAS staining revealed that the lipid content of liver was increased and the glycogen content of liver was decreased during liver development. The fatty acids biosynthesis related genes were upregulated in adult and young stages compared with juvenile stage, while lipid degradation related genes were downregulated. The genes related to glycolysis, gluconeogenesis and glycogenolysis were upregulated in juvenile or young stages compared with adult stage. Further pathway analysis indicated that the CYP450 pathway, cell cycle and amino acid metabolic pathway were induced in the process of liver maturation. Our study presents the gene expression pattern in different liver development stages of S. prenanti and may guide future studies on metabolism of S. prenanti liver.

Molecular and functional characterization of the retinol-binding protein 4 (RBP4) in hepatocytes of Schizothorax prenanti in response to palmitic acid

Retinol-binding protein 4 (RBP4) protein is a kind of adipokines synthesized and secreted by the liver, which has been verified to play important roles in liver metabolism and energy homeostasis. However, the effects of RBP4 on hepatic lipid accumulation are still elusive in fish. In the present study, we cloned and characterized the RBP4 gene in Schizothorax prenanti (S. prenanti). RBP4 gene was specifically expressed in the liver and abdominal adipose tissue. Palmitic acid (PA; 400 μM) can significantly increase lipid deposition in primary hepatocytes after 12 h of treatment. Furthermore, RBP4 knockdown can relieve the excessive lipid deposition and endoplasmic reticulum stress in the hepatocytes caused by PA. The inhibition of RBP4 abolished the ability of PA to induce the expression of genes involved in lipogenesis and endoplasmic reticulum stress. These results demonstrate that RBP4 inhibition attenuated PA-induced lipid deposition and endoplasmic reticulum stress in hepatocytes of S. prenanti. This study could contribute to improve the understanding of RBP4 functions in the PA-induced lipid deposition in hepatocytes of fish.

Decoding the roles of extremophilic microbes in the anaerobic environments: Past, Present, and Future

•

The inaccessible extreme environments harbor a large majority of anaerobic microbes which remain unknown.

•

Anaerobic microbes are used in a variety of industrial applications.

•

In the future, metagenomic-assisted techniques can be used to identify novel anaerobic microbes from the unexplored extreme environments.

•

Genetic engineering can be used to enhance the efficiency of anaerobic microbes for various processes.

The genome of an organism is directly or indirectly correlated with its environment. Consequently, different microbes have evolved to survive and sustain themselves in a variety of environments, including unusual anaerobic environments. It is believed that their genetic material could have played an important role in the early evolution of their existence in the past. Presently, out of the uncountable number of microbes found in different ecosystems we have been able to discover only one percent of the total communities. A large majority of the microbial populations exists in the most unusual and extreme environments. For instance, many anaerobic bacteria are found in the gastrointestinal tract of humans, soil, and hydrothermal vents. The recent advancements in Metagenomics and Next Generation Sequencing technologies have improved the understanding of their roles in these environments. Presently, anaerobic bacteria are used in various industries associated with biofuels, fermentation, production of enzymes, vaccines, vitamins, and dairy products. This broad applicability brings focus to the significant contribution of their genomes in these functions. Although the anaerobic microbes have become an irreplaceable component of our lives, a major and important section of such anaerobic microbes still remain unexplored. Therefore, it can be said that unlocking the role of the microbial genomes of the anaerobes can be a noteworthy discovery not just for mankind but for the entire biosystem as well.

Expression Profile Analysis of Long Non-coding RNA in OVX Models-Derived BMSCs for Postmenopausal Osteoporosis by RNA Sequencing and Bioinformatics

Osteoporosis (OP) has the characteristics of a systematically impaired bone mass, strength, and microstructure. Long non-coding RNAs (lncRNAs) are longer than 200 nt, and their functions in osteoporosis is yet not completely understood. We first harvested the bone marrow mesenchymal stem cells (BMSCs) from ovariectomy (OVX) and sham mice. Then, we systematically analyzed the differential expressions of lncRNAs and messenger RNAs (mRNAs) and constructed lncRNA–mRNA coexpression network in order to identify the function of lncRNA in osteoporosis. Totally, we screened 743 lncRNAs (461 upregulated lncRNAs and 282 downregulated lncRNAs) and 240 mRNAs (128 upregulated and 112 downregulated) with significantly differential expressions in OP compared to normal. We conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses to investigate the functions and pathways of the differential expression of messenger RNAs (mRNAs), a coexpressed network of lncRNA/mRNA. Quantitative PCR (qPCR) validated that the expressions of NONMMUT096150.1, NONMMUT083450.1, and NONMMUT029743.2 were all downregulated, whereas NONMMUT026970.2, NONMMUT051734.2, NONMMUT003617.2, and NONMMUT034049.2 were all upregulated in the OVX group. NONMMUT096150.1, as a key lncRNA in OP, was identified to modulate the adipogenesis of BMSCs. Further analysis suggested that NONMMUT096150.1 might modulate the adipogenesis of BMSCs via the peroxisome proliferator-activated receptor (PPAR) signaling pathway, AMPK signaling pathway, and the lipolysis regulation in adipocyte and adipocytokine signaling pathway. Our study expands the understanding of lncRNA in the pathogenesis of OP.