Long-Range Modulation of PAG1 Expression by 8q21 Allergy Risk Variants

Identification of the shared genetic architecture underlying seven autoimmune diseases with GWAS summary statistics

Background

The common clinical symptoms and immunopathological mechanisms have been observed among multiple autoimmune diseases (ADs), but the shared genetic etiology remains unclear.

Methods

GWAS summary statistics of seven ADs were downloaded from Open Targets Genetics and Dryad. Linkage disequilibrium score regression (LDSC) was applied to estimate overall genetic correlations, bivariate causal mixture model (MiXeR) was used to qualify the polygenic overlap, and stratified-LDSC partitioned heritability to reveal tissue and cell type specific enrichments. Ultimately, we conducted a novel adaptive association test called MTaSPUsSet for identifying pleiotropic genes.

Results

The high heritability of seven ADs ranged from 0.1228 to 0.5972, and strong genetic correlations among certain phenotypes varied between 0.185 and 0.721. There was substantial polygenic overlap, with the number of shared SNPs approximately 0.03K to 0.21K. The specificity of SNP heritability was enriched in the immune/hematopoietic related tissue and cells. Furthermore, we identified 32 pleiotropic genes associated with seven ADs, 23 genes were considered as novel genes. These genes were involved in several cell regulation pathways and immunologic signatures.

Conclusion

We comprehensively explored the shared genetic architecture across seven ADs. The findings progress the exploration of common molecular mechanisms and biological processes involved, and facilitate understanding of disease etiology.

Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity

Disease-associated single-nucleotide polymorphisms (SNPs) generally do not implicate target genes, as most disease SNPs are regulatory. Many SNP-to-gene (S2G) linking strategies have been developed to link regulatory SNPs to the genes that they regulate in cis. Here, we developed a heritability-based framework for evaluating and combining different S2G strategies to optimize their informativeness for common disease risk. Our optimal combined S2G strategy (cS2G) included seven constituent S2G strategies and achieved a precision of 0.75 and a recall of 0.33, more than doubling the recall of any individual strategy. We applied cS2G to fine-mapping results for 49 UK Biobank diseases/traits to predict 5,095 causal SNP-gene-disease triplets (with S2G-derived functional interpretation) with high confidence. We further applied cS2G to provide an empirical assessment of disease omnigenicity; we determined that the top 1% of genes explained roughly half of the SNP heritability linked to all genes and that gene-level architectures vary with variant allele frequency.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

A comparison of transcriptome analysis methods with reference genome

Background

The application of RNA-seq technology has become more extensive and the number of analysis procedures available has increased over the past years. Selecting an appropriate workflow has become an important issue for researchers in the field.

Methods

In our study, six popular analytical procedures/pipeline were compared using four RNA-seq datasets from mouse, human, rat, and macaque, respectively. The gene expression value, fold change of gene expression, and statistical significance were evaluated to compare the similarities and differences among the six procedures. qRT-PCR was performed to validate the differentially expressed genes (DEGs) from all six procedures.

Results

Cufflinks-Cuffdiff demands the highest computing resources and Kallisto-Sleuth demands the least. Gene expression values, fold change, p and q values of differential expression (DE) analysis are highly correlated among procedures using HTseq for quantification. For genes with medium expression abundance, the expression values determined using the different procedures were similar. Major differences in expression values come from genes with particularly high or low expression levels. HISAT2-StringTie-Ballgown is more sensitive to genes with low expression levels, while Kallisto-Sleuth may only be useful to evaluate genes with medium to high abundance. When the same thresholds for fold change and p value are chosen in DE analysis, StringTie-Ballgown produce the least number of DEGs, while HTseq-DESeq2, -edgeR or -limma generally produces more DEGs. The performance of Cufflinks-Cuffdiff and Kallisto-Sleuth varies in different datasets. For DEGs with medium expression levels, the biological verification rates were similar among all procedures.

Conclusion

Results are highly correlated among RNA-seq analysis procedures using HTseq for quantification. Difference in gene expression values mainly come from genes with particularly high or low expression levels. Moreover, biological validation rates of DEGs from all six procedures were similar for genes with medium expression levels. Investigators can choose analytical procedures according to their available computer resources, or whether genes of high or low expression levels are of interest. If computer resources are abundant, one can utilize multiple procedures to obtain the intersection of results to get the most reliable DEGs, or to obtain a combination of results to get a more comprehensive DE profile for transcriptomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08465-0.

Predict long-range enhancer regulation based on protein-protein interactions between transcription factors

Long-range regulation by distal enhancers plays critical roles in cell-type specific transcriptional programs. Computational predictions of genome-wide enhancer-promoter interactions are still challenging due to limited accuracy and the lack of knowledge on the molecular mechanisms. Based on recent biological investigations, the protein-protein interactions (PPIs) between transcription factors (TFs) have been found to participate in the regulation of chromatin loops. Therefore, we developed a novel predictive model for cell-type specific enhancer-promoter interactions by leveraging the information of TF PPI signatures. Evaluated by a series of rigorous performance comparisons, the new model achieves superior performance over other methods. The model also identifies specific TF PPIs that may mediate long-range regulatory interactions, revealing new mechanistic understandings of enhancer regulation. The prioritized TF PPIs are associated with genes in distinct biological pathways, and the predicted enhancer-promoter interactions are strongly enriched with cis-eQTLs. Most interestingly, the model discovers enhancer-mediated trans-regulatory links between TFs and genes, which are significantly enriched with trans-eQTLs. The new predictive model, along with the genome-wide analyses, provides a platform to systematically delineate the complex interplay among TFs, enhancers and genes in long-range regulation. The novel predictions also lead to mechanistic interpretations of eQTLs to decode the genetic associations with gene expression.

Influence of pentoxifylline on gene expression of PAG1/ miR-1206/ SNHG14 in ischemic heart disease

The regulation by immune checkpoint is able to prevent excessive tissue damage caused by ischemia reperfusion (I/R); therefore, the study aims to investigate the behavior of phosphoprotein associated with glycosphingolipid-enriched microdomains 1 (PAG1) mRNA, miR-1206 and small nucleolar RNA host gene 14 (SNHG14) during I/R and intake of pentoxifylline (PTX) as a protective drug. The relative expression level of PAG1/miR-1206/SNHG14 was determined by qRT-PCR. Cardiac tissue levels of cytotoxic T-lymphocyte associated antigen 4 (CTLA4) and PAG1 protein expression were determined by ELISA technique. The regulatory T cells achieved by the flow cytometry. The results found that the relative expression of SNHG14 was significantly upregulated in I/R, but suppressed in PTX treated groups with enhancement of the relative expression level of miR-1206. The gene and protein expression of PAG1 were downregulated with effective doses of PTX. The results showed that (30 and 40 mg/kg bwt) PTX dose suppressed the CTLA4 development significantly. The mean of the regulatory T cell in PTX protective groups is significantly reduced at (p < 0.001) in a comparison with I/R group. Spearman's correlation analysis revealed a significant negative correlation between SNHG14 and miR-1206, but a significant positive correlation between SNHG14 and PAG1 in I/R heart tissue. The results indicated that miR-1206 and SNHG14 can be used as biomarkers with perfect sensitivity and specificity. Using PTX reduced cardiac tissue damage. SNHG14 and miR-1206 can be used as a diagnostic tool in I/R.

•

Positive correlation between SNHG14 and PAG1 relative expression in I/R heart tissues.

•

SNHG14 and miR-1206 can be used as a diagnostic tool in ischemia reperfusion.

•

Negative correlation between SNHG14 and miR-1206.

•

Using pentoxifylline as a protective drug renders cardiac tissues more resistance to ischemia.

Common variants in signaling transcription-factor-binding sites drive phenotypic variability in red blood cell traits

Genome-wide association studies identify genomic variants associated with human traits and diseases. Most trait-associated variants are located within cell-type-specific enhancers, but the molecular mechanisms governing phenotypic variation are less well understood. Here, we show that many enhancer variants associated with red blood cell (RBC) traits map to enhancers that are co-bound by lineage-specific master transcription factors (MTFs) and signaling transcription factors (STFs) responsive to extracellular signals. The majority of enhancer variants reside on STF and not MTF motifs, perturbing DNA binding by various STFs (BMP/TGF-β-directed SMADs or WNT-induced TCFs) and affecting target gene expression. Analyses of engineered human blood cells and expression quantitative trait loci verify that disrupted STF binding leads to altered gene expression. Our results propose that the majority of the RBC-trait-associated variants that reside on transcription-factor-binding sequences fall in STF target sequences, suggesting that the phenotypic variation of RBC traits could stem from altered responsiveness to extracellular stimuli.

The Role of Noncoding Variants in Heritable Disease

The genetic basis of disease has largely focused on coding regions. However, it has become clear that a large proportion of the noncoding genome is functional and harbors genetic variants that contribute to disease etiology. Here, we review recent examples of inherited noncoding alterations that are responsible for Mendelian disorders or act to influence complex traits. We explore both rare and common genetic variants and discuss the wide range of mechanisms by which they affect gene regulation to promote disease. We also debate the challenges and progress associated with identifying and interpreting the functional and clinical significance of genetic variation in the context of the noncoding regulatory landscape.

PAG1 limits allergen-induced type 2 inflammation in the murine lung

BACKGROUND

Phosphoprotein associated with glycosphingolipid-enriched microdomains 1 (PAG1) is a transmembrane adaptor protein that affects immune receptor signaling in T and B cells. Evidence from genome-wide association studies of asthma suggests that genetic variants that regulate the expression of PAG1 are associated with asthma risk. However, it is not known whether PAG1 expression is causally related to asthma pathophysiology. Here, we investigated the role of PAG1 in a preclinical mouse model of house dust mite (HDM)-induced allergic sensitization and allergic airway inflammation.

METHODS

Pag1-deficient (Pag1^-/- ) and wild-type (WT) mice were sensitized or sensitized/challenged to HDM, and hallmark features of allergic inflammation were assessed. The contribution of T cells was assessed through depletion (anti-CD4 antibody) and adoptive transfer studies.

RESULTS

Type 2 inflammation (eosinophilia, eotaxin-2 expression, IL-4/IL-5/IL-13 production, mucus production) in the airways and lungs was significantly increased in HDM sensitized/challenged Pag1^-/- mice compared to WT mice. The predisposition to allergic sensitization was associated with increased airway epithelial high-mobility group box 1 (HMGB1) translocation and release, increased type 2 innate lymphoid cells (ILC2s) and monocyte-derived dendritic cell numbers in the mediastinal lymph nodes, and increased T-helper type 2 (T_H 2)-cell differentiation. CD4⁺ T-cell depletion studies or the adoptive transfer of WT OVA-specific CD4⁺ T cells to WT or Pag1^-/- recipients demonstrated that the heightened propensity for T_H 2-cell differentiation was both T cell intrinsic and extrinsic.

CONCLUSION

PAG1 deficiency increased airway epithelial activation, ILC2 expansion, and T_H 2 differentiation. As a consequence, PAG1 deficiency predisposed toward allergic sensitization and increased the severity of experimental asthma.

Identification of cancer sex-disparity in the functional integrity of p53 and its X chromosome network

The disproportionately high prevalence of male cancer is poorly understood. We tested for sex-disparity in the functional integrity of the major tumor suppressor p53 in sporadic cancers. Our bioinformatics analyses expose three novel levels of p53 impact on sex-disparity in 12 non-reproductive cancer types. First, TP53 mutation is more frequent in these cancers among US males than females, with poorest survival correlating with its mutation. Second, numerous X-linked genes are associated with p53, including vital genomic regulators. Males are at unique risk from alterations of their single copies of these genes. High expression of X-linked negative regulators of p53 in wild-type TP53 cancers corresponds with reduced survival. Third, females exhibit an exceptional incidence of non-expressed mutations among p53-associated X-linked genes. Our data indicate that poor survival in males is contributed by high frequencies of TP53 mutations and an inability to shield against deregulated X-linked genes that engage in p53 networks.

Advances in asthma and allergic disease genetics: Is bigger always better?

This review focuses on genome-wide association studies (GWASs) of asthma and allergic diseases published between January 1, 2018, and June 30, 2019. During this time period, there were 38 GWASs reported in 19 articles, including the largest performed to date for many of these conditions. Overall, we learned that childhood-onset asthma is associated with the most independent loci compared with other defined groups of asthma and allergic disease cases; adult-onset asthma and moderate-to-severe asthma are associated with fewer genes, which are largely a subset of those associated with childhood-onset asthma. There is significant genetic overlap between asthma and allergic diseases, particularly with respect to childhood-onset asthma, which involves genes that reflect the importance of barrier function biology, and to HLA region genes, which are the most frequently associated genes overall in both groups of diseases. Although the largest GWASs in African American and Latino/Hispanic populations were reported during this period, they are still significantly underpowered compared with studies reported in populations of European ancestry, highlighting the need for larger studies, particularly in patients with childhood-onset asthma and allergic diseases, in these important populations that carry the greatest burden of disease.

Histone modifications and their role in epigenetics of atopy and allergic diseases

This review covers basic aspects of histone modification and the role of posttranslational histone modifications in the development of allergic diseases, including the immune mechanisms underlying this development. Together with DNA methylation, histone modifications (including histone acetylation, methylation, phosphorylation, ubiquitination, etc.) represent the classical epigenetic mechanisms. However, much less attention has been given to histone modifications than to DNA methylation in the context of allergy. A systematic review of the literature was undertaken to provide an unbiased and comprehensive update on the involvement of histone modifications in allergy and the mechanisms underlying this development. In addition to covering the growing interest in the contribution of histone modifications in regulating the development of allergic diseases, this review summarizes some of the evidence supporting this contribution. There are at least two levels at which the role of histone modifications is manifested. One is the regulation of cells that contribute to the allergic inflammation (T cells and macrophages) and those that participate in airway remodeling [(myo-) fibroblasts]. The other is the direct association between histone modifications and allergic phenotypes. Inhibitors of histone-modifying enzymes may potentially be used as anti-allergic drugs. Furthermore, epigenetic patterns may provide novel tools in the diagnosis of allergic disorders.

The Post-GWAS Era: From Association to Function

During the past 12 years, genome-wide association studies (GWASs) have uncovered thousands of genetic variants that influence risk for complex human traits and diseases. Yet functional studies aimed at delineating the causal genetic variants and biological mechanisms underlying the observed statistical associations with disease risk have lagged. In this review, we highlight key advances in the field of functional genomics that may facilitate the derivation of biological meaning post-GWAS. We highlight the evidence suggesting that causal variants underlying disease risk often function through regulatory effects on the expression of target genes and that these expression effects might be modest and cell-type specific. We moreover discuss specific studies as proof-of-principle examples for current statistical, bioinformatic, and empirical bench-based approaches to downstream elucidation of GWAS-identified disease risk loci.

Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks

We examined common variation in asthma risk by conducting a meta-analysis of worldwide asthma genome-wide association studies (23,948 cases, 118,538 controls) from ethnically-diverse populations. We identified five new asthma loci, uncovered two additional novel associations at two known asthma loci, established asthma associations at two loci implicated previously in comorbidity of asthma plus hay fever, and confirmed nine known loci. Investigation of pleiotropy showed large overlaps in genetic variants with autoimmune and inflammatory diseases. Enrichment of asthma risk loci in enhancer marks, especially in immune cells, suggests a major role of these loci in the regulation of immune-related mechanisms.

Lessons from ten years of genome-wide association studies of asthma

Twenty-five genome-wide association studies (GWAS) of asthma were published between 2007 and 2016, the largest with a sample size of 157242 individuals. Across these studies, 39 genetic variants in low linkage disequilibrium (LD) with each other were reported to associate with disease risk at a significance threshold of P<5 × 10⁻⁸, including 31 in populations of European ancestry. Results from analyses of the UK Biobank data (n=380 503) indicate that at least 28 of the 31 associations reported in Europeans represent true-positive findings, collectively explaining 2.5% of the variation in disease liability (median of 0.06% per variant). We identified 49 transcripts as likely target genes of the published asthma risk variants, mostly based on LD with expression quantitative trait loci (eQTL). Of these genes, 16 were previously implicated in disease pathophysiology by functional studies, including TSLP, TNFSF4, ADORA1, CHIT1 and USF1. In contrast, at present, there is limited or no functional evidence directly implicating the remaining 33 likely target genes in asthma pathophysiology. Some of these genes have a known function that is relevant to allergic disease, including F11R, CD247, PGAP3, AAGAB, CAMK4 and PEX14, and so could be prioritized for functional follow-up. We conclude by highlighting three areas of research that are essential to help translate GWAS findings into clinical research or practice, namely validation of target gene predictions, understanding target gene function and their role in disease pathophysiology and genomics-guided prioritization of targets for drug development.

Functional dissection of breast cancer risk-associated TERT promoter variants

The multi-cancer susceptibility locus at 5p15.33 includes TERT, encoding the telomerase catalytic subunit. Genome-wide association studies (GWAS) have identified six single nucleotide polymorphisms (SNPs) in the TERT promoter associated with decreased breast cancer risk, although the precise causal variants and their mechanisms of action have remained elusive. Luciferase reporter assays indicated that the protective haplotype reduced TERT promoter activity in human mammary epithelial and cancer cells in an estrogen-independent manner. Using single variant constructs, we identified rs3215401 and rs2853669 as likely functional variants. Silencing of MYC decreased TERT promoter activity but neither MYC nor ETS2 silencing conferred allele-specificity. In chromatin immunoprecipitation experiments, the ETS protein GABPA, but not ETS2 or ELF1, bound rs2853669 in an allele-specific manner in mammary epithelial cells. Investigation of open chromatin in mammoplasty samples suggested involvement of three additional variants, though not rs3215401 or rs2853669. Chromosome conformation capture revealed no interaction of the TERT promoter with regulatory elements in the locus, indicating limited local impact of candidate variants on the TERT promoter. Collectively, our functional studies of the TERT-CLPTM1L breast cancer susceptibility locus describe rs2853669 as a functional variant of this association signal among three other potentially causal variants and demonstrate the versatile mechanisms by which TERT promoter variants may affect breast cancer risk.

Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling

BACKGROUND

Hundreds of genetic variants are thought to contribute to variation in asthma risk by modulating gene expression. Methods that increase the power of genome-wide association studies (GWASs) to identify risk-associated variants are needed.

OBJECTIVE

We sought to develop a method that aggregates the evidence for association with disease risk across expression quantitative trait loci (eQTLs) of a gene and use this approach to identify asthma risk genes.

METHODS

We developed a gene-based test and software package called EUGENE that (1) is applicable to GWAS summary statistics; (2) considers both cis- and trans-eQTLs; (3) incorporates eQTLs identified in different tissues; and (4) uses simulations to account for multiple testing. We applied this approach to 2 published asthma GWASs (combined n = 46,044) and used mouse studies to provide initial functional insights into 2 genes with novel genetic associations.

RESULTS

We tested the association between asthma and 17,190 genes that were found to have cis- and/or trans-eQTLs across 16 published eQTL studies. At an empirical FDR of 5%, 48 genes were associated with asthma risk. Of these, for 37, the association was driven by eQTLs located in established risk loci for allergic disease, including 6 genes not previously implicated in disease cause (eg, LIMS1, TINF2, and SAFB). The remaining 11 significant genes represent potential novel genetic associations with asthma. The association with 4 of these replicated in an independent GWAS: B4GALT3, USMG5, P2RY13, and P2RY14, which are genes involved in nucleotide synthesis or nucleotide-dependent cell activation. In mouse studies, P2ry13 and P2ry14-purinergic receptors activated by adenosine 5-diphosphate and UDP-sugars, respectively-were upregulated after allergen challenge, notably in airway epithelial cells, eosinophils, and neutrophils. Intranasal exposure with receptor agonists induced the release of IL-33 and subsequent eosinophil infiltration into the lungs.

CONCLUSION

We identified novel associations between asthma and eQTLs for 4 genes related to nucleotide synthesis/signaling and demonstrated the power of gene-based analyses of GWASs.