Decomposing the subclonal structure of tumors with two-way mixture models on copy number aberrations.

Multistage tumorigenesis is a dynamic process characterized by the accumulation of mutations. Thus, a tumor mass is composed of genetically divergent cell subclones. With the advancement of next-generation sequencing (NGS), mathematical models have been recently developed to decompose tumor subclonal architecture from a collective genome sequencing data. Most of the methods focused on single-nucleotide variants (SNVs). However, somatic copy number aberrations (CNAs) also play critical roles in carcinogenesis. Therefore, further modeling subclonal CNAs composition would hold the promise to improve the analysis of tumor heterogeneity and cancer evolution. To address this issue, we developed a two-way mixture Poisson model, named CloneDeMix for the deconvolution of read-depth information. It can infer the subclonal copy number, mutational cellular prevalence (MCP), subclone composition, and the order in which mutations occurred in the evolutionary hierarchy. The performance of CloneDeMix was systematically assessed in simulations. As a result, the accuracy of CNA inference was nearly 93% and the MCP was also accurately restored. Furthermore, we also demonstrated its applicability using head and neck cancer samples from TCGA. Our results inform about the extent of subclonal CNA diversity, and a group of candidate genes that probably initiate lymph node metastasis during tumor evolution was also discovered. Most importantly, these driver genes are located at 11q13.3 which is highly susceptible to copy number change in head and neck cancer genomes. This study successfully estimates subclonal CNAs and exhibit the evolutionary relationships of mutation events. By doing so, we can track tumor heterogeneity and identify crucial mutations during evolution process. Hence, it facilitates not only understanding the cancer development but finding potential therapeutic targets. Briefly, this framework has implications for improved modeling of tumor evolution and the importance of inclusion of subclonal CNAs.

Somatic copy number alterations detected by ultra-deep targeted sequencing predict prognosis in oral cavity squamous cell carcinoma.

BACKGROUND: Ultra-deep targeted sequencing (UDT-Seq) has advanced our knowledge on the incidence and functional significance of somatic mutations. However, the utility of UDT-Seq in detecting copy number alterations (CNAs) remains unclear. With the goal of improving molecular prognostication and identifying new therapeutic targets, we designed this study to assess whether UDT-Seq may be useful for detecting CNA in oral cavity squamous cell carcinoma (OSCC). METHODS: We sequenced a panel of clinically actionable cancer mutations in 310 formalin-fixed paraffin-embedded OSCC specimens. A linear model was developed to overcome uneven coverage across target regions and multiple samples. The 5-year rates of secondary primary tumors, local recurrence, neck recurrence, distant metastases, and survival served as the outcome measures. We confirmed the prognostic significance of the CNA signatures in an independent sample of 105 primary OSCC specimens. RESULTS: The CNA burden across 10 targeted genes was found to predict prognosis in two independent cohorts. FGFR1 and PIK3CAamplifications were associated with prognosis independent of clinical risk factors. Genes exhibiting CNA were clustered in the proteoglycan metabolism, the FOXO signaling, and the PI3K-AKT signaling pathways, for which targeted drugs are already available or currently under development. CONCLUSIONS: UDT-Seq is clinically useful to identify CNA, which significantly improve the prognostic information provided by traditional clinicopathological risk factors in OSCC patients.

MicroRNAs MiR-218, MiR-125b, and Let-7g predict prognosis in patients with oral cavity squamous cell carcinoma.

MicroRNAs (miRNAs) have a major impact on regulatory networks in human carcinogenesis. In this study, we sought to investigate the prognostic significance of miRNAs in patients with oral cavity squamous cell carcinoma (OSCC). In a discovery phase, RNA was extracted from 58 OSCC tumor samples and paired normal tissues. MiRNAs expression was evaluated with TaqMan Array Card and TaqMan MicroRNA assays. The prognostic significance of the miRNA signature identified in the discovery phase was validated by qRT-PCR in a replication set consisting of 141 formalin-fixed, paraffin-embedded (FFPE) samples. We identified a miRNA regulatory network centered on the three hub genes (SP1, MYC, and TP53) that predicted distinct clinical endpoints. Three miRNAs (miR-218, miR-125b, and let-7g) and their downstream response genes had a concordant prognostic significance on disease-free survival and disease-specific survival rates. In addition, patients with a reduced expression of miR-218, miR-125b, and let-7g have a higher risk of poor outcomes in presence of specific risk factors (p-stage III-IV, pT3-4, or pN+). Our findings indicate that specific miRNAs have prognostic significance in OSCC patients and may improve prognostic stratification over traditional risk factors.

HaplotypeCN: copy number haplotype inference with Hidden Markov Model and localized haplotype clustering.

Copy number variation (CNV) has been reported to be associated with disease and various cancers. Hence, identifying the accurate position and the type of CNV is currently a critical issue. There are many tools targeting on detecting CNV regions, constructing haplotype phases on CNV regions, or estimating the numerical copy numbers. However, none of them can do all of the three tasks at the same time. This paper presents a method based on Hidden Markov Model to detect parent specific copy number change on both chromosomes with signals from SNP arrays. A haplotype tree is constructed with dynamic branch merging to model the transition of the copy number status of the two alleles assessed at each SNP locus. The emission models are constructed for the genotypes formed with the two haplotypes. The proposed method can provide the segmentation points of the CNV regions as well as the haplotype phasing for the allelic status on each chromosome. The estimated copy numbers are provided as fractional numbers, which can accommodate the somatic mutation in cancer specimens that usually consist of heterogeneous cell populations. The algorithm is evaluated on simulated data and the previously published regions of CNV of the 270 HapMap individuals. The results were compared with five popular methods: PennCNV, genoCN, COKGEN, QuantiSNP and cnvHap. The application on oral cancer samples demonstrates how the proposed method can facilitate clinical association studies. The proposed algorithm exhibits comparable sensitivity of the CNV regions to the best algorithm in our genome-wide study and demonstrates the highest detection rate in SNP dense regions. In addition, we provide better haplotype phasing accuracy than similar approaches. The clinical association carried out with our fractional estimate of copy numbers in the cancer samples provides better detection power than that with integer copy number states.

Causal inference of gene regulation with subnetwork assembly from genetical genomics data.

Deciphering the causal networks of gene interactions is critical for identifying disease pathways and disease-causing genes. We introduce a method to reconstruct causal networks based on exploring phenotype-specific modules in the human interactome and including the expression quantitative trait loci (eQTLs) that underlie the joint expression variation of each module. Closely associated eQTLs help anchor the orientation of the network. To overcome the inherent computational complexity of causal network reconstruction, we first deduce the local causality of individual subnetworks using the selected eQTLs and module transcripts. These subnetworks are then integrated to infer a global causal network using a random-field ranking method, which was motivated by animal sociology. We demonstrate how effectively the inferred causality restores the regulatory structure of the networks that mediate lymph node metastasis in oral cancer. Network rewiring clearly characterizes the dynamic regulatory systems of distinct disease states. This study is the first to associate an RXRB-causal network with increased risks of nodal metastasis, tumor relapse, distant metastases and poor survival for oral cancer. Thus, identifying crucial upstream drivers of a signal cascade can facilitate the discovery of potential biomarkers and effective therapeutic targets.

Feature Identification of Compensatory Gene Pairs without Sequence Homology in Yeast.

Genetic robustness refers to a compensatory mechanism for buffering deleterious mutations or environmental variations. Gene duplication has been shown to provide such functional backups. However, the overall contribution of duplication-based buffering for genetic robustness is rather small. In this study, we investigated whether transcriptional compensation also exists among genes that share similar functions without sequence homology. A set of nonhomologous synthetic-lethal gene pairs was assessed by using a coexpression network, protein-protein interactions, and other types of genetic interactions in yeast. Our results are notably different from those of previous studies on buffering paralogs. The low expression similarity and the conditional coexpression alone do not play roles in identifying the functionally compensatory genes. Additional properties such as synthetic-lethal interaction, the ratio of shared common interacting partners, and the degree of coregulation were, at least in part, necessary to extract functional compensatory genes. Our network-based approach is applicable to select several well-documented cases of compensatory gene pairs and a set of new pairs. The results suggest that transcriptional reprogramming plays a limited role in functional compensation among nonhomologous genes. Our study aids in understanding the mechanism and features of functional compensation more in detail.

A novel molecular signature identified by systems genetics approach predicts prognosis in oral squamous cell carcinoma.

Molecular methods for predicting prognosis in patients with oral cavity squamous cell carcinoma (OSCC) are urgently needed, considering its high recurrence rate and tendency for metastasis. The present study investigated the genetic basis of variations in gene expression associated with poor prognosis in OSCC using Affymetrix SNP 6.0 and Affymetrix GeneChip Human Gene 1.0 ST arrays. We identified recurrent DNA amplifications scattered from 8q22.2 to 8q24.3 in 112 OSCC specimens. These amplicons demonstrated significant associations with increased incidence of extracapsular spread, development of second primary malignancies, and poor survival. Fluorescence in situ hybridization, in a validation panel consisting of 295 cases, confirmed these associations. Assessment of the effects of copy number variations (CNVs) on genome-wide variations in gene expression identified a total of 85 CNV-associated transcripts enriched in the MYC-centered regulatory network. Twenty-four transcripts associated with increased risk of second primary malignancies, tumor relapse, and poor survival. Besides MYC itself, a novel dysregulated MYC module plays a key role in OSCC carcinogenesis. This study identified a candidate molecular signature associated with poor prognosis in OSCC patients, which may ultimately facilitate patient-tailored selection of therapeutic strategies.

A novel method to identify cooperative functional modules: study of module coordination in the Saccharomyces cerevisiae cell cycle.

BACKGROUND: Identifying key components in biological processes and their associations is critical for deciphering cellular functions. Recently, numerous gene expression and molecular interaction experiments have been reported in Saccharomyces cerevisiae, and these have enabled systematic studies. Although a number of approaches have been used to predict gene functions and interactions, tools that analyze the essential coordination of functional components in cellular processes still need to be developed. RESULTS: In this work, we present a new approach to study the cooperation of functional modules (sets of functionally related genes) in a specific cellular process. A cooperative module pair is defined as two modules that significantly cooperate with certain functional genes in a cellular process. This method identifies cooperative module pairs that significantly influence a cellular process and the correlated genes and interactions that are essential to that process. Using the yeast cell cycle as an example, we identified 101 cooperative module associations among 82 modules, and importantly, we established a cell cycle-specific cooperative module network. Most of the identified module pairs cover cooperative pathways and components essential to the cell cycle. We found that 14, 36, 18, 15, and 20 cooperative module pairs significantly cooperate with genes regulated in early G1, late G1, S, G2, and M phase, respectively. Fifty-nine module pairs that correlate with Cdc28 and other essential regulators were also identified. These results are consistent with previous studies and demonstrate that our methodology is effective for studying cooperative mechanisms in the cell cycle. CONCLUSIONS: In this work, we propose a new approach to identifying condition-related cooperative interactions, and importantly, we establish a cell cycle-specific cooperation module network. These results provide a global view of the cell cycle and the method can be used to discover the dynamic coordination properties of functional components in other cellular processes.

Computational modeling with forward and reverse engineering links signaling network and genomic regulatory responses: NF-kappaB signaling-induced gene expression responses in inflammation.

BACKGROUND: Signal transduction is the major mechanism through which cells transmit external stimuli to evoke intracellular biochemical responses. Diverse cellular stimuli create a wide variety of transcription factor activities through signal transduction pathways, resulting in different gene expression patterns. Understanding the relationship between external stimuli and the corresponding cellular responses, as well as the subsequent effects on downstream genes, is a major challenge in systems biology. Thus, a systematic approach is needed to integrate experimental data and theoretical hypotheses to identify the physiological consequences of environmental stimuli. RESULTS: We proposed a systematic approach that combines forward and reverse engineering to link the signal transduction cascade with the gene responses. To demonstrate the feasibility of our strategy, we focused on linking the NF-kappaB signaling pathway with the inflammatory gene regulatory responses because NF-kappaB has long been recognized to play a crucial role in inflammation. We first utilized forward engineering (Hybrid Functional Petri Nets) to construct the NF-kappaB signaling pathway and reverse engineering (Network Components Analysis) to build a gene regulatory network (GRN). Then, we demonstrated that the corresponding IKK profiles can be identified in the GRN and are consistent with the experimental validation of the IKK kinase assay. We found that the time-lapse gene expression of several cytokines and chemokines (TNF-alpha, IL-1, IL-6, CXCL1, CXCL2 and CCL3) is concordant with the NF-kappaB activity profile, and these genes have stronger influence strength within the GRN. Such regulatory effects have highlighted the crucial roles of NF-kappaB signaling in the acute inflammatory response and enhance our understanding of the systemic inflammatory response syndrome. CONCLUSION: We successfully identified and distinguished the corresponding signaling profiles among three microarray datasets with different stimuli strengths. In our model, the crucial genes of the NF-kappaB regulatory network were also identified to reflect the biological consequences of inflammation. With the experimental validation, our strategy is thus an effective solution to decipher cross-talk effects when attempting to integrate new kinetic parameters from other signal transduction pathways. The strategy also provides new insight for systems biology modeling to link any signal transduction pathways with the responses of downstream genes of interest.

Mechanomyography is more sensitive than EMG in detecting age-related sarcopenia.

The aim of this work was to investigate the effects of age-related sarcopenia on the time and frequency domain properties of lower extremity muscles' electromyographic and mechanomyographic activities. Healthy elderly (n=10, 64.5+/-4.5yr) and young (n=10, 22.6+/-2.8yr) were recruited as participants. Participants' lean thigh volumes (LTV) and 1 RM (one repetition maximum) leg strength of quadriceps and maximum speed knee extension with different load levels (45%, 60% and 75% 1 RM) were recorded. The root mean square (RMS) and the mean frequency (MF) of the surface electromyography (EMG(RMS), EMG(MF)) and mechanomyography (MMG(RMS), MMG(MF)) signals were collected at vastus lateralis during concentric contraction with different intensity levels. Compared to the young, the elderly had significantly less LTV, absolute and relative maximal force, as well as absolute and relative maximal power (p<.05). EMG(MF) of the elderly and the young increased monotonically from 45% to 75% 1 RM testing conditions. While the MMG(RMS) of the young increased with testing intensities, the MMG(RMS) of the elderly increased only from 45% to 60% but leveled off from 60% to 75% 1 RM testing conditions. The results indicate the declines of muscle mass, force and power production capacity with aging. The observations could be explained by neuromuscular performance and change of MU activation patterns may result from age-related sarcopenia. Aging affected muscle power more than muscle strength, which could be due to fast fiber reduction. This is supported by our observations that the MMG(RMS) differences between the young and the elderly across all three intensity level where EMG(RMS) was only different at the greatest intensity. We suggest that MMG could be used as an important measurement in studying muscle contraction in age-related sarcopenia.

Identification of degenerate motifs using position restricted selection and hybrid ranking combination.

The identification of regulatory elements recognized by transcription factors and chromatin remodeling factors is essential to studying the regulation of gene expression. When no auxiliary data, such as orthologous sequences or expression profiles, are used, the accuracy of most tools for motif discovery is strongly influenced by the motif degeneracy and the lengths of sequence. Since suitable auxiliary data may not always be available, more work must be conducted to enhance tool performance to identify transcription elements in the metazoan. A non-alignment-based algorithm, MotifSeeker, is proposed to enhance the accuracy of discovering degenerate motifs. MotifSeeker utilizes the property that variable sites of transcription elements are usually position-specific to reduce exposure to noise. Consequently, the efficiency and accuracy of motif identification are improved. Using data fusion, the ranking process integrates two measures of motif significance, resulting in a more robust significance measure. Testing results for the synthetic data reveal that the accuracy of MotifSeeker is less sensitive to the motif degeneracy and the length of input sequences. Furthermore, MotifSeeker has been tested on a well-known benchmark [M. Tompa, N. Li, T.L. Bailey, G.M. Church, B. De Moor, E. Eskin, A.V. Favorov, M.C. Frith, Y. Fu, W.J. Kent, et al. (2005) Nat. Biotechnol., 23, 137-144], yielding a correlation coefficient of 0.262, which compares favorably with those of other tools. The high applicability of MotifSeeker to biological data is further demonstrated experimentally on regulons of Saccharomyces cerevisiae and liver-specific genes with experimentally verified regulatory elements.

Active leg stiffness and energy stored in the muscles during maximal counter movement jump in the aged.

The purpose of this study was to examine the effects of age on active leg stiffness adjustment, electromyogram (EMG) activities and energy stored during eccentric and concentric phases in performing a maximal functional task involving stretch-shorten cycle. Ten young (24.3+/-2 years) and 10 old (68.6+/-5 years) healthy male subjects were filmed during maximal performance of counter movement jump (CMJ) and squat jump (SJ) on force plate. Integrated EMG (IEMG), ground reaction force (GRF), active leg stiffness, energy stored/returned and active work done by the muscles were compared between two groups on eccentric (ECC) and concentric (CON) phases of CMJ. The GRF, leg stiffness and energy stored in ECC and GRF, IEMG, energy returned and active work in CON were less in the elderly (p<0.05). These results demonstrate that the neuromuscular function of adjusting active stiffness, storing elastic energy and optimizing the performance may decrease with age during CMJ.