The genetic variations among individuals are one of the notable factors determining disease severity and drug response. Nowadays, COVID-19 pandemic has been adversely affecting many aspects of human life. We used the Tehran Cardio-Metabolic Genetic Study (TCGS) data that is an ongoing genetic study including the whole-genome sequencing of 1200 individuals and chip genotyping of more than 15,000 participants. Here, the effect of ACE2 variations by focusing on the receptor-binding site of SARS-CoV-2 and ACE2 cleavage by TMPRSS2 protease were investigated through simulations study. After analyzing TCGS data, 570 genetic variations on the ACE2 gene, including single nucleotide polymorphisms (SNP) and insertion/deletion (INDEL) were detected. Interestingly, two observed missense variants, K26R and S331F, which only the first one was previously reported, can reduce the receptor affinity for the viral Spike protein. Moreover, our bioinformatics simulation of 3D structures and docking of proteins explains important details of ACE2-Spike and ACE2-TMPRSS2 interactions, especially the critical role of Arg652 of ACE2 for protease function of TMPRSS2 was uncovered. As our results show that the genetic variation of ACE2 can at least influence the affinity of this receptor to its partners, we need to consider the genetic variations on ACE2 as well as other genes in the pathways that contribute to the pathogenesis of COVID-19 for designing efficient drugs and vaccines.
Precise prognostic classification of patients and identifying survival subgroups and their associated genes can be important clinical references when designing treatment strategies for cancer patients. Multi-omics and data integration techniques are powerful tools to achieve this goal. This study aimed to introduce a machine learning method to integrate three types of biological data, and investigate the performance of two other methods, in identifying the survival dependency of patients. The data included TCGA RNA-seq gene expression, DNA methylation, and clinical data from 368 patients with colon cancer also we use an independent external validation data set, containing 232 samples. Three methods including, hyper-parameter optimized autoencoders (HPOAE), normal autoencoder, and penalized principal component analysis (PPCA) were used for simultaneous data integration and estimation under a COX hazards model. The HPOAE was thought to outperform other methods. The HPOAE had the Log Rank Mantel-Cox value of 14.27 ± 2, and a Breslow-Generalized Wilcoxon value of 13.13 ± 1. Ten miRNA, 11 methylated genes, and 28 mRNA all by (importance of marginal cutoff > 0.95) were identified. The study demonstrated that hsa-miR-485-5p targets both ZMYM1 and tp53, the latter of which has been previously associated with cancer in numerous studies. Furthermore, compared to other methods, the HPOAE exhibited a greater capacity for identifying survival subgroups and the genes associated with them in patients with colon cancer. However, all of the results were obtained by computational methods, and clinical and experimental studies are needed to validate these results.