Search PubMed⌕ Search

Biomedical subjects

Roland Eils

Publications and source records attributed to Roland Eils.

At least 19 recordsLinked to original sources

Beyond benchmarking: an expert-guided consensus approach to spatially aware clustering.

Spatial omics technologies have revolutionized the study of tissue architecture and cellular heterogeneity by integrating molecular profiles with spatial localization. In spatially resolved transcriptomics, delineating higher-order anatomical structures is critical for understanding how cellular organization affects function. However, the reliability of current benchmarks of spatially aware clustering (SAC) methods is undermined by their narrow focus on Visium and brain tissue datasets and the incorrect interpretation of manual annotation as ground truth. Here we present SACCELERATOR, a community-driven, extensible framework that standardizes data formatting, method integration and metric evaluation, enabling rapid inclusion of new methods and datasets. Our analysis revealed substantial limitations in the generalizability and reproducibility of SAC methods and shows that anatomical labels commonly used as ground truths are often biased, error prone and unsuitable for benchmarking. Rather than ranking methods, we propose a consensus-guided workflow where descriptive spatial metrics highlight high-entropy regions of method disagreement, enabling targeted feedback for tissue experts. Applied to brain and cancer datasets, this approach uncovered biologically meaningful patterns overlooked by individual SAC methods and manual annotations, highlighting the need for iterative, expert-in-the-loop evaluation.

Benchmarking↗

Managing workflow executions with WESkit.

SUMMARY: In biomedical research, managing computational workflows across numerous projects-with varying parameters, tools, and environments-creates major challenges in scalability, reproducibility, and collaboration. Here, we present WESkit, an implementation of the Global Alliance for Genomics and Health (GA4GH) Workflow Execution Service (WES) interface, designed to streamline the execution, monitoring, and documentation of data processing workflows. It addresses the complexities involved in managing numerous executions with varying parameters across diverse research projects. Supporting both Snakemake and Nextflow, the system enables consistent automation and centralized monitoring, which benefits research groups aiming for long-term reproducibility and scalable collaboration. Its suitability for larger teams and service units is further enhanced by seamless integration into cloud environments, contributing to the GA4GH cloud framework. AVAILABILITY AND IMPLEMENTATION: The software WESkit is available under MIT license at the GitLab repository (https://gitlab.com/one-touch-pipeline/weskit). The WESkit main repository is archived at Software Heritage (https://archive.softwareheritage.org/browse/origin/directory/?origin_url=https://gitlab.com/one-touch-pipeline/weskit/api.git) and can be found using "one-touch-pipeline/weskit" term in the search section.

Workflow↗

Classification of neuroblastoma patients by published gene-expression markers reveals a low sensitivity for unfavorable courses of MYCN non-amplified disease.

Currently, Pubmed lists 385 marker genes for neuroblastoma outcome. Using a customized neuroblastoma-microarray, we evaluated the prognostic impact of the gene-expression pattern of 349 of these candidates (90.6%) in 127 neuroblastoma patients with divergent outcome. By significance analysis of microarrays (SAM) and both uncorrected and Bonferroni-corrected ANOVA, 166/349 (47.5%), 218/349 (62.5%) and 128/349 (36.4%) candidates showed significant differential expression between patients with contrasting outcome. By Prediction Analysis for Microarrays (PAM), a 38-gene-classifier was derived from all markers, which classified patients outcome with an overall accuracy of 78.5%. However, patients with unfavorable outcome of MYCN non-amplified disease were largely misclassified (accuracy: 35%), suggesting that these courses are not identified by current marker genes.

Analysis of Variance↗

Customized oligonucleotide microarray gene expression-based classification of neuroblastoma patients outperforms current clinical risk stratification.

PURPOSE: To develop a gene expression-based classifier for neuroblastoma patients that reliably predicts courses of the disease. PATIENTS AND METHODS: Two hundred fifty-one neuroblastoma specimens were analyzed using a customized oligonucleotide microarray comprising 10,163 probes for transcripts with differential expression in clinical subgroups of the disease. Subsequently, the prediction analysis for microarrays (PAM) was applied to a first set of patients with maximally divergent clinical courses (n = 77). The classification accuracy was estimated by a complete 10-times-repeated 10-fold cross validation, and a 144-gene predictor was constructed from this set. This classifier's predictive power was evaluated in an independent second set (n = 174) by comparing results of the gene expression-based classification with those of risk stratification systems of current trials from Germany, Japan, and the United States. RESULTS: The first set of patients was accurately predicted by PAM (cross-validated accuracy, 99%). Within the second set, the PAM classifier significantly separated cohorts with distinct courses (3-year event-free survival [EFS] 0.86 +/- 0.03 [favorable; n = 115] v 0.52 +/- 0.07 [unfavorable; n = 59] and 3-year overall survival 0.99 +/- 0.01 v 0.84 +/- 0.05; both P < .0001) and separated risk groups of current neuroblastoma trials into subgroups with divergent outcome (NB2004: low-risk 3-year EFS 0.86 +/- 0.04 v 0.25 +/- 0.15, P < .0001; intermediate-risk 1.00 v 0.57 +/- 0.19, P = .018; high-risk 0.81 +/- 0.10 v 0.56 +/- 0.08, P = .06). In a multivariate Cox regression model, the PAM predictor classified patients of the second set more accurately than risk stratification of current trials from Germany, Japan, and the United States (P < .001; hazard ratio, 4.756 [95% CI, 2.544 to 8.893]). CONCLUSION: Integration of gene expression-based class prediction of neuroblastoma patients may improve risk estimation of current neuroblastoma trials.

Biomarkers, Tumor↗

Tropical--parameter estimation and simulation of reaction-diffusion models based on spatio-temporal microscopy images.

UNLABELLED: Tropical is a software for simulation and parameter estimation of reaction-diffusion models. Based on spatio-temporal microscopy images, Tropical estimates reaction and diffusion coefficients for user-defined models. Tropical allows the investigation of systems with an inhomogeneous distribution of molecules, making it well suited for quantitative analyses of microscopy experiments such as fluorescence recovery after photobleaching (FRAP). AVAILABILITY: Tropical is available free of charge for academic use at http://www.dkfz.de/tbi/projects/modellingAndSimulationOfCelluarSystems/tropical.jsp after signing a material transfer agreement.

Algorithms↗

Group testing for pathway analysis improves comparability of different microarray datasets.

MOTIVATION: The wide use of DNA microarrays for the investigation of the cell transcriptome triggered the invention of numerous methods for the processing of microarray data and lead to a growing number of microarray studies that examine the same biological conditions. However, comparisons made on the level of gene lists obtained by different statistical methods or from different datasets hardly converge. We aimed at examining such discrepancies on the level of apparently affected biologically related groups of genes, e.g. metabolic or signalling pathways. This can be achieved by group testing procedures, e.g. over-representation analysis, functional class scoring (FCS), or global tests. RESULTS: Three public prostate cancer datasets obtained with the same microarray platform (HGU95A/HGU95Av2) were analyzed. Each dataset was subjected to normalization by either variance stabilizing normalization (vsn) or mixed model normalization (MMN). Then, statistical analysis of microarrays was applied to the vsn-normalized data and mixed model analysis to the data normalized by MMN. For multiple testing adjustment the false discovery rate was calculated and the threshold was set to 0.05. Gene lists from the same method applied to different datasets showed overlaps between 42 and 52%, while lists from different methods applied to the same dataset had between 63 and 85% of genes in common. A number of six gene lists obtained by the two statistical methods applied to the three datasets was then subjected to group testing by Fisher's exact test. Group testing by GSEA and global test was applied to the three datasets, as well. Fisher's exact test followed by global test showed more consistent results with respect to the concordance between analyses on gene lists obtained by different methods and different datasets than the GSEA. However, all group testing methods identified pathways that had already been described to be involved in the pathogenesis of prostate cancer. Moreover, pathways recurrently identified in these analyses are more likely to be reliable than those from a single analysis on a single dataset.

Biomarkers, Tumor↗

Predicting protein subcellular locations using hierarchical ensemble of Bayesian classifiers based on Markov chains.

BACKGROUND: The subcellular location of a protein is closely related to its function. It would be worthwhile to develop a method to predict the subcellular location for a given protein when only the amino acid sequence of the protein is known. Although many efforts have been made to predict subcellular location from sequence information only, there is the need for further research to improve the accuracy of prediction. RESULTS: A novel method called HensBC is introduced to predict protein subcellular location. HensBC is a recursive algorithm which constructs a hierarchical ensemble of classifiers. The classifiers used are Bayesian classifiers based on Markov chain models. We tested our method on six various datasets; among them are Gram-negative bacteria dataset, data for discriminating outer membrane proteins and apoptosis proteins dataset. We observed that our method can predict the subcellular location with high accuracy. Another advantage of the proposed method is that it can improve the accuracy of the prediction of some classes with few sequences in training and is therefore useful for datasets with imbalanced distribution of classes. CONCLUSION: This study introduces an algorithm which uses only the primary sequence of a protein to predict its subcellular location. The proposed recursive scheme represents an interesting methodology for learning and combining classifiers. The method is computationally efficient and competitive with the previously reported approaches in terms of prediction accuracies as empirical results indicate. The code for the software is available upon request.

Amino Acid Sequence↗

Gene expression signature predicting pathologic complete response with gemcitabine, epirubicin, and docetaxel in primary breast cancer.

PURPOSE: Primary systemic therapy (PST) with gemcitabine (G), epirubicin (E), and docetaxel (Doc) has resulted in a pathologic complete response (pCR) in 26% of primary breast cancer patients. This study was aimed at the identification of a gene expression signature in diagnostic core biopsy tissue samples that predicts pCR. PATIENTS AND METHODS: Core biopsy samples from patients with operable primary breast cancer, T2-4N0-2M0, enrolled onto two phase I and II trials evaluating GEDoc (n = 48) and GE sequentially followed by Doc (GEsDoc; n = 52) as PST were snap frozen and subjected to RNA expression profiling. A signature predicting pCR was discovered in the training set (GEsDoc) applying a support vector machine algorithm, and performance of this classifier was validated on the independent test set (GEDoc) by receiver operator characteristics analysis. RESULTS: We identified a signature consisting of 512 genes, which was enriched in genes involved in transforming growth factor beta and RAS-mediated signaling pathways, that predicts pCR with a sensitivity of 78%, a specificity of 90%, and an overall accuracy of 88% (95% CI, 75% to 95%). Apart from our signature, only HER2 overexpression was an independent predictor of pCR in multivariate analysis. CONCLUSION: In conclusion, our gene expression signature allows prediction of pCR to PST containing G, E, and Doc with unprecedented high overall accuracy and robustness.

Adult↗

GOPET: a tool for automated predictions of Gene Ontology terms.

BACKGROUND: Vast progress in sequencing projects has called for annotation on a large scale. A Number of methods have been developed to address this challenging task. These methods, however, either apply to specific subsets, or their predictions are not formalised, or they do not provide precise confidence values for their predictions. DESCRIPTION: We recently established a learning system for automated annotation, trained with a broad variety of different organisms to predict the standardised annotation terms from Gene Ontology (GO). Now, this method has been made available to the public via our web-service GOPET (Gene Ontology term Prediction and Evaluation Tool). It supplies annotation for sequences of any organism. For each predicted term an appropriate confidence value is provided. The basic method had been developed for predicting molecular function GO-terms. It is now expanded to predict biological process terms. This web service is available via http://genius.embnet.dkfz-heidelberg.de/menu/biounit/open-husar CONCLUSION: Our web service gives experimental researchers as well as the bioinformatics community a valuable sequence annotation device. Additionally, GOPET also provides less significant annotation data which may serve as an extended discovery platform for the user.

Artificial Intelligence↗

Discovering functional gene expression patterns in the metabolic network of Escherichia coli with wavelets transforms.

BACKGROUND: Microarray technology produces gene expression data on a genomic scale for an endless variety of organisms and conditions. However, this vast amount of information needs to be extracted in a reasonable way and funneled into manageable and functionally meaningful patterns. Genes may be reasonably combined using knowledge about their interaction behaviour. On a proteomic level, biochemical research has elucidated an increasingly complete image of the metabolic architecture, especially for less complex organisms like the well studied bacterium Escherichia coli. RESULTS: We sought to discover central components of the metabolic network, regulated by the expression of associated genes under changing conditions. We mapped gene expression data from E. coli under aerobic and anaerobic conditions onto the enzymatic reaction nodes of its metabolic network. An adjacency matrix of the metabolites was created from this graph. A consecutive ones clustering method was used to obtain network clusters in the matrix. The wavelet method was applied on the adjacency matrices of these clusters to collect features for the classifier. With a feature extraction method the most discriminating features were selected. We yielded network sub-graphs from these top ranking features representing formate fermentation, in good agreement with the anaerobic response of hetero-fermentative bacteria. Furthermore, we found a switch in the starting point for NAD biosynthesis, and an adaptation of the l-aspartate metabolism, in accordance with its higher abundance under anaerobic conditions. CONCLUSION: We developed and tested a novel method, based on a combination of rationally chosen machine learning methods, to analyse gene expression data on the basis of interaction data, using a metabolic network of enzymes. As a case study, we applied our method to E. coli under oxygen deprived conditions and extracted physiologically relevant patterns that represent an adaptation of the cells to changing environmental conditions. In general, our concept may be transferred to network analyses on biological interaction data, when data for two comparable states of the associated nodes are made available.

Algorithms↗

Transient colocalization of X-inactivation centres accompanies the initiation of X inactivation.

The initial differential treatment of the two X chromosomes during X-chromosome inactivation is controlled by the X-inactivation centre (Xic). This locus determines how many X chromosomes are present in a cell ('counting') and which X chromosome will be inactivated in female cells ('choice'). Critical control sequences in the Xic include the non-coding RNAs Xist and Tsix, and long-range chromatin elements. However, little is known about the process that ensures that X inactivation is triggered appropriately when more than one Xic is present in a cell. Using three-dimensional fluorescence in situ hybridization (FISH) analysis, we showed that the two Xics transiently colocalize, just before X inactivation, in differentiating female embryonic stem cells. Using Xic transgenes capable of imprinted but not random X inactivation, and Xic deletions that disrupt random X inactivation, we demonstrated that Xic colocalization is linked to Xic function in random X inactivation. Both long-range sequences and the Tsix element, which generates the antisense transcript to Xist, are required for the transient interaction of Xics. We propose that transient colocalization of Xics may be necessary for a cell to determine Xic number and to ensure the correct initiation of X inactivation.

Animals↗

High-resolution genomic profiling reveals association of chromosomal aberrations on 1q and 16p with histologic and genetic subgroups of invasive breast cancer.

PURPOSE: Invasive ductal carcinoma and invasive lobular carcinoma (ILC) represent the major histologic subtypes of invasive breast cancer. They differ with regard to presentation, metastatic spread, and epidemiologic features. To elucidate the genetic basis of these differences, we analyzed copy number imbalances that differentiate the histologic subtypes. EXPERIMENTAL DESIGN: High-resolution genomic profiling of 40 invasive breast cancers using matrix-comparative genomic hybridization with an average resolution of 0.5 Mb was conducted on bacterial artificial chromosome microarrays. The data were subjected to classification and unsupervised hierarchical cluster analyses. Expression of candidate genes was analyzed in tumor samples. RESULTS: The highest discriminating power was achieved when combining the aberration patterns of chromosome arms 1q and 16p, which were significantly more often gained in ILC. These regions were further narrowed down to subregions 1q24.2-25.1, 1q25.3-q31.3, and 16p11.2. Located within the candidate gains on 1q are two genes, FMO2 and PTGS2, known to be overexpressed in ILC relative to invasive ductal carcinoma. Assessment of four candidate genes on 16p11.2 by real-time quantitative PCR revealed significant overexpression of FUS and ITGAX in ILC with 16p copy number gain. Unsupervised hierarchical cluster analysis identified three molecular subgroups that are characterized by different aberration patterns, in particular concerning gain of MYC (8q24) and the identified candidate regions on 1q24.2-25.1, 1q25.3-q31.3, and 16p11.2. These genetic subgroups differed with regard to histology, tumor grading, frequency of alterations, and estrogen receptor expression. CONCLUSIONS: Molecular profiling using bacterial artificial chromosome arrays identified DNA copy number imbalances on 1q and 16p as significant classifiers of histologic and molecular subgroups.

Biomarkers, Tumor↗

A highly standardized, robust, and cost-effective method for genome-wide transcriptome analysis of peripheral blood applicable to large-scale clinical trials.

The use of peripheral blood mononuclear cells (PBMC) for transcriptome analysis has already been proven valuable for assessing disease-associated and drug-response-related gene signatures. While these proof-of-principle studies have been critically important, the instability of RNA within PBMC prohibits their use in large-scale multicenter trials for which samples have to be transported for a prolonged time prior to RNA isolation. Therefore, a prerequisite for transcriptome analysis of peripheral blood in clinical trials will be a standardized and valid method to stabilize the RNA profile immediately after blood withdrawal. Moreover, to be able to perform such large-scale clinical studies routinely in several hundred patients more cost-effective array technologies are required. To address these critical issues, we have combined a whole-blood RNA stabilization technology with a method to reduce globin mRNA, followed by genome-wide transcriptome analysis using a newly introduced BeadChip oligonucleotide technology. We demonstrate that the globin mRNA reduction method results in significantly improved data quality of stabilized RNA samples with low intragroup variance and a detection rate of expressed genes similar to that in PBMC. More important, even small differences in gene expression such as are observed between females and males were detected and sufficient to predict gender in whole-blood samples. We therefore propose the combination of globin mRNA reduction after whole-blood RNA stabilization with a newly introduced cost-effective BeadChip array as the preferred approach for large-scale multicenter trials, especially when establishing predictive markers for disease and treatment outcome in peripheral blood.

Blood Cells↗

Argonaute--a database for gene regulation by mammalian microRNAs.

MicroRNAs (miRNAs) constitute a recently discovered class of small non-coding RNAs that regulate expression of target genes either by decreasing the stability of the target mRNA or by translational inhibition. They are involved in diverse processes, including cellular differentiation, proliferation and apoptosis. Recent evidence also suggests their importance for cancerogenesis. By far the most important model systems in cancer research are mammalian organisms. Thus, we decided to compile comprehensive information on mammalian miRNAs, their origin and regulated target genes in an exhaustive, curated database called Argonaute (http://www.ma.uni-heidelberg.de/apps/zmf/argonaute/interface). Argonaute collects latest information from both literature and other databases. In contrast to current databases on miRNAs like miRBase::Sequences, NONCODE or RNAdb, Argonaute hosts additional information on the origin of an miRNA, i.e. in which host gene it is encoded, its expression in different tissues and its known or proposed function, its potential target genes including Gene Ontology annotation, as well as miRNA families and proteins known to be involved in miRNA processing. Additionally, target genes are linked to an information retrieval system that provides comprehensive information from sequence databases and a simultaneous search of MEDLINE with all synonyms of a given gene. The web interface allows the user to get information for a single or multiple miRNAs, either selected or uploaded through a text file. Argonaute currently has information on 839 miRNAs from human, mouse and rat.

Animals↗

Prediction of clinical outcome and biological characterization of neuroblastoma by expression profiling.

Neuroblastoma is a common childhood tumor comprising cases with rapid disease progression as well as spontaneous regression. Although numerous prognostic factors have been identified, risk evaluation in individual patients remains difficult. To define a reliable prognostic predictor and gene signatures characteristic of biological subgroups, we performed mRNA expression profiling of 68 neuroblastomas of all stages. Expression data were analysed using support vector machines (SVM-rbf), prediction analysis of microarrays (PAM), k-nearest neighbors (k-NN) algorithms and multiple decision trees. SVM-rbf performed best of all methods, and predicted recurrence of neuroblastoma with an accuracy of 85% (sensitivity 77%, specificity 94%). PAM identified a classifier of 39 genes reliably predicting outcome with an accuracy of 80%. In comparison, conventional risk stratification based on stage, age and MYCN-status only reached a predictive accuracy of 64%. Kaplan-Meier analysis using the PAM classifier indicated a 5-year survival of 20 versus 78% for patients with unfavorably versus favorably predicted neuroblastomas, respectively (P = 0.0001). Significance analysis of microarrays (SAM) identified additional genes differentially expressed among subgroups. MYCN-amplification and high expression of NTRK1/TrkA demonstrated a strong association with specific gene expression patterns. Our data suggest that microarray-derived data in addition to traditional clinical factors will be useful for risk assessment and defining biological properties of neuroblastoma.

Algorithms↗

Cross-platform analysis of cancer microarray data improves gene expression based classification of phenotypes.

BACKGROUND: The extensive use of DNA microarray technology in the characterization of the cell transcriptome is leading to an ever increasing amount of microarray data from cancer studies. Although similar questions for the same type of cancer are addressed in these different studies, a comparative analysis of their results is hampered by the use of heterogeneous microarray platforms and analysis methods. RESULTS: In contrast to a meta-analysis approach where results of different studies are combined on an interpretative level, we investigate here how to directly integrate raw microarray data from different studies for the purpose of supervised classification analysis. We use median rank scores and quantile discretization to derive numerically comparable measures of gene expression from different platforms. These transformed data are then used for training of classifiers based on support vector machines. We apply this approach to six publicly available cancer microarray gene expression data sets, which consist of three pairs of studies, each examining the same type of cancer, i.e. breast cancer, prostate cancer or acute myeloid leukemia. For each pair, one study was performed by means of cDNA microarrays and the other by means of oligonucleotide microarrays. In each pair, high classification accuracies (> 85%) were achieved with training and testing on data instances randomly chosen from both data sets in a cross-validation analysis. To exemplify the potential of this cross-platform classification analysis, we use two leukemia microarray data sets to show that important genes with regard to the biology of leukemia are selected in an integrated analysis, which are missed in either single-set analysis. CONCLUSION: Cross-platform classification of multiple cancer microarray data sets yields discriminative gene expression signatures that are found and validated on a large number of microarray samples, generated by different laboratories and microarray technologies. Predictive models generated by this approach are better validated than those generated on a single data set, while showing high predictive power and improved generalization performance.

Animals↗

Local gene density predicts the spatial position of genetic loci in the interphase nucleus.

Specific chromosomal translocations are hallmarks of many human leukemias. The basis for these translocation events is poorly understood, but it has been assumed that spatial positioning of genes in the nucleus of hematopoietic cells is a contributing factor. Analysis of the nuclear 3D position of the gene MLL, frequently involved in chromosomal translocations and five of its translocation partners (AF4, AF6, AF9, ENL and ELL), and two control loci revealed a characteristic radial distribution pattern in all hematopoietic cells studied. Genes in areas of high local gene density were found positioned towards the nuclear center, whereas genes in regions of low gene density were detected closer to the nuclear periphery. The gene density within a 2 Mbp window was found to be a better predictor for the relative positioning of a genomic locus within the cell nucleus than the gene density of entire chromosomes. Analysis of the position of MLL, AF4, AF6 and AF9 in cell lines carrying chromosomal translocations involving these genes revealed that the position of the normal genes was different from that of the fusion genes, and this was again consistent with the changes in local gene density within a 2 Mbp window. Thus, alterations in gene density directly at translocation junctions could explain the change in the position of affected genes in leukemia cells.

Animals↗

Distinct molecular phenotype of malignant CD34(+) hematopoietic stem and progenitor cells in chronic myelogenous leukemia.

Chronic myelogenous leukemia (CML) is a malignant disorder of the hematopoietic stem cell characterized by the BCR-ABL oncogene. We examined gene expression profiles of highly enriched CD34(+) hematopoietic stem and progenitor cells from patients with CML in chronic phase using cDNA arrays covering 1.185 genes. Comparing CML CD34(+) cells with normal CD34(+) cells, we found 158 genes which were significantly differentially expressed. Gene expression patterns reflected BCR-ABL-induced functional alterations such as increased cell-cycle and proteasome activity. Detoxification enzymes and DNA repair proteins were downregulated in CML CD34(+) cells, which might contribute to genetic instability. Decreased expression of junction plakoglobulin and CXC chemokine receptor 4 (CXCR-4) might facilitate the release of immature precursors from bone marrow in CML. GATA-2 was upregulated in CML CD34(+) cells, suggesting an increased self-renewal in comparison with normal CD34(+) cells. Moreover, we found upregulation of the proto-oncogene SKI and of receptors for neuromediators such as opioid mu1 receptor, GABA B receptor, adenosine A1 receptor, orexin 1 and 2 receptors and corticotropine-releasing hormone receptor. Treatment of CML progenitor cells with the selective adenosine A1 receptor antagonist 8-cyclopentyl-1,3-dipropylxanthine (DPCPX) resulted in a dose-dependent significant inhibition of clonogenic growth by 40% at a concentration of 10(-5) M, which could be reversed by the equimolar addition of the receptor agonist 2-chloro-N6-cyclopentyladenosine (P<0.05). The incubation of normal progenitor cells with DPCPX resulted in an inhibition of clonogenic growth to a significantly lesser extent in comparison with CML cells (P<0.05), suggesting that the adenosine A1 receptor is of functional relevance in CML hematopoietic progenitor cells.

Adult↗