A Comparison of Methods to Harmonize Cortical Thickness Measurements Across Scanners and Sites

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants’ demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LME INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME INT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2–81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3–85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ 2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ 2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ 2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects

Remodeling of the Cortical Structural Connectome in Posttraumatic Stress Disorder:Results from the ENIGMA-PGC PTSD Consortium

BACKGROUND: Posttraumatic stress disorder (PTSD) is accompanied by disrupted cortical neuroanatomy. We investigated alteration in covariance of structural networks associated with PTSD in regions that demonstrate the case-control differences in cortical thickness (CT) and surface area (SA). METHODS: Neuroimaging and clinical data were aggregated from 29 research sites in >1,300 PTSD cases and >2,000 trauma-exposed controls (age 6.2-85.2 years) by the ENIGMA-PGC PTSD working group. Cortical regions in the network were rank-ordered by effect size of PTSD-related cortical differences in CT and SA. The top-n (n = 2 to 148) regions with the largest effect size for PTSD > non-PTSD formed hypertrophic networks, the largest effect size for PTSD < non-PTSD formed atrophic networks, and the smallest effect size of between-group differences formed stable networks. The mean structural covariance (SC) of a given n-region network was the average of all positive pairwise correlations and was compared to the mean SC of 5,000 randomly generated n-region networks. RESULTS: Patients with PTSD, relative to non-PTSD controls, exhibited lower mean SC in CT-based and SA-based atrophic networks. Comorbid depression, sex and age modulated covariance differences of PTSD-related structural networks. CONCLUSIONS: Covariance of structural networks based on CT and cortical SA are affected by PTSD and further modulated by comorbid depression, sex, and age. The structural covariance networks that are perturbed in PTSD comport with converging evidence from resting state functional connectivity networks and networks impacted by inflammatory processes, and stress hormones in PTSD

LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories

We present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples.Comment: 28 pages, 10 figures, 3 table

Measurement of the high-energy all-flavor neutrino-nucleon cross section with IceCube

The flux of high-energy neutrinos passing through the Earth is attenuated due to their interactions with matter. The interaction rate is determined by the neutrino interaction cross section and affects the flux arriving at the IceCube Neutrino Observatory, a cubic-kilometer neutrino detector embedded in the Antarctic ice sheet. We present a measurement of the neutrino cross section between 60 TeV and 10 PeV using the high-energy starting event (HESE) sample from IceCube with 7.5 years of data. The result is binned in neutrino energy and obtained using both Bayesian and frequentist statistics. We find it compatible with predictions from the Standard Model. While the cross section is expected to be flavor independent above 1 TeV, additional constraints on the measurement are included through updated experimental particle identification (PID) classifiers, proxies for the three neutrino flavors. This is the first such measurement to use a ternary PID observable and the first to account for neutrinos from tau decay

Searching for eV-scale sterile neutrinos with eight years of atmospheric neutrinos at the IceCube neutrino telescope

We report in detail on searches for eV-scale sterile neutrinos, in the context of a 3+1 model, using eight years of data from the IceCube neutrino telescope. By analyzing the reconstructed energies and zenith angles of 305,735 atmospheric $\nu_\mu$ and $\bar{\nu}_\mu$ events we construct confidence intervals in two analysis spaces: $\sin^2 (2\theta_{24})$ vs. $\Delta m^2_{41}$ under the conservative assumption $\theta_{34}=0$; and $\sin^2(2\theta_{24})$ vs. $\sin^2 (2\theta_{34})$ given sufficiently large $\Delta m^2_{41}$ that fast oscillation features are unresolvable. Detailed discussions of the event selection, systematic uncertainties, and fitting procedures are presented. No strong evidence for sterile neutrinos is found, and the best-fit likelihood is consistent with the no sterile neutrino hypothesis with a p-value of 8\% in the first analysis space and 19\% in the second.Comment: This long-form paper is a companion to the letter "An eV-scale sterile neutrino search using eight years of atmospheric muon neutrino data from the IceCube Neutrino Observatory". v2: update other experiments contours on results plo

An eV-scale sterile neutrino search using eight years of atmospheric muon neutrino data from the IceCube Neutrino Observatory

The results of a 3+1 sterile neutrino search using eight years of data from the IceCube Neutrino Observatory are presented. A total of 305,735 muon neutrino events are analyzed in reconstructed energy-zenith space to test for signatures of a matter-enhanced oscillation that would occur given a sterile neutrino state with a mass-squared differences between 0.01\,eV$^2$ and 100\,eV$^2$. The best-fit point is found to be at $\sin^2(2\theta_{24})=0.10$ and $\Delta m_{41}^2 = 4.5{\rm eV}^2$, which is consistent with the no sterile neutrino hypothesis with a p-value of 8.0\%.Comment: 11 pages, 5 figures. This letter is supported by the long-form paper "Searching for eV-scale sterile neutrinos with eight years of atmospheric neutrinos at the IceCube neutrino telescope," also appearing on arXiv. Digital data release available at: https://github.com/icecube/HE-Sterile-8year-data-releas

Paper Trails: Following the Money

In many recent elections, the candidates who raise the most money have a better shot at winning, so candidates must raise millions of dollars to win an election. A top question to consider in all elections: Where is the money coming from? Posting about the financing behind federal elections from In All Things - an online hub committed to the claim that the life, death, and resurrection of Jesus Christ has implications for the entire world. http://inallthings.org/paper-trails-following-the-money

Human subcortical brain asymmetries in 15,847 people worldwide reveal effects of age and sex

The two hemispheres of the human brain differ functionally and structurally. Despite over a century of research, the extent to which brain asymmetry is influenced by sex, handedness, age, and genetic factors is still controversial. Here we present the largest ever analysis of subcortical brain asymmetries, in a harmonized multi-site study using meta-analysis methods. Volumetric asymmetry of seven subcortical structures was assessed in 15,847 MRI scans from 52 datasets worldwide. There were sex differences in the asymmetry of the globus pallidus and putamen. Heritability estimates, derived from 1170 subjects belonging to 71 extended pedigrees, revealed that additive genetic factors influenced the asymmetry of these two structures and that of the hippocampus and thalamus. Handedness had no detectable effect on subcortical asymmetries, even in this unprecedented sample size, but the asymmetry of the putamen varied with age. Genetic drivers of asymmetry in the hippocampus, thalamus and basal ganglia may affect variability in human cognition, including susceptibility to psychiatric disorders

Cortical brain abnormalities in 4474 individuals with schizophrenia and 5098 control subjects via the enhancing neuro Imaging genetics through meta analysis (ENIGMA) Consortium

BACKGROUND: The profile of cortical neuroanatomical abnormalities in schizophrenia is not fully understood, despite hundreds of published structural brain imaging studies. This study presents the first meta-analysis of cortical thickness and surface area abnormalities in schizophrenia conducted by the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) Schizophrenia Working Group. METHODS: The study included data from 4474 individuals with schizophrenia (mean age, 32.3 years; range, 11-78 years; 66% male) and 5098 healthy volunteers (mean age, 32.8 years; range, 10-87 years; 53% male) assessed with standardized methods at 39 centers worldwide. RESULTS: Compared with healthy volunteers, individuals with schizophrenia have widespread thinner cortex (left/right hemisphere: Cohen's d = -0.530/-0.516) and smaller surface area (left/right hemisphere: Cohen's d = -0.251/-0.254), with the largest effect sizes for both in frontal and temporal lobe regions. Regional group differences in cortical thickness remained significant when statistically controlling for global cortical thickness, suggesting regional specificity. In contrast, effects for cortical surface area appear global. Case-control, negative, cortical thickness effect sizes were two to three times larger in individuals receiving antipsychotic medication relative to unmedicated individuals. Negative correlations between age and bilateral temporal pole thickness were stronger in individuals with schizophrenia than in healthy volunteers. Regional cortical thickness showed significant negative correlations with normalized medication dose, symptom severity, and duration of illness and positive correlations with age at onset. CONCLUSIONS: The findings indicate that the ENIGMA meta-analysis approach can achieve robust findings in clinical neuroscience studies; also, medication effects should be taken into account in future genetic association studies of cortical thickness in schizophrenia