Metagenomic surveillance uncovers diverse and novel viral taxa in febrile patients from Nigeria

Effective infectious disease surveillance in high-risk regions is critical for clinical care and pandemic preemption; however, few clinical diagnostics are available for the wide range of potential human pathogens. Here, we conduct unbiased metagenomic sequencing of 593 samples from febrile Nigerian patients collected in three settings: i) population-level surveillance of individuals presenting with symptoms consistent with Lassa Fever (LF); ii) real-time investigations of outbreaks with suspected infectious etiologies; and iii) undiagnosed clinically challenging cases. We identify 13 distinct viruses, including the second and third documented cases of human blood-associated dicistrovirus, and a highly divergent, unclassified dicistrovirus that we name human blood-associated dicistrovirus 2. We show that pegivirus C is a common co-infection in individuals with LF and is associated with lower Lassa viral loads and favorable outcomes. We help uncover the causes of three outbreaks as yellow fever virus, monkeypox virus, and a noninfectious cause, the latter ultimately determined to be pesticide poisoning. We demonstrate that a local, Nigerian-driven metagenomics response to complex public health scenarios generates accurate, real-time differential diagnoses, yielding insights that inform policy.

Adeno-associated virus type 2 in US children with acute severe hepatitis

As of August 2022, clusters of acute severe hepatitis of unknown aetiology in children have been reported from 35 countries, including the USA1,2. Previous studies have found human adenoviruses (HAdVs) in the blood from patients in Europe and the USA3-7, although it is unclear whether this virus is causative. Here we used PCR testing, viral enrichment-based sequencing and agnostic metagenomic sequencing to analyse samples from 16 HAdV-positive cases from 1 October 2021 to 22 May 2022, in parallel with 113 controls. In blood from 14 cases, adeno-associated virus type 2 (AAV2) sequences were detected in 93% (13 of 14), compared to 4 (3.5%) of 113 controls (P < 0.001) and to 0 of 30 patients with hepatitis of defined aetiology (P < 0.001). In controls, HAdV type 41 was detected in blood from 9 (39.1%) of the 23 patients with acute gastroenteritis (without hepatitis), including 8 of 9 patients with positive stool HAdV testing, but co-infection with AAV2 was observed in only 3 (13.0%) of these 23 patients versus 93% of cases (P < 0.001). Co-infections by Epstein-Barr virus, human herpesvirus 6 and/or enterovirus A71 were also detected in 12 (85.7%) of 14 cases, with higher herpesvirus detection in cases versus controls (P < 0.001). Our findings suggest that the severity of the disease is related to co-infections involving AAV2 and one or more helper viruses.

The Prevalence of Undiagnosed Salmonella enterica Serovar Typhi in Healthy School-Aged Children in Osun State, Nigeria

Typhoid fever remains a significant public health concern due to cases of mis-/overdiagnosis. Asymptomatic carriers play a role in the transmission and persistence of typhoid fever, especially among children, where limited data exist in Nigeria and other endemic countries. We aim to elucidate the burden of typhoid fever among healthy school-aged children using the best surveillance tool(s). In a semi-urban/urban state (Osun), 120 healthy school-aged children under 15 years were enrolled. Whole blood and fecal samples were obtained from consenting children. ELISA targeting the antigen lipopolysaccharide (LPS) and anti-LPS antibodies of Salmonella Typhi, culture, polymerase chain reaction (PCR), and next-generation sequencing (NGS) were used to analyze the samples. At least one of the immunological markers was detected in 65.8% of children, with 40.8%, 37.5%, and 39% of children testing positive for IgM, IgG, and antigen, respectively. Culture, PCR, and NGS assays did not detect the presence of Salmonella Typhi in the isolates. This study demonstrates a high seroprevalence of Salmonella Typhi in these healthy children but no carriage, indicating the inability to sustain transmission. We also demonstrate that using a single technique is insufficient for typhoid fever surveillance in healthy children living in endemic areas.

Emergence and spread of two SARS-CoV-2 variants of interest in Nigeria

Identifying the dissemination patterns and impacts of a virus of economic or health importance during a pandemic is crucial, as it informs the public on policies for containment in order to reduce the spread of the virus. In this study, we integrated genomic and travel data to investigate the emergence and spread of the SARS-CoV-2 B.1.1.318 and B.1.525 (Eta) variants of interest in Nigeria and the wider Africa region. By integrating travel data and phylogeographic reconstructions, we find that these two variants that arose during the second wave in Nigeria emerged from within Africa, with the B.1.525 from Nigeria, and then spread to other parts of the world. Data from this study show how regional connectivity of Nigeria drove the spread of these variants of interest to surrounding countries and those connected by air-traffic. Our findings demonstrate the power of genomic analysis when combined with mobility and epidemiological data to identify the drivers of transmission, as bidirectional transmission within and between African nations are grossly underestimated as seen in our import risk index estimates.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance

Investment in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing in Africa over the past year has led to a major increase in the number of sequences that have been generated and used to track the pandemic on the continent, a number that now exceeds 100,000 genomes. Our results show an increase in the number of African countries that are able to sequence domestically and highlight that local sequencing enables faster turnaround times and more-regular routine surveillance. Despite limitations of low testing proportions, findings from this genomic surveillance study underscore the heterogeneous nature of the pandemic and illuminate the distinct dispersal dynamics of variants of concern-particularly Alpha, Beta, Delta, and Omicron-on the continent. Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve while the continent faces many emerging and reemerging infectious disease threats. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century.

Future-proofing and maximizing the utility of metadata: The PHA4GE SARS-CoV-2 contextual data specification package

Background

The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard.

Results

As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories.

Conclusions

Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI’s BioSample database.

HIV-1 drug resistance and genetic diversity in a cohort of people with HIV-1 in Nigeria

OBJECTIVE

This study was designed to provide information on the genetic diversity of HIV-1 and drug resistance mutations in Nigeria, as there is limited understanding of variants circulating in the country.

METHODS

We used an advanced next-generation sequencing platform, Primer ID, to: investigate the presence of high and low abundance drug resistance mutations; characterize preexisting Integrase Strand Transfer Inhibitor (INSTI) mutations in antiretroviral therapy (ART)-experienced but dolutegravir-naive individuals; detect recent HIV-1 infections and characterize subtype diversity from a cohort of people with HIV-1 (PWH).

RESULTS

HIV-1 subtype analysis revealed the predominance of CRF02_AG and subtype G in our study population. At detection sensitivity of 30% abundance, drug resistance mutations (DRMs) were identified in 3% of samples. At a sensitivity level of 10%, DRMs were identified in 27.3% of samples. We did not detect any major INSTI mutation associated with dolutegravir-resistance. Only one recent infection was detected in our study population.

CONCLUSION

Our study suggests that dolutegravir-containing antiretroviral regimens will be effective in Nigeria. Our study also further emphasizes the high genetic diversity of HIV-1 in Nigeria and that CRF02_AG and subtype G are the dominant circulating forms of HIV-1 in Nigeria. These two circulating forms of the virus are largely driving the epidemic in the country.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa

The progression of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in Africa has so far been heterogeneous, and the full impact is not yet well understood. In this study, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations predominantly from Europe, which diminished after the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1, and C.1.1. Although distorted by low sampling numbers and blind spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a source for new variants.

VGEA: an RNA viral assembly toolkit

Next generation sequencing (NGS)-based studies have vastly increased our understanding of viral diversity. Viral sequence data obtained from NGS experiments are a rich source of information, these data can be used to study their epidemiology, evolution, transmission patterns, and can also inform drug and vaccine design. Viral genomes, however, represent a great challenge to bioinformatics due to their high mutation rate and forming quasispecies in the same infected host, bringing about the need to implement advanced bioinformatics tools to assemble consensus genomes well-representative of the viral population circulating in individual patients. Many tools have been developed to preprocess sequencing reads, carry-out de novo or reference-assisted assembly of viral genomes and assess the quality of the genomes obtained. Most of these tools however exist as standalone workflows and usually require huge computational resources. Here we present (Viral Genomes Easily Analyzed), a Snakemake workflow for analyzing RNA viral genomes. VGEA enables users to map sequencing reads to the human genome to remove human contaminants, split bam files into forward and reverse reads, carry out de novo assembly of forward and reverse reads to generate contigs, pre-process reads for quality and contamination, map reads to a reference tailored to the sample using corrected contigs supplemented by the user's choice of reference sequences and evaluate/compare genome assemblies. We designed a project with the aim of creating a flexible, easy-to-use and all-in-one pipeline from existing/stand-alone bioinformatics tools for viral genome analysis that can be deployed on a personal computer. VGEA was built on the Snakemake workflow management system and utilizes existing tools for each step: fastp (Chen et al., 2018) for read trimming and read-level quality control, BWA (Li & Durbin, 2009) for mapping sequencing reads to the human reference genome, SAMtools (Li et al., 2009) for extracting unmapped reads and also for splitting bam files into fastq files, IVA (Hunt et al., 2015) for de novo assembly to generate contigs, shiver (Wymant et al., 2018) to pre-process reads for quality and contamination, then map to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences, SeqKit (Shen et al., 2016) for cleaning shiver assembly for QUAST, QUAST (Gurevich et al., 2013) to evaluate/assess the quality of genome assemblies and MultiQC (Ewels et al., 2016) for aggregation of the results from fastp, BWA and QUAST. Our pipeline was successfully tested and validated with SARS-CoV-2 (n = 20), HIV-1 (n = 20) and Lassa Virus (n = 20) datasets all of which have been made publicly available. VGEA is freely available on GitHub at: https://github.com/pauloluniyi/VGEA under the GNU General Public License.

Microbial metagenomic approach uncovers the first rabbit haemorrhagic disease virus genome in Sub-Saharan Africa

Rabbit Haemorrhagic Disease (RHD) causes high morbidity and mortality in rabbits and hares. Here, we report the first genomic characterization of lagovirus GI.2 virus in domestic rabbits from sub-Saharan Africa. We used an unbiased microbial metagenomic Next Generation Sequencing (mNGS) approach to diagnose the pathogen causing the suspected outbreak of RHD in Ibadan, Nigeria. The liver, spleen, and lung samples of five rabbits from an outbreak in 2 farms were analyzed. The mNGS revealed one full and two partial RHDV2 genomes on both farms. Phylogenetic analysis showed close clustering with RHDV2 lineages from Europe (98.6% similarity with RHDV2 in the Netherlands, and 99.1 to 100% identity with RHDV2 in Germany), suggesting potential importation. Subsequently, all the samples were confirmed by RHDV virus-specific RT-PCR targeting the VP60 gene with the expected band size of 398 bp for the five rabbits sampled. Our findings highlight the need for increased genomic surveillance of RHDV2 to track its origin, understand its diversity and to inform public health policy in Nigeria, and Sub-Saharan Africa.

Genetic diversity and population structure of Plasmodium falciparum in Nigeria: insights from microsatellite loci analysis

BACKGROUND

Malaria remains a public health burden especially in Nigeria. To develop new malaria control and elimination strategies or refine existing ones, understanding parasite population diversity and transmission patterns is crucial.

METHODS

In this study, characterization of the parasite diversity and structure of Plasmodium falciparum isolates from 633 dried blood spot samples in Nigeria was carried out using 12 microsatellite loci of P. falciparum. These microsatellite loci were amplified via semi-nested polymerase chain reaction (PCR) and fragments were analysed using population genetic tools.

RESULTS

Estimates of parasite genetic diversity, such as mean number of different alleles (13.52), effective alleles (7.13), allelic richness (11.15) and expected heterozygosity (0.804), were high. Overall linkage disequilibrium was weak (0.006, P < 0.001). Parasite population structure was low (Fst: 0.008-0.105, AMOVA: 0.039).

CONCLUSION

The high level of parasite genetic diversity and low population structuring in this study suggests that parasite populations circulating in Nigeria are homogenous. However, higher resolution methods, such as the 24 SNP barcode and whole genome sequencing, may capture more specific parasite genetic signatures circulating in the country. The results obtained can be used as a baseline for parasite genetic diversity and structure, aiding in the formulation of appropriate therapeutic and control strategies in Nigeria.

Polymorphisms in Plasmodium falciparum chloroquine resistance transporter (Pfcrt) and multidrug-resistant gene 1 (Pfmdr-1) in Nigerian children 10 years post-adoption of artemisinin-based combination treatments

The emergence and spread of Plasmodium falciparum parasites resistant to artemisinin derivatives and their partners in southeastern Asia threatens malaria control and elimination efforts, and heightens the need for an alternative therapy. We have explored the distribution of P. falciparum chloroquine resistance transporter (Pfcrt) and multidrug-resistant gene 1 (Pfmdr-1) haplotypes 10 years following adoption of artemisinin-based combination therapies in a bid to investigate the possible re-emergence of Chloroquine-sensitive parasites in Nigeria, and investigated the effect of these P. falciparum haplotypes on treatment outcomes of patients treated with artemisinin-based combination therapies. A total of 271 children aged <5 years with uncomplicated falciparum malaria were included in this study. Polymorphisms on codons 72-76 of the Pfcrt gene and codon 86 and 184 of Pfmdr-1 were determined using the high resolution melting assay. Of 240 (88.6%) samples successfully genotyped with HRM for Pfcrt, wildtype C72M74N75K76 (42.9%) and mutant C72I74E75T76 (53.8%) were observed. Also, wildtype N86Y184 (62.9%) and mutant N86F184 (21.1%), Y86Y184 (6.4%), and Y86F184 (0.4%) haplotypes of Pfmdr-1 were observed. Measures of responsiveness to ACTs were similar in children infected with P. falciparum crt haplotypes (C72I74E75T76 and C72M74N75K76) and major mdr-1 haplotypes (N86Y184, N86F184 and Y86Y184). Despite a 10 year gap since the malaria treatment policy changed to ACTs, over 50% of the P. falciparum parasites investigated in this study harboured the Chloroquine-resistant C72I74E75T76 haplotype, however this did not compromise the efficacy of artemisinin-based combination therapies. Should complete artemisinin resistance emerge from or spread to Nigeria, chloroquine might not be a good alternative therapy.

Polymorphisms in Plasmodium falciparum dihydropteroate synthetase and dihydrofolate reductase genes in Nigerian children with uncomplicated malaria using high-resolution melting technique

In 2005, the Nigerian Federal Ministry of Health revised the treatment policy for uncomplicated malaria with the introduction of artemisinin-based combination therapies (ACTs). This policy change discouraged the use of Sulphadoxine-pyrimethamine (SP) as the second-line treatment of uncomplicated falciparum malaria. However, SP is used as an intermittent preventive treatment of malaria in pregnancy (IPTp) and seasonal malaria chemoprevention (SMC) in children aged 3–59 months. There have been increasing reports of SP resistance especially in the non-pregnant population in Nigeria, thus, the need to continually monitor the efficacy of SP as IPTp and SMC by estimating polymorphisms in dihydropteroate synthetase (dhps) and dihydrofolate reductase (dhfr) genes associated with SP resistance. The high resolution-melting (HRM) assay was used to investigate polymorphisms in codons 51, 59, 108 and 164 of the dhfr gene and codons 437, 540, 581 and 613 of the dhps gene. DNA was extracted from 271 dried bloodspot filter paper samples obtained from children (< 5 years old) with uncomplicated malaria. The dhfr triple mutant I₅₁R₅₉N₁₀₈, dhps double mutant G₄₃₇G₅₈₁ and quadruple dhfr I₅₁R₅₉N₁₀₈ + dhps G₄₃₇ mutant haplotypes were observed in 80.8%, 13.7% and 52.8% parasites, respectively. Although the quintuple dhfr I₅₁R₅₉N₁₀₈ + dhps G₄₃₇E₅₄₀ and sextuple dhfr I₅₁R₅₉N₁₀₈ + dhps G₄₃₇E₅₄₀G₅₈₁ mutant haplotypes linked with in-vivo and in-vitro SP resistance were not detected, constant surveillance of these haplotypes should be done in the country to detect any change in prevalence.

Lassa Virus Genetics

In a pattern repeated across a range of ecological niches, arenaviruses have evolved a compact four-gene genome to orchestrate a complex life cycle in a narrow range of susceptible hosts. A number of mammalian arenaviruses cross-infect humans, often causing a life-threatening viral hemorrhagic fever. Among this group of geographically bound zoonoses, Lassa virus has evolved a unique niche that leads to significant and sustained human morbidity and mortality. As a biosafety level 4 pathogen, direct study of the pathogenesis of Lassa virus is limited by the sparse availability, high operating costs, and technical restrictions of the high-level biocontainment laboratories required for safe experimentation. In this chapter, we introduce the relationship between genome structure and the life cycle of Lassa virus and outline reverse genetic approaches used to probe and describe functional elements of the Lassa virus genome. We then review the tools used to obtain viral genomic sequences used for phylogeny and molecular diagnostics, before shifting to a population perspective to assess the contributions of phylogenetic analysis in understanding the evolution and ecology of Lassa virus in West Africa. We finally consider the future outlook and clinical applications for genetic study of Lassa virus.

Real-time Metagenomic Analysis of Undiagnosed Fever Cases Unveils a Yellow Fever Outbreak in Edo State, Nigeria

Fifty patients with unexplained fever and poor outcomes presented at Irrua Specialist Teaching Hospital (ISTH) in Edo State, Nigeria, an area endemic for Lassa fever, between September 2018 - January 2019. After ruling out Lassa fever, plasma samples from these epidemiologically-linked cases were sent to the African Centre of Excellence for Genomics of Infectious Diseases (ACEGID), Redeemer's University, Ede, Osun State, Nigeria, where we carried out metagenomic sequencing which implicated yellow fever virus (YFV) as the etiology of this outbreak. Twenty-nine of the 50 samples were confirmed positive for YFV by reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), 14 of which resulted in genome assembly. Maximum likelihood phylogenetic analysis revealed that these YFV sequences formed a tightly clustered clade more closely related to sequences from Senegal than sequences from earlier Nigerian isolates, suggesting that the YFV clade responsible for this outbreak in Edo State does not descend directly from the Nigerian YFV outbreaks of the last century, but instead reflects a broader diversity and dynamics of YFV in West Africa. Here we demonstrate the power of metagenomic sequencing for identifying ongoing outbreaks and their etiologies and informing real-time public health responses, resulting in accurate and prompt disease management and control.