The Effect of Bacille Calmette-Guérin Vaccination on the Composition of the Intestinal Microbiome in Neonates From the MIS BAIR Trial

INTRODUCTION

The early-life intestinal microbiome plays an important role in the development and regulation of the immune system. It is unknown whether the administration of vaccines influences the composition of the intestinal microbiome.

OBJECTIVE

To investigate whether Bacille Calmette-Guérin (BCG) vaccine given in the first few days of life influences the abundance of bacterial taxa and metabolic pathways in the intestinal microbiome at 1 week of age.

METHODS

Healthy, term-born neonates were randomized at birth to receive BCG or no vaccine within the first few days of life. Stool samples were collected at 1 week of age from 335 neonates and analyzed using shotgun metagenomic sequencing and functional analyses.

RESULTS

The composition of the intestinal microbiome was different between neonates born by cesarean section (CS) and those born vaginally. Differences in the composition between BCG-vaccinated and BCG-naïve neonates were only minimal. CS-born BCG-vaccinated neonates had a higher abundance of Staphylococcus lugdunensis compared with CS-born BCG-naïve neonates. The latter had a higher abundance of Streptococcus infantis and Trabulsiella guamensis . Vaginally-born BCG-vaccinated neonates had a higher abundance of Clostridiaceae and Streptococcus parasanguinis compared with vaginally-born BCG-naïve neonates, and a lower abundance of Veillonella atypica and Butyricimonas faecalis. Metabolic pathways that were differently abundant between BCG-vaccinated and BCG-naïve neonates were mainly those involved in sugar degradation and nucleotide/nucleoside biosynthesis.

CONCLUSION

BCG given in the first few days of life has little effect on the composition of the intestinal microbiome at 1 week of age but does influence the abundance of certain metabolic pathways.

Haplotype-resolved genome of heterozygous African cassava cultivar TMEB117 (Manihot esculenta)

Cassava (Manihot esculenta Crantz) is a vital tropical root crop providing essential dietary energy to over 800 million people in tropical and subtropical regions. As a climate-resilient crop, its significance grows as the human population expands. However, yield improvement faces challenges from biotic and abiotic stress and limited breeding. Advanced sequencing and assembly techniques enabled the generation of a highly accurate, nearly complete, haplotype-resolved genome of the African cassava cultivar TMEB117. It is the most accurate cassava genome sequence to date with a base-level accuracy of QV > 64, N50 > 35 Mbp, and 98.9% BUSCO completeness. Over 60% of the genome comprises repetitive elements. We predicted over 45,000 gene models for both haplotypes. This achievement offers valuable insights into the heterozygosity genome organization of the cassava genome, with improved accuracy, completeness, and phased genomes. Due to its high susceptibility to African Cassava Mosaic Virus (ACMV) infections compared to other cassava varieties, TMEB117 provides an ideal reference for studying virus resistance mechanisms, including epigenetic variations and smallRNA expressions.

Structural insights into coordinating 5S RNP rotation with ITS2 pre-RNA processing during ribosome formation

The rixosome defined in Schizosaccharomyces pombe and humans performs diverse roles in pre-ribosomal RNA processing and gene silencing. Here, we isolate and describe the conserved rixosome from Chaetomium thermophilum, which consists of two sub-modules, the sphere-like Rix1-Ipi3-Ipi1 and the butterfly-like Las1-Grc3 complex, connected by a flexible linker. The Rix1 complex of the rixosome utilizes Sda1 as landing platform on nucleoplasmic pre-60S particles to wedge between the 5S rRNA tip and L1-stalk, thereby facilitating the 180° rotation of the immature 5S RNP towards its mature conformation. Upon rixosome positioning, the other sub-module with Las1 endonuclease and Grc3 polynucleotide-kinase can reach a strategic position at the pre-60S foot to cleave and 5' phosphorylate the nearby ITS2 pre-rRNA. Finally, inward movement of the L1 stalk permits the flexible Nop53 N-terminus with its AIM motif to become positioned at the base of the L1-stalk to facilitate Mtr4 helicase-exosome participation for completing ITS2 removal. Thus, the rixosome structure elucidates the coordination of two central ribosome biogenesis events, but its role in gene silencing may adapt similar strategies.

Comparative genomic and metabolomic study of three Streptomyces sp. differing in biological activity

The Streptomyces genus is known to produce many specialized metabolites of value for medicine, but the potential of these metabolites in agronomy remains largely unexplored. In this study, we investigated three phylogenetically closely related Streptomyces strains (B5, B91, and B135) isolated from three distinct soil samples in Sudan. Despite belonging to the same species, these strains exhibited different ranges of Phytophthora infestans inhibition. The objective of this work was to identify the active compound(s) responsible for the inhibition of P. infestans and of other plant pathogens by comparing the genomes and metabolomes of the three strains which showed distinct activity patterns: B5 was the strongest inhibitor of oomycetes, B5 and B91 both inhibited most fungi and B135 was the only strain showing antibacterial activity. Our comparative genomic and metabolomic analysis identified borrelidin as the bioactive compound underlying B5's strong anti-oomycete activity and highlighted a few other metabolites as putative candidates underlying the strains' antifungal and antibacterial activities. This study illustrates the power of comparative genomics and metabolomics on phylogenetically closely related strains of differing activities to highlight bioactive compounds that could contribute to new sustainable crop protection strategies.

Biological hydrogen cyanide emission globally impacts the physiology of both HCN-emitting and HCN-perceiving Pseudomonas

IMPORTANCE

Bacteria communicate by exchanging chemical signals, some of which are volatile and can remotely reach other organisms. HCN was one of the first volatiles discovered to severely impact exposed organisms by inhibiting their respiration. Using HCN-deficient mutants in two Pseudomonas strains, we demonstrate that HCN's impact goes beyond the sole inhibition of respiration and affects both emitting and receiving bacteria in a global way, modulating their motility, biofilm formation, and production of antimicrobial compounds. Our data suggest that bacteria could use HCN not only to control their own cellular functions, but also to remotely influence the behavior of other bacteria sharing the same environment. Since HCN emission occurs in both clinically and environmentally relevant Pseudomonas, these findings are important to better understand or even modulate the expression of bacterial traits involved in both virulence of opportunistic pathogens and in biocontrol efficacy of plant-beneficial strains.

Comparative genomics of Mycoplasma feriruminatoris, a fast-growing pathogen of wild Caprinae

Mycoplasma feriruminatoris is a fast-growing Mycoplasma species isolated from wild Caprinae and first described in 2013. M. feriruminatoris isolates have been associated with arthritis, kerato conjunctivitis, pneumonia and septicemia, but were also recovered from apparently healthy animals. To better understand what defines this species, we performed a genomic survey on 14 strains collected from free-ranging or zoo-housed animals between 1987 and 2017, mostly in Europe. The average chromosome size of the M. feriruminatoris strains was 1,040±0,024 kbp, with 24 % G+C and 852±31 CDS. The core genome and pan-genome of the M. feriruminatoris species contained 628 and 1312 protein families, respectively. The M. feriruminatoris strains displayed a relatively closed pan-genome, with many features and putative virulence factors shared with species from the M. mycoides cluster, including the MIB-MIP Ig cleavage system, a repertoire of DUF285 surface proteins and a complete biosynthetic pathway for galactan. M. feriruminatoris genomes were found to be mostly syntenic, although repertoires of mobile genetic elements, including Mycoplasma Integrative and Conjugative Elements, insertion sequences, and a single plasmid varied. Phylogenetic- and gene content analyses confirmed that M. feriruminatoris was closer to the M. mycoides cluster than to the ruminant species M. yeatsii and M. putrefaciens. Ancestral genome reconstruction showed that the emergence of the M. feriruminatoris species was associated with the gain of 17 gene families, some of which encode defence enzymes and surface proteins, and the loss of 25 others, some of which are involved in sugar transport and metabolism. This comparative study suggests that the M. mycoides cluster could be extended to include M. feriruminatoris. We also find evidence that the specific organization and structure of the DnaA boxes around the oriC of M. feriruminatoris may contribute to drive the remarkable fast growth of this minimal bacterium.

Advancing microbiome research with machine learning: key findings from the ML4Microbiome COST action

The rapid development of machine learning (ML) techniques has opened up the data-dense field of microbiome research for novel therapeutic, diagnostic, and prognostic applications targeting a wide range of disorders, which could substantially improve healthcare practices in the era of precision medicine. However, several challenges must be addressed to exploit the benefits of ML in this field fully. In particular, there is a need to establish "gold standard" protocols for conducting ML analysis experiments and improve interactions between microbiome researchers and ML experts. The Machine Learning Techniques in Human Microbiome Studies (ML4Microbiome) COST Action CA18131 is a European network established in 2019 to promote collaboration between discovery-oriented microbiome researchers and data-driven ML experts to optimize and standardize ML approaches for microbiome analysis. This perspective paper presents the key achievements of ML4Microbiome, which include identifying predictive and discriminatory 'omics' features, improving repeatability and comparability, developing automation procedures, and defining priority areas for the novel development of ML methods targeting the microbiome. The insights gained from ML4Microbiome will help to maximize the potential of ML in microbiome research and pave the way for new and improved healthcare practices.

Metagenomics analysis of the neonatal intestinal resistome

Introduction

The intestinal microbiome forms a major reservoir for antibiotic resistance genes (ARGs). Little is known about the neonatal intestinal resistome.

Objective

The objective of this study was to investigate the intestinal resistome and factors that influence the abundance of ARGs in a large cohort of neonates.

Methods

Shotgun metagenomics was used to analyse the resistome in stool samples collected at 1 week of age from 390 healthy, term-born neonates who did not receive antibiotics.

Results

Overall, 913 ARGs belonging to 27 classes were identified. The most abundant ARGs were those conferring resistance to tetracyclines, quaternary ammonium compounds, and macrolide-lincosamide-streptogramin-B. Phylogenetic composition was strongly associated with the resistome composition. Other factors that were associated with the abundance of ARGs were delivery mode, gestational age, birth weight, feeding method, and antibiotics in the last trimester of pregnancy. Sex, ethnicity, probiotic use during pregnancy, and intrapartum antibiotics had little effect on the abundance of ARGs.

Conclusion

Even in the absence of direct antibiotic exposure, the neonatal intestine harbours a high abundance and a variety of ARGs.

Results and lessons learned from the sbv IMPROVER metagenomics diagnostics for inflammatory bowel disease challenge

A growing body of evidence links gut microbiota changes with inflammatory bowel disease (IBD), raising the potential benefit of exploiting metagenomics data for non-invasive IBD diagnostics. The sbv IMPROVER metagenomics diagnosis for inflammatory bowel disease challenge investigated computational metagenomics methods for discriminating IBD and nonIBD subjects. Participants in this challenge were given independent training and test metagenomics data from IBD and nonIBD subjects, which could be wither either raw read data (sub-challenge 1, SC1) or processed Taxonomy- and Function-based profiles (sub-challenge 2, SC2). A total of 81 anonymized submissions were received between September 2019 and March 2020. Most participants' predictions performed better than random predictions in classifying IBD versus nonIBD, Ulcerative Colitis (UC) versus nonIBD, and Crohn's Disease (CD) versus nonIBD. However, discrimination between UC and CD remains challenging, with the classification quality similar to the set of random predictions. We analyzed the class prediction accuracy, the metagenomics features by the teams, and computational methods used. These results will be openly shared with the scientific community to help advance IBD research and illustrate the application of a range of computational methodologies for effective metagenomic classification.

Local conditions matter: Minimal and variable effects of soil disturbance on microbial communities and functions in European vineyards

Soil tillage or herbicide applications are commonly used in agriculture for weed control. These measures may also represent a disturbance for soil microbial communities and their functions. However, the generality of response patterns of microbial communities and functions to disturbance have rarely been studied at large geographical scales. We investigated how a soil disturbance gradient (low, intermediate, high), realized by either tillage or herbicide application, affects diversity and composition of soil bacterial and fungal communities as well as soil functions in vineyards across five European countries. Microbial alpha-diversity metrics responded to soil disturbance sporadically, but inconsistently across countries. Increasing soil disturbance changed soil microbial community composition at the European level. However, the effects of soil disturbance on the variation of microbial communities were smaller compared to the effects of location and soil covariates. Microbial respiration was consistently impaired by soil disturbance, while effects on decomposition of organic substrates were inconsistent and showed positive and negative responses depending on the respective country. Therefore, we conclude that it is difficult to extrapolate results from one locality to others because microbial communities and environmental conditions vary strongly over larger geographical scales.

Deeper genomic insights into tomato CLE genes repertoire identify new active peptides

Background

In eukaryotes, cell-to-cell communication relies on the activity of small signaling peptides. In plant genomes, many hundreds of genes encode for such short peptide signals. However, only few of them are functionally characterized and due to the small gene size and high sequence variability, the comprehensive identification of such peptide-encoded genes is challenging. The CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION-RELATED (CLE) gene family encodes for short peptides that have a role in plant meristem maintenance, vascular patterning and responses to environment. The full repertoire of CLE genes and the role of CLE signaling in tomato (Solanum lycopersicum)- one of the most important crop plants- has not yet been fully studied.

Results

By using a combined approach, we performed a genome-wide identification of CLE genes using the current tomato genome version SL 4.0. We identified 52 SlCLE genes, including 37 new non annotated before. By analyzing publicly available RNAseq datasets we could confirm the expression of 28 new SlCLE genes. We found that SlCLEs are often expressed in a tissue-, organ- or condition-specific manner. Our analysis shows an interesting gene diversification within the SlCLE family that seems to be a result of gene duplication events. Finally, we could show a biological activity of selected SlCLE peptides in the root growth arrest that was SlCLV2-dependent.

Conclusions

Our improved combined approach revealed 37 new SlCLE genes. These findings are crucial for better understanding of the CLE signaling in tomato. Our phylogenetic analysis pinpoints the closest homologs of Arabidopsis CLE genes in tomato genome and can give a hint about the function of newly identified SlCLEs. The strategy described here can be used to identify more precisely additional short genes in plant genomes. Finally, our work suggests that the mechanism of root-active CLE peptide perception is conserved between Arabidopsis and tomato. In conclusion, our work paves the way to further research on the CLE-dependent circuits modulating tomato development and physiological responses.

Analysis of the Hypoxic Response in a Mouse Cortical Collecting Duct-Derived Cell Line Suggests That Esrra Is Partially Involved in Hif1α-Mediated Hypoxia-Inducible Gene Expression in mCCDcl1 Cells

The kidney is strongly dependent on a continuous oxygen supply, and is conversely highly sensitive to hypoxia. Controlled oxygen gradients are essential for renal control of solutes and urine-concentrating mechanisms, which also depend on various hormones including aldosterone. The cortical collecting duct (CCD) is part of the aldosterone-sensitive distal nephron and possesses a key function in fine-tuned distal salt handling. It is well known that aldosterone is consistently decreased upon hypoxia. Furthermore, a recent study reported a hypoxia-dependent down-regulation of sodium currents within CCD cells. We thus investigated the possibility that cells from the cortical collecting duct are responsive to hypoxia, using the mouse cortical collecting duct cell line mCCDcl1 as a model. By analyzing the hypoxia-dependent transcriptome of mCCDcl1 cells, we found a large number of differentially-expressed genes (3086 in total logFC< −1 or >1) following 24 h of hypoxic conditions (0.2% O2). A gene ontology analysis of the differentially-regulated pathways revealed a strong decrease in oxygen-linked processes such as ATP metabolic functions, oxidative phosphorylation, and cellular and aerobic respiration, while pathways associated with hypoxic responses were robustly increased. The most pronounced regulated genes were confirmed by RT-qPCR. The low expression levels of Epas1 under both normoxic and hypoxic conditions suggest that Hif-1α, rather than Hif-2α, mediates the hypoxic response in mCCDcl1 cells. Accordingly, we generated shRNA-mediated Hif-1α knockdown cells and found Hif-1α to be responsible for the hypoxic induction of established hypoxically-induced genes. Interestingly, we could show that following shRNA-mediated knockdown of Esrra, Hif-1α protein levels were unaffected, but the gene expression levels of Egln3 and Serpine1 were significantly reduced, indicating that Esrra might contribute to the hypoxia-mediated expression of these and possibly other genes. Collectively, mCCDcl1 cells display a broad response to hypoxia and represent an adequate cellular model to study additional factors regulating the response to hypoxia.

Substrate recognition and cryo-EM structure of the ribosome-bound TAC toxin of Mycobacterium tuberculosis

Toxins of toxin-antitoxin systems use diverse mechanisms to control bacterial growth. Here, we focus on the deleterious toxin of the atypical tripartite toxin-antitoxin-chaperone (TAC) system of Mycobacterium tuberculosis, whose inhibition requires the concerted action of the antitoxin and its dedicated SecB-like chaperone. We show that the TAC toxin is a bona fide ribonuclease and identify exact cleavage sites in mRNA targets on a transcriptome-wide scale in vivo. mRNA cleavage by the toxin occurs after the second nucleotide of the ribosomal A-site codon during translation, with a strong preference for CCA codons in vivo. Finally, we report the cryo-EM structure of the ribosome-bound TAC toxin in the presence of native M. tuberculosis cspA mRNA, revealing the specific mechanism by which the TAC toxin interacts with the ribosome and the tRNA in the P-site to cleave its mRNA target.

Toxin-antitoxin systems are widespread in bacteria. Here the authors present structures of M. tuberculosis HigB^TAC alone and bound to the ribosome in the presence of native cspA mRNA, shedding light on its mechanism of translation inhibition.

Dedicated chaperones coordinate co-translational regulation of ribosomal protein production with ribosome assembly to preserve proteostasis

The biogenesis of eukaryotic ribosomes involves the ordered assembly of around 80 ribosomal proteins. Supplying equimolar amounts of assembly-competent ribosomal proteins is complicated by their aggregation propensity and the spatial separation of their location of synthesis and pre-ribosome incorporation. Recent evidence has highlighted that dedicated chaperones protect individual, unassembled ribosomal proteins on their path to the pre-ribosomal assembly site. Here, we show that the co-translational recognition of Rpl3 and Rpl4 by their respective dedicated chaperone, Rrb1 or Acl4, reduces the degradation of the encoding RPL3 and RPL4 mRNAs in the yeast Saccharomyces cerevisiae. In both cases, negative regulation of mRNA levels occurs when the availability of the dedicated chaperone is limited and the nascent ribosomal protein is instead accessible to a regulatory machinery consisting of the nascent-polypeptide-associated complex and the Caf130-associated Ccr4-Not complex. Notably, deregulated expression of Rpl3 and Rpl4 leads to their massive aggregation and a perturbation of overall proteostasis in cells lacking the E3 ubiquitin ligase Tom1. Taken together, we have uncovered an unprecedented regulatory mechanism that adjusts the de novo synthesis of Rpl3 and Rpl4 to their actual consumption during ribosome assembly and, thereby, protects cells from the potentially detrimental effects of their surplus production.

Living cells are packed full of molecules known as proteins, which perform many vital tasks the cells need to survive and grow. Machines called ribosomes inside the cells use template molecules called messenger RNAs (or mRNAs for short) to produce proteins. The newly-made proteins then have to travel to a specific location in the cell to perform their tasks. Some newly-made proteins are prone to forming clumps, so cells have other proteins known as chaperones that ensure these clumps do not form.

The ribosomes themselves are made up of several proteins, some of which are also prone to clumping as they are being produced. To prevent this from happening, cells control how many ribosomal proteins they make, so there are just enough to form the ribosomes the cell needs at any given time. Previous studies found that, in yeast, two ribosomal proteins called Rpl3 and Rpl4 each have their own dedicated chaperone to prevent them from clumping. However, it remained unclear whether these chaperones are also involved in regulating the levels of Rpl3 and Rpl4.

To address this question, Pillet et al. studied both of these dedicated chaperones in yeast cells. The experiments showed that the chaperones bound to their target proteins (either units of Rpl3 or Rpl4) as they were being produced on the ribosomes. This protected the template mRNAs the ribosomes were using to produce these proteins from being destroyed, thus allowing further units of Rpl3 and Rpl4 to be produced. When enough Rpl3 and Rpl4 units were made, there were not enough of the chaperones to bind them all, leaving the mRNA templates unprotected. This led to the destruction of the mRNA templates, which decreased the numbers of Rpl3 and Rpl4 units being produced.

The work of Pillet et al. reveals a feedback mechanism that allows yeast to tightly control the levels of Rpl3 and Rpl4. In the future, these findings may help us understand diseases caused by defects in ribosomal proteins, such as Diamond-Blackfan anemia, and possibly also neurodegenerative diseases caused by clumps of proteins forming in cells. The next step will be to find out whether the mechanism uncovered by Pillet et al. also exists in human and other mammalian cells.

Signaling via the FLP-14/FRPR-19 neuropeptide pathway sustains nociceptive response to repeated noxious stimuli in C. elegans

In order to thrive in constantly changing environments, animals must adaptively respond to threatening events. Noxious stimuli are not only processed according to their absolute intensity, but also to their context. Adaptation processes can cause animals to habituate at different rates and degrees in response to permanent or repeated stimuli. Here, we used a forward genetic approach in Caenorhabditis elegans to identify a neuropeptidergic pathway, essential to prevent fast habituation and maintain robust withdrawal responses to repeated noxious stimuli. This pathway involves the FRPR-19A and FRPR-19B G-protein coupled receptor isoforms produced from the frpr-19 gene by alternative splicing. Loss or overexpression of each or both isoforms can impair withdrawal responses caused by the optogenetic activation of the polymodal FLP nociceptor neuron. Furthermore, we identified FLP-8 and FLP-14 as FRPR-19 ligands in vitro. flp-14, but not flp-8, was essential to promote withdrawal response and is part of the same genetic pathway as frpr-19 in vivo. Expression and cell-specific rescue analyses suggest that FRPR-19 acts both in the FLP nociceptive neurons and downstream interneurons, whereas FLP-14 acts from interneurons. Importantly, genetic impairment of the FLP-14/FRPR-19 pathway accelerated the habituation to repeated FLP-specific optogenetic activation, as well as to repeated noxious heat and harsh touch stimuli. Collectively, our data suggest that well-adjusted neuromodulation via the FLP-14/FRPR-19 pathway contributes to promote nociceptive signals in C. elegans and counteracts habituation processes that otherwise tend to rapidly reduce aversive responses to repeated noxious stimuli.

Minimalistic mycoplasmas harbor different functional toxin-antitoxin systems

Mycoplasmas are minute bacteria controlled by very small genomes ranging from 0.6 to 1.4 Mbp. They encompass several important medical and veterinary pathogens that are often associated with a wide range of chronic diseases. The long persistence of mycoplasma cells in their hosts can exacerbate the spread of antimicrobial resistance observed for many species. However, the nature of the virulence factors driving this phenomenon in mycoplasmas is still unclear. Toxin-antitoxin systems (TA systems) are genetic elements widespread in many bacteria that were historically associated with bacterial persistence. Their presence on mycoplasma genomes has never been carefully assessed, especially for pathogenic species. Here we investigated three candidate TA systems in M. mycoides subsp. capri encoding a (i) novel AAA-ATPase/subtilisin-like serine protease module, (ii) a putative AbiEii/AbiEi pair and (iii) a putative Fic/RelB pair. We sequence analyzed fourteen genomes of M. mycoides subsp. capri and confirmed the presence of at least one TA module in each of them. Interestingly, horizontal gene transfer signatures were also found in several genomic loci containing TA systems for several mycoplasma species. Transcriptomic and proteomic data confirmed differential expression profiles of these TA systems during mycoplasma growth in vitro. While the use of heterologous expression systems based on E. coli and B. subtilis showed clear limitations, the functionality and neutralization capacities of all three candidate TA systems were successfully confirmed using M. capricolum subsp. capricolum as a host. Additionally, M. capricolum subsp. capricolum was used to confirm the presence of functional TA system homologs in mycoplasmas of the Hominis and Pneumoniae phylogenetic groups. Finally, we showed that several of these M. mycoides subsp. capri toxins tested in this study, and particularly the subtilisin-like serine protease, could be used to establish a kill switch in mycoplasmas for industrial applications.

Revealing NOTCH-dependencies in synaptic targets associated with Alzheimer's disease

Recent studies have identified NOTCH signaling as a contributor of neurodegeneration including Alzheimer's disease' (AD) pathophysiology. As part of the efforts to understand molecular mechanisms and players involved in neurodegenerative dementia, we employed transgenic mouse models with Notch1 and Rbpjk loss of function (LOF) mutation in pyramidal neurons of the CA fields. Using RNA-seq, we have investigated the differential expression of NOTCH-dependent genes either upon environmental enrichment (EE) or upon kainic acid (KA) injury. We found a substantial genetic diversity in absence of both NOTCH1 receptor or RBPJK transcriptional activator. Among differentially expressed genes, we observed a significant upregulation of Gabra2a in both knockout models, suggesting a role for NOTCH signaling in the modulation of E/I balance. Upon excitotoxic stimulation, loss of RBPJK results in decreased expression of synaptic proteins with neuroprotective effects. We confirmed Nptx2, Npy, Pdch8, TncC as direct NOTCH1/RBPJK targets and Bdnf and Scg2 as indirect targets. Finally, we translate these findings into human entorhinal cortex containing the hippocampal region from AD patients performing targeted transcripts analysis. We observe an increased trend for RBPJK and the ligand DNER starting in the mild-moderate stage of the disease with no change of NOTCH1 expression. Alongside, expression of the Notch targets Hes5 and Hey1 tend to rise in the intermediate stage of the disease and drop in severe AD. Similarly the newly discovered NOTCH targets, NPTX2, NPY, BDNF show an up-warding tendency during the mild-moderate stage, and decline in the severe phase of the disease. This study identifies NOTCH as a central signaling cascade capable of modulating synaptic transmission in response to excitatory insult through the activation of neuroprotective genes that have been associated to AD.

Genome Assembly of the Ragweed Leaf Beetle: A Step Forward to Better Predict Rapid Evolution of a Weed Biocontrol Agent to Environmental Novelties

Rapid evolution of weed biological control agents (BCAs) to new biotic and abiotic conditions is poorly understood and so far only little considered both in pre-release and post-release studies, despite potential major negative or positive implications for risks of nontargeted attacks or for colonizing yet unsuitable habitats, respectively. Provision of genetic resources, such as assembled and annotated genomes, is essential to assess potential adaptive processes by identifying underlying genetic mechanisms. Here, we provide the first sequenced genome of a phytophagous insect used as a BCA, that is, the leaf beetle Ophraella communa, a promising BCA of common ragweed, recently and accidentally introduced into Europe. A total 33.98 Gb of raw DNA sequences, representing ∼43-fold coverage, were obtained using the PacBio SMRT-Cell sequencing approach. Among the five different assemblers tested, the SMARTdenovo assembly displaying the best scores was then corrected with Illumina short reads. A final genome of 774 Mb containing 7,003 scaffolds was obtained. The reliability of the final assembly was then assessed by benchmarking universal single-copy orthologous genes (>96.0% of the 1,658 expected insect genes) and by remapping tests of Illumina short reads (average of 98.6 ± 0.7% without filtering). The number of protein-coding genes of 75,642, representing 82% of the published antennal transcriptome, and the phylogenetic analyses based on 825 orthologous genes placing O. communa in the monophyletic group of Chrysomelidae, confirm the relevance of our genome assembly. Overall, the genome provides a valuable resource for studying potential risks and benefits of this BCA facing environmental novelties.

A rationally designed oral vaccine induces immunoglobulin A in the murine gut that directs the evolution of attenuated Salmonella variants

The ability of gut bacterial pathogens to escape immunity by antigenic variation-particularly via changes to surface-exposed antigens-is a major barrier to immune clearance¹. However, not all variants are equally fit in all environments^2,3. It should therefore be possible to exploit such immune escape mechanisms to direct an evolutionary trade-off. Here, we demonstrate this phenomenon using Salmonella enterica subspecies enterica serovar Typhimurium (S.Tm). A dominant surface antigen of S.Tm is its O-antigen: a long, repetitive glycan that can be rapidly varied by mutations in biosynthetic pathways or by phase variation^4,5. We quantified the selective advantage of O-antigen variants in the presence and absence of O-antigen-specific immunoglobulin A and identified a set of evolutionary trajectories allowing immune escape without an associated fitness cost in naive mice. Through the use of rationally designed oral vaccines, we induced immunoglobulin A responses blocking all of these trajectories. This selected for Salmonella mutants carrying deletions of the O-antigen polymerase gene wzyB. Due to their short O-antigen, these evolved mutants were more susceptible to environmental stressors (detergents or complement) and predation (bacteriophages) and were impaired in gut colonization and virulence in mice. Therefore, a rationally induced cocktail of intestinal antibodies can direct an evolutionary trade-off in S.Tm. This lays the foundations for the exploration of mucosal vaccines capable of setting evolutionary traps as a prophylactic strategy.

Statistical and Machine Learning Techniques in Human Microbiome Studies: Contemporary Challenges and Solutions

The human microbiome has emerged as a central research topic in human biology and biomedicine. Current microbiome studies generate high-throughput omics data across different body sites, populations, and life stages. Many of the challenges in microbiome research are similar to other high-throughput studies, the quantitative analyses need to address the heterogeneity of data, specific statistical properties, and the remarkable variation in microbiome composition across individuals and body sites. This has led to a broad spectrum of statistical and machine learning challenges that range from study design, data processing, and standardization to analysis, modeling, cross-study comparison, prediction, data science ecosystems, and reproducible reporting. Nevertheless, although many statistics and machine learning approaches and tools have been developed, new techniques are needed to deal with emerging applications and the vast heterogeneity of microbiome data. We review and discuss emerging applications of statistical and machine learning techniques in human microbiome studies and introduce the COST Action CA18131 "ML4Microbiome" that brings together microbiome researchers and machine learning experts to address current challenges such as standardization of analysis pipelines for reproducibility of data analysis results, benchmarking, improvement, or development of existing and new tools and ontologies.

NGS-Based S. aureus Typing and Outbreak Analysis in Clinical Microbiology Laboratories: Lessons Learned From a Swiss-Wide Proficiency Test

Whole genome sequencing (WGS) enables high resolution typing of bacteria up to the single nucleotide polymorphism (SNP) level. WGS is used in clinical microbiology laboratories for infection control, molecular surveillance and outbreak analyses. Given the large palette of WGS reagents and bioinformatics tools, the Swiss clinical bacteriology community decided to conduct a ring trial (RT) to foster harmonization of NGS-based bacterial typing. The RT aimed at assessing methicillin-susceptible Staphylococcus aureus strain relatedness from WGS and epidemiological data. The RT was designed to disentangle the variability arising from differences in sample preparation, SNP calling and phylogenetic methods. Nine laboratories participated. The resulting phylogenetic tree and cluster identification were highly reproducible across the laboratories. Cluster interpretation was, however, more laboratory dependent, suggesting that an increased sharing of expertise across laboratories would contribute to further harmonization of practices. More detailed bioinformatic analyses unveiled that while similar clusters were found across laboratories, these were actually based on different sets of SNPs, differentially retained after sample preparation and SNP calling procedures. Despite this, the observed number of SNP differences between pairs of strains, an important criterion to determine strain relatedness given epidemiological information, was similar across pipelines for closely related strains when restricting SNP calls to a common core genome defined by S. aureus cgMLST schema. The lessons learned from this pilot study will serve the implementation of larger-scale RT, as a mean to have regular external quality assessments for laboratories performing WGS analyses in a clinical setting.

Genome-based characterization of two Colombian clinical Providencia rettgeri isolates co-harboring NDM-1, VIM-2, and other β-lactamases

BACKGROUND

Providencia rettgeri is a nosocomial pathogen associated with urinary tract infections and related to Healthcare-Associated Infection (HAI). In recent years isolates producing New Delhi Metallo-β-lactamase (NDM) and other β-lactamases have been reported that reduce the efficiency of clinical antimicrobial treatments. In this study, we analyzed antibiotic resistance, the presence of resistance genes and the clonal relationship of two P. rettgeri isolates obtained from male patients admitted to the same hospital in Bogotá - Colombia, 2015.

RESULTS

Antibiotic susceptibility profile evaluated by the Kirby-Bauer method revealed that both isolates were resistant to third-generation carbapenems and cephalosporins. Whole-genome sequencing (Illumina HiSeq) followed by SPAdes assembling, Prokka annotation in combination with an in-house Python program and resistance gene detection by ResFinder identified the same six β-lactamase genes in both isolates: blaNDM-1, blaVIM-2, blaCTX-M-15, blaOXA-10, blaCMY-2 and blaTEM-1. Additionally, various resistance genes associated with antibiotic target alteration (arnA, PmrE, PmrF, LpxA, LpxC, gyrB, folP, murA, rpoB, rpsL, tet34) were found and four efflux pumps (RosAB, EmrD, mdtH and cmlA). The additional resistance to gentamicin in one of the two isolates could be explained by a detected SNP in CpxA (Cys191Arg) which is involved in the stress response of the bacterial envelope. Genome BLAST comparison using CGView, the ANI value (99.99%) and the pangenome (using Roary) phylogenetic tree (same clade, small distance) showed high similarity between the isolates. The rMLST analysis indicated that both isolates were typed as rST-61,696, same as the RB151 isolate previously isolated in Bucaramanga, Colombia, 2013, and the FDAARGOS_330 isolate isolated in the USA, 2015.

CONCLUSIONS

We report the coexistence of the carbapenemase genes blaNDM-1, and blaVIM-2, together with the β-lactamase genes blaCTX-M-15, blaOXA-10, blaCMY-2 and blaTEM-1, in P. rettgeri isolates from two patients in Colombia. Whole-genome sequence analysis indicated a circulation of P. rettgeri rST-61,696 strains in America that needs to be investigated further.

Injured Axons Instruct Schwann Cells to Build Constricting Actin Spheres to Accelerate Axonal Disintegration

After a peripheral nerve lesion, distal ends of injured axons disintegrate into small fragments that are subsequently cleared by Schwann cells and later by macrophages. Axonal debris clearing is an early step of the repair process that facilitates regeneration. We show here that Schwann cells promote distal cut axon disintegration for timely clearing. By combining cell-based and in vivo models of nerve lesion with mouse genetics, we show that this mechanism is induced by distal cut axons, which signal to Schwann cells through PlGF mediating the activation and upregulation of VEGFR1 in Schwann cells. In turn, VEGFR1 activates Pak1, leading to the formation of constricting actomyosin spheres along unfragmented distal cut axons to mediate their disintegration. Interestingly, oligodendrocytes can acquire a similar behavior as Schwann cells by enforced expression of VEGFR1. These results thus identify controllable molecular cues of a neuron-glia crosstalk essential for timely clearing of damaged axons.

Study protocol for the ABERRANT study: antibiotic-induced disruption of the maternal and infant microbiome and adverse health outcomes - a prospective cohort study among children born at term

INTRODUCTION

There is compositional overlap between the maternal intestinal microbiome, the breast milk microbiome and the infant oral and intestinal microbiome. Antibiotics cause profound changes in the microbiome. However, the effect of intrapartum and early-life antibiotics on the maternal intestinal and breast milk microbiome, and the infant oral and intestinal microbiome, and whether effects are only short term or persist long term remain uncertain.

METHODS AND ANALYSES

In this prospective cohort study, we will use metagenomic sequencing to determine: (1) the effect of intrapartum antibiotics on the composition of the breast milk, and the infant oral and intestinal microbiome, including the development and persistence of antibiotic resistance; (2) the effect of antibiotic exposure in the first year of life on the composition of the infant oral and intestinal microbiome, including the development and persistence of antibiotic resistance; (3) the effect of disruption of the infant oral and intestinal microbiome on health outcomes and (4) the compositional overlap between the maternal intestinal microbiome, the breast milk microbiome and the infant oral and intestinal microbiome.

ETHICS AND DISSEMINATION

The ABERRANT study has been approved by the commission cantonale d'éthique de la recherche sur l'être humain (CER-VD) du Canton de Vaud (#2019-01567). Outcomes will be disseminated through publication and will be presented at scientific conferences.

TRIAL REGISTRATION NUMBER

NCT04091282.

Large-Scale Comparison of Toxin and Antitoxins in Listeria monocytogenes

Toxin–antitoxin systems (TASs) are widely distributed in prokaryotes and encode pairs of genes involved in many bacterial biological processes and mechanisms, including pathogenesis. The TASs have not been extensively studied in Listeria monocytogenes (Lm), a pathogenic bacterium of the Firmicutes phylum causing infections in animals and humans. Using our recently published TASmania database, we focused on the known and new putative TASs in 352 Listeria monocytogenes genomes and identified the putative core gene TASs (cgTASs) with the Pasteur BIGSdb-Lm database and, by complementarity, the putative accessory gene TAS (acTASs). We combined the cgTASs with those of an additional 227 L. monocytogenes isolates from our previous studies containing metadata information. We discovered that the differences in 14 cgTAS alleles are sufficient to separate the four main lineages of Listeria monocytogenes. Analyzing these differences in more details, we uncovered potentially co-evolving residues in some pairs of proteins in cgTASs, probably essential for protein–protein interactions within the TAS complex.

Identification of avoidance genes through neural pathway-specific forward optogenetics

Understanding how the nervous system bridges sensation and behavior requires the elucidation of complex neural and molecular networks. Forward genetic approaches, such as screens conducted in C. elegans, have successfully identified genes required to process natural sensory stimuli. However, functional redundancy within the underlying neural circuits, which are often organized with multiple parallel neural pathways, limits our ability to identify 'neural pathway-specific genes', i.e. genes that are essential for the function of some, but not all of these redundant neural pathways. To overcome this limitation, we developed a 'forward optogenetics' screening strategy in which natural stimuli are initially replaced by the selective optogenetic activation of a specific neural pathway. We used this strategy to address the function of the polymodal FLP nociceptors mediating avoidance of noxious thermal and mechanical stimuli. According to our expectations, we identified both mutations in 'general' avoidance genes that broadly impact avoidance responses to a variety of natural noxious stimuli (unc-4, unc-83, and eat-4) and mutations that produce a narrower impact, more restricted to the FLP pathway (syd-2, unc-14 and unc-68). Through a detailed follow-up analysis, we further showed that the Ryanodine receptor UNC-68 acts cell-autonomously in FLP to adjust heat-evoked calcium signals and aversive behaviors. As a whole, our work (i) reveals the importance of properly regulated ER calcium release for FLP function, (ii) provides new entry points for new nociception research and (iii) demonstrates the utility of our forward optogenetic strategy, which can easily be transposed to analyze other neural pathways.

deltaRpkm: an R package for a rapid detection of differential gene presence between related bacterial genomes

Background

Comparative genomics has seen the development of many software performing the clustering, polymorphism and gene content analysis of genomes at different phylogenetic levels (isolates, species). These tools rely on de novo assembly and/or multiple alignments that can be computationally intensive for large datasets. With a large number of similar genomes in particular, e.g., in surveillance and outbreak detection, assembling each genome can become a redundant and expensive step in the identification of genes potentially involved in a given clinical feature.

Results

We have developed deltaRpkm, an R package that performs a rapid differential gene presence evaluation between two large groups of closely related genomes. Starting from a standard gene count table, deltaRpkm computes the RPKM per gene per sample, then the inter-group δRPKM values, the corresponding median δRPKM (m) for each gene and the global standard deviation value of m (s_m). Genes with m >  = 2 ∗ s_m (standard deviation s of all the m values) are considered as “differentially present” in the reference genome group. Our simple yet effective method of differential RPKM has been successfully applied in a recent study published by our group (N = 225 genomes of Listeria monocytogenes) (Aguilar-Bultet et al. Front Cell Infect Microbiol 8:20, 2018).

Conclusions

To our knowledge, deltaRpkm is the first tool to propose a straightforward inter-group differential gene presence analysis with large datasets of related genomes, including non-coding genes, and to output directly a list of genes potentially involved in a phenotype.

Large-Scale Analysis of the Mycoplasma bovis Genome Identified Non-essential, Adhesion- and Virulence-Related Genes

Mycoplasma bovis is an important pathogen of cattle causing bovine mycoplasmosis. Clinical manifestations are numerous, but pneumonia, mastitis, and arthritis cases are mainly reported. Currently, no efficient vaccine is available and antibiotic treatments are not always satisfactory. The design of new, efficient prophylactic and therapeutic approaches requires a better understanding of the molecular mechanisms responsible for M. bovis pathogenicity. Random transposon mutagenesis has been widely used in Mycoplasma species to identify potential gene functions. Such an approach can also be used to screen genomes and search for essential and non-essential genes for growth. Here, we generated a random transposon mutant library of M. bovis strain JF4278 containing approximately 4000 independent insertion sites. We then coupled high-throughput screening of this mutant library to transposon sequencing and bioinformatic analysis to identify M. bovis non-essential, adhesion- and virulence-related genes. Three hundred and fifty-two genes of M. bovis were assigned as essential for growth in rich medium. Among the remaining non-essential genes, putative virulence-related factors were subsequently identified. The complete mutant library was screened for adhesion using primary bovine mammary gland epithelial cells. Data from this assay resulted in a list of conditional-essential genes with putative adhesion-related functions by identifying non-essential genes for growth that are essential for host cell-adhesion. By individually assessing the adhesion capacity of six selected mutants, two previously unknown factors and the adhesin TrmFO were associated with a reduced adhesion phenotype. Overall, our study (i) uncovers new, putative virulence-related genes; (ii) offers a list of putative adhesion-related factors; and (iii) provides valuable information for vaccine design and for exploring M. bovis biology, pathogenesis, and host-interaction.

TASmania: A bacterial Toxin-Antitoxin Systems database

Bacterial Toxin-Antitoxin systems (TAS) are involved in key biological functions including plasmid maintenance, defense against phages, persistence and virulence. They are found in nearly all phyla and classified into 6 different types based on the mode of inactivation of the toxin, with the type II TAS being the best characterized so far. We have herein developed a new in silico discovery pipeline named TASmania, which mines the >41K assemblies of the EnsemblBacteria database for known and uncharacterized protein components of type I to IV TAS loci. Our pipeline annotates the proteins based on a list of curated HMMs, which leads to >2^.10⁶ loci candidates, including orphan toxins and antitoxins, and organises the candidates in pseudo-operon structures in order to identify new TAS candidates based on a guilt-by-association strategy. In addition, we classify the two-component TAS with an unsupervised method on top of the pseudo-operon (pop) gene structures, leading to 1567 “popTA” models offering a more robust classification of the TAs families. These results give valuable clues in understanding the toxin/antitoxin modular structures and the TAS phylum specificities. Preliminary in vivo work confirmed six putative new hits in Mycobacterium tuberculosis as promising candidates. The TASmania database is available on the following server https://shiny.bioinformatics.unibe.ch/apps/tasmania/.

TASmania offers an extensive annotation of TA loci in a very large database of bacterial genomes, which represents a resource of crucial importance for the microbiology community. TASmania supports i) the discovery of new TA families; ii) the design of a robust experimental strategy by taking into account potential interferences in trans; iii) the comparative analysis between TA loci content, phylogeny and/or phenotypes (pathogenicity, persistence, stress resistance, associated host types) by providing a vast repertoire of annotated assemblies. Our database contains TA annotations of a given strain not only mapped to its core genome but also to its plasmids, whenever applicable.

Removal of a Subset of Non-essential Genes Fully Attenuates a Highly Virulent Mycoplasma Strain

Mycoplasmas are the smallest free-living organisms and cause a number of economically important diseases affecting humans, animals, insects, and plants. Here, we demonstrate that highly virulent Mycoplasma mycoides subspecies capri (Mmc) can be fully attenuated via targeted deletion of non-essential genes encoding, among others, potential virulence traits. Five genomic regions, representing approximately 10% of the original Mmc genome, were successively deleted using Saccharomyces cerevisiae as an engineering platform. Specifically, a total of 68 genes out of the 432 genes verified to be individually non-essential in the JCVI-Syn3.0 minimal cell, were excised from the genome. In vitro characterization showed that this mutant was similar to its parental strain in terms of its doubling time, even though 10% of the genome content were removed. A novel in vivo challenge model in goats revealed that the wild-type parental strain caused marked necrotizing inflammation at the site of inoculation, septicemia and all animals reached endpoint criteria within 6 days after experimental infection. This is in contrast to the mutant strain, which caused no clinical signs nor pathomorphological lesions. These results highlight, for the first time, the rational design, construction and complete attenuation of a Mycoplasma strain via synthetic genomics tools. Trait addition using the yeast-based genome engineering platform and subsequent in vitro or in vivo trials employing the Mycoplasma chassis will allow us to dissect the role of individual candidate Mycoplasma virulence factors and lead the way for the development of an attenuated designer vaccine.

Reproduction of contagious caprine pleuropneumonia reveals the ability of convalescent sera to reduce hydrogen peroxide production in vitro

Contagious caprine pleuropneumonia (CCPP), caused by Mycoplasma capricolum subsp. capripneumoniae is a severe disease widespread in Africa and Asia. Limited knowledge is available on the pathogenesis of this organism, mainly due to the lack of a robust in vivo challenge model and the means to do site-directed mutagenesis. This work describes the establishment of a novel caprine challenge model for CCPP that resulted in 100% morbidity using a combination of repeated intranasal spray infection followed by a single transtracheal infection employing the recent Kenyan outbreak strain ILRI181. Diseased animals displayed CCPP-related pathology and the bacteria could subsequently be isolated from pleural exudates and lung tissues in concentrations of up to 10⁹ bacteria per mL as well as in the trachea using immunohistochemistry. Reannotation of the genome sequence of ILRI181 and F38^T revealed the existence of genes encoding the complete glycerol uptake and metabolic pathways involved in hydrogen peroxide (H₂O₂) production in the phylogenetically related pathogen M. mycoides subsp. mycoides. Furthermore, the expression of L-α-glycerophosphate oxidase (GlpO) in vivo was confirmed. In addition, the function of the glycerol metabolism was verified by measurement of production of H₂O₂ in medium containing physiological serum concentrations of glycerol. Peroxide production could be inhibited with serum from convalescent animals. These results will pave the way for a better understanding of host-pathogen interactions during CCPP and subsequent vaccine development.

Suppressor mutations in Rpf2-Rrs1 or Rpl5 bypass the Cgr1 function for pre-ribosomal 5S RNP-rotation

During eukaryotic 60S biogenesis, the 5S RNP requires a large rotational movement to achieve its mature position. Cryo-EM of the Rix1-Rea1 pre-60S particle has revealed the post-rotation stage, in which a gently undulating α-helix corresponding to Cgr1 becomes wedged between Rsa4 and the relocated 5S RNP, but the purpose of this insertion was unknown. Here, we show that cgr1 deletion in yeast causes a slow-growth phenotype and reversion of the pre-60S particle to the pre-rotation stage. However, spontaneous extragenic suppressors could be isolated, which restore growth and pre-60S biogenesis in the absence of Cgr1. Whole-genome sequencing reveals that the suppressor mutations map in the Rpf2–Rrs1 module and Rpl5, which together stabilize the unrotated stage of the 5S RNP. Thus, mutations in factors stabilizing the pre-rotation stage facilitate 5S RNP relocation upon deletion of Cgr1, but Cgr1 itself could stabilize the post-rotation stage.

During biogenesis of the eukaryotic 60S ribosome, a large rotational movement of the 5S RNP is required to achieve its mature position. By analyzing extragenic suppressors of crg1—a key factor required for rotation—the authors provide mechanistic insight into a key step of ribosome biogenesis.

Genetic Separation of Listeria monocytogenes Causing Central Nervous System Infections in Animals

Listeria monocytogenes is a foodborne pathogen that causes abortion, septicemia, gastroenteritis and central nervous system (CNS) infections in ruminants and humans. L. monocytogenes strains mainly belong to two distinct phylogenetic groups, named lineages I and II. In general, clinical cases in humans and animals, in particular CNS infections, are caused by lineage I strains, while most of the environmental and food strains belong to lineage II. Little is known about why lineage I is more virulent than lineage II, even though various molecular factors and mechanisms associated with pathogenesis are known. In this study, we have used a variety of whole genome sequence analyses and comparative genomic tools in order to find characteristics that distinguish lineage I from lineage II strains and CNS infection strains from non-CNS strains. We analyzed 225 strains and identified single nucleotide variants between lineages I and II, as well as differences in the gene content. Using a novel approach based on Reads Per Kilobase per Million Mapped (RPKM), we identified 167 genes predominantly absent in lineage II but present in lineage I. These genes are mostly encoding for membrane-associated proteins. Additionally, we found 77 genes that are largely absent in the non-CNS associated strains, while 39 genes are especially lacking in our defined “non-clinical” group. Based on the RPKM analysis and the metadata linked to the L. monocytogenes strains, we identified 6 genes potentially associated with CNS cases, which include a transcriptional regulator, an ABC transporter and a non-coding RNA. Although there is not a clear separation between pathogenic and non-pathogenic strains based on phylogenetic lineages, the presence of the genes identified in our study reveals potential pathogenesis traits in ruminant L. monocytogenes strains. Ultimately, the differences that we have found in our study will help steer future studies in understanding the virulence mechanisms of the most pathogenic L. monocytogenes strains.

Malaria infected red blood cells release small regulatory RNAs through extracellular vesicles

The parasite Plasmodium falciparum causes the most severe form of malaria. Cell communication between parasites is an important mechanism to control population density and differentiation. The infected red blood cells (iRBCs) release small extracellular vesicles (EVs) that transfer cargoes between cells. The EVs synchronize the differentiation of the asexual parasites into gametocytes to initiate the transmission to the mosquito. Beside their role in parasite communication, EVs regulate vascular function. So far, the exact cargoes responsible for cellular communication remain unknown. We isolated EVs from cultured iRBCs to determine their small RNA content. We identified several types of human and plasmodial regulatory RNAs. While the miRNAs and tRNA-derived fragments were the most abundant human RNAs, we also found Y-RNAs, vault RNAs, snoRNAs and piRNAs. Interestingly, we found about 120 plasmodial RNAs, including mRNAs coding for exported proteins and proteins involved in drug resistance, as well as non-coding RNAs, such as rRNAs, small nuclear (snRNAs) and tRNAs. These data show, that iRBC-EVs carry small regulatory RNAs. A role in cellular communication is possible since the RNAs were transferred to endothelial cells. Furthermore, the presence of Plasmodium RNAs, in EVs suggests that they may be used as biomarker to track and detect disease.

Clostridium chauvoei, an Evolutionary Dead-End Pathogen

Full genome sequences of 20 strains of Clostridium chauvoei, the etiological agent of blackleg of cattle and sheep, isolated from four different continents over a period of 64 years (1951–2015) were determined and analyzed. The study reveals that the genome of the species C. chauvoei is highly homogeneous compared to the closely related species C. perfringens, a widespread pathogen that affects human and many animal species. Analysis of the CRISPR locus is sufficient to differentiate most C. chauvoei strains and is the most heterogenous region in the genome, containing in total 187 different spacer elements that are distributed as 30 – 77 copies in the various strains. Some genetic differences are found in the 3 allelic variants of fliC1, fliC2 and fliC3 genes that encode structural flagellin proteins, and certain strains do only contain one or two alleles. However, the major virulence genes including the highly toxic C.chauvoei toxin A, the sialidase and the two hyaluronidases are fully conserved as are the metabolic and structural genes of C. chauvoei. These data indicate that C. chauvoei is a strict ruminant-associated pathogen that has reached a dead end in its evolution.

Whole-Genome Sequence of a Colombian Acinetobacter baumannii Strain, a Coproducer of OXA-72 and OXA-255-Like Carbapenemases

Colombian Acinetobacter baumannii strain ST920 was isolated from the sputum of a 68-year-old male patient. This isolate possessed blaOXA-72 and blaOXA-255-like genes. The assembled genome contained 4,104,098 pb and 38.79% G+C content. This is the first case reported of the coproduction (blaOXA-72 and blaOXA-255-like) of carbapenem-hydrolyzing class D β-lactamases (CHDLs) in Acinetobacter baumannii.

The microbiome of the leaf surface of Arabidopsis protects against a fungal pathogen

We have explored the importance of the phyllosphere microbiome in plant resistance in the cuticle mutants bdg (BODYGUARD) or lacs2.3 (LONG CHAIN FATTY ACID SYNTHASE 2) that are strongly resistant to the fungal pathogen Botrytis cinerea. The study includes infection of plants under sterile conditions, 16S ribosomal DNA sequencing of the phyllosphere microbiome, and isolation and high coverage sequencing of bacteria from the phyllosphere. When inoculated under sterile conditions bdg became as susceptible as wild-type (WT) plants whereas lacs2.3 mutants retained the resistance. Adding washes of its phyllosphere microbiome could restore the resistance of bdg mutants, whereas the resistance of lacs2.3 results from endogenous mechanisms. The phyllosphere microbiome showed distinct populations in WT plants compared to cuticle mutants. One species identified as Pseudomonas sp isolated from the microbiome of bdg provided resistance to B. cinerea on Arabidopsis thaliana as well as on apple fruits. No direct activity was observed against B. cinerea and the action of the bacterium required the plant. Thus, microbes present on the plant surface contribute to the resistance to B. cinerea. These results open new perspectives on the function of the leaf microbiome in the protection of plants.

The Dedicated Chaperone Acl4 Escorts Ribosomal Protein Rpl4 to Its Nuclear Pre-60S Assembly Site

Ribosomes are the highly complex macromolecular assemblies dedicated to the synthesis of all cellular proteins from mRNA templates. The main principles underlying the making of ribosomes are conserved across eukaryotic organisms and this process has been studied in most detail in the yeast Saccharomyces cerevisiae. Yeast ribosomes are composed of four ribosomal RNAs (rRNAs) and 79 ribosomal proteins (r-proteins). Most r-proteins need to be transported from the cytoplasm to the nucleus where they get incorporated into the evolving pre-ribosomal particles. Due to the high abundance and difficult physicochemical properties of r-proteins, their correct folding and fail-safe targeting to the assembly site depends largely on general, as well as highly specialized, chaperone and transport systems. Many r-proteins contain universally conserved or eukaryote-specific internal loops and/or terminal extensions, which were shown to mediate their nuclear targeting and association with dedicated chaperones in a growing number of cases. The 60S r-protein Rpl4 is particularly interesting since it harbours a conserved long internal loop and a prominent C-terminal eukaryote-specific extension. Here we show that both the long internal loop and the C-terminal eukaryote-specific extension are strictly required for the functionality of Rpl4. While Rpl4 contains at least five distinct nuclear localization signals (NLS), the C-terminal part of the long internal loop associates with a specific binding partner, termed Acl4. Absence of Acl4 confers a severe slow-growth phenotype and a deficiency in the production of 60S subunits. Genetic and biochemical evidence indicates that Acl4 can be considered as a dedicated chaperone of Rpl4. Notably, Acl4 localizes to both the cytoplasm and nucleus and it has the capacity to capture nascent Rpl4 in a co-translational manner. Taken together, our findings indicate that the dedicated chaperone Acl4 accompanies Rpl4 from the cytoplasm to its pre-60S assembly site in the nucleus.

The Petunia GRAS Transcription Factor ATA/RAM1 Regulates Symbiotic Gene Expression and Fungal Morphogenesis in Arbuscular Mycorrhiza

Arbuscular mycorrhiza (AM) is a mutual symbiosis that involves a complex symbiotic interface over which nutrients are exchanged between the plant host and the AM fungus. Dozens of genes in the host are required for the establishment and functioning of the interaction, among them nutrient transporters that mediate the uptake of mineral nutrients delivered by the fungal arbuscules. We have isolated in a genetic mutant screen a petunia (Petunia hybrida) Gibberellic Acid Insensitive, Repressor of Gibberellic Acid Insensitive, and Scarecrow (GRAS)-type transcription factor, Atypical Arbuscule (ATA), that acts as the central regulator of AM-related genes and is required for the morphogenesis of arbuscules. Forced mycorrhizal inoculations from neighboring wild-type plants revealed an additional role of ATA in restricting mycorrhizal colonization of the root meristem. The lack of ATA, which represents the ortholog of Required For Arbuscular Mycorrhiza1 in Medicago truncatula, renders the interaction completely ineffective, hence demonstrating the central role of AM-related genes for arbuscule development and function.

SV-AUTOPILOT: optimized, automated construction of structural variation discovery and benchmarking pipelines

BACKGROUND

Many tools exist to predict structural variants (SVs), utilizing a variety of algorithms. However, they have largely been developed and tested on human germline or somatic (e.g. cancer) variation. It seems appropriate to exploit this wealth of technology available for humans also for other species. Objectives of this work included: a) Creating an automated, standardized pipeline for SV prediction. b) Identifying the best tool(s) for SV prediction through benchmarking. c) Providing a statistically sound method for merging SV calls.

RESULTS

The SV-AUTOPILOT meta-tool platform is an automated pipeline for standardization of SV prediction and SV tool development in paired-end next-generation sequencing (NGS) analysis. SV-AUTOPILOT comes in the form of a virtual machine, which includes all datasets, tools and algorithms presented here. The virtual machine easily allows one to add, replace and update genomes, SV callers and post-processing routines and therefore provides an easy, out-of-the-box environment for complex SV discovery tasks. SV-AUTOPILOT was used to make a direct comparison between 7 popular SV tools on the Arabidopsis thaliana genome using the Landsberg (Ler) ecotype as a standardized dataset. Recall and precision measurements suggest that Pindel and Clever were the most adaptable to this dataset across all size ranges while Delly performed well for SVs larger than 250 nucleotides. A novel, statistically-sound merging process, which can control the false discovery rate, reduced the false positive rate on the Arabidopsis benchmark dataset used here by >60%.

CONCLUSION

SV-AUTOPILOT provides a meta-tool platform for future SV tool development and the benchmarking of tools on other genomes using a standardized pipeline. It optimizes detection of SVs in non-human genomes using statistically robust merging. The benchmarking in this study has demonstrated the power of 7 different SV tools for analyzing different size classes and types of structural variants. The optional merge feature enriches the call set and reduces false positives providing added benefit to researchers planning to validate SVs. SV-AUTOPILOT is a powerful, new meta-tool for biologists as well as SV tool developers.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

A novel bioinformatics pipeline to discover genes related to arbuscular mycorrhizal symbiosis based on their evolutionary conservation pattern among higher plants

BACKGROUND

Genes involved in arbuscular mycorrhizal (AM) symbiosis have been identified primarily by mutant screens, followed by identification of the mutated genes (forward genetics). In addition, a number of AM-related genes has been identified by their AM-related expression patterns, and their function has subsequently been elucidated by knock-down or knock-out approaches (reverse genetics). However, genes that are members of functionally redundant gene families, or genes that have a vital function and therefore result in lethal mutant phenotypes, are difficult to identify. If such genes are constitutively expressed and therefore escape differential expression analyses, they remain elusive. The goal of this study was to systematically search for AM-related genes with a bioinformatics strategy that is insensitive to these problems. The central element of our approach is based on the fact that many AM-related genes are conserved only among AM-competent species.

RESULTS

Our approach involves genome-wide comparisons at the proteome level of AM-competent host species with non-mycorrhizal species. Using a clustering method we first established orthologous/paralogous relationships and subsequently identified protein clusters that contain members only of the AM-competent species. Proteins of these clusters were then analyzed in an extended set of 16 plant species and ranked based on their relatedness among AM-competent monocot and dicot species, relative to non-mycorrhizal species. In addition, we combined the information on the protein-coding sequence with gene expression data and with promoter analysis. As a result we present a list of yet uncharacterized proteins that show a strongly AM-related pattern of sequence conservation, indicating that the respective genes may have been under selection for a function in AM. Among the top candidates are three genes that encode a small family of similar receptor-like kinases that are related to the S-locus receptor kinases involved in sporophytic self-incompatibility.

CONCLUSIONS

We present a new systematic strategy of gene discovery based on conservation of the protein-coding sequence that complements classical forward and reverse genetics. This strategy can be applied to diverse other biological phenomena if species with established genome sequences fall into distinguished groups that differ in a defined functional trait of interest.

Patho-genetics of Clostridium chauvoei

The genomic sequence of Clostridium chauvoei, the etiological agent of blackleg, a severe disease of ruminants with high mortality specified by a myonecrosis reveals a chromosome of 2.8 million base-pairs and a cryptic plasmid of 5.5 kilo base-pairs. The chromosome contains the main pathways like glycolysis/gluconeogenesis, sugar metabolism, purine and pyrimidine metabolisms, but the notable absence of genes of the citric acid cycle and deficient or partially deficient amino acid metabolism for Histidine, Tyrosine, Phenylalanine, and Tryptophan. These essential amino acids might be acquired from host tissue damage caused by various toxins and by protein metabolism that includes 57 genes for peptidases, and several ABC transporters for amino acids import.

Complete Genome Sequences of Virulent Mycoplasma capricolum subsp. capripneumoniae Strains F38 and ILRI181

Contagious caprine pleuropneumonia (CCPP) caused by Mycoplasma capricolum subsp. capripneumoniae is a severe epidemic affecting mainly domestic Caprinae species but also affects wild Caprinae species. M. capricolum subsp. capripneumoniae belongs to the "Mycoplasma mycoides cluster." The disease features prominently in East Africa, in particular Kenya, Tanzania, and Ethiopia. CCPP also endangers wildlife and thus affects not only basic nutritional resources of large populations but also expensively built-up game resorts in affected countries. Here, we report the complete sequences of two M. capricolum subsp. capripneumoniae strains: the type strain F38 and strain ILRI181 isolated druing a recent outbreak in Kenya. Both genomes have a G+C content of 24% with sizes of 1,016,760 bp and 1,017,183 bp for strains F38 and ILRI181, respectively.

Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef)

BACKGROUND

Tef (Eragrostis tef), an indigenous cereal critical to food security in the Horn of Africa, is rich in minerals and protein, resistant to many biotic and abiotic stresses and safe for diabetics as well as sufferers of immune reactions to wheat gluten. We present the genome of tef, the first species in the grass subfamily Chloridoideae and the first allotetraploid assembled de novo. We sequenced the tef genome for marker-assisted breeding, to shed light on the molecular mechanisms conferring tef's desirable nutritional and agronomic properties, and to make its genome publicly available as a community resource.

RESULTS

The draft genome contains 672 Mbp representing 87% of the genome size estimated from flow cytometry. We also sequenced two transcriptomes, one from a normalized RNA library and another from unnormalized RNASeq data. The normalized RNA library revealed around 38000 transcripts that were then annotated by the SwissProt group. The CoGe comparative genomics platform was used to compare the tef genome to other genomes, notably sorghum. Scaffolds comprising approximately half of the genome size were ordered by syntenic alignment to sorghum producing tef pseudo-chromosomes, which were sorted into A and B genomes as well as compared to the genetic map of tef. The draft genome was used to identify novel SSR markers, investigate target genes for abiotic stress resistance studies, and understand the evolution of the prolamin family of proteins that are responsible for the immune response to gluten.

CONCLUSIONS

It is highly plausible that breeding targets previously identified in other cereal crops will also be valuable breeding targets in tef. The draft genome and transcriptome will be of great use for identifying these targets for genetic improvement of this orphan crop that is vital for feeding 50 million people in the Horn of Africa.

Comparative genome analysis of Pseudomonas knackmussii B13, the first bacterium known to degrade chloroaromatic compounds

Pseudomonas knackmussii B13 was the first strain to be isolated in 1974 that could degrade chlorinated aromatic hydrocarbons. This discovery was the prologue for subsequent characterization of numerous bacterial metabolic pathways, for genetic and biochemical studies, and which spurred ideas for pollutant bioremediation. In this study, we determined the complete genome sequence of B13 using next generation sequencing technologies and optical mapping. Genome annotation indicated that B13 has a variety of metabolic pathways for degrading monoaromatic hydrocarbons including chlorobenzoate, aminophenol, anthranilate and hydroxyquinol, but not polyaromatic compounds. Comparative genome analysis revealed that B13 is closest to Pseudomonas denitrificans and Pseudomonas aeruginosa. The B13 genome contains at least eight genomic islands [prophages and integrative conjugative elements (ICEs)], which were absent in closely related pseudomonads. We confirm that two ICEs are identical copies of the 103 kb self-transmissible element ICEclc that carries the genes for chlorocatechol metabolism. Comparison of ICEclc showed that it is composed of a variable and a 'core' region, which is very conserved among proteobacterial genomes, suggesting a widely distributed family of so far uncharacterized ICE. Resequencing of two spontaneous B13 mutants revealed a number of single nucleotide substitutions, as well as excision of a large 220 kb region and a prophage that drastically change the host metabolic capacity and survivability.

Short term evolution of a highly transmissible methicillin-resistant Staphylococcus aureus clone (ST228) in a tertiary care hospital

Staphylococcus aureus is recognized as one of the major human pathogens and is by far one of the most common nosocomial organisms. The genetic basis for the emergence of highly epidemic strains remains mysterious. Studying the microevolution of the different clones of S. aureus is essential for identifying the forces driving pathogen emergence and spread. The aim of the present study was to determine the genetic changes characterizing a lineage belonging to the South German clone (ST228) that spread over ten years in a tertiary care hospital in Switzerland. For this reason, we compared the whole genome of eight isolates recovered between 2001 and 2008 at the Lausanne hospital. The genetic comparison of these isolates revealed that their genomes are extremely closely related. Yet, a few more important genetic changes, such as the replacement of a plasmid, the loss of large fragments of DNA, or the insertion of transposases, were observed. These transfers of mobile genetic elements shaped the evolution of the ST228 lineage that spread within the Lausanne hospital. Nevertheless, although the strains analyzed differed in their dynamics, we have not been able to link a particular genetic element with spreading success. Finally, the present study showed that new sequencing technologies improve considerably the quality and quantity of information obtained for a single strain; but this information is still difficult to interpret and important investments are required for the technology to become accessible for routine investigations.

Visualization and quality assessment of de novo genome assemblies

SUMMARY

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.

CONTACT

oksana.ribagrognuz@unil.ch; yannick.wurm@unil.ch

AVAILABILITY

Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.