Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Acral lamellar ichthyosis with amino acid substitution in the C-terminus of keratin 2.

BACKGROUND: Most cases of hereditary ichthyoses present with generalized scaling and skin dryness. However, in some cases skin involvement is restricted to particular body regions as in acral lamellar ichthyosis. OBJECTIVES: We report on the genetic basis of acral ichthyosis in two families presenting with a similar phenotype. METHODS: Genetic testing was performed by targeted next generation sequencing and whole-exome sequencing. For identity-by-descent analysis, the parents were genotyped and data analysis was performed with the Chromosome Analysis Suite Software. RT-PCR with RNA extracted from skin samples was used to analyse the effect of variants on splicing. RESULTS: Genetic testing identified a few heterozygous variants, but only the variant in KRT2 c.1912 T > C, p.Phe638Leu segregated with the disease and remained the strongest candidate. Pairwise identity-by-descent analysis revealed no indication of family relationship. Phenylalanine 638 is the second last amino acid upstream of the termination codon in the tail of K2, and substitution to leucine is predicted as probably damaging. Assessment of the variant is difficult, in part due to the lack of crystal structures of this region. CONCLUSIONS: Altogether, we show that a type of autosomal dominant acral ichthyosis is most probably caused by an amino acid substitution in the C-terminus of keratin 2.

A Resilience Related Glial-Neurovascular Network Is Transcriptionally Activated after Chronic Social Defeat in Male Mice.

Upon chronic stress, a fraction of individuals shows stress resilience, which can prevent long-term mental dysfunction. The underlying molecular mechanisms are complex and have not yet been fully understood. In this study, we performed a data-driven behavioural stratification together with single-cell transcriptomics of the hippocampus in a mouse model of chronic social defeat stress. Our work revealed that in a sub-group exhibiting molecular responses upon chronic stress, the dorsal hippocampus is particularly involved in neuroimmune responses, angiogenesis, myelination, and neurogenesis, thereby enabling brain restoration and homeostasis after chronic stress. Based on these molecular insights, we applied rapamycin after the stress as a proof-of-concept pharmacological intervention and were able to substantially increase stress resilience. Our findings serve as a data resource and can open new avenues for further understanding of molecular processes underlying stress response and for targeted interventions supporting resilience.

Dichotomic Hippocampal Transcriptome After Glutamatergic vs. GABAergic Deletion of the Cannabinoid CB1 Receptor.

Brain homeostasis is the dynamic equilibrium whereby physiological parameters are kept actively within a specific range. The homeostatic range is not fixed and may change throughout the individual's lifespan, or may be transiently modified in the presence of severe perturbations. The endocannabinoid system has emerged as a safeguard of homeostasis, e.g., it modulates neurotransmission and protects neurons from prolonged or excessively strong activation. We used genetically engineered mouse lines that lack the cannabinoid type-1 receptor (CB1) either in dorsal telencephalic glutamatergic or in forebrain GABAergic neurons to create new allostatic states, resulting from alterations in the excitatory/inhibitory (E/I) balance. Previous studies with these two mouse lines have shown dichotomic results in the context of behavior, neuronal morphology, and electrophysiology. Thus, we aimed at analyzing the transcriptomic profile of the hippocampal CA region from these mice in the basal condition and after a mild behavioral stimulation (open field). Our results provide insights into the gene networks that compensate chronic E/I imbalances. Among these, there are differentially expressed genes involved in neuronal and synaptic functions, synaptic plasticity, and the regulation of behavior. Interestingly, some of these genes, e.g., Rab3b, Crhbp, and Kcnn2, and related pathways showed a dichotomic expression, i.e., they are up-regulated in one mutant line and down-regulated in the other one. Subsequent interrogation on the source of the alterations at transcript level were applied using exon-intron split analysis. However, no strong directions toward transcriptional or post-transcriptional regulation comparing both mouse lines were observed. Altogether, the dichotomic gene expression observed and their involved signaling pathways are of interest because they may act as "switches" to modulate the directionality of neural homeostasis, which then is relevant for pathologies, such as stress-related disorders and epilepsy.

Reliability of genomic variants across different next-generation sequencing platforms and bioinformatic processing pipelines.

BACKGROUND: Next Generation Sequencing (NGS) is the fundament of various studies, providing insights into questions from biology and medicine. Nevertheless, integrating data from different experimental backgrounds can introduce strong biases. In order to methodically investigate the magnitude of systematic errors in single nucleotide variant calls, we performed a cross-sectional observational study on a genomic cohort of 99 subjects each sequenced via (i) Illumina HiSeq X, (ii) Illumina HiSeq, and (iii) Complete Genomics and processed with the respective bioinformatic pipeline. We also repeated variant calling for the Illumina cohorts with GATK, which allowed us to investigate the effect of the bioinformatics analysis strategy separately from the sequencing platform's impact. RESULTS: The number of detected variants/variant classes per individual was highly dependent on the experimental setup. We observed a statistically significant overrepresentation of variants uniquely called by a single setup, indicating potential systematic biases. Insertion/deletion polymorphisms (indels) were associated with decreased concordance compared to single nucleotide polymorphisms (SNPs). The discrepancies in indel absolute numbers were particularly prominent in introns, Alu elements, simple repeats, and regions with medium GC content. Notably, reprocessing sequencing data following the best practice recommendations of GATK considerably improved concordance between the respective setups. CONCLUSION: We provide empirical evidence of systematic heterogeneity in variant calls between alternative experimental and data analysis setups. Furthermore, our results demonstrate the benefit of reprocessing genomic data with harmonized pipelines when integrating data from different studies.

Co-Inference of Data Mislabelings Reveals Improved Models in Genomics and Breast Cancer Diagnostics.

Mislabeling of cases as well as controls in case-control studies is a frequent source of strong bias in prognostic and diagnostic tests and algorithms. Common data processing methods available to the researchers in the biomedical community do not allow for consistent and robust treatment of labeled data in the situations where both, the case and the control groups, contain a non-negligible proportion of mislabeled data instances. This is an especially prominent issue in studies regarding late-onset conditions, where individuals who may convert to cases may populate the control group, and for screening studies that often have high false-positive/-negative rates. To address this problem, we propose a method for a simultaneous robust inference of Lasso reduced discriminative models and of latent group-specific mislabeling risks, not requiring any exactly labeled data. We apply it to a standard breast cancer imaging dataset and infer the mislabeling probabilities (being rates of false-negative and false-positive core-needle biopsies) together with a small set of simple diagnostic rules, outperforming the state-of-the-art BI-RADS diagnostics on these data. The inferred mislabeling rates for breast cancer biopsies agree with the published purely empirical studies. Applying the method to human genomic data from a healthy-ageing cohort reveals a previously unreported compact combination of single-nucleotide polymorphisms that are strongly associated with a healthy-ageing phenotype for Caucasians. It determines that 7.5% of Caucasians in the 1000 Genomes dataset (selected as a control group) carry a pattern characteristic of healthy ageing.

5' tRNA halves are highly expressed in the primate hippocampus and might sequence-specifically regulate gene expression.

Fragments of mature tRNAs have long been considered as mere degradation products without physiological function. However, recent reports show that tRNA-derived small RNAs (tsRNAs) play prominent roles in diverse cellular processes across a wide spectrum of species. Contrasting the situation in other small RNA pathways the mechanisms behind these effects appear more diverse, more complex, and are generally less well understood. In addition, surprisingly little is known about the expression profiles of tsRNAs across different tissues and species. Here, we provide an initial overview of tsRNA expression in different species and tissues, revealing very high levels of 5' tRNA halves (5' tRHs) particularly in the primate hippocampus. We further modulated the regulation capacity of selected 5' tRHs in human cells by transfecting synthetic tsRNA mimics ("overexpression") or antisense-RNAs ("inhibition") and identified differentially expressed transcripts based on RNA-seq. We then used a novel k-mer mapping approach to dissect the underlying targeting rules, suggesting that 5' tRHs silence genes in a sequence-specific manner, while the most efficient target sites align to the mid-region of the 5' tRH and are located within the CDS or 3' UTR of the target. This amends previous observations that tsRNAs guide Argonaute proteins to silence their targets via a miRNA-like 5' seed match and suggests a yet unknown mechanism of regulation. Finally, our data suggest that some 5' tRHs that are also able to sequence-specifically stabilize mRNAs as up-regulated mRNAs are also significantly enriched for 5' tRH target sites.

Common miRNA Patterns of Alzheimer's Disease and Parkinson's Disease and Their Putative Impact on Commensal Gut Microbiota.

With the rise of Next-Generation-Sequencing (NGS) methods, Micro-RNAs (miRNAs) have achieved an important position in the research landscape and have been found to present valuable diagnostic tools in various diseases such as multiple sclerosis or lung cancer. There is also emerging evidence that miRNAs play an important role in the pathogenesis of neurodegenerative diseases such as Alzheimer's disease (AD) or Parkinson's disease (PD). Apparently, these diseases come along with changes in miRNA expression patterns which led to attempts from researchers to use these small RNA species from several body fluids for a better diagnosis and in order to observe disease progression. Additionally, it became evident that microbial commensals might play an important role for pathology development and were shown to have a significantly different composition in patients suffering from neurodegeneration compared with healthy controls. As it could recently be shown that secreted miRNAs are able to enter microbial organisms, it is conceivable that the host's miRNA might affect the gut microbial ecosystem. As such, miRNAs may inherit a central role in shaping the "diseased microbiome" and thereby mutually act on the characteristics of these neurodegenerative diseases. We have therefore (1) compiled a list of miRNAs known to be associated with AD and/or PD, (2) performed an in silico target screen for binding sites of these miRNA on human gut metagenome sequences and (3) evaluated the hit list for interesting matches potentially relevant to the etiology of AD and or PD. The examination of protein identifiers connected to bacterial secretion system, lipopolysaccharide biosynthesis and biofilm formation revealed an overlap of 37 bacterial proteins that were targeted by human miRNAs. The identified links of miRNAs to the biological processes of bacteria connected to AD and PD have yet to be validated via in vivo experiments. However, our results show a promising new approach for understanding aspects of these neurodegenerative diseases in light of the regulation of the microbiome.

PIWI genes and piRNAs are ubiquitously expressed in mollusks and show patterns of lineage-specific adaptation.

PIWI proteins and PIWI-interacting RNAs (piRNAs) suppress transposon activity in animals, thus protecting their genomes from detrimental insertion mutagenesis. Here, we reveal that PIWI genes and piRNAs are ubiquitously expressed in mollusks, similar to the situation in arthropods. We describe lineage-specific adaptations of transposon composition in piRNA clusters in the great pond snail and the pacific oyster, likely reflecting differential transposon activity in gastropods and bivalves. We further show that different piRNA clusters with unique transposon composition are dynamically expressed during oyster development. Finally, bioinformatics analyses suggest that different populations of piRNAs presumably bound to different PIWI paralogs participate in homotypic and heterotypic ping-pong amplification loops in a tissue- and sex-specific manner. Together with recent findings from other animal species, our results support the idea that somatic piRNA expression represents the ancestral state in metazoans.

unitas: the universal tool for annotation of small RNAs.

BACKGROUND: Next generation sequencing is a key technique in small RNA biology research that has led to the discovery of functionally different classes of small non-coding RNAs in the past years. However, reliable annotation of the extensive amounts of small non-coding RNA data produced by high-throughput sequencing is time-consuming and requires robust bioinformatics expertise. Moreover, existing tools have a number of shortcomings including a lack of sensitivity under certain conditions, limited number of supported species or detectable sub-classes of small RNAs. RESULTS: Here we introduce unitas, an out-of-the-box ready software for complete annotation of small RNA sequence datasets, supporting the wide range of species for which non-coding RNA reference sequences are available in the Ensembl databases (currently more than 800). unitas combines high quality annotation and numerous analysis features in a user-friendly manner. A complete annotation can be started with one simple shell command, making unitas particularly useful for researchers not having access to a bioinformatics facility. Noteworthy, the algorithms implemented in unitas are on par or even outperform comparable existing tools for small RNA annotation that map to publicly available ncRNA databases. CONCLUSIONS: unitas brings together annotation and analysis features that hitherto required the installation of numerous different bioinformatics tools which can pose a challenge for the non-expert user. With this, unitas overcomes the problem of read normalization. Moreover, the high quality of sequence annotation and analysis, paired with the ease of use, make unitas a valuable tool for researchers in all fields connected to small RNA biology.

Temperature-responsive miRNAs in Drosophila orchestrate adaptation to different ambient temperatures.

The majority of Drosophila genes are expressed in a temperature-dependent manner, but the way in which small RNAs may contribute to this effect is completely unknown as we currently lack an idea of how small RNA transcriptomes change as a function of temperature. Applying high-throughput sequencing techniques complemented by quantitative real-time PCR experiments, we demonstrate that altered ambient temperature induces drastic but reversible changes in sequence composition and total abundance of both miRNA and piRNA populations. Further, mRNA sequencing reveals that the expression of miRNAs and their predicted target transcripts correlates inversely, suggesting that temperature-responsive miRNAs drive adaptation to different ambient temperatures on the transcriptome level. Finally, we demonstrate that shifts in temperature affect both primary and secondary piRNA pools, and the observed aberrations are consistent with altered expression levels of the involved Piwi-pathway factors. We further reason that enhanced ping-pong processing at 29°C is driven by dissolved RNA secondary structures at higher temperatures, uncovering target sites that are not accessible at low temperatures. Together, our results show that small RNAs are an important part of epigenetic regulatory mechanisms that ensure homeostasis and adaptation under fluctuating environmental conditions.