Search

Epigenomics, the study of the human genome and its interactions with proteins and other cellular elements, has become of significant interest in recent years. Such interactions have been shown to regulate essential cellular functions and are associated with multiple complex diseases. Therefore, understanding how these interactions may change across conditions is central in biomedical research. Chromatin immunoprecipitation followed by massively parallel sequencing (ChIP‐seq) is one of several techniques to detect local changes in epigenomic activity (peaks). However, existing methods for differential peak calling are not optimized for the diversity in ChIP‐seq signal profiles, are limited to the analysis of two conditions, or cannot classify specific patterns of differential change when multiple patterns exist. To address these limitations, we present a flexible and efficient method for the detection of differential epigenomic activity across multiple conditions. We utilize data from the ENCODE Consortium and show that the presented method, epigraHMM, exhibits superior performance to current tools and it is among the fastest algorithms available, while allowing the classification of combinatorial patterns of differential epigenomic activity and the characterization of chromatin regulatory states. [ABSTRACT FROM AUTHOR]

Regression calibration is the most widely used method to adjust regression parameter estimates for covariate measurement error. Yet its application in the context of a complex sampling design, for which the common bootstrap variance estimator can be less straightforward, has been less studied. We propose 2 variance estimators for a multistage probability-based sampling design, a parametric and a resampling-based multiple imputation approach, where a latent mean exposure needed for regression calibration is the target of imputation. This work was motivated by the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) data from 2008 to 2011, for which relationships between several outcomes and diet, an error-prone self-reported exposure, are of interest. We assessed the relative performance of these variance estimation strategies in an extensive simulation study built on the HCHS/SOL data. We further illustrate the proposed estimators with an analysis of the cross-sectional association of dietary sodium intake with hypertension-related outcomes in a subsample of the HCHS/SOL cohort. We have provided guidelines for the application of regression models with regression-calibrated exposures. Practical considerations for implementation of these 2 variance estimators in the setting of a large multicenter study are also discussed. Code to replicate the presented results is available online. [ABSTRACT FROM AUTHOR]

Chromatin immunoprecipitation followed by next‐generation sequencing (ChIP‐seq) is a technique to detect genomic regions containing protein‐DNA interaction, such as transcription factor binding sites or regions containing histone modifications. One goal of the analysis of ChIP‐seq experiments is to identify genomic loci enriched for sequencing reads pertaining to DNA bound to the factor of interest. The accurate identification of such regions aids in the understanding of epigenomic marks and gene regulatory mechanisms. Given the reduction of massively parallel sequencing costs, methods to detect consensus regions of enrichment across multiple samples are of interest. Here, we present a statistical model to detect broad consensus regions of enrichment from ChIP‐seq technical or biological replicates through a class of zero‐inflated mixed‐effects hidden Markov models. We show that the proposed model outperforms existing methods for consensus peak calling in common epigenomic marks by accounting for the excess zeros and sample‐specific biases. We apply our method to data from the Encyclopedia of DNA Elements and Roadmap Epigenomics projects and also from an extensive simulation study. [ABSTRACT FROM AUTHOR]

People living with HIV on antiretroviral therapy often have undetectable virus levels by standard assays, but "latent" HIV still persists in viral reservoirs. Eliminating these reservoirs is the goal of HIV cure research. The quantitative viral outgrowth assay (QVOA) is commonly used to estimate the reservoir size, that is, the infectious units per million (IUPM) of HIV-persistent resting CD4+ T cells. A new variation of the QVOA, the ultra deep sequencing assay of the outgrowth virus (UDSA), was recently developed that further quantifies the number of viral lineages within a subset of infected wells. Performing the UDSA on a subset of wells provides additional information that can improve IUPM estimation. This paper considers statistical inference about the IUPM from combined dilution assay (QVOA) and deep viral sequencing (UDSA) data, even when some deep sequencing data are missing. Methods are proposed to accommodate assays with wells sequenced at multiple dilution levels and with imperfect sensitivity and specificity, and a novel bias-corrected estimator is included for small samples. The proposed methods are evaluated in a simulation study, applied to data from the University of North Carolina HIV Cure Center, and implemented in the open-source R package SLDeepAssay. [ABSTRACT FROM AUTHOR]