Convergence of biomarkers and risk factor trait loci of coronary artery disease at 3p21.31 and HLA region

Role of the Neanderthal Genome in Genetic Susceptibility to COVID-19: 3p21.31 Locus in the Spotlight

Since the outbreak of COVID-19, genome-wide association studies have tried to discover the role of genetic predisposition in the clinical variability of this viral infection. The findings of various investigations have led to several loci for COVID-19 genetic susceptibility. Among candidate regions, the 3p21.31 locus has been in the spotlight among scientists, as it can increase the risk of severe COVID-19 by almost two fold. In addition to its substantial association with COVID-19 severity, this locus is related to some common diseases, such as diabetes, malignancies, and coronary artery disease. This locus also harbors evolutionary traces of Neanderthal genomes, which is believed to be the underlying reason for its association with COVID-19 severity. Additionally, the inheritance of this locus from Neanderthals seems to be under positive selection. This review aims to summarize a collection of evidence on the 3p21.31 locus and its impact on COVID-19 outcomes by focusing on the risk variants originated from the Neanderthal genome. Moreover, we discuss candidate genes at this locus and the possible mechanisms by which they influence the progression of COVID-19 symptoms. Better insights into human genetic susceptibility to newly emerging diseases such as COVID-19 and its evolutionary origin can provide fundamentals for risk assessment of different populations as well as the development of personalized prevention and treatments based on genomic medicine.

Inferring causal direction between two traits using R2 with application to transcriptome-wide association studies

In Mendelian randomization, two single SNP-trait correlation-based methods have been developed to infer the causal direction between an exposure (e.g., a gene) and an outcome (e.g., a trait), called MR Steiger's method and its recent extension called Causal Direction-Ratio (CD-Ratio). Here we propose an approach based on R2, the coefficient of determination, to combine information from multiple (possibly correlated) SNPs to simultaneously infer the presence and direction of a causal relationship between an exposure and an outcome. Our proposed method generalizes Steiger's method from using a single SNP to multiple SNPs as IVs. It is especially useful in transcriptome-wide association studies (TWASs) (and similar applications) with typically small sample sizes for gene expression (or another molecular trait) data, providing a more flexible and powerful approach to inferring causal directions. It can be applied to GWAS summary data with a reference panel. We also discuss the influence of invalid IVs and introduce a new approach called R2S to select and remove invalid IVs (if any) to enhance the robustness. We compared the performance of the proposed method with existing methods in simulations to demonstrate its advantages. We applied the methods to identify causal genes for high/low-density lipoprotein cholesterol (HDL/LDL) using the individual-level GTEx gene expression data and UK Biobank GWAS data. The proposed method was able to confirm some well-known causal genes while identifying some novel ones. Additionally, we illustrated an application of the proposed method to GWAS summary to infer causal relationships between HDL/LDL and stroke/coronary artery disease (CAD).

The association of rs17713054 with Neanderthal origin at 3p21.31 locus with the severity of COVID-19 in Iranian patients

Since the COVID-19 pandemic, the diversity of clinical manifestations in patients has been a tremendous challenge. It seems that genetic variations, as one of the players, contribute to the variety of symptoms. Genome-wide association studies have demonstrated the influence of certain genomic regions on the disease prognosis. Particularly, a haplotype at 3p21.31 locus, inherited from Neanderthals, showed an association with COVID-19 severity. Despite several studies regarding this haplotype, some key variants are not sufficiently addressed. In the present study, we investigated the association of rs17713054 at 3p21.31 with COVID-19 severity. We analyzed the genotype of 251 Iranian COVID-19 patients (151 patients with asymptomatic to mild form as control and 100 patients with severe to critical symptoms without any comorbidities as case group) using the ARMS-PCR method. Results demonstrated that the A allele confers an almost twofold increased risk for COVID-19 severity (P value = 0.008). The AA genotype also raises the risk by more than 11 times following the recessive model (P value = 0.013). In conclusion, the A allele in rs17713054 was a risk allele in Iranian patients and was independently associated with COVID-19 severity. More studies are beneficial to confirm these findings in other populations and to develop strategies for risk assessment, prevention, and personalized medicine.

Genome-wide search identified DNA methylation sites that regulate the metabolome

Background: Identifying DNA methylation sites that regulate the metabolome is important for several purposes. In this study, publicly available GWAS data were integrated to find methylation sites that impact metabolome through a discovery and replication scheme and by using Mendelian randomization. Results: The outcome of analyses revealed 107 methylation sites associated with 84 metabolites at the genome-wide significance level (p<5e^-8) at both the discovery and replication stages. A large percentage of the observed associations (85%) were with lipids, significantly higher than expected (p = 0.0003). A number of CpG (methylation) sites showed specificity e.g., cg20133200 within PFKP was associated with glucose only and cg10760299 within GATM impacted the level of creatinine; in contrast, there were sites associated with numerous metabolites e.g., cg20102877 on the 2p23.3 region was associated with 39 metabolites. Integrating transcriptome data enabled identifying genes (N = 82) mediating the impact of methylation sites on the metabolome and cardiometabolic traits. For example, PABPC4 mediated the impact of cg15123755-HDL on type-2 diabetes. KCNK7 mediated the impact of cg21033440-lipids on hypertension. POC5, ILRUN, FDFT1, and NEIL2 mediated the impact of CpG sites on obesity through metabolic pathways. Conclusion: This study provides a catalog of DNA methylation sites that regulate the metabolome for downstream applications.

Identification of a novel role for the immunomodulator ILRUN in the development of several T cell subsets in mice

Inflammation and lipid regulator with UBA-like and NBR1-like domains (ILRUN) is a protein-encoding gene associated with innate immune signaling, lipid metabolism and cancer. In the context of innate immunity, ILRUN inhibits IRF3-mediated transcription of antimicrobial and proinflammatory cytokines by inducing degradation of the transcriptional coactivators CBP and p300. There remains a paucity of information, however, regarding the innate immune roles of ILRUN beyond in vitro analyses. To address this, we utilize a knockout mouse model to investigate the effect of ILRUN on cytokine expression in splenocytes and on the development of immune cell populations in the spleen and thymus. We show elevated production of tumor necrosis factor and interleukin-6 cytokines in ILRUN-deficient splenocytes following stimulation with the innate immune ligands polyinosinic:polycytidylic acid or lipopolysaccharide. Differences were also observed in the populations of several T cell subsets, including regulatory, mucosal-associated invariant and natural killer. These data identify novel functions for ILRUN in the development of certain immune cell populations and support previous in vitro findings that ILRUN negatively regulates the synthesis of pathogen-stimulated cytokines. This establishes the ILRUN knockout mouse model as a valuable resource for further study of the functions of ILRUN in health and disease.

Mediation roles of neutrophils and high-density lipoprotein (HDL) on the relationship between HLA-DQB1 and rosacea

BACKGROUND

Though the previous genome-wide association studies found the association between HLA alleles and rosacea in the European populations, the data is lacking among the Asians. Moreover, neutrophils are important in the immune-related mechanism of rosacea, and dyslipidemia is closely related to rosacea. We aimed to explore the association between HLA genes and rosacea in Chinese rosacea patients, as well as the mediation effect of neutrophils, high-density lipoprotein (HDL) and low-density lipoprotein (LDL) on the relationship between HLA genes and rosacea.

METHODS

A total of 249 rosacea and 150 controls were ranked by the international investigator global rosacea severity scores. HLA genes, neutrophils, HDL, and LDL were detected. And their mediation effects on the relationship between HLA and rosacea risk or severity were analysed.

RESULTS

HLA-DQB1*03:03 allele (OR = 41.89, 95% CI: 9.80 ∼ 179.09, p = 4.7*10-7), HLA-DQB1*04:02 allele (OR = 0.16, 95% CI: 0.03 ∼ 0.81, p = 0.026) and HLA-DQB1*03:03/05:02 genotype (OR = 5.57, 95% CI: 1.13 ∼ 27.52, p = 0.0351) were significantly associated with rosacea. Moreover, HLA-DQB1*03:03 allele (b = 1.434, SE = 0.217, p = 2.0*10-10), HLA-DQB1*05:01 allele (b = 0.894, SE = 0.33520, p = 0.008) and HLA-DQB1*03:03/06:01 genotype (b = 0.998, SE = 0.472, p = 0.040) were positively associated with rosacea severity. Furthermore, we found both neutrophils and HDL, instead of LDL, have mediation effects on the relationship between HLA-DQB1*03:03 and risk or severity of rosacea.

CONCLUSIONS

We discovered novel susceptible HLA alleles for rosacea in the Chinese population, and disclosed the mediation effect of neutrophils and HDL on the relationship between HLA-DQB1 and rosacea, implying a possible correlation between rosacea and inflammatory or metabolic factors, providing hints for future studies in the mechanism of rosacea. Key messagesHLA-DQB1*03:03 allele, HLA-DQB1*04:02 allele and HLA-DQB1*03:03/05:02 genotype were significantly associated with rosacea.HLA-DQB1*03:03 allele, HLA-DQB1*05:01 allele and HLA-DQB1*03:03/06:01 genotype were positively associated with rosacea severity.Neutrophils and HDL have mediation effects on the relationship between HLA-DQB1*03:03 and risk or severity of rosacea.

Epigenome-Wide Study Identified Methylation Sites Associated with the Risk of Obesity

Here, we performed a genome-wide search for methylation sites that contribute to the risk of obesity. We integrated methylation quantitative trait locus (mQTL) data with BMI GWAS information through a SNP-based multiomics approach to identify genomic regions where mQTLs for a methylation site co-localize with obesity risk SNPs. We then tested whether the identified site contributed to BMI through Mendelian randomization. We identified multiple methylation sites causally contributing to the risk of obesity. We validated these findings through a replication stage. By integrating expression quantitative trait locus (eQTL) data, we noted that lower methylation at cg21178254 site upstream of CCNL1 contributes to obesity by increasing the expression of this gene. Higher methylation at cg02814054 increases the risk of obesity by lowering the expression of MAST3, whereas lower methylation at cg06028605 contributes to obesity by decreasing the expression of SLC5A11. Finally, we noted that rare variants within 2p23.3 impact obesity by making the cg01884057 site more susceptible to methylation, which consequently lowers the expression of POMC, ADCY3 and DNAJC27. In this study, we identify methylation sites associated with the risk of obesity and reveal the mechanism whereby a number of these sites exert their effects. This study provides a framework to perform an omics-wide association study for a phenotype and to understand the mechanism whereby a rare variant causes a disease.