BMC Medical Genomics - Most viewed articles

An integrative genomic approach reveals coordinated expression of intronic miR-335, miR-342, and miR-561 with deregulated host genes in multiple myeloma

2008-08-13

Background: The role of microRNAs (miRNAs) in multiple myeloma (MM) has yet to be fully elucidated. To identify miRNAs that are potentially deregulated in MM, we investigated those mapping within transcription units, based on evidence that intronic miRNAs are frequently coexpressed with their host genes. To this end, we monitored host transcript expression values in a panel of 20 human MM cell lines (HMCLs) and focused on transcripts whose expression varied significantly across the dataset. Methods: miRNA expression was quantified by Quantitative Real-Time PCR. Gene expression and genome profiling data were generated on Affymetrix oligonucleotide microarrays. Significant Analysis of Microarrays algorithm was used to investigate differentially expressed transcripts. Conventional statistics were used to test correlations for significance. Public libraries were queried to predict putative miRNA targets. Results: We identified transcripts specific to six miRNA host genes (CCPG1, GULP1, EVL, TACSTD1, MEST, and TNIK) whose average changes in expression varied at least 2-fold from the mean of the examined dataset. We evaluated the expression levels of the corresponding intronic miRNAs and identified a significant correlation between the expression levels of MEST, EVL, and GULP1 and those of the corresponding miRNAs miR-335, miR-342-3p, and miR-561, respectively. Genome-wide profiling of the 20 HMCLs indicated that the increased expression of the three host genes and their corresponding intronic miRNAs was not correlated with local copy number variations. Notably, miRNAs and their host genes were overexpressed in a fraction of primary tumors with respect to normal plasma cells; however, this finding was not correlated with known molecular myeloma groups. The predicted putative miRNA targets and the transcriptional profiles associated with the primary tumors suggest that MEST/miR-335 and EVL/miR-342-3p may play a role in plasma cell homing and/or interactions with the bone marrow microenvironment. Conclusion: Our data support the idea that intronic miRNAs and their host genes are regulated dependently, and may contribute to the understanding of their biological roles in cancer. To our knowledge, this is the first evidence of deregulated miRNA expression in MM, providing insights that may lead to the identification of new biomarkers and altered molecular pathways of the disease.

Role of Caveolin 1, E-Cadherin, Enolase 2 and PKCalpha on resistance to methotrexate in human HT29 colon cancer cells

2008-08-11

Background: Methotrexate is one of the earliest cytotoxic drugs used in cancer therapy, and despite the isolation of multiple other folate antagonists, methotrexate maintains its significant role as a treatment for different types of cancer and other disorders. The usefulness of treatment with methotrexate is limited by the development of drug resistance, which may be acquired through different ways. To get insights into the mechanisms associated with drug resistance and sensitization we performed a functional analysis of genes deregulated in methotrexate resistant cells, either due to its co-amplification with the dhfr gene or as a result of a transcriptome screening using microarrays. Methods: Gene expression levels were compared between triplicate samples from either HT29 sensitive cells and resistant to 10-5 M MTX by hybridization to the GeneChip® HG U133 PLUS 2.0 from Affymetrix. After normalization, a list of 3-fold differentially expressed genes with a p-value < 0.05 including multiple testing correction (Benjamini and Hochberg false discovery rate) was generated. RT-Real-time PCR was used to validate the expression levels of selected genes and copy-number was determined by qPCR. Functional validations were performed either by siRNAs or by transfection of an expression plasmid. Results: Genes adjacent to the dhfr locus and included in the 5q14 amplicon were overexpressed in HT29 MTX-resistant cells. Treatment with siRNAs against those genes caused a slight reduction in cell viability in both HT29 sensitive and resistant cells. On the other hand, microarray analysis of HT29 and HT29 MTX resistant cells unveiled overexpression of caveolin 1, enolase 2 and PKCα genes in resistant cells without concomitant copy number gain. siRNAs against these three genes effectively reduced cell viability and caused a decreased MTX resistance capacity. Moreover, overexpression of E-cadherin, which was found underexpressed in MTX-resistant cells, also sensitized the cells toward the chemotherapeutic agent. Combined treatments targeting siRNA inhibition of caveolin 1 and overexpression of E-cadherin markedly reduced cell viability in both sensitive and MTX-resistant HT29 cells. Conclusion: We provide functional evidences indicating that caveolin 1 and E-cadherin, deregulated in MTX resistant cells, may play a critical role in cell survival and may constitute potential targets for coadjuvant therapy.

Differential expression of genes mapping to recurrently abnormal chromosomal regions characterize neuroblastic tumours with distinct ploidy status

2008-08-13

Background: Neuroblastic tumours (NBTs) represent a heterogeneous spectrum of neoplastic diseases associated with multiple genetic alterations. Structural and numerical chromosomal changes are frequent and are predictive parameters of NBTs outcome. We performed a comparative analysis of the biological entities constituted by NBTs with different ploidy status. Methods: Gene expression profiling of 49 diagnostic primary NBTs with ploidy data was performed using oligonucleotide microarray. Further analyses using Quantitative Real-Time Polymerase Chain Reaction (Q-PCR); array-Comparative Genomic Hybridization (aCGH); and Fluorescent in situ Hybridization (FISH) were performed to investigate the correlation between aneuploidy, chromosomal changes and gene expression profiles. Results: Gene expression profiling of 49 primary near-triploid and near-diploid/tetraploid NBTs revealed distinct expression profiles associated with each NBT subgroup. A statistically significant portion of genes mapped to 1p36 (P = 0.01) and 17p13-q21 (P < 0.0001), described as recurrently altered in NBTs. Over 90% of these genes showed higher expression in near-triploid NBTs and the majority are involved in cell differentiation pathways. Specific chromosomal abnormalities observed in NBTs, 1p loss, 17q and whole chromosome 17 gains, were reflected in the gene expression profiles. Comparison between gene copy number and expression levels suggests that differential expression might be only partly dependent on gene copy number. Intratumoural clonal heterogeneity was observed in all NBTs, with marked interclonal variability in near-diploid/tetraploid tumours. Conclusion: NBTs with different cellular DNA content display distinct transcriptional profiles with a significant portion of differentially expressed genes mapping to specific chromosomal regions known to be associated with outcome. Furthermore, our results demonstrate that these specific genetic abnormalities are highly heterogeneous in all NBTs, and suggest that NBTs with different ploidy status may result from different mechanisms of aneuploidy driving tumourigenesis.

Tobacco use induces anti-apoptotic, proliferative patterns of gene expression in circulating leukocytes of Caucasian males

2008-08-18

Background: Strong epidemiologic evidence correlates tobacco use with a variety of serious adverse health effects, but the biological mechanisms that produce these effects remain elusive. Results: We analyzed gene transcription data to identify expression spectra related to tobacco use in circulating leukocytes of 67 Caucasian male subjects. Levels of cotinine, a nicotine metabolite, were used as a surrogate marker for tobacco exposure. Significance Analysis of Microarray and Gene Set Analysis identified 109 genes in 16 gene sets whose transcription levels were differentially regulated by nicotine exposure. We subsequently analyzed this gene set by hyperclustering, a technique that allows the data to be clustered by both expression ratio and gene annotation (e.g. Gene Ontologies). Conclusion: Our results demonstrate that tobacco use affects transcription of groups of genes that are involved in proliferation and apoptosis in circulating leukocytes. These transcriptional effects include a repertoire of transcriptional changes likely to increase the incidence of neoplasia through an altered expression of genes associated with transcription and signaling, interferon responses and repression of apoptotic pathways.

Influence of monolayer, spheroid, and tumor growth conditions on chromosome 3 gene expression in tumorigenic epithelial ovarian cancer cell lines

2008-08-07

Background: Expression microarray analyses of epithelial ovarian cancer (EOC) cell lines may be exploited to elucidate genetic and epigenetic events important in this disease. A possible variable is the influence of growth conditions on discerning candidates. The present study examined the influence of growth conditions on the expression of chromosome 3 genes in the tumorigenic EOC cell lines, OV-90, TOV-21G and TOV-112D using Affymetrix GeneChip® HG-U133A expression microarray analysis. Methods: Chromosome 3 gene expression profiles (n = 1147 probe sets, representing 735 genes) were extracted from U133A expression microarray analyses of the EOC cell lines OV-90, TOV-21G and TOV-112D that were grown as monolayers, spheroids or nude mouse xenografts and monolayers derived from these tumors. Hierarchical cluster analysis was performed to compare chromosome 3 transcriptome patterns of each growth condition. Differentially expressed genes were identified and characterized by two-way comparative analyses of fold-differences in gene expression between monolayer cultures and each of the other growth conditions, and between the maximum and minimum values of expression of all growth conditions for each EOC cell line. Results: An overall high degree of similarity (> 90%) in gene expression was observed when expression values of alternative growth conditions were compared within each EOC cell line group. Two-way comparative analysis of each EOC cell line grown in an alternative condition relative to the monolayer culture showed that overall less than 15% of probe sets exhibited at least a 3-fold difference in expression profile. Less than 23% of probe sets exhibited greater than 3-fold differences in gene expression in comparisons of the maximum and minimum value of expression of all growth conditions within each EOC cell line group. The majority of these differences were less than 5-fold. There were 17 genes in common which were differentially expressed in all EOC cell lines. However, the patterns of expression of these genes were not necessarily the same for each growth condition when one cell line was compared with another. Conclusion: The various alternative in vivo and in vitro growth conditions of tumorigenic EOC cell lines appeared to modestly influence the global chromosome 3 transcriptome supporting the notion that the in vitro cell line models are a viable option for testing gene candidates.

The gene expression profiles of primary and metastatic melanoma yields a transition point of tumor progression and metastasis

2008-04-28

Background: The process of malignant transformation, progression and metastasis of melanoma is poorly understood. Gene expression profiling of human cancer has allowed for a unique insight into the genes that are involved in these processes. Thus, we have attempted to utilize this approach through the analysis of a series of primary, non-metastatic cutaneous tumors and metastatic melanoma samples. Methods: We have utilized gene microarray analysis and a variety of molecular techniques to compare 40 metastatic melanoma (MM) samples, composed of 22 bulky, macroscopic (replaced) lymph node metastases, 16 subcutaneous and 2 distant metastases (adrenal and brain), to 42 primary cutaneous cancers, comprised of 16 melanoma, 11 squamous cell, 15 basal cell skin cancers. A Human Genome U133 Plus 2.0 array from Affymetrix, Inc. was utilized for each sample. A variety of statistical software, including the Affymetrix MAS 5.0 analysis software, was utilized to compare primary cancers to metastatic melanomas. Separate analyses were performed to directly compare only primary melanoma to metastatic melanoma samples. The expression levels of putative oncogenes and tumor suppressor genes were analyzed by semi- and real-time quantitative RT-PCR (qPCR) and Western blot analysis was performed on select genes. Results: We find that primary basal cell carcinomas, squamous cell carcinomas and thin melanomas express dramatically higher levels of many genes, including SPRR1A/B, KRT16/17, CD24, LOR, GATA3, MUC15, and TMPRSS4, than metastatic melanoma. In contrast, the metastatic melanomas express higher levels of genes such as MAGE, GPR19, BCL2A1, MMP14, SOX5, BUB1, RGS20, and more. The transition from non-metastatic expression levels to metastatic expression levels occurs as melanoma tumors thicken. We further evaluated primary melanomas of varying Breslow's tumor thickness to determine that the transition in expression occurs at different thicknesses for different genes suggesting that the "transition zone" represents a critical time for the emergence of the metastatic phenotype. Several putative tumor oncogenes (SPP-1, MITF, CITED-1, GDF-15, c-Met, HOX loci) and suppressor genes (PITX-1, CST-6, PDGFRL, DSC-3, POU2F3, CLCA2, ST7L), were identified and validated by quantitative PCR as changing expression during this transition period. These are strong candidates for genes involved in the progression or suppression of the metastatic phenotype. Conclusion: The gene expression profiling of primary, non-metastatic cutaneous tumors and metastatic melanoma has resulted in the identification of several genes that may be centrally involved in the progression and metastatic potential of melanoma. This has very important implications as we continue to develop an improved understanding of the metastatic process, allowing us to identify specific genes for prognostic markers and possibly for targeted therapeutic approaches.

Genome wide SNP comparative analysis between EGFR and KRAS mutated NSCLC and characterization of two models of oncogenic cooperation in non-small cell lung carcinoma

2008-06-12

Background: Lung cancer with EGFR mutation was shown to be a specific clinical entity. In order to better understand the biology behind this disease we used a genome wide characterization of loss of heterozygosity and amplification by Single Nucleotide Polymorphism (SNP) Array analysis to point out chromosome segments linked to EGFR mutations. To do so, we compared genetic profiles between EGFR mutated adenocarcinomas (ADC) and KRAS mutated ADC from 24 women with localized lung cancer. Results: Patterns of alterations were different between EGFR and KRAS mutated tumors and specific chromosomes alterations were linked to the EGFR mutated group. Indeed chromosome regions 14q21.3 (p = 0.027), 7p21.3-p21.2 (p = 0.032), 7p21.3 (p = 0.042) and 7p21.2-7p15.3 (p = 0.043) were found significantly amplified in EGFR mutated tumors. Within those regions 3 genes are of special interest ITGB8, HDAC9 and TWIST1. Moreover, homozygous deletions at CDKN2A and LOH at RB1 were identified in EGFR mutated tumors. We therefore tested the existence of a link between EGFR mutation, CDKN2A homozygous deletion and cyclin amplification in a larger series of tumors. Indeed, in a series of non-small-cell lung carcinoma (n = 98) we showed that homozygous deletions at CDKN2A were linked to EGFR mutations and absence of smoking whereas cyclin amplifications (CCNE1 and CCND1) were associated to TP53 mutations and smoking habit. Conclusion: All together, our results show that genome wide patterns of alteration differ between EGFR and KRAS mutated lung ADC, describe two models of oncogenic cooperation involving either EGFR mutation and CDKN2A deletion or cyclin amplification and TP53 inactivating mutations and identified new chromosome regions at 7p and 14q associated to EGFR mutations in lung cancer.

An approach to analyse the specific impact of rapamycin on mRNA-ribosome association.

2008-08-01

Background: Recent work, using both cell culture model systems and tumour derived cell lines, suggests that the differential recruitment into polysomes of mRNA populations may be sufficient to initiate and maintain tumour formation. Consequently, a major effort is underway to use high density microarray profiles to establish molecular fingerprints for cells exposed to defined drug regimes. The aim of these pharmacogenomic approaches is to provide new information on how drugs can impact on the translational read-out within a defined cellular background. Methods: We describe an approach that permits the analysis of de-novo mRNA-ribosome association in-vivo during short drug exposures. It combines hypertonic shock, polysome fractionation and high-throughput analysis to provide a molecular phenotype of translationally responsive transcripts. Compared to previous translational profiling studies, the procedure offers increased specificity due to the elimination of the drugs secondary effects (e.g. on the transcriptional read-out). For this pilot "proof-of-principle" assay we selected the drug rapamycin because of its extensively studied impact on translation initiation. Results: High throughput analysis on both the light and heavy polysomal fractions has identified mRNAs whose re-recruitment onto free ribosomes responded to short exposure to the drug rapamycin. The results of the microarray have been confirmed using real-time RT-PCR. The selective down-regulation of TOP transcripts is also consistent with previous translational profiling studies using this drug. Conclusions: The technical advance outlined in this manuscript offers the possibility of new insights into mRNA features that impact on translation initiation and provides a molecular fingerprint for transcript-ribosome association in any cell type and in the presence of a range of drugs of interest. Such molecular phenotypes defined pre-clinically may ultimately impact on the evaluation of a particular drug in a living cell.

Gene expression in periodontal tissues following treatment

2008-07-07

Background: In periodontitis, treatment aimed at controlling the periodontal biofilm infection results in a resolution of the clinical and histological signs of inflammation. Although the cell types found in periodontal tissues following treatment have been well described, information on gene expression is limited to few candidate genes. Therefore, the aim of the study was to determine the expression profiles of immune and inflammatory genes in periodontal tissues from sites with severe chronic periodontitis following periodontal therapy in order to identify genes involved in tissue homeostasis.Gingival biopsies from 12 patients with severe chronic periodontitis were taken six to eight weeks following non-surgical periodontal therapy, and from 11 healthy controls. As internal standard, RNA of an immortalized human keratinocyte line (HaCaT) was used. Total RNA was subjected to gene expression profiling using a commercially available microarray system focusing on inflammation-related genes. Post-hoc confirmation of selected genes was done by Realtime-PCR. Results: Out of the 136 genes analyzed, the 5% most strongly expressed genes compared to healthy controls were Interleukin-12A (IL-12A), Versican (CSPG-2), Matrixmetalloproteinase-1 (MMP-1), Down syndrome critical region protein-1 (DSCR-1), Macrophage inflammatory protein-2β (Cxcl-3), Inhibitor of apoptosis protein-1 (BIRC-1), Cluster of differentiation antigen 38 (CD38), Regulator of G-protein signalling-1 (RGS-1), and Finkel-Biskis-Jinkins murine osteosarcoma virus oncogene (C-FOS); the 5% least strongly expressed genes were Receptor-interacting Serine/Threonine Kinase-2 (RIP-2), Complement component 3 (C3), Prostaglandin-endoperoxide synthase-2 (COX-2), Interleukin-8 (IL-8), Endothelin-1 (EDN-1), Plasminogen activator inhibitor type-2 (PAI-2), Matrix-metalloproteinase-14 (MMP-14), and Interferon regulating factor-7 (IRF-7). Conclusion: Gene expression profiles found in periodontal tissues following therapy indicate activation of pathways that regulate tissue damage and repair.

My sister's keeper?: genomic research and the identifiability of siblings

Christopher A Cassa, Brian Schmidt, Isaac S Kohane and Kenneth D Mandl — 2008-07-25

Background: Genomic sequencing of SNPs is increasingly prevalent, though the amount of familial information these data contain has not been quantified. Methods: We provide a framework for measuring the risk to siblings of a patient's SNP genotype disclosure, and demonstrate that sibling SNP genotypes can be inferred with substantial accuracy. Results: Extending this inference technique, we determine that a very low number of matches at commonly varying SNPs is sufficient to confirm sib-ship, demonstrating that published sequence data can reliably be used to derive sibling identities. Using HapMap trio data, at SNPs where one child is homozygotic major, with a minor allele frequency ≤ 0.20, (N = 452684, 65.1%) we achieve 91.9% inference accuracy for sibling genotypes. Conclusion: These findings demonstrate that substantial discrimination and privacy risks arise from use of inferred familial genomic data.