Advances in lasso peptide discovery, biosynthesis, and function.

Lasso peptides are a large and sequence-diverse class of ribosomally synthesized and post-translationally modified peptide (RiPP) natural products characterized by their slip knot-like shape. These unique, highly stable peptides are produced by bacteria for various purposes. Their stability and sequence diversity make them a potentially useful scaffold for biomedically relevant folded peptides. However, many questions remain about lasso peptide biosynthesis, ecological function, and diversification potential for biomedical and agricultural applications. This review discusses new insights and open questions about lasso peptide biosynthesis and biological function. The role that genome mining has played in the development of new methodologies for discovering and diversifying lasso peptides is also discussed.

Power of inclusion: Enhancing polygenic prediction with admixed individuals.

Admixed individuals offer unique opportunities for addressing limited transferability in polygenic scores (PGSs), given the substantial trans-ancestry genetic correlation in many complex traits. However, they are rarely considered in PGS training, given the challenges in representing ancestry-matched linkage-disequilibrium reference panels for admixed individuals. Here we present inclusive PGS (iPGS), which captures ancestry-shared genetic effects by finding the exact solution for penalized regression on individual-level data and is thus naturally applicable to admixed individuals. We validate our approach in a simulation study across 33 configurations with varying heritability, polygenicity, and ancestry composition in the training set. When iPGS is applied to n = 237,055 ancestry-diverse individuals in the UK Biobank, it shows the greatest improvements in Africans by 48.9% on average across 60 quantitative traits and up to 50-fold improvements for some traits (neutrophil count, R2 = 0.058) over the baseline model trained on the same number of European individuals. When we allowed iPGS to use n = 284,661 individuals, we observed an average improvement of 60.8% for African, 11.6% for South Asian, 7.3% for non-British White, 4.8% for White British, and 17.8% for the other individuals. We further developed iPGS+refit to jointly model the ancestry-shared and -dependent genetic effects when heterogeneous genetic associations were present. For neutrophil count, for example, iPGS+refit showed the highest predictive performance in the African group (R2 = 0.115), which exceeds the best predictive performance for the White British group (R2 = 0.090 in the iPGS model), even though only 1.49% of individuals used in the iPGS training are of African ancestry. Our results indicate the power of including diverse individuals for developing more equitable PGS models.

Prediction of disease-free survival for precision medicine using cooperative learning on multi-omic data.

In precision medicine, both predicting the disease susceptibility of an individual and forecasting its disease-free survival are areas of key research. Besides the classical epidemiological predictor variables, data from multiple (omic) platforms are increasingly available. To integrate this wealth of information, we propose new methodology to combine both cooperative learning, a recent approach to leverage the predictive power of several datasets, and polygenic hazard score models. Polygenic hazard score models provide a practitioner with a more differentiated view of the predicted disease-free survival than the one given by merely a point estimate, for instance computed with a polygenic risk score. Our aim is to leverage the advantages of cooperative learning for the computation of polygenic hazard score models via Cox's proportional hazard model, thereby improving the prediction of the disease-free survival. In our experimental study, we apply our methodology to forecast the disease-free survival for Alzheimer's disease (AD) using three layers of data. One layer contains epidemiological variables such as sex, APOE (apolipoprotein E, a genetic risk factor for AD) status and 10 leading principal components. Another layer contains selected genomic loci, and the last layer contains methylation data for selected CpG sites. We demonstrate that the survival curves computed via cooperative learning yield an AUC of around $0.7$, above the state-of-the-art performance of its competitors. Importantly, the proposed methodology returns (1) a linear score that can be easily interpreted (in contrast to machine learning approaches), and (2) a weighting of the predictive power of the involved data layers, allowing for an assessment of the importance of each omic (or other) platform. Similarly to polygenic hazard score models, our methodology also allows one to compute individual survival curves for each patient.

Topological Twists in Nature.

The first entangled protein was observed about 30 years ago, resulting in an increased interest for uncovering the biological functions and biophysical properties of these complex topologies. Recently, the Pierced Lasso Topology (PLT) was discovered in which a covalent bond forms an intramolecular loop, leaving one or both termini free to pierce the loop. This topology is related to knots and other entanglements. PLTs exist in many well-researched systems where the PLTs have previously been unnoticed. PLTs represents 18% of all disulfide containing proteins across all kingdoms of life. In this review, we investigate the biological implications of this specific topology in which the PLT-forming disulfide may act as a molecular switch for protein function and consequently human health.

Polygenic hazard score models for the prediction of Alzheimer's free survival using the lasso for Cox's proportional hazards model.

The prediction of the susceptibility of an individual to a certain disease is an important and timely research area. An established technique is to estimate the risk of an individual with the help of an integrated risk model, that is, a polygenic risk score with added epidemiological covariates. However, integrated risk models do not capture any time dependence, and may provide a point estimate of the relative risk with respect to a reference population. The aim of this work is twofold. First, we explore and advocate the idea of predicting the time-dependent hazard and survival (defined as disease-free time) of an individual for the onset of a disease. This provides a practitioner with a much more differentiated view of absolute survival as a function of time. Second, to compute the time-dependent risk of an individual, we use published methodology to fit a Cox's proportional hazard model to data from a genetic SNP study of time to Alzheimer's disease (AD) onset, using the lasso to incorporate further epidemiological variables such as sex, APOE (apolipoprotein E, a genetic risk factor for AD) status, 10 leading principal components, and selected genomic loci. We apply the lasso for Cox's proportional hazards to a data set of 6792 AD patients (composed of 4102 cases and 2690 controls) and 87 covariates. We demonstrate that fitting a lasso model for Cox's proportional hazards allows one to obtain more accurate survival curves than with state-of-the-art (likelihood-based) methods. Moreover, the methodology allows one to obtain personalized survival curves for a patient, thus giving a much more differentiated view of the expected progression of a disease than the view offered by integrated risk models. The runtime to compute personalized survival curves is under a minute for the entire data set of AD patients, thus enabling it to handle datasets with 60,000-100,000 subjects in less than 1 h.

Multi-omics regulatory network inference in the presence of missing data.

A key problem in systems biology is the discovery of regulatory mechanisms that drive phenotypic behaviour of complex biological systems in the form of multi-level networks. Modern multi-omics profiling techniques probe these fundamental regulatory networks but are often hampered by experimental restrictions leading to missing data or partially measured omics types for subsets of individuals due to cost restrictions. In such scenarios, in which missing data is present, classical computational approaches to infer regulatory networks are limited. In recent years, approaches have been proposed to infer sparse regression models in the presence of missing information. Nevertheless, these methods have not been adopted for regulatory network inference yet. In this study, we integrated regression-based methods that can handle missingness into KiMONo, a Knowledge guided Multi-Omics Network inference approach, and benchmarked their performance on commonly encountered missing data scenarios in single- and multi-omics studies. Overall, two-step approaches that explicitly handle missingness performed best for a wide range of random- and block-missingness scenarios on imbalanced omics-layers dimensions, while methods implicitly handling missingness performed best on balanced omics-layers dimensions. Our results show that robust multi-omics network inference in the presence of missing data with KiMONo is feasible and thus allows users to leverage available multi-omics data to its full extent.

DBGRU-SE: predicting drug-drug interactions based on double BiGRU and squeeze-and-excitation attention mechanism.

The prediction of drug-drug interactions (DDIs) is essential for the development and repositioning of new drugs. Meanwhile, they play a vital role in the fields of biopharmaceuticals, disease diagnosis and pharmacological treatment. This article proposes a new method called DBGRU-SE for predicting DDIs. Firstly, FP3 fingerprints, MACCS fingerprints, Pubchem fingerprints and 1D and 2D molecular descriptors are used to extract the feature information of the drugs. Secondly, Group Lasso is used to remove redundant features. Then, SMOTE-ENN is applied to balance the data to obtain the best feature vectors. Finally, the best feature vectors are fed into the classifier combining BiGRU and squeeze-and-excitation (SE) attention mechanisms to predict DDIs. After applying five-fold cross-validation, The ACC values of DBGRU-SE model on the two datasets are 97.51 and 94.98%, and the AUC are 99.60 and 98.85%, respectively. The results showed that DBGRU-SE had good predictive performance for drug-drug interactions.

Bayesian LASSO for population stratification correction in rare haplotype association studies.

Population stratification (PS) is one major source of confounding in both single nucleotide polymorphism (SNP) and haplotype association studies. To address PS, principal component regression (PCR) and linear mixed model (LMM) are the current standards for SNP associations, which are also commonly borrowed for haplotype studies. However, the underfitting and overfitting problems introduced by PCR and LMM, respectively, have yet to be addressed. Furthermore, there have been only a few theoretical approaches proposed to address PS specifically for haplotypes. In this paper, we propose a new method under the Bayesian LASSO framework, QBLstrat, to account for PS in identifying rare and common haplotypes associated with a continuous trait of interest. QBLstrat utilizes a large number of principal components (PCs) with appropriate priors to sufficiently correct for PS, while shrinking the estimates of unassociated haplotypes and PCs. We compare the performance of QBLstrat with the Bayesian counterparts of PCR and LMM and a current method, haplo.stats. Extensive simulation studies and real data analyses show that QBLstrat is superior in controlling false positives while maintaining competitive power for identifying true positives under PS.

Comparing methods for statistical inference with model uncertainty.

Probability models are used for many statistical tasks, notably parameter estimation, interval estimation, inference about model parameters, point prediction, and interval prediction. Thus, choosing a statistical model and accounting for uncertainty about this choice are important parts of the scientific process. Here we focus on one such choice, that of variables to include in a linear regression model. Many methods have been proposed, including Bayesian and penalized likelihood methods, and it is unclear which one to use. We compared 21 of the most popular methods by carrying out an extensive set of simulation studies based closely on real datasets that span a range of situations encountered in practical data analysis. Three adaptive Bayesian model averaging (BMA) methods performed best across all statistical tasks. These used adaptive versions of Zellner's g-prior for the parameters, where the prior variance parameter g is a function of sample size or is estimated from the data. We found that for BMA methods implemented with Markov chain Monte Carlo, 10,000 iterations were enough. Computationally, we found two of the three best methods (BMA with g=√n and empirical Bayes-local) to be competitive with the least absolute shrinkage and selection operator (LASSO), which is often preferred as a variable selection technique because of its computational efficiency. BMA performed better than Bayesian model selection (in which just one model is selected).

Using phylogenetics to infer HIV-1 transmission direction between known transmission pairs.

Inferring the transmission direction between linked individuals living with HIV provides unparalleled power to understand the epidemiology that determines transmission. Phylogenetic ancestral-state reconstruction approaches infer the transmission direction by identifying the individual in whom the most recent common ancestor of the virus populations originated. While these methods vary in accuracy, it is unclear why. To evaluate the performance of phylogenetic ancestral-state reconstruction to determine the transmission direction of HIV-1 infection, we inferred the transmission direction for 112 transmission pairs where transmission direction and detailed additional information were available. We then fit a statistical model to evaluate the extent to which epidemiological, sampling, genetic, and phylogenetic factors influenced the outcome of the inference. Finally, we repeated the analysis under real-life conditions with only routinely available data. We found that whether ancestral-state reconstruction correctly infers the transmission direction depends principally on the phylogeny's topology. For example, under real-life conditions, the probability of identifying the correct transmission direction increases from 32%-when a monophyletic-monophyletic or paraphyletic-polyphyletic tree topology is observed and when the tip closest to the root does not agree with the state at the root-to 93% when a paraphyletic-monophyletic topology is observed and when the tip closest to the root agrees with the root state. Our results suggest that documenting larger differences in relative intrahost diversity increases our confidence in the transmission direction inference of linked pairs for population-level studies of HIV. These findings provide a practical starting point to determine our confidence in transmission direction inference from ancestral-state reconstruction.

CNVDeep: deep association of copy number variants with neurocognitive disorders.

BACKGROUND: Copy number variants (CNVs) have become increasingly instrumental in understanding the etiology of all diseases and phenotypes, including Neurocognitive Disorders (NDs). Among the well-established regions associated with ND are small parts of chromosome 16 deletions (16p11.2) and chromosome 15 duplications (15q3). Various methods have been developed to identify associations between CNVs and diseases of interest. The majority of methods are based on statistical inference techniques. However, due to the multi-dimensional nature of the features of the CNVs, these methods are still immature. The other aspect is that regions discovered by different methods are large, while the causative regions may be much smaller. RESULTS: In this study, we propose a regularized deep learning model to select causal regions for the target disease. With the help of the proximal [20] gradient descent algorithm, the model utilizes the group LASSO concept and embraces a deep learning model in a sparsity framework. We perform the CNV analysis for 74,811 individuals with three types of brain disorders, autism spectrum disorder (ASD), schizophrenia (SCZ), and developmental delay (DD), and also perform cumulative analysis to discover the regions that are common among the NDs. The brain expression of genes associated with diseases has increased by an average of 20 percent, and genes with homologs in mice that cause nervous system phenotypes have increased by 18 percent (on average). The DECIPHER data source also seeks other phenotypes connected to the detected regions alongside gene ontology analysis. The target diseases are correlated with some unexplored regions, such as deletions on 1q21.1 and 1q21.2 (for ASD), deletions on 20q12 (for SCZ), and duplications on 8p23.3 (for DD). Furthermore, our method is compared with other machine learning algorithms. CONCLUSIONS: Our model effectively identifies regions associated with phenotypic traits using regularized deep learning. Rather than attempting to analyze the whole genome, CNVDeep allows us to focus only on the causative regions of disease.

Exploring the folding landscape of leptin: Insights into threading pathways.

The discovery of new protein topologies with entanglements and loop-crossings have shown the impact of local amino acid arrangement and global three-dimensional structures. This phenomenon plays a crucial role in understanding how protein structure relates to folding and function, affecting the global stability, and biological activity. Protein entanglements encompassing knots and non-trivial topologies add complexity to their folding free energy landscapes. However, the initial native contacts driving the threading event for entangled proteins remains elusive. The Pierced Lasso Topology (PLT) represents an entangled topology where a covalent linker creates a loop in which the polypeptide backbone is threaded through. Compared to true knotted topologies, PLTs are simpler topologies where the covalent-loop persists in all conformations. In this work, the PLT protein leptin, is used to visualize and differentiate the preference for slipknotting over plugging transition pathways along the folding route. We utilize the Energy Landscape Visualization Method (ELViM), a multidimensional projection technique, to visualize and distinguish early threaded conformations that cannot be observed in an in vitro experiment. Critical contacts for the leptin threading mechanisms were identified where the competing pathways are determined by the formation of a hairpin loop in the unfolded basin. Thus, prohibiting the dominant slipknotting pathway. Furthermore, ELViM offers insights into distinct folding pathways associated with slipknotting and plugging providing a novel tool for de novo design and in vitro experiments with residue specific information of threading events in silico.

SpliceAPP: an interactive web server to predict splicing errors arising from human mutations.

BACKGROUND: Splicing variants are a major class of pathogenic mutations, with their severity equivalent to nonsense mutations. However, redundant and degenerate splicing signals hinder functional assessments of sequence variations within introns, particularly at branch sites. We have established a massively parallel splicing assay to assess the impact on splicing of 11,191 disease-relevant variants. Based on the experimental results, we then applied regression-based methods to identify factors determining splicing decisions and their respective weights. RESULTS: Our statistical modeling is highly sensitive, accurately annotating the splicing defects of near-exon intronic variants, outperforming state-of-the-art predictive tools. We have incorporated the algorithm and branchpoint information into a web-based tool, SpliceAPP, to provide an interactive application. This user-friendly website allows users to upload any genetic variants with genome coordinates (e.g., chr15 74,687,208 A G), and the tool will output predictions for splicing error scores and evaluate the impact on nearby splice sites. Additionally, users can query branch site information within the region of interest. CONCLUSIONS: In summary, SpliceAPP represents a pioneering approach to screening pathogenic intronic variants, contributing to the development of precision medicine. It also facilitates the annotation of splicing motifs. SpliceAPP is freely accessible using the link https://bc.imb.sinica.edu.tw/SpliceAPP . Source code can be downloaded at https://github.com/hsinnan75/SpliceAPP .

Glutamine metabolism-related genes and immunotherapy in nonspecific orbital inflammation were validated using bioinformatics and machine learning.

BACKGROUND: Nonspecific orbital inflammation (NSOI) is an idiopathic, persistent, and proliferative inflammatory condition affecting the orbit, characterized by polymorphous lymphoid infiltration. Its pathogenesis and progression have been linked to imbalances in tumor metabolic pathways, with glutamine (Gln) metabolism emerging as a critical aspect in cancer. Metabolic reprogramming is known to influence clinical outcomes in various malignancies. However, comprehensive research on glutamine metabolism's significance in NSOI is lacking. METHODS: This study conducted a bioinformatics analysis to identify and validate potential glutamine-related molecules (GlnMgs) associated with NSOI. The discovery of GlnMgs involved the intersection of differential expression analysis with a set of 42 candidate GlnMgs. The biological functions and pathways of the identified GlnMgs were analyzed using GSEA and GSVA. Lasso regression and SVM-RFE methods identified hub genes and assessed the diagnostic efficacy of fourteen GlnMgs in NSOI. The correlation between hub GlnMgs and clinical characteristics was also examined. The expression levels of the fourteen GlnMgs were validated using datasets GSE58331 and GSE105149. RESULTS: Fourteen GlnMgs related to NSOI were identified, including FTCD, CPS1, CTPS1, NAGS, DDAH2, PHGDH, GGT1, GCLM, GLUD1, ART4, AADAT, ASNSD1, SLC38A1, and GFPT2. Biological function analysis indicated their involvement in responses to extracellular stimulus, mitochondrial matrix, and lipid transport. The diagnostic performance of these GlnMgs in distinguishing NSOI showed promising results. CONCLUSIONS: This study successfully identified fourteen GlnMgs associated with NSOI, providing insights into potential novel biomarkers for NSOI and avenues for monitoring disease progression.

Targeted Learning with an Undersmoothed Lasso Propensity Score Model for Large-Scale Covariate Adjustment in Healthcare Database Studies.

Lasso regression is widely used for large-scale propensity score (PS) estimation in healthcare database studies. In these settings, previous work has shown that undersmoothing (overfitting) Lasso PS models can improve confounding control, but it can also cause problems of non-overlap in covariate distributions. It remains unclear how to select the degree of undersmoothing when fitting large-scale Lasso PS models to improve confounding control while avoiding issues that can result from reduced covariate overlap. Here, we used simulations to evaluate the performance of using collaborative-controlled targeted learning to data-adaptively select the degree of undersmoothing when fitting large-scale PS models within both singly and doubly robust frameworks to reduce bias in causal estimators. Simulations showed that collaborative learning can data-adaptively select the degree of undersmoothing to reduce bias in estimated treatment effects. Results further showed that when fitting undersmoothed Lasso PS-models, the use of cross-fitting was important for avoiding non-overlap in covariate distributions and reducing bias in causal estimates.

From bioinformatics to clinical applications: a novel prognostic model of cuproptosis-related genes based on single-cell RNA sequencing data in hepatocellular carcinoma.

OBJECTIVE AND METHODS: To ascertain the connection between cuproptosis-related genes (CRGs) and the prognosis of hepatocellular carcinoma (HCC) via single-cell RNA sequencing (scRNA-seq) and RNA sequencing (RNA-seq) data, relevant data were downloaded from the GEO and TCGA databases. The differentially expressed CRGs (DE-CRGs) were filtered by the overlaps in differentially expressed genes (DEGs) between HCC patients and normal controls (NCs) in the scRNA-seq database, DE-CRGs between high- and low-CRG-activity cells, and DEGs between HCC patients and NCs in the TCGA database. RESULTS: Thirty-three DE-CRGs in HCC were identified. A prognostic model (PM) was created employing six survival-related genes (SRGs) (NDRG2, CYB5A, SOX4, MYC, TM4SF1, and IFI27) via univariate Cox regression analysis and LASSO. The predictive ability of the model was validated via a nomogram and receiver operating characteristic curves. Research has employed tumor immune dysfunction and exclusion as a means to examine the influence of PM on immunological heterogeneity. Macrophage M0 levels were significantly different between the high-risk group (HRG) and the low-risk group (LRG), and a greater macrophage level was linked to a more unfavorable prognosis. The drug sensitivity data indicated a substantial difference in the half-maximal drug-suppressive concentrations of idarubicin and rapamycin between the HRG and the LRG. The model was verified by employing public datasets and our cohort at both the protein and mRNA levels. CONCLUSION: A PM using 6 SRGs (NDRG2, CYB5A, SOX4, MYC, TM4SF1, and IFI27) was developed via bioinformatics research. This model might provide a fresh perspective for assessing and managing HCC.

Peroxisome proliferator-activated receptors signature reveal the head and neck squamous cell carcinoma energy metabolism phenotype and clinical outcome.

BACKGROUND: Peroxisome proliferator activating receptors (PPARs) are important regulators of nuclear hormone receptor function, and they play a key role in biological processes such as lipid metabolism, inflammation and cell proliferation. However, their role in head and neck squamous cell carcinoma (HNSC) is unclear. METHODS: We used multiple datasets, including TCGA-HNSC, GSE41613, GSE139324, PRJEB23709 and IMVigor, to perform a comprehensive analysis of PPAR-related genes in HNSC. Single-cell sequencing data were preprocessed using Seurat packets, and intercellular communication was analyzed using CellChat packets. Functional enrichment analysis of PPAR-related genes was performed using ClusterProfile and GSEA. Prognostic models were constructed using LASSO and Cox regression models, and immunohistochemical analyses were performed using human protein mapping (The Human Protein Atlas). RESULTS: Our single-cell RNA sequencing analysis revealed distinct cell populations in HNSC, with T cells having the most significant transcriptome differences between tumors and normal tissues. The PPAR features were higher in most cell types in tumor tissues compared with normal tissues. We identified 17 PPAR-associated differentially expressed genes between tumors and normal tissues. A prognostic model based on seven PPAR-associated genes was constructed with high accuracy in predicting 1, 2 and 3 year survival in patients with HNSC. In addition, patients with a low risk score had a higher immune score and a higher proportion of T cells, CD8+ T cells and cytotoxic lymphocytes. They also showed higher immune checkpoint gene expression, suggesting that they might benefit from immunotherapy. PPAR-related genes were found to be closely related to energy metabolism. CONCLUSIONS: Our study provides a comprehensive understanding of the role of PPAR related genes in HNSC. The identified PPAR features and constructed prognostic models may serve as potential biomarkers for HNSC prognosis and treatment response. In addition, our study found that PPAR-related genes can differentiate energy metabolism and distinguish energy metabolic heterogeneity in HNSC, providing new insights into the molecular mechanisms of HNSC progression and therapeutic response.

The effects of noise exposure on hippocampal cognition in C57BL/6 mice via transcriptomics.

BACKGROUND: Noise is an important environmental stressor in the industrialized world and has received increasing attention in recent years. Although epidemiological research has extensively demonstrated the relationship between noise and cognitive impairment, the specific molecular mechanisms and targets remain to be fully explored and understood. METHODS: To address this issue, 5-month-old C57BL/6 mice were divided into two groups, with one group exposed to white noise at 98 dB. The effects of noise on cognition in mice were investigated through molecular biology and behavioral experiments. Subsequently, transcriptomic sequencing of the hippocampus in both groups of mice was performed and enrichment analysis of differentially expressed genes (DEGs) was conducted using KEGG and GO databases. Furthermore, LASSO analysis was used to further narrow down the relevant DEGs, followed by enrichment analysis of these genes using KEGG and GO databases. The DEGs were further validated by rt-qPCR. RESULTS: Following noise exposure, the hippocampus levels of inflammation-related factors increased, the phosphorylation of Tau protein increased, the postsynaptic density protein decreased, the number of Nissl bodies decreased, and cell shrinkage in the hippocampus increased. Moreover, the behavioral experiments manifest characteristics indicative of a decline in cognitive.A total of 472 DEGs were identified through transcriptomic analysis, and seven relevant genes were screened by the LASSO algorithm, which were further validated by PCR to confirm their consistency with the omics results. CONCLUSION: In conclusion, noise exposure affects cognitive function in mice through multiple pathways, and the omics results provide new evidence for the cognitive impairment induced by noise exposure.

A novel model for predicting prognosis and response to immunotherapy in nasopharyngeal carcinoma patients.

Blood-based biomarkers of immune checkpoint inhibitors (ICIs) response in patients with nasopharyngeal carcinoma (NPC) are lacking, so it is necessary to identify biomarkers to select NPC patients who will benefit most or least from ICIs. The absolute values of lymphocyte subpopulations, biochemical indexes, and blood routine tests were determined before ICIs-based treatments in the training cohort (n = 130). Then, the least absolute shrinkage and selection operator (Lasso) Cox regression analysis was developed to construct a prediction model. The performances of the prediction model were compared to TNM stage, treatment, and Epstein-Barr virus (EBV) DNA using the concordance index (C-index). Progression-free survival (PFS) was estimated by Kaplan-Meier (K-M) survival curve. Other 63 patients were used for validation cohort. The novel model composed of histologic subtypes, CD19+ B cells, natural killer (NK) cells, regulatory T cells, red blood cells (RBC), AST/ALT ratio (SLR), apolipoprotein B (Apo B), and lactic dehydrogenase (LDH). The C-index of this model was 0.784 in the training cohort and 0.735 in the validation cohort. K-M survival curve showed patients with high-risk scores had shorter PFS compared to the low-risk groups. For predicting immune therapy responses, the receiver operating characteristic (ROC), decision curve analysis (DCA), net reclassifcation improvement index (NRI) and integrated discrimination improvement index (IDI) of this model showed better predictive ability compared to EBV DNA. In this study, we constructed a novel model for prognostic prediction and immunotherapeutic response prediction in NPC patients, which may provide clinical assistance in selecting those patients who are likely to gain long-lasting clinical benefits to anti-PD-1 therapy.

Very short sleep duration reveals a proteomic fingerprint that is selectively associated with incident diabetes mellitus but not with incident coronary heart disease: a cohort study.

BACKGROUND: The molecular pathways linking short and long sleep duration with incident diabetes mellitus (iDM) and incident coronary heart disease (iCHD) are not known. We aimed to identify circulating protein patterns associated with sleep duration and test their impact on incident cardiometabolic disease. METHODS: We assessed sleep duration and measured 78 plasma proteins among 3336 participants aged 46-68 years, free from DM and CHD at baseline, and identified cases of iDM and iCHD using national registers. Incident events occurring in the first 3 years of follow-up were excluded from analyses. Tenfold cross-fit partialing-out lasso logistic regression adjusted for age and sex was used to identify proteins that significantly predicted sleep duration quintiles when compared with the referent quintile 3 (Q3). Predictive proteins were weighted and combined into proteomic scores (PS) for sleep duration Q1, Q2, Q4, and Q5. Combinations of PS were included in a linear regression model to identify the best predictors of habitual sleep duration. Cox proportional hazards regression models with sleep duration quintiles and sleep-predictive PS as the main exposures were related to iDM and iCHD after adjustment for known covariates. RESULTS: Sixteen unique proteomic markers, predominantly reflecting inflammation and apoptosis, predicted sleep duration quintiles. The combination of PSQ1 and PSQ5 best predicted sleep duration. Mean follow-up times for iDM (n = 522) and iCHD (n = 411) were 21.8 and 22.4 years, respectively. Compared with sleep duration Q3, all sleep duration quintiles were positively and significantly associated with iDM. Only sleep duration Q1 was positively and significantly associated with iCHD. Inclusion of PSQ1 and PSQ5 abrogated the association between sleep duration Q1 and iDM. Moreover, PSQ1 was significantly associated with iDM (HR = 1.27, 95% CI: 1.06-1.53). PSQ1 and PSQ5 were not associated with iCHD and did not markedly attenuate the association between sleep duration Q1 with iCHD. CONCLUSIONS: We here identify plasma proteomic fingerprints of sleep duration and suggest that PSQ1 could explain the association between very short sleep duration and incident DM.