On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Distinguishing species and populations of Rhipicephaline ticks with its 2 ribosomal RNA

Most populations and some species of ticks of the genera Boophilus (5 spp.) and Rhipicephalus (ca. 75 spp.) cannot be distinguished phenotypically. Moreover, there is doubt about the validity of species in these genera. I studied the entire second internal transcribed spacer (ITS 2) rRNA of 16 populations of rhipicephaline ticks to address these problems: Boophilus,microplus from Australia, Kenya, South Africa and Brazil (4 populations); Boophilus decoloratus from Kenya; Rhipicephalus appendiculatus from Kenya, Zimbabwe and Zambia (7 populations); Rhipicephalus zambesiensis from Zimbabwe (3 populations); and Rhipicephalus evertsi from Kenya. Each of the 16 populations had a unique ITS 2, but most of the nucleotide variation occurred among species and genera. ITS 2 rRNA can be used to distinguish the populations and species of Boophilus and Rhipicephalus studied here. Little support was found for the hypothesis that B. microplus from Australia and South Africa are different species. ITS 2 appears useful for phylogenetic inference in the Rhipicephalinae because in genetic distance, maximum likelihood, and maximum parsimony analyses, most branches leading to species had >95% bootstrap support. Rhipicephalus appendiculatus and R, zambeziensis are closely related, yet their ITS 2 sequences could be distinguished unambiguously. This lends weight to a previous proposal that Rhipicephalus sanguineus and Rhipicephalus turanicus, and Rhipicephalus pumlilio and Rhipicephalus camicasi, respectively, are conspecific, because each of these pairs of species had identical sequences for ca. 250 bp of ITS 2 rRNA.

V4 region of small subunit rDNA indicates polyphyly of the Fellodistomidae (Digenea) which is supported by morphology and life-cycle data

There is no morphological synapomorphy for the disparate digeneans, the Fellodistomidae Nicoll, 1909. Although all known life-cycles of the group include bivalves as first intermediate hosts, there is no convincing morphological synapomorphy that can be used to unite the group. Sequences from the V4 region of small subunit (18S) rRNA genes were used to infer phylogenetic relationships among 13 species of Fellodistomidae from four subfamilies and eight species from seven other digenean families: Bivesiculidae; Brachylaimidae; Bucephalidae; Gorgoderidae; Gymnophallidae; Opecoelidae; and Zoogonidae. Outgroup comparison was made initially with an aspidogastrean. Various species from the other digenean families were used as outgroups in subsequent analyses. Three methods of analysis indicated polyphyly of the Fellodistomidae and at least two independent radiations of the subfamilies, such that they were more closely associated with other digeneans than to each other. The Tandanicolinae was monophyletic (100% bootstrap support) and was weakly associated with the Gymnophallidae (< 50-55% bootstrap support). Monophyly of the Baccigerinae was supported with 78-87% bootstrap support, and monophyly of the Zoogonidae + Baccigerinae received 77-86% support. The remaining fellodistomid species, Fellodistomum fellis, F. agnotum and Coomera brayi (Fellodistominae) plus Proctoeces maculatus and Complexobursa sp. (Proctoecinae), formed a separate clade with 74-92% bootstrap support. On the basis of molecular, morphological and life-cycle evidence, the subfamilies Baccigerinae and Tandanicolinae are removed from the Fellodistomidae and promoted to familial status. The Baccigerinae is promoted under the senior synonym Faustulidae Poche, 1926, and the Echinobrevicecinae Dronen, Blend & McEachran, 1994 is synonymised with the Faustulidae. Consequently, species that were formerly in the Fellodistomidae are now distributed in three families: Fellodistomidae; Faustulidae (syn. Baccigerinae Yamaguti, 1954); and Tandanicolidae Johnston, 1927. We infer that the use of bivalves as intermediate hosts by this broad range of families indicates multiple host-switching events within the radiation of the Digenea.

Hemichordates and deuterostome evolution: robust molecular phylogenetic support for a hemichordate plus echinoderm clade

Hemichordates were traditionally allied to the chordates, but recent molecular analyses have suggested that hemichordates are a sister group to the echinoderms, a relationship that has important consequences for the interpretation of the evolution of deuterostome body plans. However, the molecular phylogenetic analyses to date have not provided robust support for the hemichordate + echinoderm clade. We use a maximum likelihood framework, including the parametric bootstrap, to reanalyze DNA data from complete mitochondrial genomes and nuclear 18S rRNA. This approach provides the first statistically significant support for the hemichordate + echinoderm clade from molecular data. This grouping implies that the ancestral deuterostome had features that included an adult with a pharynx and a dorsal nerve cord and an indirectly developing dipleurula-like larva.

Modelling life history strategies with capture-recapture data: Evolutionary demography of the water skink Eulamprus tympanum

Matrix population models, elasticity analysis and loop analysis can potentially provide powerful techniques for the analysis of life histories. Data from a capture-recapture study on a population of southern highland water skinks (Eulamprus tympanum) were used to construct a matrix population model. Errors in elasticities were calculated by using the parametric bootstrap technique. Elasticity and loop analyses were then conducted to identify the life history stages most important to fitness. The same techniques were used to investigate the relative importance of fast versus slow growth, and rapid versus delayed reproduction. Mature water skinks were long-lived, but there was high immature mortality. The most sensitive life history stage was the subadult stage. It is suggested that life history evolution in E. tympanum may be strongly affected by predation, particularly by birds. Because our population declined over the study, slow growth and delayed reproduction were the optimal life history strategies over this period. Although the techniques of evolutionary demography provide a powerful approach for the analysis of life histories, there are formidable logistical obstacles in gathering enough high-quality data for robust estimates of the critical parameters.

Microcystin production by a freshwater spring cyanobacterium of the genus Fischerella

We investigated the production of a hepatotoxic, cyclic heptapeptide, microcystin, by a filamentous branched cyanobacterium belonging to the order Stigonematales, genus Fischerella. The freshwater Fischerella sp. strain CENA161 was isolated from spring water in a small concrete dam in Piracicaba, Sao Paulo State, Brazil, and identified by combining a morphological description with 16S rRNA gene sequencing and phylogenetic analysis. Microcystin (MCYST) analysis performed using an ELISA assay on cultured cells gave positive results. High performance liquid chromatography-mass spectrometry (HPLC-MS) analysis detected 33.6 mu g MCYST-LR per gram dry weight of cyanobacterial cells. Microcystin profile revealed by quadrupole time-of-flight tandem mass spectrometry (Q-TOF-MS/MS) analysis confirmed the production of MCYST-LR. Furthermore, genomic DNA was analyzed by PCR for sequences similar to the ketosynthase (KS) domain of the type I polyketide synthase gene, which is involved in microcystin biosynthesis. This revealed the presence of a KS nucleotide fragment similar to the mcyD and ndaD genes of the microcystin and nodularin synthetase complexes. Phylogenetic analysis grouped the Fischerella KS sequence together with mcyD sequences of the three known microcystin synthetase operon (Microcystis, Planktothrix and Anabaena) and ndaD of the nodularin synthetase operon, with 100% bootstrap support. Our findings demonstrate that Fischerella sp. CENA161 produces MYCST-LR and for the first time identify a nucleotide sequence putatively involved in microcystin synthesis in this genus. (C) 2009 Elsevier Ltd. All rights reserved.

Ancestry informative markers in Amerindians from Brazilian Amazon

Ancestry informative markers (AIMs) are genetic loci with large frequency differences between the major ethnic groups and are very useful in admixture estimation. However, their frequencies are poorly known within South American indigenous populations, making it difficult to use them in admixture studies with Latin American populations, such as the trihybrid Brazilian population. To minimize this problem, the frequencies of the AIMs FY-null RB2300, LPL, AT3-1/1), Sb19.3, APO, and PV92 were determined via PCR and PCR-RFLP in four tribes from Brazilian Amazon (Tikuna, Kashinawa, Baniwa, and Kanamari), to evaluate their potential for discriminating indigenous populations from Europeans and Africans, as well as discriminating each tribe from the others. Although capable of differentiating tribes, as evidenced by the exact test of population differentiation, a neighbor-joining tree suggests that the AIMs are useless in obtaining reliable reconstructions of the biological relationships and evolutionary history that characterize the villages and tribes studied. The mean allele frequencies from these AIMs were very similar to those observed for North American natives. They discriminated Amerindians from Africans, but not from Europeans. On the other hand, the neighbor-joining dendrogram separated Africans and Europeans from Amerindians with a high statistical support (bootstrap = 0.989). The relatively low diversity (GST = 0.042) among North American natives and Amerindians from Brazilian Amazon agrees with the lack of intra-ethnic variation previously reported for these markers. Despite genetic drift effects, the mean allelic frequencies herein presented could be used as Amerindian parental frequencies in admixture estimates in urban Brazilian populations.

Molecular characterization of potentially zoonotic isolates of Giardia duodenalis in horses

Giardia isolates from eight horses from New York State (NY), USA and two horses from Western Australia (WA) were genetically characterized at the SSU-rDNA and triose-phosphate isomerase (TPI) genes. Phylogenetic analysis of the TPI gene provided strong support for the placement of both isolates of Giardia from horses in WA and a single isolate from a horse in NY within the assemblage AI genotype of G. duodenalis. Another two isolates from horses in NY placed within the assemblage All genotype of G. duodenalis. Phylogenetic analysis of the TPI gene also provided strong bootstrap support for the placement of four G. duodenalis isolates from horses in NY into a potentially host-specific sub-assemblage of assemblage BIV. The results of this study are consistent with previous studies showing that assemblages AI and AII of G. duodenalis provide the greatest potential zoonotic risk to humans. Horses may therefore constitute a potential source for human infection of Giardia either directly or via watersheds. (c) 2005 Elsevier B.V. All rights reserved.

Standard errors of fitted component means of normal mixtures

In this paper use consider the problem of providing standard errors of the component means in normal mixture models fitted to univariate or multivariate data by maximum likelihood via the EM algorithm. Two methods of estimation of the standard errors are considered: the standard information-based method and the computationally-intensive bootstrap method. They are compared empirically by their application to three real data sets and by a small-scale Monte Carlo experiment.

Molecular evidence of horizontal transmission of hepatitis C virus within couples

Hepatitis C virus (HCV) transmission has decreased with the adoption of universal blood donor screening and social policies to reduce the risk of infection in intravenous drug users, but remains a worldwide health problem. The objective of this study was to evaluate the phylogenetic relationships among sequences from different HCV genomic regions from sexual partners of infected patients. Nine couples with a stable relationship and without other risk factors for HCV infection and 42 control patients were selected, and the NS3 and NS5B regions were analysed. Phylogenetic analysis showed that viruses from five of the couples had a common origin, clustering in the same monophyletic group, with bootstrap values greater than 70. For the other couples, monophyletic groups were observed, but without bootstrap support. Thus, using two different viral genome regions, a common source of infection was observed in both members of five couples. These data strongly support HCV transmission within couples.

Carios mimon (Acari: Argasidae): description of adults and redescription of larva

Carios mimon is an argasid tick common on Chiroptera, originally described from larvae collected on bats Mimon crenulatum from Bolivia and Eptesicus brasiliensis from Uruguay. Later it was also registered from Argentina and recently included among the Brazilian tick fauna. In Brazil, this species is very aggressive to man, resulting in intense inflammatory response and pain. It is known only by the larval description and its morphology resembles that from other species currently included into the genus Carios, formerly classified into the subgenus Alectorobius, genus Ornithodoros. Here we describe adults and redescribe the larva of C. mimon, based on light and scanning electron microscopy. Remarks about its morphological similarity with other species of this genus are also discussed. Molecular analysis inferred from a portion of the 16S rRNA mitochondrial gene placed C. mimon in a cluster supported by maximal bootstrap value (100%) with other argasid species (mostly bat parasites in the New World), which have been classified into either the genus Ornithodoros or Carios, depending on the Argasidae classification adopted by different authors.

Detection and molecular characterization of a canine piroplasm from Brazil

In the beginning of the 20th century, a new canine disease was reported in Brazil under the name ""nambiuvu"", whose etiological agent was called Rangelia vitalii, a distinct piroplasm that was shown to parasitize not only erythrocytes, but also leucocytes and endothelial cells. In this new century, more publications on R. vitalii were reported from Brazil, including an extensive study on its ultrastructural analysis, in addition to clinical, pathological, and epidemiological data on nambiuvu. However, a molecular analysis of R. vitalii has not been performed to date. In the present study, we performed molecular phylogenetic analyses of R. vitalii based on fragments of the genes 18S rRNA and the heat shock protein 70 (hsp70), amplified by PCR performed on blood samples derived from five clinical cases of dogs presumably infected with R. vitalii in southern Brazil. In addition, we examined Giemsa-stained thin blood smears from these same dogs. DNA sequences (604-bp) of the 18S rRNA gene obtained from the five dogs were identical to each other, and by Blast analysis, this sequence shared the highest degree of sequence identity (95%) with Babesia sp. China-BQ1. DNA sequences (1056-bp) of the hsp70 gene obtained from the five dogs were identical to each other, and by Blast analysis, this sequence shared the highest degree of sequence identity (87%) with Babesia bigemina. Phylogenetic analyses inferred from either of the two genes resulted in the newly genotype being placed in the Babesia spp.sensu stricto clade with very high bootstrap support (95-100%) in three analyses (Neighbor-Joining, Maximum parsimony, and Maximum likelihood). Giemsa-stained thin blood smears from the dogs were shown to contain piroplasm organisms within erythrocytes, monocytes and neutrophils (individual forms), and schizont-like forms within neutrophils, in accordance with literature reports of R. vitalii. Based on these results, we conclude that R. vitalii, the etiological agent of ""nambiuvu"" in southern Brazil, is a valid species of piroplasm. Further studies are required to evaluate the validity of the genus Rangelia. (C) 2011 Elsevier B.V. All rights reserved.

Phylogenetic analysis of the Monocotylidae (Monogenea) inferred from 28S rDNA sequences

The current classification of the Monocotylidae (Monogenea) is based on a phylogeny generated from morphological characters. The present study tests the morphological phylogenetic hypothesis using molecular methods. Sequences from domains C2 and D1 and the partial domains C1 and D2 from the 28S rDNA gene for 26 species of monocotylids from six of the seven subfamilies were used. Trees were generated using maximum parsimony, neighbour joining and maximum likelihood algorithms. The maximum parsimony tree, with branches showing less than 70% bootstrap support collapsed, had a topology identical to that obtained using the maximum likelihood analysis. The neighbour joining tree, with branches showing less than 70% support collapsed. differed only in its placement of Heterocotyle capricornensis as the sister group to the Decacotylinae clade. The molecular tree largely supports the subfamilies established using morphological characters. Differences are primarily how the subfamilies are related to each other. The monophyly of the Calicotylinae and Merizocotylinae and their sister group relationship is supported by high bootstrap values in all three methods, but relationships within the Merizocotylinae are unclear. Merizocotyle is paraphyletic and our data suggest that Mycteronastes and Thaumatocotyle, which were synonymized with Merizocotyle after the morphological cladistic analysis, should perhaps be resurrected as valid genera. The monophyly of the Monocotylinae and Decacotylinae is also supported by high bootstrap values. The Decacotylinae, which was considered previously to be the sister group to the Calicotylinae plus Merizocotylinae, is grouped in an unresolved polychotomy with the Monocotylinae and members of the Heterocotylinae. According to our molecular data, the Heterocotylinae is paraphyletic. Molecular data support a sister group relationship between Troglocephalus rhinobatidis and Neoheterocotyle rhinobatidis to the exclusion of the other species of Neoheterocotyle and recognition of Troglocephalus renders Neoheterocotyle,le paraphyletic. We propose Troglocephalus incertae sedis. An updated classification and full species list of the Monocotylidae is provided. (C) 2001 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.