Comments on Fouchier’s Calculation of Risk and Elapsed Time for Escape of a Laboratory-Acquired Infection from His Laboratory

摘要

LETTER In a Letter to the Editor of mBio , Professor Ron Fouchier published a calculation ( 1 ) in which he finds a very low probability, P _(1), for a laboratory-acquired infection (LAI) for a single lab for a single year. Claiming numerous safety precautions in his biosafety level 3 (BSL3+) laboratory, Fouchier calculates P _(1)= 1 × 10~(–7)per person per year, and since there are 10 workers with access to his laboratory, P _(1)= 1 × 10~(–6)per lab per year. Compare this to P _(1)= 2 × 10~(–3)per lab per year for BSL3 laboratories calculated from CDC statistics for undetected or unreported LAIs ( 2 , 3 ), here called “community LAIs,” as it is assumed that an undetected or unreported LAI represents an infection that has traveled outside the lab and into the community. Recently reported escapes of LAIs from high-level biocontainment at CDC laboratories ( 4 ) and the long history of LAIs and other escapes from laboratories ( 5 ) also argue that Fouchier’s value for P _(1)is too low. Lipsitch and Inglesby ( 6 ) have supplied additional arguments as to why the Fouchier value for P _(1)is likely much too low. Fouchier uses a simplistic formula, y = 1/ P _(1), to calculate the elapsed time in years for an LAI to escape from his laboratory, y = 1/(1 × 10~(–6)) = 1 × 10~(6), that is, the million years stated in his Letter. It is not clear what this calculation tells us. Does it give us the elapsed time for a 10% chance that an LAI occurs? Does it give us elapsed time for a 50% chance, or an 80% chance? In this regard, the elapsed time for a 100% chance is infinite, as we can never be absolutely certain that an LAI will occur. I suggest attaching little weight to this elapsed time calculation and instead concentrating on risk = likelihood × consequences, starting with the P _(1)probability, specifically: potential pandemic fatalities = (probability of a community LAI) × (probability that the community LAI leads to a pandemic) × (estimated fatalities in a pandemic). My risk calculation estimates the likelihood of a community LAI for both a single laboratory and n laboratories conducting this research over y years. The total number of laboratories involved in this potential pandemic pathogen research is called here the “research enterprise.” A single, easily derived equation is used to determine the likelihood of a community LAI: (1) E = 1 ? ( 1 ? P 1 ) y n where E is the probability of at least one community LAI from n laboratories in y years. Example results are presented in Table?1 for three values of P _(1). TABLE?1? Probabilities of at least one LAI escape into the community ~(c) P _(1) Comment Estimated probability E n = 1 ~(a) n = 30 ~(b) 2.00E–03 From CDC data 0.0198 0.45 2.00E–04 10× less than CDC data 0.0020 0.058 1.00E–06 From Fouchier’s analysis 0.00001 0.0003 ~(a) A single BSL3 or BSL3+ lab. ~(b) Twice the number of NIH-funded labs studying gain-of-function pathogens with pandemic potential. ~(c) In equation 1, y is the number of years that research was carried out. In Table?1 , the number of laboratories is either n = 1 for a single laboratory, such as Fouchier’s, or n = 30, which is twice the 15 laboratories currently subject to the NIH funding pause. Picking n = 30 is a reasonable guess since there are likely many other labs throughout the world conducting this research that are not funded by NIH. y = 10?years is a reasonable time frame for this research to be completed. The rationale for picking the probabilities, P _(1), in Table?1 is as follows: P _(1)= 2 × 10~(–3)is calculated from the CDC statistics ( 2 , 3 ). P _(1)= 2 × 10~(–4)is 10-fold less and is my “guestimate” for a BSL3+ lab with rigorous safety practices. P _(1)= 1 × 10~(–6)is Fouchier’s calculated value. There are valuable observations to be gleaned from Table?1 . For instance, taking into account the whole research enterprise, not just a single lab, is important. Furthermore, even for Fouchier’s very low value for P _(1), there is an estimated probability of E = 0.0003 that there will be at least one community LAI over a 10-year period for 30 labs (likely exactly one LAI, as two is much less probable). As I will soon show, E = 0.0003, or 0.03%, is not nearly small enough to reduce risk to an acceptable level. Also, this E ?value assumes that all 30 laboratories involved in this research enterprise have the rigorous safety practices of Fouchier’s lab, a highly unlikely assumption. Summarizing the literature, Lipsitch and Inglesby ( 7 ) estimate the probability that a community LAI leads to a global spread (pandemic) to be 5 to 60%. This range is consistent with the 5 to 15% range found by Merler and coworkers ( 8 ) and with the 1 to 30% range found in a focused risk assessment ( 9 ) for infection spread beginning on crowded public transportation. As an illustration, using an intermediate value of 10% for pandemic probability, which is within the estimated ranges, the probability that a community LAI occurs and leads to a pandemic woul