RNAseq analysis of bronchial epithelial cells to identify COPD-associated genes and SNPs.

BACKGROUND: There is a need for more powerful methods to identify low-effect SNPs that contribute to hereditary COPD pathogenesis. We hypothesized that SNPs contributing to COPD risk through cis-regulatory effects are enriched in genes comprised by bronchial epithelial cell (BEC) expression patterns associated with COPD. METHODS: To test this hypothesis, normal BEC specimens were obtained by bronchoscopy from 60 subjects: 30 subjects with COPD defined by spirometry (FEV1/FVC < 0.7, FEV1% < 80%), and 30 non-COPD controls. Targeted next generation sequencing was used to measure total and allele-specific expression of 35 genes in genome maintenance (GM) genes pathways linked to COPD pathogenesis, including seven TP53 and CEBP transcription factor family members. Shrinkage linear discriminant analysis (SLDA) was used to identify COPD-classification models. COPD GWAS were queried for putative cis-regulatory SNPs in the targeted genes. RESULTS: On a network basis, TP53 and CEBP transcription factor pathway gene pair network connections, including key DNA repair gene ERCC5, were significantly different in COPD subjects (e.g., Wilcoxon rank sum test for closeness, p-value = 5.0E-11). ERCC5 SNP rs4150275 association with chronic bronchitis was identified in a set of Lung Health Study (LHS) COPD GWAS SNPs restricted to those in putative regulatory regions within the targeted genes, and this association was validated in the COPDgene non-hispanic white (NHW) GWAS. ERCC5 SNP rs4150275 is linked (D' = 1) to ERCC5 SNP rs17655 which displayed differential allelic expression (DAE) in BEC and is an expression quantitative trait locus (eQTL) in lung tissue (p = 3.2E-7). SNPs in linkage (D' = 1) with rs17655 were predicted to alter miRNA binding (rs873601). A classifier model that comprised gene features CAT, CEBPG, GPX1, KEAP1, TP73, and XPA had pooled 10-fold cross-validation receiver operator characteristic area under the curve of 75.4% (95% CI: 66.3%-89.3%). The prevalence of DAE was higher than expected (p = 0.0023) in the classifier genes. CONCLUSIONS: GM genes comprised by COPD-associated BEC expression patterns were enriched for SNPs with cis-regulatory function, including a putative cis-rSNP in ERCC5 that was associated with COPD risk. These findings support additional total and allele-specific expression analysis of gene pathways with high prior likelihood for involvement in COPD pathogenesis.

The c-myc x E2F-1/p21 interactive gene expression index augments cytomorphologic diagnosis of lung cancer in fine-needle aspirate specimens.

Morphological analysis of cytologic samples obtained by fine-needle aspirate (FNA) or bronchoscopy is an important method for diagnosing bronchogenic carcinoma. However, this approach has only about 65 to 80% diagnostic sensitivity. Based on previous studies, the c-myc x E2F-1/p21WAF1/CIP1 (p21 hereafter) gene expression index is highly sensitive and specific for distinguishing normal from malignant bronchial epithelial tissues. In an effort to improve sensitivity of diagnosing lung cancer in cytologic specimens, we used Standardized Reverse Transcriptase Polymerase Chain Reaction (StaRT-PCR) to measure the c-myc x E2F-1/p21 index in cDNA samples from 14 normal lung samples (6 normal lung parenchyma and 8 normal bronchial epithelial cell [NBEC] biopsies), and 16 FNA biopsies from 14 suspected tumors. Based on cytomorphologic criteria, 11 of the 14 suspected tumors were diagnosed as bronchogenic carcinoma and three specimens were non-diagnostic. Subsequent biopsy samples confirmed that the three non-diagnostic samples were derived from lung carcinomas. The index value for each bronchogenic carcinoma was above a cut-off value of 7000 and the index value of all but one normal sample was below 7000. Thus the c-myc x E2F-1/p21 index may augment cytomorphologic diagnosis of bronchogenic carcinoma biopsy samples, particularly those considered non-diagnostic by cytomorphologic criteria.