Human skin-resident CD8+ T cells require RUNX2 and RUNX3 for induction of cytotoxicity and expression of the integrin CD49a.

The integrin CD49a marks highly cytotoxic epidermal-tissue-resident memory (TRM) cells, but their differentiation from circulating populations remains poorly defined. We demonstrate enrichment of RUNT family transcription-factor-binding motifs in human epidermal CD8+CD103+CD49a+ TRM cells, paralleled by high RUNX2 and RUNX3 protein expression. Sequencing of paired skin and blood samples revealed clonal overlap between epidermal CD8+CD103+CD49a+ TRM cells and circulating memory CD8+CD45RA-CD62L+ T cells. In vitro stimulation of circulating CD8+CD45RA-CD62L+ T cells with IL-15 and TGF-ß induced CD49a expression and cytotoxic transcriptional profiles in a RUNX2- and RUNX3-dependent manner. We therefore identified a reservoir of circulating cells with cytotoxic TRM potential. In melanoma patients, high RUNX2, but not RUNX3, transcription correlated with a cytotoxic CD8+CD103+CD49a+ TRM cell signature and improved patient survival. Together, our results indicate that combined RUNX2 and RUNX3 activity promotes the differentiation of cytotoxic CD8+CD103+CD49a+ TRM cells, providing immunosurveillance of infected and malignant cells.

Delayed antiretroviral therapy in HIV-infected individuals leads to irreversible depletion of skin- and mucosa-resident memory T cells.

People living with HIV (PLWH) are at increased risk for developing skin and mucosal malignancies despite systemic reconstitution of CD4+ T cells upon antiretroviral therapy (ART). The underlying mechanism of chronic tissue-related immunodeficiency in HIV is unclear. We found that skin CD4+ tissue-resident memory T (Trm) cells were depleted after HIV infection and replenished only upon early ART initiation. TCR clonal analysis following early ART suggested a systemic origin for reconstituting CD4+ Trm cells. Single-cell RNA sequencing in PLWH that received late ART treatment revealed a loss of CXCR3+ Trm cells and a tolerogenic skin immune environment. Human papilloma virus-induced precancerous lesion biopsies showed reduced CXCR3+ Trm cell frequencies in the mucosa in PLWH versus HIV- individuals. These results reveal an irreversible loss of CXCR3+ Trm cells confined to skin and mucosa in PLWH who received late ART treatment, which may be a precipitating factor in the development of HPV-related cancer.

Aberrant Lipid Metabolism in Macrophages Is Associated with Granuloma Formation in Sarcoidosis.

Rationale: Chronic sarcoidosis is a complex granulomatous disease with limited treatment options that can progress over time. Understanding the molecular pathways contributing to disease would aid in new therapeutic development. Objectives: To understand whether macrophages from patients with nonresolving chronic sarcoidosis are predisposed to macrophage aggregation and granuloma formation and whether modulation of the underlying molecular pathways influence sarcoidosis granuloma formation. Methods: Macrophages were cultivated in vitro from isolated peripheral blood CD14+ monocytes and evaluated for spontaneous aggregation. Transcriptomics analyses and phenotypic and drug inhibitory experiments were performed on these monocyte-derived macrophages. Human skin biopsies from patients with sarcoidosis and a myeloid Tsc2-specific sarcoidosis mouse model were analyzed for validatory experiments. Measurements and Main Results: Monocyte-derived macrophages from patients with chronic sarcoidosis spontaneously formed extensive granulomas in vitro compared with healthy control participants. Transcriptomic analyses separated healthy and sarcoidosis macrophages and identified an enrichment in lipid metabolic processes. In vitro patient granulomas, sarcoidosis mouse model granulomas, and those directly analyzed from lesional patient skin expressed an aberrant lipid metabolism profile and contained increased neutral lipids. Conversely, a combination of statins and cholesterol-reducing agents reduced granuloma formation both in vitro and in vivo in a sarcoidosis mouse model. Conclusions: Together, our findings show that altered lipid metabolism in sarcoidosis macrophages is associated with its predisposition to granuloma formation and suggest cholesterol-reducing therapies as a treatment option in patients.

Epigenetic regulation of T cell lineages in skin and blood following hematopoietic stem cell transplantation.

Allogeneic hematopoietic stem-cell transplantation (HSCT) seeks to reconstitute the host's immune system from donor stem cells. The success of HSCT is threatened by complications including leukemia relapse or graft-versus-host-disease (GvHD). To investigate the underlying regulatory processes in central and peripheral T cell recovery, we performed sequential multi-omics analysis of T cells of the skin and blood during HSCT. We detected rapid effector T cell reconstitution, while emergence of regulatory T cells was delayed. Epigenetic and gene-regulatory programs were associated with recovering T cells and diverged greatly between skin and blood T cells. The BRG1/BRM-associated factor chromatin remodeling complex and histone deacetylases (HDACs) were epigenetic regulators involved in restoration of T cell homeostasis after transplantation. In isolated T cells of patients after HSCT, we observed class I HDAC-inhibitors to modulate their dysbalance. The present study highlights the importance of epigenetic regulation in the recovery of T cells following HSCT.

DaMold: A data-mining platform for variant annotation and visualization in molecular diagnostics research.

Next-generation sequencing (NGS) has become a powerful and efficient tool for routine mutation screening in clinical research. As each NGS test yields hundreds of variants, the current challenge is to meaningfully interpret the data and select potential candidates. Analyzing each variant while manually investigating several relevant databases to collect specific information is a cumbersome and time-consuming process, and it requires expertise and familiarity with these databases. Thus, a tool that can seamlessly annotate variants with clinically relevant databases under one common interface would be of great help for variant annotation, cross-referencing, and visualization. This tool would allow variants to be processed in an automated and high-throughput manner and facilitate the investigation of variants in several genome browsers. Several analysis tools are available for raw sequencing-read processing and variant identification, but an automated variant filtering, annotation, cross-referencing, and visualization tool is still lacking. To fulfill these requirements, we developed DaMold, a Web-based, user-friendly tool that can filter and annotate variants and can access and compile information from 37 resources. It is easy to use, provides flexible input options, and accepts variants from NGS and Sanger sequencing as well as hotspots in VCF and BED formats. DaMold is available as an online application at http://damold.platomics.com/index.html, and as a Docker container and virtual machine at https://sourceforge.net/projects/damold/.

ClinQC: a tool for quality control and cleaning of Sanger and NGS data in clinical research.

BACKGROUND: Traditional Sanger sequencing has been used as a gold standard method for genetic testing in clinic to perform single gene test, which has been a cumbersome and expensive method to test several genes in heterogeneous disease such as cancer. With the advent of Next Generation Sequencing technologies, which produce data on unprecedented speed in a cost effective manner have overcome the limitation of Sanger sequencing. Therefore, for the efficient and affordable genetic testing, Next Generation Sequencing has been used as a complementary method with Sanger sequencing for disease causing mutation identification and confirmation in clinical research. However, in order to identify the potential disease causing mutations with great sensitivity and specificity it is essential to ensure high quality sequencing data. Therefore, integrated software tools are lacking which can analyze Sanger and NGS data together and eliminate platform specific sequencing errors, low quality reads and support the analysis of several sample/patients data set in a single run. RESULTS: We have developed ClinQC, a flexible and user-friendly pipeline for format conversion, quality control, trimming and filtering of raw sequencing data generated from Sanger sequencing and three NGS sequencing platforms including Illumina, 454 and Ion Torrent. First, ClinQC convert input read files from their native formats to a common FASTQ format and remove adapters, and PCR primers. Next, it split bar-coded samples, filter duplicates, contamination and low quality sequences and generates a QC report. ClinQC output high quality reads in FASTQ format with Sanger quality encoding, which can be directly used in down-stream analysis. It can analyze hundreds of sample/patients data in a single run and generate unified output files for both Sanger and NGS sequencing data. Our tool is expected to be very useful for quality control and format conversion of Sanger and NGS data to facilitate improved downstream analysis and mutation screening. CONCLUSIONS: ClinQC is a powerful and easy to handle pipeline for quality control and trimming in clinical research. ClinQC is written in Python with multiprocessing capability, run on all major operating systems and is available at https://sourceforge.net/projects/clinqc.

Genome-wide patterns of natural variation reveal strong selective sweeps and ongoing genomic conflict in Drosophila mauritiana.

Although it is well understood that selection shapes the polymorphism pattern in Drosophila, signatures of classic selective sweeps are scarce. Here, we focus on Drosophila mauritiana, an island endemic, which is closely related to Drosophila melanogaster. Based on a new, annotated genome sequence, we characterized the genome-wide polymorphism by sequencing pooled individuals (Pool-seq). We show that the interplay between selection and recombination results in a genome-wide polymorphism pattern characteristic for D. mauritiana. Two large genomic regions (>500 kb) showed the signature of almost complete selective sweeps. We propose that the absence of population structure and limited geographic distribution could explain why such pronounced sweep patterns are restricted to D. mauritiana. Further evidence for strong adaptive evolution was detected for several nucleoporin genes, some of which were not previously identified as genes involved in genomic conflict. Since this adaptive evolution is continuing after the split of D. mauritiana and Drosophila simulans, we conclude that genomic conflict is not restricted to short episodes, but rather an ongoing process in Drosophila.

Experimental evolution reveals habitat-specific fitness dynamics among Wolbachia clades in Drosophila melanogaster.

The diversity and infection dynamics of the endosymbiont Wolbachia can be influenced by many factors, such as transmission rate, cytoplasmic incompatibility, environment, selection and genetic drift. The interplay of these factors in natural populations can result in heterogeneous infection patterns with substantial differences between populations and strains. The causes of these heterogeneities are not yet understood, partly due to the complexity of natural environments. We present experimental evolution as a new approach to study Wolbachia infection dynamics in replicate populations exposed to a controlled environment. A natural Drosophila melanogaster population infected with strains of Wolbachia belonging to different clades evolved in two laboratory environments (hot and cold) for 1.5 years. In both treatments, the rate of Wolbachia infection increased until fixation. In the hot environment, the relative frequency of different Wolbachia clades remained stable over 37 generations. In the cold environment, however, we observed marked changes in the composition of the Wolbachia population: within 15 generations, one Wolbachia clade increased more than 50% in frequency, whereas the other two clades decreased in frequency, resulting in the loss of one clade. The frequency change was highly reproducible not only among replicates, but also when flies that evolved for 42 generations in the hot environment were transferred to the cold environment. These results document how environmental factors can affect the composition of Wolbachia in D. melanogaster. The high reproducibility of the pattern suggests that experimental evolution studies can efficiently determine the functional basis of habitat-specific fitness among Wolbachia strains.

Detecting selective sweeps from pooled next-generation sequencing samples.

Due to its cost effectiveness, next-generation sequencing of pools of individuals (Pool-Seq) is becoming a popular strategy for characterizing variation in population samples. Because Pool-Seq provides genome-wide SNP frequency data, it is possible to use them for demographic inference and/or the identification of selective sweeps. Here, we introduce a statistical method that is designed to detect selective sweeps from pooled data by accounting for statistical challenges associated with Pool-Seq, namely sequencing errors and random sampling among chromosomes. This allows for an efficient use of the information: all base calls are included in the analysis, but the higher credibility of regions with higher coverage and base calls with better quality scores is accounted for. Computer simulations show that our method efficiently detects sweeps even at very low coverage (0.5× per chromosome). Indeed, the power of detecting sweeps is similar to what we could expect from sequences of individual chromosomes. Since the inference of selective sweeps is based on the allele frequency spectrum (AFS), we also provide a method to accurately estimate the AFS provided that the quality scores for the sequence reads are reliable. Applying our approach to Pool-Seq data from Drosophila melanogaster, we identify several selective sweep signatures on chromosome X that include some previously well-characterized sweeps like the wapl region.

The single-cell transcriptional landscape of innate and adaptive lymphocytes in pediatric-onset colitis.

Innate lymphoid cells (ILCs) are considered innate counterparts of adaptive T cells; however, their common and unique transcriptional signatures in pediatric inflammatory bowel disease (pIBD) are largely unknown. Here, we report a dysregulated colonic ILC composition in pIBD colitis that correlates with inflammatory activity, including accumulation of naive-like CD45RA+CD62L- ILCs. Weighted gene co-expression network analysis (WGCNA) reveals modules of genes that are shared or unique across innate and adaptive lymphocytes. Shared modules include genes associated with activation/tissue residency, naivety/quiescence, and antigen presentation. Lastly, nearest-neighbor-based analysis facilitates the identification of "most inflamed" and "least inflamed" lymphocytes in pIBD colon with unique transcriptional signatures. Our study reveals shared and unique transcriptional signatures of colonic ILCs and T cells in pIBD. We also provide insight into the transcriptional regulation of colonic inflammation, deepening our understanding of the potential mechanisms involved in pIBD.

Modular chimeric cytokine receptors with leucine zippers enhance the antitumour activity of CAR T cells via JAK/STAT signalling.

The limited availability of cytokines in solid tumours hinders maintenance of the antitumour activity of chimeric antigen receptor (CAR) T cells. Cytokine receptor signalling pathways in CAR T cells can be activated by transgenic expression or injection of cytokines in the tumour, or by engineering the activation of cognate cytokine receptors. However, these strategies are constrained by toxicity arising from the activation of bystander cells, by the suboptimal biodistribution of the cytokines and by downregulation of the cognate receptor. Here we show that replacement of the extracellular domains of heterodimeric cytokine receptors in T cells with two leucine zipper motifs provides optimal Janus kinase/signal transducer and activator of transcription signalling. Such chimeric cytokine receptors, which can be generated for common γ-chain receptors, interleukin-10 and -12 receptors, enabled T cells to survive cytokine starvation without induction of autonomous cell growth, and augmented the effector function of CAR T cells in vitro in the setting of chronic antigen exposure and in human tumour xenografts in mice. As a modular design, leucine zippers can be used to generate constitutively active cytokine receptors in effector immune cells.

PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq).

SUMMARY: Sequencing pooled DNA samples (Pool-Seq) is the most cost-effective approach for the genome-wide comparison of population samples. Here, we introduce PoPoolation2, the first software tool specifically designed for the comparison of populations with Pool-Seq data. PoPoolation2 implements a range of commonly used measures of differentiation (F(ST), Fisher's exact test and Cochran-Mantel-Haenszel test) that can be applied on different scales (windows, genes, exons, SNPs). The result may be visualized with the widely used Integrated Genomics Viewer. AVAILABILITY AND IMPLEMENTATION: PoPoolation2 is implemented in Perl and R. It is freely available on http://code.google.com/p/popoolation2/ CONTACT: christian.schloetterer@vetmeduni.ac.at SUPPLEMENTARY INFORMATION: Manual: http://code.google.com/p/popoolation2/wiki/Manual Test data and tutorial: http://code.google.com/p/popoolation2/wiki/Tutorial Validation: http://code.google.com/p/popoolation2/wiki/Validation.

Processing human skin samples for single-cell assays.

Characterizing resident immune cells in human skin using single-cell assays provides insight into their role in infections, inflammation, and cancer. We describe an optimized protocol to rapidly isolate viable cells from 6-mm skin punch-biopsies. We provide an example in which we coupled single-cell RNA sequencing (scRNA-seq) with single-cell T-cell receptor sequencing (scTCR-seq) of skin and blood cells to study transcriptional profiles and clonotypes of skin resident and peripheral circulating, memory, and effector T cells. This is an improved protocol based on Saluzzo et al. (2021). For complete details on the use and execution of this protocol, please refer to Saluzzo et al. (2021).

CD45RA+CD62L- ILCs in human tissues represent a quiescent local reservoir for the generation of differentiated ILCs.

Innate lymphoid cells (ILCs) are highly plastic and predominantly mucosal tissue-resident cells that contribute to both homeostasis and inflammation depending on the microenvironment. The discovery of naïve-like ILCs suggests an ILC differentiation process that is akin to naïve T cell differentiation. Delineating the mechanisms that underlie ILC differentiation in tissues is crucial for understanding ILC biology in health and disease. Here, we showed that tonsillar ILCs expressing CD45RA lacked proliferative activity, indicative of cellular quiescence. CD62L distinguished two subsets of CD45RA+ ILCs. CD45RA+CD62L+ ILCs (CD62L+ ILCs) resembled circulating naïve ILCs because they lacked the transcriptional, metabolic, epigenetic, and cytokine production signatures of differentiated ILCs. CD45RA+CD62L- ILCs (CD62L- ILCs) were epigenetically similar to CD62L+ ILCs but showed a transcriptional, metabolic, and cytokine production signature that was more akin to differentiated ILCs. CD62L+ and CD62L- ILCs contained uni- and multipotent precursors of ILC1s/NK cells and ILC3s. Differentiation of CD62L+ and CD62L- ILCs led to metabolic reprogramming including up-regulation of genes associated with glycolysis, which was needed for their effector functions after differentiation. CD62L- ILCs with preferential differentiation capacity toward IL-22-producing ILC3s accumulated in the inflamed mucosa of patients with inflammatory bowel disease. These data suggested distinct differentiation potential of CD62L+ and CD62L- ILCs between tissue microenvironments and identified that manipulation of these cells is a possible approach to restore tissue-immune homeostasis.

Disturbances in microbial skin recolonization and cutaneous immune response following allogeneic stem cell transfer.

The composition of the gut microbiome influences the clinical course after allogeneic hematopoietic stem cell transplantation (HSCT), but little is known about the relevance of skin microorganisms. In a single-center, observational study, we recruited a cohort of 50 patients before undergoing conditioning treatment and took both stool and skin samples up to one year after HSCT. We could confirm intestinal dysbiosis following HSCT and report that the skin microbiome is likewise perturbed in HSCT-recipients. Overall bacterial colonization of the skin was decreased after conditioning. Particularly patients that developed acute skin graft-versus-host disease (aGVHD) presented with an overabundance of Staphylococcus spp. In addition, a loss in alpha diversity was indicative of aGVHD development already before disease onset and correlated with disease severity. Further, co-localization of CD45+ leukocytes and staphylococci was observed in the skin of aGVHD patients even before disease development and paralleled with upregulated genes required for antigen-presentation in mononuclear phagocytes. Overall, our data reveal disturbances of the skin microbiome as well as cutaneous immune response in HSCT recipients with changes associated with cutaneous aGVHD.

PoPoolation DB: a user-friendly web-based database for the retrieval of natural polymorphisms in Drosophila.

BACKGROUND: The enormous potential of natural variation for the functional characterization of genes has been neglected for a long time. Only since recently, functional geneticists are starting to account for natural variation in their analyses. With the new sequencing technologies it has become feasible to collect sequence information for multiple individuals on a genomic scale. In particular sequencing pooled DNA samples has been shown to provide a cost-effective approach for characterizing variation in natural populations. While a range of software tools have been developed for mapping these reads onto a reference genome and extracting SNPs, linking this information to population genetic estimators and functional information still poses a major challenge to many researchers. RESULTS: We developed PoPoolation DB a user-friendly integrated database. Popoolation DB links variation in natural populations with functional information, allowing a wide range of researchers to take advantage of population genetic data. PoPoolation DB provides the user with population genetic parameters (Watterson's θ or Tajima's π), Tajima's D, SNPs, allele frequencies and indels in regions of interest. The database can be queried by gene name, chromosomal position, or a user-provided query sequence or GTF file. We anticipate that PoPoolation DB will be a highly versatile tool for functional geneticists as well as evolutionary biologists. CONCLUSIONS: PoPoolation DB, available at http://www.popoolation.at/pgt, provides an integrated platform for researchers to investigate natural polymorphism and associated functional annotations from UCSC and Flybase genome browsers, population genetic estimators and RNA-seq information.

Human resident memory T cells exit the skin and mediate systemic Th2-driven inflammation.

Emigration of tissue-resident memory T cells (TRMs) was recently introduced in mouse models and may drive systemic inflammation. Skin TRMs of patients undergoing allogeneic hematopoietic stem cell transplantation (HSCT) can coexist beside donor T cells, offering a unique human model system to study T cell migration. By genotyping, mathematical modeling, single-cell transcriptomics, and functional analysis of patient blood and skin T cells, we detected a small consistent population of circulating skin-derived T cells with a TRM phenotype (cTRMs) in the blood and unveil their skin origin and striking resemblance to skin TRMs. Blood from patients with active graft-versus-host disease (GVHD) contains elevated numbers of host cTRMs producing pro-inflammatory Th2/Th17 cytokines and mediating keratinocyte damage. Expression of gut-homing receptors and the occurrence of cTRMs in gastrointestinal GVHD lesions emphasize their potential to reseed and propagate inflammation in distant organs. Collectively, we describe a distinct circulating T cell population mirroring skin inflammation, which could serve as a biomarker or therapeutic target in GVHD.

The transcription factor Bcl11b promotes both canonical and adaptive NK cell differentiation.

Epigenetic landscapes can provide insight into regulation of gene expression and cellular diversity. Here, we examined the transcriptional and epigenetic profiles of seven human blood natural killer (NK) cell populations, including adaptive NK cells. The BCL11B gene, encoding a transcription factor (TF) essential for T cell development and function, was the most extensively regulated, with expression increasing throughout NK cell differentiation. Several Bcl11b-regulated genes associated with T cell signaling were specifically expressed in adaptive NK cell subsets. Regulatory networks revealed reciprocal regulation at distinct stages of NK cell differentiation, with Bcl11b repressing RUNX2 and ZBTB16 in canonical and adaptive NK cells, respectively. A critical role for Bcl11b in driving NK cell differentiation was corroborated in BCL11B-mutated patients and by ectopic Bcl11b expression. Moreover, Bcl11b was required for adaptive NK cell responses in a murine cytomegalovirus model, supporting expansion of these cells. Together, we define the TF regulatory circuitry of human NK cells and uncover a critical role for Bcl11b in promoting NK cell differentiation and function.

Contrasting seasonal niche separation between rare and abundant taxa conceals the extent of protist diversity.

With the advent of molecular methods, it became clear that microbial biodiversity had been vastly underestimated. Since then, species abundance patterns were determined for several environments, but temporal changes in species composition were not studied to the same level of resolution. Using massively parallel sequencing on the 454 GS FLX platform we identified a highly dynamic turnover of the seasonal abundance of protists in the Austrian lake Fuschlsee. We show that seasonal abundance patterns of protists closely match their biogeographic distribution. The stable predominance of few highly abundant taxa, which previously led to the suggestion of a low global protist species richness, is contrasted by a highly dynamic turnover of rare species. We suggest that differential seasonality of rare and abundant protist taxa explains the--so far--conflicting evidence in the 'everything is everywhere' dispute. Consequently temporal sampling is basic for adequate diversity and species richness estimates.

Diversity in a hidden world: potential and limitation of next-generation sequencing for surveys of molecular diversity of eukaryotic microorganisms.

With the delivery of millions of sequence reads in a single experiment, next-generation sequencing (NGS) is currently revolutionizing surveys of microorganism diversity. In particular, when applied to Eukaryotes, we are still lacking a rigorous comparison of morphological and NGS-based diversity estimates. In this report, we studied the diversity and the seasonal community turnover of alveolates (Ciliophora and Dinophyceae) in an oligotrophic freshwater lake by SSU amplicon sequencing with NGS as well as by classical morphological analysis. We complemented the morphological analysis by single-cell PCR followed by Sanger sequencing to provide an unambiguous link to the NGS data. We show that NGS and morphological analyses generally capture frequency shifts of abundant taxa over our seasonal samples. The observed incongruencies are probably largely due to rDNA copy number variation among taxa and heterogeneity in the efficiency of cell lysis. Overall, NGS-based amplicon sequencing was superior in detecting rare species. We propose that in the absence of other nuclear markers less susceptible to copy number variation, rDNA-based diversity studies need to be adjusted for confounding effects of copy number variation.