Yuriy O. Alekseyev

Characterizing the Impact of Smoking and Lung Cancer on the Airway Transcriptome Using RNA-Seq

Cigarette smoke creates a molecular field of injury in epithelial cells that line the respiratory tract. We hypothesized that transcriptome sequencing (RNA-Seq) will enhance our understanding of the field of molecular injury in response to tobacco smoke exposure and lung cancer pathogenesis by identifying gene expression differences not interrogated or accurately measured by microarrays. We sequenced the high-molecular-weight fraction of total RNA (>200 nt) from pooled bronchial airway epithelial cell brushings (n = 3 patients per pool) obtained during bronchoscopy from healthy never smoker (NS) and current smoker (S) volunteers and smokers with (C) and without (NC) lung cancer undergoing lung nodule resection surgery. RNA-Seq libraries were prepared using 2 distinct approaches, one capable of capturing non-polyadenylated RNA (the prototype NuGEN Ovation RNA-Seq protocol) and the other designed to measure only polyadenylated RNA (the standard Illumina mRNA-Seq protocol) followed by sequencing generating approximately 29 million 36 nt reads per pool and approximately 22 million 75 nt paired-end reads per pool, respectively. The NuGEN protocol captured additional transcripts not detected by the Illumina protocol at the expense of reduced coverage of polyadenylated transcripts, while longer read lengths and a paired-end sequencing strategy significantly improved the number of reads that could be aligned to the genome. The aligned reads derived from the two complementary protocols were used to define the compendium of genes expressed in the airway epithelium (n = 20,573 genes). Pathways related to the metabolism of xenobiotics by cytochrome P450, retinol metabolism, and oxidoreductase activity were enriched among genes differentially expressed in smokers, whereas chemokine signaling pathways, cytokine–cytokine receptor interactions, and cell adhesion molecules were enriched among genes differentially expressed in smokers with lung cancer. There was a significant correlation between the RNA-Seq gene expression data and Affymetrix microarray data generated from the same samples (P < 0.001); however, the RNA-Seq data detected additional smoking- and cancer-related transcripts whose expression was were either not interrogated by or was not found to be significantly altered when using microarrays, including smoking-related changes in the inflammatory genes S100A8 and S100A9 and cancer-related changes in MUC5AC and secretoglobin (SCGB3A1). Quantitative real-time PCR confirmed differential expression of select genes and non-coding RNAs within individual samples. These results demonstrate that transcriptome sequencing has the potential to provide new insights into the biology of the airway field of injury associated with smoking and lung cancer. The measurement of both coding and non-coding transcripts by RNA-Seq has the potential to help elucidate mechanisms of response to tobacco smoke and to identify additional biomarkers of lung cancer risk and novel targets for chemoprevention. Cancer Prev Res; 4(6); 803–17. ©2011 AACR.

Similarities and differences between smoking-related gene expression in nasal and bronchial epithelium.

Previous studies have shown that physiological responses to cigarette smoke can be detected via bronchial airway epithelium gene expression profiling and that heterogeneity in this gene expression response to smoking is associated with lung cancer. In this study, we sought to determine the similarity of the effects of tobacco smoke throughout the respiratory tract by determining patterns of smoking-related gene expression in paired nasal and bronchial epithelial brushings collected from 14 healthy nonsmokers and 13 healthy current smokers. Using whole genome expression arrays, we identified 119 genes whose expression was affected by smoking similarly in both bronchial and nasal epithelium, including genes related to detoxification, oxidative stress, and wound healing. While the vast majority of smoking-related gene expression changes occur in both bronchial and nasal epithelium, we also identified 27 genes whose expression was affected by smoking more dramatically in bronchial epithelium than nasal epithelium. Both common and site-specific smoking-related gene expression profiles were validated using independent microarray datasets. Differential expression of select genes was also confirmed by RT-PCR. That smoking induces largely similar gene expression changes in both nasal and bronchial epithelium suggests that the consequences of cigarette smoke exposure can be measured in tissues throughout the respiratory tract. Our findings suggest that nasal epithelial gene expression may serve as a relatively noninvasive surrogate to measure physiological responses to cigarette smoke and/or other inhaled exposures in large-scale epidemiological studies.

INACTIVATION OF AMPK ALTERS GENE EXPRESSION AND PROMOTES GROWTH OF PROSTATE CANCER CELLS

AMP-activated protein kinase (AMPK) serves as a fuel-sensing enzyme that is activated by binding of AMP and subsequent phophorylation by upstream kinases such as the tumor suppressor LKB1, when cells sense an increase in the ratio of AMP to ATP. Acute activation of AMPK stimulates fatty acid oxidation to generate more ATP and simultaneously inhibits ATP-consuming processes including fatty acid and protein syntheses, thereby preserving energy for acute cell-surviving program, whereas chronic activation leads to inhibition of cell growth. The goal of the present study is to explore the mechanisms by which AMPK regulates cell growth. Toward this end, we established stable cell lines by introducing a dominant-negative mutant of AMPK α1 subunit or its shRNA into the prostate cancer C4-2 cells and other cells, or wild type LKB1 into the lung adenocarcinoma A549 and breast MB-MDA-231 cancer cells, both of which lack functional LKB1. Our results showed that the inhibition of AMPK accelerated cell proliferation and promoted malignant behavior such as increased cell migration and anchorage-independent growth. This was associated with decreased G1 population, downregulation of p53 and p21, and upregulation of S6K, IGF-1 and IGF1R. Conversely, treatment of the C4-2 cells with 5-aminoimidazole-4-carboxamide 1-D-ribonucleoside (AICAR), a prototypical AMPK activator, caused opposite changes. In addition, our study using microarray and RT–PCR revealed that AMPK regulated gene expression involved in tumor cell growth and survival. Thus, our study provides novel insights into the mechanisms of AMPK action in cancer cells and presents AMPK as an ideal drug target for cancer therapy.

A gene expression signature of emphysema-related lung destruction and its reversal by the tripeptide GHK

BackgroundChronic obstructive pulmonary disease (COPD) is a heterogeneous disease consisting of emphysema, small airway obstruction, and/or chronic bronchitis that results in significant loss of lung function over time.MethodsIn order to gain insights into the molecular pathways underlying progression of emphysema and explore computational strategies for identifying COPD therapeutics, we profiled gene expression in lung tissue samples obtained from regions within the same lung with varying amounts of emphysematous destruction from smokers with COPD (8 regions × 8 lungs = 64 samples). Regional emphysema severity was quantified in each tissue sample using the mean linear intercept (Lm) between alveolar walls from micro-CT scans.ResultsWe identified 127 genes whose expression levels were significantly associated with regional emphysema severity while controlling for gene expression differences between individuals. Genes increasing in expression with increasing emphysematous destruction included those involved in inflammation, such as the B-cell receptor signaling pathway, while genes decreasing in expression were enriched in tissue repair processes, including the transforming growth factor beta (TGFβ) pathway, actin organization, and integrin signaling. We found concordant differential expression of these emphysema severity-associated genes in four cross-sectional studies of COPD. Using the Connectivity Map, we identified GHK as a compound that can reverse the gene-expression signature associated with emphysematous destruction and induce expression patterns consistent with TGFβ pathway activation. Treatment of human fibroblasts with GHK recapitulated TGFβ-induced gene-expression patterns, led to the organization of the actin cytoskeleton, and elevated the expression of integrin β1. Furthermore, addition of GHK or TGFβ restored collagen I contraction and remodeling by fibroblasts derived from COPD lungs compared to fibroblasts from former smokers without COPD.ConclusionsThese results demonstrate that gene-expression changes associated with regional emphysema severity within an individuals lung can provide insights into emphysema pathogenesis and identify novel therapeutic opportunities for this deadly disease. They also suggest the need for additional studies to examine the mechanisms by which TGFβ and GHK each reverse the gene-expression signature of emphysematous destruction and the effects of this reversal on disease progression.

A Dynamic Bronchial Airway Gene Expression Signature of Chronic Obstructive Pulmonary Disease and Lung Function Impairment

RATIONALE Molecular phenotyping of chronic obstructive pulmonary disease (COPD) has been impeded in part by the difficulty in obtaining lung tissue samples from individuals with impaired lung function. OBJECTIVES We sought to determine whether COPD-associated processes are reflected in gene expression profiles of bronchial airway epithelial cells obtained by bronchoscopy. METHODS Gene expression profiling of bronchial brushings obtained from 238 current and former smokers with and without COPD was performed using Affymetrix Human Gene 1.0 ST Arrays. MEASUREMENTS AND MAIN RESULTS We identified 98 genes whose expression levels were associated with COPD status, FEV1% predicted, and FEV1/FVC. In silico analysis identified activating transcription factor 4 (ATF4) as a potential transcriptional regulator of genes with COPD-associated airway expression, and ATF4 overexpression in airway epithelial cells in vitro recapitulates COPD-associated gene expression changes. Genes with COPD-associated expression in the bronchial airway epithelium had similarly altered expression profiles in prior studies performed on small-airway epithelium and lung parenchyma, suggesting that transcriptomic alterations in the bronchial airway epithelium reflect molecular events found at more distal sites of disease activity. Many of the airway COPD-associated gene expression changes revert toward baseline after therapy with the inhaled corticosteroid fluticasone in independent cohorts. CONCLUSIONS Our findings demonstrate a molecular field of injury throughout the bronchial airway of active and former smokers with COPD that may be driven in part by ATF4 and is modifiable with therapy. Bronchial airway epithelium may ultimately serve as a relatively accessible tissue in which to measure biomarkers of disease activity for guiding clinical management of COPD.

DNA polymerase V allows bypass of toxic guanine oxidation products in vivo

Reactive oxygen and nitrogen radicals produced during metabolic processes, such as respiration and inflammation, combine with DNA to form many lesions primarily at guanine sites. Understanding the roles of the polymerases responsible for the processing of these products to mutations could illuminate molecular mechanisms that correlate oxidative stress with cancer. Using M13 viral genomes engineered to contain single DNA lesions and Escherichia coli strains with specific polymerase (pol) knockouts, we show that pol V is required for efficient bypass of structurally diverse, highly mutagenic guanine oxidation products in vivo. We also find that pol IV participates in the bypass of two spiroiminodihydantoin lesions. Furthermore, we report that one lesion, 5-guanidino-4-nitroimidazole, is a substrate for multiple SOS polymerases, whereby pol II is necessary for error-free replication and pol V for error-prone replication past this lesion. The results spotlight a major role for pol V and minor roles for pol II and pol IV in the mechanism of guanine oxidation mutagenesis.

Airway gene expression in COPD is dynamic with inhaled corticosteroid treatment and reflects biological pathways associated with disease activity

Background A core feature of chronic obstructive pulmonary disease (COPD) is the accelerated decline in forced expiratory volume in one second (FEV1). The recent Groningen and Leiden Universities study of Corticosteroids in Obstructive Lung Disease (GLUCOLD) study suggested that particular phenotypes of COPD benefit from fluticasone±salmeterol by reducing the rate of FEV1 decline, yet the underlying mechanisms are unknown. Methods Whole-genome gene expression profiling using the Affymetrix Gene ST array (V.1.0) was performed on 221 bronchial biopsies available from 89 COPD patients at baseline and after 6 and 30 months of fluticasone±salmeterol and placebo treatment in GLUCOLD. Results Linear mixed effects modelling revealed that the expression of 138 genes decreased, whereas the expression of 140 genes significantly upregulated after both 6 and 30 months of treatment with fluticasone±salmeterol versus placebo. A more pronounced treatment-induced change in the expression of 50 and 55 of these 278 genes was associated with a lower rate of decline in FEV1 and Saint George Respiratory Questionnaire, respectively. Genes decreasing with treatment were involved in pathways related to cell cycle, oxidative phosphorylation, epithelial cell signalling, p53 signalling and T cell signalling. Genes increasing with treatment were involved in pathways related to focal adhesion, gap junction and extracellular matrix deposition. Finally, the fluticasone-induced gene expression changes were enriched among genes that change in the airway epithelium in smokers with versus without COPD in an independent data set. Conclusions The present study suggests that gene expression in biological pathways of COPD is dynamic with treatment and reflects disease activity. This study opens the gate to targeted and molecular phenotype-driven therapy of COPD.

Transcriptional responses of neonatal mouse lung to hyperoxia by Nrf2 status

UNLABELLED Hyperoxia exposure can inhibit alveolar growth in the neonatal lung through induction of p21/p53 pathways and is a risk factor for the development of bronchopulmonary dysplasia (BPD) in preterm infants. We previously found that activation of nuclear factor erythroid 2 p45-related factor (Nrf2) improved survival in neonatal mice exposed to hyperoxia likely due to increased expression of anti-oxidant response genes. It is not known however, whether hyperoxic induced Nrf2 activation attenuates the growth impairment caused by hyperoxia in neonatal lung. To determine if Nrf2 activation modulates cell cycle regulatory pathway genes associated with growth arrest we examined the gene expression in the lungs of Nrf2(-/-) and Nrf2(+/+) neonatal mice at one and 3days of hyperoxia exposure. METHODS Microarray analysis was performed in neonatal Nrf2(+/+) and Nrf2(-/-) lungs exposed to one and 3days of hyperoxia. Sulforaphane, an inducer of Nrf2 was given to timed pregnant mice to determine if in utero exposure attenuated p21 and IL-6 gene expression in wildtype neonatal mice exposed to hyperoxia. RESULTS Cell cycle regulatory genes were induced in Nrf2(-/-) lung at 1day of hyperoxia. At 3days of hyperoxia, induction of cell cycle regulatory genes was similar in Nrf2(+/+) and Nrf2(-/-) lungs, despite higher inflammatory gene expression in Nrf2(-/-) lung. CONCLUSION p21/p53 pathways gene expression was not attenuated by Nrf2 activation in neonatal lung. In utero SUL did not attenuate p21 expression in wildtype neonatal lung exposed to hyperoxia. These findings suggest that although Nrf2 activation induces expression of anti-oxidant genes, it does not attenuate alveolar growth arrest caused by exposure to hyperoxia.

Evaluation of commercially available RNA amplification kits for RNA sequencing using very low input amounts of total RNA

This article includes supplemental data. Please visit http://www.fasebj.org to obtain this information.Multiple recent publications on RNA sequencing (RNA-seq) have demonstrated the power of next-generation sequencing technologies in whole-transcriptome analysis. Vendor-specific protocols used for RNA library construction often require at least 100 ng total RNA. However, under certain conditions, much less RNA is available for library construction. In these cases, effective transcriptome profiling requires amplification of subnanogram amounts of RNA. Several commercial RNA amplification kits are available for amplification prior to library construction for next-generation sequencing, but these kits have not been comprehensively field evaluated for accuracy and performance of RNA-seq for picogram amounts of RNA. To address this, 4 types of amplification kits were tested with 3 different concentrations, from 5 ng to 50 pg, of a commercially available RNA. Kits were tested at multiple sites to assess reproducibility and ease of use. The human total reference RNA used was spiked with a control pool of RNA molecules in order to further evaluate quantitative recovery of input material. Additional control data sets were generated from libraries constructed following polyA selection or ribosomal depletion using established kits and protocols. cDNA was collected from the different sites, and libraries were synthesized at a single site using established protocols. Sequencing runs were carried out on the Illumina platform. Numerous metrics were compared among the kits and dilutions used. Overall, no single kit appeared to meet all the challenges of small input material. However, it is encouraging that excellent data can be recovered with even the 50 pg input total RNA.

Protein Kinase CK1αLS Promotes Vascular Cell Proliferation and Intimal Hyperplasia

Protein kinase CK1alpha regulates several fundamental cellular processes including proliferation and differentiation. Up to four forms of this kinase are expressed in vertebrates resulting from alternative splicing of exons; these exons encode either the L-insert located within the catalytic domain or the S-insert located at the C terminus of the protein. Whereas the L-insert is known to target the kinase to the nucleus, the functional significance of nuclear CK1alphaLS has been unclear. Here we demonstrate that selective L-insert-targeted short hairpin small interfering RNA-mediated knockdown of CK1alphaLS in human vascular endothelial cells and vascular smooth muscle cells impairs proliferation and abolishes hydrogen peroxide-stimulated proliferation of vascular smooth muscle cells, with the cells accumulating in G(0)/G(1). In addition, selective knockdown of CK1alphaLS in cultured human arteries inhibits vascular activation, preventing smooth muscle cell proliferation, intimal hyperplasia, and proteoglycan deposition. Knockdown of CK1alphaLS results in the harmonious down-regulation of its target substrate heterogeneous nuclear ribonucleoprotein C and results in the altered expression or alternative splicing of key genes involved in cellular activation including CXCR4, MMP3, CSF2, and SMURF1. Our results indicate that the nuclear form of CK1alpha in humans, CK1alphaLS, plays a critical role in vascular cell proliferation, cellular activation, and hydrogen peroxide-mediated mitogenic signal transduction.