The impact of measurement errors in the identification of regulatory networks

Background: There are several studies in the literature depicting measurement error in gene expression data and also, several others about regulatory network models. However, only a little fraction describes a combination of measurement error in mathematical regulatory networks and shows how to identify these networks under different rates of noise. Results: This article investigates the effects of measurement error on the estimation of the parameters in regulatory networks. Simulation studies indicate that, in both time series (dependent) and non-time series (independent) data, the measurement error strongly affects the estimated parameters of the regulatory network models, biasing them as predicted by the theory. Moreover, when testing the parameters of the regulatory network models, p-values computed by ignoring the measurement error are not reliable, since the rate of false positives are not controlled under the null hypothesis. In order to overcome these problems, we present an improved version of the Ordinary Least Square estimator in independent (regression models) and dependent (autoregressive models) data when the variables are subject to noises. Moreover, measurement error estimation procedures for microarrays are also described. Simulation results also show that both corrected methods perform better than the standard ones (i.e., ignoring measurement error). The proposed methodologies are illustrated using microarray data from lung cancer patients and mouse liver time series data. Conclusions: Measurement error dangerously affects the identification of regulatory network models, thus, they must be reduced or taken into account in order to avoid erroneous conclusions. This could be one of the reasons for high biological false positive rates identified in actual regulatory network models.

Environmental stresses inhibit splicing in the aquatic fungus Blastocladiella emersonii

Background: Exposure of cells to environmental stress conditions can lead to the interruption of several intracellular processes, in particular those performed by macromolecular complexes such as the spliceosome. Results: During nucleotide sequencing of cDNA libraries constructed using RNA isolated from B. emersonii cells submitted to heat shock and cadmium stress, a large number of ESTs with retained introns was observed. Among the 6,350 ESTs obtained through sequencing of stress cDNA libraries, 181 ESTs presented putative introns (2.9%), while sequencing of cDNA libraries from unstressed B. emersonii cells revealed only 0.2% of ESTs containing introns. These data indicate an enrichment of ESTs with introns in B. emersonii stress cDNA libraries. Among the 85 genes corresponding to the ESTs that retained introns, 19 showed more than one intron and three showed three introns, with intron length ranging from 55 to 333 nucleotides. Canonical splicing junctions were observed in most of these introns, junction sequences being very similar to those found in introns from genes previously characterized in B. emersonii, suggesting that inhibition of splicing during stress is apparently a random process. Confirming our observations, analyses of gpx3 and hsp70 mRNAs by Northern blot and S1 protection assays revealed a strong inhibition of intron splicing in cells submitted to cadmium stress. Conclusion: In conclusion, data indicate that environmental stresses, particularly cadmium treatment, inhibit intron processing in B. emersonii, revealing a new adaptive response to cellular exposure to this heavy metal.

New genes of Xanthomonas citri subsp citri involved in pathogenesis and adaptation revealed by a transposon-based mutant library

Background: Citrus canker is a disease caused by the phytopathogens Xanthomonas citri subsp. citri, Xanthomonas fuscans subsp. aurantifolli and Xanthomonas alfalfae subsp. citrumelonis. The first of the three species, which causes citrus bacterial canker type A, is the most widely spread and severe, attacking all citrus species. In Brazil, this species is the most important, being found in practically all areas where citrus canker has been detected. Like most phytobacterioses, there is no efficient way to control citrus canker. Considering the importance of the disease worldwide, investigation is needed to accurately detect which genes are related to the pathogen-host adaptation process and which are associated with pathogenesis. Results: Through transposon insertion mutagenesis, 10,000 mutants of Xanthomonas citri subsp. citri strain 306 (Xcc) were obtained, and 3,300 were inoculated in Rangpur lime (Citrus limonia) leaves. Their ability to cause citrus canker was analyzed every 3 days until 21 days after inoculation; a set of 44 mutants showed altered virulence, with 8 presenting a complete loss of causing citrus canker symptoms. Sequencing of the insertion site in all 44 mutants revealed that 35 different ORFs were hit, since some ORFs were hit in more than one mutant, with mutants for the same ORF presenting the same phenotype. An analysis of these ORFs showed that some encoded genes were previously known as related to pathogenicity in phytobacteria and, more interestingly, revealed new genes never implicated with Xanthomonas pathogenicity before, including hypothetical ORFs. Among the 8 mutants with no canker symptoms are the hrpB4 and hrpX genes, two genes that belong to type III secretion system (TTSS), two hypothetical ORFS and, surprisingly, the htrA gene, a gene reported as involved with the virulence process in animal-pathogenic bacteria but not described as involved in phytobacteria virulence. Nucleic acid hybridization using labeled cDNA probes showed that some of the mutated genes are differentially expressed when the bacterium is grown in citrus leaves. Finally, comparative genomic analysis revealed that 5 mutated ORFs are in new putative pathogenicity islands. Conclusion: The identification of these new genes related with Xcc infection and virulence is a great step towards the understanding of plant-pathogen interactions and could allow the development of strategies to control citrus canker.

Searching for Moniliophthora perniciosa pathogenicity genes

The basidiomycete Moniliophthora perniciosa is the causal agent of witches` broom disease of Theobroma cacao (cacao). Pathogenesis mechanisms of this hemibiotrophic fungus are largely unknown. An approach to identify putative pathogenicity genes is searching for sequences induced in mycelia grown under in vitro conditions. Using this approach, genes from M. perniciosa induced under limiting nitrogen and light were identified from a cDNA library enriched by suppression subtractive hybridization as potential putative pathogenicity genes. From the 159 identified unique sequences, 59 were annotated and classified by gene ontology. Two sequences were categorized as ""Defence genes, Virulence, and Cell response"" presumably coding for allergenic proteins, whose homologues from other fungi are inducers of animal or plant defences. Differential gene expression was evaluated by quantitative amplification of reversed transcripts (RT-qPCR) of the putative identified genes coding for the two allergenic proteins (Aspf13 and 88KD), and for the enzymes Arylsulfatase (AS); Aryl-Alcohol Oxidase; Aldo-Keto Reductase (AK); Cytochrome P450 (P450); Phenylalanine Ammonia-Lyase; and Peroxidase from mycelia grown under contrasting N concentrations. All genes were validated for differential expression, except for the putative Peroxidase. The same eight genes were analysed for expression in susceptible plants inoculated with M. perniciosa, and six were induced during the early asymptomatic stage of the disease. In infected host tissues, transcripts of 88KD and AS were found more abundant at the biotrophic phase, while those from Aspf13, AK, PAL, and P450 accumulated at the necrotrophic phase, enabling to suggest that mycelia transition from biotrophic to necrotrophic might occur earlier than currently considered. These sequences appeared to be virulence life-style genes, which encode factors or enzymes that enable invasion, colonization or intracellular survival, or manipulate host factors to benefit the pathogen`s own survival in the hostile environment. (C) 2010 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

No-match ORESTES explored as tumor markers

Sequencing technologies and new bioinformatics tools have led to the complete sequencing of various genomes. However, information regarding the human transcriptome and its annotation is yet to be completed. The Human Cancer Genome Project, using ORESTES (open reading frame EST sequences) methodology, contributed to this objective by generating data from about 1.2 million expressed sequence tags. Approximately 30 of these sequences did not align to ESTs in the public databases and were considered no-match ORESTES. On the basis that a set of these ESTs could represent new transcripts, we constructed a cDNA microarray. This platform was used to hybridize against 12 different normal or tumor tissues. We identified 3421 transcribed regions not associated with annotated transcripts, representing 83.3 of the platform. The total number of differentially expressed sequences was 1007. Also, 28 of analyzed sequences could represent noncoding RNAs. Our data reinforces the knowledge of the human genome being pervasively transcribed, and point out molecular marker candidates for different cancers. To reinforce our data, we confirmed, by real-time PCR, the differential expression of three out of eight potentially tumor markers in prostate tissues. Lists of 1007 differentially expressed sequences, and the 291 potentially noncoding tumor markers were provided.

Recombinant rabies virus glycoprotein synthesis in bioreactor by transfected Drosophila melanogaster S2 cells carrying a constitutive or an inducible promoter

S2 cell populations (S2AcRVGP2K and S2MtRVGP-Hy) were selected after transfection of gene expression vectors carrying the cDNA encoding the rabies virus glycoprotein (RVGP) gene under the control of the constitutive (actin) or inductive (metallothionein) promoters. These cell populations were cultivated in a 1 L bioreactor mimicking a large scale bioprocess. Cell cultures were carried out at 90 rpm and monitored/controlled for temperature (28 degrees C) and dissolved oxygen (10 or 50% air saturation). Cell growth attained similar to 1.5-3 x 10(7) cells/mL after 3-4 clays of cultivation. The constitutive synthesis of RVGP in S2AcRVGP2K cells led to values of 0.76 mu g/10(7) cells at day 4 of culture. The RVGP synthesis in S2MtRVGP-Hy cell fraction increased upon CuSO(4) induction attaining specific productivities of 1.5-2 mu g/10(7) cells at clays 4-5. RVGP values in supernatant as a result of cell lysis were always very low (<0.2 mu g/mL) indicating good integrity of cells in culture. Overall the RVGP productivity was of 1.5-3 mg/L. Our data showed an important influence of dissolved oxygen on RVGP synthesis allowing a higher and sustained productivity by S2MtRVGP-Hy cells when cultivated with a DO of 10% air saturation. The RVGP productivity in bioreactors shown here mirrors those previously observed for T-flasks and shaker bottles and allow the preparation of the large RVGP quantities required for studies of structure and function. (C) 2010 Elsevier B.V. All rights reserved.

Transcript profiling of papaya fruit reveals differentially expressed genes associated with fruit ripening

Papaya (Carica papaya L) fruit has a short shelf life due to fast ripening induced by ethylene, but little is known about the genetic control of ripening and attributes of fruit quality. Therefore, we identified ripening-related genes affected by ethylene using cDNA-AFLP (Amplified Fragment Length Polymorphism of cDNA). Transcript profiling of non-induced and ethylene-induced fruit samples was performed, and 71 differentially expressed genes were identified. Among those genes some involved in ethylene biosynthesis, regulation of transcription, and stress responses or plant defence were found (heat shock proteins, polygalacturonase-inhibiting protein, and acyl-CoA oxidases). Several transcription factors were isolated, and except for a 14-3-3 protein, an AP2 domain-containing factor, a salt-tolerant zinc finger protein, and a suppressor of PhyA-105 1, most of them were negatively affected by ethylene, including fragments of transcripts similar to VRN1, and ethylene responsive factors (ERF). With respect to fruit quality, genes related to cell wall structure or metabolism, volatiles or pigment precursors, and vitamin biosynthesis were also found. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

Effective Tityus serrulatus anti-venom produced using the Ts1 component

Scorpion stings are a public health problem in Brazil, with most incidents involving the species Tityus serrulatus. Some T serrulatus toxins may act as immunogens for the production of a specific anti-venom, but many of the component toxins remain poorly characterized. Here, we describe the immunological characteristics of the toxin Ts1 (also known as TsVII and Ts-gamma) and evaluate production of neutralizing antibodies against the crude venom of T serrulatus. Recombinant Ts1 with one copy (Ts1((1))) or two copies in tandem (Ts1((2))) was expressed in BL21 (DE3) cells. Rabbits and mice were immunized with the recombinant proteins (inclusion bodies) and then tested for production of neutralizing antibodies. Neutralization assays showed that anti-Ts1((1)) and anti-Ts1((2)) protected animals challenged with T serrulatus crude venom and native Ts1 Thus, Ts1 could be used in a mixed ""cocktail"" of immunogens for T serrulatus anti-venom production. (C) 2008 Elsevier Ltd. All rights reserved.

Genetic analysis of complement C1s deficiency associated with systemic lupus erythernatosus highlights alternative splicing of normal C1s gene

Deficiencies of complement proteins of the classical pathway are strongly associated with the development of autoimmune diseases. Deficiency of Clr has been observed to occur concomitantly with deficiency in Cls and 9 out of 15 reported cases presented systemic lupus erythernatosus (SLE). Here, we describe a family in which all four children are deficient in Cls but only two of them developed SLE. Hemolytic activity mediated by the alternative and the lectin pathways were normal, but classical pathway activation was absent in all children`s sera. Cls was undetectable, while in the parents` sera it was lower than in the normal controls. The levels of Clr observed in the siblings and parents sera were lower than in the control, while the concentrations of other complement proteins (C3, C4, MBL and MASP-2) were normal in all family members. Impairment of Cls synthesis was observed in the patients` fibroblasts when analyzed by confocal microscopy. We show that all four siblings are homozygous for a mutation at position 938 in exon 6 of the Cls cDNA that creates a premature stop codon. Our investigations led us to reveal the presence of previously uncharacterized splice variants of Cls mRNA transcripts in normal human cells. These variants are derived from the skipping of exon 3 and from the use of an alternative 3` splice site within intron I which increases the size of exon 2 by 87 nucleotides. (c) 2007 Elsevier Ltd. All rights reserved.

Isolation and structural characterization of a new fibrin(ogen)olytic metalloproteinase from Bothrops moojeni snake venom

A proteinase, named BmooMP alpha-I, from the venom of Bothrops moojeni, was purified by DEAE-Sephacel, Sephadex G-75 and heparin-agarose column chromatography. The enzyme was purified to homogeneity as judged by its migration profile in SDS-PAGE stained with coomassie blue, and showed a molecular mass of about 24.5 kDa. Its complete cDNA was obtained by RT-PCR and the 615 bp codified for a mature protein of 205 amino acid residues. The multiple alignment of its deduced amino acid sequence and those of other snake venom metalloproteinases showed a high structural similarly, mainly among class P-IB proteases. The enzyme cleaves the A alpha-chain of fibrinogen first, followed by the B beta-chain, and shows no effects on the gamma-chain. On fibrin, the enzyme hydrolyzed only the beta-chain, leaving the gamma-dimer apparently untouched. It was devoid of phospholipase A(2), hemorrhagic and thrombin-like activities. Like many venom enzymes, it is stable at pH values between 4 and 10 and stable at 70 degrees C for 15 min. The inhibitory effects of EDTA on the fibrinogenolytic activity suggest that BmooMP alpha-I is a metalloproteinase and inhibition by beta-mercaptoethanol revealed the important role of the disulfide bonds in the stabilization of the native structure. Aprotinin and benzamidine, specific serine proteinase inhibitors, had no effect on BmooMP alpha-I activity. Since the BmooMP alpha-I enzyme was found to cause defibrinogenation when administered i.p. on mice, it is expected that it may be of medical interest as a therapeutic agent in the treatment and prevention of arterial thrombosis. (C) 2007 Elsevier Ltd. All rights reserved.

Molecular characterization of BjussuSP-I, a new thrombin-like enzyme with procoagulant and kallikrein-like activity isolated from Bothrops jararacussu snake venom

A thrombin-like enzyme, named BjussuSP-I, isolated from Bothrops jararacussu snake venom, is an acidic single-chain glycoprotein with M-r = 61,000, pI similar to 3.8 and 6% sugar. BjussuSP-I shows high proteolytic activity upon synthetic substrates, such as S-2238 and S-2288. It also shows procoagulant and kallikrein-like activity, but is unable to act on platelets and plasmin. These activities are inhibited by specific inhibitors of this class of enzymes. The complete cDNA sequence of BjussuSP-I with 696 bp encodes open reading frames of 232 amino acid residues, which conserve the common domains of thrombin-like serine proteases. BjussuSP-I shows a high structural homology with other thrombin-like enzymes from snake venoms where common amino acid residues are identified as those corresponding to the catalytic site and subsites S1, S2 and S3 already reported. In this study, we also demonstrated the importance of N-linked glycans, to improve thrombin-like activity of BjussuSP-I toxin. (c) 2007 Elsevier Masson SAS. All rights reserved.

Crotalus durissus collilineatus venom gland transcriptome: Analysis of gene expression profile

Crotalus durissus rattlesnakes are responsible for the most lethal cases of snakebites in Brazil. Crotalus durissus collilineatus subspecies is related to a great number of accidents in Southeast and Central West regions, but few studies on its venom composition have been carried out to date. In an attempt to describe the transcriptional profile of the C. durissus collilineatus venom gland, we generated a cDNA library and the sequences obtained could be identified by similarity searches on existing databases. Out of 673 expressed sequence tags (ESTs) 489 produced readable sequences comprising 201 singletons and 47 clusters of two or more ESTs. One hundred and fifty reads (60.5%) produced significant hits to known sequences. The results showed a predominance of toxin-coding ESTs instead of transcripts coding for proteins involved in all cellular functions. The most frequent toxin was crotoxin, comprising 88% of toxin-coding sequences. Crotoxin B, a basic phospholipase A(2) (PLA(2)) subunit of crotoxin, was represented in more variable forms comparing to the non-enzymatic subunit (crotoxin A), and most sequences coding this molecule were identified as CB1 isoform from Crotalus durissus terrificus venom. Four percent of toxin-related sequences in this study were identified as growth factors, comprising five sequences for vascular endothelial growth factor (VEGF) and one for nerve growth factor (NGF) that showed 100% of identity with C. durissus terrificus NGF. We also identified two clusters for metalloprotease from PII class comprising 3% of the toxins, and two for serine proteases, including gyroxin (2.5%). The remaining 2.5% of toxin-coding ESTs represent singletons identified as homologue sequences to cardiotoxin, convulxin, angiotensin-converting enzyme inhibitor and C-type natriuretic peptide, Ohanin, crotamin and PLA(2) inhibitor. These results allowed the identification of the most common classes of toxins in C. durissus collilineatus snake venom, also showing some unknown classes for this subspecies and even for C. durissus species, such as cardiotoxins and VEGF. (C) 2009 Published by Elsevier Masson SAS.

Structural and functional properties of Bp-LAAO, a new L-amino acid oxidase isolated from Bothrops pauloensis snake venom

An L-amino acid oxidase (Bp-LAAO) from Bothrops pauloensis snake venom was highly purified using sequential chromatography steps on CM-Sepharose, Phenyl-Sepharose CL4B, Benzamidine Sepharose and C18 reverse-phase HPLC. Purified Bp-LAAO showed to be a homodimeric acidic glycoprotein with molecular weight around 65 kDa under reducing conditions in SDS-PAGE. The best substrates for Bp-LAAO were L-Met, L-Leu, L-Phe and L-Ile and the enzyme showed a strong reduction of its catalytic activity upon L-Met and L-Phe substrates at extreme temperatures. Bp-LAAO showed leishmanicidal, antitumoral and bactericidal activities dose dependently. Bp-LAAO induced platelet aggregation in platelet-rich plasma and this activity was inhibited by catalase. Bp-LAAC-cDNA of 1548 bp codified a mature protein with 516 amino acid residues corresponding to a theoretical isoelectric point and molecular weight of 6.3 and 58 kDa, respectively. Additionally, structural and phylogenetic studies identified residues under positive selection and their probable location in Elp-LAAO and other snake venom LAAOs (svLAAOs). Structural and functional investigations of these enzymes can contribute to the advancement of toxinology and to the elaboration of novel therapeutic agents. (C) 2009 Elsevier Masson SAS. All rights reserved.

Molecular cloning and biochemical characterization of a myotoxin inhibitor from Bothrops alternatus snake plasma

Phospholipases A(2) (PLA(2)s) are important components of Bothrops snake venoms, that can induce several effects on envenomations such as myotoxicity, inhibition or induction of platelet aggregation and edema. It is known that venomous and non-venomous snakes present PLA(2) inhibitory proteins (PLIs) in their blood plasma. An inhibitory protein that neutralizes the enzymatic and toxic activities of several PLA2s from Bothrops venoms was isolated from Bothrops alternatus snake plasma by affinity chromatography using the immobilized myotoxin BthTX-I on CNBr-activated Sepharose. Biochemical characterization of this inhibitory protein, denominated alpha BaltMIP, showed it to be a glycoprotein with Mr of similar to 24,000 for the monomeric subunit. CD spectra of the PLA(2)/inhibitor complexes are considerably different from those corresponding to the individual proteins and data deconvolution suggests that the complexes had a relative gain of helical structure elements in comparison to the individual protomers, which may indicate a more compact structure upon complexation. Theoretical and experimental structural studies performed in order to obtain insights into the structural features of aBaltMIP indicated that this molecule may potentially trimerize in solution, thus strengthening the hypothesis previously raised by other authors about snake PLIs oligomerization. (C) 2010 Elsevier Masson SAS. All rights reserved.

Enzymatic and structural characterization of a basic phospholipase A(2) from the sea anemone Condylactis gigantea

This work aimed at the isolation and structural/functional characterization of a phospholipase A(2) (CgPLA(2)) from the extract of the anemone Condylactis gigantea. CgPLA2 was isolated with a high purity level through three chromatographic steps, showing pT8.6 and molecular weights of 14,500 and 29,000 for the monomer and dimer, respectively. CgPLA2 showed a high catalytic activity upon fluorescent phospholipids inducing no direct hemolytic activity. This enzyme, which is Ca2+-dependent, showed a lower stability against temperature and pH variations when compared with snake venom enzymes. The enzymatic activity was significantly reduced or completely abolished after chemical modification of CgPLA2 with BPB. Its cDNA was then obtained, with 357 base pairs which codified for a mature protein of 119 amino acid residues. A comparative analysis of the primary structure of CgPLA2 revealed 84%, 61%, 43% and 42% similarity to the PLA2s from Adamsia carciniopados, Nematostella vectensis, Vipera russelli russelli and Both raps jararacussu, respectively. (C) 2010 Elsevier Masson SAS. All rights reserved.