The Drosophila histone methyl-transferase SET1 coordinates multiple signaling pathways in regulating male germline stem cell maintenance and differentiation.

Many tissue-specific adult stem cell lineages maintain a balance between proliferation and differentiation. Here, we study how the H3K4me3 methyltransferase, Set1, regulates early-stage male germ cells in Drosophila. Early-stage germline-specific knockdown of set1 results in temporally progressed defects, arising as germ cell loss and developing into overpopulated early-stage germ cells. These germline defects also impact the niche architecture and cyst stem cell lineage non-cell-autonomously. Additionally, wild-type Set1, but not the catalytically inactive Set1, rescues the set1 knockdown phenotypes, highlighting the functional importance of the methyl-transferase activity of Set1. Further, RNA-seq experiments reveal key signaling pathway components, such as the JAK-STAT pathway stat92E and the BMP pathway mad genes that are upregulated upon set1 knockdown. Genetic interaction assays support the functional relationships between set1 and JAK-STAT or BMP pathways, as both stat92E and mad mutations suppress the set1 knockdown phenotypes. These findings enhance our understanding of the balance between proliferation and differentiation in an adult stem cell lineage. The germ cell loss followed by over-proliferation phenotype when inhibiting a histone methyl-transferase also raise concerns about using their inhibitors in cancer therapy.

Large-scale CRISPR screen reveals context-specific genetic regulation of retinal ganglion cell regeneration.

Many genes are known to regulate retinal regeneration following widespread tissue damage. Conversely, genes controlling regeneration following limited cell loss, per degenerative diseases, are undefined. As stem/progenitor cell responses scale to injury levels, understanding how the extent and specificity of cell loss impact regenerative processes is important. Here, transgenic zebrafish enabling selective retinal ganglion cell (RGC) ablation were used to identify genes that regulate RGC regeneration. A single cell multiomics-informed screen of 101 genes identified seven knockouts that inhibited and eleven that promoted RGC regeneration. Surprisingly, 35 of 36 genes known/implicated as being required for regeneration following widespread retinal damage were not required for RGC regeneration, and seven even enhanced regeneration kinetics, including proneural factors neurog1, olig2, and ascl1a. Mechanistic analyses revealed ascl1a disruption increased the propensity of progenitor cells to produce RGCs; i.e., increased "fate bias". These data demonstrate plasticity in how Müller glia can convert to a stem-like state and context-specificity in how genes function during regeneration. Increased understanding of how the regeneration of disease-relevant cell types is specifically controlled will support the development of disease-tailored regenerative therapeutics.

New pathways to neurogenesis: Insights from injury-induced retinal regeneration.

The vertebrate retina is a tractable system for studying control of cell neurogenesis and cell fate specification. During embryonic development, retinal neurogenesis is under strict temporal regulation, with cell types generated in fixed but overlapping temporal intervals. The temporal sequence and relative numbers of retinal cell types generated during development are robust and show minimal experience-dependent variation. In many cold-blooded vertebrates, acute retinal injury induces a different form of neurogenesis, where Müller glia reprogram into retinal progenitor-like cells that selectively regenerate retinal neurons lost to injury. The extent to which the molecular mechanisms controlling developmental and injury-induced neurogenesis resemble one another has long been unclear. However, a recent study in zebrafish has shed new light on this question, using single-cell multiomic analysis to show that selective loss of different retinal cell types induces the formation of fate-restricted Müller glia-derived progenitors that differ both from one another and from progenitors in developing retina. Here, we discuss the broader implications of these findings, and their possible therapeutic relevance.

Transcriptomic comparison of two selective retinal cell ablation paradigms in zebrafish reveals shared and cell-specific regenerative responses.

Retinal Müller glia (MG) can act as stem-like cells to generate new neurons in both zebrafish and mice. In zebrafish, retinal regeneration is innate and robust, resulting in the replacement of lost neurons and restoration of visual function. In mice, exogenous stimulation of MG is required to reveal a dormant and, to date, limited regenerative capacity. Zebrafish studies have been key in revealing factors that promote regenerative responses in the mammalian eye. Increased understanding of how the regenerative potential of MG is regulated in zebrafish may therefore aid efforts to promote retinal repair therapeutically. Developmental signaling pathways are known to coordinate regeneration following widespread retinal cell loss. In contrast, less is known about how regeneration is regulated in the context of retinal degenerative disease, i.e., following the loss of specific retinal cell types. To address this knowledge gap, we compared transcriptomic responses underlying regeneration following targeted loss of rod photoreceptors or bipolar cells. In total, 2,531 differentially expressed genes (DEGs) were identified, with the majority being paradigm specific, including during early MG activation phases, suggesting the nature of the injury/cell loss informs the regenerative process from initiation onward. For example, early modulation of Notch signaling was implicated in the rod but not bipolar cell ablation paradigm and components of JAK/STAT signaling were implicated in both paradigms. To examine candidate gene roles in rod cell regeneration, including several immune-related factors, CRISPR/Cas9 was used to create G0 mutant larvae (i.e., "crispants"). Rod cell regeneration was inhibited in stat3 crispants, while mutating stat5a/b, c7b and txn accelerated rod regeneration kinetics. These data support emerging evidence that discrete responses follow from selective retinal cell loss and that the immune system plays a key role in regulating "fate-biased" regenerative processes.

Terahertz s-SNOM Imaging of a Single Cell with Nanoscale Resolution.

Terahertz scattering scanning near-field optical microscopy is a robust spectral detection technique with a nanoscale resolution. However, there are still major challenges in investigating the heterogeneity of cell membrane components in individual cells. Here, we present a novel and comprehensive analytical approach for detecting and investigating heterogeneity in cell membrane components at the single-cell level. In comparison to the resolution of the topographical atomic force microscopy image, the spatial resolution of the terahertz near-field amplitude image is 3 times that of the former. This ultrafine resolution enables the compositional distribution in the cell membrane, such as the distribution of extracellular vesicles, to be finely characterized. Furthermore, via extraction of the near-field absorption images at specific frequencies, the visualization and compositional difference analysis of cell membrane components can be presented in detail. These findings have significant implications for the intuitive and visual analysis of cell development and disease evolutionary pathways.

CellAnn: a comprehensive, super-fast, and user-friendly single-cell annotation web server.

MOTIVATION: Single-cell sequencing technology has become a routine in studying many biological problems. A core step of analyzing single-cell data is the assignment of cell clusters to specific cell types. Reference-based methods are proposed for predicting cell types for single-cell clusters. However, the scalability and lack of preprocessed reference datasets prevent them from being practical and easy to use. RESULTS: Here, we introduce a reference-based cell annotation web server, CellAnn, which is super-fast and easy to use. CellAnn contains a comprehensive reference database with 204 human and 191 mouse single-cell datasets. These reference datasets cover 32 organs. Furthermore, we developed a cluster-to-cluster alignment method to transfer cell labels from the reference to the query datasets, which is superior to the existing methods with higher accuracy and higher scalability. Finally, CellAnn is an online tool that integrates all the procedures in cell annotation, including reference searching, transferring cell labels, visualizing results, and harmonizing cell annotation labels. Through the user-friendly interface, users can identify the best annotation by cross-validating with multiple reference datasets. We believe that CellAnn can greatly facilitate single-cell sequencing data analysis. AVAILABILITY AND IMPLEMENTATION: The web server is available at www.cellann.io, and the source code is available at https://github.com/Pinlyu3/CellAnn_shinyapp.

Protein acetylation microarray reveals that NuA4 controls key metabolic target regulating gluconeogenesis.

Histone acetyltransferases (HATs) and histone deacetylases (HDACs) conduct many critical functions through nonhistone substrates in metazoans, but only chromatin-associated nonhistone substrates are known in Saccharomyces cerevisiae. Using yeast proteome microarrays, we identified and validated many nonchromatin substrates of the essential nucleosome acetyltransferase of H4 (NuA4) complex. Among these, acetylation sites (Lys19 and 514) of phosphoenolpyruvate carboxykinase (Pck1p) were determined by tandem mass spectrometry. Acetylation at Lys514 was crucial for enzymatic activity and the ability of yeast cells to grow on nonfermentable carbon sources. Furthermore, Sir2p deacetylated Pck1p both in vitro and in vivo. Loss of Pck1p activity blocked the extension of yeast chronological life span caused by water starvation. In human hepatocellular carcinoma (HepG2) cells, human Pck1 acetylation and glucose production were dependent on TIP60, the human homolog of ESA1. Our findings demonstrate a regulatory function for the NuA4 complex in glucose metabolism and life span by acetylating a critical metabolic enzyme.

Profiling the human protein-DNA interactome reveals ERK2 as a transcriptional repressor of interferon signaling.

Protein-DNA interactions (PDIs) mediate a broad range of functions essential for cellular differentiation, function, and survival. However, it is still a daunting task to comprehensively identify and profile sequence-specific PDIs in complex genomes. Here, we have used a combined bioinformatics and protein microarray-based strategy to systematically characterize the human protein-DNA interactome. We identified 17,718 PDIs between 460 DNA motifs predicted to regulate transcription and 4,191 human proteins of various functional classes. Among them, we recovered many known PDIs for transcription factors (TFs). We identified a large number of unanticipated PDIs for known TFs, as well as for previously uncharacterized TFs. We also found that over three hundred unconventional DNA-binding proteins (uDBPs)--which include RNA-binding proteins, mitochondrial proteins, and protein kinases--showed sequence-specific PDIs. One such uDBP, ERK2, acts as a transcriptional repressor for interferon gamma-induced genes, suggesting important biological roles for such proteins.

Lin28A Binds Active Promoters and Recruits Tet1 to Regulate Gene Expression.

Lin28, a well-known RNA-binding protein, regulates diverse cellular properties. All physiological functions of Lin28A characterized so far have been attributed to its repression of let-7 miRNA biogenesis or modulation of mRNA translational efficiency. Here we show that Lin28A directly binds to a consensus DNA sequence in vitro and in mouse embryonic stem cells in vivo. ChIP-seq and RNA-seq reveal enrichment of Lin28A binding around transcription start sites and a positive correlation between its genomic occupancy and expression of many associated genes. Mechanistically, Lin28A recruits 5-methylcytosine-dioxygenase Tet1 to genomic binding sites to orchestrate 5-methylcytosine and 5-hydroxymethylcytosine dynamics. Either Lin28A or Tet1 knockdown leads to dysregulated DNA methylation and expression of common target genes. These results reveal a surprising role for Lin28A in transcriptional regulation via epigenetic DNA modifications and have implications for understanding mechanisms underlying versatile functions of Lin28A in mammalian systems.

scEnhancer: a single-cell enhancer resource with annotation across hundreds of tissue/cell types in three species.

Previous studies on enhancers and their target genes were largely based on bulk samples that represent 'average' regulatory activities from a large population of millions of cells, masking the heterogeneity and important effects from the sub-populations. In recent years, single-cell sequencing technology has enabled the profiling of open chromatin accessibility at the single-cell level (scATAC-seq), which can be used to annotate the enhancers and promoters in specific cell types. A comprehensive resource is highly desirable for exploring how the enhancers regulate the target genes at the single-cell level. Hence, we designed a single-cell database scEnhancer (http://enhanceratlas.net/scenhancer/), covering 14 527 776 enhancers and 63 658 600 enhancer-gene interactions from 1 196 906 single cells across 775 tissue/cell types in three species. An unsupervised learning method was employed to sort and combine tens or hundreds of single cells in each tissue/cell type to obtain the consensus enhancers. In addition, we utilized a cis-regulatory network algorithm to identify the enhancer-gene connections. Finally, we provided a user-friendly platform with seven useful modules to search, visualize, and browse the enhancers/genes. This database will facilitate the research community towards a functional analysis of enhancers at the single-cell level.

The CombE-IDMS Assay as an Alternate Potency Method for Adjuvanted Quadrivalent Influenza Vaccines.

The potency of all currently licensed inactivated influenza viral vaccines is assayed by the single radial immunodiffusion (SRID) method. SRID relies upon antisera and reference antigen reagents which are produced, standardized, and distributed in the mass quantities needed for vaccine manufacturers only after a significant amount of time has elapsed from the seasonal strain recommendations issued by the WHO; this time delay is exacerbated under conditions of an emerging pandemic. Previously, the limited trypsin digestion isotope dilution mass spectrometry (LTD-IDMS) method, which does not require antisera or reference antigens, demonstrated comparable quantitation of immunologically active hemagglutinin, the primary viral antigen, to SRID in stressed vaccine materials. Here, we demonstrate a streamlined improvement to the LTD-IDMS method by eliminating the need for its precipitation and washing steps, saving time and labor in the sample preparation process while paving the way for plate-based high-throughput analysis. This is accomplished using dissimilar proteases in the pretreatment (a combination of chymotrypsin and elastase) and analytical (trypsin) digestion steps so that any pretreatment digests will not cause interference while monitoring analytical tryptic digests by IDMS. The combination of enzymes (CombE)-IDMS method is tested alongside LTD-IDMS and SRID for the first time on MF59 adjuvanted seasonal cell-based quadrivalent influenza vaccines (aQIVc) under stressed conditions of heating, oxidation, lowered and elevated pH, and freeze-thaw. Overall, a correlation in the degradation trend is observed between CombE-IDMS and SRID in the four strains of the quadrivalent formulation, highlighting the method's stability indicating capability as a rapid alternate potency assay in a highly complex formulation of aQIVc.

LRLoop: a method to predict feedback loops in cell-cell communication.

MOTIVATION: Intercellular communication (i.e. cell-cell communication) plays an essential role in multicellular organisms coordinating various biological processes. Previous studies discovered that feedback loops between two cell types are a widespread and vital signaling motif regulating development, regeneration and cancer progression. While many computational methods have been developed to predict cell-cell communication based on gene expression datasets, these methods often predict one-directional ligand-receptor interactions from sender to receiver cells and are not suitable to identify feedback loops. RESULTS: Here, we describe ligand-receptor loop (LRLoop), a new method for analyzing cell-cell communication based on bi-directional ligand-receptor interactions, where two pairs of ligand-receptor interactions are identified that are responsive to each other and thereby form a closed feedback loop. We first assessed LRLoop using bulk datasets and found our method significantly reduces the false positive rate seen with existing methods. Furthermore, we developed a new strategy to assess the performance of these methods in single-cell datasets. We used the between-tissue interactions as an indicator of potential false-positive prediction and found that LRLoop produced a lower fraction of between-tissue interactions than traditional methods. Finally, we applied LRLoop to the single-cell datasets obtained from retinal development. We discovered many new bi-directional ligand-receptor interactions among individual cell types that potentially control proliferation, neurogenesis and/or cell fate specification. AVAILABILITY AND IMPLEMENTATION: An R package is available at https://github.com/Pinlyu3/LRLoop. The source code can be found at figshare (https://doi.org/10.6084/m9.figshare.20126138.v1). The datasets can be found at figshare (https://doi.org/10.6084/m9.figshare.20126021.v1). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hydrogel microenvironment contributes to chemical-induced differentiation of mesenchymal stem cells: single-cell infrared microspectroscopy characterization.

Stem cell microenvironment plays vital roles in directing cell proliferation and differentiation. Due to the tiny biochemical changes in the early stage of stem cell development, technical challenges to characterize the potential effects of environmental signals remain. In this work, we have introduced synchrotron radiation-based Fourier transform infrared microspectroscopy to evaluate the synergistic effects of physical and chemical factors on stem cell differentiation at the single-cell level. By using principal component analysis and cell-cell Euclidean distance calculation, the phenotypic heterogeneity changes during stem cell osteogenesis induced by lithium chloride or Wnt5a protein loaded in the polyvinyl alcohol (PVA) hydrogel were characterized in detail. The results demonstrated that PVA hydrogel could lead to the distinct effects between low-concentration lithium and wnt5a on human mesenchymal stem cells, suggesting a vital role of niche signals in Wnt pathway. These findings highlight the importance of microenvironment to the chemical-induced effects on stem cell differentiation and also provide a label-free, noninvasive method to sensitively identify the niche function in stem cell biology.

Using 24-h intraocular pressure-related patterns to identify open-angle glaucoma in thyroid eye disease.

PURPOSE: Our study aims to develop a diagnostic model using 24-h intraocular pressure (IOP) patterns to differentiate between open-angle glaucoma (OAG) and dysthyroid optic neuropathy (DON) in thyroid eye disease (TED) patients with glaucoma-like symptoms. METHODS: TED patients with elevated IOP, abnormal optic disc, and/or visual fields were prospectively recruited. The subjects whose symptoms were relieved by DON first-line treatments were divided into the DON group, and the subjects with previous diagnosis of OAG before TED onset were divided into the OAG group. The 24-h IOP was monitored by Tono-Pen in a sitting position during awake time and in a supine position during sleep time. All subjects were divided into a training set and a testing set. The diagnostic models were generated from training set by using either IOP curve-derived parameters or principal component (PC) factors. The discrimination ability was tested in training set based on area under curve (AUC), and the calibration ability was verified in testing set by Hosmer-Lemeshow goodness-of-fit. The sensitivity and specificity were calculated by two-by-two table with the cutoff value determined by Youden's index. RESULTS: Thirty-two cases were recruited in each group. The 24-h IOP curves revealed a nocturnal pattern in both groups, with the acrophase moving slightly forward in the DON group (21:00 pm-24:00 pm) compared to the OAG group (22:00 pm-3:00 am). Several IOP curve-derived parameters differed between the two groups, with larger amplitude during sleep time (P < 0.000) and longer duration of IOP ≥ 21 mmHg at awake time (P = 0.004) in the DON group than the OAG group. However, the diagnostic model generated from IOP parameters showed poor reliability (P = 0.001) in calibration test and was rejected. The other model built on PC factors achieved good performance of discrimination (AUC = 0.943) and calibration (P = 0.139) with a sensitivity of 87.50% and a specificity of 95.83% at cutoff value of 0.538 to identify OAG cases. CONCLUSION: The diagnostic model facilitates discrimination between OAG and DON in TED patients based on 24-h IOP-related patterns. TRIAL REGISTRATION: This work was registered on Chinese Clinical Trial Registry (ChiCTR1900025394).

Multipurpose ultrasonographic characteristics of primary uveal MALT lymphoma.

PURPOSE: To investigate the ultrasonographic features in patients with primary uveal mucosa-associated lymphoid tissue (MALT) lymphoma. METHODS: Medical records of 12 patients (13 eyes) diagnosed with primary uveal MALT lymphoma between September 2014 and September 2021 were retrospectively reviewed. Ultrasonography, B-scan ultrasonography, color Doppler flow imaging, and ultrasound biomicroscopy findings were retrieved from the medical records. RESULTS: Mean age of the included patients was 59.4 ± 8.6 years. Typical ultrasonographic features of the choroidal infiltrates were flat, diffuse, and thickened, with low and homogenous internal reflectivity and with rich arterial blood flow from posterior ciliary arterioles. The mean thickness of the choroidal infiltrates was 1.34 ± 0.68 mm (n = 13). Most of the affected eyes had posterior episcleral extensions, with a mean thickness of 1.66 ± 1.21 mm (n = 12). Typical crescent-like posterior episcleral extensions were detected in nine eyes (69.2%). In six eyes, the blood flow from the choroidal infiltrates communicated with the episcleral extensions. In the ciliary body, the mean thickness of the infiltrates was 1.08 ± 0.43 mm (n = 9), and seven eyes (77.8%) had 360° ring-like infiltrations. The initial best-corrected visual acuity (BCVA) was significantly correlated with the final BCVA after treatment (p < 0.001). CONCLUSION: Multipurpose ultrasonographic imaging revealed the unique characteristics of the primary uveal MALT lymphoma and is helpful in the diagnosis of this rare disease.

Identification of Novel Serological Autoantibodies in Takayasu Arteritis Patients Using HuProt Arrays.

To identify novel autoantibodies of Takayasu arteritis (TAK) using HuProt array-based approach, a two-phase approach was adopted. In Phase I, serum samples collected from 40 TAK patients, 15 autoimmune disease patients, and 20 healthy subjects were screened to identify TAK-specific autoantibodies using human protein (HuProt) arrays. In phase II, the identified candidate autoantibodies were validated with TAK-focused arrays using an additional cohort comprised of 109 TAK patients, 110 autoimmune disease patients, and 96 healthy subjects. Subsequently, the TAK-specific autoantibodies validated in phase II were further confirmed using western blot analysis. We identified and validated eight autoantibodies as potential TAK-specific diagnostic biomarkers, including anti-SPATA7, -QDPR, -SLC25A2, -PRH2, -DIXDC1, -IL17RB, -ZFAND4, and -NOLC1 antibodies, with AUC of 0.803, 0.801, 0.780, 0.696, 0.695, 0.678, 0.635, and 0.613, respectively. SPATA7 could distinguish TAK from healthy and disease controls with 73.4% sensitivity at 85.4% specificity, while QDPR showed 71.6% sensitivity at 86.4% specificity. SLC25A22 showed the highest sensitivity of 80.7%, but at lower specificity of 67.0%. In addition, PRH2, IL17RB, and NOLC1 showed good specificities of 88.3%, 85.9%, and 86.9%, respectively, but at lower sensitivities (<50%). Finally, DIXDC1 and ZFAND4 showed moderate performance as compared with the other autoantibodies. Using a decision tree model, we could reach a specificity of 94.2% with AUC of 0.843, a significantly improved performance as compared with that by each individual biomarker. The performances of three autoantibodies, namely anti-SPATA7, -QDPR, and -PRH2, were successfully confirmed with western blot analysis. Using this two-phase strategy, we identified and validated eight novel autoantibodies as TAK-specific biomarker candidates, three of which could be readily adopted in a clinical setting.

Proteome Landscape of Epithelial-to-Mesenchymal Transition (EMT) of Retinal Pigment Epithelium Shares Commonalities With Malignancy-Associated EMT.

Stress and injury to the retinal pigment epithelium (RPE) often lead to dedifferentiation and epithelial-to-mesenchymal transition (EMT). These processes have been implicated in several retinal diseases, including proliferative vitreoretinopathy, diabetic retinopathy, and age-related macular degeneration. Despite the importance of RPE-EMT and the large body of data characterizing malignancy-related EMT, comprehensive proteomic studies to define the protein changes and pathways underlying RPE-EMT have not been reported. This study sought to investigate the temporal protein expression changes that occur in a human-induced pluripotent stem cell-based RPE-EMT model. We utilized multiplexed isobaric tandem mass tag labeling followed by high-resolution tandem MS for precise and in-depth quantification of the RPE-EMT proteome. We have identified and quantified 7937 protein groups in our tandem mass tag-based MS analysis. We observed a total of 532 proteins that are differentially regulated during RPE-EMT. Furthermore, we integrated our proteomic data with prior transcriptomic (RNA-Seq) data to provide additional insights into RPE-EMT mechanisms. To validate these results, we have performed a label-free single-shot data-independent acquisition MS study. Our integrated analysis indicates both the commonality and uniqueness of RPE-EMT compared with malignancy-associated EMT. Our comparative analysis also revealed that multiple age-related macular degeneration-associated risk factors are differentially regulated during RPE-EMT. Together, our integrated dataset provides a comprehensive RPE-EMT atlas and resource for understanding the molecular signaling events and associated biological pathways that underlie RPE-EMT onset. This resource has already facilitated the identification of chemical modulators that could inhibit RPE-EMT, and it will hopefully aid in ongoing efforts to develop EMT inhibition as an approach for the treatment of retinal disease.

URINARY METABOLOMICS OF CENTRAL SEROUS CHORIORETINOPATHY.

PURPOSE: To analyze the urinary metabolomic profile of central serous chorioretinopathy cases. METHODS: In a cross-sectional study, 80 participants with central serous chorioretinopathy were compared with 80 age-matched and sex-matched controls. Urinary metabolites were measured using Metabolon's Discovery HD4 platform. RESULTS: Of 1,031 metabolites total that were measured in urine samples, 53 were upregulated and 27 downregulated in central serous chorioretinopathy participants compared with controls. After exclusion of potentially confounding xenobiotics and bile compounds that could represent digestive processes, 14 metabolites were significantly higher and 12 metabolites were significantly lower in cases compared with controls. One upregulated metabolite (tetrahydrocortisol sulfate) is involved in the corticosteroid subpathway. The downregulated metabolites are unrelated to the identified corticosteroid subpathway. CONCLUSION: The upregulation of urinary tetrahydrocortisol sulfate in central serous chorioretinopathy cases provides a precise molecular basis to further study the role of corticosteroids in producing choroidal venous congestion.

Recursive traffic percolation on urban transportation systems.

This paper proposes a recursive traffic percolation framework to capture the dynamics of cascading failures and analyze potential overloaded bottlenecks. In particular, compared to current work, the influence of external flow is considered, providing a new perspective for the study of regional commuting. Finally, we present an empirical study to verify the accuracy and effectiveness of our framework. Further analysis indicates that external flows from different regions affect the network. Our work requires only primary data and verifies the improvement of the functional network.

SUMOylation regulates Rb hyperphosphorylation and inactivation in uveal melanoma.

Small ubiquitin-like modifier (SUMO)ylation is one of the posttranslational modifications and is implicated in many tumor types. Modulation of SUMOylation can affect tumor progression, but the underlying mechanisms remain unclear. Here, we show that, for the first time, in uveal melanoma (UM), the most common intraocular malignancy in adults, global SUMOylation is upregulated and participates in tumor growth. Inhibition of SUMOylation in UM is sufficient to reduce tumor growth both in vitro and in vivo. Furthermore, we found that retinoblastoma protein (Rb) is a target protein and a critical downstream effector of the upregulated SUMOylation activity in UM. Increased SUMOylation of the Rb protein leads to its hyperphosphorylation and inactivation in UM cells, promoting UM cell proliferation. In summary, our results provide novel insight into the mechanism underlying SUMOylation-regulated tumor growth in UM.