We present a programmable droplet-based microfluidic device that combines the reconfigurable flow-routing capabilities of integrated microvalve technology with the sample compartmentalization and dispersion-free transport that is inherent to droplets. The device allows for the execution of user-defined multistep reaction protocols in 95 individually addressable nanoliter-volume storage chambers by consecutively merging programmable sequences of picoliter-volume droplets containing reagents or cells. This functionality is enabled by "flow-controlled wetting," a droplet docking and merging mechanism that exploits the physics of droplet flow through a channel to control the precise location of droplet wetting. The device also allows for automated cross-contaminationfree recovery of reaction products from individual chambers into standard microfuge tubes for downstream analysis. The combined features of programmability, addressability, and selective recovery provide a general hardware platform that can be reprogrammed for multiple applications. We demonstrate this versatility by implementing multiple single-cell experiment types with this device: bacterial cell sorting and cultivation, taxonomic gene identification, and high-throughput single-cell whole genome amplification and sequencing using common laboratory strains. Finally, we apply the device to genome analysis of single cells and microbial consortia from diverse environmental samples including a marine enrichment culture, deep-sea sediments, and the human oral cavity. The resulting datasets capture genotypic properties of individual cells and illuminate known and potentially unique partnerships between microbial community members.two-phase flow | droplet wetting | single-cell analysis | qPCR | environmental genomics M icrofluidic devices provide numerous advantages for biological analysis including automation, enhanced sensitivity and reaction efficiency in small volumes (1, 2), favorable mass transport properties (3, 4), and the potential for scalable and costeffective small volume assays (5). Indeed, advances in microfluidics over the past decade have resulted in increasingly sophisticated functionality and the emergence of two dominant and orthogonal strategies for fluid handling, based either on the use of integrated microvalves or the transport of microdroplets, both in closed channels or over electrode surfaces.The development of soft lithography (6) and the extension of this method to the fabrication of integrated microvalves using multilayer soft lithography (5) has enabled devices with thousands of active microvalves per cm 2 . This high level of integration enables device architectures capable of executing thousands of predefined "unit cell" reactions in parallel, with applications ranging from protein structure (4) and interaction studies (7,8) to single-cell analysis and genomics (2, 9, 10). Two-phase flow systems that manipulate picoliter (pL) volume droplets in closed channels have been shown to be ideally suited to high-speed serial analysis f...

Oil in subsurface reservoirs is biodegraded by resident microbial communities. Water-mediated, anaerobic conversion of hydrocarbons to methane and CO2, catalyzed by syntrophic bacteria and methanogenic archaea, is thought to be one of the dominant processes. We compared 160 microbial community compositions in ten hydrocarbon resource environments (HREs) and sequenced twelve metagenomes to characterize their metabolic potential. Although anaerobic communities were common, cores from oil sands and coal beds had unexpectedly high proportions of aerobic hydrocarbon-degrading bacteria. Likewise, most metagenomes had high proportions of genes for enzymes involved in aerobic hydrocarbon metabolism. Hence, although HREs may have been strictly anaerobic and typically methanogenic for much of their history, this may not hold today for coal beds and for the Alberta oil sands, one of the largest remaining oil reservoirs in the world. This finding may influence strategies to recover energy or chemicals from these HREs by in situ microbial processes.

BackgroundA central challenge to understanding the ecological and biogeochemical roles of microorganisms in natural and human engineered ecosystems is the reconstruction of metabolic interaction networks from environmental sequence information. The dominant paradigm in metabolic reconstruction is to assign functional annotations using BLAST. Functional annotations are then projected onto symbolic representations of metabolism in the form of KEGG pathways or SEED subsystems.ResultsHere we present MetaPathways, an open source pipeline for pathway inference that uses the PathoLogic algorithm to map functional annotations onto the MetaCyc collection of reactions and pathways, and construct environmental Pathway/Genome Databases (ePGDBs) compatible with the editing and navigation features of Pathway Tools. The pipeline accepts assembled or unassembled nucleotide sequences, performs quality assessment and control, predicts and annotates noncoding genes and open reading frames, and produces inputs to PathoLogic. In addition to constructing ePGDBs, MetaPathways uses MLTreeMap to build phylogenetic trees for selected taxonomic anchor and functional gene markers, converts General Feature Format (GFF) files into concatenated GenBank files for ePGDB construction based on third-party annotations, and generates useful file formats including Sequin files for direct GenBank submission and gene feature tables summarizing annotations, MLTreeMap trees, and ePGDB pathway coverage summaries for statistical comparisons.ConclusionsMetaPathways provides users with a modular annotation and analysis pipeline for predicting metabolic interaction networks from environmental sequence information using an alternative to KEGG pathways and SEED subsystems mapping. It is extensible to genomic and transcriptomic datasets from a wide range of sequencing platforms, and generates useful data products for microbial community structure and function analysis. The MetaPathways software package, installation instructions, and example data can be obtained from http://hallam.microbiology.ubc.ca/MetaPathways.

BackgroundOne method for rejuvenating land polluted with anthropogenic contaminants is through phytoremediation, the reclamation of land through the cultivation of specific crops. The capacity for phytoremediation crops, such as Salix spp., to tolerate and even flourish in contaminated soils relies on a highly complex and predominantly cryptic interacting community of microbial life.MethodsHere, Illumina HiSeq 2500 sequencing and de novo transcriptome assembly were used to observe gene expression in washed Salix purpurea cv. ‘Fish Creek’ roots from trees pot grown in petroleum hydrocarbon-contaminated or non-contaminated soil. All 189,849 assembled contigs were annotated without a priori assumption as to sequence origin and differential expression was assessed.ResultsThe 839 contigs differentially expressed (DE) and annotated from S. purpurea revealed substantial increases in transcripts encoding abiotic stress response equipment, such as glutathione S-transferases, in roots of contaminated trees as well as the hallmarks of fungal interaction, such as SWEET2 (Sugars Will Eventually Be Exported Transporter). A total of 8252 DE transcripts were fungal in origin, with contamination conditions resulting in a community shift from Ascomycota to Basidiomycota genera. In response to contamination, 1745 Basidiomycota transcripts increased in abundance (the majority uniquely expressed in contaminated soil) including major monosaccharide transporter MST1, primary cell wall and lamella CAZy enzymes, and an ectomycorrhiza-upregulated exo-β-1,3-glucanase (GH5). Additionally, 639 DE polycistronic transcripts from an uncharacterised Enterobacteriaceae species were uniformly in higher abundance in contamination conditions and comprised a wide spectrum of genes cryptic under laboratory conditions but considered putatively involved in eukaryotic interaction, biofilm formation and dioxygenase hydrocarbon degradation.ConclusionsFungal gene expression, representing the majority of contigs assembled, suggests out-competition of white rot Ascomycota genera (dominated by Pyronema), a sometimes ectomycorrhizal (ECM) Ascomycota (Tuber) and ECM Basidiomycota (Hebeloma) by a poorly characterised putative ECM Basidiomycota due to contamination. Root and fungal expression involved transcripts encoding carbohydrate/amino acid (C/N) dialogue whereas bacterial gene expression included the apparatus necessary for biofilm interaction and direct reduction of contamination stress, a potential bacterial currency for a role in tripartite mutualism. Unmistakable within the metatranscriptome is the degree to which the landscape of rhizospheric biology, particularly the important but predominantly uncharacterised fungal genetics, is yet to be discovered.Electronic supplementary materialThe online version of this article (10.1186/s40168-018-0432-5) contains supplementary material, which is available to authorized users.

Thermocouple arrays were deployed on two deep-sea hydrothermal vents at Guaymas Basin (27 degrees 0.5'N, 111 degrees 24.5'W) in order to measure in situ temperatures at which microorganisms colonize the associated mineral deposits. Intact sections of three structures that formed around the arrays were collected after 4 and 72 day deployments (named BM4, BM72 and TS72). Archaeal diversity associated with discreet subsamples collected across each deposit was determined by polymerase chain reaction amplification of 16S rRNA genes. Spatial differences in archaeal diversity were observed in all deposits and appeared related to in situ temperature. In BM4, no 16S rRNA genes were detected beyond about 1.5 cm within the sample (> 200 degrees C). Phylotypes detected on the outside of this deposit belong to taxonomic groups containing mesophiles and (hyper)thermophiles, whereas only putative hyperthermophiles were detected 1.5 cm inside the structure (approximately 110 degrees C). In contrast, the more moderate thermal gradient recorded across TS72 was associated with a deeper colonization (2-3 cm inside the deposit) of putative hyperthermophilic phylotypes. Although our study does not provide a precise assessment of the highest temperature for the existence of microbial habitats inside the deposits, archaeal 16S rRNA genes were detected directly next to thermocouples that measured 110 degrees C (Methanocaldococcus spp. in BM4) and 116 degrees C (Desulfurococcaceae in TS72). The successive array deployments conducted at the Broken Mushroom (BM) site also revealed compositional differences in archaeal communities associated with immature (BM4) and mature chimneys (BM72) formed by the same fluids. These differences suggest a temporal transition in the primary carbon sources used by the archaeal communities, with potential CO(2)/H(2) methanogens prevalent in BM4 being replaced by possible methylotroph or acetoclastic methanogens and heterotrophs in BM72. This study is the first direct assessment of in situ conditions experienced by microorganisms inhabiting actively forming hydrothermal deposits at different stages of structure development.