Search

Type III Secretion Systems (T3SSs) are structurally conserved nanomachines that span the inner and outer bacterial membranes, and via a protruding needle complex contact host cell membranes and deliver type III effector proteins. T3SS are phylogenetically divided into several families based on structural basal body components. Here we have studied the evolutionary and functional conservation of four T3SS proteins from the Inv/Mxi-Spa family: a cytosolic chaperone, two hydrophobic translocators that form a plasma membrane-integral pore, and the hydrophilic 'tip complex' translocator that connects the T3SS needle to the translocon pore. Salmonella enterica serovar Typhimurium (S. Typhimurium), a common cause of food-borne gastroenteritis, possesses two T3SSs, one belonging to the Inv/Mxi-Spa family. We used invasion-deficient S. Typhimurium mutants as surrogates for expression of translocator orthologs identified from an extensive phylogenetic analysis, and type III effector translocation and host cell invasion as a readout for complementation efficiency, and identified several Inv/Mxi-Spa orthologs that can functionally substitute for the S. Typhimurium chaperone and translocator proteins. Functional complementation correlates with amino acid sequence identity between orthologs, but varies considerably between the four proteins. This is the first in-depth survey of the functional interchangeability of Inv/Mxi-Spa T3SS proteins acting directly at the host-pathogen interface., (© 2016 John Wiley & Sons Ltd.)

The Comprehensive Antibiotic Resistance Database (CARD; http://arpcard.mcmaster.ca) is a manually curated resource containing high quality reference data on the molecular basis of antimicrobial resistance (AMR), with an emphasis on the genes, proteins and mutations involved in AMR. CARD is ontologically structured, model centric, and spans the breadth of AMR drug classes and resistance mechanisms, including intrinsic, mutation-driven and acquired resistance. It is built upon the Antibiotic Resistance Ontology (ARO), a custom built, interconnected and hierarchical controlled vocabulary allowing advanced data sharing and organization. Its design allows the development of novel genome analysis tools, such as the Resistance Gene Identifier (RGI) for resistome prediction from raw genome sequence. Recent improvements include extensive curation of additional reference sequences and mutations, development of a unique Model Ontology and accompanying AMR detection models to power sequence analysis, new visualization tools, and expansion of the RGI for detection of emergent AMR threats. CARD curation is updated monthly based on an interplay of manual literature curation, computational text mining, and genome analysis., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The loss of effective antimicrobials is reducing our ability to protect the global population from infectious disease. However, the field of antibiotic drug discovery and the public health monitoring of antimicrobial resistance (AMR) is beginning to exploit the power of genome and metagenome sequencing. The creation of novel AMR bioinformatics tools and databases and their continued development will advance our understanding of the molecular mechanisms and threat severity of antibiotic resistance, while simultaneously improving our ability to accurately predict and screen for antibiotic resistance genes within environmental, agricultural, and clinical settings. To do so, efforts must be focused toward exploiting the advancements of genome sequencing and information technology. Currently, AMR bioinformatics software and databases reflect different scopes and functions, each with its own strengths and weaknesses. A review of the available tools reveals common approaches and reference data but also reveals gaps in our curated data, models, algorithms, and data-sharing tools that must be addressed to conquer the limitations and areas of unmet need within the AMR research field before DNA sequencing can be fully exploited for AMR surveillance and improved clinical outcomes., (© 2016 New York Academy of Sciences.)

Antibiotic resistance is ancient and widespread in environmental bacteria. These are therefore reservoirs of resistance elements and reflective of the natural history of antibiotics and resistance. In a previous study, we discovered that multi-drug resistance is common in bacteria isolated from Lechuguilla Cave, an underground ecosystem that has been isolated from the surface for over 4 Myr. Here we use whole-genome sequencing, functional genomics and biochemical assays to reveal the intrinsic resistome of Paenibacillus sp. LC231, a cave bacterial isolate that is resistant to most clinically used antibiotics. We systematically link resistance phenotype to genotype and in doing so, identify 18 chromosomal resistance elements, including five determinants without characterized homologues and three mechanisms not previously shown to be involved in antibiotic resistance. A resistome comparison across related surface Paenibacillus affirms the conservation of resistance over millions of years and establishes the longevity of these genes in this genus.

Development is a complex and well-defined process characterized by rapid cell proliferation and apoptosis. At this stage in life, a developmentally young organism is more sensitive to toxicants as compared to an adult. In response to pro-oxidant exposure, members of the Cap'n'Collar (CNC) basic leucine zipper (b-ZIP) transcription factor family (including Nfe2 and Nfe2-related factors, Nrfs) activate the expression of genes whose protein products contribute to reduced toxicity. Here, we studied the role of the CNC protein, Nfe2, in the developmental response to pro-oxidant exposure in the zebrafish (Danio rerio). Following acute waterborne exposures to diquat or tert-buytlhydroperoxide (tBOOH) at one of three developmental stages, wildtype (WT) and nfe2 knockout (KO) embryos and larvae were morphologically scored and their transcriptomes sequenced. Early in development, KO animals suffered from hypochromia that was made more severe through exposure to pro-oxidants; this phenotype in the KO may be linked to decreased expression of alas2, a gene involved in heme synthesis. WT and KO eleutheroembryos and larvae were phenotypically equally affected by exposure to pro-oxidants, where tBOOH caused more pronounced phenotypes as compared to diquat. Comparing diquat and tBOOH exposed embryos relative to the WT untreated control, a greater number of genes were up-regulated in the tBOOH condition as compared to diquat (tBOOH: 304 vs diquat: 148), including those commonly found to be differentially regulated in the vertebrate oxidative stress response (OSR) (e.g. hsp70.2, txn1, and gsr). When comparing WT and KO across all treatments and times, there were 1170 genes that were differentially expressed, of which 33 are known targets of the Nrf proteins Nrf1 and Nrf2. More specifically, in animals exposed to pro-oxidants a total of 968 genes were differentially expressed between WT and KO across developmental time, representing pathways involved in coagulation, embryonic organ development, body fluid level regulation, erythrocyte differentiation, and oxidation-reduction, amongst others. The greatest number of genes that changed in expression between WT and KO occurred in animals exposed to diquat at 2h post fertilization (hpf). Across time and treatment, there were six genes (dhx40, cfap70, dnajb9b, slc35f4, spi-c, and gpr19) that were significantly up-regulated in KO compared to WT and four genes (fhad1, cyp4v7, nlrp12, and slc16a6a) that were significantly down-regulated. None of these genes have been previously identified as targets of Nfe2 or the Nrf family. These results demonstrate that the zebrafish Nfe2 may be a regulator of both primitive erythropoiesis and the OSR during development., (Copyright © 2016 Elsevier B.V. All rights reserved.)

YphC and YsxC are GTPases in Bacillus subtilis that facilitate the assembly of the 50S ribosomal subunit, however their roles in this process are still uncharacterized. To explore their function, we used strains in which the only copy of the yphC or ysxC genes were under the control of an inducible promoter. Under depletion conditions, they accumulated incomplete ribosomal subunits that we named 45SYphC and 44.5SYsxC particles. Quantitative mass spectrometry analysis and the 5-6 Å resolution cryo-EM maps of the 45SYphC and 44.5SYsxC particles revealed that the two GTPases participate in the maturation of the central protuberance, GTPase associated region and key RNA helices in the A, P and E functional sites of the 50S subunit. We observed that YphC and YsxC bind specifically to the two immature particles, suggesting that they represent either on-pathway intermediates or that their structure has not significantly diverged from that of the actual substrate. These results describe the nature of these immature particles, a widely used tool to study the assembly process of the ribosome. They also provide the first insights into the function of YphC and YsxC in 50S subunit assembly and are consistent with this process occurring through multiple parallel pathways, as it has been described for the 30S subunit., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Genes of the major histocompatibility complex (MHC) exhibit heterozygote advantage in immune defence, which in turn can select for MHC-disassortative mate choice. However, many species lack this expected pattern of MHC-disassortative mating. A possible explanation lies in evolutionary processes following gene duplication: if two duplicated MHC genes become functionally diverged from each other, offspring will inherit diverse multilocus genotypes even under random mating. We used locus-specific primers for high-throughput sequencing of two expressed MHC Class II B genes in Leach's storm-petrels, Oceanodroma leucorhoa, and found that exon 2 alleles fall into two gene-specific monophyletic clades. We tested for disassortative vs. random mating at these two functionally diverged Class II B genes, using multiple metrics and different subsets of exon 2 sequence data. With good statistical power, we consistently found random assortment of mates at MHC. Despite random mating, birds had MHC genotypes with functionally diverged alleles, averaging 13 amino acid differences in pairwise comparisons of exon 2 alleles within individuals. To test whether this high MHC diversity in individuals is driven by evolutionary divergence of the two duplicated genes, we built a phylogenetic permutation model. The model showed that genotypic diversity was strongly impacted by sequence divergence between the most common allele of each gene, with a smaller additional impact of monophyly of the two genes. Divergence of allele sequences between genes may have reduced the benefits of actively seeking MHC-dissimilar mates, in which case the evolutionary history of duplicated genes is shaping the adaptive landscape of sexual selection., (© 2016 John Wiley & Sons Ltd.)

Degraded DNA from suboptimal field sampling is common in molecular ecology. However, its impact on techniques that use restriction site associated next-generation DNA sequencing (RADSeq, GBS) is unknown. We experimentally examined the effects of in situDNA degradation on data generation for a modified double-digest RADSeq approach (3RAD). We generated libraries using genomic DNA serially extracted from the muscle tissue of 8 individual lake whitefish (Coregonus clupeaformis) following 0-, 12-, 48- and 96-h incubation at room temperature posteuthanasia. This treatment of the tissue resulted in input DNA that ranged in quality from nearly intact to highly sheared. All samples were sequenced as a multiplexed pool on an Illumina MiSeq. Libraries created from low to moderately degraded DNA (12-48 h) performed well. In contrast, the number of RADtags per individual, number of variable sites, and percentage of identical RADtags retained were all dramatically reduced when libraries were made using highly degraded DNA (96-h group). This reduction in performance was largely due to a significant and unexpected loss of raw reads as a result of poor quality scores. Our findings remained consistent after changes in restriction enzymes, modified fold coverage values (2- to 16-fold), and additional read-length trimming. We conclude that starting DNA quality is an important consideration for RADSeq; however, the approach remains robust until genomic DNA is extensively degraded., (© 2015 John Wiley & Sons Ltd.)

Antimicrobial resistance is a global health challenge and has an evolutionary trajectory ranging from proto-resistance in the environment to untreatable clinical pathogens. Resistance is not static, as pathogenic strains can move among patient populations and individual resistance genes can move among pathogens. Effective treatment of resistant infections, antimicrobial stewardship, and new drug discovery increasingly rely upon genotype information, powered by decreasing costs of DNA sequencing. These new approaches will require advances in microbial informatics, particularly in development of reference databases of molecular determinants such as our Comprehensive Antibiotic Resistance Database and clinical metadata, new algorithms for prediction of resistome and resistance phenotype from genotype, and new protocols for global collection and sharing of high-throughput molecular epidemiology data., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

The International Pseudomonas aeruginosa Consortium is sequencing over 1000 genomes and building an analysis pipeline for the study of Pseudomonas genome evolution, antibiotic resistance and virulence genes. Metadata, including genomic and phenotypic data for each isolate of the collection, are available through the International Pseudomonas Consortium Database (http://ipcd.ibis.ulaval.ca/). Here, we present our strategy and the results that emerged from the analysis of the first 389 genomes. With as yet unmatched resolution, our results confirm that P. aeruginosa strains can be divided into three major groups that are further divided into subgroups, some not previously reported in the literature. We also provide the first snapshot of P. aeruginosa strain diversity with respect to antibiotic resistance. Our approach will allow us to draw potential links between environmental strains and those implicated in human and animal infections, understand how patients become infected and how the infection evolves over time as well as identify prognostic markers for better evidence-based decisions on patient care.

IslandViewer (http://pathogenomics.sfu.ca/islandviewer) is a widely used web-based resource for the prediction and analysis of genomic islands (GIs) in bacterial and archaeal genomes. GIs are clusters of genes of probable horizontal origin, and are of high interest since they disproportionately encode genes involved in medically and environmentally important adaptations, including antimicrobial resistance and virulence. We now report a major new release of IslandViewer, since the last release in 2013. IslandViewer 3 incorporates a completely new genome visualization tool, IslandPlot, enabling for the first time interactive genome analysis and gene search capabilities using synchronized circular, horizontal and vertical genome views. In addition, more curated virulence factors and antimicrobial resistance genes have been incorporated, and homologs of these genes identified in closely related genomes using strict filters. Pathogen-associated genes have been re-calculated for all pre-computed complete genomes. For user-uploaded genomes to be analysed, IslandViewer 3 can also now handle incomplete genomes, with an improved queuing system on compute nodes to handle user demand. Overall, IslandViewer 3 represents a significant new version of this GI analysis software, with features that may make it more broadly useful for general microbial genome analysis and visualization., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The major histocompatibility complex (Mhc) is subject to pathogen-mediated balancing selection and can link natural selection with mate choice. We characterized two Mhc class II B loci in Leach's storm-petrel, Oceanodroma leucorhoa, focusing on exon 2 which encodes the portion of the protein that binds pathogen peptides. We amplified and sequenced exon 2 with locus-specific nested PCR and Illumina MiSeq using individually barcoded primers. Repeat genotyping of 78 single-locus genotypes produced identical results in 77 cases (98.7%). Sequencing of messenger RNA (mRNA) from three birds confirmed expression of both loci, consistent with the observed absence of stop codons or frameshifts in all alleles. In 48 birds, we found 9 and 12 alleles at the two loci, respectively, and all 21 alleles translated to unique amino acid sequences. Unlike many studies of duplicated Mhc genes, alleles of the two loci clustered into monophyletic groups. Consistent with this phylogenetic result, interlocus gene conversion appears to have affected only two short fragments of the exon. As predicted under a paradigm of pathogen-mediated selection, comparison of synonymous and non-synonymous substitution rates found evidence of a history of positive selection at putative peptide binding sites. Overall, the results suggest that the gene duplication event leading to these two loci is not recent and that point mutations and positive selection on the peptide binding sites may be the predominant forces acting on these genes. Characterization of these loci sets the stage for population-level work on the evolutionary ecology of Mhc in this species.

Oxidative stress is an important mechanism of chemical toxicity, contributing to teratogenesis and to cardiovascular and neurodegenerative diseases. Developing animals may be especially sensitive to chemicals causing oxidative stress. The developmental expression and inducibility of anti-oxidant defenses through activation of NF-E2-related factor 2 (NRF2) affect susceptibility to oxidants, but the embryonic response to oxidants is not well understood. To assess the response to chemically mediated oxidative stress and how it may vary during development, zebrafish embryos, eleutheroembryos, or larvae at 1, 2, 3, 4, 5, and 6 days post fertilization (dpf) were exposed to DMSO (0.1%), tert-butylhydroquinone (tBHQ; 10 µM) or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD; 2 nM) for 6 hr. Transcript abundance was assessed by real-time qRT-PCR and microarray. qRT-PCR showed strong (4- to 5-fold) induction of gstp1 by tBHQ as early as 1 dpf. tBHQ also induced gclc (2 dpf), but not sod1, nqo1, or cyp1a. TCDD induced cyp1a but none of the other genes. Microarray analysis showed that 1477 probes were significantly different among the DMSO-, tBHQ-, and TCDD-treated eleutheroembryos at 4 dpf. There was substantial overlap between genes induced in developing zebrafish and a set of marker genes induced by oxidative stress in mammals. Genes induced by tBHQ in 4-dpf zebrafish included those involved in glutathione synthesis and utilization, signal transduction, and DNA damage/stress response. The strong induction of hsp70 determined by microarray was confirmed by qRT-PCR and by use of transgenic zebrafish expressing enhanced green fluorescent protein (EGFP) under control of the hsp70 promoter. Genes strongly down-regulated by tBHQ included mitfa, providing a molecular explanation for the loss of pigmentation in tBHQ-exposed embryos. These data show that zebrafish embryos are responsive to oxidative stress as early as 1 dpf, that responsiveness varies with development in a gene-specific manner, and that the oxidative stress response is substantially conserved in vertebrate animals.

Serial Analysis of Gene Expression (SAGE) was used to quantify transcriptional changes in Giardia intestinalis during its interaction with human intestinal epithelial cells (IECs, HT-29) in serum free M199 medium. Transcriptional changes were compared to those in trophozoites alone in M199 and in TYI-S-33 Giardia growth medium. In total, 90 genes were differentially expressed, mainly those involved in cellular redox homeostasis, metabolism and small molecule transport but also cysteine proteases and structural proteins of the giardin family. Only 29 genes changed their expression due to IEC interaction and the rest were due to M199 medium. Although our findings generated a small dataset, it was consistent with our earlier microarray studies performed under different interaction conditions. This study has confined the number of genes in Giardia to a small subset that specifically change their expression due to interaction with IECs., (Copyright © 2014 Elsevier B.V. All rights reserved.)

The metal responsive element-binding transcription factor-1 (MTF-1) responds to changes in cellular zinc levels caused by zinc exposure or disruption of endogenous zinc homeostasis by heavy metals or oxygen-related stress. Here we report the functional characterization of a complete zebrafish MTF-1 in comparison with the previously identified isoform lacking the highly conserved cysteine-rich motif (Cys-X-Cys-Cys-X-Cys) found in all other vertebrate MTF-1 orthologs. In an effort to develop novel molecular tools, a constitutively nuclear dominant-negative MTF-1 (dnMTF-1) was generated as tool for inhibiting endogenous MTF-1 signaling. The in vivo efficacy of the dnMTF-1 was determined by microinjecting in vitro transcribed dnMTF-1 mRNA into zebrafish embryos (1-2 cell stage) followed by transcriptomic profiling using an Agilent 4x44K array on 28- and 36-hpf embryos. A total of 594 and 560 probes were identified as differentially expressed at 28hpf and 36hpf, respectively, with interesting overlaps between timepoints. The main categories of genes affected by the inhibition of MTF-1 signaling were: nuclear receptors and genes involved in stress signaling, neurogenesis, muscle development and contraction, eye development, and metal homeostasis, including novel observations in iron and heme homeostasis. Finally, we investigate both the transcriptional activator and transcriptional repressor role of MTF-1 in potential novel target genes identified by transcriptomic profiling during early zebrafish development., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Transcription factors in the CNC-bZIP family (NFE2, NRF1, NRF2 and NRF3) regulate genes with a wide range of functions in response to both physiological and exogenous signals, including those indicating changes in cellular redox status. Given their role in helping to maintain cellular homeostasis, it is imperative to understand the expression, regulation, and function of CNC-bZIP genes during embryonic development. We explored the expression and function of six nrf genes (nfe2, nrf1a, nrf1b, nrf2a, nrf2b, and nrf3) using zebrafish embryos as a model system. Analysis by microarray and quantitative RT-PCR showed that genes in the nrf family were expressed throughout development from oocytes to larvae. The spatial expression of nrf3 suggested a role in regulating the development of the brain, brachia and pectoral fins. Knock-down by morpholino anti-sense oligonucleotides suggested that none of the genes were necessary for embryonic viability, but nfe2 was required for proper cellular organization in the pneumatic duct and subsequent swim bladder function, as well as for proper formation of the otic vesicles. nrf genes were induced by the oxidant tert-butylhydroperoxide, and some of this response was regulated through family members Nrf2a and Nrf2b. Our results provide a foundation for understanding the role of nrf genes in normal development and in regulating the response to oxidative stress in vertebrate embryos.

The field of antibiotic drug discovery and the monitoring of new antibiotic resistance elements have yet to fully exploit the power of the genome revolution. Despite the fact that the first genomes sequenced of free living organisms were those of bacteria, there have been few specialized bioinformatic tools developed to mine the growing amount of genomic data associated with pathogens. In particular, there are few tools to study the genetics and genomics of antibiotic resistance and how it impacts bacterial populations, ecology, and the clinic. We have initiated development of such tools in the form of the Comprehensive Antibiotic Research Database (CARD; http://arpcard.mcmaster.ca). The CARD integrates disparate molecular and sequence data, provides a unique organizing principle in the form of the Antibiotic Resistance Ontology (ARO), and can quickly identify putative antibiotic resistance genes in new unannotated genome sequences. This unique platform provides an informatic tool that bridges antibiotic resistance concerns in health care, agriculture, and the environment.

NF-E2-related factor 2 (NRF2; also called NFE2L2) and related NRF family members regulate antioxidant defenses by activating gene expression via antioxidant response elements (AREs), but their roles in embryonic development are not well understood. We report here that zebrafish (Danio rerio), an important developmental model species, possesses six nrf genes, including duplicated nrf1 and nrf2 genes. We cloned a novel zebrafish nrf2 paralog, nrf2b. The predicted Nrf2b protein sequence shares several domains with the original Nrf2 (now Nrf2a) but lacks the Neh4 transactivation domain. Zebrafish-human comparisons demonstrate conserved synteny involving nrf2 and hox genes, indicating that nrf2a and nrf2b are co-orthologs of human NRF2. nrf2a and nrf2b displayed distinct patterns of expression during embryonic development; nrf2b was more highly expressed at all stages. Embryos in which Nrf2a expression had been knocked down with morpholino oligonucleotides were more sensitive to tert-butylhydroperoxide but not tert-butylhydroquinone, whereas knockdown of Nrf2b did not affect sensitivity of embryos to either chemical. Gene expression profiling by microarray identified a specific role for Nrf2b as a negative regulator of several genes, including p53, cyclin G1, and heme oxygenase 1, in embryos. Nrf2a and Nrf2b exhibited different mechanisms of cross-talk with the Ahr2 signaling pathway. Together, these results demonstrate distinct roles for nrf2a and nrf2b, consistent with subfunction partitioning, and identify a novel negative regulatory role for Nrf2b during development. The identification of zebrafish nrf2 co-orthologs will facilitate new understanding of the multiple roles of NRF2 in protecting vertebrate embryos from oxidative damage.

Kinase-mediated resistance to antibiotics is a significant clinical challenge. These enzymes share a common protein fold characteristic of Ser/Thr/Tyr protein kinases. We screened 14 antibiotic resistance kinases against 80 chemically diverse protein kinase inhibitors to map resistance kinase chemical space. The screens identified molecules with both broad and narrow inhibition profiles, proving that protein kinase inhibitors offer privileged chemical matter with the potential to block antibiotic resistance. One example is the flavonol quercetin, which inhibited a number of resistance kinases in vitro and in vivo. This activity was rationalized by determination of the crystal structure of the aminoglycoside kinase APH(2″)-IVa in complex with quercetin and its antibiotic substrate kanamycin. Our data demonstrate that protein kinase inhibitors offer chemical scaffolds that can block antibiotic resistance, providing leads for co-drug design., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Proteolytic activity is important in the lifecycles of parasites and their interactions with hosts. Cysteine proteases have been best studied in Giardia, but other protease classes have been implicated in growth and/or differentiation. In this study, we employed bioinformatics to reveal the complete set of putative proteases in the Giardia genome. We identified 73 peptidase homologs distributed over 5 catalytic classes in the genome. Serial analysis of gene expression of the G. lamblia lifecycle found thirteen protease genes with significant transcriptional variation over the lifecycle, with only one serine protease transcript upregulated late in encystation. The translated gene sequence of this encystation-specific transcript was most similar to eukaryotic subtilisin-like proprotein convertases (SPC), although the typical catalytic triad was not identified. Epitope-tagged gSPC protein expressed in Giardia under its own promoter was upregulated during encystation with highest expression in cysts and it localized to encystation-specific secretory vesicles (ESV). Total gSPC from encysting cells produced proteolysis in gelatin gels that co-migrated with the epitope-tagged protease in immunoblots. Immuno-purified gSPC also had gelatinase activity. To test whether endogenous gSPC activity is involved in differentiation, trophozoites and cysts were exposed to the specific serine proteinase inhibitor 4-(2-aminoethyl)-benzenesulfonyl fluoride hydrochloride (AEBSF). After 21 h encystation, a significant decrease in ESV was observed with 1mM AEBSF and by 42 h the number of cysts was significantly reduced, but trophozoite growth was not inhibited. Concurrently, levels of cyst wall proteins 1 and 2, and AU1-tagged gSPC protein itself were decreased. Excystation of G. muris cysts was also significantly reduced in the presence of AEBSF. These results support the idea that serine protease activity is essential for Giardia encystation and excystation., (Copyright © 2010 Elsevier B.V. All rights reserved.)

Cytochrome P450 (CYP) proteins compose a highly diverse superfamily found in all domains of life. These proteins are enzymes involved in metabolism of endogenous and exogenous compounds. In vertebrates, the CYP2 family is one of the largest, most diverse and plays an important role in mammalian drug metabolism. However, there are more than 20 vertebrate CYP2 subfamilies with uncertain evolution and fairly discrete subfamily composition within vertebrate classes, hindering extrapolation of knowledge across subfamilies. To better understand CYP2 diversity, a phylogenetic analysis of 196 CYP2 protein sequences from 16 species was performed using a maximum likelihood approach and Bayesian inference. The analyses included the CYP2 compliment from human, fugu, zebrafish, stickleback, medaka, cow, and dog genomes. Additional sequences were included from rabbit, marsupial, platypus, chicken, frog, and salmonid species. Three CYP2 sequences from the tunicate Ciona intestinalis were utilized as the outgroup. Results indicate a single ancestral vertebrate CYP2 gene and monophyly of all CYP2 subfamilies. Two subfamilies (CYP2R and CYP2U) pre-date vertebrate diversification, allowing direct comparison across vertebrate classes, while all other subfamilies originated during vertebrate diversification, often within specific vertebrate lineages. Analysis of site-specific evolution indicates that some substrate recognition sites (SRS) previously proposed for CYP genes do not have elevated rates of evolution, suggesting that these regions of the protein are not necessarily important in recognition of CYP2 substrates. Type II functional divergence analysis identified multiple residues in the active site of CYP2F, CYP2A, and CYP2B proteins that have undergone radical biochemical changes and may be functionally important.

Background: Increasing use of zebrafish in drug discovery and mechanistic toxicology demands knowledge of cytochrome P450 (CYP) gene regulation and function. CYP enzymes catalyze oxidative transformation leading to activation or inactivation of many endogenous and exogenous chemicals, with consequences for normal physiology and disease processes. Many CYPs potentially have roles in developmental specification, and many chemicals that cause developmental abnormalities are substrates for CYPs. Here we identify and annotate the full suite of CYP genes in zebrafish, compare these to the human CYP gene complement, and determine the expression of CYP genes during normal development., Results: Zebrafish have a total of 94 CYP genes, distributed among 18 gene families found also in mammals. There are 32 genes in CYP families 5 to 51, most of which are direct orthologs of human CYPs that are involved in endogenous functions including synthesis or inactivation of regulatory molecules. The high degree of sequence similarity suggests conservation of enzyme activities for these CYPs, confirmed in reports for some steroidogenic enzymes (e.g. CYP19, aromatase; CYP11A, P450scc; CYP17, steroid 17a-hydroxylase), and the CYP26 retinoic acid hydroxylases. Complexity is much greater in gene families 1, 2, and 3, which include CYPs prominent in metabolism of drugs and pollutants, as well as of endogenous substrates. There are orthologous relationships for some CYP1 s and some CYP3 s between zebrafish and human. In contrast, zebrafish have 47 CYP2 genes, compared to 16 in human, with only two (CYP2R1 and CYP2U1) recognized as orthologous based on sequence. Analysis of shared synteny identified CYP2 gene clusters evolutionarily related to mammalian CYP2 s, as well as unique clusters., Conclusions: Transcript profiling by microarray and quantitative PCR revealed that the majority of zebrafish CYP genes are expressed in embryos, with waves of expression of different sets of genes over the course of development. Transcripts of some CYP occur also in oocytes. The results provide a foundation for the use of zebrafish as a model in toxicological, pharmacological and chemical disease research.

We quantified mRNA abundance from 10 stages in the Giardia lamblia life cycle in vitro using Serial Analysis of Gene Expression (SAGE). 163 abundant transcripts were expressed constitutively. 71 transcripts were upregulated specifically during excystation and 42 during encystation. Nonetheless, the transcriptomes of cysts and trophozoites showed major differences. SAGE detected co-expressed clusters of 284 transcripts differentially expressed in cysts and excyzoites and 287 transcripts in vegetative trophozoites and encysting cells. All clusters included known genes and pathways as well as proteins unique to Giardia or diplomonads. SAGE analysis of the Giardia life cycle identified a number of kinases, phosphatases, and DNA replication proteins involved in excystation and encystation, which could be important for examining the roles of cell signaling in giardial differentiation. Overall, these data pave the way for directed gene discovery and a better understanding of the biology of G. lamblia., (Copyright 2010 Elsevier B.V. All rights reserved.)

Plasmodium belongs to the phylum Apicomplexa. Within the Apicomplexa, Plasmodium, Toxoplasma and Cryptosporidium are parasites of considerable medical importance while Theileria and Eimeria are animal pathogens. P. falciparum is particularly important as it causes malaria, resulting in more than 1 million deaths each year. The malaria parasite actively invades the host cell in which it propagates and several proteins associated with the apical organelles have been implicated to be crucial in the invasion process. The biogenesis of the apical organelles is not well understood, but several studies indicate that microtubule-based vesicular transport is involved. Vesicular transport proteins are also present in Plasmodium and are presumed to be involved in transcellular transport in infected erythrocytes. Dynein is a multi-subunit motor protein involved in microtubule-based vesicular transport. In this study, we analyzed the cytoplasmic dynein light chains (Dlcs) of P. falciparum since they provide adaptor surface to the cargoes and are likely to be involved in differential transport. Dlcs consist of three different families: TcTex1/2, LC8 and LC7/roadblock. The data presented demonstrate that P. falciparum Dlcs sequences and functional domains show high sequence similarity within the species, but that only the Dlc group 1 (LC8) has a high similarity to human orthologues. TcTex1 and LC7/roadblock have low similarity to human orthologues. This sequence variation could be targeted for vaccine or drug development.

Summary: Infection of the snail, Biomphalaria glabrata, by the free-swimming miracidial stage of the human blood fluke, Schistosoma mansoni, and its subsequent development to the parasitic sporocyst stage is critical to establishment of viable infections and continued human transmission. We performed a genome-wide expression analysis of the S. mansoni miracidia and developing sporocyst using Long Serial Analysis of Gene Expression (LongSAGE). Five cDNA libraries were constructed from miracidia and in vitro cultured 6- and 20-day-old sporocysts maintained in sporocyst medium (SM) or in SM conditioned by previous cultivation with cells of the B. glabrata embryonic (Bge) cell line. We generated 21 440 SAGE tags and mapped 13 381 to the S. mansoni gene predictions (v4.0e) either by estimating theoretical 3' UTR lengths or using existing 3' EST sequence data. Overall, 432 transcripts were found to be differentially expressed amongst all 5 libraries. In total, 172 tags were differentially expressed between miracidia and 6-day conditioned sporocysts and 152 were differentially expressed between miracidia and 6-day unconditioned sporocysts. In addition, 53 and 45 tags, respectively, were differentially expressed in 6-day and 20-day cultured sporocysts, due to the effects of exposure to Bge cell-conditioned medium.

Purpose: To evaluate the activity of imatinib in treating advanced, life-threatening malignancies expressing one or more imatinib-sensitive tyrosine kinases., Experimental Design: This was a phase II, open-label, single arm study. Patients > or = 15 years old with malignancies showing histologic or molecular evidence of expression/activation of imatinib-sensitive tyrosine kinases were enrolled. Patients were treated with 400 or 800 mg/d imatinib for hematologic malignancy and solid tumors, respectively. Treatment was continued until disease progression or unacceptable toxicity. The primary objective was to identify evidence of imatinib activity with tumor response as the primary end point., Results: One hundred eighty-six patients with 40 different malignancies were enrolled (78.5% solid tumors, 21.5% hematologic malignancies). Confirmed response occurred in 8.9% of solid tumor patients (4 complete, 9 partial) and 27.5% of hematologic malignancy patients (8 complete, 3 partial). Notable activity of imatinib was observed in only five tumor types (aggressive fibromatosis, dermatofibrosarcoma protuberans, hypereosinophilic syndrome, myeloproliferative disorders, and systemic mastocytosis). A total of 106 tumors were screened for activating mutations: five KIT mutations and no platelet-derived growth factor receptor mutations were found. One patient with systemic mastocytosis and a partial response to therapy had a novel imatinib-sensitive KIT mutation (D816T). There was no clear relationship between expression or activation of wild-type imatinib-sensitive tyrosine kinases and clinical response., Conclusion: Clinical benefit was largely confined to diseases with known genomic mechanisms of activation of imatinib target kinases. Our results indicate an important role for molecular characterization of tumors to identify patients likely to benefit from imatinib treatment.

Despite the widespread use of chemotherapy and other control strategies over the past 50years, transmission rates for schistosomiasis have changed little. Regardless of the approach used, future control efforts will require a more complete understanding of fundamental parasite biology. Schistosomes undergo complex development involving an alteration of parasite generations within a mammalian and freshwater molluscan host in the completion of its lifecycle. Little is known about factors controlling schistosome development, but understanding these processes may facilitate the discovery of new control methods. Therefore, our goal in this study is to determine global developmentally regulated and stage-specific gene expression in Schistosoma mansoni using serial analysis of gene expression (SAGE). We present a preliminary analysis of genes expressed during development and sexual differentiation in the mammalian host and during early larval development in the snail host. A number of novel, differentially expressed genes have been identified, both within and between the different developmental stages found in the mammalian and snail hosts.

The genome of the eukaryotic protist Giardia lamblia, an important human intestinal parasite, is compact in structure and content, contains few introns or mitochondrial relics, and has simplified machinery for DNA replication, transcription, RNA processing, and most metabolic pathways. Protein kinases comprise the single largest protein class and reflect Giardia's requirement for a complex signal transduction network for coordinating differentiation. Lateral gene transfer from bacterial and archaeal donors has shaped Giardia's genome, and previously unknown gene families, for example, cysteine-rich structural proteins, have been discovered. Unexpectedly, the genome shows little evidence of heterozygosity, supporting recent speculations that this organism is sexual. This genome sequence will not only be valuable for investigating the evolution of eukaryotes, but will also be applied to the search for new therapeutics for this parasite.

Coatomer coated (COPI) vesicles play a pivotal role for multiple membrane trafficking steps throughout the eukaryotic cell. Our focus is on betaCOP, one of the most well known components of the COPI multi-protein complex. Amino acid differences in betaCOP may dictate functional divergence across species during the course of evolution, especially with regards to the evolutionary pressures on obligate intracellular parasites. A bioinformatic analysis of betaCOP amino acid sequences was conducted for 49 eukaryotic species. Cloning and sequence analysis of the Toxoplasma gondii betaCOP homologue revealed several amino acid insertions unique to T. gondii and one C-terminal insertion that is unique to apicomplexan parasites. These findings led us to investigate the possibility that betaCOP experienced functional divergence during the course of its evolution. Bayesian phylogenetic analysis revealed a tree consistent with pan eukaryote distribution and long-branch lengths were observed among the apicomplexans. Further analysis revealed that kinetoplast betaCOP underwent the most amount of change, leading to perhaps an overall change of function. In comparison, T. gondii exhibited subtle yet specific amino acid changes. The amino acid substitutions did not occur in the same places as other lineages, suggesting that TgbetaCOP has a role specific to the apicomplexans. Our work identifies 48 residues that are likely to be functionally important when comparing apicomplexan, kinetoplastid, and fungal betaCOP.

Matrix proteins of mitochondria, hydrogenosomes and mitosomes are typically targeted and translocated into their respective organelles using N-terminal presequences that are subsequently cleaved by a peptidase. Here we characterize a approximately 47 kDa metallopeptidase, from the hydrogenosome-bearing, unicellular eukaryote Trichomonas vaginalis, that contains the active site motif (HXXEHX(76)E) characteristic of the beta subunit of the mitochondrial processing peptidase (MPP) and localizes to hydrogenosomes. The purified recombinant protein, named hydrogenosomal processing peptidase (HPP), is capable of cleaving a hydrogenosomal presequence in vitro, in contrast to MPP which requires both an alpha and beta subunit for activity. T. vaginalis HPP forms an approximately 100 kDa homodimer in vitro and also exists in an approximately 100 kDa complex in vivo. Our phylogenetic analyses support a common origin for HPP and betaMPP and demonstrate that gene duplication gave rise to alphaMPP and betaMPP before the divergence of T. vaginalis and mitochondria-bearing lineages. These data, together with published analyses of MPPs and putative mitosomal processing peptidases, lead us to propose that the length of targeting presequences and the subunit composition of organellar processing peptidases evolved in concert. Specifically, longer mitochondrial presequences may have evolved to require an alpha/beta heterodimer for accurate cleavage, while shorter hydrogenosomal and mitosomal presequences did not.

The ability of Giardia lamblia to undergo two distinct differentiations in response to physiologic stimuli is central to its pathogenesis. The giardial cytoskeleton changes drastically during encystation and excystation. However, the signal transduction pathways mediating these transformations are poorly understood. We tested the hypothesis that PP2A, a highly conserved serine/threonine protein phosphatase, might be important in giardial differentiation. We found that in vegetatively growing trophozoites, gPP2A-C protein localizes to basal bodies/centrosomes, and to cytoskeletal structures unique to Giardia: the ventral disk, and the dense rods of the anterior, posterior-lateral, and caudal flagella. During encystation, gPP2A-C protein disappears from only the anterior flagellar dense rods. During excystation, gPP2A-C localizes to the cyst wall in excysting cysts but is not found in the wall of cysts with emerging excyzoites. Transcriptome and immunoblot analyses indicated that gPP2A-C mRNA and protein are upregulated in mature cysts and during the early stage of excystation that models passage through the host stomach. Stable expression of gPP2A-C antisense RNA did not affect vegetative growth, but strongly inhibited the formation of encystation secretory vesicles (ESV) and water-resistant cysts. Moreover, the few cysts that formed were highly defective in excystation. Thus, gPP2A-C localizes to universal cytoskeletal structures and to structures unique to Giardia. It is also important for encystation and excystation, crucial giardial transformations that entail entry into and exit from dormancy.

Since the Giardia lamblia cyst wall is necessary for survival in the environment and host infection, we tested the hypothesis that it contains proteins other than the three known cyst wall proteins. Serial analysis of gene expression during growth and encystation revealed a gene, "HCNCp" (High Cysteine Non-variant Cyst protein), that was upregulated late in encystation, and that resembled the classic Giardia variable surface proteins (VSPs) that cover the trophozoite plasmalemma. HCNCp is 13.9% cysteine, with many "CxxC" tetrapeptide motifs and a transmembrane sequence near the C-terminus. However, HCNCp has multiple "CxC" motifs rarely found in VSPs, and does not localize to the trophozoite plasmalemma. Moreover, the HCNCp C-terminus differed from the canonical VSP signature. Full-length epitope-tagged HCNCp expressed under its own promoter was upregulated during encystation with highest expression in cysts, including 42 and 21 kDa C-terminal fragments. Tagged HCNCp targeted to the nuclear envelope in trophozoites, and co-localized with cyst proteins to encystation-specific secretory vesicles during encystation. HCNCp defined a novel trafficking pathway as it localized to the wall and body of cysts, while the cyst proteins were exclusively in the wall. Unlike VSPs, HCNCp is expressed in at least five giardial strains and four WB subclones expressing different VSPs. Bioinformatics identified 60 additional large high cysteine membrane proteins (HCMp) containing > or = 20 CxxC/CxC's lacking the VSP-specific C-terminal CRGKA. HCMp were absent or rare in other model or parasite genomes, except for Tetrahymena thermophila with 30. MEME analysis classified the 61 gHCMp genes into nine groups with similar internal motifs. Our data suggest that HCNCp is a novel invariant cyst protein belonging to a new HCMp family that is abundant in the Giardia genome. HCNCp and the other HCMp provide a rich source for developing parasite-specific diagnostic reagents, vaccine candidates, and subjects for further research into Giardia biology.

Trichomonas vaginalis is a unicellular eukaryote that lacks mitochondria and contains a specialized organelle, the hydrogenosome, involved in carbohydrate metabolism and iron-sulfur cluster assembly. We report the identification of two glycine cleavage H proteins and a dihydrolipoamide dehydrogenase (L protein) of the glycine decarboxylase complex in T. vaginalis with predicted N-terminal hydrogenosomal presequences. Immunofluorescence analyses reveal that both H and L proteins are localized in hydrogenosomes, providing the first evidence for amino acid metabolism in this organelle. All three proteins were expressed in Escherichia coli and purified to homogeneity. The experimental Km of L protein for the two H proteins were 2.6 microM and 3.7 microM, consistent with both H proteins serving as substrates of L protein. Analyses using purified hydrogenosomes showed that endogenous H proteins exist as monomers and endogenous L protein as a homodimer in their native states. Phylogenetic analyses of L proteins revealed that the T. vaginalis homologue shares a common ancestry with dihydrolipoamide dehydrogenases from the firmicute bacteria, indicating its acquisition via a horizontal gene transfer event independent of the origins of mitochondria and hydrogenosomes.

The host range of African trypanosomes is influenced by innate protective molecules in the blood of primates. A subfraction of human high-density lipoprotein (HDL) containing apolipoprotein A-I, apolipoprotein L-I, and haptoglobin-related protein is toxic to Trypanosoma brucei brucei but not the human sleeping sickness parasite Trypanosoma brucei rhodesiense. It is thought that T. b. rhodesiense evolved from a T. b. brucei-like ancestor and expresses a defense protein that ablates the antitrypanosomal activity of human HDL. To directly investigate this possibility, we developed an in vitro selection to generate human HDL-resistant T. b. brucei. Here we show that conversion of T. b. brucei from human HDL sensitive to resistant correlates with changes in the expression of the variant surface glycoprotein (VSG) and abolished uptake of the cytotoxic human HDLs. Complete transcriptome analysis of the HDL-susceptible and -resistant trypanosomes confirmed that VSG switching had occurred but failed to reveal the expression of other genes specifically associated with human HDL resistance, including the serum resistance-associated gene (SRA) of T. b. rhodesiense. In addition, we found that while the original active expression site was still utilized, expression of three expression site-associated genes (ESAG) was altered in the HDL-resistant trypanosomes. These findings demonstrate that resistance to human HDLs can be acquired by T. b. brucei.

The abundant and widespread coccolithophore Emiliania huxleyi plays an important role in mediating CO2 exchange between the ocean and the atmosphere through its impact on marine photosynthesis and calcification. Here, we use long serial analysis of gene expression (SAGE) to identify E. huxleyi genes responsive to nitrogen (N) or phosphorus (P) starvation. Long SAGE is an elegant approach for examining quantitative and comprehensive gene expression patterns without a priori knowledge of gene sequences via the detection of 21-bp nucleotide sequence tags. E. huxleyi appears to have a robust transcriptional-level response to macronutrient deficiency, with 42 tags uniquely present or up-regulated twofold or greater in the N-starved library and 128 tags uniquely present or up-regulated twofold or greater in the P-starved library. The expression patterns of several tags were validated with reverse transcriptase PCR. Roughly 48% of these differentially expressed tags could be mapped to publicly available genomic or expressed sequence tag (EST) sequence data. For example, in the P-starved library a number of the tags mapped to genes with a role in P scavenging, including a putative phosphate-repressible permease and a putative polyphosphate synthetase. In short, the long SAGE analyses have (i) identified many new differentially regulated gene sequences, (ii) assigned regulation data to EST sequences with no database homology and unknown function, and (iii) highlighted previously uncharacterized aspects of E. huxleyi N and P physiology. To this end, our long SAGE libraries provide a new public resource for gene discovery and transcriptional analysis in this biogeochemically important marine organism.

Giardia lamblia is a protozoan parasite infecting the upper mammalian small intestine. Infection relies upon the ability of the parasite to attach to the intestine via a unique cytoskeletal organelle, the ventral disk. We determined the composition and structure of the disk throughout the life cycle of the parasite and identified a new disk protein, SALP-1. SALP-1 is an immunodominant protein related to striated fiber-assemblin (SFA). The disk is disassembled during encystation and stored as four fragments in the immobile cyst. Serial Analysis of Gene Expression (SAGE) showed that the mRNA levels of the disk proteins decreased in encystation but two-dimensional protein gels showed that the protein levels were more constant. The parasite emerges without a functional disk but the four disk fragments are quickly reassembled into two new disks on the dividing, early excysting form. Thus, disk proteins are stored within the cyst, ready to be used in the rapid steps of excystation.

Through a genome survey and phylogenetic analysis, we have identified and sequenced 14 new coding regions for alpha-giardins in Giardia lamblia. These proteins are related to annexins and comprise a multi-gene family with 21 members. Many alpha giardins are highly expressed proteins that are very immunogenic during acute giardiasis in humans. However, little is known about the function of these proteins. By using PCR with different combinations of gene-specific primers, we demonstrated that several of the genes localised to the same chromosomal fragment. These data point towards a molecular evolution through gene duplication and subsequent functional divergence. Semi-quantitative reverse transcriptase-PCR analysis of the Giardia life cycle revealed large differences in mRNA expression levels of the alpha giardins. Epitope tagging of the alpha-giardins localised them to different cytoskeletal components, such as the flagella and the adhesive disc, but also to the plasma membrane. These localisation experiments suggest alpha-giardins play a role in cell motility, attachment and membrane stability.

Aromatase (cytochrome P450 19, CYP19, P450arom) is the enzyme responsible for the production of estrogens, hormones critical for development and reproduction. Aromatase was sequenced from a white-sided dolphin (Lagenorhynchus acutus) ovary, transiently transfected into HEK 293 cells, and the expressed protein was characterized for aromatase activity in the presence of androstenedione and testosterone and after exposure to the aromatase inhibitor letrazole. The Kms for androstenedione and testosterone were 63.5 and 75 nM, respectively, values that are very similar to those reported for other mammalian aromatases. A Bayesian phylogenetic analysis of the vertebrate aromatases was performed on the amino acid sequences of aromatases from fish, amphibians, reptiles, birds, and mammals. Based on known species phylogeny, the cetacean aromatase showed an expected grouping with artiodactyls (cow, sheep, and goat). An analysis of functional divergence showed strong conservation of aromatase across the entire protein, which indicates that the observed sequence divergence is functionally neutral.

The Giardia lamblia genome sequencing project affords us a unique opportunity to conduct comparative analyses of core cellular systems between early and late-diverging eukaryotes on a genome-wide scale. We report a survey to identify canonical transcription components in Giardia, focusing on RNA polymerase (RNAP) subunits and transcription-initiation factors. Our survey revealed that Giardia contains homologs to 21 of the 28 polypeptides comprising eukaryal RNAPI, RNAPII, and RNAPIII; six of the seven RNAP subunits without giardial homologs are polymerase specific. Components of only four of the 12 general transcription initiation factors have giardial homologs. Surprisingly, giardial TATA-binding protein (TBP) is highly divergent with respect to archaeal and higher eukaryotic TBPs, and a giardial homolog of transcription factor IIB was not identified. We conclude that Giardia represents a transition during the evolution of eukaryal transcription systems, exhibiting a relatively complete set of RNAP subunits and a rudimentary basal initiation apparatus for each transcription system. Most class-specific RNAP subunits and basal initiation factors appear to have evolved after the divergence of Giardia from the main eukaryotic line of descent. Consequently, Giardia is predicted to be unique in many aspects of transcription initiation with respect to paradigms derived from studies in crown eukaryotes., (Copyright 2004 Cold Spring Harbor Laboratory Press ISSN)

Cytochrome P450 genes (CYP) constitute a superfamily with members known from the Bacteria, Archaea, and Eukarya. The CYP3 gene family includes the CYP3A and CYP3B subfamilies. Members of the CYP3A subfamily represent the dominant CYP forms expressed in the digestive and respiratory tracts of vertebrates. The CYP3A enzymes metabolize a wide variety of chemically diverse lipophilic organic compounds. To understand vertebrate CYP3 diversity better, we determined the killifish (Fundulus heteroclitus) CYP3A30 and CYP3A56 and the ball python (Python regius) CYP3A42 sequences. We performed phylogenetic analyses of 45 vertebrate CYP3 amino acid sequences using a Bayesian approach. Our analyses indicate that teleost, diapsid, and mammalian CYP3A genes have undergone independent diversification and that the ancestral vertebrate genome contained a single CYP3A gene. Most CYP3A diversity is the product of recent gene duplication events. There is strong support for placement of the guinea pig CYP3A genes within the rodent CYP3A diversification. The rat, mouse, and hamster CYP3A genes are mixed among several rodent CYP3A subclades, indicative of a complex history involving speciation and gene duplication.

Giardia lamblia is an early branching eukaryote, and although distinctly eukaryotic in its cell and molecular biology, transcription and translation in G. lamblia demonstrate important differences from these processes in higher eukaryotes. The cyclic octapeptide amanitin is a relatively selective inhibitor of eukaryotic RNA polymerase II (RNAP II) and is commonly used to study RNAP II transcription. Therefore, we measured the sensitivity of G. lamblia RNAP II transcription to alpha-amanitin and found that unlike most other eukaryotes, RNAP II transcription in Giardia is resistant to 1 mg/ml amanitin. In contrast, 50 microg/ml amanitin inhibits 85% of RNAP III transcription activity using leucyl-tRNA as a template. To better understand transcription in G. lamblia, we identified 10 of the 12 known eukaryotic rpb subunits, including all 10 subunits that are required for viability in Saccharomyces cerevisiae. The amanitin motif (amanitin binding site) of Rpb1 from G. lamblia has amino acid substitutions at six highly conserved sites that have been associated with amanitin resistance in other organisms. These observations of amanitin resistance of Giardia RNA polymerase II support previous proposals of the mechanism of amanitin resistance in other organisms and provide a molecular framework for the development of novel drugs with selective activity against G. lamblia.

Entamoeba histolytica and Spironucleus barkhanus have genes that encode short iron-dependent hydrogenases (Fe-hydrogenases), even though these protists lack hydrogenosomes. To understand better the biochemistry of the protist Fe-hydrogenases, we prepared a recombinant E. histolytica short Fe-hydrogenase and measured its activity in vitro. A Giardia lamblia gene encoding a short Fe-hydrogenase was identified from shotgun genomic sequences, and RT-PCR showed that cultured entamoebas and giardias transcribe short Fe-hydrogenase mRNAs. A second E. histolytica gene, which encoded a long Fe-hydrogenase, was identified from shotgun genomic sequences. Phylogenetic analyses suggested that the short Fe-hydrogenase genes of entamoeba and diplomonads share a common ancestor, while the long Fe-hydrogenase gene of entamoeba appears to have been laterally transferred from a bacterium. These results are discussed in the context of competing ideas for the origins of genes encoding fermentation enzymes of these protists.

Giardia lamblia is an important human intestinal parasite that survives outside of the host by differentiation of trophozoites into infectious cysts. Transcriptional regulation is key for encystation gene expression, but the mechanisms are unknown. Giardia genome database searches identified a myb-like gene (gmyb2) whose expression increased during encystation. Epitope-tagged gMyb2 localized to both nuclei. DNA binding and mutation analysis showed that gMyb2 binds specifically to C(T/A)ACAG, a c-Myb-like target sequence in the promoters of encystation-induced genes encoding gMyb2, three cyst wall proteins and G6PI-B, a key enzyme in cyst wall polysaccharide biosynthesis. gMyb2 binding sites were not found in the upstream regions of 31 other giardial genes. Deletion of the putative gMyb2 binding site greatly reduced encystation-specific promoter activity of g6pi-b. Fusion of gMyb2 binding sites to the constitutive ran promoter or g6pi-b promoter deletion lacking the gMyb2 binding site in-duced encystation-specific expression. gMyb2 may play an important role in transcriptional regulation of encystation genes, and may help co-ordinate synthesis of cyst wall proteins and polysaccharide. gMyb2 is the first giardial transcription factor to be functionally identified and the first that is associated with upregulation of encystation genes. This work provides a model for study of differential gene regulation in early diverging eukaryotic organisms.

Giardia lamblia and Entamoeba histolytica are amitochondriate, microaerophilic protists which use fermentation enzymes like those of bacteria to survive anaerobic conditions within the intestinal lumen. Genes encoding fermentation enzymes and related electron transport peptides (e.g., ferredoxins) in giardia organisms and amebae are hypothesized to be derived from either an ancient anaerobic eukaryote (amitochondriate fossil hypothesis), a mitochondrial endosymbiont (hydrogen hypothesis), or anaerobic bacteria (lateral transfer hypothesis). The goals here were to complete the molecular characterization of giardial and amebic fermentation enzymes and to determine the origins of the genes encoding them, when possible. A putative giardia [2Fe-2S]ferredoxin which had a hypothetical organelle-targeting sequence at its N terminus showed similarity to mitochondrial ferredoxins and the hydrogenosomal ferredoxin of Trichomonas vaginalis (another luminal protist). However, phylogenetic trees were star shaped, with weak bootstrap support, so we were unable to confirm or rule out the endosymbiotic origin of the giardia [2Fe-2S]ferredoxin gene. Putative giardial and amebic 6-kDa ferredoxins, ferredoxin-nitroreductase fusion proteins, and oxygen-insensitive nitroreductases each tentatively supported the lateral transfer hypothesis. Although there were not enough sequences to perform meaningful phylogenetic analyses, the unique common occurrence of these peptides and enzymes in giardia organisms, amebae, and the few anaerobic prokaryotes suggests the possibility of lateral transfer. In contrast, there was more robust phylogenetic evidence for the lateral transfer of G. lamblia genes encoding an NADH oxidase from a gram-positive coccus and a microbial group 3 alcohol dehydrogenase from thermoanaerobic prokaryotes. In further support of lateral transfer, the G. lamblia NADH oxidase and adh3 genes appeared to have an evolutionary history distinct from those of E. histolytica.

Short introns occur in numerous protist lineages, but there are no reports of intervening sequences in the protists Giardia lamblia and Trichomonas vaginalis, which may represent the deepest known branches in the eukaryotic line of descent. We have discovered a 35-bp spliceosomal intron in a gene encoding a putative [2Fe-2S] ferredoxin of G. lamblia. The Giardia intron contains a canonical splice site at its 3' end (AG), a noncanonical splice site at its 5' end (CT), and a branch point sequence that fits the yeast consensus sequence of TACTAAC except for the first nucleotide (AACTAAC). We have also identified several G. lamblia genes with spliceosomal peptides, including homologues of eukaryote-specific spliceosomal peptides (Prp8 and Prp11), several DExH-box RNA-helicases that have homologues in eubacteria, but serve essential functions in the splicing of introns in eukaryotes, and 11 predicted archaebacteria-like Sm and like-Sm core peptides, which coat small nuclear RNAs. Phylogenetic analyses show the Giardia Sm core peptides are the products of multiple, ancestral gene duplications followed by divergence, but they retain strong similarity to Sm and like-Sm peptides of other eukaryotes. Although we have documented only a single intron in Giardia, it likely has other introns and fully functional, spliceosomal machinery. If introns were added during eukaryotic evolution (the introns-late hypothesis), then these results push back the date of this event before the branching of G. lamblia.

A phylogenetic analysis of protein disulfide isomerase (PDI) domain evolution was performed with the inclusion of recently reported PDIs from the amitochondriate protist Giardia lamblia, yeast PDIs that contain a single thioredoxin-like domain, and PDIs from a diverse selection of protists. We additionally report and include two new giardial PDIs, each with a single thioredoxin-like domain. Inclusion of protist PDIs in our analyses revealed that the evolutionary history of the endoplasmic reticulum may not be simple. Phylogenetic analyses support common ancestry of all eukaryotic PDIs from a thioredoxin ancestor and independent duplications of thioredoxin-like domains within PDIs throughout eukaryote evolution. This was particularly evident for Acanthamoeba PDI, Dictyostelium PDI, and mammalian erp5 domains. In contrast, gene duplication, instead of domain duplication, produces PDI diversity in G. lamblia. Based on our results and the known diversity of PDIs, we present a new hypothesis that the five single-domain PDIs of G. lamblia may reflect an ancestral mechanism of protein folding in the eukaryotic endoplasmic reticulum. The PDI complement of G. lamblia and yeast suggests that a combination of PDIs may be used as a redox chain analogous to that known for bacterial Dsb proteins.

Molecular and morphological evidence points to the ancyromonad Ancyromonas as a plausible candidate for the closest relative to the common ancestor of metazoans, fungi, and choanoflagellates (the Opisthokonta). Using 18S rDNA sequences from most of the major eukaryotic lineages, maximum-likelihood, minimum-evolution, and maximum-parsimony analyses yielded congruent phylogenies supporting this hypothesis. Combined with ultrastructural similarities between Ancyromonas and opisthokonts, the evidence presented here suggests that Ancyromonas may form an independent lineage, the Ancyromonadida Cavalier-Smith 1997, closer in its relationship to the opisthokonts than is its nearest protist relatives, the Apusomonadida. However, the very low bootstrap support for deep nodes and hypothesis testing indicate that the resolving power of 18S rDNA sequences is limited for examining this aspect of eukaryotic phylogeny. Alternate branching positions for the Ancyromonas lineage cannot be robustly rejected, revealing the importance of ultrastructure when examining the origins of multicellularity. The future use of a multigene approach may additionally be needed to resolve this aspect of eukaryotic phylogeny.

The Giardia genome project database provides an online resource for Giardia lamblia (WB strain, clone C6) genome sequence information. The database includes edited single-pass reads, the results of BLASTX searches, and details of progress towards sequencing the entire 12 million-bp Giardia genome. Pre-sorted BLASTX results can be retrieved based on keyword searches and BLAST searches of the high throughput Giardia data can be initiated from the web site or through NCBI. Descriptions of the genomic DNA libraries, project protocols and summary statistics are also available. Although the Giardia genome project is ongoing, new sequences are made available on a bi-monthly basis to ensure that researchers have access to information that may assist them in the search for genes and their biological function. The current URL of the Giardia genome project database is www.mbl.edu/Giardia.