Search

Primary Sjögren’s Syndrome (pSS) is a complex autoimmune disorder characterized by exocrine gland dysfunction, leading to dry eyes and mouth. Despite growing interest in biologic therapies for pSS, FDA approval has proven challenging due to trial complications. This review addresses the absence of a molecular-target-based approach to biologic therapy development and highlights novel research on drug targets and clinical trials. A literature search identified potential pSS treatment targets and recent advances in molecular understanding. Overlooking extraglandular symptoms like fatigue and depression is a notable gap in trials. Emerging biologic agents targeting cytokines, signal pathways, and immune responses have proven efficacy. These novel therapies could complement existing methods for symptom alleviation. Improved grading systems accounting for extraglandular symptoms are needed. The future of pSS treatment may involve gene, stem-cell, and tissue-engineering therapies. This narrative review offers insights into advancing pSS management through innovative biologic interventions.

Identifying essential proteins is of utmost importance in the field of biomedical research due to their essential functions in cellular activities and their involvement in mechanisms related to diseases. In this research, a novel approach called AttentionEP for predicting essential proteins (EP) is introduced by attention mechanisms. This method leverages both cross-attention and self-attention frameworks, focusing on enhancing prediction accuracy through the integration of features across diverse scales. Spatial characteristics of proteins are obtained from the protein-protein interaction (PPI) network by employing Graph Convolutional Networks (GCN) and Graph Attention Networks (GAT). Following this, Bidirectional Long Short-Term Memory networks (BiLSTM) are employed to derive temporal features from gene expression datasets. Furthermore, spatial characteristics are derived by integrating data on subcellular localization with the application of Deep Neural Networks (DNN). In order to effectively integrate features across multiple scales, initial steps involve the application of self-attention techniques to derive essential insights from each unique data set. Following this, mechanisms involving self-attention and cross-attention are employed to enhance the interaction between diverse information sources. To identify essential proteins, a classifier based on the ResNet architecture is developed. The findings from the experiments indicate that the method introduced here shows superior performance in identifying essential proteins, recording an Area Under the Curve (AUC) value of 0.9433. This approach shows a considerable advantage over established techniques. The findings of this study provide a significant advancement in the comprehension of critical proteins, revealing promising potential for applications in the development of therapeutics and addressing various diseases.

Abstract Numerical methods are widely used to calculate the secure key rate of many quantum key distribution protocols in practice, but they consume many computing resources and are too time-consuming. In this work, we take the homodyne detection discrete-modulated continuous-variable quantum key distribution (CV-QKD) as an example, and construct a neural network that can quickly predict the secure key rate based on the experimental parameters and experimental results. Compared to traditional numerical methods, the speed of the neural network is improved by several orders of magnitude. Importantly, the predicted key rates are not only highly accurate but also highly likely to be secure. This allows the secure key rate of discrete-modulated CV-QKD to be extracted in real time on a low-power platform. Furthermore, our method is versatile and can be extended to quickly calculate the complex secure key rates of various other unstructured quantum key distribution protocols.

This paper describes an ensemble method with supervised module detection and further module prioritization for reliable network-based biomarker discovery. We design a module detection and ranking method called mRank to discover reliable network modules as cancer diagnostic biomarkers, with two procedures: (1) an iterative supervised module detection guided by phenotypic states in a specific network, (2) a block-based module ranking locally and globally via network topological centrality. We validate its effectiveness and efficiency by identifying hepatocellular carcinoma (HCC) network modules on a comprehensive gene regulatory network with specifying gene interactions by HCC RNA-seq data from the Cancer Genome Atlas (TCGA). These top-ranked modules by mRank get a mean AUC of 0.995 on TCGA HCC dataset with 371 tumor samples and 50 controls by cross-validation SVM. Based on the prior knowledge of cancer dysfunctions enriched in top-ranked modules, 69 genes are identified as HCC candidate biomarkers. They are further validated in independent cohorts with a classifier trained on TCGA HCC dataset. A mean AUC of 0.846 is achieved in distinguishing 976 disease samples from 827 controls. Moreover, some known HCC signatures such as AFP and SPP1 are also included in our identified biomarkers. mRank enables us to find more reliable network modules for cancer diagnosis. For a proof-of-concept study, we validate it in identifying HCC network biomarkers and it is generalizable to other cancers or complex disease. The overall results have demonstrated that mRank can find effective network biomarkers for cancer diagnosis which result in less false positives.

Abstract Background Mining gene regulatory network (GRN) is an important avenue for addressing cancer mechanism. Mutations in cancer genome perturb GRN and cause a rewiring in an orchestrated network. Hence, the exploration of gene regulatory network rewiring is significant to discover potential biomarkers and indicators for discriminating cancer phenotypes. Results Here, we propose a new bioinformatics method of identifying biomarkers based on network rewiring in different states. It firstly reconstructs GRN in different phenotypic conditions from gene expression data with a priori background network. We employ the algorithm based on path consistency algorithm and conditional mutual information to delete false-positive regulatory interactions between independent nodes/genes or not closely related gene pairs. And then a differential gene regulatory network (D-GRN) is constructed from the rewiring parts in the two phenotype-specific GRNs. Community detection technique is then applied for D-GRN to detect functional modules. Finally, we apply logistic regression classifier with recursive feature elimination to select biomarker genes in each module individually. The extracted feature genes result in a gene set of biomarkers with impressing ability to distinguish normal samples from controls. We verify the identified biomarkers in external independent validation datasets. For a proof-of-concept study, we apply the framework to identify diagnostic biomarkers of breast cancer. The identified biomarkers obtain a maximum AUC of 0.985 in the internal sample classification experiments. And these biomarkers achieve a maximum AUC of 0.989 in the external validations. Conclusion In conclusion, network rewiring reveals significant differences between different phenotypes, which indicating cancer dysfunctional mechanisms. With the development of sequencing technology, the amount and quality of gene expression data become available. Condition-specific gene regulatory networks that are close to the real regulations in different states will be established. Revealing the network rewiring will greatly benefit the discovery of biomarkers or signatures for phenotypes. D-GRN is a general method to meet this demand of deciphering the high-throughput data for biomarker discovery. It is also easy to be extended for identifying biomarkers of other complex diseases beyond breast cancer.

Abstract Background The successful identification of breast cancer (BRCA) prognostic biomarkers is essential for the strategic interference of BRCA patients. Recently, various methods have been proposed for exploring a small prognostic gene set that can distinguish the high-risk group from the low-risk group. Methods Regularized Cox proportional hazards (RCPH) models were proposed to discover prognostic biomarkers of BRCA from gene expression data. Firstly, the maximum connected network with 1142 genes by mapping 956 differentially expressed genes (DEGs) and 677 previously BRCA-related genes into the gene regulatory network (GRN) was constructed. Then, the 72 union genes of the four feature gene sets identified by Lasso-RCPH, Enet-RCPH, $$L_{0}$$ L 0 -RCPH and SCAD-RCPH models were recognized as the robust prognostic biomarkers. These biomarkers were validated by literature checks, BRCA-specific GRN and functional enrichment analysis. Finally, an index of prognostic risk score (PRS) for BRCA was established based on univariate and multivariate Cox regression analysis. Survival analysis was performed to investigate the PRS on 1080 BRCA patients from the internal validation. Particularly, the nomogram was constructed to express the relationship between PRS and other clinical information on the discovery dataset. The PRS was also verified on 1848 BRCA patients of ten external validation datasets or collected cohorts. Results The nomogram highlighted that the importance of PRS in guiding significance for the prognosis of BRCA patients. In addition, the PRS of 301 normal samples and 306 tumor samples from five independent datasets showed that it is significantly higher in tumors than in normal tissues ( $$P

Abstract Background Hepatocellular carcinoma (HCC) is one of the most common cancers. The discovery of specific genes severing as biomarkers is of paramount significance for cancer diagnosis and prognosis. The high-throughput omics data generated by the cancer genome atlas (TCGA) consortium provides a valuable resource for the discovery of HCC biomarker genes. Numerous methods have been proposed to select cancer biomarkers. However, these methods have not investigated the robustness of identification with different feature selection techniques. Methods We use six different recursive feature elimination methods to select the gene signiatures of HCC from TCGA liver cancer data. The genes shared in the six selected subsets are proposed as robust biomarkers. Akaike information criterion (AIC) is employed to explain the optimization process of feature selection, which provides a statistical interpretation for the feature selection in machine learning methods. And we use several methods to validate the screened biomarkers. Results In this paper, we propose a robust method for discovering biomarker genes for HCC from gene expression data. Specifically, we implement recursive feature elimination cross-validation (RFE-CV) methods based on six different classication algorithms. The overlaps in the discovered gene sets via different methods are referred as the identified biomarkers. We give an interpretation of the feature selection process based on machine learning using AIC in statistics. Furthermore, the features selected by the backward logistic stepwise regression via AIC minimum theory are completely contained in the identified biomarkers. Through the classification results, the superiority of interpretable robust biomarker discovery method is verified. Conclusions It is found that overlaps among gene subsets contain different quantitative features selected by the RFE-CV of 6 classifiers. The AIC values in the model selection provide a theoretical foundation for the feature selection process of biomarker discovery via machine learning. What’s more, genes containing in more optimally selected subsets make better biological sense and implication. The quality of feature selection is improved by the intersections of biomarkers selected from different classifiers. This is a general method suitable for screening biomarkers of complex diseases from high-throughput data.

Gene regulatory network (GRN) provides abundant information on gene interactions, which contributes to demonstrating pathology, predicting clinical outcomes, and identifying drug targets. Existing high-throughput experiments provide rich time-series gene expression data to reconstruct the GRN to further gain insights into the mechanism of organisms responding to external stimuli. Numerous machine-learning methods have been proposed to infer gene regulatory networks. Nevertheless, machine learning, especially deep learning, is generally a “black box,” which lacks interpretability. The causality has not been well recognized in GRN inference procedures. In this article, we introduce grey theory integrated with the adaptive sliding window technique to flexibly capture instant gene–gene interactions in the uncertain regulatory system. Then, we incorporate generalized multivariate Granger causality regression methods to transform the dynamic grey association into causation to generate directional regulatory links. We evaluate our model on the DREAM4 in silico benchmark dataset and real-world hepatocellular carcinoma (HCC) time-series data. We achieved competitive results on the DREAM4 compared with other state-of-the-art algorithms and gained meaningful GRN structure on HCC data respectively.

Background We showed that Beclin‐1‐dependent autophagy protects the heart in young and adult mice that underwent endotoxemia. Herein, we compared the potential therapeutic effects of Beclin‐1 activating peptide, TB‐peptide, on endotoxemia‐induced cardiac outcomes in young adult and aged mice. We further evaluated lipopolysaccharide (lipopolysaccharide)‐induced and TB‐peptide treatment‐mediated alterations in myocardial metabolism. Methods and Results C57BL/6J mice that were 10 weeks and 24 months old were challenged by lipopolysaccharide using doses at which cardiac dysfunction occurred. Following the treatment of TB‐peptide or control vehicle, heart contractility, circulating cytokines, and myocardial autophagy were evaluated. We detected that TB‐peptide boosted autophagy, attenuated cytokines, and improved cardiac performance in both young and aged mice during endotoxemia. A targeted metabolomics assay was designed to detect a pool of 361 known metabolites, of which 156 were detected in at least 1 of the heart tissue samples. Lipopolysaccharide‐induced impairments were found in glucose and amino acid metabolisms in mice of all ages, and TB‐peptide ameliorated these alterations. However, lipid metabolites were upregulated in the young group but moderately downregulated in the aged by lipopolysaccharide, suggesting an age‐dependent response. TB‐peptide mitigated lipopolysaccharide‐mediated trend of lipids in the young mice but had little effect on the aged. (Study registration: Project DOI: https://doi.org/10.21228/M8K11W). Conclusions Pharmacological activation of Beclin‐1 by TB‐peptide is cardiac protective in both young and aged population during endotoxemia, suggest a therapeutic potential for sepsis‐induced cardiomyopathy. Metabolomics analysis suggests that an age‐independent protection by TB‐peptide is associated with reprograming of energy production via glucose and amino acid metabolisms.

Abstract Background Fecundity is important for farm blue fox (Vulpes lagopus), who with asthenospermia have be a problem in some of farms in China. A key symptom of asthenospermia is decreased sperm motility. The decreased secreting beta-defensin108 (vBD108) of blue fox is speculated be related to asthenospermia. To clarify this idea, the mRNA expression of vBD108 in testis and epididymis of blue foxes with asthenospermia were detected and compared to the healthy one. The antibody was prepared and analyzed by immunohistochemistry. Results The vBD108 in testis and epididymis was found both in blue fox with asthenospermia and healthy group by the method of immunohistochemistry. The expression of vBD108 mRNA in testes (P

In this work, we provide a computational method of regularized logistic regression for discovering biomarkers of spontaneous preterm birth (SPTB) from gene expression data. The successful identification of SPTB biomarkers will greatly benefit the interference of infant gestational age for reducing the risks of pregnant women and preemies. In recent years, various approaches have been proposed for the feature selection of identifying the subset of meaningful genes that can achieve accurate classification for disease samples from controls. Here, we comprehensively summarize the regularized logistic regression with seven effective penalties developed for the selection of strongly indicative genes of SPTB from microarray data. We compare their properties and assess their classification performances in multiple datasets. It shows that elastic net, lasso, L1/2 and SCAD penalties get the better performance than others and can be successfully used to identify biomarkers of SPTB. Particularly, we make a functional enrichment analysis on these biomarkers and construct a logistic regression classifier based on them. The classifier generates an indicator of preterm risk score (PRS) for predicting SPTB. Based on the trained predictor, we verify the identified biomarkers on an independent dataset. The biomarkers achieve the AUC value of 0.933 in the SPTB classification. The results demonstrate the effectiveness and efficiency of the built-up strategy of biomarker discovery with regularized logistic regression. Obviously, the proposed method of discovering biomarkers for SPTB can be easily extended for other complex diseases.

lncRNA–protein interactions play essential roles in a variety of cellular processes. However, the experimental methods for systematically mapping of lncRNA–protein interactions remain time-consuming and expensive. Therefore, it is urgent to develop reliable computational methods for predicting lncRNA–protein interactions. In this study, we propose a computational method called LncPNet to predict potential lncRNA–protein interactions by embedding an lncRNA–protein heterogenous network. The experimental results indicate that LncPNet achieves promising performance on benchmark datasets extracted from the NPInter database with an accuracy of 0.930 and area under ROC curve (AUC) of 0.971. In addition, we further compare our method with other eight state-of-the-art methods, and the results illustrate that our method achieves superior prediction performance. LncPNet provides an effective method via a new perspective of representing lncRNA–protein heterogenous network, which will greatly benefit the prediction of lncRNA–protein interactions.

Background: Mitochondrial deficiency is a known pathology in sepsis-induced organ failure. We previously found that mitochondria-associated membranes (MAMs), a subcellular domain supporting mitochondrial status, are impaired in the heart during endotoxemia, suggesting a mechanism of mitochondrial damage occurred in sepsis. Mitophagy pathway via E3 ubiquitin ligase Parkin and PTEN-induced kinase 1 (PINK1) controls mitochondrial quality. Studies described here examined the impact of Parkin on cardiac MAMs and endotoxemia-induced cardiomyopathy. Additionally, point mutation W403A in Parkin was previously identified as a constitutively active mutation in vitro. In vivo effects of forced expression of this mutation were evaluated in the endotoxemia model.Methods: Mice of wild type (WT), Parkin-deficiency (Park2−/−), and knock-in expression of Parkin W402A (human Parkin W403A) were given lipopolysaccharide (LPS) challenge. Cardiac function was evaluated by echocardiography. In the harvested heart tissue, MAM fractions were isolated by ultracentrifugation, and their amount and function were quantified. Ultrastructure of MAMs and mitochondria was examined by electron microscopy. Mitochondrial respiratory activities were measured by enzyme assays. Myocardial inflammation was estimated by levels of pro-inflammatory cytokine IL-6. Myocardial mitophagy was assessed by levels of mitophagy factors associated with mitochondria and degrees of mitochondria-lysosome co-localization. Parkin activation, signified by phosphorylation on serine 65 of Parkin, was also evaluated.Results: Compared with WT, Park2−/− mice showed more severely impaired cardiac MAMs during endotoxemia, characterized by disrupted structure, reduced quantity, and weakened transporting function. Endotoxemia-induced cardiomyopathy was intensified in Park2−/− mice, shown by worsened cardiac contractility and higher production of IL-6. Mitochondria from the Park2−/− hearts were more deteriorated, indicated by losses in both structural integrity and respiration function. Unexpectedly, mice carrying Parkin W402A showed similar levels of cardiomyopathy and mitochondrial damage when compared with their WT counterparts. Further, Parkin W402A mutation neither enhanced mitophagy nor increased Parkin activation in myocardium under the challenge of endotoxemia.Conclusion: our results suggest that Parkin/PINK1 mitophagy participates in the regulation of cardiac MAMs during endotoxemia. Point mutation W402A (human W403A) in Parkin is not sufficient to alleviate cardiomyopathy induced by endotoxemia in vivo.

The androgen receptor (AR) plays a central role in prostate cancer. Development of castration resistant prostate cancer (CRPC) requires androgen-independent activation of AR, which involves its large N-terminal domain (NTD) and entails extensive epigenetic changes depending in part on histone lysine demethylases (KDMs) that interact with AR. The AR-NTD is rich in low-complexity sequences, including a polyQ repeat. Longer polyQ sequences were reported to decrease transcriptional activity and to protect against prostate cancer, although they can lead to muscular atrophy. However, the molecular mechanisms underlying these observations are unclear. Using NMR spectroscopy, here we identify weak interactions between the AR-NTD and the KDM4A catalytic domain, and between the AR ligand-binding domain and a central KDM4A region that also contains low-complexity sequences. We also show that the AR-NTD can undergo liquid-liquid phase separation in vitro, with longer polyQ sequences phase separating more readily. Moreover, longer polyQ sequences hinder nuclear localization in the absence of hormone and increase the propensity for formation of AR-containing puncta in the nucleus of cells treated with dihydrotestosterone. These results lead us to hypothesize that polyQ-dependent liquid-liquid phase separation may provide a mechanism to decrease the transcriptional activity of AR, potentially opening new opportunities to design effective therapies against CRPC and muscular atrophy.

AIM: To investigate the secretion levels of serum tumor markers neuron specific enolase(NSE), glycoprotein antigen CA153 and CA199 in patients with retinoblastoma(RB).METHODS: Data of 42 RB patients who received chemotherapy at Shenzhen People's Hospital between October 2017 and October 2019 were collected. The RB group was regrouped according to early and mid-term/advanced stages, and single/double eye involvement, the levels of serum tumor markers were compared between these subgroups and different genders.RESULTS: The levels of tumor markers NSE, CA153 and CA199 in the advanced group were higher than those in the early and mid-term groups(49.69±18.45ng/mL vs 36.18±14.92ng/mL, 22.38±12.03U/mL vs 15.10±8.32U/mL, 46.44±18.76U/mL vs 30.21±24.03U/mL, P0.05).CONCLUSION: Serum NSE, CA153 and CA199 are significantly increased in advanced group than early and mid-term group, and they are of great significance in the diagnosis and treatment of RB. The difference between NSE in single and double eyes is worthy of further study.

Six piperazine derivatives 6a–f containing 1,4-benzodioxan moiety have been synthesized and characterized by 1H NMR, ESI-MS and elemental analysis. The structure of 6d was further confirmed by single crystal X-ray diffraction. All these novel compounds were screened for their in vivo anti-inflammatory activity employing classical para-xylene-induced mice ear-swelling model. The results revealed that most of the target compounds showed significant anti-inflammatory activities, especially compound 6a with ortho-substituted methoxy group on the phenylpiperazine ring exhibited the best activity among the designed compounds.

Researchers have shown that adult zebrafish have the potential to regenerate 20% of the ventricular muscle within two months of apex resection, and neonatal mice have the capacity to regenerate their heart after apex resection up until day 7 after birth. The goal of this study was to determine if large mammals (porcine heart model) have the capability to fully regenerate a resected portion of the left ventricular apex during the neonatal stage, and if so, how long the regenerative potential persists. A total of 36 piglets were divided into the following groups: 0-day control and surgical groups and seven-day control and surgical groups. For the apex removal groups, each piglet was subjected to a partial wall thickness resection (~30% of the ventricular wall thickness). Heart muscle function was assessed via transthoracic echocardiograms; the seven-day surgery group experienced a decrease in ejection fraction and fractional shortening. Upon gross necropsy, for piglets euthanized four weeks post-surgery, all 0-day-old hearts showed no signs of scarring or any indication of the induced injury. Histological analysis confirmed that piglets in the 0-day surgery group exhibited various degrees of regeneration, with half of the piglets showing full regeneration and the other half showing partial regeneration. However, each piglet in the seven-day surgery group demonstrated epicardial fibrosis along with moderate to severe dissecting interstitial fibrosis, which was accompanied by an abundant collagenous extracellular matrix as the result of a scar formation in the resection site. Histology of one 0-day apex resection piglet (briefly lain on and accidentally killed by the mother sow three days post-surgery) revealed dense, proliferative mesenchymal cells bordering the fibrin and hemorrhage zone and differentiating toward immature cardiomyocytes. We further examined the heart explants at 5-days post-surgery (5D PO) and 1-week post-surgery (1W PO) to assess the repair progression. For the 0-day surgery piglets euthanized at 5D PO and 1W PO, half had abundant proliferating mesenchymal cells, suggesting active regeneration, while the other half showed increased extracellular collagen. The seven-day surgery piglets euthanized at 5D PO, and 1W PO showed evidence of greatly increased extracellular collagen, while some piglets had proliferating mesenchymal cells, suggesting a regenerative effort is ongoing while scar formation seems to predominate. In short, our qualitative findings suggest that the piglets lose the full myocardial regenerative potential by 7 days after birth, but greatly preserve the regenerative potential within 1 day post-partum.

Single-cell RNA-sequencing (scRNA-seq) is a recent high-throughput technique that can measure gene expression, reveal cell heterogeneity, rare and complex cell populations, and discover cell types and their relationships. The analysis of scRNA-seq data is challenging because of transcripts sparsity, replication noise, and outlier cell populations. A gene coexpression network (GCN) analysis effectively deciphers phenotypic differences in specific states by describing gene–gene pairwise relationships. The underlying gene modules with different coexpression patterns partially bridge the gap between genotype and phenotype. This study presents a new framework called scGENA (single-cell gene coexpression network analysis) for GCN analysis based on scRNA-seq data. Although there are several methods for scRNA-seq data analysis, we aim to build an integrative pipeline for several purposes that cover primary data preprocessing, including data exploration, quality control, normalization, imputation, and dimensionality reduction of clustering as downstream of GCN analysis. To demonstrate this integrated workflow, an scRNA-seq dataset of the human diabetic pancreas with 1600 cells and 39,851 genes was implemented to perform all these processes in practice. As a result, scGENA is demonstrated to uncover interesting gene modules behind complex diseases, which reveal biological mechanisms. scGENA provides a state-of-the-art method for gene coexpression analysis for scRNA-seq data.

Harmful algal blooms are a major type of marine disaster that endangers the marine ecological environment and human health. Since the algal bloom is a complex nonlinear process with time characteristics, traditional statistical methods often cannot provide good predictions. In this paper, we propose a method based on machine learning with the aim to predict the occurrence of algal blooms by environmental parameters. The features related to algal bloom growth have been experimented for achieving a good prediction of algal concentrations by a combination strategy. We validate the prediction performance on two real datasets from two locations in US and China, i.e., Scripps Pier, California and Sishili Bay, Shandong, respectively. The models and feature subsets have been selected to complete the missing data and predict the phytoplankton concentration. The results demonstrate the efficiency of our method in the short-term prediction of concentrations by selecting appropriate features. The comparison studies prove the advantage of our developed machine learning method. The importance of every features for the prediction performance reveals crucial factors for the outbreak of harmful algal blooms.

Education is essential to cultivate the innovations. The Nobel prize of physiology or medicine (medicine in short) has been awarded to the breakthrough achievement in medicine. It is of paramount interest to investigate the educational backgrounds of Nobel prize laureates in medicine for discovering their education characteristics. To this end, we compile an integrative and comprehensive dataset about the educational experiences of all the Nobel awardees focusing on their undergraduate and doctoral majors. Then we implement the summary statistics of numerous majors in current university disciplinary system and calculate the proportion of each major accordingly. The results demonstrate that 49.5% of the laureates of Nobel prize in medicine have multidisciplinary educational background and 69% of their undergraduate majors are not directly in medicine. After removing the laureates who graduated from the countries that do not provide undergraduate medical education, the percentage is still as high as 52%. To further check our findings, we compare the similarity and difference between the laureates in medicine and those in physics and chemistry. The percentages in the other two sciences are significantly different from those in medicine, which illustrates a distinct degree of specificity. The comparison confirms the uniqueness of academic backgrounds in major medical innovations, such as those awarded by Nobel prizes. The data-driven analysis highlights the specific pattern and inspires the urgent necessity of introducing science education in medicine for cultivating the original and cutting-edge medical innovations, especially in the undergraduate program.