Professional Documents
Culture Documents
The 2nd Yan Report
The 2nd Yan Report
Li-Meng Yan (MD, PhD)1, Shu Kang (PhD)1, Jie Guan (PhD)1, Shanchang Hu (PhD)1
1
Rule of Law Society & Rule of Law Foundation, New York, NY, USA.
Correspondence: [email protected]
Abstract
Two possibilities should be considered for the origin of SARS-CoV-2: natural evolution or laboratory
creation. In our earlier report titled “Unusual Features of the SARS-CoV-2 Genome Suggesting
Sophisticated Laboratory Modification Rather Than Natural Evolution and Delineation of Its Probable
Synthetic Route”, we disproved the possibility of SARS-CoV-2 arising naturally through evolution and
instead proved that SARS-CoV-2 must have been a product of laboratory modification. Despite this and
similar efforts, the laboratory creation theory continues to be downplayed or even diminished. This is
fundamentally because the natural origin theory remains supported by several novel coronaviruses
published after the start of the outbreak. These viruses (the RaTG13 bat coronavirus, a series of pangolin
coronaviruses, and the RmYN02 bat coronavirus) reportedly share high sequence homology with SARS-
CoV-2 and have altogether constructed a seemingly plausible pathway for the natural evolution of SARS-
CoV-2. Here, however, we use in-depth analyses of the available data and literature to prove that these
novel animal coronaviruses do not exist in nature and their sequences have been fabricated. In addition,
we also offer our insights on the hypothesis that SARS-CoV-2 may have originated naturally from a
coronavirus that infected the Mojiang miners.
Revelation of these virus fabrications renders the natural origin theory unfounded. It also strengthens
our earlier assertion that SARS-CoV-2 is a product of laboratory modification, which can be created in
approximately six months using a template virus owned by a laboratory of the People’s Liberation Army
(PLA). The fact that data fabrications were used to cover up the true origin of SARS-CoV-2 further
implicates that the laboratory modification here is beyond simple gain-of-function research.
The scale and the coordinated nature of this scientific fraud signifies the degree of corruption in the
fields of academic research and public health. As a result of such corruption, damages have been made
both to the reputation of the scientific community and to the well-being of the global community.
Importantly, while SARS-CoV-2 meets the criteria of a bioweapon specified by the PLA, its impact is
well beyond what is conceived for a typical bioweapon. In addition, records indicate that the unleashing
of this weaponized pathogen should have been intentional rather than accidental. We therefore define
SARS-CoV-2 as an Unrestricted Bioweapon and the current pandemic a result of Unrestricted Biowarfare.
We further suggest that investigations should be carried out on the suspected government and individuals
and the responsible ones be held accountable for this brutal attack on the global community.
1
Introduction
SARS-CoV-2 is a novel coronavirus and the causative agent of the COVID-19 pandemic. Despite its
tremendous impact, the origin of SARS-CoV-2, however, has been a topic of great controversy. In our
first report titled “Unusual Features of the SARS-CoV-2 Genome Suggesting Sophisticated Laboratory
Modification Rather Than Natural Evolution and Delineation of Its Probable Synthetic Route”1, we used
biological evidence and in-depth analyses to show that SARS-CoV-2 must be a laboratory product, which
was created by using a template virus (ZC45/ZXC21) owned by military research laboratories under the
control of the Chinese Communist Party (CCP) government. In addition, resources and expertise are all
in place in the Wuhan Institute of Virology (WIV) and related, other CCP-controlled institutions allowing
the creation of SARS-CoV-2 in approximately six months.
What have not been fully described in our earlier analyses are details of the novel animal coronaviruses
published by the CCP-controlled laboratories after the outbreak1. While no coronaviruses reported prior
to 2020 share more than 90% sequence identity with SARS-CoV-22,3, these recently published, novel
animal coronaviruses (the RaTG13 bat coronavirus4, a series of pangolin coronaviruses5-8, and the
RmYN02 bat coronavirus9) all share over 90% sequence identities with SARS-CoV-2. As a result, these
SARS-CoV-2-like viruses have filled an evolutionary gap and served as the founding evidence for the
theory that SARS-CoV-2 has a natural origin10-12. In this report, we provide genetic and other analyses,
which, when combined with recent findings13-21, prove that these novel animal coronaviruses do not exist
in nature and their genomic sequences are results of fabrication.
1. Evidence proving that the RaTG13 virus is fraudulent and does not exist in nature
On February 3rd, 2020, Dr. Zhengli Shi and colleagues published an article in Nature titled “A
pneumonia outbreak associated with a new coronavirus of probable bat origin” (manuscript submitted on
January 20th)4, which was one of the first publications to identify SARS-CoV-2 as the pathogen causing
the disease now widely known as COVID-19. Also reported in this article was a novel bat coronavirus
named RaTG13, the genomic sequence of which was shown to be 96.2% identical to that of SARS-CoV-
2. The close evolutionary relationship between RaTG13 and SARS-CoV-2 as suggested by the high
sequence identity had led to a conclusion that SARS-CoV-2 has a natural origin. These striking findings
have consequently made this article one of the most cited publications in the currently overwhelmed field
of coronavirus research. Interestingly, an article published by Dr. Yong-Zhen Zhang and colleagues on
the same issue of Nature, which also discovered SARS-CoV-2 as the responsible pathogen for COVID-
19, received much less citations2. This latter article made no mention of RaTG132. Instead, Zhang and
colleagues showed that, evolutionarily, SARS-CoV-2 was closest to two bat coronaviruses, ZC45 and
ZXC21, both of which were discovered and characterized by military research laboratories under the
control of the CCP government3. Immediately after the publication of this article, Dr. Zhang’s laboratory
was shut down by the CCP government with no explanations offered22.
Since its publication4, the RaTG13 virus has served as the founding evidence for the theory that SARS-
CoV-2 must have a natural origin10. However, no live virus or an intact genome of RaTG13 have ever
been isolated or recovered. Therefore, the only proof for the “existence” of RaTG13 in nature is its
genomic sequence published on GenBank.
2
1.1 The sequence of RaTG13 uploaded at GenBank can be fabricated
In order to have the sequence of a viral genome successfully uploaded onto GenBank, submitters have
to provide both the assembled genomic sequence (text only) and raw sequencing reads. The latter is used
for quality control and verification purposes. However, due to the huge amount of work involved in
assembling raw reads into complete genomes, no sufficient curation is in place to ensure the correctness
or truthfulness of the uploaded viral genomes. Therefore, an entry on GenBank, which in this case is
equivalent to the existence of an assembled viral genomic sequence and its associated sequencing reads,
is not a definitive proof that this viral genome is correct or real.
Sequencing of a viral RNA genome requires amplifying segments of it using reverse transcriptase PCR
(RT-PCR) as the first step. The products of the RT-PCR, which are double-stranded DNA, would
subsequently be sent for sequencing. The resulted sequencing reads, each ideally revealing the sequence
of a segment of the genome, are then used to assemble the genome of the virus under study (Figure 1A).
Typically, some segments of the genome may not be covered by the initial round of sequencing. Therefore,
gap filling will be carried out, where these missing segments will be amplified specifically and the DNA
products subsequently sequenced. These steps are repeated until a complete genome can be assembled,
ideally with a proper depth to ensure accuracy.
However, this process leaves room for potential fraud. If one intends to fabricate an RNA viral genome
on GenBank, he or she could do so by following these steps: create its genomic sequence on a computer,
have segments of the genome synthesized based on the sequence, amplify each DNA segment through
PCR, and then send the PCR products (may also be mixed with genetic material derived from the alleged
host of the virus to mimic an authentic sequencing sample) for sequencing (Figure 1B). The resulted raw
sequencing reads would be used, together with the created genomic sequence, for establishing an entry on
GenBank. Once accomplished, this entry would be accepted as the evidence for the natural existence of
the corresponding virus. Clearly, a viral genomic sequence and its GenBank entry can be fabricated if
well-planned.
Figure 1. Illustration of steps involved in the sequencing and assembly of coronavirus genomes. A. The normal
process. B. A possible route of fabricating a viral genome by creating a genomic sequence first and obtaining raw
sequencing reads guided by it. NGS: Next Generation Sequencing.
The complete genomic sequence of RaTG13 was first submitted to GenBank on January 27th, 2020.
The raw sequencing reads were made available on February 13th, 2020 (NCBI SRA: SRP249482).
However, the sequencing data for gap filling, which is indispensable in assembling a complete genome,
3
was only made available on May 19th, 2020 (NCBI SRA: SRX8357956). The timing and the reversed order
of events here are strange and suspicious.
The raw sequencing reads of RaTG13 have multiple abnormal features16,21. Despite the sample being
described as a fecal swab, only 0.7% of the raw sequencing reads are bacterial reads while the bacterial
abundance is typically 70~90% when other fecal swab samples were sequenced16,21. In addition, in the
identifiable region of certain sequencing reads, a vast majority of reads are eukaryotic sequences, which
is also highly unusual in the sequencing of fecal swap-derived samples16. Within these eukaryotic reads,
30% of the sequences are of non-bat origin and instead shown to be from many different types of animals
including fox, flying fox, squirrels, etc. These abnormal features are significant and indicate that the raw
sequencing reads should have been obtained via a route that is different from the normal one (Figure 1).
No independent verification of the RaTG13 sequence seems possible because, according to Dr. Zhengli
Shi, the raw sample has been exhausted and no live virus was ever isolated or recovered. Notably, this
information was known to a core circle of virologists early on and apparently accepted by them. It was
then made public, months later, by Dr. Yanyi Wang, director general of the WIV, in an TV interview on
May 23rd, 202023. Dr. Shi also confirmed this publicly in her email interview with Science in July 202024.
However, judging from Shi’s published protocol25, exhaustion of the fecal swap sample is highly
unlikely. According to this protocol, the fecal swab sample would be mixed with 1 ml of viral transport
medium and the supernatant collected. Every 140 ul of the supernatant would then yield 60 ul of extracted
RNA25. For the subsequent step, RT-PCR, 5 ul of this RNA-containing solution is required per reaction25.
Therefore, from one fecal swab sample, at least 80 RT-PCR reactions could be carried out ([1000/140] x
60/5=86). Such an amount is sufficient to support both the initial round of sequencing and the subsequent
gap filling PCR. It would be sufficient to also allow reasonable attempts to isolate live viruses, although
Dr. Shi claimed that no virus isolation was attempted24.
Therefore. the RaTG13 virus and its published sequence are suspicious and show signs of fabrication.
1.2 Other suspicions associated with RaTG13
RaTG13 was reported by Dr. Zhengli Shi from the WIV4. Dr. Shi is a fellow of the American Academy
of Microbiology and one of the most accomplished Chinese virologists. A peer-reviewed article authored
by her and published on the top journal Nature, therefore, brought a great level of comfort for the
coronavirus research community in accepting RaTG13 as a true, nature-born bat coronavirus. As a result,
RaTG13, upon its timely publication, served as the founding evidence for the natural origin theory of
SARS-CoV-2.
However, as revealed in section 1.1, the reported sequence of RaTG13, which is the only proof of the
virus’ existence in nature, is problematic and shows signs of fabrication.
Intriguingly, despite the pivotal role of RaTG13 in revealing the origin of SARS-CoV-2, the
information provided for its discovery was surprisingly scarce with key points missing (location and date
of sample collection, previous knowledge and publication of this virus, etc):
“We then found that a short region of RNA-dependent RNA polymerase (RdRp) from a bat coronavirus
(BatCoV RaTG13)—which was previously detected in Rhinolophus affinis from Yunnan province—
showed high sequence identity to 2019-nCoV. We carried out full-length sequencing on this RNA sample
(GISAID accession number EPI_ISL_402131). Simplot analysis showed that 2019-nCoV was highly
similar throughout the genome to RaTG13 (Fig. 1c), with an overall genome sequence identity of 96.2%.”4
4
Only in the source section of the NCBI entry for RaTG13 (GenBank accession code: MN996532.1),
one could find that the original sample was a “fecal swab” collected on “July 24th, 2013”. A closer look
at the sequence reveals that RaTG13 shares a 100% nucleotide sequence identity with a bat coronavirus
RaBtCoV/4991 on a short, 440-bp RNA-dependent RNA polymerase gene (RdRp) segment.
RaBtCoV/4991 was discovered by Shi and colleagues and published in 201626. As described in the 2016
publication, only a short 440-bp segment of RdRp of the RaBtCoV/4991 virus was sequenced then. Given
the 100% identity on this short gene segment between RaBtCoV/4991 and RaTG13, the field has
demanded clarification of whether or not these two names refer to the same virus. However, Dr. Shi did
not respond to the request or address this question for months. The answer finally came from Peter Daszak,
president of EcoHealth Alliance and long-term collaborator of Shi, who claimed that RaBtCoV/4991 was
RaTG1327.
RaBtCoV/4991 was discovered in the Yunnan province, China. In 2012, six miners suffered from
severe pneumonia after clearing out bat droppings in a mineshaft in Mojiang, Yunnan, and three of them
died soon afterwards28,29. Although it was initially suspected that a SARS-like bat coronavirus may be
responsible for the deaths, no coronavirus was either isolated or detected from the clinical samples30. Also,
first-hand record indicates failure of biopsy and no attempt of autopsy30, which are the gold standards in
the diagnosis of coronavirus infections30. The pathogen responsible for the miners’ deaths therefore
remained an unsolved case31. (Detailed analyses of the Mojiang Miner Passage hypothesis, which was
based on the miners’ case, are provided in section 1.6.) Despite the failed diagnosis, this unknown
pathogen nonetheless triggered immense interests in the virologists in China. Three independent teams,
including that of Dr. Shi’s, made a total of six visits to this mineshaft26,28,31. The Shi group particularly
looked for the presence of bat coronaviruses by amplifying and then sequencing a 440-bp RdRp segment29,
which is a routine procedure the Shi group follows in their surveillance studies. (As shown in section 2.1
of our first report1, this RdRp segment is also frequently used for phylogenetic analyses and is an attractive
target for antiviral drug discovery, which may have contributed to the design of incorporating a unique
RdRp into the genome of SARS-CoV-2.) Out of the many coronaviruses detected, only RaBtCoV/4991
seemed to belong to the group of SARS-related, lineage B β coronaviruses26.
The reporting of RaTG13 is suspicious in three aspects.
First, the whole genome sequencing of RaBtCoV/4991 should not have been delayed until 2020. Given
the Shi group’s consistent interests in studying SARS-like bat coronaviruses and the fact that
RaBtCoV/4991 is a SARS-like coronavirus with a possible connection to the deaths of the miners, it is
highly unlikely that the Shi group would be content with sequencing only a 440-bp segment of RdRp and
not pursue the sequencing of the receptor-binding motif (RBM)-encoding region of the spike gene. In fact,
sequencing of the spike gene is routinely attempted by the Shi group once the presence of a SARS-like
bat coronavirus is confirmed by the sequencing of the 440-bp RdRp segment25,32, although the success of
such efforts is often hindered by the poor quality of the sample.
As quoted above, in the 2020 Nature publication, Shi and colleagues strongly suggested that the
sequencing of the full genome was done in 2020 after they discovered the resemblance between RaTG13
and SARS-CoV-2 on the short RdRp segment4. This, if true, suggests that the quality of the sample should
not be poor. Therefore, there is no technical obstacle for the whole genome sequencing of RaBtCoV/4991.
Clearly, the perceivable motivation of the Shi group to study this RaBtCoV/4991 virus and the fact that
no genome sequencing of it was done for a period of seven years (2013-2020) are hard to reconcile and
explain.
5
However, an intriguing revelation took place in June 2020. Specifically, filenames of the raw
sequencing reads for RaTG13 uploaded on the database were found, which indicate that these sequencing
experiments were done in 2017 and 201833. Likely responding to this revelation, in her email interview
with Science24, Dr. Shi contradicted her own description in the Nature publication4 and admitted that the
sequencing of the full genome of RaTG13 was done in 2018.
Figure 2. Sequence alignment comparing the RBMs of SARS (top) and RaTG13 (red arrow) to the RBMs of bat
coronaviruses that Dr. Zhengli Shi published in high-profile journals between 2013 and 201825,32,34. Amino acid
residues highlighted by Shi as critical for binding the human ACE2 receptor32 are labeled in red text on top.
Alignment was done using the MultAlin webserver (http://multalin.toulouse.inra.fr/multalin/).
Second, RaTG13 has a remarkable RBM as suggested by its reported sequence, and the Shi group have
no reason to delay its publication until 2020. The most critical segment of a SARS-like β coronavirus is
the RBM in the Spike protein as it is fully responsible for binding the host ACE2 receptor and therefore
determines the virus’ potential in infecting humans. The RBM is also the most variable region because it
is under strong positive selection when the virus jumps over to a new host. Sequence alignment on this
crucial RBM motif reveals that the RaTG13 virus rivals with the most highly regarded bat coronaviruses
in terms of resemblance to SARS (Figure 2). RaTG13’s RBM not only is complete in reference to that of
SARS but also is outstanding in its preservation of five residues perceived by Dr. Shi as key in binding
human ACE2 (hACE2)32 (Figure 2, residues labeled with red texts). At position 472, RaTG13 is the only
bat coronavirus that shares a leucine (L) residue with SARS, while the other four key residues are also
largely conserved between the two viruses. Importantly, similar conservation patterns revealed in related
bat coronaviruses, Rs3367 and SHC014, had led to their publication in Nature in 201332. Furthermore,
viruses with less “attractive” RBM sequences (having large gaps and poor in the preservation of key
residues, bottom half of the sequences in Figure 2) were also published by Dr. Shi in other top virology
journals between 2013 and 201825,34. Therefore, if the genomic sequence of RaTG13 had been available
since 2018, it is unlikely that this virus, which has a possible connection to miners’ deaths in 2012 and
has an alarming SARS-like RBM, would be shelfed for two years without publication. Consistent with
this analysis, a recent study indeed proved that the RBD of RaTG13 (produced via gene synthesis based
on its published sequence) was capable of binding hACE235.
6
Third, no follow-up work on RaTG13 has been reported by the Shi group. Upon obtaining the genomic
sequence of a SARS-like bat coronavirus, the Shi group routinely investigate whether or not the virus is
capable of infecting human cells. This pattern of research activities has been shown repeatedly25,32,36-39.
However, such a pattern is not seen here despite that RaTG13 has an interesting RBM and is allegedly the
closest match evolutionarily to SARS-CoV-2.
Clearly, these three aspects deviate from normal research activities and logical thinking, which are
difficult to reconcile or explain. They should have contributed to the intentional omission of key
information in the reporting of RaTG134.
For publications of biological research, it is unethical for authors to change the name of a previously
published virus without any notice or description. It is also unethical for authors to not cite their own
publication where they had characterized and reported the same virus. The violations here by Shi and
colleagues on the reporting of RaTG13 are especially aggravating as the discovery of RaTG13 was central
to uncovering the origin of SARS-CoV-2. By the time of the publication, SARS-CoV-2 had already led
to many deaths in the city of Wuhan and had shown an alarming potential of causing a pandemic. In her
much-delayed response to Science published on July 31st, 202024, Dr. Shi finally commented on the name
change and stated that changing the name to RaTG13 was meant to better reflect the time and location of
sample collection (TG = Tongguan; 13 = 2013). However, such an intention does not seem to justify why
the previous name of RaBtCoV/4991 was never mentioned in the 2020 article4 and why they did not cite
their own 2016 publication where RaBtCoV/4991 was first reported26. Dr. Shi’s recent clarification did
not alter the fact that they have violated the reporting norms of biological research.
In summary, a range of suspicions were associated with the reporting of RaTG13, including the
violations of scientific publication principles, the inconsistency in the descriptions of the sequencing dates,
and the contradiction between the sequencing of its genome in 2018 and the publication of it in 2020 when
this virus has a striking RBM and a possible connection to pneumonia-associated deaths. Adding to these
suspicions are the exquisite timing of its publication, the problematic nature of its reported sequence and
raw sequencing reads, and the claim that no sample is left for independent verification. Collectively, these
facts justify and legitimate the concern over the true existence of the RaTG13 virus in nature and the
truthfulness of its reported genomic sequence. They also question the claim that the RaBtCoV/4991 virus
and RaTG13 are equivalent.
1.3 Genetic evidence proving the fraudulent nature of RaTG13
This evidence was revealed after a close examination of the sequences of specific genes, especially
spike, of relevant viruses. Specifically, we compared two viruses for the synonymous and non-
synonymous mutations on each gene, and we did so for two pairs of viruses. The first pair are bat
coronaviruses ZC45 and ZXC21. The second pair are SARS-CoV-2 and RaTG13. The rationale for
comparing these two pairs with each other is the following. First, ZC45 and ZXC21, each sharing an 89%
genomic sequence identity with SARS-CoV-2, are the closest relatives to SARS-CoV-2 and RaTG13.
Second, ZC45 and ZXC21 are 97% identical to each other, while SARS-CoV-2 and RaTG13 are 96%
identical. Not only the sequence identity in each case is comparable, but also the high sequence identity
indicates that, within each pair, the sequence difference should be a result of random mutations during
evolution, which ensures that synonymous and non-synonymous analyses here are appropriate and not
complicated by abrupt evolutionary events (e.g. recombination). Indeed, sequence alignment confirms
such a scenario – in both cases, the curve is smooth and the high sequence identity is maintained
throughout (Figure 3).
7
Figure 3. Simplot analyses show that high sequence identities are shared by two pairs of coronaviruses. A. the
genomic sequence of RaTG13 is plotted against that of SARS-CoV-2. B. Genomic sequence of ZXC21 is plotted
against that of ZC45.
Detailed synonymous (syn, green curve) and non-synonymous (non-syn, red curve) analyses are shown
in Figure 4. For each gene, the accumulations of syn and non-syn mutations, respectively, are illustrated
when the codons are analyzed in a sequential order. For the spike genes, between ZC45 and ZXC21, the
syn/non-syn ratio is 5.5:1 (Figure 4A left, 94 syn mutations and 17 non-syn mutations). Notably, the two
curves progress along in a roughly synchronized manner. These features reflect, to a certain extent, the
evolutionary traits resulted from random mutations during evolution in this sub-group of lineage B β
coronaviruses.
The same analysis on the spike genes of SARS-CoV-2 and RaTG13, however, revealed a different
scenario (Figure 4B right). Although the overall syn/non-syn ratio is a similar 5.4:1 (221 syn mutations
and 41 non-syn mutations), the synchronization between the two curves is non-existent. In the second half
of the sequence, which is over 700 codons (2,100 nucleotides) wide, the non-syn curve stays flat when the
syn curve climbs continuously and significantly.
Counting the syn and non-syn mutations of the S2 region (corresponding to residues 684-1273 of the
SARS-CoV-2 Spike) reveals that, between ZC45 and ZXC21, there are a total of 27 syn mutations and 5
non-syn mutations, yielding a syn/non-syn ratio of 5.4:1. In contrast, for the same S2 region, between
SARS-CoV-2 and RaTG13, there are a total of 88 syn mutations and 2 non-syn mutations, yielding a
syn/non-syn ratio of 44:1. The syn/non-syn ratios for S2, whole Spike, and other large viral proteins (Orf1a,
Orf1b, and Nucleocapsid) are summarized in Table 1. While the ratios are comparable between the two
groups for all other proteins, the ratios for the S2 protein are significantly different.
8
Figure 4. Abnormal distribution of synonymous and non-synonymous mutations in Spike revealed by the
comparison between RaTG13 and SARS-CoV-2. Synonymous and non-synonymous mutations are analyzed
between closely related coronaviruses on large viral proteins: A. Spike (S), B. Orf1a, C. Orf1b, and D. Nucleocapsid
9
(N). In each panel, the left graph is the comparison between the two bat coronaviruses ZC45 (MG772933) and
ZXC21 (MG772934), while the right graph is the comparison between SARS-CoV-2 (NC_045512) and RaTG13
(MN996532). In each graph, the accumulative growth of synonymous mutations (green curve), non-synonymous
mutations (red curve), and in-frame deletions (blue curve) are depicted, respectively. Initial sequence alignment
was done using EMBOSS Needle, which was followed by codon alignment at www.hiv.lanl.gov. Synonymous non-
synonymous analyses were performed using SNAP also at www.hiv.lanl.gov40.
The detailed syn/non-syn analyses for Orf1a, Orf1b, and N are shown in Figure 4B-D. It is also
noteworthy that, similar to that of Spike, the approximate synchronization between two curves is observed
for the Orf1a protein in the ZC45 and ZXC21 comparison (Figure 4B left) but not in the SARS-CoV-2
and RaTG13 comparison (Figure 4B right).
The S2 protein maintains trimmer formation of the Spike and, upon successive cleavages to expose the
fusion peptide, mediates membrane fusion and cell entry. Although the S2 protein is more conserved
evolutionarily than S1, the extremely high purifying pressure on S2 as suggested by the very high syn/non-
syn ratio is abnormal. In fact, Orf1b is known to be the most conserved protein in coronaviruses and yet
the syn/non-syn ratio for it is only 10.8:1 when SARS-CoV-2 and RaTG13 are compared, much lower
than the ratio of 44:1 observed for S2 (Table 1). Furthermore, since RaTG13 and SARS-CoV-2 infect
different species, no high purifying selection on S2 should be expected when these two viruses are
compared against each other.
Figure 5. Positive selection, not purifying selection, is observed for Spike in twenty randomly selected SARS-
CoV-2 sequences. GenBank accession numbers are shown in Figure 6. Collection dates of these viruses range from
December 2019 to July 2020.
10
Figure 6. Five amino acid mutations are observed in S2 (684-1273) in twenty randomly selected SARS-CoV-2
sequences. They are at positions 829, 1020, 1101, 1176, and 1191. GenBank accession number for each isolate is
shown in the sequence’s name following the country name.
11
Consistent with the above notion, a syn/non-syn analysis done for the Spike protein of twenty randomly
selected SARS-CoV-2 sequences showed that S2 was under positive selection, not purifying selection,
during the past eight months of human-to-human transmission (Figure 5). For the twenty SARS-CoV-2
isolates, amino acid mutations are observed at five different locations in S2 (Figure 6). In addition, a recent
study analyzing 2,954 genomes of SARS-CoV-2 revealed that mutations have been observed at 25
different locations in the S2 protein41, further proving that amino acid mutations are tolerated in S2 and
no high purifying pressure should be observed for S2. Evidently, the syn/non-syn ratio of 44:1 revealed
between SARS-CoV-2 and RaTG13 on the S2 region is abnormal (Table 1) and a violation of the
principles of natural evolution.
A logical interpretation of this observation is that SARS-CoV-2 and RaTG13 could not relate to
each other through natural evolution and at least one must be artificial. If one is a product of natural
evolution, then the other one must be not. It is also possible that neither of them exists naturally.
If RaTG13 is a real virus that truly exists in nature, then SARS-CoV-2 must be artificial.
However, the reality is that SARS-CoV-2 is physically present and has first appeared prior to the
reporting of RaTG134. This would then lead to the conclusion that RaTG13 is artificial, a scenario
consistent with the overwhelming suspicion that this virus does not exist in nature and its sequence has
been fabricated.
The remaining possibility is, of course, that both SARS-CoV-2 and RaTG13 are artificial: one has been
created physically and the other one exists only in the form of a fabricated sequence.
It is highly likely that the sequence of the RaTG13 genome was fabricated by lightly modifying the
SARS-CoV-2 sequence to achieve an overall 96.2% sequence identity. During this process, much editing
must have been done for the RBM region of the S1/spike because the encoded RBM determines the
interaction with ACE2 and therefore would be heavily scrutinized by others. An RBM too similar to that
of SARS-CoV-2 would be troublesome because: 1) RaTG13 could be conceived as a product of gain-of-
function research; 2) it would leave no room for an intermediate host and yet such a host is believed to
exist as the Spike/RBM needs to first adapt in an environment where the ACE2 receptor is homologous
to hACE2. In addition, modifying the sequence of the RBM is also beneficial as RaTG13 would otherwise
appear to be able to infect humans as efficiently as SARS-CoV-2 does, escalating the concern of a
laboratory leak. To eliminate such concerns, many non-syn mutations were introduced into the RBM
region.
Importantly, syn/non-syn analysis is frequently used, often at the ORF/protein level, to characterize the
evolutionary history of a virus42-44. While editing the RBM, the expert(s) carrying out this operation must
be conscious of the need to maintain a reasonable syn/non-syn ratio for the whole Spike protein. To
achieve so, however, the expert(s) must have then strictly limited the number of non-syn mutations in the
S2 half of Spike, which ended up flattening the curve (Figure 4A right).
1.4 The receptor-binding domain (RBD) of RaTG13 does not bind ACE2 of horseshoe bats
Consistent with the above conclusion that RaTG13 does not exist in nature and its sequence was
fabricated, a recent study showed that the RBD of RaTG13 could not bind the ACE2 receptors of two
different kinds of horseshoe bats, Rhinolophus macrotis and Rhinolophus pusillus45. Although the ACE2
receptor of Rhinolophus affinis (the alleged host of RaTG13) was not tested, it is unlikely that ACE2 of
R. affinis would differ significantly from those of its close relatives and be able to bind the RaTG13 RBD.
12
This result therefore implicates that RaTG13 would not be able to infect horseshoe bats, contradicting the
claim made by Shi and colleagues that the virus was detected and discovered from horseshoe bats. This is
also consistent with the above conclusion that the genomic sequence of RaTG13 is fabricated and
presumably computer-edited, which entails that the RBM/RBD suggested by the corresponding gene
sequence may not be functional in binding the ACE2 receptor of the claimed host.
1.5 Conclusion and postulation of the fabrication process
In conclusion, the evidence presented both here and from recent literature collectively prove that
RaTG13 does not exist in nature and its sequence has been fabricated.
If the RaBtCov/4991 virus is equivalent to RaTG13, then RaBtCoV/4991 must be fraudulent as well.
Apparently, in the actual process of sequence fabrication, the published sequence of the short RdRp
segment of RaBtCoV/4991 was completely inherited for RaTG13. This way, they could claim that
RaTG13 was RaBtCoV/4991, which, according to the record, was discovered in 201326. If RaTG13 had
been described as being discovered right around the time of the COVID-19 outbreak, greater suspicions
would result as tracing the evolutionary origin of a zoonotic virus is difficult and usually takes years or
decades. As described in section 2.1 of our earlier report1, the fabrication of RaTG13 should have been
planned and executed in coordination with the laboratory creation of SARS-CoV-2.
Such an approach is also safe because, except for the 440-bp RdRp segment, no other sequence
information has ever been published for the rest of the RaBtCoV/4991 genome.
It is worth noting that, due to reasons detailed in section 1.2, they still preferred to obscure the history
of RaTG13. However, they must have also anticipated that their violations of the publication norms would
invite inquiries or requests for clarifications, the number of which, however, should be limited and
manageable. RaBtCoV/4991 would then function as an additional layer of security for them in facing such
inquiries and/or requests.
Building upon the 440-bp RdRp sequence inherited from RaBtCoV/4991, the rest of the RaTG13
genome was likely fabricated by lightly editing the sequence of SARS-CoV-2. Once the genomic sequence
was finalized, DNA fragments could be synthesized individually according to the fabricated and edited
sequence and then used as templates for PCR. Amplified DNA would then be mixed with certain raw
material to give the sample a natural look (mimicking what is present in an actual RT-PCR, which is done
using RNA extracted from fecal swabs as templates). Subsequently, this sample would be sent for
sequencing. The resulted raw sequencing reads could then be uploaded together with the made-up genomic
sequence onto GenBank to create an entry for the RaTG13 genome.
1.6 The Mojiang Miner Passage (MMP) hypothesis is fatally flawed
Recently, a theory has emerged, which proposed that SARS-CoV-2 was derived from viral passaging
in the lungs of the infected Mojiang miners back in 201246. Specifically, authors believe that the
RaBtCoV/4991 virus was indeed RaTG13 and was the virus causing pneumonia in the miners in 2012.
While inside the lungs of the miners, the RaTG13 virus had evolved extensively, mimicking a viral
passage process, and eventually became SARS-CoV-2. In this process, the RBD of the virus experienced
strong positive selection, through which it became optimal in binding hACE2. Furthermore, the furin-
cleavage site at the S1/2 junction region of Spike had been acquired through recombination between the
viral spike gene and the gene encoding the human ENaC protein, which has a furin-cleavage sequence
closely resembling that of SARS-CoV-2. The end product of this passage was SARS-CoV-2, which the
13
researchers isolated from the miners’ samples and brought back to the WIV. The authors have named this
hypothesis as the Mojiang Miner Passage (MMP) hypothesis46.
However, this MMP hypothesis has fatal flaws.
First, the viral pathogen that caused the disease in the miners could not be defined or confirmed.
According to the record, which was well documented in a Master’s Thesis written by the doctor in charge,
samples from two patients (throat swabs and blood) were tested at the Center for Disease Control and
Prevention of the Chengdu Military Region between May 15th and May 20th, 2012, and yet none of the
suspected viruses, including SARS, was detected30. Furthermore, the gold standard in the clinical
diagnosis of coronavirus-caused pneumonia is biopsy and/or autopsy followed by confirmation by either
RT-PCR or isolation of the virus. However, three biopsy tests were attempted but failed30. Autopsy tests
were requested and yet all turned down by families of the deceased miners30. Due to such failure, both the
Master’s Thesis and later a PhD Dissertation, which also looked into this issue although in an indirect
manner, described the cause of the pneumonia as an unsolved case30,31.
Second, antibody tests done for the miners do not support SARS or SARS-like coronavirus infection.
According to the Master’s Thesis, samples from two miners were tested for antibodies against SARS30.
The symptoms onset date for one miner (case 3, passed away) was around April 13th, 2012. The other
miner (case 4, had severe symptoms and yet recovered) had symptoms onset around April 16th, 2012.
Antibody tests, which were recommended later by Dr. Nanshan Zhong, were done at the WIV on June
19th, 2012. However, the two samples tested were only positive for IgM30. No positive IgG or total
antibody were reported30. No antibody titer was reported either. Importantly, if the severe pneumonia was
caused by coronavirus infections, by the time of the antibody tests on June 19th, 2012, both IgM and
IgG/total antibody should be detected. In fact, IgG/total antibody should be much more abundant and
easier to detect47. On the other hand, IgM tests frequently result in false positives48. Therefore, the fact
that only IgM, and no IgG/total antibody, was tested positive suggests that the described results were most
likely false positives and the infections should not have been caused by SARS or a SARS-like coronavirus.
It is noteworthy that the later PhD Dissertation31 showed severe discrepancies with the Master’s Thesis
in the descriptions of the same clinical tests:
1. The PhD Dissertation described that samples from four miners (throat swab and blood) were sent
to the Center for Disease Control and Prevention of the Chengdu Military Region for nucleic acid
tests. However, the Master’s Thesis indicated that samples were only taken from two miners30.
2. The PhD Dissertation described samples from four miners were tested for anti-SARS antibodies
at the WIV and all were IgG positive. However, the Master’s Thesis indicated that only samples
from two miners were tested at the WIV and both were only IgM positive30.
Importantly, the Master’s Thesis was written in 2013 in Yunnan by the doctor who was in charge of
the six hospitalized miners30. The PhD dissertation, however, was written in 2016 in Beijing based only
on the clinical record. The author of the Dissertation had no direct involvement in the treatment of the
miners or in any of the described tests31. It is therefore highly likely that author of the PhD dissertation
did not verify the clinical data he presented, which makes this PhD dissertation an unreliable source of
information concerning the Mojiang miners’ case.
Third, if SARS-CoV-2 was already present in the miner’s body in 2012, it would have certainly caused
an epidemic or even pandemic then. Given the extremely high transmissibility of SARS-CoV-2, it would
be impossible for the doctors, nurses, family members of the miners, etc. to have avoided contracting the
14
virus without the protection of proper PPE. If an epidemic indeed happened in 2012, it could not have
gone unnoticed given the high transmissibility and lethality (three out of the six pneumonia patients died
despite of intense medical care provided for them).
Fourth, as shown in sections 1.1-1.5, RaTG13’s sequence is clearly fabricated and the virus does not
exist in nature. The RaBtCoV/4991 virus, which was detected in 2013, is not the RaTG13 virus that is
defined by its reported genomic sequence. No complete genomic sequence of RaBtCoV/4991 has ever
been reported likely due to the poor quality of the sample, which happens often as the RNA genome decays
easily. It is highly likely that no high homology is shared between the actual RaBtCoV/4991 virus and
SARS-CoV-2. This judgement is based on the fact that no viruses reported prior to 2020 share more than
90% sequence identity with SARS-CoV-2 despite the extensive surveillance studies of coronaviruses for
the past two decades. Therefore, even if RaBtCoV/4991 was the pathogen responsible for the pneumonia
of the miners, the theory that it has evolved in a single person’s lung into SARS-CoV-2 is far beyond
being reasonable.
Fifth, it is impossible for the Spike protein of the virus to obtain a unique furin-cleavage site at the
S1/S2 junction through recombination with the gene encoding the ENaC protein of the host cell (ENaC
carries a furin cleavage site closely resembling the one seen in SARS-CoV-2). This is because
recombination requires a significant level of sequence similarity between the two participating genes and
yet no such similarity is present between coronavirus Spike and human ENaC. The molecular basis for
recombination is non-existent. (Although recombination between ENaC and coronavirus Spike is
impossible, it is suspicious that a viral protein and a host protein would share the same sequence for their
furin-cleavage sites. It is possible, though, that the sequence of the furin-cleavage site in ENaC49, which
is known since 199750, could have been used in the design of the furin-cleavage site in the Spike of SARS-
CoV-2. Such a design may be considered sophisticated as ENaC co-expresses with ACE2 in many different
types of cells49.)
Sixth, if SARS-CoV-2 has indeed evolved from RaBtCoV/4991 in the miner’s lungs, it would look,
from every aspect, like a naturally occurring virus. In that case, there would be no need to commit sequence
fabrication for RaTG13 and for the other novel coronaviruses (parts 2 and 3) to falsify a natural origin for
SARS-CoV-2.
Finally, as revealed in our earlier report1, evidence exists in the genome of SARS-CoV-2, indicating
that genetic manipulation is part of the history of SARS-CoV-2.
2. Evidence proving that recently published pangolin coronaviruses are fraudulent and do not exist
in nature
While RaTG13 was reported to share a high sequence identity with SARS-CoV-2 and thereby hinted
a natural origin of SARS-CoV-2, significant questions remained unanswered:
• No intermediate host has been found although one was believed to exist and function as the
reservoir of the virus before it spilled over to humans.
• Despite the overall genomic resemblance of the two viruses, the RBD (particularly the RBM within
it) of RaTG13 differs significantly from that of SARS-CoV-2. The evolutionary origin of the
SARS-CoV-2 RBD, which is optimal in binding hACE2, remained unclear.
15
• A critical furin-cleavage site, which is present at the S1/S2 junction of SARS-CoV-2 Spike and
responsible for the enhanced viral infectivity and pathogenicity51-57, is absent in RaTG13 (as well
as in all known lineage B β coronaviruses58). The evolutionary origin of this furin-cleavage site
also remained mysterious.
Not long after these questions emerged, several laboratories published novel coronaviruses allegedly
found in Malayan pangolins that were smuggled from Malaysia and confiscated by the Chinese custom5-
8
. Although these novel coronaviruses share relatively lower overall sequence identities (~90%) with
SARS-CoV-2 in comparison to RaTG13 (96.2% identical to SARS-CoV-2), the RBD of the pangolin
coronaviruses resembles greatly the SARS-CoV-2 RBD (97.4% identical). In the most critical RBM
region, all amino acids except one are identical between the pangolin coronaviruses and SARS-CoV-25-8.
These observations led the authors to conclude 1) that pangolins are the likely intermediate host for the
zoonotic transfer of SARS-CoV-25,7 and 2) that a RaTG13-like ancestor coronavirus might have acquired
the RBD from a pangolin coronavirus through recombination to eventually become SARS-CoV-25-8.
Here, in part 2 of the report, we describe literature evidence and provide genetic analyses to prove that
these novel pangolin coronaviruses are results of fabrication.
2.1 A single batch of pangolin samples were used in all studies and the deposited sequencing data
showed heavy contamination and signs of fabrication
In October 2019, a team formed by three researchers from two institutions (Guangdong Institute of
Applied Biological Resources and Guangzhou Zoo) reported, for the first time, the detection of
coronavirus infections in pangolins that were allegedly smuggled from Malaysia and confiscated in the
Guangdong province in March 201959. Twenty-one pangolin samples were sequenced and five were
positive for coronavirus infections (Table 2: lung 2, 7, 8, 9, and 11), although Sendai virus infection was
also reported. However, neither the sequences of the coronaviruses nor raw sequencing data were made
available to the public for a period of three months. The raw data (NCBI BioProject PRJNA573298) was
finally released on January 22nd, 2020 after the COVID-19 outbreak started, while the article submission
date was September 30th, 2019 and the publication date was October 24th, 201959.
Between March and May 2020, four seemingly independent studies were published, all of which
reported novel pangolin coronaviruses and their assembled genomic sequences5-8. However, after a closer
look, we found that all four studies derived viral sequences from the same set of pangolin samples first
reported in the October 2019 publication59, which has been confirmed by a recent article13.
In one study6, Liu et al. (the same authors of the October 2019 publication59) re-assembled the genome
of a pangolin coronavirus by pooling two samples from the original 2019 study and one sample obtained
from another Malayan pangolin rescued in July 2019. However, although the authors stated that the more
recent raw sequencing data had been deposited at the NCBI database6, we could not find this data using
the accession number (2312773) provided. The same difficulty has been reported by others13. Therefore,
it cannot be verified whether the July 2019 dataset truly exists and has contributed to the assembly of the
reported genome.
In two other studies, Lam et al.5 and Zhang et al.8 each re-assembled the genome of a pangolin
coronavirus using only the published dataset from the October 2019 study59. Lam et al. also reported
detection of coronaviruses from smuggled Malayan pangolins that were confiscated in the Guangxi
province5, although these viruses showed lower sequence identities to SARS-CoV-2 both at the whole
genome level (~86%) and in the critical RBD region. It is noteworthy that this study was done as a
16
collaboration between Dr. Yi Guan’s group from the University of Hong Kong and Dr. Wuchun Cao’s
group from the Academy of Military Medical Sciences (AMMS), Beijing, China5. Somehow, all authors
affiliated with the AMMS were excluded from the list of authors when the article was first submitted60,
although their names eventually appeared in the final version of the publication5.
In the fourth study, Xiao et al. claimed to have examined tissue samples kept from diseased pangolins
and obtained raw sequencing data for the subsequent assembly7. However, they did not describe how the
samples were acquired. In their Extended Data Table 3, they listed the metagenome sequencing data used
in the study7, which, surprisingly, do not match with the actual data that they uploaded in the database
(Table 2). Samples M1, M5, M6, M10, and Z1 can be found in the data they deposited, but not M2, M3,
M4, and M8. Furthermore, Xiao et al. apparently were inconsistent with the reporting of these raw
sequencing reads. For samples M1, M6, pangolin3, and pangolin5, they counted paired ends numbers,
which reflect the actual number of sequenced DNA fragments in the library. For the rest of samples, the
authors counted reads numbers instead (In Illumina sequencing, there are two reads per fragment). For
samples M2, M3, M4, and M8 in this latter group7, when the reads numbers were converted to paired
ends numbers (divided by 2), they each match perfectly with lung07, lung02, lung08, and lung11,
respectively, from the October 2019 study59 (Table 2). Clearly, Xiao et al. used the data published in a
previous study but failed to disclose this necessary information in their publication7. In fact, they
intentionally presented the “number of reads” in a different format to presumably make readers overlook
the fact that the same sequencing dataset was used.
It is noteworthy that the study by Xiao et al. was also done in collaboration with the AMMS. Prior to
the publication of the manuscript, this work was first publicized in a press conference61,62. As revealed in
this conference, four principle investigators contributed to the work and one of them was Dr. Ruifu Yang
from the AMMS. However, like what happened to Dr. Cao and his AMMS colleagues in the Lam et al.
study5, Dr. Yang’s name was excluded in the submitted manuscript of Xiao et al.63. Yet, unlike the other
case, the AMMS researcher’s name did not re-appear in the final publication7. It is also noteworthy that
the two AMMS principle investigators here, Dr. Yang and Dr. Cao, are long-term collaborators and most
of their collaborative work concerned genetic analyses of SARS-CoV64-67.
Among the four studies, only two assembled complete genomes by performing gap filling using PCR6,7.
However, neither group made their gap filling sequences available13, rendering independent verification
impossible. Notably, the delayed publishing of raw sequencing reads long after the publication of genomic
sequences has occurred in the reporting of RaTG13 as well.
Adding to the above problems was the poor quality of the raw sequencing data, which has been
described recently13,14,20. We also analyzed the composition of the sequencing reads of the deposited
libraries. By performing taxonomy analysis on the NCBI SRA database, we also found that samples from
Liu et al.6 that are positive for coronavirus reads are all positive for reads that map to human genome
(Table 2). In great contrast, the rest of the samples, which are negative for viral reads, also have no human
reads detected. The same correlation is found in data presented by Xiao et al7. Although samples M5
(pangolin 6) and M6 (pangolin2) are negative for human reads, these two samples have very few viral
reads, which would hardly contribute to the viral genome assembly. Clearly, the human contamination
should not be due to sample handling as none of the coronavirus-negative samples, which must have been
handled similarly, contain such contamination. The consistent co-existence of viral reads and human reads
are highly suspicious.
17
Table 2. Analyses of the raw sequencing data deposited by Liu et al.6 and Xiao et al.7
18
These observations raise red flags not only on the credibility of the assembled sequences but also on
the authenticity of these novel pangolin coronaviruses. It is also noteworthy that the manuscript
submission dates for all four studies were between February 7th and February 18th5-8, suggesting that their
publications might have been coordinated.
2.2 No coronavirus was detected in an extensive surveillance study of Malayan pangolins
While these SARS-CoV-2-like pangolin coronaviruses were described as being detected in smuggled
Malayan pangolins59, a recent study strongly refuted the presence of such pangolin coronaviruses in nature.
A team led by Dr. Daszak examined 334 pangolin samples, which were collected in Malaysia and Sabah
from August 2009 to March 201968. Surprisingly, no coronaviridae, or any of the other families of viruses
(filoviridae, flaviviridae, orthomyxoviridae, and paramyxoviridae), were detected in any of these samples.
This is in stark contrast with the October 2019 publication where both coronavirus infection and Sendai
virus infection were reportedly detected in the smuggled Malayan pangolins59, which eventually led to the
discovery and publication of the novel pangolin coronaviruses5-8. The finding of Lee et al.68 adds
significantly to the existing suspicions and substantiates the possibility that these pangolin coronaviruses
do not exist in nature and their sequences could have been fabricated.
2.3 The RBD of the reported pangolin coronaviruses binds poorly to pangolin ACE2
If pangolin coronaviruses truly exist and have recently spilled over to infect humans, their Spike protein,
especially the RBD within Spike, should bind to pangolin ACE2 (pACE2) more efficiently than to hACE2.
However, recent findings have contradicted this theory. In an in silico study, Piplani et al. calculated,
following homology structural modeling, the binding energies involved in the association between SARS-
CoV-2 Spike and ACE2 from either human or various animals69. Interestingly, the most favorable
interaction that SARS-CoV-2 Spike makes was shown to be with hACE2, but not with ACE2 from
pangolin or any other suspected intermediate host. Furthermore, another study revealed, using a robust in
vitro binding assay, that the RBD of SARS-CoV-2 binds much tighter (greater than 9-fold) to hACE2 than
to pACE245. Although the RBD of the pangolin coronaviruses is not 100% identical to that of SARS-
CoV-2, the RBMs of the two viruses, which is the most essential segment responsible for ACE2
interactions, differ only by one amino acid5-8. Therefore, the poor binding efficiency observed between
the RBD of SARS-CoV-2 and pACE245 infers that the RBD of the reported pangolin coronaviruses must
bind to pACE2 fairly inefficiently. Indeed, a very recent study confirmed the case: the RBD of the
pangolin coronavirus binds pACE2 ten-fold weaker than to hACE270. These observations once again
refute the claim that pangolins are the probable intermediate host for SARS-CoV-2. More importantly,
the latter two studies strongly suggest that these viruses might not be able to establish infections in
pangolins, which adds significantly to the suspicion that the published sequences of the pangolin
coronaviruses may have been fabricated and these viruses do not exist in nature.
2.4 Genetic evidence proving the fraudulent nature of the pangolin coronaviruses
Evolutionarily, within the coronavirus genome, the RBD of Spike is under the strongest positive
selection as it needs to adapt for binding a new receptor whenever the virus crosses the species barrier and
enters a new host. In lineage B β coronaviruses, the most essential segment for receptor recognition is the
RBM, which fully determines the binding with ACE2. Strikingly, when the RBM sequence of the pangolin
virus MP7896 is compared to that of SARS-CoV-2, no positive selection is observed (Figure 7A). Instead,
the analysis revealed very strong purifying selection with 24 syn mutations and only one non-syn mutation.
In contrast, when two related bat coronaviruses, BM48-3171 and BtKY7272, are compared in a similar
19
manner, strong positive selection is observed as expected (Figure 7B). Here, while there are 25 syn
mutations, which is comparable to that between MP789 and SARS-CoV-2, the number of non-syn
mutations is 30 (Figure 7B). Evidently, the species difference between pangolin and human is greater than
that between the hosts of BM48-31 and BtKY72, which are two different species of bats. Therefore,
greater positive selection should be expected between MP789 and SARS-CoV-2 than that between BM48-
31 and BtKY72. The strong purifying selection observed between MP789 and SARS-CoV-2 is, therefore,
contradictory to the principles of natural evolution.
Figure 7. The extremely high purifying pressure observed for the RBM in the comparison of pangolin
coronavirus MP789 and SARS-CoV-2 contradicts the principles of natural evolution. Synonymous and non-
synonymous mutations in the RBM region are analyzed between related coronaviruses: A. pangolin coronavirus
MP789 (MT121216.1) and SARS-CoV-2 (NC_045512.2), B. bat coronaviruses BM48-31 (NC_014470.1) and
BtKY72 (KY352407.1), and C. bat coronaviruses ZC45 and ZXC21. D. Alignment of the RBM sequences from all
six viruses. The beginning and end of the RBM are labeled following the sequence of the SARS-CoV-2 Spike.
Table 3. Summary of syn/non-syn mutations in the RBM in three sets of pair-wise comparisons
Figure 8. Abnormal distribution of synonymous and non-synonymous mutations in Spike associated with
pangolin coronaviruses. A. Comparison between MP789 and P4L (MT040333.1). B. Comparison between the two
bat coronaviruses BM48-31 and BtKY72.
Table 4. Ratios of syn/non-syn mutations observed in different viral proteins as revealed by pair-wise
comparisons involving pangolin and bat coronaviruses
21
Similar syn/non-syn analyses on the overall spike further revealed the fraudulent nature of these novel
pangolin coronaviruses. Here we compared two representative pangolin coronaviruses MP7896 (a
Guangdong isolate) and P4L5 (a Guangxi isolate) as genomic sequences within each group of isolates
share very high sequence identities13. As shown in Figure 8A, similar to the abnormal pattern observed
between RaTG13 and SARS-CoV-2 (Figure 4A right), syn and non-syn curves exhibit drastically different
trajectories and the non-syn curve abruptly flattens in the S2 half of the sequence.
For comparison, we also analyzed the spike genes of two SARS-like bat coronaviruses, BM48-31 and
BtKY72. The two pangolin coronaviruses, MP789 and P4L, are 85.2% identical on the overall genome,
while bat coronaviruses BM48-31 and BtKY72 are 82.4% identical. The comparison here is therefore
appropriate. Analysis of the two bat viruses show that the two curves grow naturally in a relatively
concerted manner with no excessive flattening of the red curve observed (Figure 8B).
Counting the number of syn and non-syn mutations in each pair of comparisons further illustrated the
unnatural characteristics associated with the pangolin coronaviruses (Table 4). While the S2 protein is not
expected to be more conserved than Orf1b, the syn/non-syn ratio for S2 observed in the comparison
between MP789 and P4L is abnormally high (207 syn mutations and 9 non-syn mutations; syn/non-syn =
23:1), which is far exceeding what is observed for Orf1b (7.6:1).
As the two bat coronaviruses here were discovered in nature independently by research groups outside
of China71,72, the features displayed in Figure 8B likely represent the approximate evolutionary trait of
two coronaviruses at this level of overall divergence. According to the logic described earlier, the great
contrast between Figure 8A and 8B and the abnormal syn/non-syn ratio of 23:1 (Table 4) further prove
that, between MP789 and P4L, at least one is artificial, although we believe both groups of pangolin
coronaviruses represented by MP789 and P4L, respectively, are non-natural and fabricated.
2.5. Summary and discussion
A single source of samples was used for all studies (some spuriously independent7) reporting novel
pangolin coronaviruses. The formats of sequencing reads were manipulated with a clear intention to hide
the fact that the same dataset was used in different studies. The raw sequencing data is missing for certain
critical pieces, poor in quality, and suspicious in terms of the amounts and types of contaminations present.
The RBD encoded by the reported sequence of pangolin coronaviruses could not bind pACE2 efficiently.
As revealed by syn/non-syn analyses, sequences of the RBM and S2 regions of these pangolin
coronaviruses exhibit features that are inconsistent with natural evolution. Finally, no coronavirus was
detected in a large, decade-long surveillance study of Malayan pangolins. These observations and
evidence converge to prove that these recently reported pangolin coronaviruses do not exist in nature and
their sequences must have been fabricated.
It is noteworthy that the abnormal syn/non-syn feature revealed for S2 in the comparison between
MB789 and P4L (Figure 8A) resembles greatly that exhibited by the comparison between RaTG13 and
SARS-CoV-2 (Figure 4A right). Judging based on this reoccurring pattern, we believe that the sequence
fabrications in both cases (RaTG13 and pangolin coronaviruses) were most likely carried out by the same
person or group, whose misconception of the spike gene evolution persisted in multiple such practices and
resulted in the unnatural look of the syn/non-syn curves and numbers (Figure 4, Table 1, Figure 8, and
Table 4).
22
3. Evidence revealing the fraudulent nature of the novel bat coronavirus RmYN02
While the publications of the fabricated pangolin coronaviruses might have seemingly fulfilled the
scientific quests for an intermediate host for the zoonosis of SARS-CoV-2 as well as for an evolutionary
origin of its RBD, it had remained suspicious and unexplainable how SARS-CoV-2 could have acquired
the furin-cleavage site (-PRRAR/VS-) at the S1/2 junction through natural evolution. It is evident that,
although furin-cleavage site has been found in certain other lineages of coronaviruses at the S1/2 junction,
lineage B β coronaviruses clearly lack the ability to develop this motif at this location naturally58.
In early June, another novel bat coronavirus, RmYN02, was reported9, which shares a 93.3% sequence
identity with SARS-CoV-2 and appears to be the second closest bat coronavirus to SARS-CoV-2 (the
closest is allegedly RaTG13). This finding adds yet another member to the rapidly growing sub-lineage
of SARS-CoV-2-like coronaviruses (Figure 9), which has been completely vacant and practically non-
existent prior to the current pandemic. In addition, importantly, RmYN02 carries a unique sequence -PAA-
at the S1/S2 junction, which remotely resembles the inserted -PRRA- sequence at the same location in the
SARS-CoV-2 Spike. Despite the fact that -PAA- in RmYN02 only partially resembles the -PRRA-
insertion in SARS-CoV-2 and does not appear to be an actual insertion if properly aligned18, the authors
nonetheless claimed that the natural occurrence of -PAA- in RmYN02 proves that the -PRRA- sequence
could very likely be acquired and “inserted” into the same location in SARS-CoV-2 genome through
natural evolution9.
Figure 9. Phylogenetic analysis of SARS-CoV-2 and representative viruses from the subgenus sarbecoronavirus.
Figure redrew from Zhou et al9. Colored viruses were all reported after the COVID-19 outbreak.
23
The fact that a poor alignment was used to make a disproportional, strong argument for an evolutionary
origin of the furin-cleavage site, which appeared to be the last missing piece of the puzzle, is suspicious.
Furthermore, despite the significance of the spike sequence of RmYN02 in supporting the central
conclusion of the publication, the raw sequencing reads for spike has not been made available although
the authors stated otherwise in the article9. This is yet another repeat of the pattern that has been exhibited
in the reporting of both RaTG13 and pangolin coronaviruses, where the genomic sequence would be
published first and the raw sequencing reads would not be made available months afterwards.
Given that the CCP-controlled laboratories have repeatedly engaged in fabrication of coronaviruses to
feed the missing pieces for the puzzle, the above suspicion opens up the possibility that the RmYN02 virus
could have been fabricated as well. Judging from the fact that its sequence identity to SARS-CoV-2
(93.3%) is lower than that between RaTG13 and SARS-CoV-2 (96.2%), we suspected that the sequence
of RmYN02 might have been fabricated by modifying the sequence of RaTG13. Such an approach could
easily ensure that the evolutionary distance between RmYN02 and SARS-CoV-2 is greater than that
between RaTG13 and SARS-CoV-2. It also ensures that RmYN02 and RaTG13 would appear to be
evolutionarily close, consistent with the claim that they both infect bats although of different species.
We therefore compared the spike genes of RmYN02 and RaTG13 on the quantity and distribution of
syn and non-syn mutations. The severe divergence at the S1 portion between the two viral sequences did
not allow the S1 sequences to be properly codon-aligned. Therefore, only the S2 half was analyzed (Figure
10). For the beginning 200 codons of S2, both types of mutations accumulate steadily and gradually.
However, for the final 378 codons, once again, the non-syn curve flattens and the concerted growth of the
two curves has disappeared. In this region, there are 57 syn mutations and only one non-syn mutation. The
syn/non-syn ratio of 57:1 for a region as wide as 378 codons (1,344 nucleotides) is severely inconsistent
with what is observed naturally (Figure 4A left and Figure 8B)41.
Figure 10. Analysis of synonymous and non-synonymous mutations in S2 between RmYN02 and RaTG13. The
abrupt change of trajectory of the non-synonymous mutation (red) curve and its subsequent flattening are observed.
Logically, between RaTG13 and RmYN02, at least one must be artificial. Here, however, we are
convinced that both viruses are artificial. As shown in part 1, the sequence of RaTG13 must have been
fabricated. Therefore, the fact that the last 378 codons of RmYN02’s S2 are identical, with the exception
of one, to that of RaTG13 proves that the RmYN02 sequence must be artificial as well. This also proves
24
our earlier suspicion that the RaTG13 sequence should have been used as the template for the fabrication
of the RmYN02 sequence. RaTG13 was published in late January4, while RmYN02 was published in early
June (manuscript submitted in April)9. Therefore, enough time is in between for the sequence fabrication
to be carried out.
While introducing nucleotide changes to create the apparent divergence between the two viruses, the
expert(s) may have overly restricted amino acid changes in this part of Spike. Again, the abrupt change of
trajectory of the non-syn curve and its excessive flattening later in the sequence likely reflect their
overestimation of the purifying selection pressure on S2. The fact that this abnormal pattern has been
observed in all three cases (Figure 4A right, 8A, and 10) reiterates the point raised in section 2.5 that all
sequence fabrications may have been carried out by the same person or group.
Figure 11. All novel coronaviruses recently published by the CCP-controlled laboratories share a 100% amino
acid sequence identity on the E protein with ZC45 and ZXC21. Additional accession numbers of viruses: SARS-
CoV-2 (NC_045512.2), Pangolin-CoV (EPI_ISL_410721), P5E(MT040336.1), P3B(MT072865.1),
P2V(MT072864.1), P5L(MT040335.1), and P1E(MT040334.1).
A main goal of these fabrications was to obscure the connection between SARS-CoV-2 and
ZC45/ZXC21. Therefore, from their perspective, the fabricated viruses should resemble SARS-CoV-2
25
more than ZC45 and ZXC21 do. Because ZC45 and ZXC21 already share a 100% identity with SARS-
CoV-2 on the E protein, the fabricated viruses therefore were made to adopt this sequence completely as
well.
4.2 Important implications of this large-scale, organized scientific fraud
If SARS-CoV-2 is of a natural origin, no fabrications would be needed to suggest so. The current report,
therefore, corroborates our earlier report and further proves that SARS-CoV-2 is a laboratory product1.
As revealed1, the creation of SARS-CoV-2 is convenient by following established concepts and techniques,
some of which (for example, restriction enzyme digestion) are considered classic and yet still preferred
widely including by experts of the field35,73. A key component of the creation, the template virus
ZC45/ZXC21, is owned by military research laboratories3.
Importantly, as revealed here, multiple research laboratories and institutions have engaged in the
fabrication and cover-up4-9,59. It is clear that this was an operation orchestrated by the CCP government.
In addition, raw sequencing reads for RaTG13, which were integral parts of the fabrication, were
obtained in 2017 and 201824,33. Furthermore, manuscript reporting the falsified coronavirus infections of
Malayan pangolins was submitted for publication in September 201959. Evidently, the cover-up had been
planned and initiated before the COVID-19 outbreak. Therefore, the unleashing of the virus must be a
planned execution rather than an accident.
4.3 SARS-CoV-2 is an Unrestricted Bioweapon
Although it is not easy for the public to accept SARS-CoV-2 as a bioweapon due to its relatively low
lethality, this virus indeed meets the criteria of a bioweapon as described by Dr. Ruifu Yang. Aside from
his appointment in the AMMS, Dr. Yang is also a key member of China’s National and Military
Bioterrorism Response Consultant Group and had participated in the investigation of the Iraqi bioweapon
program as a member of the United Nations Special Commission (UNSCOM) in 1998. In 2005, Dr. Yang
specified the criteria for a pathogen to qualify as a bioweapon74:
1. It is significantly virulent and can cause large scale casualty.
2. It is highly contagious and transmits easily, often through respiratory routes in the form of aerosols.
The most dangerous scenario would be that it allows human-to-human transmission.
3. It is relatively resistant to environmental changes, can sustain transportation, and is capable of
supporting targeted release.
All of the above have been met by SARS-CoV-2: it has taken hundreds of thousands lives, led to
numerous hospitalizations, and left many with sequela and various complications; it spreads easily by
contact, droplets, and aerosols via respiratory routes and is capable of transmitting from human to
human75-77, the latter of which was initially covered up by the CCP government and the WHO and was
first revealed by Dr. Li-Meng Yan on January 19th, 2020 on Lude Press78; it is temperature-insensitive
(unlike seasonal flu) and remains viable for a long period of time on many surfaces and at 4°C (e.g. the
ice/water mixture)79,80.
Adding to the above properties is its high rate of asymptomatic transmission, which renders the control
of SARS-CoV-2 extremely challenging. In addition, the transmissibility, morbidity, and mortality of
SARS-CoV-2 also resulted in panic in the global community, disruption of social orders, and decimation
of the world’s economy. The range and destructive power of SARS-CoV-2 are both unprecedented.
26
Clearly, SARS-CoV-2 not only meets but also surpasses the standards of a traditional bioweapon.
Therefore, it should be defined as an Unrestricted Bioweapon.
4.4 The current pandemic is an attack on humanity
The scientific evidence and records indicate that the current pandemic is not a result of accidental
release of a gain-of-function product but a planned attack using an Unrestricted Bioweapon. The current
pandemic therefore should be correspondingly considered as a result of Unrestricted Biowarfare.
Under such circumstances, the infected population are being used, unconsciously, as the vectors of the
disease to facilitate the spread of the infection. The first victims of the attack were the Chinese people,
especially those in the city of Wuhan. At the initial stage, the hidden spread in Wuhan could have also
served another purpose: the final verification of the bioweapon’s functionality, an important aspect of
which is the human-to-human transmission efficiency. Upon the success of this last step, targeted release
of the pathogen might have been enabled.
Given the global presence of SARS-CoV-2 and the likelihood of its long-term persistence, it is
appropriate to say that this attack was on the humanity as a whole and has put its fate at risk.
4.5 Actions need to be taken to combat the current pandemic and save the future of humanity
Given the CCP’s role here, it is of paramount importance that the CCP is held accountable for its actions.
In addition, the world needs to find out what other variants of SARS-CoV-2 exist in the CCP-controlled
laboratories, whether or not SARS-CoV-2 or its variant(s) are still being actively released, whether or not
re-infection of SARS-CoV-2 leads to worsened outcomes due to inefficient immunity and/or antibody-
dependent enhancement (ADE)81-83, and whether other weaponized pathogens are owned by the CCP as a
result of their excessive, state-stimulated efforts in collecting novel animal pathogens and studying their
potentials in zoonosis3,25,26,28,32,36,37,84-114.
It is also of paramount importance that all the hidden knowledge of SARS-CoV-2 be brought out as
soon as possible. As illustrated in our earlier report, although a template virus was used, the creation of
SARS-CoV-2 must have involved introducing changes to the template sequence through DNA synthesis
(steps 1 and 4 in part 2 of our earlier report)1. Such a practice can be safely guided by multi-sequence
alignment of available SARS and SARS-like coronavirus sequences. The process of this practice has been
illustrated115, and both syn mutations and amino acid (non-syn) mutations at variable positions/regions
would be introduced. From the perspective of the responsible scientists, these changes are necessary
because, otherwise, the engineered nature of the virus and its connection to its template would be evident.
However, importantly, the introduced changes might have also altered the functions of the various viral
components, which could be either by design or unintended. Nonetheless, it remains to be answered
whether or how the introduced changes might be responsible for the various lasting complications that
many COVID-19 patients experience and what barriers these changes might pose to the development of
effective vaccines and other antiviral therapeutics. It is reasonable to believe that the responsible
laboratories under the control of the CCP have been engaged in this research for a long period of time and
therefore keep in possession a considerable amount of concealed knowledge of SARS-CoV-2. Some of
the knowledge may provide answers to questions that need to be addressed urgently in the global combat
against COVID-19. Such hidden knowledge ought to be made available to the world immediately.
What also need to be held accountable are the individuals and groups within certain organizations and
institutions in the fields of public health and academic research, who knowingly and collaboratively
27
facilitated the CCP’s misinformation campaign and misled the world. On January 18th and 19th, 2020, Dr.
Li-Meng Yan, then anonymously, first revealed that SARS-CoV-2 is of a laboratory origin78,116.
Immediately afterwards, on January 20th, Dr. Zhengli Shi submitted her manuscript to Nature and reported
the first fabricated virus, RaTG134. Since then, many virus fabrications have taken place and all of them
were published as peer-reviewed articles on top scientific journals4-9. Subsequently, based on such reports,
influential opinion articles promoting the natural origin theory have then been published by prominent
scientists and international organizations on such and other high-profile platforms10,117-120.
In contrast to the rigorous promotion of the natural origin theory, strict censorship has been placed by
these and other journals on manuscripts discussing a possible laboratory origin of SARS-CoV-218,121. Our
earlier report1, which was one of such manuscripts and published as a preprint article, also faced
unfounded criticisms dressed as unbiased peer reviews from two groups of scientists led by Drs. Robert
Gallo and Nancy Connell, respectively122,123 (our point-to-point responses are being prepared and will be
published soon). As a result of this collaborative efforts, the public has been largely removed from the
truth about COVID-19 and SARS-CoV-2, which has led to misjudgments, delayed actions, and greater
sufferings of the global community. It is imperative to investigate the scientists, laboratories, institutions,
and relevant collaborators responsible for the creation of SARS-CoV-2 and for the fabrications/cover-up.
It is also imperative to investigate the relevant individuals in the WHO, at the relevant scientific journals,
in the relevant funding agencies, and in other relevant bodies, which have facilitated the creation of SARS-
CoV-2 and the scientific cover-up of its true origin while under full awareness of the nature of these
operations. Finally, it also needs to be investigated which ones of the scientists engaged in the promotion
of the natural origin theory were purely misled by the scientific fraud and which ones were colluding with
the CCP government.
The time has come that the world faces the truth of COVID-19 and takes actions to save the future of
humanity.
Acknowledgements
We thank Daoyu Zhang for sharing with us the observation of abnormal distribution of non-
synonymous mutations between RaTG13 Spike and SARS-CoV-2. We thank Francisco de Asis for
revealing the filenames of the raw sequencing reads for RaTG13. We also thank other individuals,
including anonymous scientists, for uncovering various facts associated with the origin of SARS-CoV-2.
References
1. Yan, L.-M., Kang, S., Guan, J. & Hu, S. Unusual Features of the SARS-CoV-2 Genome Suggesting
Sophisticated Laboratory Modification Rather Than Natural Evolution and Delineation of Its Probable
Synthetic Route. Zenodo.org (preprint), http://doi.org/10.5281/zenodo.4028830 (2020).
2. Wu, F. et al. A new coronavirus associated with human respiratory disease in China. Nature 579, 265-269
(2020).
3. Hu, D. et al. Genomic characterization and infectivity of a novel SARS-like coronavirus in Chinese bats.
Emerg Microbes Infect 7, 154 (2018).
4. Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature
579, 270–273 (2020).
5. Lam, T.T. et al. Identifying SARS-CoV-2-related coronaviruses in Malayan pangolins. Nature (2020).
28
6. Liu, P. et al. Are pangolins the intermediate host of the 2019 novel coronavirus (SARS-CoV-2)? PLoS
Pathog 16, e1008421 (2020).
7. Xiao, K. et al. Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins. Nature (2020).
8. Zhang, T., Wu, Q. & Zhang, Z. Probable Pangolin Origin of SARS-CoV-2 Associated with the COVID-
19 Outbreak. Curr Biol 30, 1578 (2020).
9. Zhou, H. et al. A Novel Bat Coronavirus Closely Related to SARS-CoV-2 Contains Natural Insertions at
the S1/S2 Cleavage Site of the Spike Protein. Curr Biol 30, 2196-2203 e3 (2020).
10. Andersen, K.G., Rambaut, A., Lipkin, W.I., Holmes, E.C. & Garry, R.F. The proximal origin of SARS-
CoV-2. Nat Med 26, 450-452 (2020).
11. Boni, M.F. et al. Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the
COVID-19 pandemic. Nat Microbiol (2020).
12. Hu, B., Guo, H., Zhou, P. & Shi, Z. Characteristics of SARS-CoV-2 and COVID-19. Nature Reviews
Microbiology, https://doi.org/10.1038/s41579-020-00459-7 (2020).
13. Chan, Y.A. & Zhan, S.H. Single source of pangolin CoVs with a near identical Spike RBD to SARS-
CoV-2. bioRxiv, https://doi.org/10.1101/2020.07.07.184374 (2020).
14. Hassanin, A. The SARS-CoV-2-like virus found in captive pangolins from Guangdong should be better
sequenced. bioRxiv, https://doi.org/10.1101/2020.05.07.077016 (2020).
15. Lin, X. & Chen, S. Major Concerns on the Identification of Bat Coronavirus Strain RaTG13 and Quality
of Related Nature Paper. Preprints, 2020060044 (2020).
16. Rahalkar, M. & Bahulikar, R. The Abnormal Nature of the Fecal Swab Sample used for NGS Analysis of
RaTG13 Genome Sequence Imposes a Question on the Correctness of the RaTG13 Sequence.
Preprints.org, 2020080205 (2020).
17. Rahalkar, M.C. & Bahulikar, R.A. Understanding the Origin of ‘BatCoVRaTG13’, a Virus Closest to
SARS-CoV-2. Preprints, 2020050322 (2020).
18. Segreto, R. & Deigin, Y. Is considering a genetic-manipulation origin for SARS-CoV-2 a conspiracy
theory that must be censored? Preprint (Researchgate) DOI: 10.13140/RG.2.2.31358.13129/1 (2020).
19. Seyran, M. et al. Questions concerning the proximal origin of SARS-CoV-2. J Med Virol (2020).
20. Zhang, D. The Pan-SL-CoV/GD sequences may be from contamination. Preprint (zenodo.org), DOI:
10.5281/zenodo.3885333 (2020).
21. Zhang, D. Anomalies in BatCoV/RaTG13 sequencing and provenance. Preprint (zenodo.org),
https://zenodo.org/record/3987503#.Xz9GzC-z3GI (2020).
22. Lab That First Shared Novel Coronavirus Genome Still Shut Down by Chinese Government. Global
Biodefense, https://globalbiodefense.com/headlines/chinese-lab-that-first-shared-novel-coronavirus-
genome-shut-down/ (2020).
23. CGTN Exclusive: Director of Wuhan Institute of Virology says 'let science speak'. CGTN,
https://news.cgtn.com/news/2020-05-23/Exclusive-with-head-of-Wuhan-Institute-of-Virology-Let-
science-speak-QJeOjOZt4Y/index.html (2020).
24. Cohen, J. Wuhan coronavirus hunter Shi Zhengli speaks out. Science,
https://science.sciencemag.org/content/369/6503/487 (2020).
25. Hu, B. et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into
the origin of SARS coronavirus. PLoS Pathog 13, e1006698 (2017).
26. Ge, X.Y. et al. Coexistence of multiple coronaviruses in several bat colonies in an abandoned mineshaft.
Virol Sin 31, 31-40 (2016).
27. Courtney-Guy, S. Chinese scientists ‘found closest relative of coronavirus seven years ago’.
https://metro.co.uk/2020/07/05/chinese-scientists-found-closest-relative-coronavirus-seven-years-ago-
12948668/ (2020).
28. Wu, Z. et al. Novel Henipa-like virus, Mojiang Paramyxovirus, in rats, China, 2012. Emerg Infect Dis 20,
1064-6 (2014).
29. Qiu, J. How China’s ‘Bat Woman’ Hunted Down Viruses from SARS to the New Coronavirus. Scientific
American, https://www.scientificamerican.com/article/how-chinas-bat-woman-hunted-down-viruses-
from-sars-to-the-new-coronavirus1/ (2020).
30. Li, X. The Analysis of Six Patients With Severe Pneumonia Caused By Unknown Viruses. Master’s
Thesis, https://www.documentcloud.org/documents/6981198-Analysis-of-Six-Patients-With-Unknown-
Viruses.html (2013).
29
31. Huang, C. Novel Virus Discovery in Bat and the Exploration of Receptor of Bat Coronavirus HKU9. PhD
Dissertation (in Chinese), National Institute for Viral Disease Control and Prevention, Chinese Center for
Disease Control and Prevention (2016).
32. Ge, X.Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor.
Nature 503, 535-8 (2013).
33. Names of the RaTG13 Amplicon Sequences. https://graph.org/RaTG13-Amplicon-Names-07-03 (2020).
34. Zeng, L.P. et al. Bat Severe Acute Respiratory Syndrome-Like Coronavirus WIV1 Encodes an Extra
Accessory Protein, ORFX, Involved in Modulation of the Host Immune Response. J Virol 90, 6573-6582
(2016).
35. Shang, J. et al. Structural basis of receptor recognition by SARS-CoV-2. Nature 581, 221-224 (2020).
36. Ren, W. et al. Difference in receptor usage between severe acute respiratory syndrome (SARS)
coronavirus and SARS-like coronavirus of bat origin. J Virol 82, 1899-907 (2008).
37. Yang, X.L. et al. Isolation and Characterization of a Novel Bat Coronavirus Closely Related to the Direct
Progenitor of Severe Acute Respiratory Syndrome Coronavirus. J Virol 90, 3253-6 (2015).
38. Luo, C.M. et al. Discovery of Novel Bat Coronaviruses in South China That Use the Same Receptor as
Middle East Respiratory Syndrome Coronavirus. J Virol 92, DOI: 10.1128/JVI.00116-18 (2018).
39. Wang, N. et al. Characterization of a New Member of Alphacoronavirus with Unique Genomic Features
in Rhinolophus Bats. Viruses 11, https://doi.org/10.3390/v11040379 (2019).
40. Korber, B. HIV Signature and Sequence Variation Analysis. Computational Analysis of HIV Molecular
Sequences Chapter 4, 55-72 (2000).
41. Shijulal Nelson-Sathi, P.U., E Sreekumar, R Radhakrishnan Nair, Iype Joseph, Sai Ravi Chandra Nori,
Jamiema Sara Philip, Roshny Prasad, KV Navyasree, Shikha Ramesh, Heera Pillai, Sanu Ghosh, TR
Santosh Kumar, M. Radhakrishna Pillai. Structural and Functional Implications of Non-synonymous
Mutations in the Spike protein of 2,954 SARS-CoV-2 Genomes. bioRxiv,
https://doi.org/10.1101/2020.05.02.071811 (2020).
42. Li, X. et al. Emergence of SARS-CoV-2 through Recombination and Strong Purifying Selection. Science
Advances 6, eabb9153 (2020).
43. Wang, H., Pipes, L. & Nielsen, R. Synonymous mutations and the molecular evolution of SARS-Cov-2
origins. bioRxiv, https://doi.org/10.1101/2020.04.20.052019 (2020).
44. Tang, X. et al. On the origin and continuing evolution of SARS-CoV-2. National Science Review 7,
1012–1023 (2020).
45. Mou, H. et al. Mutations from bat ACE2 orthologs markedly enhance ACE2-Fc neutralization of SARS-
CoV-2. bioRxiv, https://doi.org/10.1101/2020.06.29.178459 (2020).
46. Latham, J. & Wilson, A. A Proposed Origin for SARS-CoV-2 and the COVID-19 Pandemic. Independent
Science News, https://www.independentsciencenews.org/commentaries/a-proposed-origin-for-sars-cov-2-
and-the-covid-19-pandemic/ (2020).
47. Liu, Z.L. et al. Antibody Profiles in Mild and Severe Cases of COVID-19. Clin Chem 66, 1102-1104
(2020).
48. Hicks, J. et al. Serologic cross-reactivity of SARS-CoV-2 with endemic and seasonal Betacoronaviruses.
medXriv, https://doi.org/10.1101/2020.06.22.20137695 (2020).
49. Anand, P., Puranik, A., Aravamudan, M., Venkatakrishnan, A.J. & Soundararajan, V. SARS-CoV-2
strategically mimics proteolytic activation of human ENaC. Elife 9, e58603 (2020).
50. Vallet, V., Chraibi, A., Gaeggeler, H.P., Horisberger, J.D. & Rossier, B.C. An epithelial serine protease
activates the amiloride-sensitive sodium channel. Nature 389, 607-10 (1997).
51. Claas, E.C. et al. Human influenza A H5N1 virus related to a highly pathogenic avian influenza virus.
Lancet 351, 472-7 (1998).
52. Ito, T. et al. Generation of a highly pathogenic avian influenza A virus from an avirulent field isolate by
passaging in chickens. J Virol 75, 4439-43 (2001).
53. Watanabe, R. et al. Entry from the cell surface of severe acute respiratory syndrome coronavirus with
cleaved S protein as revealed by pseudotype virus bearing cleaved S protein. J Virol 82, 11985-91 (2008).
54. Belouzard, S., Chu, V.C. & Whittaker, G.R. Activation of the SARS coronavirus spike protein via
sequential proteolytic cleavage at two distinct sites. Proc Natl Acad Sci U S A 106, 5871-6 (2009).
55. Sun, X., Tse, L.V., Ferguson, A.D. & Whittaker, G.R. Modifications to the hemagglutinin cleavage site
control the virulence of a neurotropic H1N1 influenza virus. J Virol 84, 8683-90 (2010).
30
56. Kido, H. et al. Role of host cellular proteases in the pathogenesis of influenza and influenza-induced
multiple organ failure. Biochim Biophys Acta 1824, 186-94 (2012).
57. Cheng, J. et al. The S2 Subunit of QX-type Infectious Bronchitis Coronavirus Spike Protein Is an
Essential Determinant of Neurotropism. Viruses 11, doi: 10.3390/v11100972 (2019).
58. Hoffmann, M., Kleine-Weber, H. & Pohlmann, S. A Multibasic Cleavage Site in the Spike Protein of
SARS-CoV-2 Is Essential for Infection of Human Lung Cells. Mol Cell 78, 779-784 e5 (2020).
59. Liu, P., Chen, W. & Chen, J.P. Viral Metagenomics Revealed Sendai Virus and Coronavirus Infection of
Malayan Pangolins (Manis javanica). Viruses 11, doi: 10.3390/v11110979 (2019).
60. Tommy Tsan-Yuk Lam, M.H.-H.S., Hua-Chen Zhu, Yi-Gang Tong, Xue-Bing Ni, Yun-Shi Liao, Wei
Wei, William Yiu-Man Cheung, Wen-Juan Li, Lian-Feng Li, Gabriel M Leung, Edward C. Holmes, Yan-
Ling Hu, Yi Guan. Identification of 2019-nCoV related coronaviruses in Malayan pangolins in southern
China. bioRxiv, doi.org/10.1101/2020.02.13.945485 (2020).
61. University, S.C.A. 华南农大发现穿山甲为新型冠状病毒潜在中间宿主 (South China Agricultural
University Found that Pangolins Are The Possible Intermediate Host of SARS-CoV-2). www.edu.cn,
http://www.edu.cn/ke_yan_yu_fa_zhan/gao_xiao_cheng_guo/gao_xiao_zi_xun/202002/t20200207_1710
427.shtml (2020).
62. 华南农业大学:穿山甲为新型冠状病毒潜在中间宿主 (South China Agricultural University:
Pangolins Are The Possible Intermediate Host of SARS-CoV-2). IFENG NEWS,
https://news.ifeng.com/c/7tr8u2sAQFc (2020).
63. Kangpeng Xiao, J.Z., Yaoyu Feng, Niu Zhou, Xu Zhang, Jie-Jian Zou, Na Li, Yaqiong Guo, Xiaobing Li,
Xuejuan Shen, Zhipeng Zhang, Fanfan Shu, Wanyi Huang, Yu Li, Ziding Zhang, Rui-Ai Chen, Ya-Jiang
Wu, Shi-Ming Peng, Mian Huang, Wei-Jun Xie, Qin-Hui Cai, Fang-Hui Hou, Yahong Liu, Wu Chen,
Lihua Xiao, Yongyi Shen. Isolation and Characterization of 2019-nCoV-like Coronavirus from Malayan
Pangolins. bioRxiv, doi.org/10.1101/2020.02.17.951335 (2020).
64. Bi, S. et al. Complete genome sequences of the SARS-CoV: the BJ Group (Isolates BJ01-BJ04).
Genomics Proteomics Bioinformatics 1, 180-92 (2003).
65. Qin, E. et al. A genome sequence of novel SARS-CoV isolates: the genotype, GD-Ins29, leads to a
hypothesis of viral transmission in South China. Genomics Proteomics Bioinformatics 1, 101-7 (2003).
66. Qin, E. et al. A complete sequence and comparative analysis of a SARS-associated virus (Isolate BJ01).
Chin Sci Bull 48, 941-948 (2003).
67. Zhu, X. et al. Genetic variation of the human alpha-2-Heremans-Schmid glycoprotein (AHSG) gene
associated with the risk of SARS-CoV infection. PLoS One 6, e23730 (2011).
68. Lee, J. et al. No evidence of coronaviruses or other potentially zoonotic viruses in Sunda pangolins
(Manis javanica) entering the wildlife trade via Malaysia. bioRxiv,
https://doi.org/10.1101/2020.06.19.158717 (2020).
69. Piplani, S., Singh, P.K., Winkler, D.A. & Petrovsky, N. In silico comparison of spike protein-ACE2
binding affinities across species; significance for the possible origin of the SARS-CoV-2 virus. arXiv,
arXiv:2005.06199 (2020).
70. Wrobel, A.G. et al. Structure and binding properties of Pangolin-CoV Spike glycoprotein inform the
evolution of SARS-CoV-2. Research Square, DOI: 10.21203/rs.3.rs-83072/v1 (2020).
71. Drexler, J.F. et al. Genomic characterization of severe acute respiratory syndrome-related coronavirus in
European bats and classification of coronaviruses based on partial RNA-dependent RNA polymerase gene
sequences. J Virol 84, 11336-49 (2010).
72. Tao, Y. & Tong, S. Complete Genome Sequence of a Severe Acute Respiratory Syndrome-Related
Coronavirus from Kenyan Bats. Microbiol Resour Announc 8(2019).
73. Hou, Y.J. et al. SARS-CoV-2 Reverse Genetics Reveals a Variable Infection Gradient in the Respiratory
Tract. Cell 182, 429-446 e14 (2020).
74. 生物恐怖袭击——并非杞人忧天 (Bioterrarism Attack — Not An Unfounded Concern). 保健时报
(Health Times), http://www.bjsbnet.com/Article/InArticle/BJSB200505120031 (2005).
75. COVID-19 Overview and Infection Prevention and Control Priorities in non-US Healthcare Settings.
cdc.org, https://www.cdc.gov/coronavirus/2019-ncov/hcp/non-us-
settings/overview/index.html#:~:text=COVID%2D19%20is%20primarily,inhaled%20into%20the%20lun
gs (2020).
31
76. Sia, S.F. et al. Pathogenesis and transmission of SARS-CoV-2 in golden hamsters. Nature (2020).
77. Scientific Brief: SARS-CoV-2 and Potential Airborne Transmission. cdc.org,
https://www.cdc.gov/coronavirus/2019-ncov/more/scientific-brief-sars-cov-2.html (2020).
78. The Morning Show (with An Hong and Ai Li) on Jan 19th, 2020. Lude Press (YouTube),
https://www.youtube.com/watch?v=CLTjg03CPEs (2020).
79. Chin, A.W.H. et al. Stability of SARS-CoV-2 in different environmental conditions. Lancet Microbe 1,
e10 (2020).
80. van Doremalen, N. et al. Aerosol and Surface Stability of SARS-CoV-2 as Compared with SARS-CoV-1.
N Engl J Med 382, 1564-1567 (2020).
81. Eroshenko, N. et al. Implications of antibody-dependent enhancement of infection for SARS-CoV-2
countermeasures. Nat Biotechnol 38, 789-791 (2020).
82. Bowen, T. Nevada State Public Health Lab-led team studying COVID-19 reinfection. News & Events
(University of Nevada, Reno, School of Medicine), https://med.unr.edu/news/archive/2020/covid-19-
reinfection (2020).
83. Kam, Y.W. et al. Antibodies against trimeric S glycoprotein protect hamsters against SARS-CoV
challenge despite their capacity to mediate FcgammaRII-dependent entry into B cells in vitro. Vaccine 25,
729-40 (2007).
84. Ren, W. et al. Full-length genome sequences of two SARS-like coronaviruses in horseshoe bats and
genetic variation analysis. J Gen Virol 87, 3355-9 (2006).
85. Yuan, J. et al. Intraspecies diversity of SARS-like coronaviruses in Rhinolophus sinicus and its
implications for the origin of SARS coronaviruses in humans. J Gen Virol 91, 1058-62 (2010).
86. Ge, X.Y. et al. Detection of alpha- and betacoronaviruses in rodents from Yunnan, China. Virol J 14, 98
(2017).
87. Luo, Y. et al. Longitudinal Surveillance of Betacoronaviruses in Fruit Bats in Yunnan Province, China
During 2009-2016. Virol Sin 33, 87-95 (2018).
88. Wang, Y. Preliminary investigation of viruses carried by bats on the southeast coastal area (东南沿海地
区蝙蝠携带病毒的初步调查研究). Master’s Thesis (2017).
89. Wu, Z. et al. Detection of Hantaviruses and Arenaviruzses in three-toed jerboas from the Inner Mongolia
Autonomous Region, China. Emerg Microbes Infect 7, 35 (2018).
90. Menachery, V.D. et al. A SARS-like cluster of circulating bat coronaviruses shows potential for human
emergence. Nat Med 21, 1508-13 (2015).
91. Yang, X.L. et al. Characterization of a filovirus (Mengla virus) from Rousettus bats in China. Nat
Microbiol 4, 390-395 (2019).
92. Waruhiu, C. et al. Molecular detection of viruses in Kenyan bats and discovery of novel astroviruses,
caliciviruses and rotaviruses. Virol Sin 32, 101-114 (2017).
93. Ge, X.Y. et al. Fugong virus, a novel hantavirus harbored by the small oriental vole (Eothenomys eleusis)
in China. Virol J 13, 27 (2016).
94. Yang, X.L. et al. Isolation and identification of bat viruses closely related to human, porcine and mink
orthoreoviruses. J Gen Virol 96, 3525-3531 (2015).
95. Hu, B. et al. Detection of diverse novel astroviruses from small mammals in China. J Gen Virol 95, 2442-
2449 (2014).
96. Yang, X., Zhang, Y., Ge, X., Yuan, J. & Shi, Z. A novel totivirus-like virus isolated from bat guano. Arch
Virol 157, 1093-9 (2012).
97. Zhao, L. et al. Characterization of a Novel Tanay Virus Isolated From Anopheles sinensis Mosquitoes in
Yunnan, China. Front Microbiol 10, 1963 (2019).
98. Xia, H. et al. First Isolation and Characterization of a Group C Banna Virus (BAV) from Anopheles
sinensis Mosquitoes in Hubei, China. Viruses 10, doi: 10.3390/v10100555 (2018).
99. Wang, Y., Xia, H., Zhang, B., Liu, X. & Yuan, Z. Isolation and characterization of a novel mesonivirus
from Culex mosquitoes in China. Virus Res 240, 130-139 (2017).
100. Li, L.L. et al. Detection and characterization of a novel hepacivirus in long-tailed ground squirrels
(Spermophilus undulatus) in China. Arch Virol 164, 2401-2410 (2019).
101. Zhou, Z. et al. Complete genome sequences of two crimean-congo hemorrhagic Fever viruses isolated in
china. Genome Announc 1, doi: 10.1128/genomeA.00571-13 (2013).
32
102. Lau, S.K. et al. Complete genome analysis of three novel picornaviruses from diverse bat species. J Virol
85, 8819-28 (2011).
103. Woo, P.C. et al. Discovery of seven novel Mammalian and avian coronaviruses in the genus
deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus
and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus. J Virol 86, 3995-
4008 (2012).
104. Zuo, S.Q. et al. A new hantavirus from the stripe-backed shrew (Sorex cylindricauda) in the People's
Republic of China. Virus Res 184, 82-6 (2014).
105. Zuo, S. et al. Detection of Quang Binh virus from mosquitoes in China. Virus Res 180, 31-8 (2014).
106. He, B. et al. Group A Rotaviruses in Chinese Bats: Genetic Composition, Serology, and Evidence for Bat-
to-Human Transmission and Reassortment. J Virol 91, doi: 10.1128/JVI.02493-16 (2017).
107. Feng, Y. et al. Isolation and full-length genome analysis of mosquito-borne Manzanilla virus from
Yunnan Province, China. BMC Res Notes 8, 255 (2015).
108. Xu, L. et al. Novel hantavirus identified in black-bearded tomb bats, China. Infect Genet Evol 31, 158-60
(2015).
109. He, B. et al. Hepatitis virus in long-fingered bats, Myanmar. Emerg Infect Dis 19, 638-40 (2013).
110. Sun, H. et al. Prevalence and genetic characterization of Toxoplasma gondii in bats in Myanmar. Appl
Environ Microbiol 79, 3526-8 (2013).
111. Yang, R. Progress in research concerning Yersinia pestis and its significance in military medicine. in
semanticscholar.org https://www.semanticscholar.org/paper/Progress-in-research-concerning-Yersinia-
pestis-and-Rui-fu/69c3be6d683dce8e992086d8c92c8119c039260c (2012).
112. He, B. et al. Characterization of a novel G3P[3] rotavirus isolated from a lesser horseshoe bat: a distant
relative of feline/canine rotaviruses. J Virol 87, 12357-66 (2013).
113. He, B. et al. Identification of a novel Orthohepadnavirus in pomona roundleaf bats in China. Arch Virol
160, 335-7 (2015).
114. Hu, T. et al. Characterization of a novel orthoreovirus isolated from fruit bat, China. BMC Microbiol 14,
293 (2014).
115. Becker, M.M. et al. Synthetic recombinant bat SARS-like coronavirus is infectious in cultured cells and
in mice. Proc Natl Acad Sci U S A 105, 19944-9 (2008).
116. LuDe. Twitter, https://twitter.com/ding_gang/status/1218547052084441088 (2020).
117. Calisher, C. et al. Statement in support of the scientists, public health professionals, and medical
professionals of China combatting COVID-19. Lancet 395, e42-e43 (2020).
118. Liu, S.L., Saif, L.J., Weiss, S.R. & Su, L. No credible evidence supporting claims of the laboratory
engineering of SARS-CoV-2. Emerg Microbes Infect 9, 505-507 (2020).
119. Collins, F. Genomic Study Points to Natural Origin of COVID-19. NIH Director’s Blog,
https://directorsblog.nih.gov/2020/03/26/genomic-research-points-to-natural-origin-of-covid-19/ (2020).
120. WHO. Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19).
https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-
report.pdf (2020).
121. Robinson, C. Journals censor lab origin theory for SARS-CoV-2.
(https://www.gmwatch.org/en/news/latest-news/19475-journals-censor-lab-origin-theory-for-sars-cov-2,
2020).
122. Warmbrod, K.L., West, R.M., Connell, N.D. & Gronvall, G.K. In Response: Yan et al Preprint
Examinations of the Origin of SARS-CoV-2. John Hopkins Center for Health Security,
https://www.centerforhealthsecurity.org/our-work/pubs_archive/pubs-pdfs/2020/200921-in-response-
yan.pdf (2020).
123. Koyama, T., Lauring, A., Gallo, R. & Reitz, M. Reviews of "Unusual Features of the SARS-CoV-2
Genome Suggesting Sophisticated Laboratory Modification Rather Than Natural Evolution and
Delineation of Its Probable Synthetic Route". Rapid Reviews COVID-19,
https://rapidreviewscovid19.mitpress.mit.edu/pub/78we86rp/release/2 (2020).
33