LuxRep: a technical replicate-aware method for bisulfite sequencing data analysis.

BACKGROUND: DNA methylation is commonly measured using bisulfite sequencing (BS-seq). The quality of a BS-seq library is measured by its bisulfite conversion efficiency. Libraries with low conversion rates are typically excluded from analysis resulting in reduced coverage and increased costs. RESULTS: We have developed a probabilistic method and software, LuxRep, that implements a general linear model and simultaneously accounts for technical replicates (libraries from the same biological sample) from different bisulfite-converted DNA libraries. Using simulations and actual DNA methylation data, we show that including technical replicates with low bisulfite conversion rates generates more accurate estimates of methylation levels and differentially methylated sites. Moreover, using variational inference speeds up computation time necessary for whole genome analysis. CONCLUSIONS: In this work we show that taking into account technical replicates (i.e. libraries) of BS-seq data of varying bisulfite conversion rates, with their corresponding experimental parameters, improves methylation level estimation and differential methylation detection.

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.

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.

Genome-wide profiling of interleukin-4 and STAT6 transcription factor regulation of human Th2 cell programming.

Dissecting the molecular mechanisms by which T helper (Th) cells differentiate to effector Th2 cells is important for understanding the pathogenesis of immune-mediated diseases, such as asthma and allergy. Because the STAT6 transcription factor is an upstream mediator required for interleukin-4 (IL-4)-induced Th2 cell differentiation, its targets include genes important for this process. Using primary human CD4(+) T cells, and by blocking STAT6 with RNAi, we identified a number of direct and indirect targets of STAT6 with ChIP sequencing. The integration of these data sets with detailed kinetics of IL-4-driven transcriptional changes showed that STAT6 was predominantly needed for the activation of transcription leading to the Th2 cell phenotype. This integrated genome-wide data on IL-4- and STAT6-mediated transcription provide a unique resource for studies on Th cell differentiation and, in particular, for designing interventions of human Th2 cell responses.

Genetic and epigenetic stability of human pluripotent stem cells.

Studies using high-resolution genome-wide approaches have recently reported that genomic and epigenomic alterations frequently accumulate in human pluripotent cells. Detailed characterization of these changes is crucial for understanding the impact of these alterations on self-renewal and proliferation, and particularly on the developmental and malignant potential of the cells. Such knowledge is required for the optimized and safe use of pluripotent cells for therapeutic purposes, such as regenerative cellular therapies using differentiated derivatives of pluripotent cells.In this Review, we summarize the current knowledge of the genomic and epigenomic stability of pluripotent human cells and the implications for stem cell research.

Activation of Plasmacytoid Dendritic Cells in Colon-Draining Lymph Nodes during Citrobacter rodentium Infection Involves Pathogen-Sensing and Inflammatory Pathways Distinct from Conventional Dendritic Cells.

Dendritic cells (DCs) bear the main responsibility for initiation of adaptive immune responses necessary for antimicrobial immunity. In the small intestine, afferent lymphatics convey Ags and microbial signals to mesenteric lymph nodes (LNs) to induce adaptive immune responses against microbes and food Ags derived from the small intestine. Whether the large intestine is covered by the same lymphatic system or represents its own lymphoid compartment has not been studied until very recently. We identified three small mesenteric LNs, distinct from small intestinal LNs, which drain lymph specifically from the colon, and studied DC responses to the attaching and effacing pathogen Citrobacter rodentium in these. Transcriptional profiling of conventional (CD11c(high)CD103(high)) DC and plasmacytoid (plasmacytoid DC Ag-1(high)B220(+)CD11c(int)) DC (pDC) populations during steady-state conditions revealed activity of distinct sets of genes in these two DC subsets, both in small intestinal and colon-draining LNs. C. rodentium activated DC especially in colon-draining LNs, and gene expression changed in pDC more profoundly than in conventional DC. Among the genes most upregulated in pDC were C-type lectin receptor CLEC4E, IL-1Rs (IL-1R1 and -2), proinflammatory cytokines (IL-1a and IL-6), and TLR6. Our results indicate that colon immune surveillance is distinct from that of the small intestine in terms of draining LNs, and identify pDC as active sentinels of colonic inflammation and/or microbial dysbiosis.

Placental DNA methylation marks are associated with maternal depressive symptoms during early pregnancy.

Maternal depressive symptoms during pregnancy are a significant risk factor for adverse developmental and health outcomes of the offspring. The molecular mechanisms mediating the long-term effects of this exposure are not well understood. Previous studies have found association between prenatal exposure to maternal psychological distress and placental DNA methylation of candidate genes, which can influence placental barrier function and development of the fetus. Our objective in this study was to determine epigenome wide association of maternal depressive symptoms in early pregnancy with the placental DNA methylation. For this purpose we examined DNA methylomes of 92 placental samples by using reduced representation bisulfite sequencing. The placental samples were collected after deliveries of 39 girls and 59 boys, whose mothers had Edinburgh Postnatal Depression Score ranging from 0 to 19 at gestational week 14. According to our results maternal depressive symptoms are associated with DNA methylation of 2833 CpG sites, which are particularly over-represented in genic enhancers. The genes overlapping or nearest to these sites are functionally enriched for development of neurons and show expression enrichment in several regions of developing brain. The genomic regions harboring the DNA methylation marks are enriched for single nucleotide polymorphisms associated with mental disease trait class. Potential cellular signaling cascades mediating the effects include inflammatory and hormonal pathways. As a conclusion our results suggest that maternal depressive symptoms during early pregnancy are associated with DNA methylation marks in placenta in genes, which are important for the development and long-term health of the brain. Whether similar marks can be detected in exposed children remains to be elucidated in further studies.

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.

Epigenomic Landscapes of hESC-Derived Neural Rosettes: Modeling Neural Tube Formation and Diseases.

We currently lack a comprehensive understanding of the mechanisms underlying neural tube formation and their contributions to neural tube defects (NTDs). Developing a model to study such a complex morphogenetic process, especially one that models human-specific aspects, is critical. Three-dimensional, human embryonic stem cell (hESC)-derived neural rosettes (NRs) provide a powerful resource for in vitro modeling of human neural tube formation. Epigenomic maps reveal enhancer elements unique to NRs relative to 2D systems. A master regulatory network illustrates that key NR properties are related to their epigenomic landscapes. We found that folate-associated DNA methylation changes were enriched within NR regulatory elements near genes involved in neural tube formation and metabolism. Our comprehensive regulatory maps offer insights into the mechanisms by which folate may prevent NTDs. Lastly, our distal regulatory maps provide a better understanding of the potential role of neurological-disorder-associated SNPs.

DNA methylation and Transcriptome Changes Associated with Cisplatin Resistance in Ovarian Cancer.

High-grade serous ovarian cancer is the most common ovarian cancer type. Although the combination of surgery and platinum-taxane chemotherapy provide an effective treatment, drug resistance frequently occurs leading to poor outcome. In order to clarify the molecular mechanisms of drug resistance, the DNA methylation and transcriptomic changes, associated with the development of drug resistance in high-grade serous ovarian cancer, were examined from patient derived malignant ascites cells. In parallel with large-scale transcriptome changes, cisplatin resistance was associated with loss of hypermethylation at several CpG sites primarily localized in the intergenic regions of the genome. The transcriptome and CpG methylome changes in response to cisplatin treatment of both sensitive and resistant cells were minimal, indicating the importance of post-translational mechanisms in regulating death or survival of the cells. The response of resistant cells to high concentrations of cisplatin revealed transcriptomic changes in potential key drivers of drug resistance, such as KLF4. Among the strongest changes was also induction of IL6 in resistant cells and the expression was further increased in response to cisplatin. Also, several other components of IL6 signaling were affected, further supporting previous observations on its importance in malignant transformation and development of drug resistance in ovarian cancer.

RNA Polymerase III Subunit POLR3G Regulates Specific Subsets of PolyA+ and SmallRNA Transcriptomes and Splicing in Human Pluripotent Stem Cells.

POLR3G is expressed at high levels in human pluripotent stem cells (hPSCs) and is required for maintenance of stem cell state through mechanisms not known in detail. To explore how POLR3G regulates stem cell state, we carried out deep-sequencing analysis of polyA+ and smallRNA transcriptomes present in hPSCs and regulated in POLR3G-dependent manner. Our data reveal that POLR3G regulates a specific subset of the hPSC transcriptome, including multiple transcript types, such as protein-coding genes, long intervening non-coding RNAs, microRNAs and small nucleolar RNAs, and affects RNA splicing. The primary function of POLR3G is in the maintenance rather than repression of transcription. The majority of POLR3G polyA+ transcriptome is regulated during differentiation, and the key pluripotency factors bind to the promoters of at least 30% of the POLR3G-regulated transcripts. Among the direct targets of POLR3G, POLG is potentially important in sustaining stem cell status in a POLR3G-dependent manner.

Epigenetic Silencing of the Key Antioxidant Enzyme Catalase in Karyotypically Abnormal Human Pluripotent Stem Cells.

Epigenomic regulation is likely to be important in the maintenance of genomic integrity of human pluripotent stem cells, however, the mechanisms are unknown. We explored the epigenomes and transcriptomes of human pluripotent stem cells before and after spontaneous transformation to abnormal karyotypes and in correlation to cancer cells. Our results reveal epigenetic silencing of Catalase, a key regulator of oxidative stress and DNA damage control in abnormal cells. Our findings provide novel insight into the mechanisms associated with spontaneous transformation of human pluripotent stem cells towards malignant fate. The same mechanisms may control the genomic stability of cells in somatic tissues.

Decorin in Human Colon Cancer: Localization In Vivo and Effect on Cancer Cell Behavior In Vitro.

Decorin is generally recognized as a tumor suppressing molecule. Nevertheless, although decorin has been shown to be differentially expressed in malignant tissues, it has often remained unclear whether, in addition to non-malignant stromal cells, cancer cells also express it. Here, we first used two publicly available databases to analyze the current information about decorin expression and immunoreactivity in normal and malignant human colorectal tissue samples. The analyses demonstrated that decorin expression and immunoreactivity may vary in cancer cells of human colorectal tissues. Therefore, we next examined decorin expression in normal, premalignant and malignant human colorectal tissues in more detail using both in situ hybridization and immunohistochemistry for decorin. Our results invariably demonstrate that malignant cells within human colorectal cancer tissues are devoid of both decorin mRNA and immunoreactivity. Identical results were obtained for cells of neuroendocrine tumors of human colon. Using RT-qPCR, we showed that human colon cancer cell lines are also decorin negative, in accordance with the above in vivo results. Finally, we demonstrate that decorin transduction of human colon cancer cell lines causes a significant reduction in their colony forming capability. Thus, strategies to develop decorin-based adjuvant therapies for human colorectal malignancies are highly rational.

Expression profiles of long non-coding RNAs located in autoimmune disease-associated regions reveal immune cell-type specificity.

BACKGROUND: Although genome-wide association studies (GWAS) have identified hundreds of variants associated with a risk for autoimmune and immune-related disorders (AID), our understanding of the disease mechanisms is still limited. In particular, more than 90% of the risk variants lie in non-coding regions, and almost 10% of these map to long non-coding RNA transcripts (lncRNAs). lncRNAs are known to show more cell-type specificity than protein-coding genes. METHODS: We aimed to characterize lncRNAs and protein-coding genes located in loci associated with nine AIDs which have been well-defined by Immunochip analysis and by transcriptome analysis across seven populations of peripheral blood leukocytes (granulocytes, monocytes, natural killer (NK) cells, B cells, memory T cells, naive CD4(+) and naive CD8(+) T cells) and four populations of cord blood-derived T-helper cells (precursor, primary, and polarized (Th1, Th2) T-helper cells). RESULTS: We show that lncRNAs mapping to loci shared between AID are significantly enriched in immune cell types compared to lncRNAs from the whole genome (α <0.005). We were not able to prioritize single cell types relevant for specific diseases, but we observed five different cell types enriched (α <0.005) in five AID (NK cells for inflammatory bowel disease, juvenile idiopathic arthritis, primary biliary cirrhosis, and psoriasis; memory T and CD8(+) T cells in juvenile idiopathic arthritis, primary biliary cirrhosis, psoriasis, and rheumatoid arthritis; Th0 and Th2 cells for inflammatory bowel disease, juvenile idiopathic arthritis, primary biliary cirrhosis, psoriasis, and rheumatoid arthritis). Furthermore, we show that co-expression analyses of lncRNAs and protein-coding genes can predict the signaling pathways in which these AID-associated lncRNAs are involved. CONCLUSIONS: The observed enrichment of lncRNA transcripts in AID loci implies lncRNAs play an important role in AID etiology and suggests that lncRNA genes should be studied in more detail to interpret GWAS findings correctly. The co-expression results strongly support a model in which the lncRNA and protein-coding genes function together in the same pathways.

Karyotypically abnormal human ESCs are sensitive to HDAC inhibitors and show altered regulation of genes linked to cancers and neurological diseases.

Genomic abnormalities may accumulate in human embryonic stem cells (hESCs) during in vitro maintenance. Characterization of the mechanisms enabling survival and expansion of abnormal hESCs is important due to consequences of genetic changes for the therapeutic utilization of stem cells. Furthermore, these cells provide an excellent model to study transformation in vitro. We report here that the histone deacetylase proteins, HDAC1 and HDAC2, are increased in karyotypically abnormal hESCs when compared to their normal counterparts. Importantly, similar to many cancer cell lines, we found that HDAC inhibitors repress proliferation of the karyotypically abnormal hESCs, whereas normal cells are more resistant to the treatment. The decreased proliferation correlates with downregulation of HDAC1 and HDAC2 proteins, induction of the proliferation inhibitor, cyclin-dependent kinase inhibitor 1A (CDKN1A), and altered regulation of tumor suppressor protein Retinoblastoma 1 (RB1). Through genome-wide transcriptome analysis we have identified genes with altered expression and responsiveness to HDAC inhibition in abnormal cells. Most of these genes are linked to severe developmental and neurological diseases and cancers. Our results highlight the importance of epigenetic mechanisms in the regulation of genomic stability of hESCs, and provide valuable candidates for targeted and selective growth inhibition of karyotypically abnormal cells.

High-throughput karyotyping of human pluripotent stem cells.

Genomic integrity of human pluripotent stem cell (hPSC) lines requires routine monitoring. We report here that novel karyotyping assay, utilizing bead-bound bacterial artificial chromosome probes, provides a fast and easy tool for detection of chromosomal abnormalities in hPSC lines. The analysis can be performed from low amounts of DNA isolated from whole cell pools with simple data analysis interface. The method enables routine screening of stem cell lines in a cost-efficient high-throughput manner.

Genome-wide identification of novel genes involved in early Th1 and Th2 cell differentiation.

Th cell subtypes, Th1 and Th2, are involved in the pathogenesis or progression of many immune-mediated diseases, such as type 1 diabetes and asthma, respectively. Defining the molecular networks and factors that direct Th1 and Th2 cell differentiation will help to understand the pathogenic mechanisms causing these diseases. Some of the key factors regulating this differentiation have been identified, however, they alone do not explain the process in detail. To identify novel factors directing the early differentiation, we have studied the transcriptomes of human Th1 and Th2 cells after 2, 6, and 48 h of polarization at the genome scale. Based on our current and previous studies, 288 genes or expressed sequence tags, representing approximately 1-1.5% of the human genome, are regulated in the process during the first 2 days. These transcriptional profiles revealed genes coding for components of certain pathways, such as RAS oncogene family and G protein-coupled receptor signaling, to be differentially regulated during the early Th1 and Th2 cell differentiation. Importantly, numerous novel genes with unknown functions were identified. By using short-hairpin RNA knockdown, we show that a subset of these genes is regulated by IL-4 through STAT6 signaling. Furthermore, we demonstrate that one of the IL-4 regulated genes, NDFIP2, promotes IFN-gamma production by the polarized human Th1 lymphocytes. Among the novel genes identified, there may be many factors that play a crucial role in the regulation of the differentiation process together with the previously known factors and are potential targets for developing therapeutics to modulate Th1 and Th2 responses.

Expression of human pim family genes is selectively up-regulated by cytokines promoting T helper type 1, but not T helper type 2, cell differentiation.

Cytokines are the most important inducers of T helper (Th) cell differentiation. Interleukin-12 (IL-12) and interferon-alpha (IFN-alpha) are responsible for human Th1-cell differentiation, while IL-4 is the critical cytokine promoting Th2-cell development. These two subsets of cells co-ordinate immunological responses to pathogens as well as autoimmune or allergic reactions. The pim family of proto-oncogenes encodes serine/threonine-specific kinases involved in cytokine-mediated signalling pathways in haematopoietic cells. Here we demonstrate that expression of pim-1 and pim-2 mRNAs is selectively up- or down-regulated in human cord-blood-derived CD4+ cells freshly induced to polarize towards Th1 or Th2 cells, respectively, whereas their expression is inhibited in both cell types by the immunosuppressive transforming growth factor beta (TGF-beta). Moreover, the Th1-specific cytokines IL-12 and IFN-alpha, but not the Th2-specific cytokine IL-4, transiently up-regulate pim-1 and pim-2 mRNA expression in human peripheral blood T cells and natural killer cells. In addition, the Pim-1 protein levels are strongly up-regulated by Th1-specific cytokines in all of these cell types. Taken together, our results suggest that pim genes and their protein products are involved in the early differentiation process of T helper cells.

Early target genes of IL-12 and STAT4 signaling in th cells.

IL-12 signaling through STAT4 is essential for induction of optimal levels of IFN-gamma production and commitment of Th1 cells. The molecular mechanism that controls how IL-12 and STAT4 signaling induces Th1 differentiation is poorly described. To identify the early target genes of IL-12 and STAT4 signaling, oligonucleotide arrays were used to compare the gene expression profiles of wild-type and STAT4-knockout murine Th cells during the early Th1 differentiation. According to the results, 20 genes were regulated in an IL-12- and STAT4-dependent manner. Importantly, Ifngamma was clearly the first gene induced by IL-12 in a STAT4-dependent manner. Most of the other defects in gene expression in STAT4-knockout cells were seen after 48 h of Th1 polarization. In addition to IL-12 signaling mediated by STAT4, STAT4-independent induction of a number of genes was observed immediately in response to Th1 induction. This induction was at least in part driven by IFN-gamma independently of STAT4. Importantly, addition of exogenous IFN-gamma into Th1 cell cultures of STAT4-knockout cells restored the defect in IFN-gamma production further demonstrating the critical role of IFN-gamma in early Th1 differentiation.

Kinetics and STAT4- or STAT6-mediated regulation of genes involved in lymphocyte polarization to Th1 and Th2 cells.

Many genes implicated in Th1 and Th2 differentiation have been identified in both human and mouse. However, the functional roles and hierarchy of these factors in the signaling pathways leading to either Th1 or Th2 responses are less clear. To explore at which stage of polarization the differences between Th1 and Th2 cells occur, we have studied the expression of 23 key genes implicated in the process during the first week of polarization from human precursor T helper cells using quantitative real-time reverse transcription-PCR. According to our results, 14 of the genes were clearly regulated differentially in Th1 and Th2 conditions in distinct time-dependent patterns, either during the first 2 days or after 1 week of polarization. Furthermore, 6 of these genes were identified to be targets of STAT4/6 regulation. Thus, for the first time we demonstrate expression kinetics of a number of key genes involved in Th1 and Th2 differentiation during the first week of polarization in both human and mouse. In addition, our study shows for the first time that the genes BCL-6 and TRADD are differentially regulated during the polarization of human Th1 and Th2 cells.