Web services at the European bioinformatics institute

Structure of the catalytically active APOBEC3G bound to a DNA oligonucleotide inhibitor reveals tetrahedral geometry of the transition state

APOBEC3 proteins (A3s) are enzymes that catalyze the deamination of cytidine to uridine in single-stranded DNA (ssDNA) substrates, thus playing a key role in innate antiviral immunity. However, the APOBEC3 family has also been linked to many mutational signatures in cancer cells, which has led to an intense interest to develop inhibitors of A3's catalytic activity as therapeutics as well as tools to study A3's biochemistry, structure, and cellular function. Recent studies have shown that ssDNA containing 2'-deoxy-zebularine (dZ-ssDNA) is an inhibitor of A3s such as A3A, A3B, and A3G, although the atomic determinants of this activity have remained unknown. To fill this knowledge gap, we determined a 1.5 Å resolution structure of a dZ-ssDNA inhibitor bound to active A3G. The crystal structure revealed that the activated dZ-H₂O mimics the transition state by coordinating the active site Zn²⁺ and engaging in additional stabilizing interactions, such as the one with the catalytic residue E259. Therefore, this structure allowed us to capture a snapshot of the A3's transition state and suggests that developing transition-state mimicking inhibitors may provide a new opportunity to design more targeted molecules for A3s in the future.

The IPD-IMGT/HLA Database

It is 24 years since the IPD-IMGT/HLA Database, http://www.ebi.ac.uk/ipd/imgt/hla/, was first released, providing the HLA community with a searchable repository of highly curated HLA sequences. The database now contains over 35 000 alleles of the human Major Histocompatibility Complex (MHC) named by the WHO Nomenclature Committee for Factors of the HLA System. This complex contains the most polymorphic genes in the human genome and is now considered hyperpolymorphic. The IPD-IMGT/HLA Database provides a stable and user-friendly repository for this information. Uptake of Next Generation Sequencing technology in recent years has driven an increase in the number of alleles and the length of sequences submitted. As the size of the database has grown the traditional methods of accessing and presenting this data have been challenged, in response, we have developed a suite of tools providing an enhanced user experience to our traditional web-based users while creating new programmatic access for our bioinformatics user base. This suite of tools is powered by the IPD-API, an Application Programming Interface (API), providing scalable and flexible access to the database. The IPD-API provides a stable platform for our future development allowing us to meet the future challenges of the HLA field and needs of the community.

Deciphering Human Leukocyte Antigen Susceptibility Maps From Immunopeptidomics Characterization in Oncology and Infections

Genetic variability across the three major histocompatibility complex (MHC) class I genes (human leukocyte antigen [HLA] A, B, and C) may affect susceptibility to many diseases such as cancer, auto-immune or infectious diseases. Individual genetic variation may help to explain different immune responses to microorganisms across a population. HLA typing can be fast and inexpensive; however, deciphering peptides loaded on MHC-I and II which are presented to T cells, require the design and development of high-sensitivity methodological approaches and subsequently databases. Hence, these novel strategies and databases could help in the generation of vaccines using these potential immunogenic peptides and in identifying high-risk HLA types to be prioritized for vaccination programs. Herein, the recent developments and approaches, in this field, focusing on the identification of immunogenic peptides have been reviewed and the next steps to promote their translation into biomedical and clinical practice are discussed.

Crystal Structure of a Soluble APOBEC3G Variant Suggests ssDNA to Bind in a Channel that Extends between the Two Domains

APOBEC3G (A3G) is a single-stranded DNA (ssDNA) cytosine deaminase that can restrict HIV-1 infection by mutating the viral genome. A3G consists of a non-catalytic N-terminal domain (NTD) and a catalytic C-terminal domain (CTD) connected by a short linker. While the CTD catalyzes cytosine deamination, the NTD is believed to provide additional affinity for ssDNA. Structures of both A3G domains have been solved individually; however, a full-length A3G structure has been challenging. Recently, crystal structures of full-length rhesus macaque A3G variants were solved which suggested dimerization mechanisms and RNA binding surfaces, whereas the dimerization appeared to compromise catalytic activity. We determined the crystal structure of a soluble variant of human A3G (sA3G) at 2.5 Å and from these data generated a model structure of wild-type A3G. This model demonstrated that the NTD was rotated 90° relative to the CTD along the major axis of the molecule, an orientation that forms a positively charged channel connected to the CTD catalytic site, consisting of NTD loop-1 and CTD loop-3. Structure-based mutations, in vitro deamination and DNA binding assays, and HIV-1 restriction assays identify R24, located in the NTD loop-1, as essential to a critical interaction with ssDNA. Furthermore, sA3G was shown to bind a deoxy-cytidine dinucleotide near the catalytic Zn²⁺, yet not in the catalytic position, where the interactions between deoxy-cytidines and CTD loop-1 and loop-7 residues were different from those formed with substrate. These new interactions suggest a mechanism explaining why A3G exhibits a 3' to 5' directional preference in processive deamination.

Metabolic drug targets of the cytosine metabolism pathways in the dromedary camel (Camelus dromedarius) and blood parasite Trypanosoma evansi

Trypanosomiasis is a major illness affecting camels in tropical and subtropical regions. Comparisons of camel and Trypanosoma evansi genomes can lead to the discovery of new drug targets for treating Trypanosoma infections. The synthesis pathways of cytosine, cytidine, cytidine monophosphate (CMP), cytidine diphosphate (CDP), cytidine triphosphate (CTP) deoxycytidine, deoxycytidine monophosphate (dCMP), deoxycytidine diphosphate (dCDP), and deoxycytidine triphosphate (dCTP) were compared in the dromedary camel (Camelus dromedarius) and T. evansi. None of the enzymes involved in cytosine pathway were detected in camels and T. evansi. Notably, cytidine kinase (CK) and 5'-nucleotidase, which interconverts cytidine to CMP, were not detected in T. evansi but were present in camels. UMP/CMP kinase was not predicted in T. evansi. Therefore, the presence of enzymes involved in the CTP synthesis cascade was not predicted in T. evansi. CMP synthesis might also be encoded by other enzymes, e.g., purine nucleotides kinases. Both camel and T. evansi share an efficient enzyme system for converting CDP to CTP. In conclusion, CTP synthase is important for homeostasis of cytosine nucleotides in T. evansi and could be a potential drug target against the parasite. In addition, the inhibition of UMP synthesis might contribute to parasite death as it is a shared source for CTP synthesis.

Immunoinformatics: In Silico Approaches and Computational Design of a Multi-epitope, Immunogenic Protein

Immunoinformatics is a new critical field with several tools and databases that conduct the eyesight of experimental selection and facilitate analysis of the great amount of immunologic data obtained from experimental researches and helps to design and introducing new hypothesis. Given these visages, immunoinformatics seems to be the way that develop and progress the immunological research. Bioinformatics methods and applications are successfully employed in vaccine informatics to assist different sites of the preclinical, clinical, and post-licensure vaccine enterprises. On the other hand, the progression of molecular biology and immunology caused epitope vaccines have become the focus of research on molecular vaccines. Moreover, reverse vaccinology could improve vaccine production and vaccination protocols by in silico prediction of protein-vaccine candidates from genome sequences. B- and T-cell immune epitopes could be predicted by immunoinformatics algorithms and computational methods to improve the vaccine design, protective immunity analysis, assessment of vaccine safety and efficacy, and immunization modeling. This review aims to discuss the power of computational approaches in vaccine design and their relevance to the development of effective vaccines. Furthermore, the various divisions of this field and available tools in each item are introduced and reviewed.

Characterization of gene promoters in pig: conservative elements, regulatory motifs and evolutionary trend

It is vital to understand the conservation and evolution of gene promoter sequences in order to understand environmental adaptation. The level of promoter conservation varies greatly between housekeeping (HK) and tissue-specific (TS) genes, denoting differences in the strength of the evolutionary constraints. Here, we analyzed promoter conservation and evolution to exploit differential regulation between HK and TS genes. The analysis of conserved elements showed CpG islands, short tandem repeats and G-quadruplex sequences are highly enriched in HK promoters relative to TS promoters. In addition, the type and density of regulatory motifs in TS promoters are much higher than HK promoters, indicating that TS genes show more complex regulatory patterns than HK genes. Moreover, the evolutionary dynamics of promoters showed similar evolutionary trend to coding sequences. HK promoters suffer more stringent selective pressure in the long-term evolutionary process. HK genes tend to show increased upstream sequence conservation due to stringent selection pressures acting on the promoter regions. The specificity of TS gene expression may be due to complex regulatory motifs acting in different tissues or conditions. The results from this study can be used to deepen our understanding of adaptive evolution.

Crystal structure of the catalytic domain of HIV-1 restriction factor APOBEC3G in complex with ssDNA

The human APOBEC3G protein is a cytidine deaminase that generates cytidine to deoxy-uridine mutations in single-stranded DNA (ssDNA), and capable of restricting replication of HIV-1 by generating mutations in viral genome. The mechanism by which APOBEC3G specifically deaminates 5′-CC motifs has remained elusive since structural studies have been hampered due to apparently weak ssDNA binding of the catalytic domain of APOBEC3G. We overcame the problem by generating a highly active variant with higher ssDNA affinity. Here, we present the crystal structure of this variant complexed with a ssDNA substrate at 1.86 Å resolution. This structure reveals atomic-level interactions by which APOBEC3G recognizes a functionally-relevant 5′-TCCCA sequence. This complex also reveals a key role of W211 in substrate recognition, implicating a similar recognition in activation-induced cytidine deaminase (AID) with a conserved tryptophan.

APOBEC3G (A3G) is a single-stranded DNA (ssDNA) cytidine deaminase that restricts HIV-1. Here the authors provide molecular insights into A3G substrate recognition by determining the 1.86 Å resolution crystal structure of its catalytic domain bound to ssDNA.

Occurrence of Acanthamoeba genotypes in Central West Malaysian environments

Acanthamoeba species are ubiquitous free-living protozoa that can be found worldwide. Occasionally, it can become parasitic and the causative agent of acanthamoebic keratitis (AK) and Granulomatous Amoebic Encephalitis (GAE) in man. A total of 160 environmental samples and 225 naturally-infected animal corneal swabs were collected for Acanthamoeba cultivation. Acanthamoeba was found to be high in samples collected from environments (85%, 136/160) compared to infected animal corneas (24.89%, 56/225) by microscopic examination. Analysis of nucleotide sequence of 18S rRNA gene of all the 192 cultivable Acanthamoeba isolates revealed 4 genotypes (T3, T4. T5 and T15) with T4 as the most prevalent (69.27%, 133/192) followed by T5 (20.31%), T15 (9.90%) and T3 (0.52%). Genotype T4 was from the strain of A. castellanii U07401 (44.27%), A. castellanii U07409 (20.83%) and A. polyphagaAY026243 (4.17%), but interestingly, only A. castellanii U07401 was detected in naturally infected corneal samples. In environmental samples, T4 was commonly detected in all samples including dry soil, dust, wet debris, wet soil and water. Among the T4, A. castellanii (U07409) strains were detected high occurrence in dry (45%) followed by aquatic (32.50%) and moist (22.50%) samples but however A. castellanii (U07401) strains were dominant in dry samples of soil and dust (93.10%). Subsequently, genotype T5 of A. lenticulata (U94741) strains were dominant in samples collected from aquatic environments (58.97%). In summary, A. castellanii (U07401) strains were found dominant in both environmental and corneal swab samples. Therefore, these strains are possibly the most virulent and dry soil or dusts are the most possible source of Acanthamoeba infection in cats and dogs corneas.

Cryptic gene pools in the Hypericum perforatum-H. maculatum complex: diploid persistence versus trapped polyploid melting

BACKGROUND AND AIMS

In Central Europe Hypericum perforatum and Hypericum maculatum show significant hybridization and introgression as a consequence of Pleistocene range fluctuations, and their gene pools are merging on higher ploidy levels. This paper discusses whether polyploid hybrid gene pools are trapped in the ecological climatic niche space of their diploid ancestors, and tests the idea of geographical parthenogenesis.

METHODS

DNA sequence information of nuclear ribosomal DNA and plastid loci, ploidy level estimates and ecological niche modelling are used to characterize the various diploid and polyploid gene pools and unravel spatio-temporal patterns of gene flow among them.

KEY RESULTS

On the diploid level, the three gene pools are clearly distinct between and within species of H. perforatum (two gene pools) and H. maculatum, and their divergence dates back to the first half of the Pleistocene. All polyploids in Central Europe show high levels of past and contemporary gene flow between all three gene pools. The correlation of genetic and geographical distances breaks down if the latter is larger than 250 km, indicating recent and ongoing gene flow. The two species are ecologically differentiated, but in particular hybrids among all three gene pools do not show significant niche differences compared to their parental gene pools, except for some combinations with H. maculatum.

CONCLUSIONS

Inter- and intraspecific gene flow between inter- and intra-species gene pools is limited on the diploid level, and the geographical distribution of the diploids largely reflects Pleistocene evolutionary history. Secondary contact promoted hybridization and introgression on the polyploid level, enabling offspring to escape the diploid gene pools. However, the hybrid polyploids do not show significant niche differences compared to their diploid progenitors. It is concluded that the observed absence of niche divergence has precluded further differentiation and geographical partitioning of new polyploid lineages being effectively separated from the parental lines. The predominantly apomictic reproducing polyploids are trapped in the polyploid gene pool and the ecological climatic niche space of their diploid ancestors.

The transcriptional landscape of Rhizoctonia solani AG1-IA during infection of soybean as defined by RNA-seq

Rhizoctonia solani Kühn infects most plant families and can cause significant agricultural yield losses worldwide; however, plant resistance to this disease is rare and short-lived, and therefore poorly understood, resulting in the use of chemical pesticides for its control. Understanding the functional responses of this pathogen during host infection can help elucidate the molecular mechanisms that are necessary for successful host invasion. Using the pathosystem model soybean-R. solani anastomosis group AG1-IA, we examined the global transcriptional responses of R. solani during early and late infection stages of soybean by applying an RNA-seq approach. Approximately, 148 million clean paired-end reads, representing 93% of R. solani AG1-IA genes, were obtained from the sequenced libraries. Analysis of R. solani AG1-IA transcripts during soybean invasion revealed that most genes were similarly expressed during early and late infection stages, and only 11% and 15% of the expressed genes were differentially expressed during early and late infection stages, respectively. Analyses of the differentially expressed genes (DEGs) revealed shifts in molecular pathways involved in antibiotics biosynthesis, amino acid and carbohydrate metabolism, as well as pathways involved in antioxidant production. Furthermore, several KEGG pathways were unique to each time point, particularly the up-regulation of genes related to toxin degradation (e.g., nicotinate and nicotinamid metabolism) at onset of necrosis, and those linked to synthesis of anti-microbial compounds and pyridoxine (vitamin B6) biosynthesis 24 h.p.o. of necrosis. These results suggest that particular genes or pathways are required for either invasion or disease development. Overall, this study provides the first insights into R. solani AG1-IA transcriptome responses to soybean invasion providing beneficial information for future targeted control methods of this successful pathogen.

Comparative genomics reveals convergent evolution between the bamboo-eating giant and red pandas

Phenotypic convergence between distantly related taxa often mirrors adaptation to similar selective pressures and may be driven by genetic convergence. The giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to different families in the order Carnivora, but both have evolved a specialized bamboo diet and adaptive pseudothumb, representing a classic model of convergent evolution. However, the genetic bases of these morphological and physiological convergences remain unknown. Through de novo sequencing the red panda genome and improving the giant panda genome assembly with added data, we identified genomic signatures of convergent evolution. Limb development genes DYNC2H1 and PCNT have undergone adaptive convergence and may be important candidate genes for pseudothumb development. As evolutionary responses to a bamboo diet, adaptive convergence has occurred in genes involved in the digestion and utilization of bamboo nutrients such as essential amino acids, fatty acids, and vitamins. Similarly, the umami taste receptor gene TAS1R1 has been pseudogenized in both pandas. These findings offer insights into genetic convergence mechanisms underlying phenotypic convergence and adaptation to a specialized bamboo diet.

A comparative study on the regulatory region of the PERIOD1 gene among diurnal/nocturnal primates

BACKGROUND

The circadian clock is set up around a 24-h period in humans who are awake in the daytime and sleep in the nighttime, accompanied with physiological and metabolic rhythms. Most haplorhine primates, including humans, are diurnal, while most "primitive" strepsirrhine primates are nocturnal, suggesting primates have evolved from nocturnal to diurnal habits. The mechanisms of physiological changes causing the habits and of genetic changes causing the physiological changes are, however, unknown. To reveal these mechanisms, we focus on the nucleotide sequences of the regulatory region of the PERIOD1 (PER1) gene that is known as one of the key elements of the circadian clock in mammalians.

METHODS

We determined nucleotide sequences of the regulatory region of PER1 concerning the gene expression for six primates and compared those with those of eight primates from the international DNA database. Based on the sequence data, we constructed a phylogenetic tree including both the diurnal/nocturnal species and investigated the guanine and cytosine (GC) content in the regulatory region.

RESULTS

The motif sequences regulating gene expression were evolutionary conservative in the primates examined. The phylogenetic tree simply showed phylogenetic relationship among the species and no branching pattern distinguishable between the diurnal and nocturnal groups. We found two cores showing a statistically significant difference between the diurnal and the nocturnal habits related to the GC contents of the regulatory region of PER1.

CONCLUSION

Our results suggest the possibility that the two cores in the upstream region of PER1 are related to the regulation of gene expression leading to behavioral differences between diurnal and nocturnal primates.

Sex- and tissue-specific profiles of chemosensory gene expression in a herbivorous gall-inducing fly (Diptera: Cecidomyiidae)

BACKGROUND

The chemical senses of insects mediate behaviors that are closely linked to survival and reproduction. The order Diptera contains two model organisms, the vinegar fly Drosophila melanogaster and the mosquito Anopheles gambiae, whose chemosensory genes have been extensively studied. Representing a third dipteran lineage with an interesting phylogenetic position, and being ecologically distinct by feeding on plants, the Hessian fly (Mayetiola destructor Say, Diptera: Cecidomyiidae) genome sequence has recently become available. Among plant-feeding insects, the Hessian fly is unusual in 'reprogramming' the plant to create a superior food and in being the target of plant resistance genes, a feature shared by plant pathogens. Chemoreception is essential for reproductive success, including detection of sex pheromone and plant-produced chemicals by males and females, respectively.

RESULTS

We identified genes encoding 122 odorant receptors (OR), 28 gustatory receptors (GR), 39 ionotropic receptors (IR), 32 odorant binding proteins, and 7 sensory neuron membrane proteins in the Hessian fly genome. We then mapped Illumina-sequenced transcriptome reads to the genome to explore gene expression in male and female antennae and terminal abdominal segments. Our results reveal that a large number of chemosensory genes have up-regulated expression in the antennae, and the expression is in many cases sex-specific. Sex-specific expression is particularly evident among the Or genes, consistent with the sex-divergent olfactory-mediated behaviors of the adults. In addition, the large number of Ors in the genome but the reduced set of Grs and divergent Irs suggest that the short-lived adults rely more on long-range olfaction than on short-range gustation. We also report up-regulated expression of some genes from all chemosensory gene families in the terminal segments of the abdomen, which play important roles in reproduction.

CONCLUSIONS

We show that a large number of the chemosensory genes in the Hessian fly genome have sex- and tissue-specific expression profiles. Our findings provide the first insights into the molecular basis of chemoreception in plant-feeding flies, representing an important advance toward a more complete understanding of olfaction in Diptera and its links to ecological specialization.

Gene prioritization of resistant rice gene against Xanthomas oryzae pv. oryzae by using text mining technologies

To effectively assess the possibility of the unknown rice protein resistant to Xanthomonas oryzae pv. oryzae, a hybrid strategy is proposed to enhance gene prioritization by combining text mining technologies with a sequence-based approach. The text mining technique of term frequency inverse document frequency is used to measure the importance of distinguished terms which reflect biomedical activity in rice before candidate genes are screened and vital terms are produced. Afterwards, a built-in classifier under the chaos games representation algorithm is used to sieve the best possible candidate gene. Our experiment results show that the combination of these two methods achieves enhanced gene prioritization.

Functional analysis of novel polymorphisms in the human SLCO1A2 gene that encodes the transporter OATP1A2

The solute carrier organic anion transporting polypeptide 1A2 (OATP1A2, SLCO1A2) is implicated in the cellular influx of a number of drugs. We identified five novel single nucleotide polymorphisms (SNPs) in coding exons of the SLCO1A2 gene in a cohort of subjects: G550A, G553A, G673A, A775C, and G862A, that encoded the OATP1A2 variants E184K, D185N, V255I, T259P, and D288N, respectively. The function and expression of these variant transporters were assessed in HEK-293 cells. We found that the novel variants, E184K, D185N, T259P, and D288N, were associated with impaired estrone-3-sulfate, imatinib, and methotrexate transport (∼20-50% of wild-type control); function was retained by OATP1A2-V255I. From biotinylation assays, the decreased function of these variants was due, at least in part, to impaired plasma membrane expression. The four loss-of-function variants were studied further using mutagenesis to produce variants that encode residues with different charges or steric properties. From immunoblotting, the replacement of negatively charged residues at amino acid positions 184 and 185 impaired membrane expression, while either a positive or negative charge at residue 288 supported the correct membrane targeting of OATP1A2. Replacement of T259 with bulky residues disrupted transporter stability. From molecular models, E184, D185, and D288 were located near several charged residues such that intramolecular ionic interactions may stabilize the transporter structure. Individuals who carry these novel SNPs in the SLCO1A2 gene may be at risk from impaired efficacy or enhanced toxicity during treatment with drugs that are substrates for OATP1A2.

Use of SciDBMaker as Tool for the Design of Specialized Biological Databases

The exponential growth of molecular biology research in recent decades has brought concomitant growth in the number and size of genomic and proteomic databases used to interpret experimental findings. Particularly, growth of protein sequence records created the need for smaller and manually annotated databases. Since scientists are continually developing new specific databases to enhance their understanding of biological processes, the authors created SciDBMaker to provide a tool for easy building of new specialized protein knowledge bases. This chapter also suggests best practices for specialized biological databases design, and provides examples for the implementation of these practices.

Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae)

Background

The European spruce bark beetle, Ips typographus, and the North American mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae), are severe pests of coniferous forests. Both bark beetle species utilize aggregation pheromones to coordinate mass-attacks on host trees, while odorants from host and non-host trees modulate the pheromone response. Thus, the bark beetle olfactory sense is of utmost importance for fitness. However, information on the genes underlying olfactory detection has been lacking in bark beetles and is limited in Coleoptera. We assembled antennal transcriptomes from next-generation sequencing of I. typographus and D. ponderosae to identify members of the major chemosensory multi-gene families.

Results

Gene ontology (GO) annotation indicated that the relative abundance of transcripts associated with specific GO terms was highly similar in the two species. Transcripts with terms related to olfactory function were found in both species. Focusing on the chemosensory gene families, we identified 15 putative odorant binding proteins (OBP), 6 chemosensory proteins (CSP), 3 sensory neuron membrane proteins (SNMP), 43 odorant receptors (OR), 6 gustatory receptors (GR), and 7 ionotropic receptors (IR) in I. typographus; and 31 putative OBPs, 11 CSPs, 3 SNMPs, 49 ORs, 2 GRs, and 15 IRs in D. ponderosae. Predicted protein sequences were compared with counterparts in the flour beetle, Tribolium castaneum, the cerambycid beetle, Megacyllene caryae, and the fruit fly, Drosophila melanogaster. The most notable result was found among the ORs, for which large bark beetle-specific expansions were found. However, some clades contained receptors from all four beetle species, indicating a degree of conservation among some coleopteran OR lineages. Putative GRs for carbon dioxide and orthologues for the conserved antennal IRs were included in the identified receptor sets.

Conclusions

The protein families important for chemoreception have now been identified in three coleopteran species (four species for the ORs). Thus, this study allows for improved evolutionary analyses of coleopteran olfaction. Identification of these proteins in two of the most destructive forest pests, sharing many semiochemicals, is especially important as they might represent novel targets for population control.

Expression of hsp70, hsp100 and ubiquitin in Aloe barbadensis Miller under direct heat stress and under temperature acclimation conditions

KEY MESSAGE : The study determined the tolerance of Aloe vera to high temperature, focusing on the expression of hsp70 , hsp100 and ubiquitin genes. These were highly expressed in plants acclimated at 35 °C prior to a heat shock of 45 °C. Aloe barbadensis Miller (Aloe vera), a CAM plant, was introduced into Chile in the semiarid IV and III Regions, which has summer diurnal temperature fluctuations of 25 to 40 °C and annual precipitation of 40 mm (dry years) to 170 mm (rainy years). The aim of this study was to investigate how Aloe vera responds to water and heat stress, focusing on the expression of heat shock genes (hsp70, hsp100) and ubiquitin, which not studied before in Aloe vera. The LT(50) of Aloe vera was determined as 53.2 °C. To study gene expression by semi-quantitative RT-PCR, primers were designed against conserved regions of these genes. Sequencing the cDNA fragments for hsp70 and ubiquitin showed a high identity, over 95 %, with the genes from cereals. The protein sequence of hsp70 deduced from the sequence of the cDNA encloses partial domains for binding ATP and the substrate. The protein sequence of ubiquitin deduced from the cDNA encloses a domain for interaction with the enzymes E2, UCH and CUE. The expression increased with temperature and water deficit. Hsp70 expression at 40-45 °C increased 50 % over the controls, while the expression increased by 150 % over the controls under a water deficit of 50 % FC. The expression of all three genes was also studied under 2 h of acclimation at 35 or 40 °C prior to a heat shock at 45 °C. Under these conditions, the plants showed greater expression of all genes than when they were subjected to direct heat stress.

MAPI: a software framework for distributed biomedical applications

Background

The amount of web-based resources (databases, tools etc.) in biomedicine has increased, but the integrated usage of those resources is complex due to differences in access protocols and data formats. However, distributed data processing is becoming inevitable in several domains, in particular in biomedicine, where researchers face rapidly increasing data sizes. This big data is difficult to process locally because of the large processing, memory and storage capacity required.

Results

This manuscript describes a framework, called MAPI, which provides a uniform representation of resources available over the Internet, in particular for Web Services. The framework enhances their interoperability and collaborative use by enabling a uniform and remote access. The framework functionality is organized in modules that can be combined and configured in different ways to fulfil concrete development requirements.

Conclusions

The framework has been tested in the biomedical application domain where it has been a base for developing several clients that are able to integrate different web resources. The MAPI binaries and documentation are freely available at http://www.bitlab-es.com/mapi under the Creative Commons Attribution-No Derivative Works 2.5 Spain License. The MAPI source code is available by request (GPL v3 license).

D471G mutation in LCMV-NP affects its ability to self-associate and results in a dominant negative effect in viral RNA synthesis

Arenaviruses merit significant interest because several family members are etiological agents of severe hemorrhagic fevers, representing a major burden to public health. Currently, there are no FDA-licensed vaccines against arenaviruses and the only available antiviral therapy is limited to the use of ribavirin that is partially effective. Arenavirus nucleoprotein (NP) is found associated with the genomic RNA forming the viral ribonucleoproteins (vRNPs) that together with the polymerase (L) direct viral replication and transcription. Virion formation requires the recruitment of vRNPs into budding sites, a process in which the arenavirus matrix-like protein (Z) plays a major role. Therefore, proper NP-NP and NP-Z interactions are required for the generation of infectious progeny. In this work we demonstrate the role of the amino acid residue D471 in the self-association of lymphocytic choriomeningitis virus nucleoprotein (LCMV-NP). Amino acid substitutions at this position abrogate NP oligomerization, affecting its ability to mediate replication and transcription of a minigenome reporter plasmid. However, its ability to interact with the Z protein, counteract the cellular interferon response and bind to dsRNA analogs was retained. Additionally, we also document the dominant negative effect of D471G mutation on viral infection, suggesting that NP self-association is an excellent target for the development of new antivirals against arenaviruses.

SUB1A-dependent and -independent mechanisms are involved in the flooding tolerance of wild rice species

Crop tolerance to flooding is an important agronomic trait. Although rice (Oryza sativa) is considered a flood-tolerant crop, only limited cultivars display tolerance to prolonged submergence, which is largely attributed to the presence of the SUB1A gene. Wild Oryza species have the potential to unveil adaptive mechanisms and shed light on the basis of submergence tolerance traits. In this study, we screened 109 Oryza genotypes belonging to different rice genome groups for flooding tolerance. Oryza nivara and Oryza rufipogon accessions, belonging to the A-genome group, together with Oryza sativa, showed a wide range of submergence responses, and the tolerance-related SUB1A-1 and the intolerance-related SUB1A-2 alleles were found in tolerant and sensitive accessions, respectively. Flooding-tolerant accessions of Oryza rhizomatis and Oryza eichingeri, belonging to the C-genome group, were also identified. Interestingly, SUB1A was absent in these species, which possess a SUB1 orthologue with high similarity to O. sativa SUB1C. The expression patterns of submergence-induced genes in these rice genotypes indicated limited induction of anaerobic genes, with classical anaerobic proteins poorly induced in O. rhizomatis under submergence. The results indicated that SUB1A-1 is not essential to confer submergence tolerance in the wild rice genotypes belonging to the C-genome group, which show instead a SUB1A-independent response to submergence.

MotViz: a tool for sequence motif prediction in parallel to structural visualization and analyses

Linking similar proteins structurally is a challenging task that may help in finding the novel members of a protein family. In this respect, identification of conserved sequence can facilitate understanding and classifying the exact role of proteins. However, the exact role of these conserved elements cannot be elucidated without structural and physiochemical information. In this work, we present a novel desktop application MotViz designed for searching and analyzing the conserved sequence segments within protein structure. With MotViz, the user can extract a complete list of sequence motifs from loaded 3D structures, annotate the motifs structurally and analyze their physiochemical properties. The conservation value calculated for an individual motif can be visualized graphically. To check the efficiency, predicted motifs from the data sets of 9 protein families were analyzed and MotViz algorithm was more efficient in comparison to other online motif prediction tools. Furthermore, a database was also integrated for storing, retrieving and performing the detailed functional annotation studies. In summary, MotViz effectively predicts motifs with high sensitivity and simultaneously visualizes them into 3D strucures. Moreover, MotViz is user-friendly with optimized graphical parameters and better processing speed due to the inclusion of a database at the back end. MotViz is available at http://www.fi-pk.com/motviz.html.

Comparative transcriptome analysis of two olive cultivars in response to NaCl-stress

BACKGROUND

Olive (Olea europaea L.) cultivation is rapidly expanding and low quality saline water is often used for irrigation. The molecular basis of salt tolerance in olive, though, has not yet been investigated at a system level. In this study a comparative transcriptomics approach was used as a tool to unravel gene regulatory networks underlying salinity response in olive trees by simulating as much as possible olive growing conditions in the field. Specifically, we investigated the genotype-dependent differences in the transcriptome response of two olive cultivars, a salt-tolerant and a salt-sensitive one.

METHODOLOGY/PRINCIPAL FINDINGS

A 135-day long salinity experiment was conducted using one-year old trees exposed to NaCl stress for 90 days followed by 45 days of post-stress period during the summer. A cDNA library made of olive seedling mRNAs was sequenced and an olive microarray was constructed. Total RNA was extracted from root samples after 15, 45 and 90 days of NaCl-treatment as well as after 15 and 45 days of post-treatment period and used for microarray hybridizations. SAM analysis between the NaCl-stress and the post-stress time course resulted in the identification of 209 and 36 differentially expressed transcripts in the salt-tolerant and salt-sensitive cultivar, respectively. Hierarchical clustering revealed two major, distinct clusters for each cultivar. Despite the limited number of probe sets, transcriptional regulatory networks were constructed for both cultivars while several hierarchically-clustered interacting transcription factor regulators such as JERF and bZIP homologues were identified.

CONCLUSIONS/SIGNIFICANCE

A systems biology approach was used and differentially expressed transcripts as well as regulatory interactions were identified. The comparison of the interactions among transcription factors in olive with those reported for Arabidopsis might indicate similarities in the response of a tree species with Arabidopsis at the transcriptional level under salinity stress.

Web tools for predictive toxicology model building

INTRODUCTION

The development and use of web tools in chemistry has accumulated more than 15 years of history already. Powered by the advances in the Internet technologies, the current generation of web systems are starting to expand into areas, traditional for desktop applications. The web platforms integrate data storage, cheminformatics and data analysis tools. The ease of use and the collaborative potential of the web is compelling, despite the challenges.

AREAS COVERED

The topic of this review is a set of recently published web tools that facilitate predictive toxicology model building. The focus is on software platforms, offering web access to chemical structure-based methods, although some of the frameworks could also provide bioinformatics or hybrid data analysis functionalities. A number of historical and current developments are cited. In order to provide comparable assessment, the following characteristics are considered: support for workflows, descriptor calculations, visualization, modeling algorithms, data management and data sharing capabilities, availability of GUI or programmatic access and implementation details.

EXPERT OPINION

The success of the Web is largely due to its highly decentralized, yet sufficiently interoperable model for information access. The expected future convergence between cheminformatics and bioinformatics databases provides new challenges toward management and analysis of large data sets. The web tools in predictive toxicology will likely continue to evolve toward the right mix of flexibility, performance, scalability, interoperability, sets of unique features offered, friendly user interfaces, programmatic access for advanced users, platform independence, results reproducibility, curation and crowdsourcing utilities, collaborative sharing and secure access.

Arnica (Asteraceae) phylogeny revisited using RPB2: complex patterns and multiple d-paralogues

The region coding for the second largest subunit of RNA polymerase II (RPB2) was explored for resolving interspecific relationships in Arnica and lower level taxa in general. The region between exons 17 and 23 was cloned and sequenced for 33 accessions of Arnica and four outgroup taxa. Three paralogues of the RPB2-d copy (RPB2-dA, B and C) were detected in Arnica and outgroup taxa, indicating that the duplications must have occurred before the divergence of Arnica. Parsimony and Bayesian analyses of separate alignments of the three copies reveal complex patterns in Arnica, likely reflecting a history of lineage sorting in combination with apomixis, polyploidization, and possibly hybridization. Cloned sequences of some taxa do not form monophyletic clades within paralogues, but form multiple strongly supported clades with sequences of other taxa. Some well supported groups are present in more than one paralogue and many groups are in line with earlier hypotheses regarding interspecific relationships within the genus. Low levels of homoplasy in combination with relatively high sequence variation indicates that the introns of the RPB2 region could be suitable for phylogenetic studies in low level taxonomy.

Determinants for activation of the atypical AGC kinase Greatwall during M phase entry

The atypical AGC kinase Greatwall (Gwl) mediates a pathway that prevents the precocious removal of phosphorylations added to target proteins by M phase-promoting factor (MPF); Gwl is thus essential for M phase entry and maintenance. Gwl itself is activated by M phase-specific phosphorylations that are investigated here. Many phosphorylations are nonessential, being located within a long nonconserved region, any part of which can be deleted without effect. Using mass spectrometry and mutagenesis, we have identified 3 phosphorylation sites (phosphosites) critical to Gwl activation (pT193, pT206, and pS883 in Xenopus laevis) located in evolutionarily conserved domains that differentiate Gwl from related kinases. We propose a model in which the initiating event for Gwl activation is phosphorylation by MPF of the proline-directed sites T193 and T206 in the presumptive activation loop. After this priming step, Gwl can intramolecularly phosphorylate its C-terminal tail at pS883; this site probably plays a role similar to that of the tail/Z motif of other AGC kinases. These events largely (but not completely) explain the full activation of Gwl at M phase.

A tunable zinc finger-based framework for Boolean logic computation in mammalian cells

The ability to perform molecular-level computation in mammalian cells has the potential to enable a new wave of sophisticated cell-based therapies and diagnostics. To this end, we developed a Boolean logic framework utilizing artificial Cys₂–His₂ zinc finger transcription factors (ZF-TFs) as computing elements. Artificial ZFs can be designed to specifically bind different DNA sequences and thus comprise a diverse set of components ideal for the construction of scalable networks. We generate ZF-TF activators and repressors and demonstrate a novel, general method to tune ZF-TF response by fusing ZF-TFs to leucine zipper homodimerization domains. We describe 15 transcriptional activators that display 2- to 463-fold induction and 15 transcriptional repressors that show 1.3- to 16-fold repression. Using these ZF-TFs, we compute OR, NOR, AND and NAND logic, employing hybrid promoters and split intein-mediated protein splicing to integrate signals. The split intein strategy is able to fully reconstitute the ZF-TFs, maintaining them as a uniform set of computing elements. Together, these components comprise a robust platform for building mammalian synthetic gene circuits capable of precisely modulating cellular behavior.

Dermacentor variabilis: characterization and modeling of macrophage migration inhibitory factor with phylogenetic comparisons to other ticks, insects and parasitic nematodes

We have identified and characterized the full length cDNA sequence of macrophage migration inhibitory factor (MIF) from the American dog tick, Dermacentor variabilis. The nucleotide and putative amino acid sequences from this study shared a high level of sequence conservation with other tick MIFs. The bioinformatics analysis showed across species conservation of the MIF amino acid sequence in ticks, insects and nematodes. The multiple sequence alignment identified Pro 1, 3, 55; Thr 7, 112; Asn 8, 72; Ile 64, 96; Gly 65, 110, Ser 63 and Leu 87 amino acids to be highly conserved among the sequences selected for this study. Tick MIF does not have the oxidoreductase domain as found in MIFs from other animals suggesting that tick MIF is not capable of performing as an oxidoreductase. The phylogenetic analysis revealed that tick MIFs share a closer evolutionary proximity to parasitic nematode MIFs than to insect MIFs.

C-GRAAL: common-neighbors-based global GRAph ALignment of biological networks

Networks are an invaluable framework for modeling biological systems. Analyzing protein-protein interaction (PPI) networks can provide insight into underlying cellular processes. It is expected that comparison and alignment of biological networks will have a similar impact on our understanding of evolution, biological function, and disease as did sequence comparison and alignment. Here, we introduce a novel pairwise global alignment algorithm called Common-neighbors based GRAph ALigner (C-GRAAL) that uses heuristics for maximizing the number of aligned edges between two networks and is based solely on network topology. As such, it can be applied to any type of network, such as social, transportation, or electrical networks. We apply C-GRAAL to align PPI networks of eukaryotic and prokaryotic species, as well as inter-species PPI networks, and we demonstrate that the resulting alignments expose large connected and functionally topologically aligned regions. We use the resulting alignments to transfer biological knowledge across species, successfully validating many of the predictions. Moreover, we show that C-GRAAL can be used to align human-pathogen inter-species PPI networks and that it can identify patterns of pathogen interactions with host proteins solely from network topology.

Exploiting BAC-end sequences for the mining, characterization and utility of new short sequences repeat (SSR) markers in Citrus

The aim of this study was to develop a large set of microsatellite markers based on publicly available BAC-end sequences (BESs), and to evaluate their transferability, discriminating capacity of genotypes and mapping ability in Citrus. A set of 1,281 simple sequence repeat (SSR) markers were developed from the 46,339 Citrus clementina BAC-end sequences (BES), of them 20.67% contained SSR longer than 20 bp, corresponding to roughly one perfect SSR per 2.04 kb. The most abundant motifs were di-nucleotide (16.82%) repeats. Among all repeat motifs (TA/AT)n is the most abundant (8.38%), followed by (AG/CT)n (4.51%). Most of the BES-SSR are located in the non-coding region, but 1.3% of BES-SSRs were found to be associated with transposable element (TE). A total of 400 novel SSR primer pairs were synthesized and their transferability and polymorphism tested on a set of 16 Citrus and Citrus relative's species. Among these 333 (83.25%) were successfully amplified and 260 (65.00%) showed cross-species transferability with Poncirus trifoliata and Fortunella sp. These cross-species transferable markers could be useful for cultivar identification, for genomic study of Citrus, Poncirus and Fortunella sp. Utility of the developed SSR marker was demonstrated by identifying a set of 118 markers each for construction of linkage map of Citrus reticulata and Poncirus trifoliata. Genetic diversity and phylogenetic relationship among 40 Citrus and its related species were conducted with the aid of 25 randomly selected SSR primer pairs and results revealed that citrus genomic SSRs are superior to genic SSR for genetic diversity and germplasm characterization of Citrus spp.

MAPI: towards the integrated exploitation of bioinformatics Web Services

Background

Bioinformatics is commonly featured as a well assorted list of available web resources. Although diversity of services is positive in general, the proliferation of tools, their dispersion and heterogeneity complicate the integrated exploitation of such data processing capacity.

Results

To facilitate the construction of software clients and make integrated use of this variety of tools, we present a modular programmatic application interface (MAPI) that provides the necessary functionality for uniform representation of Web Services metadata descriptors including their management and invocation protocols of the services which they represent. This document describes the main functionality of the framework and how it can be used to facilitate the deployment of new software under a unified structure of bioinformatics Web Services. A notable feature of MAPI is the modular organization of the functionality into different modules associated with specific tasks. This means that only the modules needed for the client have to be installed, and that the module functionality can be extended without the need for re-writing the software client.

Conclusions

The potential utility and versatility of the software library has been demonstrated by the implementation of several currently available clients that cover different aspects of integrated data processing, ranging from service discovery to service invocation with advanced features such as workflows composition and asynchronous services calls to multiple types of Web Services including those registered in repositories (e.g. GRID-based, SOAP, BioMOBY, R-bioconductor, and others).

Paraoxonase 1, quorum sensing, and P. aeruginosa infection: a novel model

Pseudomonas aeruginosa is a Gram-negative bacterium which exacts a heavy burden on immunocompromised patients, but is non-pathogenic in a healthy host. Using small signaling molecules called acyl-homoserine lactones (AHLs), populations of P. aeruginosa can coordinate phenotypic changes, including biofilm formation and virulence factor secretion. This concentration-dependent process is called quorum sensing (QS). Interference with QS has been identified as a potential source of new treatments for P. aeruginosa infection. The human enzyme paraoxonase 1 (PON1) degrades AHL molecules, and is a promising candidate for QS interference therapy. Although paraoxonase orthologs exist in many species, genetic redundancy in humans and other mammals has made studying the specific effects of PON1 quite difficult. Arthropods, however, do not express any PON homologs. We generated a novel model to study the specific effects of PON1 by transgenically expressing human PON1 in Drosophila melanogaster. Using this model, we showed that P. aeruginosa infection lethality is QS-dependent, and that expression of PON1 has a protective effect. This work demonstrates the value of a D. melanogaster model for investigating the specific functions of members of the paraoxonase family in vivo, and suggests that PON1 plays a role in innate immunity.

Sequence context for transcription and translation of the Arabidopsis RPL23aA and RPL23aB paralogs

The 80S cytoplasmic ribosome is responsible for translating the transcriptome into the proteome. Demand for ribosome production depends on growth rate, and both the ribosomal RNA (rRNA) and ribosomal protein (RP) components must respond coordinately and rapidly to positive and negative growth stimuli to prevent deleterious effects of excess or insufficient subunits. The 81 RPs of the Arabidopsis 80S ribosome are encoded by multigene families that often exhibit overlapping patterns of transcript accumulation; however, only one isoform of each RP family (with the exception of a small number of acidic RPs) assembles into a single ribosome. Here we dissected the regulatory regions (RRs) of both members of the RPL23a family (RPL23aA and RPL23aB) to identify salient cis-acting elements involved in transcriptional, posttranscriptional, and translational regulation of expression. Full length and truncated RRs of RPL23a paralogs were cloned upstream of a GUS reporter gene and expressed in Arabidopsis transgenic plants. High level expression in mitotically active tissues, driven by RPL23aA and RPL23aB RRs, required TATA-box, telo-box, and site II motif elements. First and second introns were found to play a minor role in posttranscriptional regulation of paralogs, and conserved transcript features (e.g., UTR base composition) may be involved in enhancing translational efficiency. Overall, our results indicate that RPL23a expression is governed by a complex network of multiple regulatory layers.

Generation, annotation and analysis of first large-scale expressed sequence tags from developing fiber of Gossypium barbadense L

Background

Cotton fiber is the world's leading natural fiber used in the manufacture of textiles. Gossypium is also the model plant in the study of polyploidization, evolution, cell elongation, cell wall development, and cellulose biosynthesis. G. barbadense L. is an ideal candidate for providing new genetic variations useful to improve fiber quality for its superior properties. However, little is known about fiber development mechanisms of G. barbadense and only a few molecular resources are available in GenBank.

Methodology and Principal Findings

In total, 10,979 high-quality expressed sequence tags (ESTs) were generated from a normalized fiber cDNA library of G. barbadense. The ESTs were clustered and assembled into 5852 unigenes, consisting of 1492 contigs and 4360 singletons. The blastx result showed 2165 unigenes with significant similarity to known genes and 2687 unigenes with significant similarity to genes of predicted proteins. Functional classification revealed that unigenes were abundant in the functions of binding, catalytic activity, and metabolic pathways of carbohydrate, amino acid, energy, and lipids. The function motif/domain-related cytoskeleton and redox homeostasis were enriched. Among the 5852 unigenes, 282 and 736 unigenes were identified as potential cell wall biosynthesis and transcription factors, respectively. Furthermore, the relationships among cotton species or between cotton and other model plant systems were analyzed. Some putative species-specific unigenes of G. barbadense were highlighted.

Conclusions/Significance

The ESTs generated in this study are from the first large-scale EST project for G. barbadense and significantly enhance the number of G. barbadense ESTs in public databases. This knowledge will contribute to cotton improvements by studying fiber development mechanisms of G. barbadense, establishing a breeding program using marker-assisted selection, and discovering candidate genes related to important agronomic traits of cotton through oligonucleotide array. Our work will also provide important resources for comparative genomics, polyploidization, and genome evolution among Gossypium species.

In silico gene expression analysis in Codonopsis lanceolata root

Expressed sequence tags (ESTs) provide valuable tools that can be used to predict the genes involved in primary and secondary metabolite synthesis. To the best of our knowledge, ESTs have not yet been developed for Codonopsis. lanceolata, and therefore, the EST referenced in this report is the first transcript for C. lanceolata. A cDNA library was constructed using the roots of C. lanceolata plants that were grown in a field. The selected 881 cDNA clones were sequenced and processed with an EST pipeline, resulting in 636 unique sequences, including 517 singletons and 119 contig sequences. Using bioinformatics tools, 81% of the EST sequence was putatively annotated. Data for unique transcripts were mined from biological databases and functionally classified using gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes Orthology, KEGG pathway maps, and protein family. The GO-based analyses were examined in terms of biotic and abiotic stress response, transport, cellular component organization, biogenesis, and secondary metabolic processes. The KEGG-based analyses of most transcripts were sorted by carbohydrate metabolism, energy metabolism, and biosynthesis of secondary metabolites. Five randomly-selected putative genes were used for an expression study using various stresses such as salt, H(2)O(2), salicylic acid, and methyl jasmonic acid. Mined data were organized in "The Codonopsis EST Database" (www.bioherbs.khu.ac.kr/Codonopsis).

Exploring the factors determining the dynamics of different protein folds

Normal mode analyses of homologous proteins at the family and superfamily level show that slow dynamics are similar and are preserved through evolution. This study investigates how the slow dynamics of proteins is affected by variation in the protein architecture and fold. For this purpose, we have used computer-generated protein models based on idealized protein structures with varying folds. These are shown to be protein-like in their behavior, and they are used to investigate the influence of architecture and fold on the slow dynamics. We compared the dynamics of models having different folds but similar architecture and found the architecture to be the dominant factor for the slow dynamics.

Gene discovery in Triatoma infestans

BACKGROUND

Triatoma infestans is the most relevant vector of Chagas disease in the southern cone of South America. Since its genome has not yet been studied, sequencing of Expressed Sequence Tags (ESTs) is one of the most powerful tools for efficiently identifying large numbers of expressed genes in this insect vector.

RESULTS

In this work, we generated 826 ESTs, resulting in an increase of 47% in the number of ESTs available for T. infestans. These ESTs were assembled in 471 unique sequences, 151 of which represent 136 new genes for the Reduviidae family.

CONCLUSIONS

Among the putative new genes for the Reduviidae family, we identified and described an interesting subset of genes involved in development and reproduction, which constitute potential targets for insecticide development.

OrchidBase: a collection of sequences of the transcriptome derived from orchids

Orchids are one of the most ecological and evolutionarily significant plants, and the Orchidaceae is one of the most abundant families of the angiosperms. Genetic databases will be useful not only for gene discovery but also for future genomic annotation. For this purpose, OrchidBase was established from 37,979,342 sequence reads collected from 11 in-house Phalaenopsis orchid cDNA libraries. Among them, 41,310 expressed sequence tags (ESTs) were obtained by using Sanger sequencing, whereas 37,908,032 reads were obtained by using next-generation sequencing (NGS) including both Roche 454 and Solexa Illumina sequencers. These reads were assembled into 8,501 contigs and 76,116 singletons, resulting in 84,617 non-redundant transcribed sequences with an average length of 459 bp. The analysis pipeline of the database is an automated system written in Perl and C#, and consists of the following components: automatic pre-processing of EST reads, assembly of raw sequences, annotation of the assembled sequences and storage of the analyzed information in SQL databases. A web application was implemented with HTML and a Microsoft .NET Framework C# program for browsing and querying the database, creating dynamic web pages on the client side, analyzing gene ontology (GO) and mapping annotated enzymes to KEGG pathways. The online resources for putative annotation can be searched either by text or by using BLAST, and the results can be explored on the website and downloaded. Consequently, the establishment of OrchidBase will provide researchers with a high-quality genetic resource for data mining and facilitate efficient experimental studies on orchid biology and biotechnology. The OrchidBase database is freely available at http://lab.fhes.tn.edu.tw/est.

Isolation and characterization of Acanthamoeba spp. from air-conditioners in Kuala Lumpur, Malaysia

During a study on the quality of the indoor environment, Acanthamoeba spp. were detected in 20 out of 87 dust samples collected from air-conditioners installed in a four-story campus building located in Kuala Lumpur, Malaysia. Twenty-one cloned Acanthamoeba isolates designated as IMU1 to IMU21 were established from the positive primary cultures. Five species were identified from the 16 isolates according to the morphological criteria of Pussard and Pons; i.e. A. castellanii, A. culbertsoni, A. griffini, A. hatchetti and A. polyphaga. Species identities for the remaining five isolates (IMU4, IMU5, IMU15, IMU20 and IMU21), however, could not be determined morphologically. At genotypic characterization, these isolates were placed into T3 (IMU14); T5 (IMU16 and IMU17) and T4 (all the remaining isolates). To predict the potential pathogenicity of these Acanthamoeba isolates, thermo- and osmotolerance tests were employed; many isolates were predicted as potential human pathogens based on the outcome of these tests. This is the first time potentially pathogenic Acanthamoeba have been isolated from air-conditioners in Malaysia.

A tandem SH2 domain in transcription elongation factor Spt6 binds the phosphorylated RNA polymerase II C-terminal repeat domain (CTD)

Spt6 is an essential transcription elongation factor and histone chaperone that binds the C-terminal repeat domain (CTD) of RNA polymerase II. We show here that Spt6 contains a tandem SH2 domain with a novel structure and CTD-binding mode. The tandem SH2 domain binds to a serine 2-phosphorylated CTD peptide in vitro, whereas its N-terminal SH2 subdomain, which we previously characterized, does not. CTD binding requires a positively charged crevice in the C-terminal SH2 subdomain, which lacks the canonical phospho-binding pocket of SH2 domains and had previously escaped detection. The tandem SH2 domain is apparently required for transcription elongation in vivo as its deletion in cells is lethal in the presence of 6-azauracil.

Topological network alignment uncovers biological function and phylogeny

Sequence comparison and alignment has had an enormous impact on our understanding of evolution, biology and disease. Comparison and alignment of biological networks will probably have a similar impact. Existing network alignments use information external to the networks, such as sequence, because no good algorithm for purely topological alignment has yet been devised. In this paper, we present a novel algorithm based solely on network topology, that can be used to align any two networks. We apply it to biological networks to produce by far the most complete topological alignments of biological networks to date. We demonstrate that both species phylogeny and detailed biological function of individual proteins can be extracted from our alignments. Topology-based alignments have the potential to provide a completely new, independent source of phylogenetic information. Our alignment of the protein-protein interaction networks of two very different species-yeast and human-indicate that even distant species share a surprising amount of network topology, suggesting broad similarities in internal cellular wiring across all life on Earth.