Clustal W and Clustal X version 2.0

Hylin a1, the first cytolytic peptide isolated from the arboreal South American frog Hypsiboas albopunctatus ("spotted treefrog")

RP-HPLC fractionation of the electrically stimulated skin secretion of the arboreal South American frog Hypsiboas albopunctatus ("spotted treefrog") led to the isolation of a cytolytic C-terminally amidated peptide. This novel peptide, named hylin a1 (Hy-a1), consists of 18 amino acid residues (IFGAILPLALGALKNLIK-NH(2)). The alpha-helical structure of the synthetic hylin a1 peptide was confirmed by CD spectroscopy in the presence of 60% (v/v) TFE. The synthetic peptide displayed broad-spectrum antimicrobial activity against Gram-negative and Gram-positive bacteria including Escherichia coli, Staphylococcus aureus, Enterococcus faecalis, Bacillus subtilis and Pseudomonas aeruginosa and also against fungi (Candida albicans, C. krusei, C. parapsilosis and Cryptococcus neoformans). Hylin a1 was also able to disrupt human erytrocytes (HC(50)=18 microM). Similarity analysis using PSI-BLAST revealed 50-44% of identity to maximins Hv, H16, H15 and H10 from Bombina maxima and also to hylins b1 and b2 (Hy-b1 and Hy-b2) from Hypsiboas lundii (synonym: Hyla biobeba).

Acetylcholinesterase genes in the basal Hexapod Orchesella villosa

Acetylcholinesterase (AChE) is a key enzyme of the cholinergic nerve system. Of the two forms found in insects, the predominant one is active in the synapses and is the target of organophosphate and carbamate insecticides, while the role of the second is currently unknown. Two acetylcholinesterase cDNAs from the basal hexapod Orchesella villosa have been characterized and compared with others reported form insects. One form conforms well to the typical structure, while the other is characterized by an unusual 3' region. No amino acid mutation could be directly associated with known resistance mutations in other insect species or to a clear signal of selection in the distribution of alleles, although the action of some population process is suggested.

Form, function and phylogenetics of NITRs in bony fish

Novel immune-type receptors (NITRs) are encoded by clusters of multigene families and have been identified in multiple bony fish species. All NITRs possess one extracellular immunoglobulin (Ig) domain of the variable (V) type and recent crystal structures of NITR V domains demonstrate their high degree of similarity to V domains of antigen receptors. Many NITRs possess a second extracellular Ig domain of the intermediate (I) type which helps differentiate NITRs from other V domain receptors. The majority of NITRs are type I transmembrane receptors; however, a small number are predicted to encode secreted proteins. Based on their sequence and structure, NITRs have been proposed to be "functional orthologs" of mammalian natural killer receptors (NKRs). Like NKRs, most NITRs possess short functional motifs permitting their classification as inhibitory or activating. NITRs lacking these motifs are functionally ambiguous. Inhibitory and activating NKRs utilize opposing signaling mechanisms to influence the response of NK cells to target cells; studies employing recombinant NITRs suggest that these signaling pathways are conserved between NKRs and NITRs. An analysis of all published NITR sequences demonstrates the conserved nature of multiple residues within the NITR Ig domains permitting the identification of NITR ESTs from salmon, cod, halibut, lake whitefish and stickleback species. Complete data sets of NITRs from the sequencing of the zebrafish and medaka genomes provide insight into the evolution of the NITRs within bony fish species. It is likely that all teleost species encode NITRs which function within innate immunity to regulate cell mediated cytotoxicity.

A missense mutation in CASK causes FG syndrome in an Italian family

First described in 1974, FG syndrome (FGS) is an X-linked multiple congenital anomaly/mental retardation (MCA/MR) disorder, characterized by high clinical variability and genetic heterogeneity. Five loci (FGS1-5) have so far been linked to this phenotype on the X chromosome, but only one gene, MED12, has been identified to date. Mutations in this gene account for a restricted number of FGS patients with a more distinctive phenotype, referred to as the Opitz-Kaveggia phenotype. We report here that a p.R28L (c.83G-->T) missense mutation in CASK causes FGS phenotype in an Italian family previously mapped to Xp11.4-p11.3 (FGS4). The identified missense mutation cosegregates with the phenotype in this family and is absent in 1000 control X chromosomes of the same ethnic origin. An extensive analysis of CASK protein functions as well as structural and dynamic studies performed by molecular dynamics (MD) simulation did not reveal significant alterations induced by the p.R28L substitution. However, we observed a partial skipping of the exon 2 of CASK, presumably a consequence of improper recognition of exonic splicing enhancers (ESEs) induced by the c.83G-->T transversion. CASK is a multidomain scaffold protein highly expressed in the central nervous system (CNS) with specific localization to the synapses, where it forms large signaling complexes regulating neurotransmission. We suggest that the observed phenotype is most likely a consequence of an altered CASK expression profile during embryogenesis, brain development, and differentiation.

Neisseria meningitidis antigen NMB0088: sequence variability, protein topology and vaccine potential

The significance of Neisseria meningitidis serogroup B membrane proteins as vaccine candidates is continually growing. Here, we studied different aspects of antigen NMB0088, a protein that is abundant in outer-membrane vesicle preparations and is thought to be a surface protein. The gene encoding protein NMB0088 was sequenced in a panel of 34 different meningococcal strains with clinical and epidemiological relevance. After this analysis, four variants of NMB0088 were identified; the variability was confined to three specific segments, designated VR1, VR2 and VR3. Secondary structure predictions, refined with alignment analysis and homology modelling using FadL of Escherichia coli, revealed that almost all the variable regions were located in extracellular loop domains. In addition, the NMB0088 antigen was expressed in E. coli and a procedure for obtaining purified recombinant NMB0088 is described. The humoral immune response elicited in BALB/c mice was measured by ELISA and Western blotting, while the functional activity of these antibodies was determined in a serum bactericidal assay and an animal protection model. After immunization in mice, the recombinant protein was capable of inducing a protective response when it was administered inserted into liposomes. According to our results, the recombinant NMB0088 protein may represent a novel antigen for a vaccine against meningococcal disease. However, results from the variability study should be considered for designing a cross-protective formulation in future studies.

Development of a novel real-time reverse-transcriptase PCR method for the detection of H275Y positive influenza A H1N1 isolates

During the 2007–2008 influenza season global strain surveillance for antiviral resistance revealed the sudden emergence of oseltamivir resistance in influenza A H1N1 isolates. Although oseltamivir resistance rates vary from region to region, 16% of isolates tested globally were found to be oseltamivir resistant by a histidine to tyrosine mutation of residue 275 of the neuraminidase gene of influenza A. In order to implement effective resistance testing locally a novel real-time reverse-transcriptase PCR (RT-PCR) assay was developed for the detection of the H275Y mutation. To evaluate this method, 40 oseltamivir resistant and 61 oseltamivir sensitive H1N1 influenza isolates were tested using Sanger sequencing, which is the reference method for detection of resistance, pyrosequencing and the novel H275Y RT-PCR assay. In comparison to Sanger sequencing, the sensitivity and specificity of the H275Y RT-PCR assay were 100% (40/40) and 100% (61/61) respectively, while the sensitivity and specificity of pyrosequencing were 100% (40/40) and 97.5% (60/61) respectively. Although all three methods were effective in detecting the H275Y mutation associated with oseltamivir resistance, the H275Y RT-PCR assay was the most rapid and could easily be incorporated into an influenza subtyping protocol.

Cold adaptation of zinc metalloproteases in the thermolysin family from deep sea and arctic sea ice bacteria revealed by catalytic and structural properties and molecular dynamics: new insights into relationship between conformational flexibility and hydrogen bonding

Increased conformational flexibility is the prevailing explanation for the high catalytic efficiency of cold-adapted enzymes at low temperatures. However, less is known about the structural determinants of flexibility. We reported two novel cold-adapted zinc metalloproteases in the thermolysin family, vibriolysin MCP-02 from a deep sea bacterium and vibriolysin E495 from an Arctic sea ice bacterium, and compared them with their mesophilic homolog, pseudolysin from a terrestrial bacterium. Their catalytic efficiencies, k(cat)/K(m) (10-40 degrees C), followed the order pseudolysin < MCP-02 < E495 with a ratio of approximately 1:2:4. MCP-02 and E495 have the same optimal temperature (T(opt), 57 degrees C, 5 degrees C lower than pseudolysin) and apparent melting temperature (T(m) = 64 degrees C, approximately 10 degrees C lower than pseudolysin). Structural analysis showed that the slightly lower stabilities resulted from a decrease in the number of salt bridges. Fluorescence quenching experiments and molecular dynamics simulations showed that the flexibilities of the proteins were pseudolysin < MCP-02 < E495, suggesting that optimization of flexibility is a strategy for cold adaptation. Molecular dynamics results showed that the ordinal increase in flexibility from pseudolysin to MCP-02 and E495, especially the increase from MCP-02 to E495, mainly resulted from the decrease of hydrogen-bond stability in the dynamic structure, which was due to the increase in asparagine, serine, and threonine residues. Finally, a model for the cold adaptation of MCP-02 and E495 was proposed. This is the first report of the optimization of hydrogen-bonding dynamics as a strategy for cold adaptation and provides new insights into the structural basis underlying conformational flexibility.

Complete genome of the broad-host-range Erwinia amylovora phage phiEa21-4 and its relationship to Salmonella phage felix O1

The first complete genome sequence for a myoviridal bacteriophage, PhiEa21-4, infecting Erwinia amylovora, Erwinia pyrifoliae, and Pantoea agglomerans strains has been determined. The unique sequence of this terminally redundant, circularly permuted genome is 84,576 bp. The PhiEa21-4 genome has a GC content of 43.8% and contains 117 putative protein-coding genes and 26 tRNA genes. PhiEa21-4 is the first phage in which a precisely conserved rho-independent terminator has been found dispersed throughout the genome, with 24 copies in all. Also notable in the PhiEa21-4 genome are the presence of tRNAs with six- and nine-base anticodon loops, the absence of a small packaging terminase subunit, and the presence of nadV, a principle component of the NAD(+) salvage pathway, which has been found in only a few phage genomes to date. PhiEa21-4 is the first reported Felix O1-like phage genome; 56% of the predicted PhiEa21-4 proteins share homology with those of the Salmonella phage. Apart from this similarity to Felix O1, the PhiEa21-4 genome appears to be substantially different, both globally and locally, from previously reported sequences. A total of 43 of the 117 genes are unique to PhiEa21-4, and 32 of the Felix O1-like genes do not appear in any phage genome sequences other than PhiEa21-4 and Felix O1. N-terminal sequencing and matrix-assisted laser desorption ionization-time of flight analysis resulted in the identification of five PhiEa21-4 genes coding for virion structural proteins, including the major capsid protein.

Crystal structure of chlorite dismutase, a detoxifying enzyme producing molecular oxygen

Chlorite dismutase (Cld) is a key enzyme of perchlorate and chlorate respiration. This heme-based protein reduces the toxic compound chlorite into the innocuous chloride anion in a very efficient way while producing molecular oxygen. A sequence comparison between Cld homologues shows a highly conserved family. The crystal structure of Azospira oryzae strain GR-1 Cld is reported to 2.1 A resolution. The structure reveals a hexameric organization of the Cld, while each monomer exhibits a ferredoxin-like fold. The six subunits are organized in a ring structure with a maximal diameter of 9 nm and an inner diameter of 2 nm. The heme active-site pocket is solvent accessible both from the inside and the outside of the ring. Moreover, a second anion binding site that could accommodate the assumed reaction intermediate ClO(-) for further transformation has been identified near the active site. The environment of the heme cofactor was investigated with electron paramagnetic resonance spectroscopy. Apart from the high-spin ferric signal of the five-coordinate resting-state enzyme, two low-spin signals were found corresponding to six-coordinate species. The current crystal structure confirms and complements a recently proposed catalytic mechanism that proceeds via a ferryl species and a ClO(-) anion. Our structural data exclude cooperativity between the iron centers.

Molecular cloning and characterization of a candidate human growth-related and time-keeping constitutive cell surface hydroquinone (NADH) oxidase

ENOX (ECTO-NOX) proteins are growth-related cell surface proteins that catalyze both hydroquinone or NADH oxidation and protein disulfide-thiol interchange and exhibit both prion-like and time-keeping (clock) properties. The two enzymatic activities they catalyze alternate to generate a regular period of 24 min in length. Here we report the cloning, expression, and characterization of a human candidate constitutive ENOX (CNOX or ENOX1) protein. The gene encoding this 643 amino acid long protein is located on chromosome 13 (13q 14.11). Functional motifs previously identified by site-directed mutagenesis in a cancer-associated ENOX (tNOX or ENOX2) as adenine nucleotide or copper binding along with essential cysteines are present, but the drug-binding motif (EEMTE) sequence of ENOX2 is absent. The activities of the recombinant protein expressed in Escherichia coli were not affected by capsaicin, EGCg, and other ENOX2-inhibiting substances. The purified recombinant protein bound ca. 2 mol of copper/mol of protein. Bound copper was necessary for activity. H260 and H579 were required for copper binding as confirmed by site-directed mutagenesis, loss of copper-binding capacity, and resultant loss of enzymatic activity. Addition of melatonin phased the 24 min period such that the next complete period began exactly 24 min after the melatonin addition as appears to be characteristic of ENOX1 activities in general. Oxidative activity was exhibited with both NAD(P)H and reduced coenzyme Q as substrate. Concentrated solutions of the purified candidate ENOX1 protein irreversibly formed insoluble aggregates, devoid of enzymatic activity, resembling amyloid.

A PAS domain with an oxygen labile [4Fe-4S](2+) cluster in the oxygen sensor kinase NreB of Staphylococcus carnosus

The cytoplasmic histidine sensor kinase NreB of Staphylococcus carnosus responds to O(2) and controls together with the response regulator NreC the expression of genes of nitrate/nitrite respiration. nreBC homologous genes were found in Staphylococcus strains and Bacillus clausii, and a modified form was found in some Lactobacillus strains. NreB contains a sensory domain with similarity to heme B binding PAS domains. Anaerobically prepared NreB of S. carnosus exhibited a (diamagnetic) [4Fe-4S](2+) cluster when assessed by Mossbauer spectroscopy. Upon reaction with air, the cluster was degraded with a half-life of approximately 2.5 min. No significant amounts of Mossbauer or EPR detectable intermediates were found during the decay, but magnetic Mossbauer spectra revealed formation of diamagnetic [2Fe-2S](2+) clusters. After extended exposure to air, NreB was devoid of a FeS cluster. Photoreduction with deazaflavin produced small amounts of [4Fe-4S](+), which were degraded subsequently. The magnetically perturbed Mossbauer spectrum of the [4Fe-4S](2+) cluster corroborated the S = 0 spin state and revealed uniform electric field gradient tensors of the iron sites, suggesting full delocalization of the valence electrons and binding of each of the Fe ions by four S ligands, including the ligand to the protein. Mutation of each of the four Cys residues inactivated NreB function in vivo in accordance with their role as ligands. [4Fe-4S](2+) cluster-containing NreB had high kinase activity. Exposure to air decreased the kinase activity and content of the [4Fe-4S](2+) cluster with similar half-lives. We conclude that the sensory domain of NreB represents a new type of PAS domain containing a [4Fe-4S](2+) cluster for sensing and function.

The A312L 5'-UTR of Chlorella virus PBCV-1 is a translational enhancer in Arabidopsis thaliana

PBCV-1 (Paramecium bursaria Chlorella virus) is a large double stranded DNA virus that replicates in certain eukaryotic chlorella like green algae. The PBCV-1 A312L gene encodes a 33-kDa protein whose function currently is unknown. The 5'-UTR of the A312L mRNA is 153 nucleotides, longer than the 5'-UTR in any other PBCV-1 gene. The sequence 5'-AAAC was repeated 17 times within 156bp 5' to the A312L gene start codon and this sequence was repeated 13 times continuously in the 5'-UTR of the mRNA. Recombinant genes were constructed in vector pBI121 that contained the A312L 5'-UTR, in both the forward and inverse-complement orientations, fused to the GUS gene under the control of the CaMV 35S promoter. These constructs were introduced into Arabidopsis thaliana and the results indicated that the A312L 5'-UTR functions as a translational enhancer only in the forward orientation. Overall, the ratio of GUS enzyme activity to GUS mRNA was 15-fold higher in constructs derived from the A312L 5'-UTR in the forward orientation as compared to constructs containing the 5'-UTR in the inverse-complement orientation or those lacking the A312L 5'-UTR.

The haloarchaeal chromosome replication machinery

The powerful combination of genetic and biochemical analysis has provided many key insights into the structure and function of the chromosomal DNA replication machineries of bacterial and eukaryotic cells. In contrast, in the archaea, biochemical studies have dominated, mainly due to the absence of efficient genetic systems for these organisms. This situation is changing, however, and, in this regard, the genetically tractable haloarchaea Haloferax volcanii and Halobacterium sp. NRC-1 are emerging as key models. In the present review, I give an overview of the components of the replication machinery in the haloarchaea, with particular emphasis on the protein factors presumed to travel with the replication fork.

PolV(PolIVb) function in RNA-directed DNA methylation requires the conserved active site and an additional plant-specific subunit

Two forms of a plant-specific RNA polymerase (Pol), PolIV(PolIVa) and PolV(PolIVb), currently defined by their respective largest subunits [NRPD1(NRPD1a) and NRPE1(NRPD1b)], have been implicated in the production and activity of 24-nt small RNAs (sRNAs) in RNA-directed DNA methylation (RdDM). Prevailing models support the view that PolIV(PolIVa) plays an upstream role in RdDM by producing the 24-nt sRNAs, whereas PolV(PolIVb) would act downstream at a structural rather than an enzymatic level to reinforce sRNA production by PolIV(PolIVa) and mediate DNA methylation. However, the composition and mechanism of action of PolIV(PolIVa)/PolV(PolIVb) remain unclear. In this work, we have identified a plant-specific PolV(PolIVb) subunit, NRPE5a, homologous to NRPB5a, a common subunit shared by PolI-III and shown here to be present in PolIV(PolIVa). Our results confirm the combinatorial diversity of PolIV(PolIVa)/PolV(PolIVb) subunit composition and indicate that these plant-specific Pols are eukaryotic-type polymerases. Moreover, we show that nrpe5a-1 mutation differentially impacts sRNAs accumulation at various PolIV(PolIVa)/PolV(PolIVb)-dependent loci, indicating a target-specific requirement for NRPE5a in the process of PolV(PolIVb)-dependent gene silencing. We then describe that the triad aspartate motif present in the catalytic center of PolV(PolIVb) is required for recapitulation of all activities associated with this Pol complex in RdDM, suggesting that RNA polymerization is important for PolV(PolIVb) to perform its regulatory functions.

Rice Pi5-mediated resistance to Magnaporthe oryzae requires the presence of two coiled-coil-nucleotide-binding-leucine-rich repeat genes

Rice blast, caused by the fungus Magnaporthe oryzae, is one of the most devastating diseases of rice. To understand the molecular basis of Pi5-mediated resistance to M. oryzae, we cloned the resistance (R) gene at this locus using a map-based cloning strategy. Genetic and phenotypic analyses of 2014 F2 progeny from a mapping population derived from a cross between IR50, a susceptible rice cultivar, and the RIL260 line carrying Pi5 enabled us to narrow down the Pi5 locus to a 130-kb interval. Sequence analysis of this genomic region identified two candidate genes, Pi5-1 and Pi5-2, which encode proteins carrying three motifs characteristic of R genes: an N-terminal coiled-coil (CC) motif, a nucleotide-binding (NB) domain, and a leucine-rich repeat (LRR) motif. In genetic transformation experiments of a susceptible rice cultivar, neither the Pi5-1 nor the Pi5-2 gene was found to confer resistance to M. oryzae. In contrast, transgenic rice plants expressing both of these genes, generated by crossing transgenic lines carrying each gene individually, conferred Pi5-mediated resistance to M. oryzae. Gene expression analysis revealed that Pi5-1 transcripts accumulate after pathogen challenge, whereas the Pi5-2 gene is constitutively expressed. These results indicate that the presence of these two genes is required for rice Pi5-mediated resistance to M. oryzae.

Development and application of versatile high density microarrays for genome-wide analysis of Streptomyces coelicolor: characterization of the HspR regulon

BACKGROUND

DNA microarrays are a key resource for global analysis of genome content, gene expression and the distribution of transcription factor binding sites. We describe the development and application of versatile high density ink-jet in situ-synthesized DNA arrays for the G+C rich bacterium Streptomyces coelicolor. High G+C content DNA probes often perform poorly on arrays, yielding either weak hybridization or non-specific signals. Thus, more than one million 60-mer oligonucleotide probes were experimentally tested for sensitivity and specificity to enable selection of optimal probe sets for the genome microarrays. The heat-shock HspR regulatory system of S. coelicolor, a well-characterized repressor with a small number of known targets, was exploited to test and validate the arrays for use in global chromatin immunoprecipitation-on-chip (ChIP-chip) and gene expression analysis.

RESULTS

In addition to confirming dnaK, clpB and lon as in vivo targets of HspR, it was revealed, using a novel ChIP-chip data clustering method, that HspR also apparently interacts with ribosomal RNA (rrnD operon) and specific transfer RNA genes (the tRNAGln/tRNAGlu cluster). It is suggested that enhanced synthesis of Glu-tRNAGlu may reflect increased demand for tetrapyrrole biosynthesis following heat-shock. Moreover, it was found that heat-shock-induced genes are significantly enriched for Gln/Glu codons relative to the whole genome, a finding that would be consistent with HspR-mediated control of the tRNA species.

CONCLUSIONS

This study suggests that HspR fulfils a broader, unprecedented role in adaptation to stresses than previously recognized -- influencing expression of key components of the translational apparatus in addition to molecular chaperone and protease-encoding genes. It is envisaged that these experimentally optimized arrays will provide a key resource for systems level studies of Streptomyces biology.

Core structure of the yeast spt4-spt5 complex: a conserved module for regulation of transcription elongation

The Spt4-Spt5 complex is an essential RNA polymerase II elongation factor found in all eukaryotes and important for gene regulation. We report here the crystal structure of Saccharomyces cerevisiae Spt4 bound to the NGN domain of Spt5. This structure reveals that Spt4-Spt5 binding is governed by an acid-dipole interaction between Spt5 and Spt4. Mutations that disrupt this interaction disrupt the complex. Residues forming this pivotal interaction are conserved in the archaeal homologs of Spt4 and Spt5, which we show also form a complex. Even though bacteria lack a Spt4 homolog, the NGN domains of Spt5 and its bacterial homologs are structurally similar. Spt4 is located at a position that may help to maintain the functional conformation of the following KOW domains in Spt5. This structural and evolutionary perspective of the Spt4-Spt5 complex and its homologs suggest that it is an ancient, core component of the transcription elongation machinery.

MYC-induced myeloid leukemogenesis is accelerated by all six members of the antiapoptotic BCL family

Signals that control the fine balance between cell death and cell survival are altered during tumorigenesis. Understanding the mechanisms by which this balance is perturbed, leading to excessive cell survival, is important for designing effective therapies. Proteins belonging to the B-cell lymphoma (BCL) family are known to regulate death responses to apoptotic signals, especially those originating within cells. A subset of BCL family members capable of inhibiting cell death is known to contribute to tumorigenesis; however, it is not known whether all six antiapoptotic BCL family members play a causal role in tumor development. Using a mouse model of MYC-driven leukemia, we showed that, in addition to the well characterized BCL2 and BCLxl (BCL2L1), the other four family members -- BCLw (BCL2L2), BCLb (BCL2L10), BFL1 (BCL2A1) and MCL1 -- also cooperate with MYC to accelerate leukemogenesis. In addition, high levels of each family member are found in either solid human tumors or cell lines derived from human leukemias or lymphomas.

Genomic survey of the non-cultivatable opportunistic human pathogen, Enterocytozoon bieneusi

Enterocytozoon bieneusi is the most common microsporidian associated with human disease, particularly in the immunocompromised population. In the setting of HIV infection, it is associated with diarrhea and wasting syndrome. Like all microsporidia, E. bieneusi is an obligate, intracellular parasite, but unlike others, it is in direct contact with the host cell cytoplasm. Studies of E. bieneusi have been greatly limited due to the absence of genomic data and lack of a robust cultivation system. Here, we present the first large-scale genomic dataset for E. bieneusi. Approximately 3.86 Mb of unique sequence was generated by paired end Sanger sequencing, representing about 64% of the estimated 6 Mb genome. A total of 3,804 genes were identified in E. bieneusi, of which 1,702 encode proteins with assigned functions. Of these, 653 are homologs of Encephalitozoon cuniculi proteins. Only one E. bieneusi protein with assigned function had no E. cuniculi homolog. The shared proteins were, in general, evenly distributed among the functional categories, with the exception of a dearth of genes encoding proteins associated with pathways for fatty acid and core carbon metabolism. Short intergenic regions, high gene density, and shortened protein-coding sequences were observed in the E. bieneusi genome, all traits consistent with genomic compaction. Our findings suggest that E. bieneusi is a likely model for extreme genome reduction and host dependence.

Enterocytozoon bieneusi is a clinically significant pathogen associated with human microsporidiosis, particularly in immunocompromised individuals. E. bieneusi is widespread in mammals, and there is no effective commercial treatment for infection. The pathogen cannot be readily cultivated, and animal models are limited. We therefore undertook a sequence survey and generated the first large-scale genomic dataset for E. bieneusi, which we used to study the organization and structure of its genome and to perform a comparative analysis with Encephalitozoon cuniculi, another microsporidian whose genome has been completely sequenced. The E. bieneusi genome showed many traits associated with genome compaction including high gene density, short intergenic regions, shortened proteins, and few introns. With one exception, all E. bieneusi proteins with assigned functions had E. cuniculi homologs. We found a paucity of genes encoding proteins associated with fatty acid and carbon metabolism. The possibility that these core functions are reduced in an intracellular parasite is intriguing, but because the genome sequence of E. bieneusi is incomplete, we cannot exclude the possibility that additional proteins associated with the various metabolic pathways would be discovered in a completed genome.

Sulfolobus tokodaii ST0053 produces a novel thermostable, NAD-dependent medium-chain alcohol dehydrogenase

An alcohol dehydrogenase (ADH) gene, ST0053, from Sulfolobus tokodaii was expressed in Escherichia coli. The purified recombinant enzyme was an NAD-dependent medium-chain ADH with high thermostability and tolerance of a wide range of pHs. This is the first step in creating an experimental functionality library of 10 genes annotated as ADHs in the S. tokodaii genome.

A novel growth hormone 1 gene-derived probe for Oncorhynchus masou formosanus distinguished from the Oncorhynchus subspecies

Traditional mitochondrial 16S rRNA is commonly used in many species identification studies. However, it is difficult to apply to the phylogenetic studies among the Oncorhynchus subspecies, which is a crucial need for management purposes for Oncorhynchus masou formosanus, Taiwan salmon. In this study, we have developed an improved species identification method for Taiwan salmon distinguished with other Oncorhynchus subspecies tested by exploiting PCR for growth hormone (GH) 1 gene. By comparing DNA sequences for GH1 from 11 species of Oncorhynchus subspecies we designed novel PCR primers that exploit differences between Taiwan salmon and other Oncorhynchus subspecies. Therefore, the technique is an important tool in the management of populations of the endangered land-locked Taiwan salmon preventing from their possible hybrids with other Oncorhynchus subspecies once tested.

Clustering of disease-causing mutations on the domain-domain interfaces of ABCC6

Mutations in ABCC6 are responsible for pseudoxanthoma elasticum (PXE), a rare genetic disease affecting the elastic tissues of the body. ABCC6 encodes a 1503 amino acid long ABC transporter, ABCC6/MRP6. The functional link between the impaired activity of the protein and the disease is not known. We have built a homology model of this transporter, and analyzed the distribution of the known 119 missense PXE-associated mutations within the predicted structure. Significant clustering of the missense mutations has been found at complex domain-domain interfaces: at the transmission interface that involves four intracellular loops and the two ABC domains as well as at the ABC-ABC interacting surfaces. The mutations affecting these regions are 2.75 and 3.53-fold more frequent than the average mutational rate along the transporter protein sequence. These data provide a genetic proof of the importance of these domain-domain interactions in the ABCC6 transporter.

Dichotomous haplotypic lineages of the immunoproteasome subunit genes, PSMB8 and PSMB10, in the MHC class I region of a Teleost Medaka, Oryzias latipes

Sequence comparison of the medaka, Oryzias latipes, major histocompatibility complex (MHC) class I region between two inbred strains, the HNI (derived from the Northern Population) and the Hd-rR (from the Southern Population), revealed a approximately 100 kb highly divergent segment encompassing two MHC class IA genes, Orla-UAA and Orla-UBA, and two immunoproteasome beta subunit genes, PSMB8 and PSMB10. To elucidate the genetic diversity of this region, we analyzed polymorphisms of the PSMB8 and PSMB10 genes using wild populations of medaka from three genetically different groups: the Northern Population, the Southern Population, and the China-West Korean Population. A total of 1,245 specimens from 10 localities were analyzed, and all the PSMB8 and PSMB10 alleles were classified into the N (fixed in the HNI strain) or the d (fixed in the Hd-rR strain) lineage. Polymerase chain reaction analysis of the region from PSMB8 to PSMB10 indicated that the two allelic lineages of these genes are segregating together constituting dichotomous haplotypic lineages. Both haplotypic lineages were identified in all three groups, although the frequency of d haplotypic lineage (73-100%) was much higher than that of N haplotypic lineage (0-27%) in all analyzed populations. The two allelic lineages of the PSMB8 gene showed curious substitutions at the 31st and 53rd residues of the mature peptide, which are likely involved in formation of the S1 pocket, suggesting that these alleles have a functional difference in cleavage specificity. These results indicate that the two medaka MHC haplotypic lineages encompassing the PSMB8 and PSMB10 genes are maintained in wild populations by a balancing selection.

Crystal structure of the conserved N-terminal domain of the peroxisomal matrix protein import receptor, Pex14p

Pex14p is a central component of the peroxisomal protein import machinery, in which the conserved N-terminal domain mediates dynamic interactions with other peroxins including Pex5p, Pex13p, and Pex19p. Here, we report the crystal structure of the conserved N-terminal domain of Pex14p with a three-helix bundle. A hydrophobic surface is composed of the conserved residues, of which two phenylalanine residues (Phe-35 and Phe-52) protrude to the solvent. Consequently, two putative binding pockets suitable for recognizing the helical WXXXF/Y motif of Pex5p are formed on the surface by the two phenylalanine residues accompanying with positively charged residues. The structural feature agrees well with our earlier findings where F35A/L36A and F52A/L53A mutants were impaired in the interactions with other peroxins such as Pex5p and Pex13p. Pex14p variants each with Phe-to-Ala mutation at positions 35, 52, and 35/52, respectively, were defective in restoring the impaired peroxisomal protein import in pex14 Chinese hamster ovary mutant ZP161 cells. Moreover, in GST pull-down assays His(6)-Pex5pL bound only to GST-Pex14p(25-70), not to any of GST-Pex14p(25-70)F35A, GST-Pex14p(25-70)F52A, and GST-Pex14p(25-70)F35A/F52A. Endogenous Pex5p was recruited to FLAG-Pex14p on peroxisomes in vivo but barely to FLAG-Pex14pF35A, FLAG-Pex14pF52A, and FLAG-Pex14pF35A/F52A. Collectively, Phe-35 and Phe-52 are essential for the Pex14p functions, including the interaction between Pex14p and Pex5p.

Properties of cytochrome bd plastoquinol oxidase from the cyanobacterium Synechocystis sp. PCC 6803

In the aerobic respiratory chain of the cyanobacterium Synechocystis sp. PCC 6803, cytochrome c oxidase serves as a major terminal oxidase while cyanide-resistant cytochrome bd serves as an alternative oxidase and evades the over-reduction of the plastoquinone pool under stress conditions. Here we expressed Synechocystis cytochrome bd in Escherichia coli and characterized enzymatic and spectroscopic properties. Cyanobacterial cytochrome bd showed the higher activity with ubiquinols than with decyl-plastoquinol and K(m) values for quinols were 2-fold smaller than those of E. coli cytochrome bd (CydAB). The dioxygen reduction site was resistant to cyanide as in E. coli oxidase while the quinol oxidation site was more sensitive to antimycin A and quinolone inhibitors. Spectroscopic analysis showed the presence of the haem b(595)-d binuclear centre but the sequence analysis indicates that cyanobacterial cytochrome bd is structurally related to cyanide-insensitive oxidase (CioAB), which does not show typical spectral changes upon reduction and ligand binding. Our data indicate that cyanobacterial cytochrome bd has unique enzymatic and structural properties and we hope that our findings will help our understanding the role and properties of CydAB and CioAB quinol oxidases in other bacterial species.

Evolutionary analysis of the Sub1 gene cluster that confers submergence tolerance to domesticated rice

BACKGROUND AND AIMS

Tolerance of complete submergence is recognized in a small number of accessions of domesticated Asian rice (Oryza sativa) and can be conferred by the Sub1A-1 gene of the polygenic Submergence-1 (Sub1) locus. In all O. sativa varieties, the Sub1 locus encodes the ethylene-responsive factor (ERF) genes Sub1B and Sub1C. A third paralogous ERF gene, Sub1A, is limited to a subset of indica accessions. It is thought that O. sativa was domesticated from the gene pools of the wild perennial species O. rufipogon Griff. and/or the annual species O. nivara Sharma et Shastry. The aim of this study was to evaluate the orthologues of the Sub1 locus in the closest relatives of O. sativa to provide insight into the origin of the gene and allelic variation of the Sub1 locus.

METHODS

Orthologues of the Sub1 genes were isolated from O. rufipogon and O. nivara by use of oligonucleotide primers corresponding to the most highly conserved regions of the Sub1 genes of domesticated rice. The phylogenetic relatedness of Sub1 genes of O. sativa and its wild relatives was evaluated.

KEY RESULTS AND CONCLUSIONS

Both O. rufipogon and O. nivara possess two Sub1 gene orthologues with strong sequence identity to the Sub1B and Sub1C alleles of cultivated rice. The phylogeny of the Sub1 genes of the domesticated and wild rice suggests that Sub1A arose from duplication of Sub1B. Variation in Sub1B alleles is correlated with the absence or presence of Sub1A. Together, the results indicate that genetic variation at the Sub1 locus is due to gene duplication and divergence that have occurred both prior to and after rice domestication.

HHP1 is involved in osmotic stress sensitivity in Arabidopsis

HHP1 (heptahelical protein 1), a protein with a predicted seven transmembrane domain structure homologous to adiponectin receptors (AdipoRs) and membrane progestin receptors (mPRs), has been characterized. Expression of HHP1 was increased in response to abscisic acid (ABA) and salt/osmotic stress as shown by quantitative real-time PCR and HHP1 promoter-controlled GUS activity. The HHP1 T-DNA insertion mutant (hhp1-1) showed a higher sensitivity to ABA and osmotic stress than the wild-type (WT), as revealed by the germination rate and post-germination growth rate. The induced expression of stress-responsive genes (RD29A, RD29B, ADH1, KIN1, COR15A, and COR47) was more sensitive to exogenous ABA and osmotic stress in hhp1-1 than in the WT. The hypersensitivity in the hhp1-1 mutant was reversed in the complementation mutant of HHP1 expressing the HHP1 gene. The data suggest that the mutation of HHP1 renders plants hypersensitive to ABA and osmotic stress and HHP1 might be a negative regulator in ABA and osmotic signalling.

Egyptian glycogen storage disease type III - identification of six novel AGL mutations, including a large 1.5 kb deletion and a missense mutation p.L620P with subtype IIId

BACKGROUND

Glycogen storage disease type III (GSD III) is caused by mutations in AGL which encodes for a single protein with two enzyme activities: oligo-1, 4-1, 4-glucantransferase (transferase) and amylo-1, 6-glucosidase. Activity of both enzymes is lost in most patients with GSD III, but in the very rare subtype IIId, transferase activity is deficient. Since the spectrum of AGL mutations is dependent on the ethnic group, we investigated the clinical and molecular characteristics in Egyptian patients with GSD III.

METHODS

Clinical features were examined in five Egyptian patients. AGL was sequenced and AGL haplotypes were determined.

RESULTS

Six novel AGL mutations were identified: a large deletion (c.3481-3588+1417del1525 bp), two insertions (c.1389insG and c.2368insA), two small deletions (c.2223-2224delGT and c.4041delT), and a missense mutation (p.L620P). p.L620P was found in a patient with IIId. Each mutation was located on a different AGL haplotype.

CONCLUSIONS

Our results suggest that there is allelic and phenotypic heterogeneity of GSD III in Egypt. This is the second description of a large deletion in AGL. p.L620P is the second mutation found in GSD IIId.

Function, regulation, and transcriptional organization of the hemin utilization locus of Bartonella quintana

Bartonella quintana is a gram-negative agent of trench fever, chronic bacteremia, endocarditis, and bacillary angiomatosis in humans. B. quintana has the highest known hemin requirement among bacteria, but the mechanisms of hemin acquisition are poorly defined. Genomic analyses revealed a potential locus dedicated to hemin utilization (hut) encoding a putative hemin receptor, HutA; a TonB-like energy transducer; an ABC transport system comprised of three proteins, HutB, HutC, and HmuV; and a hemin degradation/storage enzyme, HemS. Complementation analyses with Escherichia coli hemA show that HutA functions as a hemin receptor, and complementation analyses with E. coli hemA tonB indicate that HutA is TonB dependent. Quantitative reverse transcriptase PCR analyses show that hut locus transcription is subject to hemin-responsive regulation, which is mediated primarily by the iron response regulator (Irr). Irr functions as a transcriptional repressor of the hut locus at all hemin concentrations tested. Overexpression of the ferric uptake regulator (fur) represses transcription of tonB in the presence of excess hemin, whereas overexpression of the rhizobial iron regulator (rirA) has no effect on hut locus transcription. Reverse transcriptase PCR analyses show that hutA and tonB are divergently transcribed and that the remaining hut genes are expressed as a polycistronic mRNA. Examination of the promoter regions of hutA, tonB, and hemS reveals consensus sequence promoters that encompass an H-box element previously shown to interact with B. quintana Irr.

The basic concepts of molecular modeling

Molecular modeling techniques have made significant advances in recent years and are becoming essential components of many chemical, physical and biological studies. Here we present three widely used techniques used in the simulation of biomolecular systems: structural and homology modeling, molecular dynamics and molecular docking. For each of these topics we present a brief discussion of the underlying scientific basis of the technique, some simple examples of how the method is commonly applied, and some discussion of the limitations and caveats of which the user should be aware. References for further reading as well as an extensive list of software resources are provided.

Kinomer v. 1.0: a database of systematically classified eukaryotic protein kinases

The regulation of protein function through reversible phosphorylation by protein kinases and phosphatases is a general mechanism controlling virtually every cellular activity. Eukaryotic protein kinases can be classified into distinct, well-characterized groups based on amino acid sequence similarity and function. We recently reported a highly sensitive and accurate hidden Markov model-based method for the automatic detection and classification of protein kinases into these specific groups. The Kinomer v. 1.0 database presented here contains annotated classifications for the protein kinase complements of 43 eukaryotic genomes. These span the taxonomic range and include fungi (16 species), plants (6), diatoms (1), amoebas (2), protists (1) and animals (17). The kinomes are stored in a relational database and are accessible through a web interface on the basis of species, kinase group or a combination of both. In addition, the Kinomer v. 1.0 HMM library is made available for users to perform classification on arbitrary sequences. The Kinomer v. 1.0 database is a continually updated resource where direct comparison of kinase sequences across kinase groups and across species can give insights into kinase function and evolution. Kinomer v. 1.0 is available at http://www.compbio.dundee.ac.uk/kinomer/.

Computational analysis of the human HSPH/HSPA/DNAJ family and cloning of a human HSPH/HSPA/DNAJ expression library

In this manuscript, we describe the generation of a gene library for the expression of HSP110/HSPH, HSP70/HSPA and HSP40/DNAJ members. First, the heat shock protein (HSP) genes were collected from the gene databases and the gene families were analyzed for expression patterns, heat inducibility, subcellular localization, and protein homology using several bioinformatics approaches. These results can be used as a working draft model until data are confirmed by experimental approaches. In addition, we describe the generation of a HSPA/DNAJ overexpression library and tested the effect of different fusion tags on HSPA and DNAJ members using different techniques for measuring chaperone activity. These results show that we have cloned a high-quality heat shock protein expression library containing most members from the HSPH, HSPA, DNAJA and DNAJB families which will be useful for the chaperone community to unravel the function of the highly diverse family of human molecular chaperones.

Ethylene regulates phosphorus remobilization and expression of a phosphate transporter (PhPT1) during petunia corolla senescence

The programmed degradation of macromolecules during petal senescence allows the plant to remobilize nutrients from dying to developing tissues. Ethylene is involved in regulating the timing of nucleic acid degradation in petunia, but it is not clear if ethylene has a role in the remobilization of phosphorus during petal senescence. To investigate ethylene's role in nutrient remobilization, the P content of petals (collectively called the corolla) during early development and senescence was compared in ethylene-sensitive wild type Petunia x hybrida 'Mitchell Diploid' (MD) and transgenic petunias with reduced sensitivity to ethylene (35S::etr1-1). When compared to the total P content of corollas on the day of flower opening (the early non-senescing stage), P in MD corollas had decreased 74% by the late stage of senescence (advanced wilting). By contrast, P levels were only reduced by an average of 32% during etr1-1 corolla (lines 44568 and Z00-35-10) senescence. A high-affinity phosphate transporter, PhPT1 (PhPht1;1), was cloned from senescing petunia corollas by RT-PCR. PhPT1 expression was up-regulated during MD corolla senescence and a much smaller increase was detected during the senescence of etr1-1 petunia corollas. PhPT1 mRNA levels showed a rapid increase in detached corollas (treated at 1 d after flower opening) following treatment with low levels of ethylene (0.1 microl l(-1)). Transcripts accumulated in the presence of the protein synthesis inhibitor, cycloheximide, indicating that PhPT1 is a primary ethylene response gene. PhPT1 is a putative phosphate transporter that may function in Pi translocation during senescence.

Proteomic tools for the analysis of cytoskeleton proteins

Proteomic tools have become an essential part of the tool kit of the molecular biologist, and provide techniques for detecting homologous sequences, recognizing functional domains, modeling, and analyzing the three-dimensional structure for any given protein sequence. Although a wealth of structural and functional information is available for a large number of members of the various classes of cytoskeletal proteins, many more members remain uncharacterized. These computational tools that are freely and easily accessible to the scientific community provide an excellent starting point to predict the structural and functional properties of such partially or fully uncharacterized protein sequences, and can lead to elegantly designed experiments to probe the hypothesized function. This chapter discusses various proteomic analysis tools with a focus on protein structure and function predictions.

Are coral-dwelling crabs monophyletic? A phylogeny of the Trapezioidea (Crustacea : Decapoda : Brachyura)

Many brachyuran crab genera are known to be obligate symbionts on scleractinian corals in shallow tropical seas (e.g. Cymo, Tanaocheles, Domecia, Tetralia, Tetraloides and Trapezia). Most are currently placed in several families and superfamilies, the most prominent being the Trapezioidea, containing Domeciidae, Tetraliidae and Trapeziidae. The distinctness of the domeciids has long been recognised in Domecinae, and finally Domeciidae. Tetraliidae, however, has only recently been recognised when the classical coral obligate family Trapeziidae was split into two families, Trapeziidae sensu stricto and Tetraliidae. Recognition of Tetraliidae, however, has been controversial, and some workers have suggested that it should be regarded as a synonym of Trapeziidae, principally on the basis of carapace morphology and larval features. Using two gene markers, we evaluate the status of Trapezioidea and, in particular, that of Tetraliidae. Our analyses support the distinction between Tetraliidae and Trapeziidae and suggest that crabs have independently colonised hard corals several times.

Origins of the Xylella fastidiosa prophage-like regions and their impact in genome differentiation

Xylella fastidiosa is a Gram negative plant pathogen causing many economically important diseases, and analyses of completely sequenced X. fastidiosa genome strains allowed the identification of many prophage-like elements and possibly phage remnants, accounting for up to 15% of the genome composition. To better evaluate the recent evolution of the X. fastidiosa chromosome backbone among distinct pathovars, the number and location of prophage-like regions on two finished genomes (9a5c and Temecula1), and in two candidate molecules (Ann1 and Dixon) were assessed. Based on comparative best bidirectional hit analyses, the majority (51%) of the predicted genes in the X. fastidiosa prophage-like regions are related to structural phage genes belonging to the Siphoviridae family. Electron micrograph reveals the existence of putative viral particles with similar morphology to lambda phages in the bacterial cell in planta. Moreover, analysis of microarray data indicates that 9a5c strain cultivated under stress conditions presents enhanced expression of phage anti-repressor genes, suggesting switches from lysogenic to lytic cycle of phages under stress-induced situations. Furthermore, virulence-associated proteins and toxins are found within these prophage-like elements, thus suggesting an important role in host adaptation. Finally, clustering analyses of phage integrase genes based on multiple alignment patterns reveal they group in five lineages, all possessing a tyrosine recombinase catalytic domain, and phylogenetically close to other integrases found in phages that are genetic mosaics and able to perform generalized and specialized transduction. Integration sites and tRNA association is also evidenced. In summary, we present comparative and experimental evidence supporting the association and contribution of phage activity on the differentiation of Xylella genomes.

Analysis of ColE1 MbeC unveils an extended ribbon-helix-helix family of nicking accessory proteins

MbeC is a 13-kDa ColE1-encoded protein required for efficient mobilization of ColE1, a plasmid widely used in cloning vector technology. MbeC protein was purified and used for in vitro DNA binding, which showed that it binds specifically double-stranded DNA (dsDNA) containing the ColE1 oriT. Amino acid sequence comparison and secondary structure prediction imply that MbeC is related to the ribbon-helix-helix (RHH) protein family. Alignment with RHH members pointed to a conserved arginine (R13 in MbeC) that was mutated to alanine. The mutant MbeC(R13A) was unable to bind either single-stranded DNA or dsDNA. Limited proteolysis fragmented MbeC in two stable folding domains: the N-terminal domain, which contains the RHH motif, and the C-terminal domain, which comprises a signature shared by nicking accessory proteins. The results indicate that MbeC plays a similar role in conjugation as TraY and TrwA of plasmids F and R388, respectively. Thus, it appears that an extended, possibly universal mechanism of DNA conjugative processing exists, in which oriT-processing is carried out by relaxases assisted by homologous nicking accessory proteins. This mechanism seems to be shared by all major conjugative systems analyzed thus far.

Role of conjugative elements in the evolution of the multidrug-resistant pandemic clone Streptococcus pneumoniaeSpain23F ST81

Streptococcus pneumoniae is a human commensal and pathogen able to cause a variety of diseases that annually result in over a million deaths worldwide. The S. pneumoniae(Spain23F) sequence type 81 lineage was among the first recognized pandemic clones and was responsible for almost 40% of penicillin-resistant pneumococcal infections in the United States in the late 1990s. Analysis of the chromosome sequence of a representative strain, and comparison with other available genomes, indicates roles for integrative and conjugative elements in the evolution of pneumococci and, more particularly, the emergence of the multidrug-resistant Spain 23F ST81 lineage. A number of recently acquired loci within the chromosome appear to encode proteins involved in the production of, or immunity to, antimicrobial compounds, which may contribute to the proficiency of this strain at nasopharyngeal colonization. However, further sequencing of other pandemic clones will be required to establish whether there are any general attributes shared by these strains that are responsible for their international success.

Kalign2: high-performance multiple alignment of protein and nucleotide sequences allowing external features

In the growing field of genomics, multiple alignment programs are confronted with ever increasing amounts of data. To address this growing issue we have dramatically improved the running time and memory requirement of Kalign, while maintaining its high alignment accuracy. Kalign version 2 also supports nucleotide alignment, and a newly introduced extension allows for external sequence annotation to be included into the alignment procedure. We demonstrate that Kalign2 is exceptionally fast and memory-efficient, permitting accurate alignment of very large numbers of sequences. The accuracy of Kalign2 compares well to the best methods in the case of protein alignments while its accuracy on nucleotide alignments is generally superior. In addition, we demonstrate the potential of using known or predicted sequence annotation to improve the alignment accuracy. Kalign2 is freely available for download from the Kalign web site (http://msa.sbc.su.se/).