ABSTRACT
Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing non-homologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between "additive" and "replacing" HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting 9 method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT.
ABSTRACT
The lack of affordable and effective therapeutics against cognitive impairment has promoted research toward alternative approaches to the treatment of neurodegeneration. In recent years, a bidirectional pathway that allows the gut to communicate with the central nervous system has been recognized as the gut-brain axis. Alterations in the gut microbiota, a dynamic population of trillions of microorganisms residing in the gastrointestinal tract, have been implicated in a variety of pathological states, including neurodegenerative disorders such as Alzheimer's disease (AD). However, probiotic treatment as an affordable and accessible adjuvant therapy for the correction of dysbiosis in AD has not been thoroughly explored. Here, we sought to correct the dysbiosis in an AD mouse model with probiotic supplementation, with the intent of exploring its effects on disease progression. Transgenic 3xTg-AD mice were fed a control or a probiotic diet (Lactobacillus plantarum KY1032 and Lactobacillus curvatus HY7601) for 12 weeks, with the latter leading to a significant increase in the relative abundance of Bacteroidetes. Cognitive functions were evaluated via Barnes Maze trials and improvements in memory performance were detected in probiotic-fed AD mice. Neural tissue analysis of the entorhinal cortex and hippocampus of 10-month-old 3xTg-AD mice demonstrated that astrocytic and microglial densities were reduced in AD mice supplemented with a probiotic diet, with changes more pronounced in probiotic-fed female mice. In addition, elevated numbers of neurons in the hippocampus of probiotic-fed 3xTg-AD mice suggested neuroprotection induced by probiotic supplementation. Our results suggest that probiotic supplementation could be effective in delaying or mitigating early stages of neurodegeneration in the 3xTg-AD animal model. It is vital to explore new possibilities for palliative care for neurodegeneration, and probiotic supplementation could provide an inexpensive and easily implemented adjuvant clinical treatment for AD.
ABSTRACT
During a large-scale bacterial culturing effort of biofilms in the vicinity of a rainbow trout aquaculture facility in Idaho, USA, 10 isolates were identified as having pathogen-inhibiting activity and were characterized further. These isolates were shown to be Gram-negative, rod-shaped bacteria belonging to the genus Pseudomonas. Whole-genome comparisons and multi-locus sequence analysis using four housekeeping genes (16S rRNA, gyrA, rpoB and rpoD) showed that these 10 isolates clustered into four distinct species groups. These comparisons also indicated that these isolates were below the established species cutoffs for the genus Pseudomonas. Further phenotypic characterization using API 20NE, API ZYM and Biolog GENIII assays and chemotaxonomic analysis of cellular fatty acids were carried out. Based on the genomic, physiological and chemotaxonomic properties of these isolates, we concluded that these strains composed four novel species of the genus Pseudomonas. The proposed names are as follows: Pseudomonas aphyarum sp. nov. consisting of strains ID233, ID386T and ID387 with ID386T (=DSM 114641T=ATCC TSD-305T) as the type strain; Pseudomonas rubra sp. nov. consisting of strains ID291T, ID609 and ID1025 with ID291T (=DSM 114640T=ATCC TSD-303T) as the type strain; Pseudomonas idahonensis sp. nov. consisting of strains ID357T and ID1048 with ID357T (=DSM 114609T=ATCC TSD-304T) as the type strain; and Pseudomonas fontis sp. nov. consisting of strains ID656T and ID681 with ID656T (=DSM 114610T=ATCC TSD-306T) as the type strain.
Subject(s)
Fatty Acids , Oncorhynchus mykiss , Animals , Fatty Acids/chemistry , RNA, Ribosomal, 16S/genetics , Phylogeny , Sequence Analysis, DNA , DNA, Bacterial/genetics , Bacterial Typing Techniques , Base Composition , Pseudomonas/geneticsABSTRACT
Multisystem inflammatory syndrome in children (MIS-C) is a rare but serious condition that can develop 4-6 weeks after a school age child becomes infected by SARS-CoV-2. To date, in the United States more than 8,862 cases of MIS-C have been identified and 72 deaths have occurred. This syndrome typically affects children between the ages of 5-13; 57% are Hispanic/Latino/Black/non-Hispanic, 61% of patients are males and 100% have either tested positive for SARS-CoV-2 or had direct contact with someone with COVID-19. Unfortunately, diagnosis of MIS-C is difficult, and delayed diagnosis can lead to cardiogenic shock, intensive care admission, and prolonged hospitalization. There is no validated biomarker for the rapid diagnosis of MIS-C. In this study, we used Grating-coupled Fluorescence Plasmonic (GCFP) microarray technology to develop biomarker signatures in pediatric salvia and serum samples from patients with MIS-C in the United States and Colombia. GCFP measures antibody-antigen interactions at individual regions of interest (ROIs) on a gold-coated diffraction grating sensor chip in a sandwich immunoassay to generate a fluorescent signal based on analyte presence within a sample. Using a microarray printer, we designed a first-generation biosensor chip with the capability of capturing 33 different analytes from 80 µ L of sample (saliva or serum). Here, we show potential biomarker signatures in both saliva and serum samples in six patient cohorts. In saliva samples, we noted occasional analyte outliers on the chip within individual samples and were able to compare those samples to 16S RNA microbiome data. These comparisons indicate differences in relative abundance of oral pathogens within those patients. Microsphere Immunoassay (MIA) of immunoglobulin isotypes was also performed on serum samples and revealed MIS-C patients had several COVID antigen-specific immunoglobulins that were significantly higher than other cohorts, thus identifying potential new targets for the second-generation biosensor chip. MIA also identified additional biomarkers for our second-generation chip, verified biomarker signatures generated on the first-generation chip, and aided in second-generation chip optimization. Interestingly, MIS-C samples from the United States had a more diverse and robust signature than the Colombian samples, which was also illustrated in the MIA cytokine data. These observations identify new MIS-C biomarkers and biomarker signatures for each of the cohorts. Ultimately, these tools may represent a potential diagnostic tool for use in the rapid identification of MIS-C.
ABSTRACT
Here, we present a draft genome sequence of Plesiomonas shigelloides MD22D9, isolated from the digestive tract of the North American medicinal leech Macrobdella decora. The gut microbiome of the medicinal leech is hypothesized to be critical for maintaining host fitness. This genome can provide insights into this uncharacterized microbe-host relationship.
ABSTRACT
Periodical cicadas (Hemiptera: Magicicada) have coevolved with obligate bacteriome-inhabiting microbial symbionts, yet little is known about gut microbial symbiont composition or differences in composition among allochronic Magicicada broods (year classes) which emerge parapatrically or allopatrically in the eastern United States. Here, 16S rRNA amplicon sequencing was performed to determine gut bacterial community profiles of three periodical broods, including II (Connecticut and Virginia, 2013), VI (North Carolina, 2017), and X (Maryland, 2021, and an early emerging nymph collected in Ohio, 2017). Results showed similarities among all nymphal gut microbiomes and between morphologically distinct 17-year Magicicada, namely Magicicada septendecim (Broods II and VI) and 17-year Magicicada cassini (Brood X) providing evidence of a core microbiome, distinct from the microbiome of burrow soil inhabited by the nymphs. Generally, phyla Bacteroidetes [Bacteroidota] (> 50% relative abundance), Actinobacteria [Actinomycetota], or Proteobacteria [Pseudomonadota] represented the core. Acidobacteria and genera Cupriavidus, Mesorhizobium, and Delftia were prevalent in nymphs but less frequent in adults. The primary obligate endosymbiont, Sulcia (Bacteroidetes), was dominant amongst core genera detected. Chryseobacterium were common in Broods VI and X. Chitinophaga, Arthrobacter, and Renibacterium were common in Brood X, and Pedobacter were common to nymphs of Broods II and VI. Further taxonomic assignment of unclassified Alphaproteobacteria sequencing reads allowed for detection of multiple copies of the Hodgkinia 16S rRNA gene, distinguishable as separate operational taxonomic units present simultaneously. As major emergences of the broods examined here occur at 17-year intervals, this study will provide a valuable comparative baseline in this era of a changing climate.
Subject(s)
Gastrointestinal Microbiome , Hemiptera , Animals , Bacteria/genetics , Gastrointestinal Microbiome/genetics , Hemiptera/genetics , RNA, Ribosomal, 16S/genetics , Soil , United StatesABSTRACT
Aquaculture is an important tool for solving the growing worldwide food demand, but infectious diseases of farmed animals represent a serious roadblock to continued industry growth. Therefore, it is essential to understand the microbial communities that reside within the built environments of aquaculture facilities to identify reservoirs of bacterial pathogens and potential correlations between commensal species and specific disease agents. Here, we present the results from 3 years of sampling a commercial rainbow trout aquaculture facility. We observed that the microbial communities residing on the abiotic surfaces within the hatchery were distinct from those residing on the surfaces at the facility's water source as well as the production raceways, despite similar communities in the water column at each location. Also, a subset of the water community seeds the biofilm communities. Lastly, we detected a common fish pathogen, Flavobacterium columnare, within the hatchery, including at the source water inlet. Importantly, the relative abundance of this pathogen was correlated with clinical disease. Our results characterized the microbial communities in an aquaculture facility, established that the hatchery environment contains a unique community composition and demonstrated that a specific fish pathogen resides within abiotic surface biofilms and is seeded from the natural water source.
Subject(s)
Fish Diseases , Microbiota , Oncorhynchus mykiss , Animals , Aquaculture , Fish Diseases/microbiology , Flavobacterium/genetics , Fresh Water/microbiology , WaterABSTRACT
Nonalcoholic fatty liver disease (NAFLD) is the most common form of liver disease and is often the precursor for more serious liver conditions such as nonalcoholic steatohepatitis and cirrhosis. Although the gut microbiome has been implicated in the development of NAFLD, the strong association of obesity with NAFLD and its effect on microbiome structure has made interpreting study outcomes difficult. In the present study, we examined the taxonomic and functional differences between the microbiomes of youth with obesity and with and without NAFLD. Shotgun metagenome sequencing was performed to profile the microbiomes of 36 subjects, half of whom were diagnosed with NAFLD using abdominal magnetic resonance imaging. Beta diversity analysis showed community-wide differences between the groups (p = 0.002). Specific taxonomic differences included increased relative abundances of the species Fusicatenibacter saccharivorans (p = 0.042), Romboutsia ilealis (p = 0.046), and Actinomyces sp. ICM47 (p = 0.0009), and a decrease of Bacteroides thetaiotamicron (p = 0.0002), in the NAFLD group as compared with the non-NAFLD group. At the phylum level, Bacteroidetes (p < 0.0001) was decreased in the NAFLD group. Functionally, branched-chain amino acid (p = 0.01343) and aromatic amino acid (p = 0.01343) synthesis pathways had increased relative abundances in the NAFLD group along with numerous energy use pathways, including pyruvate fermentation to acetate (p = 0.01318). Conclusion: Community-wide differences were noted based on NAFLD status, and individual bacterial species along with specific metabolic pathways were identified as potential drivers of these differences. The results of the present study support the idea that the NAFLD phenotype displays a differentiated microbial and functional signature from the obesity phenotype.
Subject(s)
Gastrointestinal Microbiome , Non-alcoholic Fatty Liver Disease , Bacteria/genetics , Gastrointestinal Microbiome/genetics , Humans , Metabolic Networks and Pathways/genetics , Metagenomics , Non-alcoholic Fatty Liver Disease/genetics , Obesity/complicationsABSTRACT
Some of the protist species which colonize the hindguts of wood-feeding Reticulitermes termites are associated with endosymbiotic bacteria belonging to the genus Endomicrobium. In this study, we focused on the endosymbionts of three protist species from Reticulitermes flavipes, as follows: Pyrsonympha vertens, Trichonympha agilis, and Dinenympha species II. Since these protist hosts represented members of different taxa which colonize separate niches within the hindguts of their termite hosts, we investigated if these differences translated to differential gene content and expression in their endosymbionts. Following assembly and comparative genome and transcriptome analyses, we discovered that these endosymbionts differed with respect to some possible niche-specific traits, such as carbon metabolism. Our analyses suggest that species-specific genes related to carbon metabolism were acquired by horizontal gene transfer (HGT) and may have come from taxa which are common in the termite hind gut. In addition, our analyses suggested that these endosymbionts contain and express genes related to natural transformation (competence) and recombination. Taken together, the presence of genes acquired by HGT and a putative competence pathway suggest that these endosymbionts are not cut off from gene flow and that competence may be a mechanism by which members of Endomicrobium can acquire new traits. IMPORTANCE The composition and structure of wood, which contains cellulose, hemicellulose, and lignin, prevent most organisms from using this common food source. Termites are a rare exception among animals, and they rely on a complex microbiota housed in their hindguts to use wood as a source of food. The lower termite, Reticulitermes flavipes, houses a variety of protists and prokaryotes that are the key players in the disassembly of lignocellulose. Here, we describe the genomes and the gene expression profiles of five Endomicrobium endosymbionts living inside three different protist species from R. flavipes. Data from these genomes suggest that these Endomicrobium species have different mechanisms for using carbon. In addition, they harbor genes that may be used to import DNA from their environment. This process of DNA uptake may contribute to the high levels of horizontal gene transfer noted previously in Endomicrobium species.
Subject(s)
Isoptera , Animals , Bacteria , Carbon/metabolism , Eukaryota/genetics , Isoptera/microbiology , Phylogeny , Symbiosis/genetics , TranscriptomeABSTRACT
Diseases caused by the fish pathogens Flavobacterium columnare and Flavobacterium psychrophilum are major contributors of preventable losses in the aquaculture industry. The persistent and difficult-to-control infections caused by these bacteria make timely intervention and prophylactic elimination of pathogen reservoirs important measures to combat these disease-causing agents. In this study, we present two independent assays for detecting these pathogens in a range of environmental samples. Natural water samples were inoculated with F. columnare and F. psychrophilum over 5 orders of magnitude, and pathogen levels were detected using Illumina MiSeq sequencing and droplet digital PCR. Both detection methods accurately identified pathogen-positive samples and showed good agreement in quantifying each pathogen. Additionally, the real-world application of these approaches was demonstrated using environmental samples collected at a rainbow trout (Oncorhynchus mykiss) aquaculture facility. These results show that both methods can serve as useful tools for surveillance efforts in aquaculture facilities, where the early detection of these flavobacterial pathogens may direct preventative measures to reduce disease occurrence. IMPORTANCE Early detection of a deadly disease outbreak in a population can be the difference between mass mortality or mitigated effects. In the present study, we evaluated and compared two molecular techniques for detecting economically impactful aquaculture pathogens. We demonstrate that one of these techniques, 16S rRNA gene sequencing using Illumina MiSeq technology, provides the ability to accurately detect two freshwater fish pathogens, F. columnare and F. psychrophilum, while simultaneously profiling the native microbial community. The second technique, droplet digital PCR, is commonly used for pathogen detection, and the results obtained using the assays we designed with this method served to validate those obtained using the MiSeq method. These two methods offer distinct advantages. The MiSeq method pairs pathogen detection and microbial community profiling to answer immediate and long-term fish health concerns, while the droplet digital PCR method provides fast and highly sensitive detection that is useful for surveillance and rapid clinical responses.
Subject(s)
Fish Diseases , Flavobacteriaceae Infections , Oncorhynchus mykiss , Animals , Aquaculture , Fish Diseases/microbiology , Flavobacteriaceae Infections/diagnosis , Flavobacteriaceae Infections/microbiology , Flavobacteriaceae Infections/veterinary , Flavobacterium/genetics , Oncorhynchus mykiss/microbiology , RNA, Ribosomal, 16S/geneticsABSTRACT
The opportunistic pathogen Pseudomonas aeruginosa expresses a small molecular weight, cysteine-rich protein (PmtA), identified as a metallothionein (MT) protein family member. The MT family proteins have been well-characterized in eukaryotes as essential for zinc and copper homeostasis, protection against oxidative stress, and the ability to modify a variety of immune activities. Bacterial MTs share sequence homology, antioxidant chemistry, and heavy metal-binding capacity with eukaryotic MTs, however, the impact of bacterial MTs on virulence and infection have not been well-studied. In the present study, we investigated the role of PmtA in P. aeruginosa PAO1 using a PmtA-deficient strain (ΔpmtA). Here we demonstrated the virulence factor, pyocyanin, relies on the expression of PmtA. We showed that PmtA may be protective against oxidative stress, as an alternative antioxidant, glutathione, can rescue pyocyanin expression. Furthermore, the expression of phzM, which encodes a pyocyanin precursor enzyme, was decreased in the ΔpmtA mutant during early stationary phase. Upregulated pmtA expression was previously detected in confluent biofilms, which are essential for chronic infection, and we observed that the ΔpmtA mutant was disrupted for biofilm formation. As biofilms also modulate antibiotic susceptibility, we examined the ΔpmtA mutant susceptibility to antibiotics and found that the ΔpmtA mutant is more susceptible to cefepime and ciprofloxacin than the wild-type strain. Finally, we observed that the deletion of pmtA results in decreased virulence in a waxworm model. Taken together, our results support the conclusion that PmtA is necessary for the full virulence of P. aeruginosa and may represent a potential target for therapeutic intervention.
ABSTRACT
A draft genome sequence of Aeromonas popoffii ID682, isolated from a natural water source in Idaho, is presented here. A. popoffii is a relatively understudied species within a diverse, expanding freshwater genus of bacteria. Here, we describe only the second genome published for this species to date.
ABSTRACT
We present a draft genome sequence of Janthinobacterium lividum strain ID1246, isolated from within a rainbow trout hatchery raceway. Janthinobacterium spp. are well-known producers of antimicrobial compounds. Due to the unique isolation source, this genome may yield novel biosynthetic gene clusters.
ABSTRACT
Microbiota perform vital functions for their mammalian hosts, making them potential drivers of host evolution. Understanding effects of environmental factors and host characteristics on the composition and biodiversity of the microbiota may provide novel insights into the origin and maintenance of these symbiotic relationships. Our goals were to (1) characterize biodiversity of oral and rectal microbiota in bats from Puerto Rico; and (2) determine the effects of geographic location and host characteristics on that biodiversity. We collected bats and their microbiota from three sites, and used four metrics (species richness, Shannon diversity, Camargo evenness, Berger-Parker dominance) to characterize biodiversity. We quantified the relative importance of site, host sex, host species-identity, and host foraging-guild on biodiversity of the microbiota. Microbe biodiversity was highly variable among conspecifics. Geographical location exhibited consistent effects, whereas host sex did not. Within each host guild, host species exhibited consistent differences in biodiversity of oral microbiota and of rectal microbiota. Oral microbe biodiversity was indistinguishable between guilds, whereas rectal microbe biodiversity was significantly greater in carnivores than in herbivores. The high intraspecific and spatial variation in microbe biodiversity necessitate a relatively large number of samples to statistically isolate the effects of environmental or host characteristics on the microbiota. Species-specific biodiversity of oral microbiota suggests these communities are structured by direct interactions with the host immune system via epithelial receptors. In contrast, the number of microbial taxa that a host gut supports may be driven by host diet-diversity or composition.
Subject(s)
Chiroptera , Microbiota , Animals , Biodiversity , Diet , Hispanic or Latino , Humans , Puerto RicoABSTRACT
Leeches are found in terrestrial, aquatic, and marine habitats on all continents. Sanguivorous leeches have been used in medicine for millennia. Modern scientific uses include studies of neurons, anticoagulants, and gut microbial symbioses. Hirudo verbana, the European medicinal leech, maintains a gut community dominated by two bacterial symbionts, Aeromonas veronii and Mucinivorans hirudinis, which sometimes account for as much as 97% of the total crop microbiota. The highly simplified gut anatomy and microbiome of H. verbana make it an excellent model organism for studying gut microbial dynamics. The North American medicinal leech, Macrobdella decora, is a hirudinid leech native to Canada and the northern United States. In this study, we show that M. decora symbiont communities are very similar to those in H. verbana. We performed an extensive study using field-caught M. decora and purchased H. verbana from two suppliers. Deep sequencing of the V4 region of the 16S rRNA gene allowed us to determine that the core microbiome of M. decora consists of Bacteroides, Aeromonas, Proteocatella, and Butyricicoccus. The analysis revealed that the compositions of the gut microbiomes of the two leech species were significantly different at all taxonomic levels. The R2 value was highest at the genus and amplicon sequence variant (ASV) levels and much lower at the phylum, class, and order levels. The gut and bladder microbial communities were distinct. We propose that M. decora is an alternative to H. verbana for studies of wild-caught animals and provide evidence for the conservation of digestive-tract and bladder symbionts in annelid models.IMPORTANCE Building evidence implicates the gut microbiome in critical animal functions such as regulating digestion, nutrition, immune regulation, and development. Simplified, phylogenetically diverse models for hypothesis testing are necessary because of the difficulty of assigning causative relationships in complex gut microbiomes. Previous research used Hirudo verbana as a tractable animal model of digestive-tract symbioses. Our data show that Macrobdella decora may work just as well without the drawback of being an endangered organism and with the added advantage of easy access to field-caught specimens. The similarity of the microbial community structures of species from two different continents reveals the highly conserved nature of the microbial symbionts in sanguivorous leeches.
Subject(s)
Gastrointestinal Microbiome , Leeches/microbiology , Animals , Eating , High-Throughput Nucleotide Sequencing , Leeches/genetics , RNA, Ribosomal, 16S , Urinary Bladder/microbiologyABSTRACT
Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the "scale" of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.
Subject(s)
Aeromonas/genetics , Gene Transfer, Horizontal , Genomics/methods , SoftwareABSTRACT
Dysgonomonas spp. are facultative heterotrophs which colonize diverse environments, including the hindgut of the lower termite Reticulitermes flavipes Dysgonomonas genomes are enriched for genes involving oligo- and polysaccharide utilization, enabling modification of a wide array of complex glycans. Here, we report draft genome sequences for Dysgonomonas sp. strains BGC7 and HGC4.
ABSTRACT
Aeromonas veronii strain Hm21 was isolated from the medicinal leech Hirudo verbana and is used for genetic studies. We present here the 4.71-Mbp genome with a 56-kb plasmid and identify the mutations present in strains commonly used for genetic engineering.
ABSTRACT
CONTEXT: Nonalcoholic fatty liver disease (NAFLD) is the most common cause of liver disease, affecting approximately 3 in 10 obese children worldwide. OBJECTIVE: We aimed to investigate the potential relationship between gut microbiota and NAFLD in obese youth, while considering the role of PNPLA3 rs738409, a strong genetic contributor to NAFLD. DESIGN: In this cross-sectional study, participants completed an abdominal magnetic resonance imaging to measure hepatic fat fraction (HFF), oral glucose tolerance test, and PNPLA3 rs738409 genotyping. Fecal samples were collected to analyze the V4 region of the 16S rRNA gene for intestinal bacteria characterization. SETTING: Yale Pediatric Obesity Clinic. PARTICIPANTS: Obese youth (body mass indexâ >95th percentile) with NAFLD (HFFâ ≥5.5%; nâ =â 44) and without NAFLD (HFFâ <5.5%; nâ =â 29). MAIN OUTCOME MEASURE: Shannon-Wiener diversity index values and proportional bacterial abundance by NAFLD status and PNPLA3 genotype. RESULTS: Subjects with NAFLD had decreased bacterial alpha-diversity compared with those without NAFLD (Pâ =â 0.013). Subjects with NAFLD showed a higher Firmicutes to Bacteroidetes (F/B) ratio (Pâ =â 0.019) and lower abundance of Bacteroidetes (Pâ =â 0.010), Prevotella (Pâ =â 0.019), Gemmiger (Pâ =â 0.003), and Oscillospira (Pâ =â 0.036). F/B ratio, Bacteroidetes, Gemmiger, and Oscillospira were associated with HFF when controlling for group variations. We also observed an additive effect on HFF by PNPLA3 rs738409 and Gemmiger, and PNPLA3 rs738409 and Oscillospira. CONCLUSIONS: Obese youth with NAFLD have a different gut microbiota composition than those without NAFLD. These differences were still statistically significant when controlling for factors associated with NAFLD, including PNPLA3 rs738409.