Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.

A missense mutation in SNRPE linked to non-syndromal microcephaly interferes with U snRNP assembly and pre-mRNA splicing.

Malfunction of pre-mRNA processing factors are linked to several human diseases including cancer and neurodegeneration. Here we report the identification of a de novo heterozygous missense mutation in the SNRPE gene (c.65T>C (p.Phe22Ser)) in a patient with non-syndromal primary (congenital) microcephaly and intellectual disability. SNRPE encodes SmE, a basal component of pre-mRNA processing U snRNPs. We show that the microcephaly-linked SmE variant is unable to interact with the SMN complex and as a consequence fails to assemble into U snRNPs. This results in widespread mRNA splicing alterations in fibroblast cells derived from this patient. Similar alterations were observed in HEK293 cells upon SmE depletion that could be rescued by the expression of wild type but not mutant SmE. Importantly, the depletion of SmE in zebrafish causes aberrant mRNA splicing alterations and reduced brain size, reminiscent of the patient microcephaly phenotype. We identify the EMX2 mRNA, which encodes a protein required for proper brain development, as a major mis-spliced down stream target. Together, our study links defects in the SNRPE gene to microcephaly and suggests that alterations of cellular splicing of specific mRNAs such as EMX2 results in the neurological phenotype of the disease.

Combining transcription factor binding affinities with open-chromatin data for accurate gene expression prediction.

The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid costly TF ChIP-seq assays. Thus, it is important to develop computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices. TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq. Additionally, Histone-Marks (HMs) can be used to identify candidate TF binding sites. TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength. Using machine learning, we find low affinity binding sites to improve our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites. Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance. In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq data sets. Finally, these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

Discovering the mechanisms underlying serotonin (5-HT)2A and 5-HT2C receptor regulation following nicotine withdrawal in rats.

We have previously demonstrated that nicotine withdrawal produces depression-like behavior and that serotonin (5-HT)2A/2C receptor ligands modulate that mood-like state. In the present study we aimed to identify the mechanisms (changes in radioligand binding, transcription or RNA-editing) related to such a behavioral outcome. Rats received vehicle or nicotine (0.4 mg/kg, s.c.) for 5 days in home cages. Brain 5-HT2A/2C receptors were analyzed on day 3 of nicotine withdrawal. Nicotine withdrawal increased [(3)H]ketanserin binding to 5-HT2A receptors in the ventral tegmental area and ventral dentate gyrus, yet decreased binding in the nucleus accumbens shell. Reduction in [(3)H]mesulergine binding to 5-HT2C receptors was seen in the ventral dentate gyrus. Profound decrease in the 5-HT2A receptor transcript level was noted in the hippocampus and ventral tegmental area. Out of five 5-HT2C receptor mRNA editing sites, deep sequencing data showed a reduction in editing at the E site and a trend toward reduction at the C site in the hippocampus. In the ventral tegmental area, a reduction for the frequency of CD 5-HT2C receptor transcript was seen. These results show that the reduction in the 5-HT2A receptor transcript level may be an auto-regulatory response to the increased receptor density in the hippocampus and ventral tegmental area during nicotine withdrawal, while decreased 5-HT2C receptor mRNA editing may explain the reduction in receptor labeling in the hippocampus. Serotonin (5-HT)2A/2C receptor ligands alleviate depression-like state in nicotine-withdrawn rats. Here, we show that the reduction in 5-HT2A receptor transcript level may be an auto-regulatory response to the increased receptor number in the hippocampus and ventral tegmental area during nicotine withdrawal, while attenuated 5-HT2C receptor mRNA editing in the hippocampus might explain reduced inverse agonist binding to 5-HT2C receptor and suggest a shift toward a population of more active receptors. 5-HT, serotonin; 5-HT2A R, 5-HT2A receptor; 5-HT2C R, 5-HT2C receptor.

Exome sequencing helped the fine diagnosis of two siblings afflicted with atypical Timothy syndrome (TS2).

BACKGROUND: Long-QT syndrome (LQTS) causes a prolongation of the QT-interval in the ECG leading to life threatening tachyarrhythmia and ventricular fibrillation. One atypical form of LQTS, Timothy syndrome (TS), is associated with syndactyly, immune deficiency, cognitive and neurological abnormalities as well as distinct cranio-facial abnormalities. CASE PRESENTATION: On a family with both children diagnosed with clinical LQTS, we performed whole exome sequencing to comprehensively screen for causative mutations after a targeted candidate gene panel screen for Long-QT syndrome target genes failed to identify any underlying genetic defect. Using exome sequencing, we identified in both affected children, a p.402G > S mutation in exon 8 of the CACNA1C gene, a voltage-dependent Ca2+ channel. The mutation was inherited from their father, a mosaic mutation carrier. Based on this molecular finding and further more careful clinical examination, we refined the diagnosis to be Timothy syndrome (TS2) and thereby were able to present new therapeutic approaches. CONCLUSIONS: Our study highlights the difficulties in accurate diagnosis of patients with rare diseases, especially those with atypical clinical manifestation. Such challenge could be addressed with the help of comprehensive and unbiased mutation screening, such as exome sequencing.

Ethanol deprivation and central 5-HT deficiency differentially affect the mRNA editing of the 5-HT2C receptor in the mouse brain.

BACKGROUND: Serotonin (5-HT) 5-HT2C receptor mRNA editing (at five sites, A-E), implicated in neuropsychiatric disorders, including clinical depression, remains unexplored during alcohol abstinence-often accompanied by depressive symptoms. METHODS: We used deep sequencing to investigate 5-HT2C receptor editing in mice during early ethanol deprivation following prolonged alcohol exposure and mice lacking tryptophan hydroxylase (TPH)2, a key enzyme in central 5-HT production. We also examined Tph2 expression in ethanol-deprived animals using quantitative real-time PCR (qPCR). RESULTS: Cessation from chronic 10% ethanol exposure in a two-bottle choice paradigm enhanced immobility time and decreased latency in the forced swim test (FST), indicating a depression-like phenotype. In the hippocampus, ethanol-deprived "high ethanol-drinking" mice displayed reduced Tph2 expression, elevated 5-HT2C receptor editing efficiency, and decreased frequency of the D mRNA variant, encoding the less-edited INV protein isoform. Tph2-/- mice showed attenuated receptor editing in the hippocampus and elevated frequency of non-edited None and D variants. In the prefrontal cortex, Tph2 deficiency increased receptor mRNA editing at site D and reduced the frequency of AB transcript, predicting a reduction in the corresponding partially edited VNI isoform. CONCLUSIONS: Our findings reveal differential effects of 5-HT depletion and ethanol cessation on 5-HT2C receptor editing. Central 5-HT depletion attenuated editing in the prefrontal cortex and the hippocampus, whereas ethanol deprivation, coinciding with reduced Tph2 expression in the hippocampus, enhanced receptor editing efficiency specifically in this brain region. This study highlights the interplay between 5-HT synthesis, ethanol cessation, and 5-HT2C receptor editing, providing potential mechanism underlying increased ethanol consumption and deprivation.

ACCUSA--accurate SNP calling on draft genomes.

SUMMARY: Next generation sequencing technologies facilitate genome-wide analysis of several biological processes. We are interested in whole-genome genotyping. To our knowledge, none of the existing single nucleotide polymorphism (SNP) callers consider the quality of the reference genome, which is not necessary for high-quality assemblies of well-studied model organisms. However, most genome projects will remain in draft status with little to no genome assembly improvement due to time and financial constraints. Here, we present a simple yet elegant solution ('ACCUSA') that considers both the read qualities as well as the reference genome's quality using a Bayesian framework. We demonstrate that ACCUSA is as good as the current SNP calling software in detecting true SNPs. More importantly, ACCUSA does not call spurious SNPs, which originate from a poor reference sequence. AVAILABILITY: ACCUSA is available free of charge to academic users and may be obtained from ftp://bbc.mdc-berlin.de/software. ACCUSA is programmed in JAVA 6 and runs on any platform with JAVA support. CONTACT: christoph.dieterich@mdc-berlin.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

3PD: Rapid design of optimal primers for chromosome conformation capture assays.

BACKGROUND: Higher eukaryotes control the expression of their genes by mechanisms that we are just beginning to understand. A complex layer of control is the dynamic spatial organization of the nucleus. RESULTS: We present a bioinformatics solution (3PD) to support the experimentalist in detecting long-ranging intra or inter chromosomal contacts by Chromosome conformation capture (3C) assays. 3C assays take a snapshot of chromosomal contacts by a fixation step and quantify them by PCR. Our contribution is to rapidly design an optimal primer set for the crucial PCR step. Our primer design reduces the level of experimental error as primers are highly similar in terms of physical properties and amplicon length. All 3C primers are compatible with multiplex PCR reactions. Primer uniqueness is checked genome-wide with a suitable index structure. CONCLUSIONS: In summary, our software 3PD facilitates genome-wide primer design for 3C experiments in a matter of seconds. Our software is available as a web server at: http://www.pristionchus.org/3CPrimerDesign/.

Pontine Tegmental Cap Dysplasia in an Extremely Preterm Infant and Review of the Literature.

Pontine tegmental cap dysplasia is a rare hindbrain malformation syndrome with a hypoplastic pons, a tissue protrusion into the fourth ventricle, and cranial nerve dysfunction. We here report clinical, imaging, and genetic findings of the first extremely low-birth-weight preterm infant with pontine tegmental cap dysplasia born at 25 weeks of gestation and provide an overview of 29 sporadic cases. A prenatally diagnosed hypoplastic and rostrally shifted cerebellum was indicative of a hindbrain defect and later identified as an early sign of pontine tegmental cap dysplasia in our patient. The neonate exhibited severe muscle hypotonia, persistent thermolability, and clinical signs of an involvement of facial, cochlear, and hypoglossal nerves. Furthermore, paroxysmal episodes of agonizing pain with facial tics, tonic and clonic muscle contractions, blepharospasm, and singultus are highlighted as new phenotypic features of pontine tegmental cap dysplasia. With our report, we present a severe case of pontine tegmental cap dysplasia and provide a brief overview of current knowledge on this rare disease.

reChIP-seq reveals widespread bivalency of H3K4me3 and H3K27me3 in CD4(+) memory T cells.

The combinatorial action of co-localizing chromatin modifications and regulators determines chromatin structure and function. However, identifying co-localizing chromatin features in a high-throughput manner remains a technical challenge. Here we describe a novel reChIP-seq approach and tailored bioinformatic analysis tool, normR that allows for the sequential enrichment and detection of co-localizing DNA-associated proteins in an unbiased and genome-wide manner. We illustrate the utility of the reChIP-seq method and normR by identifying H3K4me3 or H3K27me3 bivalently modified nucleosomes in primary human CD4(+) memory T cells. We unravel widespread bivalency at hypomethylated CpG-islands coinciding with inactive promoters of developmental regulators. reChIP-seq additionally uncovered heterogeneous bivalency in the population, which was undetectable by intersecting H3K4me3 and H3K27me3 ChIP-seq tracks. Finally, we provide evidence that bivalency is established and stabilized by an interplay between the genome and epigenome. Our reChIP-seq approach augments conventional ChIP-seq and is broadly applicable to unravel combinatorial modes of action.

Epigenetic dynamics of monocyte-to-macrophage differentiation.

BACKGROUND: Monocyte-to-macrophage differentiation involves major biochemical and structural changes. In order to elucidate the role of gene regulatory changes during this process, we used high-throughput sequencing to analyze the complete transcriptome and epigenome of human monocytes that were differentiated in vitro by addition of colony-stimulating factor 1 in serum-free medium. RESULTS: Numerous mRNAs and miRNAs were significantly up- or down-regulated. More than 100 discrete DNA regions, most often far away from transcription start sites, were rapidly demethylated by the ten eleven translocation enzymes, became nucleosome-free and gained histone marks indicative of active enhancers. These regions were unique for macrophages and associated with genes involved in the regulation of the actin cytoskeleton, phagocytosis and innate immune response. CONCLUSIONS: In summary, we have discovered a phagocytic gene network that is repressed by DNA methylation in monocytes and rapidly de-repressed after the onset of macrophage differentiation.

An oncogenic role for alternative NF-κB signaling in DLBCL revealed upon deregulated BCL6 expression.

Diffuse large B cell lymphoma (DLBCL) is a complex disease comprising diverse subtypes and genetic profiles. Possibly because of the prevalence of genetic alterations activating canonical NF-κB activity, a role for oncogenic lesions that activate the alternative NF-κB pathway in DLBCL has remained elusive. Here, we show that deletion/mutation of TRAF3, a negative regulator of the alternative NF-κB pathway, occurs in ∼15% of DLBCLs and that it often coexists with BCL6 translocation, which prevents terminal B cell differentiation. Accordingly, in a mouse model constitutive activation of the alternative NF-κB pathway cooperates with BCL6 deregulation in DLBCL development. This work demonstrates a key oncogenic role for the alternative NF-κB pathway in DLBCL development.

Identification of PENDRIN (SLC26A4) mutations in patients with congenital hypothyroidism and "apparent" thyroid dysgenesis.

CONTEXT: Congenital hypothyroidism, the most frequent endocrine congenital disease, can occur either based on a thyroid hormone biosynthesis defect or can predominantly be due to thyroid dysgenesis. However, a genetic cause could so far only be identified in less than 10% of patients with a thyroid dysgenesis. OBJECTIVES: Exome sequencing was used for the first time to find additional genetic defects in thyroid dysgenesis. PATIENTS AND METHODS: In a consanguineous family with thyroid dysgenesis, exome sequencing was applied, and findings were further validated by Sanger sequencing in a cohort of 94 patients with thyroid dysgenesis. RESULTS: By exome sequencing we identified a homozygous missense mutation (p.Leu597Ser) in the SLC26A4 gene of a patient with hypoplastic thyroid tissue, who was otherwise healthy. In the cohort of patients with thyroid dysgenesis, we observed a second case with a homozygous missense mutation (p.Gln413Arg) in the SLC26A4 gene, who was additionally affected by severe hearing problems. Both mutations were previously described as loss-of-function mutations in patients with Pendred syndrome and nonsyndromic enlarged vestibular aqueduct. CONCLUSION: We unexpectedly identified SLC26A4 mutations that were hitherto diagnosed in thyroid dyshormonogenesis patients, now for the first time in patients with structural thyroid defects. This result resembles the historic description of thyroid atrophy in patients with the so-called myxedematous form of cretinism after severe iodine deficiency. Most likely the thyroid defect of the two homozygous SLC26A4 gene mutation carriers represents a kind of secondary thyroid atrophy, rather than a primary defect of thyroid development in the sense of thyroid agenesis. Our study extends the variable clinical spectrum of patients with SLC26A4 mutations and points out the necessity to analyze the SLC26A4 gene in patients with apparent thyroid dysgenesis in addition to the known candidate genes TSHR, PAX8, NKX2.1, NKX2.5, and FOXE1.

Combined immunodeficiency develops with age in Immunodeficiency-centromeric instability-facial anomalies syndrome 2 (ICF2).

The autosomal recessive immunodeficiency-centromeric instability-facial anomalies syndrome (ICF) is characterized by immunodeficiency, developmental delay, and facial anomalies. ICF2, caused by biallelic ZBTB24 gene mutations, is acknowledged primarily as an isolated B-cell defect. Here, we extend the phenotype spectrum by describing, in particular, for the first time the development of a combined immune defect throughout the disease course as well as putative autoimmune phenomena such as granulomatous hepatitis and nephritis. We also demonstrate impaired cell-proliferation and increased cell death of immune and non-immune cells as well as data suggesting a chromosome separation defect in addition to the known chromosome condensation defect.

Integrative analysis revealed the molecular mechanism underlying RBM10-mediated splicing regulation.

RBM10 encodes an RNA binding protein. Mutations in RBM10 are known to cause multiple congenital anomaly syndrome in male humans, the TARP syndrome. However, the molecular function of RBM10 is unknown. Here we used PAR-CLIP to identify thousands of binding sites of RBM10 and observed significant RBM10-RNA interactions in the vicinity of splice sites. Computational analyses of binding sites as well as loss-of-function and gain-of-function experiments provided evidence for the function of RBM10 in regulating exon skipping and suggested an underlying mechanistic model, which could be subsequently validated by minigene experiments. Furthermore, we demonstrated the splicing defects in a patient carrying an RBM10 mutation, which could be explained by disrupted function of RBM10 in splicing regulation. Overall, our study established RBM10 as an important regulator of alternative splicing, presented a mechanistic model for RBM10-mediated splicing regulation and provided a molecular link to understanding a human congenital disorder.