Search

The linear mixed-effects model (LME) is a versatile approach to account for dependence among observations. Many large-scale neuroimaging datasets with complex designs have increased the need for LME; however LME has seldom been used in whole-brain imaging analyses due to its heavy computational requirements. In this paper, we introduce a fast and efficient mixed-effects algorithm (FEMA) that makes whole-brain vertex-wise, voxel-wise, and connectome-wide LME analyses in large samples possible. We validate FEMA with extensive simulations, showing that the estimates of the fixed effects are equivalent to standard maximum likelihood estimates but obtained with orders of magnitude improvement in computational speed. We demonstrate the applicability of FEMA by studying the cross-sectional and longitudinal effects of age on region-of-interest level and vertex-wise cortical thickness, as well as connectome-wide functional connectivity values derived from resting state functional MRI, using longitudinal imaging data from the Adolescent Brain Cognitive DevelopmentSM Study release 4.0. Our analyses reveal distinct spatial patterns for the annualized changes in vertex-wise cortical thickness and connectome-wide connectivity values in early adolescence, highlighting a critical time of brain maturation. The simulations and application to real data show that FEMA enables advanced investigation of the relationships between large numbers of neuroimaging metrics and variables of interest while considering complex study designs, including repeated measures and family structures, in a fast and efficient manner. The source code for FEMA is available via: https://github.com/cmig-research-group/cmig_tools/.

BackgroundChildhood disadvantage is a prominent risk factor for cognitive and brain aging. Childhood disadvantage is associated with poorer episodic memory in late midlife and functional and structural brain abnormalities in the default mode network (DMN). Although age-related changes in DMN are associated with episodic memory declines in older adults, it remains unclear if childhood disadvantage has an enduring impact on this later-life brain-cognition relationship earlier in the aging process. Here, within the DMN, we examined whether its cortical microstructural integrity-an early marker of structural vulnerability that increases the risk for future cognitive decline and neurodegeneration-is associated with episodic memory in adults at ages 56-66, and whether childhood disadvantage moderates this association.MethodsCortical mean diffusivity (MD) obtained from diffusion magnetic resonance imaging was used to measure microstructural integrity in 350 community-dwelling men. We examined both visual and verbal episodic memory in relation to DMN MD and divided participants into disadvantaged and nondisadvantaged groups based on parental education and occupation.ResultsHigher DMN MD was associated with poorer visual memory but not verbal memory (β = -0.11, p = .040 vs β = -0.04, p = .535). This association was moderated by childhood disadvantage and was significant only in the disadvantaged group (β = -0.26, p = .002 vs β = -0.00, p = .957).ConclusionsLower DMN cortical microstructural integrity may reflect visual memory vulnerability in cognitively normal adults earlier in the aging process. Individuals who experienced childhood disadvantage manifested greater vulnerability to cortical microstructure-related visual memory dysfunction than their nondisadvantaged counterparts who exhibited resilience in the face of low cortical microstructural integrity.

BackgroundPlasma neurofilament light chain (NfL) is a promising biomarker of neurodegeneration with potential clinical utility in monitoring the progression of neurodegenerative diseases. However, the cross-sectional associations of plasma NfL with measures of cognition and brain have been inconsistent in community-dwelling populations.MethodsWe examined these associations in a large community-dwelling sample of early old age men (N = 969, mean age = 67.57 years, range = 61-73 years), who are either cognitively unimpaired (CU) or with mild cognitive impairment (MCI). Specifically, we investigated five cognitive domains (executive function, episodic memory, verbal fluency, processing speed, visual-spatial ability), as well as neuroimaging measures of gray and white matter.ResultsAfter adjusting for age, health status, and young adult general cognitive ability, plasma NfL level was only significantly associated with processing speed and white matter hyperintensity (WMH) volume, but not with other cognitive or neuroimaging measures. The association with processing speed was driven by individuals with MCI, as it was not detected in CU individuals.ConclusionsThese results suggest that in early old age men without dementia, plasma NfL does not appear to be sensitive to cross-sectional individual differences in most domains of cognition or neuroimaging measures of gray and white matter. The revealed plasma NfL associations were limited to WMH for all participants and processing speed only within the MCI cohort. Importantly, considering cognitive status in community-based samples will better inform the interpretation of the relationships of plasma NfL with cognition and brain and may help resolve mixed findings in the literature.

BackgroundHigh b-value diffusion-weighted images (DWI) are used for detection of clinically significant prostate cancer (csPCa). This study qualitatively and quantitatively compares synthesized DWI (sDWI) to acquired (aDWI) for detection of csPCa.MethodsOne hundred fifty-one consecutive patients who underwent prostate MRI and biopsy were included in the study. Axial DWI with b = 0, 500, 1000, and 2000 s/mm2 using a 3T clinical scanner using a 32-channel phased-array body coil were acquired. We retrospectively synthesized DWI for b = 2000 s/mm2 via extrapolation based on mono-exponential decay, using b = 0 and b = 500 s/mm2 (sDWI500) and b = 0, b = 500 s/mm2, and b = 1000 s/mm2 (sDWI1000). Differences in signal intensity between sDWI and aDWI were evaluated within different regions of interest (prostate alone, prostate plus 5 mm, 30 mm and 70 mm margin and full field of view). The maximum DWI value within each ROI was evaluated for prediction of csPCa. Classification accuracy was compared to Restriction Spectrum Imaging restriction score (RSIrs), a previously validated biomarker based on multi-exponential DWI. Discrimination of csPCa was evaluated via area under the receiver operating characteristic curve (AUC).ResultsWithin the prostate, mean ± standard deviation of percent mean differences between sDWI and aDWI signal were -46 ± 35% for sDWI1000 and -67 ± 24% for sDWI500. AUC for aDWI, sDWI500, sDWI1000, and RSIrs within the prostate 0.62[95% confidence interval: 0.53, 0.71], 0.63[0.54, 0.72], 0.65[0.56, 0.73] and 0.78[0.71, 0.86], respectively.ConclusionsDWI is qualitatively comparable to aDWI within the prostate. However, hyperintense artifacts are introduced with sDWI in the surrounding pelvic tissue that interfere with quantitative cancer detection and might mask metastases. In the prostate, RSIrs yields superior quantitative csPCa detection than sDWI or aDWI.

IntroductionThere is a pressing need for non-invasive, cost-effective tools for early detection of Alzheimer's disease (AD).MethodsUsing data from the Alzheimer's Disease Neuroimaging Initiative (ADNI), Cox proportional models were conducted to develop a multimodal hazard score (MHS) combining age, a polygenic hazard score (PHS), brain atrophy, and memory to predict conversion from mild cognitive impairment (MCI) to dementia. Power calculations estimated required clinical trial sample sizes after hypothetical enrichment using the MHS. Cox regression determined predicted age of onset for AD pathology from the PHS.ResultsThe MHS predicted conversion from MCI to dementia (hazard ratio for 80th versus 20th percentile: 27.03). Models suggest that application of the MHS could reduce clinical trial sample sizes by 67%. The PHS alone predicted age of onset of amyloid and tau.DiscussionThe MHS may improve early detection of AD for use in memory clinics or for clinical trial enrichment.HighlightsA multimodal hazard score (MHS) combined age, genetics, brain atrophy, and memory. The MHS predicted time to conversion from mild cognitive impairment to dementia. MHS reduced hypothetical Alzheimer's disease (AD) clinical trial sample sizes by 67%. A polygenic hazard score predicted age of onset of AD neuropathology.

Some evidence suggests a biphasic pattern of changes in cortical thickness wherein higher, rather than lower, thickness is associated with very early Alzheimer's disease (AD) pathology. We examined whether integrating information from AD brain signatures based on mean diffusivity (MD) can aid in the interpretation of cortical thickness/volume as a risk factor for future AD-related changes. Participants were 572 men in the Vietnam Era Twin Study of Aging who were cognitively unimpaired at baseline (mean age = 56 years; range = 51-60). Individuals with both high thickness/volume signatures and high MD signatures at baseline had lower cortical thickness/volume in AD signature regions and lower episodic memory performance 12 years later compared to those with high thickness/volume and low MD signatures at baseline. Groups did not differ in level of young adult cognitive reserve. Our findings are in line with a biphasic model in which increased cortical thickness may precede future decline and establish the value of examining cortical MD alongside cortical thickness to identify subgroups with differential risk for poorer brain and cognitive outcomes.

BackgroundComposite scores of magnetic resonance imaging-derived metrics in brain regions associated with Alzheimer's disease (AD), commonly termed AD signatures, have been developed to distinguish early AD-related atrophy from normal age-associated changes. Diffusion-based gray matter signatures may be more sensitive to early AD-related changes compared with thickness/volume-based signatures, demonstrating their potential clinical utility. The timing of early (i.e., midlife) changes in AD signatures from different modalities and whether diffusion- and thickness/volume-based signatures each capture unique AD-related phenotypic or genetic information remains unknown.MethodsOur validated thickness/volume signature, our novel mean diffusivity (MD) signature, and a magnetic resonance imaging-derived measure of brain age were used in biometrical analyses to examine genetic and environmental influences on the measures as well as phenotypic and genetic relationships between measures over 12 years. Participants were 736 men from 3 waves of the Vietnam Era Twin Study of Aging (VETSA) (baseline/wave 1: mean age [years] = 56.1, SD = 2.6, range = 51.1-60.2). Subsequent waves occurred at approximately 5.7-year intervals.ResultsMD and thickness/volume signatures were highly heritable (56%-72%). Baseline MD signatures predicted thickness/volume signatures over a decade later, but baseline thickness/volume signatures showed a significantly weaker relationship with future MD signatures. AD signatures and brain age were correlated, but each measure captured unique phenotypic and genetic variance.ConclusionsCortical MD and thickness/volume AD signatures are heritable, and each signature captures unique variance that is also not explained by brain age. Moreover, results are in line with changes in MD emerging before changes in cortical thickness, underscoring the utility of MD as a very early predictor of AD risk.

Findings in adults have shown that crystallized measures of intelligence, which are more culturally sensitive than fluid intelligence measures, have greater heritability; however, these results have not been found in children. The present study used data from 8,518 participants between 9 and 11 years old from the Adolescent Brain Cognitive Development (ABCD) Study. We found that polygenic predictors of intelligence test performance (based on genome-wide association meta-analyses of data from 269,867 individuals) and of educational attainment (based on data from 1.1 million individuals) predicted neurocognitive performance. We found that crystallized measures were more strongly associated with both polygenic predictors than were fluid measures. This mirrored heritability differences reported previously in adults and suggests similar associations in children. This may be consistent with a prominent role of gene-environment correlation in cognitive development measured by crystallized intelligence tests. Environmental and experiential mediators may represent malleable targets for improving cognitive outcomes.

Twin and family studies have historically aimed to partition phenotypic variance into components corresponding to additive genetic effects (A), common environment (C), and unique environment (E). Here we present the ACE Model and several extensions in the Adolescent Brain Cognitive Development℠ Study (ABCD Study®), employed using the new Fast Efficient Mixed Effects Analysis (FEMA) package. In the twin sub-sample (n = 924; 462 twin pairs), heritability estimates were similar to those reported by prior studies for height (twin heritability = 0.86) and cognition (twin heritability between 0.00 and 0.61), respectively. Incorporating SNP-derived genetic relatedness and using the full ABCD Study® sample (n = 9,742) led to narrower confidence intervals for all parameter estimates. By leveraging the sparse clustering method used by FEMA to handle genetic relatedness only for participants within families, we were able to take advantage of the diverse distribution of genetic relatedness within the ABCD Study® sample.

Characterizing the optimal fMRI paradigms for detecting behaviorally relevant functional connectivity (FC) patterns is a critical step to furthering our knowledge of the neural basis of behavior. Previous studies suggested that FC patterns derived from task fMRI paradigms, which we refer to as task-based FC, are better correlated with individual differences in behavior than resting-state FC, but the consistency and generalizability of this advantage across task conditions was not fully explored. Using data from resting-state fMRI and three fMRI tasks from the Adolescent Brain Cognitive Development Study ® (ABCD), we tested whether the observed improvement in behavioral prediction power of task-based FC can be attributed to changes in brain activity induced by the task design. We decomposed the task fMRI time course of each task into the task model fit (the fitted time course of the task condition regressors from the single-subject general linear model) and the task model residuals, calculated their respective FC, and compared the behavioral prediction performance of these FC estimates to resting-state FC and the original task-based FC. The FC of the task model fit was better than the FC of the task model residual and resting-state FC at predicting a measure of general cognitive ability or two measures of performance on the fMRI tasks. The superior behavioral prediction performance of the FC of the task model fit was content-specific insofar as it was only observed for fMRI tasks that probed similar cognitive constructs to the predicted behavior of interest. To our surprise, the task model parameters, the beta estimates of the task condition regressors, were equally if not more predictive of behavioral differences than all FC measures. These results showed that the observed improvement of behavioral prediction afforded by task-based FC was largely driven by the FC patterns associated with the task design. Together with previous studies, our findings highlighted the importance of task design in eliciting behaviorally meaningful brain activation and FC patterns.

Executive function encompasses effortful cognitive processes that are particularly susceptible to aging. Functional brain networks supporting executive function-such as the frontoparietal control network and the multiple demand system-have been extensively investigated. However, it remains unclear how structural networks facilitate and constrain the dynamics of functional networks to contribute to aging-related executive function declines. We examined whether changes in structural network modal controllability-a network's ability to facilitate effortful brain state transitions that support cognitive functions-are associated with changes in executive function cross-sectionally and longitudinally. Diffusion-weighted imaging and neuropsychological testing were conducted at two time points (Time 1: ages 56 to 66, N = 172; Time 2: ages 61 to 70, N = 267) in community-dwelling men from the Vietnam Era Twin Study of Aging. An executive function factor score was computed from six neuropsychological tasks. Structural networks constructed from white matter connectivity were used to estimate modal controllability in control network and multiple demand system. We showed that higher modal controllability in control network and multiple demand system was associated with better executive function at Time 2, after controlling for age, young adult general cognitive ability, and physical health status. Moreover, changes in executive function over a period of 5 to 6 years (Time 1-Time 2, N = 105) were associated with changes in modal controllability of the multiple demand system and weakly in the control network over the same time period. These findings suggest that changes in the ability of structural brain networks in facilitating effortful brain state transitions may be a key neural mechanism underlying aging-related executive function declines and cognitive aging.

PurposeDynamic contrast-enhanced MRI (DCE) and apparent diffusion coefficient (ADC) are currently used to evaluate treatment response of breast cancer. The purpose of the current study was to evaluate the three-component Restriction Spectrum Imaging model (RSI3C), a recent diffusion-weighted MRI (DWI)-based tumor classification method, combined with elastic image registration, to automatically monitor breast tumor size throughout neoadjuvant therapy.Experimental designBreast cancer patients (n=27) underwent multi-parametric 3T MRI at four time points during treatment. Elastically-registered DWI images were used to generate an automatic RSI3C response classifier, assessed against manual DCE tumor size measurements and mean ADC values. Predictions of therapy response during treatment and residual tumor post-treatment were assessed using non-pathological complete response (non-pCR) as an endpoint.ResultsTen patients experienced pCR. Prediction of non-pCR using ROC AUC (95% CI) for change in measured tumor size from pre-treatment time point to early-treatment time point was 0.65 (0.38-0.92) for the RSI3C classifier, 0.64 (0.36-0.91) for DCE, and 0.45 (0.16-0.75) for change in mean ADC. Sensitivity for detection of residual disease post-treatment was 0.71 (0.44-0.90) for the RSI3C classifier, compared to 0.88 (0.64-0.99) for DCE and 0.76 (0.50-0.93) for ADC. Specificity was 0.90 (0.56-1.00) for the RSI3C classifier, 0.70 (0.35-0.93) for DCE, and 0.50 (0.19-0.81) for ADC.ConclusionThe automatic RSI3C classifier with elastic image registration suggested prediction of response to treatment after only three weeks, and showed performance comparable to DCE for assessment of residual tumor post-therapy. RSI3C may guide clinical decision-making and enable tailored treatment regimens and cost-efficient evaluation of neoadjuvant therapy of breast cancer.

BackgroundMultiparametric magnetic resonance imaging (mpMRI) improves detection of clinically significant prostate cancer (csPCa), but the subjective Prostate Imaging Reporting and Data System (PI-RADS) system and quantitative apparent diffusion coefficient (ADC) are inconsistent. Restriction spectrum imaging (RSI) is an advanced diffusion-weighted MRI technique that yields a quantitative imaging biomarker for csPCa called the RSI restriction score (RSIrs).ObjectiveTo evaluate RSIrs for automated patient-level detection of csPCa.Design setting and participantsWe retrospectively studied all patients (n = 151) who underwent 3 T mpMRI and RSI (a 2-min sequence on a clinical scanner) for suspected prostate cancer at University of California San Diego during 2017-2019 and had prostate biopsy within 180 d of MRI.InterventionWe calculated the maximum RSIrs and minimum ADC within the prostate, and obtained PI-RADS v2.1 from medical records.Outcome measurements and statistical analysisWe compared the performance of RSIrs, ADC, and PI-RADS for the detection of csPCa (grade group ≥2) on the best available histopathology (biopsy or prostatectomy) using the area under the curve (AUC) with two-tailed α = 0.05. We also explored whether the combination of PI-RADS and RSIrs might be superior to PI-RADS alone and performed subset analyses within the peripheral and transition zones.Results and limitationsAUC values for ADC, RSIrs, and PI-RADS were 0.48 (95% confidence interval: 0.39, 0.58), 0.78 (0.70, 0.85), and 0.77 (0.70, 0.84), respectively. RSIrs and PI-RADS were each superior to ADC for patient-level detection of csPCa (p

Purpose To develop a multicompartmental signal model for whole-body diffusion-weighted imaging (DWI) and apply it to study the diffusion properties of normal tissue and metastatic prostate cancer bone lesions in vivo. Materials and Methods This prospective study (ClinicalTrials.gov: NCT03440554) included 139 men with prostate cancer (mean age, 70 years ± 9 [SD]). Multicompartmental models with two to four tissue compartments were fit to DWI data from whole-body scans to determine optimal compartmental diffusion coefficients. Bayesian information criterion (BIC) and model-fitting residuals were calculated to quantify model complexity and goodness of fit. Diffusion coefficients for the optimal model (having lowest BIC) were used to compute compartmental signal-contribution maps. The signal intensity ratio (SIR) of bone lesions to normal-appearing bone was measured on these signal-contribution maps and on conventional DWI scans and compared using paired t tests (α = .05). Two-sample t tests (α = .05) were used to compare compartmental signal fractions between lesions and normal-appearing bone. Results Lowest BIC was observed from the four-compartment model, with optimal compartmental diffusion coefficients of 0, 1.1 × 10-3, 2.8 × 10-3, and >3.0 ×10-2 mm2/sec. Fitting residuals from this model were significantly lower than from conventional apparent diffusion coefficient mapping (P < .001). Bone lesion SIR was significantly higher on signal-contribution maps of model compartments 1 and 2 than on conventional DWI scans (P < .008). The fraction of signal from compartments 2, 3, and 4 was also significantly different between metastatic bone lesions and normal-appearing bone tissue (P ≤ .02). Conclusion The four-compartment model best described whole-body diffusion properties. Compartmental signal contributions from this model can be used to examine prostate cancer bone involvement. Keywords: Whole-Body MRI, Diffusion-weighted Imaging, Restriction Spectrum Imaging, Diffusion Signal Model, Bone Metastases, Prostate Cancer Clinical trial registration no. NCT03440554 Supplemental material is available for this article. © RSNA, 2023 See also commentary by Margolis in this issue.

BackgroundEarly detection is critical for easing the rising burden of psychiatric disorders. However, the specificity of psychopathological measurements and genetic predictors is unclear among youth.MethodsWe measured associations between genetic risk for psychopathology (polygenic risk scores (PRS) and family history (FH) measures) and a wide range of behavioral measures in a large sample (n = 5,204) of early adolescent participants (9-11 years) from the Adolescent Brain and Cognitive Development StudySM . Associations were measured both with and without accounting for shared variance across measures of genetic risk.ResultsWhen controlling for genetic risk for other psychiatric disorders, polygenic risk for problematic opioid use (POU) is uniquely associated with lower behavioral inhibition. Attention deficit hyperactivity disorder (ADHD), depression (DEP), and attempted suicide (SUIC) PRS shared many significant associations with externalizing, internalizing, and psychosis-related behaviors. However, when accounting for all measures of genetic and familial risk, these PRS also showed clear, unique patterns of association. Polygenic risk for ASD, BIP, and SCZ, and attempted suicide uniquely predicted variability in cognitive performance. FH accounted for unique variability in behavior above and beyond PRS and vice versa, with FH measures explaining a greater proportion of unique variability compared to the PRS.ConclusionOur results indicate that, among youth, many behaviors show shared genetic influences; however, there is also specificity in the profile of emerging psychopathologies for individuals with high genetic risk for particular disorders. This may be useful for quantifying early, differential risk for psychopathology in development.

BackgroundPosttraumatic stress disorder (PTSD) and major depressive disorder (MDD) are commonly reported co-occurring mental health consequences of psychological trauma exposure. The disorders have high genetic overlap. Trauma is a complex phenotype but research suggests that trauma sensitivity has a heritable basis. We investigated whether sensitivity to trauma in those with MDD reflects a similar genetic component in those with PTSD.MethodsGenetic correlations between PTSD and MDD in individuals reporting trauma and MDD in individuals not reporting trauma were estimated, as well as with recurrent MDD and single-episode MDD, using genome-wide association study (GWAS) summary statistics. Genetic correlations were replicated using PTSD data from the Psychiatric Genomics Consortium and the Million Veteran Program. Polygenic risk scores were generated in UK Biobank participants who met the criteria for lifetime MDD (N = 29 471). We investigated whether genetic loading for PTSD was associated with reporting trauma in these individuals.ResultsGenetic loading for PTSD was significantly associated with reporting trauma in individuals with MDD [OR 1.04 (95% CI 1.01-1.07), Empirical-p = 0.02]. PTSD was significantly more genetically correlated with recurrent MDD than with MDD in individuals not reporting trauma (rg differences = ~0.2, p < 0.008). Participants who had experienced recurrent MDD reported significantly higher rates of trauma than participants who had experienced single-episode MDD (χ2 > 166, p < 0.001).ConclusionsOur findings point towards the existence of genetic variants associated with trauma sensitivity that might be shared between PTSD and MDD, although replication with better powered GWAS is needed. Our findings corroborate previous research highlighting trauma exposure as a key risk factor for recurrent MDD.

BackgroundProstate cancer risk stratification using single-nucleotide polymorphisms (SNPs) demonstrates considerable promise in men of European, Asian, and African genetic ancestries, but there is still need for increased accuracy. We evaluated whether including additional SNPs in a prostate cancer polygenic hazard score (PHS) would improve associations with clinically significant prostate cancer in multi-ancestry datasets.MethodsIn total, 299 SNPs previously associated with prostate cancer were evaluated for inclusion in a new PHS, using a LASSO-regularized Cox proportional hazards model in a training dataset of 72,181 men from the PRACTICAL Consortium. The PHS model was evaluated in four testing datasets: African ancestry, Asian ancestry, and two of European Ancestry-the Cohort of Swedish Men (COSM) and the ProtecT study. Hazard ratios (HRs) were estimated to compare men with high versus low PHS for association with clinically significant, with any, and with fatal prostate cancer. The impact of genetic risk stratification on the positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was also measured.ResultsThe final model (PHS290) had 290 SNPs with non-zero coefficients. Comparing, for example, the highest and lowest quintiles of PHS290, the hazard ratios (HRs) for clinically significant prostate cancer were 13.73 [95% CI: 12.43-15.16] in ProtecT, 7.07 [6.58-7.60] in African ancestry, 10.31 [9.58-11.11] in Asian ancestry, and 11.18 [10.34-12.09] in COSM. Similar results were seen for association with any and fatal prostate cancer. Without PHS stratification, the PPV of PSA testing for clinically significant prostate cancer in ProtecT was 0.12 (0.11-0.14). For the top 20% and top 5% of PHS290, the PPV of PSA testing was 0.19 (0.15-0.22) and 0.26 (0.19-0.33), respectively.ConclusionsWe demonstrate better genetic risk stratification for clinically significant prostate cancer than prior versions of PHS in multi-ancestry datasets. This is promising for implementing precision-medicine approaches to prostate cancer screening decisions in diverse populations.

Genome-Wide Association studies have typically been limited to univariate analysis in which a single outcome measure is tested against millions of variants. Recent work demonstrates that a Multivariate Omnibus Statistic Test (MOSTest) is well powered to discover genomic effects distributed across multiple phenotypes. Applied to cortical brain MRI morphology measures, MOSTest has resulted in a drastic improvement in power to discover loci when compared to established approaches (min-P). One question that arises is how well these discovered loci replicate in independent data. Here we perform 10 times cross validation within 34,973 individuals from UK Biobank for imaging measures of cortical area, thickness and sulcal depth (>1,000 dimensionality for each). By deploying a replication method that aggregates discovered effects distributed across multiple phenotypes, termed PolyVertex Score (MOSTest-PVS), we demonstrate a higher replication yield and comparable replication rate of discovered loci for MOSTest (# replicated loci: 242-496, replication rate: 96-97%) in independent data when compared with the established min-P approach (# replicated loci: 26-55, replication rate: 91-93%). An out-of-sample replication of discovered loci was conducted with a sample of 4,069 individuals from the Adolescent Brain Cognitive Development® (ABCD) study, who are on average 50 years younger than UK Biobank individuals. We observe a higher replication yield and comparable replication rate of MOSTest-PVS compared to min-P. This finding underscores the importance of using well-powered multivariate techniques for both discovery and replication of high dimensional phenotypes in Genome-Wide Association studies.

BackgroundClarifying the relationship between depression symptoms and cardiometabolic and related health could clarify risk factors and treatment targets. The objective of this study was to assess whether depression symptoms in midlife are associated with the subsequent onset of cardiometabolic health problems.MethodsThe study sample comprised 787 male twin veterans with polygenic risk score data who participated in the Harvard Twin Study of Substance Abuse ('baseline') and the longitudinal Vietnam Era Twin Study of Aging ('follow-up'). Depression symptoms were assessed at baseline [mean age 41.42 years (s.d. = 2.34)] using the Diagnostic Interview Schedule, Version III, Revised. The onset of eight cardiometabolic conditions (atrial fibrillation, diabetes, erectile dysfunction, hypercholesterolemia, hypertension, myocardial infarction, sleep apnea, and stroke) was assessed via self-reported doctor diagnosis at follow-up [mean age 67.59 years (s.d. = 2.41)].ResultsTotal depression symptoms were longitudinally associated with incident diabetes (OR 1.29, 95% CI 1.07-1.57), erectile dysfunction (OR 1.32, 95% CI 1.10-1.59), hypercholesterolemia (OR 1.26, 95% CI 1.04-1.53), and sleep apnea (OR 1.40, 95% CI 1.13-1.74) over 27 years after controlling for age, alcohol consumption, smoking, body mass index, C-reactive protein, and polygenic risk for specific health conditions. In sensitivity analyses that excluded somatic depression symptoms, only the association with sleep apnea remained significant (OR 1.32, 95% CI 1.09-1.60).ConclusionsA history of depression symptoms by early midlife is associated with an elevated risk for subsequent development of several self-reported health conditions. When isolated, non-somatic depression symptoms are associated with incident self-reported sleep apnea. Depression symptom history may be a predictor or marker of cardiometabolic risk over decades.

The locus coeruleus (LC) is one of the earliest sites of tau pathology, making it a key structure in early Alzheimer's disease (AD) progression. As the primary source of norepinephrine for the brain, reduced LC integrity may have negative consequences for brain health, yet macrostructural brain measures (e.g. cortical thickness) may not be sensitive to early stages of neurodegeneration. We therefore examined whether LC integrity was associated with differences in cortical gray matter microstructure among 435 men (mean age = 67.5; range = 62-71.7). LC structural integrity was indexed by contrast-to-noise ratio (LCCNR) from a neuromelanin-sensitive MRI scan. Restriction spectrum imaging (RSI), an advanced multi-shell diffusion technique, was used to characterize cortical microstructure, modeling total diffusion in restricted, hindered, and free water compartments. Higher LCCNR (greater integrity) was associated with higher hindered and lower free water diffusion in multiple cortical regions. In contrast, no associations between LCCNR and cortical thickness survived correction. Results suggest lower LC integrity is associated with patterns of cortical microstructure that may reflect a reduction in cytoarchitectural barriers due to broader neurodegenerative processes. These findings highlight the potential utility for LC imaging and advanced diffusion measures of cortical microstructure in assessing brain health and early identification of neurodegenerative processes.

BackgroundWe previously developed an African-ancestry-specific polygenic hazard score (PHS46+African) that substantially improved prostate cancer risk stratification in men with African ancestry. The model consists of 46 SNPs identified in Europeans and 3 SNPs from 8q24 shown to improve model performance in Africans. Herein, we used principal component (PC) analysis to uncover subpopulations of men with African ancestry for whom the utility of PHS46+African may differ.Materials and methodsGenotypic data were obtained from the PRACTICAL consortium for 6253 men with African genetic ancestry. Genetic variation in a window spanning 3 African-specific 8q24 SNPs was estimated using 93 PCs. A Cox proportional hazards framework was used to identify the pair of PCs most strongly associated with the performance of PHS46+African. A calibration factor (CF) was formulated using Cox coefficients to quantify the extent to which the performance of PHS46+African varies with PC.ResultsCF of PHS46+African was strongly associated with the first and twentieth PCs. Predicted CF ranged from 0.41 to 2.94, suggesting that PHS46+African may be up to 7 times more beneficial to some African men than others. The explained relative risk for PHS46+African varied from 3.6% to 9.9% for individuals with low and high CF values, respectively. By cross-referencing our data set with 1000 Genomes, we identified significant associations between continental and calibration groupings.ConclusionWe identified PCs within 8q24 that were strongly associated with the performance of PHS46+African. Further research to improve the clinical utility of polygenic risk scores (or models) is needed to improve health outcomes for men of African ancestry.

Background and aimsSmoking is associated with increased risk for brain aging/atrophy and dementia. Few studies have examined early associations with brain aging. This study aimed to measure whether adult men with a history of heavier smoking in early mid-life would have older than predicted brain age 16-28 years later.DesignProspective cohort observational study, utilizing smoking pack years data from average age 40 (early mid-life) predicting predicted brain age difference scores (PBAD) at average ages 56, 62 (later mid-life) and 68 years (early old age). Early mid-life alcohol use was also evaluated.SettingPopulation-based United States sample.Participants/casesParticipants were male twins of predominantly European ancestry who served in the United States military between 1965 and 1975. Structural magnetic resonance imaging (MRI) began at average age 56. Subsequent study waves included most baseline participants; attrition replacement subjects were added at later waves.MeasurementsSelf-reported smoking information was used to calculate pack years smoked at ages 40, 56, 62, and 68. MRIs were processed with the Brain-Age Regression Analysis and Computation Utility software (BARACUS) program to create PBAD scores (chronological age-predicted brain age) acquired at average ages 56 (n = 493; 2002-08), 62 (n = 408; 2009-14) and 68 (n = 499; 2016-19).FindingsIn structural equation modeling, age 40 pack years predicted more advanced age 56 PBAD [β = -0.144, P = 0.012, 95% confidence interval (CI) = -0.257, -0.032]. Age 40 pack years did not additionally predict PBAD at later ages. Age 40 alcohol consumption, but not a smoking × alcohol interaction, predicted more advanced PBAD at age 56 (β = -0.166, P = 0.001, 95% CI = -0.261, -0.070) with additional influences at age 62 (β = -0.115, P = 0.005, 95% CI = -0.195, -0.036). Age 40 alcohol did not predict age 68 PBAD. Within-twin-pair analyses suggested some genetic mechanism partially underlying effects of alcohol, but not smoking, on PBAD.ConclusionsHeavier smoking and alcohol consumption by age 40 appears to predict advanced brain aging by age 56 in men.

BackgroundPosttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).MethodsA GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.ResultsGWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.ConclusionsThrough using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

PurposeRestriction spectrum imaging (RSI) decomposes the diffusion-weighted MRI signal into separate components of known apparent diffusion coefficients (ADCs). The number of diffusion components and optimal ADCs for RSI are organ-specific and determined empirically. The purpose of this work was to determine the RSI model for breast tissues.MethodsThe diffusion-weighted MRI signal was described using a linear combination of multiple exponential components. A set of ADC values was estimated to fit voxels in cancer and control ROIs. Later, the signal contributions of each diffusion component were estimated using these fixed ADC values. Relative-fitting residuals and Bayesian information criterion were assessed. Contrast-to-noise ratio between cancer and fibroglandular tissue in RSI-derived signal contribution maps was compared to DCE imaging.ResultsA total of 74 women with breast cancer were scanned at 3.0 Tesla MRI. The fitting residuals of conventional ADC and Bayesian information criterion suggest that a 3-component model improves the characterization of the diffusion signal over a biexponential model. Estimated ADCs of triexponential model were D1,3 = 0, D2,3 = 1.5 × 10-3 , and D3,3 = 10.8 × 10-3 mm2 /s. The RSI-derived signal contributions of the slower diffusion components were larger in tumors than in fibroglandular tissues. Further, the contrast-to-noise and specificity at 80% sensitivity of DCE and a subset of RSI-derived maps were equivalent.ConclusionBreast diffusion-weighted MRI signal was best described using a triexponential model. Tumor conspicuity in breast RSI model is comparable to that of DCE without the use of exogenous contrast. These data may be used as differential features between healthy and malignant breast tissues.

To determine the impact of genetic variants on the brain, we used genetically informed brain atlases in genome-wide association studies of regional cortical surface area and thickness in 39,898 adults and 9136 children. We uncovered 440 genome-wide significant loci in the discovery cohort and 800 from a post hoc combined meta-analysis. Loci in adulthood were largely captured in childhood, showing signatures of negative selection, and were linked to early neurodevelopment and pathways associated with neuropsychiatric risk. Opposing gradations of decreased surface area and increased thickness were associated with common inversion polymorphisms. Inferior frontal regions, encompassing Broca's area, which is important for speech, were enriched for human-specific genomic elements. Thus, a mixed genetic landscape of conserved and human-specific features is concordant with brain hierarchy and morphogenetic gradients.

During late childhood behavioral changes, such as increased risk-taking and emotional reactivity, have been associated with the maturation of cortico-cortico and cortico-subcortical circuits. Understanding microstructural changes in both white matter and subcortical regions may aid our understanding of how individual differences in these behaviors emerge. Restriction spectrum imaging (RSI) is a framework for modelling diffusion-weighted imaging that decomposes the diffusion signal from a voxel into hindered, restricted, and free compartments. This yields greater specificity than conventional methods of characterizing diffusion. Using RSI, we quantified voxelwise restricted diffusion across the brain and measured age associations in a large sample (n = 8086) from the Adolescent Brain and Cognitive Development (ABCD) study aged 9-14 years. Older participants showed a higher restricted signal fraction across the brain, with the largest associations in subcortical regions, particularly the basal ganglia and ventral diencephalon. Importantly, age associations varied with respect to the cytoarchitecture within white matter fiber tracts and subcortical structures, for example age associations differed across thalamic nuclei. This suggests that age-related changes may map onto specific cell populations or circuits and highlights the utility of voxelwise compared to ROI-wise analyses. Future analyses will aim to understand the relevance of this microstructural developmental for behavioral outcomes.

Because longitudinal studies of aging typically lack cognitive data from earlier ages, it is unclear how general cognitive ability (GCA) changes throughout the life course. In 1173 Vietnam Era Twin Study of Aging (VETSA) participants, we assessed young adult GCA at average age 20 and current GCA at 3 VETSA assessments beginning at average age 56. The same GCA index was used throughout. Higher young adult GCA and better GCA maintenance were associated with stronger specific cognitive abilities from age 51 to 73. Given equivalent GCA at age 56, individuals who had higher age 20 GCA outperformed those whose GCA remained stable in terms of memory, executive function, and working memory abilities from age 51 to 73. Thus, paradoxically, despite poorer maintenance of GCA, high young adult GCA still conferred benefits. Advanced predicted brain age and the combination of elevated vascular burden and APOE-ε4 status were associated with poorer maintenance of GCA. These findings highlight the importance of distinguishing between peak and current GCA for greater understanding of cognitive aging.

BackgroundParkinson's disease (PD) is a highly age-related disorder, where common genetic risk variants affect both disease risk and age at onset. A statistical approach that integrates these effects across all common variants may be clinically useful for individual risk stratification. A polygenic hazard score methodology, leveraging a time-to-event framework, has recently been successfully applied in other age-related disorders.ObjectivesWe aimed to develop and validate a polygenic hazard score model in sporadic PD.MethodsUsing a Cox regression framework, we modeled the polygenic hazard score in a training data set of 11,693 PD patients and 9841 controls. The score was then validated in an independent test data set of 5112 PD patients and 5372 controls and a small single-study sample of 360 patients and 160 controls.ResultsA polygenic hazard score predicts the onset of PD with a hazard ratio of 3.78 (95% confidence interval 3.49-4.10) when comparing the highest to the lowest risk decile. Combined with epidemiological data on incidence rate, we apply the score to estimate genetically stratified instantaneous PD risk across age groups.ConclusionsWe demonstrate the feasibility of a polygenic hazard approach in PD, integrating the genetic effects on disease risk and age at onset in a single model. In combination with other predictive biomarkers, the approach may hold promise for risk stratification in future clinical trials of disease-modifying therapies, which aim at postponing the onset of PD. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Diffusion-weighted MRI (DW-MRI) offers a potential adjunct to dynamic contrast-enhanced MRI to discriminate benign from malignant breast lesions by yielding quantitative information about tissue microstructure. Multi-component modeling of the DW-MRI signal over an extended b-value range (up to 3000 s/mm2) theoretically isolates the slowly diffusing (restricted) water component in tissues. Previously, a three-component restriction spectrum imaging (RSI) model demonstrated the ability to distinguish malignant lesions from healthy breast tissue. We further evaluated the utility of this three-component model to differentiate malignant from benign lesions and healthy tissue in 12 patients with known malignancy and synchronous pathology-proven benign lesions. The signal contributions from three distinct diffusion compartments were measured to generate parametric maps corresponding to diffusivity on a voxel-wise basis. The three-component model discriminated malignant from benign and healthy tissue, particularly using the restricted diffusion C1 compartment and product of the restricted and intermediate diffusion compartments (C1 and C2). However, benign lesions and healthy tissue did not significantly differ in diffusion characteristics. Quantitative discrimination of these three tissue types (malignant, benign, and healthy) in non-pre-defined lesions may enhance the clinical utility of DW-MRI in reducing excessive biopsies and aiding in surveillance and surgical evaluation without repeated exposure to gadolinium contrast.

The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discover 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% are validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development are found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders.

Age has a major effect on brain volume. However, the normative studies available are constrained by small sample sizes, restricted age coverage and significant methodological variability. These limitations introduce inconsistencies and may obscure or distort the lifespan trajectories of brain morphometry. In response, we capitalized on the resources of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to examine age-related trajectories inferred from cross-sectional measures of the ventricles, the basal ganglia (caudate, putamen, pallidum, and nucleus accumbens), the thalamus, hippocampus and amygdala using magnetic resonance imaging data obtained from 18,605 individuals aged 3-90 years. All subcortical structure volumes were at their maximum value early in life. The volume of the basal ganglia showed a monotonic negative association with age thereafter; there was no significant association between age and the volumes of the thalamus, amygdala and the hippocampus (with some degree of decline in thalamus) until the sixth decade of life after which they also showed a steep negative association with age. The lateral ventricles showed continuous enlargement throughout the lifespan. Age was positively associated with inter-individual variability in the hippocampus and amygdala and the lateral ventricles. These results were robust to potential confounders and could be used to examine the functional significance of deviations from typical age-related morphometric patterns.

Delineating the association of age and cortical thickness in healthy individuals is critical given the association of cortical thickness with cognition and behavior. Previous research has shown that robust estimates of the association between age and brain morphometry require large-scale studies. In response, we used cross-sectional data from 17,075 individuals aged 3-90 years from the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to infer age-related changes in cortical thickness. We used fractional polynomial (FP) regression to quantify the association between age and cortical thickness, and we computed normalized growth centiles using the parametric Lambda, Mu, and Sigma method. Interindividual variability was estimated using meta-analysis and one-way analysis of variance. For most regions, their highest cortical thickness value was observed in childhood. Age and cortical thickness showed a negative association; the slope was steeper up to the third decade of life and more gradual thereafter; notable exceptions to this general pattern were entorhinal, temporopolar, and anterior cingulate cortices. Interindividual variability was largest in temporal and frontal regions across the lifespan. Age and its FP combinations explained up to 59% variance in cortical thickness. These results may form the basis of further investigation on normative deviation in cortical thickness and its significance for behavioral and cognitive outcomes.