<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1d3 20150301//EN" "http://jats.nlm.nih.gov/publishing/1.1d3/JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.1d3" 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="nlm-ta">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group>
<journal-title>PLOS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, CA USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338</article-id>
<article-id pub-id-type="publisher-id">PONE-D-20-15330</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v3">
<subject>People and places</subject><subj-group><subject>Geographical locations</subject><subj-group><subject>Asia</subject><subj-group><subject>Japan</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Social sciences</subject><subj-group><subject>Economics</subject><subj-group><subject>Health economics</subject><subj-group><subject>Health insurance</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Medicine and health sciences</subject><subj-group><subject>Health care</subject><subj-group><subject>Health economics</subject><subj-group><subject>Health insurance</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Earth sciences</subject><subj-group><subject>Geography</subject><subj-group><subject>Human geography</subject><subj-group><subject>Urban geography</subject><subj-group><subject>Urban areas</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Social sciences</subject><subj-group><subject>Human geography</subject><subj-group><subject>Urban geography</subject><subj-group><subject>Urban areas</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Earth sciences</subject><subj-group><subject>Geography</subject><subj-group><subject>Geographic areas</subject><subj-group><subject>Urban areas</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>People and places</subject><subj-group><subject>Population groupings</subject><subj-group><subject>Age groups</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Medicine and health sciences</subject><subj-group><subject>Health care</subject><subj-group><subject>Health care policy</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>People and places</subject><subj-group><subject>Population groupings</subject><subj-group><subject>Age groups</subject><subj-group><subject>Adults</subject><subj-group><subject>Elderly</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Social sciences</subject><subj-group><subject>Sociology</subject><subj-group><subject>Population mobility</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3">
<subject>Medicine and health sciences</subject><subj-group><subject>Public and occupational health</subject><subj-group><subject>Global health</subject></subj-group></subj-group></subj-group></article-categories>
<title-group>
<article-title>Effect of population inflow and outflow between rural and urban areas on regional antimicrobial use surveillance</article-title>
<alt-title alt-title-type="running-head">Effect of rural–urban population mobility on regional antimicrobial use surveillance</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Koizumi</surname>
<given-names>Ryuji</given-names>
</name>
<role content-type="https://casrai.org/credit/">Conceptualization</role>
<role content-type="https://casrai.org/credit/">Data curation</role>
<role content-type="https://casrai.org/credit/">Formal analysis</role>
<role content-type="https://casrai.org/credit/">Investigation</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Writing – original draft</role>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes" xlink:type="simple">
<contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0002-7072-2338</contrib-id>
<name name-style="western">
<surname>Kusama</surname>
<given-names>Yoshiki</given-names>
</name>
<role content-type="https://casrai.org/credit/">Conceptualization</role>
<role content-type="https://casrai.org/credit/">Data curation</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Supervision</role>
<role content-type="https://casrai.org/credit/">Writing – review &amp; editing</role>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff002"><sup>2</sup></xref>
<xref ref-type="corresp" rid="cor001">*</xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<contrib-id authenticated="true" contrib-id-type="orcid">https://orcid.org/0000-0001-6951-5676</contrib-id>
<name name-style="western">
<surname>Muraki</surname>
<given-names>Yuichi</given-names>
</name>
<role content-type="https://casrai.org/credit/">Conceptualization</role>
<role content-type="https://casrai.org/credit/">Data curation</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Supervision</role>
<role content-type="https://casrai.org/credit/">Writing – review &amp; editing</role>
<xref ref-type="aff" rid="aff003"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Ishikane</surname>
<given-names>Masahiro</given-names>
</name>
<role content-type="https://casrai.org/credit/">Conceptualization</role>
<role content-type="https://casrai.org/credit/">Data curation</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Supervision</role>
<role content-type="https://casrai.org/credit/">Writing – review &amp; editing</role>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff004"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Yamasaki</surname>
<given-names>Daisuke</given-names>
</name>
<role content-type="https://casrai.org/credit/">Conceptualization</role>
<role content-type="https://casrai.org/credit/">Data curation</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Supervision</role>
<role content-type="https://casrai.org/credit/">Writing – review &amp; editing</role>
<xref ref-type="aff" rid="aff005"><sup>5</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Tanabe</surname>
<given-names>Masaki</given-names>
</name>
<role content-type="https://casrai.org/credit/">Conceptualization</role>
<role content-type="https://casrai.org/credit/">Data curation</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Supervision</role>
<role content-type="https://casrai.org/credit/">Writing – review &amp; editing</role>
<xref ref-type="aff" rid="aff005"><sup>5</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Ohmagari</surname>
<given-names>Norio</given-names>
</name>
<role content-type="https://casrai.org/credit/">Funding acquisition</role>
<role content-type="https://casrai.org/credit/">Methodology</role>
<role content-type="https://casrai.org/credit/">Supervision</role>
<role content-type="https://casrai.org/credit/">Writing – review &amp; editing</role>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff002"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff004"><sup>4</sup></xref>
</contrib>
</contrib-group>
<aff id="aff001"><label>1</label> <addr-line>AMR Clinical Reference Center, National Center for Global Health and Medicine, Tokyo, Japan</addr-line></aff>
<aff id="aff002"><label>2</label> <addr-line>Department of Emerging and Re-emerging Infectious Diseases, Tohoku University School of Medicine, Miyagi, Japan</addr-line></aff>
<aff id="aff003"><label>3</label> <addr-line>Department of Clinical Pharmacoepidemiology, Kyoto Pharmaceutical University, Kyoto, Japan</addr-line></aff>
<aff id="aff004"><label>4</label> <addr-line>Disease Control and Prevention Center, National Center for Global Health and Medicine, Tokyo, Japan</addr-line></aff>
<aff id="aff005"><label>5</label> <addr-line>Department of Infection Control and Prevention, Mie University Hospital, Mie, Japan</addr-line></aff>
<contrib-group>
<contrib contrib-type="editor" xlink:type="simple">
<name name-style="western">
<surname>Suppiah</surname>
<given-names>Vijayaprakash</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"/>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>University of South Australia, AUSTRALIA</addr-line></aff>
<author-notes>
<fn fn-type="conflict" id="coi001">
<p>The authors have declared that no competing interests exist.</p>
</fn>
<corresp id="cor001">* E-mail: <email xlink:type="simple">stone.bagle@gmail.com</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>3</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>16</volume>
<issue>3</issue>
<elocation-id>e0248338</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>5</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>24</day>
<month>2</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-year>2021</copyright-year>
<copyright-holder>Koizumi et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">
<license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="info:doi/10.1371/journal.pone.0248338"/>
<abstract>
<sec id="sec001">
<title>Purpose</title>
<p>Regional-level measures can complement national antimicrobial stewardship programs. In Japan, sub-prefectural regions called secondary medical areas (SMAs) provide general inpatient care within their borders, and regional antimicrobial stewardship measures are frequently implemented at this level. There is therefore a need to conduct antimicrobial use (AMU) surveillance at this level to ascertain antimicrobial consumption. However, AMU estimates are generally standardized to residence-based nighttime populations, which do not account for population mobility across regional borders. We examined the impact of population in/outflow on SMA-level AMU estimates by comparing the differences between standardization using daytime and nighttime populations.</p>
</sec>
<sec id="sec002">
<title>Methods</title>
<p>We obtained AMU information from the National Database of Health Insurance Claims and Specific Health Checkups of Japan. AMU was quantified at the prefectural and SMA levels using the number of defined daily doses (DDDs) divided by (a) 1,000 nighttime population per day or (b) 1,000 daytime population per day. We identified and characterized the discrepancies between the two types of estimates at the prefectural and SMA levels.</p>
</sec>
<sec id="sec003">
<title>Results</title>
<p>The national AMU was 17.21 DDDs per 1,000 population per day. The mean (95% confidence interval) prefectural-level DDDs per 1,000 nighttime and daytime population per day were 17.27 (14.10, 20.44) and 17.41 (14.30, 20.53), respectively. The mean (95% confidence interval) SMA-level DDDs per 1,000 nighttime and daytime population per day were 16.12 (9.84, 22.41) and 16.41 (10.57, 22.26), respectively. The nighttime population-standardized estimates were generally higher than the daytime population-standardized estimates in urban areas, but lower in the adjacent suburbs. Large differences were observed in the main metropolitan hubs in eastern and western Japan.</p>
</sec>
<sec id="sec004">
<title>Conclusion</title>
<p>Regional-level AMU estimates, especially of smaller regions such as SMAs, are susceptible to the use of different populations for standardization. This finding indicates that AMU standardization based on population values is not suitable for AMU estimates in small regions.</p>
</sec>
</abstract>
<funding-group>
<award-group id="award001">
<funding-source>
<institution-wrap>
<institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100003478</institution-id>
<institution>Ministry of Health, Labour and Welfare</institution>
</institution-wrap>
</funding-source>
<award-id>20HA2003</award-id>
<principal-award-recipient>
<name name-style="western">
<surname>Ohmagari</surname>
<given-names>Norio</given-names>
</name>
</principal-award-recipient>
</award-group>
<funding-statement>This work was supported by the Ministry of Health, Labor and Welfare (MHLW) research grant of Japan (20HA2003).</funding-statement>
</funding-group>
<counts>
<fig-count count="4"/>
<table-count count="1"/>
<page-count count="11"/>
</counts>
<custom-meta-group>
<custom-meta id="data-availability">
<meta-name>Data Availability</meta-name>
<meta-value>All relevant data are within the manuscript and its <xref ref-type="sec" rid="sec016">Supporting information</xref> files.</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec id="sec005" sec-type="intro">
<title>Introduction</title>
<p>Antimicrobial resistance is a major global health concern, and there is an urgent need to reduce inappropriate antimicrobial use (AMU) as a countermeasure [<xref ref-type="bibr" rid="pone.0248338.ref001">1</xref>]. In 2016, the Japanese government published the National Action Plan on Antimicrobial Resistance, which highlighted the integral role that regional cooperation plays in complementing national antimicrobial stewardship programs [<xref ref-type="bibr" rid="pone.0248338.ref002">2</xref>]. In order to improve AMU within a specific region, it is first necessary to ascertain the region’s actual drug consumption patterns. However, there is a lack of information on best practice methodologies for regional-level AMU surveillance.</p>
<p>Under Japan’s universal healthcare system, patients are free to seek care at any medical institution across regional borders without restriction. Furthermore, regions can experience daily fluxes in population size due to the inflow and outflow of people across borders for work and schooling. Attempts to characterize AMU at the regional level may therefore be hindered by differences between the locations where antimicrobials are prescribed and where patients reside. Accordingly, there are fundamental difficulties in accurately ascertaining the actual AMU of a region. In regions characterized by construction and air pollution, population migration has been reported to affect urban development and exposure to pollutants [<xref ref-type="bibr" rid="pone.0248338.ref003">3</xref>, <xref ref-type="bibr" rid="pone.0248338.ref004">4</xref>]. We also posit that the impact of population inflow and outflow on AMU estimates would increase as the size of the target regional unit decreases, but the extent of such an effect has yet to be explored.</p>
<p>To date, AMU in Japan has been examined at the prefectural level [<xref ref-type="bibr" rid="pone.0248338.ref005">5</xref>], but there is a lack of information on sub-prefectural regions. Because healthcare policies are frequently implemented at the sub-prefectural level, understanding the trends in AMU at this level can help to identify region-specific problems and support the development of more precise and effective antimicrobial stewardship programs.</p>
<p>National-level AMU is generally indicated using the number of defined daily doses (DDDs) per 1,000 inhabitants per day (DID) based on population statistics. Regional-level DID estimates are dependent on the definition of each region’s population, which in turn is affected by population mobility across borders. However, the effects of different population definitions on DID estimates at the regional level are unknown. To improve our understanding of regional-level AMU surveillance methodologies, this study aimed to elucidate the impact of population inflow and outflow on sub-prefectural DID estimates in Japan.</p>
</sec>
<sec id="sec006" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="sec007">
<title>Japan’s health insurance system and secondary medical areas</title>
<p>In Japan, all residents are required to enroll in health insurance, which enables them to receive healthcare at any medical institution throughout the country. Each enrollee’s insurance plan is dependent on his/her age and occupation. Enrollees pay monthly premiums to their insurers, and also pay a portion of the medical charges (i.e., copayments) at the point of care when receiving insurance-covered treatments and medications. These copayments range from 10% to 30% depending on the enrollee’s insurance plan and income level. The healthcare providers send insurance claims to the applicable insurers through a claims processing agency in order to be reimbursed for the remaining charges.</p>
<p>Japan’s healthcare provision infrastructure involves three increasing levels of geographical divisions—designated primary, secondary, and tertiary medical areas—that serve as units for the implementation of healthcare policies. Primary medical areas comprise the nation’s municipalities, and are equipped to provide basic primary outpatient care. Secondary medical areas (SMAs) are sub-prefectural regions comprising several primary medical areas, and are designed to meet each region’s need for general inpatient care (including emergency care). Tertiary medical areas are mostly represented by prefectures, and provide specialized care that requires advanced technology and equipment. Although primary and tertiary medical areas generally use existing regional borders, SMAs are delineated based on the presence of healthcare facilities that enable them to fulfill their designated functions. SMAs, which are most frequently used as the basic unit for healthcare planning, include designated core hospitals that treat critically ill inpatients, specialized outpatient clinics, and hospitals that provide routine in-hospital treatment. As of December 2020, there are 344 SMAs located throughout Japan (see <xref ref-type="supplementary-material" rid="pone.0248338.s001">S1 Table</xref> for the list of SMAs and their constituent municipalities).</p>
</sec>
<sec id="sec008">
<title>Data source</title>
<p>For this retrospective study, data were obtained from the National Database of Health Insurance Claims and Specific Health Checkups of Japan (NDB). The NDB has collected and maintained insurance claims data provided by the Ministry of Health, Labour and Welfare since April 2009, and these data can be used for research purposes through the submission and approval of an application [<xref ref-type="bibr" rid="pone.0248338.ref006">6</xref>]. Because insurance-covered care accounts for the majority of medical treatments provided in Japan, the NDB represents a near-comprehensive database of all treatments performed throughout the country. However, the database does not include claims data from patients with fully publicly funded healthcare (e.g., patients with intractable diseases, atomic bomb survivors, patients on welfare, patients with tuberculosis, and patients with human immunodeficiency virus infections) and patients who personally pay for all of their medical expenses (e.g., foreign travelers and cosmetic surgery patients). In this study, we calculated the number of prescriptions generated for systemic antimicrobial drugs (both oral and injection) in each SMA in 2015. We also identified the SMAs where each prescribing healthcare facility and dispensing pharmacy were located. The NDB data were accessed in January 2020.</p>
</sec>
<sec id="sec009">
<title>Data processing</title>
<p>Antimicrobial drugs were identified using the J01 classification in the Anatomical Therapeutic Chemical/Defined Daily Doses system established by the World Health Organization’s Collaborating Centre for Drug Statistics Methodology [<xref ref-type="bibr" rid="pone.0248338.ref007">7</xref>]. The populations used for analyses were the national population stratified by municipality (i.e., cities, towns, villages, and wards) and the daytime population published by the Statistics Bureau of the Ministry of Internal Affairs and Communications [<xref ref-type="bibr" rid="pone.0248338.ref008">8</xref>]. The national population estimate is the estimated size of the population on every October 1st that reflects the natural population growth, social dynamics, and nationwide migration of Japanese nationals. These parameters are reported for each year (from October 1st of the previous year to September 30th of the index year) based on population data obtained from a national census of all households (conducted every five years) and the annual population inflow and outflow estimates. These population estimates have conventionally been used to calculate the national DID for AMU surveillance. As these estimates are based on residences, they represent the nighttime population. The daytime and nighttime populations at the national, prefectural, and SMA levels are presented in <xref ref-type="supplementary-material" rid="pone.0248338.s002">S2 Table</xref>.</p>
<p>The daytime population of a region accounts for the number of people at work or school, and was calculated based on the national census using the formula shown below.</p>
<disp-formula id="pone.0248338.e001">
<alternatives>
<graphic id="pone.0248338.e001g" mimetype="image" position="anchor" xlink:href="info:doi/10.1371/journal.pone.0248338.e001" xlink:type="simple"/>
<mml:math display="block" id="M1">
<mml:mrow><mml:mi>D</mml:mi><mml:mi>a</mml:mi><mml:mi>y</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mspace width="4pt"/><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>p</mml:mi><mml:mi>u</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi><mml:mi>i</mml:mi><mml:mi>g</mml:mi><mml:mi>h</mml:mi><mml:mi>t</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mspace width="4pt"/><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>p</mml:mi><mml:mi>u</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mo>−</mml:mo><mml:mo>(</mml:mo><mml:mrow><mml:mi>O</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi><mml:mi>f</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mspace width="4pt"/><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>p</mml:mi><mml:mi>u</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mi>I</mml:mi><mml:mi>n</mml:mi><mml:mi>f</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mspace width="4pt"/><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>p</mml:mi><mml:mi>u</mml:mi><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow>
</mml:math>
</alternatives>
</disp-formula>
<p>The nighttime and daytime populations of each SMA were calculated by totaling the respective populations of its municipalities.</p>
</sec>
<sec id="sec010">
<title>Analysis</title>
<p>The study period was 2015, which was the most recent year with available statistics on the daytime population. For this study, AMU was quantified using DID estimates. National-level, prefectural-level, and SMA-level DDDs were shown in <xref ref-type="supplementary-material" rid="pone.0248338.s003">S3 Table</xref>. We calculated and compared the national-level, prefectural-level, and SMA-level DID values that were standardized to either the nighttime population or daytime population. To calculate the AMU at the various levels, the DID values of their constituent regions were totaled. Next, we evaluated the distributions of daytime and nighttime population–standardized DID values at the prefectural and SMA levels using the Kolmogorov–Smirnov test. The mean DID values were compared with the national AMU using the one-sample <italic>t</italic>-test. The population-standardized DID at the prefectural and SMA levels were used to generate violin plots, which were examined to identify regions with notable discrepancies between the two types of populations. Correlations between nighttime and daytime population–standardized DID values were examined using Pearson’s correlation coefficients.</p>
<p>We then calculated the difference between the nighttime population–standardized DID values and the daytime population–standardized DID values for each prefecture and SMA. Choropleth maps were created based on these differences, and the distributions of regions with substantial differences were examined. The Tokyo/Ku-chuoubu SMA was excluded from the choropleth map as it was an extreme outlier. Furthermore, we identified the SMAs with the largest positive and negative differences in DID values, as well as the SMAs with the smallest absolute differences.</p>
<p>Finally, we used the population number for each age category (children: &lt;15 years, working-age adults: 15–64 years, and older adults: ≥65 years) for each SMA, and analyzed the correlation between the absolute difference in DID values and each age category. For this analysis, the Tokyo/Ku-chuoubu SMA was excluded as it was an extreme outlier.</p>
</sec>
<sec id="sec011">
<title>Data management</title>
<p>The mapping of prefectures was performed using Tableau version 2019.1.0 (Tableau Software, Washington, USA). For the visualization of the SMAs, geocoding was performed using the geographical information in Tableau based on the National Land Numerical Information published by the National Spatial Planning and Regional Policy Bureau of the Ministry of Land, Infrastructure, Transport and Tourism. Finally, correlations between the absolute difference in DID values in each SMA and the population age categories were examined using Pearson’s correlation coefficients. Statistical analyses were performed using R ver 4.0.0 (R Core Team, Vienna, Austria), and <italic>P</italic> values below 0.05 were considered statistically significant.</p>
</sec>
<sec id="sec012">
<title>Ethics</title>
<p>The study did not involve any interventions in human subjects, and the NDB data were anonymized before being received by the authors. This study was approved by the institutional review board of the National Center for Global Health and Medicine (Approval Number: NCGM-G-002505-00).</p>
</sec>
</sec>
<sec id="sec013" sec-type="results">
<title>Results</title>
<p><xref ref-type="fig" rid="pone.0248338.g001">Fig 1</xref> shows the national-level, prefectural-level, and SMA-level DID values that were standardized to either the nighttime population or daytime population. When standardized with the nighttime population, the mean DID was 17.27 (95% confidence interval [14.10, 20.44]) at the prefectural level and 16.12 (95% confidence interval [9.84, 22.41]) at the SMA level. When standardized with the daytime population, the mean DID was 17.41 (95% confidence interval [14.30, 20.53]) at the prefectural level and 16.41 (95% confidence interval [10.57, 22.26]) at the SMA level. Both the prefectural-level and SMA-level DID values were normally distributed regardless of the population used for standardization. The daytime and nighttime population–standardized mean DID values at the prefectural level were not significantly different from the national-level DID values (nighttime: <italic>P</italic> = .385; daytime: <italic>P</italic> = .811); however, the corresponding mean DID values at the SMA level were significantly different from the national-level DID values (nighttime: <italic>P</italic> &lt; .001; daytime: <italic>P</italic> &lt; .001). As shown in <xref ref-type="fig" rid="pone.0248338.g001">Fig 1</xref>, daytime population–standardized DID at the SMA level had a narrower dispersion and a median value (center of the widest section in the violin plot) that was closer to the national-level DID than the nighttime population–standardized DID. In contrast, the daytime population and nighttime population–standardized DID at the prefectural level exhibited similar shapes in the violin plot.</p>
<fig id="pone.0248338.g001" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0248338.g001</object-id>
<label>Fig 1</label>
<caption>
<title>Violin plots of prefectural-level and SMA-level DDDs per 1,000 nighttime population per day and DDDs per 1,000 daytime population per day.</title>
<p>The blue broken line represents the national-level antimicrobial use. DDD, defined daily dose; SMA, secondary medical area.</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.g001" xlink:type="simple"/>
</fig>
<p><xref ref-type="fig" rid="pone.0248338.g002">Fig 2</xref> shows scatter plots of nighttime population–standardized DID values against daytime population–standardized DID values. The correlation coefficient between nighttime and daytime population–standardized DID values was higher at the prefectural level (0.90; <italic>P</italic> &lt; .001) than at the SMA level (0.80; <italic>P</italic> &lt; .001). At the prefectural level, the nighttime population–standardized DID values were higher than the daytime population–standardized DID values in central urban areas such as Tokyo and Osaka, but lower in surrounding prefectures such as Saitama, Kanagawa, Chiba, Gifu, and Mie. At the SMA level, the nighttime population–standardized DID value for Tokyo/Ku-chuoubu—which has a high concentration of companies and schools—was extremely high. In addition, the nighttime population–standardized DID values were higher in the urban SMAs of Tokyo/Ku-seibu and Osaka/Osaka-shi, but lower in surrounding SMAs such as Aichi/Ama, Wakayama/Naga, Fukuoka/Munakata, Tokyo/Kitatama-hokubu, and Kyoto/Yamashiro-minami.</p>
<fig id="pone.0248338.g002" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0248338.g002</object-id>
<label>Fig 2</label>
<caption>
<title>Scatter plots of DDDs per 1,000 daytime population per day vs. DDDs per 1,000 nighttime population per day.</title>
<p>(A) Prefectural level and (B) Secondary medical area level. The dotted lines represent y = x. DDD, defined daily dose.</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.g002" xlink:type="simple"/>
</fig>
<p><xref ref-type="fig" rid="pone.0248338.g003">Fig 3</xref> shows choropleth maps of the differences between the nighttime population–standardized and daytime population–standardized DID values at the prefectural and SMA levels. Among the prefectures, a positive difference (indicating that the nighttime population–standardized DID values were higher than the daytime population–standardized DID values) was observed for prefectures that contained more urban areas, such as Tokyo, Osaka, Aichi, and Kyoto. In contrast, a negative difference (indicating that the nighttime population–standardized DID values were lower than the daytime population–standardized DID values) was observed for prefectures adjacent to the urban areas. Prefectures showing substantial differences in DID values were limited to the eastern Kanto (which includes the Greater Tokyo Area) and western Kansai (which includes Osaka, Hyogo, Kyoto, and surrounding prefectures) metropolitan regions, with no notable differences observed in the other regions. Similarly, a positive difference in DID values was generally observed in urban SMAs, whereas a negative difference was observed in the surrounding SMAs.</p>
<fig id="pone.0248338.g003" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0248338.g003</object-id>
<label>Fig 3</label>
<caption>
<title>Choropleth maps of Japan showing the differences between DDDs per 1,000 nighttime population per day and DDDs per 1,000 daytime population per day.</title>
<p>(A) Prefectural level and (B) Secondary medical area level. The numbers represent the differences between DDDs per 1,000 nighttime population per day and DDDs per 1,000 daytime population per day. The red and blue colors represent positive and negative differences, respectively. Data from Tokyo/Ku-chuoubu are not shown in Fig 3B because it was an extreme outlier (43.08). DDD, defined daily dose.</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.g003" xlink:type="simple"/>
</fig>
<p>The SMAs with the largest positive, largest negative, and smallest absolute differences between the nighttime population–standardized and daytime population–standardized DID values are presented in <xref ref-type="table" rid="pone.0248338.t001">Table 1</xref>. Highly urbanized SMAs in Tokyo, Osaka, Aichi, and Fukuoka prefectures showed large positive differences. In contrast, SMAs in Kyoto, Tokyo, Fukuoka, Kanagawa, and Saitama prefectures (which are adjacent to urbanized areas) showed large negative differences. SMAs with the smallest absolute values were mostly found in rural prefectures such as Hokkaido, Shimane, and Kagoshima.</p>
<table-wrap id="pone.0248338.t001" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0248338.t001</object-id>
<label>Table 1</label> <caption><title>Difference between nighttime population–standardized and daytime population–standardized defined daily doses per 1,000 population per day.</title></caption>
<alternatives>
<graphic id="pone.0248338.t001g" mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.t001" xlink:type="simple"/>
<table>
<colgroup>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
</colgroup>
<thead>
<tr>
<th align="left">Rank</th>
<th align="left">Prefecture/Secondary medical area</th>
<th align="left">Difference</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" colspan="3">Difference in AMU between nighttime population–standardized and daytime population–standardized (ranked in descending order)</td>
</tr>
<tr>
<td align="left">1</td>
<td align="left">Tokyo/Ku-chuoubu</td>
<td align="char" char=".">30.94</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">Osaka/Osaka-shi</td>
<td align="char" char=".">4.94</td>
</tr>
<tr>
<td align="left">3</td>
<td align="left">Tokyo/Ku-seibu</td>
<td align="char" char=".">4.50</td>
</tr>
<tr>
<td align="left">4</td>
<td align="left">Tokyo/Ku-seinanbu</td>
<td align="char" char=".">3.15</td>
</tr>
<tr>
<td align="left">5</td>
<td align="left">Fukushima/Soso</td>
<td align="char" char=".">2.30</td>
</tr>
<tr>
<td align="left">6</td>
<td align="left">Aichi/Nagoya</td>
<td align="char" char=".">2.24</td>
</tr>
<tr>
<td align="left">7</td>
<td align="left">Tokyo/Ku-nanbu</td>
<td align="char" char=".">1.99</td>
</tr>
<tr>
<td align="left">8</td>
<td align="left">Fukuoka/Itoshima</td>
<td align="char" char=".">1.49</td>
</tr>
<tr>
<td align="left">9</td>
<td align="left">Fukuoka/Noogata, Kurate</td>
<td align="char" char=".">1.34</td>
</tr>
<tr>
<td align="left">10</td>
<td align="left">Aichi/Nishimikawa-hokubu</td>
<td align="char" char=".">1.25</td>
</tr>
<tr>
<td align="left" colspan="3">Difference in AMU between nighttime population–standardized and daytime population–standardized (ranked in ascending order)</td>
</tr>
<tr>
<td align="left">1</td>
<td align="left">Kyoto/Yamashiro-minami</td>
<td align="char" char=".">-4.42</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">Tokyo/Kitatama-hokubu</td>
<td align="char" char=".">-3.78</td>
</tr>
<tr>
<td align="left">3</td>
<td align="left">Fukuoka/Munakata</td>
<td align="char" char=".">-3.70</td>
</tr>
<tr>
<td align="left">4</td>
<td align="left">Kanagawa/Kawasaki-hokubu</td>
<td align="char" char=".">-3.69</td>
</tr>
<tr>
<td align="left">5</td>
<td align="left">Wakayama/Naga</td>
<td align="char" char=".">-3.28</td>
</tr>
<tr>
<td align="left">6</td>
<td align="left">Saitama/Nanbu</td>
<td align="char" char=".">-2.91</td>
</tr>
<tr>
<td align="left">7</td>
<td align="left">Aichi/Ama</td>
<td align="char" char=".">-2.86</td>
</tr>
<tr>
<td align="left">8</td>
<td align="left">Saitama/Keno</td>
<td align="char" char=".">-2.85</td>
</tr>
<tr>
<td align="left">9</td>
<td align="left">Nara/Seiwa</td>
<td align="char" char=".">-2.82</td>
</tr>
<tr>
<td align="left">10</td>
<td align="left">Chiba/Toukatsu-hokubu</td>
<td align="char" char=".">-2.82</td>
</tr>
<tr>
<td align="left" colspan="3">Absolute value of the AMU difference between nighttime population–-standardized and daytime population–-standardized (ranked in ascending order)</td>
</tr>
<tr>
<td align="left">1</td>
<td align="left">Shimane/Masuda</td>
<td align="char" char=".">0.00</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">Hokkaido/Tokachi</td>
<td align="char" char=".">0.00</td>
</tr>
<tr>
<td align="left">3</td>
<td align="left">Yamagata/Murayama</td>
<td align="char" char=".">0.00</td>
</tr>
<tr>
<td align="left">4</td>
<td align="left">Hokkaido/Kushiro</td>
<td align="char" char=".">0.00</td>
</tr>
<tr>
<td align="left">5</td>
<td align="left">Gifu/Hida</td>
<td align="char" char=".">0.01</td>
</tr>
<tr>
<td align="left">6</td>
<td align="left">Fukushima/Kennaka</td>
<td align="char" char=".">0.01</td>
</tr>
<tr>
<td align="left">7</td>
<td align="left">Kouchi/Chuou</td>
<td align="char" char=".">0.01</td>
</tr>
<tr>
<td align="left">8</td>
<td align="left">Kagoshima/Nansatsu</td>
<td align="char" char=".">0.01</td>
</tr>
<tr>
<td align="left">9</td>
<td align="left">Kagoshima/Amami</td>
<td align="char" char=".">0.01</td>
</tr>
<tr>
<td align="left">10</td>
<td align="left">Ehime/Matsuyama</td>
<td align="char" char=".">0.01</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>The nighttime population–standardized and daytime population–standardized DID values according to age group were calculated. <xref ref-type="fig" rid="pone.0248338.g004">Fig 4</xref> shows the correlations between the absolute difference in DID values and the proportion of each age group in the population. The correlation coefficients were 0.14 (<italic>P</italic> = 0.0082) for children, 0.49 (<italic>P</italic> &lt; 0.001) for working-age adults, and -0.44 (<italic>P</italic> &lt; 0.001) for older adults.</p>
<fig id="pone.0248338.g004" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0248338.g004</object-id>
<label>Fig 4</label>
<caption>
<title>Scatter plots of the proportion of age groups vs. absolute differences between DDDs per 1,000 nighttime population per day and DDDs per 1,000 daytime population per day.</title>
<p>(A) &lt;15 years, (B) 15–64 years, and (C) ≥65 years. Each dot represents a secondary medical area. The correlation coefficients were 0.14 (<italic>P</italic> = 0.0082) for children (&lt;15 years), 0.49 (<italic>P</italic> &lt;0.001) for working-age adults (15–64 years), and -0.44 (<italic>P</italic> &lt;0.001) for older adults (≥65 years). Data from Tokyo/Ku-chuoubu are not shown in the plots because this region was an extreme outlier (43.08).</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.g004" xlink:type="simple"/>
</fig>
</sec>
<sec id="sec014" sec-type="conclusions">
<title>Discussion</title>
<p>In this retrospective nationwide study, we comparatively examined the effects of using the nighttime population and daytime population to adjust regional-level AMU in Japan. Previous studies have examined the effects of using different denominator values when calculating AMU in hospitals [<xref ref-type="bibr" rid="pone.0248338.ref009">9</xref>–<xref ref-type="bibr" rid="pone.0248338.ref012">12</xref>]. However, these effects on regional AMU surveillance have not been explored. Even the World Health Organization’s AMU surveillance methodology does not address the appropriate methods for analyzing sub-national regions [<xref ref-type="bibr" rid="pone.0248338.ref013">13</xref>]. In our analysis, SMA-level AMU (standardized using either the daytime or nighttime population) was found to be significantly different from the national-level AMU. This suggests that the population-standardized DID values of smaller regions are susceptible to the effects of population inflow and outflow, which can lead to erroneous results.</p>
<p>As more medical examinations and prescriptions are received during the day than at night, the calculation of regional AMU indices should account for the effects of population inflow and outflow. At the SMA level, daytime population–standardization produced fewer outliers, narrower 95% confidence intervals, and mean DID values that were closer to the national DID than nighttime population–standardization. These findings showed that when analyzing smaller regional units such as SMAs, the use of different populations for standardization has a considerable effect on DID estimates in central urban areas and their surrounding regions. Our insights indicate that AMU standardization based on population values is not suitable for AMU estimates in small regions.</p>
<p>When comparing the daytime population–standardized DID values with the nighttime population–standardized values at the prefectural level, the former tended to be higher in bedroom communities such as Gifu, Nara, Kanagawa, Chiba, and Saitama, but lower in the central urban areas of Tokyo and Osaka. A similar trend was observed at the SMA level. These observations may be explained by the higher concentration of people in the city centers during the day for work or schooling (i.e., population outflow from bedroom communities during the day). In particular, there was a large difference between the DID values standardized for nighttime population (43.08) and daytime population (12.15) in the Tokyo/Ku-chuoubu SMA, which experiences high population inflow during the day. However, substantial differences were mainly observed in the Kanto and Kansai regions that are centered around large metropolitan areas, with many other regions unaffected. An exception was the Fukushima/Soso region, which underwent evacuations due to the nuclear power plant disaster in 2011. The SMAs with small absolute differences tended to be located in rural regions and remote islands with higher levels of medical self-sufficiency (e.g., Hokkaido/Tokachi, Yamagata/Murayama, Kagoshima/Amami, and Hokkaido/Kushiro), which would therefore have low population outflow.</p>
<p>Based on a hypothesis that population age structure may affect mobility across regional borders, we examined the correlations between the proportion of each age group in the population and the absolute difference in DID values. There was a significant and positive correlation between the proportion of working-age adults (who would be the most mobile among the age groups) and the absolute difference in DID values. Regions with a high proportion of working-age adults would experience higher population mobility, resulting in a larger difference in DID values between the different populations. In contrast, a significant and negative correlation was observed for older adults, who would have the least mobility among the age groups. Therefore, regions with many older adults would experience lower population mobility, resulting in a smaller difference in DID values between the different populations.</p>
<p>This study has several limitations. First, the NDB does not include information on the treatment of patients who pay their own expenses and those for whom the municipality shares the cost. However, almost all necessary medical services in Japan are covered by health insurance. Therefore, the NDB covers the majority of healthcare provided throughout the country. Next, daytime population statistics are only published once every five years in Japan. As these statistics are based on weekday estimates, weekends and holidays (accounting for approximately 30% of the year) are overlooked. Finally, the SMAs may differ from year to year due to the merging of municipalities. Although not examined in this study, AMU surveillance should consider such regional changes over time. Despite these limitations, our findings demonstrated the effects of population inflow and outflow on the population-standardized AMU.</p>
<p>Although regional AMU estimates would help to inform AMU-related policymaking, we recommend avoiding AMU evaluation in small regions standardized by the population number. If we wish to monitor AMU evaluations in small regions, only the temporal change in a region should be considered, for example.</p>
</sec>
<sec id="sec015" sec-type="conclusions">
<title>Conclusion</title>
<p>AMU surveillance has conventionally used the nighttime population for standardization. However, regional-level AMU estimates, especially of smaller regions such as SMAs, are more susceptible to the influence of whatever population is used for standardization, which can lead to erroneous estimates. Therefore, new approaches are required to monitor AMU evaluations in small regions, for example, observing only the temporal changes in a region.</p>
</sec>
<sec id="sec016" sec-type="supplementary-material">
<title>Supporting information</title>
<supplementary-material id="pone.0248338.s001" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.s001" xlink:type="simple">
<label>S1 Table</label>
<caption>
<title>Prefectures, secondary medical areas, and municipalities in Japan.</title>
<p>(DOCX)</p>
</caption>
</supplementary-material>
<supplementary-material id="pone.0248338.s002" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.s002" xlink:type="simple">
<label>S2 Table</label>
<caption>
<title>Daytime and nighttime populations at the national level, prefectural level, and secondary medical area level in Japan.</title>
<p>(DOCX)</p>
</caption>
</supplementary-material>
<supplementary-material id="pone.0248338.s003" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.s003" xlink:type="simple">
<label>S3 Table</label>
<caption>
<title>Defined daily doses of antimicrobials at the national level, prefectural level, and secondary medical area level in Japan.</title>
<p>(DOCX)</p>
</caption>
</supplementary-material>
</sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="pone.0248338.ref001"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">World Health Organization. Global action plan on antimicrobial resistance. 2015 [Cited 2020 Mar 12] <ext-link ext-link-type="uri" xlink:href="https://apps.who.int/iris/bitstream/handle/10665/193736/9789241509763_eng.pdf?sequence=1" xlink:type="simple">https://apps.who.int/iris/bitstream/handle/10665/193736/9789241509763_eng.pdf?sequence=1</ext-link></mixed-citation></ref>
<ref id="pone.0248338.ref002"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">The Government of Japan. National action plan on antimicrobial resistance (AMR) 2016–2020 [Internet]. 2016. <ext-link ext-link-type="uri" xlink:href="https://www.mhlw.go.jp/file/06-Seisakujouhou-10900000-Kenkoukyoku/0000138942.pdf" xlink:type="simple">https://www.mhlw.go.jp/file/06-Seisakujouhou-10900000-Kenkoukyoku/0000138942.pdf</ext-link></mixed-citation></ref>
<ref id="pone.0248338.ref003"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Qi</surname> <given-names>W</given-names></name>, <name name-style="western"><surname>Liu</surname> <given-names>S</given-names></name>, <name name-style="western"><surname>Gao</surname> <given-names>X</given-names></name>, <name name-style="western"><surname>Zhao</surname> <given-names>M</given-names></name>. <article-title>Modeling the spatial distribution of urban population during the daytime and at night based on land use: A case study in Beijing, China</article-title>. <source>J Geogr Sci</source>. <year>2015</year>;<volume>25</volume>: <fpage>756</fpage>–<lpage>68</lpage>.</mixed-citation></ref>
<ref id="pone.0248338.ref004"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Park</surname> <given-names>YM</given-names></name>, <name name-style="western"><surname>Kwan</surname> <given-names>MP</given-names></name>. <article-title>Understanding racial disparities in exposure to traffic-related air pollution: Considering the spatiotemporal dynamics of population distribution</article-title>. <source>Int J Environ Res Public Health</source>. <year>2020</year>;<volume>17</volume>(<issue>3</issue>): <fpage>908</fpage>. Published 2020 Feb 1. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3390/ijerph17030908" xlink:type="simple">10.3390/ijerph17030908</ext-link></comment> <object-id pub-id-type="pmid">32024171</object-id></mixed-citation></ref>
<ref id="pone.0248338.ref005"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kusama</surname> <given-names>Y</given-names></name>, <name name-style="western"><surname>Ishikane</surname> <given-names>M</given-names></name>, <name name-style="western"><surname>Tanaka</surname> <given-names>C</given-names></name>, <name name-style="western"><surname>Kimura</surname> <given-names>Y</given-names></name>, <name name-style="western"><surname>Yumura</surname> <given-names>E</given-names></name>, <name name-style="western"><surname>Hayakawa</surname> <given-names>K</given-names></name>, <etal>et al</etal>. <article-title>Regional variation of antimicrobial use in Japan from 2013–2016, as estimated by the sales data</article-title>. <source>Jpn J Infect Dis</source>. <year>2019</year>;<volume>72</volume>: <fpage>326</fpage>–<lpage>9</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.7883/yoken.JJID.2018.417" xlink:type="simple">10.7883/yoken.JJID.2018.417</ext-link></comment> <object-id pub-id-type="pmid">31061356</object-id></mixed-citation></ref>
<ref id="pone.0248338.ref006"><label>6</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Matsuda</surname> <given-names>S</given-names></name>, <name name-style="western"><surname>Fujimori</surname> <given-names>K</given-names></name>. <article-title>The claim database in Japan</article-title>. <source>Asian Pac J Dis Manag</source>. <year>2012</year>;<volume>6</volume>: <fpage>55</fpage>–<lpage>9</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.7223/apjdm.6.55" xlink:type="simple">10.7223/apjdm.6.55</ext-link></comment></mixed-citation></ref>
<ref id="pone.0248338.ref007"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">WHO Collaborating Centre for Drug Statistics Methodology. ATC/DDD Index 2019. <ext-link ext-link-type="uri" xlink:href="https://www.whocc.no/atc_ddd_index/" xlink:type="simple">https://www.whocc.no/atc_ddd_index/</ext-link></mixed-citation></ref>
<ref id="pone.0248338.ref008"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Ministry of Internal Affairs and Communications. Portal Site of Official Statistics of Japan -e-Stat. <ext-link ext-link-type="uri" xlink:href="https://www.e-stat.go.jp/en" xlink:type="simple">https://www.e-stat.go.jp/en</ext-link></mixed-citation></ref>
<ref id="pone.0248338.ref009"><label>9</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Avedissian</surname> <given-names>SN</given-names></name>, <name name-style="western"><surname>Scheetz</surname> <given-names>MH</given-names></name>, <name name-style="western"><surname>Zembower</surname> <given-names>TR</given-names></name>, <name name-style="western"><surname>Silkaitis</surname> <given-names>C</given-names></name>, <name name-style="western"><surname>Maxwell</surname> <given-names>R</given-names></name>, <name name-style="western"><surname>Jenkins</surname> <given-names>C</given-names></name>, <etal>et al</etal>. <article-title>Measuring the impact of varying denominator definitions on standardized antibiotic consumption rates: Implications for antimicrobial stewardship programmes</article-title>. <source>J Antimicrob Chemother</source>. <year>2018</year>;<volume>73</volume>: <fpage>2876</fpage>–<lpage>82</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1093/jac/dky275" xlink:type="simple">10.1093/jac/dky275</ext-link></comment> <object-id pub-id-type="pmid">30085084</object-id></mixed-citation></ref>
<ref id="pone.0248338.ref010"><label>10</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kuster</surname> <given-names>SP</given-names></name>, <name name-style="western"><surname>Ruef</surname> <given-names>C</given-names></name>, <name name-style="western"><surname>Ledergerber</surname> <given-names>B</given-names></name>, <name name-style="western"><surname>Hintermann</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Deplazes</surname> <given-names>C</given-names></name>, <name name-style="western"><surname>Neuber</surname> <given-names>L</given-names></name>, <etal>et al</etal>. <article-title>Quantitative antibiotic use in hospitals: Comparison of measurements, literature review, and recommendations for a standard of reporting</article-title>. <source>Infection</source>. <year>2008</year>;<volume>36</volume>: <fpage>549</fpage>–<lpage>59</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1007/s15010-008-7462-z" xlink:type="simple">10.1007/s15010-008-7462-z</ext-link></comment> <object-id pub-id-type="pmid">19011740</object-id></mixed-citation></ref>
<ref id="pone.0248338.ref011"><label>11</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Stanić Benić</surname> <given-names>M</given-names></name>, <name name-style="western"><surname>Milanič</surname> <given-names>R</given-names></name>, <name name-style="western"><surname>Monnier</surname> <given-names>AA</given-names></name>, <name name-style="western"><surname>Gyssens</surname> <given-names>IC</given-names></name>, <name name-style="western"><surname>Adriaenssens</surname> <given-names>N</given-names></name>, <name name-style="western"><surname>Versporten</surname> <given-names>A</given-names></name>, <etal>et al</etal>. <article-title>Metrics for quantifying antibiotic use in the hospital setting: Results from a systematic review and international multidisciplinary consensus procedure</article-title>. <source>J Antimicrob Chemother</source>. <year>2018</year>;<volume>73</volume>: <fpage>vi50</fpage>–<lpage>8</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1093/jac/dky118" xlink:type="simple">10.1093/jac/dky118</ext-link></comment> <object-id pub-id-type="pmid">29878222</object-id></mixed-citation></ref>
<ref id="pone.0248338.ref012"><label>12</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Yu</surname> <given-names>KC</given-names></name>, <name name-style="western"><surname>Moisan</surname> <given-names>E</given-names></name>, <name name-style="western"><surname>Tartof</surname> <given-names>SY</given-names></name>, <name name-style="western"><surname>Nguyen</surname> <given-names>HM</given-names></name>, <name name-style="western"><surname>Rieg</surname> <given-names>G</given-names></name>, <name name-style="western"><surname>Ramaprasad</surname> <given-names>C</given-names></name>, <etal>et al</etal>. <article-title>Benchmarking inpatient antimicrobial use: A comparison of risk-adjusted observed-to-expected ratios</article-title>. <source>Clin Infect Dis</source>. <year>2018</year>;<volume>67</volume>: <fpage>1677</fpage>–<lpage>85</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1093/cid/ciy354" xlink:type="simple">10.1093/cid/ciy354</ext-link></comment> <object-id pub-id-type="pmid">29688279</object-id></mixed-citation></ref>
<ref id="pone.0248338.ref013"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">World Health Organization. WHO methodology for a global programme on surveillance of antimicrobial consumption. <ext-link ext-link-type="uri" xlink:href="https://www.who.int/medicines/areas/rational_use/WHO_AMCsurveillance_1.0.pdf?ua=1" xlink:type="simple">https://www.who.int/medicines/areas/rational_use/WHO_AMCsurveillance_1.0.pdf?ua=1</ext-link></mixed-citation></ref>
</ref-list>
</back>
<sub-article article-type="aggregated-review-documents" id="pone.0248338.r001" specific-use="decision-letter">
<front-stub>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338.r001</article-id>
<title-group>
<article-title>Decision Letter 0</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name name-style="western">
<surname>Hepp</surname>
<given-names>Susan</given-names>
</name>
<role>Staff Editor</role>
</contrib>
</contrib-group>
<permissions>
<copyright-year>2021</copyright-year>
<copyright-holder>Susan Hepp</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<related-object document-id="10.1371/journal.pone.0248338" document-id-type="doi" document-type="article" id="rel-obj001" link-type="peer-reviewed-article"/>
<custom-meta-group>
<custom-meta>
<meta-name>Submission Version</meta-name>
<meta-value>0</meta-value>
</custom-meta>
</custom-meta-group>
</front-stub>
<body>
<p>
<named-content content-type="letter-date">17 Nov 2020</named-content>
</p>
<p>PONE-D-20-15330</p>
<p>Effect of population inflow and outflow between rural and urban areas on regional antimicrobial use surveillance</p>
<p>PLOS ONE</p>
<p>Dear Dr. Kusama,</p>
<p>Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.</p>
<p>Specifically, two experts in the field reviewed your manuscript. Their major concerns relate to statistical analyses and elaboration of the discussion.</p>
<p>Please submit your revised manuscript by Dec 31 2020 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at <email xlink:type="simple">plosone@plos.org</email>. When you're ready to submit your revision, log on to <ext-link ext-link-type="uri" xlink:href="https://www.editorialmanager.com/pone/" xlink:type="simple">https://www.editorialmanager.com/pone/</ext-link> and select the 'Submissions Needing Revision' folder to locate your manuscript file.</p>
<p>Please include the following items when submitting your revised manuscript:</p>
<p><list list-type="bullet"><list-item><p>A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.</p></list-item><list-item><p>A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.</p></list-item><list-item><p>An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.</p></list-item></list></p>
<p>If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.</p>
<p>If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: <ext-link ext-link-type="uri" xlink:href="http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols" xlink:type="simple">http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols</ext-link></p>
<p>We look forward to receiving your revised manuscript.</p>
<p>Kind regards,</p>
<p>Susan Hepp</p>
<p>Academic Editor</p>
<p>PLOS ONE</p>
<p>Journal Requirements:</p>
<p>When submitting your revision, we need you to address these additional requirements.</p>
<p>1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at</p>
<p><ext-link ext-link-type="uri" xlink:href="https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf" xlink:type="simple">https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf</ext-link> and</p>
<p><ext-link ext-link-type="uri" xlink:href="https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf" xlink:type="simple">https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf</ext-link></p>
<p>2. In your ethics statement in the Methods section and in the online submission form, please clarify whether all data were fully anonymized before you accessed them.</p>
<p>3. Please state the date(s) when the data was obtained/accessed from the NDB.</p>
<p>4. Thank you for stating in your manuscript text</p>
<p>"This study was approved by the institutional review board of the National Center for Global Health and Medicine (Approval Number: NCGM-G-002505-00)."</p>
<p>Please also include this in your ethics statement in the online submission form.</p>
<p>5. Please include a table of relevant demographic details summarizing the studied population.</p>
<p>6. Please upload a copy of Figure 4, to which you refer in your text on page 16. If the figure is no longer to be included as part of the submission please remove all reference to it within the text.</p>
<p>7. Please include captions for your Supporting Information files at the end of your manuscript, and update any in-text citations to match accordingly. Please see our Supporting Information guidelines for more information: <ext-link ext-link-type="uri" xlink:href="about:blank" xlink:type="simple">http://journals.plos.org/plosone/s/supporting-information</ext-link></p>
<p>[Note: HTML markup is below. Please do not edit.]</p>
<p>Reviewers' comments:</p>
<p>Reviewer's Responses to Questions</p>
<p><!-- <font color="black"> --><bold>Comments to the Author</bold></p>
<p>1. Is the manuscript technically sound, and do the data support the conclusions?</p>
<p>The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. <!-- </font> --></p>
<p>Reviewer #1: Yes</p>
<p>Reviewer #2: Partly</p>
<p>**********</p>
<p><!-- <font color="black"> -->2. Has the statistical analysis been performed appropriately and rigorously? <!-- </font> --></p>
<p>Reviewer #1: Yes</p>
<p>Reviewer #2: No</p>
<p>**********</p>
<p><!-- <font color="black"> -->3. Have the authors made all data underlying the findings in their manuscript fully available?</p>
<p>The <ext-link ext-link-type="uri" xlink:href="http://www.plosone.org/static/policies.action#sharing" xlink:type="simple">PLOS Data policy</ext-link> requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.<!-- </font> --></p>
<p>Reviewer #1: Yes</p>
<p>Reviewer #2: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->4. Is the manuscript presented in an intelligible fashion and written in standard English?</p>
<p>PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.<!-- </font> --></p>
<p>Reviewer #1: Yes</p>
<p>Reviewer #2: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->5. Review Comments to the Author</p>
<p>Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)<!-- </font> --></p>
<p>Reviewer #1: Thank you for the opportunity to review this educational manuscript by Koizumi et al. They describe the effects of regional population inflow and outflow on estimates of antimicrobial use in Japan, and the implications this may have for future policies regarding antimicrobial stewardship. I have a few comments and questions that I hope will serve to improve the manuscript further:</p>
<p>1) The discussion is concise, but would benefit from a further elaboration of how this recommended shift from nighttime to daytime population usage will change antimicrobial use policy and management. What implications does this have for future regulation? What steps will be taken next? A more in-depth discussion would be helpful.</p>
<p>2) I did not receive Figure 4 – the description is provided, but not the figure itself. This should be provided and reviewed.</p>
<p>3) Figure 3 displays at a quite small size, making the color of the map and differences between region\\SMA difficult to see. Would enlarge.</p>
<p>Reviewer #2: Thank you for the opportunity to review this paper. This is an interesting manuscript presenting a results of a retrospective study that examine the impact of population inflow and outflow on sub- prefectural DID estimates in Japan.</p>
<p>My review mainly concerns only the statistical aspects of the study. Some questions reported below were raised and in my view, it is not acceptable in this form for the publication in this journal.</p>
<p>• The authors don’t’ compare median value of DID standardized between nighttime and daytime population. Are the difference statistically significant?</p>
<p>• In Figure 1 I think that the crosses represent the median and not the mean value as described in legend of figure 1.</p>
<p>• I suggest to the authors to estimate an indicator of agreement to evaluate the relationship between nighttime population–standardized DID values respect to daytime population–standardized DID values. (figure 2)</p>
<p>• Figure 4 is not included in the submission material</p>
<p>**********</p>
<p><!-- <font color="black"> -->6. PLOS authors have the option to publish the peer review history of their article (<ext-link ext-link-type="uri" xlink:href="https://journals.plos.org/plosone/s/editorial-and-peer-review-process#loc-peer-review-history" xlink:type="simple">what does this mean?</ext-link>). If published, this will include your full peer review and any attached files.</p>
<p>If you choose “no”, your identity will remain anonymous but your review may still be made public.</p>
<p><bold>Do you want your identity to be public for this peer review?</bold> For information about this choice, including consent withdrawal, please see our <ext-link ext-link-type="uri" xlink:href="https://www.plos.org/privacy-policy" xlink:type="simple">Privacy Policy</ext-link>.<!-- </font> --></p>
<p>Reviewer #1: No</p>
<p>Reviewer #2: No</p>
<p>[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]</p>
<p>While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, <ext-link ext-link-type="uri" xlink:href="https://pacev2.apexcovantage.com/" xlink:type="simple">https://pacev2.apexcovantage.com/</ext-link>. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at <email xlink:type="simple">figures@plos.org</email>. Please note that Supporting Information files do not need this step.</p>
</body>
</sub-article>
<sub-article article-type="author-comment" id="pone.0248338.r002">
<front-stub>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338.r002</article-id>
<title-group>
<article-title>Author response to Decision Letter 0</article-title>
</title-group>
<related-object document-id="10.1371/journal.pone.0248338" document-id-type="doi" document-type="peer-reviewed-article" id="rel-obj002" link-type="rebutted-decision-letter" object-id="10.1371/journal.pone.0248338.r001" object-id-type="doi" object-type="decision-letter"/>
<custom-meta-group>
<custom-meta>
<meta-name>Submission Version</meta-name>
<meta-value>1</meta-value>
</custom-meta>
</custom-meta-group>
</front-stub>
<body>
<p>
<named-content content-type="author-response-date">14 Dec 2020</named-content>
</p>
<p>Response to Reviewers</p>
<p>Reviewer #1:</p>
<p>Response: Thank you for your time and effort in reviewing our manuscript. We have revised the manuscript in accordance with your suggestions. </p>
<p>We have added a statistical analysis of the differences between the national-level DID values and the prefectural- and SMA-level DID values. This analysis showed that the SMA-level DID was significantly different from the national-level DID, irrespective of whether standardization was performed with the nighttime or daytime population. Therefore, as a major revision, we have softened the tone of our conclusion that we should use the daytime standardized–population as the SMA-level AMU’s denominator, and indicated that it may be better to avoid using SMA-level DID. Finally, we apologize for omitting Figure 4, and have included it in this resubmission. Thank you for giving us the chance to revise our manuscript, and we look forward to your evaluation.</p>
<p>Conclusions (Page 20, Lines 334-338)</p>
<p>Before</p>
<p>AMU surveillance has conventionally used the nighttime population for standardization. However, it is evident that nighttime population standardization can lead to erroneous estimates due to its omission of population inflow and outflow, especially for smaller regions. In Japan, the estimation of SMA-level DID values should be based on daytime population standardization to enable more impartial comparisons.</p>
<p>After</p>
<p>AMU surveillance has conventionally used the nighttime population for standardization. However, regional-level AMU estimates, especially of smaller regions such as SMAs, are more susceptible to the influence of whatever population is used for standardization, which can lead to erroneous estimates. Nevertheless, daytime populations that account for population inflow and outflow may be preferable to nighttime populations to minimize such errors and enable more precise comparisons.</p>
<p>1. The discussion is concise, but would benefit from a further elaboration of how this recommended shift from nighttime to daytime population usage will change antimicrobial use policy and management. What implications does this have for future regulation? What steps will be taken next? A more in-depth discussion would be helpful.</p>
<p>Thank you for your advice. As suggested, we have addressed the potential implications of our findings, although we prefer to refrain from a more in-depth discussion at this stage (we foresee that the shift from nighttime to daytime population could help produce more accurate AMU estimates, but anything downstream of that would be speculative). In this study, we identified population in/out flow to be one of the factors that can introduce bias into regional AMU estimates. As a next step, we aim to further elucidate the other factors that affect regional AMU, and to develop a standardization formula similar to the Standardized Antimicrobial Administration Ratio, which was developed by the US CDC and is used in hospital AMU comparisons. We have added this perspective to the end of the Discussion section.</p>
<p>Page 20, Lines 322-331</p>
<p>The shift from nighttime population to daytime population for standardization could provide more accurate regional-level AMU estimates, which would help to inform AMU-related policymaking, guide the efficient allocation of resources for region-specific antimicrobial stewardship programs, and contextualize the decision-making process for AMU regulations. As a next step, we aim to develop a formula to standardize regional-level AMU that can account for various sources of biases (e.g., population inflow and outflow) similar to the Standardized Antimicrobial Administration Ratio, which was developed by the US Centers for Disease Control and Prevention and is used in hospital-level AMU comparisons. In this way, our study represents the first step for developing a standardization method for regional-level AMU surveillance in Japan. Further studies are needed to identify the other factors of regional AMU surveillance.</p>
<p>Reference</p>
<p>14. van Santen K L, Edwards JR, Webb A K, Pollack L A, O’Leary E, Neuhauser M N, et al. The Standardized Antimicrobial Administration Ratio: A New Metric for Measuring and Comparing Antibiotic Use. Clin Infect Dis. 2018;67: 179-185. doi: 10.1093/cid/ciy075</p>
<p>2. I did not receive Figure 4 – the description is provided, but not the figure itself. This should be provided and reviewed.</p>
<p>We apologize for omitting Figure 4. This figure has been included in the resubmitted manuscript, and we would appreciate any comments or advice.</p>
<p>3. Figure 3 displays at a quite small size, making the color of the map and differences between region SMA difficult to see. Would enlarge.</p>
<p>Thank you for pointing this out. We have reformatted Figure 3 with a higher resolution (1200 dpi) and enlarged. We hope that it is more readable. </p>
<p>Reviewer #2:</p>
<p>Response: Thank you for your time and effort in reviewing our manuscript. We have revised the manuscript in accordance with your suggestions. </p>
<p>We have added a statistical analysis of the differences between the national-level DID values and the prefectural- and SMA-level DID values. This analysis showed that the SMA-level DID was significantly different from the national-level DID, irrespective of whether standardization was performed with the nighttime or daytime population. Therefore, as a major revision, we have softened the tone of our conclusion that we should use the daytime standardized–population as the SMA-level AMU’s denominator, and indicated that it may be better to avoid using SMA-level DID. Finally, we apologize for omitting Figure 4, and have included it in this resubmission. Thank you for giving us the chance to revise our manuscript, and we look forward to your evaluation.</p>
<p>Conclusions (Page 20, Lines 334-338)</p>
<p>Before</p>
<p>AMU surveillance has conventionally used the nighttime population for standardization. However, it is evident that nighttime population standardization can lead to erroneous estimates due to its omission of population inflow and outflow, especially for smaller regions. In Japan, the estimation of SMA-level DID values should be based on daytime population standardization to enable more impartial comparisons.</p>
<p>After</p>
<p>AMU surveillance has conventionally used the nighttime population for standardization. However, regional-level AMU estimates, especially of smaller regions such as SMAs, are more susceptible to the influence of whatever population is used for standardization, which can lead to erroneous estimates. Nevertheless, daytime populations that account for population inflow and outflow may be preferable to nighttime populations to minimize such errors and enable more precise comparisons.</p>
<p>1. The authors don’t’ compare median value of DID standardized between nighttime and daytime population. Are the differences statistically significant?</p>
<p>Thank you for your important advice. We have consulted with a statistician regarding this concern. Because prefectural- and SMA-level DID values do not have a “correct value” or gold standard, it would be difficult to interpret a statistically significant difference in median values between the nighttime and daytime populations. Instead, we have calculated the correlations between daytime and nighttime population–standardized DID in Figure 2. Furthermore, in Figure 1, we show the mean and 95% confidence interval for each indicator, and created violin plots to visualize their dispersions. Additionally, we compared each DID value to the national-level AMU using the one-sample t-test. The results showed that neither the daytime nor nighttime population–standardized DID at the prefectural level had any significant difference with the national-level AMU. However, both the daytime and nighttime population–standardized DID values at the SMA level were significantly different from the national-level AMU. Accordingly, we concluded that SMA-level DID values standardized using either of these populations would be difficult to evaluate. However, the 95% confidence interval was narrower and the mean (center of the widest section in the violin plot) was closer to the national AMU in the daytime population–standardized DID than in the nighttime population–standardized DID at the SMA level. Therefore, if SMA-level DID values must be calculated, the daytime population may represent the better option for standardization. We have added these points to the manuscript as follows:</p>
<p>Abstract (Page 3, Lines 38-42)</p>
<p>The national AMU was 17.21 DDDs per 1,000 population per day. The mean (95% confidence interval) prefectural-level DDDs per 1,000 nighttime and daytime population per day were 17.27 (14.10, 20.44) and 17.41 (14.30, 20.53), respectively. The mean (95% confidence interval) SMA-level DDDs per 1,000 nighttime and daytime population per day were 16.12 (9.84, 22.41) and 16.41 (10.57, 22.26), respectively.</p>
<p>Methods (Page 9, Lines 147 -150)</p>
<p>Next, we evaluated the distributions of daytime and nighttime population–standardized DID values at the prefectural and SMA levels using the Kolmogorov–Smirnov test. The mean DID values were compared with the national AMU using the one-sample t-test.</p>
<p>Results (Page 11, Lines 183-197)</p>
<p>When standardized with the nighttime population, the mean DID was 17.27 (95% confidence interval [14.10, 20.44]) at the prefectural level and 16.12 (95% confidence interval [9.84, 22.41]) at the SMA level. When standardized with the daytime population, the median DID was 17.41 (95% confidence interval [14.30, 20.53]) at the prefectural level and 16.41 (95% confidence interval [10.57, 22.26]) at the SMA level. Both the prefectural-level and SMA-level DID values were normally distributed regardless of the population used for standardization. The daytime and nighttime population–standardized mean DID values at the prefectural level were not significantly different from the national-level DID values (nighttime: P = .385; daytime: P = .811); however, the corresponding mean DID values at the SMA level were significantly different from the national-level DID values (nighttime: P &lt; .001; daytime: P &lt; .001). As shown in Fig 1, daytime population–standardized DID at the SMA level had a narrower dispersion and a mean value (center of the widest section in the violin plot) that was closer to the national-level DID than the nighttime population–standardized DID. In contrast, the daytime population and nighttime population–standardized DID at the prefectural level exhibited similar shapes in the violin plot.</p>
<p>Discussion (Page 17, Lines 269-279)</p>
<p>Previous studies have examined the effects of using different denominator values when calculating AMU in hospitals.9–12 However, these effects on regional AMU surveillance have not been explored. Even the World Health Organization’s AMU surveillance methodology does not address the appropriate methods for analyzing sub-national regions.13 In our analysis, SMA-level AMU (standardized using either the daytime or nighttime population) was found to be significantly different from the national-level AMU. This suggests that the population-standardized DID values of smaller regions are susceptible to the effects of population inflow and outflow, which can lead to erroneous results. </p>
<p>At the SMA level, daytime population–standardization produced fewer outliers, narrower 95% confidence intervals, and mean DID values that were closer to the national-level DID than nighttime population–standardization.</p>
<p>Discussion (Page 18, Lines 283-284)</p>
<p>Although standardization using either population may be problematic for SMAs, our findings indicate that the use of the daytime population is preferable to the nighttime population for such adjustments.</p>
<p>2. In Figure 1 I think that the crosses represent the median and not the mean value as described in legend of figure 1.</p>
<p>Thank you for the suggestion. As mentioned in the response to your previous comment, we have replaced the boxplots with violin plots.</p>
<p>3. I suggest to the authors to estimate an indicator of agreement to evaluate the relationship between nighttime population–standardized DID values respect to daytime population–standardized DID values. (figure 2)</p>
<p>As advised, we calculated Pearson’s correlation coefficients and P values for Figure 2. The following sentences have been added:</p>
<p>Methods (Page 9, Lines 152-153)</p>
<p>Correlations between nighttime and daytime population–standardized DID values were examined using Pearson’s correlation coefficients.</p>
<p>Results (Page 12, Lines 205-207)</p>
<p>The correlation coefficient between nighttime and daytime population–standardized DID values was higher at the prefectural level (0.90; P &lt; .001) than at the SMA level (0.80; P &lt; .001).</p>
<p>4. Figure 4 is not included in the submission material</p>
<p>We apologize for omitting Figure 4. This figure has been included in the resubmitted manuscript, and we would appreciate any comments or advice.</p>
<supplementary-material id="pone.0248338.s004" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.s004" xlink:type="simple">
<label>Attachment</label>
<caption>
<p>Submitted filename: <named-content content-type="submitted-filename">PONE-D-20-15330_R1_Response_to_Editors_and_Reviewers.docx</named-content></p>
</caption>
</supplementary-material>
</body>
</sub-article>
<sub-article article-type="aggregated-review-documents" id="pone.0248338.r003" specific-use="decision-letter">
<front-stub>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338.r003</article-id>
<title-group>
<article-title>Decision Letter 1</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name name-style="western">
<surname>Suppiah</surname>
<given-names>Vijayaprakash</given-names>
</name>
<role>Academic Editor</role>
</contrib>
</contrib-group>
<permissions>
<copyright-year>2021</copyright-year>
<copyright-holder>Vijayaprakash Suppiah</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<related-object document-id="10.1371/journal.pone.0248338" document-id-type="doi" document-type="article" id="rel-obj003" link-type="peer-reviewed-article"/>
<custom-meta-group>
<custom-meta>
<meta-name>Submission Version</meta-name>
<meta-value>1</meta-value>
</custom-meta>
</custom-meta-group>
</front-stub>
<body>
<p>
<named-content content-type="letter-date">4 Feb 2021</named-content>
</p>
<p>PONE-D-20-15330R1</p>
<p>Effect of population inflow and outflow between rural and urban areas on regional antimicrobial use surveillance</p>
<p>PLOS ONE</p>
<p>Dear Dr. Kusama,</p>
<p>Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.</p>
<p>Please submit your revised manuscript by Mar 21 2021 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at <email xlink:type="simple">plosone@plos.org</email>. When you're ready to submit your revision, log on to <ext-link ext-link-type="uri" xlink:href="https://www.editorialmanager.com/pone/" xlink:type="simple">https://www.editorialmanager.com/pone/</ext-link> and select the 'Submissions Needing Revision' folder to locate your manuscript file.</p>
<p>Please include the following items when submitting your revised manuscript:</p>
<p><list list-type="bullet"><list-item><p>A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.</p></list-item><list-item><p>A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.</p></list-item><list-item><p>An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.</p></list-item></list></p>
<p>If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.</p>
<p>If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: <ext-link ext-link-type="uri" xlink:href="http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols" xlink:type="simple">http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols</ext-link></p>
<p>We look forward to receiving your revised manuscript.</p>
<p>Kind regards,</p>
<p>Vijayaprakash Suppiah, PhD</p>
<p>Academic Editor</p>
<p>PLOS ONE</p>
<p>[Note: HTML markup is below. Please do not edit.]</p>
<p>Reviewers' comments:</p>
<p>Reviewer's Responses to Questions</p>
<p><!-- <font color="black"> --><bold>Comments to the Author</bold></p>
<p>1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.<!-- </font> --></p>
<p>Reviewer #3: (No Response)</p>
<p>**********</p>
<p><!-- <font color="black"> -->2. Is the manuscript technically sound, and do the data support the conclusions?</p>
<p>The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. <!-- </font> --></p>
<p>Reviewer #3: Partly</p>
<p>**********</p>
<p><!-- <font color="black"> -->3. Has the statistical analysis been performed appropriately and rigorously? <!-- </font> --></p>
<p>Reviewer #3: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->4. Have the authors made all data underlying the findings in their manuscript fully available?</p>
<p>The <ext-link ext-link-type="uri" xlink:href="http://www.plosone.org/static/policies.action#sharing" xlink:type="simple">PLOS Data policy</ext-link> requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.<!-- </font> --></p>
<p>Reviewer #3: No</p>
<p>**********</p>
<p><!-- <font color="black"> -->5. Is the manuscript presented in an intelligible fashion and written in standard English?</p>
<p>PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.<!-- </font> --></p>
<p>Reviewer #3: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->6. Review Comments to the Author</p>
<p>Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)<!-- </font> --></p>
<p>Reviewer #3: Summary: Understanding rates of anti-microbial use is an important part of health system surveillance, however, these rates may be sensitive to the population size used in the denominator, specifically whether it represents the nighttime or daytime population. Here the authors calculate AMU rates standardied for both populations at the prefecture and secondary medical area across the country of Japan, finding that these two types of metrics differ, particularly for urban regions with high levels of daily commuting.</p>
<p>Strengths:</p>
<p>The analysis seems to be done correctly, and the writing is clear and succinct.</p>
<p>The authors did a very comprehensive job of responding to previous comments from reviewers, however I have a couple of small comments.</p>
<p>Figure 2 is very clear and does an excellent job conveying the results.</p>
<p>Major Comments:</p>
<p>The authors make claims about which population value (nighttime vs. daytime) should be used to calculate AMU, based on which is more correct. However, the datasets used include no information about the "true" AMU, so it is unclear to me what this claim in based on. The authors do rightly make claims as to where this difference in calculation could lead to the widest discrepancies between the two metrics, but I don't understand how this leads to the conclusion that it is better to use daytime populations. Perhaps there is a well-supported assumption that people are prescribed medicine near their work and school places, but I am not familiar with the Japanese healthcare system and so do not know if this is the case. Other readers may also not be familiar with Japanese healthcare, and it would be helpful if the authors could provide more reasoning for the case they make.</p>
<p>Minor Comments:</p>
<p>Line 162: Do you mean the populations size of each age category for each SMA?</p>
<p>Line 186: Why do you compare median here and mean on line 184? Is one of these a typo, since you do compare them later in this paragraph?</p>
<p>Table 1: What do the three different tables in Table 1 represent? Is it possible to add a sub-title, such as a row that seperates each table that lables what each one is?</p>
<p>Line 284: It is not clear to me why the daytime-population standardized DID is better (see comment above). Could you explain this further?</p>
<p>Fig. 3: Please label the color bar with what it represents.</p>
<p>Data Availability:</p>
<p>The authors have made the population data available in the supplement, but the data on AMU is not available, which is required to conduct the analyses.</p>
<p>**********</p>
<p><!-- <font color="black"> -->7. PLOS authors have the option to publish the peer review history of their article (<ext-link ext-link-type="uri" xlink:href="https://journals.plos.org/plosone/s/editorial-and-peer-review-process#loc-peer-review-history" xlink:type="simple">what does this mean?</ext-link>). If published, this will include your full peer review and any attached files.</p>
<p>If you choose “no”, your identity will remain anonymous but your review may still be made public.</p>
<p><bold>Do you want your identity to be public for this peer review?</bold> For information about this choice, including consent withdrawal, please see our <ext-link ext-link-type="uri" xlink:href="https://www.plos.org/privacy-policy" xlink:type="simple">Privacy Policy</ext-link>.<!-- </font> --></p>
<p>Reviewer #3: No</p>
<p>[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]</p>
<p>While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, <ext-link ext-link-type="uri" xlink:href="https://pacev2.apexcovantage.com/" xlink:type="simple">https://pacev2.apexcovantage.com/</ext-link>. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at <email xlink:type="simple">figures@plos.org</email>. Please note that Supporting Information files do not need this step.</p>
</body>
</sub-article>
<sub-article article-type="author-comment" id="pone.0248338.r004">
<front-stub>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338.r004</article-id>
<title-group>
<article-title>Author response to Decision Letter 1</article-title>
</title-group>
<related-object document-id="10.1371/journal.pone.0248338" document-id-type="doi" document-type="peer-reviewed-article" id="rel-obj004" link-type="rebutted-decision-letter" object-id="10.1371/journal.pone.0248338.r003" object-id-type="doi" object-type="decision-letter"/>
<custom-meta-group>
<custom-meta>
<meta-name>Submission Version</meta-name>
<meta-value>2</meta-value>
</custom-meta>
</custom-meta-group>
</front-stub>
<body>
<p>
<named-content content-type="author-response-date">17 Feb 2021</named-content>
</p>
<p>Response to Reviewers</p>
<p>Reviewer #3:</p>
<p>Thank you for your useful advice. We have corrected our manuscript according to your suggestions. We hope that our revisions meet your requirements.</p>
<p>Major Comments:</p>
<p>The authors make claims about which population value (nighttime vs. daytime) should be used to calculate AMU, based on which is more correct. However, the datasets used include no information about the "true" AMU, so it is unclear to me what this claim in based on. The authors do rightly make claims as to where this difference in calculation could lead to the widest discrepancies between the two metrics, but I don't understand how this leads to the conclusion that it is better to use daytime populations. Perhaps there is a well-supported assumption that people are prescribed medicine near their work and school places, but I am not familiar with the Japanese healthcare system and so do not know if this is the case. Other readers may also not be familiar with Japanese healthcare, and it would be helpful if the authors could provide more reasoning for the case they make.</p>
<p>Thank you for pointing this out. We agree with your opinion. Although we described that daytime population standardization is a superior method for estimating AMU in small regions than nighttime population standardization, this is not justified because of the absence of a standard AMU, as you pointed out. Therefore, we have corrected the Discussion and Conclusions as follows: </p>
<p>Abstract (Lines 47-49)</p>
<p>Before</p>
<p>Regional-level AMU estimates, especially of smaller regions such as SMAs, are susceptible to the use of different populations for standardization. Daytime populations that account for population in/outflow may be preferable to nighttime populations for such adjustments.</p>
<p>After</p>
<p>Regional-level AMU estimates, especially of smaller regions such as SMAs, are susceptible to the use of different populations for standardization. This finding indicates that AMU standardization based on population values is not suitable for AMU estimates in small regions. </p>
<p>Discussion (Line 277-284)</p>
<p>Before</p>
<p>At the SMA level, daytime population–standardization produced fewer outliers, narrower 95% confidence intervals, and mean DID values that were closer to the national-level DID than nighttime population–standardization. These findings showed that when analyzing smaller regional units such as SMAs, the use of different populations for standardization had a considerable effect on DID estimates in central urban areas and their surrounding regions. As more medical examinations and prescriptions are received during the day than at night, the calculation of regional AMU indices should account for the effects of population inflow and outflow. Although standardization using either population may be problematic for SMAs, our findings indicate that the use of the daytime population is preferable to the nighttime population for such adjustments.</p>
<p>After</p>
<p>As more medical examinations and prescriptions are received during the day than at night, the calculation of regional AMU indices should account for the effects of population inflow and outflow. At the SMA level, daytime population standardization produced fewer outliers, narrower 95% confidence intervals, and mean DID values that were closer to the national DID than nighttime population standardization. These findings showed that when analyzing smaller regional units, such as SMAs, the use of different populations for standardization has a considerable effect on DID estimates in central urban areas and their surrounding regions. Our insights indicate that AMU standardization based on population values is not suitable for AMU estimates in small regions.</p>
<p>Discussion (Line 319-322)</p>
<p>Before</p>
<p>The shift from nighttime population to daytime population for standardization could provide more accurate regional-level AMU estimates, which would help to inform AMU-related policymaking, guide the efficient allocation of resources for region-specific antimicrobial stewardship programs, and contextualize the decision-making process for AMU regulations. As a next step, we aim to develop a formula to standardize regional-level AMU that can account for various sources of biases (e.g., population inflow and outflow) similar to the Standardized Antimicrobial Administration Ratio, which was developed by the US Centers for Disease Control and Prevention and is used in hospital-level AMU comparisons. In this way, our study represents the first step for developing a standardization method for regional-level AMU surveillance in Japan. Further studies are needed to identify the other factors of regional AMU surveillance.</p>
<p>After</p>
<p>Although regional AMU estimates would help inform AMU-related policymaking, we recommend avoiding AMU evaluation in small regions standardized by the population number. If we wish to monitor AMU evaluations in small regions, only the temporal change in a region should be considered, for example. </p>
<p>Conclusion (Line 328-329)</p>
<p>Before</p>
<p>Nevertheless, daytime populations that account for population inflow and outflow may be preferable to nighttime populations to minimize such errors and enable more precise comparisons.</p>
<p>After</p>
<p>Therefore, new approaches are required to monitor AMU evaluations in small regions, for example, observing only the temporal changes in a region.</p>
<p>Minor Comments:</p>
<p>1. Line 162: Do you mean the populations size of each age category for each SMA?</p>
<p>Yes, we used population numbers according to the age categories. We have corrected the description as you indicated.</p>
<p>Methods (Line 160-161)</p>
<p>Before</p>
<p>Finally, we divided each SMA’s population into three categories based on age (children: &lt;15 years, working-age adults: 15–64 years, and older adults: ≥65 years), and analyzed the correlation between the absolute difference in DID values and each age category. For this analysis, the Tokyo/Ku-chuoubu SMA was again excluded as it was an extreme outlier.</p>
<p>After</p>
<p>Finally, we used the population number for each age category (children: &lt;15 years, working-age adults: 15–64 years, and older adults: ≥65 years) for each SMA, and analyzed the correlation between the absolute difference in DID values and each age category. For this analysis, the Tokyo/Ku-chuoubu SMA was excluded as it was an extreme outlier.</p>
<p>2. Line 186: Why do you compare median here and mean on line 184? Is one of these a typo, since you do compare them later in this paragraph?</p>
<p>Thank you for pointing this out. We corrected “median” to “mean.”</p>
<p>3. Table 1: What do the three different tables in Table 1 represent? Is it possible to add a sub-title, such as a row that seperates each table that lables what each one is?</p>
<p>We added descriptions to each table. </p>
<p>4. Line 284: It is not clear to me why the daytime-population standardized DID is better (see comment above). Could you explain this further?</p>
<p>Thank you. We corrected our discussion and conclusions as presented above.</p>
<p>5. Fig. 3: Please label the color bar with what it represents.</p>
<p>Thank you. We have added explanations of the color bars. </p>
<p>6. The authors have made the population data available in the supplement, but the data on AMU is not available, which is required to conduct the analyses.</p>
<p>Thank you. We have created Supplementary Table 3 to display AMU data.</p>
<p>Methods (Line 155-156)</p>
<p>each DID value was shown in Table S3</p>
<supplementary-material id="pone.0248338.s005" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" position="float" xlink:href="info:doi/10.1371/journal.pone.0248338.s005" xlink:type="simple">
<label>Attachment</label>
<caption>
<p>Submitted filename: <named-content content-type="submitted-filename">PONE-D-20-15330_R2_Response_to_reviewer.docx</named-content></p>
</caption>
</supplementary-material>
</body>
</sub-article>
<sub-article article-type="aggregated-review-documents" id="pone.0248338.r005" specific-use="decision-letter">
<front-stub>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338.r005</article-id>
<title-group>
<article-title>Decision Letter 2</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name name-style="western">
<surname>Suppiah</surname>
<given-names>Vijayaprakash</given-names>
</name>
<role>Academic Editor</role>
</contrib>
</contrib-group>
<permissions>
<copyright-year>2021</copyright-year>
<copyright-holder>Vijayaprakash Suppiah</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<related-object document-id="10.1371/journal.pone.0248338" document-id-type="doi" document-type="article" id="rel-obj005" link-type="peer-reviewed-article"/>
<custom-meta-group>
<custom-meta>
<meta-name>Submission Version</meta-name>
<meta-value>2</meta-value>
</custom-meta>
</custom-meta-group>
</front-stub>
<body>
<p>
<named-content content-type="letter-date">25 Feb 2021</named-content>
</p>
<p>Effect of population inflow and outflow between rural and urban areas on regional antimicrobial use surveillance</p>
<p>PONE-D-20-15330R2</p>
<p>Dear Dr. Kusama,</p>
<p>We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.</p>
<p>Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.</p>
<p>An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at <ext-link ext-link-type="uri" xlink:href="http://www.editorialmanager.com/pone/" xlink:type="simple">http://www.editorialmanager.com/pone/</ext-link>, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at <email xlink:type="simple">authorbilling@plos.org</email>.</p>
<p>If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact <email xlink:type="simple">onepress@plos.org</email>.</p>
<p>Kind regards,</p>
<p>Vijayaprakash Suppiah, PhD</p>
<p>Academic Editor</p>
<p>PLOS ONE</p>
<p>Reviewers' comments:</p>
<p>Reviewer's Responses to Questions</p>
<p><!-- <font color="black"> --><bold>Comments to the Author</bold></p>
<p>1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.<!-- </font> --></p>
<p>Reviewer #3: All comments have been addressed</p>
<p>**********</p>
<p><!-- <font color="black"> -->2. Is the manuscript technically sound, and do the data support the conclusions?</p>
<p>The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. <!-- </font> --></p>
<p>Reviewer #3: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->3. Has the statistical analysis been performed appropriately and rigorously? <!-- </font> --></p>
<p>Reviewer #3: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->4. Have the authors made all data underlying the findings in their manuscript fully available?</p>
<p>The <ext-link ext-link-type="uri" xlink:href="http://www.plosone.org/static/policies.action#sharing" xlink:type="simple">PLOS Data policy</ext-link> requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.<!-- </font> --></p>
<p>Reviewer #3: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->5. Is the manuscript presented in an intelligible fashion and written in standard English?</p>
<p>PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.<!-- </font> --></p>
<p>Reviewer #3: Yes</p>
<p>**********</p>
<p><!-- <font color="black"> -->6. Review Comments to the Author</p>
<p>Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)<!-- </font> --></p>
<p>Reviewer #3: The authors have addressed all of my comments. I find the new paragraph in the discussion regarding how daytime population numbers affect the AMU indices especially helpful and I think it really ties the manuscript together. Great work!</p>
<p>**********</p>
<p><!-- <font color="black"> -->7. PLOS authors have the option to publish the peer review history of their article (<ext-link ext-link-type="uri" xlink:href="https://journals.plos.org/plosone/s/editorial-and-peer-review-process#loc-peer-review-history" xlink:type="simple">what does this mean?</ext-link>). If published, this will include your full peer review and any attached files.</p>
<p>If you choose “no”, your identity will remain anonymous but your review may still be made public.</p>
<p><bold>Do you want your identity to be public for this peer review?</bold> For information about this choice, including consent withdrawal, please see our <ext-link ext-link-type="uri" xlink:href="https://www.plos.org/privacy-policy" xlink:type="simple">Privacy Policy</ext-link>.<!-- </font> --></p>
<p>Reviewer #3: No</p>
</body>
</sub-article>
<sub-article article-type="editor-report" id="pone.0248338.r006" specific-use="acceptance-letter">
<front-stub>
<article-id pub-id-type="doi">10.1371/journal.pone.0248338.r006</article-id>
<title-group>
<article-title>Acceptance letter</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name name-style="western">
<surname>Suppiah</surname>
<given-names>Vijayaprakash</given-names>
</name>
<role>Academic Editor</role>
</contrib>
</contrib-group>
<permissions>
<copyright-year>2021</copyright-year>
<copyright-holder>Vijayaprakash Suppiah</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<related-object document-id="10.1371/journal.pone.0248338" document-id-type="doi" document-type="article" id="rel-obj006" link-type="peer-reviewed-article"/>
</front-stub>
<body>
<p>
<named-content content-type="letter-date">2 Mar 2021</named-content>
</p>
<p>PONE-D-20-15330R2 </p>
<p>Effect of population inflow and outflow between rural and urban areas on regional antimicrobial use surveillance </p>
<p>Dear Dr. Kusama:</p>
<p>I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department. </p>
<p>If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact <email xlink:type="simple">onepress@plos.org</email>.</p>
<p>If we can help with anything else, please email us at <email xlink:type="simple">plosone@plos.org</email>. </p>
<p>Thank you for submitting your work to PLOS ONE and supporting open access. </p>
<p>Kind regards, </p>
<p>PLOS ONE Editorial Office Staff</p>
<p>on behalf of</p>
<p>Dr. Vijayaprakash Suppiah  </p>
<p>Academic Editor</p>
<p>PLOS ONE</p>
</body>
</sub-article>
</article>