Non-coding Variants Research Articles

Sudden Cardiac Death and Channelopathies: What Lies behind the Clinical Significance of Rare Splice-Site Alterations in the Genes Involved?

Background and objectives: Sudden cardiac death (SCD) is a natural and unexpected death of cardiac origin that occurs within 1 h from the onset of acute symptoms. The major leading causes of SCD are cardiomyopathies and channelopathies. In this review, we focus on channelopathies, inherited diseases caused by mutations affecting genes encoding membrane ion channels (sodium, potassium or calcium channels) or cellular structures that affect Ca2+ availability. The diagnosis of diseases such as long QT, Brugada syndrome, short QT and catecholaminergic polymorphic ventricular tachycardia (CPVT) is still challenging. Currently, genetic testing and next-generation sequencing allow us to identify many rare alterations. However, some non-coding variants, e.g., splice-site variants, are usually difficult to interpret and to classify. Methods: In our review, we searched for splice-site variants of genes involved in channelopathies, focusing on variants of unknown significance (VUSs) registered on ClinVar up to now. Results: The research led to a high number of splice-site VUSs of genes involved in channelopathies, suggesting the performance of deeper studies. Conclusions: In order to interpret the correlation between variants and pathologies, we discuss experimental studies, such as RNA sequencing and functional analysis of proteins. Unfortunately, as these in vitro analyses cannot always be performed, we draw attention to in silico studies as future perspectives in genetics. This review has the aim of discussing the potential methods of detection and interpretation of VUSs, bringing out the need for a future reclassification of variants with currently unknown significance.

Characterization of non-coding variants associated with transcription factor binding through ATAC-seq-defined footprint QTLs in liver.

Non-coding variants discovered by genome-wide association studies (GWAS) are enriched in regulatory elements harboring transcription factor (TF) binding motifs, strongly suggesting a connection between disease association and the disruption of cis-regulatory sequences. Occupancy of a TF inside a region of open chromatin can be detected in ATAC-seq where bound TFs block the transposase Tn5, leaving a pattern of relatively depleted Tn5 insertions known as a "footprint". Here, we sought to identify variants associated with TF-binding, or "footprint quantitative trait loci" (fpQTLs) in ATAC-seq data generated from 170 human liver samples. We used computational tools to scan the ATAC-seq reads to quantify TF binding likelihood as "footprint scores" at variants derived from whole genome sequencing generated in the same samples. We tested for association between genotype and footprint score and observed 693 fpQTLs associated with footprint-inferred TF binding (FDR < 5%). Given that Tn5 insertion sites are measured with base-pair resolution, we show that fpQTLs can aid GWAS and QTL fine-mapping by precisely pinpointing TF activity within broad trait-associated loci where the underlying causal variant is unknown. Liver fpQTLs were strongly enriched across ChIP-seq peaks, liver expression QTLs (eQTLs), and liver-related GWAS loci, and their inferred effect on TF binding was concordant with their effect on underlying sequence motifs in 80% of cases. We conclude that fpQTLs can reveal causal GWAS variants, define the role of TF binding site disruption in disease and provide functional insights into non-coding variants, ultimately informing novel treatments for common diseases.

PSAP-Genomic-Regions: A Method Leveraging Population Data to Prioritize Coding and Non-Coding Variants in Whole Genome Sequencing for Rare Disease Diagnosis.

The introduction of Next-Generation Sequencing technologies in the clinics has improved rare disease diagnosis. Nonetheless, for very heterogeneous or very rare diseases, more than half of cases still lack molecular diagnosis. Novel strategies are needed to prioritize variants within a single individual. The Population Sampling Probability (PSAP) method was developed to meet this aim but only for coding variants in exome data. Here, we propose an extension of the PSAP method to the non-coding genome called PSAP-genomic-regions. In this extension, instead of considering genes as testing units (PSAP-genes strategy), we use genomic regions defined over the whole genome that pinpoint potential functional constraints. We conceived an evaluation protocol for our method using artificially generated disease exomes and genomes, by inserting coding and non-coding pathogenic ClinVar variants in large data sets of exomes and genomes from the general population. PSAP-genomic-regions significantly improves the ranking of these variants compared to using a pathogenicity score alone. Using PSAP-genomic-regions, more than 50% of non-coding ClinVar variants were among the top 10 variants of the genome. On real sequencing data from six patients with Cerebral Small Vessel Disease and nine patients with male infertility, all causal variants were ranked in the top 100 variants with PSAP-genomic-regions. By revisiting the testing units used in the PSAP method to include non-coding variants, we have developed PSAP-genomic-regions, an efficient whole-genome prioritization tool which offers promising results for the diagnosis of unresolved rare diseases.

LncRNA Pvt1 aggravates cardiomyocyte apoptosis via the microRNA-216/Ccnd3 axis

ObjectiveOur study aims to evaluate the role of long non-coding RNA variant translocation gene Pvt1 in cardiomyocyte apoptosis, as well as the potential targets and mechanisms involved in Pvt1-miRNA-mRNA axis. Methods1.Pvt1 knockdown in cells by transfection with small interfering RNA (si-Pvt1), HL-1 cells were randomly divided into control group, hypoxia group, hypoxia + negative control group and hypoxia + si-Pvt1 group. Apoptosis-related genes expression was detected by Western blot assay, RT-qPCR and Flow cytometry assay. 2.Pvt1 knockdown model (sh-Pvt1) was established by injecting adeno-associated virus (AAV) vector shRNA-Pvt1 into the caudal vein 7 days before myocardial infarction, and echocardiography was used to measure cardiac function 7 days after myocardial infarction induced by ligation of the left anterior descending branch. HE staining was used to evaluate the pathological injury of mouse heart tissue, and the apoptotic protein expression was detected by Western blot. 3.lncRNA-related microRNAs were predicted by bioinformatics tools and further verified by dual luciferase experiment. Western blot analysis was used to identify the expression of apoptotic genes following the simultaneous transfection of si-Pvt1 and miR-216 mimics. Genes differentially expressed in hypoxia + si-NC and hypoxia + si-Pvt1 groups were identified by RNA sequencing. These genes were then compared with the target genes of miR-216 predicted by bioinformatics tools. The gene of interest Ccnd3 was excluded from the analysis. Western blot analysis was used to assess the expression of Apoptosis-related proteins in HL-1 cells co-transfected with miR-216 mimics and overexpressed Ccnd3. Results1. Pvt1 was highly expressed in HL-1 induced by hypoxia, and Pvt1 knockdown can reduce cell apoptosis in hypoxia cells. 2. MI causes myocardial injury in mice, and inhibition of Pvt 11 can improve the cardiac function of mice with myocardial infarction, prevent some inflammatory cell infiltration, and reduce myocardial cell apoptosis. 3. Pvt1 acts as a sponge for miR-216 and promotes the expression of Ccnd3. ConclusionPvt1 may promote myocardial infarct-induced apoptosis through the miR-216/Ccnd3 axis.

Selection on the vascular-remodeling BMPER gene is associated with altitudinal adaptation in an insular lizard

Abstract High altitude imposes several extreme constraints on life, such as low oxygen pressure and high levels of ultraviolet radiation, which require specialized adaptations. Many studies have focused on how endothermic vertebrates respond to these challenging environments, but there is still uncertainty on how ectotherms adapt to these conditions. Here, we used whole-genome sequencing of low-altitude (100–600 m) and high-altitude (3,550 m) populations of the wide-ranging Tenerife lizard Gallotia galloti to uncover signatures of selection for altitudinal adaptation. The studied populations show reduced differentiation, sharing similar patterns of genetic variation. Selective sweep mapping suggests that signatures of adaptation to high altitude are not widespread across the genome, clustering in a relatively small number of genomic regions. One of these regions contains BMPER, a gene involved with vascular remodeling, and that has been associated with hypoxia-induced angiogenic response. By genotyping samples across 2 altitudinal transects, we show that allele frequency changes at this locus are not gradual, but rather show a well-defined shift above ca. 1,900 m. Transcript and protein structure analyses on this gene suggest that putative selection likely acts on noncoding variation. These results underline how low oxygen pressure generates the most consistent selective constraint in high-altitude environments, to which vertebrates with vastly contrasting physiological profiles need to adapt in the context of ongoing climate change.

Retracted: Long Non-Coding RNA Plasmacytoma Variant Translocation 1 (PVT1) Enhances Proliferation, Migration, and Epithelial-Mesenchymal Transition (EMT) of Pituitary Adenoma Cells by Activating β-Catenin, c-Myc, and Cyclin D1 Expression.

The Editors of Medical Science Monitor wish to inform you that the above manuscript has been retracted from publication due to concerns with the credibility and originality of the study, the manuscript content, and the Figure images. Reference: Yihua Zhang, Yang Tan, Hao Wang, Minhui Xu, Lunshan Xu. Long Non-Coding RNA Plasmacytoma Variant Translocation 1 (PVT1) Enhances Proliferation, Migration, and Epithelial-Mesenchymal Transition (EMT) of Pituitary Adenoma Cells by Activating ß-Catenin, c-Myc, and Cyclin D1 Expression. Med Sci Monit, 2019; 25: 7652-7659. DOI: 10.12659/MSM.917110.

Variation within the non-coding genome influences genetic and epigenetic regulation of the human leukocyte antigen genes.

Variation within the non-coding genome may influence the regulation and expression of important genes involved in immune control such as the human leukocyte antigen (HLA) system. Class I and Class II HLA molecules are essential for peptide presentation which is required for T lymphocyte activation. Single nucleotide polymorphisms within non-coding regions of HLA Class I and Class II genes may influence the expression of these genes by affecting the binding of transcription factors and chromatin modeling molecules. Furthermore, an interplay between genetic and epigenetic factors may also influence HLA expression. Epigenetic factors such as DNA methylation and non-coding RNA, regulate gene expression without changing the DNA sequence. However, genetic variation may promote or allow genes to escape regulation by epigenetic factors, resulting in altered expression. The HLA system is central to most diseases, therefore, understanding the role of genetics and epigenetics on HLA regulation will tremendously impact healthcare. The knowledge gained from these studies may lead to novel and cost-effective diagnostic approaches and therapeutic interventions. This review discusses the role of non-coding variants on HLA regulation. Furthermore, we discuss the interplay between genetic and epigenetic factors on the regulation of HLA by evaluating literature based on polymorphisms within DNA methylation and miRNA regulatory sites within class I and Class II HLA genes. We also provide insight into the importance of the HLA non-coding genome on disease, discuss ethnic-specific differences across the HLA region and provide guidelines for future HLA studies.

DNA breathing integration with deep learning foundational model advances genome-wide binding prediction of human transcription factors.

It was previously shown that DNA breathing, thermodynamic stability, as well as transcriptional activity and transcription factor (TF) bindings are functionally correlated. To ascertain the precise relationship between TF binding and DNA breathing, we developed the multi-modal deep learning model EPBDxDNABERT-2, which is based on the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA dynamics model. To train our EPBDxDNABERT-2, we used chromatin immunoprecipitation sequencing (ChIP-Seq) data comprising 690 ChIP-seq experimental results encompassing 161 distinct TFs and 91 human cell types. EPBDxDNABERT-2 significantly improves the prediction of over 660 TF-DNA, with an increase in the area under the receiver operating characteristic (AUROC) metric of up to 9.6% when compared to the baseline model that does not leverage DNA biophysical properties. We expanded our analysis to in vitro high-throughput Systematic Evolution of Ligands by Exponential enrichment (HT-SELEX) dataset of 215 TFs from 27 families, comparing EPBD with established frameworks. The integration of the DNA breathing features with DNABERT-2 foundational model, greatly enhanced TF-binding predictions. Notably, EPBDxDNABERT-2, trained on a large-scale multi-species genomes, with a cross-attention mechanism, improved predictive power shedding light on the mechanisms underlying disease-related non-coding variants discovered in genome-wide association studies.

Nutrient starvation activates ECM remodeling gene enhancers associated with inflammatory bowel disease risk in fibroblasts.

Nutrient deprivation induces a reversible cell cycle arrest state termed quiescence, which often accompanies transcriptional silencing and chromatin compaction. Paradoxically, nutrient deprivation is associated with activated fibroblast states in pathological microenvironments in which fibroblasts drive extracellular matrix (ECM) remodeling to alter tissue environments. The relationship between nutrient deprivation and fibroblast activation remains unclear. Here, we report that serum deprivation extensively activates transcription of ECM remodeling genes in cultured fibroblasts, despite the induction of quiescence. Starvation-induced transcriptional activation accompanied large-scale histone acetylation of putative distal enhancers, but not promoters. The starvation-activated putative enhancers were enriched for non-coding genetic risk variants associated with inflammatory bowel disease (IBD), suggesting that the starvation-activated gene regulatory network may contribute to fibroblast activation in IBD. Indeed, the starvation-activated gene PLAU, encoding uPA serine protease for plasminogen and ECM, was upregulated in inflammatory fibroblasts in the intestines of IBD patients. Furthermore, the starvation-activated putative enhancer at PLAU, which harbors an IBD risk variant, gained chromatin accessibility in IBD patient fibroblasts. This study implicates nutrient deprivation in transcriptional activation of ECM remodeling genes in fibroblasts and suggests nutrient deprivation as a potential mechanism for pathological fibroblast activation in IBD.

Targeted long-read sequencing enriches disease-relevant genomic regions of interest to provide complete Mendelian disease diagnostics.

Despite advances in sequencing technologies, a molecular diagnosis remains elusive in many patients with Mendelian disease. Current short-read clinical sequencing approaches cannot provide chromosomal phase information or epigenetic information without further sample processing, which is not routinely done and can result in an incomplete molecular diagnosis in patients. The ability to provide phased genetic and epigenetic information from a single sequencing run would improve the diagnostic rate of Mendelian conditions. Here, we describe targeted long-read sequencing of Mendelian disease genes (TaLon-SeqMD) using a real-time adaptive sequencing approach. Optimization of bioinformatic targeting enabled selective enrichment of multiple disease-causing regions of the human genome. Haplotype-resolved variant calling and simultaneous resolution of epigenetic base modification could be achieved in a single sequencing run. The TaLon-SeqMD approach was validated in a cohort of 18 individuals with previous genetic testing targeting 373 inherited retinal disease (IRD) genes, yielding the complete molecular diagnosis in each case. This approach was then applied in 2 IRD cases with inconclusive testing, which uncovered noncoding and structural variants that were difficult to characterize by standard short-read sequencing. Overall, these results demonstrate TaLon-SeqMD as an approach to provide rapid phased-variant calling to provide the molecular basis of Mendelian diseases.

Non-coding cause of congenital heart defects: Abnormal RNA splicing with multiple isoforms as a mechanism for heterotaxy

Heterotaxy is a disorder characterized by severe congenital heart defects (CHDs) and abnormal left-right patterning in other thoracic or abdominal organs. Clinical and research-based genetic testing has previously focused on evaluation of coding variants to identify causes of CHDs, leaving non-coding causes of CHDs largely unknown. Variants in the transcription factor zinc finger of the cerebellum 3 (ZIC3) cause X-linked heterotaxy. We identified an X-linked heterotaxy pedigree without a coding variant in ZIC3. Whole-genome sequencing revealed a deep intronic variant (ZIC3 c.1224+3286A>G) predicted to alter RNA splicing. An invitro minigene splicing assay confirmed the variant acts as a cryptic splice acceptor. CRISPR-Cas9 served to introduce the ZIC3 c.1224+3286A>G variant into human embryonic stem cells demonstrating pseudoexon inclusion caused by the variant. Surprisingly, Sanger sequencing of the resulting ZIC3 c.1224+3286A>G amplicons revealed several isoforms, many of which bypass the normal coding sequence of the third exon of ZIC3, causing a disruption of a DNA-binding domain and a nuclear localization signal. Short- and long-read mRNA sequencing confirmed these initial results and identified additional splicing patterns. Assessment of four isoforms determined abnormal functions invitro and invivo while treatment with a splice-blocking morpholino partially rescued ZIC3. These results demonstrate that pseudoexon inclusion in ZIC3 can cause heterotaxy and provide functional validation of non-coding disease causation. Our results suggest the importance of non-coding variants in heterotaxy and the need for improved methods to identify and classify non-coding variation that may contribute to CHDs.

Whole-genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles

AbstractGenetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole-genome sequencing (WGS) data provide better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the National Heart, Lung, and Blood Institute’s Trans-Omics for Precision Medicine (TOPMed) program (n = 32 572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n = 131 340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, 4 are driven by common variants of small effect with reported minor allele frequency (MAF) at least 10 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor 2 conditionally distinct, noncoding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA) contains 7 distinct signals, including 1 novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR = 0.180; MAFEUR = 0.008). Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation.

Genetic advances in neurodevelopmental disorders

Abstract Neurodevelopmental disorders (NDDs) are a group of highly heterogeneous diseases that affect children’s social, cognitive, and emotional functioning. The etiology is complicated with genetic factors playing an important role. During the past decade, large-scale whole exome sequencing (WES) and whole genome sequencing (WGS) have vastly advanced the genetic findings of NDDs. Various forms of variants have been reported to contribute to NDDs, such as de novo mutations (DNMs), copy number variations (CNVs), rare inherited variants (RIVs), and common variation. By far, over 200 high-risk NDD genes have been identified, which are involved in biological processes including synaptic function, transcriptional and epigenetic regulation. In addition, monogenic, oligogenic, polygenetic, and omnigenic models have been proposed to explain the genetic architecture of NDDs. However, the majority of NDD patients still do not have a definitive genetic diagnosis. In the future, more types of risk factors, as well as noncoding variants, are await to be identified, and including their interplay mechanisms are key to resolving the etiology and heterogeneity of NDDs.

The roles and mechanisms of coding and noncoding RNA variations in cancer.

Functional variations in coding and noncoding RNAs are crucial in tumorigenesis, with cancer-specific alterations often resulting from chemical modifications and posttranscriptional processes mediated by enzymes. These RNA variations have been linked to tumor cell proliferation, growth, metastasis, and drug resistance and are valuable for identifying diagnostic or prognostic cancer biomarkers. The diversity of posttranscriptional RNA modifications, such as splicing, polyadenylation, methylation, and editing, is particularly significant due to their prevalence and impact on cancer progression. Additionally, other modifications, including RNA acetylation, circularization, miRNA isomerization, and pseudouridination, are recognized as key contributors to cancer development. Understanding the mechanisms underlying these RNA modifications in cancer can enhance our knowledge of cancer biology and facilitate the development of innovative therapeutic strategies. Targeting these RNA modifications and their regulatory enzymes may pave the way for novel RNA-based therapies, enabling tailored interventions for specific cancer subtypes. This review provides a comprehensive overview of the roles and mechanisms of various coding and noncoding RNA modifications in cancer progression and highlights recent advancements in RNA-based therapeutic applications.

Identification of Noncoding Functional Regulatory Variants of STIM1 Gene in Idiopathic Pulmonary Arterial Hypertension.

STIM1 (stromal interaction molecule 1) regulates store-operated calcium entry and is involved in pulmonary artery vasoconstriction and pulmonary artery smooth muscle cell proliferation, leading to pulmonary arterial hypertension (PAH). Bioinformatics analysis and a 2-stage matched case-control study were conducted to screen for noncoding variants that may potentially affect STIM1 transcriptional regulation in 242 patients with idiopathic PAH and 414 healthy controls. Luciferase reporter assay, real-time quantitative polymerase chain reaction, western blot, 5-ethynyl-2'-deoxyuridine (EdU) assay, and intracellular Ca2+ measurement were performed to study the mechanistic roles of those STIM1 noncoding variants in PAH. Five noncoding variants (rs3794050, rs7934581, rs3750996, rs1561876, and rs3750994) were identified and genotyped using Sanger sequencing. Rs3794050, rs7934581, and rs1561876 were associated with idiopathic PAH (recessive model, all P<0.05). Bioinformatics analysis showed that these 3 noncoding variants possibly affect the enhancer function of STIM1 or the microRNA (miRNA) binding to STIM1. Functional validation performed in HEK293 and pulmonary artery smooth muscle cells demonstrated that the noncoding variant rs1561876-G (STIM1 mutant) had significantly stronger transcriptional activity than the wild-type counterpart, rs1561876-A, by affecting the transcriptional regulatory function of both hsa-miRNA-3140-5p and hsa-miRNA-4766-5p. rs1561876-G enhanced intracellular Ca2+ signaling in human pulmonary artery smooth muscle cells secondary to calcium-sensing receptor activation and promoted proliferation of pulmonary artery smooth muscle cells under both normoxia and hypoxia conditions, suggesting a possible contribution to PAH development. The potential clinical implications of the 3 noncoding variants of STIM1, rs3794050, rs7934581, and rs1561876, are 2-fold, as they may help predict the risk and prognosis of idiopathic PAH and guide investigations on novel therapeutic pathway(s).

Noncoding variants are a rare cause of recessive developmental disorders in trans with coding variants

PurposeIdentifying pathogenic noncoding variants is challenging. A single protein-altering variant is often identified in a recessive gene in individuals with developmental disorders (DD), but the prevalence of pathogenic noncoding “second hits” in trans with these is unknown. MethodsIn 4073 genetically undiagnosed rare-disease trio probands from the 100,000 Genomes project, we identified rare heterozygous protein-altering variants in recessive DD-associated genes. We identified rare noncoding variants on the other haplotype in introns, untranslated regions, promoters, and candidate enhancer regions. We clinically evaluated the top candidates for phenotypic fit and performed functional testing where possible. ResultsWe identified 3761 rare heterozygous loss-of-function or ClinVar pathogenic variants in recessive DD-associated genes in 2430 probands. For 1366 (36.3%) of these, we identified at least 1 rare noncoding variant in trans. Bioinformatic filtering and clinical review, revealed 7 to be a good clinical fit. After detailed characterization, we identified likely diagnoses for 3 probands (in GAA, NPHP3, and PKHD1) and candidate diagnoses in a further 3 (PAH, LAMA2, and IGHMBP2). ConclusionWe developed a systematic approach to uncover new diagnoses involving compound heterozygous coding/noncoding variants and conclude that this mechanism is likely to be a rare cause of DDs.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic risk to childhood-onset asthma

Asthma is a complex disease caused by genetic and environmental factors. Studies show that wheezing during rhinovirus infection correlates with childhood asthma development. Over 150 non-coding risk variants for asthma have been identified, many affecting gene regulation in Tcells, but the effects of most risk variants remain unknown. We hypothesized that airway epithelial cells could also mediate genetic susceptibility to asthma given they are the first line of defense against respiratory viruses and allergens. We integrated genetic data with transcriptomics of airway epithelial cells subject to different stimuli. We demonstrate that rhinovirus infection significantly upregulates childhood-onset asthma-associated genes, particularly in non-ciliated cells. This enrichment is also observed with influenza infection but not with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or cytokine activation. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Rare germline genetic variation in PAX8 transcription factor binding sites and susceptibility to epithelial ovarian cancer.

Common genetic variation throughout the genome together with rare coding variants identified to date explain about a half of the inherited genetic component of epithelial ovarian cancer risk. It is likely that rare variation in the non-coding genome will explain some of the unexplained heritability, but identifying such variants is challenging. The primary problem is lack of statistical power to identifying individual risk variants by association as power is a function of sample size, effect size and allele frequency. Power can be increased by using burden tests which test for association of carriers of any variant in a specified genomic region. This has the effect of increasing the putative effect allele frequency. PAX8 is a transcription factor that plays a critical role in tumour progression, migration and invasion. Furthermore, regulatory elements proximal to target genes of PAX8 are enriched for common ovarian cancer risk variants. We hypothesised that rare variation in PAX8 binding sites are also associated with ovarian cancer risk, but unlikely to be associated with risk of breast, colorectal or endometrial cancer. We have used publicly available, whole-genome sequencing data from the UK 100,000 Genomes Project to evaluate the burden of rare variation in PAX8 binding sites across the genome. Data were available for 522 ovarian cancers, 2,984 breast cancers, 2,696 colorectal cancers, 836 endometrial cancers and 2253 non-cancer controls. Active binding sites were defined using data from multiple PAX8 and H3K27 ChIPseq experiments. We found no association between the burden of rare variation in PAX8 binding sites (defined in several ways) and risk of ovarian, breast or endometrial cancer. An apparent association with colorectal cancer was likely to be a technical artefact as a similar association was also detected for rare variation in random regions of the genome. Despite the null result this study provides a proof-of -principle for using burden testing to identify rare, non-coding germline genetic variation associated with disease. Larger sample sizes available from large-scale sequencing projects together with improved understanding of the function of the non-coding genome will increase the potential of similar studies in the future.

Deep Learning Sequence Models for Transcriptional Regulation.

Deciphering the regulatory code of gene expression and interpreting the transcriptional effects of genome variation are critical challenges in human genetics. Modern experimental technologies have resulted in an abundance of data, enabling the development of sequence-based deep learning models that link patterns embedded in DNA to the biochemical and regulatory properties contributing to transcriptional regulation, including modeling epigenetic marks, 3D genome organization, and gene expression, with tissue and cell-type specificity. Such methods can predict the functional consequences of any noncoding variant in the human genome, even rare or never-before-observed variants, and systematically characterize their consequences beyond what is tractable from experiments or quantitative genetics studies alone. Recently, the development and application of interpretability approaches have led to the identification of key sequence patterns contributing to the predicted tasks, providing insights into the underlying biological mechanisms learned and revealing opportunities for improvement in future models.

G protein–coupled receptor kinase 5 regulates thrombin signaling in platelets

BackgroundOur prior genome-wide association study of thrombin-induced platelet aggregation identified a G protein–coupled receptor kinase 5 (GRK5) noncoding variant (rs10886430-G) that is strongly associated with increased platelet reactivity to thrombin. This variant predisposes to increased risk of stroke, pulmonary embolism, and venous thromboembolism. ObjectivesTo determine role of platelet specific GRK5 in platelet responses to agonists and injury. MethodsPlatelets from GRK5 mutant mice have been shown to have increased thrombin sensitivity, indicating that GRK5 may be a negative regulator of platelet activation. However, this has not been studied in a platelet-specific manner. We therefore used platelet-specific GRK5 mutant mice and models of thrombosis and pulmonary embolism. ResultsWe now demonstrate that mice lacking GRK5 specifically in platelets had a mild increase in thrombin responses in vitro and a shortened time to arterial thrombosis in vivo. In addition, platelet GRK5 mutant mice had increased thrombin but not collagen-induced thrombus burden in a mouse model of pulmonary embolism. ConclusionThese data indicate that platelet GRK5 has a significant role in limiting platelet responses to thrombin.