Total infectome characterization of respiratory infections in pre-COVID-19 Wuhan,China

At the end of 2019 Wuhan witnessed an outbreak of “atypical pneumonia” that later developed into a global pandemic. Metagenomic sequencing rapidly revealed the causative agent of this outbreak to be a novel coronavirus denoted SARS-CoV-2. To provide a snapshot of the pathogens in pneumonia-associated respiratory samples from Wuhan prior to the emergence of SARS-CoV-2, we collected bronchoalveolar lavage fluid samples from 408 patients presenting with pneumonia and acute respiratory infections at the Central Hospital of Wuhan between 2016 and 2017. Unbiased total RNA sequencing was performed to reveal their “total infectome”, including viruses, bacteria and fungi. We identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections (13.5%). SARS-CoV-2 was not present. These data depict a stable core infectome comprising common respiratory pathogens such as rhinoviruses and influenza viruses, an atypical respiratory virus (EV-D68), and a single case of a sporadic zoonotic pathogen–Chlamydia psittaci. Samples from patients experiencing respiratory disease on average had higher pathogen abundance than healthy controls. Phylogenetic analyses of individual pathogens revealed multiple origins and global transmission histories, highlighting the connectedness of the Wuhan population. This study provides a comprehensive overview of the pathogens associated with acute respiratory infections and pneumonia, which were more diverse and complex than obtained using targeted PCR or qPCR approaches. These data also suggest that SARS-CoV-2 or closely related viruses were absent from Wuhan in 2016–2017.     

Isolation and Growth Characteristics of SARS-CoV-2 in Vero Cell

The coronavirus disease 2019 (COVID-19) broke out in early December 2019 in Wuhan, China and escalated into a global pandemic. There is an urgent need to understand the biology of SARS-CoV-2. In this letter, we report the isolation and characterization of seven isolates of SARS-CoV-2. Results show that our viruses have 99% sequence identity with published virus sequences. In addition, all viruses grew well in Vero cells, and one of the viruses had a deletion mutation after short passage. These results shall facilitate the understanding of the characteristics of SARS-CoV-2 in vitro.     

Contrasting Patterns in the Early Stage of SARS-CoV-2 Evolution between Humans and Minks

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor-binding domain and receptor-binding motif. An excess of high-frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low-frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies that the virus may not be preadapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.     

Host-Specific and Segment-Specific Evolutionary Dynamics of Avian and Human Influenza A Viruses: A Systematic Review

Understanding the evolutionary dynamics of influenza viruses is essential to control both avian and human influenza. Here, we analyze host-specific and segment-specific Tajima’s D trends of influenza A virus through a systematic review using viral sequences registered in the National Center for Biotechnology Information. To avoid bias from viral population subdivision, viral sequences were stratified according to their sampling locations and sampling years. As a result, we obtained a total of 580 datasets each of which consists of nucleotide sequences of influenza A viruses isolated from a single population of hosts at a single sampling site within a single year. By analyzing nucleotide sequences in the datasets, we found that Tajima’s D values of viral sequences were different depending on hosts and gene segments. Tajima’s D values of viruses isolated from chicken and human samples showed negative, suggesting purifying selection or a rapid population growth of the viruses. The negative Tajima’s D values in rapidly growing viral population were also observed in computer simulations. Tajima’s D values of PB2, PB1, PA, NP, and M genes of the viruses circulating in wild mallards were close to zero, suggesting that these genes have undergone neutral selection in constant-sized population. On the other hand, Tajima’s D values of HA and NA genes of these viruses were positive, indicating HA and NA have undergone balancing selection in wild mallards. Taken together, these results indicated the existence of unknown factors that maintain viral subtypes in wild mallards.     

Genomic analysis of SARS-CoV-2 reveals local viral evolution in Ghana

The confirmed case fatality rate for the coronavirus disease 2019 (COVID-19) in Ghana has dropped from a peak of 2% in March to be consistently below 1% since May 2020. Globally, case fatality rates have been linked to the strains/clades of circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) within a specific country. Here we present 46 whole genomes of SARS-CoV-2 circulating in Ghana, from two separate sequencing batches: 15 isolates from the early epidemic (March 12–April 1 2020) and 31 from later time-points ( 25–27 May 2020). Sequencing was carried out on an Illumina MiSeq system following an amplicon-based enrichment for SARS-CoV-2 cDNA. After genome assembly and quality control processes, phylogenetic analysis showed that the first batch of 15 genomes clustered into five clades: 19A, 19B, 20A, 20B, and 20C, whereas the second batch of 31 genomes clustered to only three clades 19B, 20A, and 20B. The imported cases (6/46) mapped to circulating viruses in their countries of origin, namely, India, Hungary, Norway, the United Kingdom, and the United States of America. All genomes mapped to the original Wuhan strain with high similarity (99.5–99.8%). All imported strains mapped to the European superclade A, whereas 5/9 locally infected individuals harbored the B4 clade, from the East Asian superclade B. Ghana appears to have 19B and 20B as the two largest circulating clades based on our sequence analyses. In line with global reports, the D614G linked viruses seem to be predominating. Comparison of Ghanaian SARS-CoV-2 genomes with global genomes indicates that Ghanaian strains have not diverged significantly from circulating strains commonly imported into Africa. The low level of diversity in our genomes may indicate lower levels of transmission, even for D614G viruses, which is consistent with the relatively low levels of infection reported in Ghana.     

SARS-CoV-2 Serological Survey of Cats in China before and after the Pandemic

The high percentage of seropositivity of cats in Wuhan, the epicentre of the SARS-CoV-2 epidemic, could be due to the large number of infected human cases (more than 14 folds than that in other cities) (Table 1) where cats may have been more frequently exposed to SARS-CoV-2 patients or contaminated environment. As yet, SARS-CoV-2 serological prevalence of cats in other Chinese cities remains unknown. Therefore, a serological survey including more cities with numbers of SARS-CoV-2 human cases in China will be valuable for elucidating the role of cats in transmission of the viruses and relieve public concerns. In this study, 630 cat serum samples collected before November 2019 and 423 cat serum samples collected during SARS-CoV-2 outbreak (from February 2020 to April 2020) in 20 cities in China for detecting the prevalence of SARS-CoV-2 specific antibodies.     

Characteristics of SARS-CoV-2 transmission in a medium-sized city with traditional communities during the early COVID-19 epidemic in China

The nationwide COVID-19 epidemic ended in 2020, a few months after its outbreak in Wuhan, China at the end of 2019. Most COVID-19 cases occurred in Hubei Province, with a few local outbreaks in other provinces of China. A few studies have reported the early SARS-CoV-2 epidemics in several large cities or provinces of China. However, information regarding the early epidemics in small and medium-sized cities, where there are still traditionally large families and community culture is more strongly maintained and thus, transmission profiles may differ, is limited. In this study, we characterized 60 newly sequenced SARS-CoV-2 genomes from Anyang as a representative of small and medium-sized Chinese cities, compared them with more than 400 reference genomes from the early outbreak, and studied the SARS-CoV-2 transmission profiles. Genomic epidemiology revealed multiple SARS-CoV-2 introductions in Anyang and a large-scale expansion of the epidemic because of the large family size. Moreover, our study revealed two transmission patterns in a single outbreak, which were attributed to different social activities. We observed the complete dynamic process of single-nucleotide polymorphism development during community transmission and found that intrahost variant analysis was an effective approach to studying cluster infections. In summary, our study provided new SARS-CoV-2 transmission profiles representative of small and medium-sized Chinese cities as well as information on the evolution of SARS-CoV-2 strains during the early COVID-19 epidemic in China.     

Serosurvey for SARS-CoV-2 among blood donors in Wuhan,China from September to December 2019

The emerging of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused COVID-19 pandemic. The first case of COVID-19 was reported at early December in 2019 in Wuhan City, China. To examine specific antibodies against SARS-CoV-2 in biological samples before December 2019 would give clues when the epidemic of SARS-CoV-2 might start to circulate in populations. We obtained all 88,517 plasmas from 76,844 blood donors in Wuhan between 1 September and 31 December 2019. We first evaluated the pan-immunoglobin (pan-Ig) against SARS-CoV-2 in 43,850 samples from 32,484 blood donors with suitable sample quality and enough volume. Two hundred and sixty-four samples from 213 donors were pan-Ig reactive, then further tested IgG and IgM, and validated by neutralizing antibodies against SARS-CoV-2. Two hundred and thirteen samples (from 175 donors) were only pan-Ig reactive, 8 (from 4 donors) were pan-Ig and IgG reactive, and 43 (from 34 donors) were pan-Ig and IgM reactive. Microneutralization assay showed all negative results. In addition, 213 screened reactive donors were analyzed and did not show obviously temporal or regional tendency, but the distribution of age showed a difference compared with all tested donors. Then we reviewed SARS-CoV-2 antibody results from these donors who donated several times from September 2019 to June 2020, partly tested in a previous published study, no one was found a significant increase in S/CO of antibodies against SARS-CoV-2. Our findings showed no SARS-CoV-2-specific antibodies existing among blood donors in Wuhan, China before 2020, indicating no evidence of transmission of COVID-19 before December 2019 in Wuhan, China.     

TMPRSS2 promotes SARS-CoV-2 evasion from NCOA7-mediated restriction

Interferons play a critical role in regulating host immune responses to SARS-CoV-2, but the interferon (IFN)-stimulated gene (ISG) effectors that inhibit SARS-CoV-2 are not well characterized. The IFN-inducible short isoform of human nuclear receptor coactivator 7 (NCOA7) inhibits endocytic virus entry, interacts with the vacuolar ATPase, and promotes endo-lysosomal vesicle acidification and lysosomal protease activity. Here, we used ectopic expression and gene knockout to demonstrate that NCOA7 inhibits infection by SARS-CoV-2 as well as by lentivirus particles pseudotyped with SARS-CoV-2 Spike in lung epithelial cells. Infection with the highly pathogenic, SARS-CoV-1 and MERS-CoV, or seasonal, HCoV-229E and HCoV-NL63, coronavirus Spike-pseudotyped viruses was also inhibited by NCOA7. Importantly, either overexpression of TMPRSS2, which promotes plasma membrane fusion versus endosomal fusion of SARS-CoV-2, or removal of Spike’s polybasic furin cleavage site rendered SARS-CoV-2 less sensitive to NCOA7 restriction. Collectively, our data indicate that furin cleavage sensitizes SARS-CoV-2 Spike to the antiviral consequences of endosomal acidification by NCOA7, and suggest that the acquisition of furin cleavage may have favoured the co-option of cell surface TMPRSS proteases as a strategy to evade the suppressive effects of IFN-induced endo-lysosomal dysregulation on virus infection.     

Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen

Virus host shifts are generally associated with novel adaptations to exploit the cells of the new host species optimally. Surprisingly, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has apparently required little to no significant adaptation to humans since the start of the Coronavirus Disease 2019 (COVID-19) pandemic and to October 2020. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus the early SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects evidence for significant positive episodic diversifying selection acting at the base of the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in these ancestral bat hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor about 1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. While an undiscovered “facilitating” intermediate species cannot be discounted, collectively, our results support the progenitor of SARS-CoV-2 being capable of efficient human–human transmission as a consequence of its adaptive evolutionary history in bats, not humans, which created a relatively generalist virus.

A study of the natural origins of SARS-CoV-2 reveals very little adaptive evolution occurring since it emerged in humans, but strong evolutionary signals in the bat virus lineage from which SARS-CoV-2 arose. Evolution in bats involved lineage-specific depletion of CpG nucleotides (linked to host anti-viral molecules), and clear evidence of recombination across these lineages, supporting bat host species’ influence on the ancestral viruses.     

The SARS-CoV-2 Delta variant induces an antibody response largely focused on class 1 and 2 antibody epitopes

Exposure histories to SARS-CoV-2 variants and vaccinations will shape the specificity of antibody responses. To understand the specificity of Delta-elicited antibody immunity, we characterize the polyclonal antibody response elicited by primary or mRNA vaccine-breakthrough Delta infections. Both types of infection elicit a neutralizing antibody response focused heavily on the receptor-binding domain (RBD). We use deep mutational scanning to show that mutations to the RBD’s class 1 and class 2 epitopes, including sites 417, 478, and 484–486 often reduce binding of these Delta-elicited antibodies. The anti-Delta antibody response is more similar to that elicited by early 2020 viruses than the Beta variant, with mutations to the class 1 and 2, but not class 3 epitopes, having the largest effects on polyclonal antibody binding. In addition, mutations to the class 1 epitope (e.g., K417N) tend to have larger effects on antibody binding and neutralization in the Delta spike than in the D614G spike, both for vaccine- and Delta-infection-elicited antibodies. These results help elucidate how the antigenic impacts of SARS-CoV-2 mutations depend on exposure history.     

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203–205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral ‘RG’ motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2’s continued adaptation to human infection.     

SARS-CoV-2 suppresses IFNβ production mediated by NSP1, 5, 6, 15, ORF6 and ORF7b but does not suppress the effects of added interferon

Type I Interferons (IFN-Is) are a family of cytokines which play a major role in inhibiting viral infection. Resultantly, many viruses have evolved mechanisms in which to evade the IFN-I response. Here we tested the impact of expression of 27 different SARS-CoV-2 genes in relation to their effect on IFN production and activity using three independent experimental methods. We identified six gene products; NSP6, ORF6, ORF7b, NSP1, NSP5 and NSP15, which strongly (>10-fold) blocked MAVS-induced (but not TRIF-induced) IFNβ production. Expression of the first three of these SARS-CoV-2 genes specifically blocked MAVS-induced IFNβ-promoter activity, whereas all six genes induced a collapse in IFNβ mRNA levels, corresponding with suppressed IFNβ protein secretion. Five of these six genes furthermore suppressed MAVS-induced activation of IFNλs, however with no effect on IFNα or IFNγ production. In sharp contrast, SARS-CoV-2 infected cells remained extremely sensitive to anti-viral activity exerted by added IFN-Is. None of the SARS-CoV-2 genes were able to block IFN-I signaling, as demonstrated by robust activation of Interferon Stimulated Genes (ISGs) by added interferon. This, despite the reduced levels of STAT1 and phospho-STAT1, was likely caused by broad translation inhibition mediated by NSP1. Finally, we found that a truncated ORF7b variant that has arisen from a mutant SARS-CoV-2 strain harboring a 382-nucleotide deletion associating with mild disease (Δ382 strain identified in Singapore & Taiwan in 2020) lost its ability to suppress type I and type III IFN production. In summary, our findings support a multi-gene process in which SARS-CoV-2 blocks IFN-production, with ORF7b as a major player, presumably facilitating evasion of host detection during early infection. However, SARS-CoV-2 fails to suppress IFN-I signaling thus providing an opportunity to exploit IFN-Is as potential therapeutic antiviral drugs.     

The Human Cytomegalovirus UL133-138 Gene Locus Attenuates the Lytic Viral Cycle in Fibroblasts

The genomes of HCMV clinical strains (e.g. FIX, TR, PH, etc) contain a 15 kb region that encodes 20 putative ORFs. The region, termed ULb’, is lost after serial passage of virus in human foreskin fibroblast (HFF) cell culture. Compared to clinical strains, laboratory strains replicate faster and to higher titers of infectious virus. We made recombinant viruses with 22, 14, or 7 ORFs deleted from the ULb’ region using FIX and TR as model clinical strains. We also introduced a stop codon into single ORFs between UL133 and UL138 to prevent protein expression. All deletions within ULb’ and all stop codon mutants within the UL133 to UL138 region increased to varying degrees, viral major immediate early RNA and protein, DNA, and cell-free infectious virus compared to the wild type viruses. The wild type viral proteins slowed down the viral replication process along with cell-free infectious virus release from human fibroblast cells.     

Signal-regulatory protein alpha is an anti-viral entry factor targeting viruses using endocytic pathways

Signal-regulatory protein alpha (SIRPA) is a well-known inhibitor of phagocytosis when it complexes with CD47 expressed on target cells. Here we show that SIRPA decreased in vitro infection by a number of pathogenic viruses, including New World and Old World arenaviruses, Zika virus, vesicular stomatitis virus and pseudoviruses bearing the Machupo virus, Ebola virus and SARS-CoV-2 glycoproteins, but not HSV-1, MLV or mNoV. Moreover, mice with targeted mutation of the Sirpa gene that renders it non-functional were more susceptible to infection with the New World arenaviruses Junín virus vaccine strain Candid 1 and Tacaribe virus, but not MLV or mNoV. All SIRPA-inhibited viruses have in common the requirement for trafficking to a low pH endosomal compartment. This was clearly demonstrated with SARS-CoV-2 pseudovirus, which was only inhibited by SIRPA in cells in which it required trafficking to the endosome. Similar to its role in phagocytosis inhibition, SIRPA decreased virus internalization but not binding to cell surface receptors. We also found that increasing SIRPA levels via treatment with IL-4 led to even greater anti-viral activity. These data suggest that enhancing SIRPA’s activity could be a target for anti-viral therapies.     

Extensive C->U transition biases in the genomes of a wide range of mammalian RNA viruses; potential associations with transcriptional mutations,damage- or host-mediated editing of viral RNA

The rapid evolution of RNA viruses has been long considered to result from a combination of high copying error frequencies during RNA replication, short generation times and the consequent extensive fixation of neutral or adaptive changes over short periods. While both the identities and sites of mutations are typically modelled as being random, recent investigations of sequence diversity of SARS coronavirus 2 (SARS-CoV-2) have identified a preponderance of C->U transitions, proposed to be driven by an APOBEC-like RNA editing process. The current study investigated whether this phenomenon could be observed in datasets of other RNA viruses. Using a 5% divergence filter to infer directionality, 18 from 36 datasets of aligned coding region sequences from a diverse range of mammalian RNA viruses (including Picornaviridae, Flaviviridae, Matonaviridae, Caliciviridae and Coronaviridae) showed a >2-fold base composition normalised excess of C->U transitions compared to U->C (range 2.1x–7.5x), with a consistently observed favoured 5’ U upstream context. The presence of genome scale RNA secondary structure (GORS) was the only other genomic or structural parameter significantly associated with C->U/U->C transition asymmetries by multivariable analysis (ANOVA), potentially reflecting RNA structure dependence of sites targeted for C->U mutations. Using the association index metric, C->U changes were specifically over-represented at phylogenetically uninformative sites, potentially paralleling extensive homoplasy of this transition reported in SARS-CoV-2. Although mechanisms remain to be functionally characterised, excess C->U substitutions accounted for 11–14% of standing sequence variability of structured viruses and may therefore represent a potent driver of their sequence diversification and longer-term evolution.     

Modeling the onset of symptoms of COVID-19: Effects of SARS-CoV-2 variant

Identifying order of symptom onset of infectious diseases might aid in differentiating symptomatic infections earlier in a population thereby enabling non-pharmaceutical interventions and reducing disease spread. Previously, we developed a mathematical model predicting the order of symptoms based on data from the initial outbreak of SARS-CoV-2 in China using symptom occurrence at diagnosis and found that the order of COVID-19 symptoms differed from that of other infectious diseases including influenza. Whether this order of COVID-19 symptoms holds in the USA under changing conditions is unclear. Here, we use modeling to predict the order of symptoms using data from both the initial outbreaks in China and in the USA. Whereas patients in China were more likely to have fever before cough and then nausea/vomiting before diarrhea, patients in the USA were more likely to have cough before fever and then diarrhea before nausea/vomiting. Given that the D614G SARS-CoV-2 variant that rapidly spread from Europe to predominate in the USA during the first wave of the outbreak was not present in the initial China outbreak, we hypothesized that this mutation might affect symptom order. Supporting this notion, we found that as SARS-CoV-2 in Japan shifted from the original Wuhan reference strain to the D614G variant, symptom order shifted to the USA pattern. Google Trends analyses supported these findings, while weather, age, and comorbidities did not affect our model’s predictions of symptom order. These findings indicate that symptom order can change with mutation in viral disease and raise the possibility that D614G variant is more transmissible because infected people are more likely to cough in public before being incapacitated with fever.     

Phylogenetic and Evolutionary History of Influenza B Viruses,which Caused a Large Epidemic in 2011–2012, Taiwan

The annual recurrence of the influenza epidemic is considered to be primarily associated with immune escape due to changes to the virus. In 2011–2012, the influenza B epidemic in Taiwan was unusually large, and influenza B was predominant for a long time. To investigate the genetic dynamics of influenza B viruses during the 2011–2012 epidemic, we analyzed the sequences of 4,386 influenza B viruses collected in Taiwan from 2004 to 2012. The data provided detailed insight into the flux patterns of multiple genotypes. We found that a re-emergent TW08-I virus, which was the major genotype and had co-circulated with the two others, TW08-II and TW08-III, from 2007 to 2009 in Taiwan, successively overtook TW08-II in March and then underwent a lineage switch in July 2011. This lineage switch was followed by the large epidemic in Taiwan. The whole-genome compositions and phylogenetic relationships of the representative viruses of various genotypes were compared to determine the viral evolutionary histories. We demonstrated that the large influenza B epidemic of 2011–2012 was caused by Yamagata lineage TW08-I viruses that were derived from TW04-II viruses in 2004–2005 through genetic drifts without detectable reassortments. The TW08-I viruses isolated in both 2011–2012 and 2007–2009 were antigenically similar, indicating that an influenza B virus have persisted for 5 years in antigenic stasis before causing a large epidemic. The results suggest that in addition to the emergence of new variants with mutations or reassortments, other factors, including the interference of multi-types or lineages of influenza viruses and the accumulation of susceptible hosts, can also affect the scale and time of an influenza B epidemic.     

SARS-CoV-2 lineage B.6 was the major contributor to early pandemic transmission in Malaysia

Malaysia had 10,219 confirmed cases of COVID-19 as of September 20, 2020. About 33% were associated with a Tablighi Jamaat religious mass gathering held in Kuala Lumpur between February 27 and March 3, 2020, which drove community transmission during Malaysia’s second wave. We analysed genome sequences of SARS-CoV-2 from Malaysia to better understand the molecular epidemiology and spread. We obtained 58 SARS-CoV-2 whole genome sequences from patients in Kuala Lumpur and performed phylogenetic analyses on these and a further 57 Malaysian sequences available in the GISAID database. Nine different SARS-CoV-2 lineages (A, B, B.1, B.1.1, B.1.1.1, B.1.36, B.2, B.3 and B.6) were detected in Malaysia. The B.6 lineage was first reported a week after the Tablighi mass gathering and became predominant (65.2%) despite being relatively rare (1.4%) globally. Direct epidemiological links between lineage B.6 viruses and the mass gathering were identified. Increases in reported total cases, Tablighi-associated cases, and community-acquired B.6 lineage strains were temporally linked. Non-B.6 lineages were mainly travel-associated and showed limited onward transmission. There were also temporally correlated increases in B.6 sequences in other Southeast Asian countries, India and Australia, linked to participants returning from this event. Over 95% of global B.6 sequences originated from Asia Pacific. We also report a nsp3-C6310A substitution found in 47.3% of global B.6 sequences which was associated with reduced sensitivity using a commercial diagnostic real-time PCR assay. Lineage B.6 became the predominant cause of community transmission in Malaysia after likely introduction during a religious mass gathering. This event also contributed to spikes of lineage B.6 in other countries in the Asia-Pacific. Mass gatherings can be significant causes of local and global spread of COVID-19. Shared genomic surveillance can be used to identify SARS-CoV-2 transmission chains to aid prevention and control, and to monitor diagnostic molecular assays.Clinical Trial Registration: COVID-19 paper.