Epidemiology of Human Seasonal Coronaviruses Among People With Mild and Severe Acute Respiratory Illness in Blantyre, Malawi, 2011-2017.

BACKGROUND: The aim of this study was to characterize the epidemiology of human seasonal coronaviruses (HCoVs) in southern Malawi. METHODS: We tested for HCoVs 229E, OC43, NL63, and HKU1 using real-time polymerase chain reaction (PCR) on upper respiratory specimens from asymptomatic controls and individuals of all ages recruited through severe acute respiratory illness (SARI) surveillance at Queen Elizabeth Central Hospital, Blantyre, and a prospective influenza-like illness (ILI) observational study between 2011 and 2017. We modeled the probability of having a positive PCR for each HCoV using negative binomial models, and calculated pathogen-attributable fractions (PAFs). RESULTS: Overall, 8.8% (539/6107) of specimens were positive for ≥1 HCoV. OC43 was the most frequently detected HCoV (3.1% [191/6107]). NL63 was more frequently detected in ILI patients (adjusted incidence rate ratio [aIRR], 9.60 [95% confidence interval {CI}, 3.25-28.30]), while 229E (aIRR, 8.99 [95% CI, 1.81-44.70]) was more frequent in SARI patients than asymptomatic controls. In adults, 229E and OC43 were associated with SARI (PAF, 86.5% and 89.4%, respectively), while NL63 was associated with ILI (PAF, 85.1%). The prevalence of HCoVs was similar between children with SARI and controls. All HCoVs had bimodal peaks but distinct seasonality. CONCLUSIONS: OC43 was the most prevalent HCoV in acute respiratory illness of all ages. Individual HCoVs had distinct seasonality that differed from temperate settings.

Defective Monocyte Enzymatic Function and an Inhibitory Immune Phenotype in Human Immunodeficiency Virus-Exposed Uninfected African Infants in the Era of Antiretroviral Therapy.

BACKGROUND: Human immunodeficiency virus-exposed uninfected (HEU) infants are a rapidly expanding population in sub-Saharan Africa and are highly susceptible to encapsulated bacterial disease in the first year of life. The mechanism of this increased risk is still poorly understood. We investigated whether human immunodeficiency virus (HIV)-exposure dysregulates HEU immunity, vaccine-antibody production, and human herpes virus amplify this effect. METHODS: Thirty-four HIV-infected and 44 HIV-uninfected pregnant women were recruited into the birth cohort and observed up to 6 weeks of age; and then a subsequent 43 HIV-infected and 61 HIV-uninfected mother-infant pairs were recruited into a longitudinal infant cohort at either: 5-7 to 14-15; or 14-15 to 18-23 weeks of age. We compared monocyte function, innate and adaptive immune cell phenotype, and vaccine-induced antibody responses between HEU and HIV-unexposed uninfected (HU) infants. RESULTS: We demonstrate (1) altered monocyte phagosomal function and B-cell subset homeostasis and (2) lower vaccine-induced anti-Haemophilus influenzae type b (Hib) and anti-tetanus toxoid immunoglobulin G titers in HEU compared with HU infants. Human herpes virus infection was similar between HEU and HU infants. CONCLUSIONS: In the era of antiretroviral therapy-mediated viral suppression, HIV exposure may dysregulate monocyte and B-cell function, during the vulnerable period of immune maturation. This may contribute to the high rates of invasive bacterial disease and pneumonia in HEU infants.

Evolutionary Genomics at the Human-Environment Interface in Africa.

We report on the first meeting of SMBE in Africa. SMBE Malawi was initiated to bring together African and international researchers who use genetics or genomics to study natural systems impacted by human activities. The goals of this conference were 1) to reach a world-class standard of science with a large number of contributions from Africa, 2) to initiate exchange between African and international researchers, and 3) to identify challenges and opportunities for evolutionary genomics research in Africa. As repored, we think that we have achieved these goals and make suggestions on the way forward for African evolutionary genomics research.

Epidemiological and Molecular Characterization of an Invasive Group A Streptococcus emm32.2 Outbreak.

An emm32.2 invasive group A streptococcus (iGAS) outbreak occurred in Liverpool from January 2010 to September 2012. This genotype had not previously been identified in Liverpool, but was responsible for 32% (14/44) of all iGAS cases reported during this time period. We performed a case-case comparison of emm32.2 iGAS cases with non-emm32.2 control iGAS cases identified in the Liverpool population over the same time period to assess patient risk factors for emm32.2 iGAS infection. The emm32.2 iGAS cases were confined to the adult population. We show that homelessness, intravenous drug use, and alcohol abuse predisposed patients to emm32.2 iGAS disease; however, no obvious epidemiological linkage between the patients with emm32.2 iGAS could be identified. Comparative whole-genome sequencing analysis of emm32.2 iGAS and non-emm32.2 control isolates was also performed to identify pathogen factors which might have driven the outbreak. We identified 19 genes, five of which had previously been implicated in virulence, which were present in all of the emm32.2 iGAS isolates but not present in any of the non-emm32.2 control isolates. We report that a novel emm32.2 genotype emerged in Liverpool in 2010 and identified a specific subset of genes, which could have allowed this novel emm32.2 genotype to persist in a disadvantaged population in the region over a 3-year period.

Understanding pneumococcal serotype 1 biology through population genomic analysis.

BACKGROUND: Pneumococcus kills over one million children annually and over 90 % of these deaths occur in low-income countries especially in Sub-Saharan Africa (SSA) where HIV exacerbates the disease burden. In SSA, serotype 1 pneumococci particularly the endemic ST217 clone, causes majority of the pneumococcal disease burden. To understand the evolution of the virulent ST217 clone, we analysed ST217 whole genomes from isolates sampled from African and Asian countries. METHODS: We analysed 226 whole genome sequences from the ST217 lineage sampled from 9 African and 4 Asian countries. We constructed a whole genome alignment and used it for phylogenetic and coalescent analyses. We also screened the genomes to determine presence of antibiotic resistance conferring genes. RESULTS: Population structure analysis grouped the ST217 isolates into five sequence clusters (SCs), which were highly associated with different geographical regions and showed limited intracontinental and intercontinental spread. The SCs showed lower than expected genomic sequence, which suggested strong purifying selection and small population sizes caused by bottlenecks. Recombination rates varied between the SCs but were lower than in other successful clones such as PMEN1. African isolates showed higher prevalence of antibiotic resistance genes than Asian isolates. Interestingly, certain West African isolates harbored a defective chloramphenicol and tetracycline resistance-conferring element (Tn5253) with a deletion in the loci encoding the chloramphenicol resistance gene (cat pC194), which caused lower chloramphenicol than tetracycline resistance. Furthermore, certain genes that promote colonisation were absent in the isolates, which may contribute to serotype 1's rarity in carriage and consequently its lower recombination rates. CONCLUSIONS: The high phylogeographic diversity of the ST217 clone shows that this clone has been in circulation globally for a long time, which allowed its diversification and adaptation in different geographical regions. Such geographic adaptation reflects local variations in selection pressures in different locales. Further studies will be required to fully understand the biological mechanisms which makes the ST217 clone highly invasive but unable to successfully colonise the human nasopharynx for long durations which results in lower recombination rates.

Evidence of efficient stop codon readthrough in four mammalian genes.

Stop codon readthrough is used extensively by viruses to expand their gene expression. Until recent discoveries in Drosophila, only a very limited number of readthrough cases in chromosomal genes had been reported. Analysis of conserved protein coding signatures that extend beyond annotated stop codons identified potential stop codon readthrough of four mammalian genes. Here we use a modified targeted bioinformatic approach to identify a further three mammalian readthrough candidates. All seven genes were tested experimentally using reporter constructs transfected into HEK-293T cells. Four displayed efficient stop codon readthrough, and these have UGA immediately followed by CUAG. Comparative genomic analysis revealed that in the four readthrough candidates containing UGA-CUAG, this motif is conserved not only in mammals but throughout vertebrates with the first six of the seven nucleotides being universally conserved. The importance of the CUAG motif was confirmed using a systematic mutagenesis approach. One gene, OPRL1, encoding an opiate receptor, displayed extremely efficient levels of readthrough (∼31%) in HEK-293T cells. Signals both 5' and 3' of the OPRL1 stop codon contribute to this high level of readthrough. The sequence UGA-CUA alone can support 1.5% readthrough, underlying its importance.

GWIPS-viz: development of a ribo-seq genome browser.

We describe the development of GWIPS-viz (http://gwips.ucc.ie), an online genome browser for viewing ribosome profiling data. Ribosome profiling (ribo-seq) is a recently developed technique that provides genome-wide information on protein synthesis (GWIPS) in vivo. It is based on the deep sequencing of ribosome-protected messenger RNA (mRNA) fragments, which allows the ribosome density along all mRNA transcripts present in the cell to be quantified. Since its inception, ribo-seq has been carried out in a number of eukaryotic and prokaryotic organisms. Owing to the increasing interest in ribo-seq, there is a pertinent demand for a dedicated ribo-seq genome browser. GWIPS-viz is based on The University of California Santa Cruz (UCSC) Genome Browser. Ribo-seq tracks, coupled with mRNA-seq tracks, are currently available for several genomes: human, mouse, zebrafish, nematode, yeast, bacteria (Escherichia coli K12, Bacillus subtilis), human cytomegalovirus and bacteriophage lambda. Our objective is to continue incorporating published ribo-seq data sets so that the wider community can readily view ribosome profiling information from multiple studies without the need to carry out computational processing.

Comparative Genomic Analysis of Meningitis- and Bacteremia-Causing Pneumococci Identifies a Common Core Genome.

Streptococcus pneumoniae is a nasopharyngeal commensal that occasionally invades normally sterile sites to cause bloodstream infection and meningitis. Although the pneumococcal population structure and evolutionary genetics are well defined, it is not clear whether pneumococci that cause meningitis are genetically distinct from those that do not. Here, we used whole-genome sequencing of 140 isolates of S. pneumoniae recovered from bloodstream infection (n = 70) and meningitis (n = 70) to compare their genetic contents. By fitting a double-exponential decaying-function model, we show that these isolates share a core of 1,427 genes (95% confidence interval [CI], 1,425 to 1,435 genes) and that there is no difference in the core genome or accessory gene content from these disease manifestations. Gene presence/absence alone therefore does not explain the virulence behavior of pneumococci that reach the meninges. Our analysis, however, supports the requirement of a range of previously described virulence factors and vaccine candidates for both meningitis- and bacteremia-causing pneumococci. This high-resolution view suggests that, despite considerable competency for genetic exchange, all pneumococci are under considerable pressure to retain key components advantageous for colonization and transmission and that these components are essential for access to and survival in sterile sites.

Darned in 2013: inclusion of model organisms and linking with Wikipedia.

DARNED (DAtabase of RNa EDiting, available at http://darned.ucc.ie) is a centralized repository of reference genome coordinates corresponding to RNA nucleotides having altered templated identities in the process of RNA editing. The data in DARNED are derived from published datasets of RNA editing events. RNA editing instances have been identified with various methods, such as bioinformatics screenings, deep sequencing and/or biochemical techniques. Here we report our current progress in the development and expansion of the DARNED. In addition to novel database features the DARNED update describes inclusion of Drosophila melanogaster and Mus musculus RNA editing events and the launch of a community-based annotation in the RNA WikiProject.

The African Human Microbiome Portal: a public web portal of curated metagenomic metadata.

There is growing evidence that comprehensive and harmonized metadata are fundamental for effective public data reusability. However, it is often challenging to extract accurate metadata from public repositories. Of particular concern is the metagenomic data related to African individuals, which often omit important information about the particular features of these populations. As part of a collaborative consortium, H3ABioNet, we created a web portal, namely the African Human Microbiome Portal (AHMP), exclusively dedicated to metadata related to African human microbiome samples. Metadata were collected from various public repositories prior to cleaning, curation and harmonization according to a pre-established guideline and using ontology terms. These metadata sets can be accessed at https://microbiome.h3abionet.org/. This web portal is open access and offers an interactive visualization of 14 889 records from 70 bioprojects associated with 72 peer reviewed research articles. It also offers the ability to download harmonized metadata according to the user's applied filters. The AHMP thereby supports metadata search and retrieve operations, facilitating, thus, access to relevant studies linked to the African Human microbiome. Database URL:  https://microbiome.h3abionet.org/.

Gene expression data visualization tool on the o²S²PARC platform.

Background: The identification of differentially expressed genes and their associated biological processes, molecular function, and cellular components are important for genetic diseases studies because they present potential biomarkers and therapeutic targets. Methods: In this study, we developed an o²S²PARC template representing an interactive pipeline for the gene expression data visualization and ontologies data analysis and visualization.  To demonstrate the usefulness of the tool, we performed a case study on a publicly available dataset. Results: The tool enables users to identify the differentially expressed genes (DEGs) and visualize them in a volcano plot format. The ontologies associated with the DEGs are determined and visualized in barplots. Conclusions: The "Expression data visualization" template is publicly available on the o²S²PARC platform.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

African Genomic Medicine Portal: A Web Portal for Biomedical Applications.

Genomics data are currently being produced at unprecedented rates, resulting in increased knowledge discovery and submission to public data repositories. Despite these advances, genomic information on African-ancestry populations remains significantly low compared with European- and Asian-ancestry populations. This information is typically segmented across several different biomedical data repositories, which often lack sufficient fine-grained structure and annotation to account for the diversity of African populations, leading to many challenges related to the retrieval, representation and findability of such information. To overcome these challenges, we developed the African Genomic Medicine Portal (AGMP), a database that contains metadata on genomic medicine studies conducted on African-ancestry populations. The metadata is curated from two public databases related to genomic medicine, PharmGKB and DisGeNET. The metadata retrieved from these source databases were limited to genomic variants that were associated with disease aetiology or treatment in the context of African-ancestry populations. Over 2000 variants relevant to populations of African ancestry were retrieved. Subsequently, domain experts curated and annotated additional information associated with the studies that reported the variants, including geographical origin, ethnolinguistic group, level of association significance and other relevant study information, such as study design and sample size, where available. The AGMP functions as a dedicated resource through which to access African-specific information on genomics as applied to health research, through querying variants, genes, diseases and drugs. The portal and its corresponding technical documentation, implementation code and content are publicly available.

DARNED: a DAtabase of RNa EDiting in humans.

MOTIVATION: RNA editing is a phenomenon, which is responsible for the alteration of particular nucleotides in RNA sequences relative to their genomic templates. Recently, a large number of RNA editing instances in humans have been identified using bioinformatic screens and high-throughput experimental investigations utilizing next-generation sequencing technologies. However, the available data on RNA editing are not uniform and difficult to access. RESULTS: Here, we describe a new database DARNED (DAtabase of RNa EDiting) that provides centralized access to available published data related to RNA editing. RNA editing locations are mapped on the reference human genome. The current release of the database contains information on approximately 42 000 human genome coordinates corresponding to RNA locations that undergo RNA editing, mostly involving adenosine-to-inosine (A-to-I) substitutions. The data can be queried using a range of genomic coordinates, their corresponding functional localization in RNA molecules [Exons, Introns, CoDing Sequence (CDS) and UnTranslated Regions (UTRs)] and information regarding tissue/organ/cell sources where RNA editing has been observed. It is also possible to obtain RNA editing information for a specific gene or an RNA molecule using corresponding accession numbers. Search results provide information on the number of expressed sequence tags (ESTs) supporting edited and genomic bases, functional localization of RNA editing and existence of known single nucleotide polymorphisms (SNPs). Editing data can be explored in UCSC and Ensembl genome browsers, in conjunction with additional data provided by these popular genome browsers. DARNED has been designed for researchers seeking information on RNA editing and for the developers of novel algorithms for its prediction. AVAILABILITY: DARNED is accessible at http://darned.ucc.ie

Transmission dynamics of SARS-CoV-2 within-host diversity in two major hospital outbreaks in South Africa.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes acute, highly transmissible respiratory infection in humans and a wide range of animal species. Its rapid global spread has resulted in a major public health emergency, necessitating commensurately rapid research to improve control strategies. In particular, the ability to effectively retrace transmission chains in outbreaks remains a major challenge, partly due to our limited understanding of the virus' underlying evolutionary dynamics within and between hosts. We used high-throughput sequencing whole-genome data coupled with bottleneck analysis to retrace the pathways of viral transmission in two nosocomial outbreaks that were previously characterised by epidemiological and phylogenetic methods. Additionally, we assessed the mutational landscape, selection pressures, and diversity at the within-host level for both outbreaks. Our findings show evidence of within-host selection and transmission of variants between samples. Both bottleneck and diversity analyses highlight within-host and consensus-level variants shared by putative source-recipient pairs in both outbreaks, suggesting that certain within-host variants in these outbreaks may have been transmitted upon infection rather than arising de novo independently within multiple hosts. Overall, our findings demonstrate the utility of combining within-host diversity and bottleneck estimations for elucidating transmission events in SARS-CoV-2 outbreaks, provide insight into the maintenance of viral genetic diversity, provide a list of candidate targets of positive selection for further investigation, and demonstrate that within-host variants can be transferred between patients. Together these results will help in developing strategies to understand the nature of transmission events and curtail the spread of SARS-CoV-2.

H3ABioNet genomic medicine and microbiome data portals hackathon proceedings.

African genomic medicine and microbiome datasets are usually not well characterized in terms of their origin, making it difficult to find and extract data for specific African ethnic groups or even countries. The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for developing data portals for African genomic medicine and African microbiomes to address this and ran a hackathon to initiate their development. The two portals were designed and significant progress was made in their development during the hackathon. All the participants worked in a very synergistic and collaborative atmosphere in order to achieve the hackathon's goals. The participants were divided into content and technical teams and worked over a period of 6 days. In response to one of the survey questions of what the participants liked the most during the hackathon, 55% of the hackathon participants highlighted the familial and friendly atmosphere, the team work and the diversity of team members and their expertise. This paper describes the preparations for the portals hackathon and the interaction between the participants and reflects upon the lessons learned about its impact on successfully developing the two data portals as well as building scientific expertise of younger African researchers. Database URL: The code for developing the two portals was made publicly available in GitHub repositories: [https://github.com/codemeleon/Database; https://github.com/codemeleon/AfricanMicrobiomePortal].

Genomic Characteristics of Invasive Streptococcus pneumoniae Serotype 1 in New Caledonia Prior to the Introduction of PCV13.

Streptococcus pneumoniae serotype 1 is a common cause of global invasive pneumococcal disease. In New Caledonia, serotype 1 is the most prevalent serotype and led to two major outbreaks reported in the 2000s. The pneumococcal conjugate vaccine 13 (PCV13) was introduced into the vaccination routine, intending to prevent the expansion of serotype 1 in New Caledonia. Aiming to provide a baseline for monitoring the post-PCV13 changes, we performed a whole-genome sequence analysis on 67 serotype 1 isolates collected prior to the PCV13 introduction. To highlight the S. pneumoniae serotype 1 population structure, we performed a multilocus sequence typing (MLST) analysis revealing that NC serotype 1 consisted of 2 sequence types: ST3717 and the highly dominant ST306. Both sequence types harbored the same resistance genes to beta-lactams, macrolide, streptogramin B, fluoroquinolone, and lincosamide antibiotics. We have also identified 36 virulence genes that were ubiquitous to all the isolates. Among these virulence genes, the pneumolysin sequence presented an allelic profile associated with disease outbreaks and reduced hemolytic activity. Moreover, recombination hotspots were identified in 4 virulence genes and more notably in the cps locus (cps2L), potentially leading to capsular switching, a major mechanism of the emergence of nonvaccine types. In summary, this study represents the first overview of the genomic characteristics of S. pneumoniae serotype 1 in New Caledonia prior to the introduction of PCV13. This preliminary description represents a baseline to assess the impact of PCV13 on serotype 1 population structure and genomic diversity.

Characterization of DNA methylation in Malawian Mycobacterium tuberculosis clinical isolates.

BACKGROUND: Although Mycobacterium tuberculosis (Mtb) strains exhibit genomic homology of >99%, there is considerable variation in the phenotype. The underlying mechanisms of phenotypic heterogeneity in Mtb are not well understood but epigenetic variation is thought to contribute. At present the methylome of Mtb has not been completely characterized. METHODS: We completed methylomes of 18 Mycobacterium tuberculosis (Mtb) clinical isolates from Malawi representing the largest number of Mtb genomes to be completed in a single study using Single Molecule Real Time (SMRT) sequencing to date. RESULTS: We replicate and confirm four methylation disrupting mutations in 4 lineages of Mtb. For the first time we report complete loss of methylation courtesy of C758T (S253L) mutation in the MamB gene of Indo-oceanic lineage of Mtb. Additionally, we report a novel missense mutation G454A (G152S) in the MamA gene of the Euro-American lineage which could potentially be attributed to total disruption of methylation in the CCCAG motif but partial loss in a partner motif. Through a genomic and methylome comparative analysis with a global sample of sixteen, we report previously unknown mutations affecting the pks15/1 locus in L6 isolates. We confirm that methylation in Mtb is lineage specific although some unresolved issues still remain.