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.

The craniosacral progression of muscle development influences the emergence of neuromuscular junction alterations in a severe murine model for spinal muscular atrophy.

AIMS As 4-day-old mice of the severe spinal muscular atrophy (SMA) model (dying at 5-8 days) display pronounced neuromuscular changes in the diaphragm but not the soleus muscle, we wanted to gain more insight into the relationship between muscle development and the emergence of pathological changes and additionally to analyse intercostal muscles which are affected in human SMA. METHODS Structures of muscle fibres and neuromuscular junctions (NMJs) of the diaphragm, intercostal and calf muscles of prenatal (E21) and postnatal (P0 and P4) healthy and SMA mice were analysed by light and transmission electron microscopy. NMJ innervation was studied by whole mount immunofluorescence in diaphragms of P4 mice. RESULTS During this period, the investigated muscles still show a significant neck-to-tail developmental gradient. The diaphragm and calf muscles are most and least advanced, respectively, with respect to muscle fibre fusion and differentiation. The number and depth of subsynaptic folds increases, and perisynaptic Schwann cells (PSCs) acquire a basal lamina on their outer surface. Subsynaptic folds are connected to an extensive network of tubules and beaded caveolae, reminiscent of the T system in adult muscle. Interestingly, intercostal muscles from P4 SMA mice show weaker pathological involvement (that is, vacuolization of PSCs and perineurial cells) than those previously described by us for the diaphragm, whereas calf muscles show no pathological changes. CONCLUSION SMA-related alterations appear to occur only when the muscles have reached a certain developmental maturity. Moreover, glial cells, in particular PSCs, play an important role in SMA pathogenesis.

Development of a high- versus low-pathogenicity model of the free-living amoeba Naegleria fowleri

Species in the genus Naegleria are free-living amoebae of the soil and warm fresh water. Although around 30 species have been recognized, Naegleria fowleri is the only one that causes primary amoebic meningoencephalitis (PAM) in humans. PAM is an acute and fast progressing disease affecting the central nervous system. Most of the patients die within 1-2 weeks of exposure to the infectious water source. The fact that N. fowleri causes such fast progressing and highly lethal infections has opened many questions regarding the relevant pathogenicity factors of the amoeba. In order to investigate the pathogenesis of N. fowleri under defined experimental conditions, we developed a novel high- versus low-pathogenicity model for this pathogen. We showed that the composition of the axenic growth media influenced growth behaviour and morphology, as well as in vitro cytotoxicity and in vivo pathogenicity of N. fowleri. Trophozoites maintained in Nelson's medium were highly pathogenic for mice, demonstrated rapid in vitro proliferation, characteristic expression of surface membrane vesicles and a small cell diameter, and killed target mouse fibroblasts by both contact-dependent and -independent destruction. In contrast, N. fowleri cultured in PYNFH medium exhibited a low pathogenicity, slower growth, increased cell size and contact-dependent target cell destruction. However, cultivation of the amoeba in PYNFH medium supplemented with liver hydrolysate (LH) resulted in trophozoites that were highly pathogenic in mice, and demonstrated an intermediate proliferation rate in vitro, diminished cell diameter and contact-dependent target cell destruction. Thus, in this model, the presence of LH resulted in increased proliferation of trophozoites in vitro and enhanced pathogenicity of N. fowleri in mice. However, neither in vitro cytotoxicity mechanisms nor the presence of membrane vesicles on the surface correlated with the pathologic potential of the amoeba. This indicated that the pathogenicity of N. fowleri remains a complex interaction between as-yet-unidentified cellular mechanisms.

CX3CR1+ cells facilitate the activation of CD4 T cells in the colonic lamina propria during antigen-driven colitis

Dendritic cells (DCs) and macrophages populate the intestinal lamina propria to initiate immune responses required for the maintenance of intestinal homeostasis. To investigate whether CX3CR1(+) phagocytes communicate with CD4 T cells during the development of transfer colitis, we established an antigen-driven colitis model induced by the adoptive transfer of DsRed OT-II cells in CX3CR1(GFP/+) × RAG(-/-) recipients challenged with Escherichia coli expressing ovalbumin (OVA) fused to a cyan fluorescent protein (CFP). After colonization of CX3CR1(GFP/+) × RAG(-/-) animals with red fluorescent E. coli pCherry-OVA, colonic CX3CR1(+) cells but not CD103(+) DCs phagocytosed E. coli pCherry-OVA. Degraded bacterial-derived antigens are transported by CD103(+) DCs to mesenteric lymph nodes (MLNs), where CD103(+) DCs prime naive T cells. In RAG(-/-) recipients reconstituted with OT II cells and gavaged with OVA-expressing E. coli, colonic CX3CR1(+) phagocytes are in close contact with CD4 T cells and presented bacterial-derived antigens to CD4 T cells to activate and expand effector T cells.

The coupled atmosphere–chemistry–ocean model SOCOL-MPIOM

The newly developed atmosphere–ocean-chemistry-climate model SOCOL-MPIOM is presented by demonstrating the influence of the interactive chemistry module on the climate state and the variability. Therefore, we compare pre-industrial control simulations with (CHEM) and without (NOCHEM) interactive chemistry. In general, the influence of the chemistry on the mean state and the variability is small and mainly restricted to the stratosphere and mesosphere. The largest differences are found for the atmospheric dynamics in the polar regions, with slightly stronger northern and southern winter polar vortices in CHEM. The strengthening of the vortex is related to larger stratospheric temperature gradients, which are attributed to a parametrization of the absorption of ozone and oxygen in the Lyman-alpha, Schumann–Runge, Hartley, and Higgins bands. This effect is parametrized in the version with interactive chemistry only. A second reason for the temperature differences between CHEM and NOCHEM is related to diurnal variations in the ozone concentrations in the higher atmosphere, which are missing in NOCHEM. Furthermore, stratospheric water vapour concentrations differ substantially between the two experiments, but their effect on the temperatures is small. In both setups, the simulated intensity and variability of the northern polar vortex is inside the range of present day observations. Sudden stratospheric warming events are well reproduced in terms of their frequency, but the distribution amongst the winter months is too uniform. Additionally, the performance of SOCOL-MPIOM under changing external forcings is assessed for the period 1600–2000 using an ensemble of simulations driven by a spectral solar forcing reconstruction. The amplitude of the reconstruction is large in comparison to other state-of-the-art reconstructions, providing an upper limit for the importance of the solar signal. In the pre-industrial period (1600–1850) the simulated surface temperature trends are in reasonable agreement with temperature reconstructions, although the multi-decadal variability is more pronounced. This enhanced variability can be attributed to the variability in the solar forcing. The simulated temperature reductions during the Maunder Minimum are in the lowest probability range of the proxy records. During the Dalton Minimum, when also volcanic forcing is an important driver of temperature variations, the agreement is better. In the industrial period from 1850 onward SOCOL-MPIOM overestimates the temperature increase in comparison to observational data sets. Sensitivity simulations show that this overestimation can be attributed to the increasing trend in the solar forcing reconstruction that is used in this study and an additional warming induced by the simulated ozone changes.

A Drosophila XPD model links cell cycle coordination with neuro-development and suggests links to cancer

XPD functions in transcription, DNA repair and in cell cycle control. Mutations in human XPD (also known as ERCC2) mainly cause three clinical phenotypes: xeroderma pigmentosum (XP), Cockayne syndrome (XP/CS) and trichothiodystrophy (TTD), and only XP patients have a high predisposition to developing cancer. Hence, we developed a fly model to obtain novel insights into the defects caused by individual hypomorphic alleles identified in human XP-D patients. This model revealed that the mutations that displayed the greatest in vivo UV sensitivity in Drosophila did not correlate with those that led to tumor formation in humans. Immunoprecipitations followed by targeted quantitative MS/MS analysis showed how different xpd mutations affected the formation or stability of different transcription factor IIH (TFIIH) subcomplexes. The XP mutants most clearly linked to high cancer risk, Xpd R683W and R601L, showed a reduced interaction with the core TFIIH and also an abnormal interaction with the Cdk-activating kinase (CAK) complex. Interestingly, these two XP alleles additionally displayed high levels of chromatin loss and free centrosomes during the rapid nuclear division phase of the Drosophila embryo. Finally, the xpd mutations showing defects in the coordination of cell cycle timing during the Drosophila embryonic divisions correlated with those human mutations that cause the neurodevelopmental abnormalities and developmental growth defects observed in XP/CS and TTD patients.

Development and characterization of a scaffold-free 3D spheroid model of iPSC-derived human cardiomyocytes

Purpose: Cardiomyocytes are terminally differentiated cells in the adult heart and ischemia and cardiotoxic compounds can lead to cell death and irreversible decline of cardiac function. As testing platforms, isolated organs and primary cells from rodents have been the standard in research and toxicology, but there is a need for better models that more faithfully recapitulate native human biology. Hence, a new in vitro model comprising the advantages of 3D cell culture and the availability of induced pluripotent stem cells (iPSC) from human origin was developed and characterized. Methods: Human cardiomyocytes (CMs) derived from induced pluripotent stem cells (iPSCs) were studied in standard 2D culture and as cardiac microtissues (MTs) formed in hanging drops. 2D cultures were examined using immunofluorescence microscopy and Western blotting while the cardiac MTs were subjected to immunofluorescence, contractility, and pharmacological investigations. Results: iPSC-derived CMs in 2D culture showed well-formed myofibrils, cell-cell contacts positive for connexin-43, and other typical cardiac proteins. The cells reacted to pro-hypertrophic growth factors with a substantial increase in myofibrils and sarcomeric proteins. In hanging drop cultures, iPSC-derived cardiomyocytes formed spheroidal MTs within 4 days showing a homogeneous tissue structure with well-developed myofibrils extending throughout the whole spheroid without a necrotic core. MTs showed spontaneous contractions for more than 4 weeks that were recorded by optical motion tracking, sensitive to temperature, and responsive to electrical pacing. Contractile pharmacology was tested with several agents known to modulate cardiac rate and viability. Calcium-transients underlay the contractile activity and were also responsive to electrical stimulation, caffeine-induced Ca2+-release, extracellular calcium levels. Conclusions: 3D culture using iPSC-derived human cardiomyocytes provides an organoid human-based cellular platform that is free of necrosis and recapitulates vital cardiac functionality, thereby providing new and promising relevant model for the evaluation and development of new therapies and detection of cardiotoxicity.

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.

Solvation-Driven Charge Transfer and Localization in Metal Complexes

In any physicochemical process in liquids, the dynamical response of the solvent to the solutes out of equilibrium plays a crucial role in the rates and products: the solvent molecules react to the changes in volume and electron density of the solutes to minimize the free energy of the solution, thus modulating the activation barriers and stabilizing (or destabilizing) intermediate states. In charge transfer (CT) processes in polar solvents, the response of the solvent always assists the formation of charge separation states by stabilizing the energy of the localized charges. A deep understanding of the solvation mechanisms and time scales is therefore essential for a correct description of any photochemical process in dense phase and for designing molecular devices based on photosensitizers with CT excited states. In the last two decades, with the advent of ultrafast time-resolved spectroscopies, microscopic models describing the relevant case of polar solvation (where both the solvent and the solute molecules have a permanent electric dipole and the mutual interaction is mainly dipole−dipole) have dramatically progressed. Regardless of the details of each model, they all assume that the effect of the electrostatic fields of the solvent molecules on the internal electronic dynamics of the solute are perturbative and that the solvent−solute coupling is mainly an electrostatic interaction between the constant permanent dipoles of the solute and the solvent molecules. This well-established picture has proven to quantitatively rationalize spectroscopic effects of environmental and electric dynamics (time-resolved Stokes shifts, inhomogeneous broadening, etc.). However, recent computational and experimental studies, including ours, have shown that further improvement is required. Indeed, in the last years we investigated several molecular complexes exhibiting photoexcited CT states, and we found that the current description of the formation and stabilization of CT states in an important group of molecules such as transition metal complexes is inaccurate. In particular, we proved that the solvent molecules are not just spectators of intramolecular electron density redistribution but significantly modulate it. Our results solicit further development of quantum mechanics computational methods to treat the solute and (at least) the closest solvent molecules including the nonperturbative treatment of the effects of local electrostatics and direct solvent−solute interactions to describe the dynamical changes of the solute excited states during the solvent response.

Development of a Tumour Growth Inhibition Model to Elucidate the Effects of Ritonavir on Intratumoural Metabolism and Anti-tumour Effect of Docetaxel in a Mouse Model for Hereditary Breast Cancer

In a mouse tumour model for hereditary breast cancer, we previously explored the anti-cancer effects of docetaxel, ritonavir and the combination of both and studied the effect of ritonavir on the intratumoural concentration of docetaxel. The objective of the current study was to apply pharmacokinetic (PK)-pharmacodynamic (PD) modelling on this previous study to further elucidate and quantify the effects of docetaxel when co-administered with ritonavir. PK models of docetaxel and ritonavir in plasma and in tumour were developed. The effect of ritonavir on docetaxel concentration in the systemic circulation of Cyp3a knock-out mice and in the implanted tumour (with inherent Cyp3a expression) was studied, respectively. Subsequently, we designed a tumour growth inhibition model that included the inhibitory effects of both docetaxel and ritonavir. Ritonavir decreased docetaxel systemic clearance with 8% (relative standard error 0.4%) in the co-treated group compared to that in the docetaxel only-treated group. The docetaxel concentration in tumour tissues was significantly increased by ritonavir with mean area under the concentration-time curve 2.5-fold higher when combined with ritonavir. Observed tumour volume profiles in mice could be properly described by the PK/PD model. In the co-treated group, the enhanced anti-tumour effect was mainly due to increased docetaxel tumour concentration; however, we demonstrated a small but significant anti-tumour effect of ritonavir addition (p value <0.001). In conclusion, we showed that the increased anti-tumour effect observed when docetaxel is combined with ritonavir is mainly caused by enhanced docetaxel tumour concentration and to a minor extent by a direct anti-tumour effect of ritonavir.

Development of a Novel Rabies Simulation Model for Application in a Non-endemic Environment.

Domestic dog rabies is an endemic disease in large parts of the developing world and also epidemic in previously free regions. For example, it continues to spread in eastern Indonesia and currently threatens adjacent rabies-free regions with high densities of free-roaming dogs, including remote northern Australia. Mathematical and simulation disease models are useful tools to provide insights on the most effective control strategies and to inform policy decisions. Existing rabies models typically focus on long-term control programs in endemic countries. However, simulation models describing the dog rabies incursion scenario in regions where rabies is still exotic are lacking. We here describe such a stochastic, spatially explicit rabies simulation model that is based on individual dog information collected in two remote regions in northern Australia. Illustrative simulations produced plausible results with epidemic characteristics expected for rabies outbreaks in disease free regions (mean R0 1.7, epidemic peak 97 days post-incursion, vaccination as the most effective response strategy). Systematic sensitivity analysis identified that model outcomes were most sensitive to seven of the 30 model parameters tested. This model is suitable for exploring rabies spread and control before an incursion in populations of largely free-roaming dogs that live close together with their owners. It can be used for ad-hoc contingency or response planning prior to and shortly after incursion of dog rabies in previously free regions. One challenge that remains is model parameterisation, particularly how dogs' roaming and contacts and biting behaviours change following a rabies incursion in a previously rabies free population.

Gene regulation during placenta development: Murine adenosine deaminase as a model to study placenta specific expression

The formation of the placenta is one of the first and most important developmental events that occur in early mammalian embryogenesis. Even given this importance of the placenta, the academic community has largely ignored studying gene regulation during the development and maturation of the placenta. For this reason, an in-depth study of gene regulation in the trophoblast layer of the placenta using murine Adenosine Deaminase (Ada) as a model system has been undertaken. It has been determined that Ada is highly expressed in the placenta and is critical for embryo development. Dr. Kellems' laboratory has previously described a 1.8 kb fragment of the Ada 5 ′ flanking region that is capable of directing trophoblast specific expression in a transgenic model system. Preliminary studies have demonstrated several critical portions of this fragment that are necessary for the correct tissue specific expression in the placenta. My first specific aim was to elucidate the trans factor binding to one of these sequences, the FP3. Through electromobility shift assays (EMSA), the 30 bp FP3 was narrowed to a 5 bp sequence which computer databases predicted bound to Acute Myeloid Leukemia 1 (AML-1). This was confirmed by supershift analysis. The functional importance of this binding was demonstrated by a transgenic approach. A significant difference in expression of the reporter in the placenta was seen when the 5 bp sequence was mutated. This finding is a novel use for the AML-1 transcription factor which is the DNA binding portion of the heterodimer Core Binding Protein (CBP). The 5′ 240 bp region has also been demonstrated to contain functionally significant sequence. Through EMSA assays and computer predictions, the area has been narrowed to two pertinent regions that are predicted to contain GATA binding motifs. ^

Development and preliminary validation of the cancer family impact scale for colorectal cancer

This study aimed to develop and validate The Cancer Family Impact Scale (CFIS), an instrument for use in studies investigating relationships among family factors and colorectal cancer (CRC) screening when family history is a risk factor. We used existing data to develop the measure from 1,285 participants (637 families) across the United States who were in the Johns Hopkins Colon Cancer Genetic Testing study. Participants were 94% white with an average age of 50.1 years, and 60% were women. None had a personal CRC history, and eighty percent had 1 FDR with CRC and 20% had more than one FDR with CRC. The study had three aims: (1) to identify the latent factors underlying the CFIS via exploratory factor analysis (EFA); (2) to confirm the findings of the EFA via confirmatory factor analysis (CFA); and (3) to assess the reliability of the scale via Cronbach's alpha. Exploratory analyses were performed on a split half of the sample, and the final model was confirmed on the other half. The EFA suggested the CFIS was an 18-item measure with 5 latent constructs: (1) NEGATIVE: negative effects of cancer on the family; (2) POSITIVE: positive effects of cancer on the family; (3) COMMUNICATE: how families communicate about cancer; (4) FLOW: how information about cancer is conveyed in families; and (5) NORM: how individuals react to family norms about cancer. CFA on the holdout sample showed the CFIS to have a reasonably good fit (Chi-square = 389.977, df = 122, RMSEA= 0.058 (.052-.065), CFI=.902, TLI=.877, GF1=.939). The overall reliability of the scale was α=0.65. The reliability of the subscales was: (1) NEGATIVE α = 0.682; (2) POSITIVE α = 0.686; (3) COMMUNICATE α = 0.723; (4) FLOW α = 0.467; and (5) NORM α = 0.732. ^ We concluded the CFIS to be a good measure with most fit levels over 0.90. The CFIS could be used to compare theoretically driven hypotheses about the pathways through which family factors could influence health behavior among unaffected individuals at risk due to family history, and also aid in the development and evaluation of cancer prevention interventions including a family component. ^