key: cord-0688476-97ywyhaj
authors: Mazzotta, Valentina; Cozzi-Lepri, Alessandro; Colavita, Francesca; Lanini, Simone; Rosati, Silvia; Lalle, Eleonora; Mastrorosa, Ilaria; Cimaglia, Claudia; Vergori, Alessandra; Bevilacqua, Nazario; Lapa, Daniele; Mariano, Andrea; Bettini, Aurora; Agrati, Chiara; Piselli, Pierluca; Girardi, Enrico; Castilletti, Concetta; Garbuglia, Anna Rosa; Vaia, Francesco; Nicastri, Emanuele; Antinori, Andrea
title: Emulation of a Target Trial From Observational Data to Compare Effectiveness of Casirivimab/Imdevimab and Bamlanivimab/Etesevimab for Early Treatment of Non-Hospitalized Patients With COVID-19
date: 2022-04-20
journal: Front Immunol
DOI: 10.3389/fimmu.2022.868020
sha: 27c230ab1aa726b006c23cf79b2b7d28928a2445
doc_id: 688476
cord_uid: 97ywyhaj

OBJECTIVES: Comparative analysis between different monoclonal antibodies (mAbs) against SARS-CoV-2 are lacking. We present an emulation trial from observational data to compare effectiveness of Bamlanivimab/Etesevimab (BAM/ETE) and Casirivimab/Imdevimab (CAS/IMD) in outpatients with early mild-to-moderate COVID-19 in a real-world scenario of variants of concern (VoCs) from Alpha to Delta. METHODS: Allocation to treatment was subject to mAbs availability, and the measured factors were not used to determine which combination to use. Patients were followed through day 30. Viral load was measured by cycle threshold (CT) on D1 (baseline) and D7. Primary outcome was time to COVID-19-related hospitalization or death from any cause over days 0-30. Weighted pooled logistic regression and marginal structural Cox model by inverse probability weights were used to compare BAM/ETE vs. CAS/IMD. ANCOVA was used to compare mean D7 CT values by intervention. Models were adjusted for calendar month, MASS score and VoCs. We evaluated effect measure modification by VoCs, vaccination, D1 CT levels and enrolment period. RESULTS: COVID19-related hospitalization or death from any cause occurred in 15 of 237 patients in the BAM/ETE group (6.3%) and in 4 of 196 patients in the CAS/IMD group (2.0%) (relative risk reduction [1 minus the relative risk] 72%; p=0.024). Subset analysis carried no evidence that the effect of the intervention was different across stratification factors. There was no evidence in viral load reduction from baseline through day 7 across the two groups (+0.17, 95% -1.41;+1.74, p=0.83). Among patients who experienced primary outcome, none showed a negative RT-PCR test in nasopharyngeal swab (p=0.009) and 82.4% showed still high viral load (p<0.001) on D7. CONCLUSIONS: In a pre-Omicron epidemiologic scenario, CAS/IMD reduced risk of clinical progression of COVID-19 compared to BAM/ETE. This effect was not associated with a concomitant difference in virological response.

Objectives: Comparative analysis between different monoclonal antibodies (mAbs) against SARS-CoV-2 are lacking. We present an emulation trial from observational data to compare effectiveness of Bamlanivimab/Etesevimab (BAM/ETE) and Casirivimab/ Imdevimab (CAS/IMD) in outpatients with early mild-to-moderate COVID-19 in a realworld scenario of variants of concern (VoCs) from Alpha to Delta.

Methods: Allocation to treatment was subject to mAbs availability, and the measured factors were not used to determine which combination to use. Patients were followed through day 30. Viral load was measured by cycle threshold (CT) on D1 (baseline) and D7.Primary outcome was time to COVID-19-related hospitalization or death from any cause over days 0-30. Weighted pooled logistic regression and marginal structural Cox model by inverse probability weights were used to compare BAM/ETE vs. CAS/IMD. ANCOVA was used to compare mean D7 CT values by intervention. Models were adjusted for calendar month, MASS score and VoCs. We evaluated effect measure modification by VoCs, vaccination, D1 CT levels and enrolment period.

The widespread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, causing coronavirus disease-2019 (COVID- 19) , continues to be a challenge for global public health.

Monoclonal antibodies (mAbs) against SARS-CoV-2 have emerged as the strategy of choice for the treatment of early mildto-moderate COVID-19 in outpatients at increased risk of clinical progression (1) and, based on data from randomized clinical trials (RCT), Bamlanivimab/Etesevimab (2, 3), Casirivimab/Imdevimab (4, 5) and Sotrovimab (6) received emergency use authorizations by the Italian Medicines Agency (AIFA) (7) .

In addition to RCT, data from observational cohorts (8) (9) (10) (11) have also confirmed the effectiveness of mAbs, but comparative analysis between different options available are lacking (12, 13) . Considering the rapid epidemiological evolution, with the emergence of new Variants of Concern (VoCs) (14) (15) (16) that have been shown to escape (17) the action of mAbs in vitro (18, 19) and in vivo (20) , real-life data about clinical impact and mAb comparison are useful to better clarify the scenario of currently existing drugs.

The aim of this analysis was to compare the clinical effectiveness of two mAb combinations, Bamlanivimab/ Etesevimab (BAM/ETE) and Casirivimab/Imdevimab (CAS/ IMD), in a real-life setting during a period in which the prevalent lineages in Italy were B.1.1.7 (Alpha) and B.1.617.2 (Delta), and B.1.529 (Omicron) was not yet circulating.

On March 2021, mAb administration program started at the National Institute for Infectious Diseases Lazzaro Spallanzani IRCCS in Rome, ruled by AIFA. Eligibility criteria included outpatients with a confirmed diagnosis of SARS-CoV-2 infection by an antigenic or molecular nasopharyngeal swab (NPS), mild or moderate symptoms of COVID-19 for 10 days or less and at least one of following risk factors for progression to severe disease: body mass index (BMI) ≥35, chronical peritoneal dialysis or hemodialysis, uncontrolled or complicated diabetes mellitus, primary or secondary immunodeficiency. Subjects who were 55 years and older were also eligible if they had any cardiocerebrovascular diseases, Chronic Obstructive Pulmonary Disease (COPD) or other chronic respiratory diseases. Patients requiring hospitalization for COVID-19 or supplemental oxygen therapy were excluded. In June 2021 (21) , AIFA expanded use of mAbs including all patients with one of the following: 65 years and older, BMI >30, any chronic renal impairment (including subjects undergoing peritoneal dialysis or hemodialysis), uncontrolled or complicated diabetes mellitus, any immunocompromising condition [including primary and secondary immunodeficiencies (1)], cardio-cerebrovascular diseases (including hypertension with concomitant organ damage), COPD or other chronic respiratory diseases, chronic liver disease, hemoglobinopathies, neurodevelopmental and neurodegenerative disease.

Comorbidity burden was assessed using Monoclonal Antibodies Screening Score [MASS (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) ] that assigned points, as follows: age ≥65 (2 points), BMI ≥35 (1 point), diabetes mellitus (2 points), chronic kidney disease (3 points), cardiovascular disease in a patient ≥55 years (2 points), chronic respiratory disease in a patient ≥55 years (2 points), hypertension in a patient ≥55 years (1 point), and immunocompromised status (3 points).

All consecutive adult patient (age>18) who provided a written informed consent were included in study population. The study was approved by AIFA and National Ethics Committee.

Bamlanivimab/Etesevimab (700 mg/1400 mg) or Casirivimab/ Imdevimab (1200 mg/1200mg) were administered by one-hour intravenous infusion and patients were observed for one hour after infusion. Allocation to treatments was pseudo-random, as a criterion of daily alternation (subject to drugs availability) was adopted, and as not many of the measured factors were used to determine which combination to infuse.

Patients who received Bamlanivimab as monotherapy or Sotrovimab were excluded.

Outpatients visits were scheduled at baseline (D1) and at 7 (D7) and 30 (D30) days after infusion. Medical evaluation, vital signs recording, laboratory tests and reports on adverse effects were performed at each visit. If patients missed person visits, they were called by telephone to assess clinical conditions. As a real-life study, due to the high diagnostic demands related to the COVID-19 pandemic, different methods were used to investigate serology and virological parameters according to the laboratory workflow and tests availability.

SARS-COV-2 serology was performed by ELISA detecting anti-SARS-CoV-2 IgG, IgM, and IgA (ENZY-WELL SARS-CoV-2; DIESSE, Diagnostica Senese, Siena, Italy; positive index values ≥1.1), or by two chemiluminescence microparticle assays (CMIA) detecting anti-Nucleoprotein and anti-Spike/RBD IgG (ARCHITECT SARS-CoV-2 IgG, and ARCHITECT SARS-CoV-2 IgG II Quantitative; Abbott Laboratories, Wiesbaden, Germany, respectively). According to the to manufacturer's instructions, for the two CMIA, Index >1.4 and Binding Antibody Units (BAU)/mL ≥7.1 are considered positive for anti-N and anti-Spike/RBD IgG, respectively.

Semi-quantitative estimation of viral load in NPS was assessed by RT-PCR using DiaSorin Simplexa ® COVID-19 Direct platform (DiaSorin, Saluggia, Italy), based on cycle threshold (CT) values of S and ORF1ab genes amplification. Other RT-PCR methods used to verify the presence of SARS-CoV-2 were the Abbott m2000 RealTime System (Abbott Laboratories, Wiesbaden, Germany) and the Cobas ® SARS-CoV-2 Test on the fully-automated cobas ® 6800 Systems (Roche Diagnostics, Rotkreuz, Switzerland).

Identification of VoCs was conducted by Sanger sequencing of the Spike coding gene on the D1 samples. During the period of Delta variant wave, the RT-PCR Simplexa ® SARS-CoV-2 Variants Direct kit (DiaSorin, Saluggia, Italy). was included in the study as rapid method for the qualitative detection and differentiation of the N501Y, E484K, E484Q and L452R mutations.

The primary outcome was defined as time to hospitalization due to development of severe COVID-19 or death from any cause over days 0-30. Secondary outcomes were a) time to hospitalization or death from any cause by day 30, b) time to hospitalization or death from any cause over days 3-30, c) the impact of intervention on CT values change from baseline to D7. The proportion of participants reporting adverse events were also shown as a safety endpoint.

Main characteristics of the participants, assessed at D1, were compared by treatment strategy using Mann-Whitney U test for continuous variables, expressed as median (IQR) or c 2 test or Fisher's exact test as appropriate for Categorical variables expressed as numbers and percentages.

The effectiveness of the two strategies on the three outcomes was estimated and compared using a weighted pooled logistic regression model which approximates the parameters of a marginal structural Cox model by mean of inverse probability weights. Participants' follow-up accrued from the date of infusion until the date of hospitalisation, death or date of discharge. Administrative censoring was also applied at 31/12/ 2021 the date at which the database was frozen. Weights have been calculated using the predicted values from the pooled logistic models for the probabilities of starting BAM/ETE vs. CAS/IMD and those of censoring, respectively. Treatment was fitted as a time-fixed variable and there was no need to account for immortal time bias. Potential informative censoring was controlled for using inverse probability of censoring weights. Unweighted and weighted hazard ratios (HRs) with 95% confidence intervals (CI) were shown. Assumptions regarding the underlying causal link between measured factors are pictured in Figure 1 . Unweighted and weighting Kaplan-Meier estimates of the primary outcome stratified by treatment strategy were fitted. Interactions between the intervention and study population strata (type of VoC, vaccination status, baseline SARS-COV-2 serology, level of D1 CT values and period of enrolment) were formally tested by including a multiplicative term in the marginal Cox regression model and adjusted hazard ratios (HRs) with 95% confidence intervals (CI) were shown in a forest plot.

Positive serology was defined as detection of IgG and/or IgA with ELISA or detection of IgG anti-S and/or anti-N by CMIA.

We compared mean CT values at D7 by treatment in analysis of covariance (ANCOVA-model adjusted for D1 CT value, month of enrolment and type of VoC) and we described, using box-plots, median (IQR) of D1 and D7 CT values and its variation.

We also evaluated the associations between D7 CT response (using the cut-offs of 40 for negativity and 25 for reduced viral load) with both the intervention and the primary outcome using a chi-square and Fisher exact test as appropriate.

A descriptive analysis of self-reported side effects was also performed.

A two-sided test of less than 0.05 was considered statistically significant. All statistical analyses were performed using the SAS software, version 9.4 (Carey USA). Of the 61 vaccinated patients tested with ELISA, 41 (67%) showed a positive SARS-COV-2 serology; of the 138 vaccinated patients tested by CMIA, 130 (94%) were positive to anti-Spike antibodies, and only 7 were positive to anti-N antibodies. Among vaccinated patients, mRNA-1273 was used more frequently in patients receiving BAM/ETE, and ChAdOx1 less frequently (p=0.05). Overall, the participants receiving the two treatment strategies appeared to be balanced with respect to key predictors of outcomes as expected under our pseudorandom allocation design. Main characteristics according to the two treatment groups are reported in Table 1 .

COVID19-related hospitalization or death from any cause occurred in 19 participants: 15 patients in the BAM/ETE group (6.3%) and 4 patients in CAS/IMD group (2%) ( Table 2) . Two deaths were observed, both in patients treated with BAM/ETE experiencing COVI-19 clinical progression.

The majority of the events occurred before D7. In the weighted Kaplan-Meier analysis there was greater evidence for a difference in risk by intervention group (>10% for BAM/ETE) than in the unweighted (Figures 3A, B) .

The relative hazard of COVID-19-related hospitalization or death for CAS/IMD compared to BAM/ETE was 0.28 (95%CI 0.09-0.85; p=0.024, Table 3A ).

Hospitalization or death from any cause by day 30 occurred in 29 participants (9% in BAM/ETE vs. 4% in CAS/IMD). Overall, 18 of these secondary events (62%) occurred within day three after infusion ( Table 2) .

Tables 3B,C, showed the weighted relative hazard ratio for secondary endpoints.

The analysis by subsets carried no evidence that the effect of the intervention was different across a number of stratification variables for the primary endpoint (p-values for interaction >0.18, Figure 4 ). The effect of the intervention appeared attenuated in participants with the Delta VoC but with large uncertainty around these subgroup estimates. 

Thirty-one participants reported adverse events; a breakdown of the most common adverse events reported is described in Table 4 . Overall, 19 events (10%) in the CAS/IMD group and 6 events (3%) in the BAM/ETE group were considered related to mAb infusion (p=0.001). One patient only (0.2%) reported severe dyspnea requiring hospitalization: this event was not considered by the investigators to be related to treatment.

To our knowledge, this analysis is the first real-life comparison of two routinely used anti-SARS-CoV-2 mAbs by a trial emulation methodology using observational data collected from outpatients with early mild-to-moderate COVID-19 and at high-risk for disease progression. In an evolving scenario of SARS-CoV-2 variants from Alpha to Delta, we found that patients receiving Patients with hospitalization or death for any cause over days 0-30. c Patients with hospitalization or death for any cause over days 3-30. **adjusted for calendar month, MASS score and VoCs. FIGURE 4 | Forest plot of the effect of the intervention by subsets -Primary endpoint. *p-value corresponds to the test for interaction between intervention and each subgroup unadjusted for multiplicity. **aHR1-(vaccination, serology and CT strata) adjusted for calendar month and VoCs; **aHR2 -(VoCs stratum) adjusted for MASS score because of a positivity issue leading to abnormally high weights in the stratum of participants infected with the Delta VoC; **aHR3 -(enrolment period stratum) adjusted for MASS score and VoCs. PDFU, person days of follow-up.

Casirivimab/Imdevimab had a 72% lower risk of COVID-19 related-hospitalization or death from any cause than patients receiving Bamlanivimab/Etesevimab. The greater benefit of Casirivimab/Imdevimab was evident also excluding patients who experienced failure by day2 (a time window in which events cannot be ascribed to lack of treatment effect), with a confirmed 78% risk reduction.

Overall, we found no evidence that the magnitude of the difference of the effect between mAb interventions on the risk of clinical outcomes was different in specific subsets of the study population. Nevertheless, the effect of Casirivimab/Imdevimab appeared to be larger in subgroups with the highest known risk of disease progression: patients with higher baseline viral load (23) , unvaccinated (24, 25) and those enrolled before June 2021, when the target population included people with higher risk of severe outcome (26) .

On the other hand, there was larger uncertainty around the hazard ratio comparing interventions in the subgroup of patients infected with the Delta VoC with a 95% CI not excluding superiority of Bamlanivimab/Etesevimab vs. Casirivimab/ Imdevimab. This finding is consistent with in vitro observed retained activity of Bamlanivimab/Etesevimab (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) against Delta variant (B.1617.2, non-AY.1/AY.2), and supports the recommended use of both Bamlanivimab/Etesevimab and Casirivimab/Imdevimab in settings of elevated Delta VoC prevalence.

To our knowledge, this is the first analysis from real-life data evaluating virological response to mAb treatment. It has been suggested that mAbs may act as antiviral neutralizing agents through multiple mechanisms, such as targeting free virus and virally infected cells (28) . Significant decrease in viral load was described in RCT (2-4) for the two mAbs combinations analyzed in this study, but relation between virological and clinical outcomes remains uncertain (29). Interestingly, our analysis showed no difference in terms of viral load reduction from D1 to D7 between Casirivimab/Imdevimab and Bamlanivimab/ Etesevimab, but also displayed a very strong association between clinical and virological outcomes, suggesting that patients developing severe disease also failed in viral clearing. However, viral load was measured after the occurrence of all clinical events, so it is difficult to determine how much clinical outcome was mediated by the virological response or whether lack of virological clearance was actually a consequence of the clinical picture. Finally, our data confirmed safety and tolerability of these two mAb combinations in a real-life unselected population.

Our analysis has some limitations. First, due to the observational nature of the study conducted in a single COVID health care center and to the lack of a randomized design, confounding bias cannot be ruled out. Further, eligibility criteria changed over time concurrently with the advent of Delta wave and with a wider use of Bamlanivimab/Etesevimab, due to available supplies. However, results were similar after controlling for MASS score in the regression models. Moreover, the lack of an early measure of CT (e.g. at D3) prevented us from investigating viral load as a potential mediator. Finally, the study was conducted before the emergence of Omicron B.1.1.529 VoC, which is going to subvert previous assessment about mAbs treatment as several in vitro studies suggest that both Bamlanivimab/Etesevimab and Casirivimab/Imdevimab did not retain a remarkable activity against Omicron (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) . Despite this, even today a proportion of illnesses and consequent hospitalizations are still due to the VoCs different from Omicron, and so knowledge of comparative data between available mAbs is still crucial for optimizing treatment in pandemic times.

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

The studies involving human participants were reviewed and approved by National Ethics Committee. The patients/ participants provided their written informed consent to participate in this study.

AA, VM conceptualized and designed the study. VM and IM wrote the protocol. VM, FC and AC wrote the first draft of the manuscript and referred to appropriate literature. AC was also the main responsible person for formal data analysis. AA, AC, VM, CA, FC, CCa and SL conceived, supervised the study and contributed to data interpretation. CCi and PP were responsible for data curation. AV, IM, AG, SC, EG, EN, FV revised the manuscript content, reviewed and edited the manuscript. LE, CF, EL, FC, CCa and AG performed all virological test. IM, AV, SL and SR enrolled participants. All authors agreed with and approved the final version of the manuscript. 

The study was performed in the framework of the SARS-CoV-2 surveillance and response program implemented by the Lazio Region Health Authority. This study was supported by funds to the Istituto Nazionale per le Malattie Infettive Lazzaro Spallanzani IRCCS, Rome (Italy), from Ministero della Salute (Ricerca Finalizzata COVID-2020-12371675 and Ricerca Corrente linea 1 on emerging and reemerging infections").

COVID-19) Treatment Guidelines

Effect of Bamlanivimab as Monotherapy or in Combination With Etesevimab on Viral Load in Patients With Mild to Moderate COVID-19: A Randomized Clinical Trial

Bamlanivimab Plus Etesevimab in Mild or Moderate Covid-19

REGN-COV2, a Neutralizing Antibody Cocktail, in Outpatients With Covid-19

REGEN-COV Antibody Combination and Outcomes in Outpatients With Covid-19

Early Treatment for Covid-19 With SARS-Cov-2 Neutralizing Antibody Sotrovimab

Authorisation for the Temporary Distribution of Monoclonal Antibody-Based Medicines for the Treatment of COVID-19

Monoclonal Antibody Treatment Program. Real-World Clinical Outcomes of Bamlanivimab and Casirivimab-Imdevimab Among High-Risk Patients With Mild to Moderate Coronavirus Disease

Real-World Effect of Monoclonal Antibody Treatment in COVID-19 Patients in a Diverse Population in the United States

Real-World Effectiveness and Tolerability of Monoclonal Antibody Therapy for Ambulatory Patients With Early COVID-19

Casirivimab-Imdevimab Treatment Is Associated With Reduced Rates of Hospitalization Among High-Risk Patients With Mild to Moderate Coronavirus Disease-19

A Learning Health System Randomized Trial of Monoclonal Antibodies for

Effectiveness of Casirivimab and Imdevimab, and Sotrovimab During Delta Variant Surge: A Prospective Cohort Study and Comparative Effectiveness Randomized

SARS-CoV-2 Variants B.1.351 and P.1 Escape From Neutralizing Antibodies

Reduced Sensitivity of SARS-CoV-2 Variant Delta to Antibody Neutralization

Broadly Neutralizing Antibodies Overcome SARS-Cov-2 Omicron Antigenic Shift

Complete Map of SARS-CoV-2 RBD Mutations That Escape the Monoclonal Antibody LY-Cov555 and Its Cocktail With LY-Cov016

Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7

Escape From Neutralizing Antibodies by SARS-Cov-2 Spike Protein Variants

In Vivo Monoclonal Antibody Efficacy Against SARS-CoV-2 Variant Strains

A Framework for Outpatient Infusion of Anti-Spike Monoclonal Antibodies to High-Risk Patients With Mild to Moderate Coronavirus Disease-19: The Mayo Clinic Model

A Randomized, Placebo-Controlled Clinical Trial of Bamlanivimab and Etesevimab Together in High-Risk Ambulatory Patients With COVID-19 and Validation of the Prognostic Value of Persistently High Viral Load

Estimating the Early Impact of the US COVID-19 Vaccination Programme on COVID-19 Cases, Emergency Department Visits, Hospital Admissions, and Deaths Among Adults Aged 65 Years and Older: An Ecological Analysis of National Surveillance Data

Association Between Mrna Vaccination and COVID-19 Hospitalization and Disease Severity

Monoclonal Antibody Treatment of Breakthrough COVID-19 in Fully Vaccinated Individuals With High-Risk Comorbidities

SARS-CoV-2 Variant B.1.617 Is Resistant to Bamlanivimab and Evades Antibodies Induced by Infection and Vaccination

Bamlanivimab Reduces Nasopharyngeal SARS-CoV-2 RNA Levels But Not Symptom Duration in Non-Hospitalized Adults With

Considerable Escape of SARS-CoV-2 Variant Omicron to Antibody Neutralization

Omicron Escapes the Majority of Existing SARS-CoV-2 Neutralizing Antibodies

SARS-CoV-2 Omicron-B.1.1.529 Leads to Widespread Escape From Neutralizing Antibody Responses

The authors gratefully acknowledge nurse staff, all the patients and all members of the Outpatient clinic for MAbs administration and of laboratory: Chiara Agrati, Andrea Antinori, Tommaso Ascoli-Bartoli, Francesco Baldini, Bartolini

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.Copyright © 2022 Mazzotta, Cozzi-Lepri, Colavita, Lanini, Rosati, Lalle, Mastrorosa, Cimaglia, Vergori, Bevilacqua, Lapa, Mariano, Bettini, Agrati, Piselli, Girardi, Castilletti, Garbuglia, Vaia, Nicastri and Antinori. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.