Shannon McWeeney

Defining the CREB Regulon: A Genome-Wide Analysis of Transcription Factor Regulatory Regions

The CREB transcription factor regulates differentiation, survival, and synaptic plasticity. The complement of CREB targets responsible for these responses has not been identified, however. We developed a novel approach to identify CREB targets, termed serial analysis of chromatin occupancy (SACO), by combining chromatin immunoprecipitation (ChIP) with a modification of SAGE. Using a SACO library derived from rat PC12 cells, we identified approximately 41,000 genomic signature tags (GSTs) that mapped to unique genomic loci. CREB binding was confirmed for all loci supported by multiple GSTs. Of the 6302 loci identified by multiple GSTs, 40% were within 2 kb of the transcriptional start of an annotated gene, 49% were within 1 kb of a CpG island, and 72% were within 1 kb of a putative cAMP-response element (CRE). A large fraction of the SACO loci delineated bidirectional promoters and novel antisense transcripts. This study represents the most comprehensive definition of transcription factor binding sites in a metazoan species.

Oncogenic CSF3R mutations in chronic neutrophilic leukemia and atypical CML

BACKGROUND The molecular causes of many hematologic cancers remain unclear. Among these cancers are chronic neutrophilic leukemia (CNL) and atypical (BCR-ABL1-negative) chronic myeloid leukemia (CML), both of which are diagnosed on the basis of neoplastic expansion of granulocytic cells and exclusion of genetic drivers that are known to occur in other myeloproliferative neoplasms and myeloproliferative-myelodysplastic overlap neoplasms. METHODS To identify potential genetic drivers in these disorders, we used an integrated approach of deep sequencing coupled with the screening of primary leukemia cells obtained from patients with CNL or atypical CML against panels of tyrosine kinase-specific small interfering RNAs or small-molecule kinase inhibitors. We validated candidate oncogenes using in vitro transformation assays, and drug sensitivities were validated with the use of assays of primary-cell colonies. RESULTS We identified activating mutations in the gene encoding the receptor for colony-stimulating factor 3 (CSF3R) in 16 of 27 patients (59%) with CNL or atypical CML. These mutations segregate within two distinct regions of CSF3R and lead to preferential downstream kinase signaling through SRC family-TNK2 or JAK kinases and differential sensitivity to kinase inhibitors. A patient with CNL carrying a JAK-activating CSF3R mutation had marked clinical improvement after the administration of the JAK1/2 inhibitor ruxolitinib. CONCLUSIONS Mutations in CSF3R are common in patients with CNL or atypical CML and represent a potentially useful criterion for diagnosing these neoplasms. (Funded by the Leukemia and Lymphoma Society and others.).

MIFlowCyt: The Minimum Information About a Flow Cytometry Experiment

A fundamental tenet of scientific research is that published results are open to independent validation and refutation. Minimum data standards aid data providers, users, and publishers by providing a specification of what is required to unambiguously interpret experimental findings. Here, we present the Minimum Information about a Flow Cytometry Experiment (MIFlowCyt) standard, stating the minimum information required to report flow cytometry (FCM) experiments. We brought together a cross‐disciplinary international collaborative group of bioinformaticians, computational statisticians, software developers, instrument manufacturers, and clinical and basic research scientists to develop the standard. The standard was subsequently vetted by the International Society for Advancement of Cytometry (ISAC) Data Standards Task Force, Standards Committee, membership, and Council. The MIFlowCyt standard includes recommendations about descriptions of the specimens and reagents included in the FCM experiment, the configuration of the instrument used to perform the assays, and the data processing approaches used to interpret the primary output data. MIFlowCyt has been adopted as a standard by ISAC, representing the FCM scientific community including scientists as well as software and hardware manufacturers. Adoptionof MIFlowCyt by the scientific and publishing communities will facilitate third‐party understanding and reuse of FCM data.

A Viral microRNA Down-Regulates Multiple Cell Cycle Genes through mRNA 5′UTRs

Global gene expression data combined with bioinformatic analysis provides strong evidence that mammalian miRNAs mediate repression of gene expression primarily through binding sites within the 3′ untranslated region (UTR). Using RNA induced silencing complex immunoprecipitation (RISC-IP) techniques we have identified multiple cellular targets for a human cytomegalovirus (HCMV) miRNA, miR-US25-1. Strikingly, this miRNA binds target sites primarily within 5′UTRs, mediating significant reduction in gene expression. Intriguingly, many of the genes targeted by miR-US25-1 are associated with cell cycle control, including cyclin E2, BRCC3, EID1, MAPRE2, and CD147, suggesting that miR-US25-1 is targeting genes within a related pathway. Deletion of miR-US25-1 from HCMV results in over expression of cyclin E2 in the context of viral infection. Our studies demonstrate that a viral miRNA mediates translational repression of multiple cellular genes by targeting mRNA 5′UTRs.

Evaluating gene expression in C57BL/6J and DBA/2J mouse striatum using RNA-Seq and microarrays.

C57BL/6J (B6) and DBA/2J (D2) are two of the most commonly used inbred mouse strains in neuroscience research. However, the only currently available mouse genome is based entirely on the B6 strain sequence. Subsequently, oligonucleotide microarray probes are based solely on this B6 reference sequence, making their application for gene expression profiling comparisons across mouse strains dubious due to their allelic sequence differences, including single nucleotide polymorphisms (SNPs). The emergence of next-generation sequencing (NGS) and the RNA-Seq application provides a clear alternative to oligonucleotide arrays for detecting differential gene expression without the problems inherent to hybridization-based technologies. Using RNA-Seq, an average of 22 million short sequencing reads were generated per sample for 21 samples (10 B6 and 11 D2), and these reads were aligned to the mouse reference genome, allowing 16,183 Ensembl genes to be queried in striatum for both strains. To determine differential expression, ‘digital mRNA counting’ is applied based on reads that map to exons. The current study compares RNA-Seq (Illumina GA IIx) with two microarray platforms (Illumina MouseRef-8 v2.0 and Affymetrix MOE 430 2.0) to detect differential striatal gene expression between the B6 and D2 inbred mouse strains. We show that by using stringent data processing requirements differential expression as determined by RNA-Seq is concordant with both the Affymetrix and Illumina platforms in more instances than it is concordant with only a single platform, and that instances of discordance with respect to direction of fold change were rare. Finally, we show that additional information is gained from RNA-Seq compared to hybridization-based techniques as RNA-Seq detects more genes than either microarray platform. The majority of genes differentially expressed in RNA-Seq were only detected as present in RNA-Seq, which is important for studies with smaller effect sizes where the sensitivity of hybridization-based techniques could bias interpretation.

Identification of β-catenin binding regions in colon cancer cells using ChIP-Seq

Deregulation of the Wnt/β-catenin signaling pathway is a hallmark of colon cancer. Mutations in the adenomatous polyposis coli (APC) gene occur in the vast majority of colorectal cancers and are an initiating event in cellular transformation. Cells harboring mutant APC contain elevated levels of the β-catenin transcription coactivator in the nucleus which leads to abnormal expression of genes controlled by β-catenin/T-cell factor 4 (TCF4) complexes. Here, we use chromatin immunoprecipitation coupled with massively parallel sequencing (ChIP-Seq) to identify β-catenin binding regions in HCT116 human colon cancer cells. We localized 2168 β-catenin enriched regions using a concordance approach for integrating the output from multiple peak alignment algorithms. Motif discovery algorithms found a core TCF4 motif (T/A–T/A–C–A–A–A–G), an extended TCF4 motif (A/T/G–C/G–T/A–T/A–C–A–A–A–G) and an AP-1 motif (T–G–A–C/T–T–C–A) to be significantly represented in β-catenin enriched regions. Furthermore, 417 regions contained both TCF4 and AP-1 motifs. Genes associated with TCF4 and AP-1 motifs bound β-catenin, TCF4 and c-Jun in vivo and were activated by Wnt signaling and serum growth factors. Our work provides evidence that Wnt/β-catenin and mitogen signaling pathways intersect directly to regulate a defined set of target genes.

Cancer Stem Cell Tumor Model Reveals Invasive Morphology and Increased Phenotypical Heterogeneity

The recently developed concept of cancer stem cells (CSC) sheds new light on various aspects of tumor growth and progression. Here, we present a mathematical model of malignancies to investigate how a hierarchical organized cancer cell population affects the fundamental properties of solid malignancies. We establish that tumors modeled in a CSC context more faithfully resemble human malignancies and show invasive behavior, whereas tumors without a CSC hierarchy do not. These findings are corroborated by in vitro studies. In addition, we provide evidence that the CSC model is accompanied by highly altered evolutionary dynamics compared with the ones predicted to exist in a stochastic, nonhierarchical tumor model. Our main findings indicate that the CSC model allows for significantly higher tumor heterogeneity, which may affect therapy resistance. Moreover, we show that therapy which fails to target the CSC population is not only unsuccessful in curing the patient, but also promotes malignant features in the recurring tumor. These include rapid expansion, increased invasion, and enhanced heterogeneity.

Modeling Host Genetic Regulation of Influenza Pathogenesis in the Collaborative Cross

Genetic variation contributes to host responses and outcomes following infection by influenza A virus or other viral infections. Yet narrow windows of disease symptoms and confounding environmental factors have made it difficult to identify polymorphic genes that contribute to differential disease outcomes in human populations. Therefore, to control for these confounding environmental variables in a system that models the levels of genetic diversity found in outbred populations such as humans, we used incipient lines of the highly genetically diverse Collaborative Cross (CC) recombinant inbred (RI) panel (the pre-CC population) to study how genetic variation impacts influenza associated disease across a genetically diverse population. A wide range of variation in influenza disease related phenotypes including virus replication, virus-induced inflammation, and weight loss was observed. Many of the disease associated phenotypes were correlated, with viral replication and virus-induced inflammation being predictors of virus-induced weight loss. Despite these correlations, pre-CC mice with unique and novel disease phenotype combinations were observed. We also identified sets of transcripts (modules) that were correlated with aspects of disease. In order to identify how host genetic polymorphisms contribute to the observed variation in disease, we conducted quantitative trait loci (QTL) mapping. We identified several QTL contributing to specific aspects of the host response including virus-induced weight loss, titer, pulmonary edema, neutrophil recruitment to the airways, and transcriptional expression. Existing whole-genome sequence data was applied to identify high priority candidate genes within QTL regions. A key host response QTL was located at the site of the known anti-influenza Mx1 gene. We sequenced the coding regions of Mx1 in the eight CC founder strains, and identified a novel Mx1 allele that showed reduced ability to inhibit viral replication, while maintaining protection from weight loss.

A gene expression signature of CD34+ cells to predict major cytogenetic response in chronic-phase chronic myeloid leukemia patients treated with imatinib

In chronic-phase chronic myeloid leukemia (CML) patients, the lack of a major cytogenetic response (< 36% Ph(+) metaphases) to imatinib within 12 months indicates failure and mandates a change of therapy. To identify biomarkers predictive of imatinib failure, we performed gene expression array profiling of CD34(+) cells from 2 independent cohorts of imatinib-naive chronic-phase CML patients. The learning set consisted of retrospectively selected patients with a complete cytogenetic response or more than 65% Ph(+) metaphases within 12 months of imatinib therapy. Based on analysis of variance P less than .1 and fold difference 1.5 or more, we identified 885 probe sets with differential expression between responders and nonresponders, from which we extracted a 75-probe set minimal signature (classifier) that separated the 2 groups. On application to a prospectively accrued validation set, the classifier correctly predicted 88% of responders and 83% of nonresponders. Bioinformatics analysis and comparison with published studies revealed overlap of classifier genes with CML progression signatures and implicated beta-catenin in their regulation, suggesting that chronic-phase CML patients destined to fail imatinib have more advanced disease than evident by morphologic criteria. Our classifier may allow directing more aggressive therapy upfront to the patients most likely to benefit while sparing good-risk patients from unnecessary toxicity.

Serial analysis of chromatin occupancy identifies β-catenin target genes in colorectal carcinoma cells

Most instances of colorectal cancer are due to abnormalities in the Wnt signaling pathway, resulting in nuclear accumulation of β-catenin. β-Catenin activates transcription of target genes primarily by associating with the T cell factor/lymphoid enhancer-binding factor (TCF/Lef) family of transcription factors. In this report, we use serial analysis of chromatin occupancy (SACO) to identify 412 high-confidence β-catenin targets in HCT116 colorectal carcinoma cells. Of these targets, 84% contained a consensus TCF motif and were occupied by TCF4 in vivo. Examination of the flanking 5-bp residues in each consensus revealed motif-specific enrichment at neighboring sites. β-Catenin binding was localized to the 5′ promoters, internal regions, and 3′ UTRs of protein-coding genes. Furthermore, 15 components of the canonical Wnt pathway were identified as β-catenin target genes, suggesting that feed-forward and feedback mechanisms exist to modulate the Wnt signal in colon cancer cells.