Early-to-mid stage idiopathic Parkinson's disease shows enhanced cytotoxicity and differentiation in CD8 T-cells in females.

Neuroinflammation in the brain contributes to the pathogenesis of Parkinson's disease (PD), but the potential dysregulation of peripheral immunity has not been systematically investigated for idiopathic PD (iPD). Here we showed an elevated peripheral cytotoxic immune milieu, with more terminally-differentiated effector memory (TEMRA) CD8 T, CD8+ NKT cells and circulating cytotoxic molecules in fresh blood of patients with early-to-mid iPD, especially females, after analyzing > 700 innate and adaptive immune features. This profile, also reflected by fewer CD8+FOXP3+ T cells, was confirmed in another subcohort. Co-expression between cytotoxic molecules was selectively enhanced in CD8 TEMRA and effector memory (TEM) cells. Single-cell RNA-sequencing analysis demonstrated the accelerated differentiation within CD8 compartments, enhanced cytotoxic pathways in CD8 TEMRA and TEM cells, while CD8 central memory (TCM) and naïve cells were already more-active and transcriptionally-reprogrammed. Our work provides a comprehensive map of dysregulated peripheral immunity in iPD, proposing candidates for early diagnosis and treatments.

Akkermansia muciniphila exacerbates food allergy in fibre-deprived mice.

Alterations in the gut microbiome, including diet-driven changes, are linked to the rising prevalence of food allergy. However, little is known about how specific gut bacteria trigger the breakdown of oral tolerance. Here we show that depriving specific-pathogen-free mice of dietary fibre leads to a gut microbiota signature with increases in the mucin-degrading bacterium Akkermansia muciniphila. This signature is associated with intestinal barrier dysfunction, increased expression of type 1 and 2 cytokines and IgE-coated commensals in the colon, which result in an exacerbated allergic reaction to food allergens, ovalbumin and peanut. To demonstrate the causal role of A. muciniphila, we employed a tractable synthetic human gut microbiota in gnotobiotic mice. The presence of A. muciniphila within the microbiota, combined with fibre deprivation, resulted in stronger anti-commensal IgE coating and innate type-2 immune responses, which worsened symptoms of food allergy. Our study provides important insights into how gut microbes can regulate immune pathways of food allergy in a diet-dependent manner.

Novel, computational IgE-clustering in a population-based cross-sectional study: Mapping the allergy burden.

BACKGROUND: Even though the prevalence of allergies is increasing, population-based data are still scarce. As a read-out for chronic inflammatory information, new methods are needed to integrate individual biological measurements and lifestyle parameters to mitigate the consequences and costs of allergic burden for society. METHODS: More than 480.000 data points were collected from 1462 Luxembourg adults during the representative, cross-sectional European Health Examination Survey, spanning health and lifestyle reports. Deep IgE-profiles based on unsupervised clustering were correlated with data of the health survey. FINDINGS: 42.6% of the participants reported a physician-diagnosed allergy and 44% were found to be IgE-positive to at least one allergen or extract. The main sensitization sources were tree pollens followed by grass pollens and mites (52.4%, 51.8% and 40.3% of sensitized participants respectively), suggesting seasonal as well as perennial burden. The youngest group of participants (25-34 years old) showed the highest burden of sensitization, with 18.2% of them having IgE to 10 or more allergen groups. Unsupervised clustering revealed that the biggest cluster of 24.4% of participants was also the one with the highest medical need, marked by their multi-sensitization to respiratory sources. INTERPRETATION: Our novel approach to analyzing large biosample datasets together with health information allows the measurement of the chronic inflammatory disease burden in the general population and led to the identification of the most vulnerable groups in need of better medical care.

High-dimensional immune profiles correlate with phenotypes of peanut allergy during food-allergic reactions.

BACKGROUND: Food challenges carry a burden of safety, effort and resources. Clinical reactivity and presentation, such as thresholds and symptoms, are considered challenging to predict ex vivo. AIMS: To identify changes of peripheral immune signatures during oral food challenges (OFC) that correlate with the clinical outcome in patients with peanut allergy (PA). METHODS: Children with a positive (OFC+ , n = 16) or a negative (OFC- , n = 10) OFC-outcome were included (controls, n = 7). Single-cell mass cytometry/unsupervised analysis allowed unbiased immunophenotyping during OFC. RESULTS: Peripheral immune profiles correlated with OFC outcome. OFC+ -profiles revealed mainly decreased Th2 cells, memory Treg and activated NK cells, which had an increased homing marker expression signifying immune cell migration into effector tissues along with symptom onset. OFC- -profiles had also signs of ongoing inflammation, but with a signature of a controlled response, lacking homing marker expression and featuring a concomitant increase of Th2-shifted CD4+ T cells and Treg cells. Low versus high threshold reactivity-groups had differential frequencies of intermediate monocytes and myeloid dendritic cells at baseline. Low threshold was associated with increased CD8+ T cells and reduced memory cells (central memory [CM] CD4+ [Th2] T cells, CM CD8+ T cells, Treg). Immune signatures also discriminated patients with preferential skin versus gastrointestinal symptoms, whereby skin signs correlated with increased expression of CCR4, a molecule enabling skin trafficking, on various immune cell types. CONCLUSION: We showed that peripheral immune signatures reflected dynamics of clinical outcome during OFC with peanut. Those immune alterations hold promise as a basis for predictive OFC biomarker discovery to monitor disease outcome and therapy of PA.

Allergic airway inflammation delays glioblastoma progression and reinvigorates systemic and local immunity in mice.

BACKGROUND: Numerous patient-based studies have highlighted the protective role of immunoglobulin E-mediated allergic diseases on glioblastoma (GBM) susceptibility and prognosis. However, the mechanisms behind this observation remain elusive. Our objective was to establish a preclinical model able to recapitulate this phenomenon and investigate the role of immunity underlying such protection. METHODS: An immunocompetent mouse model of allergic airway inflammation (AAI) was initiated before intracranial implantation of mouse GBM cells (GL261). RAG1-KO mice served to assess tumor growth in a model deficient for adaptive immunity. Tumor development was monitored by MRI. Microglia were isolated for functional analyses and RNA-sequencing. Peripheral as well as tumor-associated immune cells were characterized by flow cytometry. The impact of allergy-related microglial genes on patient survival was analyzed by Cox regression using publicly available datasets. RESULTS: We found that allergy establishment in mice delayed tumor engraftment in the brain and reduced tumor growth resulting in increased mouse survival. AAI induced a transcriptional reprogramming of microglia towards a pro-inflammatory-like state, uncovering a microglia gene signature, which correlated with limited local immunosuppression in glioma patients. AAI increased effector memory T-cells in the circulation as well as tumor-infiltrating CD4+ T-cells. The survival benefit conferred by AAI was lost in mice devoid of adaptive immunity. CONCLUSION: Our results demonstrate that AAI limits both tumor take and progression in mice, providing a preclinical model to study the impact of allergy on GBM susceptibility and prognosis, respectively. We identify a potentiation of local and adaptive systemic immunity, suggesting a reciprocal crosstalk that orchestrates allergy-induced immune protection against GBM.

Sustained high expression of multiple APOBEC3 cytidine deaminases in systemic lupus erythematosus.

APOBEC3 (A3) enzymes are best known for their role as antiviral restriction factors and as mutagens in cancer. Although four of them, A3A, A3B, A3F and A3G, are induced by type-1-interferon (IFN-I), their role in inflammatory conditions is unknown. We thus investigated the expression of A3, and particularly A3A and A3B because of their ability to edit cellular DNA, in Systemic Lupus Erythematosus (SLE), a chronic inflammatory disease characterized by high IFN-α serum levels. In a cohort of 57 SLE patients, A3A and A3B, but also A3C and A3G, were upregulated ~ 10 to 15-fold (> 1000-fold for A3B) compared to healthy controls, particularly in patients with flares and elevated serum IFN-α levels. Hydroxychloroquine, corticosteroids and immunosuppressive treatment did not reverse A3 levels. The A3AΔ3B polymorphism, which potentiates A3A, was detected in 14.9% of patients and in 10% of controls, and was associated with higher A3A mRNA expression. A3A and A3B mRNA levels, but not A3C or A3G, were correlated positively with dsDNA breaks and negatively with lymphopenia. Exposure of SLE PBMCs to IFN-α in culture induced massive and sustained A3A levels by 4 h and led to massive cell death. Furthermore, the rs2853669 A > G polymorphism in the telomerase reverse transcriptase (TERT) promoter, which disrupts an Ets-TCF-binding site and influences certain cancers, was highly prevalent in SLE patients, possibly contributing to lymphopenia. Taken together, these findings suggest that high baseline A3A and A3B levels may contribute to cell frailty, lymphopenia and to the generation of neoantigens in SLE patients. Targeting A3 expression could be a strategy to reverse cell death and the generation of neoantigens.

Induction of IL-10-producing type 2 innate lymphoid cells by allergen immunotherapy is associated with clinical response.

The role of innate immune cells in allergen immunotherapy that confers immune tolerance to the sensitizing allergen is unclear. Here, we report a role of interleukin-10-producing type 2 innate lymphoid cells (IL-10+ ILC2s) in modulating grass-pollen allergy. We demonstrate that KLRG1+ but not KLRG1- ILC2 produced IL-10 upon activation with IL-33 and retinoic acid. These cells attenuated Th responses and maintained epithelial cell integrity. IL-10+ KLRG1+ ILC2s were lower in patients with grass-pollen allergy when compared to healthy subjects. In a prospective, double-blind, placebo-controlled trial, we demonstrated that the competence of ILC2 to produce IL-10 was restored in patients who received grass-pollen sublingual immunotherapy. The underpinning mechanisms were associated with the modification of retinol metabolic pathway, cytokine-cytokine receptor interaction, and JAK-STAT signaling pathways in the ILCs. Altogether, our findings underscore the contribution of IL-10+ ILC2s in the disease-modifying effect by allergen immunotherapy.

Comprehensive mapping of immune tolerance yields a regulatory TNF receptor 2 signature in a murine model of successful Fel d 1-specific immunotherapy using high-dose CpG adjuvant.

BACKGROUND: The prevalence of allergy to cat is expanding worldwide. Allergen-specific immunotherapy (AIT) has advantages over symptomatic pharmacotherapy and promises long-lasting disease control in allergic patients. However, there is still a need to improve cat AIT regarding efficacy, safety, and adherence to the treatment. Here, we aim to boost immune tolerance to the major cat allergen Fel d 1 by increasing the anti-inflammatory activity of AIT with the established immunomodulatory adjuvant CpG, but at a higher dose than previously used in AIT. METHODS: Together with CpG, we used endotoxin-free Fel d 1 as therapeutic allergen throughout the study in a BALB/c model of allergy to Fel d 1, thus mimicking the conditions of human AIT trials. Multidimensional immune phenotyping including mass cytometry (CyTOF) was applied to analyze AIT-specific immune signatures. RESULTS: We show that AIT with high-dose CpG in combination with endotoxin-free Fel d 1 reverts all major hallmarks of allergy. High-dimensional CyTOF analysis of the immune cell signatures initiating and sustaining the AIT effect indicates the simultaneous engagement of both, the pDC-Treg and B-cell axis, with the emergence of a systemic GATA3+ FoxP3hi biTreg population. The regulatory immune signature also suggests the involvement of the anti-inflammatory TNF/TNFR2 signaling cascade in NK and B cells at an early stage and in Tregs later during AIT. CONCLUSION: Our results highlight the potential of CpG adjuvant in a novel formulation to be further exploited for inducing allergen-specific tolerance in patients with cat allergy or other allergic diseases.

GigaSOM.jl: High-performance clustering and visualization of huge cytometry datasets.

BACKGROUND: The amount of data generated in large clinical and phenotyping studies that use single-cell cytometry is constantly growing. Recent technological advances allow the easy generation of data with hundreds of millions of single-cell data points with >40 parameters, originating from thousands of individual samples. The analysis of that amount of high-dimensional data becomes demanding in both hardware and software of high-performance computational resources. Current software tools often do not scale to the datasets of such size; users are thus forced to downsample the data to bearable sizes, in turn losing accuracy and ability to detect many underlying complex phenomena. RESULTS: We present GigaSOM.jl, a fast and scalable implementation of clustering and dimensionality reduction for flow and mass cytometry data. The implementation of GigaSOM.jl in the high-level and high-performance programming language Julia makes it accessible to the scientific community and allows for efficient handling and processing of datasets with billions of data points using distributed computing infrastructures. We describe the design of GigaSOM.jl, measure its performance and horizontal scaling capability, and showcase the functionality on a large dataset from a recent study. CONCLUSIONS: GigaSOM.jl facilitates the use of commonly available high-performance computing resources to process the largest available datasets within minutes, while producing results of the same quality as the current state-of-art software. Measurements indicate that the performance scales to much larger datasets. The example use on the data from a massive mouse phenotyping effort confirms the applicability of GigaSOM.jl to huge-scale studies.

Glutathione Restricts Serine Metabolism to Preserve Regulatory T Cell Function.

Regulatory T cells (Tregs) maintain immune homeostasis and prevent autoimmunity. Serine stimulates glutathione (GSH) synthesis and feeds into the one-carbon metabolic network (1CMet) essential for effector T cell (Teff) responses. However, serine's functions, linkage to GSH, and role in stress responses in Tregs are unknown. Here, we show, using mice with Treg-specific ablation of the catalytic subunit of glutamate cysteine ligase (Gclc), that GSH loss in Tregs alters serine import and synthesis and that the integrity of this feedback loop is critical for Treg suppressive capacity. Although Gclc ablation does not impair Treg differentiation, mutant mice exhibit severe autoimmunity and enhanced anti-tumor responses. Gclc-deficient Tregs show increased serine metabolism, mTOR activation, and proliferation but downregulated FoxP3. Limitation of cellular serine in vitro and in vivo restores FoxP3 expression and suppressive capacity of Gclc-deficient Tregs. Our work reveals an unexpected role for GSH in restricting serine availability to preserve Treg functionality.

High-throughput sequencing of murine immunoglobulin heavy chain repertoires using single side unique molecular identifiers on an Ion Torrent PGM.

With the advent of high-throughput sequencing (HTS), profiling immunoglobulin (IG) repertoires has become an essential part of immunological research. Advances in sequencing technology enable the IonTorrent Personal Genome Machine (PGM) to cover the full-length of IG mRNA transcripts. Nucleotide insertions and deletions (indels) are the dominant errors of the PGM sequencing platform and can critically influence IG repertoire assessments. Here, we present a PGM-tailored IG repertoire sequencing approach combining error correction through unique molecular identifier (UID) barcoding and indel detection through ImMunoGeneTics (IMGT), the most commonly used sequence alignment database for IG sequences. Using artificially falsified sequences for benchmarking, we found that IMGT's underlying algorithms efficiently detect 98% of the introduced indels. Undetected indels are either located at the end of the sequences or produce masked frameshifts with an insertion and deletion in close proximity. The complementary determining regions 3 (CDR3s) are returned correct for up to 3 insertions or 3 deletions through conservative culling. We further show, that our PGM-tailored unique molecular identifiers result in highly accurate HTS data if combined with the presented processing strategy. In this regard, considering sequences with at least two copies from datasets with UID families of minimum 3 reads result in correct sequences with over 99% confidence. Finally, we show that the protocol can readily be used to generate homogenous datasets for bulk sequencing of murine bone marrow samples. Taken together, this approach will help to establish benchtop-scale sequencing of IG heavy chain transcripts in the field of IG repertoire research.

Genome scaffolding and annotation for the pathogen vector Ixodes ricinus by ultra-long single molecule sequencing.

BACKGROUND: Global warming and other ecological changes have facilitated the expansion of Ixodes ricinus tick populations. Ixodes ricinus is the most important carrier of vector-borne pathogens in Europe, transmitting viruses, protozoa and bacteria, in particular Borrelia burgdorferi (sensu lato), the causative agent of Lyme borreliosis, the most prevalent vector-borne disease in humans in the Northern hemisphere. To faster control this disease vector, a better understanding of the I. ricinus tick is necessary. To facilitate such studies, we recently published the first reference genome of this highly prevalent pathogen vector. Here, we further extend these studies by scaffolding and annotating the first reference genome by using ultra-long sequencing reads from third generation single molecule sequencing. In addition, we present the first genome size estimation for I. ricinus ticks and the embryo-derived cell line IRE/CTVM19. RESULTS: 235,953 contigs were integrated into 204,904 scaffolds, extending the currently known genome lengths by more than 30% from 393 to 516 Mb and the N50 contig value by 87% from 1643 bp to a N50 scaffold value of 3067 bp. In addition, 25,263 sequences were annotated by comparison to the tick's North American relative Ixodes scapularis. After (conserved) hypothetical proteins, zinc finger proteins, secreted proteins and P450 coding proteins were the most prevalent protein categories annotated. Interestingly, more than 50% of the amino acid sequences matching the homology threshold had 95-100% identity to the corresponding I. scapularis gene models. The sequence information was complemented by the first genome size estimation for this species. Flow cytometry-based genome size analysis revealed a haploid genome size of 2.65Gb for I. ricinus ticks and 3.80 Gb for the cell line. CONCLUSIONS: We present a first draft sequence map of the I. ricinus genome based on a PacBio-Illumina assembly. The I. ricinus genome was shown to be 26% (500 Mb) larger than the genome of its American relative I. scapularis. Based on the genome size of 2.65 Gb we estimated that we covered about 67% of the non-repetitive sequences. Genome annotation will facilitate screening for specific molecular pathways in I. ricinus cells and provides an overview of characteristics and functions.

Focussing reduced representation CpG sequencing through judicious restriction enzyme choice.

Current restriction enzyme based reduced representation methylation analyses aim for limited, but unbiased, methylome coverage. As the current best estimate suggests that only ~20% of CpGs are dynamically regulated, we characterised the CpG and genomic context surrounding all suitable restriction enzyme sites to identify those that were located in regions rich in dynamically methylated CpGs. The restriction-site distributions for MspI, BstUI, and HhaI were non-random. CpGs in CGI and shelf+shore could be enriched, particularly in gene bodies for all genomic regions, promoters (TSS1500, TSS200), intra- (1st exon, gene body, 3'UTR, 5'UTR) and inter-genic regions. HpyCH4IV enriched CpG elements in the open sea for all genomic elements. Judicious restriction enzyme choice improves the focus of reduced representation approaches by avoiding the monopolization of read coverage by genomic regions that are irrelevant, unwanted or difficult to map, and only sequencing the most informative fraction of CpGs.

Where does transcription start? 5'-RACE adapted to next-generation sequencing.

The variability and complexity of the transcription initiation process was examined by adapting RNA ligase-mediated rapid amplification of 5' cDNA ends (5'-RACE) to Next-Generation Sequencing (NGS). We oligo-labelled 5'-m(7)G-capped mRNA from two genes, the simple mono-exonic Beta-2-Adrenoceptor (ADRB2R)and the complex multi-exonic Glucocorticoid Receptor (GR, NR3C1), and detected a variability in TSS location that has received little attention up to now. Transcription was not initiated at a fixed TSS, but from loci of 4 to 10 adjacent nucleotides. Individual TSSs had frequencies from <0.001% to 38.5% of the total gene-specific 5' m(7)G-capped transcripts. ADRB2R used a single locus consisting of 4 adjacent TSSs. Unstimulated, the GR used a total of 358 TSSs distributed throughout 38 loci, that were principally in the 5' UTRs and were spliced using established donor and acceptor sites. Complete demethylation of the epigenetically sensitive GR promoter with 5-azacytidine induced one new locus and 127 TSSs, 12 of which were unique. We induced GR transcription with dexamethasone and Interferon-γ, adding one new locus and 185 additional TSSs distributed throughout the promoter region. In-vitro the TSS microvariability regulated mRNA translation efficiency and the relative abundance of the different GRN-terminal protein isoform levels.

Integration of Ixodes ricinus genome sequencing with transcriptome and proteome annotation of the naïve midgut.

BACKGROUND: In Europe, Ixodes ricinus ticks are the most important vectors of diseases threatening humans, livestock, wildlife and companion animals. Nevertheless, genomic sequence information is missing and functional annotation of transcripts and proteins is limited. This lack of information is restricting studies of the vector and its interactions with pathogens and hosts. Here we present and integrate the first analysis of the I. ricinus genome with the transcriptome and proteome of the unfed I. ricinus midgut. METHODS: Whole genome sequencing was performed on I. ricinus ticks and the sequences were de novo assembled. In parallel, I. ricinus ticks were dissected and the midgut transcriptome sequenced. Both datasets were integrated by transcript discovery analysis to identify putative genes and genome contigs were screened for homology. An alignment-based and a motif-search-based approach were combined for the annotation of the midgut transcriptome. Additionally, midgut proteins were identified and annotated by mass spectrometry with public databases and the in-house built transcriptome database as references and results were cross-validated. RESULTS: The de novo assembly of 1 billion DNA sequences to a reference genome of 393 Mb length provides an unprecedented insight into the I. ricinus genome. A homology search revealed sequences in the assembled genome contigs homologous to 89% of the I. scapularis genome scaffolds indicating coverage of most genome regions. We identified moreover 6,415 putative genes. More than 10,000 transcripts from naïve midgut were annotated with respect of predicted function and/or cellular localization. By combining an alignment-based with a motif-search-based annotation approach, we doubled the number of annotations throughout all functional categories. In addition, 574 gel spots were significantly identified by mass spectrometry (p<0.05) and 285 distinct proteins expressed in the naïve midgut were annotated functionally and/or for cellular localization. Our systems approach reveals a midgut metabolism of the unfed tick that is prepared to sense and process an anticipated blood meal. CONCLUSIONS: This multiple-omics study vastly extends the publicly available DNA and RNA databases for I. ricinus, paving the way for further in-depth analysis of the most important European disease vector and its interactions with pathogens and hosts.

Bayesian inference of the evolution of HBV/E.

Despite its wide spread and high prevalence in sub-Saharan Africa, hepatitis B virus genotype E (HBV/E) has a surprisingly low genetic diversity, indicating an only recent emergence of this genotype in the general African population. Here, we performed extensive phylogeographic analyses, including Bayesian MCMC modeling. Our results indicate a mutation rate of 1.9 × 10(-4) substitutions per site and year (s/s/y) and confirm a recent emergence of HBV/E, most likely within the last 130 years, and only after the transatlantic slave-trade had come to an end. Our analyses suggest that HBV/E originated from the area of Nigeria, before rapidly spreading throughout sub-Saharan Africa. Interestingly, viral strains found in Haiti seem to be the result of multiple introductions only in the second half of the 20th century, corroborating an absence of a significant number of HBV/E strains in West Africa several centuries ago. Our results confirm that the hyperendemicity of HBV(E) in today's Africa is a recent phenomenon and likely the result of dramatic changes in the routes of viral transmission in a relatively recent past.