<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Physiol.</journal-id>
<journal-title>Frontiers in Physiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Physiol.</abbrev-journal-title>
<issn pub-type="epub">1664-042X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1112115</article-id>
<article-id pub-id-type="doi">10.3389/fphys.2023.1112115</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Unstable SpO<sub>2</sub> in preterm infants: The key role of reduced ventilation to perfusion ratio</article-title>
<alt-title alt-title-type="left-running-head">Stoecklin et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fphys.2023.1112115">10.3389/fphys.2023.1112115</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Stoecklin</surname>
<given-names>Benjamin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1663595/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Choi</surname>
<given-names>Y. Jane</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2158905/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Dassios</surname>
<given-names>Theodore</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1294928/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jones</surname>
<given-names>J. Gareth</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2121307/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Lockwood</surname>
<given-names>Geoffrey G.</given-names>
</name>
<xref ref-type="aff" rid="aff6">
<sup>6</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Pillow</surname>
<given-names>J. Jane</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/459734/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Neonatology</institution>, <institution>University Children&#x2019;s Hospital Basel (UKBB)</institution>, <addr-line>Basel</addr-line>, <country>Switzerland</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Human Sciences</institution>, <institution>The University of Western Australia</institution>, <addr-line>Crawley</addr-line>, <addr-line>WA</addr-line>, <country>Australia</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Wal-Yan Respiratory Research Centre</institution>, <institution>Telethon Kids Institute</institution>, <addr-line>Nedlands</addr-line>, <addr-line>WA</addr-line>, <country>Australia</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Neonatal Intensive Care Unit</institution>, <institution>King&#x2019;s College Hospital NHS Foundation Trust Denmark Hill</institution>, <addr-line>London</addr-line>, <country>United Kingdom</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Cambridge University Clinical School</institution>, <addr-line>Cambridge</addr-line>, <country>United Kingdom</country>
</aff>
<aff id="aff6">
<sup>6</sup>
<institution>Anaesthetic Department</institution>, <institution>Hammersmith Hospital</institution>, <addr-line>London</addr-line>, <country>United Kingdom</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1218440/overview">Rajasvaran Logeswaran</ext-link>, Asia Pacific University of Technology &#x26; Innovation, Malaysia</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1251821/overview">Kevin Dysart</ext-link>, Alfred I. duPont Hospital for Children, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2133996/overview">Juliann Di Fiore</ext-link>, Case Western Reserve University, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Benjamin Stoecklin, <email>benjamin.stoecklin@ukbb.ch</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Respiratory Physiology and Pathophysiology, a section of the journal Frontiers in Physiology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>07</day>
<month>02</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1112115</elocation-id>
<history>
<date date-type="received">
<day>30</day>
<month>11</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>25</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Stoecklin, Choi, Dassios, Jones, Lockwood and Pillow.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Stoecklin, Choi, Dassios, Jones, Lockwood and Pillow</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>
<bold>Introduction:</bold> Instability of peripheral oxyhemoglobin saturation (SpO<sub>2</sub>) in preterm infants is correlated with late disability and is poorly understood. We hypothesised that a reduced ventilation to perfusion ratio (V<sub>A</sub>/Q) is the key predisposing factor for SpO<sub>2</sub> instability.</p>
<p>
<bold>Methods:</bold> We first used a mathematical model to compare the effects of reduced V<sub>A</sub>/Q or shunt on SaO<sub>2</sub> stability (SaO<sub>2</sub> and SpO<sub>2</sub> are used for model and clinical studies respectively). Stability was inferred from the slope of the SaO<sub>2</sub> vs. inspired oxygen pressure (<italic>P</italic>
<sub>I</sub>O<sub>2</sub>) curve as it intersects the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line (breathing air). Then, in a tertiary neonatal intensive care unit, paired hourly readings of SpO<sub>2</sub> and <italic>P</italic>
<sub>I</sub>O<sub>2</sub> were recorded over a 24&#xa0;h period in week old extremely preterm infants. We noted SpO<sub>2</sub> variability and used an algorithm to derive V<sub>A</sub>/Q and shunt from the paired SpO<sub>2</sub> and <italic>P</italic>
<sub>I</sub>O<sub>2</sub> measurements.</p>
<p>
<bold>Results:</bold> Our model predicted that when V<sub>A</sub>/Q &#x3c; 0.4, a 1% change in <italic>P</italic>
<sub>I</sub>O<sub>2</sub> results in &#x3e;8% fluctuation in SaO<sub>2</sub> at 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub>. In contrast, when a 20% intrapulmonary shunt was included in the model, a 1% change in <italic>P</italic>
<sub>I</sub>O<sub>2</sub> results in &#x3c;1% fluctuation in the SaO<sub>2</sub>. Moreover, further reducing the V<sub>A</sub>/Q from 0.4 to 0.3&#xa0;at 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> resulted in a 24% fall in SaO<sub>2</sub>. All 31 preterm infants [mean gestation (&#xb1;standard deviation) 26.2 (&#xb1;1)&#xa0;week] had V<sub>A</sub>/Q &#x3c; 0.74 (normal &#x3e;0.85) but only two infants had increased shunt at 1.1 (&#xb1;0.5)&#xa0;weeks&#x2019; postnatal age. Median (IQR) SpO<sub>2</sub> fluctuation was 8 (7)%. The greatest SpO<sub>2</sub> fluctuations were seen in infants with V<sub>A</sub>/Q &#x3c; 0.52 (<italic>n</italic> &#x3d; 10): SpO<sub>2</sub> fluctuations ranged from 11%&#x2013;17% at a constant <italic>P</italic>
<sub>I</sub>O<sub>2</sub> when V<sub>A</sub>/Q &#x3c; 0.52. Two infants had reduced V<sub>A</sub>/Q and increased shunt (21% and 27%) which resolved into low V<sub>A</sub>/Q after 3&#x2013;6&#xa0;h.</p>
<p>
<bold>Discussion:</bold> Routine monitoring of <italic>P</italic>
<sub>I</sub>O<sub>2</sub> and SpO<sub>2</sub> can be used to derive a hitherto elusive measure of V<sub>A</sub>/Q. Predisposition to SpO<sub>2</sub> instability results from reduced V<sub>A</sub>/Q rather than increased intrapulmonary shunt in preterm infants with cardiorespiratory disease. SpO<sub>2</sub> instability can be prevented by a small increase in <italic>P</italic>
<sub>I</sub>O<sub>2</sub>.</p>
</abstract>
<kwd-group>
<kwd>infant</kwd>
<kwd>premature</kwd>
<kwd>neonatal intensive care unit</kwd>
<kwd>pulmonary gas exchange</kwd>
<kwd>oxygen inhalation therapy</kwd>
</kwd-group>
<contract-num rid="cn001">GNT1057514 GNT1047689 RF1077691 GNT1196188</contract-num>
<contract-sponsor id="cn001">National Health and Medical Research Council<named-content content-type="fundref-id">10.13039/501100000925</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Schweizerischer Nationalfonds zur F&#xf6;rderung der Wissenschaftlichen Forschung<named-content content-type="fundref-id">10.13039/501100001711</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Peripheral oxyhemoglobin saturation (SpO<sub>2</sub>) instability may result in more than 100 hypoxemic events per day within the first 8&#xa0;weeks of life in preterm infants (<xref ref-type="bibr" rid="B9">Di Fiore et al., 2010</xref>). Moreover, SpO<sub>2</sub> instability is poorly documented despite the correlation with an increased rate of late death or disability at 18&#xa0;months of age (<xref ref-type="bibr" rid="B18">Poets et al., 2015</xref>). SpO<sub>2</sub> instability increases when the infant is in supine compared to the prone position, independent of the mode of respiratory support (<xref ref-type="bibr" rid="B16">Miller-Barmak et al., 2020</xref>). A contributing factor for hypoxemic episodes in preterm infants includes sudden decrease in lung volume leading to small airway collapse and intrapulmonary shunt (<xref ref-type="bibr" rid="B3">Bolivar et al., 1995</xref>). We have previously shown in adults that SpO<sub>2</sub> was unstable when ventilation to perfusion ratio (V<sub>A</sub>/Q) was reduced, but was stable with increased shunt (<xref ref-type="bibr" rid="B13">Jones and Jones, 2000</xref>).</p>
<p>The concept of V<sub>A</sub>/Q is well established but infrequently used in neonatal practice, as the technics for measuring V<sub>A</sub>/Q in infants are technically difficult e.g. A-a nitrogen difference, or impossible e.g. Multiple Inert Gas Elimination Technique (MIGET) (<xref ref-type="bibr" rid="B5">Corbet et al., 1974</xref>; <xref ref-type="bibr" rid="B11">Hand et al., 1990</xref>; <xref ref-type="bibr" rid="B19">Roca and Wagner, 1994</xref>). Our non-invasive method for measuring V<sub>A</sub>/Q in adults is based on the different effects on SpO<sub>2</sub> of changing inspired oxygen pressure (<italic>P</italic>
<sub>I</sub>O<sub>2</sub>) when either V<sub>A</sub>/Q is reduced or shunt increased (<xref ref-type="bibr" rid="B13">Jones and Jones, 2000</xref>). The method depends on the shape of the oxyhaemoglobin dissociation curve (ODC). We adapted this method for preterm infants by using the neonatal rather than the adult ODC as the reference and plotting SpO<sub>2</sub> at different <italic>P</italic>
<sub>I</sub>O<sub>2</sub> (<xref ref-type="bibr" rid="B15">Lockwood et al., 2014</xref>). A computer algorithm derived a model of V<sub>A</sub>/Q and shunt from the paired SpO<sub>2</sub> and <italic>P</italic>
<sub>I</sub>O<sub>2</sub> dataset. Decreasing V<sub>A</sub>/Q is reflected by a right shift of the SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve. The right shift leads to a steeper slope of the curve as it intersects the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line (<italic>P</italic>
<sub>I</sub>O<sub>2</sub> &#x3d; F<sub>I</sub>O<sub>2</sub> x (barometric pressure&#x2014;saturated water vapour pressure), which means 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> is equal to room air at sea level (<xref ref-type="fig" rid="F1">Figure 1</xref>). Any paired SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> measurement located on the steep section of the ODC may predispose the infant to SpO<sub>2</sub> instability with small changes in alveolar oxygen.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>The neonatal oxyhaemoglobin dissociation curve is derived by plotting oxygen saturation (SaO<sub>2</sub>) against arterial oxygen pressure (PaO<sub>2</sub>) and determines the shape of the normal infant SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve (N), which is displaced to the right proportional to PCO<sub>2</sub>. Reducing Ventilation to Perfusion ratio (V<sub>A</sub>/Q) shifts the curve further to the right (dashed line). Its slope is now steep as it intersects the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line where a 1% change in <italic>P</italic>
<sub>I</sub>O<sub>2</sub> results in large SpO<sub>2</sub> instability. Shunt displaces the plateau downwards with trivial effect on stability.</p>
</caption>
<graphic xlink:href="fphys-14-1112115-g001.tif"/>
</fig>
<p>We hypothesised that SpO<sub>2</sub> instability in preterm infants with cardiorespiratory disease is not an epiphenomenon but an important clinical sign of a reduced V<sub>A</sub>/Q rather than right to left intrapulmonary shunt. We used a new algorithm and routine SpO<sub>2</sub> monitoring to derive V<sub>A</sub>/Q and shunt.</p>
</sec>
<sec sec-type="methods" id="s2">
<title>2 Methods</title>
<p>We conducted a prospective observational study in two phases. Firstly, we used a mathematical model of pulmonary gas exchange to examine changes in V<sub>A</sub>/Q or increasing shunt on the slope of the SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve.</p>
<p>Secondly, we recorded 24 hourly SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> measurements in extremely preterm infants requiring continuous SpO<sub>2</sub> monitoring. From these measurements, we derived V<sub>A</sub>/Q and shunt using a pulmonary gas exchange algorithm (<xref ref-type="bibr" rid="B15">Lockwood et al., 2014</xref>).</p>
<sec id="s2-1">
<title>2.1 Gas exchange model</title>
<p>We explored the effects of reducing V<sub>A</sub>/Q or increasing shunt on the arterial oxygen saturation (SaO<sub>2</sub>) vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve using a mathematical model of pulmonary gas exchange described by Olszowka and Wagner (<xref ref-type="bibr" rid="B17">Olszowka et al., 1980</xref>). Datasets were generated using the equations implemented on a spreadsheet supplied by Dr AJ Olszowka. The model allowed calculation of exact values of SaO<sub>2</sub> for a given <italic>P</italic>
<sub>I</sub>O<sub>2</sub>. The model lung is subdivided into three compartments: a shunt and two ventilated regions with different alveolar ventilation-perfusion ratios (V<sub>A</sub>/Q). The values of cardiac output, oxygen consumption, hemoglobin concentration, shunt fraction, and the distribution of blood flow and alveolar ventilation to the ventilated compartments can be set. The perfusion of one compartment was set at 90% of non-shunt flow while V<sub>A</sub>/Q was reduced stepwise from 0.85 to 0.3. Shunt was fixed at 2%, <italic>P</italic>
<sub>I</sub>O<sub>2</sub> was varied between 15&#xa0;kPa and 30&#xa0;kPa (F<sub>I</sub>O<sub>2</sub> &#x3d; 0.15&#x2013;0.3) and the SaO<sub>2</sub> was derived at the corresponding <italic>P</italic>
<sub>I</sub>O<sub>2</sub>. In the next step, the V<sub>A</sub>/Q was kept constant at 0.85 and the shunt was increased stepwise from 2%&#x2013;25%. The <italic>P</italic>
<sub>I</sub>O<sub>2</sub> was again varied between 15 and 30&#xa0;kPa (F<sub>I</sub>O<sub>2</sub> &#x3d; 0.15&#x2013;0.3) and the corresponding SpO<sub>2</sub> recorded. The slope of the SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve was calculated as it intersected the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line (F<sub>I</sub>O<sub>2</sub> &#x3d; 0.21).</p>
</sec>
<sec id="s2-2">
<title>2.2 Clinical study</title>
<sec id="s2-2-1">
<title>2.2.1 Study design</title>
<p>We conducted a prospective observational study at King Edward Memorial Hospital in Perth in Western Australia. The study was approved by the Women and Newborn Health Service Human Research Ethics Committee (HREC:1883EW and 20130193EW) in Perth.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Recruitment period, inclusion and exclusion criteria</title>
<p>Preterm infants born &#x2264;28&#xa0;weeks&#x2019; gestation without major congenital malformations were recruited from the Neonatal Clinical Care Unit at King Edward Memorial Hospital for Women in Perth, Western Australia (KEMH) between 21st August 2017 and the 1st February 2018. All included infants were part of the Preterm Infant Functional and Clinical Outcomes (PIFCO) study (ACTRN12613001062718). We started recruitment for this substudy based on the findings from the main PIFCO cohort, hence the shorter recruitment period and the much smaller number of infants included in the study (<xref ref-type="bibr" rid="B22">Svedenkrans et al., 2019</xref>). Informed consent was obtained from parents before the first measurement.</p>
</sec>
<sec id="s2-2-3">
<title>2.2.3 Conduct of study</title>
<p>Infants were assessed at 1&#xa0;week of age. SpO<sub>2</sub> measurements were recorded at hourly intervals for 24&#xa0;h (Masimo Infant Pulse Oximeter Adhesive Sensor RD SET<sup>&#xae;</sup> Inf; Philips Monitor. IntelliVue MP50 or MP70 Neonatal). Measurements were postponed until the following day in infants with changing respiratory support on the day of measurement. F<sub>I</sub>O<sub>2</sub> was adjusted by the bedside nurses to achieve a SpO<sub>2</sub> within the target range SpO<sub>2</sub> 90%&#x2013;94%. F<sub>I</sub>O<sub>2</sub> was later converted into <italic>P</italic>
<sub>I</sub>O<sub>2</sub> (<italic>P</italic>
<sub>I</sub>O<sub>2</sub> &#x3d; F<sub>I</sub>O<sub>2</sub> x (barometric pressure &#x2013; saturated water vapour pressure).</p>
</sec>
<sec id="s2-2-4">
<title>2.2.4 Analysis of results</title>
<p>The slope of the SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve in preterm infants was analysed using the pulmonary gas exchange algorithm with three lung compartments (<xref ref-type="bibr" rid="B15">Lockwood et al., 2014</xref>). Reference normative data of V<sub>A</sub>/Q, right shift of the oxyhemoglobin dissociation curve and right to left shunt in healthy term infants studied in the first week of life were used to quantify the magnitude of the abnormalities in our population of extremely preterm infants (<xref ref-type="bibr" rid="B7">Dassios et al., 2017</xref>; <xref ref-type="bibr" rid="B6">Dassios et al., 2019</xref>).</p>
</sec>
<sec id="s2-2-5">
<title>2.2.5 Background data</title>
<p>Patient data information including duration of respiratory support and oxygen therapy were collected from the medical charts or the discharge summaries.</p>
</sec>
<sec id="s2-2-6">
<title>2.2.6 Statistical analyses</title>
<p>Study data were collected and managed using Research Electronic Data Capture (REDCap) software hosted at The University of Western Australia. REDCap is a secure, web-based application designed to support data capture for research studies (<xref ref-type="bibr" rid="B12">Harris et al., 2009</xref>).</p>
<p>Parametric data are reported as mean and standard deviation (SD) and non-parametric data as median and variance. Statistical analyses included Student&#x2019;s <italic>t</italic>-test for the comparison of parametric and Mann-Whitney-U test for the comparison of non-parametric data. Data were analysed within SPSS (v25&#xb7;0&#xb7;0&#xb7;1; IBM Corp, United States).</p>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Gas exchange model</title>
<p>Reducing V<sub>A</sub>/Q shifted the SaO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve to the right so that its slope increased as it intercepted the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line (<xref ref-type="fig" rid="F2">Figure 2A</xref>). SaO<sub>2</sub> at the intercept fell in increasingly large increments particularly with V<sub>A</sub>/Q below 0.5: e.g., a 24% SaO<sub>2</sub> fall when V<sub>A</sub>/Q dropped from 0.4 to 0.3. The effect of V<sub>A</sub>/Q on the slope of the curve (<xref ref-type="fig" rid="F2">Figure 2B</xref>) suggested an unstable SaO<sub>2</sub> once V<sub>A</sub>/Q fell below 0.5: a 1% change in <italic>P</italic>
<sub>I</sub>O<sub>2</sub> at 0.3 V<sub>A</sub>/Q resulted in a &#x3e;8% change in SaO<sub>2</sub>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>The gas exchange model derived <italic>P</italic>
<sub>I</sub>O<sub>2</sub> vs. SaO<sub>2</sub> curves when V<sub>A</sub>/Q was reduced from 0.85 to 0.3. (Figure 2A). The slope of the curve as it intercepted the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line is shown for V<sub>A</sub>/Q (Figure 2B). The <italic>P</italic>
<sub>I</sub>O<sub>2</sub> vs. SaO<sub>2</sub> curves when shunt increased from 2%&#x2013;25% (Figure 2C) and their slopes (Figure 2D).</p>
</caption>
<graphic xlink:href="fphys-14-1112115-g002.tif"/>
</fig>
<p>In contrast, increasing shunt (<xref ref-type="fig" rid="F2">Figure 2C</xref>) displaced the curve downwards but even large changes in <italic>P</italic>
<sub>I</sub>O<sub>2</sub> caused only small changes in both SaO<sub>2</sub> and slope (<xref ref-type="fig" rid="F2">Figure 2D</xref>): for example in the setting of a 20% intrapulmonary shunt, a 1% change in <italic>P</italic>
<sub>I</sub>O<sub>2</sub> resulted in &#x3c;1% change in SaO<sub>2</sub>.</p>
</sec>
<sec id="s3-2">
<title>3.2 Clinical study</title>
<p>We studied 31 extremely preterm infants (mean (SD) 26.2 (1.0)&#xa0;weeks&#x2019; gestation) at median (IQR) 1.1 (0.5)&#xa0;weeks&#x2019; postnatal age. Eight infants were mechanically ventilated and 23 infants were on non-invasive respiratory support during the measurements. Demographics for this infant cohort are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Demographics of infants studied (<italic>n</italic> &#x3d; 31).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Male (n, %)</th>
<th align="left">20 (64.5)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">GA (w)</td>
<td align="left">26.2 &#xb1; 1.0</td>
</tr>
<tr>
<td align="left">Age at Test (w)</td>
<td align="left">1.1 (0.5)</td>
</tr>
<tr>
<td align="left">Weight at birth (g)</td>
<td align="left">890 &#xb1; 196</td>
</tr>
<tr>
<td align="left">Mechanical Ventilation at test (n, %)</td>
<td align="left">8 (26.0)</td>
</tr>
<tr>
<td align="left">Non-invasive respiratory support at test (n, %)</td>
<td align="left">23 (74.0)</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>GA, gestational age; w, week; g, gram. Values are reported as mean &#xb1; SD, median (IQR) or n (%).</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Ten infants received 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> throughout the 24&#xa0;h monitoring period. Median (IQR) <italic>P</italic>
<sub>I</sub>O<sub>2</sub> in the remainder was 28 (<xref ref-type="bibr" rid="B13">Jones and Jones, 2000</xref>) % to maintain target SpO<sub>2</sub> 90%&#x2013;94%. The number of consecutive hourly readings at the same <italic>P</italic>
<sub>I</sub>O<sub>2</sub> in each infant are shown in <xref ref-type="table" rid="T2">Table 2</xref>. There were sufficient data to describe the steep part of the SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve in all infants in terms of V<sub>A</sub>/Q and shift, but complete characterization of the top end of the curve was constrained by the upper limit of the SpO<sub>2</sub> target range and therefore restricted <italic>P</italic>
<sub>I</sub>O<sub>2</sub>. The V<sub>A</sub>/Q was less than 0.74 in every infant and only two infants (Nr. 3 and 6) had a large shunt resolving after a few hours into predominantly low V<sub>A</sub>/Q compartments. The median (IQR) SpO<sub>2</sub> of all preterm infants was 8.0 (7.0) % at a constant <italic>P</italic>
<sub>I</sub>O<sub>2</sub> over at least 5&#xa0;h. At a V<sub>A</sub>/Q &#x3c; 0.52, in 10 out of 31 infants the variability of the SpO<sub>2</sub> ranged from 11%&#x2013;17%, whereas in the four infants with V<sub>A</sub>/Q &#x3e; 0.52, the range of SpO<sub>2</sub> recorded was &#x2264;4% (V<sub>A</sub>/Q 0.73, 0.56, 0.57, 0.54). Seven infants had a profound desaturation with a SpO<sub>2</sub> &#x2264; 84%. All seven infants had a V<sub>A</sub>/Q &#x2264; 0.49 and four of them V<sub>A</sub>/Q &#x3c; 0.31. Examples of the pattern of SpO<sub>2</sub> in three infants are shown in <xref ref-type="fig" rid="F3">Figure 3</xref> using the gas exchange algorithm (<xref ref-type="bibr" rid="B15">Lockwood et al., 2014</xref>) and a format similar to <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Measurements of included infants.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Infant number</th>
<th align="center">MV</th>
<th align="center">V<sub>A</sub>/Q</th>
<th align="center">Shift kPa</th>
<th align="center">
<italic>P</italic>
<sub>I</sub>O<sub>2</sub>kPa</th>
<th align="center">Number of SpO<sub>2</sub> at that <italic>P</italic>
<sub>I</sub>O<sub>2</sub>
</th>
<th align="center">Min SpO<sub>2</sub>%</th>
<th align="center">Med SpO<sub>2</sub>%</th>
<th align="center">Range SpO<sub>2</sub>%</th>
<th align="center">Variance SpO<sub>2</sub>%</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">1</td>
<td align="center">N</td>
<td align="center">0.52</td>
<td align="center">11.6</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">93</td>
<td align="center">98</td>
<td align="center">14</td>
<td align="center">6.8</td>
</tr>
<tr>
<td align="center">2</td>
<td align="center">N</td>
<td align="center">0.49</td>
<td align="center">12.4</td>
<td align="center">22</td>
<td align="center">15</td>
<td align="center">83</td>
<td align="center">94</td>
<td align="center">15</td>
<td align="center">11.0</td>
</tr>
<tr>
<td align="center">3</td>
<td align="center">Y</td>
<td align="center">0.22</td>
<td align="center">27.5</td>
<td align="center">35</td>
<td align="center">14</td>
<td align="center">90</td>
<td align="center">94</td>
<td align="center">8</td>
<td align="center">6.0</td>
</tr>
<tr>
<td align="center">4</td>
<td align="center">N</td>
<td align="center">0.21</td>
<td align="center">28.8</td>
<td align="center">38</td>
<td align="center">5</td>
<td align="center">88</td>
<td align="center">95</td>
<td align="center">9</td>
<td align="center">28.8</td>
</tr>
<tr>
<td align="center">5</td>
<td align="center">N</td>
<td align="center">0.27</td>
<td align="center">22.8</td>
<td align="center">30</td>
<td align="center">14</td>
<td align="center">87</td>
<td align="center">93</td>
<td align="center">9</td>
<td align="center">5.9</td>
</tr>
<tr>
<td align="center">6</td>
<td align="center">Y</td>
<td align="center">0.31</td>
<td align="center">19.2</td>
<td align="center">27</td>
<td align="center">13</td>
<td align="center">84</td>
<td align="center">96</td>
<td align="center">13</td>
<td align="center">16.9</td>
</tr>
<tr>
<td align="center">7</td>
<td align="center">N</td>
<td align="center">0.5</td>
<td align="center">12.0</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">93</td>
<td align="center">97</td>
<td align="center">6</td>
<td align="center">3.3</td>
</tr>
<tr>
<td align="center">8</td>
<td align="center">N</td>
<td align="center">0.73</td>
<td align="center">8.3</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">96</td>
<td align="center">98</td>
<td align="center">4</td>
<td align="center">1.5</td>
</tr>
<tr>
<td align="center">9</td>
<td align="center">Y</td>
<td align="center">0.46</td>
<td align="center">13.3</td>
<td align="center">21</td>
<td align="center">19</td>
<td align="center">91</td>
<td align="center">95</td>
<td align="center">6</td>
<td align="center">3.3</td>
</tr>
<tr>
<td align="center">10</td>
<td align="center">N</td>
<td align="center">0.56</td>
<td align="center">10.8</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">97</td>
<td align="center">98</td>
<td align="center">3</td>
<td align="center">0.6</td>
</tr>
<tr>
<td align="center">11</td>
<td align="center">N</td>
<td align="center">0.57</td>
<td align="center">10.6</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">96</td>
<td align="center">99</td>
<td align="center">4</td>
<td align="center">1.7</td>
</tr>
<tr>
<td align="center">12</td>
<td align="center">N</td>
<td align="center">0.54</td>
<td align="center">11.3</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">96</td>
<td align="center">98</td>
<td align="center">4</td>
<td align="center">1.6</td>
</tr>
<tr>
<td align="center">13</td>
<td align="center">Y</td>
<td align="center">0.35</td>
<td align="center">17.5</td>
<td align="center">25</td>
<td align="center">15</td>
<td align="center">91</td>
<td align="center">94</td>
<td align="center">5</td>
<td align="center">2.1</td>
</tr>
<tr>
<td align="center">14</td>
<td align="center">N</td>
<td align="center">0.3</td>
<td align="center">20.0</td>
<td align="center">28</td>
<td align="center">17</td>
<td align="center">90</td>
<td align="center">96</td>
<td align="center">8</td>
<td align="center">4.0</td>
</tr>
<tr>
<td align="center">15</td>
<td align="center">N</td>
<td align="center">0.39</td>
<td align="center">15.4</td>
<td align="center">24</td>
<td align="center">8</td>
<td align="center">84</td>
<td align="center">91</td>
<td align="center">14</td>
<td align="center">22.3</td>
</tr>
<tr>
<td align="center">16</td>
<td align="center">N</td>
<td align="center">0.49</td>
<td align="center">12.3</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">93</td>
<td align="center">96</td>
<td align="center">7</td>
<td align="center">2.8</td>
</tr>
<tr>
<td align="center">17</td>
<td align="center">N</td>
<td align="center">0.3</td>
<td align="center">20.2</td>
<td align="center">28</td>
<td align="center">16</td>
<td align="center">91</td>
<td align="center">95</td>
<td align="center">5</td>
<td align="center">2.7</td>
</tr>
<tr>
<td align="center">18</td>
<td align="center">N</td>
<td align="center">0.49</td>
<td align="center">12.4</td>
<td align="center">21</td>
<td align="center">19</td>
<td align="center">92</td>
<td align="center">96</td>
<td align="center">7</td>
<td align="center">4.1</td>
</tr>
<tr>
<td align="center">19</td>
<td align="center">Y</td>
<td align="center">0.24</td>
<td align="center">25.0</td>
<td align="center">32</td>
<td align="center">10</td>
<td align="center">84</td>
<td align="center">93</td>
<td align="center">14</td>
<td align="center">17.5</td>
</tr>
<tr>
<td align="center">20</td>
<td align="center">Y</td>
<td align="center">0.27</td>
<td align="center">22.6</td>
<td align="center">30</td>
<td align="center">12</td>
<td align="center">88</td>
<td align="center">94</td>
<td align="center">8</td>
<td align="center">7.7</td>
</tr>
<tr>
<td align="center">21</td>
<td align="center">N</td>
<td align="center">0.28</td>
<td align="center">22.7</td>
<td align="center">30</td>
<td align="center">12</td>
<td align="center">88</td>
<td align="center">93.5</td>
<td align="center">9</td>
<td align="center">8.4</td>
</tr>
<tr>
<td align="center">22</td>
<td align="center">Y</td>
<td align="center">0.3</td>
<td align="center">20.1</td>
<td align="center">30</td>
<td align="center">9</td>
<td align="center">89</td>
<td align="center">95</td>
<td align="center">8</td>
<td align="center">6.3</td>
</tr>
<tr>
<td align="center">23</td>
<td align="center">N</td>
<td align="center">0.44</td>
<td align="center">13.8</td>
<td align="center">21</td>
<td align="center">14</td>
<td align="center">81</td>
<td align="center">93</td>
<td align="center">16</td>
<td align="center">16.8</td>
</tr>
<tr>
<td align="center">24</td>
<td align="center">N</td>
<td align="center">0.27</td>
<td align="center">22.8</td>
<td align="center">30</td>
<td align="center">10</td>
<td align="center">81</td>
<td align="center">94</td>
<td align="center">17</td>
<td align="center">22</td>
</tr>
<tr>
<td align="center">25</td>
<td align="center">N</td>
<td align="center">0.47</td>
<td align="center">12.8</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">91</td>
<td align="center">95</td>
<td align="center">8</td>
<td align="center">3.7</td>
</tr>
<tr>
<td align="center">26</td>
<td align="center">N</td>
<td align="center">0.52</td>
<td align="center">11.6</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">95</td>
<td align="center">97.5</td>
<td align="center">4</td>
<td align="center">11.6</td>
</tr>
<tr>
<td align="center">27</td>
<td align="center">N</td>
<td align="center">0.35</td>
<td align="center">17.3</td>
<td align="center">25</td>
<td align="center">13</td>
<td align="center">92</td>
<td align="center">94</td>
<td align="center">6</td>
<td align="center">17.3</td>
</tr>
<tr>
<td align="center">28</td>
<td align="center">N</td>
<td align="center">0.39</td>
<td align="center">15.5</td>
<td align="center">23</td>
<td align="center">10</td>
<td align="center">86</td>
<td align="center">94</td>
<td align="center">12</td>
<td align="center">15.5</td>
</tr>
<tr>
<td align="center">29</td>
<td align="center">N</td>
<td align="center">0.47</td>
<td align="center">12.8</td>
<td align="center">21</td>
<td align="center">24</td>
<td align="center">92</td>
<td align="center">96</td>
<td align="center">5</td>
<td align="center">12.6</td>
</tr>
<tr>
<td align="center">30</td>
<td align="center">N</td>
<td align="center">0.27</td>
<td align="center">22.8</td>
<td align="center">30</td>
<td align="center">13</td>
<td align="center">84</td>
<td align="center">93</td>
<td align="center">12</td>
<td align="center">22.8</td>
</tr>
<tr>
<td align="center">31</td>
<td align="center">Y</td>
<td align="center">0.29</td>
<td align="center">20.8</td>
<td align="center">28</td>
<td align="center">9</td>
<td align="center">87</td>
<td align="center">94</td>
<td align="center">11</td>
<td align="center">20.8</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Spot reading of SpO<sub>2</sub> and <italic>P</italic>
<sub>I</sub>O<sub>2</sub> at 1&#xa0;h intervals within a 24&#xa0;h monitoring period. Infant 3 had four readings at <italic>P</italic>
<sub>I</sub>O<sub>2</sub> 78, 92, 92 &#x26; 91&#xa0;kPa from 3&#x2013;6&#xa0;h to derive a 27% shunt, V<sub>A</sub>/Q 0.24 and shift of 34.1&#xa0;kPa. From 10&#x2013;24&#xa0;h the 27% shunt resolved at constant 35&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> to a single compartment V<sub>A</sub>/Q of 0.22 and 27.4&#xa0;kPa shift. MV, mechanical ventilation.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>The normal <italic>P</italic>
<sub>I</sub>O<sub>2</sub> vs. SpO<sub>2</sub> curve is shown on the left. SpO<sub>2</sub> readings (closed boxes) with the calculated gas exchange model curves for infants 3, 6 and 10. Infant 3 had 27% shunt and V<sub>A</sub>/Q of 0.24 (3a), which resolved during the day to a single V<sub>A</sub>/Q compartment of 0.22 and a SpO<sub>2</sub> range of 8% (3b) and subsequently into a final point (open box) at the end of the day to almost overlying the V<sub>A</sub>/Q 0.31 line of infant 6 (the changes are indicated by the dashed arrow). Infant 6 initially had a 21% shunt and V<sub>A</sub>/Q of 0.43 (not shown in Figure), which resolved into a V<sub>A</sub>/Q of 0.31 in the 9&#x2013;24&#xa0;h study period. Infant 10 had a V<sub>A</sub>/Q of 0.56 and a SpO<sub>2</sub> range of 3% from 24 SpO<sub>2</sub> values.</p>
</caption>
<graphic xlink:href="fphys-14-1112115-g003.tif"/>
</fig>
<p>One infant (Nr. 3) required a <italic>P</italic>
<sub>I</sub>O<sub>2</sub> &#x3e; 50&#xa0;kPa to maintain SpO<sub>2</sub> within the target range. This infant had an initial large shunt and greatly reduced V<sub>A</sub>/Q, which later resolved into a single low V<sub>A</sub>/Q compartment. Infant Nr. 3, 6 and 10 are described in detail in the legend to <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<p>The 10 infants exhibiting the most variable SpO<sub>2</sub> (range &#x3e;10% in SpO<sub>2</sub>) were equally distributed between those breathing spontaneously (<italic>n</italic> &#x3d; 5) and those requiring mechanical ventilation (<italic>n</italic> &#x3d; 5) (<xref ref-type="table" rid="T1">Table 1</xref>). The increased variance in SpO<sub>2</sub> as V<sub>A</sub>/Q falls below 0.5 is illustrated in <xref ref-type="fig" rid="F4">Figure 4A</xref> shift increases &#x3e;11&#xa0;kPa in <xref ref-type="fig" rid="F4">Figure 4B</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The fall in V<sub>A</sub>/Q below 0.52 (Figure 4A) and an increase in shift &#x3e;11&#xa0;kPa (Figure 4B), were associated with an abrupt widening of the variance of SpO<sub>2</sub>. The two uppermost points have a range of 16% and 17% SpO<sub>2</sub> in infants 23 and 24. The normal values for V<sub>A</sub>/Q and shift are shown by the vertical line. The normal variance for SpO<sub>2</sub> is 2.2 shown by the contrasting boxes on the vertical line (<xref ref-type="bibr" rid="B7">Dassios et al., 2017</xref>; <xref ref-type="bibr" rid="B6">Dassios et al., 2019</xref>).</p>
</caption>
<graphic xlink:href="fphys-14-1112115-g004.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>We suggest that a reduced V<sub>A</sub>/Q predisposes preterm infants to SpO<sub>2</sub> instability. The pulmonary gas exchange model showed that reducing V<sub>A</sub>/Q shifted the <italic>P</italic>
<sub>I</sub>O<sub>2</sub> vs. SaO<sub>2</sub> curve to the right, whereas shunt displaced the curve downwards. Right shift increased the effective curve slope up to 16 times where it intercepted the 21&#xa0;kPa <italic>P</italic>
<sub>I</sub>O<sub>2</sub> line e.g., breathing room air. The right shift predicted large fluctuations in SaO<sub>2</sub> with small changes in alveolar oxygen. Reducing V<sub>A</sub>/Q from 0.4 to 0.3 predicted a 24% fall in SaO<sub>2</sub> breathing air. In contrast a large shunt had little effect on curve slope and a small fall in SaO<sub>2</sub>.</p>
<p>These predictions were confirmed in extremely preterm infants by routine monitoring of SpO<sub>2</sub> and <italic>P</italic>
<sub>I</sub>O<sub>2</sub>. The paired SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> measurements were used to derive V<sub>A</sub>/Q, right shift and shunt using Lockwood&#x2019;s algorithm (<xref ref-type="bibr" rid="B15">Lockwood et al., 2014</xref>). All infants had a reduced V<sub>A</sub>/Q. The SpO<sub>2</sub> was unstable with a median (IQR) 8.0 (7.0) % at a constant <italic>P</italic>
<sub>I</sub>O<sub>2</sub> compared to a mean (&#xb1;SD) of 3% (&#xb1;1.5%) in healthy newborn infants breathing room air (<xref ref-type="bibr" rid="B7">Dassios et al., 2017</xref>; <xref ref-type="bibr" rid="B6">Dassios et al., 2019</xref>). Infants with V<sub>A</sub>/Q &#x3c; 0.52 had considerable increase in SpO<sub>2</sub> variability, seven infants had episodes with SpO<sub>2</sub> &#x2264; 84%, whilst four of these infants were dependent on supplemental oxygen &#x2265;27%. Two infants had transient shunt lasting for a few hours, which was not associated with SpO<sub>2</sub> instability. In the remaining infants, the <italic>P</italic>
<sub>I</sub>O<sub>2</sub> was insufficiently high to fully characterize shunt due to the SpO<sub>2</sub> target range (90%&#x2013;94%).</p>
<p>The non-invasive method for deriving V<sub>A</sub>/Q and shunt normally depends on varying <italic>P</italic>
<sub>I</sub>O<sub>2</sub> (<xref ref-type="bibr" rid="B20">Rowe et al., 2010</xref>; <xref ref-type="bibr" rid="B1">Bamat et al., 2015</xref>; <xref ref-type="bibr" rid="B7">Dassios et al., 2017</xref>; <xref ref-type="bibr" rid="B6">Dassios et al., 2019</xref>; <xref ref-type="bibr" rid="B22">Svedenkrans et al., 2019</xref>; <xref ref-type="bibr" rid="B2">Bamat et al., 2022</xref>). More recently, we showed that V<sub>A</sub>/Q and/or shift can be derived from serial SpO<sub>2</sub> measurements at a fixed <italic>P</italic>
<sub>I</sub>O<sub>2</sub> (<xref ref-type="bibr" rid="B22">Svedenkrans et al., 2019</xref>; <xref ref-type="bibr" rid="B21">Stoecklin et al., 2021</xref>). A reduced V<sub>A</sub>/Q lowers SpO<sub>2</sub> at any given <italic>P</italic>
<sub>I</sub>O<sub>2</sub> and increases the effective slope of the SpO<sub>2</sub> vs. <italic>P</italic>
<sub>I</sub>O<sub>2</sub> curve, but does not by itself make SpO<sub>2</sub> unstable unless V<sub>A</sub>/Q itself is unstable (<xref ref-type="bibr" rid="B14">Jones, 2021</xref>). The reduced V<sub>A</sub>/Q amplifies the destabilizing effects on SpO<sub>2</sub> during changes in cardiac output, ventilation, central or obstructive apnea and posture. Interestingly, the 10 infants with the most unstable SpO<sub>2</sub> were equally distributed between those breathing spontaneously and those on mechanical ventilation. Reduced V<sub>A</sub>/Q can lead to SpO<sub>2</sub> instability independent of the mode of respiratory support used in an infant.</p>
<p>In agreement with our results, a reduced V<sub>A</sub>/Q in adults leads to SpO<sub>2</sub> instability breathing room air. By adding a time dimension, the different effects of V<sub>A</sub>/Q or shunt in individual adult patients breathing air is more clearly illustrated by dynamic waterfall plots of continuously monitored SpO<sub>2</sub> (<xref ref-type="bibr" rid="B10">Entwistle et al., 1991</xref>; <xref ref-type="bibr" rid="B13">Jones and Jones, 2000</xref>; <xref ref-type="bibr" rid="B14">Jones, 2021</xref>) than do &#x201c;static&#x201d; histograms (<xref ref-type="bibr" rid="B4">Borenstein-Levin et al., 2020</xref>). A typical example of unstable SpO2 in a patient with reduced V<sub>A</sub>/Q breathing air is shown in <xref ref-type="fig" rid="F3">Figure 3</xref> in Ref (<xref ref-type="bibr" rid="B14">Jones, 2021</xref>). This patient had blunt curves with SpO<sub>2</sub> fluctuating from 77% to 95%. In contrast, a patient with increased shunt had superimposed stable SpO<sub>2</sub> peaks within a much narrower SpO<sub>2</sub> range.</p>
<p>We note that the right to left shunt in our study might include some element of cardiac shunting occurring <italic>via</italic> an open ductus arteriosus. Our method does not allow differentiation between intrapulmonary or cardiac shunt. However, in the majority of preterm infants the ductus arteriosus would have functionally closed on day seven of life. In the few infants with a persistent ductus arteriosus, at day seven of life a predominantly left to right shunt would be expected which would not affect our calculations (<xref ref-type="bibr" rid="B8">de Klerk et al., 2020</xref>).</p>
<p>The clinical applicability of our findings is that routine measurements of SpO<sub>2</sub> and <italic>P</italic>
<sub>I</sub>O<sub>2</sub> can be analysed with our algorithm to derive V<sub>A</sub>/Q and shunt, thereby characterizing respiratory disease in preterm infants. V<sub>A</sub>/Q and shunt incorporate information on the mechanisms of hypoxemia and explain SpO<sub>2</sub> instability. Consequently, SpO2 instability can be prevented by a small increase in <italic>P</italic>
<sub>I</sub>O<sub>2</sub>.</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>In conclusion, we reported that predisposition to oxygen saturation instability in preterm infants results from a reduced ventilation to perfusion ratio rather than from increased intrapulmonary shunt. We have highlighted how routine monitoring can be used to derive non-invasive measurements of oxygenation impairment in preterm infants. Future research should aim at strategies to not only detect SpO<sub>2</sub> instability, but to also prevent such instability.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Ethics statement</title>
<p>The study was approved by the Women and Newborn Health Service Human Research Ethics Committee (HREC:1883EW and 20130193EW) in Perth. Written informed consent to participate in this study was provided by the participants&#x2019; legal guardian/next of kin.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>BS collected and analysed the data, interpreted the data, assisted with the REDCap database design and development, drafted the initial manuscript together with JGJ, performed literature search, and approved the final manuscript as submitted. YC collected the data, reviewed and revised the manuscript, and approved the final manuscript as submitted. JGJ developed the algorithms used for the calculation of shunt, shift and VA/Q, drafted reviewed and revised the manuscript, drafted the Figures and approved the final manuscript as submitted. GL developed the algorithms used for the calculation of shunt, shift and VA/Q, reviewed and revised the manuscript. TD drafted the initial manuscript, revised and approved the final manuscript. JP was the principal investigator obtaining funding, leading study design including development of the REDCap database, verified all statistical calculations, interpreted the data, critically reviewed and revised the manuscript, and approved the final manuscript as submitted.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>All phases of this study were supported by the University of Western Australia and the Women and Newborn Health Service of Western Australia. Funded by National Health and Medical Research Council (NHMRC) of Australia (GNT1047689, GNT1057514) and the Metropolitan Health Research Infrastructure Fund (MHRIF). BS was supported by the Swiss National Science Foundation (P2BSP3_158837) and a Research Training Program scholarship from The University of Western Australia. JP was supported by a NHMRC Fellowships (RF1077691, GNT1196188).</p>
</sec>
<ack>
<p>The authors thank Dr A. Olszowka, Department of Physiology, University of Buffalo, New York, United States for providing his pulmonary gas exchange program.</p>
</ack>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bamat</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Ghavam</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>DeMauro</surname>
<given-names>S. B.</given-names>
</name>
<name>
<surname>Jensen</surname>
<given-names>E. A.</given-names>
</name>
<name>
<surname>Roberts</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Reliability of a noninvasive measure of V./Q. Mismatch for bronchopulmonary dysplasia</article-title>. <source>Ann. Am. Thorac. Soc.</source> <volume>12</volume> (<issue>5</issue>), <fpage>727</fpage>&#x2013;<lpage>733</lpage>. <pub-id pub-id-type="doi">10.1513/AnnalsATS.201410-462OC</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Bamat</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Orians</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>Abbasi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Morley</surname>
<given-names>C. J.</given-names>
</name>
<name>
<surname>Ross Russell</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Panitch</surname>
<given-names>H. B.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <source>Use of ventilation/perfusion mismatch to guide individualised CPAP level selection in preterm infants: A feasibility trial</source>. <publisher-name>Arch Dis Child Fetal Neonatal Ed</publisher-name>. <pub-id pub-id-type="doi">10.1136/archdischild-2022-324474</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolivar</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Gerhardt</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Gonzalez</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hummler</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Claure</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Everett</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>1995</year>). <article-title>Mechanisms for episodes of hypoxemia in preterm infants undergoing mechanical ventilation</article-title>. <source>J. Pediatr.</source> <volume>127</volume> (<issue>5</issue>), <fpage>767</fpage>&#x2013;<lpage>773</lpage>. <pub-id pub-id-type="doi">10.1016/s0022-3476(95)70171-0</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Borenstein-Levin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Konikoff</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Solimano</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Clinical quantification of SpO2 instability using a new histogram classification system: A clinical study</article-title>. <source>Pediatr. Res.</source> <volume>87</volume> (<issue>4</issue>), <fpage>716</fpage>&#x2013;<lpage>720</lpage>. <pub-id pub-id-type="doi">10.1038/s41390-019-0566-6</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Corbet</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Ross</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Beaudry</surname>
<given-names>P. H.</given-names>
</name>
<name>
<surname>Stern</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>1974</year>). <article-title>Ventilation-perfusion relationships as assessed by a ADN2 in hyaline membrane disease</article-title>. <source>J. Appl. Physiol.</source> <volume>36</volume> (<issue>1</issue>), <fpage>74</fpage>&#x2013;<lpage>81</lpage>. <pub-id pub-id-type="doi">10.1152/jappl.1974.36.1.74</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dassios</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Rossor</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Greenough</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Using the fetal oxyhaemoglobin dissociation curve to calculate the ventilation/perfusion ratio and right to left shunt in healthy newborn infants</article-title>. <source>J. Clin. Monit. Comput.</source> <volume>33</volume> (<issue>3</issue>), <fpage>545</fpage>&#x2013;<lpage>546</lpage>. <pub-id pub-id-type="doi">10.1007/s10877-018-0168-6</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dassios</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Rossor</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Greenough</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Ventilation/perfusion ratio and right to left shunt in healthy newborn infants</article-title>. <source>J. Clin. Monit. Comput.</source> <volume>31</volume> (<issue>6</issue>), <fpage>1229</fpage>&#x2013;<lpage>1234</lpage>. <pub-id pub-id-type="doi">10.1007/s10877-016-9969-7</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>de Klerk</surname>
<given-names>J. C. A.</given-names>
</name>
<name>
<surname>Engbers</surname>
<given-names>A. G. J.</given-names>
</name>
<name>
<surname>van Beek</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Flint</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>Reiss</surname>
<given-names>I. K. M.</given-names>
</name>
<name>
<surname>Voller</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Spontaneous closure of the ductus arteriosus in preterm infants: A systematic review</article-title>. <source>Front. Pediatr.</source> <volume>8</volume>, <fpage>541</fpage>. <pub-id pub-id-type="doi">10.3389/fped.2020.00541</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Di Fiore</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Bloom</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Orge</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Schutt</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Schluchter</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Cheruvu</surname>
<given-names>V. K.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>A higher incidence of intermittent hypoxemic episodes is associated with severe retinopathy of prematurity</article-title>. <source>J. Pediatr.</source> <volume>157</volume> (<issue>1</issue>), <fpage>69</fpage>&#x2013;<lpage>73</lpage>. <pub-id pub-id-type="doi">10.1016/j.jpeds.2010.01.046</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Entwistle</surname>
<given-names>M. D.</given-names>
</name>
<name>
<surname>Roe</surname>
<given-names>P. G.</given-names>
</name>
<name>
<surname>Sapsford</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Berrisford</surname>
<given-names>R. G.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Patterns of oxygenation after thoracotomy</article-title>. <source>Br. J. Anaesth.</source> <volume>67</volume> (<issue>6</issue>), <fpage>704</fpage>&#x2013;<lpage>711</lpage>. <pub-id pub-id-type="doi">10.1093/bja/67.6.704</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hand</surname>
<given-names>I. L.</given-names>
</name>
<name>
<surname>Shepard</surname>
<given-names>E. K.</given-names>
</name>
<name>
<surname>Krauss</surname>
<given-names>A. N.</given-names>
</name>
<name>
<surname>Auld</surname>
<given-names>P. A.</given-names>
</name>
</person-group> (<year>1990</year>). <article-title>Ventilation-perfusion abnormalities in the preterm infant with hyaline membrane disease: A two-compartment model of the neonatal lung</article-title>. <source>Pediatr. Pulmonol.</source> <volume>9</volume> (<issue>4</issue>), <fpage>206</fpage>&#x2013;<lpage>213</lpage>. <pub-id pub-id-type="doi">10.1002/ppul.1950090404</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harris</surname>
<given-names>P. A.</given-names>
</name>
<name>
<surname>Taylor</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Thielke</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Payne</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Gonzalez</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Conde</surname>
<given-names>J. G.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Research electronic data capture (REDCap)&#x2014;a metadata-driven methodology and workflow process for providing translational research informatics support</article-title>. <source>J. Biomed. Inf.</source> <volume>42</volume> (<issue>2</issue>), <fpage>377</fpage>&#x2013;<lpage>381</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbi.2008.08.010</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>S. E.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Discriminating between the effect of shunt and reduced VA/Q on arterial oxygen saturation is particularly useful in clinical practice</article-title>. <source>J. Clin. Monit. Comput.</source> <volume>16</volume> (<issue>5-6</issue>), <fpage>337</fpage>&#x2013;<lpage>350</lpage>. <pub-id pub-id-type="doi">10.1023/a:1011495416005</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
</person-group> (<year>2021</year>). <source>Oxygen saturation instability in suspected covid-19 patients; contrasting effects of reduced VA/Q and shunt</source>. <publisher-name>medRxiv</publisher-name>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lockwood</surname>
<given-names>G. G.</given-names>
</name>
<name>
<surname>Fung</surname>
<given-names>N. L.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Evaluation of a computer program for non-invasive determination of pulmonary shunt and ventilation-perfusion mismatch</article-title>. <source>J. Clin. Monit. Comput.</source> <volume>28</volume> (<issue>6</issue>), <fpage>581</fpage>&#x2013;<lpage>590</lpage>. <pub-id pub-id-type="doi">10.1007/s10877-014-9554-x</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miller-Barmak</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Riskin</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hochwald</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Haddad</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Dinur</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Vortman</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Oxygenation instability assessed by oxygen saturation histograms during supine vs prone position in very low birthweight infants receiving noninvasive respiratory support</article-title>. <source>J. Pediatr.</source> <volume>226</volume>, <fpage>123</fpage>&#x2013;<lpage>128</lpage>. <pub-id pub-id-type="doi">10.1016/j.jpeds.2020.06.066</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Olszowka</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Wagner</surname>
<given-names>P. D.</given-names>
</name>
</person-group> (<year>1980</year>). &#x201c;<article-title>8 - numerical analysis of gas exchange</article-title>,&#x201d; in <source>Ventilation, blood flow, and diffusion</source>. Editor <person-group person-group-type="editor">
<name>
<surname>West</surname>
<given-names>J. B.</given-names>
</name>
</person-group> (<publisher-name>Academic Press</publisher-name>), <fpage>263</fpage>&#x2013;<lpage>306</lpage>.</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Poets</surname>
<given-names>C. F.</given-names>
</name>
<name>
<surname>Roberts</surname>
<given-names>R. S.</given-names>
</name>
<name>
<surname>Schmidt</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Whyte</surname>
<given-names>R. K.</given-names>
</name>
<name>
<surname>Asztalos</surname>
<given-names>E. V.</given-names>
</name>
<name>
<surname>Bader</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Association between intermittent hypoxemia or bradycardia and late death or disability in extremely preterm infants</article-title>. <source>Jama</source> <volume>314</volume> (<issue>6</issue>), <fpage>595</fpage>&#x2013;<lpage>603</lpage>. <pub-id pub-id-type="doi">10.1001/jama.2015.8841</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roca</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wagner</surname>
<given-names>P. D.</given-names>
</name>
</person-group> (<year>1994</year>). <article-title>Contribution of multiple inert gas elimination technique to pulmonary medicine. 1. Principles and information content of the multiple inert gas elimination technique</article-title>. <source>Thorax</source> <volume>49</volume> (<issue>8</issue>), <fpage>815</fpage>&#x2013;<lpage>824</lpage>. <pub-id pub-id-type="doi">10.1136/thx.49.8.815</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rowe</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Quine</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Bhushan</surname>
<given-names>S. S.</given-names>
</name>
<name>
<surname>Stenson</surname>
<given-names>B. J.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>A simplified method for deriving shunt and reduced VA/Q in infants</article-title>. <source>Archives Dis. Child. Fetal neonatal Ed.</source> <volume>95</volume> (<issue>1</issue>), <fpage>F47</fpage>&#x2013;<lpage>F52</lpage>. <pub-id pub-id-type="doi">10.1136/adc.2009.160010</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stoecklin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Rakshasbhuvankar</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Svedenkrans</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Pillow</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Simplified bedside assessment of pulmonary gas exchange in very preterm infants at 36 weeks&#x2019; postmenstrual age</article-title>. <source>Thorax</source> <volume>76</volume>, <fpage>689</fpage>&#x2013;<lpage>695</lpage>. <pub-id pub-id-type="doi">10.1136/thoraxjnl-2020-214659</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Svedenkrans</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Stoecklin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Doherty</surname>
<given-names>D. A.</given-names>
</name>
<name>
<surname>Pillow</surname>
<given-names>J. J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Physiology and predictors of impaired gas exchange in infants with bronchopulmonary dysplasia</article-title>. <source>Am. J. Respir. Crit. Care Med.</source> <volume>200</volume> (<issue>4</issue>), <fpage>471</fpage>&#x2013;<lpage>480</lpage>. <pub-id pub-id-type="doi">10.1164/rccm.201810-2037OC</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>