Decoding complexity in biomolecular recognition of DNA i-motifs with microarrays.

DNA i-motifs (iMs) are non-canonical C-rich secondary structures implicated in numerous cellular processes. Though iMs exist throughout the genome, our understanding of iM recognition by proteins or small molecules is limited to a few examples. We designed a DNA microarray containing 10976 genomic iM sequences to examine the binding profiles of four iM-binding proteins, mitoxantrone and the iMab antibody. iMab microarray screens demonstrated that pH 6.5, 5% BSA buffer was optimal, and fluorescence was correlated with iM C-tract length. hnRNP K broadly recognizes diverse iM sequences, favoring 3-5 cytosine repeats flanked by thymine-rich loops of 1-3 nucleotides. Array binding mirrored public ChIP-Seq datasets, in which 35% of well-bound array iMs are enriched in hnRNP K peaks. In contrast, other reported iM-binding proteins had weaker binding or preferred G-quadruplex (G4) sequences instead. Mitoxantrone broadly binds both shorter iMs and G4s, consistent with an intercalation mechanism. These results suggest that hnRNP K may play a role in iM-mediated regulation of gene expression in vivo, whereas hnRNP A1 and ASF/SF2 are possibly more selective in their binding preferences. This powerful approach represents the most comprehensive investigation of how biomolecules selectively recognize genomic iMs to date.

bZIP Dimers CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun Prefer to Bind to Some Double-Stranded DNA Sequences Containing 5-Formylcytosine and 5-Carboxylcytosine.

In mammalian cells, 5-methylcytosine (5mC) occurs in genomic double-stranded DNA (dsDNA) and is enzymatically oxidized to 5-hydroxymethylcytosine (5hmC), then to 5-formylcytosine (5fC), and finally to 5-carboxylcytosine (5caC). These cytosine modifications are enriched in regulatory regions of the genome. The effect of these oxidative products on five bZIP dimers (CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun) binding to five types of dsDNA was measured using protein binding microarrays. The five dsDNAs contain either cytosine in both DNA strands or cytosine in one strand and either 5mC, 5hmC, 5fC, or 5caC in the second strand. Some sequences containing the CEBP half-site GCAA are bound more strongly by all five bZIP domains when dsDNA contains 5mC, 5hmC, or 5fC. dsDNA containing 5caC in some TRE (AP-1)-like sequences, e.g., TGACTAA, is better bound by Zta, ATF3|cJun, and cFos|cJun.

Custom G4 Microarrays Reveal Selective G-Quadruplex Recognition of Small Molecule BMVC: A Large-Scale Assessment of Ligand Binding Selectivity.

G-quadruplexes (G4) are considered new drug targets for human diseases such as cancer. More than 10,000 G4s have been discovered in human chromatin, posing challenges for assessing the selectivity of a G4-interactive ligand. 3,6-bis(1-Methyl-4-vinylpyridinium) carbazole diiodide (BMVC) is the first fluorescent small molecule for G4 detection in vivo. Our previous structural study shows that BMVC binds to the MYC promoter G4 (MycG4) with high specificity. Here, we utilize high-throughput, large-scale custom DNA G4 microarrays to analyze the G4-binding selectivity of BMVC. BMVC preferentially binds to the parallel MycG4 and selectively recognizes flanking sequences of parallel G4s, especially the 3'-flanking thymine. Importantly, the microarray results are confirmed by orthogonal NMR and fluorescence binding analyses. Our study demonstrates the potential of custom G4 microarrays as a platform to broadly and unbiasedly assess the binding selectivity of G4-interactive ligands, and to help understand the properties that govern molecular recognition.

Replacing C189 in the bZIP domain of Zta with S, T, V, or A changes DNA binding specificity to four types of double-stranded DNA.

Zta is a bZIP transcription factor (TF) in the Epstein-Barr virus that binds unmethylated and methylated DNA sequences. Substitution of cysteine 189 of Zta to serine (Zta(C189S)) results in a virus that is unable to execute the lytic cycle, which was attributed to a change in binding to methylated DNA sequences. To learn more about the role of this position in defining sequence-specific DNA binding, we mutated cysteine 189 to four other amino acids, producing Zta(C189S), Zta(C189T), Zta(C189A), and Zta(C189V) mutants. Zta and mutants were used in protein binding microarray (PBM) experiments to evaluate sequence-specific DNA binding to four types of double-stranded DNA (dsDNA): 1) with cytosine in both strands (DNA(C|C)), 2) with 5-methylcytosine (5mC) in one strand and cytosine in the second strand (DNA(5mC|C)), 3) with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand (DNA(5hmC|C)), and 4) with both cytosines in all CG dinucleotides containing 5mC (DNA(5mCG)). Zta(C189S) and Zta(C189T) bound the TRE (AP-1) motif (TGAG/CTCA) more strongly than wild-type Zta, while binding to other sequences, including the C/EBP half site GCAA was reduced. Binding of Zta(C189S) and Zta(C189T) to DNA containing modified cytosines (DNA(5mC|C), DNA(5hmC|C), and DNA(5mCG)) was reduced compared to Zta. Zta(C189A) and Zta(C189V) had higher non-specific binding to all four types of DNA. Our data suggests that position C189 in Zta impacts sequence-specific binding to DNA containing modified and unmodified cytosine.

The Epstein-Barr Virus B-ZIP Protein Zta Recognizes Specific DNA Sequences Containing 5-Methylcytosine and 5-Hydroxymethylcytosine.

The Epstein-Barr virus (EBV) B-ZIP transcription factor Zta binds to many DNA sequences containing methylated CG dinucleotides. Using protein binding microarrays (PBMs), we analyzed the sequence specific DNA binding of Zta to four kinds of double-stranded DNA (dsDNA): (1) DNA containing cytosine in both strands, (2) DNA with 5-methylcytosine (5mC) in one strand and cytosine in the second strand, (3) DNA with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand, and (4) DNA in which both cytosines in all CG dinucleotides contain 5mC. We compared these data to PBM data for three additional B-ZIP proteins (CREB1 and CEBPB homodimers and cJun|cFos heterodimers). With cytosine, Zta binds the TRE motif TGAC/GTCA as previously reported. With CG dinucleotides containing 5mC on both strands, many TRE motif variants containing a methylated CG dinucleotide at two positions in the motif, such as MGAGTCA and TGAGMGA (where M = 5mC), were preferentially bound. 5mC inhibits binding of Zta to both TRE motif half-sites GTCA and CTCA. Like the CREB1 homodimer, the Zta homodimer and the cJun|cFos heterodimer more strongly bind the C/EBP half-site tetranucleotide GCAA when it contains 5mC. Zta also binds dsDNA sequences containing 5hmC in one strand, although the effect is less dramatic than that observed for 5mC. Our results identify new DNA sequences that are well-bound by the viral B-ZIP protein Zta only when they contain 5mC or 5hmC, uncovering the potential for discovery of new viral and host regulatory programs controlled by EBV.

5-Methylcytosine (5mC) and 5-Hydroxymethylcytosine (5hmC) Enhance the DNA Binding of CREB1 to the C/EBP Half-Site Tetranucleotide GCAA.

In human and mouse stem cells and brain, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) can occur outside of CG dinucleotides. Using protein binding microarrays (PBMs) containing 60-mer DNA probes, we evaluated the effect of 5mC and 5hmC on one DNA strand on the double-stranded DNA binding of the mouse B-ZIP transcription factors (TFs) CREB1, ATF1, and JUND. 5mC inhibited binding of CREB1 to the canonical CRE half-site |GTCA but enhanced binding to the C/EBP half-site |GCAA. 5hmC inhibited binding of CREB1 to all 8-mers except TGAT|GCAA, where binding is enhanced. We observed similar DNA binding patterns with ATF1, a closely related B-ZIP domain. In contrast, both 5mC and 5hmC inhibited binding of JUND. These results identify new DNA sequences that are well-bound by CREB1 and ATF1 only when they contain 5mC or 5hmC. Analysis of two X-ray structures examines the consequences of 5mC and 5hmC on DNA binding by CREB and FOS|JUN.

Inheritance of Cytosine Methylation.

There are numerous examples of parental transgenerational inheritance that is epigenetic. The informational molecules include RNA, chromatin modifications, and cytosine methylation. With advances in DNA sequencing technologies, the molecular and epigenetic mechanisms mediating these effects are now starting to be uncovered. This mini-review will highlight some of the examples of epigenetic inheritance, the establishment of cytosine methylation in sperm, and recent genomic studies linking sperm cytosine methylation to epigenetic effects on offspring. A recent paper examining changes in diet and sperm cytosine methylation from pools of eight animals each, found differences between a normal diet, a high fat diet, and a low protein diet. However, epivariation between individuals within a group was greater than the differences between groups obscuring any potential methylation changes linked to diet. Learning more about epivariation may help unravel the mechanisms that regulate cytosine methylation. In addition, other experimental and genetic systems may also produce more dramatic changes in the sperm methylome, making it easier to unravel potential transgenerational phenomena. J. Cell. Physiol. 231: 2346-2352, 2016. © 2016 Wiley Periodicals, Inc.

A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters.

The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially approximately 100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.

Nucleosome maps of the human cytomegalovirus genome reveal a temporal switch in chromatin organization linked to a major IE protein.

Human CMV (hCMV) establishes lifelong infections in most of us, causing developmental defects in human embryos and life-threatening disease in immunocompromised individuals. During productive infection, the viral >230,000-bp dsDNA genome is expressed widely and in a temporal cascade. The hCMV genome does not carry histones when encapsidated but has been proposed to form nucleosomes after release into the host cell nucleus. Here, we present hCMV genome-wide nucleosome occupancy and nascent transcript maps during infection of permissive human primary cells. We show that nucleosomes occupy nuclear viral DNA in a nonrandom and highly predictable fashion. At early times of infection, nucleosomes associate with the hCMV genome largely according to their intrinsic DNA sequence preferences, indicating that initial nucleosome formation is genetically encoded in the virus. However, as infection proceeds to the late phase, nucleosomes redistribute extensively to establish patterns mostly determined by nongenetic factors. We propose that these factors include key regulators of viral gene expression encoded at the hCMV major immediate-early (IE) locus. Indeed, mutant virus genomes deficient for IE1 expression exhibit globally increased nucleosome loads and reduced nucleosome dynamics compared with WT genomes. The temporal nucleosome occupancy differences between IE1-deficient and WT viruses correlate inversely with changes in the pattern of viral nascent and total transcript accumulation. These results provide a framework of spatial and temporal nucleosome organization across the genome of a major human pathogen and suggest that an hCMV major IE protein governs overall viral chromatin structure and function.

The DNA-encoded nucleosome organization of a eukaryotic genome.

Nucleosome organization is critical for gene regulation. In living cells this organization is determined by multiple factors, including the action of chromatin remodellers, competition with site-specific DNA-binding proteins, and the DNA sequence preferences of the nucleosomes themselves. However, it has been difficult to estimate the relative importance of each of these mechanisms in vivo, because in vivo nucleosome maps reflect the combined action of all influencing factors. Here we determine the importance of nucleosome DNA sequence preferences experimentally by measuring the genome-wide occupancy of nucleosomes assembled on purified yeast genomic DNA. The resulting map, in which nucleosome occupancy is governed only by the intrinsic sequence preferences of nucleosomes, is similar to in vivo nucleosome maps generated in three different growth conditions. In vitro, nucleosome depletion is evident at many transcription factor binding sites and around gene start and end sites, indicating that nucleosome depletion at these sites in vivo is partly encoded in the genome. We confirm these results with a micrococcal nuclease-independent experiment that measures the relative affinity of nucleosomes for approximately 40,000 double-stranded 150-base-pair oligonucleotides. Using our in vitro data, we devise a computational model of nucleosome sequence preferences that is significantly correlated with in vivo nucleosome occupancy in Caenorhabditis elegans. Our results indicate that the intrinsic DNA sequence preferences of nucleosomes have a central role in determining the organization of nucleosomes in vivo.

The homeobox transcription factor DUXBL controls exit from totipotency.

In mice, exit from the totipotent two-cell (2C) stage embryo requires silencing of the 2C-associated transcriptional program. However, the molecular mechanisms involved in this process remain poorly understood. Here we demonstrate that the 2C-specific transcription factor double homeobox protein (DUX) mediates an essential negative feedback loop by inducing the expression of DUXBL to promote this silencing. We show that DUXBL gains accessibility to DUX-bound regions specifically upon DUX expression. Furthermore, we determine that DUXBL interacts with TRIM24 and TRIM33, members of the TRIM superfamily involved in gene silencing, and colocalizes with them in nuclear foci upon DUX expression. Importantly, DUXBL overexpression impairs 2C-associated transcription, whereas Duxbl inactivation in mouse embryonic stem cells increases DUX-dependent induction of the 2C-transcriptional program. Consequently, DUXBL deficiency in embryos results in sustained expression of 2C-associated transcripts leading to early developmental arrest. Our study identifies DUXBL as an essential regulator of totipotency exit enabling the first divergence of cell fates.

Microenvironment shapes small-cell lung cancer neuroendocrine states and presents therapeutic opportunities.

Small-cell lung cancer (SCLC) is the most fatal form of lung cancer. Intratumoral heterogeneity, marked by neuroendocrine (NE) and non-neuroendocrine (non-NE) cell states, defines SCLC, but the cell-extrinsic drivers of SCLC plasticity are poorly understood. To map the landscape of SCLC tumor microenvironment (TME), we apply spatially resolved transcriptomics and quantitative mass spectrometry-based proteomics to metastatic SCLC tumors obtained via rapid autopsy. The phenotype and overall composition of non-malignant cells in the TME exhibit substantial variability, closely mirroring the tumor phenotype, suggesting TME-driven reprogramming of NE cell states. We identify cancer-associated fibroblasts (CAFs) as a crucial element of SCLC TME heterogeneity, contributing to immune exclusion, and predicting exceptionally poor prognosis. Our work provides a comprehensive map of SCLC tumor and TME ecosystems, emphasizing their pivotal role in SCLC's adaptable nature, opening possibilities for reprogramming the TME-tumor communications that shape SCLC tumor states.

DECODING COMPLEXITY IN BIOMOLECULAR RECOGNITION OF DNA I-MOTIFS.

DNA i-motifs (iMs) are non-canonical C-rich secondary structures implicated in numerous cellular processes. Though iMs exist throughout the genome, our understanding of iM recognition by proteins or small molecules is limited to a few examples. We designed a DNA microarray containing 10,976 genomic iM sequences to examine the binding profiles of four iM-binding proteins, mitoxantrone, and the iMab antibody. iMab microarray screens demonstrated that pH 6.5, 5% BSA buffer was optimal, and fluorescence was correlated with iM C-tract length. hnRNP K broadly recognizes diverse iM sequences, favoring 3-5 cytosine repeats flanked by thymine-rich loops of 1-3 nucleotides. Array binding mirrored public ChIP-Seq datasets, in which 35% of well-bound array iMs are enriched in hnRNP K peaks. In contrast, other reported iM-binding proteins had weaker binding or preferred G-quadruplex (G4) sequences instead. Mitoxantrone broadly binds both shorter iMs and G4s, consistent with an intercalation mechanism. These results suggest that hnRNP K may play a role in iM-mediated regulation of gene expression in vivo, whereas hnRNP A1 and ASF/SF2 are possibly more selective in their binding preferences. This powerful approach represents the most comprehensive investigation of how biomolecules selectively recognize genomic iMs to date.

Direct ionic stress sensing and mitigation by the transcription factor NFAT5.

Homeostatic control of intracellular ionic strength is essential for protein, organelle and genome function, yet mechanisms that sense and enable adaptation to ionic stress remain poorly understood in animals. We find that the transcription factor NFAT5 directly senses solution ionic strength using a C-terminal intrinsically disordered region. Both in intact cells and in a purified system, NFAT5 forms dynamic, reversible biomolecular condensates in response to increasing ionic strength. This self-associative property, conserved from insects to mammals, allows NFAT5 to accumulate in the nucleus and activate genes that restore cellular ion content. Mutations that reduce condensation or those that promote aggregation both reduce NFAT5 activity, highlighting the importance of optimally tuned associative interactions. Remarkably, human NFAT5 alone is sufficient to reconstitute a mammalian transcriptional response to ionic or hypertonic stress in yeast. Thus NFAT5 is both the sensor and effector of a cell-autonomous ionic stress response pathway in animal cells.

A novel GATA2 distal enhancer mutation results in MonoMAC syndrome in 2 second cousins.

Mutations in the transcription factor GATA2 can cause MonoMAC syndrome, a GATA2 deficiency disease characterized by several findings, including disseminated nontuberculous mycobacterial infections, severe deficiencies of monocytes, natural killer cells, and B lymphocytes, and myelodysplastic syndrome. GATA2 mutations are found in âˆ¼90% of patients with a GATA2 deficiency phenotype and are largely missense mutations in the conserved second zinc-finger domain. Mutations in an intron 5 regulatory enhancer element are also well described in GATA2 deficiency. Here, we present a multigeneration kindred with the clinical features of GATA2 deficiency but lacking an apparent GATA2 mutation. Whole genome sequencing revealed a unique adenine-to-thymine variant in the GATA2 -110 enhancer 116,855 bp upstream of the GATA2 ATG start site. The mutation creates a new E-box consensus in position with an existing GATA-box to generate a new hematopoietic regulatory composite element. The mutation segregates with the disease in several generations of the family. Cell type-specific allelic imbalance of GATA2 expression was observed in the bone marrow of a patient with higher expression from the mutant-linked allele. Allele-specific overexpression of GATA2 was observed in CRISPR/Cas9-modified HL-60 cells and in luciferase assays with the enhancer mutation. This study demonstrates overexpression of GATA2 resulting from a single nucleotide change in an upstream enhancer element in patients with MonoMAC syndrome. Patients in this study were enrolled in the National Institute of Allergy and Infectious Diseases clinical trial and the National Cancer Institute clinical trial (both trials were registered at www.clinicaltrials.gov as #NCT01905826 and #NCT01861106, respectively).

Altering the Double-Stranded DNA Specificity of the bZIP Domain of Zta with Site-Directed Mutagenesis at N182.

Zta, the Epstein-Barr virus bZIP transcription factor (TF), binds both unmethylated and methylated double-stranded DNA (dsDNA) in a sequence-specific manner. We studied the contribution of a conserved asparagine (N182) to sequence-specific dsDNA binding to four types of dsDNA: (i) dsDNA with cytosine in both strands ((DNA(C|C)), (ii, iii) dsDNA with 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and cytosine in the second strand ((DNA(5mC|C) and DNA(5hmC|C)), and (iv) dsDNA with methylated cytosine in both strands in all CG dinucleotides ((DNA(5mCG)). We replaced asparagine with five similarly sized amino acids (glutamine (Q), serine (S), threonine (T), isoleucine (I), or valine (V)) and used protein binding microarrays to evaluate sequence-specific dsDNA binding. Zta preferentially binds the pseudo-palindrome TRE (AP1) motif (T-4G-3A-2G/C 0T2C3A4 ). Zta (N182Q) changes binding to A3 in only one half-site. Zta(N182S) changes binding to G3 in one or both halves of the motif. Zta(N182S) and Zta(N182Q) have 34- and 17-fold weaker median dsDNA binding, respectively. Zta(N182V) and Zta(N182I) have increased binding to dsDNA(5mC|C). Molecular dynamics simulations rationalize some of these results, identifying hydrogen bonds between glutamine and A3 , but do not reveal why serine preferentially binds G3 , suggesting that entropic interactions may mediate this new binding specificity.

ASXL1 and STAG2 are common mutations in GATA2 deficiency patients with bone marrow disease and myelodysplastic syndrome.

Patients with GATA2 deficiencyharbor de novo or inherited germline mutations in the GATA2 transcription factor gene, predisposing them to myeloid malignancies. There is considerable variation in disease progression, even among family members with the same mutation in GATA2. We investigated somatic mutations in 106 patients with GATA2 deficiency to identify acquired mutations that are associated with myeloid malignancies. Myelodysplastic syndrome (MDS) was the most common diagnosis (∼44%), followed by GATA2 bone marrow immunodeficiency disorder (G2BMID; ∼37%). Thirteen percent of the cohort had GATA2 mutations but displayed no disease manifestations. There were no correlations between age or sex with disease progression or survival. Cytogenetic analyses showed a high incidence of abnormalities (∼43%), notably trisomy 8 (∼23%) and monosomy 7 (∼12%), but the changes did not correlate with lower survival. Somatic mutations in ASXL1 and STAG2 were detected in ∼25% of patients, although the mutations were rarely concomitant. Mutations in DNMT3A were found in ∼10% of patients. These somatic mutations were found similarly in G2BMID and MDS, suggesting clonal hematopoiesis in early stages of disease, before the onset of MDS. ASXL1 mutations conferred a lower survival probability and were more prevalent in female patients. STAG2 mutations also conferred a lower survival probability, but did not show a statistically significant sex bias. There was a conspicuous absence of many commonly mutated genes associated with myeloid malignancies, including TET2, IDH1/2, and the splicing factor genes. Notably, somatic mutations in chromatin-related genes and cohesin genes characterized disease progression in GATA2 deficiency.

Spatial meta-transcriptomics reveal associations of intratumor bacteria burden with lung cancer cells showing a distinct oncogenic signature.

BACKGROUND: The lung intratumor microbiome influences lung cancer tumorigenesis and treatment responses, but detailed data on the extent, location, and effects of microbes within lung tumors are missing, information needed for improved prognosis and treatment. METHODS: To address this gap, we developed a novel spatial meta-transcriptomic method simultaneously detecting the expression level of 1,811 host genes and 3 microbe targets (bacteria, fungi, and cytomegalovirus). After rigorous validation, we analyzed the spatial meta-transcriptomic profiles of tumor cells, T cells, macrophages, other immune cells, and stroma in surgically resected tumor samples from 12 patients with early-stage lung cancer. RESULTS: Bacterial burden was significantly higher in tumor cells compared with T cells, macrophages, other immune cells, and stroma. This burden increased from tumor-adjacent normal lung and tertiary lymphoid structures to tumor cells to the airways, suggesting that lung intratumor bacteria derive from the latter route of entry. Expression of oncogenic ß-catenin was strongly correlated with bacterial burden, as were tumor histological subtypes and environmental factors. CONCLUSIONS: Intratumor bacteria were enriched with tumor cells and associated with multiple oncogenic pathways, supporting a rationale for reducing the local intratumor microbiome in lung cancer for patient benefit. TRIAL REGISTRATION NUMBER: NCT00242723, NCT02146170.

A Knock-In Mouse Model of Thymoma With the GTF2I L424H Mutation.

INTRODUCTION: The pathogenesis of thymic epithelial tumors remains largely unknown. We previously identified GTF2I L424H as the most frequently recurrent mutation in thymic epithelial tumors. Nevertheless, the precise role of this mutation in tumorigenesis of thymic epithelial cells is unclear. METHODS: To investigate the role of GTF2I L424H mutation in thymic epithelial cells in vivo, we generated and characterized a mouse model in which the Gtf2i L424H mutation was conditionally knocked-in in the Foxn1+ thymic epithelial cells. Digital spatial profiling was performed on thymomas and normal thymic tissues with GeoMx-mouse whole transcriptome atlas. Immunohistochemistry staining was performed using both mouse tissues and human thymic epithelial tumors. RESULTS: We observed that the Gtf2i mutation impairs development of the thymic medulla and maturation of medullary thymic epithelial cells in young mice and causes tumor formation in the thymus of aged mice. Cell cycle-related pathways, such as E2F targets and MYC targets, are enriched in the tumor epithelial cells. Results of gene set variation assay analysis revealed that gene signatures of cortical thymic epithelial cells and thymic epithelial progenitor cells are also enriched in the thymomas of the knock-in mice, which mirrors the human counterparts in The Cancer Genome Atlas database. Immunohistochemistry results revealed similar expression pattern of epithelial cell markers between mouse and human thymomas. CONCLUSIONS: We have developed and characterized a novel thymoma mouse model. This study improves knowledge of the molecular drivers in thymic epithelial cells and provides a tool for further study of the biology of thymic epithelial tumors and for development of novel therapies.

REL Domain of NFATc2 Binding to Five Types of DNA Using Protein Binding Microarrays.

NFATc2 is a DNA binding protein in the Rel family transcription factors, which binds a CGGAA motif better when both cytosines in the CG dinucleotide are methylated. Using protein binding microarrays (PBMs), we examined the DNA binding of NFATc2 to three additional types of DNA: single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) with either 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and a cytosine in the second strand. ATTTCCAC, the complement of the core GGAA motif, is better bound as ssDNA compared to dsDNA. dsDNA containing the 5-mer CGGAA with either 5mC or 5hmC in one DNA strand is bound stronger than CGGAA. In contrast, the reverse complement TTCCG is bound weaker when it contains 5mC. Analysis of the available NFATc2:dsDNA complexes rationalizes these PBM data.