Hecatomb: an integrated software platform for viral metagenomics.

BACKGROUND: Modern sequencing technologies offer extraordinary opportunities for virus discovery and virome analysis. Annotation of viral sequences from metagenomic data requires a complex series of steps to ensure accurate annotation of individual reads and assembled contigs. In addition, varying study designs will require project-specific statistical analyses. FINDINGS: Here we introduce Hecatomb, a bioinformatic platform coordinating commonly used tasks required for virome analysis. Hecatomb means "a great sacrifice." In this setting, Hecatomb is "sacrificing" false-positive viral annotations using extensive quality control and tiered-database searches. Hecatomb processes metagenomic data obtained from both short- and long-read sequencing technologies, providing annotations to individual sequences and assembled contigs. Results are provided in commonly used data formats useful for downstream analysis. Here we demonstrate the functionality of Hecatomb through the reanalysis of a primate enteric and a novel coral reef virome. CONCLUSION: Hecatomb provides an integrated platform to manage many commonly used steps for virome characterization, including rigorous quality control, host removal, and both read- and contig-based analysis. Each step is managed using the Snakemake workflow manager with dependency management using Conda. Hecatomb outputs several tables properly formatted for immediate use within popular data analysis and visualization tools, enabling effective data interpretation for a variety of study designs. Hecatomb is hosted on GitHub (github.com/shandley/hecatomb) and is available for installation from Bioconda and PyPI.

Gemykibivirus detection in acute encephalitis patients from Nepal.

Acute Encephalitis Syndrome (AES) causes significant morbidity and mortality worldwide. In Nepal, Japanese encephalitis virus (JEV) accounts for ~ 5-20% of AES cases, but ~75% of AES cases are of unknown etiology. We identified a gemykibivirus in CSF collected in 2020 from a male child with AES using metagenomic next-generation sequencing. Gemykibiviruses are single stranded, circular DNA viruses in the family Genomoviridae. The complete genome of 2211 nucleotides was sequenced which shared 98.69% nucleotide identity to its closest relative, Human associated gemykibivirus 2 isolate SAfia-449D. Two real-time PCR assays were designed, and screening of 337 CSF and 164 serum samples from AES patients in Nepal collected in 2020 and 2022 yielded 11 CSF and 1 serum sample that were positive in both PCR assays. Complete genomes of 7 of the positives were sequenced. These results identify a candidate etiologic agent of encephalitis in Nepal.

Draft genomes of 12 Bifidobacterium isolates from human IBD fecal samples.

Twelve Bifidobacterium strains were isolated from fecal samples of inflammatory bowel disease patients and matched "household control" individuals. These include the species Bifidobacterium adolescentis, Bifidobacterium animalis, Bifidobacterium breve, Bifidobacterium catenulatum, Bifidobacterium longum, and Bifidobacterium pseudocatenulatum.

Human enteric viruses autonomously shape inflammatory bowel disease phenotype through divergent innate immunomodulation.

Altered enteric microorganisms in concert with host genetics shape inflammatory bowel disease (IBD) phenotypes. However, insight is limited to bacteria and fungi. We found that eukaryotic viruses and bacteriophages (collectively, the virome), enriched from non-IBD, noninflamed human colon resections, actively elicited atypical anti-inflammatory innate immune programs. Conversely, ulcerative colitis or Crohn's disease colon resection viromes provoked inflammation, which was successfully dampened by non-IBD viromes. The IBD colon tissue virome was perturbed, including an increase in the enterovirus B species of eukaryotic picornaviruses, not previously detected in fecal virome studies. Mice humanized with non-IBD colon tissue viromes were protected from intestinal inflammation, whereas IBD virome mice exhibited exacerbated inflammation in a nucleic acid sensing-dependent fashion. Furthermore, there were detrimental consequences for IBD patient-derived intestinal epithelial cells bearing loss-of-function mutations within virus sensor MDA5 when exposed to viromes. Our results demonstrate that innate recognition of IBD or non-IBD human viromes autonomously influences intestinal homeostasis and disease phenotypes. Thus, perturbations in the intestinal virome, or an altered ability to sense the virome due to genetic variation, contribute to the induction of IBD. Harnessing the virome may offer therapeutic and biomarker potential.

Divergent Enteroviruses from Macaques with Chronic Diarrhea.

We report the draft genome sequences of five novel members of the family Picornaviridae that were isolated from the stool of rhesus macaques (Macaca mulatta) with chronic diarrhea. The strains were named NOLA-1 through NOLA-5 because the macaques were residents of the Tulane National Primate Research Center.

Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders.

Understanding the tissue-specific genetic controls of protein levels is essential to uncover mechanisms of post-transcriptional gene regulation. In this study, we generated a genomic atlas of protein levels in three tissues relevant to neurological disorders (brain, cerebrospinal fluid and plasma) by profiling thousands of proteins from participants with and without Alzheimer's disease. We identified 274, 127 and 32 protein quantitative trait loci (pQTLs) for cerebrospinal fluid, plasma and brain, respectively. cis-pQTLs were more likely to be tissue shared, but trans-pQTLs tended to be tissue specific. Between 48.0% and 76.6% of pQTLs did not co-localize with expression, splicing, DNA methylation or histone acetylation QTLs. Using Mendelian randomization, we nominated proteins implicated in neurological diseases, including Alzheimer's disease, Parkinson's disease and stroke. This first multi-tissue study will be instrumental to map signals from genome-wide association studies onto functional genes, to discover pathways and to identify drug targets for neurological diseases.

Multi-Omics Analyses Show Disease, Diet, and Transcriptome Interactions With the Virome.

BACKGROUND & AIMS: The gut virome includes eukaryotic viruses and bacteriophages that can shape the gut bacterial community and elicit host responses. The virome can be implicated in diseases, such as irritable bowel syndrome (IBS), where gut bacteria play an important role in pathogenesis. We provide a comprehensive and longitudinal characterization of the virome, including DNA and RNA viruses and paired multi-omics data in a cohort of healthy subjects and patients with IBS. METHODS: We selected 2 consecutive stool samples per subject from a longitudinal study cohort and performed metagenomic sequencing on DNA and RNA viruses after enriching for viral-like particles. Viral sequence abundance was evaluated over time, as well as in the context of diet, bacterial composition and function, metabolite levels, colonic gene expression, host genetics, and IBS subsets. RESULTS: We found that the gut virome was temporally stable and correlated with the colonic transcriptome. We identified IBS-subset-specific changes in phage populations; Microviridae, Myoviridae, and Podoviridae species were elevated in diarrhea-predominant IBS, and other Microviridae and Myoviridae species were elevated in constipation-predominant IBS compared to healthy controls. We identified correlations between subsets of the virome and bacterial composition (unclassifiable "dark matter" and phages) and diet (eukaryotic viruses). CONCLUSIONS: We found that the gut virome is stable over time but varies among subsets of patients with IBS. It can be affected by diet and potentially influences host function via interactions with gut bacteria and/or altering host gene expression.

Alzheimer's disease alters oligodendrocytic glycolytic and ketolytic gene expression.

INTRODUCTION: Sporadic Alzheimer's disease (AD) is strongly correlated with impaired brain glucose metabolism, which may affect AD onset and progression. Ketolysis has been suggested as an alternative pathway to fuel the brain. METHODS: RNA-seq profiles of post mortem AD brains were used to determine whether dysfunctional AD brain metabolism can be determined by impairments in glycolytic and ketolytic gene expression. Data were obtained from the Knight Alzheimer's Disease Research Center (62 cases; 13 controls), Mount Sinai Brain Bank (110 cases; 44 controls), and the Mayo Clinic Brain Bank (80 cases; 76 controls), and were normalized to cell type: astrocytes, microglia, neurons, oligodendrocytes. RESULTS: In oligodendrocytes, both glycolytic and ketolytic pathways were significantly impaired in AD brains. Ketolytic gene expression was not significantly altered in neurons, astrocytes, and microglia. DISCUSSION: Oligodendrocytes may contribute to brain hypometabolism observed in AD. These results are suggestive of a potential link between hypometabolism and dysmyelination in disease physiology. Additionally, ketones may be therapeutic in AD due to their ability to fuel neurons despite impaired glycolytic metabolism.

TREM2 activation on microglia promotes myelin debris clearance and remyelination in a model of multiple sclerosis.

Multiple sclerosis (MS) is an inflammatory, demyelinating, and neurodegenerative disease of the central nervous system (CNS) triggered by autoimmune mechanisms. Microglia are critical for the clearance of myelin debris in areas of demyelination, a key step to allow remyelination. TREM2 is expressed by microglia and promotes microglial survival, proliferation, and phagocytic activity. Herein we demonstrate that TREM2 was highly expressed on myelin-laden phagocytes in active demyelinating lesions in the CNS of subjects with MS. In gene expression studies, macrophages from subjects with TREM2 genetic deficiency displayed a defect in phagocytic pathways. Treatment with a new TREM2 agonistic antibody promoted the clearance of myelin debris in the cuprizone model of CNS demyelination. Effects included enhancement of myelin uptake and degradation, resulting in accelerated myelin debris removal by microglia. Most importantly, antibody-dependent TREM2 activation on microglia increased density of oligodendrocyte precursors in areas of demyelination, as well as the formation of mature oligodendrocytes thus enhancing remyelination and axonal integrity. These results are relevant as they propose TREM2 on microglia as a potential new target to promote remyelination.

Improving Characterization of Understudied Human Microbiomes Using Targeted Phylogenetics.

Whole-genome bacterial sequences are required to better understand microbial functions, niche-specific bacterial metabolism, and disease states. Although genomic sequences are available for many of the human-associated bacteria from commonly tested body habitats (e.g., feces), as few as 13% of bacterium-derived reads from other sites such as the skin map to known bacterial genomes. To facilitate a better characterization of metagenomic shotgun reads from underrepresented body sites, we collected over 10,000 bacterial isolates originating from 14 human body habitats, identified novel taxonomic groups based on full-length 16S rRNA gene sequences, clustered the sequences to ensure that no individual taxonomic group was overselected for sequencing, prioritized bacteria from underrepresented body sites (such as skin and respiratory and urinary tracts), and sequenced and assembled genomes for 665 new bacterial strains. Here, we show that addition of these genomes improved read mapping rates of Human Microbiome Project (HMP) metagenomic samples by nearly 30% for the previously underrepresented phylum Fusobacteria, and 27.5% of the novel genomes generated here had high representation in at least one of the tested HMP samples, compared to 12.5% of the sequences in the public databases, indicating an enrichment of useful novel genomic sequences resulting from the prioritization procedure. As our understanding of the human microbiome continues to improve and to enter the realm of therapy developments, targeted approaches such as this to improve genomic databases will increase in importance from both an academic and a clinical perspective.IMPORTANCE The human microbiome plays a critically important role in health and disease, but current understanding of the mechanisms underlying the interactions between the varying microbiome and the different host environments is lacking. Having access to a database of fully sequenced bacterial genomes provides invaluable insights into microbial functions, but currently sequenced genomes for the human microbiome have largely come from a limited number of body sites (primarily feces), while other sites such as the skin, respiratory tract, and urinary tract are underrepresented, resulting in as little as 13% of bacterium-derived reads mapping to known bacterial genomes. Here, we sequenced and assembled 665 new bacterial genomes, prioritized from a larger database to select underrepresented body sites and bacterial taxa in the existing databases. As a result, we substantially improve mapping rates for samples from the Human Microbiome Project and provide an important contribution to human bacterial genomic databases for future studies.

The TMEM106B FTLD-protective variant, rs1990621, is also associated with increased neuronal proportion.

Apart from amyloid ß deposition and tau neurofibrillary tangles, Alzheimer's disease (AD) is a neurodegenerative disorder characterized by neuronal loss and astrocytosis in the cerebral cortex. The goal of this study is to investigate genetic factors associated with the neuronal proportion in health and disease. To identify cell-autonomous genetic variants associated with neuronal proportion in cortical tissues, we inferred cellular population structure from bulk RNA-Seq derived from 1536 individuals. We identified the variant rs1990621 located in the TMEM106B gene region as significantly associated with neuronal proportion (p value = 6.40 × 10-07) and replicated this finding in an independent dataset (p value = 7.41 × 10-04) surpassing the genome-wide threshold in the meta-analysis (p value = 9.42 × 10-09). This variant is in high LD with the TMEM106B non-synonymous variant p.T185S (rs3173615; r2 = 0.98) which was previously identified as a protective variant for frontotemporal lobar degeneration (FTLD). We stratified the samples by disease status, and discovered that this variant modulates neuronal proportion not only in AD cases, but also several neurodegenerative diseases and in elderly cognitively healthy controls. Furthermore, we did not find a significant association in younger controls or schizophrenia patients, suggesting that this variant might increase neuronal survival or confer resilience to the neurodegenerative process. The single variant and gene-based analyses also identified an overall genetic association between neuronal proportion, AD and FTLD risk. These results suggest that common pathways are implicated in these neurodegenerative diseases, that implicate neuronal survival. In summary, we identified a protective variant in the TMEM106B gene that may have a neuronal protection effect against general aging, independent of disease status, which could help elucidate the relationship between aging and neuronal survival in the presence or absence of neurodegenerative disorders. Our findings suggest that TMEM106B could be a potential target for neuronal protection therapies to ameliorate cognitive and functional deficits.

A single-nuclei RNA sequencing study of Mendelian and sporadic AD in the human brain.

BACKGROUND: Alzheimer's disease (AD) is the most common form of dementia. This neurodegenerative disorder is associated with neuronal death and gliosis heavily impacting the cerebral cortex. AD has a substantial but heterogeneous genetic component, presenting both Mendelian and complex genetic architectures. Using bulk RNA-seq from the parietal lobes and deconvolution methods, we previously reported that brains exhibiting different AD genetic architecture exhibit different cellular proportions. Here, we sought to directly investigate AD brain changes in cell proportion and gene expression using single-cell resolution. METHODS: We generated unsorted single-nuclei RNA sequencing data from brain tissue. We leveraged the tissue donated from a carrier of a Mendelian genetic mutation, PSEN1 p.A79V, and two family members who suffer from sporadic AD, but do not carry any autosomal mutations. We evaluated alternative alignment approaches to maximize the titer of reads, genes, and cells with high quality. In addition, we employed distinct clustering strategies to determine the best approach to identify cell clusters that reveal neuronal and glial cell types and avoid artifacts such as sample and batch effects. We propose an approach to cluster cells that reduces biases and enable further analyses. RESULTS: We identified distinct types of neurons, both excitatory and inhibitory, and glial cells, including astrocytes, oligodendrocytes, and microglia, among others. In particular, we identified a reduced proportion of excitatory neurons in the Mendelian mutation carrier, but a similar distribution of inhibitory neurons. Furthermore, we investigated whether single-nuclei RNA-seq from the human brains recapitulate the expression profile of disease-associated microglia (DAM) discovered in mouse models. We also determined that when analyzing human single-nuclei data, it is critical to control for biases introduced by donor-specific expression profiles. CONCLUSION: We propose a collection of best practices to generate a highly detailed molecular cell atlas of highly informative frozen tissue stored in brain banks. Importantly, we have developed a new web application to make this unique single-nuclei molecular atlas publicly available.

TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers.

BACKGROUND: Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript. METHODS: We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million. RESULTS: The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007). CONCLUSION: Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

Polycomb repressive complex 1 modifies transcription of active genes.

This study examines the role of Polycomb repressive complex 1 (PRC1) at active genes. The PRC1 and PRC2 complexes are crucial for epigenetic silencing during development of an organism. They are recruited to Polycomb response elements (PREs) and establish silenced domains over several kilobases. Recent studies show that PRC1 is also directly recruited to active genes by the cohesin complex. Cohesin participates broadly in control of gene transcription, but it is unknown whether cohesin-recruited PRC1 also plays a role in transcriptional control of active genes. We address this question using genome-wide RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq). The results show that PRC1 influences transcription of active genes, and a significant fraction of its effects are likely direct. The roles of different PRC1 subunits can also vary depending on the gene. Depletion of PRC1 subunits by RNA interference alters phosphorylation of RNA polymerase II (Pol II) and occupancy by the Spt5 pausing-elongation factor at most active genes. These effects on Pol II phosphorylation and Spt5 are likely linked to changes in elongation and RNA processing detected by nascent RNA-seq, although the mechanisms remain unresolved. The experiments also reveal that PRC1 facilitates association of Spt5 with enhancers and PREs. Reduced Spt5 levels at these regulatory sequences upon PRC1 depletion coincide with changes in Pol II occupancy and phosphorylation. Our findings indicate that, in addition to its repressive roles in epigenetic gene silencing, PRC1 broadly influences transcription of active genes and may suppress transcription of nonpromoter regulatory sequences.

Drosophila TDP-43 RNA-Binding Protein Facilitates Association of Sister Chromatid Cohesion Proteins with Genes, Enhancers and Polycomb Response Elements.

The cohesin protein complex mediates sister chromatid cohesion and participates in transcriptional control of genes that regulate growth and development. Substantial reduction of cohesin activity alters transcription of many genes without disrupting chromosome segregation. Drosophila Nipped-B protein loads cohesin onto chromosomes, and together Nipped-B and cohesin occupy essentially all active transcriptional enhancers and a large fraction of active genes. It is unknown why some active genes bind high levels of cohesin and some do not. Here we show that the TBPH and Lark RNA-binding proteins influence association of Nipped-B and cohesin with genes and gene regulatory sequences. In vitro, TBPH and Lark proteins specifically bind RNAs produced by genes occupied by Nipped-B and cohesin. By genomic chromatin immunoprecipitation these RNA-binding proteins also bind to chromosomes at cohesin-binding genes, enhancers, and Polycomb response elements (PREs). RNAi depletion reveals that TBPH facilitates association of Nipped-B and cohesin with genes and regulatory sequences. Lark reduces binding of Nipped-B and cohesin at many promoters and aids their association with several large enhancers. Conversely, Nipped-B facilitates TBPH and Lark association with genes and regulatory sequences, and interacts with TBPH and Lark in affinity chromatography and immunoprecipitation experiments. Blocking transcription does not ablate binding of Nipped-B and the RNA-binding proteins to chromosomes, indicating transcription is not required to maintain binding once established. These findings demonstrate that RNA-binding proteins help govern association of sister chromatid cohesion proteins with genes and enhancers.

Metagenomic Approach for Identification of the Pathogens Associated with Diarrhea in Stool Specimens.

The potential to rapidly capture the entire microbial community structure and/or gene content makes metagenomic sequencing an attractive tool for pathogen identification and the detection of resistance/virulence genes in clinical settings. Here, we assessed the consistency between PCR from a diagnostic laboratory, quantitative PCR (qPCR) from a research laboratory, 16S rRNA gene sequencing, and metagenomic shotgun sequencing (MSS) for Clostridium difficile identification in diarrhea stool samples. Twenty-two C. difficile-positive diarrhea samples identified by PCR and qPCR and five C. difficile-negative diarrhea controls were studied. C. difficile was detected in 90.9% of C. difficile-positive samples using 16S rRNA gene sequencing, and C. difficile was detected in 86.3% of C. difficile-positive samples using MSS. CFU inferred from qPCR analysis were positively correlated with the relative abundance of C. difficile from 16S rRNA gene sequencing (r(2) = -0.60) and MSS (r(2) = -0.55). C. difficile was codetected with Clostridium perfringens, norovirus, sapovirus, parechovirus, and anellovirus in 3.7% to 27.3% of the samples. A high load of Candida spp. was found in a symptomatic control sample in which no causative agents for diarrhea were identified in routine clinical testing. Beta-lactamase and tetracycline resistance genes were the most prevalent (25.9%) antibiotic resistance genes in these samples. In summary, the proof-of-concept study demonstrated that next-generation sequencing (NGS) in pathogen detection is moderately correlated with laboratory testing and is advantageous in detecting pathogens without a priori knowledge.

The effect of dietary resistant starch type 2 on the microbiota and markers of gut inflammation in rural Malawi children.

BACKGROUND: Resistant starch (RS) decreases intestinal inflammation in some settings. We tested the hypothesis that gut inflammation will be reduced with dietary supplementation with RS in rural Malawian children. Eighteen stunted 3-5-year-old children were supplemented with 8.5 g/day of RS type 2 for 4 weeks. The fecal samples were analyzed for the microbiota, the microbiome, short chain fatty acids, metabolome, and proteins indicative of inflammation before and after the intervention. Subjects served as their own controls. RESULTS: The consumption of RS changed the composition of the microbiota; at the phylum level Actinobacteria increased, while Firmicutes decreased. Among the most prevalent genera, Lactobacillus was increased and Roseburia, Blautia, and Lachnospiracea incertae sedis were decreased. The Shannon H index at the genus level decreased from 2.02 on the habitual diet and 1.76 after the introduction of RS (P < 0.01). Fecal acetate concentration decreased, and fecal propionate concentration increased after RS administration (-5.2 and 2.0 µmol/g, respectively). Fecal calprotectin increased from 29 ± 69 to 89 ± 49 µg/g (P = 0.003) after RS was given. The lipopolysaccharide biosynthesis pathway was upregulated. CONCLUSIONS: Our findings do not support the hypothesis that RS reduces gut inflammation in rural Malawian children.

Metagenomic analysis of double-stranded DNA viruses in healthy adults.

BACKGROUND: The Human Microbiome Project (HMP) was undertaken with the goal of defining microbial communities in and on the bodies of healthy individuals using high-throughput, metagenomic sequencing analysis. The viruses present in these microbial communities, the 'human virome', are an important aspect of the human microbiome that is particularly understudied in the absence of overt disease. We analyzed eukaryotic double-stranded DNA (dsDNA) viruses, together with dsDNA replicative intermediates of single-stranded DNA viruses, in metagenomic sequence data generated by the HMP. 706 samples from 102 subjects were studied, with each subject sampled at up to five major body habitats: nose, skin, mouth, vagina, and stool. Fifty-one individuals had samples taken at two or three time points 30 to 359 days apart from at least one of the body habitats. RESULTS: We detected an average of 5.5 viral genera in each individual. At least 1 virus was detected in 92% of the individuals sampled. These viruses included herpesviruses, papillomaviruses, polyomaviruses, adenoviruses, anelloviruses, parvoviruses, and circoviruses. Each individual had a distinct viral profile, demonstrating the high interpersonal diversity of the virome. Some components of the virome were stable over time. CONCLUSIONS: This study is the first to use high-throughput DNA sequencing to describe the diversity of eukaryotic dsDNA viruses in a large cohort of normal individuals who were sampled at multiple body sites. Our results show that the human virome is a complex component of the microbial flora. Some viruses establish long-term infections that may be associated with increased risk or possibly with protection from disease. A better understanding of the composition and dynamics of the virome may hold important keys to human health.

Exploration of bacterial community classes in major human habitats.

BACKGROUND: Determining bacterial abundance variation is the first step in understanding bacterial similarity between individuals. Categorization of bacterial communities into groups or community classes is the subsequent step in describing microbial distribution based on abundance patterns. Here, we present an analysis of the groupings of bacterial communities in stool, nasal, skin, vaginal and oral habitats in a healthy cohort of 236 subjects from the Human Microbiome Project. RESULTS: We identify distinct community group patterns in the anterior nares, four skin sites, and vagina at the genus level. We also confirm three enterotypes previously identified in stools. We identify two clusters with low silhouette values in most oral sites, in which bacterial communities are more homogeneous. Subjects sharing a community class in one habitat do not necessarily share a community class in another, except in the three vaginal sites and the symmetric habitats of the left and right retroauricular creases. Demographic factors, including gender, age, and ethnicity, significantly influence community composition in several habitats. Community classes in the vagina, retroauricular crease and stool are stable over approximately 200 days. CONCLUSION: The community composition, association of demographic factors with community classes, and demonstration of community stability deepen our understanding of the variability and dynamics of human microbiomes. This also has significant implications for experimental designs that seek microbial correlations with clinical phenotypes.

Nyamiviridae: proposal for a new family in the order Mononegavirales.

Nyamanini virus (NYMV) and Midway virus (MIDWV) are unclassified tick-borne agents that infect land birds and seabirds, respectively. The recent molecular characterization of both viruses confirmed their already known close serological relationship and revealed them to be nonsegmented, single- and negative-stranded RNA viruses that are clearly related to, but quite distinct from, members of the order Mononegavirales (bornaviruses, filoviruses, paramyxoviruses, and rhabdoviruses). A third agent, soybean cyst nematode virus 1 (SbCNV-1, previously named soybean cyst nematode nyavirus), was recently found to be an additional member of this new virus group. Here, we review the current knowledge about all three viruses and propose classifying them as members of a new mononegaviral family, Nyamiviridae.