Abbreviations
-
- ACE2
-
- angiotensin‐converting enzyme 2
-
- ADAM
-
- a disintegrin and metallopeptidase
-
- AEC
-
- airway epithelial cell
-
- ANOVA
-
- analysis of variance
-
- ARDS
-
- acute respiratory distress syndrome
-
- BDP
-
- beclomethasone dipropionate
-
- COPD
-
- chronic obstructive pulmonary disease
-
- COVID
-
- coronavirus disease
-
- CTL
-
- cathepsin‐L
-
- FEV1
-
- forced expiratory volume in 1 s
-
- GINA
-
- Global initiative for Asthma
-
- GOLD
-
- Global Initiative for Chronic Obstructive Lung Disease
-
- HC
-
- healthy control
-
- HSPA
-
- heat shock protein
-
- ICS
-
- inhaled corticosteroid
-
- IL
-
- interleukin
-
- PCR
-
- polymerase chain reaction
-
- PI4KB
-
- phosphatidylinositol 4‐kinase beta
-
- qPCR
-
- quantitative PCR
-
- SARS‐CoV
-
- severe acute respiratory syndrome‐associated coronavirus
-
- TMPRSS
-
- transmembrane serine protease
INTRODUCTION
In December 2019, an outbreak of a novel respiratory illness resulting in severe disease and respiratory failure emerged in Wuhan, China, and has since spread to become the major pandemic in living memory, now known as coronavirus disease (COVID)‐19.1 The virus responsible was rapidly identified as a previously unknown beta‐coronavirus, closely related to severe acute respiratory syndrome‐associated coronavirus (SARS‐CoV) and is now identified as SARS‐CoV‐2.2, 3 The virus binds to the angiotensin‐converting enzyme 2 (ACE2) receptor and requires the transmembrane serine protease 2 (TMPRSS2) to cleave the viral spike protein in order to enter a cell.3 This step appears to be facilitated by endosomal proteases such as cathepsin‐L (CTL) and enhanced by the protein furin4, 5; then, the virus enters the host cell by endocytosis.
COVID‐19 has had a devastating impact, now infecting more than 30 million and is responsible for at least 100 000 000 deaths worldwide (World Health Organization (WHO) reports).6 Infection begins in the upper respiratory tract and while the majority of people experience a mild disease course, 10–15% have more severe disease with progression to infection of the lower airways, resulting in pneumonia, acute respiratory distress syndrome (ARDS) and death.7, 8 People with chronic respiratory disease, especially asthma and chronic obstructive pulmonary disease (COPD), are usually at heightened risk of complications from acute respiratory viral infections.9 This is, however, not clearly the case with COVID‐19, while there appears to be a heightened risk for active smokers,10, 11 those with asthma do not appear to be over‐represented. Age and male sex appear to be associated with worse outcomes, and children generally experience only mild illness.12–14 One factor that may predispose to more severe disease and pneumonia is the level of expression of the ACE2 receptor in airway epithelial cells (AEC) of the lower respiratory tract allowing infection to spread more easily from the upper airway, resulting in a more intense lower airway immune response and more severe clinical disease.
We sought to determine the expression of ACE2 and other genes that are crucial for viral entry in primary human AEC in a cohort of healthy individuals and in people with asthma and COPD.
METHODS
Participants
Subjects for this study had all been previously and independently recruited to provide lower airway samples for studies of chronic asthma or COPD and comprise two independent cohorts based in Newcastle and Perth, Australia (for more details, see Appendix S1 in Supplementary Information). The study was approved by the Hunter New England Human Research Ethics Committee (Reference No. 05/08/10/3.09), Perth Children’s Hospital Ethics Committee (RGS1470) and St John of Gods Human Ethics Committee (SJOG#901), and all participants provided written informed consent.
Gene expression by qPCR
Detailed description of the quantitative PCR (qPCR) is provided in Table S1 (Supplementary Information). Results were calculated using 2−ΔΔCt (where Ct is the threshold cycle) relative to the mean ΔCt of the healthy control group as described previously.15
Immunohistochemistry and immunofluorescence
Formalin‐fixed paraffin‐embedded endobronchial biopsies from subjects of the Newcastle cohort were taken from the third‐ to fourth‐generation airways. They were sectioned and rehydrated for immunofluorescence and immunohistochemistry analysis. Detailed methods are provided in Appendix S1 (Supplementary Information).
Statistical analysis
Data were analysed using Stata software version 15 (StataCorp., College Station, TX, USA) or with GraphPad Prism 8 (San Diego, CA, USA). Results are reported as mean (SD) or median (interquartile range), unless otherwise stated. Continuous measures were analysed using Student’s t‐test or one‐way analysis of variance (ANOVA) as appropriate. Spearman’s correlation coefficients were calculated for univariate association between continuous variables. Linear regression models were used to assess the associations between the variables of interest and ACE2 expression. We initially examined demographic variables alone (bivariate associations and then a multivariable model), and then the clinical characteristics (presence of asthma and COPD) for which the sample was enriched; we also assessed models that included cardiovascular disease, hypertension and diabetes as simple bivariate regression models as well as a model including demographic variables. Due to the relatively low sample size, rather than fitting a full multivariable model and risk over‐fitting, we used a Random Forest to identify the important variables from the collective set. Non‐parametric bootstrapping was used to assess the reliability of results by fitting the model on each of the 50 different samples taken with replacement from the original sample. Regression analyses were performed using R version 4.0.0 (R Core Team (2020), Vienna, Austria).
RESULTS
We assessed gene expression in primary bronchial epithelial cells (pBEC) from participants of the Newcastle cohort (Table 1) and details provided by disease group are summarized separately (Table 2).
Number | 116 |
Age | 62 (13.5) |
Sex | Male 49 (42%), female 67 (58%) |
Smoking status | Never 48 (41%), ex‐smoker 59 (51%), current 9 (8%) |
BMI | 29.6 (6.6) |
Asthma | 39 (33.6%) |
COPD | 38 (33%) |
Atopy | 23/92 (33%) |
Cardiovascular disease | 8 (7%) |
Hypertension | 22 (19%) |
Diabetes | 12 (10%) |
Malignancy | 9 (8%) |
Gastro‐oesophageal reflux | 55 (47%) |
Prescribed regular inhaled corticosteroids | 55 (47%) |
Prescribed regular oral corticosteroids | 5 (4%) |
Prescribed regular ACE or ARB inhibitor | 12 (10%) |
Prescribed regular proton pump inhibitor | 23 (20%) |
- Fourteen subjects had no documented allergy by either skin prick test or radioimmunosorbent assay.
- ACE, angiotensin‐converting enzyme; ARB, angiotensin II receptor blocker; BMI, body mass index; COPD, chronic obstructive pulmonary disease.
Healthy | Asthma | COPD | Analysis | |
---|---|---|---|---|
Number | 39 | 39 | 38 | NA |
Age years (mean and SD) | 64 (55.1, 71.5) | 53.2 (14.8) | 68.2 (8.2) | P = 0.8 |
Sex, male:female | 11:28 | 16:23 | 16:22 | P = 0.7 |
Pack‐years smoked (mean and SD) | 0 | 0 | 38.6 (34.9) | P < 0.001 |
FEV1pp (median and interquartile range) | 90 (83.5, 130.1) | 64.5 (49.8, 86.3) | 46.5 (40, 64) | P < 0.001† |
Regular ICS (%) | 0 | 39 (100%) | 17 (42%) | P = 0.001 |
ICS, median BDP/day, (median and interquartile range) | 0 | 1600 (800, 2000) | 1600 (800, 2000) | P = 0.8‡ |
Atopy | 0 | 19 (49%) | 4 (10.1%) | P < 0.001 |
Airway eosinophils, % total cell count (median and interquartile range) | 1.0 (0.25, 1.6) | 2.5 (1, 6.9) | 0.9 (0.5, 10.1) | P = 0.002§ |
GINA | ||||
GINA step 3 | 16 | |||
GINA step 4 | 3 | |||
GINA step 5 | 20 | |||
TH2 high disease | NA | 28 (72%) | NA | |
GOLD | ||||
GOLD 2 | 15 | |||
GOLD 3 | 19 | |||
GOLD 4 | 3 |
- Pack‐years smoked, 1 packet of 25 cigarettes or equivalent per day for 1 year. ICS dose is expressed as equivalent BDP in mcg per day, where 1 μg of BDP is equivalent to 1 mcg budesonide, 0.5 mcg fluticasone propionate and 0.4 mcg ciclosenide. TH2 high disease was defined as evidence of either atopy or elevated airway eosinophils (2.75% of the total cell count).
- †
FEV1 was significantly lower in participants with COPD compared to both healthy controls (P < 0.001) and participants with asthma (P = 0.01). FEV1 was lower in participants with asthma compared to healthy controls (P < 0.001) (Kruskal–Wallis with Dunn’s comparison).
- ‡
ICS dose was different between those with asthma and COPD.
- §
Participants with asthma had higher airway eosinophils compared to healthy controls (Kruskal–Wallis with Dunn’s comparison).
- BDP, beclomethasone dipropionate; COPD, chronic obstructive pulmonary disease; FEV1, forced expiratory volume in 1 s; FEV1pp, FEV1 percentage of the predicted value; ICS, inhaled corticosteroid; GINA, Global initiative for Asthma; GOLD, Global Initiative for Chronic Obstructive Lung Disease.
Gene expression of ACE2 and related genes
Expression of the following genes was assessed: ACE2, the human receptor that SARS‐CoV‐2 is known to bind, TMPRSS2, CTL and furin, proteases known to be crucial in cleaving the SARS‐CoV‐2 spike protein and facilitating viral entry.3–5, 16 In addition, we examined the expression of genes that may be involved in facilitating SARS‐CoV‐2 entry: TMPRSS11A/D that may also be involved in spike protein cleavage17; a disintegrin and metallopeptidase (ADAM)‐10 and ADAM‐17, which mediate shedding of ACE2 from the cell surface and promote uptake of SARS‐CoV18, 19; and phosphatidylinositol 4‐kinase beta (PI4KB) that creates a lipid microenvironment at the cell surface required for the entry of SARS‐CoV‐1 as well as other RNA viruses.20 We also assessed the gene expression of proposed alternate receptors for the SARS‐CoV‐2 notably CD147, heat shock protein (HSPA) 5 or GRP78 as potential alternative binding sites for the virus.21, 22
We first determined if there were differences in gene expression between subjects who were healthy controls or had a diagnosis of asthma or COPD (Fig. 1). ACE2 expression was significantly lower in those with asthma (mean difference fold change ddCT with healthy controls, 0.23, P = 0.01). Furin was also lower in subjects with asthma (mean difference: 1.1, P = 0.02), but there were no differences between groups in the expression of TMPRSS2 or CTSL. There were no differences in TMPRSS11A/D. ADAM‐10 was higher in COPD (mean difference: 0.15, P = 0.02) and ADAM‐17 was higher in subjects with asthma (mean difference: 0.18, P = 0.02). The groups were not different in CD147 or PI4KB expression. HSPA5 expression was higher in subjects with asthma (mean difference: 0.3, P = 0.02).
Age and male sex are associated with higher gene expression
Older age was associated with increased ACE2 (Fig. 2), although the correlation was not strong (Pearson’s rho = 0.2, P = 0.02). There was no correlation with body mass index (BMI), pack‐years smoked or dose of inhaled corticosteroids (ICS) with any of the genes assessed (data not shown). There were no differences seen in gene expression when analysed separately by obesity (BMI > 30); atopy; history of hypertension, cardiovascular disease or diabetes; use of inhaled or oral corticosteroids; use of ACE inhibitors, ACE2 or proton pump inhibitors (data not shown); and evidence of atopy or airway eosinophil count (data not shown).
When we investigated samples from the Perth cohort, there were no differences in ACE2 gene expression between children and adults, but the negative correlation with age still held. We compared ACE2 gene expression between the Perth participants, the Newcastle healthy controls and participants with asthma. While there were no differences between the healthy participant groups, despite the younger age of the Perth healthy controls compared to those from Newcastle, those with asthma still had significantly less ACE2 expression (Fig. 3).
We performed a bivariate regression analysis for each variable, adjusting for age, sex, smoking, BMI, asthma, COPD, atopy, hypertension, diabetes, cardiac disease, ICS, prednisone and use of ACE2 inhibitors in all subjects. Age and sex had significant associations with ACE2 expression, and collectively the demographics explained 7.5% of the variation in ACE2 (Table 3). Asthma was also associated with lower ACE2 expression; however, this association was not preserved when adjusting for age, sex and smoking.
Combined model coefficients | Estimate | SE | Probability > [T] |
---|---|---|---|
Intercept | 0.788732 | 0.238309 | <0.01 |
Age | 0.004731 | 0.003588 | 0.04 |
Sex | −0.18289 | 0.087999 | 0.04 |
Pack‐years smoked | 0.002361 | 0.001627 | 0.15 |
Asthma | −0.158238 | 0.103480 | 0.13 |
- ACE2 expression was adjusted for age, sex, smoking, pack‐years smoked and weight. Multiple R‐squared: 0.1185, adjusted R‐squared: 0.08619, degrees of freedom, P = 0.008.
- ACE2, angiotensin‐converting enzyme 2.
Expression of ACE2 in the airway epithelium
We stained endobronchial biopsies from five healthy controls, five subjects with asthma and five subjects with COPD, with representative slides displayed (Fig. 4). In comparing slides between healthy controls aged less than 40 years (n = 5) and greater than 65 years (n = 5), there appeared to be a marked difference in ACE2 expression in the epithelium in older donors (Fig. 4A,B). Comparing healthy controls to asthma, there again appeared to be less ACE2 in subjects with asthma (Fig. 4A,C). In healthy controls and subjects with COPD, there was appreciably greater ACE2 observed in the epithelium and subepithelium compared to subjects with asthma (Fig. 4D–F). Most ACE2 expression appeared to be associated with ciliated as well as basal AEC. It was concentrated at the apical end of the AEC, although moderate nuclear and cytoplasmic staining was also observed. The quantification of ACE2 expression (%) in the epithelium with immunohistochemistry staining (IHC) staining indicated less ACE2 expression in subjects with asthma compared to healthy controls and subjects with COPD (Fig. 4G). The percentage (%) ACE2 expression in epithelium was significantly lower in subjects with asthma (median (range): 6.654 (3.955–8.667)) than in subjects with COPD (median (range): 11.084 (9.934–12.136); P = 0.008) and in healthy control (median (range): 9.063 (7.075–12.213); P = 0.05). ACE2 expression (%) in subjects with COPD appeared higher than in controls but failed to reach significance.
DISCUSSION
The COVID‐19 pandemic represents a unique and critical challenge for the world and a greater understanding of the mechanisms of disease overall and in different patient groups is crucial in determining effective therapeutic responses that so far remain elusive. We found ACE2 gene expression is greater in lower airway cells in older individuals and males. We determined that ACE2 gene and protein levels were lower in subjects with asthma. This was also associated with reduced furin expression and increased ADAM‐17 expression in the epithelium from subjects with asthma, demonstrating a potential mechanism through which IL‐13‐induced airway inflammation may downregulate ACE2 expression.
Viral respiratory tract infections are common triggers of acute asthma,23 and initially there was great concern that people with asthma would be highly susceptible to SARS‐CoV‐2. Early reports from China7 and Italy14 did not suggest that people with asthma were over‐represented in those hospitalized from SARS‐CoV‐2. More recently, the largest cohort review of more than 17 000 hospitalized cases of COVID‐19 from the UK has been pre‐published demonstrating that older age, being male and social deprivation were strongly associated with mortality.24 While chronic lung disease, including COPD, was associated with mortality, only those with asthma who had required the recent use of oral corticosteroids had increased mortality.24 We found increased ACE2 was associated with older age and male sex in two independent cohorts in keeping with the epidemiological studies. This suggests that increased lower airway ACE2 expression may be associated with increased risk of more severe disease, and this is consistent with increased ACE2 expression in the small airway epithelium from current smoker and those with COPD.11, 25 Asthma alone is not over‐represented in people with COVID‐19 nor is it associated with greater disease severity. This raises the question as to why people with asthma are not at increased susceptibility to SARS‐CoV‐2.
SARS‐CoV‐2 is a highly infectious virus transmitted by airway droplets.26 Cell entry is the first step necessary and cell tropism is key to infectivity. SARS‐CoV‐2 uses the ACE2 receptor as its point of entry and this receptor has been found to be highly expressed in the upper airway.27, 28 While ACE2 is expressed in the lower respiratory tract, in keeping with our results, there are relatively lower levels of expression compared to the upper airway. SARS‐CoV‐2 is able to first infect ciliated columnar epithelial cells and then appears to spread to surrounding secretory cells.29 ACE2 expression is critical for lower airway infection to occur, and mice deficient in the receptor had greatly reduced SARS‐CoV infection and replication.30 Thus, it is likely that the level of ACE2 expression in lower airway cells may be an important factor in determining more severe infection.
The entry of SARS‐CoV and SARS‐CoV‐2 into cells is facilitated by the interaction between viral S‐protein with the extracellular domain of the transmembrane ACE2 protein. Endocytosis of the virus then results in downregulation of surface ACE2 expression.3, 31, 32 Once infection is established, SARS‐CoV has been shown to upregulate ADAM‐17 that facilitates further viral endocytosis with ACE2, while knockdown of ADAM‐17 by siRNA prior to infection substantially attenuates SARS‐CoV entry.33 In 12 selected patients with COVID‐19 and the 6 of these who developed ARDS, circulating angiotensin II levels were markedly elevated, which correlated with viral load as measured by PCR and the degree of pathology. This suggests a link in SARS‐CoV‐2 infection with tissue ACE2 downregulation facilitating the development of ARDS.34 The spread of infection to the lower airways is reflected in the clinical presentation of the illness, with the first 4–5 days often asymptomatic, followed by dry cough and fever, with viral load peaking at day 10, with those who go onto to develop severe disease, ARDS and severe systemic immune responses 7–10 days after infection.7, 35
The biologic importance of reduced ACE2 gene expression and protein levels that we find in asthma is also supported by the observation of reduced expression of furin, a protease that together with TMPRSS2 has been found to be critical in facilitating SARS‐CoV‐2 endocytosis.16 We also found enhanced expression of ADAM‐17. Phosphorylation of ACE2 enhances its catalytic activity and increases its shedding from the epithelial surface.36 The role of ADAM‐17 is not well defined in asthma but its presence in AEC has been proposed to be anti‐inflammatory and linked to increased airway remodelling and goblet cell hyperplasia.37 Enhanced ADAM‐17 activity is a potential mechanism by which ACE2 receptor expression could be downregulated specifically in asthma.
Three other studies have examined ACE2 mRNA expression in epithelial cells from people with asthma. Saheb Sharif‐Askari et al. examined Biobanked Affymetrix (Santa Clara CA, USA) libraries of bronchial epithelial brushings but did not find reduced ACE2 expression in asthma, although they did find an increase in COPD samples.38 Similar to us, Jackson et al. showed reduced ACE2 expression in nasal samples from children and adults and bronchial epithelial cells from steroid‐naïve mild allergic asthma patients.39 They demonstrated a negative correlation of ACE2 expression with type 2 airway inflammation, exhaled nitric oxide and nasal epithelial IL‐13 expression. We demonstrate ACE2 expression is seen equally in those with non‐allergic and more severe asthma on regular ICS. Interestingly, it has been shown that IL‐13‐induced epithelial proliferation is mediated by ADAM‐17.40 There is now evidence that treatment of AEC cultures with IL‐13 reduces the expression of ACE2.41 This could provide a unifying mechanism by which the effect of IL‐13 on tissue remodelling of the airway epithelium as it occurs in severe allergic disease and asthma may lead to downregulation of ACE2 expression in the lower airways.
In contrast to our findings, Bradding et al. examined bronchial brushings and endobronchial biopsies for the expression of ACE2, TMPRSS2 and furin.42 Results were available from 356 patients (88 healthy volunteers and 268 patients with asthma (mild to severe)). They did not see any difference between participants with asthma or healthy controls but did find weak, although statistically significant positive, correlations with a TH‐17 gene signature and negative correlation with a TH‐2 gene signature. In contrast, we saw reduced ACE2 expression in patients with asthma but did not find a difference in ACE2 expression with atopy or bronchial lavage eosinophil count. The cells we have measured ACE2 expression in had been cultured without stimulation, as opposed to endobronchial brushings. Our endobronchial biopsies were collected in a similar way.
An important possibility remains that ICS may reduce ACE2 expression. Increased expression of ACE2 and TMPRSS2 was also found in over 300 induced sputum samples from a cohort with severe asthma, although the use of inhaled steroids was associated with lower ACE2 expression.43 These results however are from induced sputum and the RNA content will have largely come from immune cells and not the airway epithelium as is the case for us. Belgian researchers have seen a similar trend with reduced ACE2 expression in people with COPD or asthma/COPD overlap on regular ICS.44 While a pre‐print report using in vitro and animal models of ICS in COPD also suggests that they may reduce ACE2 expression.45 All our participants with asthma were on ICS and 42% of those with COPD. We did not see a negative correlation between ICS dose and ACE2 expression and subjects with airways disease on ICS independent of asthma did not have lower ACE2 expression. Leung et al. demonstrated increased ACE2 expression from resected lung tissue and endobronchial brushings from 26 subjects who were current or former smokers as well as subjects with COPD, suggesting increased expression in COPD.25 This relationship is likely to be quite complex and observational data in asthma and COPD will inevitably be confounded by ICS use. It is important to note that neither asthma nor ICS use are associated with increased ACE2 expression and on balance the evidence suggests they may be protective. It will take further mechanistic studies to answer these questions. There are now four actively recruiting clinical trials that are assessing if ICS are protective in COVID‐19.
Our study has some weaknesses. Subjects were recruited to study chronic airways disease and were enriched for subjects with moderate to severe asthma and COPD. Our numbers overall were not large enough to give a clear answer with multivariate analysis and a large proportion with both asthma and COPD were on ICS.
We present the largest pathology study that examined ACE2 gene and protein expression from cultured lower airway bronchial epithelial cells enriched with people with asthma and COPD. We have confirmed the presence of ACE2 in samples from children and adults in two independent cohorts and determined with bivariate analysis that age and male sex are associated with increased lower airway ACE2 expression. Our data provide further evidence for the differences observed in RNA expression confirming reduced ACE2 protein levels in subjects with asthma by immunohistochemistry. We are now the second group to find that people with asthma have reduced ACE2 expression in lower airway cells. We propose that this may provide a degree of protection from the complications of SARS‐CoV‐2. The next step is to determine from current epidemiological studies of COVID transmission what is the risk of infection in people with asthma and whether they are relatively protected. We find reduced ACE2 expression associated with reduced furin expression in subjects with asthma and observe increased expression of ADAM‐17. This provides a potential mechanism by which chronic airway inflammation associated with asthma and the action of IL‐13 on the airway epithelium could downregulate ACE2 expression via ADAM‐17 prior to infection. Further studies should now be performed to define this mechanism and determine whether this can be manipulated as a therapeutic target to reduce SARS‐CoV‐2 infection of the airway epithelium.
Acknowledgements
S.S.S. is supported by the Clifford Craig Foundation Launceston General Hospital and the Rebecca L. Cooper Medical Research Foundation. P.M.H. is funded by the Clifford Craig Foundation Launceston General Hospital, a Rebecca L. Cooper Fellowship and grants from the National Health and Medical Research Foundation Council (NHMRC) of Australia (1079187 and 1175134) and SPHERE.
Author contributions
Conceptualization: P.A.B.W., P.S.P., G.K., S.S.S., S.M.S., P.M.H., A.K. Data curation: P.A.B.W., P.S.P., G.K., A.A., E.N.S., L.W.G., M.S.E., A.R., P.V., A.C.‐Y.H., K.L., A.K. Formal analysis: P.A.B.W., L.C., E.N.S., M.S.E., C.O., A.K. Investigation: K.N., A.A., L.C., L.W.G., S.S.S., W.L., N.B., A.R., P.V., A.C.‐Y.H., K.L., T.I. Methodology: P.S.P., G.K., K.N., E.N.S., L.W.G., S.S.S., W.L., M.S.E., A.R., A.C.‐Y.H., K.L., T.I. Project administration: P.A.B.W., K.N., S.M.S., P.M.H., A.K. Resources: P.A.B.W., E.N.S., S.S.S., N.B., S.M.S., A.K. Supervision: P.A.B.W. Software: L.W.G. Validation: K.N., A.A., E.N.S., L.W.G., M.S.E., A.K. Visualization: K.N., L.C. Writing—original draft: P.A.B.W., P.S.P., G.K., K.N., L.C., E.N.S., L.W.G., S.S.S., M.S.E., C.O., N.B., A.R., K.L., T.I., S.M.S., A.K. Writing—review and editing: P.A.B.W., P.S.P., G.K., K.N., A.A., L.W.G., S.S.S., M.S.E., C.O., A.R., T.I., S.M.S., P.M.H., A.K.
REFERENCES