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

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

Low birthweight and the risk of developing subsequent high blood pressure and/or dyslipidemia in childhood

The association between birthweight and blood pressure (BP), and birthweight and serum lipid concentrations at age 7 through 11 years was examined in 1446 black and white children. The prevalence ratio (with 95% confidence interval) for being in the race-, sex- and age-specific upper decile of diastolic BP in children born with low birthweight (LBW, $<$2500 grams) versus children with birthweight $\geq$2500 grams was for black boys, 2.66 (1.24-5.70). In the other race-sex groups for diastolic BP, and in all race-sex groups for systolic BP this ratio did not differ from one. Among white boys with LBW, but not in the other race-sex groups, higher than expected percentages of subjects were in the highest decile group of triglyceride concentrations (0.01 $<$ p $<$ 0.05). The prevalence ratio was 2.42 (1.19-4.91). When prematures were excluded only more than expected white girls with LBW were in the highest decile group of triglyceride concentrations. The prevalence ratio was 3.23 (1.16-9.00). Prevalence ratios for triglyceride concentrations in black boys and girls, and for LDL/HDL-C ratio, cholesterol and VLDL-C concentrations in all race-sex groups were not different from one in analyses including and in those excluding prematures. Mean triglyceride concentrations stratified by tertiles of Quetelet Index, race and sex showed a strongly positive association between triglyceride concentrations and Quetelet Index, and in the upper tertile of the Quetelet Index an association between LBW and raised triglyceride concentrations. Multiple linear regression analyses showed that after adjusting for sex, race and age present Quetelet Index (p $<$ 0.001) is a much stronger predictor of systolic and diastolic BP, and also of LDL-C/HDL-C ratio and triglyceride concentrations in this age group than birthweight (p $>$ 0.05). Thus, an association between LBW and subsequent risk for elevated BP was confirmed for diastolic BP in black boys, but not for the other race-sex groups, and not for systolic BP in any group. This is the first study finding an association between LBW and elevated triglyceride concentrations in boys (white and black) and girls (white). A follow-up study to assess whether the findings can be confirmed at adult age is recommended. ^

The role of serotonin in depression and clotting in the coronary artery disease population

Background: The mechanisms underlying the relationship between depression and acute coronary syndrome (ACS) remain unclear. Platelet serotonin has been associated with both depression and coronary artery disease in stable outpatients. Understanding the association between depression and platelet serotonin, during ACS, may explain some of the acute cardiovascular events seen in some individuals with depression. ^ Objectives: This study was designed to evaluate whether levels of platelet serotonin, during ACS, differ between individuals who screen positive for depression and individuals who screen negative for depression and to determine if a dose-response relationship exists between depressive symptoms and platelet serotonin levels. ^ Methods: In this cross-sectional study, data was collected on 51 patients hospitalized for ACS. Multiple linear regression models were used to determine if a relationship exists between depression and platelet serotonin levels. ^ Results: Of the 51 ACS patients, 24 screened positive for depression and 27 screened negative for depression. Platelet serotonin levels were not significantly different between the depressed group (942.10 ± 461.3) and the non-depressed group (1192.41 ± 764.3) (p= .293 and β= -4.093) and a dose-response relationship between depressive symptoms and platelet serotonin levels was not found (p= .250 and β= -.254). ^ Discussion: In this study, a relationship between depression and platelet serotonin levels was not found. Future research should focus on gaining a better understanding of the variables that may influence platelet serotonin levels in the ACS population. ^

Identification of genetic variants that influence atherosclerosis in young persons

Atherosclerosis is widely accepted as a complex genetic phenotype and is the usual cause of cardiovascular disease, the world’s leading killer. Genetic factors have been proven to be important risk contributors for atherosclerosis and much work has been done to identify promising candidates that might play a role in the development of atherosclerosis. It is well known that many independent replications are needed to unequivocally establish a valid genotype-phenotype association across different populations before the findings are extended to clinical settings and to the expensive follow-up studies designed to identify causal genetic variants. Aiming to replicate the association with atherosclerosis in the Pathobiological Determinants of Atherosclerosis in Youth (PDAY) study, we assessed the relationship of 32 atherosclerosis candidate SNPs to atherosclerosis in the PDAY cohort, consisting of AA and EA young people aged 15-34 years who died of non-medical causes. Two association studies, a whole sample study and a 1:1 matched case control study were performed by use of multiple linear regression and logistic regression analyses, respectively. For the whole sample association study, 32 SNPs among 2,650 individuals (1,369 AA and 1,281 EA) were tested for the association with six early atherosclerosis phenotypes: abdominal aorta fatty streaks, abdominal aorta raised lesions, right coronary artery fatty streaks, right coronary artery raised lesions, thoracic aorta fatty streaks, and thoracic aorta raised lesions. For the matched case-control association study, 337 case-control paired samples were included; cases were chosen with the highest total raised lesion scores from the studied population, while controls were randomly selected from individuals that had no raised lesions and matched to cases by age, gender and race. Sixteen SNPs in 13 genes were found to be significantly associated with atherosclerosis in at least one of the PDAY association studies. Among these 16 findings: eight SNPs (rs9579646, rs6053733, rs3849150, rs10499903, rs2148079, rs5073691, rs10116277, and rs17228212) successfully replicated previous results, six SNPs (rs17222814, rs10811661, rs7028570, rs7291467, rs16996148 and rs10401969) were reported as new findings exclusive to our study, the last two of the 16 SNPs, rs501120 and rs6922269, showed either intriguing or conflicting result. SNP rs17222814 in ALOX5AP and SNP rs3849150 in LRRC18 were consistently associated with atherosclerosis in both prior and the two PDAY association studies. SNP rs3849150 was also identified to be highly correlated with a non-synonymous coding SNP, rs17772611, which may damage the protein (polyphen score = 0.996), suggesting that SNP rs17772611 may be the causal functional variant.^ In conclusion, our study added more support for the association of these candidate genes with atherosclerosis. SNPs rs3849150 and rs17772611 of LRRC18, as well as SNP rs17222814 of ALOX5AP, were the most significant findings from our study, and may be ranked among the best for further study.^

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

Evaluation of racial and ethnic disparities of glycemic control in the Health and Retirement Study. 2003 Diabetes Survey

Objective. To examine and evaluate racial and ethnic disparities in glycemic control among HRS respondents with diabetes aged 55-94 years. ^ Methods. The HRS Diabetes 2003 database provides data on blood-drawn glycemic control and self-reported demographics, socioeconomic status, clinical, health access and self-care characteristics. 1,141 non-Hispanic White, non-Hispanic Black, and Hispanic respondents were included in multiple logistic regression of glycemic control. ^ Results. The rate of poor control was significantly higher among Blacks (61.5%, 105/171) and Hispanics (65.3% 72/110) than among Whites (45.0% 387/860) (p < 0.01). After controlling for influential covariates and interactions, Blacks and Hispanics had a three-fold increased risk for poor control compared to Whites when duration was five years or less. ^ Conclusions. Clinical and self-perception variables, like duration, medication, and self-rated poor diabetes control affected glycemic control independent of race and ethnicity, but there remains unexplained racial and ethnic disparities for newly-diagnosed individuals. This is the first study to find an interaction between duration and race and ethnicity on glycemic control. Future research should incorporate cultural beliefs and attitudes about diabetes control that may explain the racial and ethnic disparity. ^

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

Dietary intake patterns and relationships to polybrominated diphenyl ether (PBDE) and phthalate body burden

Polybrominated diphenyl ethers (PBDEs) and phthalates are chemicals of concern because of high levels measured in people and the environment as well as the demonstrated toxicity in animal studies and limited epidemiological studies. Exposure to these chemicals has been associated with a range of toxicological outcomes, including developmental effects, behavioral changes, endocrine disruption, effects on sexual health, and cancer. Previous research has shown that both of these classes of chemicals contaminate food in the United States and worldwide. However, how large a role diet plays in exposure to these chemicals is currently unknown. To address this question, an exploratory analysis of data collected as part of the 2003-04 National Health and Nutrition Examination Survey (NHANES) was conducted. Associations between dietary intake (assessed by 24-hour dietary recalls) for a range of food types (meat, poultry, fish, and dairy) and levels PBDEs and phthalate metabolites were analyzed using multiple linear regression modeling. Levels of individual PBDE congeners 28, 47, 99, 100 as well as total PBDEs were found to be significantly associated with the consumption of poultry. Metabolites of di-(2-ethylhexyl) phthalate (DEHP) were found to be associated with the consumption of poultry, as well as with an increased consumption of fat of animal origin. These results, combined with results from previous studies, suggest that diet is an important route of intake for both PBDEs and phthalates. Further research needs to be conducted to determine the sources of food contamination with these toxic chemicals as well as to describe the levels of contamination of US food in a large, representative sample.^