Abstract
Keywords
Introduction
Pulmonary atresia and ventricular septal defect (PA-VSD) is a complex and heterogeneous congenital cardiopulmonary malformation characterized by a deficiency of luminal continuity from the right ventricle (RV) to the pulmonary arteries (PAs), combined with a large malaligned VSD resulting from anterior deviation of the conal septum. 1
There are two concurrent treatment strategies, primary one-stage biventricular repair (p-BVR) and staged surgical repair (SR), with a greater survival rate in p-BVR patients than in those treated with SR but also a significantly greater incidence of right ventricular outflow tract (RVOT) obstruction.2–5 In contrast, patients treated with SR have a decreased probability of requiring reoperation or reintervention during follow-up compared with patients treated with p-BVR but also have a higher mortality after shunt palliation and during the interim period before the Rastelli-type procedure.3,4 Furthermore, many studies have shown that a substantial proportion of patients have early distal conduit or PA branch stenosis after RVOT reconstruction, especially infants and young children.6–8
In Vietnam, in the conditions of a developing country, we have no homograft material, as well as limited fetal diagnosis, and late detection is rather common. Before 2016, almost all patients in our center received Blalock-Taussig-Thomas (BTT) shunt surgery followed by Rastelli surgery at 1.5 to 2 years of age due to limited access to the conduit and concerns that the PA was too small to undergo total repair. The interruption of Gore product distribution in Vietnam poses significant challenges for patients with congenital heart defects with ductal dependencies, particularly those who need an aortopulmonary shunt or Fontan operation. Since 2016, the availability of the Contegra conduit has shifted our treatment strategy for this complex group of disorders, resulting in early complete Rastelli SR for young infants or neonates with lower body weight. However, we are faced with the question of whether early p-BVR or SR will increase the burden of postoperative intervention as in previous studies.
The primary objective of this study was to describe our experience with operative repair of ductal-dependent PA-VSD patients via both p-BVR and SR approaches and to determine the impact of branch PA arterioplasty on clinical outcomes, especially related to RVOT reintervention.
Materials and Methods
Patients and Definitions
This study was approved by the Vietnam National Children's Hospital Research Ethics Council (Ref No: 1076/BVNTU-HDDD in June 2022). A retrospective review was performed for all data collected, and informed consent was waived according to our National Legislation and Institutional Requirements.
Between April 2014 and March 2023, a total of 205 patients diagnosed with PA-VSD underwent SR at Vietnam National Children's Hospital. Patients who were diagnosed with PA-VSD-MAPCA (61 patients) were excluded, and 144 patients who were PA-VSD and ductal-dependent were analyzed—70 patients who underwent p-BVR and 74 patients who underwent SR.
In accordance with our national legislation and institutional requirements, informed consent was waived because of the observational nature of the study. Preoperative, intraoperative, postoperative, and follow-up data were acquired via retrospective review of medical records.
Patient follow-up was completed via regular checkups at 1 month, 3 months, 6 months, and every year after operation at the outpatient clinic department. Contact was maintained with the patient's family via telephone and postal or through direct contact with the referring hospital to obtain external medical records when the patient lived far from Hanoi.
The criteria for p-BVR were a Nakata index greater than 130 mm2/m2 and echocardiogram demonstrating balanced ventricles. The BTT shunt was indicated as the first palliation if the patients had a Nakata Index smaller than 130 mm2/m2 or patients admitted to our hospital with severe cyanosis combined with multiple organ failure due to late presentation. A central shunt was indicated as the primary operation if severe hypoplasia of both PA branches occurred. Ductal stenting and right ventricle-to-pulmonary artery (RV-PA) shunts were rarely performed in selected patients.
Early mortality was defined as death occurring in the hospital before discharge or within 30 days after the initial operation. Interim mortality was defined as death occurring after discharge from the hospital after shunt palliation and before the Rastelli-type procedure in the SR group. Patients who survived after the shunt operation but completed the Rastelli operation at other institutions (four patients) were excluded from the statistical analysis of the Rastelli procedure. All other deaths were considered late deaths. RVOT reintervention was defined as reoperation or cardiac catheterization related to RVOTO, distal conduit or PA branch stenosis, or severe conduit regurgitation.
Pulmonary Artery Diameter Measurement
A pediatric cardiologist reviewed all echocardiograms, multislice computed tomography (CT), and cardiac catheterization data. The right PA (RPA) and left PA (LPA) diameters measured immediately at the hilum during systole were used to determine the Z scores, and the Nakata index was calculated as the sum of the left and right pulmonary cross-sectional areas indexed to the body surface area. The diameter of the branch PAs, each branch PA Z score, and the Nakata index were compared before the first operation.
Surgical Techniques
The Rastelli-type procedure was performed via median sternotomy with cardiopulmonary bypass and moderate hypothermia. The PA branches were dissected up to the hilum to perform pulmonary arterioplasty with a beating heart to reduce the ischemic time of the aortic cross-clamp. If the patient had a small native PA, irrespective of the diameter, the proximal PA was ligated and divided, with the PA stump preserved as a PA bifurcation. Pulmonary arterioplasty was performed on individual PA branches with separate triangular untreated autologous pericardial patches. If the native PA did not exist, a single rectangular pericardial patch was used to perform pulmonary arterioplasty from the hilum to the hilum, and the central part of the patch was perforated to create the distal RV-PA connection. The distal RV-PA conduit was then sutured to the PA bifurcation before cardioplegia arrest (Figure 1). In neonates or small infants with a body weight less than 5 kg, we always consider reducing the RV-PA conduit size by performing resection of one cusp of the Contegra conduit 12 mm in diameter to create a non-oversized bicuspidized conduit 8 to 9 mm (Supplemental Video 1).

Surgical techniques in pulmonary atresia and ventricular septal defect (PA-VSD) patients who underwent the Rastelli operation via individual autologous pericardium patches to increase the size of the PA branches and create a PA bifurcation that was as close to normal as possible. MAPCA, major aortopulmonary collateral arteries.
Our institutional protocol for the BTT shunt was as follows: (1) The size of the shunt was calculated according to the body weight (<3 kg: 3 mm; 3-4 kg: 3.5 mm and >4 kg: 4 mm polytetrafluoroethylene graft). (2) The shunt direction was contralateral to the direction of the aortic arch, which was the left side shunt in a patient with a right aortic arch. A central shunt was placed from the ascending aorta to the main PA in patients with severe PA hypoplasia. Patent ductus arteriosus (PDA) stenting was performed in patients whose left PDA and right aortic arch were affected. The median time from the first palliation to the Rastelli-type procedure was 1.02 years (IQR, 0.78-1.34 years).
As per our protocol, all patients received anticoagulant/antiplatelet medications following the BTT shunt and the Rastelli operation. Heparin was administered at 10 IU/kg/h for the first three to five days following the procedure, and aspirin was administered at a dose of 3 to 5 mg/kg/day for the rest of the patient's life.
Data Analysis
Categorical variables are presented as absolute frequencies and percentages. Continuous variables are presented as medians (interquartile ranges) for skewed data. Categorical variables were compared via Pearson χ2 test. Continuous variables were compared via Wilcoxon rank-sum test. Fisher exact test was used for cases in which the expected frequency was <5. Predictors for death and RVOT reintervention were determined via a multivariable Cox proportional hazards regression model with a baseline model of multivariable models that included all variables that had a P value <.15. A multivariate variable with a P value <.05 with a hazard ratio and 95% confidence interval (CI) was defined as a multivariate risk factor. The Fine-Gray competing risk regression model was used to plot the cumulative incidence function curve, which compared the two groups’ reoperation hazards. Aalen Johansen methods were used to estimate the survival functions while accounting for competitive risk. Time-dependent endpoints for survival and freedom from reintervention were analyzed via the Kaplan-Meier method, while the log-rank test was used to compare Kaplan-Meier curves, with a P value <.05 indicating statistical significance.
Results
Patient Demographic Characteristics
Among 144 patients who underwent surgical intervention for PA-VSD, 84 (58.3%) were males. These 144 patients were classified into two groups according to the primary operation of each group: the p-BVR group (n = 70, 48.6%) and the SR group (n = 74, 51.4%). There were 10 premature babies (14.3%, 10/70) in the p-BVR group and 11 (14.9%, 11/74) in the SR group.
The median age and median weight at the Rastelli operation in the p-BVR group were 3.4 months (IQR, 1.67-7.47 months) and 4.65 kg (IQR, 3.70-6.85 kg), respectively. Eleven of 70 patients (15.7%) underwent surgery during the neonatal period, and a total of 31 of 70 patients (44.3%) were less than 3 months of age in the p-BVR group. In the SR group, the median age and weight at palliation were 1.4 months (IQR, 0.85-3.1 months) and 3.5 kg (IQR, 3-4.5 kg), respectively, whereas 24 of 74 patients (32.4%) underwent shunt palliation in the neonatal period, and 51 patients (68.9%) were less than three months of age. The Nakata index in the p-BVR group was 174 mm2/m2 (IQR, 139-226 mm2/m2), whereas it was 109 mm2/m2 (IQR, 88-138 mm2/m2) in the SR group (P < .001). Figure 2 illustrates the patients’ progression diagram for the study's completion, and Table 1 demonstrates the patients characteristics before the first operation.

Flowchart for patients’ progression for the study completion.
Patient Characteristics Between Groups Before the First Operation.
Abbreviations: IQR, interquartile range; LPA, left pulmonary artery; PA, pulmonary artery; p-BVR, primary biventricular repair; RPA, right pulmonary artery; SR, staged repair.
Perioperative Variables and Mortality
In the p-BVR group, RV-PA reconstruction involved a Contegra conduit in 64 of 70 patients (91.4%), a transannular autologous pericardial patch in four of 70 patients (5.7%), a Hancock valve conduit in 1 of 70 patients (1.4%), and a polytetrafluoroethylene conduit in 1 of 70 patients (1.4%). Sixty-four patients (91.4%, 64/70) underwent LPA patch arterioplasty, whereas 58 of 70 (82.9%) underwent RPA patch arterioplasty, all with autologous pericardium. The median size of the RV-PA conduit was 12 mm (IQR, 12-14 mm) in the p-BVR group, whereas the median size of the RV-PA conduit was 16 mm (IQR, 14-16 mm) in the SR group (P < .001).
In the SR group, the main palliative procedure was a BTT shunt (85.1%, 63/74), followed by a central shunt (9.5%, 7/74), ductal stenting (4.1%, 3/74), and an RV-PA shunt (1.3%, 1/74). Fifty patients (67.5%, 50/74) in the SR group completed the Rastelli-type procedure at our institution at a median of 12 months (IQR, 9.5-16.3 months) after shunt palliation, whereas four patients (5.4%, 4/74) underwent definitive repair at other institutions (these four patients had all survived at the last follow-up). In 50 patients who underwent the Rastelli procedure at our institution, the RV-PA connection involved a Contegra conduit in 44 patients (88%, 44/50), a Hancock valve conduit in three of 50 patients (6%), a polytetrafluoroethylene valve conduit in two of 50 patients (4.0%), and a transannular patch using autologous pericardium in one of 50 patients (2%). Forty-six of 50 patients (92%) underwent LPA patch arterioplasty, and 41 of 50 (82%) underwent RPA patch arterioplasty concomitant with the Rastelli procedure using autologous pericardium as its only material. The perioperative characteristics of the patients who underwent the Rastelli operation are described in Table 2.
Patient Characteristics at Rastelli Operation.
Abbreviations: ECMO, extracorporeal membrane oxygenation; IQR, interquartile range; LPA, left pulmonary artery; PA, pulmonary artery; p-BVR, primary biventricular repair; RPA, right pulmonary artery; RVOT, right ventricular outflow tract; RV/LV, right ventricle/left ventricle; SR, staged repair.
Overall, 11 in-hospital deaths (11/144, 7.6%) occurred, with five deaths (5/70, 7.1%) in the p-BVR group and six (6/74, 8.1%) in the SR group. There were four late deaths (4/70, 5.7%) in the p-BVR group, but eight interstage deaths and two deaths (10/74, 13.5%) after the Rastelli procedure in the SR group. The causes of death and patient characteristics are described in Supplemental Table 1. The overall survival rates at one year and six years after initial surgery were 86% (95% CI, 81%-92%) and 83% (95% CI, 76%-89%), respectively. The survival rates for each group at one year and six years after initial surgery were 91% (95% CI, 85%-98%) and 89% (95% CI, 82%-97%), respectively, in the p-BVR group and 81% (95% CI, 73%-91%) and 76% (95% CI, 66%-87%), respectively, in the SR group (P = .048) (Figure 3). Multivariate Cox analysis revealed that the risk factor for mortality after the Rastelli procedure was an aortic cross-clamp time of more than 80 min (HR, 5.72; 95% CI, 1.26-25.9; P = .024).

Kaplan-Meier curves for survival stratified according to the type of repair.
Right Ventricular Outflow Tract Reintervention and Risk Factors for RVOT Reintervention
Twenty patients (13.8%, 20/144) underwent 27 RVOT reinterventions during a follow-up of 2.6 years (IQR, 0.9-5.1 years). In the p-BVR group, 15 patients (21.4%, 15/70) underwent reintervention: cardiac RVOT reoperation in five patients (7.1%, 5/70), cardiac catheterization for branch PA dilatation in six patients (8.6%, 6/70), and both reinterventions in four patients (5.7%, 4/70). Five patients (10%, 5/50) required reintervention in the SR group: RVOT reoperations were performed in two patients (4%, 2/50), and both RVOT percutaneous interventions combined with reoperation were performed in three patients (6%, 3/50). The overall rates of freedom from RVOT reintervention after the Rastelli operation at one year and six years were 94% (95% CI, 90%-99%) and 79% (95% CI, 69%-91%), respectively, at a median of 1.9 years (IQR, 0.9-5.0 years) from the Rastelli operation to the first reintervention. The rate of freedom from reintervention for the overall type of RVOT, including stenosis or severe regurgitation at six years, was 69% (95% CI, 53%-90%) in the p-BVR group compared with 91% (95% CI, 82%-100%) in the SR group (P = .059) (Figure 4).

Kaplan-Meier curves for freedom from right ventricular outflow tract (RVOT) reintervention. (A) Freedom from overall RVOT reintervention; (B) freedom from RVOT reintervention between groups; (C) freedom from the distal conduit and PA branch reintervention; and (D) freedom from reintervention of RVOT stratified by LPAs Z-score. Abbreviations: p-BVR, primary biventricular repair; LPA, left pulmonary artery.
There were 12 patients with distal conduit and/or PA branch stenosis among the 20 patients who needed RVOT reintervention, resulting in rates of freedom from distal conduit and PA branch reintervention at one year and six years of 94% (95% CI, 90%-99%) and 86% (95% CI, 78%-94%), respectively. The rate of freedom from reintervention for distal conduit and PA branch stenosis at six years was 82% (95% CI, 72%-94%) in the p-BVR group and 91% (95% CI, 82%-100%) in the SR group (P = .2). The risk factors for RVOT reintervention after definitive repair were patients who required PGE1 transfusion before the first operation (HR, 4.4; 95% CI, 1.73-11.2; P = .002) and a Z score of less than −1.4 for the LPA (HR, 4.78; 95% CI, 1.82-12.6; P = .002).
A competing risk analysis was performed, with the events being survival without RVOT reintervention, death without RVOT reintervention, and required RVOT reintervention. The cumulative incidence of survival without RVOT reintervention in the p-BVR group at seven years after initial surgery was 68.6% (95% CI, 55.9%-78.3%), that of death without RVOT reintervention was 8.6% (95% CI, 4.1%-18.9%), and that of RVOT reintervention was 20% (95% CI, 11.4%-30.2%); those in the SR group were 85.2% (95% CI, 71.7%-92.6%), 3.7% (95% CI, 6.7%-11.4%), and 9.3% (95% CI, 3.1%-19.7%), respectively (Figure 5). Fine and Gray's model for cumulative incidence across groups revealed that being in the p-BVR group was a high-risk factor for RVOT reintervention but not significant relative to being in the SR group (P = .12, SHR = 2.22, 95% CI, 0.82-6.03).

The cumulative incidence between groups of survival and freedom from right ventricular outflow tract (RVOT) reintervention.
Discussion
The adverse effects of the Rastelli operation include a heavy burden of RVOT reintervention in patients with ductal-dependent PA-VSD, especially in neonates or small children who undergo p-BVR.3–5 Given the circumstances in developing countries, the rise in reoperations and reinterventions shortly after surgery will impose additional pressure on the system's limited resources, particularly for high-volume centers. Our study suggested that concomitant pulmonary arterioplasty, especially with the primary Rastelli procedure, created a PA bifurcation that was as close to normal as possible and that integration with a nonoversize conduit might have the potential to limit distal anastomotic stenosis, allowing the natural growth of PA branches after definitive repair of the PA-VSD and potentially reducing the development of PA branch stenosis. In contrast with previous studies,2–5 our Fine and Gray's model for cumulative incidence across groups demonstrated that the rate of freedom from RVOT reintervention in the p-BVR group is lower but did not differ substantially (P = .12) from that in the SR group, although the weight and age of patients at the time of the Rastelli operation in the p-BVR group were significantly lower than those in the SR group (Figure 5). Moreover, for significantly smaller conduits in the p-BVR group than in the SR group (P < .001), the freedom for reintervention for distal PA and PA branch stenosis did not differ between the groups (P = .2), at least at the six-year follow-up after the Rastelli operation.
Impact of Concomitant Pulmonary Arterioplasty on Right Ventricular Outflow Tract Reintervention
As shown in previous studies, the management of RVOTO after the Rastelli operation in patients with PA-VSD requires multiple reinterventions but combines multiple approaches, including balloon angioplasty catheterization, stent implantation, and reoperation, especially with primary definitive repair.2–5 Kim et al 5 reported that RVOTO occurred much earlier in the p-BVR group, with 42.9% and 28.6% freedom from RVOT reintervention at one year and two years, respectively. Lee et al 4 also demonstrated a sharp decline in freedom from RVOT reintervention during the first two years of life after the Rastelli operation. Data from previous studies show that after RVOT reconstruction, even with a reasonable preoperative Nakata index and equivalent diameter of the RPA and LPA, recurrent RVOTO and hypoplastic PA branches quickly develop.
The mechanism of recurrent RVOTO is complex, but we hypothesized that the main reason is that direct conduit implantation to small PA branches, without a natural PA bifurcation, will create shear stress flow directly on the PA to induce quick stenosis at the distal conduit and PA branches. Another issue in previous studies is that a significant proportion of patients underwent p-BVR with the transannular patch technique. Given the tiny size of the PA branches and the persistently high PA resistance, one of the potential hazards associated with surgery in neonates is the development of severe pulmonary regurgitation. Therefore, establishing a competent pulmonary valve was essential to the success of postoperative resuscitation in our situation, when iNO was frequently inaccessible. In ductal-dependent PA-VSD patients, our surgical approach of enlarging individual PA branches and forming a natural-like PA bifurcation shows promise in reducing recurrent RVOTO, although the use of a small conduit has the unavoidable consequence of requiring RVOT reoperation with conduit replacement. The median time in our study from Rastelli operation to the first RVOTO reintervention was 1.9 years. This demonstrates the potential for delayed RVOT reintervention using our strategy of concomitant PA arterioplasty.
Additionally, previous studies revealed that distal conduit stenosis and younger age were risk factors for RVOT reintervention after Contegra insertion, with freedom from RVOT reintervention ranging from 24.5% to 52.9% at 5 to 6 years.6–9 Our results modestly suggest that the potential midterm benefit of concomitant PA arterioplasty by the creation of a PA bifurcation that is as close to normal as possible and PA branch enlargement can optimize the outcomes of the Rastelli operation by not only reducing distal PA and PA branch reintervention but also delaying conduit explantation as much as possible for children aged up to six years to implant the largest RV-PA conduit (Table 3). Furthermore, the presence of neonates and small infants was not a risk factor for RVOT reintervention, and a superior rate of freedom from distal conduit and PA branch stenosis (82% at six years) in our p-BVR group, in contrast with previous studies, potentially provides new insight into improving outcomes after primary Rastelli surgery for neonates and young infants. Although our results were encouraging for performing p-BVR without an increased risk of reintervention, a smaller LPA appeared to be a multivariate risk factor for RVOT reintervention, suggesting that patients with LPA coarctation who requested close follow-up after the Rastelli operation and LPA diameter might not benefit in terms of growth, as described by previous studies,10,11 irrespective of extensive arterioplasty at the time of the Rastelli operation.
Baseline Characteristics of Outcomes Within Studies of Ductal-Dependent PA-VSD.
Abbreviations: PA-VSD, pulmonary atresia and ventricular septal defect; p-BVR, primary biventricular repair; RVOTO, right ventricular outflow tract obstruction; SR, staged repair.
Risk Factors for Death and Survival Between Groups
As seen in previous studies,2,5 we have seen the superiority of p-BVR over SR in patients with ductal-dependent PA-VSD. Obviously, the interim deaths after shunt palliation contributed substantially to the overall mortality from SRs, especially in low- to middle-income countries with limited resources to follow-up patients with a vulnerable balance of pulmonary circulation. Consequently, p-BVR for ductal-dependent PA-VSDs has become a vital option for midterm outcomes because it can achieve immediate restoration of normal circulation, resulting in comparable overall survival to that reported in previous studies.2–4 Additionally, it does not impose a heavy burden on early RVOT catheterization or reoperation in our healthcare system owing to resource constraints. However, our patient characteristics might reveal the potential for natural selection with the most severe group of genetic/extracardiac anomalies, which is not a multivariate predictor of death but rather of the complexity of the Rastelli procedure with a longer ischemic time. Nonetheless, we do not present any patient exclusion criteria in our center except for a limited number of patients whose families refuse medical treatment for their children with multiple extracardiac anomalies.
Limitations
This study is limited by its retrospective nature, nonrandomized design, and selection bias. The primary limitation of our study is the lack of a control group of PA-VSD patients who underwent the Rastelli operation without concomitant PA angioplasty.
Although echocardiography and multislice CT results were carefully examined by one cardiologist, accurate measurements could be obtained via cardiac catheterization. The statistical analyses were limited to a small number of patients and outcomes.
Conclusions
Concomitant pulmonary arterioplasty may be associated with increased freedom from RVOT reintervention and delayed RVOT reintervention, provide good midterm outcomes, and should be considered in PA-VSD patients who undergo the Rastelli operation. Primary BVR combined with pulmonary arterioplasty for PA-VSD patients is correlated with better survival than SR (with high interim-stage deaths) and it does not increase the burden of RVOT catheterization or reoperation.
Supplemental Material
sj-docx-1-pch-10.1177_21501351251363158 - Supplemental material for Outcome of Infants With Ductal-Dependent Pulmonary Atresia and Ventricular Septal Defect in a Developing Nation: A Single-Center Experience
Supplemental material, sj-docx-1-pch-10.1177_21501351251363158 for Outcome of Infants With Ductal-Dependent Pulmonary Atresia and Ventricular Septal Defect in a Developing Nation: A Single-Center Experience by Vinh Quang Tran, Anh Vuong Doan, Van-Anh Thi Nguyen, Thai Quoc Ngo, Mai Tuan Nguyen, Doan Vu Tu Quyen, Duyen Dinh Mai, Nhat Huu Nguyen, Doan Quoc Hung, James St. Louis and Truong Ly Thinh Nguyen in World Journal for Pediatric and Congenital Heart Surgery
Footnotes
Abbreviations
Acknowledgments
The authors would like to thank Dr Chang-Ha Lee for his instructions and valuable advice on how to perform p-BRV for ductal-dependent PA-VSD in our hospital.
Author Contributions
Truong Nguyen Ly Thinh contributed conceptualization, project administration, investigation, visualization, writing—original draft, and writing—review and editing. Vinh Quang Tran contributed data curation, validation, and writing—review and editing. Anh Vuong Doan contributed data curation, formal analysis, validation, and writing—review and editing. Van—Anh Thi Nguyen contributed data curation, formal analysis, validation, and writing—review and editing. Thai Quoc Ngo contributed investigation, visualization, and writing—review and editing. Mai Tuan Nguyen contributed investigation, visualization, and writing—review and editing. Doan Vu Tu Quyen contributed investigation, visualization, and writing—review and editing. Duyen Dinh Mai contributed investigation, visualization, and writing—review and editing. Nhat Huu Nguyen contributed investigation, visualization, writing—review and editing. Doan Quoc Hung contributed investigation, visualization, and writing—review and editing. James St. Louis contributed conceptualization, project administration, investigation, visualization, writing—original draft, and writing—review and editing.
Data Availability Statement
The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
