Integration of multiple epigenomic marks improves prediction of variant impact in saturation mutagenesis reporter assay

The integrative analysis of highâ throughput reporter assays, machine learning, and profiles of epigenomic chromatin state in a broad array of cells and tissues has the potential to significantly improve our understanding of noncoding regulatory element function and its contribution to human disease. Here, we report results from the CAGI 5 regulation saturation challenge where participants were asked to predict the impact of nucleotide substitution at every base pair within five diseaseâ associated human enhancers and nine diseaseâ associated promoters. A library of mutations covering all bases was generated by saturation mutagenesis and altered activity was assessed in a massively parallel reporter assay (MPRA) in relevant cell lines. Reporter expression was measured relative to plasmid DNA to determine the impact of variants. The challenge was to predict the functional effects of variants on reporter expression. Comparative analysis of the full range of submitted prediction results identifies the most successful models of transcription factor binding sites, machine learning algorithms, and ways to choose among or incorporate diverse datatypes and cellâ types for training computational models. These results have the potential to improve the design of future studies on more diverse sets of regulatory elements and aid the interpretation of diseaseâ associated genetic variation.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151884/1/humu23797_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151884/2/humu23797.pd

Genome-Wide Meta-Analysis Validates a Role for S1PR1 in Microtubule Targeting Agent-Induced Sensory Peripheral Neuropathy

Microtubule targeting agents (MTAs) are anticancer therapies commonly prescribed for breast cancer and other solid tumors. Sensory peripheral neuropathy (PN) is the major dose-limiting toxicity for MTAs and can limit clinical efficacy. The current pharmacogenomic study aimed to identify genetic variations that explain patient susceptibility and drive mechanisms underlying development of MTA-induced PN. A meta-analysis of genome-wide association studies (GWAS) from two clinical cohorts treated with MTAs (CALGB 40502 and CALGB 40101) was conducted using a Cox regression model with cumulative dose to first instance of grade 2 or higher PN. Summary statistics from a GWAS of European subjects (n = 469) in CALGB 40502 that estimated cause-specific risk of PN were meta-analyzed with those from a previously published GWAS of European ancestry (n = 855) from CALGB 40101 that estimated the risk of PN. Novel single nucleotide polymorphisms in an enhancer region downstream of sphingosine-1-phosphate receptor 1 (S1PR1 encoding S1PR1; e.g., rs74497159, βCALGB 40101 per allele log hazard ratio (95% CI) = 0.591 (0.254 – 0.928), βCALGB 40502 per allele log hazard ratio (95% CI) = 0.693 (0.334 – 1.053); PMETA = 3.62×10−7) were the most highly ranked associations based on P-values with risk of developing grade 2 and higher PN. In silico functional analysis identified multiple regulatory elements and potential enhancer activity for S1PR1 within this genomic region. Inhibition of S1PR1 function in iPSC-derived human sensory neurons shows partial protection against paclitaxel-induced neurite damage. These pharmacogenetic findings further support ongoing clinical evaluations to target S1PR1 as a therapeutic strategy for prevention and/or treatment of MTA-induced neuropathy

Risk-Association of DNMT1 Gene Polymorphisms with Coronary Artery Disease in Chinese Han Population

Recently, a significant epigenetic component in the pathogenesis of Coronary Artery Disease (CAD) has been realized. Here, we evaluated the possible association of candidate Single Nucleotide Polymorphisms (SNPs) in the epigenetic-regulatory gene, DNA methyltransferase 1 (DNMT1), with CAD in Chinese Han population. Five tag SNPs (rs16999593, rs2336691, rs2228611, rs4804494, rs7253062) were analyzed by High Resolution Melt (HRM) method in 476 CAD patients and 478 controls. Overall, there were significant differences in the genotype and allele distributions of rs2228611 and rs2336691, between patients and controls. The minor A allele of rs2228611 was associated with a lower risk of CAD (p = 0.034); modest effect in the additive analysis but also marginal significance was found in the recessive model [ORadditive = 0.404 (0.184, 0.884), p = 0.023 and ORrecessive = 0.452 (0.213, 0.963), p = 0.040] after adjusting for confounders. While the rs2336691 A allele were associated with a higher risk of developing CAD (p = 0.037); borderline significant association in both additive and dominant models [ORadditive = 1.632 (1.030, 2.583), p = 0.037 and ORdominant = 1.599 (1.020, 2.507), p = 0.040]. In conclusion, these data provide the first evidence that occurrence of CAD may be moderated by genetic variation in the gene involved in the epigenetic machinery

Deletion of Pr130 Interrupts Cardiac Development in Zebrafish

Protein phosphatase 2 regulatory subunit B, alpha (PPP2R3A), a regulatory subunit of protein phosphatase 2A (PP2A), is a major serine/threonine phosphatase that regulates crucial function in development and growth. Previous research has implied that PPP2R3A was involved in heart failure, and PR130, the largest transcription of PPP2R3A, functioning in the calcium release of sarcoplasmic reticulum (SR), plays an important role in the excitation-contraction (EC) coupling. To obtain a better understanding of PR130 functions in myocardium and cardiac development, two pr130-deletion zebrafish lines were generated using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated proteins (Cas) system. Pr130-knockout zebrafish exhibited cardiac looping defects and decreased cardiac function (decreased fractional area and fractional shortening). Hematoxylin and eosin (H&E) staining demonstrated reduced cardiomyocytes. Subsequent transmission electron microscopy revealed that the bright and dark bands were narrowed and blurred, the Z- and M-lines were fogged, and the gaps between longitudinal myocardial fibers were increased. Additionally, increased apoptosis was observed in cardiomyocyte in pr130-knockout zebrafish compared to wild-type (WT). Taken together, our results suggest that pr130 is required for normal myocardium formation and efficient cardiac contractile function

Saturation mutagenesis of twenty disease-associated regulatory elements at single base-pair resolution.

The majority of common variants associated with common diseases, as well as an unknown proportion of causal mutations for rare diseases, fall in noncoding regions of the genome. Although catalogs of noncoding regulatory elements are steadily improving, we have a limited understanding of the functional effects of mutations within them. Here, we perform saturation mutagenesis in conjunction with massively parallel reporter assays on 20 disease-associated gene promoters and enhancers, generating functional measurements for over 30,000 single nucleotide substitutions and deletions. We find that the density of putative transcription factor binding sites varies widely between regulatory elements, as does the extent to which evolutionary conservation or integrative scores predict functional effects. These data provide a powerful resource for interpreting the pathogenicity of clinically observed mutations in these disease-associated regulatory elements, and comprise a rich dataset for the further development of algorithms that aim to predict the regulatory effects of noncoding mutations

P‐Glycoprotein Inhibition Exacerbates Paclitaxel Neurotoxicity in Neurons and Patients With Cancer

Paclitaxel-induced peripheral neuropathy (PIPN) is a common and dose-limiting adverse event. The role of P-glycoprotein (P-gp) in the neuronal efflux of paclitaxel was assessed using a translational approach. SH-SY5Y cells were differentiated to neurons and paclitaxel toxicity in the absence and presence of a P-gp inhibitor was determined. Paclitaxel caused marked dose-dependent toxicity in SH-SY5Y-derived neurons. Paclitaxel neurotoxicity was exacerbated with concomitant P-gp inhibition by valspodar and verapamil, consistent with increased intracellular accumulation of paclitaxel. Patients with cancer treated with paclitaxel and P-gp inhibitors had a 2.4-fold (95% confidence interval (CI) 1.3-4.3) increased risk of peripheral neuropathy-induced dose modification while a 4.7-fold (95% CI 1.9-11.9) increased risk for patients treated with strong P-gp inhibitors was observed, and a 7.0-fold (95% CI 2.3-21.5) increased risk in patients treated with atorvastatin. Atorvastatin also increased neurotoxicity by paclitaxel in SH-SY5Y-derived neurons. Clinicians should be aware that comedication with P-gp inhibitors may lead to increased risk of PIPN

Deletion of CTCF sites in the SHH locus alters enhancer-promoter interactions and leads to acheiropodia.

Acheiropodia, congenital limb truncation, is associated with homozygous deletions in the LMBR1 gene around ZRS, an enhancer regulating SHH during limb development. How these deletions lead to this phenotype is unknown. Using whole-genome sequencing, we fine-mapped the acheiropodia-associated region to 12 kb and show that it does not function as an enhancer. CTCF and RAD21 ChIP-seq together with 4C-seq and DNA FISH identify three CTCF sites within the acheiropodia-deleted region that mediate the interaction between the ZRS and the SHH promoter. This interaction is substituted with other CTCF sites centromeric to the ZRS in the disease state. Mouse knockouts of the orthologous 12 kb sequence have no apparent abnormalities, showcasing the challenges in modelling CTCF alterations in animal models due to inherent motif differences between species. Our results show that alterations in CTCF motifs can lead to a Mendelian condition due to altered enhancer-promoter interactions

Integration of multiple epigenomic marks improves prediction of variant impact in saturation mutagenesis reporter assay

The integrative analysis of high-throughput reporter assays, machine learning, and profiles of epigenomic chromatin state in a broad array of cells and tissues has the potential to significantly improve our understanding of noncoding regulatory element function and its contribution to human disease. Here, we report results from the CAGI 5 regulation saturation challenge where participants were asked to predict the impact of nucleotide substitution at every base pair within five disease-associated human enhancers and nine disease-associated promoters. A library of mutations covering all bases was generated by saturation mutagenesis and altered activity was assessed in a massively parallel reporter assay (MPRA) in relevant cell lines. Reporter expression was measured relative to plasmid DNA to determine the impact of variants. The challenge was to predict the functional effects of variants on reporter expression. Comparative analysis of the full range of submitted prediction results identifies the most successful models of transcription factor binding sites, machine learning algorithms, and ways to choose among or incorporate diverse datatypes and cell-types for training computational models. These results have the potential to improve the design of future studies on more diverse sets of regulatory elements and aid the interpretation of disease-associated genetic variation