RUNX transcription factors are essential in maintaining epididymal epithelial differentiation.

Apart from the androgen receptor, transcription factors (TFs) that are required for the development and formation of the different segments of the epididymis have remained unknown. We identified TF families expressed in the developing epididymides, of which many showed segment specificity. From these TFs, down-regulation of runt related transcription factors (RUNXs) 1 and 2 expression coincides with epithelial regression in Dicer1 cKO mice. Concomitant deletion of both Runx1 and Runx2 in a mouse epididymal epithelial cell line affected cell morphology, adhesion and mobility in vitro. Furthermore, lack of functional RUNXs severely disturbed the formation of 3D epididymal organoid-like structures. Transcriptomic analysis of the epididymal cell organoid-like structures indicated that RUNX1 and RUNX2 are involved in the regulation of MAPK signaling, NOTCH pathway activity, and EMT-related gene expression. This suggests that RUNXs are master regulators of several essential signaling pathways, and necessary for the maintenance of proper differentiation of the epididymal epithelium.

PIM kinases regulate early human Th17 cell differentiation.

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

Profiling of peripheral blood B-cell transcriptome in children who developed coeliac disease in a prospective study.

Background: In coeliac disease (CoD), the role of B-cells has mainly been considered to be production of antibodies. The functional role of B-cells has not been analysed extensively in CoD. Methods: We conducted a study to characterize gene expression in B-cells from children developing CoD early in life using samples collected before and at the diagnosis of the disease. Blood samples were collected from children at risk at 12, 18, 24 and 36 months of age. RNA from peripheral blood CD19+ cells was sequenced and differential gene expression was analysed using R package Limma. Findings: Overall, we found one gene, HNRNPL, modestly downregulated in all patients (logFC -0·7; q = 0·09), and several others downregulated in those diagnosed with CoD already by the age of 2 years. Interpretation: The data highlight the role of B-cells in CoD development. The role of HNRPL in suppressing enteroviral replication suggests that the predisposing factor for both CoD and enteroviral infections is the low level of HNRNPL expression. Funding: EU FP7 grant no. 202063, EU Regional Developmental Fund and research grant PRG712, The Academy of Finland Centre of Excellence in Molecular Systems Immunology and Physiology Research (SyMMyS) 2012-2017, grant no. 250114) and, AoF Personalized Medicine Program (grant no. 292482), AoF grants 292335, 294337, 319280, 31444, 319280, 329277, 331790) and grants from the Sigrid Jusélius Foundation (SJF).

SMG6 localizes to the chromatoid body and shapes the male germ cell transcriptome to drive spermatogenesis.

Nonsense-mediated RNA decay (NMD) is a highly conserved and selective RNA turnover pathway that depends on the endonuclease SMG6. Here, we show that SMG6 is essential for male germ cell differentiation in mice. Germ-cell conditional knockout (cKO) of Smg6 induces extensive transcriptome misregulation, including a failure to eliminate meiotically expressed transcripts in early haploid cells, and accumulation of NMD target mRNAs with long 3' untranslated regions (UTRs). Loss of SMG6 in the male germline results in complete arrest of spermatogenesis at the early haploid cell stage. We find that SMG6 is strikingly enriched in the chromatoid body (CB), a specialized cytoplasmic granule in male germ cells also harboring PIWI-interacting RNAs (piRNAs) and the piRNA-binding protein PIWIL1. This raises the possibility that SMG6 and the piRNA pathway function together, which is supported by several findings, including that Piwil1-KO mice phenocopy Smg6-cKO mice and that SMG6 and PIWIL1 co-regulate many genes in round spermatids. Together, our results demonstrate that SMG6 is an essential regulator of the male germline transcriptome, and highlight the CB as a molecular platform coordinating RNA regulatory pathways to control sperm production and fertility.

MYO10-filopodia support basement membranes at pre-invasive tumor boundaries.

Ductal carcinoma in situ (DCIS) is a pre-invasive stage of breast cancer. During invasion, the encapsulating DCIS basement membrane (BM) is compromised, and tumor cells invade the surrounding stroma. The mechanisms that regulate functional epithelial BMs in vivo are poorly understood. Myosin-X (MYO10) is a filopodia-inducing protein associated with metastasis and poor clinical outcome in invasive breast cancer (IBC). We identify elevated MYO10 expression in human DCIS and IBC, and this suggests links with disease progression. MYO10 promotes filopodia formation and cell invasion in vitro and cancer-cell dissemination from progressively invasive human DCIS xenografts. However, MYO10-depleted xenografts are more invasive. These lesions exhibit compromised BMs, poorly defined borders, and increased cancer-cell dispersal and EMT-marker-positive cells. In addition, cancer spheroids are dependent on MYO10-filopodia to generate a near-continuous extracellular matrix boundary. Thus, MYO10 is protective in early-stage breast cancer, correlating with tumor-limiting BMs, and pro-invasive at later stages, facilitating cancer-cell dissemination.

Umbilical cord blood DNA methylation in children who later develop type 1 diabetes.

AIMS/HYPOTHESIS: Distinct DNA methylation patterns have recently been observed to precede type 1 diabetes in whole blood collected from young children. Our aim was to determine whether perinatal DNA methylation is associated with later progression to type 1 diabetes. METHODS: Reduced representation bisulphite sequencing (RRBS) analysis was performed on umbilical cord blood samples collected within the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) Study. Children later diagnosed with type 1 diabetes and/or who tested positive for multiple islet autoantibodies (n = 43) were compared with control individuals (n = 79) who remained autoantibody-negative throughout the DIPP follow-up until 15 years of age. Potential confounding factors related to the pregnancy and the mother were included in the analysis. RESULTS: No differences in the umbilical cord blood methylation patterns were observed between the cases and controls at a false discovery rate <0.05. CONCLUSIONS/INTERPRETATION: Based on our results, differences between children who progress to type 1 diabetes and those who remain healthy throughout childhood are not yet present in the perinatal DNA methylome. However, we cannot exclude the possibility that such differences would be found in a larger dataset.

High intratumoral dihydrotestosterone is associated with antiandrogen resistance in VCaP prostate cancer xenografts in castrated mice.

Antiandrogen treatment resistance is a major clinical concern in castration-resistant prostate cancer (CRPC) treatment. Using xenografts of VCaP cells we showed that growth of antiandrogen resistant CRPC tumors were characterized by a higher intratumor dihydrotestosterone (DHT) concentration than that of treatment responsive tumors. Furthermore, the slow tumor growth after adrenalectomy was associated with a low intratumor DHT concentration. Reactivation of androgen signaling in enzalutamide-resistant tumors was further shown by the expression of several androgen-dependent genes. The data indicate that intratumor DHT concentration and expression of several androgen-dependent genes in CRPC lesions is an indication of enzalutamide treatment resistance and an indication of the need for further androgen blockade. The presence of an androgen synthesis, independent of CYP17A1 activity, has been shown to exist in prostate cancer cells, and thus, novel androgen synthesis inhibitors are needed for the treatment of enzalutamide-resistant CRPC tumors that do not respond to abiraterone.

A systematic comparison of FOSL1, FOSL2 and BATF-mediated transcriptional regulation during early human Th17 differentiation.

Th17 cells are essential for protection against extracellular pathogens, but their aberrant activity can cause autoimmunity. Molecular mechanisms that dictate Th17 cell-differentiation have been extensively studied using mouse models. However, species-specific differences underscore the need to validate these findings in human. Here, we characterized the human-specific roles of three AP-1 transcription factors, FOSL1, FOSL2 and BATF, during early stages of Th17 differentiation. Our results demonstrate that FOSL1 and FOSL2 co-repress Th17 fate-specification, whereas BATF promotes the Th17 lineage. Strikingly, FOSL1 was found to play different roles in human and mouse. Genome-wide binding analysis indicated that FOSL1, FOSL2 and BATF share occupancy over regulatory regions of genes involved in Th17 lineage commitment. These AP-1 factors also share their protein interacting partners, which suggests mechanisms for their functional interplay. Our study further reveals that the genomic binding sites of FOSL1, FOSL2 and BATF harbour hundreds of autoimmune disease-linked SNPs. We show that many of these SNPs alter the ability of these transcription factors to bind DNA. Our findings thus provide critical insights into AP-1-mediated regulation of human Th17-fate and associated pathologies.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age.

AIMS/HYPOTHESIS: Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies. METHODS: Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4+ T cell, CD8+ T cell and CD4-CD8- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing. RESULTS: We identified 79, 56 and 45 differentially methylated regions in CD4+ T cells, CD8+ T cells and CD4-CD8- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4+ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4+ T cells. CONCLUSIONS/INTERPRETATION: These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

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.

Evaluation of tools for identifying large copy number variations from ultra-low-coverage whole-genome sequencing data.

BACKGROUND: Detection of copy number variations (CNVs) from high-throughput next-generation whole-genome sequencing (WGS) data has become a widely used research method during the recent years. However, only a little is known about the applicability of the developed algorithms to ultra-low-coverage (0.0005-0.8×) data that is used in various research and clinical applications, such as digital karyotyping and single-cell CNV detection. RESULT: Here, the performance of six popular read-depth based CNV detection algorithms (BIC-seq2, Canvas, CNVnator, FREEC, HMMcopy, and QDNAseq) was studied using ultra-low-coverage WGS data. Real-world array- and karyotyping kit-based validation were used as a benchmark in the evaluation. Additionally, ultra-low-coverage WGS data was simulated to investigate the ability of the algorithms to identify CNVs in the sex chromosomes and the theoretical minimum coverage at which these tools can accurately function. Our results suggest that while all the methods were able to detect large CNVs, many methods were susceptible to producing false positives when smaller CNVs (< 2 Mbp) were detected. There was also significant variability in their ability to identify CNVs in the sex chromosomes. Overall, BIC-seq2 was found to be the best method in terms of statistical performance. However, its significant drawback was by far the slowest runtime among the methods (> 3 h) compared with FREEC (~ 3 min), which we considered the second-best method. CONCLUSIONS: Our comparative analysis demonstrates that CNV detection from ultra-low-coverage WGS data can be a highly accurate method for the detection of large copy number variations when their length is in millions of base pairs. These findings facilitate applications that utilize ultra-low-coverage CNV detection.

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.

Exon-level estimates improve the detection of differentially expressed genes in RNA-seq studies.

Detection of differentially expressed genes (DEGs) between different biological conditions is a key data analysis step of most RNA-sequencing studies. Conventionally, computational tools have used gene-level read counts as input to test for differential gene expression between sample condition groups. Recently, it has been suggested that statistical testing could be performed with increased power at a lower feature level prior to aggregating the results to the gene level. In this study, we systematically compared the performance of calling the DEGs when using read count data at different levels (gene, transcript, and exon) as input, in the context of two publicly available data sets. Additionally, we tested two different methods for aggregating the lower feature-level p-values to gene-level: Lancaster and empirical Brown's method. Our results show that detection of DEGs is improved compared to the conventional gene-level approach regardless of the lower feature-level used for statistical testing. The overall best balance between accuracy and false discovery rate was obtained using the exon-level approach with empirical Brown's aggregation method, which we provide as a freely available Bioconductor package EBSEA (https://bioconductor.org/packages/release/bioc/html/EBSEA.html).

Reproducibility-optimized detection of differential DNA methylation.

Aim: DNA methylation is a key epigenetic mechanism regulating gene expression. Identifying differentially methylated regions is integral to DNA methylation analysis and there is a need for robust tools reliably detecting regions with significant differences in their methylation status. Materials & methods: We present here a reproducibility-optimized test statistic (ROTS) for detection of differential DNA methylation from high-throughput sequencing or array-based data. Results: Using both simulated and real data, we demonstrate the ability of ROTS to identify differential methylation between sample groups. Conclusion: Compared with state-of-the-art methods, ROTS shows competitive sensitivity and specificity in detecting consistently differentially methylated regions.

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.

Urine microRNA Profiling Displays miR-125a Dysregulation in Children with Fragile X Syndrome.

A triplet repeat expansion leading to transcriptional silencing of the FMR1 gene results in fragile X syndrome (FXS), which is a common cause of inherited intellectual disability and autism. Phenotypic variation requires personalized treatment approaches and hampers clinical trials in FXS. We searched for microRNA (miRNA) biomarkers for FXS using deep sequencing of urine and identified 28 differentially regulated miRNAs when 219 reliably identified miRNAs were compared in dizygotic twin boys who shared the same environment, but one had an FXS full mutation, and the other carried a premutation allele. The largest increase was found in miR-125a in the FXS sample, and the miR-125a levels were increased in two independent sets of urine samples from a total of 19 FXS children. Urine miR-125a levels appeared to increase with age in control subjects, but varied widely in FXS subjects. Should the results be generalized, it could suggest that two FXS subgroups existed. Predicted gene targets of the differentially regulated miRNAs are involved in molecular pathways that regulate developmental processes, homeostasis, and neuronal function. Regulation of miR-125a has been associated with type I metabotropic glutamate receptor signaling (mGluR), which has been explored as a treatment target for FXS, reinforcing the possibility that urine miR-125a may provide a novel biomarker for FXS.

Systematic evaluation of differential splicing tools for RNA-seq studies.

Differential splicing (DS) is a post-transcriptional biological process with critical, wide-ranging effects on a plethora of cellular activities and disease processes. To date, a number of computational approaches have been developed to identify and quantify differentially spliced genes from RNA-seq data, but a comprehensive intercomparison and appraisal of these approaches is currently lacking. In this study, we systematically evaluated 10 DS analysis tools for consistency and reproducibility, precision, recall and false discovery rate, agreement upon reported differentially spliced genes and functional enrichment. The tools were selected to represent the three different methodological categories: exon-based (DEXSeq, edgeR, JunctionSeq, limma), isoform-based (cuffdiff2, DiffSplice) and event-based methods (dSpliceType, MAJIQ, rMATS, SUPPA). Overall, all the exon-based methods and two event-based methods (MAJIQ and rMATS) scored well on the selected measures. Of the 10 tools tested, the exon-based methods performed generally better than the isoform-based and event-based methods. However, overall, the different data analysis tools performed strikingly differently across different data sets or numbers of samples.

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.

Peripheral blood DNA methylation differences in twin pairs discordant for Alzheimer's disease.

BACKGROUND: Alzheimer's disease results from a neurodegenerative process that starts well before the diagnosis can be made. New prognostic or diagnostic markers enabling early intervention into the disease process would be highly valuable. Environmental and lifestyle factors largely modulate the disease risk and may influence the pathogenesis through epigenetic mechanisms, such as DNA methylation. As environmental and lifestyle factors may affect multiple tissues of the body, we hypothesized that the disease-associated DNA methylation signatures are detectable in the peripheral blood of discordant twin pairs. RESULTS: Comparison of 23 disease discordant Finnish twin pairs with reduced representation bisulfite sequencing revealed peripheral blood DNA methylation differences in 11 genomic regions with at least 15.0% median methylation difference and FDR adjusted p value ≤ 0.05. Several of the affected genes are primarily associated with neuronal functions and pathologies and do not display disease-associated differences in gene expression in blood. The DNA methylation mark in ADARB2 gene was found to be differentially methylated also in the anterior hippocampus, including entorhinal cortex, of non-twin cases and controls. Targeted bisulfite pyrosequencing of the DNA methylation mark in ADARB2 gene in 62 Finnish and Swedish twin pairs revealed that, in addition to the disease status, DNA methylation of this region is influenced by gender, age, zygosity, APOE genotype, and smoking. Further analysis of 120 Swedish twin pairs indicated that this specific DNA methylation mark is not predictive for Alzheimer's disease and becomes differentially methylated after disease onset. CONCLUSIONS: DNA methylation differences can be detected in the peripheral blood of twin pairs discordant for Alzheimer's disease. These DNA methylation signatures may have value as disease markers and provide insights into the molecular mechanisms of pathogenesis. We found no evidence that the DNA methylation marks would be associated with gene expression in blood. Further studies are needed to elucidate the potential importance of the associated genes in neuronal functions and to validate the prognostic or diagnostic value of the individual marks or marker panels.