Search

Abstract Nanoplasmas induced by intense laser fields have attracted enormous attention due to their accompanied spectacular physical phenomena which are vigorously expected by the community of science and industry. For instance, the energetic electrons and ions produced in laser-driven nanoplasmas are significant for the development of compact beam sources. Nevertheless, effective confinement on the propagating charged particles, which was realized through magnetic field modulation and target structure design in big facilities, are largely absent in the microscopic regime. Here we introduce a reliable scheme to provide control on the emission direction of protons generated from surface ionization in gold nanoparticles driven by intense femtosecond laser fields. The ionization level of the nanosystem provides us a knob to manipulate the characteristics of the collective proton emission. The most probable emission direction can be precisely steered by tuning the excitation strength of the laser pulses. This work opens new avenue for controlling the ion emission in nanoplasmas and can vigorously promote the fields such as development of on-chip beam sources at micro-/nano-scales.

Abstract 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; I 2 = 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 (I 2 = 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 (Pro state cancer A ge-based P heWAS; 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.

Abstract The development of an effective survival prediction tool is key for reducing colorectal cancer mortality. Here, we apply a three-stage study to devise a polygenic prognostic score (PPS) for stratifying colorectal cancer overall survival. Leveraging two cohorts of 3703 patients, we first perform a genome-wide survival association analysis to develop eight candidate PPSs. Further using an independent cohort with 470 patients, we identify the 287 variants-derived PPS (i.e., PPS287) achieving an optimal prediction performance [hazard ratio (HR) per SD = 1.99, P = 1.76 × 10−8], accompanied by additional tests in two external cohorts, with HRs per SD of 1.90 (P = 3.21 × 10−14; 543 patients) and 1.80 (P = 1.11 × 10−9; 713 patients). Notably, the detrimental impact of pathologic characteristics and genetic risk could be attenuated by a healthy lifestyle, yielding a 7.62% improvement in the 5-year overall survival rate. Therefore, our findings demonstrate the integrated contribution of pathologic characteristics, germline variants, and lifestyle exposure to the prognosis of colorectal cancer patients.

Exposure to persistent organic pollutants (POPs) may contribute to colorectal cancer risk, but the underlying mechanisms of crucial POPs exposure remain unclear. Hence, we systematically investigated the associations among POPs exposure, genetics and epigenetics and their effects on colorectal cancer. A case-control study was conducted in the Chinese population for detecting POPs levels. We measured the concentrations of 24 POPs in the plasma using gas chromatography-tandem mass spectrometry (GC–MS/MS) and evaluated the clinical significance of POPs by calculating the area under the receiver operating characteristic curve (AUC). To assess the associations between candidate genetic variants and colorectal cancer risk, unconditional logistic regression was used. Compared with healthy control individuals, individuals with colorectal cancer exhibited higher concentrations of the majority of POPs. Exposure to PCB153 was positively associated with colorectal cancer risk, and PCB153 demonstrated superior accuracy (AUC=0.72) for predicting colorectal cancer compared to other analytes. On PCB153-related genes, the rs67734009 C allele was significantly associated with reduced colorectal cancer risk and lower plasma levels of PCB153. Moreover, rs67734009 exhibited an expression quantitative trait locus (eQTL) effect on ESR1, of which the expression level was negatively related to PCB153 concentration. Mechanistically, the risk allele of rs67734009 increased ESR1 expression via miR-3492 binding and m6A modification. Collectively, this study sheds light on potential genetic and epigenetic mechanisms linking PCB153 exposure and colorectal cancer risk, thereby providing insight into the accurate protection against POPs exposure.

Summary: 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.

Environmental pollutants known as polycyclic aromatic hydrocarbons (PAHs) are produced through the incomplete combustion of organic material. While PAHs have been investigated as genotoxicants, they can also operate through nongenotoxic pathways in estrogen-dependent malignancies, such as breast, cervical and ovarian cancer. However, whether PAHs induce colorectal cancer (CRC) risk through estrogenic effects is still illusive. Here, we systematically investigated the abnormal expression and activation of estrogen receptor beta (ERβ) regulated by PAHs in CRC as well as the underlying mechanisms of ERβ-mediated CRC risk. Based on the 300 plasma samples from CRC patients and healthy controls detected by GC–MS/MS, we found that the plasma concentrations of benzo[a]pyrene (BaP) were significantly higher in CRC cases than in healthy controls, with significant estrogenic effects. Moreover, histone deacetylase 2 (HDAC2)-induced deacetylation of the promoter decreases ERβ expression, which is associated with poor overall survival and advanced tumor stage. The study also revealed that BaP and estradiol (E2) had different carcinogenic effects, with BaP promoting cell proliferation and inhibiting apoptosis, while E2 had the opposite effects. Additionally, this study mapped ERβ genomic binding regions by performing ChIP-seq and ATAC-seq and identified genetic variants of rs1411680 and its high linkage disequilibrium SNP rs6477937, which were significantly associated with CRC risk through meta-analysis of two independent Chinese population genome-wide association studies comprising 2,248 cases and 3,173 controls and then validation in a large-scale European population. By integrating data from functional genomics, we validated the regulatory effect of rs6477937 as an ERβ binding-disrupting SNP that mediated allele-specific expression of LINC02977 in a long-range chromosomal interaction manner, which was found to be highly expressed in CRC tissues. Overall, this study suggests that the different active effects on ERβ by PAHs and endogenous E2 may play a crucial role in the development and progression of CRC and highlights the potential of targeting ERβ and its downstream targets for CRC prevention and treatment.

Background: Environmental pollutant measurement is essential for accurate health risk assessment. However, the detection of humans’ internal exposure to pollutants is cost-intensive and consumes time and energy. Polygenic risk scores (PRSs) have been widely applied in genetic studies of complex trait diseases. It is important to construct a genetically relevant environmental surrogate for pollutant exposure and to explore its utility for disease prediction and risk assessment. Objectives: This study enrolled 714 individuals with complete genomic data and exposomic data on 22 plasma-persistent organic pollutants (POPs). Methods: We first conducted 22 POP genome-wide association studies (GWAS) and constructed the corresponding environmental pollutant-based PRS (EpPRS) by clumping and P value thresholding (C + T), lassosum, and PRS-CS methods. The best-fit EpPRS was chosen by its regression R2. An adverse outcome pathway (AOP) framework was developed to assess the effects of contaminants on candidate diseases. Furthermore, Mendelian randomization (MR) analysis was performed to explore the causal association between POPs and cancer risk. Results: The C + T method produced the best-performing EpPRSs for 7 PCBs and 4 PBDEs. EpPRSs replicated the correlations of environmental exposure measurements based on consistent patterns. The diagnostic performance of type 2 diabetes mellitus (T2DM) PRS was improved by the combined model of T2DM-EpPRS of PCB126/BDE153. Finally, the AKT1-mediated AOP framework illustrated that PCB126 and BDE153 may increase the risk of T2DM by decreasing AKT1 phosphorylation through the cGMP-PKG pathway and promoting abnormal glucose homeostasis. MR analysis showed that digestive system tumors, such as colorectal cancer and biliary tract cancer, are more sensitive to POP exposure. Conclusions: EpPRSs can serve as a proxy for assessing pollutant internal exposure. The application of the EpPRS to disease risk assessment can reveal the toxic pathway and mode of action linking exposure and disease in detail, providing a basis for the development of environmental pollutant control strategies.

Abstract Background Colorectal cancer is one of the most common malignant digestive tract tumors with a poor prognosis. RNA 5‐methylcytosine (m5C) is an important posttranscriptional widespread modification involved in many biological processes. However, the association between genetic variations of m5C modification genes and the prognostic value of colorectal cancer remains unclear. Methods We investigated the association between candidate single nucleotide polymorphisms (SNPs) in 13 m5C modification genes and colorectal cancer overall survival (OS) after chemotherapy by the Cox regression model. The combined effect of selected SNPs on OS, progression‐free survival (PFS), and disease control rate (DCR) was assessed by the number of risk alleles (NRA). The GTEx and TCGA database were used to perform expression qualitative trait locus (eQTL) analysis. Results We identified that two SNPs in YBX1 were associated with OS after chemotherapy (HR = 1.43, p = 0.001 for rs10890208; HR = 1.36, p = 0.025 for rs3862218). A striking dose–response effect between NRA and OS after chemotherapy was found (ptrend = 0.002). The DCR of patients receiving oxaliplatin chemotherapy in the 3–4 NRA group was markedly reduced in comparison to that in the 0–2 NRA group (OR = 1.49, p = 0.036). Moreover, YBX1 mRNA expression was significantly overexpressed in tumor tissues (p

Abstract Background The genetic architectures of colorectal cancer are distinct across different populations. To date, the majority of polygenic risk scores (PRSs) are derived from European (EUR) populations, which limits their accurate extrapolation to other populations. Here, we aimed to generate a PRS by incorporating East Asian (EAS) and EUR ancestry groups and validate its utility for colorectal cancer risk assessment among different populations. Methods A large-scale colorectal cancer genome-wide association study (GWAS), harboring 35,145 cases and 288,934 controls from EAS and EUR populations, was used for the EAS-EUR GWAS meta-analysis and the construction of candidate EAS-EUR PRSs via different approaches. The performance of each PRS was then validated in external GWAS datasets of EAS (727 cases and 1452 controls) and EUR (1289 cases and 1284 controls) ancestries, respectively. The optimal PRS was further tested using the UK Biobank longitudinal cohort of 355,543 individuals and ultimately applied to stratify individual risk attached by healthy lifestyle. Results In the meta-analysis across EAS and EUR populations, we identified 48 independent variants beyond genome-wide significance (P < 5 × 10−8) at previously reported loci. Among 26 candidate EAS-EUR PRSs, the PRS-CSx approach-derived PRS (defined as PRSCSx) that harbored genome-wide variants achieved the optimal discriminatory ability in both validation datasets, as well as better performance in the EAS population compared to the PRS derived from known variants. Using the UK Biobank cohort, we further validated a significant dose-response effect of PRSCSx on incident colorectal cancer, in which the risk was 2.11- and 3.88-fold higher in individuals with intermediate and high PRSCSx than in the low score subgroup (P trend = 8.15 × 10−53). Notably, the detrimental effect of being at a high genetic risk could be largely attenuated by adherence to a favorable lifestyle, with a 0.53% reduction in 5-year absolute risk. Conclusions In summary, we systemically constructed an EAS-EUR PRS to effectively stratify colorectal cancer risk, which highlighted its clinical implication among diverse ancestries. Importantly, these findings also supported that a healthy lifestyle could reduce the genetic impact on incident colorectal cancer.

Field-induced tunneling is one of the most fundamental quantum phenomena. In atoms this process has been successfully explored by attosecond interferometry based on high harmonic generation. Adapting this method to solids, the reconstruction of the subcycle tunneling dynamics calls for rigorous theoretical models that are able to properly map the experimental observables to the character of the tunneling process. Unlike in atomic gases, for crystalline solid-state systems the validity and applicability of the semiclassical trajectory-based model are still highly debated. Here we present a saddle-point analysis for solid-state systems which includes the quantum dynamics during tunneling. This allows us to quantify the initial conditions of electrons and holes when they emerge after the tunneling process in the classically allowed region. Our quantum trajectory simulations clarify the crucial role of the tunneling dynamics for the subsequent evolution and the harmonic emission from solids. Besides a nonzero initial electron-hole separation, a nonzero initial velocity of electrons/holes at the tunneling exits is revealed which arises from nonadiabatic tunneling. We find that depending on the ionization time, both inward and outward movements of the electron and its associated hole at the tunneling exit can occur. Our results provide intuitive insight into the nonadiabatic tunneling dynamics in solids and have direct implications for revealing fundamental quantum mechanical phenomena in solid-state systems with attosecond spectroscopic techniques.

Abstract 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, P IVW = 3.05×10-3) but decreased breast cancer risk (ORIVW = 0.68, P IVW = 8.19×10-3); however, the significance of breast cancer was dampened (P multivariable 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.

We theoretically investigate the non-sequential double ionization of Ar atoms in the combined fields of linearly polarized laser and circularly polarized laser through 3D semiclassical simulations. By partially overlapping the two time-delayed multicycle laser pulses, we construct an optical waveform whose polarization ellipticity increase slowly for consecutive optical cycles. This composite laser pulses with the time-dependent ellipticity can tunnel-ionize atoms and steer the first tunneling electron to recollision with the second bound electron through different trajectories, in which the recollision occurs with different return times of the first ionized electron. Through tuning delay time between the two laser pulses, the double ionization yields and recollision trajectories with different return times can be controlled. The time-dependent ellipticity with different delay time can enhance or suppress the probability of different return times. This work provides a scheme exploring electron dynamics in few optical cycle or even subcycle time scale in a multicycle laser field without having to be limited to near-single-cycle laser pulses.

Autophagy-related genes have a vital effect on colorectal cancer (CRC) by affecting genomic stability and regulating immune responses. However, the associations between genetic variants in autophagy-related genes and CRC outcomes for chemotherapy therapy remain unclear. The Cox regression model was used to evaluate the associations between single-nucleotide polymorphisms (SNPs) in autophagy-related genes and overall survival (OS) and progression-free survival (PFS) of CRC patients. The results were corrected by the false discovery rate (FDR) correction. We used the logistic regression model to investigate the associations of SNPs with the disease control rate (DCR) of patients. Gene expression analysis was explored based on an in-house dataset and other databases. The associations between gene expression and infiltrating immune cells were evaluated using the Tumor Immune Estimation Resource (TIMER) database. We observed that ATG2B rs17094017 A > T was significantly associated with increased OS (HR = 0.65, 95% CI = 0.50-0.86, P = 2.54×10-3), PFS (HR = 0.76, 95% CI = 0.62-0.93, P = 7.34×10-3), and DCR (OR = 0.60, 95% CI = 0.37-0.96, P = 3.31×10-2) of CRC patients after chemotherapy. The expression of ATG2B was down-expressed in CRC tissues than in adjacent normal tissues. Moreover, ATG2B expression influenced the infiltration of CD8+ T cells, CD4+ T cells, B cells, and T cell receptor signaling pathways, which may inhibit the occurrence of CRC by affecting the immune system. This study suggests that genetic variants in the autophagy-related gene ATG2B play a critical role in predicting the prognosis of CRC prognosis undergoing chemotherapy.

Abstract Small non-coding RNAs (ncRNAs) are vital regulators of biological activities, and aberrant levels of small ncRNAs are commonly found in precancerous lesions and cancer. PIWI-interacting RNAs (piRNAs) are a novel type of small ncRNA initially discovered in germ cells that have a specific length (24–31 nucleotides), bind to PIWI proteins, and show 2′-O-methyl modification at the 3′-end. Numerous studies have revealed that piRNAs can play important roles in tumorigenesis via multiple biological regulatory mechanisms, including silencing transcriptional and posttranscriptional gene processes and accelerating multiprotein interactions. piRNAs are emerging players in the malignant transformation of normal cells and participate in the regulation of cancer hallmarks. Most of the specific cancer hallmarks regulated by piRNAs are involved in sustaining proliferative signaling, resistance to cell death or apoptosis, and activation of invasion and metastasis. Additionally, piRNAs have been used as biomarkers for cancer diagnosis and prognosis and have great potential for clinical utility. However, research on the underlying mechanisms of piRNAs in cancer is limited. Here, we systematically reviewed recent advances in the biogenesis and biological functions of piRNAs and relevant bioinformatics databases with the aim of providing insights into cancer diagnosis and clinical applications. We also focused on some cancer hallmarks rarely reported to be related to piRNAs, which can promote in-depth research of piRNAs in molecular biology and facilitate their clinical translation into cancer treatment.

In two-color strong field ionization of atoms, dynamical characteristics such as ionization rate and electron trajectory depend heavily on the relative phase of the strong fundamental field and its second harmonic. The phase-of-the-phase spectroscopy method reveals the oscillation of the photoelectron momentum distribution with the relative phase of the two-color field, and the relative phase contrast and its phase are usually obtained by the fitting process. Instead, we apply the generalized Bessel function to the strong field approximation and derive analytical expressions for the first- and second-order relative phase contrast and phase-of-the-phase spectra. Our analytical results are in good agreement with fitting-based ones, but calculating them is much less computationally intensive. We illustrate the advantage of using generalized Bessel functions in two-color phase-of-the-phase spectroscopy of argon.

Abstract The association between air pollution and childhood respiratory disease is inconsistent. In the present study, we investigated a short-term effect of ambient air pollutants and daily childhood lower respiratory diseases (CLRD). Daily air pollutants, weather data, and CLRD data were collected from January 2014 to April 2015 (452 days) in Nanjing, China. Time-series regression and generalized additive models were used to assess the effects of air pollutants (PM10, PM2.5, NO2, SO2, O3, and CO) on CLRD. We observed that an interquartile range (IQR) increase in concentrations of PM10, NO2, and SO2 significantly increased the daily CLRD with 6 days cumulative effects (difference of estimates: 2.8%, 95% CI: 0.6–5.0%; 4.1%, 1.2–7.0%; 5.6%, 2.6–8.6%, respectively). However, no significant association was found in IQR concentrations of PM2.5, O3, and CO. Specifically, elevated PM10, PM2.5, NO2, and SO2 significantly increased the numbers of CLRD in cool season (3.6%, 1.5–5.7%; 2.4%, 0.3–4.5%; 4.9%, 2.9–7.0%; 6.3%, 3.7–9.0%, respectively). Additionally, the effect estimates of PM10, NO2, and SO2 in female and age >27 months were more pronounced than in male and age ≤27 months. This study suggested that short-term exposure to ambient PM10, NO2, and SO2 were associated with the increased CLRD numbers.

The current study aimed to investigate associations of circRNAs and related genetic variants with the risk of prostate cancer (PCa) as well as to elucidate biological mechanisms underlying the associations. We first compared expression levels of circRNAs between 25 paired PCa and adjacent normal tissues to identify riskassociated circRNAs by using the MiOncoCirc database. We then used logistic regression models to evaluate associations between genetic variants in candidate circRNAs and PCa risk among 4 662 prostate cancer patients and 3 114 healthy controls, and identified circHIBADH rs11973492 T>C as a significant risk-associated variant (odds ratio = 1.20, 95% confidence interval: 1.08–1.34, P = 7.06 × 10−4) in a dominant genetic model, which altered the secondary structure of the corresponding RNA chain. In the in silico analysis, we found that circHIBADH sponged and silenced 21 RNA-binding proteins (RBPs) enriched in the RNA splicing pathway, among which HNRNPA1 was identified and validated as a hub RBP using an external RNA-sequencing data as well as the in-house (four tissue samples) and publicly available single-cell transcriptomes. Additionally, we demonstrated that HNRNPA1 influenced hallmarks including MYC target, DNA repair, and E2F target signaling pathways, thereby promoting carcinogenesis. In conclusion, genetic variants in circHIBADH may act as sponges and inhibitors of RNA splicing-associated RBPs including HNRNPA1, playing an oncogenic role in PCa. [ABSTRACT FROM AUTHOR]

Genetic variants in super-enhancers (SEs) are increasingly implicated as a disease risk-driving mechanism. Previous studies have reported an associations between benzo[a]pyrene (BaP) exposure and some malignant tumor risk. Currently, it is unclear whether BaP is involved in the effect of genetic variants in SEs on prostate cancer risk, nor the associated intrinsic molecular mechanisms. In the current study, by using logistic regression analysis, we found that rs5750581T>C in 22q-SE was significantly associated with prostate cancer risk (odds ratio = 1.26, P = 7.61 × 10-5). We also have found that the rs6001092T>G, in a high linkage disequilibrium with rs5750581T>C (r² = 0.98), is located in a regulatory aryl hydrocarbon receptor (AhR) motif and may interact with the FAM227A promoter in further bioinformatics analysis. We then performed a series of functional and BaP acute exposure experiments to assess biological function of the genetic variant and the target gene. Biologically, the rs6001092-G allele strengthened the transcription factor binding affinity to AhR, thereby upregulating FAM227A, especially upon exposure to BaP, which induced the malignant phenotypes of prostate cancer. The current study highlights that AhR acts as an environmental sensor of BaP and is involved in the SE-mediated prostate cancer risk, which may provide new insights into the etiology of prostate cancer associated with the inherited SE variants under environmental carcinogen stressors. [ABSTRACT FROM AUTHOR]

The interplay between DNA replication stress and immune microenvironment alterations is known to play a crucial role in colorectal tumorigenesis, but a comprehensive understanding of their association with and relevant biomarkers involved in colorectal tumorigenesis is lacking. To address this gap, we conducted a study aiming to investigate this association and identify relevant biomarkers. We analyzed transcriptomic and proteomic profiles of 904 colorectal tumor tissues and 342 normal tissues to examine pathway enrichment, biological activity, and the immune microenvironment. Additionally, we evaluated genetic effects of single variants and genes on colorectal cancer susceptibility using data from genome-wide association studies (GWASs) involving both East Asian (7 062 cases and 195 745 controls) and European (24 476 cases and 23 073 controls) populations. We employed mediation analysis to infer the causal pathway, and applied multiplex immunofluorescence to visualize colocalized biomarkers in colorectal tumors and immune cells. Our findings revealed that both DNA replication activity and the flap structure-specific endonuclease 1 (FEN1) gene were significantly enriched in colorectal tumor tissues, compared with normal tissues. Moreover, a genetic variant rs4246215 G>T in FEN1 was associated with a decreased risk of colorectal cancer (odds ratio = 0.94, 95% confidence interval: 0.90-0.97, Pmeta = 4.70 × 10-9). Importantly, we identified basophils and eosinophils that both exhibited a significantly decreased infiltration in colorectal tumors, and were regulated by rs4246215 through causal pathways involving both FEN1 and DNA replication. In conclusion, this trans-omics incorporating GWAS data provides insights into a plausible pathway connecting DNA replication and immunity, expanding biological knowledge of colorectal tumorigenesis and therapeutic targets. [ABSTRACT FROM AUTHOR]