1. Introduction

Coronaviridae or coronaviruses are a large family of viruses that can cause disease in mammals and birds. In humans, the viruses cause illness, ranging from common cold to severe respiratory diseases. In 2002, severe acute respiratory syndrome (SARS) emerged, and in 2012, Middle East respiratory syndrome (MERS) emerged—both are betacoronoviruses transmitted from animal to human that result in severe respiratory diseases in affected individuals [1]. In December 2019, a novel SARS-coronavirus (CoV)-2, also known as 2019 novel corona virus (2019-nCoV) or coronavirus disease-2019 (COVID-19), caused a pneumonia outbreak in China and subsequently expanded worldwide, leading to a pandemic. Recently, the first genomic sequence of COVID-19 was released, and through comparing to the genomes of SARS-CoV and MERS-CoV, researchers found that COVID-19 has better genomic sequence homology with SARS-CoV than that of MERS-CoV [2,3].

During the first SARS-CoV outbreak, Li et al. [4] identified angiotensin-converting enzyme 2 (ACE2) as the human host factor or cell entry receptor for SARS-CoV. Overexpression of ACE2 and injection of SARS-CoV spike protein developed severe acute lung failure in mice, which can be attenuated by blocking the renin-angiotensin pathway [5]. Recent studies have revealed that COVID-19 spike protein has strong affinity with ACE2 on host cells, which is significantly higher than that of SARS-CoV [1,6]. These studies also pointed out that treatment with transmembrane protease serine 2 (TMPRSS2) inhibitor significantly blocked SARS-CoV cell entry; therefore, either ACE2 or TMPRSS2 blockers can be a potential targets for anti-viral intervention. Another study reported that the COVID-19 receptor binding domain was capable of entering cells expressing human ACE2, while other receptors are ineffective, confirming that human ACE2 is the prime receptor for COVID-19 [7]. Since the host cell receptor plays a crucial role in virus entry, targeting the precise receptor ACE2 is a promising preventive strategy for COVID-19 infection. Recent studies have demonstrated that ACE2 overexpression was frequently observed in gastrointestinal tissues and colon cell lines [1], which is comparatively higher than that of other tissues, including lung tissues. Therefore, in this study, HT-29, a colon adenocarcinoma cell line, was employed to investigate the ACE2 inhibitory effect of test samples in vitro.

At present, there is no definite treatment or vaccine developed for the coronavirus that causes COVID-19. Nevertheless, many possible treatments for COVID-19 have been thrust into the spotlight by scientists and health industries. For example, the anti-malarial drug combination of chloroquine and hydroxychloroquine as well as anti-HIV drugs ritonavir and lopinavir have been recommended [8,9]. Since the drugs directly target the pathogen, the effectiveness of these drugs are largely anecdotal. Additionally, the development of new drugs for targeting ACE2 and treating COVID-19 could be time-consuming. Hence, the safety efficacy of the new drugs are a prime concern, which requires a long time for testing, while the infection is growing fast. The traditional medicine systems from many geographical areas use herbs as the primary treatment of viral infections, including those caused by SARS-CoV. For example, leaf extracts of Toona sinensis inhibit SARS-CoV replication [10]. Licorice has been suggested as a promising treatment for SARS-CoV [11]. In addition, natural products including diterpenoids, sesquiterpenoids, triterpenoids, lignoids, curcumin, and ginsenosside-Rb1 have been shown to inhibit SARS-CoV [12,13]. A recent in silico study showed that baicalin, scutellarin, hesperetin, nicotianamine, and glycyrrhizin were capable of inhibiting ACE2 [3].

There are a wide-range of essential oils, and their components have been clinically proven to possess antiviral properties [14,15]. A study by Jackwood et al. [16] found that treatment with a mixture of oleoresins and essential oils from botanicals decreased the severity of clinical signs and lesions in chickens that carried the avian infectious bronchitis virus (IBV-CoV). However, the effects of plant essential oils on human coronaviruses are yet to be explored. In particular, their effects on host cell receptors have barely been investigated. Geranium essential oil is derived from the leaves of Pelargonium graveolens, which is widely utilized in the cosmetic industry, perfumery, and aromatherapy [17]. Traditionally, geranium essential oil has been used in both physiological and psychological complications, including anxiety, insomnia, high blood pressure, worry, anger, frustration, restlessness, nervousness, weight loss, hypercholesterolemia, gastrointestinal disorders, and respiratory tract infection [18,19]. Geranium oil represents a potent immune modulator, and stimulates and cleans lymphatic system [20]. The oil is clinically used for treating diarrhea, jaundice, diabetes, hepatitis, ulcerative colitis, cholecystitis, and renal stone [21]. The present study aimed to investigate the ACE2 inhibitory effects of plant essential oils and their major components in vitro.

2. Results

2.1. Cytotoxic Effect of Essential Oils on HT-29 Cells

Prior to investigation of the ACE2 inhibitory effect of essential oils, we determined the cytotoxicity of essential oils on HT-29 cells by 3-(4,5-dimethyl-thiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) assay. As shown in Table 1, we treated HT-29 cells with increasing concentrations of essential oils (25–200 mg/mL) for 48 h. We found that geranium, cypress, eucalyptus, juniper berry, marjoram, and tea tree oils were not cytotoxic up to a concentration of 200 mg/mL, whereas fennel, frankincense, myrtle, neroli, petitgrain, and tangerine oils displayed moderate cytotoxicity to HT-29 cells, with IC₅₀ values of 143.8, 155.91, 144.22, 119.52, 151.63, and 109.31 mg/mL, respectively. Treatment with citronella, may chang, kunzea, palmarosa, black pepper, and elemi oils exhibited strong cytotoxicity towards HT-29 cells, with IC₅₀ values of 0.07, 0.11, 0.44, 1.23, 7.62, and 7.93 mg/mL. The histogram analysis of individual essential oils on cell viability is shown in Figure S1. On the basis of the MTT assay, we selected the optimum non-cytotoxic concentration of essential oils for our further ACE2 assay.

Table 1. Cytotoxic and ACE2 inhibitory effects of essential oils.

Values represent the mean ± SD of three independent experiments. p-values of less than 0.05*, 0.01**, 0.001***, and 0.0001**** were considered statistically significant for the sample treatment group vs. the control group

2.2. Essential Oils Downregulate ACE2 Activity in HT-29 Cells

Next, we examined the effect of selected essential oils on ACE2 activity in HT-29 cells. All the treated essential oils exhibited significant inhibition on ACE2 activity (Table 1). Among them, geranium and lemon essential oils strongly reduced the human ACE2 levels in HT-29 cells from 17.68 ng/mL (control) to 1.43 ng/mL and 4.34 ng/mL, respectively (Figure 1A). In addition, we found that geranium and lemon essential oils dose-dependently reduced ACE2 activity in HT-29 cells (Figure 1B,C). Further, to confirm this effect, we examined the protein expression levels of ACE2 by immunoblotting. Similar to that of the ACE2 ELISA assay, the ACE2 protein level was significantly reduced by geranium and lemon essential oils (Figure 2A). Further qPCR analysis also confirmed that geranium and lemon essential oils significantly downregulated ACE2 (Figure 2B) and TMPRSS2 (Figure 2C) mRNA levels in HT-29 cells.

Figure 1. Essential oils reduced human angiotensin-converting enzyme 2 (ACE2) levels in HT-29 cells. (A) The ACE2 inhibitory effects of selected essential oils. HT-29 cells were incubated with indicated concentration of selected essential oils for 48 h. Cell lysates were subjected to determine ACE2 levels using a commercially available ELISA kit. (B,C) The dose-dependent inhibitory effects of geranium and lemon essential oils were determined in the same course of time. Values represent the mean ± SD of three independent experiments. p-values of less than 0.05*, 0.01**, 0.001***, and 0.0001**** were considered statistically significant for the sample treatment group vs. the control group.

Figure 2. Geranium and lemon essential oils downregulated ACE2 expression in HT-29 cells. (A) Protein expression levels of ACE2 were determined by immunoblotting. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as an internal loading control. Relative density of one representative experiment is shown, where ACE2 signal was normalized to GAPDH signal. (B,C) Relative expression of ACE2 and TMPRSS2 mRNAs in HT-29 cells. Total RNA was extracted from cells treated with indicated concentration of geranium and lemon essential oils for 48 h. The transcription levels of ACE2 and TMPRSS2 were quantified by qPCR and a representative experiment is shown. The D^ct values of ACE2 and TMPRSS2 mRNAs were normalized to GAPDH mRNA. Values represent the mean ± SD of three independent experiments. p-values of less than 0.05*, 0.01**, and 0.001*** were considered statistically significant for the sample treatment group vs. the control group.

2.3. Chemical Compositions of Geranium and Lemon Essential Oils

The major chemical constituents of geranium and lemon oils and their relative amounts were determined by GC–MS analysis. The chemical composition of tested essential oils are listed in Table S1–S8). The GC–MS profiles and the major compounds of geranium and lemon oils are shown in Figure 3A,B. The relative contents (%) in geranium and lemon essential oils are shown in Tables 2 and 3, respectively. A total of 22 compounds were identified in geranium oil, accounting for 90.27% of the whole oil. The major component in geranium oils were citronellol (27.1%), geraniol (21.4%), and neryl acetate (10.5 %), which made up around 60% of the content of the whole geranium oil (Table 2). Lemon oil contains nine identified compounds, accounting for 98.1% of the whole oil. The main component of lemon oil was found to be limonene, which represents 73.0% of the whole oil. The second and third most abundant compounds were γ-terpinene (9.2%) and β-pinene (8.6%) (Table 3).

2.4. ACE2 Inhibitory Effects of Major Constituents in Geranium and Lemon Essential Oils

The cytotoxic effect of citronellol, geraniol, limonene, and neryl acetate were determined by MTT assay. Treatment with either citronellol, geraniol, limonene, or neryl acetate did not display cytotoxicity in HT-29 cells up to a concentration of 100 mM for 48 h (Figure S2). Next, we examined the ACE2 inhibitory effects of these compounds. Interestingly, the ACE2 activity was significantly inhibited by citronellol (50 mM), geraniol (50 mM), limonene (50 mM), and neryl acetate (50 mM), as these compound reduced the ACE2 levels from 18.0 ng/mL (control) to 7.67 ng/mL, 10.44 ng/mL, 12.92 ng/mL, and 16.63 ng/mL, respectively (Figure 4A). Similar to the ELISA assay, the immunobloting also confirmed that treatment with citronellol, geraniol, limonene, and neryl acetate significantly inhibited ACE2 protein expression in HT-29 cells (Figure 4B). qPCR analysis further supported this notion that citronellol, geraniol, and limonene significantly downregulate ACE2 mRNA levels, whereas the inhibition of ACE2 mRNA by neryl acetate was not statistically significant (Figure 4C). Interestingly, treatment with major compounds significantly downregulated the TMPRSS2 mRNA levels (Figure 4D), which were involved in S-protein priming during SARS-CoV-2 cell entry.

3. Discussion

Essential oils have been used in folk medicine throughout history as principal ingredients in aromatherapy and psychotherapy. However, modern preclinical studies have revealed that essential oils possess a broad spectrum of pharmacological activities, including antiseptic, diuretic, choleretic, spasmolytic, hyperemic, expectorant, anti-anxiety, antioxidant, anti-inflammatory, and anticancer activities [22]. The pharmacological properties of essential oils are attributed by their unique major constituents, such as monoterpenoids, sesquiterpenoids, and phenylpropanoids. In addition, there is considerable evidence emerging from in vitro studies and controlled trails of the potential for essential oils as antiviral agents for the treatment of human viral infections, including SARS coronaviruses [23]. These antiviral essential oils were mostly tested against enveloped RNA or DNA viruses such as herpes simplex virus type 1, Junin virus, influenza virus, dengue virus type 2, and SARS coronaviruses, as well as non-enveloped viruses such as coxsackie virus B1, poliovirus, and adenovirus type 3 [22]. Most of these clinically useful antiviral agents are substances that act on specific steps of the viral biosynthesis, inhibiting viral replication in particular. On the other hand, virucidal agents denature the structure or glycoproteins of the virus, thus reducing or completely blocking the infectivity of virus particles [24].

Figure 3. GC–MS analysis of geranium and lemon essential oils. GC profiles and chemical structures of the major compounds in geranium (A) and lemon (B) essential oils.

Recently, ACE2 was identified as a functional SARS coronavirus receptor in epithelial cells. In particular, the SARS-CoV spike protein engages ACE2 as the cell entry receptor and employs the cellular serine protease TMPRSS2 for spike protein priming [1]. Thus, either inhibition of ACE2 or TMPRSS2 receptors can be potential targets for SARS-CoV prevention, as ACE2 displays distinct modes of enzyme action and tissue distribution. Abundant endogenous expression of ACE2 has been observed in cardiovascular and colon-specific cell lines, while fainted expression has been noted in kidney cells [25]. Recent studies have also confirmed that the endogenous expression of ACE2 was higher in gastrointestinal tissues [26] and colon cell lines [1], with these reports pointing out that ACE2 expression is comparatively higher in intestinal tissues than in lung tissues. Therefore, in this study, we employed HT-29, a colon adenocarcinoma cell line, which endogenously expresses ACE2, in order to investigate the ACE2 inhibitory effect of essential oils. Initially, we examined the cytotoxic effects of essential oils in HT-29 cells. Assessment of cytotoxicity is clearly an important aspect of the evaluation of a potent antiviral agent, as a potential agent should be selective for virus- or cell-specific processes with no or limited effects on cellular metabolism [27]. Our results indicated that each essential oil exhibits differential cytotoxicity, some of them are strongly cytotoxic and others of them do not exhibit cytotoxicity up to a concentration of 200 mg/mL.

 Figure 4. ACE2 inhibitory effects of major compounds of geranium and lemon essential oils. (A) ACE2 activity in HT-29 cells were determined by ELISA assay. (B) Protein expression levels of ACE2 were determined by immunoblotting. Relative density of one representative experiment is shown, where ACE2 signal was normalized to GAPDH signal. (C,D) Relative expression of ACE2 and TMPRSS2 mRNAs in HT-29 cells. Values represent the mean ± SD of three independent experiments. p-values of less than 0.05*, 0.01**, 0.001***, and 0.0001**** were considered statistically significant for the sample treatment group vs. the control group.

Table 2. The main components and their relative contents (%) of geranium oil.

RT: retention time; * MS: National Institute of Standards and Technology (NIST) and Wiley libraries and the literature; KI: Kovats index on a DB-5MS column in reference to n-alkanes; ST: authentic standard compounds.

Table 3. Main components and their relative contents (%) of lemon oil.

RT: retention time; * MS: NIST and Wiley libraries and the literature; KI: Kovats index on a DB-5MS column in reference to n-alkanes; ST: authentic standard compounds.

Recent in silico studies reported that organosulfur compounds in garlic essential oil and other natural products, such as baicalin, scutellarin, hesperetin, nicotianamine, glycyrrhizin, (E,E)-α-farnesene, (E)-β-farnesene, and (E,E)−farnesol, have the potential to bind the human ACE2 receptor, thereby possibly blocking SARS-CoV-2 cell entry [3,28,29]. One study evaluated the in vitro antiviral effects against SARS-CoV of seven essential oils from Lebanese species that included Laurus nobilis T., Juniperus oxycedrus spp. oxycedrus, Thuja orientalis, Cupressus sempervirens spp. pyramidalis, Pistacia palaestina, Salvia officinalis, and Satureja thymbra. Among them, the fruit essential oils of L. nobilis, T. orientalis, and J. oxycedrus spp. oxycedrus displayed inhibition against SARS-CoV-induced cytopathogenic effect with IC₅₀ values of 120 mg/mL, 130 mg/mL, and 270 mg/mL, respectively [27]. However, there are no studies directly exhibiting the ACE2 inhibitory effects of essential oils or their major compounds in vitro or in vivo. In this study, we screened the ACE2 inhibitory effect of 10 essential oils; among them, 8 essential oils exhibited significant ACE2 inhibition. Interestingly, geranium and lemon essential oils strongly inhibited ACE2 activity without displaying cytotoxicity. To the best of our knowledge, this is the first report indicating that geranium and lemon essential oils and their major components citronellol, geraniol, limonene, linalool, and neryl acetate downregulate ACE2 receptor activity in virus–host epithelial cells.

4. Materials and Methods

4.1. Chemicals and Reagents

All the essential oils were provided by Bio-Jourdeness International Groups Co. Ltd. (Taichung Taiwan). Citronellol, geraniol, neryl acetate, and limonene were obtained from Tokyo Chemical Industry Co., Ltd. (Tokya, Japan). Purity of these compounds were above 99%, according to the GC and ¹H-NMR analysis. McCoy’s medium, fetal bovine serum (FBS), sodium pyruvate, penicillin, and streptomycin were obtained from Invitrogen (Carlsbad, CA, USA). 3-(4,5-Dimethyl-thiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) and dimethylsulfoxide (DMSO) were purchased from Sigma-Aldrich. Antibody against ACE2 was obtained from Arigo Biolaboratories (Hsinchu, Taiwan). Antibody against GAPDH was obtained from Cell Signaling Technology Inc. All other chemicals were reagent grade or HPLC grade and were supplied by either Merck or Sigma-Aldrich.

4.2. Cell Culture and Cell Viability Assay

Human colorectal adenocarcinoma cell line (HT-29) was cultured and maintained in McCoy’s medium supplemented with 10% fetal bovine serum, 1% L-glutamine, and 1% penicillin/streptomycin. Cell viability was assessed by MTT colorimetric assay. Briefly, HT-29 cells were seeded into a 48-well culture plate with a density of 5 × 10⁴ cells per well. After 24 h incubation, we treated cells with increasing concentrations of essential oils for 48 h. The culture supernatant was removed and 1 mg/mL of MTT in 0.1 mL fresh medium was added. The MTT formazon crystals were dissolved in 0.4 mL of DMSO and the samples were measured at 570 nm (A₅₇₀) using ELISA micro-plate reader (Bio-Tek Instruments, Winooski, VT). The percentage of cell viability (%) was calculated as (A₅₇₀ of treated cells/A₅₇₀ of untreated cells) ´ 100.

4.3. Determination of ACE2 Activity

Cellular ACE2 activity was measured using a commercially available human ACE2 ELISA kit (Elabscience Biotechnology Inc, Hubei, China). Briefly, HT-29 cells were seeded into a 10 cm cell culture dish with a density of 2 × 10⁶ cells per well. After 48 h, we treated cells with a selective dose (non-cytotoxic high concentration) of essential oils or derived compounds for a further 48 h. Then, the cells were gently washed with pre-cooled phosphate buffered saline (PBS) and we dissociated the cells using trypsin. Cell suspension was centrifuged for 5 min at 1000× g. Then, the cells were lysed with a mixture of radioimmuno precipitation assay (RIPA) buffer (Thermo Fisher Scientific, Waltham, MA) and transmembrane protein extraction reagent (Fivephoton Biochemical, San Diego, CA, USA) at a ratio of 1:1, followed by centrifugation for 10 min at 1500× g at 4°C. Protein concentration was determined by Bio-Rad protein assay reagent (Bio-Rad Laboratories, Hercules, CA, USA). Equal amounts of protein samples (100 µg/well) were assayed according to the manufacturer’s protocol. The optical density (OD value) of each sample was determined at 450 nm (A₄₅₀) using an ELISA micro-plate reader (Bio-Tek Instruments).The percentage of ACE2 activity (%) was calculated as (A₄₅₀ of treated cells/A₄₅₀ of untreated cells) ´ 100.

4.4. Determination of ACE2 Protein

HT-29 cells were seeded in a 10 cm dish at a density of 2 ´ 10⁶ cells per dish. After 24 h, we exposed the cells to test samples for 48 h. Cellular and membrane-bound proteins were isolated by a mixture of RIPA lysis buffer and transmembrane protein extraction reagent. Protein concentration was determined by Bio-Rad protein assay reagent. Equal amounts of denatured protein samples (60 µg) were separated by 10% SDS-PAGE, and the separated proteins were transferred onto polyvinylidene difluoride (PVDF) membrane overnight. Then, the membranes were blocked with 5% non-fat dried milk for 30 min, followed by incubation with ACE2 or GAPDH antibodies overnight, and then incubated with either horseradish peroxidase-conjugated goat anti-rabbit or anti-mouse antibodies for 1 h. Immunoblots were developed with the enhanced chemiluminescence (ECL) reagents (Advansta Inc., San Jose, CA, USA), images were captured by CheniDoc XRS⁺ docking system, and the protein bands were quantified by using Imagelab software (Bio-Rad laboratories, Hercules, CA, USA). 

4.5. Quantitative Real-Time PCR

Total RNA was extracted by total RNA purification kit (GeneMark, New Taipei City, Taiwan). RNA concentration was quantified with a NanoVue Plus spectrophotometer (GE Health Care Life Sciences, Chicago, IL, USA). Quantitative PCR (qPCR) was performed on a real-time PCR detection system and software (Applied Biosystems, Foster City, CA, USA). First-strand complementary DNA (cDNA) was generated by SuperScript III reverse transcriptase kit (Invitrogen). Quantification of mRNA expression for genes of interest was performed by qPCR reactions with equal volume of cDNA, forward and reverse primers (10 µM), and power SYBR Green Master Mix (Applied Biosystems, Foster city, CA, USA). The sequence of the PCR primers were as follows: ACE2: forward 5¢-GCTGCTCAGTCCACCATTGAG-3¢, reverse 5¢-GCTTCGTGGTTAAACTTGTCCAA-3¢; TMPRSS2: forward 5¢-AATCGGTGTGTTCGCCTCTAC-3¢, reverse 5¢-GCGGCTGTCACGATCC-3¢; GAPDH: forward 5'-TCCTGGTATGACAACGAAT-3', reverse 5'-GGTCTCTCTCTTCCTCTTG-3' [30]. The copy number of each transcript was calculated as the relative copy number normalized by GAPDH copy number.

4.6. GC–MS Analysis

To determine the essential oil composition, we carried out analyses using an ITQ 900 mass spectrometer coupled with a DB-5MS column. The temperature program was as follows: 45 °C for 3 min, then increased to 3 °C/min to 180 °C, and then increased to 10 °C/min to 280 °C hold for 5 min. The other parameters were injection temperature, 240 °C; ion source temperature, 200 °C; EI, 70 eV; carrier gas, He 1 mL/min; and mass scan range, 40–600 m/z. The volatile compounds were identified by Wiley/NBS Registry of mass spectral databases (V. 8.0, Hoboken, NJ, USA), National Institute of Standards and Technology (NIST) Ver. 2.0 GC/MS libraries and the Kovats indices were calculated for all volatile constituents using a homologous series of n-alkanes C₉-C₂₄. The major components were identified by co-injection with standards (wherever possible).

4.7. Statistical Analysis

Data are expressed as mean ± SD. All data were analyzed using the statistical software Graphpad Prism version 6.0 for Windows (GraphPad Software, San Diego, CA, USA). Statistical analysis was performed using one-way ANOVA followed by Dunnett’s test for multiple comparison. p-values of less than 0.05*, 0.01**, 0.001***, and 0.0001**** were considered statistically significant for the sample treatment group vs. the control group.

5. Conclusions

The recent COVID-19/SARS-CoV-2 pneumonia pandemic constitutes the largest global public health crisis and has created international anxiety due to its relatively high infectious, rapid progression, as well as its relatively high death rate. Drug development for treating COVID-19 is currently important due to its rapid progression. The conventional antiviral vaccine development is time-consuming and its safety needs to be verified. The fact is that there is no conventional medicine developed for the treatment of COVID-19. However, there is evidence that suggests that essential oils and their major components have displayed potent antiviral activity to other coronaviruses, such as SARS-CoV-1, although the mechanism of action of these oils and their components were found to be mainly through inhibition of viral replication [27]. Recently, ACE2, a receptor in host cell favoring virus cell entry, was recognized as one of the prime targets to minimize the infection. In this study, we presented the first piece of evidence that geranium and lemon essential oils and their major compounds, citronellol, geraniol, limonene, linalool, and neryl acetate, could downregulate ACE2 expression in epithelial cells, thereby blocking virus entry into host cells, and eventually preventing viral infection. However, further studies are highly warranted to unveil the underlying molecular mechanisms of this inhibitory effect.

Supplementary Materials: The following are available online www.mdpi.com/xxx/s1. Figure S1: Cytotoxic effect of essential oils on HT-29 cells. Figure S2: Cytotoxic effect of major compounds of geranium and lemon essential oils. Table S1. Main components and their relative contents (%) of Bergamot oil derived from citrus rind (peel) of Citrus bergamia. Table S2. Main components and their relative contents (%) of cypress oil derived from needles of Cupressus sempervirens. Table S3. Main components and their relative contents (%) of Eucalyptus oil derived from leaves of Eucalyptus globulus. Table S4. Main components and their relative contents (%) of juniper berry oil derived from berries of Juniperus communis. Table S5. Main components and their relative contents (%) of neroli oil derived from flowers of Citrus aurantium. Table S6. Main components and their relative contents (%) of petitgrain oil derived from leaves and twigs of Citrus aurantium. Table S7. Main components and their relative contents (%) of tangerine oil derived from fruit peel of Citrus reticulata. Table S8. Main components and their relative contents (%) of tea tree oil derived from leaves of Mellaleuca alternifolia.

Author ContributionsConceptualization, S.-Y.W. and K.J.S.K.; methodology and investigation, K.J.S.K. and M.G.V., C.-S.L; formal analysis, K.J.S.K., C.-S.L., C.-S.W., L.-P.L., and S.-S.H.; resources, C.-C.C., C.-H.H., and Y.-C.C.; writing—original draft preparation, K.J.S.K. and S.-Y.W.; supervision, S.-Y.W.

Funding: This research received no external funding.

Conflicts of Interest: The authors declare that there is no conflict of interest. The company provides no financial support for the research of this article. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

1. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Kruger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 2020, 181, 271–280e.
2. Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N.; et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding. Lancet 2020, 395, 565–574.
3. Chen, H.; Du, Q. Potential natural compounds for preventing 2019-nCoV infection. Preprints 2020, 2020010358.
4. Li, W.; Moore, M.J.; Vasilieva, N.; Sui, J.; Wong, S.K.; Berne, M.A.; Somasundaran, M.; Sullivan, J.L.; Luzuriaga, K.; Greenough, T.C.; et al. Angiotensin-converting enzyme 2 is a functional receptor for the sars coronavirus. Nature 2003, 426, 450–454.
5. Kuba, K.; Imai, Y.; Rao, S.; Gao, H.; Guo, F.; Guan, B.; Huan, Y.; Yang, P.; Zhang, Y.; Deng, W.; et al. A crucial role of angiotensin converting enzyme 2 (ACE2) in sars coronavirus-induced lung injury. Med. 2005, 11, 875–879.
6. Ziegler, C.G.K.; Allon, S.J.; Nyquist, S.K.; Mbano, I.M.; Miao, V.N.; Tzouanas, C.N.; Cao, Y.; Yousif, A.S.; Bals, J.; Hauser, B.M.; et al. SARS-CoV-2 receptor ACE2 is an interferon-stimulated gene in human airway epithelial cells and is detected in specific cell subsets across tissues. Cell 2020, 181, 1016–1035.e19.
7. Letko, M.; Munster, V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Microbiol. 2020, 5, 562–569.
8. Pastick, K.A.; Okafor, E.C.; Wang, F.; Lofgren, S.M.; Skipper, C.P.; Nicol, M.R.; Pullen, M.F.; Rajasingham, R.; McDonald, E.G.; Lee, T.C.; et al. Review: Hydroxychloroquine and chloroquine for treatment of sars-cov-2 (covid-19). Open Forum Infect. Dis. 2020, 7, ofaa130.
9. Sanders, J.M.; Monogue, M.L.; Jodlowski, T.Z.; Cutrell, J.B. Pharmacologic treatments for coronavirus disease 2019 (COVID-19): A review. JAMA 2020, 323, 1824–1836.
10. Chen, C.-J.; Michaelis, M.; Hsu, H.-K.; Tsai, C.-C.; Yang, K.D.; Wu, Y.-C.; Cinatl, J.; Doerr, H.W. Toona sinensis roem tender leaf extract inhibits sars coronavirus replication. Ethnopharmacol. 2008, 120, 108–111.
11. Pilcher, H. Liquorice may tackle SARS. Nature 2003, doi:10.1038/news030609-16.
12. Wen, C.C.; Kuo, Y.H.; Jan, J.T.; Liang, P.H.; Wang, S.Y.; Liu, H.G.; Lee, C.K.; Chang, S.T.; Kuo, C.J.; Lee, S.S.; et al. Specific plant terpenoids and lignoids possess potent antiviral activities against severe acute respiratory syndrome coronavirus. Med. Chem. 2007, 50, 4087–4095.
13. Wu, C.Y.; Jan, J.T.; Ma, S.H.; Kuo, C.J.; Juan, H.F.; Cheng, Y.S.; Hsu, H.H.; Huang, H.C.; Wu, D.; Brik, A.; et al. Small molecules targeting severe acute respiratory syndrome human coronavirus. Natl. Acad. Sci. USA 2004, 101, 10012–10017.
14. Astani, A.; Reichling, J.; Schnitzler, P. Comparative study on the antiviral activity of selected monoterpenes derived from essential oils. Res. 2010, 24, 673–679.
15. Swamy, M.K.; Akhtar, M.S.; Sinniah, U.R. Antimicrobial properties of plant essential oils against human pathogens and their mode of action: An updated review. Based Complement. Alternat. Med. 2016, 2016, 3012462–3012462.
16. Jackwood, M.W.; Rosenbloom, R.; Petteruti, M.; Hilt, D.A.; McCall, A.W.; Williams, S.M. Avian coronavirus infectious bronchitis virus susceptibility to botanical oleoresins and essential oils in vitro and in vivo. Virus Res. 2010, 149, 86–94.
17. Fayed, S.A. Antioxidant and anticancer activities of citrus reticulate (petitgrain mandarin) and pelargonium graveolens (geranium) essential oils. J. Agric. Biol. Sci. 2009, 5, 740–747.
18. Lis-Balchin, M. A chemotaxonomic study of the pelargonium (geraniaceae) species and their modern cultivars. Hort. Sci. 1997, 72, 791–795.
19. Asgarpanah, J.; Ramezanloo, F. An overview on phytopharmacology of pelargonium graveolens l. J. Trad. Knowl. 2015, 14, 558–563.
20. Standen, M.D.; Connellan, P.A.; Leach, D.N. Natural killer cell activity and lymphocyte activation: Investigating the effects of a selection of essential oils and components in vitro. J. Aromather. 2006, 16, 133–139.
21. Peterson, A.; Machmudah, S.; Roy, B.C.; Goto, M.; Sasaki, M.; Hirose, T. Extraction of essential oil from geranium (Pelargonium graveolens) with supercritical carbon dioxide. Chem. Technol. Biotechnol. 2006, 81, 167–172.
22. Reichling, J.; Schnitzler, P.; Suschke, U.; Saller, R. Essential oils of aromatic plants with antibacterial, antifungal, antiviral, and cytotoxic properties—An overview. Med. Res. 2009, 16, 79–90.
23. Nadjib, B.M. Effective antiviral activity of essential oils and their characteristic terpenes against coronaviruses: An update. Pharmacol. Clin. Toxicol. 2020, 8, 1138.
24. Zhu, J.D.; Meng, W.; Wang, X.J.; Wang, H.C.R. Broad-spectrum antiviral agents. Microbiol. 2015, 6, 517–517.
25. Warner, F.J.; Lew, R.A.; Smith, A.I.; Lambert, D.W.; Hooper, N.M.; Turner, A.J. Angiotensin-converting enzyme 2 (ACE2), but not ace, is preferentially localized to the apical surface of polarized kidney cells. Biol. Chem. 2005, 280, 39353–39362.
26. Xu, H.; Zhong, L.; Deng, J.; Peng, J.; Dan, H.; Zeng, X.; Li, T.; Chen, Q. High expression of ace2 receptor of 2019-ncov on the epithelial cells of oral mucosa. J. Oral Sci. 2020, 12, 8.
27. Loizzo, M.R.; Saab, A.M.; Tundis, R.; Statti, G.A.; Menichini, F.; Lampronti, I.; Gambari, R.; Cinatl, J.; Doerr, H.W. Phytochemical analysis and in vitro antiviral activities of the essential oils of seven lebanon species. Biodivers. 2008, 5, 461–470.
28. Thuy, B.T.P.; My, T.T.A.; Hai, N.T.T.; Hieu, L.T.; Hoa, T.T.; Thi Phuong Loan, H.; Triet, N.T.; Anh, T.T.V.; Quy, P.T.; Tat, P.V.; et al. Investigation into SARS-CoV-2 resistance of compounds in garlic essential oil. ACS Omega 2020, 5, 8312–8320.
29. da Silva, J.K.R.; Figueiredo, P.L.B.; Byler, K.G.; Setzer, W.N. Essential oils as antiviral agents, potential of essential oils to treat sars-cov-2 infection: An in-silico investigation. J. Mol. Sci. 2020, 21, 3426.
30. Bernardi, S.; Zennaro, C.; Palmisano, S.; Velkoska, E.; Sabato, N.; Toffoli, B.; Giacomel, G.; Buri, L.; Zanconati, F.; Bellini, G.; et al. Characterization and significance of ace2 and mas receptor in human colon adenocarcinoma. Renin Angiotensin Aldosterone Syst. 2012, 13, 202–209.