Basic reproduction number

In epidemiology, the basic reproduction number, or basic reproductive number (sometimes called basic reproduction ratio or basic reproductive rate), denoted (pronounced R nought or R zero),[1] of an infection is the expected number of cases directly generated by one case in a population where all individuals are susceptible to infection.[2] The definition assumes that no other individuals are infected or immunized (naturally or through vaccination). Some definitions, such as that of the Australian Department of Health, add the absence of "any deliberate intervention in disease transmission".[3] The basic reproduction number is not necessarily the same as the effective reproduction number (usually written [t for time], sometimes ),[4] which is the number of cases generated in the current state of a population, which does not have to be the uninfected state. is a dimensionless number (persons infected per person infecting) and not a time rate, which would have units of time−1,[5] or units of time like doubling time.[6]

is the average number of people infected from one other person. For example, Ebola has an of two, so on average, a person who has Ebola will pass it on to two other people.

is not a biological constant for a pathogen as it is also affected by other factors such as environmental conditions and the behaviour of the infected population. values are usually estimated from mathematical models, and the estimated values are dependent on the model used and values of other parameters. Thus values given in the literature only make sense in the given context and it is not recommended to compare values based on different models.[7] does not by itself give an estimate of how fast an infection spreads in the population.

The most important uses of are determining if an emerging infectious disease can spread in a population and determining what proportion of the population should be immunized through vaccination to eradicate a disease. In commonly used infection models, when the infection will be able to start spreading in a population, but not if . Generally, the larger the value of , the harder it is to control the epidemic. For simple models, the proportion of the population that needs to be effectively immunized (meaning not susceptible to infection) to prevent sustained spread of the infection has to be larger than .[8] This is the so-called Herd immunity threshold or herd immunity level. Here, herd immunity means that the disease cannot spread in the population because each infected person, on average, can only transmit the infection to less than one other contact.[9] Conversely, the proportion of the population that remains susceptible to infection in the endemic equilibrium is . However, this threshold is based on simple models that assume a fully mixed population with no structured relations between the individuals. For example, if there is some correlation between people's immunization (e.g., vaccination) status, then the formula may underestimate the herd immunity threshold.[9]

Graph of herd immunity threshold vs basic reproduction number with selected diseases

The basic reproduction number is affected by several factors, including the duration of infectivity of affected people, the contagiousness of the microorganism, and the number of susceptible people in the population that the infected people contact.[10]

History

edit

The roots of the basic reproduction concept can be traced through the work of Ronald Ross, Alfred Lotka and others,[11] but its first modern application in epidemiology was by George Macdonald in 1952,[12] who constructed population models of the spread of malaria. In his work he called the quantity basic reproduction rate and denoted it by  .

Overview of estimation methods

edit

Compartmental models

edit

Compartmental models are a general modeling technique often applied to the mathematical modeling of infectious diseases. In these models, population members are assigned to 'compartments' with labels – for example, S, I, or R, (Susceptible, Infectious, or Recovered). These models can be used to estimate  .

Epidemic models on networks

edit

Epidemics can be modeled as diseases spreading over networks of contact and disease transmission between people.[13] Nodes in these networks represent individuals and links (edges) between nodes represent the contact or disease transmission between them. If such a network is a locally tree-like network, then the basic reproduction can be written in terms of the average excess degree of the transmission network such that:

 

where   is the mean-degree (average degree) of the network and   is the second moment of the transmission network degree distribution.

Heterogeneous populations

edit

In populations that are not homogeneous, the definition of   is more subtle. The definition must account for the fact that a typical infected individual may not be an average individual. As an extreme example, consider a population in which a small portion of the individuals mix fully with one another while the remaining individuals are all isolated. A disease may be able to spread in the fully mixed portion even though a randomly selected individual would lead to fewer than one secondary case. This is because the typical infected individual is in the fully mixed portion and thus is able to successfully cause infections. In general, if the individuals infected early in an epidemic are on average either more likely or less likely to transmit the infection than individuals infected late in the epidemic, then the computation of   must account for this difference. An appropriate definition for   in this case is "the expected number of secondary cases produced, in a completely susceptible population, produced by a typical infected individual".[14]

The basic reproduction number can be computed as a ratio of known rates over time: if a contagious individual contacts   other people per unit time, if all of those people are assumed to contract the disease, and if the disease has a mean infectious period of  , then the basic reproduction number is just  . Some diseases have multiple possible latency periods, in which case the reproduction number for the disease overall is the sum of the reproduction number for each transition time into the disease.

Effective reproduction number

edit
An explanation of the   number in simple terms from the Welsh Government.

In reality, varying proportions of the population are immune to any given disease at any given time. To account for this, the effective reproduction number   or   is used.   is the average number of new infections caused by a single infected individual at time t in the partially susceptible population. It can be found by multiplying   by the fraction S of the population that is susceptible. When the fraction of the population that is immune increases (i. e. the susceptible population S decreases) so much that   drops below, herd immunity has been achieved and the number of cases occurring in the population will gradually decrease to zero.[15][16][17]

Limitations of

edit

Use of   in the popular press has led to misunderstandings and distortions of its meaning.   can be calculated from many different mathematical models. Each of these can give a different estimate of  , which needs to be interpreted in the context of that model.[10] Therefore, the contagiousness of different infectious agents cannot be compared without recalculating   with invariant assumptions.   values for past outbreaks might not be valid for current outbreaks of the same disease. Generally speaking,   can be used as a threshold, even if calculated with different methods: if  , the outbreak will die out, and if  , the outbreak will expand. In some cases, for some models, values of   can still lead to self-perpetuating outbreaks. This is particularly problematic if there are intermediate vectors between hosts (as is the case for zoonoses), such as malaria.[18] Therefore, comparisons between values from the "Values of   of well-known contagious diseases" table should be conducted with caution.

Although   cannot be modified through vaccination or other changes in population susceptibility, it can vary based on a number of biological, sociobehavioral, and environmental factors.[7] It can also be modified by physical distancing and other public policy or social interventions,[19][7] although some historical definitions exclude any deliberate intervention in reducing disease transmission, including nonpharmacological interventions.[3] And indeed, whether nonpharmacological interventions are included in   often depends on the paper, disease, and what if any intervention is being studied.[7] This creates some confusion, because   is not a constant; whereas most mathematical parameters with "nought" subscripts are constants.

  depends on many factors, many of which need to be estimated. Each of these factors adds to uncertainty in estimates of  . Many of these factors are not important for informing public policy. Therefore, public policy may be better served by metrics similar to  , but which are more straightforward to estimate, such as doubling time or half-life ( ).[20][21]

Methods used to calculate   include the survival function, rearranging the largest eigenvalue of the Jacobian matrix, the next-generation method,[22] calculations from the intrinsic growth rate,[23] existence of the endemic equilibrium, the number of susceptibles at the endemic equilibrium, the average age of infection[24] and the final size equation.[25] Few of these methods agree with one another, even when starting with the same system of differential equations.[18] Even fewer actually calculate the average number of secondary infections. Since   is rarely observed in the field and is usually calculated via a mathematical model, this severely limits its usefulness.[26]

Sample values for various contagious diseases

edit

Despite the difficulties in estimating  mentioned in the previous section, estimates have been made for a number of genera, and are shown in this table. Each genus may be composed of many species, strains, or variants. Estimations of   for species, strains, and variants are typically less accurate than for genera, and so are provided in separate tables below for diseases of particular interest (influenza and COVID-19).

Values of R0 and herd immunity thresholds (HITs) of contagious diseases prior to intervention
Disease Transmission R0 HIT[a]
Measles Aerosol 12–18[27][7] 92–94%
Chickenpox (varicella) Aerosol 10–12[28] 90–92%
Mumps Respiratory droplets 10–12[29] 90–92%
COVID-19 (see values for specific strains below) Respiratory droplets and aerosol 2.9-9.5[30] 65–89%
Rubella Respiratory droplets 6–7[b] 83–86%
Polio Fecal–oral route 5–7[b] 80–86%
Pertussis Respiratory droplets 5.5[35] 82%
Smallpox Respiratory droplets 3.5–6.0[36] 71–83%
HIV/AIDS Body fluids 2–5[37] 50–80%
SARS Respiratory droplets 2–4[38] 50–75%
Diphtheria Saliva 2.6 (1.74.3)[39] 62% (4177%)
Common cold (e.g., rhinovirus) Respiratory droplets 2–3[40][medical citation needed] 50–67%
Mpox Physical contact, body fluids, respiratory droplets, sexual (MSM) 2.1 (1.12.7)[41][42] 53% (2263%)
Ebola (2014 outbreak) Body fluids 1.8 (1.41.8)[43] 44% (3144%)
Influenza (seasonal strains) Respiratory droplets 1.3 (1.21.4)[44] 23% (1729%)
Andes hantavirus Respiratory droplets and body fluids 1.2 (0.81.6)[45] 16% (036%)[c]
Nipah virus Body fluids 0.5[46] 0%[c]
MERS Respiratory droplets 0.5 (0.30.8)[47] 0%[c]

Estimates for strains of influenza.

Values of R0 and herd immunity thresholds (HITs) for specific influenza strains
Disease Transmission R0 HIT[a]
Influenza (1918 pandemic strain) Respiratory droplets 2[48] 50%
Influenza (2009 pandemic strain) Respiratory droplets 1.6 (1.32.0)[2] 37% (2551%)
Influenza (seasonal strains) Respiratory droplets 1.3 (1.21.4)[44] 23% (1729%)

Estimates for variants of SARS-CoV-2.

Values of R0 and herd immunity thresholds (HITs) for variants of SARS-CoV-2
Disease Transmission R0 HIT[a]
COVID-19 (Omicron variant) Respiratory droplets and aerosol 9.5[30] 89%
COVID-19 (Delta variant) Respiratory droplets and aerosol 5.1[49] 80%
COVID-19 (Alpha variant) Respiratory droplets and aerosol 4–5[50][medical citation needed] 75–80%
COVID-19 (ancestral strain) Respiratory droplets and aerosol[51] 2.9 (2.43.4)[52] 65% (5871%)


edit

In the 2011 film Contagion, a fictional medical disaster thriller, a blogger's calculations for   are presented to reflect the progression of a fatal viral infection from isolated cases to a pandemic.[19]

See also

edit

Notes

edit
  1. ^ a b c Calculated using p = 1 − 1/R0.
  2. ^ a b From a module of a training course[31] with data modified from other sources.[32][33][34]
  3. ^ a b c When R0 < 1.0, the disease naturally disappears.

References

edit
  1. ^ Milligan GN, Barrett AD (2015). Vaccinology : an essential guide. Chichester, West Sussex: Wiley Blackwell. p. 310. ISBN 978-1-118-63652-7. OCLC 881386962.
  2. ^ a b Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, et al. (June 2009). "Pandemic potential of a strain of influenza A (H1N1): early findings". Science. 324 (5934): 1557–61. Bibcode:2009Sci...324.1557F. doi:10.1126/science.1176062. PMC 3735127. PMID 19433588.
  3. ^ a b Becker NG, Glass K, Barnes B, Caley P, Philp D, McCaw JM, et al. (April 2006). "The reproduction number". Using Mathematical Models to Assess Responses to an Outbreak of an Emerged Viral Respiratory Disease. National Centre for Epidemiology and Population Health. ISBN 1-74186-357-0. Archived from the original on February 1, 2020. Retrieved February 1, 2020.
  4. ^ Adam D (July 2020). "A guide to R - the pandemic's misunderstood metric". Nature. 583 (7816): 346–348. Bibcode:2020Natur.583..346A. doi:10.1038/d41586-020-02009-w. PMID 32620883.
  5. ^ Jones J. "Notes On R0" (PDF). Stanford University.
  6. ^ Siegel E. "Why 'Exponential Growth' Is So Scary For The COVID-19 Coronavirus". Forbes. Retrieved March 19, 2020.
  7. ^ a b c d e Delamater PL, Street EJ, Leslie TF, Yang YT, Jacobsen KH (January 2019). "Complexity of the Basic Reproduction Number (R0)". Emerging Infectious Diseases. 25 (1): 1–4. doi:10.3201/eid2501.171901. PMC 6302597. PMID 30560777.
  8. ^ Fine, P.; Eames, K.; Heymann, D. L. (April 1, 2011). "'Herd Immunity': A Rough Guide". Clinical Infectious Diseases. 52 (7): 911–916. doi:10.1093/cid/cir007. PMID 21427399.
  9. ^ a b Hiraoka, Takayuki; K. Rizi, Abbas; Kivelä, Mikko; Saramäki, Jari (May 12, 2022). "Herd immunity and epidemic size in networks with vaccination homophily". Physical Review E. 105 (5): L052301. arXiv:2112.07538. Bibcode:2022PhRvE.105e2301H. doi:10.1103/PhysRevE.105.L052301. PMID 35706197. S2CID 245130970.
  10. ^ a b Vegvari C. "Commentary on the use of the reproduction number R during the COVID-19 pandemic". Stat Methods Med Res. PMID 34569883.
  11. ^ Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE (April 5, 2012). "Ross, macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens". PLOS Pathogens. 8 (4): e1002588. doi:10.1371/journal.ppat.1002588. PMC 3320609. PMID 22496640.
  12. ^ Macdonald G (September 1952). "The analysis of equilibrium in malaria". Tropical Diseases Bulletin. 49 (9): 813–29. PMID 12995455.
  13. ^ Network Science by Albert-László Barabási.
  14. ^ Diekmann O, Heesterbeek JA, Metz JA (1990). "On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations". Journal of Mathematical Biology. 28 (4): 365–82. doi:10.1007/BF00178324. hdl:1874/8051. PMID 2117040. S2CID 22275430.
  15. ^ Garnett GP (February 2005). "Role of herd immunity in determining the effect of vaccines against sexually transmitted disease". The Journal of Infectious Diseases. 191 (Suppl 1): S97-106. doi:10.1086/425271. PMID 15627236.
  16. ^ Rodpothong P, Auewarakul P (October 2012). "Viral evolution and transmission effectiveness". World Journal of Virology. 1 (5): 131–4. doi:10.5501/wjv.v1.i5.131. PMC 3782273. PMID 24175217.
  17. ^ Dabbaghian V, Mago VK (2013). Theories and Simulations of Complex Social Systems. Springer. pp. 134–35. ISBN 978-3642391491. Retrieved March 29, 2015.
  18. ^ a b Li J, Blakeley D, Smith RJ (2011). "The failure of R0". Computational and Mathematical Methods in Medicine. 2011 (527610): 527610. doi:10.1155/2011/527610. PMC 3157160. PMID 21860658.
  19. ^ a b Byrne M (October 6, 2014), "The Misunderstood Number That Predicts Epidemics", vice.com, retrieved March 23, 2020
  20. ^ Balkew TM (December 2010). The SIR Model When S(t) is a Multi-Exponential Function (Thesis). East Tennessee State University.
  21. ^ Ireland MW, ed. (1928). The Medical Department of the United States Army in the World War, vol. IX: Communicable and Other Diseases. Washington: U.S.: U.S. Government Printing Office. pp. 116–7.
  22. ^ Diekmann O, Heesterbeek JA (2000). "The Basic Reproduction Ratio". Mathematical Epidemiology of Infectious Diseases : Model Building, Analysis and Interpretation. New York: Wiley. pp. 73–98. ISBN 0-471-49241-8.
  23. ^ Chowell G, Hengartner NW, Castillo-Chavez C, Fenimore PW, Hyman JM (July 2004). "The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda". Journal of Theoretical Biology. 229 (1): 119–26. arXiv:q-bio/0503006. Bibcode:2004JThBi.229..119C. doi:10.1016/j.jtbi.2004.03.006. PMID 15178190. S2CID 7298792.
  24. ^ Ajelli M, Iannelli M, Manfredi P, Ciofi degli Atti ML (March 2008). "Basic mathematical models for the temporal dynamics of HAV in medium-endemicity Italian areas". Vaccine. 26 (13): 1697–707. doi:10.1016/j.vaccine.2007.12.058. PMID 18314231.
  25. ^ von Csefalvay, Chris (January 1, 2023), von Csefalvay, Chris (ed.), "2 - Simple compartmental models: The bedrock of mathematical epidemiology", Computational Modeling of Infectious Disease, Academic Press, pp. 19–91, doi:10.1016/b978-0-32-395389-4.00011-6, ISBN 978-0-323-95389-4, retrieved March 2, 2023
  26. ^ Heffernan JM, Smith RJ, Wahl LM (September 2005). "Perspectives on the basic reproductive ratio". Journal of the Royal Society, Interface. 2 (4): 281–93. doi:10.1098/rsif.2005.0042. PMC 1578275. PMID 16849186.
  27. ^ Guerra FM, Bolotin S, Lim G, Heffernan J, Deeks SL, Li Y, Crowcroft NS (December 2017). "The basic reproduction number (R0) of measles: a systematic review". The Lancet. Infectious Diseases. 17 (12): e420–e428. doi:10.1016/S1473-3099(17)30307-9. PMID 28757186.
  28. ^ Ireland's Health Services. Health Care Worker Information (PDF). Retrieved March 27, 2020.
  29. ^ Australian government Department of Health Mumps Laboratory Case Definition (LCD)
  30. ^ a b Liu, Y (March 9, 2022). "The effective reproductive number of the Omicron variant of SARS-CoV-2 is several times relative to Delta". Journal of Travel Medicine. 29 (3). Table 1. doi:10.1093/jtm/taac037. ISSN 1708-8305. PMC 8992231. PMID 35262737.
  31. ^ Centers for Disease Control and Prevention; World Health Organization (2001). "History and epidemiology of global smallpox eradication". Smallpox: disease, prevention, and intervention (training course) (Presentation). Atlanta: Centers for Disease Control and Prevention (published August 25, 2014). cdc:27929. Archived (PDF) from the original on March 17, 2017. Retrieved June 17, 2021.
  32. ^ Fine, Paul E. M. (1993). "Herd Immunity: History, Theory, Practice". Epidemiologic Reviews. 15 (2): 265–302. doi:10.1093/oxfordjournals.epirev.a036121. PMID 8174658.
  33. ^ Luman, ET; Barker, LE; Simpson, DM; Rodewald, LE; Szilagyi, PG; Zhao, Z (May 2001). "National, state, and urban-area vaccination-coverage levels among children aged 19–35 months, United States, 1999". American Journal of Preventive Medicine. 20 (4): 88–153. doi:10.1016/s0749-3797(01)00274-4. PMID 12174806.
  34. ^ Jiles, RB; Fuchs, C; Klevens, RM (September 22, 2000). "Vaccination coverage among children enrolled in Head Start programs or day care facilities or entering school". Morbidity and Mortality Weekly Report. 49 (9): 27–38. PMID 11016876.
  35. ^ Kretzschmar M, Teunis PF, Pebody RG (June 2010). "Incidence and reproduction numbers of pertussis: estimates from serological and social contact data in five European countries". PLOS Medicine. 7 (6): e1000291. doi:10.1371/journal.pmed.1000291. PMC 2889930. PMID 20585374.
  36. ^ Gani R, Leach S (December 2001). "Transmission potential of smallpox in contemporary populations". Nature. 414 (6865): 748–51. Bibcode:2001Natur.414..748G. doi:10.1038/414748a. PMID 11742399. S2CID 52799168. Retrieved March 18, 2020.
  37. ^ "Playing the Numbers Game: R0". National Emerging Special Pathogen Training and Education Center. January 30, 2020. Archived from the original on May 12, 2020. Retrieved December 27, 2020. [...] while infections that require sexual contact like HIV have a lower R0 (2-5).
  38. ^ Consensus document on the epidemiology of severe acute respiratory syndrome (SARS). Department of Communicable Disease Surveillance and Response (Technical report). World Health Organization. p. 26. hdl:10665/70863. WHO/CDS/CSR/GAR/2003.11. A number of researchers have estimated the basic reproduction number by fitting models to the initial growth of epidemics in a number of countries. Their observations indicate that the SARS-CoV is less transmissible than initially thought with estimates of Ro in the range of 2-4.
  39. ^ Truelove SA, Keegan LT, Moss WJ, Chaisson LH, Macher E, Azman AS, Lessler J (June 2020). "Clinical and Epidemiological Aspects of Diphtheria: A Systematic Review and Pooled Analysis". Clinical Infectious Diseases. 71 (1): 89–97. doi:10.1093/cid/ciz808. PMC 7312233. PMID 31425581.
  40. ^ Freeman C (November 6, 2014). "Magic formula that will determine whether Ebola is beaten". The Telegraph. Telegraph.Co.Uk. Archived from the original on January 12, 2022. Retrieved March 30, 2020.
  41. ^ Grant R, Nguyen LL, Breban R (September 1, 2020). "Modelling human-to-human transmission of monkeypox" (PDF). Bulletin of the World Health Organization. 98 (9): 638–640. doi:10.2471/BLT.19.242347. ISSN 0042-9686. PMC 7463189. PMID 33012864. Archived from the original (PDF) on December 11, 2020.
  42. ^ Al-Raeei M (February 2023). "The study of human monkeypox disease in 2022 using the epidemic models: herd immunity and the basic reproduction number case". Annals of Medicine & Surgery. 85 (2): 316–321. doi:10.1097/MS9.0000000000000229. ISSN 2049-0801. PMC 9949786. PMID 36845803.
  43. ^ Wong ZS, Bui CM, Chughtai AA, Macintyre CR (April 2017). "A systematic review of early modelling studies of Ebola virus disease in West Africa". Epidemiology and Infection. 145 (6): 1069–1094. doi:10.1017/S0950268817000164. PMC 9507849. PMID 28166851. The median of the R0 mean estimate for the ongoing epidemic (overall) is 1.78 (interquartile range: 1.44, 1.80)
  44. ^ a b Chowell G, Miller MA, Viboud C (June 2008). "Seasonal influenza in the United States, France, and Australia: transmission and prospects for control". Epidemiology and Infection. 136 (6). Cambridge University Press: 852–64. doi:10.1017/S0950268807009144. PMC 2680121. PMID 17634159. The reproduction number across influenza seasons and countries lied in the range 0.9–2.0 with an overall mean of 1.3, and 95% confidence interval (CI) 1.2–1.4.
  45. ^ Martínez, Valeria P.; Di Paola, Nicholas; Alonso, Daniel O.; Pérez-Sautu, Unai; Bellomo, Carla M.; Iglesias, Ayelén A.; et al. (December 3, 2020). "'Super-Spreaders' and Person-to-Person Transmission of Andes Virus in Argentina". New England Journal of Medicine. 383 (23): 2230–2241. doi:10.1056/NEJMoa2009040. PMID 33264545. S2CID 227259435.
  46. ^ Luby SP (October 2013). "The pandemic potential of Nipah virus". Antiviral Research. 100 (1): 38–43. doi:10.1016/j.antiviral.2013.07.011. PMID 23911335.
  47. ^ Kucharski AJ, Althaus CL (June 2015). "The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission". Euro Surveillance. 20 (25): 14–8. doi:10.2807/1560-7917.ES2015.20.25.21167. PMID 26132768.
  48. ^ "Omicron transmission: how contagious diseases spread". Nebraska Medicine. December 21, 2021. Retrieved January 25, 2022.
  49. ^ Liu, Ying; Rocklöv, Joacim (October 1, 2021). "The reproductive number of the Delta variant of SARS-CoV-2 is far higher compared to the ancestral SARS-CoV-2 virus". Journal of Travel Medicine. 28 (7). doi:10.1093/jtm/taab124. ISSN 1708-8305. PMC 8436367. PMID 34369565.
  50. ^ Gallagher, James (June 12, 2021). "Covid: Is there a limit to how much worse variants can get?". BBC News. Retrieved July 21, 2021.
  51. ^ Prather, Kimberly A.; Marr, Linsey C.; Schooley, Robert T.; McDiarmid, Melissa A.; Wilson, Mary E.; Milton, Donald K. (October 16, 2020). "Airborne transmission of SARS-CoV-2". Science. 370 (6514): 303.2–304. Bibcode:2020Sci...370..303P. doi:10.1126/science.abf0521. PMID 33020250. S2CID 222145689.
  52. ^ Billah, Arif; Miah, Mamun; Khan, Nuruzzaman (November 11, 2020). "Reproductive number of coronavirus: A systematic review and meta-analysis based on global level evidence". PLOS ONE. 15 (11): e0242128. Bibcode:2020PLoSO..1542128B. doi:10.1371/journal.pone.0242128. PMC 7657547. PMID 33175914.

Further reading

edit