Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Genome assembly, Full-length transcriptome, and isoform diversity of Red Snapper, Lutjanus argentimaculatus

The mangrove red snapper, Lutjanus argentimaculatus, is a marine food fish of economic and aquaculture importance. The application of genomic selection-based breeding programs for this species is limited by the absence of a reference genome and transcriptome profiles. The current study attempted to fill this void by generating genomic and transcriptomic resources for red snapper. Using PacBio long reads, and Arima Hi-C linked reads, a scaffold-level genome assembly was generated for L. argentimaculatus. The assembly is of 1.03 Gb comprising of 400 scaffolds with N50 of 33.8 Mb and was assessed to be 97.2% complete upon benchmarking with BUSCO. Full-length transcriptome generated with PacBio Iso-Sequencing strategy using six tissues (muscle, gills, liver, kidney, stomach, and gonad) contained 56,515 isoforms belonging to 18,108 unique genes with N50 length of 3,973 bp. The resources generated will have potential applications in the functional studies, conservation, broodstock management and selective breeding programmes of L. argentimaculatus.

Identification, diversity, and evolution analysis of thioester-containing protein family in Pacific oyster (Crassostrea gigas) and immune response to biotic and abiotic stresses

Thioester-containing proteins (TEPs) play a vital role in the innate immune response to biotic and abiotic stresses. In this study, the TEPs in C. gigas were identified, and their gene structure, phylogenetic relationships, collinearity relationships, expression profiles, sequence diversity, and alternative splicing were analyzed. Eight Tep genes were identified in C. gigas genome. Functional analysis and evolutionary relationships indicated a high level of homology to other mollusks TEPs. The transcriptome quantitative analysis results showed that the Tep genes in C. gigas respond to heat stress and Vibrio stress. Alternative splicing analysis revealed four Tep genes (designated A2M_1, CD109_3, CD109_5, complement C3) encode multiple alternative splice variants. Analysis of gene structure and multiple alignments revealed that seven CD109_5 variants are produced through the alternative splicing of the 19th exon, which encodes the highly variable central region. Sequence diversity analysis revealed thirteen missense variants within the 19th exon region of these seven CD109_5 alternative splice variants. Furthermore, the differential alternative splicing analysis showed significant induction of CD109_5, A2M_1 and A2M_2 variants after infection with V. parahaemolyticus. This study explores the Tep genes of C. gigas, providing insights into the molecular mechanisms underlying the involvement of C. gigas TEPs in innate immunity.

Transcriptome analyses reveal key roles of alternative splicing regulation in atlantic salmon during the infectious process of Piscirickettsiosis disease

In the Chilean salmon farming industry, infection by Piscirickettsia salmonis is the primary cause of the main bacterial disease known as Piscirickettsiosis, which has an overwhelming economic impact. Although it has been demonstrated that Piscirickettsiosis modifies the expression of numerous salmonids genes, it is yet unknown how alternative splicing (AS) contributes to salmonids bacterial infection. AS, has the potential to create heterogeneity at the protein and RNA levels and has been associated as a relevant molecular mechanism in the immune response of eukaryotes to several diseases. In this study, we used RNA data to survey P. salmonis-induced modifications in the AS of Atlantic salmon and found that P. salmonis infection promoted a substantial number (158,668) of AS events. Differentially spliced genes (DSG) sensitive to Piscirickettsiosis were predominantly enriched in genes involved in RNA processing, splicing and spliceosome processes (e.g., hnRNPm, hnRPc, SRSF7, SRSF45), whereas among the DSG of resistant and susceptible to Piscirickettsiosis, several metabolic and immune processes were found, most notably associated to the regulation of GTPase, lysosome and telomere organization-maintenance. Furthermore, we found that DSG were mostly not differentially expressed (5-7 %) and were implicated in distinct biological pathways. Therefore, our results underpin AS achieving a significant regulatory performance in the response of salmonids to Piscirickettsiosis.

Physiological, metabolomic, and transcriptomic reveal metabolic pathway alterations in Gymnocypris przewalskii due to cold exposure

Teleost fish have evolved various adaptations that allow them to tolerate cold water conditions. However, the underlying mechanism of this adaptation is poorly understood in Tibetan Plateau fish. RNA-seq combined with liquid chromatography‒mass spectrometry (LC‒MS/MS) metabolomics was used to investigate the physiological responses of a Tibetan Plateau-specific teleost, Gymnocypris przewalskii, under cold conditions. The 8-month G. przewalskii juvenile fish were exposed to cold (4 ℃, cold acclimation, CA) and warm (17 ℃, normal temperature, NT) temperature water for 15 days. Then, the transcript profiles of eight tissues, including the brain, gill, heart, intestine, hepatopancreas, kidney, muscle, and skin, were evaluated by transcriptome sequencing. The metabolites of the intestine, hepatopancreas, and muscle were identified by LC‒MS/MS. A total of 5,745 differentially expressed genes (DEGs) were obtained in the CA group. The key DEGs were annotated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The DEGs from the eight tissues were significantly enriched in spliceosome pathways, indicating that activated alternative splicing is a critical biological process that occurs in the tissues to help fish cope with cold stress. Additionally, 82, 97, and 66 differentially expressed metabolites were identified in the intestine, hepatopancreas, and muscle, respectively. Glutathione metabolism was the only overlapping significant pathway between the transcriptome and metabolome analyses in these three tissues, indicating that an activated antioxidative process was triggered during cold stress. In combination with the multitissue transcriptome and metabolome, we established a physiology-gene‒metabolite interaction network related to energy metabolism during cold stress and found that gluconeogenesis and long-chain fatty acid metabolism played critical roles in glucose homeostasis and energy supply.

The Responses of Alternative Splicing during Heat Stress in the Pacific White Shrimp Litopenaeus vannamei

Heat tolerance is increasingly becoming a crucial trait for aquaculture species in the face of rapidly changing climate conditions. Alternative splicing (AS) is a vital mechanism within cells that modulates gene abundance and functional diversity, enabling organisms to effectively respond to diverse stressful conditions, including thermal stress. However, it is still uncertain whether AS contributes to heat tolerance in shrimp. In this study, we conducted an extensive transcriptome analysis on the Pacific white shrimp, Litopenaeus vannamei, revealing a total of 1267, 987, and 130 differential AS events (DAS) in the gill, hepatopancreas, and muscle, respectively, following exposure to heat stress. Among all of the DAS events, exon skipping (ES) was the predominant form of splicing modification observed. Interestingly, a minor portion of DAS genes exhibited overlap across the three tissues, implying that heat stress exerts unique effects on various tissue types. Moreover, the functional enrichment analysis demonstrated that commonly identified DAS genes were primarily associated with the "spliceosome" pathway, indicating that the AS of splicing-related genes played a crucial role in the response to heat stress. Our findings also revealed that heat stress tended to induce longer mRNA isoforms through differential alternative 3' splice site (A3SS) events. Notably, A3SS events exhibited the highest proportion of maintained open reading frames (ORFs) under heat stress. Interestingly, we observed a limited overlap between the genes exhibiting DAS and those showing differential gene expression (DEG), indicating that AS may function as a distinct regulatory mechanism independent of transcriptional regulation in response to heat stress. This is the first comprehensive study on AS in crustacea species under heat stress, which broadens our understanding of the regulatory mechanisms governing the crustaceans' response to environmental stress, providing valuable insights for the aquaculture breeding of shrimp and other aquatic animals.

Unveiling the Role of Dynamic Alternative Splicing Modulation After Infestation with Sea Lice (Caligus rogercresseyi) in Atlantic Salmon

Sea lice are pathogenic marine ectoparasite copepods that represent a severe risk to the worldwide salmon industry. Several transcriptomic investigations have characterized the regulation of gene expression response of Atlantic salmon to sea lice infestation. These studies have focused on the levels of transcript, overlooking the potentially relevant role of alternative splicing (AS), which corresponds to an essential control mechanism of gene expression through RNA processing. In the present study, we performed a genome-wide bioinformatics characterization of differential AS event dynamics in control and infested C. rogercresseyi Atlantic salmon and in resistant and susceptible phenotypes. We identified a significant rise of alternative splicing events and AS genes after infestation and 176 differential alternative splicing events (DASE) from 133 genes. In addition, a higher number of DASE and AS genes were observed among resistant and susceptible phenotypes. Functional annotation of AS genes shows several terms and pathways associated with behavior, RNA splicing, immune response, and RNA binding. Furthermore, three protein-coding genes were identified undergoing differential transcript usage events, among resistant and susceptible phenotypes. Our findings support AS performing a relevant regulatory role in the response of salmonids to sea lice infestation.

Possible Contribution of Alternative Transcript Isoforms in Mature Biofilm Growth Phase of Candida glabrata

Expression of genome-wide alternative transcript isoforms and differential transcript isoform usage in different biological conditions (isoform switching) are responsible for the varied proteomic functional diversity in higher eukaryotic organisms. However, these mechanisms have not been studied in Candida glabrata, which is a potent eukaryotic opportunistic pathogen. Biofilm formation is an important virulence factor of C. glabrata that withstands antifungal drug stress and overcomes the host-immune response. Here, we present the genome-wide differential transcript isoform expression (DTE) and differential transcript isoform usage (DTU) in a mature biofilm growth phase of C. glabrata (clinical isolate; NCCPF 100,037) using the RNA sequencing approach. The DTE analysis generated 7837 transcript isoforms from the C. glabrata genome (5293 genes in total), and revealed that transcript isoforms generated from 292 genes showed significant DTU in the mature biofilm cells. Gene ontology, pathway analysis and protein-protein interactions of significant transcript isoforms, further substantiated that their specific expression and differential usage is required for transitioning the planktonic cells to biofilm in C. glabrata. The present study reported the possible role of expression of alternative transcript isoforms and differential transcript isoform usage in the mature biofilms of C. glabrata. The observation derived from the study may prove to be beneficial for making future antifungal therapeutic strategies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-022-01036-7.

High-Altitude Stress Orchestrates mRNA Expression and Alternative Splicing of Ovarian Follicle Development Genes in Tibetan Sheep

High-altitude stress threatens the survival rate of Tibetan sheep and reduces their fertility. However, the molecular basis of this phenomenon remains elusive. Here, we used RNA-seq to elucidate the transcriptome dynamics of high-altitude stress in Tibetan sheep ovaries. In total, 104 genes were characterized as high-altitude stress-related differentially expressed genes (DEGs). In addition, 36 DEGs contributed to ovarian follicle development, and 28 of them were downregulated under high-altitude stress. In particular, high-altitude stress significantly suppressed the expression of two ovarian lymphatic system marker genes: LYVE1 and ADAMTS-1. Network analysis revealed that luteinizing hormone (LH)/follicle-stimulating hormone (FSH) signaling-related genes, such as EGR1, FKBP5, DUSP1, and FOS, were central regulators in the DEG network, and these genes were also suppressed under high-altitude stress. As a post-transcriptional regulation mechanism, alternative splicing (AS) is ubiquitous in Tibetan sheep. High-altitude stress induced 917 differentially alternative splicing (DAS) events. High-altitude stress modulated DAS in an AS-type-specific manner: suppressing skipped exon events but increasing retained intron events. C2H2-type zinc finger transcription factors and RNA processing factors were mainly enriched in DAS. These findings revealed high-altitude stress repressed ovarian development by suppressing the gene expression of LH/FSH hormone signaling genes and inducing intron retention of C2H2-type zinc finger transcription factors.

Understanding the Function and Mechanism of Zebrafish Tmem39b in Regulating Cold Resistance

Autophagy and endoplasmic reticulum (ER) stress response are among the key pathways regulating cold resistance of fish through eliminating damaged cellular components and facilitating the restoration of cell homeostasis upon exposure to acute cold stress. The transmembrane protein 39A (TMEM39A) was reported to regulate both autophagy and ER stress response, but its vertebrate-specific paralog, the transmembrane protein 39B (TMEM39B), has not been characterized. In the current study, we generate tmem39b-knockout zebrafish lines and characterize their survival ability under acute cold stress. We observed that the dysfunction of Tmem39b remarkably decreased the cold resilience of both the larval and adult zebrafish. Gene transcription in the larvae exposed to cold stress and rewarming were characterized by RNA sequencing (RNA-seq) to explore the mechanisms underlying functions of Tmem39b in regulating cold resistance. The results indicate that the deficiency of Tmem39b attenuates the up-regulation of both cold- and rewarming-induced genes. The cold-induced transcription factor genes bif1.2, fosab, and egr1, and the rewarming-activated immune genes c3a.3, il11a, and sting1 are the representatives influenced by Tmem39b dysfunction. However, the loss of tmem39b has little effect on the transcription of the ER stress response- and autophagy-related genes. The measurements of the phosphorylated H2A histone family member X (at Ser 139, abbreviated as γH2AX) demonstrate that zebrafish Tmem39b protects the cells against DNA damage caused by exposure to the cold-warming stress and facilitates tissue damage repair during the recovery phase. The gene modules underlying the functions of Tmem39b in zebrafish are highly enriched in biological processes associated with immune response. The dysfunction of Tmem39b also attenuates the up-regulation of tissue C-reactive protein (CRP) content upon rewarming. Together, our data shed new light on the function and mechanism of Tmem39b in regulating the cold resistance of fish.

Mitochondria Dysfunction and Cell Apoptosis Limit Resistance of Nile Tilapia (Oreochromis niloticus) to Lethal Cold Stress

Simple Summary

Sensitivity of Nile tilapia (Oreochromis niloticus) to cold stress represents a major concern for both aquaculture and theoretical study; however, the cellular and molecular mechanisms determining cold susceptibility of it remain largely unknown. In this study, we first estimated the median survival time of juvenile Nile tilapia under exposure to lethal cold stress (12 °C). The fish were classified as cold-sensitive or cold-tolerant based on their behavioral manifestation after exposed to 12 °C for 3 days. Subsequently, histological, biochemical and gene expression analyses were performed for the fish with different cold resistance to explore the cellular and molecular events underlying cold susceptibility of Nile tilapia. We found that exposure of Nile tilapia to lethal cold stress caused systemic tissue structure changes, mitochondrial swelling and dysfunction, induction of apoptosis and endoplasmic reticulum (ER) stress-related genes and cell apoptosis. The extent of these adverse cellular and molecular events determines an individual’s ability to survive cold stress. Our data indicate that mitochondria dysfunction and mitochondria-mediated cell apoptosis are the main factors limiting Nile tilapia’s cold resistance.

Abstract

Inability of Nile tilapia (Oreochromis niloticus) to withstand cold stress represents a major economic concern, which restricts the culture area, limits the growing period and even results in mass mortality in cold seasons. However, the cellular and molecular mechanisms determining cold susceptibility of Nile tilapia remain largely unknown. In this study, we characterized the ability of juvenile Nile tilapia to survive lethal cold stress (12 °C) and the median survival time (LT50) of the experimental fish under exposure to 12 °C cold stress was estimated as 3.14 d. After being exposed to 12 °C for 3 d, the survivors that lost equilibrium (LE) and those that swam normally (NO) were regarded as cold-sensitive and cold-tolerant, respectively. The untreated (Ctrl), NO and LE fish were subjected to histological, biochemical and gene expression analyses to explore the cellular and molecular events underlying cold susceptibility of Nile tilapia. Exposure of Nile tilapia to lethal cold stress caused systemic tissue structure changes, mitochondrial swelling and dysfunction, induction of apoptosis and endoplasmic reticulum (ER) stress-related genes and cell apoptosis. The extent of these adverse cellular and molecular events determines an individual’s ability to survive cold stress. Our data indicate that mitochondria dysfunction and mitochondria-mediated cell apoptosis are the main factors limiting Nile tilapia’s cold resistance.

Tissue-Specific Analysis of Alternative Splicing Events and Differential Isoform Expression in Large Yellow Croaker (Larimichthys crocea) After Cryptocaryon irritans Infection

The large yellow croaker (Larimichthys crocea) is one of the most important mariculture fish in China. Recently, cryptocaryonosis caused by Cryptocryon irritans infection has brought huge economic losses and threatens the healthy and sustainable development of the L. crocea industry. However, the molecular mechanism and regulation process for L. crocea resistance to C. irritans infection has not been fully researched. Alternative splicing (AS) is an important post-transcriptional regulatory mechanism that allows cells to produce transcriptional and proteomic diversity. The results of AS are tissue dependent, and the expression of tissue-specific transcription subtype genes is determined by AS and transcriptional regulation. However, studies on the tissue specificity of AS events in L. crocea following infection with C. irritans have not been performed. In this study, the L. crocea were artificially infected with C. irritans; their skin and gill were collected at 0 h, 24 h, 48 h, 72 h, and 96 h post infection. After sequencing and differential expression analysis, a set of 452, 692, 934, 711, 534, and 297 differential alternative splicing (DAS) events were identified in 0 h, 12 h, 24 h, 48 h, 72 h, and 96 h post infection respectively. Furthermore, 4160 differentially expressed isoforms (DEIs) and 4209 DEI genes were identified from all time point groups. GO enrichment and pathway analysis indicated that many genes of DAS and DEIs were rich in immune-related GO terms and KEGG pathways, such as the Toll and Imd signaling pathway, NOD-like receptor signaling pathway, TNF signaling pathway, and TNF signaling pathway. Among hub DEI genes, alternative splicing-related genes (cwc25, prpf8, and sf3a3), skin function-related gene (fa2h), and oxygen deprivation-related gene (hyo1) were found in DEI genes. This study provided insight into the temporal change of DAS and DEIs between skin and gill of L. crocea against C. irritans infection and revealed that these differences might play immune-related roles in the infection process.

Non-targeted Metabolomics Analysis Based on LC–MS to Assess the Effects of Different Cold Exposure Times on Piglets

Pigs are susceptible to low temperature conditions, and cold stress causes metabolic changes in the body to increase heat production as an adaption to adverse environments. To characterize and validate different metabolites in piglet livers at different cold exposure times, sixteen 30-day-old male weaned piglets with similar weights were randomly divided into four groups: the normal temperature group (24 ± 2°C, NT) and cold exposure (4 ± 2°C) 2-h group (CS2), 6-h group (CS6), and 12-h group (CS12). At the end of the experiment, the liver samples were analyzed using systemic non-targeted metabolomics. Eight known differentially abundant metabolites (farnesyl pyrophosphate, isocitrate, triethanolamine, phenylethylamine, deoxynosine, citric acid, maltotriose, and epinephrine) were observed between the CS groups and the control group in positive and negative ion modes. The eight main differentially abundant metabolites involved in seven metabolite classifications. Metabolic pathways and enrichment analyses revealed that the pathways involved three KEGG pathway classifications. Most of the pathways were related to amino acid or energy metabolism. Moreover, the metabolic pathways were not identical under different cold exposure times, with those following 2 and 6 h of cold exposure more related to carbohydrates and energy production and those following 12 h of cold exposure more related to the metabolism connected with epinephrine. Thus, under different cold exposure times, the metabolite profiles and metabolic pathways differed.

The alternative splicing landscape of a coral reef fish during a marine heatwave

Alternative splicing is a molecular mechanism that enables a single gene to encode multiple transcripts and proteins by post-transcriptional modification of pre-RNA molecules. Changes in the splicing scheme of genes can lead to modifications of the transcriptome and the proteome. This mechanism can enable organisms to respond to environmental fluctuations. In this study, we investigated patterns of alternative splicing in the liver of the coral reef fish Acanthochromis polyacanthus in response to the 2016 marine heatwave on the Great Barrier Reef. The differentially spliced (DS; n = 40) genes during the onset of the heatwave (i.e., 29.49°C or +1°C from average) were related to essential cellular functions such as the MAPK signaling system, Ca(2+) binding, and homeostasis. With the persistence of the heatwave for a period of one month (February to March), 21 DS genes were detected, suggesting that acute warming during the onset of the heatwave is more influential on alternative splicing than the continued exposure to elevated temperatures. After the heatwave, the water temperature cooled to ~24.96°C, and fish showed differential splicing of genes related to cyto-protection and post-damage recovery (n = 26). Two-thirds of the DS genes detected across the heatwave were also differentially expressed, revealing that the two molecular mechanisms act together in A. polyacanthus to cope with the acute thermal change. This study exemplifies how splicing patterns of a coral reef fish can be modified by marine heatwaves. Alternative splicing could therefore be a potential mechanism to adjust cellular physiological states under thermal stress and aid coral reef fishes in their response to more frequent acute thermal fluctuations in upcoming decades.

Genome-Wide Analysis of Alternative Splicing (AS) Mechanism Provides Insights into Salinity Adaptation in the Livers of Three Euryhaline Teleosts, including Scophthalmus maximus, Cynoglossus semilaevis and Oncorhynchus mykiss

Simple Summary

Alternative splicing (AS) is a key post-transcriptional regulatory mechanism that acts an important regulator in response to environmental stimuli in organisms. In the present study, 18 RNA-Seq datasets were utilized to investigate the potential roles of AS in response to different salinity environments in the livers of three euryhaline teleosts, including turbot (Scophthalmus maximus), tongue sole (Cynoglossus semilaevis) and steelhead trout (Oncorhynchus mykiss). The results indicated that different salinity environments changed the splicing patterns of numerous RNA splicing regulators, which might affect the splicing decisions of many downstream target genes in response to salinity changes. This study provides preliminary evidence for the important roles of AS events in salinity adaptation in teleosts.

Abstract

Salinity is an important environmental factor that directly affects the survival of aquatic organisms, including fish. However, the underlying molecular mechanism of salinity adaptation at post-transcriptional regulation levels is still poorly understood in fish. In the present study, 18 RNA-Seq datasets were utilized to investigate the potential roles of alternative splicing (AS) in response to different salinity environments in the livers of three euryhaline teleosts, including turbot (Scophthalmus maximus), tongue sole (Cynoglossus semilaevis) and steelhead trout (Oncorhynchus mykiss). A total of 10,826, 10,741 and 10,112 AS events were identified in the livers of the three species. The characteristics of these AS events were systematically investigated. Furthermore, a total of 940, 590 and 553 differentially alternative splicing (DAS) events were determined and characterized in the livers of turbot, tongue sole and steelhead trout, respectively, between low- and high-salinity environments. Functional enrichment analysis indicated that these DAS genes in the livers of three species were commonly enriched in some GO terms and KEGG pathways associated with RNA processing. The most common DAS genes work as RNA-binding proteins and play crucial roles in the regulation of RNA splicing. The study provides new insights into uncovering the molecular mechanisms of salinity adaptation in teleosts.

RNA binding protein NKAP protects glioblastoma cells from ferroptosis by promoting SLC7A11 mRNA splicing in an m6A-dependent manner

Ferroptosis is a form of cell death characterized by lipid peroxidation. Previous studies have reported that knockout of NF-κB activating protein (NKAP), an RNA-binding protein, increased lipid peroxidation level in naive T cells and induced cell death in colon cancer cells. However, there was no literature reported the relationship between NKAP and ferroptosis in glioblastoma cells. Notably, the mechanism of NKAP modulating ferroptosis is still unknown. Here, we found NKAP knockdown induced cell death in glioblastoma cells. Silencing NKAP increased the cell sensitivity to ferroptosis inducers both in vitro and in vivo. Exogenous overexpression of NKAP promoted cell resistance to ferroptosis inducers by positively regulating a ferroptosis defense protein, namely cystine/glutamate antiporter (SLC7A11). The regulation of SLC7A11 by NKAP can be weakened by the m⁶A methylation inhibitor cycloleucine and knockdown of the m⁶A writer METTL3. NKAP combined the “RGAC” motif which was exactly in line with the m⁶A motif “RGACH” (R = A/G, H = A/U/C) uncovered by the m⁶A-sequence. RNA Immunoprecipitation (RIP) and Co-Immunoprecipitation (Co-IP) proved the interaction between NKAP and m⁶A on SLC7A11 transcript. Following its binding to m⁶A, NKAP recruited the splicing factor proline and glutamine-rich (SFPQ) to recognize the splice site and then conducted transcription termination site (TTS) splicing event on SLC7A11 transcript and the retention of the last exon, screened by RNA-sequence and Mass Spectrometry (MS). In conclusion, NKAP acted as a new ferroptosis suppressor by binding to m⁶A and then promoting SLC7A11 mRNA splicing and maturation.

Development of Disease-Resistance-Associated Microsatellite DNA Markers for Selective Breeding of Tilapia (Oreochromis spp.) Farmed in Taiwan

There are numerous means to improve the tilapia aquaculture industry, and one is to develop disease resistance through selective breeding using molecular markers. In this study, 11 disease-resistance-associated microsatellite markers including 3 markers linked to hamp2, 4 linked to hamp1, 1 linked to pgrn2, 2 linked to pgrn1, and 1 linked to piscidin 4 (TP4) genes were established for tilapia strains farmed in Taiwan after challenge with Streptococcus inae. The correlation analysis of genotypes and survival revealed a total of 55 genotypes related to survival by the chi-square and Z-test. Although fewer markers were found in B and N2 strains compared with A strain, they performed well in terms of disease resistance. It suggested that this may be due to the low potency of some genotypes and the combinatorial arrangement between them. Therefore, a predictive model was built by the genotypes of the parental generation and the mortality rate of different combinations was calculated. The results show the same trend of predicted mortality in the offspring of three new disease-resistant strains as in the challenge experiment. The present findings is a nonkilling method without requiring the selection by challenge with bacteria or viruses and might increase the possibility of utilization of selective breeding using SSR markers in farms.

Characterization of Biological Pathways Regulating Acute Cold Resistance of Zebrafish

Low temperature stress represents a major threat to the lives of both farmed and wild fish species. However, biological pathways determining the development of cold resistance in fish remain largely unknown. Zebrafish larvae at 96 hpf were exposed to lethal cold stress (10 °C) for different time periods to evaluate the adverse effects at organism, tissue and cell levels. Time series RNA sequencing (RNA-seq) experiments were performed to delineate the transcriptomic landscape of zebrafish larvae under cold stress and during the subsequent rewarming phase. The genes regulated by cold stress were characterized by progressively enhanced or decreased expression, whereas the genes associated with rewarming were characterized by rapid upregulation upon return to normal temperature (28 °C). Genes such as trib3, dusp5 and otud1 were identified as the representative molecular markers of cold-induced damages through network analysis. Biological pathways involved in cold stress responses were mined from the transcriptomic data and their functions in regulating cold resistance were validated using specific inhibitors. The autophagy, FoxO and MAPK (mitogen-activated protein kinase) signaling pathways were revealed to be survival pathways for enhancing cold resistance, while apoptosis and necroptosis were the death pathways responsible for cold-induced mortality. Functional mechanisms of the survival-enhancing factors Foxo1, ERK (extracellular signal-regulated kinase) and p38 MAPK were further characterized by inhibiting their activities upon cold stress and analyzing gene expression though RNA-seq. These factors were demonstrated to determine the cold resistance of zebrafish through regulating apoptosis and p53 signaling pathway. These findings have provided novel insights into the stress responses elicited by lethal cold and shed new light on the molecular mechanisms underlying cold resistance of fish.

Characterization and Comparative Transcriptomic Analysis of Skeletal Muscle in Pekin Duck at Different Growth Stages Using RNA-Seq

Simple Summary

Skeletal muscle is an important tissue and its development is strictly regulated by genes. In this study, in order to understand the muscle-related gene expression in Pekin duck, RNA-seq was performed to analyze and compare skeletal muscle at different growth stages. Alternative splicing, single nucleotide polymorphisms and insertion–deletions were detected, and 299 novel genes were discovered. MYL4, IGF2BP1, CSRP3, SPP1, KLHL31, LAMB2, LAMA2, ITGB1 and OPN played crucial roles in skeletal muscle development. Oxidative phosphorylation, ECM-receptor interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids participated in the regulation of skeletal muscle development in Pekin duck. This study provides an important reference for revealing the developmental mechanisms of pectoral and leg muscles in duck.

Abstract

Skeletal muscle, accounting for approximately 50% of body weight, is the largest and most important tissue. In this study, the gene expression profiles and pathways in skeletal muscle of Pekin duck were investigated and compared at embryonic day 17, 21, and 27 and postnatally at 6 months of age. An average of 49,555,936 reads in each sample was obtained from the transcriptome libraries. Over 70.0% of alternative splicing (AS) in each sample was mainly alternative 5′ first exon (transcription start site)—the first exon splicing (TSS) and alternative 3′ last exon (transcription terminal site)—the last exon splicing (TTS), indicating that TSS and TTS were the most common AS event in Pekin ducks, and these AS events were closely related to the regulation of muscle development at different growth stages. The results provided a valuable genomic resource for selective breeding and functional studies of genes. A total of 299 novel genes with ≥2 exons were obtained. There were 294 to 2806 differentially expressed genes (DEGs) in each pairwise comparison of Pekin duck. Notably, 90 DEGs in breast muscle and 9 DEGs in leg muscle were co-expressed at all developmental points. DEGs were validated by qPCR analysis, which confirmed the tendency of the expression. DEGs related to muscle development were involved in biological processes such as “endodermal cell differentiation”, “muscle cell cellular homeostasis”, “skeletal muscle tissue growth” and “skeletal muscle cell differentiation”, and were involved in pathways such as oxidative phosphorylation, ECM-receptor (extracellular matrix receptor) interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids. Some DEGs, including MYL4, IGF2BP1, CSRP3, SPP1 and KLHL31, as well as LAMB2, LAMA2, ITGB1 and OPN, played crucial roles in muscle growth and development. This study provides valuable information about the expression profile of mRNAs and pathways from duck skeletal muscle at different growth stages, and further functional study of these mRNAs and pathways could provide new ideas for studying the molecular networks of growth and development in duck skeletal muscle.

Establishment and characterization of a cold-sensitive neural cell line from the brain of tilapia (Oreochromis niloticus)

The aquaculture of tilapia (Oreochromis sp.) is adversely affected by the sensitivity to cold stress. A large number of genes in tilapia were found to be regulated by cold stress, but their functions and mechanisms in cold tolerance remain largely unknown, partially due to the lack of a suitable in vitro model. An immortal neural cell line designated as tilapia brain neural (TBN) was established from brain tissue of the genetically improved farmed tilapia strain of Nile tilapia (Oreochromis niloticus). The TBN cells show a neuron-like morphology at low density and form a fibroblast-like monolayer at high density. Transcriptome profiling through RNA-sequencing revealed that a total of 15,011 genes were expressed in the TBN cells. The TBN cells express a wide array of marker genes for neural cells. A comparative analysis of the featured genes among the 17 cell clusters isolated from the subventricular zone of mouse brain revealed the highest transcriptome similarity between the TBN cells and the transient amplifying progenitors (TAPs). The TBN cells tolerate relatively high culture temperatures, and the highest growth rate was observed for the cells cultured at 32°C compared with those at 30°C, 28°C and 26°C. Nonetheless, this cell line is cold sensitive. Exposure of the cells to 16°C or lower temperatures significantly decreased cell confluences and induced apoptosis. The TBN cells were more sensitive to cold stress than the ZF4 cells (embryonic zebrafish fibroblasts). Moreover, the TBN cells can be efficiently transfected through electroporation. This study provides an invaluable research tool to understand the nature of cold sensitivity of tilapia and to dissect the function and mechanism of genes in regulating cold tolerance of fish.

Molecular link between circadian clocks and cardiac function: a network of core clock, slave clock, and effectors

The circadian rhythm has a strong influence on both cardiac physiology and disease in humans. Preclinical studies primarily using tissue-specific transgenic mouse models have contributed to our understanding of the molecular mechanism of the circadian clock in the cardiovascular system. The core clock driven by CLOCK:BMAL1 complex functions as a universal timing machinery that primarily sets the pace in all mammalian cell types. In one specific cell or tissue type, core clock may control a secondary transcriptional oscillator, conceptualized as slave clock, which confers the oscillatory expression of tissue-specific effectors. Here, we discuss a core clock-slave clock-effectors network, which links the molecular clock to cardiac function.

Temperature-Dependent Alternative Splicing of Precursor mRNAs and Its Biological Significance: A Review Focused on Post-Transcriptional Regulation of a Cold Shock Protein Gene in Hibernating Mammals

Multiple mRNA isoforms are often generated during processing such as alternative splicing of precursor mRNAs (pre-mRNA), resulting in a diversity of generated proteins. Alternative splicing is an essential mechanism for the functional complexity of eukaryotes. Temperature, which is involved in all life activities at various levels, is one of regulatory factors for controlling patterns of alternative splicing. Temperature-dependent alternative splicing is associated with various phenotypes such as flowering and circadian clock in plants and sex determination in poikilothermic animals. In some specific situations, temperature-dependent alternative splicing can be evoked even in homothermal animals. For example, the splicing pattern of mRNA for a cold shock protein, cold-inducible RNA-binding protein (CIRP or CIRBP), is changed in response to a marked drop in body temperature during hibernation of hamsters. In this review, we describe the current knowledge about mechanisms and functions of temperature-dependent alternative splicing in plants and animals. Then we discuss the physiological significance of hypothermia-induced alternative splicing of a cold shock protein gene in hibernating and non-hibernating animals.