Microbial genomes are available at an ever-increasing pace, as cultivation and sequencing become cheaper and obtaining metagenome-assembled genomes (MAGs) becomes more effective. Phylogenetic ...placement methods to contextualize hundreds of thousands of genomes must thus be efficiently scalable and sensitive from closely related strains to divergent phyla. We present PhyloPhlAn 3.0, an accurate, rapid, and easy-to-use method for large-scale microbial genome characterization and phylogenetic analysis at multiple levels of resolution. PhyloPhlAn 3.0 can assign genomes from isolate sequencing or MAGs to species-level genome bins built from >230,000 publically available sequences. For individual clades of interest, it reconstructs strain-level phylogenies from among the closest species using clade-specific maximally informative markers. At the other extreme of resolution, it scales to large phylogenies comprising >17,000 microbial species. Examples including Staphylococcus aureus isolates, gut metagenomes, and meta-analyses demonstrate the ability of PhyloPhlAn 3.0 to support genomic and metagenomic analyses.
Although microbial communities are associated with human, environmental, plant, and animal health, there exists no cost-effective method for precisely characterizing species and genes in such ...communities. While deep whole-metagenome shotgun (WMS) sequencing provides high taxonomic and functional resolution, it is often prohibitively expensive for large-scale studies. The prevailing alternative, 16S rRNA gene amplicon (16S) sequencing, often does not resolve taxonomy past the genus level and provides only moderately accurate predictions of the functional profile; thus, there is currently no widely accepted approach to affordable, high-resolution, taxonomic, and functional microbiome analysis. To address this technology gap, we evaluated the information content of shallow shotgun sequencing with as low as 0.5 million sequences per sample as an alternative to 16S sequencing for large human microbiome studies. We describe a library preparation protocol enabling shallow shotgun sequencing at approximately the same per-sample cost as 16S sequencing. We analyzed multiple real and simulated biological data sets, including two novel human stool samples with ultradeep sequencing of 2.5 billion sequences per sample, and found that shallow shotgun sequencing recovers more-accurate species-level taxonomic and functional profiles of the human microbiome than 16S sequencing. We discuss the inherent limitations of shallow shotgun sequencing and note that 16S sequencing remains a valuable and important method for taxonomic profiling of novel environments. Although deep WMS sequencing remains the gold standard for high-resolution microbiome analysis, we recommend that researchers consider shallow shotgun sequencing as a useful alternative to 16S sequencing for large-scale human microbiome research studies where WMS sequencing may be cost-prohibitive.
A common refrain in recent microbiome-related academic meetings is that the field needs to move away from broad taxonomic surveys using 16S sequencing and toward more powerful longitudinal studies using shotgun sequencing. However, performing deep shotgun sequencing in large longitudinal studies remains prohibitively expensive for all but the most well-funded research labs and consortia, which leads many researchers to choose 16S sequencing for large studies, followed by deep shotgun sequencing on a subset of targeted samples. Here, we show that shallow- or moderate-depth shotgun sequencing may be used by researchers to obtain species-level taxonomic and functional data at approximately the same cost as amplicon sequencing. While shallow shotgun sequencing is not intended to replace deep shotgun sequencing for strain-level characterization, we recommend that microbiome scientists consider using shallow shotgun sequencing instead of 16S sequencing for large-scale human microbiome studies.
High-throughput amplicon sequencing of large genomic regions remains challenging for short-read technologies. Here, we report a high-throughput amplicon sequencing approach combining unique molecular ...identifiers (UMIs) with Oxford Nanopore Technologies (ONT) or Pacific Biosciences circular consensus sequencing, yielding high-accuracy single-molecule consensus sequences of large genomic regions. We applied our approach to sequence ribosomal RNA operon amplicons (~4,500 bp) and genomic sequences (>10,000 bp) of reference microbial communities in which we observed a chimera rate <0.02%. To reach a mean UMI consensus error rate <0.01%, a UMI read coverage of 15× (ONT R10.3), 25× (ONT R9.4.1) and 3× (Pacific Biosciences circular consensus sequencing) is needed, which provides a mean error rate of 0.0042%, 0.0041% and 0.0007%, respectively.
STING (also known as MITA) mediates the innate antiviral signaling and ubiquitination of STING is key to its function. However, the deubiquitination process of STING is unclear. Here we report that ...USP18 recruits USP20 to deconjugate K48-1inked ubiquitination chains from STING and promotes the stability of STING and the expression of type I IFNs and proinflammatory cytokines after DNA virus infection. USP18 deficiency or knockdown of USP20 resulted in enhanced K48-1inked ubiquitination and accelerated degradation of STING, and impaired activation of IRF3 and NF-κB as well as induction of downstream genes after infection with DNA virus HSV-1 or transfeetion of various DNA ligands. In addition, Uspl8-/- mice were more susceptible to HSV-1 infection compared with the wildtype littermates. USP18 did not deubiquitinate STING in vitro but facilitated USP20 to catalyze deubiquitination of STING in a manner independent of the enzymatic activity of USP18. In addition, reconstitution of STING into Uspl8-/- MEFs restored HSV-1-induced expression of downstream genes and cellular antiviral responses. Our findings thus uncover previously uncharacterized roles of USPI8 and USP20 in mediating virus-triggered signaling and contribute to the understanding of the complicated regulatory system of the innate antiviral responses.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
First pass methods based on BLAST match are commonly used as an initial step to separate the different phylogenetic histories of genes in microbial genomes, and target putative horizontal gene ...transfer (HGT) events. This will continue to be necessary given the rapid growth of genomic data and the technical difficulties in conducting large-scale explicit phylogenetic analyses. However, these methods often produce misleading results due to their inability to resolve indirect phylogenetic links and their vulnerability to stochastic events.
A new computational method of rapid, exhaustive and genome-wide detection of HGT was developed, featuring the systematic analysis of BLAST hit distribution patterns in the context of a priori defined hierarchical evolutionary categories. Genes that fall beyond a series of statistically determined thresholds are identified as not adhering to the typical vertical history of the organisms in question, but instead having a putative horizontal origin. Tests on simulated genomic data suggest that this approach effectively targets atypically distributed genes that are highly likely to be HGT-derived, and exhibits robust performance compared to conventional BLAST-based approaches. This method was further tested on real genomic datasets, including Rickettsia genomes, and was compared to previous studies. Results show consistency with currently employed categories of HGT prediction methods. In-depth analysis of both simulated and real genomic data suggests that the method is notably insensitive to stochastic events such as gene loss, rate variation and database error, which are common challenges to the current methodology. An automated pipeline was created to implement this approach and was made publicly available at: https://github.com/DittmarLab/HGTector. The program is versatile, easily deployed, has a low requirement for computational resources.
HGTector is an effective tool for initial or standalone large-scale discovery of candidate HGT-derived genes.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Influenza A viruses (IAVs) continuously challenge the poultry industry and human health. Elucidation of the host factors that modulate the IAV lifecycle is vital for developing antiviral drugs and ...vaccines. In this study, we infected A549 cells with IAVs and found that host protein contactin-1 (CNTN1), a member of the immunoglobulin superfamily, enhanced viral replication. Bioinformatic prediction and experimental validation indicated that the expression of CNTN1 was reduced by microRNA-200c (miR-200c) through directly targeting. We further showed that CNTN1-modulated viral replication in A549 cells is dependent on type I interferon signaling. Co-immunoprecipitation experiments revealed that CNTN1 specifically interacts with MAVS and promotes its proteasomal degradation by removing its K63-linked ubiquitination. Moreover, we discovered that the deubiquitinase USP25 is recruited by CNTN1 to catalyze the deubiquitination of K63-linked MAVS. Consequently, the CNTN1-induced degradation cascade of MAVS blocked RIG-I-MAVS-mediated interferon signaling, leading to enhanced viral replication. Taken together, our data reveal novel roles of CNTN1 in the type I interferon pathway and regulatory mechanism of IAV replication.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Accurate microbial identification and abundance estimation are crucial for metagenomics analysis. Various methods for classification of metagenomic data and estimation of taxonomic profiles, broadly ...referred to as metagenomic profilers, have been developed. Nevertheless, benchmarking of metagenomic profilers remains challenging because some tools are designed to report relative sequence abundance while others report relative taxonomic abundance. Here we show how misleading conclusions can be drawn by neglecting this distinction between relative abundance types when benchmarking metagenomic profilers. Moreover, we show compelling evidence that interchanging sequence abundance and taxonomic abundance will influence both per-sample summary statistics and cross-sample comparisons. We suggest that the microbiome research community pay attention to potentially misleading biological conclusions arising from this issue when benchmarking metagenomic profilers, by carefully considering the type of abundance data that were analyzed and interpreted and clearly stating the strategy used for metagenomic profiling.
Influenza A virus (IAV) has evolved various strategies to counteract the innate immune response using different viral proteins. However, the mechanism is not fully elucidated. In this study, we ...identified the PB1 protein of H7N9 virus as a new negative regulator of virus- or poly(I:C)-stimulated IFN induction and specifically interacted with and destabilized MAVS. A subsequent study revealed that PB1 promoted E3 ligase RNF5 to catalyze K27-linked polyubiquitination of MAVS at Lys362 and Lys461. Moreover, we found that PB1 preferentially associated with a selective autophagic receptor neighbor of BRCA1 (NBR1) that recognizes ubiquitinated MAVS and delivers it to autophagosomes for degradation. The degradation cascade mediated by PB1 facilitates H7N9 virus infection by blocking the RIG-I-MAVS-mediated innate signaling pathway. Taken together, these data uncover a negative regulatory mechanism involving the PB1-RNF5-MAVS-NBR1 axis and provide insights into an evasion strategy employed by influenza virus that involves selective autophagy and innate signaling pathways.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The activity and stability of the adapter protein MAVS (also known as VISA, Cardif and IPS-1), which critically mediates cellular antiviral responses, are extensively regulated by ubiquitination. ...However, the process whereby MAVS is deubiquitinated is unclear. Here, we report that the ovarian tumor family deubiquitinase 4 (OTUD4) targets MAVS for deubiquitination. Viral infection leads to the IRF3/7-dependent upregulation of OTUD4 which interacts with MAVS to remove K48-linked polyubiquitin chains, thereby maintaining MAVS stability and promoting innate antiviral signaling. Knockout or knockdown of OTUD4 impairs RNA virus-triggered activation of IRF3 and NF-κB, expression of their downstream target genes, and potentiates VSV replication in vitro and in vivo. Consistently, Cre-ER Otud4
or Lyz2-Cre Otud4
mice produce decreased levels of type I interferons and proinflammatory cytokines and exhibit increased sensitivity to VSV infection compared to their control littermates. In addition, reconstitution of MAVS into OTUD4-deficient cells restores virus-induced expression of downstream genes and cellular antiviral responses. Together, our findings uncover an essential role of OTUD4 in virus-triggered signaling and contribute to the understanding of deubiquitination-mediated regulation of innate antiviral responses.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Rapid growth of genome data provides opportunities for updating microbial evolutionary relationships, but this is challenged by the discordant evolution of individual genes. Here we build a reference ...phylogeny of 10,575 evenly-sampled bacterial and archaeal genomes, based on a comprehensive set of 381 markers, using multiple strategies. Our trees indicate remarkably closer evolutionary proximity between Archaea and Bacteria than previous estimates that were limited to fewer "core" genes, such as the ribosomal proteins. The robustness of the results was tested with respect to several variables, including taxon and site sampling, amino acid substitution heterogeneity and saturation, non-vertical evolution, and the impact of exclusion of candidate phyla radiation (CPR) taxa. Our results provide an updated view of domain-level relationships.