Nigrostriatal degeneration determines dynamics of glial inflammatory and phagocytic activity.

Glial cells are key players in the initiation of innate immunity in neurodegeneration. Upon damage, they switch their basal activation state and acquire new functions in a context and time-dependent manner. Since modulation of neuroinflammation is becoming an interesting approach for the treatment of neurodegenerative diseases, it is crucial to understand the specific contribution of these cells to the inflammatory reaction and to select experimental models that recapitulate what occurs in the human disease. Previously, we have characterized a region-specific activation pattern of CD11b+ cells and astrocytes in the α-synuclein overexpression mouse model of Parkinson´s disease (PD). In this study we hypothesized that the time and the intensity of dopaminergic neuronal death would promote different glial activation states. Dopaminergic degeneration was induced with two administration regimens of the neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP), subacute (sMPTP) and chronic (cMPTP). Our results show that in the sMPTP mouse model, the pro-inflammatory phenotype of striatal CD11b+ cells was counteracted by an anti-inflammatory astrocytic profile. In the midbrain the roles were inverted, CD11b+ cells exhibited an anti-inflammatory profile and astrocytes were pro-inflammatory. The overall response generated resulted in decreased CD4 T cell infiltration in both regions. Chronic MPTP exposure resulted in a mild and prolonged neuronal degeneration that generated a pro-inflammatory response and increased CD4 T cell infiltration in both regions. At the onset of the neurodegenerative process, microglia and astrocytes cooperated in the removal of dopaminergic terminals. With time, only microglia maintained the phagocytic activity. In the ventral midbrain, astrocytes were the main phagocytic mediators at early stages of degeneration while microglia were the major phagocytic cells in the chronic state. In this scenario, we questioned which activation pattern recapitulates better the features of glial activation in PD. Glial activation in the cMPTP mouse model reflects many pathways of their corresponding counterparts in the human brain with advanced PD. Altogether, our results point toward a context-dependent cooperativity of microglia/myeloid cells and astrocytes in response to neuronal damage and the relevance of selecting the right experimental models for the study of neuroinflammation.

Generation of heart and vascular system in rodents by blastocyst complementation.

Generating organs from stem cells through blastocyst complementation is a promising approach to meet the clinical need for transplants. In order to generate rejection-free organs, complementation of both parenchymal and vascular cells must be achieved, as endothelial cells play a key role in graft rejection. Here, we used a lineage-specific cell ablation system to produce mouse embryos unable to form both the cardiac and vascular systems. By mouse intraspecies blastocyst complementation, we rescued heart and vascular system development separately and in combination, obtaining complemented hearts with cardiomyocytes and endothelial cells of exogenous origin. Complemented chimeras were viable and reached adult stage, showing normal cardiac function and no signs of histopathological defects in the heart. Furthermore, we implemented the cell ablation system for rat-to-mouse blastocyst complementation, obtaining xenogeneic hearts whose cardiomyocytes were completely of rat origin. These results represent an advance in the experimentation towards the in vivo generation of transplantable organs.

Optimization of universal allogeneic CAR-T cells combining CRISPR and transposon-based technologies for treatment of acute myeloid leukemia.

Despite the potential of CAR-T therapies for hematological malignancies, their efficacy in patients with relapse and refractory Acute Myeloid Leukemia has been limited. The aim of our study has been to develop and manufacture a CAR-T cell product that addresses some of the current limitations. We initially compared the phenotype of T cells from AML patients and healthy young and elderly controls. This analysis showed that T cells from AML patients displayed a predominantly effector phenotype, with increased expression of activation (CD69 and HLA-DR) and exhaustion markers (PD1 and LAG3), in contrast to the enriched memory phenotype observed in healthy donors. This differentiated and more exhausted phenotype was also observed, and corroborated by transcriptomic analyses, in CAR-T cells from AML patients engineered with an optimized CAR construct targeting CD33, resulting in a decreased in vivo antitumoral efficacy evaluated in xenograft AML models. To overcome some of these limitations we have combined CRISPR-based genome editing technologies with virus-free gene-transfer strategies using Sleeping Beauty transposons, to generate CAR-T cells depleted of HLA-I and TCR complexes (HLA-IKO/TCRKO CAR-T cells) for allogeneic approaches. Our optimized protocol allows one-step generation of edited CAR-T cells that show a similar phenotypic profile to non-edited CAR-T cells, with equivalent in vitro and in vivo antitumoral efficacy. Moreover, genomic analysis of edited CAR-T cells revealed a safe integration profile of the vector, with no preferences for specific genomic regions, with highly specific editing of the HLA-I and TCR, without significant off-target sites. Finally, the production of edited CAR-T cells at a larger scale allowed the generation and selection of enough HLA-IKO/TCRKO CAR-T cells that would be compatible with clinical applications. In summary, our results demonstrate that CAR-T cells from AML patients, although functional, present phenotypic and functional features that could compromise their antitumoral efficacy, compared to CAR-T cells from healthy donors. The combination of CRISPR technologies with transposon-based delivery strategies allows the generation of HLA-IKO/TCRKO CAR-T cells, compatible with allogeneic approaches, that would represent a promising option for AML treatment.

Uncovering perturbations in human hematopoiesis associated with healthy aging and myeloid malignancies at single-cell resolution.

Early hematopoiesis is a continuous process in which hematopoietic stem and progenitor cells (HSPCs) gradually differentiate toward specific lineages. Aging and myeloid malignant transformation are characterized by changes in the composition and regulation of HSPCs. In this study, we used single-cell RNA sequencing (scRNA-seq) to characterize an enriched population of human HSPCs obtained from young and elderly healthy individuals.Based on their transcriptional profile, we identified changes in the proportions of progenitor compartments during aging, and differences in their functionality, as evidenced by gene set enrichment analysis. Trajectory inference revealed that altered gene expression dynamics accompanied cell differentiation, which could explain aging-associated changes in hematopoiesis. Next, we focused on key regulators of transcription by constructing gene regulatory networks (GRNs) and detected regulons that were specifically active in elderly individuals. Using previous findings in healthy cells as a reference, we analyzed scRNA-seq data obtained from patients with myelodysplastic syndrome (MDS) and detected specific alterations of the expression dynamics of genes involved in erythroid differentiation in all patients with MDS such as TRIB2. In addition, the comparison between transcriptional programs and GRNs regulating normal HSPCs and MDS HSPCs allowed identification of regulons that were specifically active in MDS cases such as SMAD1, HOXA6, POU2F2, and RUNX1 suggesting a role of these transcription factors (TFs) in the pathogenesis of the disease.In summary, we demonstrate that the combination of single-cell technologies with computational analysis tools enable the study of a variety of cellular mechanisms involved in complex biological systems such as early hematopoiesis and can be used to dissect perturbed differentiation trajectories associated with perturbations such as aging and malignant transformation. Furthermore, the identification of abnormal regulatory mechanisms associated with myeloid malignancies could be exploited for personalized therapeutic approaches in individual patients.

Our blood contains many different types of cells; red blood cells carry oxygen through the body, platelets help to stop bleeding and a variety of white blood cells fight infections. All of these critical components come from a pool of immature cells in bone marrow, which can develop and specialise into any of these. However, as we get older, these immature cells can accumulate damage, including mutations in specific genes. This increases the risk of diseases such as myelodysplastic syndromes (MDS), a type of cancer in which the cells cannot develop and the patient does not have enough healthy mature blood cells. The changes in gene activity in the immature cells have previously been studied using samples from young and elderly people, as well as individuals with MDS. These studies examined large numbers of cells together, revealing differences between young and elderly people, and individuals with MDS. However, this does not describe how the different types alter their behaviour. To address this, Ainciburu, Ezponda et al. used a technique called single-cell RNA sequencing to study the gene activity in individual immature blood cells. This revealed changes associated with maturation that may account for the different combinations of cell populations in younger and older people. The results confirmed findings from previous studies and suggested new genes involved in ageing or MDS. Ainciburu, Ezponda et al. used these results to create an analytical system that highlights gene activity differences in individual MDS patients that are independent of age-related changes. These results provide new insights that could help further research into the development of MDS and the ageing process. In addition, scientists could study other diseases using this approach of analysing individual patients' gene activity. In future, this could help to personalise clinical decisions on diagnosis and treatment.

The transcription factor DDIT3 is a potential driver of dyserythropoiesis in myelodysplastic syndromes.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in older individuals. Here we analyze the transcriptome of human HSCs purified from young and older healthy adults, as well as MDS patients, identifying transcriptional alterations following different patterns of expression. While aging-associated lesions seem to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detect the upregulation of the transcription factor DNA Damage Inducible Transcript 3 (DDIT3). Overexpression of DDIT3 in human healthy HSCs induces an MDS-like transcriptional state, and dyserythropoiesis, an effect associated with a failure in the activation of transcriptional programs required for normal erythroid differentiation. Moreover, DDIT3 knockdown in CD34+ cells from MDS patients with anemia is able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythroid differentiation characterizing MDS patients.

CAR density influences antitumoral efficacy of BCMA CAR T cells and correlates with clinical outcome.

Identification of new markers associated with long-term efficacy in patients treated with CAR T cells is a current medical need, particularly in diseases such as multiple myeloma. In this study, we address the impact of CAR density on the functionality of BCMA CAR T cells. Functional and transcriptional studies demonstrate that CAR T cells with high expression of the CAR construct show an increased tonic signaling with up-regulation of exhaustion markers and increased in vitro cytotoxicity but a decrease in in vivo BM infiltration. Characterization of gene regulatory networks using scRNA-seq identified regulons associated to activation and exhaustion up-regulated in CARHigh T cells, providing mechanistic insights behind differential functionality of these cells. Last, we demonstrate that patients treated with CAR T cell products enriched in CARHigh T cells show a significantly worse clinical response in several hematological malignancies. In summary, our work demonstrates that CAR density plays an important role in CAR T activity with notable impact on clinical response.

Deconvolution of the hematopoietic stem cell microenvironment reveals a high degree of specialization and conservation.

Understanding the regulation of normal and malignant human hematopoiesis requires comprehensive cell atlas of the hematopoietic stem cell (HSC) regulatory microenvironment. Here, we develop a tailored bioinformatic pipeline to integrate public and proprietary single-cell RNA sequencing (scRNA-seq) datasets. As a result, we robustly identify for the first time 14 intermediate cell states and 11 stages of differentiation in the endothelial and mesenchymal BM compartments, respectively. Our data provide the most comprehensive description to date of the murine HSC-regulatory microenvironment and suggest a higher level of specialization of the cellular circuits than previously anticipated. Furthermore, this deep characterization allows inferring conserved features in human, suggesting that the layers of microenvironmental regulation of hematopoiesis may also be shared between species. Our resource and methodology is a stepping-stone toward a comprehensive cell atlas of the BM microenvironment.

The bone marrow niche regulates redox and energy balance in MLL::AF9 leukemia stem cells.

Eradicating leukemia requires a deep understanding of the interaction between leukemic cells and their protective microenvironment. The CXCL12/CXCR4 axis has been postulated as a critical pathway dictating leukemia stem cell (LSC) chemoresistance in AML due to its role in controlling cellular egress from the marrow. Nevertheless, the cellular source of CXCL12 in the acute myeloid leukemia (AML) microenvironment and the mechanism by which CXCL12 exerts its protective role in vivo remain unresolved. Here, we show that CXCL12 produced by Prx1+ mesenchymal cells but not by mature osteolineage cells provide the necessary cues for the maintenance of LSCs in the marrow of an MLL::AF9-induced AML model. Prx1+ cells promote survival of LSCs by modulating energy metabolism and the REDOX balance in LSCs. Deletion of Cxcl12 leads to the accumulation of reactive oxygen species and DNA damage in LSCs, impairing their ability to perpetuate leukemia in transplantation experiments, a defect that can be attenuated by antioxidant therapy. Importantly, our data suggest that this phenomenon appears to be conserved in human patients. Hence, we have identified Prx1+ mesenchymal cells as an integral part of the complex niche-AML metabolic intertwining, pointing towards CXCL12/CXCR4 as a target to eradicate parenchymal LSCs in AML.

In vivo CRISPR-Cas9 inhibition of hepatic LDH as treatment of primary hyperoxaluria.

Genome-editing strategies, especially CRISPR-Cas9 systems, have substantially increased the efficiency of innovative therapeutic approaches for monogenic diseases such as primary hyperoxalurias (PHs). We have previously demonstrated that inhibition of glycolate oxidase using CRISPR-Cas9 systems represents a promising therapeutic option for PH type I (PH1). Here, we extended our work evaluating the efficacy of liver-specific inhibition of lactate dehydrogenase (LDH), a key enzyme responsible for converting glyoxylate to oxalate; this strategy would not be limited to PH1, being applicable to other PH subtypes. In this work, we demonstrate a liver-specific inhibition of LDH that resulted in a drastic reduction of LDH levels in the liver of PH1 and PH3 mice after a single-dose delivery of AAV8 vectors expressing the CRISPR-Cas9 system, resulting in reduced urine oxalate levels and kidney damage without signs of toxicity. Deep sequencing analysis revealed that this approach was safe and specific, with no off-targets detected in the liver of treated animals and no on-target/off-tissue events. Altogether, our data provide evidence that in vivo genome editing using CRISPR-Cas9 systems would represent a valuable tool for improved therapeutic approaches for PH.

Conserved and distinct roles of H3K27me3 demethylases regulating flowering time in Brassica rapa.

Epigenetic regulation is necessary for optimal organism development and preservation of gene expression profiles in the cell. In plants, the trimethylation of histone H3 lysine 27 (H3K27me3) is a silencing epigenetic mark relevant for developmental transitions like flowering. The floral transition is a key agronomic trait; however, the epigenetic mechanisms of flowering time regulation in crops remain poorly understood. Here we study the Jumonji H3K27me3 demethylases BraA.REF6 and BraA.ELF6 in Brassica rapa. Phenotypic characterization of novel mutant lines and genome-wide H3K27me3 chromatin immunoprecipitation and transcriptomic analyses indicated that BraA.REF6 plays a greater role than BraA.ELF6 in fine-tuning H3K27me3 levels. In addition, we found that braA.elf6 mutants were early flowering due to high H3K27me3 levels at B. rapa homologs of the floral repressor FLC. Unlike mutations in Arabidopsis thaliana, braA.ref6 mutants were late flowering without altering the expression of B. rapa FLC genes. Remarkably, we found that BraA.REF6 regulated a number of gibberellic acid (GA) biosynthetic genes, including a homolog of GA1, and that GA-treatment complemented the late flowering mutant phenotype. This study increases our understanding of the epigenetic regulation of flowering time in B. rapa, highlighting conserved and distinct regulatory mechanisms between model and crop species.

Single-Cell RNA Sequencing Analysis Reveals a Crucial Role for CTHRC1 (Collagen Triple Helix Repeat Containing 1) Cardiac Fibroblasts After Myocardial Infarction.

BACKGROUND: Cardiac fibroblasts (CFs) have a central role in the ventricular remodeling process associated with different types of fibrosis. Recent studies have shown that fibroblasts do not respond homogeneously to heart injury. Because of the limited set of bona fide fibroblast markers, a proper characterization of fibroblast population heterogeneity in response to cardiac damage is lacking. The purpose of this study was to define CF heterogeneity during ventricular remodeling and the underlying mechanisms that regulate CF function. METHODS: Collagen1α1-GFP (green fluorescent protein)-positive CFs were characterized after myocardial infarction (MI) by single-cell and bulk RNA sequencing, assay for transposase-accessible chromatin sequencing, and functional assays. Swine and patient samples were studied using bulk RNA sequencing. RESULTS: We identified and characterized a unique CF subpopulation that emerges after MI in mice. These activated fibroblasts exhibit a clear profibrotic signature, express high levels of Cthrc1 (collagen triple helix repeat containing 1), and localize into the scar. Noncanonical transforming growth factor-ß signaling and different transcription factors including SOX9 are important regulators mediating their response to cardiac injury. Absence of CTHRC1 results in pronounced lethality attributable to ventricular rupture. A population of CFs with a similar transcriptome was identified in a swine model of MI and in heart tissue from patients with MI and dilated cardiomyopathy. CONCLUSIONS: We report CF heterogeneity and their dynamics during the course of MI and redefine the CFs that respond to cardiac injury and participate in myocardial remodeling. Our study identifies CTHRC1 as a novel regulator of the healing scar process and a target for future translational studies.

Midbrain microglia mediate a specific immunosuppressive response under inflammatory conditions.

BACKGROUND: Inflammation is a critical process for the progression of neuronal death in neurodegenerative disorders. Microglia play a central role in neuroinflammation and may affect neuron vulnerability. Next generation sequencing has shown the molecular heterogeneity of microglial cells; however, the variability in their response to pathological inputs remains unknown. METHODS: To determine the effect of an inflammatory stimulus on microglial cells, lipopolysaccharide (LPS) was administered peripherally to mice and the inflammatory status of the cortex, hippocampus, midbrain, and striatum was assessed. Microglial activation and interaction with the immune system were analyzed in single cell suspensions obtained from the different brain regions by fluorescence-activated cell sorting, next generation RNA sequencing, real-time PCR, and immunohistochemical techniques. Antigen-presenting properties of microglia were evaluated by the ability of isolated cells to induce a clonal expansion of CD4+ T cells purified from OT-II transgenic mice. RESULTS: Under steady-state conditions, the midbrain presented a high immune-alert state characterized by the presence of two unique microglial subpopulations, one expressing the major histocompatibility complex class II (MHC-II) and acting as antigen-presenting cells and another expressing the toll-like receptor 4 (TLR4), and by the presence of a higher proportion of infiltrating CD4+ T cells. This state was not detected in the cortex, hippocampus, or striatum. Systemic LPS administration induced a general increase in classic pro-inflammatory cytokines, in co-inhibitory programmed death ligand 1 (PD-L1), and in cytotoxic T lymphocyte antigen 4 (CTLA-4) receptors, as well as a decrease in infiltrating effector T cells in all brain regions. Interestingly, a specific immune-suppressive response was observed in the midbrain which was characterized by the downregulation of MHC-II microglial expression, the upregulation of the anti-inflammatory cytokines IL10 and TGFß, and the increase in infiltrating regulatory T cells. CONCLUSIONS: These data show that the midbrain presents a high immune-alert state under steady-state conditions that elicits a specific immune-suppressive response when exposed to an inflammatory stimulus. This specific inflammatory tone and response may have an impact in neuronal viability.

Prokaryotic diversity and community composition in the Salar de Uyuni, a large scale, chaotropic salt flat.

Salar de Uyuni (SdU), with a geological history that reflects 50 000 years of climate change, is the largest hypersaline salt flat on Earth and is estimated to be the biggest lithium reservoir in the world. Its salinity reaches saturation levels for NaCl, a kosmotropic salt, and high concentrations of MgCL2 and LiCl, both salts considered important chaotrophic stressors. In addition, extreme temperatures, anoxic conditions, high UV irradiance, high albedo and extremely low concentrations of phosphorous, make SdU a unique natural extreme environment in which to contrast hypotheses about limiting factors of life diversification. Geophysical studies of brines from different sampling stations show that water activity is rather constant along SdU. Geochemical measurements show significant differences in magnesium concentration, ranging from 0.2 to 2M. This work analyses the prokaryotic diversity and community structure at four SdU sampling stations, selected according to their location and ionic composition. Prokaryotic communities were composed of both Archaea (with members of the classes Halobacteria, Thermoplasmata and Nanohaloarchaea, from the Euryarchaeota and Nanohaloarcheota phyla respectively) and Bacteria (mainly belonging to Bacteroidetes and Proteobacteria phyla). The important differences in composition of microbial communities inversely correlate with Mg2+ concentration, suggesting that prokaryotic diversity at SdU is chaotropic dependent.

Transposase interaction with the ß sliding clamp: effects on insertion sequence proliferation and transposition rate.

Insertion sequences (ISs) are ubiquitous and abundant mobile genetic elements in prokaryotic genomes. ISs often encode only one protein, the transposase, which catalyzes their transposition. Recent studies have shown that transposases of many different IS families interact with the ß sliding clamp, a DNA replication factor of the host. However, it was unclear to what extent this interaction limits or favors the ability of ISs to colonize a chromosome from a phylogenetically-distant organism, or if the strength of this interaction affects the transposition rate. Here we describe the proliferation of a member of the IS1634 family in Acidiphilium over ~600 generations of cultured growth. We demonstrate that the purified transposase binds to the ß sliding clamp of Acidiphilium, Leptospirillum and E. coli. Further, we also demonstrate that the Acidiphilium IS1634 transposase binds to the archaeal sliding clamp (PCNA) from Methanosarcina, and that the transposase encoded by Methanosarcina IS1634 binds to Acidiphilium ß. Finally, we demonstrate that increasing the strength of the interaction between ß and transposase results in a higher transposition rate in vivo. Our results suggest that the interaction could determine the potential of ISs to be mobilized in bacterial populations and also their ability to proliferate within chromosomes.

Nickel-resistance determinants in Acidiphilium sp. PM identified by genome-wide functional screening.

Acidiphilium spp. are conspicuous dwellers of acidic, metal-rich environments. Indeed, they are among the most metal-resistant organisms; yet little is known about the mechanisms behind the metal tolerance in this genus. Acidiphilium sp. PM is an environmental isolate from Rio Tinto, an acidic, metal-laden river located in southwestern Spain. The characterization of its metal resistance revealed a remarkable ability to tolerate high Ni concentrations. Here we report the screening of a genomic library of Acidiphilium sp. PM to identify genes involved in Ni resistance. This approach revealed seven different genes conferring Ni resistance to E. coli, two of which form an operon encoding the ATP-dependent protease HslVU (ClpQY). This protease was found to enhance resistance to both Ni and Co in E. coli, a function not previously reported. Other Ni-resistance determinants include genes involved in lipopolysaccharide biosynthesis and the synthesis of branched amino acids. The diversity of molecular functions of the genes recovered in the screening suggests that Ni resistance in Acidiphilium sp. PM probably relies on different molecular mechanisms.

Microbial mediated formation of Fe-carbonate minerals under extreme acidic conditions.

Discovery of Fe-carbonate precipitation in Rio Tinto, a shallow river with very acidic waters, situated in Huelva, South-western Spain, adds a new dimension to our understanding of carbonate formation. Sediment samples from this low-pH system indicate that carbonates are formed in physico-chemical conditions ranging from acid to neutral pH. Evidence for microbial mediation is observed in secondary electron images (Fig. 1), which reveal rod-shaped bacteria embedded in the surface of siderite nanocrystals. The formation of carbonates in Rio Tinto is related to the microbial reduction of ferric iron coupled to the oxidation of organic compounds. Herein, we demonstrate for the first time, that Acidiphilium sp. PM, an iron-reducing bacterium isolated from Rio Tinto, mediates the precipitation of siderite (FeCO3) under acidic conditions and at a low temperature (30°C). We describe nucleation of siderite on nanoglobules in intimate association with the bacteria cell surface. This study has major implications for understanding carbonate formation on the ancient Earth or extraterrestrial planets.

Draft genome sequence of the electricigen Acidiphilium sp. strain PM (DSM 24941).

Acidiphilium sp. strain PM (DSM 24941) was isolated from Rio Tinto's acidic, heavy metal-rich waters. Voltammetry experiments revealed that this strain is capable of electricity production even under aerobic conditions. Here we report the draft genome sequence of Acidiphilium sp. PM and a preliminary genome analysis that reveals a versatile respiratory metabolism.

Extraction of extracellular polymeric substances from extreme acidic microbial biofilms.

The efficiency of five extraction methods for extracellular polymeric substances (EPS) was compared on three benthic eukaryotic biofilms isolated from an extreme acidic river, Río Tinto (SW, Spain). Three chemical methods (MilliQ water, NaCl, and ethylenediamine tetraacetic acid [EDTA]) and two physical methods (Dowex 50.8 and Crown Ether cation exchange resins) were tested. The quality and quantity of the EPS extracted from acidic biofilms varied according to which EPS extraction protocol was used. Higher amounts were obtained when NaCl and Crown Ether resins were used as extractant agents, followed by EDTA, Dowex, and MilliQ. EPS amounts varied from approximately 155 to 478 mg g(-1) of dry weight depending on the extraction method and biofilm analyzed. EPS were primarily composed of carbohydrate, heavy metals, and humic acid, plus small quantities of proteins and DNA. Neutral hexose concentrations corresponded to more than 90% of the total EPS dry weight. The proportions of each metals in the EPS extracted with EDTA are similar to the proportions present in the water from each locality where the biofilms were collected except for Al, Cu, Zn, and Pb. In this study, the extracellular matrix heavy metal sorption efficiencies of five methods for extracting EPS from eukaryotic acidic biofilms were compared.