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.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer.

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

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.

Cross-oncopanel study reveals high sensitivity and accuracy with overall analytical performance depending on genomic regions.

BACKGROUND: Targeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing. RESULTS: All panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5-20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden. CONCLUSION: This comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.

Evaluating the analytical validity of circulating tumor DNA sequencing assays for precision oncology.

Circulating tumor DNA (ctDNA) sequencing is being rapidly adopted in precision oncology, but the accuracy, sensitivity and reproducibility of ctDNA assays is poorly understood. Here we report the findings of a multi-site, cross-platform evaluation of the analytical performance of five industry-leading ctDNA assays. We evaluated each stage of the ctDNA sequencing workflow with simulations, synthetic DNA spike-in experiments and proficiency testing on standardized, cell-line-derived reference samples. Above 0.5% variant allele frequency, ctDNA mutations were detected with high sensitivity, precision and reproducibility by all five assays, whereas, below this limit, detection became unreliable and varied widely between assays, especially when input material was limited. Missed mutations (false negatives) were more common than erroneous candidates (false positives), indicating that the reliable sampling of rare ctDNA fragments is the key challenge for ctDNA assays. This comprehensive evaluation of the analytical performance of ctDNA assays serves to inform best practice guidelines and provides a resource for precision oncology.

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.

Evaluation of DNA microarray results with quantitative gene expression platforms.

We have evaluated the performance characteristics of three quantitative gene expression technologies and correlated their expression measurements to those of five commercial microarray platforms, based on the MicroArray Quality Control (MAQC) data set. The limit of detection, assay range, precision, accuracy and fold-change correlations were assessed for 997 TaqMan Gene Expression Assays, 205 Standardized RT (Sta)RT-PCR assays and 244 QuantiGene assays. TaqMan is a registered trademark of Roche Molecular Systems, Inc. We observed high correlation between quantitative gene expression values and microarray platform results and found few discordant measurements among all platforms. The main cause of variability was differences in probe sequence and thus target location. A second source of variability was the limited and variable sensitivity of the different microarray platforms for detecting weakly expressed genes, which affected interplatform and intersite reproducibility of differentially expressed genes. From this analysis, we conclude that the MAQC microarray data set has been validated by alternative quantitative gene expression platforms thus supporting the use of microarray platforms for the quantitative characterization of gene expression.

Using RNA sample titrations to assess microarray platform performance and normalization techniques.

We have assessed the utility of RNA titration samples for evaluating microarray platform performance and the impact of different normalization methods on the results obtained. As part of the MicroArray Quality Control project, we investigated the performance of five commercial microarray platforms using two independent RNA samples and two titration mixtures of these samples. Focusing on 12,091 genes common across all platforms, we determined the ability of each platform to detect the correct titration response across the samples. Global deviations from the response predicted by the titration ratios were observed. These differences could be explained by variations in relative amounts of messenger RNA as a fraction of total RNA between the two independent samples. Overall, both the qualitative and quantitative correspondence across platforms was high. In summary, titration samples may be regarded as a valuable tool, not only for assessing microarray platform performance and different analysis methods, but also for determining some underlying biological features of the samples.

The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements.

Over the last decade, the introduction of microarray technology has had a profound impact on gene expression research. The publication of studies with dissimilar or altogether contradictory results, obtained using different microarray platforms to analyze identical RNA samples, has raised concerns about the reliability of this technology. The MicroArray Quality Control (MAQC) project was initiated to address these concerns, as well as other performance and data analysis issues. Expression data on four titration pools from two distinct reference RNA samples were generated at multiple test sites using a variety of microarray-based and alternative technology platforms. Here we describe the experimental design and probe mapping efforts behind the MAQC project. We show intraplatform consistency across test sites as well as a high level of interplatform concordance in terms of genes identified as differentially expressed. This study provides a resource that represents an important first step toward establishing a framework for the use of microarrays in clinical and regulatory settings.

The balance of reproducibility, sensitivity, and specificity of lists of differentially expressed genes in microarray studies.

BACKGROUND: Reproducibility is a fundamental requirement in scientific experiments. Some recent publications have claimed that microarrays are unreliable because lists of differentially expressed genes (DEGs) are not reproducible in similar experiments. Meanwhile, new statistical methods for identifying DEGs continue to appear in the scientific literature. The resultant variety of existing and emerging methods exacerbates confusion and continuing debate in the microarray community on the appropriate choice of methods for identifying reliable DEG lists. RESULTS: Using the data sets generated by the MicroArray Quality Control (MAQC) project, we investigated the impact on the reproducibility of DEG lists of a few widely used gene selection procedures. We present comprehensive results from inter-site comparisons using the same microarray platform, cross-platform comparisons using multiple microarray platforms, and comparisons between microarray results and those from TaqMan - the widely regarded "standard" gene expression platform. Our results demonstrate that (1) previously reported discordance between DEG lists could simply result from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion with a non-stringent P-value cutoff filtering, the DEG lists become much more reproducible, especially when fewer genes are selected as differentially expressed, as is the case in most microarray studies; and (3) the instability of short DEG lists solely based on P-value ranking is an expected mathematical consequence of the high variability of the t-values; the more stringent the P-value threshold, the less reproducible the DEG list is. These observations are also consistent with results from extensive simulation calculations. CONCLUSION: We recommend the use of FC-ranking plus a non-stringent P cutoff as a straightforward and baseline practice in order to generate more reproducible DEG lists. Specifically, the P-value cutoff should not be stringent (too small) and FC should be as large as possible. Our results provide practical guidance to choose the appropriate FC and P-value cutoffs when selecting a given number of DEGs. The FC criterion enhances reproducibility, whereas the P criterion balances sensitivity and specificity.

Fetal-juvenile origins of point mutations in the adult human tracheal-bronchial epithelium: absence of detectable effects of age, gender or smoking status.

Allele-specific mismatch amplification mutation assays (MAMA) of anatomically distinct sectors of the upper bronchial tracts of nine nonsmokers revealed many numerically dispersed clusters of the point mutations C742T, G746T, G747T of the TP53 gene, G35T of the KRAS gene and G508A of the HPRT1 gene. Assays of these five mutations in six smokers have yielded quantitatively similar results. One hundred and eighty four micro-anatomical sectors of 0.5-6x10(6) tracheal-bronchial epithelial cells represented en toto the equivalent of approximately 1.7 human smokers' bronchial trees to the fifth bifurcation. Statistically significant mutant copy numbers above the 95% upper confidence limits of historical background controls were found in 198 of 425 sector assays. No significant differences (P=0.1) for negative sector fractions, mutant fractions, distributions of mutant cluster size or anatomical positions were observed for smoking status, gender or age (38-76 year). Based on the modal cluster size of mitochondrial point mutants, the size of the adult bronchial epithelial maintenance turnover unit was estimated to be about 32 cells. When data from all 15 lungs were combined the log2 of nuclear mutant cluster size plotted against log2 of the number of clusters of a given cluster size displayed a slope of approximately 1.1 over a range of cluster sizes from approximately 2(6) to 2(15) mutant copies. A parsimonious interpretation of these nuclear and previously reported data for lung epithelial mitochondrial point mutant clusters is that they arose from mutations in stem cells at a high but constant rate per stem cell doubling during at least ten stem cell doublings of the later fetal-juvenile period. The upper and lower decile range of summed point mutant fractions among lungs was about 7.5-fold, suggesting an important source of stratification in the population with regard to risk of tumor initiation.

The interactive transcript abundance index [c-myc*p73alpha]/[p21*Bcl-2] correlates with baseline level of apoptosis and response to CPT-11 in human bronchogenic carcinoma cell lines.

Currently available cytotoxic chemotherapy is ineffective for treating bronchogenic carcinoma, and this is partly due to unpredictable inter-tumor variation in resistance. For example, tumors with inactivating p53 mutations or deletions are less likely to respond to certain chemotherapeutics. However, even if p53 is intact, a tumor may be unresponsive if defects in other p53 pathway genes compromise apoptosis. In an effort to identify biomarkers that better predict response to camptothecin, we investigated the association of CPT-11 (irinotecan)-induced cytotoxicity (IC50) with apoptosis or expression of genes upstream or downstream of p53 in cell lines that retain wild-type p53. CPT-11 response was greater in cell lines with higher baseline apoptosis (p<0.05). In addition, the interactive transcript abundance index (ITAI) [c-myc*p73alpha]/[p21*Bcl-2] was directly correlated with baseline apoptosis (p<0.01) and CPT-11 response (p<0.05). The ITAI was also correlated with CPT-11 response among cell lines derived from a variety of tissues that had inactivating p53 mutations or deletions, supporting its applicability for predicting response to camptothecins in other tissues regardless of p53 status.

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.