Micro RNA expression alterations in patients with glioblastoma

Glioblastoma, also known as glioblastoma multiform or GBM, is the most common and most malignant primary brain tumor. The clinical history of patients with glioblastoma is short, usually less than 3 months in more than 50% of cases after diagnosis. Currently, the methods of glioblastoma treatment are chemotherapy, radiotherapy and surgery. Even with the more effective treatment options, patients with glioblastoma most likely have a median survival time of 10 to 12 months. It is necessary to seek other treatment methods, including gene-targeted treatment. The success of gene-targeted treatment depends critically on the knowledge of genes that may be the cause of, or contribute to disease. To establish a correlate between glioblastoma survival timeline and micro RNA expression alteration, a study of 91 glioblastoma patients was conducted at the University of Texas M. D. Anderson Cancer Center. These 91 glioblastoma patients were newly diagnosed from 2002 to 2007. Statistical analysis was conducted to test the association of miRNA expression alteration between long-term survival and short-term survival glioblastoma. The completion of this proposed study will provide a better understanding of the regulatory role of miRNA in glioblastoma progression.^

Performance evaluation in Bayesian adaptive randomization

Bayesian adaptive randomization (BAR) is an attractive approach to allocate more patients to the putatively superior arm based on the interim data while maintains good statistical properties attributed to randomization. Under this approach, patients are adaptively assigned to a treatment group based on the probability that the treatment is better. The basic randomization scheme can be modified by introducing a tuning parameter, replacing the posterior estimated response probability, setting a boundary to randomization probabilities. Under randomization settings comprised of the above modifications, operating characteristics, including type I error, power, sample size, imbalance of sample size, interim success rate, and overall success rate, were evaluated through simulation. All randomization settings have low and comparable type I errors. Increasing tuning parameter decreases power, but increases imbalance of sample size and interim success rate. Compared with settings using the posterior probability, settings using the estimated response rates have higher power and overall success rate, but less imbalance of sample size and lower interim success rate. Bounded settings have higher power but less imbalance of sample size than unbounded settings. All settings have better performance in the Bayesian design than in the frequentist design. This simulation study provided practical guidance on the choice of how to implement the adaptive design. ^

Gastric cancer survival amongst Hispanics- A Hispanic paradox? A case series study using the SEER 1973--2005 data registry

Stomach cancer is the fourth most common cancer in the world, and ranked 16th in the US in 2008. The age-adjusted rates among Hispanics were 2.8 times that of non-Hispanic Whites in 1998-2002. In spite of that, previous research has found that Hispanics with non-cardia adenocarcinoma of the stomach have a slightly better survival than non-Hispanic Whites. However, such previous research did not include a comparison with African-Americans, and it was limited to data released for the years 1973-2000 in the nine original Surveillance, Epidemiology, and End Results Cancer Registries. This finding was interpreted as related to the Hispanic Paradox, a phenomenon that refers to the fact that Hispanics in the USA tend to paradoxically have substantially better health than other ethnic groups in spite of what their aggregate socio-economic indicators would predict. We extended such research to the SEER 17 Registry, 1973-2005, with varying years of diagnosis per registry, and compared the survival of non-cardia adenocarcinoma of the stomach according to ethnicity (Hispanics, non-Hispanic Whites and African-Americans), while controlling for age, gender, marital status, stage of disease and treatment using Cox regression survival analysis. We found that Hispanic ethnicity by itself did not confer an advantage on survival from non-cardia adenocarcinoma of the stomach, but that being born abroad was independently associated with the apparent 'Hispanic Paradox' previously reported, and that such advantage was seen among foreign born persons across all race/ethnic groups.^

Genome-wide gene-gene interaction analysis for cardiovascular disease

Numerous studies have been carried out to try to better understand the genetic predisposition for cardiovascular disease. Although it is widely believed that multifactorial diseases such as cardiovascular disease is the result from effects of many genes which working alone or interact with other genes, most genetic studies have been focused on identifying of cardiovascular disease susceptibility genes and usually ignore the effects of gene-gene interactions in the analysis. The current study applies a novel linkage disequilibrium based statistic for testing interactions between two linked loci using data from a genome-wide study of cardiovascular disease. A total of 53,394 single nucleotide polymorphisms (SNPs) are tested for pair-wise interactions, and 8,644 interactions are found to be significant with p-values less than 3.5×10-11. Results indicate that known cardiovascular disease susceptibility genes tend not to have many significantly interactions. One SNP in the CACNG1 (calcium channel, voltage-dependent, gamma subunit 1) gene and one SNP in the IL3RA (interleukin 3 receptor, alpha) gene are found to have the most significant pair-wise interactions. Findings from the current study should be replicated in other independent cohort to eliminate potential false positive results.^

An ordinal logistic regression model with misclassification of the outcome variable and categorical covariate

Ordinal outcomes are frequently employed in diagnosis and clinical trials. Clinical trials of Alzheimer's disease (AD) treatments are a case in point using the status of mild, moderate or severe disease as outcome measures. As in many other outcome oriented studies, the disease status may be misclassified. This study estimates the extent of misclassification in an ordinal outcome such as disease status. Also, this study estimates the extent of misclassification of a predictor variable such as genotype status. An ordinal logistic regression model is commonly used to model the relationship between disease status, the effect of treatment, and other predictive factors. A simulation study was done. First, data based on a set of hypothetical parameters and hypothetical rates of misclassification was created. Next, the maximum likelihood method was employed to generate likelihood equations accounting for misclassification. The Nelder-Mead Simplex method was used to solve for the misclassification and model parameters. Finally, this method was applied to an AD dataset to detect the amount of misclassification present. The estimates of the ordinal regression model parameters were close to the hypothetical parameters. β1 was hypothesized at 0.50 and the mean estimate was 0.488, β2 was hypothesized at 0.04 and the mean of the estimates was 0.04. Although the estimates for the rates of misclassification of X1 were not as close as β1 and β2, they validate this method. X 1 0-1 misclassification was hypothesized as 2.98% and the mean of the simulated estimates was 1.54% and, in the best case, the misclassification of k from high to medium was hypothesized at 4.87% and had a sample mean of 3.62%. In the AD dataset, the estimate for the odds ratio of X 1 of having both copies of the APOE 4 allele changed from an estimate of 1.377 to an estimate 1.418, demonstrating that the estimates of the odds ratio changed when the analysis includes adjustment for misclassification. ^

Survival of gastrointestinal and stromal tumor patients at UT M.D. Anderson Cancer Center

Objective. Gastrointestinal Stromal Tumors (GISTs) are rare mesenchymal tumors of the gastrointestinal (GI) tract with spindled cell, epithelioid, or occasionally pleomorphic morphology. The primary objective of this paper is to describe the demographic and clinical characteristics and survival among GIST patients registered at the University of Texas M.D. Anderson Cancer Center (MDACC). ^ Methods. This cohort study includes 783 consecutive patients diagnosed with GIST from 1995 to 2007. Demographic, clinical and survival information were obtained from the MDACC cancer registry. ^ Statistical Analysis. Kaplan-Meier survival curves, univariate and multivariate Cox proportional hazards analysis were conducted to estimate survival and identify prognostic clinical factors associated with survival. Results. The age at diagnosis of MDACC GIST cases ranged from 17 to 91 with a mean of 57 years and a male-to-female ratio of 1.3:1. The racial distribution was whites 77%, African-Americans 9.5%, Hispanics 9.3% and other races 4.2%. Fifty per cent of the GISTs arose from stomach, 35% small intestine, 7% retroperitoneal space, 6% colorectal and 2% were omentum and mesentery. About half of the tumors were less than 10 cm in size. Fifty eight per cent of the tumors were localized whereas 36% were metastatic. MDACC GIST patients were generally comparable to SEER patients, but, on the average, were 7 years younger than SEER patients and were predominantly whites. ^ Stratification of 783 GIST cases by year of diagnosis based on the introduction of imatinib treatment in 2000 revealed that 60% of the GIST cases were first diagnosed between 2000 and 2007 whereas, 40% were first diagnosed between 1995 and 1999. There was a significant difference between the two cohorts in the distribution of race, GIST symptom, tumor size, tumor site, and stage of the tumor at diagnosis. The 1- and 5-year survival was 93% and 59% in the 1995–2007 cohort. Multivariate Cox regression analysis identified age at diagnosis (p<0.001), female sex (p=0.047), tumor size (p=0.07), multiple cancers (p=0.002), and GIST diagnosed between 2000 and 2007 (p<0.001) were significantly associated with survival. Approximately, 58% of the cases were treated with imatinib whereas 42% did not receive imatinib in 2000–2005 cohort. There was a significant difference in survival between imatinib and non-imatinib groups and in the distribution of tumor size categories, stage of the tumor at diagnosis and cancers before the diagnosis of GIST. The 1- and 5-year survival for imatinib patients was 99% and 73% and was 91% and 63% for non-imatinib patients. Multivariate Cox regression analysis of the 2000–2007 cohort identified, age at diagnosis and tumor stage as possible prognostic factors associated with survival.^

A multivariate frailty model for disease recurrences and survival

A multivariate frailty hazard model is developed for joint-modeling of three correlated time-to-event outcomes: (1) local recurrence, (2) distant recurrence, and (3) overall survival. The term frailty is introduced to model population heterogeneity. The dependence is modeled by conditioning on a shared frailty that is included in the three hazard functions. Independent variables can be included in the model as covariates. The Markov chain Monte Carlo methods are used to estimate the posterior distributions of model parameters. The algorithm used in present application is the hybrid Metropolis-Hastings algorithm, which simultaneously updates all parameters with evaluations of gradient of log posterior density. The performance of this approach is examined based on simulation studies using Exponential and Weibull distributions. We apply the proposed methods to a study of patients with soft tissue sarcoma, which motivated this research. Our results indicate that patients with chemotherapy had better overall survival with hazard ratio of 0.242 (95% CI: 0.094 - 0.564) and lower risk of distant recurrence with hazard ratio of 0.636 (95% CI: 0.487 - 0.860), but not significantly better in local recurrence with hazard ratio of 0.799 (95% CI: 0.575 - 1.054). The advantages and limitations of the proposed models, and future research directions are discussed. ^

Ethnic differences in tumor characteristics, treatment and survival in glioma patients

Purpose. A descriptive analysis of glioma patients by race was carried out in order to better elucidate potential differences between races in demographics, treatment, characteristics, prognosis and survival. ^ Patients and Methods. Among 1,967 patients ≥ 18 years diagnosed with glioma seen between July 2000 and September 2006 at The University of Texas M.D. Anderson Cancer Center (UTMDACC). Data were collated from the UTMDACC Patient History Database (PHDB) and the UTMDACC Tumor Registry Database (TRDB). Chi-square analysis, uni- /multivariate Cox proportional hazards modeling and survival analysis were used to analyze differences by race. ^ Results. Demographic, treatment and histologic differences exist between races. Though risk differences were seen between races, race was not found to be a significant predictor in multivariate regression analysis after accounting for age, surgery, chemotherapy, radiation, tumor type as stratified by WHO tumor grade. Age was the most consistent predictor in risk for death. Overall survival by race was significantly different (p=0.0049) only in low-grade gliomas after adjustment for age although survival differences were very slight. ^ Conclusion. Among this cohort of glioma patients, age was the strongest predictor for survival. It is likely that survival is more influenced by age, time to treatment, tumor grade and surgical expertise rather than racial differences. However, age at diagnosis, gender ratios, histology and history of cancer differed significantly between race and genetic differences to this effect cannot be excluded. ^

Involvement of lipocalin 2 in leukemia and breast cancer

Human lipocalin 2 is described as the neutrophil gelatinase-associated lipocalin (NGAL). The lipocalin 2 gene encodes a small, secreted glycoprotein that possesses a variety of functions, of which the best characterized function is organic iron binding activity. Elevated NGAL expression has been observed in many human cancers including breast, colorectal, pancreatic and ovarian cancers. I focused on the characterization of NGAL function in chronic myelogenous leukemia (CML) and breast cancer. Using the leukemic xenograft mouse model, we demonstrated that over-expression of NGAL in K562 cells, a leukemic cell line, led to a higher apoptotic rate and an atrophy phenotype in the spleen of inoculated mice compared to K562 cells alone. These results indicate that NGAL plays a primary role in suppressing hematopoiesis by inducing apoptosis within normal hematopoietic cells. In the breast cancer project, we analyzed two microarray data sets of breast cancer cell lines ( n = 54) and primary breast cancer samples (n = 318), and demonstrated that high NGAL expression is significantly correlated with several tumor characteristics, including negative estrogen receptor (ER) status, positive HER2 status, high tumor grade, and lymph node metastasis. Ectopic NGAL expression in non-aggressive (ZR75.1 and MCF7) cells led to aggressive tumor phenotypes in vitro and in vivo. Conversely, knockdown of NGAL expression in various breast cancer cell lines by shRNA lentiviral infection significantly decreased migration, invasion, and metastasis activities of tumor cells both in vitro and in vivo . It has been previously reported that transgenic mice with a mutation in the region of trans-membrane domain (V664E) of HER2 develop mammary tumors that progress to lung metastasis. However, we observed that genetic deletion of the 24p3 gene, a mouse homolog of NGAL, in HER2 transgenic mice by breeding with 24p3-null mice resulted in a significant delay of mammary tumor formation and reduction of lung metastasis. Strikingly, we also found that treatment with affinity purified 24p3 antibodies in the 4T1 breast cancer mice strongly reduced lung metastasis. Our studies provide evidence that NGAL plays a critical role in breast cancer development and progression, and thus NGAL has potential as a new therapeutic target in breast cancer.^

Systematic review of genetic risk score in coronary heart disease and other diseases

In order to better take advantage of the abundant results from large-scale genomic association studies, investigators are turning to a genetic risk score (GRS) method in order to combine the information from common modest-effect risk alleles into an efficient risk assessment statistic. The statistical properties of these GRSs are poorly understood. As a first step toward a better understanding of GRSs, a systematic analysis of recent investigations using a GRS was undertaken. GRS studies were searched in the areas of coronary heart disease (CHD), cancer, and other common diseases using bibliographic databases and by hand-searching reference lists and journals. Twenty-one independent case-control studies, cohort studies, and simulation studies (12 in CHD, 9 in other diseases) were identified. The underlying statistical assumptions of the GRS using the experience of the Framingham risk score were investigated. Improvements in the construction of a GRS guided by the concept of composite indicators are discussed. The GRS will be a promising risk assessment tool to improve prediction and diagnosis of common diseases.^

Variance of familial dietary intake

Although many family-based genetic studies have collected dietary data, very few have used the dietary information in published findings. No single solution has been presented or discussed in the literature to deal with the problem of using factor analyses for the analyses of dietary data from several related individuals from a given household. The standard statistical approach of factor analysis cannot be applied to the VIVA LA FAMILIA Study diet data to ascertain dietary patterns since this population consists of three children from each family, thus the dietary patterns of the related children may be correlated and non-independent. Addressing this problem in this project will enable us to describe the dietary patterns in Hispanic families and to explore the relationships between dietary patterns and childhood obesity. ^ In the VIVA LA FAMILIA Study, an overweight child was first identified and then his/her siblings and parents were brought in for data collection which included 24 hour recalls and food frequency questionnaire (FFQ). Dietary intake data were collected using FFQ and 24 hour recalls on 1030 Hispanic children from 319 families. ^ The design of the VIVA LA FAMILIA Study has important and unique statistical considerations since its participants are related to each other, the majority form distinct nuclear families. Thus, the standard approach of factor analysis cannot be applied to these diet data to ascertain dietary patterns. In this project we propose to investigate whether the determinants of the correlation matrix of each family unit will allow us to adjust the original correlation matrix of the dietary intake data prior to ascertaining dietary intake patterns. If these methods are appropriate, then in the future the dietary patterns among related individuals could be assessed by standard orthogonal principal component factor analysis.^

The combined effect of pravastatin and aspirin on clinical cardiovascular events

The 3-hydroxy-3methylglutaryl coenzyme A (HMG-CoA) reductase inhibitors, or statins, can achieve significant reductions in plasma low-density lipoprotein (LDL)-cholesterol levels. Experimental and clinical evidence now shows that some statins interfere with formation of atherosclerotic lesions independent of their hypolipidemic properties. Vulnerable plaque rupture can result in thrombus formation and artery occlusion; this plaque deterioration is responsible for most acute coronary syndromes, including myocardial infarction (MI), unstable angina, and coronary death, as well as coronary heart diseaseequivalent non-hemorrhagic stroke. Inhibition of HMG-CoA reductase has potential pleiotropic effects other than lipid-lowering, as statins block mevalonic acid production, a precursor to cholesterol and numerous other metabolites. Statins' beneficial effects on clinical events may also thus involve nonlipid-related mechanisms that modify endothelial function, inflammatory responses, plaque stability, and thrombus formation. Aspirin, routinely prescribed to post-MI patients as adjunct therapy, may potentiate statins beneficial effects, as aspirin does not compete metabolically with statins but acts similarly on atherosclerotic lesions. Common functions of both medications include inhibition of platelet activity and aggregation, reduction in atherosclerotic plaque macrophage cell count, and prevention of atherosclerotic vessel endothelial dysfunction. The Cholesterol and Recurrent Events (CARE) trial provides an ideal population in which to examine the combined effects of pravastatin and aspirin. Lipid levels, intermediate outcomes, are examined by pravastatin and aspirin status, and differences between the two pravastatin groups are found. A modified Cox proportional-hazards model with aspirin as a time-dependent covariate was used to determine the effect of aspirin and pravastatin on the clinical cardiovascular composite endpoint of coronary heart disease death, recurrent MI or stroke. Among those assigned to pravastatin, use of aspirin reduced the composite primary endpoint by 35%; this result was similar by gender, race, and diabetic status. Older patients demonstrated a nonsignificant 21% reduction in the primary outcome, whereas the younger had a significant reduction of 43% in the composite primary outcome. Secondary outcomes examined include coronary artery bypass graft (38% reduction), nonsurgical bypass, peripheral vascular disease, and unstable angina. Pravastatin and aspirin in a post-MI population was found to be a beneficial combination that seems to work through lipid and nonlipid, anti-inflammatory mechanisms. ^

Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

Assessing the improved discriminatory power of a new biomarker in prognostic models

Although the area under the receiver operating characteristic (AUC) is the most popular measure of the performance of prediction models, it has limitations, especially when it is used to evaluate the added discrimination of a new biomarker in the model. Pencina et al. (2008) proposed two indices, the net reclassification improvement (NRI) and integrated discrimination improvement (IDI), to supplement the improvement in the AUC (IAUC). Their NRI and IDI are based on binary outcomes in case-control settings, which do not involve time-to-event outcome. However, many disease outcomes are time-dependent and the onset time can be censored. Measuring discrimination potential of a prognostic marker without considering time to event can lead to biased estimates. In this dissertation, we have extended the NRI and IDI to survival analysis settings and derived the corresponding sample estimators and asymptotic tests. Simulation studies were conducted to compare the performance of the time-dependent NRI and IDI with Pencina’s NRI and IDI. For illustration, we have applied the proposed method to a breast cancer study.^ Key words: Prognostic model, Discrimination, Time-dependent NRI and IDI ^

Performance recovery from exercise: A review from a statistical point of view

Background. Research into methods for recovery from fatigue due to exercise is a popular topic among sport medicine, kinesiology and physical therapy. However, both the quantity and quality of studies and a clear solution of recovery are lacking. An analysis of the statistical methods in the existing literature of performance recovery can enhance the quality of research and provide some guidance for future studies. Methods: A literature review was performed using SCOPUS, SPORTDiscus, MEDLINE, CINAHL, Cochrane Library and Science Citation Index Expanded databases to extract the studies related to performance recovery from exercise of human beings. Original studies and their statistical analysis for recovery methods including Active Recovery, Cryotherapy/Contrast Therapy, Massage Therapy, Diet/Ergogenics, and Rehydration were examined. Results: The review produces a Research Design and Statistical Method Analysis Summary. Conclusion: Research design and statistical methods can be improved by using the guideline from the Research Design and Statistical Method Analysis Summary. This summary table lists the potential issues and suggested solutions, such as, sample size calculation, sports specific and research design issues consideration, population and measure markers selection, statistical methods for different analytical requirements, equality of variance and normality of data, post hoc analyses and effect size calculation.^