miR-142-3p prevents macrophage differentiation during cancer-induced myelopoiesis.

Tumor progression is accompanied by an altered myelopoiesis causing the accumulation of immunosuppressive cells. Here, we showed that miR-142-3p downregulation promoted macrophage differentiation and determined the acquisition of their immunosuppressive function in tumor. Tumor-released cytokines signaling through gp130, the common subunit of the interleukin-6 cytokine receptor family, induced the LAP∗ isoform of C/EBPß transcription factor, promoting macrophage generation. miR-142-3p downregulated gp130 by canonical binding to its messenger RNA (mRNA) 3' UTR and repressed C/EBPß LAP∗ by noncanonical binding to its 5' mRNA coding sequence. Enforced miR expression impaired macrophage differentiation both in vitro and in vivo. Mice constitutively expressing miR-142-3p in the bone marrow showed a marked increase in survival following immunotherapy with tumor-specific T lymphocytes. By modulating a specific miR in bone marrow precursors, we thus demonstrated the feasibility of altering tumor-induced macrophage differentiation as a potent tool to improve the efficacy of cancer immunotherapy.

MiR&moRe2: A Bioinformatics Tool to Characterize microRNAs and microRNA-Offset RNAs from Small RNA-Seq Data.

MicroRNA-offset RNAs (moRNAs) are microRNA-like small RNAs generated by microRNA precursors. To date, little is known about moRNAs and bioinformatics tools to inspect their expression are still missing. We developed miR&moRe2, the first bioinformatics method to consistently characterize microRNAs, moRNAs, and their isoforms from small RNA sequencing data. To illustrate miR&moRe2 discovery power, we applied it to several published datasets. MoRNAs identified by miR&moRe2 were in agreement with previous research findings. Moreover, we observed that moRNAs and new microRNAs predicted by miR&moRe2 were downregulated upon the silencing of the microRNA-biogenesis pathway. Further, in a sizeable dataset of human blood cell populations, tens of novel miRNAs and moRNAs were discovered, some of them with significantly varied expression levels among the cell types. Results demonstrate that miR&moRe2 is a valid tool for a comprehensive study of small RNAs generated from microRNA precursors and could help to investigate their biogenesis and function.

miRNA-mRNA integrative analysis in primary myelofibrosis CD34+ cells: role of miR-155/JARID2 axis in abnormal megakaryopoiesis.

Primary myelofibrosis (PMF) is a myeloproliferative neoplasm characterized by megakaryocyte (MK) hyperplasia, bone marrow fibrosis, and abnormal stem cell trafficking. PMF may be associated with somatic mutations in JAK2, MPL, or CALR. Previous studies have shown that abnormal MKs play a central role in the pathophysiology of PMF. In this work, we studied both gene and microRNA (miRNA) expression profiles in CD34(+) cells from PMF patients. We identified several biomarkers and putative molecular targets such as FGR, LCN2, and OLFM4. By means of miRNA-gene expression integrative analysis, we found different regulatory networks involved in the dysregulation of transcriptional control and chromatin remodeling. In particular, we identified a network gathering several miRNAs with oncogenic potential (eg, miR-155-5p) and targeted genes whose abnormal function has been previously associated with myeloid neoplasms, including JARID2, NR4A3, CDC42, and HMGB3. Because the validation of miRNA-target interactions unveiled JARID2/miR-155-5p as the strongest relationship in the network, we studied the function of this axis in normal and PMF CD34(+) cells. We showed that JARID2 downregulation mediated by miR-155-5p overexpression leads to increased in vitro formation of CD41(+) MK precursors. These findings suggest that overexpression of miR-155-5p and the resulting downregulation of JARID2 may contribute to MK hyperplasia in PMF.

Characterization and discovery of novel miRNAs and moRNAs in JAK2V617F-mutated SET2 cells.

To gain insights into a possible role of microRNAs in myeloproliferative neoplasms, we performed short RNA massive sequencing and extensive bioinformatic analysis in the JAK2V617F-mutated SET2 cell line. Overall, 652 known mature miRNAs were detected, of which 21 were highly expressed, thus being responsible of most of miRNA-mediated gene repression. microRNA putative targets were enriched in specific signaling pathways, providing information about cell activities under massive posttranscriptional regulation. The majority of miRNAs were mixtures of sequence variants, called isomiRs, mainly because of alternative, noncanonical processing of hairpin precursors. We also identified 78 novel miRNAs (miRNA*) derived from known hairpin precursors. Both major and minor (*) forms of miRNAs were expressed concurrently from half of expressed hairpins, highlighting the relevance of miRNA* and the complexity of strand selection bias regulation. Finally, we discovered that SET2 cells express a number of miRNA-offset RNAs (moRNAs), short RNAs derived from genomic regions flanking mature miRNAs. We provide novel data about the possible origin of moRNAs, although their functional role remains to be elucidated. Overall, this study shed light on the complexity of microRNA-mediated gene regulation in SET2 cells and represents the basis for future studies in JAK2V617F-mutated cellular models.

miRNome of Italian Large White pig subcutaneous fat tissue: new miRNAs, isomiRs and moRNAs.

Small RNAs, such as micro-RNAs (miRNAs), are decisive regulators of gene expression, and they could determine adipose tissue traits. A better knowledge of porcine fat genomics is relevant given that the pig is a biomedical model for metabolic and cardiovascular human pathologies. Adipose tissue is particularly important for the meat industry. We explored the miRNome of two adult Italian Large White pig backfat samples by Illumina RNA-Seq. Using custom bioinformatic methods, the expressed miRNAs were identified and quantified and the nucleotide sequence variability of miRNA isoforms were analysed. We detected 222 known miRNAs, 68 new miRNAs and 17 miRNA-offset RNAs (moRNAs) expressed from known hairpins, and 312 new miRNAs expressed from 253 new hairpins. Porcine transcripts targeted by the most expressed miRNAs were predicted, showing that these miRNAs may have an impact on Wnt, insulin signalling and axon guidance pathways. The expression of five small RNAs, including moRNA ssc-5'-moR-21 and a miRNA from a new hairpin, was validated by a qRT-PCR assay, thus confirming the robustness of our results. The depicted miRNome complexity suggests that quantitative and qualitative features of miRNAs and non-canonical products of their precursors are worthy of further investigation to clarify their roles in the adipose tissue biology.

MAGIA²: from miRNA and genes expression data integrative analysis to microRNA-transcription factor mixed regulatory circuits (2012 update).

MAGIA(2) (http://gencomp.bio.unipd.it/magia2) is an update, extension and evolution of the MAGIA web tool. It is dedicated to the integrated analysis of in silico target prediction, microRNA (miRNA) and gene expression data for the reconstruction of post-transcriptional regulatory networks. miRNAs are fundamental post-transcriptional regulators of several key biological and pathological processes. As miRNAs act prevalently through target degradation, their expression profiles are expected to be inversely correlated to those of the target genes. Low specificity of target prediction algorithms makes integration approaches an interesting solution for target prediction refinement. MAGIA(2) performs this integrative approach supporting different association measures, multiple organisms and almost all target predictions algorithms. Nevertheless, miRNAs activity should be viewed as part of a more complex scenario where regulatory elements and their interactors generate a highly connected network and where gene expression profiles are the result of different levels of regulation. The updated MAGIA(2) tries to dissect this complexity by reconstructing mixed regulatory circuits involving either miRNA or transcription factor (TF) as regulators. Two types of circuits are identified: (i) a TF that regulates both a miRNA and its target and (ii) a miRNA that regulates both a TF and its target.

Impact of microRNAs on regulatory networks and pathways in human colorectal carcinogenesis and development of metastasis.

BACKGROUND: Qualitative alterations or abnormal expression of microRNAs (miRNAs) in colon cancer have mainly been demonstrated in primary tumors. Poorly overlapping sets of oncomiRs, tumor suppressor miRNAs and metastamiRs have been linked with distinct stages in the progression of colorectal cancer. To identify changes in both miRNA and gene expression levels among normal colon mucosa, primary tumor and liver metastasis samples, and to classify miRNAs into functional networks, in this work miRNA and gene expression profiles in 158 samples from 46 patients were analysed. RESULTS: Most changes in miRNA and gene expression levels had already manifested in the primary tumors while these levels were almost stably maintained in the subsequent primary tumor-to-metastasis transition. In addition, comparing normal tissue, tumor and metastasis, we did not observe general impairment or any rise in miRNA biogenesis. While only few mRNAs were found to be differentially expressed between primary colorectal carcinoma and liver metastases, miRNA expression profiles can classify primary tumors and metastases well, including differential expression of miR-10b, miR-210 and miR-708. Of 82 miRNAs that were modulated during tumor progression, 22 were involved in EMT. qRT-PCR confirmed the down-regulation of miR-150 and miR-10b in both primary tumor and metastasis compared to normal mucosa and of miR-146a in metastases compared to primary tumor. The upregulation of miR-201 in metastasis compared both with normal and primary tumour was also confirmed. A preliminary survival analysis considering differentially expressed miRNAs suggested a possible link between miR-10b expression in metastasis and patient survival. By integrating miRNA and target gene expression data, we identified a combination of interconnected miRNAs, which are organized into sub-networks, including several regulatory relationships with differentially expressed genes. Key regulatory interactions were validated experimentally. Specific mixed circuits involving miRNAs and transcription factors were identified and deserve further investigation. The suppressor activity of miR-182 on ENTPD5 gene was identified for the first time and confirmed in an independent set of samples. CONCLUSIONS: Using a large dataset of CRC miRNA and gene expression profiles, we describe the interplay of miRNA groups in regulating gene expression, which in turn affects modulated pathways that are important for tumor development.

Improving biomarker list stability by integration of biological knowledge in the learning process.

BACKGROUND: The identification of robust lists of molecular biomarkers related to a disease is a fundamental step for early diagnosis and treatment. However, methodologies for biomarker discovery using microarray data often provide results with limited overlap. It has been suggested that one reason for these inconsistencies may be that in complex diseases, such as cancer, multiple genes belonging to one or more physiological pathways are associated with the outcomes. Thus, a possible approach to improve list stability is to integrate biological information from genomic databases in the learning process; however, a comprehensive assessment based on different types of biological information is still lacking in the literature. In this work we have compared the effect of using different biological information in the learning process like functional annotations, protein-protein interactions and expression correlation among genes. RESULTS: Biological knowledge has been codified by means of gene similarity matrices and expression data linearly transformed in such a way that the more similar two features are, the more closely they are mapped. Two semantic similarity matrices, based on Biological Process and Molecular Function Gene Ontology annotation, and geodesic distance applied on protein-protein interaction networks, are the best performers in improving list stability maintaining almost equal prediction accuracy. CONCLUSIONS: The performed analysis supports the idea that when some features are strongly correlated to each other, for example because are close in the protein-protein interaction network, then they might have similar importance and are equally relevant for the task at hand. Obtained results can be a starting point for additional experiments on combining similarity matrices in order to obtain even more stable lists of biomarkers. The implementation of the classification algorithm is available at the link: http://www.math.unipd.it/~dasan/biomarkers.html.

MAGIA, a web-based tool for miRNA and Genes Integrated Analysis.

MAGIA (miRNA and genes integrated analysis) is a novel web tool for the integrative analysis of target predictions, miRNA and gene expression data. MAGIA is divided into two parts: the query section allows the user to retrieve and browse updated miRNA target predictions computed with a number of different algorithms (PITA, miRanda and Target Scan) and Boolean combinations thereof. The analysis section comprises a multistep procedure for (i) direct integration through different functional measures (parametric and non-parametric correlation indexes, a variational Bayesian model, mutual information and a meta-analysis approach based on P-value combination) of mRNA and miRNA expression data, (ii) construction of bipartite regulatory network of the best miRNA and mRNA putative interactions and (iii) retrieval of information available in several public databases of genes, miRNAs and diseases and via scientific literature text-mining. MAGIA is freely available for Academic users at http://gencomp.bio.unipd.it/magia.

Motif discovery in promoters of genes co-localized and co-expressed during myeloid cells differentiation.

Genes co-expressed may be under similar promoter-based and/or position-based regulation. Although data on expression, position and function of human genes are available, their true integration still represents a challenge for computational biology, hampering the identification of regulatory mechanisms. We carried out an integrative analysis of genomic position, functional annotation and promoters of genes expressed in myeloid cells. Promoter analysis was conducted by a novel multi-step method for discovering putative regulatory elements, i.e. over-represented motifs, in a selected set of promoters, as compared with a background model. The combination of transcriptional, structural and functional data allowed the identification of sets of promoters pertaining to groups of genes co-expressed and co-localized in regions of the human genome. The application of motif discovery to 26 groups of genes co-expressed in myeloid cells differentiation and co-localized in the genome showed that there are more over-represented motifs in promoters of co-expressed and co-localized genes than in promoters of simply co-expressed genes (CEG). Motifs, which are similar to the binding sequences of known transcription factors, non-uniformly distributed along promoter sequences and/or occurring in highly co-expressed subset of genes were identified. Co-expressed and co-localized gene sets were grouped in two co-expressed genomic meta-regions, putatively representing functional domains of a high-level expression regulation.

A-MADMAN: annotation-based microarray data meta-analysis tool.

BACKGROUND: Publicly available datasets of microarray gene expression signals represent an unprecedented opportunity for extracting genomic relevant information and validating biological hypotheses. However, the exploitation of this exceptionally rich mine of information is still hampered by the lack of appropriate computational tools, able to overcome the critical issues raised by meta-analysis. RESULTS: This work presents A-MADMAN, an open source web application which allows the retrieval, annotation, organization and meta-analysis of gene expression datasets obtained from Gene Expression Omnibus. A-MADMAN addresses and resolves several open issues in the meta-analysis of gene expression data. CONCLUSION: A-MADMAN allows i) the batch retrieval from Gene Expression Omnibus and the local organization of raw data files and of any related meta-information, ii) the re-annotation of samples to fix incomplete, or otherwise inadequate, metadata and to create user-defined batches of data, iii) the integrative analysis of data obtained from different Affymetrix platforms through custom chip definition files and meta-normalization. Software and documentation are available on-line at http://compgen.bio.unipd.it/bioinfo/amadman/.

Expanding the repertoire of miRNAs and miRNA-offset RNAs expressed in multiple myeloma by small RNA deep sequencing.

Microarray analysis of the multiple myeloma (MM) miRNome has unraveled the differential expression of miRNAs in cytogenetic subgroups, their involvement in the tumor biology and their effectiveness in prognostic models. Herein, the small RNA transcriptional landscape in MM has been investigated exploiting the possibilities offered by small RNA-seq, including accurate quantification of known mature species, discovery and characterization of isomiRs, and miRNA-offset RNAs (moRNAs). Matched small RNA-seq and miRNA GeneChip® microarray expression profiles were obtained in a representative panel of 30 primary MM tumors, fully characterized for genomic aberrations and mutations. RNA-seq and microarray gave concordant estimations of known species. Enhanced analysis of RNA-seq data with the miR&moRe pipeline led to the characterization of 655 known and 17 new mature miRNAs and of 74 moRNAs expressed in the considered cohort, 5 of which (moR-150-3p, moR-24-2-5p, moR-421-5p, moR-21-5p, and moR-6724-5p) at high level. Ectopic expression of miR-135a-3p in t(4;14) patients, upregulation of moR-150-3p and moR-21-5p in t(14;16)/t(14;20) samples, and of moR-6724-1-5p in patients overexpressing CCND1 were uncovered and validated by qRT-PCR. Overall, RNA-seq offered a more complete overview of small non-coding RNA in MM tumors, indicating specific moRNAs that demand further investigations to explore their role in MM biology.

Detecting seeded motifs in DNA sequences.

The problem of detecting DNA motifs with functional relevance in real biological sequences is difficult due to a number of biological, statistical and computational issues and also because of the lack of knowledge about the structure of searched patterns. Many algorithms are implemented in fully automated processes, which are often based upon a guess of input parameters from the user at the very first step. In this paper, we present a novel method for the detection of seeded DNA motifs, composed by regions with a different extent of variability. The method is based on a multi-step approach, which was implemented in a motif searching web tool (MOST). Overrepresented exact patterns are extracted from input sequences and clustered to produce motifs core regions, which are then extended and scored to generate seeded motifs. The combination of automated pattern discovery algorithms and different display tools for the evaluation and selection of results at several analysis steps can potentially lead to much more meaningful results than complete automation can produce. Experimental results on different yeast and human real datasets proved the methodology to be a promising solution for finding seeded motifs. MOST web tool is freely available at http://telethon.bio.unipd.it/bioinfo/MOST.

Disentangling the microRNA regulatory milieu in multiple myeloma: integrative genomics analysis outlines mixed miRNA-TF circuits and pathway-derived networks modulated in t(4;14) patients.

The identification of overexpressed miRNAs in multiple myeloma (MM) has progressively added a further level of complexity to MM biology. miRNA and gene expression profiles of two large representative MM datasets, available from retrospective and prospective series and encompassing a total of 249 patients at diagnosis, were analyzed by means of in silico integrative genomics methods, based on MAGIA2 and Micrographite computational procedures. We first identified relevant miRNA/transcription factors/target gene regulation circuits in the disease and linked them to biological processes. Members of the miR-99b/let-7e/miR-125a cluster, or of its paralog, upregulated in t(4;14), were connected with the specific transcription factors PBX1 and CEBPA and several target genes. These results were validated in two additional independent plasma cell tumor datasets. Then, we reconstructed a non-redundant miRNA-gene regulatory network in MM, linking miRNAs, such as let-7g, miR-19a, mirR-20a, mir-21, miR-29 family, miR-34 family, miR-125b, miR-155, miR-221 to pathways associated with MM subtypes, in particular the ErbB, the Hippo, and the Acute myeloid leukemia associated pathways.

A multistep bioinformatic approach detects putative regulatory elements in gene promoters.

BACKGROUND: Searching for approximate patterns in large promoter sequences frequently produces an exceedingly high numbers of results. Our aim was to exploit biological knowledge for definition of a sheltered search space and of appropriate search parameters, in order to develop a method for identification of a tractable number of sequence motifs. RESULTS: Novel software (COOP) was developed for extraction of sequence motifs, based on clustering of exact or approximate patterns according to the frequency of their overlapping occurrences. Genomic sequences of 1 Kb upstream of 91 genes differentially expressed and/or encoding proteins with relevant function in adult human retina were analyzed. Methodology and results were tested by analysing 1,000 groups of putatively unrelated sequences, randomly selected among 17,156 human gene promoters. When applied to a sample of human promoters, the method identified 279 putative motifs frequently occurring in retina promoters sequences. Most of them are localized in the proximal portion of promoters, less variable in central region than in lateral regions and similar to known regulatory sequences. COOP software and reference manual are freely available upon request to the Authors. CONCLUSION: The approach described in this paper seems effective for identifying a tractable number of sequence motifs with putative regulatory role.

Novel genes, possibly relevant for molecular diagnosis or therapy of human rhabdomyosarcoma, detected by genomic expression profiling.

Transcriptional profiles of an alveolar rhabdomyosarcoma (RMS) and of a RMS cell line were reconstructed by a computational and statistical approach. Expression data of 29,963 genes in 11 adult human healthy tissues and in 37 tumour tissues were analysed for comparison. We identified 202 genes differentially expressed in at least one RMS sample, as compared with normal skeletal muscle. Among them, 107 resulted specifically overexpressed in RMS, but in no tumour affecting other tissues. Cluster analysis applied to expression data detected a series of genes presumably co-expressed with genes encoding known tumour markers and/or reportedly involved in genesis or development of rhabdomyosarcoma. This study succeeded in identifying a number of genes, which become candidates for in vitro study, thus facilitating discovery of novel tumour markers or targets for drug therapy.

Small RNA Sequencing Uncovers New miRNAs and moRNAs Differentially Expressed in Normal and Primary Myelofibrosis CD34+ Cells.

Myeloproliferative neoplasms (MPN) are chronic myeloid cancers thought to arise at the level of CD34+ hematopoietic stem/progenitor cells. They include essential thrombocythemia (ET), polycythemia vera (PV) and primary myelofibrosis (PMF). All can progress to acute leukemia, but PMF carries the worst prognosis. Increasing evidences indicate that deregulation of microRNAs (miRNAs) might plays an important role in hematologic malignancies, including MPN. To attain deeper knowledge of short RNAs (sRNAs) expression pattern in CD34+ cells and of their possible role in mediating post-transcriptional regulation in PMF, we sequenced with Illumina HiSeq2000 technology CD34+ cells from healthy subjects and PMF patients. We detected the expression of 784 known miRNAs, with a prevalence of miRNA up-regulation in PMF samples, and discovered 34 new miRNAs and 99 new miRNA-offset RNAs (moRNAs), in CD34+ cells. Thirty-seven small RNAs were differentially expressed in PMF patients compared with healthy subjects, according to microRNA sequencing data. Five miRNAs (miR-10b-5p, miR-19b-3p, miR-29a-3p, miR-379-5p, and miR-543) were deregulated also in PMF granulocytes. Moreover, 3'-moR-128-2 resulted consistently downregulated in PMF according to RNA-seq and qRT-PCR data both in CD34+ cells and granulocytes. Target predictions of these validated small RNAs de-regulated in PMF and functional enrichment analyses highlighted many interesting pathways involved in tumor development and progression, such as signaling by FGFR and DAP12 and Oncogene Induced Senescence. As a whole, data obtained in this study deepened the knowledge of miRNAs and moRNAs altered expression in PMF CD34+ cells and allowed to identify and validate a specific small RNA profile that distinguishes PMF granulocytes from those of normal subjects. We thus provided new information regarding the possible role of miRNAs and, specifically, of new moRNAs in this disease.

H-ferritin-regulated microRNAs modulate gene expression in K562 cells.

In a previous study, we showed that the silencing of the heavy subunit (FHC) offerritin, the central iron storage molecule in the cell, is accompanied by a modification in global gene expression. In this work, we explored whether different FHC amounts might modulate miRNA expression levels in K562 cells and studied the impact of miRNAs in gene expression profile modifications. To this aim, we performed a miRNA-mRNA integrative analysis in K562 silenced for FHC (K562shFHC) comparing it with K562 transduced with scrambled RNA (K562shRNA). Four miRNAs, namely hsa-let-7g, hsa-let-7f, hsa-let-7i and hsa-miR-125b, were significantly up-regulated in silenced cells. The remarkable down-regulation of these miRNAs, following FHC expression rescue, supports a specific relation between FHC silencing and miRNA-modulation. The integration of target predictions with miRNA and gene expression profiles led to the identification of a regulatory network which includes the miRNAs up-regulated by FHC silencing, as well as91 down-regulated putative target genes. These genes were further classified in 9 networks; the highest scoring network, "Cell Death and Survival, Hematological System Development and Function, Hematopoiesis", is composed by 18 focus molecules including RAF1 and ERK1/2. We confirmed that, following FHC silencing, ERK1/2 phosphorylation is severely impaired and that RAF1 mRNA is significantly down-regulated. Taken all together, our data indicate that, in our experimental model, FHC silencing may affect RAF1/pERK1/2 levels through the modulation of a specific set of miRNAs and add new insights in to the relationship among iron homeostasis and miRNAs.

A constitutive active MAPK/ERK pathway due to BRAFV600E positively regulates AHR pathway in PTC.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor mediating the toxicity and tumor-promoting properties of dioxin. AHR has been reported to be overexpressed and constitutively active in a variety of solid tumors, but few data are currently available concerning its role in thyroid cancer. In this study we quantitatively explored a series of 51 paired-normal and papillary thyroid carcinoma (PTC) tissues for AHR-related genes. We identified an increased AHR expression/activity in PTC, independently from its nuclear dimerization partner and repressor but strictly related to a constitutive active MAPK/ERK pathway. The AHR up-regulation followed by an increased expression of AHR target genes was confirmed by a meta-analysis of published microarray data, suggesting a ligand-independent active AHR pathway in PTC. In-vitro studies using a PTC-derived cell line (BCPAP) and HEK293 cells showed that BRAFV600E may directly modulate AHR localization, induce AHR expression and activity in an exogenous ligand-independent manner. The AHR pathway might represent a potential novel therapeutic target for PTC in the clinical practice.

Detection of chromosomal regions showing differential gene expression in human skeletal muscle and in alveolar rhabdomyosarcoma.

BACKGROUND: Rhabdomyosarcoma is a relatively common tumour of the soft tissue, probably due to regulatory disruption of growth and differentiation of skeletal muscle stem cells. Identification of genes differentially expressed in normal skeletal muscle and in rhabdomyosarcoma may help in understanding mechanisms of tumour development, in discovering diagnostic and prognostic markers and in identifying novel targets for drug therapy. RESULTS: A Perl-code web client was developed to automatically obtain genome map positions of large sets of genes. The software, based on automatic search on Human Genome Browser by sequence alignment, only requires availability of a single transcribed sequence for each gene. In this way, we obtained tissue-specific chromosomal maps of genes expressed in rhabdomyosarcoma or skeletal muscle. Subsequently, Perl software was developed to calculate gene density along chromosomes, by using a sliding window. Thirty-three chromosomal regions harbouring genes mostly expressed in rhabdomyosarcoma were identified. Similarly, 48 chromosomal regions were detected including genes possibly related to function of differentiated skeletal muscle, but silenced in rhabdomyosarcoma. CONCLUSION: In this study we developed a method and the associated software for the comparative analysis of genomic expression in tissues and we identified chromosomal segments showing differential gene expression in human skeletal muscle and in alveolar rhabdomyosarcoma, appearing as candidate regions for harbouring genes involved in origin of alveolar rhabdomyosarcoma representing possible targets for drug treatment and/or development of tumor markers.