Statistical mechanics of protein complexed and condensed DNA

In this thesis I treat various biophysical questions arising in the context of complexed / ”protein-packed” DNA and DNA in confined geometries (like in viruses or toroidal DNA condensates). Using diverse theoretical methods I consider the statistical mechanics as well as the dynamics of DNA under these conditions. In the first part of the thesis (chapter 2) I derive for the first time the single molecule ”equation of state”, i.e. the force-extension relation of a looped DNA (Eq. 2.94) by using the path integral formalism. Generalizing these results I show that the presence of elastic substructures like loops or deflections caused by anchoring boundary conditions (e.g. at the AFM tip or the mica substrate) gives rise to a significant renormalization of the apparent persistence length as extracted from single molecule experiments (Eqs. 2.39 and 2.98). As I show the experimentally observed apparent persistence length reduction by a factor of 10 or more is naturally explained by this theory. In chapter 3 I theoretically consider the thermal motion of nucleosomes along a DNA template. After an extensive analysis of available experimental data and theoretical modelling of two possible mechanisms I conclude that the ”corkscrew-motion” mechanism most consistently explains this biologically important process. In chapter 4 I demonstrate that DNA-spools (architectures in which DNA circumferentially winds on a cylindrical surface, or onto itself) show a remarkable ”kinetic inertness” that protects them from tension-induced disruption on experimentally and biologically relevant timescales (cf. Fig. 4.1 and Eq. 4.18). I show that the underlying model establishes a connection between the seemingly unrelated and previously unexplained force peaks in single molecule nucleosome and DNA-toroid stretching experiments. Finally in chapter 5 I show that toroidally confined DNA (found in viruses, DNAcondensates or sperm chromatin) undergoes a transition to a twisted, highly entangled state provided that the aspect ratio of the underlying torus crosses a certain critical value (cf. Eq. 5.6 and the phase diagram in Fig. 5.4). The presented mechanism could rationalize several experimental mysteries, ranging from entangled and supercoiled toroids released from virus capsids to the unexpectedly short cholesteric pitch in the (toroidaly wound) sperm chromatin. I propose that the ”topological encapsulation” resulting from our model may have some practical implications for the gene-therapeutic DNA delivery process.

Computational modelling of human stem cell-derived cardiomyocytes - Texture descriptors for biological image processing

This thesis investigates two distinct research topics. The main topic (Part I) is the computational modelling of cardiomyocytes derived from human stem cells, both embryonic (hESC-CM) and induced-pluripotent (hiPSC-CM). The aim of this research line lies in developing models of the electrophysiology of hESC-CM and hiPSC-CM in order to integrate the available experimental data and getting in-silico models to be used for studying/making new hypotheses/planning experiments on aspects not fully understood yet, such as the maturation process, the functionality of the Ca2+ hangling or why the hESC-CM/hiPSC-CM action potentials (APs) show some differences with respect to APs from adult cardiomyocytes. Chapter I.1 introduces the main concepts about hESC-CMs/hiPSC-CMs, the cardiac AP, and computational modelling. Chapter I.2 presents the hESC-CM AP model, able to simulate the maturation process through two developmental stages, Early and Late, based on experimental and literature data. Chapter I.3 describes the hiPSC-CM AP model, able to simulate the ventricular-like and atrial-like phenotypes. This model was used to assess which currents are responsible for the differences between the ventricular-like AP and the adult ventricular AP. The secondary topic (Part II) consists in the study of texture descriptors for biological image processing. Chapter II.1 provides an overview on important texture descriptors such as Local Binary Pattern or Local Phase Quantization. Moreover the non-binary coding and the multi-threshold approach are here introduced. Chapter II.2 shows that the non-binary coding and the multi-threshold approach improve the classification performance of cellular/sub-cellular part images, taken from six datasets. Chapter II.3 describes the case study of the classification of indirect immunofluorescence images of HEp2 cells, used for the antinuclear antibody clinical test. Finally the general conclusions are reported.

Early steps in ventral nerve cord development in chelicerates and myriapods and formation of brain compartments in spiders

The central point of this work is the investigation of neurogenesis in chelicerates and myriapods. By comparing decisive mechanisms in neurogenesis in the four arthropod groups (Chelicerata, Crustacea, Insecta, Myriapoda) I was able to show which of these mechanisms are conserved and which developmental modules have diverged. Thereby two processes of embryonic development of the central nervous system were brought into focus. On the one hand I studied early neurogenesis in the ventral nerve cord of the spiders Cupiennius salei and Achaearanea tepidariorum and the millipede Glomeris marginata and on the other hand the development of the brain in Cupiennius salei.rnWhile the nervous system of insects and crustaceans is formed by the progeny of single neural stem cells (neuroblasts), in chelicerates and myriapods whole groups of cells adopt the neural cell fate and give rise to the ventral nerve cord after their invagination. The detailed comparison of the positions and the number of the neural precursor groups within the neuromeres in chelicerates and myriapods showed that the pattern is almost identical which suggests that the neural precursors groups in these arthropod groups are homologous. This pattern is also very similar to the neuroblast pattern in insects. This raises the question if the mechanisms that confer regional identity to the neural precursors is conserved in arthropods although the mode of neural precursor formation is different. The analysis of the functions and expression patterns of genes which are known to be involved in this mechanism in Drosophila melanogaster showed that neural patterning is highly conserved in arthropods. But I also discovered differences in early neurogenesis which reflect modifications and adaptations in the development of the nervous systems in the different arthropod groups.rnThe embryonic development of the brain in chelicerates which was investigated for the first time in this work shows similarities but also some modifications to insects. In vertebrates and arthropods the adult brain is composed of distinct centres with different functions. Investigating how these centres, which are organised in smaller compartments, develop during embryogenesis was part of this work. By tracing the morphogenetic movements and analysing marker gene expressions I could show the formation of the visual brain centres from the single-layered precheliceral neuroectoderm. The optic ganglia, the mushroom bodies and the arcuate body (central body) are formed by large invaginations in the peripheral precheliceral neuroectoderm. This epithelium itself contains neural precursor groups which are assigned to the respective centres and thereby build the three-dimensional optical centres. The single neural precursor groups are distinguishable during this process leading to the assumption that they carry positional information which might subdivide the individual brain centres into smaller functional compartments.rn

Relationship between bispectral index values and volatile anesthetic concentrations during the maintenance phase of anesthesia in the B-Unaware trial

Hypnotic depth during anesthesia affects electroencephalography waveforms and electroencephalogram-derived indices, such as the bispectral index (BIS). Titrating anesthetic administration against the BIS assumes reliable relationships between BIS values, electroencephalogram waveforms, and effect site concentration, beyond loss of responsiveness. Associations among BIS, end-tidal anesthetic concentrations (ETAC), and patient characteristics were examined during anesthetic maintenance, using B-Unaware trial data.

Extensive clonal spread and extreme longevity in saw palmetto, a foundation clonal plant

The lack of effective tools have hampered our ability to assess the size, growth and ages of clonal plants. With Serenoa repens (saw palmetto) as a model, we introduce a novel analytical framework that integrates DNA fingerprinting and mathematical modelling to simulate growth and estimate ages of clonal plants. We also demonstrate the application of such life-history information of clonal plants to provide insight into management plans. Serenoa is an ecologically important foundation species in many Southeastern United States ecosystems; yet, many land managers consider Serenoa a troublesome invasive plant. Accordingly, management plans have been developed to reduce or eliminate Serenoa with little understanding of its life history. Using Amplified Fragment Length Polymorphisms, we genotyped 263 Serenoa and 134 Sabal etonia (a sympatric non-clonal palmetto) samples collected from a 20 X 20 m study plot in Florida scrub. Sabal samples were used to assign small field-unidentifiable palmettos to Serenoa or Sabal and also as a negative control for clone detection. We then mathematically modelled clonal networks to estimate genet ages. Our results suggest that Serenoa predominantly propagate via vegetative sprouts and 10000-year-old genets may be common, while showing no evidence of clone formation by Sabal. The results of this and our previous studies suggest that: (i) Serenoa has been part of scrub associations for thousands of years, (ii) Serenoa invasion are unlikely and (ii) once Serenoa is eliminated from local communities, its restoration will be difficult. Reevaluation of the current management tools and plans is an urgent task.

Extensive clonal spread and extreme longevity in saw palmetto, a foundation clonal plant

The lack of effective tools has hampered our ability to assess the size, growth and ages of clonal plants. With Serenoa repens (saw palmetto) as a model, we introduce a novel analytical frame work that integrates DNA fingerprinting and mathematical modelling to simulate growth and estimate ages of clonal plants. We also demonstrate the application of such life-history information of clonal plants to provide insight into management plans. Serenoa is an ecologically important foundation species in many Southeastern United States ecosystems; yet, many land managers consider Serenoa a troublesome invasive plant. Accordingly, management plans have been developed to reduce or eliminate Serenoa with little understanding of its life history. Using Amplified Fragment Length Polymorphisms, we genotyped 263 Serenoa and 134 Sabal etonia (a sympatric non-clonal palmetto) samples collected from a 20 x 20 m study plot in Florida scrub. Sabal samples were used to assign small field-unidentifiable palmettos to Serenoa or Sabal and also as a negative control for clone detection. We then mathematically modelled clonal networks to estimate genet ages. Our results suggest that Serenoa predominantly propagate via vegetative sprouts and 10000-year-old genets maybe common, while showing no evidence of clone formation by Sabal. The results of this and our previous studies suggest that: (i) Serenoa has been part of scrub associations for thousands of years, (ii) Serenoa invasions are unlikely and (ii) once Serenoa is eliminated from local communities, its restoration will be difficult. Reevaluation of the current management tools and plans is an urgent task.

The structural architecture of adult mammalian articular cartilage evolves by a synchronized process of tissue resorption and neoformation during postnatal development

OBJECTIVE: During postnatal development, mammalian articular cartilage acts as a surface growth plate for the underlying epiphyseal bone. Concomitantly, it undergoes a fundamental process of structural reorganization from an immature isotropic to a mature (adult) anisotropic architecture. However, the mechanism underlying this structural transformation is unknown. It could involve either an internal remodelling process, or complete resorption followed by tissue neoformation. The aim of this study was to establish which of these two alternative tissue reorganization mechanisms is physiologically operative. We also wished to pinpoint the articular cartilage source of the stem cells for clonal expansion and the zonal location of the chondrocyte pool with high proliferative activity. METHODS: The New Zealand white rabbit served as our animal model. The analysis was confined to the high-weight-bearing (central) areas of the medial and lateral femoral condyles. After birth, the articular cartilage layer was evaluated morphologically at monthly intervals from the first to the eighth postnatal month, when this species attains skeletal maturity. The overall height of the articular cartilage layer at each juncture was measured. The growth performance of the articular cartilage layer was assessed by calcein labelling, which permitted an estimation of the daily growth rate of the epiphyseal bone and its monthly length-gain. The slowly proliferating stem-cell pool was identified immunohistochemically (after labelling with bromodeoxyuridine), and the rapidly proliferating chondrocyte population by autoradiography (after labelling with (3)H-thymidine). RESULTS: The growth activity of the articular cartilage layer was highest 1 month after birth. It declined precipitously between the first and third months, and ceased between the third and fourth months, when the animal enters puberty. The structural maturation of the articular cartilage layer followed a corresponding temporal trend. During the first 3 months, when the articular cartilage layer is undergoing structural reorganization, the net length-gain in the epiphyseal bone exceeded the height of the articular cartilage layer. This finding indicates that the postnatal reorganization of articular cartilage from an immature isotropic to a mature anisotropic structure is not achieved by a process of internal remodelling, but by the resorption and neoformation of all zones except the most superficial (stem-cell) one. The superficial zone was found to consist of slowly dividing stem cells with bidirectional mitotic activity. In the horizontal direction, this zone furnishes new stem cells that replenish the pool and effect a lateral expansion of the articular cartilage layer. In the vertical direction, the superficial zone supplies the rapidly dividing, transit-amplifying daughter-cell pool that feeds the transitional and upper radial zones during the postnatal growth phase of the articular cartilage layer. CONCLUSIONS: During postnatal development, mammalian articular cartilage fulfils a dual function, viz., it acts not only as an articulating layer but also as a surface growth plate. In the lapine model, this growth activity ceases at puberty (3-4 months of age), whereas that of the true (metaphyseal) growth plate continues until the time of skeletal maturity (8 months). Hence, the two structures are regulated independently. The structural maturation of the articular cartilage layer coincides temporally with the cessation of its growth activity - for the radial expansion and remodelling of the epiphyseal bone - and with sexual maturation. That articular cartilage is physiologically reorganized by a process of tissue resorption and neoformation, rather than by one of internal remodelling, has important implications for the functional engineering and repair of articular cartilage tissue.

Impact of weather and climate on the incidence of acute coronary syndromes

BACKGROUND: First investigations of the interactions between weather and the incidence of acute myocardial infarctions date back to 1938. The early observation of a higher incidence of myocardial infarctions in the cold season could be confirmed in very different geographical regions and cohorts. While the influence of seasonal variations on the incidence of myocardial infarctions has been extensively documented, the impact of individual meteorological parameters on the disease has so far not been investigated systematically. Hence the present study intended to assess the impact of the essential variables of weather and climate on the incidence of myocardial infarctions. METHODS: The daily incidence of myocardial infarctions was calculated from a national hospitalization survey. The hourly weather and climate data were provided by the database of the national weather forecast. The epidemiological and meteorological data were correlated by multivariate analysis based on a generalized linear model assuming a log-link-function and a Poisson distribution. RESULTS: High ambient pressure, high pressure gradients, and heavy wind activity were associated with an increase in the incidence of the totally 6560 hospitalizations for myocardial infarction irrespective of the geographical region. Snow- and rainfall had inconsistent effects. Temperature, Foehn, and lightning showed no statistically significant impact. CONCLUSIONS: Ambient pressure, pressure gradient, and wind activity had a statistical impact on the incidence of myocardial infarctions in Switzerland from 1990 to 1994. To establish a cause-and-effect relationship more data are needed on the interaction between the pathophysiological mechanisms of the acute coronary syndrome and weather and climate variables.

Association between sleep apnea severity and blood coagulability: Treatment effects of nasal continuous positive airway pressure

A prothrombotic state may contribute to the elevated cardiovascular risk in patients with obstructive sleep apnea (OSA). We investigated the relationship between apnea severity and hemostasis factors and effect of continuous positive airway pressure (CPAP) treatment on hemostatic activity. We performed full overnight polysomnography in 44 OSA patients (mean age 47+/-10 years), yielding apnea-hypopnea index (AHI) and mean nighttime oxyhemoglobin saturation (SpO2) as indices of apnea severity. For treatment, subjects were double-blind randomized to 2 weeks of either therapeutic CPAP (n = 18), 3 l/min supplemental nocturnal oxygen (n = 16) or placebo-CPAP (<1 cm H2O) (n = 10). Levels of von Willebrand factor antigen (VWF:Ag), soluble tissue factor (sTF), D-dimer, and plasminogen activator inhibitor (PAI)-1 antigen were measured in plasma pre- and posttreatment. Before treatment, PAI-1 was significantly correlated with AHI (r = 0.47, p = 0.001) and mean nighttime SpO2 (r = -0.32, p = 0.035), but these OSA measures were not significantly related with VWF:Ag, sTF, and D-dimer. AHI was a significant predictor of PAI-1 (R2 = 0.219, standardized beta = 0.47, p = 0.001), independent of mean nighttime SpO2, body mass index (BMI), and age. A weak time-by-treatment interaction for PAI-1 was observed (p = 0.041), even after adjusting for age, BMI, pre-treatment AHI, and mean SpO2 (p = 0.046). Post hoc analyses suggested that only CPAP treatment was associated with a decrease in PAI-1 (p = 0.039); there were no changes in VWF:Ag, sTF, and D-dimer associated with treatment with placebo-CPAP or with nocturnal oxygen. Apnea severity may be associated with impairment in the fibrinolytic capacity. To the extent that our sample size was limited, the observation that CPAP treatment led to a decrease in PAI-1 in OSA must be regarded as tentative.

A unifying approach for haplotype analysis of quantitative traits in family-based association studies: Testing and estimating gene-environment interactions with complex exposure variables

We propose robust and e±cient tests and estimators for gene-environment/gene-drug interactions in family-based association studies. The methodology is designed for studies in which haplotypes, quantitative pheno- types and complex exposure/treatment variables are analyzed. Using causal inference methodology, we derive family-based association tests and estimators for the genetic main effects and the interactions. The tests and estimators are robust against population admixture and strati¯cation without requiring adjustment for confounding variables. We illustrate the practical relevance of our approach by an application to a COPD study. The data analysis suggests a gene-environment interaction between a SNP in the Serpine gene and smok- ing status/pack years of smoking that reduces the FEV1 volume by about 0.02 liter per pack year of smoking. Simulation studies show that the pro- posed methodology is su±ciently powered for realistic sample sizes and that it provides valid tests and effect size estimators in the presence of admixture and stratification.

Invasion patterns of emerald ash borer and European earthworms in forested ecosystems

Invasive and exotic species present a serious threat to the health and sustainability of natural ecosystems. These species often benefit from anthropogenic activities that aid their introduction and dispersal. This dissertation focuses on invasion dynamics of the emerald ash borer, native to Asia, and European earthworms. These species have shown detrimental impacts in invaded forest ecosystems across the Great Lakes region, and continue to spread via human-assisted long distance dispersal and by natural modes of dispersal into interior forests from areas of introduction. Successful forest management requires that the impact and effect of invasive species be considered and incorporated into management plans. Understanding patterns and constraints of introduction, establishment, and spread will aid in this effort. To assist in efforts to locate introduction points of emerald ash borer, a multicriteria risk model was developed to predict the highest risk areas. Important parameters in the model were road proximity, land cover type, and campground proximity. The model correctly predicted 85% of known emerald ash borer invasion sites to be at high risk. The model’s predictions across northern Michigan can be used to focus and guide future monitoring efforts. Similar modeling efforts were applied to the prediction of European earthworm invasion in northern Michigan forests. Field sampling provided a means to improve upon modeling efforts for earthworms to create current and future predictions of earthworm invasion. Those sites with high soil pH and high basal area of earthworm preferred overstory species (such as basswood and maples) had the highest likelihood of European earthworm invasion. Expanding beyond Michigan into the Upper Great Lakes region, earthworm populations were sampled across six National Wildlife Refuges to identify potential correlates and deduce specific drivers and constraints of earthworm invasion. Earthworm communities across all refuges were influenced by patterns of anthropogenic activity both within refuges and in surrounding ecoregions of study. Forest composition, soil pH, soil organic matter, anthropogenic cover, and agriculture proximity also proved to be important drivers of earthworm abundance and community composition. While there are few management options to remove either emerald ash borer or European earthworms from forests after they have become well established, prevention and early detection are important and can be beneficial. An improved understanding the factors controlling the distribution and invasion patterns of exotic species across the landscape will aid efforts to determine their consequences and generate appropriate forest management solutions to sustain ecosystem health in the presence of these invaders.

Critical evaluation of diameter increment modelling in the Great Lakes region

Simulations of forest stand dynamics in a modelling framework including Forest Vegetation Simulator (FVS) are diameter driven, thus the diameter or basal area increment model needs a special attention. This dissertation critically evaluates diameter or basal area increment models and modelling approaches in the context of the Great Lakes region of the United States and Canada. A set of related studies are presented that critically evaluate the sub-model for change in individual tree basal diameter used in the Forest Vegetation Simulator (FVS), a dominant forestry model in the Great Lakes region. Various historical implementations of the STEMS (Stand and Tree Evaluation and Modeling System) family of diameter increment models, including the current public release of the Lake States variant of FVS (LS-FVS), were tested for the 30 most common tree species using data from the Michigan Forest Inventory and Analysis (FIA) program. The results showed that current public release of the LS-FVS diameter increment model over-predicts 10-year diameter increment by 17% on average. Also the study affirms that a simple adjustment factor as a function of a single predictor, dbh (diameter at breast height) used in the past versions, provides an inadequate correction of model prediction bias. In order to re-engineer the basal diameter increment model, the historical, conceptual and philosophical differences among the individual tree increment model families and their modelling approaches were analyzed and discussed. Two underlying conceptual approaches toward diameter or basal area increment modelling have been often used: the potential-modifier (POTMOD) and composite (COMP) approaches, which are exemplified by the STEMS/TWIGS and Prognosis models, respectively. It is argued that both approaches essentially use a similar base function and neither is conceptually different from a biological perspective, even though they look different in their model forms. No matter what modelling approach is used, the base function is the foundation of an increment model. Two base functions – gamma and Box-Lucas – were identified as candidate base functions for forestry applications. The results of a comparative analysis of empirical fits showed that quality of fit is essentially similar, and both are sufficiently detailed and flexible for forestry applications. The choice of either base function in order to model diameter or basal area increment is dependent upon personal preference; however, the gamma base function may be preferred over the Box-Lucas, as it fits the periodic increment data in both a linear and nonlinear composite model form. Finally, the utility of site index as a predictor variable has been criticized, as it has been widely used in models for complex, mixed species forest stands though not well suited for this purpose. An alternative to site index in an increment model was explored, using site index and a combination of climate variables and Forest Ecosystem Classification (FEC) ecosites and data from the Province of Ontario, Canada. The results showed that a combination of climate and FEC ecosites variables can replace site index in the diameter increment model.