BIRC6 modifies risk of invasive bacterial infection in Kenyan children
Abstract
Invasive bacterial disease is a major cause of morbidity and mortality in African children. Despite being caused by diverse pathogens, children with sepsis are clinically indistinguishable from one another. In spite of this, most genetic susceptibility loci for invasive infection that have been discovered to date are pathogen specific and are not therefore suggestive of a shared genetic architecture of bacterial sepsis. Here we utilise probabilistic diagnostic models to identify children with a high probability of invasive bacterial disease among critically unwell Kenyan children with P. falciparum parasitaemia. We construct a joint dataset including 1,445 bacteraemia cases and 1,143 severe malaria cases, and population controls, among critically unwell Kenyan children that have previously been genotyped for human genetic variation. Using these data we perform a cross-trait genome-wide association study of invasive bacterial infection, weighting cases according to their probability of bacterial disease. In doing so we identify and validate a novel risk locus for invasive infection secondary to multiple bacterial pathogens, that has no apparent effect on malaria risk. The locus identified modifies splicing of BIRC6 in stimulated monocytes, implicating regulation of apoptosis and autophagy in the pathogenesis of sepsis in Kenyan children.
Data availability
Patient level genotype and phenotype data are available via the European Genome-Phenome Archive, with accession codes EGAD00010000950 (WTCCC2: bacteraemia cases and controls) and EGAD00010000904 (MalariaGEN Consortium: severe malaria cases and controls).Full GWAS summary statistics have been deposited with the GWAS Catalog with accession code GCST90094632.Code and source data underlying each figure (and supplementary figure) are available at: https://github.com/jjgilchrist/Kenya_bacteraemia_malaria
Article and author information
Author details
Funding
Wellcome Trust (202800)
- Thomas N Williams
Wellcome Trust (098532)
- J Anthony G Scott
National Institute for Health and Care Research
- James Gilchrist
National Institute for Health and Care Research
- Alexander J Mentzer
Wellcome Trust (223253/Z/21/Z)
- James A Watson
Wellcome Trust (209265/Z/17/Z)
- Kathryn Maitland
- Thomas N Williams
Wellcome Trust (HCUZZ0)
- Adrian VS Hill
European Research Council (294557)
- Adrian VS Hill
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Human subjects: Following explanation of the study, written informed consent was obtained from the parent or guardian of each child included in the study. Ethical approval was obtained from the Kenya Medical Research Institute (KEMRI) National Scientific Steering and Research Committees (approval numbers; SCC1192 and SCC967) and the Oxford Tropical Research Ethics Committee (OxTREC, approval numbers; 020-06 and 014-01).
Copyright
© 2022, Gilchrist et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 628
- views
-
- 159
- downloads
-
- 6
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
- Microbiology and Infectious Disease
Evolution of gene expression frequently drives antibiotic resistance in bacteria. We had previously (Patel and Matange, eLife, 2021) shown that, in Escherichia coli, mutations at the mgrB locus were beneficial under trimethoprim exposure and led to overexpression of dihydrofolate reductase (DHFR), encoded by the folA gene. Here, we show that DHFR levels are further enhanced by spontaneous duplication of a genomic segment encompassing folA and spanning hundreds of kilobases. This duplication was rare in wild-type E. coli. However, its frequency was elevated in a lon-knockout strain, altering the mutational landscape early during trimethoprim adaptation. We then exploit this system to investigate the relationship between trimethoprim pressure and folA copy number. During long-term evolution, folA duplications were frequently reversed. Reversal was slower under antibiotic pressure, first requiring the acquisition of point mutations in DHFR or its promoter. Unexpectedly, despite resistance-conferring point mutations, some populations under high trimethoprim pressure maintained folA duplication to compensate for low abundance DHFR mutants. We find that evolution of gene dosage depends on expression demand, which is generated by antibiotic and exacerbated by proteolysis of drug-resistant mutants of DHFR. We propose a novel role for proteostasis as a determinant of copy number evolution in antibiotic-resistant bacteria.
-
- Evolutionary Biology
- Genetics and Genomics
Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.