Modeling the dynamics of human hair cycles by a follicular automaton

The hair follicle cycle successively goes through the anagen, catagen, telogen, and latency phases, which correspond, respectively, to hair growth, arrest, shedding, and absence before a new anagen phase is initiated. Experimental observations collected over a period of 14 years in a group of 10 male volunteers, alopecic and nonalopecic, allowed us to determine the characteristics of scalp hair follicle cycles. On the basis of these observations, we propose a follicular automaton model to simulate the dynamics of human hair cycles. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states anagen, telogen, and latency, and the subsequent return to anagen. The transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a variance. The follicular automaton model accounts both for the dynamical transitions observed in a single follicle and for the behavior of an ensemble of independently cycling follicles. Thus, the model successfully reproduces the evolution of the fractions of follicle populations in each of the three phases, which fluctuate around steady-state or slowly drifting values. We apply the follicular automaton model to the study of spatial patterns of follicular growth that result from a spatially heterogeneous distribution of parameters such as the mean duration of anagen phase. When considering that follicles die or miniaturize after going through a critical number of successive cycles, the model can reproduce the evolution to hair patterns similar to well known types of diffuse or androgenetic alopecia.

Assessment and discrimination of odor stimuli in rat olfactory bulb by dynamic functional MRI

Dynamic blood oxygenation level-dependent functional MRI was applied at 7 T in the rat olfactory bulb (OB) with pulsed delivery of iso-amyl acetate (IAA) and limonene. Acquisition times for single-slice and whole OB data were 8 and 32 s, respectively, with spatial resolution of 220 × 220 × 250 μm. On an intrasubject basis, short IAA exposures of 0.6 min separated by 3.5-min intervals induced reproducible spatial activity patterns (SAPs) in the olfactory nerve layer, glomerular layer, and external plexiform layer. During long exposures (≈10 min), the initially dominant dorsal SAPs declined in intensity and area, whereas in some OB regions, the initially weak ventral/lateral SAPs increased first and then decreased. The SAPs of different concentrations were topologically similar, which implies that whereas an odor at various concentrations activates the same subsets of receptor cells, different concentrations are assessed and discriminated by variable magnitudes of laminarspecific activations. IAA and limonene reproducibly activated different subsets of receptor cells with some overlaps. Whereas qualitative topographical agreement was observed with results from other methods, the current dynamic blood oxygenation level-dependent functional MRI results can provide quantitative SAPs of the entire OB.

Dynamic localization of a cytoplasmic signal transduction response regulator controls morphogenesis during the Caulobacter cell cycle

We present evidence that a bacterial signal transduction cascade that couples morphogenesis with cell cycle progression is regulated by dynamic localization of its components. Previous studies have implicated two histidine kinases, DivJ and PleC, and the response regulator, DivK, in the regulation of morphogenesis in the dimorphic bacterium Caulobacter crescentus. Here, we show that the cytoplasmic response regulator, DivK, exhibits a dynamic, cyclical localization that culminates in asymmetric distribution of DivK within the two cell types that are characteristic of the Caulobacter cell cycle; DivK is dispersed throughout the cytoplasm of the progeny swarmer cell and is localized to the pole of the stalked cell. The membrane-bound DivJ and PleC histidine kinases, which are asymmetrically localized at the opposite poles of the predivisional cell, control the temporal and spatial localization of DivK. DivJ mediates DivK targeting to the poles whereas PleC controls its release from one of the poles at times and places that are consistent with the activities and location of DivJ and PleC in the late predivisional cell. Thus, dynamic changes in subcellular location of multiple components of a signal transduction cascade may constitute a novel mode of prokaryotic regulation to generate and maintain cellular asymmetry.

The neural development and organization of letter recognition: Evidence from functional neuroimaging, computational modeling, and behavioral studies

Although much of the brain’s functional organization is genetically predetermined, it appears that some noninnate functions can come to depend on dedicated and segregated neural tissue. In this paper, we describe a series of experiments that have investigated the neural development and organization of one such noninnate function: letter recognition. Functional neuroimaging demonstrates that letter and digit recognition depend on different neural substrates in some literate adults. How could the processing of two stimulus categories that are distinguished solely by cultural conventions become segregated in the brain? One possibility is that correlation-based learning in the brain leads to a spatial organization in cortex that reflects the temporal and spatial clustering of letters with letters in the environment. Simulations confirm that environmental co-occurrence does indeed lead to spatial localization in a neural network that uses correlation-based learning. Furthermore, behavioral studies confirm one critical prediction of this co-occurrence hypothesis, namely, that subjects exposed to a visual environment in which letters and digits occur together rather than separately (postal workers who process letters and digits together in Canadian postal codes) do indeed show less behavioral evidence for segregated letter and digit processing.

Dynamic contrast-enhanced magnetic resonance imaging reveals stress-induced angiogenesis in MCF7 human breast tumors.

The mechanism of contrast enhancement of tumors using magnetic resonance imaging was investigated in MCF7 human breast cancer implanted in nude mice. Dynamic contrast-enhanced images recorded at high spatial resolution were analyzed by an image analysis method based on a physiological model, which included the blood circulation, the tumor, the remaining tissues, and clearance via the kidneys. This analysis enabled us to map in rapidly enhancing regions within the tumor, the capillary permeability factor (capillary permeability times surface area per voxel volume) and the fraction of leakage space. Correlation of these maps with T2-weighted spin echo images, with histopathology, and with immunohistochemical staining of endothelial cells demonstrated the presence of dense permeable microcapillaries in the tumor periphery and in intratumoral regions that surrounded necrotic loci. The high leakage from the intratumoral permeable capillaries indicated an induction of a specific angiogenic process associated with stress conditions that cause necrosis. This induction was augmented in tumors responding to tamoxifen treatment. Determination of the distribution and extent of this stress-induced angiogenic activity by contrast-enhanced MRI might be of diagnostic and of prognostic value.

Spatial integration and cortical dynamics.

Cells in adult primary visual cortex are capable of integrating information over much larger portions of the visual field than was originally thought. Moreover, their receptive field properties can be altered by the context within which local features are presented and by changes in visual experience. The substrate for both spatial integration and cortical plasticity is likely to be found in a plexus of long-range horizontal connections, formed by cortical pyramidal cells, which link cells within each cortical area over distances of 6-8 mm. The relationship between horizontal connections and cortical functional architecture suggests a role in visual segmentation and spatial integration. The distribution of lateral interactions within striate cortex was visualized with optical recording, and their functional consequences were explored by using comparable stimuli in human psychophysical experiments and in recordings from alert monkeys. They may represent the substrate for perceptual phenomena such as illusory contours, surface fill-in, and contour saliency. The dynamic nature of receptive field properties and cortical architecture has been seen over time scales ranging from seconds to months. One can induce a remapping of the topography of visual cortex by making focal binocular retinal lesions. Shorter-term plasticity of cortical receptive fields was observed following brief periods of visual stimulation. The mechanisms involved entailed, for the short-term changes, altering the effectiveness of existing cortical connections, and for the long-term changes, sprouting of axon collaterals and synaptogenesis. The mutability of cortical function implies a continual process of calibration and normalization of the perception of visual attributes that is dependent on sensory experience throughout adulthood and might further represent the mechanism of perceptual learning.

Modeling and optimization of populations subject to time-dependent mutation.

It has become clear that many organisms possess the ability to regulate their mutation rate in response to environmental conditions. So the question of finding an optimal mutation rate must be replaced by that of finding an optimal mutation schedule. We show that this task cannot be accomplished with standard population-dynamic models. We then develop a "hybrid" model for populations experiencing time-dependent mutation that treats population growth as deterministic but the time of first appearance of new variants as stochastic. We show that the hybrid model agrees well with a Monte Carlo simulation. From this model, we derive a deterministic approximation, a "threshold" model, that is similar to standard population dynamic models but differs in the initial rate of generation of new mutants. We use these techniques to model antibody affinity maturation by somatic hypermutation. We had previously shown that the optimal mutation schedule for the deterministic threshold model is phasic, with periods of mutation between intervals of mutation-free growth. To establish the validity of this schedule, we now show that the phasic schedule that optimizes the deterministic threshold model significantly improves upon the best constant-rate schedule for the hybrid and Monte Carlo models.

Comparison of the defensive elicitation activity of cerato-populin and cerato-platanin: a structural model

Plants can defend themselves from potential pathogenic microorganisms relying on a complex interplay of signaling pathways: activation of the MAPK cascade, transcription of defense related genes, production of reactive oxygen species, nitric oxide and synthesis of other defensive compounds such as phytoalexins. These events are triggered by the recognition of pathogen’s effectors (effector-triggered immunity) or PAMPs (PAMP-triggered immunity). The Cerato Platanin Family (CPF) members are Cys-rich proteins secreted and localized on fungal cell walls, involved in several aspects of fungal development and pathogen-host interactions. Although more than hundred genes of the CPF have been identified and analyzed, the structural and functional characterization of the expressed proteins has been restricted only to few members of the family. Interestingly, those proteins have been shown to bind chitin with diverse affinity and after foliar treatment they elicit defensive mechanisms in host and non-host plants. This property turns cerato platanins into interesting candidates, worth to be studied to develop new fungal elicitors with applications in sustainable agriculture. This study focus on cerato-platanin (CP), core member of the family and on the orthologous cerato-populin (Pop1). The latter shows an identity of 62% and an overall homology of 73% with respect to CP. Both proteins are able to induce MAPKs phosphorylation, production of reactive oxygen species and nitric oxide, overexpression of defense’s related genes, programmed cell death and synthesis of phytoalexins. CP, however, when compared to Pop1, induces a faster response and, in some cases, a stronger activity on plane leaves. Aim of the present research is to verify if the dissimilarities observed in the defense elicitation activity of these proteins can be associated to their structural and dynamic features. Taking advantage of the available CP NMR structure, Pop1’s 3D one was obtained by homology modeling. Experimental residual dipolar couplings and 1H, 15N, 13C resonance assignments were used to validate the model. Previous works on CPF members, addressed the highly conserved random coil regions (loops b1-b2 and b2-b3) as sufficient and necessary to induce necrosis in plants’ leaves: that region was investigated in both Pop1 and CP. In the two proteins the loops differ, in their primary sequence, for few mutations and an insertion with a consequent diversification of the proteins’ electrostatic surface. A set of 2D and 3D NMR experiments was performed to characterize both the spatial arrangement and the dynamic features of the loops. NOE data revealed a more extended network of interactions between the loops in Pop1 than in CP. In addition, in Pop1 we identified a salt bridge Lys25/Asp52 and a strong hydrophobic interaction between Phe26/Trp53. These structural features were expected not only to affect the loops’ spatial arrangement, but also to reduce the degree of their conformational freedom. Relaxation data and the order parameter S2 indeed highlighted reduced flexibility, in particular for loop b1-b2 of Pop1. In vitro NMR experiments, where Pop1 and CP were titrated with oligosaccharides, supported the hypothesis that the loops structural and dynamic differences may be responsible for the different chitin-binding properties of the two proteins: CP selectively binds tetramers of chitin in a shallow groove on one side of the barrel defined by loops b1-b2, b2-b3 and b4-b5, Pop1, instead, interacts in a non-specific fashion with oligosaccharides. Because the region involved in chitin-binding is also responsible for the defense elicitation activity, possibly being recognized by plant's receptors, it is reasonable to expect that those structural and dynamic modifications may also justify the different extent of defense elicitation. To test that hypothesis, the initial steps of a protocol aimed to the identify a receptor for CP, in silico, are presented.

Logic-Based Outer Approximation for the Design of Discrete-Continuous Dynamic Systems with Implicit Discontinuities

We address the optimization of discrete-continuous dynamic optimization problems using a disjunctive multistage modeling framework, with implicit discontinuities, which increases the problem complexity since the number of continuous phases and discrete events is not known a-priori. After setting a fixed alternative sequence of modes, we convert the infinite-dimensional continuous mixed-logic dynamic (MLDO) problem into a finite dimensional discretized GDP problem by orthogonal collocation on finite elements. We use the Logic-based Outer Approximation algorithm to fully exploit the structure of the GDP representation of the problem. This modelling framework is illustrated with an optimization problem with implicit discontinuities (diver problem).

Amphibian Hotspots and Conservation Priorities in Eastern Cuba Identified by Species Distribution Modeling

The high rate of amphibian endemism and the severe habitat modification in the Caribbean islands make them an ideal place to test if the current protected areas network might protect this group. In this study, we model distribution and map species richness of the 40 amphibian species from eastern Cuba with the objectives of identify hotspots, detect gaps in species representation in protected areas, and select additional areas to fill these gaps. We used two modeling methods, Maxent and Habitat Suitability Models, to reach a consensus distribution map for each species, then calculate species richness by combining specific models and finally performed gap analyses for species and hotspots. Our results showed that the models were robust enough to predict species distributions and that most of the amphibian hotspots were represented in reserves, but 50 percent of the species were incompletely covered and Eleutherodactylus rivularis was totally uncovered by the protected areas. We identified 1441 additional km2 (9.9% of the study area) that could be added to the current protected areas, allowing the representation of every species and all hotspots. Our results are relevant for the conservation planning in other Caribbean islands, since studies like this could contribute to fill the gaps in the existing protected areas and to design a future network. Both cases would benefit from modeling amphibian species distribution using available data, even if they are incomplete, rather than relying only in the protection of known or suspected hotspots.

Modélisation de l’activité électrique des oreillettes avant et après ablation par cathéter

- Réalisé au centre de recherche de l'hospital du Sacré-Coeur de Montréal. - Programme conjoint entre Université de Montréal et École Polytechnique de Montréal.

Modeling and validation of the thermal behavior of buildings for the development of demand response methods

The representation of the thermal behaviour of the building is achieved through a relatively simple dynamic model that takes into account the effects due to the thermal mass of the building components. The model of a intra-floor apartment has been built in the Matlab-Simulink environment and considers the heat transmission through the external envelope, wall and windows, the internal thermal masses, (i.e. furniture, internal wall and floor slabs) and the sun gain due to opaque and see-through surfaces of the external envelope. The simulations results for the entire year have been compared and the model validated, with the one obtained with the dynamic building simulation software Energyplus.

A 20 m spatial resolution seamless multisource Digital Elevation/Depth Model for Lizard Island, northern Great Barrier Reef

Hydrographers have traditionally referred to the nearshore area as the "white ribbon" area due to the challenges associated with the collection of elevation data in this highly dynamic transitional zone between terrestrial and marine environments. Accordingly, available information in this zone is typically characterised by a range of datasets from disparate sources. In this paper we propose a framework to 'fill' the white ribbon area of a coral reef system by integrating multiple elevation and bathymetric datasets acquired by a suite of remote-sensing technologies into a seamless digital elevation model (DEM). A range of datasets are integrated, including field-collected GPS elevation points, terrestrial and bathymetric LiDAR, single and multibeam bathymetry, nautical chart depths and empirically derived bathymetry estimations from optical remote sensing imagery. The proposed framework ranks data reliability internally, thereby avoiding the requirements to quantify absolute error and results in a high resolution, seamless product. Nested within this approach is an effective spatially explicit technique for improving the accuracy of bathymetry estimates derived empirically from optical satellite imagery through modelling the spatial structure of residuals. The approach was applied to data collected on and around Lizard Island in northern Australia. Collectively, the framework holds promise for filling the white ribbon zone in coastal areas characterised by similar data availability scenarios. The seamless DEM is referenced to the horizontal coordinate system MGA Zone 55 - GDA 1994, mean sea level (MSL) vertical datum and has a spatial resolution of 20 m.

The effects of age, spatial ability, and navigational information on navigation performance.

Federal Highway Administration, Washington, D.C.

Transformation of Nine-Accelerometer-Package (NAP) data for replicating headpart kinematics and dynamic loading. Final report.

National Highway Traffic Safety Administration, Washington, D.C.