Identification and Phylogenetic Analysis of Heme Synthesis Genes in Trypanosomatids and Their Bacterial Endosymbionts

It has been known for decades that some insect-infecting trypanosomatids can survive in culture without heme supplementation while others cannot, and that this capability is associated with the presence of a betaproteobacterial endosymbiont in the flagellate's cytoplasm. However, the specific mechanisms involved in this process remained obscure. In this work, we sequence and phylogenetically analyze the heme pathway genes from the symbionts and from their hosts, as well as from a number of heme synthesis-deficient Kinetoplastida. Our results show that the enzymes responsible for synthesis of heme are encoded on the symbiont genomes and produced in close cooperation with the flagellate host. Our evidence suggests that this synergistic relationship is the end result of a history of extensive gene loss and multiple lateral gene transfer events in different branches of the phylogeny of the Trypanosomatidae.

Var transcription profiling of Plasmodium falciparum 3D7: assignment of cytoadherent phenotypes to dominant transcripts

Background: Cytoadherence of Plasmodium falciparum-infected red blood cells is mediated by var gene-encoded P. falciparum erythrocyte membrane protein-1 and host receptor preference depends in most cases on which of the 50-60 var genes per genome is expressed. Enrichment of phenotypically homogenous parasites by panning on receptor expressing cells is fundamental for the identification of the corresponding var transcript. Methods: P. falciparum 3D7 parasites were panned on several transfected CHO-cell lines and their var transcripts analysed by i) reverse transcription/PCR/cloning/sequencing using a universal DBL alpha specific oligonucleotide pair and ii) by reverse transcription followed by quantitative PCR using 57 different oligonucleotide pairs. Results: Each cytoadherence selected parasite line also adhered to untransfected CHO-745 cells and upregulation of the var gene PFD995/PFD1000c was consistently associated with cytoadherence to all but one CHO cell line. In addition, parasites panned on different CHO cell lines revealed candidate var genes which reproducibly associated to the respective cytoadherent phenotype. The transcription profile obtained by RT-PCR/cloning/sequencing differed significantly from that of RT-quantitative PCR. Conclusion: Transfected CHO cell lines are of limited use for the creation of monophenotypic cytoadherent parasite lines. Nevertheless, 3D7 parasites can be reproducibly selected for the transcription of different determined var genes without genetic manipulation. Most importantly, var transcription analysis by RT-PCR/cloning/sequencing may lead to erroneous interpretation of var transcription profiles.

Leishmania-Specific Surface Antigens Show Sub-Genus Sequence Variation and Immune Recognition

Background: A family of hydrophilic acylated surface (HASP) proteins, containing extensive and variant amino acid repeats, is expressed at the plasma membrane in infective extracellular (metacyclic) and intracellular (amastigote) stages of Old World Leishmania species. While HASPs are antigenic in the host and can induce protective immune responses, the biological functions of these Leishmania-specific proteins remain unresolved. Previous genome analysis has suggested that parasites of the sub-genus Leishmania (Viannia) have lost HASP genes from their genomes. Methods/Principal Findings: We have used molecular and cellular methods to analyse HASP expression in New World Leishmania mexicana complex species and show that, unlike in L. major, these proteins are expressed predominantly following differentiation into amastigotes within macrophages. Further genome analysis has revealed that the L. (Viannia) species, L. (V.) braziliensis, does express HASP-like proteins of low amino acid similarity but with similar biochemical characteristics, from genes present on a region of chromosome 23 that is syntenic with the HASP/SHERP locus in Old World Leishmania species and the L. (L.) mexicana complex. A related gene is also present in Leptomonas seymouri and this may represent the ancestral copy of these Leishmania-genus specific sequences. The L. braziliensis HASP-like proteins (named the orthologous (o) HASPs) are predominantly expressed on the plasma membrane in amastigotes and are recognised by immune sera taken from 4 out of 6 leishmaniasis patients tested in an endemic region of Brazil. Analysis of the repetitive domains of the oHASPs has shown considerable genetic variation in parasite isolates taken from the same patients, suggesting that antigenic change may play a role in immune recognition of this protein family. Conclusions/Significance: These findings confirm that antigenic hydrophilic acylated proteins are expressed from genes in the same chromosomal region in species across the genus Leishmania. These proteins are surface-exposed on amastigotes (although L. (L.) major parasites also express HASPB on the metacyclic plasma membrane). The central repetitive domains of the HASPs are highly variant in their amino acid sequences, both within and between species, consistent with a role in immune recognition in the host.

Supramolecular polymorphism of DNA in non-cationic L(alpha) lipid phases

The structure of a complex between hydrated DNA and a non-cationic lipid is studied, including its phase diagram. The complex is spontaneously formed by adding DNA fragments (ca. 150 base pairs in length) to non-cationic lipids and water. The self-assembly process often leads to highly ordered structures. The structures were studied by combining X-ray scattering, fluorescence and polarized microscopy, as well as freeze-fracture experiments with transmission electron microscopy. We observe a significant increase of the smectic order as DNA is incorporated into the water layers of the lamellar host phase, and stabilization of single phase domains for large amounts of DNA. The effect of confinement on DNA ordering is investigated by varying the water content, following three dilution lines. A rich polymorphism is found, ranging from weakly correlated DNA-DNA in-plane organizations to highly ordered structures, where transmembrane correlations lead to the formation of columnar rectangular and columnar hexagonal superlattices of nucleotides embedded between lipid lamellae. From these observations, we suggest that addition of DNA to the lamellar phase significantly restricts membrane fluctuations above a certain concentration and helps the formation of the lipoplex. The alteration of membrane steric interactions, together with the appearance of interfacial interactions between membranes and DNA molecules may be a relevant mechanism for the emergence of highly ordered structures in the concentrated regime.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Identification of COL6A1 as a differentially expressed gene in human astrocytomas

Diffuse infiltrating gliomas are the most common tumors of the central nervous system. Gliomas are classified by the WHO according to their histopathological and clinical characteristics into four classes: grade I (pilocytic astrocytoma), grade II (diffuse astrocytoma), grade III (anaplastic astrocytoma), and grade IV (glioblastoma multiforme). Several genes have already been correlated with astrocytomas, but many others are yet to be uncovered. By analyzing the public SAGE data from 21 patients, comprising low malignant grade astrocytomas and glioblastomas, we found COL6A1 to be differentially expressed, confirming this finding by real time RT-PCR in 66 surgical samples. To the best of our knowledge, COL6A1 has never been described in gliomas. The expression of this gene has significantly different means when normal glia is compared with low-grade astrocytomas (grades I and II) and high-grade astrocytomas (grades III and IV), with a tendency to be greater in higher grade samples, thus rendering it a powerful tumor marker.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Construct and Compare Gene Coexpression Networks with DAPfinder and DAPview

Background: DAPfinder and DAPview are novel BRB-ArrayTools plug-ins to construct gene coexpression networks and identify significant differences in pairwise gene-gene coexpression between two phenotypes. Results: Each significant difference in gene-gene association represents a Differentially Associated Pair (DAP). Our tools include several choices of filtering methods, gene-gene association metrics, statistical testing methods and multiple comparison adjustments. Network results are easily displayed in Cytoscape. Analyses of glioma experiments and microarray simulations demonstrate the utility of these tools. Conclusions: DAPfinder is a new friendly-user tool for reconstruction and comparison of biological networks.

Beta-2 adrenergic receptor gene polymorphisms Gln27Glu, Arg16Gly in patients with heart failure

Background -: Beta-2 adrenergic receptor gene polymorphisms Gln27Glu, Arg16Gly and Thr164Ile were suggested to have an effect in heart failure. We evaluated these polymorphisms relative to clinical characteristics and prognosis of alarge cohort of patients with heart failure of different etiologies. Methods -: We studied 501 patients with heart failure of different etiologies. Mean age was 58 years (standard deviation 14.4 years), 298 (60%) were men. Polymorphisms were identified by polymerase chain reaction-restriction fragment length polymorphism. Results -: During the mean follow-up of 12.6 months (standard deviation 10.3 months), 188 (38%) patients died. Distribution of genotypes of polymorphism Arg16Gly was different relative to body mass index (chi(2) = 9.797; p = 0.04). Overall the probability of survival was not significantly predicted by genotypes of Gln27Glu, Arg16Gly, or Thr164Ile. Allele and haplotype analysis also did not disclose any significant difference regarding mortality. Exploratory analysis through classification trees pointed towards a potential association between the Gln27Glu polymorphism and mortality in older individuals. Conclusion -: In this study sample, we were not able to demonstrate an overall influence of polymorphisms Gln27Glu and Arg16Gly of beta-2 receptor gene on prognosis. Nevertheless, Gln27Glu polymorphism may have a potential predictive value in older individuals.

Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp aurantifolii

Background: Citrus canker is a disease that has severe economic impact on the citrus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri, whose genome sequence was made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C. Results: We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on type III secreted effector repertoires. In addition to pthA, already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI, are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content differences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, type IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP, which codes for a natriuretic protein. Conclusion: We have uncovered numerous commonalities and differences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. The gained knowledge will be instrumental for improving citrus canker control.

Structure-Function Analysis of the HrpB2-HrcU Interaction in the Xanthomonas citri Type III Secretion System

Bacterial type III secretion systems deliver protein virulence factors to host cells. Here we characterize the interaction between HrpB2, a small protein secreted by the Xanthomonas citri subsp. citri type III secretion system, and the cytosolic domain of the inner membrane protein HrcU, a paralog of the flagellar protein FlhB. We show that a recombinant fragment corresponding to the C-terminal cytosolic domain of HrcU produced in E. coli suffers cleavage within a conserved Asn264-Pro265-Thr266-His267 (NPTH) sequence. A recombinant HrcU cytosolic domain with N264A, P265A, T266A mutations at the cleavage site (HrcU(AAAH)) was not cleaved and interacted with HrpB2. Furthermore, a polypeptide corresponding to the sequence following the NPTH cleavage site also interacted with HrpB2 indicating that the site for interaction is located after the NPTH site. Non-polar deletion mutants of the hrcU and hrpB2 genes resulted in a total loss of pathogenicity in susceptible citrus plants and disease symptoms could be recovered by expression of HrpB2 and HrcU from extrachromossomal plasmids. Complementation of the Delta hrcU mutant with HrcU(AAAH) produced canker lesions similar to those observed when complemented with wild-type HrcU. HrpB2 secretion however, was significantly reduced in the Delta hrcU mutant complemented with HrcU(AAAH), suggesting that an intact and cleavable NPTH site in HrcU is necessary for total functionally of T3SS in X. citri subsp. citri. Complementation of the Delta hrpB2 X. citri subsp. citri strain with a series of hrpB2 gene mutants revealed that the highly conserved HrpB2 C-terminus is essential for T3SS-dependent development of citrus canker symptoms in planta.

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.

A NOVEL EPIPHYTIC CYANOBACTERIAL SPECIES FROM THE GENUS BRASILONEMA CAUSING DAMAGE TO EUCALYPTUS LEAVES

A cyanobacterial mat colonizing the leaves of Eucalyptus grandis was determined to be responsible for serious damage affecting the growth and development of whole plants under the clonal hybrid nursery conditions. The dominant cyanobacterial species was isolated in BG-11 medium lacking a source of combined nitrogen and identified by cell morphology characters and molecular phylogenetic analysis (16S rRNA gene and cpcBA-IGS sequences). The isolated strain represents a novel species of the genus Brasilonema and is designated Brasilonema octagenarum strain UFV-E1. Thin sections of E. grandis leaves analyzed by light and electron microscopy showed that the B. octagenarum UFV-E1 filaments penetrate into the leaf mesophyll. The depth of infection and the mechanism by which the cyanobacterium invades leaf tissue were not determined. A major consequence of colonization by this cyanobacterium is a reduction in photosynthesis in the host since the cyanobacterial mats decrease the amount of light incident on leaf surfaces. Moreover, the cyanobacteria also interfere with stomatal gas exchange, decreasing CO2 assimilation. To our knowledge, this is the first report of an epiphytic cyanobacterial species causing damage to E. grandis leaves.

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.