Integrating Fréchet distance and AI reveals the evolutionary trajectory and origin of SARS-CoV-2

A genome, composed of a precisely ordered sequence of four nucleotides (ATCG), encompasses a multitude of specific genome features like AAA motif. Mutations occurring within a genome disrupt the sequential order and composition of these features, thereby influencing the evolutionary trajectories and yielding variants. The evolutionary relatedness between a variant and its ancestor can be estimated by assessing evolutionary distances across a spectrum of genome features. This study develops a novel, alignment-free algorithm that considers both the sequential order and composition of genome features, enabling computation of the Fréchet distance (Fr) across multiple genome features to quantify the evolutionary status of a variant. Integrating this algorithm with an artificial recurrent neural network (RNN) reveals the quantitative evolutionary trajectory and origin of SARS-CoV-2, a puzzle unsolved by alignment-based phylogenetics. The RNN generates the evolutionary trajectory from Fr data at two levels: genome sequence mutations and organism variants. At the genome sequence level, SARS-CoV-2 evolutionarily shortens its genome to enhance its infectious capacity. Mutating signature features, such as TTA and GCT, increases its infectious potential and drives its evolution. At the organism level, variants mutating a single biomarker possess low infectious potential. However, mutating multiple markers dramatically increases their infectious capacity, propelling the COVID-19 pandemic. SARS-CoV-2 likely originates from mink coronavirus variants, with its origin trajectory traced as follows: mink, cat, tiger, mouse, hamster, dog, lion, gorilla, leopard, bat, and pangolin. Together, mutating multiple signature features and biomarkers delineates the evolutionary trajectory of mink-origin SARS-CoV-2, leading to the COVID-19 pandemic.

The signature of SARS-CoV-2 evolution reflects selective pressures within human guts

In somatic cells, microRNAs (miRNAs) bind to the genomes of RNA viruses and influence their translation and replication. In London and Berlin samples represented in GISAID database, we traced severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages and divided these sequenced in two groups, "Ancestral variants" and "Omicrons," and analyzed them through the prism of the tissue-specific binding between host miRNAs and viral messenger RNAs. We demonstrate a significant number of miRNA-binding sites in the NSP4 region of the SARS-CoV-2 genome, with evidence of evolutionary pressure within this region exerted by human intestinal miRNAs. Notably, in infected cells, NSP4 promotes the formation of double-membrane vesicles, which serve as the scaffolds for replication-transcriptional complexes and protect viral RNA from intracellular destruction. In 3 years of selection, the loss of many miRNA-binding sites in general and those within the NSP4 in particular has shaped the SARS-CoV-2 genomes. With that, the descendants of the BA.2 variants were promoted as dominant strains, which define current momentum of the pandemics.

Virome analysis provides new insights into the association between viruses and Parkinson's disease

Parkinson's disease (PD) is a kind of neurodegenerative disease that causes a huge burden to society. Previous studies have suggested the association between PD and multiple viruses. However, there is still a lack of a virome study about PD. This study systematically identified viruses from the public RNA-sequencing data of more than 700 samples from both PD patients and the control group (most were healthy people). Only nine viruses such as human betaherpesvirus 5 and Merkel cell polyomavirus have been detected in several human brain tissues of the central nervous system, the appendix, and blood of PD patients, and all of these viruses were also detected in the control group. Most viruses were observed to have low abundance in no more than three tissues. No statistically significant differences were observed between the virus abundance in the PD patients and the control group for all viruses. The positive rates of most viruses in PD patients were higher or similar to that in the control group, although those were less than 5% for most viruses. Overall, this is the first study to systematically investigate the virome in PD patients, and provides new insights into the association between viruses and PD.

Short-term forecasts of Monkeypox cases in multiple countries: keep calm and don't panic

We propose a parametric regression model for incidence indicators based on the use of the Richards' curve (a generalized logistic function) to analyse the current Monkeypox epidemic data for the most 10 affected countries worldwide. At present, results show that the outbreak is under control in most countries.

Evaluation of diagnostic test procedures for SARS-CoV-2 using latent class models

Polymerase chain reaction (PCR) and antigen tests have been used extensively for screening during the severe acute respiratory syndrome coronavirus 2 pandemics. However, the real‐world sensitivity and specificity of the two testing procedures in the field have not yet been estimated without assuming that the PCR constitutes a gold standard test. We use latent class models to estimate the in situ performance of both tests using data from the Danish national registries. We find that the specificity of both tests is very high (>99.7%), while the sensitivities are 95.7% (95% confidence interval [CI]: 92.8%–98.4%) and 53.8% (95% CI: 49.8%–57.9%) for the PCR and antigen tests, respectively. These findings have implications for the use of confirmatory PCR tests following a positive antigen test result: we estimate that serial testing is counterproductive at higher prevalence levels.

Case Report: Genomic Characteristics of the First Known Case of SARS-CoV-2 Imported From Spain to Sichuan, China

The pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been basically under control in China since March 2020, but the import of domestic SARS-CoV-2 has begun to increase. This study reported the first case of asymptomatic SARS-CoV-2 infection imported from Spain into Sichuan Province, China, on March 11, 2020. The infected male had a body temperature of 37.5°C, normal blood oxygen saturation levels, and a computed tomography (CT) examination showed that his lungs had no shadows. However, a throat swab from the subject tested positive for SARS-CoV-2 using qPCR assay. In this study, we conducted transcriptome sequencing on respiratory throat swabs from the subject and found that the dominant SARS-CoV-2 sequence (Gene Bank ID: MW301121) was a spike protein D614G mutant strain, which is currently popular throughout world. We downloaded and analyzed SARS-CoV-2 sequences collected from cases in China and Spain for comparison and tracing purposes. After March 11, 2020, the Chinese domestic clade was naturally divided into the imported SARS-CoV-2 D614G mutant strain and evolutionarily-related similar sequences and that of sequences collected in the original Wuhan area. The sequence reported in this study was located on a small branch, far from the evolution of Wuhan sequences. As expected, the identified sequence was closely related to the evolution of the SARS-CoV-2 D614G mutant strain circulating in Spain.

Understanding epidemic data and statistics: A case study of COVID-19

The 2019 novel-coronavirus (COVID-19) has affected 181 countries with approximately 1197405 confirmed cases (by 5th April). Understanding the transmission dynamics of the infection in each country which got affected on a daily basis and evaluating the effectiveness of control policies are critical for our further actions. To date, the statistics of COVID-19 reported cases show that more than 80% of infected are mild cases of disease, around 14% of infected have severe complications, and about 5% are categorized as critical disease victims. Today's report (5th April 2020; daily updates in the prepared website) shows that the confirmed cases of COVID-19 in the United States, Spain, Italy, and Germany are 308850, 126168, 124632, and 96092, respectively. Calculating the total case fatality rate (CFR) of Italy (4th April 2020), about 13.3% of confirmed cases have passed away. Compared with South Korea's rate of 1.8% (seven times lower than Italy) and China's 4% (69% lower than Italy), the CFR of Italy is too high. Some effective policies that yielded significant changes in the trend of cases were the lockdown policy in China, Italy, and Spain (the effect observed after some days), the shutdown of all nonessential companies in Hubei (the effect observed after 5 days), combined policy in South Korea, and reducing working hours in Iran.

Validity of administrative database detection of previously resolved hepatitis B virus in Japan

The risk of hepatitis B virus (HBV) reactivation has increased owing to advances in the immunosuppressive therapy field. However, the HBV reactivation incidence among patients with previously resolved HBV (prHBV) infection during immunosuppressive therapy for rheumatoid arthritis (RA) remains unclear. The objective of this work is to describe the validity of detecting prHBV infection from administrative data through comparisons with chart abstraction and determine the incidence of HBV reactivation during immunosuppressive therapy for RA in Japan. In this retrospective cohort study, data on selected patients were extracted from administrative claims data. To identify patients with prHBV infection and de novo hepatitis, and HBsAg carriers, we conducted chart abstraction. The incidence rate of de novo hepatitis was 1.23 of 100 person‐years. The positive predictive value (PPV) and its 95% confidence interval (CI) of administrative data for the identification of suspected prHBV infections was 85.8% (95% CI: 81.7%‐89.3%). This study evaluated the PPV of the algorithm of HBV‐DNA testing with immunosuppressive therapy performed four times or more per year for the detection of prHBV infection from administrative data. Additionally, we determined the incidence rate of HBV reactivation among preHBV infections during immunosuppressive therapy for RA to be 1.23 of 100 person‐years.

Genotype I hepatitis A virus introduction in Italy: Bayesian phylogenetic analysis to date different epidemics

Despite a significant decrease in acute hepatitis A in the last 2 decades in Italy, outbreaks were observed occurring mostly in southern Italy. In this study, Bayesian phylogenetic analysis was used to analyze the origin of these epidemics. With this aim, 5 different data sets of hepatitis A virus sequences were built to perform genotyping by the neighbor-joining method to estimate the evolutionary rates by using a Bayesian Markov chain Monte Carlo approach and to investigate the demographic history by independent Markov chain Monte Carlo runs enforcing both a strict and relaxed clock. The estimated mean value of the evolutionary rate, representing Ia and Ib strains, was 1.21 × 10^-3 and 2.0 × 10^-3 substitutions/site/year, respectively. The Bayesian maximum clade credibility tree of hepatitis A virus (HAV) Ia and Ib strains showed that Italian sequences mostly formed separate clusters. The root of the time for the most recent common ancestor (tMRCA) for HAV Ia and Ib strains dated back to 1981 and to 1988, respectively, showing in both cases different epidemic entrances. Phylodynamic analysis showed that genotype Ia increased in 1997, when the Apulia epidemic started, then suffered a bottleneck, probably consequent to vaccination and to the herd immunity, followed by a new increase in virus population in the years 2013-2014 consequent to the epidemic caused by the ingestion of mixed frozen berries. A similar trend without an evident bottleneck was observed also in the case of genotype Ib. In conclusion, the Bayesian phylogenetic analysis represents a good tool to measure the effectiveness of the public health plans used for HAV control.

Evaluation of nasal and nasopharyngeal swab collection for the detection of Epstein-Barr virus in nasopharyngeal carcinoma

Epstein-Barr virus detection using nasopharyngeal swabs has been suggested as a potential screening test that could improve the specificity of current EBV-based serological assays. However, application requires insertion of the swab deep into the nasopharynx, a procedure not amenable to non-clinic screening. We reasoned that swabbing the more easily accessible nasal cavity might provide an appealing alternative for NPC detection. Patients > 18 years of age diagnosed with histologically confirmed NPC were recruited from the Otolaryngology Department at the National Taiwan University Hospital. ENT clinicians collected both nasal and nasopharyngeal swabs. EBV DNA and cellular beta-globulin DNA were quantified using quantitative PCR targeting a highly-conserved region of the BKRF1 gene. EBV DNA was detectable (non-zero) in all 34 nasopharyngeal swabs and above the positivity threshold of 1666 EBV copies in 30 (88.2%) patients. EBV DNA was detectable in 50% of 34 nasal swabs and above the positivity threshold in four (11.8%) patients. Average EBV DNA levels were >3-fold higher (P < 0.001) in nasopharyngeal compared to nasal swabs. Among the 17 NPC patients with detectable EBV DNA in both swab types, we observed correlation (P < 0.01) between EBV DNA measurements. Our data represent the first evaluation of EBV DNA collected from nasal swabs. Given current EBV DNA amplification techniques, nasopharyngeal swabs remain more sensitive than nasal swabs for NPC detection.

Human papillomaviruses genotyping in plantar warts

Plantar warts are caused by human papillomaviruses (HPVs) and have been associated with several HPV genotypes. However, there are few studies focused exclusively on plantar warts. In this work, we aim to identify the HPV genotypes of plantar warts and explore their relation to demographic and clinical characteristics of patients. A total of 72 patients diagnosed with plantar warts were recruited at the Laser unit at Podiatric Hospital, University of Barcelona, Spain. Inner hyperkeratosis laminar sections of warts were collected and DNA of samples were extracted. Amplification of a conserved region of the HPV L1 gene was performed with the SK-Polymerase chain reaction method. DNA amplicons were sequenced and HPV types identified. The most prevalent genotypes detected among the 105 analyzed plantar warts were HPV-57 (37.1%), HPV-27 (23.8%), HPV-1a (20.9%), HPV-2 (15.2%), and HPV-65 (2.8%). The majority of patients (78%) presented one single plantar wart, whereas multiple warts were detected in 22.2% of patients. One patient with multiple warts presented HPV types from two different genera, suggesting the spread of warts by self-inoculation as well as by de novo infection. No significant differences between the number of warts in toes, midfoot and heel were found. The most prevalent HPV types detected in all areas belonged to the alpha genus. This work provides new insight on plantar warts and their associated HPV genotypes, and evidences the usefulness and reliability of both the sample collection procedure and the PCR method used for HPV detection and typing. J. Med. Virol. 89:902-907, 2017. © 2016 Wiley Periodicals, Inc.