Managing and monitoring a pandemic: showcasing a practical approach for the genomic surveillance of SARS-CoV-2.

With the rapidly growing amount of biological data, powerful but also flexible data management and visualization systems are of increasingly crucial importance. The COVID-19 pandemic has more than highlighted this need and the challenges scientists are facing. Here, we provide an example and a step-by-step template for non-IT personnel to easily implement an intuitive, interactive data management solution to manage and visualize the high influx of biological samples and associated metadata in a laboratory setting. Our approach is illustrated with the genomic surveillance for SARS-CoV-2 in Germany, covering over 11 600 internal and 130 000 external samples from multiple datasets. We compare three data management options used in laboratories: (i) simple, yet error-prone and inefficient spreadsheets, (ii) complex and long-to-implement laboratory information management systems and (iii) high-performance database management systems. We highlight the advantages and pitfalls of each option and outline why a document-oriented NoSQL option via MongoDB Atlas can be a suitable solution for many labs. Our example can be treated as a template and easily adapted to allow scientists to focus on their core work and not on complex data administration.

Faltwerk: a library for spatial exploratory data analysis of protein structures.

Summary: Proteins are fundamental building blocks of life and are investigated in a broad range of scientific fields, especially in the context of recent progress using in silico structure prediction models and the surge of resulting protein structures in public databases. However, exploratory data analysis of these proteins can be slow because of the need for several methods, ranging from geometric and spatial analysis to visualization. The Python library faltwerk provides an integrated toolkit to perform explorative work with rapid feedback. This toolkit includes support for protein complexes, spatial analysis (point density or spatial autocorrelation), ligand binding site prediction and an intuitive visualization interface based on the grammar of graphics. Availability and implementation: faltwerk is distributed under the permissive BSD-3 open source license. Source code and documentation, including an extensive common-use case tutorial, can be found at github.com/phiweger/faltwerk; binaries are available from the pypi repository.

Nanopore-based enrichment of antimicrobial resistance genes - a case-based study.

Rapid screening of hospital admissions to detect asymptomatic carriers of resistant bacteria can prevent pathogen outbreaks. However, the resulting isolates rarely have their genome sequenced due to cost constraints and long turn-around times to get and process the data, limiting their usefulness to the practitioner. Here we used real-time, on-device target enrichment ("adaptive") sequencing as a highly multiplexed assay covering 1,147 antimicrobial resistance genes. We compared its utility against standard and metagenomic sequencing, focusing on an isolate of Raoultella ornithinolytica harbouring three carbapenemases (NDM, KPC, VIM). Based on this experimental data, we then modelled the influence of several variables on the enrichment results and predicted the large effect of nucleotide identity (higher is better) and read length (shorter is better). Lastly, we showed how all relevant resistance genes are detected using adaptive sequencing on a miniature ("Flongle") flow cell, motivating its use in a clinical setting to monitor similar cases and their surroundings.

Enhanced glycerol assimilation and lipid production in Rhodotorula toruloides CBS14 upon addition of hemicellulose primarily correlates with early transcription of energy-metabolism-related genes.

BACKGROUND: Lipid formation from glycerol was previously found to be activated in Rhodotorula toruloides when the yeast was cultivated in a mixture of crude glycerol (CG) and hemicellulose hydrolysate (CGHH) compared to CG as the only carbon source. RNA samples from R. toruloides CBS14 cell cultures grown on either CG or CGHH were collected at different timepoints of cultivation, and a differential gene expression analysis was performed between cells grown at a similar physiological situation. RESULTS: We observed enhanced transcription of genes involved in oxidative phosphorylation and enzymes localized in mitochondria in CGHH compared to CG. Genes involved in protein turnover, including those encoding ribosomal proteins, translation elongation factors, and genes involved in building the proteasome also showed an enhanced transcription in CGHH compared to CG. At 10 h cultivation, another group of activated genes in CGHH was involved in ß-oxidation, handling oxidative stress and degradation of xylose and aromatic compounds. Potential bypasses of the standard GUT1 and GUT2-glycerol assimilation pathway were also expressed and upregulated in CGHH 10 h. When the additional carbon sources from HH were completely consumed, at CGHH 36 h, their transcription decreased and NAD+-dependent glycerol-3-phosphate dehydrogenase was upregulated compared to CG 60 h, generating NADH instead of NADPH with glycerol catabolism. TPI1 was upregulated in CGHH compared to cells grown on CG in all physiological situations, potentially channeling the DHAP formed through glycerol catabolism into glycolysis. The highest number of upregulated genes encoding glycolytic enzymes was found after 36 h in CGHH, when all additional carbon sources were already consumed. CONCLUSIONS: We suspect that the physiological reason for the accelerated glycerol assimilation and faster lipid production, was primarily the activation of enzymes that provide energy.

What the Phage: a scalable workflow for the identification and analysis of phage sequences.

Phages are among the most abundant and diverse biological entities on earth. Phage prediction from sequence data is a crucial first step to understanding their impact on the environment. A variety of bacteriophage prediction tools have been developed over the years. They differ in algorithmic approach, results, and ease of use. We, therefore, developed "What the Phage" (WtP), an easy-to-use and parallel multitool approach for phage prediction combined with an annotation and classification downstream strategy, thus supporting the user's decision-making process by summarizing the results of the different prediction tools in charts and tables. WtP is reproducible and scales to thousands of datasets through a workflow manager (Nextflow). WtP is freely available under a GPL-3.0 license (https://github.com/replikation/What_the_Phage).

Near Chromosome-Level Genome Assembly and Annotation of Rhodotorula babjevae Strains Reveals High Intraspecific Divergence.

The genus Rhodotorula includes basidiomycetous oleaginous yeast species. Rhodotorula babjevae can produce compounds of biotechnological interest such as lipids, carotenoids, and biosurfactants from low value substrates such as lignocellulose hydrolysate. High-quality genome assemblies are needed to develop genetic tools and to understand fungal evolution and genetics. Here, we combined short- and long-read sequencing to resolve the genomes of two R. babjevae strains, CBS 7808 (type strain) and DBVPG 8058, at chromosomal level. Both genomes are 21 Mbp in size and have a GC content of 68.2%. Allele frequency analysis indicates that both strains are tetraploid. The genomes consist of a maximum of 21 chromosomes with a size of 0.4 to 2.4 Mbp. In both assemblies, the mitochondrial genome was recovered in a single contig, that shared 97% pairwise identity. Pairwise identity between most chromosomes ranges from 82 to 87%. We also found indications for strain-specific extrachromosomal endogenous DNA. A total of 7591 and 7481 protein-coding genes were annotated in CBS 7808 and DBVPG 8058, respectively. CBS 7808 accumulated a higher number of tandem duplications than DBVPG 8058. We identified large translocation events between putative chromosomes. Genome divergence values between the two strains indicate that they may belong to different species.

Evaluation of microbiome enrichment and host DNA depletion in human vaginal samples using Oxford Nanopore's adaptive sequencing.

Metagenomic sequencing is promising for clinical applications to study microbial composition concerning disease or patient outcomes. Alterations of the vaginal microbiome are associated with adverse pregnancy outcomes, like preterm premature rupture of membranes and preterm birth. Methodologically these samples often have to deal with low relative amounts of prokaryotic DNA and high amounts of host DNA (> 90%), decreasing the overall microbial resolution. Nanopore's adaptive sampling method offers selective DNA depletion or target enrichment to directly reject or accept DNA molecules during sequencing without specialized sample preparation. Here, we demonstrate how selective 'human host depletion' resulted in a 1.70 fold (± 0.27 fold) increase in total sequencing depth, providing higher taxonomic profiling sensitivity. At the same time, the microbial composition remains consistent with the control experiments. The complete removal of all human host sequences is not yet possible and should be considered as an ethical approval statement might still be necessary. Adaptive sampling increased microbial sequencing yield in all 15 sequenced clinical routine vaginal samples, making it a valuable tool for clinical surveillance and medical-based research, which can be used in addition to other host depletion methods before sequencing.

Leclercia pneumoniae sp. nov., a bacterium isolated from clinical specimen in Leipzig, Germany.

Strain 49125T was isolated from an infant with pneumonia and septicaemia at the Leipzig University Hospital. Phenotypic and genomic traits were investigated. The strain's biochemical profile and its MALDI-TOF spectrogram did not differ from comparative samples of Leclercia adecarboxylata, thus far the sole member of the Leclercia species. A circular genome with a size of 4.4 Mbp and a G+C content of 55.0 mol% was reconstructed using hybrid Illumina and Nanopore sequencing. Phylogenetic analysis was based on 172 marker genes and validated using a k-mer-based search against a large genome collection including subsequent in silico DNA-DNA hybridization. Whole genome average nucleotide identity to any described species was below 95%, suggesting that strain 49125T represents a new species, for which we propose the name Leclercia pneumoniae sp. nov. with the type strain 49125T (=LMG 32245T=DSM 112336T).

Context-aware genomic surveillance reveals hidden transmission of a carbapenemase-producing Klebsiella pneumoniae.

Genomic surveillance can inform effective public health responses to pathogen outbreaks. However, integration of non-local data is rarely done. We investigate two large hospital outbreaks of a carbapenemase-carrying Klebsiella pneumoniae strain in Germany and show the value of contextual data. By screening about 10 000 genomes, over 400 000 metagenomes and two culture collections using in silico and in vitro methods, we identify a total of 415 closely related genomes reported in 28 studies. We identify the relationship between the two outbreaks through time-dated phylogeny, including their respective origin. One of the outbreaks presents extensive hidden transmission, with descendant isolates only identified in other studies. We then leverage the genome collection from this meta-analysis to identify genes under positive selection. We thereby identify an inner membrane transporter (ynjC) with a putative role in colistin resistance. Contextual data from other sources can thus enhance local genomic surveillance at multiple levels and should be integrated by default when available.

Chromosome-level genome assembly and transcriptome-based annotation of the oleaginous yeast Rhodotorula toruloides CBS 14.

Rhodotorula toruloides is an oleaginous yeast with high biotechnological potential. In order to understand the molecular physiology of lipid synthesis in R. toruloides and to advance metabolic engineering, a high-resolution genome is required. We constructed a genome draft of R. toruloides CBS 14, using a hybrid assembly approach, consisting of short and long reads generated by Illumina and Nanopore sequencing, respectively. The genome draft consists of 23 contigs and 3 scaffolds, with a N50 length of 1,529,952 bp, thus largely representing chromosomal organization. The total size of the genome is 20,534,857 bp and the overall GC content is 61.83%. Transcriptomic data from different growth conditions was used to aid species-specific gene annotation. We annotated 9464 genes and identified 11,691 transcripts. Furthermore, we demonstrated the presence of a potential plasmid, an extrachromosomal circular structure of about 11 kb with a copy number about three times as high as the other chromosomes.

Metagenomics workflow for hybrid assembly, differential coverage binning, metatranscriptomics and pathway analysis (MUFFIN).

Metagenomics has redefined many areas of microbiology. However, metagenome-assembled genomes (MAGs) are often fragmented, primarily when sequencing was performed with short reads. Recent long-read sequencing technologies promise to improve genome reconstruction. However, the integration of two different sequencing modalities makes downstream analyses complex. We, therefore, developed MUFFIN, a complete metagenomic workflow that uses short and long reads to produce high-quality bins and their annotations. The workflow is written by using Nextflow, a workflow orchestration software, to achieve high reproducibility and fast and straightforward use. This workflow also produces the taxonomic classification and KEGG pathways of the bins and can be further used for quantification and annotation by providing RNA-Seq data (optionally). We tested the workflow using twenty biogas reactor samples and assessed the capacity of MUFFIN to process and output relevant files needed to analyze the microbial community and their function. MUFFIN produces functional pathway predictions and, if provided de novo metatranscript annotations across the metagenomic sample and for each bin. MUFFIN is available on github under GNUv3 licence: https://github.com/RVanDamme/MUFFIN.

VIDHOP, viral host prediction with deep learning.

MOTIVATION: Zoonosis, the natural transmission of infections from animals to humans, is a far-reaching global problem. The recent outbreaks of Zikavirus, Ebolavirus and Coronavirus are examples of viral zoonosis, which occur more frequently due to globalization. In case of a virus outbreak, it is helpful to know which host organism was the original carrier of the virus to prevent further spreading of viral infection. Recent approaches aim to predict a viral host based on the viral genome, often in combination with the potential host genome and arbitrarily selected features. These methods are limited in the number of different hosts they can predict or the accuracy of the prediction. RESULTS: Here, we present a fast and accurate deep learning approach for viral host prediction, which is based on the viral genome sequence only. We tested our deep neural network (DNN) on three different virus species (influenza A virus, rabies lyssavirus and rotavirus A). We achieved for each virus species an AUC between 0.93 and 0.98, allowing highly accurate predictions while using only fractions (100-400 bp) of the viral genome sequences. We show that deep neural networks are suitable to predict the host of a virus, even with a limited amount of sequences and highly unbalanced available data. The trained DNNs are the core of our virus-host prediction tool VIrus Deep learning HOst Prediction (VIDHOP). VIDHOP also allows the user to train and use models for other viruses. AVAILABILITY AND IMPLEMENTATION: VIDHOP is freely available under https://github.com/flomock/vidhop. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hypoxia compromises the mitochondrial metabolism of Alzheimer's disease microglia via HIF1.

Genetic Alzheimer's disease (AD) risk factors associate with reduced defensive amyloid ß plaque-associated microglia (AßAM), but the contribution of modifiable AD risk factors to microglial dysfunction is unknown. In AD mouse models, we observe concomitant activation of the hypoxia-inducible factor 1 (HIF1) pathway and transcription of mitochondrial-related genes in AßAM, and elongation of mitochondria, a cellular response to maintain aerobic respiration under low nutrient and oxygen conditions. Overactivation of HIF1 induces microglial quiescence in cellulo, with lower mitochondrial respiration and proliferation. In vivo, overstabilization of HIF1, either genetically or by exposure to systemic hypoxia, reduces AßAM clustering and proliferation and increases Aß neuropathology. In the human AD hippocampus, upregulation of HIF1α and HIF1 target genes correlates with reduced Aß plaque microglial coverage and an increase of Aß plaque-associated neuropathology. Thus, hypoxia (a modifiable AD risk factor) hijacks microglial mitochondrial metabolism and converges with genetic susceptibility to cause AD microglial dysfunction.

Assessing genetic diversity and similarity of 435 KPC-carrying plasmids.

The global spread and diversification of multidrug-resistant Gram-negative (MRGN) bacteria poses major challenges to healthcare. In particular, carbapenem-resistant Klebsiella pneumoniae strains have been frequently identified in infections and hospital-wide outbreaks. The most frequently underlying resistance gene (blaKPC) has been spreading over the last decade in the health care setting. blaKPC seems to have rapidly diversified and has been found in various species and on different plasmid types. To review the progress and dynamics of this diversification, all currently available KPC plasmids in the NCBI database were analysed in this work. Plasmids were grouped into 257 different representative KPC plasmids, of which 79.4% could be clearly assigned to incompatibility (Inc) group or groups. In almost half of all representative plasmids, the KPC gene is located on Tn4401 variants, emphasizing the importance of this transposon type for the transmission of KPC genes to other plasmids. The transposons also seem to be responsible for the occurrence of altered or uncommon fused plasmid types probably due to incomplete transposition. Moreover, many KPC plasmids contain genes that encode proteins promoting recombinant processes and mutagenesis; in consequence accelerating the diversification of KPC genes and other colocalized resistance genes.

Direct RNA nanopore sequencing of full-length coronavirus genomes provides novel insights into structural variants and enables modification analysis.

Sequence analyses of RNA virus genomes remain challenging owing to the exceptional genetic plasticity of these viruses. Because of high mutation and recombination rates, genome replication by viral RNA-dependent RNA polymerases leads to populations of closely related viruses, so-called "quasispecies." Standard (short-read) sequencing technologies are ill-suited to reconstruct large numbers of full-length haplotypes of (1) RNA virus genomes and (2) subgenome-length (sg) RNAs composed of noncontiguous genome regions. Here, we used a full-length, direct RNA sequencing (DRS) approach based on nanopores to characterize viral RNAs produced in cells infected with a human coronavirus. By using DRS, we were able to map the longest (∼26-kb) contiguous read to the viral reference genome. By combining Illumina and Oxford Nanopore sequencing, we reconstructed a highly accurate consensus sequence of the human coronavirus (HCoV)-229E genome (27.3 kb). Furthermore, by using long reads that did not require an assembly step, we were able to identify, in infected cells, diverse and novel HCoV-229E sg RNAs that remain to be characterized. Also, the DRS approach, which circumvents reverse transcription and amplification of RNA, allowed us to detect methylation sites in viral RNAs. Our work paves the way for haplotype-based analyses of viral quasispecies by showing the feasibility of intra-sample haplotype separation. Even though several technical challenges remain to be addressed to exploit the potential of the nanopore technology fully, our work illustrates that DRS may significantly advance genomic studies of complex virus populations, including predictions on long-range interactions in individual full-length viral RNA haplotypes.

Retroclival arachnoid cyst in a preterm infant after ventriculitis and intraventricular hemorrhage-a case report.

DISCUSSION: We present a case of a preterm boy of 29 gestational weeks with complicated posthemorrhagic hydrocephalus who developed a retroclival cyst extending to the upper cervical spinal canal and leading to brainstem compression. At the age of 2 months, third ventriculocisternocystostomy (ETVCC) resulted in temporary relief but was followed by reclosure of the stoma within weeks. At the age of 4 months, navigated endoscopic ventriculocisternocystostomy and endoscopy-guided stent catheter placement was performed with connection to a ventriculoperitoneal (VP) shunt system through a burr hole reservoir. CONCLUSION: This strategy improved the multiloculated hydrocephalus by establishing a communication between theretroclival subarachnoid cyst and the ventricular system, leading to uneventful further development of the child.

The Gini coefficient: a methodological pilot study to assess fetal brain development employing postmortem diffusion MRI.

BACKGROUND: Diffusion-weighted imaging (DWI) is important in the assessment of fetal brain development. However, it is clinically challenging and time-consuming to prepare neuromorphological examinations to assess real brain age and to detect abnormalities. OBJECTIVE: To demonstrate that the Gini coefficient can be a simple, intuitive parameter for modelling fetal brain development. MATERIALS AND METHODS: Postmortem fetal specimens(n = 28) were evaluated by diffusion-weighted imaging (DWI) on a 3-T MRI scanner using 60 directions, 0.7-mm isotropic voxels and b-values of 0, 150, 1,600 s/mm(2). Constrained spherical deconvolution (CSD) was used as the local diffusion model. Fractional anisotropy (FA), apparent diffusion coefficient (ADC) and complexity (CX) maps were generated. CX was defined as a novel diffusion metric. On the basis of those three parameters, the Gini coefficient was calculated. RESULTS: Study of fetal brain development in postmortem specimens was feasible using DWI. The Gini coefficient could be calculated for the combination of the three diffusion parameters. This multidimensional Gini coefficient correlated well with age (Adjusted R(2) = 0.59) between the ages of 17 and 26 gestational weeks. CONCLUSIONS: We propose a new method that uses an economics concept, the Gini coefficient, to describe the whole brain with one simple and intuitive measure, which can be used to assess the brain's developmental state.

Potential Pediatric Applications of PET/MR.

Medical imaging with multimodality and whole-body technologies has continuously improved in recent years. The advent of combined modalities such as PET/CT and PET/MR offers new tools with an exact fusion of molecular imaging and high-resolution anatomic imaging. For noninvasive pediatric diagnostics, molecular imaging and whole-body MR have become important, especially in pediatric oncology. Because it has a lower radiation exposure than PET/CT, combined PET/MR is expected to be of special use in pediatric diagnostics. This review focuses on possible pediatric applications of PET/MR hybrid imaging, particularly pediatric oncology and neurology but also the diagnosis of infectious or inflammatory diseases.

Sirolimus treatment of severe PTEN hamartoma tumor syndrome: case report and in vitro studies.

BACKGROUND: Phosphatase and tensin homolog (PTEN) hamartoma tumor syndrome (PHTS) is caused by germ line mutations in the PTEN gene. Symptoms include cancer predisposition, immune deviations, and lipomas/lipomatosis. No causal standard therapy is available. We describe a therapeutic attempt with the mammalian target of rapamycin (mTOR) inhibitor sirolimus for a PHTS patient suffering from thymus hyperplasia and lipomatosis. We furthermore assessed the in vitro effects of sirolimus and other inhibitors on lipoma cells of the patient. METHODS: The patient underwent clinical and blood examinations and whole-body magnetic resonance imaging to assess tumor sizes. Lipoma cells of the patient were incubated with inhibitors of the phosphoinositide-3-kinase (PI3K)/AKT/mTOR signaling pathway to analyze the effects on proliferation, adipocyte differentiation, and survival in vitro. RESULTS: Sirolimus treatment improved somatic growth and reduced thymus volume. These effects diminished over the treatment period of 19 mo. Sirolimus decreased lipoma cell proliferation and adipocyte differentiation in vitro but did not cause apoptosis. PI3K and AKT inhibitors induced apoptosis significantly. CONCLUSION: Sirolimus treatment led to an improvement of the patient's clinical status and a transient reduction of the thymus. Our in vitro findings point to PI3K and AKT inhibitors as potential treatment options for patients with severe forms of PHTS.