Genetic and potential antigenic evolution of influenza A(H1N1)pdm09 viruses circulating in Kenya during 2009-2018 influenza seasons.

Influenza viruses undergo rapid evolutionary changes, which requires continuous surveillance to monitor for genetic and potential antigenic changes in circulating viruses that can guide control and prevention decision making. We sequenced and phylogenetically analyzed A(H1N1)pdm09 virus genome sequences obtained from specimens collected from hospitalized patients of all ages with or without pneumonia between 2009 and 2018 from seven sentinel surveillance sites across Kenya. We compared these sequences with recommended vaccine strains during the study period to infer genetic and potential antigenic changes in circulating viruses and associations of clinical outcome. We generated and analyzed a total of 383 A(H1N1)pdm09 virus genome sequences. Phylogenetic analyses of HA protein revealed that multiple genetic groups (clades, subclades, and subgroups) of A(H1N1)pdm09 virus circulated in Kenya over the study period; these evolved away from their vaccine strain, forming clades 7 and 6, subclades 6C, 6B, and 6B.1, and subgroups 6B.1A and 6B.1A1 through acquisition of additional substitutions. Several amino acid substitutions among circulating viruses were associated with continued evolution of the viruses, especially in antigenic epitopes and receptor binding sites (RBS) of circulating viruses. Disease severity declined with an increase in age among children aged < 5 years. Our study highlights the necessity of timely genomic surveillance to monitor the evolutionary changes of influenza viruses. Routine influenza surveillance with broad geographic representation and whole genome sequencing capacity to inform on prioritization of antigenic analysis and the severity of circulating strains are critical to improved selection of influenza strains for inclusion in vaccines.

Differences in proteome perturbations caused by the Wolbachia strain wAu suggest multiple mechanisms of Wolbachia-mediated antiviral activity.

Some strains of the inherited bacterium Wolbachia have been shown to be effective at reducing the transmission of dengue virus (DENV) and other RNA viruses by Aedes aegypti in both laboratory and field settings and are being deployed for DENV control. The degree of virus inhibition varies between Wolbachia strains. Density and tissue tropism can contribute to these differences but there are also indications that this is not the only factor involved: for example, strains wAu and wAlbA are maintained at similar intracellular densities but only wAu produces strong DENV inhibition. We previously reported perturbations in lipid transport dynamics, including sequestration of cholesterol in lipid droplets, with strains wMel/wMelPop in Ae. aegypti. To further investigate the cellular basis underlying these differences, proteomic analysis of midguts was carried out on Ae. aegypti lines carrying strains wAu and wAlbA: with the hypothesis that differences in perturbations may underline Wolbachia-mediated antiviral activity. Surprisingly, wAu-carrying midguts not only showed distinct proteome perturbations when compared to non-Wolbachia carrying and wAlbA-carrying midguts but also wMel-carrying midguts. There are changes in RNA processing pathways and upregulation of a specific set of RNA-binding proteins in the wAu-carrying line, including genes with known antiviral activity. Lipid transport and metabolism proteome changes also differ between strains, and we show that strain wAu does not produce the same cholesterol sequestration phenotype as wMel. Moreover, in contrast to wMel, wAu antiviral activity was not rescued by cyclodextrin treatment. Together these results suggest that wAu could show unique features in its inhibition of arboviruses compared to previously characterized Wolbachia strains.

Genomic epidemiology of the rotavirus G2P[4] strains in coastal Kenya pre- and post-rotavirus vaccine introduction, 2012-8.

The introduction of rotavirus vaccines into the national immunization programme in many countries has led to a decline in childhood diarrhoea disease burden. Coincidentally, the incidence of some rotavirus group A (RVA) genotypes has increased, which may result from non-vaccine-type replacement. Here, we investigate the evolutionary genomics of rotavirus G2P[4] which has shown an increase in countries that introduced the monovalent Rotarix® vaccine. We examined sixty-three RVA G2P[4] strains sampled from children (aged below 13 years) admitted to Kilifi County Hospital, coastal Kenya, pre- (2012 to June 2014) and post-(July 2014 to 2018) rotavirus vaccine introduction. All the sixty-three genome sequences showed a typical DS-1-like genome constellation (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Pre-vaccine G2 sequences predominantly classified as sub-lineage IVa-3 and co-circulated with low numbers of sub-lineage IVa-1 strains, whereas post-vaccine G2 sequences mainly classified into sub-lineage IVa-3. In addition, in the pre-vaccine period, P[4] sub-lineage IVa strains co-circulated with low numbers of P[4] lineage II strains, but P[4] sub-lineage IVa strains predominated in the post-vaccine period. On the global phylogeny, the Kenyan pre- and post-vaccine G2P[4] strains clustered separately, suggesting that different virus populations circulated in the two periods. However, the strains from both periods exhibited conserved amino acid changes in the known antigenic epitopes, suggesting that replacement of the predominant G2P[4] cluster was unlikely a result of immune escape. Our findings demonstrate that the pre- and post-vaccine G2P[4] strains circulating in Kilifi, coastal Kenya, differed genetically but likely were antigenically similar. This information informs the discussion on the consequences of rotavirus vaccination on rotavirus diversity.

Phylogenomic analysis uncovers a 9-year variation of Uganda influenza type-A strains from the WHO-recommended vaccines and other Africa strains.

Genetic characterisation of circulating influenza viruses directs annual vaccine strain selection and mitigation of infection spread. We used next-generation sequencing to locally generate whole genomes from 116 A(H1N1)pdm09 and 118 A(H3N2) positive patient swabs collected across Uganda between 2010 and 2018. We recovered sequences from 92% (215/234) of the swabs, 90% (193/215) of which were whole genomes. The newly-generated sequences were genetically and phylogenetically compared to the WHO-recommended vaccines and other Africa strains sampled since 1994. Uganda strain hemagglutinin (n = 206), neuraminidase (n = 207), and matrix protein (MP, n = 213) sequences had 95.23-99.65%, 95.31-99.79%, and 95.46-100% amino acid similarity to the 2010-2020 season vaccines, respectively, with several mutated hemagglutinin antigenic, receptor binding, and N-linked glycosylation sites. Uganda influenza type-A virus strains sequenced before 2016 clustered uniquely while later strains mixed with other Africa and global strains. We are the first to report novel A(H1N1)pdm09 subclades 6B.1A.3, 6B.1A.5(a,b), and 6B.1A.6 (± T120A) that circulated in Eastern, Western, and Southern Africa in 2017-2019. Africa forms part of the global influenza ecology with high viral genetic diversity, progressive antigenic drift, and local transmissions. For a continent with inadequate health resources and where social distancing is unsustainable, vaccination is the best option. Hence, African stakeholders should prioritise routine genome sequencing and analysis to direct vaccine selection and virus control.

Anti-merozoite antibodies induce natural killer cell effector function and are associated with immunity against malaria.

Natural killer (NK) cells are potent immune effectors that can be activated via antibody-mediated Fc receptor engagement. Using multiparameter flow cytometry, we found that NK cells degranulate and release IFN-γ upon stimulation with antibody-opsonized Plasmodium falciparum merozoites. Antibody-dependent NK (Ab-NK) activity was largely strain transcending and enhanced invasion inhibition into erythrocytes. Ab-NK was associated with the successful control of parasitemia after experimental malaria challenge in African adults. In an independent cohort study in children, Ab-NK increased with age, was boosted by concurrent P. falciparum infections, and was associated with a lower risk of clinical episodes of malaria. Nine of the 14 vaccine candidates tested induced Ab-NK, including some less well-characterized antigens: P41, P113, MSP11, RHOPH3, and Pf_11363200. These data highlight an important role of Ab-NK activity in immunity against malaria and provide a potential mechanism for evaluating vaccine candidates.

Maintenance of high temporal Plasmodium falciparum genetic diversity and complexity of infection in asymptomatic and symptomatic infections in Kilifi, Kenya from 2007 to 2018.

BACKGROUND: High levels of genetic diversity are common characteristics of Plasmodium falciparum parasite populations in high malaria transmission regions. There has been a decline in malaria transmission intensity over 12 years of surveillance in the community in Kilifi, Kenya. This study sought to investigate whether there was a corresponding reduction in P. falciparum genetic diversity, using msp2 as a genetic marker. METHODS: Blood samples were obtained from children (< 15 years) enrolled into a cohort with active weekly surveillance between 2007 and 2018 in Kilifi, Kenya. Asymptomatic infections were defined during the annual cross-sectional blood survey and the first-febrile malaria episode was detected during the weekly follow-up. Parasite DNA was extracted and successfully genotyped using allele-specific nested polymerase chain reactions for msp2 and capillary electrophoresis fragment analysis. RESULTS: Based on cross-sectional surveys conducted in 2007-2018, there was a significant reduction in malaria prevalence (16.2-5.5%: P-value < 0.001), however msp2 genetic diversity remained high. A high heterozygosity index (He) (> 0.95) was observed in both asymptomatic infections and febrile malaria over time. About 281 (68.5%) asymptomatic infections were polyclonal (> 2 variants per infection) compared to 46 (56%) polyclonal first-febrile infections. There was significant difference in complexity of infection (COI) between asymptomatic 2.3 [95% confidence interval (CI) 2.2-2.5] and febrile infections 2.0 (95% CI 1.7-2.3) (P = 0.016). Majority of asymptomatic infections (44.2%) carried mixed alleles (i.e., both FC27 and IC/3D7), while FC27 alleles were more frequent (53.3%) among the first-febrile infections. CONCLUSIONS: Plasmodium falciparum infections in Kilifi are still highly diverse and polyclonal, despite the reduction in malaria transmission in the community.

Genomic Epidemiology of SARS-CoV-2 in Seychelles, 2020-2021.

Seychelles, an archipelago of 155 islands in the Indian Ocean, had confirmed 24,788 cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by the 31st of December 2021. The first SARS-CoV-2 cases in Seychelles were reported on the 14th of March 2020, but cases remained low until January 2021, when a surge was observed. Here, we investigated the potential drivers of the surge by genomic analysis of 1056 SARS-CoV-2 positive samples collected in Seychelles between 14 March 2020 and 31 December 2021. The Seychelles genomes were classified into 32 Pango lineages, 1042 of which fell within four variants of concern, i.e., Alpha, Beta, Delta and Omicron. Sporadic cases of SARS-CoV-2 detected in Seychelles in 2020 were mainly of lineage B.1 (lineage predominantly observed in Europe) but this lineage was rapidly replaced by Beta variant starting January 2021, and which was also subsequently replaced by the Delta variant in May 2021 that dominated till November 2021 when Omicron cases were identified. Using the ancestral state reconstruction approach, we estimated that at least 78 independent SARS-CoV-2 introduction events occurred in Seychelles during the study period. The majority of viral introductions into Seychelles occurred in 2021, despite substantial COVID-19 restrictions in place during this period. We conclude that the surge of SARS-CoV-2 cases in Seychelles in January 2021 was primarily due to the introduction of more transmissible SARS-CoV-2 variants into the islands.

Transmission networks of SARS-CoV-2 in Coastal Kenya during the first two waves: A retrospective genomic study.

Background: Detailed understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) regional transmission networks within sub-Saharan Africa is key for guiding local public health interventions against the pandemic. Methods: Here, we analysed 1139 SARS-CoV-2 genomes from positive samples collected between March 2020 and February 2021 across six counties of Coastal Kenya (Mombasa, Kilifi, Taita Taveta, Kwale, Tana River, and Lamu) to infer virus introductions and local transmission patterns during the first two waves of infections. Virus importations were inferred using ancestral state reconstruction, and virus dispersal between counties was estimated using discrete phylogeographic analysis. Results: During Wave 1, 23 distinct Pango lineages were detected across the six counties, while during Wave 2, 29 lineages were detected; 9 of which occurred in both waves and 4 seemed to be Kenya specific (B.1.530, B.1.549, B.1.596.1, and N.8). Most of the sequenced infections belonged to lineage B.1 (n = 723, 63%), which predominated in both Wave 1 (73%, followed by lineages N.8 [6%] and B.1.1 [6%]) and Wave 2 (56%, followed by lineages B.1.549 [21%] and B.1.530 [5%]). Over the study period, we estimated 280 SARS-CoV-2 virus importations into Coastal Kenya. Mombasa City, a vital tourist and commercial centre for the region, was a major route for virus imports, most of which occurred during Wave 1, when many Coronavirus Disease 2019 (COVID-19) government restrictions were still in force. In Wave 2, inter-county transmission predominated, resulting in the emergence of local transmission chains and diversity. Conclusions: Our analysis supports moving COVID-19 control strategies in the region from a focus on international travel to strategies that will reduce local transmission. Funding: This work was funded by The Wellcome (grant numbers: 220985, 203077/Z/16/Z, 220977/Z/20/Z, and 222574/Z/21/Z) and the National Institute for Health and Care Research (NIHR), project references: 17/63/and 16/136/33 using UK Aid from the UK government to support global health research, The UK Foreign, Commonwealth and Development Office. The views expressed in this publication are those of the author(s) and not necessarily those of the funding agencies.

Incidence of chikungunya virus infections among Kenyan children with neurological disease, 2014-2018: A cohort study.

BACKGROUND: Neurological complications due to chikungunya virus (CHIKV) infection have been described in different parts of the world, with children being disproportionately affected. However, the burden of CHIKV-associated neurological disease in Africa is currently unknown and given the lack of diagnostic facilities in routine care it is possible that CHIKV is an unrecognized etiology among children with encephalitis or other neurological illness. METHODS AND FINDINGS: We estimated the incidence of CHIKV infection among children hospitalized with neurological disease in Kilifi County, coastal Kenya. We used reverse transcriptase polymerase chain reaction (RT-PCR) to systematically test for CHIKV in cerebrospinal fluid (CSF) samples from children aged <16 years hospitalized with symptoms of neurological disease at Kilifi County Hospital between January 2014 and December 2018. Clinical records were linked to the Kilifi Health and Demographic Surveillance System and population incidence rates of CHIKV infection estimated. There were 18,341 pediatric admissions for any reason during the 5-year study period, of which 4,332 (24%) had CSF collected. The most common clinical reasons for CSF collection were impaired consciousness, seizures, and coma (47%, 22%, and 21% of all collections, respectively). After acute investigations done for immediate clinical care, CSF samples were available for 3,980 admissions, of which 367 (9.2%) were CHIKV RT-PCR positive. Case fatality among CHIKV-positive children was 1.4% (95% CI 0.4, 3.2). The annual incidence of CHIKV-associated neurological disease varied between 13 to 58 episodes per 100,000 person-years among all children <16 years old. Among children aged <5 years, the incidence of CHIKV-associated neurological disease was 77 per 100,000 person-years, compared with 20 per 100,000 for cerebral malaria and 7 per 100,000 for bacterial meningitis during the study period. Because of incomplete case ascertainment due to children not presenting to hospital, or not having CSF collected, these are likely minimum estimates. Study limitations include reliance on hospital-based surveillance and limited CSF sampling in children in coma or other contraindications to lumbar puncture, both of which lead to under-ascertainment of incidence and of case fatality. CONCLUSIONS: In this study, we observed that CHIKV infections are relatively more common than cerebral malaria and bacterial meningitis among children hospitalized with neurological disease in coastal Kenya. Given the wide distribution of CHIKV mosquito vectors, studies to determine the geographic extent of CHIKV-associated neurological disease in Africa are essential.

Cessation of exclusive breastfeeding and seasonality, but not small intestinal bacterial overgrowth, are associated with environmental enteric dysfunction: A birth cohort study amongst infants in rural Kenya.

Background: Environmental Enteric Dysfunction (EED) is a chronic intestinal inflammatory disorder of unclear aetiology prevalent amongst children in low-income settings and associated with stunting. We aimed to characterise development of EED and its putative risk factors amongst rural Kenyan infants. Methods: In a birth cohort study in Junju, rural coastal Kenya, between August 2015 and January 2017, 100 infants were each followed for nine months. Breastfeeding status was recorded weekly and anthropometry monthly. Acute illnesses and antibiotics were captured by active and passive surveillance. Intestinal function and small intestinal bacterial overgrowth (SIBO) were assessed by monthly urinary lactulose mannitol (LM) and breath hydrogen tests. Faecal alpha-1-antitrypsin, myeloperoxidase and neopterin were measured as EED biomarkers, and microbiota composition assessed by 16S sequencing. Findings: Twenty nine of the 88 participants (33%) that underwent length measurement at nine months of age were stunted (length-for-age Z score <-2). During the rainy season, linear growth was slower and LM ratio was higher. In multivariable models, LM ratio, myeloperoxidase and neopterin increased after cessation of continuous-since-birth exclusive breastfeeding. For LM ratio this only occurred during the rainy season. EED markers were not associated with antibiotics, acute illnesses, SIBO, or gut microbiota diversity. Microbiota diversified with age and was not strongly associated with complementary food introduction or linear growth impairment. Interpretation: Our data suggest that intensified promotion of uninterrupted exclusive breastfeeding amongst infants under six months during the rainy season, where rainfall is seasonal, may help prevent EED. Our findings also suggest that therapeutic strategies directed towards SIBO are unlikely to impact on EED in this setting. However, further development of non-invasive diagnostic methods for SIBO is required. Funding: This research was funded in part by the Wellcome Trust (Research Training Fellowship to RJC (103376/Z/13/Z)). EPKP was supported by the MRC/DfID Newton Fund (MR/N006259/1). JAB was supported by the MRC/DFiD/Wellcome Trust Joint Global Health Trials scheme (MR/M007367/1) and the Bill & Melinda Gates Foundation (OPP1131320). HHU was supported by the NIHR Oxford Biomedical Research Centre (IS-BRC-1215-20008).

Amplicon Sequencing as a Potential Surveillance Tool for Complexity of Infection and Drug Resistance Markers in Plasmodium falciparum Asymptomatic Infections.

BACKGROUND: Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers. METHODS: We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker. RESULTS: We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65-2.77] vs 2.22 [95% CI, 2.17-2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased. CONCLUSIONS: This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.

Targeted Amplicon Deep Sequencing for Monitoring Antimalarial Resistance Markers in Western Kenya.

Molecular surveillance of Plasmodium falciparum parasites is important to track emerging and new mutations and trends in established mutations and should serve as an early warning system for antimalarial resistance. Dried blood spots were obtained from a Plasmodium falciparum malaria survey in school children conducted across eight counties in western Kenya in 2019. Real-time PCR identified 500 P. falciparum-positive samples that were amplified at five drug resistance loci for targeted amplicon deep sequencing (TADS). The absence of important kelch 13 mutations was similar to previous findings in Kenya pre-2019, and low-frequency mutations were observed in codons 569 and 578. The chloroquine resistance transporter gene codons 76 and 145 were wild type, indicating that the parasites were chloroquine and piperaquine sensitive, respectively. The multidrug resistance gene 1 haplotypes based on codons 86, 184, and 199 were predominantly present in mixed infections with haplotypes NYT and NFT, driven by the absence of chloroquine pressure and the use of lumefantrine, respectively. The sulfadoxine-pyrimethamine resistance profile was a "superresistant" combination of triple mutations in both Pfdhfr (51I 59R 108N) and Pfdhps (436H 437G 540E), rendering sulfadoxine-pyrimethamine ineffective. TADS highlighted the low-frequency variants, allowing the early identification of new mutations, Pfmdr1 codon 199S and Pfdhfr codon 85I and emerging 164L mutations. The added value of TADS is its accuracy in identifying mixed-genotype infections and for high-throughput monitoring of antimalarial resistance markers.

Targeted amplicon deep sequencing of ama1 and mdr1 to track within-host P. falciparum diversity throughout treatment in a clinical drug trial.

Introduction: Antimalarial therapeutic efficacy studies are routinely conducted in malaria-endemic countries to assess the effectiveness of antimalarial treatment strategies. Targeted amplicon sequencing (AmpSeq) uniquely identifies and quantifies genetically distinct parasites within an infection. In this study, AmpSeq of Plasmodium falciparum apical membrane antigen 1 ( ama1), and multidrug resistance gene 1 ( mdr1), were used to characterise the complexity of infection (COI) and drug-resistance genotypes, respectively. Methods: P. falciparum-positive samples were obtained from a triple artemisinin combination therapy clinical trial conducted in 30 children under 13 years of age between 2018 and 2019 in Kilifi, Kenya. Nine of the 30 participants presented with recurrent parasitemia from day 26 (624h) onwards. The ama1 and mdr1 genes were amplified and sequenced, while msp1, msp2 and glurp data were obtained from the original clinical study. Results: The COI was comparable between ama1 and msp1, msp2 and glurp; overall, ama1 detected more microhaplotypes. Based on ama1, a stable number of microhaplotypes were detected throughout treatment until day 3. Additionally, a recrudescent infection was identified with an ama1 microhaplotype initially observed at 30h and later in an unscheduled follow-up visit. Using the relative frequencies of ama1 microhaplotypes and parasitemia, we identified a fast (<1h) and slow (>5h) clearing microhaplotype. As expected, only two mdr1 microhaplotypes (NF and NY) were identified based on the combination of amino acid polymorphisms at codons 86 and 184. Conclusions: This study highlights AmpSeq as a tool for highly-resolution tracking of parasite microhaplotypes throughout treatment and can detect variation in microhaplotype clearance estimates. AmpSeq can also identify slow-clearing microhaplotypes, a potential early sign of selection during treatment. Consequently, AmpSeq has the capability of improving the discriminatory power to distinguish recrudescences from reinfections accurately.

Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya, 2016-2018.

Background: Virus genome sequencing is increasingly utilized in epidemiological surveillance. Genomic data allows comprehensive evaluation of underlying viral diversity and epidemiology to inform control. For human rhinovirus (HRV), genomic amplification and sequencing is challenging due to numerous types, high genetic diversity and inadequate reference sequences. Methods: We developed a tiled amplicon type-specific protocol for genome amplification and sequencing on the Illumina MiSeq platform of two HRV types, A15 and A101. We then assessed added value in analyzing whole genomes relative to the VP4/2 region only in the investigation of HRV molecular epidemiology within the community in Kilifi, coastal Kenya. Results: We processed 73 nasopharyngeal swabs collected between 2016-2018, and 48 yielded at least 70% HRV genome coverage. These included all A101 samples (n=10) and 38 (60.3%) A15 samples.  Phylogenetic analysis revealed that the Kilifi A101 sequences interspersed with global A101 genomes available in GenBank collected between 1999-2016. On the other hand, our A15 sequences formed a monophyletic group separate from the global genomes collected in 2008 and 2019. An improved phylogenetic resolution was observed with the genome phylogenies compared to the VP4/2 phylogenies. Conclusions: We present a type-specific full genome sequencing approach for obtaining HRV genomic data and characterizing infections.

Characterizing the Countrywide Epidemic Spread of Influenza A(H1N1)pdm09 Virus in Kenya between 2009 and 2018.

The spatiotemporal patterns of spread of influenza A(H1N1)pdm09 viruses on a countrywide scale are unclear in many tropical/subtropical regions mainly because spatiotemporally representative sequence data are lacking. We isolated, sequenced, and analyzed 383 A(H1N1)pdm09 viral genomes from hospitalized patients between 2009 and 2018 from seven locations across Kenya. Using these genomes and contemporaneously sampled global sequences, we characterized the spread of the virus in Kenya over several seasons using phylodynamic methods. The transmission dynamics of A(H1N1)pdm09 virus in Kenya were characterized by (i) multiple virus introductions into Kenya over the study period, although only a few of those introductions instigated local seasonal epidemics that then established local transmission clusters, (ii) persistence of transmission clusters over several epidemic seasons across the country, (iii) seasonal fluctuations in effective reproduction number (Re) associated with lower number of infections and seasonal fluctuations in relative genetic diversity after an initial rapid increase during the early pandemic phase, which broadly corresponded to epidemic peaks in the northern and southern hemispheres, (iv) high virus genetic diversity with greater frequency of seasonal fluctuations in 2009-2011 and 2018 and low virus genetic diversity with relatively weaker seasonal fluctuations in 2012-2017, and (v) virus spread across Kenya. Considerable influenza virus diversity circulated within Kenya, including persistent viral lineages that were unique to the country, which may have been capable of dissemination to other continents through a globally migrating virus population. Further knowledge of the viral lineages that circulate within understudied low-to-middle-income tropical and subtropical regions is required to understand the full diversity and global ecology of influenza viruses in humans and to inform vaccination strategies within these regions.

Tracking the introduction and spread of SARS-CoV-2 in coastal Kenya.

Genomic surveillance of SARS-CoV-2 is important for understanding both the evolution and the patterns of local and global transmission. Here, we generated 311 SARS-CoV-2 genomes from samples collected in coastal Kenya between 17th March and 31st July 2020. We estimated multiple independent SARS-CoV-2 introductions into the region were primarily of European origin, although introductions could have come through neighbouring countries. Lineage B.1 accounted for 74% of sequenced cases. Lineages A, B and B.4 were detected in screened individuals at the Kenya-Tanzania border or returning travellers. Though multiple lineages were introduced into coastal Kenya following the initial confirmed case, none showed extensive local expansion other than lineage B.1. International points of entry were important conduits of SARS-CoV-2 importations into coastal Kenya and early public health responses prevented established transmission of some lineages. Undetected introductions through points of entry including imports from elsewhere in the country gave rise to the local epidemic at the Kenyan coast.

Near-Complete Genome Sequences of Eight Human Astroviruses Recovered from Diarrheal Stool Samples of Hospitalized Children in Coastal Kenya in 2019.

Here, using a sequence-independent sequencing approach (M. V. Phan, P. Hong Anh, N. Van Cuong, B. Oude Munnink, et al., Virus Evol 2:vew027, 2016, https://doi.org/10.1093/ve/vew027), we determined human astrovirus (HAstV) genome sequences from eight diarrheal stool samples collected in coastal Kenya in 2019. Phylogenetic analysis identified the following 4 genotypes: HAstV-1 (n = 4), HAstV-2 (n = 1), HAstV-3 (n = 1), and HAstV-5 (n = 2).

Endemic chikungunya fever in Kenyan children: a prospective cohort study.

BACKGROUND: Chikungunya fever (CHIKF) was first described in Tanzania in 1952. Several epidemics including East Africa have occurred, but there are no descriptions of longitudinal surveillance of endemic disease. Here, we estimate the incidence of CHIKF in coastal Kenya and describe the associated viral phylogeny. METHODS: We monitored acute febrile illnesses among 3500 children visiting two primary healthcare facilities in coastal Kenya over a 5-year period (2014-2018). Episodes were linked to a demographic surveillance system and blood samples obtained. Cross-sectional sampling in a community survey of a different group of 435 asymptomatic children in the same study location was done in 2016. Reverse-transcriptase PCR was used for chikungunya virus (CHIKV) screening, and viral genomes sequenced for phylogenetic analyses. RESULTS: We found CHIKF to be endemic in this setting, associated with 12.7% (95% CI 11.60, 13.80) of all febrile presentations to primary healthcare. The prevalence of CHIKV infections among asymptomatic children in the community survey was 0.7% (95% CI 0.22, 2.12). CHIKF incidence among children < 1 year of age was 1190 cases/100,000-person years and 63 cases/100,000-person years among children aged ≥10 years. Recurrent CHIKF episodes, associated with fever and viraemia, were observed among 19 of 170 children with multiple febrile episodes during the study period. All sequenced viral genomes mapped to the ECSA genotype albeit distinct from CHIKV strains associated with the 2004 East African epidemic. CONCLUSIONS: CHIKF may be a substantial public health burden in primary healthcare on the East African coast outside epidemic years, and recurrent infections are common.

Detection of SARS-CoV-2 variant 501Y.V2 in Comoros Islands in January 2021.

Background. Genomic data is key in understanding the spread and evolution of SARS-CoV-2 pandemic and informing the design and evaluation of interventions. However, SARS-CoV-2 genomic data remains scarce across Africa, with no reports yet from the Indian Ocean islands. Methods. We genome sequenced six SARS-CoV-2 positive samples from the first major infection wave in the Union of Comoros in January 2021 and undertook detailed phylogenetic analysis. Results. All the recovered six genomes classified within the 501Y.V2 variant of concern (also known as lineage B.1.351) and appeared to be from 2 sub-clusters with the most recent common ancestor dated 30 th Oct-2020 (95% Credibility Interval: 06 th Sep-2020 to 10 th Dec-2020). Comparison of the Comoros genomes with those of 501Y.V2 variant of concern from other countries deposited into the GISAID database revealed their close association with viruses identified in France and Mayotte (part of the Comoros archipelago and a France, Overseas Department). Conclusions. The recovered genomes, albeit few, confirmed local transmission following probably multiple introductions of the SARS-CoV-2 501Y.V2 variant of concern during the Comoros's first major COVID-19 wave. These findings demonstrate the importance of genomic surveillance and have implications for ongoing control strategies on the islands.

Enrichment approach for unbiased sequencing of respiratory syncytial virus directly from clinical samples.

Background: Nasopharyngeal samples contain higher quantities of bacterial and host nucleic acids relative to viruses; presenting challenges during virus metagenomics sequencing, which underpins agnostic sequencing protocols. We aimed to develop a viral enrichment protocol for unbiased whole-genome sequencing of respiratory syncytial virus (RSV) from nasopharyngeal samples using the Oxford Nanopore Technology (ONT) MinION platform. Methods: We assessed two protocols using RSV positive samples. Protocol 1 involved physical pre-treatment of samples by centrifugal processing before RNA extraction, while Protocol 2 entailed direct RNA extraction without prior enrichment. Concentrates from Protocol 1 and RNA extracts from Protocol 2 were each divided into two fractions; one was DNase treated while the other was not. RNA was then extracted from both concentrate fractions per sample and RNA from both protocols converted to cDNA, which was then amplified using the tagged Endoh primers through Sequence-Independent Single-Primer Amplification (SISPA) approach, a library prepared, and sequencing done. Statistical significance during analysis was tested using the Wilcoxon signed-rank test. Results: DNase-treated fractions from both protocols recorded significantly reduced host and bacterial contamination unlike the untreated fractions (in each protocol p<0.01). Additionally, DNase treatment after RNA extraction (Protocol 2) enhanced host and bacterial read reduction compared to when done before (Protocol 1). However, neither protocol yielded whole RSV genomes. Sequenced reads mapped to parts of the nucleoprotein (N gene) and polymerase complex (L gene) from Protocol 1 and 2, respectively. Conclusions: DNase treatment was most effective in reducing host and bacterial contamination, but its effectiveness improved if done after RNA extraction than before. We attribute the incomplete genome segments to amplification biases resulting from the use of short length random sequence (6 bases) in tagged Endoh primers. Increasing the length of the random nucleotides from six hexamers to nine or 12 in future studies may reduce the coverage biases.