Bothrops insularis venomics: A proteomic analysis supported by transcriptomic-generated sequence data

A joint transcriptomic and proteomic approach employing two-dimensional electrophoresis, liquid chromatography and mass spectrometry was carried out to identify peptides and proteins expressed by the venom gland of the snake Bothrops insularis, an endemic species of Queimada Grande Island, Brazil. Four protein families were mainly represented in processed spots, namely metalloproteinase, serine proteinase, phospholipase A(2) and lectin. Other represented families were growth factors, the developmental protein G10, a disintegrin and putative novel bradykinin-potentiating peptides. The enzymes were present in several isoforms. Most of the experimental data agreed with predicted values for isoelectric point and M(r) of proteins found in the transcriptome of the venom gland. The results also support the existence of posttranslational modifications and of proteolytic processing of precursor molecules which could lead to diverse multifunctional proteins. This study provides a preliminary reference map for proteins and peptides present in Bothrops insularis whole venom establishing the basis for comparative studies of other venom proteomes which could help the search for new drugs and the improvement of venom therapeutics. Altogether, our data point to the influence of transcriptional and post-translational events on the final venom composition and stress the need for a multivariate approach to snake venomics studies. (c) 2009 Elsevier B.V. All rights reserved.

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

The R2 mobile element of Rhynchosciara americana: Molecular, cytological and dynamic aspects

Ribosomal RNA genes are encoded by large units clustered (18S, 5S, and 28S) in the nucleolar organizer region in several organisms. Sometimes additional insertions are present in the coding region for the 28S rDNA. These insertions are specific non-long terminal repeat retrotransposons that have very restricted integration targets within the genome. The retrotransposon present in the genome of Rhynchosciara americana, RaR2, was isolated by the screening of a genomic library. Sequence analysis showed the presence of conserved regions, such as a reverse transcriptase domain and a zinc finger motif in the amino terminal region. The insertion site was highly conserved in R. americana and a phylogenetic analysis showed that this element belongs to the R2 clade. The chromosomal localization confirmed that the RaR2 mobile element was inserted into a specific site in the rDNA gene. The expression level of RaR2 in salivary glands during larval development was determined by quantitative RT-PCR, and the increase of relative expression in the 3P of the fourth instar larval could be related to intense gene activity characteristic of this stage. 5`-Truncated elements were identified in different DNA samples. Additionally, in three other Rhynchosciara species, the R2 element was present as a full-length element.

Structural and phylogenetic relationship of ORF 31 from the Anticarsia gemmatalis MNPV to poly (ADP-ribose) polymerases (PARP)

ORF 31 is a unique baculovirus gene in the genome of Anticarsia gemmatalis multiple nucleopolyhedrovirus isolate 2D (AgMNPV-2D). It encodes a putative polypeptide of 369 aa homologous to poly (ADP-ribose) polymerase (PARP) found in the genomes of several organisms. Moreover, we found a phylogenetic association with Group I PARP proteins and a 3D homology model of its conserved PARP C-terminal catalytic domain indicating that had almost an exact spatial superimposition of < 1 angstrom with other PARP available structures. The 5` end of ORF 31 mRNA was located at the first nucleotide of a CATT motif at position -27. Using real-time PCR we detected transcripts at 3 h post-infection (p.i.) increasing until 24 h p.i., which coincides with the onset of DNA replication, suggestive of a possible role in DNA metabolism.

Phylogenetic analysis of a near-full-length sequence of an erythrovirus genotype 3 strain isolated in Brazil

Human parvovirus B19 is the only member of the genus Erythrovirus that causes human disease. Recent findings of several strains with considerable sequence divergence from B19 have suggested a new classification for parvovirus genotypes as 1 (B19), 2 (A-6 and LaLi) and 3 (V9). In their overall DNA sequence, the three genotypes differ by similar to 10%. Here, we report the isolation of a genotype-3-related strain named BR543 during a prospective study conducted in Sao Paulo, Brazil. Analysis of the nearly full-length genome sequence of BR543 indicates that this B19 variant sequence clusters with Gh2768, a strain from Ghana belonging to subtype 3b, and showed mostly synonymous substitutions.

Paracoccidioides brasiliensis Vaccine Formulations Based on the gp43-Derived P10 Sequence and the Salmonella enterica FliC Flagellin

Paracoccidioidomycosis (PCM) is a systemic granulomatous disease caused by the dimorphic fungus Paracoccidioides brasiliensis. Anti-PCM vaccine formulations based on the secreted fungal cell wall protein (gp43) or the derived P10 sequence containing a CD4(+) T-cell-specific epitope have shown promising results. In the present study, we evaluated new anti-PCM vaccine formulations based on the intranasal administration of P. brasiliensis gp43 or the P10 peptide in combination with the Salmonella enterica FliC flagellin, an innate immunity agonist binding specifically to the Toll-like receptor 5, in a murine model. BALB/c mice immunized with gp43 developed high-specific-serum immunoglobulin G1 responses and enhanced interleukin-4 (IL-4) and IL-10 levels. On the other hand, mice immunized with recombinant purified flagellins genetically fused with P10 at the central hypervariable domain, either flanked or not by two lysine residues, or the synthetic P10 peptide admixed with purified FliC elicited a prevailing Th1-type immune response based on lung cell-secreted type 1 cytokines. Mice immunized with gp43 and FliC and intratracheally challenged with P. brasiliensis yeast cells had increased fungal proliferation and lung tissue damage. In contrast, mice immunized with the chimeric flagellins and particularly those immunized with P10 admixed with FliC reduced P. brasiliensis growth and lung damage. Altogether, these results indicate that S. enterica FliC flagellin modulates the immune response to P. brasiliensis P10 antigen and represents a promising alternative for the generation of anti-PCM vaccines.

Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58 Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates. (C) 2009 Elsevier B.V. All rights reserved.

Worldwide sequence conservation of transmission-blocking vaccine candidate Pvs230 in Plasmodium vivax

Pfs230, surface protein of gametocyte/gamete of the human malaria parasite, Plasmodium falciparum, is a prime candidate of malaria transmission-blocking vaccine. Plasmodium vivax has an ortholog of Pfs230 (Pvs230), however, there has been no study in any aspects on Pvs230 to date. To investigate whether Pvs230 can be a vivax malaria transmission-blocking vaccine, we performed evolutionary and population genetic analysis of the Pvs230 gene (pvs230: PVX_003905). Our analysis of Pvs230 and its orthologs in eight Plasmodium species revealed two distinctive parts: an interspecies variable part (IVP) containing species-specific oligopeptide repeats at the N-terminus and a 7.5 kb interspecies conserved part (ICP) containing 14 cysteine-rich domains. Pvs230 was closely related to its orthologs, Pks230 and Pcys230, in monkey malaria parasites. Analysis of 113 pvs230 sequences obtained from worldwide, showed that nucleotide diversity is remarkably low in the non-repeat 8-kb region of pvs230 (theta pi = 0.00118) with 77 polymorphic nucleotide sites, 40 of which results in amino acid replacements. A signature of purifying selection but not of balancing selection was seen on pvs230. Functional and/or structural constraints may limit the level of polymorphism in pvs230. The observed limited polymorphism in pvs230 should ground for utilization of Pvs230 as an effective transmission-blocking vaccine. (C) 2011 Elsevier Ltd. All rights reserved.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7 Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5` and 3` segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition. we uncovered two surprising patterns of evolution in repetitive sequence. Firsts in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi. (C) 2009 Elsevier B.V. All rights reserved.

Sequence-specific reconstruction from fragmentary databases using seed sequences: implementation and validation on SAGE, proteome and generic sequencing data

Motivation: DNA assembly programs classically perform an all-against-all comparison of reads to identify overlaps, followed by a multiple sequence alignment and generation of a consensus sequence. If the aim is to assemble a particular segment, instead of a whole genome or transcriptome, a target-specific assembly is a more sensible approach. GenSeed is a Perl program that implements a seed-driven recursive assembly consisting of cycles comprising a similarity search, read selection and assembly. The iterative process results in a progressive extension of the original seed sequence. GenSeed was tested and validated on many applications, including the reconstruction of nuclear genes or segments, full-length transcripts, and extrachromosomal genomes. The robustness of the method was confirmed through the use of a variety of DNA and protein seeds, including short sequences derived from SAGE and proteome projects.

Intraspecific sequence variation in 16S rRNA gene of Ureaplasma diversum isolates

Ureaplasma diversum infection in bulls may result in seminal vesiculitis, balanoposthitis and alterations in spermatozoids. In cows, it can cause placentitis, fetal alveolitis, abortion and the birth of weak calves. U. diversum ATCC 49782 (serogroups A), ATCC 49783 (serogroup C) and 34 field isolates were used for this study. These microorganisms were submitted to Polymerase Chain Reaction for 16S gene sequence determination using Tact High Fidelity and the products were purified and bi-directionally sequenced. Using the sequence obtained, a fragment containing four hypervariable regions was selected and nucleotide polymorphisms were identified based on their position within the 16S rRNA gene. Forty-four single nucleotide polymorphisms (SNP) were detected. The genotypic variability of the 16S rRNA gene of U. diversum isolates shows that the taxonomy classification of these organisms is likely much more complex than previously described and that 16S rRNA gene sequencing may be used to suggest an epidemiologic pattern of different origin strains. (c) 2011 Elsevier B.V. All rights reserved.

Identification of a novel ligand binding motif in the transthyretin channel

The design of therapeutic compounds targeting transthyretin (TTR) is challenging due to the low specificity of interaction in the hormone binding site. Such feature is highlighted by the interactions of TTR with diclofenac, a compound with high affinity for TTR, in two dissimilar modes, as evidenced by crystal structure of the complex. We report here structural analysis of the interactions of TTR with two small molecules, 1-amino-5-naphthalene sulfonate (1,5-AmNS) and 1-anilino-8-naphthalene sulfonate (1,8-ANS). Crystal structure of TTR: 1,8-ANS complex reveals a peculiar interaction, through the stacking of the naphthalene ring between the side-chain of Lys15 and Leu17. The sulfonate moiety provides additional interaction with Lys15` and a water-mediated hydrogen bond with Thr119`. The uniqueness of this mode of ligand recognition is corroborated by the crystal structure of TTR in complex with the weak analogue 1,5-AmNS, the binding of which is driven mainly by hydrophobic partition and one electrostatic interaction between the sulfonate group and the Lys15. The ligand binding motif unraveled by 1,8-ANS may open new possibilities to treat TTR amyloid diseases by the elucidation of novel candidates for a more specific pharmacophoric pattern. (C) 2009 Published by Elsevier Ltd.

The three-dimensional structure of bothropasin, the main hemorrhagic factor from Bothrops jararaca venom: Insights for a new classification of snake venom metalloprotease subgroups

Bothropasin is a 48 kDa hemorrhagic PIII snake venom metalloprotease (SVMP) isolated from Bothrops jararaca, containing disintegrin/cysteine-rich adhesive domains. Here we present the crystal structure of bothropasin complexed with the inhibitor POL647. The catalytic domain consists of a scaffold of two subdomains organized similarly to those described for other SVMPs, including the zinc and calcium-binding sites. The free cysteine residue Cys(189) is located within a hydrophobic core and it is not available for disulfide bonding or other interactions. There is no identifiable secondary structure for the disintegrin domain, but instead it is composed mostly of loops stabilized by seven disulfide bonds and by two calcium ions. The ECD region is in a loop and is structurally related to the RGD region of RGD disintegrins, which are derived from I`ll SVMPs. The ECD motif is stabilized by the Cys(117)_Cys(310) disulfide bond (between the disintegrin and cysteine-rich domains) and by one calcium ion. The side chain of Glu(276) of the ECD motif is exposed to solvent and free to make interactions. In bothropasin, the HVR (hyper-variable region) described for other Pill SVMPs in the cysteine-rich domain, presents a well-conserved sequence with respect to several other Pill members from different species. We propose that this subset be referred to as PIII-HCR (highly conserved region) SVMPs. The differences in the disintegrin-like, cysteine-rich or disintegrin-like cysteine-rich domains may be involved in selecting target binding, which in turn could generate substrate diversity or specificity for the catalytic domain. (C) 2008 Elsevier Ltd. All rights reserved.

Sequence and Structural Analysis of the 5 ` Noncoding Region of Hepatitis C Virus in Patients With Chronic Infection

Hepatitis C virus (HCV), exhibits considerable genetic diversity, but presents a relatively well conserved 5 ` noncoding region (5 ` NCR) among all genotypes. In this study, the structural features and translational efficiency of the HCV 5 ` NCR sequences were analyzed using the programs RNAfold, RNAshapes and RNApdist and with a bicistronic dual luciferase expression system, respectively. RNA structure prediction software indicated that base substitutions will alter potentially the 5 ` NCR structure. The heterogeneous sequence observed on 5 ` NCR led to important changes in their translation efficiency in different cell culture lines. Interactions of the viral RNA with cellular transacting factors may vary according to the cell type and viral genome polymorphisms that may result in the translational efficiency observed. J. Med. Virol. 81: 1212-1219, 2009. (C) 2009 Wiley-Liss, Inc.

Braid groups of non-orientable surfaces and the Fadell-Neuwirth short exact sequence

Let M be a compact, connected non-orientable surface without boundary and of genus g >= 3. We investigate the pure braid groups P,(M) of M, and in particular the possible splitting of the Fadell-Neuwirth short exact sequence 1 -> P(m)(M \ {x(1), ..., x(n)}) hooked right arrow P(n+m)(M) (P*) under right arrow P(n)(M) -> 1, where m, n >= 1, and p* is the homomorphism which corresponds geometrically to forgetting the last m strings. This problem is equivalent to that of the existence of a section for the associated fibration p: F(n+m)(M) -> F(n)(M) of configuration spaces, defined by p((x(1), ..., x(n), x(n+1), ..., x(n+m))) = (x(1), ..., x(n)). We show that p and p* admit a section if and only if n = 1. Together with previous results, this completes the resolution of the splitting problem for surface pure braid groups. (C) 2009 Elsevier B.V. All rights reserved.