Highly divergent genes for methanopterin-linked C1 transfer reactions in Lake Washington, assessed via metagenomic analysis and mRNA detection

The origins and the evolutionary history of tetrahydromethanopterin-linked C1 transfer reactions that are part of two environmentally important biotransformations, methylotrophy and methanogenesis, are still not well understood. In previous studies, we have expanded the known phylogenetic diversity of these reactions by identifying genes highly diverging from the ones associated with cultivated Proteobacteria, Planctomycetes, or Archaea (M. G. Kalyuzhnaya, M. E. Lidstrom, and L. Chistoserdova, Microb. Ecol. 48:463-472, 2004; M. G. Kalyuzhnaya, O. Nercessian, M. E. Lidstrom, and L. Chistoserdova, Environ. Microbiol. 7:1269-1274, 2005). Here we used a metagenomic approach to demonstrate that these divergent genes are present with high abundance in the microbial community inhabiting Lake Washington sediment. We also gained preliminary insights into the genomic composition of the organisms possessing these genes by sequencing genomic fragments from three uncultured microbes possessing the genes of interest. Phylogenetic analyses suggested that, although distantly related to each other, these organisms deeply diverge from known Bacteria and Archaea, with more relation to the former, suggesting their affiliation with a new bacterial phylum. We also demonstrate, via specific mRNA detection, that these divergent genes are expressed in the environment, pointing toward their potential role in local carbon cycling.     

Poly-beta-hydroxybutyrate biosynthesis in the facultative methylotroph methylobacterium extorquens AM1: identification and mutation of gap11, gap20, and phaR

Methylobacterium extorquens AM1, a serine cycle facultative methylotroph, accumulates poly-beta-hydroxybutyrate (PHB) as a carbon and energy reserve material during growth on both multicarbon- and single-carbon substrates. Recently, the identification and mutation of the genes involved in the biosynthesis and degradation of PHB have been described for this bacterium, demonstrating that two of the genes of the PHB cycle (phaA and phaB) are also involved in C(1) and C(2) metabolism, as part of a novel pathway for glyoxylate regeneration in the serine cycle (N. Korotkova and M. E. Lidstrom, J. Bacteriol. 183:1038-1046, 2001; N. Korotkova, L. Chistoserdova, V. Kuksa, and M. E. Lidstrom, J. Bacteriol. 184:1750-1758, 2002). In this work, three new genes involved in PHB biosynthesis in this bacterium have been investigated via mutation and phenotypic analysis: gap11, gap20, and phaR. We demonstrate that gap11 and gap20 encode two major granule-associated proteins (phasins) and that mutants with mutations in these genes are defective in PHB production and also in growth on C(2) compounds, while they show wild-type growth characteristics on C(1) or multicarbon compounds. The phaR mutant shows defects in both PHB accumulation and growth characteristics when grown on C(1) compounds and has defects in PHB accumulation but grows normally on C(3) and C(4) compounds, while both PHB accumulation and growth rate are at wild-type levels during growth on C(2) compounds. Our results suggest that this phenotype is due to altered fluxes of acetyl coenzyme A (CoA), a major intermediate in C(1), C(2), and heterotrophic metabolism in M. extorquens AM1, as well as the entry metabolite for the PHB cycle. Therefore, it seems likely that PhaR acts to control acetyl-CoA flux to PHB in this methylotrophic bacterium.     

Identification of a fourth formate dehydrogenase in Methylobacterium extorquens AM1 and confirmation of the essential role of formate oxidation in methylotrophy

A mutant of Methylobacterium extorquens AM1 with lesions in genes for three formate dehydrogenase (FDH) enzymes was previously described by us (L. Chistoserdova, M. Laukel, J.-C. Portais, J. A. Vorholt, and M. E. Lidstrom, J. Bacteriol. 186:22-28, 2004). This mutant had lost its ability to grow on formate but still maintained the ability to grow on methanol. In this work, we further investigated the phenotype of this mutant. Nuclear magnetic resonance experiments with [¹³C]formate, as well as ¹⁴C-labeling experiments, demonstrated production of labeled CO₂ in the mutant, pointing to the presence of an additional enzyme or a pathway for formate oxidation. The tungsten-sensitive phenotype of the mutant suggested the involvement of a molybdenum-dependent enzyme. Whole-genome array experiments were conducted to test for genes overexpressed in the triple-FDH mutant compared to the wild type, and a gene (fdh4A) was identified whose translated product carried similarity to an uncharacterized putative molybdopterin-binding oxidoreductase-like protein sharing relatively low similarity with known formate dehydrogenase alpha subunits. Mutation of this gene in the triple-FDH mutant background resulted in a methanol-negative phenotype. When the gene was deleted in the wild-type background, the mutant revealed diminished growth on methanol with accumulation of high levels of formate in the medium, pointing to an important role of FDH4 in methanol metabolism. The identity of FDH4 as a novel FDH was also confirmed by labeling experiments that revealed strongly reduced CO₂ formation in growing cultures. Mutation of a small open reading frame (fdh4B) downstream of fdh4A resulted in mutant phenotypes similar to the phenotypes of fdh4A mutants, suggesting that fdh4B is also involved in formate oxidation.     

The Exopolyphosphatase Gene from Sulfolobus solfataricus: Characterization of the First Gene Found To Be Involved in Polyphosphate Metabolism in Archaea

Inorganic polyphosphate (polyP) polymers are widely distributed in all kinds of organisms. Although the presence of polyP in members of the domain Archaea has been described, at present nothing is known about the enzymology of polyP metabolism or the genes involved in this domain. We have cloned, sequenced, and overexpressed an exopolyphosphatase (PPX) gene (ppx) from thermophilic Sulfolobus solfataricus. The gene codes for a functional PPX and possesses an open reading frame for 417 amino acids (calculated mass, 47.9 kDa). The purified recombinant PPX was highly active, degrading long-chain polyP (700 to 800 residues) in vitro at 50 to 60°C. The putative PPXs present in known archaeal genomes showed the highest similarity to yeast PPXs. In contrast, informatic analysis revealed that the deduced amino acid sequence of S. solfataricus PPX showed the highest similarity (25 to 45%) to sequences of members of the bacterial PPXs, possessing all of their conserved motifs. To our knowledge, this is the first report of an enzyme characterized to be involved in polyP metabolism in members of the Archaea.     

A survey of small RNA-encoding genes in Escherichia coli

Small RNA (sRNA) molecules have gained much interest lately, as recent genome-wide studies have shown that they are widespread in a variety of organisms. The relatively small family of 10 known sRNA-encoding genes in Escherichia coli has been significantly expanded during the past two years with the discovery of 45 novel genes. Most of these genes are still uncharacterized and their cellular roles are unknown. In this survey we examined the sequence and genomic features of the 55 currently known sRNA-encoding genes in E.coli, attempting to identify their common characteristics. Such characterization is important for both expanding our understanding of this unique gene family and for improving the methods to predict and identify sRNA-encoding genes based on genomic information.     

Gut Symbionts from Distinct Hosts Exhibit Genotoxic Activity via Divergent Colibactin Biosynthesis Pathways

Secondary metabolites produced by nonribosomal peptide synthetase (NRPS) or polyketide synthase (PKS) pathways are chemical mediators of microbial interactions in diverse environments. However, little is known about their distribution, evolution, and functional roles in bacterial symbionts associated with animals. A prominent example is colibactin, a largely unknown family of secondary metabolites produced by Escherichia coli via a hybrid NRPS-PKS biosynthetic pathway that inflicts DNA damage upon eukaryotic cells and contributes to colorectal cancer and tumor formation in the mammalian gut. Thus far, homologs of this pathway have only been found in closely related Enterobacteriaceae, while a divergent variant of this gene cluster was recently discovered in a marine alphaproteobacterial Pseudovibrio strain. Herein, we sequenced the genome of Frischella perrara PEB0191, a bacterial gut symbiont of honey bees and identified a homologous colibactin biosynthetic pathway related to those found in Enterobacteriaceae. We show that the colibactin genomic island (GI) has conserved gene synteny and biosynthetic module architecture across F. perrara, Enterobacteriaceae, and the Pseudovibrio strain. Comparative metabolomics analyses of F. perrara and E. coli further reveal that these two bacteria produce related colibactin pathway-dependent metabolites. Finally, we demonstrate that F. perrara, like E. coli, causes DNA damage in eukaryotic cells in vitro in a colibactin pathway-dependent manner. Together, these results support that divergent variants of the colibactin biosynthetic pathway are widely distributed among bacterial symbionts, producing related secondary metabolites and likely endowing its producer with functional capabilities important for diverse symbiotic associations.     

Genetics of the serine cycle in Methylobacterium extorquens AM1: identification, sequence, and mutation of three new genes involved in C1 assimilation, orf4, mtkA, and mtkB.

In a recent paper we reported the sequence of the beginning of a serine cycle gene cluster on the Methylobacterium extorquens AM1 chromosome, containing the genes encoding serine glyoxylate aminotransferase (sgaA), hydroxypyruvate reductase (hprA), and 5,10-methylenetetrahydrofolate dehydrogenase (mtdA) (L. V. Chistoserdova and M. E. Lidstrom J. Bacteriol. 176:1957-1968, 1994). Here we present the sequence of the adjacent downstream region containing three full and one partial open reading frames. The first of the full open reading frames (orf4) remains unidentified, while the other two (mtkA and mtkB) code for the two subunits of malate thiokinase, and the fourth, a partial open reading frame (ppcA), apparently encodes phosphoenolpyruvate carboxylase. Mutants containing insertion mutations in orf4, mtdA, and mtdB all were unable to grow on C1 compounds, showing that these three newly identified genes are indispensable for the operation of the serine cycle. Mutants in orf4 were also unable to grow on C2 compounds, but growth was restored by glyoxylate, suggesting that orf4 might be required for the conversion of acetyl coenzyme A to glyoxylate.     

Development of Antibiotic-Overproducing Strains by Site-Directed Mutagenesis of the rpsL Gene in Streptomyces lividans

Certain rpsL (which encodes the ribosomal protein S12) mutations that confer resistance to streptomycin markedly activate the production of antibiotics in Streptomyces spp. These rpsL mutations are known to be located in the two conserved regions within the S12 protein. To understand the roles of these two regions in the activation of silent genes, we used site-directed mutagenesis to generate eight novel mutations in addition to an already known (K88E) mutation that is capable of activating antibiotic production in Streptomyces lividans. Of these mutants, two (L90K and R94G) activated antibiotic production much more than the K88E mutant. Neither the L90K nor the R94G mutation conferred an increase in the level of resistance to streptomycin and paromomycin. Our results demonstrate the efficacy of the site-directed mutagenesis technique for strain improvement.     

Whole-Genome Reciprocal BLAST Analysis Reveals that Planctomycetes Do Not Share an Unusually Large Number of Genes with Eukarya and Archaea

The genome sequences of Rhodopirellula baltica, formerly Pirellula sp. strain 1, Blastopirellula marina, Gemmata obscuriglobus, and Kuenenia stuttgartiensis were used in a series of pairwise reciprocal best-hit analyses to evaluate the contested evolutionary position of Planctomycetes. Contrary to previous reports which suggested that R. baltica had a high percentage of genes with closest matches to Archaea and Eukarya, we show here that these Planctomycetes do not share an unusually large number of genes with the Archaea or Eukarya, compared with other Bacteria. Thus, best-hit analyses may assign phylogenetic affinities incorrectly if close relatives are absent from the sequence database.     

Comparative Genomic Hybridization (CGH) Reveals a Neo-X Chromosome and Biased Gene Movement in Stalk-Eyed Flies (Genus Teleopsis)

Chromosomal location has a significant effect on the evolutionary dynamics of genes involved in sexual dimorphism, impacting both the pattern of sex-specific gene expression and the rate of duplication and protein evolution for these genes. For nearly all non-model organisms, however, knowledge of chromosomal gene content is minimal and difficult to obtain on a genomic scale. In this study, we utilized Comparative Genomic Hybridization (CGH), using probes designed from EST sequence, to identify genes located on the X chromosome of four species in the stalk-eyed fly genus Teleopsis. Analysis of log₂ ratio values of female-to-male hybridization intensities from the CGH microarrays for over 3,400 genes reveals a strongly bimodal distribution that clearly differentiates autosomal from X-linked genes for all four species. Genotyping of 33 and linkage mapping of 28 of these genes in Teleopsis dalmanni indicate the CGH results correctly identified chromosomal location in all cases. Syntenic comparison with Drosophila indicates that 90% of the X-linked genes in Teleopsis are homologous to genes located on chromosome 2L in Drosophila melanogaster, suggesting the formation of a nearly complete neo-X chromosome from Muller element B in the dipteran lineage leading to Teleopsis. Analysis of gene movement both relative to Drosophila and within Teleopsis indicates that gene movement is significantly associated with 1) rates of protein evolution, 2) the pattern of gene duplication, and 3) the evolution of eyespan sexual dimorphism. Overall, this study reveals that diopsids are a critical group for understanding the evolution of sex chromosomes within Diptera. In addition, we demonstrate that CGH is a useful technique for identifying chromosomal sex-linkage and should be applicable to other organisms with EST or partial genomic information.     

Methane- and Sulfur-Metabolizing Microbial Communities Dominate the Lost City Hydrothermal Field Ecosystem

Hydrothermal venting and the formation of carbonate chimneys in the Lost City hydrothermal field (LCHF) are driven predominantly by serpentinization reactions and cooling of mantle rocks, resulting in a highly reducing, high-pH environment with abundant dissolved hydrogen and methane. Phylogenetic and terminal restriction fragment length polymorphism analyses of 16S rRNA genes in fluids and carbonate material from this site indicate the presence of organisms similar to sulfur-oxidizing, sulfate-reducing, and methane-oxidizing Bacteria as well as methanogenic and anaerobic methane-oxidizing Archaea. The presence of these metabolic groups indicates that microbial cycling of sulfur and methane may be the dominant biogeochemical processes active within this ultramafic rock-hosted environment. 16S rRNA gene sequences grouping within the Methylobacter and Thiomicrospira clades were recovered from a chemically diverse suite of carbonate chimney and fluid samples. In contrast, 16S rRNA genes corresponding to the Lost City Methanosarcinales phylotype were found exclusively in high-temperature chimneys, while a phylotype of anaerobic methanotrophic Archaea (ANME-1) was restricted to lower-temperature, less vigorously venting sites. A hyperthermophilic habitat beneath the LCHF may be reflected by 16S rRNA gene sequences belonging to Thermococcales and uncultured Crenarchaeota identified in vent fluids. The finding of a diverse microbial ecosystem supported by the interaction of high-temperature, high-pH fluids resulting from serpentinization reactions in the subsurface provides insight into the biogeochemistry of what may be a pervasive process in ultramafic subseafloor environments.     

Major Soybean Maturity Gene Haplotypes Revealed by SNPViz Analysis of 72 Sequenced Soybean Genomes

In this Genomics Era, vast amounts of next-generation sequencing data have become publicly available for multiple genomes across hundreds of species. Analyses of these large-scale datasets can become cumbersome, especially when comparing nucleotide polymorphisms across many samples within a dataset and among different datasets or organisms. To facilitate the exploration of allelic variation and diversity, we have developed and deployed an in-house computer software to categorize and visualize these haplotypes. The SNPViz software enables users to analyze region-specific haplotypes from single nucleotide polymorphism (SNP) datasets for different sequenced genomes. The examination of allelic variation and diversity of important soybean [Glycine max (L.) Merr.] flowering time and maturity genes may provide additional insight into flowering time regulation and enhance researchers'' ability to target soybean breeding for particular environments. For this study, we utilized two available soybean genomic datasets for a total of 72 soybean genotypes encompassing cultivars, landraces, and the wild species Glycine soja. The major soybean maturity genes E1, E2, E3, and E4 along with the Dt1 gene for plant growth architecture were analyzed in an effort to determine the number of major haplotypes for each gene, to evaluate the consistency of the haplotypes with characterized variant alleles, and to identify evidence of artificial selection. The results indicated classification of a small number of predominant haplogroups for each gene and important insights into possible allelic diversity for each gene within the context of known causative mutations. The software has both a stand-alone and web-based version and can be used to analyze other genes, examine additional soybean datasets, and view similar genome sequence and SNP datasets from other species.     

Correlation of glucocorticoid-mediated E4BP4 upregulation with altered expression of pro- and anti-apoptotic genes in CEM human lymphoblastic leukemia cells

In Caenorhabditiselegans, motorneuron apoptosis is regulated via a ces-2–ces-1–egl-1 pathway. We tested whether human CEM lymphoblastic leukemia cells undergo apoptosis via an analogous pathway. We have previously shown that E4BP4, a ces-2 ortholog, mediates glucocorticoid (GC)-dependent upregulation of BIM, an egl-1 ortholog, in GC-sensitive CEM C7-14 cells and in CEM C1-15mE#3 cells, which are sensitized to GCs by ectopic expression of E4BP4. In the present study, we demonstrate that the human ces-1 orthologs, SLUG and SNAIL, are not significantly repressed in correlation with E4BP4 expression. Expression of E4BP4 homologs, the PAR family genes, especially HLF, encoding a known anti-apoptotic factor, was inverse to that of E4BP4 and BIM. Expression of pro- and anti-apoptotic genes in CEM cells was analyzed via an apoptosis PCR Array. We identified BIRC3 and BIM as genes whose expression paralleled that of E4BP4, while FASLG, TRAF4, BCL2A1, BCL2L1, BCL2L2 and CD40LG as genes whose expression was opposite to that of E4BP4.     

Distribution and diversity of rhizobia associated with wild soybean (Glycine soja Sieb. & Zucc.) in Northwest China

A total of 155 nodule isolates that originated from seven sites in Northwest China were characterized by PCR-RFLP of the 16S rRNA gene and sequence analysis of multiple core genes (16S rRNA, recA, atpD, and glnII) in order to investigate the diversity and biogeography of Glycine soja-nodulating rhizobia. Among the isolates, 80 were Ensifer fredii, 19 were Ensifer morelense, 49 were Rhizobium radiobacter, and 7 were putative novel Rhizobium species. The phylogenies of E. fredii and E. morelense isolates in a concatenate tree (assembly of all housekeeping genes) were generally consistent with those in individual gene trees. However, incongruence was found in the phylogenies of the different genes of Rhizobium isolates, indicating that lateral transfer or recombination possibly occurred in these gene loci. Despite their species identity, all the isolates in this study formed a single lineage related to E. fredii in nodAand nifH gene phylogenies, which also indicated that the symbiotic genes were laterally transferred between different species. Biogeographic patterns were found at the species and strain genomic type levels, as revealed by BOXA1R fingerprinting, demonstrating that the evolution of rhizobial populations in different geographic locations was related to soil types, altitude and spatial effects. This study is the first to report that E. morelense, R. radiobacter, and Rhizobium sp. are microsymbionts of G. soja, as well as showing that the diversity of G. soja rhizobia is enhanced and new rhizobia have evolved in Northwest China.     

Comparative Genomic Analysis Reveals 2-Oxoacid Dehydrogenase Complex Lipoylation Correlation with Aerobiosis in Archaea

Metagenomic analyses have advanced our understanding of ecological microbial diversity, but to what extent can metagenomic data be used to predict the metabolic capacity of difficult-to-study organisms and their abiotic environmental interactions? We tackle this question, using a comparative genomic approach, by considering the molecular basis of aerobiosis within archaea. Lipoylation, the covalent attachment of lipoic acid to 2-oxoacid dehydrogenase multienzyme complexes (OADHCs), is essential for metabolism in aerobic bacteria and eukarya. Lipoylation is catalysed either by lipoate protein ligase (LplA), which in archaea is typically encoded by two genes (LplA-N and LplA-C), or by a lipoyl(octanoyl) transferase (LipB or LipM) plus a lipoic acid synthetase (LipA). Does the genomic presence of lipoylation and OADHC genes across archaea from diverse habitats correlate with aerobiosis? First, analyses of 11,826 biotin protein ligase (BPL)-LplA-LipB transferase family members and 147 archaeal genomes identified 85 species with lipoylation capabilities and provided support for multiple ancestral acquisitions of lipoylation pathways during archaeal evolution. Second, with the exception of the Sulfolobales order, the majority of species possessing lipoylation systems exclusively retain LplA, or either LipB or LipM, consistent with archaeal genome streamlining. Third, obligate anaerobic archaea display widespread loss of lipoylation and OADHC genes. Conversely, a high level of correspondence is observed between aerobiosis and the presence of LplA/LipB/LipM, LipA and OADHC E2, consistent with the role of lipoylation in aerobic metabolism. This correspondence between OADHC lipoylation capacity and aerobiosis indicates that genomic pathway profiling in archaea is informative and that well characterized pathways may be predictive in relation to abiotic conditions in difficult-to-study extremophiles. Given the highly variable retention of gene repertoires across the archaea, the extension of comparative genomic pathway profiling to broader metabolic and homeostasis networks should be useful in revealing characteristics from metagenomic datasets related to adaptations to diverse environments.     

External Loops at the C Terminus of Erwinia chrysanthemi Pectate Lyase C Are Required for Species-Specific Secretion through the Out Type II Pathway

The type II secretion system (main terminal branch of the general secretion pathway) is used by diverse gram-negative bacteria to secrete extracellular proteins. Proteins secreted by this pathway are synthesized with an N-terminal signal peptide which is removed upon translocation across the inner membrane, but the signals which target the mature proteins for secretion across the outer membrane are unknown. The plant pathogens Erwinia chrysanthemi and Erwinia carotovora secrete several isozymes of pectate lyase (Pel) by the out-encoded type II pathway. However, these two bacteria cannot secrete Pels encoded by heterologously expressed pel genes from the other species, suggesting the existence of species-specific secretion signals within these proteins. The functional cluster of E. chrysanthemi out genes carried on cosmid pCPP2006 enables Escherichia coli to secrete E. chrysanthemi, but not E. carotovora, Pels. We exploited the high sequence similarity between E. chrysanthemi PelC and E. carotovora Pel1 to construct 15 hybrid proteins in which different regions of PelC were replaced with homologous sequences from Pel1. The differential secretion of these hybrid proteins by E. coli(pCPP2006) revealed M118 to D175 and V215 to C329 as regions required for species-specific secretion of PelC. We propose that the primary targeting signal is contained within the external loops formed by G274 to C329 but is dependent on residues in M118 to D170 and V215 to G274 for proper positioning.     

Phylogeny of Dissimilatory Sulfite Reductases Supports an Early Origin of Sulfate Respiration

Microorganisms that use sulfate as a terminal electron acceptor for anaerobic respiration play a central role in the global sulfur cycle. Here, we report the results of comparative sequence analysis of dissimilatory sulfite reductase (DSR) genes from closely and distantly related sulfate-reducing organisms to infer the evolutionary history of DSR. A 1.9-kb DNA region encoding most of the α and β subunits of DSR could be recovered only from organisms capable of dissimilatory sulfate reduction with a PCR primer set targeting highly conserved regions in these genes. All DNA sequences obtained were highly similar to one another (49 to 89% identity), and their inferred evolutionary relationships were nearly identical to those inferred on the basis of 16S rRNA. We conclude that the high similarity of bacterial and archaeal DSRs reflects their common origin from a conserved DSR. This ancestral DSR was either present before the split between the domains Bacteria, Archaea, and Eucarya or laterally transferred between Bacteria and Archaea soon after domain divergence. Thus, if the physiological role of the DSR was constant over time, then early ancestors of Bacteria and Archaea already possessed a key enzyme of sulfate and sulfite respiration.     

Genetics of the serine cycle in Methylobacterium extorquens AM1: identification of sgaA and mtdA and sequences of sgaA, hprA, and mtdA.

In a previous paper, we reported identification of the 5' part of hprA of Methylobacterium extorquens AM1, which encodes the serine cycle enzyme hydroxypyruvate reductase (L. V. Chistoserdova and M. E. Lidstrom, J. Bacteriol. 174:71-77, 1992). Here we present the complete sequence of hprA and partial sequence of genes adjacent to hprA. Upstream of hprA, the 3' part of an open reading frame was discovered, separated from hprA by 263 bp. This open reading frame was identified as the gene encoding another serine cycle enzyme, serine glyoxylate aminotransferase (sgaA). Cells containing an insertion mutation into sgaA were unable to grow on C1 compounds, demonstrating that the gene is required for C1 metabolism. Sequencing downstream of hprA has revealed the presence of another open reading frame (mtdA), which is probably cotranscribed with hprA. This open reading frame was identified as the gene required for the synthesis of 5,10-methylenetetrahydrofolate dehydrogenase. Our data suggest that this enzyme plays an integral role in methylotrophic metabolism in M. extorquens AM1, either in formaldehyde oxidation or as part of the serine cycle.     

Activity-based cell sorting reveals responses of uncultured archaea and bacteria to substrate amendment

Metagenomic studies have revolutionized our understanding of the metabolic potential of uncultured microorganisms in various ecosystems. However, many of these genomic predictions have yet to be experimentally tested, and the functional expression of genomic potential often remains unaddressed. In order to obtain a more thorough understanding of cell physiology, novel techniques capable of testing microbial metabolism under close to in situ conditions must be developed. Here, we provide a benchmark study to demonstrate that bioorthogonal non-canonical amino acid tagging (BONCAT) in combination with fluorescence-activated cell sorting (FACS) and 16S rRNA gene sequencing can be used to identify anabolically active members of a microbial community incubated in the presence of various growth substrates or under changing physicochemical conditions. We applied this approach to a hot spring sediment microbiome from Yellowstone National Park (Wyoming, USA) and identified several microbes that changed their activity levels in response to substrate addition, including uncultured members of the phyla Thaumarchaeota, Acidobacteria, and Fervidibacteria. Because shifts in activity in response to substrate amendment or headspace changes are indicative of microbial preferences for particular growth conditions, results from this and future BONCAT-FACS studies could inform the development of cultivation media to specifically enrich uncultured microbes. Most importantly, BONCAT-FACS is capable of providing information on the physiology of uncultured organisms at as close to in situ conditions as experimentally possible.Subject terms: Environmental microbiology, Microbial communities, Microbial ecology, Archaea     

Novel Conserved Genotypes Correspond to Antibiotic Resistance Phenotypes of E. coli Clinical Isolates

Current efforts to understand antibiotic resistance on the whole genome scale tend to focus on known genes even as high throughput sequencing strategies uncover novel mechanisms. To identify genomic variations associated with antibiotic resistance, we employed a modified genome-wide association study; we sequenced genomic DNA from pools of E. coli clinical isolates with similar antibiotic resistance phenotypes using SOLiD technology to uncover single nucleotide polymorphisms (SNPs) unanimously conserved in each pool. The multidrug-resistant pools were genotypically similar to SMS-3-5, a previously sequenced multidrug-resistant isolate from a polluted environment. The similarity was evenly spread across the entire genome and not limited to plasmid or pathogenicity island loci. Among the pools of clinical isolates, genomic variation was concentrated adjacent to previously reported inversion and duplication differences between the SMS-3-5 isolate and the drug-susceptible laboratory strain, DH10B. SNPs that result in non-synonymous changes in gyrA (encoding the well-known S83L allele associated with fluoroquinolone resistance), mutM, ligB, and recG were unanimously conserved in every fluoroquinolone-resistant pool. Alleles of the latter three genes are tightly linked among most sequenced E. coli genomes, and had not been implicated in antibiotic resistance previously. The changes in these genes map to amino acid positions in alpha helices that are involved in DNA binding. Plasmid-encoded complementation of null strains with either allelic variant of mutM or ligB resulted in variable responses to ultraviolet light or hydrogen peroxide treatment as markers of induced DNA damage, indicating their importance in DNA metabolism and revealing a potential mechanism for fluoroquinolone resistance. Our approach uncovered evidence that additional DNA binding enzymes may contribute to fluoroquinolone resistance and further implicate environmental bacteria as a reservoir for antibiotic resistance.