
==== Front
Medicine (Baltimore)
Medicine (Baltimore)
MD
Medicine
0025-7974
1536-5964
Lippincott Williams & Wilkins Hagerstown, MD

MD-D-24-00500
00062
10.1097/MD.0000000000039404
3
4900
Research Article
Observational Study
Causal effects of sepsis on structural changes in cerebral cortex: A Mendelian randomization investigation
https://orcid.org/0009-0005-5654-3973
Zhou Dengfeng BD a*
Wang Weina MM weinawang1225@163.com
a
Gu Jiaying MM 474426139@qq.com
a
Lu Qiaofa MM luqiaofa@163.com
a
a Department of Respiratory and Critical Care Medicine, Wuhan Fourth Hospital, Wuhan, Hubei Province, China.
* Correspondence: Dengfeng Zhou, Department of Respiratory and Critical Care Medicine, Wuhan Fourth Hospital, 76, Jiefang Avenue, Wuhan 430030, China (e-mail: 18971406617@189.cn).
06 9 2024
06 9 2024
103 36 e3940415 1 2024
23 4 2024
01 8 2024
Copyright © 2024 the Author(s). Published by Wolters Kluwer Health, Inc.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the Creative Commons Attribution License 4.0 (CCBY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous research has shown a strong correlation between sepsis and brain structure. However, whether this relationship represents a causality remains elusive. In this study, we employed Mendelian randomization (MR) to probe the associations of genetically predicted sepsis and sepsis-related death with structural changes in specific brain regions. Genome-wide association study (GWAS) data for sepsis phenotypes (sepsis and sepsis-related death) were obtained from the IEU OpenGWAS. Correspondingly, GWAS data for brain structural traits (volume of the subcortical structure, cortical thickness, and surface area) were derived from the ENIGMA consortium. Inverse variance weighted was mainly utilized to assess the causal effects, while weighted median and MR-Egger regression served as complementary methods. Sensitivity analyses were implemented with Cochran Q test, MR-Egger regression, and MR-PRESSO. In addition, a reverse MR analysis was carried out to assess the possibility of reverse causation. We identified that genetic liability to sepsis was normally significantly associated with a reduced surface area of the postcentral gyrus (β = −35.5280, SE = 13.7465, P = .0096). The genetic liability to sepsis-related death showed a suggestive positive correlation with the surface area of fusiform gyrus (β = 11.0920, SE = 3.6412, P = .0023) and posterior cingulate gyrus (β = 3.6530, SE = 1.6684, P = .0286), While it presented a suggestive negative correlation with surface area of the caudal middle frontal gyrus (β = −11.4586, SE = 5.1501, P = .0261) and frontal pole (β = −1.0024, SE = 0.4329, P = .0206). We also indicated a possible bidirectional causal association between genetic liability to sepsis-related death and the thickness of the transverse temporal gyrus. Sensitivity analyses verified the robustness of the above associations. These findings suggested that genetically determined liability to sepsis might influence the specific brain structure in a causal way, offering new perspectives to investigate the mechanism of sepsis-related neuropsychiatric disorders.

brain structures
causality
Mendelian randomization
sepsis
OPEN-ACCESSTRUE
SDCT
==== Body
pmc1. Introduction

Sepsis is one of the leading causes of death worldwide, posing a serious threat to public health.[1] According to epidemiological surveys, sepsis affected nearly 48.9 million people worldwide in 2017, with 11 million deaths, accounting for 19.7% of global fatalities.[2] Sepsis is an intrinsically complex and heterogeneous disease, characterized by multi-organ dysfunction including the brain.[3] Acute brain dysfunction caused by sepsis, known as sepsis-associated encephalopathy, is a common complication during sepsis with an incidence rate of up to 70%.[4] It is manifested as delirium and coma closely, which is associated with increased mortality rates.[4] More concerning, sepsis is a major risk factor for long-term cognitive disorders. Hospitalization for severe sepsis increases the incidence of moderate to severe cognitive dysfunction by over 10%.[5,6] Also, the incidence of mental disorders such as depression, anxiety, and posttraumatic stress disorder among sepsis survivors is significantly higher than that in the general population.[7] These findings suggested that sepsis severely impacts the quality of daily life for sepsis survivors, imposing a substantial healthcare burden on families and society. Therefore, it is of great significance to fully explore the neuropsychiatric diseases associated with sepsis to reduce the mortality and improve the survival quality of sepsis.

The brain, including the cortex and subcortical structures, is the human’s highest neural center, responsible for various functions such as movement, sensation, and cognition. Brain structural changes are thought to be the basis for functional changes, with specific brain structure patterns always observed in various neuropsychiatric diseases. A growing number of studies have shown a correlation between sepsis and brain structures.[8–11] For example, a study based on brain magnetic resonance imaging (MRI) revealed that atrophy occurred in 60% of sepsis patients with acute brain dysfunction, especially in areas like cerebral and cerebellar white matter, cortex, hippocampus, and amygdala.[8] Reductions in the volume of the caudate, putamen, and thalamus were associated with poor prognosis of sepsis.[8] Furthermore, Gunther et al showed that sepsis survivors experienced general or partial brain atrophy after discharge, with the degree of atrophy in specific brain regions being associated with reduced cognitive function, executive function, and visual attention disorders.[11] However, the causal link between sepsis and brain structures remains uncertain due to traditional observational studies are prone to biases from confounding factors and reverse causation.

Mendelian randomization (MR) analysis is a new epidemiological method utilizing genetic variants as instrumental variables (IVs) for investigating causality between exposures and outcomes.[12] Based on the law of segregation, it is assumed that genetic variants are randomly distributed to offspring before birth, akin to the random assignment in randomized clinical trials, thereby mitigating the bias of confounding factors and reverse causation.[13] So far, MR has been employed to investigate the causal effect of various exposures (such as disease exposures including schizophrenia and other neuropsychiatric disorders, lifestyle behaviors, and circulating biomarkers) on brain structures.[14–16] However, the causal relationship between sepsis and cortical brain structures has not yet been examined. In this study, we employed a MR analysis to explore the causal effects of genetic liability to sepsis and sepsis-related death on cortical and subcortical structures (34 regional cortical thickness and surface area, 7 subcortical volume). Findings from this study may provide new perspectives on the relationship between sepsis and neuropsychiatric disorders.

2. Materials and methods

2.1. Study design

This study followed the guidelines of the Strengthening the Reporting of Observational Studies in Epidemiology Using MR (STROBE-MR), a checklist of which can be found in the Supporting Information (Table S1, Supplemental Digital Content, http://links.lww.com/MD/N414).[17] Employing summary statistics derived from genome-wide association studies (GWAS), we conducted 2-sample MR analyses to evaluate the causality of sepsis, sepsis (28-day death), and brain cortical structure. IVs are chosen to meet 3 assumptions: IVs are required to be strongly associated with the exposure of interest; IVs must be independent of unmeasured confounders; and IVs influence outcomes only through the exposure of interest.[18] In this MR study, sepsis, sepsis (28-day death), and brain cortical structure, and were used as exposure and outcome, respectively. A flowchart depicting the overall design is presented in Figure 1.

Figure 1. Overall design of our research.

2.2. Sepsis samples

Summary statistics for sepsis phenotypes were extracted from the IEU OpenGWAS, while summary-level data from the UK Biobank was also utilized (https://gwas.mrcieu.ac.uk/). The UK Biobank is an extensive cohort of UK adult participants, details were found in other sections.[19] Sepsis phenotypes included sepsis (total cases: 11,643; total controls: 474,841) and sepsis-related 28-day mortality (total cases: 1896; total controls: 484,588). In the hospital episode statistics sourced from the UK Biobank, cases were identified when the relevant code was listed as the primary or secondary diagnosis. The analysis was adjusted for age, gender, microchip, and first 10 principal components using regenie v2.2.4.[20] Sepsis admissions identified by ICD codes from the UK Biobank linking secondary care data. In line with existing literature, ICD-10 codes A02, A39, A40, and A41 were used to identify sepsis.[21] The participants of the studies were all of European ancestry.

2.3. Brain cortical structure samples

Brain structural characterization data were selected from the GWAS results of MRI-derived brain morphometry performed by the ENIGMA consortium. Regarding cortical thickness and surface area, 51,665 individuals from 60 cohorts worldwide were included. Thirty-four brain regions and the entire cortex were defined using the Desikan–Killiany Cortical Atlas, and estimates were weighted for the entire brain. Regarding the volume of subcortical structures, 7 brain regions were measured for 30,717 participants. All data were adjusted for intracranial volume. Phenotype was defined as the mean estimate of the left and right hemispheres (thickness in millimeters, surface area in square millimeters, and volume of subcortical structures in cubic centimeters).[22,23]

2.4. Selection of IVs

As the eligible number of IVs was extremely small (P < 5 × 10−8), a relatively high threshold (P < 5 × 10−6) was used, which was in line with Zhou et al.[24] In parallel, we performed quality control according to the following steps: First, the linkage disequilibrium threshold was set as 0.001, and the clumping window was 10 Mb. Second, we computed the F-statistic to quantify the genetic variation strength and then abandoned single nucleotide polymorphisms (SNPs) with an F-statistic of <10, indicating insufficient strength. Third, harmonizing processes were conducted to exclude ambiguous and palindromic SNPs. Ultimately, the brain cortical structure associated with the outcome exhibited a significance level of P < 5.0 × 10−6, and brain cortical structures containing fewer than 3 SNPs were excluded from the analysis. Subsequently, the remaining SNPs were employed as IVs following the aforementioned steps.

2.5. Primary analysis

For the primary analysis, we employed random effects inverse variance weighted (IVW) estimation, which combines the Wald ratio of each SNP to the outcome, assuming all genetic variations are valid. This approach offers the highest power for MR estimation, yet it is susceptible to multidirectional bias.[25] Hence, IVW was used as the primary method to estimate the causal effect of sepsis, sepsis-related 28-day death on structural changes in specific brain regions.

2.6. Sensitivity analysis

Sensitivity analyses were subsequently performed to assess the bias of the MR assumptions for the identified significant estimates (PIVW < .05). Some other MR analyses, such as weighted median and MR-Egger regression, were also used as complementary methods. MR-Egger regressions can test for multiplicity and considerable heterogeneity of imbalances, whereas for the same change in underexposure, larger sample sizes are required.[26] In situations where at least half of the weighted variance provided by the horizontal pleated product effect is valid, the weighted median method offers consistent estimates of the effect.[27] In addition, sensitivity analyses are also critical in MR studies to evaluate any bias of the MR assumptions. Consequently, Cochran Q statistic, and MR-Egger intercept tests are employed to identify the presence of heterogeneity and pleiotropy, as well as evaluate the robustness of the obtained results. Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) tests were conducted to test for outliers with potential horizontal pleiotropy. Should any outliers be discovered, they are eliminated to produce unbiased causal estimates from the outlier-corrected MR analysis.[28] Hence, the underlying brain structure impacted by sepsis was recognized as follows: PIVW＜.05; the direction and amplitude of the 3 MR methods were consistent; no pleiotropy was observed. P-values <.05 were considered nominally significant, whereas high-confidence results were those that survived adjustment for multiple testing (Bonferroni correction threshold of 0.05/75). All analyses were run by using the R package TwoSampleMR (version 0.5.6) in R (version 4.2.3).

2.7. Confounding analysis and reverse MR analysis

Despite employing various statistical approaches in sensitivity analyses to investigate potential violations of MR assumptions, we also utilized the Phenoscanner V2 website (http://www.phenoscanner.medchsl.cam.ac.uk/) to examine whether sepsis-related SNPs were concurrently associated with multiple common risk factors that might influence MR estimates, including all neurological and psychiatric disorders, hypoxemia, obesity, intelligence, and other potential confounders.[29,30] If the correlation between SNPs and these potential confounders reached a threshold value of P < 5 × 10−6, IVW was reiterated following the removal of these SNPs to validate the robustness of the findings. Additionally, reverse MR analysis was performed on the sepsis which was found to be causally associated with brain cortical structure in forward MR analysis. The procedure for reverse MR analysis is the same as the MR analysis above.

3. Results

After a series of quality control steps, the characteristics of selected SNPs were shown in Table S2, Supplemental Digital Content, http://links.lww.com/MD/N414. F statistics of all IVs were more than 10, demonstrating the absence of weak instrument bias.

3.1. Causal effect of genetic liability to sepsis on brain structure

Our preliminary IVW analysis showed that genetic liability to sepsis had no causal effect on subcortical structures’ volume or cortical thickness (Figs. 2 and 3), but was causally associated with increased surface area of the lingual gyrus (β = 26.9884, SE = 11.2470, P = .0164) and pericalcarine (β = 19.3815, SE = 7.9255, P = .0145) (Fig. 4). By performing sensitivity analysis, 3 MR analysis methods identified inconsistent directions of the above effect (Tables S3–S5, Supplemental Digital Content, http://links.lww.com/MD/N414). Interestingly, IVW demonstrated genetic liability to sepsis was associated with reduced surface area of the postcentral gyrus (β = −35.5280, SE = 13.7465, P = .0096) after removing SNP outliers identified by the MR-PRESSO method. Cochran Q-derived P values (P > .05) suggested that there was no heterogeneity. Besides, neither the MR-Egger regression test nor the MR-PRESSO test exhibited any evidence of horizontal pleiotropy (all P > .05).

Figure 2. IVW estimates of the effect of sepsis and sepsis-related death on the volume of subcortical structures.

Figure 3. IVW estimates of the effect of sepsis and sepsis-related death on cortical thickness.

Figure 4. IVW estimates of the effect of sepsis and sepsis-related death on surface area.

3.2. Causal effect of genetic liability to sepsis-related death on brain structure

The results of IVW analysis demonstrated that there was no causal relationship between genetic liability to sepsis-related death and subcortical structures (Fig. 2). For cortical thickness, the genetic liability to sepsis-related death was associated with increased thickness of the precentral gyrus (β = 0.0039, SE = 0.0015, P = .0075) and transverse temporal gyrus (β = 0.0063, SE = 0.0026, P = .0144) as well as reduced thickness of lingual gyrus (β = −0.0028, SE = 0.0014, P = .0448) (Fig. 3). However, 3 MR analysis methods showed inconsistency in the direction of the magnitude of the effect of sepsis-related death in relation to the thickness of the precentral gyrus and lingual gyrus (Table S6, Supplemental Digital Content, http://links.lww.com/MD/N414). For cortical surface area, the genetic liability to sepsis-related death exhibited positive correlations with surface area of the fusiform gyrus (β = 11.0920, SE = 3.6412, P = .0023) and posterior cingulate gyrus (β = 3.6530, SE = 1.6684, P = .0286), While it exhibited negative correlations with surface area of the caudal middle frontal gyrus (β = −11.4586, SE = 5.1501, P = .0261) and frontal pole (β = −1.0024, SE = 0.4329, P = .0206) (Fig. 4). After Bonferroni correction, all of these results were nominally significant. Moreover, we conducted several sensitivity analyses, including Cochran Q test, MR-Egger regression, and MR-PRESSO global test, confirming the rigidity of the above results (Tables S3, S6, and S7, Supplemental Digital Content, http://links.lww.com/MD/N414).

3.3. Confounding analysis and Reverse MR analysis

Regarding postcentral cortex surface area, 1 SNP (rs4841254) in sepsis phenotype was associated with irritability-related phenotypes. Causality remained significant after removal of this SNP (βIVW = −34.0216, P = .02631). Regarding cortex surface area, SNPs of the sepsis-related death phenotype were discovered to be causally related to caudal middle frontal gyrus, frontal pole, fusiform gyrus, and posterior cingulate gyrus cortex surface area, independent of any confounding factors (Fig. 5). No confounders were also found for SNPs related to a causal relationship between genetic liability to sepsis-related death and transverse temporal gyrus cortical thickness. In addition, reverse MR showed a positive correlation between the thickness of the transverse temporal gyrus (β = 3.1983, SE = 1.4157, P = .0239) and the risk of sepsis-related death (Table S8, Supplemental Digital Content, http://links.lww.com/MD/N414).

Figure 5. The result of confounding analysis.

4. Discussion

To our knowledge, this is the first study using MR analysis to examine the causal effect of genetic liability to sepsis and sepsis-related death on specific brain region structural changes. Our results showed suggestive associations between genetic liability to sepsis and the postcentral gyrus surface area and sepsis-related death with regional cortical surface area (caudal middle frontal gyrus, frontal pole, fusiform gyrus, and posterior cingulate gyrus). We also found a possible bidirectional causal association between genetic liability to sepsis-related death and the thickness of the transverse temporal gyrus.

For cortical thickness, our study showed no evident associations between genetic liability to sepsis and cortical thickness, but a bidirectional causal association between genetic liability to sepsis-related death and thickness of transverse temporal gyrus. The transverse temporal gyrus contains the primary auditory cortex, which can be damaged by infection, trauma, ischemia, and tumors.[31,32] A cohort study by Kamkwalala et al showed elevated plasma soluble CD14 was associated with volume reductions in several regions of the temporal lobe in women with HIV.[33] This is inconsistent with our findings, suggesting that structural changes in the brain may be caused by pathological processes subsequent to sepsis, such as systemic inflammation, in addition to infection. Furthermore, brain anatomy is closely related to the brain’s ability to resist external pathogens, which may explain the poor prognosis of sepsis due to the increased cortical thickness of transverse temporal gyrus.[34] However, due to a paucity of studies explicitly linking sepsis to the transverse temporal gyrus, we are unable to corroborate the causal relationship discovered.

For cortical surface area, genetic liability to sepsis was suggested to be negatively correlated with surface area of the postcentral gyrus. Whereas the genetic liability to sepsis-related death was normally significantly negatively correlated with surface area of the caudal middle frontal gyrus and frontal pole, and normally significantly positively correlated with surface area of the fusiform gyrus and posterior cingulate gyrus. The postcentral gyrus is the location of the main somatosensory cortex, a significant brain region in charge of proprioception. We hypothesize that sensory dysfunction in patients with sepsis may be mediated by affecting the structure and function of the posterior central gyrus.[35] The caudal middle frontal gyrus and the frontal pole are both part of the frontal cortex. An LPS-induced sepsis rats model by Semmler et al found a reduced density of neurons in the prefrontal cortex during acute sepsis.[36] Additionally, sepsis typically affects prefrontal cortex-related cognitive abilities such as verbal fluency and memory.[37] The posterior cingulate gyrus is a component of the cingulate gyrus, which has been demonstrated to have considerable neuronal loss in SIBD patients.[38] As a crucial node in the salience network of the brain, cingulate gyrus failure can result in apathy and impaired consciousness that are common in sepsis.[38] And there has been little investigation into sepsis and the fusiform gyrus. Moreover, an interesting finding is that these brain structures have been demonstrated in previous studies to be related with Alzheimer disease. A MR study suggested that sepsis significantly increased the risk of Alzheimer disease.[39] We hypothesize that sepsis may influence the development of Alzheimer disease with mediators of specific brain structures.[40] In sum, research related to sepsis and these brain structures is lacking, but there is some indirect evidence of their relevance. Therefore, further research on the structural and functional alterations, as well as related mechanisms of the specific brain structures during sepsis is required.

For subcortical structures, neither genetic liability to sepsis-related death nor the sepsis-related death had a causal relationship with subcortical structure volume. However, contrary to our findings, most previous studies reported reduced volume of subcortical structures such as amygdala, caudate, hippocampus, putamen, and thalamus in sepsis.[8,38,41] The inconsistencies could be attributed to the insufficient sample size as well as the bias of confounding factors and reverse causal associations that are difficult to avoid in observational studies. Based on large sample size GWAS data and minimal confounding factor bias, our study can offer strong support for the inference of causation.

There are several possible mechanisms that might explain how sepsis affects brain structure. Firstly, sepsis may trigger microglia activation via various mechanisms, and notably, the sustained activation of these cells creates a cytotoxic environment, causing neuronal dysfunction and cell death.[42] Secondly, damage to the blood–brain barrier involving intercellular junctions of the endothelial cells may allow peripheral toxic substances to enter the central nervous system, leading to additional damage.[43] Thirdly, microcirculatory dysfunctions in sepsis could lead to insufficient perfusion, resulting in cerebral ischemic necrosis. Post‐mortem study of patients who died of septic shock frequently reveal cerebral ischemic lesions.[6,44] Lastly, the prevalent gut microbiota dysbiosis in sepsis might lead to brain injury via the microbiota-gut-brain axis.[45]

The principal strength of our study lies in its utilization of comprehensive MR analysis to elucidate the causal relationship between sepsis and brain structures, effectively circumventing the confounding bias and reverse causal associations inherent in traditional observational studies. Nevertheless, some limitations needed to be considered in the current study. Firstly, due to our reliance on summary-level statistics instead of individual-level data, the study was unable to conduct subgroup analyses. Secondly, the predominance of participants with European ancestry limits the generalizability of our findings to other ethnic groups. Finally, while our study proposes a causal link between sepsis and specific brain structures, it underscores the need for further research to explore potential functional alterations and the mechanisms underpinning these changes.

In summary, our study suggested a potential causal link between genetic liability to sepsis and cortical/subcortical structural changes, adding new evidence for a relationship between sepsis and brain structure changes.

Acknowledgments

We thank all the authors and participants of GWASs, the online scientific mapping site RAWGraphs 2.0 beta (https://app.rawgraphs.io).

Author contributions

Conceptualization: Dengfeng Zhou, Weina Wang, Qiaofa Lu.

Data curation: Dengfeng Zhou, Weina Wang.

Formal analysis: Dengfeng Zhou, Weina Wang.

Investigation: Dengfeng Zhou, Jiaying Gu.

Methodology: Dengfeng Zhou.

Validation: Weina Wang.

Visualization: Weina Wang.

Writing – original draft: Dengfeng Zhou, Weina Wang.

Writing – review & editing: Qiaofa Lu.

Supplementary Material

Abbreviations:

GWAS genome-wide association study

IVs instrumental variables

IVW inverse variance weighted

MR Mendelian randomization

SE standard error

SNPs single nucleotide polymorphisms

This study did not require separate ethical approval.

The authors have no funding and conflicts of interest to disclose.

The datasets generated during and/or analyzed during the current study are publicly available.

Supplemental Digital Content is available for this article.

How to cite this article: Zhou D, Wang W, Gu J, Lu Q. Causal effects of sepsis on structural changes in cerebral cortex: A Mendelian randomization investigation. Medicine 2024;103:36(e39404).
==== Refs
References

[1] Fleischmann C Scherag A Adhikari NK . International Forum of Acute Care Trialists. Assessment of global incidence and mortality of hospital-treated sepsis. Current estimates and limitations. Am J Respir Crit Care Med. 2016;193 :259–72.26414292
[2] Rudd KE Johnson SC Agesa KM . Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study. Lancet. 2020;395 :200–11.31954465
[3] Singer M Deutschman CS Seymour CW . The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315 :801–10.26903338
[4] Widmann CN Heneka MT . Long-term cerebral consequences of sepsis. Lancet Neurol. 2014;13 :630–6.24849863
[5] Sakusic A O’Horo JC Dziadzko M . Potentially modifiable risk factors for long-term cognitive impairment after critical illness: a systematic review. Mayo Clin Proc. 2018;93 :68–82.29304923
[6] Heming N Mazeraud A Verdonk F Bozza FA Chrétien F Sharshar T . Neuroanatomy of sepsis-associated encephalopathy. Crit Care. 2017;21 :65.28320461
[7] Wintermann GB Brunkhorst FM Petrowski K . Stress disorders following prolonged critical illness in survivors of severe sepsis. Crit Care Med. 2015;43 :1213–22.25760659
[8] Orhun G Tüzün E Bilgiç B . Brain volume changes in patients with acute brain dysfunction due to sepsis. Neurocrit Care. 2020;32 :459–68.31187433
[9] Hosokawa T Kinoshita K Ihara S . Relationship between brain volume reduction during the acute phase of sepsis and activities of daily living in elderly patients: a prospective cohort study. PLoS One. 2023;18 :e0284886.37192211
[10] Nakae R Sekine T Tagami T . Rapidly progressive brain atrophy in septic ICU patients: a retrospective descriptive study using semiautomatic CT volumetry. Crit Care. 2021;25 :411.34844648
[11] Gunther ML Morandi A Krauskopf E . VISIONS Investigation, VISualizing Icu SurvivOrs Neuroradiological Sequelae. The association between brain volumes, delirium duration, and cognitive outcomes in intensive care unit survivors: the VISIONS cohort magnetic resonance imaging study*. Crit Care Med. 2012;40 :2022–32.22710202
[12] Smith GD Ebrahim S . “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32 :1–22.12689998
[13] Davey Smith G Hemani G . Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23 :R89–98.25064373
[14] Williams JA Burgess S Suckling J . PIMS Collaboration. Inflammation and brain structure in schizophrenia and other neuropsychiatric disorders: a Mendelian randomization study. JAMA Psychiatry. 2022;79 :498–507.35353173
[15] Mavromatis LA Rosoff DB Cupertino RB Garavan H Mackey S Lohoff FW . Association between brain structure and alcohol use behaviors in adults: a Mendelian randomization and multiomics study. JAMA Psychiatry. 2022;79 :869–78.35947372
[16] Zeng Y Guo R Cao S Yang H . Causal associations between blood lipids and brain structures: a Mendelian randomization study. Cereb Cortex. 2023;33 :10901–8.37718242
[17] Skrivankova VW Richmond RC Woolf BAR . Strengthening the reporting of observational studies in epidemiology using Mendelian randomization: the STROBE-MR statement. JAMA. 2021;326 :1614–21.34698778
[18] Boef AG Dekkers OM le Cessie S . Mendelian randomization studies: a review of the approaches used and the quality of reporting. Int J Epidemiol. 2015;44 :496–511.25953784
[19] Bycroft C Freeman C Petkova D . The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562 :203–9.30305743
[20] Hamilton FW Thomas M Arnold D . Therapeutic potential of IL6R blockade for the treatment of sepsis and sepsis-related death: a Mendelian randomisation study. PLoS Med. 2023;20 :e1004174.36716318
[21] Zekavat SM Lin SH Bick AG . Biobank Japan Project. Hematopoietic mosaic chromosomal alterations increase the risk for diverse types of infection. Nat Med. 2021;27 :1012–24.34099924
[22] Grasby KL Jahanshad N Painter JN . Alzheimer’s Disease Neuroimaging Initiative. The genetic architecture of the human cerebral cortex. Science. 2020;367 :eaay6690.32193296
[23] Hibar DP Stein JL Renteria ME . Alzheimer’s Disease Neuroimaging Initiative. Common genetic variants influence human subcortical brain structures. Nature. 2015;520 :224–9.25607358
[24] Zhou S Wei T Liu X . Causal effects of COVID-19 on structural changes in specific brain regions: a Mendelian randomization study. BMC Med. 2023;21 :261.37468885
[25] Pierce BL Burgess S . Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. 2013;178 :1177–84.23863760
[26] Bowden J Davey Smith G Burgess S . Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44 :512–25.26050253
[27] Bowden J Davey Smith G Haycock PC Burgess S . Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40 :304–14.27061298
[28] Verbanck M Chen CY Neale B Do R . Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50 :693–8.29686387
[29] Mulugeta A Lumsden A Hyppönen E . Unlocking the causal link of metabolically different adiposity subtypes with brain volumes and the risks of dementia and stroke: a Mendelian randomization study. Neurobiol Aging. 2021;102 :161–9.33770530
[30] Song W Qian W Wang W Yu S Lin GN . Mendelian randomization studies of brain MRI yield insights into the pathogenesis of neuropsychiatric disorders. BMC Genomics. 2021;22 (Suppl 3 ):342.34078268
[31] Qi R Cao Z Surento W . RORA rs8042149 polymorphism moderates the association between PTSD symptom severity and transverse temporal gyrus thickness in Han Chinese adults who lost their only child. J Affect Disord. 2022;314 :318–24.35878841
[32] Mirsattari SM Blume WT . Update on temporal lobe epilepsy. Epilepsy Res Treat. 2013;2013 :579487.23533737
[33] Kamkwalala AR Wang X Maki PM . Brief report: higher peripheral monocyte activation markers are associated with smaller frontal and temporal cortical volumes in women with HIV. J Acquir Immune Defic Syndr. 2020;84 :54–9.31914004
[34] Alves de Lima K Rustenhoven J Kipnis J . Meningeal immunity and its function in maintenance of the central nervous system in health and disease. Annu Rev Immunol. 2020;38 :597–620.32340575
[35] Zauner C Gendo A Kramer L . Impaired subcortical and cortical sensory evoked potential pathways in septic patients. Crit Care Med. 2002;30 :1136–9.12006815
[36] Semmler A Frisch C Debeir T . Long-term cognitive impairment, neuronal loss and reduced cortical cholinergic innervation after recovery from sepsis in a rodent model. Exp Neurol. 2007;204 :733–40.17306796
[37] Iwashyna TJ Ely EW Smith DM Langa KM . Long-term cognitive impairment and functional disability among survivors of severe sepsis. JAMA. 2010;304 :1787–94.20978258
[38] Orhun G Esen F Özcan PE . Neuroimaging findings in sepsis-induced brain dysfunction: association with clinical and laboratory findings. Neurocrit Care. 2019;30 :106–17.30027347
[39] Yang H Xu H Li Q . Study of brain morphology change in Alzheimer’s disease and amnestic mild cognitive impairment compared with normal controls. Gen Psychiatr. 2019;32 :e100005.31179429
[40] Zeng Y Cao S Pang K Tang J Lin G . Causal association between sepsis and neurodegenerative diseases: a bidirectional two-sample Mendelian randomization study. J Alzheimers Dis. 2024;97 :229–37.38189756
[41] Semmler A Widmann CN Okulla T . Persistent cognitive impairment, hippocampal atrophy and EEG changes in sepsis survivors. J Neurol Neurosurg Psychiatry. 2013;84 :62–9.23134661
[42] Zrzavy T Höftberger R Berger T . Pro-inflammatory activation of microglia in the brain of patients with sepsis. Neuropathol Appl Neurobiol. 2019;45 :278–90.29804289
[43] Erikson K Tuominen H Vakkala M . Brain tight junction protein expression in sepsis in an autopsy series. Crit Care. 2020;24 :385.32600371
[44] Sharshar T Annane D de la Grandmaison GL Brouland JP Hopkinson NS Françoise G . The neuropathology of septic shock. Brain Pathol. 2004;14 :21–33.14997934
[45] Zhang H Xu J Wu Q . Gut microbiota mediates the susceptibility of mice to sepsis-associated encephalopathy by butyric acid. J Inflamm Res. 2022;15 :2103–19.35386224
