Environmental change and predator diversity drive alpha and beta diversity in freshwater macro and microorganisms

Silica aggravates pulmonary fibrosis through disrupting lung microbiota and amino acid metabolites

Silicosis, recognized as a severe global public health issue, is an irreversible pulmonary fibrosis caused by the long-term inhalation of silica particles. Given the intricate pathogenesis of silicosis, there is no effective intervention measure, which poses a severe threat to public health. Our previous study reported that dysbiosis of lung microbiota is associated with the development of pulmonary fibrosis, potentially involving the lipopolysaccharides/toll-like receptor 4 pathway. Similarly, the process of pulmonary fibrosis is accompanied by alterations in metabolic pathways. This study employed a combined approach of 16S rDNA sequencing and metabolomic analysis to investigate further the role of lung microbiota in silicosis delving deeper into the potential pathogenesis of silicosis. Silica exposure can lead to dysbiosis of the lung microbiota and the occurrence of pulmonary fibrosis, which was alleviated by a combination antibiotic intervention. Additionally, significant metabolic disturbances were found in silicosis, involving 85 differential metabolites among the three groups, which are mainly focused on amino acid metabolic pathways. The changed lung metabolites showed a substantial correlation with lung microbiota. The relative abundance of Pseudomonas negatively correlated with L-Aspartic acid, L-Glutamic acid, and L-Threonine levels. These results indicate that dysbiosis in pulmonary microbiota exacerbates silica-induced fibrosis through impacts on amino acid metabolism, providing new insights into the potential mechanisms and interventions of silicosis.

Brownification increases the abundance of microorganisms related to carbon and nitrogen cycling in shallow lakes

Brownification in aquatic ecosystems under global change has attracted attention. The composition and quantity of dissolved organic matter transported from various land use types to lakes differ significantly, causing varying ecological effects of lake brownification by region. Bacterial communities make a significant contribution to the material cycle of ecosystems and are sensitive to environmental changes. In this study, a series of mesocosm systems were used to simulate forest lakes and urban lakes with different degrees of brownification, and a high-throughput amplicon sequencing technique was used to explore the changes in the composition, structure, and function of bacterial communities in shallow lakes undergoing brownification. Principal coordinate analysis (PCoA) and Jensen‒Shannon distance typing analysis both indicated significant differences in bacterial communities between forest lakes and urban lakes. The α diversity of bacterial communities in urban lakes increased with the degree of brownification. However, whether forest lakes or urban lakes, brownification increased the abundance of carbon cycling-related bacterial phyla (Proteobacteria, Poribacteria, and Chloroflexi) and nitrogen cycling-related bacterial genera (Microbacteriaceae, Limnohabitans, Comamonadaceae, Bacillus, and Rhizobiales_Incertae_Sedis). Additionally, the carbon and nitrogen cycling functions of bacterial communities in forest lakes are dominant, while those in urban lakes are dominated by functions related to light. Our study has preliminarily revealed that lake brownification promotes the growth of carbon and nitrogen cycling microorganisms, providing a new paradigm for understanding the response of lake ecosystems in different catchment areas to environmental changes and the carbon and nitrogen cycling processes in shallow lake ecosystems.

Elucidating biogeochemical characterization of nitrogen in the vadose zone integrating geochemistry, microorganism, and numerical simulation

A thorough comprehension of nitrogen biogeochemical processes in the vadose zone is crucial for the effective prevention and remediation of soil-groundwater system contamination. Despite the growing research on this subject, the full scope of nitrogen biogeochemical characterization in different geological environments remains poorly understood. This study addresses this knowledge gap by integrating geochemical, microbiological and numerical simulation approaches to gain a deeper insight into nitrogen biogeochemistry in agriculture. Our findings indicate the biogeochemical behavior of nitrogen in the vadose zone is mediated by microorganisms, driven by hydraulics, influenced by geological conditions and environmental factors. Along the groundwater flow, NH4+-N was found to be heavily accumulated in the topsoil of 0-40 cm, while NO3--N was transported and driven by hydrodynamics from both vertical and horizontal directions. Microbial diversity, species composition and functional microorganisms were significantly influenced by soil depth, rather than geomorphological types. Oxidation-reduction potential (ORP), total organic carbon (TOC), soil moisture (MOI), bicarbonate (HCO3-), and ferrous (Fe2+) were identified as the principal environmental factors that regulate nitrogen metabolism and the dominant biochemical processes, encompassing nitrogen fixation, nitrification, and denitrification. Driven by hydrodynamics, NH4+-N, NO2--N and NO3--N tend to form distinct biochemical reaction zones in the vertical vadose zone. These areas are dynamic and subject to geomorphologies. It should be noted that NO3--N can migrate towards groundwater from the clayey sand in the Alluvial Plain, which presents a potential risk of groundwater contamination. The fissure structure of loess may serve as the major transport pathway for groundwater nitrogen contamination in the Loess Tableland. This finding highlights the importance of integrating microbiology, geochemistry and hydraulics to elucidate the biogeochemical processes of nitrogen in the vadose zone with a dynamic mindset.

Variations in the Odonata Assemblages: How Do the Dry Season and Water Bodies Influence Them?

Diverse abiotic and biotic factors drive the ecological variation of communities across spatial and temporal dimensions. Within the Amazonian landscape, various freshwater environments exhibit distinct physicochemical characteristics. Thus, our study delved into the fluctuations of Odonata assemblages amidst distinct water bodies within Amazonia, encompassing two distinct climatic seasons. Comparative analysis was conducted on Odonata species diversity and assemblage composition across a blackwater pond, a lake, and a stream, spanning the initiation and culmination of the dry season in the southwestern Amazon region in Peru. Our methodology involved capturing adult Odonata using entomological nets on three separate occasions between 11:00 and 14:00 h for each water body in May (beginning of the dry season) and October (end of the dry season) of 2018. We also evaluated the influence of temperature, precipitation, and percent cloud cover on the abundance and richness of adult Odonata. Species richness and composition differed among the three water bodies in both periods of the dry season. No effect of the dry season periods on species richness and abundance was observed. However, except in the oxbow lake, the more abundant species were substituted to the end of the dry season. Our study highlights the influence of water body types on Odonata species diversity and composition. The effects of the sampling period during the dry season may not be immediately apparent in conventional diversity metrics, such as species richness and abundance. Instead, its effects manifest predominantly in the relative abundance of the species that compose these assemblages.

Investigation of the mechanism of silica-induced pulmonary fibrosis: The role of lung microbiota dysbiosis and the LPS/TLR4 signaling pathway

The widespread manufacture of silica and its extensive use, and potential release of silica into the environment pose a serious human health hazard. Silicosis, a severe global public health issue, is caused by exposure to silica, leading to persistent inflammation and fibrosis of the lungs. The underlying pathogenic mechanisms of silicosis remain elusive. Lung microbiota dysbiosis is associated with the development of inflammation and fibrosis. However, limited information is currently available regarding the role of lung microbiota in silicosis. The study therefore is designed to conduct a comprehensive analysis of the role of lung microbiota dysbiosis and establish a basis for future investigations into the potential mechanisms underlying silicosis. Here, the pathological and biochemical parameters were used to systematically assessed the degree of inflammation and fibrosis following silica exposure and treatment with combined antibiotics. The underlying mechanisms were studied via integrative multi-omics analyses of the transcriptome and microbiome. Analysis of 16S ribosomal DNA revealed dysbiosis of the microbial community in silicosis, characterized by a predominance of gram-negative bacteria. Exposure to silica has been shown to trigger lung inflammation and fibrosis, leading to an increased concentration of lipopolysaccharides in the bronchoalveolar lavage fluid. Furthermore, Toll-like receptor 4 was identified as a key molecule in the lung microbiota dysbiosis associated with silica-induced lung fibrosis. All of these outcomes can be partially controlled through combined antibiotic administration. The study findings demonstrate that the dysbiosis of lung microbiota enhances silica-induced fibrosis associated with the lipopolysaccharides/Toll-like receptor 4 pathway and provided a promising target for therapeutic intervention of silicosis.

Influence of Ecological Multiparameters on Facets of β-Diversity of Freshwater Plankton Ciliates

Understanding the relative importance of the factors that drive global patterns of biodiversity is among the major topics of ecological and biogeographic research. In freshwater bodies, spatial, temporal, abiotic, and biotic factors are important structurers of these ecosystems and can trigger distinct responses according to the facet of biodiversity considered. The objective was to evaluate how different facets of β-diversity (taxonomic, functional, and phylogenetic) based on data from the planktonic ciliate community of a Neotropical floodplain, are influenced by temporal, spatial, abiotic, and biotic factors. The research was conducted in the upper Paraná River floodplain between the years 2010 and 2020 in different water bodies. All predictors showed significant importance on the facets of β-diversity, except the abiotic predictors on species composition data, for the taxonomic facet. The functional and phylogenetic facets were mostly influenced by abiotic, biotic, and spatial factors. For temporal predictors, results showed influence on taxonomic (structure and composition data) and functional (structure data) facets. Also, a fraction of shared explanation between the temporal and abiotic components was observed for the distinct facets. Significant declines in β-diversity in continental ecosystems have been evidenced, especially those with drastic implications for ecosystemic services. Therefore, the preservation of a high level of diversity in water bodies, also involving phylogenetic and functional facets, should be a priority in conservation plans and goals, to ensure the maintenance of important ecological processes involving ciliates.

Biotic homogenisation and differentiation as directional change in beta diversity: synthesising driver-response relationships to develop conceptual models across ecosystems

Biotic homogenisation is defined as decreasing dissimilarity among ecological assemblages sampled within a given spatial area over time. Biotic differentiation, in turn, is defined as increasing dissimilarity over time. Overall, changes in the spatial dissimilarities among assemblages (termed 'beta diversity') is an increasingly recognised feature of broader biodiversity change in the Anthropocene. Empirical evidence of biotic homogenisation and biotic differentiation remains scattered across different ecosystems. Most meta-analyses quantify the prevalence and direction of change in beta diversity, rather than attempting to identify underlying ecological drivers of such changes. By conceptualising the mechanisms that contribute to decreasing or increasing dissimilarity in the composition of ecological assemblages across space, environmental managers and conservation practitioners can make informed decisions about what interventions may be required to sustain biodiversity and can predict potential biodiversity outcomes of future disturbances. We systematically reviewed and synthesised published empirical evidence for ecological drivers of biotic homogenisation and differentiation across terrestrial, marine, and freshwater realms to derive conceptual models that explain changes in spatial beta diversity. We pursued five key themes in our review: (i) temporal environmental change; (ii) disturbance regime; (iii) connectivity alteration and species redistribution; (iv) habitat change; and (v) biotic and trophic interactions. Our first conceptual model highlights how biotic homogenisation and differentiation can occur as a function of changes in local (alpha) diversity or regional (gamma) diversity, independently of species invasions and losses due to changes in species occurrence among assemblages. Second, the direction and magnitude of change in beta diversity depends on the interaction between spatial variation (patchiness) and temporal variation (synchronicity) of disturbance events. Third, in the context of connectivity and species redistribution, divergent beta diversity outcomes occur as different species have different dispersal characteristics, and the magnitude of beta diversity change associated with species invasions also depends strongly on alpha and gamma diversity prior to species invasion. Fourth, beta diversity is positively linked with spatial environmental variability, such that biotic homogenisation and differentiation occur when environmental heterogeneity decreases or increases, respectively. Fifth, species interactions can influence beta diversity via habitat modification, disease, consumption (trophic dynamics), competition, and by altering ecosystem productivity. Our synthesis highlights the multitude of mechanisms that cause assemblages to be more or less spatially similar in composition (taxonomically, functionally, phylogenetically) through time. We consider that future studies should aim to enhance our collective understanding of ecological systems by clarifying the underlying mechanisms driving homogenisation or differentiation, rather than focusing only on reporting the prevalence and direction of change in beta diversity, per se.

Effects of warming on the structure of aquatic communities in tropical bromeliad microecosystems

Freshwaters are among the most vulnerable ecosystems to climate warming, with projected temperature increases over the coming decades leading to significant losses of aquatic biodiversity. Experimental studies that directly warm entire natural ecosystems in the tropics are needed, for understanding the disturbances on aquatic communities. Therefore, we conducted an experiment to test the impacts of predicted future warming on density, alpha diversity, and beta diversity of freshwater aquatic communities, inhabiting natural microecosystems-Neotropical tank bromeliads. Aquatic communities within the tanks bromeliads were experimentally exposed to warming, with temperatures ranging from 23.58 to 31.72°C. Linear regression analysis was used to test the impacts of warming. Next, distance-based redundancy analysis was performed to assess how warming might alter total beta diversity and its components. This experiment was conducted across a gradient of habitat size (bromeliad water volume) and availability of detrital basal resources. A combination of the highest detritus biomass and higher experimental temperatures resulted in the greatest density of flagellates. However, the density of flagellates declined in bromeliads with higher water volume and lower detritus biomass. Moreover, the combination of the highest water volume and high temperature reduced density of copepods. Finally, warming changed microfauna species composition, mostly through species substitution (β _repl component of total beta-diversity). These findings indicate that warming strongly structures freshwater communities by reducing or increasing densities of different aquatic communities groups. It also enhances beta-diversity, and many of these effects are modulated by habitat size or detrital resources.

Revealing the driving synergistic degradation mechanism of Rhodococcus sp. B2 on the bioremediation of pretilachlor-contaminated soil

The pretilachlor has been widely used worldwide and has contaminated the environment for many years. The environmental fate of pretilachlor and its residues removal from the contaminated environment have attracted great concern. Reportedly, pretilachlor could partly be transformed to HECDEPA by Rhodococcus sp. B2. However, the effects of pretilachlor on soil bacterial communities and its complete metabolic pathway remain unknown. Herein, we investigated the mechanism of driving synergistic degradation of pretilachlor by strain B2 in the soil. The results revealed that pretilachlor showed a negative effect on bacterial communities and caused significant variations in the community structure. Strain B2 showed the ability to remediate the pretilachlor-contaminated soils and network analysis revealed that it may drive the enrichment of potential pretilachlor-degrading bacteria from the soil. The soil pretilachlor degradation may be facilitated by the members of the keystone families Comamonadaceae, Caulobacteraceae, Rhodospirillaceae, Chitinophagaceae, and Sphingomonadaceae. Meanwhile, Sphingomonas sp. M6, a member of the Sphingomonadaceae family, has been isolated from the strain B2 inoculation sample soil. The co-culture, comprising strain M6 and B2, could synergistic degrade pretilachlor within 30 h, which is the highest degradation rate. Strain M6 could completely degrade the HECDEPA via CDEPA and DEA. In the soil, a comparable pretilachlor degradation pathway may exist. This study suggested that strain B2 had the potential to drive the remediation of pretilachlor-contaminated soils.

Lung microbiome and transcriptome reveal mechanisms underlying PM2.5 induced pulmonary fibrosis

Airborne fine particulate matter (PM_2.5) is considered to be a risk factor for lung fibrosis, and therefore, it has attracted public attention due to its various physicochemical features and its adverse effects on health. However, little remains to be known regarding the mechanism of PM_2.5-induced pulmonary fibrosis. The lung microbiota may be a potential factor involved in the adverse outcomes of pulmonary fibrosis. Meanwhile, miRNAs are thought to be key regulators that participate in the complex interplay between the host and the microbiota. Hence, to investigate the potential mechanisms of pulmonary fibrosis, and to explore the impact of PM_2.5-induced alterations in miRNAs and the lung microbiota and possible interaction patterns in mice models, we took advantage of 16S rDNA gene sequencing, miRNAs sequencing (miRNAs-Seq), and mining of public databases profiling. The results of 16S rDNA analysis showed that PM_2.5 interfered with the microbial community composition, resulting in Proteobacteria becoming an additional dominant phylum. In addition, differentially expressed miRNAs were enriched in HIF-1 signaling, the IL-17 signaling, as well as Th17 cell differentiation pathways, which are closely related to microbial functional pathways. Significantly, a target miRNA, miR-149-5p, may be a key factor triggering the MAPK signal pathway related to pulmonary fibrosis and disturbing the homeostasis of lung bacterial flora. These results indicate that PM_2.5 may lead to interaction between lung microbiota dysbiosis and an imbalance of miRNA levels to form a vicious cycle that promotes lung fibrogenesis. The current study provides new insights into the progression of pulmonary fibrosis.

Alterations of the rhizosphere soil microbial community composition and metabolite profiles of Zea mays by polyethylene-particles of different molecular weights

Polyethylene film is the most widely used plastic film in agricultural production activities, and its depolymerization products are mainly polyethylene-particles (PE-particles) of different molecular weights. However, the impact of the molecular weights of the PE-particles on soil-crop microenvironment has not been elucidated. In this study, a potted microcosmic simulation system was used to evaluate the impact of low, medium and high molecular weight PE-particles on soil metabolism, microbial community structure, and crop growth. There were obvious differences in the shape and surface microstructure of PE-particles with different molecular weights. Soil sucrase and peroxidase had significant responses to PE-particles of different molecular weights. In the rhizosphere, the number of microorganisms and the microbial alpha diversity index increased with increasing PE-particles molecular weight. Rhizobacter, Nitrospira, and Sphingomonas were the dominant microorganisms induced by PE-particles to regulate the metabolism of elements. Carbohydrate and amino acid contents in rhizosphere soils were the key factors affecting the species abundance of Lysobacter, Polyclovorans, Rhizobacter, and Sphingomonas. In plants, PE-particles treatment reduced the plant biomass and photosynthetic rate and disrupted normal mineral nutrient metabolism. Different molecular weight PE-particles may therefore have adverse effects on the soil-plant system.

Warming and top predator loss drive direct and indirect effects on multiple trophic groups within and across ecosystems

The interspecific interactions within and between adjacent ecosystems strongly depend on the changes in their abiotic and biotic components. However, little is known about how climate change and biodiversity loss in a specific ecosystem can impact the multiple trophic interactions of different biological groups within and across ecosystems. We used natural microecosystems (tank-bromeliads) as a model system to investigate the main and interactive effects of aquatic warming and aquatic top predator loss (i.e. trophic downgrading) on trophic relationships in three integrated food web compartments: (a) aquatic micro-organisms, (b) aquatic macro-organisms and (c) terrestrial predators (i.e. via cross-ecosystem effects). The aquatic top predator loss substantially impacted the three food web compartments. In the aquatic macrofauna compartment, trophic downgrading increased the filter feeder richness and abundance directly and indirectly via an increase in detritivore richness, likely through a facilitative interaction. For the microbiota compartment, aquatic top predator loss had a negative effect on algae richness, probably via decreasing the input of nutrients from predator biological activities. Furthermore, the more active terrestrial predators responded more to aquatic top predator loss, via an increase in some components of aquatic macrofauna, than more stationary terrestrial predators. The aquatic trophic downgrading indirectly altered the richness and abundance of cursorial terrestrial predators, but these effects had different direction according to the aquatic functional group, filter feeder or other detritivores. The web-building predators were indirectly affected by aquatic trophic downgrading due to increased filter feeder richness. Aquatic warming did not affect the aquatic micro- or macro-organisms but did positively affect the abundance of web-building terrestrial predators. These results allow us to raise a predictive framework of how different anthropogenic changes predicted for the next decades, such as aquatic warming and top predator loss, could differentially affect multiple biological groups through interactions within and across ecosystems.

Trophic downgrading decreases species asynchrony and community stability regardless of climate warming

Theory and some evidence suggest that biodiversity promotes stability. However, evidence of how trophic interactions and environmental changes modulate this relationship in multitrophic communities is lacking. Given the current scenario of biodiversity loss and climate changes, where top predators are disproportionately more affected, filling these knowledge gaps is crucial. We simulated climate warming and top predator loss in natural microcosms to investigate their direct and indirect effects on temporal stability of microbial communities and the role of underlying stabilising mechanisms. Community stability was insensitive to warming, but indirectly decreased due to top predator loss via increased mesopredator abundance and consequent reduction of species asynchrony and species stability. The magnitude of destabilising effects differed among trophic levels, being disproportionally higher at lower trophic levels (e.g. producers). Our study unravels major patterns and causal mechanisms by which trophic downgrading destabilises large food webs, regardless of climate warming scenarios.

Cascading effects of moth outbreaks on subarctic soil food webs

The increasing severity and frequency of natural disturbances requires a better understanding of their effects on all compartments of biodiversity. In Northern Fennoscandia, recent large-scale moth outbreaks have led to an abrupt change in plant communities from birch forests dominated by dwarf shrubs to grass-dominated systems. However, the indirect effects on the belowground compartment remained unclear. Here, we combined eDNA surveys of multiple trophic groups with network analyses to demonstrate that moth defoliation has far-reaching consequences on soil food webs. Following this disturbance, diversity and relative abundance of certain trophic groups declined (e.g., ectomycorrhizal fungi), while many others expanded (e.g., bacterivores and omnivores) making soil food webs more diverse and structurally different. Overall, the direct and indirect consequences of moth outbreaks increased belowground diversity at different trophic levels. Our results highlight that a holistic view of ecosystems improves our understanding of cascading effects of major disturbances on soil food webs.

Warming and leaf litter functional diversity, not litter quality, drive decomposition in a freshwater ecosystem

Environment, litter composition and decomposer community are known to be the main drivers of litter decomposition in aquatic ecosystems. However, it remains unclear whether litter quality or functional diversity prevails under warming conditions. Using tank bromeliad ecosystems, we evaluated the combined effects of warming, litter quality and litter functional diversity on the decomposition process. We also assessed the contribution of macroinvertebrates and microorganisms in explaining litter decomposition patterns using litter bags made with different mesh sizes. Our results showed that litter decomposition was driven by litter functional diversity and was increasingly higher under warming, in both mesh sizes. Decomposition was explained by increasing litter dissimilarities in C and N. Our results highlight the importance of considering different aspects of litter characteristics (e.g., quality and functional diversity) in order to predict the decomposition process in freshwater ecosystems. Considering the joint effect of warming and litter traits aspects allow a more refined understanding of the underlying mechanisms of climate change and biodiversity shifts effects on ecosystem functioning.

The rise of a marine generalist predator and the fall of beta diversity

Determining the importance of physical and biological drivers in shaping biodiversity in diverse ecosystems remains a global challenge. Advancements have been made towards this end in large marine ecosystems with several studies suggesting environmental forcing as the primary driver. However, both empirical and theoretical studies point to additional drivers of changes in diversity involving trophic interactions and, in particular, predation. Moreover, a more integrated but less common approach to the assessment of biodiversity changes involves analyses of spatial β diversity, whereas most studies to date assess only changes in species richness (α diversity). Recent research has established that when cod, a dominant generalist predator, was overfished and collapsed in a northwest Atlantic food web, spatial β diversity increased; that is, the spatial structure of the fish assemblage became increasingly heterogeneous. If cod were to recover, would this situation be reversible, given the inherent complexity and non-linear dynamics that typify such systems? A dramatic increase of cod in an ecologically similar large marine ecosystem may provide an answer. Here we show that spatial β diversity of fish assemblages in the Barents Sea decreased with increasing cod abundance, while decadal scale changes in temperature did not play a significant role. These findings indicate a reversibility of the fish assemblage structure in response to changing levels of an apex predator and highlight the frequently overlooked importance of trophic interactions in determining large-scale biodiversity patterns. As increased cod abundance was largely driven by changes in fisheries management, our study also shows that management policies and practices, particularly those involving apex predators, can have a strong effect in shaping spatial diversity patterns, and one should not restrict the focus to effects of climate change alone.

Comparative Transcriptome and Proteome Analysis of Heat Acclimation in Predatory Mite Neoseiulus barkeri

In our previous study, we reported a high temperature adapted strain (HTAS) of the predatory mite Neoseiulus barkeri was artificially selected via a long-term heat acclimation (35°C) and frequent heat hardenings. To understand the molecular basis of heat acclimation, 'omics' analyses were performed to compare the differences between HTAS female adults to conventional strain (CS) at transcriptional and translational levels. We obtained a total of 5,374 differentially expressed genes and 500 differentially expressed proteins. Among them, 119 transcripts had concurrent transcription and translation profiles. It's conserved that some processes, such as high expression of heat shock protein (HSP) genes, involved in heat tolerance of transcriptome analyses, while many protective enzymes including glutathione S-transferase, superoxide dismutase, peroxidase, and cytochrome P450 displayed down-regulated expression. KEGG analysis mapped 4,979 and 348 differentially expressed genes and proteins, to 299 and 253 pathways, respectively. The mitogen-activated protein kinases (MAPK) signaling pathway may provide new insights for the investigation of the molecular mechanisms of heat tolerance. Correlation enriched pathways indicated that there were four pathways associated with heat acclimation involving in energy metabolism and immunity. In addition, the expression patterns of ten randomly selected genes including HSP were consistent with the transcriptome results obtained through quantitative real-time PCR. Comparisons between transcriptome and proteome results indicated the upregulation of HSPs and genes participated in ATP production, immunity and energy metabolism process. A majority of antioxidant-related genes and detoxication-related genes were down-regulated suggesting a fitness cost of heat acclimation. Our results demonstrated that heat tolerance during a long-time acclimation of N. barkeri is a fairly complicated process of physiological regulations. These findings also contribute to a better understanding of the mechanisms of thermal responses of phytoseiid mites which could provide useful information for biological control through natural enemies.