Metabolism and evolution of Haemophilus influenzae deduced from a whole-genome comparison with Escherichia coli

OrthoRefine: automated enhancement of prior ortholog identification via synteny

BACKGROUND

Identifying orthologs continues to be an early and imperative step in genome analysis but remains a challenging problem. While synteny (conservation of gene order) has previously been used independently and in combination with other methods to identify orthologs, applying synteny in ortholog identification has yet to be automated in a user-friendly manner. This desire for automation and ease-of-use led us to develop OrthoRefine, a standalone program that uses synteny to refine ortholog identification.

RESULTS

We developed OrthoRefine to improve the detection of orthologous genes by implementing a look-around window approach to detect synteny. We tested OrthoRefine in tandem with OrthoFinder, one of the most used software for identification of orthologs in recent years. We evaluated improvements provided by OrthoRefine in several bacterial and a eukaryotic dataset. OrthoRefine efficiently eliminates paralogs from orthologous groups detected by OrthoFinder. Using synteny increased specificity and functional ortholog identification; additionally, analysis of BLAST e-value, phylogenetics, and operon occurrence further supported using synteny for ortholog identification. A comparison of several window sizes suggested that smaller window sizes (eight genes) were generally the most suitable for identifying orthologs via synteny. However, larger windows (30 genes) performed better in datasets containing less closely related genomes. A typical run of OrthoRefine with ~ 10 bacterial genomes can be completed in a few minutes on a regular desktop PC.

CONCLUSION

OrthoRefine is a simple-to-use, standalone tool that automates the application of synteny to improve ortholog detection. OrthoRefine is particularly efficient in eliminating paralogs from orthologous groups delineated by standard methods.

The DmsABC S-oxide reductase is an essential component of a novel, hypochlorite-inducible system of extracellular stress defense in Haemophilus influenzae

Defenses against oxidative damage to cell components are essential for survival of bacterial pathogens during infection, and here we have uncovered that the DmsABC S-/N-oxide reductase is essential for virulence and in-host survival of the human-adapted pathogen, Haemophilus influenzae. In several different infection models, H. influenzae ΔdmsA strains showed reduced immunogenicity as well as lower levels of survival in contact with host cells. Expression of DmsABC was induced in the presence of hypochlorite and paraquat, closely linking this enzyme to defense against host-produced antimicrobials. In addition to methionine sulfoxide, DmsABC converted nicotinamide- and pyrimidine-N-oxide, precursors of NAD and pyrimidine for which H. influenzae is an auxotroph, at physiologically relevant concentrations, suggesting that these compounds could be natural substrates for DmsABC. Our data show that DmsABC forms part of a novel, periplasmic system for defense against host-induced S- and N-oxide stress that also comprises the functionally related MtsZ S-oxide reductase and the MsrAB peptide methionine sulfoxide reductase. All three enzymes are induced following exposure of the bacteria to hypochlorite. MsrAB is required for physical resistance to HOCl and protein repair. In contrast, DmsABC was required for intracellular colonization of host cells and, together with MtsZ, contributed to resistance to N-Chlorotaurine. Our work expands and redefines the physiological role of DmsABC and highlights the importance of different types of S-oxide reductases for bacterial virulence.

The environmentally-regulated interplay between local three-dimensional chromatin organisation and transcription of proVWX in E. coli

Nucleoid associated proteins (NAPs) maintain the architecture of bacterial chromosomes and regulate gene expression. Thus, their role as transcription factors may involve three-dimensional chromosome re-organisation. While this model is supported by in vitro studies, direct in vivo evidence is lacking. Here, we use RT-qPCR and 3C-qPCR to study the transcriptional and architectural profiles of the H-NS (histone-like nucleoid structuring protein)-regulated, osmoresponsive proVWX operon of Escherichia coli at different osmolarities and provide in vivo evidence for transcription regulation by NAP-mediated chromosome re-modelling in bacteria. By consolidating our in vivo investigations with earlier in vitro and in silico studies that provide mechanistic details of how H-NS re-models DNA in response to osmolarity, we report that activation of proVWX in response to a hyperosmotic shock involves the destabilization of H-NS-mediated bridges anchored between the proVWX downstream and upstream regulatory elements (DRE and URE), and between the DRE and ygaY that lies immediately downstream of proVWX. The re-establishment of these bridges upon adaptation to hyperosmolarity represses the operon. Our results also reveal additional structural features associated with changes in proVWX transcript levels such as the decompaction of local chromatin upstream of the operon, highlighting that further complexity underlies the regulation of this model operon. H-NS and H-NS-like proteins are wide-spread amongst bacteria, suggesting that chromosome re-modelling may be a typical feature of transcriptional control in bacteria.

Ribosomal Protein Cluster Organization in Asgard Archaea

It has been proposed that the superphylum of Asgard Archaea may represent a historical link between the Archaea and Eukarya. Following the discovery of the Archaea, it was soon appreciated that archaeal ribosomes were more similar to those of Eukarya rather than Bacteria. Coupled with other eukaryotic-like features, it has been suggested that the Asgard Archaea may be directly linked to eukaryotes. However, the genomes of Bacteria and non-Asgard Archaea generally organize ribosome-related genes into clusters that likely function as operons. In contrast, eukaryotes typically do not employ an operon strategy. To gain further insight into conservation of the r-protein genes, the genome order of conserved ribosomal protein (r-protein) coding genes was identified in 17 Asgard genomes (thirteen complete genomes and four genomes with less than 20 contigs) and compared with those found previously in non-Asgard archaeal and bacterial genomes. A universal core of two clusters of 14 and 4 cooccurring r-proteins, respectively, was identified in both the Asgard and non-Asgard Archaea. The equivalent genes in the E. coli version of the cluster are found in the S10 and spc operons. The large cluster of 14 r-protein genes (uS19-uL22-uS3-uL29-uS17 from the S10 operon and uL14-uL24-uL5-uS14-uS8-uL6-uL18-uS5-uL30-uL15 from the spc operon) occurs as a complete set in the genomes of thirteen Asgard genomes (five Lokiarchaeotes, three Heimdallarchaeotes, one Odinarchaeote, and four Thorarchaeotes). Four less conserved clusters with partial bacterial equivalents were found in the Asgard. These were the L30e (str operon in Bacteria) cluster, the L18e (alpha operon in Bacteria) cluster, the S24e-S27ae-rpoE1 cluster, and the L31e, L12..L1 cluster. Finally, a new cluster referred to as L7ae was identified. In many cases, r-protein gene clusters/operons are less conserved in their organization in the Asgard group than in other Archaea. If this is generally true for nonribosomal gene clusters, the results may have implications for the history of genome organization. In particular, there may have been an early transition to or from the operon approach to genome organization. Other nonribosomal cellular features may support different relationships. For this reason, it may be important to consider ribosome features separately.

In vivo gene expression profile of Haemophilus influenzae during human pneumonia

Haemophilus influenzae is a major cause of community-acquired pneumonia. While studied extensively in various laboratory models, less is known about the cell function while inside the human lung. We present the first analysis of the global gene expression of H. influenzae while the bacteria are in the lung during pneumonia (in vivo conditions) and contrast it with bacterial isolates that have been cultured under standard laboratory conditions (in vitro conditions). Patients with pneumonia were recruited from emergency departments and intensive care units during 2018-2020 (n = 102). Lower respiratory samples were collected for bacterial culture and RNA extraction. Patient samples with H. influenzae (n = 8) and colonies from bacterial cultures (n = 6) underwent RNA sequencing. The reads were then pseudo-aligned to core and pan genomes created from 15 reference strains. While bacteria cultured in vitro clustered tightly by principal component analysis of core genome (n = 1067) gene expression, bacteria in the patient samples had more diverse transcriptomic signatures and did not group with their lab-cultured counterparts. In total, 328 core genes were significantly differentially expressed between in vitro and in vivo conditions. The most highly upregulated genes in vivo included tbpA and fbpA, which are involved in the acquisition of iron from transferrin, and the stress response gene msrAB. The biosynthesis of nucleotides/purines and molybdopterin-scavenging processes were also significantly enriched in vivo. In contrast, major metabolic pathways and iron-sequestering genes were downregulated under this condition. In conclusion, extensive transcriptomic differences were found between bacteria while in the human lung and bacteria that were cultured in vitro. IMPORTANCE The human-specific pathogen Haemophilus influenzae is generally not well suited for studying in animal models, and most laboratory models are unlikely to approximate the diverse environments encountered by bacteria in the human airways accurately. Thus, we have examined the global gene expression of H. influenzae during pneumonia. Extensive differences in the global gene expression profiles were found in H. influenzae while in the human lung compared to bacteria that were grown in the laboratory. In contrast, the gene expression profiles of isolates collected from different patients were found to cluster together when grown under the same laboratory conditions. Interesting observations were made of how H. influenzae acquires and uses iron and molybdate, endures oxidative stress, and regulates central metabolism while in the lung. Our results indicate important processes during infection and can guide future research on genes and pathways that are relevant in the pathogenesis of H. influenzae pneumonia.

Multiplatform Metabolomics Characterization Reveals Novel Metabolites and Phospholipid Compositional Rules of Haemophilus influenzae Rd KW20

Haemophilus influenzae is a gram-negative bacterium of relevant clinical interest. H. influenzae Rd KW20 was the first organism to be sequenced and for which a genome-scale metabolic model (GEM) was developed. However, current H. influenzae GEMs are unable to capture several aspects of metabolome nature related to metabolite pools. To directly and comprehensively characterize the endometabolome of H. influenzae Rd KW20, we performed a multiplatform MS-based metabolomics approach combining LC-MS, GC-MS and CE-MS. We obtained direct evidence of 15-20% of the endometabolome present in current H. influenzae GEMs and showed that polar metabolite pools are interconnected through correlating metabolite islands. Notably, we obtained high-quality evidence of 18 metabolites not previously included in H. influenzae GEMs, including the antimicrobial metabolite cyclo(Leu-Pro). Additionally, we comprehensively characterized and evaluated the quantitative composition of the phospholipidome of H. influenzae, revealing that the fatty acyl chain composition is largely independent of the lipid class, as well as that the probability distribution of phospholipids is mostly related to the conditional probability distribution of individual acyl chains. This finding enabled us to provide a rationale for the observed phospholipid profiles and estimate the abundance of low-level species, permitting the expansion of the phospholipidome characterization through predictive probabilistic modelling.

Disturbances of purine and lipid metabolism in the microbiota-gut-brain axis in male adolescent nonhuman primates with depressive-like behaviors

INTRODUCTION

Major depressive disorder (MDD) in adolescents is a widespread and growing global public health concern with unique characteristics and pathophysiological mechanisms that are distinct from MDD in adults.

OBJECTIVE

The purpose of our work was to address this knowledge gap about the unique characteristics and pathophysiological mechanisms of adolescent depression from a microbiota-gut-brain (MGB) axis perspective.

METHOD

Ten healthy male cynomolgus macaques (Macaca fascicularis) were paired into five pairs based on age and body weight, and two cynomolgus macaques from each pair were randomly allocated to chronic unpredictable mild stress group, or unstressed control group. At endpoint, microbe composition from cecum, ascending colon, transverse colon, and descending colon were analyzed by metagenome sequencing, and the metabolite profiles of MGB axis including central (prefrontal cortex, hippocampus and amygdala) and peripheral (plasma, gut and feces of cecum, ascending colon, transverse colon and descending colon) samples were analyzed by metabolomic profiling. Then, we compare the gut microbiome and metabolic signatures in MGB axis between adolescent and adult depressed macaques.

RESULTS

The microbial composition and gut-brain metabolic signatures were widely divergent between adolescent and adult depressed macaques, though the phylum Firmicutes and lipid metabolism pathways were persistently altered in both populations. Purine and arginine biosynthesis metabolism were a specific hallmark of adolescent depressed macaques, while fatty acyl metabolism was specially altered in adult. These differential metabolic pathways in adolescent and adult depressed macaques were mainly mapped into the prefrontal cortex and hippocampus, respectively. Notably, the genus Clostridium and Haemophilus, characteristically disturbed in adolescent depressed macaques but not in adult, were also significantly associated with the majority of purine metabolites in MGB axis.

CONCLUSION

These findings provide a new framework describing divergent pathophysiological mechanisms between adolescent and adult depression, and may open new windows for more effective treatment strategies of adolescent depression.

The central role of arginine in Haemophilus influenzae survival in a polymicrobial environment with Streptococcus pneumoniae and Moraxella catarrhalis

Haemophilus influenzae, Streptococcus pneumoniae and Moraxella catarrhalis are bacterial species which frequently co-colonise the nasopharynx, but can also transit to the middle ear to cause otitis media. Chronic otitis media is often associated with a polymicrobial infection by these bacteria. However, despite being present in polymicrobial infections, the molecular interactions between these bacterial species remain poorly understood. We have previously reported competitive interactions driven by pH and growth phase between H. influenzae and S. pneumoniae. In this study, we have revealed competitive interactions between the three otopathogens, which resulted in reduction of H. influenzae viability in co-culture with S. pneumoniae and in triple-species culture. Transcriptomic analysis by mRNA sequencing identified a central role of arginine in mediating these interactions. Arginine supplementation was able to increase H. influenzae survival in a dual-species environment with S. pneumoniae, and in a triple-species environment. Arginine was used by H. influenzae for ATP production, and levels of ATP generated in dual- and triple-species co-culture at early stages of growth were significantly higher than the combined ATP levels of single-species cultures. These results indicate a central role for arginine-mediated ATP production by H. influenzae in the polymicrobial community.

Antibacterial Effects of Bacteriocin PLNC8 against Helicobacter pylori and Its Potential Mechanism of Action

Helicobacter pylori (H. pylori) is a bacterium that can cause a variety of gastric diseases. Most bacteriocins have gained popularity due to their non-toxic effects on cells and antibacterial effects against a wide range of pathogenic bacteria. In this study, the chemical synthesis of the bipeptide bacteriocin PLNC8 was used to investigate its possible action mechanism against H. pylori ZJC03 in vitro. Results showed that PLNC8 had significant anti-H. pylori ZJC03 potential, which resulted in a significant reduction in urease activity and a minimum inhibitory concentration (MIC) of 80 μM. PLNC8 inhibited the growth of H. pylori ZJC03, disrupting its structure as observed by scanning electron microscopy (SEM) and transmission electron microscopy (TEM). In addition, PLNC8 decreased the ATP level and hydrogen peroxide sensitivity of H. pylori ZJC03. In conclusion, PLNC8 disrupts the ability of H. pylori ZJC03 to alter the host environment, providing a new avenue for the prevention and control of H. pylori infection, providing a theoretical foundation for further elucidation of its regulatory mechanism.

NCMW: A Python Package to Analyze Metabolic Interactions in the Nasal Microbiome

The human upper respiratory tract is the reservoir of a diverse community of commensals and potential pathogens (pathobionts), including Streptococcus pneumoniae (pneumococcus), Haemophilus influenzae, Moraxella catarrhalis, and Staphylococcus aureus, which occasionally turn into pathogens causing infectious diseases, while the contribution of many nasal microorganisms to human health remains undiscovered. To better understand the composition of the nasal microbiome community, we create a workflow of the community model, which mimics the human nasal environment. To address this challenge, constraint-based reconstruction of biochemically accurate genome-scale metabolic models (GEMs) networks of microorganisms is mandatory. Our workflow applies constraint-based modeling (CBM), simulates the metabolism between species in a given microbiome, and facilitates generating novel hypotheses on microbial interactions. Utilizing this workflow, we hope to gain a better understanding of interactions from the metabolic modeling perspective. This article presents nasal community modeling workflow (NCMW)—a python package based on GEMs of species as a starting point for understanding the composition of the nasal microbiome community. The package is constructed as a step-by-step mathematical framework for metabolic modeling and analysis of the nasal microbial community. Using constraint-based models reduces the need for culturing species in vitro, a process that is not convenient in the environment of human noses.

Availability: NCMW is freely available on the Python Package Index (PIP) via pip install NCMW. The source code, documentation, and usage examples (Jupyter Notebook and example files) are available at https://github.com/manuelgloeckler/ncmw.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

Haemophilus influenzae Glucose Catabolism Leading to Production of the Immunometabolite Acetate Has a Key Contribution to the Host Airway-Pathogen Interplay

Chronic obstructive pulmonary disease (COPD) is characterized by abnormal inflammatory responses and impaired airway immunity, which provides an opportunistic platform for nontypeable Haemophilus influenzae (NTHi) infection. Clinical evidence supports that the COPD airways present increased concentrations of glucose, which may facilitate proliferation of pathogenic bacteria able to use glucose as a carbon source. NTHi metabolizes glucose through respiration-assisted fermentation, leading to the excretion of acetate, formate, and succinate. We hypothesized that such specialized glucose catabolism may be a pathoadaptive trait playing a pivotal role in the NTHi airway infection. To find out whether this is true, we engineered and characterized bacterial mutant strains impaired to produce acetate, formate, or succinate by inactivating the ackA, pflA, and frdA genes, respectively. While the inactivation of the pflA and frdA genes only had minimal physiological effects, the inactivation of the ackA gene affected acetate production and led to reduced bacterial growth, production of lactate under low oxygen tension, and bacterial attenuation in vivo. Moreover, bacterially produced acetate was able to stimulate the expression of inflammatory genes by cultured airway epithelial cells. These results back the notion that the COPD lung supports NTHi growth on glucose, enabling production of fermentative end products acting as immunometabolites at the site of infection. Thus, glucose catabolism may contribute not only to NTHi growth but also to bacterially driven airway inflammation. This information has important implications for developing nonantibiotic antimicrobials, given that airway glucose homeostasis modifying drugs could help prevent microbial infections associated with chronic lung disease.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Metabolic analyses reveal common adaptations in two invasive Haemophilus influenzae strains

Non-typeable Haemophilus influenzae (NTHi) is a major pathogen in upper and lower respiratory tract infections in humans, and is increasingly also associated with invasive disease. We have examined two unrelated NTHi invasive disease isolates, R2866 and C188, in order to identify metabolic and physiological properties that distinguish them from respiratory tract disease isolates such as Hi2019. While the general use of the Hi metabolic network was similar across all three strains, the two invasive isolates secreted increased amounts of succinate, which can have anti-inflammatory properties. In addition, they showed a common shift in their carbon source utilization patterns, with strongly enhanced metabolism of nucleoside substrates, glucose and sialic acid. The latter two are major compounds present in blood and cerebrospinal fluid (CSF). Interestingly, C188 and R2866 also shared a reduced ability to invade or survive intracellularly in 16HBE14 bronchial epithelial cells relative to Hi2019 (4-fold (4 h), 25-fold (24 h) reduction). Altered metabolic properties, such as the ones observed here, could arise from genomic adaptations that NTHi undergo during infection. Together these data indicate that shifts in substrate preferences in otherwise conserved metabolic pathways may underlie strain niche specificity and thus have the potential to alter the outcomes of host-NTHi interactions.

Operon Concatenation Is an Ancient Feature That Restricts the Potential to Rearrange Bacterial Chromosomes

The last common ancestor of the Gammaproteobacteria carried an important 40-kb chromosome section encoding 51 proteins of the transcriptional and translational machinery. These genes were organized into eight contiguous operons (rrnB-tufB-secE-rpoBC-str-S10-spc-alpha). Over 2 Gy of evolution, in different lineages, some of the operons became separated by multigene insertions. Surprisingly, in many Enterobacteriaceae, much of the ancient organization is conserved, indicating a strong selective force on the operons to remain colinear. Here, we show for one operon pair, tufB-secE in Salmonella, that an interruption of contiguity significantly reduces growth rate. Our data show that the tufB-secE operons are concatenated by an interoperon terminator–promoter overlap that plays a significant role regulating gene expression. Interrupting operon contiguity interferes with this regulation, reducing cellular fitness. Six operons of the ancestral chromosome section remain contiguous in Salmonella (tufB-secE-rpoBC and S10-spc-alpha) and, strikingly, each of these operon pairs is also connected by an interoperon terminator–promoter overlap. Accordingly, we propose that operon concatenation is an ancient feature that restricts the potential to rearrange bacterial chromosomes and can select for the maintenance of a colinear operon organization over billions of years.

Genomic characterization of Haemophilus influenzae: a focus on the capsule locus

Background

Haemophilus influenzae (Hi) can cause invasive diseases such as meningitis, pneumonia, or sepsis. Typeable Hi includes six serotypes (a through f), each expressing a unique capsular polysaccharide. The capsule, encoded by the genes within the capsule locus, is a major virulence factor of typeable Hi. Non-typeable (NTHi) does not express capsule and is associated with invasive and non-invasive diseases.

Methods

A total of 395 typeable and 293 NTHi isolates were characterized by whole genome sequencing (WGS). Phylogenetic analysis and multilocus sequence typing were used to characterize the overall genetic diversity. Pair-wise comparisons were used to evaluate the capsule loci. A WGS serotyping method was developed to predict the Hi serotype. WGS serotyping results were compared to slide agglutination (SAST) or real-time PCR (rt-PCR) serotyping.

Results

Isolates of each Hi serotype clustered into one or two subclades, with each subclade being associated with a distinct sequence type (ST). NTHi isolates were genetically diverse, with seven subclades and 125 STs being detected. Regions I and III of the capsule locus were conserved among the six serotypes (≥82% nucleotide identity). In contrast, genes in Region II were less conserved, with only six gene pairs from all serotypes showing ≥56% nucleotide identity. The WGS serotyping method was 99.9% concordant with SAST and 100% concordant with rt-PCR in determining the Hi serotype.

Conclusions

Genomic analysis revealed a higher degree of genetic diversity among NTHi compared to typeable Hi. The WGS serotyping method accurately predicted the Hi capsule type and can serve as an alternative method for Hi serotyping.

Assembly of Schizosaccharomyces cryophilus chromosomes and their comparative genomic analyses revealed principles of genome evolution of the haploid fission yeasts

The fission yeast clade, which has a distinct life history from other yeasts, can provide important clues about evolutionary changes. To reveal these changes the large S. cryophilus supercontigs were assembled into chromosomes using synteny relationships and the conserved pericentromeric, subtelomeric genes. Togetherness of the supercontigs was confirmed by PCR. Investigation of the gene order revealed localisation of the rDNA arrays, more than 300 new conserved orthologues and proved that S. cryophilus supercontigs were mosaics of collinear blocks. PFGE analysis showed that size of the S. cryophilus chromosomes differ from the S. pombe chromosomes. Comparative genomic analyses of the newly assembled chromosomes confirmed that the closest relative of S. cryophilus was S. octosporus not just in sequence similarity but also in a structural way, and revealed that preservation of the conserved regions did not arise from the lower number of chromosomal rearrangements. Translocations were more typical in the closely related species, while the number of inversions increased with the phylogenetic distances. Our data suggested that sites of the chromosomal rearrangements were not random and often associated with repetitive sequences, structural- and nucleotide evolution might correlate. Chromosomal rearrangements of the fission yeasts compared to other lineages were also discussed.

Experimental design and metabolic flux analysis tools to optimize industrially relevant Haemophilus influenzae type b growth medium

Haemophilus influenzae type b (Hib), a Gram-negative capsulated bacterium, is a causative agent of meningitis worldwide. The capsular polysaccharide, a high molecular mass polymer consisting of the repeated units of the polyribosyl-ribitol-phosphate, is considered the main virulence factor and it is used as an antigen to vaccines, conjugated to a carrier protein. The industrial production of the polysaccharide requires the cultivation of Hib in rich medium, which impacts process costs and product recovery. In this study, a central composite rotational experimental design strategy was used to access the influence of key components of culture medium (soy peptone, yeast extract and glucose) on biomass formation and polysaccharide production in shake-flasks. The optimized medium formulation, containing half of the usual yeast extract and soytone concentrations, was further validated in batch bioreactor cultivations. High polysaccharide production (∼500 mg/L) was obtained in a cheaper and more competitive production process for use in Hib vaccine production. In addition, simulations of a metabolic model describing Hib central metabolism were used to assess the role of key amino acids on growth. A chemically defined medium supplemented only with amino acids from α-ketoglutarate and oxaloacetate families as well as phenylalanine was suggested as a promising alternative for reduced acetate accumulation and enhanced polysaccharide production in Hib cultures. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:1508-1519, 2017.

DEF: an automated dead-end filling approach based on quasi-endosymbiosis

Motivation

Gap filling for the reconstruction of metabolic networks is to restore the connectivity of metabolites via finding high-confidence reactions that could be missed in target organism. Current methods for gap filling either fall into the network topology or have limited capability in finding missing reactions that are indirectly related to dead-end metabolites but of biological importance to the target model.

Results

We present an automated dead-end filling (DEF) approach, which is derived from the wisdom of endosymbiosis theory, to fill gaps by finding the most efficient dead-end utilization paths in a constructed quasi-endosymbiosis model. The recalls of reactions and dead ends of DEF reach around 73% and 86%, respectively. This method is capable of finding indirectly dead-end-related reactions with biological importance for the target organism and is applicable to any given metabolic model. In the E. coli iJR904 model, for instance, about 42% of the dead-end metabolites were fixed by our proposed method.

Availabilty and Implementaion

DEF is publicly available at http://bis.zju.edu.cn/DEF/.

Contact

mchen@zju.edu.cn

Supplimentary Information

Supplementary data are available at Bioinformatics online.

Metabolic versatility in Haemophilus influenzae: a metabolomic and genomic analysis

Haemophilus influenzae is a host adapted human pathogen known to contribute to a variety of acute and chronic diseases of the upper and lower respiratory tract as well as the middle ear. At the sites of infection as well as during growth as a commensal the environmental conditions encountered by H. influenzae will vary significantly, especially in terms of oxygen availability, however, the mechanisms by which the bacteria can adapt their metabolism to cope with such changes have not been studied in detail. Using targeted metabolomics the spectrum of metabolites produced during growth of H. influenzae on glucose in RPMI-based medium was found to change from acetate as the main product during aerobic growth to formate as the major product during anaerobic growth. This change in end-product is likely caused by a switch in the major route of pyruvate degradation. Neither lactate nor succinate or fumarate were major products of H. influenzae growth under any condition studied. Gene expression studies and enzyme activity data revealed that despite an identical genetic makeup and very similar metabolite production profiles, H. influenzae strain Rd appeared to favor glucose degradation via the pentose phosphate pathway, while strain 2019, a clinical isolate, showed higher expression of enzymes involved in glycolysis. Components of the respiratory chain were most highly expressed during microaerophilic and anaerobic growth in both strains, but again clear differences existed in the expression of genes associated e.g., with NADH oxidation, nitrate and nitrite reduction in the two strains studied. Together our results indicate that H. influenzae uses a specialized type of metabolism that could be termed “respiration assisted fermentation” where the respiratory chain likely serves to alleviate redox imbalances caused by incomplete glucose oxidation, and at the same time provides a means of converting a variety of compounds including nitrite and nitrate that arise as part of the host defence mechanisms.

Gene identification in prokaryotic genomes, phages, metagenomes, and EST sequences with GeneMarkS suite

This unit describes how to use several gene-finding programs from the GeneMark line developed for finding protein-coding ORFs in genomic DNA of prokaryotic species, in genomic DNA of eukaryotic species with intronless genes, in genomes of viruses and phages, and in prokaryotic metagenomic sequences, as well as in EST sequences with spliced-out introns. These bioinformatics tools were demonstrated to have state-of-the-art accuracy, and have been frequently used for gene annotation in novel nucleotide sequences. An additional advantage of these sequence-analysis tools is that the problem of algorithm parameterization is solved automatically, with parameters estimated by iterative self-training (unsupervised training).

Comparative genomic analysis reveals distinct genotypic features of the emerging pathogen Haemophilus influenzae type f

BACKGROUND

The incidence of invasive disease caused by encapsulated Haemophilus influenzae type f (Hif) has increased in the post-H. influenzae type b (Hib) vaccine era. We previously annotated the first complete Hif genome from a clinical isolate (KR494) that caused septic shock and necrotizing myositis. Here, the full genome of Hif KR494 was compared to sequenced reference strains Hib 10810, capsule type d (Hid) Rd Kw20, and finally nontypeable H. influenzae 3655. The goal was to identify possible genomic characteristics that may shed light upon the pathogenesis of Hif.

RESULTS

The Hif KR494 genome exhibited large regions of synteny with other H. influenzae, but also distinct genome rearrangements. A predicted Hif core genome of 1390 genes was shared with the reference strains, and 6 unique genomic regions comprising half of the 191 unique coding sequences were revealed. The majority of these regions were inserted genetic fragments, most likely derived from the closely-related Haemophilus spp. including H. aegyptius, H. haemolyticus and H. parainfluenzae. Importantly, the KR494 genome possessed several putative virulence genes that were distinct from non-type f strains. These included the sap2 operon, aef3 fimbriae, and genes for kanamycin nucleotidyltranserase, iron-utilization proteins, and putative YadA-like trimeric autotransporters that may increase the bacterial virulence. Furthermore, Hif KR494 lacked a hisABCDEFGH operon for de novo histidine biosynthesis, hmg locus for lipooligosaccharide biosynthesis and biofilm formation, the Haemophilus antibiotic resistance island and a Haemophilus secondary molybdate transport system. We confirmed the histidine auxotrophy and kanamycin resistance in Hif by functional experiments. Moreover, the pattern of unique or missing genes of Hif KR494 was similar in 20 Hif clinical isolates obtained from different years and geographical areas. A cross-species comparison revealed that the Hif genome shared more characteristics with H. aegyptius than Hid and NTHi.

CONCLUSIONS

The genomic comparative analyses facilitated identification of genotypic characteristics that may be related to the specific virulence of Hif. In relation to non-type f H. influenzae strains, the Hif genome contains differences in components involved in metabolism and survival that may contribute to its invasiveness.

Genome-derived vaccines

Vaccine research entered a new era when the complete genome of a pathogenic bacterium was published in 1995. Since then, more than 97 bacterial pathogens have been sequenced and at least 110 additional projects are now in progress. Genome sequencing has also dramatically accelerated: high-throughput facilities can draft the sequence of an entire microbe (two to four megabases) in 1 to 2 days. Vaccine developers are using microarrays, immunoinformatics, proteomics and high-throughput immunology assays to reduce the truly unmanageable volume of information available in genome databases to a manageable size. Vaccines composed by novel antigens discovered from genome mining are already in clinical trials. Within 5 years we can expect to see a novel class of vaccines composed by genome-predicted, assembled and engineered T- and Bcell epitopes. This article addresses the convergence of three forces--microbial genome sequencing, computational immunology and new vaccine technologies--that are shifting genome mining for vaccines onto the forefront of immunology research.

"Affect of anaerobiosis on the antibiotic susceptibility of H. influenzae"

Background

Haemophilus influenzae is a human-restricted facultative anaerobe which resides mostly in the oropharynx. The majority of isolates recovered from the throat are unencapsulated commensals (NTHi), but depending on host susceptibility they cause bronchitis, otitis media and on occasion bacteremia and meningitis. Because of the variable oxygen availability in the various niche permitting bacterium replication, the organism must thrive in well oxygenated surfaces, such as pharyngeal epithelium to anoxic environments like the bottom of a Biofilm and in airway mucus. Other reports indicate that H. influenzae use aerobic respiration, anaerobic respiration and fermentation to generate ATP. To gain insight in to the activity of several classes of antibiotics against five well-characterized unencapsulated H. influenzae in room air, in 5% CO₂ and under strict anaerobiosis. We also tested for the role of oxidative killing by all cidal antibiotics.

Results

In comparison to room air, testing in 5% CO₂ had minimal effects on the susceptibility to aminoglycosides, cephalosporins, tetracycline and chloramphenicol: the MIC of rifampin and ciprofloxacin increased eight fold with certain strains in 5% CO₂. All antibiotics, except trimethoprim were cidal under both growth conditions. Aminoglycosides remained bactericidal in a strict anaerobic environment, while a reliable MBC was obtained with trimethoprim only under anaerobic conditions. Kinetic analysis of the cidal action of spectinomycin and tetracycline indicated slower killing anaerobically. An oxidative mechanism for aerobic killing could not be demonstrated.

Conclusions

We conclude that β-lactams, cephalosporins, macrolides, tetracycline’s, aminoglycosides, chloramphenicol, rifampin and ciprofloxacin are bactericidal against five well-characterizes H. influenzae in an aerobic and anaerobic environment. The activity of trimethoprim was increased in anaerobic conditions.

Gene identification in prokaryotic genomes, phages, metagenomes, and EST sequences with GeneMarkS suite

This unit describes how to use several gene-finding programs from the GeneMark line developed for finding protein-coding ORFs in genomic DNA of prokaryotic species, in genomic DNA of eukaryotic species with intronless genes, in genomes of viruses and phages, and in prokaryotic metagenomic sequences, as well as in EST sequences with spliced-out introns. These bioinformatics tools were demonstrated to have state-of-the-art accuracy and have been frequently used for gene annotation in novel nucleotide sequences. An additional advantage of these sequence-analysis tools is that the problem of algorithm parameterization is solved automatically, with parameters estimated by iterative self-training (unsupervised training).

A novel esterase gene cloned from a metagenomic library from neritic sediments of the South China Sea

BACKGROUND

Marine microbes are a large and diverse group, which are exposed to a wide variety of pressure, temperature, salinity, nutrient availability and other environmental conditions. They provide a huge potential source of novel enzymes with unique properties that may be useful in industry and biotechnology. To explore the lipolytic genetic resources in the South China Sea, 23 sediment samples were collected in the depth < 100 m marine areas.

RESULTS

A metagenomic library of South China Sea sediments assemblage in plasmid vector containing about 194 Mb of community DNA was prepared. Screening of a part of the unamplified library resulted in isolation of 15 unique lipolytic clones with the ability to hydrolyze tributyrin. A positive recombinant clone (pNLE1), containing a novel esterase (Est_p1), was successfully expressed in E. coli and purified. In a series of assays, Est_p1 displayed maximal activity at pH 8.57, 40°C, with ρ-Nitrophenyl butyrate (C4) as substrate. Compared to other metagenomic esterases, Est_p1 played a notable role in specificity for substrate C4 (kcat/Km value 11,500 S-1m M-1) and showed no inhibited by phenylmethylsulfonyl fluoride, suggested that the substrate binding pocket was suitable for substrate C4 and the serine active-site residue was buried at the bottom of substrate binding pocket which sheltered by a lid structure.

CONCLUSIONS

Esterase, which specificity towards short chain fatty acids, especially butanoic acid, is commercially available as potent flavoring tools. According the outstanding activity and specificity for substrate C4, Est_p1 has potential application in flavor industries requiring hydrolysis of short chain esters.

A novel nickel responsive MerR-like regulator, NimR, from Haemophilus influenzae

We have identified a novel regulator from the MerR family of transcription factors in the bacterial pathogen Haemophilus influenzae (HI1623; nickel-associated merR-like Regulator--NimR). NimR regulates the expression of a Ni(2+) uptake transporter (NikKLMQO). The promoters for nimR and the nik operon are divergent and overlapping and NimR binds at a site between the promoter elements for nikKLMQO. Expression of this operon requires NimR and depends on Ni(2+). Growth rates of the H. influenzae nimR and nikQ mutants were reduced in chemically defined media compared to the wild type and the mutants were unable to grow in the presence of EDTA. The mutant strains were less tolerant of acidic pH and the wild type Rd KW20 could not tolerate low pH in the presence of fluoramide, a urease specific inhibitor, confirming that both nickel transport and urea hydrolysis are a central process in pH control. H. influenzae nimR and nikQ strains were deficient in urease activity, but this could be specifically restored by the addition of excess Ni(2+). NimR did not directly regulate the expression of urease genes but the activity of urease requires both nimR and nikQ. Purified NimR is a dimer that binds 1 Ni(2+)ion. NimR is the first example of a Ni-dependent regulator from the MerR family and targeting a metal ion uptake system; it is distinct from NikR the Ni-responsive regulators of the ribbon-helix-helix family.

Evaluating ortholog prediction algorithms in a yeast model clade

BACKGROUND

Accurate identification of orthologs is crucial for evolutionary studies and for functional annotation. Several algorithms have been developed for ortholog delineation, but so far, manually curated genome-scale biological databases of orthologous genes for algorithm evaluation have been lacking. We evaluated four popular ortholog prediction algorithms (MultiParanoid; and OrthoMCL; RBH: Reciprocal Best Hit; RSD: Reciprocal Smallest Distance; the last two extended into clustering algorithms cRBH and cRSD, respectively, so that they can predict orthologs across multiple taxa) against a set of 2,723 groups of high-quality curated orthologs from 6 Saccharomycete yeasts in the Yeast Gene Order Browser.

RESULTS

Examination of sensitivity [TP/(TP+FN)], specificity [TN/(TN+FP)], and accuracy [(TP+TN)/(TP+TN+FP+FN)] across a broad parameter range showed that cRBH was the most accurate and specific algorithm, whereas OrthoMCL was the most sensitive. Evaluation of the algorithms across a varying number of species showed that cRBH had the highest accuracy and lowest false discovery rate [FP/(FP+TP)], followed by cRSD. Of the six species in our set, three descended from an ancestor that underwent whole genome duplication. Subsequent differential duplicate loss events in the three descendants resulted in distinct classes of gene loss patterns, including cases where the genes retained in the three descendants are paralogs, constituting 'traps' for ortholog prediction algorithms. We found that the false discovery rate of all algorithms dramatically increased in these traps.

CONCLUSIONS

These results suggest that simple algorithms, like cRBH, may be better ortholog predictors than more complex ones (e.g., OrthoMCL and MultiParanoid) for evolutionary and functional genomics studies where the objective is the accurate inference of single-copy orthologs (e.g., molecular phylogenetics), but that all algorithms fail to accurately predict orthologs when paralogy is rampant.

Systems-level approaches for identifying and analyzing genetic interaction networks in Escherichia coli and extensions to other prokaryotes

Molecular interactions define the functional organization of the cell. Epistatic (genetic, or gene-gene) interactions, one of the most informative and commonly encountered forms of functional relationships, are increasingly being used to map process architecture in model eukaryotic organisms. In particular, 'systems-level' screens in yeast and worm aimed at elucidating genetic interaction networks have led to the generation of models describing the global modular organization of gene products and protein complexes within a cell. However, comparable data for prokaryotic organisms have not been available. Given its ease of growth and genetic manipulation, the Gram-negative bacterium Escherichia coli appears to be an ideal model system for performing comprehensive genome-scale examinations of genetic redundancy in bacteria. In this review, we highlight emerging experimental and computational techniques that have been developed recently to examine functional relationships and redundancy in E. coli at a systems-level, and their potential application to prokaryotes in general. Additionally, we have scanned PubMed abstracts and full-text published articles to manually curate a list of approximately 200 previously reported synthetic sick or lethal genetic interactions in E. coli derived from small-scale experimental studies.

Comparative transcriptomics in the Triticeae

Background

Barley and particularly wheat are two grass species of immense agricultural importance. In spite of polyploidization events within the latter, studies have shown that genotypically and phenotypically these species are very closely related and, indeed, fertile hybrids can be created by interbreeding. The advent of two genome-scale Affymetrix GeneChips now allows studies of the comparison of their transcriptomes.

Results

We have used the Wheat GeneChip to create a "gene expression atlas" for the wheat transcriptome (cv. Chinese Spring). For this, we chose mRNA from a range of tissues and developmental stages closely mirroring a comparable study carried out for barley (cv. Morex) using the Barley1 GeneChip. This, together with large-scale clustering of the probesets from the two GeneChips into "homologous groups", has allowed us to perform a genomic-scale comparative study of expression patterns in these two species. We explore the influence of the polyploidy of wheat on the results obtained with the Wheat GeneChip and quantify the correlation between conservation in gene sequence and gene expression in wheat and barley. In addition, we show how the conservation of expression patterns can be used to elucidate, probeset by probeset, the reliability of the Wheat GeneChip.

Conclusion

While there are many differences in expression on the level of individual genes and tissues, we demonstrate that the wheat and barley transcriptomes appear highly correlated. This finding is significant not only because given small evolutionary distance between the two species it is widely expected, but also because it demonstrates that it is possible to use the two GeneChips for comparative studies. This is the case even though their probeset composition reflects rather different design principles as well as, of course, the present incomplete knowledge of the gene content of the two species. We also show that, in general, the Wheat GeneChip is not able to distinguish contributions from individual homoeologs. Furthermore, the comparison between the two species leads us to conclude that the conservation of both gene sequence as well as gene expression is positively correlated with absolute expression levels, presumably reflecting increased selection pressure on genes coding for proteins present at high levels. In addition, the results indicate the presence of a correlation between sequence and expression conservation within the Triticeae.

Uncovering metabolic pathways relevant to phenotypic traits of microbial genomes

Identifying the biochemical basis of microbial phenotypes is a main objective of comparative genomics. Here we present a novel method using multivariate machine learning techniques for comparing automatically derived metabolic reconstructions of sequenced genomes on a large scale. Applying our method to 266 genomes directly led to testable hypotheses such as the link between the potential of microorganisms to cause periodontal disease and their ability to degrade histidine, a link also supported by clinical studies.

Trends in prokaryotic evolution revealed by comparison of closely related bacterial and archaeal genomes

In order to explore microevolutionary trends in bacteria and archaea, we constructed a data set of 41 alignable tight genome clusters (ATGCs). We show that the ratio of the medians of nonsynonymous to synonymous substitution rates (dN/dS) that is used as a measure of the purifying selection pressure on protein sequences is a stable characteristic of the ATGCs. In agreement with previous findings, parasitic bacteria, notwithstanding the sometimes dramatic genome shrinkage caused by gene loss, are typically subjected to relatively weak purifying selection, presumably owing to relatively small effective population sizes and frequent bottlenecks. However, no evidence of genome streamlining caused by strong selective pressure was found in any of the ATGCs. On the contrary, a significant positive correlation between the genome size, as well as gene size, and selective pressure was observed, although a variety of free-living prokaryotes with very close selective pressures span nearly the entire range of genome sizes. In addition, we examined the connections between the sequence evolution rate and other genomic features. Although gene order changes much faster than protein sequences during the evolution of prokaryotes, a strong positive correlation was observed between the "rearrangement distance" and the amino acid distance, suggesting that at least some of the events leading to genome rearrangement are subjected to the same type of selective constraints as the evolution of amino acid sequences.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

Prokaryotic gene prediction using GeneMark and GeneMark.hmm

In this unit, the GeneMark and GeneMark.hmm programs are presented as two different methods for the in silico prediction of genes in prokaryotes. GeneMark can be used for whole genome analysis as well as for the local analysis of a particular gene and its surrounding regions. GeneMark.hmm makes use of Hidden Markov models to find the transition points (boundaries) between protein coding states and noncoding states and can be efficiently used for larger genome sequences. These methods can be used in conjunction with each other for a higher sensitivity of gene detection.

Whole genome-based assessment of the taxonomic position of the arthropod pathogenic bacterium Rickettsiella grylli

Rickettsiella grylli is an intracellular bacterial pathogen of aquatic and terrestrial arthropods. Previous determination of its 16S rRNA-encoding sequence has led to the taxonomic classification of the genus Rickettsiella in the class Gammaproteobacteria, order Legionellales, family Coxiellaceae, i.e. in close vicinity to vertebrate pathogenic bacteria of the genera Coxiella and Legionella. Here we use the additional information available from the recently published first whole genome sequence from this genus to evaluate critically the taxonomic classification of R. grylli beyond the 16S rRNA gene level. Using phylogenetic reconstruction, together with significance testing on a data basis defined by a core set of 211 previously identified families of protein-encoding genes, together with a reanalysis of 16S rRNA gene data, the present study firmly corroborates the assignment of this species to both the class Gammaproteobacteria and the order Legionellales. However, the results obtained from concatenated and single protein, single protein-encoding gene, and 16S rRNA gene data demonstrate a similar phylogenetic distance of R. grylli to both the Coxiellaceae and the Legionellaceae and are, therefore, inconsistent with its current family-level classification. Consequently, a respective reorganization of the order Legionellales is proposed.

Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea

Background

An evolutionary classification of genes from sequenced genomes that distinguishes between orthologs and paralogs is indispensable for genome annotation and evolutionary reconstruction. Shortly after multiple genome sequences of bacteria, archaea, and unicellular eukaryotes became available, an attempt on such a classification was implemented in Clusters of Orthologous Groups of proteins (COGs). Rapid accumulation of genome sequences creates opportunities for refining COGs but also represents a challenge because of error amplification. One of the practical strategies involves construction of refined COGs for phylogenetically compact subsets of genomes.

Results

New Archaeal Clusters of Orthologous Genes (arCOGs) were constructed for 41 archaeal genomes (13 Crenarchaeota, 27 Euryarchaeota and one Nanoarchaeon) using an improved procedure that employs a similarity tree between smaller, group-specific clusters, semi-automatically partitions orthology domains in multidomain proteins, and uses profile searches for identification of remote orthologs. The annotation of arCOGs is a consensus between three assignments based on the COGs, the CDD database, and the annotations of homologs in the NR database. The 7538 arCOGs, on average, cover ~88% of the genes in a genome compared to a ~76% coverage in COGs. The finer granularity of ortholog identification in the arCOGs is apparent from the fact that 4538 arCOGs correspond to 2362 COGs; ~40% of the arCOGs are new. The archaeal gene core (protein-coding genes found in all 41 genome) consists of 166 arCOGs. The arCOGs were used to reconstruct gene loss and gene gain events during archaeal evolution and gene sets of ancestral forms. The Last Archaeal Common Ancestor (LACA) is conservatively estimated to possess 996 genes compared to 1245 and 1335 genes for the last common ancestors of Crenarchaeota and Euryarchaeota, respectively. It is inferred that LACA was a chemoautotrophic hyperthermophile that, in addition to the core archaeal functions, encoded more idiosyncratic systems, e.g., the CASS systems of antivirus defense and some toxin-antitoxin systems.

Conclusion

The arCOGs provide a convenient, flexible framework for functional annotation of archaeal genomes, comparative genomics and evolutionary reconstructions. Genomic reconstructions suggest that the last common ancestor of archaea might have been (nearly) as advanced as the modern archaeal hyperthermophiles. ArCOGs and related information are available at: .

Reviewers

This article was reviewed by Peer Bork, Patrick Forterre, and Purificacion Lopez-Garcia.

Metabolic analysis of Moraxella catarrhalis and the effect of selected in vitro growth conditions on global gene expression

The nucleotide sequence from the genome of Moraxella catarrhalis ATCC 43617 was annotated and used both to assess the metabolic capabilities and limitations of this bacterium and to design probes for a DNA microarray. An absence of gene products for utilization of exogenous carbohydrates was noteworthy and could be correlated with published phenotypic data. Gene products necessary for aerobic energy generation were present, as were a few gene products generally ascribed to anaerobic systems. Enzymes for synthesis of all amino acids except proline and arginine were present. M. catarrhalis DNA microarrays containing 70-mer oligonucleotide probes were designed from the genome-derived nucleotide sequence data. Analysis of total RNA extracted from M. catarrhalis ATCC 43617 cells grown under iron-replete and iron-restricted conditions was used to establish the utility of these DNA microarrays. These DNA microarrays were then used to analyze total RNA from M. catarrhalis cells grown in a continuous-flow biofilm system and in the planktonic state. The genes whose expression was most dramatically increased by growth in the biofilm state included those encoding a nitrate reductase, a nitrite reductase, and a nitric oxide reductase. Real-time reverse transcriptase PCR analysis was used to validate these DNA microarray results. These results indicate that growth of M. catarrhalis in a biofilm results in increased expression of gene products which can function not only in energy generation but also in resisting certain elements of the innate immune response.

Optimization based automated curation of metabolic reconstructions

Background

Currently, there exists tens of different microbial and eukaryotic metabolic reconstructions (e.g., Escherichia coli, Saccharomyces cerevisiae, Bacillus subtilis) with many more under development. All of these reconstructions are inherently incomplete with some functionalities missing due to the lack of experimental and/or homology information. A key challenge in the automated generation of genome-scale reconstructions is the elucidation of these gaps and the subsequent generation of hypotheses to bridge them.

Results

In this work, an optimization based procedure is proposed to identify and eliminate network gaps in these reconstructions. First we identify the metabolites in the metabolic network reconstruction which cannot be produced under any uptake conditions and subsequently we identify the reactions from a customized multi-organism database that restores the connectivity of these metabolites to the parent network using four mechanisms. This connectivity restoration is hypothesized to take place through four mechanisms: a) reversing the directionality of one or more reactions in the existing model, b) adding reaction from another organism to provide functionality absent in the existing model, c) adding external transport mechanisms to allow for importation of metabolites in the existing model and d) restore flow by adding intracellular transport reactions in multi-compartment models. We demonstrate this procedure for the genome- scale reconstruction of Escherichia coli and also Saccharomyces cerevisiae wherein compartmentalization of intra-cellular reactions results in a more complex topology of the metabolic network. We determine that about 10% of metabolites in E. coli and 30% of metabolites in S. cerevisiae cannot carry any flux. Interestingly, the dominant flow restoration mechanism is directionality reversals of existing reactions in the respective models.

Conclusion

We have proposed systematic methods to identify and fill gaps in genome-scale metabolic reconstructions. The identified gaps can be filled both by making modifications in the existing model and by adding missing reactions by reconciling multi-organism databases of reactions with existing genome-scale models. Computational results provide a list of hypotheses to be queried further and tested experimentally.

Multiple consecutive lavage samplings reveal greater burden of disease and provide direct access to the nontypeable Haemophilus influenzae biofilm in experimental otitis media

The typically recovered quantity of nontypeable Haemophilus influenzae (NTHi) bacteria in an ex vivo middle ear (ME) aspirate from the chinchilla model of experimental otitis media is insufficient for direct analysis of gene expression by microarray or of lipopolysaccharide glycoforms by mass spectrometry. This prompted us to investigate a strategy of multiple consecutive lavage samplings to increase ex vivo bacterial recovery. As multiple consecutive lavage samples significantly increased the total number of bacterial CFU collected during nasopharyngeal colonization or ME infection, this led us to evaluate whether bacteria sequentially acquired from consecutive lavages were similar. Comparative observation of complete ex vivo sample series by microscopy initially revealed ME inflammatory fluid consisting solely of planktonic-phase NTHi. In contrast, subsequent lavage samplings of the same infected ear revealed the existence of bacteria in two additional growth states, filamentous and biofilm encased. Gene expression analysis of such ex vivo samples was in accord with different bacterial growth phases in sequential lavage specimens. The existence of morphologically distinct NTHi subpopulations with varying levels of gene expression indicates that the pooling of specimens requires caution until methods for their separation are developed. This study based on multiple consecutive lavages is consistent with prior reports that NTHi forms a biofilm in vivo, describes the means to directly acquire ex vivo biofilm samples without sacrificing the animal, and has broad applicability for a study of mucosal infections. Moreover, this approach revealed that the actual burden of bacteria in experimental otitis media is significantly greater than was previously reported. Such findings may have direct implications for antibiotic treatment and vaccine development against NTHi.

Histidine auxotrophy in commensal and disease-causing nontypeable Haemophilus influenzae

Histidine biosynthesis is one of the best studied metabolic pathways in bacteria. Although this pathway is thought to be highly conserved within and between bacterial species, a previous study identified a genetic region within the histidine operon (his) of nontypeable strains of Haemophilus influenzae (NTHI) that was more prevalent among otitis media strains than among throat commensal NTHI strains. In the present study, we further characterized this region and showed that genes in the complete his operon (hisG, -D, -C, -NB, -H, -A, -F, and -IE) are >99% conserved among four fully sequenced NTHI strains, are present in the same location in these four genomes, and are situated in the same gene order. Using PCR and dot blot hybridization, we determined that the his operon was significantly more prevalent in otitis media NTHI strains (106/121; 87.7%) than in throat strains (74/137; 54%) (prevalence ratio, 1.62; P<0.0001), suggesting a possible role in middle ear survival and/or acute otitis media. NTHI strains lacking the his operon showed attenuated growth in histidine-restricted media, confirming them as his-negative auxotrophs. Our results suggest that the ability to make histidine is an important factor in bacterial growth and survival in the middle ear, where nutrients such as histidine may be found in limited amounts. Those isolates lacking the histidine pathway were still able to survive well in the throat, which suggests that histidine is readily available in the throat environment.

Rebuilding microbial genomes

Engineered microbes are of great potential utility in biotechnology and basic research. In principle, a cell can be built from scratch by assembling small molecule sets with auto-catalytic properties. Alternatively, DNA can be isolated or directly synthesized and molded into a synthetic genome using existing genomic blueprints and molecular biology tools. Activating such a synthetic genome will yield a synthetic cell. Here we examine obstacles associated with this latter approach using a model system whereby a donor genome from H. influenzae is fragmented, and the pieces are then modified and reassembled stepwise in an E. coli host cell. There are obstacles associated with this strategy related to DNA transfer, DNA replication, cross-talk in gene regulation and compatibility of gene products between donor and host. Encouragingly, analysis of gene expression indicates widespread transcription of H. influenzae genes in E. coli, and analysis of gap locations in H. influenzae and other microbial genome assemblies reveals few genes routinely incompatible with E. coli. In conclusion, rebuilding and booting a genome remains a feasible and pragmatic approach to creating a synthetic microbial cell.

Complete genome sequence of Haemophilus somnus (Histophilus somni) strain 129Pt and comparison to Haemophilus ducreyi 35000HP and Haemophilus influenzae Rd

Haemophilus somnus can be either a commensal of bovine mucosal surfaces or an opportunistic pathogen. Pathogenic strains of H. somnus are a significant cause of systemic disease in cattle. We report the genome sequence of H. somnus 129Pt, a nonpathogenic commensal preputial isolate, and the results of a genome-wide comparative analysis of H. somnus 129Pt, Haemophilus influenzae Rd, and Haemophilus ducreyi 35000HP. We found unique genes in H. somnus 129Pt involved in lipooligosaccharide biosynthesis, carbohydrate uptake and metabolism, cation transport, amino acid metabolism, ubiquinone and menaquinone biosynthesis, cell surface adhesion, biosynthesis of cofactors, energy metabolism, and electron transport. There were also many genes in common among the three organisms. Our comparative analyses of H. somnus 129Pt, H. influenzae Rd, and H. ducreyi 35000HP revealed similarities and differences in the numbers and compositions of genes involved in metabolism, host colonization, and persistence. These results lay a foundation for research on the host specificities and niche preferences of these organisms. Future comparisons between H. somnus 129Pt and virulent strains will aid in the development of protective strategies and vaccines to protect cattle against H. somnus disease.

Identifying metabolic enzymes with multiple types of association evidence

Background

Existing large-scale metabolic models of sequenced organisms commonly include enzymatic functions which can not be attributed to any gene in that organism. Existing computational strategies for identifying such missing genes rely primarily on sequence homology to known enzyme-encoding genes.

Results

We present a novel method for identifying genes encoding for a specific metabolic function based on a local structure of metabolic network and multiple types of functional association evidence, including clustering of genes on the chromosome, similarity of phylogenetic profiles, gene expression, protein fusion events and others. Using E. coli and S. cerevisiae metabolic networks, we illustrate predictive ability of each individual type of association evidence and show that significantly better predictions can be obtained based on the combination of all data. In this way our method is able to predict 60% of enzyme-encoding genes of E. coli metabolism within the top 10 (out of 3551) candidates for their enzymatic function, and as a top candidate within 43% of the cases.

Conclusion

We illustrate that a combination of genome context and other functional association evidence is effective in predicting genes encoding metabolic enzymes. Our approach does not rely on direct sequence homology to known enzyme-encoding genes, and can be used in conjunction with traditional homology-based metabolic reconstruction methods. The method can also be used to target orphan metabolic activities.

Gene fusion in Helicobacter pylori: making the ends meet

Fusion genes have been reported as a means of enabling the development of novel or enhanced functions. In this report, we analyzed fusion genes in the genomes of two Helicobacter pylori strains (26695 and J99) and identified 32 fusion genes that are present as neighbours in one strain (components) and are fused in the second (composite), and vice-versa. The mechanism for each case of gene fusion is explored. 28 out of 32 genes identified as fusion products in this analysis were reported as essential genes in the previously documented transposon mutagenesis of H. pylori strain G27. This observation suggests the potential of the products of fusion genes as putative microbial drug targets. These results underscore the utility of bacterial genomic sequence comparisons for understanding gene evolution and for in silico drug target identification in the post-genomic era.

Transcriptional tradeoff between metabolic and stress-response programs in Pseudomonas putida KT2440 cells exposed to toluene

When Pseudomonas putida KT2440 cells encounter toluene in the growth medium, they perceive it simultaneously as a potential nutrient to be metabolized, as a membrane-damaging toxic drug to be extruded, and as a macromolecule-disrupting agent from which to protect proteins. Each of these inputs requires a dedicated transcriptional response that involves a large number of genes. We used DNA array technology to decipher the interplay between these responses in P. putida KT2440 subjected to a short challenge (15 min) with toluene. We then compared the results with those in cells exposed to o-xylene (a non-biodegradable toluene counterpart) and 3-methylbenzoate (a specific substrate of the lower TOL pathway of the P. putida pWW0 plasmid). The resulting expression profiles suggest that the bulk of the available transcriptional machinery is reassigned to endure general stress, whereas only a small share of the available machinery is redirected to the degradation of the aromatic compounds. Specifically, both toluene and o-xylene induce the TOL pathways and a dedicated but not always productive metabolic program. Similarly, 3-methylbenzoate induces the expression not only of the lower meta pathway but also of the non-productive and potentially deleterious genes for the metabolism of (nonsubstituted) benzoate. In addition, toluene (and to a lesser extent o-xylene) inhibit motility functions as an unequivocal response to aromatic toxicity. We argue that toluene is sensed by P. putida KT2440 as a stressor rather than as a nutrient and that the inhibition by the aromatic compounds of many functions we tested is the tradeoff for activating stress tolerance genes at a minimal cost in terms of energy.

Orthologs, paralogs, and evolutionary genomics

Orthologs and paralogs are two fundamentally different types of homologous genes that evolved, respectively, by vertical descent from a single ancestral gene and by duplication. Orthology and paralogy are key concepts of evolutionary genomics. A clear distinction between orthologs and paralogs is critical for the construction of a robust evolutionary classification of genes and reliable functional annotation of newly sequenced genomes. Genome comparisons show that orthologous relationships with genes from taxonomically distant species can be established for the majority of the genes from each sequenced genome. This review examines in depth the definitions and subtypes of orthologs and paralogs, outlines the principal methodological approaches employed for identification of orthology and paralogy, and considers evolutionary and functional implications of these concepts.

Bacterial colonization of nasal mucosa induces expression of siderocalin, an iron-sequestering component of innate immunity

Host-microbe interactions often begin with colonization of mucosal surfaces. These relationships are highly specific, as certain microbial species are found only in particular microenvironments. Transcriptional microarrays were used to screen host genes whose expression in the murine nasal mucosa was affected by colonization with the Gram-positive bacterium Streptococcus pneumoniae. Siderocalin (Scn, or lipocalin 2 or neutrophil gelatinase-associated lipocalin) expression was increased up to 65-fold during colonization by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR). Western analysis showed that Scn was secreted into airway surface fluid in colonized animals. Immunohistochemical analysis localized Scn expression primarily to secretory Bowman's glands. Similar results were observed during colonization with the Gram-negative bacterium Haemophilus influenzae, suggesting that Scn secretion is a general response. Western analysis of human nasal secretions also demonstrated secretion of Scn at potentially bacteriostatic levels. This is a previously unrecognized response that may have a role in determining the establishment or maintenance of mucosal colonization. Scn contributes to antimicrobial defence by sequestration of a subset of microbial siderophores. As neither S. pneumoniae nor H. influenzae are known to produce or utilize siderophores, successful colonizers of the nasal passages may have evolved siderophore-independent mechanisms to acquire essential iron and to evade the inhibitory effects of Scn.