Functional dissection of complex and molecular trait variants at single nucleotide resolution

Identifying the causal variants and mechanisms that drive complex traits and diseases remains a core problem in human genetics. The majority of these variants have individually weak effects and lie in non-coding gene-regulatory elements where we lack a complete understanding of how single nucleotide alterations modulate transcriptional processes to affect human phenotypes. To address this, we measured the activity of 221,412 trait-associated variants that had been statistically fine-mapped using a Massively Parallel Reporter Assay (MPRA) in 5 diverse cell-types. We show that MPRA is able to discriminate between likely causal variants and controls, identifying 12,025 regulatory variants with high precision. Although the effects of these variants largely agree with orthogonal measures of function, only 69% can plausibly be explained by the disruption of a known transcription factor (TF) binding motif. We dissect the mechanisms of 136 variants using saturation mutagenesis and assign impacted TFs for 91% of variants without a clear canonical mechanism. Finally, we provide evidence that epistasis is prevalent for variants in close proximity and identify multiple functional variants on the same haplotype at a small, but important, subset of trait-associated loci. Overall, our study provides a systematic functional characterization of likely causal common variants underlying complex and molecular human traits, enabling new insights into the regulatory grammar underlying disease risk.

Myocardial Recovery in Recent Onset Dilated Cardiomyopathy: Role of CDCP1 and Cardiac Fibrosis

BACKGROUND

Dilated cardiomyopathy (DCM) is a major cause of heart failure and carries a high mortality rate. Myocardial recovery in DCM-related heart failure patients is highly variable, with some patients having little or no response to standard drug therapy. A genome-wide association study may agnostically identify biomarkers and provide novel insight into the biology of myocardial recovery in DCM.

METHODS

A genome-wide association study for change in left ventricular ejection fraction was performed in 686 White subjects with recent-onset DCM who received standard pharmacotherapy. Genome-wide association study signals were subsequently functionally validated and studied in relevant cellular models to understand molecular mechanisms that may have contributed to the change in left ventricular ejection fraction.

RESULTS

The genome-wide association study identified a highly suggestive locus that mapped to the 5'-flanking region of the CDCP1 (CUB [complement C1r/C1s, Uegf, and Bmp1] domain containing protein 1) gene (rs6773435; P=7.12×10-7). The variant allele was associated with improved cardiac function and decreased CDCP1 transcription. CDCP1 expression was significantly upregulated in human cardiac fibroblasts (HCFs) in response to the PDGF (platelet-derived growth factor) signaling, and knockdown of CDCP1 significantly repressed HCF proliferation and decreased AKT (protein kinase B) phosphorylation. Transcriptomic profiling after CDCP1 knockdown in HCFs supported the conclusion that CDCP1 regulates HCF proliferation and mitosis. In addition, CDCP1 knockdown in HCFs resulted in significantly decreased expression of soluble ST2 (suppression of tumorigenicity-2), a prognostic biomarker for heart failure and inductor of cardiac fibrosis.

CONCLUSIONS

CDCP1 may play an important role in myocardial recovery in recent-onset DCM and mediates its effect primarily by attenuating cardiac fibrosis.

Pharmacological Targeting of Androgen Receptor Elicits Context-Specific Effects in Estrogen Receptor-Positive Breast Cancer

Androgen receptor (AR) is expressed in 80% to 90% of estrogen receptor α-positive (ER+) breast cancers. Accumulated evidence has shown that AR is a tumor suppressor and that its expression is associated with improved prognosis in ER+ breast cancer. However, both a selective AR agonist (RAD140) and an AR inhibitor (enzalutamide, ENZ) have shown a therapeutic effect on ER+ breast cancer, so the potential for clinical application of AR-targeting therapy for ER+ breast cancer is still in dispute. In this study, we evaluated the efficacy of ENZ and RAD140 in vivo and in vitro in AR+/ER+ breast cancer models, characterizing the relationship of AR and ER levels to response to AR-targeting drugs and investigating the alterations of global gene expression and chromatin binding of AR and ERα after ENZ treatment. In the AR-low setting, ENZ directly functioned as an ERα antagonist. Cell growth inhibition by ENZ in breast cancer with low AR expression was independent of AR and instead dependent on ER. In AR-high breast cancer models, AR repressed ERα signaling and ENZ promoted ERα signaling by antagonizing AR. In contrast, RAD140 activated AR signaling and suppressed AR-high tumor growth by deregulating ERα expression and blocking ERα function. Overall, analysis of the dynamic efficacies and outcomes of AR agonist, and antagonist in the presence of different AR and ERα levels reveals regulators of response and supports the clinical investigation of ENZ in selected ER+ tumors with a low AR/ER ratio and AR agonists in tumors with a high AR/ER ratio.

SIGNIFICANCE

The ratio of androgen receptor to estrogen receptor in breast cancer dictates the response to AR-targeted therapies, providing guidelines for developing AR-directed treatment strategies for patients with breast cancer.

Glucocorticoids mediate transcriptome-wide alternative polyadenylation: Potential mechanistic and clinical implications

Alternative polyadenylation (APA) is a common genetic regulatory mechanism that generates distinct 3' ends for RNA transcripts. Changes in APA have been associated with multiple biological processes and disease phenotypes. However, the role of hormones and their drug analogs in APA remains largely unknown. In this study, we investigated transcriptome-wide the impact of glucocorticoids on APA in 30 human B-lymphoblastoid cell lines. We found that glucocorticoids could regulate APA for a subset of genes, possibly by changing the expression of 142 RNA-binding proteins, some with known APA-regulating properties. Interestingly, genes with glucocorticoid-mediated APA were enriched in viral translation-related pathways, while genes with glucocorticoid-mediated expression were enriched in interferon and interleukin pathways, suggesting that glucocorticoid-mediated APA might result in functional consequences distinct from gene expression. For example, glucocorticoids, a pharmacotherapy for severe COVID-19, were found to change the APA but not the expression of LY6E, an important antiviral inhibitor in coronavirus diseases. Glucocorticoid-mediated APA was also cell-type-specific, suggesting an action of glucocorticoids that may be unique to immune regulation. We also observed evidence for genotype-dependent glucocorticoid-mediated APA (referred to as pharmacogenomic-alterative polyadenylation quantitative trait loci), providing potential functional mechanisms for a series of common genetic variants that had previously been associated with immune disorders, but without a clear mechanism. In summary, this study reports a series of observations regarding the impact of glucocorticoids on APA, raising the possibility that this mechanism might have implications for both disease pathophysiology and drug therapy.

Mood Regulatory Actions of Active and Sham Nucleus Accumbens Deep Brain Stimulation in Antidepressant Resistant Rats

The antidepressant actions of deep brain stimulation (DBS) are associated with progressive neuroadaptations within the mood network, modulated in part, by neurotrophic mechanisms. We investigated the antidepressant-like effects of chronic nucleus accumbens (NAc) DBS and its association with change in glycogen synthase kinase 3 (GSK3) and mammalian target of rapamycin (mTOR) expression in the infralimbic cortex (IL), and the dorsal (dHIP) and ventral (vHIP) subregions of the hippocampus of antidepressant resistant rats. Antidepressant resistance was induced via daily injection of adrenocorticotropic hormone (ACTH; 100 μg/day; 15 days) and confirmed by non-response to tricyclic antidepressant treatment (imipramine, 10 mg/kg). Portable microdevices provided continuous bilateral NAc DBS (130 Hz, 200 μA, 90 μs) for 7 days. A control sham electrode group was included, together with ACTH- and saline-treated control groups. Home cage monitoring, open field, sucrose preference, and, forced swim behavioral tests were performed. Post-mortem levels of GSK3 and mTOR, total and phosphorylated, were determined with Western blot. As previously reported, ACTH treatment blocked the immobility-reducing effects of imipramine in the forced swim test. In contrast, treatment with either active DBS or sham electrode placement in the NAc significantly reduced forced swim immobility time in ACTH-treated animals. This was associated with increased homecage activity in the DBS and sham groups relative to ACTH and saline groups, however, no differences in locomotor activity were observed in the open field test, nor were any group differences seen for sucrose consumption across groups. The antidepressant-like actions of NAc DBS and sham electrode placements were associated with an increase in levels of IL and vHIP phospho-GSK3β and phospho-mTOR, however, no differences in these protein levels were observed in the dHIP region. These data suggest that early response to electrode placement in the NAc, irrespective of whether active DBS or sham, has antidepressant-like effects in the ACTH-model of antidepressant resistance associated with distal upregulation of phospho-GSK3β and phospho-mTOR in the IL and vHIP regions of the mood network.

Selective Serotonin Reuptake Inhibitor Pharmaco-Omics: Mechanisms and Prediction

Selective serotonin reuptake inhibitors (SSRIs) are a standard of care for the pharmacotherapy of patients suffering from Major Depressive Disorder (MDD). However, only one-half to two-thirds of MDD patients respond to SSRI therapy. Recently, a "multiple omics" research strategy was applied to identify genetic differences between patients who did and did not respond to SSRI therapy. As a first step, plasma metabolites were assayed using samples from the 803 patients in the PGRN-AMPS SSRI MDD trial. The metabolomics data were then used to "inform" genomics by performing a genome-wide association study (GWAS) for plasma concentrations of the metabolite most highly associated with clinical response, serotonin (5-HT). Two genome-wide or near genome-wide significant single nucleotide polymorphism (SNP) signals were identified, one that mapped near the TSPAN5 gene and another across the ERICH3 gene, both genes that are highly expressed in the brain. Knocking down TSPAN5 and ERICH3 resulted in decreased 5-HT concentrations in neuroblastoma cell culture media and decreased expression of enzymes involved in 5-HT biosynthesis and metabolism. Functional genomic studies demonstrated that ERICH3 was involved in clathrin-mediated vesicle formation and TSPAN5 was an ethanol-responsive gene that may be a marker for response to acamprosate pharmacotherapy of alcohol use disorder (AUD), a neuropsychiatric disorder highly co-morbid with MDD. In parallel studies, kynurenine was the plasma metabolite most highly associated with MDD symptom severity and application of a metabolomics-informed pharmacogenomics approach identified DEFB1 and AHR as genes associated with variation in plasma kynurenine levels. Both genes also contributed to kynurenine-related inflammatory pathways. Finally, a multiply replicated predictive algorithm for SSRI clinical response with a balanced predictive accuracy of 76% (compared with 56% for clinical data alone) was developed by including the SNPs in TSPAN5, ERICH3, DEFB1 and AHR. In summary, application of a multiple omics research strategy that used metabolomics to inform genomics, followed by functional genomic studies, identified novel genes that influenced monoamine biology and made it possible to develop a predictive algorithm for SSRI clinical outcomes in MDD. A similar pharmaco-omic research strategy might be broadly applicable for the study of other neuropsychiatric diseases and their drug therapy.

Mood-Stabilizing Antiepileptic Treatment Response in Bipolar Disorder: A Genome-Wide Association Study

Several antiepileptic drugs (AEDs) have US Food and Drug Administration (FDA) approval for use as mood stabilizers in bipolar disorder (BD), but not all BD patients respond to these AED mood stabilizers (AED‐MSs). To identify genetic polymorphisms that contribute to the variability in AED‐MS response, we performed a discovery genome‐wide association study (GWAS) of 199 BD patients from the Mayo Clinic Bipolar Disorder Biobank. Most of these patients had been treated with the AED‐MS valproate/divalproex and/or lamotrigine. AED‐MS response was assessed using the Alda scale, which quantifies clinical improvement while accounting for potential confounding factors. We identified two genome‐wide significant single‐nucleotide polymorphism (SNP) signals that mapped to the THSD7A (rs78835388, P = 7.1E‐09) and SLC35F3 (rs114872993, P = 3.2E‐08) genes. We also identified two genes with statistically significant gene‐level associations: ABCC1 (P = 6.7E‐07; top SNP rs875740, P = 2.0E‐6), and DISP1 (P = 8.9E‐07; top SNP rs34701716, P = 8.9E‐07). THSD7A SNPs were previously found to be associated with risk for several psychiatric disorders, including BD. Both THSD7A and SLC35F3 are expressed in excitatory/glutamatergic and inhibitory/γ‐aminobutyric acidergic (GABAergic) neurons, which are targets of AED‐MSs. ABCC1 is involved in the transport of valproate and lamotrigine metabolites, and the SNPs in ABCC1 and DISP1 with the strongest evidence of association in our GWAS are strong splicing quantitative trait loci in the human gut, suggesting a possible influence on drug absorption. In conclusion, our pharmacogenomic study identified novel genetic loci that appear to contribute to AED‐MS treatment response, and may facilitate precision medicine in BD.

A review of brain insulin signaling in mood disorders: From biomarker to clinical target

Patients with mood disorders are at increased risk for metabolic dysfunction. Co-occurrence of the two conditions is typically associated with a more severe disease course and poorer treatment outcomes. The specific pathophysiological mechanisms underlying this bidirectional relationship between mood and metabolic dysfunction remains poorly understood. However, it is likely that impairment of metabolic processes within the brain play a critical role. The insulin signaling pathway mediates metabolic homeostasis and is important in the regulation of neurotrophic and synaptic plasticity processes, including those involved in neurodegenerative diseases like Alzheimer's. Thus, insulin signaling in the brain may serve to link metabolic function and mood. Central insulin signaling is mediated through locally secreted insulin and widespread insulin receptor expression. Here we review the preclinical and clinical data addressing the relationships between central insulin signaling, cellular metabolism, neurotrophic processes, and mood regulation, including key points of mechanistic overlap. These relationships have important implications for developing biomarker-based diagnostics and precision medicine approaches to treat severe mood disorders.

Corifollitropin alfa vs recombinant FSH for controlled ovarian stimulation in women aged 35-42 years with a body weight ≥50 kg: a randomized controlled trial

STUDY QUESTION

Is corifollitropin alfa 150 μg equivalent to follitropin beta 300 IU/day for controlled ovarian hyperstimulation (COS) in older women weighing ≥50 kg undergoing IVF and/or ICSI in Vietnam?

SUMMARY ANSWER

Corifollitropin alfa 150 μg was equivalent to follitropin beta 300 IU/day with respect to the number of oocytes retrieved, the ongoing, cumulative and live birth rates and obstetric outcomes.

WHAT IS KNOWN ALREADY

Corifollitropin alfa is a recombinant FSH (rFSH) preparation with slow absorption and a long half-life allowing administration of a single dose for COS lasting 7 days. Several randomized, controlled clinical trials have reported that COS with corifollitropin alfa is associated with similar outcomes compared with COS using daily rFSH. However, limited data are available in Asian patients.

STUDY DESIGN, SIZE, DURATION

This randomized controlled trial was conducted at a single large IVF centre in Vietnam from June 2015 to August 2016. A total of 400 patients were included, 200 in each treatment group. The primary outcome measure was the number of oocytes retrieved. Patients were followed for 1 year after randomization.

PARTICIPANTS /MATERIALS, SETTING, METHODS

Participants aged 35–42 years with a body weight ≥50 kg who were undergoing an IVF cycle were randomized to undergo COS with a single dose of corifollitropin alfa 150 μg on Day 2 or 3 of the menstrual cycle, or follitropin beta 300 IU/day for 7 days starting on Day 2 or 3 of the menstrual cycle. All underwent ICSI according to standard institutional protocols. A beta hCG test was performed 17 days after ovum pick-up, and positive tests were confirmed on vaginal and/or abdominal ultrasound at 5–6 weeks after embryo transfer (clinical pregnancy) and at ≥10 weeks (ongoing pregnancy). Rates of ovarian hyperstimulation syndrome, and maternal and foetal outcomes after one cycle of ICSI were monitored over 12 months.

MAIN RESULTS AND THE ROLE OF CHANCE

Patients in the corifollitropin alfa and follitropin beta groups were well matched at baseline (mean age 37.5 ± 1.9 vs 37.7 ± 2.0 years, mean body weight 53.7 ± 5.4 vs 52.5 ± 4.8 kg). There was no significant difference between the corifollitropin alfa and follitropin beta groups in the number of oocytes retrieved (11.4 ± 5.9 vs 10.8 ± 5.8; P = 0.338). The ongoing pregnancy rate (31.5 vs 32.0%; P = 0.99) and live birth rate (30.5 vs 32.0%; P = 0.83) (both per initiated cycle at 12 months after randomization) were also similar in the two treatment groups. Complication rates were low and similar in the corifollitropin alfa and follitropin beta groups, and there were no significant between-group differences in obstetric outcomes.

LIMITATIONS, REASONS FOR CAUTION

This study had an open-label design, and therefore, the potential for bias cannot be excluded. The findings are only applicable to patient populations with similar characteristics to those enroled in the study.

WIDER IMPLICATIONS OF THE FINDINGS

This study adds to the body of evidence supporting the equivalence of corifollitropin alfa and follitropin beta for COS in a variety of patients undergoing IVF and/or ICSI. The ability to provide COS with corifollitropin alfa has the potential to reduce the burden of treatment for patients.

STUDY FUNDING/COMPETING INTERESTS

This study was supported by Merck Sharp and Dohme. The authors state that they have no financial or commercial conflicts of interest.

TRIAL REGISTRATION NUMBER

The trial was registered with clinicaltrials.gov (NCT02466204).

TRIAL REGISTRATION DATE

2 June 2015.

DATE OF FIRST PATIENT’S ENROLMENT

19 June 2015.