The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

JBrowse 2: a modular genome browser with views of synteny and structural variation.

We present JBrowse 2, a general-purpose genome annotation browser offering enhanced visualization of complex structural variation and evolutionary relationships. It retains core features of JBrowse while adding new views for synteny, dotplots, breakpoints, gene fusions, and whole-genome overviews. It allows users to share sessions, open multiple genomes, and navigate between views. It can be embedded in a web page, used as a standalone application, or run from Jupyter notebooks or R sessions. These improvements are enabled by a ground-up redesign using modern web technology. We describe application functionality, use cases, performance benchmarks, and implementation notes for web administrators and developers.

JBrowse Jupyter: a Python interface to JBrowse 2.

MOTIVATION: JBrowse Jupyter is a package that aims to close the gap between Python programming and genomic visualization. Web-based genome browsers are routinely used for publishing and inspecting genome annotations. Historically they have been deployed at the end of bioinformatics pipelines, typically decoupled from the analysis itself. However, emerging technologies such as Jupyter notebooks enable a more rapid iterative cycle of development, analysis and visualization. RESULTS: We have developed a package that provides a Python interface to JBrowse 2's suite of embeddable components, including the primary Linear Genome View. The package enables users to quickly set up, launch and customize JBrowse views from Jupyter notebooks. In addition, users can share their data via Google's Colab notebooks, providing reproducible interactive views. AVAILABILITY AND IMPLEMENTATION: JBrowse Jupyter is released under the Apache License and is available for download on PyPI. Source code and demos are available on GitHub at https://github.com/GMOD/jbrowse-jupyter.

Oligonucleotide-directed proximity-interactome mapping (O-MAP): A unified method for discovering RNA-interacting proteins, transcripts and genomic loci in situ.

Throughout biology, RNA molecules form complex networks of molecular interactions that are central to their function, but remain challenging to investigate. Here, we introduce Oligonucleotide-mediated proximity-interactome MAPping (O-MAP), a straightforward method for elucidating the biomolecules near an RNA of interest, within its native cellular context. O-MAP uses programmable oligonucleotide probes to deliver proximity-biotinylating enzymes to a target RNA, enabling nearby molecules to be enriched by streptavidin pulldown. O-MAP induces exceptionally precise RNA-localized in situ biotinylation, and unlike alternative methods it enables straightforward optimization of its targeting accuracy. Using the 47S pre-ribosomal RNA and long noncoding RNA Xist as models, we develop O-MAP workflows for unbiased discovery of RNA-proximal proteins, transcripts, and genomic loci. This revealed unexpected co-compartmentalization of Xist and other chromatin-regulatory RNAs and enabled systematic characterization of nucleolar-chromatin interactions across multiple cell lines. O-MAP is portable to cultured cells, organoids, and tissues, and to RNAs of various lengths, abundances, and sequence composition. And, O-MAP requires no genetic manipulation and uses exclusively off-the-shelf parts. We therefore anticipate its application to a broad array of RNA phenomena.

Cohesin controls X chromosome structure remodeling and X-reactivation during mouse iPSC-reprogramming.

Reactivation of the inactive X chromosome is a hallmark epigenetic event during reprogramming of mouse female somatic cells to induced pluripotent stem cells (iPSCs). This involves global structural remodeling from a condensed, heterochromatic into an open, euchromatic state, thereby changing a transcriptionally inactive into an active chromosome. Despite recent advances, very little is currently known about the molecular players mediating this process and how this relates to iPSC-reprogramming in general. To gain more insight, here we perform a RNAi-based knockdown screen during iPSC-reprogramming of mouse fibroblasts. We discover factors important for X chromosome reactivation (XCR) and iPSC-reprogramming. Among those, we identify the cohesin complex member SMC1a as a key molecule with a specific function in XCR, as its knockdown greatly affects XCR without interfering with iPSC-reprogramming. Using super-resolution microscopy, we find SMC1a to be preferentially enriched on the active compared with the inactive X chromosome and that SMC1a is critical for the decompacted state of the active X. Specifically, depletion of SMC1a leads to contraction of the active X both in differentiated and in pluripotent cells, where it normally is in its most open state. In summary, we reveal cohesin as a key factor for remodeling of the X chromosome from an inactive to an active structure and that this is a critical step for XCR during iPSC-reprogramming.

PaintSHOP enables the interactive design of transcriptome- and genome-scale oligonucleotide FISH experiments.

Fluorescence in situ hybridization (FISH) allows researchers to visualize the spatial position and quantity of nucleic acids in fixed samples. Recently, considerable progress has been made in developing oligonucleotide (oligo)-based FISH methods that have enabled researchers to study the three-dimensional organization of the genome at super-resolution and visualize the spatial patterns of gene expression for thousands of genes in individual cells. However, there are few existing computational tools to support the bioinformatics workflows necessary to carry out these experiments using oligo FISH probes. Here, we introduce paint server and homology optimization pipeline (PaintSHOP), an interactive platform for the design of oligo FISH experiments. PaintSHOP enables researchers to identify probes for their experimental targets efficiently, to incorporate additional necessary sequences such as primer pairs and to easily generate files documenting library design. PaintSHOP democratizes and standardizes the process of designing complex probe sets for the oligo FISH community.

JBrowseR: an R interface to the JBrowse 2 genome browser.

MOTIVATION: Genome browsers are an essential tool in genome analysis. Modern genome browsers enable complex and interactive visualization of a wide variety of genomic data modalities. While such browsers are very powerful, they can be challenging to configure and program for bioinformaticians lacking expertise in web development. RESULTS: We have developed an R package that provides an interface to the JBrowse 2 genome browser. The package can be used to configure and customize the browser entirely with R code. The browser can be deployed from the R console, or embedded in Shiny applications or R Markdown documents. AVAILABILITY AND IMPLEMENTATION: JBrowseR is available for download from CRAN, and the source code is openly available from the Github repository at https://github.com/GMOD/JBrowseR/.

3D mapping and accelerated super-resolution imaging of the human genome using in situ sequencing.

There is a need for methods that can image chromosomes with genome-wide coverage, as well as greater genomic and optical resolution. We introduce OligoFISSEQ, a suite of three methods that leverage fluorescence in situ sequencing (FISSEQ) of barcoded Oligopaint probes to enable the rapid visualization of many targeted genomic regions. Applying OligoFISSEQ to human diploid fibroblast cells, we show how four rounds of sequencing are sufficient to produce 3D maps of 36 genomic targets across six chromosomes in hundreds to thousands of cells, implying a potential to image thousands of targets in only five to eight rounds of sequencing. We also use OligoFISSEQ to trace chromosomes at finer resolution, following the path of the X chromosome through 46 regions, with separate studies showing compatibility of OligoFISSEQ with immunocytochemistry. Finally, we combined OligoFISSEQ with OligoSTORM, laying the foundation for accelerated single-molecule super-resolution imaging of large swaths of, if not entire, human genomes.

OligoMinerApp: a web-server application for the design of genome-scale oligonucleotide in situ hybridization probes through the flexible OligoMiner environment.

Fluorescence in situ hybridization (FISH) is a powerful single-cell technique that harnesses nucleic acid base pairing to detect the abundance and positioning of cellular RNA and DNA molecules in fixed samples. Recent technology development has paved the way to the construction of FISH probes entirely from synthetic oligonucleotides (oligos), allowing the optimization of thermodynamic properties together with the opportunity to design probes against any sequenced genome. However, comparatively little progress has been made in the development of computational tools to facilitate the oligos design, and even less has been done to extend their accessibility. OligoMiner is an open-source and modular pipeline written in Python that introduces a novel method of assessing probe specificity that employs supervised machine learning to predict probe binding specificity from genome-scale sequence alignment information. However, its use is restricted to only those people who are confident with command line interfaces because it lacks a Graphical User Interface (GUI), potentially cutting out many researchers from this technology. Here, we present OligoMinerApp (http://oligominerapp.org), a web-based application that aims to extend the OligoMiner framework through the implementation of a smart and easy-to-use GUI and the introduction of new functionalities specially designed to make effective probe mining available to everyone.

The Effect of Mouse Strain, Sex, and Carcinogen Dose on Toxicity and the Development of Lung Dysplasia and Squamous Cell Carcinomas in Mice.

In order to translate new treatments to the clinic, it is necessary to use animal models that closely recapitulate human disease. Lung cancer develops after extended exposure to carcinogens. It has one of the highest mutation rates of all cancer and is highly heterogenic. Topical treatment with N-nitrosotris-(2-chloroethyl)urea (NTCU) induces lung squamous cell carcinoma (SCC) with nonsynonymous mutation rates similar to those reported for human non-small cell lung cancer. However, NTCU induces lung cancer with variable efficacy and toxicity depending on the mouse strain. A detailed characterization of the NTCU model is needed. We have compared the effect of three different NTCU doses (20, 30, and 40 mmol/L) in female and male of NIH Swiss, Black Swiss, and FVB mice on tumor incidence, survival, and toxicity. The main findings in this study are (1) NIH Swiss mice present with a higher incidence of SCC and lower mortality compared with Black Swiss and FVB mice; (2) 30 mmol/L NTCU dose induces SCC at the same rate and incidence as the 40 mmol/L dose with lower mortality; (3) female mice present higher grade and incidence of preinvasive lesions and SCC compared with males; (4) NTCU-induced transformation is principally within the respiratory system; and (5) NTCU treatment does not affect the ability to elicit a specific adaptive immune response. This study provides a reference point for experimental designs to evaluate either preventive or therapeutic treatments for lung SCC, including immunotherapies, before initiating human clinical trials.