Combinatorial Anti-arenaviral Therapy with the Small Molecule SKI-1/S1P Inhibitor PF-429242 and Ribavirin

Arenaviruses are enveloped negative strand viruses that cause acute and chronic infections. Several Arenaviruses can cause severe hemorrhagic fever in humans. In West Africa Lassa virus causes several hundred thousand infections per year, while Junin, Machupo, Guanarito, and Sabia virus have emerged in South America. So far, only one drug is licensed against arenaviruses, the nucleoside analogue Ribavirin (Rib), which is effective when given early in disease, but shows only minor therapeutic effects in late stages of the infection. Previous works demonstrated that processing of the arenavirus glycoprotein precursor (GPC) by the cellular proprotein convertase site 1 protease (S1P), also known as subtilisin-kexinisozyme 1 (SKI-1), is crucial for cell-to-cell propagation of infectionand production of infectious virus. Recently, the SKI-1/S1P inhibitor PF-429242wasshownto inhibit Old World arenavirusGPCprocessing, cell-to-cell propagation, and infectious virus production. In the present study, we assessed the activity of PF-429242 against processing of the GPCs of the genetically and structurally more distant New World arenaviruses and found potent inhibition of processing of the GPCs of Junin, Machupo, and Guanarito virus. Using the prototypic arenavirus lymphocytic choriomeningitis virus (LCMV), we studied the potency of PF-429242 in the context of acute and chronic infection. In line with published data, PF-429242 potently inhibited acute LCMV infection. PF-429242 was also highly active against chronic infection and drug treatment resulted in rapid extinction of the virus without emergence of drug-resistant variants. In a combinatorial drug approach, we found that PF-429242 potentiated the anti-viral effect of Rib in treatment of acute andchronic infection. Taken together, we showed that the SKI-1/S1P inhibitor PF-429242 is broadly active against GPC processing of all major human pathogenic arenaviruses. Apart from being potent in acute infection, the drug is remarkably active in clearing chronic infection and potentiated the anti-arenaviral activity of Rib.

Role and regulation of polyphosphate in leishmania parasites

Le protozoaire unicellulaire Leishmania est l'agent responsable de la leishmaniose, une maladie parasitaire humaine qui se manifeste par des lésions de la peau, se résolvant le plus souvent spontanément, jusqu'à des lésions viscérales fatales. Le parasite est transmis de l'insecte à l'hôte mammifère lors d'un repas sanguin de la mouche des sables et y réside respectivement sous formes extra- et intracellulaires. On estime que cette maladie touche environ 12 millions de personnes dans 98 pays. Etant donné que les médicaments disponibles à ce jour sont faiblement efficaces et/ou hautement toxiques, il est indispensable de consolider les connaissances sur le fonctionnement et la survie du parasite pour pouvoir développer de nouvelles stratégies de traitements et de préventions. Tous les organismes vivants, dont Leishmania, contiennent du polyphopshate (polyP). Cette molécule chargée négativement est constituée de trois jusqu'à plusieurs centaines de résidus de phosphates reliées par des liaisons à haute énergie. Le polyP sert donc de source d'énergie et de réservoir de phosphate; dans certaines espèces, il joue aussi un rôle dans l'adaptation au stress et la virulence de pathogènes. Ceci nous a amené à étudier le rôle du polyP dans le parasite Leishmania. L'enzyme responsable de la synthèse de polyP a été identifié récemment dans la levure : il s'agit de la chaperone de transport vacuolaire 4 (Vtc4). Nous avons identifié un homologue de Vtc4 chez les Trypanosomatidae, et avons donc décidé d'examiner sa fonction dans le métabolisme du polyP chez Leishmania. En éliminant l'expression de Vtc4 chez L. major et L. guyanensis, nous avons pu démontrer qu'il est indispensable pour la production de polyP chez Leishmania. De plus, nous avons constaté que ces parasites possèdent des chaînes de polyP allant de trois jusqu'à environ 300 résidus de phosphate. Le taux de polyP dans la cellule est précisément régulé et varie entre un très haut niveau durant la phase proliférative des promastigotes à un niveau bas en phase stationnaire tardive, alors que l'expression de Vtc4p reste stable. Dans les amastigotes intracellulaires, seulement des petites quantités de polyP et de Vtc4p sont détectées. En outre, l'absence de Vc4p et de polyP n'a pas d'effet significatif sur les infections in vivo de souris, ce qui indique que le polyP n'est pas nécessaire au développement de la leishmaniose. Ceci suggère que Vtc4p n 'est pas une bonne cible pour le développement de nouveaux traitements contre Leishmania. "Néanmoins, la présence du polyP favorise fortement la survie du parasite suite à un choc de température (37°C) et aide ainsi à sa persistance intracellulaire pendant les premiers jours d'infection de macrophages. En résumé, nos résultats indiquent que si le polyP a peu d'importance pendant l'infection et le développement de la leishmaniose chez la souris, il est par contre crucial pour l'adaptation à des situations de stress comme l'augmentation de la température. Le fait que le polyP a été conservé dans tous les organismes durant l'évolution suggère toutefois que cette molécule joue un rôle fondamental. Etant donné que l'absence de polyP n'a pas d'effet sur la survie des amastigotes, il pourrait être plus important dans la forme promastigote infectant la mouche des sables. - The unicellular protozoan parasite Leishmania is the causative agent of the human disease leishmaniasis, which can range from self-healing skin lesions to fatal visceral lesions. The parasite is transmitted from the insect vector to the mammalian host when the sand fly takes its blood meal and exists in an extra- and an intracellular form, respectively. The disease is estimated to affect 12 million people in 98 countries and currently available drug treatments are of relatively low potency and/or high toxicity. Thus, investigating parasite survival mechanisms and parasite adaptation to the two host environments contributes to the general understanding of Leishmania propagation and might therefore help to develop future treatments or preventions. All living cells, including Leishmania, contain a negatively charged polymer of a few up to several hundred phosphate residues. These so-called polyphosphates (polyPs) serve as an energy source and phosphate reservoir. In some organisms, polyP is also involved in adaptation to stresses and virulence of pathogens. Therefore we were interested in investigating the importance of polyP in Leishmania parasites. Recently, an eukaryotic enzyme responsible for polyP synthesis has been identified as the vacuolar transporter chaperone 4 (Vtc4) in yeast. We, and others, found a Vtc4 homologue in trypanosomatids and decided to examine its potential function in polyP metabolism. By generating VTC4 knock-out cell lines in L. major and Vtc4 knock-down cell lines in L. guyanensis, we were able to demonstrate that Vtc4p is responsible for the total amount of cellular polyP. We also observed that Leishmania polyP chain length ranges from a few up to around 300 residues and that its level is tightly regulated. PolyP abundance is highest during the logarithmic proliferating phase of promastigotes and decreases in the stationary phase, while Vtc4 protein expression remains stable during both phases. In the intracellular amastigote form, only low amounts of polyP and Vtc4p were detectable. Furthermore, absence of Vtc4p and polyP did not have a significant effect on in vivo mouse infections, indicating that polyP is not necessary for Leishmania disease progression. This suggests that Vtc4p would be a poor drug target against Leishmania infection. However, presence of the polymer strongly supported parasite survival during heat shock (37°C) and thereby promoted intracellular persistence during the first days of macrophage infections. Taken together, we found that polyP has little importance in Leishmania {in vivo) infection but that it plays a crucial role during adaptation to stress, such as heat shock. Given that polyP has been preciously conserved in all organisms during evolution it seems to play a fundamental role. Since absence of polyP does not affect amastigote survival, it might be significant for promastigote existence in the sand fly vector.

Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations.

Brain fluctuations at rest are not random but are structured in spatial patterns of correlated activity across different brain areas. The question of how resting-state functional connectivity (FC) emerges from the brain's anatomical connections has motivated several experimental and computational studies to understand structure-function relationships. However, the mechanistic origin of resting state is obscured by large-scale models' complexity, and a close structure-function relation is still an open problem. Thus, a realistic but simple enough description of relevant brain dynamics is needed. Here, we derived a dynamic mean field model that consistently summarizes the realistic dynamics of a detailed spiking and conductance-based synaptic large-scale network, in which connectivity is constrained by diffusion imaging data from human subjects. The dynamic mean field approximates the ensemble dynamics, whose temporal evolution is dominated by the longest time scale of the system. With this reduction, we demonstrated that FC emerges as structured linear fluctuations around a stable low firing activity state close to destabilization. Moreover, the model can be further and crucially simplified into a set of motion equations for statistical moments, providing a direct analytical link between anatomical structure, neural network dynamics, and FC. Our study suggests that FC arises from noise propagation and dynamical slowing down of fluctuations in an anatomically constrained dynamical system. Altogether, the reduction from spiking models to statistical moments presented here provides a new framework to explicitly understand the building up of FC through neuronal dynamics underpinned by anatomical connections and to drive hypotheses in task-evoked studies and for clinical applications.

Comparison of Energy Minimization Methods for 3-D Brain Tissue Classification

This paper presents 3-D brain tissue classificationschemes using three recent promising energy minimizationmethods for Markov random fields: graph cuts, loopybelief propagation and tree-reweighted message passing.The classification is performed using the well knownfinite Gaussian mixture Markov Random Field model.Results from the above methods are compared with widelyused iterative conditional modes algorithm. Theevaluation is performed on a dataset containing simulatedT1-weighted MR brain volumes with varying noise andintensity non-uniformities. The comparisons are performedin terms of energies as well as based on ground truthsegmentations, using various quantitative metrics.

Neuronal activity in the hub of extrasynaptic Schwann cell-axon interactions.

The integrity and function of neurons depend on their continuous interactions with glial cells. In the peripheral nervous system glial functions are exerted by Schwann cells (SCs). SCs sense synaptic and extrasynaptic manifestations of action potential propagation and adapt their physiology to support neuronal activity. We review here existing literature data on extrasynaptic bidirectional axon-SC communication, focusing particularly on neuronal activity implications. To shed light on underlying mechanisms, we conduct a thorough analysis of microarray data from SC-rich mouse sciatic nerve at different developmental stages and in neuropathic models. We identify molecules that are potentially involved in SC detection of neuronal activity signals inducing subsequent glial responses. We further suggest that alterations in the activity-dependent axon-SC crosstalk impact on peripheral neuropathies. Together with previously reported data, these observations open new perspectives for deciphering glial mechanisms of neuronal function support.

Post-Plutonic Magmatism: from the Source to the Emplacement of an Upper Crustal Dyke Swarm, Adamello Massif, Italy

Magmas of the arc-tholeiitic and calc-alkaline differentiation suites contribute substantially to the formation of continental crust in subduction zones. Different geochemical-petrological models have been put forward to achieve evolved magmas forming large volumes of tonalitic to granitic plutons, building an important part of the continental crust. Primary magmas produced in the mantle wedge overlying the subducted slab migrate through the mantle and the crust. During the transfer, magma can accumulate in intermediate reservoirs at different levels where crystallization leads to differentiation and the heat transfer from the magma, together with gained heat from solidification, lead to partial melting of the crust. Partial melts can be assimilated and mix with more primitive magma. Moreover, already formed crystal cumulates or crystal mushes can be recycled and reactivated to transfer to higher crustal levels. Magma transport in the crust involves fow through fractures within a brittle elastic rock. The solidified magma filled crack, a dyke, can crosscut previously formed geological structures and thus serves as a relative or absolute time marker. The study area is situated in the Adamello massif. The Adamello massif is a composite of plutons that were emplaced between 42 and 29 million years. A later dyke swarm intruded into the southern part of the Adamello Batholith. A fractionation model covering dyke compositions from picrobasalts to dacites results in the cummulative crystallization of 17% olivine, 2% Cr-rich spinel, 18% clinopyroxene, 41% amphibole, 4% plagioclase and 0.1% magnetite to achieve an andesitic composition out of a hydrous primitive picrobasalt. These rocks show a similar geochemical evolution as experimental data simulating fractional crystallization and associated magma differentiation at lower crustal depth (7-10 kbar). The peraluminous, corundum normative composition is one characteristic of more evolved dacitic magmas, which has been explained in a long lasting debate with two di_erent models. Melting of mafic crust or politic material provides one model, whereas an alternative is fractionation from primary mantle derived melts. Amphibole occurring in basaltic-andesitic and andesitic dyke rocks as fractionating cumulate phase extracted from lower crustal depth (6-7.5 kbar) is driving the magmas to peraluminous, corundum normative compositions, which are represented by tonalites forming most of the Adamello Batholith. Most primitive picrobasaltic dykes have a slightly steepened chondrite normalized rare earth elements (REE) pattern and the increased enrichment of light-REE (LREE) for andesites and dacites can be explained by the fractional crystallization model originating from a picrobasalt, taking the changing fractionating phase assemblage and temperature into account. The injection of hot basaltic magma (~1050°C) in a closely spaced dyke swarm increases the surface of the contact to the mainly tonalitic wallrock. Such a setting induces partial melting of the wall rock and selective assimilation. Partial melting of the tonalite host is further expressed through intrusion breccias from basaltic dykes. Heat conduction models with instantaneous magma injection for such a dyke swarm geometry can explain features of partial melting observed in the field. Geochemical data of minerals and bulk rock further underline the selective or bulk assimilation of the tonalite host rock at upper crustal levels (~2-3 kbar), in particular with regard to light ion lithophile elements (LILE) such as Sr, Ba and Rb. Primitive picrobasalts carry an immiscible felsic assimilant as enclaves that bring along refractory rutile and zircon with textures typically found in oceanic plagiogranites or high pressure/low-temperature metamorphic rocks in general. U-Pb data implies a lower Cretaceous age for zircon not yet described as assimilant in Eocene to Oligocene magmatic rocks of the Central Southern Alps. The distribution of post-plutonic dykes in large batholiths such as the Adamello is one of the key features for understanding the regional stress field during the post-batholith emplacement cooling history. The emplacement of the regional dyke swarm covering the southern part of the Adamello massif was associated with consistent left lateral strike-slip movement along magma dilatation planes, leading to en echelon segmentation of dykes. Through the dilation by magma of pre-existing weaknesses and cracks in an otherwise uniform host rock, the dyke propagation and according orientation in the horizontal plane adjusted continuously perpendicular to least compressive remote stress σ3, resulting in an inferred rotation of the remote principal stress field. Les magmas issus des zones de subduction contribuent substantiellement à la formation de la croûte continentale. Les plutons tonalitiques et granitiques représentent, en effet, une partie importante de la croûte continentale. Des magmas primaires produits dans le 'mantle wedge ', partie du manteau se trouvant au-dessus de la plaque plongeante dans des zones de subduction, migrent à travers le manteau puis la croûte. Pendant ce transfert, le magma peut s'accumuler dans des réservoirs intermédiaires à différentes profondeurs. Le stockage de magma dans ces réservoirs engendre, d'une part, la différentiation des magmas par cristallisation fractionnée et, d'autre part, une fusion partielle la croûte continentale préexistante associée au transfert de la chaleur des magmas vers l'encaissant. Ces liquides magmatiques issus de la croûte peuvent, ensuite, se mélanger avec des magmas primaires. Le transport du magma dans la croûte implique notamment un flux de magma à travers différentes fractures recoupant les roches encaissantes élastiques. Au cours de ce processus de migration, des cumulats de cristaux ou des agrégats de cristaux encore non-solidifiés, peuvent être recyclés et réactivés pour être transportés à des niveaux supérieures de la croûte. Le terrain d'étude est situé dans le massif d'Adamello. Celui-ci est composé de plusieurs plutons mis en place entre 42 et 29 millions d'années. Dans une phase tardive de l'activité magmatique liée à ce batholite, une série de filons de composition variable allant de picrobasalte à des compositions dacitiques s'est mise en place la partie sud du massif. Deux modèles sont proposés dans la littérature, pour expliquer la formation des magmas dacitiques caractérisés par des compositions peralumineux (i.e. à corindon normatif). Le premier modèle propose que ces magmas soient issus de la fusion de matériel mafique et pélitique présent dans la partie inférieur de la croûte, alors que le deuxième modèle suggère une évolution par cristallisation fractionnée à partir de liquides primaires issus du manteau. Un modèle de cristallisation fractionnée a pu être développé pour expliquer l'évolution des filons de l'Adamello. Ce modèle explique la formation des filons dacitiques par la cristallisation fractionnée de 17% olivine, 2% spinelle riche en Cr, 18% clinopyroxène, 41% amphibole, 4% plagioclase et 0.1% magnetite à partir de liquide de compositions picrobasaltiques. Ce modèle prend en considération les contraintes pétrologiques déduites de l'observation des différents filons ainsi que du champ de stabilité des différentes phases en fonction de la température. Ces roches montrent une évolution géochimique similaire aux données expérimentales simulant la cristallisation fractionnée de magmas évoluant à des niveaux inférieurs de la croûte (7-10 kbar). Le modèle montre, en particulier, le rôle prépondérant de l'amphibole, une phase qui contrôle en particulier le caractère peralumineux des magmas différentiés ainsi que leurs compositions en éléments en traces. Des phénomènes de fusion partielle de l'encaissant tonalitique lors de la mise en place de _lons mafiques sont observée sur le terrain. L'injection du magma basaltique chaud (~1050°C) sous forme de filons rapprochés augmente la surface du contact avec l'encaissante tonalitique. Une telle situation produit la fusion partielle des roches encaissantes nécessaire à l'incorporation d'enclaves mafiques observés au sein des tonalites. Pour comprendre les conditions nécessaires pour la fusion partielle des roches encaissantes, des modèles de conduction thermique pour une injection simultanée d'une série de filons ont été développées. Des données géochimiques sur les minéraux et sur les roches totales soulignent qu'au niveau supérieur de la croûte, l'assimilation sélective ou totale de l'encaissante tonalitique modifie la composition du liquide primaire pour les éléments lithophiles tel que le Sr, Ba et Rb. Un autre aspect important concernant la pétrologie des filons de l'Adamello est la présence d'enclaves felsiques dans les filons les plus primitifs. Ces enclaves montrent, en particulier, des textures proches de celles rencontrées dans des plagiogranites océaniques ou dans des roches métamorphiques de haute pression/basse température. Ces enclaves contiennent du zircon et du rutile. La datations de ces zircons à l'aide du géochronomètre U-Pb indique un âge Crétacé inférieur. Cet âge est important, car aucune roche de cet âge n'a été considérée comme un assimilant potentiel pour des roches magmatiques d'âge Eocène à Oligocène dans les Alpes Sud Centrales. La réparation spatiale des filons post-plutoniques dans des grands batholites tel que l'Adamello, est une caractéristique clé pour la compréhension des champs de contraintes lors du refroidissement du batholite. L'orientation des filons va, en particulier, indiqué la contrainte minimal au sein des roches encaissante. La mise en place de la série de filon recoupant la partie Sud du massif de l'Adamello est associée à un décrochement senestre, un décrochement que l'on peut lié aux contraintes tectoniques régionales auxquelles s'ajoutent l'effet de la dilatation produite par la mise en place du batholite lui-même. Ce décrochement senestre produit une segmentation en échelon des filons.

Analyse des modifications de la chromatine impliquées dans l'effet barrière de CTF1 aux télomères humains

Eukaryotic genomes are compartmentalized in different structural domains that can affect positively or negatively gene expression. These regions of euchromatin and heterochromatin are characterized by distinct histones marks which can facilitate or repress gene transcription. The chromatin environment represents thus one of the main problems to control gene expression in biotechnological applications or gene therapy, since its expression is affected by the chromatin neighboring its locus of insertion. Some chromatin regions like telomeres are composed of constitutive heterochromatin which leads to the telomeric position effect (TPE) that silences genes adjacent to the telomere. TPE is known to spread by the selfrecruitment of the SIR histone deacetylase complex from the telomere in S.cerevisiae, but the histone marks that are associated to telomeric chromatin in mammalian cells remain mostly unknown. The transcription factor CTF1 has shown antisilencing properties in mammalian cells and also a boundary activity against TPE in yeast cells when fused to the yeast Gal4 DNA binding domain. In the work presented here, we describe a dual-reporter system to assess the boundary activity of proteins such as CTF1 at human telomeres. When located between the two reporter genes, CTF1 shields the telomere distal gene from TPE, while the telomereproximal gene remains silenced by telomeric heterochromatin. The boundary activity of CTF1 is shown to act regardless its function of transcriptional activator, by opposition to the transcriptional activator VP16 which activates indifferently both transgenes. Moreover, this study shows that CTF1 boundary activity is linked to its H3 binding function, as expected from a chromatin remodeler. ChIP experiments showed that histone deacetylation is the main histone modification involved in gene silencing at mammalian cell telomeres. Distinctly to yeast cells, the histone deacetylation signal in human cells extented over a short range along the chromosome. CTF1 may help to block this propagation and therefore to restore histones acetylation level on telomere protected locus. Surprisingly, other histone marks such as trimethyl-H3K9 or trimethyl-H4K20 were found on telomere protected locus, while in another clone, unsilencing of telomere distal transgene was associated with recruitment of the histone variant H2A.Z. Thus, I conclude that CTF1 displays a chromatin boundary function which is independent of its transcriptional activity and therefore exhibit features required for use as chromatin insulator in biotechnological applications. RESUME Les génomes eucaryotes sont compartementalisés en domaines structurels qui peuvent affecter positivement ou négativement l'expression des gènes avoisinants. Ces régions dites d'euchromatine ou d'hétérochromatine sont caractérisées par des modifications posttraductionnelles des histones qui peuvent faciliter ou au contraire inhiber la transcription des gènes qui s'y trouvent. Ainsi, isoler un gène de son environnement chromatinien est problème fréquent lorsqu'il s'agit de contrôler son expression dans le cadre d'applications en biotechnologie ou encore en thérapie génique. Certaines régions de chromatine telles que les télomères sont composées d'hétérochromatine constitutive qui mène au silençage des gènes avoisinants. Cet effet de position télomérique (TPE) est connu dans la levure S.cerevisiae comme se propageant par auto-recrutement du complexe de déacétylation d'histone SIR, alors que peu de modifications de chromatine ont pu être associées à ce phénomène dans les cellules de mammifères. Le facteur de transcription CTF1 a montré des propriétés d'anti-silençage dans les cellules de mammifères, ainsi qu'une activité barrière contre le silençage télomérique dans les cellules de levures lorsqu'il est fusionné au domaine de liaison à l'ADN de la protéine de levure Gal4. Dans le travail présenté ci-après est décrit un système à deux gènes rapporteurs permettant de mesurer l'activité barrière de protéines telles que CTF1 aux télomères humains, et les modifications de chromatine qui y sont associées. Lorsque CTF1 est placé entre les deux gènes rapporteurs, le gène distant du télomère est protégé du silençage qui lui est associé, alors que le gène proche du télomère reste soumis à ce silençage induit par l'hétérochromatine télomérique. L'activité barrière de CTF1 est montrée ici comme agissant indépendamment de son activité transcriptionnelle, par opposition à l'activateur transcriptionnel VP16 qui active indifféremment les deux transgènes. En outre, cette étude appuie l'hypothèse stipulant que CTF1 agisse comme remodeleur chromatinien puisqu'elle démontre que son activité barrière est directement dépendante de son activité de liaison avec l'histone H3. De plus, des expériences d'immuno-précipitation de la chromatine démontrent que la déacétylation des histones est le majeur phénomène intervenant dans le silençage télomérique. Par opposition à la levure, ce signal de déacétylation ne se propage dans les cellules humaines que sur une courte distance le long du chromosome. CTF1 agit ainsi en bloquant cette propagation et en restaurant le niveau d'acétylation des histones sur le locus protégé du télomère. De manière surprenante et inattendue, d'autres modifications d'histones telles que 4 les H3K9 et H4K20 triméthylées sont aussi observées à ce locus, tandis le recrutement du variant H2A.Z peut aussi être suffisant à restaurer l'expression du gène distant du télomère. En terme de cette analyse, CTF1 exhibe ainsi une fonction de barrière chromatinienne qui exclue une activité transcriptionnelle non désirée - propriété qui est requise dans l'établissement des isolateurs visant à permettre le contrôle d'un transgène dans le cadre d'applications en biotechnologies.

An active strain electromechanical model for cardiac tissue

We propose a finite element approximation of a system of partial differential equations describing the coupling between the propagation of electrical potential and large deformations of the cardiac tissue. The underlying mathematical model is based on the active strain assumption, in which it is assumed that a multiplicative decomposition of the deformation tensor into a passive and active part holds, the latter carrying the information of the electrical potential propagation and anisotropy of the cardiac tissue into the equations of either incompressible or compressible nonlinear elasticity, governing the mechanical response of the biological material. In addition, by changing from an Eulerian to a Lagrangian configuration, the bidomain or monodomain equations modeling the evolution of the electrical propagation exhibit a nonlinear diffusion term. Piecewise quadratic finite elements are employed to approximate the displacements field, whereas for pressure, electrical potentials and ionic variables are approximated by piecewise linear elements. Various numerical tests performed with a parallel finite element code illustrate that the proposed model can capture some important features of the electromechanical coupling, and show that our numerical scheme is efficient and accurate.

Propagation velocity kinetics and repolarization alternans in a free-behaving sheep model of pacing-induced atrial fibrillation.

AIMS: Experimental models have reported conflicting results regarding the role of dispersion of repolarization in promoting atrial fibrillation (AF). Repolarization alternans, a beat-to-beat alternation in action potential duration, enhances dispersion of repolarization when propagation velocity is involved. METHODS AND RESULTS: In this work, original electrophysiological parameters were analysed to study AF susceptibility in a chronic sheep model of pacing-induced AF. Two pacemakers were implanted, each with a single right atrial lead. Right atrial depolarization and repolarization waves were documented at 2-week intervals. A significant and gradual decrease in the propagation velocity at all pacing rates and in the right atrial effective refractory period (ERP) was observed during the weeks of burst pacing before sustained AF developed when compared with baseline conditions. Right atrial repolarization alternans was observed, but because of the development of 2/1 atrioventricular block with far-field ventricular interference, its threshold could not be precisely measured. Non-sustained AF was not observed at baseline, but appeared during the electrical remodelling in association with a decrease in both ERP and propagation velocity. CONCLUSION: We report here on the feasibility of measuring ERP, atrial repolarization alternans, and propagation velocity kinetics and their potential in predicting susceptibility to AF in a free-behaving model of pacing-induced AF using the standard pacemaker technology.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Simulations in evolution. II. Relative fitness and the propagation of mutants.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous article presented a histogram model [1] consisting in populations of individuals whose number changed under the influence of variation and/or fitness, the total population remaining constant. Individuals are classified into bins, and the content of each bin is calculated generation after generation by an Excel spreadsheet. Here, we apply the histogram model to a stable population with fitness F(1)=1.00 in which one or two fitter mutants emerge. In a first scenario, a single mutant emerged in the population whose fitness was greater than 1.00. The simulations ended when the original population was reduced to a single individual. The histogram model was validated by excellent agreement between its predictions and those of a classical continuous function (Eqn. 1) which predicts the number of generations needed for a favorable mutation to spread throughout a population. But in contrast to Eqn. 1, our histogram model is adaptable to more complex scenarios, as demonstrated here. In the second and third scenarios, the original population was present at time zero together with two mutants which differed from the original population by two higher and distinct fitness values. In the fourth scenario, the large original population was present at time zero together with one fitter mutant. After a number of generations, when the mutant offspring had multiplied, a second mutant was introduced whose fitness was even greater. The histogram model also allows Shannon entropy (SE) to be monitored continuously as the information content of the total population decreases or increases. The results of these simulations illustrate, in a graphically didactic manner, the influence of natural selection, operating through relative fitness, in the emergence and dominance of a fitter mutant.

AltitudeOmics: The Integrative Physiology of Human Acclimatization to Hypobaric Hypoxia and Its Retention upon Reascent.

An understanding of human responses to hypoxia is important for the health of millions of people worldwide who visit, live, or work in the hypoxic environment encountered at high altitudes. In spite of dozens of studies over the last 100 years, the basic mechanisms controlling acclimatization to hypoxia remain largely unknown. The AltitudeOmics project aimed to bridge this gap. Our goals were 1) to describe a phenotype for successful acclimatization and assess its retention and 2) use these findings as a foundation for companion mechanistic studies. Our approach was to characterize acclimatization by measuring changes in arterial oxygenation and hemoglobin concentration [Hb], acute mountain sickness (AMS), cognitive function, and exercise performance in 21 subjects as they acclimatized to 5260 m over 16 days. We then focused on the retention of acclimatization by having subjects reascend to 5260 m after either 7 (n = 14) or 21 (n = 7) days at 1525 m. At 16 days at 5260 m we observed: 1) increases in arterial oxygenation and [Hb] (compared to acute hypoxia: PaO2 rose 9±4 mmHg to 45±4 while PaCO2 dropped a further 6±3 mmHg to 21±3, and [Hb] rose 1.8±0.7 g/dL to 16±2 g/dL; 2) no AMS; 3) improved cognitive function; and 4) improved exercise performance by 8±8% (all changes p<0.01). Upon reascent, we observed retention of arterial oxygenation but not [Hb], protection from AMS, retention of exercise performance, less retention of cognitive function; and noted that some of these effects lasted for 21 days. Taken together, these findings reveal new information about retention of acclimatization, and can be used as a physiological foundation to explore the molecular mechanisms of acclimatization and its retention.

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

AltitudeOmics: enhanced cerebrovascular reactivity and ventilatory response to CO2 with high-altitude acclimatization and reexposure.

The present study is the first to examine the effect of high-altitude acclimatization and reexposure on the responses of cerebral blood flow and ventilation to CO2. We also compared the steady-state estimates of these parameters during acclimatization with the modified rebreathing method. We assessed changes in steady-state responses of middle cerebral artery velocity (MCAv), cerebrovascular conductance index (CVCi), and ventilation (V(E)) to varied levels of CO2 in 21 lowlanders (9 women; 21 ± 1 years of age) at sea level (SL), during initial exposure to 5,260 m (ALT1), after 16 days of acclimatization (ALT16), and upon reexposure to altitude following either 7 (POST7) or 21 days (POST21) at low altitude (1,525 m). In the nonacclimatized state (ALT1), MCAv and V(E) responses to CO2 were elevated compared with those at SL (by 79 ± 75% and 14.8 ± 12.3 l/min, respectively; P = 0.004 and P = 0.011). Acclimatization at ALT16 further elevated both MCAv and Ve responses to CO2 compared with ALT1 (by 89 ± 70% and 48.3 ± 32.0 l/min, respectively; P < 0.001). The acclimatization gained for V(E) responses to CO2 at ALT16 was retained by 38% upon reexposure to altitude at POST7 (P = 0.004 vs. ALT1), whereas no retention was observed for the MCAv responses (P > 0.05). We found good agreement between steady-state and modified rebreathing estimates of MCAv and V(E) responses to CO2 across all three time points (P < 0.001, pooled data). Regardless of the method of assessment, altitude acclimatization elevates both the cerebrovascular and ventilatory responsiveness to CO2. Our data further demonstrate that this enhanced ventilatory CO2 response is partly retained after 7 days at low altitude.