Investigation of One Familial Cluster of COVID-19 in Taiwan: Differentiation of Genetic Variation Among Isolates and Implications for Epidemiological Investigation and Surveillance by Genomic Assay
Authors Jian MJ, Chung HY, Chang CK, Hsieh SS, Lin JC, Yeh KM, Chen CW, Chang FY, Chiu SK, Hung KS, Liu MT, Yang JR, Perng CL, Shang HS
Received 22 December 2020
Accepted for publication 16 February 2021
Published 11 March 2021 Volume 2021:14 Pages 971—977
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 3
Editor who approved publication: Dr Sahil Khanna
Ming-Jr Jian,1 Hsing-Yi Chung,1 Chih-Kai Chang,1 Shan-Shan Hsieh,1 Jung-Chung Lin,2 Kuo-Ming Yeh,2 Chien-Wen Chen,3 Feng-Yee Chang,2 Sheng-Kang Chiu,2 Kuo-Sheng Hung,4 Ming-Tsan Liu,5 Ji-Rong Yang,5 Cherng-Lih Perng,1 Hung-Sheng Shang1
1Division of Clinical Pathology, Department of Pathology, Tri-Service General Hospital, National Defense Medical Center, Taipei, Taiwan, Republic of China; 2Division of Infectious Diseases and Tropical Medicine, Department of Medicine, Tri-Service General Hospital, National Defense Medical Center, Taipei, Taiwan, Republic of China; 3Division of Pulmonary and Critical Care Medicine, Department of Medicine, Tri-Service General Hospital, National Defense Medical Center, Taipei, Taiwan, Republic of China; 4Center for Precision Medicine and Genomics, Tri-Service General Hospital, National Defense Medical Center, Taipei, Taiwan, Republic of China; 5Centers for Disease Control, Taipei, Taiwan, Republic of China
Correspondence: Hung-Sheng Shang
Division of Clinical Pathology, Department of Pathology, Tri-Service General Hospital, National Defense Medical Center, 325, Section 2, Cheng-Kung Road, Neihu, Taipei, 114, Taiwan, Republic of China
Email [email protected]
Abstract: The COVID-19 pandemic has caused a global public health crisis. Taiwan experienced two waves of imported cases of coronavirus disease 2019 (COVID-19), first from China in January to late February, 2020 then from other countries starting in early March. As of Dec 14, 2020, 733 cases have been reported in Taiwan, with cases of entire families being infected. This study aimed to investigate the clinical characteristics and differentiation of genetic variation among isolates from a cluster of familial COVID-19 infection. The parents had pneumonia (Case 14, father, and Case 15, mother), the elder son (Case 17) had mild cough, and the younger son (Case 18) was asymptomatic. In this study, four full viral genomes were sequenced by Illumina sequencing directly from specimens. Phylogenetic tree analysis revealed that these sequences came from Italy, not China, indicating that no major strain has been circulating in Taiwan. Several novel mutations were observed in the asymptomatic patient, such as nsp2, nsp12, and nsp14. These mutations may be associated with the severity of COVID-19 infection.
Keywords: COVID-19, SARS-CoV-2, imported family infections, whole-genome sequencing, phylogenetic analysis
Coronavirus disease (COVID-19) has spread globally, mainly via person–person transmission; it poses a major public health concern.1 An unknown etiology of novel coronavirus (SARS-CoV-2) infection was reported in December 2019.2 The disease caused by this virus has been called novel coronavirus disease 2019 (COVID-19). On March 11, 2020, the World Health Organization declared the worldwide spread of COVID-19 a pandemic.3 As of Dec 14, 733 cases have been reported in Taiwan. Here, we report the first four-member familial cluster of COVID-19 in Taiwan. Two patients presented with pneumonia following admission to our hospital (Tri-Service General Hospital). The other two patients were admitted based on positive results of RT-PCR tests for SARS-CoV-2 viral nucleic acid. All patients were successfully treated with oseltamivir, piperacillin/tazobactam, and clarithromycin as empirical therapy, and were eventually discharged. Previous studies have identified positive COVID-19 cases via nucleic acid tests, where some patients remained afebrile and asymptomatic.4,5 However, the viral shedding pattern and genomic characteristics of COVID-19 remain poorly understood. Up-to-date information on viral genomics is crucial in understanding the global dispersion characteristics of the virus and providing insights into its viral pathogenicity and transmission. In this study, we performed whole-genome sequencing of four clinical specimens of SARS-CoV-2 taken from the same infected family. We compared the genomes obtained with those of other strains from the GISAID database to understand their evolutionary trajectory. We also assessed the number of nucleotide substitutions in one family member who was asymptomatic. Through genome-wide association analysis on SARS-CoV-2 genomes, we sought to identify genetic variations that might be associated with the severity of COVID-19 infection. SARS-CoV-2 was confirmed by RT-PCR as described previously6 and re-confirmed by the Taiwan Centers for Disease Control. Extracted viral RNA was stored at −80°C for further analysis. Data such as demographic information, underlying illnesses, clinical manifestations, and outcomes were collected from medical records. The case definition for COVID-19 was in accordance with Taiwan Centers for Disease Control’s recommendations.
Confirmed case number information came from the content of the Taiwan’s Central Epidemic Command Center (CECC) press conference and a press release from the Taiwan Centers for Disease Control . A 53-year-old Taiwanese man (Case 14), his wife (Case 15, 53 years old), and two sons (Cases 17 and 18, 24 and 22 years old, respectively) departed by air on January 22, 2020 from Taiwan, transferred at Hong Kong (where they stayed at the airport for 1 h), and arrived in Italy on January 23; they returned to Taiwan on February 1 (Figure 1). Case 17 had a mild cough since January 27, and case 18 reported no discomfort. Cases 17 and 18 were quarantined, admitted to the hospital on February 7 and 8, and discharged from the hospital on February 27 and March 6, respectively. Case 18 was the first asymptomatic diagnosis in Taiwan and the viral load detected in the body was high. Although the younger son had no symptoms, the other three family members were diagnosed successively (laboratory findings shown in Table 1). Exposure to the same source of infection could have been from contact with another infected person while on an airplane, in an airport, or in Italy.
Table 1 Clinical and Laboratory Characteristics of the Familial Cluster Infection Cases
Figure 1 Genogram of family members infected with COVID-19 who returned to Taiwan from Italy, including their clinical courses.
Whole-Genome Sequencing of SARS-CoV-2
Whole genome sequences of SARS-CoV-2 were obtained using the Illumina TruSeq Stranded mRNA Library Prep Kit protocol to enrich SARS-CoV-2 cDNA using multiplex RT-PCR amplicons. The NovaSeq 6000 platform (Illumina, San Diego, USA) was used for sequencing of the indexed libraries. To determine the full viral genome sequences, the Next Generation Sequencing reads were first mapped to the human genome (hg38), and unmapped reads were compared to the SARS-CoV-2 Wuhan-Hu-1 reference genome sequence (29.9kb ss-RNA; GenBank ID: NC_045512.2) using SPAdes v3.13.2 with the parameter “–trusted-contigs” to assemble the specimen-specific SARS-CoV- 2 genome sequence. All viral assemblies were uploaded to GISAID (Supplementary Table S1).
Phylogenetic Relationship Analysis
SARS-CoV-2 genome sequences from GISAID and those of the four cases were aligned using Clustal Omega.7 Aligned nucleic acid sequences were used to construct phylogenetic trees based on the neighbor-joining tree algorithms provided by PHYLIP.8 The tree file created by PHYLIP was drawn and managed in MEGA X.9 Different COVID-19 sample collection times and phylogenetic trees created by viral genome sequencing with branch lengths were analyzed using the RelTime with Dated Tips (RTDT) method10,11 supported by MEGA X.9 Bootstrap analysis was calculated with 1000 replicates and analysis values were shown next to the branch nodes.12
Four full viral genomes were resolved in this study directly from specimens. Phylogenetic tree analysis showed that these four cases were imported from Italy, not China (Figure 2), indicating that the outbreak of COVID-19 in Italy may have occurred earlier than the end of January, 2020. One of the sequences (case 14, strain name: TSGH-35) lies near the Wuhan SARS-CoV-2 reference sequence representing the early virus strains as an index patient, also mapping the divergence time of late-December 2019 to mid-January 2020 as the origin of the Wuhan strain (Wuhan-Hu-1). To decipher the time point when the family could have been infected, another 75 strains were aligned; the closest area found was Wuhan (Figure 3), indicating that the early Wuhan virus strain spread to Italy. None of the four genomes harbored the D614G mutation, which has been linked to higher mortality in COVID-19.13 A previous study reported that COVID-19 cases were more likely to have comorbidities than non-severe cases, with host contributing factors including hypertension, diabetes mellitus, etc.14 Davies et al also showed that the age disparities in observed cases could be explained by children having lower susceptibility to infection and a lower propensity to showing clinical symptoms.15 This may explain why Case 17 and 18 in our study were paucisymptomatic or asymptomatic. We further investigated the asymptomatic case, in which some novel nucleotide substitutions were observed (Figure 4), such as nsp12 (RdRp) with the P378L mutation, nsp14 with the A23S mutation, nsp2 with the M141L mutation, and N protein with the P309L and F363V mutations. It has been suggested that the number of asymptomatic infections of COVID-19 has been underestimated due to an increasing number of infected people being found who have not traveled to epidemic hotspots or been linked to known COVID-19 patients.16 Previous studies have identified COVID-19 transmission caused by asymptomatic carriers who had normal chest CT findings, but were found to have a viral load similar to that of a symptomatic patient.17,18 RNA viruses show variation in their genomes due to nucleotide substitutions generated by the low fidelity of RNA-dependent RNA polymerase during replication. Further studies are required to investigate the mutation sites correlated with the function of SARS-CoV-2 in asymptomatic cases.
Figure 3 Phylogenetic tree with 75 viral genome strains linking the four familial sequences.
Figure 4 Novel nucleotide substitutions found via whole-genome analysis of the asymptomatic case, Case 18.
Timely sequencing and sharing of whole genomes of SARS-CoV-2 from different locations is important for monitoring genetic changes in the virus that may be associated with viral spreading and clinical manifestations. Here, we determined the sequences of SARS-CoV-2 in one familial cluster of COVID-19 in Taiwan and showed that sequencing is a useful tool for tracing the infection source for this type of RNA virus. Moreover, the existence of fewer nucleotide changes in viral genomes from asymptomatic patients in clustered infections might help us elucidate the viral shedding pattern and replication model with SARS-CoV-2 infection.
This study was approved by the Institutional Review Board, Tri-Service General Hospital (TSGHIRB No.: C202005041), registered on 20 March 2020. Written informed consent was obtained from the participants for the publication of the case report.
Ming-Jr Jian: Concepts, Design, Data analysis, Manuscript preparation; Hsing-Yi Chung, Chih-Kai Chang, Shan-Shan Hsieh, Kuo-Sheng Hung, Cherng-Lih Perng: Data analysis, Definition of intellectual content, Literature search; Jung-Chung Lin, Kuo-Ming Yeh, Chien-Wen Chen, Feng-Yee Chang, Sheng-Kang Chiu, Ming-Tsan Liu, Ji-Rong Yang: Concepts, Design, Definition of intellectual content; Hung-Sheng Shang: Manuscript editing, Manuscript review, Guarantor. All authors made substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; agreed on submission to the current journal; gave final approval of the version to be published; and agreed to be accountable for all aspects of the work.
This study was supported by Tri-Service General Hospital, Taipei, Taiwan, ROC, Grant Numbers: TSGH-D-109142 and TSGH-D-110100. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors report no conflicts of interest in this work.
1. Chan JF-W, Yuan S, Kok K-H, et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet. 2020;395(10223):514–523. doi:10.1016/S0140-6736(20)30154-9
2. Zhu N, Zhang D, Wang W, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. 2020;382(8):727–733. doi:10.1056/NEJMoa2001017
3. Chang FY, Chen HC, Chen PJ, et al. Immunologic aspects of characteristics, diagnosis, and treatment of coronavirus disease 2019 (COVID-19). J Biomed Sci. 2020;27(1):72. doi:10.1186/s12929-020-00663-w
4. Gao Z, Xu Y, Sun C, et al. A systematic review of asymptomatic infections with COVID-19. J Microbiol Immunol Infect. 2020. doi:10.1016/j.jmii.2020.05.001
5. Tabata S, Imai K, Kawano S, et al. Clinical characteristics of COVID-19 in 104 people with SARS-CoV-2 infection on the Diamond Princess cruise ship: a retrospective analysis. Lancet Infect Dis. 2020;20(9):1043–1050. doi:10.1016/S1473-3099(20)30482-5
6. Perng CL, Jian MJ, Chang CK, et al. Novel rapid identification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by real-time RT-PCR using BD Max Open System in Taiwan. PeerJ. 2020;8:e9318. doi:10.7717/peerj.9318
7. Sievers F, Wilm A, Dineen D, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539. doi:10.1038/msb.2011.75
8. Felsenstein J. PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by author. Seattle: Department of Genome Sciences, University of Washington; 2005. Available from: http://evolution.genetics.washington.edu/phylip.html.
9. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–1549. doi:10.1093/molbev/msy096
10. Tamura K, Battistuzzi FU, Billing-Ross P, Murillo O, Filipski A, Kumar S. Estimating divergence times in large molecular phylogenies. Proc Natl Acad Sci U S A. 2012;109(47):19333–19338.
11. Tamura K, Tao Q, Kumar S. Theoretical foundation of the RelTime method for estimating divergence times from variable evolutionary rates. Mol Biol Evol. 2018;35(7):1770–1782.
12. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39(4):783–791. doi:10.2307/2408678
13. Eaaswarkhanth M, Al Madhoun A, Al-Mulla F. Could the D614G substitution in the SARS-CoV-2 spike (S) protein be associated with higher COVID-19 mortality? Int J Infect Dis. 2020;96:459–460. doi:10.1016/j.ijid.2020.05.071
14. Lai -C-C, Liu YH, Wang C-Y, et al. Asymptomatic carrier state, acute respiratory disease, and pneumonia due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2): facts and myths. J Microbiol Immunol Infect. 2020;53(3):404–412. doi:10.1016/j.jmii.2020.02.012
15. Davies NG, Klepac P, Liu Y, et al. Age-dependent effects in the transmission and control of COVID-19 epidemics. Nat Med. 2020;26(8):1205–1211. doi:10.1038/s41591-020-0962-9
16. Qiu J. Covert coronavirus infections could be seeding new outbreaks. Nature. 2020. doi:10.1038/d41586-020-00822-x
17. Bai Y, Yao L, Wei T, et al. Presumed asymptomatic carrier transmission of COVID-19. JAMA. 2020.
18. Zou L, Ruan F, Huang M, et al. SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N Engl J Med. 2020;382(12):1177–1179. doi:10.1056/NEJMc2001737
This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.Download Article [PDF]