Changes in spatial chromatin interactions are now emerging as a unifying mechanism orchestrating the regulation of gene expression. Hi-C sequencing technology allows insight into chromatin ...interactions on a genome-wide scale. However, Hi-C data contains many DNA sequence- and technology-driven biases. These biases prevent effective comparison of chromatin interactions aimed at identifying genomic regions differentially interacting between, e.g., disease-normal states or different cell types. Several methods have been developed for normalizing individual Hi-C datasets. However, they fail to account for biases between two or more Hi-C datasets, hindering comparative analysis of chromatin interactions.
We developed a simple and effective method, HiCcompare, for the joint normalization and differential analysis of multiple Hi-C datasets. The method introduces a distance-centric analysis and visualization of the differences between two Hi-C datasets on a single plot that allows for a data-driven normalization of biases using locally weighted linear regression (loess). HiCcompare outperforms methods for normalizing individual Hi-C datasets and methods for differential analysis (diffHiC, FIND) in detecting a priori known chromatin interaction differences while preserving the detection of genomic structures, such as A/B compartments.
HiCcompare is able to remove between-dataset bias present in Hi-C matrices. It also provides a user-friendly tool to allow the scientific community to perform direct comparisons between the growing number of pre-processed Hi-C datasets available at online repositories. HiCcompare is freely available as a Bioconductor R package https://bioconductor.org/packages/HiCcompare/ .
Chronic alcohol abuse has been linked to the disruption of executive function and allostatic conditioning of reward response dysregulation in the mesocorticolimbic pathway (MCL). Here, we analyzed ...genome-wide mRNA and miRNA expression from matched cases with alcohol dependence (AD) and controls (n = 35) via gene network analysis to identify unique and shared biological processes dysregulated in the prefrontal cortex (PFC) and nucleus accumbens (NAc). We further investigated potential mRNA/miRNA interactions at the network and individual gene expression levels to identify the neurobiological mechanisms underlying AD in the brain. By using genotyped and imputed SNP data, we identified expression quantitative trait loci (eQTL) uncovering potential genetic regulatory elements for gene networks associated with AD. At a Bonferroni corrected p≤0.05, we identified significant mRNA (NAc = 6; PFC = 3) and miRNA (NAc = 3; PFC = 2) AD modules. The gene-set enrichment analyses revealed modules preserved between PFC and NAc to be enriched for immune response processes, whereas genes involved in cellular morphogenesis/localization and cilia-based cell projection were enriched in NAc modules only. At a Bonferroni corrected p≤0.05, we identified significant mRNA/miRNA network module correlations (NAc = 6; PFC = 4), which at an individual transcript level implicated miR-449a/b as potential regulators for cellular morphogenesis/localization in NAc. Finally, we identified eQTLs (NAc: mRNA = 37, miRNA = 9; PFC: mRNA = 17, miRNA = 16) which potentially mediate alcohol's effect in a brain region-specific manner. Our study highlights the neurotoxic effects of chronic alcohol abuse as well as brain region specific molecular changes that may impact the development of alcohol addiction.
Although epigenome-wide association studies (EWAS) have been successful in identifying DNA methylation (DNAm) patterns associated with disease states, any further characterization of etiologic ...mechanisms underlying disease remains elusive. This knowledge gap does not originate from a lack of DNAm-trait associations, but rather stems from study design issues that affect the interpretability of EWAS results. Despite known limitations in predicting the function of a particular CpG site, most EWAS maintain the broad assumption that altered DNAm results in a concomitant change of transcription at the most proximal gene. This study integrated DNAm and gene expression (GE) measurements in two cohorts, the Adolescent and Young Adult Twin Study (AYATS) and the Pregnancy, Race, Environment, Genes (PREG) study, to improve the understanding of epigenomic regulatory mechanisms. CpG sites associated with GE in cis were enriched in areas of transcription factor binding and areas of intermediate-to-low CpG density. CpG sites associated with trans GE were also enriched in areas of known regulatory significance, including enhancer regions. These results highlight issues with restricting DNAm-transcript annotations to small genomic intervals and question the validity of assuming a cis DNAm-GE pathway. Based on these findings, the interpretation of EWAS results is limited in studies without multi-omic support and further research should identify genomic regions in which GE-associated DNAm is overrepresented. An in-depth characterization of GE-associated CpG sites could improve predictions of the downstream functional impact of altered DNAm and inform best practices for interpreting DNAm-trait associations generated by EWAS.
Neuropsychiatric and substance use disorders (NPSUDs) have a complex etiology that includes environmental and polygenic risk factors with significant cross-trait genetic correlations. Genome-wide ...association studies (GWAS) of NPSUDs yield numerous association signals. However, for most of these regions, we do not yet have a firm understanding of either the specific risk variants or the effects of these variants. Post-GWAS methods allow researchers to use GWAS summary statistics and molecular mediators (transcript, protein, and methylation abundances) infer the effect of these mediators on risk for disorders. One group of post-GWAS approaches is commonly referred to as transcriptome/proteome/methylome-wide association studies, which are abbreviated as T/P/MWAS (or collectively as XWAS). Since these approaches use biological mediators, the multiple testing burden is reduced to the number of genes (∼20,000) instead of millions of GWAS SNPs, which leads to increased signal detection. In this work, our aim is to uncover likely risk genes for NPSUDs by performing XWAS analyses in two tissues-blood and brain. First, to identify putative causal risk genes, we performed an XWAS using the Summary-data-based Mendelian randomization, which uses GWAS summary statistics, reference xQTL data, and a reference LD panel. Second, given the large comorbidities among NPSUDs and the shared cis-xQTLs between blood and the brain, we improved XWAS signal detection for underpowered analyses by performing joint concordance analyses between XWAS results i) across the two tissues and ii) across NPSUDs. All XWAS signals i) were adjusted for heterogeneity in dependent instruments (HEIDI) (non-causality)
-values and ii) used to test for pathway enrichment. The results suggest that there were widely shared gene/protein signals within the major histocompatibility complex region on chromosome 6 (
and
) and elsewhere in the genome (
and
). The identification of putative molecular genes and pathways underlying risk may offer new targets for therapeutic development. Our study revealed an enrichment of XWAS signals in vitamin D and omega-3 gene sets. So, including vitamin D and omega-3 in treatment plans may have a modest but beneficial effect on patients with bipolar disorder.
Summary Background Treatments with survival benefit are greatly needed for women with heavily pretreated metastatic breast cancer. Eribulin mesilate is a non-taxane microtubule dynamics inhibitor ...with a novel mode of action. We aimed to compare overall survival of heavily pretreated patients receiving eribulin versus currently available treatments. Methods In this phase 3 open-label study, women with locally recurrent or metastatic breast cancer were randomly allocated (2:1) to eribulin mesilate (1·4 mg/m2 administered intravenously during 2–5 min on days 1 and 8 of a 21-day cycle) or treatment of physician's choice (TPC). Patients had received between two and five previous chemotherapy regimens (two or more for advanced disease), including an anthracycline and a taxane, unless contraindicated. Randomisation was stratified by geographical region, previous capecitabine treatment, and human epidermal growth factor receptor 2 status. Patients and investigators were not masked to treatment allocation. The primary endpoint was overall survival in the intention-to-treat population. This study is registered at ClinicalTrials.gov , number NCT00388726. Findings 762 women were randomly allocated to treatment groups (508 eribulin, 254 TPC). Overall survival was significantly improved in women assigned to eribulin (median 13·1 months, 95% CI 11·8–14·3) compared with TPC (10·6 months, 9·3–12·5; hazard ratio 0·81, 95% CI 0·66–0·99; p=0·041). The most common adverse events in both groups were asthenia or fatigue (270 54% of 503 patients on eribulin and 98 40% of 247 patients on TPC at all grades) and neutropenia (260 52% patients receiving eribulin and 73 30% of those on TPC at all grades). Peripheral neuropathy was the most common adverse event leading to discontinuation from eribulin, occurring in 24 (5%) of 503 patients. Interpretation Eribulin showed a significant and clinically meaningful improvement in overall survival compared with TPC in women with heavily pretreated metastatic breast cancer. This finding challenges the notion that improved overall survival is an unrealistic expectation during evaluation of new anticancer therapies in the refractory setting. Funding Eisai.
Genotype imputation across populations of mixed ancestry is critical for optimal discovery in large‐scale genome‐wide association studies (GWAS). Methods for direct imputation of GWAS ...summary‐statistics were previously shown to be practically as accurate as summary statistics produced after raw genotype imputation, while incurring orders of magnitude lower computational burden. Given that direct imputation needs a precise estimation of linkage‐disequilibrium (LD) and that most of the methods using a small reference panel for example, ~2,500‐subject coming from the 1000 Genome‐Project, there is a great need for much larger and more diverse reference panels. To accurately estimate the LD needed for an exhaustive analysis of any cosmopolitan cohort, we developed DISTMIX2. DISTMIX2: (a) uses a much larger and more diverse reference panel compared to traditional reference panels, and (b) can estimate weights of ethnic‐mixture based solely on Z‐scores, when allele frequencies are not available. We applied DISTMIX2 to GWAS summary‐statistics from the psychiatric genetic consortium (PGC). DISTMIX2 uncovered signals in numerous new regions, with most of these findings coming from the rarer variants. Rarer variants provide much sharper location for the signals compared with common variants, as the LD for rare variants extends over a lower distance than for common ones. For example, while the original PGC post‐traumatic stress disorder GWAS found only 3 marginal signals for common variants, we now uncover a very strong signal for a rare variant in PKN2, a gene associated with neuronal and hippocampal development. Thus, DISTMIX2 provides a robust and fast (re)imputation approach for most psychiatric GWAS‐studies.
Purpose
To assess the safety, pharmacokinetics, pharmacodynamics, and preliminary efficacy of cetrelimab (JNJ-63723283), a monoclonal antibody programmed cell death protein-1 (PD-1) inhibitor, in ...patients with advanced/refractory solid tumors in the phase 1/2 LUC1001 study.
Methods
In phase 1, patients with advanced solid tumors received intravenous cetrelimab 80, 240, 460, or 800 mg every 2 weeks (Q2W) or 480 mg Q4W. In phase 2, patients with melanoma, non-small-cell lung cancer (NSCLC), and microsatellite instability–high (MSI-H)/DNA mismatch repair-deficient colorectal cancer (CRC) received cetrelimab 240 mg Q2W. Response was assessed Q8W until Week 24 and Q12W thereafter.
Results
In phase 1, 58 patients received cetrelimab. Two dose-limiting toxicities were reported and two recommended phase 2 doses (RP2D) were defined (240 mg Q2W or 480 mg Q4W). After a first dose, mean maximum serum concentrations (
C
max
) ranged from 24.7 to 227.0 µg/mL; median time to
C
max
ranged from 2.0 to 3.2 h. Pharmacodynamic effect was maintained throughout the dosing period across doses. In phase 2, 146 patients received cetrelimab 240 mg Q2W. Grade ≥ 3 adverse events (AEs) occurred in 53.9% of patients. Immune-related AEs (any grade) occurred in 35.3% of patients (grade ≥ 3 in 6.9%). Overall response rate was 18.6% across tumor types, 34.3% in NSCLC, 52.6% in programmed death ligand 1–high (≥ 50% by immunohistochemistry) NSCLC, 28.0% in melanoma, and 23.8% in centrally confirmed MSI-H CRC.
Conclusions
The RP2D for cetrelimab was established. Pharmacokinetic/pharmacodynamic characteristics, safety profile, and clinical activity of cetrelimab in immune-sensitive advanced cancers were consistent with known PD-1 inhibitors.
Trial registrations
NCT02908906 at ClinicalTrials.gov, September 21, 2016; EudraCT 2016–002,017-22 at clinicaltrialsregister.eu, Jan 11, 2017.
Abstract MicroRNAs (miRNAs) are small non-coding RNAs that mainly function as negative regulators of gene expression ( Lai, 2002 ) and have been shown to be involved in schizophrenia etiology through ...genetic and expression studies ( Burmistrova et al., 2007; Hansen et al., 2007a; Perkins et al., 2007; Beveridge et al., 2010; Kim et al., 2010 ). In a mega analysis of genome-wide association study (GWAS) of schizophrenia (SZ) and bipolar disorders (BP), a polymorphism (rs1625579) located in the primary transcript of a miRNA gene, hsa-miR-137, was reported to be strongly associated with SZ. Four SZ loci ( CACNA1C , TCF4 , CSMD1 , C10orf26 ) achieving genome-wide significance in the same study were predicted and later experimentally validated (Kwon et al., 2011) as hsa-miR-137 targets. Here, using in silico , cellular and luciferase based approaches we also provide evidence that another well replicated candidate schizophrenia gene, ZNF804A , is also target for hsa-miR-137.
Rare variant-based analyses are beginning to identify risk genes for neuropsychiatric disorders and other diseases. However, the identified genes only account for a fraction of predicted causal ...genes. Recent studies have shown that rare damaging variants are significantly enriched in specific gene-sets. Methods which are able to jointly model rare variants and gene-sets to identify enriched gene-sets and use these enriched gene-sets to prioritize additional risk genes could improve understanding of the genetic architecture of diseases.
We propose DECO (Integrated analysis of de novo mutations, rare case/control variants and omics information via gene-sets), an integrated method for rare-variant and gene-set analysis. The method can (i) test the enrichment of gene-sets directly within the statistical model, and (ii) use enriched gene-sets to rank existing genes and prioritize additional risk genes for tested disorders. In simulations, DECO performs better than a homologous method that uses only variant data. To demonstrate the application of the proposed protocol, we have applied this approach to rare-variant datasets of schizophrenia. Compared with a method which only uses variant information, DECO is able to prioritize additional risk genes.
DECO can be used to analyze rare-variants and biological pathways or cell types for any disease. The package is available on Github https://github.com/hoangtn/DECO.
Alcohol consumption is known to lead to gene expression changes in the brain. After performing weighted gene co-expression network analyses (WGCNA) on genome-wide mRNA and microRNA (miRNA) expression ...in Nucleus Accumbens (NAc) of subjects with alcohol dependence (AD; N = 18) and of matched controls (N = 18), six mRNA and three miRNA modules significantly correlated with AD were identified (Bonferoni-adj. p≤ 0.05). Cell-type-specific transcriptome analyses revealed two of the mRNA modules to be enriched for neuronal specific marker genes and downregulated in AD, whereas the remaining four mRNA modules were enriched for astrocyte and microglial specific marker genes and upregulated in AD. Gene set enrichment analysis demonstrated that neuronal specific modules were enriched for genes involved in oxidative phosphorylation, mitochondrial dysfunction and MAPK signaling. Glial-specific modules were predominantly enriched for genes involved in processes related to immune functions, i.e. cytokine signaling (all adj. p≤ 0.05). In mRNA and miRNA modules, 461 and 25 candidate hub genes were identified, respectively. In contrast to the expected biological functions of miRNAs, correlation analyses between mRNA and miRNA hub genes revealed a higher number of positive than negative correlations (χ2 test p≤ 0.0001). Integration of hub gene expression with genome-wide genotypic data resulted in 591 mRNA cis-eQTLs and 62 miRNA cis-eQTLs. mRNA cis-eQTLs were significantly enriched for AD diagnosis and AD symptom counts (adj. p = 0.014 and p = 0.024, respectively) in AD GWAS signals in a large, independent genetic sample from the Collaborative Study on Genetics of Alcohol (COGA). In conclusion, our study identified putative gene network hubs coordinating mRNA and miRNA co-expression changes in the NAc of AD subjects, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.