Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Time-resolved single-cell transcriptomics defines immune trajectories in glioblastoma.

Deciphering the cell-state transitions underlying immune adaptation across time is fundamental for advancing biology. Empirical in vivo genomic technologies that capture cellular dynamics are currently lacking. We present Zman-seq, a single-cell technology recording transcriptomic dynamics across time by introducing time stamps into circulating immune cells, tracking them in tissues for days. Applying Zman-seq resolved cell-state and molecular trajectories of the dysfunctional immune microenvironment in glioblastoma. Within 24 hours of tumor infiltration, cytotoxic natural killer cells transitioned to a dysfunctional program regulated by TGFB1 signaling. Infiltrating monocytes differentiated into immunosuppressive macrophages, characterized by the upregulation of suppressive myeloid checkpoints Trem2, Il18bp, and Arg1, over 36 to 48 hours. Treatment with an antagonistic anti-TREM2 antibody reshaped the tumor microenvironment by redirecting the monocyte trajectory toward pro-inflammatory macrophages. Zman-seq is a broadly applicable technology, enabling empirical measurements of differentiation trajectories, which can enhance the development of more efficacious immunotherapies.

LGR5 expressing skin fibroblasts define a major cellular hub perturbed in scleroderma.

Systemic sclerosis (scleroderma, SSc) is an incurable autoimmune disease with high morbidity and mortality rates. Here, we conducted a population-scale single-cell genomic analysis of skin and blood samples of 56 healthy controls and 97 SSc patients at different stages of the disease. We found immune compartment dysfunction only in a specific subtype of diffuse SSc patients but global dysregulation of the stromal compartment, particularly in a previously undefined subset of LGR5+-scleroderma-associated fibroblasts (ScAFs). ScAFs are perturbed morphologically and molecularly in SSc patients. Single-cell multiome profiling of stromal cells revealed ScAF-specific markers, pathways, regulatory elements, and transcription factors underlining disease development. Systematic analysis of these molecular features with clinical metadata associates specific ScAF targets with disease pathogenesis and SSc clinical traits. Our high-resolution atlas of the sclerodermatous skin spectrum will enable a paradigm shift in the understanding of SSc disease and facilitate the development of biomarkers and therapeutic strategies.

Coupled scRNA-Seq and Intracellular Protein Activity Reveal an Immunosuppressive Role of TREM2 in Cancer.

Cell function and activity are regulated through integration of signaling, epigenetic, transcriptional, and metabolic pathways. Here, we introduce INs-seq, an integrated technology for massively parallel recording of single-cell RNA sequencing (scRNA-seq) and intracellular protein activity. We demonstrate the broad utility of INs-seq for discovering new immune subsets by profiling different intracellular signatures of immune signaling, transcription factor combinations, and metabolic activity. Comprehensive mapping of Arginase 1-expressing cells within tumor models, a metabolic immune signature of suppressive activity, discovers novel Arg1+ Trem2+ regulatory myeloid (Mreg) cells and identifies markers, metabolic activity, and pathways associated with these cells. Genetic ablation of Trem2 in mice inhibits accumulation of intra-tumoral Mreg cells, leading to a marked decrease in dysfunctional CD8+ T cells and reduced tumor growth. This study establishes INs-seq as a broadly applicable technology for elucidating integrated transcriptional and intra-cellular maps and identifies the molecular signature of myeloid suppressive cells in tumors.

Lipid-Associated Macrophages Control Metabolic Homeostasis in a Trem2-Dependent Manner.

Immune cells residing in white adipose tissue have been highlighted as important factors contributing to the pathogenesis of metabolic diseases, but the molecular regulators that drive adipose tissue immune cell remodeling during obesity remain largely unknown. Using index and transcriptional single-cell sorting, we comprehensively map all adipose tissue immune populations in both mice and humans during obesity. We describe a novel and conserved Trem2+ lipid-associated macrophage (LAM) subset and identify markers, spatial localization, origin, and functional pathways associated with these cells. Genetic ablation of Trem2 in mice globally inhibits the downstream molecular LAM program, leading to adipocyte hypertrophy as well as systemic hypercholesterolemia, body fat accumulation, and glucose intolerance. These findings identify Trem2 signaling as a major pathway by which macrophages respond to loss of tissue-level lipid homeostasis, highlighting Trem2 as a key sensor of metabolic pathologies across multiple tissues and a potential therapeutic target in metabolic diseases.

Understanding m6A changes in chromophobe renal cell carcinoma and predicting patient outcomes survival.

N6-methyladenosine (m6A) is a prevalent mRNA modification known for its implications in various cancer types, yet its role in chromophobe renal cell carcinoma (chRCC) remains largely unexplored. In this study, we performed m6A-SEAL-seq and RNA-seq analyses on tissues from three chRCC subjects, aiming to uncover m6A alterations in chRCC. Our findings revealed reduced expression levels of four m6A regulators in chRCC tissues and highlighted differences in m6A levels compared to normal tissues. Furthermore, we identified specific genes and cancer-related pathways affected by these differences, including notable candidates like NOTCH1 and FGFR1, implicated in chRCC development. Additionally, we developed a predictive model based on the expression level of m6A associated genes, demonstrating promising prognostic capabilities for patient survival prediction. Overall, our study provides valuable insights into the role of m6A in chRCC and its potential as a prognostic indicator.

Advancing primate surveillance with image recognition techniques from unmanned aerial vehicles.

Using unmanned aerial vehicles (UAVs) for surveys on thermostatic animals has gained prominence due to their ability to provide practical and precise dynamic censuses, contributing to developing and refining conservation strategies. However, the practical application of UAVs for animal monitoring necessitates the automation of image interpretation to enhance their effectiveness. Based on our past experiences, we present the Sichuan snub-nosed monkey (Rhinopithecus roxellana) as a case study to illustrate the effective use of thermal cameras mounted on UAVs for monitoring monkey populations in Qinling, a region characterized by magnificent biodiversity. We used the local contrast method for a small infrared target detection algorithm to collect the total population size. Through the experimental group, we determined the average optimal grayscale threshold, while the validation group confirmed that this threshold enables automatic detection and counting of target animals in similar datasets. The precision rate obtained from the experiments ranged from 85.14% to 97.60%. Our findings reveal a negative correlation between the minimum average distance between thermal spots and the count of detected individuals, indicating higher interference in images with closer thermal spots. We propose a formula for adjusting primate population estimates based on detection rates obtained from UAV surveys. Our results demonstrate the practical application of UAV-based thermal imagery and automated detection algorithms for primate monitoring, albeit with consideration of environmental factors and the need for data preprocessing. This study contributes to advancing the application of UAV technology in wildlife monitoring, with implications for conservation management and research.

Advancements in enhancing resistant starch type 3 (RS3) content in starchy food and its impact on gut microbiota: A review.

Resistant starch type 3 (RS3), often found in cooked starchy food, has various health benefits due to its indigestible properties and physiological functions such as promoting the abundance of gut beneficial microbial flora and inhibiting the growth of intestinal pathogenic bacteria. However, it is challenging to develop starchy food with high RS3 content. This review aims to provide a detailed overview of current advancements to enhance RS3 content in starchy food and its effects of RS3 on gut microbiota. These approaches include breeding high-amylose cereals through gene editing techniques, processing, enzyme treatments, storage, formation of RS3 nanoparticles, and the incorporation of bioactive compounds. The mechanisms, specific conditions, advantages, and disadvantages associated with each approach and the potential effects of RS3 prepared by different methods on gut microbiota are summarized. In conclusion, this review contains important information that aims to provide guidelines for developing an efficient RS3 preparation process and promote the consumption of RS3-enriched starchy foods to improve overall health outcomes.

Social and paternal female choice for male MHC genes in golden snub-nosed monkeys (Rhinopithecus roxellana).

Genes of the major histocompatibility complex (MHC) play important roles in vertebrate immunocompetence. MHC genes thus offer females indirect benefits to mate choice through the production of offspring of an optimal MHC genotype. Females may choose males with specific MHC haplotypes, dissimilar MHC genotypes, MHC heterozygous males or MHC-diverse males. We tested these four alternatives for both female social and paternal choice in wild golden snub-nosed monkeys (Rhinopithecus roxellana) by examining overall genetic variability (via microsatellites) and four MHC-genes (DRB1, DRB2, DQA1 and DQB1). Monte Carlo randomization tests showed that MHC dissimilarity was favoured for social choice (males to which females were socially affiliated) and intermediate MHC dissimilarity was favoured in paternal choice (fathers of offspring). No evidence of inbreeding avoidance was found for either social or paternal mates. We found that MHC heterozygotes, higher microsatellite multilocus heterozygosity and higher microsatellites diversity were favoured for social mates, and higher microsatellite diversity was favoured for paternal mates. Independent of male age, we found that the formation of male-female social pairings is significantly predicted by compatibility based on the sharing of MHC haplotypes. However, we found no evidence of independent genetic effects on the duration of male-female social pairings, male social status (achieving OMU leader male status or not), the number of females with which individual leader males paired, the likelihood of potential male-female pairings producing offspring, or whether males fathered offspring or not. Overall, our findings suggest different genetic factors are involved in social and paternal choice in R. roxellana.

Inducing Inorganic Carbon Accrual in Subsoil through Biochar Application on Calcareous Topsoil.

Biochar amendments add persistent organic carbon to soil and can stabilize rhizodeposits and existing soil organic carbon (SOC), but effects of biochar on subsoil carbon stocks have been overlooked. We quantified changes in soil inorganic carbon (SIC) and SOC to 2 m depth 10 years after biochar application to calcareous soil. The total soil carbon (i.e., existing SOC, SIC, and biochar-C) increased by 71, 182, and 210% for B30, B60, and B90, respectively. Biochar application at 30, 60, and 90 t ha-1 rates significantly increased SIC by 10, 38, and 68 t ha-1, respectively, with accumulation mainly occurring in the subsoil (below 1 m). This huge increase of SIC (mainly CaCO3) is ∼100 times larger than the inorganic carbon present in the added biochar (0.3, 0.6, or 0.9 t ha-1). The benzene polycarboxylic acid method showed that the biochar-amended soil contained more black carbon particles (6.8 times higher than control soil) in the depth of 1.4-1.6 m, which provided the direct quantitative evidence for biochar migration into subsoil after a decade. Spectral and energy spectrum analysis also showed an obvious biochar structure in the biochar-amended subsoil, accompanied by a Ca/Mg carbonate cluster, which provided further evidence for downward migration of biochar after a decade. To explain SIC accumulation in subsoil with biochar amendment, the interacting mechanisms are proposed: (1) biochar amendment significantly increases subsoil pH (0.3-0.5 units) 10 years after biochar application, thus forming a favorable pH environment in the subsoil to precipitate HCO3-; and (2) the transported biochar in subsoil can act as nuclei to precipitate SIC. Biochar amendment enhanced SIC by up to 80%; thus, the effects on carbon stocks in subsoil must be understood to inform strategies for carbon dioxide removal through biochar application. Our study provided critical knowledge on the impact of biochar application to topsoil on carbon stocks in subsoil in the long term.

Consistency of limb preference across unimanual feeding, bipedal locomotion, and social grooming in golden snub-nosed monkeys (Rhinopithecus roxellana).

The golden snub-nosed monkey (Rhinopithecus roxellana) is a typical arboreal group-living Old World primate. While limb preference has been extensively studied in this species, limb preference consistency has not yet been explored. Here, based on 26 R. roxellana adults, we investigated whether individuals exhibit consistent motor biases in manual- (e.g., unimanual feeding and social grooming) and foot-related (e.g., bipedal locomotion) tasks and whether limb preference consistency is influenced by increased social interactions during social grooming. Results showed no consistency in the direction or strength of limb preference among tasks, except for lateral strength in handedness for unimanual feeding and footedness in the initiation of locomotion. Population-level foot preference was only found among right-handers. Marked lateral bias was found in unimanual feeding, indicating that it may be a sensitive behavioural measure for assessing manual preference, especially for provisioned populations. This study not only improves our understanding of the relationship between hand and foot preference in R. roxellana but also reveals potential differential hemispheric regulation of limb preference and the influence of increased social interaction on handedness consistency.

Comparison of Depth-Specific Prediction of Soil Properties: MIR vs. Vis-NIR Spectroscopy.

The prediction of soil properties at different depths is an important research topic for promoting the conservation of black soils and the development of precision agriculture. Mid-infrared spectroscopy (MIR, 2500-25000 nm) has shown great potential in predicting soil properties. This study aimed to explore the ability of MIR to predict soil organic matter (OM) and total nitrogen (TN) at five different depths with the calibration from the whole depth (0-100 cm) or the shallow layers (0-40 cm) and compare its performance with visible and near-infrared spectroscopy (vis-NIR, 350-2500 nm). A total of 90 soil samples containing 450 subsamples (0-10 cm, 10-20 cm, 20-40 cm, 40-70 cm, and 70-100 cm depths) and their corresponding MIR and vis-NIR spectra were collected from a field of black soil in Northeast China. Multivariate adaptive regression splines (MARS) were used to build prediction models. The results showed that prediction models based on MIR (OM: RMSEp = 1.07-3.82 g/kg, RPD = 1.10-5.80; TN: RMSEp = 0.11-0.15 g/kg, RPD = 1.70-4.39) outperformed those based on vis-NIR (OM: RMSEp = 1.75-8.95 g/kg, RPD = 0.50-3.61; TN: RMSEp = 0.12-0.27 g/kg; RPD = 1.00-3.11) because of the higher number of characteristic bands. Prediction models based on the whole depth calibration (OM: RMSEp = 1.09-2.97 g/kg, RPD = 2.13-5.80; TN: RMSEp = 0.08-0.19 g/kg, RPD = 1.86-4.39) outperformed those based on the shallow layers (OM: RMSEp = 1.07-8.95 g/kg, RPD = 0.50-3.93; TN: RMSEp = 0.11-0.27 g/kg, RPD = 1.00-2.24) because the soil sample data of the whole depth had a larger and more representative sample size and a wider distribution. However, prediction models based on the whole depth calibration might provide lower accuracy in some shallow layers. Accordingly, it is suggested that the methods pertaining to soil property prediction based on the spectral library should be considered in future studies for an optimal approach to predicting soil properties at specific depths. This study verified the superiority of MIR for soil property prediction at specific depths and confirmed the advantage of modeling with the whole depth calibration, pointing out a possible optimal approach and providing a reference for predicting soil properties at specific depths.

Linkage disequilibrium under polysomic inheritance.

Linkage disequilibrium (LD) is the non-random association of alleles at different loci. Squared LD coefficients r2 (for phased genotypes) and [Formula: see text] (for unphased genotypes) will converge to constants that are determined by the sample size, the recombination frequency, the effective population size and the mating system. LD can therefore be used for gene mapping and the estimation of effective population size. However, current methods work only with diploids. To resolve this problem, we here extend the linkage disequilibrium measures to include polysomic inheritance. We derive the values of r2 and [Formula: see text] at equilibrium state for various mating systems and different ploidy levels. For unlinked loci, [Formula: see text] for monoecious and dioecious (with random pairing) mating systems or [Formula: see text] for dioecious mating systems (with lifetime pairing), where f is the number of females in a half-sib family and η is a constant related to the ploidy level. We simulate the application of estimating Ne using unphased genotypes. We find that estimating Ne in polyploids requires similar sample sizes and numbers of loci as in diploids, with the main source of bias due to using 0.5 as the recombination frequency.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs.

MAIN CONCLUSION: The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNAVal-UAC and tRNAIle-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNAMet seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

Estimating the contribution of plant traits to light partitioning in simultaneous maize/soybean intercropping.

Spatial configuration and plant phenotypic plasticity contribute to increased light capture in relay intercropping, but there is little information on whether these factors also increase light capture in simultaneous intercropping. We developed and validated a three-dimensional functional-structural plant model to simulate light capture in maize and soybean sole crops and intercrop scenarios, using species traits observed in sole crops and intercrops. The intercrop maize phenotype had 2% greater light capture than the sole crop phenotype in a pure stand. The soybean intercrop phenotype had 5-10% lower light capture than the sole crop phenotype in a pure stand. The intercrop configuration increased the light capture of maize by 29% and reduced the light capture of soybean by 42%, compared with the light capture expected from sole crops. However, intercrop configuration only marginally affected total light capture by the intercrop system (+1%). Testing of individual soybean plant traits revealed that plasticity in leaf dimensions was the main reason for differences in light capture by soybean in simulated sole crops and intercrops. The results of this study illustrate a major shift of light capture from shorter species (soybean) to the taller component (maize) in a simultaneous strip intercrop. Plastic plant traits modulate this overall effect, but only marginally.

Corrigendum: A Severe Lack of Evidence Limits Effective Conservation of the World's Primates.

[This corrects the article DOI: 10.1093/biosci/biaa082.].

The geometry of resource constraint: An empirical study of the golden snub-nosed monkey.

Apposite conceptualization and measurement of resource variation is critical for understanding many issues in ecology, including ecological niches, persistence and distribution of populations, the structure of communities and population resilience to perturbations. We apply the nutritional geometry framework to conceptualize and quantify the responses of a temperate-living primate, the golden snub-nosed monkey Rhinopithecus roxellana to variation in resource quality and quantity and in nutrient requirements associated with seasonal environments. We present a geometric model distinguishing qualitative constraint, quantitative constraint and 'pseudo-constraint' whereby nutrient intakes resemble response to qualitative resource constraint but are in fact driven by variation in nutrient requirements. The model is applied to analyse nutrient intakes recorded in 164 full-day observations of monkeys from two populations, one wild and the other captive, across seasons. Additionally, we recorded the diet of a single animal over 32 consecutive days in the wild. Despite considerable differences in available resources, the captive and wild populations showed marked similarities in nutrient intakes, including indistinguishable amounts and ratios of ingested macronutrients during summer and autumn and strong year-round maintenance of protein compared to seasonally variable fat and carbohydrate intakes. These similarities suggest homeostatically regulated nutritional targets and provide reference points to identify factors driving population differences in macronutrient intake in winter and spring. Our framework enabled us to distinguish examples of quantitative, qualitative and 'pseudo-constraint'. We suggest that this approach can increase the resolution at which resource constraint is conceptualized and measured in ecological studies.

Effects of small-molecule compounds on fibroblast properties in golden snub-nosed monkey (Rhinopithecus roxellana).

BACKGROUND: Golden snub-nosed monkey (Rhinopithecus roxellana) is an endangered primate species, whose molecular material for conservation purposes has not yet been maintained. Although small-molecule compounds (SMCs) have been reported to improve induced pluripotent stem cells (iPSCs), their efficiency in the interspecies-transferred nucleus is still unknown. METHODS: We thus used the fibroblasts from the golden snub-nosed monkey treated with SMC as donor cells, injected into the enucleated oocytes of goats, to test such efficiency. Gene expression profiles in the cell-constructed embryos with and without SMCs were compared by qPCR. RESULTS: The results show that cell morphology undergoes remarkable changes (volume is smaller than normal cells, and many black spots in the cytoplasm were found); pluripotent genes (Oct4, Sox2, and Nanog) significantly increased with SMC treatment. CONCLUSIONS: This study demonstrates that SMCs alter the properties of donor cells and promote the expression of pluripotent genes in hybrid embryos.

Use of historical and contemporary distribution of mammals in China to inform conservation.

A systematic understanding of dynamic animal extinction trajectories for different regions in a nation like China is critically important to developing practical conservation strategies. We explored historical and contemporary changes in terrestrial mammalian diversity to determine how diversity in each of the 5 regions in China has changed over time and to examine the conservation potential of these regions. We used records from databases on Pleistocene mammalian fossils and historical distribution records (1175-2020) for Primates (as a case study) to reconstruct evolutionary and historical distribution trajectories of the 11 orders of terrestrial mammals and to predict their prospective survival based on the national conservation strategy applied. The results indicated that since the Pleistocene, 4-5 mammalian orders have been lost in the northeast, 3 in central China, 2 along the coast, and 1 in the northwest. In the southwest, all 11 orders were maintained. Contemporarily, the coast and southwest had the highest and second-highest species densities. The southwest region and southeastern sections of the northwest region were the most historically and contemporarily diverse areas, which suggests that they should be the first priority for protected area (PA) designation. The central and coastal areas should be secondarily prioritized. In these 2 regions, conservation should focus on human coexistence with nature. Less attention should be paid to the PA in the northeast and western northwest because in these areas ecosystems are depauperate and the climate is harsh. Conservation in these areas should focus principally on avoiding further human encroachment on natural areas. Article impact statement: Historical and contemporary patterns of extinction can be a basis for mammalian conservation strategies.

Uso de la Distribución Histórica y Contemporánea de los Mamíferos en China para Orientar a la Conservación Resumen El entendimiento sistemático de las trayectorias de extinción dinámica de los animales para diferentes regiones en un país como en China es de importancia crítica para poder desarrollar estrategias de conservación prácticas. Exploramos los cambios históricos y contemporáneos en la diversidad de mamíferos terrestres para determinar cómo ha cambiado con el tiempo la diversidad en cada una de las cinco regiones de China y para examinar el potencial de conservación de aquellas regiones. Usamos registros tomados de bases de datos con información sobre los fósiles de mamíferos del Pleistoceno y registros de la distribución histórica (1175 - 2020) de los primates (como estudio de caso) para reconstruir las trayectorias de distribución históricas y contemporáneas de los once órdenes de mamíferos terrestres presentes en China. También utilizamos esta información para predecir la supervivencia potencial de estos órdenes con base en la estrategia nacional de conservación aplicada. Los resultados indicaron que, desde el Pleistoceno, entre cuatro y cinco órdenes de mamíferos se han perdido en el noreste de China, tres en la zona central, dos a lo largo de la costa y uno en el noroeste. En el suroeste del país se han mantenido los once órdenes. Actualmente, la costa y el suroeste tienen la primera y segunda densidad más alta de especies. La región del suroeste y las secciones al sureste de la región del noroeste son las áreas con la mayor diversidad histórica y contemporánea, lo que sugiere que deberían ser de primera prioridad para la designación de áreas protegidas (AP). Las áreas de la costa y el centro deberían ser las de segunda prioridad por la misma razón. En estas dos últimas regiones, la conservación debería enfocarse en la coexistencia entre los humanos y la naturaleza. Se les debería prestar una menor atención a las AP en el noreste y en la zona oeste del noroeste porque en estas áreas los ecosistemas son paupérrimos y el clima es hostil. La conservación en estas áreas debería enfocarse principalmente en evitar más invasiones humanas dentro de las áreas naturales.

Extracellular polymeric substances induced cell-surface interactions facilitate bacteria transport in saturated porous media.

Bacteria often respond to dynamic soil environment through the secretion of extracellular polymeric substances (EPS). The EPS modifies cell surface properties and soil pore-scale hydration status, which in turn, influences bacteria transport in soil. However, the effect of soil particle size and EPS-mediated surface properties on bacterial transport in the soil is not well understood. In this study, the simultaneous impacts of EPS and collector size on Escherichia coli (E. coli) transport and deposition in a sand column were investigated. E. coli transport experiments were carried out under steady-state flow in saturated columns packed with quartz sand with different size ranges, including 0.300-0.425 mm (sand-I), 0.212-0.300 mm (sand-II), 0.106-0.150 mm (sand-III) and 0.075-0.106 mm (sand-IV). Bacterial retention increased with decreasing sand collector size, suggesting that straining played an important role in fine-textured media. Both experiment and simulation results showed a clear drop in the retention rate of the bacterial population with the presence of additional EPS (200 mg L-1) (EPS+). The inhibited retention of cells in sand columns under EPS+ scenario was likely attributed to enhanced bacteria hydrophilicity and electrostatic repulsion between cells and sand particles as well as reduced straining. Calculations of the extended Derjaguin-Landau-Verwey-Overbeek (XDLVO) interactions energies revealed that high repulsive energy barrier existed between bacterial cells and sand particles in EPS+ environment, primarily due to high repulsive electrostatic force and Lewis acid-base force, as well as low attractive Lifshitz-van der Waals force, which retarded bacterial population deposition. Steric stabilization of EPS would also prevent the approaching of cells close to the quartz surface and thereby hinder cell attachment. This study was the first to show that EPS reduced bacterial straining in saturated porous media. These findings provide new insight into the functional effects of extrinsic EPS on bacterial transport behavior in the saturated soil environment, e.g., aquifers.