Bacterial whole genome-based phylogeny: construction of a new benchmarking dataset and assessment of some existing methods.

BACKGROUND: Whole genome sequencing (WGS) is increasingly used in diagnostics and surveillance of infectious diseases. A major application for WGS is to use the data for identifying outbreak clusters, and there is therefore a need for methods that can accurately and efficiently infer phylogenies from sequencing reads. In the present study we describe a new dataset that we have created for the purpose of benchmarking such WGS-based methods for epidemiological data, and also present an analysis where we use the data to compare the performance of some current methods. RESULTS: Our aim was to create a benchmark data set that mimics sequencing data of the sort that might be collected during an outbreak of an infectious disease. This was achieved by letting an E. coli hypermutator strain grow in the lab for 8 consecutive days, each day splitting the culture in two while also collecting samples for sequencing. The result is a data set consisting of 101 whole genome sequences with known phylogenetic relationship. Among the sequenced samples 51 correspond to internal nodes in the phylogeny because they are ancestral, while the remaining 50 correspond to leaves. We also used the newly created data set to compare three different online available methods that infer phylogenies from whole-genome sequencing reads: NDtree, CSI Phylogeny and REALPHY. One complication when comparing the output of these methods with the known phylogeny is that phylogenetic methods typically build trees where all observed sequences are placed as leafs, even though some of them are in fact ancestral. We therefore devised a method for post processing the inferred trees by collapsing short branches (thus relocating some leafs to internal nodes), and also present two new measures of tree similarity that takes into account the identity of both internal and leaf nodes. CONCLUSIONS: Based on this analysis we find that, among the investigated methods, CSI Phylogeny had the best performance, correctly identifying 73% of all branches in the tree and 71% of all clades. We have made all data from this experiment (raw sequencing reads, consensus whole-genome sequences, as well as descriptions of the known phylogeny in a variety of formats) publicly available, with the hope that other groups may find this data useful for benchmarking and exploring the performance of epidemiological methods. All data is freely available at: https://cge.cbs.dtu.dk/services/evolution_data.php .

Geographical Distribution and Pattern of Pesticides in Danish Drinking Water 2002-2018: Reducing Data Complexity.

Pesticides are a large and heterogenous group of chemicals with a complex geographic distribution in the environment. The purpose of this study was to explore the geographic distribution of pesticides in Danish drinking water and identify potential patterns in the grouping of pesticides. Our data included 899,169 analyses of 167 pesticides and metabolites, of which 55 were identified above the detection limit. Pesticide patterns were defined by (1) pesticide groups based on chemical structure and pesticide-metabolite relations and (2) an exploratory factor analysis identifying underlying patterns of related pesticides within waterworks. The geographic distribution was evaluated by mapping the pesticide categories for groups and factor components, namely those detected, quantified, above quality standards, and not analysed. We identified five and seven factor components for the periods 2002-2011 and 2012-2018, respectively. In total, 16 pesticide groups were identified, of which six were representative in space and time with regards to the number of waterworks and analyses, namely benzothiazinone, benzonitriles, organophosphates, phenoxy herbicides, triazines, and triazinones. Pesticide mapping identified areas where multiple pesticides were detected, indicating areas with a higher pesticide burden. The results contribute to a better understanding of the pesticide pattern in Danish drinking water and may contribute to exposure assessments for future epidemiological studies.

An interactive database for the investigation of high-density peptide microarray guided interaction patterns and antivenom cross-reactivity.

Snakebite envenoming is a major neglected tropical disease that affects millions of people every year. The only effective treatment against snakebite envenoming consists of unspecified cocktails of polyclonal antibodies purified from the plasma of immunized production animals. Currently, little data exists on the molecular interactions between venom-toxin epitopes and antivenom-antibody paratopes. To address this issue, high-density peptide microarray (hdpm) technology has recently been adapted to the field of toxinology. However, analysis of such valuable datasets requires expert understanding and, thus, complicates its broad application within the field. In the present study, we developed a user-friendly, and high-throughput web application named "Snake Toxin and Antivenom Binding Profiles" (STAB Profiles), to allow straight-forward analysis of hdpm datasets. To test our tool and evaluate its performance with a large dataset, we conducted hdpm assays using all African snake toxin protein sequences available in the UniProt database at the time of study design, together with eight commercial antivenoms in clinical use in Africa, thus representing the largest venom-antivenom dataset to date. Furthermore, we introduced a novel method for evaluating raw signals from a peptide microarray experiment and a data normalization protocol enabling intra-microarray and even inter-microarray chip comparisons. Finally, these data, alongside all the data from previous similar studies by Engmark et al., were preprocessed according to our newly developed protocol and made publicly available for download through the STAB Profiles web application (http://tropicalpharmacology.com/tools/stab-profiles/). With these data and our tool, we were able to gain key insights into toxin-antivenom interactions and were able to differentiate the ability of different antivenoms to interact with certain toxins of interest. The data, as well as the web application, we present in this article should be of significant value to the venom-antivenom research community. Knowledge gained from our current and future analyses of this dataset carry the potential to guide the improvement and optimization of current antivenoms for maximum patient benefit, as well as aid the development of next-generation antivenoms.

Using Polygenic Scores in Social Science Research: Unraveling Childlessness.

Biological, genetic, and socio-demographic factors are all important in explaining reproductive behavior, yet these factors are typically studied in isolation. In this study, we explore an innovative sociogenomic approach, which entails including key socio-demographic (marriage, education, occupation, religion, cohort) and genetic factors related to both behavioral [age at first birth (AFB), number of children ever born (NEB)] and biological fecundity-related outcomes (endometriosis, age at menopause and menarche, polycystic ovary syndrome, azoospermia, testicular dysgenesis syndrome) to explain childlessness. We examine the association of all sets of factors with childlessness as well as the interplay between them. We derive polygenic scores (PGS) from recent genome-wide association studies (GWAS) and apply these in the Health and Retirement Study (N = 10,686) and Wisconsin Longitudinal Study (N = 8,284). Both socio-demographic and genetic factors were associated with childlessness. Whilst socio-demographic factors explain 19-46% in childlessness, the current PGS explains <1% of the variance, and only PGSs from large GWASs are related to childlessness. Our findings also indicate that genetic and socio-demographic factors are not independent, with PGSs for AFB and NEB related to education and age at marriage. The explained variance by polygenic scores on childlessness is limited since it is largely a behavioral trait, with genetic explanations expected to increase somewhat in the future with better-powered GWASs. As genotyping of individuals in social science surveys becomes more prevalent, the method described in this study can be applied to other outcomes.