KLK7, KLK10, and KLK11 in Papillary Thyroid Cancer: Bioinformatic Analysis and Experimental Validation

Despite many studies on papillary thyroid carcinoma (PTC) in the past few decades, some critical and significant genes remain undiscovered. To explore genes that may play crucial roles in PTC, a detailed analysis of the expression levels, mutations, and clinical significance of Kallikrein-related peptidases (KLKs) family genes in PTC was undertaken to provide new targets for the precise treatment of the disease. A comprehensive analysis of KLK family genes was performed using various online tools, such as GEPIA, Kaplan-Meier Plotter, LinkedOmics, GSCA, TIMER, and Cluego. KLK7, KLK10, and KLK11 were critical factors of KLK family genes. Then, functional assays were carried out on KLK7/10/11 to determine their proliferation, migration, and invasion capabilities in PTC. The mRNA expression levels of KLK7, KLK10, KLK11, and KLK13 were significantly elevated in thyroid carcinoma, while KLK1, KLK2, KLK3 and KLK4 mRNA levels were decreased compared to normal tissues. Correlations between KLK2/7-12/15 expression levels and tumor stage were also observed in thyroid carcinoma. Survival analysis demonstrated that KLK4/5/7/9-12/14 was associated with overall survival in patients with thyroid cancer. Not only were KLK genes strongly associated with cancer-related pathways, but also KLK7/10/11 was associated with immune-cell infiltration. Finally, silencing KLK7/10/11 impaired human papillary thyroid carcinoma cells' growth, migration ability, and invasiveness. The increased expression of KLK7, KLK10, and KLK11 may serve as molecular markers to identify PTC patients. KLK7, KLK10, and KLK11 could be potential prognostic indicators and targets for precision therapy against PTC.

KLK6/PAR1 Axis Promotes Tumor Growth and Metastasis by Regulating Cross-Talk between Tumor Cells and Macrophages

Kallikrein-related peptidase (KLK)6 is associated with inflammatory diseases and neoplastic progression. KLK6 is aberrantly expressed in several solid tumors and regulates cancer development, metastatic progression, and drug resistance. However, the function of KLK6 in the tumor microenvironment remains unclear. This study aimed to determine the role of KLK6 in the tumor microenvironment. Here, we uncovered the mechanism underlying KLK6-mediated cross-talk between cancer cells and macrophages. Compared with wild-type mice, KLK6-/- mice showed less tumor growth and metastasis in the B16F10 melanoma and Lewis lung carcinoma (LLC) xenograft model. Mechanistically, KLK6 promoted the secretion of tumor necrosis factor-alpha (TNF-α) from macrophages via the activation of protease-activated receptor-1 (PAR1) in an autocrine manner. TNF-α secreted from macrophages induced the release of the C-X-C motif chemokine ligand 1 (CXCL1) from melanoma and lung carcinoma cells in a paracrine manner. The introduction of recombinant KLK6 protein in KLK6-/- mice rescued the production of TNF-α and CXCL1, tumor growth, and metastasis. Inhibition of PAR1 activity suppressed these malignant phenotypes rescued by rKLK6 in vitro and in vivo. Our findings suggest that KLK6 functions as an important molecular link between macrophages and cancer cells during malignant progression, thereby providing opportunities for therapeutic intervention.

Activity-Based Probes for Proteases Pave the Way to Theranostic Applications

Proteases are important enzymes in health and disease. Their activities are regulated at multiple levels. In fact, proteases are synthesized as inactive proenzymes (zymogens) that are activated by proteolytic removal of their pro-peptide sequence and can remain active or their activity can be attenuated by complex formation with specific endogenous inhibitors or by limited proteolysis or degradation. Consequently, quite often, only a fraction of the protease molecules is in the active/functional form, thus, the abundance of a protease is not always linearly proportional to the (patho)physiological function(s). Therefore, assays to determine the active forms of proteases are needed, not only in research but also in molecular diagnosis and therapy. Activity-based probes (ABPs) are chemical entities that bind covalently to the active enzyme/protease. ABPs carry a detection tag to enable localization and quantification of specific enzymatic/proteolytic activities with applications in molecular imaging and diagnosis. Moreover, ABPs act as suicide inhibitors of proteases, which can be exploited for delineation of the functional role(s) of a given protease in (patho) biological context and as potential therapeutics. In this sense, ABPs represent new theranostic agents. We outline recent developments pertaining to ABPs for proteases with potential therapeutic applications, with the aim to highlight their importance in theranostics.

Human Tissue Kallikreins-Related Peptidases Are Targets for the Treatment of Skin Desquamation Diseases

Human tissue Kallikrein-related peptidases (hKLKs) are serine proteases distributed in several tissues that are involved in several biological processes. In skin, many are responsible for skin desquamation in the Stratum Corneum (SC) of the epidermis, specially hKLK5, hKLK7, hKLK6, hKLK8, and hKLK14. In SC, hKLKs cleave proteins of corneodesmosomes, an important structure responsible to maintain corneocytes attached. As part of skin desquamation, hKLKs are also involved in skin diseases with abnormal desquamation and inflammation, such as Atopic Dermatitis (AD), psoriasis, and the rare disease Netherton Syndrome (NS). Many studies point to hKLK overexpression or overactive in skin diseases, and they are also part of the natural skin inflammation process, through the PAR2 cleavage pathway. Therefore, the control of hKLK activity may offer successful treatments for skin diseases, improving the quality of life in patients. Diseases like AD, Psoriasis, and NS have an impact on social life, causing pain, itchy and mental disorders. In this review, we address the molecular mechanisms of skin desquamation, emphasizing the roles of human tissue Kallikrein-related peptidases, and the promising therapies targeting the inhibition of hKLKs.

Reconstructing the epidermal proteolytic cascades in health and disease

The epidermis is the outer stratified epithelium of the skin, forming the physical barrier that is indispensable for homeostasis. Epidermal proteolysis, mainly but not exclusively executed by kallikrein-related peptidases (KLKs), is tightly regulated to ensure maintenance of physiological skin renewal and an intact skin barrier. Perturbation of epidermal proteolytic networks is implicated in a wide array of rare and common skin pathologies of diverse genetic backgrounds. Recent studies of monogenic human skin diseases and newly developed animal models have revealed new mechanisms of regulation of proteolytic pathways in epidermal physiology and in disease states. These new data have challenged some accepted views, for example the role of matriptase in epidermal desquamation, which turned out to be restricted to mouse skin. The significance of PAR2 signaling in skin inflammation should also be reconsidered in the face of recent findings. Cumulatively, recent studies necessitate a sophisticated redefinition of the proteolytic and signaling pathways that operate in human skin. We elaborate how epidermal proteolysis is finely regulated at multiple levels, and in a spatial manner that was not taken into consideration so far, in which specific proteases are confined to distinct epidermal sublayers. Of interest, transglutaminases have emerged as regulators of epidermal proteolysis and desquamation by spatially fixing endogenous protease inhibitors, constituting regulatory factors that were not recognized before. Furthermore, new evidence suggests a link between proteolysis and lipid metabolism. By synthesis of established notions and recent discoveries, we provide an up-to-date critical parathesis of current knowledge and the extended complexity of proteolysis regulation and signaling pathways in skin. This article is protected by copyright. All rights reserved.

Kallikreins emerge as new regulators of viral infections

Kallikrein-related peptidases (KLKs) or kallikreins have been linked to diverse (patho) physiological processes, such as the epidermal desquamation and inflammation, seminal clot liquefaction, neurodegeneration, and cancer. Recent mounting evidence suggests that KLKs also represent important regulators of viral infections. It is well-established that certain enveloped viruses, including influenza and coronaviruses, require proteolytic processing of their hemagglutinin or spike proteins, respectively, to infect host cells. Similarly, the capsid protein of the non-enveloped papillomavirus L1 should be proteolytically cleaved for viral uncoating. Consequently, extracellular or membrane-bound proteases of the host cells are instrumental for viral infections and represent potential targets for drug development. Here, we summarize how extracellular proteolysis mediated by the kallikreins is implicated in the process of influenza (and potentially coronavirus and papillomavirus) entry into host cells. Besides direct proteolytic activation of viruses, KLK5 and 12 promote viral entry indirectly through proteolytic cascade events, like the activation of thrombolytic enzymes that also can process hemagglutinin, while additional functions of KLKs in infection cannot be excluded. In the light of recent evidence, KLKs represent potential host targets for the development of new antivirals. Humanized animal models to validate their key functions in viral infections will be valuable.

Research in practice: Towards deciphering the role of epidermal proteases in recessive dystrophic epidermolysis bullosa progression

Recessive dystrophic epidermolysis bullosa (RDEB) is an incurable severe skin disease caused by loss of collagen VII, an extracellular protein that ensures skin cohesion. It manifests in skin blistering and unresolved cycles of wounding and healing that progressively lead to dermal stiffening and early development of aggressive cutaneous squamous cell carcinomas. Inflammation and subsequent tissue fibrosis highly contribute to RDEB pathogenicity and targeting them could provide new therapeutic options. Kallikreins (KLKs) are epidermal secreted proteases, which contribute to skin desquamation and inflammation. Kallikreins are involved in the pathogenesis of several inflammatory skin disorders, but interestingly also in the initiation and progression of different cancers. Our project aims at deciphering the role of KLKs in inflammation, fibrosis, and tumor development in RDEB.

Tristetraprolin expression by keratinocytes protects against skin carcinogenesis

Cancer is caused primarily by genomic alterations resulting in deregulation of gene regulatory circuits in key growth, apoptosis, or DNA repair pathways. Multiple genes associated with the initiation and development of tumors are also regulated at the level of mRNA decay, through the recruitment of RNA-binding proteins to AU-rich elements (AREs) located in their 3'-untranslated regions. One of these ARE-binding proteins, tristetraprolin (TTP; encoded by Zfp36), is consistently dysregulated in many human malignancies. Herein, using regulated overexpression or conditional ablation in the context of cutaneous chemical carcinogenesis, we show that TTP represents a critical regulator of skin tumorigenesis. We provide evidence that TTP controlled both tumor-associated inflammation and key oncogenic pathways in neoplastic epidermal cells. We identify Areg as a direct target of TTP in keratinocytes and show that EGFR signaling potentially contributed to exacerbated tumor formation. Finally, single-cell RNA-Seq analysis indicated that ZFP36 was downregulated in human malignant keratinocytes. We conclude that TTP expression by epidermal cells played a major role in the control of skin tumorigenesis.

Redirecting drug repositioning to discover innovative cosmeceuticals

Skin appearance is essential for self-esteem and quality of life; consequently, skin care products represent a huge market. In particular, cosmeceuticals constitute a hybrid category of skin care formulations, at the interphase of cosmetics and pharmaceuticals, rationally designed to target (patho) physiological mechanisms aiming to enhance skin health and appearance. Cosmeceuticals are marketed as anti-ageing, anti-wrinkle, hair regrowth, skin whitening and wound healing agents with special emphasis on scar-free healing. An overview on recent cutting-edge advances concerning the discovery and development of enhanced performance cosmeceuticals by drug repositioning approaches is presented here. In this context, we propose "target repositioning," a new term, to highlight that druggable protein targets implicated in multiple diseases (hubs in the diseasome) can be exploited to accelerate the discovery of molecularly targeted cosmeceuticals that can promote skin health as an added benefit, which is a novel concept not described before. In this direction, emphasis is placed on the role of mouse models, for often untreatable skin diseases, as well as recent breakthroughs on monogenic rare skin syndromes, in promoting compound repositioning to innovative cosmeceuticals.

Ectopic expression of KLK6 in MDA-MB-435 melanoma cells reduces tumorigenicity in vivo

Melanoma is an aggressive form of cancer with poor prognosis therefore, identification of associated pathophysiological mechanisms is imperative towards the development of new therapeutic strategies. The KLK6 is a serine protease normally expressed in the epidermis. Recently, we found that elimination of Klk6 in mice results in enhanced resistance to chemically induced non-melanoma skin cancer. To delineate putative roles of KLK6 in melanoma, the invasive KLK6-non-expressing MDA-MB-435 melanoma cell line was stably transfected with the full-length KLK6 cDNA and expression of the corresponding RNA and protein were confirmed. Interestingly, restoration of KLK6 expression resulted in markedly suppressed growth of primary tumors when orthotopically implanted in SCID mice. Analysis of data retrieved from the human protein atlas revealed that melanomas with high KLK6 expression have a trend for longer survival. Collectively, we suggest that KLK6 inhibits growth of melanomas.

Comprehensive Identification of the Human Secretome as Potential Indicators in Treatment Outcome of HPV-Positive and -Negative Cervical Cancer Patients

BACKGROUND

The aim of this work was to explore the novel and promising biomarkers for the diagnosis and prognosis of cervical cancer patients.

METHODS

The secretome of primary cervical tissues was extracted and then determined by the LC-MS/MS assay. The level of screened targets was confirmed using the RT-PCR and ELISA in cervical cancer tissue samples. The median expression level of certain targets was used as a cutoff value to divide the patients into 2 groups, and then the patients were followed up. The predictive abilities of the targets on the prognosis were further studied.

RESULTS

LC-MS/MS, together with bioinformatic analysis, demonstrated that totally 95 targets were dysregulated in cervical cancer. Among them, ECM2, KLK6, and MASP1 were increased in cervical cancer in a stage-dependent manner, whereas FGA was negatively associated with the stage of cervical cancers. Overall survival (OS) and disease-free survival (DFS) rates were significantly decreased in the KLK6 high group, whereas little difference was found between the high and low groups of other 3 cases. Univariate analysis of the 5-year OS and DFS revealed a significantly worse outcome for patients with KLK6 high tumors. In multivariate analysis, KLK6 remained a highly significant prognostic marker for OS and DFS. Combined survival analysis of KLK6 expression and the HPV infection revealed that KLK6highHPV(-) predicted the most poor OS rate and the KLK6lowHPV(+) group showed the best prognosis.

CONCLUSION

Through the secretome analysis, we identified a series of secreted proteins differentially expressed in the clinical cancer, among which KLK6 has the potential to become a promising biomarker for the diagnosis and prognosis of cervical cancer patients.

Deficiency of the serine peptidase Kallikrein 6 does not affect the levels and the pathological accumulation of a-synuclein in mouse brain

Several lines of evidence indicate that the propagation of misfolded α-synuclein (α-syn) plays a central role in the progression and manifestation of Parkinson's disease. Pathogenic α-syn species can be present in the extracellular space. Thus, the identification and modulation of the key enzymes implicated in extracellular α-syn turnover becomes vital. Kallikrein peptidase 6 has been identified as one of the major α-syn degrading enzymes and has been implicated in the clearance of extracellular α-syn. However, the physiological role of this enzyme in regulating α-syn, in vivo, still remains elusive. Here, by utilizing Klk6 knock-out (Klk6^-/- ) mice as our experimental model, we provide insight into the physiologic relevance of endogenous KLK6 expression on α-syn processing. Behavioral phenotyping showed that Klk6^-/- mice display no gross behavioral abnormalities. Further in vivo characterization of this mouse model, in the context of α-syn accumulation, showed that KLK6 deletion had no impact on the protein levels of intracellular or extracellular α-syn. Upon in vivo administration of α-syn pre-formed fibrils (PFF), α-syn pathologic accumulations were evident both in the brains of Klk6^-/- mice and wt mice without significant differences. Intrastriatal delivery of active KLK6, did not affect secreted α-syn levels observed in the A53T α-syn over-expressing mice. These findings suggest that in the in vivo setting of PFF pathology induction, KLK6 alone is not able to modulate pathology transmission. Our study raises implications for the use of recombinant α-syn fibrils in α-syn turnover studies.

The role of kallikreins in inflammatory skin disorders and their potential as therapeutic targets

The skin is a vital organ of the human body, serving numerous protective and functional roles that are essential for survival. Residing in the epidermis are various epidermal proteases responsible for the establishment and regulation of barrier function. The human tissue kallikrein-related peptidase family conserves homeostasis of the skin barrier through their roles in desquamation, antimicrobial defense, innate immune response, and barrier maintenance. The activity of kallikreins is tightly regulated and dysregulation of kallikrein activity is seen to contribute to the formation of several inflammatory skin disorders. This review highlights the roles of kallikreins in skin homeostasis and pathologies. Due to their part in these skin disorders, inhibitors of the skin kallikreins have become attractive therapeutics. Over the past few years, both natural and synthetic inhibitors of several kallikreins have been identified and are undergoing further development as treatments to restore compromised barrier function. This review summarizes the kallikrein inhibitors under development for this purpose. These inhibitors remain promising therapeutics in cases of severe skin inflammation not well managed by current therapies.

Cathelicidin represents a new target for manipulation of skin inflammation in Netherton syndrome

Netherton syndrome (NS) is a severe ichthyosis caused by inactivating mutations in the SPINK5 gene encoding the serine protease inhibitor LEKTI. Spink5^-/- mice recapitulate NS and die perinatally from extensive dehydration as a result of a severe defect of the epidermal barrier. We showed that deletion of Klk5 in Spink5^-/- rescues neonatal lethality (Furio et al., 2015). However, Spink5^-/-Klk5^-/- mice developed skin shedding and inflammation during the first week from birth and the majority (70%) succumbed on P7. The remaining mice lived short (i.e. mean survival was 5 months) indicating alternative inflammatory pathways. Since cathelicidin is increased in Spink5^-/- epidermis, we investigated whether it could be implicated in NS pathology. Ablation of Camp in Spink5^-/- suppressed epidermal inflammation and restored abnormal epidermal differentiation, nevertheless, it failed to inhibit overdesquamation and Spink5^-/-Camp^-/- succumbed perinatally due to skin barrier defect, similarly to Spink5^-/-. Joint invalidation of Klk5 and Camp significantly extended survival of Spink5^-/-Klk5^-/-Camp^-/- mice. We provide evidence that cathelicidin is implicated in NS-associated skin inflammation in vivo. Therefore, marketed products that are known to reduce cathelicidin expression could be repurposed for the management of NS.

Kallikreins: Essential epidermal messengers for regulation of the skin microenvironment during homeostasis, repair and disease

As the outermost layer of the skin, the epidermis is playing a major role in organism homeostasis providing the first barrier against external aggressions. Although considered as an extracellular matrix (ECM)-poor subtissue, the epidermal microenvironment is a key regulator of skin homeostasis and functionality. Among the proteins essential for upholding the epidermal microenvironment are the members of the kallikrein (KLK) family composed of 15 secreted serine proteases. Most of the members of these epithelial-specific proteins are present in skin and regulate skin desquamation and inflammation. However, although epidermal products, the consequences of KLK activities are not confined to the epidermis but widespread in the skin. In this review starting with the location and proteolytic activation cascade of KLKs, we present KLKs involvement in skin homeostasis, regeneration and pathology. KLKs have a large variety of substrates including ECM proteins, and evidence suggests that they are involved in the different steps of skin wound healing as discussed here. KLKs are also used as prognosis/diagnosis markers for many cancer types and we are focusing later on KLKs in cutaneous cancers, although their pathogenicity remains to be fully elucidated. Dysregulation of the KLK cascade is directly responsible for skin diseases with heavy inflammatory aspects, highlighting their involvement in skin immune homeostasis. Future studies will be needed to support the therapeutic potential of adjusting KLK activities for treatment of inflammatory skin diseases and wound healing pathologies.

•

Regulation of the microenvironment even in an extracellular matrix-poor tissue can heavily impact organ function.

•

Extracellular activities of kallikreins maintain skin homeostasis by regulating desquamation and inflammation.

•

The activation of skin epidermal-specific kallikrein family of proteases is regulated by an intricate proteolytic cascade.

•

Kallikreins are emerging as players during skin wound healing.

•

Dysregulated kallikrein expression and activity occur in cancers and inflammatory skin diseases.

Biochemical pathways mediated by KLK6 protease in breast cancer

Kallikrein-related peptidase 6 (KLK6) is a serine protease normally expressed in mammary tissue and aberrantly regulated in breast cancer. At physiological levels, KLK6 functions as a suppressor of breast cancer, while its aberrant overexpression (> 50-fold higher than normal) is characteristic of a subset of breast cancers and has been linked to accelerated growth of primary breast tumors in severe combined immunodeficiency mice (Pampalakis et al. Cancer Res 2009, 69, 3779). Here, we investigated the molecular mechanisms underlying the concentration-dependent functions of KLK6 by comparing MDA-MB-231 stable transfectants expressing increasing levels of KLK6 in in vitro and in vivo tumorigenicity assays (soft agar, xenograft growth, tail vein metastasis). Quantitative proteomics was applied to identify proteins that are altered upon re-expression of KLK6 in MDA-MB-231 at normal or constitutive levels. Overexpression of KLK6 is associated with increased metastatic ability of breast cancer cells into lungs, increased expression of certain S100 proteins (S100A4, S100A11) and keratins (KRT), and downregulation of the apoptosis-related proteases CASP7 and CASP8, and RABs. On the other hand, KLK6 re-expression at physiological levels leads to inhibition of lung metastases associated with suppression of S100 proteins (S100A4, S100A10, S100A13, S100A16) and induced CASP7 and CASP8 expression. As this is the first report that KLK6 expression is associated with S100 proteins, caspases, RABs, and KRTs, we validated this finding in clinical datasets. By integrating proteomics and microarray data from breast cancer patients, we generated two composite scores, KLK6 + S100B-S100A7 and KLK6 + S100B-S100A14-S100A16, to predict long-term survival of breast cancer patients. We present previously unknown pathways implicating KLK6 in breast cancer. The findings promise to aid our understanding of the functional roles of KLK6 in breast cancer and may yield new biomarkers for the cancer types in which KLK6 is known to be aberrantly upregulated.

A proinflammatory role of KLK6 protease in Netherton syndrome

BACKGROUND

Netherton syndrome (NS) is a rare but severe type of ichthyosis characterized by atopy, allergies, and potentially lethal skin overdesquamation associated with highly elevated proteolytic activities in LEKTI-deficient epidermis. NS symptoms are recapitulated in Spink5^-/- mouse where the gene encoding Lekti has been invalidated. Spink5^-/- mice die within 5h from birth due to their severe skin barrier defect leading to dehydration. Spink5^-/- mice also serve as a model for atopic dermatitis. The KLK6 protease is expressed by epidermal keratinocytes and shown in vitro to cleave desmosomal components.

OBJECTIVE

To investigate in vivo whether KLK6 is implicated in epidermal overdesquamation and/or inflammation associated with NS.

METHODS

The role of KLK6 was evaluated by generating Spink5^-/-Klk6^-/- double knockout mice. The phenotype was assessed by macroscopic observation, immunohistochemistry for differentiation markers, in situ zymography for proteolysis, and quantification of proinflammatory cytokines.

RESULTS

Elimination of Klk6 in Spink5^-/- remarkably suppresses the expression of Tslp, a major itching-inducing factor and driver of allergic reactions. Tnfα and the Th17 promoting cytokine Il-23 were also suppressed. Spink5^-/-Klk6^-/- mice display normalized keratinocyte differentiation, nevertheless, epidermal proteolytic activities and the associated overdesquamation were not ameliorated, and Spink5^-/-Klk6^-/- still died from a severe epidermal barrier defect as the Spink5^-/-.

CONCLUSIONS

Ablation of Klk6 largely suppresses epidermal inflammation but cannot rescue overdesquamation leading to the lethal NS phenotype. Nonetheless, our findings demonstrate for the first time that KLK6 is implicated in skin inflammation and may represent a novel druggable target for NS and other inflammatory conditions e.g. atopic dermatitis.

OMICfpp: a fuzzy approach for paired RNA-Seq counts

BACKGROUND

RNA sequencing is a widely used technology for differential expression analysis. However, the RNA-Seq do not provide accurate absolute measurements and the results can be different for each pipeline used. The major problem in statistical analysis of RNA-Seq and in the omics data in general, is the small sample size with respect to the large number of variables. In addition, experimental design must be taken into account and few tools consider it.

RESULTS

We propose OMICfpp, a method for the statistical analysis of RNA-Seq paired design data. First, we obtain a p-value for each case-control pair using a binomial test. These p-values are aggregated using an ordered weighted average (OWA) with a given orness previously chosen. The aggregated p-value from the original data is compared with the aggregated p-value obtained using the same method applied to random pairs. These new pairs are generated using between-pairs and complete randomization distributions. This randomization p-value is used as a raw p-value to test the differential expression of each gene. The OMICfpp method is evaluated using public data sets of 68 sample pairs from patients with colorectal cancer. We validate our results through bibliographic search of the reported genes and using simulated data set. Furthermore, we compared our results with those obtained by the methods edgeR and DESeq2 for paired samples. Finally, we propose new target genes to validate these as gene expression signatures in colorectal cancer. OMICfpp is available at http://www.uv.es/ayala/software/OMICfpp_0.2.tar.gz .

CONCLUSIONS

Our study shows that OMICfpp is an accurate method for differential expression analysis in RNA-Seq data with paired design. In addition, we propose the use of randomized p-values pattern graphic as a powerful and robust method to select the target genes for experimental validation.