{'fraction_non_alphanumeric': 0.029423076923076923, 'fraction_numerical': 0.013685897435897435, 'mean_word_length': 5.164987156688402, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 3, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Drug repositioning for esophageal squamous cell carcinoma

Esophageal cancer (EC) remains a significant challenge globally, having the 8th highest incidence and 6th highest mortality worldwide. Esophageal squamous cell carcinoma (ESCC) is the most common form of EC in Asia. Crucially, more than 90% of EC cases in China are ESCC. The high mortality rate of EC is likely due to the limited number of effective therapeutic options. To increase patient survival, novel therapeutic strategies for EC patients must be devised. Unfortunately, the development of novel drugs also presents its own significant challenges as most novel drugs do not make it to market due to lack of efficacy or safety concerns. A more time and cost-effective strategy is to identify existing drugs, that have already been approved for treatment of other diseases, which can be repurposed to treat EC patients, with drug repositioning. This can be achieved by comparing the gene expression profiles of disease-states with the effect on gene-expression by a given drug. In our analysis, we used previously published microarray data and identified 167 differentially expressed genes (DEGs). Using weighted key driver analysis, 39 key driver genes were then identified. These driver genes were then used in Overlap Analysis and Network Analysis in Pharmomics. By extracting drugs common to both analyses, 24 drugs are predicted to demonstrate therapeutic effect in EC patients. Several of which have already been shown to demonstrate a therapeutic effect in EC, most notably Doxorubicin, which is commonly used to treat EC patients, and Ixazomib, which was recently shown to induce apoptosis and supress growth of EC cell lines. Additionally, our analysis predicts multiple psychiatric drugs, including Venlafaxine, as repositioned drugs. This is in line with recent research which suggests that psychiatric drugs should be investigated for use in gastrointestinal cancers such as EC. Our study shows that a drug repositioning approach is a feasible strategy for identifying novel ESCC therapies and can also improve the understanding of the mechanisms underlying the drug targets.

Esophageal cancer (EC) remains a significant challenge globally, having the 8th highest incidence and 6th highest mortality worldwide. Esophageal squamous cell carcinoma (ESCC) is the most common form of EC in Asia. Crucially, more than 90% of EC cases in China are ESCC. The high mortality rate of EC is likely due to the limited number of effective therapeutic options. To increase patient survival, novel therapeutic strategies for EC patients must be devised. Unfortunately, the development of novel drugs also presents its own significant challenges as most novel drugs do not make it to market due to lack of efficacy or safety concerns. A more time and cost-effective strategy is to identify existing drugs, that have already been approved for treatment of other diseases, which can be repurposed to treat EC patients, with drug repositioning. This can be achieved by comparing the gene expression profiles of disease-states with the effect on gene-expression by a given drug. In our analysis, we used previously published microarray data and identified 167 differentially expressed genes (DEGs). Using weighted key driver analysis, 39 key driver genes were then identified. These driver genes were then used in Overlap Analysis and Network Analysis in Pharmomics. By extracting drugs common to both analyses, 24 drugs are predicted to demonstrate therapeutic effect in EC patients. Several of which have already been shown to demonstrate a therapeutic effect in EC, most notably Doxorubicin, which is commonly used to treat EC patients, and Ixazomib, which was recently shown to induce apoptosis and supress growth of EC cell lines. Additionally, our analysis predicts multiple psychiatric drugs, including Venlafaxine, as repositioned drugs. This is in line with recent research which suggests that psychiatric drugs should be investigated for use in gastrointestinal cancers such as EC. Our study shows that a drug repositioning approach is a feasible strategy for identifying novel ESCC therapies and can also improve the understanding of the mechanisms underlying the drug targets.

Introduction
There are two major subtypes of Esophageal cancer (EC), esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC) (Ye et al., 2021). In China, more than 90% of esophageal cancer cases are ESCC . EC as a whole remains a significant challenge globally, having the 8 th highest incidence and the 6th highest mortality worldwide killing over 500,000 people in 2020 (Sung et al., 2021). A major driver of the high mortality rate is likely due to the fact that there are very few effective therapeutic options for EC patients. In recent years, there has been a significant increase in survival for many cancers, largely due to the availability of targeted therapies. For EC, however, targeted therapies are yet to make a significant impact on patient survival. Consequently, patients are often relying on more traditional therapies such as chemotherapy and surgical resection. In-order-to increase patient survival, novel therapeutic strategies for EC patients must be devised. Unfortunately, the development of novel drugs also presents its own significant challenges as most novel drugs do not make it to market due to lack of efficacy or safety concerns. Therefore, it is more time and cost effective to identify existing drugs, that have already been approved for treatment of other diseases, which can be repurposed to treat EC patients. This can be achieved using a drugs gene signature, the alterations in gene expression as a result of exposure to the drug. The gene signature of a drug indicates the underlying biological pathways and mechanisms that are involved in the therapeutic effect of the drug. With this knowledge, we can then identify candidate drugs which have gene signatures capable of reversing aberrant gene expression patterns observed in disease-states to those observed in normal cells. This gene signature-based approach has been adopted by previous research to identify drugs that can be repositioned to treat a variety of diseases including, but not limited to, cancer, Alzheimer's, hyperlipidaemia, hypertension, and inflammatory disease (Corbett et al., 2012;Hall et al., 2014;Guney et al., 2016;Subramanian et al., 2017;Cheng et al., 2018;Carvalho et al., 2021;Wu et al., 2022). To date, drug repositioning to target gene signatures has primarily involved identifying directly overlapping drug genes and disease genes (herein referred to as overlap analysis) (Subramanian et al., 2017;Wang et al., 2018;Chen et al., 2022). More recently, network analysis has been greatly employed in this area as it offers distinct advantages over more traditional statistical methods. This is due to the fact that the models that can be built with this methodology are an excellent way to capture a molecules relationship with other molecules. In particular, nodes can be used to represent multiple entities such as genes, molecules, proteins, etc, and the edges can also represent a vast array of information such as mode-of-actions (MoAs), underlying mechanisms, or functional similarities (Jarada et al., 2020) Hence, network-based methods can accurately represent the biological mechanisms which are driving diseases (Barabási et al., 2011). As a result, network-based drug repositioning can identify drugs which target the underlying biology of the disease. It is worth noting, however, that other methods of computational drug repositioning have also been adopted, such as Data Mining and Machine Learning. An excellent review of the different methodologies, as well as their advantages and disadvantages has recently been published (Jarada et al., 2020). Due to the success of drug repositioning overall, and the absence of effective treatments for ESCC, it has been proposed that this method be used to identify novel treatment strategies for ESCC. However,

FIGURE 1
Drug Repositioning Analysis Methodology. The initial step of the analysis included differential gene expression on previously published array data from GEO (accession: GSE23400). DEGs were then used to identify key driver genes in a weighted key driver analysis. The key driver genes were then used as input in 2 arms; overlap analysis and network analysis. Qualify control was performed to filter out erroneous results and identify candidate drugs. Drugs which were common to both arms were considered robust and considered ESCC Repositioned Drugs.
Frontiers in Genetics frontiersin.org 02 these studies have largely, though not completely, been limited to testing existing cancer drugs in vitro with drug screening methods . Herein, we adopt both a network-based and overlap-based drug repositioning methodology, to identify existing drugs that can specifically target the aberrant expression profile of ESCC and impede oncogenesis. To do this, we used previously published data for in-silico drug repositioning analysis utilising the PharmOmics webserver . The repositioning analysis consisted of two arms, the 'overlap analysis' arm and the 'network analysis' arm ( Figure 1), which utilise two methods of drug repositioning.

Identification of DEGs
The dataset GSE23400 was downloaded using the GEOquery R package function getGEO. In total, 167 DEGs were identified between ESC and normal samples (Details of the differential gene expression analysis can be found in the methods section). Of which, 65 were upregulated and 102 were downregulated. The top 5 most upregulated genes are MMP1, SPP1, POSTN, COL1A1, and JUP. The top 5 most downregulated genes are CRISP3, MAL, CRNN, SCEL, CLCA4. The top 25 up-regulated genes can be observed in Table 1, whereas the top 25 downregulated genes can be observed in Table 2.

Functional and pathway enrichment analyses
Functional and pathway enrichment analyses were performed used the 'clusterProfiler' R package. Gene Set Enrichment Analysis (GSEA) was performed with Gene Ontology (GO) (hereafter referred to as GSEA-GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway (hereafter referred to as GSEA-KEGG). GSEA-GO analysis was performed with the gene set categories Biological Process (BP), Cellular Component (CC), and Molecular Function (MF), which identified 253, 36, 25 enriched gene  Frontiers in Genetics frontiersin.org 03 sets, respectively. Numerous BP gene sets identified by the analysis are related to extracellular matrix and cell differentiation. Ranking BP analysis by adjusted p-value, the top 5 most enriched gene sets are cellular component organization, cellular component organization or biogenesis, extracellular matrix organization, extracellular structure organization, multicellular organism development. According to adjusted p-value, the top 5 most enriched CC category are endoplasmic reticulum lumen, external encapsulating structure, extracellular matrix, fibrillar collagen trimer, banded collagen fibril. Furthermore, the top 5 categories in the MF analysis identified extracellular matrix structural constituent, extracellular matrix structural constituent conferring tensile strength, protein-containing complex binding, cell adhesion molecule binding, glycosaminoglycan binding. The full results for BP, CC, MF can be observed in Supplementary Tables S1-S3, respectively. Gene Set Enrichment Analysis of KEGG (GSEA-KEGG) identified 12 enriched gene sets (Supplementary Table  S4), including those previously identified as ESCC-related, such as Focal adhesion, ECM-receptor interaction, PI3K-Akt signalling pathway.

Weighted key driver analysis
Weighted Key Driver Analysis (wKDA) was performed using Mergeomics webserver. In this analysis, genes which possess a local network neighbourhood that have a significant enrichment of genes that are ESCC-associated are considered key drivers (KDs) (Ding et al., 2021). The analysis identified 89 key driver genes which were then filtered to select those which possessed an FDR <0.05, ensuring that only the strongest KDs are used in subsequent analyses. This resulted in 39 key driver genes (Supplementary Table S5). The top 10 key driver genes are NCAPG, PLG, NUSAP1, COL17A1, ASPM, TOP2A, ITGB3, P4HB, TTK, and COL7A1.

Repositioned drugs
Drug repositioning analysis was performed using both the Overlap Drug Repositioning and the Network Drug Repositioning modules from PharmOmics . The potential drugs from the analysis, were then filtered to

Drug validation
To validate our findings, we performed a literature search to determine whether any of the drugs identified by our analysis are currently used in ESCC treatment (Table 4). We found that 7 of the top 10 repositioned drugs, according to z-score, are already used to treat ESCC or have been shown to demonstrate efficacy in clinical trials. Candidate drugs were then validated using Binding DB (Gilson et al., 2016). Each drug was searched in the database to ascertain whether they bind to proteins known to be involved in ESCC. We found that 21 out of 25 repositioned drugs have a strong binding affinity to proteins that have been associated with ESCC in some manner previously (Table 5). Additionally, we performed a literature search to assess whether there is any biological evidence (in vitro or in vivo) that demonstrates efficacy or establishes a plausible mechanism by which the novel repositioned drugs could be beneficial for ESCC patients (Table 6). We found that all of our novel ESCC drugs, except for Venlafaxine, target pathways or proteins which have been demonstrated to drive oncogenesis in several cancers, including ESCC. Therefore, these drugs should be able to target the underlying biological processes driving oncogenesis Tamoxifen No -- Frontiers in Genetics frontiersin.org 05 in ESCC and inhibit proliferation and/or initiate apoptosis in ESCC.

Discussion
ESCC is one of the most common malignancies and possess a significant mortality rate worldwide. This is largely due to late diagnosis and scarcity of efficacious treatment strategies upon being diagnosed (Feng et al., 2021). To address this, we performed a disease-based drug repositioning analysis with previously published ESCC gene expression data from paired patient samples. Differential gene expression analysis data identified 167 differentially expressed genes (DEGs) which were then used in wKDA and identified 39 key driver genes (KDGs). The genes with the highest absolute logFC identified by our differential gene expression analysis are MMP1, CRISP3, MAL, CRNN, SCEL. The most upregulated gene, MMP1, encodes a protein involved in the breakdown of the extracellular matrix (ECM) by cleaving collagens and other molecules. The most downregulated gene, CRISP3, encodes a protein located in the ECM and thought to be involved in cellular matrix remodelling (Ribeiro et al., 2011). The wKDA identified 39 significant driver genes for ESCC. Amongst the top 10 most significant KDGs, NUSAP1, COL17A1, ITGB3 and COL7A1 are involved in ECM maintenance. For example, the 4 th most significant key driver gene, COL17A1, encodes a protein involved in cell-matrix adhesion (Jones et al., 2020). Taken together, these results suggest that alterations in ECM are an important driver of ESCC oncogenesis . Furthermore, KEGG analysis found both Focal adhesion and ECM-receptor interaction to be the 3rd and 4th most enriched term, respectively. This is in line with previous research that indicates higher levels of Serum human relaxin 2 (H2 RLN), a protein involved in ECM, collagen, and matrix metalloproteinase is associated with worse prognosis, including higher clinical stage and poorer survival (Ren et al., 2013;Napier et al., 2014).
Using 39 KDGs in an ESCC drug repositioning analysis, we identified 25 drugs that are predicted to have therapeutic effect in ESCC. Of which, 7 are either currently used in the clinic or have been used in clinical trials and 2 have shown efficacy in vitro or in vivo. Importantly, those which have been used in clinical trials have demonstrated efficacy particularly when used in combination with other drugs, such as chemotherapy. This is not surprising, however, as combination therapy has long been a standard practice in cancer therapy, including for ESCC where the current first-line treatment regimen is a combination of 5-fluorouracil and cisplatin (Hiramoto et al., 2018;Hirano and Kato, 2019). Each repositioned drug was validated in-silico using the drug binding database BindingDB, to identify which drug targets have previously been associated with ESCC (Table 5). Significantly, 21 of the 25 repositioned drugs have targets that have previously been associated with ESCC in some manner, which demonstrate the  Kim et al. (2015) Anastrozole Aromatase Inhibition Has been used with Anti-Fibroblast growth factor receptor 1 (FGFR1) drug in breast cancer. Evidence that FGFR1 can be used as a independent prognosis marker in ESCC and anti-FGFR1 decreases proliferation via MEK-ERK downstream pathways Milani et al., 2009;Chen et al., 2017 Frontiers in Genetics frontiersin.org robustness of our findings. To further validate our findings, we performed a literature search on the novel repositioned drugs to examine whether there is an underlying biological mechanism which would justify the drugs appearance in the results (Table 6). We found that almost all of the repositioned drugs have been shown to demonstrate anti-cancer effects in multiple cancers, most notably breast cancers and non-small cell lung carcinoma (NSCLC). Interestingly, many of the repositioned drugs target specific pathways; EGFR (Erlotinib, Crizotinib, and Lapatinib), estrogen signalling (Tamoxifen, Fulvestrant, Hydrocortisone, and Anastrozole) and TRAIL-mediated apoptosis (Azithromycin and Anastrozole) pathways, suggesting that these pathways are key drivers of ESCC. This is in line with previous research which identified the EGFR AND ER pathways as drivers of ESCC oncogenesis and metastasis and have also been associated with patient outcome (Maron et al., 2020). Crucially, some of these drugs have been shown to have therapeutic potential in vitro. For example, Lapatinib, which acts through EGFR and HER2 has been shown to be efficacious in ESCC patient-derived xenografts (Rong et al., 2020). The potential mechanisms by which novel drugs identified by our study can be observed in Table 6. It is also worth noting that there are 2 anti-depressants present in our results, Venlafaxine and Nefazodone. These results are particularly interesting as it has recently been shown that psychiatric drugs offer potential as anti-cancer therapeutics (Loehr et al., 2021). Moreover, a recent review has specifically addressed the need to investigate psychiatric drugs for treatment of gastrointestinal cancers (Avendaño-Félix et al., 2020). We hypothesise that Crizotinib, Lapatinib, and Dasatinib are amongst the drugs with the most potential. Particularly Crizotinib and Lapatinib are of note as they target the EGFR pathway which is already targeted in ESCC treatment with Erlotinib. Dasatinib also has high potential due to targeting Tyrosine-and threonine-specific cdc2-inhibitory kinase and CDK1, proteins known to be involved in ESCC, and also due to displaying efficacy in cell lines (Chen et al., 2015;Zhang et al., 2019). There are several limitations to our study, however, most notably that due to limited data availability, the sample size of patient samples was relatively small. We were unable to stratify patients according to subtype of ESCC. This means that the analysis is focussed on ESCC as whole and does not take into consideration specific subtypes. Moreover, as multiple drugs identified in our study are more efficacious when in combination with another drug, it would be beneficial to know what other drugs should be used in combination with the novel therapeutics identified. However, this analysis does not predict drug combinations that would be effective in treating ESCC.
On the other hand, our study has several strengths. To our knowledge, this is the first study to adopt a primarily computational approach to perform drug repositioning analysis in ESCC. Particularly, there are studies that have a computational component, but they do not use patient samples to identify drugs based on network analysis of differentially expressed genes . Moreover, this is also the first to adopt Pharmomics unique network analysis to perform the analysis on ESCC. Furthermore, as Pharmomics contains >18000 species/tissue-specific gene signatures for 941 drugs and chemicals, it provides a larger scope of potential drugs compared to other studies in ESCC. Another strength of the study is that we used a two armed approach to ensure robust findings as each repositioning methodology has its own strengths. The overlap-based repositioning allows us to identify drugs which target the KDGs whereas the network-based repositioning allows for insights into the molecular and mechanistic therapeutic effects of the drugs. As this specific form of network-based repositioning is unique to PharmOmics, our study can provide valuable insights into the underlying molecular mechanisms driving ESCC. Another strength of our study is the consistency of our results with previously published literature. DEGs which displayed the highest absolute logFC were consistent with previously published literature including MMP1, SPP1, COL1A2, and COL1A1 amongst the top upregulated genes and CRISPR3, MAL, TMPRSS11E, and CRNN amongst the top downregulated genes (Feng et al., 2021;Song et al., 2021). Indeed, the gene with the highest absolute logFC, MMP1, is already known to be associated with ESCC oncogenesis . Additionally, higher MMP1 is associated with poorer prognosis (Feng et al., 2021). Moreover, the most significant key driver genes (KDGs) identified by our wKDA are consistent with previously published studies Yu-jing et al., 2020;Wang M. et al., 2021). Significantly, multiple drugs identified by our analysis target key pathways known to be involved in ESCC oncogenesis and metastasis.

Conclusion
Herein we utilised in silico disease-based drug repositioning to identify novel therapeutics for esophageal squamous cell carcinoma. Amongst 25 potential repositioned drugs identified in our study, 9 are currently used in the clinic or have shown promising results in clinical trials in combination with other treatments. Crucially, we identified 16 novel therapeutic strategies which possess a strong biological rationale for use in ESCC patients. Our study shows that drug repositioning approach is a feasible strategy in ESCC therapies and can improve the understanding of the mechanisms of the drug targets.

Materials and methodology Data acquisition and identification of DEGs
Dataset was acquired from the Gene Expression Omnibus under accession code GSE23400 using the GEOquery R package function getGEO. The dataset consists of 53 paired patient samples from Esophageal squamous cell carcinoma (ESCC) patients. Additionally, 14,335 genes were in the dataset. Differentially expressed genes (DEGs) between paired tumour and non-tumour samples were identified using the limma package. The log-fold change (logFC) was calculated for DEGs. Genes with absolute logFC >1.5 and adjusted p-value < 0.01 were considered significant and used in subsequent analyses.

Functional and pathway enrichment analyses
The DEGs identified above were analysed using the clusterProfiler R package in order to identify biological annotations from the Gene ontology (GO) functional enrichment and Kyoto Encyclopaedia of Genes and Genomes (KEGG). The GO analysis was performed for biological process (BP), cellular component (CC) and molecular function (MF). An adjusted p-value < 0.05 was considered as statically significant for all analyses.

Weighted key driver analysis
Weighted key driver analysis (wKDA) was performed on DEGs using Mergeomics webserver to identify key driver genes (KDGs). wKDA has higher accuracy than standard key driver analysis as it considers edge weight information. The network used in the analysis was STRING PPI Network and default parameters were used (Search depth of 1, Undirected Edges, Min Hub Overlap of 0.33, and edge factor of 0.0). Genes which had an FDR <0.05 were considered as significant KDGs and used in subsequent analyses.

Drug repositioning analysis
By the Pharmomics webserver, Drug Repositioning Analyses was performed using genes obtained from the wKDA analysis. The analysis consisted of two arms: the Overlap Drug Repositioning and ADR Analysis (Overlap-DR) arm and the Meta-Signature Network Drug Repositioning and ADR Analysis (Meta-Net-DR) arm. The network analysis adopted by Pharmomics uses a network proximity measure between drug DEGs and disease-related genes that has been adopted previously for protein-network-based analysis. Specifically, tissue-specific Bayesian gene regulatory networks (BNs) are used and then the mean shortest distance between drug DEGs and disease genes are tested. Hence, it combines species and tissue specific in vivo drug signatures with gene networks to identify connections between disease genes and known drug targets. On the other hand, overlap analysis adopted by Pharmomics is largely similar to that which has been adopted previously, and assesses direct overlap between input genes and drug gene signatures. To do so, the Jaccard score, gene overlap fold enrichment, and Fisher's exact test p values as measures of direct gene overlap are calculated. This analysis is based upon the premise that if disease and drug signatures target similar pathways then they would more than likely have gene overlaps and/or connect extensively in a gene network. Meta-Signature Network Drug Repositioning and ADR Analysis was performed using the multi-tissue network. In Overlap-DR, Jaccard score was used to measure the similarity between the 39 KDG's gene networks and the drug target gene networks. In Meta-Net-DR, the connectivity of the gene network between drug signatures from PharmOmics and the KDs is used. The z-score of each drug is calculated which represents the distance between the KD network and the PharmOmics drug network. The smaller the z-score, the closer the distance between the networks. The output from these analyses were considered as possible repositioned drugs and were then filtered in Drug Candidate Selection to identify ESCC repositioned drugs.

Drug candidate selection
Repositioned drugs from both Pharmomics analyses were used as candidates to identify potential drugs for ESCC. Candidate drugs from Overlap-DR results were filtered to keep drugs with an adjusted P-value < 0.05, species equal to Homo sapiens, and a within species rank >0. The mean Jaccard score was then calculated and drugs with a Jaccard score less than the mean were removed. Subsequently, the drugs were sorted according to 'Drug Name', 'Within Species Rank', 'Jaccard Score', and 'P-value' and duplicate drugs were removed, keeping only the highest-ranking occurrence of each drug. Candidate drugs from Meta-Net-DR were filtered to keep drugs with adjusted p-value < 0.05. Candidate drugs were then sorted according to 'Drug Name' and 'Rank' and then duplicate drugs were removed, keeping only the highest-ranking occurrence of each drug. The filtered results from Overlap-DR and Meta-Net-DR were then compared to extract candidate drugs common to both arms of analysis. Drugs common to both arms were considered ESCC Repositioned Drugs. 'Study', 'Jaccard Score', 'Odds Ratio', 'Adj. P-Value', 'Within Species Rank' data from the Overlap-DR analysis and 'z-score' from Meta-Net-DR was used to construct the final ESCC Repositioned Drugs Repositioned Drugs were then sorted according to z-score (Table 3).

Drug candidate validation
In order to validate the repositioned drugs that were identified by the analysis, we performed a literature search to ascertain whether the drugs have previously been used in ESCC treatment (Table 4). Each drug was then investigated using the drug binding database Binding DB. For each ESCC Repositioned Drug, we identified which proteins they display a high binding affinity to. We then performed a literature search on these proteins, using Google Scholar and PubMed, to ascertain whether or not they have previously been shown to be ESCCrelated in vitro or in vivo (Table 5). Finally, we performed a literature search on novel drugs identified by our analysis to elucidate the underlying biological processes and causal mechanisms which would explain why it is predicted to have therapeutic utility (Table 6).

Data availability Statement
Publicly available datasets were analyzed in this study. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement
The studies involving human participants were reviewed and approved by Human Subjects Ethics Sub-Committee at City University of Hong Kong (Reference number 2-11-201810_02). The patients/participants provided their written informed consent to participate in this study.

Author contributions
NL, W-KS, and KC envisioned and directed the project, assisted in writing manuscript, and performed data interpretation. AB and RH developed the analytical pipeline and interpreted data. AB performed the analysis and wrote the manuscript. QH assisted in analytical pipeline development and contributed towards the manuscript. KC coordinated the development of the analysis and revised the manuscript.

Funding
This project was funded by The Jockey Club College of Veterinary Medicine and Life Sciences (JCC) Interdisciplinary PhD Programme at City University of Hong Kong in Collaboration with Cornell University. This work was supported by the City University of Hong Kong New Research Initiatives/Infrastructure Support from Central (APRC; grant number 9610401).

Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene. 2022.991842/full#supplementary-material SUPPLEMENTARY FIGURE S1 This figure shows the enriched gene set categories for the Gene Set Enrichment Analysis for Biological Processes.

SUPPLEMENTARY FIGURE S2
This figure shows the enriched gene set categories for the Gene Set Enrichment Analysis for Cellular Component.

SUPPLEMENTARY FIGURE S3
This figure shows the enriched gene set categories for the Gene Set Enrichment Analysis for Molecular Function.
Frontiers in Genetics frontiersin.org