Extension of Murray's law using a non-Newtonian model of blood flow.

BACKGROUND: So far, none of the existing methods on Murray's law deal with the non-Newtonian behavior of blood flow although the non-Newtonian approach for blood flow modelling looks more accurate. MODELING: In the present paper, Murray's law which is applicable to an arterial bifurcation, is generalized to a non-Newtonian blood flow model (power-law model). When the vessel size reaches the capillary limitation, blood can be modeled using a non-Newtonian constitutive equation. It is assumed two different constraints in addition to the pumping power: the volume constraint or the surface constraint (related to the internal surface of the vessel). For a seek of generality, the relationships are given for an arbitrary number of daughter vessels. It is shown that for a cost function including the volume constraint, classical Murray's law remains valid (i.e. SigmaR(c) = cste with c = 3 is verified and is independent of n, the dimensionless index in the viscosity equation; R being the radius of the vessel). On the contrary, for a cost function including the surface constraint, different values of c may be calculated depending on the value of n. RESULTS: We find that c varies for blood from 2.42 to 3 depending on the constraint and the fluid properties. For the Newtonian model, the surface constraint leads to c = 2.5. The cost function (based on the surface constraint) can be related to entropy generation, by dividing it by the temperature. CONCLUSION: It is demonstrated that the entropy generated in all the daughter vessels is greater than the entropy generated in the parent vessel. Furthermore, it is shown that the difference of entropy generation between the parent and daughter vessels is smaller for a non-Newtonian fluid than for a Newtonian fluid.

In vitro characterization of PlySK1249, a novel phage lysin, and assessment of its antibacterial activity in a mouse model of Streptococcus agalactiae bacteremia.

Beta-hemolytic Streptococcus agalactiae is the leading cause of bacteremia and invasive infections. These diseases are treated with β-lactams or macrolides, but the emergence of less susceptible and even fully resistant strains is a cause for concern. New bacteriophage lysins could be promising alternatives against such organisms. They hydrolyze the bacterial peptidoglycan at the end of the phage cycle, in order to release the phage progeny. By using a bioinformatic approach to screen several beta-hemolytic streptococci, a gene coding for a lysin was identified on a prophage carried by Streptococcus dysgalactiae subsp. equisimilis SK1249. The gene product, named PlySK1249, harbored an original three-domain structure with a central cell wall-binding domain surrounded by an N-terminal amidase and a C-terminal CHAP domain. Purified PlySK1249 was highly lytic and bactericidal for S. dysgalactiae (2-log10 CFU/ml decrease within 15 min). Moreover, it also efficiently killed S. agalactiae (1.5-log10 CFU/ml decrease within 15 min) but not several streptococcal commensal species. We further investigated the activity of PlySK1249 in a mouse model of S. agalactiae bacteremia. Eighty percent of the animals (n = 10) challenged intraperitoneally with 10(6) CFU of S. agalactiae died within 72 h, whereas repeated injections of PlySK1249 (45 mg/kg 3 times within 24 h) significantly protected the mice (P < 0.01). Thus, PlySK1249, which was isolated from S. dysgalactiae, demonstrated high cross-lytic activity against S. agalactiae both in vitro and in vivo. These encouraging results indicated that PlySK1249 might represent a good candidate to be developed as a new enzybiotic for the treatment of systemic S. agalactiae infections.

Simulation of stromatolite growth using diffusion-limited aggregation and a discrete model of biogeochemical exchanges in microbial mat

Resume : Mieux comprendre les stromatolithes et les tapis microbiens est un sujet important en biogéosciences puisque cela aide à l'étude des premières formes de vie sur Terre, a mieux cerner l'écologie des communautés microbiennes et la contribution des microorganismes a la biominéralisation, et même à poser certains fondements dans les recherches en exobiologie. D'autre part, la modélisation est un outil puissant utilisé dans les sciences naturelles pour appréhender différents phénomènes de façon théorique. Les modèles sont généralement construits sur un système d'équations différentielles et les résultats sont obtenus en résolvant ce système. Les logiciels disponibles pour implémenter les modèles incluent les logiciels mathématiques et les logiciels généraux de simulation. L'objectif principal de cette thèse est de développer des modèles et des logiciels pour aider a comprendre, via la simulation, le fonctionnement des stromatolithes et des tapis microbiens. Ces logiciels ont été développés en C++ en ne partant d'aucun pré-requis de façon a privilégier performance et flexibilité maximales. Cette démarche permet de construire des modèles bien plus spécifiques et plus appropriés aux phénomènes a modéliser. Premièrement, nous avons étudié la croissance et la morphologie des stromatolithes. Nous avons construit un modèle tridimensionnel fondé sur l'agrégation par diffusion limitée. Le modèle a été implémenté en deux applications C++: un moteur de simulation capable d'exécuter un batch de simulations et de produire des fichiers de résultats, et un outil de visualisation qui permet d'analyser les résultats en trois dimensions. Après avoir vérifié que ce modèle peut en effet reproduire la croissance et la morphologie de plusieurs types de stromatolithes, nous avons introduit un processus de sédimentation comme facteur externe. Ceci nous a mené a des résultats intéressants, et permis de soutenir l'hypothèse que la morphologie des stromatolithes pourrait être le résultat de facteurs externes autant que de facteurs internes. Ceci est important car la classification des stromatolithes est généralement fondée sur leur morphologie, imposant que la forme d'un stromatolithe est dépendante de facteurs internes uniquement (c'est-à-dire les tapis microbiens). Les résultats avancés dans ce mémoire contredisent donc ces assertions communément admises. Ensuite, nous avons décidé de mener des recherches plus en profondeur sur les aspects fonctionnels des tapis microbiens. Nous avons construit un modèle bidimensionnel de réaction-diffusion fondé sur la simulation discrète. Ce modèle a été implémenté dans une application C++ qui permet de paramétrer et exécuter des simulations. Nous avons ensuite pu comparer les résultats de simulation avec des données du monde réel et vérifier que le modèle peut en effet imiter le comportement de certains tapis microbiens. Ainsi, nous avons pu émettre et vérifier des hypothèses sur le fonctionnement de certains tapis microbiens pour nous aider à mieux en comprendre certains aspects, comme la dynamique des éléments, en particulier le soufre et l'oxygène. En conclusion, ce travail a abouti à l'écriture de logiciels dédiés à la simulation de tapis microbiens d'un point de vue tant morphologique que fonctionnel, suivant deux approches différentes, l'une holistique, l'autre plus analytique. Ces logiciels sont gratuits et diffusés sous licence GPL (General Public License). Abstract : Better understanding of stromatolites and microbial mats is an important topic in biogeosciences as it helps studying the early forms of life on Earth, provides clues re- garding the ecology of microbial ecosystems and their contribution to biomineralization, and gives basis to a new science, exobiology. On the other hand, modelling is a powerful tool used in natural sciences for the theoretical approach of various phenomena. Models are usually built on a system of differential equations and results are obtained by solving that system. Available software to implement models includes mathematical solvers and general simulation software. The main objective of this thesis is to develop models and software able to help to understand the functioning of stromatolites and microbial mats. Software was developed in C++ from scratch for maximum performance and flexibility. This allows to build models much more specific to a phenomenon rather than general software. First, we studied stromatolite growth and morphology. We built a three-dimensional model based on diffusion-limited aggregation. The model was implemented in two C++ applications: a simulator engine, which can run a batch of simulations and produce result files, and a Visualization tool, which allows results to be analysed in three dimensions. After verifying that our model can indeed reproduce the growth and morphology of several types of stromatolites, we introduced a sedimentation process as an external factor. This lead to interesting results, and allowed to emit the hypothesis that stromatolite morphology may be the result of external factors as much as internal factors. This is important as stromatolite classification is usually based on their morphology, imposing that a stromatolite shape is dependant on internal factors only (i.e. the microbial mat). This statement is contradicted by our findings, Second, we decided to investigate deeper the functioning of microbial mats, We built a two-dimensional reaction-diffusion model based on discrete simulation, The model was implemented in a C++ application that allows setting and running simulations. We could then compare simulation results with real world data and verify that our model can indeed mimic the behaviour of some microbial mats. Thus, we have proposed and verified hypotheses regarding microbial mats functioning in order to help to better understand them, e.g. the cycle of some elements such as oxygen or sulfur. ln conclusion, this PhD provides a simulation software, dealing with two different approaches. This software is free and available under a GPL licence.

Limited role of the c-Jun N-terminal kinase pathway in a neonatal rat model of cerebral hypoxia-ischemia

D-JNKI1, a cell-permeable peptide inhibitor of the c-Jun N-terminal kinase (JNK) pathway, has been shown to be a powerful neuroprotective agent after focal cerebral ischemia in adult mice and young rats. We have investigated the potential neuroprotective effect of D-JNKI1 and the involvement of the JNK pathway in a neonatal rat model of cerebral hypoxia-ischemia. Seven-day-old rats underwent a permanent ligation of the right common carotid artery followed by 2h of hypoxia (8% oxygen). Treatment with D-JNKI1 (0.3mg/kg intraperitoneally) significantly reduced early calpain activation, late caspase-3 activation and, in the thalamus, autophagosome formation, indicating an involvement of JNK in different types of cell death: necrotic, apoptotic and autophagic. However the size of the lesion was unchanged. Further analysis showed that neonatal hypoxia-ischemia induced an immediate decrease in JNK phosphorylation (reflecting mainly P-JNK1) followed by a slow progressive increase (including P-JNK3 54kDa), whereas c-jun and c-fos expression were both strongly activated immediately after hypoxia-ischemia. In conclusion, unlike in adult ischemic models, JNK is only moderately activated after severe cerebral hypoxia-ischemia in neonatal rats and the observed positive effects of D-JNKI1 are insufficient to give neuroprotection. Thus, for perinatal asphyxia, D-JNKI1 can only be considered in association with other therapies.

Respiratory failure in a mouse model of myotonic dystrophy does not correlate with the CTG repeat length.

Myotonic dystrophy (DM1) is a multisystemic disease caused by an expansion of CTG repeats in the region of DMPK, the gene encoding DM protein kinase. The severity of muscle disability in DM1 correlates with the size of CTG expansion. As respiratory failure is one of the main causes of death in DM1, we investigated the correlation between respiratory impairment and size of the (CTG)n repeat in DM1 animal models. Using pressure plethysmography the respiratory function was assessed in control and transgenic mice carrying either 600 (DM600) or >1300 CTG repeats (DMSXL). The statistical analysis of respiratory parameters revealed that both DM1 transgenic mice sub-lines show respiratory impairment compared to control mice. In addition, there is no significant difference in breathing functions between the DM600 and DMSXL mice. In conclusion, these results indicate that respiratory impairment is present in both transgenic mice sub-lines, but the severity of respiratory failure is not related to the size of the (CTG)n expansion.

Continuous positive airway pressure causes lung injury in a model of sepsis.

Continuous positive airway pressure, aimed at preventing pulmonary atelectasis, has been used for decades to reduce lung injury in critically ill patients. In neonatal practice, it is increasingly used worldwide as a primary form of respiratory support due to its low cost and because it reduces the need for endotracheal intubation and conventional mechanical ventilation. We studied the anesthetized in vivo rat and determined the optimal circuit design for delivery of continuous positive airway pressure. We investigated the effects of continuous positive airway pressure following lipopolysaccharide administration in the anesthetized rat. Whereas neither continuous positive airway pressure nor lipopolysaccharide alone caused lung injury, continuous positive airway pressure applied following intravenous lipopolysaccharide resulted in increased microvascular permeability, elevated cytokine protein and mRNA production, and impaired static compliance. A dose-response relationship was demonstrated whereby higher levels of continuous positive airway pressure (up to 6 cmH(2)O) caused greater lung injury. Lung injury was attenuated by pretreatment with dexamethasone. These data demonstrate that despite optimal circuit design, continuous positive airway pressure causes significant lung injury (proportional to the airway pressure) in the setting of circulating lipopolysaccharide. Although we would currently avoid direct extrapolation of these findings to clinical practice, we believe that in the context of increasing clinical use, these data are grounds for concern and warrant further investigation.

The anticancer drug tamoxifen counteracts the pathology in a mouse model of duchenne muscular dystrophy.

Duchenne muscular dystrophy (DMD) is a severe disorder characterized by progressive muscle wasting,respiratory and cardiac impairments, and premature death. No treatment exists so far, and the identification of active substances to fight DMD is urgently needed. We found that tamoxifen, a drug used to treat estrogen-dependent breast cancer, caused remarkable improvements of muscle force and of diaphragm and cardiac structure in the mdx(5Cv) mouse model of DMD. Oral tamoxifen treatment from 3 weeks of age for 15 months at a dose of 10 mg/kg/day stabilized myofiber membranes, normalized whole body force, and increased force production and resistance to repeated contractions of the triceps muscle above normal values. Tamoxifen improved the structure of leg muscles and diminished cardiac fibrosis by~ 50%. Tamoxifen also reduced fibrosis in the diaphragm, while increasing its thickness,myofiber count, and myofiber diameter, thereby augmenting by 72% the amount of contractile tissue available for respiratory function. Tamoxifen conferred a markedly slower phenotype to the muscles.Tamoxifen and its metabolites were present in nanomolar concentrations in plasma and muscles,suggesting signaling through high-affinity targets. Interestingly, the estrogen receptors ERa and ERb were several times more abundant in dystrophic than in normal muscles, and tamoxifen normalized the relative abundance of ERb isoforms. Our findings suggest that tamoxifen might be a useful therapy for DMD.

Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Torque vibration model of axial-flux surface-mounted permanent magnet synchronous machine

In order that the radius and thus ununiform structure of the teeth and otherelectrical and magnetic parts of the machine may be taken into consideration the calculation of an axial flux permanent magnet machine is, conventionally, doneby means of 3D FEM-methods. This calculation procedure, however, requires a lotof time and computer recourses. This study proves that also analytical methods can be applied to perform the calculation successfully. The procedure of the analytical calculation can be summarized into following steps: first the magnet is divided into slices, which makes the calculation for each section individually, and then the parts are submitted to calculation of the final results. It is obvious that using this method can save a lot of designing and calculating time. Thecalculation program is designed to model the magnetic and electrical circuits of surface mounted axial flux permanent magnet synchronous machines in such a way, that it takes into account possible magnetic saturation of the iron parts. Theresult of the calculation is the torque of the motor including the vibrations. The motor geometry and the materials and either the torque or pole angle are defined and the motor can be fed with an arbitrary shape and amplitude of three-phase currents. There are no limits for the size and number of the pole pairs nor for many other factors. The calculation steps and the number of different sections of the magnet are selectable, but the calculation time is strongly depending on this. The results are compared to the measurements of real prototypes. The permanent magnet creates part of the flux in the magnetic circuit. The form and amplitude of the flux density in the air-gap depends on the geometry and material of the magnetic circuit, on the length of the air-gap and remanence flux density of the magnet. Slotting is taken into account by using the Carter factor in the slot opening area. The calculation is simple and fast if the shape of the magnetis a square and has no skew in relation to the stator slots. With a more complicated magnet shape the calculation has to be done in several sections. It is clear that according to the increasing number of sections also the result will become more accurate. In a radial flux motor all sections of the magnets create force with a same radius. In the case of an axial flux motor, each radial section creates force with a different radius and the torque is the sum of these. The magnetic circuit of the motor, consisting of the stator iron, rotor iron, air-gap, magnet and the slot, is modelled with a reluctance net, which considers the saturation of the iron. This means, that several iterations, in which the permeability is updated, has to be done in order to get final results. The motor torque is calculated using the instantaneous linkage flux and stator currents. Flux linkage is called the part of the flux that is created by the permanent magnets and the stator currents passing through the coils in stator teeth. The angle between this flux and the phase currents define the torque created by the magnetic circuit. Due to the winding structure of the stator and in order to limit the leakage flux the slot openings of the stator are normally not made of ferromagnetic material even though, in some cases, semimagnetic slot wedges are used. In the slot opening faces the flux enters the iron almost normally (tangentially with respect to the rotor flux) creating tangential forces in the rotor. This phenomenon iscalled cogging. The flux in the slot opening area on the different sides of theopening and in the different slot openings is not equal and so these forces do not compensate each other. In the calculation it is assumed that the flux entering the left side of the opening is the component left from the geometrical centre of the slot. This torque component together with the torque component calculated using the Lorenz force make the total torque of the motor. It is easy to assume that when all the magnet edges, where the derivative component of the magnet flux density is at its highest, enter the slot openings at the same time, this will have as a result a considerable cogging torque. To reduce the cogging torquethe magnet edges can be shaped so that they are not parallel to the stator slots, which is the common way to solve the problem. In doing so, the edge may be spread along the whole slot pitch and thus also the high derivative component willbe spread to occur equally along the rotation. Besides forming the magnets theymay also be placed somewhat asymmetric on the rotor surface. The asymmetric distribution can be made in many different ways. All the magnets may have a different deflection of the symmetrical centre point or they can be for example shiftedin pairs. There are some factors that limit the deflection. The first is that the magnets cannot overlap. The magnet shape and the relative width compared to the pole define the deflection in this case. The other factor is that a shifting of the poles limits the maximum torque of the motor. If the edges of adjacent magnets are very close to each other the leakage flux from one pole to the other increases reducing thus the air-gap magnetization. The asymmetric model needs some assumptions and simplifications in order to limit the size of the model and calculation time. The reluctance net is made for symmetric distribution. If the magnets are distributed asymmetrically the flux in the different pole pairs will not be exactly the same. Therefore, the assumption that the flux flows from the edges of the model to the next pole pairs, in the calculation model from one edgeto the other, is not correct. If it were wished for that this fact should be considered in multi-pole pair machines, this would mean that all the poles, in other words the whole machine, should be modelled in reluctance net. The error resulting from this wrong assumption is, nevertheless, irrelevant.

Dynamic model of the left ventricle for use in simulation of myocardial perfusion SPECT and gated SPECT

Simulation is a useful tool in cardiac SPECT to assess quantification algorithms. However, simple equation-based models are limited in their ability to simulate realistic heart motion and perfusion. We present a numerical dynamic model of the left ventricle, which allows us to simulate normal and anomalous cardiac cycles, as well as perfusion defects. Bicubic splines were fitted to a number of control points to represent endocardial and epicardial surfaces of the left ventricle. A transformation from each point on the surface to a template of activity was made to represent the myocardial perfusion. Geometry-based and patient-based simulations were performed to illustrate this model. Geometry-based simulations modeled ~1! a normal patient, ~2! a well-perfused patient with abnormal regional function, ~3! an ischaemic patient with abnormal regional function, and ~4! a patient study including tracer kinetics. Patient-based simulation consisted of a left ventricle including a realistic shape and motion obtained from a magnetic resonance study. We conclude that this model has the potential to study the influence of several physical parameters and the left ventricle contraction in myocardial perfusion SPECT and gated-SPECT studies.

Murine model of pneumonia caused by Parachlamydia acanthamoebae

The role of Parachlamydia acanthamoebae as an agent of pneumonia is suggested by sero-epidemiological studies, molecular surveys and by the permissivity of macrophages, lung fibroblasts and pneumocytes to this obligate intracellular bacteria. We thus developed a murine model of pneumonia due to Parachlamydia. Mice were inoculated intratracheally with Parachlamydia acanthamoebae. Pneumonia-associated mortality was of 50% 5 days post-inoculation. Lungs histopathology was characterized by purulent and interstitial pneumonia. The presence of Parachlamydia in the lesions was demonstrated by PCR, immunohistochemistry and electron microscopy. Moreover, living Parachlamydia could be recovered from the lungs of infected mice using amoebal co-culture. All control mice inoculated with heat-inactivated bacteria were free of symptoms and survived. Thus, we demonstrated that Parachlamydia induce a severe pneumonia in mice. This animal model, which confirms the third and fourth Koch postulates, may be suitable to test in vivo efficient therapeutic regimens against Parachlamydia.

A bovine model of a respiratory Parachlamydia acanthamoebae infection.

The aim of this study was to evaluate the pathogenicity of Parachlamydia (P.) acanthamoebae as a potential agent of lower respiratory tract disease in a bovine model of induced lung infection. Intrabronchial inoculation with P. acanthamoebae was performed in healthy calves aged 2-3 months using two challenge doses: 10(8) and 10(10) bacteria per animal. Controls received 10(8) heat-inactivated bacteria. Challenge with 10(8) viable Parachlamydia resulted in a mild degree of general indisposition, whereas 10(10) bacteria induced a more severe respiratory illness becoming apparent 1-2 days post inoculation (dpi), affecting 9/9 (100%) animals and lasting for 6 days. The extent of macroscopic pulmonary lesions was as high as 6.6 (6.0)% [median (range)] of lung tissue at 2-4 dpi and correlated with parachlamydial genomic copy numbers detected by PCR, and with bacterial load estimated by immunohistochemistry in lung tissue. Clinical outcome, acute phase reactants, pathological findings and bacterial load exhibited an initial dose-dependent effect on severity. Animals fully recovered from clinical signs of respiratory disease within 5 days. The bovine lung was shown to be moderately susceptible to P. acanthamoebae, exhibiting a transient pneumonic inflammation after intrabronchial challenge. Further studies are warranted to determine the precise pathophysiologic pathways of host-pathogen interaction.

Model of black hole evolution

We investigate a model where the quantum dynamics of black hole evaporation is determined by imposing a boundary on the apparent horizon with suitable boundary conditions. An unconventional scenario for the evolution emerges: only an insignificant fraction of energy of order (mG)-1 is radiated out; the outgoing wave carries a very small part of the quantum-mechanical information of the collapsed body, the bulk of the information remaining in the final stable black hole geometry.

Lack of GDAP1 induces neuronal calcium and mitochondrial defects in a knockout mouse model of charcot-marie-tooth neuropathy.

Mutations in GDAP1, which encodes protein located in the mitochondrial outer membrane, cause axonal recessive (AR-CMT2), axonal dominant (CMT2K) and demyelinating recessive (CMT4A) forms of Charcot-Marie-Tooth (CMT) neuropathy. Loss of function recessive mutations in GDAP1 are associated with decreased mitochondrial fission activity, while dominant mutations result in impairment of mitochondrial fusion with increased production of reactive oxygen species and susceptibility to apoptotic stimuli. GDAP1 silencing in vitro reduces Ca2+ inflow through store-operated Ca2+ entry (SOCE) upon mobilization of endoplasmic reticulum (ER) Ca2+, likely in association with an abnormal distribution of the mitochondrial network. To investigate the functional consequences of lack of GDAP1 in vivo, we generated a Gdap1 knockout mouse. The affected animals presented abnormal motor behavior starting at the age of 3 months. Electrophysiological and biochemical studies confirmed the axonal nature of the neuropathy whereas histopathological studies over time showed progressive loss of motor neurons (MNs) in the anterior horn of the spinal cord and defects in neuromuscular junctions. Analyses of cultured embryonic MNs and adult dorsal root ganglia neurons from affected animals demonstrated large and defective mitochondria, changes in the ER cisternae, reduced acetylation of cytoskeletal α-tubulin and increased autophagy vesicles. Importantly, MNs showed reduced cytosolic calcium and SOCE response. The development and characterization of the GDAP1 neuropathy mice model thus revealed that some of the pathophysiological changes present in axonal recessive form of the GDAP1-related CMT might be the consequence of changes in the mitochondrial network biology and mitochondria-endoplasmic reticulum interaction leading to abnormalities in calcium homeostasis.