An immune cell atlas reveals the dynamics of human macrophage specification during prenatal development.

Macrophages are heterogeneous and play critical roles in development and disease, but their diversity, function, and specification remain inadequately understood during human development. We generated a single-cell RNA sequencing map of the dynamics of human macrophage specification from PCW 4-26 across 19 tissues. We identified a microglia-like population and a proangiogenic population in 15 macrophage subtypes. Microglia-like cells, molecularly and morphologically similar to microglia in the CNS, are present in the fetal epidermis, testicle, and heart. They are the major immune population in the early epidermis, exhibit a polarized distribution along the dorsal-lateral-ventral axis, and interact with neural crest cells, modulating their differentiation along the melanocyte lineage. Through spatial and differentiation trajectory analysis, we also showed that proangiogenic macrophages are perivascular across fetal organs and likely yolk-sac-derived as microglia. Our study provides a comprehensive map of the heterogeneity and developmental dynamics of human macrophages and unravels their diverse functions during development.

Accurate transcriptome-wide identification and quantification of alternative polyadenylation from RNA-seq data with APAIQ.

Alternative polyadenylation (APA) enables a gene to generate multiple transcripts with different 3' ends, which is dynamic across different cell types or conditions. Many computational methods have been developed to characterize sample-specific APA using the corresponding RNA-seq data, but suffered from high error rate on both polyadenylation site (PAS) identification and quantification of PAS usage (PAU), and bias toward 3' untranslated regions. Here we developed a tool for APA identification and quantification (APAIQ) from RNA-seq data, which can accurately identify PAS and quantify PAU in a transcriptome-wide manner. Using 3' end-seq data as the benchmark, we showed that APAIQ outperforms current methods on PAS identification and PAU quantification, including DaPars2, Aptardi, mountainClimber, SANPolyA, and QAPA. Finally, applying APAIQ on 421 RNA-seq samples from liver cancer patients, we identified >540 tumor-associated APA events and experimentally validated two intronic polyadenylation candidates, demonstrating its capacity to unveil cancer-related APA with a large-scale RNA-seq data set.

CRISPR-iPAS: a novel dCAS13-based method for alternative polyadenylation interference.

Alternative polyadenylation (APA) plays an important role in gene regulation. With the recent application of novel sequencing technology in APA profiling, an ever-increasing number of APA genes/sites have been identified. However, the phenotypic relevance of most of these APA isoforms remains elusive, which is largely due to the lack of a convenient genetics tool for APA interference. To address this issue, herein, an efficient method is developed based on the CRISPR-dCas13 system, termed as CRISPR-iPAS. Out of eight different dCas13 proteins, Porphyromonas gulae (Pgu) dCas13b, is identified as the most effective one in blocking the usage of the polyadenylation site (PAS). With guide RNAs targeting at core regulatory elements, dPguCas13b enabled APA regulation of endogenous genes with different APA types, including tandem 3'UTR, alternative terminal exon, as well as intronic PAS. Finally, we demonstrated that the proposed APA perturbation tool could be used to investigate the functional relevance of APA isoforms.

Mammalian splicing divergence is shaped by drift, buffering in trans, and a scaling law.

Alternative splicing is ubiquitous, but the mechanisms underlying its pattern of evolutionary divergence across mammalian tissues are still underexplored. Here, we investigated the cis-regulatory divergences and their relationship with tissue-dependent trans-regulation in multiple tissues of an F1 hybrid between two mouse species. Large splicing changes between tissues are highly conserved and likely reflect functional tissue-dependent regulation. In particular, micro-exons frequently exhibit this pattern with high inclusion levels in the brain. Cis-divergence of splicing appears to be largely non-adaptive. Although divergence is in general associated with higher densities of sequence variants in regulatory regions, events with high usage of the dominant isoform apparently tolerate more mutations, explaining why their exon sequences are highly conserved but their intronic splicing site flanking regions are not. Moreover, we demonstrate that non-adaptive mutations are often masked in tissues where accurate splicing likely is more important, and experimentally attribute such buffering effect to trans-regulatory splicing efficiency.

QAUST: Protein Function Prediction Using Structure Similarity, Protein Interaction, and Functional Motifs.

The number of available protein sequences in public databases is increasing exponentially. However, a significant percentage of these sequences lack functional annotation, which is essential for the understanding of how biological systems operate. Here, we propose a novel method, Quantitative Annotation of Unknown STructure (QAUST), to infer protein functions, specifically Gene Ontology (GO) terms and Enzyme Commission (EC) numbers. QAUST uses three sources of information: structure information encoded by global and local structure similarity search, biological network information inferred by protein-protein interaction data, and sequence information extracted from functionally discriminative sequence motifs. These three pieces of information are combined by consensus averaging to make the final prediction. Our approach has been tested on 500 protein targets from the Critical Assessment of Functional Annotation (CAFA) benchmark set. The results show that our method provides accurate functional annotation and outperforms other prediction methods based on sequence similarity search or threading. We further demonstrate that a previously unknown function of human tripartite motif-containing 22 (TRIM22) protein predicted by QAUST can be experimentally validated.

Reciprocal signaling between mTORC1 and MNK2 controls cell growth and oncogenesis.

eIF4E plays key roles in protein synthesis and tumorigenesis. It is phosphorylated by the kinases MNK1 and MNK2. Binding of MNKs to eIF4G enhances their ability to phosphorylate eIF4E. Here, we show that mTORC1, a key regulator of mRNA translation and oncogenesis, directly phosphorylates MNK2 on Ser74. This suppresses MNK2 activity and impairs binding of MNK2 to eIF4G. These effects provide a novel mechanism by which mTORC1 signaling impairs the function of MNK2 and thereby decreases eIF4E phosphorylation. MNK2[S74A] knock-in cells show enhanced phosphorylation of eIF4E and S6K1 (i.e., increased mTORC1 signaling), enlarged cell size, and increased invasive and transformative capacities. MNK2[Ser74] phosphorylation was inversely correlated with disease progression in human prostate tumors. MNK inhibition exerted anti-proliferative effects in prostate cancer cells in vitro. These findings define a novel feedback loop whereby mTORC1 represses MNK2 activity and oncogenic signaling through eIF4E phosphorylation, allowing reciprocal regulation of these two oncogenic pathways.

Integrative multi-omics analysis of a colon cancer cell line with heterogeneous Wnt activity revealed RUNX2 as an epigenetic regulator of EMT.

Epithelial-mesenchymal transition (EMT) program, which facilitates tumor metastasis, stemness and therapy resistance, is a reversible biological process that is largely orchestrated at the epigenetic level under the regulation of different cell signaling pathways. EMT state is often heterogeneous within individual tumors, though the epigenetic drivers underlying such heterogeneity remain elusive. In colon cancer, hyperactivation of the Wnt/ß-catenin signaling not only drives tumor initiation, but also promotes metastasis in late stage by promoting EMT program. However, it is unknown whether the intratumorally heterogeneous Wnt activity could directly drive EMT heterogeneity, and, if so, what are the underlying epigenetic driver(s). Here, by analyzing a phenotypically and molecularly heterogeneous colon cancer cell line using single-cell RNA sequencing, we identified two distinct cell populations with positively correlated Wnt activity and EMT state. Integrative multi-omics analysis of these two cell populations revealed RUNX2 as a critical transcription factor epigenetically driving the EMT heterogeneity. Both in vitro and in vivo genetic perturbation assays validated the EMT-enhancing effect of RUNX2, which remodeled chromatin landscape and activated a panel of EMT-associated genes through binding to their promoters and/or potential enhancers. Finally, by exploring the clinical data, we showed that RUNX2 expression is positively correlated with metastasis development and poor survival of colon cancer patients, as well as patients afflicted with other types of cancer. Taken together, our work revealed RUNX2 as a new EMT-promoting epigenetic regulator in colon cancer, which may potentially serve as a prognostic marker for tumor metastasis.

Design, synthesis and activity of Mnk1 and Mnk2 selective inhibitors containing thieno[2,3-d]pyrimidine scaffold.

The mitogen-activated protein kinase-interacting kinases 1 and 2 (MNK1 and MNK2) phosphorylate eukaryotic initiation factor 4E (eIF4E) and play important roles in promoting tumorigenesis and metabolic disease. Thus, inhibiting these enzymes might be valuable in the treatment of such conditions. We designed and synthesized a series of 4-((4-fluoro-2-isopropoxyphenyl)amino)-5-methylthieno[2,3-d]pyrimidine derivatives, and evaluated their inhibitory activity against the MNKs. We found 15 compounds that were active as MNK inhibitors and that one in particular, designated MNK-7g, which was potent against MNK1 and substantially more potent against MNK2. The compound MNK-7g did not affect other signaling pathways tested and had no adverse effects on cell viability. As expected from earlier studies, MNK-7g also inhibited cell migration. Therefore, the compound MNK-7g, which forms an ionic bond with Asp226 in MNK2 and possesses a substituted aniline in a thieno[2,3-d] pyrimidine structure, is a promising starting point for the future development of novel drugs for treating or managing cancer and metabolic disease.

mlDEEPre: Multi-Functional Enzyme Function Prediction With Hierarchical Multi-Label Deep Learning.

As a great challenge in bioinformatics, enzyme function prediction is a significant step toward designing novel enzymes and diagnosing enzyme-related diseases. Existing studies mainly focus on the mono-functional enzyme function prediction. However, the number of multi-functional enzymes is growing rapidly, which requires novel computational methods to be developed. In this paper, following our previous work, DEEPre, which uses deep learning to annotate mono-functional enzyme's function, we propose a novel method, mlDEEPre, which is designed specifically for predicting the functionalities of multi-functional enzymes. By adopting a novel loss function, associated with the relationship between different labels, and a self-adapted label assigning threshold, mlDEEPre can accurately and efficiently perform multi-functional enzyme prediction. Extensive experiments also show that mlDEEPre can outperform the other methods in predicting whether an enzyme is a mono-functional or a multi-functional enzyme (mono-functional vs. multi-functional), as well as the main class prediction across different criteria. Furthermore, due to the flexibility of mlDEEPre and DEEPre, mlDEEPre can be incorporated into DEEPre seamlessly, which enables the updated DEEPre to handle both mono-functional and multi-functional predictions without human intervention.

Mycobacterium tuberculosis subverts negative regulatory pathways in human macrophages to drive immunopathology.

Tuberculosis remains a global pandemic and drives lung matrix destruction to transmit. Whilst pathways driving inflammatory responses in macrophages have been relatively well described, negative regulatory pathways are less well defined. We hypothesised that Mycobacterium tuberculosis (Mtb) specifically targets negative regulatory pathways to augment immunopathology. Inhibition of signalling through the PI3K/AKT/mTORC1 pathway increased matrix metalloproteinase-1 (MMP-1) gene expression and secretion, a collagenase central to TB pathogenesis, and multiple pro-inflammatory cytokines. In patients with confirmed pulmonary TB, PI3Kδ expression was absent within granulomas. Furthermore, Mtb infection suppressed PI3Kδ gene expression in macrophages. Interestingly, inhibition of the MNK pathway, downstream of pro-inflammatory p38 and ERK MAPKs, also increased MMP-1 secretion, whilst suppressing secretion of TH1 cytokines. Cross-talk between the PI3K and MNK pathways was demonstrated at the level of eIF4E phosphorylation. Mtb globally suppressed the MMP-inhibitory pathways in macrophages, reducing levels of mRNAs encoding PI3Kδ, mTORC-1 and MNK-1 via upregulation of miRNAs. Therefore, Mtb disrupts negative regulatory pathways at multiple levels in macrophages to drive a tissue-destructive phenotype that facilitates transmission.

Oncogenic MNK signalling regulates the metastasis suppressor NDRG1.

The protein N-myc down-regulated gene 1 (NDRG1) represses tumour metastasis. It is phosphorylated at several sites by serum and glucocorticoid-regulated kinase 1 (SGK1). Here we show that NDRG1 is also regulated by the oncogenic MAP kinase-interacting kinase (MNK) pathway, a target for cancer therapy.Inhibiting MNKs increases the expression of NDRG1 protein and mRNA in breast cancer cells. MNK inhibition also decreases the phosphorylation of NDRG1. Phosphorylation of NDRG1 is reduced in cells lacking MNK1, but not MNK2-knockout cells, indicating that NDRG1 phosphorylation is a specific target for MNK1. However, MNK1 cannot directly phosphorylate NDRG1 in vitro, indicating that additional signalling connections are involved. Taken together, our data indicate that MNK signaling regulates NDRG1 at transcriptional and post-translational levels.We show that SGK1 phosphorylates MNK1 at a conserved site, which represses its activity. NDRG1, SGK1 and the MNKs are implicated in cell migration and metastasis. As expected, knocking-down NDRG1 promoted cell migration. However, whereas MNK inhibition impairs these processes irrespective of NDRG1 levels, SGK inhibition only did so in NDRG1-depleted cells. Thus, MNKs and SGK affect migration/invasion through distinct mechanisms.Our data reveal several novel connections between signalling pathways important for tumour biology.

The MAP kinase-interacting kinases regulate cell migration, vimentin expression and eIF4E/CYFIP1 binding.

The MAP kinase-interacting kinases (Mnk1 and Mnk2) are activated by ERK and are best known for phosphorylating the translation initiation factor eIF4E. Genetic knockout of the Mnks impaired the migration of embryonic fibroblasts both in two-dimensional wound-healing experiments and in three-dimensional migration assays. Furthermore, a novel and selective Mnk inhibitor, Mnk-I1, which potently blocks eIF4E phosphorylation, blocked the migration of fibroblasts and cancer cells, without exerting 'off-target' effects on other signalling pathways such as Erk. Mnk-I1 or genetic knockout of the Mnks decreased the expression of vimentin, a marker of mesenchymal cells, without affecting vimentin mRNA levels. Vimentin protein levels were much lower in Mnk1/2-knockout cells than in controls, although mRNA levels were similar. Our data suggest that the Mnks regulate the translation of the vimentin mRNA and the stability of the vimentin protein. Inhibition or genetic knockout of the Mnks increased the binding of eIF4E to the cytoplasmic FMRP-interacting protein 1 (CYFIP1), which binds the fragile-X mental retardation protein, FMRP, a translational repressor. Since FMRP binds mRNAs for proteins involved in metastasis, the Mnk-dependent release of CYFIP1 from eIF4E is expected to release the repression of translation of FMRP-bound mRNAs, potentially providing a molecular mechanism for the control of cell migration by the Mnks. As Mnk1/2 are not essential for viability, inhibition of the Mnks may be a useful approach to tackling cancer metastasis, a key process contributing to mortality in cancer patients.

pelo is required for high efficiency viral replication.

Viruses hijack host factors for their high speed protein synthesis, but information about these factors is largely unknown. In searching for genes that are involved in viral replication, we carried out a forward genetic screen for Drosophila mutants that are more resistant or sensitive to Drosophila C virus (DCV) infection-caused death, and found a virus-resistant line in which the expression of pelo gene was deficient. Our mechanistic studies excluded the viral resistance of pelo deficient flies resulting from the known Drosophila anti-viral pathways, and revealed that pelo deficiency limits the high level synthesis of the DCV capsid proteins but has no or very little effect on the expression of some other viral proteins, bulk cellular proteins, and transfected exogenous genes. The restriction of replication of other types of viruses in pelo deficient flies was also observed, suggesting pelo is required for high level production of capsids of all kinds of viruses. We show that both pelo deficiency and high level DCV protein synthesis increase aberrant 80S ribosomes, and propose that the preferential requirement of pelo for high level synthesis of viral capsids is at least partly due to the role of pelo in dissociation of stalled 80S ribosomes and clearance of aberrant viral RNA and proteins. Our data demonstrated that pelo is a host factor that is required for high efficiency translation of viral capsids and targeting pelo could be a strategy for general inhibition of viral infection.