New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

The Role of Acidic Fibroblast Growth Factor and Fibroblast Growth Factor Receptor 4 in High Grade Serous Carcinoma

Background: High grade serous carcinoma whether ovarian, tubal or primary peritoneal, continues to be the most lethal gynecologic malignancy in the USA. Although combination chemotherapy and aggressive surgical resection has improved survival in the past decade the majority of patients still succumb to chemo-resistant disease recurrence. It has recently been reported that amplification of 5q31-5q35.3 is associated with poor prognosis in patients with high grade serous ovarian carcinoma. Although the amplicon contains over 50 genes, it is notable for the presence of several members of the fibroblast growth factor signaling axis. In particular acidic fibroblast growth factor (FGF1) has been demonstrated to be one of the driving genes in mediating the observed prognostic effect of the amplicon in ovarian cancer patients. This study seeks to further validate the prognostic value of fibroblast growth receptor 4 (FGFR4), another candidate gene of the FGF/FGFR axis located in the same amplicon. The emphasis will be delineating the role the FGF1/FGFR4 signaling axis plays in high grade serous ovarian carcinoma; and test the feasibility of targeting the FGF1/FGFR4 axis therapeutically. Materials and Methods: Spearman and Pearson correlation studies on data generated from array CGH and transcriptome profiling analyses on 51 microdissected tumor samples were used to identify genes located on chromosome 5q31-35.3 that showed significant correlation between DNA and mRNA copy numbers. Significant correlation between FGF1 and FGFR4 DNA copy numbers was further validated by qPCR analysis on DNA isolated from 51 microdissected tumor samples. Immunolocalization and quantification of FGFR4 expression were performed on paraffin embedded tissue samples from 183 cases of high-grade serous ovarian carcinoma. The expression was then correlated with clinical data to assess impact on survival. The expression of FGF1 and FGFR4 in vitro was quantified by real-time PCR and western blotting in six high-grade serous ovarian carcinoma cell lines and compared to those in human ovarian surface epithelial cells to identify overexpression. The effect of FGF1 on these cell lines after serum starvation was quantified for in vitro cellular proliferation, migration/invasion, chemoresistance and survival utilizing a combination of commercially available colorimetric, fluorometric and electrical impedance assays. FGFR4 expression was then transiently silenced via siRNA transfection and the effects on response to FGF1, cellular proliferation, and migration were quantified. To identify relevant cellular pathways involved, responsive cell lines were transduced with different transcription response elements using the Cignal-Lenti reporter system and treated with FGF1 with and without transient FGFR4 knock down. This was followed by western blot confirmation for the relevant phosphoproteins. Anti-FGF1 antibodies and FGFR trap proteins were used to attempt inhibition of FGF mediated phenotypic changes and relevant signaling in vitro. Orthotopic intraperitoneal tumors were established in nude mice using serous cell lines that have been previously transfected with luciferase expressing constructs. The mice were then treated with FGFR trap protein. Tumor progression was then followed via bioluminescent imaging. The FGFR4 gene from 52 clinical samples was sequenced to screen for mutations. Results: FGFR4 DNA and mRNA copy numbers were significantly correlated and FGFR4 DNA copy number was significantly correlated with that of FGF1. Survival of patients with high FGFR4 expressing tumors was significantly shorter that those with low expression(median survival 28 vs 55 month p< 0.001) In a multivariate cox regression model FGFR expression significantly increased risk of death (HR 2.1, p<0.001). FGFR4 expression was significantly higher in all cell lines tested compared to HOSE, OVCA432 cell line in particular had very high expression suggesting amplification. FGF1 was also particularly overexpressed in OVCA432. FGF1 significantly increased cell survival after serum deprivation in all cell lines. Transient knock down of FGFR4 caused significant reduction in cell migration and proliferation in vitro and significantly decreased the proliferative effects of FGF1 in vitro. FGFR1, FGFR4 traps and anti-FGF1 antibodies did not show activity in vitro. OVCA432 transfected with the cignal lenti reporter system revealed significant activation of MAPK, NFkB and WNT pathways, western blotting confirmed the results. Reverse phase protein array (RPPA) analysis also showed activation of MAPK, AKT, WNT pathways and down regulation of E Cadherin. FGFR trap protein significantly reduced tumor growth in vivo in an orthotopic mouse model. Conclusions: Overexpression and amplification of several members of the FGF signaling axis present on the amplicon 5q31-35.3 is a negative prognostic indicator in high grade serous ovarian carcinoma and may drive poor survival associated with that amplicon. Activation of The FGF signaling pathway leads to downstream activation of MAPK, AKT, WNT and NFkB pathways leading to a more aggressive cancer phenotype with increased tumor growth, evasion of apoptosis and increased migration and invasion. Inhibition of FGF pathway in vivo via FGFR trap protein leads to significantly decreased tumor growth in an orthotopic mouse model.

Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

Response to chemotherapy, recurrence and survival in advanced-stage ovarian, fallopian tube and primary peritoneal cancer patients with non-Ashkenazi Jewish BRCA mutations, compared to those without

Objectives. Previous studies have shown a survival advantage in ovarian cancer patients with Ashkenazi-Jewish (AJ) BRCA founder mutations, compared to sporadic ovarian cancer patients. The purpose of this study was to determine if this association exists in ovarian cancer patients with non-Ashkenazi Jewish BRCA mutations. In addition, we sought to account for possible "survival bias" by minimizing any lead time that may exist between diagnosis and genetic testing. ^ Methods. Patients with stage III/IV ovarian, fallopian tube, or primary peritoneal cancer and a non-Ashkenazi Jewish BRCA1 or 2 mutation, seen for genetic testing January 1996-July 2007, were identified from genetics and institutional databases. Medical records were reviewed for clinical factors, including response to initial chemotherapy. Patients with sporadic (non-hereditary) ovarian, fallopian tube, or primary peritoneal cancer, without family history of breast or ovarian cancer, were compared to similar cases, matched by age, stage, year of diagnosis, and vital status at time interval to BRCA testing. When possible, 2 sporadic patients were matched to each BRCA patient. An additional group of unmatched, sporadic ovarian, fallopian tube and primary peritoneal cancer patients was included for a separate analysis. Progression-free (PFS) & overall survival (OS) were calculated by the Kaplan-Meier method. Multivariate Cox proportional hazards models were calculated for variables of interest. Matched pairs were treated as clusters. Stratified log rank test was used to calculate survival data for matched pairs using paired event times. Fisher's exact test, chi-square, and univariate logistic regression were also used for analysis. ^ Results. Forty five advanced-stage ovarian, fallopian tube and primary peritoneal cancer patients with non-Ashkenazi Jewish (non-AJ) BRCA mutations, 86 sporadic-matched and 414 sporadic-unmatched patients were analyzed. Compared to the sporadic-matched and sporadic-unmatched ovarian cancer patients, non-AJ BRCA mutation carriers had longer PFS (17.9 & 13.8 mos. vs. 32.0 mos., HR 1.76 [95% CI 1.13–2.75] & 2.61 [95% CI 1.70–4.00]). In relation to the sporadic- unmatched patients, non-AJ BRCA patients had greater odds of complete response to initial chemotherapy (OR 2.25 [95% CI 1.17–5.41]) and improved OS (37.6 mos. vs. 101.4 mos., HR 2.64 [95% CI 1.49–4.67]). ^ Conclusions. This study demonstrates a significant survival advantage in advanced-stage ovarian cancer patients with non-AJ BRCA mutations, confirming the previous studies in the Jewish population. Our efforts to account for "survival bias," by matching, will continue with collaborative studies. ^

Integrating sequence information in microarray data analysis by free energy modeling

Microarray technology is a high-throughput method for genotyping and gene expression profiling. Limited sensitivity and specificity are one of the essential problems for this technology. Most of existing methods of microarray data analysis have an apparent limitation for they merely deal with the numerical part of microarray data and have made little use of gene sequence information. Because it's the gene sequences that precisely define the physical objects being measured by a microarray, it is natural to make the gene sequences an essential part of the data analysis. This dissertation focused on the development of free energy models to integrate sequence information in microarray data analysis. The models were used to characterize the mechanism of hybridization on microarrays and enhance sensitivity and specificity of microarray measurements. ^ Cross-hybridization is a major obstacle factor for the sensitivity and specificity of microarray measurements. In this dissertation, we evaluated the scope of cross-hybridization problem on short-oligo microarrays. The results showed that cross hybridization on arrays is mostly caused by oligo fragments with a run of 10 to 16 nucleotides complementary to the probes. Furthermore, a free-energy based model was proposed to quantify the amount of cross-hybridization signal on each probe. This model treats cross-hybridization as an integral effect of the interactions between a probe and various off-target oligo fragments. Using public spike-in datasets, the model showed high accuracy in predicting the cross-hybridization signals on those probes whose intended targets are absent in the sample. ^ Several prospective models were proposed to improve Positional Dependent Nearest-Neighbor (PDNN) model for better quantification of gene expression and cross-hybridization. ^ The problem addressed in this dissertation is fundamental to the microarray technology. We expect that this study will help us to understand the detailed mechanism that determines sensitivity and specificity on the microarrays. Consequently, this research will have a wide impact on how microarrays are designed and how the data are interpreted. ^

Improving the health and well-being of women at risk for neural tube defect recurrence

There is growing interest in providing women with internatal care, a package of healthcare and ancillary services that can improve their health during the period after the termination of one pregnancy but before the conception of the next pregnancy. Women who have had a pregnancy affected by a neural tube defect can especially benefit from internatal care because they are at increased risk for recurrence and improvements to their health during the inter-pregnancy period can prevent future negative birth outcomes. The dissertation provides three papers that inform the content of internatal care for women at risk for recurrence by examining descriptive epidemiology to develop an accurate risk profile of the population, assessing whether women at risk for recurrence would benefit from a psychosocial intervention, and determining how to improve health promotion efforts targeting folic acid use.^ Paper one identifies information relevant for developing risk profiles and conducting risk assessments. A number of investigations have found that the risk for neural tube defects differs between non-Hispanic Whites and Hispanics. To understand the risk difference, the descriptive epidemiology of spina bifida and anencephaly was examined for Hispanics and non-Hispanic Whites based on data from the Texas Birth Defects Registry for the years 1999 through 2004. Crude and adjusted birth prevalence ratios and corresponding 95% confidence intervals were calculated between descriptive epidemiologic characteristics and anencephaly and spina bifida for non-Hispanic Whites and for Hispanics. In both race/ethnic groups, anencephaly expressed an inverse relationship with maternal age and a positive linear relationship with parity. Both relationships were stronger in non-Hispanic Whites. Female infants had a higher risk for anencephaly in non-Hispanic Whites. Lower maternal education was associated with increased risk for spina bifida in Hispanics.^ Paper two assesses the need for a psychosocial intervention. For mothers who have children with spina bifida, the transition to motherhood can be stressful. This qualitative study explored the process of becoming a mother to a child with spina bifida focusing particularly on stress and coping in the immediate postnatal environment. Semi-structured interviews were conducted with six mothers who have children with spina bifida. Mothers were asked about their initial emotional and problem-based coping efforts, the quality and kind of support provided by health providers, and the characteristics of their meaning-based coping efforts; questions matched Transactional Model of Stress and Coping (TMSC) constructs. Analysis of the responses revealed a number of modifiable stress and coping transactions, the most salient being: health providers are in a position to address beliefs about self-causality and prevent mothers from experiencing the repercussions that stem from maintaining these beliefs. ^ Paper three identifies considerations when creating health promotion materials targeting folic acid use. A brochure was designed using concepts from the Precaution Adoption Process Model (PAPM). Three focus groups comprising 26 mothers of children with spina bifida evaluated the brochure. One focus group was conducted in Spanish-only, the other two focus groups were conducted in English and Spanish combined. Qualitative analysis of coded transcripts revealed that a brochure is a helpful adjunct. Questions about folic acid support the inclusion of an insert with basic information. There may be a need to develop different educational material for Hispanics so the importance of folic acid is provided in a situational context. Some participants blamed themselves for their pregnancy outcome which may affect their receptivity to messages in the brochure. The women's desire for photographs that affect their perception of threat and their identification with the second role model indicate they belong to PAPM Stage 2 and 3. Participants preferred colorful envelopes, high quality paper, intimidating photographs, simple words, conversational style sentences, and positive messages.^ These papers develop the content of risk assessment, psychosocial intervention, and health promotion components of internatal care as they apply to women at risk for recurrence. The findings provided evidence for considering parity and maternal age when assessing nutritional risk. The two dissimilarities between the two race/ethnic groups, infant sex and maternal education lent support to creating separate risk profiles. Interviews with mothers of children with spina bifida revealed the existence of unmet needs-suggesting that a psychosocial intervention provided as part of internatal care can strengthen and support women's well-being. Segmenting the audience according to race/ethnicity and PAPM stage can improve the relevance of print materials promoting folic acid use.^

Enlarged tracheoesophageal puncture: A systematic review and retrospective analysis

Background. An enlarged tracheoesophageal puncture (TEP) results in aspiration around the voice prosthesis (VP) and may lead to pneumonia. The aims of this research were: (1) to conduct a systematic review and meta-analysis on enlarged TEP; (2) to analyze preoperative, perioperative, and postoperative risk factors for enlarged TEP; and (3) to evaluate control of leakage around the VP using conservative treatments and adverse events in patients with enlarged TEP.^ Methods. A systematic review was conducted (1978-2008). A summary risk estimate was calculated using a random-effects meta-analysis model. A retrospective cohort study was completed. Patients who underwent total laryngectomy and TEP at The University of Texas M. D. Anderson Cancer Center (MDACC) were included. Multiple logistic regression methods were used to assess risk factors for enlargement. Descriptive and bivariate statistics were calculated to evaluate outcomes and adverse events. Results: Twenty-seven manuscripts were included in the systematic review. The summary risk estimate of enlarged TEP/leakage around the VP was 7.2% (95% CI: 4.8%-9.6%). Temporary VP removal and TEP-site injections were the most commonly reported treatments. Neither prosthetic diameter (p=0.076) nor timing of TEP (p=0.297) significantly increased risk of enlargement per stratified analyses of published outcomes. The cumulative incidence of enlarged TEP was 18.6% (36/194, 95% CI: 13.0%-24.1%) in the MDACC cohort. Enlarged TEP occurred exclusively in irradiated patients. Adjusting for length of follow-up and timing of TEP, advanced nodal disease (ORadjusted: 4.3, 95% CI: 1.0-19.1), stricture (ORadjusted : 3.2, 95% CI: 1.2-8.6), and locoregional recurrence/distant metastasis after laryngectomy (ORadjusted: 6.2, 95% CI: 2.3-16.4) increased risk of enlarged TEP. At last follow-up, conservative methods controlled leakage around the VP in 81% (29/36) of patients. Unresolved leakage was associated with recurrent cancer (p=0.081) and TEP-site irregularity (p=0.003). Relative to those without enlargement, enlarged TEP patients had significantly higher risk of pneumonia (RR: 3.4, 95% CI: 1.9-6.2).^ Conclusions. These data establish that enlarged TEP poses serious health risks, and provide insight into medical and oncologic factors that may contribute to development of this complication. In addition, this research supports the use of conservative treatments to address leakage after enlarged TEP in lieu of complete TEP closure.^

Venous thromboembolism in patients with acute leukemia: Prevalence, incidence, recurrence of thrombosis, risk factors, and impact on survival

Venous thromboembolism (VTE), including deep vein thrombosis (DVT) and pulmonary embolism (PE), is the third most preventable cardiovascular disease and a growing public health problem in the United States. The incidence of VTE remains high with an annual estimate of more than 600,000 symptomatic events. DVT affects an estimated 2 million American each year with a death toll of 300,000 persons per year from DVT-related PE. Leukemia patients are at high risk for both hemorrhage and thrombosis; however, little is known about thrombosis among acute leukemia patients. The ultimate goal of this dissertation was to obtain deep understanding of thrombotic issue among acute leukemia patients. The dissertation was presented in a format of three papers. First paper mainly looked at distribution and risk factors associated with development of VTE among patients with acute leukemia prior to leukemia treatment. Second paper looked at incidence, risk factors, and impact of VTE on survival of patients with acute lymphoblastic leukemia during treatment. Third paper looked at recurrence and risk factors for VTE recurrence among acute leukemia patients with an initial episode of VTE. Descriptive statistics, Chi-squared or Fisher's exact test, median test, Mann-Whitney test, logistic regression analysis, Nonparametric Estimation Kaplan-Meier with a log-rank test or Cox model were used when appropriate. Results from analyses indicated that acute leukemia patients had a high prevalence, incidence, and recurrent rate of VTE. Prior history of VTE, obesity, older age, low platelet account, presence of Philadelphia positive ALL, use of oral contraceptives or hormone replacement therapy, presence of malignancies, and co-morbidities may place leukemia patients at an increased risk for VTE development or recurrence. Interestingly, development of VTE was not associated with a higher risk of death among hospitalized acute leukemia patients.^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Renal Cryoablation: Investigation of Periprocedural Visualization Tools and Treatment Response Quantification

Cryoablation for small renal tumors has demonstrated sufficient clinical efficacy over the past decade as a non-surgical nephron-sparing approach for treating renal masses for patients who are not surgical candidates. Minimally invasive percutaneous cryoablations have been performed with image guidance from CT, ultrasound, and MRI. During the MRI-guided cryoablation procedure, the interventional radiologist visually compares the iceball size on monitoring images with respect to the original tumor on separate planning images. The comparisons made during the monitoring step are time consuming, inefficient and sometimes lack the precision needed for decision making, requiring the radiologist to make further changes later in the procedure. This study sought to mitigate uncertainty in these visual comparisons by quantifying tissue response to cryoablation and providing visualization of the response during the procedure. Based on retrospective analysis of MR-guided cryoablation patient data, registration and segmentation algorithms were investigated and implemented for periprocedural visualization to deliver iceball position/size with respect to planning images registered within 3.3mm with at least 70% overlap and a quantitative logit model was developed to relate perfusion deficit in renal parenchyma visualized in verification images as a result of iceball size visualized in monitoring images. Through retrospective study of 20 patient cases, the relationship between likelihood of perfusion loss in renal parenchyma and distance within iceball was quantified and iteratively fit to a logit curve. Using the parameters from the logit fit, the margin for 95% perfusion loss likelihood was found to be 4.28 mm within the iceball. The observed margin corresponds well with the clinically accepted margin of 3-5mm within the iceball. In order to display the iceball position and perfusion loss likelihood to the radiologist, algorithms were implemented to create a fast segmentation and registration module which executed in under 2 minutes, within the clinically-relevant 3 minute monitoring period. Using 16 patient cases, the average Hausdorff distance was reduced from 10.1mm to 3.21 mm with average DSC increased from 46.6% to 82.6% before and after registration.