Development of a Bayesian Joint Logistic Model to Better Study the Association between Haplotypes and Disease

In 2011, there will be an estimated 1,596,670 new cancer cases and 571,950 cancer-related deaths in the US. With the ever-increasing applications of cancer genetics in epidemiology, there is great potential to identify genetic risk factors that would help identify individuals with increased genetic susceptibility to cancer, which could be used to develop interventions or targeted therapies that could hopefully reduce cancer risk and mortality. In this dissertation, I propose to develop a new statistical method to evaluate the role of haplotypes in cancer susceptibility and development. This model will be flexible enough to handle not only haplotypes of any size, but also a variety of covariates. I will then apply this method to three cancer-related data sets (Hodgkin Disease, Glioma, and Lung Cancer). I hypothesize that there is substantial improvement in the estimation of association between haplotypes and disease, with the use of a Bayesian mathematical method to infer haplotypes that uses prior information from known genetics sources. Analysis based on haplotypes using information from publically available genetic sources generally show increased odds ratios and smaller p-values in both the Hodgkin, Glioma, and Lung data sets. For instance, the Bayesian Joint Logistic Model (BJLM) inferred haplotype TC had a substantially higher estimated effect size (OR=12.16, 95% CI = 2.47-90.1 vs. 9.24, 95% CI = 1.81-47.2) and more significant p-value (0.00044 vs. 0.008) for Hodgkin Disease compared to a traditional logistic regression approach. Also, the effect sizes of haplotypes modeled with recessive genetic effects were higher (and had more significant p-values) when analyzed with the BJLM. Full genetic models with haplotype information developed with the BJLM resulted in significantly higher discriminatory power and a significantly higher Net Reclassification Index compared to those developed with haplo.stats for lung cancer. Future analysis for this work could be to incorporate the 1000 Genomes project, which offers a larger selection of SNPs can be incorporated into the information from known genetic sources as well. Other future analysis include testing non-binary outcomes, like the levels of biomarkers that are present in lung cancer (NNK), and extending this analysis to full GWAS studies.

Evidence of an association between sleep and levodopa-induced dyskinesia in an animal model of Parkinson's disease.

Levodopa-induced dyskinesia (LID) represents a major challenge for clinicians treating patients affected by Parkinson's disease (PD). Although levodopa is the most effective treatment for PD, the remodeling effects induced by disease progression and the pharmacologic treatment itself cause a narrowing of the therapeutic window because of the development of LID. Although animal models of PD provide strong evidence that striatal plasticity underlies the development of dyskinetic movements, the pathogenesis of LID is not entirely understood. In recent years, slow homeostatic adjustment of intrinsic excitability occurring during sleep has been considered fundamental for network stabilization by gradually modifying plasticity thresholds. So far, how sleep affects on LID has not been investigated. Therefore, we measured synaptic downscaling across sleep episodes in a parkinsonian animal model showing dyskinetic movements similar to LID. Our electrophysiological, molecular, and behavioral results are consistent with an impaired synaptic homeostasis during sleep in animals showing dyskinesia. Accordingly, sleep deprivation causes an anticipation and worsening of LID supporting a link between sleep and the development of LID.

Experimental infection of the pig with Mycobacterium ulcerans: a novel model for studying the pathogenesis of Buruli ulcer disease.

BACKGROUND Buruli ulcer (BU) is a slowly progressing, necrotising disease of the skin caused by infection with Mycobacterium ulcerans. Non-ulcerative manifestations are nodules, plaques and oedema, which may progress to ulceration of large parts of the skin. Histopathologically, BU is characterized by coagulative necrosis, fat cell ghosts, epidermal hyperplasia, clusters of extracellular acid fast bacilli (AFB) in the subcutaneous tissue and lack of major inflammatory infiltration. The mode of transmission of BU is not clear and there is only limited information on the early pathogenesis of the disease available. METHODOLOGY/PRINCIPAL FINDINGS For evaluating the potential of the pig as experimental infection model for BU, we infected pigs subcutaneously with different doses of M. ulcerans. The infected skin sites were excised 2.5 or 6.5 weeks after infection and processed for histopathological analysis. With doses of 2 × 10(7) and 2 × 10(6) colony forming units (CFU) we observed the development of nodular lesions that subsequently progressed to ulcerative or plaque-like lesions. At lower inoculation doses signs of infection found after 2.5 weeks had spontaneously resolved at 6.5 weeks. The observed macroscopic and histopathological changes closely resembled those found in M. ulcerans disease in humans. CONCLUSION/SIGNIFICANCE Our results demonstrate that the pig can be infected with M. ulcerans. Productive infection leads to the development of lesions that closely resemble human BU lesions. The pig infection model therefore has great potential for studying the early pathogenesis of BU and for the development of new therapeutic and prophylactic interventions.

Identification of a missense mutation in the bovine ATP2A1 gene in congenital pseudomyotonia of Chianina cattle: an animal model of human Brody disease

Congenital pseudomyotonia in Chianina cattle is a muscle function disorder very similar to that of Brody disease in humans. Mutations in the human ATP2A1 gene, encoding SERCA1, cause Brody myopathy. The analysis of the collected Chianina pedigree data suggested monogenic autosomal recessive inheritance and revealed that all 17 affected individuals traced back to a single founder. A deficiency of SERCA1 function in skeletal muscle of pseudomyotonia affected Chianina cattle was observed as SERCA1 activity in affected animals was decreased by about 70%. Linkage analysis showed that the mutation was located in the ATP2A1 gene region on BTA25 and subsequent mutation analysis of the ATP2A1 exons revealed a perfectly associated missense mutation in exon 6 (c.491G>A) leading to a p.Arg164His substitution. Arg164 represents a functionally important and strongly conserved residue of SERCA1. This study provides a suitable large animal model for human Brody disease.

Genome profiling of sterol synthesis shows convergent evolution in parasites and guides chemotherapeutic attack

Sterols are an essential class of lipids in eukaryotes, where they serve as structural components of membranes and play important roles as signaling molecules. Sterols are also of high pharmacological significance: cholesterol-lowering drugs are blockbusters in human health, and inhibitors of ergosterol biosynthesis are widely used as antifungals. Inhibitors of ergosterol synthesis are also being developed for Chagas's disease, caused by Trypanosoma cruzi. Here we develop an in silico pipeline to globally evaluate sterol metabolism and perform comparative genomics. We generate a library of hidden Markov model-based profiles for 42 sterol biosynthetic enzymes, which allows expressing the genomic makeup of a given species as a numerical vector. Hierarchical clustering of these vectors functionally groups eukaryote proteomes and reveals convergent evolution, in particular metabolic reduction in obligate endoparasites. We experimentally explore sterol metabolism by testing a set of sterol biosynthesis inhibitors against trypanosomatids, Plasmodium falciparum, Giardia, and mammalian cells, and by quantifying the expression levels of sterol biosynthetic genes during the different life stages of T. cruzi and Trypanosoma brucei. The phenotypic data correlate with genomic makeup for simvastatin, which showed activity against trypanosomatids. Other findings, such as the activity of terbinafine against Giardia, are not in agreement with the genotypic profile.

CD4 T cells are required for both development and maintenance of disease in a new mouse model of reversible colitis

Current therapies to treat inflammatory bowel diseases have limited efficacy, significant side effects, and often wane over time. Little is known about the cellular and molecular mechanisms operative in the process of mucosal healing from colitis. To study such events, we developed a new model of reversible colitis in which adoptive transfer of CD4(+)CD45RB(hi) T cells into Helicobacter typhlonius-colonized lymphopenic mice resulted in a rapid onset of colonic inflammation that was reversible through depletion of colitogenic T cells. Remission was associated with an improved clinical and histopathological score, reduced immune cell infiltration to the intestinal mucosa, altered intestinal gene expression profiles, regeneration of the colonic mucus layer, and the restoration of epithelial barrier integrity. Notably, colitogenic T cells were not only critical for induction of colitis but also for maintenance of disease. Depletion of colitogenic T cells resulted in a rapid drop in tumor necrosis factor α (TNFα) levels associated with reduced infiltration of inflammatory immune cells to sites of inflammation. Although neutralization of TNFα prevented the onset of colitis, anti-TNFα treatment of mice with established disease failed to resolve colonic inflammation. Collectively, this new model of reversible colitis provides an important research tool to study the dynamics of mucosal healing in chronic intestinal remitting-relapsing disorders.Mucosal Immunology advance online publication 16 September 2015; doi:10.1038/mi.2015.93.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Past and future evolution of Abies alba forests in Europe - comparison of a dynamic vegetation model with palaeo data and observations

Information on how species distributions and ecosystem services are impacted by anthropogenic climate change is important for adaptation planning. Palaeo data suggest that Abies alba formed forests under significantly warmer-than-present conditions in Europe and might be a native substitute for widespread drought-sensitive temperate and boreal tree species such as beech (Fagus sylvatica) and spruce (Picea abies) under future global warming conditions. Here, we combine pollen and macrofossil data, modern observations, and results from transient simulations with the LPX-Bern dynamic global vegetation model to assess past and future distributions of A. alba in Europe. LPX-Bern is forced with climate anomalies from a run over the past 21 000 years with the Community Earth System Model, modern climatology, and with 21st-century multimodel ensemble results for the high-emission RCP8.5 and the stringent mitigation RCP2.6 pathway. The simulated distribution for present climate encompasses the modern range of A. alba, with the model exceeding the present distribution in north-western and southern Europe. Mid-Holocene pollen data and model results agree for southern Europe, suggesting that at present, human impacts suppress the distribution in southern Europe. Pollen and model results both show range expansion starting during the Bølling–Allerød warm period, interrupted by the Younger Dryas cold, and resuming during the Holocene. The distribution of A. alba expands to the north-east in all future scenarios, whereas the potential (currently unrealized) range would be substantially reduced in southern Europe under RCP8.5. A. alba maintains its current range in central Europe despite competition by other thermophilous tree species. Our combined palaeoecological and model evidence suggest that A. alba may ensure important ecosystem services including stand and slope stability, infrastructure protection, and carbon sequestration under significantly warmer-than-present conditions in central Europe.

Sleep-Wake Cycle Dysfunction in the TgCRND8 Mouse Model of Alzheimer's Disease: From Early to Advanced Pathological Stages.

In addition to cognitive decline, individuals affected by Alzheimer's disease (AD) can experience important neuropsychiatric symptoms including sleep disturbances. We characterized the sleep-wake cycle in the TgCRND8 mouse model of AD, which overexpresses a mutant human form of amyloid precursor protein resulting in high levels of β-amyloid and plaque formation by 3 months of age. Polysomnographic recordings in freely-moving mice were conducted to study sleep-wake cycle architecture at 3, 7 and 11 months of age and corresponding levels of β-amyloid in brain regions regulating sleep-wake states were measured. At all ages, TgCRND8 mice showed increased wakefulness and reduced non-rapid eye movement (NREM) sleep during the resting and active phases. Increased wakefulness in TgCRND8 mice was accompanied by a shift in the waking power spectrum towards fast frequency oscillations in the beta (14-20 Hz) and low gamma range (20-50 Hz). Given the phenotype of hyperarousal observed in TgCRND8 mice, the role of noradrenergic transmission in the promotion of arousal, and previous work reporting an early disruption of the noradrenergic system in TgCRND8, we tested the effects of the alpha-1-adrenoreceptor antagonist, prazosin, on sleep-wake patterns in TgCRND8 and non-transgenic (NTg) mice. We found that a lower dose (2 mg/kg) of prazosin increased NREM sleep in NTg but not in TgCRND8 mice, whereas a higher dose (5 mg/kg) increased NREM sleep in both genotypes, suggesting altered sensitivity to noradrenergic blockade in TgCRND8 mice. Collectively our results demonstrate that amyloidosis in TgCRND8 mice is associated with sleep-wake cycle dysfunction, characterized by hyperarousal, validating this model as a tool towards understanding the relationship between β-amyloid overproduction and disrupted sleep-wake patterns in AD.

Loss of Nogo-A-expressing neurons in a rat model of Parkinson's disease

The myelin-associated protein Nogo-A is among the most potent neurite growth inhibitors in the adult CNS. Recently, Nogo-A expression was demonstrated in a number of neuronal subpopulations of the adult and developing CNS but at present, little is known about the expression of Nogo-A in the nigrostriatal system, a brain structure severely affected in Parkinson's disease (PD). The present study sought to characterize the expression pattern of Nogo-A immunoreactive (ir) cells in the adult ventral mesencephalon of control rats and in the 6-hydroxydopamine (6-OHDA) rat model of PD. Immunohistochemical analyses of normal adult rat brain showed a distinct expression of Nogo-A in the ventral mesencephalon, with the highest level in the substantia nigra pars compacta (SNc) where it co-localized with dopaminergic neurons. Analyses conducted 1week and 1 month after unilateral striatal injections of 6-OHDA disclosed a severe loss of the number of Nogo-A-ir cells in the SNc. Notably, at 1week after treatment, more dopaminergic neurons expressing Nogo-A were affected by the 6-OHDA toxicity than Nogo-A-negative dopaminergic neurons. However, at later time points more of the surviving dopaminergic neurons expressed Nogo-A. In the striatum, both small and large Nogo-A-positive cells were detected. The large cells were identified as cholinergic interneurons. Our results suggest yet unidentified functions of Nogo-A in the CNS beyond the inhibition of axonal regeneration and plasticity, and may indicate a role for Nogo-A in PD.

A multivariate frailty model for disease recurrences and survival

A multivariate frailty hazard model is developed for joint-modeling of three correlated time-to-event outcomes: (1) local recurrence, (2) distant recurrence, and (3) overall survival. The term frailty is introduced to model population heterogeneity. The dependence is modeled by conditioning on a shared frailty that is included in the three hazard functions. Independent variables can be included in the model as covariates. The Markov chain Monte Carlo methods are used to estimate the posterior distributions of model parameters. The algorithm used in present application is the hybrid Metropolis-Hastings algorithm, which simultaneously updates all parameters with evaluations of gradient of log posterior density. The performance of this approach is examined based on simulation studies using Exponential and Weibull distributions. We apply the proposed methods to a study of patients with soft tissue sarcoma, which motivated this research. Our results indicate that patients with chemotherapy had better overall survival with hazard ratio of 0.242 (95% CI: 0.094 - 0.564) and lower risk of distant recurrence with hazard ratio of 0.636 (95% CI: 0.487 - 0.860), but not significantly better in local recurrence with hazard ratio of 0.799 (95% CI: 0.575 - 1.054). The advantages and limitations of the proposed models, and future research directions are discussed. ^

Disease risk factors and humoral antibody response to periodontal-associated pathogens: A model to evaluate periodontal disease in older adults

Periodontal diseases (PD) are infectious, inflammatory, and tissue destructive events which affect the periodontal ligament that surround and support the teeth. Periodontal diseases are the major cause of tooth loss after age 35, with gingivitis and periodontitis affecting 75% of the adult population. A select group of bacterial organisms are associated with periodontal pathogenesis. There is a direct association between oral hygiene and prevention of PD. The importance of genetic differences and host immune response capabilities in determining host, susceptibility or resistance to PD has not been established. This study examined the risk factors and serum (humoral) immune response to periodontal diseased-associated pathogens in a 55 to 80+ year old South Texas study sample with PD. This study sample was described by: age, sex, ethnicity, the socioeconomic factors marital status, income and occupation, IgG, IgA, IgM immunoglobulin status, and the autoimmune response markers rheumatoid factor (RF) and antinuclear antibody (ANA). These variables were used to determine the risk factors associated with development of PD. Serum IgG, IgA, IgM antibodies to bacterial antigens provided evidence for disease exposure.^ A causal model for PD was constructed from associations for risk factors (ethnicity, marital status, income, and occupation) with dental exam and periodontitis. The multiple correlation between PD and ethnicity, income and dental exam was significant. Hispanics of low income were least likely to have had a dental exam in the last year and most likely to have PD. The etiologic agents for PD, as evidenced by elevated humoral antibody responses, were the Gram negative microorganisms Bacteroides gingivalis, serotypes FDC381 and SUNYaBA7A1-28, and Wolinella recta. Recommendation for a PD prevention and control program are provided. ^