TiO2 nanoparticles' library toxicity (UV and non-UV exposure) - High-throughput in vivo transcriptomics reveals mechanisms.

The hazards of nanomaterials/nanoparticles (NMs/NPs) are mostly assessed using individual NMs, and a more systematic approach, using many NMs, is needed to evaluate its risks in the environment. Libraries of NMs, with a range of identified different but related characters/descriptors allow the comparison of effects across many NMs. The effects of a custom designed Fe-doped TiO2 NMs library containing 11 NMs was assessed on the soil model Enchytraeus crypticus (Oligochaeta), both with and without UV (standard fluorescent) radiation. Effects were analyzed at organism (phenotypic, survival and reproduction) and gene expression level (transcriptomics, high-throughput 4x44K microarray) to understand the underlying mechanisms. A total of 48 microarrays (20 test conditions) were done plus controls (UV and non-UV). Unique mechanisms induced by TiO2 NPs exposure included the impairment in RNA processing for TiO2_10nm, or deregulated apoptosis for 2%FeTiO2_10nm. Strikingly apparent was the size dependent effects such as induction of reproductive effects via smaller TiO2 NPs (≤12 nm) - embryo interaction, while larger particles (27 nm) caused reproductive effects through different mechanisms. Also, phagocytosis was affected by 12 and 27 nm NPs, but not by ≤11 nm. The organism level study shows the integrated response, i.e. the result after a cascade of events. While uni-cell models offer key mechanistic information, we here deliver a combined biological system level (phenotype and genotype), seldom available, especially for environmental models.

A cross entropy test allows quantitative statistical comparison of t-SNE and UMAP representations.

The advent of high-dimensional single-cell data has necessitated the development of dimensionality-reduction tools. t-Distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) are the two most frequently used approaches, allowing clear visualization of complex single-cell datasets. Despite the need for quantitative comparison, t-SNE and UMAP have largely remained visualization tools due to the lack of robust statistical approaches. Here, we have derived a statistical test for evaluating the difference between dimensionality-reduced datasets using the Kolmogorov-Smirnov test on the distributions of cross entropy of single cells within each dataset. As the approach uses the inter-relationship of single cells for comparison, the resulting statistic is robust and capable of identifying true biological variation. Further, the test provides a valid distance between single-cell datasets, allowing the organization of multiple samples into a dendrogram for quantitative comparison of complex datasets. These results demonstrate the largely untapped potential of dimensionality-reduction tools for biomedical data analysis beyond visualization.

Context-dependent effects of IL-2 rewire immunity into distinct cellular circuits.

Interleukin 2 (IL-2) is a key homeostatic cytokine, with therapeutic applications in both immunogenic and tolerogenic immune modulation. Clinical use has been hampered by pleiotropic functionality and widespread receptor expression, with unexpected adverse events. Here, we developed a novel mouse strain to divert IL-2 production, allowing identification of contextual outcomes. Network analysis identified priority access for Tregs and a competitive fitness cost of IL-2 production among both Tregs and conventional CD4 T cells. CD8 T and NK cells, by contrast, exhibited a preference for autocrine IL-2 production. IL-2 sourced from dendritic cells amplified Tregs, whereas IL-2 produced by B cells induced two context-dependent circuits: dramatic expansion of CD8+ Tregs and ILC2 cells, the latter driving a downstream, IL-5-mediated, eosinophilic circuit. The source-specific effects demonstrate the contextual influence of IL-2 function and potentially explain adverse effects observed during clinical trials. Targeted IL-2 production therefore has the potential to amplify or quench particular circuits in the IL-2 network, based on clinical desirability.

Astrocyte-targeted gene delivery of interleukin 2 specifically increases brain-resident regulatory T cell numbers and protects against pathological neuroinflammation.

The ability of immune-modulating biologics to prevent and reverse pathology has transformed recent clinical practice. Full utility in the neuroinflammation space, however, requires identification of both effective targets for local immune modulation and a delivery system capable of crossing the blood-brain barrier. The recent identification and characterization of a small population of regulatory T (Treg) cells resident in the brain presents one such potential therapeutic target. Here, we identified brain interleukin 2 (IL-2) levels as a limiting factor for brain-resident Treg cells. We developed a gene-delivery approach for astrocytes, with a small-molecule on-switch to allow temporal control, and enhanced production in reactive astrocytes to spatially direct delivery to inflammatory sites. Mice with brain-specific IL-2 delivery were protected in traumatic brain injury, stroke and multiple sclerosis models, without impacting the peripheral immune system. These results validate brain-specific IL-2 gene delivery as effective protection against neuroinflammation, and provide a versatile platform for delivery of diverse biologics to neuroinflammatory patients.

AutoSpill is a principled framework that simplifies the analysis of multichromatic flow cytometry data.

Compensating in flow cytometry is an unavoidable challenge in the data analysis of fluorescence-based flow cytometry. Even the advent of spectral cytometry cannot circumvent the spillover problem, with spectral unmixing an intrinsic part of such systems. The calculation of spillover coefficients from single-color controls has remained essentially unchanged since its inception, and is increasingly limited in its ability to deal with high-parameter flow cytometry. Here, we present AutoSpill, an alternative method for calculating spillover coefficients. The approach combines automated gating of cells, calculation of an initial spillover matrix based on robust linear regression, and iterative refinement to reduce error. Moreover, autofluorescence can be compensated out, by processing it as an endogenous dye in an unstained control. AutoSpill uses single-color controls and is compatible with common flow cytometry software. AutoSpill allows simpler and more robust workflows, while reducing the magnitude of compensation errors in high-parameter flow cytometry.

Increased IL-10-producing regulatory T cells are characteristic of severe cases of COVID-19.

OBJECTIVES: The pandemic spread of the coronavirus SARS-CoV-2 is due, in part, to the immunological properties of the host-virus interaction. The clinical presentation varies from individual to individual, with asymptomatic carriers, mild-to-moderate-presenting patients and severely affected patients. Variation in immune response to SARS-CoV-2 may underlie this clinical variation. METHODS: Using a high-dimensional systems immunology platform, we have analysed the peripheral blood compartment of 6 healthy individuals, 23 mild-to-moderate and 20 severe COVID-19 patients. RESULTS: We identify distinct immunological signatures in the peripheral blood of the mild-to-moderate and severe COVID-19 patients, including T-cell lymphopenia, more consistent with peripheral hypo- than hyper-immune activation. Unique to the severe COVID-19 cases was a large increase in the proportion of IL-10-secreting regulatory T cells, a lineage known to possess anti-inflammatory properties in the lung. CONCLUSION: As IL-10-secreting regulatory T cells are known to possess anti-inflammatory properties in the lung, their proportional increase could contribute to a more severe COVID-19 phenotype. We openly provide annotated data (https://flowrepository.org/experiments/2713) with clinical correlates as a systems immunology resource for the COVID-19 research community.

Microglia Require CD4 T Cells to Complete the Fetal-to-Adult Transition.

The brain is a site of relative immune privilege. Although CD4 T cells have been reported in the central nervous system, their presence in the healthy brain remains controversial, and their function remains largely unknown. We used a combination of imaging, single cell, and surgical approaches to identify a CD69+ CD4 T cell population in both the mouse and human brain, distinct from circulating CD4 T cells. The brain-resident population was derived through in situ differentiation from activated circulatory cells and was shaped by self-antigen and the peripheral microbiome. Single-cell sequencing revealed that in the absence of murine CD4 T cells, resident microglia remained suspended between the fetal and adult states. This maturation defect resulted in excess immature neuronal synapses and behavioral abnormalities. These results illuminate a role for CD4 T cells in brain development and a potential interconnected dynamic between the evolution of the immunological and neurological systems. VIDEO ABSTRACT.

High-throughput transcriptomics: Insights into the pathways involved in (nano) nickel toxicity in a key invertebrate test species.

Nickel nanoparticles (NiNPs) have an estimated production of ca. 20 tons per year in the US. Nickel has been risk-assessed for long in Europe, but not NiNPs, hence the concern for the environment. In the present study, we focused on investigating the mechanisms of toxicity of NiNPs and the comparison to NiNO3. The high-throughput microarray for the soil ecotox model Enchytraeus crypticus (Oligochaeta) was used. To anchor gene to phenotype effect level, organisms were exposed to reproduction effect concentrations EC20 and EC50, for 3 and 7 days. Results showed commonly affected pathways between NiNPs and NiNO3, including increase in proteolysis, apoptosis and inflammatory response, and interference with the nervous system. Mechanisms unique to NiNO3 were also observed (e.g. glutathione synthesis). No specific mechanisms for NiNPs were found, which could indicate that longer exposure period (>7 days) is required to capture the peak response to NiNPs. A mechanisms scheme is assembled, showing both common and unique mechanisms to NiNO3 and NiNPs, providing an important framework for further, more targeted, studies.

Murine myeloproliferative disorder as a consequence of impaired collaboration between dendritic cells and CD4 T cells.

Dendritic cells (DCs) are a key cell type in the initiation of the adaptive immune response. Recently, an additional role for DCs in suppressing myeloproliferation was discovered. Myeloproliferative disorder (MPD) was observed in murine studies with constitutive depletion of DCs, as well as in patients with congenital deficiency in DCs caused by mutations in GATA2 or IRF8 The mechanistic link between DC deficiency and MPD was not predicted through the known biology and has remained an enigma. Prevailing models suggest numerical DC deficiency leads to MPD through compensatory myeloid differentiation. Here, we formally tested whether MPD can also arise through a loss of DC function without numerical deficiency. Using mice whose DCs are deficient in antigen presentation, we find spontaneous MPD that is characterized by splenomegaly, neutrophilia, and extramedullary hematopoiesis, despite normal numbers of DCs. Disease development was dependent on loss of the MHC class II (MHCII) antigen-presenting complex on DCs and was eliminated in mice deficient in total lymphocytes. Mice lacking MHCII and CD4 T cells did not develop disease. Thus, MPD was paradoxically contingent on the presence of CD4 T cells and on a failure of DCs to activate CD4 T cells, trapping the cells in a naive Flt3 ligand-expressing state. These results identify a novel requirement for intercellular collaboration between DCs and CD4 T cells to regulate myeloid differentiation. Our findings support a new conceptual framework of DC biology in preventing MPD in mice and humans.

NFIL3 mutations alter immune homeostasis and sensitise for arthritis pathology.

OBJECTIVES: NFIL3 is a key immunological transcription factor, with knockout mice studies identifying functional roles in multiple immune cell types. Despite the importance of NFIL3, little is known about its function in humans. METHODS: Here, we characterised a kindred of two monozygotic twin girls with juvenile idiopathic arthritis at the genetic and immunological level, using whole exome sequencing, single cell sequencing and flow cytometry. Parallel studies were performed in a mouse model. RESULTS: The patients inherited a novel p.M170I in NFIL3 from each of the parents. The mutant form of NFIL3 demonstrated reduced stability in vitro. The potential contribution of this mutation to arthritis susceptibility was demonstrated through a preclinical model, where Nfil3-deficient mice upregulated IL-1ß production, with more severe arthritis symptoms on disease induction. Single cell sequencing of patient blood quantified the transcriptional dysfunctions present across the peripheral immune system, converging on IL-1ß as a pivotal cytokine. CONCLUSIONS: NFIL3 mutation can sensitise for arthritis development, in mice and humans, and rewires the innate immune system for IL-1ß over-production.

Mechanisms of (photo)toxicity of TiO2 nanomaterials (NM103, NM104, NM105): using high-throughput gene expression in Enchytraeus crypticus.

Titanium dioxide (TiO2) based nanomaterials (NMs) are among the most produced NMs worldwide. When irradiated with light, particularly UV, TiO2 is photoactive, a property that is explored for several purposes. There are an increasing number of reports on the negative effects of photoactivated TiO2 on non-target organisms. We have here studied the effect of a suite of reference type TiO2 NMs (i.e. NM103, NM104, and NM105 and compared these to the bulk) with and without UV radiation to the oligochaete Enchytraeus crypticus. High-throughput gene expression was used to assess the molecular mechanisms, while also anchoring it to the known effects at the organism level (i.e., reproduction). Results showed that the photoactivity of TiO2 (UV exposed) played a major role in enhancing TiO2 toxicity, activating the transcription of oxidative stress, lysosome damage and apoptosis mechanisms. For non-UV activated TiO2, where toxicity at the organism level (reproduction) was lower, results showed potential for long-term effects (i.e., mutagenic and epigenetic). NM specific mechanisms were identified: NM103 affected transcription and translation, NM104_UV negatively affected the reproductive system/organs, and NM105_UV activated superoxide anion response. Results provided mechanistic information on UV-related phototoxicity of TiO2 materials and evidence for the potential long-term effects.

High-throughput gene expression in soil invertebrate embryos - Mechanisms of Cd toxicity in Enchytraeus crypticus.

Gene expression can vary with the organisms' life stage. It is known that embryos can be more sensitive to toxicant exposure, as previously demonstrated for Enchytraeus crypticus (Oligochaeta) exposed to cadmium (Cd), known to cause embryotoxicity and hatching delay. It was shown that Ca enters embryos via the L-type Ca channels in the cocoon membrane, this being affected in Cd exposed embryos (Cd-Ca competition is well-known). In the present study, the embryotoxic mechanisms of Cd were studied via high-throughput gene expression for E. crypticus. Cocoons (1-2 days old), instead of the adult organism, were exposed in Cd spiked LUFA 2.2 soil during 1 day. Results showed that Cd affected Ca homeostasis which is implicated in several other molecular processes. Several of the major modulators of Cd toxicity (e.g., impaired gene expression, cell cycle arrest, DNA and mitochondrial damage) were identified in the embryos showing its relevancy as a model in ecotoxicogenomics. The draft Adverse Outcome Pathway was improved. Previously was hypothesized that gene regulation mechanisms were activated to synthesize more Ca channel proteins - this was confirmed here. Further, novel evidences were that, besides the extracellular competition, Cd competes intracellularly which causes a reduction in Ca efflux, and potentiates Cd embryotoxicity.

Identifying conserved UV exposure genes and mechanisms.

Studies have been showing how changes in ultraviolet (UV) affect the terrestrial system, mostly focusing on higher plants and indirect effects, e.g. UV changed food quality/decomposition. Much less attention has been given to direct effect on terrestrial species, although the negative effects have been recognized for some earthworms. Further, the actual mechanisms of UV toxicity to soil invertebrates are even less understood. We here studied the effect of UV on the soil oligochaete Enchytraeus crypticus, and attempted to identify the possible mechanisms of toxicity using high-throughput gene expression. Applying a UV dose equivalent to UV during the winter months in northern Europe we observed an 80% decrease in reproduction. For these organisms, approximately 5% of the genes were differentially expressed. Among the observations was an activation of the DNA repair mechanisms, nucleotide excision repair, which correlated with survival of the organisms. An observed repressing of apoptosis seems to have deleterious effects (e.g. because it may lead to the accumulation of aberrant cells) leading to a decline in reproduction. The mechanisms activated by UV were similar to those mechanisms activated in humans, showing conservation across species.

High-throughput tool to discriminate effects of NMs (Cu-NPs, Cu-nanowires, CuNO3, and Cu salt aged): transcriptomics in Enchytraeus crypticus.

The current testing of nanomaterials (NMs) via standard toxicity tests does not cover many of the NMs specificities. One of the recommendations lays on understanding the mechanisms of action, as these can help predicting long-term effects and safe-by-design production. In the present study, we used the high-throughput gene expression tool, developed for Enchytraeus crypticus (4 × 44k Agilent microarray), to study the effects of exposure to several copper (Cu) forms. The Cu treatments included two NMs (spherical and wires) and two copper-salt treatments (CuNO3 spiked and Cu salt field historical contamination). To relate gene expression with higher effect level, testing was done with reproduction effect concentrations (EC20, EC50), using 3 and 7 days as exposure periods. Results showed that time plays a major role in the transcriptomic response, most of it occurring after 3 days. Analysis of gene expression profiles showed that Cu-salt-aged and Cu-nanowires (Nwires) differed from CuNO3 and Cu-nanoparticles (NPs). Functional analysis revealed specific mechanisms: Cu-NPs uniquely affected senescence and cuticle pattern formation, which can result from the contact of the NPs with the worms' tegument. Cu-Nwires affected reproduction via male gamete generation and hermaphrodite genitalia development. CuNO3 affected neurotransmission and locomotory behavior, both of which can be related with avoidance response. Cu salt-aged uniquely affected phagocytosis and reproductive system development (via different mechanisms than Cu-Nwires). For the first time for Cu (nano)materials, the adverse outcome pathways (AOPs) drafted here provide an overview for common and unique effects per material and linkage with apical effects.

Corrigendum: Variation-preserving normalization unveils blind spots in gene expression profiling.

This corrects the article DOI: 10.1038/srep42460.

JAK/STAT and TGF-ß activation as potential adverse outcome pathway of TiO2NPs phototoxicity in Caenorhabditis elegans.

Titanium dioxide nanoparticles (TiO2NPs) are widely used nanoparticles, whose catalytic activity is mainly due to photoactivation. In this study, the toxicity of TiO2NPs was investigated on the nematode Caenorhabditis elegans, with and without UV activation. Comparative analyses across the four treatments revealed that UV-activated TiO2NPs led to significant reproductive toxicity through oxidative stress. To understand the underlying molecular mechanism, transcriptomics and metabolomics analyses were conducted, followed by whole-genome network-based pathway analyses. Differential expression analysis from microarray data revealed only 4 DEGs by exposure to TiO2NPs alone, compared to 3,625 and 3,286 DEGs by UV alone and UV-activated TiO2NPs, respectively. Pathway analyses suggested the possible involvement of the JAK/STAT and TGF-ß pathways in the phototoxicity of TiO2NPs, which correlated with the observation of increased gene expression of those pathways. Comparative analysis of C. elegans response across UV activation and TiO2NPs exposure was performed using loss-of-function mutants of genes in these pathways. Results indicated that the JAK/STAT pathway was specific to TiO2NPs, whereas the TGF-ß pathway was specific to UV. Interestingly, crosstalk between these pathways was confirmed by further mutant analysis. We consider that these findings will contribute to understand the molecular mechanisms of toxicity of TiO2NPs in the natural environment.

Variation-preserving normalization unveils blind spots in gene expression profiling.

RNA-Seq and gene expression microarrays provide comprehensive profiles of gene activity, but lack of reproducibility has hindered their application. A key challenge in the data analysis is the normalization of gene expression levels, which is currently performed following the implicit assumption that most genes are not differentially expressed. Here, we present a mathematical approach to normalization that makes no assumption of this sort. We have found that variation in gene expression is much larger than currently believed, and that it can be measured with available assays. Our results also explain, at least partially, the reproducibility problems encountered in transcriptomics studies. We expect that this improvement in detection will help efforts to realize the full potential of gene expression profiling, especially in analyses of cellular processes involving complex modulations of gene expression.

An Integrated Data-Driven Strategy for Safe-by-Design Nanoparticles: The FP7 MODERN Project.

The development and implementation of safe-by-design strategies is key for the safe development of future generations of nanotechnology enabled products. The safety testing of the huge variety of nanomaterials that can be synthetized is unfeasible due to time and cost constraints. Computational modeling facilitates the implementation of alternative testing strategies in a time and cost effective way. The development of predictive nanotoxicology models requires the use of high quality experimental data on the structure, physicochemical properties and bioactivity of nanomaterials. The FP7 Project MODERN has developed and evaluated the main components of a computational framework for the evaluation of the environmental and health impacts of nanoparticles. This chapter describes each of the elements of the framework including aspects related to data generation, management and integration; development of nanodescriptors; establishment of nanostructure-activity relationships; identification of nanoparticle categories; hazard ranking and risk assessment.

A systems toxicology approach reveals the Wnt-MAPK crosstalk pathway mediated reproductive failure in Caenorhabditis elegans exposed to graphene oxide (GO) but not to reduced graphene oxide (rGO).

The potential hazards of graphene nanomaterials were investigated by exposing the nematode Caenorhabditis elegans to graphene oxide (GO) and reduced graphene oxide (rGO). The underlying mechanisms of the nano-bio interaction were addressed with an integrated systems toxicology approach using global transcriptomics, network-based pathway analysis, and experimental validation of the in-silico-derived hypotheses. Graphene oxide was found to reduce the worms' reproductive health to a greater degree than rGO, but it did not affect survival (24 h endpoint). Comparative analysis of GO vs. rGO effects found that the wingless-type MMTV integration site family (Wnt) pathway and the mitogen-activated protein kinase (MAPK) pathway were evoked in GO- but not in rGO-exposed worms. We therefore hypothesized that crosstalk between the Wnt and MAPK pathways is responsible for C. elegans' reproductive sensitivity to GO exposure. By targeting the individual components of the Wnt-MAPK crosstalk pathway (with qPCR gene expression and mutant reproduction analysis), we found a signaling cascade of MOM-2 → MOM-5 → MOM-4 → LIT-1 → POP-1 → EGL-5. Specifically, the activation of POP-1 (the TCF protein homolog) and subsequent repression of the Wnt/ß-catenin target gene (EGL-5), analyzed with target-gene-specific RNAi in POP-1 mutant [pop-1(q645)] worms, were the central mechanisms of reduced reproductive potential in the worms exposed to GO. Our results highlight the distinct biological and molecular mechanisms of GO and rGO exposure and the role of Wnt-MAPK pathway crosstalk in regulating GO-induced reproductive failure in in vivo systems, and they will contribute to the development of efficient and innocuous graphene applications as well to improvements in mechanism-based risk assessment.