Allele-specific Binding Research Articles

Genetic variation in histone modifications and gene expression identifies regulatory variants in the mammary gland of cattle.

Causal variants for complex traits, such as eQTL are often found in non-coding regions of the genome, where they are hypothesised to influence phenotypes by regulating gene expression. Many regulatory regions are marked by histone modifications, which can be assayed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). Sequence reads from ChIP-seq form peaks at putative regulatory regions, which may reflect the amount of regulatory activity at this region. Therefore, eQTL which are also associated with differences in histone modifications are excellent candidate causal variants. We assayed the histone modifications H3K4Me3, H3K4Me1 and H3K27ac and mRNA in the mammary gland of up to 400 animals. We identified QTL for peak height (histone QTL), exon expression (eeQTL), allele specific expression (aseQTL) and allele specific binding (asbQTL). By intersecting these results, we identify variants which may influence gene expression by altering regulatory regions of the genome, and may be causal variants for other traits. Lastly, we find that these variants are found in putative transcription factor binding sites, identifying a mechanism for the effect of many eQTL. We find that allele specific and traditional QTL analysis often identify the same genetic variants and provide evidence that many eQTL are regulatory variants which alter activity at regulatory regions of the bovine genome. Our work provides methodological and biological updates on how regulatory mechanisms interplay at multi-omics levels.

Characterization of sequence determinants of enhancer function using natural genetic variation.

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F1-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.

Abstract P1108: Functional Genomics Analysis Of Sox17 Locus To Define The Pathogenic Role Of FUBP1 In Pulmonary Arterial Hypertension

Rationale: Pulmonary arterial hypertension (PAH) is an enigmatic and morbid disease where insights are emerging regarding genetic susceptibility. Genome-wide Association Studies (GWAS) have identified SOX17 as the most significant PAH-associated genomic locus. The allele of tag SNP in this locus is associated with 1.8-fold higher PAH risk. It has been challenging to define the mechanisms underlying the contribution of the PAH-associated functional SNPs (fSNPs) to pathogenesis of the disease. Methods: We developed a post-GWAS functional genomics strategy to define the causative fSNPs and identify the associated biological mechanisms. This analysis includes Reel-seq (Regulatory element-sequencing), an EMSA-based high-throughput technique to identify fSNPs in a synthetic DNA library containing PAH-associated SNPs; SDCP-MS (SNP-specific DNA competition pulldown-mass spectrometry) to identify proteins that specifically bind to fSNPs; and AIDP-Wb (allele-imbalanced DNA pulldown-Western blot) to show allele-imbalanced binding of these proteins to fSNPs. The regulation of risk gene expression by these fSNP-binding proteins and their pathogenicity were determined in human pulmonary arterial endothelial cells (PAECs) and confirmed in PAH animal models and patients. Results: By using high-throughput Reel-seq and subsequent validation with EMSA, intergenic SNP rs4738801 in SOX17 locus was identified as a fSNP from a library containing 254 PAH-associated haplotype SNPs. This fSNP resides in a remote upstream enhancer region of SOX17, an endothelial effector increasingly associated with PAH pathogenesis. Using SDCP-MS and AIDP-Wb, we found that the transcription factor FUBP1 binds to rs4738801 risk allele C with lower affinity than non-risk allele G, resulting in a decrease in SOX17 expression. FUBP1 and target gene SOX17 controlled PAH-associated pathophenotypes in PAECs, including proliferation, apoptosis, and angiogenesis. Downregulated by the major acquired PAH trigger hypoxia, FUBP1 and SOX17 were decreased in lungs and pulmonary ECs isolated from PAH patients and mouse models. A 3.77-fold enrichment of fSNP rs4738801 risk allele C was found in patients with PAH induced by hypoxia, but not in PAH associated with connective tissue disease or congenital heart disease. Conclusions: FUBP1 controls SOX17 expression via allele-specific binding to PAH-associated fSNP rs4738801. The reduced binding of FUBP1 to risk allele C defines the genomic architecture contributing to the SOX17-dependent genetic susceptibility of PAH. The downregulation of FUBP1-SOX17 by hypoxia results in endothelial dysfunction, contributing to the acquired pathogenesis of PAH. These findings identify a novel role of FUBP1 in the functional regulation of SOX17 locus and elucidating a pathogenic mechanism that combines the acquired PAH-triggering factors and altered genetic susceptibility.

NFIA determines the cis-effect of genetic variation on Ucp1 expression in murine thermogenic adipocytes

SummaryThermogenic brown and beige adipocytes counteract obesity by enhancing energy dissipation via uncoupling protein-1 (Ucp1). However, the effect of genetic variation on these cells, a major source of disease susceptibility, has been less well studied. Here we examined beige adipocytes from obesity-prone C57BL/6J (B6) and obesity-resistant 129X1/SvJ (129) mouse strains and identified a cis-regulatory variant rs47238345 that is responsible for differential Ucp1 expression. The alternative T allele of rs47238345 at the Ucp1 -12kb enhancer in 129 facilitates the allele-specific binding of nuclear factor I-A (NFIA) to mediate allele-specific enhancer-promoter interaction and Ucp1 transcription. Furthermore, CRISPR-Cas9/Cpf1-mediated single nucleotide polymorphism (SNP) editing of rs47238345 resulted in increased Ucp1 expression. We also identified Lim homeobox protein 8 (Lhx8), whose expression is higher in 129 than in B6, as a trans-acting regulator of Ucp1 in mice and humans. These results demonstrate the cis- and trans-acting effects of genetic variation on Ucp1 expression that underlie phenotypic diversity.

Abstract 1475: Uncovering the functional variants and target genes of the 7q32 pancreatic cancer risk locus

Abstract In western nations, pancreatic ductal adenocarcinoma (PDAC) is set to become second leading cause of cancer mortality over the next decade. Despite advances in treatment, PDAC has a 5-year survival rate of only ~9%, necessitating a better understanding of its etiology. Risk factors for PDAC are both environmental and genetic, with heritability estimated at around 20%. However, highly penetrant coding variants are rare within the population and have been estimated to account for as little as 1% of the phenotypic variance for PDAC. In an effort to identify common PDAC susceptibility alleles, genome-wide association studies (GWAS) have identified 17 loci at which germline variants associate with PDAC risk. Mostly localizing to non-coding loci, many such risk variants are expected to confer risk via allele-specific changes in transcription factor binding, leading to altered target gene expression over a lifetime. The identification of functional variants from GWAS is complicated by the correlation of alleles through linkage disequilibrium (LD). Here, we have leveraged Bayesian fine mapping to identify a credible set of seven single nucleotide polymorphisms (SNPs) in high LD at the 7q32 risk locus. Fine-mapped variants were annotated using epigenomic datasets we previously generated from a range of pancreatic cell lines, organoids and primary patient samples, encompassing histone modification marks (ChIP-seq), open chromatin (ATAC-seq) and chromatin conformation capture (Capture-HiC, HiChIP). Two variants from the credible set (rs6971499 and rs6970779) were found to consistently lie within regions annotated as super-enhancers and marked by open chromatin within our data. Both SNPs exhibit chromatin looping to the promoters of nearby plausible target genes LINC-PINT and miR29a. Electrophoretic mobility shift assays (EMSA) illustrate that rs6971499-T and rs6970779-G exhibit allele-specific protein binding in PANC-1 and MIA PaCa-2 cell lines, while luciferase reporter assays in the same cell lines demonstrate that both alleles drive significant changes in transcriptional activity. Ongoing proteomic experiments aim to identify the transcription factors differentially bound to these alleles, while CRISPR methods are being employed to validate enhancer variants in their native chromatin context. Future experiments may further characterise these transcription factors and their target genes with an eye toward better understanding their roles in PDAC carcinogenesis. Citation Format: Aidan O'Brien, Jason W. Hoskins, Daina R. Eiser, Katelyn E. Connelly, Jun Zhong, Thorkell Andresson, Irene Collins, Pancreatic Cancer Cohort Consortium, Pancreatic Cancer Case-Control Consortium, Stephen J. Chanock, Alison P. Klein, H. Efsun Arda, Laufey T. Amundadottir. Uncovering the functional variants and target genes of the 7q32 pancreatic cancer risk locus [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1475.

Discovery and Functional Characterization of Two Regulatory Variants Underlying Lupus Susceptibility at 2p13.1.

Genome-wide association studies have identified 2p13.1 as a prominent susceptibility locus for systemic lupus erythematosus (SLE)—a complex, multisystem autoimmune disease. However, the identity of underlying causal variant (s) and molecular mechanisms for increasing disease susceptibility are poorly understood. Using meta-analysis (cases = 10,252, controls = 21,604) followed by conditional analysis, bioinformatic annotation, and eQTL and 3D-chromatin interaction analyses, we computationally prioritized potential functional variants and subsequently experimentally validated their effects. Ethnicity-specific meta-analysis revealed striking allele frequency differences between Asian and European ancestries, but with similar odds ratios. We identified 20 genome-wide significant (p < 5 × 10−8) variants, and conditional analysis pinpointed two potential functional variants, rs6705628 and rs2272165, likely to explain the association. The two SNPs are near DGUOK, mitochondrial deoxyguanosine kinase, and its associated antisense RNA DGUOK-AS1. Using luciferase reporter gene assays, we found significant cell type- and allele-specific promoter activity at rs6705628 and enhancer activity at rs2272165. This is supported by ChIP-qPCR showing allele-specific binding with three histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), transcriptional coactivator p300, CCCTC-binding factor (CTCF), and transcription factor ARID3A. Transcriptome data across 28 immune cell types from Asians showed both SNPs are cell-type-specific but only in B-cells. Splicing QTLs showed strong regulation of DGUOK-AS1. Genotype-specific DGOUK protein levels are supported by Western blots. Promoter capture Hi-C data revealed long-range chromatin interactions between rs2272165 and several nearby promoters, including DGUOK. Taken together, we provide mechanistic insights into how two noncoding variants underlie SLE risk at the 2p13.1 locus.

Evolution of transcription factor binding through sequence variations and turnover of binding sites.

Variations in noncoding regulatory sequences play a central role in evolution. Interpreting such variations, however, remains difficult even in the context of defined attributes such as transcription factor (TF) binding sites. Here, we systematically link variations in cis-regulatory sequences to TF binding by profiling the allele-specific binding of 27 TFs expressed in a yeast hybrid, in which two related genomes are present within the same nucleus. TFs localize preferentially to sites containing their known consensus motifs but occupy only a small fraction of the motif-containing sites available within the genomes. Differential binding of TFs to the orthologous alleles was well explained by variations that alter motif sequence, whereas differences in chromatin accessibility between alleles were of little apparent effect. Motif variations that abolished binding when present in only one allele were still bound when present in both alleles, suggesting evolutionary compensation, with a potential role for sequence conservation at the motif's vicinity. At the level of the full promoter, we identify cases of binding-site turnover, in which binding sites are reciprocally gained and lost, yet most interspecific differences remained uncompensated. Our results show the flexibility of TFs to bind imprecise motifs and the fast evolution of TF binding sites between related species.

Multiplex Fragment Analysis for Flexible Detection of All SARS-CoV-2 Variants of Concern.

BackgroundSevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to emerge, and effective tracking requires rapid return of results. Surveillance of variants is typically performed by whole genome sequencing (WGS), which can be financially prohibitive and requires specialized equipment and bioinformatic expertise. Genotyping approaches are rapid methods for monitoring SARS-CoV-2 variants but require continuous adaptation. Fragment analysis may represent an approach for improved SARS-CoV-2 variant detection.MethodsA multiplex fragment analysis approach (CoVarScan) was validated using PCR targeting variants by size and fluorescent color. Eight SARS-CoV-2 mutational hot spots in variants of concern (VOCs) were targeted. Three primer pairs (recurrently deleted region [RDR] 1, RDR2, and RDR3–4) flank RDRs in the S-gene. Three allele-specific primers target recurrent spike receptor binding domain mutants. Lastly, 2 primer pairs target recurrent deletions or insertions in ORF1A and ORF8. Fragments were resolved and analyzed by capillary electrophoresis (ABI 3730XL), and mutational signatures were compared to WGS results.ResultsWe validated CoVarScan using 3544 clinical respiratory specimens. The assay exhibited 96% sensitivity and 99% specificity compared to WGS. The limit of detection for the core targets (RDR1, RDR2, and ORF1A) was 5 copies/reaction. Variants were identified in 95% of samples with cycle threshold (CT) <30 and 75% of samples with a CT 34 to 35. Assay design was frozen April 2021, but all subsequent VOCs have been detected including Delta (n = 2820), Mu, (n = 6), Lambda (n = 6), and Omicron (n = 309). Genotyping results are available in as little as 4 h.ConclusionsMultiplex fragment analysis is adaptable and rapid and has similar accuracy to WGS to classify SARS-CoV-2 variants.

ANANASTRA: annotation and enrichment analysis of allele-specific transcription factor binding at SNPs.

We present ANANASTRA, https://ananastra.autosome.org, a web server for the identification and annotation of regulatory single-nucleotide polymorphisms (SNPs) with allele-specific binding events. ANANASTRA accepts a list of dbSNP IDs or a VCF file and reports allele-specific binding (ASB) sites of particular transcription factors or in specific cell types, highlighting those with ASBs significantly enriched at SNPs in the query list. ANANASTRA is built on top of a systematic analysis of allelic imbalance in ChIP-Seq experiments and performs the ASB enrichment test against background sets of SNPs found in the same source experiments as ASB sites but not displaying significant allelic imbalance. We illustrate ANANASTRA usage with selected case studies and expect that ANANASTRA will help to conduct the follow-up of GWAS in terms of establishing functional hypotheses and designing experimental verification.

A general framework for predicting the transcriptomic consequences of non-coding variation and small molecules.

Genome wide association studies (GWASs) for complex traits have implicated thousands of genetic loci. Most GWAS-nominated variants lie in noncoding regions, complicating the systematic translation of these findings into functional understanding. Here, we leverage convolutional neural networks to assist in this challenge. Our computational framework, peaBrain, models the transcriptional machinery of a tissue as a two-stage process: first, predicting the mean tissue specific abundance of all genes and second, incorporating the transcriptomic consequences of genotype variation to predict individual abundance on a subject-by-subject basis. We demonstrate that peaBrain accounts for the majority (>50%) of variance observed in mean transcript abundance across most tissues and outperforms regularized linear models in predicting the consequences of individual genotype variation. We highlight the validity of the peaBrain model by calculating non-coding impact scores that correlate with nucleotide evolutionary constraint that are also predictive of disease-associated variation and allele-specific transcription factor binding. We further show how these tissue-specific peaBrain scores can be leveraged to pinpoint functional tissues underlying complex traits, outperforming methods that depend on colocalization of eQTL and GWAS signals. We subsequently: (a) derive continuous dense embeddings of genes for downstream applications; (b) highlight the utility of the model in predicting transcriptomic impact of small molecules and shRNA (on par with in vitro experimental replication of external test sets); (c) explore how peaBrain can be used to model difficult-to-study processes (such as neural induction); and (d) identify putatively functional eQTLs that are missed by high-throughput experimental approaches.

Regulation of the Late Onset alzheimer's Disease Associated HLA-DQA1/DRB1 Expression.

(Genome-wide Association Studies) GWAS have identified ∼42 late-onset Alzheimer's disease (LOAD)-associated loci, each of which contains multiple single nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) and most of these SNPs are in the non-coding region of human genome. However, how these SNPs regulate risk gene expression remains unknown. In this work, by using a set of novel techniques, we identified 6 functional SNPs (fSNPs) rs9271198, rs9271200, rs9281945, rs9271243, and rs9271247 on the LOAD-associated HLA-DRB1/DQA1 locus and 42 proteins specifically binding to five of these 6 fSNPs. As a proof of evidence, we verified the allele-specific binding of GATA2 and GATA3, ELAVL1 and HNRNPA0, ILF2 and ILF3, NFIB and NFIC, as well as CUX1 to these five fSNPs, respectively. Moreover, we demonstrate that all these nine proteins regulate the expression of both HLA-DQA1 and HLA-DRB1 in human microglial cells. The contribution of HLA class II to the susceptibility of LOAD is discussed.

Unsupervised Mining of HLA-I Peptidomes Reveals New Binding Motifs and Potential False Positives in the Community Database.

Modern vaccine designs and studies of human leukocyte antigen (HLA)-mediated immune responses rely heavily on the knowledge of HLA allele-specific binding motifs and computational prediction of HLA-peptide binding affinity. Breakthroughs in HLA peptidomics have considerably expanded the databases of natural HLA ligands and enabled detailed characterizations of HLA-peptide binding specificity. However, cautions must be made when analyzing HLA peptidomics data because identified peptides may be contaminants in mass spectrometry or may weakly bind to the HLA molecules. Here, a hybrid de novo peptide sequencing approach was applied to large-scale mono-allelic HLA peptidomics datasets to uncover new ligands and refine current knowledge of HLA binding motifs. Up to 12-40% of the peptidomics data were low-binding affinity peptides with an arginine or a lysine at the C-terminus and likely to be tryptic peptide contaminants. Thousands of these peptides have been reported in a community database as legitimate ligands and might be erroneously used for training prediction models. Furthermore, unsupervised clustering of identified ligands revealed additional binding motifs for several HLA class I alleles and effectively isolated outliers that were experimentally confirmed to be false positives. Overall, our findings expanded the knowledge of HLA binding specificity and advocated for more rigorous interpretation of HLA peptidomics data that will ensure the high validity of community HLA ligandome databases.

The renal cancer risk allele at 14q24.2 activates a novel hypoxia-inducible transcription factor-binding enhancer of DPF3 expression

Evolution of clear cell renal cell carcinoma is guided by dysregulation of hypoxia-inducible transcription factor (HIF) pathways following loss of the von Hippel-Lindau tumor suppressor protein. Renal cell carcinoma (RCC)-associated polymorphisms influence HIF–DNA interactions at enhancers of important oncogenes thereby modulating the risk of developing renal cancer. A strong signal of genome-wide association with RCC was determined for the single nucleotide polymorphism (SNP) rs4903064, located on chr14q.24.2 within an intron of DPF3, encoding for Double PHD Fingers 3, a member of chromatin remodeling complexes; however, it is unclear how the risk allele operates in renal cells. In this study, we used tissue specimens and primary renal cells from a large cohort of RCC patients to examine the function of this polymorphism. In clear cell renal cell carcinoma tissue, isolated tumor cells as well as in primary renal tubular cells, in which HIF was stabilized, we determined genotype-specific increases of DPF3 mRNA levels and identified that the risk SNP resides in an active enhancer region, creating a novel HIF-binding motif. We then confirmed allele-specific HIF binding to this locus using chromatin immunoprecipitation of HIF subunits. Consequentially, HIF-mediated DPF3 regulation was dependent on the presence of the risk allele. Finally, we show that DPF3 deletion in proximal tubular cells retarded cell growth, indicating potential roles for DPF3 in cell proliferation. Our analyses suggest that the HIF pathway differentially operates on a SNP-induced hypoxia-response element at 14q24.2, thereby affecting DPF3 expression, which increases the risk of developing renal cancer.

Allele-specific transcription factor binding in a cellular model of orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a frequent congenital malformation with multifactorial etiology. While recent genome-wide association studies (GWAS) have identified several nsCL/P risk loci, the functional effects of the associated non-coding variants are largely unknown. Furthermore, additional risk loci remain undetected due to lack of power. As genetic variants might alter binding of transcription factors (TF), we here hypothesized that the integration of data from TF binding sites, expression analyses and nsCL/P GWAS might help to (i) identify functionally relevant variants at GWAS loci, and (ii) highlight novel risk variants that have been previously undetected. Analysing the craniofacial TF TFAP2A in human embryonic palatal mesenchyme (HEPM) cells, we identified 2845 TFAP2A ChIP-seq peaks, several of which were located near nsCL/P candidate genes (e.g. MSX1 and SPRY2). Comparison with independent data suggest that 802 of them might be specific to craniofacial development, and genes near these peaks are enriched in processes relevant to nsCL/P. Integration with nsCL/P GWAS data, however, did not show robust evidence for co-localization of common nsCL/P risk variants with TFAP2A ChIP-seq peaks. This data set represents a new resource for the analyses of craniofacial processes, and similar approaches with additional cell lines and TFs could be applied to generate further insights into nsCL/P etiology.

Comprehensive Statistical and Bioinformatics Analysis in the Deciphering of Putative Mechanisms by Which Lipid-Associated GWAS Loci Contribute to Coronary Artery Disease.

The study was designed to evaluate putative mechanisms by which lipid-associated loci identified by genome-wide association studies (GWAS) are involved in the molecular pathogenesis of coronary artery disease (CAD) using a comprehensive statistical and bioinformatics analysis. A total of 1700 unrelated individuals of Slavic origin from the Central Russia, including 991 CAD patients and 709 healthy controls were examined. Sixteen lipid-associated GWAS loci were selected from European studies and genotyped using the MassArray-4 system. The polymorphisms were associated with plasma lipids such as total cholesterol (rs12328675, rs4846914, rs55730499, and rs838880), LDL-cholesterol (rs3764261, rs55730499, rs1689800, and rs838880), HDL-cholesterol (rs3764261) as well as carotid intima-media thickness/CIMT (rs12328675, rs11220463, and rs1689800). Polymorphisms such as rs4420638 of APOC1 (p = 0.009), rs55730499 of LPA (p = 0.0007), rs3136441 of F2 (p < 0.0001), and rs6065906 of PLTP (p = 0.002) showed significant associations with the risk of CAD, regardless of sex, age, and body mass index. A majority of the observed associations were successfully replicated in large independent cohorts. Bioinformatics analysis allowed establishing (1) phenotype-specific and shared epistatic gene–gene and gene–smoking interactions contributing to all studied cardiovascular phenotypes; (2) lipid-associated GWAS loci might be allele-specific binding sites for transcription factors from gene regulatory networks controlling multifaceted molecular mechanisms of atherosclerosis.

Fuchs Endothelial Corneal Dystrophy associated risk variant, rs3768617 in LAMC1 shows allele specific binding of GFI1B

AimsTo investigate the genetic and functional association of an intronic variant of LAMC1, rs3768617 with Fuchs endothelial corneal dystrophy (FECD) in the Indian population. MethodsBlood samples were collected from age and sex matched 356 controls and 120 FECD patients after a detailed assessment via specular microscopy. Genomic DNA was extracted and genotyping was done by fluorescence based capillary electrophoresis. The genetic association of rs3768617 polymorphisms was computed by the chi-square (χ2) test. Bioinformatics studies were performed to find the allele specific binding of different transcription factors in the region of rs3768617 and functional evaluation assessed by luciferase assay followed by Electrophoretic Mobility Shift Assay (EMSA) and Chromatin Immunoprecipitation assay (ChIP). Immunofluorescence assay was carried out to check for any differential expression of GFI1B between control and FECD endothelium samples. ResultsSNP rs3768617 {chr1:183123365 (GRCh38.p13)} was found to be genetically associated with FECD in Indian population (p = 2.646 × 10−8). Luciferase assay suggested that the rs3768617 locus has a regulatory role. In silico analysis showed that the transcription factor, GFI1B binds to the risk allele ‘G' of rs3768617, but not to the protective allele ‘A’ which was also experimentally validated by EMSA. High enrichment of DNA flanking the surrounding region of rs3768617 was also found in presence of GFI1B specific antibody in ChIP assay. There was a 0.63 fold decrease in GFI1B expression in FECD affected corneal endothelium compared to control endothelium. ConclusionsThe genetic association of rs3768617 in LAMC1 with FECD pathogenesis is mediated by GFI1B, thus finding the functional role of LAMC1 in FECD pathogenesis.

Functional analysis of the 1p34.3 risk locus implicates GNL2 in high-grade serous ovarian cancer

The high-grade serous ovarian cancer (HGSOC) risk locus at chromosome 1p34.3 resides within a frequently amplified genomic region signifying the presence of an oncogene. Here, we integrate in silico variant-to-function analysis with functional studies to characterize the oncogenic potential of candidate genes in the 1p34.3 locus. Fine mapping of genome-wide association statistics identified candidate causal SNPs local to H3K27ac-demarcated enhancer regions that exhibit allele-specific binding for CTCF in HGSOC and normal fallopian tube secretory epithelium cells (FTSECs). SNP risk associations colocalized with eQTL for six genes (DNALI1, GNL2, RSPO1, SNIP1, MEAF6, and LINC01137) that are more highly expressed in carriers of the risk allele, and three (DNALI1, GNL2, and RSPO1) were upregulated in HGSOC compared to normal ovarian surface epithelium cells and/or FTSECs. Increased expression of GNL2 and MEAF6 was associated with shorter survival in HGSOC with 1p34.3 amplifications. Despite its activation of β-catenin signaling, RSPO1 overexpression exerted no effects on proliferation or colony formation in our study of ovarian cancer and FTSECs. Instead, GNL2, MEAF6, and SNIP1 silencing impaired in vitro ovarian cancer cell growth. Additionally, GNL2 silencing diminished xenograft tumor formation, whereas overexpression stimulated proliferation and colony formation in FTSECs. GNL2 influences 60S ribosomal subunit maturation and global protein synthesis in ovarian cancer and FTSECs, providing a potential mechanism of how GNL2 upregulation might promote ovarian cancer development and mediate genetic susceptibility of HGSOC.

Polymorphic Variants in the GRK5 Gene Promoter Are Associated With Diastolic Dysfunction in Coronary Artery Bypass Graft Surgery Patients.

The G-protein-coupled receptor kinase 5 (GRK5) is a mediator of cardiovascular homeostasis and participates in inflammation and cardiac fibrosis, both being involved in the development of diastolic dysfunction (DD). While mechanisms of transcriptional regulation of the GRK5 promoter are unclear, we tested the hypotheses, that (1) GRK5 expression varies depending on functional single nucleotide polymorphisms (SNPs) in the GRK5 promoter and (2) this is associated with DD in patients undergoing coronary artery bypass graft (CABG) surgery. We amplified and sequenced the GRK5 promoter followed by cloning, reporter assays, and electrophoretic mobility shift assays (EMSA). GRK5 messenger ribonucleic acid (mRNA) expression was determined in right atrial tissue sampled from 50 patients undergoing CABG surgery. In another prospective study, GRK5 genotypes were associated with determinants of diastolic function using transesophageal echocardiography in 255 patients with CABG with normal systolic left ventricular (LV) function. Specifically, we measured ejection fraction (EF), transmitral Doppler early filling velocity (E), tissue Doppler early diastolic lateral mitral annular velocity (E' lateral), and calculated E/E', E' norm and the difference of E' lateral and E' norm to account for age-related changes in diastolic function. We identified 6 SNPs creating 3 novel haplotypes with the greatest promoter activation in haplotype tagging (ht) SNP T(-678)C T-allele constructs (P < .001). EMSAs showed allele-specific transcription factor binding proving functional activity. GRK5 mRNA expression was greatest in TT genotypes (TT: 131 fg/µg [95% CI, 108-154]; CT: 109 [95% confidence interval {CI}, 93-124]; CC: 83 [95% CI, 54-112]; P = .012). Moreover, GRK5 genotypes were significantly associated with determinants of diastolic function. Grading of DD revealed more grade 3 patients in TT compared to CT and CC genotypes (58% vs 38% vs 4%; P = .023). E´ lateral was lowest in TT genotypes (P = .007) and corresponding E/E' measurements showed 1.27-fold increased values in TT versus CC genotypes (P = .01), respectively. While E' norm values were not different between genotypes (P = .182), the difference between E' lateral and E' norm was significantly higher in TT genotypes compared to CC and CT genotypes (-1.2 [interquartile range {IQR}, 2.7], -0.5 [IQR, 3.4], and -0.4 [IQR, 4.2; P = .035], respectively). A functional GRK5 SNP results in allele-dependent differences in GRK5 promoter activity and mRNA expression. This is associated with altered echocardiographic determinants of diastolic function. Thus, SNPs in the GRK5 promoter are associated with altered perioperative diastolic cardiac function. In the future, preoperative testing for these and other SNPs might allow to initiate more specific diagnostic and perioperative pathways to benefit patients at risk.

Genome-wide association and functional interrogation identified a variant at 3p26.1 modulating ovarian cancer survival among Chinese women

Ovarian cancer survival varies considerably among patients, to which germline variation may also contribute in addition to mutational signatures. To identify genetic markers modulating ovarian cancer outcome, we performed a genome-wide association study in 2130 Chinese ovarian cancer patients and found a hitherto unrecognized locus at 3p26.1 to be associated with the overall survival (Pcombined = 8.90 × 10−10). Subsequent statistical fine-mapping, functional annotation, and eQTL mapping prioritized a likely casual SNP rs9311399 in the non-coding regulatory region. Mechanistically, rs9311399 altered its enhancer activity through an allele-specific transcription factor binding and a long-range interaction with the promoter of a lncRNA BHLHE40-AS1. Deletion of the rs9311399-associated enhancer resulted in expression changes in several oncogenic signaling pathway genes and a decrease in tumor growth. Thus, we have identified a novel genetic locus that is associated with ovarian cancer survival possibly through a long-range gene regulation of oncogenic pathways.

Abstract 11643: ACTN4 Regulates HLA-DPA1-Related Pulmonary Endothelial Pathophenotypes in Pulmonary Arterial Hypertension Through Allele-Specific Binding to Functional SNP rs9277336

Introduction: Pulmonary arterial hypertension (PAH) is a progressive disease characterized by increased pulmonary arterial pressure and vascular resistance. This vascular injury promotes immune dysregulation, further worsening the vasculopathy. A genome-wide association study identified a single nucleotide polymorphism (SNP) associated with PAH risk (rs2856830) at Human Leukocyte Antigen DPA1 (HLA-DPA1) locus, a gene known to be involved in adaptive immunity. However, little is known about the mechanism underlying this clinical association. Hypothesis: Allele-specific binding of transcription factors to a disease-associated functional SNP (fSNP) controls HLA-DPA1 expression, driving pulmonary vascular cell dysfunction and PAH. Methods: Oligonucleotides containing alleles of rs2856830 and eight SNPs in linkage disequilibrium (LD) were subjected to electrophoretic mobility shift assays (EMSA) to identify fSNPs characterized by allele-imbalanced binding to human pulmonary arterial endothelial cells (PAEC) nuclear proteins. SNP genotyping was done using blood from PAH (n=134) and control patients (n=131) at UPMC. Transcriptomics, angiogenesis, and migration were assessed in siRNA mediated HLA-DPA1 and Alpha Actinin 4 (ACTN4) knockdown (KD) in PAECs. Results: Through EMSA, rs9277336 (LD r 2 =0.87) was identified as a fSNP in HLA-DPA1 locus. SNP pulldown and immunoblotting confirmed allele-specific binding of nuclear ACTN4 protein to rs9277336 risk allele A (3.26 ± 0.49, p=0.001). SNP genotyping revealed a 2.04 ± 0.67 fold enrichment of risk allele A in PAH patients (χ 2 test, p=0.025). ACTN4 KD in cultured human PAECs decreased HLA-DPA1 transcript (0.30 ± 0.10, p&lt;0.001) . Transcriptomics revealed overlap of several cell structure and immune system-related pathways by HLA-DPA1 and ACTN4 KD. HLA-DPA1 and ACTN4 KD affected PAH associated pathophenotypes, including angiogenesis and migration. Conclusions: Our functional analysis identified ACTN4 as a novel regulator to disease-associated gene HLA-DPA1, affecting the endothelial pathophenotypes in PAH. The allele-specific binding to fSNP rs9277336 provides a mechanistic explanation for the clinical association between HLA-DPA1 locus and PAH outcomes.