A BAYESIAN HIERARCHICAL FRAMEWORK FOR SPATIAL MODELING OF fMRI DATA

Functional neuroimaging techniques enable investigations into the neural basis of human cognition, emotions, and behaviors. In practice, applications of functional magnetic resonance imaging (fMRI) have provided novel insights into the neuropathophysiology of major psychiatric,neurological, and substance abuse disorders, as well as into the neural responses to their treatments. Modern activation studies often compare localized task-induced changes in brain activity between experimental groups. One may also extend voxel-level analyses by simultaneously considering the ensemble of voxels constituting an anatomically defined region of interest (ROI) or by considering means or quantiles of the ROI. In this work we present a Bayesian extension of voxel-level analyses that offers several notable benefits. First, it combines whole-brain voxel-by-voxel modeling and ROI analyses within a unified framework. Secondly, an unstructured variance/covariance for regional mean parameters allows for the study of inter-regional functional connectivity, provided enough subjects are available to allow for accurate estimation. Finally, an exchangeable correlation structure within regions allows for the consideration of intra-regional functional connectivity. We perform estimation for our model using Markov Chain Monte Carlo (MCMC) techniques implemented via Gibbs sampling which, despite the high throughput nature of the data, can be executed quickly (less than 30 minutes). We apply our Bayesian hierarchical model to two novel fMRI data sets: one considering inhibitory control in cocaine-dependent men and the second considering verbal memory in subjects at high risk for Alzheimer’s disease. The unifying hierarchical model presented in this manuscript is shown to enhance the interpretation content of these data sets.

Functional relevance of DNA polymorphisms within the promoter region of the prion protein gene and their association to BSE infection

Transmissible spongiform encephalopathies (TSEs) are a group of neurodegenerative diseases that can occur spontaneously or can be caused by infection or mutations within the prion protein gene PRNP. Nonsynonymous DNA polymorphisms within the PRNP gene have been shown to influence susceptibility/resistance to infection in sheep and humans. Analysis of DNA polymorphisms within the core promoter region of the PRNP gene in four major German bovine breeds resulted in the identification of both SNPs and insertion/deletion (indel) polymorphisms. Comparative genotyping of both controls and animals that tested positive for bovine spongiform encephalopathy (BSE) revealed a significantly different distribution of two indel polymorphisms and two SNPs within Braunvieh animals, suggesting an association of these polymorphisms with BSE susceptibility. The functional relevance of these polymorphisms was analyzed using reporter gene constructs in neuronal cells. A specific haplotype near exon 1 was identified that exhibited a significantly lower expression level. Genotyping of nine polymorphisms within the promoter region and haplotype calculation revealed that the haplotype associated with the lowest expression level was underrepresented in the BSE group of all breeds compared to control animals, indicating a correlation of reduced PRNP expression and increased resistance to BSE.

Short communication: duplication in the 5'-flanking region of the beta-lactoglobulin gene is linked to the BLG A allele

beta-Lactoglobulin (beta-LG) is the major whey protein in the milk of cows and other ruminants. It is well established that the predominant genetic variants beta-LG A and B are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants. However, the genetic basis for the differentially expressed BLG A and B alleles is still elusive. In this study additional genetic variation further upstream in the 5'-flanking region of the BLG gene was identified, including 6 single nucleotide substitutions, a single nucleotide deletion, and a 7-bp duplication. Comparison of DNA sequences showed that the investigated 5'-flanking region is highly conserved between ruminants, and the duplication g.-1885_-1879dupCTCTCGC and the substitution g.-1888A>G are only found in the BLG A and D alleles in cattle. The cytosine at position g.-1957 and the thymines at positions g.-2008 and g.-2049 are only found in BLG B alleles of cattle. It is suggested that the described genetic variability contributes to the differential allelic expression of the BLG gene.

A homologue of aliB is found in the capsule region of nonencapsulated Streptococcus pneumoniae

The epidemiology, phylogeny, and biology of nonencapsulated Streptococcus pneumoniae are largely unknown. Increased colonization capacity and transformability are, however, intriguing features of these pneumococci and play an important role. Twenty-seven nonencapsulated pneumococci were identified in a nationwide collection of 1,980 nasopharyngeal samples and 215 blood samples obtained between 1998 and 2002. On the basis of multilocus sequence typing and capsule region analysis we divided the nonencapsulated pneumococci into two groups. Group I was closely related to encapsulated strains. Group II had a clonal population structure, including two geographically widespread clones able to cause epidemic conjunctivitis and invasive diseases. Group II strains also carried a 1,959-bp homologue of aliB (aliB-like ORF 2) in the capsule region, which was highly homologous to a sequence in the capsule region of Streptococcus mitis. In addition, strains of the two major clones in group II had an additional sequence, aliB-like ORF 1 (1,968 to 2,004 bp), upstream of aliB-like ORF 2. Expression of aliB-like ORF 1 was detected by reverse transcription-PCR, and the corresponding RNA was visualized by Northern blotting. A gene fragment homologous to capN of serotypes 33 and 37 suggests that group II strains were derived from encapsulated pneumococci some time ago. Therefore, loss of capsule expression in vivo was found to be associated with the importation of one or two aliB homologues in some nonencapsulated pneumococci.

Adequacy of an early arterial phase low-volume contrast protocol in 64-detector computed tomography angiography for aortoiliac aneurysms

PURPOSE: This retrospective study was conducted to determine whether a low-volume contrast medium protocol provides sufficient enhancement for 64-detector computed tomography angiography (CTA) in patients with aortoiliac aneurysms. METHODS: Evaluated were 45 consecutive patients (6 women; mean age, 72 +/- 6 years) who were referred for aortoiliac computed tomography angiography between October 2005 and January 2007. Group A (22 patients; creatinine clearance, 64.2 +/- 8.1 mL/min) received 50 mL of the contrast agent. Group B (23 patients; creatinine clearance, 89.4 +/- 7.3 mL/min) received 100 mL of the contrast agent. The injection rate was 3.5 mL/s, followed by 30 mL of saline at 3.5 mL/s. Studies were performed on the same 64-detector computed tomography scanner using a real-time bolus-tracking technique. Quantitative analysis was performed by determination of mean vascular attenuation at 10 regions of interest from the suprarenal aorta to the common femoral artery by one reader blinded to type and amount of contrast agent and compared using the Student t test. Image quality according to a 4-point scale was assessed in consensus by two readers blinded to type and amount of contrast medium and compared using the Mann-Whitney test. Multivariable adjustments were performed using ordinal regression analysis. RESULTS: Mean total attenuation did not differ significantly between both groups (196.5 +/- 33.0 Hounsfield unit [HU] in group A and 203.1 +/- 44.2 HU in group B; P = .57 by univariate and P > .05 by multivariable analysis). Accordingly, attenuation at each region of interest was not significantly different (P > .35). Image quality was excellent or good in all patients. No significant differences in visual assessment were found comparing both contrast medium protocols (P > .05 by univariate and by multivariable analysis). CONCLUSIONS: Aortoiliac aneurysm imaging can be performed with substantially reduced amounts of contrast medium using 64-detector computed tomography angiography technology.

Measuring energy spectra of TeV gamma-ray emission from the Cygnus region of our galaxy with Milagro

High energy gamma rays can provide fundamental clues to the origins of cosmic rays. In this thesis, TeV gamma-ray emission from the Cygnus region is studied. Previously the Milagro experiment detected five TeV gamma-ray sources in this region and a significant excess of TeV gamma rays whose origin is still unclear. To better understand the diffuse excess the separation of sources and diffuse emission is studied using the latest and most sensitive data set of the Milagro experiment. In addition, a newly developed technique is applied that allows the energy spectrum of the TeV gamma rays to be reconstructed using Milagro data. No conclusive statement can be made about the spectrum of the diffuse emission from the Cygnus region because of its low significance of 2.2 σ above the background in the studied data sample. The entire Cygnus region emission is best fit with a power law with a spectral index of α=2.40 (68% confidence interval: 1.35-2.92) and a exponential cutoff energy of 31.6 TeV (10.0-251.2 TeV). In the case of a simple power law assumption without a cutoff energy the best fit yields a spectral index of α=2.97 (68% confidence interval: 2.83-3.10). Neither of these best fits are in good agreement with the data. The best spectral fit to the TeV emission from MGRO J2019+37, the brightest source in the Cygnus region, yields a spectral index of α=2.30 (68% confidence interval: 1.40-2.70) with a cutoff energy of 50.1 TeV (68% confidence interval: 17.8-251.2 TeV) and a spectral index of α=2.75 (68% confidence interval: 2.65-2.85) when no exponential cutoff energy is assumed. According to the present analysis, MGRO J2019+37 contributes 25% to the differential flux from the entire Cygnus at 15 TeV.

Performance evaluation of VIP latrines in the Upper West Region of Ghana

The prevalence of Ventilated Improved Pit (VIP) latrines in Ghana suggests that the design must have a high user acceptance. The two key factors attributed to user acceptance of a VIP latrine over an alternative latrine design, such as the basic pit latrine, are its ability to remove foul odors and maintain low fly populations; both of which are a direct result of an adequate ventilation flow rate. Adequate ventilation for odorless conditions in a VIP latrine has been defined by the United Nations Development Program (UNDP) and the World Bank, as an air flow rate equivalent to 6 air changes per hour (6 ACH) of the superstructure’s air volume. Additionally, the UNDP determined that the three primary factors that affect ventilation are: 1) wind passing over the mouth of the vent pipe, 2) wind passing into the superstructure, and 3) solar radiation on to the vent pipe. Previous studies also indicate that vent pipes with larger diameters increase flow rates, and the application of carbonaceous materials to the pit sludge reduces odor and insect prevalence. Furthermore, proper design and construction is critical for the correct functioning of VIP latrines. Under-designing could cause problems with odor and insect control; over-designing would increase costs unnecessarily, thereby making it potentially unaffordable for benefactors to independently construct, repair or replace a VIP latrine. The present study evaluated the design of VIP latrines used by rural communities in the Upper West Region of Ghana with the focus of assessing adequate ventilation for odor removal and insect control. Thirty VIP latrines from six communities in the Upper West Region of Ghana were sampled. Each VIP latrine’s ventilation flow rate and micro-environment was measured using a hot-wire anemometer probe and portable weather station for a minimum of four hours. To capture any temporal or seasonal variations in ventilation, ten of the latrines were sampled monthly over the course of three months for a minimum of 12 hours. A latrine usage survey and a cost analysis were also conducted to further assess the VIP latrine as an appropriated technology for sustainable development in the Upper West Region. It was found that the average air flow rate over the entire sample set was 11.3 m3/hr. The minimum and maximum air flow rates were 0.0 m3/hr and 48.0 m3/hr respectively. Only 1 of the 30 VIP latrines (3%) was found to have an air flow rate greater than the UNDP-defined odorless condition of 6 ACH. Furthermore, 19 VIP latrines (63%) were found to have an average air flow rate of less than half the flow rate required to achieve 6 ACH. The dominant factors affecting ventilation flow rate were wind passing over the mouth of the vent pipe and air buoyancy forces, which were the effect of differences in temperature between the substructure and the ambient environment. Of 76 usable VIP latrines found in one community, 68.4% were in actual use. The cost of a VIP latrine was found to be equivalent to approximately 12% of the mean annual household income for Upper West Region inhabitants.

Magnetization transfer contrast and T2 mapping in the evaluation of cartilage repair tissue with 3T MRI

PURPOSE: To use magnetization transfer (MT) imaging in the visualization of healthy articular cartilage and cartilage repair tissue after different cartilage repair procedures, and to assess global as well as zonal values and compare the results to T2-relaxation. MATERIALS AND METHODS: Thirty-four patients (17 after microfracture [MFX] and 17 after matrix-associated autologous cartilage transplantation [MACT]) were examined with 3T MRI. The MT ratio (MTR) was calculated from measurements with and without MT contrast. T2-values were evaluated using a multiecho, spin-echo approach. Global (full thickness of cartilage) and zonal (deep and superficial aspect) region-of-interest assessment of cartilage repair tissue and normal-appearing cartilage was performed. RESULTS: In patients after MFX and MACT, the global MTR of cartilage repair tissue was significantly lower compared to healthy cartilage. In contrast, using T2, cartilage repair tissue showed significantly lower T2 values only after MFX, whereas after MACT, global T2 values were comparable to healthy cartilage. For zonal evaluation, MTR and T2 showed a significant stratification within healthy cartilage, and T2 additionally within cartilage repair tissue after MACT. CONCLUSION: MT imaging is capable and sensitive in the detection of differences between healthy cartilage and areas of cartilage repair and might be an additional tool in biochemical cartilage imaging. For both MTR and T2 mapping, zonal assessment is desirable.

In vivo biochemical 7.0 Tesla magnetic resonance: preliminary results of dGEMRIC, zonal T2, and T2* mapping of articular cartilage

INTRODUCTION: Ultra-high-field whole-body systems (7.0 T) have a high potential for future human in vivo magnetic resonance imaging (MRI). In musculoskeletal MRI, biochemical imaging of articular cartilage may benefit, in particular. Delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) and T2 mapping have shown potential at 3.0 T. Although dGEMRIC, allows the determination of the glycosaminoglycan content of articular cartilage, T2 mapping is a promising tool for the evaluation of water and collagen content. In addition, the evaluation of zonal variation, based on tissue anisotropy, provides an indicator of the nature of cartilage ie, hyaline or hyaline-like articular cartilage.Thus, the aim of our study was to show the feasibility of in vivo dGEMRIC, and T2 and T2* relaxation measurements, at 7.0 T MRI; and to evaluate the potential of T2 and T2* measurements in an initial patient study after matrix-associated autologous chondrocyte transplantation (MACT) in the knee. MATERIALS AND METHODS: MRI was performed on a whole-body 7.0 T MR scanner using a dedicated circular polarization knee coil. The protocol consisted of an inversion recovery sequence for dGEMRIC, a multiecho spin-echo sequence for standard T2 mapping, a gradient-echo sequence for T2* mapping and a morphologic PD SPACE sequence. Twelve healthy volunteers (mean age, 26.7 +/- 3.4 years) and 4 patients (mean age, 38.0 +/- 14.0 years) were enrolled 29.5 +/- 15.1 months after MACT. For dGEMRIC, 5 healthy volunteers (mean age, 32.4 +/- 11.2 years) were included. T1 maps were calculated using a nonlinear, 2-parameter, least squares fit analysis. Using a region-of-interest analysis, mean cartilage relaxation rate was determined as T1 (0) for precontrast measurements and T1 (Gd) for postcontrast gadopentate dimeglumine [Gd-DTPA(2-)] measurements. T2 and T2* maps were obtained using a pixelwise, monoexponential, non-negative least squares fit analysis; region-of-interest analysis was carried out for deep and superficial cartilage aspects. Statistical evaluation was performed by analyses of variance. RESULTS: Mean T1 (dGEMRIC) values for healthy volunteers showed slightly different results for femoral [T1 (0): 1259 +/- 277 ms; T1 (Gd): 683 +/- 141 ms] compared with tibial cartilage [T1 (0): 1093 +/- 281 ms; T1 (Gd): 769 +/- 150 ms]. Global mean T2 relaxation for healthy volunteers showed comparable results for femoral (T2: 56.3 +/- 15.2 ms; T2*: 19.7 +/- 6.4 ms) and patellar (T2: 54.6 +/- 13.0 ms; T2*: 19.6 +/- 5.2 ms) cartilage, but lower values for tibial cartilage (T2: 43.6 +/- 8.5 ms; T2*: 16.6 +/- 5.6 ms). All healthy cartilage sites showed a significant increase from deep to superficial cartilage (P < 0.001). Within healthy cartilage sites in MACT patients, adequate values could be found for T2 (56.6 +/- 13.2 ms) and T2* (18.6 +/- 5.3 ms), which also showed a significant stratification. Within cartilage repair tissue, global mean values showed no difference, with 55.9 +/- 4.9 ms for T2 and 16.2 +/- 6.3 ms for T2*. However, zonal assessment showed only a slight and not significant increase from deep to superficial cartilage (T2: P = 0.174; T2*: P = 0.150). CONCLUSION: In vivo T1 dGEMRIC assessment in healthy cartilage, and T2 and T2* mapping in healthy and reparative articular cartilage, seems to be possible at 7.0 T MRI. For T2 and T2*, zonal variation of articular cartilage could also be evaluated at 7.0 T. This zonal assessment of deep and superficial cartilage aspects shows promising results for the differentiation of healthy and affected articular cartilage. In future studies, optimized protocol selection, and sophisticated coil technology, together with increased signal at ultra-high-field MRI, may lead to advanced biochemical cartilage imaging.

High-resolution morphological and biochemical imaging of articular cartilage of the ankle joint at 3.0 T using a new dedicated phased array coil: in vivo reproducibility study

OBJECTIVE: The objective of this study was to evaluate the feasibility and reproducibility of high-resolution magnetic resonance imaging (MRI) and quantitative T2 mapping of the talocrural cartilage within a clinically applicable scan time using a new dedicated ankle coil and high-field MRI. MATERIALS AND METHODS: Ten healthy volunteers (mean age 32.4 years) underwent MRI of the ankle. As morphological sequences, proton density fat-suppressed turbo spin echo (PD-FS-TSE), as a reference, was compared with 3D true fast imaging with steady-state precession (TrueFISP). Furthermore, biochemical quantitative T2 imaging was prepared using a multi-echo spin-echo T2 approach. Data analysis was performed three times each by three different observers on sagittal slices, planned on the isotropic 3D-TrueFISP; as a morphological parameter, cartilage thickness was assessed and for T2 relaxation times, region-of-interest (ROI) evaluation was done. Reproducibility was determined as a coefficient of variation (CV) for each volunteer; averaged as root mean square (RMSA) given as a percentage; statistical evaluation was done using analysis of variance. RESULTS: Cartilage thickness of the talocrural joint showed significantly higher values for the 3D-TrueFISP (ranging from 1.07 to 1.14 mm) compared with the PD-FS-TSE (ranging from 0.74 to 0.99 mm); however, both morphological sequences showed comparable good results with RMSA of 7.1 to 8.5%. Regarding quantitative T2 mapping, measurements showed T2 relaxation times of about 54 ms with an excellent reproducibility (RMSA) ranging from 3.2 to 4.7%. CONCLUSION: In our study the assessment of cartilage thickness and T2 relaxation times could be performed with high reproducibility in a clinically realizable scan time, demonstrating new possibilities for further investigations into patient groups.

T2 mapping in the knee after microfracture at 3.0 T: correlation of global T2 values and clinical outcome - preliminary results

OBJECTIVE: The aim of our study was to correlate global T2 values of microfracture repair tissue (RT) with clinical outcome in the knee joint. METHODS: We assessed 24 patients treated with microfracture in the knee joint. Magnetic resonance (MR) examinations were performed on a 3T MR unit, T2 relaxation times were obtained with a multi-echo spin-echo technique. T2 maps were obtained using a pixel wise, mono-exponential non-negative least squares fit analysis. Slices covering the cartilage RT were selected and region of interest analysis was done. An individual T2 index was calculated with global mean T2 of the RT and global mean T2 of normal, hyaline cartilage. The Lysholm score and the International Knee Documentation Committee (IKDC) knee evaluation forms were used for the assessment of clinical outcome. Bivariate correlation analysis and a paired, two tailed t test were used for statistics. RESULTS: Global T2 values of the RT [mean 49.8ms, standards deviation (SD) 7.5] differed significantly (P<0.001) from global T2 values of normal, hyaline cartilage (mean 58.5ms, SD 7.0). The T2 index ranged from 61.3 to 101.5. We found the T2 index to correlate with outcome of the Lysholm score (r(s)=0.641, P<0.001) and the IKDC subjective knee evaluation form (r(s)=0.549, P=0.005), whereas there was no correlation with the IKDC knee form (r(s)=-0.284, P=0.179). CONCLUSION: These findings indicate that T2 mapping is sensitive to assess RT function and provides additional information to morphologic MRI in the monitoring of microfracture.

An ethnobotanical survey of the economic and cultural significance of non-timber forest products in the southwest Rhodope Mountain region of Bulgaria

The people of the southwestern Rhodope Mountains of Bulgaria live in small, mountainous villages and rural areas. They rely on berries, herbs, and mushrooms provided by the forest and maintain a lifestyle and culture of gathering them. This study determined the economic and landscape concentration of Non-Timber Forest Products (NTFPs) and how this has changed in the past twenty years in the region of Garmen. The objective was to gauge the cultural and economic significance of NTFPs in the lives of the people who live there. Data was collected using informal, open-ended interviews and through participant observation. Results indicate that ethnicity influence how resources are utilized. Roma people collect mushrooms for income generation; Orthodox Bulgarians gather herbs, berries, and mushrooms for medicinal purposes, to supplement their diets, and to carry on traditions. Bulgarian Muslims collect for a combination of the aforementioned reasons. Changes that occur in the forests affect each of the ethnic groups in different ways and forest management practices should include people’s knowledge and uses of NTFPs.