<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Vet. Sci.</journal-id>
<journal-title>Frontiers in Veterinary Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Vet. Sci.</abbrev-journal-title>
<issn pub-type="epub">2297-1769</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fvets.2023.1208804</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Veterinary Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Factors associated with heartworm preventative use in the golden retriever lifetime study</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes"><name><surname>Wisnieski</surname><given-names>Lauren</given-names></name><xref rid="c001" ref-type="corresp"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1415461/overview"/>
</contrib>
<contrib contrib-type="author"><name><surname>Faulkner</surname><given-names>Vina</given-names></name>
</contrib>
<contrib contrib-type="author"><name><surname>Faulkner</surname><given-names>Charles</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/2316111/overview"/>
</contrib>
</contrib-group>
<aff><institution>Richard A. Gillespie College of Veterinary Medicine, Lincoln Memorial University</institution>, <addr-line>Harrogate, TN</addr-line>, <country>United States</country></aff>
<author-notes>
<fn id="fn0001" fn-type="edited-by"><p>Edited by: Javier Sanchez, University of Prince Edward Island, Canada</p></fn>
<fn id="fn0002" fn-type="edited-by"><p>Reviewed by: Margaret Ruth Slater, Texas A and M University, United States; Kennedy Kapala Mwacalimba, Zoetis, United States</p></fn>
<corresp id="c001">&#x002A;Correspondence: Lauren Wisnieski, <email>lauren.wisnieski@lmunet.edu</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>09</day>
<month>06</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1208804</elocation-id>
<history>
<date date-type="received">
<day>19</day>
<month>04</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>19</day>
<month>05</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Wisnieski, Faulkner and Faulkner.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Wisnieski, Faulkner and Faulkner</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>Heartworm disease is preventable with use of heartworm preventatives, but the reported prevalence of heartworm preventative use in the United States is low, some estimates falling around 50% of dogs. However, there are very few estimates of prevalence and its associated factors.</p>
</sec>
<sec>
<title>Methods</title>
<p>We aimed to estimate prevalence and evaluate factors, including vaccination status, demographics, lifestyle, physical conditions, medications and supplements, and environment and living conditions, for their association with heartworm preventative use in a large dataset from the Golden Retriever Lifetime Study (<italic>N</italic> = 2,998). Due to the large number of predictors evaluated, we built a bootstrapped elastic net logistic regression model, which is robust to overfitting and multicollinearity. Variables were evaluated by calculating covariate stability (&#x003E;80%) and statistical significance (<italic>p</italic>&#x003C;0.02).</p>
</sec>
<sec>
<title>Results</title>
<p>In our sample, the prevalence of heartworm use was 39.5%. In our elastic net model, receiving vaccinations (rabies, Bordetella, or any other vaccine), being located in the Southern U.S., being altered, having an infectious disease or ear/ nose/throat system disease diagnosis, being on heartworm preventatives in the past, currently being on tick preventative, having sun exposure in an area with concrete flooring, living in a house with more rooms with carpeted floors, and spending time on hardwood flooring inside were associated with greater odds of heartworm preventative use. Supplementation use and being in the top quartile of height were associated with lower odds of heartworm preventative use.</p>
</sec>
<sec>
<title>Discussion</title>
<p>The explanatory factors we identified can be used to improve client communication. In addition, target populations for educational interventions and outreach can be identified. Future studies can validate the findings in a more diverse population of dogs.</p>
</sec>
</abstract>
<kwd-group>
<kwd>elastic net regression</kwd>
<kwd>heartworm disease</kwd>
<kwd>prophylactic</kwd>
<kwd>bootstrapping</kwd>
<kwd>canine</kwd>
<kwd>preventative</kwd>
<kwd>dog</kwd>
<kwd>big data</kwd>
</kwd-group>
<counts>
<fig-count count="1"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="19"/>
<page-count count="8"/>
<word-count count="5225"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Veterinary Epidemiology and Economics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="sec1" sec-type="intro"><label>1.</label>
<title>Introduction</title>
<p>Canine heartworm disease (CHD) incidence and geographic range is increasing in the United States (<xref ref-type="bibr" rid="ref1">1</xref>, <xref ref-type="bibr" rid="ref2">2</xref>). From 2013 to 2016, Drake and Wiseman (<xref ref-type="bibr" rid="ref1">1</xref>) estimated that the relative change in CHD was +15.3% in the US, with highest increases in the Southeast region. Heartworm disease is preventable with the use of prophylactic macrolide medications, which are safe and effective (<xref ref-type="bibr" rid="ref3">3</xref>). Previous studies have estimated that adherence to the recommended 12-month use of heartworm preventative is low (24.4% for monthly HWD preventatives) (<xref ref-type="bibr" rid="ref4">4</xref>, <xref ref-type="bibr" rid="ref5">5</xref>). One study found that use of heartworm preventatives was estimated to be approximately 74&#x2013;79% of dogs presenting to a veterinary teaching hospital, however it was reduced to 50% in winter (<xref ref-type="bibr" rid="ref5">5</xref>). Very few studies have evaluated what factors are associated with heartworm preventative use in dogs. Bir et al. (<xref ref-type="bibr" rid="ref6">6</xref>) found greater heartworm preventative use among dogs with owners that were male, younger, lived in the South, and had higher socioeconomic status. Gates and Nolan (<xref ref-type="bibr" rid="ref5">5</xref>) evaluated multiple factors included signalment, presenting complaint, and vaccination status and found that only month of presentation and neuter status were associated with heartworm preventative use. Due to the sparse research on this topic, there is a need to identify additional factors associated with heartworm preventative use, including lifestyle, the strength of the human-animal bond, the environment, and living conditions. Therefore, the objective of this study was to evaluate the association of multiple factors with heartworm preventative use utilizing data available from the Golden Retriever Lifetime Study, which is a large prospective study that aims to identify risk factors for disease in golden retrievers (<xref ref-type="bibr" rid="ref7">7</xref>). Due to the large number of predictors evaluated, we utilized elastic net regression, which is a regularized regression method that reduces risk of overfitting and multicollinearity (<xref ref-type="bibr" rid="ref8">8</xref>).</p>
</sec>
<sec id="sec2" sec-type="materials|methods"><label>2.</label>
<title>Materials and methods</title>
<p>The Institutional Review Board at Lincoln Memorial University approved this study (1,068&#x2009;V.0). All data cleaning and analyses were conducted in Stata version 17.0 (College Station, TX).</p>
<sec id="sec3"><label>2.1.</label>
<title>Dataset and data cleaning</title>
<p>The Golden Retriever Lifetime study is a large, prospective cohort study that was established in 2012, which investigates risk factors for a variety of health conditions, including the primary endpoints of lymphoma, hemangiosarcoma, high-grade mast cell tumors, and osteosarcoma. Over each dog&#x2019;s lifetime, owners and veterinarians complete an annual questionnaire and clinical exam and collect samples. All dogs in the study visit their veterinarian yearly to complete a physical exam and sample collection (<xref ref-type="bibr" rid="ref7">7</xref>). For this analysis, all data from the owner and veterinarian questionnaires were requested, except for diet. The full list of variables screened can be obtained from the Annual Veterinarian Questionnaire at<xref rid="fn0003" ref-type="fn"><sup>1</sup></xref> and the Annual Owner Questionnaire at.<xref rid="fn0004" ref-type="fn"><sup>2</sup></xref> After receiving the data, variables were further screened for feasibility of a potential relationship with heartworm prophylaxis use by our study team. The remaining variables underwent data cleaning. Numerical variables were evaluated for linearity. Variables that violated linearity requirements and that were unable to be transformed (ex. Logarithm transformation) were converted to categorical variables based on quartiles of the original variable. Categorical variables with more than two categories were converted to dummy variables as required for elastic net regression. Variables with a large amount of missing data or lack of variation were removed from the analysis. Complete case analysis was used, so observations with any missing data of the remaining variables were removed. The final list of predictor variables and how they were coded are presented in <xref rid="tab1" ref-type="table">Table 1</xref>.</p>
<table-wrap position="float" id="tab1"><label>Table 1</label>
<caption>
<p>Variable names, descriptions, and coding used for a bootstrapped elastic net model with binary outcome variable &#x201C;heartworm preventative use.&#x201D;</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Category</th>
<th align="left" valign="top">Variable</th>
<th align="left" valign="top">Description</th>
<th align="left" valign="top">Coding</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top" rowspan="3">Demographics</td>
<td align="left" valign="top">Coat color</td>
<td align="left" valign="top">Dog&#x2019;s reported coat color</td>
<td align="left" valign="top">Converted to 4 dummy variables for blonde, dark gold or red, medium gold, and light gold (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Region</td>
<td align="left" valign="top">Region of residence reported on the AOQ</td>
<td align="left" valign="top">Converted to 5 dummy variables for South, Midwest, Mountain, Northeast, and Pacific (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Sex</td>
<td align="left" valign="top">Dog&#x2019;s reported sex status</td>
<td align="left" valign="top">Converted to 4 dummy variables for male neutered, male intact, female neutered, and female intact (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top" rowspan="14">Lifestyle</td>
<td align="left" valign="top">Insurance</td>
<td align="left" valign="top">Dog health insurance status</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Primary activity</td>
<td align="left" valign="top">Primacy activity/lifestyle</td>
<td align="left" valign="top">Converted to 1 dummy variables for companion animal (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Brushes teeth</td>
<td align="left" valign="top">Dog gets teeth brushed</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Dental treat</td>
<td align="left" valign="top">Dog gets dental treats</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Dental food</td>
<td align="left" valign="top">Dog gets dental food</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Mouthwash</td>
<td align="left" valign="top">Dog gets mouthwash/water additive</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Homegroom frequency</td>
<td align="left" valign="top">Homegroom &#x003E;1x a month</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Professional groom</td>
<td align="left" valign="top">Ever professional groom</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Sun exposure</td>
<td align="left" valign="top">Dog&#x2019;s sun exposure duration</td>
<td align="left" valign="top">1&#x2009;=&#x2009;3 or more hours, 0&#x2009;=&#x2009;less than 3&#x2009;h</td>
</tr>
<tr>
<td align="left" valign="top">Sun exposure location</td>
<td align="left" valign="top">Dog&#x2019;s sun exposure location</td>
<td align="left" valign="top">Converted to 3 dummy variables for concrete, grass, and dirt (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Swim location</td>
<td align="left" valign="top">Dog&#x2019;s swim location</td>
<td align="left" valign="top">Converted to 4 dummy variables for swimming pool, pond/lake, river/stream/agricultural ditch, and ocean (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Attachment 1</td>
<td align="left" valign="top">Dogs tends to sit close to or in contact with owner (or others) when owner (or others) are sitting down</td>
<td align="left" valign="top">Left as numerical scale (1&#x2009;=&#x2009;never, 2&#x2009;=&#x2009;seldom, 3&#x2009;=&#x2009;sometimes, 4&#x2009;=&#x2009;usually, 5&#x2009;=&#x2009;always)</td>
</tr>
<tr>
<td align="left" valign="top">Attachment 2</td>
<td align="left" valign="top">Dogs tends to nudge, nuzzle, or paw owner (or others) for attention when owner (or others) are sitting down</td>
<td align="left" valign="top">Left as numerical scale (1&#x2009;=&#x2009;never, 2&#x2009;=&#x2009;seldom, 3&#x2009;=&#x2009;sometimes, 4&#x2009;=&#x2009;usually, 5&#x2009;=&#x2009;always)</td>
</tr>
<tr>
<td align="left" valign="top">Dog activity level</td>
<td align="left" valign="top">Dog is active, energetic, always on the go</td>
<td align="left" valign="top">Left as numerical scale (1&#x2009;=&#x2009;never, 2&#x2009;=&#x2009;seldom, 3&#x2009;=&#x2009;sometimes, 4&#x2009;=&#x2009;usually, 5&#x2009;=&#x2009;always)</td>
</tr>
<tr>
<td align="left" valign="top" rowspan="5">AOQ medications and supplements</td>
<td align="left" valign="top">Supplements</td>
<td align="left" valign="top">Supplements given</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Medication</td>
<td align="left" valign="top">Medications given</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Tick</td>
<td align="left" valign="top">Dog is on tick preventative</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Prior heartworm</td>
<td align="left" valign="top">Dog has been on heartworm preventative in past</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Heartworm</td>
<td align="left" valign="top">Dog is currently on heartworm preventative</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top" rowspan="22">Environment and living conditions</td>
<td align="left" valign="top">Wood flooring</td>
<td align="left" valign="top">Number of rooms with wood flooring</td>
<td align="left" valign="top">0, 1, 2&#x2026;</td>
</tr>
<tr>
<td align="left" valign="top">Carpeted flooring</td>
<td align="left" valign="top">Number of rooms with carpeted flooring</td>
<td align="left" valign="top">0, 1, 2&#x2026;</td>
</tr>
<tr>
<td align="left" valign="top">Tile flooring</td>
<td align="left" valign="top">Number of rooms with tile flooring</td>
<td align="left" valign="top">0, 1, 2&#x2026;</td>
</tr>
<tr>
<td align="left" valign="top">Linoleum flooring</td>
<td align="left" valign="top">Rooms with linoleum flooring</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends most time</td>
<td align="left" valign="top">Where dog spends most time</td>
<td align="left" valign="top">Converted to 2 dummy variables for indoors, indoors/outdoors and outdoors (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Crate</td>
<td align="left" valign="top">Dog spends time indoors in indoor crate</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends time on carpeted flooring</td>
<td align="left" valign="top">Dog spends time indoors in a room with carpeted flooring</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends time on hardwood flooring</td>
<td align="left" valign="top">Dog spends time indoors in a room with hardwood flooring</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends time on tile flooring</td>
<td align="left" valign="top">Dog spends time indoors in a room with tile flooring</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends time on linoleum flooring</td>
<td align="left" valign="top">Dog spends time indoors in a room with linoleum flooring</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends time in room with furniture</td>
<td align="left" valign="top">Dog spends time indoors in a room with furniture</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Spends time in another room</td>
<td align="left" valign="top">Dog spends time indoors in another room</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Kennel</td>
<td align="left" valign="top">Dog spends time in kennel</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Fenced area</td>
<td align="left" valign="top">Dog spends time in fenced area</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Chain or lead</td>
<td align="left" valign="top">Dog spends time on a chain or lead</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Poison control</td>
<td align="left" valign="top">Number of poison control or veterinarian calls due to dog ingesting poison or other hazardous material in the past 12&#x2009;months</td>
<td align="left" valign="top">Converted to 1 dummy variable (0&#x2009;=&#x2009;never, 1&#x2009;=&#x2009;at least once)</td>
</tr>
<tr>
<td align="left" valign="top">Smoke</td>
<td align="left" valign="top">Dog exposed to secondhand smoke in the last 12&#x2009;months</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Neighborhood type</td>
<td align="left" valign="top">Primary neighborhood type dog has lived in the past 12&#x2009;months</td>
<td align="left" valign="top">Converted to 3 dummy variables for urban, suburban, and rural (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Single family home</td>
<td align="left" valign="top">Primary home dog has lived in the past 12&#x2009;months is a single family home</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Water source</td>
<td align="left" valign="top">Primary address water source</td>
<td align="left" valign="top">Converted to 2 dummy variables for municipal water and well water (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Central A/C</td>
<td align="left" valign="top">Primary address central A/C status</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Secondary home</td>
<td align="left" valign="top">Dog has a secondary address</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top" rowspan="5">Physical conditions</td>
<td align="left" valign="top">Height</td>
<td align="left" valign="top">Dog&#x2019;s height in cm</td>
<td align="left" valign="top">Converted to 4 dummy variables based on quartiles: &#x003C; 54.6, 54.6 - &#x003C;57.2, 57.2 - &#x003C;63.5, 63.5 and higher (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">Weight</td>
<td align="left" valign="top">Dog&#x2019;s weight in kg</td>
<td align="left" valign="top">Converted to 4 dummy variables based on quartiles: &#x003C;24.4, 24.4 - &#x003C;27.7, 27.7 - &#x003C;31.5, 31.5 and greater (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top">BCS</td>
<td align="left" valign="top">Purina BCS</td>
<td align="left" valign="top">1, 2, 3,&#x2026;9</td>
</tr>
<tr>
<td align="left" valign="top">Parasite</td>
<td align="left" valign="top">Dog diagnosed with gastrointestinal parasite past 12&#x2009;months</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Body system diagnosis</td>
<td align="left" valign="top">Whether the dog had any diagnoses in body system category</td>
<td align="left" valign="top">Converted to 9 dummy variables for ear/nose/throat, eye, gastrointestinal, infectious disease, musculoskeletal, reproductive, skin, urinary, and other (1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no)</td>
</tr>
<tr>
<td align="left" valign="top" rowspan="4">Vaccines</td>
<td align="left" valign="top">Bordetella</td>
<td align="left" valign="top">Dog received Bordetella vaccine</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Rabies</td>
<td align="left" valign="top">Dog received rabies vaccine</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Other vaccine</td>
<td align="left" valign="top">Dog received a different vaccine than Bordetella and rabies</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
<tr>
<td align="left" valign="top">Vaccine</td>
<td align="left" valign="top">Dog received any vaccine</td>
<td align="left" valign="top">1&#x2009;=&#x2009;yes, 0&#x2009;=&#x2009;no</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="sec4"><label>2.2.</label>
<title>Statistical analysis</title>
<p>Traditional regression models are prone to overfitting and inflated standard errors when a large number of predictors are used. Regularization methods, also called penalized regression models, reduce risk of overfitting by constraining the estimated coefficients and decreasing out of sample error. Elastic net models are a type of penalized regression technique that combines ridge regression and lasso regression. Elastic net models incorporate the ridge regression penalty which can handle correlated predictors and the feature selection characteristic of lasso regression (<xref ref-type="bibr" rid="ref8">8</xref>). Elastic net logistic regression models were built with the &#x2018;<italic>elasticnet</italic>&#x2019; package.</p>
<p>A bootstrapping procedure proposed by Bunea et al. (<xref ref-type="bibr" rid="ref9">9</xref>) was used to calculate covariate stability. Bootstrapped value of ps and confidence intervals were calculated as previously reported for robust estimation of model parameters (<xref ref-type="bibr" rid="ref10 ref11 ref12">10&#x2013;12</xref>). This procedure involved performing 10-fold cross-validation 10 times in 500 bootstrapped samples, created using the <italic>&#x2018;bsample&#x2019;</italic> command. The optimal tuning parameters, alpha and lambda, were selected through this process (<xref ref-type="bibr" rid="ref8">8</xref>). A large range of alpha values was tested (0.1&#x2013;1.0 at 0.1 intervals). A large range of lambdas was also tested using the &#x2018;<italic>selection</italic> (cv.<italic>, alllambdas</italic>)&#x2019; option. Covariate stability was defined as the percent of time the covariate was selected to be in the model in the 500 bootstrapped samples. Bootstrapped value of ps were calculated as one minus the percent of time the coefficient estimates were on the majority side of zero. For example, if the coefficient estimate was greater than zero 450 out of 500 times, the bootstrapped value of <italic>p</italic> would be 1 &#x2013; (450/500)&#x2009;=&#x2009;0.10. Lastly, 90% bootstrapped confidence intervals were calculated using the &#x2018;<italic>bootstrap&#x2019;</italic> command. A plot of covariates stabilities against significance was used to determine thresholds for including variables in the final model. Variables were included in the final model if the bootstrapped value of <italic>p</italic> was &#x003C;0.02 and the covariate stability was &#x003E;80% (<xref rid="fig1" ref-type="fig">Figure 1</xref>) based on separation of the datapoints on the plot.</p>
<fig position="float" id="fig1"><label>Figure 1</label>
<caption>
<p><italic>p</italic>-value versus covariate stability of variables in a bootstrapped elastic net model.</p>
</caption>
<graphic xlink:href="fvets-10-1208804-g001.tif"/>
</fig>
</sec>
</sec>
<sec id="sec5" sec-type="results"><label>3.</label>
<title>Results</title>
<sec id="sec6"><label>3.1.</label>
<title>Dataset description</title>
<p>The raw baseline dataset included data from 3,044 golden retrievers. After removing 46 records that had missing data in the final set of variables, 2,998 golden retrievers remained for analysis. Less than half (39.5%) of dogs were currently on heartworm preventative. The sample was approximately half female (49.5%) and male (50.5%). Less than half of males (32.9%) and females (42.0%) were neutered and spayed, respectively. Although age was not used in the analysis due to a high degree of missingness (29.7%), the mean age of the dogs in the 70.3% that reported age was 1.3&#x2009;years, with a minimum and maximum value of approximately 5&#x2009;months and 3&#x2009;years, respectively. The American Heartworm Society recommends starting puppies on preventative as early as possible, no later than 8&#x2009;weeks of age, so ideally all dogs in the sample would be on preventative.<xref rid="fn0005" ref-type="fn"><sup>3</sup></xref> The sample represented the entire U.S., with 14.3% from the Pacific region, 27.8% from the South, 23.7% from the Midwest, 20.7% from the Northeast, and 13.5% from the Mountain region. The vast majority (83.1%) of dogs were considered companion animals. <xref rid="sec14" ref-type="sec">Supplementary Tables S1 and S2</xref> describe heartworm preventative use stratified by predictor variables.</p>
</sec>
<sec id="sec7"><label>3.2.</label>
<title>Elastic net model results</title>
<p>The final bootstrapped elastic net model is presented in <xref rid="tab2" ref-type="table">Table 2</xref>. Dogs from the Southern U.S. had higher odds of being on heartworm preventative [OR (95% CI): 1.34 (1.33&#x2013;1.35)] compared to other regions. Dogs from the Pacific region had lower odds [0.78 (0.77&#x2013;0.79)] compared to other regions. There were no differences in heartworm preventative use in males (39.5%) versus females (39.7%), but altered dogs had higher heartworm preventative use (43.7%) compared to intact dogs (37.0%). In the final model, binary variables for neutered and intact males were retained. Neutered male dogs had higher odds of being on heartworm preventative [1.19 (1.18&#x2013;1.20)] compared to the other sex categories, while intact male dogs had lower odds [0.859 (0.854&#x2013;0.865)]. Dogs in the highest height quartile (63.5&#x2009;cm and taller) had lower odds [0.742 (0.735&#x2013;0.750)] of being on heartworm preventative compared to shorter dogs.</p>
<table-wrap position="float" id="tab2"><label>Table 2</label>
<caption>
<p>Factors associated with heartworm preventative use in a sample of golden retrievers using a bootstrapped elastic net model (<italic>N</italic>&#x2009;=&#x2009;2,998 dogs).</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Variable<xref rid="tfn1" ref-type="table-fn"><sup>1</sup></xref></th>
<th align="center" valign="top"><italic>n</italic> (number in variable category)</th>
<th align="center" valign="top">Covariate Stability<xref rid="tfn2" ref-type="table-fn"><sup>2</sup></xref></th>
<th align="center" valign="top">Bootstrapped OR<xref rid="tfn3" ref-type="table-fn"><sup>3</sup></xref></th>
<th align="center" valign="top">Lower 90% bootstrapped CI</th>
<th align="center" valign="top">Upper 90% bootstrapped CI</th>
<th align="center" valign="top">Bootstrapped <italic>p</italic>-value<xref rid="tfn4" ref-type="table-fn"><sup>4</sup></xref></th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="bottom">Intercept</td>
<td align="center" valign="middle">&#x2013;</td>
<td align="center" valign="middle">1.00</td>
<td align="center" valign="middle">0.167</td>
<td align="center" valign="middle">0.162</td>
<td align="center" valign="middle">0.170</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Region (South)</td>
<td align="center" valign="middle">832</td>
<td align="center" valign="middle">0.99</td>
<td align="center" valign="middle">1.34</td>
<td align="center" valign="middle">1.33</td>
<td align="center" valign="middle">1.35</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Region (Pacific)</td>
<td align="center" valign="middle">428</td>
<td align="center" valign="middle">0.93</td>
<td align="center" valign="middle">0.78</td>
<td align="center" valign="middle">0.77</td>
<td align="center" valign="middle">0.79</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Sex (male neutered)</td>
<td align="center" valign="middle">498</td>
<td align="center" valign="middle">0.88</td>
<td align="center" valign="middle">1.19</td>
<td align="center" valign="middle">1.18</td>
<td align="center" valign="middle">1.20</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Sex (male intact)</td>
<td align="center" valign="middle">1,016</td>
<td align="center" valign="middle">0.88</td>
<td align="center" valign="middle">0.859</td>
<td align="center" valign="middle">0.854</td>
<td align="center" valign="middle">0.865</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Height (63.5&#x2009;cm and taller)</td>
<td align="center" valign="middle">326</td>
<td align="center" valign="middle">0.97</td>
<td align="center" valign="middle">0.742</td>
<td align="center" valign="middle">0.735</td>
<td align="center" valign="middle">0.750</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Infectious disease system diagnosis</td>
<td align="center" valign="middle">714</td>
<td align="center" valign="middle">0.84</td>
<td align="center" valign="middle">1.15</td>
<td align="center" valign="middle">1.14</td>
<td align="center" valign="middle">1.16</td>
<td align="center" valign="middle">0.0048</td>
</tr>
<tr>
<td align="left" valign="bottom">Ear/nose/throat system diagnosis</td>
<td align="center" valign="middle">594</td>
<td align="center" valign="middle">0.89</td>
<td align="center" valign="middle">1.17</td>
<td align="center" valign="middle">1.16</td>
<td align="center" valign="middle">1.18</td>
<td align="center" valign="middle">0.0067</td>
</tr>
<tr>
<td align="left" valign="bottom">Sun exposure location (concrete)</td>
<td align="center" valign="middle">442</td>
<td align="center" valign="middle">0.86</td>
<td align="center" valign="middle">1.19</td>
<td align="center" valign="middle">1.18</td>
<td align="center" valign="middle">1.20</td>
<td align="center" valign="middle">0.0047</td>
</tr>
<tr>
<td align="left" valign="bottom">Prior heartworm preventative use</td>
<td align="center" valign="middle">160</td>
<td align="center" valign="middle">0.96</td>
<td align="center" valign="middle">1.47</td>
<td align="center" valign="middle">1.45</td>
<td align="center" valign="middle">1.49</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Supplements given</td>
<td align="center" valign="middle">1,460</td>
<td align="center" valign="middle">0.87</td>
<td align="center" valign="middle">0.889</td>
<td align="center" valign="middle">0.884</td>
<td align="center" valign="middle">0.894</td>
<td align="center" valign="middle">0.0162</td>
</tr>
<tr>
<td align="left" valign="bottom">On tick preventative</td>
<td align="center" valign="middle">351</td>
<td align="center" valign="middle">1.00</td>
<td align="center" valign="middle">20.7</td>
<td align="center" valign="middle">20.4</td>
<td align="center" valign="middle">21.1</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom"># of rooms with carpeted flooring</td>
<td align="center" valign="middle">N/A<xref rid="tfn5" ref-type="table-fn"><sup>5</sup></xref></td>
<td align="center" valign="middle">0.86</td>
<td align="center" valign="middle">1.026</td>
<td align="center" valign="middle">1.025</td>
<td align="center" valign="middle">1.027</td>
<td align="center" valign="middle">0.0116</td>
</tr>
<tr>
<td align="left" valign="bottom">Spends time indoors in a room with hardwood floor</td>
<td align="center" valign="middle">1,620</td>
<td align="center" valign="middle">0.96</td>
<td align="center" valign="middle">1.20</td>
<td align="center" valign="middle">1.19</td>
<td align="center" valign="middle">1.21</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Received Bordetella vaccine</td>
<td align="center" valign="middle">1,439</td>
<td align="center" valign="middle">1.00</td>
<td align="center" valign="middle">1.40</td>
<td align="center" valign="middle">1.39</td>
<td align="center" valign="middle">1.41</td>
<td align="center" valign="middle">0.0000</td>
</tr>
<tr>
<td align="left" valign="bottom">Received rabies vaccine</td>
<td align="center" valign="middle">2,505</td>
<td align="center" valign="middle">0.83</td>
<td align="center" valign="middle">1.24</td>
<td align="center" valign="middle">1.23</td>
<td align="center" valign="middle">1.26</td>
<td align="center" valign="middle">0.0048</td>
</tr>
<tr>
<td align="left" valign="bottom">Received a different vaccine than rabies and Bordetella</td>
<td align="center" valign="middle">2,563</td>
<td align="center" valign="middle">0.99</td>
<td align="center" valign="middle">1.60</td>
<td align="center" valign="middle">1.58</td>
<td align="center" valign="middle">1.62</td>
<td align="center" valign="middle">0.0000</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="tfn1"><label>1</label>
<p>Variables were included in the final model if the bootstrapped <italic>p</italic>-value was&#x2009;&#x003C;&#x2009;0.02 and the covariate stability was&#x2009;&#x003E;&#x2009;80% (<xref rid="fig1" ref-type="fig">Figure 1</xref>).</p></fn>
<fn id="tfn2"><label>2</label>
<p>Covariate stability was defined as the percent of time the covariate was selected to be in the model in 500 bootstrapped samples.</p></fn>
<fn id="tfn3"><label>3</label>
<p>Bootstrapped ORs and 90% CI calculated from 500 bootstrapped samples. Each variable that is a level of a categorical variable is being compared to all other levels of that same categorical variable.</p></fn>
<fn id="tfn4"><label>4</label>
<p>Bootstrapped value of ps were calculated as one minus the percent of time the coefficient estimates were on the majority side of zero.</p></fn>
<fn id="tfn5"><label>5</label>
<p>Numerical predictor. Median 3, minimum 0, maximum 14.</p></fn>
</table-wrap-foot>
</table-wrap>
<p>Dogs with an infectious disease diagnosis or an ear/nose/throat system disease diagnosis had higher odds of being on heartworm preventative [1.15 (1.14&#x2013;1.16) and 1.17 (1.16&#x2013;1.18), respectively]. Those that had prior heartworm preventative use had almost 50% greater odds [1.47 (1.45&#x2013;1.49)] of currently being on heartworm preventatives. Interestingly, those that were on supplements had lower odds [0.889 (0.884&#x2013;0.894)] of being on heartworm preventative. The variable most strongly associated with heartworm preventative use was being on tick preventative. Those that were on tick preventative had 20.7 times greater odds (95% CI: 20.4&#x2013;21.1) of being on heartworm preventative compared to those not on tick preventative.</p>
<p>Being exposed to sun on concrete flooring was associated with higher heartworm preventative use [1.19 (1.18&#x2013;1.20)]. The odds of being on heartworm preventative increased with the number of carpeted rooms in the house [1.026 (1.025&#x2013;1.027)]. Spending time indoors on a room with hardwood flooring was associated with greater odds [1.20 (1.19&#x2013;1.21)] of being on heartworm preventative.</p>
<p>Dogs that received vaccines had higher odds of being on heartworm preventative, with the highest odds [1.60 (1.58&#x2013;1.62)] among those that received a different vaccine other than rabies and Bordetella (dogs in this category could also have gotten rabies and/or Bordetella vaccinations).</p>
</sec>
</sec>
<sec id="sec8" sec-type="discussions"><label>4.</label>
<title>Discussion</title>
<p>There are very few studies that have investigated what factors influence heartworm preventative use in dogs (<xref ref-type="bibr" rid="ref5">5</xref>, <xref ref-type="bibr" rid="ref6">6</xref>). To address this knowledge gap, we examined multiple factors associated with heartworm preventative use in a large sample of golden retrievers from the Golden Retriever Lifetime Study (<xref ref-type="bibr" rid="ref7">7</xref>). Due to the high number of predictors, we implemented a robust statistical analysis technique (bootstrapped elastic net regression) that can handle multicollinearity unlike traditional regression techniques. Regularization methods, also called penalized regression models, reduce risk of overfitting by constraining the estimated coefficients and decreasing out of sample error (<xref ref-type="bibr" rid="ref13">13</xref>). We identified multiple variables associated with heartworm preventative use, including vaccine use, U.S. region, sex, height, other disease conditions, medication and supplementation history, and characteristics of the spaces where the dog spends time.</p>
<p>Overall, the prevalence of heartworm preventative use was 39.5%. This was much lower than an estimate from a sample of dogs presenting at the veterinary teaching hospital at the University of Pennsylvania, which reported that 79.8% of dogs were currently on preventatives at the time of their admittance (<xref ref-type="bibr" rid="ref5">5</xref>). In our study, even among dogs located in the Southern U.S., which has the highest prevalence of heartworm preventative use and rate of heartworm disease, prevalence was only 45.9% (<xref ref-type="bibr" rid="ref2">2</xref>, <xref ref-type="bibr" rid="ref6">6</xref>). Correspondingly, use of heartworm preventatives was lowest in the Pacific region which has some of the lowest rates of heartworm disease (<xref ref-type="bibr" rid="ref2">2</xref>). Our estimates are similar to results reported in Bir et al., 2018 which found that 54% used heartworm prevention in a nationally representative sample. We also found that dogs in the highest quartile of height were less likely to be on heartworm preventative. Height is likely correlated with weight and cost could be prohibitive (<xref ref-type="bibr" rid="ref14">14</xref>). However, weight did not make it into the final model, so the mechanism is unclear. It could be due to the distribution of the data and the quartile ranges may not correspond directly to the weight ranges of heartworm medications. We also found that sex was associated with heartworm preventative use. The &#x201C;sex&#x201D; variable was a combination of biological sex and spay/neuter status. While there were no differences between males versus females, intact males and females had lower rates of heartworm preventative use compared to their altered counterparts. Spay/neuter status likely represents owner investment in health and veterinary care. This corroborates other research from Gates and Nolan (<xref ref-type="bibr" rid="ref5">5</xref>) who found that intact dogs were less likely to be on preventative compared to neutered dogs and that overall, there were no differences between males and females.</p>
<p>Dogs that received a Bordetella vaccine, rabies, or other vaccine had greater heartworm preventative use. This could be due to greater trust in veterinarian recommendations (<xref ref-type="bibr" rid="ref15">15</xref>). Interestingly, supplements were associated with lower heartworm preventative use. Pet owners that use supplements may prefer a holistic, naturopathic approach to pet care. They may use supplements such as black walnut or home remedies such as brewer&#x2019;s yeast or essential oils instead of heartworm preventatives (<xref ref-type="bibr" rid="ref15">15</xref>, <xref ref-type="bibr" rid="ref16">16</xref>). Prior use of heartworm preventatives and current use of tick preventatives were strongly associated with current heartworm preventative use which is not surprising since some heartworm preventatives also protect against infections from ticks. Offering a combination drug may increase adherence to heartworm preventatives, especially if owners do not understand and perceive the risk of heartworm disease (<xref ref-type="bibr" rid="ref17">17</xref>). Dogs that were diagnosed with an infectious disease system or an ear/nose/throat system disease were more likely to be on preventative. The infectious disease system includes gastrointestinal parasites; therefore, owners may be more likely to be on medications for parasite control after a diagnosis and many of these medications include heartworm control. The mechanism behind why an ear/nose/throat system diagnosis is associated with greater heartworm preventative use is unclear. One possible mechanism is that some of the symptoms of heartworm disease involve the ear/nose/throat system, such as coughing, so the dog may have been tested for heartworm disease and then subsequently put on preventative. In a survey among members of a hunting club, only 70% tested annually, which is expected to be lower in the general population (<xref ref-type="bibr" rid="ref17">17</xref>). Some of the dogs that do not get tested annually may only be tested when ill and presenting with symptoms that could be attributed to heartworm disease.</p>
<p>Dogs that were exposed to the sun on concrete flooring were more likely to be on preventative. The other sun exposure flooring categories were grass and dirt but were not retained in the final model during the variable selection process. Concrete flooring may be significant as it may indicate a closer proximity of the dog to the house, whereas dogs kept on grass and dirt may include dogs that free roam and that receive poorer general care compared to confined dogs (<xref ref-type="bibr" rid="ref18">18</xref>). Heartworm preventative use increased with the number of carpeted rooms in the house and was higher among dogs that spent time on hardwood floors. Both variables could be proxies for socioeconomic status. Owners of higher socioeconomic status may be more likely to spend money on heartworm preventative as they have more disposable income. However, Gates and Nolan (<xref ref-type="bibr" rid="ref5">5</xref>) did not find that income was significantly associated with heartworm preventative use in multivariable analyses, although Evason et al. (<xref ref-type="bibr" rid="ref15">15</xref>) reported that 33% of clients said cost was a barrier and downside to using heartworm preventatives. It is of note that variables that represent attachment and the human-animal bond were not significant in the final model. This is similar to results reported in Shore et al. (<xref ref-type="bibr" rid="ref19">19</xref>) who did not find a clear trend between attachment levels and heartworm preventative use.</p>
<p>A major limitation of this study is the limited generalizability outside the study population. The study population consisted of all golden retrievers that were less than 6&#x2009;years of old. The results of the study may not be generalizable to other breeds, older dogs, and to the general dog population. In addition, owners that participate in the GRLS may be inherently different than owners that do not participate. Owners that participate may be more invested in their dog therefore more likely to adopt veterinarian recommendations for heartworm prevention. However, this does not seem to be the case, as the prevalence of heartworm preventative use was low in our sample and was close to other national estimates (<xref ref-type="bibr" rid="ref6">6</xref>). Another limitation is that we did not have owner demographic data. We were only able to infer socioeconomic status from some of the variables. Owner factors including gender, age, and socioeconomic status have been shown to be associated with heartworm preventative use (<xref ref-type="bibr" rid="ref14">14</xref>). This study also had many strengths. We were able to leverage a large dataset representative of the U.S. that collected multiple variables in a standardized questionnaire. In addition, we were able to evaluate a number of factors that have not been tested for their relationship with heartworm preventative use prior to this study. Still, the findings of this study can be strengthened by conducting additional studies outside of the study population.</p>
<p>Overall, we found multiple factors associated with heartworm preventive use. Results of this study can be used to help veterinarians understand motivations and factors behind heartworm preventative use in dogs so that client communication can be improved. Results can also be used to identify populations at risk of low heartworm preventative use so that education and outreach efforts can be better targeted. Future studies can confirm the findings in a more diverse study population and use the results of this study to develop and test client communication and educational interventions.</p>
</sec>
<sec id="sec9" sec-type="data-availability">
<title>Data availability statement</title>
<p>The data analyzed in this study is subject to the following licenses/restrictions: Data can be requested through the Morris Animal Foundation Golden Retriever Lifetime Study RFP. Requests to access these datasets should be directed to <ext-link xlink:href="https://www.morrisanimalfoundation.org/golden-retriever-lifetime-study-rfp" ext-link-type="uri">https://www.morrisanimalfoundation.org/golden-retriever-lifetime-study-rfp</ext-link>.</p>
</sec>
<sec id="sec10">
<title>Ethics statement</title>
<p>The secondary data analysis conducted in this study was approved by the Lincoln Memorial University Institutional Review Board.</p>
</sec>
<sec id="sec11">
<title>Author contributions</title>
<p>LW was involved in developing the project proposal, study design, data analysis, and writing the first draft of the manuscript. CF and VF were involved with study design and editing the manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="sec12" sec-type="funding-information">
<title>Funding</title>
<p>This study was supported by the Golden Retriever Lifetime Study and this manuscript were made possible through financial support provided by the Morris Family Foundation, the Mark &#x0026; Bette Morris Family Foundation, VCA, the V Foundation, Blue Buffalo Company, Petco Love, Zoetis, Antech Inc., Elanco, the Purina Institute, Orvis, the Golden Retriever Foundation, the Hadley and Marion Stuart Foundation, Mars Veterinary, generous private donors and the Flint Animal Cancer Center at Colorado State University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<ack>
<p>We would like to thank Morris Animal Foundation, its staff members and all participants in the Golden Retriever Lifetime Study, including the dog owners, their golden retrievers and the Study veterinarians who made this work possible.</p>
</ack>
<sec id="sec14" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fvets.2023.1208804/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fvets.2023.1208804/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.xlsx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_2.xlsx" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><label>1.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Drake</surname><given-names>J</given-names></name> <name><surname>Wiseman</surname><given-names>S</given-names></name></person-group>. <article-title>Increasing incidence of Dirofilaria immitis in dogs in USA with focus on the southeast region 2013&#x2013;2016</article-title>. <source>Parasites Vectors</source>. (<year>2018</year>) <volume>11</volume>:<fpage>39</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13071-018-2631-0</pub-id>, PMID: <pub-id pub-id-type="pmid">29343304</pub-id></citation></ref>
<ref id="ref2"><label>2.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Self</surname><given-names>SW</given-names></name> <name><surname>Pulaski</surname><given-names>CN</given-names></name> <name><surname>McMahan</surname><given-names>CS</given-names></name> <name><surname>Brown</surname><given-names>DA</given-names></name> <name><surname>Yabsley</surname><given-names>MJ</given-names></name> <name><surname>Gettings</surname><given-names>JR</given-names></name></person-group>. <article-title>Regional and local temporal trends in the prevalence of canine heartworm infection in the contiguous United States: 2012&#x2013;2018</article-title>. <source>Parasites Vectors</source>. (<year>2019</year>) <volume>12</volume>:<fpage>380</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13071-019-3633-2</pub-id>, PMID: <pub-id pub-id-type="pmid">31362754</pub-id></citation></ref>
<ref id="ref3"><label>3.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Noack</surname><given-names>S</given-names></name> <name><surname>Harrington</surname><given-names>J</given-names></name> <name><surname>Carithers</surname><given-names>DS</given-names></name> <name><surname>Kaminsky</surname><given-names>R</given-names></name> <name><surname>Selzer</surname><given-names>PM</given-names></name></person-group>. <article-title>Heartworm disease &#x2013; overview, intervention, and industry perspective</article-title>. <source>Int J Parasitol Drugs Drug Resist</source>. (<year>2021</year>) <volume>16</volume>:<fpage>65</fpage>&#x2013;<lpage>89</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijpddr.2021.03.004</pub-id>, PMID: <pub-id pub-id-type="pmid">34030109</pub-id></citation></ref>
<ref id="ref4"><label>4.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mwacalimba</surname><given-names>K</given-names></name> <name><surname>Amodie</surname><given-names>D</given-names></name> <name><surname>Swisher</surname><given-names>L</given-names></name> <name><surname>Moldavchuk</surname><given-names>M</given-names></name> <name><surname>Brennan</surname><given-names>C</given-names></name> <name><surname>Walther</surname><given-names>C</given-names></name> <etal/></person-group>. <article-title>Pharmacoeconomic analysis of heartworm preventive compliance and revenue in veterinary practices in the United States</article-title>. <source>Front Vet Sci</source>. (<year>2021</year>) <volume>8</volume>:<fpage>602622</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fvets.2021.602622</pub-id></citation></ref>
<ref id="ref5"><label>5.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gates</surname><given-names>MC</given-names></name> <name><surname>Nolan</surname><given-names>TJ</given-names></name></person-group>. <article-title>Factors influencing heartworm, flea, and tick preventative use in patients presenting to a veterinary teaching hospital</article-title>. <source>Prev Vet Med</source>. (<year>2010</year>) <volume>93</volume>:<fpage>193</fpage>&#x2013;<lpage>200</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.prevetmed.2009.10.012</pub-id>, PMID: <pub-id pub-id-type="pmid">19931925</pub-id></citation></ref>
<ref id="ref6"><label>6.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bir</surname><given-names>C</given-names></name> <name><surname>Olynk Widmar</surname><given-names>N</given-names></name> <name><surname>Croney</surname><given-names>C</given-names></name></person-group>. <article-title>Exploring social desirability bias in perceptions of dog adoption: all&#x2019;s well that ends well? Or does the method of adoption matter?</article-title> <source>Animals</source>. (<year>2018</year>) <volume>8</volume>:<fpage>154</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ani8090154</pub-id>, PMID: <pub-id pub-id-type="pmid">30217039</pub-id></citation></ref>
<ref id="ref7"><label>7.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guy</surname><given-names>MK</given-names></name> <name><surname>Page</surname><given-names>RL</given-names></name> <name><surname>Jensen</surname><given-names>WA</given-names></name> <name><surname>Olson</surname><given-names>PN</given-names></name> <name><surname>Haworth</surname><given-names>JD</given-names></name> <name><surname>Searfoss</surname><given-names>EE</given-names></name> <etal/></person-group>. <article-title>The Golden retriever lifetime study: establishing an observational cohort study with translational relevance for human health</article-title>. <source>Phil Trans R Soc B</source>. (<year>2015</year>) <volume>370</volume>:<fpage>20140230</fpage>. doi: <pub-id pub-id-type="doi">10.1098/rstb.2014.0230</pub-id>, PMID: <pub-id pub-id-type="pmid">26056371</pub-id></citation></ref>
<ref id="ref8"><label>8.</label><citation citation-type="book"><person-group person-group-type="author"><name><surname>Kuhn</surname><given-names>M</given-names></name> <name><surname>Johnson</surname><given-names>K</given-names></name></person-group>. <source>Applied predictive modeling</source>. <publisher-loc>Germany</publisher-loc>: <publisher-name>Springer New York</publisher-name> (<year>2013</year>). p. <fpage>122</fpage>&#x2013;<lpage>128</lpage>.</citation></ref>
<ref id="ref9"><label>9.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bunea</surname><given-names>F</given-names></name> <name><surname>She</surname><given-names>Y</given-names></name> <name><surname>Ombao</surname><given-names>H</given-names></name> <name><surname>Gongvatana</surname><given-names>A</given-names></name> <name><surname>Devlin</surname><given-names>K</given-names></name> <name><surname>Cohen</surname><given-names>R</given-names></name></person-group>. <article-title>Penalized least squares regression methods and applications to neuroimaging</article-title>. <source>Neuroimage</source>. (<year>2011</year>) <volume>55</volume>:<fpage>1519</fpage>&#x2013;<lpage>27</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.neuroimage.2010.12.028</pub-id>, PMID: <pub-id pub-id-type="pmid">21167288</pub-id></citation></ref>
<ref id="ref10"><label>10.</label><citation citation-type="book"><person-group person-group-type="author"><name><surname>Hastie</surname><given-names>T</given-names></name> <name><surname>Tibshirani</surname><given-names>R</given-names></name> <name><surname>Mainwright</surname><given-names>M</given-names></name></person-group>. <source>Statistical learning with sparsity: the lasso and generalisations</source>. <publisher-name>Chapman and Hall/CRC</publisher-name> (<year>2015</year>) p. <fpage>142</fpage>&#x2013;<lpage>147</lpage>.</citation></ref>
<ref id="ref11"><label>11.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lima</surname><given-names>E</given-names></name> <name><surname>Green</surname><given-names>M</given-names></name> <name><surname>Lovatt</surname><given-names>F</given-names></name> <name><surname>Davies</surname><given-names>P</given-names></name> <name><surname>King</surname><given-names>L</given-names></name> <name><surname>Kaler</surname><given-names>J</given-names></name></person-group>. <article-title>Use of bootstrapped, regularised regression to identify factors associated with lamb-derived revenue on commercial sheep farms</article-title>. <source>Prev Vet Med</source>. (<year>2020</year>) <volume>174</volume>:<fpage>104851</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.prevetmed.2019.104851</pub-id>, PMID: <pub-id pub-id-type="pmid">31778947</pub-id></citation></ref>
<ref id="ref12"><label>12.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>McLaughlin</surname><given-names>D</given-names></name> <name><surname>Bradley</surname><given-names>A</given-names></name> <name><surname>Dottorini</surname><given-names>T</given-names></name> <name><surname>Giebel</surname><given-names>K</given-names></name> <name><surname>Leach</surname><given-names>K</given-names></name> <name><surname>Hyde</surname><given-names>R</given-names></name> <etal/></person-group>. <article-title>Identifying associations between management practices and antimicrobial resistances of sentinel bacteria recovered from bulk tank milk on dairy farms</article-title>. <source>Prev Vet Med</source>. (<year>2022</year>) <volume>204</volume>:<fpage>105666</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.prevetmed.2022.105666</pub-id>, PMID: <pub-id pub-id-type="pmid">35594608</pub-id></citation></ref>
<ref id="ref13"><label>13.</label><citation citation-type="book"><person-group person-group-type="author"><name><surname>Boehmke</surname><given-names>B</given-names></name> <name><surname>Greenwell</surname><given-names>BM</given-names></name></person-group>. <source>Hands-on machine learning with R</source>. <publisher-loc>United States</publisher-loc>: <publisher-name>CRC Press</publisher-name> (<year>2019</year>). p. <fpage>122</fpage>&#x2013;<lpage>126</lpage>.</citation></ref>
<ref id="ref14"><label>14.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bir</surname><given-names>C</given-names></name> <name><surname>Ortez</surname><given-names>M</given-names></name> <name><surname>Olynk Widmar</surname><given-names>NJ</given-names></name> <name><surname>Wolf</surname><given-names>CA</given-names></name> <name><surname>Hansen</surname><given-names>C</given-names></name> <name><surname>Ouedraogo</surname><given-names>FB</given-names></name></person-group>. <article-title>Familiarity and use of veterinary services by US resident dog and cat owners</article-title>. <source>Animals</source>. (<year>2020</year>) <volume>10</volume>:<fpage>483</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ani10030483</pub-id>, PMID: <pub-id pub-id-type="pmid">32183120</pub-id></citation></ref>
<ref id="ref15"><label>15.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Evason</surname><given-names>M</given-names></name> <name><surname>McGrath</surname><given-names>M</given-names></name> <name><surname>Stull</surname><given-names>J</given-names></name></person-group>. <article-title>Companion animal preventive care at a veterinary teaching hospital &#x2013; knowledge, attitudes, and practices of clients</article-title>. <source>Can Vet J</source>. (<year>2021</year>) <volume>62</volume>:<fpage>484</fpage>&#x2013;<lpage>90</lpage>. PMID: <pub-id pub-id-type="pmid">33967287</pub-id></citation></ref>
<ref id="ref16"><label>16.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Means</surname><given-names>C</given-names></name></person-group>. <article-title>Selected herbal hazards</article-title>. <source>Vet Clin Small Anim</source>. (<year>2002</year>) <volume>32</volume>:<fpage>367</fpage>&#x2013;<lpage>82</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0195-5616(01)00009-2</pub-id></citation></ref>
<ref id="ref17"><label>17.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rohrbach</surname><given-names>BW</given-names></name> <name><surname>Lutzy</surname><given-names>A</given-names></name> <name><surname>Patton</surname><given-names>S</given-names></name></person-group>. <article-title>Attributes, knowledge, beliefs, and behaviors relating to prevention of heartworm in dogs among members of a national hunting dog club</article-title>. <source>Vet Parasitol</source>. (<year>2011</year>) <volume>176</volume>:<fpage>324</fpage>&#x2013;<lpage>32</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.vetpar.2011.01.017</pub-id>, PMID: <pub-id pub-id-type="pmid">21310533</pub-id></citation></ref>
<ref id="ref18"><label>18.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Astorga</surname><given-names>F</given-names></name> <name><surname>Poo-Mu&#x00F1;oz</surname><given-names>DA</given-names></name> <name><surname>Organ</surname><given-names>J</given-names></name> <name><surname>Medina-Vogel</surname><given-names>G</given-names></name></person-group>. <article-title>Why let the dogs out? Exploring variables associated with dog confinement and general characteristics of the free-ranging owned-dog population in a peri-urban area</article-title>. <source>J Appl Anim Welf Sci</source>. (<year>2022</year>) <volume>25</volume>:<fpage>311</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1080/10888705.2020.1820334</pub-id></citation></ref>
<ref id="ref19"><label>19.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shore</surname><given-names>ER</given-names></name> <name><surname>Douglas</surname><given-names>DK</given-names></name> <name><surname>Riley</surname><given-names>ML</given-names></name></person-group>. <article-title>What&#x2019;s in it for the companion animal? Pet attachment and college students&#x2019; behaviors towards pets</article-title>. <source>J Appl Anim Welf Sci</source>. (<year>2005</year>) <volume>8</volume>:<fpage>1</fpage>&#x2013;<lpage>11</lpage>. doi: <pub-id pub-id-type="doi">10.1207/s15327604jaws0801_1</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn0003"><p><sup>1</sup><ext-link xlink:href="https://www.morrisanimalfoundation.org/sites/default/files/filesync/GRLS-Annual-Veterinarian-Questionnaire.pdf" ext-link-type="uri">https://www.morrisanimalfoundation.org/sites/default/files/filesync/GRLS-Annual-Veterinarian-Questionnaire.pdf</ext-link></p></fn>
<fn id="fn0004"><p><sup>2</sup><ext-link xlink:href="https://www.morrisanimalfoundation.org/sites/default/files/filesync/GRLS-Annual-Owner-Questionnaire.pdf" ext-link-type="uri">https://www.morrisanimalfoundation.org/sites/default/files/filesync/GRLS-Annual-Owner-Questionnaire.pdf</ext-link></p></fn>
<fn id="fn0005"><p><sup>3</sup><ext-link xlink:href="https://www.heartwormsociety.org/images/pdf/2018-AHS-Canine-Guidelines.pdf" ext-link-type="uri">https://www.heartwormsociety.org/images/pdf/2018-AHS-Canine-Guidelines.pdf</ext-link></p></fn>
</fn-group>
</back>
</article>