SUMMER: a Mendelian randomization interactive server to systematically evaluate the causal effects of risk factors and circulating biomarkers on pan-cancer survival.

Genome-wide association studies (GWASs) underlying case-control design have uncovered hundreds of genetic loci involved in tumorigenesis and provided rich resources for identifying risk factors and biomarkers associated with cancer susceptibility. However, the application of GWAS in determining the genetic architecture of cancer survival remains unestablished. Here, we systematically evaluated genetic effects at the genome-wide level on cancer survival that included overall survival (OS) and cancer-specific survival (CSS), leveraging data deposited in the UK Biobank cohort of a total of 19 628 incident patients across 17 cancer types. Furthermore, we assessed the causal effects of risk factors and circulating biomarkers on cancer prognosis via a Mendelian randomization (MR) analytic framework, which integrated cancer survival GWAS dataset, along with phenome-wide association study (PheWAS) and blood genome-wide gene expression/DNA methylation quantitative trait loci (eQTL/meQTL) datasets. On average, more than 10 traits, 700 genes, and 4,500 CpG sites were prone to cancer prognosis. Finally, we developed a user-friendly online database, SUrvival related cancer Multi-omics database via MEndelian Randomization (SUMMER; http://njmu-edu.cn:3838/SUMMER/), to help users query, browse, and download cancer survival results. In conclusion, SUMMER provides an important resource to assist the research community in understanding the genetic mechanisms of cancer survival.

Genetic variants in hypoxia-inducible factor pathway are associated with colorectal cancer risk and immune infiltration.

Hypoxia-inducible factor (HIF) pathway genes influence tumorigenesis and immune status. However, the associations between genetic variants in hypoxia-related genes and colorectal cancer risk and the immune status of hypoxia-associated genes in colorectal cancer have not been systematically characterized. The associations between genetic variants and colorectal cancer risk were evaluated in Chinese, Japanese and European populations using logistic regression analysis. The relationships between target genes and tumour immune infiltration were predicted by Tumour Immune Estimation Resource (TIMER). We found that rs34533650 in EPAS1 was associated with colorectal cancer risk (OR = 1.43, 95% CI = 1.20-1.70, P(FDR) = 8.35 × 10-4 ), and this finding was validated in two independent populations (Japanese: OR = 1.07, 95% CI = 1.01-1.15, p = 3.38 × 10-2 ; European: OR = 1.11, 95% CI = 1.03-1.19, p = 6.04 × 10-3 ). EPAS1-associated genes were enriched in immune-related pathways. In addition, we found that EPAS1 copy number variation (CNV) was associated with the degree of infiltration of immune cells and observed correlations between EPAS1 expression and immune cell infiltration levels in colorectal cancer. These results highlight that genetic variants of hypoxia-related genes play roles in colorectal cancer risk and provide new insight that EPAS1 might be a promising predictor of colorectal cancer susceptibility and immune status.

An early-onset specific polygenic risk score optimizes age-based risk estimate and stratification of prostate cancer: population-based cohort study.

BACKGROUND: Early-onset prostate cancer (EOPC, ≤ 55 years) has a unique clinical entity harboring high genetic risk, but the majority of EOPC patients still substantial opportunity to be early-detected thus suffering an unfavorable prognosis. A refined understanding of age-based polygenic risk score (PRS) for prostate cancer (PCa) would be essential for personalized risk stratification. METHODS: We included 167,517 male participants [4882 cases including 205 EOPC and 4677 late-onset PCa (LOPC)] from UK Biobank. A General-, an EOPC- and an LOPC-PRS were derived from age-specific genome-wide association studies. Weighted Cox proportional hazard models were applied to estimate the risk of PCa associated with PRSs. The discriminatory capability of PRSs were validated using time-dependent receiver operating characteristic (ROC) curves with additional 4238 males from PLCO and TCGA. Phenome-wide association studies underlying Mendelian Randomization were conducted to discover EOPC linking phenotypes. RESULTS: The 269-PRS calculated via well-established risk variants was more strongly associated with risk of EOPC [hazard ratio (HR) = 2.35, 95% confidence interval (CI) 1.99-2.78] than LOPC (HR = 1.95, 95% CI 1.89-2.01; I2 = 79%). EOPC-PRS was dramatically related to EOPC risk (HR = 4.70, 95% CI 3.98-5.54) but not to LOPC (HR = 0.98, 95% CI 0.96-1.01), while LOPC-PRS had similar risk estimates for EOPC and LOPC (I2 = 0%). Particularly, EOPC-PRS performed optimal discriminatory capability for EOPC (area under the ROC = 0.613). Among the phenomic factors to PCa deposited in the platform of ProAP (Prostate cancer Age-based PheWAS; https://mulongdu.shinyapps.io/proap ), EOPC was preferentially associated with PCa family history while LOPC was prone to environmental and lifestyles exposures. CONCLUSIONS: This study comprehensively profiled the distinct genetic and phenotypic architecture of EOPC. The EOPC-PRS may optimize risk estimate of PCa in young males, particularly those without family history, thus providing guidance for precision population stratification.

Metal Exposure Promotes Colorectal Tumorigenesis via the Aberrant N6-Methyladenosine Modification of ATP13A3.

Element contamination, including that from heavy metals, is associated with gastrointestinal tumorigenesis, but the effects and mechanisms of crucial element exposure associated with colorectal cancer remain unclear. We profiled 56 elements by ICP-MS and used logistic regression, LASSO, BKMR, and GAM to identify colorectal cancer-relevant elements. A series of biochemical experiments were performed to demonstrate the cytotoxicity and the mechanisms of malignant transformation after metal exposure. Using an elementomics approach, we first found that the metal thallium (Tl) was positively correlated with many toxic metals and was associated with a significantly increased risk of colorectal cancer. Acute exposure to Tl induced cytotoxicity and cell death by accelerating the generation of reactive oxygen species and DNA damage. Chronic exposure to Tl led to the inhibition of cell death and thereby induced the malignant transformation of normal colon cells and xenograft tumor formation in nude mice. Furthermore, we describe the first identification of a significant metal quantitative trait locus for the novel colorectal cancer susceptibility locus rs1511625 near ATP13A3. Mechanistically, Tl increased the level of aberrant N6-methyladenosine (m6A) modification of ATP13A3 via the METLL3/METTL14/ALKBH5-ATP13A3 axis to promote colorectal tumorigenesis. This study provides a basis for the development of public health strategies for reducing metal exposure among populations at high risk for colorectal cancer.

Systematic evaluation of the effects of genetic variants on PIWI-interacting RNA expression across 33 cancer types.

PIWI-interacting RNAs (piRNAs) are an emerging class of non-coding RNAs involved in tumorigenesis. Expression quantitative trait locus (eQTL) analysis has been demonstrated to help reveal the genetic mechanism of single nucleotide polymorphisms (SNPs) in cancer etiology. However, there are no databases that have been constructed to provide an eQTL analysis between SNPs and piRNA expression. In this study, we collected genotyping and piRNA expression data for 10 997 samples across 33 cancer types from The Cancer Genome Atlas (TCGA). Using linear regression cis-eQTL analysis with adjustment of appropriate covariates, we identified millions of SNP-piRNA pairs in tumor (76 924 831) and normal (24 431 061) tissues. Further, we performed differential expression and survival analyses, and linked the eQTLs to genome-wide association study (GWAS) data to comprehensively decipher the functional roles of identified cis-piRNA eQTLs. Finally, we developed a user-friendly database, piRNA-eQTL (http://njmu-edu.cn:3838/piRNA-eQTL/), to help users query, browse and download corresponding eQTL results. In summary, piRNA-eQTL could serve as an important resource to assist the research community in understanding the roles of genetic variants and piRNAs in the development of cancers.

Effect of polycyclic aromatic hydrocarbons on cancer risk causally mediated via vitamin D levels.

Polycyclic aromatic hydrocarbons (PAHs) widely exist in environmental substrates and are closely related to individual circulating vitamin D levels and tumorigenesis. Therefore, we proposed to evaluate the relationship between PAH exposure, vitamin D, and the risks for 14 cancer types via a causal inference framework underlying the mediation analysis. We evaluated seven urine monohydroxylated PAH (OH-PAH) and serum vitamin D concentrations of 3306 participants from the National Health and Nutrition Examination Survey between the 2013 and 2016 survey cycles and measured PAH concentrations in 150 subjects from the Nanjing cohort. We observed a significant negative dose-response relationship between increased OH-PAH levels and vitamin D deficiency. Each unit increase in ∑OH-PAHs could lead to a decrease in vitamin D levels (ßadj = -0.98, Padj = 2.05 × 10-4 ). Body mass index could have interaction effects with ∑OH-PAHs and affect vitamin D levels. Coexposure to naphthalene and fluorene metabolites mutually affected vitamin D levels. Notably, vitamin D could causally mediate the relationship between OH-PAHs and nine types of cancer (e.g., colorectal cancer, liver cancers, etc.). This study first emphasizes the causal cascade of individual OH-PAHs, vitamin D, and cancer risk, providing insights into prevention via the environment.

Association between circulating vitamin E and ten common cancers: evidence from large-scale Mendelian randomization analysis and a longitudinal cohort study.

BACKGROUND: The association between vitamin E and cancer risk has been widely investigated by observational studies, but the findings remain inconclusive. Here, we aimed to evaluate the causal effect of circulating vitamin E on the risk of ten common cancers, including bladder, breast, colorectal, esophagus, lung, oral and pharynx, ovarian, pancreatic, prostate, and kidney cancer. METHODS: A Mendelian randomization (MR) analytic framework was applied to data from a cancer-specific genome-wide association study (GWAS) comprising a total of 297,699 cancer cases and 304,736 controls of European ancestry. Three genetic instrumental variables associated with circulating vitamin E were selected. Summary statistic-based methods of inverse variance weighting (IVW) and likelihood-based approach, as well as the individual genotyping-based method of genetic risk score (GRS) were used. Multivariable IVW analysis was further performed to control for potential confounding effects. Furthermore, the UK Biobank cohort was used as external validation, supporting 355,543 European participants (incident cases ranged from 437 for ovarian cancer to 4882 for prostate cancer) for GRS-based estimation of circulating vitamin E, accompanied by a one-sample MR analysis of dietary vitamin E intake underlying the time-to-event analytic framework. RESULTS: Specific to cancer GWAS, we found that circulating vitamin E was significantly associated with increased bladder cancer risk (odds ratios [OR]IVW = 6.23, PIVW = 3.05×10-3) but decreased breast cancer risk (ORIVW = 0.68, PIVW = 8.19×10-3); however, the significance of breast cancer was dampened (Pmultivariable IVW > 0.05) in the subsequent multivariable MR analysis. In the validation stage of the UK Biobank cohort, we did not replicate convincing causal effects of genetically predicted circulating vitamin E concentrations and dietary vitamin E intake on the risk of ten cancers. CONCLUSIONS: This large-scale population study upon data from cancer-specific GWAS and a longitudinal biobank cohort indicates plausible non-causal associations between circulating vitamin E and ten common cancers in the European populations. Further studies regarding ancestral diversity are warranted to validate such causal associations.

CoSMeD: a user-friendly web server to estimate 5-year survival probability of left-sided and right-sided colorectal cancer patients using molecular data.

SUMMARY: Colorectal cancer is a heterogeneous disease with diverse prognoses between left-sided and right-sided patients; therefore, it is necessary to precisely evaluate the survival probability of side-specific colorectal cancer patients. Here, we collected multi-omics data from The Cancer Genome Atlas program, including gene expression, DNA methylation and microRNA expression. Specificity measure and robust likelihood-based survival analysis were used to identify 6 left-sided and 28 right-sided prognostic biomarkers. Compared to the performance of clinical prognostic models, the addition of these biomarkers could significantly improve the discriminatory ability and calibration in predicting side-specific 5-year survival for colorectal cancer. Additional dataset derived from Gene Expression Omnibus was used to validate the prognostic value of side-specific genes. Finally, we constructed colorectal cancer side-specific molecular database (CoSMeD), a user-friendly interface for estimating side-specific colorectal cancer 5-year survival probability, which can lay the basis for personalized management of left-sided and right-sided colorectal cancer patients. AVAILABILITY AND IMPLEMENTATION: CoSMeD is freely available at https://mulongdu.shinyapps.io/cosmed. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Photoelectron interference of He atoms in the attosecond ionization gating.

We theoretically investigate the photoelectron momentum distribution of He atoms by numerically solving the time-dependent Schro¨dinger equation (TDSE) in few-cycle ionization gating, which is synthesized by two linearly polarized laser pulses. When applying the TDSE, we can clearly see the spider-like structures in the photoelectron momentum spectra. We also find that the spider-like structures can be isolated by changing the relative phase. The directionality of the spider-like structure is changed from right-side to left-side and the ring-like interference structure gradually appears in the photoelectron momentum spectra when increasing the relative phase. The interference patterns observed in TDSE are recaptured well by the quantum-trajectory Monte Carlo (QTMC) model. We separate the ionization time window of the tunneling electron by analyzing the ionization rate. With the help of QTMC simulation, we illustrate the change of the interference structure and its directionality in the photoelectron momentum spectra. By changing the relative phase, the forward-backward asymmetry of the momentum distribution of the emitted electrons can also be controlled. Moreover, we find that the relative contribution of the nonrescattering and the rescattering trajectories can be controlled. These properties are beneficial for the application of photoelecron holography in probing atomic and molecular structures and dynamics.

Genetic variations in the CTLA-4 immune checkpoint pathway are associated with colon cancer risk, prognosis, and immune infiltration via regulation of IQCB1 expression.

The programmed cell death-1 (PD-1)/cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) immune checkpoint pathways serve as targets of immunotherapy for colorectal cancer. However, the associations between genetic variations in these pathways and colorectal cancer risk, prognosis, and immune status remain unclear. The associations between single-nucleotide polymorphisms (SNPs) and colorectal cancer risk and survival were evaluated in a case-control study comprising 1150 cases and 1342 controls along with 287 cases with overall survival information. We found that individuals with the A allele of B7-2 rs2681416 in CTLA-4 immune checkpoint pathway had a significantly increased risk of colorectal cancer [odds ratio (OR) = 1.37, P = 3.17 × 10-4] than those with G allele under the dominant model, which had a predominant site-specific effect in colon cancer (OR = 1.55, P = 3.11 × 10-5). In addition, rs2681416 significantly decreased the overall survival time of patients with colon cancer [hazard ratio (HR) = 1.96, P = 1.10 × 10-2], but not of patients with rectal cancer (P = 0.271). Moreover, rs2681416 had an expression quantitative trait locus effect on the B7-2 flanking gene IQCB1 in colon tissues, which contributed to colon cancer risk by regulating genome organization and influenced the expression of IQCB1 in an allele-specific manner. IQCB1 expression was upregulated in colorectal cancer tissues compared with normal tissues, accounting for various critical carcinogenic states in colon cancer and promoting immune infiltration of Th17 cells in the tumor microenvironment. Our study highlights the important roles of genetic variations in immune checkpoint pathways and provides new insight into potential site-specific independent biomarkers for colorectal cancer susceptibility, prognosis, and tumor immune status.

Evaluation of common genetic variants in vitamin E-related pathway genes and colorectal cancer susceptibility.

Vitamin E is effective for preventing the risk of cancer. However, few studies have elucidated the mechanism of vitamin E in cancer occurrence. Herein, we aimed to identify the genetic variants in vitamin E-related pathway genes associated with colorectal cancer risk. We applied logistic regression models to assess the association between single-nucleotide polymorphisms (SNPs) in vitamin E-related pathway genes and colorectal cancer risk in the Chinese and European population. The false discovery rate (FDR) method was used to correct multiple comparisons. The mRNA and protein expression analysis were evaluated in public database and in-house RNA-Seq data. SCARB1 rs73227586 was identified significantly increased risk of colorectal cancer in the Chinese population (odd ratio (OR) = 1.46, 95% confidence interval (CI) = 1.22-1.75, P = 2.99 × 10-5). This finding was further validated in the European population (OR = 1.11, 95% CI = 1.02-1.20, P = 1.44 × 10-2). Additionally, the mRNA and protein expression of SCARB1 were markedly up-regulated in colorectal tumor tissues. Moreover, rs73227586 T allele could increase the minimum free energy (MFE) and weaken binding ability to transcription factor ELL2. Our findings indicated that SCARB1 may play a carcinogenic role in colorectal cancer. Genetic variants in vitamin E-related pathway genes may concern to be predictors of colorectal cancer risk.

Investigation of electron vortices in time-delayed circularly polarized laser pulses with a semiclassical perspective.

We theoretically investigate strong-filed electron vortices in time-delayed circularly polarized laser pulses by a generalized quantum-trajectory Monte Carlo (GQTMC) model. Vortex interference patterns in photoelectron momentum distributions (PMDs) with various laser parameters can be well reproduced by the semiclassical simulation. The phase difference responsible for the interference structures is analytically identified through trajectory-based analysis and simple-man theory, which reveal the underlying mechanism of electron vortex phenomena for both co-rotating and counter-rotating component. This semiclassical analysis can also demonstrate the influences of laser intensity and wavelength on the number of arms of vortices. Furthermore, we show the influence of the Coulomb effect on the PMDs. Finally, the controlling of the ionization time intervals in the tens to hundreds of attosecond magnitude is qualitatively discussed.

Genetic variants in circTUBB interacting with smoking can enhance colorectal cancer risk.

Circular RNAs (circRNAs) are an intriguing class of regulatory RNAs involved in the tumorigenesis of many cancers, including colorectal cancer. Mechanistically, circRNAs sponge microRNAs (miRNAs) with specific miRNA response elements (MREs) and compete for regulatory target genes downstream. However, the genetic effects of MREs on colorectal cancer susceptibility remain unclear. The MREs of colorectal cancer-associated circRNAs (CRC-circRNAs) and corresponding single nucleotide polymorphisms (SNPs) were identified by the transcriptome of cancer cells and the 1000 Genomes Project. The association between candidate SNPs and colorectal cancer risk was evaluated in a Chinese population (1150 cases and 1342 controls) and two European populations (9023 cases and 386,896 controls) using logistic regression analysis. Among the 197 candidate SNPs within MREs of 186 CRC-circRNAs, rs25497 in circTUBB was significantly associated with colorectal cancer risk in a Chinese population after false discovery rate (FDR) correction [odds ratio (OR) = 1.78, 95% confidence interval (CI) = 1.44-2.21, P = 1.42 × 10-7, PFDR = 2.80 × 10-5] and even reached significance in a genome-wide association study (GWAS) under the dominant model (P = 1.28 × 10-8). Similar results were found in the European populations (ORmeta = 1.30, 95% CI = 1.10-1.53). Both stratification and interaction analyses revealed that rs25497 interacting with smoking affected the colorectal cancer risk (Pinteraction = 1.48 × 10-2). Here, we first comprehensively identified genetic variants in MREs of CRC-circRNAs and evaluated their effects on colorectal cancer risk in both Chinese and European populations. These findings provide basis for a comprehensive understanding of colorectal cancer etiology.

Vitamin B2 intake reduces the risk for colorectal cancer: a dose-response analysis.

PURPOSE: Several epidemiological studies have assessed the ability of vitamin B2 to prevent colorectal cancer (CRC), but the results are controversial results. We conducted a dose-response meta-analysis to investigate the association between vitamin B2 and CRC risk. METHODS: We searched the PubMed and EMBASE database until January 3, 2018 to identify relevant studies. The pooled relative risks (RRs) with the corresponding 95% confidence intervals (CIs) were calculated using a random-effects model or fixed-effects model. The dose-response relationship was assessed by restricted cubic splines. RESULTS: A total of 14 studies reporting vitamin B2 intake and two studies reporting blood vitamin B2 concentration, comprising 14,934 cases and 1593 cases, respectively, were included in the meta-analysis. Vitamin B2 intake was inversely associated with CRC risk (RR = 0.87; 95% CI 0.81-0.93). Similar results were found for total vitamin B2 intake from diet and supplements (RR = 0.86; 95% CI 0.78-0.94) and dietary vitamin B2 intake (RR = 0.89, 95% CI 0.82-0.98) in subgroup analyses. The dose-response model indicated a non-linear trend, and CRC risk was reduced by 10% when vitamin B2 intake increased to 5 mg/day. In addition, high blood concentrations of vitamin B2 could also reduce the CRC risk (RR = 0.74; 95% CI 0.59-0.92). CONCLUSIONS: This dose-response analysis indicates that vitamin B2 intake is inversely associated with CRC risk. The inverse association may also exist between blood vitamin B2 concentration and CRC risk. These results suggest the importance of vitamin B2 intake in the prevention of CRC.

Multiple recollisions in nonsequential double ionization by counter-rotating two-color circularly polarized laser fields.

With the three-dimensional (3D) classical ensemble method, we theoretically investigate the recollision dynamics in strong-field nonsequential double ionization (NSDI) of Ar by counter-rotating two-color circularly polarized laser fields. With the analysis of the NSDI trajectories, we find that not only multiple-recollision but also single-recollision processes occur in the double ionization events. Furthermore, the multiple-recollision and single-recollision processes both undergo the recollision-induced excitation with subsequent ionization (RESI) and recollision-induced ionization (RII). The angle between the momentum and the force of the laser field at the recollision moment can affect the times of the recollision.

Circadian clock pathway genes associated with colorectal cancer risk and prognosis.

Circadian clock genes influence biological processes and may be involved in tumorigenesis. We systematically evaluated genetic variants in the circadian clock pathway genes associated with colorectal cancer risk and survival. We evaluated the association of 119 single nucleotide polymorphisms (SNPs) in 27 circadian clock pathway genes with the risk of colorectal cancer in a case-control study (1150 cases and 1342 controls). The false discovery rate (FDR) method was applied to correct for multiple comparisons. Gene-based analysis was performed by the sequence kernel association test (SKAT). Cox proportional hazards regression was used to calculate the effects of SNPs on the overall survival of patients. We identified that compared to those with the G allele, individuals with the rs76436997 A allele in RORA had a significant 1.33-fold increased risk of colorectal cancer (P = 3.83 × 10- 4). Specifically, the GA/AA genotypes were related to an enhanced risk of colorectal cancer compared with that associated with the GG genotype, which was more common in patients with well and moderately differentiated tumors and Dukes A/B stages. The SNP rs76436997 significantly increased the overall survival time of colorectal cancer patients (P = 0.044). Furthermore, RNA-seq data showed that the mRNA levels of RORA were significantly lower in colorectal tumors than the paired normal tissues. Gene-based analysis revealed a significant association between RORA and colorectal cancer risk. These findings highlight the important roles of genetic variations in circadian clock pathway genes play in colorectal cancer risk and suggest that RORA is potentially related to colorectal carcinogenesis.

Nonsequential double ionization of Mg from a doubly excited complex driven by circularly polarized laser field.

With the classical ensemble method, the correlated-electron dynamics of Mg atom from a doubly excited, transition Coulomb complex in few-cycle circularly polarized (CP) laser field at low laser intensity is theoretically investigated. The low energy transfer during the recollision process indicates that the two electrons cannot release directly, but it can pass through a doubly excited state, and then escape with the ionization time difference. The numerical results show that the feature of the sequential double ionization (SDI) can be observed in the nonsequential double ionization (NSDI) process. The SDI-like results demonstrate that the intermediate state has lost any memory of its formation dynamics. The distribution of the angle between the two release directions of the two electrons also depends on the ionization time difference. Finally, the influence of e-e Coulomb repulsion is discussed.

Nonsequential double ionization channels control of Ar with few-cycle elliptically polarized laser pulse by carrier-envelope-phase.

The carrier-envelop-phase (CEP) dependence of nonsequential double ionization (NSDI) of atomic Ar with few-cycle elliptically polarized laser pulse is investigated using 2D classical ensemble method. We distinguish two particular recollision channels in NSDI, which are recollision-impact ionization (RII) and recollision-induced excitation with subsequent ionization (RESI). We separate the RII and RESI channels according to the delay time between recollision and final double ionization. By tracing the history of the trajectories, we find the electron correlation spectra as well as the competition between the two channels are sensitively dependent on the laser field CEP. Finally, control can be achieved between the two channels by varying the CEP.

Microglia-endothelial cross-talk regulates diabetes-induced retinal vascular dysfunction through remodeling inflammatory microenvironment.

Inflammation-mediated crosstalk between neuroglial cells and endothelial cells (ECs) is a fundamental feature of many vascular diseases. Nevertheless, the landscape of inflammatory processes during diabetes-induced microvascular dysfunction remains elusive. Here, we applied single-cell RNA sequencing to elucidate the transcriptional landscape of diabetic retinopathy (DR). The transcriptome characteristics of microglia and ECs revealed two microglial subpopulations and three EC populations. Exploration of intercellular crosstalk between microglia and ECs showed that diabetes-induced interactions mainly participated in the inflammatory response and vessel development, with colony-stimulating factor 1 (CSF1) and CSF1 receptor (CSF1R) playing important roles in early cell differentiation. Clinically, we found that CSF1/CSF1R crosstalk dysregulation was associated with proliferative DR. Mechanistically, ECs secrete CSF1 and activate CSF1R endocytosis and the CSF1R phosphorylation-mediated MAPK signaling pathway, which elicits the differentiation of microglia and triggers the secretion of inflammatory factors, and subsequently foster angiogenesis by remodeling the inflammatory microenvironment through a positive feedback mechanism.