Models for Natural Killer Cell Repertoire Formation

Natural killer (NK) cells lyse only cells that do not express sufficient levels of self class I MHC molecules. Inhibition of lysis is mediated by inhibitory receptors expressed by NK cells, such as the murine Ly49 receptors, that bind to MHC class I molecules. Since inhibitory receptor genes and MHC class I genes are located on different chromosomes, and are hence not automatically co-inherited, NK cells apparently adapt to the MHC environment during their development. Two models have been proposed to account for this “education” process of NK cells. The two-step selection model postulates that developing NK cells initiate the stable expression of a random set of Ly49 genes, and then undergo two selection steps, one for cells that express a sufficient number of self-MHC receptors, and one against cells that express too many inhibitory receptors. The sequential model postulates that a cell keeps initiating the stable expression of additional inhibitory receptors until a sufficient expression level of self-MHC specific receptors is reached, and the cell matures. In this study we implement both models in computer simulations, and compare simulation results to experimental data, in order to evaluate the relative plausibility of the two models

Probing Natural Killer Cell Education by Ly49 Receptor Expression Analysis and Computational Modelling in Single MHC Class I Mice

Murine natural killer (NK) cells express inhibitory Ly49 receptors for MHC class I molecules, which allows for “missing self” recognition of cells that downregulate MHC class I expression. During murine NK cell development, host MHC class I molecules impose an “educating impact” on the NK cell pool. As a result, mice with different MHC class I expression display different frequency distributions of Ly49 receptor combinations on NK cells. Two models have been put forward to explain this impact. The two-step selection model proposes a stochastic Ly49 receptor expression followed by selection for NK cells expressing appropriate receptor combinations. The sequential model, on the other hand, proposes that each NK cell sequentially expresses Ly49 receptors until an interaction of sufficient magnitude with self-class I MHC is reached for the NK cell to mature. With the aim to clarify which one of these models is most likely to reflect the actual biological process, we simulated the two educational schemes by mathematical modelling, and fitted the results to Ly49 expression patterns, which were analyzed in mice expressing single MHC class I molecules. Our results favour the two-step selection model over the sequential model. Furthermore, the MHC class I environment favoured maturation of NK cells expressing one or a few self receptors, suggesting a possible step of positive selection in NK cell education. Based on the predicted Ly49 binding preferences revealed by the model, we also propose, that Ly49 receptors are more promiscuous than previously thought in their interactions with MHC class I molecules, which was supported by functional studies of NK cell subsets expressing individual Ly49 receptors

Natural killer cell education in mice with single or multiple major histocompatibility complex class I molecules

The ability of murine NK cells to reject cells lacking self MHC class I expression results from an in vivo education process. To study the impact of individual MHC class I alleles on this process, we generated mice expressing single MHC class I alleles (Kb, Db, Dd, or Ld) or combinations of two or more alleles. All single MHC class I mice rejected MHC class I–deficient cells in an NK cell–dependent way. Expression of Kb or Dd conveyed strong rejection of MHC class I–deficient cells, whereas the expression of Db or Ld resulted in weaker responses. The educating impact of weak ligands (Db and Ld) was further attenuated by the introduction of additional MHC class I alleles, whereas strong ligands (Kb and Dd) maintained their educating impact under such conditions. An analysis of activating and inhibitory receptors in single MHC class I mice suggested that the educating impact of a given MHC class I molecule was controlled both by the number of NK cells affected and by the strength of each MHC class I–Ly49 receptor interaction, indicating that NK cell education may be regulated by a combination of qualitative and quantitative events

Classifying Medulloblastoma Subgroups Based on Small, Clinically Achievable Gene Sets

As treatment protocols for medulloblastoma (MB) are becoming subgroup-specific, means for reliably distinguishing between its subgroups are a timely need. Currently available methods include immunohistochemical stains, which are subjective and often inconclusive, and molecular techniques—e.g., NanoString, microarrays, or DNA methylation assays—which are time-consuming, expensive and not widely available. Quantitative PCR (qPCR) provides a good alternative for these methods, but the current NanoString panel which includes 22 genes is impractical for qPCR. Here, we applied machine-learning–based classifiers to extract reliable, concise gene sets for distinguishing between the four MB subgroups, and we compared the accuracy of these gene sets to that of the known NanoString 22-gene set. We validated our results using an independent microarray-based dataset of 92 samples of all four subgroups. In addition, we performed a qPCR validation on a cohort of 18 patients diagnosed with SHH, Group 3 and Group 4 MB. We found that the 22-gene set can be reduced to only six genes (IMPG2, NPR3, KHDRBS2, RBM24, WIF1, and EMX2) without compromising accuracy. The identified gene set is sufficiently small to make a qPCR-based MB subgroup classification easily accessible to clinicians, even in developing, poorly equipped countries

Reactive stroma and trastuzumab resistance in HER2-positive early breast cancer

We investigated the value of reactive stroma as a predictor for trastuzumab resistance in patients with early HER2-positive breast cancer receiving adjuvant therapy. The pathological reactive stroma and the mRNA gene signatures that reflect reactive stroma in 209 HER2-positive breast cancer samples from the FinHer adjuvant trial were evaluated. Levels of stromal gene signatures were determined as a continuous parameter, and pathological reactive stromal findings were defined as stromal predominant breast cancer (SPBC; >= 50% stromal) and correlated with distant disease-free survival. Gene signatures associated with reactive stroma in HER2-positive early breast cancer (N = 209) were significantly associated with trastuzumab resistance in estrogen receptor (ER)-negative tumors (hazard ratio [HR] = 1.27 p interaction = 0.014 [DCN], HR = 1.58, p interaction = 0.027 [PLAU], HR = 1.71, p interaction = 0.019 [HER2STROMA, novel HER2 stromal signature]), but not in ER-positive tumors (HR = 0.73 p interaction = 0.47 [DCN], HR = 0.71, p interaction = 0.73 [PLAU], HR = 0.84; p interaction = 0.36 [HER2STROMA]). Pathological evaluation of HER2-positive/ER-negative tumors suggested an association between SPBC and trastuzumab resistance. Reactive stroma did not correlate with tumor-infiltrating lymphocytes (TILs), and the expected benefit from trastuzumab in patients with high levels of TILs was pronounced only in tumors with low stromal reactivity (SPBCPeer reviewe

NuRD-mediated deacetylation of H3K27 facilitates recruitment of Polycomb Repressive Complex 2 to direct gene repression

The NURD and Polycomb complexes PRC1 and PRC2 have been implicated in stem cell differentiation although their molecular roles are unclear. This study identifies a common set of genes that are sequentially regulated by these complexes, the NURD complex deacetylates H3K27 as a prerequisite for subsequent methylation by PRC2

A regression analysis of gene expression in ES cells reveals two gene classes that are significantly different in epigenetic patterns

<p>Abstract</p> <p>Background</p> <p>To understand the gene regulatory system that governs the self-renewal and pluripotency of embryonic stem cells (ESCs) is an important step for promoting regenerative medicine. In it, the role of several core transcription factors (TFs), such as Oct4, Sox2 and Nanog, has been intensively investigated, details of their involvement in the genome-wide gene regulation are still not well clarified.</p> <p>Methods</p> <p>We constructed a predictive model of genome-wide gene expression in mouse ESCs from publicly available ChIP-seq data of 12 core TFs. The tag sequences were remapped on the genome by various alignment tools. Then, the binding density of each TF is calculated from the genome-wide bona fide TF binding sites. The TF-binding data was combined with the data of several epigenetic states (DNA methylation, several histone modifications, and CpG island) of promoter regions. These data as well as the ordinary peak intensity data were used as predictors of a simple linear regression model that predicts absolute gene expression. We also developed a pipeline for analyzing the effects of predictors and their interactions.</p> <p>Results</p> <p>Through our analysis, we identified two classes of genes that are either well explained or inefficiently explained by our model. The latter class seems to be genes that are not directly regulated by the core TFs. The regulatory regions of these gene classes show apparently distinct patterns of DNA methylation, histone modifications, existence of CpG islands, and gene ontology terms, suggesting the relative importance of epigenetic effects. Furthermore, we identified statistically significant TF interactions correlated with the epigenetic modification patterns.</p> <p>Conclusions</p> <p>Here, we proposed an improved prediction method in explaining the ESC-specific gene expression. Our study implies that the majority of genes are more or less directly regulated by the core TFs. In addition, our result is consistent with the general idea of relative importance of epigenetic effects in ESCs.</p

Nuclear hormone 1α,25-dihydroxyvitamin D3 elicits a genome-wide shift in the locations of VDR chromatin occupancy

A global understanding of the actions of the nuclear hormone 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3) and its vitamin D receptor (VDR) requires a genome-wide analysis of VDR binding sites. In THP-1 human monocytic leukemia cells we identified by ChIP-seq 2340 VDR binding locations, of which 1171 and 520 occurred uniquely with and without 1α,25(OH)2D3 treatment, respectively, while 649 were common. De novo identified direct repeat spaced by 3 nucleotides (DR3)-type response elements (REs) were strongly associated with the ligand-responsiveness of VDR occupation. Only 20% of the VDR peaks diminishing most after ligand treatment have a DR3-type RE, in contrast to 90% for the most growing peaks. Ligand treatment revealed 638 1α,25(OH)2D3 target genes enriched in gene ontology categories associated with immunity and signaling. From the 408 upregulated genes, 72% showed VDR binding within 400 kb of their transcription start sites (TSSs), while this applied only for 43% of the 230 downregulated genes. The VDR loci showed considerable variation in gene regulatory scenarios ranging from a single VDR location near the target gene TSS to very complex clusters of multiple VDR locations and target genes. In conclusion, ligand binding shifts the locations of VDR occupation to DR3-type REs that surround its target genes and occur in a large variety of regulatory constellations

Zeb1 controls neuron differentiation and germinal zone exit by a mesenchymal-epithelial-like transition

In the developing mammalian brain, differentiating neurons mature morphologically via neuronal polarity programs. Despite discovery of polarity pathways acting concurrently with differentiation, it's unclear how neurons traverse complex polarity transitions or how neuronal progenitors delay polarization during development. We report that zinc finger and homeobox transcription factor-1 (Zeb1), a master regulator of epithelial polarity, controls neuronal differentiation by transcriptionally repressing polarity genes in neuronal progenitors. Necessity-sufficiency testing and functional target screening in cerebellar granule neuron progenitors (GNPs) reveal that Zeb1 inhibits polarization and retains progenitors in their germinal zone (GZ). Zeb1 expression is elevated in the Sonic Hedgehog (SHH) medulloblastoma subgroup originating from GNPs with persistent SHH activation. Restored polarity signaling promotes differentiation and rescues GZ exit, suggesting a model for future differentiative therapies. These results reveal unexpected parallels between neuronal differentiation and mesenchymal-to-epithelial transition and suggest that active polarity inhibition contributes to altered GZ exit in pediatric brain cancers.National Institute of Neurological Disorders and Stroke grant: (1R01NS066936); March of Dimes Foundation grant: (#1-FY12-455).info:eu-repo/semantics/publishedVersio

Target Gene Analysis by Microarrays and Chromatin Immunoprecipitation Identifies HEY Proteins as Highly Redundant bHLH Repressors

HEY bHLH transcription factors have been shown to regulate multiple key steps in cardiovascular development. They can be induced by activated NOTCH receptors, but other upstream stimuli mediated by TGFß and BMP receptors may elicit a similar response. While the basic and helix-loop-helix domains exhibit strong similarity, large parts of the proteins are still unique and may serve divergent functions. The striking overlap of cardiac defects in HEY2 and combined HEY1/HEYL knockout mice suggested that all three HEY genes fulfill overlapping function in target cells. We therefore sought to identify target genes for HEY proteins by microarray expression and ChIPseq analyses in HEK293 cells, cardiomyocytes, and murine hearts. HEY proteins were found to modulate expression of their target gene to a rather limited extent, but with striking functional interchangeability between HEY factors. Chromatin immunoprecipitation revealed a much greater number of potential binding sites that again largely overlap between HEY factors. Binding sites are clustered in the proximal promoter region especially of transcriptional regulators or developmental control genes. Multiple lines of evidence suggest that HEY proteins primarily act as direct transcriptional repressors, while gene activation seems to be due to secondary or indirect effects. Mutagenesis of putative DNA binding residues supports the notion of direct DNA binding. While class B E-box sequences (CACGYG) clearly represent preferred target sequences, there must be additional and more loosely defined modes of DNA binding since many of the target promoters that are efficiently bound by HEY proteins do not contain an E-box motif. These data clearly establish the three HEY bHLH factors as highly redundant transcriptional repressors in vitro and in vivo, which explains the combinatorial action observed in different tissues with overlapping expression