Precision transcriptomics of viral foci reveals the spatial regulation of immune-signaling genes and identifies RBOHD as an important player in the incompatible interaction between potato virus Y and potato.

Whereas the activation of resistance (R) proteins has been intensively studied, the downstream signaling mechanisms leading to the restriction of the pathogen remain mostly unknown. We studied the immunity network response conditioned by the potato Ny-1 gene against potato virus Y. We analyzed the processes in the cell death zone and surrounding tissue on the biochemical and gene expression levels in order to reveal the spatiotemporal regulation of the immune response. We show that the transcriptional response in the cell death zone and surrounding tissue is dependent on salicylic acid (SA). For some genes the spatiotemporal regulation is completely lost in the SA-deficient line, whereas other genes show a different response, indicating multiple connections between hormonal signaling modules. The induction of NADPH oxidase RBOHD expression occurs specifically on the lesion border during the resistance response. In plants with silenced RBOHD, the functionality of the resistance response is perturbed and the spread of the virus is not arrested at the site of infection. RBOHD is required for the spatial accumulation of SA, and conversely RBOHD is under the transcriptional regulation of SA. Using spatially resolved RNA-seq, we also identified spatial regulation of an UDP-glucosyltransferase, another component in feedback activation of SA biosynthesis, thus deciphering a novel aspect of resistance signaling.

Novel algorithm including CA-125, HE4 and body mass index in the diagnosis of endometrial cancer.

OBJECTIVES: To evaluate the diagnostic and prognostic potential of preoperative serum CA-125 and HE4 levels in patients with endometrial cancer. METHODS: Prospective case-control study of 133 women who underwent surgical treatment at the University Medical Centre Ljubljana (64 patients with endometrial cancer, 69 control patients with prolapsed uterus or myoma). Serum CA-125 and HE4 levels were determined using electrochemiluminescent assays. RESULTS: Serum CA-125 and HE4 levels were significantly higher in patients with endometrial cancer, compared to the controls (p=2.67×10-4, 1.36×10-7, respectively). A diagnostic model that combines serum CA-125 and HE4 levels and body mass index separated patients with endometrial cancer from controls, with AUC of 0.804, sensitivity of 66.7%, and specificity of 84.6%. Serum HE4 levels showed good prognostic potential and stratified the patients according to presence/absence of deep myometrial invasion (p=0.001) or lymphovascular invasion (p=0.003), with AUCs of 0.78 and 0.81, respectively. In low-risk patients with grade 1 and 2 endometrioid cancer for whom lymphadenectomy can be avoided, HE4 allowed stratification according to deep myometrial invasion (p=3.39×10-4), with AUC of 0.84. Although median HE4 levels were higher in patients with lymphovascular invasion, this difference did not reach significance (p=0.06). CONCLUSIONS: A model based on preoperative serum CA-125 and HE4 levels and body mass index has good diagnostic accuracy for separation of patients with endometrial cancer and control patients. In patients with endometrial cancer, serum HE4 levels allow prediction of deep myometrial and lymphovascular invasion.

Inter-laboratory assessment of different digital PCR platforms for quantification of human cytomegalovirus DNA.

Quantitative PCR (qPCR) is an important tool in pathogen detection. However, the use of different qPCR components, calibration materials and DNA extraction methods reduces comparability between laboratories, which can result in false diagnosis and discrepancies in patient care. The wider establishment of a metrological framework for nucleic acid tests could improve the degree of standardisation of pathogen detection and the quantification methods applied in the clinical context. To achieve this, accurate methods need to be developed and implemented as reference measurement procedures, and to facilitate characterisation of suitable certified reference materials. Digital PCR (dPCR) has already been used for pathogen quantification by analysing nucleic acids. Although dPCR has the potential to provide robust and accurate quantification of nucleic acids, further assessment of its actual performance characteristics is needed before it can be implemented in a metrological framework, and to allow adequate estimation of measurement uncertainties. Here, four laboratories demonstrated reproducibility (expanded measurement uncertainties below 15%) of dPCR for quantification of DNA from human cytomegalovirus, with no calibration to a common reference material. Using whole-virus material and extracted DNA, an intermediate precision (coefficients of variation below 25%) between three consecutive experiments was noted. Furthermore, discrepancies in estimated mean DNA copy number concentrations between laboratories were less than twofold, with DNA extraction as the main source of variability. These data demonstrate that dPCR offers a repeatable and reproducible method for quantification of viral DNA, and due to its satisfactory performance should be considered as candidate for reference methods for implementation in a metrological framework.

Applicability of human osteoarthritic chondrocytes for in vitro efficacy testing of anti-TNFα drugs.

In vitro cell-based models are important tools for assessing efficacies of new leads in early phases of drug development. Human osteoarthritic chondrocytes (OACs), obtained from biomedical waste material, represent a valuable, relatively accessible cellular source that could be used for this purpose. By employing reverse transcription-polymerase chain reaction (qRT-PCR) we compared gene expression profiles of key anabolic, catabolic and inflammatory genes of freshly isolated vs. monolayer cultured OACs (passages P0-P2) and non-stimulated vs. tumor necrosis factor alpha (TNF-α) stimulated P2 OACs. After expansion of OACs in monolayer cultures, the expression of almost all analyzed genes significantly decreased. The subsequent addition of TNF-α to OACs at P2 significantly increased expressions of all catabolic and inflammatory genes, leaving the anabolic profile almost unchanged. TNF-α-treated OACs were later utilized for efficacy testing of anti-TNF-α drugs infliximab and etanercept and both significantly reduced the expressions of all catabolic and inflammatory genes tested.

Antibacterial Activity of Wild Mushroom Extracts on Bacterial Wilt Pathogen Ralstonia solanacearum.

In total, 150 protein extracts from 94 different basidiomycete and ascomycete wild mushroom species were tested for antibacterial activity against the quarantine plant-pathogen bacterium Ralstonia solanacearum. In in vitro microtiter plate assays, 15 extracts with moderate to high antibacterial activities were identified: 11 completely inhibited bacterial growth and 4 showed partial inhibition. Of these 15 extracts, 5 were further tested and 3 extracts slowed disease progression and reduced disease severity in artificially inoculated tomato and potato plants. However, the in vitro activities of the extracts did not always correlate with their in vivo activities, which emphasizes the importance of performing early screening tests also in vivo. Testing of selected extracts against 12 R. solanacearum strains identified 6 with potential for broader applicability. Further analysis of extracts from Amanita phalloides and Clitocybe geotropa showed that the active substances are proteins with an approximate size of 180 kDa. To our knowledge, this is the first in vitro and in vivo study that demonstrates that mushroom protein extracts can be promising for treatment of bacterial wilt caused by R. solanacearum.

Bimodal dynamics of primary metabolism-related responses in tolerant potato-Potato virus Y interaction.

BACKGROUND: Potato virus Y (PVY) is a major pathogen that causes substantial economic losses in worldwide potato production. Different potato cultivars differ in resistance to PVY, from severe susceptibility, through tolerance, to complete resistance. The aim of this study was to better define the mechanisms underlying tolerant responses of potato to infection by the particularly aggressive PVY(NTN) strain. We focused on the dynamics of the primary metabolism-related processes during PVY(NTN) infection. RESULTS: A comprehensive analysis of the dynamic changes in primary metabolism was performed, which included whole transcriptome analysis, nontargeted proteomics, and photosynthetic activity measurements in potato cv. Désirée and its transgenic counterpart depleted for accumulation of salicylic acid (NahG-Désirée). Faster multiplication of virus occurred in the NahG-Désirée, with these plants developing strong disease symptoms. We show that while the dynamics of responses at the transcriptional level are extensive and bimodal, this is only partially translated to the protein level, and to the final functional outcome. Photosynthesis-related genes are transiently induced before viral multiplication is detected and it is down-regulated later on. This is reflected as a deficiency of the photosynthetic apparatus at the onset of viral multiplication only. Interestingly, specific and constant up-regulation of some RuBisCO transcripts was detected in Désirée plants, which might be important, as these proteins have been shown to interact with viral proteins. In SA-deficient and more sensitive NahG-Désirée plants, consistent down-regulation of photosynthesis-related genes was detected. A constant reduction in the photochemical efficiency from the onset of viral multiplication was identified; in nontransgenic plants this decrease was only transient. The transient reduction in net photosynthetic rate occurred in both genotypes with the same timing, and coincided with changes in stomatal conductivity. CONCLUSIONS: Down-regulation of photosynthesis-related gene expression and decreased photosynthetic activity is in line with other studies that have reported the effects of biotic stress on photosynthesis. Here, we additionally detected induction of light-reaction components in the early stages of PVY(NTN) infection of tolerant interaction. As some of these components have already been shown to interact with viral proteins, their overproduction might contribute to the absence of symptoms in cv. Désirée.

Transcriptome study and identification of potential marker genes related to the stable expression of recombinant proteins in CHO clones.

BACKGROUND: Chinese hamster ovary (CHO) cells have become the host of choice for the production of recombinant proteins, due to their capacity for correct protein folding, assembly, and posttranslational modifications. The most widely used system for recombinant proteins is the gene amplification procedure that uses the CHO-Dhfr expression system. However, CHO cells are known to have a very unstable karyotype. This is due to chromosome rearrangements that can arise from translocations and homologous recombination, especially when cells with the CHO-Dhfr expression system are treated with methotrexate hydrate. The present method used in the industry for testing clones for their long-term stability of recombinant protein production is empirical, and it involves their cultivation over extended periods of time prior to the selection of the most suitable clone for further bioprocess development. The aim of the present study was the identification of marker genes that can predict stable expression of recombinant genes in particular clones early in the development stage. RESULTS: The transcriptome profiles of CHO clones with stable and unstable recombinant protein production were investigated over 10-weeks of cultivation, using a DNA microarray. We identified 14 genes that were differentially expressed between the stable and unstable clones already at 2 weeks from the beginning of the cultivation. Their expression was validated by reverse-transcription quantitative real-time PCR (RT-qPCR). Furthermore, the k-nearest neighbour algorithm approach shows that the combination of the gene expression patterns of only five of these 14 genes is sufficient to predict stable recombinant protein production in clones in the early phases of cell-line development. CONCLUSIONS: The exact molecular mechanisms that cause unstable recombinant protein production are not fully understood. However, the expression profiles of some genes in clones with stable and unstable recombinant protein production allow prediction of such instability early in the cell-line development stage. We have thus developed a proof-of-concept for a novel approach to eliminate unstable clones in the CHO-Dhfr expression system, which saves time and labour-intensive work in cell-line development.

Gene expression of cultured human chondrocytes as a model for assessing neutralization efficacy of soluble TNFα by TNFα antagonists.

Tumor necrosis factor-alpha (TNFα) antagonists are efficacious in the treatment of various immune-mediated inflammatory diseases. Because of rapidly growing demand for developing new or biosimilar versions of these biologicals, the need to create in vitro testing models that best represent physiological conditions is increasing. Primary human chondrocytes were used for potency evaluation and comparison between the molecular effects of anti-TNFα biologicals. Infliximab and etanercept were chosen to assess the suitability of chondrocyte cell culture for determination of anti-TNFα neutralization efficacy employing quantitative reverse transcription-polymerase chain reaction (qRT-PCR) technology. Use of both anti-TNFα biologics resulted in decrease of TNFα-stimulated expression of various matrix metalloproteinases, interleukins and other inflammation-related genes in our cell model. Significant differences in inhibition efficacy of etanercept and infliximab were observed, which were confirmed also on protein level. To evaluate the potency of anti-TNFα biologicals, a selection of TNFα-responsive target genes was made from the gene array data. The selected genes were employed in development of statistical model, which enables comparability of anti-TNFα biologicals. The presented analytical approach is suitable for assessment of the neutralization efficacy of various anti-TNFα biologicals. As such, it can be used for additional comprehensive characterization and comparability of TNF antagonists in preclinical drug testing.

Transcriptomic variation between different Chinese hamster ovary cell lines.

OBJECTIVES: To identify transcription markers that uniquely determine specific Chinese hamster ovary (CHO) cell lines and can be used for the identification of cell lines in the process of biopharmaceutical cell-line development. RESULTS: Five CHO cell lines with different origins were extensively characterised at the transcriptomic level and the results were compared to their karyotype characterisation. The analysed cell lines differ in their karyotype but, due to the genome instability observed during parental and recombinant cell-line establishment, karyotyping is not the preferred method for accurate identification of the various CHO cell lines. Marker genes unique to a specific cell line were identified by microarrays, and their expression was validated by reverse-transcription quantitative real-time PCR. The analysed cell lines can be differentiated by the presence/absence of detectable marker gene expression. Additionally, the similarity of the transcriptional profiles is dependent on cell-line history but independent of the manipulation steps involved in the recombinant cell-line development process. CONCLUSIONS: Certain transcripts can be used as markers for the identification of a CHO cell line undergoing recombinant development and thus represent a powerful tool for ensuring the maintenance of high quality standards.

Identification of plasma biomarker candidates in glioblastoma using an antibody-array-based proteomic approach.

BACKGROUND: Glioblastoma multiforme (GBM) is a brain tumour with a very high patient mortality rate, with a median survival of 47 weeks. This might be improved by the identification of novel diagnostic, prognostic and predictive therapy-response biomarkers, preferentially through the monitoring of the patient blood. The aim of this study was to define the impact of GBM in terms of alterations of the plasma protein levels in these patients. MATERIALS AND METHODS: We used a commercially available antibody array that includes 656 antibodies to analyse blood plasma samples from 17 healthy volunteers in comparison with 17 blood plasma samples from patients with GBM. RESULTS: We identified 11 plasma proteins that are statistically most strongly associated with the presence of GBM. These proteins belong to three functional signalling pathways: T-cell signalling and immune responses; cell adhesion and migration; and cell-cycle control and apoptosis. Thus, we can consider this identified set of proteins as potential diagnostic biomarker candidates for GBM. In addition, a set of 16 plasma proteins were significantly associated with the overall survival of these patients with GBM. Guanine nucleotide binding protein alpha (GNAO1) was associated with both GBM presence and survival of patients with GBM. CONCLUSIONS: Antibody array analysis represents a useful tool for the screening of plasma samples for potential cancer biomarker candidates in small-scale exploratory experiments; however, clinical validation of these candidates requires their further evaluation in a larger study on an independent cohort of patients.

Niche-based mechanisms operating within extreme habitats: a case study of subterranean amphipod communities.

It has been suggested that both niche-based and neutral mechanisms are important for biological communities to evolve and persist. For communities in extreme and isolated environments such as caves, theoretical and empirical considerations (low species turnover, high stress, strong convergence owing to strong directional selection) predict neutral mechanisms and functional equivalence of species. We tested this prediction using subterranean amphipod communities from caves and interstitial groundwater. Contrary to expectations, functional morphological diversity within communities in both habitats turned out to be significantly higher than the null model of randomly assembled communities. This suggests that even the most extreme, energy-poor environments still maintain the potential for diversification via differentiation of niches.

Gene Expression Profiling of Recombinant Protein Producing E. coli at Suboptimal Growth Temperature.

Recent studies have revealed that at lower cultivation temperatures (25 °C) much higher percentage of correctly folded recombinant hG-CSF protein can be extracted from inclusion bodies. Hence, the goal of our research was to investigate mechanisms determining characteristics of non-classical inclusion bodies production using gene expression profiling, focusing on proteases and chaperones gene expression. Statistical analysis of microarray data showed prominent changes in energy metabolism, in metabolism of amino acids and nucleotides, as well as in biosynthesis of cofactors and secondary metabolites if the culture was grown below its optimal temperature. Moreover, 24 differentially expressed up to now known genes classified among proteases, chaperones and other heat or stress related genes. Among chaperones UspE and among proteases YaeL and YeaZ might play an important role in accumulation of correctly folded recombinant proteins. Membrane localized protease yaeL gene was found to have higher activity at 25 °C and is thus potentially functionally related to the more efficient recombinant protein production at lower temperatures. The results of this study represent advance in the understanding of recombinant protein production in E. coli. Genes potentially influencing production of recombinant protein at lower growth temperature represent basis for further research towards improvement of E. coli production strains as well as fermentation process.

pISA-tree - a data management framework for life science research projects using a standardised directory tree.

We developed pISA-tree, a straightforward and flexible data management solution for organisation of life science project-associated research data and metadata. pISA-tree was initiated by end-user requirements thus its strong points are practicality and low maintenance cost. It enables on-the-fly creation of enriched directory tree structure (project/Investigation/Study/Assay) based on the ISA model, in a standardised manner via consecutive batch files. Templates-based metadata is generated in parallel at each level enabling guided submission of experiment metadata. pISA-tree is complemented by two R packages, pisar and seekr. pisar facilitates integration of pISA-tree datasets into bioinformatic pipelines and generation of ISA-Tab exports. seekr enables synchronisation with the FAIRDOMHub repository. Applicability of pISA-tree was demonstrated in several national and international multi-partner projects. The system thus supports findable, accessible, interoperable and reusable (FAIR) research and is in accordance with the Open Science initiative. Source code and documentation of pISA-tree are available at https://github.com/NIB-SI/pISA-tree .

Hay meadow vibroscape and interactions within insect vibrational community.

Our experiences shape our knowledge and understanding of the world around us. The natural vibrational environment (vibroscape) is hidden to human senses but is nevertheless perceived and exploited by the majority of animals. Here, we show that the vibroscape recorded on plants in a temperate hay meadow is a dynamic low-frequency world, rich in species-specific vibrational signals. The overall vibroscape composition changed throughout the season and also depended on the plant species, as well as on the spatial position of individual plants within the meadow. Within the studied community, vibrationally signaling species sharing this communication channel avoided interference primarily by partitioning vibrational space on a fine temporal scale. The vibroscape is a reliable source of information in the environment and expands our understanding of ecological and evolutionary processes.

Development of sampling approaches for the determination of the presence of genetically modified organisms at the field level.

In order to comply with the European Union regulatory threshold for the adventitious presence of genetically modified organisms (GMOs) in food and feed, it is important to trace GMOs from the field. Appropriate sampling methods are needed to accurately predict the presence of GMOs at the field level. A 2-year field experiment with two maize varieties differing in kernel colour was conducted in Slovenia. Based on the results of data mining analyses and modelling, it was concluded that spatial relations between the donor and receptor field were the most important factors influencing the distribution of outcrossing rate (OCR) in the field. The approach for estimation fitting function parameters in the receptor (non-GM) field at two distances from the donor (GM) field (10 and 25 m) for estimation of the OCR (GMO content) in the whole receptor field was developed. Different sampling schemes were tested; a systematic random scheme in rows was proposed to be applied for sampling at the two distances for the estimation of fitting function parameters for determination of OCR. The sampling approach had already been validated with some other OCR data and was practically applied in the 2009 harvest in Poland. The developed approach can be used for determination of the GMO presence at the field level and for making appropriate labelling decisions. The importance of this approach lies in its possibility to also address other threshold levels beside the currently prescribed labelling threshold of 0.9% for food and feed.

Evaluation of the reliability of maize reference assays for GMO quantification.

A reliable PCR reference assay for relative genetically modified organism (GMO) quantification must be specific for the target taxon and amplify uniformly along the commercialised varieties within the considered taxon. Different reference assays for maize (Zea mays L.) are used in official methods for GMO quantification. In this study, we evaluated the reliability of eight existing maize reference assays, four of which are used in combination with an event-specific polymerase chain reaction (PCR) assay validated and published by the Community Reference Laboratory (CRL). We analysed the nucleotide sequence variation in the target genomic regions in a broad range of transgenic and conventional varieties and lines: MON 810 varieties cultivated in Spain and conventional varieties from various geographical origins and breeding history. In addition, the reliability of the assays was evaluated based on their PCR amplification performance. A single base pair substitution, corresponding to a single nucleotide polymorphism (SNP) reported in an earlier study, was observed in the forward primer of one of the studied alcohol dehydrogenase 1 (Adh1) (70) assays in a large number of varieties. The SNP presence is consistent with a poor PCR performance observed for this assay along the tested varieties. The obtained data show that the Adh1 (70) assay used in the official CRL NK603 assay is unreliable. Based on our results from both the nucleotide stability study and the PCR performance test, we can conclude that the Adh1 (136) reference assay (T25 and Bt11 assays) as well as the tested high mobility group protein gene assay, which also form parts of CRL methods for quantification, are highly reliable. Despite the observed uniformity in the nucleotide sequence of the invertase gene assay, the PCR performance test reveals that this target sequence might occur in more than one copy. Finally, although currently not forming a part of official quantification methods, zein and SSIIb assays are found to be highly reliable in terms of nucleotide stability and PCR performance and are proposed as good alternative targets for a reference assay for maize.

'Bois noir' phytoplasma induces significant reprogramming of the leaf transcriptome in the field grown grapevine.

BACKGROUND: Phytoplasmas are bacteria without cell walls from the class Mollicutes. They are obligate intracellular plant pathogens which cause diseases in hundreds of economically important plants including the grapevine (Vitis vinifera). Knowledge of their biology and the mechanisms of their interactions with hosts is largely unknown because they are uncultivable and experimentally inaccessible in their hosts. We detail here the global transcriptional profiling in grapevine responses to phytoplasmas. The gene expression patterns were followed in leaf midribs of grapevine cv. 'Chardonnay' naturally infected with a phytoplasma from the stolbur group 16SrXII-A, which is associated with the grapevine yellows disease 'Bois noir'. RESULTS: We established an on field experimental system in a productive vineyard that allowed application of molecular tools in a plant natural environment. Global transcription profiles of infected samples were compared with the healthy ones using microarray datasets and metabolic pathway analysis software (MapMan). The two-year-long experiment revealed that plant genes involved in primary and secondary metabolic pathways were changed in response to infection and that these changes might support phytoplasma nutrition. A hypothesis that phytoplasmas interact with the plant carbohydrate metabolism was proven and some possibilities how the products of this pathway might be utilized by phytoplasmas are discussed. In addition, several photosynthetic genes were largely down-regulated in infected plants, whereas defense genes from the metabolic pathway leading to formation of flavonoids and some PR proteins were significantly induced. Few other genes involved in defense-signaling were differentially expressed in healthy and infected plants. A set of 17 selected genes from several differentially expressed pathways was additionally analyzed with quantitative real-time PCR and confirmed to be suitable for a reliable classification of infected plants and for the characterization of susceptibility features in the field conditions. CONCLUSION: This study revealed some fundamental aspects of grapevine interactions with the stolbur 'Bois noir' phytoplasma in particular and some plant interactions with phytoplasmas in general. In addition, the results of the study will likely have an impact on grape improvement by yielding marker genes that can be used in new diagnostic assays for phytoplasmas or by identifying candidate genes that contribute to the improved properties of grape.

DiNAR: revealing hidden patterns of plant signalling dynamics using Differential Network Analysis in R.

BACKGROUND: Progress in high-throughput molecular methods accompanied by more complex experimental designs demands novel data visualisation solutions. To specifically answer the question which parts of the specifical biological system are responding in particular perturbation, integrative approach in which experimental data are superimposed on a prior knowledge network is shown to be advantageous. RESULTS: We have developed DiNAR, Differential Network Analysis in R, a user-friendly application with dynamic visualisation that integrates multiple condition high-throughput data and extensive biological prior knowledge. Implemented differential network approach and embedded network analysis allow users to analyse condition-specific responses in the context of topology of interest (e.g. immune signalling network) and extract knowledge concerning patterns of signalling dynamics (i.e. rewiring in network structure between two or more biological conditions). We validated the usability of software on the Arabidopsis thaliana and Solanum tuberosum datasets, but it is set to handle any biological instances. CONCLUSIONS: DiNAR facilitates detection of network-rewiring events, gene prioritisation for future experimental design and allows capturing dynamics of complex biological system. The fully cross-platform Shiny App is hosted and freely available at https://nib-si.shinyapps.io/DiNAR. The most recent version of the source code is available at https://github.com/NIB-SI/DiNAR/ with a DOI 10.5281/zenodo.1230523 of the archived version in Zenodo.

Localization patterns of cathepsins K and X and their predictive value in glioblastoma.

Background Glioblastoma is a highly aggressive central nervous system neoplasm characterized by extensive infiltration of malignant cells into brain parenchyma, thus preventing complete tumor eradication. Cysteine cathepsins B, S, L and K are involved in cancer progression and are overexpressed in glioblastoma. We report here for the first time that cathepsin X mRNA and protein are also abundantly present in malignant glioma. Materials and methods Gene expression of cathepsins K and X was analyzed using publically-available tran-scriptomic datasets and correlated with glioma grade and glioblastoma subtype. Kaplan-Maier survival analysis was performed to evaluate the predictive value of cathepsin K and X mRNA expression. Cathepsin protein expression was localized and semi-quantified in tumor tissues by immunohistochemistry. Results Highest gene expression of cathepsins K and X was found in glioblastoma, in particular in the mesenchymal subtype. Overall, high mRNA expression of cathepsin X, but not that of cathepsin K, correlated with poor patients' survival. Cathepsin K and X proteins were abundantly and heterogeneously expressed in glioblastoma tissue. Immuno-labeling of cathepsins K and X was observed in areas of CD133-positive glioblastoma stem cells, localized around arterioles in their niches that also expressed SDF-1α and CD68. mRNA levels of both cathepsins K and X correlated with mRNA levels of markers of glioblastoma stem cells and their niches. Conclusions The presence of both cathepsins in glioblastoma stem cell niche regions indicates their possible role in regulation of glioblastoma stem cell homing in their niches. The clinical relevance of this data needs to be elaborated in further prospective studies.