Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Functional analysis of histones H3 and H4 in Tup1 repression and in GCN5 activation

Tup1 forms a complex with Ssn6 in yeast. Ssn6-Tup1 complex is recruited via direct interactions with specific DNA binding proteins to a specific promoter region and mediates repression of several sets of genes including a-cell specific genes (asg) in $\alpha$ cells. It has been shown that repression of asgs also requires histone H4 and that Tup1 can directly interact with H3 and H4 in vitro. To address whether histone H3 is required for the repression of asgs, I have examined the effect of H3 and H4 mutations on the expression of a $\alpha$2-controlled LacZ reporter. Assay of $\beta$-glactosidase shows that mutations in either H3 or H4 cause a weak derepression of the reporter gene. Some double mutations result in a stronger derepression, while others do not. The H3 N-terminal deletion also leads to a slightly decreased expression of the reporter gene in $\alpha$ cells. Our data suggest that the N-termini of both H3 and H4 are cooperatively involved in the repression of a-cell specific genes in $\alpha$ cells, possibly through their interaction with Tup1.^ GCN5 was originally identified as a transcriptional regulator required to activate a subset of genes in yeast. Recently, it has been shown that GCN5 encodes the catalytic subunit of a nuclear histone acetyltransferase, providing the first direct link between histone acetylation and gene transcription. Recombinant Gcn5p (rGcn5p) exhibits a limited substrate specificity in vitro. However, neither the specificity of this enzyme in vivo nor the importance of particular acetylated residues to transcription or cell growth are well defined. In order to define the sites of histone acetylation mediated by Gcn5p in vivo and assess the significance of histone acetylation, more than 30 yeast strains have been constructed to bear specific H3 and/or H4 mutations in the presence or absence of GCN5 function. Our genetic data suggest that Gcn5p may have additional targets in vivo that are not identified as the targets of rGcn5p by previous studies. Western analysis using antibodies specifically recognizing particular acetylated isoforms of H3 and H4 led us to conclude that Gcn5p is necessary for full acetylation of multiple sites in both H3 and H4 in vivo. Consistent with these observations, rGcn5p still acetylates histones H3 and H4 bearing mutations either in H3 K14 or H4 K8,16, sites previously identified as the targets of acetylation by rGcn5p in H3 and H4. Our data also demonstrated that Gcn5p-mediated acetylation events are important for normal progression of the cell cycle and for transcriptional activation. Furthermore, a critical overall level of acetylation is essential for cell viability. ^

Analysis of primary Wilms tumors for WT1 RNA processing alterations

Wilms tumor (WT) is an embryonal renal tumor with a heterogeneous genetic etiology that serves as a valuable model for studying tumorigenesis. Biallelic inactivation of the tumor suppressor gene WT1, a zinc-finger transcriptional regulator located at 11p13, is critical for the development of some Wilms tumors. Interestingly, WT1 genomic analysis has demonstrated mutations in less than 20% of WT cases. This suggests either other genes play a more major role in Wilms tumorigenesis or WT1 is functionally altered by mechanisms other than DNA mutation. Previous observations in rat and in WT xenograft cell lines have suggested that abnormal WT1 RNA processing (exon 6 RNA editing and aberrant exon 2 splicing, respectively) is a potential mechanism of altering WT1 function in the absence of a WT1 DNA mutation. However, the role of this abnormal RNA processing has not previously been assessed in primary Wilms tumors. ^ To test the hypothesis that abnormal WT1 RNA processing is a mechanism of WT1alteration during tumor development, WT1 RNA from 85 primary tumors was analyzed using reverse transcription and polymerase chain reaction amplification (RT-PCR). Although no evidence for WT1 RNA editing was observed, variable levels (5% to 50%) of aberrant WT1 exon 2 splicing were detected for 11 tumors in the absence of a detectable WT1 DNA mutation. Also, alteration of normal WT1 alternative splicing, observed as RNA isoform loss, was detected in five tumors with no apparent WT1 genomic alteration, although no consistent pattern of RNA isoform loss was detected. This abnormal WT1 splicing, detected by either loss of exon 2 from some of the transcripts or loss of RNA isoforms, is statistically correlated with relapse (p = 0.005). These studies demonstrate that abnormal WT1 RNA processing is not a common mechanism of abrogating normal WT1 function in primary tumors. However, in those cases in which abnormal WTI splicing is present, these data indicate that it may serve as a useful prognostic marker for relapse in WT patients. ^

IMMUNOLOGICAL AND CYTOLOGICAL INVESTIGATION OF GAMMA INTERFERON PRODUCTION IN HUMAN T-LYMPHOCYTES INDUCED BY PMA AND PHA (IMMUNOGOLD, CELL CYCLE ANALYSIS)

The present study examined cellular mechanisms involved in the production and secretion of human (gamma)IFN. The hypothesis of this investigation was that (gamma)IFN is an export glycoprotein whose synthesis in human T lymphocytes is dependent on membrane stimulation, polypeptide synthesis in the rough endoplasmic reticulum, packaging in the Golgi complex, and release from the cell by exocytosis.^ The model system for this examination utilized T lymphocytes from normal donors and patients with chronic lymphocytic leukemia (CLL) induced in vitro with the tumor promoter, phorbol 12-myristate 13-acetate (PMA) and the lectin, phytohemagglutinin (PHA) to produce (gamma)IFN. This study reconfirmed the ability of PMA and PHA to synergistically induce (gamma)IFN production in normal T lymphocytes, as measured by viral inhibition assays and radio-immunoassays for (gamma)IFN. The leukemic T cells were demonstrated to produce (gamma)IFN in response to treatment with PHA. PMA treatment also induced (gamma)IFN production in the leukemic T cells, which was much greater than that observed in similarly treated normal T cells. In these same cells, however, combined treatment of the agents was shown to be ineffective at inducing (gamma)IFN production beyond the levels stimulated by the individual agents. In addition, the present study reiterated the synergistic effect of PMA/PHA on the stimulation of growth kinetics in normal T cells. The cell cycle of the leukemic T cells was also responsive to treatment with the agents, particularly with PMA treatment. A number of morphological alterations were attributed to PMA treatment including the acquisition of an elongated configuration, nuclear folds, and large cytoplasmic vacuoles. Many of the effects were observed to be reversible with dilution of the agents, and reversion to this state occurred more rapidly in the leukemic T cells. Most importantly, utilization of a thin section immuno-colloidal gold labelling technique for electron microscopy provided, for the first time, direct evidence of the cellular mechanism of (gamma)IFN production and secretion. The results of this latter study support the idea that (gamma)IFN is produced in the rough endoplasmic reticulum, transferred to the Golgi complex for accumulation and packaging, and released from the T cells by exocytosis. ^

ANALYSIS OF THE CORYNEBACTERIUM PARVUM-INDUCED DECLINE OF MURINE NATURAL KILLER CELL-MEDIATED CYTOTOXICITY

Treatment of mice with the immunomodulating agent, Corynebacterium parvum (C. parvum), was shown to result in a severe and long-lasting depression of splenic natural killer (NK) cell-mediated cytotoxicity 5-21 days post-inoculation. Because NK cells have been implicated in immunosurveillance against malignancy (due to their spontaneous occurrence and rapid reactivity to a variety of histological types of tumors), as well as in resistance to established tumors, this decreased activity was of particular concern, since this effect is contrary to that which would be considered therapeutically desirable in cancer treatment (i.e. a potentiation of antitumor effector functions, including NK cell activity, would be expected to lead to a more effective destruction of malignant cells). Therefore, an analysis of the mechanism of this decline of splenic NK cell activity in C.parvum treated mice was undertaken.^ From in vitro co-culturing experiments, it was found that low NK-responsive C. parvum splenocytes were capable of reducing the normally high-reactivity of cells from untreated syngeneic mice to YAC-1 lymphoma, suggesting the presence of NK-directed suppressor cells in C. parvum treated animals. This was further supported by the demonstration of normal levels of cytotoxicity in C. parvum splenocyte preparations following Ficoll-Hypaque separation, which coincided with removal of the NK-suppressive capabilities of these cells. The T cell nature of these regulatory cells was indicated by (1) the failure of C. parvum to cause a reduction of NK cell activity, or the generation of NK-directed suppressor cells in T cell-deficient athymic mice, (2) the removal of C. parvum-induced suppression by T cell-depleting fractionation procedures or treatments, and (3) demonstration of suppression of NK cell activity by T cell-enriched C. parvum splenocytes. These studies suggest, therefore, that the eventual reduction of suppression by T cell elimination and/or inhibition, may result in a promotion of the antitumor effectiveness of C. parvum due to the contribution of "freed" NK effector cell activity.^ However, the temporary suppression of NK cell activity induced by C. parvum (reactivity of treated mice returns to normal levels within 28 days after C. parvum injection), may in fact be favorable in some situations, e.g. in bone marrow transplantation cases, since NK cells have been suggested to play a role also in the process of bone marrow graft rejection.^ Therefore, the discriminate use of agents such as C. parvum may allow for the controlled regulation of NK cell activity suggested to be necessary for the optimalization of therapeutic regimens. ^

CYTOGENETIC AND DNA CONTENT ANALYSIS IN MULTIPLE MYELOMA

In this dissertation, the cytogenetic characteristics of bone marrow cells from 41 multiple myeloma patients were investigated. These cytogenetic data were correlated with the total DNA content as measured by flow cytometry. Both the cytogenetic information and DNA content were then correlated with clinical data to determine if diagnosis and prognosis of multiple myeloma could be improved.^ One hundred percent of the patients demonstrated abnormal chromosome numbers per metaphase. The average chromosome number per metaphase ranged from 42 to 49.9, with a mean of 44.99. The percent hypodiploidy ranged from 0-100% and the percent hyperdiploidy from 0-53%. Detailed cytogenetic analyses were very difficult to perform because of the paucity of mitotic figures and the poor chromosome morphology. Thus, detailed chromosome banding analysis on these patients was impossible.^ Thirty seven percent of the patients had normal total DNA content, whereas 63% had abnormal amounts of DNA (one patient with less than normal amounts and 25 patients with greater than normal amounts of DNA).^ Several clinical parameters were used in the statistical analyses: tumor burden, patient status at biopsy, patient response status, past therapy, type of treatment and percent plasma cells. Only among these clinical parameters were any statistically significant correlations found: pretreatment tumor burden versus patient response, patient biopsy status versus patient response and past therapy versus patient response.^ No correlations were found between percent hypodiploid, diploid, hyperdiploid or DNA content, and the patient response status, nor were any found between those patients with: (a) normal plasma cells, low pretreatment tumor mass burden and more than 50% of the analyzed metaphases with 46 chromosomes; (b) normal amounts of DNA, low pretreatment tumor mass burden and more than 50% of the metaphases with 46 chromosomes; (c) normal amounts of DNA and normal quantities of plasma cells; (d) abnormal amounts of DNA, abnormal amounts of plasma cells, high pretreatment tumor mass burden and less than 50% of the metaphases with 46 chromosomes.^ Technical drawbacks of both cytogenetic and DNA content analysis in these multiple myeloma patients are discussed along with the lack of correlations between DNA content and chromosome number. Refined chromosome banding analysis awaits technical improvements before we can understand which chromosome material (if any) makes up the "extra" amounts of DNA in these patients. None of the correlations tested can be used as diagnostic or prognostic aids for multiple myeloma. ^

POU domain transcription factor Brn3b is critical for the normal development and survival of retinal ganglion cells

The POU domain transcription factor Brn3b/POU4F2 plays a critical role regulating gene expression in mouse retinal ganglion cells (RGCs). Previous investigations have shown that Brn3b is not required for initial cell fate specification or migration; however, it is essential for normal RGC differentiation. In contrast to wild type axons, the mutant neurites were phenotypically different: shorter, rougher, disorganized, and poorly fasciculated. Wild type axons stained intensely with axon specific marker tau-1, while mutant projections were weakly stained and the mutant projections showed strong labeling with dendrite specific marker MAP2. Brn-3b mutant axonal projections contained more microtubules and fewer neurofilaments, a dendritic characteristic, than the wild type. The mutant neurites also exhibited significantly weaker staining of neurofilament low-molecular-weight (NF-L) in the axon when compared to the wild type, and NF-L accumulation in the neuron cell body. The absence of Brn-3b results in an inability to form normal axons and enhanced apoptosis in RGCs, suggesting that Brn-3b may control a set of genes involved in axon formation. ^ Brn3b contains several distinct sequence motifs: a glycine/serine rich region, two histidine rich regions, and a fifteen amino acid conserved sequence shared by all Brn3 family members in the N-terminus and a POU specific and POU homeodomain in the C-terminus. Brn3b activates a Luciferase reporter over 25 fold in cell culture when binding to native brn3 binding sites upstream of a minimal promoter. When fused to the Gal4 DNA Binding domain (DBD) and driven by either a strong (CMV) or weaker (pAHD) promoter, the N-terminal of Brn3b is capable of similar activation when binding to Gal4 UAS sites, indicating a presumptive activator of transcription. Both full length Brn3b or the C-terminus fused to the Gal4DBD and driven by pCMV repressed a Luciferase reporter downstream of UAS binding sites. Lower levels of expression of the fusion protein driven by pADH resulted in an alleviation of repression. This repression appears to be a limitation of this system of transcriptional analysis and a potential pitfall in conventional pCMV based transfection assays. ^

Bayesian generalized linear models for meta-analysis of diagnostic tests

With the recognition of the importance of evidence-based medicine, there is an emerging need for methods to systematically synthesize available data. Specifically, methods to provide accurate estimates of test characteristics for diagnostic tests are needed to help physicians make better clinical decisions. To provide more flexible approaches for meta-analysis of diagnostic tests, we developed three Bayesian generalized linear models. Two of these models, a bivariate normal and a binomial model, analyzed pairs of sensitivity and specificity values while incorporating the correlation between these two outcome variables. Noninformative independent uniform priors were used for the variance of sensitivity, specificity and correlation. We also applied an inverse Wishart prior to check the sensitivity of the results. The third model was a multinomial model where the test results were modeled as multinomial random variables. All three models can include specific imaging techniques as covariates in order to compare performance. Vague normal priors were assigned to the coefficients of the covariates. The computations were carried out using the 'Bayesian inference using Gibbs sampling' implementation of Markov chain Monte Carlo techniques. We investigated the properties of the three proposed models through extensive simulation studies. We also applied these models to a previously published meta-analysis dataset on cervical cancer as well as to an unpublished melanoma dataset. In general, our findings show that the point estimates of sensitivity and specificity were consistent among Bayesian and frequentist bivariate normal and binomial models. However, in the simulation studies, the estimates of the correlation coefficient from Bayesian bivariate models are not as good as those obtained from frequentist estimation regardless of which prior distribution was used for the covariance matrix. The Bayesian multinomial model consistently underestimated the sensitivity and specificity regardless of the sample size and correlation coefficient. In conclusion, the Bayesian bivariate binomial model provides the most flexible framework for future applications because of its following strengths: (1) it facilitates direct comparison between different tests; (2) it captures the variability in both sensitivity and specificity simultaneously as well as the intercorrelation between the two; and (3) it can be directly applied to sparse data without ad hoc correction. ^

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

A cross-sectional analysis of physical activity and obesity in inner city Houston Texas Mexican American youth

Obesity prevalence among children and adolescents is rising. It is one of the most attributable causes of hospitalization and death. Overweight and obese children are more likely to suffer from associated conditions such as hypertension, dyslipidemia, chronic inflammation, increased blood clotting tendency, endothelial dysfunction, hyperinsulinemia, and asthma. These children and adolescents are also more likely to be overweight and obese in adulthood. Interestingly, rates of obesity and overweight are not evenly distributed across racial and ethnic groups. Mexican American youth have higher rates of obesity and are at higher risk of becoming obese than non-Hispanic black and non-Hispanic white children. ^ Methods. This cross-sectional study describes the association between rates of obesity and physical activity in a sample of 1313 inner-city Mexican American children and adolescents (5-19 years of age) in Houston, Texas. This study is important because it will contribute to our understanding of childhood and adolescent obesity in this at-risk population. ^ Data from the Mexican American Feasibility Cohort using the Mano a Mano questionnaire are used to describe this population's status of obesity and physical activity. An initial sample taken from 5000 households in inner city Houston Texas was used as the baseline for this prospective cohort. The questionnaire was given in person to the participants to complete (or to parents for younger children) at a home visit by two specially trained bilingual interviewers. Analysis comprised prevalence estimates of obesity represented as percentile rank (<85%= normal weight, >85%= at risk, >95%= obese) by age and gender. The association between light, moderate, strenuous activity, and obesity was also examined using linear regression. ^ Results. Overall, 46% of this Mexican American Feasibility cohort is overweight or obese. The prevalence for children in the 6-11 age range (53.2%) was significantly greater than that reported from NHANES, 1999–2002 data (39.4%). Although the percentage of overweight and obese among the 12-19 year olds was greater than that reported in NHANES (38.5% versus 38.6%) this difference was not statistically significant. ^ A significant association between BMI and sit time and moderate physical activity (both p < 0.05) found in this sample. For males, this association was significant for moderate physical activity (p < 0.01). For the females, this association was significant for BMI and sit time (p < 0.05). These results need to be interpreted in the light of design and measurement limitations. ^ Conclusion. This study supports observations that the inner city Houston Texas Mexican American child and adolescent population is more overweight and obese than nationally reported figures, and that there are positive relationships between BMI, activity levels, and sit time in this population. This study supports the need for public health initiatives within the Houston Hispanic community. ^

A meta-analysis: Obesity and colorectal cancer screening

Of cancer death, colorectal cancer death ranks second in the United States. Obesity is an important risk factor for colorectal cancer (1). Early detection of colorectal cancer when it is localized can effectively reduce mortality of colorectal cancer and increase survival time of patients if they are treated. Also, previous studies showed that obese women were more likely to delay breast cancer screening and cervical cancer screening than normal weight women (2-5). However, results from prior studies demonstrating the relationship between obesity and colorectal cancer screening are not consistent. This research was done to conduct a meta-analysis of previous cross-sectional studies selected from the Medline database and to evaluate the association between obesity and colorectal cancer screening. While the odds ratio was not statistically different from one, the results from this meta-analysis under the random effects model showed that obese people are slightly less likely to have colorectal cancer screening compared to normal weight individuals (OR,0.93;95% CI 0.75-1.15). This meta-analysis was particularly sensitive to one individual study (6) and the effect of obesity on colorectal cancer screening was statistically significant (OR, 0.87; 95% CI, 0.81-0.92) after removing Heo's study. Further systematic studies focused on whether the effect of obesity on colorectal cancer screening is limited to women only are suggested. ^

Identification and analysis of novel mitotic inhibitors of the Caenorhabditis elegans Aurora B kinase

Proper execution of mitosis requires the accurate segregation of replicated DNA into each daughter cell. The highly conserved mitotic kinase AIR-2/Aurora B is a dynamic protein that interacts with subsets of cofactors and substrates to coordinate chromosome segregation and cytokinesis in Caenorhabdiris elegans. To identify components of the AIR-2 regulatory pathway, a genome-wide RNAi-based screen for suppressors of air-2 temperature-sensitive mutant lethality was conducted. Here, I present evidence that two classes of suppressors identified in this screen are bona fide regulators of the AIR-2 kinase. The strongest suppressor cdc-48.3, encodes an Afg2/Spaf-related Cdc48-like AAA+ ATPase that regulates AIR-2 kinase activity and stability during C. elegans embryogenesis. Loss of CDC-48.3 suppresses the lethality of air-2 mutant embryos, marked by the restoration of the dynamic behavior of AIR-2 and rescue of chromosome segregation and cytokinesis defects. Loss of CDC-48.3 leads to mitotic delays and abnormal accumulation of AIR-2 during late telophase/mitotic exit. In addition, AIR-2 kinase activity is significantly upregulated from metaphase through mitotic exit in CDC-48.3 depleted embryos. Inhibition of the AIR-2 kinase is dependent on (1) a direct physical interaction between CDC-48.3 and AIR-2, and (2) CDC-48.3 ATPase activity. Importantly, the increase in AIR-2 kinase activity does not correlate with the stabilization of AIR-2 in late mitosis. Hence, CDC-48.3 is a bi-functional inhibitor of AIR-2 that is likely to act via distinct mechanisms. The second class of suppressors consists of psy-2/smk-1 and pph-4.1, which encode two components of the conserved PP4 phosphatase complex that is essential for spindle assembly, chromosome segregation, and overall mitotic progression. AIR-2 and its substrates are likely to be targets of this complex since mitotic AIR-2 kinase activity is significantly increased during mitosis when either PSY-2/SMK-1 or PPH-4.l is depleted. Altogether, this study demonstrates that during the C. elegans embryonic cell cycle, regulators including the CDC-48.3 ATPase and PP4 phosphatase complex interact with and control the kinase activity, targeting behavior and protein stability of the Aurora B kinase to ensure accurate and timely progression of mitosis. ^

Transcriptional regulation of lung tumor suppressor GPRC5A in malignant and normal cells

Chronic exposure of the airways to cigarette smoke induces inflammatory response and genomic instability that play important roles in lung cancer development. Nuclear factor kappa B (NF-κB), the major intracellular mediator of inflammatory signals, is frequently activated in preneoplastic and malignant lung lesions. ^ Previously, we had shown that a lung tumor suppressor GPRC5A is frequently repressed in human non-small cell lung cancers (NSCLC) cells and lung tumor specimens. Recently, other groups have shown that human GPRC5A transcript levels are higher in bronchial samples of former than of current smokers. These results suggested that smoking represses GPRC5A expression and thus promotes the occurrence of lung cancer. We hypothesized that cigarette smoking or associated inflammatory response repressed GPRC5A expression through NF-κB signaling. ^ To determine the effect of inflammation, we examined GPRC5A protein expression in several lung cell lines following by TNF-α treatment. TNF-α significantly suppressed GPRC5A expression in normal small airway epithelial cells (SAEC) as well as in Calu-1 cells. Real-time PCR analysis indicated that TNF-α inhibits GPRC5A expression at the transcriptional level. NF-κB, the major downstream effectors of TNF-α signaling, mediates TNF-α-induced repression of GPRC5A because over-expression of NF-κB suppressed GPRC5A. To determine the region in the GPRC5A promoter through which NF-κB acts, we examined the ability of TNF-α to inhibit a series of reporter constructs with different deletions of GPRC5A promoter. The luciferase assay showed that the potential NF-κB binding sites containing region are irresponsible for TNF-α-induced suppression. Further analysis using constructs with different deletions in p65 revealed that NF-κB-mediated repression of GPRC5A is transcription-independent. Co-immunoprecipitation assays revealed that NF-κB could form a complex with RAR/RXR heterodimer. Moreover, the inhibitory effect of NF-κB has been found to be proportional to NF-κB/RAR ratio in luciferase assay. Finally, Chromatin IP demonstrated that NF-κB/p65 bound to GPRC5A promoter as well as RAR/RXR and suppressed transcription. Taken together, we propose that inflammation-induced NF-κB activation disrupts the RA signaling and suppresses GPRC5A expression and thus contributes to the oncogenesis of lung cancer. Our studies shed new light on the pathogenesis of lung cancer and potentially provide novel interventions for preventing and treating this disease. ^

DEVELOPMENT OF A BEAM-SPECIFIC PLANNING TARGET VOLUME AND A ROBUST PLAN ANALYSIS TOOL FOR PROTON THERAPY

Proton therapy is growing increasingly popular due to its superior dose characteristics compared to conventional photon therapy. Protons travel a finite range in the patient body and stop, thereby delivering no dose beyond their range. However, because the range of a proton beam is heavily dependent on the tissue density along its beam path, uncertainties in patient setup position and inherent range calculation can degrade thedose distribution significantly. Despite these challenges that are unique to proton therapy, current management of the uncertainties during treatment planning of proton therapy has been similar to that of conventional photon therapy. The goal of this dissertation research was to develop a treatment planning method and a planevaluation method that address proton-specific issues regarding setup and range uncertainties. Treatment plan designing method adapted to proton therapy: Currently, for proton therapy using a scanning beam delivery system, setup uncertainties are largely accounted for by geometrically expanding a clinical target volume (CTV) to a planning target volume (PTV). However, a PTV alone cannot adequately account for range uncertainties coupled to misaligned patient anatomy in the beam path since it does not account for the change in tissue density. In order to remedy this problem, we proposed a beam-specific PTV (bsPTV) that accounts for the change in tissue density along the beam path due to the uncertainties. Our proposed method was successfully implemented, and its superiority over the conventional PTV was shown through a controlled experiment.. Furthermore, we have shown that the bsPTV concept can be incorporated into beam angle optimization for better target coverage and normal tissue sparing for a selected lung cancer patient. Treatment plan evaluation method adapted to proton therapy: The dose-volume histogram of the clinical target volume (CTV) or any other volumes of interest at the time of planning does not represent the most probable dosimetric outcome of a given plan as it does not include the uncertainties mentioned earlier. Currently, the PTV is used as a surrogate of the CTV’s worst case scenario for target dose estimation. However, because proton dose distributions are subject to change under these uncertainties, the validity of the PTV analysis method is questionable. In order to remedy this problem, we proposed the use of statistical parameters to quantify uncertainties on both the dose-volume histogram and dose distribution directly. The robust plan analysis tool was successfully implemented to compute both the expectation value and its standard deviation of dosimetric parameters of a treatment plan under the uncertainties. For 15 lung cancer patients, the proposed method was used to quantify the dosimetric difference between the nominal situation and its expected value under the uncertainties.

DEVELOPMENT OF NOVEL METHODS TO MINIMIZE THE IMPACT OF SEQUENCING ERRORS IN THE NEXT-GENERATION SEQUENCING DATA ANALYSIS

Next-generation sequencing (NGS) technology has become a prominent tool in biological and biomedical research. However, NGS data analysis, such as de novo assembly, mapping and variants detection is far from maturity, and the high sequencing error-rate is one of the major problems. . To minimize the impact of sequencing errors, we developed a highly robust and efficient method, MTM, to correct the errors in NGS reads. We demonstrated the effectiveness of MTM on both single-cell data with highly non-uniform coverage and normal data with uniformly high coverage, reflecting that MTM’s performance does not rely on the coverage of the sequencing reads. MTM was also compared with Hammer and Quake, the best methods for correcting non-uniform and uniform data respectively. For non-uniform data, MTM outperformed both Hammer and Quake. For uniform data, MTM showed better performance than Quake and comparable results to Hammer. By making better error correction with MTM, the quality of downstream analysis, such as mapping and SNP detection, was improved. SNP calling is a major application of NGS technologies. However, the existence of sequencing errors complicates this process, especially for the low coverage (