MEDTEC Students against Coronavirus: Investigating the Role of Hemostatic Genes in the Predisposition to COVID-19 Severity

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of the coronavirus disease 2019 (COVID-19) pandemic. Besides virus intrinsic characteristics, the host genetic makeup is predicted to account for the extreme clinical heterogeneity of the disease, which is characterized, among other manifestations, by a derangement of hemostasis associated with thromboembolic events. To date, large-scale studies confirmed that genetic predisposition plays a role in COVID-19 severity, pinpointing several susceptibility genes, often characterized by immunologic functions. With these premises, we performed an association study of common variants in 32 hemostatic genes with COVID-19 severity. We investigated 49,845 single-nucleotide polymorphism in a cohort of 332 Italian severe COVID-19 patients and 1668 controls from the general population. The study was conducted engaging a class of students attending the second year of the MEDTEC school (a six-year program, held in collaboration between Humanitas University and the Politecnico of Milan, allowing students to gain an MD in Medicine and a Bachelor’s Degree in Biomedical Engineering). Thanks to their willingness to participate in the fight against the pandemic, we evidenced several suggestive hits (p < 0.001), involving the PROC, MTHFR, MTR, ADAMTS13, and THBS2 genes (top signal in PROC: chr2:127192625:G:A, OR = 2.23, 95%CI = 1.50–3.34, p = 8.77 × 10−5). The top signals in PROC, MTHFR, MTR, ADAMTS13 were instrumental for the construction of a polygenic risk score, whose distribution was significantly different between cases and controls (p = 1.62 × 10−8 for difference in median levels). Finally, a meta-analysis performed using data from the Regeneron database confirmed the contribution of the MTHFR variant chr1:11753033:G:A to the predisposition to severe COVID-19 (pooled OR = 1.21, 95%CI = 1.09–1.33, p = 4.34 × 10−14 in the weighted analysis).

Comparison of DNA Methylation Profiles of Hemostatic Genes between Liver Tissue and Peripheral Blood within Individuals

DNA methylation has become increasingly recognized in the etiology of complex diseases, including thrombotic disorders. Blood is often collected in epidemiological studies for genotyping and has recently also been used to examine DNA methylation in epigenome-wide association studies. DNA methylation patterns are often tissue-specific, thus, peripheral blood may not accurately reflect the methylation pattern in the tissue of relevance. Here, we collected paired liver and blood samples concurrently from 27 individuals undergoing liver surgery. We performed targeted bisulfite sequencing for a set of 35 hemostatic genes primarily expressed in liver to analyze DNA methylation levels of >10,000 cytosine-phosphate-guanine (CpG) dinucleotides. We evaluated whether DNA methylation in blood could serve as a proxy for DNA methylation in liver at individual CpGs. Approximately 30% of CpGs were nonvariable and were predominantly hypo- (<25%) or hypermethylated (>70%) in both tissues. While blood can serve as a proxy for liver at these CpGs, the low variability renders these unlikely to explain phenotypic differences. We therefore focused on CpG sites with variable methylation levels in liver. The level of blood–liver tissue correlation varied widely across these variable CpGs; moderate correlations (0.5 ≤ r < 0.75) were detected for 6% and strong correlations (r ≥ 0.75) for a further 4%. Our findings indicate that it is essential to study the concordance of DNA methylation between blood and liver at individual CpGs. This paired blood–liver dataset is intended as a resource to aid interpretation of blood-based DNA methylation results.

A Comprehensive Sequencing-Based Analysis of Allelic Methylation Patterns in Hemostatic Genes in Human Liver

Characterizing the relationship between genetic, epigenetic (e.g., deoxyribonucleic acid [DNA] methylation), and transcript variation could provide insights into mechanisms regulating hemostasis and potentially identify new drug targets. Several hemostatic factors are synthesized in the liver, yet high-resolution DNA methylation data from human liver tissue is currently lacking for these genes. Single-nucleotide polymorphisms (SNPs) can influence DNA methylation in cis which can affect gene expression. This can be analyzed through allele-specific methylation (ASM) experiments. We performed targeted genomic DNA- and bisulfite-sequencing of 35 hemostatic genes in human liver samples for SNP and DNA methylation analysis, respectively, and integrated the data for ASM determination. ASM-associated SNPs (ASM-SNPs) were tested for association to gene expression in liver using in-house generated ribonucleic acid-sequencing data. We then assessed whether ASM-SNPs associated with gene expression, plasma proteins, or other traits relevant for hemostasis using publicly available data. We identified 112 candidate ASM-SNPs. Of these, 68% were associated with expression of their respective genes in human liver or in other human tissues and 54% were associated with the respective plasma protein levels, activity, or other relevant hemostatic genome-wide association study traits such as venous thromboembolism, coronary artery disease, stroke, and warfarin dose maintenance. Our study provides the first detailed map of the DNA methylation landscape and ASM analysis of hemostatic genes in human liver tissue, and suggests that methylation regulated by genetic variants in cis may provide a mechanistic link between noncoding SNPs and variation observed in circulating hemostatic proteins, prothrombotic diseases, and drug response.

Prothrombogenic polymorphic variants of hemostatic and folate metabolism genes In patients with aseptic cerebral venous thrombosis

Cerebral venous sinus thrombosis (CVT) becomes the cause of stroke in less than 1% of patients. In 20-30% of patients, the cause of thrombosis remains unclear, and thrombosis is considered idiopathic. Inherited hypercoagulable conditions significantly increase the risk of CVT. The aim of the study was to evaluate the frequency of prothrombogenic polymorphic variants of hemostatic and methionine-homocysteine metabolism genes alleles and genotypes in patients with aseptic CVT. Fifty one patients aged 18–75 with aseptic CVT were examined. The control group included 36 healthy volunteers. Neuroimaging methods included brain MRI in standard modes (T1, T2, T2 d-f (FLAIR), DWI) and MR venosinusography. All patients were surveyed to identify carriers of prothrombogenic polymorphic variants of hemostatic and folate metabolism genes alleles and genotypes. Prothrombogenic polymorphic variants of hemostatic genes were detected in 94% of patients, and the variants of the methionine-homocysteine metabolism genes were observed in 86% of patients. The differences between distributions of alleles and genotypes 5G6754G of the PAI-1 gene, G103T of the FXIIIA1 gene, A66G of the MTRR gene, A2756G of the MTR gene in the group of patients with CVT and in the control group were significant. Allele 4G, genotypes 4G/4G and 5G/4G of 5G6754G polymorphism of the PAI-1 gene; allele T of G103Т polymorphism of the FXIIIA1 gene; allele G and genotype A/G of A66G polymorphism of the MTRR gene; allele G and genotype A/G of A2756G polymorphism of the MTR gene correlated with aseptic CVT. It was concluded that the gene polymorphisms 5G6754G (PAI-1), G103T (FXIIIA1), A66G (MTRR) and A2756G (MTR) carriage increased the risk of aseptic CVT and did not affect the thrombosis clinical manifestations.

Hemostatic Genes Exhibit a High Degree of Allele-Specific Regulation in Liver

Objective Elucidating the genetic basis underlying hepatic hemostatic gene expression variability may contribute to unraveling genetic factors contributing to thrombotic or bleeding disorders. We aimed to identify novel cis-regulatory variants involved in regulating hemostatic genes by analyzing allele-specific expression (ASE) in human liver samples. Study Design Biopsies of human liver tissue and blood were collected from adults undergoing liver surgery at the Sahlgrenska University Hospital (n = 20). Genomic deoxyribonucleic acid (gDNA) and total ribonucleic acid (RNA) were isolated. A targeted approach was used to enrich and sequence 35 hemostatic genes for single nucleotide polymorphism (SNP) analysis (gDNAseq) and construct individualized genomes for transcript alignment. The allelic ratio of transcripts from targeted RNAseq was determined via ASE analysis. Public expression quantitative trait loci (eQTL) and genome-wide association study (GWAS) data were used to assess novelty and importance of the ASE SNPs (and proxies, r 2 ≥ 0.8) for relevant traits/diseases. Results Sixty percent of the genes studied showed allelic imbalance across 53 SNPs. Of these, 7 SNPs were previously validated in liver eQTL studies. For 32 with eQTLs in other cell/tissue types, this is the first time genotype-specific expression is demonstrated in liver, and for 14 ASE SNPs, this is the first ever reported genotype–expression association. A total of 29 ASE SNPs were previously associated with the respective plasma protein levels and 17 ASE SNPs to other relevant GWAS traits including venous thromboembolism, coronary artery disease, and stroke. Conclusion Our study provides a comprehensive ASE analysis of hemostatic genes and insights into the regulation of hemostatic genes in human liver.