Genome-Wide Role of HSF1 in Transcriptional Regulation of Desiccation Tolerance in the Anhydrobiotic Cell Line, Pv11.

The Pv11, an insect cell line established from the midge Polypedilum vanderplanki, is capable of extreme hypometabolic desiccation tolerance, so-called anhydrobiosis. We previously discovered that heat shock factor 1 (HSF1) contributes to the acquisition of desiccation tolerance by Pv11 cells, but the mechanistic details have yet to be elucidated. Here, by analyzing the gene expression profiles of newly established HSF1-knockout and -rescue cell lines, we show that HSF1 has a genome-wide effect on gene regulation in Pv11. The HSF1-knockout cells exhibit a reduced desiccation survival rate, but this is completely restored in HSF1-rescue cells. By comparing mRNA profiles of the two cell lines, we reveal that HSF1 induces anhydrobiosis-related genes, especially genes encoding late embryogenesis abundant proteins and thioredoxins, but represses a group of genes involved in basal cellular processes, thus promoting an extreme hypometabolism state in the cell. In addition, HSF1 binding motifs are enriched in the promoters of anhydrobiosis-related genes and we demonstrate binding of HSF1 to these promoters by ChIP-qPCR. Thus, HSF1 directly regulates the transcription of anhydrobiosis-related genes and consequently plays a pivotal role in the induction of anhydrobiotic ability in Pv11 cells.

Systematic analysis of transcription start sites in avian development.

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds.

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

The level of protein in the maternal murine diet modulates the facial appearance of the offspring via mTORC1 signaling.

The development of craniofacial skeletal structures is fascinatingly complex and elucidation of the underlying mechanisms will not only provide novel scientific insights, but also help develop more effective clinical approaches to the treatment and/or prevention of the numerous congenital craniofacial malformations. To this end, we performed a genome-wide analysis of RNA transcription from non-coding regulatory elements by CAGE-sequencing of the facial mesenchyme of human embryos and cross-checked the active enhancers thus identified against genes, identified by GWAS for the normal range human facial appearance. Among the identified active cis-enhancers, several belonged to the components of the PI3/AKT/mTORC1/autophagy pathway. To assess the functional role of this pathway, we manipulated it both genetically and pharmacologically in mice and zebrafish. These experiments revealed that mTORC1 signaling modulates craniofacial shaping at the stage of skeletal mesenchymal condensations, with subsequent fine-tuning during clonal intercalation. This ability of mTORC1 pathway to modulate facial shaping, along with its evolutionary conservation and ability to sense external stimuli, in particular dietary amino acids, indicate that the mTORC1 pathway may play a role in facial phenotypic plasticity. Indeed, the level of protein in the diet of pregnant female mice influenced the activity of mTORC1 in fetal craniofacial structures and altered the size of skeletogenic clones, thus exerting an impact on the local geometry and craniofacial shaping. Overall, our findings indicate that the mTORC1 signaling pathway is involved in the effect of environmental conditions on the shaping of craniofacial structures.

ZENBU-Reports: a graphical web-portal builder for interactive visualization and dissemination of genome-scale data.

In the genomic era, data dissemination and visualization is an integral part of scientific publications and research projects involving international consortia producing massive genome-wide data sets, intra-organizational collaborations, or individual labs. However, creating custom supporting websites is oftentimes impractical due to the required programming effort, web server infrastructure, and data storage facilities, as well as the long-term maintenance burden. ZENBU-Reports (https://fantom.gsc.riken.jp/zenbu/reports) is a web application to create interactive scientific web portals by using graphical interfaces while providing storage and secured collaborative sharing for data uploaded by users. ZENBU-Reports provides the scientific visualization elements commonly used in supplementary websites, publications and presentations, presenting a complete solution for the interactive display and dissemination of data and analysis results during the full lifespan of a scientific project both during the active research phase and after publication of the results.

Dosage compensation of Z sex chromosome genes in avian fibroblast cells.

In birds, sex is genetically determined; however, the molecular mechanism is not well-understood. The avian Z sex chromosome (chrZ) lacks whole chromosome inactivation, in contrast to the mammalian chrX. To investigate chrZ dosage compensation and its role in sex specification, we use a highly quantitative method and analyze transcriptional activities of male and female fibroblast cells from seven bird species. Our data indicate that three fourths of chrZ genes are strictly compensated across Aves, similar to mammalian chrX. We also present a complete list of non-compensated chrZ genes and identify Ribosomal Protein S6 (RPS6) as a conserved sex-dimorphic gene in birds.

Dissection of gastric homeostasis in vivo facilitates permanent capture of isthmus-like stem cells in vitro.

The glandular stomach is composed of two regenerative compartments termed corpus and antrum, and our understanding of the transcriptional networks that maintain these tissues is incomplete. Here we show that cell types with equivalent functional roles in the corpus and antrum share similar transcriptional states including the poorly characterized stem cells of the isthmus region. To further study the isthmus, we developed a monolayer two-dimensional (2D) culture system that is continually maintained by Wnt-responsive isthmus-like cells capable of differentiating into several gastric cell types. Importantly, 2D cultures can be converted into conventional three-dimensional organoids, modelling the plasticity of gastric epithelial cells in vivo. Finally, we utilized the 2D culture system to show that Sox2 is both necessary and sufficient to generate enterochromaffin cells. Together, our data provide important insights into gastric homeostasis, establish a tractable culture system to capture isthmus cells and uncover a role for Sox2 in enterochromaffin cells.

Unraveling the Mechanisms of Sensitivity to Anti-FGF Therapies in Imatinib-Resistant Gastrointestinal Stromal Tumors (GIST) Lacking Secondary KIT Mutations.

We showed previously that inhibition of KIT signaling in GISTs activates FGFR-signaling pathway rendering cancer cells resistant to receptor tyrosine kinase inhibitor (RTKi) imatinib mesylate (IM) (Gleevec) despite of absence of secondary KIT mutations and thereby illustrating a rationale for the combined (e.g., KIT- and FGFR-targeted) therapies. We show here that long-term culture of IM-resistant GISTs (GIST-R1) with IM substantially down-regulates KIT expression and induces activation of the FGFR-signaling cascade, evidenced by increased expression of total and phosphorylated forms of FGFR1 and 2, FGF-2, and FRS-2, the well-known adaptor protein of the FGF-signaling cascade. This resulted in activation of both AKT- and MAPK-signaling pathways shown on mRNA and protein levels, and rendered cancer cells highly sensitive to pan-FGFR-inhibitors (BGJ 398, AZD 4547, and TAS-120). Indeed, we observed a significant decrease of IC50 values for BGJ 398 in the GIST subclone (GIST-R2) derived from GIST-R1 cells continuously treated with IM for up to 12 months. An increased sensitivity of GIST-R2 cells to FGFR inhibition was also revealed on the xenograft models, illustrating a substantial (>70%) decrease in tumor size in BGJ 398-treated animals when treated with this pan-FGFR inhibitor. Similarly, an increased intra-tumoral apoptosis as detected by immunohistochemical (IHC)-staining for cleaved caspase-3 on day 5 of the treatment was found. As expected, both BGJ 398 and IM used alone lacked the pro-apoptotic and growth-inhibitory activities on GIST-R1 xenografts, thereby revealing their resistance to these TKis when used alone. Important, the knockdown of FGFR2, and, in much less content, FGF-2, abrogated BGJ 398's activity against GIST-R2 cells both in vitro and in vivo, thereby illustrating the FGF-2/FGFR2-signaling axis in IM-resistant GISTs as a primary molecular target for this RTKi. Collectively, our data illustrates that continuous inhibition of KIT signaling in IM-resistant GISTs lacking secondary KIT mutations induced clonal heterogeneity of GISTs and resulted in accumulation of cancer cells with overexpressed FGF-2 and FGFR1/2, thereby leading to activation of FGFR-signaling. This in turn rendered these cells extremely sensitive to the pan-FGFR inhibitors used in combination with IM, or even alone, and suggests a rationale to re-evaluate the effectiveness of FGFR-inhibitors in order to improve the second-line therapeutic strategies for selected subgroups of GIST patients (e.g., IM-resistant GISTs lacking secondary KIT mutations and exhibiting the activation of the FGFR-signaling pathway).

Population Genomics of Two Closely Related Anhydrobiotic Midges Reveals Differences in Adaptation to Extreme Desiccation.

The sleeping chironomid Polypedilum vanderplanki is capable of anhydrobiosis, a striking example of adaptation to extreme desiccation. Tolerance to complete desiccation in this species is associated with emergence of multiple paralogs of protective genes. One of the gene families highly expressed under anhydrobiosis and involved in this process is protein-L-isoaspartate (D-aspartate) O-methyltransferases (PIMTs). Recently, another closely related midge was discovered, Polypedilum pembai, which is able not only to tolerate desiccation but also to survive multiple desiccation-rehydration cycles. To investigate the evolution of anhydrobiosis in these species, we sequenced and assembled the genome of P. pembai and compared it with P. vanderplanki and also performed a population genomics analysis of several populations of P. vanderplanki and one population of P. pembai. We observe positive selection and radical changes in the genetic architecture of the PIMT locus between the two species, including its amplification in the P. pembai lineage. In particular, PIMT-4, the most highly expressed of these PIMTs, is present in six copies in the P. pembai; these copies differ in expression profiles, suggesting possible sub- or neofunctionalization. The nucleotide diversity of the genomic region carrying these new genes is decreased in P. pembai, but not in the orthologous region carrying the ancestral gene in P. vanderplanki, providing evidence for a selective sweep associated with postduplication adaptation in the former. Overall, our results suggest an extensive relatively recent and likely ongoing adaptation of the mechanisms of anhydrobiosis.

A previously uncharacterized Factor Associated with Metabolism and Energy (FAME/C14orf105/CCDC198/1700011H14Rik) is related to evolutionary adaptation, energy balance, and kidney physiology.

In this study we use comparative genomics to uncover a gene with uncharacterized function (1700011H14Rik/C14orf105/CCDC198), which we hereby name FAME (Factor Associated with Metabolism and Energy). We observe that FAME shows an unusually high evolutionary divergence in birds and mammals. Through the comparison of single nucleotide polymorphisms, we identify gene flow of FAME from Neandertals into modern humans. We conduct knockout experiments on animals and observe altered body weight and decreased energy expenditure in Fame knockout animals, corresponding to genome-wide association studies linking FAME with higher body mass index in humans. Gene expression and subcellular localization analyses reveal that FAME is a membrane-bound protein enriched in the kidneys. Although the gene knockout results in structurally normal kidneys, we detect higher albumin in urine and lowered ferritin in the blood. Through experimental validation, we confirm interactions between FAME and ferritin and show co-localization in vesicular and plasma membranes.

Identification of Genomic Safe Harbors in the Anhydrobiotic Cell Line, Pv11.

Genomic safe harbors (GSHs) provide ideal integration sites for generating transgenic organisms and cells and can be of great benefit in advancing the basic and applied biology of a particular species. Here we report the identification of GSHs in a dry-preservable insect cell line, Pv11, which derives from the sleeping chironomid, Polypedilum vanderplanki, and similar to the larvae of its progenitor species exhibits extreme desiccation tolerance. To identify GSHs, we carried out genome analysis of transgenic cell lines established by random integration of exogenous genes and found four candidate loci. Targeted knock-in was performed into these sites and the phenotypes of the resulting transgenic cell lines were examined. Precise integration was achieved for three candidate GSHs, and in all three cases integration did not alter the anhydrobiotic ability or the proliferation rate of the cell lines. We therefore suggest these genomic loci represent GSHs in Pv11 cells. Indeed, we successfully constructed a knock-in system and introduced an expression unit into one of these GSHs. We therefore identified several GSHs in Pv11 cells and developed a new technique for producing transgenic Pv11 cells without affecting the phenotype.

Characteristics and Resistance to Cisplatin of Human Neuroblastoma Cells Co-Cultivated with Immune and Stromal Cells.

We investigated the features of the morphology and cytokine profiles of neuroblastoma SH-SY5Y cells, bone marrow-derived mesenchymal stromal/stem cells (BM-MSCs), and peripheral blood mononuclear cells (PBMCs) in double (BM-MSCs + SH-SY5Y cells) and triple (BM-MSCs + SH-SY5Y cells + PBMCs) co-cultures incubated on plastic and Matrigel. Cells in the co-cultures communicated by vesicular transport and by exchanging membrane and cytoplasmic components. The cytokine profile of double and triple co-cultures incubated on Matrigel and plastic had differences and showed the highest concentration of a number of chemokines/cytokines, such as CXCL8/IL-8, I-TAC/CXCL11, IP10/CXCL10, MDC/CCL22, MIP-1α/CCL3, IL-1ß, ENA-78/CXCL5, Gro-α/CXCL1, MCP-1/CCL2, TERC/CCL25, CXCL8/IL-8, and IL-6. High concentrations of inflammatory chemokines/cytokines in the conditioned medium of triple co-culture form a chronic inflammation, which brings the presented co-cultivation system closer to a natural tumor. Triple co-cultures were more resistant to cisplatin (CDDP) than the double- and monoculture of SH-SY5Y. The mRNA levels of BCL2, BCL2L1, RAC1, CAV1, CASP3, and BAX genes were changed in cells after co-culturing and CDDP treatment in double and triple co-cultures. The expression of the BCL2, BAX, CAV1, and CASP3 proteins in SH-SY5Y cells after the triple co-culture and CAV1 and BAX protein expression in SH-SY5Y cells after the double co-culture were determined. This study demonstrated the nature of the cellular interactions between components of tumor niche and the intercellular influence on chemoresistance observed in our tumor model, which should enable the development of novel test systems for anti-tumor agents.

Upper cortical layer-driven network impairment in schizophrenia.

Schizophrenia is one of the most widespread and complex mental disorders. To characterize the impact of schizophrenia, we performed single-nucleus RNA sequencing (snRNA-seq) of >220,000 neurons from the dorsolateral prefrontal cortex of patients with schizophrenia and matched controls. In addition, >115,000 neurons were analyzed topographically by immunohistochemistry. Compositional analysis of snRNA-seq data revealed a reduction in abundance of GABAergic neurons and a concomitant increase in principal neurons, most pronounced for upper cortical layer subtypes, which was substantiated by histological analysis. Many neuronal subtypes showed extensive transcriptomic changes, the most marked in upper-layer GABAergic neurons, including down-regulation in energy metabolism and up-regulation in neurotransmission. Transcription factor network analysis demonstrated a developmental origin of transcriptomic changes. Last, Visium spatial transcriptomics further corroborated upper-layer neuron vulnerability in schizophrenia. Overall, our results point toward general network impairment within upper cortical layers as a core substrate associated with schizophrenia symptomatology.

High quality genome assembly of the anhydrobiotic midge provides insights on a single chromosome-based emergence of extreme desiccation tolerance.

Non-biting midges (Chironomidae) are known to inhabit a wide range of environments, and certain species can tolerate extreme conditions, where the rest of insects cannot survive. In particular, the sleeping chironomid Polypedilum vanderplanki is known for the remarkable ability of its larvae to withstand almost complete desiccation by entering a state called anhydrobiosis. Chromosome numbers in chironomids are higher than in other dipterans and this extra genomic resource might facilitate rapid adaptation to novel environments. We used improved sequencing strategies to assemble a chromosome-level genome sequence for P. vanderplanki for deep comparative analysis of genomic location of genes associated with desiccation tolerance. Using whole genome-based cross-species and intra-species analysis, we provide evidence for the unique functional specialization of Chromosome 4 through extensive acquisition of novel genes. In contrast to other insect genomes, in the sleeping chironomid a uniquely high degree of subfunctionalization in paralogous anhydrobiosis genes occurs in this chromosome, as well as pseudogenization in a highly duplicated gene family. Our findings suggest that the Chromosome 4 in Polypedilum is a site of high genetic turnover, allowing it to act as a 'sandbox' for evolutionary experiments, thus facilitating the rapid adaptation of midges to harsh environments.

Integrative transcription start site analysis and physiological phenotyping reveal torpor-specific expression program in mouse skeletal muscle.

Mice enter an active hypometabolic state, called daily torpor when they experience a lowered caloric intake under cold ambient temperature. During torpor, the oxygen consumption rate in some animals drops to less than 30% of the normal rate without harming the body. This safe but severe reduction in metabolism is attractive for various clinical applications; however, the mechanism and molecules involved are unclear. Therefore, here we systematically analyzed the gene expression landscape on the level of the RNA transcription start sites in mouse skeletal muscles under various metabolic states to identify torpor-specific transcribed regulatory patterns. We analyzed the soleus muscles from 38 mice in torpid and non-torpid conditions and identified 287 torpor-specific promoters out of 12,862 detected promoters. Furthermore, we found that the transcription factor ATF3 is highly expressed during torpor deprivation and its binding motif is enriched in torpor-specific promoters. Atf3 was also highly expressed in the heart and brown adipose tissue during torpor and systemically knocking out Atf3 affected the torpor phenotype. Our results demonstrate that mouse torpor combined with powerful genetic tools is useful for studying active hypometabolism.

Cas9-mediated genome editing reveals a significant contribution of calcium signaling pathways to anhydrobiosis in Pv11 cells.

Pv11 is an insect cell line established from the midge Polypedilum vanderplanki, whose larval form exhibits an extreme desiccation tolerance known as anhydrobiosis. Pv11 itself is also capable of anhydrobiosis, which is induced by trehalose treatment. Here we report the successful construction of a genome editing system for Pv11 cells and its application to the identification of signaling pathways involved in anhydrobiosis. Using the Cas9-mediated gene knock-in system, we established Pv11 cells that stably expressed GCaMP3 to monitor intracellular Ca2+ mobilization. Intriguingly, trehalose treatment evoked a transient increase in cytosolic Ca2+ concentration, and further experiments revealed that the calmodulin-calcineurin-NFAT pathway contributes to tolerance of trehalose treatment as well as desiccation tolerance, while the calmodulin-calmodulin kinase-CREB pathway conferred only desiccation tolerance on Pv11 cells. Thus, our results show a critical contribution of the trehalose-induced Ca2+ surge to anhydrobiosis and demonstrate temporally different roles for each signaling pathway.

Diversity and Regulation of S-Adenosylmethionine Dependent Methyltransferases in the Anhydrobiotic Midge.

Multiple co-localized paralogs of genes in Polypedilum vanderplanki's genome have strong transcriptional response to dehydration and considered to be a part of adaptation machinery at the larvae stage. One group of such genes represented by L-isoaspartate O-methyltransferases (PIMT). In order to highlight specific role of PIMT paralogization in desiccation tolerance of the larvae we annotated and compared S-adenosylmethionine (SAM) dependent methyltransferases of four insect species. From another side we applied co-expression analysis in desiccation/rehydration time course and showed that PIMT coding genes could be separated into five clusters by expression profile. We found that among Polypedilum vanderplanki's PIMTs only PIMT1 and PIMT2 have enzymatic activity in normal physiological conditions. From in silico analysis of the protein structures we found two highly variable regions outside of the active center, but also amino acid substitutions which may affect SAM stabilization. Overall, in this study we demonstrated features of Polypedilum vanderplanki's PIMT coding paralogs related to different roles in desiccation tolerance of the larvae. Our results also suggest a role of different SAM-methyltransferases in the adaptation, including GSMT, JHAMT, and candidates from other classes, which could be considered in future studies.

Correction: Genetic algorithm-based personalized models of human cardiac action potential.

[This corrects the article DOI: 10.1371/journal.pone.0231695.].

Genetic algorithm-based personalized models of human cardiac action potential.

We present a novel modification of genetic algorithm (GA) which determines personalized parameters of cardiomyocyte electrophysiology model based on set of experimental human action potential (AP) recorded at different heart rates. In order to find the steady state solution, the optimized algorithm performs simultaneous search in the parametric and slow variables spaces. We demonstrate that several GA modifications are required for effective convergence. Firstly, we used Cauchy mutation along a random direction in the parametric space. Secondly, relatively large number of elite organisms (6-10% of the population passed on to new generation) was required for effective convergence. Test runs with synthetic AP as input data indicate that algorithm error is low for high amplitude ionic currents (1.6±1.6% for IKr, 3.2±3.5% for IK1, 3.9±3.5% for INa, 8.2±6.3% for ICaL). Experimental signal-to-noise ratio above 28 dB was required for high quality GA performance. GA was validated against optical mapping recordings of human ventricular AP and mRNA expression profile of donor hearts. In particular, GA output parameters were rescaled proportionally to mRNA levels ratio between patients. We have demonstrated that mRNA-based models predict the AP waveform dependence on heart rate with high precision. The latter also provides a novel technique of model personalization that makes it possible to map gene expression profile to cardiac function.

New group of transmembrane proteins associated with desiccation tolerance in the anhydrobiotic midge Polypedilum vanderplanki.

Larvae of the sleeping chironomid Polypedilum vanderplanki are known for their extraordinary ability to survive complete desiccation in an ametabolic state called "anhydrobiosis". The unique feature of P. vanderplanki genome is the presence of expanded gene clusters associated with anhydrobiosis. While several such clusters represent orthologues of known genes, there is a distinct set of genes unique for P. vanderplanki. These include Lea-Island-Located (LIL) genes with no known orthologues except two of LEA genes of P. vanderplanki, PvLea1 and PvLea3. However, PvLIL proteins lack typical features of LEA such as the state of intrinsic disorder, hydrophilicity and characteristic LEA_4 motif. They possess four to five transmembrane domains each and we confirmed membrane targeting for three PvLILs. Conserved amino acids in PvLIL are located in transmembrane domains or nearby. PvLEA1 and PvLEA3 proteins are chimeras combining LEA-like parts and transmembrane domains, shared with PvLIL proteins. We have found that PvLil genes are highly upregulated during anhydrobiosis induction both in larvae of P. vanderplanki and P. vanderplanki-derived cultured cell line, Pv11. Thus, PvLil are a new intriguing group of genes that are likely to be associated with anhydrobiosis due to their common origin with some LEA genes and their induction during anhydrobiosis.