{'fraction_non_alphanumeric': 0.032134563484591726, 'fraction_numerical': 0.01575346764576665, 'mean_word_length': 4.6381811748053785, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mycoplasma hyopneumoniae Surveillance in Pig Populations: Establishing Sampling Guidelines for Detection in Growing Pigs

Antemortem detection of Mycoplasma hyopneumoniae infection in swine production systems has relied on antibody testing, but the availability of tests based on DNA detection and novel diagnostic specimens, e.g., tracheal swabs and oral fluids, has the potential to improve M. hyopneumoniae surveillance. A field study was performed over a 14-week period during which 10 pigs in one pen at the center of a room with 1,250 6-week-old pigs housed in 46 pens were intratracheally inoculated with M. hyopneumoniae.

primarily through nose-to-nose contact, such that one infected pig can infect 1.16 pigs during a 6-week nursery period (3). This low rate of transmission is one of the most challenging aspects of M. hyopneumoniae disease management and surveillance, especially for early detection in naive populations.
M. hyopneumoniae surveillance may be done using various sampling strategies and testing protocols. Selection of specimen and test depends on the level of disease in the herd, accuracy of the test, and cost. Most commonly, antemortem monitoring of M. hyopneumoniae is done using enzyme-linked immunosorbent assay (ELISA) because of the ease of sample collection, lower cost, and diagnostic performance (4). However, the utility of ELISAs is compromised by the highly variable period between infection and antibody production (;3 to 8 weeks), the ambiguous relationship between positive results and clinical disease, and the inability to differentiate natural infection from vaccination (3)(4)(5).
Alternatively, PCR may be used, but the process of sampling must provide the highest likelihood of the presence of M. hyopneumoniae DNA in the sample. M. hyopneumoniae establishes itself within the lower respiratory tract, mainly the trachea and bronchi, and thus sampling needs to target these sites (6). Several specimens have been evaluated for the detection of M. hyopneumoniae by PCR, including nasal, tonsil, laryngeal, and tracheal samples. A recent study carried out by Sponheim et al. (7) showed a significant increase in diagnostic sensitivity of tracheal samples compared to that of laryngeal swabs. However, collection of tracheal samples requires more time and animal restraint than serum collection. Furthermore, for an individual animal sample, the cost per PCR is roughly six times that of an ELISA.
In contrast to sera and tracheal samples, aggregate specimens, such as oral fluids, require less labor, are less stressful for the pigs, and have demonstrated greater sensitivity versus individual pig samples for the detection of other swine diseases, such as porcine reproductive and respiratory syndrome virus (PRRSV) (8). Previous research suggested that oral fluids could serve as a suitable diagnostic specimen for M. hyopneumoniae, especially during active clinical infection (9).
Surveillance protocols for introduction of negative replacement gilts into naive herds have traditionally relied on the testing of serum samples to validate M. hyopneumoniae-negative herd status. This approach provides M. hyopneumoniae surveillance with minimal labor and expense, but the probability of M. hyopneumoniae detection using this approach is uncertain. Likewise, estimates of the probability of detection for alternative sampling approaches, e.g., tracheal samples and oral fluids, are not available. Therefore, the objective of this study was to estimate the probability of M. hyopneumoniae detection in tracheal samples (DNA), oral fluids (DNA), and sera (antibodies) as a function of M. hyopneumoniae prevalence and sample size.

MATERIALS AND METHODS
Animal care and housing. This study was reviewed and approved by Iowa State University (ISU)'s Institutional Animal Care and Use Committee (IACUC-18-141). The study was conducted on a wean-tofinish site with two connected double-wide barns (2 rooms per barn) housing ;5,000 pigs. Barns had automatic tunnel ventilation with liquid propane gas brooders for weaned pigs and water misters for heat dispersion. Manure was collected in deep pits, and the site was managed all-in-all-out. The study was carried out in one room (1,250 6-week-old pigs) with 23 pens (;28 pigs per pen) on either side of a central alleyway. Two pens were kept empty, were half the size of regular pens, and were strictly used as recovery or hospital pens by the production system. The room was stocked with 21-day-old barrows confirmed negative for M. hyopneumoniae, PRRSV, and influenza A virus (IAV), utilizing molecular (VetMax PRRSV NA & EU, VetMax, and Gold SIV detection kits; Thermo Fisher Scientific) and serological assays (IDEXX M. hyo Ab test, IDEXX PRRS OF Ab ELISA test, and IDEXX influenza Ab test). During the study, pigs did not receive M. hyopneumoniae vaccine and were not treated with M. hyopneumoniae-susceptible medication. All animal veterinary care, housing, handling, and feeding were under the supervision of production system veterinarians.
M. hyopneumoniae inoculation and sampling of seeder pigs (inoculated pen). Seven days prior to the initiation of the experiment, one tracheal, one serum, and one oral fluid sample were collected from each pen from randomly selected pigs to confirm the M. hyopneumoniae-negative status of the room.
To initiate M. hyopneumoniae infection in the room, 10 pigs conveniently selected from a centrally located pen were ear-tagged and intratracheally inoculated with an M. hyopneumoniae lung homogenate (1 Â 10 5 CCU/ml of strain 232, provided by Iowa State University, Ames, IA). Briefly, pigs were immobilized using a snare, and the airway was visualized with a mouth speculum and laryngoscope. A total of 10 ml of lung homogenate was deposited into the trachea via the use of a feeding catheter (10). To confirm M. hyopneumoniae infection, serum and tracheal samples were collected on a weekly basis from every pig in the inoculated pen until each pig was confirmed positive by M. hyopneumoniae ELISA and M. hyopneumoniae PCR, respectively. In addition, one pen-based oral fluid sample was collected weekly from the inoculated pen until the sample was confirmed positive by M. hyopneumoniae PCR.
Sampling of uninoculated pens. Beginning at 14 days postinoculation (DPI) and continuing every 14 days over a 14-week period, one tracheal sample, four sera, and one oral fluid sample were collected from each uninoculated pen in the room (n = 45). Tracheal samples were collected by restraining the pig with a snare and mouth speculum and introducing a single-use catheter into the trachea, as previously described (11). Serum samples were collected via venipuncture from the jugular vein. Oral fluid samples were collected by suspending a rope from the pen gate such that the end of the rope was level with the pigs' shoulders. The rope was suspended in the pen for 30 min (12). During each sampling event, samples were collected by first starting at the southeast section of the barn and then moving north until the last pen from that row was sampled. At this point, samplers then moved to the southwest section, moving north until the last pen from that row was sampled. All samplers were blind to test pen results over the course of the study. The first randomly selected pig in each pen was selected for both tracheal and serum samples. This protocol was executed for a total of seven sampling events over a 14-week period.
Diagnostic testing. All tracheal samples, sera, and oral fluids were submitted to the Iowa State University Veterinary Diagnostic Laboratory (VDL) for testing by either PCR or ELISA. M. hyopneumoniae DNA from tracheal samples and oral fluids was extracted using the MagMax-96 pathogen RNA/DNA kit (Applied Biosystems, Carlsbad, CA). Extracted DNA was amplified using a PCR TaqMan Fast virus 1-step mastermix (Life Technologies, Carlsbad, CA) with a previously described protocol (13). PCR runs were performed on an ABI Prism 7500 machine (Applied Biosystems, Waltham, MA). A sample with cycle threshold (C T ) values lower than 37 were considered positive; otherwise, samples were considered negative.
M. hyopneumoniae-specific antibodies in all sera were measured utilizing a commercially available enzyme-linked immunosorbent assay (ELISA) test (IDEXX M. hyo Ab test). All samples were tested using the same ELISA equipment, namely, plate washer (ELx405; Biotek Instruments, Inc., Winooski, VT), ELISA reader (EMax Plus microplate reader; Molecular Devices, San Jose, CA), and reader software (SoftMax Pro 7.0; Molecular Devices). Following manufacturer's instructions, antibody concentration was expressed as the ratio of optical densities (OD) from the sample and mean positive (sample to positive [S/P] ratio). Samples with an S/P ratio equal to or greater than 0.3 were considered positive; otherwise, samples were considered negative.
A total of 26 samples from randomly selected PCR-positive tracheal samples were submitted to the ISU VDL for p146 M. hyopneumoniae sequencing to confirm circulation of the challenge strain (M. hyopneumoniae 232) (14).
Respiratory distress index (SOMO devices). Five SOMO devices (SoundTalks NV, Precision Livestock Farming, Belgium) were placed equidistant from each other in the alleyway of the room. SOMO 1 was placed at the far north end of the room. A second device, SOMO 2, was placed north of the center of the room. SOMO 3 was placed in the center of the room. SOMO 4 was placed south of the center of the room. SOMO 5 was placed at the south end of the room. These devices continuously recorded sound throughout the trial. Data from these devices provided measurements for temperature, humidity, and a respiratory distress index (RDI). RDI was determined through a proprietary algorithm. The algorithm automatically detected and classified individual coughs. The algorithm differentiated cough from other sounds such as grunts, squeals, and background noise such as fans, gates, and feeders. A threshold was incorporated into the algorithm to signal an RDI alert. In addition, researchers recorded the first observation of coughing during sampling events.
Statistical analysis. A hierarchical Bayesian model based on a latent spatial piecewise exponential model (15) was constructed to estimate the spread of M. hyopneumoniae and the diagnostic sensitivities of tracheal and oral fluid PCRs and serum ELISA. Modeling and graphs were performed in R (R program version 3.6.0, package rjags 4.8; R core team 2019).
Let i i ¼ 1; 2; . . . ; 46 ð Þ denote the i th pen in the barn and n i denote the number of pigs in the i th pen. The sampling time points were 0 ¼ t 0 ,t 1 ,t 2 , . . . ,t K ,1, and the time unit was 2 weeks. Let Þdenote the diagnostic test methods, 1 for tracheal PCR test, 2 for serum ELISA, and 3 for oral fluid PCR test. Let u ijkl be the test outcome of pig j in pen i at the k th sampling time using test method l. Note that the oral fluid PCR test was a pen-level method, where the term u i*k3 represented the oral fluid PCR test result in the pen i at the k th sampling time. To any of the tests, the outcome was binned to 0 or 1 (negative or positive result). Let y ijk denote the true infection status of pig j in pen i at the k th sampling time; y ijk = 1 if it was infected and 0 otherwise. Once a pig was considered infected, it remained infected until the end of the study. If y ijk0 ¼ 1, then y ijk0 ¼ 1y ijk1 ¼ 1 for all k 1 $ k 0 . Let t ij denote the unknown time to event of pig j in pen i, so that y ijk ¼ I t ij , t k ð Þ , where I(·) is an indicator function that takes a value of 1 when the condition in the parentheses is true and 0 otherwise. The distance between pens was defined to estimate the transmission of M. hyopneumoniae. The barn comprised two rows of pens. Let d ii9 be the distance between pen i and i9. If pens i and i9 were adjacent, then d ii9 5 1; otherwise, d ii9 was calculated by Euclidean distance between pens i and i9. Note that d ii 5 0. Let g 1l and g 0l denote the diagnostic sensitivity and specificity of the l th test and p ik be the within-pen prevalence of the i th pen at the k th sampling time. In this model, it was assumed that g 0l ¼ 1 for estimating g 1l ; l ¼ 1; 2; 3: (i) Model to represent misclassification. Given the unknown true infection status, the distribution of diagnostic test results was estimated as follows.
For l 5 1 (tracheal PCR test), let To account for delayed detection of M. hyopneumoniae antibodies that are typically detected around 21 to 56 days after exposure (16), a time lag of 21 days (1.5 time units) was implemented; hence, for l 5 2 (serum ELISA), let A lag parameter, q, was introduced to account for the delayed detection in oral fluids. Hence, for l 5 3 (oral fluid PCR test), let For each sampling event, four serum ELISA and one tracheal sample PCR result were obtained for each uninoculated pen. The tracheal sample and the first serum sample were collected from the same pig. A new variable was introduced, namely z ikm ; discrete uniform 1; Modified spatial piecewise exponential model. The Cox proportional hazards model (17) for the time to event is as follows: where l 0 (t) is the baseline hazard function, x is a vector of explanatory variables, and bis a vector of fixed effect parameters. In the study, the duration of the study was partitioned by (K 1 1) sampling time points, 0 ¼ t 0 ,t 1 ,t 2 , . . . ,t K ,1. Let the k th interval be t k21 ; t k ð , and l ijk be the hazard of pig j in pen i being infected in the interval t k21 ; t k ð , given that pig was not infected by time t k21 . The conditional distribution of t ij jy ij;k21 ¼ 0 followed an exponential distribution.
Assume that the pens' prevalence at time t k21 can affect the hazard of pens at time t k . To account for the prevalence effect, the prevalence and the spatial distance were introduced as the covariates into the following hazard function: where b 0 is the regression coefficient for the baseline hazard function, b 1 is a regression coefficient about the effect of pen prevalence but not associated with the spatial distance, and b 2 is a coefficient considering both the spatial distance and pen prevalence effect. After transformation, the hazard function for the interval was t k21 ; t k ð : Parameter C was introduced to control whether a pen was infected or not. Let s ijk be a threshold parameter where It was assumed a pig in a pen can be infected only if the pen was infected. Pen i was infected at time point t k if s iÃk , C. Then the conditional probability of pig j in pen i being infected at t k ðk $ 1) is The model parameters g 1,1 , g 1,2 , g 1,3 , b 0 , b 1 , b 2 , and C were estimated through a hierarchical Bayesian model. The empirical priors of g 1,1 , g 1,2 , and g 1,3 were calculated using a beta-binomial model. The priors for other parameters were noninformative.

RESULTS
Confirmation of inoculated pen infection. All tracheal samples, oral fluids and serum samples collected before inoculation tested negative for M. hyopneumoniae by PCR and ELISA. All seeder pigs were positive by tracheal sample PCR on DPI 7 and by serum ELISA on DPI 42. All seeder pig pen mates were positive by tracheal sample PCR on DPI 21 and by serum ELISA on DPI 70. Pen-level oral fluids were positive by PCR on DPI 56. Two of the seeder pigs died during the course of the study, and their lungs were harvested and submitted to the ISU VDL for analysis. PCR testing, histopathology, and p146 sequencing (14) Table 1). The parameter estimates, standard errors, and 95% credible intervals are shown in Table 1. The probability of M. hyopneumoniae detection in tracheal samples, oral fluids, and serum samples over time was simulated, at each iteration of which the estimated parameters b 0 , b 1 , b 2 , and C were used to simulate the spread of infection and calculate the within-barn prevalence at sampling point t k ; k 2 f0; 1; . . . ; Kg. The simulation was initiated with one randomly selected infected pen containing one infected pig. Sampling time (units of 7 days) and sample sizes were specified. The following sampling rules were implemented in the model: A. If the number of samples was less than the number of pens, then each pen had one sample at most.
B. If the number of samples was more than the number of pens, then each pen had at least one sample for each diagnostic test method.
C. The same pig cannot be sampled twice for the same diagnostic test, i.e., for the ELISA; if the number of samples was 50 (exceeding the number of pens), then some pens may have two samples and the two samples are taken from different pigs in that pen.
D. Sampling was memoryless. For example, if a pig was sampled at time t k in a simulation iteration, then that pig could also be sampled after t k in that simulation iteration.
E. The same pig could be sampled twice for different types of diagnostic test methods at the same time. For example, if a pig was sampled for the serology test method at time t k , then this pig could also be sampled for a tracheal PCR test at time t k .
Given a sample size, samples using the fixed spatial sampling strategy at each sampling time t k for each diagnostic test method were selected. Fixed spatial sampling was defined as allocating samples equidistantly across the sampling space. The probability of detection was estimated as the proportion of detection out of 10,000 simulations. This same approach was used to simulate probability of detection over time if five pigs were infected in one pen on DPI 0. These estimates are shown in Table 2 and Fig. 2.
To estimate the point in time probability of detection for tracheal samples and sera, the number of positive pens and the within-pen prevalence was defined. To carry out simulated sampling for the defined number of positive pens and within-pen prevalence, the following rules were applied: A. If the number of samples was less than the number of pens, then each pen had one sample at most.
B. If the number of samples was more than the number of pens, then each pen had at least one sample for each diagnostic test method.
C. The same pig could not be sampled twice for the same diagnostic test, i.e., for the ELISA; if the number of samples was 50 (exceeding the number of pens), then some pens may have two samples and the two samples are taken from different pigs in that pen.
For each iteration, given the number of positive pens and the within-pen prevalence, the positive pens and the infected pigs were simulated. For example, if there were five positive pens and the within-pen prevalence was 0.2, then five pens were randomly selected to be positive and the infection status of the pigs in those pens followed a Bernoulli distribution (P 5 0.2). After simulating positive pens and infected pigs, samples were selected based on the sample size. Pens were selected using fixed spatial sampling, while pigs within pens were selected using random sampling. The probability of detection was calculated as the proportion of simulations (out of 10,000) with $1 positive pig among the total pigs sampled. Probability of detection estimates are shown in Tables 2 to 4.
Respiratory distress index. Figure 3 showed the RDI values recorded during the study for each of the five SOMO devices. The first RDI alert was recorded on DPI 55 by SOMO 3, located in the center of the room nearest to the inoculated pen. On DPI 78, SOMO 1, located at the north end of the barn recorded an RDI alert. Both SOMO 1 and 2, set up on the north end of the room, recorded two RDI alerts on DPI 86. In total, 8 RDI alerts were recorded throughout the study.

DISCUSSION
Effective surveillance must achieve timely detection of the disease of interest in the target population. A single standardized approach for effective surveillance of all pathogens does not exist; rather, each surveillance program must adopt a pathogenspecific approach that optimizes detection in a cost-effective manner. Due to the elusiveness of M. hyopneumoniae and the limited sensitivity of diagnostic tools, designing an effective M. hyopneumoniae surveillance program is particularly challenging. The foundation of any surveillance program is the estimation of the probability of detection associated with a given combination of sample type, diagnostic test, sample size, sample allocation, and prevalence (8,18). This is a prerequisite for understanding the effectiveness of a herd surveillance program. Therefore, the objective of this study was to estimate the probability of M. hyopneumoniae detection in tracheal samples, oral fluids, and sera as a function of prevalence, time since initial introduction, and sample size. Overall, the timing of M. hyopneumoniae detection varied among sample types. Tracheal samples provided the most consistent and earliest detection throughout the study, followed by serum and oral fluids. The M. hyopneumoniae pattern of detection in tracheal samples and sera followed the direction of the airflow, beginning north of the inoculated pen and spreading throughout the south end of the room before spreading to the north end of the room. This observed movement of M. hyopneumoniae may suggest an important role of aerosol transmission within an airspace and its potential implications for early detection. A previous study demonstrated airborne spread of M. hyopneumoniae as a critical risk factor for interfarm transmission (19). However, fomites or pig handlers may have also played a role in the transmission of M. hyopneumoniae in this study.
While the first detection of M. hyopneumoniae was obtained by serum ELISA on DPI 14, M. hyopneumoniae was detected by PCR in tracheal samples in more pens at DPI 42, 56, and 70 compared to detection in serum. These results were expected given the documented lag between infection and production of antibodies (16). Compared to tracheal samples and serum, M. hyopneumoniae detection in oral fluids was delayed, more variable, and less sensitive until DPI 84. Previous research suggested that detection of M. hyopneumoniae in oral fluids was associated with clinical signs, e.g., cough, which may serve to increase the antigen in the oral cavity (9). In this study, cough activity was monitored using SOMO devices. SOMO 3, located in the center of the room, recorded the first RDI alert on DPI 55. The highest average RDI across all SOMO devices was recorded on DPI 95, which coincided with the highest detection of positive pens by oral fluid samples, an observation consistent with the hypothesis proposed by Hernandez-Garcia et al. (9) (Fig. 3).  The field data were modeled to estimate the probability of detection for each sample type, which had not been previously evaluated. M. hyopneumoniae infection was solely determined by tracheal sample and oral fluid PCR and serum ELISA results; the model also provided diagnostic sensitivity estimates for each sample type. In this study, the diagnostic sensitivity estimated for tracheal samples was estimated at 0.965. These results are in accordance with those of a recent study (7). Serum followed, with a sensitivity of 0.818. Recently, diagnostic sensitivity for the IDEXX ELISA was estimated to be 0.56 (20). The noted differences between sensitivities are likely due to the differences in definition of false positives between studies. In this study, a time lag parameter of 21 days was included to account for delayed detection of M. hyopneumoniae antibodies (16). Oral fluids provided the lowest diagnostic sensitivity at 0.396 (Table 1). This number was less than half the diagnostic sensitivity of tracheal and serum samples, supporting the use of these sample types over oral fluid considering current available diagnostic tools.
Two different scenarios (1 versus 5 initially infected animals) were constructed to evaluate the effect of time, sample size, and sample type on the probability of M. hyopneumoniae detection. The number of initially infected animals was chosen to mimic  a As a function of the within-pen prevalence (5% or 50%), number of pens sampled using a fixed spatial approach (n), and number of positive pens in the barn. For serum results, the model assumed that the agent had been in the population at least 21 days. b TS, tracheal sample; SS, serum sample. scenarios of very low initial prevalence (0.08% versus 0.4%), and thus the resulting estimates provided in the study were conservative. In scenario 1, 120 tracheal samples collected on DPI 7 provided an ;50% probability of detection. Furthermore, 120 tracheal samples would provide a probability of detection of ;96% on DPI 35. In contrast, probability of detection estimates for serum and oral fluid samples could not be generated until DPI 21 and DPI 35, respectively. The earliest detection for serum and oral fluid samples at ;50% probability of detection occurred with sample sizes of 90 and 60, respectively, on DPI 35. To achieve higher probabilities of detection using serum samples or oral fluids, sample sizes would need to be increased, but time to detection would still be delayed. For example, 120 serum samples would need to be collected at DPI 56 to have a 95% probability of detection. Likewise, 120 oral fluids collected on DPI 49 provided a 95% probability of detection. Interestingly, at DPI 35, 30 oral fluid samples provided almost double the probability of detection (22% versus 39%) as the same number of serum samples (Fig. 2). Reasonably, at an increased initial prevalence of 0.4% (scenario 2), the required sample size and time of detection decreases across sample types; however, tracheal samples remain the optimal sample type. For example, 60 tracheal samples collected on DPI 14 would provide a 57% probability of detection. To achieve 95% probability of detection, 90 tracheal samples on DPI 35 or 120 on DPI 28, respectively, would need to be collected ( Table 2). One of the highest-risk events for a swine production system is the introduction of replacement animals into a herd (21). In most production systems, these animals will undergo a quarantine period (;30 days) in a separate facility prior to being introduced into the population. During quarantine, it is imperative that diagnostics are performed to determine the health status of the replacement animals. The results of this study support the collection of a high number of tracheal samples for systems that implement a #30-day quarantine period (Table 2). In the most conservative scenario, a traditional protocol consisting of 30 serum samples at 28 days post quarantine start would yield a probability of detection of 13% at a cost of $165 (assuming $5.50/ELISA). Given the high risk of introducing potentially positive animals, this traditional approach does not provide a high level of confidence. While doubling the sample size and changing the sample type (i.e., tracheal samples) increases the likelihood of detection by 50%, it also increases the cost 10-fold ($1,800; assuming $30/PCR). For example, collecting $60 tracheal samples at 28 days increases the probability of detection to 75%. Despite the significant increase, this sampling approach (i.e., 60 tracheal samples) still does not reach the industry's traditional expectation of 95% probability of detection (22). Thus, producers should consider extending the quarantine period to 60 days to achieve probabilities of detection closer to 100% using more convenient and economical sample types (i. e., serum and oral fluids). Selection of sample type and size depends on the level of risk, the impact of the disease, the possibility for extended duration of quarantine, and the cost the operation is willing to pay. For instance, if a production site is experiencing clinical signs suggestive of M. hyopneumoniae or had a recent biosecurity breach, the herd veterinarian will need to decide between collecting and testing a large set of tracheal samples immediately or collecting a smaller set of tracheal, serum, or oral fluid samples at a later time point (Table 2). In the immediate response scenario, the cost of collecting and testing a large set of tracheal samples is high; however, the probability of detection also increases. In the delayed-response scenario, the cost associated with labor and testing is reduced, but this response implies that the risk of delaying diagnosis is minimal.
Probability of detection estimates are necessary for implementation of routine surveillance programs to demonstrate freedom from disease in negative populations. In these cases, veterinarians select a sample size given the desired level of confidence and the cost of sample collection and testing. Testing all animals in a population would provide high confidence in a herd's negative status, but this approach is impractical. Traditionally, sample sizes have been selected based on estimates provided by Cannon and Roe (22) and Cannon (23), which are based on the approximation of a hypergeometric distribution. The results of this study indicated that 95% probability of detection would require 90 tracheal samples to be collected using a fixed spatial sampling approach from a population of 1,250 in which 50% of the animals across four pens were M. hyopneumoniae positive (Table 3). In contrast, using Cannon's (23) approximation under the same assumptions would require 65 tracheal samples to be collected using a simple random sampling approach. The underestimation of sample size by Cannon (23) highlights the challenge found in the approximation of the hypergeometric distribution when applied to swine disease surveillance. The differences between the sample sizes are due to the use of a simple hypergeometric model versus the spatial piecewise exponential model. The spatial piecewise exponential model accounts for spatial distance and its relationship with transmission, and is supported by field data. Hypergeometric distribution assumes that the variable of interest, in this case, disease, is randomly distributed with less variation than a spatial model. Furthermore, as demonstrated in previous work, estimates found by Cannon (23) cannot be applied to oral fluid samples, as these are an aggregate sample (18).
Modeling the field data supported previously documented trends related to probability of detection (18,22). At larger sample size and higher disease prevalence, probability of detection increases for all three sample types. Tracheal samples provided the highest diagnostic sensitivity, earliest detection, and highest probability of detection estimates, but at a significantly higher cost. While the cost associated with increased labor and collection materials incurred in tracheal sampling cannot be reduced, pooling may be an effective approach to decreasing the cost of testing. Future studies should evaluate the effect of pooling schemes on probability of detection. While diagnostic sensitivity of sera was higher than oral fluids, after DPI 35 (scenario 1), oral fluid samples provided a higher probability of detection at every sampling point for the same sample sizes. It is important to point out that increasing the sample size of serum samples at the same time point will achieve similar probability of detection, at a lower testing cost. For example, 60 serum samples, would provide the same 39% probability of detection at a third of the cost of 30 oral fluid samples. While testing penbased oral fluid samples allows for an increased number of pigs represented in one sample and can be tested for multiple pathogens, future studies should explore potential improvements in M. hyopneumoniae diagnostic performance of this sample type. Probability of detection estimates for all M. hyopneumoniae-specific sample types have been generated by this study and these will be critical in the refinement of current and future M. hyopneumoniae surveillance programs.