Wen-Ping Hsieh

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.

Quantitative trait transcripts for nicotine resistance in Drosophila melanogaster

Although most genetic association studies are performed with the intention of detecting nucleotide polymorphisms that are correlated with a complex trait, transcript abundance should also be expected to associate with diseases or phenotypes. We performed a scan for such quantitative trait transcripts in adult female heads of the fruit fly (Drosophila melanogaster) that might explain variation for nicotine resistance. The strongest association was seen for abundance of ornithine aminotransferase transcripts, implicating detoxification and neurotransmitter biosynthesis as mediators of the quantitative response to the drug. Subsequently, genetic analysis and metabolite profiling confirmed a complex role for ornithine and GABA levels in modification of survival time upon chronic nicotine exposure. Differences between populations from North Carolina and California suggest that the resistance mechanism may be an evolved response to environmental exposure.

Detection of transcriptional difference of porcine imprinted genes using different microarray platforms

BackgroundPresently, multiple options exist for conducting gene expression profiling studies in swine. In order to determine the performance of some of the existing microarrays, Affymetrix Porcine, Affymetrix Human U133+2.0, and the U.S. Pig Genome Coordination Program spotted glass oligonucleotide microarrays were compared for their reproducibility, coverage, platform independent and dependent sensitivity using fibroblast cell lines derived from control and parthenogenic porcine embryos.ResultsArray group correlations between technical replicates demonstrated comparable reproducibility in both Affymetrix arrays. Glass oligonucleotide arrays showed greater variability and, in addition, approximately 10% of probes had to be discarded due to slide printing defects. Probe level analysis of Affymetrix Human arrays revealed significant variability within probe sets due to the effects of cross-species hybridization. Affymetrix Porcine arrays identified the greatest number of differentially expressed genes amongst probes common to all arrays, a measure of platform sensitivity. Affymetrix Porcine arrays also identified the greatest number of differentially expressed known imprinted genes using all probes on each array, an ad hoc measure of realistic performance for this particular experiment.ConclusionWe conclude that of the platforms currently available and tested, the Affymetrix Porcine array is the most sensitive and reproducible microarray for swine genomic studies.

Global screening of potential Candida albicans biofilm-related transcription factors via network comparison.

BackgroundCandida albicans is a commonly encountered fungal pathogen in humans. The formation of biofilm is a major virulence factor in C. albicans pathogenesis and is related to antidrug resistance of this organism. Although many factors affecting biofilm have been analyzed, molecular mechanisms that regulate biofilm formation still await to be elucidated.ResultsIn this study, from the gene regulatory network perspective, we developed an efficient computational framework, which integrates different kinds of data from genome-scale analysis, for global screening of potential transcription factors (TFs) controlling C. albicans biofilm formation. S. cerevisiae information and ortholog data were used to infer the possible TF-gene regulatory associations in C. albicans. Based on TF-gene regulatory associations and gene expression profiles, a stochastic dynamic model was employed to reconstruct the gene regulatory networks of C. albicans biofilm and planktonic cells. The two networks were then compared and a score of relevance value (RV) was proposed to determine and assign the quantity of correlation of each potential TF with biofilm formation. A total of twenty-three TFs are identified to be related to the biofilm formation; ten of them are previously reported by literature evidences.ConclusionsThe results indicate that the proposed screening method can successfully identify most known biofilm-related TFs and also identify many others that have not been previously reported. Together, this method can be employed as a pre-experiment screening approach that reveals new target genes for further characterization to understand the regulatory mechanisms in biofilm formation, which can serve as the starting point for therapeutic intervention of C. albicans infections.

Responses of Candida albicans to the human antimicrobial peptide LL-37

Candida albicans is amajor fungal pathogen in humans. Antimicrobial peptides (AMPs) are critical components of the innate immune response in vertebrates and represent the first line of defense against microbial infection. LL-37 is the only member of the human family of cathelicidin AMPs and is commonly expressed by various tissues and cells, including surfaces of epithelia. The candidacidal effects of LL-37 have been well documented, but the mechanisms by which LL-37 kills C. albicans are not completely understood. In this study, we examined the effects of LL-37 on cell wall and cellular responses in C. albicans. Using transmission electron microscopy, carbohydrate analyses, and staining for β-1,3-glucan, changing of C. albicans cell wall integrity was detected upon LL-37 treatment. In addition, LL-37 also affected cell wall architecture of the pathogen. Finally, DNA microarray analysis and quantitative PCR demonstrated that sub-lethal concentrations of LL-37 modulated the expression of genes with a variety of functions, including transporters, regulators for biological processes, response to stress or chemical stimulus, and pathogenesis. Together, LL-37 induces complex responses in C. albicans, making LL-37 a promising candidate for use as a therapeutic agent against fungal infections.

Dynamic Transcript Profiling of Candida albicans Infection in Zebrafish: A Pathogen-Host Interaction Study

Candida albicans is responsible for a number of life-threatening infections and causes considerable morbidity and mortality in immunocompromised patients. Previous studies of C. albicans pathogenesis have suggested several steps must occur before virulent infection, including early adhesion, invasion, and late tissue damage. However, the mechanism that triggers C. albicans transformation from yeast to hyphae form during infection has yet to be fully elucidated. This study used a systems biology approach to investigate C. albicans infection in zebrafish. The surviving fish were sampled at different post-infection time points to obtain time-lapsed, genome-wide transcriptomic data from both organisms, which were accompanied with in sync histological analyses. Principal component analysis (PCA) was used to analyze the dynamic gene expression profiles of significant variations in both C. albicans and zebrafish. The results categorized C. albicans infection into three progressing phases: adhesion, invasion, and damage. Such findings were highly supported by the corresponding histological analysis. Furthermore, the dynamic interspecies transcript profiling revealed that C. albicans activated its filamentous formation during invasion and the iron scavenging functions during the damage phases, whereas zebrafish ceased its iron homeostasis function following massive hemorrhage during the later stages of infection. Most of the immune related genes were expressed as the infection progressed from invasion to the damage phase. Such global, inter-species evidence of virulence-immune and iron competition dynamics during C. albicans infection could be crucial in understanding control fungal pathogenesis.

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.

Interspecies protein-protein interaction network construction for characterization of host-pathogen interactions: a Candida albicans-zebrafish interaction study

BackgroundDespite clinical research and development in the last decades, infectious diseases remain a top global problem in public health today, being responsible for millions of morbidities and mortalities each year. Therefore, many studies have sought to investigate host-pathogen interactions from various viewpoints in attempts to understand pathogenic and defensive mechanisms, which could help control pathogenic infections. However, most of these efforts have focused predominately on the host or the pathogen individually rather than on a simultaneous analysis of both interaction partners.ResultsIn this study, with the help of simultaneously quantified time-course Candida albicans-zebrafish interaction transcriptomics and other omics data, a computational framework was developed to construct the interspecies protein-protein interaction (PPI) network for C. albicans-zebrafish interactions based on the inference of ortholog-based PPIs and the dynamic modeling of regulatory responses. The identified C. albicans-zebrafish interspecies PPI network highlights the association between C. albicans pathogenesis and the zebrafish redox process, indicating that redox status is critical in the battle between the host and pathogen.ConclusionsAdvancing from the single-species network construction method, the interspecies network construction approach allows further characterization and elucidation of the host-pathogen interactions. With continued accumulation of interspecies transcriptomics data, the proposed method could be used to explore progressive network rewiring over time, which could benefit the development of network medicine for infectious diseases.

The chromatin modification by SUMO-2/3 but not SUMO-1 prevents the epigenetic activation of key immune-related genes during Kaposi’s sarcoma associated herpesvirus reactivation

BackgroundSUMOylation, as part of the epigenetic regulation of transcription, has been intensively studied in lower eukaryotes that contain only a single SUMO protein; however, the functions of SUMOylation during mammalian epigenetic transcriptional regulation are largely uncharacterized. Mammals express three major SUMO paralogues: SUMO-1, SUMO-2, and SUMO-3 (normally referred to as SUMO-1 and SUMO-2/3). Herpesviruses, including Kaposi’s sarcoma associated herpesvirus (KSHV), seem to have evolved mechanisms that directly or indirectly modulate the SUMO machinery in order to evade host immune surveillance, thus advancing their survival. Interestingly, KSHV encodes a SUMO E3 ligase, K-bZIP, with specificity toward SUMO-2/3 and is an excellent model for investigating the global functional differences between SUMO paralogues.ResultsWe investigated the effect of experimental herpesvirus reactivation in a KSHV infected B lymphoma cell line on genomic SUMO-1 and SUMO-2/3 binding profiles together with the potential role of chromatin SUMOylation in transcription regulation. This was carried out via high-throughput sequencing analysis. Interestingly, chromatin immunoprecipitation sequencing (ChIP-seq) experiments showed that KSHV reactivation is accompanied by a significant increase in SUMO-2/3 modification around promoter regions, but SUMO-1 enrichment was absent. Expression profiling revealed that the SUMO-2/3 targeted genes are primarily highly transcribed genes that show no expression changes during viral reactivation. Gene ontology analysis further showed that these genes are involved in cellular immune responses and cytokine signaling. High-throughput annotation of SUMO occupancy of transcription factor binding sites (TFBS) pinpointed the presence of three master regulators of immune responses, IRF-1, IRF-2, and IRF-7, as potential SUMO-2/3 targeted transcriptional factors after KSHV reactivation.ConclusionOur study is the first to identify differential genome-wide SUMO modifications between SUMO paralogues during herpesvirus reactivation. Our findings indicate that SUMO-2/3 modification near protein-coding gene promoters occurs in order to maintain host immune-related gene unaltered during viral reactivation.

Mixture modeling of transcript abundance classes in natural populations

BackgroundPopulations diverge in genotype and phenotype under the influence of such evolutionary processes as genetic drift, mutation accumulation, and natural selection. Because genotype maps onto phenotype by way of transcription, it is of interest to evaluate how these evolutionary factors influence the structure of variation at the level of transcription. Here, we explore the distributions of cis-acting and trans-acting factors and their relative contributions to expression of transcripts that exhibit two or more classes of abundance among individuals within populations.ResultsExpression profiling using cDNA microarrays was conducted in Drosophila melanogaster adult female heads for 58 nearly isogenic lines from a North Carolina population and 50 from a California population. Using a mixture modeling approach, transcripts were identified that exhibit more than one mode of transcript abundance across the samples. Power studies indicate that sample sizes of 50 individuals will generally be sufficient to detect divergent transcript abundance classes. The distribution of transcript abundance classes is skewed toward low frequency minor classes, which is reminiscent of the typical skew in genotype frequencies. Similar results are observed in reported data on gene expression in human lymphoblast cell lines, in which analysis of association with linked polymorphisms implies that cis-acting single nucleotide polymorphisms make only a modest contribution to bimodal distributions of transcript abundance.ConclusionPopulation surveys of gene expression may complement genetical genomics as a general approach to quantifying sources of transcriptional variation. Differential expression of transcripts among individuals is due to a complex interplay of cis-acting and trans-acting factors.