High-resolution phylogeographical surveillance of Hantaan orthohantavirus using rapid amplicon-based Flongle sequencing, Republic of Korea.

Orthohantaviruses, etiological agents of hemorrhagic fever with renal syndrome (HFRS) and hantavirus cardiopulmonary syndrome, pose a critical public health threat worldwide. Hantaan orthohantavirus (HTNV) outbreaks are particularly endemic in Gyeonggi Province in northern area of the Republic of Korea (ROK). Small mammals were collected from three regions in the Gyeonggi Province during 2017 and 2018. Serological and molecular prevalence of HTNV was 25/201 (12.4%) and 10/25 (40%), respectively. A novel nanopore-based diagnostic assay using a cost-efficient Flongle chip was developed to rapidly and sensitively detect HTNV infection in rodent specimens within 3 h. A rapid phylogeographical surveillance of HTNV at high-resolution phylogeny was established using the amplicon-based Flongle sequencing. In total, seven whole-genome sequences of HTNV were newly obtained from wild rodents collected in Paju-si (Gaekhyeon-ri) and Yeoncheon-gun (Hyeonga-ri and Wangnim-ri), Gyeonggi Province. Phylogenetic analyses revealed well-supported evolutionary divergence and genetic diversity, enhancing the resolution of the phylogeographic map of orthohantaviruses in the ROK. Incongruences in phylogenetic patterns were identified among HTNV tripartite genomes, suggesting differential evolution for each segment. These findings provide crucial insights into on-site diagnostics, genome-based surveillance, and the evolutionary dynamics of orthohantaviruses to mitigate hantaviral outbreaks in HFRS-endemic areas in the ROK.

A multipathogen DNA vaccine elicits protective immune responses against two class A bioterrorism agents, anthrax and botulism.

The potential use of biological agents has become a major public health concern worldwide. According to the CDC classification, Bacillus anthracis and Clostridium botulinum, the bacterial pathogens that cause anthrax and botulism, respectively, are considered to be the most dangerous potential biological agents. Currently, there is no licensed vaccine that is well suited for mass immunization in the event of an anthrax or botulism epidemic. In the present study, we developed a dual-expression system-based multipathogen DNA vaccine that encodes the PA-D4 gene of B. anthracis and the HCt gene of C. botulinum. When the multipathogen DNA vaccine was administered to mice and guinea pigs, high level antibody responses were elicited against both PA-D4 and HCt. Analysis of the serum IgG subtype implied a combined Th1/Th2 response to both antigens, but one that was Th2 skewed. In addition, immunization with the multipathogen DNA vaccine induced effective neutralizing antibody activity against both PA-D4 and HCt. Finally, the protection efficiency of the multipathogen DNA vaccine was determined by sequential challenge with 10 LD50 of B. anthracis spores and 10 LD50 of botulinum toxin, or vice versa, and the multipathogen DNA vaccine provided higher than 50% protection against lethal challenge with both high-risk biothreat agents. Our studies suggest the strategy used for this anthrax-botulinum multipathogen DNA vaccine as a prospective approach for developing emergency vaccines that can be immediately distributed on a massive scale in response to a biothreat emergency or infectious disease outbreak. Key points • A novel multipathogen DNA vaccine was constructed against anthrax and botulism. • Robust immune responses were induced following vaccination. • Suggests a potential vaccine development strategy against biothreat agents.

Active Targeted Surveillance to Identify Sites of Emergence of Hantavirus.

BACKGROUND: Endemic outbreaks of hantaviruses pose a critical public health threat worldwide. Hantaan orthohantavirus (HTNV) causes hemorrhagic fever with renal syndrome (HFRS) in humans. Using comparative genomic analyses of partial and nearly complete sequences of HTNV from humans and rodents, we were able to localize, with limitations, the putative infection locations for HFRS patients. Partial sequences might not reflect precise phylogenetic positions over the whole-genome sequences; finer granularity of rodent sampling reflects more precisely the circulation of strains. METHODS: Five HFRS specimens were collected. Epidemiological surveys were conducted with the patients during hospitalization. We conducted active surveillance at suspected HFRS outbreak areas. We performed multiplex polymerase chain reaction-based next-generation sequencing to obtain the genomic sequence of HTNV from patients and rodents. The phylogeny of human- and rodent-derived HTNV was generated using the maximum likelihood method. For phylogeographic analyses, the tracing of HTNV genomes from HFRS patients was defined on the bases of epidemiological interviews, phylogenetic patterns of the viruses, and geographic locations of HTNV-positive rodents. RESULTS: The phylogeographic analyses demonstrated genetic clusters of HTNV strains from clinical specimens, with HTNV circulating in rodents at suspected sites of patient infections. CONCLUSIONS: This study demonstrates a major shift in molecular epidemiological surveillance of HTNV. Active targeted surveillance was performed at sites of suspected infections, allowing the high-resolution phylogeographic analysis to reveal the site of emergence of HTNV. We posit that this novel approach will make it possible to identify infectious sources, perform disease risk assessment, and implement preparedness against vector-borne viruses.

Prolonged shedding of type 55 human adenovirus in immunocompetent adults with adenoviral respiratory infections.

Human adenovirus (HAdV) is a common pathogen causing respiratory infections with outbreaks reported in the military and community. However, little information is available on the shedding kinetics. We performed a prospective study of immunocompetent adults confirmed with HAdV respiratory infection by multiplex real-time PCR during an outbreak of HAdV-55. Consecutive respiratory specimens of sputum or nasopharyngeal swab were collected from each patient every 2 days. Viral load was measured by real-time quantitative PCR. Of 32 enrolled patients, 27 (84.4%) had pneumonia. Five patients (15.6%) received cidofovir. Viral load was highest in the earliest samples at 8.69 log10 copies/mL. In a linear regression model, viral load declined consistently in a log-linear fashion at the rate of - 0.15 log10 copies/mL per day (95% confidence interval (CI): - 0.18, - 0.12; R2 = 0.32). However, the regression model estimated the viral shedding duration to be 55 days. The rate of decline in viral load did not differ between patients who received cidofovir and who did not. Patients with prominent respiratory symptoms or extensive involvement on chest radiograph had higher volume of viral excretion. Prolonged viral shedding was observed in otherwise healthy adults with HAdV-55 respiratory infection. This finding should be considered in the establishment of infection control and prevention strategies.

Multiplex PCR-Based Next-Generation Sequencing and Global Diversity of Seoul Virus in Humans and Rats.

Seoul virus (SEOV) poses a worldwide public health threat. This virus, which is harbored by Rattus norvegicus and R. rattus rats, is the causative agent of hemorrhagic fever with renal syndrome (HFRS) in humans, which has been reported in Asia, Europe, the Americas, and Africa. Defining SEOV genome sequences plays a critical role in development of preventive and therapeutic strategies against the unique worldwide hantavirus. We applied multiplex PCR-based next-generation sequencing to obtain SEOV genome sequences from clinical and reservoir host specimens. Epidemiologic surveillance of R. norvegicus rats in South Korea during 2000-2016 demonstrated that the serologic prevalence of enzootic SEOV infections was not significant on the basis of sex, weight (age), and season. Viral loads of SEOV in rats showed wide dissemination in tissues and dynamic circulation among populations. Phylogenetic analyses showed the global diversity of SEOV and possible genomic configuration of genetic exchanges.

Febrile Respiratory Illness Associated with Human Adenovirus Type 55 in South Korea Military, 2014-2016

An outbreak of febrile respiratory illness associated with human adenovirus (HAdV) occurred in the South Korea military during the 2014-15 influenza season and thereafter. Molecular typing and phylogenetic analysis of patient samples identified HAdV type 55 as the causative agent. Emergence of this novel HAdV necessitates continued surveillance in military and civilian populations.

Construction of Bacillus thuringiensis Simulant Strains Suitable for Environmental Release.

For a surrogate bacterium to be used in outdoor studies, it is important to consider environmental and human safety and ease of detection. Recently, Bacillus thuringiensis, a popular bioinsecticide bacterium, has been gaining attention as a surrogate bacterium for use in biodefense. In this study, we constructed simulant strains of B. thuringiensis with enhanced characteristics for environmental studies. Through transposon mutagenesis, pigment genes were inserted into the chromosome, producing yellow-colored colonies for easy detection. To prevent persistence of spores in the environment, a genetic circuit was designed to produce a spore without sporulation capability. Two loxP sites were inserted, one on each side of the spo0A gene, which encodes a sporulation master regulator, and a sporulation-dependent Cre expression cassette was inserted into the chromosome. This genetic circuit successfully deleted spo0A during sporulation, producing spores that lacked the spo0A gene. In addition, two major α/ß-type small acid-soluble spore protein (SASP) genes, predicted by synteny analysis, were deleted. The spores of the mutant strain showed increased UV-C sensitivity and quickly lost viability when tested in a solar simulator. When the spores of the mutant strain were administered to the lungs of BALB/c mice, cells were quickly removed from the body, suggesting enhanced in vivo safety. All strains constructed in this study contain no antibiotic resistance markers and all heterologous genes were inserted into the chromosome, which are useful features for simulants to be released into the environment.IMPORTANCEB. thuringiensis has recently been receiving increasing attention as a good spore simulant in biodefense research. However, few studies were done to properly address many important features of B. thuringiensis as a simulant in environmental studies. Since spores can persist in the environment for years after release, environmental contamination is a big problem, especially when genetically engineered strains are used. To solve these problems, we report here the development of B. thuringiensis simulant strains that are capable of forming yellow colonies for easy detection, incapable of forming spores more than once due to a genetic circuit, and lacking in two major SASP genes. The genetic circuit to produce a spore without sporulation capability, together with the deletion of SASP genes, ensures the environmental and human safety of the simulant strains developed in this study. All of these features will allow wider use of B. thuringiensis as a simulant for Bacillus anthracis in environmental release studies.

Whole-genome sequence analysis of Zika virus, amplified from urine of traveler from the Philippines.

Zika virus (ZIKV) (genus Flavivirus, family Flaviviridae) is an emerging pathogen associated with microcephaly and Guillain-Barré syndrome. The rapid spread of ZIKV disease in over 60 countries and the large numbers of travel-associated cases have caused worldwide concern. Thus, intensified surveillance of cases among immigrants and tourists from ZIKV-endemic areas is important for disease control and prevention. In this study, using Next Generation Sequencing, we reported the first whole-genome sequence of ZIKV strain AFMC-U, amplified from the urine of a traveler returning to Korea from the Philippines. Phylogenetic analysis showed geographic-specific clustering. Our results underscore the importance of examining urine in the diagnosis of ZIKV infection.

Detection of Hantaan virus RNA from anti-Hantaan virus IgG seronegative rodents in an area of high endemicity in Republic of Korea.

Hantaan virus (HTNV), of the family Bunyaviridae, causes hemorrhagic fever with renal syndrome (HFRS) in humans. Although the majority of epidemiologic studies have found that rodents are seropositive for hantavirus-specific immunoglobulin, the discovery of hantavirus RNA in seronegative hosts has led to an investigation of the presence of HTNV RNA in rodents captured in HFRS endemic areas. HTNV RNA was detected in seven (3.8%) of 186 anti-HTNV IgG seronegative rodents in Republic of Korea (ROK) during 2013-2014. RT-qPCR for HTNV RNA revealed dynamic virus-host interactions of HTNV in areas of high endemicity, providing important insights into the epidemiology of hantaviruses.

Sensing of microbial molecular patterns by Toll-like receptors.

Toll-like receptors (TLRs) sense structural patterns in microbial molecules and initiate immune defense mechanisms. The structures of many extracellular and intracellular domains of TLRs have been studied in the last 10 years. These structures reveal the extraordinary diversity of TLR-ligand interactions. Some TLRs use internal hydrophobic pockets to bind bacterial ligands and others use solvent-exposed surfaces to bind hydrophilic ligands. The structures suggest a common activation mechanism for TLRs: ligand binding to extracellular domains induces dimerization of the intracellular domains and so activates intracellular signaling pathways. Recently, the structure of the death domain complex of one of the signaling adapters, myeloid differentiation factor 88 (MyD88), has been determined. This structure shows how aggregation of signaling adapters recruits downstream kinases. However, we are still far from a complete understanding of TLR activation. We need to study the structures of TLR7-10 in complex with their ligands. We also need to determine the structures of TLR-adapter aggregates to understand activation mechanisms and the specificity of the signaling pathways. Ultimately, we will have to study the structures of the complete TLR signaling complexes containing full-length receptors, ligands, signaling, and bridging adapters, and some of the downstream kinases to understand how TLRs sense microbial infections and activate immune responses against them.

The structural basis of lipopolysaccharide recognition by the TLR4-MD-2 complex.

The lipopolysaccharide (LPS) of Gram negative bacteria is a well-known inducer of the innate immune response. Toll-like receptor (TLR) 4 and myeloid differentiation factor 2 (MD-2) form a heterodimer that recognizes a common 'pattern' in structurally diverse LPS molecules. To understand the ligand specificity and receptor activation mechanism of the TLR4-MD-2-LPS complex we determined its crystal structure. LPS binding induced the formation of an m-shaped receptor multimer composed of two copies of the TLR4-MD-2-LPS complex arranged symmetrically. LPS interacts with a large hydrophobic pocket in MD-2 and directly bridges the two components of the multimer. Five of the six lipid chains of LPS are buried deep inside the pocket and the remaining chain is exposed to the surface of MD-2, forming a hydrophobic interaction with the conserved phenylalanines of TLR4. The F126 loop of MD-2 undergoes localized structural change and supports this core hydrophobic interface by making hydrophilic interactions with TLR4. Comparison with the structures of tetra-acylated antagonists bound to MD-2 indicates that two other lipid chains in LPS displace the phosphorylated glucosamine backbone by approximately 5 A towards the solvent area. This structural shift allows phosphate groups of LPS to contribute to receptor multimerization by forming ionic interactions with a cluster of positively charged residues in TLR4 and MD-2. The TLR4-MD-2-LPS structure illustrates the remarkable versatility of the ligand recognition mechanisms employed by the TLR family, which is essential for defence against diverse microbial infection.

Post-Vaccination Delivery of CpG ODNs Enhances the Th2-Associated Protective Immunity of the Smallpox DNA Vaccine.

Potential threat of smallpox bioterrorism and concerns related to the adverse effects of currently licensed live-virus vaccines suggest the need to develop novel vaccines with better efficacy against smallpox. Use of DNA vaccines containing specific antigen-encoding plasmids prevents the risks associated with live-virus vaccines, offering a promising alternative to conventional smallpox vaccines. In this study, we investigated the efficiency of toll-like receptor (TLR) ligands in enhancing the immunogenicity of smallpox DNA vaccines. BALB/c mice were immunized with a DNA vaccine encoding the vaccinia virus L1R protein, along with the cytosine-phosphate-guanine (CpG) motif as a vaccine adjuvant, and their immune response was analyzed. Administration of B-type CpG oligodeoxynucleotides (ODNs) as TLR9 ligands 24 h after DNA vaccination enhanced the Th2-biased L1R-specific antibody immunity in mice. Moreover, B-type CpG ODNs improved the protective effects of the DNA vaccine against the lethal Orthopoxvirus challenge. Therefore, use of L1R DNA vaccines with CpG ODNs as adjuvants is a promising approach to achieve effective immunogenicity against smallpox infection.

Crystallographic and mutational analysis of the CD40-CD154 complex and its implications for receptor activation.

CD40 is a tumor necrosis factor receptor (TNFR) family protein that plays an important role in B cell development. CD154/CD40L is the physiological ligand of CD40. We have determined the crystal structure of the CD40-CD154 complex at 3.5 Å resolution. The binding site of CD40 is located in a crevice formed between two CD154 subunits. Charge complementarity plays a critical role in the CD40-CD154 interaction. Some of the missense mutations found in hereditary hyper-IgM syndrome can be mapped to the CD40-CD154 interface. The CD40 interaction area of one of the CD154 subunits is twice as large as that of the other subunit forming the binding crevice. This is because cysteine-rich domain 3 (CRD3) of CD40 has a disulfide bridge in an unusual position that alters the direction of the ladder-like structure of CD40. The Ser(132) loop of CD154 is not involved in CD40 binding but its substitution significantly reduces p38- and ERK-dependent signaling by CD40, whereas JNK-dependent signaling is not affected. These findings suggest that ligand-induced di- or trimerization is necessary but not sufficient for complete activation of CD40.

Molecular characterization of SARS-CoV-2 from the saliva of patients in the Republic of Korea in 2020.

Background and aims: Despite global vaccination efforts, the number of confirmed cases of coronavirus disease 2019 (COVID-19) remains high. To overcome the crisis precipitated by the ongoing pandemic, characteristic studies such as virus diagnosis, isolation, and genome analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are necessary. Herein, we report the isolation and molecular characterization of SARS-CoV-2 from the saliva of patients who had tested positive for COVID-19 at Proving Ground in Taean County, Republic of Korea, in 2020. Methods: We analyzed the whole-genome sequence of SARS-CoV-2 isolated from the saliva samples of patients through next-generation sequencing. We also successfully isolated SARS-CoV-2 from the saliva samples of two patients by using cell culture, which was used to study the cytopathic effects and viral replication in Vero E6 cells. Results: Whole-genome sequences of the isolates, SARS-CoV-2 ADD-2 and ADD-4, obtained from saliva were identical, and phylogenetic analysis using Bayesian inference methods showed SARS-CoV-2 GH clade (B.1.497) genome-specific clustering. Typical coronavirus-like particles, with diameters of 70-120 nm, were observed in the SARS-CoV-2 infected Vero E6 cells using transmission electron microscopy. Conclusion: In conclusion, this report provides insights into the molecular diagnosis, isolation, genetic characteristics, and diversity of SARS-CoV-2 isolated from the saliva of patients. Further studies are needed to explore and monitor the evolution and characteristics of SARS-CoV-2 variants.

A Sensitive Immunodetection Assay Using Antibodies Specific to Staphylococcal Enterotoxin B Produced by Baculovirus Expression.

Staphylococcal enterotoxin B (SEB) is a potent bacterial toxin that causes inflammatory stimulation and toxic shock, thus it is necessary to detect SEB in food and environmental samples. Here, we developed a sensitive immunodetection system using monoclonal antibodies (mAbs). Our study is the first to employ a baculovirus expression vector system (BEVS) to produce recombinant wild-type SEB. BEVS facilitated high-quantity and pure SEB production from suspension-cultured insect cells, and the SEB produced was characterized by mass spectrometry analysis. The SEB was stable at 4 °C for at least 2 years, maintaining its purity, and was further utilized for mouse immunization to generate mAbs. An optimal pair of mAbs non-competitive to SEB was selected for sandwich enzyme-linked immunosorbent assay-based immunodetection. The limit of detection of the immunodetection method was 0.38 ng/mL. Moreover, it displayed higher sensitivity in detecting SEB than commercially available immunodetection kits and retained detectability in various matrices and S. aureus culture supernatants. Thus, the results indicate that BEVS is useful for producing pure recombinant SEB with its natural immunogenic property in high yield, and that the developed immunodetection assay is reliable and sensitive for routine identification of SEB in various samples, including foods.

A Portable Diagnostic Assay, Genetic Diversity, and Isolation of Seoul Virus from Rattus norvegicus Collected in Gangwon Province, Republic of Korea.

Seoul virus (SEOV), an etiological agent for hemorrhagic fever with renal syndrome, poses a significant public health threat worldwide. This study evaluated the feasibility of a mobile Biomeme platform for facilitating rapid decision making of SEOV infection. A total of 27 Rattus norvegicus were collected from Seoul Metropolitan City and Gangwon Province in Republic of Korea (ROK), during 2016-2020. The serological and molecular prevalence of SEOV was 5/27 (18.5%) and 2/27 (7.4%), respectively. SEOV RNA was detected in multiple tissues of rodents using the Biomeme device, with differences in Ct values ranging from 0.6 to 2.1 cycles compared to a laboratory benchtop system. Using amplicon-based next-generation sequencing, whole-genome sequences of SEOV were acquired from lung tissues of Rn18-1 and Rn19-5 collected in Gangwon Province. Phylogenetic analysis showed a phylogeographical diversity of rat-borne orthohantavirus collected in Gangwon Province. We report a novel isolate of SEOV Rn19-5 from Gangwon Province. Our findings demonstrated that the Biomeme system can be applied for the molecular diagnosis of SEOV comparably to the laboratory-based platform. Whole-genome sequencing of SEOV revealed the phylogeographical diversity of orthohantavirus in the ROK. This study provides important insights into the field-deployable diagnostic assays and genetic diversity of orthohantaviruses for the rapid response to hantaviral outbreaks in the ROK.

Whole-genome sequencing and genetic diversity of severe fever with thrombocytopenia syndrome virus using multiplex PCR-based nanopore sequencing, Republic of Korea.

BACKGROUND: Whole-genome sequencing plays a critical role in the genomic epidemiology intended to improve understanding the spread of emerging viruses. Dabie bandavirus, causing severe fever with thrombocytopenia syndrome (SFTS), is a zoonotic tick-borne virus that poses a significant public health threat. We aimed to evaluate a novel amplicon-based nanopore sequencing tool to obtain whole-genome sequences of Dabie bandavirus, also known as SFTS virus (SFTSV), and investigate the molecular prevalence in wild ticks, Republic of Korea (ROK). PRINCIPAL FINDINGS: A total of 6,593 ticks were collected from Gyeonggi and Gangwon Provinces, ROK in 2019 and 2020. Quantitative polymerase chain reaction revealed the presence of SFSTV RNA in three Haemaphysalis longicornis ticks. Two SFTSV strains were isolated from H. longicornis captured from Pocheon and Cheorwon. Multiplex polymerase chain reaction-based nanopore sequencing provided nearly full-length tripartite genome sequences of SFTSV within one hour running. Phylogenetic and reassortment analyses were performed to infer evolutionary relationships among SFTSVs. Phylogenetic analysis grouped SFTSV Hl19-31-4 and Hl19-31-13 from Pocheon with sub-genotype B-1 in all segments. SFTSV Hl20-8 was found to be a genomic organization compatible with B-1 (for L segment) and B-2 (for M and S segments) sub-genotypes, indicating a natural reassortment between sub-genotypes. CONCLUSION/SIGNIFICANCE: Amplicon-based next-generation sequencing is a robust tool for whole-genome sequencing of SFTSV using the nanopore platform. The molecular prevalence and geographical distribution of SFTSV enhanced the phylogeographic map at high resolution for sophisticated prevention of emerging SFTS in endemic areas. Our findings provide important insights into the rapid whole-genome sequencing and genetic diversity for the genome-based diagnosis of SFTSV in the endemic outbreak.

Discovery of ß-aminoacyl containing thiazolidine derivatives as potent and selective dipeptidyl peptidase IV inhibitors.

A series of ß-aminoacyl containing thiazolidine derivatives was synthesized and evaluated for their ability to inhibit DPP-IV. Several thiazolidine derivatives with an acid moiety were found to be potent DPP-IV inhibitors. Among them, compound 2da is the most active in this series with an IC(50) value of 1 nM, and it showed excellent selectivity over DPP-IV related enzymes including DPP-2, DPP-8, and DPP-9. Compound 2da is chemically and metabolically stable, and showed no CYP inhibition, hERG binding or cytotoxicity. Compound 2db, an ester prodrug of 2da, showed good in vivo DPP-IV inhibition after oral administration in rat and dog models.

Urinary genome detection and tracking of Hantaan virus from hemorrhagic fever with renal syndrome patients using multiplex PCR-based next-generation sequencing.

BACKGROUND: Hantavirus infection occurs through the inhalation of aerosolized excreta, including urine, feces, and saliva of infected rodents. The presence of Hantaan virus (HTNV) RNA or infectious particles in urine specimens of patient with hemorrhagic fever with renal syndrome (HFRS) remains to be investigated. METHODOLOGY/PRINCIPAL FINDINGS: We collected four urine and serum specimens of Republic of Korea Army (ROKA) patients with HFRS. We performed multiplex PCR-based next-generation sequencing (NGS) to obtain the genome sequences of clinical HTNV in urine specimens containing ultra-low amounts of viral genomes. The epidemiological and phylogenetic analyses of HTNV demonstrated geographically homogenous clustering with those in Apodemus agrarius captured in highly endemic areas, indicating that phylogeographic tracing of HTNV genomes reveals the potential infection sites of patients with HFRS. Genetic exchange analyses showed a genetic configuration compatible with HTNV L segment exchange in nature. CONCLUSION/SIGNIFICANCE: Our results suggest that whole or partial genome sequences of HTNV from the urine enabled to track the putative infection sites of patients with HFRS by phylogeographically linking to the zoonotic HTNV from the reservoir host captured at endemic regions. This report raises awareness among physicians for the presence of HTNV in the urine of patients with HFRS.

Genomic investigation of the coronavirus disease-2019 outbreak in the Republic of Korea.

The South Korean government effectively contained the coronavirus disease-2019 (COVID-19) outbreak primarily associated with a religious group. We conducted SARS-CoV-2 whole genome sequencing of 66 cases to investigate connections among the initial South Korean cases and the religious group outbreak. We assessed the accuracy of genomic investigation by comparing the whole genome sequences with comprehensive contact tracing records. Five transmission clusters were estimated among the 15 initial cases. The six close-contact cases and two potential exposure pairs identified by contact tracing showed two or fewer nucleotide base differences. Additionally, we identified two transmission clusters that were phylogenetically distinct from the initial clusters, sharing common G11083T, G26144T, and C14805T markers. The strain closest to the two additional clusters was identified from a pair of identical sequences isolated from individuals who traveled from Wuhan to Italy. Our findings provide insights into the origins of community spread of COVID-19.