Current concepts, advances, and challenges in deciphering the human microbiota with metatranscriptomics.

Metatranscriptomics refers to the analysis of the collective microbial transcriptome of a sample. Its increased utilization for the characterization of human-associated microbial communities has enabled the discovery of many disease-state related microbial activities. Here, we review the principles of metatranscriptomics-based analysis of human-associated microbial samples. We describe strengths and weaknesses of popular sample preparation, sequencing, and bioinformatics approaches and summarize strategies for their use. We then discuss how human-associated microbial communities have recently been examined and how their characterization may change. We conclude that metatranscriptomics insights into human microbiotas under health and disease have not only expanded our knowledge on human health, but also opened avenues for rational antimicrobial drug use and disease management.

Comparative genomics of Lactobacillus crispatus suggests novel mechanisms for the competitive exclusion of Gardnerella vaginalis.

BACKGROUND: Lactobacillus crispatus is a ubiquitous micro-organism encountered in a wide range of host-associated habitats. It can be recovered from the gastrointestinal tract of animals and it is a common constituent of the vaginal microbiota of humans. Moreover, L. crispatus can contribute to the urogenital health of the host through competitive exclusion and the production of antimicrobial agents. In order to investigate the genetic diversity of this important urogenital species, we performed a comparative genomic analysis of L. crispatus. RESULTS: Utilizing the completed genome sequence of a strain ST1 and the draft genome sequences of nine other L. crispatus isolates, we defined the scale and scope of the pan- and core genomic potential of L. crispatus. Our comparative analysis identified 1,224 and 2,705 ortholog groups present in all or only some of the ten strains, respectively. Based on mathematical modeling, sequencing of additional L. crispatus isolates would result in the identification of new genes and functions, whereas the conserved core of the ten strains was a good representation of the final L. crispatus core genome, estimated to level at about 1,116 ortholog groups. Importantly, the current core was observed to encode bacterial components potentially promoting urogenital health. Using antibody fragments specific for one of the conserved L. crispatus adhesins, we demonstrated that the L. crispatus core proteins have a potential to reduce the ability of Gardnerella vaginalis to adhere to epithelial cells. These findings thereby suggest that L. crispatus core proteins could protect the vagina from G. vaginalis and bacterial vaginosis. CONCLUSIONS: Our pan-genome analysis provides insights into the intraspecific genome variability and the collective molecular mechanisms of the species L. crispatus. Using this approach, we described the differences and similarities between the genomes and identified features likely to be important for urogenital health. Notably, the conserved genetic backbone of L. crispatus accounted for close to 60% of the ortholog groups of an average L. crispatus strain and included factors for the competitive exclusion of G. vaginalis, providing an explanation on how this urogenital species could improve vaginal health.

Understanding human health through metatranscriptomics.

Metatranscriptomics has revolutionized our ability to explore and understand transcriptional programs in microbial communities. Moreover, it has enabled us to gain deeper and more specific insight into the microbial activities in human gut, respiratory, oral, and vaginal communities. Perhaps the most important contribution of metatranscriptomics arises, however, from the analyses of disease-associated communities. We review the advantages and disadvantages of metatranscriptomics analyses in understanding human health and disease. We focus on human tissues low in microbial biomass and conditions associated with dysbiotic microbiota. We conclude that a more widespread use of metatranscriptomics and increased knowledge on microbe activities will uncover critical interactions between microbes and host in human health and provide diagnostic basis for culturing-independent, direct functional pathogen identification.

BLANNOTATOR: enhanced homology-based function prediction of bacterial proteins.

BACKGROUND: Automated function prediction has played a central role in determining the biological functions of bacterial proteins. Typically, protein function annotation relies on homology, and function is inferred from other proteins with similar sequences. This approach has become popular in bacterial genomics because it is one of the few methods that is practical for large datasets and because it does not require additional functional genomics experiments. However, the existing solutions produce erroneous predictions in many cases, especially when query sequences have low levels of identity with the annotated source protein. This problem has created a pressing need for improvements in homology-based annotation. RESULTS: We present an automated method for the functional annotation of bacterial protein sequences. Based on sequence similarity searches, BLANNOTATOR accurately annotates query sequences with one-line summary descriptions of protein function. It groups sequences identified by BLAST into subsets according to their annotation and bases its prediction on a set of sequences with consistent functional information. We show the results of BLANNOTATOR's performance in sets of bacterial proteins with known functions. We simulated the annotation process for 3090 SWISS-PROT proteins using a database in its state preceding the functional characterisation of the query protein. For this dataset, our method outperformed the five others that we tested, and the improved performance was maintained even in the absence of highly related sequence hits. We further demonstrate the value of our tool by analysing the putative proteome of Lactobacillus crispatus strain ST1. CONCLUSIONS: BLANNOTATOR is an accurate method for bacterial protein function prediction. It is practical for genome-scale data and does not require pre-existing sequence clustering; thus, this method suits the needs of bacterial genome and metagenome researchers. The method and a web-server are available at http://ekhidna.biocenter.helsinki.fi/poxo/blannotator/.

Ischemic Heart Disease Selectively Modifies the Right Atrial Appendage Transcriptome.

Background: Although many pathological changes have been associated with ischemic heart disease (IHD), molecular-level alterations specific to the ischemic myocardium and their potential to reflect disease severity or therapeutic outcome remain unclear. Currently, diagnosis occurs relatively late and evaluating disease severity is largely based on clinical symptoms, various imaging modalities, or the determination of risk factors. This study aims to identify IHD-associated signature RNAs from the atrial myocardium and evaluate their ability to reflect disease severity or cardiac surgery outcomes. Methods and Results: We collected right atrial appendage (RAA) biopsies from 40 patients with invasive coronary angiography (ICA)-positive IHD undergoing coronary artery bypass surgery and from 8 patients ICA-negative for IHD (non-IHD) undergoing valvular surgery. Following RNA sequencing, RAA transcriptomes were analyzed against 429 donors from the GTEx project without cardiac disease. The IHD transcriptome was characterized by repressed RNA expression in pathways for cell-cell contacts and mitochondrial dysfunction. Increased expressions of the CSRNP3, FUT10, SHD, NAV2-AS4, and hsa-mir-181 genes resulted in significance with the complexity of coronary artery obstructions or correlated with a functional cardiac benefit from bypass surgery. Conclusions: Our results provide an atrial myocardium-focused insight into IHD signature RNAs. The specific gene expression changes characterized here, pave the way for future disease mechanism-based identification of biomarkers for early detection and treatment of IHD.

Genome sequence of Lactobacillus crispatus ST1.

Lactobacillus crispatus is a common member of the beneficial microbiota present in the vertebrate gastrointestinal and human genitourinary tracts. Here, we report the genome sequence of L. crispatus ST1, a chicken isolate displaying strong adherence to vaginal epithelial cells.

Mutation accumulation in cancer genes relates to nonoptimal outcome in chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm accounting for ∼15% of all leukemia. Progress of the disease from an indolent chronic phase to the more aggressive accelerated phase or blast phase (BP) occurs in a minority of cases and is associated with an accumulation of somatic mutations. We performed genetic profiling of 85 samples and transcriptome profiling of 12 samples from 59 CML patients. We identified recurrent somatic mutations in ABL1 (37%), ASXL1 (26%), RUNX1 (16%), and BCOR (16%) in the BP and observed that mutation signatures in the BP resembled those of acute myeloid leukemia (AML). We found that mutation load differed between the indolent and aggressive phases and that nonoptimal responders had more nonsilent mutations than did optimal responders at the time of diagnosis, as well as in follow-up. Using RNA sequencing, we identified other than BCR-ABL1 cancer-associated hybrid genes in 6 of the 7 BP samples. Uncovered expression alterations were in turn associated with mechanisms and pathways that could be targeted in CML management and by which somatic alterations may emerge in CML. Last, we showed the value of genetic data in CML management in a personalized medicine setting.

Aggressive natural killer-cell leukemia mutational landscape and drug profiling highlight JAK-STAT signaling as therapeutic target.

Aggressive natural killer-cell (NK-cell) leukemia (ANKL) is an extremely aggressive malignancy with dismal prognosis and lack of targeted therapies. Here, we elucidate the molecular pathogenesis of ANKL using a combination of genomic and drug sensitivity profiling. We study 14 ANKL patients using whole-exome sequencing (WES) and identify mutations in STAT3 (21%) and RAS-MAPK pathway genes (21%) as well as in DDX3X (29%) and epigenetic modifiers (50%). Additional alterations include JAK-STAT copy gains and tyrosine phosphatase mutations, which we show recurrent also in extranodal NK/T-cell lymphoma, nasal type (NKTCL) through integration of public genomic data. Drug sensitivity profiling further demonstrates the role of the JAK-STAT pathway in the pathogenesis of NK-cell malignancies, identifying NK cells to be highly sensitive to JAK and BCL2 inhibition compared to other hematopoietic cell lineages. Our results provide insight into ANKL genetics and a framework for application of targeted therapies in NK-cell malignancies.

Functional genomics provides insights into the role of Propionibacterium freudenreichii ssp. shermanii JS in cheese ripening.

Propionibacterium freudenreichii is a commercially important bacterium that is essential for the development of the characteristic eyes and flavor of Swiss-type cheeses. These bacteria grow actively and produce large quantities of flavor compounds during cheese ripening at warm temperatures but also appear to contribute to the aroma development during the subsequent cold storage of cheese. Here, we advance our understanding of the role of P. freudenreichii in cheese ripening by presenting the 2.68-Mbp annotated genome sequence of P. freudenreichii ssp. shermanii JS and determining its global transcriptional profiles during industrial cheese-making using transcriptome sequencing. The annotation of the genome identified a total of 2377 protein-coding genes and revealed the presence of enzymes and pathways for formation of several flavor compounds. Based on transcriptome profiling, the expression of 348 protein-coding genes was altered between the warm and cold room ripening of cheese. Several propionate, acetate, and diacetyl/acetoin production related genes had higher expression levels in the warm room, whereas a general slowing down of the metabolism and an activation of mobile genetic elements was seen in the cold room. A few ripening-related and amino acid catabolism involved genes were induced or remained active in cold room, indicating that strain JS contributes to the aroma development also during cold room ripening. In addition, we performed a comparative genomic analysis of strain JS and 29 other Propionibacterium strains of 10 different species, including an isolate of both P. freudenreichii subspecies freudenreichii and shermanii. Ortholog grouping of the predicted protein sequences revealed that close to 86% of the ortholog groups of strain JS, including a variety of ripening-related ortholog groups, were conserved across the P. freudenreichii isolates. Taken together, this study contributes to the understanding of the genomic basis of P. freudenreichii and sheds light on its activities during cheese ripening.

Genome-Wide Analysis of Evolutionary Markers of Human Influenza A(H1N1)pdm09 and A(H3N2) Viruses May Guide Selection of Vaccine Strain Candidates.

Here we analyzed whole-genome sequences of 3,969 influenza A(H1N1)pdm09 and 4,774 A(H3N2) strains that circulated during 2009-2015 in the world. The analysis revealed changes at 481 and 533 amino acid sites in proteins of influenza A(H1N1)pdm09 and A(H3N2) strains, respectively. Many of these changes were introduced as a result of random drift. However, there were 61 and 68 changes that were present in relatively large number of A(H1N1)pdm09 and A(H3N2) strains, respectively, that circulated during relatively long time. We named these amino acid substitutions evolutionary markers, as they seemed to contain valuable information regarding the viral evolution. Interestingly, influenza A(H1N1)pdm09 and A(H3N2) viruses acquired non-overlapping sets of evolutionary markers. We next analyzed these characteristic markers in vaccine strains recommended by the World Health Organization for the past five years. Our analysis revealed that vaccine strains carried only few evolutionary markers at antigenic sites of viral hemagglutinin (HA) and neuraminidase (NA). The absence of these markers at antigenic sites could affect the recognition of HA and NA by human antibodies generated in response to vaccinations. This could, in part, explain moderate efficacy of influenza vaccines during 2009-2014. Finally, we identified influenza A(H1N1)pdm09 and A(H3N2) strains, which contain all the evolutionary markers of influenza A strains circulated in 2015, and which could be used as vaccine candidates for the 2015/2016 season. Thus, genome-wide analysis of evolutionary markers of influenza A(H1N1)pdm09 and A(H3N2) viruses may guide selection of vaccine strain candidates.

Comparative Analysis of Whole-Genome Sequences of Influenza A(H1N1)pdm09 Viruses Isolated from Hospitalized and Nonhospitalized Patients Identifies Missense Mutations That Might Be Associated with Patient Hospital Admissions in Finland during 2009 to 2014.

Here, we report 40 new whole-genome sequences of influenza A(H1N1)pdm09 viruses isolated from Finnish patients during 2009 to 2014. A preliminary analysis of these and 186 other whole genomes of influenza A(H1N1)pdm09 viruses isolated from hospitalized and nonhospitalized patients during 2009 to 2014 in Finland revealed several viral mutations that might be associated with patient hospitalizations.

Invited review: Breaking barriers--attack on innate immune defences by omptin surface proteases of enterobacterial pathogens.

The omptin family of Gram-negative bacterial transmembrane aspartic proteases comprises surface proteins with a highly conserved beta-barrel fold but differing biological functions. The omptins OmpT of Escherichia coli, PgtE of Salmonella enterica, and Pla of Yersinia pestis differ in their substrate specificity as well as in control of their expression. Their functional differences are in accordance with the differing pathogenesis of the infections caused by E. coli, Salmonella, and Y. pestis, which suggests that the omptins have adapted to the life-styles of their host species. The omptins Pla and PgtE attack on innate immunity by affecting the plasminogen/plasmin, complement, coagulation, fibrinolysis, and matrix metalloproteinase systems, by inactivating antimicrobial peptides, and by enhancing bacterial adhesiveness and invasiveness. Although the mechanistic details of the functions of Pla and PgtE differ, the outcome is the same: enhanced spread and multiplication of Y. pestis and S. enterica in the host. The omptin OmpT is basically a housekeeping protease but it also degrades cationic antimicrobial peptides and may enhance colonization of E. coli at uroepithelia. The catalytic residues in the omptin molecules are spatially conserved, and the differing polypeptide substrate specificities are dictated by minor sequence variations at regions surrounding the catalytic cleft. For enzymatic activity, omptins require association with lipopolysaccharide on the outer membrane. Modification of lipopolysaccharide by in vivo conditions or by bacterial gene loss has an impact on omptin function. Creation of bacterial surface proteolysis is thus a coordinated function involving several surface structures.