TP53 mutation prevalence in normal airway epithelium as a biomarker for lung cancer risk.

BACKGROUND: There is a need for biomarkers that improve accuracy compared with current demographic risk indices to detect individuals at the highest lung cancer risk. Improved risk determination will enable more effective lung cancer screening and better stratification of lung nodules into high or low-risk category. We previously reported discovery of a biomarker for lung cancer risk characterized by increased prevalence of TP53 somatic mutations in airway epithelial cells (AEC). Here we present results from a validation study in an independent retrospective case-control cohort. METHODS: Targeted next generation sequencing was used to identify mutations within three TP53 exons spanning 193 base pairs in AEC genomic DNA. RESULTS: TP53 mutation prevalence was associated with cancer status (P < 0.001). The lung cancer detection receiver operator characteristic (ROC) area under the curve (AUC) for the TP53 biomarker was 0.845 (95% confidence limits 0.749-0.942). In contrast, TP53 mutation prevalence was not significantly associated with age or smoking pack-years. The combination of TP53 mutation prevalence with PLCOM2012 risk score had an ROC AUC of 0.916 (0.846-0.986) and this was significantly higher than that for either factor alone (P < 0.03). CONCLUSIONS: These results support the validity of the TP53 mutation prevalence biomarker and justify taking additional steps to assess this biomarker in AEC specimens from a prospective cohort and in matched nasal brushing specimens as a potential non-invasive surrogate specimen.

Technical advance in targeted NGS analysis enables identification of lung cancer risk-associated low frequency TP53, PIK3CA, and BRAF mutations in airway epithelial cells.

BACKGROUND: Standardized Nucleic Acid Quantification for SEQuencing (SNAQ-SEQ) is a novel method that utilizes synthetic DNA internal standards spiked into each sample prior to next generation sequencing (NGS) library preparation. This method was applied to analysis of normal appearing airway epithelial cells (AEC) obtained by bronchoscopy in an effort to define a somatic mutation field effect associated with lung cancer risk. There is a need for biomarkers that reliably detect those at highest lung cancer risk, thereby enabling more effective screening by annual low dose CT. The purpose of this study was to test the hypothesis that lung cancer risk is characterized by increased prevalence of low variant allele frequency (VAF) somatic mutations in lung cancer driver genes in AEC. METHODS: Synthetic DNA internal standards (IS) were prepared for 11 lung cancer driver genes and mixed with each AEC genomic (g) DNA specimen prior to competitive multiplex PCR amplicon NGS library preparation. A custom Perl script was developed to separate IS reads and respective specimen gDNA reads from each target into separate files for parallel variant frequency analysis. This approach identified nucleotide-specific sequencing error and enabled reliable detection of specimen mutations with VAF as low as 5 × 10- 4 (0.05%). This method was applied in a retrospective case-control study of AEC specimens collected by bronchoscopic brush biopsy from the normal airways of 19 subjects, including eleven lung cancer cases and eight non-cancer controls, and the association of lung cancer risk with AEC driver gene mutations was tested. RESULTS: TP53 mutations with 0.05-1.0% VAF were more prevalent (p < 0.05) and also enriched for tobacco smoke and age-associated mutation signatures in normal AEC from lung cancer cases compared to non-cancer controls matched for smoking and age. Further, PIK3CA and BRAF mutations in this VAF range were identified in AEC from cases but not controls. CONCLUSIONS: Application of SNAQ-SEQ to measure mutations in the 0.05-1.0% VAF range enabled identification of an AEC somatic mutation field of injury associated with lung cancer risk. A biomarker comprising TP53, PIK3CA, and BRAF somatic mutations may better stratify individuals for optimal lung cancer screening and prevention outcomes.

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.

A lung cancer risk classifier comprising genome maintenance genes measured in normal bronchial epithelial cells.

BACKGROUND: Annual low dose CT (LDCT) screening of individuals at high demographic risk reduces lung cancer mortality by more than 20%. However, subjects selected for screening based on demographic criteria typically have less than a 10% lifetime risk for lung cancer. Thus, there is need for a biomarker that better stratifies subjects for LDCT screening. Toward this goal, we previously reported a lung cancer risk test (LCRT) biomarker comprising 14 genome-maintenance (GM) pathway genes measured in normal bronchial epithelial cells (NBEC) that accurately classified cancer (CA) from non-cancer (NC) subjects. The primary goal of the studies reported here was to optimize the LCRT biomarker for high specificity and ease of clinical implementation. METHODS: Targeted competitive multiplex PCR amplicon libraries were prepared for next generation sequencing (NGS) analysis of transcript abundance at 68 sites among 33 GM target genes in NBEC specimens collected from a retrospective cohort of 120 subjects, including 61 CA cases and 59 NC controls. Genes were selected for analysis based on contribution to the previously reported LCRT biomarker and/or prior evidence for association with lung cancer risk. Linear discriminant analysis was used to identify the most accurate classifier suitable to stratify subjects for screening. RESULTS: After cross-validation, a model comprising expression values from 12 genes (CDKN1A, E2F1, ERCC1, ERCC4, ERCC5, GPX1, GSTP1, KEAP1, RB1, TP53, TP63, and XRCC1) and demographic factors age, gender, and pack-years smoking, had Receiver Operator Characteristic area under the curve (ROC AUC) of 0.975 (95% CI: 0.96-0.99). The overall classification accuracy was 93% (95% CI 88%-98%) with sensitivity 93.1%, specificity 92.9%, positive predictive value 93.1% and negative predictive value 93%. The ROC AUC for this classifier was significantly better (p < 0.0001) than the best model comprising demographic features alone. CONCLUSIONS: The LCRT biomarker reported here displayed high accuracy and ease of implementation on a high throughput, quality-controlled targeted NGS platform. As such, it is optimized for clinical validation in specimens from the ongoing LCRT blinded prospective cohort study. Following validation, the biomarker is expected to have clinical utility by better stratifying subjects for annual lung cancer screening compared to current demographic criteria alone.

Haplotype and diplotype analyses of variation in ERCC5 transcription cis-regulation in normal bronchial epithelial cells.

Excision repair cross-complementation group 5 (ERCC5) gene plays an important role in nucleotide excision repair, and dysregulation of ERCC5 is associated with increased lung cancer risk. Haplotype and diplotype analyses were conducted in normal bronchial epithelial cells (NBEC) to better understand mechanisms responsible for interindividual variation in transcript abundance regulation of ERCC5 We determined genotypes at putative ERCC5 cis-regulatory SNPs (cis-rSNP) rs751402 and rs2296147, and marker SNPs rs1047768 and rs17655. ERCC5 allele-specific transcript abundance was assessed by a recently developed targeted sequencing method. Syntenic relationships among alleles at rs751402, rs2296147, and rs1047768 were assessed by allele-specific PCR followed by Sanger sequencing. We then assessed association of ERCC5 allele-specific expression at rs1047768 with haplotype and diplotype structure at cis-rSNPs rs751402 and rs2296147. Genotype analysis revealed significantly (P < 0.005) higher interindividual variation in allelic ratios in cDNA samples relative to matched gDNA samples at both rs1047768 and rs17655. By diplotype analysis, mean expression was higher at the rs1047768 alleles syntenic with rs2296147 T allele compared with rs2296147 C allele. Furthermore, mean expression was lower at rs17655 C allele, which is syntenic with G allele at a linked SNP rs873601 (D' = 0.95). These data support the conclusions that in NBEC, T allele at SNP rs2296147 upregulates ERCC5, variation at rs751402 does not alter ERCC5 regulation, and that C allele at SNP rs17655 downregulates ERCC5 Variation in ERCC5 transcript abundance associated with allelic variation at these SNPs could result in variation in NER function in NBEC and lung cancer risk.

Microcystin-LR aerosol induces inflammatory responses in healthy human primary airway epithelium.

Harmful algal blooms plague bodies of freshwater globally. These blooms are often composed of outgrowths of cyanobacteria capable of producing the heptapeptide Microcystin-LR (MC-LR) which is a well-known hepatotoxin. Recently, MC-LR has been detected in aerosols generated from lake water. However, the risk for human health effects due to MC-LR inhalation exposure have not been extensively investigated. In this study, we exposed a fully differentiated 3D human airway epithelium derived from 14 healthy donors to MC-LR-containing aerosol once a day for 3 days. Concentrations of MC-LR ranged from 100 pM to 1 µM. Although there were little to no detrimental alterations in measures of the airway epithelial function (i.e. cell survival, tissue integrity, mucociliary clearance, or cilia beating frequency), a distinct shift in the transcriptional activity was found. Genes related to inflammation were found to be upregulated such as C-C motif chemokine 5 (CCL5; log2FC = 0.57, p = 0.03) and C-C chemokine receptor type 7 (CCR7; log2FC = 0.84, p = 0.03). Functionally, conditioned media from MC-LR exposed airway epithelium was also found to have significant chemo-attractive properties for primary human neutrophils. Additionally, increases were found in the concentration of secreted chemokine proteins in the conditioned media such as CCL1 (log2FC = 5.07, p = 0.0001) and CCL5 (log2FC = 1.02, p = 0.046). These results suggest that MC-LR exposure to the human airway epithelium is capable of inducing an inflammatory response that may potentiate acute or chronic disease.

Deep oncopanel sequencing reveals within block position-dependent quality degradation in FFPE processed samples.

BACKGROUND: Clinical laboratories routinely use formalin-fixed paraffin-embedded (FFPE) tissue or cell block cytology samples in oncology panel sequencing to identify mutations that can predict patient response to targeted therapy. To understand the technical error due to FFPE processing, a robustly characterized diploid cell line was used to create FFPE samples with four different pre-tissue processing formalin fixation times. A total of 96 FFPE sections were then distributed to different laboratories for targeted sequencing analysis by four oncopanels, and variants resulting from technical error were identified. RESULTS: Tissue sections that fail more frequently show low cellularity, lower than recommended library preparation DNA input, or target sequencing depth. Importantly, sections from block surfaces are more likely to show FFPE-specific errors, akin to "edge effects" seen in histology, while the inner samples display no quality degradation related to fixation time. CONCLUSIONS: To assure reliable results, we recommend avoiding the block surface portion and restricting mutation detection to genomic regions of high confidence.

Advancing NGS quality control to enable measurement of actionable mutations in circulating tumor DNA.

The primary objective of the FDA-led Sequencing and Quality Control Phase 2 (SEQC2) project is to develop standard analysis protocols and quality control metrics for use in DNA testing to enhance scientific research and precision medicine. This study reports a targeted next-generation sequencing (NGS) method that will enable more accurate detection of actionable mutations in circulating tumor DNA (ctDNA) clinical specimens. To accomplish this, a synthetic internal standard spike-in was designed for each actionable mutation target, suitable for use in NGS following hybrid capture enrichment and unique molecular index (UMI) or non-UMI library preparation. When mixed with contrived ctDNA reference samples, internal standards enabled calculation of technical error rate, limit of blank, and limit of detection for each variant at each nucleotide position in each sample. True-positive mutations with variant allele fraction too low for detection by current practice were detected with this method, thereby increasing sensitivity.

Cis-acting genetic variation at an E2F1/YY1 response site and putative p53 site is associated with altered allele-specific expression of ERCC5 (XPG) transcript in normal human bronchial epithelium.

ERCC5 (XPG) is a key component of the nucleotide excision DNA repair pathway. In two recent case-control studies, we determined that ERCC5 transcript expression pattern in grossly normal human bronchial epithelial cells (NBEC) was different in individuals diagnosed with lung cancer compared with non-lung cancer controls. In this study, we tested the hypothesis that variation at cis-acting sites contributed to observed variation in ERCC5 transcript expression in NBEC. Allele-specific expression (ASE) was measured at transcribed polymorphic site rs1047768 in exon 2 of ERCC5 in NBEC complementary DNA (cDNA) of 22 individuals using allele-specific competitive polymerase chain reaction. ASE at rs1047768 was then assessed for association with allelotype at polymorphic sites rs751402 (E2F1 and YY1 recognition and response site) and rs2296147 (putative P53 recognition site) in the proximal promoter and 5' untranslated region, respectively, of ERCC5. Interindividual variation in recombination between rs751402, rs2296147 and rs1047768 in poly-heterozygotes was controlled for by allele-specific sequencing. Measured rs1047768 T:C allelic ratio was (i) significantly higher in NBEC cDNA compared with genomic DNA controls (P < 0.001) among samples heterozygous at both rs751402 and rs2296147; (ii) less high (P = 0.02) for samples homozygous at rs751402 but heterozygous at rs2296147 and (iii) not significantly different (P = 0.18) for doubly homozygous individuals. Here, we demonstrate that rs751402 A allele and rs2296147 T allele are associated with higher ASE of ERCC5 T allele transcript at rs1047768 in NBEC.

Platelet-T cell aggregates in lung cancer patients: Implications for thrombosis.

Platelet-leukocyte aggregates (PLAs) are associated with increased thrombosis risk. The influence of PLA formation is especially important for cancer patients, since thrombosis accounts for approximately 10% of cancer-associated deaths. Our objective was to characterize and quantify PLAs in whole blood samples from lung cancer patients compared to healthy volunteers with the intent to analyze PLA formation in the context of lung cancer-associated thrombosis. Consenting lung cancer patients (57) and healthy volunteers (56) were enrolled at the Dana Cancer Center at the University of Toledo Health Science Campus. Peripheral blood samples were analyzed by flow cytometry. Patient medical history was reviewed through electronic medical records. Most importantly, we found lung cancer patients to have higher percentages of platelet-T cell aggregates (PTCAs) than healthy volunteers among both CD4+ T lymphocyte and CD8+ T lymphocyte populations. Our findings demonstrate that characterization of PTCAs may have clinical utility in differentiating lung cancer patients from healthy volunteers and stratifying lung cancer patients by history of thrombosis.

Circulating Lactonase Activity but Not Protein Level of PON-1 Predicts Adverse Outcomes in Subjects with Chronic Kidney Disease.

The burden of cardiovascular disease and death in chronic kidney disease (CKD) outpaces that of the other diseases and is not adequately described by traditional risk factors alone. Diminished activity of paraoxonase (PON)-1 is associated with increased oxidant stress, a common feature underlying the pathogenesis of CKD. We aimed to assess the prognostic value of circulating PON-1 protein and PON lactonase activity on adverse clinical outcomes across various stages and etiologies of CKD. Circulating PON-1 protein levels and PON lactonase activity were measured simultaneously in patients with CKD as well as a cohort of apparently healthy non-CKD subjects. Both circulating PON-1 protein levels and PON lactonase activity were significantly lower in CKD patients compared to the non-CKD subjects. Similarly, across all stages of CKD, circulating PON-1 protein and PON lactonase activity were significantly lower in patients with CKD compared to the non-CKD controls. Circulating PON lactonase activity, but not protein levels, predicted future adverse clinical outcomes, even after adjustment for traditional risk factors. The combination of lower circulating protein levels and higher activity within the CKD subjects were associated with the best survival outcomes. These findings demonstrate that diminished circulating PON lactonase activity, but not protein levels, predicts higher risk of future adverse clinical outcomes in patients with CKD.

Chronic Low Dose Oral Exposure to Microcystin-LR Exacerbates Hepatic Injury in a Murine Model of Non-Alcoholic Fatty Liver Disease.

Microcystins are potent hepatotoxins that have become a global health concern in recent years. Their actions in at-risk populations with pre-existing liver disease is unknown. We tested the hypothesis that the No Observed Adverse Effect Level (NOAEL) of Microcystin-LR (MC-LR) established in healthy mice would cause exacerbation of hepatic injury in a murine model (Leprdb/J) of Non-alcoholic Fatty Liver Disease (NAFLD). Ten-week-old male Leprdb/J mice were gavaged with 50 µg/kg, 100 µg/kg MC-LR or vehicle every 48 h for 4 weeks (n = 15-17 mice/group). Early mortality was observed in both the 50 µg/kg (1/17, 6%), and 100 µg/kg (3/17, 18%) MC-LR exposed mice. MC-LR exposure resulted in significant increases in circulating alkaline phosphatase levels, and histopathological markers of hepatic injury as well as significant upregulation of genes associated with hepatotoxicity, necrosis, nongenotoxic hepatocarcinogenicity and oxidative stress response. In addition, we observed exposure dependent changes in protein phosphorylation sites in pathways involved in inflammation, immune function, and response to oxidative stress. These results demonstrate that exposure to MC-LR at levels that are below the NOAEL established in healthy animals results in significant exacerbation of hepatic injury that is accompanied by genetic and phosphoproteomic dysregulation in key signaling pathways in the livers of NAFLD mice.

Stable low-level expression of p21WAF1/CIP1 in A549 human bronchogenic carcinoma cell line-derived clones down-regulates E2F1 mRNA and restores cell proliferation control.

BACKGROUND: Deregulated cell cycle progression and loss of proliferation control are key properties of malignant cells. In previous studies, an interactive transcript abundance index (ITAI) comprising three cell cycle control genes, [MYC x E2F1]/p21 accurately distinguished normal from malignant bronchial epithelial cells (BEC), using a cut-off threshold of 7,000. This cut-off is represented by a line with a slope of 7,000 on a bivariate plot of p21 versus [MYC x E2F1], with malignant BEC above the line and normal BEC below the line. This study was an effort to better quantify, at the transcript abundance level, the difference between normal and malignant BEC. The hypothesis was tested that experimental elevation of p21 in a malignant BEC line would decrease the value of the [MYC x E2F1]/p21 ITAI to a level below this line, resulting in loss of immortality and limited cell population doubling capacity. In order to test the hypothesis, a p21 expression vector was transfected into the A549 human bronchogenic carcinoma cell line, which has low constitutive p21 TA expression relative to normal BEC. RESULTS: Following transfection of p21, four A549/p21 clones with stable two-fold up-regulated p21 expression were isolated and expanded. For each clone, the increase in p21 transcript abundance (TA) was associated with increased total p21 protein level, more than 5-fold reduction in E2F1 TA, and 10-fold reduction in the [MYC x E2F1]/p21 ITAI to a value below the cut-off threshold. These changes in regulation of cell cycle control genes were associated with restoration of cell proliferation control. Specifically, each transfectant was capable of only 15 population doublings compared with unlimited population doublings for parental A549. This change was associated with an approximate 2-fold increase in population doubling time to 38.4 hours (from 22.3 hrs), resumption of contact-inhibition, and reduced dividing cell fraction as measured by flow cytometric DNA analysis. CONCLUSION: These results, likely due to increased p21-mediated down-regulation of E2F1 TA at the G1/S phase transition, are consistent with our hypothesis. Specifically, they provide experimental confirmation that a line with slope of 7,000 on the p21 versus [MYC x E2F1] bivariate plot quantifies the difference between normal and malignant BEC at the level of transcript abundance.

ABCC5, ERCC2, XPA and XRCC1 transcript abundance levels correlate with cisplatin chemoresistance in non-small cell lung cancer cell lines.

BACKGROUND: Although 40-50% of non-small cell lung cancer (NSCLC) tumors respond to cisplatin chemotherapy, there currently is no way to prospectively identify potential responders. The purpose of this study was to determine whether transcript abundance (TA) levels of twelve selected DNA repair or multi-drug resistance genes (LIG1, ERCC2, ERCC3, DDIT3, ABCC1, ABCC4, ABCC5, ABCC10, GTF2H2, XPA, XPC and XRCC1) were associated with cisplatin chemoresistance and could therefore contribute to the development of a predictive marker. Standardized RT (StaRT)-PCR, was employed to assess these genes in a set of NSCLC cell lines with a previously published range of sensitivity to cisplatin. Data were obtained in the form of target gene molecules relative to 106 beta-actin (ACTB) molecules. To cancel the effect of ACTB variation among the different cell lines individual gene expression values were incorporated into ratios of one gene to another. Each two-gene ratio was compared as a single variable to chemoresistance for each of eight NSCLC cell lines using multiple regression. In an effort to validate these results, six additional lines then were evaluated. RESULTS: Following validation, single variable models best correlated with chemoresistance (p < 0.001), were ERCC2/XPC, ABCC5/GTF2H2, ERCC2/GTF2H2, XPA/XPC and XRCC1/XPC. All single variable models were examined hierarchically to achieve two variable models. The two variable model with the highest correlation was (ABCC5/GTF2H2, ERCC2/GTF2H2) with an R2 value of 0.96 (p < 0.001). CONCLUSION: These results provide markers suitable for assessment of small fine needle aspirate biopsies in an effort to prospectively identify cisplatin resistant tumors.

Variation in transcriptional regulation of cyclin dependent kinase inhibitor p21waf1/cip1 among human bronchogenic carcinomas.

BACKGROUND: Cell proliferation control depends in part on the carefully ordered regulation of transcription factors. The p53 homolog p73, contributes to this control by directly upregulating the cyclin dependent kinase inhibitor, p21waf1/cip1. E2F1, an inducer of cell proliferation, directly upregulates p73 and in some systems upregulates p21 directly. Because of its central role in controlling cell proliferation, upregulation of p21 has been explored as a modality for treating bronchogenic carcinoma (BC). Improved understanding of p21 transcriptional regulation will facilitate identification of BC tissues that are responsive to p21-directed therapies. Toward this goal, we investigated the role that E2F1 and p73 each play in the transcriptional regulation of p21. RESULTS: Among BC samples (N = 21) p21 transcript abundance (TA) levels varied over two orders of magnitude with values ranging from 400 to 120,000 (in units of molecules/106 molecules beta-actin). The p21 values in many BC were high compared to those observed in normal bronchial epithelial cells (BEC) (N = 18). Among all BC samples, there was no correlation between E2F1 and p21 TA but there was positive correlation between E2F1 and p73alpha (p < 0.001) TA. Among BC cell lines with inactivated p53 and wild type p73 (N = 7) there was positive correlation between p73alpha and p21 TA (p < 0.05). Additionally, in a BC cell line in which both p53 and p73 were inactivated (H1155), E2F1 TA level was high (50,000), but p21 TA level was low (470). Transiently expressed exogenous p73alpha in the BC cell line Calu-1, was associated with a significant (p < 0.05) 90% increase in p21 TA and a 20% reduction in E2F1 TA. siRNA mediated reduction of p73 TA in the N417 BC cell line was associated with a significant reduction in p21 TA level (p < 0.01). CONCLUSION: p21 TA levels vary considerably among BC patients which may be attributable to 1) genetic alterations in Rb and p53 and 2) variation in TA levels of upstream transcription factors E2F1 and p73. Here we provide evidence that p73 upregulates p21 TA in BC tissues and upregulated p21 TA may result from E2F1 upregulation of p73 but not from E2F1 directly.

CEBPG transcription factor correlates with antioxidant and DNA repair genes in normal bronchial epithelial cells but not in individuals with bronchogenic carcinoma.

BACKGROUND: Cigarette smoking is the primary cause of bronchogenic carcinoma (BC), yet only 10-15% of heavy smokers develop BC and it is likely that this variation in risk is, in part, genetically determined. We previously reported a set of antioxidant genes for which transcript abundance was lower in normal bronchial epithelial cells (NBEC) of BC individuals compared to non-BC individuals. In unpublished studies of the same NBEC samples, transcript abundance values for several DNA repair genes were correlated with these antioxidant genes. From these data, we hypothesized that antioxidant and DNA repair genes are co-regulated by one or more transcription factors and that inter-individual variation in expression and/or function of one or more of these transcription factors is responsible for inter-individual variation in risk for BC. METHODS: The putative transcription factor recognition sites common to six of the antioxidant genes were identified through in silico DNA sequence analysis. The transcript abundance values of these transcription factors (n = 6) and an expanded group of antioxidant and DNA repair genes (n = 16) were measured simultaneously by quantitative PCR in NBEC of 24 non-BC and 25 BC individuals. RESULTS: CEBPG transcription factor was significantly (p < 0.01) correlated with eight of the antioxidant or DNA repair genes in non-BC individuals but not in BC individuals. In BC individuals the correlation with CEBPG was significantly (p < 0.01) lower than that of non-BC individuals for four of the genes (XRCC1, ERCC5, GSTP1, and SOD1) and the difference was nearly significant for GPX1. The only other transcription factor correlated with any of these five target genes in non-BC individuals was E2F1. E2F1 was correlated with GSTP1 among non-BC individuals, but in contrast to CEBPG, there was no significant difference in this correlation in non-BC individuals compared to BC individuals. CONCLUSION: We conclude that CEBPG is the transcription factor primarily responsible for regulating transcription of key antioxidant and DNA repair genes in non-BC individuals. Further, we conclude that the heavy smokers selected for development of BC are those who have sub-optimal regulation of antioxidant and DNA repair genes by CEBPG.

Control for stochastic sampling variation and qualitative sequencing error in next generation sequencing.

BACKGROUND: Clinical implementation of Next-Generation Sequencing (NGS) is challenged by poor control for stochastic sampling, library preparation biases and qualitative sequencing error. To address these challenges we developed and tested two hypotheses. METHODS: Hypothesis 1: Analytical variation in quantification is predicted by stochastic sampling effects at input of a) amplifiable nucleic acid target molecules into the library preparation, b) amplicons from library into sequencer, or c) both. We derived equations using Monte Carlo simulation to predict assay coefficient of variation (CV) based on these three working models and tested them against NGS data from specimens with well characterized molecule inputs and sequence counts prepared using competitive multiplex-PCR amplicon-based NGS library preparation method comprising synthetic internal standards (IS). Hypothesis 2: Frequencies of technically-derived qualitative sequencing errors (i.e., base substitution, insertion and deletion) observed at each base position in each target native template (NT) are concordant with those observed in respective competitive synthetic IS present in the same reaction. We measured error frequencies at each base position within amplicons from each of 30 target NT, then tested whether they correspond to those within the 30 respective IS. RESULTS: For hypothesis 1, the Monte Carlo model derived from both sampling events best predicted CV and explained 74% of observed assay variance. For hypothesis 2, observed frequency and type of sequence variation at each base position within each IS was concordant with that observed in respective NTs (R2 = 0.93). CONCLUSION: In targeted NGS, synthetic competitive IS control for stochastic sampling at input of both target into library preparation and of target library product into sequencer, and control for qualitative errors generated during library preparation and sequencing. These controls enable accurate clinical diagnostic reporting of confidence limits and limit of detection for copy number measurement, and of frequency for each actionable mutation.

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.

Standardized RT-PCR and the standardized expression measurement center.

Standardized reverse transcriptase polymerase chain reaction (StaRT-PCR) is a modification of the competitive template (CT) RT method described by Gilliland et al. StaRT-PCR allows rapid, reproducible, standardized, quantitative measurement of data for many genes simultaneously. An internal standard CT is prepared for each gene, cloned to generate enough for 10(9) assays and CTs for up to 1,000 genes are mixed together. Each target gene is normalized to a reference gene to control for cDNA loaded in a standardized mixture of internal standards (SMIS) into the reaction. Each target gene and reference gene is measured relative to its respective internal standard within the SMIS. Because each target gene and reference gene is simultaneously measured relative to a known number of internal standard molecules in the SMIS, it is possible to report each gene expression measurement as a numerical value in units of target gene cDNA molecules/ 10(6) reference gene cDNA molecules. Calculation of data in this format allows for entry into a common databank, direct interexperimental comparison, and combination of values into interactive gene expression indices.

A multiplex two-color real-time PCR method for quality-controlled molecular diagnostic testing of FFPE samples.

BACKGROUND: Reverse transcription quantitative real-time PCR (RT-qPCR) tests support personalized cancer treatment through more clinically meaningful diagnosis. However, samples obtained through standard clinical pathology procedures are formalin-fixed, paraffin-embedded (FFPE) and yield small samples with low integrity RNA containing PCR interfering substances. RT-qPCR tests able to assess FFPE samples with quality control and inter-laboratory reproducibility are needed. METHODS: We developed an RT-qPCR method by which 1) each gene was measured relative to a known number of its respective competitive internal standard molecules to control for interfering substances, 2) two-color fluorometric hydrolysis probes enabled analysis on a real-time platform, 3) external standards controlled for variation in probe fluorescence intensity, and 4) pre-amplification maximized signal from FFPE RNA samples. Reagents were developed for four genes comprised by a previously reported lung cancer diagnostic test (LCDT) then subjected to analytical validation using synthetic native templates as test articles to assess linearity, signal-to-analyte response, lower detection threshold, imprecision and accuracy. Fitness of this method and these reagents for clinical testing was assessed in FFPE normal (N = 10) and malignant (N = 10) lung samples. RESULTS: Reagents for each of four genes, MYC, E2F1, CDKN1A and ACTB comprised by the LCDT had acceptable linearity (R(2)>0.99), signal-to-analyte response (slope 1.0 ± 0.05), lower detection threshold (<10 molecules) and imprecision (CV <20%). Poisson analysis confirmed accuracy of internal standard concentrations. Internal standards controlled for experimentally introduced interference, prevented false-negatives and enabled pre-amplification to increase signal without altering measured values. In the fitness for purpose testing of this two-color fluorometric LCDT using surgical FFPE samples, the diagnostic accuracy was 93% which was similar to that previously reported for analysis of fresh samples. CONCLUSIONS: This quality-controlled two-color fluorometric RT-qPCR approach will facilitate the development of reliable, robust RT-qPCR-based molecular diagnostic tests in FFPE clinical samples.