Analysis of spatial-temporal gene expression patterns reveals dynamics and regionalization in developing mouse brain.

Allen Brain Atlas (ABA) provides a valuable resource of spatial/temporal gene expressions in mammalian brains. Despite rich information extracted from this database, current analyses suffer from several limitations. First, most studies are either gene-centric or region-centric, thus are inadequate to capture the superposition of multiple spatial-temporal patterns. Second, standard tools of expression analysis such as matrix factorization can capture those patterns but do not explicitly incorporate spatial dependency. To overcome those limitations, we proposed a computational method to detect recurrent patterns in the spatial-temporal gene expression data of developing mouse brains. We demonstrated that regional distinction in brain development could be revealed by localized gene expression patterns. The patterns expressed in the forebrain, medullary and pontomedullary, and basal ganglia are enriched with genes involved in forebrain development, locomotory behavior, and dopamine metabolism respectively. In addition, the timing of global gene expression patterns reflects the general trends of molecular events in mouse brain development. Furthermore, we validated functional implications of the inferred patterns by showing genes sharing similar spatial-temporal expression patterns with Lhx2 exhibited differential expression in the embryonic forebrains of Lhx2 mutant mice. These analysis outcomes confirm the utility of recurrent expression patterns in studying brain development.

Focal amplification of HOXD-harboring chromosome region is implicated in multiple-walled carbon nanotubes-induced carcinogenicity.

Multiple-walled carbon nanotubes (MWCNTs) may cause carcinogenesis. We found that long-term exposure to MWCNTs can induce irreversible oncogenic transformation of human bronchial epithelial cells and tumorigenicity in vivo. A genome-wide array-comparative genomic hybridization (aCGH) analysis revealed global chromosomal aberration in MWCNTs-treated clones, predominantly at chromosome 2q31-32, where the potential oncogenes HOXD9 and HOXD13 are located. Functional assays confirmed that this variation can modulate oncogenic signaling and plays a part in MWCNTs-induced tumorigenesis, suggesting that MWCNTs are carcinogens that act by altering genomic stability and oncogenic copy numbers.

Tissue-specific gene expression templates for accurate molecular characterization of the normal physiological states of multiple human tissues with implication in development and cancer studies.

BACKGROUND: To elucidate the molecular complications in many complex diseases, we argue for the priority to construct a model representing the normal physiological state of a cell/tissue. RESULTS: By analyzing three independent microarray datasets on normal human tissues, we established a quantitative molecular model GET, which consists of 24 tissue-specific Gene Expression Templates constructed from a set of 56 genes, for predicting 24 distinct tissue types under disease-free condition. 99.2% correctness was reached when a large-scale validation was performed on 61 new datasets to test the tissue-prediction power of GET. Network analysis based on molecular interactions suggests a potential role of these 56 genes in tissue differentiation and carcinogenesis.Applying GET to transcriptomic datasets produced from tissue development studies the results correlated well with developmental stages. Cancerous tissues and cell lines yielded significantly lower correlation with GET than the normal tissues. GET distinguished melanoma from normal skin tissue or benign skin tumor with 96% sensitivity and 89% specificity. CONCLUSIONS: These results strongly suggest that a normal tissue or cell may uphold its normal functioning and morphology by maintaining specific chemical stoichiometry among genes. The state of stoichiometry can be depicted by a compact set of representative genes such as the 56 genes obtained here. A significant deviation from normal stoichiometry may result in malfunction or abnormal growth of the cells.

Clustered genomic alterations in chromosome 7p dictate outcomes and targeted treatment responses of lung adenocarcinoma with EGFR-activating mutations.

PURPOSE: Although epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs) have been proven more effective for patients with lung adenocarcinoma with EGFR-activating mutation rather than wild type, the former group still includes approximately 30% nonresponders. The molecular basis of this substantial response heterogeneity is unknown. Our purpose was to seek molecular aberrations contributing to disease progression at the genome-wide level and identify the prognostic signature unique to patients with EGFR-activating mutation. PATIENTS AND METHODS: We first investigated the molecular differences between tumors with EGFR-activating mutation and wild-type tumors by conducting high-density array comparative genomic hybridization on a collection of 138 adenocarcinoma tissues. We then used an independent group of 114 patients to validate the clinical relevance of copy-number alterations (CNAs) in predicting overall and disease-free survival. Finally, focusing on 23 patients with EGFR mutation receiving EGFR-TKI treatment, we investigated the association between CNAs and response to EGFR-TKIs. RESULTS: We identified chromosome regions with differential CNAs between tumors with EGFR-activating mutation and wild-type tumors and found the aberration sites to cluster highly on chromosome 7p. A cluster of six representative chromosome 7p genes predicted overall and disease-free survival for patients with EGFR-activating mutation but not for those with wild type. Importantly, simultaneous presence of more genes with increased CNAs in this cluster correlated with less favorable response to EGFR-TKIs in patients with EGFR-activating mutation. CONCLUSION: Our results shed light on why responses to EGFR-TKIs are heterogeneous among patients with EGFR-activating mutation. They may lead to better patient management in this population.

An optimized procedure greatly improves EST vector contamination removal.

BACKGROUND: The enormous amount of sequence data available in the public domain database has been a gold mine for researchers exploring various themes in life sciences, and hence the quality of such data is of serious concern to researchers. Removal of vector contamination is one of the most significant operations to obtain accurate sequence data containing only a cDNA insert from the basecalls output by an automatic DNA sequencer. Popular bioinformatics programs to accomplish vector trimming include LUCY, cross_match and SeqClean. RESULTS: In a recent study, where the program SeqClean was used to remove vector contamination from our test set of EST data compiled through various library construction systems, however, a significant number of errors remained after preliminary trimming. These errors were later almost completely corrected by simply using a re-linearized form of the cloning vector to compare against the target ESTs. The modified trimming procedure for SeqClean was also compared with the trimming efficiency of the other two popular programs, LUCY2, and cross_match. Using SeqClean with a re-linearized form of the cloning vector significantly surpassed the other two programs in all tested conditions, while the performance of the other two programs was not influenced by the modified procedure. Vector contamination in dbEST was also investigated in this study: 2203 out of the 48212 ESTs sampled from dbEST (2007-04-18 freeze) were found to match sequences in UNIVEC. CONCLUSION: Vector contamination remains a serious concern to the data quality in the public sequence database nowadays. Based on the results presented here, we feel that our modified procedure with SeqClean should be recommended to all researchers for the task of vector removal from EST or genomic sequences.