Large Time-Varying Parameter VARs

In this paper we develop methods for estimation and forecasting in large timevarying parameter vector autoregressive models (TVP-VARs). To overcome computational constraints with likelihood-based estimation of large systems, we rely on Kalman filter estimation with forgetting factors. We also draw on ideas from the dynamic model averaging literature and extend the TVP-VAR so that its dimension can change over time. A final extension lies in the development of a new method for estimating, in a time-varying manner, the parameter(s) of the shrinkage priors commonly-used with large VARs. These extensions are operationalized through the use of forgetting factor methods and are, thus, computationally simple. An empirical application involving forecasting inflation, real output, and interest rates demonstrates the feasibility and usefulness of our approach.

Large Bayesian VARMAs

Vector Autoregressive Moving Average (VARMA) models have many theoretical properties which should make them popular among empirical macroeconomists. However, they are rarely used in practice due to over-parameterization concerns, difficulties in ensuring identification and computational challenges. With the growing interest in multivariate time series models of high dimension, these problems with VARMAs become even more acute, accounting for the dominance of VARs in this field. In this paper, we develop a Bayesian approach for inference in VARMAs which surmounts these problems. It jointly ensures identification and parsimony in the context of an efficient Markov chain Monte Carlo (MCMC) algorithm. We use this approach in a macroeconomic application involving up to twelve dependent variables. We find our algorithm to work successfully and provide insights beyond those provided by VARs.

Secondary contact zones and hybridizations: the case of the lesser white-toothed shrew (Crocidura suaveolens group, Soricidae)

In the present study, we analyzed 58 samples of the lesser white-toothed shrew group (Crocidura suaveolens) from eastern Europe and Turkey, where, according to previous publications, three different mitochondrial and nuclear lineages are present. We sequenced 799 bp of the nuclear BRCA1 gene and 400 bp of the mitochondrial cytochrome b gene to: (1) determine a potential contact zone between the lineages; (2) detect hybridizations and introgressions between them; and (3) comment on the level of reproductive isolation of the different lineages. We revealed two zones of hybridization in Turkey, of which the first occurred west of the Bosphorus Straits (three hybrids) and the second in Anatolia (twelve hybrids). In the latter, the nuclear markers revealed a large zone of hybridization, of approximately 600 km. It also revealed that hybrids of first, second, and later generations are present within the populations, and therefore that the reproductive isolation between the different lineages is weak.

Species identification in mammals from mixed biological samples based on mitochondrial DNA control region length polymorphism.

The ability to identify the species origin of an unknown biological sample is relevant in the fields of human and wildlife forensics. However, the detection of several species mixed in the same sample still remains a challenge. We developed and tested a new approach for mammal DNA identification in mixtures of two or three species, based on the analysis of mitochondrial DNA control region interspecific length polymorphism followed by direct sequencing. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and is not based on a pre-defined panel of species. Amplicons can be separated either on agarose gels or using CE. The advantages and limitations of the assay are discussed under different conditions, such as variable template concentration, amplicon sizes and size difference among the amplicons present in the mixture. For the first time, this protocol provides a simple, reliable and flexible method for simultaneous identification of multiple mammalian species from mixtures, without any prior knowledge of the species involved.

Description of Leishmania equatorensis sp. n. (Kinetoplastida: Trypanosomatidae), a new parasite infecting arboreal mammals in Ecuador

Characterization is given of a new parasite, Leishmania equatoriensis sp.n. wich was isolated from the viscera of a sloth (Choloepus hoffmanni) and a squirrel (Sciurus granatensis), captured in humid tropical forest onthe Pacific Coast of Ecuador. Data based on biological and molecular criteria, as well as numerical zymotaxonomical analysis, indicate that this parasite is a new species of the L. brasiliensis complex. L. equatoriensis is cleary distinguishable form all other known species within this complex, using the following molecular criteria: reactivity patterns with specific monoclonal antibodies, isoenzyme electrophoresis, and restriction-endonuclease fragment patterns of kinetoplast DNA (k-DNA).

Ectoparasites Ixodida Leach, 1817 on wild mammals in the state of Paraná, Brazil

A taxonomical and ecological study was made on 264 samples of 12 species of ectoparasite ixodides collected on wild mammals from several natural regions of Paraná state, Brazil. These species of ticks as their hosts are listed by their identified evolutive forms and capture locations. A review of the early studies on the Ixodida from South Brazil is made considering the frequency of the parasite species on each host concerning the geographic distribution of such species. This paper is a contribution to the knowledge of the Ixodidae ectoparasites species of wild mammals in Brazil.

Large-scale imatinib dose-concentration-effect study in CML patients under routine care conditions.

This observational study analyzed imatinib pharmacokinetics and response in 2478 chronic myeloid leukemia (CML) patients. Data were obtained through centralized therapeutic drug monitoring (TDM) at median treatment duration of ≥2 years. First, individual initial trough concentrations under 400mg/day imatinib starting dose were estimated. Second, their correlation (C^min(400mg)) with reported treatment response was verified. Low imatinib levels were predicted in young male patients and those receiving P-gp/CYP3A4 inducers. These patients had also lower response rates (7% lower 18-months MMR in male, 17% lower 1-year CCyR in young patients, Kaplan-Meier estimates). Time-point independent multivariate regression confirmed a correlation of individual C^min(400mg) with response and adverse events. Possibly due to confounding factors (e.g. dose modifications, patient selection bias), the relationship seemed however flatter than previously reported from prospective controlled studies. Nonetheless, these observational results strongly suggest that a subgroup of patients could benefit from early dosage optimization assisted by TDM, because of lower imatinib concentrations and lower response rates.

The current molecular phylogeny of Eutherian mammals challenges previous interpretations of placental evolution.

Based on histology, the placentae of eutherians are currently grouped in epitheliochorial, endotheliochorial and haemochorial placentae. In a haeckelian sense, the epitheliochorial contact with marked histiotrophic feeding by uterine milk is generally considered as primitive, especially since similar contacts exist in Marsupials. In contrast, the more intimate endotheliochorial and haemochorial contact, facilitating haemotrophic nutrition, is interpreted as a derived state. A cladistic analysis based on the phylogenetic relationships established by molecular analyses reveals that the basic clades are all characterized by an endotheliochorial or haemochorial placenta, and that the epitheliochorial placenta evolved at least three times in a convergent manner. This evolution may be explained by the fact that the epitheliochorial placenta in eutherians is more efficient in nutritional transfer (flow rate by exchange surface). Moreover, this arrangement may confer an advantage to the mother who can probably reduce the degree of manipulation by a genetically imprinted embryo.

Immunoglobulin classes in aquatic mammals. Characterization by serologic cross-reactivity, molecular size and binding of human free secretory component.

On the basis of serologic cross-reactivity, three immunoglobulin classes homologous to human IgG, IgM and IgA were identified in two species of acquatic mammal representing the orders Cetacea (dolphin) and Pinnipedea (sea lion). Molecular size was estimated by sucrose density gradient ultracentrifugation and Sephadex G-200 chromatography, indicating a 7S IgG, 19S IgM and heterogeneous serum IgA. Human secretory component was readily bound to the IgM of both species and to an apparently lesser extent to the larger molecular size populations of IgA. No binding was observed with IgG. Several antisera specific for human γ-chains gave a single precipitin line with the sea lion IgG but when made to react with dolphin serum produced two lines, suggesting the presence of two different subclasses of IgG in this species.

Study of humoral immune response in mammals immunized with Plasmodium falciparum antigenic preparations

Six Plasmodium falciparum protein fractions, isolated under reducing conditions, were used to immunize mice, rabbits and the squirrel monkey Saimiri sciureus. Five or seven subcutaneous injections of each antigenic preparation, in conjunction with Freund's complete or incomplete adjuvants, were administered. This led to the development of specific antibodies detected by IFAT, ELISA or immunobloting which inhibited merozoite reinvasion in in vitro P. falciparum cultures. This activity seems to be associated with rhoptry proteins contained in fractions Pf F2 and Pf F4.

Superinfection with drug-resistant HIV is rare and does not contribute substantially to therapy failure in a large European cohort.

BACKGROUND: Superinfection with drug resistant HIV strains could potentially contribute to compromised therapy in patients initially infected with drug-sensitive virus and receiving antiretroviral therapy. To investigate the importance of this potential route to drug resistance, we developed a bioinformatics pipeline to detect superinfection from routinely collected genotyping data, and assessed whether superinfection contributed to increased drug resistance in a large European cohort of viremic, drug treated patients. METHODS: We used sequence data from routine genotypic tests spanning the protease and partial reverse transcriptase regions in the Virolab and EuResist databases that collated data from five European countries. Superinfection was indicated when sequences of a patient failed to cluster together in phylogenetic trees constructed with selected sets of control sequences. A subset of the indicated cases was validated by re-sequencing pol and env regions from the original samples. RESULTS: 4425 patients had at least two sequences in the database, with a total of 13816 distinct sequence entries (of which 86% belonged to subtype B). We identified 107 patients with phylogenetic evidence for superinfection. In 14 of these cases, we analyzed newly amplified sequences from the original samples for validation purposes: only 2 cases were verified as superinfections in the repeated analyses, the other 12 cases turned out to involve sample or sequence misidentification. Resistance to drugs used at the time of strain replacement did not change in these two patients. A third case could not be validated by re-sequencing, but was supported as superinfection by an intermediate sequence with high degenerate base pair count within the time frame of strain switching. Drug resistance increased in this single patient. CONCLUSIONS: Routine genotyping data are informative for the detection of HIV superinfection; however, most cases of non-monophyletic clustering in patient phylogenies arise from sample or sequence mix-up rather than from superinfection, which emphasizes the importance of validation. Non-transient superinfection was rare in our mainly treatment experienced cohort, and we found a single case of possible transmitted drug resistance by this route. We therefore conclude that in our large cohort, superinfection with drug resistant HIV did not compromise the efficiency of antiretroviral treatment.

Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

Large multi-gene phylogenetic trees of the grasses (Poaceae): progress towards complete tribal and generic level sampling.

In this paper we included a very broad representation of grass family diversity (84% of tribes and 42% of genera). Phylogenetic inference was based on three plastid DNA regions rbcL, matK and trnL-F, using maximum parsimony and Bayesian methods. Our results resolved most of the subfamily relationships within the major clades (BEP and PACCMAD), which had previously been unclear, such as, among others the: (i) BEP and PACCMAD sister relationship, (ii) composition of clades and the sister-relationship of Ehrhartoideae and Bambusoideae + Pooideae, (iii) paraphyly of tribe Bambuseae, (iv) position of Gynerium as sister to Panicoideae, (v) phylogenetic position of Micrairoideae. With the presence of a relatively large amount of missing data, we were able to increase taxon sampling substantially in our analyses from 107 to 295 taxa. However, bootstrap support and to a lesser extent Bayesian inference posterior probabilities were generally lower in analyses involving missing data than those not including them. We produced a fully resolved phylogenetic summary tree for the grass family at subfamily level and indicated the most likely relationships of all included tribes in our analysis.