Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

Oxygen uptake kinetics and middle distance swimming performance.

OBJECTIVE: The aim of this study was to determine whether V˙O(2) kinetics and specifically, the time constant of transitions from rest to heavy (τ(p)H) and severe (τ(p)S) exercise intensities, are related to middle distance swimming performance. DESIGN: Fourteen highly trained male swimmers (mean ± SD: 20.5 ± 3.0 yr; 75.4 ± 12.4 kg; 1.80 ± 0.07 m) performed an discontinuous incremental test, as well as square wave transitions for heavy and severe swimming intensities, to determine V˙O(2) kinetics parameters using two exponential functions. METHODS: All the tests involved front-crawl swimming with breath-by-breath analysis using the Aquatrainer swimming snorkel. Endurance performance was recorded as the time taken to complete a 400 m freestyle swim within an official competition (T400), one month from the date of the other tests. RESULTS: T400 (Mean ± SD) (251.4 ± 12.4 s) was significantly correlated with τ(p)H (15.8 ± 4.8s; r=0.62; p=0.02) and τ(p)S (15.8 ± 4.7s; r=0.61; p=0.02). The best single predictor of 400 m freestyle time, out of the variables that were assessed, was the velocity at V˙O(2max)vV˙O(2max), which accounted for 80% of the variation in performance between swimmers. However, τ(p)H and V˙O(2max) were also found to influence the prediction of T400 when they were included in a regression model that involved respiratory parameters only. CONCLUSIONS: Faster kinetics during the primary phase of the V˙O(2) response is associated with better performance during middle-distance swimming. However, vV˙O(2max) appears to be a better predictor of T400.

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings.

Mathematical methods combined with measurements of single-cell dynamics provide a means to reconstruct intracellular processes that are only partly or indirectly accessible experimentally. To obtain reliable reconstructions, the pooling of measurements from several cells of a clonal population is mandatory. However, cell-to-cell variability originating from diverse sources poses computational challenges for such process reconstruction. We introduce a scalable Bayesian inference framework that properly accounts for population heterogeneity. The method allows inference of inaccessible molecular states and kinetic parameters; computation of Bayes factors for model selection; and dissection of intrinsic, extrinsic and technical noise. We show how additional single-cell readouts such as morphological features can be included in the analysis. We use the method to reconstruct the expression dynamics of a gene under an inducible promoter in yeast from time-lapse microscopy data.

Fat oxidation kinetics during exercise in obese individuals

Introduction An impaired ability to oxidize fat may be a factor in the obesity's aetiology (3). Moreover, the exercise intensity (Fatmax) eliciting the maximal fat oxidation rate (MFO) was lower in obese (O) compared with lean (L) individuals (4). However, difference in fat oxidation rate (FOR) during exercise between O and L remains equivocal and little is known about FORs during high intensities (>60% ) in O compared with L. This study aimed to characterize fat oxidation kinetics over a large range of intensities in L and O. Methods 12 healthy L [body mass index (BMI): 22.8±0.4] and 16 healthy O men (BMI: 38.9±1.4) performed submaximal incremental test (Incr) to determine whole-body fat oxidation kinetics using indirect calorimetry. After a 15-min resting period (Rest) and 10-min warm-up at 20% of maximal power output (MPO, determined by a maximal incremental test), the power output was increased by 7.5% MPO every 6-min until respiratory exchange ratio reached 1.0. Venous lactate and glucose and plasma concentration of epinephrine (E), norepinephrine (NE), insulin and non-esterified fatty acid (NEFA) were assessed at each step. A mathematical model (SIN) (1), including three variables (dilatation, symmetry, translation), was used to characterize fat oxidation (normalized by fat-free mass) kinetics and to determine Fatmax and MFO. Results FOR at Rest and MFO were not significantly different between groups (p≥0.1). FORs were similar from 20-60% (p≥0.1) and significantly lower from 65-85% in O than in L (p≤0.04). Fatmax was significantly lower in O than in L (46.5±2.5 vs 56.7±1.9 % respectively; p=0.005). Fat oxidation kinetics was characterized by similar translation (p=0.2), significantly lower dilatation (p=0.001) and tended to a left-shift symmetry in O compared with L (p=0.09). Plasma E, insulin and NEFA were significantly higher in L compared to O (p≤0.04). There were no significant differences in glucose, lactate and plasma NE between groups (p≥0.2). Conclusion The study showed that O presented a lower Fatmax and a lower reliance on fat oxidation at high, but not at moderate, intensities. This may be linked to a: i) higher levels of insulin and lower E concentrations in O, which may induce blunted lipolysis; ii) higher percentage of type II and a lower percentage of type I fibres (5), and iii) decreased mitochondrial content (2), which may reduce FORs at high intensities and Fatmax. These findings may have implications for an appropriate exercise intensity prescription for optimize fat oxidation in O. References 1. Cheneviere et al. Med Sci Sports Exerc. 2009 2. Holloway et al. Am J Clin Nutr. 2009 3. Kelley et al. Am J Physiol. 1999 4. Perez-Martin et al. Diabetes Metab. 2001 5. Tanner et al. Am J Physiol Endocrinol Metab. 2002

Viral-mediated overexpression of mutant huntingtin to model HD in various species.

Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by an expansion of CAG repeats in the huntingtin (Htt) gene. Despite intensive efforts devoted to investigating the mechanisms of its pathogenesis, effective treatments for this devastating disease remain unavailable. The lack of suitable models recapitulating the entire spectrum of the degenerative process has severely hindered the identification and validation of therapeutic strategies. The discovery that the degeneration in HD is caused by a mutation in a single gene has offered new opportunities to develop experimental models of HD, ranging from in vitro models to transgenic primates. However, recent advances in viral-vector technology provide promising alternatives based on the direct transfer of genes to selected sub-regions of the brain. Rodent studies have shown that overexpression of mutant human Htt in the striatum using adeno-associated virus or lentivirus vectors induces progressive neurodegeneration, which resembles that seen in HD. This article highlights progress made in modeling HD using viral vector gene transfer. We describe data obtained with of this highly flexible approach for the targeted overexpression of a disease-causing gene. The ability to deliver mutant Htt to specific tissues has opened pathological processes to experimental analysis and allowed targeted therapeutic development in rodent and primate pre-clinical models.

Viral superantigen drives extrafollicular and follicular B cell differentiation leading to virus-specific antibody production.

Mouse mammary tumor virus (MMTV[SW]) encodes a superantigen expressed by infected B cells. It evokes an antibody response specific for viral envelope protein, indicating selective activation of antigen-specific B cells. The response to MMTV(SW) in draining lymph nodes was compared with the response to haptenated chicken gamma globulin (NP-CGG) using flow cytometry and immunohistology. T cell priming occurs in both responses, with T cells proliferating in association with interdigitating dendritic cells in the T zone. T cell proliferation continues in the presence of B cells in the outer T zone, and B blasts then undergo exponential growth and differentiation into plasma cells in the medullary cords. Germinal centers develop in both responses, but those induced by MMTV(SW) appear later and are smaller. Most T cells activated in the T zone and germinal centers in the MMTV(SW) response are superantigen specific and these persist for weeks in lymph nodes draining the site MMTV(SW) injection: this contrasts with the selective loss of superantigen-specific T cells from other secondary lymphoid tissues. The results indicate that this viral superantigen, when expressed by professional antigen-presenting cells, drives extrafollicular and follicular B cell differentiation leading to virus-specific antibody production.

Purification and ligand binding of a soluble class I major histocompatibility complex molecule consisting of the first three domains of H-2Kd fused to beta 2-microglobulin expressed in the baculovirus-insect cell system.

A recombinant baculovirus encoding a single-chain murine major histocompatibility complex class I molecule in which the first three domains of H-2Kd are fused to beta 2-microglobulin (beta 2-m) via a 15-amino acid linker has been isolated and used to infect lepidopteran cells. A soluble, 391-amino acid single-chain H-2Kd (SC-Kd) molecule of 48 kDa was synthesized and glycosylated in insect cells and could be purified in the absence of detergents by affinity chromatography using the anti-H-2Kd monoclonal antibody SF1.1.1.1. We tested the ability of SC-Kd to bind antigenic peptides using a direct binding assay based on photoaffinity labeling. The photoreactive derivative was prepared from the H-2Kd-restricted Plasmodium berghei circumsporozoite protein (P.b. CS) peptide 253-260 (YIPSAEKI), a probe that we had previously shown to be unable to bind to the H-2Kd heavy chain in infected cells in the absence of co-expressed beta 2-microglobulin. SC-Kd expressed in insect cells did not require additional mouse beta 2-m to bind the photoprobe, indicating that the covalently attached beta 2-m could substitute for the free molecule. Similarly, binding of the P.b. CS photoaffinity probe to the purified SC-Kd molecule was unaffected by the addition of exogenous beta 2-m. This is in contrast to H-2KdQ10, a soluble H-2Kd molecule in which beta 2-m is noncovalently bound to the soluble heavy chain, whose ability to bind the photoaffinity probe is greatly enhanced in the presence of an excess of exogenous beta 2-m. The binding of the probe to SC-Kd was allele-specific, since labeling was selectively inhibited only by antigenic peptides known to be presented by the H-2Kd molecule.

Encéphalites virales [Viral encephalitis: update]

Suspicion of viral encephalitis should always be considered as a medical emergency and the prognosis depend on both the immune status of the host and the virulence of the virus. Among them, the herpes simplex virus is by far the most important one since it can be associated with severe encephalitis in immunocompetent host, and because a good response to acyclovir can be expected when rapidly initiated. Nevertheless, confirmation of the diagnosis requires exclusion of both metabolic or toxic encephalopathy and inflammatory encephalitis of non-infectious origin. In addition, other germs than viruses can mimic viral encephalitis and must be taken into consideration. The purpose of this review is to update the investigation that should be performed in clinical practice for any patient with suspicion of acute viral encephalitis.

Siglec-1 is a novel dendritic cell receptor that mediates HIV-1 trans-infection through recognition of viral membrane gangliosides.

Dendritic cells (DCs) are essential antigen-presenting cells for the induction of immunity against pathogens. However, HIV-1 spread is strongly enhanced in clusters of DCs and CD4(+) T cells. Uninfected DCs capture HIV-1 and mediate viral transfer to bystander CD4(+) T cells through a process termed trans-infection. Initial studies identified the C-type lectin DC-SIGN as the HIV-1 binding factor on DCs, which interacts with the viral envelope glycoproteins. Upon DC maturation, however, DC-SIGN is down-regulated, while HIV-1 capture and trans-infection is strongly enhanced via a glycoprotein-independent capture pathway that recognizes sialyllactose-containing membrane gangliosides. Here we show that the sialic acid-binding Ig-like lectin 1 (Siglec-1, CD169), which is highly expressed on mature DCs, specifically binds HIV-1 and vesicles carrying sialyllactose. Furthermore, Siglec-1 is essential for trans-infection by mature DCs. These findings identify Siglec-1 as a key factor for HIV-1 spread via infectious DC/T-cell synapses, highlighting a novel mechanism that mediates HIV-1 dissemination in activated tissues.

Cloning of a cDNA encoding a plasma membrane-associated, uronide binding phosphoprotein with physical properties similar to viral movement proteins.

Oligogalacturonides are structural and regulatory homopolymers from the extracellular pectic matrix of plants. In vitro micromolar concentrations of oligogalacturonates and polygalacturonates were shown previously to stimulate the phosphorylation of a small plasma membrane-associated protein in potato. Immunologically cross-reactive proteins were detected in plasma membrane-enriched fractions from all angiosperm subclasses in the Cronquist system. Polygalacturonate-enhanced phosphorylation of the protein was observed in four of the six dicotyledon subclasses but not in any of the five monocotyledon subclasses. A cDNA for the protein was cloned from potato. The deduced protein is extremely hydrophilic and has a proline-rich N terminus. The C-terminal half of the protein was predicted to be a coiled coil, suggesting that the protein interacts with other macromolecules. The recombinant protein was found to bind both simple and complex galacturonides. The behavior of the protein suggests several parallels with viral proteins involved in intercellular communication.