Genetic Predictors of Clinical Outcomes in Non-Small Cell Lung Cancer Patients

Lung cancer is the leading cause of cancer-related mortality in the US. Emerging evidence has shown that host genetic factors can interact with environmental exposures to influence patient susceptibility to the diseases as well as clinical outcomes, such as survival and recurrence. We aimed to identify genetic prognostic markers for non-small cell lung cancer (NSCLC), a major (85%) subtype of lung cancer, and also in other subgroups. With the fast evolution of genotyping technology, genetic association studies have went through candidate gene approach, to pathway-based approach, to the genome wide association study (GWAS). Even in the era of GWAS, pathway-based approach has its own advantages on studying cancer clinical outcomes: it is cost-effective, requiring a smaller sample size than GWAS easier to identify a validation population and explore gene-gene interactions. In the current study, we adopted pathway-based approach focusing on two critical pathways - miRNA and inflammation pathways. MicroRNAs (miRNA) post-transcriptionally regulate around 30% of human genes. Polymorphisms within miRNA processing pathways and binding sites may influence patients’ prognosis through altered gene regulation. Inflammation plays an important role in cancer initiation and progression, and also has shown to impact patients’ clinical outcomes. We first evaluated 240 single nucleotide polymorphisms (SNPs) in miRNA biogenesis genes and predicted binding sites in NSCLC patients to determine associations with clinical outcomes in early-stage (stage I and II) and late-stage (stage III and IV) lung cancer patients, respectively. First, in 535 early-stage patients, after correcting multiple comparisons, FZD4:rs713065 (hazard ratio [HR]:0.46, 95% confidence interval [CI]:0.32-0.65) showed a significant inverse association with survival in early stage surgery-only patients. SP1:rs17695156 (HR:2.22, 95% CI:1.44-3.41) and DROSHA:rs6886834 (HR:6.38, 95% CI:2.49-16.31) conferred increased risk of progression in the all patients and surgery-only populations, respectively. FAS:rs2234978 was significantly associated with improved survival in all patients (HR:0.59, 95% CI:0.44-0.77) and in the surgery plus chemotherapy populations (HR:0.19, 95% CI:0.07-0.46).. Functional genomics analysis demonstrated that this variant creates a miR-651 binding site resulting in altered miRNA regulation of FAS, providing biological plausibility for the observed association. We then analyzed these associations in 598 late-stage patients. After multiple comparison corrections, no SNPs remained significant in the late stage group, while the top SNP NAT1:rs15561 (HR=1.98, 96%CI=1.32-2.94) conferred a significantly increased risk of death in the chemotherapy subgroup. To test the hypothesis that genetic variants in the inflammation-related pathways may be associated with survival in NSCLC patients, we first conducted a three-stage study. In the discovery phase, we investigated a comprehensive panel of 11,930 inflammation-related SNPs in three independent lung cancer populations. A missense SNP (rs2071554) in HLA-DOB was significantly associated with poor survival in the discovery population (HR: 1.46, 95% CI: 1.02-2.09), internal validation population (HR: 1.51, 95% CI: 1.02-2.25), and external validation (HR: 1.52, 95% CI: 1.01-2.29) population. Rs2900420 in KLRK1 was significantly associated with a reduced risk for death in the discovery (HR: 0.76, 95% CI: 0.60-0.96) and internal validation (HR: 0.77, 95% CI: 0.61-0.99) populations, and the association reached borderline significance in the external validation population (HR: 0.80, 95% CI: 0.63-1.02). We also evaluated these inflammation-related SNPs in NSCLC patients in never smokers. Lung cancer in never smokers has been increasingly recognized as distinct disease from that in ever-smokers. A two-stage study was performed using a discovery population from MD Anderson (411 patients) and a validation population from Mayo Clinic (311 patients). Three SNPs (IL17RA:rs879576, BMP8A:rs698141, and STK:rs290229) that were significantly associated with survival were validated (pCD74:rs1056400 and CD38:rs10805347) were borderline significant (p=0.08) in the Mayo Clinic population. In the combined analysis, IL17RA:rs879576 resulted in a 40% reduction in the risk for death (p=4.1 × 10-5 [p=0.61, heterogeneity test]). We also validated a survival tree created in MD Anderson population in the Mayo Clinic population. In conclusion, our results provided strong evidence that genetic variations in specific pathways that examined (miRNA and inflammation pathways) influenced clinical outcomes in NSCLC patients, and with further functional studies, the novel loci have potential to be translated into clinical use.

Critical determinants of estrogen receptor -mediated agonist activity

Exogenous ligands that bind to the estrogen receptor (ER) exhibit unique pharmacologies distinct from that observed with the endogenous hormone, 17β-estradiol (ED. Differential activity among ER ligands has been observed at the level of receptor binding, promoter interaction and transcriptional activation. Furthermore, xenoestrogens can display tissue-specific agonist activity on the cellular level, functioning as an agonist in one tissue and as an antagonist in another. That the same ligand, functioning through the same receptor, can produce differing agonist responses on the cellular level indicates that there are tissue-specific determinants of agonist activity. In these studies critical molecular determinants of agonist activity were characterized for several cell types. In the normal and neoplastic myometrium a proliferative response was dependent upon activation of AF2 of the ER, functioning as a determinant of agonism in this cell type. Progesterone receptor (PR) ligands transdominantly suppressed ER-mediated transcription and proliferation in uterine leiomyoma cells, indicating that ER/PR cross-talk can modulate agonist activity in a myometrial cell background. In the breast, the agonist response to ER ligands was investigated by employing a functional genomics approach to generate gene expression profiles. Treatment of breast cancer cells with the selective estrogen receptor modulator tamoxifen largely recapitulated the expression profile induced by treatment with the agonist E2, despite the well-characterized antiproliferative effects produced by tamoxifen in this cell type. While the expression of many genes involved in regulating cell cycle progression, including fos, myc, cdc25a, stk15 and cyclin A, were induced by both E2 and tamoxifen in breast cells, treatment with the agonist E2 specifically induced the expression of cyclin D1, fra-1 , and uracil DNA glycosylase. These results suggest that the inability of tamoxifen to transactivate expression of only a few key genes, functioning as cellular gatekeepers, prevent tamoxifen-treated breast cells from entering the cell cycle. Thus, the expression of these agonist-specific marker genes is a potential determinant of agonist activity at the cellular level in the breast. Collectively, studies in the breast and uterine myometrium have identified several mechanisms whereby ER ligands modulate ER-mediated signaling and provide insights into the biology of tissue-specific agonist activity in hormone-responsive tissues. ^

Thermographic imaging: assessment of drought and heat tolerance in Spanish germplasm of Brachypodium distachyon.

The annual grass Brachypodium distachyon has been recently recognized as the model plant for functional genomics of temperate grasses, including cereals of economic relevance like wheat and barley. Sixty-two lines of B. distachyon were assessed for response to drought stress and heat tolerance. All these lines, except the reference genotype BD21, derive from specimens collected in 32 distinct locations of the Iberian Peninsula, covering a wide range of geo- climatic conditions. Sixteen lines of Brachypodium hybridum, an allotetraploid closely related to B. distachyon were used as reference of abiotic-stress well-adapted genotypes. Drought tolerance was assessed in a green-house trial. At the rosette-stage, no irrigation was applied to treated plants whereas their replicates at the control were maintained well watered during all the experiment. Thermographic images of treated and control plants were taken after 2 and 3 weeks of drought treatment, when stressed plants showed medium and extreme wilting symptoms. The mean leaf temperature of stressed (LTs) and control (LTc) plants was estimated based upon thermographic records from selected pixels (183 per image) that strictly correspond to leaf tissue. The response to drought was based on the analysis of two parameters: LTs and the thermal difference (TD) between stressed and control plants (LTs – LTc). The response to heat stress was based on LTc. Comparison of the mean values of these parameters showed that: 1) Genotypes better adapted to drought (B. hybridum lines) presented a higher LTs and TD than B. distachyon lines. 2) Under high temperature conditions, watered plants of B. hybridum lines maintained lower LTc than those of B. distachyon. Those results suggest that in these species adaptation to drought is linked to a more efficient stomata regulation: under water stress stomata are closed, increasing foliar temperature but also water use efficiency by reducing transpiration. With high temperature and water availability the results are less definite, but still seems that opening stomata allow plants to increase transpiration and therefore to diminish foliar temperature.

Toward controlling gene expression at will: Specific regulation of the erbB-2/HER-2 promoter by using polydactyl zinc finger proteins constructed from modular building blocks

To create a universal system for the control of gene expression, we have studied methods for the construction of novel polydactyl zinc finger proteins that recognize extended DNA sequences. Elsewhere we have described the generation of zinc finger domains recognizing sequences of the 5′-GNN-3′ subset of a 64-member zinc finger alphabet. Here we report on the use of these domains as modular building blocks for the construction of polydactyl proteins specifically recognizing 9- or 18-bp sequences. A rapid PCR assembly method was developed that, together with this predefined set of zinc finger domains, provides ready access to 17 million novel proteins that bind the 5′-(GNN)6-3′ family of 18-bp DNA sites. To examine the efficacy of this strategy in gene control, the human erbB-2 gene was chosen as a model. A polydactyl protein specifically recognizing an 18-bp sequence in the 5′-untranslated region of this gene was converted into a transcriptional repressor by fusion with Krüppel-associated box (KRAB), ERD, or SID repressor domains. Transcriptional activators were generated by fusion with the herpes simplex VP16 activation domain or with a tetrameric repeat of VP16’s minimal activation domain, termed VP64. We demonstrate that both gene repression and activation can be achieved by targeting designed proteins to a single site within the transcribed region of a gene. We anticipate that gene-specific transcriptional regulators of the type described here will find diverse applications in gene therapy, functional genomics, and the generation of transgenic organisms.

Genetic footprinting with mariner-based transposition in Pseudomonas aeruginosa

The complete DNA sequence of Pseudomonas aeruginosa provides an opportunity to apply functional genomics to a major human pathogen. A comparative genomics approach combined with genetic footprinting was used as a strategy to identify genes required for viability in P. aeruginosa. Use of a highly efficient in vivo mariner transposition system in P. aeruginosa facilitated the analysis of candidate genes of this class. We have developed a rapid and efficient allelic exchange system by using the I-SceI homing endonuclease in conjunction with in vitro mariner mutagenesis to generate mutants within targeted regions of the P. aeruginosa chromosome for genetic footprinting analyses. This technique for generating transposon insertion mutants should be widely applicable to other organisms that are not naturally transformable or may lack well developed in vivo transposition systems. We tested this system with three genes in P. aeruginosa that have putative essential homologs in Haemophilus influenzae. We show that one of three H. influenzae essential gene homologs is needed for growth in P. aeruginosa, validating the practicality of this comparative genomics strategy to identify essential genes in P. aeruginosa.

Chromosomal transposition of a Tc1/mariner-like element in mouse embryonic stem cells

Mouse has become an increasingly important organism for modeling human diseases and for determining gene function in a mammalian context. Unfortunately, transposon-tagged mutagenesis, one of the most valuable tools for functional genomics, still is not available in this organism. On the other hand, it has long been speculated that members of the Tc1/mariner-like elements may be less dependent on host factors and, hence, can be introduced into heterologous organisms. However, this prediction has not been realized in mice. We report here the chromosomal transposition of the Sleeping Beauty (SB) element in mouse embryonic stem cells, providing evidence that it can be used as an in vivo mutagen in mice.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

The Medicago Genome Initiative: a model legume database

The Medicago Genome Initiative (MGI) is a database of EST sequences of the model legume Medicago truncatula. The database is available to the public and has resulted from a collaborative research effort between the Samuel Roberts Noble Foundation and the National Center for Genome Resources to investigate the genome of M.truncatula. MGI is part of the greater integrated Medicago functional genomics program at the Noble Foundation (http://www.noble .org), which is taking a global approach in studying the genetic and biochemical events associated with the growth, development and environmental interactions of this model legume. Our approach will include: large-scale EST sequencing, gene expression profiling, the generation of M.truncatula activation-tagged and promoter trap insertion mutants, high-throughput metabolic profiling, and proteome studies. These multidisciplinary information pools will be interfaced with one another to provide scientists with an integrated, holistic set of tools to address fundamental questions pertaining to legume biology. The public interface to the MGI database can be accessed at http://www.ncgr.org/research/mgi.

Analysis of RNA-protein interactions by flow cytometry

Flow cytometry, in combination with advances in bead coding technologies, is maturing as a powerful high-throughput approach for analyzing molecular interactions. Applications of this technology include antibody assays and single nucleotide polymorphism mapping. This review describes the recent development of a microbead flow cytometric approach to analyze RNA-protein interactions and discusses emerging bead coding strategies that together will allow genome-wide identification of RNA-protein complexes. The microbead flow cytometric approach is flexible and provides new opportunities for functional genomic studies and small-molecule screening.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Systematic expression profiling of the mouse transcriptome using RIKEN cDNA mircoarrays

The number of known mRNA transcripts in the mouse has been greatly expanded by the RIKEN Mouse Gene Encyclopedia project. Validation of their reproducible expression in a tissue is an important contribution to the study of functional genomics. In this report, we determine the expression profile of 57,931 clones on 20 mouse tissues using cDNA microarrays. Of these 57,931 clones, 22,928 clones correspond to the FANTOM2 clone set. The set represents 20,234 transcriptional units (TUs) out of 33,409 TUs in the FANTOM2 set. We identified 7206 separate clones that satisfied stringent criteria for tissue-specific expression. Gene Ontology terms were assigned for these 7206 clones, and the proportion of 'molecular function' ontology for each tissue-specific clone was examined. These data will provide insights into the function of each tissue. Tissue-specific gene expression profiles obtained using our cDNA microarrays were also compared with the data extracted from the GNF Expression Atlas based on Affymetrix microarrays. One major outcome of the RIKEN transcriptome analysis is the identification of numerous nonprotein-coding mRNAs. The expression profile was also used to obtain evidence of expression for putative noncoding RNAs. In addition, 1926 clones (70%) of 2768 clones that were categorized as unknown EST, and 1969 (58%) clones of 3388 clones that were categorized as unclassifiable were also shown to be reproducibly expressed.

Expression profiling in alcoholism research

This article represents the proceedings of a symposium at the 2004 International Society for Biomedical Research on Alcoholism in Mannheim, Germany, organized and co-chaired by Susan E. Bergeson and Wolfgang Sommer. The presentations and presenter were (1) Gene Expression in Brains of AlcoholPreferring and Non-Preferring Rats, by Howard J. Edenberg (2) Candidate Treatment Targets for Alcoholism: Leads from Functional Genomics Approaches, by Wolfgang Sommer (3) Microarray Analysis of Acute and Chronic Alcohol Response in Brain, by Susan E. Bergeson (4) On the Integration of QTL and Gene Expression Analysis, by Robert J. Hitzemann (5) Microarray and Proteomic Analysis of the Human Alcoholic Brain, by Peter R. Dodd.

Streptolysin-O reversible permeabilisation is an effective method to transfect siRNAs into myeloma cells

RNA interference (RNAi) has been shown to be a valuable tool to specifically target gene expression in a number of organisms becoming an indispensable weapon in the arsenal in functional genomics. In this study, we demonstrate that streptolysin-O (SLO) reversible permeabilisation is an efficient method to deliver small interfering RNAs (siRNAs) to hard-to-transfect human myeloma cell lines. We used published, pre-validated siRNAs for ERK2 and non-silencing siRNA control. We transfected siRNAs into human myeloma cell lines using SLO reversible permeabilisation method. Flow cytometry and western blot analysis were performed to assess the effect of SLO on transfection efficiency and ERK2 knockdown. These experiments demonstrate that SLO reversible permeabilisation method is an efficient and easy-to-use method to deliver siRNAs into human myeloma cell lines. Optimised SLO permeabilisation method showed to transfect >80% of JIM-3, H929, RPM18226 and U266 cells, with minimal effect on cell viability (<10%) and cell cycle. Equally important, SLO permeabilisation induced a substantial knockdown of ERK2 at the protein level. These studies demonstrate that reversible SLO permeabilisation can successfully be applied to hard-to-transfect human myeloma cell lines to effectively silence genes. (C) 2008 Published by Elsevier B.V.

Stimulating Neoblast-Like Cell Proliferation in Juvenile Fasciola hepatica Supports Growth and Progression towards the Adult Phenotype In Vitro

Fascioliasis (or fasciolosis) is a socioeconomically important parasitic disease caused by liver flukes of the genus Fasciola. Flukicide resistance has exposed the need for new drugs and/or a vaccine for liver fluke control. A rapidly improving 'molecular toolbox' for liver fluke encompasses quality genomic/transcriptomic datasets and an RNA interference platform that facilitates functional genomics approaches to drug/vaccine target validation. The exploitation of these resources is undermined by the absence of effective culture/maintenance systems that would support in vitro studies on juvenile fluke development/biology. Here we report markedly improved in vitro maintenance methods for Fasciola hepatica that achieved 65% survival of juvenile fluke after 6 months in standard cell culture medium supplemented with 50% chicken serum. We discovered that this long-term maintenance was dependent upon fluke growth, which was supported by increased proliferation of cells resembling the "neoblast" stem cells described in other flatworms. Growth led to dramatic morphological changes in juveniles, including the development of the digestive tract, reproductive organs and the tegument, towards more adult-like forms. The inhibition of DNA synthesis prevented neoblast-like cell proliferation and inhibited growth/development. Supporting our assertion that we have triggered the development of juveniles towards adult-like fluke, mass spectrometric analyses showed that growing fluke have an excretory/secretory protein profile that is distinct from that of newly-excysted juveniles and more closely resembles that of ex vivo immature and adult fluke. Further, in vitro maintained fluke displayed a transition in their movement from the probing behaviour associated with migrating stage worms to a slower wave-like motility seen in adults. Our ability to stimulate neoblast-like cell proliferation and growth in F. hepatica underpins the first simple platform for their long-term in vitro study, complementing the recent expansion in liver fluke resources and facilitating in vitro target validation studies of the developmental biology of liver fluke.