Vertical stratification of the air microbiome in the lower troposphere.

The troposphere constitutes the final frontier of global ecosystem research due to technical challenges arising from its size, low biomass, and gaseous state. Using a vertical testing array comprising a meteorological tower and a research aircraft, we conducted synchronized measurements of meteorological parameters and airborne biomass (n = 480) in the vertical air column up to 3,500 m. The taxonomic analysis of metagenomic data revealed differing patterns of airborne microbial community composition with respect to time of day and height above ground. The temporal and spatial resolution of our study demonstrated that the diel cycle of airborne microorganisms is a ground-based phenomenon that is entirely absent at heights >1,000 m. In an integrated analysis combining meteorological and biological data, we demonstrate that atmospheric turbulence, identified by potential temperature and high-frequency three-component wind measurements, is the key driver of bioaerosol dynamics in the lower troposphere. Multivariate regression analysis shows that at least 50% of identified airborne microbial taxa (n = ∼10,000) are associated with either ground or height, allowing for an understanding of dispersal patterns of microbial taxa in the vertical air column. Due to the interconnectedness of atmospheric turbulence and temperature, the dynamics of microbial dispersal are likely to be impacted by rising global temperatures, thereby also affecting ecosystems on the planetary surface.

Microbial communities in the tropical air ecosystem follow a precise diel cycle.

The atmosphere is vastly underexplored as a habitable ecosystem for microbial organisms. In this study, we investigated 795 time-resolved metagenomes from tropical air, generating 2.27 terabases of data. Despite only 9 to 17% of the generated sequence data currently being assignable to taxa, the air harbored a microbial diversity that rivals the complexity of other planetary ecosystems. The airborne microbial organisms followed a clear diel cycle, possibly driven by environmental factors. Interday taxonomic diversity exceeded day-to-day and month-to-month variation. Environmental time series revealed the existence of a large core of microbial taxa that remained invariable over 13 mo, thereby underlining the long-term robustness of the airborne community structure. Unlike terrestrial or aquatic environments, where prokaryotes are prevalent, the tropical airborne biomass was dominated by DNA from eukaryotic phyla. Specific fungal and bacterial species were strongly correlated with temperature, humidity, and CO2 concentration, making them suitable biomarkers for studying the bioaerosol dynamics of the atmosphere.

Metagenomics Reveals a Core Macrolide Resistome Related to Microbiota in Chronic Respiratory Disease.

Rationale: Long-term antibiotic use for managing chronic respiratory disease is increasing; however, the role of the airway resistome and its relationship to host microbiomes remains unknown.Objectives: To evaluate airway resistomes and relate them to host and environmental microbiomes using ultradeep metagenomic shotgun sequencing.Methods: Airway specimens from 85 individuals with and without chronic respiratory disease (severe asthma, chronic obstructive pulmonary disease, and bronchiectasis) were subjected to metagenomic sequencing to an average depth exceeding 20 million reads. Respiratory and device-associated microbiomes were evaluated on the basis of taxonomical classification and functional annotation including the Comprehensive Antibiotic Resistance Database to determine airway resistomes. Co-occurrence networks of gene-microbe association were constructed to determine potential microbial sources of the airway resistome. Paired patient-inhaler metagenomes were compared (n = 31) to assess for the presence of airway-environment overlap in microbiomes and/or resistomes.Measurements and Main Results: Airway metagenomes exhibit taxonomic and metabolic diversity and distinct antimicrobial resistance patterns. A "core" airway resistome dominated by macrolide but with high prevalence of ß-lactam, fluoroquinolone, and tetracycline resistance genes exists and is independent of disease status or antibiotic exposure. Streptococcus and Actinomyces are key potential microbial reservoirs of macrolide resistance including the ermX, ermF, and msrD genes. Significant patient-inhaler overlap in airway microbiomes and their resistomes is identified where the latter may be a proxy for airway microbiome assessment in chronic respiratory disease.Conclusions: Metagenomic analysis of the airway reveals a core macrolide resistome harbored by the host microbiome.

Longitudinal assessment of antibiotic resistance gene profiles in gut microbiomes of infants at risk of eczema.

BACKGROUND: While there is increasing knowledge about the gut microbiome, the factors influencing and the significance of the gut resistome are still not well understood. Infant gut commensals risk transferring multidrug-resistant antibiotic resistance genes (ARGs) to pathogenic bacteria. The rapid spread of multidrug-resistant pathogenic bacteria is a worldwide public health concern. Better understanding of the naïve infant gut resistome may build the evidence base for antimicrobial stewardship in both humans and in the food industry. Given the high carriage rate of extended spectrum beta-lactamase (ESBL)-producing Enterobacteriaceae in Asia, we aimed to evaluate community prevalence, dynamics, and longitudinal changes in antibiotic resistance gene (ARG) profiles and prevalence of ESBL-producing E. coli and K. pneumoniae in the intestinal microbiome of infants participating in the Growing Up in Singapore Towards Healthy Outcomes (GUSTO) study, a longitudinal cohort study of pregnant women and their infants. METHODS: We analysed ARGs in the first year of life among 75 infants at risk of eczema who had stool samples collected at multiple timepoints using metagenomics. RESULTS: The mean number of ARGs per infant increased with age. The most common ARGs identified confer resistance to aminoglycoside, beta-lactam, macrolide and tetracycline antibiotics; all infants harboured these antibiotic resistance genes at some point in the first year of life. Few ARGs persisted throughout the first year of life. Beta-lactam resistant Escherichia coli and Klebsiella pneumoniae were detected in 4 (5.3%) and 32 (42.7%) of subjects respectively. CONCLUSION: In this longitudinal cohort study of infants living in a region with high endemic antibacterial resistance, we demonstrate that majority of the infants harboured several antibiotic resistance genes in their gut and showed that the infant gut resistome is diverse and dynamic over the first year of life.

Long-read sequencing uncovers the adaptive topography of a carnivorous plant genome.

Utricularia gibba, the humped bladderwort, is a carnivorous plant that retains a tiny nuclear genome despite at least two rounds of whole genome duplication (WGD) since common ancestry with grapevine and other species. We used a third-generation genome assembly with several complete chromosomes to reconstruct the two most recent lineage-specific ancestral genomes that led to the modern U. gibba genome structure. Patterns of subgenome dominance in the most recent WGD, both architectural and transcriptional, are suggestive of allopolyploidization, which may have generated genomic novelty and led to instantaneous speciation. Syntenic duplicates retained in polyploid blocks are enriched for transcription factor functions, whereas gene copies derived from ongoing tandem duplication events are enriched in metabolic functions potentially important for a carnivorous plant. Among these are tandem arrays of cysteine protease genes with trap-specific expression that evolved within a protein family known to be useful in the digestion of animal prey. Further enriched functions among tandem duplicates (also with trap-enhanced expression) include peptide transport (intercellular movement of broken-down prey proteins), ATPase activities (bladder-trap acidification and transmembrane nutrient transport), hydrolase and chitinase activities (breakdown of prey polysaccharides), and cell-wall dynamic components possibly associated with active bladder movements. Whereas independently polyploid Arabidopsis syntenic gene duplicates are similarly enriched for transcriptional regulatory activities, Arabidopsis tandems are distinct from those of U. gibba, while still metabolic and likely reflecting unique adaptations of that species. Taken together, these findings highlight the special importance of tandem duplications in the adaptive landscapes of a carnivorous plant genome.

Complete Genome Sequence of Penicillium oxalicum Strain SGAir0226 Isolated from Outdoor Tropical Air in Singapore.

Penicillium oxalicum strain SGAir0226 was isolated from a tropical air sample collected in Singapore. The complete genome was assembled from long reads obtained from single-molecule real-time sequencing and was further polished and error corrected using short read sequencing data. The assembly comprises 20 contigs with a total length of 30.7 Mb. The genome was predicted to contain 8310 protein-coding genes, 237 tRNAs and 83 rRNAs.

Whole-Genome Sequencing of Aspergillus terreus Species Complex.

Aspergillus terreus species complex is an opportunistic fungal pathogen increasingly implicated in invasive infection, as well as chronic respiratory disease. Currently, an understanding of A. terreus pathogenicity is impeded by a limited number of whole-genome sequences of this fungal pathogen. We here describe a high-quality whole-genome assembly of European A. terreus clinical isolate M6925, derived by single-molecule real-time sequencing with short-read polishing.

BatMis: a fast algorithm for k-mismatch mapping.

MOTIVATION: Second-generation sequencing (SGS) generates millions of reads that need to be aligned to a reference genome allowing errors. Although current aligners can efficiently map reads allowing a small number of mismatches, they are not well suited for handling a large number of mismatches. The efficiency of aligners can be improved using various heuristics, but the sensitivity and accuracy of the alignments are sacrificed. In this article, we introduce Basic Alignment tool for Mismatches (BatMis)--an efficient method to align short reads to a reference allowing k mismatches. BatMis is a Burrows-Wheeler transformation based aligner that uses a seed and extend approach, and it is an exact method. RESULTS: Benchmark tests show that BatMis performs better than competing aligners in solving the k-mismatch problem. Furthermore, it can compete favorably even when compared with the heuristic modes of the other aligners. BatMis is a useful alternative for applications where fast k-mismatch mappings, unique mappings or multiple mappings of SGS data are required. AVAILABILITY AND IMPLEMENTATION: BatMis is written in C/C++ and is freely available from http://code.google.com/p/batmis/

The genomic characterisation and comparison of Bacillus cereus strains isolated from indoor air.

BACKGROUND: Bacillus cereus is ubiquitous in nature, found in environments such as soil, plants, air, and part of the insect and human gut microbiome. The ability to produce endospores and biofilms contribute to their pathogenicity, classified in two types of food poisoning: diarrheal and emetic syndromes. Here we report gap-free, whole-genome sequences of two B. cereus strains isolated from air samples and analyse their emetic and diarrheal potential. RESULTS: Genome assemblies of the B. cereus strains consist of one chromosome and seven plasmids each. The genome size of strain SGAir0260 is 6.30-Mb with 6590 predicted coding sequences (CDS) and strain SGAir0263 is 6.47-Mb with 6811 predicted CDS. Macrosynteny analysis showed 99% collinearity between the strains isolated from air and 90.2% with the reference genome. Comparative genomics with 57 complete B. cereus genomes suggests these strains from air are closely associated with strains isolated from foodborne illnesses outbreaks. Due to virulence potential of B. cereus and its reported involvement in nosocomial infections, antibiotic resistance analyses were performed and confirmed resistance to ampicillin and fosfomycin, with susceptibility to ciprofloxacin, tetracycline and vancomycin in both strains. CONCLUSION: Phylogenetic analysis combined with detection of haemolytic (hblA, hblC, and hblD) and non-haemolytic (nheA, nheB, and nheC) enterotoxin genes in both air-isolated strains point to the diarrheic potential of the air isolates, though not emetic. Characterization of these airborne strains and investigation of their potential disease-causing genes could facilitate identification of environmental sources of contamination leading to foodborne illnesses and nosocomial infections transported by air.

Experimental parameters defining ultra-low biomass bioaerosol analysis.

Investigation of the microbial ecology of terrestrial, aquatic and atmospheric ecosystems requires specific sampling and analytical technologies, owing to vastly different biomass densities typically encountered. In particular, the ultra-low biomass nature of air presents an inherent analytical challenge that is confounded by temporal fluctuations in community structure. Our ultra-low biomass pipeline advances the field of bioaerosol research by significantly reducing sampling times from days/weeks/months to minutes/hours, while maintaining the ability to perform species-level identification through direct metagenomic sequencing. The study further addresses all experimental factors contributing to analysis outcome, such as amassment, storage and extraction, as well as factors that impact on nucleic acid analysis. Quantity and quality of nucleic acid extracts from each optimisation step are evaluated using fluorometry, qPCR and sequencing. Both metagenomics and marker gene amplification-based (16S and ITS) sequencing are assessed with regard to their taxonomic resolution and inter-comparability. The pipeline is robust across a wide range of climatic settings, ranging from arctic to desert to tropical environments. Ultimately, the pipeline can be adapted to environmental settings, such as dust and surfaces, which also require ultra-low biomass analytics.

Taxonomic composition and seasonal dynamics of the air microbiome in West Siberia.

Here, we describe taxonomical composition, as well as seasonal and diel dynamics of airborne microbial communities in West Siberia. A total of 78 airborne biomass samples from 39 time intervals were analysed, within a temperature range of 48 °C (26 °C to - 22 °C). We observed a 5-170-fold decrease in DNA yield extracted from the airborne biomass in winter compared to summer, nevertheless, yielding sufficient material for metagenomic analysis. The airborne microbial communities included Actinobacteria and Proteobacteria, Ascomycota and Basidiomycota fungi as major components, as well as some Streptophyta plants. In summer, bacterial and fungal plant pathogens, and wood-rotting saprophytes were predominant. In winter, Ascomycota moulds and cold-related or stress environment bacterial species were enriched, while the fraction of wood-rotting and mushroom-forming Basidiomycota fungi was largely reduced. As recently reported for the tropical climate, the airborne microbial communities performed a diel cycle in summer, however, in winter diel dynamics were not observed.

Complete genome of Enterobacter sichuanensis strain SGAir0282 isolated from air in Singapore.

BACKGROUND: Enterobacter cloacae complex (ECC) bacteria, such as E. cloacae, E. sichuanensis, E. kobei, and E. roggenkampii, have been emerging as nosocomial pathogens. Many strains isolated from medical clinics were found to be resistant to antibiotics, and in the worst cases, acquired multidrug resistance. We present the whole genome sequence of SGAir0282, isolated from the outdoor air in Singapore, and its relevance to other ECC bacteria by in silico genomic analysis. RESULTS: Complete genome assembly of E. sichuanensis strain SGAir0282 was generated using PacBio RSII and Illumina MiSeq platforms, and the datasets were used for de novo assembly using Hierarchical Genome Assembly Process (HGAP) and error corrected with Pilon. The genome assembly consisted of a single contig of 4.71 Mb and with a G+C content of 55.5%. No plasmid was detected in the assembly. The genome contained 4371 coding genes, 83 tRNA and 25 rRNA genes, as predicted by NCBI's Prokaryotic Genome Annotation Pipeline (PGAP). Among the genes, the antibiotic resistance related genes were included: Streptothricin acetdyltransferase (SatA), fosfomycin resistance protein (FosA) and metal-dependent hydrolases of the beta-lactamase superfamily I (BLI). CONCLUSION: Based on whole genome alignment and phylogenetic analysis, the strain SGAir0282 was identified to be Enterobacter sichuanensis. The strain possesses gene clusters for virulence, disease and defence, that can also be found in other multidrug resistant ECC type strains.

A compromised developmental trajectory of the infant gut microbiome and metabolome in atopic eczema.

Evidence is accumulating that the establishment of the gut microbiome in early life influences the development of atopic eczema. In this longitudinal study, we used integrated multi-omics analyses to infer functional mechanisms by which the microbiome modulates atopic eczema risk. We measured the functionality of the gut microbiome and metabolome of 63 infants between ages 3 weeks and 12 months with well-defined eczema cases and controls in a sub-cohort from the Growing Up in Singapore Toward healthy Outcomes (GUSTO) mother-offspring cohort. At 3 weeks, the microbiome and metabolome of allergen-sensitized atopic eczema infants were characterized by an enrichment of Escherichia coli and Klebsiella pneumoniae, associated with increased stool D-glucose concentration and increased gene expression of associated virulence factors. A delayed colonization by beneficial Bacteroides fragilis and subsequent delayed accumulation of butyrate and propionate producers after 3 months was also observed. Here, we describe an aberrant developmental trajectory of the gut microbiome and stool metabolome in allergen sensitized atopic eczema infants. The infographic describes an impaired developmental trajectory of the gut microbiome and metabolome in allergen-sensitized atopic eczema (AE) infants and infer its contribution in modulating allergy risk in the Singaporean mother-offspring GUSTO cohort. The key microbial signature of AE is characterized by (1) an enrichment of Escherichia coli and Klebsiella pneumoniae which are associated with accumulation of pre-glycolysis intermediates (D-glucose) via the trehalose metabolic pathway, increased gene expression of associated virulence factors (invasin, adhesin, flagellin and lipopolysaccharides) by utilizing ATP from oxidative phosphorylation and delayed production of butyrate and propionate, (2) depletion of Bacteroides fragilis which resulted in lower expression of immunostimulatory bacterial cell envelope structure and folate (vitamin B9) biosynthesis pathway, and (3) accompanied depletion of bacterial groups with the ability to derive butyrate and propionate through direct or indirect pathways which collectively resulted in reduced glycolysis, butyrate and propionate biosynthesis.

Complete Genome Sequence of Brevundimonas sp. Strain SGAir0440, Isolated from Indoor Air in Singapore.

Brevundimonas sp. strain SGAir0440 was isolated from indoor air samples collected in Singapore. Its genome was assembled using single-molecule real-time sequencing data, resulting in one circular chromosome with a length of 3.1 Mbp. The genome consists of 3,033 protein-coding genes, 48 tRNAs, and 6 rRNA operons.

Complete Genome Sequence of Streptomyces sp. Strain SGAir0924, an Actinobacterium Isolated from Outdoor Air in Singapore.

Streptomyces sp. strain SGAir0924 was isolated from outdoor air collected in Singapore. Its genome was assembled using long reads generated by single-molecule real-time sequencing. The final assembly had one chromosome of 7.65 Mb and three plasmids with an average length of 142 kb. The genome contained 6,825 protein-coding genes, 68 tRNAs, and 18 rRNAs.

Complete Genome Sequence of Enterococcus faecalis Strain SGAir0397, Isolated from a Tropical Air Sample Collected in Singapore.

Enterococcus faecalis strain SGAir0397 was isolated from a tropical air sample collected in Singapore. Its genome was assembled using single-molecule real-time sequencing data and comprises one circular chromosome with a length of 2.69 Mbp. The genome contains 2,595 protein-coding genes, 59 tRNAs, and 12 rRNAs.

Complete Genome Sequence of Citricoccus sp. Strain SGAir0253, Isolated from Indoor Air in Singapore.

Citricoccus sp. strain SGAir0253 was isolated from indoor air collected in Singapore. Its genome sequence was assembled using single-molecule real-time sequencing. It comprises one chromosome of 3.32 Mb and two plasmids of 137 kb and 99 kb. The genome consists of 2,950 protein-coding genes, 49 tRNAs, and 9 rRNAs.

Complete Genome Sequence of Micrococcus luteus Strain SGAir0127, Isolated from Indoor Air Samples from Singapore.

Micrococcus luteus strain SGAir0127 was isolated from indoor air samples collected in Singapore. The assembly, based on single-molecule real-time sequencing reads, resulted in two contigs, one chromosomal contig with a length of 2.57 Mbp and one nonchromosomal contig of 8.68 kbp. The genome has a total of 2,564 genes.

Complete Genome Sequence of Rhodococcus sp. Strain SGAir0479, Isolated from Indoor Air Collected in Singapore.

The complete genome sequence of Rhodococcus sp. strain SGAir0479 is presented here. This organism was isolated from an air sample collected in an indoor location in Singapore. The consensus assembly generated one chromosome of 4.86 Mb (G+C content of 69.8%) and one plasmid of 104,493 bp.

Complete Genome Sequence of Nissabacter sp. Strain SGAir0207, Isolated from an Air Sample Collected in Singapore.

Nissabacter sp. strain SGAir0207 was isolated from a tropical air sample collected in Singapore. Its genome was assembled using a hybrid approach with long and short reads, resulting in one chromosome of 3.9 Mb and 7 plasmids. The complete genome consists of 4,403 protein-coding, 84 tRNA, and 22 rRNA genes.