Efficient Genome Editing in Clostridium cellulolyticum via CRISPR-Cas9 Nickase.

The CRISPR-Cas9 system is a powerful and revolutionary genome-editing tool for eukaryotic genomes, but its use in bacterial genomes is very limited. Here, we investigated the use of the Streptococcus pyogenes CRISPR-Cas9 system in editing the genome of Clostridium cellulolyticum, a model microorganism for bioenergy research. Wild-type Cas9-induced double-strand breaks were lethal to C. cellulolyticum due to the minimal expression of nonhomologous end joining (NHEJ) components in this strain. To circumvent this lethality, Cas9 nickase was applied to develop a single-nick-triggered homologous recombination strategy, which allows precise one-step editing at intended genomic loci by transforming a single vector. This strategy has a high editing efficiency (>95%) even using short homologous arms (0.2 kb), is able to deliver foreign genes into the genome in a single step without a marker, enables precise editing even at two very similar target sites differing by two bases preceding the seed region, and has a very high target site density (median interval distance of 9 bp and 95.7% gene coverage in C. cellulolyticum). Together, these results establish a simple and robust methodology for genome editing in NHEJ-ineffective prokaryotes.

Continuous cellulosic bioethanol fermentation by cyclic fed-batch cocultivation.

Cocultivation of cellulolytic and saccharolytic microbial populations is a promising strategy to improve bioethanol production from the fermentation of recalcitrant cellulosic materials. Earlier studies have demonstrated the effectiveness of cocultivation in enhancing ethanolic fermentation of cellulose in batch fermentation. To further enhance process efficiency, a semicontinuous cyclic fed-batch fermentor configuration was evaluated for its potential in enhancing the efficiency of cellulose fermentation using cocultivation. Cocultures of cellulolytic Clostridium thermocellum LQRI and saccharolytic Thermoanaerobacter pseudethanolicus strain X514 were tested in the semicontinuous fermentor as a model system. Initial cellulose concentration and pH were identified as the key process parameters controlling cellulose fermentation performance in the fixed-volume cyclic fed-batch coculture system. At an initial cellulose concentration of 40 g liter(-1), the concentration of ethanol produced with pH control was 4.5-fold higher than that without pH control. It was also found that efficient cellulosic bioethanol production by cocultivation was sustained in the semicontinuous configuration, with bioethanol production reaching 474 mM in 96 h with an initial cellulose concentration of 80 g liter(-1) and pH controlled at 6.5 to 6.8. These results suggested the advantages of the cyclic fed-batch process for cellulosic bioethanol fermentation by the cocultures.

StressChip as a high-throughput tool for assessing microbial community responses to environmental stresses.

Microbial community responses to environmental stresses are critical for microbial growth, survival, and adaptation. To fill major gaps in our ability to discern the influence of environmental changes on microbial communities from engineered and natural environments, a functional gene-based microarray, termed StressChip, has been developed. First, 46 functional genes involved in microbial responses to environmental stresses such as changes to temperature, osmolarity, oxidative status, nutrient limitation, or general stress response were selected and curated. A total of 22,855 probes were designed, covering 79,628 coding sequences from 985 bacterial, 76 archaeal, and 59 eukaryotic species/strains. Probe specificity was computationally verified. Second, the usefulness of functional genes as indicators of stress response was examined by surveying their distribution in metagenome data sets. The abundance of individual stress response genes is consistent with expected distributions based on respective habitats. Third, the StressChip was used to analyze marine microbial communities from the Deepwater Horizon oil spill. That functional stress response genes were detected in higher abundance (p < 0.05) in oil plume compared to nonplume samples indicated shifts in community composition and structure, consistent with previous results. In summary, StressChip provides a new tool for accessing microbial community functional structure and responses to environmental changes.

Developing virtual and augmented reality applications for science, technology, engineering and math education.

The Rhode Island IDeA Network of Biomedical Research Excellence Molecular Informatics Core at the University of Rhode Island Information Technology Services Innovative Learning Technologies developed virtual and augmented reality applications to teach concepts in biomedical science, including pharmacology, medicinal chemistry, cell culture and nanotechnology. The apps were developed as full virtual reality/augmented reality and 3D gaming versions, which do not require virtual reality headsets. Development challenges included creating intuitive user interfaces, text-to-voice functionality, visualization of molecules and implementing complex science concepts. In-app quizzes are used to assess the user's understanding of topics, and user feedback was collected for several apps to improve the experience. The apps were positively reviewed by users and are being implemented into the curriculum at the University of Rhode Island.

Genome-wide analysis of acetivibrio cellulolyticus provides a blueprint of an elaborate cellulosome system.

BACKGROUND: Microbial degradation of plant cell walls and its conversion to sugars and other byproducts is a key step in the carbon cycle on Earth. In order to process heterogeneous plant-derived biomass, specialized anaerobic bacteria use an elaborate multi-enzyme cellulosome complex to synergistically deconstruct cellulosic substrates. The cellulosome was first discovered in the cellulolytic thermophile, Clostridium thermocellum, and much of our knowledge of this intriguing type of protein composite is based on the cellulosome of this environmentally and biotechnologically important bacterium. The recently sequenced genome of the cellulolytic mesophile, Acetivibrio cellulolyticus, allows detailed comparison of the cellulosomes of these two select cellulosome-producing bacteria. RESULTS: Comprehensive analysis of the A. cellulolyticus draft genome sequence revealed a very sophisticated cellulosome system. Compared to C. thermocellum, the cellulosomal architecture of A. cellulolyticus is much more extensive, whereby the genome encodes for twice the number of cohesin- and dockerin-containing proteins. The A. cellulolyticus genome has thus evolved an inflated number of 143 dockerin-containing genes, coding for multimodular proteins with distinctive catalytic and carbohydrate-binding modules that play critical roles in biomass degradation. Additionally, 41 putative cohesin modules distributed in 16 different scaffoldin proteins were identified in the genome, representing a broader diversity and modularity than those of Clostridium thermocellum. Although many of the A. cellulolyticus scaffoldins appear in unconventional modular combinations, elements of the basic structural scaffoldins are maintained in both species. In addition, both species exhibit similarly elaborate cell-anchoring and cellulosome-related gene- regulatory elements. CONCLUSIONS: This work portrays a particularly intricate, cell-surface cellulosome system in A. cellulolyticus and provides a blueprint for examining the specific roles of the various cellulosomal components in the degradation of complex carbohydrate substrates of the plant cell wall by the bacterium.

Functional characterization of Crp/Fnr-type global transcriptional regulators in Desulfovibrio vulgaris Hildenborough.

Crp/Fnr-type global transcriptional regulators regulate various metabolic pathways in bacteria and typically function in response to environmental changes. However, little is known about the function of four annotated Crp/Fnr homologs (DVU0379, DVU2097, DVU2547, and DVU3111) in Desulfovibrio vulgaris Hildenborough. A systematic study using bioinformatic, transcriptomic, genetic, and physiological approaches was conducted to characterize their roles in stress responses. Similar growth phenotypes were observed for the crp/fnr deletion mutants under multiple stress conditions. Nevertheless, the idea of distinct functions of Crp/Fnr-type regulators in stress responses was supported by phylogeny, gene transcription changes, fitness changes, and physiological differences. The four D. vulgaris Crp/Fnr homologs are localized in three subfamilies (HcpR, CooA, and cc). The crp/fnr knockout mutants were well separated by transcriptional profiling using detrended correspondence analysis (DCA), and more genes significantly changed in expression in a ΔDVU3111 mutant (JW9013) than in the other three paralogs. In fitness studies, strain JW9013 showed the lowest fitness under standard growth conditions (i.e., sulfate reduction) and the highest fitness under NaCl or chromate stress conditions; better fitness was observed for a ΔDVU2547 mutant (JW9011) under nitrite stress conditions and a ΔDVU2097 mutant (JW9009) under air stress conditions. A higher Cr(VI) reduction rate was observed for strain JW9013 in experiments with washed cells. These results suggested that the four Crp/Fnr-type global regulators play distinct roles in stress responses of D. vulgaris. DVU3111 is implicated in responses to NaCl and chromate stresses, DVU2547 in nitrite stress responses, and DVU2097 in air stress responses.

Isolation and characterization of Shigella flexneri G3, capable of effective cellulosic saccharification under mesophilic conditions.

A novel Shigella strain (Shigella flexneri G3) showing high cellulolytic activity under mesophilic, anaerobic conditions was isolated and characterized. The bacterium is Gram negative, short rod shaped, and nonmotile and displays effective production of glucose, cellobiose, and other oligosaccharides from cellulose (Avicel PH-101) under optimal conditions (40°C and pH 6.5). Approximately 75% of the cellulose was hydrolyzed in modified ATCC 1191 medium containing 0.3% cellulose, and the oligosaccharide production yield and specific production rate reached 375 mg g Avicel(-1) and 6.25 mg g Avicel(-1) h(-1), respectively, after a 60-hour incubation. To our knowledge, this represents the highest oligosaccharide yield and specific rate from cellulose for mesophilic bacterial monocultures reported so far. The results demonstrate that S. flexneri G3 is capable of rapid conversion of cellulose to oligosaccharides, with potential biofuel applications under mesophilic conditions.

Correlation of genomic and physiological traits of thermoanaerobacter species with biofuel yields.

Thermophilic anaerobic noncellulolytic Thermoanaerobacter species are of great biotechnological importance in cellulosic ethanol production due to their ability to produce high ethanol yields by simultaneous fermentation of hexose and pentose. Understanding the genome structure of these species is critical to improving and implementing these bacteria for possible biotechnological use in consolidated bioprocessing schemes (CBP) for cellulosic ethanol production. Here we describe a comparative genome analysis of two ethanologenic bacteria, Thermoanaerobacter sp. X514 and Thermoanaerobacter pseudethanolicus 39E. Compared to 39E, X514 has several unique key characteristics important to cellulosic biotechnology, including additional alcohol dehydrogenases and xylose transporters, modifications to pentose metabolism, and a complete vitamin B12 biosynthesis pathway. Experimental results from growth, metabolic flux, and microarray gene expression analyses support genome sequencing-based predictions which help to explain the distinct differences in ethanol production between these strains. The availability of whole-genome sequence and comparative genomic analyses will aid in engineering and optimizing Thermoanaerobacter strains for viable CBP strategies.

RI-INBRE: A Statewide NIH Program Grant to Improve Institutional Biomedical Research Capacity in Rhode Island.

The overarching goal of the Rhode Island-IDeA Network of Biomedical Research Excellence (RI-INBRE) is to improve institutional capacity for biomedical research excellence and expand student experiential training opportunities in the State of Rhode Island. RI-INBRE comprises five major core components: The Administrative Core, the Bioinformatics Core, the Centralized Research Core Facility, the Training Core, and the Developmental Research Project Program Core. Since its inception in 2001, RI-INBRE has made significant investments and marked advancements in the biomedical research infrastructure of Rhode Island. RI-INBRE funding has increased the scale and quality of faculty research and engaged undergraduate students, graduate students, and postdoctoral fellows in structured and mentored research training experiences. Over the last 19 years, RI-INBRE has supported 212 faculty researchers and over 533 projects and has provided research-training opportunities for nearly 2,000 students, resulting in 757 publications. Through its student-training program, RI-INBRE has contributed to regional workforce development by engaging students and encouraging them to pursue careers in biomedical fields. Many of these students have been admitted to graduate or medical schools and obtained biomedical industry jobs following graduation. RI-INBRE has been particularly influential in building the research infrastructure at primarily undergraduate institutions, which have seen significant improvements in research quality and output, student training, and research infrastructure.

Sequencing of multiple clostridial genomes related to biomass conversion and biofuel production.

Modern methods to develop microbe-based biomass conversion processes require a system-level understanding of the microbes involved. Clostridium species have long been recognized as ideal candidates for processes involving biomass conversion and production of various biofuels and other industrial products. To expand the knowledge base for clostridial species relevant to current biofuel production efforts, we have sequenced the genomes of 20 species spanning multiple genera. The majority of species sequenced fall within the class III cellulosome-encoding Clostridium and the class V saccharolytic Thermoanaerobacteraceae. Species were chosen based on representation in the experimental literature as model organisms, ability to degrade cellulosic biomass either by free enzymes or by cellulosomes, ability to rapidly ferment hexose and pentose sugars to ethanol, and ability to ferment synthesis gas to ethanol. The sequenced strains significantly increase the number of noncommensal/nonpathogenic clostridial species and provide a key foundation for future studies of biomass conversion, cellulosome composition, and clostridial systems biology.

Hydrogen peroxide-induced oxidative stress responses in Desulfovibrio vulgaris Hildenborough.

To understand how sulphate-reducing bacteria respond to oxidative stresses, the responses of Desulfovibrio vulgaris Hildenborough to H(2)O(2)-induced stresses were investigated with transcriptomic, proteomic and genetic approaches. H(2)O(2) and induced chemical species (e.g. polysulfide, ROS) and redox potential shift increased the expressions of the genes involved in detoxification, thioredoxin-dependent reduction system, protein and DNA repair, and decreased those involved in sulfate reduction, lactate oxidation and protein synthesis. A gene coexpression network analysis revealed complicated network interactions among differentially expressed genes, and suggested possible importance of several hypothetical genes in H(2)O(2) stress. Also, most of the genes in PerR and Fur regulons were highly induced, and the abundance of a Fur regulon protein increased. Mutant analysis suggested that PerR and Fur are functionally overlapped in response to stresses induced by H(2)O(2) and reaction products, and the upregulation of thioredoxin-dependent reduction genes was independent of PerR or Fur. It appears that induction of those stress response genes could contribute to the increased resistance of deletion mutants to H(2)O(2)-induced stresses. In addition, a conceptual cellular model of D. vulgaris responses to H(2)O(2) stress was constructed to illustrate that this bacterium may employ a complicated molecular mechanism to defend against the H(2)O(2)-induced stresses.

Succession of the bacterial community and dynamics of hydrogen producers in a hydrogen-producing bioreactor.

Variation in the hydrogen production rate was consistent with the succession of dominant bacteria during the batch fermentation process. Denaturing gradient gel electrophoresis (DGGE) of 16S rRNA genes and quantitative analysis of the hydA genes at both the DNA and mRNA levels confirmed that Clostridium perfringens was the most dominant hydrogen producer in the bioreactor.

Global transcriptional, physiological, and metabolite analyses of the responses of Desulfovibrio vulgaris hildenborough to salt adaptation.

The response of Desulfovibrio vulgaris Hildenborough to salt adaptation (long-term NaCl exposure) was examined by performing physiological, global transcriptional, and metabolite analyses. Salt adaptation was reflected by increased expression of genes involved in amino acid biosynthesis and transport, electron transfer, hydrogen oxidation, and general stress responses (e.g., heat shock proteins, phage shock proteins, and oxidative stress response proteins). The expression of genes involved in carbon metabolism, cell growth, and phage structures was decreased. Transcriptome profiles of D. vulgaris responses to salt adaptation were compared with transcriptome profiles of D. vulgaris responses to salt shock (short-term NaCl exposure). Metabolite assays showed that glutamate and alanine accumulated under salt adaptation conditions, suggesting that these amino acids may be used as osmoprotectants in D. vulgaris. Addition of amino acids (glutamate, alanine, and tryptophan) or yeast extract to the growth medium relieved salt-related growth inhibition. A conceptual model that links the observed results to currently available knowledge is proposed to increase our understanding of the mechanisms of D. vulgaris adaptation to elevated NaCl levels.

Fenofibrate Improves Liver Function and Reduces the Toxicity of the Bile Acid Pool in Patients With Primary Biliary Cholangitis and Primary Sclerosing Cholangitis Who Are Partial Responders to Ursodiol.

Cholestatic liver diseases result in the hepatic retention of bile acids, causing subsequent liver toxicity. Peroxisome proliferator-activated receptor alpha (PPARα) regulates bile acid metabolism. In this retrospective observational study, we assessed the effects of fenofibrate (a PPARα agonist) therapy on bile acid metabolism when given to patients with primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC) who have had an incomplete response to Ursodiol monotherapy. When fenofibrate was added to Ursodiol therapy there was a significant reduction and in some cases normalization of serum alkaline phosphatase, alanine aminotransferase, and aspartate aminotransferase abnormalities, as well as pro-inflammatory cytokines. Combination fenofibrate treatment also reduced 7α-hydroxy-4-cholesten-3-one (C4), the bile acid precursor, as well as total, primary, and conjugated bile acids. In addition, principal components analysis and heatmap analysis show that bile acid metabolites trended closer to that of healthy control subjects. These favorable effects of fenofibrate on bile acid metabolism may contribute to its beneficial clinical effects in patients with PBC and PSC experiencing a subtherapeutic response to Ursodiol monotherapy.

Characterization of the central metabolic pathways in Thermoanaerobacter sp. strain X514 via isotopomer-assisted metabolite analysis.

Thermoanaerobacter sp. strain X514 has great potential in biotechnology due to its capacity to ferment a range of C(5) and C(6) sugars to ethanol and other metabolites under thermophilic conditions. This study investigated the central metabolism of strain X514 via (13)C-labeled tracer experiments using either glucose or pyruvate as both carbon and energy sources. X514 grew on minimal medium and thus contains complete biosynthesis pathways for all macromolecule building blocks. Based on genome annotation and isotopic analysis of amino acids, three observations can be obtained about the central metabolic pathways in X514. First, the oxidative pentose phosphate pathway in X514 is not functional, and the tricarboxylic acid cycle is incomplete under fermentative growth conditions. Second, X514 contains (Re)-type citrate synthase activity, although no gene homologous to the recently characterized (Re)-type citrate synthase of Clostridium kluyveri was found. Third, the isoleucine in X514 is derived from acetyl coenzyme A and pyruvate via the citramalate pathway rather than being synthesized from threonine via threonine ammonia-lyase. The functionality of the citramalate synthase gene (cimA [Teth514_1204]) has been confirmed by enzymatic activity assays, while the presence of intracellular citramalate has been detected by mass spectrometry. This study demonstrates the merits of combining (13)C-assisted metabolite analysis, enzyme assays, and metabolite detection not only to examine genome sequence annotations but also to discover novel enzyme activities.

The genome sequence of the anaerobic, sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough.

Desulfovibrio vulgaris Hildenborough is a model organism for studying the energy metabolism of sulfate-reducing bacteria (SRB) and for understanding the economic impacts of SRB, including biocorrosion of metal infrastructure and bioremediation of toxic metal ions. The 3,570,858 base pair (bp) genome sequence reveals a network of novel c-type cytochromes, connecting multiple periplasmic hydrogenases and formate dehydrogenases, as a key feature of its energy metabolism. The relative arrangement of genes encoding enzymes for energy transduction, together with inferred cellular location of the enzymes, provides a basis for proposing an expansion to the 'hydrogen-cycling' model for increasing energy efficiency in this bacterium. Plasmid-encoded functions include modification of cell surface components, nitrogen fixation and a type-III protein secretion system. This genome sequence represents a substantial step toward the elucidation of pathways for reduction (and bioremediation) of pollutants such as uranium and chromium and offers a new starting point for defining this organism's complex anaerobic respiration.

NifH-Harboring Bacterial Community Composition across an Alaskan Permafrost Thaw Gradient.

Since nitrogen (N) is often limiting in permafrost soils, we investigated the N2-fixing genetic potential and the inferred taxa harboring those genes by sequencing nifH gene fragments in samples taken along a permafrost thaw gradient in an Alaskan boreal soil. Samples from minimally, moderately and extensively thawed sites were taken to a depth of 79 cm to encompass zones above and below the depth of the water table. NifH reads were translated with frameshift correction and 112,476 sequences were clustered at 5% amino acid dissimilarity resulting in 1,631 OTUs. Sample depth in relation to water table depth was correlated to differences in the NifH sequence classes with those most closely related to group I nifH-harboring Alpha- and Beta-Proteobacteria in higher abundance above water table depth while those related to group III nifH-harboring Delta Proteobacteria more abundant below. The most dominant below water table depth NifH sequences, comprising 1/3 of the total, were distantly related to Verrucomicrobia-Opitutaceae. Overall, these results suggest that permafrost thaw alters the class-level composition of N2-fixing communities in the thawed soil layers and that this distinction corresponds to the depth of the water table. These nifH data were also compared to nifH sequences obtained from a study at an Alaskan taiga site, and to those of other geographically distant, non-permafrost sites. The two Alaska sites were differentiated largely by changes in relative abundances of the same OTUs, whereas the non-Alaska sites were differentiated by the lack of many Alaskan OTUs, and the presence of unique halophilic, sulfate- and iron-reducing taxa in the Alaska sites.

Lateral Gene Transfer in a Heavy Metal-Contaminated-Groundwater Microbial Community.

UNLABELLED: Unraveling the drivers controlling the response and adaptation of biological communities to environmental change, especially anthropogenic activities, is a central but poorly understood issue in ecology and evolution. Comparative genomics studies suggest that lateral gene transfer (LGT) is a major force driving microbial genome evolution, but its role in the evolution of microbial communities remains elusive. To delineate the importance of LGT in mediating the response of a groundwater microbial community to heavy metal contamination, representative Rhodanobacter reference genomes were sequenced and compared to shotgun metagenome sequences. 16S rRNA gene-based amplicon sequence analysis indicated that Rhodanobacter populations were highly abundant in contaminated wells with low pHs and high levels of nitrate and heavy metals but remained rare in the uncontaminated wells. Sequence comparisons revealed that multiple geochemically important genes, including genes encoding Fe(2+)/Pb(2+) permeases, most denitrification enzymes, and cytochrome c553, were native to Rhodanobacter and not subjected to LGT. In contrast, the Rhodanobacter pangenome contained a recombinational hot spot in which numerous metal resistance genes were subjected to LGT and/or duplication. In particular, Co(2+)/Zn(2+)/Cd(2+) efflux and mercuric resistance operon genes appeared to be highly mobile within Rhodanobacter populations. Evidence of multiple duplications of a mercuric resistance operon common to most Rhodanobacter strains was also observed. Collectively, our analyses indicated the importance of LGT during the evolution of groundwater microbial communities in response to heavy metal contamination, and a conceptual model was developed to display such adaptive evolutionary processes for explaining the extreme dominance of Rhodanobacter populations in the contaminated groundwater microbiome. IMPORTANCE: Lateral gene transfer (LGT), along with positive selection and gene duplication, are the three main mechanisms that drive adaptive evolution of microbial genomes and communities, but their relative importance is unclear. Some recent studies suggested that LGT is a major adaptive mechanism for microbial populations in response to changing environments, and hence, it could also be critical in shaping microbial community structure. However, direct evidence of LGT and its rates in extant natural microbial communities in response to changing environments is still lacking. Our results presented in this study provide explicit evidence that LGT played a crucial role in driving the evolution of a groundwater microbial community in response to extreme heavy metal contamination. It appears that acquisition of genes critical for survival, growth, and reproduction via LGT is the most rapid and effective way to enable microorganisms and associated microbial communities to quickly adapt to abrupt harsh environmental stresses.

Comparative metagenomics reveals impact of contaminants on groundwater microbiomes.

To understand patterns of geochemical cycling in pristine versus contaminated groundwater ecosystems, pristine shallow groundwater (FW301) and contaminated groundwater (FW106) samples from the Oak Ridge Integrated Field Research Center (OR-IFRC) were sequenced and compared to each other to determine phylogenetic and metabolic difference between the communities. Proteobacteria (e.g., Burkholderia, Pseudomonas) are the most abundant lineages in the pristine community, though a significant proportion ( >55%) of the community is composed of poorly characterized low abundance (individually <1%) lineages. The phylogenetic diversity of the pristine community contributed to a broader diversity of metabolic networks than the contaminated community. In addition, the pristine community encodes redundant and mostly complete geochemical cycles distributed over multiple lineages and appears capable of a wide range of metabolic activities. In contrast, many geochemical cycles in the contaminated community appear truncated or minimized due to decreased biodiversity and dominance by Rhodanobacter populations capable of surviving the combination of stresses at the site. These results indicate that the pristine site contains more robust and encodes more functional redundancy than the stressed community, which contributes to more efficient nutrient cycling and adaptability than the stressed community.

Genomic insights into gene regulation of Desulfovibrio vulgaris Hildenborough.

Traditional laboratory studies of the sulfate-reducing bacteria have focused primarily on the biochemistry of the organisms. As genomic sequences of sulfate-reducing species have become available, insights have been gained into the metabolic and regulatory networks of these organisms. A computational analysis is reported of the transcriptional regulatory networks of Desulfovibrio vulgaris Hildenborough, the first mesophilic gram-negative sulfate-reducing bacterium for which a genome sequence is available. A set of conserved DNA motifs were derived from libraries of potential promoter regions of putative D. vulgaris regulons with the AlignACE program suite. Although one motif showed similarity to the Escherichia coli GlpR binding site, most of the motifs returned were apparently unique. A number of expected orthologs for regulatory proteins have not yet been recognized in D. vulgaris.