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.^

The effect of yoga on cardiovascular health as measured by heart rate variability in a sample of breast cancer patients

The study is a three-armed randomized controlled trial comparing values for heart rate variability (HRV), a measure of cardiovascular health, throughout a yoga intervention of breast cancer patients undergoing radiotherapy. Patients attended either a yoga (n=45), stretch, (n=46), or control (n=42) condition 3 times per week for 6 weeks of radiation. Electrocardiograms (ECGs) were conducted on each participant to provide the values necessary for HRV analysis. Analyses focused on examining scores for those participants with HRV baseline values considered to be below the cutoff point for healthy HRV levels, defined by the authors as below the cutpoint of 68 ms. From the entire sample of 133 with available baselines, 26 yogis, 26 stretchers, and 23 controls were determined to be “pathologic” in terms of HRV, and selected for follow-up analysis at 3 weeks and then again at 6 weeks. Though no statistically significant differences were found between either group means at each timepoint or group change score means, the yoga group had consistently higher mean score and mean change scores. These findings are suggestive and indicate the need to refine the use of ECGs and HRV analysis programs to more accurately and comprehensively assess the effects of yoga on cardiovascular health in cancer patients.^

HPV vaccine coverage in Texas counties: An application of multilevel, small area estimation

Health departments, research institutions, policy-makers, and healthcare providers are often interested in knowing the health status of their clients/constituents. Without the resources, financially or administratively, to go out into the community and conduct health assessments directly, these entities frequently rely on data from population-based surveys to supply the information they need. Unfortunately, these surveys are ill-equipped for the job due to sample size and privacy concerns. Small area estimation (SAE) techniques have excellent potential in such circumstances, but have been underutilized in public health due to lack of awareness and confidence in applying its methods. The goal of this research is to make model-based SAE accessible to a broad readership using clear, example-based learning. Specifically, we applied the principles of multilevel, unit-level SAE to describe the geographic distribution of HPV vaccine coverage among females aged 11-26 in Texas.^ Multilevel (3 level: individual, county, public health region) random-intercept logit models of HPV vaccination (receipt of ≥ 1 dose Gardasil® ) were fit to data from the 2008 Behavioral Risk Factor Surveillance System (outcome and level 1 covariates) and a number of secondary sources (group-level covariates). Sampling weights were scaled (level 1) or constructed (levels 2 & 3), and incorporated at every level. Using the regression coefficients (and standard errors) from the final models, I simulated 10,000 datasets for each regression coefficient from the normal distribution and applied them to the logit model to estimate HPV vaccine coverage in each county and respective demographic subgroup. For simplicity, I only provide coverage estimates (and 95% confidence intervals) for counties.^ County-level coverage among females aged 11-17 varied from 6.8-29.0%. For females aged 18-26, coverage varied from 1.9%-23.8%. Aggregated to the state level, these values translate to indirect state estimates of 15.5% and 11.4%, respectively; both of which fall within the confidence intervals for the direct estimates of HPV vaccine coverage in Texas (Females 11-17: 17.7%, 95% CI: 13.6, 21.9; Females 18-26: 12.0%, 95% CI: 6.2, 17.7).^ Small area estimation has great potential for informing policy, program development and evaluation, and the provision of health services. Harnessing the flexibility of multilevel, unit-level SAE to estimate HPV vaccine coverage among females aged 11-26 in Texas counties, I have provided (1) practical guidance on how to conceptualize and conduct modelbased SAE, (2) a robust framework that can be applied to other health outcomes or geographic levels of aggregation, and (3) HPV vaccine coverage data that may inform the development of health education programs, the provision of health services, the planning of additional research studies, and the creation of local health policies.^

Genetic epidemiology of childhood brain tumors

A rare familial cancer syndrome involving childhood brain tumors (CBT), breast cancer, sarcomas and an array of other tumors has been described (Li and Fraumeni 1969, 1975, 1982, 1987). A survey of CBT identified through the Connnecticut Tumor Registry in 1984 revealed a high frequency of CBT, leukemia and other childhood cancer in siblings of CBT patients (Farwell and Flannery, 1984). Other syndromes such as neurofibromatosis and nevoid basal cell carcinoma syndrome have also been associated with CBT; however, no systematic family studies have been conducted to determine the extent to which cancer aggregates in family members of CBT patients. This family study was designed to determine the frequency of cancer aggregation overall or at specific sites, to determine the frequency of known or potentially hereditary syndromes in families of CBT patients, and to determine a genetic model to characterize familial cancer syndromes and to identify specific kindreds to which such a model(s) might apply. This study includes 244 confirmed CBT patients referred to the University of Texas M. D. Anderson Cancer Center between the years 1944 and 1983, diagnosed under the age of 15 years and resident in the U.S. or Canada. Family histories were obtained on the proband's first (parents, siblings and offspring) and second degree (proband's aunts, uncles and grandparents) relatives following sequential sampling scheme rules. To determine if cancer aggregates in families, we compared the cancer experience in the population to that expected in the general population using Connecticut Tumor Registry calendar year, age, race and sex-specific rates. The standardized incidence ratio (SIR) for cancer overall was 0.91 (41 observed (O) and 44.94 expected (E); 95% Confidence Interval (CI) = 0.65-1.24). We observed a significant excess of colon cancer among the proband's first degree relatives (O/E = 5/1.64; 95% CI = 1.01-7.65), in particular those under age 45 year. Segregation analysis showed evidence for multifactorial inheritance in the small percentage (N = 5) of the families. ^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

A Bayesian analysis for spatial processes with application to mapping

In geographical epidemiology, maps of disease rates and disease risk provide a spatial perspective for researching disease etiology. For rare diseases or when the population base is small, the rate and risk estimates may be unstable. Empirical Bayesian (EB) methods have been used to spatially smooth the estimates by permitting an area estimate to "borrow strength" from its neighbors. Such EB methods include the use of a Gamma model, of a James-Stein estimator, and of a conditional autoregressive (CAR) process. A fully Bayesian analysis of the CAR process is proposed. One advantage of this fully Bayesian analysis is that it can be implemented simply by using repeated sampling from the posterior densities. Use of a Markov chain Monte Carlo technique such as Gibbs sampler was not necessary. Direct resampling from the posterior densities provides exact small sample inferences instead of the approximate asymptotic analyses of maximum likelihood methods (Clayton & Kaldor, 1987). Further, the proposed CAR model provides for covariates to be included in the model. A simulation demonstrates the effect of sample size on the fully Bayesian analysis of the CAR process. The methods are applied to lip cancer data from Scotland, and the results are compared. ^

Genetic variants in the mammalian target of rapamycin (mTOR) signaling pathway as predictors of survival and clinical response in women with ovarian cancer

Background. The mTOR pathway is commonly altered in human tumors and promotes cell survival and proliferation. Preliminary evidence suggests this pathway's involvement in chemoresistance to platinum and taxanes, first line therapy for epithelial ovarian cancer. A pathway-based approach was used to identify individual germline single nucleotide polymorphisms (SNPs) and cumulative effects of multiple genetic variants in mTOR pathway genes and their association with clinical outcome in women with ovarian cancer. ^ Methods. The case-series was restricted to 319 non-Hispanic white women with high grade ovarian cancer treated with surgery and platinum-based chemotherapy. 135 SNPs in 20 representative genes in the mTOR pathway were genotyped. Hazard ratios (HRs) for death and Odds ratios (ORs) for failure to respond to primary therapy were estimated for each SNP using the multivariate Cox proportional hazards model and multivariate logistic regression model, respectively, while adjusting for age, stage, histology and treatment sequence. A survival tree analysis of SNPs with a statistically significant association (p<0.05) was performed to identify higher order gene-gene interactions and their association with overall survival. ^ Results. There was no statistically significant difference in survival by tumor histology or treatment regimen. The median survival for the cohort was 48.3 months. Seven SNPs were significantly associated with decreased survival. Compared to those with no unfavorable genotypes, the HR for death increased significantly with the increasing number of unfavorable genotypes and women in the highest risk category had HR of 4.06 (95% CI 2.29–7.21). The survival tree analysis also identified patients with different survival patterns based on their genetic profiles. 13 SNPs on five different genes were found to be significantly associated with a treatment response, defined as no evidence of disease after completion of primary therapy. Rare homozygous genotype of SNP rs6973428 showed a 5.5-fold increased risk compared to the wild type carrying genotypes. In the cumulative effect analysis, the highest risk group (individuals with ≥8 unfavorable genotypes) was significantly less likely to respond to chemotherapy (OR=8.40, 95% CI 3.10–22.75) compared to the low risk group (≤4 unfavorable genotypes). ^ Conclusions. A pathway-based approach can demonstrate cumulative effects of multiple genetic variants on clinical response to chemotherapy and survival. Therapy targeting the mTOR pathway may modify outcome in select patients.^

Stressful life events and their contributions to symptoms of anxiety and depression in adolescents

There are several innovative aspects to this thesis that extend our current knowledge of the relations between stress and psychiatric symptoms in adolescents. First, distal and proximal stressors are differentiated. This approach allows one to specifically examine the role of early childhood stressors as well as stressors experienced more recently as they impact the expression of depression and anxiety during adolescence. Second, a state-of-the-art assessment instrument was used to examine proximal stressors, helping to distinguish several aspects of stress, including objective stress and subjective stress. Third, the parent study from which these data were derived was designed to examine the role of familial risk for depression and related risk factors for the initial development of depression and alcohol use disorders. This allowed for a very thorough collection of demographic characteristics of the study population. Accordingly, this thesis examines the initial prodromal expression of anxiety and depressive symptoms as they are originally expressed prior to the development, if any, of a full-blown psychiatric disorder.^

Identification of micro-RNAs targeting genes of PI3K/AKT pathway in ovarian cancer

Ovarian cancer is the leading cause of cancer-related death for females due to lack of specific early detection method. It is of great interest to find molecular-based biomarkers which are sensitive and specific to ovarian cancer for early diagnosis, prognosis and therapeutics. miRNAs have been proposed to be potential biomarkers that could be used in cancer prevention and therapeutics. The current study analyzed the miRNA and mRNA expression data extracted from the Cancer Genome Atlas (TCGA) database. Using simple linear regression and multiple regression models, we found 71 miRNA-mRNA pairs which were negatively associated between 56 miRNAs and 24 genes of PI3K/AKT pathway. Among these miRNA and mRNA target pairs, 9 of them were in agreement with the predictions from the most commonly used target prediction programs including miRGen, miRDB, miRTarbase and miR2Disease. These shared miRNA-mRNA pairs were considered to be the most potential genes that were involved in ovarian cancer. Furthermore, 4 of the 9 target genes encode cell cycle or apoptosis related proteins including Cyclin D1, p21, FOXO1 and Bcl2, suggesting that their regulator miRNAs including miR-16, miR-96 and miR-21 most likely played important roles in promoting tumor growth through dysregulated cell cycle or apoptosis. miR-96 was also found to directly target IRS-1. In addition, the results showed that miR-17 and miR-9 may be involved in ovarian cancer through targeting JAK1. This study might provide evidence for using miRNA or miRNA profile as biomarker.^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

Pathway analysis of genome-wide association study data in primary biliary cirrhosisa

Genome-wide association studies (GWAS) have successfully identified several genetic loci associated with inherited predisposition to primary biliary cirrhosis (PBC), the most common autoimmune disease of the liver. Pathway-based tests constitute a novel paradigm for GWAS analysis. By evaluating genetic variation across a biological pathway (gene set), these tests have the potential to determine the collective impact of variants with subtle effects that are individually too weak to be detected in traditional single variant GWAS analysis. To identify biological pathways associated with the risk of development of PBC, GWAS of PBC from Italy (449 cases and 940 controls) and Canada (530 cases and 398 controls) were independently analyzed. The linear combination test (LCT), a recently developed pathway-level statistical method was used for this analysis. For additional validation, pathways that were replicated at the P <0.05 level of significance in both GWAS on LCT analysis were also tested for association with PBC in each dataset using two complementary GWAS pathway approaches. The complementary approaches included a modification of the gene set enrichment analysis algorithm (i-GSEA4GWAS) and Fisher's exact test for pathway enrichment ratios. Twenty-five pathways were associated with PBC risk on LCT analysis in the Italian dataset at P<0.05, of which eight had an FDR<0.25. The top pathway in the Italian dataset was the TNF/stress related signaling pathway (p=7.38×10 -4, FDR=0.18). Twenty-six pathways were associated with PBC at the P<0.05 level using the LCT in the Canadian dataset with the regulation and function of ChREBP in liver pathway (p=5.68×10-4, FDR=0.285) emerging as the most significant pathway. Two pathways, phosphatidylinositol signaling system (Italian: p=0.016, FDR=0.436; Canadian: p=0.034, FDR=0.693) and hedgehog signaling (Italian: p=0.044, FDR=0.636; Canadian: p=0.041, FDR=0.693), were replicated at LCT P<0.05 in both datasets. Statistically significant association of both pathways with PBC genetic susceptibility was confirmed in the Italian dataset on i-GSEA4GWAS. Results for the phosphatidylinositol signaling system were also significant in both datasets on applying Fisher's exact test for pathway enrichment ratios. This study identified a combination of known and novel pathway-level associations with PBC risk. If functionally validated, the findings may yield fresh insights into the etiology of this complex autoimmune disease with possible preventive and therapeutic application.^

Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

Survival analysis of circulating tumor cells in triple-negative breast cancer

Background. Breast cancer is the most frequently diagnosed cancer and the leading cause of cancer death among females, accounting for 23% (1.38 million) of the total new cancer cases and 14% (458,400) of the total cancer deaths in 2008. [1] Triple-negative breast cancer (TNBC) is an aggressive phenotype comprising 10–20% of all breast cancers (BCs). [2-4] TNBCs show absence of estrogen, progesterone and HER2/neu receptors on the tumor cells. Because of the absence of these receptors, TNBCs are not candidates for targeted therapies. Circulating tumor cells (CTCs) are observed in blood of breast cancer patients even at early stages (Stage I & II) of the disease. Immunological and molecular analysis can be used to detect the presence of tumor cells in the blood (Circulating tumor cells; CTCs) of many breast cancer patients. These cells may explain relapses in early stage breast cancer patients even after adequate local control. CTC detection may be useful in identifying patients at risk for disease progression, and therapies targeting CTCs may improve outcome in patients harboring them. Methods . In this study we evaluated 80 patients with TNBC who are enrolled in a larger prospective study conducted at M D Anderson Cancer Center in order to determine whether the presence of circulating tumor cells is a significant prognostic factor in relapse free and overall survival . Patients with metastatic disease at the time of presentation were excluded from the study. CTCs were assessed using CellSearch System™ (Veridex, Raritan, NJ). CTCs were defined as nucleated cells lacking the presence of CD45 but expressing cytokeratins 8, 18 or 19. The distribution of patient and tumor characteristics was analyzed using chi square test and Fisher's exact test. Log rank test and Cox regression analysis was applied to establish the association of circulating tumor cells with relapse free and overall survival. Results. The median age of the study participants was 53years. The median duration of follow-up was 40 months. Eighty-eight percent (88%) of patients were newly diagnosed (without a previous history of breast cancer), and (60%) of patients were chemo naïve (had not received chemotherapy at the time of their blood draw for CTC analysis). Tumor characteristics such as stage (P=0.40), tumor size (P=69), sentinel nodal involvement (P=0.87), axillary lymph node involvement (P=0.13), adjuvant therapy (P=0.83), and high histological grade of tumor (P=0.26) did not predict the presence of CTCs. However, CTCs predicted worse relapse free survival (1 or more CTCs log rank P value = 0.04, at 2 or more CTCs P = 0.02 and at 3 or more CTCs P < 0.0001) and overall survival (at 1 or more CTCs log rank P value = 0.08, at 2 or more CTCs P = 0.01 and at 3 or more CTCs P = 0.0001. Conclusions. The number of circulating tumor cells predicted worse relapse free survival and overall survival in TNBC patients.^

Identification of genetic risk factors for cerebellar mutism in pediatric brain tumor patients

Following posterior fossa surgery for resection of childhood medulloblastoma and primitive neuroectodermal tumor (M/PNET), cerebellar mutism (CM) may develop. This is a condition of absent or diminished speech in a conscious patient with no evidence of oral apraxia, which can be accompanied by other symptoms of the posterior fossa syndrome complex, which includes ataxia and hypotonia. Little is known about the etiology. Therefore, we conducted a SNP, gene, and pathway-level analysis to assess the role of host genetic variation on the risk of CM in M/PNET subjects following treatment. Cases (n= 20) and controls (n= 53) were recruited from the Childhood Cancer Epidemiology and Prevention Center, in Houston, TX. DNA samples were genotyped using the Illumina Human 1M Quad SNP chip. Ten pathways were identified from logistic regression used to identify the marginal effect of each SNP on CM risk. The minP test was conducted to identify associations between SNPs categorized to genes and CM risk. Pathways were assessed to determine if there was a significant enrichment of genes in the pathway compared to all other pathways. There were 78 genes that reached the threshold of min P ≤0.05 in 948 genes. The Neurotoxicity pathway was the most significant pathway after adjusting for multiple comparisons (q=0.040 and q=0.005, using Fisher's exact test and a test of proportions, respectively). Most genes within the Neurotoxicity pathway that reached a threshold of minP ≤0.05 were known to have an apoptosis function, possibly inducing neuronal apoptosis in the dentatothalamocortical pathway, and may be important in CM etiology in this population. This is the first study to assess the potential role of genetic risk factors on CM. As an exploratory study, these results should be replicated in a larger sample. ^