Search

Obtaining a complete good-quality sequence and annotation for the long double-stranded DNA genome of the African swine fever virus (ASFV) from next-generation sequencing (NGS) technology has proven difficult, despite the increasing availability of reference genome sequences and the increasing affordability of NGS. A gap analysis conducted by the global African swine fever research alliance (GARA) partners identified that a standardized, automatic pipeline for NGS analysis was urgently needed, particularly for new outbreak strains. Whilst there are several diagnostic and research labs worldwide that collect isolates of the ASFV from outbreaks, many do not have the capability to analyze, annotate, and format NGS data from outbreaks for submission to NCBI, and some publicly available ASFV genomes have missing or incorrect annotations. We developed an automated, standardized pipeline for the analysis of NGS reads that directly provides users with assemblies and annotations formatted for their submission to NCBI. This pipeline is freely available on GitHub and has been tested through the GARA partners by examining two previously sequenced ASFV genomes; this study also aimed to assess the accuracy and limitations of two strategies present within the pipeline: reference-based (Illumina reads) and de novo assembly (Illumina and Nanopore reads) strategies.

African swine fever (ASF) is a devastating disease of domestic pigs that has spread across the globe since its introduction into Georgia in 2007. The etiological agent is a large double-stranded DNA virus with a genome of 170 to 180 kb in length depending on the isolate. Much of the differences in genome length between isolates are due to variations in the copy number of five different multigene families that are encoded in repetitive regions that are towards the termini of the covalently closed ends of the genome. Molecular epidemiology of African swine fever virus (ASFV) is primarily based on Sanger sequencing of a few conserved and variable regions, but due to the stability of the dsDNA genome changes in the variable regions occur relatively slowly. Observations in Europe and Asia have shown that changes in other genetic loci can occur and that this could be useful in molecular tracking. ASFV has been circulating in Western Africa for at least forty years. It is therefore reasonable to assume that changes may have accumulated in regions of the genome other than the standard targets over the years. At present only one full genome sequence is available for an isolate from Western Africa, that of a highly virulent isolate collected from Benin during an outbreak in 1997. In Cameroon, ASFV was first reported in 1981 and outbreaks have been reported to the present day and is considered endemic. Here we report three full genome sequences from Cameroon isolates of 1982, 1994 and 2018 outbreaks and identify novel single nucleotide polymorphisms and insertion-deletions that may prove useful for molecular epidemiology studies in Western Africa and beyond.

Lumpy skin disease virus (LSDV) is a member of the capripoxvirus (CPPV) genus of the Poxviridae family. LSDV is a rapidly emerging, high-consequence pathogen of cattle, recently spreading from Africa and the Middle East into Europe and Asia. We have sequenced the whole genome of historical LSDV isolates from the Pirbright Institute virus archive, and field isolates from recent disease outbreaks in Sri Lanka, Mongolia, Nigeria and Ethiopia. These genome sequences were compared to published genomes and classified into different subgroups. Two subgroups contained vaccine or vaccine-like samples (“Neethling-like” clade 1.1 and “Kenya-like” subgroup, clade 1.2.2). One subgroup was associated with outbreaks of LSD in the Middle East/Europe (clade 1.2.1) and a previously unreported subgroup originated from cases of LSD in west and central Africa (clade 1.2.3). Isolates were also identified that contained a mix of genes from both wildtype and vaccine samples (vaccine-like recombinants, grouped in clade 2). Whole genome sequencing and analysis of LSDV strains isolated from different regions of Africa, Europe and Asia have provided new knowledge of the drivers of LSDV emergence, and will inform future disease control strategies.

Abstract Background Non-targeted whole genome sequencing is a powerful tool to comprehensively identify constituents of microbial communities in a sample. There is no need to direct the analysis to any identification before sequencing which can decrease the introduction of bias and false negatives results. It also allows the assessment of genetic aberrations in the genome (e.g., single nucleotide variants, deletions, insertions and copy number variants) including in noncoding protein regions. Methods The performance of four different random priming amplification methods to recover RNA viral genetic material of SARS-CoV-2 were compared in this study. In method 1 (H-P) the reverse transcriptase (RT) step was performed with random hexamers whereas in methods 2–4 RT incorporating an octamer primer with a known tag. In methods 1 and 2 (K-P) sequencing was applied on material derived from the RT-PCR step, whereas in methods 3 (SISPA) and 4 (S-P) an additional amplification was incorporated before sequencing. Results The SISPA method was the most effective and efficient method for non-targeted/random priming whole genome sequencing of SARS-CoV-2 that we tested. The SISPA method described in this study allowed for whole genome assembly of SARS-CoV-2 and influenza A(H1N1)pdm09 in mixed samples. We determined the limit of detection and characterization of SARS-CoV-2 virus which was 103 pfu/ml (Ct, 22.4) for whole genome assembly and 101 pfu/ml (Ct, 30) for metagenomics detection. Conclusions The SISPA method is predominantly useful for obtaining genome sequences from RNA viruses or investigating complex clinical samples as no prior sequence information is needed. It might be applied to monitor genomic virus changes, virus evolution and can be used for fast metagenomics detection or to assess the general picture of different pathogens within the sample.

Abstract Background Hepatocellular carcinoma is the second most deadly cancer with late presentation and limited treatment options, highlighting an urgent need to better understand HCC to facilitate the identification of early-stage biomarkers and uncover therapeutic targets for the development of novel therapies for HCC. Methods Deep transcriptome sequencing of tumor and paired non-tumor liver tissues was performed to comprehensively evaluate the profiles of both the host and HBV transcripts in HCC patients. Differential gene expression patterns and the dys-regulated genes associated with clinical outcomes were analyzed. Somatic mutations were identified from the sequencing data and the deleterious mutations were predicted. Lastly, human-HBV chimeric transcripts were identified, and their distribution, potential function and expression association were analyzed. Results Expression profiling identified the significantly upregulated TP73 as a nodal molecule modulating expression of apoptotic genes. Approximately 2.5% of dysregulated genes significantly correlated with HCC clinical characteristics. Of the 110 identified genes, those involved in post-translational modification, cell division and/or transcriptional regulation were upregulated, while those involved in redox reactions were downregulated in tumors of patients with poor prognosis. Mutation signature analysis identified that somatic mutations in HCC tumors were mainly non-synonymous, frequently affecting genes in the micro-environment and cancer pathways. Recurrent mutations occur mainly in ribosomal genes. The most frequently mutated genes were generally associated with a poorer clinical prognosis. Lastly, transcriptome sequencing suggest that HBV replication in the tumors of HCC patients is rare. HBV-human fusion transcripts are a common observation, with favored HBV and host insertion sites being the HBx C-terminus and gene introns (in tumors) and introns/intergenic-regions (in non-tumors), respectively. HBV-fused genes in tumors were mainly involved in RNA binding while those in non-tumors tissues varied widely. These observations suggest that while HBV may integrate randomly during chronic infection, selective expression of functional chimeric transcripts may occur during tumorigenesis. Conclusions Transcriptome sequencing of HCC patients reveals key cancer molecules and clinically relevant pathways deregulated/mutated in HCC patients and suggests that while HBV may integrate randomly during chronic infection, selective expression of functional chimeric transcripts likely occur during the process of tumorigenesis.

Foot-and-mouth disease (FMD) is endemic in large parts of sub-Saharan Africa, Asia and South America, where outbreaks in cloven-hooved livestock threaten food security and have severe economic impacts. Vaccination in endemic regions remains the most effective control strategy. Current FMD vaccines are produced from chemically inactivated foot-and-mouth disease virus (FMDV) grown in suspension cultures of baby hamster kidney 21 cells (BHK-21). Strain diversity means vaccines produced from one subtype may not fully protect against circulating disparate subtypes, necessitating the development of new vaccine strains that “antigenically match”. However, some viruses have proven difficult to adapt to cell culture, slowing the manufacturing process, reducing vaccine yield and limiting the availability of effective vaccines, as well as potentiating the selection of undesired antigenic changes. To circumvent the need to cell culture adapt FMDV, we have used a systematic approach to develop recombinant suspension BHK-21 that stably express the key FMDV receptor integrin αvβ6. We show that αvβ6 expression is retained at consistently high levels as a mixed cell population and as a clonal cell line. Following exposure to field strains of FMDV, these recombinant BHK-21 facilitated higher virus yields compared to both parental and control BHK-21, whilst demonstrating comparable growth kinetics. The presented data supports the application of these recombinant αvβ6-expressing BHK-21 in future FMD vaccine production.

Peste des petits ruminants virus (PPRV) causes a highly devastating disease of sheep and goats that threatens food security, small ruminant production and susceptible endangered wild ruminants. With policy directed towards achieving global PPR eradication, the establishment of cost-effective genomic surveillance tools is critical where PPR is endemic. Genomic data can provide sufficient in-depth information to identify the pockets of endemicity responsible for PPRV persistence and viral evolution, and direct an appropriate vaccination response. Yet, access to the required sequencing technology is low in resource-limited settings and is compounded by the difficulty of transporting clinical samples from wildlife across international borders due to the Convention on International Trade in Endangered Species (CITES) of Wild Fauna and Flora, and Nagoya Protocol regulations. Oxford nanopore MinION sequencing technology has recently demonstrated an extraordinary performance in the sequencing of PPRV due to its rapidity, utility in endemic countries and comparatively low cost per sample when compared to other whole-genome (WGS) sequencing platforms. In the present study, Oxford nanopore MinION sequencing was utilised to generate complete genomes of PPRV isolates collected from infected goats in Ngorongoro and Momba districts in the northern and southern highlands of Tanzania during 2016 and 2018, respectively. The tiling multiplex polymerase chain reaction (PCR) was carried out with twenty-five pairs of long-read primers. The resulting PCR amplicons were used for nanopore library preparation and sequencing. The analysis of output data was complete genomes of PPRV, produced within four hours of sequencing (accession numbers: MW960272 and MZ322753). Phylogenetic analysis of the complete genomes revealed a high nucleotide identity, between 96.19 and 99.24% with lineage III PPRV currently circulating in East Africa, indicating a common origin. The Oxford nanopore MinION sequencer can be deployed to overcome diagnostic and surveillance challenges in the PPR Global Control and Eradication program. However, the coverage depth was uneven across the genome and amplicon dropout was observed mainly in the GC-rich region between the matrix (M) and fusion (F) genes of PPRV. Thus, larger field studies are needed to allow the collection of sufficient data to assess the robustness of nanopore sequencing technology.

Many viruses that cause serious diseases in humans and animals, including the betacoronaviruses (beta-CoVs), such as SARS-CoV, MERS-CoV, and the recently identified SARS-CoV-2, have natural reservoirs in bats. Because these viruses rely entirely on the host cellular machinery for survival, their evolution is likely to be guided by the link between the codon usage of the virus and that of its host. As a result, specific cellular microenvironments of the diverse hosts and/or host tissues imprint peculiar molecular signatures in virus genomes. Our study is aimed at deciphering some of these signatures. Using a variety of genetic methods we demonstrated that trends in codon usage across chiroptera-hosted CoVs are collaboratively driven by geographically different host-species and temporal-spatial distribution. We not only found that chiroptera-hosted CoVs are the ancestors of SARS-CoV-2, but we also revealed that SARS-CoV-2 has the codon usage characteristics similar to those seen in CoVs infecting the Rhinolophus sp. Surprisingly, the envelope gene of beta-CoVs infecting Rhinolophus sp., including SARS-CoV-2, had extremely high CpG levels, which appears to be an evolutionarily conserved trait. The dissection of the furin cleavage site of various CoVs infecting hosts revealed host-specific preferences for arginine codons; however, arginine is encoded by a wider variety of synonymous codons in the murine CoV (MHV-A59) furin cleavage site. Our findings also highlight the latent diversity of CoVs in mammals that has yet to be fully explored.

Abstract Background There are a large number of biological databases publicly available for scientists in the web. Also, there are many private databases generated in the course of research projects. These databases are in a wide variety of formats. Web standards have evolved in the recent times and semantic web technologies are now available to interconnect diverse and heterogeneous sources of data. Therefore, integration and querying of biological databases can be facilitated by techniques used in semantic web. Heterogeneous databases can be converted into Resource Description Format (RDF) and queried using SPARQL language. Searching for exact queries in these databases is trivial. However, exploratory searches need customized solutions, especially when multiple databases are involved. This process is cumbersome and time consuming for those without a sufficient background in computer science. In this context, a search engine facilitating exploratory searches of databases would be of great help to the scientific community. Results We present BioCarian, an efficient and user-friendly search engine for performing exploratory searches on biological databases. The search engine is an interface for SPARQL queries over RDF databases. We note that many of the databases can be converted to tabular form. We first convert the tabular databases to RDF. The search engine provides a graphical interface based on facets to explore the converted databases. The facet interface is more advanced than conventional facets. It allows complex queries to be constructed, and have additional features like ranking of facet values based on several criteria, visually indicating the relevance of a facet value and presenting the most important facet values when a large number of choices are available. For the advanced users, SPARQL queries can be run directly on the databases. Using this feature, users will be able to incorporate federated searches of SPARQL endpoints. We used the search engine to do an exploratory search on previously published viral integration data and were able to deduce the main conclusions of the original publication. BioCarian is accessible via http://www.biocarian.com . Conclusions We have developed a search engine to explore RDF databases that can be used by both novice and advanced users.

ALK is an established causative oncogenic driver in neuroblastoma, and is likely to emerge as a routine biomarker in neuroblastoma diagnostics. At present, the optimal strategy for clinical diagnostic evaluation of ALK protein, genomic and hotspot mutation status is not well-studied. We evaluated ALK immunohistochemical (IHC) protein expression using three different antibodies (ALK1, 5A4 and D5F3 clones), ALK genomic status using single-color chromogenic in situ hybridization (CISH), and ALK hotspot mutation status using conventional Sanger sequencing and a next-generation sequencing platform (Ion Torrent Personal Genome Machine (IT-PGM)), in archival formalin-fixed, paraffin-embedded neuroblastoma samples. We found a significant difference in IHC results using the three different antibodies, with the highest percentage of positive cases seen on D5F3 immunohistochemistry. Correlation with ALK genomic and hotspot mutational status revealed that the majority of D5F3 ALK-positive cases did not possess either ALK genomic amplification or hotspot mutations. Comparison of sequencing platforms showed a perfect correlation between conventional Sanger and IT-PGM sequencing. Our findings suggest that D5F3 immunohistochemistry, single-color CISH and IT-PGM sequencing are suitable assays for evaluation of ALK status in future neuroblastoma clinical trials.