RNA sequencing in human HepG2 hepatocytes reveals PPAR-α mediates transcriptome responsiveness of bilirubin.

Bilirubin is a potent antioxidant that reduces inflammation and the accumulation of fat. There have been reports of gene responses to bilirubin, which was mostly attributed to its antioxidant function. Using RNA sequencing, we found that biliverdin, which is rapidly reduced to bilirubin, induced transcriptome responses in human HepG2 hepatocytes in a peroxisome proliferator-activated receptor (PPAR)-α-dependent fashion (398 genes with >2-fold change; false discovery rate P < 0.05). For comparison, a much narrower set of genes demonstrated differential expression when PPAR-α was suppressed via lentiviral shRNA knockdown (23 genes). Gene set enrichment analysis revealed the bilirubin-PPAR-α transcriptome mediates pathways for oxidation-reduction processes, mitochondrial function, response to nutrients, fatty acid oxidation, and lipid homeostasis. Together, these findings suggest that transcriptome responses from the generation of bilirubin are mostly PPAR-α dependent, and its antioxidant function regulates a smaller set of genes.

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.

Antioxidant Therapy Significantly Attenuates Hepatotoxicity following Low Dose Exposure to Microcystin-LR in a Murine Model of Diet-Induced Non-Alcoholic Fatty Liver Disease.

We have previously shown in a murine model of Non-alcoholic Fatty Liver Disease (NAFLD) that chronic, low-dose exposure to the Harmful Algal Bloom cyanotoxin microcystin-LR (MC-LR), resulted in significant hepatotoxicity including micro-vesicular lipid accumulation, impaired toxin metabolism as well as dysregulation of the key signaling pathways involved in inflammation, immune response and oxidative stress. On this background we hypothesized that augmentation of hepatic drug metabolism pathways with targeted antioxidant therapies would improve MC-LR metabolism and reduce hepatic injury in NAFLD mice exposed to MC-LR. We chose N-acetylcysteine (NAC, 40 mM), a known antioxidant that augments the glutathione detoxification pathway and a novel peptide (pNaKtide, 25 mg/kg) which is targeted to interrupting a specific Src-kinase mediated pro-oxidant amplification mechanism. Histological analysis showed significant increase in hepatic inflammation in NAFLD mice exposed to MC-LR which was attenuated on treatment with both NAC and pNaKtide (both p ≤ 0.05). Oxidative stress, as measured by 8-OHDG levels in urine and protein carbonylation in liver sections, was also significantly downregulated upon treatment with both antioxidants after MC-LR exposure. Genetic analysis of key drug transporters including Abcb1a, Phase I enzyme-Cyp3a11 and Phase II metabolic enzymes-Pkm (Pyruvate kinase, muscle), Pklr (Pyruvate kinase, liver, and red blood cell) and Gad1 (Glutamic acid decarboxylase) was significantly altered by MC-LR exposure as compared to the non-exposed control group (all p ≤ 0.05). These changes were significantly attenuated with both pNaKtide and NAC treatment. These results suggest that MC-LR metabolism and detoxification is significantly impaired in the setting of NAFLD, and that these pathways can potentially be reversed with targeted antioxidant treatment.

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.

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.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing.

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

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.

Whole-Genome Sequencing of Bacterial Isolates That Degrade the Cyanobacterial Toxin Microcystin-LR.

We previously demonstrated that 13 bacterial isolates from Lake Erie, when grown in groups of four to five isolates per group, degraded the cyanobacterial toxin microcystin-LR (MC-LR) into nontoxic fragments. Whole-genome sequencing of these bacteria was performed to provide genus and species information and to predict putative MC-LR-degrading genes.

CD40 Receptor Knockout Protects against Microcystin-LR (MC-LR) Prolongation and Exacerbation of Dextran Sulfate Sodium (DSS)-Induced Colitis.

Inflammatory Bowel Disease (IBD) is one of the most common gastrointestinal (GI) disorders around the world, and includes diagnoses such as Crohn's disease and ulcerative colitis. The etiology of IBD is influenced by genetic and environmental factors. One environmental perturbagen that is not well studied within the intestines is microcystin-leucine arginine (MC-LR), which is a toxin produced by cyanobacteria in freshwater environments around the world. We recently reported that MC-LR has limited effects within the intestines of healthy mice, yet interestingly has significant toxicity within the intestines of mice with pre-existing colitis induced by dextran sulfate sodium (DSS). MC-LR was found to prolong DSS-induced weight loss, prolong DSS-induced bloody stools, exacerbate DSS-induced colonic shortening, exacerbate DSS-induced colonic ulceration, and exacerbate DSS-induced inflammatory cytokine upregulation. In addition, we previously reported a significant increase in expression of the pro-inflammatory receptor CD40 in the colons of these mice, along with downstream products of CD40 activation, including plasminogen activator inhibitor-1 (PAI-1) and monocyte chemoattractant protein-1 (MCP-1). In the current study, we demonstrate that knocking out CD40 attenuates the effects of MC-LR in mice with pre-existing colitis by decreasing the severity of weight loss, allowing a full recovery in bloody stools, preventing the exacerbation of colonic shortening, preventing the exacerbation of colonic ulceration, and preventing the upregulation of the pro-inflammatory and pro-fibrotic cytokines IL-1ß, MCP-1, and PAI-1. We also demonstrate the promising efficacy of a CD40 receptor blocking peptide to ameliorate the effects of MC-LR exposure in a proof-of-concept study. Our findings suggest for the first time that MC-LR acts through a CD40-dependent mechanism to exacerbate colitis.

Hyperglycemia induces key genetic and phenotypic changes in human liver epithelial HepG2 cells which parallel the Leprdb/J mouse model of non-alcoholic fatty liver disease (NAFLD).

Non-alcoholic fatty liver disease (NAFLD) is a growing global health concern. With a propensity to progress towards non-alcoholic steatohepatitis (NASH), cirrhosis, and hepatocellular carcinoma, NAFLD is an important link amongst a multitude of comorbidities including obesity, diabetes, and cardiovascular and kidney disease. As several in vivo models of hyperglycemia and NAFLD are employed to investigate the pathophysiology of this disease process, we aimed to characterize an in vitro model of hyperglycemia that was amenable to address molecular mechanisms and therapeutic targets at the cellular level. Utilizing hyperglycemic cell culturing conditions, we induced steatosis within a human hepatocyte cell line (HepG2 cells), as confirmed by electron microscopy. The deposition and accumulation of lipids within hyperglycemic HepG2 cells is significantly greater than in normoglycemic cells, as visualized and quantified by Nile red staining. Alanine aminotransferase (ALT) and alkaline phosphatase (ALP), diagnostic biomarkers for liver damage and disease, were found to be upregulated in hyperglycemic HepG2 cells as compared with normoglycemic cells. Suppression of CEACAM1, GLUT2, and PON1, and elevation of CD36, PCK1, and G6PK were also found to be characteristic in hyperglycemic HepG2 cells compared with normoglycemic cells, suggesting insulin resistance and NAFLD. These in vitro findings mirror the characteristic genetic and phenotypic profile seen in Leprdb/J mice, a well-established in vivo model of NAFLD. In conclusion, we characterize an in vitro model displaying several key genetic and phenotypic characteristics in common with NAFLD that may assist future studies in addressing the molecular mechanisms and therapeutic targets to combat this disease.

Exposure to the Harmful Algal Bloom (HAB) Toxin Microcystin-LR (MC-LR) Prolongs and Increases Severity of Dextran Sulfate Sodium (DSS)-Induced Colitis.

Inflammatory Bowel Disease (IBD) represents a collection of gastrointestinal disorders resulting from genetic and environmental factors. Microcystin-leucine arginine (MC-LR) is a toxin produced by cyanobacteria during algal blooms and demonstrates bioaccumulation in the intestinal tract following ingestion. Little is known about the impact of MC-LR ingestion in individuals with IBD. In this study, we sought to investigate MC-LR's effects in a dextran sulfate sodium (DSS)-induced colitis model. Mice were separated into four groups: (a) water only (control), (b) DSS followed by water (DSS), (c) water followed by MC-LR (MC-LR), and (d) DSS followed by MC-LR (DSS + MC-LR). DSS resulted in weight loss, splenomegaly, and severe colitis marked by transmural acute inflammation, ulceration, shortened colon length, and bloody stools. DSS + MC-LR mice experienced prolonged weight loss and bloody stools, increased ulceration of colonic mucosa, and shorter colon length as compared with DSS mice. DSS + MC-LR also resulted in greater increases in pro-inflammatory transcripts within colonic tissue (TNF-α, IL-1ß, CD40, MCP-1) and the pro-fibrotic marker, PAI-1, as compared to DSS-only ingestion. These findings demonstrate that MC-LR exposure not only prolongs, but also worsens the severity of pre-existing colitis, strengthening evidence of MC-LR as an under-recognized environmental toxin in vulnerable populations, such as those with IBD.

Postural orthostatic tachycardia syndrome is associated with platelet storage pool deficiency.

Mechanisms have been postulated to explain postural orthostatic tachycardia syndrome (POTS), however, the etiology of this often debilitating disorder remains unknown. We conducted a retrospective case-control study of 181 POTS patients who exhibited/reported bleeding symptoms for a specific platelet (PL) dysfunction disorder, delta granule storage pool deficiency (δ-SPD).Patients were included only if results of blood tests for δ-SPD were available. Electron microscopy was utilized to diagnose δ-SPD. An ELISA assay was used to determine serotonin (5HT) concentration in PLs and medical record review was employed to collect patients' clinical symptoms.The most common bleeding symptom was easy bruising (71%) but frequent nose bleeds, heavy menstrual bleeding, and a family history of bleeding were also commonly reported. Of the patients studied, 81% were diagnosed with δ-SPD. Our investigation of 5HT concentration extracted from PLs revealed significantly lower levels of 5HT in POTS patients when compared to that of control subjects. Our data suggest that patients with POTS have significant comorbidities including bleeding symptoms and/or family bleeding histories, and have diminished PL 5HT levels supporting the hypothesis that POTS is a low 5HT level disorder. While we describe a significant relationship with POTS and δ-SPD, this finding does not constitute an etiology for POTS.Our results establish an additional comorbidity frequently seen in POTS that could explain a number of disparate symptoms often affecting the severity of POTS.

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.

Quality control methods for optimal BCR-ABL1 clinical testing in human whole blood samples.

Reliable breakpoint cluster region (BCR)--Abelson (ABL) 1 measurement is essential for optimal management of chronic myelogenous leukemia. There is a need to optimize quality control, sensitivity, and reliability of methods used to measure a major molecular response and/or treatment failure. The effects of room temperature storage time, different primers, and RNA input in the reverse transcription (RT) reaction on BCR-ABL1 and ß-glucuronidase (GUSB) cDNA yield were assessed in whole blood samples mixed with K562 cells. BCR-ABL1 was measured relative to GUSB to control for sample loading, and each gene was measured relative to known numbers of respective internal standard molecules to control for variation in quality and quantity of reagents, thermal cycler conditions, and presence of PCR inhibitors. Clinical sample and reference material measurements with this test were concordant with results reported by other laboratories. BCR-ABL1 per 10(3) GUSB values were significantly reduced (P = 0.004) after 48-hour storage. Gene-specific primers yielded more BCR-ABL1 cDNA than random hexamers at each RNA input. In addition, increasing RNA inhibited the RT reaction with random hexamers but not with gene-specific primers. Consequently, the yield of BCR-ABL1 was higher with gene-specific RT primers at all RNA inputs tested, increasing to as much as 158-fold. We conclude that optimal measurement of BCR-ABL1 per 10(3) GUSB in whole blood is obtained when gene-specific primers are used in RT and samples are analyzed within 24 hours after blood collection.

CEBPG Exhibits Allele-Specific Expression in Human Bronchial Epithelial Cells.

Inter-individual variation in CCAAT/enhancer binding protein gamma (CEBPG) transcript expression in normal human bronchial epithelial cells (NBEC) is associated with predisposition to lung cancer. We hypothesize that this inter-individual variation is in part explained by cis-acting genetic variation in CEBPG. To test this hypothesis we measured transcript expression derived from each parental copy of CEBPG (ie, allele-specific expression; ASE). There was a significant 2.9-fold higher cell cycle-specific variation in ASE of CEBPG rs2772 A compared to C allele (P < 0.001). In 20% of NBEC samples, CEBPG rs2772 A allele was expressed on average 2.10 fold greater than rs2772 C allele. These data support the hypothesis that genetic variation in linkage disequilibrium with rs2772 influences regulation of CEBPG transcript expression through a trans-effect downstream of RNA polymerase II transcription and confirm that cis-acting genetic variation contributes to inter-individual variation in CEBPG transcript expression in NBEC, which is associated with variation in lung cancer risk.

Targeted RNA-sequencing with competitive multiplex-PCR amplicon libraries.

Whole transcriptome RNA-sequencing is a powerful tool, but is costly and yields complex data sets that limit its utility in molecular diagnostic testing. A targeted quantitative RNA-sequencing method that is reproducible and reduces the number of sequencing reads required to measure transcripts over the full range of expression would be better suited to diagnostic testing. Toward this goal, we developed a competitive multiplex PCR-based amplicon sequencing library preparation method that a) targets only the sequences of interest and b) controls for inter-target variation in PCR amplification during library preparation by measuring each transcript native template relative to a known number of synthetic competitive template internal standard copies. To determine the utility of this method, we intentionally selected PCR conditions that would cause transcript amplification products (amplicons) to converge toward equimolar concentrations (normalization) during library preparation. We then tested whether this approach would enable accurate and reproducible quantification of each transcript across multiple library preparations, and at the same time reduce (through normalization) total sequencing reads required for quantification of transcript targets across a large range of expression. We demonstrate excellent reproducibility (R²â€Š= 0.997) with 97% accuracy to detect 2-fold change using External RNA Controls Consortium (ERCC) reference materials; high inter-day, inter-site and inter-library concordance (R²â€Š= 0.97-0.99) using FDA Sequencing Quality Control (SEQC) reference materials; and cross-platform concordance with both TaqMan qPCR (R²â€Šâ€Š= 0.96) and whole transcriptome RNA-sequencing following "traditional" library preparation using Illumina NGS kits (R²â€Š= 0.94). Using this method, sequencing reads required to accurately quantify more than 100 targeted transcripts expressed over a 107-fold range was reduced more than 10,000-fold, from 2.3×109 to 1.4×105 sequencing reads. These studies demonstrate that the competitive multiplex-PCR amplicon library preparation method presented here provides the quality control, reproducibility, and reduced sequencing reads necessary for development and implementation of targeted quantitative RNA-sequencing biomarkers in molecular diagnostic testing.