A multiplicative model for standardizing vital rates without using external standards

Standardization is a common method for adjusting confounding factors when comparing two or more exposure category to assess excess risk. Arbitrary choice of standard population in standardization introduces selection bias due to healthy worker effect. Small sample in specific groups also poses problems in estimating relative risk and the statistical significance is problematic. As an alternative, statistical models were proposed to overcome such limitations and find adjusted rates. In this dissertation, a multiplicative model is considered to address the issues related to standardized index namely: Standardized Mortality Ratio (SMR) and Comparative Mortality Factor (CMF). The model provides an alternative to conventional standardized technique. Maximum likelihood estimates of parameters of the model are used to construct an index similar to the SMR for estimating relative risk of exposure groups under comparison. Parametric Bootstrap resampling method is used to evaluate the goodness of fit of the model, behavior of estimated parameters and variability in relative risk on generated sample. The model provides an alternative to both direct and indirect standardization method. ^

Multifactor effects and evidence of potential interaction between complement factor H Y402H and LOC387715 A69S in age-related macular degeneration.

BACKGROUND: Variants in the complement cascade genes and the LOC387715/HTRA1, have been widely reported to associate with age-related macular degeneration (AMD), the most common cause of visual impairment in industrialized countries. METHODS/PRINCIPAL FINDINGS: We investigated the association between the LOC387715 A69S and complement component C3 R102G risk alleles in the Finnish case-control material and found a significant association with both variants (OR 2.98, p = 3.75 x 10(-9); non-AMD controls and OR 2.79, p = 2.78 x 10(-19), blood donor controls and OR 1.83, p = 0.008; non-AMD controls and OR 1.39, p = 0.039; blood donor controls), respectively. Previously, we have shown a strong association between complement factor H (CFH) Y402H and AMD in the Finnish population. A carrier of at least one risk allele in each of the three susceptibility loci (LOC387715, C3, CFH) had an 18-fold risk of AMD when compared to a non-carrier homozygote in all three loci. A tentative gene-gene interaction between the two major AMD-associated loci, LOC387715 and CFH, was found in this study using a multiplicative (logistic regression) model, a synergy index (departure-from-additivity model) and the mutual information method (MI), suggesting that a common causative pathway may exist for these genes. Smoking (ever vs. never) exerted an extra risk for AMD, but somewhat surprisingly, only in connection with other factors such as sex and the C3 genotype. Population attributable risks (PAR) for the CFH, LOC387715 and C3 variants were 58.2%, 51.4% and 5.8%, respectively, the summary PAR for the three variants being 65.4%. CONCLUSIONS/SIGNIFICANCE: Evidence for gene-gene interaction between two major AMD associated loci CFH and LOC387715 was obtained using three methods, logistic regression, a synergy index and the mutual information (MI) index.

Eye Movement Measures of Cognitive Control in Children with Tourette Syndrome

Tourette Syndrome begins in childhood and is characterized by uncontrollable repetitive actions like neck craning or hopping and noises such as sniffing or chirping. Worst in early adolescence, these tics wax and wane in severity and occur in bouts unpredictably, often drawing unwanted attention from bystanders. Making matters worse, over half of children with Tourette Syndrome also suffer from comorbid, or concurrent, disorders such as attention deficit hyperactivity disorder (ADHD) and obsessive-compulsive disorder (OCD). These disorders introduce anxious thoughts, impulsivity, inattention, and mood variability that further disrupt children with Tourette Syndrome from focusing and performing well at school and home. Thus, deficits in the cognitive control functions of response inhibition, response generation, and working memory have long been ascribed to Tourette Syndrome. Yet, without considering the effect of medication, age, and comorbidity, this is a premature attribution. This study used an infrared eye tracking camera and various computer tasks requiring eye movement responses to evaluate response inhibition, response generation, and working memory in Tourette Syndrome. This study, the first to control for medication, age, and comorbidity, enrolled 39 unmedicated children with Tourette Syndrome and 29 typically developing peers aged 10-16 years who completed reflexive and voluntary eye movement tasks and diagnostic rating scales to assess symptom severities of Tourette Syndrome, ADHD, and OCD. Children with Tourette Syndrome and comorbid ADHD and/or OCD, but not children with Tourette Syndrome only, took longer to respond and made more errors and distracted eye movements compared to typically-developing children, displaying cognitive control deficits. However, increasing symptom severities of Tourette Syndrome, ADHD, and OCD correlated with one another. Thus, cognitive control deficits were not specific to Tourette Syndrome patients with comorbid conditions, but rather increase with increasing tic severity, suggesting that a majority of Tourette Syndrome patients, regardless of a clinical diagnosis of ADHD and/or OCD, have symptoms of cognitive control deficits at some level. Therefore, clinicians should evaluate and counsel all families of children with Tourette Syndrome, with or without currently diagnosed ADHD and/or OCD, about the functional ramifications of comorbid symptoms and that they may wax and wane with tic severity.

NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

Mechanisms of perceptual learning of depth discrimination in random dot stereograms.

Perceptual learning is a training induced improvement in performance. Mechanisms underlying the perceptual learning of depth discrimination in dynamic random dot stereograms were examined by assessing stereothresholds as a function of decorrelation. The inflection point of the decorrelation function was defined as the level of decorrelation corresponding to 1.4 times the threshold when decorrelation is 0%. In general, stereothresholds increased with increasing decorrelation. Following training, stereothresholds and standard errors of measurement decreased systematically for all tested decorrelation values. Post training decorrelation functions were reduced by a multiplicative constant (approximately 5), exhibiting changes in stereothresholds without changes in the inflection points. Disparity energy model simulations indicate that a post-training reduction in neuronal noise can sufficiently account for the perceptual learning effects. In two subjects, learning effects were retained over a period of six months, which may have application for training stereo deficient subjects.

Primary lung cancer after breast cancer: The role of radiation therapy and cigarette smoking

The magnitude of the interaction between cigarette smoking, radiation therapy, and primary lung cancer after breast cancer remains unresolved. This case control study further examines the main and joint effects of cigarette smoking and radiation therapy (XRT) among breast cancer patients who subsequently developed primary lung cancer, at The University of Texas M. D. Anderson Cancer Center (MDACC) in Houston, Texas. Cases (n = 280) were women diagnosed with primary lung cancer between 1955 and 1970, between 30–89 years of age, who had a prior history of breast cancer, and were U.S. residents. Controls (n = 300) were randomly selected from 37,000 breast cancer patients at MDACC and frequency matched to cases on age at diagnosis (in 5-year strata), ethnicity, year of breast cancer diagnosis (in 5-year strata), and had survived at least as long as the time interval for lung cancer diagnosis in the cases. Stratified analysis and unconditional logistic regression modeling were used to calculate the main and joint effects of cigarette smoking and radiation treatment on lung cancer risk. Medical record review yielded smoking information on 93% of cases and 84% of controls, and among cases 45% received XRT versus 44% of controls. Smoking increased the odds of lung cancer in women who did not receive XRT (OR = 6.0, 95%CI, 3.5–10.1) whereas XRT was not associated with increased odds (OR = 0.5, 95%CI, 0.2–1.1) in women who did not smoke. Overall the odds ratio for both XRT and smoking together compared with neither exposure was 9.00 (9 5% CI, 5.1–15.9). Similarly, when stratifying on laterality of the lung cancer in relation to the breast cancer, and when the time interval between breast and lung cancers was >10 years, there was an increased odds for both smoking and XRT together for lung cancers on the same side as the breast cancer (ipsilateral) (OR = 11.5, 95% CI, 4.9–27.8) and lung cancers on the opposite side of the breast cancer (contralateral) (OR= 9.6, 95% CI, 2.9–0.9). After 20 years the odds for the ipsilateral lung were even more pronounced (OR = 19.2, 95% CI, 4.2–88.4) compared to the contralateral lung (OR = 2.6, 95% CI, 0.2–2.1). In conclusion, smoking was a significant independent risk factor for lung cancer after breast cancer. Moreover, a greater than multiplicative effect was observed with smoking and XRT combined being especially evident after 10 years for both the ipsilateral and contralateral lung and after 20 years for the ipsilateral lung. ^

Survival analysis of longevity in siblings

Do siblings of centenarians tend to have longer life spans? To answer this question, life spans of 184 siblings for 42 centenarians have been evaluated. Two important questions have been addressed in analyzing the sibling data. First, a standard needs to be established, to which the life spans of 184 siblings are compared. In this report, an external reference population is constructed from the U.S. life tables. Its estimated mortality rates are treated as baseline hazards from which the relative mortality of the siblings are estimated. Second, the standard survival models which assume independent observations are invalid when correlation within family exists, underestimating the true variance. Methods that allow correlations are illustrated by three different methods. First, the cumulative relative excess mortality between siblings and their comparison group is calculated and used as an effective graphic tool, along with the Product Limit estimator of the survival function. The variance estimator of the cumulative relative excess mortality is adjusted for the potential within family correlation using Taylor linearization approach. Second, approaches that adjust for the inflated variance are examined. They are adjusted one-sample log-rank test using design effect originally proposed by Rao and Scott in the correlated binomial or Poisson distribution setting and the robust variance estimator derived from the log-likelihood function of a multiplicative model. Nether of these two approaches provide correlation estimate within families, but the comparison with the comparison with the standard remains valid under dependence. Last, using the frailty model concept, the multiplicative model, where the baseline hazards are known, is extended by adding a random frailty term that is based on the positive stable or the gamma distribution. Comparisons between the two frailty distributions are performed by simulation. Based on the results from various approaches, it is concluded that the siblings of centenarians had significant lower mortality rates as compared to their cohorts. The frailty models also indicate significant correlations between the life spans of the siblings. ^

Diabetes genes and risk of prostate cancer in the Atherosclerosis Risk in Communities study

Prostate cancer (PrCa) is a leading cause of morbidity and mortality, yet the etiology remains uncertain. Meta-analyses show that PrCa risk is reduced by 16% in men with type 2 diabetes (T2D), but the mechanism is unknown. Recent genome-wide association studies and meta-analyses have found single nucleotide polymorphisms (SNPs) that consistently predict T2D risk. We evaluated associations of incident PrCa with 14 T2D SNPs in the Atherosclerosis Risk in Communities (ARIC) study. From 1987-2000, there were 397 incident PrCa cases ascertained from state or local cancer registries among 6,642 men (1,560 blacks and 5,082 whites) aged 45-64 years at baseline. Genotypes were determined by TaqMan assay. Cox proportional hazards models were used to assess the association between PrCa and increasing number of T2D risk-raising alleles for individual SNPs and for genetic risk scores (GRS) comprised of the number of T2D risk-raising alleles across SNPs. Two-way gene-gene interactions were evaluated with likelihood ratio tests. Using additive genetic models, the T2D risk-raising allele was associated with significantly reduced risk of PrCa for IGF2BP2 rs4402960 (hazard ratio [HR]=0.79; P=0.07 among blacks only), SLC2A2 rs5400 (race-adjusted HR=0.85; P=0.05) and UCP2 rs660339 (race-adjusted HR=0.84; P=0.02), but significantly increased risk of PrCa for CAPN10 rs3792267 (race-adjusted HR=1.20; P=0.05). No other SNPs were associated with PrCa using an additive genetic model. However, at least one copy of the T2D risk-raising allele for TCF7L2 rs7903146 was associated with reduced PrCa risk using a dominant genetic model (race-adjusted HR=0.79; P=0.03). These results imply that the T2D-PrCa association may be partly due to shared genetic variation, but these results should be verified since multiple tests were performed. When the combined, additive effects of these SNPs were tested using a GRS, there was nearly a 10% reduction in risk of PrCa per T2D risk-raising allele (race-adjusted HR=0.92; P=0.02). SNPs in IGF2BP2, KCNJ11 and SLC2A2 were also involved in multiple synergistic gene-gene interactions on a multiplicative scale. In conclusion, it appears that the T2D-PrCa association may be due, in part, to common genetic variation. Further knowledge of T2D gene-PrCa mechanisms may improve understanding of PrCa etiology and may inform PrCa prevention and treatment.^

Intensive care noise and mean arterial blood pressure in ELBW neonates

Objective. Loud noises in neonatal intensive care units (NICUs) may impede growth and development for extremely low birthweight (ELBW, < 1000 grams) newborns. The objective of this study was to measure the association between NICU sound levels and ELBW neonates' arterial blood pressure to determine whether these newborns experience noise-induced stress. ^ Methods. Noise and arterial blood pressure recordings were collected for 9 ELBW neonates during the first week of life. Sound levels were measured inside the incubator, and each subject's arterial blood pressures were simultaneously recorded for 15 minutes (at 1 sec intervals). Time series cross-correlation functions were calculated for NICU noise and mean arterial blood pressure (MABP) recordings for each subject. The grand mean noise-MABP cross-correlation was calculated for all subjects and for lower and higher birthweight groups for comparison. ^ Results. The grand mean noise-MABP cross-correlation for all subjects was mostly negative (through 300 sec lag time) and nearly reached significance at the 95% level at 111 sec lag (mean r = -0.062). Lower birthweight newborns (454-709 g) experienced significant decreases in blood pressure with increasing NICU noise after 145 sec lag (peak r = -0.074). Higher birthweight newborns had an immediate negative correlation with NICU sound levels (at 3 sec lag, r = -0.071), but arterial blood pressures increased to a positive correlation with noise levels at 197 sec lag (r = 0.075). ^ Conclusions. ELBW newborns' arterial blood pressure was influenced by NICU noise levels during the first week of life. Lower birthweight newborns may have experienced an orienting reflex to NICU sounds. Higher birthweight newborns experienced an immediate orienting reflex to increasing sound levels, but arterial blood pressure increased approximately 3 minutes after increases in noise levels. Increases in arterial blood pressure following increased NICU sound levels may result from a stress response to noise. ^

Information bias and lifetime mortality risks of radiation-induced cancer

Additive and multiplicative models of relative risk were used to measure the effect of cancer misclassification and DS86 random errors on lifetime risk projections in the Life Span Study (LSS) of Hiroshima and Nagasaki atomic bomb survivors. The true number of cancer deaths in each stratum of the cancer mortality cross-classification was estimated using sufficient statistics from the EM algorithm. Average survivor doses in the strata were corrected for DS86 random error ($\sigma$ = 0.45) by use of reduction factors. Poisson regression was used to model the corrected and uncorrected mortality rates with covariates for age at-time-of-bombing, age at-time-of-death and gender. Excess risks were in good agreement with risks in RERF Report 11 (Part 2) and the BEIR-V report. Bias due to DS86 random error typically ranged from $-$15% to $-$30% for both sexes, and all sites and models. The total bias, including diagnostic misclassification, of excess risk of nonleukemia for exposure to 1 Sv from age 18 to 65 under the non-constant relative projection model was $-$37.1% for males and $-$23.3% for females. Total excess risks of leukemia under the relative projection model were biased $-$27.1% for males and $-$43.4% for females. Thus, nonleukemia risks for 1 Sv from ages 18 to 85 (DRREF = 2) increased from 1.91%/Sv to 2.68%/Sv among males and from 3.23%/Sv to 4.02%/Sv among females. Leukemia excess risks increased from 0.87%/Sv to 1.10%/Sv among males and from 0.73%/Sv to 1.04%/Sv among females. Bias was dependent on the gender, site, correction method, exposure profile and projection model considered. Future studies that use LSS data for U.S. nuclear workers may be downwardly biased if lifetime risk projections are not adjusted for random and systematic errors. (Supported by U.S. NRC Grant NRC-04-091-02.) ^

Genetic susceptibility, cigarette smoking, and wood dust exposure in lung cancer: Case control analyses

This case control study was conducted to assess the association between lung cancer risk, mutagen sensitivity (a marker of cancer susceptibility), and a putative lung carcinogen, wood dust exposure. There were 165 cases (98 African-Americans, 67 Mexican-Americans) with newly diagnosed, previously untreated lung cancer, and 239 controls, frequency-matched on age, sex, and ethnicity.^ Mutagen sensitivity ($\ge$1 break/cell) was associated with a statistically significant elevated risk for lung cancer (odds ratio (OR) = 4.1, 95% confidence limits (CL) = 2.3,7.2). Wood dust exposure was also a significant predictor of risk (OR = 2.8, 95% CL = 1.2,6.6) after controlling for smoking and mutagen sensitivity. When stratified by ethnicity, wood dust exposure was a significant risk factor for African-Americans (OR = 4.0, 95% CL = 1.4,11.5), but not for Mexican-Americans (OR = 1.5, 95% CL = 0.3,7.1). Stratified analysis suggested a greater than multiplicative interaction between wood dust exposure and both mutagen sensitivity and smoking.^ The cases had significantly more breaks on chromosomes 4 and 5 than the controls did with ORs of 4.9 (95% CL = 2.0, 11.7) and 3.9 (95% CL = 1.6, 9.3), respectively. Breaks at 4p14, 4q27, 4q31, 5q21-22, 5q31, and 5q33 were significantly more common in lung cancer patients than in controls. Lung cancer risk had a dose-response relationship with breaks on chromosomes 4 and 5. Cigarette smoking had a strong interaction with breaks on chromosomes 2, 4, and 5.^ In a molecular cytogenetic study, using chromosome painting and G-banding, we showed that: (1) the proportion of chromosome 5 abnormalities surviving as chromosome-type aberrations remained significantly higher in cells of lung cancer cases (14%) than in controls (5%) (P $<$ 0.001). However, no significant differences were detected in chromosome 4 abnormalities between cases and controls; (2) the proportion of chromosome 5q13-22 abnormalities was 5.3% in the cases and 0.7% in the controls (P $<$ 0.001). 5q13-22 regions represented 40% of all abnormalities on chromosome 5 in the cases and only 14% in the controls.^ This study suggests that mutagen sensitivity, wood dust exposure, and cigarette smoking were independent risk factors for lung cancer, and the susceptibility of particular chromosome loci to mutagenic damage may be a genetic marker for specific types of lung cancer. ^