Toward a Stochastic Complete Active Space Second-Order Perturbation Theory.

In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.

Identifying vulnerable brain networks associated with Alzheimer's disease risk.

The selective vulnerability of brain networks in individuals at risk for Alzheimer's disease (AD) may help differentiate pathological from normal aging at asymptomatic stages, allowing the implementation of more effective interventions. We used a sample of 72 people across the age span, enriched for the APOE4 genotype to reveal vulnerable networks associated with a composite AD risk factor including age, genotype, and sex. Sparse canonical correlation analysis (CCA) revealed a high weight associated with genotype, and subgraphs involving the cuneus, temporal, cingulate cortices, and cerebellum. Adding cognitive metrics to the risk factor revealed the highest cumulative degree of connectivity for the pericalcarine cortex, insula, banks of the superior sulcus, and the cerebellum. To enable scaling up our approach, we extended tensor network principal component analysis, introducing CCA components. We developed sparse regression predictive models with errors of 17% for genotype, 24% for family risk factor for AD, and 5 years for age. Age prediction in groups including cognitively impaired subjects revealed regions not found using only normal subjects, i.e. middle and transverse temporal, paracentral and superior banks of temporal sulcus, as well as the amygdala and parahippocampal gyrus. These modeling approaches represent stepping stones towards single subject prediction.

Remote Sensing Low Signal-to-Noise-Ratio Target Detection Enhancement.

In real-time remote sensing application, frames of data are continuously flowing into the processing system. The capability of detecting objects of interest and tracking them as they move is crucial to many critical surveillance and monitoring missions. Detecting small objects using remote sensors is an ongoing, challenging problem. Since object(s) are located far away from the sensor, the target's Signal-to-Noise-Ratio (SNR) is low. The Limit of Detection (LOD) for remote sensors is bounded by what is observable on each image frame. In this paper, we present a new method, a "Multi-frame Moving Object Detection System (MMODS)", to detect small, low SNR objects that are beyond what a human can observe in a single video frame. This is demonstrated by using simulated data where our technology-detected objects are as small as one pixel with a targeted SNR, close to 1:1. We also demonstrate a similar improvement using live data collected with a remote camera. The MMODS technology fills a major technology gap in remote sensing surveillance applications for small target detection. Our method does not require prior knowledge about the environment, pre-labeled targets, or training data to effectively detect and track slow- and fast-moving targets, regardless of the size or the distance.

A multicontrast MR atlas of the Wistar rat brain.

We describe a multi-contrast, multi-dimensional atlas of the Wistar rat acquired at microscopic spatial resolution using magnetic resonance histology (MRH). Diffusion weighted images, and associated scalar images were acquired of a single specimen with a fully sampled Fourier reconstruction, 61 angles and b=3000 s/mm2 yielding 50 um isotropic spatial resolution. The higher angular sampling allows use of the GQI algorithm improving the angular invariance of the scalar images and yielding an orientation distribution function to assist in delineating subtle boundaries where there are crossing fibers  and track density images providing insight into local fiber architecture.  A multigradient echo image of the same specimen was acquired at 25 um isotropic spatial resolution. A quantitative susceptibility map enhances fiber architecture relative to the magnitude images.  An accompanying multi-specimen atlas (n=6) was acquired with compressed sensing with the same diffusion protocol as used for the single specimen atlas.  An average was created using diffeomorphic mapping. Scalar volumes from the diffusion data, a T2* weighted volume, a quantitative susceptibility map, and a track density volume, all registered to the same space provide multiple contrasts to assist in anatomic delineation. The new template  provides significantly increased contrast in the scalar DTI images when compared to previous atlases. A compact interactive viewer based on 3D Slicer is provided to facilitate comparison among the contrasts in the multiple volumes. The single volume and average atlas with multiple 3D volumes provide an improved template for anatomic interrogation of the Wistar rat brain. The improved contrast to noise in the scalar DTI images and the addition of other volumes (eg. QA,QSM,TDI ) will facilitate automated label registration for MR histology and preclinical imaging.

Variability and heritability of mouse brain structure: Microscopic MRI atlases and connectomes for diverse strains.

Genome-wide association studies have demonstrated significant links between human brain structure and common DNA variants. Similar studies with rodents have been challenging because of smaller brain volumes. Using high field MRI (9.4 T) and compressed sensing, we have achieved microscopic resolution and sufficiently high throughput for rodent population studies. We generated whole brain structural MRI and diffusion connectomes for four diverse isogenic lines of mice (C57BL/6J, DBA/2J, CAST/EiJ, and BTBR) at spatial resolution 20,000 times higher than human connectomes. We measured narrow sense heritability (h2) I.e. the fraction of variance explained by strains in a simple ANOVA model for volumes and scalar diffusion metrics, and estimates of residual technical error for 166 regions in each hemisphere and connectivity between the regions. Volumes of discrete brain regions had the highest mean heritability (0.71 ± 0.23 SD, n = 332), followed by fractional anisotropy (0.54 ± 0.26), radial diffusivity (0.34 ± 0.022), and axial diffusivity (0.28 ± 0.19). Connection profiles were statistically different in 280 of 322 nodes across all four strains. Nearly 150 of the connection profiles were statistically different between the C57BL/6J, DBA/2J, and CAST/EiJ lines. Microscopic whole brain MRI/DTI has allowed us to identify significant heritable phenotypes in brain volume, scalar DTI metrics, and quantitative connectomes.

Qualitative and Quantitative Neuropathology Approaches Using Magnetic Resonance Microscopy (Diffusion Tensor Imaging) and Stereology in a Hexachlorophene Model of Myelinopathy in Sprague-Dawley Rats.

It is well established that hexachlorophene, which is used as an antibacterial agent, causes intramyelinic edema in humans and animal models. The hexachlorophene myelinopathy model, in which male Sprague-Dawley rats received 25 to 30 mg/kg hexachlorophene by gavage for up to 5 days, provided an opportunity to compare traditional neuropathology evaluations with magnetic resonance microscopy (MRM) findings. In addition, stereology assessments of 3 neuroanatomical sites were compared to quantitative measurements of similar structures by MRM. There were positive correlations between hematoxylin and eosin and luxol fast blue stains and MRM for identifying intramyelinic edema in the cingulum of corpus callosum, optic chiasm, anterior commissure (aca), lateral olfactory tracts, pyramidal tracts (py), and white matter tracts in the cerebellum. Stereology assessments were focused on the aca, longitudinal fasciculus of the pons, and py and demonstrated differences between control and treated rats, as was observed using MRM. The added value of MRM assessments was the ability to acquire qualitative 3-dimensional (3-D) images and obtain quantitative measurements of intramyelinic edema in 26 neuroanatomical sites in the intact brain. Also, diffusion tensor imaging (fractional anisotropy [FA]) indicated that there were changes in the cytoarchitecture of the white matter as detected by decreases in the FA in the treated compared to the control rats. This study demonstrates creative strategies that are possible using qualitative and quantitative assessments of potential white matter neurotoxicants in nonclinical toxicity studies. Our results lead us to the conclusion that volumetric analysis by MRM and stereology adds significant value to the standard 2-D microscopic evaluations.

Four-component full configuration interaction quantum Monte Carlo for relativistic correlated electron problems.

An adaptation of the full configuration interaction quantum Monte Carlo (FCIQMC) method is presented for correlated electron problems containing heavy elements and the presence of significant relativistic effects. The modified algorithm allows for the sampling of the four-component spinors of the Dirac-Coulomb(-Breit) Hamiltonian within the relativistic no-pair approximation. The loss of spin symmetry and the general requirement for complex-valued Hamiltonian matrix elements are the most immediate considerations in expanding the scope of FCIQMC into the relativistic domain, and the alternatives for their efficient implementation are motivated and demonstrated. For the canonical correlated four-component chemical benchmark application of thallium hydride, we show that the necessary modifications do not particularly adversely affect the convergence of the systematic (initiator) error to the exact correlation energy for FCIQMC calculations, which is primarily dictated by the sparsity of the wavefunction, allowing the computational effort to somewhat bypass the formal increases in Hilbert space dimension for these problems. We apply the method to the larger problem of the spectroscopic constants of tin oxide, correlating 28 electrons in 122 Kramers-paired spinors, finding good agreement with experimental and prior theoretical relativistic studies.

Efficient and stochastic multireference perturbation theory for large active spaces within a full configuration interaction quantum Monte Carlo framework.

Full Configuration Interaction Quantum Monte Carlo (FCIQMC) has been effectively applied to very large configuration interaction (CI) problems and was recently adapted for use as an active space solver and combined with orbital optimization. In this work, we detail an approach within FCIQMC to allow for efficient sampling of fully internally contracted multireference perturbation theories within the same stochastic framework. Schemes are described to allow for the close control over the resolution of stochastic sampling of the effective higher-body intermediates within the active space. It is found that while complete active space second-order perturbation theory seems less amenable to a stochastic reformulation, strongly contracted N-Electron Valence second-order Perturbation Theory (NEVPT2) is far more stable, requiring a similar number of walkers to converge the sc-NEVPT2 expectation values as to converge the underlying CI problem. We demonstrate the application of the stochastic approach to the computation of sc-NEVPT2 within a (24, 24) active space in a biologically relevant system and show that small numbers of walkers are sufficient for a faithful sampling of the sc-NEVPT2 energy to chemical accuracy, despite the active space already exceeding the limits of practicality for traditional approaches. This raises prospects of an efficient stochastic solver for multireference chemical problems requiring large active spaces, with an accurate treatment of external orbitals.

NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods.

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24 000 central processing unit cores. In this paper, we describe the core functionalities of NECI and its recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing us to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs, and it is licensed under GPL-3.0.

Global patterns of kelp forest change over the past half-century.

Kelp forests (Order Laminariales) form key biogenic habitats in coastal regions of temperate and Arctic seas worldwide, providing ecosystem services valued in the range of billions of dollars annually. Although local evidence suggests that kelp forests are increasingly threatened by a variety of stressors, no comprehensive global analysis of change in kelp abundances currently exists. Here, we build and analyze a global database of kelp time series spanning the past half-century to assess regional and global trends in kelp abundances. We detected a high degree of geographic variation in trends, with regional variability in the direction and magnitude of change far exceeding a small global average decline (instantaneous rate of change = -0.018 y-1). Our analysis identified declines in 38% of ecoregions for which there are data (-0.015 to -0.18 y-1), increases in 27% of ecoregions (0.015 to 0.11 y-1), and no detectable change in 35% of ecoregions. These spatially variable trajectories reflected regional differences in the drivers of change, uncertainty in some regions owing to poor spatial and temporal data coverage, and the dynamic nature of kelp populations. We conclude that although global drivers could be affecting kelp forests at multiple scales, local stressors and regional variation in the effects of these drivers dominate kelp dynamics, in contrast to many other marine and terrestrial foundation species.

PITUITARY EVALUATION IN PATIENTS WITH LOW PROSTATE-SPECIFIC ANTIGEN.

OBJECTIVE: To evaluate pituitary function in men with a low screening prostate-specific antigen (PSA) of ≤0.1 ng/mL and test the hypothesis that low PSA is associated with hypogonadism alone or other hormone deficiency. METHODS: This was a case-control study evaluating the rates of hypogonadism and low insulin-like growth factor (IGF)-1 in a cohort of men with low or normal screening PSA level. Sixty-four men >40 years old without known prostate disease were divided into a low-PSA group (PSA ≤0.1 ng/mL) and normal-PSA group (PSA 1 to 4 ng/mL). Hormonal evaluation included total testosterone, prolactin, luteinizing hormone, follicle-stimulating hormone, IGF-1, growth hormone, thyroid-stimulating hormone, free thyroxine, morning cortisol, and adrenocorticotropic hormone. The difference between each patient's observed IGF-1 and the IGF-1 age-specific lower limit was calculated. The odds ratios (ORs) for having hypogonadism and associated 95% confidence intervals (CIs) were calculated using the Cochran-Mantel-Haenszel test. RESULTS: The rate of hypogonadism was significantly higher in the low-PSA group (n = 44) compared with the normal-PSA control group (n = 20) (45.5% vs. 15.0%; OR, 4.7; 95% CI, 1.2 to 18.4; P = .027). The total testosterone in the low-PSA group was significantly lower compared with the control group (181.7 ng/dL vs. 263.7 ng/dL; P = .008). IGF-1 values were below their lower bound in 18.6% of subjects in the low-PSA group, compared with 0% in the control group. CONCLUSION: Men with low PSA have significantly higher rates of hypogonadism and low IGF-1 compared with those with normal PSA. In such men, we recommend hormonal evaluation to exclude associated pituitary dysfunction. ABBREVIATIONS: BMI = body mass index; GH = growth hormone; IGF-1 = insulin-like growth factor 1; MRI = magnetic resonance imaging; PSA = prostate-specific antigen; T2DM = type 2 diabetes mellitus; VA-NWIHCS = VA-Nebraska Western Iowa Health Care System.

A phylogeographic investigation of the kelp genus Laminaria (Laminariales, Phaeophyceae), with emphasis on the South Atlantic Ocean.

The genus Laminaria has a wide distribution range compared with other kelp genera because it is found in both the North and the South Atlantic, on both sides of the North Pacific, as well as in the Mediterranean. Hypotheses behind this biogeographical pattern have been discussed by several authors but have not yet been fully evaluated with time-calibrated phylogenies. Based on the analysis of four molecular markers (ITS2, rbcL, atp8 and trnWI), our goal was to reassess the Laminaria species diversity in South Africa, assess its relationship with the other species distributed in the South Atlantic and reconstruct the historical biogeography of the genus. Our results confirm the occurrence of a single species, L. pallida, in southern Africa, and its sister relationship with the North Atlantic L. ochroleuca. Both species belonged to a clade containing the other South Atlantic species: L. abyssalis from Brazil, and the Mediterranean L. rodriguezii. Our time-calibrated phylogenies suggest that Laminaria originated in the northern Pacific around 25 mya, followed by at least two migration events through the Bering Strait after its opening (~5.32 mya). Today, the first is represented by L. solidungula in the Arctic, while the second gave rise to the rest of the Atlantic species. The colonization of the North Atlantic was followed by a gradual colonization southward along the west coast of Europe, into the Mediterranean (~2.07 mya) and two recent, but disconnected, migrations (~1.34 and 0.87 mya) across the equator, giving rise to L. abyssalis in Brazil and L. pallida in southern Africa, respectively.

Molecular systematics reveals increased diversity within the South African Laurencia complex (Rhodomelaceae, Rhodophyta).

Previous publications list ten species in the Laurencia complex from South Africa with all ascribed to the genus Laurencia sensu stricto. However, the diversity of the complex in South Africa has not yet been re-assessed following the numerous recent taxonomic changes. This study investigated the phylogenetic relationships and taxonomy of this group in South Africa using recent collections. Methods included molecular phylogenetic analyses of plastid rbcL gene sequences (a total of 146; including eleven outgroup taxa) using Maximum Likelihood and Bayesian Inference, and the examination of morphological and anatomical characters, including the number of corps en cerise when present. The seven genera of the Laurencia complex formed monophyletic clades with high posterior probabilities. Seventeen morphotypes were identified: 14 in the genus Laurencia sensu stricto, among which eight corresponded to Laurencia species currently recognized from South Africa and one each to species of Palisada, Chondrophycus, and Laurenciella. The six remaining morphotypes in Laurencia sensu stricto did not match any descriptions and are described here as five new species: Laurencia alfredensis sp. nov., Laurencia dichotoma sp. nov., Laurencia digitata sp. nov., Laurencia multiclavata sp. nov. and Laurencia sodwaniensis sp. nov. and a new variety: Laurencia pumila var. dehoopiensis var. nov. Laurencia stegengae nom. nov. is established to replace Laurencia peninsularis Stegenga, Bolton and Anderson nom. illeg. The diversity is likely greater, with six additional unidentified specimens found in this molecular investigation. These findings place South Africa alongside Australia in having one of the most diverse floras of this group in the world.

USE OF A CARBONIC ANHYDRASE INHIBITOR IN X-LINKED RETINOSCHISIS: Effect on Cystic-Appearing Macular Lesions and Visual Acuity.

PURPOSE: To evaluate changes in cystic-appearing macular lesions and visual acuity in patients with X-linked retinoschisis while being treated with a carbonic anhydrase inhibitor. METHODS: A retrospective analysis of 68 eyes from 36 patients between the ages of 5 years and 61 years with X-linked retinoschisis were monitored while on a carbonic anhydrase inhibitor. Macular cystic-appearing lesions were monitored with optical coherence tomography. Snellen visual acuity measurements were converted to logarithm of the minimum angle of resolution equivalent Early Treatment Diabetic Retinopathy Study letters for analysis. Analyses for changes in both visual acuity and macular cysts included comparisons between treatment and pretreatment segments. RESULTS: Forty-five eyes (66%) had a reduction of their cysts while on a carbonic anhydrase inhibitor. Twenty eyes (29%) showed no cystic change, whereas 3 eyes (4%) demonstrated worsening of their cysts with treatment when compared with pretreatment. There was a statistically significant improvement in logarithm of the minimum angle of resolution visual acuity while on treatment relative to pretreatment (P < 0.0001). The estimated average Early Treatment Diabetic Retinopathy Study equivalent improvement was 0.09 (slightly less than one line on the Early Treatment Diabetic Retinopathy Study chart) with a 95% confidence interval of 0.08 to 0.11. CONCLUSION: Considering the entire 36 patients in this cohort, while statistically significant, the average improvement in visual acuity was modest. Nonetheless, in individual patients, the improvement was more substantial. Improvement in the extent of cystic macular lesions was observed in a high percentage of cases.

The fornix provides multiple biomarkers to characterize circuit disruption in a mouse model of Alzheimer's disease.

Multivariate biomarkers are needed for detecting Alzheimer's disease (AD), understanding its etiology, and quantifying the effect of therapies. Mouse models provide opportunities to study characteristics of AD in well-controlled environments that can help facilitate development of early interventions. The CVN-AD mouse model replicates multiple AD hallmark pathologies, and we identified multivariate biomarkers characterizing a brain circuit disruption predictive of cognitive decline. In vivo and ex vivo magnetic resonance imaging (MRI) revealed that CVN-AD mice replicate the hippocampal atrophy (6%), characteristic of humans with AD, and also present changes in subcortical areas. The largest effect was in the fornix (23% smaller), which connects the septum, hippocampus, and hypothalamus. In characterizing the fornix with diffusion tensor imaging, fractional anisotropy was most sensitive (20% reduction), followed by radial (15%) and axial diffusivity (2%), in detecting pathological changes. These findings were strengthened by optical microscopy and ultrastructural analyses. Ultrastructual analysis provided estimates of axonal density, diameters, and myelination-through the g-ratio, defined as the ratio between the axonal diameter, and the diameter of the axon plus the myelin sheath. The fornix had reduced axonal density (47% fewer), axonal degeneration (13% larger axons), and abnormal myelination (1.5% smaller g-ratios). CD68 staining showed that white matter pathology could be secondary to neuronal degeneration, or due to direct microglial attack. In conclusion, these findings strengthen the hypothesis that the fornix plays a role in AD, and can be used as a disease biomarker and as a target for therapy.

Canopy interactions and physical stress gradients in subtidal communities.

Species interactions are integral drivers of community structure and can change from competitive to facilitative with increasing environmental stress. In subtidal marine ecosystems, however, interactions along physical stress gradients have seldom been tested. We observed seaweed canopy interactions across depth and latitudinal gradients to test whether light and temperature stress structured interaction patterns. We also quantified interspecific and intraspecific interactions among nine subtidal canopy seaweed species across three continents to examine the general nature of interactions in subtidal systems under low consumer pressure. We reveal that positive and neutral interactions are widespread throughout global seaweed communities and the nature of interactions can change from competitive to facilitative with increasing light stress in shallow marine systems. These findings provide support for the stress gradient hypothesis within subtidal seaweed communities and highlight the importance of canopy interactions for the maintenance of subtidal marine habitats experiencing environmental stress.

Update on blood pressure control and renal outcomes in diabetes mellitus.

Diabetic nephropathy (DN) is associated with a high incidence of cardiovascular (CV) morbidity and mortality. Although relationships between hypertension and diabetic nephropathy are complex, blood pressure (BP) control is an important management strategy in the prevention of onset and progression of DN in patients with diabetes mellitus (DM). Recent guidelines recommend less stringent BP targets among patients with type 2 DM and chronic kidney disease. These recommendations are based mostly on lack of benefit in CV outcomes with a low BP target. We review the current information on efficacy of BP control in improving renal outcomes in patients with type 2 DM. Presently, although intensive BP control has been was beneficial in decreasing albuminuria, it has not translated into reductions in risks of hard renal endpoints, such as progression to end-stage renal disease, the need for renal replacement therapy, and mortality from renal causes.

A new research strategy for integrating studies of paleoclimate, paleoenvironment, and paleoanthropology.

Paleoanthropologists (scientists studying human origins) universally recognize the evolutionary significance of ancient climates and environments for understanding human origins. Even those scientists working in recent phases of human evolution, when modern humans evolved, agree that hunter-gatherer adaptations are tied to the way that climate and environment shape the food and technological resource base. The result is a long tradition of paleoanthropologists engaging with climate and environmental scientists in an effort to understand if and how hominin bio-behavioral evolution responded to climate and environmental change. Despite this unusual consonance, the anticipated rewards of this synergy are unrealized and, in our opinion, will not reach potential until there are some fundamental changes in the way the research model is constructed. Discovering the relation between climate and environmental change to human origins must be grounded in a theoretical framework and a causal understanding of the connection between climate, environment, resource patterning, behavior, and morphology, then move beyond the strict correlative research that continues to dominate the field.

A molecular investigation of the genus Ecklonia (Phaeophyceae, Laminariales) with special focus on the Southern Hemisphere.

Brown algae of the order Laminariales, commonly referred to as kelps, are the largest and most productive primary producers in the coastal inshore environment. The genus Ecklonia (Lessoniaceae, Phaeophyceae) consists of seven species with four species in the Northern Hemisphere and three in the Southern Hemisphere. It was recently transferred to the family Lessoniaceae based on phylogenetic analyses of nuclear and chloroplastic markers, though the type of the genus was not included and its relationship with allied genera Eckloniopsis and Eisenia remained unresolved. The present study is the first to produce a phylogeny focussed on the genus Ecklonia. It included sequences from nuclear, mitochondrial, and chloroplastic DNA, for most of the distribution range of the three current Southern Hemisphere species (Ecklonia radiata, Ecklonia maxima, and a sample of a putative Ecklonia brevipes specimen), sequences for East Asiatic species (Ecklonia cava, Ecklonia kurome, and Ecklonia stolonifera), as well as the closely related genera Eckloniopsis and Eisenia. Results confirmed E. radiata and E. maxima as two distinct species in South Africa, E. radiata as a single species throughout the Southern Hemisphere (in South Africa, Australia, and New Zealand) and East Asiatic species as a distinct lineage from the Southern Hemisphere clade. Results further pointed out a close sister relationship between Eckloniopsis radicosa and two Eisenia species (including the type species: Eisenia arborea) to the genus Ecklonia suggesting that the genera Eckloniopsis and Eisenia are superfluous.

Simultaneous multislice spectral-spatial excitations for reduced signal loss susceptibility artifact in BOLD functional MRI.

PURPOSE: Simultaneous multislice (SMS) imaging can significantly increase image acquisition rates and improve temporal resolution and contrast in gradient-echo blood oxygen level-dependent (BOLD) functional MRI (fMRI) experiments. Through-plane signal loss due to B(0) inhomogeneities at air-tissue interfaces limits fMRI of structures near the nasal cavity and ear canals. This study implemented spectral-spatial (SPSP) radiofrequency pulses for reduced through-plane signal loss across multiple simultaneously excited slices. THEORY AND METHODS: Multiband (MB) and power independent of number of slices (PINS) methods are combined with SPSP excitation for signal loss compensation in slice-accelerated human brain imaging. Nine simultaneous slices of 5-mm thickness and 20 mm apart were excited using standard MB radiofrequency pulses and the proposed SPSP-SMS pulses, yielding coverage of 36 slices in four shots with 350-ms volume pulse repetition time. The pulses were compared in breath-hold fMRI at 3T. RESULTS: The SPSP-SMS pulses recovered ∼45% of voxels with signal loss in standard SMS images. Activation in areas of signal recovery increased by 26.4% using a 12.6-ms SPSP-MB pulse and 20.3% using a 12.1-ms SPSP-PINS pulse. CONCLUSIONS: It is demonstrated that SPSP-SMS pulses can improve BOLD sensitivity in areas of signal loss across simultaneous multiple slices.