Protective neutralizing influenza antibody response in the absence of T follicular helper cells.

Virus infection induces the development of T follicular helper (TFH) and T helper 1 (TH1) cells. Although TFH cells are important in anti-viral humoral immunity, the contribution of TH1 cells to a protective antibody response remains unknown. We found that IgG2 antibodies predominated in the response to vaccination with inactivated influenza A virus (IAV) and were responsible for protective immunity to lethal challenge with pathogenic H5N1 and pandemic H1N1 IAV strains, even in mice that lacked TFH cells and germinal centers. The cytokines interleukin-21 and interferon-γ, which are secreted from TH1 cells, were essential for the observed greater persistence and higher titers of IgG2 protective antibodies. Our results suggest that TH1 induction could be a promising strategy for producing effective neutralizing antibodies against emerging influenza viruses.

The transcription factor ATF7 mediates lipopolysaccharide-induced epigenetic changes in macrophages involved in innate immunological memory.

Immunological memory is thought to be mediated exclusively by lymphocytes. However, enhanced innate immune responses caused by a previous infection increase protection against reinfection, which suggests the presence of innate immunological memory. Here we identified an important role for the stress-response transcription factor ATF7 in innate immunological memory. ATF7 suppressed a group of genes encoding factors involved in innate immunity in macrophages by recruiting the histone H3K9 dimethyltransferase G9a. Treatment with lipopolysaccharide, which mimics bacterial infection, induced phosphorylation of ATF7 via the kinase p38, which led to the release of ATF7 from chromatin and a decrease in repressive histone H3K9me2 marks. A partially disrupted chromatin structure and increased basal expression of target genes were maintained for long periods, which enhanced resistance to pathogens. ATF7 might therefore be important in controlling memory in cells of the innate immune system.

Clonal Origin and Lineage Ambiguity in Mixed Neuroendocrine Carcinoma of the Uterine Cervix.

Small-cell neuroendocrine carcinoma (SCNEC) of the cervix is a rare disease characterized by a high incidence of mixed tumors with other types of cancer. The mechanism underlying this mixed phenotype is not well understood. This study established a panel of organoid lines from patients with SCNEC of the cervix and ultimately focused on one line, which retained a mixed tumor phenotype, both in vitro and in vivo. Histologically, both organoids and xenograft tumors showed distinct differentiation into either SCNEC or adenocarcinoma in some regions and ambiguous differentiation in others. Tracking single cells indicated the existence of cells with bipotential differentiation toward SCNEC and adenocarcinomas. Single-cell transcriptional analysis identified three distinct clusters: SCNEC-like, adenocarcinoma-like, and a cluster lacking specific differentiation markers. The expression of neuroendocrine markers was enriched in the SCNEC-like cluster but not exclusively. Human papillomavirus 18 E6 was enriched in the SCNEC-like cluster, which showed higher proliferation and lower levels of the p53 pathway. After treatment with anticancer drugs, the expression of adenocarcinoma markers increased, whereas that of SCNEC decreased. Using a reporter system for keratin 19 expression, changes in the differentiation of each cell were shown to be associated with the shift in differentiation induced by drug treatment. These data suggest that mixed SCNEC/cervical tumors have a clonal origin and are characterized by an ambiguous and flexible differentiation state.

Ligand-specific c-Fos expression emerges from the spatiotemporal control of ErbB network dynamics.

Activation of ErbB receptors by epidermal growth factor (EGF) or heregulin (HRG) determines distinct cell-fate decisions, although signals propagate through shared pathways. Using mathematical modeling and experimental approaches, we unravel how HRG and EGF generate distinct, all-or-none responses of the phosphorylated transcription factor c-Fos. In the cytosol, EGF induces transient and HRG induces sustained ERK activation. In the nucleus, however, ERK activity and c-fos mRNA expression are transient for both ligands. Knockdown of dual-specificity phosphatases extends HRG-stimulated nuclear ERK activation, but not c-fos mRNA expression, implying the existence of a HRG-induced repressor of c-fos transcription. Further experiments confirmed that this repressor is mainly induced by HRG, but not EGF, and requires new protein synthesis. We show how a spatially distributed, signaling-transcription cascade robustly discriminates between transient and sustained ERK activities at the c-Fos system level. The proposed control mechanisms are general and operate in different cell types, stimulated by various ligands.

Hepatic posttranscriptional network comprised of CCR4-NOT deadenylase and FGF21 maintains systemic metabolic homeostasis.

Whole-body metabolic homeostasis is tightly controlled by hormone-like factors with systemic or paracrine effects that are derived from nonendocrine organs, including adipose tissue (adipokines) and liver (hepatokines). Fibroblast growth factor 21 (FGF21) is a hormone-like protein, which is emerging as a major regulator of whole-body metabolism and has therapeutic potential for treating metabolic syndrome. However, the mechanisms that control FGF21 levels are not fully understood. Herein, we demonstrate that FGF21 production in the liver is regulated via a posttranscriptional network consisting of the CCR4-NOT deadenylase complex and RNA-binding protein tristetraprolin (TTP). In response to nutrient uptake, CCR4-NOT cooperates with TTP to degrade AU-rich mRNAs that encode pivotal metabolic regulators, including FGF21. Disruption of CCR4-NOT activity in the liver, by deletion of the catalytic subunit CNOT6L, increases serum FGF21 levels, which ameliorates diet-induced metabolic disorders and enhances energy expenditure without disrupting bone homeostasis. Taken together, our study describes a hepatic CCR4-NOT/FGF21 axis as a hitherto unrecognized systemic regulator of metabolism and suggests that hepatic CCR4-NOT may serve as a target for devising therapeutic strategies in metabolic syndrome and related morbidities.

Tamoxifen mechanically reprograms the tumor microenvironment via HIF-1A and reduces cancer cell survival.

The tumor microenvironment is fundamental to cancer progression, and the influence of its mechanical properties is increasingly being appreciated. Tamoxifen has been used for many years to treat estrogen-positive breast cancer. Here we report that tamoxifen regulates the level and activity of collagen cross-linking and degradative enzymes, and hence the organization of the extracellular matrix, via a mechanism involving both the G protein-coupled estrogen receptor (GPER) and hypoxia-inducible factor-1 alpha (HIF-1A). We show that tamoxifen reduces HIF-1A levels by suppressing myosin-dependent contractility and matrix stiffness mechanosensing. Tamoxifen also downregulates hypoxia-regulated genes and increases vascularization in PDAC tissues. Our findings implicate the GPER/HIF-1A axis as a master regulator of peri-tumoral stromal remodeling and the fibrovascular tumor microenvironment and offer a paradigm shift for tamoxifen from a well-established drug in breast cancer hormonal therapy to an alternative candidate for stromal targeting strategies in PDAC and possibly other cancers.

Translation-dependent unwinding of stem-loops by UPF1 licenses Regnase-1 to degrade inflammatory mRNAs.

Regnase-1-mediated mRNA decay (RMD), in which inflammatory mRNAs harboring specific stem-loop structures are degraded, is a critical part of proper immune homeostasis. Prior to initial translation, Regnase-1 associates with target stem-loops but does not carry out endoribonucleolytic cleavage. Single molecule imaging revealed that UPF1 is required to first unwind the stem-loops, thus licensing Regnase-1 to proceed with RNA degradation. Following translation, Regnase-1 physically associates with UPF1 using two distinct points of interaction: The Regnase-1 RNase domain binds to SMG1-phosphorylated residue T28 in UPF1; in addition, an intrinsically disordered segment in Regnase-1 binds to the UPF1 RecA domain, enhancing the helicase activity of UPF1. The SMG1-UPF1-Regnase-1 axis targets pioneer rounds of translation and is critical for rapid resolution of inflammation through restriction of the number of proteins translated by a given mRNA. Furthermore, small-molecule inhibition of SMG1 prevents RNA unwinding in dendritic cells, allowing post-transcriptional control of innate immune responses.

Transcriptionally inducible Pleckstrin homology-like domain, family A, member 1, attenuates ErbB receptor activity by inhibiting receptor oligomerization.

Feedback control is a key mechanism in signal transduction, intimately involved in regulating the outcome of the cellular response. Here, we report a novel mechanism by which PHLDA1, Pleckstrin homology-like domain, family A, member 1, negatively regulates ErbB receptor signaling by inhibition of receptor oligomerization. We have found that the ErbB3 ligand, heregulin, induces PHILDA1 expression in MCF-7 cells. Transcriptionally-induced PHLDA1 protein directly binds to ErbB3, whereas knockdown of PHLDA1 increases complex formation between ErbB3 and ErbB2. To provide insight into the mechanism for our time-course and single-cell experimental observations, we performed a systematic computational search of network topologies of the mathematical models based on receptor dimer-tetramer formation in the ErbB activation processes. Our results indicate that only a model in which PHLDA1 inhibits formation of both dimers and tetramer can explain the experimental data. Predictions made from this model were further validated by single-molecule imaging experiments. Our studies suggest a unique regulatory feature of PHLDA1 to inhibit the ErbB receptor oligomerization process and thereby control the activity of receptor signaling network.

ASBEL-TCF3 complex is required for the tumorigenicity of colorectal cancer cells.

Wnt/ß-catenin signaling plays a key role in the tumorigenicity of colon cancer. Furthermore, it has been reported that lncRNAs are dysregulated in several steps of cancer development. Here we show that ß-catenin directly activates the transcription of the long noncoding RNA (lncRNA) ASBEL [antisense ncRNA in the ANA (Abundant in neuroepithelium area)/BTG3 (B-cell translocation gene 3) locus] and transcription factor 3 (TCF3), both of which are required for the survival and tumorigenicity of colorectal cancer cells. ASBEL interacts with and recruits TCF3 to the activating transcription factor 3 (ATF3) locus, where it represses the expression of ATF3. Furthermore, we demonstrate that ASBEL-TCF3-mediated down-regulation of ATF3 expression is required for the proliferation and tumorigenicity of colon tumor cells. ATF3, in turn, represses the expression of ASBEL Our results reveal a pathway involving an lncRNA and two transcription factors that plays a key role in Wnt/ß-catenin-mediated tumorigenesis. These results may provide insights into the variety of biological and pathological processes regulated by Wnt/ß-catenin signaling.

Long noncoding RNA UPAT promotes colon tumorigenesis by inhibiting degradation of UHRF1.

Many long noncoding RNAs (lncRNAs) are reported to be dysregulated in human cancers and play critical roles in tumor development and progression. Furthermore, it has been reported that many lncRNAs regulate gene expression by recruiting chromatin remodeling complexes to specific genomic loci or by controlling transcriptional or posttranscriptional processes. Here we show that an lncRNA termed UPAT [ubiquitin-like plant homeodomain (PHD) and really interesting new gene (RING) finger domain-containing protein 1 (UHRF1) Protein Associated Transcript] is required for the survival and tumorigenicity of colorectal cancer cells. UPAT interacts with and stabilizes the epigenetic factor UHRF1 by interfering with its ß-transducin repeat-containing protein (TrCP)-mediated ubiquitination. Furthermore, we demonstrate that UHRF1 up-regulates Stearoyl-CoA desaturase 1 and Sprouty 4, which are required for the survival of colon tumor cells. Our study provides evidence for an lncRNA that regulates protein ubiquitination and degradation and thereby plays a critical role in the survival and tumorigenicity of tumor cells. Our results suggest that UPAT and UHRF1 may be promising molecular targets for the therapy of colon cancer.

Genome-scale regression analysis reveals a linear relationship for promoters and enhancers after combinatorial drug treatment.

MOTIVATION: Drug combination therapy for treatment of cancers and other multifactorial diseases has the potential of increasing the therapeutic effect, while reducing the likelihood of drug resistance. In order to reduce time and cost spent in comprehensive screens, methods are needed which can model additive effects of possible drug combinations. RESULTS: We here show that the transcriptional response to combinatorial drug treatment at promoters, as measured by single molecule CAGE technology, is accurately described by a linear combination of the responses of the individual drugs at a genome wide scale. We also find that the same linear relationship holds for transcription at enhancer elements. We conclude that the described approach is promising for eliciting the transcriptional response to multidrug treatment at promoters and enhancers in an unbiased genome wide way, which may minimize the need for exhaustive combinatorial screens. AVAILABILITY AND IMPLEMENTATION: The CAGE sequence data used in this study is available in the DDBJ Sequence Read Archive (http://trace.ddbj.nig.ac.jp/index_e.html), accession number DRP001113. CONTACT: xin.gao@kaust.edu.sa or erik.arner@riken.jp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Base-pairing probability in the microRNA stem region affects the binding and editing specificity of human A-to-I editing enzymes ADAR1-p110 and ADAR2.

Adenosine deaminases acting on RNA (ADARs) catalyze the deamination of adenosine (A) to inosine (I). A-to-I RNA editing targets double-stranded RNA (dsRNA), and increases the complexity of gene regulation by modulating base pairing-dependent processes such as splicing, translation, and microRNA (miRNA)-mediated gene silencing. This study investigates the genome-wide binding preferences of the nuclear constitutive isoforms ADAR1-p110 and ADAR2 on human miRNA species by RNA immunoprecipitation of ADAR-bound small RNAs (RIP-seq). Our results suggest that secondary structure predicted by base-pairing probability in the mainly double-stranded region of a pre-miRNA or mature miRNA duplex may determine ADAR isoform preference for binding distinct subpopulations of miRNAs. Furthermore, we identify 31 unique editing sites with statistical significance, 19 sites of which are novel editing sites. Editing sites are enriched in the seed region responsible for target recognition by miRNAs, and isoform-specific nucleotide motifs in the immediate vicinity and opposite of editing sites are consistent with previous studies, and further reveal that ADAR2 may edit A/C bulges more frequently than ADAR1-p110 in the context of miRNA.

Mathematical modeling of atopic dermatitis reveals "double-switch" mechanisms underlying 4 common disease phenotypes.

BACKGROUND: The skin barrier acts as the first line of defense against constant exposure to biological, microbial, physical, and chemical environmental stressors. Dynamic interplay between defects in the skin barrier, dysfunctional immune responses, and environmental stressors are major factors in the development of atopic dermatitis (AD). A systems biology modeling approach can yield significant insights into these complex and dynamic processes through integration of prior biological data. OBJECTIVE: We sought to develop a multiscale mathematical model of AD pathogenesis that describes the dynamic interplay between the skin barrier, environmental stress, and immune dysregulation and use it to achieve a coherent mechanistic understanding of the onset, progression, and prevention of AD. METHODS: We mathematically investigated synergistic effects of known genetic and environmental risk factors on the dynamic onset and progression of the AD phenotype, from a mostly asymptomatic mild phenotype to a severe treatment-resistant form. RESULTS: Our model analysis identified a "double switch," with 2 concatenated bistable switches, as a key network motif that dictates AD pathogenesis: the first switch is responsible for the reversible onset of inflammation, and the second switch is triggered by long-lasting or frequent activation of the first switch, causing irreversible onset of systemic TH2 sensitization and worsening of AD symptoms. CONCLUSIONS: Our mathematical analysis of the bistable switch predicts that genetic risk factors decrease the threshold of environmental stressors to trigger systemic TH2 sensitization. This analysis predicts and explains 4 common clinical AD phenotypes from a mild and reversible phenotype through to severe and recalcitrant disease and provides a mechanistic explanation for clinically demonstrated preventive effects of emollient treatments against development of AD.

Fluctuation of Rac1 activity is associated with the phenotypic and transcriptional heterogeneity of glioma cells.

Phenotypic heterogeneity of cancer cells is caused not only by genetic and epigenetic alterations but also by stochastic variation of intracellular signaling molecules. Using cells that stably express Förster resonance energy transfer (FRET) biosensors, we show here a correlation between a temporal fluctuation in the activity of Rac1 and the invasive properties of C6 glioma cells. By using long-term time-lapse imaging, we found that Rac1 activity in C6 glioma cells fluctuated over a timescale that was substantially longer than that of the replication cycle. Because the relative level of Rac1 activity in each cell was unaffected by a suspension-adhesion procedure, we were able to sort C6 glioma cells according to the levels of Rac1 activity, yielding Rac1(high) and Rac1(low) cells. The Rac1(high) cells invaded more efficiently than did Rac1(low) cells in a Matrigel invasion assay. We assessed the transcriptional profiles of Rac1(high) and Rac1(low) cells and performed gene ontology analysis. Among the 14 genes that were most associated with the term 'membrane' (membrane-related genes) in Rac1(high) cells, we identified four genes that were associated with glioma invasion and Rac1 activity by using siRNA knockdown experiments. Among the transcription factors upregulated in Rac1(high) cells, Egr2 was found to positively regulate expression of the four membrane-related invasion-associated genes. The identified signaling network might cause the fluctuations in Rac1 activity and the heterogeneity in the invasive capacity of glioma cells.

Transcriptional dynamics reveal critical roles for non-coding RNAs in the immediate-early response.

The immediate-early response mediates cell fate in response to a variety of extracellular stimuli and is dysregulated in many cancers. However, the specificity of the response across stimuli and cell types, and the roles of non-coding RNAs are not well understood. Using a large collection of densely-sampled time series expression data we have examined the induction of the immediate-early response in unparalleled detail, across cell types and stimuli. We exploit cap analysis of gene expression (CAGE) time series datasets to directly measure promoter activities over time. Using a novel analysis method for time series data we identify transcripts with expression patterns that closely resemble the dynamics of known immediate-early genes (IEGs) and this enables a comprehensive comparative study of these genes and their chromatin state. Surprisingly, these data suggest that the earliest transcriptional responses often involve promoters generating non-coding RNAs, many of which are produced in advance of canonical protein-coding IEGs. IEGs are known to be capable of induction without de novo protein synthesis. Consistent with this, we find that the response of both protein-coding and non-coding RNA IEGs can be explained by their transcriptionally poised, permissive chromatin state prior to stimulation. We also explore the function of non-coding RNAs in the attenuation of the immediate early response in a small RNA sequencing dataset matched to the CAGE data: We identify a novel set of microRNAs responsible for the attenuation of the IEG response in an estrogen receptor positive cancer cell line. Our computational statistical method is well suited to meta-analyses as there is no requirement for transcripts to pass thresholds for significant differential expression between time points, and it is agnostic to the number of time points per dataset.

Statistical significance of combinatorial regulations.

More than three transcription factors often work together to enable cells to respond to various signals. The detection of combinatorial regulation by multiple transcription factors, however, is not only computationally nontrivial but also extremely unlikely because of multiple testing correction. The exponential growth in the number of tests forces us to set a strict limit on the maximum arity. Here, we propose an efficient branch-and-bound algorithm called the "limitless arity multiple-testing procedure" (LAMP) to count the exact number of testable combinations and calibrate the Bonferroni factor to the smallest possible value. LAMP lists significant combinations without any limit, whereas the family-wise error rate is rigorously controlled under the threshold. In the human breast cancer transcriptome, LAMP discovered statistically significant combinations of as many as eight binding motifs. This method may contribute to uncover pathways regulated in a coordinated fashion and find hidden associations in heterogeneous data.

Dynamically varying interactions between heregulin and ErbB proteins detected by single-molecule analysis in living cells.

Heregulin (HRG) belongs to the family of EGFs and activates the receptor proteins ErbB3 and ErbB4 in a variety of cell types to regulate cell fate. The interactions between HRG and ErbB3/B4 are important to the pathological mechanisms underlying schizophrenia and some cancers. Here, we observed the reaction kinetics between fluorescently labeled single HRG molecules and ErbB3/B4 on the surfaces of MCF-7 human breast cancer cells. The equilibrium association and the dissociation from equilibrium were also measured using single-molecule imaging techniques. The unitary association processes mirrored the EGF and ErbB1 interactions in HeLa cells [Teramura Y, et al. (2006) EMBO J 25:4215-4222], suggesting that the predimerization of the receptors, followed by intermediate formation (between the first and second ligand-binding events to a receptor dimer), accelerated the formation of doubly liganded signaling dimers of the receptor molecules. However, the dissociation analysis suggested that the first HRG dissociation from the doubly liganded dimer was rapid, but the second dissociation from the singly liganded dimer was slow. The dissociation rate constant from the liganded monomer was intermediate. The dynamic changes in the association and dissociation kinetics in relation to the dimerization of ErbB displayed negative cooperativity, which resulted in apparent low- and high-affinity sites of HRG association on the cell surface.

Individual hematopoietic stem cells in human bone marrow of patients with aplastic anemia or myelodysplastic syndrome stably give rise to limited cell lineages.

Mutation of the phosphatidylinositol N-acetylglucosaminyltransferase subunit A (PIG-A) gene in hematopoietic stem cells (HSCs) results in the loss of glycosylphosphatidylinositol-anchored proteins (GPI-APs) on HSCs, but minimally affects their development, and thus can be used as a clonal maker of HSCs. We analyzed GPI-APs expression on six major lineage cells in a total of 574 patients with bone marrow (BM) failure in which microenvironment itself is thought to be unaffected, including aplastic anemia (AA) or myelodysplastic syndrome (MDS). GPI-APs-deficient (GPI-APs(-) ) cells were detected in 250 patients. Whereas the GPI-APs(-) cells were seen in all six lineages in a majority of patients who had higher proportion ([dbmtequ]3%) of GPI-APs(-) cells, they were detected in only limited lineages in 92.9% of cases in the lower proportion (<3%) group. In all 250 cases, the same lineages of GPI-APs(-) cells were detected even after 6-18-month intervals, indicating that the GPI-APs(-) cells reflect hematopoiesis maintained by a self-renewing HSC in most of cases. The frequency of clones with limited lineages seen in mild cases of AA was similar to that in severe cases, and clones with limited lineages were seen even in two health volunteer cases. These results strongly suggest most individual HSCs produce only restricted lineages even in a steady state. While this restriction could reflect heterogeneity in the developmental potential of HSCs, we propose an alternative model in which the BM microenvironment is mosaic in supporting commitment of progenitors toward distinct lineages. Our computer simulation based on this model successfully recapitulated the observed clinical data.

LncRNA ZNNT1 induces p53 degradation by interfering with the interaction between p53 and the SART3-USP15 complex.

Mammalian genomes encode large number of long noncoding RNAs (lncRNAs) that play key roles in various biological processes, including proliferation, differentiation, and stem cell pluripotency. Recent studies have addressed that some lncRNAs are dysregulated in human cancers and may play crucial roles in tumor development and progression. Here, we show that the lncRNA ZNNT1 is required for the proliferation and tumorigenicity of colon cancer cells with wild-type p53. ZNNT1 knockdown leads to decreased ubiquitination and stabilization of p53 protein. Moreover, we demonstrate that ZNNT1 needs to interact with SART3 to destabilize p53 and to promote the proliferation and tumorigenicity of colon cancer cells. We further show that SART3 is associated with the ubiquitin-specific peptidase USP15 and that ZNNT1 may induce p53 destabilization by inhibiting this interaction. These results suggest that ZNNT1 interferes with the SART3-USP15 complex-mediated stabilization of p53 protein and thereby plays important roles in the proliferation and tumorigenicity of colon cancer cells. Our findings suggest that ZNNT1 may be a promising molecular target for the therapy of colon cancer.

Integrated quantitative analysis of the phosphoproteome and transcriptome in tamoxifen-resistant breast cancer.

Quantitative phosphoproteome and transcriptome analysis of ligand-stimulated MCF-7 human breast cancer cells was performed to understand the mechanisms of tamoxifen resistance at a system level. Phosphoproteome data revealed that WT cells were more enriched with phospho-proteins than tamoxifen-resistant cells after stimulation with ligands. Surprisingly, decreased phosphorylation after ligand perturbation was more common than increased phosphorylation. In particular, 17ß-estradiol induced down-regulation in WT cells at a very high rate. 17ß-Estradiol and the ErbB ligand heregulin induced almost equal numbers of up-regulated phospho-proteins in WT cells. Pathway and motif activity analyses using transcriptome data additionally suggested that deregulated activation of GSK3ß (glycogen-synthase kinase 3ß) and MAPK1/3 signaling might be associated with altered activation of cAMP-responsive element-binding protein and AP-1 transcription factors in tamoxifen-resistant cells, and this hypothesis was validated by reporter assays. An examination of clinical samples revealed that inhibitory phosphorylation of GSK3ß at serine 9 was significantly lower in tamoxifen-treated breast cancer patients that eventually had relapses, implying that activation of GSK3ß may be associated with the tamoxifen-resistant phenotype. Thus, the combined phosphoproteome and transcriptome data set analyses revealed distinct signal transcription programs in tumor cells and provided a novel molecular target to understand tamoxifen resistance.