Assessment of methods for prediction of human West Nile virus (WNV) disease from WNV-infected dead birds
© Veksler et al; licensee BioMed Central Ltd. 2009
Received: 03 June 2008
Accepted: 05 June 2009
Published: 05 June 2009
West Nile virus (WNV) is currently the leading cause of arboviral-associated encephalitis in the U.S., and can lead to long-term neurologic sequelae. Improvements in dead bird specimen processing time, including the availability of rapid field laboratory tests, allows reassessment of the effectiveness of using WNV-positive birds in forecasting human WNV disease.
Using New York State integrated WNV surveillance data from transmissions seasons in 2001–2003, this study determined which factors associated with WNV-positive dead birds are most closely associated with human disease. The study also addressed the 'delay' period between the distribution of the dead bird variable and the distribution of the human cases. In the last step, the study assessed the relative risk of contracting WNV disease for people who lived in counties with a 'signal' value of the predictor variable versus people who lived in counties with no 'signal' value of the predictor variable.
The variable based on WNV-positive dead birds [(Positive/Tested)*(Population/Area)] was identified as the optimum variable for predicting WNV human disease at a county level. The delay period between distribution of the variable and human cases was determined to be approximately two weeks. For all 3 years combined, the risk of becoming a WNV case for people who lived in 'exposed' counties (those with levels of the positive dead bird variable above the signal value) was about 2 times higher than the risk for people who lived in 'unexposed' counties, but risk varied by year.
This analysis develops a new variable based on WNV-positive dead birds, [(Positive/Tested)*(Population/Area)] to be assessed in future real-time studies for forecasting the number of human cases in a county. A delay period of approximately two weeks between increases in this variable and the human case onset was identified. Several threshold 'signal' values were assessed and found effective at indicating human case risk, although specific thresholds are likely to vary by region and surveillance system differences.
West Nile virus (WNV) was first recognized in the Northeast United States in 1999 . Since 1999 the virus has spread across the country, resulting in 28,943 human cases and 1130 deaths through 2008 (reported as of February 13, 2009) . WNV is now found throughout the western hemisphere . Infection can lead to long-term neurologic sequelae in people , and is currently the leading cause of arboviral-type encephalitis in the U.S. .
The number of human cases reported to the Centers for Disease Control and Prevention (CDC) currently varies widely in the U.S., from a low of 1 case in North Carolina, South Carolina and West Virginia in 2006 to a high of 996 cases in Idaho . Some of this variation can be attributed to different human population sizes and surveillance systems for infection, with some states more aggressive at testing and reporting milder, non-neuroinvasive disease. Although most WNV exposures do not result in clinical disease and most disease is mild, with case fatality rates ranging from 3% to 15%  disease prevention should be a priority. Personal protective measures and mosquito control measures all have a resource cost, and thus it is important to try and determine when the risk of human infection is high or low so that individuals and governmental agencies can make appropriate decisions about prevention.
Dead bird surveillance has served as one key method for tracking WNV activity in the U.S. [8–10]. Previous studies have documented that dead crow sightings can serve as a valuable index for forecasting human cases before or without laboratory confirmation of WNV infection, because crows have had a case-fatality rate close to 100% and are reasonably easy for the public to recognize and report [11, 12]. Multiple laboratory diagnostic methods may be used for serologic and virologic diagnosis of WNV, although collection, submission, processing, and testing can be lengthy [13, 14]. Thus, dead crow sightings, which can be used for surveillance immediately upon reporting, provide a more immediate indicator of WNV activity than WNV-positive birds. However, more rapid processing methods are in use in many areas, including VecTest and RAMP, that can be used to provide a laboratory result from swabs taken where the bird is found [13, 15]. Development of dead bird indicators using only WNV-positive birds eliminates the possibility of misclassification inherent in a dead crow sighting index, and it allows use of other species, which is particularly important in areas with few crows or other corvids (blue jays, ravens), or where WNV has reduced the crow population. Previous studies have indicated that WNV is sometimes first detected in an area in a bird species other than a corvid .
This study is an exploratory study focusing on determining which factors, using WNV-positive and tested dead birds, are most closely associated with the number of human WNV cases at a county level. The study also addresses the issue of "delay" between the distribution of a predictor variable based on laboratory-tested dead birds and the distribution of human disease. Previous studies have noted that the time between mosquito bites and human disease onset (incubation period) is 2–14 days . This study expands upon prior analyses by considering daily distributions of the predictor variable and human disease and applying statistical methods to assess the "delay" period between these variables. Finally, this study assesses the predictive value for the risk of human disease of using a weekly county level predictor variable based on WNV-positive and tested dead birds with several signal levels.
To identify measures of WNV activity in birds that might provide an indication of increased risk in humans, analyses focused on NYS integrated WNV surveillance data collected during the transmission season, defined as a 16-week period using weeks as defined for reporting by the CDC as the 26th week to the 41st week for the years 2001–2003 (June 24–Oct. 13 in 2001, June 23–Oct. 12 in 2002, and June 22–Oct. 11 in 2003). This time frame was optimal for inclusion of all human disease in New York State.
The NYS integrated surveillance system includes real-time surveillance components for humans, mammals, birds, and mosquitoes. The surveillance system relies on the public to report sightings of dead birds to local health departments (LHDs). Some of the reported dead birds were collected and tissue specimens were tested for WNV infection at the NYSDOH Wadsworth Center's Arthropod-borne Disease Laboratory, according to laboratory protocols previously described .
For WNV surveillance in humans, healthcare providers were asked to report patients with encephalitis and aseptic meningitis to the LHDs. Thus, the NYSDOH surveillance system emphasizes reporting of human cases with neuroinvasive disease, although cases without neuroinvasive disease are not excluded. Patients were tested for WNV infection at the NYSDOH Wadsworth Center's Diagnostic Immunology Laboratory (serology), Arbovirus Laboratory (PRNT), and Viral Encephalitis Laboratory (PCR), according to laboratory protocols previously described .
According to the CDC case definitions, New York (excluding New York City) had six (one excluded from study) confirmed or probable human WNV cases disease in 2001, 52 (three excluded) in 2002 and 40 (three excluded) in 2003. Two counties had human WNV cases in 2001, 12 counties had human WNV cases in 2002, and ten counties reported WNV human cases in 2003. Data from 57 NYS counties were included in this study. New York City data were excluded from analysis because New York City developed its own WNV surveillance system to monitor dead bird reports, with different priorities for reporting.
These analyses were conducted with county as the unit of analysis for several reasons. The small number of human cases made further subdivision challenging. In addition, the actual location of infection for humans and birds is unknown, so county can serve as a reasonable surrogate, acknowledging the potential for misclassification and reduced statistical power if infection actually occurred in a different county. Finally, these analyses were conducted to aid in prevention and control decision-making, which occurs at the county level.
Variables measuring WNV activity in dead birds
Variables constructed empirically from positive bird and tested bird surveillance variables
Surveillance variables standardized by population of county, area of county, or population density a
Proportion of positive birds and proportion of positive birds standardized by population or population density
A1 = [Positive/Area]
A6 = [Positive/Tested]
A2 = [Tested/Area]
A7 = [(Positive/Tested)*Population]
A3 = [Positive/Population]
A8 = [(Positive/Tested)*(Population/Area)]
A4 = [Positive*(Population/Area)]
A5 = [Tested*(Population/Area)]
To determine the significant variables among the eight variables considered as potential predictor variables using WNV-positive and/or tested dead birds, a correlation matrix between each of the variables for every county with at least one human case and the human cases in each county across all three years was generated using SAS software. The variables were considered as highly correlated with human disease if p ≤ 0.05. Multiple regression was then used to determine the best model predicting the number of human disease cases using the adjusted R2 statistic (adjusted for the number of parameters in the model) .
Previous studies have noted that WNV-positive dead birds are usually found before the onset of the human infection [21, 22]. Studies in the eastern U.S. found that the delay period between the first WNV-positive bird and the first human case varied from 15 days to 92 days . Possible explanations for the delay include the amplification cycle between birds and mosquitoes and the transmission cycle from mosquitoes to humans through mosquito bites. The time between the exposure (mosquito bite) and onset of human disease (incubation period) is known to be two to 14 days . To determine the delay period for our study, the distribution of the variable with the highest correlation with human disease for every day of the transmission period across all counties and for all three years pooled together was considered.
To obtain the best estimate of the delay period between the occurrence of WNV-positive birds and onset of human disease, we calculated the mean value of the optimal predictor dead bird variable for every day of the transmission period across all counties with human cases for all three years. The study period each year was 16 weeks (113 days). For each of the 113 days for every county with human cases we found the values of the predictor variable identified in the previous modeling for the particular day, then summarized them across all counties and divided the sum by 24, the number of counties with human cases in 2001–2003. If no birds were tested on that day, then the predictor variable was considered as zero.
To determine the delay period, two methods, non-parametric and parametric, were considered. The non-parametric two-sample Kolmogorov-Smirnov test was used to test whether the two underlying probability distributions of the predictor variable and human disease differed significantly . This test was performed using SAS software with the proc npar1way, edf option. By shifting the values of our variable against the daily distribution of human cases, the best shift in days that maximized goodness of fit between two distributions was determined. The second (parametric) method was applied to smoothed daily distributions of our predictor variable and human cases. The goal of this method was to maximize correlation between two distributions over different shifts and select the shift that provides the maximum. This test was performed using SAS proc corr for every value of shift.
Estimation of WNV disease risk
The Cochran-Mantel-Haenszel (CMH) test was used to calculate point and interval estimates for the relative risk of becoming a human WNV case depending on the value of the predictor variable in a person's county of residence. For each week, a table was constructed comparing the week's human cases per population in counties with a "signal" value of the predictor variable (greater than a threshold pre-determined for county use based on preliminary estimates) and the human cases per population in counties with no signal. Analyses were conducted with a weekly distribution of the predictor variable, to avoid small numbers or zeroes in daily distributions. The CMH chi-square statistic was used to compare the incidence (risk) of WNV disease in signal areas with the incidence in the non-signal areas over all the weeks of study. The CMH test was performed for several threshold values of the predictor variable using SAS Proc freq to obtain the estimate of odds ratio and relative risk [SAS System for Windows V8, SAS Institute, Cary, NC, USA].
The county was included in the CMH analysis if it tested at least 10 birds during the transmission period. In 2001, only 16 counties tested at least 10 birds. In 2002, 53 counties tested at least 10 birds. In 2003, 48 counties tested at least 10 birds. All counties with human cases, except one, were retained in the analysis because they all tested at least 10 birds during the transmission period. The exception was Schuyler County with a human case in 2003 but only five birds tested. This county was excluded from analysis in 2003.
Sample CMH data table for assessing relative risk of becoming a human WNV case
Week 38, 2003
Threshold = defined levela
In counties with predictor variableb ≥ threshold at defined time periodc before case onset
In counties with predictor variable < threshold at defined time period before case onset
Number of persons with disease onset
Number of persons without disease onset
Association of WNV-positive dead birds with human WNV disease
Correlation matrix, number of human cases with positive and tested dead bird variablesa
Pearson Correlation Coefficients, N = 24 (Counties) Prob > |r| under H0: Rho = 0
Multiple regression models for prediction of the number of human WNV cases in a county
Using three variables most highly correlated with number of human cases
A6, A7, A8
Using three variables with lowest intercorrelations
A3, Pos, A8
We also examined the smoothed plot of our two distributions (Figure 1c). We found correlations between human cases at time t and our predictor variable at time t+p. We determined the optimal shift p when the correlation reached its maximum. Although the correlation reached its maximum at 14 days, values for shifts 12, 13, 14 and 15 days were very close to each other, so the optimal shift can be 12–15 days. Figure 1d represents the plot of values of the correlation coefficient for the predictor variable against shifts in days.
Estimation of WNV disease risk
Association of weekly predictor variablea values and number of human cases, by signal value
2.79 (0.57, 13.53)
1.43 (0.79, 2.58)
4.03 (1.99, 8.12)
2.21 (1.43, 3.41)
1.54 (0.79, 2.97)
2.34 (1.51, 3.62)
4.87 (2.38, 9.90)
2.43 (1.57, 3.77)
1.55 (0.84, 2.84)
2.45 (1.58, 3.45)
1.36 (0.76, 2.42)
4.73 (2.37, 9.46)
2.22 (1.45, 3.40)
Previous studies have assessed the value of other WNV variables for forecasting potential increases in human disease. In one study that used dead bird surveillance information (WNV positive birds and tested birds) to construct the predictor variables, the study focused on the early transmission season, identified as a six-week period in June-July . Another study utilized dead crow sightings, defining the appropriate period for forecasting as up to two weeks before human case onset based on incubation period, but did not conduct an analysis of delay period . Many studies that tracked WNV-positive birds reported that they were found before onset of human disease. Our study is the first statistical comparison of the daily distributions of the dead bird predictor variable and the number of human cases to confirm a delay period between these two distributions. Using either raw or smoothed data, the delay period was about two weeks, which is consistent with the previous conclusions based on incubation period.
The estimates of relative risk will be affected by the choice of predictor variable. For example, most of the NY human WNV cases in 2001–2003 were reported from counties with a high human population density. Counties with low population density usually had only single WNV cases (Clinton, Orleans, and Yates counties in 2002; Warren, Cattaraugus, Yates, and Dutchess counties in 2003). A notable exception is Broome County in 2002, with seven human WNV cases. Although the overall county human population density in Broome County is relatively low (283.64 persons per square mile), the cases were clustered around the higher population density city of Binghamton. Because the analyses for this study were done at the county level, Broome County was classified as not 'exposed' for the dead bird predictor variable thresholds we considered. This may be one factor affecting the low relative risk estimates for 2002.
To address the issue of the large variation in population density among counties, future studies may benefit from combining data from several states and grouping counties with similar population densities. Different threshold values may be applicable for each group. Another approach would be to assign weights to each county according to its population density, to account for the large variation among counties. The weights can be determined, for example, as a ratio of the population densities of the most populated county to the given county. When the dead bird predictor variable is multiplied by this weight, the counties will become more comparable.
In interpreting these findings, several limitations of the study methodology must be considered. Reporting of dead birds was a passive system, depending solely on the public. Variability in reporting interest between counties could lead to misclassification on the predictor variable, reducing the chance of detecting an association. Not all reported dead birds were collected and tested, which could also affect estimates of the predictor variable. Only reported human cases were used in the analysis, and thus human infection (either asymptomatic or mild, and thus unreported) was likely rarely recognized and reported. This would also reduce the statistical power of detecting an association. The number of NY human cases was relatively small in each year compared to some other states, especially in 2001, which was one reason for conducting our analyses at the county level. With a larger number of cases, analyses could be done at the sub-county level allowing more precise interpretation of risk for vector control decisions. Infection in other counties due to movement of humans and birds will lead to misclassification and reduction of statistical power. Although the relative risks were elevated using the threshold values in the study for all three years, they were only statistically significant for 2003 and all three years combined. More years of data, or combining data across states, would increase the power of our analysis.
Estimation of WNV risk using the CMH test has limitations based on the parameters chosen for this study. This study restricted county entry into the study by having at least 10 birds tested during the transmission season, and the effect of using other numbers of birds for determining inclusion was not examined. Because the predictor variable is [(Positive/Tested)*(Population/Area)], higher levels of birds tested per week could be considered in order to justify that a county has sufficient surveillance to warrant inclusion in the study. Different threshold values could also be considered, and may be appropriate in areas with different WNV disease dynamics and different surveillance systems. The predictor variable consists of two parts: proportion of positive birds and human population density. The proportion can take values between zero and one; therefore the maximum value of our predictor variable is the population density of the particular county. Counties with small population densities may never be able to reach specific threshold values to provide a 'signal' of risk. On the other hand, their small human population size indicates an inherent lower risk of those counties having a human case, in comparison to counties with higher human populations even with the same proportion of infected mosquitoes and birds. Thus, regional rather than county analyses may be more useful in more sparsely populated areas, particularly if such analyses can compare urban and rural areas. However, if using the analyses to determine specific prevention and control activities based on having a "signal" of risk, analyses in smaller geographic units can be more helpful in determining where to target those interventions. Finally, these "signals" were generated based on bird surveillance alone, and it is possible that "signals" developed that also incorporated mosquito surveillance data might offer even better prediction of human case risk, but mosquito surveillance data is often not widely available with rapid test results for real-time analyses.
In New York State, decisions about WNV prevention and control are primarily made at the county level. These decisions have resource consequences, and use of chemicals for mosquito control can raise concerns. Thus, decision-making is usually based on a wide variety of inputs to help determine level of risk and need for control, including more focal information such as dead bird clusters or areas with high mosquito infection rates. Previous studies have indicated the value of dead bird indicators [7–16, 21, 24], and this study provides a novel statistical approach to demonstrate the value of using WNV-infected dead birds as initial indicator of WNV disease risk for a county. Signals values of our weekly predictor variable above all of the threshold values evaluated in this study were associated with increased risk of WNV disease, with significant associations in one year and across the three years of the study. The North American WNV strain appears to be a phenotype highly virulent to American crows . However, dead bird indicators using American crows may be of less value in the future if there are reductions in the case fatality rate or reduced numbers of birds for surveillance due to WNV die-offs.
The authors' work on this publication was partially supported through the Health Research Institute and Cooperative Agreement Number U90/CCU216988 from the Centers for Disease Control and Prevention (CDC), and National Institute of Allergy and Infectious Diseases-NIH contract NO1-AI-25490 to Laura Kramer. The publication's contents are solely the responsibility of the authors and do not necessarily represent the official views of CDC or NIH. The authors thank the following staff who were instrumental in developing the data used in this study: the NYSDOH Center for Community Health's Zoonoses Program for bird surveillance data and the Arthropod-borne Disease Program for human surveillance data; the NYSDOH Wadsworth Center's Diagnostic Immunology Laboratory and Viral Encephalitis Laboratory for human sample laboratory testing and the Arbovirus Laboratory for bird testing; the NYS Department of Environmental Conservation's Wildlife Pathology Unit for dead bird processing; the United States Department of Agriculture's Wildlife Services New York Office for maintaining the dead bird reporting hotline; and the local health departments, health care providers, and others submitting specimens and reports to the NYS WNV surveillance system.
- Lanciotti RS, Roehrig JT, Deubel V, Smith J, Parker M, Steele K, Crise B, Volpe KE, Crabtree MB, Scherret JH, Hall RA, MacKenzie JS, Cropp CB, Panigraphy B, Ostlund E, Schmitt B, Malkinson M, Banet C, Weissman J, Komar N, Savage HM, Stone W, McNamara T, Gubler DJ: Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States. Science. 1999, 286 (5448): 2333-7. 10.1126/science.286.5448.2333View ArticlePubMedGoogle Scholar
- CDC (Centers for Disease Control and Prevention). West Nile Virus: Statistics, Surveillance and Control. [http://www.cdc.gov/ncidod/dvbid/westnile/surv%26;control.htm]
- Kramer LD, Li J, Shi PY: West Nile virus. Lancet Neurol. 2007, 6: 171-181. 10.1016/S1474-4422(07)70030-3View ArticlePubMedGoogle Scholar
- Sejvar JJ: The Long-Term Outcomes of Human West Nile Virus Infection. Clin Infect Dis. 2007, 44: 1617-24. 10.1086/518281View ArticlePubMedGoogle Scholar
- Bitto A, Hakim JA, Pula MF: What Seniors Say About West Nile Virus – The Threat is Not Over. Californian J Health Promot. 2005, 3 (3): 127-36.Google Scholar
- CDC (Centers for Disease Control and Prevention). 2006 Human Case Counts. West Nile Virus Statistics, Surveillance and Control. 2005, [http://www.cdc.gov/ncidod/dvbid/westnile/surv%26;controlCaseCount06_detailed.htm]
- CDC (Centers for Disease Control and Prevention). Cases of West Nile Encephalitis: Questions and Answers. [http://www.cdc.gov/ncidod/dvbid/westnile/qa/cases.htm]
- Eidson M, Komar N, Sorhage F, Nelson R, Talbot T, Mostashari F, McLean R, : Crow deaths as a sentinel surveillance system for West Nile virus in the northeastern United States, 1999. Emerg Infect Dis. 2001, 7 (4): 615-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Blackmore CG, Stark LM, Jeter WC, Oliver RL, Brooks RG, Conti LA, Wiersma ST: Surveillance results from the first West Nile virus transmission season in Florida, 2001. Am J Trop Med Hyg. 2003, 69 (2): 141-50.PubMedGoogle Scholar
- Watson JT, Jones RC, Gibbs K, Paul W: Dead crow reports and location of human West Nile virus cases, Chicago, 2002. Emerg Infect Dis. 2004, 10 (5): 938-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Eidson M: "Neon needles" in a haystack: the advantages of passive surveillance for West Nile virus. Ann N Y Acad Sci. 2001, 951: 38-53.View ArticlePubMedGoogle Scholar
- Eidson M, Schmit K, Hagiwara Y, Anand M, Backenson PB, Gotham I, Kramer L: Dead crow density and West Nile virus monitoring, New York. Emerg Infect Dis. 2005, 11 (9): 1370-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Stone WB, Okoniewski JC, Therrien JE, Kramer LD, Kauffman EB, Eidson M: VecTest as diagnostic and surveillance tool for West Nile virus in dead birds. Emerg Infect Dis. 2004, 10 (12): 2175-81.PubMed CentralView ArticlePubMedGoogle Scholar
- Trevejo RT, Eidson M: West Nile virus. J Am Vet Med Assoc. 2008, 232 (9): 1302-1309. 10.2460/javma.232.9.1302View ArticlePubMedGoogle Scholar
- Stone WB, Therrien JE, Benson R, Kramer L, Kauffman EB, Eidson M, Campbell S: Assays to detect West Nile virus in dead birds. Emerg Infect Dis. 2005, 11 (11): 1770-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Eidson M, Kramer L, Stone W, Hagiwara Y, Schmit K, : Dead bird surveillance as an early warning system for West Nile virus. Emerg Infect Dis. 2001, 7 (4): 615-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen LR, Marfin AA: West Nile virus: a primer for the clinician. Ann Intern Med. 2002, 137: 173-9.View ArticlePubMedGoogle Scholar
- Kauffman EB, Jones SA, Dupuis AP, Ngo KA, Bernard KA, Kramer LD: Virus detection protocols for West Nile virus in vertebrate and mosquito specimens. J Clin Microbiol. 2003, 41 (8): 3661-3667. 10.1128/JCM.41.8.3661-3667.2003PubMed CentralView ArticlePubMedGoogle Scholar
- Wong SJ, Demarest VL, Boyle RH, Wang T, Ledizet M, Kar K, Kramer LD, Fikrig E, Koski RA: Detection of human anti-flavivirus antibodies with a West Nile virus recombinant antigen microsphere immunoassay. J Clin Microbiol. 2004, 42 (1): 65-72. 10.1128/JCM.42.1.65-72.2004PubMed CentralView ArticlePubMedGoogle Scholar
- Der G, Brian SE: A Handbook of Statistical Analyses using SAS Second edition. Boca Raton, FL: Chapman and Hall/CRC; 2002.Google Scholar
- Julian KG, Eidson M, Kipp AM, Weiss E, Peterson LR, Miller JR, Hinten SR, Marfin A: Early season crow mortality as a sentinel for West Nile virus disease in humans, northeastern United States. Vector Borne Zoonotic Dis. 2002, 2: 145-55. 10.1089/15303660260613710View ArticlePubMedGoogle Scholar
- Marfin AA, Peterson LR, Eidson M, Miller J, Hadler J, Farello C, Werner B, Campbell GL, Layton M, Smith P, Bresnitz E, Cartter M, Scaletta J, Obiri G, Bunning M, Craven RC, Roehrig JT, Julian KG, Hinten SR, Gubler DJ, : Widespread West Nile Virus Activity, Eastern United States, 2000. Emerg Infect Dis. 2001, 7 (4): 730-35.PubMed CentralView ArticlePubMedGoogle Scholar
- SAS online documentation. Tests based on the Empirical Distribution Function. [http://support.sas.com/documentation/onlinedoc/91pdf/sasdoc_91/stat_ug_7313.pdf]
- Eidson M, Miller J, Kramer L, Cherry B, Hagiwara Y, : Dead Crow densities and human cases of West Nile Virus, New York State, 2002. Emerg Infect Dis. 2001, 10: 938-40.Google Scholar
- Brault AC, Huang CY-H, Langevin SA, Kinney RM, Bown RA, Ramey WN, Panella NA, Holmes ED, Powers AM, Miller BR: A single positively selected West Nile viral mutation confers incrased virogenesis in American crows. Nat Genet. 2007, 39 (9): 1162-6. 10.1038/ng2097PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ([http://creativecommons.org/licenses/by/2.0]), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.