Research article | Open | Open Peer Review | Published:
Methacholine bronchial provocation measured by spirometry versus wheeze detection in preschool children
BMC Pediatricsvolume 5, Article number: 19 (2005)
Determination of PC20-FEV1 during Methacholine bronchial provocation test (MCT) is considered to be impossible in preschool children, as it requires repetitive spirometry sets. The aim of this study was to assess the feasibility of determining PC20-FEV1 in preschool age children and compares the results to the wheeze detection (PCW) method.
55 preschool children (ages 2.8–6.4 years) with recurrent respiratory symptoms were recruited. Baseline spirometry and MCT were performed according to ATS/ERS guidelines and the following parameters were determined at baseline and after each inhalation: spirometry-indices, lung auscultation at tidal breathing, oxygen saturation, respiratory and heart rate. Comparison between PCW and PC20-FEV1 and clinical parameters at these end-points was done by paired Student's t-tests.
Results and discussion
Thirty-six of 55 children (65.4%) successfully performed spirometry-sets up to the point of PCW. PC20-FEV1 occurred at a mean concentration of 1.70+/-2.01 while PCW occurred at a mean concentration of 4.37+/-3.40 mg/ml (p < 0.05). At PCW, all spirometry-parameters were markedly reduced: FVC by 41.3+/-16.4% (mean +/-SD); FEV1 by 44.7+/-14.5%; PEFR by 40.5+/-14.5 and FEF25–75 by 54.7+/-14.4% (P < 0.01 for all parameters). This reduction was accompanied by de-saturation, hyperpnoea, tachycardia and a response to bronchodilators.
Determination of PC20-FEV1 by spirometry is feasible in many preschool children. PC20-FEV1 often appears at lower provocation dose than PCW. The lower dose may shorten the test and encourage participation. Significant decrease in spirometry indices at PCW suggests that PC20-FEV1 determination may be safer.
Measurements of bronchial hyper-reactivity (BHR) have provided insight into the physiological basis of asthma, and provide a tool for asthma diagnosis, assessment of asthma severity and response to treatment [1, 2]. The bronchial provocation tests require an objective outcome measurement that reflects airway function. Forced expiratory volume in 1 second (FEV1) has been standardized to measure changes in airway caliber that occur with bronchial provocation . In the Methacholine challenge test (MCT), the provocative concentration reducing FEV1 by 20% from baseline (PC20-FEV1) is considered the end point of the test. Traditionally, spirometry in young children has been difficult to achieve. Therefore, techniques that do not require cooperation (i.e., detection of wheeze during normal breathing, a fall of 5% in O2-saturation (SaO2), or an increase of 50% in respiratory rate and/or heart rate) have been used as alternative end points in bronchial provocation tests in the preschool age [4–7]. Recently it has been shown that young children can be taught to perform reliable forced expiratory maneuvers [8–11]. Yet, it is unclear whether these young children have the drive to perform and tolerate repetitive reproducible spirometry-sets that are measured during the interval between inhalations. Concentration of methacholine (MCH) causing wheeze, a fall of 5% in O2- Saturation, an increase of 50% in respiratory rate and/or heart rate (PCW) and PC20-FEV1 were compared in school children and a good correlation was found between the two methods [7, 12–14].
This study assesses the ability of young asthmatic preschool children to cooperate with repetitive spirometry-sets during MCT, and thereby allow determination of PC20-FEV1 in comparison with PCW.
Consecutive preschool children referred to the Pediatric Pulmonary Clinic, Meyer Children's Hospital, Rambam Medical Center, Haifa, over a 6-month period were recruited. Of 62 families offered participation in the study, seven refused. None of the children had experienced spirometry previously. Inclusion criteria were: 2.5–6.5 year-old children who were asthmatic according to GINA guidelines  with recurrent episodes of wheeze, cough and/or shortness of breath with clinical response to bronchodilator; normal chest auscultation and FEV1 >75% of predicted for healthy preschool children  after saline inhalation. Exclusion criteria were: presence of other chronic respiratory conditions; emergency room visit in the past three months; respiratory infection in the past month; oral or inhaled steroids or other anti-inflammatory medication taken in the last week; bronchodilator taken within 24 hours prior to the test.
The Rambam Medical Center Ethics Board approved the study. Parental consent was obtained for each child.
Tests were performed in a designated room at the Pediatric Pulmonary Unit, Meyer Children's Hospital, Haifa, Israel. A parent and the investigating team (a pediatric pulmonary physician, respiratory physiologist and technician) were present throughout the test. MCT was performed according to published guidelines, , with doubling doses of fresh Methacholine solutions (0.06 to 8.00 mg/ml) dissolved in saline. Solutions were driven by compressed air of 5 l/min flow (giving a mean output of 0.4 ml/min), and nebulized using a Hudson nebulizer (Hudson RCI, Temecula, CA, USA). Inhalations were performed using a facemask while the child was sitting up straight and breathing normally. Nebulized Methacholine was inhaled for 2 minutes, with 5-minute intervals between doses, until the maximal concentration or the end point was reached. To ensure safety in light of the risk of airway closer, the MCH increment was only half the usual amount when transient wheeze or cough was noted, keeping in mind that the accumulative dose is affected by this manipulation. Oxygen saturation and heart rate were monitored continuously by pulse oximetry (Biox 3700e; Ohmeda). A single observer (LB) performed auscultation for 20 seconds over the trachea and two zones of both lungs (upper front and lower back) according to Springer et al. 
The following indices were considered "end of test": appearance of audible wheeze, a fall of ≥5% in O2-saturation, or an increase of ≥50% in respiratory rate and/or heart rate . At the "end of test", spirometric measurements were performed, followed by administration of nebulized Albuterol (2.5 mg).
Forced expiratory flow volume (FEFV) curves were measured with a ZAN100 commercial spirometer (ZAN Messgeraete GmbH, Oberthulba, Germany). Calibration was performed before the testing sessions. The curves were monitored on the computer screen to ensure best effort. Results were corrected to BTPS conditions. The software included an interactive animated computer game (SpiroGame®) set by targets of the FEFV maneuver, combining forced inhalation preceding forced expiration, peak expiratory flow rate (PEFR) and forced vital capacity (FVC) with emphasis on prolonged expiration.  The targets were the extrapolated values derived from comparative data from older children, corrected for height.  An experienced pulmonary technician instructed each child how to operate the game. Teaching time was limited to 15 minutes. On-line rejection of curves was based on visual inspection for "non-cooperation" errors and included: poor effort; incomplete expiration; cough; glottis closure. Curves had to show a rapid rise to peak flow, and gradual, smooth decline of flow down to residual volume. Baseline maneuvers were repeated to visually obtain best possible efforts on at least 3 technically acceptable FEFV curves. After obtaining baseline spirometry, MCT was performed. A duplicate spirometry set was performed immediately after auscultation. PC20-FEV1 was determined off line by the provocative concentration that reduced FEV1 by 20% from baseline. PC values were log-transformed before statistical analyses. Spirometry indices included FVC, FEV1, PEFR, forced expiratory flow at 50% FVC (FEF50), FEV1/FVC ratio.
Analysis and statistics
Three baseline spirometry curves were analyzed for acceptability criteria according to ATS/ERS guidelines [17, 18] and in comparison with similar data for preschool children [11, 19]. These included: a) "Start of test" criteria: time to peak expiratory flow and backward extrapolated volume (Vbe) b) "End of test criteria": described by "total expiratory time" and the ratio of "no change in expiratory volume" to "total expiratory time" c) reproducibility (coefficient of variation) of the three baseline curves, calculated as SD/mean*100.
After inhalations, the curves were inspected visually online, and were analyzed offline in relation to baseline using paired t-test. Differences were considered significant when p < 0.05. The level of agreement between the dose at end of test and the dose of PC20 were compared by Bland and Altman analysis (20).
A total of 55 children (28F/27M, age range 2.8–6.4 years) were recruited. Eleven children failed spirometry and underwent MCT by auscultation only. Failure to perform spirometry was due to lack of comprehension (4 children) or failure to repeat spirometry after baseline measurements (7 children). Failure was not age dependent. Eight children refused to cooperate with either test. Thirty-six of 55 (65.5%) children performed the MCT with spirometry tests and with auscultation. Of these 36 children, eleven were 2.5–3.9 years old, 15 were 4–5 years old, and 10 were >5 years old. Three children failed to produce FEV1 on the baseline measurements but were able to produce it after saline administration. In these children, post saline FEV1 measurements were considered as baseline. FEV1 at that point was >75% predicted. The anthropometric data and baseline lung function of the 36 patients are presented in Table 1 and clinical characteristics in Table 2.
The 36 children participating in both tests had a previous response to bronchodilators as judged by clinical observation. The average duration of respiratory symptoms was 18 ± 14 weeks. Five children were not receiving any medication for a period of weeks. Nine children were receiving bronchodilators as needed, and 22 were using both inhaled steroids and bronchodilators as needed.
Quality of baseline maneuvers: Start of test: Peak expiratory flow rates were reached within a mean of 98 ± 7 ms (range 89–115 ms) and mean Vbe was 3.4 ± 1.5% of FVC (range 1.2–5.7). Intra-subject reproducibility for the baseline triple maneuvers was: for FVC, 4.1 ± 2.3% (range 1.8–6.3); for FEV1, 3.8 ± 2.3% (range 0.4–7.3); for PEFR, 4.4 ± 2.8% (range 0.3–8.6) and for FEF25–75, 7.9 ± 3.5% (range 2.7–13.2). End of test: Mean expiratory time was 1.48 ± 0.47 seconds and the ratio of "no change in expiratory-volume" to "total expiration time" was 0.20 ± 0.06.
Children's response to MCT (n = 36) is summarized in Figure 1 and Table 3. Average test time to reach PC20-FEV1 was 29 ± 11 minutes, while for PCW it was 41 ± 10 minutes (not including bronchodilator administration) (p < 0.001). The end point of the challenge was determined by the pediatric pulmonologist as positive in 35/36 children. One child did not display any of the determined criteria for PCW up to 8 mg/ml and was considered to have no BHR. The mean (± SD) concentration at PCW for the 36 children was 4.26± 3.31 mg/ml. Wheezing at the end point was observed in 26/36 children and in 9/36 the test was ended before the appearance of wheeze due to either oxygen desaturation or tachypnea accompanied by audible long expiration. Mean increase in heart rate at PCW was 25.5 ± 11% (range 10–42%); respiratory rate increased by 30.0 ± 21.1% (range 0–42%) and SaO2 decreased by 6.3 ± 2.7% (range 2.3–10.3%).
PC20-FEV1 occurred at a mean concentration value of 1.96- ± 1.83 mg/ml. The one child who did not respond to MCH of up to 8 mg/ml by PCW (negative BHR) did not show a fall of 20% from baseline FEV1 value either. The other 35 children exhibited a fall of 20% in FEV1 from baseline values in response to MCH ≤8 mg/ml (Figure 1 and Table 3). A representative set of FEFV curves from a single patient that includes the predicted curve, baseline, PC20-FEV1 and end of test curves is shown in Figure 2.
At PC20-FEV1 there was a mean increase in heart rate of 13.5 ± 11.0%, respiratory rate increased by 15.4 ± 15.8% and SaO2 decreased by 2.4 ± 2.1% from baseline level. These changes were significantly lower than those found at PCW (p < 0.01 for three parameters). The appearance of PC20-FEV1 occurred 2 concentrations earlier than PCW in 5 children, 1.5-concentrations earlier in 3 children, one concentration earlier in 17 children, 0.5 concentrations earlier in 3 children and at the same concentration as PCW in 7 children (Figure 1). The effects of MCH on the spirometry parameters are presented in Table 4. At PC20-FEV1, parameters were moderately decreased, while at end point, test parameters were markedly reduced. The severity of FEV1 reduction at PCW was variable, ranging from 30.8 to 68.2% of baseline. The level of agreement between the dose at end of test (PCW) and the dose at PC20 is presented in Figure 3. Dotted lines represent 95% coefficient of variation values.
Bronchodilators improved FEV1 by 43 ± 29% from PCW values and all respiratory symptoms disappeared shortly after bronchodilator administration.
In this study we assessed the feasibility of determining PC20-FEV1 during Methacholine bronchial provocation testing in asthmatic preschool children. We found MCT was feasible in 65% of this group of wheezy preschool children. Children as young as 3 years old complied and cooperated with what seems to be a most fatiguing procedure. Baseline measurements met most of the ATS criteria for older children and adults [17, 18] and quality control studies on spirometry in preschool children [11, 19]. We found that PC20-FEV1 correlates with PCW. However, PC20-FEV1 frequently precedes PCW. All spirometry parameters at PC20-FEV1 were significantly higher than those measured at PCW.
In this study, we used interactive spirometry games  with multiple spirometry targets, since single targeted games (usually peak expiratory flow targeted) have not fulfilled expectations [21, 22]. Our teaching method is supported by the findings that 65% of the children fully cooperated not only with baseline measurements but also with spirometry sets. Of note, 26 of the 36 children were younger than 5 years. Conforming quality control was necessary to proceed with the test. The quality control of baseline spirometry in our study met most ATS/ERS criteria concerning reproducibility and start of test criteria [17, 18] and matched those reported for preschool children [11, 19], encouraging us to continue with the MCT test. Vbe = 5%FVC found in our study is narrower than reported , as we have rejected in advance curves with Vbe >5%FVC at the expense of success rate. It should be stressed that our work did not compare verbal coaching  or other spirometry games [11, 23] as the preferable methodology for keeping the child going and performing repetitive spirometry sets.
The mean PC20-FEV1 of 1.96 ± 1.83 mg/ml found in our group reflects a mild degree of BHR, as we recruited children with mild asthmatic symptoms. Our findings for PC20-FEV1 are comparable to those of Hayden et al , who found a mean PC20-FEV1 at FEV0.5 of 2.49 ± 2.55 mg/ml in infants. Adinoff et al  reported a mean provocative dose of 3.0 mg/ml Methacholine in their preschool children and infants. Tepper  et al. reported that infants with asthma-like respiratory symptoms might respond to MCH concentrations as low as 1.25 mg/ml. In that respect we found that PC20-FEV0.5 occurred at a mean concentration value of 1.29- ± 1.47 mg/ml, meaning that the responsiveness of the airways in the preschool age may be similar to that of infants, despite differences in the measurement techniques. It is important to note that PC20-FEV0.5 occurred at a significant mean lower concentration than PC20-FEV1 (1.96- ± 1.83 mg/ml; p < 0.01), however, standardization is needed to accept the PC20-FEV0.5 value for the determination of hyper-reactive airways.
We found that PCW occurred in our group at a mean concentration of 4.26 ± 3.31 mg/ml. PCW values in our study were much higher than the PCW (0.4 mg/ml) reported by Springer et al . The difference may be attributed to inclusion of more severe asthmatics in their study group.
Spirometry at PCW
We found that PC20-FEV1 occurred at a lower concentration than PCW in most subjects. This finding is in agreement with several other studies comparing PCW detection to PC20-FEV1 in school age children. [4–7, 14]. However, in none of these studies were spirometry measurements carried out to the point of wheeze. We expected to find a good correlation between the two tests (PCW = 1.2195*PC20-FEV1 + 0.0288; R2 = 0.9733; p < 0.005), yet the Blant and Altman analysis revealed that in children with higher mean provocation PCW dose (≥6 mg), the level of agreement between the methods was low, reflecting higher sensitivity of the PC20 method, especially in mild airway reactivity (Figure 3).
We found that at PCW, FEFV curves visually seemed to be smaller and all parameters were reduced simultaneously (Figure 2), with a highly significant reduction in flows and volume parameters. The reduction in curve was gradual in most children, accompanied by an increase in respiratory symptoms (Table 2), and responded to bronchodilators, and hence was not considered to reflect fatigue. To further strengthen this point a representative curve of one child illustrating, a poor effort performed at teaching process vs. end of test curve is shown in figure 4. The poor-effort curve did not fulfill start of test criteria and is round while the "end of test curve" has an obstructed shape.
The reduced FVC and flows are most likely due to a severe degree of airway narrowing involving small to medium airways that may be accompanied by air trapping, partial closure of airways and elevation in FRC. Reduced FVC may also be due to increased glottic narrowing due to MCH irritation [26, 27], but the flow volume curve was not suggestive of upper airway obstruction (trimmed PEFR). Alternatively, the upper airways response to methacholine may contribute to the increase in total respiratory resistance . This pattern occurred in some cases before appearance of wheeze or other clinical end-points. Indeed, in 9/36 subjects, the test was terminated due to oxygen desaturation or tachypnea rather than wheeze. Similar to our results, Sprikkelman et al  reported that wheeze was detected in only 33% of 15 school-age asthmatic children at PC20-FEV1, and Springer et al  terminated the test without the presence of wheeze in 19.2% of young children. In this respect we would argue that FEV1 does make a contribution beyond simply asking the subject if they wheeze. Novitzki et al (4) found in 5–8 year-old children that FEV1 is decreased by 33.3 ± 7.4% at PCW. Spence et al  reported a mean fall of 51 ± 14% from baseline FEV1 when wheeze appeared in their asthmatic older subjects. Our results strengthen these prior findings, and suggest that spirometric PC20-FEV1 may be achieved with inhalation of lower MCH concentrations than those used to achieve wheeze.
Measuring PCW during tidal volume breathing has the advantage that no active cooperation on the child's part is needed. Therefore the success rate of PCW is higher than spirometry (44/55 children). However, using PC20-FEV1 (or PC20-FEV0.5) can preclude inhalation of higher concentrations of MCH used to achieve wheeze, leading to alarmingly diminished flows found at PCW and a significant shortening of test time relative to PCW.
We conclude that PC20-FEV1 is feasible in preschool asthmatic children when using respiratory games teaching techniques and that the children tolerate repetitive duplicate sets of spirometry maneuvers. PC20-FEV1 in preschool children appears to be as sensitive as in adults and school children. Yet, many questions remain open as to the usefulness of this test in a random sample of young children and/or how discriminating this test is as a diagnostic tool. It would also be necessary to assess the sensitivity of this test to various severities of disease. Further studies are needed for standardization and definition of methodological criteria.
Cockcroft DW: Bronchoprovocation methods: direct challenges. Clin Rev Allergy Immunol. 2003, 24: 19-26. 10.1385/CRIAI:24:1:19.
Spiropoulos K, Stevens J, Eigen H, Spiropoulos A: Specificity and sensitivity of methacholine challenge test in children with normal and hyper-reactive airways. Acta Paediatr Scand. 1986, 75: 737-743.
Guidelines for Methacholine and exercise challenge testing - 1999. The Official Statement of the American Thoracic Society. Am J Respir Crit Care Med. 2000, 161: 309-329.
Noviski N, Cohen L, Springer C, Bar-Yishay E, Avital A, Godfrey S: Bronchial provocation determined by breath sounds compared with lung function. Arch Dis Child. 1991, 66: 952-955.
Wilson NM, Bridge P, Silverman M: The measurement of methacholine responsiveness in 5-year-old children: three methods compared. Eur Respir J. 1995, 8: 364-370. 10.1183/09031936.95.08030364.
Wilts M, Hop WC, van der Heyden GH, Kerrebijn KF, de Jongste JC: Measurement of bronchial responsiveness in young children: comparison of transcutaneous oxygen tension and functional residual capacity during induced bronchoconstriction and dilatation. Pediatr Pulmonol. 1992, 12: 181-185.
Springer C, Godfrey S, Picard E, Uwyyed K, Rotschild M, Hanania S, Noviski N, Avital A: Efficacy and Safety of Methacholine Bronchial Challenge performed by auscultation in young asthmatic children. Am J Respir Crit Care Med. 2000, 162: 857-860.
Vilozni D, Barker M, Jellouschek H, Heimann G, Blau H: An interactive computer-animated system (SpiroGame) facilitates spirometry in pre-school children. Am J Respir Crit Care Med. 2001, 164: 2200-2205.
Eigen H, Bieler H, Grant D, Christoph K, Terrill D, Heilman DK, Ambrosius WT, Tepper RS: Spirometric pulmonary function in healthy preschool children. Am J Respir Crit Care Med. 2001, 163: 619-623.
Zapletal A, Chalupova J: Forced expiratory parameters in healthy preschool children (3–6 years of age). Pediatr Pulmono. 2003, 35: 200-207. 10.1002/ppul.10265.
Aurora P, Stocks J, Oliver C, Saunders C, Castle R, Chaziparasidis G, Bush A: Quality control for spirometry in preschool children with and without lung disease. Am J Respir Crit Care Med. 2004, 169: 1152-1159. 10.1164/rccm.200310-1453OC.
Sprikkelman AB, Grol MH, Lourens MS, Gerritsen J, Heymans HS, van Aalderen WM: Use of tracheal auscultation for the assessment of bronchial responsiveness in asthmatic children. Thorax. 1996, 51: 317-319.
Hayden MJ, Devadason SG, Sly PD, Wildhaber JH, LeSouef PN: Methacholine responsiveness using the raised volume forced expiration technique in infants. Am J Respir Crit Care Med. 1997, 155: 1670-1675.
Martinez FC, Yarza GPE, Ruiz AA, Knorr EJI, Blecua CM, Aramburu MJ: Agreement between tracheal auscultation and pulmonary function in methacholine bronchial inhalation challenge in asthmatic children. An Esp Pediatr. 2002, 56: 304-309.
Global Strategy for Asthma Management and Prevention. 1995, National Institutes of Health, National Heart, Lung and Blood Institute, Bethesda
Zapletal A, Samanek M, Prague TP: Lung function in children and adolescents: Methods, reference values. Basel, Karger. 1987
Standardization of spirometry-1994 update. Statement of the American Thoracic Society. Am J Respir Crit Care Med. 1995, 153: 1107-1236.
Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, Peslin R, Yernault JC: Lung volumes and forced ventilatory flows. Report Working Party Standardization of Lung Function Tests, European Community for Steel and Coal. Official Statement of the European Respiratory Society. Eur Respir J. 1993, 5-40.
Crenesse D, Berlioz M, Bourrier T, Albertini M: Spirometry in children aged 3 to 5 years: reliability of forced expiratory maneuvers. Pediatr Pulmonol. 2001, 32: 56-61. 10.1002/ppul.1089.
Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 1: 307-310.
Nystad W, Samuelsen SO, Nafstad P, Edvardsen E, Stensrud T, Jaakkola JJ: Feasibility of measuring lung function in preschool children. Thorax. 2002, 57: 1021-1027. 10.1136/thorax.57.12.1021.
Gracchi V, Boel M, Van der Laag J, Van der Ent CK: Spirometry in young children: should computer-animation programs be used during testing?. Eur Respir J. 2003, 2: 872-875. 10.1183/09031936.03.00059902.
Kozlowska W, Aurora P, Stocks J: The use of computer-animation programs during spirometry in preschool children. Eur Respir J. 2004, 2: 494-495.
Adinoff AD, Schlosberg RT, Strunk RC: Methacholine inhalation challenge in young children: results of testing and follow-up. Ann Allergy. 1988, 61: 282-286.
Tepper RS: Airway reactivity in infants: a positive response to methacholine and metaproterenol. J Appl Physiol. 1987, 62: 1155-1159.
England SJ, Ho V, Zamel N: Laryngeal constriction in normal humans during experimentally induced bronchoconstriction. J Appl Physiol. 1985, 58: 352-356.
Marchal F, Loos N, Monin P, Peslin R: Methacholine-induced volume dependence of respiratory resistance in preschool children. Eur Respir J. 1999, 14 (5): 1167-74. 10.1183/09031936.99.14511679.
Sprikkelman AB, Schouten JP, Lourens MS, Heymans HS, van Aalderen WM: Agreement between spirometry and tracheal auscultation in assessing bronchial responsiveness in asthmatic children. Respir Med. 1999, 93: 102-107. 10.1016/S0954-6111(99)90298-6.
Spence DP, Graham DR, Jamieson G, Cheetham BM, Calverley PM, Earis JE: The relationship between wheezing and lung mechanics during methacholine induced broncho-constriction in asthmatic subjects. Am J Respir Crit Care Med. 1996, 154: 290-294.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2431/5/19/prepub
The Study was funded by the Israel Lung Association, Tel-Aviv, Israel
The author(s) declare that they have no competing interest.
The SpiroGame program is privately patented in USA, granted to Dr. Vilozni.
Dr. Vilozni does not foresee any financial gain or loss, now or in the future from publishing this manuscript. The patent is not commercialized.
Dr. Lea Bentur and Dr. Daphna Vilozni had primary responsibility for protocol development, outcome assessment, data analysis and writing of the manuscript.
Dr. Raphael Beck, Dr. Nael Elias, Dr. Asher Barak, Dr. Ori Efrati and Prof. Yaacov Yahav contributed to this study by patients screening, patient enrollment, analysis of the data and quality control of the data.