Genomic Characterization of Multidrug-Resistant Pathogenic Enteric Bacteria from Healthy Children in Osun State, Nigeria.

Antimicrobial resistance (AMR) is responsible for the spread and persistence of bacterial infections. Surveillance of AMR in healthy individuals is usually not considered, though these individuals serve as reservoirs for continuous disease transmission. Therefore, it is essential to conduct epidemiological surveillance of AMR in healthy individuals to fully understand the dynamics of AMR transmission in Nigeria. Thirteen multidrug-resistant Citrobacter spp., Enterobacter spp., Klebsiella pneumoniae, and Escherichia coli isolated from stool samples of healthy children were subjected to whole genome sequencing (WGS) using Illumina and Oxford nanopore sequencing platforms. A bioinformatics analysis revealed antimicrobial resistance genes such as the pmrB_Y358N gene responsible for colistin resistance detected in E. coli ST219, virulence genes such as senB, and ybtP&Q, and plasmids in the isolates sequenced. All isolates harbored more than three plasmid replicons of either the Col and/or Inc type. Plasmid reconstruction revealed an integrated tetA gene, a toxin production caa gene in two E. coli isolates, and a cusC gene in K. quasivariicola ST3879, which induces neonatal meningitis. The global spread of AMR pathogenic enteric bacteria is of concern, and surveillance should be extended to healthy individuals, especially children. WGS for epidemiological surveillance will improve the detection of AMR pathogens for management and control.

Molecular characterization of non-aureus staphylococci and Mammaliicoccus from Hipposideros bats in Southwest Nigeria.

Bats are not only ecologically valuable mammals but also reservoirs of zoonotic pathogens. Their vast population, ability to fly, and inhabit diverse ecological niches could play some role in the spread of antibiotic resistance. This study investigated non-aureus staphylococci and Mammaliicoccus colonization in the Hipposideros bats at Obafemi Awolowo University, Ile-Ife, Nigeria. Pharyngeal samples (n = 23) of the insectivorous bats were analyzed, and the presumptive non-aureus staphylococcal and Mammaliicoccus isolates were confirmed by matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS). The isolates were characterized based on antibiotic susceptibility testing and whole-genome sequencing (WGS). Six bacterial genomes were assembled, and three species were identified, including Mammaliicoccus sciuri (n = 4), Staphylococcus gallinarum (n = 1), and Staphylococcus nepalensis (n = 1). All the isolates were resistant to clindamycin, while the M. sciuri and S. gallinarum isolates were also resistant to fusidic acid. WGS analysis revealed that the M. sciuri and S. gallinarum isolates were mecA-positive. In addition, the M. sciuri isolates possessed some virulence (icaA, icaB, icaC, and sspA) genes. Multi-locus sequence typing identified two new M. sciuri sequence types (STs) 233 and ST234. The identification of these new STs in a migratory mammal deserves close monitoring because previously known ST57, ST60, and ST65 sharing ack (8), ftsZ (13), glpK (14), gmk (6), and tpiA (10) alleles with ST233 and ST234 have been linked to mastitis in animals. Moreover, the broad host range of M. sciuri could facilitate the dispersal of antibiotic resistance genes. This study provides evidence of the importance of including migratory animals in monitoring the development and spread of antibiotic resistance.

Genomic Characterisation of Multidrug-Resistant Pathogenic Enteric Bacteria from healthy children in Osun State, Nigeria.

Antimicrobial resistance (AMR) has been established to be a significant driver for the persistence and spread of bacterial infections. It is, therefore, essential to conduct epidemiological surveillance of AMR in healthy individuals to understand the actual dynamics of AMR in Nigeria. Multi-drug resistant Klebsiella quasivariicola (n=1), Enterobacter hormaechei (n=1), and Escherichia coli (n=3) from stool samples of healthy children were subjected to whole genome sequencing using Illumina Nextseq1000/2000 and Oxford nanopore. Bioinformatics analysis reveals antimicrobial resistance, virulence genes, and plasmids. This pathogenic enteric bacteria harbored more than three plasmid replicons of either Col and/or Inc type associated with outbreaks and AMR resistant gene pmrB responsible for colistin resistance. Plasmid reconstruction revealed an integrated tetA gene responsible for tetracycline resistance, and caa gene responsible for toxin production in two of the E.coli isolates, and a cusC gene known to induce neonatal meningitis in the K. quasivariicola ST3879. The global spread of MDR pathogenic enteric bacteria is a worrying phenomenon, and close surveillance of healthy individuals, especially children, is strongly recommended to prevent the continuous spread and achieve the elimination and eradication of these infections. Molecular epidemiological surveillance using whole genome sequencing (WGS) will improve the detection of MDR pathogens in Nigeria.

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.

Emergence and Genomic Characterization of Multidrug Resistant Candida auris in Nigeria, West Africa.

Candida auris is an emerging multidrug-resistant fungal pathogen that has become a worldwide public health threat due to the limitations of treatment options, difficulty in diagnosis, and its potential for clonal transmission. Four ICU patients from three different healthcare facilities in Southern Nigeria presented features suggestive of severe sepsis and the blood cultures yielded the growth of Candida spp., which was identified using VITEK 2 as C. auris. Further confirmation was performed using whole genome sequencing (WGS). From the genomic analysis, two had mutations that conferred resistance to the antifungal azole group and other non-synonymous mutations in hotspot genes, such as ERG2, ERG11, and FKS1. From the phylogenetic analysis, cases 2 and 4 had a confirmed mutation (ERG11:Y132F) that conferred drug resistance to azoles clustered with clade 1, whilst cases 1 and 3 clustered with clade 4. Three of the patients died, and the fourth was most likely a case of colonization since he received no antifungals and was discharged home. These first cases of C. auris reported from Nigeria were most likely introduced from different sources. It is of public health importance as it highlights diagnostic gaps in our setting and the need for active disease surveillance in the region.

Simplified Cas13-based assays for the fast identification of SARS-CoV-2 and its variants.

The widespread transmission and evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) call for rapid nucleic acid diagnostics that are easy to use outside of centralized clinical laboratories. Here we report the development and performance benchmarking of Cas13-based nucleic acid assays leveraging lyophilised reagents and fast sample inactivation at ambient temperature. The assays, which we named SHINEv.2 (for 'streamlined highlighting of infections to navigate epidemics, version 2'), simplify the previously reported RNA-extraction-free SHINEv.1 technology by eliminating heating steps and the need for cold storage of the reagents. SHINEv.2 detected SARS-CoV-2 in nasopharyngeal samples with 90.5% sensitivity and 100% specificity (benchmarked against the reverse transcription quantitative polymerase chain reaction) in less than 90 min, using lateral-flow technology and incubation in a heat block at 37 °C. SHINEv.2 also allows for the visual discrimination of the Alpha, Beta, Gamma, Delta and Omicron SARS-CoV-2 variants, and can be run without performance losses by using body heat. Accurate, easy-to-use and equipment-free nucleic acid assays could facilitate wider testing for SARS-CoV-2 and other pathogens in point-of-care and at-home settings.

Equipment-free detection of SARS-CoV-2 and Variants of Concern using Cas13.

The COVID-19 pandemic, and the recent rise and widespread transmission of SARS-CoV-2 Variants of Concern (VOCs), have demonstrated the need for ubiquitous nucleic acid testing outside of centralized clinical laboratories. Here, we develop SHINEv2, a Cas13-based nucleic acid diagnostic that combines quick and ambient temperature sample processing and lyophilized reagents to greatly simplify the test procedure and assay distribution. We benchmarked a SHINEv2 assay for SARS-CoV-2 detection against state-of-the-art antigen-capture tests using 96 patient samples, demonstrating 50-fold greater sensitivity and 100% specificity. We designed SHINEv2 assays for discriminating the Alpha, Beta, Gamma and Delta VOCs, which can be read out visually using lateral flow technology. We further demonstrate that our assays can be performed without any equipment in less than 90 minutes. SHINEv2 represents an important advance towards rapid nucleic acid tests that can be performed in any location.

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.

Whole genome sequencing of clinical samples reveals extensively drug resistant tuberculosis (XDR TB) strains from the Beijing lineage in Nigeria, West Africa.

Multi-drug (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) continues to be a global public health problem especially in high TB burden countries like Nigeria. Many of these cases are undetected and go on to infect high risk individuals. Clinical samples from positive rifampicin resistant Xpert®MTB/Rif assay were subjected to direct whole genome sequencing and bioinformatics analysis to identify the full antibiotics resistance and lineage profile. We report two (2) XDR TB samples also belonging to the East-Asian/Beijing family of lineage 2 Mycobacterium tuberculosis complex from clinical samples in Nigeria. Our findings further reveal the presence of mutations that confer resistance to first-line drugs (rifampicin, isoniazid, ethambutol and pyrazanimide), second-line injectables (capreomycin, streptomycin, kanamycin and/or amikacin) and at least one of the fluoroquinolones (ofloxacin, moxifloxacin, levofloxacin and/or ciprofloxacin) in both samples. The genomic sequence data from this study not only provide the first evidence of XDR TB in Nigeria and West Africa, but also emphasize the importance of WGS in accurately detecting MDR and XDR TB, to ensure adequate and proper management treatment regimens for affected individuals. This will greatly aid in preventing the spread of drug resistance TB in high burden countries.

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 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 I51R59N108, dhps double mutant G437G581 and quadruple dhfr I51R59N108 + dhps G437 mutant haplotypes were observed in 80.8%, 13.7% and 52.8% parasites, respectively. Although the quintuple dhfr I51R59N108 + dhps G437E540 and sextuple dhfr I51R59N108 + dhps G437E540G581 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.

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.

Streamlined inactivation, amplification, and Cas13-based detection of SARS-CoV-2.

The COVID-19 pandemic has highlighted that new diagnostic technologies are essential for controlling disease transmission. Here, we develop SHINE (Streamlined Highlighting of Infections to Navigate Epidemics), a sensitive and specific diagnostic tool that can detect SARS-CoV-2 RNA from unextracted samples. We identify the optimal conditions to allow RPA-based amplification and Cas13-based detection to occur in a single step, simplifying assay preparation and reducing run-time. We improve HUDSON to rapidly inactivate viruses in nasopharyngeal swabs and saliva in 10 min. SHINE's results can be visualized with an in-tube fluorescent readout - reducing contamination risk as amplification reaction tubes remain sealed - and interpreted by a companion smartphone application. We validate SHINE on 50 nasopharyngeal patient samples, demonstrating 90% sensitivity and 100% specificity compared to RT-qPCR with a sample-to-answer time of 50 min. SHINE has the potential to be used outside of hospitals and clinical laboratories, greatly enhancing diagnostic capabilities.

Deployable CRISPR-Cas13a diagnostic tools to detect and report Ebola and Lassa virus cases in real-time.

Recent outbreaks of viral hemorrhagic fevers (VHFs), including Ebola virus disease (EVD) and Lassa fever (LF), highlight the urgent need for sensitive, deployable tests to diagnose these devastating human diseases. Here we develop CRISPR-Cas13a-based (SHERLOCK) diagnostics targeting Ebola virus (EBOV) and Lassa virus (LASV), with both fluorescent and lateral flow readouts. We demonstrate on laboratory and clinical samples the sensitivity of these assays and the capacity of the SHERLOCK platform to handle virus-specific diagnostic challenges. We perform safety testing to demonstrate the efficacy of our HUDSON protocol in heat-inactivating VHF viruses before SHERLOCK testing, eliminating the need for an extraction. We develop a user-friendly protocol and mobile application (HandLens) to report results, facilitating SHERLOCK's use in endemic regions. Finally, we successfully deploy our tests in Sierra Leone and Nigeria in response to recent outbreaks.

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.

Genomic Analysis of Lassa Virus during an Increase in Cases in Nigeria in 2018.

During 2018, an unusual increase in Lassa fever cases occurred in Nigeria, raising concern among national and international public health agencies. We analyzed 220 Lassa virus genomes from infected patients, including 129 from the 2017-2018 transmission season, to understand the viral populations underpinning the increase. A total of 14 initial genomes from 2018 samples were generated at Redeemer's University in Nigeria, and the findings were shared with the Nigerian Center for Disease Control in real time. We found that the increase in cases was not attributable to a particular Lassa virus strain or sustained by human-to-human transmission. Instead, the data were consistent with ongoing cross-species transmission from local rodent populations. Phylogenetic analysis also revealed extensive viral diversity that was structured according to geography, with major rivers appearing to act as barriers to migration of the rodent reservoir.