Candidate historical events for the emergence of Human Coronavirus OC43: A critical reassessment of the molecular evidence.

The "Russian Influenza"-coronavirus theory (RICT) proposes that the pandemic of 1889-1892, conventionally regarded as an influenza pandemic, was caused by the emergence of human coronavirus OC43 (HCoV-OC43) as a zoonosis of bovine coronavirus (BCoV). RICT is based on a Bayesian phylogenetic calculation of the date of the most recent common ancestor (MRCA) of HCoV-OC43 and BCoV. The theory also draws on comparison of both symptoms and some epidemiological parameters of the best studied coronavirus pandemic, i.e. COVID-19, with those reported in 1889-1892. The case is completed with circumstantial evidence involving a panzoonotic among cattle in the decade prior to the "Russian Influenza", with characteristics suggesting it may have been caused by BCoV. In this paper, we review the Bayesian phylogenetic evidence for RICT, replicating previous studies and adding our own, in each case critically reviewing the suitability of the datasets used and the parameters applied. We conclude that the most probable date for the MRCA of HCoV-OC43 and BCoV is 1898-1902. This is a decade too late for compatibility with RICT but happens to coincide with another serious outbreak of respiratory illness, reported in both the USA and the UK, in the winter of 1899-1900.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World.

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) procaviensis Isolate 253, Strain LV425.

Leishmania (Mundinia) procaviensis is a parasitic kinetoplastid that was first isolated from a rock hyrax in Namibia in 1975. We present the complete genome sequence of Leishmania (Mundinia) procaviensis isolate 253, strain LV425, sequenced using combined short- and long-read technologies. This genome will contribute to our understanding of hyraxes as a Leishmania reservoir.

Chromosome-Scale Assembly of the Complete Genome Sequence of Porcisia hertigi, Isolate C119, Strain LV43.

Porcisia hertigi is a parasitic kinetoplastid first isolated from porcupines (Coendou rothschildi) in central Panama in 1965. We present the complete genome sequence of P. hertigi, isolate C119, strain LV43, sequenced using combined short- and long-read technologies. This complete genome sequence will contribute to our knowledge of the parasitic genus Porcisia.

ICTV Virus Taxonomy Profile: Herpesviridae 2021.

Members of the family Herpesviridae have enveloped, spherical virions with characteristic complex structures consisting of symmetrical and non-symmetrical components. The linear, double-stranded DNA genomes of 125-241 kbp contain 70-170 genes, of which 43 have been inherited from an ancestral herpesvirus. In general, herpesviruses have coevolved with and are highly adapted to their hosts, which comprise many mammalian, avian and reptilian species. Following primary infection, they are able to establish lifelong latent infection, during which there is limited viral gene expression. Severe disease is usually observed only in the foetus, the very young, the immunocompromised or following infection of an alternative host. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Herpesviridae, which is available at ictv.global/report/herpesviridae.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) sp. Ghana, Isolate GH5, Strain LV757.

Leishmania (Mundinia) sp. Ghana is a kinetoplastid parasite isolated in 2015 in Ghana. We report the complete genome sequence of L. (M.) sp. Ghana, sequenced using combined short-read and long-read technologies. This will facilitate greater understanding of this novel pathogen and its relationships within the subgenus Mundinia.

Chromosome-scale genome sequencing, assembly and annotation of six genomes from subfamily Leishmaniinae.

We provide the raw and processed data produced during the genome sequencing of isolates from six species of parasites from the sub-family Leishmaniinae: Leishmania martiniquensis (Thailand), Leishmania orientalis (Thailand), Leishmania enriettii (Brazil), Leishmania sp. Ghana, Leishmania sp. Namibia and Porcisia hertigi (Panama). De novo assembly was performed using Nanopore long reads to construct chromosome backbone scaffolds. We then corrected erroneous base calling by mapping short Illumina paired-end reads onto the initial assembly. Data has been deposited at NCBI as follows: raw sequencing output in the Sequence Read Archive, finished genomes in GenBank, and ancillary data in BioSample and BioProject. Derived data such as quality scoring, SAM files, genome annotations and repeat sequence lists have been deposited in Lancaster University's electronic data archive with DOIs provided for each item. Our coding workflow has been deposited in GitHub and Zenodo repositories. This data constitutes a resource for the comparative genomics of parasites and for further applications in general and clinical parasitology.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) enriettii, Isolate CUR178, Strain LV763.

Leishmania (Mundinia) enriettii is a parasitic kinetoplastid first isolated from a guinea pig in Brazil in 1946. We present the complete genome sequence of L. (M.) enriettii, isolate CUR178, strain LV763, sequenced using combined short-read and long-read technologies. This will facilitate a greater understanding of the genome diversity within L. (M.) enriettii.

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) orientalis, Isolate LSCM4, Strain LV768.

Leishmania (Mundinia) orientalis is a kinetoplastid parasite first isolated in 2014 in Thailand. We report the complete genome sequence of L. (M.) orientalis, sequenced using combined short-read and long-read technologies. This will facilitate greater understanding of this novel pathogen and its relationship to other members of the subgenus Mundinia.

LGAAP: Leishmaniinae Genome Assembly and Annotation Pipeline.

We present the LGAAP computational pipeline, which was successfully used to assemble six genomes of the parasite subfamily Leishmaniinae to chromosome-scale completeness from a combination of long- and short-read sequencing data. LGAAP is open source, and we suggest that it may easily be ported for assembly of any genome of comparable size (∼35 Mb).

Chromosome-Scale Assembly of the Complete Genome Sequence of Leishmania (Mundinia) martiniquensis, Isolate LSCM1, Strain LV760.

Leishmania (Mundinia) martiniquensis is a kinetoplastid parasite that was first isolated in 1995 on Martinique. We report the first complete genome for Leishmania martiniquensis from Asia, isolate LSCM1, strain LV760, which was sequenced using combined short-read and long-read technologies. This will facilitate greater understanding of the evolution of the geographically dispersed subgenus Mundinia.

Mitochondrial DNA Profiling Reveals Two Lineages of Sun Bears in East and West Malaysia.

Sun bear populations are fragmented and at risk from habitat loss and exploitation for body parts. These threats are made worse by significant gaps in knowledge of sun bear population genetic diversity, population connectivity, and taxonomically significant management units. Using a complete sun bear mitochondrial genome, we developed a set of mitochondrial markers to assess haplotype variation and the evolutionary history of sun bears from Peninsular (West) Malaysia and Sabah (East Malaysia). Genetic samples from 28 sun bears from Peninsular Malaysia, 36 from Sabah, and 18 from Thailand were amplified with primers targeting a 1800 bp region of the mitochondrial genome including the complete mitochondrial control region and adjacent genes. Sequences were analyzed using phylogenetic methods. We identified 51 mitochondrial haplotypes among 82 sun bears. Phylogenetic and network analyses provided strong support for a deep split between Malaysian sun bears and sun bears in East Thailand and Yunnan province in China. The Malaysian lineage was further subdivided into two clades: Peninsular Malaysian and Malaysian Borneo (Sabah). Sun bears from Thailand occurred in both Sabah and Peninsular Malaysian clades. Our study supports recent findings that sun bears from Sundaland form a distinct clade from those in China and Indochina with Thailand possessing lineages from the three clades. Importantly, we demonstrate a more recent and clear genetic delineation between sun bears from the Malay Peninsula and Sabah indicating historical barriers to gene flow within the Sundaic region.

Insecticide-impregnated netting: A surface treatment for killing Lutzomyia longipalpis (Diptera: Psychodidae), the vector of Leishmania infantum.

The sand fly Lutzomyia longipalpis is the main vector of Leishmania infantum in Brazil. Synthetic male-produced sex/aggregation pheromone co-located with micro-encapsulated λ-cyhalothrin in chicken sheds can significantly reduce canine infection and sand fly densities in a lure-and-kill strategy. In this study, we determined if insecticide-impregnated netting (IN) could replace insecticide residual spraying (IRS). We compared numbers of Lu. longipalpis attracted and killed in experimental and real chicken sheds baited with pheromone and treated with a 1 m2 area of either insecticide spray or netting. First, we compared both treatments in experimental sheds to control mortality established from light trap captures. We then compared the long-term killing effect of insecticide spray and netting, without renewal, in experimental sheds over a period of 16 weeks. Finally, a longitudinal intervention study in real chicken sheds compared the numbers and proportions of Lu. longipalpis collected and killed before and after application of both treatments. In Experiment 1, a higher proportion of males and females captured in IRS- and IN-treated sheds were dead at 24 h compared to controls (P < 0.05). No difference was found in the proportion of females killed in sheds treated with IN or IRS (P = 0.15). A slightly higher proportion of males were killed by IRS (100%) compared to IN (98.6%; P < 0.05). In Experiment 2, IN- and IRS-treated traps were equally effective at killing females (P = 0.21) and males (P = 0.08). However, IRS killed a significantly higher proportion of females and males after 8 (P < 0.05) and 16 (P < 0.05) weeks. In Experiment 3, there was no significant difference between treatments in the proportion of females killed before (P = 0.88) or after (P = 0.29) or males killed before (P = 0.76) or after (P = 0.73) intervention. Overall, initially the IN was as effective as IRS at killing female and male Lu. longipalpis in both experimental and real chicken sheds. However, the relative lethal effect of the IN deteriorated over time when stored under prevailing environmental conditions.

Reflections on integrating bioinformatics into the undergraduate curriculum: The Lancaster experience.

Bioinformatics is an essential discipline for biologists. It also has a reputation of being difficult for those without a strong quantitative and computer science background. At Lancaster University, we have developed modules for the integration of bioinformatics skills training into our undergraduate biology degree portfolio. This article describes those modules, situating them in the context of the accumulated quarter century of literature on bioinformatics education. The constant evolution of bioinformatics as a discipline is emphasized, drawing attention to the continual necessity to revise and upgrade those skills being taught, even at undergraduate level. Our overarching aim is to equip students both with a portfolio of skills in the currently most essential bioinformatics tools and with the confidence to continue their own bioinformatics skills development at postgraduate or professional level.

Visualization of protein sequence space with force-directed graphs, and their application to the choice of target-template pairs for homology modelling.

The protein sequence-structure gap results from the contrast between rapid, low-cost deep sequencing, and slow, expensive experimental structure determination techniques. Comparative homology modelling may have the potential to close this gap by predicting protein structure in target sequences using existing experimentally solved structures as templates. This paper presents the first use of force-directed graphs for the visualization of sequence space in two dimensions, and applies them to the choice of suitable RNA-dependent RNA polymerase (RdRP) target-template pairs within human-infective RNA virus genera. Measures of centrality in protein sequence space for each genus were also derived and used to identify centroid nearest-neighbour sequences (CNNs) potentially useful for production of homology models most representative of their genera. Homology modelling was then carried out for target-template pairs in different species, different genera and different families, and model quality assessed using several metrics. Reconstructed ancestral RdRP sequences for individual genera were also used as templates for the production of ancestral RdRP homology models. High quality ancestral RdRP models were consistently produced, as were good quality models for target-template pairs in the same genus. Homology modelling between genera in the same family produced mixed results and inter-family modelling was unreliable. We present a protocol for the production of optimal RdRP homology models for use in further experiments, e.g. docking to discover novel anti-viral compounds. (219 words).

Delineating the Plausible Molecular Vaccine Candidates and Drug Targets of Multidrug-Resistant Acinetobacter baumannii.

Nosocomial infections have become alarming with the increase of multidrug-resistant bacterial strains of Acinetobacter baumannii. Being the causative agent in ~80% of the cases, these pathogenic gram-negative species could be deadly for hospitalized patients, especially in intensive care units utilizing ventilators, urinary catheters, and nasogastric tubes. Primarily infecting an immuno-compromised system, they are resistant to most antibiotics and are the root cause of various types of opportunistic infections including but not limited to septicemia, endocarditis, meningitis, pneumonia, skin, and wound sepsis and even urinary tract infections. Conventional experimental methods including typing, computational methods encompassing comparative genomics, and combined methods of reverse vaccinology and proteomics had been proposed to differentiate and develop vaccines and/or drugs for several outbreak strains. However, identifying proteins suitable enough to be posed as drug targets and/or molecular vaccines against the multidrug-resistant pathogenic bacterial strains has probably remained an open issue to address. In these cases of novel protein identification, the targets either are uncharacterized or have been unable to confer the most coveted protection either in the form of molecular vaccine candidates or as drug targets. Here, we report a strategic approach with the 3,766 proteins from the whole genome of A. baumannii ATCC19606 (AB) to rationally identify plausible candidates and propose them as future molecular vaccine candidates and/or drug targets. Essentially, we started with mapping the vaccine candidates (VaC) and virulence factors (ViF) of A. baumannii strain AYE onto strain ATCC19606 to identify them in the latter. We move on to build small networks of VaC and ViF to conceptualize their position in the network space of the whole genomic protein interactome (GPIN) and rationalize their candidature for drugs and/or molecular vaccines. To this end, we propose new sets of known proteins unearthed from interactome built using key factors, KeF, potent enough to compete with VaC and ViF. Our method is the first of its kind to propose, albeit theoretically, a rational approach to identify crucial proteins and pose them for candidates of vaccines and/or drugs effective enough to combat the deadly pathogenic threats of A. baumannii.

A decade of sustained selection pressure on two surface sites of the VP1 protein of Enterovirus A71 suggests that immune evasion may be an indirect driver for virulence.

Enterovirus A71 (EV-A71) is an emerging pathogen in the Enterovirus A species group. EV-A71 causes hand, foot and mouth disease (HFMD), with virulent variants exhibiting polio-like acute flaccid paralysis and other central nervous system manifestations. We analysed all enterovirus A71 complete genomes with collection dates from 2008 to mid-2018. All sub-genotypes exhibit a strong molecular clock with omega (dN/dS) suggesting strong purifying selection. In sub-genotypes B5 and C4, positive selection can be detected at two surface sites on the VP1 protein, also detected in positive selection studies performed prior to 2008. Toggling of a limited repertoire of amino acids at these positively selected residues over the last decade suggests that EV-A71 may be undergoing a sustained frequency-dependent selection process for immune evasion, raising issues for vaccine development. These same sites have also been previously implicated in virus-host binding and strain-associated severity of HFMD, suggesting that immune evasion may be an indirect driver for virulence (154 words).

Sequence Variation in Multidrug-Resistant Plasmid pLUH01, Isolated from Human Nasopharyngeal Swabs.

Three variants of the multidrug-resistant plasmid pLUH01 were assembled by deep sequencing from nasopharyngeal swabs. All have a 21-bp deletion in the RS14515 hypothetical gene. Variants 1 through 3 have 2, 6, and 3 nucleotide substitutions, respectively, compared to the pLUH01 reference genome. We named the new plasmid variants pLUH01/Lancaster/2015/1 to pLUH01/Lancaster/2015/3.

Paradigm Shift in Drug Re-purposing From Phenalenone to Phenaleno-Furanone to Combat Multi-Drug Resistant Salmonella enterica Serovar Typhi.

Over recent years, typhoid fever has gained increasing attention with several cases reporting treatment failure due to multidrug resistant (MDR) strains of Salmonella enterica serovar Typhi. While new drug development strategies are being devised to combat the threat posed by these MDR pathogens, drug repurposing or repositioning has become a good alternative. The latter is considered mainly due to its capacity for saving sufficient time and effort for pre-clinical and optimization studies. Owing to the possibility of an unsuccessful repositioning, due to the mismatch in the optimization of the drug ligand for the changed biochemical properties of "old" and "new" targets, we have chosen a "targeted" approach of adopting a combined chemical moiety-based drug repurposing. Using small molecules selected from a combination of earlier approved drugs having phenalenone and furanone moieties, we have computationally delineated a step-wise approach to drug design against MDR Salmonella. We utilized our network analysis-based pre-identified, essential chaperone protein, SicA, which regulates the folding and quality of several secretory proteins including the Hsp70 chaperone, SigE. To this end, another crucial chaperone protein, Hsp70 DnaK, was also considered due to its importance for pathogen survival under the stress conditions typically encountered during antibiotic therapies. These were docked with the 19 marketed anti-typhoid drugs along with two phenalenone-furanone derivatives, 15 non-related drugs which showed 70% similarity to phenalenone and furanone derivatives and other analogous small molecules. Furthermore, molecular dynamics simulation studies were performed to check the stability of the protein-drug complexes. Our results showed the best binding interaction and stability, under the parameters of a virtual human body environment, with XR770, a phenaleno-furanone moiety based derivative. We therefore propose XR770, for repurposing for therapeutic intervention against emerging and significant drug resistance conferred by pathogenic Salmonella strains.