Assessing the emergence time of SARS-CoV-2 zoonotic spillover.

Understanding the evolution of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) and its relationship to other coronaviruses in the wild is crucial for preventing future virus outbreaks. While the origin of the SARS-CoV-2 pandemic remains uncertain, mounting evidence suggests the direct involvement of the bat and pangolin coronaviruses in the evolution of the SARS-CoV-2 genome. To unravel the early days of a probable zoonotic spillover event, we analyzed genomic data from various coronavirus strains from both human and wild hosts. Bayesian phylogenetic analysis was performed using multiple datasets, using strict and relaxed clock evolutionary models to estimate the occurrence times of key speciation, gene transfer, and recombination events affecting the evolution of SARS-CoV-2 and its closest relatives. We found strong evidence supporting the presence of temporal structure in datasets containing SARS-CoV-2 variants, enabling us to estimate the time of SARS-CoV-2 zoonotic spillover between August and early October 2019. In contrast, datasets without SARS-CoV-2 variants provided mixed results in terms of temporal structure. However, they allowed us to establish that the presence of a statistically robust clade in the phylogenies of gene S and its receptor-binding (RBD) domain, including two bat (BANAL) and two Guangdong pangolin coronaviruses (CoVs), is due to the horizontal gene transfer of this gene from the bat CoV to the pangolin CoV that occurred in the middle of 2018. Importantly, this clade is closely located to SARS-CoV-2 in both phylogenies. This phylogenetic proximity had been explained by an RBD gene transfer from the Guangdong pangolin CoV to a very recent ancestor of SARS-CoV-2 in some earlier works in the field before the BANAL coronaviruses were discovered. Overall, our study provides valuable insights into the timeline and evolutionary dynamics of the SARS-CoV-2 pandemic.

Evolution and conservation genetics of pangolins.

Pangolins (Pholidota, Manidae) are classified as an evolutionarily distinct and globally endangered mammal due to their unique morphology (nail-like scales and a myrmecophagous diet) and being the victim of heavy poaching and worldwide trafficking. As such, pangolins serve as a textbook example for studying the special phenotypic evolutionary adaptations and conservation genetics of an endangered species. Recent years have demonstrated significant advancements in the fields of molecular genetics and genomics, which have translated to a series of important research achievements and breakthroughs concerning the evolution and conservation genetics of pangolins. This review comprehensively presents the hitherto advances in phylogeny, adaptive evolution, conservation genetics, and conservation genomics that are related to pangolins, which will provide an ample understanding of their diversity, molecular adaptation mechanisms, and evolutionary potentials. In addition, we highlight the priority of investigating species/population diversity among pangolins and suggest several avenues of research that are highly relevant for future pangolin conservation.

Visualisation and detection of latent DNA deposited by pangolin scales onto plastic packaging materials.

We report on the detection and visualisation of latent DNA from pangolin scales deposited onto a plastic packaging material through the use of a nucleic acid staining dye. This latent DNA deposited by pangolin scales was subsequently isolated and analysed using DNA barcoding method. Pangolins are the most illegally traded mammalian species due to the demand for their scales and meat. The demand for their scales were mostly fuelled by its use in traditional medicines. The scales are usually packed into bags and transported globally via sea routes. This is the first report detailing the detection of trace latent DNA from processed wildlife products, on surfaces of bags that they were packaged in. Prior to this report, it was not known if the dried pangolin scales contained transferable quantities of biological material for DNA analyses. To address this, scales were removed from a roadkill Sunda pangolin (Manis javanica), processed by drying and packaged into one of five plastic bags. The presence of pangolin latent DNA was detected on the surface of the plastic bags and visualised using Diamond™ nucleic acid dye. Swabs were then used to recover the stained biological material from various locations in the five bags. The DNA was isolated and quantified using a newly designed quantitative PCR (qPCR) specific to M. javanica to amplify a fragment of the mitochondrial DNA cytochrome b gene. There was a positive correlation between the number of stained particles and DNA quantity, and a greater number of stained particles were found at the bottom of the bag than were found at the top. Conventional PCR targeting part of the cyt b gene amplified a product from all 30 samples taken from the bags and in all cases, sequence data generated matched that of the Sunda pangolin, as expected. All negative controls yielded no results. The method described here is the very first use of a nucleic acid staining dye to detect latent DNA from a mammalian species, other than humans, and highlights the opportunity for further use of Diamond™ nucleic acid dye in wildlife forensic science. It is anticipated that this method will be invaluable in retrieving latent DNA deposited by wildlife products from the environment in which they were contained, to determine the presence of these illegal wildlife products even when previously hidden, inaccessible, or no longer present physically. Further research is required to understand if the use on non-human mammalian wildlife species is feasible.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays.

BACKGROUND: Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes. RESULTS: We benchmark eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compare experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, is lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieve the best overall performance at distinguishing disruptive and neutral variants, and controlling for overall call rate genome-wide, SpliceAI and Pangolin have superior sensitivity. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues. CONCLUSION: SpliceAI and Pangolin show the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals.

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

Genome sequence analysis of Malayan pangolin (Manis javanica) forensic samples reveals the presence of Paraburkholderia fungorum sequences.

Background: The Malayan pangolin (Manis javanica) is a placental mammal and is listed as Critically Endangered on the IUCN Red List of Threatened Species. Most previous attempts to breed pangolins in captivity have met with little success because of dietary issues, infections, and other complications, although a previous study reported breeding pangolins in captivity to the third generation. In our previous pangolin genome sequencing data analysis, we obtained a considerable amount of bacterial DNA from a pregnant female Malayan pangolin (named "UM3"), which was likely infected by Paraburkholderia fungorum-an agent of biodegradation and bioremediation in agriculture. Methodology: Here, we further confirmed and characterized this bacterial species using PCR, histological staining, whole-genome sequencing, and bioinformatics approaches. PCR assays with in-house designed primer sets and 16S universal primers showed clear positive bands in the cerebrum, cerebellum, lung, and blood of UM3 suggesting that UM3 might have developed septicaemia. Histological staining showed the presence of Gram-negative rod-shaped bacteria in the pangolin brain and lungs, indicating the colonization of the bacteria in these two organs. In addition, PCR screening of UM3's fetal tissues revealed the presence of P. fungorum in the gastrocnemius muscle, but not in other tissues that we examined. We also sequenced and reconstructed the genome of pangolin P. fungorum, which has a genome size of 7.7 Mbps. Conclusion: Our study is the first to present detailed evidence of the presence of P. fungorum in a pangolin and her fetus (although preliminary results were presented in our previous article). Here, we raise the concern that P. fungorum may potentially infect humans, especially YOPI (young, old, pregnant, and immunocompromised) people. Therefore, caution should be exercised when using this bacterial species as biodegradation or bioremediation agents in agriculture.

Genomic analysis reveals a cryptic pangolin species.

Eight extant species of pangolins are currently recognized. Recent studies found that two mitochondrial haplotypes identified in confiscations in Hong Kong could not be assigned to any known pangolin species, implying the existence of a species. Here, we report that two additional mitochondrial haplotypes identified in independent confiscations from Yunnan align with the putative species haplotypes supporting the existence of this mysterious species/population. To verify the new species scenario we performed a comprehensive analysis of scale characteristics and 138 whole genomes representing all recognized pangolin species and the cryptic new species, 98 of which were generated here. Our morphometric results clearly attributed this cryptic species to Asian pangolins (Manis sp.) and the genomic data provide robust and compelling evidence that it is a pangolin species distinct from those recognized previously, which separated from the Philippine pangolin and Malayan pangolin over 5 Mya. Our study provides a solid genomic basis for its formal recognition as the ninth pangolin species or the fifth Asian one, supporting a new taxonomic classification of pangolins. The effects of glacial climate changes and recent anthropogenic activities driven by illegal trade are inferred to have caused its population decline with the genomic signatures showing low genetic diversity, a high level of inbreeding, and high genetic load. Our finding greatly expands current knowledge of pangolin diversity and evolution and has vital implications for conservation efforts to prevent the extinction of this enigmatic and endangered species from the wild.

Pathogenicity, tissue tropism and potential vertical transmission of SARSr-CoV-2 in Malayan pangolins.

Malayan pangolin SARS-CoV-2-related coronavirus (SARSr-CoV-2) is closely related to SARS-CoV-2. However, little is known about its pathogenicity in pangolins. Using CT scans we show that SARSr-CoV-2 positive Malayan pangolins are characterized by bilateral ground-glass opacities in lungs in a similar manner to COVID-19 patients. Histological examination and blood gas tests are indicative of dyspnea. SARSr-CoV-2 infected multiple organs in pangolins, with the lungs the major target, and histological expression data revealed that ACE2 and TMPRSS2 were co-expressed with viral RNA. Transcriptome analysis indicated that virus-positive pangolins were likely to have inadequate interferon responses, with relative greater cytokine and chemokine activity in the lung and spleen. Notably, both viral RNA and viral proteins were detected in three pangolin fetuses, providing initial evidence for vertical virus transmission. In sum, our study outlines the biological framework of SARSr-CoV-2 in pangolins, revealing striking similarities to COVID-19 in humans.

Decay of Skin-Specific Gene Modules in Pangolins.

The mammalian skin exhibits a rich spectrum of evolutionary adaptations. The pilosebaceous unit, composed of the hair shaft, follicle, and the sebaceous gland, is the most striking synapomorphy. The evolutionary diversification of mammals across different ecological niches was paralleled by the appearance of an ample variety of skin modifications. Pangolins, order Pholidota, exhibit keratin-derived scales, one of the most iconic skin appendages. This formidable armor is intended to serve as a deterrent against predators. Surprisingly, while pangolins have hair on their abdomens, the occurrence of sebaceous and sweat glands is contentious. Here, we explore various molecular modules of skin physiology in four pangolin genomes, including that of sebum production. We show that genes driving wax monoester formation, Awat1/2, show patterns of inactivation in the stem pangolin branch, while the triacylglycerol synthesis gene Dgat2l6 seems independently eroded in the African and Asian clades. In contrast, Elovl3 implicated in the formation of specific neutral lipids required for skin barrier function is intact and expressed in the pangolin skin. An extended comparative analysis shows that genes involved in skin pathogen defense and structural integrity of keratinocyte layers also show inactivating mutations: associated with both ancestral and independent pseudogenization events. Finally, we deduce that the suggested absence of sweat glands is not paralleled by the inactivation of the ATP-binding cassette transporter Abcc11, as previously described in Cetacea. Our findings reveal the sophisticated and complex history of gene retention and loss as key mechanisms in the evolution of the highly modified mammalian skin phenotypes.

Comparative analyses and phylogenetic relationships of thirteen Pholidota species (Orchidaceae) inferred from complete chloroplast genomes.

BACKGROUND: The orchid genus Pholidota Lindl. ex Hook. is economically important as some species has long been used in traditional medicine. However, the systematic status of the genus and intergeneric relationships inferred from previous molecular studies are unclear due to insufficient sampling and lack of informative sites. So far, only limited genomic information has been available. The taxonomy of Pholidota remains unresolved and somewhat controversial. In this study, the complete chloroplast (cp.) genomes of thirteen Pholidota species were sequenced and analyzed to gain insight into the phylogeny of Pholidota and mutation patterns in their cp. genomes. RESULTS: All examined thirteen Pholidota cp. genomes exhibited typical quadripartite circular structures, with the size ranging from 158,786 to 159,781 bp. The annotation contained a total of 135 genes in each cp. genome, i.e., 89 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The codon usage analysis indicated the preference of A/U-ending codons. Repeat sequence analysis identified 444 tandem repeats, 322 palindromic repeats and 189 dispersed repeats. A total of 525 SSRs, 13,834 SNPs and 8,630 InDels were detected. Six mutational hotspots were identified as potential molecular markers. These molecular markers and highly variable regions are expected to facilitate future genetic and genomic studies. Our phylogenetic analyses confirmed the polyphyletic status of the genus Pholidota, with species grouped into four main clades: Pholidota s.s. was resolved as the sister to a clade containing species of Coelogyne; the other two clades clustered together with species of Bulleyia and Panisea, respectively; species P. ventricosa was placed at the basal position, deviated from all other species. CONCLUSION: This is the first study to comprehensively examine the genetic variations and systematically analyze the phylogeny and evolution of Pholidota based on plastid genomic data. These findings contribute to a better understanding of plastid genome evolution of Pholidota and provide new insights into the phylogeny of Pholidota and its closely related genera within the subtribe Coelogyninae. Our research has laid the foundation for future studies on the evolutionary mechanisms and classification of this economically and medicinally important genus.

Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity.

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.

Characteristics of the multiple replicon plasmid IncX1-X1 in multidrug-resistant Escherichia coli from Malayan pangolin (Manis javanica).

Potential zoonotic pathogens may be transmitted from wildlife to humans through the illegal wild meat trade, which has become a pressing issue. However, research on the antimicrobial resistance genes (ARGs) of Malayan pangolin (Manis javanica) intestinal bacteria is limited. Here, multidrug-resistant Escherichia coli M172-1 (ST354) isolated from Malayan pangolin feces in 2019 was found to be resistant to 13 antibiotics. BGWAS analysis revealed 4 plasmids, namely, pM172-1.1, pM172-1.2, pM172-1.3, and pM172-1.4, in the isolate. The pM172-1.2, pM172-1.3, and pM172-1.4 plasmids carried ARGs, namely, IncHI2-HI2A, IncX1-X1, and IncX1, respectively. pM172-1.3 and pM172-1.4 contained intact IntI1 integrons (Is26/IntI1/arr2/cmlA5/blaOXA-10 /ant(3″)-IIA/dfrA14/Is26). Notably, pM172-1.3 resulted from the fusion of 2 pM172-1.4 copies and carried many more ARGs. In addition to pM172-1.3 from the same host, other drug-resistant bacteria (E. coli M159-1 (ST48), E. coli S171-1 (ST206), and Klebsiella pneumoniae S174-1 (ST2354)) in the same Malayan pangolin fecal samples also carried 3 plasmids with 100% gene coverage of pM172-1.4 and 99.98% identity. Therefore, ARGs in IncX1 might spread in the intestinal flora of Malayan pangolin and between species via the illegal food chain, posing a potential threat to public health and safety.

Role of the Pangolin in Origin of SARS-CoV-2: An Evolutionary Perspective.

After the recent emergence of SARS-CoV-2 infection, unanswered questions remain related to its evolutionary history, path of transmission or divergence and role of recombination. There is emerging evidence on amino acid substitutions occurring in key residues of the receptor-binding domain of the spike glycoprotein in coronavirus isolates from bat and pangolins. In this article, we summarize our current knowledge on the origin of SARS-CoV-2. We also analyze the host ACE2-interacting residues of the receptor-binding domain of spike glycoprotein in SARS-CoV-2 isolates from bats, and compare it to pangolin SARS-CoV-2 isolates collected from Guangdong province (GD Pangolin-CoV) and Guangxi autonomous regions (GX Pangolin-CoV) of South China. Based on our comparative analysis, we support the view that the Guangdong Pangolins are the intermediate hosts that adapted the SARS-CoV-2 and represented a significant evolutionary link in the path of transmission of SARS-CoV-2 virus. We also discuss the role of intermediate hosts in the origin of Omicron.

Molecular Genetics Supports a Potential Fifth Asian Pangolin Species (Mammalia, Pholidota, Manis).

Recently, two mitochondrial haplotypes, H4 and H8, of Manis sp. were found in two seizures in Hong Kong that do not correspond to Manis javanica, Manis pentadactyla or Manis crassicaudata of Asian pangolin species or any African pangolin species. It was proposed that both haplotypes derived from Manis culionensis, an unknown lineage of M. javanica, or a thus far unidentified Asian pangolin species (Manis sp.). To further investigate these three hypotheses, we used two mitochondrial genes of all eight known extant pangolin species and conducted phylogenetic tree reconstructions, divergence time estimation, and species delimitation analyses. All analyses consistently confirmed that these two haplotypes of Manis sp. constitute a distinct lineage, potentially representing a fifth Asian pangolin species, which originated around the Late Miocene to Early Pliocene (6.95 [4.64-9.85] million years ago). Our study provides genetic support for a potential fifth Asian pangolin species and helps to better understand species diversity of Asian pangolins, which is urgently needed for effective conservation work.

Combined proteomics and transcriptomics reveal the genetic basis underlying the differentiation of skin appendages and immunity in pangolin.

Pangolin (Mains javanica) is an interesting endangered mammal with special morphological characteristics. Here, we applied proteomics and transcriptomics to explore the differentiation of pangolin skin appendages at two developmental stages and to compare gene expression profiles between abdomen hair and dorsal scale tissues. We identified 4,311 genes and 91 proteins differentially expressed between scale-type and hair-type tissue, of which 6 genes were shared by the transcriptome and proteome. Differentiation altered the abundance of hundreds of proteins and mRNA in the two types of skin appendages, many of which are involved in keratinocyte differentiation, epidermal cell differentiation, and multicellular organism development based on GO enrichment analysis, and FoxO, MAPK, and p53 signalling pathways based on KEGG enrichment analysis. DEGs in scale-type tissues were also significantly enriched in immune-related terms and pathways compared with that in hair-type tissues. Thus, we propose that pangolins have a normal skin innate immune system. Compared with the abdomen, the back skin of pangolins had more genes involved in the regulation of immune function, which may be an adaptive adjustment for the vulnerability of scaly skin to infection and injury. This investigation provides a scientific basis for the study of development and immunity of pangolin skin, which may be helpful in the protection of wild pangolin in China.

Pangolin Indexing System: implications in forensic surveillance of large seizures.

Demand for pangolin scales in East Asia has increased dramatically in the past two decades, raising concern to the pangolin survival and bringing them to the brink of local extinction. Enumerating the number of individuals from the seized pangolin scales primarily goes undocumented, mostly due to the unavailability of the appropriate methods. In this study, we developed a Pangolin Indexing System, a multi-locus STR panel of eight dinucleotide microsatellites that showed promising results in individualization and assignment of scales into Chinese and Indian pangolins. The combined power of exclusion was 0.83 and 0.99 for Chinese and Indian pangolin. The select panel of eight polymorphic STRs exhibited the cumulative probability of identity 3.7 × 10-9 for Indian pangolin and 3.6 × 10-7 for Chinese pangolin and identified 51 unique genotypes from the 74 scales selected from the four pangolin seizures. The study demonstrated the first report of cross-species validation of STRs developed from Malayan pangolin to Indian pangolin and showed the potential application of Pangolin Indexing System in screening of large seizures through DNA profiling from the scales of Indian and Chinese pangolin.

Cytosolic DNA sensing through cGAS and STING is inactivated by gene mutations in pangolins.

The release of DNA into the cytoplasm upon damage to the nucleus or during viral infection triggers an interferon-mediated defense response, inflammation and cell death. In human cells cytoplasmic DNA is sensed by cyclic GMP-AMP Synthase (cGAS) and Absent In Melanoma 2 (AIM2). Here, we report the identification of a "natural knockout" model of cGAS. Comparative genomics of phylogenetically diverse mammalian species showed that cGAS and its interaction partner Stimulator of Interferon Genes (STING) have been inactivated by mutations in the Malayan pangolin whereas other mammals retained intact copies of these genes. The coding sequences of CGAS and STING1 are also disrupted by premature stop codons and frame-shift mutations in Chinese and tree pangolins, suggesting that expression of these genes was lost in a common ancestor of all pangolins that lived more than 20 million years ago. AIM2 is retained in a functional form in pangolins whereas it is inactivated by mutations in carnivorans, the phylogenetic sister group of pangolins. The deficiency of cGAS and STING points to the existence of alternative mechanisms of controlling cytoplasmic DNA-associated cell damage and viral infections in pangolins.

Molecular evolution and phylogenetic analysis of SARS-CoV-2 and hosts ACE2 protein suggest Malayan pangolin as intermediary host.

An emergence of a novel coronavirus, causative agent of COVID19, named as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), occurred due to cross-species transmission. Coronaviruses are a large family of viruses able to infect a great number of hosts. Entrance of SARS-CoV-2 depends on the surface (S) protein interaction with host ACE2 protein and cleavage by TMPRSS2. ACE2 could be a species-specific barrier that interferes with bat-to-human coronavirus cross-species transmission. Molecular analysis supported bats as natural hosts for SARS-CoV and involved them in MERS-CoV origin. The genomic similarity between bat RaTG13 CoV strain and SARS-CoV-2 implicates bats in the origin of the new outbreak. Additionally, there is a hypothesis for the zoonotic transmission based on contact with Malayan pangolins by humans in Huanan seafood market in Wuhan, China. To investigate bats and pangolin as hosts in SARS-CoV-2 cross-species transmission, we perform an evolutionary analysis combining viral and host phylogenies and divergence of ACE2 and TMPRSS2 amino acid sequences between CoV hosts. Phylogeny showed SARS-like-CoV-2 strains that infected pangolin and bats are close to SARS-CoV-2. In contrast to TMPRSS2, pangolin ACE2 amino acid sequence has low evolutionary divergence compared with humans and is more divergent from bats. Comparing SARS-CoV with SARS-CoV-2 origins, pangolin has yet lower ACE2 evolutionary divergence with humans than civet-the main intermediary host of SARS-CoV. Thus, pangolin has become an opportune host to intermediates bat-to-human SARS-CoV-2 jump and entry.

Isolation and characterization of 30 STRs in Temminck's ground pangolin (Smutsia temminckii) and potential for cross amplification in other African species.

Temminck's ground pangolin (Smutsia temminckii) is one of four species of pangolin, endemic to Africa. Two of the African pangolin species are listed as vulnerable and two are listed as endangered on the International Union for Conservation of Nature Red List of Threatened Species due to their ongoing exploitation for traditional medicine and bushmeat. In this study, we developed 30 species-specific short-tandem repeats (STRs) in Temminck's ground pangolin using next-generation sequencing. The markers were also optimized for crossamplification in other African species. All the markers amplified successfully in Temminck's ground pangolin with allelic polymorphisms observed in 87% of the markers in giant pangolin (S. gigantea) whereas 60% of the markers were amplified polymorphic loci in both whitebellied pangolin (Phataginus tricuspis) and black-bellied pangolin (P. tetradactyla). Analysis of diversity estimates showed moderate levels of variability in Temminck's ground pangolin (Na = 5; Ho = 0.559), giant pangolin (Na = 4.909; Ho = 0.514) and white-bellied pangolin (Na= 2.686; Ho = 0.541) with lower values being observed in black-bellied pangolin (Na = 3; Ho = 0.242). This study provides data of the first available STR markers which was amplified in all four African pangolin species that can now be used in conservation genetic and evolutionary aspects of population histories.