Leveraging the adolescent brain cognitive development study to improve behavioral prediction from neuroimaging in smaller replication samples.

Neuroimaging is a popular method to map brain structural and functional patterns to complex human traits. Recently published observations cast doubt upon these prospects, particularly for prediction of cognitive traits from structural and resting state functional magnetic resonance imaging (MRI). We leverage baseline data from thousands of children in the Adolescent Brain Cognitive DevelopmentSM Study to inform the replication sample size required with univariate and multivariate methods across different imaging modalities to detect reproducible brain-behavior associations. We demonstrate that by applying multivariate methods to high-dimensional brain imaging data, we can capture lower dimensional patterns of structural and functional brain architecture that correlate robustly with cognitive phenotypes and are reproducible with only 41 individuals in the replication sample for working memory-related functional MRI, and ~ 100 subjects for structural and resting state MRI. Even with 100 random re-samplings of 100 subjects in discovery, prediction can be adequately powered with 66 subjects in replication for multivariate prediction of cognition with working memory task functional MRI. These results point to an important role for neuroimaging in translational neurodevelopmental research and showcase how findings in large samples can inform reproducible brain-behavior associations in small sample sizes that are at the heart of many research programs and grants.

Early cortical microstructural changes in aging are linked to vulnerability to Alzheimer's disease pathology.

BACKGROUND: Early identification of Alzheimer's disease (AD) risk is critical for improving treatment success. Cortical thickness is a macrostructural measure used to assess AD neurodegeneration. However, cortical microstructural changes appear to precede macrostructural atrophy and may improve early risk identification. Currently, whether cortical microstructural changes in aging are linked to vulnerability to AD pathophysiology remains unclear in non-clinical populations, who are precisely the target for early risk identification. METHODS: In 194 adults, we calculated MRI-derived maps of changes in cortical mean diffusivity (microstructure) and cortical thickness (macrostructure) over 5-6 years (meanage: Time1=61.82; Time2=67.48). Episodic memory was assessed using three well-established tests. We obtained PET-derived maps of AD pathology deposition (beta-amyloid, tau) and neurotransmitter receptors (cholinergic, glutamatergic) implicated in AD pathophysiology. Spatial correlational analyses were used to compare pattern similarity among maps. RESULTS: Spatial patterns of cortical macrostructural changes resembled patterns of cortical organization sensitive to age-related processes (r=-0.31, p<0.05), whereas microstructural changes resembled the patterns of tau (r=0.39, p=0.015) deposition in AD. Individuals with patterns of microstructural changes that more closely resembled stereotypical tau deposition exhibited greater memory decline (ß=0.21, p=0.036). Microstructural changes and AD pathology deposition were enriched in areas with greater densities of cholinergic and glutamatergic receptors (ps<0.05). CONCLUSIONS: Patterns of cortical microstructural changes were more AD-like than patterns of macrostructural changes, which appeared to reflect more general aging processes. Microstructural changes may better inform early risk prediction efforts as a sensitive measure of vulnerability to pathological processes prior to overt atrophy and cognitive decline.

Partitioning variance in cortical morphometry into genetic, environmental, and subject-specific components.

The relative contributions of genetic variation and experience in shaping the morphology of the adolescent brain are not fully understood. Using longitudinal data from 11,665 subjects in the ABCD Study, we fit vertex-wise variance components including family effects, genetic effects, and subject-level effects using a computationally efficient framework. Variance in cortical thickness and surface area is largely attributable to genetic influence, whereas sulcal depth is primarily explained by subject-level effects. Our results identify areas with heterogeneous distributions of heritability estimates that have not been seen in previous work using data from cortical regions. We discuss the biological importance of subject-specific variance and its implications for environmental influences on cortical development and maturation.

Improved functional mapping of complex trait heritability with GSA-MiXeR implicates biologically specific gene sets.

While genome-wide association studies are increasingly successful in discovering genomic loci associated with complex human traits and disorders, the biological interpretation of these findings remains challenging. Here we developed the GSA-MiXeR analytical tool for gene set analysis (GSA), which fits a model for the heritability of individual genes, accounting for linkage disequilibrium across variants and allowing the quantification of partitioned heritability and fold enrichment for small gene sets. We validated the method using extensive simulations and sensitivity analyses. When applied to a diverse selection of complex traits and disorders, including schizophrenia, GSA-MiXeR prioritizes gene sets with greater biological specificity compared to standard GSA approaches, implicating voltage-gated calcium channel function and dopaminergic signaling for schizophrenia. Such biologically relevant gene sets, often with fewer than ten genes, are more likely to provide insights into the pathobiology of complex diseases and highlight potential drug targets.

Quality over quantity: powering neuroimaging samples in psychiatry.

Neuroimaging has been widely adopted in psychiatric research, with hopes that these non-invasive methods will provide important clues to the underpinnings and prediction of various mental health symptoms and outcomes. However, the translational impact of neuroimaging has not yet reached its promise, despite the plethora of computational methods, tools, and datasets at our disposal. Some have lamented that too many psychiatric neuroimaging studies have been underpowered with respect to sample size. In this review, we encourage this discourse to shift from a focus on sheer increases in sample size to more thoughtful choices surrounding experimental study designs. We propose considerations at multiple decision points throughout the study design, data modeling and analysis process that may help researchers working in psychiatric neuroimaging boost power for their research questions of interest without necessarily increasing sample size. We also provide suggestions for leveraging multiple datasets to inform each other and strengthen our confidence in the generalization of findings to both population-level and clinical samples. Through a greater emphasis on improving the quality of brain-based and clinical measures rather than merely quantity, meaningful and potentially translational clinical associations with neuroimaging measures can be achieved with more modest sample sizes in psychiatry.

Clinical Impact of Contouring Variability for Prostate Cancer Tumor Boost.

PURPOSE: The focal radiotherapy (RT) boost technique was shown in a phase III randomized controlled trial (RCT) to improve prostate cancer outcomes without increasing toxicity. This technique relies on the accurate delineation of prostate tumors on MRI. A recent prospective study evaluated radiation oncologists' accuracy when asked to delineate prostate tumors on MRI and demonstrated high variability in tumor contours. We sought to evaluate the impact of contour variability and inaccuracy on predicted clinical outcomes. We hypothesized that radiation oncologists' contour inaccuracies would yield meaningfully worse clinical outcomes. MATERIALS & METHODS: 45 radiation oncologists and 2 expert radiologists contoured prostate tumors on 30 patient cases. Of these cases, those with CT simulation or diagnostic CT available were selected for analysis. A knowledge-based planning model was developed to generate focal RT boost plans for each contour per the RCT protocol. Probability of biochemical failure (BF) was determined using a model from the RCT. The primary metric evaluated was delta BF (ΔBF = Participant BF - Expert BF). An absolute increase in BF ≥5% was considered clinically meaningful. RESULTS: 8 patient cases and 394 target volumes for focal RT boost planning were included in this analysis. In general, participant plans were associated with worse predicted clinical outcomes compared to the expert plan, with an average absolute increase in BF of 4.3%. 37% of participant plans were noted to have an absolute increase in BF of 5% or more. CONCLUSION: Radiation oncologists' attempts to contour tumor targets for focal RT boost are frequently inaccurate enough to yield meaningfully inferior clinical outcomes for patients.

Dissecting the Shared Genetic Architecture of Common Epilepsies With Cortical Brain Morphology.

Background and Objectives: Epilepsies are associated with differences in cortical thickness (TH) and surface area (SA). However, the mechanisms underlying these relationships remain elusive. We investigated the extent to which these phenotypes share genetic influences. Methods: We analyzed genome-wide association study data on common epilepsies (n = 69,995) and TH and SA (n = 32,877) using Gaussian mixture modeling MiXeR and conjunctional false discovery rate (conjFDR) analysis to quantify their shared genetic architecture and identify overlapping loci. We biologically interrogated the loci using a variety of resources and validated in independent samples. Results: The epilepsies (2.4 k-2.9 k variants) were more polygenic than both SA (1.8 k variants) and TH (1.3 k variants). Despite absent genome-wide genetic correlations, there was a substantial genetic overlap between SA and genetic generalized epilepsy (GGE) (1.1 k), all epilepsies (1.1 k), and juvenile myoclonic epilepsy (JME) (0.7 k), as well as between TH and GGE (0.8 k), all epilepsies (0.7 k), and JME (0.8 k), estimated with MiXeR. Furthermore, conjFDR analysis identified 15 GGE loci jointly associated with SA and 15 with TH, 3 loci shared between SA and childhood absence epilepsy, and 6 loci overlapping between SA and JME. 23 loci were novel for epilepsies and 11 for cortical morphology. We observed a high degree of sign concordance in the independent samples. Discussion: Our findings show extensive genetic overlap between generalized epilepsies and cortical morphology, indicating a complex genetic relationship with mixed-effect directions. The results suggest that shared genetic influences may contribute to cortical abnormalities in epilepsies.

Unraveling the shared genetics of common epilepsies and general cognitive ability.

Objective: Cognitive impairment is prevalent among individuals with epilepsy, and it is possible that genetic factors can underlie this relationship. Here, we investigated the potential shared genetic basis of common epilepsies and general cognitive ability (COG). Methods: We applied linkage disequilibrium score (LDSC) regression, MiXeR and conjunctional false discovery rate (conjFDR) to analyze different aspects of genetic overlap between COG and epilepsies. We used the largest available genome-wide association study data on COG (n = 269,867) and common epilepsies (n = 27,559 cases, 42,436 controls), including the broad phenotypes 'all epilepsy', focal epilepsies and genetic generalized epilepsies (GGE), and as well as specific subtypes. We functionally annotated the identified loci using a variety of biological resources and validated the results in independent samples. Results: Using MiXeR, COG (11.2k variants) was estimated to be almost four times more polygenic than 'all epilepsy', GGE, juvenile myoclonic epilepsy (JME), and childhood absence epilepsy (CAE) (2.5k - 2.9k variants). The other epilepsy phenotypes were insufficiently powered for analysis. We show extensive genetic overlap between COG and epilepsies with significant negative genetic correlations (-0.23 to -0.04). COG was estimated to share 2.9k variants with both GGE and 'all epilepsy', and 2.3k variants with both JME and CAE. Using conjFDR, we identified 66 distinct loci shared between COG and epilepsies, including novel associations for GGE (27), 'all epilepsy' (5), JME (5) and CAE (5). The implicated genes were significantly expressed in multiple brain regions. The results were validated in independent samples (COG: p = 1.0 × 10-14; 'all epilepsy': p = 5.6 × 10-3). Significance: Our study demonstrates a substantial genetic basis shared between epilepsies and COG and identifies novel overlapping genomic loci. Enhancing our understanding of the relationship between epilepsies and COG may lead to the development of novel comorbidity-targeted epilepsy treatments.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Associations of plasma neurofilament light chain with cognition and neuroimaging measures in community-dwelling early old age men.

BACKGROUND: Plasma 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. METHODS: We 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. RESULTS: After 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. CONCLUSIONS: These 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.

Dissecting the genetic overlap between three complex phenotypes with trivariate MiXeR.

Comorbidities are an increasing global health challenge. Accumulating evidence suggests overlapping genetic architectures underlying comorbid complex human traits and disorders. The bivariate causal mixture model (MiXeR) can quantify the polygenic overlap between complex phenotypes beyond global genetic correlation. Still, the pattern of genetic overlap between three distinct phenotypes, which is important to better characterize multimorbidities, has previously not been possible to quantify. Here, we present and validate the trivariate MiXeR tool, which disentangles the pattern of genetic overlap between three phenotypes using summary statistics from genome-wide association studies (GWAS). Our simulations show that the trivariate MiXeR can reliably reconstruct different patterns of genetic overlap. We further demonstrate how the tool can be used to estimate the proportions of genetic overlap between three phenotypes using real GWAS data, providing examples of complex patterns of genetic overlap between diverse human traits and diseases that could not be deduced from bivariate analyses. This contributes to a better understanding of the etiology of complex phenotypes and the nature of their relationship, which may aid in dissecting comorbidity patterns and their biological underpinnings.

Quantitative MRI biomarker for classification of clinically significant prostate cancer: calibration for reproducibility across echo times.

Background: Restriction Spectrum Imaging restriction score (RSIrs) is a quantitative biomarker for detecting clinically significant prostate cancer (csPCa). However, the quantitative value of the RSIrs is affected by imaging parameters such as echo time (TE). Purpose: The purpose of the present study is to develop a calibration method to account for differences in echo times and facilitate use of RSIrs as a quantitative biomarker for the detection of csPCa. Methods: This study included 197 consecutive patients who underwent MRI and biopsy examination; 97 were diagnosed with csPCa (grade group ≥ 2). RSI data were acquired three times during the same session: twice at minimum TE∼75ms and once at TE=90ms (TEmin 1 , TEmin 2 , and TE90, respectively). A proposed calibration method, trained on patients without csPCa, estimated a linear scaling factor (f) for each of the four diffusion compartments (C) of the RSI signal model. A linear regression model was determined to match C-maps of TE90 to the reference C-maps of TEmin 1 within the interval ranging from 95 th to 99 th percentile of signal intensity within the prostate. RSIrs comparisons were made at 98 th percentile within each patient's prostate. We compared RSIrs from calibrated TE90 (RSIrs TE90corr ) and uncorrected TE90 (RSIrs TE90 ) to RSIrs from reference TEmin 1 (RSIrs TEmin1 ) and repeated TEmin 2 (RSIrs TEmin2 ). Calibration performance was evaluated with sensitivity, specificity, area under the ROC curve, positive predicted value, negative predicted value, and F1-score. Results: Scaling factors for C 1 , C 2 , C 3 , and C 4 were estimated as 1.70, 1.38, 1.03, and 1.19, respectively. In non-csPCa cases, the 98 th percentile of RSIrs TEmin2 and RSIrs TEmin1 differed by 0.27±0.86SI (mean±standard deviation), whereas RSIrs TE90 differed from RSIrs TEmin1 by 1.81±1.20SI. After calibration, this bias was reduced to -0.41±1.20SI, representing a 78% reduction in absolute error. For patients with csPCa, the difference was 0.54±1.98SI between RSIrs TEmin2 and RSIrs TEmin1 and 2.28±2.06SI between RSIrs TE90 and RSIrs TEmin1 . After calibration, the mean difference decreased to -0.86SI, a 38% reduction in absolute error. At the Youden index for patient-level classification of csPCa (8.94SI), RSIrs TEmin1 has a sensitivity of 66% and a specificity of 72%. Prior to calibration, RSIrs TE90 at the same threshold tended to over-diagnose benign cases (sensitivity 44%, specificity 88%). Post-calibration, RSIrs TE90corr performs more similarly to the reference (sensitivity 71%, specificity 62%). Conclusion: The proposed linear calibration method produces similar quantitative biomarker values for acquisitions with different TE, reducing TE-induced error by 78% and 38% for non-csPCa and csPCa, respectively.

FEMA: Fast and efficient mixed-effects algorithm for large sample whole-brain imaging data.

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/.

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis.

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

How Real-World Data Can Facilitate the Development of Precision Medicine Treatment in Psychiatry.

Precision medicine has the ambition to improve treatment response and clinical outcomes through patient stratification and holds great potential for the treatment of mental disorders. However, several important factors are needed to transform current practice into a precision psychiatry framework. Most important are 1) the generation of accessible large real-world training and test data including genomic data integrated from multiple sources, 2) the development and validation of advanced analytical tools for stratification and prediction, and 3) the development of clinically useful management platforms for patient monitoring that can be integrated into health care systems in real-life settings. This narrative review summarizes strategies for obtaining the key elements-well-powered samples from large biobanks integrated with electronic health records and health registry data using novel artificial intelligence algorithms-to predict outcomes in severe mental disorders and translate these models into clinical management and treatment approaches. Key elements are massive mental health data and novel artificial intelligence algorithms. For the clinical translation of these strategies, we discuss a precision medicine platform for improved management of mental disorders. We use cases to illustrate how precision medicine interventions could be brought into psychiatry to improve the clinical outcomes of mental disorders.

Genome-wide analyses reveal shared genetic architecture and novel risk loci between opioid use disorder and general cognitive ability.

BACKGROUND: Opioid use disorder (OUD), a serious health burden worldwide, is associated with lower cognitive function. Recent studies have demonstrated a negative genetic correlation between OUD and general cognitive ability (COG), indicating a shared genetic basis. However, the specific genetic variants involved, and the underlying molecular mechanisms remain poorly understood. Here, we aimed to quantify and identify the genetic basis underlying OUD and COG. METHODS: We quantified the extent of genetic overlap between OUD and COG using a bivariate causal mixture model (MiXeR) and identified specific genetic loci applying conditional/conjunctional FDR. Finally, we investigated biological function and expression of implicated genes using available resources. RESULTS: We estimated that ~94% of OUD variants (4.8k out of 5.1k variants) also influence COG. We identified three novel OUD risk loci and one locus shared between OUD and COG. Loci identified implicated biological substrates in the basal ganglia. CONCLUSION: We provide new insights into the complex genetic risk architecture of OUD and its genetic relationship with COG.

Childhood Disadvantage Moderates Late Midlife Default Mode Network Cortical Microstructure and Visual Memory Association.

BACKGROUND: Childhood 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. METHODS: Cortical 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. RESULTS: Higher 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). CONCLUSIONS: Lower 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.

Identification of Novel Genomic Loci for Anxiety and Extensive Genetic Overlap with Psychiatric Disorders.

Background: Anxiety disorders are prevalent and anxiety symptoms co-occur with many psychiatric disorders. We aimed to identify genomic risk loci associated with anxiety, characterize its genetic architecture, and genetic overlap with psychiatric disorders. Methods: We used the GWAS of anxiety symptoms, schizophrenia, bipolar disorder, major depression, and attention deficit hyperactivity disorder (ADHD). We employed MiXeR and LAVA to characterize the genetic architecture and genetic overlap between the phenotypes. Conditional and conjunctional false discovery rate analyses were performed to boost the identification of genomic loci associated with anxiety and those shared with psychiatric disorders. Gene annotation and gene set analyses were conducted using OpenTargets and FUMA, respectively. Results: Anxiety was polygenic with 12.9k estimated genetic risk variants and overlapped extensively with psychiatric disorders (4.1-11.4k variants). MiXeR and LAVA revealed predominantly positive genetic correlations between anxiety and psychiatric disorders. We identified 114 novel loci for anxiety by conditioning on the psychiatric disorders. We also identified loci shared between anxiety and major depression (n = 47), bipolar disorder (n = 33), schizophrenia (n = 71), and ADHD (n = 20). Genes annotated to anxiety loci exhibit enrichment for a broader range of biological pathways and differential tissue expression in more diverse tissues than those annotated to the shared loci. Conclusions: Anxiety is a highly polygenic phenotype with extensive genetic overlap with psychiatric disorders. These genetic overlaps enabled the identification of novel loci for anxiety. The shared genetic architecture may underlie the extensive cross-disorder comorbidity of anxiety, and the identified genetic loci implicate molecular pathways that may lead to potential drug targets.

Restriction spectrum imaging with elastic image registration for automated evaluation of response to neoadjuvant therapy in breast cancer.

Purpose: Dynamic 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 design: Breast 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. Results: Ten 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. Conclusion: The 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.