Packaging and containerization of computational methods.

Methods for analyzing the full complement of a biomolecule type, e.g., proteomics or metabolomics, generate large amounts of complex data. The software tools used to analyze omics data have reshaped the landscape of modern biology and become an essential component of biomedical research. These tools are themselves quite complex and often require the installation of other supporting software, libraries and/or databases. A researcher may also be using multiple different tools that require different versions of the same supporting materials. The increasing dependence of biomedical scientists on these powerful tools creates a need for easier installation and greater usability. Packaging and containerization are different approaches to satisfy this need by delivering omics tools already wrapped in additional software that makes the tools easier to install and use. In this systematic review, we describe and compare the features of prominent packaging and containerization platforms. We outline the challenges, advantages and limitations of each approach and some of the most widely used platforms from the perspectives of users, software developers and system administrators. We also propose principles to make the distribution of omics software more sustainable and robust to increase the reproducibility of biomedical and life science research.

Integration of 168,000 samples reveals global patterns of the human gut microbiome.

Understanding the factors that shape variation in the human microbiome is a major goal of research in biology. While other genomics fields have used large, pre-compiled compendia to extract systematic insights requiring otherwise impractical sample sizes, there has been no comparable resource for the 16S rRNA sequencing data commonly used to quantify microbiome composition. To help close this gap, we have assembled a set of 168,484 publicly available human gut microbiome samples, processed with a single pipeline and combined into the largest unified microbiome dataset to date. We use this resource, which is freely available at microbiomap.org, to shed light on global variation in the human gut microbiome. We find that Firmicutes, particularly Bacilli and Clostridia, are almost universally present in the human gut. At the same time, the relative abundance of the 65 most common microbial genera differ between at least two world regions. We also show that gut microbiomes in undersampled world regions, such as Central and Southern Asia, differ significantly from the more thoroughly characterized microbiomes of Europe and Northern America. Moreover, humans in these overlooked regions likely harbor hundreds of taxa that have not yet been discovered due to this undersampling, highlighting the need for diversity in microbiome studies. We anticipate that this new compendium can serve the community and enable advanced applied and methodological research.

High-throughput Oligopaint screen identifies druggable 3D genome regulators.

The human genome functions as a three-dimensional chromatin polymer, driven by a complex collection of chromosome interactions1-3. Although the molecular rules governing these interactions are being quickly elucidated, relatively few proteins regulating this process have been identified. Here, to address this gap, we developed high-throughput DNA or RNA labelling with optimized Oligopaints (HiDRO)-an automated imaging pipeline that enables the quantitative measurement of chromatin interactions in single cells across thousands of samples. By screening the human druggable genome, we identified more than 300 factors that influence genome folding during interphase. Among these, 43 genes were validated as either increasing or decreasing interactions between topologically associating domains. Our findings show that genetic or chemical inhibition of the ubiquitous kinase GSK3A leads to increased long-range chromatin looping interactions in a genome-wide and cohesin-dependent manner. These results demonstrate the importance of GSK3A signalling in nuclear architecture and the use of HiDRO for identifying mechanisms of spatial genome organization.

An atlas of lamina-associated chromatin across twelve human cell types reveals an intermediate chromatin subtype.

BACKGROUND: Association of chromatin with lamin proteins at the nuclear periphery has emerged as a potential mechanism to coordinate cell type-specific gene expression and maintain cellular identity via gene silencing. Unlike many histone modifications and chromatin-associated proteins, lamina-associated domains (LADs) are mapped genome-wide in relatively few genetically normal human cell types, which limits our understanding of the role peripheral chromatin plays in development and disease. RESULTS: To address this gap, we map LAMIN B1 occupancy across twelve human cell types encompassing pluripotent stem cells, intermediate progenitors, and differentiated cells from all three germ layers. Integrative analyses of this atlas with gene expression and repressive histone modification maps reveal that lamina-associated chromatin in all twelve cell types is organized into at least two subtypes defined by differences in LAMIN B1 occupancy, gene expression, chromatin accessibility, transposable elements, replication timing, and radial positioning. Imaging of fluorescently labeled DNA in single cells validates these subtypes and shows radial positioning of LADs with higher LAMIN B1 occupancy and heterochromatic histone modifications primarily embedded within the lamina. In contrast, the second subtype of lamina-associated chromatin is relatively gene dense, accessible, dynamic across development, and positioned adjacent to the lamina. Most genes gain or lose LAMIN B1 occupancy consistent with cell types along developmental trajectories; however, we also identify examples where the enhancer, but not the gene body and promoter, changes LAD state. CONCLUSIONS: Altogether, this atlas represents the largest resource to date for peripheral chromatin organization studies and reveals an intermediate chromatin subtype.

Co-depletion of NIPBL and WAPL balance cohesin activity to correct gene misexpression.

The relationship between cohesin-mediated chromatin looping and gene expression remains unclear. NIPBL and WAPL are two opposing regulators of cohesin activity; depletion of either is associated with changes in both chromatin folding and transcription across a wide range of cell types. However, a direct comparison of their individual and combined effects on gene expression in the same cell type is lacking. We find that NIPBL or WAPL depletion in human HCT116 cells each alter the expression of ~2,000 genes, with only ~30% of the genes shared between the conditions. We find that clusters of differentially expressed genes within the same topologically associated domain (TAD) show coordinated misexpression, suggesting some genomic domains are especially sensitive to both more or less cohesin. Finally, co-depletion of NIPBL and WAPL restores the majority of gene misexpression as compared to either knockdown alone. A similar set of NIPBL-sensitive genes are rescued following CTCF co-depletion. Together, this indicates that altered transcription due to reduced cohesin activity can be functionally offset by removal of either its negative regulator (WAPL) or the physical barriers (CTCF) that restrict loop-extrusion events.

BiomeHorizon: Visualizing Microbiome Time Series Data in R.

Despite playing a key role in the health of their hosts, host-associated microbial communities demonstrate considerable variation over time. These communities comprise thousands of temporally dynamic taxa, which makes visualizing microbial time series data challenging. As such, a method to visualize both the proportional and absolute change in the relative abundance of multiple taxa across multiple subjects over time is needed. To address this gap, we developed BiomeHorizon, the first automated, open-source R package that visualizes longitudinal compositional microbiome data using horizon plots. BiomeHorizon is available at https://github.com/blekhmanlab/biomehorizon/ and a guide with step-by-step instructions for using the package is provided at https://blekhmanlab.github.io/biomehorizon/. IMPORTANCE Host-associated microbiota (i.e., the number and types of bacteria in the body) can have profound impacts on an animal's day-to-day functioning as well as their long-term health. Recent work has shown that these microbial communities change substantially over time, so it is important to be able to link changes in the abundance of certain microbes with host health outcomes. However, visualizing such changes is difficult because the microbiome comprises thousands of different microbes. To address this issue, we developed BiomeHorizon, an R package for visualizing longitudinal microbiome data using horizon plots. BiomeHorizon accepts a range of data formats and was developed with two common microbiome study designs in mind: human health studies, where the microbiome is sampled at set time points, and observational wildlife studies, where samples may be collected at irregular time intervals. BiomeHorizon thus provides a flexible, user-friendly approach to microbiome time series data visualization and analysis.

Public human microbiome data are dominated by highly developed countries.

The importance of sampling from globally representative populations has been well established in human genomics. In human microbiome research, however, we lack a full understanding of the global distribution of sampling in research studies. This information is crucial to better understand global patterns of microbiome-associated diseases and to extend the health benefits of this research to all populations. Here, we analyze the country of origin of all 444,829 human microbiome samples that are available from the world's 3 largest genomic data repositories, including the Sequence Read Archive (SRA). The samples are from 2,592 studies of 19 body sites, including 220,017 samples of the gut microbiome. We show that more than 71% of samples with a known origin come from Europe, the United States, and Canada, including 46.8% from the US alone, despite the country representing only 4.3% of the global population. We also find that central and southern Asia is the most underrepresented region: Countries such as India, Pakistan, and Bangladesh account for more than a quarter of the world population but make up only 1.8% of human microbiome samples. These results demonstrate a critical need to ensure more global representation of participants in microbiome studies.

Comparing quality of reporting between preprints and peer-reviewed articles in the biomedical literature.

BACKGROUND: Preprint usage is growing rapidly in the life sciences; however, questions remain on the relative quality of preprints when compared to published articles. An objective dimension of quality that is readily measurable is completeness of reporting, as transparency can improve the reader's ability to independently interpret data and reproduce findings. METHODS: In this observational study, we initially compared independent samples of articles published in bioRxiv and in PubMed-indexed journals in 2016 using a quality of reporting questionnaire. After that, we performed paired comparisons between preprints from bioRxiv to their own peer-reviewed versions in journals. RESULTS: Peer-reviewed articles had, on average, higher quality of reporting than preprints, although the difference was small, with absolute differences of 5.0% [95% CI 1.4, 8.6] and 4.7% [95% CI 2.4, 7.0] of reported items in the independent samples and paired sample comparison, respectively. There were larger differences favoring peer-reviewed articles in subjective ratings of how clearly titles and abstracts presented the main findings and how easy it was to locate relevant reporting information. Changes in reporting from preprints to peer-reviewed versions did not correlate with the impact factor of the publication venue or with the time lag from bioRxiv to journal publication. CONCLUSIONS: Our results suggest that, on average, publication in a peer-reviewed journal is associated with improvement in quality of reporting. They also show that quality of reporting in preprints in the life sciences is within a similar range as that of peer-reviewed articles, albeit slightly lower on average, supporting the idea that preprints should be considered valid scientific contributions.

International authorship and collaboration across bioRxiv preprints.

Preprints are becoming well established in the life sciences, but relatively little is known about the demographics of the researchers who post preprints and those who do not, or about the collaborations between preprint authors. Here, based on an analysis of 67,885 preprints posted on bioRxiv, we find that some countries, notably the United States and the United Kingdom, are overrepresented on bioRxiv relative to their overall scientific output, while other countries (including China, Russia, and Turkey) show lower levels of bioRxiv adoption. We also describe a set of 'contributor countries' (including Uganda, Croatia and Thailand): researchers from these countries appear almost exclusively as non-senior authors on international collaborations. Lastly, we find multiple journals that publish a disproportionate number of preprints from some countries, a dynamic that almost always benefits manuscripts from the US.

Challenges and recommendations to improve the installability and archival stability of omics computational tools.

Developing new software tools for analysis of large-scale biological data is a key component of advancing modern biomedical research. Scientific reproduction of published findings requires running computational tools on data generated by such studies, yet little attention is presently allocated to the installability and archival stability of computational software tools. Scientific journals require data and code sharing, but none currently require authors to guarantee the continuing functionality of newly published tools. We have estimated the archival stability of computational biology software tools by performing an empirical analysis of the internet presence for 36,702 omics software resources published from 2005 to 2017. We found that almost 28% of all resources are currently not accessible through uniform resource locators (URLs) published in the paper they first appeared in. Among the 98 software tools selected for our installability test, 51% were deemed "easy to install," and 28% of the tools failed to be installed at all because of problems in the implementation. Moreover, for papers introducing new software, we found that the number of citations significantly increased when authors provided an easy installation process. We propose for incorporation into journal policy several practical solutions for increasing the widespread installability and archival stability of published bioinformatics software.

Rxivist.org: Sorting biology preprints using social media and readership metrics.

Preprints have arrived. In increasing numbers, researchers across the life sciences are embracing the once-niche practice, shaking off decades of reluctance and posting hundreds of papers per week to preprint servers, sharing their findings with the community before embarking on the weary march through peer review. However, there are limited methods for individuals sifting through this avalanche of research to identify the preprints that are most relevant to their interests. Here, we describe Rxivist.org, a website that indexes all preprints posted to bioRxiv.org, the largest preprint server in the life sciences, and allows users to filter and sort papers based on download metrics and Twitter activity over a variety of categories and time periods. In this work, we hope to make it easier for readers to find relevant research on bioRxiv and to improve the visibility of preprints currently being read and discussed online.

Tracking the popularity and outcomes of all bioRxiv preprints.

The growth of preprints in the life sciences has been reported widely and is driving policy changes for journals and funders, but little quantitative information has been published about preprint usage. Here, we report how we collected and analyzed data on all 37,648 preprints uploaded to bioRxiv.org, the largest biology-focused preprint server, in its first five years. The rate of preprint uploads to bioRxiv continues to grow (exceeding 2,100 in October 2018), as does the number of downloads (1.1 million in October 2018). We also find that two-thirds of preprints posted before 2017 were later published in peer-reviewed journals, and find a relationship between the number of downloads a preprint has received and the impact factor of the journal in which it is published. We also describe Rxivist.org, a web application that provides multiple ways to interact with preprint metadata.