Differential regulation of distinct Vps34 complexes by AMPK in nutrient stress and autophagy.

Autophagy is a stress response protecting cells from unfavorable conditions, such as nutrient starvation. The class III phosphatidylinositol-3 kinase, Vps34, forms multiple complexes and regulates both intracellular vesicle trafficking and autophagy induction. Here, we show that AMPK plays a key role in regulating different Vps34 complexes. AMPK inhibits the nonautophagy Vps34 complex by phosphorylating T163/S165 in Vps34 and therefore suppresses overall PI(3)P production and protects cells from starvation. In parallel, AMPK activates the proautophagy Vps34 complex by phosphorylating S91/S94 in Beclin1 to induce autophagy. Atg14L, an autophagy-essential gene present only in the proautophagy Vps34 complex, inhibits Vps34 phosphorylation but increases Beclin1 phosphorylation by AMPK. As such, Atg14L dictates the differential regulation (either inhibition or activation) of different Vps34 complexes in response to glucose starvation. Our study reveals an intricate molecular regulation of Vps34 complexes by AMPK in nutrient stress response and autophagy.

LeafNet: a tool for segmenting and quantifying stomata and pavement cells.

Stomata play important roles in gas and water exchange in leaves. The morphological features of stomata and pavement cells are highly plastic and are regulated during development. However, it is very laborious and time-consuming to collect accurate quantitative data from the leaf surface by manual phenotyping. Here, we introduce LeafNet, a tool that automatically localizes stomata, segments pavement cells (to prepare them for quantification), and reports multiple morphological parameters for a variety of leaf epidermal images, especially bright-field microscopy images. LeafNet employs a hierarchical strategy to identify stomata using a deep convolutional network and then segments pavement cells on stomata-masked images using a region merging method. LeafNet achieved promising performance on test images for quantifying different phenotypes of individual stomata and pavement cells compared with six currently available tools, including StomataCounter, Cellpose, PlantSeg, and PaCeQuant. LeafNet shows great flexibility, and we improved its ability to analyze bright-field images from a broad range of species as well as confocal images using transfer learning. Large-scale images of leaves can be efficiently processed in batch mode and interactively inspected with a graphic user interface or a web server (https://leafnet.whu.edu.cn/). The functionalities of LeafNet could easily be extended and will enhance the efficiency and productivity of leaf phenotyping for many plant biologists.

The polyA tail facilitates splicing of last introns with weak 3' splice sites via PABPN1.

The polyA tail of mRNAs is important for many aspects of RNA metabolism. However, whether and how it regulates pre-mRNA splicing is still unknown. Here, we report that the polyA tail acts as a splicing enhancer for the last intron via the nuclear polyA binding protein PABPN1 in HeLa cells. PABPN1-depletion induces the retention of a group of introns with a weaker 3' splice site, and they show a strong 3'-end bias and mainly locate in nuclear speckles. The polyA tail is essential for PABPN1-enhanced last intron splicing and functions in a length-dependent manner. Tethering PABPN1 to nonpolyadenylated transcripts also promotes splicing, suggesting a direct role for PABPN1 in splicing regulation. Using TurboID-MS, we construct the PABPN1 interactome, including many spliceosomal and RNA-binding proteins. Specifically, PABPN1 can recruit RBM26&27 to promote splicing by interacting with the coiled-coil and RRM domain of RBM27. PABPN1-regulated terminal intron splicing is conserved in mice. Together, our study establishes a novel mode of post-transcriptional splicing regulation via the polyA tail and PABPN1.

MVIP: multi-omics portal of viral infection.

Virus infections are huge threats to living organisms and cause many diseases, such as COVID-19 caused by SARS-CoV-2, which has led to millions of deaths. To develop effective strategies to control viral infection, we need to understand its molecular events in host cells. Virus related functional genomic datasets are growing rapidly, however, an integrative platform for systematically investigating host responses to viruses is missing. Here, we developed a user-friendly multi-omics portal of viral infection named as MVIP (https://mvip.whu.edu.cn/). We manually collected available high-throughput sequencing data under viral infection, and unified their detailed metadata including virus, host species, infection time, assay, and target, etc. We processed multi-layered omics data of more than 4900 viral infected samples from 77 viruses and 33 host species with standard pipelines, including RNA-seq, ChIP-seq, and CLIP-seq, etc. In addition, we integrated these genome-wide signals into customized genome browsers, and developed multiple dynamic charts to exhibit the information, such as time-course dynamic and differential gene expression profiles, alternative splicing changes and enriched GO/KEGG terms. Furthermore, we implemented several tools for efficiently mining the virus-host interactions by virus, host and genes. MVIP would help users to retrieve large-scale functional information and promote the understanding of virus-host interactions.

MaGenDB: a functional genomics hub for Malvaceae plants.

Malvaceae is a family of flowering plants containing many economically important plant species including cotton, cacao and durian. Recently, the genomes of several Malvaceae species have been decoded, and many omics data were generated for individual species. However, no integrative database of multiple species, enabling users to jointly compare and analyse relevant data, is available for Malvaceae. Thus, we developed a user-friendly database named MaGenDB (http://magen.whu.edu.cn) as a functional genomics hub for the plant community. We collected the genomes of 13 Malvaceae species, and comprehensively annotated genes from different perspectives including functional RNA/protein element, gene ontology, KEGG orthology, and gene family. We processed 374 sets of diverse omics data with the ENCODE pipelines and integrated them into a customised genome browser, and designed multiple dynamic charts to present gene/RNA/protein-level knowledge such as dynamic expression profiles and functional elements. We also implemented a smart search system for efficiently mining genes. In addition, we constructed a functional comparison system to help comparative analysis between genes on multiple features in one species or across closely related species. This database and associated tools will allow users to quickly retrieve large-scale functional information for biological discovery.

A mammalian autophagosome maturation mechanism mediated by TECPR1 and the Atg12-Atg5 conjugate.

Autophagy is a major catabolic pathway in eukaryotes associated with a broad spectrum of human diseases. In autophagy, autophagosomes carrying cellular cargoes fuse with lysosomes for degradation. However, the molecular mechanism underlying autophagosome maturation is largely unknown. Here we report that TECPR1 binds to the Atg12-Atg5 conjugate and phosphatidylinositol 3-phosphate (PtdIns[3]P) to promote autophagosome-lysosome fusion. TECPR1 and Atg16 form mutually exclusive complexes with the Atg12-Atg5 conjugate, and TECPR1 binds PtdIns(3)P upon association with the Atg12-Atg5 conjugate. Strikingly, TECPR1 localizes to and recruits Atg5 to autolysosome membrane. Consequently, elimination of TECPR1 leads to accumulation of autophagosomes and blocks autophagic degradation of LC3-II and p62. Finally, autophagosome maturation marked by GFP-mRFP-LC3 is defective in TECPR1-deficient cells. Thus, we propose that the concerted interactions among TECPR1, Atg12-Atg5, and PtdIns(3)P provide the fusion specificity between autophagosomes and lysosomes and that the assembly of this complex initiates the autophagosome maturation process.

Long-term trend in vegetation gross primary production, phenology and their relationships inferred from the FLUXNET data.

Climate-induced changes in plant phenology and physiology plays an important role in control of carbon exchange between terrestrial ecosystems and the atmosphere. Based on dataset during 1997-2014 from 41 flux tower sites (440 site-years) across the northern hemisphere, relationships between long-term trends in start of growing season (SOS), end of growing season (EOS), length of growing season (LOS), maximal gross primary production (GPPmax), and seasonal and annual gross primary production (GPP) were analyzed. Statistical Models of Integrated Phenology and Physiology (SMIPP) were built for predicting the long-term trends in annual GPP. Results showed that SOS advanced and EOS delayed for forest sites, while both SOS and EOS for grassland (GRA) sites delayed. Long-term trends in SOS and EOS of evergreen needle-leaf forests (ENF) sites were greater than those of deciduous broadleaf forests (DBF) sites. Seasonal and annual GPP for forest sites increased, among which long-term trend in annual GPP of ENF sites was the largest. Spring GPP of GRA sites decreased, but annual GPP increased. Strong relationships between long-term trends in phenological and physiological indicators and seasonal GPP were found. Long-term trend in GPPmax had the highest relationship with long-term trend in annual GPP for forest sites, but long-term trend in SOS was the most related to long-term trend in annual GPP for GRA sites. Increases in spring and autumn GPP due to a one-day advance in SOS and delay in EOS for DBF sites were greater than ENF sites. Delay in EOS resulted in more carbon sequestration than advance in SOS for forest sites, while advance in SOS significantly increased spring GPP for GRA sites. The SMIPP model driven by long-term trends in LOS and GPPmax had stronger explanatory power for predicting long-term trend in annual GPP than the SMIPP model driven by long-term trends in SOS, EOS, and GPPmax. Long-term trend in annual GPP was accurately predicted by using the SMIPP model, while long-term trend in annual GPP for GRA sites was more difficult to be captured than the forest sites. Drought and disturbance effects on phenology and physiology were major factors for model uncertainty. This study is helpful to understand changes in phenology and carbon uptake and their differences among different vegetation types and provides a potential way for predicting annual rate of change in carbon uptake through vegetation photosynthesis at a global scale.

Autophagosome targeting and membrane curvature sensing by Barkor/Atg14(L).

The class III phosphatidylinositol 3-kinase (PI3KC3) is crucial for autophagosome biogenesis. It has been long speculated to nucleate the autophagosome membrane, but the biochemical mechanism of such nucleation activity remains unsolved. We recently identified Barkor/Atg14(L) as the targeting factor for PI3KC3 to autophagosome membrane. Here, we show that we have characterized the region of Barkor/Atg14(L) required for autophagosome targeting and identified the BATS [Barkor/Atg14(L) autophagosome targeting sequence] domain at the carboxyl terminus of Barkor. Bioinformatics and mutagenesis analyses revealed that the BATS domain binds to autophagosome membrane via the hydrophobic surface of an intrinsic amphipathic alpha helix. BATS puncta overlap with Atg16 and LC3, and partially with DFCP1, in a stress-inducible manner. Ectopically expressed BATS accumulates on highly curved tubules that likely represent intermediate autophagic structures. PI3KC3 recruitment and autophagy stimulation by Barkor/Atg14(L) require the BATS domain. Furthermore, our biochemical analyses indicate that the BATS domain directly binds to the membrane, and it favors membrane composed of phosphatidylinositol 3-phosphate [PtdIns(3)P] and phosphatidylinositol 4,5-biphosphate [PtdIns(4,5)P2]. By binding preferentially to curved membranes incorporated with PtdIns(3)P but not PtdIns(4,5)P2, the BATS domain is capable of sensing membrane curvature. Thus, we propose a novel model of PI3KC3 autophagosome membrane nucleation in which its autophagosome-specific adaptor, Barkor, accumulates on highly curved PtdIns(3)P enriched autophagic membrane via its BATS domain to sense and maintain membrane curvature.

Atlas of interactions between SARS-CoV-2 macromolecules and host proteins.

The proteins and RNAs of viruses extensively interact with host proteins after infection. We collected and reanalyzed all available datasets of protein-protein and RNA-protein interactions related to SARS-CoV-2. We investigated the reproducibility of those interactions and made strict filters to identify highly confident interactions. We systematically analyzed the interaction network and identified preferred subcellular localizations of viral proteins, some of which such as ORF8 in ER and ORF7A/B in ER membrane were validated using dual fluorescence imaging. Moreover, we showed that viral proteins frequently interact with host machinery related to protein processing in ER and vesicle-associated processes. Integrating the protein- and RNA-interactomes, we found that SARS-CoV-2 RNA and its N protein closely interacted with stress granules including 40 core factors, of which we specifically validated G3BP1, IGF2BP1, and MOV10 using RIP and Co-IP assays. Combining CRISPR screening results, we further identified 86 antiviral and 62 proviral factors and associated drugs. Using network diffusion, we found additional 44 interacting proteins including two proviral factors previously validated. Furthermore, we showed that this atlas could be applied to identify the complications associated with COVID-19. All data are available in the AIMaP database (https://mvip.whu.edu.cn/aimap/) for users to easily explore the interaction map.

Ligand-induced activation and G protein coupling of prostaglandin F2α receptor.

Prostaglandin F2α (PGF2α), an endogenous arachidonic acid metabolite, regulates diverse physiological functions in many tissues and cell types through binding and activation of a G-protein-coupled receptor (GPCR), the PGF2α receptor (FP), which also is the primary therapeutic target for glaucoma and several other diseases. Here, we report cryo-electron microscopy (cryo-EM) structures of the human FP bound to endogenous ligand PGF2α and anti-glaucoma drugs LTPA and TFPA at global resolutions of 2.67 Å, 2.78 Å, and 3.14 Å. These structures reveal distinct features of FP within the lipid receptor family in terms of ligand binding selectivity, its receptor activation, and G protein coupling mechanisms, including activation in the absence of canonical PIF and ERY motifs and Gq coupling through direct interactions with receptor transmembrane helix 1 and intracellular loop 1. Together with mutagenesis and functional studies, our structures reveal mechanisms of ligand recognition, receptor activation, and G protein coupling by FP, which could facilitate rational design of FP-targeting drugs.

The RUN domain of rubicon is important for hVps34 binding, lipid kinase inhibition, and autophagy suppression.

The class III phosphatidylinositol 3-kinase (PI3KC3) plays a central role in autophagy. Rubicon, a RUN domain-containing protein, is newly identified as a PI3KC3 subunit through its association with Beclin 1. Rubicon serves as a negative regulator of PI3KC3 and autophagosome maturation. The molecular mechanism underlying the PI3KC3 and autophagy inhibition by Rubicon is largely unknown. Here, we demonstrate that Rubicon interacts with the PI3KC3 catalytic subunit hVps34 via its RUN domain. The RUN domain contributes to the efficient inhibition of PI3KC3 lipid kinase activity by Rubicon. Furthermore, a Rubicon RUN domain deletion mutant fails to complement the autophagy deficiency in Rubicon-depleted cells. Hence, these results reveal a critical role of the Rubicon RUN domain in PI3KC3 and autophagy regulation.

Osmotic stress activates RIPK3/MLKL-mediated necroptosis by increasing cytosolic pH through a plasma membrane Na+/H+ exchanger.

Necroptosis is a form of cell death triggered by stimuli such as the tumor necrosis factor family of cytokines, which induce necrotic cell death through the RIPK1-RIPK3-MLKL pathway. We report here that necroptosis is also activated by extracellular osmotic stresses. Unlike the previously identified inducers of necroptosis, osmotic stress stimulated necroptosis through the direct activation of the kinase activity of RIPK3 by an increase in cytosolic pH mediated by the Na+/H+ exchanger SLC9A1. Knockout, knockdown, or chemical inhibition of SLC9A1 blocked necroptosis induced by osmotic stresses. Moreover, setting intracellular pH at above-physiological values directly activated RIPK3 and necroptosis. The activation of RIPK3 by osmotic stresses did not require its RHIM domain, the protein-interacting domain required for the activation of RIPK3 when cells respond to other previously identified necroptotic stimuli. These results thus delineate a pathway that activates necroptosis in response to osmotic stresses.

Identification of Barkor as a mammalian autophagy-specific factor for Beclin 1 and class III phosphatidylinositol 3-kinase.

Autophagy mediates the cellular response to nutrient deprivation, protein aggregation, and pathogen invasion in human. Dysfunction of autophagy has been implicated in multiple human diseases including cancer. The identification of novel autophagy factors in mammalian cells will provide critical mechanistic insights into how this complicated cellular pathway responds to a broad range of challenges. Here, we report the cloning of an autophagy-specific protein that we called Barkor (Beclin 1-associated autophagy-related key regulator) through direct interaction with Beclin 1 in the human phosphatidylinositol 3-kinase class III complex. Barkor shares 18% sequence identity and 32% sequence similarity with yeast Atg14. Elimination of Barkor expression by RNA interference compromises starvation- and rapamycin-induced LC3 lipidation and autophagosome formation. Overexpression of Barkor leads to autophagy activation and increased number and enlarged volume of autophagosomes. Tellingly, Barkor is also required for suppression of the autophagy-mediated intracellular survival of Salmonella typhimurium in mammalian cells. Mechanistically, Barkor competes with UV radiation resistance associated gene product (UVRAG) for interaction with Beclin 1, and the complex formation of Barkor and Beclin1 is required for their localizations to autophagosomes. Therefore, we define a regulatory signaling pathway mediated by Barkor that positively controls autophagy through Beclin 1 and represents a potential target for drug development in the treatment of human diseases implicated in autophagic dysfunction.

NEPdb: A Database of T-Cell Experimentally-Validated Neoantigens and Pan-Cancer Predicted Neoepitopes for Cancer Immunotherapy.

T-cell recognition of somatic mutation-derived cancer neoepitopes can lead to tumor regression. Due to the difficulty to identify effective neoepitopes, constructing a database for sharing experimentally validated cancer neoantigens will be beneficial to precise cancer immunotherapy. Meanwhile, the routine neoepitope prediction in silico is important but laborious for clinical use. Here we present NEPdb, a database that contains more than 17,000 validated human immunogenic neoantigens and ineffective neoepitopes within human leukocyte antigens (HLAs) via curating published literature with our semi-automatic pipeline. Furthermore, NEPdb also provides pan-cancer level predicted HLA-I neoepitopes derived from 16,745 shared cancer somatic mutations, using state-of-the-art predictors. With a well-designed search engine and visualization modes, this database would enhance the efficiency of neoantigen-based cancer studies and treatments. NEPdb is freely available at http://nep.whu.edu.cn/.

ASER: Animal Sex Reversal Database.

Sex reversal, representing extraordinary sexual plasticity during the life cycle, not only triggers reproduction in animals but also affects reproductive and endocrine system-related diseases and cancers in humans. Sex reversal has been broadly reported in animals; however, an integrated resource hub of sex reversal information is still lacking. Here, we constructed a comprehensive database named ASER (Animal Sex Reversal) by integrating sex reversal-related data of 18 species from teleostei to mammalia. We systematically collected 40,018 published papers and mined the sex reversal-associated genes (SRGs), including their regulatory networks, from 1611 core papers. We annotated homologous genes and computed conservation scores for whole genomes across the 18 species. Furthermore, we collected available RNA-seq datasets and investigated the expression dynamics of SRGs during sex reversal or sex determination processes. In addition, we manually annotated 550 in situ hybridization (ISH), fluorescence in situ hybridization (FISH), and immunohistochemistry (IHC) images of SRGs from the literature and described their spatial expression in the gonads. Collectively, ASER provides a unique and integrated resource for researchers to query and reuse organized data to explore the mechanisms and applications of SRGs in animal breeding and human health. The ASER database is publicly available at http://aser.ihb.ac.cn/.

Flotillin-mediated endocytosis and ALIX-syntenin-1-mediated exocytosis protect the cell membrane from damage caused by necroptosis.

Necroptosis is a form of regulated necrosis that is implicated in various human diseases including Alzheimer's disease. Necroptosis requires the translocation of the pseudokinase MLKL from the cytosol to the plasma membrane after its phosphorylation by the kinase RIPK3. Using protein cross-linking followed by affinity purification, we detected the lipid raft-associated proteins flotillin-1 and flotillin-2 and the ESCRT-associated proteins ALIX and syntenin-1 in membrane-localized MLKL immunoprecipitates. Phosphorylated MLKL was removed from membranes through either flotillin-mediated endocytosis followed by lysosomal degradation or ALIX-syntenin-1-mediated exocytosis. Thus, cells undergoing necroptosis need to overcome these independent suppressive mechanisms before plasma membrane disruption can occur.

MTORC1-mediated NRBF2 phosphorylation functions as a switch for the class III PtdIns3K and autophagy.

NRBF2/Atg38 has been identified as the fifth subunit of the macroautophagic/autophagic class III phosphatidylinositol 3-kinase (PtdIns3K) complex, along with ATG14/Barkor, BECN1/Vps30, PIK3R4/p150/Vps15 and PIK3C3/Vps34. However, its functional mechanism and regulation are not fully understood. Here, we report that NRBF2 is a fine tuning regulator of PtdIns3K controlled by phosphorylation. Human NRBF2 is phosphorylated by MTORC1 at S113 and S120. Upon nutrient starvation or MTORC1 inhibition, NRBF2 phosphorylation is diminished. Phosphorylated NRBF2 preferentially interacts with PIK3C3/PIK3R4. Suppression of NRBF2 phosphorylation by MTORC1 inhibition alters its binding preference from PIK3C3/PIK3R4 to ATG14/BECN1, leading to increased autophagic PtdIns3K complex assembly, as well as enhancement of ULK1 protein complex association. Consequently, NRBF2 in its unphosphorylated form promotes PtdIns3K lipid kinase activity and autophagy flux, whereas its phosphorylated form blocks them. This study reveals NRBF2 as a critical molecular switch of PtdIns3K and autophagy activation, and its on/off state is precisely controlled by MTORC1 through phosphorylation.

Membrane curvature response in autophagy.

Membrane trafficking is key for signal transduction, cargo transportation, and in the case of autophagy, delivering cytoplasmic substrates to the lysosome for degradation. Autophagy requires the formation of a unique double membrane vesicle, the autophagosome. However, the mechanism by which the autophagosome forms is unknown. Our recent study focused on the role of Barkor/Atg14(L) in targeting the autophagy-specific class III phosphatidylinositol-3-kinase (PtdIns3KC3) complex to the early autophagosome has implicated this complex in autophagosome formation. This study found that the BATS domain of Barkor targets the PtdIns3KC3 complex to early autophagic structures and senses highly curved membranes enriched in phosphatidylinositol-3-phosphate (PtdIns(3)P). Consequently, this study uncovered an exciting new role for the PtdIns3KC3 complex as a potential inducer of autophagosome formation.

Keap1 facilitates p62-mediated ubiquitin aggregate clearance via autophagy.

The accumulation of ubiquitin-positive protein aggregates has been implicated in the pathogenesis of neurodegenerative diseases, heart disease and diabetes. Emerging evidence indicates that the autophagy lysosomal pathway plays a critical role in the clearance of ubiquitin aggregates, a process that is mediated by the ubiquitin binding protein p62. In addition to binding ubiquitin, p62 also interacts with LC3 and transports ubiquitin conjugates to autophagosomes for degradation. The exact regulatory mechanism of this process is still largely unknown. Here we report the identification of Keap1 as a binding partner for p62 and LC3. Keap1 inhibits Nrf2 by sequestering it in the cytosol and preventing its translocation to the nucleus and activation of genes involved in the oxidative stress response. In this study, we found that Keap1 interacts with p62 and LC3 in a stress-inducible manner, and that Keap1 colocalizes with LC3 and p62 in puromycin-induced ubiquitin aggregates. Moreover, p62 serves as a bridge between Keap1 and ubiquitin aggregates and autophagosomes. Finally, genetic ablation of Keap1 leads to the accumulation of ubiquitin aggregates, increased cytotoxicity of misfolded protein aggregates, and defective activation of autophagy. Therefore, this study assigns a novel positive role of Keap1 in upregulating p62-mediated autophagic clearance of ubiquitin aggregates.