Respiratory viruses dynamics and interactions: ten years of surveillance in central Europe.

BACKGROUND: Lower respiratory tract infections are among the main causes of death. Although there are many respiratory viruses, diagnostic efforts are focused mainly on influenza. The Respiratory Viruses Network (RespVir) collects infection data, primarily from German university hospitals, for a high diversity of infections by respiratory pathogens. In this study, we computationally analysed a subset of the RespVir database, covering 217,150 samples tested for 17 different viral pathogens in the time span from 2010 to 2019. METHODS: We calculated the prevalence of 17 respiratory viruses, analysed their seasonality patterns using information-theoretic measures and agglomerative clustering, and analysed their propensity for dual infection using a new metric dubbed average coinfection exclusion score (ACES). RESULTS: After initial data pre-processing, we retained 206,814 samples, corresponding to 1,408,657 performed tests. We found that Influenza viruses were reported for almost the half of all infections and that they exhibited the highest degree of seasonality. Coinfections of viruses are frequent; the most prevalent coinfection was rhinovirus/bocavirus and most of the virus pairs had a positive ACES indicating a tendency to exclude each other regarding infection. CONCLUSIONS: The analysis of respiratory viruses dynamics in monoinfection and coinfection contributes to the prevention, diagnostic, treatment, and development of new therapeutics. Data obtained from multiplex testing is fundamental for this analysis and should be prioritized over single pathogen testing.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq.

KEY MESSAGE: We studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode. MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

Evaluation of a Rapid Antigen Test To Detect SARS-CoV-2 Infection and Identify Potentially Infectious Individuals.

The identification and isolation of highly infectious SARS-CoV-2-infected individuals is an important public health strategy. Rapid antigen detection tests (RADT) are promising tools for large-scale screenings due to timely results and feasibility for on-site testing. Nonetheless, the diagnostic performance of RADT in detecting infectious individuals is not yet fully determined. In this study, RT-qPCR and virus culture of RT-qPCR-positive samples were used to evaluate and compare the performance of the Standard Q COVID-19 Ag test in detecting SARS-CoV-2-infected and possibly infectious individuals. To this end, two combined oro- and nasopharyngeal swabs were collected at a routine SARS-CoV-2 diagnostic center. A total of 2,028 samples were tested, and 118 virus cultures were inoculated. SARS-CoV-2 infection was detected in 210 samples by RT-qPCR, representing a positive rate of 10.36%. The Standard Q COVID-19 Ag test yielded a positive result in 92 (4.54%) samples resulting in an overall sensitivity and specificity of 42.86 and 99.89%, respectively. For adjusted CT values of <20 (n = 14), <25 (n = 57), and <30 (n = 88), the RADT reached sensitivities of 100, 98.25, and 88.64%, respectively. All 29 culture-positive samples were detected by the RADT. Although the overall sensitivity was low, the Standard Q COVID-19 Ag test reliably detected patients with high RNA loads. In addition, negative RADT results fully corresponded with the lack of viral cultivability in Vero E6 cells. These results indicate that RADT can be a valuable tool for the detection of individuals with high RNA loads that are likely to transmit SARS-CoV-2.

Weighted elastic net for unsupervised domain adaptation with application to age prediction from DNA methylation data.

MOTIVATION: Predictive models are a powerful tool for solving complex problems in computational biology. They are typically designed to predict or classify data coming from the same unknown distribution as the training data. In many real-world settings, however, uncontrolled biological or technical factors can lead to a distribution mismatch between datasets acquired at different times, causing model performance to deteriorate on new data. A common additional obstacle in computational biology is scarce data with many more features than samples. To address these problems, we propose a method for unsupervised domain adaptation that is based on a weighted elastic net. The key idea of our approach is to compare dependencies between inputs in training and test data and to increase the cost of differently behaving features in the elastic net regularization term. In doing so, we encourage the model to assign a higher importance to features that are robust and behave similarly across domains. RESULTS: We evaluate our method both on simulated data with varying degrees of distribution mismatch and on real data, considering the problem of age prediction based on DNA methylation data across multiple tissues. Compared with a non-adaptive standard model, our approach substantially reduces errors on samples with a mismatched distribution. On real data, we achieve far lower errors on cerebellum samples, a tissue which is not part of the training data and poorly predicted by standard models. Our results demonstrate that unsupervised domain adaptation is possible for applications in computational biology, even with many more features than samples. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/PfeiferLabTue/wenda. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Disentangling transcription factor binding site complexity.

The binding motifs of many transcription factors (TFs) comprise a higher degree of complexity than a single position weight matrix model permits. Additional complexity is typically taken into account either as intra-motif dependencies via more sophisticated probabilistic models or as heterogeneities via multiple weight matrices. However, both orthogonal approaches have limitations when learning from in vivo data where binding sites of other factors in close proximity can interfere with motif discovery for the protein of interest. In this work, we demonstrate how intra-motif complexity can, purely by analyzing the statistical properties of a given set of TF-binding sites, be distinguished from complexity arising from an intermix with motifs of co-binding TFs or other artifacts. In addition, we study the related question whether intra-motif complexity is represented more effectively by dependencies, heterogeneities or variants in between. Benchmarks demonstrate the effectiveness of both methods for their respective tasks and applications on motif discovery output from recent tools detect and correct many undesirable artifacts. These results further suggest that the prevalence of intra-motif dependencies may have been overestimated in previous studies on in vivo data and should thus be reassessed.

InMoDe: tools for learning and visualizing intra-motif dependencies of DNA binding sites.

Summary: Recent studies have shown that the traditional position weight matrix model is often insufficient for modeling transcription factor binding sites, as intra-motif dependencies play a significant role for an accurate description of binding motifs. Here, we present the Java application InMoDe, a collection of tools for learning, leveraging and visualizing such dependencies of putative higher order. The distinguishing feature of InMoDe is a robust model selection from a class of parsimonious models, taking into account dependencies only if justified by the data while choosing for simplicity otherwise. Availability and Implementation: InMoDe is implemented in Java and is available as command line application, as application with a graphical user-interface, and as an integration into Galaxy on the project website at http://www.jstacs.de/index.php/InMoDe . Contact: ralf.eggeling@cs.helsinki.fi.

Inferring intra-motif dependencies of DNA binding sites from ChIP-seq data.

BACKGROUND: Statistical modeling of transcription factor binding sites is one of the classical fields in bioinformatics. The position weight matrix (PWM) model, which assumes statistical independence among all nucleotides in a binding site, used to be the standard model for this task for more than three decades but its simple assumptions are increasingly put into question. Recent high-throughput sequencing methods have provided data sets of sufficient size and quality for studying the benefits of more complex models. However, learning more complex models typically entails the danger of overfitting, and while model classes that dynamically adapt the model complexity to data have been developed, effective model selection is to date only possible for fully observable data, but not, e.g., within de novo motif discovery. RESULTS: To address this issue, we propose a stochastic algorithm for performing robust model selection in a latent variable setting. This algorithm yields a solution without relying on hyperparameter-tuning via massive cross-validation or other computationally expensive resampling techniques. Using this algorithm for learning inhomogeneous parsimonious Markov models, we study the degree of putative higher-order intra-motif dependencies for transcription factor binding sites inferred via de novo motif discovery from ChIP-seq data. We find that intra-motif dependencies are prevalent and not limited to first-order dependencies among directly adjacent nucleotides, but that second-order models appear to be the significantly better choice. CONCLUSIONS: The traditional PWM model appears to be indeed insufficient to infer realistic sequence motifs, as it is on average outperformed by more complex models that take into account intra-motif dependencies. Moreover, using such models together with an appropriate model selection procedure does not lead to a significant performance loss in comparison with the PWM model for any of the studied transcription factors. Hence, we find it worthwhile to recommend that any modern motif discovery algorithm should attempt to take into account intra-motif dependencies.

Deep Learning in the Detection of Rare Fractures - Development of a "Deep Learning Convolutional Network" Model for Detecting Acetabular Fractures. / Deep Learning in der Erkennung seltener Frakturen ­ Entwicklung eines "Deep Convolutional Neural Network" (DCNN) am Beispiel der Azetabulumfraktur.

BACKGROUND: Fracture detection by artificial intelligence and especially Deep Convolutional Neural Networks (DCNN) is a topic of growing interest in current orthopaedic and radiological research. As learning a DCNN usually needs a large amount of training data, mostly frequent fractures as well as conventional X-ray are used. Therefore, less common fractures like acetabular fractures (AF) are underrepresented in the literature. The aim of this pilot study was to establish a DCNN for detection of AF using computer tomography (CT) scans. METHODS: Patients with an acetabular fracture were identified from the monocentric consecutive pelvic injury registry at the BG Trauma Center XXX from 01/2003 - 12/2019. All patients with unilateral AF and CT scans available in DICOM-format were included for further processing. All datasets were automatically anonymised and digitally post-processed. Extraction of the relevant region of interests was performed and the technique of data augmentation (DA) was implemented to artificially increase the number of training samples. A DCNN based on Med3D was used for autonomous fracture detection, using global average pooling (GAP) to reduce overfitting. RESULTS: From a total of 2,340 patients with a pelvic fracture, 654 patients suffered from an AF. After screening and post-processing of the datasets, a total of 159 datasets were enrolled for training of the algorithm. A random assignment into training datasets (80%) and test datasets (20%) was performed. The technique of bone area extraction, DA and GAP increased the accuracy of fracture detection from 58.8% (native DCNN) up to an accuracy of 82.8% despite the low number of datasets. CONCLUSION: The accuracy of fracture detection of our trained DCNN is comparable to published values despite the low number of training datasets. The techniques of bone extraction, DA and GAP are useful for increasing the detection rates of rare fractures by a DCNN. Based on the used DCNN in combination with the described techniques from this pilot study, the possibility of an automatic fracture classification of AF is under investigation in a multicentre study.

Dynamics and durability of HIV-1 neutralization are determined by viral replication.

Human immunodeficiency virus type 1 (HIV-1)-neutralizing antibodies (nAbs) that prevent infection are the main goal of HIV vaccine discovery. But as no nAb-eliciting vaccines are yet available, only data from HIV-1 neutralizers-persons with HIV-1 who naturally develop broad and potent nAbs-can inform about the dynamics and durability of nAb responses in humans, knowledge which is crucial for the design of future HIV-1 vaccine regimens. To address this, we assessed HIV-1-neutralizing immunoglobulin G (IgG) from 2,354 persons with HIV-1 on or off antiretroviral therapy (ART). Infection with non-clade B viruses, CD4+ T cell counts <200 µl-1, being off ART and a longer time off ART were independent predictors of a more potent and broad neutralization. In longitudinal analyses, we found nAb half-lives of 9.3 and 16.9 years in individuals with no- or low-level viremia, respectively, and 4.0 years in persons who newly initiated ART. Finally, in a potent HIV-1 neutralizer, we identified lower fractions of serum nAbs and of nAb-encoding memory B cells after ART initiation, suggesting that a decreasing neutralizing serum activity after antigen withdrawal is due to lower levels of nAbs. These results collectively show that HIV-1-neutralizing responses can persist for several years, even at low antigen levels, suggesting that an HIV-1 vaccine may elicit a durable nAb response.

Kinetics and correlates of the neutralizing antibody response to SARS-CoV-2 infection in humans.

Understanding antibody-based SARS-CoV-2 immunity is critical for overcoming the COVID-19 pandemic and informing vaccination strategies. We evaluated SARS-CoV-2 antibody dynamics over 10 months in 963 individuals who predominantly experienced mild COVID-19. Investigating 2,146 samples, we initially detected SARS-CoV-2 antibodies in 94.4% of individuals, with 82% and 79% exhibiting serum and IgG neutralization, respectively. Approximately 3% of individuals demonstrated exceptional SARS-CoV-2 neutralization, with these "elite neutralizers" also possessing SARS-CoV-1 cross-neutralizing IgG. Multivariate statistical modeling revealed age, symptomatic infection, disease severity, and gender as key factors predicting SARS-CoV-2-neutralizing activity. A loss of reactivity to the virus spike protein was observed in 13% of individuals 10 months after infection. Neutralizing activity had half-lives of 14.7 weeks in serum versus 31.4 weeks in purified IgG, indicating a rather long-term IgG antibody response. Our results demonstrate a broad spectrum in the initial SARS-CoV-2-neutralizing antibody response, with sustained antibodies in most individuals for 10 months after mild COVID-19.

Polyclonal and convergent antibody response to Ebola virus vaccine rVSV-ZEBOV.

Recombinant vesicular stomatitis virus-Zaire Ebola virus (rVSV-ZEBOV) is the most advanced Ebola virus vaccine candidate and is currently being used to combat the outbreak of Ebola virus disease (EVD) in the Democratic Republic of the Congo (DRC). Here we examine the humoral immune response in a subset of human volunteers enrolled in a phase 1 rVSV-ZEBOV vaccination trial by performing comprehensive single B cell and electron microscopy structure analyses. Four studied vaccinees show polyclonal, yet reproducible and convergent B cell responses with shared sequence characteristics. EBOV-targeting antibodies cross-react with other Ebolavirus species, and detailed epitope mapping revealed overlapping target epitopes with antibodies isolated from EVD survivors. Moreover, in all vaccinees, we detected highly potent EBOV-neutralizing antibodies with activities comparable or superior to the monoclonal antibodies currently used in clinical trials. These include antibodies combining the IGHV3-15/IGLV1-40 immunoglobulin gene segments that were identified in all investigated individuals. Our findings will help to evaluate and direct current and future vaccination strategies and offer opportunities for novel EVD therapies.

On the value of intra-motif dependencies of human insulator protein CTCF.

The binding affinity of DNA-binding proteins such as transcription factors is mainly determined by the base composition of the corresponding binding site on the DNA strand. Most proteins do not bind only a single sequence, but rather a set of sequences, which may be modeled by a sequence motif. Algorithms for de novo motif discovery differ in their promoter models, learning approaches, and other aspects, but typically use the statistically simple position weight matrix model for the motif, which assumes statistical independence among all nucleotides. However, there is no clear justification for that assumption, leading to an ongoing debate about the importance of modeling dependencies between nucleotides within binding sites. In the past, modeling statistical dependencies within binding sites has been hampered by the problem of limited data. With the rise of high-throughput technologies such as ChIP-seq, this situation has now changed, making it possible to make use of statistical dependencies effectively. In this work, we investigate the presence of statistical dependencies in binding sites of the human enhancer-blocking insulator protein CTCF by using the recently developed model class of inhomogeneous parsimonious Markov models, which is capable of modeling complex dependencies while avoiding overfitting. These findings lead to a more detailed characterization of the CTCF binding motif, which is only poorly represented by independent nucleotide frequencies at several positions, predominantly at the 3' end.