ANSWER: Annotation Software for Electronic Reporting

PURPOSE To better use genetic testing, which is used by clinicians to explain the molecular mechanism of disease and to suggest clinical actionability and new treatment options, clinical next-generation sequencing (NGS) laboratories must send the results into reports in PDF and discrete data element format (HL7). Although most clinical diagnostic tests have set molecular markers tested and have a set range of values or a binary result (positive or negative), the NGS genetic test could examine hundreds or thousands of genes with no predefined list of variants. Although there are some commercial and open-source tools for clinically reporting genomics results for oncology testing, they often lack necessary features. METHODS Using several available software tools for data storage including MySQL and MongoDB, database querying with Python, and a web-based user application using JAVA and JAVA script, we have developed a tool to store and query complex genomics and demographics data, which can be manually curated and reported by the user. RESULTS We have developed a tool, Annotation SoftWare for Electronic Reporting (ANSWER), that can allow molecular pathologists to (1) filter variants to find those meeting quality control metrics in the genes that are clinically actionable by diagnosis; (2) visualize variants using data generated in the bioinformatics analysis; (3) create annotations that can be reused in future reports with association specific to the gene, variant, or diagnosis; (4) select variants and annotations that should be reported to match the details of the case; and (5) generate a report that includes demographics, reported variants, clinical actionability annotation, and references that can be exported into PDF or HL7 format, which can be electronically sent to an electronic health record. CONCLUSION ANSWER is a tool that can be installed locally and is designed to meet the clinical reporting needs of a clinical oncology NGS laboratory for reporting.

Evaluating denoising strategies in resting-state fMRI in traumatic brain injury (EpiBioS4Rx)

Objective: Resting-state functional MRI is increasingly used in the clinical setting and is now included in some diagnostic guidelines for severe brain injury patients. However, to ensure high-quality data, one should mitigate fMRI-related noise typical of this population. Therefore, we aimed to evaluate the ability of different preprocessing strategies to mitigate noise-related signal (i.e., in-scanner movement and physiological noise) in functional connectivity of traumatic brain injury patients. Methods: We applied nine commonly used denoising strategies, combined into 17 pipelines, to 88 traumatic brain injury patients from the Epilepsy Bioinformatics Study for Anti-epileptogenic Therapy clinical trial (EpiBioS4Rx). Pipelines were evaluated by three quality control metrics across three exclusion regimes based on the participant's head movement profile. Results: While no pipeline eliminated noise effects on functional connectivity, some pipelines exhibited relatively high effectiveness depending on the exclusion regime. Once high-motion participants were excluded, the choice of denoising pipeline becomes secondary - although this strategy leads to substantial data loss. Pipelines combining spike regression with physiological regressors were the best performers, whereas pipelines that used automated data driven methods performed comparatively worse. Conclusion: In this study, we report the first large-scale evaluation of denoising pipelines aimed at reducing noise-related functional connectivity in a clinical population known to be highly susceptible to in-scanner motion and significant anatomical abnormalities. If resting-state functional magnetic resonance is to be a successful clinical technique, it is crucial that procedures mitigating the effect of noise be systematically evaluated in the most challenging populations, such as traumatic brain injury datasets.

Angiotensin converting enzyme 2 (ACE2) is expressed in murine cutaneous under single-cell transcriptome resolution

Angiotensin converting enzyme 2 (Ace2) is widely distributed in human organs, which was identified as a functional receptor for severe acute respiratory syndrome (SARS) coronavirus in human beings. It was also confirmed that SARS-CoV-2 uses the same cell entry receptor, ACE2, as SARS-CoV. However, related research still not discover the expression data associated with murine skin under single cell RNA resolution. In this study, we performed single-cell RNA sequencing (scRNA-seq) on unsorted cells from mouse dorsal skin after 7 days post-wounding. 8312 sequenced cells from four skin samples met quality control metrics and were analyzed.

Characterization of the consensus mucosal microbiome of colorectal cancer

Dysbioisis is an imbalance of an organ's microbiome and plays a role in colorectal cancer pathogenesis. Characterizing the bacteria in the microenvironment of a cancer through genome sequencing has advantages compared to culture-based profiling. However, there are notable technical and analytical challenges in characterizing universal features of tumor microbiomes. Colorectal tumors demonstrate microbiome variation among different studies and across individual patients. To address these issues, we conducted a computational study to determine a consensus microbiome for colorectal cancer, analyzing 924 tumors from eight independent RNA-Seq data sets. A standardized meta-transcriptomic analysis pipeline was established with quality control metrics. Microbiome profiles across different cohorts were compared and recurrently altered microbial shifts specific to colorectal cancer were determined. We identified cancer-specific set of 114 microbial species associated with tumors that were found among all investigated studies. Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria were among the four most abundant phyla for the colorectal cancer microbiome. Member species of Clostridia were depleted and Fusobacterium nucleatum was one of the most enriched bacterial species in tumors. Associations between the consensus species and specific immune cell types were noted. Our results are available as a web data resource for other researchers to explore (https://crc-microbiome.stanford.edu).

PEPATAC: an optimized pipeline for ATAC-seq data analysis with serial alignments

As chromatin accessibility data from ATAC-seq experiments continues to expand, there is continuing need for standardized analysis pipelines. Here, we present PEPATAC, an ATAC-seq pipeline that is easily applied to ATAC-seq projects of any size, from one-off experiments to large-scale sequencing projects. PEPATAC leverages unique features of ATAC-seq data to optimize for speed and accuracy, and it provides several unique analytical approaches. Output includes convenient quality control plots, summary statistics, and a variety of generally useful data formats to set the groundwork for subsequent project-specific data analysis. Downstream analysis is simplified by a standard definition format, modularity of components, and metadata APIs in R and Python. It is restartable, fault-tolerant, and can be run on local hardware, using any cluster resource manager, or in provided Linux containers. We also demonstrate the advantage of aligning to the mitochondrial genome serially, which improves the accuracy of alignment statistics and quality control metrics. PEPATAC is a robust and portable first step for any ATAC-seq project. BSD2-licensed code and documentation are available at https://pepatac.databio.org.

Proteomic profiling dataset of chemical perturbations in multiple biological backgrounds

While gene expression profiling has traditionally been the method of choice for large-scale perturbational profiling studies, proteomics has emerged as an effective tool in this context for directly monitoring cellular responses to perturbations. We previously reported a pilot library containing 3400 profiles of multiple perturbations across diverse cellular backgrounds in the reduced-representation phosphoproteome (P100) and chromatin space (Global Chromatin Profiling, GCP). Here, we expand our original dataset to include profiles from a new set of cardiotoxic compounds and from astrocytes, an additional neural cell model, totaling 5300 proteomic signatures. We describe filtering criteria and quality control metrics used to assess and validate the technical quality and reproducibility of our data. To demonstrate the power of the library, we present two case studies where data is queried using the concept of “connectivity” to obtain biological insight. All data presented in this study have been deposited to the ProteomeXchange Consortium with identifiers PXD017458 (P100) and PXD017459 (GCP) and can be queried at https://clue.io/proteomics.

CHIPS: A Snakemake pipeline for quality control and reproducible processing of chromatin profiling data

Motivation: The chromatin profile measured by ATAC-seq, ChIP-seq, or DNase-seq experiments can identify genomic regions critical in regulating gene expression and provide insights on biological processes such as diseases and development. However, quality control and processing chromatin profiling data involves many steps, and different bioinformatics tools are used at each step. It can be challenging to manage the analysis. Results: We developed a Snakemake pipeline called CHIPS (CHromatin enrIchment ProcesSor) to streamline the processing of ChIP-seq, ATAC-seq, and DNase-seq data. The pipeline supports single- and paired-end data and is flexible to start with FASTQ or BAM files. It includes basic steps such as read trimming, mapping, and peak calling. In addition, it calculates quality control metrics such as contamination profiles, polymerase chain reaction bottleneck coefficient, the fraction of reads in peaks, percentage of peaks overlapping with the union of public DNaseI hypersensitivity sites, and conservation profile of the peaks. For downstream analysis, it carries out peak annotations, motif finding, and regulatory potential calculation for all genes. The pipeline ensures that the processing is robust and reproducible. Availability: CHIPS is available at https://github.com/liulab-dfci/CHIPS.

CHIPS: A Snakemake pipeline for quality control and reproducible processing of chromatin profiling data

MotivationThe chromatin profile measured by ATAC-seq, ChIP-seq, or DNase-seq experiments can identify genomic regions critical in regulating gene expression and provide insights on biological processes such as diseases and development. However, quality control and processing chromatin profiling data involve many steps, and different bioinformatics tools are used at each step. It can be challenging to manage the analysis.ResultsWe developed a Snakemake pipeline called CHIPS (CHromatin enrichment Processor) to streamline the processing of ChIP-seq, ATAC-seq, and DNase-seq data. The pipeline supports single- and paired-end data and is flexible to start with FASTQ or BAM files. It includes basic steps such as read trimming, mapping, and peak calling. In addition, it calculates quality control metrics such as contamination profiles, PCR bottleneck coefficient, the fraction of reads in peaks, percentage of peaks overlapping with the union of public DNaseI hypersensitivity sites, and conservation profile of the peaks. For downstream analysis, it carries out peak annotations, motif finding, and regulatory potential calculation for all genes. The pipeline ensures that the processing is robust and reproducible.AvailabilityCHIPS is available at https://github.com/liulab-dfci/CHIPS

A multiplexed, next generation sequencing platform for high-throughput detection of SARS-CoV-2

Population scale sweeps of viral pathogens, such as SARS-CoV-2, require high intensity testing for effective management. Here, we describe “Systematic Parallel Analysis of RNA coupled to Sequencing for Covid-19 screening” (C19-SPAR-Seq), a multiplexed, scalable, readily automated platform for SARS-CoV-2 detection that is capable of analyzing tens of thousands of patient samples in a single run. To address strict requirements for control of assay parameters and output demanded by clinical diagnostics, we employ a control-based Precision-Recall and Receiver Operator Characteristics (coPR) analysis to assign run-specific quality control metrics. C19-SPAR-Seq coupled to coPR on a trial cohort of several hundred patients performs with a specificity of 100% and sensitivity of 91% on samples with low viral loads, and a sensitivity of >95% on high viral loads associated with disease onset and peak transmissibility. This study establishes the feasibility of employing C19-SPAR-Seq for the large-scale monitoring of SARS-CoV-2 and other pathogens.