Species-level classification for 16S rRNA gene sequences remains a serious challenge for microbiome researchers, because existing taxonomic classification tools for 16S rRNA gene sequences either do not provide species-level classification, or their classification results are unreliable. The unreliable results are due to the limitations in the existing methods which either lack solid probabilistic-based criteria to evaluate the confidence of their taxonomic assignments, or use nucleotide k-mer frequency as the proxy for sequence similarity measurement. We have developed a method that shows significantly improved species-level classification results over existing methods. Our method calculates true sequence similarity between query sequences and database hits using pairwise sequence alignment. Taxonomic classifications are assigned from the species to the phylum levels based on the lowest common ancestors of multiple database hits for each query sequence, and further classification reliabilities are evaluated by bootstrap confidence scores. The novelty of our method is that the contribution of each database hit to the taxonomic assignment of the query sequence is weighted by a Bayesian posterior probability based upon the degree of sequence similarity of the database hit to the query sequence. Our method does not need any training datasets specific for different taxonomic groups. Instead only a reference database is required for aligning to the query sequences, making our method easily applicable for different regions of the 16S rRNA gene or other phylogenetic marker genes. Reliable species-level classification for 16S rRNA or other phylogenetic marker genes is critical for microbiome research. Our software shows significantly higher classification accuracy than the existing tools and we provide probabilistic-based confidence scores to evaluate the reliability of our taxonomic classification assignments based on multiple database matches to query sequences. Despite its higher computational costs, our method is still suitable for analyzing large-scale microbiome datasets for practical purposes. Furthermore, our method can be applied for taxonomic classification of any phylogenetic marker gene sequences. Our software, called BLCA, is freely available at https://github.com/qunfengdong/BLCA .

Domestic combustion of biomass fuels, such as wood, charcoal, crop residue and dung causes Household Air Pollution (HAP). These inhaled particulates affect more than half of the world’s population, causing respiratory problems such as infection and inflammatory lung disease. We examined whether the presence of black carbon in alveolar macrophages was associated with alterations in the lung microbiome in a Malawi population. Bronchoalveolar lavage samples from 44 healthy adults were sequenced using 16S rDNA amplification to assess microbial diversity, richness and relative taxa abundance. Individuals were classified as high or low particulate exposure as determined by questionnaire and the percentage of black carbon within their alveolar macrophages. Subjects in the low and high particulate groups did not differ in terms of source of fuels used for cooking or lighting. There was no difference in alpha or beta diversity by particulate group. Neisseria and Streptococcus were significantly more abundant in samples from high particulate exposed individuals, and Tropheryma was found less abundant. Petrobacter abundance was higher in people using biomass fuel for household cooking and lighting, compared with exclusive use of electricity. Healthy adults in Malawi exposed to higher levels of particulates have higher abundances of potentially pathogenic bacteria (Streptococcus, Neisseria) within their lung microbiome. Domestic biomass fuel use was associated with an uncommon environmental bacterium (Petrobacter) associated with oil-rich niches.