A sequential methodology for the rapid identification and characterization of breast cancer-associated functional SNPs

Integration of multi-omics transcriptome-wide analysis for the identification of novel therapeutic drug targets in diabetic retinopathy

BACKGROUND

Diabetic retinopathy (DR) is the most important complication of Type 2 Diabetes (T2D) in eyes. Despite its prevalence, the early detection and management of DR continue to pose considerable challenges. Our research aims to elucidate potent drug targets that could facilitate the identification of DR and propel advancements in its therapeutic strategies.

METHODS

A broad multi-omics exploration of DR was presented to decipher the drug targets of DR and proliferative diabetic retinopathy (PDR). Transcriptome-Wide Association Studies (TWAS), fine-mapping and conditional analysis were applied to unearth potential tissue-specific gene associations with DR. Summary Data-based Mendelian Randomization (SMR) provided secondary analysis of high confidence genes. Cis-instrument of druggable genes were extracted from the eQTLGen Consortium and PsychENCODE, facilitating drug-target MR supported by colocalization analysis. Phenome-Wide Association Studies (PheWAS) was conducted on the high confidence genes. Metabolomic and immunomic MR-profiling further augmented our research as complement.

RESULTS

TWAS identified multiple robust genetic loci in both DR and PDR (WFS1, RPS26, and SRPK1) through genetic associations across different tissues. Meanwhile, we have delineated both the commonalities and discrepancies between DR and PDR at the transcriptomic level, represented by DCLRE1B as the hub gene that DR progressed into PDR. SMR revealed 92 key DR-related genes and 55 PDR-related genes. HLA-DQ family genes have a frequent occurrence, while RPS26, WFS1 and SRPK1 were validated as the genetic network's linchpins. Drug-target MR casted ERBB3 and SRPK1 as candidate effector genes for DR and PDR susceptibility. In addition, metabolomics and immunomics analyses also revealed multifaceted pathogenic factors for DR.

CONCLUSIONS

Our research offers targeted therapeutic insights for early-stage DR and facilitates multi-omic comparisons of it and PDR.

Methods for Functional Characterization of Genetic Polymorphisms of Non-Coding Regulatory Regions of the Human Genome

Currently, numerous associations between genetic polymorphisms and various diseases have been characterized through the Genome-Wide Association Studies. Majority of the clinically significant polymorphisms are localized in non-coding regions of the genome. While modern bioinformatic resources make it possible to predict molecular mechanisms that explain influence of the non-coding polymorphisms on gene expression, such hypotheses require experimental verification. This review discusses the methods for elucidating molecular mechanisms underlying dependence of the disease pathogenesis on specific genetic variants within the non-coding sequences. A particular focus is on the methods for identification of transcription factors with binding efficiency dependent on polymorphic variations. Despite remarkable progress in bioinformatic resources enabling prediction of the impact of polymorphisms on the disease pathogenesis, there is still the need for experimental approaches to investigate this issue.

Allele-specific control of rodent and human lncRNA KMT2E-AS1 promotes hypoxic endothelial pathology in pulmonary hypertension

Hypoxic reprogramming of vasculature relies on genetic, epigenetic, and metabolic circuitry, but the control points are unknown. In pulmonary arterial hypertension (PAH), a disease driven by hypoxia inducible factor (HIF)-dependent vascular dysfunction, HIF-2α promoted expression of neighboring genes, long noncoding RNA (lncRNA) histone lysine N-methyltransferase 2E-antisense 1 (KMT2E-AS1) and histone lysine N-methyltransferase 2E (KMT2E). KMT2E-AS1 stabilized KMT2E protein to increase epigenetic histone 3 lysine 4 trimethylation (H3K4me3), driving HIF-2α-dependent metabolic and pathogenic endothelial activity. This lncRNA axis also increased HIF-2α expression across epigenetic, transcriptional, and posttranscriptional contexts, thus promoting a positive feedback loop to further augment HIF-2α activity. We identified a genetic association between rs73184087, a single-nucleotide variant (SNV) within a KMT2E intron, and disease risk in PAH discovery and replication patient cohorts and in a global meta-analysis. This SNV displayed allele (G)-specific association with HIF-2α, engaged in long-range chromatin interactions, and induced the lncRNA-KMT2E tandem in hypoxic (G/G) cells. In vivo, KMT2E-AS1 deficiency protected against PAH in mice, as did pharmacologic inhibition of histone methylation in rats. Conversely, forced lncRNA expression promoted more severe PH. Thus, the KMT2E-AS1/KMT2E pair orchestrates across convergent multi-ome landscapes to mediate HIF-2α pathobiology and represents a key clinical target in pulmonary hypertension.

UTP23 Is a Promising Prognostic Biomarker and Is Associated with Immune Infiltration in Breast Cancer

Breast cancer is one of the malignant tumors with a high incidence and mortality rate among women worldwide, and its prevalence is increasing year by year, posing a serious health risk to women. UTP23 (UTP23 Small Subunit Processome Component) is a nucleolar protein that is essential for ribosome production. As we all know, disruption of ribosome structure and function results in improper protein function, affecting the body's normal physiological processes and promoting cancer growth. However, little research has shown a connection between UTP23 and cancer. We analyzed the mRNA expression of UTP23 in normal tissue and breast cancer using The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) database, and the protein expression of UTP23 using The Human Protein Atlas (HPA) database. Next, we examined the relationship between UTP23 high expression and Overall Survival (OS) using Kaplan-Meier Plotters and enriched 980 differentially expressed genes in UTP23 high and low expression samples using GO/KEGG and GSEA to identify potential biological functions of UTP23 and signaling pathways that it might influence. Finally, we also investigated the relationship between UTP23 and immune infiltration and examined the effect of UTP23 on the proliferation of human breast cancer cell lines by knocking down UTP23. We found that UTP23 levels in breast cancer patient samples were noticeably greater than those in healthy individuals and that high UTP23 levels were strongly linked with poor prognoses (P = 0.008). Functional enrichment analysis revealed that UTP23 expression was connected to the humoral immune response. Besides, UTP23 expression was found to be positively correlated with immune cell infiltration. Furthermore, UTP23 knockdown has been shown to inhibit the proliferation of human breast cancer cells MDA-MB-231 and HCC-1806. Taken together, our study demonstrated that UTP23 is a promising target in detecting and treating breast cancer and is intimately linked to immune infiltration.

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

Single Nucleotide Polymorphism rs9277336 Controls the Nuclear Alpha Actinin 4-Human Leukocyte Antigen-DPA1 Axis and Pulmonary Endothelial Pathophenotypes in Pulmonary Arterial Hypertension

Background Pulmonary arterial hypertension (PAH) is a complex, fatal disease where disease severity has been associated with the single nucleotide polymorphism (SNP) rs2856830, located near the human leukocyte antigen DPA1 (HLA-DPA1) gene. We aimed to define the genetic architecture of functional variants associated with PAH disease severity by identifying allele-specific binding transcription factors and downstream targets that control endothelial pathophenotypes and PAH. Methods and Results Electrophoretic mobility shift assays of oligonucleotides containing SNP rs2856830 and 8 SNPs in linkage disequilibrium revealed functional SNPs via allele-imbalanced binding to human pulmonary arterial endothelial cell nuclear proteins. DNA pulldown proteomics identified SNP-binding proteins. SNP genotyping and clinical correlation analysis were performed in 84 patients with PAH at University of Pittsburgh Medical Center and in 679 patients with PAH in the All of Us database. SNP rs9277336 was identified as a functional SNP in linkage disequilibrium (r2>0.8) defined by rs2856830, and the minor allele was associated with decreased hospitalizations and improved cardiac output in patients with PAH, an index of disease severity. SNP pulldown proteomics showed allele-specific binding of nuclear ACTN4 (alpha actinin 4) protein to rs9277336 minor allele. Both ACTN4 and HLA-DPA1 were downregulated in pulmonary endothelium in human patients and rodent models of PAH. Via transcriptomic and phenotypic analyses, knockdown of HLA-DPA1 phenocopied knockdown of ACTN4, both similarly controlling cell structure pathways, immune pathways, and endothelial dysfunction. Conclusions We defined the pathogenic activity of functional SNP rs9277336, entailing the allele-specific binding of ACTN4 and controlling expression of the neighboring HLA-DPA1 gene. Through inflammatory or genetic means, downregulation of this ACTN4-HLA-DPA1 regulatory axis promotes endothelial pathophenotypes, providing a mechanistic explanation for the association between this SNP and PAH outcomes.

Multiplexed functional genomic assays to decipher the noncoding genome

Linkage disequilibrium and the incomplete regulatory annotation of the noncoding genome complicates the identification of functional noncoding genetic variants and their causal association with disease. Current computational methods for variant prioritization have limited predictive value, necessitating the application of highly parallelized experimental assays to efficiently identify functional noncoding variation. Here, we summarize two distinct approaches, massively parallel reporter assays and CRISPR-based pooled screens and describe their flexible implementation to characterize human noncoding genetic variation at unprecedented scale. Each approach provides unique advantages and limitations, highlighting the importance of multimodal methodological integration. These multiplexed assays of variant effects are undoubtedly poised to play a key role in the experimental characterization of noncoding genetic risk, informing our understanding of the underlying mechanisms of disease-associated loci and the development of more robust predictive classification algorithms.

Genetic susceptibility loci for Chlamydia trachomatis endometrial infection influence expression of genes involved in T cell function, tryptophan metabolism and epithelial integrity

Objectives

Identify genetic loci of enhanced susceptibility to Chlamydial trachomatis (Ct) upper genital tract infection in women.

Methods

We performed an integrated analysis of DNA genotypes and blood-derived mRNA profiles from 200 Ct-exposed women to identify expression quantitative trait loci (eQTL) and determine their association with endometrial chlamydial infection using a mediation test. We further evaluated the effect of a lead eQTL on the expression of CD151 by immune cells from women with genotypes associated with low and high whole blood expression of CD151, respectively.

Results

We identified cis-eQTLs modulating mRNA expression of 81 genes (eGenes) associated with altered risk of ascending infection. In women with endometrial infection, eGenes involved in proinflammatory signaling were upregulated. Downregulated eGenes included genes involved in T cell functions pivotal for chlamydial control. eGenes encoding molecules linked to metabolism of tryptophan, an essential chlamydial nutrient, and formation of epithelial tight junctions were also downregulated in women with endometrial infection. A lead eSNP rs10902226 was identified regulating CD151, a tetrospanin molecule important for immune cell adhesion and migration and T cell proliferation. Further in vitro experiments showed that women with a CC genotype at rs10902226 had reduced rates of endometrial infection with increased CD151 expression in whole blood and T cells when compared to women with a GG genotype.

Conclusions

We discovered genetic variants associated with altered risk for Ct ascension. A lead eSNP for CD151 is a candidate genetic marker for enhanced CD4 T cell function and reduced susceptibility.

Post-GWAS functional analysis identifies CUX1 as a regulator of p16INK4a and cellular senescence

Accumulation of senescent cells with age is an important driver of aging and age-related diseases. However, the mechanisms and signaling pathways that regulate senescence remain elusive. In this report, we performed post-genome-wide association studies (GWAS) functional studies on the CDKN2A/B locus, a locus known to be associated with multiple age-related diseases and overall human lifespan. We demonstrate that transcription factor CUX1 (Cut-Like Homeobox 1) specifically binds to an atherosclerosis-associated functional single-nucleotide polymorphism (fSNP) (rs1537371) within the locus and regulates the CDKN2A/B-encoded proteins p14ARF, p15INK4b and p16INK4a and the antisense noncoding RNA in the CDK4 (INK4) locus (ANRIL) in endothelial cells (ECs). Endothelial CUX1 expression correlates with telomeric length and is induced by both DNA-damaging agents and oxidative stress. Moreover, induction of CUX1 expression triggers both replicative and stress-induced senescence via activation of p16INK4a expression. Thus, our studies identify CUX1 as a regulator of p16INK4a-dependent endothelial senescence and a potential therapeutic target for atherosclerosis and other age-related diseases.

Coupling high-throughput mapping with proteomics analysis delineates cis-regulatory elements at high resolution

Growing evidence suggests that functional cis-regulatory elements (cis-REs) not only exist in epigenetically marked but also in unmarked sites of the human genome. While it is already difficult to identify cis-REs in the epigenetically marked sites, interrogating cis-REs residing within the unmarked sites is even more challenging. Here, we report adapting Reel-seq, an in vitro high-throughput (HTP) technique, to fine-map cis-REs at high resolution over a large region of the human genome in a systematic and continuous manner. Using Reel-seq, as a proof-of-principle, we identified 408 candidate cis-REs by mapping a 58 kb core region on the aging-related CDKN2A/B locus that harbors p16INK4a. By coupling Reel-seq with FREP-MS, a proteomics analysis technique, we characterized two cis-REs, one in an epigenetically marked site and the other in an epigenetically unmarked site. These elements are shown to regulate the p16INK4a expression over an ∼100 kb distance by recruiting the poly(A) binding protein PABPC1 and the transcription factor FOXC2. Downregulation of either PABPC1 or FOXC2 in human endothelial cells (ECs) can induce the p16INK4a-dependent cellular senescence. Thus, we confirmed the utility of Reel-seq and FREP-MS analyses for the systematic identification of cis-REs at high resolution over a large region of the human genome.

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.

A Panel of rSNPs Demonstrating Allelic Asymmetry in Both ChIP-seq and RNA-seq Data and the Search for Their Phenotypic Outcomes through Analysis of DEGs

Currently, the detection of the allele asymmetry of gene expression from RNA-seq data or the transcription factor binding from ChIP-seq data is one of the approaches used to identify the functional genetic variants that can affect gene expression (regulatory SNPs or rSNPs). In this study, we searched for rSNPs using the data for human pulmonary arterial endothelial cells (PAECs) available from the Sequence Read Archive (SRA). Allele-asymmetric binding and expression events are analyzed in paired ChIP-seq data for H3K4me3 mark and RNA-seq data obtained for 19 individuals. Two statistical approaches, weighted z-scores and predicted probabilities, were used to improve the efficiency of finding rSNPs. In total, we identified 14,266 rSNPs associated with both allele-specific binding and expression. Among them, 645 rSNPs were associated with GWAS phenotypes; 4746 rSNPs were reported as eQTLs by GTEx, and 11,536 rSNPs were located in 374 candidate transcription factor binding motifs. Additionally, we searched for the rSNPs associated with gene expression using an SRA RNA-seq dataset for 281 clinically annotated human postmortem brain samples and detected eQTLs for 2505 rSNPs. Based on these results, we conducted Gene Ontology (GO), Disease Ontology (DO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and constructed the protein-protein interaction networks to represent the top-ranked biological processes with a possible contribution to the phenotypic outcome.

Comparative Proteomic Characterization of Ventral Hippocampus in Susceptible and Resilient Rats Subjected to Chronic Unpredictable Stress

Chronic stress is an essential factor leading to depression. However, there exist individual differences in people exposed to the same stressful stimuli. Some people display negative psychology and behavior, while others are normal. Given the importance of individual difference, finding differentially expressed proteins in stress-resistant and stress-susceptible groups has great significance for the study of pathogenesis and treatment of depression. In this study, stress-susceptible rats and stress-resilient rats were first distinguished by sucrose preference test. These stress-susceptible rats also displayed depression-like behaviors in forced swimming test and open field test. Then, we employed label-free quantitative proteomics to analyze proteins in the ventral hippocampus. There were 4,848 proteins totally identified. Based on statistical analysis, we found 276 differentially expressed proteins. Bioinformatics analysis revealed that the biological processes of these differential proteins were related to mitochondrion organization, protein localization, coenzyme metabolic process, cerebral cortex tangential migration, vesicle-mediated transport, and so on. The KEGG pathways were mainly involved in metabolic pathways, axon guidance, autophagy, and tight junction. Furthermore, we ultimately found 20 stress-susceptible proteins and two stress-resilient proteins. These stress-related proteins could not only be potential biomarkers for depression diagnosis but also contribute to finding new therapeutic targets and providing personalized medicine.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

Post-GWAS functional studies reveal an RA-associated CD40-induced NF-kB signal transduction and transcriptional regulation network targeted by class II HDAC inhibitors

Currently, it remains difficult to identify which single nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) are functional and how various functional SNPs (fSNPs) interact and contribute to disease susceptibility. GWAS have identified a CD40 locus that is associated with rheumatoid arthritis (RA). We previously used two techniques developed in our laboratory, single nucleotide polymorphism-next-generation sequencing (SNP-seq) and flanking restriction enhanced DNA pulldown-mass spectrometry (FREP-MS), to determine that the RA risk gene RBPJ regulates CD40 expression via a fSNP at the RA-associated CD40 locus. In the present work, by applying the same approach, we report the identification of six proteins that regulate RBPJ expression via binding to two fSNPs on the RA-associated RBPJ locus. Using these findings, together with the published data, we constructed an RA-associated signal transduction and transcriptional regulation network (STTRN) that functionally connects multiple RA-associated risk genes via transcriptional regulation networks (TRNs) linked by CD40-induced nuclear factor kappa B (NF-kB) signaling. Remarkably, this STTRN provides insight into the potential mechanism of action for the histone deacetylase inhibitor givinostat, an approved therapy for systemic juvenile idiopathic arthritis. Thus, the generation of disease-associated STTRNs based on post-GWAS functional studies is demonstrated as a novel and effective approach to apply GWAS for mechanistic studies and target identification.

Transcriptional Regulation of CD40 Expression by 4 Ribosomal Proteins via a Functional SNP on a Disease-Associated CD40 Locus

Previously, using FREP-MS, we identified a protein complex including eight proteins that specifically bind to the functional SNP (fSNP) rs6032664 at a CD40 locus associated with autoimmune diseases. Among these eight proteins, four are ribosomal proteins RPL26, RPL4, RPL8, and RPS9 that normally make up the ribosomal subunits involved in the cellular process of protein translation. So far, no publication has shown these ribosomal proteins function as transcriptional regulators. In this work, we demonstrate that four ribosomal proteins: RPL26, RPL4, RPL8, and RPS9 are bona fide CD40 transcriptional regulators via binding to rs6032664. In addition, we show that suppression of CD40 expression by RPL26 RNAi knockdown inactivates NF-κB p65 by dephosphorylation via NF-κB signaling pathway in fibroblast-like synoviocytes (FLS), which further reduces the transcription of disease-associated risk genes such as STAT4, CD86, TRAF1 and ICAM1 as the direct targets of NF-κB p65. Based on these findings, a disease-associated risk gene transcriptional regulation network (TRN) is generated, in which decreased expression of, at least, RPL26 results in the downregulation of risk genes: STAT4, CD86, TRAF1 and ICAM1, as well as the two proinflammatory cytokines: IL1β and IL6 via CD40-induced NF-κB signaling. We believe that further characterization of this disease-associated TRN in the CD40-induced NF-κB signaling by identifying both the upstream and downstream regulators will potentially enable us to identify the best targets for drug development.