Advances in Bioinspired Artificial Systems Enabling Biomarker-Driven Therapy.

Bioinspired molecular engineering strategies have emerged as powerful tools that significantly enhance the development of novel therapeutics, improving efficacy, specificity, and safety in disease treatment. Recent advancements have focused on identifying and utilizing disease-associated biomarkers to optimize drug activity and address challenges inherent in traditional therapeutics, such as frequent drug administrations, poor patient adherence, and increased risk of adverse effects. In this review, we provide a comprehensive overview of the latest developments in bioinspired artificial systems (BAS) that use molecular engineering to tailor therapeutic responses to drugs in the presence of disease-specific biomarkers. We examine the transition from open-loop systems, which rely on external cues, to closed-loop feedback systems capable of autonomous self-regulation in response to disease-associated biomarkers. We detail various BAS modalities designed to achieve biomarker-driven therapy, including activatable prodrug molecules, smart drug delivery platforms, autonomous artificial cells, and synthetic receptor-based cell therapies, elucidating their operational principles and practical in vivo applications. Finally, we discuss the current challenges and future perspectives in the advancement of BAS-enabled technology and envision that ongoing advancements toward more programmable and customizable BAS-based therapeutics will significantly enhance precision medicine.

Use of NAD tagSeq II to identify growth phase-dependent alterations in E. coli RNA NAD+ capping.

Recent findings regarding nicotinamide adenine dinucleotide (NAD+)-capped RNAs (NAD-RNAs) indicate that prokaryotes and eukaryotes employ noncanonical RNA capping to regulate gene expression. Two methods for transcriptome-wide analysis of NAD-RNAs, NAD captureSeq and NAD tagSeq, are based on copper-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry to label NAD-RNAs. However, copper ions can fragment/degrade RNA, interfering with the analyses. Here we report development of NAD tagSeq II, which uses copper-free, strain-promoted azide-alkyne cycloaddition (SPAAC) for labeling NAD-RNAs, followed by identification of tagged RNA by single-molecule direct RNA sequencing. We used this method to compare NAD-RNA and total transcript profiles of Escherichia coli cells in the exponential and stationary phases. We identified hundreds of NAD-RNA species in E. coli and revealed genome-wide alterations of NAD-RNA profiles in the different growth phases. Although no or few NAD-RNAs were detected from some of the most highly expressed genes, the transcripts of some genes were found to be primarily NAD-RNAs. Our study suggests that NAD-RNAs play roles in linking nutrient cues with gene regulation in E. coli.

SPAAC-NAD-seq, a sensitive and accurate method to profile NAD+-capped transcripts.

Nicotinamide adenine diphosphate (NAD+) is a novel messenger RNA 5' cap in Escherichia coli, yeast, mammals, and Arabidopsis Transcriptome-wide identification of NAD+-capped RNAs (NAD-RNAs) was accomplished through NAD captureSeq, which combines chemoenzymatic RNA enrichment with high-throughput sequencing. NAD-RNAs are enzymatically converted to alkyne-RNAs that are then biotinylated using a copper-catalyzed azide-alkyne cycloaddition (CuAAC) reaction. Originally applied to E. coli RNA, which lacks the m7G cap, NAD captureSeq was then applied to eukaryotes without extensive verification of its specificity for NAD-RNAs vs. m7G-capped RNAs (m7G-RNAs). In addition, the Cu2+ ion in the CuAAC reaction causes RNA fragmentation, leading to greatly reduced yield and loss of full-length sequence information. We developed an NAD-RNA capture scheme utilizing the copper-free, strain-promoted azide-alkyne cycloaddition reaction (SPAAC). We examined the specificity of CuAAC and SPAAC reactions toward NAD-RNAs and m7G-RNAs and found that both prefer the former, but also act on the latter. We demonstrated that SPAAC-NAD sequencing (SPAAC-NAD-seq), when combined with immunodepletion of m7G-RNAs, enables NAD-RNA identification with accuracy and sensitivity, leading to the discovery of new NAD-RNA profiles in Arabidopsis Furthermore, SPAAC-NAD-seq retained full-length sequence information. Therefore, SPAAC-NAD-seq would enable specific and efficient discovery of NAD-RNAs in prokaryotes and, when combined with m7G-RNA depletion, in eukaryotes.

Transcriptome analysis of the Nile tilapia (Oreochromis niloticus) reveals altered expression of immune genes by cadmium.

Nile tilapia (Oreochromis niloticus) is a worldwide farmed fish and has been widely used for the study on comparative immunology in teleosts. It is well known that cadmium (Cd) can cause a variety of adverse effects in fish. However, data on the effects of Cd in fish liver and the defensive mechanisms of these effects using transcriptome approach are relatively scarce to date. In this study, by using an RNA sequencing approach, the gene expression profiling was performed in livers of tilapia exposed to 0 (control), 50, 100, and 200 µg/L of Cd for 2 months. The results showed that exposure to 50 µg/L Cd altered the expressions of 911 genes, while exposure to 100 and 200 µg/L Cd resulted in 4318 and 3737 differentially expressed genes compared to the control. Weighted correlation network analysis (WGCNA) and gene ontology (GO) analysis identified a 14-gene network linked to the immune system development. Further, in a fuzzy analysis, the GO term immune system development was enriched in cluster 3, and gene expression decreased with increasing Cd levels in a concentration-dependent manner. The qPCR and RNA-seq results identified 4 genes, i.e., dnmt3bb.1, sf3b1, SMARCAL1, and zap70, as convenient potential biological indicators for detecting waterborne Cd. The present results help systematically understand the effects of Cd on the hepatic transcriptome in tilapia.

Using blood transcriptome analysis for Alzheimer's disease diagnosis and patient stratification.

INTRODUCTION: Blood protein biomarkers demonstrate potential for Alzheimer's disease (AD) diagnosis. Limited studies examine the molecular changes in AD blood cells. METHODS: Bulk RNA-sequencing of blood cells was performed on AD patients of Chinese descent (n = 214 and 26 in the discovery and validation cohorts, respectively) with normal controls (n = 208 and 38 in the discovery and validation cohorts, respectively). Weighted gene co-expression network analysis (WGCNA) and deconvolution analysis identified AD-associated gene modules and blood cell types. Regression and unsupervised clustering analysis identified AD-associated genes, gene modules, cell types, and established AD classification models. RESULTS: WGCNA on differentially expressed genes revealed 15 gene modules, with 6 accurately classifying AD (areas under the receiver operating characteristics curve [auROCs] > 0.90). These modules stratified AD patients into subgroups with distinct disease states. Cell-type deconvolution analysis identified specific blood cell types potentially associated with AD pathogenesis. DISCUSSION: This study highlights the potential of blood transcriptome for AD diagnosis, patient stratification, and mechanistic studies. HIGHLIGHTS: We comprehensively analyze the blood transcriptomes of a well-characterized Alzheimer's disease cohort to identify genes, gene modules, pathways, and specific blood cells associated with the disease. Blood transcriptome analysis accurately classifies and stratifies patients with Alzheimer's disease, with some gene modules achieving classification accuracy comparable to that of the plasma ATN biomarkers. Immune-associated pathways and immune cells, such as neutrophils, have potential roles in the pathogenesis and progression of Alzheimer's disease.

Contamination and Carcinogenic Risks of Lead in Lip Cosmetics in China.

There are growing concerns about elevated lead (Pb) levels in lip cosmetics, yet in China, the largest lip cosmetic market, recent Pb contamination in lip cosmetics and associated Pb exposure remain unclear. Here, we measured Pb levels of 29 popular lip cosmetics in China and conducted the bioaccessibility-corrected carcinogenic risk assessments and sensitivity analysis regarding Pb exposure for consumers using Monte Carlo simulation. The Pb concentrations of collected samples ranged from undetectable (< 0.05 µg/kg) to 0.21 mg/kg, all of which were well below the Pb concentration limit set for cosmetics in China (10 mg/kg). The 50th percentile incremental lifetime cancer risk (ILCR) of Pb in Chinese cosmetics (1.20E-07) was below the acceptable level (1E-06), indicating that the application of lip cosmetics and subsequent Pb exposure does not pose carcinogenic risks to consumers in most cases. The results of this study provide new insights into understanding the Pb risk in lip cosmetics.

AtHDA6 functions as an H3K18ac eraser to maintain pericentromeric CHG methylation in Arabidopsis thaliana.

Pericentromeric DNA, consisting of high-copy-number tandem repeats and transposable elements, is normally silenced through DNA methylation and histone modifications to maintain chromosomal integrity and stability. Although histone deacetylase 6 (HDA6) has been known to participate in pericentromeric silencing, the mechanism is still yet unclear. Here, using whole genome bisulfite sequencing (WGBS) and chromatin immunoprecipitation-sequencing (ChIP-Seq), we mapped the genome-wide patterns of differential DNA methylation and histone H3 lysine 18 acetylation (H3K18ac) in wild-type and hda6 mutant strains. Results show pericentromeric CHG hypomethylation in hda6 mutants was mediated by DNA demethylases, not by DNA methyltransferases as previously thought. DNA demethylases can recognize H3K18ac mark and then be recruited to the chromatin. Using biochemical assays, we found that HDA6 could function as an 'eraser' enzyme for H3K18ac mark to prevent DNA demethylation. Oxford Nanopore Technology Direct RNA Sequencing (ONT DRS) also revealed that hda6 mutants with H3K18ac accumulation and CHG hypomethylation were shown to have transcriptionally active pericentromeric DNA.

Transcriptomics integrated with metabolomics reveals the effect of Lycium barbarum polysaccharide on apoptosis in Nile tilapia (Oreochromis niloticus).

Lycium barbarum polysaccharide (LBP) is one of the main active ingredients in the fruit of L. barbarum L. It has been used as herbal medicine for thousands of years in China. In this study, Nile tilapia (Oreochromis niloticus) was taken as the research object. After feeding tilapia with 5 different doses of LBP (0 mg/kg, 500 mg/kg, 1000 mg/kg, 1500 mg/kg, 2000 mg/kg) for 55 d, it was found that LBP could promote the growth of tilapia, and this effect was the strongest at Group 1500 mg/kg. Apoptosis analysis in the liver and spleen showed that dietary supplementation with 1000 mg/kg LBP had the best protective effect on the spleen and liver in tilapia. Combined transcriptomics and metabolomics of the spleen in tilapia at Group 0 mg/kg and 1000 mg/kg showed that the differentially expressed genes (DEGs) such as NT5C2L1, pmm1, FasL and the differentially metabolites such as xanthine, dGMP, guanine and glutamate were mainly concentrated in signaling pathways such as Purine metabolism and FoxO signaling pathway. In conclusion, LBP regulates the metabolic waste levels of tilapia mainly through Purine metabolism and the FoxO signaling pathway, thereby inhibiting cell apoptosis, improving the utilization of nutrients, and promoting the growth of tilapia. This study not only provides a theoretical basis for the application of LBP in aquatic animals but also provides useful information for the healthy development of the aquaculture.

Arabidopsis EXTRA-LARGE G PROTEIN 1 (XLG1) functions together with XLG2 and XLG3 in PAMP-triggered MAPK activation and immunity.

Pattern-triggered immunity (PTI) is an essential strategy used by plants to deploy broad-spectrum resistance against pathogen attacks. Heterotrimeric G proteins have been reported to contribute to PTI. Of the three non-canonical EXTRA-LARGE G PROTEINs (XLGs) in Arabidopsis thaliana, XLG2 and XLG3 were shown to positively regulate immunity, but XLG1 was not considered to function in defense, based on the analysis of a weak xlg1 allele. In this study, we characterized the xlg1 xlg2 xlg3 triple knockout mutants generated from an xlg1 knockout allele. The strong xlg1 xlg2 xlg3 triple mutants compromised pathogen-associated molecular pattern (PAMP)-triggered activation of mitogen-activated protein kinases (MAPKs) and resistance to pathogen infection. The three XLGs interacted with MAPK cascade proteins involved in defense signaling, including the MAPK kinase kinases MAPKKK3 and MAPKKK5, the MAPK kinases MKK4 and MKK5, and the MAPKs MPK3 and MPK6. Expressing a constitutively active form of MKK4 restored MAPK activation and partially recovered the compromised disease resistance seen in the strong xlg1 xlg2 xlg3 triple mutant. Furthermore, mutations of all three XLGs largely restored the phenotype of the autoimmunity mutant bak1-interacting receptor-like kinase 1. Our study reveals that all three XLGs function redundantly in PAMP-triggered MAPK activation and plant immunity.

Arabidopsis PUB2 and PUB4 connect signaling components of pattern-triggered immunity.

Plants use pattern recognition receptors (PRRs) to detect pathogen-associated molecular patterns (PAMPs) and activate pattern-triggered immunity (PTI). Precise regulation of information from PRRs to downstream signaling components is vital to mounting an appropriate immune response and requires dynamic interactions of these PTI components. We used transcriptome profiling, phenotypic analysis, molecular genetics, and protein-protein interaction analysis to understand the roles of the Arabidopsis plant U-box (PUB) proteins PUB2 and PUB4 in disease resistance and PTI signaling. Loss of function of both PUB2 and PUB4 diminishes the PAMP-triggered oxidative bursts and dampens mitogen-activated protein kinase signaling, resulting in a severe compromise in resistance to not only pathogenic but also nonpathogenic strains of Pseudomonas syringae. Within PUB4, the E3 ligase activity is dispensable, but the armadillo repeat region is essential and sufficient for its function in immunity. PUB2 and PUB4 interact with PTI signaling components, including FLS2, BIK1, PBL27, and RbohD, and enhance FLS2-BIK1 and BIK1-RbohD interactions. Our study reveals that PUB2 and PUB4 are critical components of plant immunity and connect PTI components to positively regulate defense responses.

Taurine inhibits hydrogen peroxide-induced oxidative stress, inflammatory response and apoptosis in liver of Monopterus albus.

Fish are extremely vulnerable to environmental stimulation and produce oxidative stress. Among them, hydrogen peroxide is an oxidative stress source that cannot be ignored in fish, which can cause physical disorders, inflammation and even death. Taurine was revealed to reduce oxidative damage and inflammation caused by toxic substances, but whether it can reduce toxicity of rice field eel caused by H2O2 has not been determined. Thus, the intervention effects of taurine on H2O2-induced oxidative stress, inflammation, apoptosis, and autophagy in rice field eel. The results showed that oxidative injury in the liver was determined after H2O2 injection, as indicated by enhanced serum AST and ALT activities, inhibited the antioxidant function (increased MDA and ROS contents, decreased antioxidant enzymes, inhibited nrf2 transcription level), and induced inflammatory response (upregulated il-1ß, il-6, il-8, and il-12ß gene expression, downregulated tgf-ß1 gene expression, activated the transcription level of nf-κb, tlr-3, and tlr-7). In addition, bax, caspase3, beclin1, and Lc3B gene expression were significantly upregulated after H2O2 injection, while bcl2 and p62 gene expression were downregulated, leading to the occurrence of apoptosis and autophagy. In contrast, adding 0.2 and 0.5% taurine to feed significantly alleviated this damage, as indicated by the recovery of the aforementioned bioindicators, and the effect of 0.5% taurine addition is better than 0.2%. Overall, these results suggested that taurine can relieve the liver toxicity induced by H2O2, which enriched the toxic mechanism of H2O2 on fish and provided evidence for the protective effect of taurine on liver.

A novel interleukin-1 receptor-associated kinase 4 from blunt snout bream (Megalobrama amblycephala) is involved in inflammatory response via MyD88-mediated NF-κB signal pathway.

Interleukin-1 receptor-associated kinase 4 (IRAK4) plays a crucial role in the Toll-like receptor/IL-1R signal pathway, which mediates the downstream signal transduction involved in innate and adaptive immunity. In the present study, an IRAK4 homologue (named as MaIRAK4) from blunt snout bream (Megalobrama amblycephala) was cloned and characterized. The open reading frame (ORF) of MaIRAK4 contains 1422 nucleotides, encoding a putative protein of 473 amino acids. Protein structural analysis revealed that MaIRAK4 has an N-terminal death domain (DD) and a central kinase domain (S_TKc), similar to those of mammals and other fishes. Multiple sequence alignment demonstrated that MaIRAK4 is highly homologous with that of grass carp (97.67%). The qRT-PCR analysis showed that MaIRAK4 expressed widely in all examined tissues, including heart, liver, spleen, kidney, head-kidney, gill, intestine and muscle, with the highest expression in the liver and spleen. After stimulation with LPS, MaIRAK4 expression upregulated significantly and reached a peak at 6 h and 12 h post LPS stimulation in the spleen and head-kidney, respectively. After challenge with Aeromonas hydrophila, MaIRAK4 expression peaked at 48 h and 72 h in spleen/head-kidney and liver, respectively. These results implied that MaIRAK4 is involved in the host defense against bacterial infection. Subcellular localization analysis indicated that MaIRAK4 distributed in the cytoplasm. Co-immunoprecipitation and subcellular co-localization assay revealed that MaIRAK4 can combine with MaMyD88 through DD domain. MaIRAK4 overexpression can induce slightly the NF-κB promoter activity in HEK 293 cells. However, the activity of NF-κB promoter was dramatically enhanced after co-transfection with MaIRAK4 and MaMyD88 plasmids. The results showed that MaIRAK4 was involved in NF-κB signal pathway mediated by maMyD88. The expression level of pro-inflammatory cytokines (IL-1ß, IL-6, IL-8 and TNF-α) decreased significantly after the siRNA-mediated knockdown of MaIRAK4. Together, these results suggest that MaIRAK4 plays an important function in the innate immunity of M. amblycephala by inducing cytokines expression.

Looping Mercury Cycle in Global Environmental-Economic System Modeling.

The Minamata Convention on Mercury calls for Hg control actions to protect the environment and human beings from the adverse impacts of Hg pollution. It aims at the entire life cycle of Hg. Existing studies on the Hg cycle in the global environmental-economic system have characterized the emission-to-impact pathway of Hg pollution. That is, Hg emissions/releases from the economic system can have adverse impacts on human health and ecosystems. However, current modeling of the Hg cycle is not fully looped. It ignores the feedback of Hg-related environmental impacts (including human health impacts and ecosystem impacts) to the economic system. This would impede the development of more comprehensive Hg control actions. By synthesizing recent information on Hg cycle modeling, this critical review found that Hg-related environmental impacts would have feedbacks to the economic system via the labor force and biodiversity loss. However, the interactions between Hg-related activities in the environmental and economic systems are not completely clear. The cascading effects of Hg-related environmental impacts to the economic system throughout global supply chains have not been revealed. Here, we emphasize the knowledge gaps and propose possible approaches for looping the Hg cycle in global environmental-economic system modeling. This progress is crucial for formulating more dynamic and flexible Hg control measures. It provides new perspectives for the implementation of the Minamata Convention on Mercury.

NAD tagSeq reveals that NAD+-capped RNAs are mostly produced from a large number of protein-coding genes in Arabidopsis.

The 5' end of a eukaryotic mRNA transcript generally has a 7-methylguanosine (m7G) cap that protects mRNA from degradation and mediates almost all other aspects of gene expression. Some RNAs in Escherichia coli, yeast, and mammals were recently found to contain an NAD+ cap. Here, we report the development of the method NAD tagSeq for transcriptome-wide identification and quantification of NAD+-capped RNAs (NAD-RNAs). The method uses an enzymatic reaction and then a click chemistry reaction to label NAD-RNAs with a synthetic RNA tag. The tagged RNA molecules can be enriched and directly sequenced using the Oxford Nanopore sequencing technology. NAD tagSeq can allow more accurate identification and quantification of NAD-RNAs, as well as reveal the sequences of whole NAD-RNA transcripts using single-molecule RNA sequencing. Using NAD tagSeq, we found that NAD-RNAs in Arabidopsis were produced by at least several thousand genes, most of which are protein-coding genes, with the majority of these transcripts coming from <200 genes. For some Arabidopsis genes, over 5% of their transcripts were NAD capped. Gene ontology terms overrepresented in the 2,000 genes that produced the highest numbers of NAD-RNAs are related to photosynthesis, protein synthesis, and responses to cytokinin and stresses. The NAD-RNAs in Arabidopsis generally have the same overall sequence structures as the canonical m7G-capped mRNAs, although most of them appear to have a shorter 5' untranslated region (5' UTR). The identification and quantification of NAD-RNAs and revelation of their sequence features can provide essential steps toward understanding the functions of NAD-RNAs.

Transcriptomic evidence of luteinizing hormone-releasing hormone agonist (LHRH-A) regulation on lipid metabolism in grass carp (Ctenopharyngodon idella).

In this study, RNA sequencing was used to identify the hepatic gene expression profile in grass carp associated with luteinizing hormone-releasing hormone agonist (LHRH-A) treatment. A total of 93,912,172 reads were generated by HiSeq 4000 sequencing platform. After filtering, 83,450,860 clean reads were mapped to the reference genome. By calculating the FPKM of genes, 1475 differentially expressed genes were identified. PPAR signaling pathway was enriched with upregulated genes in LHRH-A injection group showing the regulation of the lipid metabolism by LHRH-A. The expression of eight key genes in PPAR signaling pathway was confirmed by qPCR and the results suggested that ACSL4A, ACSL4B, ANGPTL4, LPL, RXRBA and SLC27A1B were significantly stimulated by LHRH-A injection. This investigation provides the evidence that LHRH-A could play a role in lipid metabolism.

Chromosome-level assembly of the Hypophthalmichthys molitrix (Cypriniformes: Cyprinidae) genome provides insights into its ecological adaptation.

Hypophthalmichthys molitrix (silver carp) is phytoplanktivorous and is an economically and ecologically important fish species. As a well-known invasive species, a number of factors associated with the ecological adaptations of this species are largely unknown. Here, we present a chromosomal-level assembly of the species based on the PacBio Sequel II platform and Hi-C scaffolding technology. Based on the high-quality genome sequences and previous genome sequencing projects, a number of genes that were probably subject to positive selection reside in the genome of H. molitrix, and the last common ancestors of H. molitrix and H. nobilis were identified. Some of these genes may partially explain the mechanisms of H. molitrix for surviving damage due to toxic algae. Demographic history estimation suggests that the effective population size (EPS) of the species may have constantly increased along with the uplift of the Qinghai-Tibet Plateau, started to decline when quaternary glaciation started, and further declined during the Younger Dryas Period. Moreover, the introgression from H. nobilis to H. molitrix in North America was corroborated based on the whole-genome sequencing data, and the proportion of introgressed regions was estimated to be approximately 5.8%. Based on the high-quality assembly, the possible mechanisms by which H. molitrix adapts to its endemic and invaded locations were profiled.

Low expression of miR-19a-5p is associated with high mRNA expression of diacylglycerol O-acyltransferase 2 (DGAT2) in hybrid tilapia.

DGAT2 (acyl CoA:diacylglycerol acyltransferase 2) is a key and rate-limiting enzyme that catalyzes the final step of triglyceride (TG) synthesis. In this study, hybrid tilapia were generated from Nile tilapia (♀) and blue tilapia (♂) crossing. The TG content levels in the liver of these tilapia were measured. The results showed that the TG content was higher in the hybrid tilapia. In addition, protein and mRNA expression levels in the tilapia livers were determined. Higher hepatic mRNA and protein expression of DGAT2 in the hybrid fish was found. A luciferase reporter assay with HEK293T cells revealed that miRNA-19a-5p targeted the 3'UTR of DGAT2, suggesting a direct regulatory mechanism. Using qRT-PCR, we found that DGAT2 mRNA levels had a negative correlation with miRNA-19a-5p expression in Nile tilapia and hybrid. Taken together, these findings provide evidence that miRNA-19a-5p is involved in TG synthesis in the regulation of lipid metabolism in tilapia.

Reply to Comments on "A Scientometric Analysis of Recent Literature on Arsenic Bioaccumulation and Biotransformation in Marine Ecosystems".

In this letter, we formally replied to the comments on our recently published paper "A scientometric analysis of recent literature on arsenic bioaccumulation and biotransformation in marine ecosystems" with evidence and facts.

Assessment of the Bioavailability of Mercury Sulfides in Paddy Soils Using Sodium Thiosulfate Extraction - Results from Microcosm Experiments.

Mercury sulfides (HgS), one of the largest Hg sinks in the lithosphere, has long been considered to be highly inert. Recently, several HgS speciation (e.g., nano- or micro-sized HgS particles) in paddy soils have been found to be reactive and bioavailable, increasing the possibility of methylation and bioaccumulation and posing a potential risk to humans. However, a simple and uniform method for investigating HgS bioavailability is still lacking. To address this issue, we extracted dissolved Hg from HgS particles by sodium thiosulfate (Na2S2O3) in paddy soils and analyzed the correlation between extracted Hg and soil methylmercury (MeHg). Results showed that the amounts of Hg extracted by Na2S2O3 had a strong positive correlation with the levels of soil MeHg (R 2 adj = 0.893, p < 0.05). It is suggested that Na2S2O3 extraction may be a good method of predicting Hg bioavailability in paddy soils. Our results would help to give clues in better predicting Hg risk in natural environments.