A novel method for evaluating an interaction effect of correlated continuous covariates in a linear model

Interaction effect is an important scientific interest for many areas of research. Common approach for investigating the interaction effect of two continuous covariates on a response variable is through a cross-product term in multiple linear regression. In epidemiological studies, the two-way analysis of variance (ANOVA) type of method has also been utilized to examine the interaction effect by replacing the continuous covariates with their discretized levels. However, the implications of model assumptions of either approach have not been examined and the statistical validation has only focused on the general method, not specifically for the interaction effect.^ In this dissertation, we investigated the validity of both approaches based on the mathematical assumptions for non-skewed data. We showed that linear regression may not be an appropriate model when the interaction effect exists because it implies a highly skewed distribution for the response variable. We also showed that the normality and constant variance assumptions required by ANOVA are not satisfied in the model where the continuous covariates are replaced with their discretized levels. Therefore, naïve application of ANOVA method may lead to an incorrect conclusion. ^ Given the problems identified above, we proposed a novel method modifying from the traditional ANOVA approach to rigorously evaluate the interaction effect. The analytical expression of the interaction effect was derived based on the conditional distribution of the response variable given the discretized continuous covariates. A testing procedure that combines the p-values from each level of the discretized covariates was developed to test the overall significance of the interaction effect. According to the simulation study, the proposed method is more powerful then the least squares regression and the ANOVA method in detecting the interaction effect when data comes from a trivariate normal distribution. The proposed method was applied to a dataset from the National Institute of Neurological Disorders and Stroke (NINDS) tissue plasminogen activator (t-PA) stroke trial, and baseline age-by-weight interaction effect was found significant in predicting the change from baseline in NIHSS at Month-3 among patients received t-PA therapy.^

Correlates of tobacco cessation counseling among Hispanic physicians in the US: a cross-sectional survey study.

BACKGROUND: Physician advice is an important motivator for attempting to stop smoking. However, physicians' lack of intervention with smokers has only modestly improved in the last decade. Although the literature includes extensive research in the area of the smoking intervention practices of clinicians, few studies have focused on Hispanic physicians. The purpose of this study was to explore the correlates of tobacco cessation counseling practices among Hispanic physicians in the US. METHODS: Data were collected through a validated survey instrument among a cross-sectional sample of self-reported Hispanic physicians practicing in New Mexico, and who were members of the New Mexico Hispanic Medical Society in the year 2001. Domains of interest included counseling practices, self-efficacy, attitudes/responsibility, and knowledge/skills. Returned surveys were analyzed to obtain frequencies and descriptive statistics for each survey item. Other analyses included: bivariate Pearson's correlation, factorial ANOVAs, and multiple linear regressions. RESULTS: Respondents (n = 45) reported a low level of compliance with tobacco control guidelines and recommendations. Results indicate that physicians' familiarity with standard cessation protocols has a significant effect on their tobacco-related practices (r = .35, variance shared = 12%). Self-efficacy and gender were both significantly correlated to tobacco related practices (r = .42, variance shared = 17%). A significant correlation was also found between self-efficacy and knowledge/skills (r = .60, variance shared = 36%). Attitudes/responsibility was not significantly correlated with any of the other measures. CONCLUSION: More resources should be dedicated to training Hispanic physicians in tobacco intervention. Training may facilitate practice by increasing knowledge, developing skills and, ultimately, enhancing feelings of self-efficacy.

Understanding causes of hospitalization and access to care among newly diagnosed HIV patients in Houston, TX

More than a quarter of patients with HIV in the United States are diagnosed in hospital settings most often with advanced HIV related conditions.(1) There has been little research done on the causes of hospitalization when the patients are first diagnosed with HIV. The aim of this study was to determine if the patients are hospitalized due to an HIV related cause or due to some other co-morbidity. Reduced access to care could be one possible reason why patients are diagnosed late in the course of the disease. This study compared the access to care of patients diagnosed with HIV in hospital and outpatient setting. The data used for the study was a part of the ongoing study “Attitudes and Beliefs and Steps of HIV Care”. The participants in the study were newly diagnosed with HIV and recruited from both inpatient and outpatient settings. The primary and the secondary diagnoses from hospital discharge reports were extracted and a primary reason for hospitalization was ascertained. These were classified as HIV-related, other infectious causes, non–infectious causes, other systemic causes, and miscellaneous causes. Access to care was determined by a score based on responses to a set of questions derived from the HIV Cost and Services Utilization Study (HCSUS) on a 6 point scale. The mean score of the hospitalized patients and mean score of the patients diagnosed in an outpatient setting was compared. We used multiple linear regressions to compare mean differences in the two groups after adjusting for age, sex, race, household income educational level and health insurance at the time of diagnosis. There were 185 participants in the study, including 78 who were diagnosed in hospital settings and 107 who were diagnosed in outpatient settings. We found that HIV-related conditions were the leading cause of hospitalization, accounting for 60% of admissions, followed by non-infectious causes (20%) and then other infectious causes (17%). The inpatient diagnosed group did not have greater perceived access-to-care as compared to the outpatient group. Regression analysis demonstrated a statistically significant improvement in access-to-care with advancing education level (p=0.04) and with better health insurance (p=0.004). HIV-related causes account for many hospitalizations when patients are first diagnosed with HIV. Many of these HIV-related hospitalizations could have been prevented if patients were diagnosed early and linked to medical care. Programs to increase HIV awareness need to be an integral part of activities aimed at control of spread of HIV in the community. Routine testing for HIV infection to promote early HIV diagnosis can prevent significant morbidity and mortality.^

3D MPRAGE improves classification of cortical lesions in multiple sclerosis.

BACKGROUND: Gray matter lesions are known to be common in multiple sclerosis (MS) and are suspected to play an important role in disease progression and clinical disability. A combination of magnetic resonance imaging (MRI) techniques, double-inversion recovery (DIR), and phase-sensitive inversion recovery (PSIR), has been used for detection and classification of cortical lesions. This study shows that high-resolution three-dimensional (3D) magnetization-prepared rapid acquisition with gradient echo (MPRAGE) improves the classification of cortical lesions by allowing more accurate anatomic localization of lesion morphology. METHODS: 11 patients with MS with previously identified cortical lesions were scanned using DIR, PSIR, and 3D MPRAGE. Lesions were identified on DIR and PSIR and classified as purely intracortical or mixed. MPRAGE images were then examined, and lesions were re-classified based on the new information. RESULTS: The high signal-to-noise ratio, fine anatomic detail, and clear gray-white matter tissue contrast seen in the MPRAGE images provided superior delineation of lesion borders and surrounding gray-white matter junction, improving classification accuracy. 119 lesions were identified as either intracortical or mixed on DIR/PSIR. In 89 cases, MPRAGE confirmed the classification by DIR/PSIR. In 30 cases, MPRAGE overturned the original classification. CONCLUSION: Improved classification of cortical lesions was realized by inclusion of high-spatial resolution 3D MPRAGE. This sequence provides unique detail on lesion morphology that is necessary for accurate classification.

Deep gray matter atrophy in multiple sclerosis: a tensor based morphometry.

Tensor based morphometry (TBM) was applied to determine the atrophy of deep gray matter (DGM) structures in 88 relapsing multiple sclerosis (MS) patients. For group analysis of atrophy, an unbiased atlas was constructed from 20 normal brains. The MS brain images were co-registered with the unbiased atlas using a symmetric inverse consistent nonlinear registration. These studies demonstrate significant atrophy of thalamus, caudate nucleus, and putamen even at a modest clinical disability, as assessed by the expanded disability status score (EDSS). A significant correlation between atrophy and EDSS was observed for different DGM structures: (thalamus: r=-0.51, p=3.85 x 10(-7); caudate nucleus: r=-0.43, p=2.35 x 10(-5); putamen: r=-0.36, p=6.12 x 10(-6)). Atrophy of these structures also correlated with 1) T2 hyperintense lesion volumes (thalamus: r=-0.56, p=9.96 x 10(-9); caudate nucleus: r=-0.31, p=3.10 x 10(-3); putamen: r=-0.50, p=6.06 x 10(-7)), 2) T1 hypointense lesion volumes (thalamus: r=-0.61, p=2.29 x 10(-10); caudate nucleus: r=-0.35, p=9.51 x 10(-4); putamen: r=-0.43, p=3.51 x 10(-5)), and 3) normalized CSF volume (thalamus: r=-0.66, p=3.55 x 10(-12); caudate nucleus: r=-0.52, p=2.31 x 10(-7), and putamen: r=-0.66, r=2.13 x 10(-12)). More severe atrophy was observed mainly in thalamus at higher EDSS. These studies appear to suggest a link between the white matter damage and DGM atrophy in MS.

Proton magnetic resonance spectroscopy in multiple sclerosis.

Proton magnetic resonance spectroscopy ((1)H-MRS) provides tissue metabolic information in vivo. This article reviews the role of MRS-determined metabolic alterations in lesions, normal-appearing white matter, gray matter, and spinal cord in advancing our knowledge of pathologic changes in multiple sclerosis (MS). In addition, the role of MRS in objectively evaluating therapeutic efficacy is reviewed. This potential metabolic information makes MRS a unique tool to follow MS disease evolution, understand its pathogenesis, evaluate the disease severity, establish a prognosis, and objectively evaluate the efficacy of therapeutic interventions.

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

14-3-3zeta overexpression in multiple human cancers plays a critical role in cancer progression

14-3-3 is a family of highly conserved and ubiquitously expressed proteins in eukaryotic organisms. 14-3-3 isoforms bind in a phospho-serine/threonine-dependent manner to a host of proteins involved in essential cellular processes including cell cycle, signal transduction and apoptosis. We fortuitously discovered 14-3-3 zeta overexpression in many human primary cancers, such as breast, lung, and sarcoma, and in a majority of cancer cell lines. To determine 14-3-3 zeta involvement in breast cancer progression, we used immunohistochemical analysis to examine 14-3-3 zeta expression in human primary invasive breast carcinomas. High 14-3-3 zeta expression was significantly correlated with poor prognosis of breast cancer patients. Increased expression of 14-3-3 zeta was also significantly correlated with elevated PKB/Akt activation in patient samples. Thus, 14-3-3 zeta is a marker of poor prognosis in breast cancers. Furthermore, up-regulation of 14-3-3 zeta enhanced malignant transformation of cancer cells in vitro. ^ To determine the biological significance of 14-3-3 zeta in human cancers, small interfering RNAs (siRNA) were used to specifically block 14-3-3 zeta expression in cancer cells. 14-3-3 zeta siRNA inhibited cellular proliferation by inducing a G1 arrest associated with up-regulation of p27 KIP1 and p21CIP1 cyclin dependent kinase inhibitors. Reduced 14-3-3 zeta inhibited PKB/Akt activation while stimulating the p38 signaling pathway. Silencing 14-3-3 zeta expression also increased stress-induced apoptosis by caspase activation. Notably, 14-3-3 zeta siRNA inhibited transformation related properties of breast cancer cells in vitro and inhibited tumor progression of breast cancer cells in vivo. 14-3-3 zeta may be a key regulatory factor controlling multiple signaling pathways leading to tumor progression. ^ The data indicate 14-3-3 zeta is a major regulator of cell growth and apoptosis and may play a critical role in the development of multiple cancer types. Hence, blocking 14-3-3 zeta may be a promising therapeutic approach for numerous cancers. ^

CARMA3 serves as a critical mediator of NF-kappaB in multiple signal transduction pathways

The central dogma of molecular biology dictates that DNA is transcribed into RNA, which is later translated into protein. One of the early activators in this process is the transcription factor NF-κB. We have determined that an NF-κB inducer, CARMA3, is required for proper neural tube closure, similar to other NF-κB inducers. Using a genetic knockout of CARMA3, we demonstrated that it is required for Gαq-coupled GPCR-induced NF-κB activation. This is facilitated through a MAPK and IKK phosphorylation-independent mechanism, most likely by controlling NEMO-associated ubiquitination. We have also shown that CARMA3 is required for EGF and HRG-induced NF-κB activation. This activation requires the activity of both EGFR and HER2, as well as PKC. Again, we observed no defect in IKK phosphorylation, although we determined a clear defect in IKK activation. Finally, we have begun to determine the role of CARMA3 to both EGFR and HER2-induced tumorigenicity. By overexpressing a constitutive active mutant of HER2 in our CARMA3 WT and KO MEF cells, we have shown CARMA3 is important for HER2-driven soft agar colony growth. We have also shown that knockdown of endogenous CARMA3 in the EGFR-overexpressing A431 cell line abolishes EGF-induced NF-κB activation. These same cells have a dramatically reduced capacity to form colonies in soft agar as well. Using both mouse xenografts and a transgenic model of HER2-induced breast cancer, we have initiated studies which will help to determine the role of CARMA3 to in vivo tumorigenesis. Collectively, this work reveals novel roles for the CARMA3 protein in development, GPCR and EGFR/HER2 signaling. It also suggests that CARMA3 is involved in EGFR/HER2 mediated tumorigenesis, possibly indicating a novel therapeutic target for use in treatment of cancer. ^

The survival role and the potential therapeutic targeting of MET in multiple myeloma

Multiple myeloma (MM) is a debilitating and incurable B-cell malignancy. Previous studies have documented that the hepatocyte growth factor (HGF) plays a role in the pathobiology of MM. The receptor tyrosine kinase MET induced signaling initiates when its ligand HGF binds to the MET receptor. However, the direct importance of MET in MM has not been elucidated. The present work used three different but complementary approaches to reduce MET protein levels or its activity to demonstrate the importance of MET in MM. ^ In the first approach, MET transcript and protein levels were reduced by directly targeting the cellular MET transcripts using shRNA retroviral infection techniques. This direct reduction of MET mRNA leads to a reduction of MET protein levels, which caused an inhibition of growth and induction of cell death. ^ In the second approach, a global transcription inhibitor flavopiridol was used as a potential pharmacological tool to reduce MET levels. MET has a short half-life of 30 min for mRNA and 4 hours for protein; therefore using a RNA pol II inhibitor such as flavopiridol would be a viable option to reduce MET levels. When using flavopiridol in MM cell lines, there was a reduction of MET transcript and protein levels, which was associated with the induction of cell death. ^ Finally in the last strategy, MET kinase activity was suppressed by MP470, a small molecule inhibitor that binds to the ATP binding pocket in the kinase domain. At concentrations where phosphorylation of MET was inhibited there was induction of cell death in MM cell lines and primary cells from patients. In addition, in MM cell lines there was a decrease in phosphorylation of AKT (ser473) and caspase-9 (ser196); downstream of MET, suggesting that the mechanism of action for survival may be through these cascade of events. ^ Overall, this study provides a proof-of-principle that MET is important for the survival of MM cell lines as well as primary plasma cells obtained from patients. Therefore, targeting MET therapeutically may be a possible strategy to treat patients with this debilitating disease of MM. ^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^

Robust effect sizes and their confidence intervals for group difference between trajectories in hierarchical linear growth model

Hierarchical linear growth model (HLGM), as a flexible and powerful analytic method, has played an increased important role in psychology, public health and medical sciences in recent decades. Mostly, researchers who conduct HLGM are interested in the treatment effect on individual trajectories, which can be indicated by the cross-level interaction effects. However, the statistical hypothesis test for the effect of cross-level interaction in HLGM only show us whether there is a significant group difference in the average rate of change, rate of acceleration or higher polynomial effect; it fails to convey information about the magnitude of the difference between the group trajectories at specific time point. Thus, reporting and interpreting effect sizes have been increased emphases in HLGM in recent years, due to the limitations and increased criticisms for statistical hypothesis testing. However, most researchers fail to report these model-implied effect sizes for group trajectories comparison and their corresponding confidence intervals in HLGM analysis, since lack of appropriate and standard functions to estimate effect sizes associated with the model-implied difference between grouping trajectories in HLGM, and also lack of computing packages in the popular statistical software to automatically calculate them. ^ The present project is the first to establish the appropriate computing functions to assess the standard difference between grouping trajectories in HLGM. We proposed the two functions to estimate effect sizes on model-based grouping trajectories difference at specific time, we also suggested the robust effect sizes to reduce the bias of estimated effect sizes. Then, we applied the proposed functions to estimate the population effect sizes (d ) and robust effect sizes (du) on the cross-level interaction in HLGM by using the three simulated datasets, and also we compared the three methods of constructing confidence intervals around d and du recommended the best one for application. At the end, we constructed 95% confidence intervals with the suitable method for the effect sizes what we obtained with the three simulated datasets. ^ The effect sizes between grouping trajectories for the three simulated longitudinal datasets indicated that even though the statistical hypothesis test shows no significant difference between grouping trajectories, effect sizes between these grouping trajectories can still be large at some time points. Therefore, effect sizes between grouping trajectories in HLGM analysis provide us additional and meaningful information to assess group effect on individual trajectories. In addition, we also compared the three methods to construct 95% confident intervals around corresponding effect sizes in this project, which handled with the uncertainty of effect sizes to population parameter. We suggested the noncentral t-distribution based method when the assumptions held, and the bootstrap bias-corrected and accelerated method when the assumptions are not met.^

Impact of medications prescribed for treatment of attention-deficit hyperactivity disorder on physical growth in children and adolescents with HIV.

OBJECTIVE: To examine the relationships between physical growth and medications prescribed for symptoms of attention-deficit hyperactivity disorder in children with HIV. METHODS: Analysis of data from children with perinatally acquired HIV (N = 2251; age 3-19 years), with and without prescriptions for stimulant and nonstimulant medications used to treat attention-deficit hyperactivity disorder, in a long-term observational study. Height and weight measurements were transformed to z scores and compared across medication groups. Changes in z scores during a 2-year interval were compared using multiple linear regression models adjusting for selected covariates. RESULTS: Participants with (n = 215) and without (n = 2036) prescriptions were shorter than expected based on US age and gender norms (p < .001). Children without prescriptions weighed less at baseline than children in the general population (p < .001) but gained height and weight at a faster rate (p < .001). Children prescribed stimulants were similar to population norms in baseline weight; their height and weight growth velocities were comparable with the general population and children without prescriptions (for weight, p = .511 and .100, respectively). Children prescribed nonstimulants had the lowest baseline height but were similar to population norms in baseline weight. Their height and weight growth velocities were comparable with the general population but significantly slower than children without prescriptions (p = .01 and .02, respectively). CONCLUSION: The use of stimulants to treat symptoms of attention-deficit hyperactivity disorder does not significantly exacerbate the potential for growth delay in children with HIV and may afford opportunities for interventions that promote physical growth. Prospective studies are needed to confirm these findings.

Genes to blood pressure: The effects of variation in genes of the renin-angiotensin system on the hormonal and renal control of blood pressure

Objective. Essential hypertension affects 25% of the US adult population and is a leading contributor to morbidity and mortality. Because BP is a multifactorial phenotype that resists simple genetic analysis, intermediate phenotypes within the complex network of BP regulatory systems may be more accessible to genetic dissection. The Renin-Angiotensin System (RAS) is known to influence intermediate and long-term blood pressure regulation through alterations in vascular tone and renal sodium and fluid resorption. This dissertation examines associations between renin (REN), angiotensinogen (AGT), angiotensin-converting enzyme (ACE) and angiotensin II type 1 receptor (AT1) gene variation and interindividual differences in plasma hormone levels, renal hemodynamics, and BP homeostasis.^ Methods. A total of 150 unrelated men and 150 unrelated women, between 20.0 and 49.9 years of age and free of acute or chronic illness except for a history of hypertension (11 men and 7 women, all off medications), were studied after one week on a controlled sodium diet. RAS plasma hormone levels, renal hemodynamics and BP were determined prior to and during angiotensin II (Ang II) infusion. Individuals were genotyped by PCR for a variable number tandem repeat (VNTR) polymorphism in REN, and for the following restriction fragment length polymorphisms (RFLP): AGT M235T, ACE I/D, and AT1 A1166C. Associations between clinical measurements and allelic variation were examined using multiple linear regression statistical models.^ Results. Women homozygous for the AT1 1166C allele demonstrated higher intracellular levels of sodium (p = 0.044). Men homozygous for the AGT T235 allele demonstrated a blunted decrement in renal plasma flow in response to Ang II infusion (p = 0.0002). There were no significant associations between RAS gene variation and interindividual variation in RAS plasma hormone levels or BP.^ Conclusions. Rather than identifying new BP controlling genes or alleles, the study paradigm employed in this thesis (i.e., measured genes, controlled environments and interventions) may provide mechanistic insight into how candidate genes affect BP homeostasis. ^

Association between temperament and physical activity in a population of minority women

Physical activity has been, and remains, a significant public health issue. Thus, increasing physical activity has been identified as a top priority according to Healthy People 2010. Various behavioral variables have been associated with participation in physical activity, including the Type A behavior pattern (TABP). This study was a secondary data analysis of the Women On The Move pilot study data and examined the relationship between Type A behavior with physical activity. The study population consisted of fifty-six (56) adult minority women 40 years of age and above. The Thurstone Activity Scale was adapted for use in this study to measure TABP. Physical activity behavior was measured using an accelerometer (Computer Science Application, [CSA]) and a physical activity diary. All study questions were examined using multiple linear regression analysis. In all analyses age, household income, and level of education were entered as covariates. The results found no association with TABP and exercise or physical activity. More research involving a larger, more active study population is recommended in order to more precisely determine the relationship of TABP and physical activity. ^