Abstract
Understanding the factors that mediate overdispersion in SARS-CoV-2 transmissibility is crucial towards mitigating the COVID-19 pandemic. Using a systematically developed dataset, we meta-analyze respiratory viral loads (rVLs) of SARS-CoV-2, SARS-CoV-1 and influenza A(H1N1)pdm09 from 15 countries and model infectiousness by shedding viable virus via droplets and aerosols. Our results indicate heterogeneity in rVL as an intrinsic virological factor facilitating greater overdispersion in the COVID-19 pandemic than in the 2009 H1N1 pandemic. For COVID-19, case heterogeneity remains broad throughout the infectious period, including for pediatric and asymptomatic infections. Many cases inherently present minimal transmission risk, whereas highly infectious individuals shed tens to thousands of SARS-CoV-2 virions/min via droplets and aerosols while breathing, talking and singing. Coughing increases the contagiousness, especially in close contact, of symptomatic cases relative to asymptomatic ones. Our findings show how individual case variations influence SARS-CoV-2 transmissibility and present considerations for disease control.
One Sentence Summary Intrinsic case variation in respiratory viral load facilitates broad overdispersion in SARS-CoV-2 transmissibility.
Main Text
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread globally, causing the coronavirus disease 2019 (COVID-19) pandemic with more than 50.9 million infections and 1.2 million deaths (as of 10 November 2020) (1). For respiratory virus transmission, airway epithelial cells shed virions to the extracellular fluid before atomization (from breathing, talking, singing, coughing and aerosol-generating procedures) partitions them into a polydisperse mixture of particles that are expelled to the ambient environment. Aerosols (≤100 μm) can be inhaled nasally, whereas droplets (>100 μm) tend to be excluded (2, 3). For direct transmission, droplets must be sprayed ballistically onto susceptible tissue (4). Hence, droplets predominantly deposit on nearby surfaces, potentiating indirect transmission. Aerosols can be further categorized based on typical travel characteristics: short-range aerosols (50-100 μm) tend to settle within 2 m; long-range ones (10-50 μm) often travel beyond 2 m based on emission force; and buoyant aerosols (≤10 μm) remain suspended and travel based on airflow profiles for minutes to many hours (4, 5). Although proximity has been associated with infection risk for COVID-19 (6), studies have also suggested that long-range airborne transmission occurs conditionally (7-9).
While the basic reproductive number has been estimated to be 2.0-3.6 (10, 11), transmissibility of SARS-CoV-2 is highly overdispersed (dispersion parameter k, 0.10-0.58), with numerous instances of superspreading (7-9) and few cases (10-20%) causing many secondary infections (80%) (12-14). Similarly, few cases drive the transmission of SARS-CoV-1 (k, 0.16-0.17) (15), whereas influenza A(H1N1)pdm09 transmits more homogeneously (k, 7.4-14.4) (16, 17), despite both viruses spreading by contact, droplets and aerosols (18, 19). Although understanding the determinants of overdispersion is crucial towards characterizing transmissibility and developing effective strategies to limit infection (20), mechanistic associations for k remain unclear. As an empirical estimate, k depends on myriad extrinsic (behavioral, environmental and invention) and host factors. Nonetheless, k remains similar across distinct outbreaks of a virus (15), suggesting that intrinsic virological factors mediate overdispersion for emerging infections.
Here, we investigated how intrinsic case variation in rVL facilitates overdispersion in SARS-CoV-2 transmissibility. By systematic review, we developed a comprehensive dataset of rVLs from cases of COVID-19, SARS and A(H1N1)pdm09. Using comparative meta-analyses, we found that heterogeneity in rVL was associated with overdispersion among these emerging infections. To assess potential sources of case heterogeneity, we analyzed SARS-CoV-2 rVLs across age and symptomatology subgroups as well as disease course. To interpret the influence of heterogeneity in rVL on individual infectiousness, we modelled likelihoods of shedding viable virus via respiratory droplets and aerosols.
Systematic review
We conducted a systematic review on virus quantitation in respiratory specimens taken during the infectious periods of SARS-CoV-2 (−3 to 10 days from symptom onset [DFSO]) (21-23), SARS-CoV-1 (0-20 DFSO) (24) and A(H1N1)pdm09 (−2 to 9 DFSO) (25) (Methods). The systematic search (Tables S1 to S5) identified 4,274 results. After screening and full-text review, 63 studies met the inclusion criteria and contributed to the systematic dataset (Fig. 1) (N = 9,616 total specimens), which included adult (N = 5,030) and pediatric (N = 1,594) cases from 15 countries and specimen measurements for asymptomatic (N = 2,387), presymptomatic (N = 28) and symptomatic (N = 7,146) infections. According to a hybrid Joanna Briggs Institute critical appraisal checklist, risk of bias was low for most contributing studies (Table S6).
Association of overdispersion with heterogeneity in rVL
We hypothesized that individual case variation in rVL facilitates the distinctions in k among COVID-19, SARS and A(H1N1)pdm09. For each study in the systematic dataset, we used specimen measurements (viral RNA concentration in a respiratory specimen) to estimate rVLs (viral RNA concentration in the respiratory tract) (Methods). To investigate the relationship between k and heterogeneity in rVL, we performed a meta-regression using each contributing study (fig. S1), which showed a weak, negative association between the two variables (meta-regression slope t-test: P = 0.040, Pearson’s r = −0.26).
Using contributing studies with low risk of bias, meta-regression (Fig. 2) showed a strong, negative association between k and heterogeneity in rVL (meta-regression slope t-test: P < 10−8, Pearson’s r = −0.73). In this case, each unit increase (1 log10 copies/ml) in the SD of rVL decreased log(k) by a factor of −1.41 (95% confidence interval [CI]: −1.80 to −1.03), suggesting that broader heterogeneity in rVL facilitates greater overdispersion in the transmissibility of SARS-CoV-2 than of A(H1N1)pmd09. To investigate mechanistic aspects of this association, we conducted a series of analyses on rVL and then modelled the influence of heterogeneity in rVL on individual infectiousness.
Meta-analysis and subgroup analyses of rVL
We first compared rVLs among the emerging infections. We performed a random-effects meta-analysis (fig. S2), which assessed the expected rVL when encountering a COVID-19, SARS or A(H1N1)pdm09 case during the infectious period. This showed that the expected rVL of SARS-CoV-2 was not different from that of SARS-CoV-1 (one-sided Welch’s t-test: P = 0.102) but lesser than that of A(H1N1)pdm09 (P = 0.044). We also performed random-effects subgroup analyses for COVID-19 (Fig. 3), which showed that expected SARS-CoV-2 rVLs were similar between pediatric and adult cases (P = 0.409) and comparable between symptomatic/presymptomatic and asymptomatic infections (P = 0.107). Since these meta-analyses had significant between-study heterogeneity among the mean estimates (Cochran’s Q test: P < 0.001 for each meta-analysis), we conducted risk-of-bias sensitivity analyses; meta-analyses of low-risk-of-bias studies continued to show significant heterogeneity (figs. S3 to S7).
Distributions of rVL
We next analyzed rVL distributions. For all three viruses, rVLs best conformed to Weibull distributions (fig. S8), and we fitted the entirety of individual sample data for each virus in the systematic dataset (Fig. 4A and fig. S8N). While COVID-19 and SARS cases tended to shed lesser virus than those with A(H1N1)pdm09 (fig. S2), broad heterogeneity in SARS-CoV-2 and SARS-CoV-1 rVLs inverted this relationship for highly infectious individuals (Fig. 4A and fig. S9, A to C). At the 90th case percentile (cp) throughout the infectious period, the estimated rVL was 8.91 (95% CI: 8.83-9.00) log10 copies/ml for SARS-CoV-2, whereas it was 8.62 (8.47-8.76) log10 copies/ml for A(H1N1)pdm09 (Table S7). The standard deviation (SD) of the overall rVL distribution for SARS-CoV-2 was 2.04 log10 copies/ml, while it was 1.45 log10 copies/ml for A(H1N1)pdm09, showing that heterogeneity in rVL was indeed broader for SARS-CoV-2.
To assess potential sources of heterogeneity in SARS-CoV-2 rVL, we compared rVL distributions among COVID-19 subgroups. In addition to comparable mean estimates (Fig. 3), adult, pediatric, symptomatic/presymptomatic and asymptomatic COVID-19 cases showed similar rVL distributions (Fig. 4, B and C), with SDs of 2.04, 2.06, 2.01 and 2.02 log10 copies/ml, respectively. Thus, age and symptomatology minimally influenced case variation in SARS-CoV-2 rVL during the infectious period.
SARS-CoV-2 kinetics during respiratory infection
To analyze the influences of disease course, we delineated individual SARS-CoV-2 rVLs by DFSO and fitted the mean estimates to a mechanistic model for viral kinetics (Fig. 4D, Methods). The outputs indicated that, on average, each productively infected cell in the airway epithelium shed SARS-CoV-2 at 1.83 (95% CI: 0.69-2.97) copies/ml day-1 and infected up to 9.48 susceptible cells (Table S8). The turnover rate for infected epithelial cells was 0.55 (0.30-0.79) days-1, while the half-life of SARS-CoV-2 in the respiratory tract was 2.65 (1.39-28.2) hours. By extrapolating the model to an initial rVL of 0 log10 copies/ml, the estimated incubation period was 5.09 days, which agrees with epidemiological findings (10). Conversely, the expected duration of shedding was 30.4 DFSO. Thus, SARS-CoV-2 rVL increased exponentially after infection, peaked around 1 DFSO along with the proportion of infected epithelial cells (fig. S10) and then diminished exponentially.
To evaluate case heterogeneity across the infectious period, we fitted distributions for each DFSO (Fig. 4E), which showed that high SARS-CoV-2 rVLs also increased from the presymptomatic period, peaked at 1 DFSO and then decreased towards the end of the first week of illness. For the 90th cp at 1 DFSO, the rVL was 9.83 (95% CI: 9.12-10.61) log10 copies/ml, an order of magnitude greater than the overall 90th-cp estimate. High rVLs between 1-5 DFSO were elevated above the expected values from the overall rVL distribution (Table S7). At −1 DFSO, the 90th-cp rVL was 8.30 (6.88-10.02) log10 copies/ml, while it was 7.96 (7.37-8.58) log10 copies/ml at 10 DFSO. Moreover, heterogeneity in rVL remained broad across the infectious period, with SDs of 1.85-2.46 log10 copies/ml between −1 to 10 DFSO (fig. S9, H to S).
Likelihood that droplets and aerosols contain virions
Towards analyzing the influence of heterogeneity in rVL on individual infectiousness, we first modelled the likelihood of respiratory particles containing viable SARS-CoV-2. Since rVL is an intensive quantity, the volume fraction of virions is low and viral partitioning coincides with atomization, we used Poisson statistics to model likelihood profiles. To calculate an unbiased estimator of partitioning (the expected number of viable copies per particle), our method multiplied rVL estimates with particle volumes during atomization and an assumed viability proportion of 0.1% in equilibrated particles (Methods).
When expelled by the mean COVID-19 case during the infectious period, respiratory particles showed low likelihoods of carrying viable SARS-CoV-2 (fig. S11B). Aerosols (equilibrium aerodynamic diameter [da] ≤ 100 µm) were ≤0.69% (95% CI: 0.43-0.95%) likely to contain a virion. Droplets also had low likelihoods: at da = 330 µm, they were 19.4% (18.5-20.3%), 2.42% (2.15-2.69%) and 0.20% (0.16-0.24%) likely to contain one, two or three virions, respectively.
COVID-19 cases with high rVLs, however, expelled particles with considerably greater likelihoods of carrying viable copies (Fig. 5, A and B, and fig. S11, D and E). For the 80th cp during the infectious period, aerosols (da ≤ 100 µm) were ≤36.6% (95% CI: 35.3-38.0%) likely to carry at least one SARS-CoV-2 virion. For the 90th cp, this likelihood was ≤96.7% (96.5-96.9%), with larger aerosols tending to contain multiple virions (fig. S11E). At 1 DFSO, these estimates were greatest, and ≤65.2% (56.1-74.3%) of buoyant aerosols (da ≤ 10 µm) contained at least one viable copy of SARS-CoV-2 for the 98th cp. When expelled by high cps, droplets (da > 100 µm) tended to contain tens to thousands of SARS-CoV-2 virions (fig. S11E and Fig. 5B)
Shedding SARS-CoV-2 via respiratory droplets and aerosols
Using the partitioning estimates in conjunction with published profiles of the particles expelled by respiratory activities (fig. S12), we next modelled the rates at which talking, singing, breathing and coughing shed viable SARS-CoV-2 across da (Fig. 5, C to F). Singing shed virions more rapidly than talking based on the increased emission of aerosols. Voice amplitude, however, had a significant effect on aerosol production, and talking loudly emitted aerosols at similar rates to singing (fig. S12E). Based on the generation of larger aerosols and droplets, talking and singing shed virions more rapidly than breathing (Fig. 5, C to E). Each cough shed similar quantities of virions as in a minute of talking (Fig. 5, C and F).
Each of these respiratory activities expelled aerosols at greater rates than droplets, but particle size correlated with the likelihood of containing virions. Talking, singing and coughing expelled virions at greater proportions via droplets (80.6-86.0%) than aerosols (14.0-19.4%) (Fig. 5G). Moreover, short-range aerosols predominantly mediated the virions (90.8-92.6%) shed via aerosols while talking normally and coughing. In comparison, while singing, talking loudly and breathing, buoyant (29.7-68.4%) and long-range (10.3-31.6%) aerosols carried a larger proportion of the virions shed via aerosols (Fig. 5G).
Influence of heterogeneity in rVL on individual infectiousness
To interpret how heterogeneity in rVL influences individual infectiousness, we modelled total SARS-CoV-2 shedding rates (over all particle sizes) for each respiratory activity (Fig. 5H and fig. S13). Between the 1st and the 99th cps, the estimates for a respiratory activity spanned ≥8.55 orders of magnitude on each DFSO; cumulatively from −1 to 10 DFSO, they spanned 11.0 orders of magnitude. Hence, many COVID-19 cases inherently presented minimal transmission risk, whereas highly infectious individuals shed considerable quantities of SARS-CoV-2. For the 98th cp at 1 DFSO, singing expelled 350 (95% CI: 36.2-4,213) virions/min to the ambient environment, talking emitted 328 (33.9-3,945) virions/min, breathing exhaled 1.73 (0.18-20.7) virions/min and coughing discharged 279 (28.8-3,350) virions/cough; these estimates were approximately two orders of magnitude greater than those for the 86th cp. For the 98th cp at −1 DFSO, singing shed 14.5 (0.15-4,515) virions/min and breathing exhaled 7.13×10−2 (7.20×10−4-22.2) virions/min. The estimates at 9-10 DFSO were similar to these presymptomatic ones (Fig. 5H and fig. S13B). As indicated by comparable mean rVLs (Fig. 3) and heterogeneities in rVL (Fig. 4, B and C), adult, pediatric, symptomatic/presymptomatic and asymptomatic COVID-19 subgroups presented similar distributions for shedding virions through these activities.
We also compared the influence of case variation on individual infectiousness between A(H1N1)pdm09 and COVID-19. Aerosol spread accounted for approximately half of A(H1N1)pdm09 transmission events (19), and the 50% human infectious dose for aerosolized influenza A virus is approximately 1-3 virions in the absence of neutralizing antibodies (26). Based on the model, 62.9% of A(H1N1)pdm09 cases were infectious (shed ≥1 virion) via aerosols within 24 h of talking loudly or singing (fig. S14A). The estimate was 58.6% within 24 h of talking normally and 22.3% within 24 h of breathing. In comparison, 48.0% of COVID-19 cases shed ≥1 virion via aerosols in 24 h of talking loudly or singing (fig. S14C). Notably, only 61.4% of COVID-19 cases shed ≥1 virion via either droplets or aerosols in 24 h of talking loudly or singing (fig. S14D). While the human infectious dose of SARS-CoV-2 by any exposure route remains unelucidated, it must be at least one viable copy. Thus, at least 38.6% of COVID-19 cases were expected to present negligible risk to spread SARS-CoV-2 through either droplets or aerosols in 24 h. The proportion of inherently infectious cases further decreased as the infectious dose increased: 55.8, 42.5 and 25.0% of COVID-19 cases were expected to shed ≥2, ≥10 and ≥100 virions, respectively, in 24 h of talking loudly or singing.
While these analyses indicated that a greater proportion of A(H1N1)pdm09 cases were inherently infectious, 18.8% of COVID-19 cases shed virions more rapidly than those infected with A(H1N1)pdm09 (Fig. 4A). At the 98th cp for A(H1N1)pdm09, singing expelled 4.38 (2.85-6.78) virions/min and breathing exhaled 2.15×10−2 (1.40×10−2-3.34×10−2) virions/min. Highly infectious COVID-19 cases expelled virions at rates that were up to 1-2 orders of magnitude greater than their A(H1N1)pdm09 counterparts (Fig. 5H and fig. S15).
Discussion
This study provided systematic analyses of several factors characterizing SARS-CoV-2 transmissibility. First, we found that broader heterogeneity in rVL facilitates greater overdispersion for SARS-CoV-2 than A(H1N1)pdm09. Our results suggest that many COVID-19 cases infect no one (12-14) because they inherently present minimal transmission risk via droplets or aerosols, although behavioral and environmental factors can further abate risk. Meanwhile, highly infectious cases can shed tens to thousands of SARS-CoV-2 virions/min. The model estimates, when corrected to copies rather than virions, align with recent clinical findings for exhalation rates of SARS-CoV-2 (27). In comparison, a greater proportion of A(H1N1)pdm09 cases are infectious but shed virions at low rates, which concurs with more uniform transmission and few superspreading events observed during the 2009 H1N1 pandemic (16, 17). Moreover, our analyses suggest that heterogeneity in rVL may be generally associated with overdispersion for viral respiratory infections. In this case, rVL distribution can serve as an early correlate for transmission patterns, including superspreading, during outbreaks of novel respiratory viruses, providing insight for disease control before large-scale epidemiological studies empirically characterize k. When transmission is highly overdispersed, targeted interventions may disproportionately mitigate infection (20), with models estimating that focusing half of control efforts on the most infectious 20% of cases outperforms random control policies threefold (15).
Second, we analyzed SARS-CoV-2 kinetics during respiratory infection. While heterogeneity remains broad throughout the infectious period, rVL tends to peak at 1 DFSO and be elevated for 1-5 DFSO, coinciding with the period of highest attack rates observed among close contacts (28). These results indicate that transmission risk tends to be greatest soon after illness rather than in the presymptomatic period, which concurs with large tracing studies (6.4-12.6% of secondary infections from presymptomatic transmission) (29, 30) rather than early temporal models (∼44%) (23). Furthermore, our kinetic analysis suggests that, on average, SARS-CoV-2 reaches diagnostic concentrations 1.58-3.03 days after respiratory infection (−3.51 to −2.06 DFSO), assuming assay detection limits of 1-3 log10 copies/ml, respectively, for nasopharyngeal swabs immersed in 1 ml of transport media.
Third, we assessed the relative infectiousness of COVID-19 subgroups. As a common symptom of COVID-19 (31), coughing sheds considerable numbers of virions via droplets and short-range aerosols. Thus, symptomatic infections tend to be more contagious than asymptomatic ones, providing one reason as to why asymptomatic cases transmit SARS-CoV-2 at lower relative rates (32), especially in close contact (33), despite similar rVLs and increased contact patterns. Accordingly, children (48-54% of symptomatic cases present with cough) (34, 35) tend to be less contagious than adults (68-80%) (31, 35) based on tendencies of symptomatology rather than rVL. Conversely, coughing sheds few virions via smaller aerosols. Our analyses suggest that asymptomatic and symptomatic infections present comparable risks for airborne spread, as do adult and pediatric cases. While singing and talking loudly, highly infectious individuals shed tens to hundreds of SARS-CoV-2 virions/min via long-range and buoyant aerosols.
Our study has limitations. The systematic search found a limited number of studies reporting quantitative specimen measurements from the presymptomatic period, meaning these estimates may be sensitive to sampling bias. Although additional studies have reported semiquantitative metrics (cycle thresholds), these data were excluded because they cannot be compared on an absolute scale due to batch effects (36), limiting use in compound analyses. Furthermore, this study considered population-level estimates of the infectious periods, viability proportions and rate profiles for respiratory particles, which omit individual or environmental variation. Distinctions in phonetic tendencies and, especially for young children, respiratory capacity lead to variation in particle emission rates (37). Some patients shed SARS-CoV-2 with diminishing viability soon after symptom onset (21), whereas others produce replication-competent virus for weeks (38). It remains unclear how case characteristics and environmental factors affect the viability dynamics of SARS-CoV-2.
Taken together, our findings provide a potential path forward for disease control. They support aerosol spread as a transmission mode for SARS-CoV-2, including for conditional superspreading by highly infectious cases. However, with short durations of stay in well-ventilated areas, the exposure risk for aerosols, including long-range and buoyant ones, remains correlated with proximity to infectious cases (2, 4). Strategies to abate infection should limit crowd numbers and duration of stay while reinforcing distancing, low voice amplitudes and widespread mask usage; well-ventilated settings can be recognized as lower risk venues. Coughing can shed considerable quantities of virions, while rVL tends to peak at 1 DFSO and can be high throughout the infectious period. Thus, immediate, sustained self-isolation upon illness is crucial to curb transmission from symptomatic cases. Collectively, our analyses highlight the role of infections with high rVLs in propelling the COVID-19 pandemic. While diagnosing COVID-19, qRT-PCR can also triage contact tracing, prioritizing these patients: for nasopharyngeal swabs immersed in 1 ml of transport media, ≥7.14 (95% CI: 7.07-7.22) log10 copies/ml corresponds to the top 20% of COVID-19 cases. Doing so may identify asymptomatic and presymptomatic infections more efficiently, a key step towards mitigation as the pandemic continues.
Data Availability
Data will be made available upon request. All raw data, code and model outputs from this study will be made publicly available in online repositories after peer review. Search strategies for the systematic review are shown in Supplementary Tables 1-5. The systematic review protocol was prospectively registered on PROSPERO (registration number, CRD42020204637).
Funding
This work was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Toronto COVID-19 Action Fund. P.Z.C. was supported by the NSERC Vanier Scholarship (608544). D.N.F. was supported by the Canadian Institutes of Health Research (Canadian COVID-19 Rapid Research Fund, OV4-170360). F.X.G. was supported by the NSERC Senior Industrial Research Chair.
Author contributions
Conceptualization, P.Z.C.; Methodology, P.Z.C., N.B., and Z.P.; Formal Analysis, P.Z.C.; Investigation, P.Z.C. and N.B.; Resources, Z.P.; Writing – Original Draft, P.Z.C.; Writing – Review & Editing, P.Z.C., N.B., Z.P., M.K., D.N.F., and F.X.G.; Visualization, P.Z.C.; Supervision, M.K., D.N.F., and F.X.G; Funding Acquisition, F.X.G.
Data and materials availability
The systematic dataset and model outputs from this study are uploaded to Zenodo (https://zenodo.org/record/4266560). The code generated during this study is available at GitHub (https://github.com/paulzchen/sars2-heterogeneity). Search strategies for the systematic review are shown in Tables S1 to S5. The systematic review protocol was prospectively registered on PROSPERO (registration number, CRD42020204637).
Competing interests
The authors declare no competing interests.
Acknowledgments
We thank T. Alba (Toronto) for discussion on statistical methods. We thank J. Jimenez (Colorado) for discussion on the characteristics of aerosols and droplets. We thank E. Lavezzo and A. Chrisanti (Padova) and A. Wyllie and N. Grubaugh (Yale) for responses to data inquiries.
References and Notes
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.
- 40.
- 41.
- 42.
- 43.
- 44.
- 45.
- 46.
- 47.
- 48.
- 49.
- 50.
- 51.
- 52.
- 53.
- 54.
- 55.
- 56.
- 57.
- 58.
- 59.
- 60.
- 61.
- 62.
- 63.
- 64.
- 65.
- 66.
- 67.
- 68.
- 69.
- 70.
- 71.
- 72.
- 73.
- 74.
- 75.
- 76.
- 77.
- 78.
- 79.
- 80.
- 81.
- 82.
- 83.
- 84.
- 85.
- 86.
- 87.
- 88.
- 89.
- 90.
- 91.
- 92.
- 93.
- 94.
- 95.
- 96.
- 97.
- 98.
- 99.
- 100.
- 101.
- 102.
- 103.
- 104.
- 105.
- 106.
- 107.
- 108.
- 109.
- 110.
- 111.
- 112.
- 113.
- 114.
- 115.
- 116.
- 117.
- 118.
- 119.