Presence and Persistence of ESKAPEE Bacteria before and after Hospital Wastewater Treatment.

The metagenomic surveillance of antimicrobial resistance in wastewater has been suggested as a methodological tool to characterize the distribution, status, and trends of antibiotic-resistant bacteria. In this study, a cross-sectional collection of samples of hospital-associated raw and treated wastewater were obtained from February to March 2020. Shotgun metagenomic sequencing and bioinformatic analysis were performed to characterize bacterial abundance and antimicrobial resistance gene analysis. The main bacterial phyla found in all the samples were as follows: Proteobacteria, Bacteroides, Firmicutes, and Actinobacteria. At the species level, ESKAPEE bacteria such as E. coli relative abundance decreased between raw and treated wastewater, but S. aureus, A. baumannii, and P. aeruginosa increased, as did the persistence of K. pneumoniae in both raw and treated wastewater. A total of 172 different ARGs were detected; blaOXA, blaVEB, blaKPC, blaGES, mphE, mef, erm, msrE, AAC(6'), ant(3″), aadS, lnu, PBP-2, dfrA, vanA-G, tet, and sul were found at the highest abundance and persistence. This study demonstrates the ability of ESKAPEE bacteria to survive tertiary treatment processes of hospital wastewater, as well as the persistence of clinically important antimicrobial resistance genes that are spreading in the environment.

Characterization of immunoglobulin loci in the gigantic genome of Ambystoma mexicanum.

Background: The axolotl, Ambystoma mexicanum is a unique biological model for complete tissue regeneration. Is a neotenic endangered species and is highly susceptible to environmental stress, including infectious disease. In contrast to other amphibians, the axolotl is particularly vulnerable to certain viral infections. Like other salamanders, the axolotl genome is one of the largest (32 Gb) and the impact of genome size on Ig loci architecture is unknown. To better understand the immune response in axolotl, we aimed to characterize the immunoglobulin loci of A. mexicanum and compare it with other model vertebrates. Methods: The most recently published genome sequence of A. mexicanum (V6) was used for alignment-based annotation and manual curation using previously described axolotl Ig sequences or reference sequences from other vertebrates. Gene models were further curated using A. mexicanum spleen RNA-seq data. Human, Xenopus tropicalis, Danio rerio (zebrafish), and eight tetrapod reference genomes were used for comparison. Results: Canonical A. mexicanum heavy chain (IGH), lambda (IGL), sigma (IGS), and the putative surrogate light chain (SLC) loci were identified. No kappa locus was found. More than half of the IGHV genes and the IGHF gene are pseudogenes and there is no clan I IGHV genes. Although the IGH locus size is proportional to genome size, we found local size restriction in the IGHM gene and the V gene intergenic distances. In addition, there were V genes with abnormally large V-intron sizes, which correlated with loss of gene functionality. Conclusion: The A. mexicanum immunoglobulin loci share the same general genome architecture as most studied tetrapods. Consistent with its large genome, Ig loci are larger; however, local size restrictions indicate evolutionary constraints likely to be imposed by high transcriptional demand of certain Ig genes, as well as the V(D)J recombination over very long genomic distance ranges. The A. mexicanum has undergone an extensive process of Ig gene loss which partially explains a reduced potential repertoire diversity that may contribute to its impaired antibody response.

Association of Gut Microbiota with Dietary-dependent Childhood Obesity.

AIM: To evaluate the taxonomic profile of the gut microbiota using metagenomics and the association with diet-dependent childhood obesity. METHODS: A cross-sectional study of a subsample of 46 children was conducted. The children were classified as normal-weight, overweight, and obese according to their age and sex and the World Health Organization (WHO) guidelines. Dietary patterns were determined through principal component analysis. The profile of the human gut microbiota was determined by bioinformatic analysis using whole metagenome shotgun sequencing. The association of gut microbiota and z-BMI, waist circumference and hip circumference, and the possible modifying effect of diet were analyzed using multiple regression models. RESULTS: Children with an abundance of Holdemania spp. and high protein and complex carbohydrate consumption had a lower z-BMI (ß -19.06, p = 0.011), waist circumference (ß -171.92, p = 0.003), and hip circumference (ß -157.57, p = 0.004). In contrast, observed a positive association between Coprococcus catus and the low intake of this dietary pattern with hip circumference (ß 147.87, p = 0.025). Furthermore, the presence of Bilophila spp. and Paraprevotella xylaniphila with high saturated fat and simple carbohydrate consumption we observed a positive association between z-BMI (ß 47.5, p = 0.002), hip circumference (ß 44.54, p = 0.025), and waist circumference (ß 44.34, p = 0.004). CONCLUSION: We suggest that the synergism between diet and the profile of children's gut microbiota can be a factor that could be associated with the development of obesity and its complications in childhood.

Author Correction: Analysis of sequencing strategies and tools for taxonomic annotation: Defining standards for progressive metagenomics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Analysis of sequencing strategies and tools for taxonomic annotation: Defining standards for progressive metagenomics.

Metagenomics research has recently thrived due to DNA sequencing technologies improvement, driving the emergence of new analysis tools and the growth of taxonomic databases. However, there is no all-purpose strategy that can guarantee the best result for a given project and there are several combinations of software, parameters and databases that can be tested. Therefore, we performed an impartial comparison, using statistical measures of classification for eight bioinformatic tools and four taxonomic databases, defining a benchmark framework to evaluate each tool in a standardized context. Using in silico simulated data for 16S rRNA amplicons and whole metagenome shotgun data, we compared the results from different software and database combinations to detect biases related to algorithms or database annotation. Using our benchmark framework, researchers can define cut-off values to evaluate the expected error rate and coverage for their results, regardless the score used by each software. A quick guide to select the best tool, all datasets and scripts to reproduce our results and benchmark any new method are available at https://github.com/Ales-ibt/Metagenomic-benchmark . Finally, we stress out the importance of gold standards, database curation and manual inspection of taxonomic profiling results, for a better and more accurate microbial diversity description.

Lower IgG somatic hypermutation rates during acute dengue virus infection is compatible with a germinal center-independent B cell response.

BACKGROUND: The study of human B cell response to dengue virus (DENV) infection is critical to understand serotype-specific protection and the cross-reactive sub-neutralizing response. Whereas the first is beneficial and thus represents the ultimate goal of vaccination, the latter has been implicated in the development of severe disease, which occurs in a small, albeit significant, fraction of secondary DENV infections. Both primary and secondary infections are associated with the production of poly-reactive and cross-reactive IgG antibodies. METHODS: To gain insight into the effect of DENV infection on the B cell repertoire, we used VH region high-throughput cDNA sequencing of the peripheral blood IgG B cell compartment of 19 individuals during the acute phase of infection. For 11 individuals, a second sample obtained 6 months later was analyzed for comparison. Probabilities of sequencing antibody secreting cells or memory B cells were estimated using second-order Monte Carlo simulation. RESULTS: We found that in acute disease there is an increase in IgG B cell diversity and changes in the relative use of segments IGHV1-2, IGHV1-18, and IGHV1-69. Somewhat unexpectedly, an overall low proportion of somatic hypermutated antibody genes was observed during the acute phase plasmablasts, particularly in secondary infections and those cases with more severe disease. CONCLUSIONS: Our data are consistent with an innate-like antiviral recognition system mediated by B cells using defined germ-line coded B cell receptors, which could provide a rapid germinal center-independent antibody response during the early phase of infection. A model describing concurrent T-dependent and T-independent B cell responses in the context of DENV infection is proposed, which incorporates the selection of B cells using hypomutated IGHV segments and their potential role in poly/cross-reactivity. Its formal demonstration could lead to a definition of its potential implication in antibody-dependent enhancement, and may contribute to rational vaccine development efforts.

Longitudinal analysis of the peripheral B cell repertoire reveals unique effects of immunization with a new influenza virus strain.

BACKGROUND: Despite the potential to produce antibodies that can neutralize different virus (heterotypic neutralization), there is no knowledge of why vaccination against influenza induces protection predominantly against the utilized viral strains (homotypic response). Identification of structural patterns of the B cell repertoire associated to heterotypic neutralization may contribute to identify relevant epitopes for a universal vaccine against influenza. METHODS: Blood samples were collected from volunteers immunized with 2008/2009 trivalent inactivated vaccine (TIV), pandemic H1N1 (pdmH1N1) monovalent inactivated vaccine (MIV) and the 2014/2015 TIV. Neutralization was assessed by hemagglutination and microneutralization test. IgG V(H) amplicons derived from peripheral blood RNA from pre-immune and 7 days post vaccination were subjected to 454-Roche sequencing. Full reconstruction of the sampled repertoires was done with ImmunediveRsity. RESULTS: The TIV induced a predominantly homotypic neutralizing serologic response, while the 09 MIV induced a heterotypic neutralizing seroconversion in 17% of the individuals. Both the 08/09 and the 14/15 TIV were associated with a reduction in clonotypic diversity, whereas 09 MIV was the opposite. Moreover, TIV and MIV induced distinctive patterns of IGHV segment use that are consistent with B cell selection by conserved antigenic determinants shared by the pre-pandemic and the pandemic strains. However, low somatic hypermutation rates in IgG after 09 MIV immunization, but not after 08/09 and 14/15 TIV immunization were observed. Furthermore, no evidence of the original antigenic sin was found in the same individuals after vaccination with the three vaccines. CONCLUSIONS: Immunization with a new influenza virus strain (2009 pdmH1N1) induced unique effects in the peripheral B cell repertoire clonal structure, a stereotyped response involving distinctive IGHV segment use and low somatic hypermutation levels. These parameters were contrastingly different to those observed in response to pre-pandemic and post-pandemic vaccination, and may be the result of clonal selection of common antigenic determinants, as well as germinal center-independent responses that wane as the pandemic strain becomes seasonal. Our findings may contribute in the understanding of the structural and cellular basis required to develop a universal influenza vaccine.

Reconstructing and mining the B cell repertoire with ImmunediveRsity.

The B cell antigen receptor repertoire is highly diverse and constantly modified by clonal selection. High-throughput DNA sequencing (HTS) of the lymphocyte repertoire (Rep-Seq) represents a promising technology to explore such diversity ex-vivo and assist in the identification of antigen-specific antibodies based on molecular signatures of clonal selection. Therefore, integrative tools for repertoire reconstruction and analysis from antibody sequences are needed. We developed ImmunediveRity, a stand-alone pipeline primarily based in R programming for the integral analysis of B cell repertoire data generated by HTS. The pipeline integrates GNU software and in house scripts to perform quality filtering, sequencing noise correction and repertoire reconstruction based on V, D and J segment assignment, clonal origin and unique heavy chain identification. Post-analysis scripts generate a wealth of repertoire metrics that in conjunction with a rich graphical output facilitates sample comparison and repertoire mining. Its performance was tested with raw and curated human and mouse 454-Roche sequencing benchmarks providing good approximations of repertoire structure. Furthermore, ImmunediveRsity was used to mine the B cell repertoire of immunized mice with a model antigen, allowing the identification of previously validated antigen-specific antibodies, and revealing different and unexpected clonal diversity patterns in the post-immunization IgM and IgG compartments. Although ImmunediveRsity is similar to other recently developed tools, it offers significant advantages that facilitate repertoire analysis and repertoire mining. ImmunediveRsity is open source and free for academic purposes and it runs on 64 bit GNU/Linux and MacOS. Available at: https://bitbucket.org/ImmunediveRsity/immunediversity/.