Mammarenaviruses of Rodents, South Africa and Zimbabwe.

We conducted a survey for group-specific indirect immunofluorescence antibody to mammarenaviruses by using Lassa fever and Mopeia virus antigens on serum specimens of 5,363 rodents of 33 species collected in South Africa and Zimbabwe during 1964-1994. Rodents were collected for unrelated purposes or for this study and stored at -70°C. We found antibody to be widely distributed in the 2 countries; antibody was detected in serum specimens of 1.2%-31.8% of 14 species of myomorph rodents, whereas 19 mammarenavirus isolates were obtained from serum specimens and viscera of 4 seropositive species. Phylogenetic analysis on the basis of partial nucleoprotein sequences indicates that 14 isolates from Mastomys natalensis, the Natal multimammate mouse, were Mopeia virus, whereas Merino Walk virus was characterized as a novel virus in a separate study. The remaining 4 isolates from 3 rodent species potentially constitute novel viruses pending full characterization.

Persistence of Crimean-Congo Hemorrhagic Fever Virus RNA.

Crimean-Congo hemorrhagic fever virus (CCHFV) causes severe disease with fatalities. Awareness of potential sources of infection is important to reduce risk to healthcare workers and contacts. We detected CCHFV RNA in formalin-fixed, paraffin-embedded tissues from a spontaneous abortion that were submitted for histology 9 weeks after a suspected CCHFV infection in the mother.

Basic insights into Zika virus infection of neuroglial and brain endothelial cells.

Zika virus (ZIKV) has recently emerged as an important human pathogen due to the strong evidence that it causes disease of the central nervous system, particularly microcephaly and Guillain-Barré syndrome. The pathogenesis of disease, including mechanisms of neuroinvasion, may include both invasion via the blood-brain barrier and via peripheral (including cranial) nerves. Cellular responses to infection are also poorly understood. This study characterizes the in vitro infection of laboratory-adapted ZIKV African MR766 and two Asian strains of (1) brain endothelial cells (hCMEC/D3 cell line) and (2) olfactory ensheathing cells (OECs) (the neuroglia populating cranial nerve I and the olfactory bulb; both human and mouse OEC lines) in comparison to kidney epithelial cells (Vero cells, in which ZIKV infection is well characterized). Readouts included infection kinetics, intracellular virus localization, viral persistence and cytokine responses. Although not as high as in Vero cells, viral titres exceeded 104 plaque-forming units (p.f.u.) ml-1 in the endothelial/neuroglial cell types, except hOECs. Despite these substantial titres, a relatively small proportion of neuroglial cells were primarily infected. Immunolabelling of infected cells revealed localization of the ZIKV envelope and NS3 proteins in the cytoplasm; NS3 staining overlapped with that of dsRNA replication intermediate and the endoplasmic reticulum (ER). Infected OECs and endothelial cells produced high levels of pro-inflammatory chemokines. Nevertheless, ZIKV was also able to establish persistent infection in hOEC and hCMEC/D3 cells. Taken together, these results provide basic insights into ZIKV infection of endothelial and neuroglial cells and will form the basis for further study of ZIKV disease mechanisms.

Taxonomy of the order Bunyavirales: second update 2018.

In October 2018, the order Bunyavirales was amended by inclusion of the family Arenaviridae, abolishment of three families, creation of three new families, 19 new genera, and 14 new species, and renaming of three genera and 22 species. This article presents the updated taxonomy of the order Bunyavirales as now accepted by the International Committee on Taxonomy of Viruses (ICTV).

Taxonomy of the order Bunyavirales: update 2019.

In February 2019, following the annual taxon ratification vote, the order Bunyavirales was amended by creation of two new families, four new subfamilies, 11 new genera and 77 new species, merging of two species, and deletion of one species. This article presents the updated taxonomy of the order Bunyavirales now accepted by the International Committee on Taxonomy of Viruses (ICTV).

Global Genomic Diversity of Human Papillomavirus 11 Based on 433 Isolates and 78 Complete Genome Sequences.

UNLABELLED: Human papillomavirus 11 (HPV11) is an etiological agent of anogenital warts and laryngeal papillomas and is included in the 4-valent and 9-valent prophylactic HPV vaccines. We established the largest collection of globally circulating HPV11 isolates to date and examined the genomic diversity of 433 isolates and 78 complete genomes (CGs) from six continents. The genomic variation within the 2,800-bp E5a-E5b-L1-upstream regulatory region was initially studied in 181/207 (87.4%) HPV11 isolates collected for this study. Of these, the CGs of 30 HPV11 variants containing unique single nucleotide polymorphisms (SNPs), indels (insertions or deletions), or amino acid changes were fully sequenced. A maximum likelihood tree based on the global alignment of 78 HPV11 CGs (30 CGs from our study and 48 CGs from GenBank) revealed two HPV11 lineages (lineages A and B) and four sublineages (sublineages A1, A2, A3, and A4). HPV11 (sub)lineage-specific SNPs within the CG were identified, as well as the 208-bp representative region for CG-based phylogenetic clustering within the partial E2 open reading frame and noncoding region 2. Globally, sublineage A2 was the most prevalent, followed by sublineages A1, A3, and A4 and lineage B. IMPORTANCE: This collaborative international study defined the global heterogeneity of HPV11 and established the largest collection of globally circulating HPV11 genomic variants to date. Thirty novel complete HPV11 genomes were determined and submitted to the available sequence repositories. Global phylogenetic analysis revealed two HPV11 variant lineages and four sublineages. The HPV11 (sub)lineage-specific SNPs and the representative region identified within the partial genomic region E2/noncoding region 2 (NCR2) will enable the simpler identification and comparison of HPV11 variants worldwide. This study provides an important knowledge base for HPV11 for future studies in HPV epidemiology, evolution, pathogenicity, prevention, and molecular assay development.

Global genomic diversity of human papillomavirus 6 based on 724 isolates and 190 complete genome sequences.

UNLABELLED: Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Comparative analysis of the L, M, and S RNA segments of Crimean-Congo haemorrhagic fever virus isolates from southern Africa.

Crimean-Congo haemorrhagic fever virus (CCHFV) is a member of the Bunyaviridae family with a tripartite, negative sense RNA genome. This study used predictive software to analyse the L (large), M (medium), and S (small) segments of 14 southern African CCHFV isolates. The OTU-like cysteine protease domain and the RdRp domain of the L segment are highly conserved among southern African CCHFV isolates. The M segment encodes the structural glycoproteins, GN and GC, and the non-structural glycoproteins which are post-translationally cleaved at highly conserved furin and subtilase SKI-1 cleavage sites. All of the sites previously identified were shown to be conserved among southern African CCHFV isolates. The heavily O-glycosylated N-terminal variable mucin-like domain of the M segment shows the highest sequence variability of the CCHFV proteins. Five transmembrane domains are predicted in the M segment polyprotein resulting in three regions internal to and three regions external to the membrane across the G(N), NS(M) and G(C) glycoproteins. The corroboration of conserved genome domains and sequence identity among geographically diverse isolates may assist in the identification of protein function and pathogenic mechanisms, as well as the identification of potential targets for antiviral therapy and vaccine design. As detailed functional studies are lacking for many of the CCHFV proteins, identification of functional domains by prediction of protein structure, and identification of amino acid level similarity to functionally characterised proteins of related viruses or viruses with similar pathogenic mechanisms are a necessary step for selection of areas for further study.

Molecular evolution of viruses of the family Filoviridae based on 97 whole-genome sequences.

Viruses in the Ebolavirus and Marburgvirus genera (family Filoviridae) have been associated with large outbreaks of hemorrhagic fever in human and nonhuman primates. The first documented cases occurred in primates over 45 years ago, but the amount of virus genetic diversity detected within bat populations, which have recently been identified as potential reservoir hosts, suggests that the filoviruses are much older. Here, detailed Bayesian coalescent phylogenetic analyses are performed on 97 whole-genome sequences, 55 of which are newly reported, to comprehensively examine molecular evolutionary rates and estimate dates of common ancestry for viruses within the family Filoviridae. Molecular evolutionary rates for viruses belonging to different species range from 0.46 × 10(-4) nucleotide substitutions/site/year for Sudan ebolavirus to 8.21 × 10(-4) nucleotide substitutions/site/year for Reston ebolavirus. Most recent common ancestry can be traced back only within the last 50 years for Reston ebolavirus and Zaire ebolavirus species and suggests that viruses within these species may have undergone recent genetic bottlenecks. Viruses within Marburg marburgvirus and Sudan ebolavirus species can be traced back further and share most recent common ancestors approximately 700 and 850 years before the present, respectively. Examination of the whole family suggests that members of the Filoviridae, including the recently described Lloviu virus, shared a most recent common ancestor approximately 10,000 years ago. These data will be valuable for understanding the evolution of filoviruses in the context of natural history as new reservoir hosts are identified and, further, for determining mechanisms of emergence, pathogenicity, and the ongoing threat to public health.

Chikungunya: a re-emerging virus.

In the past decade, chikungunya--a virus transmitted by Aedes spp mosquitoes--has re-emerged in Africa, southern and southeastern Asia, and the Indian Ocean Islands as the cause of large outbreaks of human disease. The disease is characterised by fever, headache, myalgia, rash, and both acute and persistent arthralgia. The disease can cause severe morbidity and, since 2005, fatality. The virus is endemic to tropical regions, but the spread of Aedes albopictus into Europe and the Americas coupled with high viraemia in infected travellers returning from endemic areas increases the risk that this virus could establish itself in new endemic regions. This Seminar focuses on the re-emergence of this disease, the clinical manifestations, pathogenesis of virus-induced arthralgia, diagnostic techniques, and various treatment modalities.

Prevalence and occupational exposure to zoonotic diseases in high-risk populations in the Free State Province, South Africa.

Introduction: Zoonotic diseases are responsible for 2.5 billion human cases globally and approximately 2.7 million deaths annually. Surveillance of animal handlers and livestock for zoonotic pathogens contributes to understanding the true disease burden and risk factors within a community. This study investigated the prevalence of selected zoonoses in cattle, farm workers and occupational exposure to endemic zoonotic diseases and their associated risk factors. Methods: Sputum samples from farmworkers were screened for Mycobacterium bovis. Blood specimens from farmworkers and archived sera were tested for serological evidence of Brucella sp., hantaviruses, and Leptospira sp. Communal and commercial cattle herds were tested for bovine tuberculosis and brucellosis. Results: Mycobacterium bovis was not isolated from human samples. A total of 327 human sera were screened, and 35/327 (10.7%) were Brucella sp. IgG positive, 17/327 (5.2%) Leptospira sp. IgM positive, and 38/327 (11.6%) hantavirus IgG positive (95% CI). A higher proportion of Brucella sp. IgG-positive samples were detected among veterinarians (value of p = 0.0006). Additionally, two cattle from a commercial dairy farm were bovine tuberculosis (bTB) positive using the bTB skin test and confirmatory interferon-gamma assay. A higher percentage of confirmed brucellosis-positive animals were from communal herds (8.7%) compared to commercial herds (1.1%). Discussion: These findings highlight the brucellosis and M. bovis prevalence in commercial and communal herds, the zoonotic disease risk in commercial and subsistence farming in developing countries, and the occupational and rural exposure risk to zoonotic pathogens.

Factors affecting the use of biosecurity measures for the protection of ruminant livestock and farm workers against infectious diseases in central South Africa.

Biosecurity measures have been introduced to limit economic losses and zoonotic exposures to humans by preventing and controlling animal diseases. However, they are implemented on individual farms with varying frequency. The goal of this study was to evaluate which biosecurity measures were used by farmers to prevent infectious diseases in ruminant livestock and to identify factors that influenced these decisions. We conducted a survey in 264 ruminant livestock farmers in a 40,000 km2 area in the Free State and Northern Cape provinces of South Africa. We used descriptive statistics, to characterize biosecurity measures and farm attributes, then multivariable binomial regression to assess the strength of the association between the attributes and the implementation of biosecurity measures including property fencing, separate equipment use on different species, separate rearing of species, isolation of sick animals, isolation of pregnant animals, quarantine of new animals, animal transport cleaning, vaccination, tick control and insect control. Ninety-nine percent of farmers reported using at least one of the 10 biosecurity measures investigated (median [M]: 6; range: 0-10). The most frequently used biosecurity measures were tick control (81%, 214 out of 264), vaccination (80%, 211 out of 264) and isolation of sick animals (72%, 190 out of 264). More biosecurity measures were used on farms with 65-282 animals (M: 6; odds ratio [OR]: 1.52) or farms with 283-12,030 animals (M: 7; OR: 1.87) than on farms with fewer than 65 animals (M: 4). Furthermore, farmers who kept two animal species (M: 7; OR: 1.41) or three or more species (M: 7) used more biosecurity measures than single-species operations (M: 4). Farmers with privately owned land used more biosecurity measures (M: 6; OR: 1.51) than those grazing their animals on communal land (M: 3.5). Farms that reported previous Rift Valley fever (RVF) outbreaks used more biosecurity measures (M: 7; OR: 1.25) compared with farms without RVF reports (M: 6) and those that purchased animals in the 12 months prior to the survey (M: 7; OR: 1.19) compared with those that did not (M: 6). When introducing new animals into their herds (n = 122), most farmers used fewer biosecurity measures than they did for their existing herd: 34% (41 out of 122) used multiple biosecurity measures like those of vaccination, tick control, quarantine or antibiotic use, whereas 36% (44 out of 122) used only one and 30% (37 out of 122) used none. Certain farm features, primarily those related to size and commercialization, were associated with more frequent use of biosecurity measures. Given the variation in the application of biosecurity measures, more awareness and technical assistance are needed to support the implementation of a biosecurity management plan appropriate for the type of farm operation and available resources.

Complete genome sequence of a HPV31 isolate from laryngeal squamous cell carcinoma and biological consequences for p97 promoter activity.

Human papillomavirus type 31, although detected less frequently than HPV types 16 and 18, is associated with head and neck squamous cell carcinomas. Previous studies suggest that polymorphisms in the long control region (LCR) may alter the oncogenic potential of the virus. This study reports the first complete genome of a South African HPV31 isolate from a laryngeal squamous cell carcinoma. Sequence variations relative to the HPV31 prototype sequence were identified. The pBlue-Topo® vector, a reporter gene system was used to investigate the possible influence of these variations on the LCR promoter activity in vitro. Using mutagenesis to create two different fragments, ß-galactosidase assays were used to monitor the effect of nucleotide variations on the p97 promoter. Increased ß-galactosidase expression was observed in mutants when compared to the South African HPV31 LCR isolate. Enhanced transcriptional activity was observed with the mutant that possessed a single nucleotide change within the YY1 transcription factor binding site. In conclusion, sequence variation within the LCR of HPV31 isolates may have a functional effect on viral p97 promoter activity.

Whole-Genome Sequence and Comparative Analysis of Human Papillomavirus Type 18 Isolated from a Nasopharyngeal Carcinoma from South Africa.

We report the complete genome sequence of human papillomavirus type 18 isolated from a nasopharyngeal carcinoma in South Africa.

A systematic review and meta-analysis of the sensitivity of antibody tests for the laboratory confirmation of COVID-19.

Aim: The aim of this study was to investigate the utility of serological tests for the diagnosis of COVID-19 during the first week of symptom onset in patients confirmed with the real-time RT-PCR. Materials & methods: A systematic review and meta-analysis of 58 publications were performed using data obtained from Academic Search Ultimate, Africa-wide, Scopus, Web of Science and MEDLINE. Results: We found that the highest pooled sensitivities were obtained with ELISA IgM-IgG and chemiluminescence immunoassay IgM tests. Conclusion: Serological tests have low sensitivity within the first week of symptom onset and cannot replace nucleic acid amplification tests. However, serological assays can be used to support nucleic acid amplification tests.

Arthritogenic alphaviruses: epidemiological and clinical perspective on emerging arboviruses.

Mosquito-borne viruses, or arboviruses, have been part of the infectious disease landscape for centuries, and are often, but not exclusively, endemic to equatorial and subtropical regions of the world. The past two decades saw the re-emergence of arthritogenic alphaviruses, a genus of arboviruses that includes several members that cause severe arthritic disease. Recent outbreaks further highlight the substantial public health burden caused by these viruses. Arthritogenic alphaviruses are often reported in the context of focused outbreaks in specific regions (eg, Caribbean, southeast Asia, and Indian Ocean) and cause debilitating acute disease that can extend to chronic manifestations for years after infection. These viruses are classified among several antigenic complexes, span a range of hosts and mosquito vectors, and can be distributed along specific geographical locations. In this Review, we highlight key features of alphaviruses that are known to cause arthritic disease in humans and outline the present findings pertaining to classification, immunogenicity, pathogenesis, and experimental approaches aimed at limiting disease manifestations. Although the most prominent alphavirus outbreaks in the past 15 years featured chikungunya virus, and a large body of work has been dedicated to understanding chikungunya disease mechanisms, this Review will instead focus on other arthritogenic alphaviruses that have been identified globally and provide a comprehensive appraisal of present and future research directions.

Risk factors associated with exposure to Crimean-Congo haemorrhagic fever virus in animal workers and cattle, and molecular detection in ticks, South Africa.

Crimean-Congo haemorrhagic fever (CCHF) is a severe tick-borne viral zoonosis endemic to parts of Africa, Europe, the Middle East and Central Asia. Human cases are reported annually in South Africa, with a 25% case fatality rate since the first case was recognized in 1981. We investigated CCHF virus (CCHFV) seroprevalence and risk factors associated with infection in cattle and humans, and the presence of CCHFV in Hyalomma spp. ticks in central South Africa in 2017-18. CCHFV IgG seroprevalence was 74.2% (95%CI: 64.2-82.1%) in 700 cattle and 3.9% (95%CI: 2.6-5.8%) in 541 farm and wildlife workers. No veterinary personnel (117) or abattoir workers (382) were seropositive. The prevalence of CCHFV RNA was significantly higher in Hyalomma truncatum (1.6%) than in H. rufipes (0.2%) (P = 0.002). Seroprevalence in cattle increased with age and was greater in animals on which ticks were found. Seroprevalence in cattle also showed significant geographic variation. Seroprevalence in humans increased with age and was greater in workers who handled livestock for injection and collection of samples. Our findings support previous evidence of widespread high CCHFV seroprevalence in cattle and show significant occupational exposure amongst farm and wildlife workers. Our seroprevalence estimate suggests that CCHFV infections are five times more frequent than the 215 confirmed CCHF cases diagnosed in South Africa in the last four decades (1981-2019). With many cases undiagnosed, the potential seriousness of CCHF in people, and the lack of an effective vaccine or treatment, there is a need to improve public health awareness, prevention and disease control.

Comorbidities in SARS-CoV-2 Patients: a Systematic Review and Meta-Analysis.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly spread across the globe at unprecedented speed and is showing no signs of slowing down. The outbreak of coronavirus disease 2019 (COVID-19) has led to significant health burden in infected patients especially in those with underlying comorbidities. The aim of this study was to evaluate the correlation between comorbidities and their role in the exacerbation of disease in COVID-19 patients leading to fatal outcomes. A systematic review was conducted using data from MEDLINE, Scopus, Web of Science, and EMBASE databases published from 1 December 2019 to 15 September 2020. Fifty-three articles were included in the systematic review. Of those 53 articles, 8 articles were eligible for meta-analysis. Hypertension, obesity, and diabetes mellitus were identified to be the most prevalent comorbidities in COVID-19 patients. Our meta-analysis showed that cancer, chronic kidney diseases, diabetes mellitus, and hypertension were independently associated with mortality in COVID-19 patients. Chronic kidney disease was statistically the most prominent comorbidity leading to death. However, despite having high prevalence, obesity was not associated with mortality in COVID-19 patients.IMPORTANCE COVID-19 has plagued the world since it was first identified in December 2019. Previous systematic reviews and meta-analysis were limited by various factors such as the usage of non-peer reviewed data and were also limited by the lack of clinical data on a global scale. Comorbidities are frequently cited as risk factors for severe COVID-19 outcomes. However, the degree to which specific comorbidities impact the disease is debatable. Our study selection involves a global reach and covers all comorbidities that were reported to be involved in the exacerbation of COVID-19 leading to fatal outcomes, which allows us to identify the specific comorbidities that have higher risk in patients. The study highlights COVID-19 high-risk groups. However, further research should focus on the status of comorbidities and prognosis in COVID-19 patients.

A simple and rapid approach to prepare Sindbis and West Nile viral RNA controls for differentiation between positive samples and laboratory contamination.

Reverse transcription-polymerase chain reaction (RT-PCR) is frequently used for surveillance and diagnosis of arboviruses and emerging viruses. A disadvantage of RT-PCR assays, especially nested assays, is the potential for false-positive results caused by laboratory contamination from either positive controls or positive samples. Positive reactors usually require sequence determination for confirmation which delay timeous reporting of a result. Thus, the aim of the study was to use a simple technique to prepare a positive control allowing true positives to be differentiated from laboratory contamination based on size differentiation for conventional PCR, or melt temperatures for real time assays. A flavivirus positive control and an alphavirus positive control were prepared for two RT-PCR assays that we are currently using for arbovirus surveillance in South Africa. Primers targeting a region of the partial genes of interest cloned in pGEM®T-easy were modified at the 5' ends with non-viral nucleotides. The resulting amplicons were circularised, resulting in pGEM®T-easy constructs with 51 and 65 non-viral bases inserted into the partial flaviviral and alphaviral genes respectively and used as template for transcribing RNA. Sequence analysis was used to confirm the manipulation of the partial genes. Using virus specific primer pairs, viral RNA could be readily differentiated from the modified positive controls either by size differentiation, or melt temperature in a SYBR®Green real time RT-PCR. This study demonstrates how simple recombinant technology can be used to produce a positive control that has application in the laboratory for surveillance studies or as a diagnostic tool using synthetic genes to abrogate the requirement for handling infectious virus.

Marburg hemorrhagic fever associated with multiple genetic lineages of virus.

BACKGROUND: An outbreak of Marburg hemorrhagic fever was first observed in a gold-mining village in northeastern Democratic Republic of the Congo in October 1998. METHODS: We investigated the outbreak of Marburg hemorrhagic fever most intensively in May and October 1999. Sporadic cases and short chains of human-to-human transmission continued to occur until September 2000. Suspected cases were identified on the basis of a case definition; cases were confirmed by the detection of virus antigen and nucleic acid in blood, cell culture, antibody responses, and immunohistochemical analysis. RESULTS: A total of 154 cases (48 laboratory-confirmed and 106 suspected) were identified (case fatality rate, 83 percent); 52 percent of cases were in young male miners. Only 27 percent of these men reported having had contact with other affected persons, whereas 67 percent of patients who were not miners reported such contact (P<0.001). Most of the affected miners (94 percent) worked in an underground mine. Cessation of the outbreak coincided with flooding of the mine. Epidemiologic evidence of multiple introductions of infection into the population was substantiated by the detection of at least nine genetically distinct lineages of virus in circulation during the outbreak. CONCLUSIONS: Marburg hemorrhagic fever can have a very high case fatality rate. Since multiple genetic variants of virus were identified, ongoing introduction of virus into the population helped perpetuate this outbreak. The findings imply that reservoir hosts of Marburg virus inhabit caves, mines, or similar habitats.