Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

Pulse wave analysis of aortic pressure: diastole should also be considered.

BACKGROUND: : The systolic augmentation index (sAix), calculated from the central aortic pulse wave (reconstructed from the noninvasive recording of the radial pulse with applanation tonometry), is widely used as a simple index of central arterial stiffness, but has the disadvantage of also being influenced by the timing of the reflected with respect to the forward pressure wave, as shown by its inverse dependence on heart rate (HR). During diastole, the central aortic pulse also contains reflected waves, but their relationship to arterial stiffness and HR has not been studied. METHODS: : In 48 men and 45 women, all healthy, with ages ranging from 19 to 70 years, we measured pulse wave velocity (PWV, patients supine), a standard evaluator of arterial stiffness, and carried out radial applanation tonometry (patients sitting and supine). The impact of reflected waves on the diastolic part of the aortic pressure waveform was quantified in the form of a diastolic augmentation index (dAix). RESULTS: : Across ages, sexes, and body position, there was an inverse relationship between the sAix and the dAix. When PWV and HR were added as covariates to a prediction model including age, sex and body position as main factors, the sAix was directly related to PWV (P < 0.0001) and inversely to HR (P < 0.0001). With the same analysis, the dAix was inversely related to PWV (P < 0.0001) and independent of HR (P = 0.52). CONCLUSION: : The dAix has the same degree of linkage to arterial stiffness as the more conventional sAix, while being immune to the confounding effect of HR. The quantification of diastolic aortic pressure augmentation by reflected waves could be a useful adjunct to pulse wave analysis.

Cellular distribution of glucose and monocarboxylate transporters in human brain white matter and multiple sclerosis lesions.

To ensure efficient energy supply to the high demanding brain, nutrients are transported into brain cells via specific glucose (GLUT) and monocarboxylate transporters (MCT). Mitochondrial dysfunction and altered glucose metabolism are thought to play an important role in the progression of neurodegenerative diseases, including multiple sclerosis (MS). Here, we investigated the cellular localization of key GLUT and MCT proteins in human brain tissue of non-neurological controls and MS patients. We show that in control brain tissue GLUT and MCT proteins were abundantly expressed in a variety of central nervous system cells, particularly in microglia and endothelial cells. In active MS lesions, GLUTs and MCTs were highly expressed in infiltrating leukocytes and reactive astrocytes. Astrocytes manifest increased MCT1 staining and maintain GLUT expression in inactive lesions, whereas demyelinated axons exhibit significantly reduced GLUT3 and MCT2 immunoreactivity in inactive lesions. Finally, we demonstrated that the co-transcription factor peroxisome proliferator-activated receptor gamma co-activator 1-alpha (PGC-1α), an important protein involved in energy metabolism, is highly expressed in reactive astrocytes in active MS lesions. Overexpression of PGC-1α in astrocyte-like cells resulted in increased production of several GLUT and MCT proteins. In conclusion, we provide for the first time a comprehensive overview of key nutrient transporters in white matter brain samples. Moreover, our data demonstrate an altered expression of these nutrient transporters in MS brain tissue, including a marked reduction of axonal GLUT3 and MCT2 expression in chronic lesions, which may impede efficient nutrient supply to the hypoxic demyelinated axons thereby contributing to the ongoing neurodegeneration in MS. GLIA 2014;62:1125-1141.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Twitch and M-wave potentiation induced by intermittent maximal voluntary quadriceps contractions: Differences between direct quadriceps and femoral nerve stimulation.

Introduction: To investigate differences in twitch and M-wave potentiation in the quadriceps femoris when electrical stimulation is applied over the quadriceps muscle belly versus the femoral nerve trunk. Methods: M-waves and mechanical twitches were evoked using direct quadriceps muscle and femoral nerve stimulation between 48 successive isometric maximal voluntary contractions (MVC) from 10 young, healthy subjects. Potentiation was investigated by analyzing the changes in M-wave amplitude recorded from the vastus medialis (VM) and vastus lateralis (VL) muscles and in quadriceps peak twitch force. Results: Potentiation of twitch, VM M-wave, and VL M-wave were greater for femoral nerve than for direct quadriceps stimulation (P<0.05). Despite a 50% decrease in MVC force, the amplitude of the M-waves increased significantly during exercise. Conclusions: In addition to enhanced electrogenic Na(+) -K(+) pumping, other factors (such as synchronization in activation of muscle fibers and muscle architectural properties) might significantly influence the magnitude of M-wave enlargement. © 2013 Wiley Periodicals, Inc.

Grey matter changes in motor conversion disorder.

OBJECTIVE: To detect anatomical differences in areas related to motor processing between patients with motor conversion disorder (CD) and controls. METHODS: T1-weighted 3T brain MRI data of 15 patients suffering from motor CD (nine with hemiparesis and six with paraparesis) and 25 age- and gender-matched healthy volunteers were compared using voxel-based morphometry (VBM) and voxel-based cortical thickness (VBCT) analysis. RESULTS: We report significant cortical thickness (VBCT) increases in the bilateral premotor cortex of hemiparetic patients relative to controls and a trend towards increased grey matter volume (VBM) in the same region. Regression analyses showed a non-significant positive correlation between cortical thickness changes and symptom severity as well as illness duration in CD patients. CONCLUSIONS: Cortical thickness increases in premotor cortical areas of patients with hemiparetic CD provide evidence for altered brain structure in a condition with presumed normal brain anatomy. These may either represent premorbid vulnerability or a plasticity phenomenon related to the disease with the trends towards correlations with clinical variables supporting the latter.

Proton T1 relaxation times of metabolites in human occipital white and gray matter at 7 T.

Proton T1 relaxation times of metabolites in the human brain have not previously been published at 7 T. In this study, T1 values of CH3 and CH2 group of N-acetylaspartate and total creatine as well as nine other brain metabolites were measured in occipital white matter and gray matter at 7 T using an inversion-recovery technique combined with a newly implemented semi-adiabatic spin-echo full-intensity acquired localized spectroscopy sequence (echo time = 12 ms). The mean T1 values of metabolites in occipital white matter and gray matter ranged from 0.9 to 2.2 s. Among them, the T1 of glutathione, scyllo-inositol, taurine, phosphorylethanolamine, and N-acetylaspartylglutamate were determined for the first time in the human brain. Significant differences in T1 between white matter and gray matter were found for water (-28%), total choline (-14%), N-acetylaspartylglutamate (-29%), N-acetylaspartate (+4%), and glutamate (+8%). An increasing trend in T1 was observed when compared with previously reported values of N-acetylaspartate (CH3 ), total creatine (CH3 ), and total choline at 3 T. However, for N-acetylaspartate (CH3 ), total creatine, and total choline, no substantial differences compared to previously reported values at 9.4 T were discernible. The T1 values reported here will be useful for the quantification of metabolites and signal-to-noise optimization in human brain at 7 T. Magn Reson Med 69:931-936, 2013. © 2012 Wiley Periodicals, Inc.

Structural covariance of superficial white matter in mild Alzheimer's disease compared to normal aging.

INTRODUCTION: Interindividual variations in regional structural properties covary across the brain, thus forming networks that change as a result of aging and accompanying neurological conditions. The alterations of superficial white matter (SWM) in Alzheimer's disease (AD) are of special interest, since they follow the AD-specific pattern characterized by the strongest neurodegeneration of the medial temporal lobe and association cortices. METHODS: Here, we present an SWM network analysis in comparison with SWM topography based on the myelin content quantified with magnetization transfer ratio (MTR) for 39 areas in each hemisphere in 15 AD patients and 15 controls. The networks are represented by graphs, in which nodes correspond to the areas, and edges denote statistical associations between them. RESULTS: In both groups, the networks were characterized by asymmetrically distributed edges (predominantly in the left hemisphere). The AD-related differences were also leftward. The edges lost due to AD tended to connect nodes in the temporal lobe to other lobes or nodes within or between the latter lobes. The newly gained edges were mostly confined to the temporal and paralimbic regions, which manifest demyelination of SWM already in mild AD. CONCLUSION: This pattern suggests that the AD pathological process coordinates SWM demyelination in the temporal and paralimbic regions, but not elsewhere. A comparison of the MTR maps with MTR-based networks shows that although, in general, the changes in network architecture in AD recapitulate the topography of (de)myelination, some aspects of structural covariance (including the interhemispheric asymmetry of networks) have no immediate reflection in the myelination pattern.

B1-insensitive T2 preparation for improved coronary magnetic resonance angiography at 3 T.

At 3 T, the effective wavelength of the RF field is comparable to the dimension of the human body, resulting in B1 standing wave effects and extra variations in phase. This effect is accompanied by an increase in B0 field inhomogeneity compared to 1.5 T. This combination results in nonuniform magnetization preparation by the composite MLEV weighted T2 preparation (T2 Prep) sequence used for coronary magnetic resonance angiography (MRA). A new adiabatic refocusing T2 Prep sequence is presented in which the magnetization is tipped into the transverse plane with a hard RF pulse and refocused using a pair of adiabatic fast-passage RF pulses. The isochromats are subsequently returned to the longitudinal axis using a hard RF pulse. Numerical simulations predict an excellent suppression of artifacts originating from B1 inhomogeneity while achieving good contrast enhancement between coronary arteries and surrounding tissue. This was confirmed by an in vivo study, in which coronary MR angiograms were obtained without a T2 Prep, with an MLEV weighted T2 Prep and the proposed adiabatic T2 Prep. Improved quantitative and qualitative coronary MRA image measurement was achieved using the adiabatic T2 Prep at 3 T.

Prediction of functional outcome 18 months after a first psychotic episode: a proton magnetic resonance spectroscopy study.

CONTEXT: Recent magnetic resonance imaging studies have attempted to relate volumetric brain measurements in early schizophrenia to clinical and functional outcome some years later. These studies have generally been negative, perhaps because gray and white matter volumes inaccurately assess the underlying dysfunction that might be predictive of outcome. OBJECTIVE: To investigate the predictive value of frontal and temporal spectroscopy measures for outcome in patients with first-episode psychoses. DESIGN: Left prefrontal cortex and left mediotemporal lobe voxels were assessed using proton magnetic resonance spectroscopy to provide the ratio of N-acetylaspartate (NAA) and choline-containing compounds to creatine and phosphocreatine (Cr) (NAA/Cr ratio). These data were used to predict outcome at 18 months after admission, as assessed by a systematic medical record audit. SETTING: Early psychosis clinic. PARTICIPANTS: Forty-six patients with first-episode psychosis. MAIN OUTCOME MEASURES: We used regression models that included age at imaging and duration of untreated psychosis to predict outcome scores on the Global Assessment of Functioning Scale, Clinical Global Impression scales, and Social and Occupational Functional Assessment Scale, as well as the number of admissions during the treatment period. We then further considered the contributions of premorbid function and baseline level of negative symptoms. RESULTS: The only spectroscopic predictor of outcome was the NAA/Cr ratio in the prefrontal cortex. Low scores on this variable were related to poorer outcome on all measures. In addition, the frontal NAA/Cr ratio explained 17% to 30% of the variance in outcome. CONCLUSIONS: Prefrontal neuronal dysfunction is an inconsistent feature of early psychosis; rather, it is an early marker of poor prognosis across the first years of illness. The extent to which this can be used to guide treatment and whether it predicts outcome some years after first presentation are questions for further research.