Determination of whole-genome expression differences in larynx cancers
Expression differences in laryngeal cancer
Authors
Abstract
AimMore than 200.000 new laryngeal cancer cases have been diagnosed worldwide, annually. The incidence and prevalence have increased during the past decades. It also has a high mortality rate. To have a better knowledge about this particular cancer type, we aimed to conduct a comparative analysis of whole genome expression differences between laryngeal squamous cell carcinoma and adjacent normal tissues
MethodsTissue samples were collected from specimens of laryngeal carcinoma and normal epithelium tissues adjacent to the carcinoma. RNAs isolated from these specimens were used for genome-wide gene expression analysis using microarrays. Genes that were expressed significantly differently in the tumor sample compared to normal tissues were identified. Pathway analysis of these genes was also performed.
ResultsIt was observed that 24 genes were significantly differentially expressed in cancer tissue. Expressions of MMP1, MMP12, S100A2, S100A3, CK14, CK16, SLC2A1, ITGA6, CEP55, KLK6, LAMC2, IL1F9, TP63, VSNL1, CXCL1, COL4A5, COL4A6, FSCN1, KRT6B, KRT17, WDR66 and ACOT7 genes were increased, while expressions of ITM2A and CFD genes were decreased.
ConclusionSome of the genes mentioned above are known to be involved in the etiology of laryngeal cancer, but to the best of our knowledge, 10 genes have been associated with laryngeal cancer for the first time in our study. These genes might be useful as biomarkers in the clinic, for early diagnosis, prognosis and personalized targeted therapy.
Keywords
Introduction
laryngeal carcinoma remains a considerable cause of morbidity and mortality with ~180,000 new cases and 99,000 deaths, accounting for 50% of mortality worldwide in 2020.1 several genetic and environmental factors have been identified to be associated with the development and progression of laryngeal squamous cell carcinoma (lscc). among them, the major etiological agents are tobacco and alcohol consumption. in recent years, these agents have been shown to cause molecular changes and a sequence of events in cells, including oncogene activation, inhibition of tumor suppressor genes, and changes in gene expression profiling. clinical outcomes in laryngeal cancer patients are thought to be based on clinicopathological features; however, few studies have indicated that gene expression alterations, as well as genomic and epigenetic changes, can also be used to predict prognosis. due to this reason, molecular characterization of laryngeal cancers seems to be extremely important. certain success has been achieved in molecular studies on the etiology of larynx cancer, but previous studies on larynx cancers were uncertain and incomprehensive.2,3
lately, gene expression profiling (gep) using microarrays has been identified as a hopeful method to identify variations occurring in the pathophysiology of complex diseases such as cancers. identification of gene expression changes associated with laryngeal cancer gives us information about the mechanism of the disease and the pathways associated with this disease. thus, new biological and prognostic markers can be identified.4
there are many studies in the literature aimed at clarifying the etiology of the disease, but the pathogenesis of this neoplasm and associated pathways are still not fully understood. the lack of significant improvement in survival rates over the past 25 years is proof of this.5 that is why we aimed to conduct a comparative analysis of whole genome expression differences between laryngeal squamous cell carcinoma and adjacent normal tissues. these efforts will ensure a better understanding of molecular mechanisms underlying laryngeal carcinoma, as well as enable identification of novel biomarkers and therapeutic targets.
Materials and Methods
Tissue SamplesFresh tumor and adjacent normal tissue biopsy specimens were taken from n = 12 patients diagnosed with larynx cancer between March 2015 and April 2016. Clinical characteristics of the patients are shown in Table 1. These tissues were placed in RNA Later solution and stocked at -80 C until use. Informed consent was obtained from patients or their families. This study was approved by the ethics committee of Necmettin Erbakan University (decision number 2015/146) and financial support was provided by the Scientific Research Projects of Necmettin Erbakan University (Project Number:151518008).
RNA IsolationRNA isolation was performed by Norgen Total RNA Purf. Kit (cat.no. 25700). RNA concentration was detected by Thermo Scientific Nonodrop 2000c. Then, RNA quality and quantity were identified by Agilent RNA 6000 Nano Kit (reorder-no 5067-1511). Thus, RNA was purified from other molecules with the same absorbance.
Synthesis of CDNA and In Vitro Transcription of Biotin-CRNAThe obtained RNA was converted to cDNA using the TargetAmp™-Nano Labeling Kit (Illumina® Expression BeadChip®) before the hybridization step. First-strand cDNA synthesis was catalyzed by Super Script III Reverse Transcriptase. The produced cDNA was converted to double-stranded cDNA for use in the subsequent in vitro transcription reaction. Then, in vitro transcription was performed to obtain biotinylated cRNA. The synthesized biotin-labeled RNA was cleaned up again (Norgen RNA Clean-Up and Conc. Micro Kit 23600,43200) and cRNA levels of all samples were equalized before being sent to the microarray platform.
Microarray AnalysisThe biotin-labeled RNA was processed through hybridization, washing, labeling with streptavidin, and drying before installation. The labeled RNA strand was hybridized to the bead on the BeadChip (HumanHT-12 v4.0 Exp. BeadChip Kit) containing the complementary gene-specific sequence. This kit presents a comprehensive analysis of genome-wide expression. Each channel on the BeadChip contains approximately 47,000 different bead types distributed throughout the genome. Contents of the kit were created with reference to the NCBI RefSeq database. Hybridization was performed at 58 C for 20 h. BeadChips were washed after overnight hybridization. Unbound RNAs were removed. BeadChips were labeled with streptavidin, centrifuged immediately after washing, and stored in a dark, ozone-free environment until scanning. Finally, BeadChips were read in a microarray laser reader (Illumina® Whole-Genome Gene Expression Direct Hybridization Assay system).
Bioinformatics AnalysisPreliminary analyses and quality controls of the obtained data were performed using the Genome Studio software of the Illumina iScan Microarray system. Fold change (FC) is a strong marker indicating statistical significance of genes observed in the gene expression profile. By FC analysis, changes larger than two-fold were accepted as increase or decrease for each probe. Heat-map analysis for differentially expressed genes was generated and GO (gene ontology) analysis was performed.
Ethical ApprovalEthics Committee approval for the study was obtained.
Results
In this study, we identified that 14,294 genes were differentially expressed between laryngeal squamous cancer tissues and non-neoplastic tissues. When the FC value is considered to be 2 for the significance of the expression differences, increased expression of 22 genes (MMP1, MMP12, S100A2, S100A3, CK14, CK16, SLC2A1, ITGA6, CEP55, KLK6, KRT17, LAMC2, IL1F9, TP63, VSNL1, CXCL1, COL4A5, COL4A6, FSCN1, KRT6B, WDR66 and ACOT7) and decreased expression of two genes (ITM2A and CFD) were observed. Fourteen of 24 genes were known to be involved in the etiology of laryngeal tumors in the literature, but as far as we know, the effect on laryngeal cancer of the remaining 10 genes (VSNL1, CXCL1, COL4A5, COL4A6, FSCN1, KRT6B, WDR66, ACOT7, ITM2A, CFD) was identified for the first time in our study. Differentially expressed genes between LSCC and adjacent normal tissues are shown in Table 2. The hierarchical clustering of the expression of these genes is shown in Figure 1.
Also, to determine the biological functions of differentially expressed genes, WikiPathway functional enrichment analysis was carried out using the Gene Spring gx (Agilent) program. WikiPathways was established to contribute to the information about the pathway and to facilitate the access by the biology community. WikiPathways is a new database including and complementing databases such as KEGG, Reactome and Pathway Commons.
The pathways such as cell cycle, G1 to S cell cycle control, cell cycle checkpoints, regulation of DNA replication, DNA replication, degradation of the extracellular matrix were upregulated; however, the pathways such as electron transport chain, oxidative phosphorylation, integrin mediated cell adhesion, the citric acid cycle and respiratory electron transport, adipogenesis, mitochondrial translation, fatty acid triacylglycerol and ketone body metabolism were downregulated. The WikiPathway analysis of the most affected pathways is shown in Table 3.
Discussion
Cancers occur as a consequence of accumulation of genetic aberrations. Almost every neoplasm has its own unique molecular changes. Characterization of these changes is extremely important in matters of clinical care and the development of personalized treatment strategies.
Laryngeal carcinomas have been identified as aggressive tumors because of their high recurrence and metastasis rate. The first step for regional recurrence and distant metastasis is the deterioration of the extracellular matrix components. Matrix metalloproteinases (MMP) are a family of enzymes that play a key role in this process. Overexpression or dysregulation of these enzymes is known to be related to various types of cancer.6 Krecicki et al. demonstrated strong immunoreactivity of MMP1 in 36 (72%) of 50 laryngeal squamous cell carcinoma cases using immunohistochemical methods.7 Liu et al. reported a difference in MMP12 expression between tumor and normal tissues, also metastatic and non-metastatic tumors.8 Kallikrein-related peptidase 6 (KLK6), a member of the serine protease family, is involved in cellular processes such as degradation of the extracellular matrix in a similar way to MMP. In a study of n = 162 patients with head and neck tumors, KLK6 overexpression was detected in 42.6% of the cases.9,10 In addition, overexpression of the ITGA6 gene, which is involved in cell-cell and cell-matrix interactions, has been shown to be effective in tumor invasion and metastasis in many cancers, including LSCC.11,12 Similarly, laminin encoded by the LAMC2 gene, which plays a role in cell differentiation, migration, and metastasis, was found to be overexpressed in squamous laryngeal cancer cases by microarray technique.13 We also detected upregulated expression in MMP1, MMP12, KLK6, LAMC2, and ITGA6 genes in all of our cases.
S100 proteins are components of the Epidermal Differentiation Complex and are known to be related with skin diseases and human cancers, including squamous cell laryngeal cancer, lung, ovarian, renal, colorectal, skin, and gastric carcinomas. In the study by Tyszkiewicz et al., while S100A1 and S100A4 genes were down-regulated, S100A2, S100A3, S100A4, and S100A11 genes were significantly up-regulated in n = 93 cases of head and neck cancer. In our study, S100A2 and S100A3 genes were upregulated in all cases.14
CK14, CK16, KRT17, and KRT6 genes, encoding members of the keratin family, are responsible for the structural integrity of epithelial cells. Lauriola et al. reported that CK14 expression difference was detected in 54 of 62 laryngeal carcinoma cases by immunocytochemistry, and this expression was strictly associated with S100A2.15 Elazezy showed that keratin 16 overexpression might be related to more aggressive breast cancer.16 Khanom R. et al. demonstrated that KRT6 and KRT17 were overexpressed in oral squamous cell carcinoma (OSCC) cell lines.17
CEP55 (centrosomal protein 55kDa) has been identified as a tumor-associated antigen functioning as a regulator in the PI3K/AKT pathway. Its upregulation has been associated with poor prognosis in various neoplasms.18 In laryngeal squamous cell carcinoma cases, CEP55 expression was increased in early and advanced stages.19 Our results also support the oncogenic role of CEP55 in LSCC. IL1F9 (interleukin 1 family, member 9) activates NF-kappa B and PI3K pathways, increasing cell growth, differentiation, and survival while suppressing apoptosis.20 PI3K pathway is upregulated in over 90% of head and neck carcinoma cases.21
In recent years, Glut-1 protein encoded by the SLC2A1 gene has become popular in cancer research. Starska et al. observed SLC2A1 gene expression upregulated in 83% of 106 laryngeal carcinoma cases, associated with poor prognosis.22 Another study showed increased SLC2A1 expression in 30 of 38 head and neck carcinoma patients.23 In our study, SLC2A1 expression was significantly increased, consistent with the literature. TP63, P53 tumor suppressor gene homolog, is highly expressed in various cancers, including LSCC.24
Although genes such as VSNL1, CXCL1, COL4A5, COL4A6, FSCN1, KRT6B, KRT17, WDR66, ACOT7, ITM2A, and CFD are known to be altered in other cancers, they have not been described before in LSCC. Type IV collagen protein plays an important role in migration and adhesion. Degradation of type IV collagen is related to cancer progression, invasion, and metastasis. In colorectal cancer patients, loss of COL4A5 and COL4A6 expression was observed due to promoter hypermethylation.25 In esophageal squamous cell carcinoma specimens, increased COL4A5 and COL4A6 expression was found.26 FSCN1 encodes a protein involved in cell migration, motility, adhesion, and interactions. Overexpression of FSCN1 increases metastasis via enhanced motility.27 Increased FSCN1 expression was detected in n = 129 cases of oral and oropharyngeal carcinoma.28
VSNL1 has been identified as a tumor suppressor in esophageal cancer, but as an oncogene in neuroblastoma. It can also be used as a marker for Alzheimer’s disease and other neurodegenerative diseases.29,30 CXCL1 gene encodes a protein involved in inflammation and neutrophil chemoattraction. Abnormal expression is associated with melanomas.31 ACOT7 gene participates in several neoplasms including melanomas, and its upregulation is linked to poor prognosis in acute myeloid leukemia.32,34
WDR66 is a large family of proteins involved in signal transduction, apoptosis, cell cycle control, transcription regulation, and autophagy.35 Overexpression has been found in gastric, thyroid, lung, and esophageal carcinomas, suggested as tumor markers.36,38 ITM2A’s role in tumorigenesis is unknown, but decreased expression was observed in ovarian cancer tissues, suggesting a tumor suppressor role.39 CFD has anti-inflammatory effects via the complement system. Ye H et al. detected decreased CFD expression in tongue squamous cell carcinoma cases.40
In this study, analysis of expression-altering genes and pathway analysis were performed. The cell cycle and extracellular matrix degradation pathways were over-represented, while the electron transport chain pathway was under-represented. All affected pathways are shown in Table 3.
Conclusion
We have identified several differentially expressed genes and described their functions in the cell and the affected pathways in laryngeal squamous cell carcinoma. These genes may be used as new biomarkers for diagnosis, prognosis and therapy of LSCC. It is also believed that these genes will shed light on the studies that we wish to do in the future.
Declarations
Animal and Human Rights Statement
All procedures performed in this study were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. No animal or human studies were carried out by the authors for this article.
Data Availability
The datasets used and/or analyzed during the current study are not publicly available due to patient privacy reasons but are available from the corresponding author on reasonable request.
Conflict of Interest
None of the authors received any type of financial support that could be considered potential conflict of interest regarding the manuscript or its submission.
Funding
None.
References
-
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7-30. doi:10.3322/caac.21590
-
Ma LJ, Li W, Zhang X, et al. Differential gene expression profiling of laryngeal squamous cell carcinoma by laser capture microdissection and complementary DNA microarrays. Arch Med Res. 2009;40(2):114-123. doi:10.1016/j.arcmed.2008.12.005
-
Ni RS, Shen X, Qian X, Yu C, Wu H, Gao X. Detection of differentially expressed genes and association with clinicopathological features in laryngeal squamous cell carcinoma. Oncol Lett. 2012;4(6):1354-1360. doi:10.3892/ol.2012.920
-
Liu W, Jiao Y, Tian C, et al. Gene expression profiling studies using microarray in osteoarthritis: genes in common and different conditions. Arch Immunol Ther Exp (Warsz). 2020;68(5):28. doi:10.1007/s00005-020-00592-4
-
Jemal A, Murray T, Ward E, et al. Cancer statistics. CA Cancer J Clin. 2005;55(1):10-30. doi:10.3322/canjclin.55.1.10
-
Mondal S, Adhikari N, Banerjee S, Amin SA, Jha T. Matrix metalloproteinase-9 and its inhibitors in cancer: a minireview. Eur J Med Chem. 2020;194:112260. doi:10.1016/j.ejmech.2020.112260
-
Krecicki T, Fraczek M, Jelen M, Podhorska M, Szkudlarek T, Zatonski T. Expression of collagenase-1, collagenase-3 and tissue inhibitor of matrix metalloproteinase-1 in laryngeal squamous cell carcinomas. Eur Arch Otorhinolaryngol. 2003;260(9):494-497. doi:10.1007/s00405-003-0610-2
-
Liu W, Xiang C, Jia S. Study on the expression of matrix metalloproteinase-12 and NM_002426 in squamous laryngeal carcinoma. Lin Chuang Er Bi Yan Hou Ke Za Zhi. 2006;20(24):1120-1123.
-
Borgoño CA, Diamandis EP. The emerging roles of human tissue kallikreins in cancer. Nat Rev Cancer. 2004;4(11):876-890. doi:10.1038/nrc1474
-
Schrader CH, Kolb M, Zaoui K, et al. Kallikrein-related peptidase 6 regulates epithelial-to-mesenchymal transition and serves as prognostic biomarker for head and neck squamous cell carcinoma patients. Mol Cancer. 2015;14:107. doi:10.1186/s12943-015-0381-6
-
Stewart RL, O’Connor KL. Clinical significance of the integrin α6β4 in human malignancies. Lab Invest. 2015;95(9):976-986. doi:10.1038/labinvest.2015.82
-
Zhong F, Lu HP, Chen G, et al. The clinical significance and potential molecular mechanism of integrin subunit beta 4 in laryngeal squamous cell carcinoma. Pathol Res Pract. 2020;216(2):152785. doi:10.1016/j.prp.2019.152785
-
Hui L, Yang N, Yang H, Guo X, Jang X. Identification of biomarkers with a tumor stage-dependent expression and exploration of the mechanism involved in laryngeal squamous cell carcinoma. Oncol Rep. 2015;34(5):2627-2635. doi:10.3892/or.2015.4230
-
Tyszkiewicz T, Jarzab M, Szymczyk C, et al. Epidermal differentiation complex locus 1q21 gene expression in head and neck cancer and normal mucosa. Folia Histochem Cytobiol. 2014;52(2):79-89. doi:10.5603/fhc.2014.0018
-
Lauriola L, Michetti F, Maggiano N, et al. Prognostic significance of the Ca(2+) binding protein S100A2 in laryngeal squamous-cell carcinoma. Int J Cancer. 2000;89(4):345-349. doi:10.1002/1097-0215(20000720)89:43.0.co;2-t
-
Elazezy M, Schwentesius S, Stegat L, et al. Emerging insights into keratin 16 expression during metastatic progression of breast cancer. Cancers (Basel). 2021;13(15):3869. doi:10.3390/cancers13153869
-
Khanom R, Nguyen CT, Kayamori K, et al. Keratin 17 is induced in oral cancer and facilitates tumor growth. PLoS One. 2016;11(8):e0161163. doi:10.1371/journal.pone.0161163
-
Yang L, He Y, Zhang Z, Wang W. Upregulation of CEP55 predicts dismal prognosis in patients with liver cancer. Biomed Res Int. 2020;2020:4139320. doi:10.1155/2020/4139320
-
Hui L, Yang N, Yang H, Guo X, Jang X. Identification of biomarkers with a tumor stage-dependent expression and exploration of the mechanism involved in laryngeal squamous cell carcinoma. Oncol Rep. 2015;34(5):2627-2635. doi:10.3892/or.2015.4230
-
Madonna S, Girolomoni G, Dinarello CA, Albanesi C. The significance of IL-36 hyperactivation and IL-36R targeting in psoriasis. Int J Mol Sci. 2019;20(13):3318. doi:10.3390/ijms20133318
-
Molinolo AA, Hewitt SM, Amornphimoltham P, et al. Dissecting the Akt/mammalian target of rapamycin signaling network: emerging results from the head and neck cancer tissue array initiative. Clin Cancer Res. 2007;13(17):4964-4973. doi:10.1158/1078-0432.ccr-07-1041
-
Starska K, Forma E, Jóźwiak P, et al. Gene and protein expression of glucose transporter 1 and glucose transporter 3 in human laryngeal cancer: the relationship with regulatory hypoxia-inducible factor-1α expression, tumor invasiveness, and patient prognosis. Tumour Biol. 2015;36(4):2309-2321. doi:10.1007/s13277-014-2838-4
-
Zhou S, Wang S, Wu Q, Fan J, Wang Q. Expression of glucose transporter-1 and -3 in the head and neck carcinoma: the correlation of the expression with the biological behaviors. ORL J Otorhinolaryngol Relat Spec. 2008;70(3):189-194. doi:10.1159/000124293
-
Jia H, Wang J, Shan Y, Hou J. The clinical significance of the expression of p63 gene in laryngeal squamous cell carcinoma. Lin Chuang Er Bi Yan Hou Ke Za Zhi. 2006;20(3):119-121. Chinese.
-
Ikeda K, Iyama K, Ishikawa N, et al. Loss of expression of type IV collagen alpha5 and alpha6 chains in colorectal cancer associated with the hypermethylation of their promoter region. Am J Pathol. 2006;168(3):856-865. doi:10.2353/ajpath.2006.050384
-
Baba Y, Iyama K, Ikeda K, et al. The expression of type IV collagen alpha6 chain is related to the prognosis in patients with esophageal squamous cell carcinoma. Ann Surg Oncol. 2008;15(2):555-565. doi:10.1245/s10434-007-9592-4
-
Shi S, Zheng HC, Zhang ZG. Roles of fascin mRNA expression in colorectal cancer: meta-analysis and bioinformatics analysis. Mol Clin Oncol. 2020;13(2):119-128. doi:10.3892/mco.2020.2069
-
Chen SF, Yang SF, Li JW, et al. Expression of fascin in oral and oropharyngeal squamous cell carcinomas has prognostic significance: a tissue microarray study of 129 cases. Histopathology. 2007;51(2):173-183. doi:10.1111/j.1365-2559.2007.02755.x
-
Mahloogi H, González-Guerrico AM, Lopez De Cicco R, et al. Overexpression of the calcium sensor visinin-like protein-1 leads to a cAMP-mediated decrease of in vivo and in vitro growth and invasiveness of squamous cell carcinoma cells. Cancer Res. 2003;63(16):4997-5004.
-
Xie Y, Chan H, Fan J, et al. Involvement of visinin-like protein-1 in regulating proliferative and invasive properties of neuroblastoma. Carcinogenesis. 2007;28(10):2122-2130. doi:10.1093/carcin/bgm147
-
Haghnegahdar H, Du J, Wang D, et al. The tumorigenic and angiogenic effects of MGSA/GRO proteins in melanoma. J Leukoc Biol. 2000;67(1):53-62. doi:10.1002/jlb.67.1.53
-
Feng H, Liu X. Interaction between ACOT7 and LncRNA NMRAL2P via methylation regulates gastric cancer progression. Yonsei Med J. 2020;61(6):471-481. doi:10.3349/ymj.2020.61.6.471
-
Sumantran VN, Mishra P, Sudhakar N. Microarray analysis of differentially expressed genes regulating lipid metabolism during melanoma progression. Indian J Biochem Biophys. 2015;52(2):125-131.
-
Zhang X, Liu B, Zhang J, et al. Expression level of ACOT7 influences the prognosis in acute myeloid leukemia patients. Cancer Biomark. 2019;26(4):441-449. doi:10.3233/cbm-182287
-
Wang Q, Ma C, Kemmner W. Wdr66 is a novel marker for risk stratification and involved in epithelial-mesenchymal transition of esophageal squamous cell carcinoma. BMC Cancer. 2013;13:137. doi:10.1186/1471-2407-13-137
-
Zhang Y, Yu S, Zhang Z, Zhao G, Xu J. Long non-coding RNA AK096174 promotes cell proliferation and invasion in gastric cancer by regulating WDR66 expression. Biosci Rep. 2018;38(4):BSR20180277.
-
Akdi A, Giménez EM, García-Quispes W, et al. WDR3 gene haplotype is associated with thyroid cancer risk in a Spanish population. Thyroid. 2010;20(7):803-809. doi:10.1089/thy.2010.0072
-
Sato N, Koinuma J, Fujita M, et al. Activation of WD repeat and high-mobility group box DNA binding protein 1 in pulmonary and esophageal carcinogenesis. Clin Cancer Res. 2010;16(1):226-239. doi:10.1158/1078-0432.ccr-09-1405
-
Nguyen TM, Shin IW, Lee TJ, et al. Loss of ITM2A, a novel tumor suppressor of ovarian cancer through G2/M cell cycle arrest, is a poor prognostic factor of epithelial ovarian cancer. Gynecol Oncol. 2016;140(3):545-553. doi:10.1016/j.ygyno.2015.12.006
-
Ye H, Yu T, Temam S, et al. Transcriptomic dissection of tongue squamous cell carcinoma. BMC Genomics. 2008;9:69. doi:10.1186/1471-2164-9-69
Additional Information
Publisher’s Note
Bayrakol MP remains neutral with regard to jurisdictional and institutional claims.
Rights and Permissions
About This Article
How to Cite This Article
Emine Göktaş, Kayhan Öztürk, Ayşe Gül Zamani, Mahmut Selman Yıldırım. Determination of whole-genome expression differences in larynx cancers. Ann Clin Anal Med 2023;14(Suppl 1):S31-36. doi:10.4328/ACAM.21602
- Received:
- January 20, 2023
- Accepted:
- February 25, 2023
- Published Online:
- March 11, 2023
- Printed:
- March 25, 2023
