Geographic origin shapes the adaptive divergences of Rotaria rotatoria (Rotifera, Bdelloidea) to thermal stress: Insights from ecology and transcriptomics

Global warming has raised concerns regarding the potential impact on aquatic biosafety and health. To illuminate the adaptive mechanisms of bdelloid rotifers in response to global warming, the ecological and transcriptomic characteristics of two strains (HX and ZJ) of Rotaria rotatoria were investigated at 25°C and 35°C. Our results showed an obvious genetic divergence between the two geographic populations. Thermal stress significantly reduced the average lifespan of R. rotatoria in both strains, but increased the offspring production in the ZJ strain. Furthermore, the expression levels of genes Hsp70 were significantly upregulated in the HX strain, while GSTo1 and Cu/Zn-SOD were on the contrary. In the ZJ strain, the expression levels of genes Hsp70, CAT2, and GSTo1 were upregulated under thermal stress. Conversely, a significant decrease in the expression level of the Mn-SOD gene was observed in the ZJ strain under thermal stress. Transcriptomic profiling analysis revealed a total of 105 and 5288 differentially expressed genes (DEGs) in the HX and ZJ strains under thermal stress, respectively. The PCA results showed clear differences in gene expression pattern between HX and ZJ strains under thermal stress. Interestingly, compared to the HX strain, numerous downregulated DEGs in the ZJ strain were enriched into pathways related to metabolism under thermal stress, suggesting that rotifers from the ZJ strain prioritize resource allocation to reproduction by suppressing costly metabolic processes. This finding is consistent with the life table results. This study provides new insights into the adaptive evolution of aquatic animals in the context of global climate change.

A genome-wide screening of the 70 kDa heat shock protein (HSP70) genes in the rotifer Brachionus plicatilis sensu stricto with a characterization of two heat-inducible HSP70 genes

The 70 kDa heat shock proteins (HSP70s) and the constitutive members of the HSP70 family (heat shock cognates; HSC70s) play essential roles in various biological processes. The number of hsp70/hsc70 in the database is rapidly increasing because of their importance and the automatic annotation of newly sequenced genomes. However, accumulating evidence indicates that neither hsp70 nor hsc70 forms a monophyletic gene family, raising the need to reconsider the annotation strategy based on the traditional concept of the inducible HSP70 and constitutive HSC70s. The main aim of this study is to establish a systematic scheme to annotate hsp70-like genes taking the latest phylogenetic insights into account. We cloned two hsp70s from the rotifer Brachionus plicatilis sensu stricto (s.s.), an emerging model in evolutionary genetics, and conducted a genome-wide screening of B. plicatilis s.s. hsp70s using the two sequences as queries. A total of 15 hsp70-like genes were found, and 7 of them encoded distant members of the HSP70 family, the function of which largely remains unknown. Eight canonical hsp70s were annotated according to a recently proposed nomenclature based on the molecular evolution: e.g., HSP70cA1/B1 for the cytosolic lineage, HSP70er1 for the endoplasmic reticulum lineage, and HSP70m1 for the mitochondrial lineage. The two cloned hsp70s, HSP70cB1 and HSP70cB2, ubiquitously increased their mRNA levels up to 7.5 times after heat treatment as assessed by semi-quantitative PCR, real-time PCR, and in situ hybridization. This systematic screening incorporating a reasonable update to the annotation strategy would provide a useful example for future HSP70 studies, especially those in non-traditional model organisms.