Beyond linearity in neuroimaging: Capturing nonlinear relationships with application to longitudinal studies.

The ubiquitous adoption of linearity for quantitative predictors in statistical modeling is likely attributable to its advantages of straightforward interpretation and computational feasibility. The linearity assumption may be a reasonable approximation especially when the variable is confined within a narrow range, but it can be problematic when the variable's effect is non-monotonic or complex. Furthermore, visualization and model assessment of a linear fit are usually omitted because of challenges at the whole brain level in neuroimaging. By adopting a principle of learning from the data in the presence of uncertainty to resolve the problematic aspects of conventional polynomial fitting, we introduce a flexible and adaptive approach of multilevel smoothing splines (MSS) to capture any nonlinearity of a quantitative predictor for population-level neuroimaging data analysis. With no prior knowledge regarding the underlying relationship other than a parsimonious assumption about the extent of smoothness (e.g., no sharp corners), we express the unknown relationship with a sufficient number of smoothing splines and use the data to adaptively determine the specifics of the nonlinearity. In addition to introducing the theoretical framework of MSS as an efficient approach with a counterbalance between flexibility and stability, we strive to (a) lay out the specific schemes for population-level nonlinear analyses that may involve task (e.g., contrasting conditions) and subject-grouping (e.g., patients vs controls) factors; (b) provide modeling accommodations to adaptively reveal, estimate and compare any nonlinear effects of a predictor across the brain, or to more accurately account for the effects (including nonlinear effects) of a quantitative confound; (c) offer the associated program 3dMSS to the neuroimaging community for whole-brain voxel-wise analysis as part of the AFNI suite; and (d) demonstrate the modeling approach and visualization processes with a longitudinal dataset of structural MRI scans.

The NIMH Intramural Longitudinal Study of the Endocrine and Neurobiological Events Accompanying Puberty: Protocol and rationale for methods and measures.

Delineating the relationship between human neurodevelopment and the maturation of the hypothalamic-pituitary-gonadal (HPG) axis during puberty is critical for investigating the increase in vulnerability to neuropsychiatric disorders that is well documented during this period. Preclinical research demonstrates a clear association between gonadal production of sex steroids and neurodevelopment; however, identifying similar associations in humans has been complicated by confounding variables (such as age) and the coactivation of two additional endocrine systems (the adrenal androgenic system and the somatotropic growth axis) and requires further elucidation. In this paper, we present the design of, and preliminary observations from, the ongoing NIMH Intramural Longitudinal Study of the Endocrine and Neurobiological Events Accompanying Puberty. The aim of this study is to directly examine how the increase in sex steroid hormone production following activation of the HPG-axis (i.e., gonadarche) impacts neurodevelopment, and, additionally, to determine how gonadal development and maturation is associated with longitudinal changes in brain structure and function in boys and girls. To disentangle the effects of sex steroids from those of age and other endocrine events on brain development, our study design includes 1) selection criteria that establish a well-characterized baseline cohort of healthy 8-year-old children prior to the onset of puberty (e.g., prior to puberty-related sex steroid hormone production); 2) temporally dense longitudinal, repeated-measures sampling of typically developing children at 8-10 month intervals over a 10-year period between the ages of eight and 18; 3) contemporaneous collection of endocrine and other measures of gonadal, adrenal, and growth axis function at each timepoint; and 4) collection of multimodal neuroimaging measures at these same timepoints, including brain structure (gray and white matter volume, cortical thickness and area, white matter integrity, myelination) and function (reward processing, emotional processing, inhibition/impulsivity, working memory, resting-state network connectivity, regional cerebral blood flow). This report of our ongoing longitudinal study 1) provides a comprehensive review of the endocrine events of puberty; 2) details our overall study design; 3) presents our selection criteria for study entry (e.g., well-characterized prepubertal baseline) along with the endocrinological considerations and guiding principles that underlie these criteria; 4) describes our longitudinal outcome measures and how they specifically relate to investigating the effects of gonadal development on brain development; and 5) documents patterns of fMRI activation and resting-state networks from an early, representative subsample of our cohort of prepubertal 8-year-old children.

Brain-Derived Neurotrophic Factor Val66Met Polymorphism Affects the Relationship Between an Anxiety-Related Personality Trait and Resting Regional Cerebral Blood Flow.

Brain-derived neurotrophic factor (BDNF) is an important modulator of constitutive stress responses mediated by limbic frontotemporal circuits, and its gene contains a functional polymorphism (Val66Met) that may influence trait stress sensitivity. Reports of an association of this polymorphism with anxiety-related personality traits have been controversial and without clear neurophysiological support. We, therefore, determined the relationship between resting regional cerebral blood flow (rCBF) and a well-validated measure of anxiety-related personality, the TPQ Harm Avoidance (HA) scale, as a function of BDNF Val66Met genotype. Sixty-four healthy participants of European ancestry underwent resting H215O positron emission tomography scans. For each genotype group separately, we first determined the relationship between participants' HA scores and their resting rCBF values in each voxel across the entire brain, and then directly compared these HA-rCBF relationships between Val66Met genotype groups. HA-rCBF relationships differed between Val homozygotes and Met carriers in several regions relevant to stress regulation: subgenual cingulate, orbital frontal cortex, and the hippocampal/parahippocampal region. In each of these areas, the relationship was positive in Val homozygotes and negative in Met carriers. These data demonstrate a coupling between trait anxiety and basal resting blood flow in frontolimbic neurocircuitry that may be determined in part by genetically mediated BDNF signaling.

Altered pubertal timing in 7q11.23 copy number variations and associated genetic mechanisms.

Pubertal timing, including age at menarche (AAM), is a heritable trait linked to lifetime health outcomes. Here, we investigate genetic mechanisms underlying AAM by combining genome-wide association study (GWAS) data with investigations of two rare genetic conditions clinically associated with altered AAM: Williams syndrome (WS), a 7q11.23 hemideletion characterized by early puberty; and duplication of the same genes (7q11.23 Duplication syndrome [Dup7]) characterized by delayed puberty. First, we confirm that AAM-derived polygenic scores in typically developing children (TD) explain a modest amount of variance in AAM (R2 = 0.09; p = 0.04). Next, we demonstrate that 7q11.23 copy number impacts AAM (WS < TD < Dup7; p = 1.2x10-8, η2 = 0.45) and pituitary volume (WS < TD < Dup7; p = 3x10-5, ηp2 = 0.2) with greater effect sizes. Finally, we relate an AAM-GWAS signal in 7q11.23 to altered expression in postmortem brains of STAG3L2 (p = 1.7x10-17), a gene we also find differentially expressed with 7q11.23 copy number (p = 0.03). Collectively, these data explicate the role of 7q11.23 in pubertal onset, with STAG3L2 and pituitary development as potential mediators.

Brain-derived neurotrophic factor Val66Met polymorphism affects resting regional cerebral blood flow and functional connectivity differentially in women versus men.

The human Val66Met single nucleotide polymorphism in the brain-derived neurotrophic factor (BDNF) gene impacts BDNF signaling at the cellular level. At the neural-systems level, it is associated with differences in prefrontal cortex (PFC) and hippocampal function during performance of cognitive and affective tasks. Because the impact of this variant on basal prefrontal and hippocampal activity is not known but may be relevant to understanding the function of this gene in health and disease, we studied 94 healthy individuals with H2 ¹5O PET to assess regional cerebral blood flow (rCBF) during rest and tested for between-genotype differences. Because BDNF and gonadal steroid hormones conjointly influence neuronal growth, survival, and plasticity in hippocampus and PFC, we also tested for sex × genotype interactions. Finally, in light of the known impact of BDNF on plasticity and dendritic arborization, we complimented direct rCBF comparisons with connectivity analyses to determine how activity in hippocampal and prefrontal regions showing between-genotype group differences covaries with rCBF in other nodes throughout the brain in a genotype- or sex-dependent manner. Compared with Val homozygotes, Met carriers had higher rCBF in prefrontal (BA25 extending into BA10) and hippocampal/parahippocampal regions. Moreover, there were significant sex × genotype interactions in regions (including frontal, parahippocampal, and lateral temporal cortex) in which Val homozygotes showed higher rCBF in females than males, but Met carriers showed the opposite relationship. Functional connectivity analysis demonstrated that correlations of BA25, hippocampus, and parahippocampus with frontal and temporal networks were positive for Val homozygotes and negative for Met carriers. In addition, sex × genotype analysis of functional connectivity revealed that genotype affected directionality of the inter-regional correlations differentially in men versus women. Our data indicate that BDNF allelic variation and sex interactively affect basal prefrontal and hippocampal function.

The short-term effects of estradiol, raloxifene, and a phytoestrogen in women with perimenopausal depression.

OBJECTIVE: We examined the short-term efficacies of three estrogen-like compounds under placebo-controlled conditions in women with perimenopause-related depression (PMD). METHODS: Women with PMD were randomized in a double-blind parallel design to one of four treatments: transdermal 17-beta estradiol (TE) (100 mcg/d); oral raloxifene (60 mg/d); a proprietary phytoestrogen compound, Rimostil (1,000 mg twice/d); or placebo for 8 weeks. The main outcome measures were the Center for Epidemiology Studies Depression Scale, 17-item Hamilton Rating Scale for Depression (HRSD), and the Beck Depression Inventory completed at each clinic visit. Secondary outcomes included a visual analogue self-rating completed at each clinic visit, and daily self-ratings of hot flush severity. Cognitive tests were performed at pretreatment baseline and at the end of the trial. In the primary analysis, we obtained four repeated measures in each woman in the four treatment arms. Analyses were done with SAS Version 9.4 software (SAS Institute, Inc, Cary, NC), using PROC MIXED (for mixed models). All models included the following four explanatory variables, regardless of whether they were statistically significant: 1) treatment group (TE, raloxifene, Rimostil, placebo); 2) week (W2, W4, W6, W8); 3) treatment group-by-week interaction; and 4) baseline value of the measure being analyzed. The inclusion of additional variables was evaluated individually for each outcome measure. RESULTS: Sixty-six women were randomized into the trial, four women dropped out of the trial, and 62 women were included in the final data analysis. No effect of treatment group was observed in either the Center for Epidemiology Studies Depression Scale (P = 0.34) or Beck Depression Inventory (P = 0.27) scores; however, there was a difference in HRSD scores between treatment groups (P = 0.0037) that pair-wise comparisons of the combined weekly scores in each treatment demonstrated TE's beneficial effects on HRSD scores compared with Rimostil (P = 0.0005), and less consistently with placebo (P = 0.099). The average (SD) of the baseline scores for each treatment group on the HRSD was as follows: TE-15.3 (4.5), raloxifene-16.0 (3.7), Rimostil-14.0 (2.7), and placebo-15.2 (3.0). Whereas the HRSD scores after 8 weeks of treatment (least-square means) were TE-5.2(1.1), raloxifene-5.8(1.2), Rimostil-11.2(1.4), and placebo-7.8(1.1). No differences were observed between raloxifene and either TE or placebo in any scale score. HRSD scores in women assigned to TE were improved compared with those on Rimostil during weeks 6 and 8 (P values = 0.0008, 0.0011, respectively). Cognitive testing at week 8 showed that none of the three active treatment groups performed better than placebo. CONCLUSIONS: This study did not identify significant therapeutic benefits of TE, Rimostil, or raloxifene compared with placebo in PMD. However, improvements in depression ratings were observed between TE compared with Rimostil. Thus, our findings do not support the role of ERbeta compounds in the treatment of PMD (and indeed could suggest a more important role of ERalpha).

The Cortisol and ACTH Response to Dex/CRH Testing in Women With and Without Perimenopausal Depression.

CONTEXT: Abnormalities in the hypothalamic-pituitary-adrenal (HPA) axis are frequent accompaniments of depression, and studies have documented the role of stress and stressful life events in the ontogeny of perimenopausal depressions (PMD). Because HPA axis function in women is further modulated both by aging and ovarian steroids, it is possible that a dysregulated HPA axis contributes to the increased risk of PMD. OBJECTIVE: We examined HPA axis function in perimenopausal women with and without depression using the combined dexamethasone-corticotropin-releasing hormone (Dex/CRH) test. METHODS: Dex/CRH tests were performed on 20 women with PMD and 20 women who were also perimenopausal but without current or past depression (control women). Main outcome measures were plasma levels of cortisol and adrenocorticotropin (ACTH) and 24-hour urinary free cortisol (UFC). Five women took chronic stable medications, otherwise all women were medically healthy, and both groups were comparable with respect to reproductive stage and age. Standardized symptom rating scales were administered to each woman prior to Dex/CRH testing. RESULTS: No group differences were present in either baseline or stimulated ACTH and cortisol secretion. Baseline plasma measures of estradiol, progesterone, and 24-hour UFC levels similarly did not differ in PMD and control women. CONCLUSION: Despite reports of increased stress responsiveness in PMD, we observed no abnormalities of HPA axis activity associated with PMD compared with women without depression. These findings suggest that PMD is not uniformly associated with HPA dysregulation and could reflect underlying pathophysiologic processes that are distinct from women with nonreproductive-related depressions.

Subgenual cingulate resting regional cerebral blood flow in premenstrual dysphoric disorder: differential regulation by ovarian steroids and preliminary evidence for an association with expression of ESC/E(Z) complex genes.

Substantial evidence suggests that circulating ovarian steroids modulate behavior differently in women with PMDD than in those without this condition. However, hormonal state-related abnormalities of neural functioning in PMDD remain to be better characterized. In addition, while altered neural function in PMDD likely co-exists with alterations in intrinsic cellular function, such a relationship has not been explored. Here, we investigated the effects of ovarian steroids on basal, resting regional cerebral blood flow (rCBF) in PMDD, and, in an exploratory analysis, we tested whether the rCBF findings were linked to the expression of ESC/E(Z) genes, which form an essential ovarian steroid-regulated gene-silencing complex. Resting rCBF was measured with oxygen-15 water PET (189 PET sessions in 43 healthy women and 20 women with PMDD) during three self-as-own-control conditions: GnRH agonist (Lupron)-induced ovarian suppression, estradiol add-back, and progesterone add-back. ESC/E(Z) gene expression data were obtained from RNA-sequencing of lymphoblastoid cell lines performed in a previous study and were examined in relation to hormone-induced changes in rCBF. In the rCBF PET data, there was a significant diagnosis-by-hormone interaction in the subgenual cingulate (PFDR = 0.05), an important neuroanatomical hub for regulating affective state. Whereas control women showed no hormonally-related changes in resting rCBF, those with PMDD showed decreased resting rCBF during both estradiol (P = 0.02) and progesterone (P = 0.0002) add-back conditions. In addition, in PMDD, ESC/E(Z) gene expression correlated with the change in resting rCBF between Lupron-alone and progesterone conditions (Pearson r = -0.807, P = 0.016). This work offers a formulation of PMDD that integrates behavioral, neural circuit, and cellular mechanisms, and may provide new targets for future therapeutic interventions.

Abnormalities of dorsolateral prefrontal function in women with premenstrual dysphoric disorder: a multimodal neuroimaging study.

OBJECTIVE: To investigate the neural substrate of premenstrual dysphoric disorder (PMDD), the authors used [15O]H2O positron emission tomography (PET) regional cerebral blood flow (rCBF) and blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) signal measurements during working memory in conjunction with a 6-month hormone manipulation protocol. METHOD: PET and fMRI scans were obtained from women with prospectively confirmed PMDD and asymptomatic comparison subjects while they completed the n-back task during three hormone conditions: ovarian suppression induced by the gonadotropin-releasing hormone agonist leuprolide acetate, leuprolide plus estradiol, and leuprolide plus progesterone. Fifteen patients and 15 matched comparison subjects underwent PET imaging. Fourteen patients and 14 comparison subjects underwent fMRI. For each hormone condition, rCBF was measured with [15O]H2O PET, and BOLD signal was measured with fMRI, both during an n-back working memory paradigm. Global Assessment of Functioning Scale (GAF) scores and clinical characteristics were obtained for each patient before hormone manipulation, and symptoms were measured before and during the protocol. RESULTS: In both the PET and fMRI studies, a main effect of diagnosis was observed, with PMDD patients showing greater prefrontal activation than comparison subjects. In the patient group, the degree to which dorsolateral prefrontal cortex activation was abnormally increased correlated with several dimensions of disease: disability as indicated by GAF scores, age at symptom onset, duration of PMDD, and differences in pre- and postmenses PMDD symptoms. CONCLUSIONS: Abnormal working memory activation in PMDD, specifically in the dorsolateral prefrontal cortex, is related to PMDD severity, symptoms, age at onset, and disease burden. These results support the clinical relevance of the findings and the proposal that dorsolateral prefrontal cortex dysfunction represents a substrate of risk for PMDD. The concordance of the fMRI and PET data attests to the neurobiological validity of the results.

Functional neuroanatomy of the parahippocampal region in the rat: the perirhinal and postrhinal cortices.

The parahippocampal region in the rodent brain includes the perirhinal, postrhinal, and entorhinal cortices, the presubiculum, and the parasubiculum. In recent years, the perirhinal and postrhinal cortices have been a focus in memory research because they supply highly processed, polymodal sensory information to the hippocampus, both directly and via the entorhinal cortex. Available evidence indicates that these cortices receive different complements of cortical information, which are then forwarded to the hippocampus via parallel pathways. Here we have summarized the cortical, subcortical, and hippocampal connections of the perirhinal and postrhinal cortices in order to provide further insight into the nature of the information that is processed by these regions prior to arriving in the hippocampus. As has been previously described, the cortical afferents of the rodent postrhinal cortex are dominated by structures known to be involved in the processing of visual and spatial information, whereas the cortical afferents of the perirhinal cortex result in remarkable convergence of polymodal sensory information. The two regions are also differentiated by their cortical efferents. The perirhinal cortex projects more strongly to piriform, frontal, and insular regions, whereas the postrhinal cortex projects preferentially to visual and visuospatial regions. The subcortical connections of the two regions provide further evidence that they have different functions. For example, the perirhinal cortex has strong reciprocal connections with the amygdala, which suggest involvement in processing affective stimuli. Subcortical input to the postrhinal cortex is dominated by projections from dorsal thalamic structures, particularly the lateral posterior nucleus. Although the perirhinal and postrhinal cortices are considered to contribute to the episodic memory system, many questions remain about their particular roles. A detailed description of the anatomical connections of the perirhinal and postrhinal cortices will permit the generation of new, anatomically guided, hypotheses about their role in episodic memory and other cognitive processes.