Temperature-dependent life history and transcriptomic responses in heat-tolerant versus heat-sensitive Brachionus rotifers

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.

Gene and Allele-Specific Expression Underlying the Electric Signal Divergence in African Weakly Electric Fish

In the African weakly electric fish genus Campylomormyrus, electric organ discharge signals are strikingly different in shape and duration among closely related species, contribute to prezygotic isolation, and may have triggered an adaptive radiation. We performed mRNA sequencing on electric organs and skeletal muscles (from which the electric organs derive) from 3 species with short (0.4 ms), medium (5 ms), and long (40 ms) electric organ discharges and 2 different cross-species hybrids. We identified 1,444 upregulated genes in electric organ shared by all 5 species/hybrid cohorts, rendering them candidate genes for electric organ-specific properties in Campylomormyrus. We further identified several candidate genes, including KCNJ2 and KLF5, and their upregulation may contribute to increased electric organ discharge duration. Hybrids between a short (Campylomormyrus compressirostris) and a long (Campylomormyrus rhynchophorus) discharging species exhibit electric organ discharges of intermediate duration and showed imbalanced expression of KCNJ2 alleles, pointing toward a cis-regulatory difference at this locus, relative to electric organ discharge duration. KLF5 is a transcription factor potentially balancing potassium channel gene expression, a crucial process for the formation of an electric organ discharge. Unraveling the genetic basis of the species-specific modulation of the electric organ discharge in Campylomormyrus is crucial for understanding the adaptive radiation of this emerging model taxon of ecological (perhaps even sympatric) speciation.

Temporal Distribution Patterns of Cryptic Brachionus calyciflorus (Rotifera) Species in Relation to Biogeographical Gradient Associated with Latitude

Sympatric distribution and temporal overlap of cryptic zooplankton species pose a challenge to the framework of the niche differentiation theory and the mechanisms allowing competitor coexistence. We applied the methods of phylogenetic analysis, DNA taxonomy, and statistical analysis to study the temporal distribution patterns of the cryptic B. calyciflorus species, an excellent model, in three lakes, and to explore the putative mechanisms for their seasonal succession and temporal overlap. The results showed that in the warm-temperate Lake Yunlong, B. fernandoi and B. calyciflorus s.s. underwent a seasonal succession, which was largely attributed to their differential adaptation to water temperature. In the subtropical Lake Jinghu, B. fernandoi, B. calyciflorus s.s., and B. dorcas exhibited both seasonal succession and temporal overlap. Seasonal successions were largely attributed to their differential adaptation to temperature, and temporal overlap resulted from their differential responses to algal food concentration. In the tropical Lake Jinniu, B. calyciflorus s.s. persisted throughout the year and overlapped with B. dorcas for 5 months. The temporal overlap resulted from their differential responses to copepod predation. These results indicated that the temporal distribution pattern of the cryptic B. calyciforus species and the mechanism that allows competitor coexistence vary with different climate zones.

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.

Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

Differential transcriptomic responses to heat stress in surface and subterranean diving beetles

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CT_max 38 °C cf. 42–46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.

The genome of the freshwater monogonont rotifer Brachionus rubens: Identification of phase I, II, and III detoxification genes

Monogonont rotifers are common species in aquatic environments and make model species for ecotoxicology studies. Whole genomes of several species of the genus Brachionus have been assembled, but no information on the freshwater rotifer Brachionus rubens has been reported. In this study, the whole-genome sequence of B. rubens was successfully assembled using NextDenovo. The total length of the genome was 132.7 Mb (N50 = 2.51 Mb), including 122 contigs. The GC contents accounted for 29.96% of the genome. Aquatic organisms are always exposed to various external stresses, and a comprehensive genomic analysis is needed to better understand the adverse effects on organisms. This paper focuses on the ecotoxicological aspect and conducted genome analysis of representative gene families involved in detoxification mechanisms against environmental stressors. Specifically, we identified cytochrome P450 genes (CYPs) of phase I, glutathione S-transferase genes (GSTs) of phase II, and ATP-binding cassette transporter genes (ABCs) of phase III in the genome of B. rubens. Gene duplications were found in CYP, GST, and ABC genes, as is the case for other Brachionus rotifers. Our results suggest that these detoxification-related gene families have evolved in a species-specific and/or lineage-specific manner. This paper improves our understanding of how the freshwater Brachionus rotifers respond to environmental stressors in a molecular ecotoxicology context.