Sex-biased dispersal in a migratory bat: a characterization using sex-specific demographic parameters.

We studied the noctule bat (Nyctalus noctula), in which the mitochondrial F(ST) is about 10 times that revealed by nuclear markers, to address two questions. We first verified whether random dispersal of one sex is compatible with highly contrasted mitochondrial and nuclear population structures. Using computer simulations, we then assessed the power of multilocus population differentiation tests when the expected population structure departs only slightly from panmixia. Using an island model with sex-specific demographic parameters, we found that random male dispersal is consistent with the population structure observed in the noctule. However, other parameter combinations are also compatible with the data. We computed the minimum sex bias in dispersal (at least 69% of the dispersing individuals are males), a result that would not be available if we had used more classical population genetic models. The power of multilocus population differentiation tests was unexpectedly high, the tests being significant in almost 100% of the replicates, although the observed population structure infered from nuclear markers was extremely low (F(ST) = 0.6%).

Robust synchronization of coupled circadian and cell cycle oscillators in single mammalian cells.

Circadian cycles and cell cycles are two fundamental periodic processes with a period in the range of 1 day. Consequently, coupling between such cycles can lead to synchronization. Here, we estimated the mutual interactions between the two oscillators by time-lapse imaging of single mammalian NIH3T3 fibroblasts during several days. The analysis of thousands of circadian cycles in dividing cells clearly indicated that both oscillators tick in a 1:1 mode-locked state, with cell divisions occurring tightly 5 h before the peak in circadian Rev-Erbα-YFP reporter expression. In principle, such synchrony may be caused by either unidirectional or bidirectional coupling. While gating of cell division by the circadian cycle has been most studied, our data combined with stochastic modeling unambiguously show that the reverse coupling is predominant in NIH3T3 cells. Moreover, temperature, genetic, and pharmacological perturbations showed that the two interacting cellular oscillators adopt a synchronized state that is highly robust over a wide range of parameters. These findings have implications for circadian function in proliferative tissues, including epidermis, immune cells, and cancer.

Interregional compensatory mechanisms of motor functioning in progressing preclinical neurodegeneration.

Understanding brain reserve in preclinical stages of neurodegenerative disorders allows determination of which brain regions contribute to normal functioning despite accelerated neuronal loss. Besides the recruitment of additional regions, a reorganisation and shift of relevance between normally engaged regions are a suggested key mechanism. Thus, network analysis methods seem critical for investigation of changes in directed causal interactions between such candidate brain regions. To identify core compensatory regions, fifteen preclinical patients carrying the genetic mutation leading to Huntington's disease and twelve controls underwent fMRI scanning. They accomplished an auditory paced finger sequence tapping task, which challenged cognitive as well as executive aspects of motor functioning by varying speed and complexity of movements. To investigate causal interactions among brain regions a single Dynamic Causal Model (DCM) was constructed and fitted to the data from each subject. The DCM parameters were analysed using statistical methods to assess group differences in connectivity, and the relationship between connectivity patterns and predicted years to clinical onset was assessed in gene carriers. In preclinical patients, we found indications for neural reserve mechanisms predominantly driven by bilateral dorsal premotor cortex, which increasingly activated superior parietal cortices the closer individuals were to estimated clinical onset. This compensatory mechanism was restricted to complex movements characterised by high cognitive demand. Additionally, we identified task-induced connectivity changes in both groups of subjects towards pre- and caudal supplementary motor areas, which were linked to either faster or more complex task conditions. Interestingly, coupling of dorsal premotor cortex and supplementary motor area was more negative in controls compared to gene mutation carriers. Furthermore, changes in the connectivity pattern of gene carriers allowed prediction of the years to estimated disease onset in individuals. Our study characterises the connectivity pattern of core cortical regions maintaining motor function in relation to varying task demand. We identified connections of bilateral dorsal premotor cortex as critical for compensation as well as task-dependent recruitment of pre- and caudal supplementary motor area. The latter finding nicely mirrors a previously published general linear model-based analysis of the same data. Such knowledge about disease specific inter-regional effective connectivity may help identify foci for interventions based on transcranial magnetic stimulation designed to stimulate functioning and also to predict their impact on other regions in motor-associated networks.

Perioperative nutrition is still a surgical orphan

Objective: To assess if screening programs and treatment of preoperative malnutrition have been implemented into surgical practice to decrease morbidity. There is strong evidence that postoperative morbidity can be minimized by early identifying and treating patients at nutritional risk before major surgery.The validated nutritional risk score (NRS) is recommended by the European Society of Parenteral and Enteral Nutrition for nutritional screening. It remains unclear whether routine preoperative nutritional assessment and perioperative nutrition is widely implemented.Methods: A survey was conducted in 173 Swiss and Austrian surgical departments. Implementation of nutritional screening, perioperative nutrition, and estimated impact on clinical outcome were assessed. Non-responders were repeatedly contacted by the authors.Results: The overall response rate was 55%, whereby 69% (54/78) of Swiss and 44% (42/95) of Austrian centers responded. Despite 80% and 59% of the responding centers are aware of a reduced complication rate and shortened hospital stay, respectively, only 20% of them implemented routine nutritional screening. Financial (49%) and logistic restrictions (33%) are the predominant reasons against the routine clinical use. Screening is mainly performed either in the outpatient's clinic (52%) or during admission (54%). The NRS is only used by 14%. Instead, various clinical (78%), e.g. BMI and laboratory findings (56%), e.g. albumine, are used. Indication for perioperative nutrition is based on preoperative screening in 49%.While 23% use preoperative nutrition, 68% apply nutritional support pre- and postoperatively. Preoperative nutritional treatment ranged from three days (33%), to five days (31%) and even seven days (20%).Conclusion: Despite malnutrition is well recognized as major risk factor for increased postoperative morbidity, the majority of surgeons are reluctant to implement routine screening and nutritional support. If nutritional assessment is performed, local institutional screening parameters are still preferred. It remains difficult to overcome traditions, and to change surgeon's mind.

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

Effect of oxidized cholesterol on age-associated changes to immune parameters in spleen lymphocytes and peritoneal exudate cells derived from rats

The effects of oxidized cholesterol on immune parameters were examined by using spleen lymphocytes and peritoneal exudate cells (PEC) derived from 5-week- (Young) and 9-month-old (Adult) rats. The immunoglobulin (Ig) G and IgM production was inhibited by oxidized cholesterol in the rats of both ages when lymphocytes were exposed to 30 micrograms/ml of oxidized cholesterol for 24 hr. The intracellular IgA level was also lowered by 30 micrograms/ml of oxidized cholesterol, irrespective of age. In contrast, IgE production was significantly increased by the addition of 30 micrograms/ml of oxidized cholesterol in only young lymphocytes. Moreover, oxidized cholesterol enhanced the intracellular histamine accumulation in only adult PEC, although the total histamine level produced by PEC was similar in the rats of both ages. These results thus suggest the possibility that oxidized cholesterol can have different effects on the age-related modulation of immune functions such as Igs production and histamine release.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Contribution pollinique à l'étude du pléistocène de la région lémanique (Suisse) : paléoenvironnement des derniers 800'000 ans

Abstact : Pollinic contribution to a Pleistocene study of the Lake Geneva region (Switzerland) Palaeoenvironment of the last 800,000 years. Although the "Plateau Romand" (tableland of French-speaking Switzerland) has been marked by the repetitive action of different Quaternary glacial thrusts, numerous preserved outcrops here and there, thanks to favourable topographical conditions, provide sufficient pollinic content to reconstitute the significant botanical fluctuations and allow bio- and chrono-stratigraphical correlations with the long European sequences. The present study examines around ten Pleistocene deposits rich in pollens, whose presumed ages cover a range of 800,000 years. - The Early and Early-Middle Pleistocene are examined at Ecoteaux (VD), undeniably the oldest site in this study. The "Lower Formation of Ecoteaux", whose low pollinic content reveals an environment of the cold desert type, shows an inverse palaeomagnetic remanence. An age earlier than 780,000 years BP (Matuyama period or earlier) is proposed. The "Upper Formation of Ecoteaux" contains over eight major fluctuations in the vegetal covering, of which at least four temporal phases and one climatic optimum. The presence of botanical genera, relics of the Tertiary in particular Pterocarya and Carya, associated with the vegetal dynamic made up of short oscillations, contrives to link this formation to the Cromerian Complex. The thermomeres recorded could correspond to one or several interglacial periods from Middle Pleistocene. Possible correlations with marine isotopic stages 15, 17 and 19 are discussed. - The Middle and Late Pleistocene is approached through the sediments of Onnens (VD), Creux d'Enfer (FR), Port-Valais (VS), La Dénériaz (VD) and Cortaillod (NE). Each site provides a temperate interglacial flora, whose pollinic content is evaluated according to the different conceivable bio-stratigraphic correlations. - Finally the Würm is studied thanks to the outcrop of Marly (FR), together with the excavations at Villars-sous-Yens (VD), Versoix (GE) and Ollon (VD). Together, these observations tend to show that the west of the Swiss Plateau is a zone of transition between the meridian vegetation of the lower Rhone valley and that of beyond the Jura. At the outlet of the Rhone glacier, and on the borders of the prealpine and Jurassian glaciers, the region is situated at the crossroads of post-glacial botanical migrations, between the major axis of the Rhone valley and the Eastern refuge zones; this situation confers particular importance on it. Although believed to be deeply scraped by ice, one notices today that the "Plateau Romand" harbours overtwenty interglacial and interstadial sites. The abrasive action of the Rhone glacier has thus been more modest than previously estimated. The effective extension of glacial tongues, together with atmospheric circulation, are doubtless varying from one interglacial to another. Finally, the speeds of glacial movements, and the unbelievable vegetal vitality, are important parameters, capable of leading to deposits, very early enriched in pollens. RÉSUMÉS : Contribution pollinique à l'étude du Pléistocène de la région lémanique (Suisse) - Paléoenvironnement des derniers 800'000 ans Bien que le Plateau romand soit marqué par l'action répétitive des différentes poussées glaciaires quaternaires, de nombreux affleurements préservés ça et là, à la faveur de conditions topographiques favorables, livrent un contenu pollinique suffisant pour reconstituer des fluctuations botaniques significatives et autoriser des corrélations bioet chronostratigraphiques avec les longues séquences européennes. Le présent travail étudie une dizaine de dépôts pléistocènes riches en pollens, dont les âges présumés balaient un éventail de plus de 800'000 ans. - Le Pléistocène Ancien et Moyen inférieur, sont examinés à Ecoteaux (VD); indéniablement le site le plus ancien de cette étude. La "Formation inférieure d'Ecoteaux", dont le maigre contenu pollinique traduit un environnement de type désertique froid, présente une rémanence magnétique inverse. Un âge antérieur à 780'000 ans BP (période de Matuyama ou antérieure) est proposé. La "Formation supérieure d'Ecoteaux" comprend plus de huit fluctuations majeures du couvert végétal, parmi lesquelles on compte au moins quatre épisodes tempérés et un optimum climatique. La présence de genres botaniques reliques du Tertiaire, en particulier Pterocarya et Carya, associée à la dynamique végétale faite de courtes oscillations, concourent à lier cette formation au Complexe Cromérien. Les thermomères enregistrées pourraient correspondre à un- ou plusieurs interglaciaires du Pléistocène moyen. Les corrélations possibles avec les stades isotopiques marins, MIS 15, 17 et 19 sont discutées. - Le Pléistocène Moyen et Récent est abordé à travers les sédiments d'Onnens (VD), du Creux d'Enfer (FR), de Port-Valais (VS), de La Dénériaz (VD) et de Cortaillod (NE). Chaque site livre une flore tempérée interglaciaire, dont le contenu pollinique est évalué en fonction des différentes corrélations bio-stratigraphiques envisageables. - Enfin le Würm est étudié grâce à l'affleurement de Marly (FR), accompagné des forages de Villars-sous-Yens (VD), Versoix (GE) et Ollon (VD). L'ensemble de ces observations tend à montrer que l'ouest du Plateau suisse est une zone charnière entre la végétation méridionale de la basse vallée du Rhône et celle d'au-delà du Jura. Au débouché du glacier du Rhône et aux confins des glaciers préalpins et jurassiens, cette région est située au carrefour des migrations botaniques post-glaciaires, entre l'axe majeur de la vallée du Rhône et les zones refuges de l'Est. Cette situation lui confère une importance toute particulière. Alors qu'on pensait le Plateau romand profondément raboté par les glaces, on constate aujourd'hui qu'il recèle plus d'une vingtaine de sites interglaciaires et interstadiaires. L'action abrasive du glacier rhodanien a donc été plus modeste qu'envisagée jusqu'ici. L'extension effective des langues glaciaires, tout comme la circulation atmosphérique, diffèrent sans doute d'un interglaciaire à l'autre. Enfin la vitesse des mouvements glaciaires et l'incroyable vitalité végétale sont des paramètres importants, susceptibles d'avoir entraîné des dépôts, très précocement enrichis en pollens.

Swain: Stata module to correct the SEM chi-square overidentification test in small sample sizes or complex models. Statistical Software Components S457617, Boston College Department of Economics.

Swain corrects the chi-square overidentification test (i.e., likelihood ratio test of fit) for structural equation models whethr with or without latent variables. The chi-square statistic is asymptotically correct; however, it does not behave as expected in small samples and/or when the model is complex (cf. Herzog, Boomsma, & Reinecke, 2007). Thus, particularly in situations where the ratio of sample size (n) to the number of parameters estimated (p) is relatively small (i.e., the p to n ratio is large), the chi-square test will tend to overreject correctly specified models. To obtain a closer approximation to the distribution of the chi-square statistic, Swain (1975) developed a correction; this scaling factor, which converges to 1 asymptotically, is multiplied with the chi-square statistic. The correction better approximates the chi-square distribution resulting in more appropriate Type 1 reject error rates (see Herzog & Boomsma, 2009; Herzog, et al., 2007).

Prediction of infant drug exposure through breastfeeding: population PK modeling and simulation of fluoxetine

BACKGROUND: Risks of significant infant drug exposurethrough breastmilk are poorly defined for many drugs, and largescalepopulation data are lacking. We used population pharmacokinetics(PK) modeling to predict fluoxetine exposure levels ofinfants via mother's milk in a simulated population of 1000 motherinfantpairs.METHODS: Using our original data on fluoxetine PK of 25breastfeeding women, a population PK model was developed withNONMEM and parameters, including milk concentrations, wereestimated. An exponential distribution model was used to account forindividual variation. Simulation random and distribution-constrainedassignment of doses, dosing time, feeding intervals and milk volumewas conducted to generate 1000 mother-infant pairs with characteristicssuch as the steady-state serum concentrations (Css) and infantdose relative to the maternal weight-adjusted dose (relative infantdose: RID). Full bioavailability and a conservative point estimate of1-month-old infant CYP2D6 activity to be 20% of the adult value(adjusted by weigth) according to a recent study, were assumed forinfant Css calculations.RESULTS: A linear 2-compartment model was selected as thebest model. Derived parameters, including milk-to-plasma ratios(mean: 0.66; SD: 0.34; range, 0 - 1.1) were consistent with the valuesreported in the literature. The estimated RID was below 10% in >95%of infants. The model predicted median infant-mother Css ratio was0.096 (range 0.035 - 0.25); literature reported mean was 0.07 (range0-0.59). Moreover, the predicted incidence of infant-mother Css ratioof >0.2 was less than 1%.CONCLUSION: Our in silico model prediction is consistent withclinical observations, suggesting that substantial systemic fluoxetineexposure in infants through human milk is rare, but further analysisshould include active metabolites. Our approach may be valid forother drugs. [supported by CIHR and Swiss National Science Foundation(SNSF)]

The contribution of patch topology and demographic parameters to PVA predictions: the case of the European tree frog

Population viability analyses (PVA) are increasingly used in metapopulation conservation plans. Two major types of models are commonly used to assess vulnerability and to rank management options: population-based stochastic simulation models (PSM such as RAMAS or VORTEX) and stochastic patch occupancy models (SPOM). While the first set of models relies on explicit intrapatch dynamics and interpatch dispersal to predict population levels in space and time, the latter is based on spatially explicit metapopulation theory where the probability of patch occupation is predicted given the patch area and isolation (patch topology). We applied both approaches to a European tree frog (Hyla arborea) metapopulation in western Switzerland in order to evaluate the concordances of both models and their applications to conservation. Although some quantitative discrepancies appeared in terms of network occupancy and equilibrium population size, the two approaches were largely concordant regarding the ranking of patch values and sensitivities to parameters, which is encouraging given the differences in the underlying paradigms and input data.

Functional MRI evaluation of liver tumour response after radiofrequency: short- and mid-term evolution of diffusion parameters

Purpose: To evaluate the short- and mid-term evolutions of the apparent diffusion coefficient of lesions treated with RF, in order to determine if the ADC can be used as a marker of tumour response. Methods and Materials: Twenty patients were treated for a liver malignancy with RF and were examined on a 1.5 T/3.0 T machine with T2, gadolinium-enhanced T1 and diffusion sequences: before treatment (< 1 month), just after treatment (< 1 month) and midterm (3-6 months). The ADC was measured in the whole lesion and in the area with the most restricted diffusion (MRDA). The ROI size was also measured on the diffusion map. The Pearson/ANOVA tests were used. Results: All patients were successfully treated with complete disappearance of CE. The lesional size on T2 showed a negative evolution in time (p < 0.002). The ADC in the whole lesion showed a bell-shaped evolution (increasing just after RF, then decreasing, p = 0.02). The ROI size on the diffusion map followed a similar course (p = 0.01). For the MRDA, such evolutions were also found, but they were not significant. There was a negative correlation between CE and the ADC (p < 0.02) and between the lesional size on T2 and ADC (p = 0.03) in the whole lesion. There were also positive correlations between the ROI size and ADC (p = 0.0008) and between CE and the size on T2 (p = 0.0002). The ADC in MRDA showed some non-significant correlations with other variables. Conclusion: The lesions successfully treated with RF have a clear and predictable evolution in terms of T2 size, CE and ADC.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Aerodynamic drag modeling of alpine skiers performing giant slalom turns.

PURPOSE: Aerodynamic drag plays an important role in performance for athletes practicing sports that involve high-velocity motions. In giant slalom, the skier is continuously changing his/her body posture, and this affects the energy dissipated in aerodynamic drag. It is therefore important to quantify this energy to understand the dynamic behavior of the skier. The aims of this study were to model the aerodynamic drag of alpine skiers in giant slalom simulated conditions and to apply these models in a field experiment to estimate energy dissipated through aerodynamic drag. METHODS: The aerodynamic characteristics of 15 recreational male and female skiers were measured in a wind tunnel while holding nine different skiing-specific postures. The drag and the frontal area were recorded simultaneously for each posture. Four generalized and two individualized models of the drag coefficient were built, using different sets of parameters. These models were subsequently applied in a field study designed to compare the aerodynamic energy losses between a dynamic and a compact skiing technique. RESULTS: The generalized models estimated aerodynamic drag with an accuracy of between 11.00% and 14.28%, and the individualized models estimated aerodynamic drag with an accuracy between 4.52% and 5.30%. The individualized model used for the field study showed that using a dynamic technique led to 10% more aerodynamic drag energy loss than using a compact technique. DISCUSSION: The individualized models were capable of discriminating different techniques performed by advanced skiers and seemed more accurate than the generalized models. The models presented here offer a simple yet accurate method to estimate the aerodynamic drag acting upon alpine skiers while rapidly moving through the range of positions typical to turning technique.

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.