Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Genetic study of late-onset neonatal sepsis reveals significant differences by sex

Timothy H. Ciesielski, Xueyi Zhang, Alessandra Tacconelli, Irja Lutsar, Vincent Meiffredy de Cabre, Emmanuel Roilides, Cinzia Ciccacci, Paola Borgiani, View ORCID ProfileWilliam K. Scott, NeoMero Consortium, Scott M. Williams, Giorgio Sirugo
doi: https://doi.org/10.1101/2021.08.12.21261209
Timothy H. Ciesielski
1The Department of Population and Quantitative Health Sciences at Case Western Reserve University School of Medicine, Cleveland, OH
2Mary Ann Swetland Center for Environmental Health at Case Western Reserve University School of Medicine, Cleveland, OH
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xueyi Zhang
1The Department of Population and Quantitative Health Sciences at Case Western Reserve University School of Medicine, Cleveland, OH
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alessandra Tacconelli
3Ospedale San Pietro Fatebenefratelli, Rome, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Irja Lutsar
4Department of Microbiology, School of Medicine, University of Tartu, Estonia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vincent Meiffredy de Cabre
5SC10-US019, INSERM, Villejuif, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emmanuel Roilides
6Laboratory of Infectious Diseases. 3rd Department of Paediatrics, School of Medicine. Aristotle University, Greece
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cinzia Ciccacci
7Unicamillus, Saint Camillus International University of Health Sciences, Rome, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paola Borgiani
8Dipartimento di Biomedicina e Prevenzione, Facolta’ di Medicina e Chirurgia, Universita’ di Tor Vergata, Rome, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
William K. Scott
9John P. Hussman Institute for Human Genomics, University of Miami, Miami, Florida
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for William K. Scott
10Consortium members listed in acknowledgements
Scott M. Williams
1The Department of Population and Quantitative Health Sciences at Case Western Reserve University School of Medicine, Cleveland, OH
11Department of Genetics and Genome Sciences, Case Western Reserve University, Cleveland, Ohio
12Cleveland Institute for Computational Biology, Cleveland OH, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: smw154@case.edu giorgio.sirugo@pennmedicine.upenn.edu
Giorgio Sirugo
13Institute of Systems Pharmacology and Translational Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
14Division of Translational Medicine and Human Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: smw154@case.edu giorgio.sirugo@pennmedicine.upenn.edu
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Late-Onset Neonatal Sepsis (LOS) is a rare, but life-threatening condition, involving widespread infection, immune disruption, organ dysfunction, and often death. Because exposure to pathogens is not 100% preventable, identifying susceptibility factors is critical to defining neonates at greater risk. Prior studies demonstrated that both genetics and infant sex influence susceptibility. We, therefore, performed a genome wide association study (GWAS) with 224 cases and 273 controls from six European countries to identify genetic risk factors. We tested for association with both autosomal and X-chromosome variants in the total sample and in the samples stratified by sex. In total 71 SNPs associated with neonatal sepsis at p<1×10−4 in at least one analysis. Most importantly, the sex stratified analyses revealed associations with multiple SNPs (28 SNPs in males and 16 in females), but none of 44 SNPs from single-sex analyses associated with sepsis in the other sex at p<0.05, and the vast majority showed significant SNP X sex interactions (17 of 28 from the male-only analyses and all 16 from the female-only study). Pathway analyses showed that NOTCH signaling is over-represented among loci identified by the analyses. Our results indicate that susceptibility to LOS is genetically sexually dimorphic, and suggest that NOTCH signaling is likely to play a role in it.

Introduction

Late-Onset Neonatal Sepsis (LOS) is a serious infection response that develops after the third day of life.1,2 In this rare disorder, the neonate develops widespread infection, diffuse and ineffective immune responses, and cardiovascular instability that leads to organ failure and death in approximately 15-30% of cases.1-3 In high income countries this condition affects approximately 2-4 of 1000 live births4, but in low and middle income countries the picture is far less clear due to worse surveillance.5,6 There are many reasons that neonatal infections may be more common in poor resource settings; hospital hygiene is less stringent and effective treatments are less available than in upper-middle and high income countries.5,6 LOS can only occur when infants are exposed to infectious agents in the first few days of life, and a variety of procedures and clinical practices can prevent these contacts. 5,7 However as virulent pathogens are not completely avoidable, it is critical to identify host susceptibility factors, learn more about the pathophysiology, and develop interventions that can help the neonates at greatest risk.

Neonatal sepsis is a phenotypically complex and heterogeneous disorder but there is evidence that both sex8-12 and genetics13 modulate risk. Males appear to be more susceptible to LOS as immune function and infection susceptibility differs between the sexes at all ages.8-12 Thus, it is possible that at least some of the etiologic processes and genetic risk factors differ by sex. Therefore, we performed a genome wide association study (GWAS) to identity variants that associate with neonatal sepsis in sex-combined and sex-stratified analyses. This approach is critical in the context of traits with sex differences 8,10,14-16, and is in line with recent consensus statements and published analytic guidelines16-20 that recommend considering the sexes independently in biomedical research. Failing to stratify by sex might impair the recognition of context dependent risk factors because the biology may differ.

Overall, we designed our analyses to embrace the etiologic heterogeneity of neonatal sepsis with the goal of better understanding pathophysiology to develop effective new interventions.

Methods

Data Source

Neonates with LOS were identified in 21 hospitals across Europe as previously described as part of the NeoMero studies.1,21 The details of Neomero1 are presented elsewhere 22. The recruitment centers were in 6 different countries (Italy, Greece, Estonia, Spain, Lithuania, and UK), and local ethics committees provided approval. Informed consent was given by parents prior to enrollment and a separate consent was given for genetic studies for this work.

Detailed information is listed on EudraCT (2011-001515-31 and 2011-001521-25) and clinicaltrials.gov (NCT01551394 and NCT01554124).21 LOS was defined as sepsis in infants between 3 and 90 days of age.21 The diagnosis was established in one of two ways: 1) culture confirmed or 2) clinical sepsis.2,3 Culture confirmed cases had both an abnormal clinical or laboratory finding and a positive bacterial culture. Clinical Sepsis cases met clinical and laboratory criteria of Goldstein et al 20053 or the 2010 EMA report2, but had a negative bacterial culture. All exclusion and inclusion criteria have been described previously22.

Controls were randomly selected from existing DNA banks in 3 of the European countries: Greece, Italy, and Spain from whom cases were recruited. Italian and Spanish controls were previously described adult population-based samples23,24. Greek controls were obtained from the Laboratory of Medical Genetics, Medical School, National & Kapodistrian, University of Athens, Greece, and were originally sent to the laboratory as part of a screening protocol for cystic fibrosis mutations. Although phenotypes were not available for the controls, these population-based controls should represent an appropriate comparator group as these individuals for the most part will not have had LOS due to its rarity and the likelihood of mortality for those who did contract it in these cohorts. Even at a prevalence of 4 per 10004 only ∼1 control would be misclassified, causing a very small drop in power. All control samples were de-identified.

Genotyping and Genomic Quality Control

DNA was extracted as previously described from 233 case blood samples (500 ul/sample)25. All samples were genotyped using the Illumina MEGA Consortium V2 Beadchip at the Hussman Institute for Human Genomics (University of Miami). Genome Studio V2011.1 (Illumina: San Diego, California) was used to generate the raw genotype calls. Standard quality control (QC) steps were performed26-28, and all of the following were removed prior to the initial association analyses: 1) SNPs missing calls in > 5% of participants, 2) participants with genotyping calls for less than 90% of the measured SNPs, 3) SNPs with minor allele frequencies < 5%, 4) participants with an excessively high or low number of heterozygous calls (> 3 SD away from the mean number of heterozygous calls for an individual) as individuals with a very low number of heterozygous calls may indicate inbreeding while a very high number may indicate sample contamination, and 5) participants who were cryptically related (pairs of samples with pihat > 0.2 were evaluated and one member of each pair was removed; this removed sample duplicates, full sibling, parent-child, and identical twin relationships from the data29. No pairs with pihat >0.2 remained after this step). It was important to eliminate relatives from the data because logistic regression analyses assume independent observations. After these steps there were 497 participants (224 cases, 273 controls) and 601876 autosomal/pseudo-autosomal SNPs available for association analyses (Figure S1). X chromosome SNPs not from the pseudo-autosomal region were handled separately as described below. After initial association analyses were performed, three additional QC checks were conducted on the top hits. SNPs were removed when there was evidence of: 1) deviation from Hardy-Weinberg Equilibrium (p < 1×10−4 among the controls), 2) differential patterns of missing SNPs between cases and controls (p <0.05 in a Fisher’s exact test), 3) evidence of SNP call errors in the raw data intensity plot from GenomeStudio (cluster separation score > 0.2).28

X chromosome SNPs were separated from the autosomal and pseudo-autosomal SNPs and QC steps described above were conducted separately in males and females. Subsequent association analyses were also conducted separately in males and females. 16664 SNPs passed the initial QC steps among the females (78 cases, 76 controls), and 15732 SNPs passed these steps among the males (146 cases, 197 controls). For the three post association QC steps: 1) the cluster separation scores were assessed as presented in the raw data (which included all participants), 2) differential patterns of missing SNPs between cases and controls was performed separately in the males and females, and 3) Hardy-Weinberg Equilibrium analyses were only performed for females.

Statistical Analyses

We conducted logistic regression analyses on the 601876 autosomal/pseudo-autosomal SNPs and because none of the associations reached the canonical threshold for genome wide significance (p < 5×10−8) we present SNPs with p < 1×10−4 in each analysis. A similar approach has been utilized in prior GWAS of Neonatal Sepsis30, and may be required in GWAS of low incidence, severe/acute illness where the collection of large sample sizes is unlikely. Additionally, this method allows for hypothesis generation in the early stages of research, where type 2 error is most costly.31 Overall, the canonical p-value threshold (p < 5×10−8) has multiple limitations in genomic discovery settings 32-36 and we present results to facilitate further evaluation. We also performed logistic regression analyses on the top hits adjusting for sex and the first 4 principal components (to adjust for population substructure). We selected this number of principal components (PCs) because only the first 4 PCs had eigenvalues above 1.5 (Figure S2).

Neonatal sepsis is more common in males and infection susceptibility differs by sex8-11, thereby providing a priori reason to stratify analyses. This approach is critical from a biological perspective8,10,14-16 in the context of a trait with sex differences, but recent analytic guidelines strongly recommend considering the sexes independently, even when there is no a priori reason to suspect etiologic differences 16-20. Thus, we separated the males and females and reanalyzed the autosomal/pseudo-autosomal SNPs within sexes. Because top hits differed by sex, we evaluated a SNP-sex interaction term for each of the top hits (logistic regression model with terms for sex, the first 4 PCs, and SNP-sex interaction term) in the sex-combined data. This allowed us to determine if there was a significant difference in patterns of SNP-outcome association between the sexes.

Descriptive statistics were calculated using SAS (Cary, NC). QC and logistic regression analyses were conducted in Plink 1.07. PCs were obtained using Plink 1.9. Sex status was assigned with genetic confirmation using Plink 1.07. SNPs were assigned a gene annotation if they were intragenic or within 10Kb of a gene. Locus zoom (http://locuszoom.org/) was used to visualize regions around the top hits. Reactome (version 69, https://reactome.org) 37,38 was used to determine if the gene lists from our analyses were enriched for genes in known physiologic pathways.

Validation of Prior Findings

Several published associations were reassessed. Srinivasan et al conducted a similar GWAS among extremely premature infants.30 The top ten SNPs from their work identified 4 regions of interest and we evaluated SNPs within 10 Kb of these loci. Additionally, a recent meta-analysis evaluated three candidate genes (TNFα, IL6, and IL10) 39, and we also evaluated SNPs that were intragenic or within 10Kb of these genes.

Results

Altogether DNA samples from 233 patients with LOS and 288 controls were obtained for genetic analysis. After QC, 224 LOS cases and 273 adult controls remained (Figure S1); these included 197 (72.2%) male controls and 146 (65.2%) male cases (Table 1). Most LOS cases were born preterm (70.5%), 44.2% were culture confirmed, and the mean birth weight was 1.9 kg (Table 2).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1. Characteristics of the 497 study participants
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2. Characteristics of the neonatal sepsis cases

Analysis of male and female participants combined

In unadjusted logistic regression analyses, 37 SNPs from the autosomes and the pseudo-autosomal regions of the X chromosome were associated with LOS at p < 1×10−4 (Table 3), but none reached the canonical genome-wide significance level (p < 5×10−8). The p-values and odds ratios were similar after adjusting for sex and the first four principal components (Table 3). In all analyses the reference allele for OR calculations was the major allele. The 37 SNPs were annotated at the gene level (Table 3) and a pathway analysis was conducted using the curated Reactome Database37,38 (https://reactome.org/). A hypergeometric overrepresentation test indicated that this gene list was significantly enriched for components of 12 physiologic pathways (FDR <0.05); all these pathways involved NOTCH signaling (Table S1).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3: Autosomal SNPs associated with neonatal sepsis among males and females combined a

Analysis stratified by sex

In the sex-stratified logistic regression analyses 28 SNPs were detected in males (Table 4) and 16 SNPs in females (Table 5)(p < 1×10−4). Results remained similar after adjusting for the first four principal components. When evaluating the top male SNPs in females, none associated with LOS at p < 0.05 (Table 4), and none of the top female SNPs associated in males (Table 5). In combined data set analyses, 17 of the 28 male SNPs had a significant sex-SNP interaction term (p < 0.05; Table 4). All 16 of the top female SNPS had a significant sex-SNP interaction term (Table 5). The sex stratified analyses identified fewer SNPs and genes of interest at p < 1×10−4, and many of the loci identified in the sex-combined analysis were also top SNPs in the males-only analysis. No pathways were significantly overrepresented among the top hit genes for the males (FDR <0.05), and pathway analysis was uninformative among the females because only one of the genes is in the Reactome database. Top SNPs from the sex stratified analyses, and SNPs in LD with these top SNPs, had sex-dependent significance patterns (Figure S3-S6). Top SNPs from the sex-combined analysis, and SNPs in LD with these top SNPs, were more significant in the males-only analyses than in the females-only analyses (Figure S7-S8).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 4: Autosomal SNPs associated with neonatal sepsis among the males a
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 5: Autosomal SNPs associated with neonatal sepsis among the females a

Analysis of X chromosome stratified by sex

No SNPs on the X chromosome significantly associated with LOS at p < 1×10−4. However, there were 13 SNPs in males and 9 SNPs in females that associated at p < 1×10−3. Similar to the pattern seen in the autosomal analysis, none of the top SNPs detected on the X chromosome in one sex replicated in the other sex (Table S2 and Table S3).

Replication of Prior Findings

Four regions of interest on chromosome 2, 4, 6, and 16 were identified in a prior GWAS of neonatal sepsis.30 From these regions, we identified 3 SNPs that associated with neonatal sepsis (p < 0.05) on chromosome 16, and 1 SNP on chromosome 6 (Table 6). Three candidate genes of interest have been evaluated in prior meta-analyses.39 In evaluation of these three genes, we identified five SNPs from IL10, and one SNP from TNFα, that associated with LOS (Table 7, Figure S9). The five associating SNPs near IL10 (intragenic +/-10kb) were in two linkage disequilibrium (LD) blocks (Figure S9). The minor alleles for rs1518111 and rs1800871, located near the transcription start site, were in strong LD (r2 = 0.9) and associated with decreased odds of neonatal sepsis (ORs 0.74 and 0.69, respectively). In contrast, the minor alleles for the 3’ SNPs, rs55705316, rs3024505, and rs3024493, were in LD (0.51 < r2 < 0.98) and were associated with increased odds of LOS (ORs 1.47 to 1.60).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 6: Reassessment of regions of interest from a prior GWAS a
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 7: Reassessment of regions of interest from prior candidate gene studies a

Discussion

We identified 71 SNPs that associated with LOS at p < 1×10−4 in at least one of five regression analyses: 1) sex-combined-autosomal, 2) males-only-autosomal, 3) females-only-autosomal, 4) males-only-X-chromosomal, and 5) females-only-X-chromosomal. However, none were significant at the canonical genome-wide significance level. Overall, this was not unexpected as the sample size was modest, the outcome was complex, and exposure to infectious agents varied. Nonetheless, many biologically meaningful association signals will not reach this canonical threshold 40,41 and requiring this level of significance can be problematic for discovery.31,36 Thus, the top hits in these analyses may still be of interest if they are corroborated with diverse convergent evidence.36 In other words, statistical evidence from high dimensional observational epidemiology (e.g. GWAS based associations) is of greater interest when informatic evidence (e.g. pathway analyses in curated biological databases) and laboratory experimentation (e.g. transgenic knockout models or in vitro biochemical experiments) can validate the findings. Here we use the DiCE metric36 to characterize the diversity and strength of validation for the biological signals. Overall, we observed three key patterns that have implications for further research: 1) links between LOS and IL10 and two uncharacterized genes on chromosomes 16 and 6; 2) the pattern of genes identified in the sex-combined-autosomes analysis implicates NOTCH signaling; and 3) many of the observed associations differed compellingly by sex.

Comparisons to prior genetic analyses of neonatal sepsis

Although our analyses did not reveal genome-wide significant associations, they were still useful for replicating prior findings. However, the existing studies evaluated males and females together, and this approach should be expected to detect variants that modify susceptibility in both sexes in the same direction, but not variants that modify susceptibility only in one sex. Hence, we compared the existing literature to only the findings from our sex-combined analysis because our results indicate that a precision medicine approach to LOS will necessitate stratifying by sex.

Most relevant literature in LOS comes from candidate gene studies, and three of the most studied genes (TNFα, IL6, and IL10) were recently evaluated in a meta-analysis. Evidence indicated trends for TNFα and IL10, but only IL10 reached significance.39 We identified a similar pattern in our combined male/female analysis: no IL6 SNPs associated with neonatal sepsis but 5 IL10 SNPs, and 1 TNFα SNP associated at p <0.05 (Table 7). This convergence with prior meta-analyses serves to validate both our approach and novel findings. Importantly, we identified two LD blocks among the five significant SNPs in the IL10 region. The minor alleles for rs1518111 and rs1800871 are located near the transcription start site and were associated with decreased odds of neonatal sepsis, while the minor alleles for the three SNPs near the 3’ end of gene (rs55705316, rs3024505, and rs3024493) associated with increased odds of neonatal sepsis. These findings reinforce the putative importance of IL10 in LOS and they may be used to develop testable molecular hypotheses for animal or in vitro models.

In addition to prior work with candidate genes, at least one GWAS has been conducted on neonatal sepsis.30 This study was restricted to extremely premature infants and it did not include sex stratified analyses, so the comparisons to our study are limited. Yet, we found corroborating evidence in our sex-combined analysis for 2 of the 4 loci identified (Table 6). In short, we identified 3 significant SNPs (p<0.05) in LOC105371393 on chromosome 16, overlapping with the prior study. Given the variation in the phenotypes between our and the previous GWAS, the results indicate that the genes in common are likely to be strong modulators of sepsis risk operating across sub-phenotypes. While this locus remains uncharacterized molecularly, the association with neonatal sepsis in an independent GWAS is intriguing. Thus, this locus currently has a DiCE score of 3 based on omic/observational evidence alone.36 Because no informatic or laboratory data yet exists for this locus, 4 is the maximum possible score at this stage in the research process; thus, additional research on this region is warranted. We also identified a significant SNP (p<0.05) in LOC107986666 of chromosome 6. Again, this uncharacterized locus has a DiCE score of 3 in a setting where a maximum of 4 points are possible, thus further investigation of this region is also warranted. Details of the DiCE metric calculation are available36, but in brief, the current DiCE score of 3 was achieved as follows: 3 points for statistical associations observed in two separate GWAS, 0 points for informatic evidence, and 0 points for experimental evidence.

NOTCH signaling in Neonatal Sepsis

In pathway analyses the gene list from the sex-combined autosome analysis was significantly enriched for components of 12 physiologic pathways all of which implicated NOTCH signaling (Table S1). Overall, this supports further research into the role of NOTCH signaling in sepsis, as well as, the existing small-molecule modulators of NOTCH signaling.42-49

Sex Differences in Neonatal Sepsis

None of the top SNPs in males were associated with sepsis in the females (at p<0.05), and none of the top SNPs in females associated with sepsis in males (at p<0.05). These statements hold for both the autosomal and X chromosomal variants. This pattern is quite striking and was corroborated with formal statistical interaction analyses in the sex-combined-autosomes data. Seventeen of the 28 hits among the males, and 16 of the 16 hits among the females, significantly interacted with sex. Furthermore, all 16 of the sex*SNP interaction terms among the females were highly significant (p<1×10−3). This pattern of interactions is quite remarkable given the modest size of our dataset; the probability of having all these interaction terms be significant is of the order of 10−48. Additionally, the locus zoom plots reveal that the SNPs in LD with the top hit SNPs from the stratified analyses have sex dependent significance patterns, corroborating patterns observed for the top hits. We note that we could not combine males and females to evaluate interaction terms for the X-chromosome SNPs as such analyses are biologically and analytically incompatible. Lastly, we found that associations among the males had excess influence in driving the results of the sex-combined analysis. Specifically, we observed that the identified genes of interest, and the SNP significance patterns were similar in the sex-combined and males-only analyses (Figure S5-S8). This finding is not surprising as males are more susceptible to sepsis8-11 and there were more males in our study. However, males may have also been overrepresented in prior studies of neonatal sepsis that did not stratify by sex.

Overall, the prominent sex differences we observed indicate that some components of the etiology of LOS are different in males and females, and future research in this area should always consider the sexes separately. The existence of sex specific pathophysiologies in neonatal sepsis 8-11 was proposed over 50 years ago, but this is often ignored in neonatal sepsis research. Failing to attend to this may have hampered prior studies due to etiological heterogeneity. However, the sex of the neonates is nearly always recorded, and thus this heterogeneity is retrospectively addressable. In summary, our findings suggest that sex specific analyses should be conducted routinely in neonatal sepsis research.

Weaknesses and Strengths

Several study weaknesses and strengths deserve mention. One weakness was that we were not able to obtain controls from neonates in the six European countries, but we obtained adult controls from three of the six countries (Greece, Italy, and Spain). However, this control selection strategy is not likely to generate meaningful misclassification or bias because the incidence of neonatal sepsis in high income countries is 2-4 per 1000 live births.4 Thus, we would expect that no more than 1 of the 273 controls to be misclassified (273 × 0.004 = 1.092). A second weakness of this study is its small sample size (497 participants and 224 cases) and we would expect this to reduce power to both replicate prior findings and to find novel associations. However, despite limited power, we have identified consistent, significant patterns that corroborate and extend prior findings. A strength and weakness is that all of our samples were from Europe. The relatively homogenous study populations did not show substantial ancestry variation and this should increase the power and validity of our analyses among Europeans but may reduce generalizability to other ancestry groups50. Finally, we emphasize that our sex stratification approach allowed us to address a critical but underappreciated source of etiologic heterogeneity that has hindered neonatal sepsis research.

Conclusions

We studied the role of genetics in LOS using a GWAS strategy and made several key observations. First, our sex-combined association analysis identified 37 SNPs of interest and corroborated previously identified links between neonatal sepsis and IL10 as well for two uncharacterized loci from chromosomes 6 and 16. Second, our pathway analyses motivates exploration of interventions that utilize small molecule modulators of NOTCH signaling. Third, our sex stratified analyses revealed 28 SNPs of interest among only males and 16 non-overlapping SNPs of interest among only females. Most of the top male SNPs and all the top female SNPs demonstrated statistically significant interactions with sex. Overall, these findings add to the growing evidence that male and female infants differ in their physiologic response to infectious disease, and they should routinely be considered separately in future studies of neonatal sepsis.

Data Availability

Data is not generally available due to ethical considerations but may be requested from the authors.

Funding

This work was supported by the European Commission under the FP7 (ID: 242146).

NeoMero Consortium

Jean Pierre Aboulker, Oguz Akbas, Antonella Allegro, Cinzia Auriti, Chiara Bertaina, Davide Bilardi, Giulia Bonatti, Fuat Emre Canpolat, Francesca Ippolita Calo Carducci, Corine Chazallon, Nijole Drazdiene, Susanna Esposito, Silvia Faggion, Isabelle Fournier, Eva Germovsek, Carlo Giaquinto, Genny Gottardi, Hayriye Gözde, Tiziana Grossele, Maarja Hallik, Cristina Haass, Paul Heath, Tatiana Munera Huertas, Valentina Ierardi, Mari-Liis Ilmoja, Elias Iosifidis, Sandrine Kahi, Paraskevi Karagianni, Aspasia Katragkou, Eve Kaur, Birgit Kiilaspää, Karin Kipper, Aggeliki Kontou, Victoria Kougia, Kanmaz Kutman, Jelena Kuznetsova, Elisabetta Lolli, Tuuli Metsvaht, Laurence Meyer, George Mitsiakos, Valentina Montinaro, Fabio Mosca, Makis Mylonas, Kader Ben Abdelkader Emmanuelle Netzer, Clarissa Oeser, Felix Omenaca, Zoi Dorothea Pana, Maria Luisa Paoloni, Simona Perniciaro, Laura Picault, Carlo Pietrasanta, Lorenza Pugni, Andrea Ronchi, Paolo Rossi, Suzan Şahin, Yacine Saidi, Laura Sanchez, Kosmas Sarafidis, Michael Sharland, Marina Spinelli, Joseph Standing, Claudia Tagliabue, Tuuli Tammekunn, Nina Tiburzi, Ursula Trafojer, Vytautas Usonis, Adilia Warris.

Supplemental Materials

Figure S1:
  • Download figure
  • Open in new tab
Figure S1: Quality Control Procedures Flow Chart
Figure S2:
  • Download figure
  • Open in new tab
Figure S2: Scree Plot of Eigenvalues for the Principal Components
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S1: Reactome based pathway analysis for genes identified in the combined sex autosome analysis (Table 1)

12 physiologic pathways are better represented in the gene list than would be expected by chance (FDR <0.05). All of these pathways involve NOTCH signaling. Note that BRINP3, CSMD1, SASH1 were not in the curated Reactome database and this analysis was therefore based on the remaining gene hits: KCNH6, TRAPPC9, PIK3C2G, MAP2K4, KIF13A, HEY2, SCOC, SLC32A1, ARFGAP3, TNRC6C. https://reactome.org/PathwayBrowser/#TOOL=AT

Figure S3.
  • Download figure
  • Open in new tab
Figure S3. Locus Zoom Plots displaying p-values and LD structure around the top hit in the females-only analysis
Figure S4.
  • Download figure
  • Open in new tab
Figure S4. Locus Zoom Plots displaying p-values and LD structure around the second hit in the females-only analysis
Figure S5.
  • Download figure
  • Open in new tab
Figure S5. Locus Zoom Plots displaying p-values and LD structure around the top hit in the males-only analysis
Figure S6.
  • Download figure
  • Open in new tab
Figure S6. Locus Zoom Plots displaying p-values and LD structure around the second hit in the males-only analysis
Figure S7.
  • Download figure
  • Open in new tab
Figure S7. Locus Zoom Plots displaying p-values and LD structure around the top hit in the sex-combined analysis
Figure S8.
  • Download figure
  • Open in new tab
Figure S8. Locus Zoom Plots displaying p-values and LD structure around the second hit in the sex-combined analysis
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S2: X chromosome SNPs associated with neonatal sepsis among the males a
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S3: X chromosome SNPs associated with neonatal sepsis among the females a
Figure S9:
  • Download figure
  • Open in new tab
Figure S9: Haploview plot of the LD structure for the IL-10 SNPs with the relative genomic location overlaid from the Locus Zoom Plot.

LD structure is displayed with r2 values and ORs for SNPs with significant associations in the sex-combined autosomal analysis (p<0.05) are presented above their location on the plot. In all analyses the reference allele for OR calculations is the major allele. Thus, rs55705316, rs3024505, and rs3024493 are in LD, and their minor alleles are associated with increased odds of neonatal sepsis. Additionally, rs1518111 and rs1800871 are in LD, and their minor alleles are associated with decreased odds of neonatal sepsis.

References

  1. 1.↵
    Lutsar I, Trafojer UM, Heath PT, et al. Meropenem vs standard of care for treatment of late onset sepsis in children of less than 90 days of age: study protocol for a randomised controlled trial. Trials 2011;12:215. doi: 10.1186/1745-6215-12-215 [published Online First: 2011/10/01]
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    Rossi P, Botgros R, Tibby S, et al. Report on the Expert Meeting on Neonatal and Paediatric Sepsis European Medicines Agency 2010.
  3. 3.↵
    Goldstein B, Giroir B, Randolph A. International pediatric sepsis consensus conference: definitions for sepsis and organ dysfunction in pediatrics. Pediatric critical care medicine: a journal of the Society of Critical Care Medicine and the World Federation of Pediatric Intensive and Critical Care Societies 2005;6(1):2–8. doi: 10.1097/01.Pcc.0000149131.72248.E6 [published Online First: 2005/01/08]
    OpenUrlCrossRef
  4. 4.↵
    Baltimore RS. Neonatal sepsis: epidemiology and management. Paediatric drugs 2003;5(11):723–40. doi: 10.2165/00148581-200305110-00002 [published Online First: 2003/10/29]
    OpenUrlCrossRefPubMed
  5. 5.↵
    Zaidi AK, Huskins WC, Thaver D, et al. Hospital-acquired neonatal infections in developing countries. Lancet 2005;365(9465):1175–88. doi: 10.1016/s0140-6736(05)71881-x [published Online First: 2005/03/30]
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    Qazi SA, Stoll BJ. Neonatal sepsis: a major global public health challenge. The Pediatric infectious disease journal 2009;28(1 Suppl):S1–2. doi: 10.1097/INF.0b013e31819587a9 [published Online First: 2009/01/10]
    OpenUrlCrossRefPubMed
  7. 7.↵
    Garland JS, Uhing MR. Strategies to prevent bacterial and fungal infection in the neonatal intensive care unit. Clinics in perinatology 2009;36(1):1–13. doi: 10.1016/j.clp.2008.09.005 [published Online First: 2009/01/24]
    OpenUrlCrossRefPubMed
  8. 8.↵
    Washburn TC, Medearis DN, Jr.., Childs B. SEX DIFFERENCES IN SUSCEPTIBILITY TO INFECTIONS. Pediatrics 1965;35:57–64. [published Online First: 1965/01/01]
    OpenUrlAbstract/FREE Full Text
  9. 9.
    Dutta S, Reddy R, Sheikh S, et al. Intrapartum antibiotics and risk factors for early onset sepsis. Archives of disease in childhood Fetal and neonatal edition 2010;95(2):F99–103. doi: 10.1136/adc.2009.163220 [published Online First: 2009/12/10]
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    Jaillon S, Berthenet K, Garlanda C. Sexual Dimorphism in Innate Immunity. Clinical reviews in allergy & immunology 2017 doi: 10.1007/s12016-017-8648-x [published Online First: 2017/10/01]
    OpenUrlCrossRefPubMed
  11. 11.↵
    van Lunzen J, Altfeld M. Sex differences in infectious diseases-common but neglected. J Infect Dis 2014;209 Suppl 3:S79–80. doi: 10.1093/infdis/jiu159 [published Online First: 2014/06/27]
    OpenUrlCrossRefPubMed
  12. 12.↵
    Muenchhoff M, Goulder PJ. Sex differences in pediatric infectious diseases. J Infect Dis 2014;209 Suppl 3:S120–6. doi: 10.1093/infdis/jiu232 [published Online First: 2014/06/27]
    OpenUrlCrossRefPubMed
  13. 13.↵
    Bizzarro MJ, Jiang Y, Hussain N, et al. The impact of environmental and genetic factors on neonatal late-onset sepsis. The Journal of pediatrics 2011;158(2):234-8.e1. doi: 10.1016/j.jpeds.2010.07.060 [published Online First: 2010/09/21]
    OpenUrlCrossRefPubMed
  14. 14.↵
    Klein SL, Flanagan KL. Sex differences in immune responses. Nature reviews Immunology 2016;16(10):626–38. doi: 10.1038/nri.2016.90 [published Online First: 2016/08/23]
    OpenUrlCrossRefPubMed
  15. 15.
    Koch MA. Sex Bias in Sepsis. Cell host & microbe 2018;24(5):613–15. doi: 10.1016/j.chom.2018.10.014 [published Online First: 2018/11/16]
    OpenUrlCrossRef
  16. 16.↵
    Powers MS, Smith PH, McKee SA, et al. From sexless to sexy: Why it is time for human genetics to consider and report analyses of sex. Biology of sex differences 2017;8:15. doi: 10.1186/s13293-017-0136-8 [published Online First: 2017/05/06]
    OpenUrlCrossRef
  17. 17.
    Accounting for sex in the genome. Nat Med 2017;23(11):1243. doi: 10.1038/nm.4445 [published Online First: 2017/11/09]
    OpenUrlCrossRef
  18. 18.
    Heidari S, Babor TF, De Castro P, et al. Sex and Gender Equity in Research: rationale for the SAGER guidelines and recommended use. Research integrity and peer review 2016;1:2. doi: 10.1186/s41073-016-0007-6 [published Online First: 2016/05/03]
    OpenUrlCrossRefPubMed
  19. 19.
    Khramtsova EA, Davis LK, Stranger BE. The role of sex in the genomics of human complex traits. Nat Rev Genet 2019;20(3):173–90. doi: 10.1038/s41576-018-0083-1 [published Online First: 2018/12/26]
    OpenUrlCrossRefPubMed
  20. 20.↵
    Clayton JA. Applying the new SABV (sex as a biological variable) policy to research and clinical care. Physiology & behavior 2018;187:2–5. doi: 10.1016/j.physbeh.2017.08.012 [published Online First: 2017/08/22]
    OpenUrlCrossRefPubMed
  21. 21.↵
    Germovsek E, Lutsar I, Kipper K, et al. Plasma and CSF pharmacokinetics of meropenem in neonates and young infants: results from the NeoMero studies. The Journal of antimicrobial chemotherapy 2018;73(7):1908–16. doi: 10.1093/jac/dky128 [published Online First: 2018/04/24]
    OpenUrlCrossRef
  22. 22.↵
    Lutsar I, Chazallon C, Trafojer U, et al. Meropenem vs standard of care for treatment of neonatal late onset sepsis (NeoMero1): A randomised controlled trial. PloS one 2020;15(3):e0229380. doi: 10.1371/journal.pone.0229380 [published Online First: 2020/03/05]
    OpenUrlCrossRef
  23. 23.↵
    Corte-Real HB, Macaulay VA, Richards MB, et al. Genetic diversity in the Iberian Peninsula determined from mitochondrial sequence analysis. Ann Hum Genet 1996;60(4):331–50. doi: 10.1111/j.1469-1809.1996.tb01196.x [published Online First: 1996/07/01]
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    Raveane A, Aneli S, Montinaro F, et al. Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe. Sci Adv 2019;5(9):eaaw3492. doi: 10.1126/sciadv.aaw3492 [published Online First: 2019/09/14]
    OpenUrlFREE Full Text
  25. 25.↵
    Malferrari G, Monferini E, DeBlasio P, et al. High-quality genomic DNA from human whole blood and mononuclear cells. Biotechniques 2002;33(6):1228–30. doi: 10.2144/02336bm09 [published Online First: 2002/12/31]
    OpenUrlCrossRefPubMed
  26. 26.↵
    Marees AT, de Kluiver H, Stringer S, et al. A tutorial on conducting genome-wide association studies: Quality control and statistical analysis. International journal of methods in psychiatric research 2018;27(2):e1608. doi: 10.1002/mpr.1608 [published Online First: 2018/02/28]
    OpenUrlCrossRefPubMed
  27. 27.
    Turner S, Armstrong LL, Bradford Y, et al. Quality control procedures for genome-wide association studies. Current protocols in human genetics 2011;Chapter 1:Unit1.19. doi: 10.1002/0471142905.hg0119s68 [published Online First: 2011/01/15]
    OpenUrlCrossRef
  28. 28.↵
    Igo RP, Jr.., Cooke Bailey JN, Romm J, et al. Quality Control for the Illumina HumanExome BeadChip. Curr Protoc Hum Genet 2016;90:2.14.1-2.14.16. doi: 10.1002/cphg.15 [published Online First: 2016/07/02]
    OpenUrlCrossRef
  29. 29.↵
    Malinowski J, Goodloe R, Brown-Gentry K, et al. Cryptic relatedness in epidemiologic collections accessed for genetic association studies: experiences from the Epidemiologic Architecture for Genes Linked to Environment (EAGLE) study and the National Health and Nutrition Examination Surveys (NHANES). Front Genet 2015;6:317. doi: 10.3389/fgene.2015.00317 [published Online First: 2015/11/19]
    OpenUrlCrossRef
  30. 30.↵
    Srinivasan L, Page G, Kirpalani H, et al. Genome-wide association study of sepsis in extremely premature infants. Archives of disease in childhood Fetal and neonatal edition 2017;102(5):F439–f45. doi: 10.1136/archdischild-2016-311545 [published Online First: 2017/03/12]
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    Williams SM, Haines JL. Correcting away the hidden heritability. Ann Hum Genet 2011;75(3):348–50. doi: 10.1111/j.1469-1809.2011.00640.x [published Online First: 2011/04/15]
    OpenUrlCrossRefPubMed
  32. 32.↵
    Bender R, Lange S. Multiple test procedures other than Bonferroni’s deserve wider use. Bmj 1999;318(7183):600–1. doi: 10.1136/bmj.318.7183.600a [published Online First: 1999/02/26]
    OpenUrlFREE Full Text
  33. 33.
    Kraft P. Curses--winner’s and otherwise--in genetic epidemiology. Epidemiology 2008;19(5):649–51; discussion 57-8. doi: 10.1097/EDE.0b013e318181b865 [published Online First: 2008/08/16]
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.
    Sobota RS, Shriner D, Kodaman N, et al. Addressing population-specific multiple testing burdens in genetic association studies. Annals of human genetics 2015;79(2):136–47. doi: 10.1111/ahg.12095 [published Online First: 2015/02/04]
    OpenUrlCrossRefPubMed
  35. 35.
    Rothman KJ. No adjustments are needed for multiple comparisons. Epidemiology 1990;1(1):43–6. [published Online First: 1990/01/01]
    OpenUrlCrossRefPubMed
  36. 36.↵
    Ciesielski TH, Pendergrass SA, White MJ, et al. Diverse convergent evidence in the genetic analysis of complex disease: coordinating omic, informatic, and experimental evidence to better identify and validate risk factors. BioData Min 2014;7:10. doi: 10.1186/1756-0381-7-10
    OpenUrlCrossRef
  37. 37.↵
    Fabregat A, Sidiropoulos K, Garapati P, et al. The Reactome pathway Knowledgebase. Nucleic acids research 2016;44(D1):D481–7. doi: 10.1093/nar/gkv1351 [published Online First: 2015/12/15]
    OpenUrlCrossRefPubMed
  38. 38.↵
    Fabregat A, Sidiropoulos K, Viteri G, et al. Reactome pathway analysis: a high-performance in-memory approach. BMC Bioinformatics 2017;18(1):142. doi: 10.1186/s12859-017-1559-2 [published Online First: 2017/03/03]
    OpenUrlCrossRefPubMed
  39. 39.↵
    Srinivasan L, Swarr DT, Sharma M, et al. Systematic Review and Meta-analysis: Gene Association Studies in Neonatal Sepsis. Am J Perinatol 2017;34(7):684–92. doi: 10.1055/s-0036-1597132 [published Online First: 2016/12/14]
    OpenUrlCrossRef
  40. 40.↵
    Zaykin DV, Zhivotovsky LA. Ranks of genuine associations in whole-genome scans. Genetics 2005;171(2):813–23. doi: 10.1534/genetics.105.044206 [published Online First: 2005/07/16]
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    Panagiotou OA, Ioannidis JP. What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations. Int J Epidemiol 2012;41(1):273–86. doi: 10.1093/ije/dyr178 [published Online First: 2012/01/19]
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Saleem SJ, Conrad DH. Hematopoietic cytokine-induced transcriptional regulation and Notch signaling as modulators of MDSC expansion. International immunopharmacology 2011;11(7):808–15. doi: 10.1016/j.intimp.2011.03.010 [published Online First: 2011/03/24]
    OpenUrlCrossRefPubMed
  43. 43.
    Tsao PN, Wei SC, Huang MT, et al. Lipopolysaccharide-induced Notch signaling activation through JNK-dependent pathway regulates inflammatory response. Journal of biomedical science 2011;18:56. doi: 10.1186/1423-0127-18-56 [published Online First: 2011/08/17]
    OpenUrlCrossRefPubMed
  44. 44.
    Pan T, Liu Z, Yin J, et al. Notch Signaling Pathway Was Involved in Regulating Programmed Cell Death 1 Expression during Sepsis-Induced Immunosuppression. Mediators of inflammation 2015;2015:539841. doi: 10.1155/2015/539841 [published Online First: 2015/06/13]
    OpenUrlCrossRefPubMed
  45. 45.
    Skirecki T, Kawiak J, Machaj E, et al. Early severe impairment of hematopoietic stem and progenitor cells from the bone marrow caused by CLP sepsis and endotoxemia in a humanized mice model. Stem cell research & therapy 2015;6:142. doi: 10.1186/s13287-015-0135-9 [published Online First: 2015/08/15]
    OpenUrlCrossRef
  46. 46.
    Netzer C, Knape T, Kuchler L, et al. Apoptotic Diminution of Immature Single and Double Positive Thymocyte Subpopulations Contributes to Thymus Involution During Murine Polymicrobial Sepsis. Shock (Augusta, Ga) 2017;48(2):215–26. doi: 10.1097/shk.0000000000000842 [published Online First: 2017/07/15]
    OpenUrlCrossRef
  47. 47.
    Ge Y, Huang M, Ma YF. The effects of microRNA-34a regulating Notch-1/NF-kappaB signaling pathway on lipopolysaccharide-induced human umbilical vein endothelial cells. World journal of emergency medicine 2017;8(4):292–96. doi: 10.5847/wjem.j.1920-8642.2017.04.008 [published Online First: 2017/11/11]
    OpenUrlCrossRef
  48. 48.
    Yuyun X, Xi C, Qing Y, et al. Asiatic acid attenuates lipopolysaccharide-induced injury by suppressing activation of the Notch signaling pathway. Oncotarget 2018;9(19):15036–46. doi: 10.18632/oncotarget.24542 [published Online First: 2018/03/31]
    OpenUrlCrossRef
  49. 49.↵
    Bai X, He T, Liu Y, et al. Acetylation-Dependent Regulation of Notch Signaling in Macrophages by SIRT1 Affects Sepsis Development. Front Immunol 2018;9:762. doi: 10.3389/fimmu.2018.00762 [published Online First: 2018/06/06]
    OpenUrlCrossRef
  50. 50.↵
    Sirugo G, Williams SM, Tishkoff SA. The Missing Diversity in Human Genetic Studies. Cell 2019;177(1):26–31. doi: 10.1016/j.cell.2019.02.048 [published Online First: 2019/03/23]
    OpenUrlCrossRef
  51. 51.
    Quinn JA, Munoz FM, Gonik B, et al. Preterm birth: Case definition & guidelines for data collection, analysis, and presentation of immunisation safety data. Vaccine 2016;34(49):6047–56. doi: 10.1016/j.vaccine.2016.03.045 [published Online First: 2016/10/17]
    OpenUrlCrossRef
Back to top
PreviousNext
Posted August 13, 2021.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Genetic study of late-onset neonatal sepsis reveals significant differences by sex
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Genetic study of late-onset neonatal sepsis reveals significant differences by sex
Timothy H. Ciesielski, Xueyi Zhang, Alessandra Tacconelli, Irja Lutsar, Vincent Meiffredy de Cabre, Emmanuel Roilides, Cinzia Ciccacci, Paola Borgiani, William K. Scott, NeoMero Consortium, Scott M. Williams, Giorgio Sirugo
medRxiv 2021.08.12.21261209; doi: https://doi.org/10.1101/2021.08.12.21261209
Reddit logo Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Genetic study of late-onset neonatal sepsis reveals significant differences by sex
Timothy H. Ciesielski, Xueyi Zhang, Alessandra Tacconelli, Irja Lutsar, Vincent Meiffredy de Cabre, Emmanuel Roilides, Cinzia Ciccacci, Paola Borgiani, William K. Scott, NeoMero Consortium, Scott M. Williams, Giorgio Sirugo
medRxiv 2021.08.12.21261209; doi: https://doi.org/10.1101/2021.08.12.21261209

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Genetic and Genomic Medicine
Subject Areas
All Articles
  • Addiction Medicine (228)
  • Allergy and Immunology (506)
  • Anesthesia (110)
  • Cardiovascular Medicine (1245)
  • Dentistry and Oral Medicine (206)
  • Dermatology (147)
  • Emergency Medicine (282)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (534)
  • Epidemiology (10032)
  • Forensic Medicine (5)
  • Gastroenterology (500)
  • Genetic and Genomic Medicine (2467)
  • Geriatric Medicine (238)
  • Health Economics (480)
  • Health Informatics (1647)
  • Health Policy (754)
  • Health Systems and Quality Improvement (637)
  • Hematology (250)
  • HIV/AIDS (536)
  • Infectious Diseases (except HIV/AIDS) (11872)
  • Intensive Care and Critical Care Medicine (626)
  • Medical Education (253)
  • Medical Ethics (75)
  • Nephrology (268)
  • Neurology (2290)
  • Nursing (139)
  • Nutrition (352)
  • Obstetrics and Gynecology (454)
  • Occupational and Environmental Health (537)
  • Oncology (1249)
  • Ophthalmology (377)
  • Orthopedics (134)
  • Otolaryngology (226)
  • Pain Medicine (158)
  • Palliative Medicine (50)
  • Pathology (325)
  • Pediatrics (734)
  • Pharmacology and Therapeutics (315)
  • Primary Care Research (282)
  • Psychiatry and Clinical Psychology (2281)
  • Public and Global Health (4844)
  • Radiology and Imaging (843)
  • Rehabilitation Medicine and Physical Therapy (492)
  • Respiratory Medicine (652)
  • Rheumatology (286)
  • Sexual and Reproductive Health (241)
  • Sports Medicine (227)
  • Surgery (269)
  • Toxicology (44)
  • Transplantation (125)
  • Urology (99)