A statistical method for measuring activation of gene regulatory networks.

MOTIVATION: Gene expression data analysis is of great importance for modern molecular biology, given our ability to measure the expression profiles of thousands of genes and enabling studies rooted in systems biology. In this work, we propose a simple statistical model for the activation measuring of gene regulatory networks, instead of the traditional gene co-expression networks. RESULTS: We present the mathematical construction of a statistical procedure for testing hypothesis regarding gene regulatory network activation. The real probability distribution for the test statistic is evaluated by a permutation based study. To illustrate the functionality of the proposed methodology, we also present a simple example based on a small hypothetical network and the activation measuring of two KEGG networks, both based on gene expression data collected from gastric and esophageal samples. The two KEGG networks were also analyzed for a public database, available through NCBI-GEO, presented as Supplementary Material. AVAILABILITY: This method was implemented in an R package that is available at the BioConductor project website under the name maigesPack.

Gene network analyses point to the importance of human tissue kallikreins in melanoma progression.

BACKGROUND: A wide variety of high-throughput microarray platforms have been used to identify molecular targets associated with biological and clinical tumor phenotypes by comparing samples representing distinct pathological states. METHODS: The gene expression profiles of human cutaneous melanomas were determined by cDNA microarray analysis. Next, a robust analysis to determine functional classifications and make predictions based on data-oriented hypotheses was performed. Relevant networks that may be implicated in melanoma progression were also considered. RESULTS: In this study we aimed to analyze coordinated gene expression changes to find molecular pathways involved in melanoma progression. To achieve this goal, ontologically-linked modules with coordinated expression changes in melanoma samples were identified. With this approach, we detected several gene networks related to different modules that were induced or repressed during melanoma progression. Among them we observed high coordinated expression levels of genes involved in a) cell communication (KRT4, VWF and COMP); b) epidermal development (KLK7, LAMA3 and EVPL); and c) functionally related to kallikreins (EVPL, KLK6, KLK7, KLK8, SERPINB13, SERPING1 and SLPI). Our data also indicated that hKLK7 protein expression was significantly associated with good prognosis and survival. CONCLUSIONS: Our findings, derived from a different type of analysis of microarray data, highlight the importance of analyzing coordinated gene expression to find molecular pathways involved in melanoma progression.

Early changes in gene expression induced by tobacco smoke: Evidence for the importance of estrogen within lung tissue.

Lung cancer is the leading cause of cancer deaths in the United States, surpassing breast cancer as the primary cause of cancer-related mortality in women. The goal of the present study was to identify early molecular changes in the lung induced by exposure to tobacco smoke and thus identify potential targets for chemoprevention. Female A/J mice were exposed to either tobacco smoke or HEPA-filtered air via a whole-body exposure chamber (6 h/d, 5 d/wk for 3, 8, and 20 weeks). Gene expression profiles of lung tissue from control and smoke-exposed animals were established using a 15K cDNA microarray. Cytochrome P450 1b1, a phase I enzyme involved in both the metabolism of xenobiotics and the 4-hydroxylation of 17beta-estradiol (E(2)), was modulated to the greatest extent following smoke exposure. A panel of 10 genes were found to be differentially expressed in control and smoke-exposed lung tissues at 3, 8, and 20 weeks (P < 0.001). The interaction network of these differentially expressed genes revealed new pathways modulated by short-term smoke exposure, including estrogen metabolism. In addition, E(2) was detected within murine lung tissue by gas chromatography-coupled mass spectrometry and immunohistochemistry. Identification of the early molecular events that contribute to lung tumor formation is anticipated to lead to the development of promising targeted chemopreventive therapies. In conclusion, the presence of E(2) within lung tissue when combined with the modulation of cytochrome P450 1b1 and other estrogen metabolism genes by tobacco smoke provides novel insight into a possible role for estrogens in lung cancer.

Characterization of global transcription profile of normal and HPV-immortalized keratinocytes and their response to TNF treatment.

BACKGROUND: Persistent infection by high risk HPV types (e.g. HPV-16, -18, -31, and -45) is the main risk factor for development of cervical intraepithelial neoplasia and cervical cancer. Tumor necrosis factor (TNF) is a key mediator of epithelial cell inflammatory response and exerts a potent cytostatic effect on normal or HPV16, but not on HPV18 immortalized keratinocytes. Moreover, several cervical carcinoma-derived cell lines are resistant to TNF anti-proliferative effect suggesting that the acquisition of TNF-resistance may constitute an important step in HPV-mediated carcinogenesis. In the present study, we compared the gene expression profiles of normal and HPV16 or 18 immortalized human keratinocytes before and after treatment with TNF for 3 or 60 hours. METHODS: In this study, we determined the transcriptional changes 3 and 60 hours after TNF treatment of normal, HPV16 and HPV18 immortalized keratinocytes by microarray analysis. The expression pattern of two genes observed by microarray was confirmed by Northern Blot. NF-kappaB activation was also determined by electrophoretic mobility shift assay (EMSA) using specific oligonucleotides and nuclear protein extracts. RESULTS: We observed the differential expression of a common set of genes in two TNF-sensitive cell lines that differs from those modulated in TNF-resistant ones. This information was used to define genes whose differential expression could be associated with the differential response to TNF, such as: KLK7 (kallikrein 7), SOD2 (superoxide dismutase 2), 100P (S100 calcium binding protein P), PI3 (protease inhibitor 3, skin-derived), CSTA (cystatin A), RARRES1 (retinoic acid receptor responder 1), and LXN (latexin). The differential expression of the KLK7 and SOD2 transcripts was confirmed by Northern blot. Moreover, we observed that SOD2 expression correlates with the differential NF-kappaB activation exhibited by TNF-sensitive and TNF-resistant cells. CONCLUSION: This is the first in depth analysis of the differential effect of TNF on normal and HPV16 or HPV18 immortalized keratinocytes. Our findings may be useful for the identification of genes involved in TNF resistance acquisition and candidate genes which deregulated expression may be associated with cervical disease establishment and/or progression.

New insights about host response to smallpox using microarray data.

BACKGROUND: Smallpox is a lethal disease that was endemic in many parts of the world until eradicated by massive immunization. Due to its lethality, there are serious concerns about its use as a bioweapon. Here we analyze publicly available microarray data to further understand survival of smallpox infected macaques, using systems biology approaches. Our goal is to improve the knowledge about the progression of this disease. RESULTS: We used KEGG pathways annotations to define groups of genes (or modules), and subsequently compared them to macaque survival times. This technique provided additional insights about the host response to this disease, such as increased expression of the cytokines and ECM receptors in the individuals with higher survival times. These results could indicate that these gene groups could influence an effective response from the host to smallpox. CONCLUSION: Macaques with higher survival times clearly express some specific pathways previously unidentified using regular gene-by-gene approaches. Our work also shows how third party analysis of public datasets can be important to support new hypotheses to relevant biological problems.

Expression profile of malignant and nonmalignant lesions of esophagus and stomach: differential activity of functional modules related to inflammation and lipid metabolism.

Adenocarcinomas of stomach and esophagus are frequently associated with preceding inflammatory alterations of the normal mucosa. Whereas intestinal metaplasia of the gastric mucosa is associated with higher risk of malignization, Barrett's disease is a risk factor for adenocarcinoma of the esophagus. Barrett's disease is characterized by the substitution of the squamous mucosa of the esophagus by a columnar tissue classified histopathologically as intestinal metaplasia. Using cDNA microarrays, we determined the expression profile of normal gastric and esophageal mucosa as well as intestinal metaplasia and adenocarcinomas from both organs. Data were explored to define functional alterations related to the transformation from squamous to columnar epithelium and the malignant transformation from intestinal metaplasia to adenocarcinomas. Based on their expression profile, adenocarcinomas of the esophagus showed stronger correlation with intestinal metaplasia of the stomach than with Barrett's mucosa. Second, we identified two functional modules, lipid metabolism and cytokine, as being altered with higher statistical significance. Whereas the lipid metabolism module is active in samples representing intestinal metaplasia and inactive in adenocarcinomas, the cytokine module is inactive in samples representing normal esophagus and esophagitis. Using the concept of relevance networks, we determined the changes in linear correlation of genes pertaining to these two functional modules. Exploitation of the data presented herein will help in the precise molecular characterization of adenocarcinoma from the distal esophagus, avoiding the topographical and descriptive classification that is currently adopted, and help with the proper management of patients with Barrett's disease.