Transcriptomic responses to heat stress in gill and liver of endangered Brachymystax lenok tsinlingensis

Application of Environmental DNA for Assessing the Distribution and Biomass of Brachymystax lenok Tsinlingensis in the Zhouzhi Heihe River

In recent years, the population of Brachymystax lenok tsinlingensis has experienced a rapid decline due to habitat degradation caused by overfishing and dam construction. In the current study, water samples were collected from four tanks containing B. lenok tsinlingensis individuals with wet weights of 5.8 g, 12.3 g, 34.7 g, and 81.7 g, respectively, and the cycle threshold (Ct) values of B. lenok tsinlingensis in each tank were measured using quantitative PCR. Additionally, water samples were collected from six sites along the Zhouzhi Heihe River (Yu Dongquan, Ban Fangzi, Ba Mugou, two sites along Wang Jia, and Qingshui Rivers), and the Ct values of B. lenok tsinlingensis at these sites were quantified. The results show that the relationship between the Ct values measured by the quantitative test and the fish biomass in the corresponding tanks followed an exponential function, with a regression coefficient of 0.9987. The biomass for the six sites was determined using the established relationship, and the results indicate that the biomass of B. lenok tsinlingensis was highest in Yu Dongquan, followed by Ban Fangzi, Ba Mugou, Wang Jia River 2, and Qingshui, with the lowest biomass observed in Wang Jia River 1. These findings are consistent with the species' adaptive survival distribution across the six sites. This study presents a standardized approach for using an eDNA-based quantitative PCR assay to monitor B. lenok tsinlingensis. It provides an advanced technical method applicable to the protection and monitoring of other aquatic species in future studies.

Analysis of Differential Alternative Splicing in Largemouth Bass After High Temperature Exposure

Temperature is one of the critical factors affecting the physiological functions of fish. With ongoing global warming, changes in water temperature have a profound impact on fish species. Alternative splicing, being a significant mechanism for gene expression regulation, facilitates fish to adapt and thrive in dynamic and varied aquatic environments. Our study used transcriptome sequencing to analyze alternative splicing in largemouth bass gills at 34 °C for 24 h. The findings indicated an increase in both alternative splicing events and alternative splicing genes after high temperature treatment. Specifically, the comparative analysis revealed a total of 674 differential alternative splicing events and 517 differential alternative splicing genes. Enrichment analysis of differential alternative splicing genes revealed significant associations with various gene ontology (GO) terms and KEGG pathways, particularly in immune-related pathways like necroptosis, apoptosis, and the C-type lectin receptor signaling pathway. These results emphasize that some RNA splicing-related genes are involved in the response of largemouth bass to high temperatures.

High-Temperature-Induced Differential Expression of miRNA Mediates Liver Inflammatory Response in Tsinling Lenok Trout (Brachymystax lenok tsinlingensis)

High-temperature stress poses a significant environmental challenge for aquatic organisms, including tsinling lenok trout (Brachymystax lenok tsinlingensis). This study aimed to investigate the role of microRNAs (miRNAs) in inducing liver inflammation in tsinling lenok trout under high-temperature stress. Tsinling lenok trout were exposed to high-temperature conditions (24 °C) for 8 h, and liver samples were collected for analysis. Through small RNA sequencing, we identified differentially expressed miRNAs in the liver of high-temperature-stressed tsinling lenok trout compared to the control group (maintained at 16 °C). Several miRNAs, including novel-m0105-5p and miR-8159-x, showed significant changes in expression levels. Additionally, we conducted bioinformatics analysis to explore the potential target genes of these differentially expressed miRNAs. Our findings revealed that these miRNA target genes are involved in inflammatory response pathways, such as NFKB1 and MAP3K5. The downregulation of novel-m0105-5p and miR-8159-x in the liver of high-temperature-stressed tsinling lenok trout suggests their role in regulating liver inflammatory responses. To validate this, we performed a dual-luciferase reporter assay to confirm the regulatory relationship between miRNAs and target genes. Our results demonstrated that novel-m0105-5p and miR-8159-x enhance the inflammatory response of hepatocytes by promoting the expression of NFKB1 and MAP3K5, respectively. In conclusion, our study provides evidence that high-temperature stress induces liver inflammation in tsinling lenok trout through dysregulation of miRNAs. Understanding the molecular mechanisms underlying the inflammatory response in tsinling lenok trout under high-temperature stress is crucial for developing strategies to mitigate the negative impacts of environmental stressors on fish health and aquaculture production.

Tissue-specific transcriptional response of post-larval clownfish to ocean warming

Anthropogenically driven climate change is predicted to increase average sea surface temperatures, as well as the frequency and intensity of marine heatwaves in the future. This increasing temperature is predicted to have a range of negative physiological impacts on multiple life-stages of coral reef fish. Nevertheless, studies of early-life stages remain limited, and tissue-specific transcriptomic studies of post-larval coral reef fish are yet to be conducted. Here, in an aquaria-based study we investigate the tissue-specific (brain, liver, muscle, and digestive tract) transcriptomic response of post-larval (20 dph) Amphiprion ocellaris to temperatures associated with future climate change (+3 °C). Additionally, we utilized metatranscriptomic sequencing to investigate how the microbiome of the digestive tract changes at +3 °C. Our results show that the transcriptional response to elevated temperatures is highly tissue-specific, as the number of differentially expressed genes (DEGs) and gene functions varied amongst the brain (102), liver (1785), digestive tract (380), and muscle (447). All tissues displayed DEGs associated with thermal stress, as 23 heat-shock protein genes were upregulated in all tissues. Our results indicate that post-larval clownfish may experience liver fibrosis-like symptoms at +3 °C as genes associated with extracellular matrix structure, oxidative stress, inflammation, glucose transport, and metabolism were all upregulated. We also observe a shift in the digestive tract microbiome community structure, as Vibrio sp. replace Escherichia coli as the dominant bacteria. This shift is coupled with the dysregulation of various genes involved in immune response in the digestive tract. Overall, this study highlights post-larval clownfish will display tissue-specific transcriptomic responses to future increases in temperature, with many potentially harmful pathways activated at +3 °C.

Redox state and metabolic responses to severe heat stress in lenok Brachymystax lenok (Salmonidae)

In order to provide new insights into the physiological responses of lenok (Brachymystax lenok: Salmonidae) to acute and severe heat stress (25°C, 48 h), dynamic changes in redox state and metabolic responses are studied combined biochemical index and non-targeted metabolome. Nicotinamide adenine dinucleotide (NAD⁺) consumption causes significant increases in ratio of reduced NADH to NAD⁺ and ratio of reduced nicotinamide adenine dinucleotide phosphate (NADPH) to NADP⁺, which induced the redox imbalance in heat stressed lenok. Lowered reduced glutathione/oxidized glutathione (GSH/GSSG) ratios suggested that more oxidized conditions occurred in heat-stressed lenok, leading to membrane lipid oxidation. The first few hours of heat stress promoted the activity of enzymes involved in anaerobic glycolysis (hexokinase, pyruvate kinase, lactic dehydrogenase) and glutamicpyruvic transaminase and glutamic oxaloacetic transaminase, which might lead to consumption of many carbohydrates and amino acid catabolism. These enzyme activities decreased with time in a possible compensatory strategy to manage anabolic and catabolic metabolism, maintaining the redox homeostasis. After 48 h of recovery, NAD⁺, carbohydrate levels and enzyme activities had returned to control levels, whereas many amino acids were consumed for repair and new synthesis. GSH remained at levels lower than controls, and the more oxidized conditions had not recovered, aggravating oxidative damage. Glutamic acid, glutamine, lysine and arginine may play important roles in survival of heat-stressed lenok.

Clownfish larvae exhibit faster growth, higher metabolic rates and altered gene expression under future ocean warming

Increasing ocean temperatures have been demonstrated to have a range of negative impacts on coral reef fishes. However, despite a wealth of studies of juvenile/adult reef fish, studies of how early developmental stages respond to ocean warming are limited. As overall population persistence is influenced by the development of early life stages, detailed studies of larval responses to ocean warming are essential. Here, in an aquaria-based study we investigate how temperatures associated with future warming and present-day marine heatwaves (+3 °C) impact the growth, metabolic rate, and transcriptome of 6 discrete developmental stages of clownfish larvae (Amphiprion ocellaris). A total of 6 clutches of larvae were assessed, with 897 larvae imaged, 262 larvae undergoing metabolic testing and 108 larvae subject to transcriptome sequencing. Our results show that larvae reared at +3 °C grow and develop significantly faster and exhibit higher metabolic rates than those in control conditions. Finally, we highlight the molecular mechanisms underpinning the response of larvae from different developmental stages to higher temperatures, with genes associated with metabolism, neurotransmission, heat stress and epigenetic reprogramming differentially expressed at +3 °C. Overall, these results indicate that clownfish development could be altered under future warming, with developmental rate, metabolic rate, and gene expression all affected. Such changes may lead to altered larval dispersal, changes in settlement time and increased energetic costs.

Combined effects of high temperature and pesticide mixture exposure on free-swimming behaviors and hepatic cytochrome P450 1A expression in goldfish, Carassius auratus

The synergy between multiple compounds and other stressors, including heat, creates volatility and greater unpredictability than standard single-chemical toxicity testing, especially in the case of pesticides and metabolites which might contain several noxious ingredients resulting in adverse ecological effects. To address this, the aim of this study was to examine the dose- and time-dependent effects of low- and high-dose pesticide mixture (metalachlor, linuron, isoproturon, tebucanazole, aclonifen, atrazine, pendimethalin, azinphos-methyl) and heat stress co-exposure (22°C control/32°C treatment for 4-week) on free-swimming behaviors and cumulative actionless time (CAT) of goldfish. Behavioral analysis showed a dose- and time-dependent decrease in distance swam, as well as a subsequent increase in CAT. Vertical and horizontal spatial behavioral use were affected under heat and pesticides co-exposure conditions. In 3- and 4-week(s) exposure groups, horizontal spatial behavioral use demonstrated elevated time spent in the lower third of the aquarium. Similarly, during 3- and 4-week(s) exposure (32°C control and 32°C high doses) vertical spatial behavioral use was found to increase time spent in the outermost edges of the aquarium. In all treatment groups, the final condition factor (KM) showed significant attenuation when compared to the initial KM. However, there was an unclear relationship between heat/pesticide co-exposure and growth most notably in 32°C high-dose groups. In addition, the expression of hepatic cytochrome P450 1A mRNA was significantly higher in pesticide-exposed groups. Taken together, data demonstrated that co-exposure with low- or high-dose pesticide mixture and heat stress significantly impacted natural swimming patterns, which over time might result in the broader population and ecological effects.

Sardines at a junction: Seascape genomics reveals ecological and oceanographic drivers of variation in the NW Mediterranean Sea

By evaluating genetic variation across the entire genome, one can address existing questions in a novel way while raising new ones. The latter includes how different local environments influence adaptive and neutral genomic variation within and among populations, providing insights into local adaptation of natural populations and their responses to global change. Here, under a seascape genomic approach, ddRAD data of 4609 single nucleotide polymorphisms (SNPs) from 398 sardines (Sardina pilchardus) collected in 11 Mediterranean and one Atlantic site were generated. These were used along with oceanographic and ecological information to detect signals of adaptive divergence with gene flow across environmental gradients. The studied sardines constitute two clusters (F_ST  = 0.07), a pattern attributed to outlier loci, highlighting putative local adaptation. The trend in the number of days with sea surface temperature above 19°C, a critical threshold for successful sardine spawning, was crucial at all levels of population structuring with implications on the species' key biological processes. Outliers link candidate SNPs to the region's environmental heterogeneity. Our findings provide evidence for a dynamic equilibrium in which population structure is maintained by physical and ecological factors under the opposing influences of migration and selection. This dynamic in a natural system warrants continuous monitoring under a seascape genomic approach that might benefit from a temporal and more detailed spatial dimension. Our results may contribute to complementary studies aimed at providing deeper insights into the mechanistic processes underlying population structuring. Those are key to understanding and predicting future changes and responses of this highly exploited species in the face of climate change.

Thermal Stress Induces Metabolic Responses in Juvenile Qingtian Paddy Field Carp Cyprinus carpio var qingtianensis

Extreme fluctuations in water temperature lead to significant economic losses for the aquaculture industry. Cyprinus carpio var qingtianensis (locally called Qingtian paddy field carp), is a local variety commonly found in Zhejiang province, China. Unlike traditional aquaculture environments, the water temperature range between day and night in the rice field environment is much larger, and the high temperature in summer may exceed the growth threshold of fish because there is no manual intervention; therefore, the study of how the Qingtian paddy field carp (PF carp) adapts to high-temperature conditions can shed light how the species adapt to the rice field environment. To investigate the molecular mechanisms of this fish under thermal stress, the liver metabolomics of Qiangtian paddy field carp (PF carp) were analyzed. In this study, metabolomics was used to examine the metabolic reaction of PF carp (102 days old, 104.69 ± 3.08 g in weight, 14.65 ± 0.46 cm in length) at water temperatures of 28 °C (control group, CG), 34 °C (experimental group (EG) 34), and 38 °C (EG38). The results show that 175 expression profile metabolites (DEMs), including 115 upregulated and 60 downregulated metabolites, were found in the CG vs. EG34. A total of 354 DEMs were inspected in CG vs. EG38, with 85 metabolites downregulated and 269 metabolites upregulated. According to the pathway enrichment study, various pathways were altered by thermal stress, including those of lipid, amino-acid, and carbohydrate metabolism. Our study presents a potential metabolic profile for PF carp under thermal stress. It also demonstrates how the host responds to thermal stress on a metabolic and molecular level.

Taxifolin Modulates Transcriptomic Response to Heat Stress in Rainbow Trout, Oncorhynchus mykiss

Taxifolin is a natural flavonoid known for its antioxidant, anti-inflammatory, and antiproliferative effects on animals. In this work, we have studied the effect of this compound on rainbow trout, Oncorhynchus mykiss, a major object of aquaculture, under slowly increasing ambient temperature and Gyrodactylus flatworm infection. Transcriptomic profiling of liver samples performed by using the Illumina HiSeq 2500 sequencing platform shows that a combined taxifolin/heat treatment, unlike heat treatment alone, downregulates the production of isopentenyl diphosphate, likely affecting the production of cholesterol and other sterols. Taxifolin treatment also modulates multiple apoptosis regulators and affects the expression of HSPs in response to increasing temperature. On the other hand, the expression of antioxidant enzymes in response to heat is not significantly affected by taxifolin. As for the Gyrodactylus infection, the parasite load is not affected by taxifolin treatment, although it was lower in the high-temperature group. Parasite load also did not induce a statistically significant transcriptomic response within the no heat/no taxifolin group.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6–26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.