Anaerobically Grown Escherichia coli Has an Enhanced Mutation Rate and Distinct Mutational Spectra.

Oxidative stress is a major cause of mutation but little is known about how growth in the absence of oxygen impacts the rate and spectrum of mutations. We employed long-term mutation accumulation experiments to directly measure the rates and spectra of spontaneous mutation events in Escherichia coli populations propagated under aerobic and anaerobic conditions. To detect mutations, whole genome sequencing was coupled with methods of analysis sufficient to identify a broad range of mutational classes, including structural variants (SVs) generated by movement of repetitive elements. The anaerobically grown populations displayed a mutation rate nearly twice that of the aerobic populations, showed distinct asymmetric mutational strand biases, and greater insertion element activity. Consistent with mutation rate and spectra observations, genes for transposition and recombination repair associated with SVs were up-regulated during anaerobic growth. Together, these results define differences in mutational spectra affecting the evolution of facultative anaerobes.

Comparative Genomics of Rumen Butyrivibrio spp. Uncovers a Continuum of Polysaccharide-Degrading Capabilities.

Plant polysaccharide breakdown by microbes in the rumen is fundamental to digestion in ruminant livestock. Bacterial species belonging to the rumen genera Butyrivibrio and Pseudobutyrivibrio are important degraders and utilizers of lignocellulosic plant material. These bacteria degrade polysaccharides and ferment the released monosaccharides to yield short-chain fatty acids that are used by the ruminant for growth and the production of meat, milk, and fiber products. Although rumen Butyrivibrio and Pseudobutyrivibrio species are regarded as common rumen inhabitants, their polysaccharide-degrading and carbohydrate-utilizing enzymes are not well understood. In this study, we analyzed the genomes of 40 Butyrivibrio and 6 Pseudobutyrivibrio strains isolated from the plant-adherent fraction of New Zealand dairy cows to explore the polysaccharide-degrading potential of these important rumen bacteria. Comparative genome analyses combined with phylogenetic analysis of their 16S rRNA genes and short-chain fatty acid production patterns provide insight into the genomic diversity and physiology of these bacteria and divide Butyrivibrio into 3 species clusters. Rumen Butyrivibrio bacteria were found to encode a large and diverse spectrum of degradative carbohydrate-active enzymes (CAZymes) and binding proteins. In total, 4,421 glycoside hydrolases (GHs), 1,283 carbohydrate esterases (CEs), 110 polysaccharide lyases (PLs), 3,605 glycosyltransferases (GTs), and 1,706 carbohydrate-binding protein modules (CBM) with predicted activities involved in the depolymerization and transport of the insoluble plant polysaccharides were identified. Butyrivibrio genomes had similar patterns of CAZyme families but varied greatly in the number of genes within each category in the Carbohydrate-Active Enzymes database (CAZy), suggesting some level of functional redundancy. These results suggest that rumen Butyrivibrio species occupy similar niches but apply different degradation strategies to be able to coexist in the rumen.IMPORTANCE Feeding a global population of 8 billion people and climate change are the primary challenges facing agriculture today. Ruminant livestock are important food-producing animals, and maximizing their productivity requires an understanding of their digestive systems and the roles played by rumen microbes in plant polysaccharide degradation. Members of the genera Butyrivibrio and Pseudobutyrivibrio are a phylogenetically diverse group of bacteria and are commonly found in the rumen, where they are a substantial source of polysaccharide-degrading enzymes for the depolymerization of lignocellulosic material. Our findings have highlighted the immense enzymatic machinery of Butyrivibrio and Pseudobutyrivibrio species for the degradation of plant fiber, suggesting that these bacteria occupy similar niches but apply different degradation strategies in order to coexist in the competitive rumen environment.

Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome.

Ruminant livestock represent the single largest anthropogenic source of the potent greenhouse gas methane, which is generated by methanogenic archaea residing in ruminant digestive tracts. While differences between individual animals of the same breed in the amount of methane produced have been observed, the basis for this variation remains to be elucidated. To explore the mechanistic basis of this methane production, we measured methane yields from 22 sheep, which revealed that methane yields are a reproducible, quantitative trait. Deep metagenomic and metatranscriptomic sequencing demonstrated a similar abundance of methanogens and methanogenesis pathway genes in high and low methane emitters. However, transcription of methanogenesis pathway genes was substantially increased in sheep with high methane yields. These results identify a discrete set of rumen methanogens whose methanogenesis pathway transcription profiles correlate with methane yields and provide new targets for CH4 mitigation at the levels of microbiota composition and transcriptional regulation.

Metasecretome-selective phage display approach for mining the functional potential of a rumen microbial community.

BACKGROUND: In silico, secretome proteins can be predicted from completely sequenced genomes using various available algorithms that identify membrane-targeting sequences. For metasecretome (collection of surface, secreted and transmembrane proteins from environmental microbial communities) this approach is impractical, considering that the metasecretome open reading frames (ORFs) comprise only 10% to 30% of total metagenome, and are poorly represented in the dataset due to overall low coverage of metagenomic gene pool, even in large-scale projects. RESULTS: By combining secretome-selective phage display and next-generation sequencing, we focused the sequence analysis of complex rumen microbial community on the metasecretome component of the metagenome. This approach achieved high enrichment (29 fold) of secreted fibrolytic enzymes from the plant-adherent microbial community of the bovine rumen. In particular, we identified hundreds of heretofore rare modules belonging to cellulosomes, cell-surface complexes specialised for recognition and degradation of the plant fibre. CONCLUSIONS: As a method, metasecretome phage display combined with next-generation sequencing has a power to sample the diversity of low-abundance surface and secreted proteins that would otherwise require exceptionally large metagenomic sequencing projects. As a resource, metasecretome display library backed by the dataset obtained by next-generation sequencing is ready for i) affinity selection by standard phage display methodology and ii) easy purification of displayed proteins as part of the virion for individual functional analysis.

Functional genome analysis of Bifidobacterium breve UCC2003 reveals type IVb tight adherence (Tad) pili as an essential and conserved host-colonization factor.

Development of the human gut microbiota commences at birth, with bifidobacteria being among the first colonizers of the sterile newborn gastrointestinal tract. To date, the genetic basis of Bifidobacterium colonization and persistence remains poorly understood. Transcriptome analysis of the Bifidobacterium breve UCC2003 2.42-Mb genome in a murine colonization model revealed differential expression of a type IVb tight adherence (Tad) pilus-encoding gene cluster designated "tad(2003)." Mutational analysis demonstrated that the tad(2003) gene cluster is essential for efficient in vivo murine gut colonization, and immunogold transmission electron microscopy confirmed the presence of Tad pili at the poles of B. breve UCC2003 cells. Conservation of the Tad pilus-encoding locus among other B. breve strains and among sequenced Bifidobacterium genomes supports the notion of a ubiquitous pili-mediated host colonization and persistence mechanism for bifidobacteria.

Complete genome sequence of Methanosphaera sp. ISO3-F5, a rumen methylotrophic methanogen.

Methanosphaera spp. are methylotrophic methanogenic archaea and members of the order Methanobacteriales with few cultured representatives. Methanosphaera sp. ISO3-F5 was isolated from sheep rumen contents in New Zealand. Here, we report its complete genome, consisting of a large chromosome and a megaplasmid (GenBank accession numbers CP118753 and CP118754, respectively).

Methanobrevibacter boviskoreani JH1T growth on alcohols allows development of a high throughput bioassay to detect methanogen inhibition.

Rumen methanogenic archaea use by-products of fermentation to carry out methanogenesis for energy generation. A key fermentation by-product is hydrogen (H2), which acts as the source of reducing potential for methane (CH4) formation in hydrogenotrophic methanogens. The in vitro cultivation of hydrogenotrophic rumen methanogens requires pressurised H2 which limits the ability to conduct high-throughput screening experiments with these organisms. The genome of the hydrogenotrophic methanogen Methanobrevibacter boviskoreani JH1T harbors genes encoding an NADP-dependent alcohol dehydrogenase and a F420-dependent NADP reductase, which may facilitate the transfer of reducing potential from ethanol to F420 via NADP. The aim of this study was to explore the anaerobic culturing of JH1T without pressurised H2, using a variety of short chain alcohols. The results demonstrate that in the absence of H2, JHIT can use ethanol, 1-propanol, and 1-butanol but not methanol, as a source of reducing potential for methanogenesis. The ability to use ethanol to drive CH4 formation in JH1T makes it possible to develop a high throughput culture-based bioassay enabling screening of potential anti-methanogen compounds. The development of this resource will help researchers globally to accelerate the search for methane mitigation technologies for ruminant animals. Global emissions pathways that are consistent with the temperature goal of the Paris Agreement, rely on substantial reductions of agricultural greenhouse gasses, particularly from ruminant animals.

Extracellular polysaccharide-degrading proteome of Butyrivibrio proteoclasticus.

Plant polysaccharide-degrading rumen microbes are fundamental to the health and productivity of ruminant animals. Butyrivibrio proteoclasticus B316(T) is a gram-positive, butyrate-producing anaerobic bacterium with a key role in hemicellulose degradation in the rumen. Gel-based proteomics was used to examine the growth-phase-dependent abundance patterns of secreted proteins recovered from cells grown in vitro with xylan or xylose provided as the sole supplementary carbon source. Five polysaccharidases and two carbohydrate-binding proteins (CBPs) were among 30 identified secreted proteins. The endo-1,4-ß-xylanase Xyn10B was 17.5-fold more abundant in the culture medium of xylan-grown cells, which suggests it plays an important role in hemicellulose degradation. The secretion of three nonxylanolytic enzymes and two CBPs implies they augment hemicellulose degradation by hydrolysis or disruption of associated structural polysaccharides. Sixteen ATP-binding cassette (ABC) transporter substrate-binding proteins were identified, several of which had altered relative abundance levels between growth conditions, which suggests they are important for oligosaccharide uptake. This study demonstrates that B. proteoclasticus modulates the secretion of hemicellulose-degrading enzymes and ATP-dependent sugar uptake systems in response to growth substrate and supports the notion that this organism makes an important contribution to polysaccharide degradation in the rumen.

Hydrogen and formate production and utilisation in the rumen and the human colon.

Molecular hydrogen (H2) and formate (HCOO-) are metabolic end products of many primary fermenters in the mammalian gut. Both play a vital role in fermentation where they are electron sinks for individual microbes in an anaerobic environment that lacks external electron acceptors. If H2 and/or formate accumulate within the gut ecosystem, the ability of primary fermenters to regenerate electron carriers may be inhibited and microbial metabolism and growth disrupted. Consequently, H2- and/or formate-consuming microbes such as methanogens and homoacetogens play a key role in maintaining the metabolic efficiency of primary fermenters. There is increasing interest in identifying approaches to manipulate mammalian gut environments for the benefit of the host and the environment. As H2 and formate are important mediators of interspecies interactions, an understanding of their production and utilisation could be a significant entry point for the development of successful interventions. Ruminant methane mitigation approaches are discussed as a model to help understand the fate of H2 and formate in gut systems.

Electron flow: key to mitigating ruminant methanogenesis.

Disposal of electrons generated during the fermentation of ingested feed is a fundamental feature of anaerobic microbial gut ecosystems. Here, we focus on the well-studied rumen environment to highlight how electrons are transferred through anaerobic fermentation pathways and how manipulating this electron flow is important to reducing methane emissions from ruminants. Priorities for research that can accelerate understanding in this area are highlighted.

Genomic insights into the physiology of Quinella, an iconic uncultured rumen bacterium.

Quinella is a genus of iconic rumen bacteria first reported in 1913. There are no cultures of these bacteria, and information on their physiology is scarce and contradictory. Increased abundance of Quinella was previously found in the rumens of some sheep that emit low amounts of methane (CH4) relative to their feed intake, but whether Quinella contributes to low CH4 emissions is not known. Here, we concentrate Quinella cells from sheep rumen contents, extract and sequence DNA, and reconstruct Quinella genomes that are >90% complete with as little as 0.20% contamination. Bioinformatic analyses of the encoded proteins indicate that lactate and propionate formation are major fermentation pathways. The presence of a gene encoding a potential uptake hydrogenase suggests that Quinella might be able to use free hydrogen (H2). None of the inferred metabolic pathways is predicted to produce H2, a major precursor of CH4, which is consistent with the lower CH4 emissions from those sheep with high abundances of this bacterium.

The large episomes of Butyrivibrio proteoclasticus B316T have arisen through intragenomic gene shuttling from the chromosome to smaller Butyrivibrio-specific plasmids.

The genome of Butyrivibrio proteoclasticus B316(T) contains three large episomes including a 302 kb chromid (BPc2) and two large plasmids of 361 (pCY360) and 186 kb (pCY186). The two plasmids are largely cryptic and it is therefore difficult to gauge their contributions or importance to the biology of B. proteoclasticus. Here, we provide evidence that at least BPc2 and pCY360 are essential as neither could be cured using several previously described curing techniques. We show that BPc2 exists at a copy number of 1, while pCY360 and pCY186 exist at copy numbers of 4 and 0.9, respectively. Yet the transcriptional activities of each episome are much less than that of the 3.5 Mb chromosome. Codon usage analyses did not support the hypothesis that the genes of all three episomes were acquired horizontally. Instead our analyses suggest that the vast majority of genes on each episome were transferred from the 3.5 Mb B. proteoclasticus chromosome. Analysis of their replication origins, however, suggests the plasmid backbones share an evolutionary lineage with the smaller Butyrivibrio specific plasmids, pRJF1 and pRJF2. A survey of 13 species of the Butyrivibrio/Pseudobutyrivibrio assemblage identified similar large episomes in nine strains. DNA hybridization experiments revealed none contained an rRNA operon and only a 145 kb episome from Pseudobutyrivibrioruminis possessed an ortholog of the pCY360 plasmid replication initiation protein. The size and distribution of episomes within the nine strains of Butyrivibrio/Pseudobutyrivibrio showed no correlation with 16S rRNA based phylogeny, leading to a hypothesis that the large episomes of Butyrivibrio spp., have arisen through intragenomic gene transfer events from the chromosome to small horizontally acquired elements.

Presence of novel, potentially homoacetogenic bacteria in the rumen as determined by analysis of formyltetrahydrofolate synthetase sequences from ruminants.

Homoacetogens produce acetate from H(2) and CO(2) via the Wood-Ljungdahl pathway. Some homoacetogens have been isolated from the rumen, but these organisms are expected to be only part of the full diversity present. To survey the presence of rumen homoacetogens, we analyzed sequences of formyltetrahydrofolate synthetase (FTHFS), a key enzyme of the Wood-Ljungdahl pathway. A total of 275 partial sequences of genes encoding FTHFS were PCR amplified from rumen contents of a cow, two sheep, and a deer. Phylogenetic trees were constructed using these FTHFS gene sequences and the translated amino acid sequences, together with other sequences from public databases and from novel nonhomoacetogenic bacteria isolated from the rumen. Over 90% of the FTHFS sequences fell into 34 clusters defined with good bootstrap support. Few rumen-derived FTHFS sequences clustered with sequences of known homoacetogens. Conserved residues were identified in the deduced FTHFS amino acid sequences from known homoacetogens, and their presence in the other sequences was used to determine a "homoacetogen similarity" (HS) score. A homoacetogen FTHFS profile hidden Markov model (HoF-HMM) was used to assess the homology of rumen and homoacetogen FTHFS sequences. Many clusters had low HS scores and HoF-HMM matches, raising doubts about whether the sequences originated from homoacetogens. In keeping with these findings, FTHFS sequences from nonhomoacetogenic bacterial isolates grouped in these clusters with low scores. However, sequences that formed 10 clusters containing no known isolates but representing 15% of our FTHFS sequences from rumen samples had high HS scores and HoF-HMM matches and so could represent novel homoacetogens.

Complete Genome Sequence of the Polysaccharide-Degrading Rumen Bacterium Pseudobutyrivibrio xylanivorans MA3014 Reveals an Incomplete Glycolytic Pathway.

Bacterial species belonging to the genus Pseudobutyrivibrio are important members of the rumen microbiome contributing to the degradation of complex plant polysaccharides. Pseudobutyrivibrio xylanivorans MA3014 was selected for genome sequencing to examine its ability to breakdown and utilize plant polysaccharides. The complete genome sequence of MA3014 is 3.58 Mb, consists of three replicons (a chromosome, chromid, and plasmid), has an overall G + C content of 39.6%, and encodes 3,265 putative protein-coding genes (CDS). Comparative pan-genomic analysis of all cultivated and currently available P. xylanivorans genomes has revealed a strong correlation of orthologous genes within this rumen bacterial species. MA3014 is metabolically versatile and capable of growing on a range of simple mono- or oligosaccharides derived from complex plant polysaccharides such as pectins, mannans, starch, and hemicelluloses, with lactate, butyrate, and formate as the principal fermentation end products. The genes encoding these metabolic pathways have been identified and MA3014 is predicted to encode an extensive range of Carbohydrate-Active enZYmes with 78 glycoside hydrolases, 13 carbohydrate esterases, and 54 glycosyl transferases, suggesting an important role in solubilization of plant matter in the rumen.

Occurrence and expression of genes encoding methyl-compound production in rumen bacteria.

BACKGROUND: Digestive processes in the rumen lead to the release of methyl-compounds, mainly methanol and methylamines, which are used by methyltrophic methanogens to form methane, an important agricultural greenhouse gas. Methylamines are produced from plant phosphatidylcholine degradation, by choline trimethylamine lyase, while methanol comes from demethoxylation of dietary pectins via pectin methylesterase activity. We have screened rumen metagenomic and metatranscriptomic datasets, metagenome assembled genomes, and the Hungate1000 genomes to identify organisms capable of producing methyl-compounds. We also describe the enrichment of pectin-degrading and methane-forming microbes from sheep rumen contents and the analysis of their genomes via metagenomic assembly. RESULTS: Screens of metagenomic data using the protein domains of choline trimethylamine lyase (CutC), and activator protein (CutD) found good matches only to Olsenella umbonata and to Caecibacter, while the Hungate1000 genomes and metagenome assembled genomes from the cattle rumen found bacteria within the phyla Actinobacteria, Firmicutes and Proteobacteria. The cutC and cutD genes clustered with genes that encode structural components of bacterial microcompartment proteins. Prevotella was the dominant genus encoding pectin methyl esterases, with smaller numbers of sequences identified from other fibre-degrading rumen bacteria. Some large pectin methyl esterases (> 2100 aa) were found to be encoded in Butyrivibrio genomes. The pectin-utilising, methane-producing consortium was composed of (i) a putative pectin-degrading bacterium (phylum Tenericutes, class Mollicutes), (ii) a galacturonate-using Sphaerochaeta sp. predicted to produce acetate, lactate, and ethanol, and (iii) a methylotrophic methanogen, Methanosphaera sp., with the ability to form methane via a primary ethanol-dependent, hydrogen-independent, methanogenesis pathway. CONCLUSIONS: The main bacteria that produce methyl-compounds have been identified in ruminants. Their enzymatic activities can now be targeted with the aim of finding ways to reduce the supply of methyl-compound substrates to methanogens, and thereby limit methylotrophic methanogenesis in the rumen.

Improved taxonomic assignment of rumen bacterial 16S rRNA sequences using a revised SILVA taxonomic framework.

The taxonomy and associated nomenclature of many taxa of rumen bacteria are poorly defined within databases of 16S rRNA genes. This lack of resolution results in inadequate definition of microbial community structures, with large parts of the community designated as incertae sedis, unclassified, or uncultured within families, orders, or even classes. We have begun resolving these poorly-defined groups of rumen bacteria, based on our desire to name these for use in microbial community profiling. We used the previously-reported global rumen census (GRC) dataset consisting of >4.5 million partial bacterial 16S rRNA gene sequences amplified from 684 rumen samples and representing a wide range of animal hosts and diets. Representative sequences from the 8,985 largest operational units (groups of sequence sharing >97% sequence similarity, and covering 97.8% of all sequences in the GRC dataset) were used to identify 241 pre-defined clusters (mainly at genus or family level) of abundant rumen bacteria in the ARB SILVA 119 framework. A total of 99 of these clusters (containing 63.8% of all GRC sequences) had no unique or had inadequate taxonomic identifiers, and each was given a unique nomenclature. We assessed this improved framework by comparing taxonomic assignments of bacterial 16S rRNA gene sequence data in the GRC dataset with those made using the original SILVA 119 framework, and three other frameworks. The two SILVA frameworks performed best at assigning sequences to genus-level taxa. The SILVA 119 framework allowed 55.4% of the sequence data to be assigned to 751 uniquely identifiable genus-level groups. The improved framework increased this to 87.1% of all sequences being assigned to one of 871 uniquely identifiable genus-level groups. The new designations were included in the SILVA 123 release (https://www.arb-silva.de/documentation/release-123/) and will be perpetuated in future releases.

Applications of the Soil, Plant and Rumen Microbiomes in Pastoral Agriculture.

The production of dairy, meat, and fiber by ruminant animals relies on the biological processes occurring in soils, forage plants, and the animals' rumens. Each of these components has an associated microbiome, and these have traditionally been viewed as distinct ecosystems. However, these microbiomes operate under similar ecological principles and are connected via water, energy flows, and the carbon and nitrogen nutrient cycles. Here, we summarize the microbiome research that has been done in each of these three environments (soils, forage plants, animals' rumen) and investigate what additional benefits may be possible through understanding the interactions between the various microbiomes. The challenge for future research is to enhance microbiome function by appropriate matching of plant and animal genotypes with the environment to improve the output and environmental sustainability of pastoral agriculture.

Diverse hydrogen production and consumption pathways influence methane production in ruminants.

Farmed ruminants are the largest source of anthropogenic methane emissions globally. The methanogenic archaea responsible for these emissions use molecular hydrogen (H2), produced during bacterial and eukaryotic carbohydrate fermentation, as their primary energy source. In this work, we used comparative genomic, metatranscriptomic and co-culture-based approaches to gain a system-wide understanding of the organisms and pathways responsible for ruminal H2 metabolism. Two-thirds of sequenced rumen bacterial and archaeal genomes encode enzymes that catalyse H2 production or consumption, including 26 distinct hydrogenase subgroups. Metatranscriptomic analysis confirmed that these hydrogenases are differentially expressed in sheep rumen. Electron-bifurcating [FeFe]-hydrogenases from carbohydrate-fermenting Clostridia (e.g., Ruminococcus) accounted for half of all hydrogenase transcripts. Various H2 uptake pathways were also expressed, including methanogenesis (Methanobrevibacter), fumarate and nitrite reduction (Selenomonas), and acetogenesis (Blautia). Whereas methanogenesis-related transcripts predominated in high methane yield sheep, alternative uptake pathways were significantly upregulated in low methane yield sheep. Complementing these findings, we observed significant differential expression and activity of the hydrogenases of the hydrogenogenic cellulose fermenter Ruminococcus albus and the hydrogenotrophic fumarate reducer Wolinella succinogenes in co-culture compared with pure culture. We conclude that H2 metabolism is a more complex and widespread trait among rumen microorganisms than previously recognised. There is evidence that alternative hydrogenotrophs, including acetogenic and respiratory bacteria, can prosper in the rumen and effectively compete with methanogens for H2. These findings may help to inform ongoing strategies to mitigate methane emissions by increasing flux through alternative H2 uptake pathways, including through animal selection, dietary supplementation and methanogenesis inhibitors.

Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection.

Productivity of ruminant livestock depends on the rumen microbiota, which ferment indigestible plant polysaccharides into nutrients used for growth. Understanding the functions carried out by the rumen microbiota is important for reducing greenhouse gas production by ruminants and for developing biofuels from lignocellulose. We present 410 cultured bacteria and archaea, together with their reference genomes, representing every cultivated rumen-associated archaeal and bacterial family. We evaluate polysaccharide degradation, short-chain fatty acid production and methanogenesis pathways, and assign specific taxa to functions. A total of 336 organisms were present in available rumen metagenomic data sets, and 134 were present in human gut microbiome data sets. Comparison with the human microbiome revealed rumen-specific enrichment for genes encoding de novo synthesis of vitamin B12, ongoing evolution by gene loss and potential vertical inheritance of the rumen microbiome based on underrepresentation of markers of environmental stress. We estimate that our Hungate genome resource represents ∼75% of the genus-level bacterial and archaeal taxa present in the rumen.

Dynamics and genetic diversification of Escherichia coli during experimental adaptation to an anaerobic environment.

BACKGROUND: Many bacteria are facultative anaerobes, and can proliferate in both anoxic and oxic environments. Under anaerobic conditions, fermentation is the primary means of energy generation in contrast to respiration. Furthermore, the rates and spectra of spontaneous mutations that arise during anaerobic growth differ to those under aerobic growth. A long-term selection experiment was undertaken to investigate the genetic changes that underpin how the facultative anaerobe, Escherichia coli, adapts to anaerobic environments. METHODS: Twenty-one populations of E. coli REL4536, an aerobically evolved 10,000th generation descendent of the E. coli B strain, REL606, were established from a clonal ancestral culture. These were serially sub-cultured for 2,000 generations in a defined minimal glucose medium in strict aerobic and strict anaerobic environments, as well as in a treatment that fluctuated between the two environments. The competitive fitness of the evolving lineages was assessed at approximately 0, 1,000 and 2,000 generations, in both the environment of selection and the alternative environment. Whole genome re-sequencing was performed on random colonies from all lineages after 2,000-generations. Mutations were identified relative to the ancestral genome, and based on the extent of parallelism, traits that were likely to have contributed towards adaptation were inferred. RESULTS: There were increases in fitness relative to the ancestor among anaerobically evolved lineages when tested in the anaerobic environment, but no increases were found in the aerobic environment. For lineages that had evolved under the fluctuating regime, relative fitness increased significantly in the anaerobic environment, but did not increase in the aerobic environment. The aerobically-evolved lineages did not increase in fitness when tested in either the aerobic or anaerobic environments. The strictly anaerobic lineages adapted more rapidly to the anaerobic environment than did the fluctuating lineages. Two main strategies appeared to predominate during adaptation to the anaerobic environment: modification of energy generation pathways, and inactivation of non-essential functions. Fermentation pathways appeared to alter through selection for mutations in genes such as nadR, adhE, dcuS/R, and pflB. Mutations were frequently identified in genes for presumably dispensable functions such as toxin-antitoxin systems, prophages, virulence and amino acid transport. Adaptation of the fluctuating lineages to the anaerobic environments involved mutations affecting traits similar to those observed in the anaerobically evolved lineages. DISCUSSION: There appeared to be strong selective pressure for activities that conferred cell yield advantages during anaerobic growth, which include restoring activities that had previously been inactivated under long-term continuous aerobic evolution of the ancestor.