Divergent composition and transposon-silencing activity of small RNAs in mammalian oocytes.

BACKGROUND: Small RNAs are essential for germ cell development and fertilization. However, fundamental questions remain, such as the level of conservation in small RNA composition between species and whether small RNAs control transposable elements in mammalian oocytes. RESULTS: Here, we use high-throughput sequencing to profile small RNAs and poly(A)-bearing long RNAs in oocytes of 12 representative vertebrate species (including 11 mammals). The results show that miRNAs are generally expressed in the oocytes of each representative species (although at low levels), whereas endo-siRNAs are specific to mice. Notably, piRNAs are predominant in oocytes of all species (except mice) and vary widely in length. We find PIWIL3-associated piRNAs are widespread in mammals and generally lack 3'-2'-O-methylation. Additionally, sequence identity is low between homologous piRNAs in different species, even among those present in syntenic piRNA clusters. Despite the species-specific divergence, piRNAs retain the capacity to silence younger TE subfamilies in oocytes. CONCLUSIONS: Collectively, our findings illustrate a high level of diversity in the small RNA populations of mammalian oocytes. Furthermore, we identify sequence features related to conserved roles of small RNAs in silencing TEs, providing a large-scale reference for future in-depth study of small RNA functions in oocytes.

Characterizing the cellular and molecular variabilities of peripheral immune cells in healthy recipients of BBIBP-CorV inactivated SARS-CoV-2 vaccine by single-cell RNA sequencing.

Over 3 billion doses of inactivated vaccines for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been administered globally. However, our understanding of the immune cell functional transcription and T cell receptor (TCR)/B cell receptor (BCR) repertoire dynamics following inactivated SARS-CoV-2 vaccination remains poorly understood. Here, we performed single-cell RNA and TCR/BCR sequencing on peripheral blood mononuclear cells at four time points after immunization with the inactivated SARS-CoV-2 vaccine BBIBP-CorV. Our analysis revealed an enrichment of monocytes, central memory CD4+ T cells, type 2 helper T cells and memory B cells following vaccination. Single-cell TCR-seq and RNA-seq comminating analysis identified a clonal expansion of CD4+ T cells (but not CD8+ T cells) following a booster vaccination that corresponded to a decrease in the TCR diversity of central memory CD4+ T cells and type 2 helper T cells. Importantly, these TCR repertoire changes and CD4+ T cell differentiation were correlated with the biased VJ gene usage of BCR and the antibody-producing function of B cells post-vaccination. Finally, we compared the functional transcription and repertoire dynamics in immune cells elicited by vaccination and SARS-CoV-2 infection to explore the immune responses under different stimuli. Our data provide novel molecular and cellular evidence for the CD4+ T cell-dependent antibody response induced by inactivated vaccine BBIBP-CorV. This information is urgently needed to develop new prevention and control strategies for SARS-CoV-2 infection. (ClinicalTrials.gov Identifier: NCT04871932).

Moderate static magnetic fields enhance antitumor CD8+ T cell function by promoting mitochondrial respiration.

With the discovery of magnetoreceptor mechanisms in animals, it materialized the novel applications of controlling cell and animal behaviors using magnetic fields. T cells have shown to be sensitive to magnetic fields. Here, we reported that exposure to moderate SMFs (static magnetic fields) led to increased granule and cytokine secretion as well as ATP production and mitochondrial respiration from CD8+ T cells. These effects were inhibited by knocking down the Uqcrb and Ndufs6 genes of mitochondrial respiratory chain, whose transcriptions were regulated by candidate magnetoreceptor genes Isca1 and Cry1/Cry2. SMF exposure also promoted CD8+ T cell granule and cytokine secretion and repressed tumor growth in vivo. SMFs enhanced CD8+ T cell cytotoxicity, and the adoptive transfer into tumor-bearing mice resulted in enhanced antitumor effects. Collectively, our study suggests that moderate SMFs enhance CD8+ T cell cytotoxicity by promoting mitochondrial respiration and promoted the antitumor function of CD8+ T cells.

Transcriptomic profile of early zebrafish PGCs by single cell sequencing.

Single cell RNA-seq is a powerful and sensitive way to capture the genome-wide gene expression. Here, single cell RNA-seq was utilized to study the transcriptomic profile of early zebrafish PGCs (primordial germ cells) at three different developmental stages. The three stages were 6, 11 and 24 hpf (hours post fertilization). For each developmental stage, three zebrafish PGCs from one embryo were collected, and 9 samples in total were used in this experiment. Single cell RNA-seq results showed that 5099-7376 genes were detected among the 9 samples, and the number of expressed genes decreased as development progressed. Based on the gene expression pattern, samples from 6 and 11 hpf clustered closely, while samples from 24 hpf were more dispersed. By WGCNA (weighted gene co-expression network analysis), the two biggest modules that had inverse gene expression patterns were found to be related to PGC formation or migration. Functional enrichment analysis for these two modules showed that PGCs mainly conducted migration and cell division in early development (6/11 hpf) and translation activity became active in late development (24 hpf). Differentially expressed gene analyses showed that more genes were downregulated than upregulated between two adjacent stages, and genes related to PGC formation or migration reported by previous studies decreased significantly from 11 to 24 hpf. Our results provide base knowledge about zebrafish PGC development at the single cell level and can be further studied by other researchers interested in biological development.

Single-cell CAS-seq reveals a class of short PIWI-interacting RNAs in human oocytes.

Small RNAs have important functions. However, small RNAs in primate oocytes remain unexplored. Herein, we develop CAS-seq, a single-cell small RNA sequencing method, and profile the small RNAs in human oocytes and embryos. We discover a class of ~20-nt small RNAs that are predominantly expressed in human and monkey oocytes, but not in mouse oocytes. They are specifically associated with HIWI3 (PIWIL3), whereas significantly shorter than the commonly known PIWI-interacting RNAs (piRNAs), designated as oocyte short piRNAs (os-piRNAs). Notably, the os-piRNAs in human oocytes lack 2'-O-methylation at the 3' end, a hallmark of the classic piRNAs. In addition, the os-piRNAs have a strong 1U/10 A bias and are enriched on the antisense strands of recently evolved transposable elements (TEs), indicating the potential function of silencing TEs by cleavage. Therefore, our study has identified an oocyte-specific piRNA family with distinct features and provides valuable resources for studying small RNAs in primate oocytes.

High throughput small RNA and transcriptome sequencing reveal capacitation-related microRNAs and mRNA in boar sperm.

BACKGROUND: Capacitation, a prerequisite for oocyte fertilization, is a complex process involving series of structural and functional changes in sperms such as membrane modifications, modulation of enzyme activities, and protein phosphorylation. In order to penetrate and fertilize an oocyte, mammalian sperms must undergo capacitation. Nevertheless, the process of sperm capacitation remains poorly understood and requires further elucidation. In the current study, via high throughput sequencing, we identified and explored the differentially expressed microRNAs (miRNAs) and mRNAs involved in boar sperm capacitation. RESULTS: We identified a total of 5342 mRNAs and 204 miRNAs that were differentially expressed in fresh and capacitated boar sperms. From these, 12 miRNAs (8 known and 4 newly identified miRNAs) and their differentially expressed target mRNAs were found to be involved in sperm capacitation-related PI3K-Akt, MAPK, cAMP-PKA and Ca2+signaling pathways. CONCLUSIONS: Our study is first to provide the complete miRNA and transcriptome profiles of boar sperm. Our findings provide important insights for the understanding of the RNA profile in boar sperm and future elucidation of the underlying molecular mechanism relevant to mammalian sperm capacitation.

Regulation of the terminal maturation of iNKT cells by mediator complex subunit 23.

Invariant natural killer T cells (iNKT cells) are a specific subset of T cells that recognize glycolipid antigens and upon activation rapidly exert effector functions. This unique function is established during iNKT cell development; the detailed mechanisms of this process, however, remain to be elucidated. Here the authors show that deletion of the mediator subunit Med23 in CD4+CD8+ double positive (DP) thymocytes completely blocks iNKT cell development at stage 2. This dysregulation is accompanied by a bias in the expression of genes related to the regulation of transcription and metabolism, and functional impairment of the cells including the loss of NK cell characteristics, reduced ability to secrete cytokines and attenuated recruitment capacity upon activation. Moreover, Med23-deficient iNKT cells exhibit impaired anti-tumor activity. Our study identifies Med23 as an essential transcriptional regulator that controls iNKT cell differentiation and terminal maturation.

Structural insights into the sequence-specific recognition of Piwi by Drosophila Papi.

The Tudor domain-containing (Tdrd) family proteins play a critical role in transposon silencing in animal gonads by recognizing the symmetrically dimethylated arginine (sDMA) on the (G/A)R motif of the N-terminal of PIWI family proteins via the eTud domains. Papi, also known as "Tdrd2," is involved in Zucchini-mediated PIWI-interacting RNA (piRNA) 3'-end maturation. Intriguingly, a recent study showed that, in papi mutant flies, only Piwi-bound piRNAs increased in length, and not Ago3-bound or Aub-bound piRNAs. However, the molecular and structural basis of the Papi-Piwi complex is still not fully understood, which limits mechanistic understanding of the function of Papi in piRNA biogenesis. In the present study, we determined the crystal structures of Papi-eTud in the apo form and in complex with a peptide containing unmethylated or dimethylated R10 residues. Structural and biochemical analysis showed that the Papi interaction region on the Drosophila Piwi contains an RGRRR motif (R7-R11) distinct from the consensus (G/A)R motif recognized by canonical eTud. Mass spectrometry results indicated that Piwi is the major binding partner of Papi in vivo. The papi mutant flies suffered from both fertility and transposon-silencing defects, supporting the important role conferred to Papi in piRNA 3' processing through direct interaction with Piwi proteins.

Highly sensitive sequencing reveals dynamic modifications and activities of small RNAs in mouse oocytes and early embryos.

Small RNAs play important roles in early embryonic development. However, their expression dynamics and modifications are poorly understood because of the scarcity of RNA that is obtainable for sequencing analysis. Using an improved deep sequencing method that requires as little as 10 ng of total RNA or 50 oocytes, we profile small RNAs in mouse oocytes and early embryos. We find that microRNA (miRNA) expression starts soon after fertilization, and the mature miRNAs carried into the zygote by sperm during fertilization are relatively rare compared to the oocyte miRNAs. Intriguingly, the zygotic miRNAs display a marked increase in 3' mono- and oligoadenylation in one- to two-cell embryos, which may protect the miRNAs from the massive degradation taking place during that time. Moreover, bioinformatics analyses show that the function of miRNA is suppressed from the oocyte to the two-cell stage and appears to be reactivated after the two-cell stage to regulate genes important in embryonic development. Our study thus provides a highly sensitive profiling method and valuable data sets for further examination of small RNAs in early embryos.

MyCompoundID MS/MS Search: Metabolite Identification Using a Library of Predicted Fragment-Ion-Spectra of 383,830 Possible Human Metabolites.

We report an analytical tool to facilitate metabolite identification based on an MS/MS spectral match of an unknown to a library of predicted MS/MS spectra of possible human metabolites. To construct the spectral library, the known endogenous human metabolites in the Human Metabolome Database (HMDB) (8,021 metabolites) and their predicted metabolic products via one metabolic reaction in the Evidence-based Metabolome Library (EML) (375,809 predicted metabolites) were subjected to in silico fragmentation to produce the predicted MS/MS spectra. This spectral library is hosted at the public MCID Web site ( www.MyCompoundID.org ), and a spectral search program, MCID MS/MS, has been developed to allow a user to search one or a batch of experimental MS/MS spectra against the library spectra for possible match(s). Using MS/MS spectra generated from standard metabolites and a human urine sample, we demonstrate that this tool is very useful for putative metabolite identification. It allows a user to narrow down many possible structures initially found by using an accurate mass search of an unknown metabolite to only one or a few candidates, thereby saving time and effort in selecting or synthesizing metabolite standard(s) for eventual positive metabolite identification.

Identification of lncRNA-associated competing triplets reveals global patterns and prognostic markers for cancer.

Recent studies have suggested that long non-coding RNAs (lncRNAs) can interact with microRNAs (miRNAs) and indirectly regulate miRNA targets though competing interactions. However, the molecular mechanisms underlying these interactions are still largely unknown. In this study, these lncRNA-miRNA-gene interactions were defined as lncRNA-associated competing triplets (LncACTs), and an integrated pipeline was developed to identify lncACTs that are active in cancer. Competing lncRNAs had sponge features distinct from non-competing lncRNAs. In the lncACT cross-talk network, disease-associated lncRNAs, miRNAs and coding-genes showed specific topological patterns indicative of their competence and control of communication within the network. The construction of global competing activity profiles revealed that lncACTs had high activity specific to cancers. Analyses of clustered lncACTs revealed that they were enriched in various cancer-related biological processes. Based on the global cross-talk network and cluster analyses, nine cancer-specific sub-networks were constructed. H19- and BRCA1/2-associated lncACTs were able to discriminate between two groups of patients with different clinical outcomes. Disease-associated lncACTs also showed variable competing patterns across normal and cancer patient samples. In summary, this study uncovered and systematically characterized global properties of human lncACTs that may have prognostic value for predicting clinical outcome in cancer patients.

Identifying a polymorphic 'switch' that influences miRNAs' regulation of a myasthenia gravis risk pathway.

The significant roles of genetic variants in myasthenia gravis (MG) pathogenesis have been demonstrated in many studies, and recently it has been revealed that aberrant level/regulation of microRNAs (miRNAs) might contribute to the initiation and progression of MG. However, the dysfunction of miRNA associated with single nucleotide polymorphisms (miRSNPs) has not been well investigated in MG. In this study, we created a contemporary catalog of 89 MG risk genes via manual literature-mining. Based on this risk gene catalog, we obtained 18 MG risk pathways. Furthermore, we identified 93 miRNAs that target MG risk pathways and revealed the miRSNPs 'switches' in miRNA regulation in the MG risk pathways by integrating the database information of miRSNPs. We also constructed a miRNA-mediated SNP switching pathway network (MSSPN) to intuitively analyze miRNA regulation of MG risk pathways and the relationship of the polymorphism 'switch' with these changes in regulation. Moreover, we carried out in-depth dissection on the correlation between hsa05200 (pathway in cancer) and MG development, and elaborated the significance of 4 high-risk genes. By network analysis and literature mining, we proposed a potential mechanism of miRSNPs→gene→pathway effects on MG pathogenesis, especially for rs28457673 (miR-15/16/195/424/497 family)→IGF1R→hsa05200 (pathway in cancer). Therefore, our studies have revealed a functional role for genetic modulators in MG pathogenesis at a systemic level, which could be informative for further miRNA and miRSNPs studies in MG.

SNP@lincTFBS: an integrated database of polymorphisms in human LincRNA transcription factor binding sites.

Large intergenic non-coding RNAs (lincRNAs) are a new class of functional transcripts, and aberrant expression of lincRNAs was associated with several human diseases. The genetic variants in lincRNA transcription factor binding sites (TFBSs) can change lincRNA expression, thereby affecting the susceptibility to human diseases. To identify and annotate these functional candidates, we have developed a database SNP@lincTFBS, which is devoted to the exploration and annotation of single nucleotide polymorphisms (SNPs) in potential TFBSs of human lincRNAs. We identified 6,665 SNPs in 6,614 conserved TFBSs of 2,423 human lincRNAs. In addition, with ChIPSeq dataset, we identified 139,576 SNPs in 304,517 transcription factor peaks of 4,813 lincRNAs. We also performed comprehensive annotation for these SNPs using 1000 Genomes Project datasets across 11 populations. Moreover, one of the distinctive features of SNP@lincTFBS is the collection of disease-associated SNPs in the lincRNA TFBSs and SNPs in the TFBSs of disease-associated lincRNAs. The web interface enables both flexible data searches and downloads. Quick search can be query of lincRNA name, SNP identifier, or transcription factor name. SNP@lincTFBS provides significant advances in identification of disease-associated lincRNA variants and improved convenience to interpret the discrepant expression of lincRNAs. The SNP@lincTFBS database is available at http://bioinfo.hrbmu.edu.cn/SNP_lincTFBS.

Identification of a core miRNA-pathway regulatory network in glioma by therapeutically targeting miR-181d, miR-21, miR-23b, ß-Catenin, CBP, and STAT3.

The application of microRNAs (miRNAs) in the therapeutics of glioma and other human diseases is an area of intense interest. However, it's still a great challenge to interpret the functional consequences of using miRNAs in glioma therapy. Here, we examined paired deep sequencing expression profiles of miRNAs and mRNAs from human glioma cell lines after manipulating the levels of miRNAs miR-181d, -21, and -23b, as well as transcriptional regulators ß-catenin, CBP, and STAT3. An integrated approach was used to identify functional miRNA-pathway regulatory networks (MPRNs) responding to each manipulation. MiRNAs were identified to regulate glioma related biological pathways collaboratively after manipulating the level of either post-transcriptional or transcriptional regulators, and functional synergy and crosstalk was observed between different MPRNs. MPRNs responsive to multiple interventions were found to occupy central positions in the comprehensive MPRN (cMPRN) generated by integrating all the six MPRNs. Finally, we identified a core module comprising 14 miRNAs and five pathways that could predict the survival of glioma patients and represent potential targets for glioma therapy. Our results provided novel insight into miRNA regulatory mechanisms implicated in therapeutic interventions and could offer more inspiration to miRNA-based glioma therapy.

LincSNP: a database of linking disease-associated SNPs to human large intergenic non-coding RNAs.

BACKGROUND: Genome-wide association studies (GWAS) have successfully identified a large number of single nucleotide polymorphisms (SNPs) that are associated with a wide range of human diseases. However, many of these disease-associated SNPs are located in non-coding regions and have remained largely unexplained. Recent findings indicate that disease-associated SNPs in human large intergenic non-coding RNA (lincRNA) may lead to susceptibility to diseases through their effects on lincRNA expression. There is, therefore, a need to specifically record these SNPs and annotate them as potential candidates for disease. DESCRIPTION: We have built LincSNP, an integrated database, to identify and annotate disease-associated SNPs in human lincRNAs. The current release of LincSNP contains approximately 140,000 disease-associated SNPs (or linkage disequilibrium SNPs), which can be mapped to around 5,000 human lincRNAs, together with their comprehensive functional annotations. The database also contains annotated, experimentally supported SNP-lincRNA-disease associations and disease-associated lincRNAs. It provides flexible search options for data extraction and searches can be performed by disease/phenotype name, SNP ID, lincRNA name and chromosome region. In addition, we provide users with a link to download all the data from LincSNP and have developed a web interface for the submission of novel identified SNP-lincRNA-disease associations. CONCLUSIONS: The LincSNP database aims to integrate disease-associated SNPs and human lincRNAs, which will be an important resource for the investigation of the functions and mechanisms of lincRNAs in human disease. The database is available at http://bioinfo.hrbmu.edu.cn/LincSNP.

PEP search in MyCompoundID: detection and identification of dipeptides and tripeptides using dimethyl labeling and hydrophilic interaction liquid chromatography tandem mass spectrometry.

Small peptides, such as dipeptides and tripeptides, are naturally present in many biological samples (e.g., human biofluids and cell extracts). They have attracted great attention in many research fields because of their important biological functions as well as potential roles as disease biomarkers. Tandem mass spectrometry (MS/MS) can be used to profile these small peptides. However, the type and number of fragment ions generated in MS/MS are often limited for unambiguous identification. Herein we report a novel database-search strategy based on the use of MS/MS spectra of both unlabeled and dimethyl labeled peptides to identify and confirm amino acid sequences of di/tripeptides that are separated using hydrophilic interaction (HILIC) liquid chromatography (LC). To facilitate the di/tripeptide identification, a database consisting of all the predicted MS/MS spectra from 400 dipeptides and 8000 tripeptides was created, and a search tool, PEP Search, was developed and housed at the MyCompoundID website ( www.mycompoundid.org/PEP). To evaluate the identification specificity of this method, we used acid hydrolysis to degrade a standard protein, cytochrome c, to produce many di/tripeptides with known sequences for LC/MS/MS. The resultant MS/MS spectra were searched against the database to generate a list of matches which were compared to the known sequences. We correctly identified the di/tripeptides in the protein hydrolysate. We then applied this method to detect and identify di/tripeptides naturally present in human urine samples with high confidence. We envisage the use of this method as a complementary tool to various LC/MS techniques currently available for small molecule or metabolome profiling with an added benefit of covering all di/tripeptide chemical space.

A global map for dissecting phenotypic variants in human lincRNAs.

Large intergenic noncoding RNAs (lincRNAs) are emerging as key factors of multiple cellular processes. Cumulative evidence has linked lincRNA polymorphisms to diverse diseases. However, the global properties of lincRNA polymorphisms and their implications for human disease remain largely unknown. Here we performed a systematic analysis of naturally occurring variants in human lincRNAs, with a particular focus on lincRNA polymorphism as novel risk factor of disease etiology. We found that lincRNAs exhibited a relatively low level of polymorphisms, and low single-nucleotide polymorphism (SNP) density lincRNAs might have a broad range of functions. We also found that some polymorphisms in evolutionarily conserved regions of lincRNAs had significant effects on predicted RNA secondary structures, indicating their potential contribution to diseases. We mapped currently available phenotype-associated SNPs to lincRNAs and found that lincRNAs were associated with a wide range of human diseases. Some lincRNAs could be responsible for particular diseases. Our results provided not only a global perspective on genetic variants in human lincRNAs but also novel insights into the function and etiology of lincRNA. All the data in this study can be accessed and retrieved freely via a web server at http://bioinfo.hrbmu.edu.cn/lincPoly.

MyCompoundID: using an evidence-based metabolome library for metabolite identification.

Identification of unknown metabolites is a major challenge in metabolomics. Without the identities of the metabolites, the metabolome data generated from a biological sample cannot be readily linked with the proteomic and genomic information for studies in systems biology and medicine. We have developed a web-based metabolite identification tool ( http://www.mycompoundid.org ) that allows searching and interpreting mass spectrometry (MS) data against a newly constructed metabolome library composed of 8,021 known human endogenous metabolites and their predicted metabolic products (375,809 compounds from one metabolic reaction and 10,583,901 from two reactions). As an example, in the analysis of a simple extract of human urine or plasma and the whole human urine by liquid chromatography-mass spectrometry and MS/MS, we are able to identify at least two times more metabolites in these samples than by using a standard human metabolome library. In addition, it is shown that the evidence-based metabolome library (EML) provides a much superior performance in identifying putative metabolites from a human urine sample, compared to the use of the ChemPub and KEGG libraries.

mirTarPri: improved prioritization of microRNA targets through incorporation of functional genomics data.

MicroRNAs (miRNAs) are a class of small (19-25 nt) non-coding RNAs. This important class of gene regulator downregulates gene expression through sequence-specific binding to the 3'untranslated regions (3'UTRs) of target mRNAs. Several computational target prediction approaches have been developed for predicting miRNA targets. However, the predicted target lists often have high false positive rates. To construct a workable target list for subsequent experimental studies, we need novel approaches to properly rank the candidate targets from traditional methods. We performed a systematic analysis of experimentally validated miRNA targets using functional genomics data, and found significant functional associations between genes that were targeted by the same miRNA. Based on this finding, we developed a miRNA target prioritization method named mirTarPri to rank the predicted target lists from commonly used target prediction methods. Leave-one-out cross validation has proved to be successful in identifying known targets, achieving an AUC score up to 0. 84. Validation in high-throughput data proved that mirTarPri was an unbiased method. Applying mirTarPri to prioritize results of six commonly used target prediction methods allowed us to find more positive targets at the top of the prioritized candidate list. In comparison with other methods, mirTarPri had an outstanding performance in gold standard and CLIP data. mirTarPri was a valuable method to improve the efficacy of current miRNA target prediction methods. We have also developed a web-based server for implementing mirTarPri method, which is freely accessible at http://bioinfo.hrbmu.edu.cn/mirTarPri.

Prioritizing cancer therapeutic small molecules by integrating multiple OMICS datasets.

Drug design is crucial for the effective discovery of anti-cancer drugs. The success or failure of drug design often depends on the leading compounds screened in pre-clinical studies. Many efforts, such as in vivo animal experiments and in vitro drug screening, have improved this process, but these methods are usually expensive and laborious. In the post-genomics era, it is possible to seek leading compounds for large-scale candidate small-molecule screening with multiple OMICS datasets. In the present study, we developed a computational method of prioritizing small molecules as leading compounds by integrating transcriptomics and toxicogenomics data. This method provides priority lists for the selection of leading compounds, thereby reducing the time required for drug design. We found 11 known therapeutic small molecules for breast cancer in the top 100 candidates in our list, 2 of which were in the top 10. Furthermore, another 3 of the top 10 small molecules were recorded as closely related to cancer treatment in the DrugBank database. A comparison of the results of our approach with permutation tests and shared gene methods demonstrated that our OMICS data-based method is quite competitive. In addition, we applied our method to a prostate cancer dataset. The results of this analysis indicated that our method surpasses both the shared gene method and random selection. These analyses suggest that our method may be a valuable tool for directing experimental studies in cancer drug design, and we believe this time- and cost-effective computational strategy will be helpful in future studies in cancer therapy.