Inter and intra-host diversity of RSV in hematopoietic stem cell transplant adults with normal and delayed viral clearance

Respiratory syncytial virus (RSV) infection in immunocompromised individuals often leads to prolonged illness, progression to severe lower respiratory tract infection, and even death. How the host immune environment of the hematopoietic stem cell transplant (HCT) adults can affect viral genetic variation during an acute infection is not understood well. In the present study, we performed whole genome sequencing of RSV/A or RSV/B from samples collected longitudinally from HCT adults with normal (<14 days) and delayed (≥14 days) RSV clearance who were enrolled in a ribavirin trial. We determined the inter-host and intra-host genetic variation of RSV and the effect of mutations on putative glycosylation sites. The inter-host variation of RSV is centered in the attachment (G) and fusion (F) glycoprotein genes followed by polymerase (L) and matrix (M) genes. Interestingly, the overall genetic variation was constant between normal and delayed clearance groups for both RSV/A and RSV/B. Intra-host variation primarily occurred in the G gene followed by non-structural protein (NS1) and L genes; however, gain or loss of stop codons and frameshift mutations appeared only in the G gene and only in the delayed viral clearance group. Potential gain or loss of O-linked glycosylation sites in the G gene occurred both in RSV/A and RSV/B isolates. For RSV F gene, loss of N-linked glycosylation site occurred in three RSV/B isolates within an antigenic epitope. Both oral and aerosolized ribavirin did not cause any mutations in the L gene. In summary, prolonged viral shedding and immune deficiency resulted in RSV variation, especially in structural mutations in the G gene, possibly associated with immune evasion. Therefore, sequencing and monitoring of RSV isolates from immunocompromised patients are crucial as they can create escape mutants that can impact the effectiveness of upcoming vaccines and treatments.

Modulating a prebiotic food source influences inflammation and immune-regulating gut microbes and metabolites: insights from the BE GONE trial

BACKGROUND

Accessible prebiotic foods hold strong potential to jointly target gut health and metabolic health in high-risk patients. The BE GONE trial targeted the gut microbiota of obese surveillance patients with a history of colorectal neoplasia through a straightforward bean intervention.

METHODS

This low-risk, non-invasive dietary intervention trial was conducted at MD Anderson Cancer Center (Houston, TX, USA). Following a 4-week equilibration, patients were randomized to continue their usual diet without beans (control) or to add a daily cup of study beans to their usual diet (intervention) with immediate crossover at 8-weeks. Stool and fasting blood were collected every 4 weeks to assess the primary outcome of intra and inter-individual changes in the gut microbiome and in circulating markers and metabolites within 8 weeks. This study was registered on ClinicalTrials.gov as NCT02843425, recruitment is complete and long-term follow-up continues.

FINDINGS

Of the 55 patients randomized by intervention sequence, 87% completed the 16-week trial, demonstrating an increase on-intervention in diversity [n = 48; linear mixed effect and 95% CI for inverse Simpson index: 0.16 (0.02, 0.30); p = 0.02] and shifts in multiple bacteria indicative of prebiotic efficacy, including increased Faecalibacterium, Eubacterium and Bifidobacterium (all p < 0.05). The circulating metabolome showed parallel shifts in nutrient and microbiome-derived metabolites, including increased pipecolic acid and decreased indole (all p < 0.002) that regressed upon returning to the usual diet. No significant changes were observed in circulating lipoproteins within 8 weeks; however, proteomic biomarkers of intestinal and systemic inflammatory response, fibroblast-growth factor-19 increased, and interleukin-10 receptor-α decreased (p = 0.01).

INTERPRETATION

These findings underscore the prebiotic and potential therapeutic role of beans to enhance the gut microbiome and to regulate host markers associated with metabolic obesity and colorectal cancer, while further emphasizing the need for consistent and sustainable dietary adjustments in high-risk patients.

FUNDING

This study was funded by the American Cancer Society.

Distinct transcriptomic profiles in children prior to the appearance of type 1 diabetes-linked islet autoantibodies and following enterovirus infection

Although the genetic basis and pathogenesis of type 1 diabetes have been studied extensively, how host responses to environmental factors might contribute to autoantibody development remains largely unknown. Here, we use longitudinal blood transcriptome sequencing data to characterize host responses in children within 12 months prior to the appearance of type 1 diabetes-linked islet autoantibodies, as well as matched control children. We report that children who present with insulin-specific autoantibodies first have distinct transcriptional profiles from those who develop GADA autoantibodies first. In particular, gene dosage-driven expression of GSTM1 is associated with GADA autoantibody positivity. Moreover, compared with controls, we observe increased monocyte and decreased B cell proportions 9-12 months prior to autoantibody positivity, especially in children who developed antibodies against insulin first. Lastly, we show that control children present transcriptional signatures consistent with robust immune responses to enterovirus infection, whereas children who later developed islet autoimmunity do not. These findings highlight distinct immune-related transcriptomic differences between case and control children prior to case progression to islet autoimmunity and uncover deficient antiviral response in children who later develop islet autoimmunity.

Tumor-resident Lactobacillus iners confer chemoradiation resistance through lactate-induced metabolic rewiring

Tumor microbiota can produce active metabolites that affect cancer and immune cell signaling, metabolism, and proliferation. Here, we explore tumor and gut microbiome features that affect chemoradiation response in patients with cervical cancer using a combined approach of deep microbiome sequencing, targeted bacterial culture, and in vitro assays. We identify that an obligate L-lactate-producing lactic acid bacterium found in tumors, Lactobacillus iners, is associated with decreased survival in patients, induces chemotherapy and radiation resistance in cervical cancer cells, and leads to metabolic rewiring, or alterations in multiple metabolic pathways, in tumors. Genomically similar L-lactate-producing lactic acid bacteria commensal to other body sites are also significantly associated with survival in colorectal, lung, head and neck, and skin cancers. Our findings demonstrate that lactic acid bacteria in the tumor microenvironment can alter tumor metabolism and lactate signaling pathways, causing therapeutic resistance. Lactic acid bacteria could be promising therapeutic targets across cancer types.

Gastrointestinal Infections Modulate the Risk for Insulin Autoantibodies as the First-Appearing Autoantibody in the TEDDY Study

OBJECTIVE

To investigate gastrointestinal infection episodes (GIEs) in relation to the appearance of islet autoantibodies in The Environmental Determinants of Diabetes in the Young (TEDDY) cohort.

RESEARCH DESIGN AND METHODS

GIEs on risk of autoantibodies against either insulin (IAA) or GAD (GADA) as the first-appearing autoantibody were assessed in a 10-year follow-up of 7,867 children. Stool virome was characterized in a nested case-control study.

RESULTS

GIE reports (odds ratio [OR] 2.17 [95% CI 1.39-3.39]) as well as Norwalk viruses found in stool (OR 5.69 [1.36-23.7]) at <1 year of age were associated with an increased IAA risk at 2-4 years of age. GIEs reported at age 1 to <2 years correlated with a lower risk of IAA up to 10 years of age (OR 0.48 [0.35-0.68]). GIE reports at any other age were associated with an increase in IAA risk (OR 2.04 for IAA when GIE was observed 12-23 months prior [1.41-2.96]). Impacts on GADA risk were limited to GIEs <6 months prior to autoantibody development in children <4 years of age (OR 2.16 [1.54-3.02]).

CONCLUSIONS

Bidirectional associations were observed. GIEs were associated with increased IAA risk when reported before 1 year of age or 12-23 months prior to IAA. Norwalk virus was identified as one possible candidate factor. GIEs reported during the 2nd year of life were associated with a decreased IAA risk.

Wastewater sequencing reveals community and variant dynamics of the collective human virome

Wastewater is a discarded human by-product, but its analysis may help us understand the health of populations. Epidemiologists first analyzed wastewater to track outbreaks of poliovirus decades ago, but so-called wastewater-based epidemiology was reinvigorated to monitor SARS-CoV-2 levels while bypassing the difficulties and pit falls of individual testing. Current approaches overlook the activity of most human viruses and preclude a deeper understanding of human virome community dynamics. Here, we conduct a comprehensive sequencing-based analysis of 363 longitudinal wastewater samples from ten distinct sites in two major cities. Critical to detection is the use of a viral probe capture set targeting thousands of viral species or variants. Over 450 distinct pathogenic viruses from 28 viral families are observed, most of which have never been detected in such samples. Sequencing reads of established pathogens and emerging viruses correlate to clinical data sets of SARS-CoV-2, influenza virus, and monkeypox viruses, outlining the public health utility of this approach. Viral communities are tightly organized by space and time. Finally, the most abundant human viruses yield sequence variant information consistent with regional spread and evolution. We reveal the viral landscape of human wastewater and its potential to improve our understanding of outbreaks, transmission, and its effects on overall population health.

Phage-bacteria dynamics during the first years of life revealed by trans-kingdom marker gene analysis

Humans are colonized with commensal bacteria soon after birth, and, while this colonization is affected by lifestyle and other factors, bacterial colonization proceeds through well-studied phases. However, less is known about phage communities in early human development due to small study sizes, inability to leverage large databases, and lack of appropriate bioinformatics tools. In this study, whole genome shotgun sequencing data from the TEDDY study, composed of 12,262 longitudinal samples from 887 children in 4 countries, is reanalyzed to assess phage and bacterial dynamics simultaneously. Reads from these samples were mapped to marker genes from both bacteria and a new database of tens of thousands of phage taxa from human microbiomes. We uncover that each child is colonized by hundreds of different phages during the early years, and phages are more transitory than bacteria. Participants' samples continually harbor new phage species over time whereas the diversification of bacterial species begins to saturate. Phage data improves the ability for machine learning models to discriminate samples by country. Finally, while phage populations were individual-specific, striking patterns arose from the larger dataset, showing clear trends of ecological succession amongst phages, which correlated well with putative host bacteria. Improved understanding of phage-bacterial relationships may reveal new means by which to shape and modulate the microbiome and its constituents to improve health and reduce disease, particularly in vulnerable populations where antibiotic use and/or other more drastic measures may not be advised.

Functional Genomics of Gastrointestinal Escherichia coli Isolated from Patients with Cancer and Diarrhea

We describe the epidemiology and clinical characteristics of 29 patients with cancer and diarrhea in whom Enteroaggregative Escherichia coli (EAEC) was initially identified by GI BioFire panel multiplex. E. coli strains were successfully isolated from fecal cultures in 14 of 29 patients. Six of the 14 strains were identified as EAEC and 8 belonged to other diverse E. coli groups of unknown pathogenesis. We investigated these strains by their adherence to human intestinal organoids, cytotoxic responses, antibiotic resistance profile, full sequencing of their genomes, and annotation of their functional virulome. Interestingly, we discovered novel and enhanced adherence and aggregative patterns for several diarrheagenic pathotypes that were not previously seen when co-cultured with immortalized cell lines. EAEC isolates displayed exceptional adherence and aggregation to human colonoids compared not only to diverse GI E. coli , but also compared to prototype strains of other diarrheagenic E. coli . Some of the diverse E. coli strains that could not be classified as a conventional pathotype also showed an enhanced aggregative and cytotoxic response. Notably, we found a high carriage rate of antibiotic resistance genes in both EAEC strains and diverse GI E. coli isolates and observed a positive correlation between adherence to colonoids and the number of metal acquisition genes carried in both EAEC and the diverse E. coli strains. This work indicates that E. coli from cancer patients constitute strains of remarkable pathotypic and genomic divergence, including strains of unknown disease etiology with unique virulomes. Future studies will allow for the opportunity to re-define E. coli pathotypes with greater diagnostic accuracy and into more clinically relevant groupings.

Abstract 649: Comparative effects of combustible cigarette versus electronic cigarette exposures on K-ras mutant lung cancer

Combustible cigarette smoking (CCS) is causally related to ~90% of all lung cancers and induces various tumor-initiating effects, some of which (e.g., inflammation) are not reversible even after smoking cessation. More recently, use of non-combustible smoking vectors, such as electronic cigarette vapors (ECV), have rapidly increased, especially among youth. Despite this emerging public health concern, the long-term impact of ECV exposure is poorly understood. Our lab sought to assess the safety of ECV and determine their influence on immune dysfunction and lung cancer progression by comparing the effects of CCS and ECV in an airway epithelial specific K-ras mutant mouse model of lung adenocarcinoma (CC-LR). Three cohorts of 6-week-old CC-LR mice were exposed to either room air, CCS, or ECV 5 days per week for 8 weeks. For CCS exposure, mice were exposed to CCS from 3R4F research cigarettes at a rate of 3 puffs/min for 2 hours each day using the SCIREQ integrated Cigarette Smoking Robot. For ECV exposure, mice were exposed to 72 mg/ml of liquid nicotine in 50% propylene glycol/vegetable glycerin solution at a rate of 3 puffs/min for 2 hours using the SCIREQ Electronic Nicotine Delivery System. While CCS led to a significant increase in tumor burden (25%), ECV showed no significant changes compared to room air-exposed control. Bronchoalveolar lavage fluid had significantly increased total lung immune cell infiltrates in both CCS (30%) and ECV (25%), particularly in macrophage and lymphocyte populations. Immunophenotyping of CCS and ECV exposed lungs displayed pronounced pro-tumor immunosuppressive phenotypes, characterized by significantly decreased CD4+ IFNγ+ and CD8+ GZMB+ cytotoxic T cells along with a significantly elevated CD4+ FOXP3+ regulatory T cells. Furthermore, the cytotoxic cytokine IFNγ was significantly reduced in lungs of CCS- (7-fold) and ECV-exposed (2-fold) mice compared to room air-exposed controls. We also found increased levels of the immunosuppressive cytokine IL-10 in lungs of CCS- and ECV-exposed mice. Previous studies have shown IL-10 to play an important role in microbiome-associated immune modulation, hence we performed a preliminary gut and lung microbiome study to assess microbiome dissimilarity in these cohorts. Taxonomic profiling via 16S rRNA gene sequencing of matched stool and lung samples showed differences in the relative abundance of several lung Proteobacteria spp. while gut Firmicutes, particularly Turicibacter and Ileibacterium, were increased by CCS and ECV. In conclusion, while both CCS and ECV increased immune suppression, ECV did not significantly promote tumorigenesis during this timeframe of intervention. Future studies probing differences in microbiome-modulated CCS and ECV immune dysfunction can help pave way for identification of new targets that foster new interception and early treatment strategies for K-ras mutant lung cancer.

Citation Format: Walter V. Velasco Torrez, Maria T. Grimaldo, Michael J. Clowers, Bo Yuan, Segundo del Aguila Soto, Iman Bouchelkia, Javier Eduardo Moreno Barragan, Umesh C. Karandikar, Joseph F. Petrosino, Florencia McAllister, Humam Kadara, Kristi Louise Hoffman, Seyed Javad Moghaddam. Comparative effects of combustible cigarette versus electronic cigarette exposures on K-ras mutant lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 649.

The Association between Caffeine Intake and the Colonic Mucosa-Associated Gut Microbiota in Humans-A Preliminary Investigation

We examined the association between caffeine and coffee intake and the community composition and structure of colonic microbiota. A total of 34 polyp-free adults donated 97 colonic biopsies. Microbial DNA was sequenced for the 16S rRNA gene V4 region. The amplicon sequence variant was assigned using DADA2 and SILVA. Food consumption was ascertained using a food frequency questionnaire. We compared the relative abundance of taxonomies by low (<82.9 mg) vs. high (≥82.9 mg) caffeine intake and by never or <2 cups vs. 2 cups vs. ≥3 cups coffee intake. False discovery rate-adjusted p values (q values) <0.05 indicated statistical significance. Multivariable negative binomial regression models were used to estimate the incidence rate ratio and its 95% confidence interval of having a non-zero count of certain bacteria by intake level. Higher caffeine and coffee intake was related to higher alpha diversity (Shannon index p < 0.001), higher relative abundance of Faecalibacterium and Alistipes, and lower relative abundance of Erysipelatoclostridium (q values < 0.05). After adjustment of vitamin B2 in multivariate analysis, the significant inverse association between Erysipelatoclostridium count and caffeine intake remained statistically significant. Our preliminary study could not evaluate other prebiotics in coffee.

Modeling nonsegmented negative-strand RNA virus (NNSV) transcription with ejective polymerase collisions and biased diffusion

Infections by non-segmented negative-strand RNA viruses (NNSV) are widely thought to entail gradient gene expression from the well-established existence of a single promoter at the 3' end of the viral genome and the assumption of constant transcriptional attenuation between genes. But multiple recent studies show viral mRNA levels in infections by respiratory syncytial virus (RSV), a major human pathogen and member of NNSV, that are inconsistent with a simple gradient. Here we integrate known and newly predicted phenomena into a biophysically reasonable model of NNSV transcription. Our model succeeds in capturing published observations of respiratory syncytial virus and vesicular stomatitis virus (VSV) mRNA levels. We therefore propose a novel understanding of NNSV transcription based on the possibility of ejective polymerase-polymerase collisions and, in the case of RSV, biased polymerase diffusion.

Butyrate-Producing Bacteria and Insulin Homeostasis: The Microbiome and Insulin Longitudinal Evaluation Study (MILES)

Gut microbiome studies have documented depletion of butyrate-producing taxa in type 2 diabetes. We analyzed associations between butyrate-producing taxa and detailed measures of insulin homeostasis, whose dysfunction underlies diabetes in 224 non-Hispanic Whites and 129 African Americans, all of whom completed an oral glucose tolerance test. Stool microbiome was assessed by whole-metagenome shotgun sequencing with taxonomic profiling. We examined associations among 36 butyrate-producing taxa (n = 7 genera and 29 species) and insulin sensitivity, insulin secretion, disposition index, insulin clearance, and prevalence of dysglycemia (prediabetes plus diabetes, 46% of cohort), adjusting for age, sex, BMI, and race. The genus Coprococcus was associated with higher insulin sensitivity (β = 0.14; P = 0.002) and disposition index (β = 0.12; P = 0.012) and a lower rate of dysglycemia (odds ratio [OR] 0.91; 95% CI 0.85-0.97; P = 0.0025). In contrast, Flavonifractor was associated with lower insulin sensitivity (β = -0.13; P = 0.004) and disposition index (β = -0.11; P = 0.04) and higher prevalence of dysglycemia (OR 1.22; 95% CI 1.08-1.38; P = 0.0013). Species-level analyses found 10 bacteria associated with beneficial directions of effects and two bacteria with adverse associations on insulin homeostasis and dysglycemia. Although most butyrate producers analyzed appear to be metabolically beneficial, this is not the case for all such bacteria, suggesting that microbiome-directed therapeutic measures to prevent or treat diabetes should be targeted to specific butyrate-producing taxa rather than all butyrate producers.

Maternal gut microbiota mediate intergenerational effects of high-fat diet on descendant social behavior

Dysbiosis of the maternal gut microbiome during pregnancy is associated with adverse neurodevelopmental outcomes. We previously showed that maternal high-fat diet (MHFD) in mice induces gut dysbiosis, social dysfunction, and underlying synaptic plasticity deficits in male offspring (F₁). Here, we reason that, if HFD-mediated changes in maternal gut microbiota drive offspring social deficits, then MHFD-induced dysbiosis in F₁ female MHFD offspring would likewise impair F₂ social behavior. Metataxonomic sequencing reveals reduced microbial richness among female F₁ MHFD offspring. Despite recovery of microbial richness among MHFD-descendant F₂ mice, they display social dysfunction. Post-weaning Limosilactobacillus reuteri treatment increases the abundance of short-chain fatty acid-producing taxa and rescues MHFD-descendant F₂ social deficits. L. reuteri exerts a sexually dimorphic impact on gut microbiota configuration, increasing discriminant taxa between female cohorts. Collectively, these results show multigenerational impacts of HFD-induced dysbiosis in the maternal lineage and highlight the potential of maternal microbiome-targeted interventions for neurodevelopmental disorders.

Metagenomes of rectal swabs in larger, advanced stage cervical cancers have enhanced mucus degrading functionalities and distinct taxonomic structure

Background

Gut microbiome community composition differs between cervical cancer (CC) patients and healthy controls, and increased gut diversity is associated with improved outcomes after treatment. We proposed that functions of specific microbial species adjoining the mucus layer may directly impact the biology of CC.

Method

Metagenomes of rectal swabs in 41 CC patients were examined by whole-genome shotgun sequencing to link taxonomic structures, molecular functions, and metabolic pathway to patient’s clinical characteristics.

Results

Significant association of molecular functions encoded by the metagenomes was found with initial tumor size and stage. Profiling of the molecular function abundances and their distributions identified 2 microbial communities co-existing in each metagenome but having distinct metabolism and taxonomic structures. Community A (Clostridia and Proteobacteria predominant) was characterized by high activity of pathways involved in stress response, mucus glycan degradation and utilization of degradation byproducts. This community was prevalent in patients with larger, advanced stage tumors. Conversely, community B (Bacteroidia predominant) was characterized by fast growth, active oxidative phosphorylation, and production of vitamins. This community was prevalent in patients with smaller, early-stage tumors.

Conclusions

In this study, enrichment of mucus degrading microbial communities in rectal metagenomes of CC patients was associated with larger, more advanced stage tumors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09997-0.

Temporal changes in gastrointestinal fungi and the risk of autoimmunity during early childhood: the TEDDY study

Fungal infections are a major health problem that often begin in the gastrointestinal tract. Gut microbe interactions in early childhood are critical for proper immune responses, yet there is little known about the development of the fungal population from infancy into childhood. Here, as part of the TEDDY (The Environmental Determinants of Diabetes in the Young) study, we examine stool samples of 888 children from 3 to 48 months and find considerable differences between fungi and bacteria. The metagenomic relative abundance of fungi was extremely low but increased while weaning from milk and formula. Overall fungal diversity remained constant over time, in contrast with the increase in bacterial diversity. Fungal profiles had high temporal variation, but there was less variation from month-to-month in an individual than among different children of the same age. Fungal composition varied with geography, diet, and the use of probiotics. Multiple Candida spp. were at higher relative abundance in children than adults, while Malassezia and certain food-associated fungi were lower in children. There were only subtle fungal differences associated with the subset of children that developed islet autoimmunity or type 1 diabetes. Having proper fungal exposures may be crucial for children to establish appropriate responses to fungi and limit the risk of infection: the data here suggests those gastrointestinal exposures are limited and variable.

Alterations of the gut microbial community structure and function with aging in the spontaneously hypertensive stroke prone rat

Gut dysbiosis, a pathological imbalance of bacteria, has been shown to contribute to the development of hypertension (HT), systemic- and neuro-inflammation, and blood-brain barrier (BBB) disruption in spontaneously hypertensive stroke prone rats (SHRSP). However, to date individual species that contribute to HT in the SHRSP model have not been identified. One potential reason, is that nearly all studies of the SHRSP gut microbiota have analyzed samples from rats with established HT. The goal of this study was to examine the SHRSP gut microbiota before, during, and after the onset of hypertension, and in normotensive WKY control rats over the same age range. We hypothesized that we could identify key microbes involved in the development of HT by comparing WKY and SHRSP microbiota during the pre-hypertensive state and longitudinally. Systolic blood pressure (SBP) was measured by tail-cuff plethysmography and fecal microbiota analyzed by16S rRNA gene sequencing. SHRSP showed significant elevations in SBP, as compared to WKY, beginning at 8 weeks of age (p < 0.05 at each time point). Bacterial community structure was significantly different between WKY and SHRSP as early as 4 weeks of age, and remained different throughout the study (p = 0.001-0.01). At the phylum level we observed significantly reduced Firmicutes and Deferribacterota, and elevated Bacteroidota, Verrucomicrobiota, and Proteobacteria, in pre-hypertensive SHRSP, as compared to WKY. At the genus level we identified 18 bacteria whose relative abundance was significantly different in SHRSP versus WKY at the pre-hypertensive ages of 4 or 6 weeks. In an attempt to further refine bacterial candidates that might contribute to the SHRSP phenotype, we compared the functional capacity of WKY versus SHRSP microbial communities. We identified significant differences in amino acid metabolism. Using untargeted metabolomics we found significant reductions in metabolites of the tryptophan-kynurenine pathway and increased indole metabolites in SHRSP versus WKY plasma. Overall, we provide further evidence that gut dysbiosis contributes to hypertension in the SHRSP model, and suggest for the first time the potential involvement of tryptophan metabolizing microbes.

Spatial Characteristics of Colonic Mucosa-Associated Gut Microbiota in Humans

Limited data exist on the spatial distribution of the colonic bacteria in humans. We collected the colonic biopsies from five segments of 27 polyp-free adults and collected feces from 13 of them. We sequenced the V4 region of the bacterial 16S rRNA gene using the MiSeq platform. The sequencing data were assigned to the amplicon sequence variant (ASV) using SILVA. Biodiversity and the relative abundance of the ASV were compared across the colonic segments and between the rectal and fecal samples. Bacterial functional capacity was assessed using Tax4fun. Each individual had a unique bacterial community composition (Weighted Bray-Curtis P value = 0.001). There were no significant differences in richness, evenness, community composition, and the taxonomic structure across the colon segments in all the samples. Firmicutes (47%), Bacteroidetes (39%), and Proteobacteria (6%) were the major phyla in all segments, followed by Verrucomicrobia, Fusobacteria, Desulfobacterota, and Actinobacteria. There were 15 genera with relative abundance > 1%, including Bacteroides, Faecalibacterium, Escherichia/Shigella, Sutterella, Akkermansia, Parabacteroides, Prevotella, Lachnoclostridium, Alistipes, Fusobacterium, Erysipelatoclostridium, and four Lachnospiraceae family members. Intra-individually, the community compositional dissimilarity was the greatest between the cecum and the rectum. There were significant differences in biodiversity and the taxonomic structure between the rectal and fecal bacteria. The bacterial community composition and structure were homogeneous across the large intestine in adults. The inter-individual variability of the bacteria was greater than inter-segment variability. The rectal and fecal bacteria differed in the community composition and structure.

Fully resolved assembly of Cryptosporidium parvum

Background

Cryptosporidium parvum is an apicomplexan parasite commonly found across many host species with a global infection prevalence in human populations of 7.6%. Understanding its diversity and genomic makeup can help in fighting established infections and prohibiting further transmission. The basis of every genomic study is a high-quality reference genome that has continuity and completeness, thus enabling comprehensive comparative studies.

Findings

Here, we provide a highly accurate and complete reference genome of Cryptosporidium parvum. The assembly is based on Oxford Nanopore reads and was improved using Illumina reads for error correction. We also outline how to evaluate and choose from different assembly methods based on 2 main approaches that can be applied to other Cryptosporidium species. The assembly encompasses 8 chromosomes and includes 13 telomeres that were resolved. Overall, the assembly shows a high completion rate with 98.4% single-copy BUSCO genes.

Conclusions

This high-quality reference genome of a zoonotic IIaA17G2R1 C. parvum subtype isolate provides the basis for subsequent comparative genomic studies across the Cryptosporidium clade. This will enable improved understanding of diversity, functional, and association studies.

Urinary Tract Infections in Children with Vesicoureteral Reflux Are Accompanied by Alterations in Urinary Microbiota and Metabolome Profiles

Children with vesicoureteral reflux (VUR) are at an increased risk of recurrent urinary tract infections (UTIs) and renal scarring. Gut microbiota are associated with disease phenotypes, but there has been no study that associates urinary microbiota (uMB) and metabolic profiles with VUR pathology. To identify dominant uMB genera and metabolites associated with UTIs in VUR, urine samples collected under sterile conditions underwent 16S ribosomal RNA sequencing (n = 49) and metabolomic analysis by mass spectrometry (n = 96). Alterations in uMB and metabolomic profiles in VUR patients suggest remodeling of urinary bacterial communities after UTIs: Dorea- and Escherichia-dominant uMB profiles were more frequently identified in participants with VUR. Prevotella- and Lactobacillus-dominant uMB profiles were more prevalent in controls (p < 0.001). Microbial composition varied based on recurrent febrile UTI status (p = 0.001). A total of 243 urinary metabolites involved in energy, amino acid, nucleotide, and lipid metabolism were altered in VUR patients with UTIs (p < 0.05). Importantly, VUR specimens revealed changes in the bacteria-associated metabolic pathways such as glutamate degradation, methyl-citrate cycle, and bile acid metabolism. PATIENT SUMMARY: Differences in urinary commensal bacteria and metabolites exist between children with and without vesicoureteral reflux (VUR). These changes may be utilized to identify patients at risk of VUR-associated kidney damage.

Dietary fiber and probiotics influence the gut microbiome and melanoma immunotherapy response

Gut bacteria modulate the response to immune checkpoint blockade (ICB) treatment in cancer, but the effect of diet and supplements on this interaction is not well studied. We assessed fecal microbiota profiles, dietary habits, and commercially available probiotic supplement use in melanoma patients and performed parallel preclinical studies. Higher dietary fiber was associated with significantly improved progression-free survival in 128 patients on ICB, with the most pronounced benefit observed in patients with sufficient dietary fiber intake and no probiotic use. Findings were recapitulated in preclinical models, which demonstrated impaired treatment response to anti–programmed cell death 1 (anti–PD-1)–based therapy in mice receiving a low-fiber diet or probiotics, with a lower frequency of interferon-γ–positive cytotoxic T cells in the tumor microenvironment. Together, these data have clinical implications for patients receiving ICB for cancer.

Human milk oligosaccharide DSLNT and gut microbiome in preterm infants predicts necrotising enterocolitis

OBJECTIVE

Necrotising enterocolitis (NEC) is a devastating intestinal disease primarily affecting preterm infants. The underlying mechanisms are poorly understood: mother's own breast milk (MOM) is protective, possibly relating to human milk oligosaccharide (HMO) and infant gut microbiome interplay. We investigated the interaction between HMO profiles and infant gut microbiome development and its association with NEC.

DESIGN

We performed HMO profiling of MOM in a large cohort of infants with NEC (n=33) with matched controls (n=37). In a subset of 48 infants (14 with NEC), we also performed longitudinal metagenomic sequencing of infant stool (n=644).

RESULTS

Concentration of a single HMO, disialyllacto-N-tetraose (DSLNT), was significantly lower in MOM received by infants with NEC compared with controls. A MOM threshold level of 241 nmol/mL had a sensitivity and specificity of 0.9 for NEC. Metagenomic sequencing before NEC onset showed significantly lower relative abundance of Bifidobacterium longum and higher relative abundance of Enterobacter cloacae in infants with NEC. Longitudinal development of the microbiome was also impacted by low MOM DSLNT associated with reduced transition into preterm gut community types dominated by Bifidobacterium spp and typically observed in older infants. Random forest analysis combining HMO and metagenome data before disease accurately classified 87.5% of infants as healthy or having NEC.

CONCLUSION

These results demonstrate the importance of HMOs and gut microbiome in preterm infant health and disease. The findings offer potential targets for biomarker development, disease risk stratification and novel avenues for supplements that may prevent life-threatening disease.

Drivers of transcriptional variance in human intestinal epithelial organoids

Human intestinal epithelial organoids (enteroids and colonoids) are tissue cultures used for understanding the physiology of the human intestinal epithelium. Here, we explored the effect on the transcriptome of common variations in culture methods, including extracellular matrix substrate, format, tissue segment, differentiation status, and patient heterogeneity. RNA-sequencing datasets from 276 experiments performed on 37 human enteroid and colonoid lines from 29 patients were aggregated from several groups in the Texas Medical Center. DESeq2 and gene set enrichment analysis (GSEA) were used to identify differentially expressed genes and enriched pathways. PERMANOVA, Pearson's correlation, and dendrogram analysis of the data originally indicated three tiers of influence of culture methods on transcriptomic variation: substrate (collagen vs. Matrigel) and format (3-D, transwell, and monolayer) had the largest effect; segment of origin (duodenum, jejunum, ileum, colon) and differentiation status had a moderate effect; and patient heterogeneity and specific experimental manipulations (e.g., pathogen infection) had the smallest effect. GSEA identified hundreds of pathways that varied between culture methods, such as IL1 cytokine signaling enriched in transwell versus monolayer cultures and E2F target genes enriched in collagen versus Matrigel cultures. The transcriptional influence of the format was furthermore validated in a synchronized experiment performed with various format-substrate combinations. Surprisingly, large differences in organoid transcriptome were driven by variations in culture methods such as format, whereas experimental manipulations such as infection had modest effects. These results show that common variations in culture conditions can have large effects on intestinal organoids and should be accounted for when designing experiments and comparing results between laboratories. Our data constitute the largest RNA-seq dataset interrogating human intestinal epithelial organoids.

Oral Health and the Altered Colonic Mucosa-Associated Gut Microbiota

BACKGROUND

Systemic diseases have been associated with oral health and gut microbiota. We examined the association between oral health and the community composition and structure of the adherent colonic gut microbiota.

METHODS

We obtained 197 snap-frozen colonic biopsies from 62 colonoscopy-confirmed polyp-free individuals. Microbial DNA was sequenced for the 16S rRNA V4 region using the Illumina MiSeq, and the sequences were assigned to the operational taxonomic unit based on SILVA. We used a questionnaire to ascertain tooth loss, gum disease, and lifestyle factors. We compared biodiversity and relative abundance of bacterial taxa based on the amount of tooth loss and the presence of gum disease. The multivariable negative binomial regression model for panel data was used to estimate the association between the bacterial count and oral health. False discovery rate-adjusted P value (q value) < .05 indicated statistical significance.

RESULTS

More tooth loss and gum disease were associated with lower bacterial alpha diversity. The relative abundance of Faecalibacterium was lower (q values < .05) with more tooth loss. The association was significant after adjusting for age, ethnicity, obesity, smoking, alcohol use, hypertension, diabetes, and the colon segment. The relative abundance of Bacteroides was higher in those with gum disease.

CONCLUSIONS

Oral health was associated with alteration in the community composition and structure of the adherent gut bacteria in the colon. The reduced anti-inflammatory Faecalibacterium in participants with more tooth loss may indicate systemic inflammation. Future studies are warranted to confirm our findings and investigate the systemic role of Faecalibacterium.

Oligonucleotide capture sequencing of the SARS-CoV-2 genome and subgenomic fragments from COVID-19 individuals

The newly emerged and rapidly spreading SARS-CoV-2 causes coronavirus disease 2019 (COVID-19). To facilitate a deeper understanding of the viral biology we developed a capture sequencing methodology to generate SARS-CoV-2 genomic and transcriptome sequences from infected patients. We utilized an oligonucleotide probe-set representing the full-length genome to obtain both genomic and transcriptome (subgenomic open reading frames [ORFs]) sequences from 45 SARS-CoV-2 clinical samples with varying viral titers. For samples with higher viral loads (cycle threshold value under 33, based on the CDC qPCR assay) complete genomes were generated. Analysis of junction reads revealed regions of differential transcriptional activity among samples. Mixed allelic frequencies along the 20kb ORF1ab gene in one sample, suggested the presence of a defective viral RNA species subpopulation maintained in mixture with functional RNA in one sample. The associated workflow is straightforward, and hybridization-based capture offers an effective and scalable approach for sequencing SARS-CoV-2 from patient samples.

Oral insulin immunotherapy in children at risk for type 1 diabetes in a randomised controlled trial

AIMS/HYPOTHESIS

Oral administration of antigen can induce immunological tolerance. Insulin is a key autoantigen in childhood type 1 diabetes. Here, oral insulin was given as antigen-specific immunotherapy before the onset of autoimmunity in children from age 6 months to assess its safety and immune response actions on immunity and the gut microbiome.

METHODS

A phase I/II randomised controlled trial was performed in a single clinical study centre in Germany. Participants were 44 islet autoantibody-negative children aged 6 months to 2.99 years who had a first-degree relative with type 1 diabetes and a susceptible HLA DR4-DQ8-containing genotype. Children were randomised 1:1 to daily oral insulin (7.5 mg with dose escalation to 67.5 mg) or placebo for 12 months using a web-based computer system. The primary outcome was immune efficacy pre-specified as induction of antibody or T cell responses to insulin and measured in a central treatment-blinded laboratory.

RESULTS

Randomisation was performed in 44 children. One child in the placebo group was withdrawn after the first study visit and data from 22 insulin-treated and 21 placebo-treated children were analysed. Oral insulin was well tolerated with no changes in metabolic variables. Immune responses to insulin were observed in children who received both insulin (54.5%) and placebo (66.7%), and the trial did not demonstrate an effect on its primary outcome (p = 0.54). In exploratory analyses, there was preliminary evidence that the immune response and gut microbiome were modified by the INS genotype Among children with the type 1 diabetes-susceptible INS genotype (n = 22), antibody responses to insulin were more frequent in insulin-treated (72.7%) as compared with placebo-treated children (18.2%; p = 0.03). T cell responses to insulin were modified by treatment-independent inflammatory episodes.

CONCLUSIONS/INTERPRETATION

The study demonstrated that oral insulin immunotherapy in young genetically at-risk children was safe, but was not associated with an immune response as predefined in the trial primary outcome. Exploratory analyses suggested that antibody responses to oral insulin may occur in children with a susceptible INS genotype, and that inflammatory episodes may promote the activation of insulin-responsive T cells.

TRIAL REGISTRATION

Clinicaltrials.gov NCT02547519 FUNDING: The main funding source was the German Center for Diabetes Research (DZD e.V.).

Restructuring the Gut Microbiota by Intermittent Fasting Lowers Blood Pressure

[Figure: see text].

Gut microbiome diversity is an independent predictor of survival in cervical cancer patients receiving chemoradiation

Diversity of the gut microbiome is associated with higher response rates for cancer patients receiving immunotherapy but has not been investigated in patients receiving radiation therapy. Additionally, current studies investigating the gut microbiome and outcomes in cancer patients may not have adjusted for established risk factors. Here, we sought to determine if diversity and composition of the gut microbiome was independently associated with survival in cervical cancer patients receiving chemoradiation. Our study demonstrates that the diversity of gut microbiota is associated with a favorable response to chemoradiation. Additionally, compositional variation among patients correlated with short term and long-term survival. Short term survivor fecal samples were significantly enriched in Porphyromonas, Porphyromonadaceae, and Dialister, whereas long term survivor samples were significantly enriched in Escherichia Shigella, Enterobacteriaceae, and Enterobacteriales. Moreover, analysis of immune cells from cervical tumor brush samples by flow cytometry revealed that patients with a high microbiome diversity had increased tumor infiltration of CD4+ lymphocytes as well as activated subsets of CD4 cells expressing ki67+ and CD69+ over the course of radiation therapy. Modulation of the gut microbiota before chemoradiation might provide an alternative way to enhance treatment efficacy and improve treatment outcomes in cervical cancer patients.

The gut microbiome contributes to blood-brain barrier disruption in spontaneously hypertensive stroke prone rats

In recent years, it has become apparent that the gut microbiome can influence the functioning and pathological states of organs and systems throughout the body. In this study, we tested the hypothesis that the gut microbiome has a major role in the disruption of the blood-brain barrier (BBB) in the spontaneously hypertensive stroke prone rats (SHRSP), an animal model for hypertensive cerebral small vessel disease (CSVD). Loss of BBB is thought to be an early and initiating component to the full expression of CSVD in animal models and humans. To test this hypothesis, newly born SHRSP pups were placed with foster dams of the SHRSP strain or dams of the WKY strain, the control strain that does not demonstrate BBB dysfunction or develop hypertensive CSVD. Similarly, WKY pups were placed with foster dams of the same or opposite strain. The rationale for cross fostering is that the gut microbiomes are shaped by environmental bacteria of the foster dam and the nesting surroundings. Analysis of the bacterial genera in feces, using 16S rRNA analysis, demonstrated that the gut microbiome in the rat pups was influenced by the foster dam. SHRSP offspring fostered on WKY dams had systolic blood pressures (SBPs) that were significantly decreased by 26 mmHg (P < .001) from 16-20 weeks, compared to SHRSP offspring fostered on SHRSP dams. Similarly WKY offspring fostered on SHRSP dams had significantly increased SBP compared to WKY offspring fostered on WKY dams, although the magnitude of SBP change was not as robust. At ~20 weeks of age, rats fostered on SHRSP dams showed enhanced inflammation in distal ileum regardless of the strain of the offspring. Disruption of BBB integrity, an early marker of CSVD onset, was improved in SHRSPs that were fostered on WKY dams when compared to the SHRSP rats fostered on SHRSP dams. Although SHRSP is a genetic model for CSVD, environmental factors such as the gut microbiota of the foster dam have a major influence in the loss of BBB integrity.

Fructan-sensitive children with irritable bowel syndrome have distinct gut microbiome signatures

BACKGROUND

Dietary fructans may worsen gastrointestinal symptoms in children with irritable bowel syndrome (IBS).

AIM

To determine whether gut microbiome composition and function are associated with childhood IBS fructan-induced symptoms.

METHODS

Faecal samples were collected from 38 children aged 7-17 years with paediatric Rome III IBS, who previously completied a double-blind, randomised, placebo-controlled crossover (fructan vs maltodextrin) trial. Fructan sensitivity was defined as an increase of ≥30% in abdominal pain frequency during the fructan diet. Gut microbial composition was determined via 16Sv4 rDNA sequencing. LEfSe evaluated taxonomic composition differences. Tax4Fun2 predicted microbial fructan metabolic pathways.

RESULTS

At baseline, 17 fructan-sensitive (vs 21 fructan-tolerant) subjects had lower alpha diversity (q < 0.05) and were enriched in the genus Holdermania. In contrast, fructan-tolerant subjects were enriched in 14 genera from the class Clostridia. During the fructan diet, fructan-sensitive (vs tolerant) subjects were enriched in both Agathobacter (P = 0.02) and Cyanobacteria (P = 0.0001). In contrast, fructan-tolerant subjects were enriched in three genera from the Clostridia class. Comparing the fructan vs maltodextrin diet, fructan-sensitive subjects had a significantly increased relative abundance of Bifidobacterium (P = 0.02) while fructan-tolerant subjects had increased Anaerostipes (P = 0.03) during the fructan diet. Only fructan-sensitive subjects had a trend towards increased predicted β-fructofuranosidase during the fructan vs maltodextrin diet.

CONCLUSIONS

Fructan-sensitive children with IBS have distinct gut microbiome signatures. These microbiome signatures differ both at baseline and in response to a fructan challenge.

A study of innate immune kinetics reveals a role for a chloride transporter in a virulent Francisella tularensis type B strain

Tularemia is a zoonotic disease of global proportions. Francisella tularensis subspecies tularensis (type A) and holarctica (type B) cause disease in healthy humans, with type A infections resulting in higher mortality. Repeated passage of a type B strain in the mid-20th century generated the Live Vaccine Strain (LVS). LVS remains unlicensed, does not protect against high inhalational doses of type A, and its exact mechanisms of attenuation are poorly understood. Recent data suggest that live attenuated vaccines derived from type B may cross-protect against type A. However, there is a dearth of knowledge regarding virulent type B pathogenesis and its capacity to stimulate the host's innate immune response. We therefore sought to increase our understanding of virulent type B in vitro characteristics using strain OR96-0246 as a model. Adding to our knowledge of innate immune kinetics in macrophages following infection with virulent type B, we observed robust replication of strain OR96-0246 in murine and human macrophages, reduced expression of pro-inflammatory cytokine genes from "wild type" type B-infected macrophages compared to LVS, and delayed macrophage cell death suggesting that virulent type B may suppress macrophage activation. One disruption in LVS is in the gene encoding the chloride transporter ClcA. We investigated the role of ClcA in macrophage infection and observed a replication delay in a clcA mutant. Here, we propose its role in acid tolerance. A greater understanding of LVS attenuation may reveal new mechanisms of pathogenesis and inform strategies toward the development of an improved vaccine against tularemia.

Houston hurricane Harvey health (Houston-3H) study: assessment of allergic symptoms and stress after hurricane Harvey flooding

BACKGROUND

In August 2017, Hurricane Harvey caused unprecedented flooding across the greater Houston area. Given the potential for widespread flood-related exposures, including mold and sewage, and the emotional and mental toll caused by the flooding, we sought to evaluate the short- and long-term impact of flood-related exposures on the health of Houstonians. Our objectives were to assess the association of flood-related exposures with allergic symptoms and stress among Houston-area residents at two time points: within approximately 30 days (T1) and 12 months (T2) after Hurricane Harvey's landfall.

METHODS

The Houston Hurricane Harvey Health (Houston-3H) Study enrolled a total of 347 unique participants from four sites across Harris County at two times: within approximately 1-month of Harvey (T1, n = 206) and approximately 12-months after Harvey (T2, n = 266), including 125 individuals who participated at both time points. Using a self-administered questionnaire, participants reported details on demographics, flood-related exposures, and health outcomes, including allergic symptoms and stress.

RESULTS

The majority of participants reported hurricane-related flooding in their homes at T1 (79.1%) and T2 (87.2%) and experienced at least one allergic symptom after the hurricane (79.4% at T1 and 68.4% at T2). In general, flood-exposed individuals were at increased risk of upper respiratory tract allergic symptoms, reported at both the T1 and T2 time points, with exposures to dirty water and mold associated with increased risk of multiple allergic symptoms. The mean stress score of study participants at T1 was 8.0 ± 2.1 and at T2, 5.1 ± 3.2, on a 0-10 scale. Participants who experienced specific flood-related exposures reported higher stress scores when compared with their counterparts, especially 1 year after Harvey. Also, a supplementary paired-samples analysis showed that reports of wheezing, shortness of breath, and skin rash did not change between T1 and T2, though other conditions were less commonly reported at T2.

CONCLUSION

These initial Houston-3H findings demonstrate that flooding experiences that occurred as a consequence of Hurricane Harvey had lasting impacts on the health of Houstonians up to 1 year after the hurricane.

Microbiome for Mars: surveying microbiome connections to healthcare with implications for long-duration human spaceflight, virtual workshop, July 13, 2020

The inaugural "Microbiome for Mars" virtual workshop took place on July 13, 2020. This event assembled leaders in microbiome research and development to discuss their work and how it may relate to long-duration human space travel. The conference focused on surveying current microbiome research, future endeavors, and how this growing field could broadly impact human health and space exploration. This report summarizes each speaker's presentation in the order presented at the workshop.

Oligonucleotide Capture Sequencing of the SARS-CoV-2 Genome and Subgenomic Fragments from COVID-19 Individuals

The newly emerged and rapidly spreading SARS-CoV-2 causes coronavirus disease 2019 (COVID-19). To facilitate a deeper understanding of the viral biology we developed a capture sequencing methodology to generate SARS-CoV-2 genomic and transcriptome sequences from infected patients. We utilized an oligonucleotide probe-set representing the full-length genome to obtain both genomic and transcriptome (subgenomic open reading frames [ORFs]) sequences from 45 SARS-CoV-2 clinical samples with varying viral titers. For samples with higher viral loads (cycle threshold value under 33, based on the CDC qPCR assay) complete genomes were generated. Analysis of junction reads revealed regions of differential transcriptional activity and provided evidence of expression of ORF10. Heterogeneous allelic frequencies along the 20kb ORF1ab gene suggested the presence of a defective interfering viral RNA species subpopulation in one sample. The associated workflow is straightforward, and hybridization-based capture offers an effective and scalable approach for sequencing SARS-CoV-2 from patient samples.

Baseline Oral Microbiome and All-cancer Incidence in a Cohort of Nonsmoking Mexican American Women

Given the increasing evidence that the oral microbiome is involved in obesity, diabetes, and cancer risk, we investigated baseline oral microbiota profiles in relation to all-cancer incidence among nonsmoking women enrolled in a Texas cohort of first- and second-generation immigrants of Mexican origin. We characterized the 16Sv4 rDNA microbiome in oral mouthwash samples collected at baseline from a representative subset of 305 nonsmoking women, ages 20-75 years. We evaluated within- (alpha) and between-sample (beta) diversity by incident cancer status and applied linear discriminant analysis (LDA) effect size analysis to assess differentially abundant taxa. Diversity and candidate taxa in relation to all-cancer incidence were evaluated in multivariable-adjusted Cox regression models. Over 8.8 median years of follow-up, 31 incident cancer cases were identified and verified. Advanced age, greater acculturation, and cardiometabolic risk factors were associated with all-cancer incidence. Higher alpha diversity (age-adjusted P _difference < 0.01) and distinct biological communities (P _difference = 0.002) were observed by incident cancer status. Each unit increase in the Shannon diversity index yielded >8-fold increase in all-cancer and obesity-related cancer risk [multivariable-adjusted HR (95% confidence interval), 8.11 (3.14-20.94) and 10.72 (3.30-34.84), respectively] with similar findings for the inverse Simpson index. Streptococcus was enriched among women who did not develop cancer, while Fusobacterium, Prevotella, Mogibacterium, Campylobacter, Lachnoanaerobaculum, Dialister, and Atopobium were higher among women who developed cancer (LDA score ≥ 3; q-value < 0.01). This initial study of oral microbiota and overall cancer risk in nonsmoking Mexican American women suggests the readily accessible oral microbiota as a promising biomarker. PREVENTION RELEVANCE: Mexican American women suffer a disproportionate burden of chronic health conditions that increase cancer risk. Few investigations of the microbiome, a key determinant of host health, have been conducted among this group. Oral microbiota profiles may provide early and accessible cancer biomarker data on invasive bacteria or community disruptions.

Effect of Antibiotics on Gut and Vaginal Microbiomes Associated with Cervical Cancer Development in Mice

Antibiotics affect microbial diversity in the gut, leading to dysbiosis and impaired immunity. However, the impact of antibiotics on microbial communities at other sites, such as vagina is less understood. It is also not clear whether changes induced by antibiotics in both microbiomes affect the development of cervical cancer. In this study, we utilized the murine model to evaluate these questions. We show that oral application of broad-spectrum antibiotics in mice changed not only diversity, but composition and sharing of gut and vaginal microbiomes in mice and influenced cervical cancer development in an orthotopic tumor model. Antibiotics decreased richness and diversity indexes in the gut but increased them in the vagina. Some beneficial taxa, such as Bacteroides, Ruminococcaceae, and Lachnospiraceae increased their abundance in the vagina while other pathogenic species, such as Proteobacteria, were decreased. As a result of the changes, mice with greater richness and diversity of the vaginal microbiome after antibiotics exposure were less likely developed tumors. No association between richness and diversity of the gut microbiome and tumor development was identified.

Metagenomic Information Recovery from Human Stool Samples Is Influenced by Sequencing Depth and Profiling Method

Sequencing of the 16S rRNA gene (16S) has long been a go-to method for microbiome characterization due to its accessibility and lower cost compared to shotgun metagenomic sequencing (SMS). However, 16S sequencing rarely provides species-level resolution and cannot provide direct assessment of other taxa (e.g., viruses and fungi) or functional gene content. Shallow shotgun metagenomic sequencing (SSMS) has emerged as an approach to bridge the gap between 16S sequencing and deep metagenomic sequencing. SSMS is cost-competitive with 16S sequencing, while also providing species-level resolution and functional gene content insights. In the present study, we evaluated the effects of sequencing depth on marker gene-mapping- and alignment-based annotation of bacteria in healthy human stool samples. The number of identified taxa decreased with lower sequencing depths, particularly with the marker gene-mapping-based approach. Other annotations, including viruses and pathways, also showed a depth-dependent effect on feature recovery. These results refine the understanding of the suitability and shortcomings of SSMS, as well as annotation tools for metagenomic analyses in human stool samples. Results may also translate to other sample types and may open the opportunity to explore the effect of sequencing depth and annotation method.

Dietary inflammatory potential in relation to the gut microbiome: results from a cross-sectional study

Diet has direct and indirect effects on health through inflammation and the gut microbiome. We investigated total dietary inflammatory potential via the literature-derived index (Dietary Inflammatory Index (DII®)) with gut microbiota diversity, composition and function. In cancer-free patient volunteers initially approached at colonoscopy and healthy volunteers recruited from the medical centre community, we assessed 16S ribosomal DNA in all subjects who provided dietary assessments and stool samples (n 101) and the gut metagenome in a subset of patients with residual fasting blood samples (n 34). Associations of energy-adjusted DII scores with microbial diversity and composition were examined using linear regression, permutational multivariate ANOVA and linear discriminant analysis. Spearman correlation was used to evaluate associations of species and pathways with DII and circulating inflammatory markers. Across DII levels, α- and β-diversity did not significantly differ; however, Ruminococcus torques, Eubacterium nodatum, Acidaminococcus intestini and Clostridium leptum were more abundant in the most pro-inflammatory diet group, while Akkermansia muciniphila was enriched in the most anti-inflammatory diet group. With adjustment for age and BMI, R. torques, E. nodatum and A. intestini remained significantly associated with a more pro-inflammatory diet. In the metagenomic and fasting blood subset, A. intestini was correlated with circulating plasminogen activator inhibitor-1, a pro-inflammatory marker (rho = 0·40), but no associations remained significant upon correction for multiple testing. An index reflecting overall inflammatory potential of the diet was associated with specific microbes, but not overall diversity of the gut microbiome in our study. Findings from this preliminary study warrant further research in larger samples and prospective cohorts.

Rationale, design and baseline characteristics of the Microbiome and Insulin Longitudinal Evaluation Study (MILES)

AIM

To investigate the role of the gut microbiome in regulating key insulin homeostasis traits (insulin sensitivity, insulin secretion and insulin clearance) whose dysfunction leads to type 2 diabetes (T2D).

MATERIALS AND METHODS

The Microbiome and Insulin Longitudinal Evaluation Study (MILES) focuses on African American and non-Hispanic white participants aged 40-80 years without diabetes. Three study visits are planned (at baseline, 15 and 30 months). Baseline measurements include assessment of the stool microbiome and administration of an oral glucose tolerance test, which will yield indexes of insulin sensitivity, insulin secretion and insulin clearance. The gut microbiome profile (composition and function) will be determined using whole metagenome shotgun sequencing along with analyses of plasma short chain fatty acids. Additional data collected include dietary history, sociodemographic factors, health habits, anthropometry, medical history, medications and family history. Most assessments are repeated 15 and 30 months following baseline.

RESULTS

After screening 875 individuals, 129 African American and 224 non-Hispanic white participants were enrolled. At baseline, African American participants have higher blood pressure, weight, body mass index, waist and hip circumferences but similar waist-hip ratio compared with the non-Hispanic white participants. On average, African American participants are less insulin-sensitive and have higher acute insulin secretion and lower insulin clearance.

CONCLUSIONS

The longitudinal design and robust characterization of potential mediators will allow for the assessment of glucose and insulin homeostasis and gut microbiota as they change over time, improving our ability to discern causal relationships between the microbiome and the insulin homeostasis traits whose deterioration determines T2D, setting the stage for future microbiome-directed therapies to prevent and treat T2D.

Longitudinal Changes in Early Nasal Microbiota and the Risk of Childhood Asthma

OBJECTIVES

Although the airway microbiota is a highly dynamic ecology, the role of longitudinal changes in airway microbiota during early childhood in asthma development is unclear. We aimed to investigate the association of longitudinal changes in early nasal microbiota with the risk of developing asthma.

METHODS

In this prospective, population-based birth cohort study, we followed children from birth to age 7 years. The nasal microbiota was tested by using 16S ribosomal RNA gene sequencing at ages 2, 13, and 24 months. We applied an unsupervised machine learning approach to identify longitudinal nasal microbiota profiles during age 2 to 13 months (the primary exposure) and during age 2 to 24 months (the secondary exposure) and examined the association of these profiles with the risk of physician-diagnosed asthma at age 7 years.

RESULTS

Of the analytic cohort of 704 children, 57 (8%) later developed asthma. We identified 4 distinct longitudinal nasal microbiota profiles during age 2 to 13 months. In the multivariable analysis, compared with the persistent Moraxella dominance profile during age 2 to 13 months, the persistent Moraxella sparsity profile was associated with a significantly higher risk of asthma (adjusted odds ratio, 2.74; 95% confidence interval, 1.20-6.27). Similar associations were observed between the longitudinal changes in nasal microbiota during age 2 to 24 months and risk of asthma.

CONCLUSIONS

Children with an altered longitudinal pattern in the nasal microbiota during early childhood had a high risk of developing asthma. Our data guide the development of primary prevention strategies (eg, early identification of children at high risk and modification of microbiota) for childhood asthma. These observations present a new avenue for risk modification for asthma (eg, microbiota modification).

Alcohol use alters the colonic mucosa-associated gut microbiota in humans

Alcohol misuse is a risk factor for many adverse health outcomes. Alcohol misuse has been associated with an imbalance of gut microbiota in preclinical models and alcoholic diseases. We hypothesized that daily alcohol use would change the community composition and structure of the human colonic gut microbiota. Thirty-four polyp-free individuals donated 97 snap-frozen colonic biopsies. Microbial DNA was sequenced for the 16S ribosomal RNA gene hypervariable region 4. The SILVA database was used for operational taxonomic unit classification. Alcohol use was assessed using a food frequency questionnaire. We compared the biodiversity and relative abundance of the taxa among never drinkers (ND, n = 9), former drinkers (FD, n = 10), current light drinkers (LD, <2 drinks daily, n = 9), and current heavy drinkers (HD, ≥2 drinks daily, n = 6). False discovery rate-adjusted P values (q values) < .05 indicated statistical significance. HD had the lowest α diversity (Shannon index q value < 0.001), and HD's microbial composition differed the most from the other groups (P value = .002). LD had the highest relative abundance of Akkermansia (q values < 0.001). HD had the lowest relative abundance of Subdoligranulum, Roseburia, and Lachnospiraceaeunc91005 but the highest relative abundance of Lachnospiraceaeunc8895 (all q values < 0.05). The multivariable negative binomial regression model supported these observations. ND and FD had a similar microbial profile. Heavy alcohol use was associated with impaired gut microbiota that may partially mediate its effect on health outcomes.

Rotavirus infection induces glycan availability to promote ileum-specific changes in the microbiome aiding rotavirus virulence

Multiple studies have identified changes within the gut microbiome in response to diarrheal-inducing bacterial pathogens. However, examination of the microbiome in response to viral pathogens remains understudied. Compounding this, many studies use fecal samples to assess microbiome composition; which may not accurately mirror changes within the small intestine, the primary site for most enteric virus infections. As a result, the functional significance of small intestinal microbiome shifts during infection is not well defined. To address these gaps, rotavirus-infected neonatal mice were examined for changes in bacterial community dynamics, host gene expression, and tissue recovery during infection. Profiling bacterial communities using 16S rRNA sequencing suggested significant and distinct changes in ileal communities in response to rotavirus infection, with no significant changes for other gastrointestinal (GI) compartments. At 1-d post-infection, we observed a loss in Lactobacillus species from the ileum, but an increase in Bacteroides and Akkermansia, both of which exhibit mucin-digesting capabilities. Concomitant with the bacterial community shifts, we observed a loss of mucin-filled goblet cells in the small intestine at d 1, with recovery occurring by d 3. Rotavirus infection of mucin-producing cell lines and human intestinal enteroids (HIEs) stimulated release of stored mucin granules, similar to in vivo findings. In vitro, incubation of mucins with Bacteroides or Akkermansia members resulted in significant glycan degradation, which altered the binding capacity of rotavirus in silico and in vitro. Taken together, these data suggest that the response to and recovery from rotavirus-diarrhea is unique between sub-compartments of the GI tract and may be influenced by mucin-degrading microbes.

Exploration of the Vaginal and Gut Microbiome in African American Women by Body Mass Index, Class of Obesity, and Gestational Weight Gain: A Pilot Study

OBJECTIVE

This study determines the differences in the distal gut and vaginal microbiome in African American (AA) women by prepregnancy body mass index and gestational weight gain (GWG) comparing women with and without obesity and by obesity class.

STUDY DESIGN

We prospectively sampled the vaginal and distal gut microbiome in pregnant AA women at two time points during pregnancy. Samples were analyzed using high-throughput sequencing of the V4 region of the 16S ribosomal ribonucleic acid gene.

RESULTS

Distinct differences in vaginal and distal gut α-diversity were observed at time point 1 between women with and without obesity by total GWG. Significant differences in distal gut β-diversity were also found at time point 1 in obese women by GWG. Within the Bacteroides genus, a significant association was observed by total GWG among obese women which was absent in nonobese women. Women with class III obesity who experienced low GWG had the lowest abundance of distal gut Bacteroides and appreciably higher relative abundance of a consortia of vaginal taxa including Atopobium, Gardnerella, Prevotella, and Sneathia.

CONCLUSION

These results contribute new evidence showing that GWG in combination with obesity and obesity class is associated with an altered distal gut and vaginal composition early in pregnancy among AA women.

Fungal cutaneous microbiome and host determinants in preterm and term neonates

BACKGROUND

The neonatal cutaneous mycobiome has not been characterized in preterm infants. Invasive fungal infections in preterm neonates are associated with high mortality. The immaturity of the preterm skin predisposes neonates to invasive infection by skin colonizers. We report the clinical and host determinants that influence the skin mycobiome.

METHODS

Skin swabs from the antecubital fossa, forehead, and gluteal region of 15 preterm and 15 term neonates were obtained during the first 5 weeks of life. The mycobiome was sequenced using the conserved pan-fungal ITS2 region. Blood samples were used to genotype immune modulating genes. Clinical metadata was collected to determine the clinical predictors of the abundance and diversity of the skin mycobiome.

RESULTS

The neonatal mycobiome is characterized by few taxa. Alpha diversity of the mycobiome is influenced by antibiotic exposure, the forehead body site, and the neonatal intensive care unit (NICU) environment. Beta diversity varies with mode of delivery, diet, and body site. The host determinants of the cutaneous microbiome include single-nucleotide polymorphisms in TLR4, NLRP3,CARD8, and NOD2.

CONCLUSION

The neonatal cutaneous mycobiome is composed of few genera and is influenced by clinical factors and host genetics, the understanding of which will inform preventive strategies against invasive fungal infections.

Metagenomics of the faecal virome indicate a cumulative effect of enterovirus and gluten amount on the risk of coeliac disease autoimmunity in genetically at risk children: the TEDDY study

OBJECTIVE

Higher gluten intake, frequent gastrointestinal infections and adenovirus, enterovirus, rotavirus and reovirus have been proposed as environmental triggers for coeliac disease. However, it is not known whether an interaction exists between the ingested gluten amount and viral exposures in the development of coeliac disease. This study investigated whether distinct viral exposures alone or together with gluten increase the risk of coeliac disease autoimmunity (CDA) in genetically predisposed children.

DESIGN

The Environmental Determinants of Diabetes in the Young study prospectively followed children carrying the HLA risk haplotypes DQ2 and/or DQ8 and constructed a nested case-control design. From this design, 83 CDA case-control pairs were identified. Median age of CDA was 31 months. Stool samples collected monthly up to the age of 2 years were analysed for virome composition by Illumina next-generation sequencing followed by comprehensive computational virus profiling.

RESULTS

The cumulative number of stool enteroviral exposures between 1 and 2 years of age was associated with an increased risk for CDA. In addition, there was a significant interaction between cumulative stool enteroviral exposures and gluten consumption. The risk conferred by stool enteroviruses was increased in cases reporting higher gluten intake.

CONCLUSIONS

Frequent exposure to enterovirus between 1 and 2 years of age was associated with increased risk of CDA. The increased risk conferred by the interaction between enteroviruses and higher gluten intake indicate a cumulative effect of these factors in the development of CDA.

Oligonucleotide capture sequencing of the SARS-CoV-2 genome and subgenomic fragments from COVID-19 individuals

The newly emerged and rapidly spreading SARS-CoV-2 causes coronavirus disease 2019 (COVID-19). To facilitate a deeper understanding of the viral biology we developed a capture sequencing methodology to generate SARS-CoV-2 genomic and transcriptome sequences from infected patients. We utilized an oligonucleotide probe-set representing the full-length genome to obtain both genomic and transcriptome (subgenomic open reading frames [ORFs]) sequences from 45 SARS-CoV-2 clinical samples with varying viral titers. For samples with higher viral loads (cycle threshold value under 33, based on the CDC qPCR assay) complete genomes were generated. Analysis of junction reads revealed regions of differential transcriptional activity and provided evidence of expression of ORF10. Heterogeneous allelic frequencies along the 20kb ORF1ab gene suggested the presence of a defective interfering viral RNA species subpopulation in one sample. The associated workflow is straightforward, and hybridization-based capture offers an effective and scalable approach for sequencing SARS-CoV-2 from patient samples.

SMRT Sequencing of Paramecium Bursaria Chlorella Virus-1 Reveals Diverse Methylation Stability in Adenines Targeted by Restriction Modification Systems

Chloroviruses (family Phycodnaviridae) infect eukaryotic, freshwater, unicellular green algae. A unique feature of these viruses is an abundance of DNA methyltransferases, with isolates dedicating up to 4.5% of their protein coding potential to these genes. This diversity highlights just one of the long-standing values of the chlorovirus model system; where group-wide epigenomic characterization might begin to elucidate the function(s) of DNA methylation in large dsDNA viruses. We characterized DNA modifications in the prototype chlorovirus, PBCV-1, using single-molecule real time (SMRT) sequencing (aka PacBio). Results were compared to total available sites predicted in silico based on DNA sequence alone. SMRT-software detected N6-methyl-adenine (m6A) at GATC and CATG recognition sites, motifs previously shown to be targeted by PBCV-1 DNA methyltransferases M.CviAI and M. CviAII, respectively. At the same time, PacBio analyses indicated that 10.9% of the PBCV-1 genome had large interpulse duration ratio (ipdRatio) values, the primary metric for DNA modification identification. These events represent 20.6x more sites than can be accounted for by all available adenines in GATC and CATG motifs, suggesting base or backbone modifications other than methylation might be present. To define methylation stability, we cross-compared methylation status of each GATC and CATG sequence in three biological replicates and found ∼81% of sites were stably methylated, while ∼2% consistently lack methylation. The remaining 17% of sites were stochastically methylated. When methylation status was analyzed for both strands of each target, we show that palindromes existed in completely non-methylated states, fully-methylated states, or hemi-methylated states, though GATC sites more often lack methylation than CATG sequences. Given that both sequences are targeted by not just methyltransferases, but by restriction endonucleases that are together encoded by PBCV-1 as virus-originating restriction modification (RM) systems, there is strong selective pressure to modify all target sites. The finding that most instances of non-methylation are associated with hemi-methylation is congruent with observations that hemi-methylated palindromes are resistant to cleavage by restriction endonucleases. However, sites where hemi-methylation is conserved might represent a unique regulatory function for PBCV-1. This study serves as a baseline for future investigation into the epigenomics of chloroviruses and their giant virus relatives.