ROTS: An R package for reproducibility-optimized statistical testing.

Differential expression analysis is one of the most common types of analyses performed on various biological data (e.g. RNA-seq or mass spectrometry proteomics). It is the process that detects features, such as genes or proteins, showing statistically significant differences between the sample groups under comparison. A major challenge in the analysis is the choice of an appropriate test statistic, as different statistics have been shown to perform well in different datasets. To this end, the reproducibility-optimized test statistic (ROTS) adjusts a modified t-statistic according to the inherent properties of the data and provides a ranking of the features based on their statistical evidence for differential expression between two groups. ROTS has already been successfully applied in a range of different studies from transcriptomics to proteomics, showing competitive performance against other state-of-the-art methods. To promote its widespread use, we introduce here a Bioconductor R package for performing ROTS analysis conveniently on different types of omics data. To illustrate the benefits of ROTS in various applications, we present three case studies, involving proteomics and RNA-seq data from public repositories, including both bulk and single cell data. The package is freely available from Bioconductor (https://www.bioconductor.org/packages/ROTS).

Histone H3K4me3 breadth in hypoxia reveals endometrial core functions and stress adaptation linked to endometriosis.

Trimethylation of histone H3 at lysine 4 (H3K4me3) is a marker of active promoters. Broad H3K4me3 promoter domains have been associated with cell type identity, but H3K4me3 dynamics upon cellular stress have not been well characterized. We assessed this by exposing endometrial stromal cells to hypoxia, which is a major cellular stress condition. We observed that hypoxia modifies the existing H3K4me3 marks and that promoter H3K4me3 breadth rather than height correlates with transcription. Broad H3K4me3 domains mark genes for endometrial core functions and are maintained or selectively extended upon hypoxia. Hypoxic extension of H3K4me3 breadth associates with stress adaptation genes relevant for the survival of endometrial cells including transcription factor KLF4, for which we found increased protein expression in the stroma of endometriosis lesions. These results substantiate the view on broad H3K4me3 as a marker of cell identity genes and reveal participation of H3K4me3 extension in cellular stress adaptation.

Differential ATAC-seq and ChIP-seq peak detection using ROTS.

Changes in cellular chromatin states fine-tune transcriptional output and ultimately lead to phenotypic changes. Here we propose a novel application of our reproducibility-optimized test statistics (ROTS) to detect differential chromatin states (ATAC-seq) or differential chromatin modification states (ChIP-seq) between conditions. We compare the performance of ROTS to existing and widely used methods for ATAC-seq and ChIP-seq data using both synthetic and real datasets. Our results show that ROTS outperformed other commonly used methods when analyzing ATAC-seq data. ROTS also displayed the most accurate detection of small differences when modeling with synthetic data. We observed that two-step methods that require the use of a separate peak caller often more accurately called enrichment borders, whereas one-step methods without a separate peak calling step were more versatile in calling sub-peaks. The top ranked differential regions detected by the methods had marked correlation with transcriptional differences of the closest genes. Overall, our study provides evidence that ROTS is a useful addition to the available differential peak detection methods to study chromatin and performs especially well when applied to study differential chromatin states in ATAC-seq data.

Complex Interplay Between MAZR and Runx3 Regulates the Generation of Cytotoxic T Lymphocyte and Memory T Cells.

The BTB zinc finger transcription factor MAZR (also known as PATZ1) controls, partially in synergy with the transcription factor Runx3, the development of CD8 lineage T cells. Here we explored the role of MAZR as well as combined activities of MAZR/Runx3 during cytotoxic T lymphocyte (CTL) and memory CD8+ T cell differentiation. In contrast to the essential role of Runx3 for CTL effector function, the deletion of MAZR had a mild effect on the generation of CTLs in vitro. However, a transcriptome analysis demonstrated that the combined deletion of MAZR and Runx3 resulted in much more widespread downregulation of CTL signature genes compared to single Runx3 deletion, indicating that MAZR partially compensates for loss of Runx3 in CTLs. Moreover, in line with the findings made in vitro, the analysis of CTL responses to LCMV infection revealed that MAZR and Runx3 cooperatively regulate the expression of CD8α, Granzyme B and perforin in vivo. Interestingly, while memory T cell differentiation is severely impaired in Runx3-deficient mice, the deletion of MAZR leads to an enlargement of the long-lived memory subset and also partially restored the differentiation defect caused by loss of Runx3. This indicates distinct functions of MAZR and Runx3 in the generation of memory T cell subsets, which is in contrast to their cooperative roles in CTLs. Together, our study demonstrates complex interplay between MAZR and Runx3 during CTL and memory T cell differentiation, and provides further insight into the molecular mechanisms underlying the establishment of CTL and memory T cell pools.

Histone deacetylases 1 and 2 restrain CD4+ cytotoxic T lymphocyte differentiation.

Some effector CD4+ T cell subsets display cytotoxic activity, thus breaking the functional dichotomy of CD4+ helper and CD8+ cytotoxic T lymphocytes. However, molecular mechanisms regulating CD4+ cytotoxic T lymphocyte (CD4+ CTL) differentiation are poorly understood. Here we show that levels of histone deacetylases 1 and 2 (HDAC1-HDAC2) are key determinants of CD4+ CTL differentiation. Deletions of both Hdac1 and 1 Hdac2 alleles (HDAC1cKO-HDAC2HET) in CD4+ T cells induced a T helper cytotoxic program that was controlled by IFN-γ-JAK1/2-STAT1 signaling. In vitro, activated HDAC1cKO-HDAC2HET CD4+ T cells acquired cytolytic activity and displayed enrichment of gene signatures characteristic of effector CD8+ T cells and human CD4+ CTLs. In vivo, murine cytomegalovirus-infected HDAC1cKO-HDAC2HET mice displayed a stronger induction of CD4+ CTL features compared with infected WT mice. Finally, murine and human CD4+ T cells treated with short-chain fatty acids, which are commensal-produced metabolites acting as HDAC inhibitors, upregulated CTL genes. Our data demonstrate that HDAC1-HDAC2 restrain CD4+ CTL differentiation. Thus, HDAC1-HDAC2 might be targets for the therapeutic induction of CD4+ CTLs.

RepViz: a replicate-driven R tool for visualizing genomic regions.

OBJECTIVE: Visualization of sequencing data is an integral part of genomic data analysis. Although there are several tools to visualize sequencing data on genomic regions, they do not offer user-friendly ways to view simultaneously different groups of replicates. To address this need, we developed a tool that allows efficient viewing of both intra- and intergroup variation of sequencing counts on a genomic region, as well as their comparison to the output of user selected analysis methods, such as peak calling. RESULTS: We present an R package RepViz for replicate-driven visualization of genomic regions. With ChIP-seq and ATAC-seq data we demonstrate its potential to aid visual inspection involved in the evaluation of normalization, outlier behavior, detected features from differential peak calling analysis, and combined analysis of multiple data types. RepViz is readily available on Bioconductor ( https://www.bioconductor.org/packages/devel/bioc/html/RepViz.html ) and on Github ( https://github.com/elolab/RepViz ).

The Transcription Factor MAZR/PATZ1 Regulates the Development of FOXP3+ Regulatory T Cells.

Forkhead box protein P3+ (FOXP3+) regulatory T cells (Treg cells) play a key role in maintaining tolerance and immune homeostasis. Here, we report that a T cell-specific deletion of the transcription factor MAZR (also known as PATZ1) leads to an increased frequency of Treg cells, while enforced MAZR expression impairs Treg cell differentiation. Further, MAZR expression levels are progressively downregulated during thymic Treg cell development and during in-vitro-induced human Treg cell differentiation, suggesting that MAZR protein levels are critical for controlling Treg cell development. However, MAZR-deficient Treg cells show only minor transcriptional changes ex vivo, indicating that MAZR is not essential for establishing the transcriptional program of peripheral Treg cells. Finally, the loss of MAZR reduces the clinical score in dextran-sodium sulfate (DSS)-induced colitis, suggesting that MAZR activity in T cells controls the extent of intestinal inflammation. Together, these data indicate that MAZR is part of a Treg cell-intrinsic transcriptional network that modulates Treg cell development.