Seth Frietze

Genetic analysis of Ikaros target genes and tumor suppressor function in BCR-ABL1+ pre–B ALL

Inactivation of the tumor suppressor gene encoding the transcriptional regulator Ikaros (IKZF1) is a hallmark of BCR-ABL1+ precursor B cell acute lymphoblastic leukemia (pre–B ALL). However, the mechanisms by which Ikaros functions as a tumor suppressor in pre–B ALL remain poorly understood. Here, we analyzed a mouse model of BCR-ABL1+ pre–B ALL together with a new model of inducible expression of wild-type Ikaros in IKZF1 mutant human BCR-ABL1+ pre–B ALL. We performed integrated genome-wide chromatin and expression analyses and identified Ikaros target genes in mouse and human BCR-ABL1+ pre–B ALL, revealing novel conserved gene pathways associated with Ikaros tumor suppressor function. Notably, genetic depletion of different Ikaros targets, including CTNND1 and the early hematopoietic cell surface marker CD34, resulted in reduced leukemic growth. Our results suggest that Ikaros mediates tumor suppressor function by enforcing proper developmental stage–specific expression of multiple genes through chromatin compaction at its target genes.

Intranuclear and higher-order chromatin organization of the major histone gene cluster in breast cancer†

Alterations in nuclear morphology are common in cancer progression. However, the degree to which gross morphological abnormalities translate into compromised higher‐order chromatin organization is poorly understood. To explore the functional links between gene expression and chromatin structure in breast cancer, we performed RNA‐seq gene expression analysis on the basal breast cancer progression model based on human MCF10A cells. Positional gene enrichment identified the major histone gene cluster at chromosome 6p22 as one of the most significantly upregulated (and not amplified) clusters of genes from the normal‐like MCF10A to premalignant MCF10AT1 and metastatic MCF10CA1a cells. This cluster is subdivided into three sub‐clusters of histone genes that are organized into hierarchical topologically associating domains (TADs). Interestingly, the sub‐clusters of histone genes are located at TAD boundaries and interact more frequently with each other than the regions in‐between them, suggesting that the histone sub‐clusters form an active chromatin hub. The anchor sites of loops within this hub are occupied by CTCF, a known chromatin organizer. These histone genes are transcribed and processed at a specific sub‐nuclear microenvironment termed the major histone locus body (HLB). While the overall chromatin structure of the major HLB is maintained across breast cancer progression, we detected alterations in its structure that may relate to gene expression. Importantly, breast tumor specimens also exhibit a coordinate pattern of upregulation across the major histone gene cluster. Our results provide a novel insight into the connection between the higher‐order chromatin organization of the major HLB and its regulation during breast cancer progression.

Antagonism of B cell enhancer networks by STAT5 drives leukemia and poor patient survival

The transcription factor STAT5 has a critical role in B cell acute lymphoblastic leukemia (B-ALL). How STAT5 mediates this effect is unclear. Here we found that activation of STAT5 worked together with defects in signaling components of the precursor to the B cell antigen receptor (pre-BCR), including defects in BLNK, BTK, PKCβ, NF-κB1 and IKAROS, to initiate B-ALL. STAT5 antagonized the transcription factors NF-κB and IKAROS by opposing regulation of shared target genes. Super-enhancers showed enrichment for STAT5 binding and were associated with an opposing network of transcription factors, including PAX5, EBF1, PU.1, IRF4 and IKAROS. Patients with a high ratio of active STAT5 to NF-κB or IKAROS had more-aggressive disease. Our studies indicate that an imbalance of two opposing transcriptional programs drives B-ALL and suggest that restoring the balance of these pathways might inhibit B-ALL.

Conserved IKAROS-regulated genes associated with B-progenitor acute lymphoblastic leukemia outcome

Genetic alterations disrupting the transcription factor IKZF1 (encoding IKAROS) are associated with poor outcome in B lineage acute lymphoblastic leukemia (B-ALL) and occur in >70% of the high-risk BCR-ABL1+ (Ph+) and Ph-like disease subtypes. To examine IKAROS function in this context, we have developed novel mouse models allowing reversible RNAi-based control of Ikaros expression in established B-ALL in vivo. Notably, leukemias driven by combined BCR-ABL1 expression and Ikaros suppression rapidly regress when endogenous Ikaros is restored, causing sustained disease remission or ablation. Comparison of transcriptional profiles accompanying dynamic Ikaros perturbation in murine B-ALL in vivo with two independent human B-ALL cohorts identified nine evolutionarily conserved IKAROS-repressed genes. Notably, high expression of six of these genes is associated with inferior event–free survival in both patient cohorts. Among them are EMP1, which was recently implicated in B-ALL proliferation and prednisolone resistance, and the novel target CTNND1, encoding P120-catenin. We demonstrate that elevated Ctnnd1 expression contributes to maintenance of murine B-ALL cells with compromised Ikaros function. These results suggest that IKZF1 alterations in B-ALL leads to induction of multiple genes associated with proliferation and treatment resistance, identifying potential new therapeutic targets for high-risk disease.

Salivary microbiomes of indigenous Tsimane mothers and infants are distinct despite frequent premastication

Background Premastication, the transfer of pre-chewed food, is a common infant and young child feeding practice among the Tsimane, forager-horticulturalists living in the Bolivian Amazon. Research conducted primarily with Western populations has shown that infants harbor distinct oral microbiota from their mothers. Premastication, which is less common in these populations, may influence the colonization and maturation of infant oral microbiota, including via transmission of oral pathogens. We collected premasticated food and saliva samples from Tsimane mothers and infants (9–24 months of age) to test for evidence of bacterial transmission in premasticated foods and overlap in maternal and infant salivary microbiota. We extracted bacterial DNA from two premasticated food samples and 12 matched salivary samples from maternal-infant pairs. DNA sequencing was performed with MiSeq (Illumina). We evaluated maternal and infant microbial composition in terms of relative abundance of specific taxa, alpha and beta diversity, and dissimilarity distances. Results The bacteria in saliva and premasticated food were mapped to 19 phyla and 400 genera and were dominated by Firmicutes, Proteobacteria, Actinobacteria, and Bacteroidetes. The oral microbial communities of Tsimane mothers and infants who frequently share premasticated food were well-separated in a non-metric multi-dimensional scaling ordination (NMDS) plot. Infant microbiotas clustered together, with weighted Unifrac distances significantly differing between mothers and infants. Infant saliva contained more Firmicutes (p < 0.01) and fewer Proteobacteria (p < 0.05) than did maternal saliva. Many genera previously associated with dental and periodontal infections, e.g. Neisseria, Gemella, Rothia, Actinomyces, Fusobacterium, and Leptotrichia, were more abundant in mothers than in infants. Conclusions Salivary microbiota of Tsimane infants and young children up to two years of age do not appear closely related to those of their mothers, despite frequent premastication and preliminary evidence that maternal bacteria is transmitted to premasticated foods. Infant physiology and diet may constrain colonization by maternal bacteria, including several oral pathogens.

Bivalent Epigenetic Control of Oncofetal Gene Expression in Cancer

ABSTRACT Multiple mechanisms of epigenetic control that include DNA methylation, histone modification, noncoding RNAs, and mitotic gene bookmarking play pivotal roles in stringent gene regulation during lineage commitment and maintenance. Experimental evidence indicates that bivalent chromatin domains, i.e., genome regions that are marked by both H3K4me3 (activating) and H3K27me3 (repressive) histone modifications, are a key property of pluripotent stem cells. Bivalency of developmental genes during the G1 phase of the pluripotent stem cell cycle contributes to cell fate decisions. Recently, some cancer types have been shown to exhibit partial recapitulation of bivalent chromatin modifications that are lost along with pluripotency, suggesting a mechanism by which cancer cells reacquire properties that are characteristic of undifferentiated, multipotent cells. This bivalent epigenetic control of oncofetal gene expression in cancer cells may offer novel insights into the onset and progression of cancer and may provide specific and selective options for diagnosis as well as for therapeutic intervention.

Transcriptomic Signatures of Tacaribe Virus-Infected Jamaican Fruit Bats

As reservoir hosts of viruses associated with human disease, little is known about the interactions between bats and viruses. Using Jamaican fruit bats infected with Tacaribe virus (TCRV) as a model, we characterized the gene expression responses to infection in different tissues and identified pathways involved with the response to infection. This report is the most detailed gene discovery work in the species to date and the first to describe immune gene expression responses in bats during a pathogenic viral infection. ABSTRACT Tacaribe virus (TCRV) is a mammalian arenavirus that was first isolated from artibeus bats in the 1950s. Subsequent experimental infection of Jamaican fruit bats (Artibeus jamaicensis) caused a disease similar to that of naturally infected bats. Although substantial attention has focused on bats as reservoir hosts of viruses that cause human disease, little is known about the interactions between bats and their pathogens. We performed a transcriptome-wide study to illuminate the response of Jamaican fruit bats experimentally infected with TCRV. Differential gene expression analysis of multiple tissues revealed global and organ-specific responses associated with innate antiviral responses, including interferon alpha/beta and Toll-like receptor signaling, activation of complement cascades, and cytokine signaling, among others. Genes encoding proteins involved in adaptive immune responses, such as gamma interferon signaling and costimulation of T cells by the CD28 family, were also altered in response to TCRV infection. Immunoglobulin gene expression was also elevated in the spleens of infected bats, including IgG, IgA, and IgE isotypes. These results indicate an active innate and adaptive immune response to TCRV infection occurred but did not prevent fatal disease. This de novo assembly provides a high-throughput data set of the Jamaican fruit bat and its host response to TCRV infection, which remains a valuable tool to understand the molecular signatures involved in antiviral responses in bats. IMPORTANCE As reservoir hosts of viruses associated with human disease, little is known about the interactions between bats and viruses. Using Jamaican fruit bats infected with Tacaribe virus (TCRV) as a model, we characterized the gene expression responses to infection in different tissues and identified pathways involved with the response to infection. This report is the most detailed gene discovery work in the species to date and the first to describe immune gene expression responses in bats during a pathogenic viral infection.

Genome-wide analysis reveals positional-nucleosome-oriented binding pattern of pioneer factor FOXA1

The compaction of nucleosomal structures creates a barrier for DNA-binding transcription factors (TFs) to access their cognate cis-regulatory elements. Pioneer factors (PFs) such as FOXA1 are able to directly access these cis-targets within compact chromatin. However, how these PFs interplay with nucleosomes remains to be elucidated, and is critical for us to understand the underlying mechanism of gene regulation. Here, we have conducted a computational analysis on a strand-specific paired-end ChIP-exo (termed as ChIP-ePENS) data of FOXA1 in LNCaP cells by our novel algorithm ePEST. We find that FOXA1 chromatin binding occurs via four distinct border modes (or footprint boundary patterns), with a preferential footprint boundary patterns relative to FOXA1 motif orientation. In addition, from this analysis three fundamental nucleotide positions (oG, oS and oH) emerged as major determinants for blocking exo-digestion and forming these four distinct border modes. By integrating histone MNase-seq data, we found an astonishingly consistent, ‘well-positioned’ configuration occurs between FOXA1 motifs and dyads of nucleosomes genome-wide. We further performed ChIP-seq of eight chromatin remodelers and found an increased occupancy of these remodelers on FOXA1 motifs for all four border modes (or footprint boundary patterns), indicating the full occupancy of FOXA1 complex on the three blocking sites (oG, oS and oH) likely produces an active regulatory status with well-positioned phasing for protein binding events. Together, our results suggest a positional-nucleosome-oriented accessing model for PFs seeking target motifs, in which FOXA1 can examine each underlying DNA nucleotide and is able to sense all potential motifs regardless of whether they face inward or outward from histone octamers along the DNA helix axis.

Targeted Genome Sequencing Reveals Varicella-Zoster Virus Open Reading Frame 12 Deletion

ABSTRACT The neurotropic herpesvirus varicella-zoster virus (VZV) establishes a lifelong latent infection in humans following primary infection. The low abundance of VZV nucleic acids in human neurons has hindered an understanding of the mechanisms that regulate viral gene transcription during latency. To overcome this critical barrier, we optimized a targeted capture protocol to enrich VZV DNA and cDNA prior to whole-genome/transcriptome sequence analysis. Since the VZV genome is remarkably stable, it was surprising to detect that VZV32, a VZV laboratory strain with no discernible growth defect in tissue culture, contained a 2,158-bp deletion in open reading frame (ORF) 12. Consequently, ORF 12 and 13 protein expression was abolished and Akt phosphorylation was inhibited. The discovery of the ORF 12 deletion, revealed through targeted genome sequencing analysis, points to the need to authenticate the VZV genome when the virus is propagated in tissue culture. IMPORTANCE Viruses isolated from clinical samples often undergo genetic modifications when cultured in the laboratory. Historically, VZV is among the most genetically stable herpesviruses, a notion supported by more than 60 complete genome sequences from multiple isolates and following multiple in vitro passages. However, application of enrichment protocols to targeted genome sequencing revealed the unexpected deletion of a significant portion of VZV ORF 12 following propagation in cultured human fibroblast cells. While the enrichment protocol did not introduce bias in either the virus genome or transcriptome, the findings indicate the need for authentication of VZV by sequencing when the virus is propagated in tissue culture.

Varicella Zoster Virus Induces Nuclear Translocation of the Neurokinin-1 Receptor, Promoting Lamellipodia Formation and Viral Spread in Spinal Astrocytes

Background Varicella zoster virus (VZV) can present as a myelopathy with spinal astrocyte infection. Recent studies support a role for the neurokinin-1 receptor (NK-1R) in virus infections, as well as for cytoskeletal alterations that may promote viral spread. Thus, we examined the role of NK-1R in VZV-infected primary human spinal astrocytes (HA-sps) to shed light on the pathogenesis of VZV myelopathy. Methods Mock- and VZV-infected HA-sps were examined for substance P (subP) production, NK-1R localization, morphological changes, and viral spread in the presence or absence of the NK-1R antagonists aprepitant and rolapitant. Results VZV infection of HA-sps induced nuclear localization of full-length and truncated NK-1R in the absence of the endogenous ligand, subP, and was associated with extensive lamellipodia formation and viral spread that was inhibited by NK-1R antagonists. Conclusions We have identified a novel, subP-independent, proviral function of nuclear NK-1R associated with lamellipodia formation and viral spread that is distinct from subP-induced NK-1R cell membrane/cytoplasmic localization without lamellipodia formation. These results suggest that binding of a putative viral ligand to NK-1R produces a dramatically different NK-1R downstream effect than binding of subP. Finally, the NK-1R antagonists aprepitant and rolapitant provide promising alternatives to nucleoside analogs in treating VZV infections, including myelopathy.