Regulatory Noncoding and Predicted Pathogenic Coding Variants of CCR5 Predispose to Severe COVID-19

Genome-wide association studies (GWAS) found locus 3p21.31 associated with severe COVID-19. CCR5 resides at the same locus and, given its known biological role in other infection diseases, we investigated if common noncoding and rare coding variants, affecting CCR5, can predispose to severe COVID-19. We combined single nucleotide polymorphisms (SNPs) that met the suggestive significance level (P ≤ 1 × 10⁻⁵) at the 3p21.31 locus in public GWAS datasets (6406 COVID-19 hospitalized patients and 902,088 controls) with gene expression data from 208 lung tissues, Hi-C, and Chip-seq data. Through whole exome sequencing (WES), we explored rare coding variants in 147 severe COVID-19 patients. We identified three SNPs (rs9845542, rs12639314, and rs35951367) associated with severe COVID-19 whose risk alleles correlated with low CCR5 expression in lung tissues. The rs35951367 resided in a CTFC binding site that interacts with CCR5 gene in lung tissues and was confirmed to be associated with severe COVID-19 in two independent datasets. We also identified a rare coding variant (rs34418657) associated with the risk of developing severe COVID-19. Our results suggest a biological role of CCR5 in the progression of COVID-19 as common and rare genetic variants can increase the risk of developing severe COVID-19 by affecting the functions of CCR5.     

Genetic Predisposition to the Mortality in Septic Shock Patients: From GWAS to the Identification of a Regulatory Variant Modulating the Activity of a CISH Enhancer

The high mortality rate in septic shock patients is likely due to environmental and genetic factors, which influence the host response to infection. Two genome-wide association studies (GWAS) on 832 septic shock patients were performed. We used integrative bioinformatic approaches to annotate and prioritize the sepsis-associated single nucleotide polymorphisms (SNPs). An association of 139 SNPs with death based on a false discovery rate of 5% was detected. The most significant SNPs were within the CISH gene involved in cytokine regulation. Among the 139 SNPs associated with death and the 1311 SNPs in strong linkage disequilibrium with them, we investigated 1439 SNPs within non-coding regions to identify regulatory variants. The highest integrative weighted score (IW-score) was obtained for rs143356980, indicating that this SNP is a robust regulatory candidate. The rs143356980 region is located in a non-coding region close to the CISH gene. A CRISPR-Cas9-mediated deletion of this region and specific luciferase assays in K562 cells showed that rs143356980 modulates the enhancer activity in K562 cells. These analyses allowed us to identify several genes associated with death in patients with septic shock. They suggest that genetic variations in key genes, such as CISH, perturb relevant pathways, increasing the risk of death in sepsis patients.     

Influence of Insulin Receptor Single Nucleotide Polymorphisms on Glycaemic Control and Formation of Anti-Insulin Antibodies in Diabetes Mellitus

Single nucleotide polymorphisms (SNPs) in insulin and insulin receptor genes may influence the interaction between the two molecules, as may anti-insulin antibodies (IAs), commonly found in patients with type 1 diabetes mellitus (T1D) or type 2 diabetes mellitus (T2D) treated with exogenous insulin. We examined the impact of two SNPs in the human insulin gene (INS), rs3842752 and rs689, and two in the insulin receptor gene (INSR) rs2245649 and rs2229429, on disease susceptibility, glycaemic control, and IAs formation in 100 T1D patients and 101 T2D patients treated with insulin. 79 individuals without diabetes were typed as healthy controls. The minor alleles of rs3842752 and rs689 in INS protected against T1D (OR: 0.50, p = 0.01 and OR: 0.44; p = 0.002, respectively). The minor alleles of both rs2245649 and rs2229429 in INSR were risk factors for poor glycaemic control (HbA1c ≥ 80 mmol/mol) in T1D (OR: 5.35, p = 0.009 and OR: 3.10, p = 0.01, respectively). Surprisingly, the minor alleles of rs2245649 and rs2229429 in INSR associated strongly with the absence of IAs in T1D (OR = 0.28, p = 0.008 and OR = 0.30, p = 0.002, respectively). In conclusion, the minor alleles of the investigated INS SNPs protect against T1D, and the minor alleles of the investigated INSR SNPs are associated with poor glycaemic control and the absence of IAs in T1D.     

The Analysis of a Genome-Wide Association Study (GWAS) of Overweight and Obesity in Psoriasis

There is evidence that the concomitance of psoriasis and obesity may originate from the interplay between multiple genetic pathways and involve gene–gene interactions. The aim of this study was to compare the genetic background related to obesity among psoriatic patients versus healthy controls by means of a Genome-Wide Association Study (GWAS). A total of 972 psoriatic patients and a total of 5878 healthy donors were enrolled in this study. DNA samples were genotyped for over 500,000 single nucleotide polymorphisms (SNPs) using Infinium CoreExome BeadChips (Illumina, San Diego, CA, USA). Statistical analysis identified eleven signals (p < 1 × 10⁻⁵) associated with BMI across the study groups and revealed a varying effect size in each sub-cohort. Seven of the alternative alleles (rs1558902 in the FTO gene, rs696574 in the CALCRL gene, as well as rs10968110, rs4551082, rs4609724, rs9320269, and rs2338833,) are associated with increased BMI among all psoriatic patients and four (rs1556519 in the ITLN2 gene, rs12972098 in the AC003006.7 gene, rs12676670 in the PAG1 gene, and rs1321529) are associated with lower BMI. The results of our study may lead to further insights into the understanding of the pathogenesis of obesity among psoriatic patients.     

The Role of Nonshivering Thermogenesis Genes on Leptin Levels Regulation in Residents of the Coldest Region of Siberia

Leptin plays an important role in thermoregulation and is possibly associated with the microevolutionary processes of human adaptation to a cold climate. In this study, based on the Yakut population (n = 281 individuals) living in the coldest region of Siberia (t°minimum −71.2 °C), we analyze the serum leptin levels and data of 14 single nucleotide polymorphisms (SNPs) of 10 genes (UCP1, UCP2, UCP3, FNDC5, PPARGC1A, CIDEA, PTGS2, TRPV1, LEPR, BDNF) that are possibly involved in nonshivering thermogenesis processes. Our results demonstrate that from 14 studied SNPs of 10 genes, 2 SNPs (the TT rs3811787 genotype of the UCP1 gene and the GG rs6265 genotype of the BDNF gene) were associated with the elevated leptin levels in Yakut females (p < 0.05). Furthermore, of these two SNPs, the rs3811787 of the UCP1 gene demonstrated more indications of natural selection for cold climate adaptation. The prevalence gradient of the T-allele (rs3811787) of UCP1 increased from the south to the north across Eurasia, along the shore of the Arctic Ocean. Thereby, our study suggests the potential involvement of the UCP1 gene in the leptin-mediated thermoregulation mechanism, while the distribution of its allelic variants is probably related to human adaptation to a cold climate.     

Implication of Melanocortin Receptor Genes in the Familial Comorbidity of Type 2 Diabetes and Depression

The melanocortin receptors are G-protein-coupled receptors, which are essential components of the hypothalamic–pituitary–adrenal axis, and they mediate the actions of melanocortins (melanocyte-stimulating hormones: α-MSH, β-MSH, and γ-MSH) as well as the adrenocorticotropin hormone (ACTH) in skin pigmentation, adrenal steroidogenesis, and stress response. Three melanocortin receptor genes (MC1R, MC2R, and MC5R) contribute to the risk of major depressive disorder (MDD), and one melanocortin receptor gene (MC4R) contributes to the risk of type 2 diabetes (T2D). MDD increases T2D risk in drug-naïve patients; thus, MDD and T2D commonly coexist. The five melanocortin receptor genes might confer risk for both disorders. However, they have never been investigated jointly to evaluate their potential contributing roles in the MDD-T2D comorbidity, specifically within families. In 212 Italian families with T2D and MDD, we tested 11 single nucleotide polymorphisms (SNPs) in the MC1R gene, 9 SNPs in MC2R, 3 SNPs in MC3R, 4 SNPs in MC4R, and 2 SNPs in MC5R. The testing used 2-point parametric linkage and linkage disequilibrium (LD) (i.e., association) analysis with four models (dominant with complete penetrance (D1), dominant with incomplete penetrance (D2), recessive with complete penetrance (R1), and recessive with incomplete penetrance (R2)). We detected significant (p ≤ 0.05) linkage and/or LD (i.e., association) to/with MDD for one SNP in MC2R (rs111734014) and one SNP in MC5R (rs2236700), and to/with T2D for three SNPs in MC1R (rs1805007 and rs201192930, and rs2228479), one SNP in MC2R (rs104894660), two SNPs in MC3R (rs3746619 and rs3827103), and one SNP in MC4R genes (Chr18-60372302). The linkage/LD/association was significant across different linkage patterns and different modes of inheritance. All reported variants are novel in MDD and T2D. This is the first study to report risk variants in MC1R, MC2R, and MC3R genes in T2D. MC2R and MC5R genes are replicated in MDD, with one novel variant each. Within our dataset, only the MC2R gene appears to confer risk for both MDD and T2D, albeit with different risk variants. To further clarity the role of the melanocortin receptor genes in MDD-T2D, these findings should be sought among other ethnicities as well.     

Matrix Metalloproteinase Gene Polymorphisms Are Associated with Breast Cancer in the Caucasian Women of Russia

We conducted this study to explore the association between matrix metalloproteinase (MMP) gene polymorphisms and breast cancer (BC) risk in the Caucasian women of Russia. In total, 358 affected (BC) and 746 unaffected (cancer-free) women were included in this case-control retrospective study. From BC-related genes in previous studies, ten single nucleotide polymorphisms (SNPs) in five MMP genes (MMP1, 2, 3, 8, 9) were genotyped. The BC risk was calculated by logistic regression (to evaluate the SNPs’ independent effects) and model-based multifactor dimensionality reduction (MB-MDR) (to identify SNP–SNP interactions) methods. The allelic variants’ distribution of c.836 A > G (rs17576) and c. 1721 C > G (rs2250889) MMP9 was significantly different between BC and cancer-free women: for G minor alleles, these SNPs manifested disorder protective effects (OR 0.82 and OR 0.67–0.71, respectively, p_perm ≤ 0.035). Eleven haplotypes of six SNPs MMP9 were involved in BC risk (nine haplotypes) and protective (two haplotypes) effects. All 10 SNPs of the MMP genes examined were associated with BC within the 13 SNP–SNP interaction simulated models, with a pivotal role of the two-locus (rs17577 × rs3918242) MMP9 epistatic interaction (defined as 1.81% BC entropy within more than 60% of the genetic models). Under in silico bioinformatics, BC susceptibility MMP polymorphic loci are located in functionally active genome regions and impact genes expression and splicing “regulators” in the mammary gland. The biological pathways of BC MMP candidate genes are mainly realized due to metalloendopeptidase activity and extracellular matrix organization (structure, disassembly, metabolic process, etc.). In conclusion, our data show that MMP gene polymorphisms are related to BC susceptibility in the Caucasian women of Russia.     

Variants of the Low Oxygen Sensors EGLN1 and HIF-1AN Associated with Acute Mountain Sickness

Two low oxygen sensors, Egl nine homolog 1 (EGLN1) and hypoxia-inducible factor 1-α inhibitor (HIF-1AN), play pivotal roles in the regulation of HIF-1α, and high altitude adaption may be involved in the pathology of acute mountain sickness (AMS). Here, we aimed to analyze single nucleotide polymorphisms (SNPs) in the untranslated regions of the EGLN1 and HIF-1AN genes and SNPs chosen from a genome-wide adaptation study of the Han Chinese population. To assess the association between EGLN1 and HIF-1AN SNPs and AMS in a Han Chinese population, a case–control study was performed including 190 patients and 190 controls. In total, thirteen SNPs were genotyped using the MassARRAY^® MALDI-TOF system. Multiple genetic models were tested; The Akaike’s information criterion (AIC) and Bayesian information criterion (BIC) values indicated that the dominant model may serve as the best-fit model for rs12406290 and rs2153364 of significant difference. However, these data were not significant after Bonferroni correction. No significant association was noted between AMS and rs12757362, rs1339894, rs1361384, rs2009873, rs2739513 or rs2486729 before and after Bonferroni correction. Further haplotype analyses indicated the presence of two blocks in EGLN1; one block consists of rs12406290-rs2153364, located upstream of the EGLN1 gene. Carriers of the “GG” haplotype of rs12406290-rs2153364 exhibited an increased risk of AMS after adjustments for age and smoking status. However, no significant association was observed among HIF-1AN 3''-untranslated region (3''-UTR) polymorphisms, haplotype and AMS. Our study indicates that variants in the EGLN1 5''-UTR influence the susceptibility to AMS in a Han Chinese population.     

Association of LTA and SOD Gene Polymorphisms with Cerebral White Matter Hyperintensities in Migraine Patients

White matter hyperintensities (WMHs) in migraine could be related to inflammatory and antioxidant events. The aim of this study is to verify whether migraine patients with WMHs carry a genetic pro-inflammatory/pro-oxidative status. To test this hypothesis, we analyzed lymphotoxin alpha (LTA; rs2071590T and rs2844482G) and superoxide dismutase 1 (SOD1; rs2234694C) and 2 (SOD2; rs4880T) gene polymorphisms (SNPs) in 370 consecutive patients affected by episodic (EM; n = 251) and chronic (CM; n = 119) migraine and in unrelated healthy controls (n = 100). Brain magnetic resonance was available in 183/370 patients. The results obtained show that genotypes and allele frequencies for all tested SNPs did not differ between patients and controls. No association was found between single SNPs or haplotypes and sex, migraine type, cardiovascular risk factors or disorders. Conversely, the LTA rs2071590T (OR = 2.2) and the SOD1 rs2234694C (OR = 4.9) alleles were both associated with WMHs. A four-loci haplotype (TGCT haplotype: rs2071590T/rs2844482G/rs2234694C/rs4880T) was significantly more frequent in migraineurs with WMHs (7 of 38) compared to those without WMHs (4 of 134; OR = 8.7). We may, therefore, conclude by suggesting that that an imbalance between pro-inflammatory/pro-oxidative and antioxidant events in genetically predisposed individuals may influence the development of WMHs.     

Interleukin-6 Receptor rs7529229 T/C Polymorphism Is Associated with Left Main Coronary Artery Disease Phenotype in a Chinese Population

Left main coronary artery disease (LMCAD) is a particular severe phenotype of coronary artery disease (CAD) and heritability. Interleukin (IL) may play important roles in the pathogenesis of CAD. Although several single nucleotide polymorphisms (SNPs) identified in IL related genes have been evaluated for their roles in inflammatory diseases and CAD predisposition, the investigations between genetic variants and CAD phenotype are limited. We hypothesized that some of these gene SNPs may contribute to LMCAD phenotype susceptibility compared with more peripheral coronary artery disease (MPCAD). In a hospital-based case-only study, we studied IL-1A rs1800587 C/T, IL-1B rs16944 G/A, IL-6 rs1800796 C/G, IL-6R rs7529229 T/C, IL-8 rs4073 T/A, IL-10 rs1800872 A/C, and IL-10 rs1800896 A/G SNPs in 402 LMCAD patients and 804 MPCAD patients in a Chinese population. Genotyping was done using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) and ligation detection reaction (LDR) method. When the IL-6R rs7529229 TT homozygote genotype was used as the reference group, the CC or TC/CC genotypes were associated with the increased risk for LMCAD (CC vs. TT, adjusted odds ratio(OR) = 1.46, 95% confidence interval (CI) = 1.02–2.11, p = 0.042; CC + TC vs. TT, adjusted OR = 1.31, 95% CI = 1.02–1.69, p = 0.037). None of the other six SNPs achieved any significant differences between LMCAD and MPCAD. The present study suggests that IL-6R rs7529229 T/C functional SNP may contribute to the risk of LMCAD in a Chinese population. However, our results were limited. Validation by a larger study from a more diverse ethnic population is needed.     

Peptidylarginine Deiminase 2 Gene Polymorphisms in Subjects with Periodontitis Predispose to Rheumatoid Arthritis

Epidemiologic studies have shown associations between periodontitis and rheumatoid arthritis (RA), but a causal relationship has not been established. Citrullination of gingival proteins by human peptidylarginine deiminases (PADs) or PAD from Porphyromonas gingivalis has been proposed to generate autoantigens in anti-CCP-positive RA. This study investigated whether the association between periodontitis and RA is influenced by single nucleotide polymorphisms (SNPs) in the genes encoding PAD2 and PAD4 that catalyze aberrant citrullination in RA and often are overexpressed in inflamed gingival connective tissue in subjects with periodontitis. The study included 137 RA patients and 161 controls with self-reported periodontitis. Periodontitis onset preceded RA onset by 13 years on average and was not associated with any of the SNPs investigated. In subjects with periodontitis, carriage of the minor alleles of rs2057094 and rs2235912 in PADI2 significantly increased the risk of RA (odds ratios 1.42 [p = 0.03] and 1.48 [p = 0.02], respectively), and this effect was driven by the anti-CCP-negative RA patients. The minor alleles of these SNPs only increased risk of anti-CCP-positive RA in individuals with periodontitis and a history of smoking. These data suggest that individuals with periodontitis carrying the minor alleles of SNPs rs2057094, rs2076616 and rs2235912 in PADI2 may be at increased risk of RA.     

A Polymorphism at the 3′-UTR Region of the Aromatase Gene Is Associated with the Efficacy of the Aromatase Inhibitor,Anastrozole, in Metastatic Breast Carcinoma

Estrogen-related genes and the fat mass and obesity-associated (FTO) gene play a critical role in estrogen metabolism, and those polymorphisms are associated with a poor prognosis in breast cancer. However, little is known about the association between these polymorphisms and the efficacy of anastrozole. The aim was to investigate the impact of the genetic polymorphisms, CYP19A1, 17-β-HSD-1 and FTO, on the response to anastrozole in metastatic breast carcinoma (MBC) and to evaluate the impact of those polymorphisms on various clinicopathologic features. Two-hundred seventy-two women with hormone receptor-positive MBC treated with anastrozole were identified retrospectively. DNA was extracted from peripheral blood and genotyped for five variants in three candidate genes. Time to progression was improved in patients carrying the variant alleles of rs4646 when compared to patients with the wild-type allele (16.40 months versus 13.52 months; p = 0.049). The rs4646 variant alleles were significantly associated with longer overall survival (37.3 months versus 31.6 months; p = 0.007). This relationship was not observed with the rs10046, rs2830, rs9926298 and rs9939609 polymorphisms. The findings of this study indicate that rs4646 polymorphism in the CYP19A1 gene may serve as a prognostic maker of the response to anastrozole in patients with MBC who are treated with anastrozole.     

The Influence of PSCA Gene Variation on Its Expression and Gastric Adenocarcinoma Susceptibility in the Northwest Chinese Population

Gastric adenocarcinoma (GAC) imposes a considerable health burden around the world. Gene variation in prostate stem cell antigen gene (PSCA) has been identified to be associated with GAC risk, while the results showed regional variation. To explore the influence of PSCA gene variation on its expression and GAC risk in the Northwest Chinese population, four single nucleotide polymorphisms (SNPs) of PSCA were genotyped in 476 GAC cases and 481 controls using MassARRAY system. Two SNPs of rs2294008 (C>T) and rs2976392 (G>A) were identified to be associated with GAC risk. rs2294008, rs2976392 and rs10216533 made up two statistically significant haplotypes (Hap-CGG and Hap-TAG). Additionally, PSCA expression was analyzed by quantitative real time PCR, immunohistochemistry and tissue microarray. The results showed that PSCA expression was decreased in GAC tissues compared with adjacent normal tissues. For normal tissues, PSCA expression was higher with Hap-TA than that with Hap-CG. For GAC tissues, the differentiation degree of Hap-TA was higher than that of Hap-CG. The expression distribution of PSCA in multiple human organs showed disparity. These results suggest that PSCA gene variation has a potential effect on its expression and GAC risk in the Northwest Chinese population.     

Alzheimer’s Disease-Associated SNP rs708727 in SLC41A1 May Increase Risk for Parkinson’s Disease: Report from Enlarged Slovak Study

SLC41A1 (A1) SNPs rs11240569 and rs823156 are associated with altered risk for Parkinson’s disease (PD), predominantly in Asian populations, and rs708727 has been linked to Alzheimer’s disease (AD). In this study, we have examined a potential association of the three aforementioned SNPs and of rs9438393, rs56152218, and rs61822602 (all three lying in the A1 promoter region) with PD in the Slovak population. Out of the six tested SNPs, we have identified only rs708727 as being associated with an increased risk for PD onset in Slovaks. The minor allele (A) in rs708727 is associated with PD in dominant and completely over-dominant genetic models (OR_D = 1.36 (1.05–1.77), p = 0.02, and OR_COD = 1.34 (1.04–1.72), p = 0.02). Furthermore, the genotypic triplet GG_(rs708727) + AG_(rs823156) + CC_(rs61822602) might be clinically relevant despite showing a medium (h ≥ 0.5) size difference (h = 0.522) between the PD and the control populations. RandomForest modeling has identified the power of the tested SNPs for discriminating between PD-patients and the controls to be essentially zero. The identified association of rs708727 with PD in the Slovak population leads us to hypothesize that this A1 polymorphism, which is involved in the epigenetic regulation of the expression of the AD-linked gene PM20D1, is also involved in the pathoetiology of PD (or universally in neurodegeneration) through the same or similar mechanism as in AD.     

Two Single Nucleotide Polymorphisms in the Purinergic Receptor P2X7 Gene Are Associated with Disease Severity in Multiple Sclerosis

Multiple sclerosis (MS) is a chronic autoimmune inflammatory disease of the central nervous system (CNS) that leads to progressive physical disability. Recent evidence has suggested that P2X7 receptor (P2X7R)-mediated purinergic signalling pathways play a role in MS-associated neuroinflammation, possibly contributing to disease pathogenesis. To evaluate possible associations between P2X7R polymorphisms and MS disease severity, we performed an association study of five non-synonymous SNPs coding variants of the P2X7R gene: rs1718119 Ala348Thr, rs2230911 Thr357Ser, rs2230912 Gln460Arg, rs3751143 Glu496Ala, and rs28360457 Arg307Gln, modulating P2X7R expression in 128 MS patients (relapsing remitting MS, RRMS: n = 94; secondary progressive, SPMS: n = 34). All patients were genotyped, and multiple sclerosis severity score (MSSS) was evaluated in every case; 189 healthy subjects were enrolled as well as controls. Results showed that P2X7R rs1718119(A) 348Thr and rs22390912(G) 464Arg, two SNPs of minor allele frequency (MAF) known to confer gain of function to the P2X7R protein, were associated with significantly higher MSSS in RRMS patients alone (SMRR (p < 0.001, p = 0.01, respectively)). Interestingly, two whole haplotypes resulted in having significant association with MSSS in these same patients. Thus: (1) the P2X7R-4 “ACGAG” haplotype, characterized by the co-presence of the rs1718119-rs2230912 AG MAF alleles, was associated with higher MSSS (Beta: 1.11 p = 0.04), and (2) the P2X7R-1 “GCAAG” complementary haplotype, which contains the rs1718119 and rs2230912 GA wild-type alleles, was more frequently carried by patients with lower MSSS and less severe disease (Beta: −1.54 p < 0.001). Although being preliminary and needing confirmation in an ampler cohort, these results suggest that 348Thr and 464Arg variants have a role as modulators of disease severity in RRMS patients.     

Association of Vitamin D Metabolism Gene Polymorphisms with Autoimmunity: Evidence in Population Genetic Studies

A high prevalence of vitamin D (calcidiol) serum deficiency has been described in several autoimmune diseases, including multiple sclerosis (MS), rheumatoid arthritis (AR), and systemic lupus erythematosus (SLE). Vitamin D is a potent immunonutrient that through its main metabolite calcitriol, regulates the immunomodulation of macrophages, dendritic cells, T and B lymphocytes, which express the vitamin D receptor (VDR), and they produce and respond to calcitriol. Genetic association studies have shown that up to 65% of vitamin D serum variance may be explained due to genetic background. The 90% of genetic variability takes place in the form of single nucleotide polymorphisms (SNPs), and SNPs in genes related to vitamin D metabolism have been linked to influence the calcidiol serum levels, such as in the vitamin D binding protein (VDBP; rs2282679 GC), 25-hydroxylase (rs10751657 CYP2R1), 1α-hydroxylase (rs10877012, CYP27B1) and the vitamin D receptor (FokI (rs2228570), BsmI (rs1544410), ApaI (rs7975232), and TaqI (rs731236) VDR). Therefore, the aim of this comprehensive literature review was to discuss the current findings of functional SNPs in GC, CYP2R1, CYP27B1, and VDR associated to genetic risk, and the most common clinical features of MS, RA, and SLE.