Alterations of oral microbiota and impact on the gut microbiome in type 1 diabetes mellitus revealed by integrated multi-omic analyses.

BACKGROUND: Alterations to the gut microbiome have been linked to multiple chronic diseases. However, the drivers of such changes remain largely unknown. The oral cavity acts as a major route of exposure to exogenous factors including pathogens, and processes therein may affect the communities in the subsequent compartments of the gastrointestinal tract. Here, we perform strain-resolved, integrated meta-genomic, transcriptomic, and proteomic analyses of paired saliva and stool samples collected from 35 individuals from eight families with multiple cases of type 1 diabetes mellitus (T1DM). RESULTS: We identified distinct oral microbiota mostly reflecting competition between streptococcal species. More specifically, we found a decreased abundance of the commensal Streptococcus salivarius in the oral cavity of T1DM individuals, which is linked to its apparent competition with the pathobiont Streptococcus mutans. The decrease in S. salivarius in the oral cavity was also associated with its decrease in the gut as well as higher abundances in facultative anaerobes including Enterobacteria. In addition, we found evidence of gut inflammation in T1DM as reflected in the expression profiles of the Enterobacteria as well as in the human gut proteome. Finally, we were able to follow transmitted strain-variants from the oral cavity to the gut at the individual omic levels, highlighting not only the transfer, but also the activity of the transmitted taxa along the gastrointestinal tract. CONCLUSIONS: Alterations of the oral microbiome in the context of T1DM impact the microbial communities in the lower gut, in particular through the reduction of "mouth-to-gut" transfer of Streptococcus salivarius. Our results indicate that the observed oral-cavity-driven gut microbiome changes may contribute towards the inflammatory processes involved in T1DM. Through the integration of multi-omic analyses, we resolve strain-variant "mouth-to-gut" transfer in a disease context. Video Abstract.

[The gut microbiome in Parkinson's disease]. / Das Darmmikrobiom bei der Parkinson-Krankheit.

The vast majority of Parkinson's disease (PD) cases are of sporadic origin and despite extensive research in recent years, the etiology still remains unclear. Several current case control studies are aiming to characterize a putative PD-specific composition of the gut microbiome, reflecting the potential relevance of microbiota in the pathogenesis of PD. Although methodologies and cohort sizes differed, the currently available studies showed reproducible or consistent results in terms of PD-specific alterations to the intestinal bacteria. By applying metagenomic sequencing procedures, it is even possible to distinguish PD cases from healthy individuals at a very early disease stage by means of individually modified microbiota. Among others, microbiota that are associated with an altered intestinal barrier or immune function, such as Akkermansia, Lactobacillus, Faecalibacterium and Prevotella were significantly over-represented or under-represented. There may even be a prodromal microbiome, as a comparable microbial shift is also found in patients with rapid eye movement (REM) sleep behavior disorder (RBD), a risk factor for the later development of synucleinopathies, such as PD.

The gut microbiome is associated with behavioural task in honey bees.

The gut microbiome is recognised as playing an integral role in the health and ecology of a wide variety of animal taxa. However, the relationship between social behavioural traits and the microbial community has received little attention. Honey bees are highly social and the workers perform different behavioural tasks in the colony that cause them to be exposed to different local environments. Here we examined whether the gut microbial community composition of worker honey bees is associated with the behavioural tasks they perform, and therefore also the local environment they are exposed to. We set up five observation hives, in which all workers were matched in age and observed the behaviour of marked bees in each colony over 4 days. The gut bacterial communities of bees seen performing predominantly foraging or predominantly in nest tasks were then characterised and compared based on amplicon sequencing of the 16S rRNA gene. Our results show that some core members of the unique honey bee gut bacterial community are represented in different relative abundances in bees performing different behavioural tasks. The differentially represented bacterial taxa include some thought to be important in carbohydrate metabolism and transport, and also linked to bee health. The results suggest an influence of task-related local environment exposure and diet on the honey bee gut microbial community and identify focal core taxa for further functional analyses.

Functional implications of microbial and viral gut metagenome changes in early stage L-DOPA-naïve Parkinson's disease patients.

BACKGROUND: Parkinson's disease (PD) presently is conceptualized as a protein aggregation disease in which pathology involves both the enteric and the central nervous system, possibly spreading from one to another via the vagus nerves. As gastrointestinal dysfunction often precedes or parallels motor symptoms, the enteric system with its vast diversity of microorganisms may be involved in PD pathogenesis. Alterations in the enteric microbial taxonomic level of L-DOPA-naïve PD patients might also serve as a biomarker. METHODS: We performed metagenomic shotgun analyses and compared the fecal microbiomes of 31 early stage, L-DOPA-naïve PD patients to 28 age-matched controls. RESULTS: We found increased Verrucomicrobiaceae (Akkermansia muciniphila) and unclassified Firmicutes, whereas Prevotellaceae (Prevotella copri) and Erysipelotrichaceae (Eubacterium biforme) were markedly lowered in PD samples. The observed differences could reliably separate PD from control with a ROC-AUC of 0.84. Functional analyses of the metagenomes revealed differences in microbiota metabolism in PD involving the ẞ-glucuronate and tryptophan metabolism. While the abundances of prophages and plasmids did not differ between PD and controls, total virus abundance was decreased in PD participants. Based on our analyses, the intake of either a MAO inhibitor, amantadine, or a dopamine agonist (which in summary relates to 90% of PD patients) had no overall influence on taxa abundance or microbial functions. CONCLUSIONS: Our data revealed differences of colonic microbiota and of microbiota metabolism between PD patients and controls at an unprecedented detail not achievable through 16S sequencing. The findings point to a yet unappreciated aspect of PD, possibly involving the intestinal barrier function and immune function in PD patients. The influence of the parkinsonian medication should be further investigated in the future in larger cohorts.

eggNOG v2.0: extending the evolutionary genealogy of genes with enhanced non-supervised orthologous groups, species and functional annotations.

The identification of orthologous relationships forms the basis for most comparative genomics studies. Here, we present the second version of the eggNOG database, which contains orthologous groups (OGs) constructed through identification of reciprocal best BLAST matches and triangular linkage clustering. We applied this procedure to 630 complete genomes (529 bacteria, 46 archaea and 55 eukaryotes), which is a 2-fold increase relative to the previous version. The pipeline yielded 224,847 OGs, including 9724 extended versions of the original COG and KOG. We computed OGs for different levels of the tree of life; in addition to the species groups included in our first release (i.e. fungi, metazoa, insects, vertebrates and mammals), we have now constructed OGs for archaea, fishes, rodents and primates. We automatically annotate the non-supervised orthologous groups (NOGs) with functional descriptions, protein domains, and functional categories as defined initially for the COG/KOG database. In-depth analysis is facilitated by precomputed high-quality multiple sequence alignments and maximum-likelihood trees for each of the available OGs. Altogether, eggNOG covers 2,242 035 proteins (built from 2,590,259 proteins) and provides a broad functional description for at least 1,966,709 (88%) of them. Users can access the complete set of orthologous groups via a web interface at: http://eggnog.embl.de.

Quantitative assessment of protein function prediction from metagenomics shotgun sequences.

To assess the potential of protein function prediction in environmental genomics data, we analyzed shotgun sequences from four diverse and complex habitats. Using homology searches as well as customized gene neighborhood methods that incorporate intergenic and evolutionary distances, we inferred specific functions for 76% of the 1.4 million predicted ORFs in these samples (83% when nonspecific functions are considered). Surprisingly, these fractions are only slightly smaller than the corresponding ones in completely sequenced genomes (83% and 86%, respectively, by using the same methodology) and considerably higher than previously thought. For as many as 75,448 ORFs (5% of the total), only neighborhood methods can assign functions, illustrated here by a previously undescribed gene associated with the well characterized heme biosynthesis operon and a potential transcription factor that might regulate a coupling between fatty acid biosynthesis and degradation. Our results further suggest that, although functions can be inferred for most proteins on earth, many functions remain to be discovered in numerous small, rare protein families.

Quantitative phylogenetic assessment of microbial communities in diverse environments.

The taxonomic composition of environmental communities is an important indicator of their ecology and function. We used a set of protein-coding marker genes, extracted from large-scale environmental shotgun sequencing data, to provide a more direct, quantitative, and accurate picture of community composition than that provided by traditional ribosomal RNA-based approaches depending on the polymerase chain reaction. Mapping marker genes from four diverse environmental data sets onto a reference species phylogeny shows that certain communities evolve faster than others. The method also enables determination of preferred habitats for entire microbial clades and provides evidence that such habitat preferences are often remarkably stable over time.

The Helmholtz Network for Bioinformatics: an integrative web portal for bioinformatics resources.

SUMMARY: The Helmholtz Network for Bioinformatics (HNB) is a joint venture of eleven German bioinformatics research groups that offers convenient access to numerous bioinformatics resources through a single web portal. The 'Guided Solution Finder' which is available through the HNB portal helps users to locate the appropriate resources to answer their queries by employing a detailed, tree-like questionnaire. Furthermore, automated complex tool cascades ('tasks'), involving resources located on different servers, have been implemented, allowing users to perform comprehensive data analyses without the requirement of further manual intervention for data transfer and re-formatting. Currently, automated cascades for the analysis of regulatory DNA segments as well as for the prediction of protein functional properties are provided. AVAILABILITY: The HNB portal is available at http://www.hnbioinfo.de

HGVbase: a human sequence variation database emphasizing data quality and a broad spectrum of data sources.

HGVbase (Human Genome Variation database; http://hgvbase.cgb.ki.se, formerly known as HGBASE) is an academic effort to provide a high quality and non-redundant database of available genomic variation data of all types, mostly comprising single nucleotide polymorphisms (SNPs). Records include neutral polymorphisms as well as disease-related mutations. Online search tools facilitate data interrogation by sequence similarity and keyword queries, and searching by genome coordinates is now being implemented. Downloads are freely available in XML, Fasta, SRS, SQL and tagged-text file formats. Each entry is presented in the context of its surrounding sequence and many records are related to neighboring human genes and affected features therein. Population allele frequencies are included wherever available. Thorough semi-automated data checking ensures internal consistency and addresses common errors in the source information. To keep pace with recent growth in the field, we have developed tools for fully automated annotation. All variants have been uniquely mapped to the draft genome sequence and are referenced to positions in EMBL/GenBank files. Data utility is enhanced by provision of genotyping assays and functional predictions. Recent data structure extensions allow the capture of haplotype and genotype information, and a new initiative (along with BiSC and HUGO-MDI) aims to create a central repository for the broad collection of clinical mutations and associated disease phenotypes of interest.

The Spir actin organizers are involved in vesicle transport processes.

The p150-Spir protein, which was discovered as a phosphorylation target of the Jun N-terminal kinase, is an essential regulator of the polarization of the Drosophila oocyte. Spir proteins are highly conserved between species and belong to the family of Wiskott-Aldrich homology region 2 (WH2) proteins involved in actin organization. The C-terminal region of Spir encodes a zinc finger structure highly homologous to FYVE motifs. A region with high homology between the Spir family proteins is located adjacent (N-terminal) to the modified FYVE domain and is designated as "Spir-box." The Spir-box has sequence similarity to a region of rabphilin-3A, which mediates interaction with the small GTPase Rab3A. Coexpression of p150-Spir and green fluorescent protein-tagged Rab GTPases in NIH 3T3 cells revealed that the Spir protein colocalized specifically with the Rab11 GTPase, which is localized at the trans-Golgi network (TGN), post-Golgi vesicles, and the recycling endosome. The distinct Spir localization pattern was dependent on the integrity of the modified FYVE finger motif and the Spir-box. Overexpression of a mouse Spir-1 dominant interfering mutant strongly inhibited the transport of the vesicular stomatitis virus G (VSV G) protein to the plasma membrane. The viral protein was arrested in membrane structures, largely colocalizing with the TGN marker TGN46. Our findings that the Spir actin organizer is targeted to intracellular membrane structures by its modified FYVE zinc finger and is involved in vesicle transport processes provide a novel link between actin organization and intracellular transport.

Novel protein domains and repeats in Drosophila melanogaster: insights into structure, function, and evolution.

Sequence database searching methods such as BLAST, are invaluable for predicting molecular function on the basis of sequence similarities among single regions of proteins. Searches of whole databases however, are not optimized to detect multiple homologous regions within a single polypeptide. Here we have used the prospero algorithm to perform self-comparisons of all predicted Drosophila melanogaster gene products. Predicted repeats, and their homologs from all species, were analyzed further to detect hitherto unappreciated evolutionary relationships. Results included the identification of novel tandem repeats in the human X-linked retinitis pigmentosa type-2 gene product, repeated segments in cystinosin, associated with a defect in cystine transport, and 'nested' homologous domains in dysferlin, whose gene is mutated in limb girdle muscular dystrophy. Novel signaling domain families were found that may regulate the microtubule-based cytoskeleton and ubiquitin-mediated proteolysis, respectively. Two families of glycosyl hydrolases were shown to contain internal repetitions that hint at their evolution via a piecemeal, modular approach. In addition, three examples of fruit fly genes were detected with tandem exons that appear to have arisen via internal duplication. These findings demonstrate how completely sequenced genomes can be exploited to further understand the relationships between molecular structure, function, and evolution.

The phylogenetic distribution of frataxin indicates a role in iron-sulfur cluster protein assembly.

Much has been learned about the cellular pathology of Friedreich's ataxia, a recessive neurodegenerative disease resulting from insufficient expression of the mitochondrial protein frataxin. However, the biochemical function of frataxin has remained obscure, hampering attempts at therapeutic intervention. To predict functional interactions of frataxin with other proteins we investigated whether its gene specifically co-occurs with any other genes in sequenced genomes. In 56 available genomes we identified two genes with identical phylogenetic distributions to the frataxin/cyaY gene: hscA and hscB/JAC1. These genes have not only emerged in the same evolutionary lineage as the frataxin gene, they have also been lost at least twice with it, and they have been horizontally transferred with it in the evolution of the mitochondria. The proteins encoded by hscA and hscB, the chaperone HSP66 and the co-chaperone HSP20, have been shown to be required for the synthesis of 2Fe-2S clusters on ferredoxin in proteobacteria. JAC1, an ortholog of hscB, and SSQ1, a paralog of hscA, have been shown to be required for iron-sulfur cluster assembly in mitochondria of Saccharomyces cerevisiae. Combining data on the co-occurrence of genes in genomes with experimental and predicted cellular localization data of their proteins supports the hypothesis that frataxin is directly involved in iron-sulfur cluster protein assembly. They indicate that frataxin is specifically involved in the same sub-process as HSP20/Jac1p.

Evolution of tuf genes: ancient duplication, differential loss and gene conversion.

The tuf gene of eubacteria, encoding the EF-tu elongation factor, was duplicated early in the evolution of the taxon. Phylogenetic and genomic location analysis of 20 complete eubacterial genomes suggests that this ancient duplication has been differentially lost and maintained in eubacteria.

XplorMed: a tool for exploring MEDLINE abstracts.

The most frequent access to the MEDLINE database of scientific abstracts is by keyword search. However, this is often not sufficient because although the user might find all the useful abstracts, these are buried in hundreds that are irrelevant. The exploratory tool XplorMed has been developed to analyse the result of any MEDLINE query. It suggests main groups of related topics and documents, sparing the user the need of reading all abstracts.

Comparison of ARM and HEAT protein repeats.

ARM and HEAT motifs are tandemly repeated sequences of approximately 50 amino acid residues that occur in a wide variety of eukaryotic proteins. An exhaustive search of sequence databases detected new family members and revealed that at least 1 in 500 eukaryotic protein sequences contain such repeats. It also rendered the similarity between ARM and HEAT repeats, believed to be evolutionarily related, readily apparent. All the proteins identified in the database searches could be clustered by sequence similarity into four groups: canonical ARM-repeat proteins and three groups of the more divergent HEAT-repeat proteins. This allowed us to build improved sequence profiles for the automatic detection of repeat motifs. Inspection of these profiles indicated that the individual repeat motifs of all four classes share a common set of seven highly conserved hydrophobic residues, which in proteins of known three-dimensional structure are buried within or between repeats. However, the motifs differ at several specific residue positions, suggesting important structural or functional differences among the classes. Our results illustrate that ARM and HEAT-repeat proteins, while having a common phylogenetic origin, have since diverged significantly. We discuss evolutionary scenarios that could account for the great diversity of repeats observed.

Frameshift peptide-derived T-cell epitopes: a source of novel tumor-specific antigens.

Microsatellite instability (MSI) caused by defective DNA mismatch repair (MMR) is a hallmark of hereditary nonpolyposis colorectal cancers (HNPCC) but also occurs in about 15% of sporadic tumors. If instability affects microsatellites in coding regions, translational frameshifts lead to truncated proteins often marked by unique frameshift peptide sequences at their C-terminus. Since MSI tumors show enhanced lymphocytic infiltration and our previous analysis identified numerous coding mono- and dinucleotide repeat-bearing candidate genes as targets of genetic instability, we examined the role of frameshift peptides in triggering cellular immune responses. Using peptide pulsed autologous CD40-activated B cells, we have generated cytotoxic T lymphocytes (CTL) that specifically recognize HLA-A2.1-restricted peptides derived from frameshift sequences. Among 16 frameshift peptides predicted from mutations in 8 different genes, 3 peptides conferred specific lysis of target cells exogenously loaded with cognate peptide. One peptide derived from a (-1) frameshift mutation in the TGFbetaIIR gene gave rise to a CTL bulk culture capable of lysing the MSI colorectal cancer cell line HCT116 carrying this frameshift mutation. Given the huge number of human coding microsatellites and assuming only a fraction being mutated and encoding immunologically relevant peptides in MSI tumors, frameshift protein sequences represent a novel subclass of tumor-specific antigens. It is tempting to speculate that a frameshift peptide-directed vaccination approach not only could offer new treatment modalities for existing MSI tumors but also might benefit asymptomatic at-risk individuals in HNPCC families by a prophylactic vaccination strategy.

Systematic identification of genes with coding microsatellites mutated in DNA mismatch repair-deficient cancer cells.

Microsatellite instability (MSI) caused by deficient DNA mismatch-repair functions is a hallmark of cancers associated with the hereditary nonpolyposis colorectal cancer (HNPCC) syndrome but is also found in about 15% of all sporadic tumors. Most affected microsatellites reside in untranslated intergenic or intronic sequences. However, recently few genes with coding microsatellites were also shown to be mutational targets in MSI-positive cancers and might represent important mutation targets in their pathogenesis. The systematic identification of such genes and the analysis of their mutation frequency in MSI-positive cancers might thus reveal major clues to their functional role in MSI-associated carcinogenesis. We therefore initiated a systematic database search in 33,595 distinctly annotated human genes and identified 17,654 potentially coding mononucleotide repeats (cMNRs) and 2,028 coding dinucleotide repeats (cDNRs), which consist of n > or = 6 and n > or = 4 repeat units, respectively. Expression pattern and mutation frequency of 19 of these genes with the longest repeats were compared between DNA mismatch repair-deficient (MSI(+)) and proficient (MSS) cancer cells. Instability frequencies in these coding microsatellite genes ranged from 10% to 100% in MSI-H tumor cells, whereas MSS cancer cells did not show mutations. RT-PCR analysis further showed that most of the affected genes (10/15) were highly expressed in tumor cells. The approach outlined here identified a new set of genes frequently affected by mutations in MSI-positive tumor cells. It will lead to novel and highly specific diagnostic and therapeutic targets for microsatellite unstable cancers.