Intensity warping for multisite MRI harmonization.

In multisite neuroimaging studies there is often unwanted technical variation across scanners and sites. These "scanner effects" can hinder detection of biological features of interest, produce inconsistent results, and lead to spurious associations. We propose mica (multisite image harmonization by cumulative distribution function alignment), a tool to harmonize images taken on different scanners by identifying and removing within-subject scanner effects. Our goals in the present study were to (1) establish a method that removes scanner effects by leveraging multiple scans collected on the same subject, and, building on this, (2) develop a technique to quantify scanner effects in large multisite studies so these can be reduced as a preprocessing step. We illustrate scanner effects in a brain MRI study in which the same subject was measured twice on seven scanners, and assess our method's performance in a second study in which ten subjects were scanned on two machines. We found that unharmonized images were highly variable across site and scanner type, and our method effectively removed this variability by aligning intensity distributions. We further studied the ability to predict image harmonization results for a scan taken on an existing subject at a new site using cross-validation.

Distinguishing Extravascular from Intravascular Ferumoxytol Pools within the Brain: Proof of Concept in Patients with Treated Glioblastoma.

BACKGROUND AND PURPOSE: Glioblastoma-associated macrophages are a major constituent of the immune response to therapy and are known to engulf the iron-based MR imaging contrast agent, ferumoxytol. Current ferumoxytol MR imaging techniques for localizing macrophages are confounded by contaminating intravascular signal. The aim of this study was to assess the utility of a newly developed MR imaging technique, segregation and extravascular localization of ferumoxytol imaging, for differentiating extravascular-from-intravascular ferumoxytol contrast signal at a delayed 24-hour imaging time point. MATERIALS AND METHODS: Twenty-three patients with suspected post-chemoradiotherapy glioblastoma progression underwent ferumoxytol-enhanced SWI. Segregation and extravascular localization of ferumoxytol imaging maps were generated as the voxelwise difference of the delayed (24 hours) from the early (immediately after administration) time point SWI maps. Continuous segregation and extravascular localization of ferumoxytol imaging map values were separated into positive and negative components. Image-guided biologic correlation was performed. RESULTS: Negative segregation and extravascular localization of ferumoxytol imaging values correlated with early and delayed time point SWI values, demonstrating that intravascular signal detected in the early time point persists into the delayed time point. Positive segregation and extravascular localization of ferumoxytol imaging values correlated only with delayed time point SWI values, suggesting successful detection of the newly developed extravascular signal. CONCLUSIONS: Segregation and extravascular localization of ferumoxytol MR imaging improves on current techniques by eliminating intrinsic tissue and intravascular ferumoxytol signal and may inform glioblastoma outcomes by serving as a more specific metric of macrophage content compared with uncorrected T1 and SWI techniques.

An Automated Statistical Technique for Counting Distinct Multiple Sclerosis Lesions.

BACKGROUND AND PURPOSE: Lesion load is a common biomarker in multiple sclerosis, yet it has historically shown modest association with clinical outcome. Lesion count, which encapsulates the natural history of lesion formation and is thought to provide complementary information, is difficult to assess in patients with confluent (ie, spatially overlapping) lesions. We introduce a statistical technique for cross-sectionally counting pathologically distinct lesions. MATERIALS AND METHODS: MR imaging was used to assess the probability of a lesion at each location. The texture of this map was quantified using a novel technique, and clusters resembling the center of a lesion were counted. Validity compared with a criterion standard count was demonstrated in 60 subjects observed longitudinally, and reliability was determined using 14 scans of a clinically stable subject acquired at 7 sites. RESULTS: The proposed count and the criterion standard count were highly correlated (r = 0.97, P < .001) and not significantly different (t59 = -.83, P = .41), and the variability of the proposed count across repeat scans was equivalent to that of lesion load. After accounting for lesion load and age, lesion count was negatively associated (t58 = -2.73, P < .01) with the Expanded Disability Status Scale. Average lesion size had a higher association with the Expanded Disability Status Scale (r = 0.35, P < .01) than lesion load (r = 0.10, P = .44) or lesion count (r = -.12, P = .36) alone. CONCLUSIONS: This study introduces a novel technique for counting pathologically distinct lesions using cross-sectional data and demonstrates its ability to recover obscured longitudinal information. The proposed count allows more accurate estimation of lesion size, which correlated more closely with disability scores than either lesion load or lesion count alone.

Volumetric Analysis from a Harmonized Multisite Brain MRI Study of a Single Subject with Multiple Sclerosis.

BACKGROUND AND PURPOSE: MR imaging can be used to measure structural changes in the brains of individuals with multiple sclerosis and is essential for diagnosis, longitudinal monitoring, and therapy evaluation. The North American Imaging in Multiple Sclerosis Cooperative steering committee developed a uniform high-resolution 3T MR imaging protocol relevant to the quantification of cerebral lesions and atrophy and implemented it at 7 sites across the United States. To assess intersite variability in scan data, we imaged a volunteer with relapsing-remitting MS with a scan-rescan at each site. MATERIALS AND METHODS: All imaging was acquired on Siemens scanners (4 Skyra, 2 Tim Trio, and 1 Verio). Expert segmentations were manually obtained for T1-hypointense and T2 (FLAIR) hyperintense lesions. Several automated lesion-detection and whole-brain, cortical, and deep gray matter volumetric pipelines were applied. Statistical analyses were conducted to assess variability across sites, as well as systematic biases in the volumetric measurements that were site-related. RESULTS: Systematic biases due to site differences in expert-traced lesion measurements were significant (P < .01 for both T1 and T2 lesion volumes), with site explaining >90% of the variation (range, 13.0-16.4 mL in T1 and 15.9-20.1 mL in T2) in lesion volumes. Site also explained >80% of the variation in most automated volumetric measurements. Output measures clustered according to scanner models, with similar results from the Skyra versus the other 2 units. CONCLUSIONS: Even in multicenter studies with consistent scanner field strength and manufacturer after protocol harmonization, systematic differences can lead to severe biases in volumetric analyses.

Magnetic resonance imaging of the proximal upper extremity musculature in boys with Duchenne muscular dystrophy.

There is a pressing need for biomarkers and outcomes that can be used across disease stages in Duchenne muscular dystrophy (DMD), to facilitate the inclusion of a wider range of participants in clinical trials and to improve our understanding of the natural history of DMD. Quantitative magnetic resonance imaging (qMRI) and spectroscopy (MRS) biomarkers show considerable promise in both the legs and forearms of individuals with DMD, but have not yet been examined in functionally important proximal upper extremity muscles such as the biceps brachii and deltoid. The primary objective of this study was to examine the feasibility of implementing qMRI and MRS biomarkers in the proximal upper extremity musculature, and the secondary objective was to examine the relationship between MR measures of arm muscle pathology and upper extremity functional endpoints. Biomarkers included MRS and MRI measures of fat fraction and transverse relaxation time (T 2). The MR exam was well tolerated in both ambulatory and non-ambulatory boys. qMR biomarkers differentiated affected and unaffected participants and correlated strongly with upper extremity function (r = 0.91 for biceps brachii T 2 versus performance of upper limb score). These qMR outcome measures could be highly beneficial to the neuromuscular disease community, allowing measurement of the quality of functionally important muscles across disease stages to understand the natural history of DMD and particularly to broaden the opportunity for clinical trial participation.

Characterizing the white matter hyperintensity penumbra with cerebral blood flow measures.

OBJECTIVE: White matter hyperintensities (WMHs) are common with age, grow over time, and are associated with cognitive and motor impairments. Mechanisms underlying WMH growth are unclear. We aimed to determine the presence and extent of decreased normal appearing white matter (NAWM) cerebral blood flow (CBF) surrounding WMHs to identify 'WM at risk', or the WMH CBF penumbra. We aimed to further validate cross-sectional finding by determining whether the baseline WMH penumbra CBF predicts the development of new WMHs at follow-up. METHODS: Sixty-one cognitively intact elderly subjects received 3 T MPRAGE, FLAIR, and pulsed arterial spin labeling (PASL). Twenty-four subjects returned for follow-up MRI. The inter-scan interval was 18 months. A NAWM layer mask, comprised of fifteen layers, 1 mm thick each surrounding WMHs, was generated for periventricular (PVWMH) and deep (DWMH) WMHs. Mean CBF for each layer was computed. New WMH and persistent NAWM voxels for each penumbra layer were defined from follow-up MRI. RESULTS: CBF in the area surrounding WMHs was significantly lower than the total brain NAWM, extending approximately 12 mm from both the established PVWMH and DWMH. Voxels with new WMH at follow-up had significantly lower baseline CBF than voxels that maintained NAWM, suggesting that baseline CBF can predict the development of new WMHs over time. CONCLUSIONS: A CBF penumbra exists surrounding WMHs, which is associated with future WMH expansion. ASL MRI can be used to monitor interventions to increase white matter blood flow for the prevention of further WM damage and its cognitive and motor consequences.

Loss of pons-to-hypothalamic white matter tracks in brainstem obesity.

Hyperphagia and obesity have been reported following damage to the hypothalamus in humans. Other brain sites are also postulated to be involved in the control of food intake and body weight regulation, such as the amygdala and brainstem. The brainstem, however, is thought to primarily integrate short-term meal-related signals but not affect long-term alterations in body weight, which is controlled by higher centers. The objective of this study was to identify structural pathways damaged in a patient with a brainstem cavernoma who experienced sudden onset of hyperphagia and >50 kg weight gain in <1 year following surgical drainage via a midline suboccipital craniotomy. Diffusion tensor imaging revealed loss of nerve fiber connections between her brainstem, hypothalamus and higher brain centers with preservation of motor tracks. Imaging and endocrine testing confirmed normal hypothalamic structure and function. Gastric bypass surgery restored normal appetite and body weight to baseline. This is the first report of 'brainstem obesity' and adds to the brain regions that can determine the long-term body weight set point in humans.

Comparison of 3 T MRI and CT for the measurement of visceral and subcutaneous adipose tissue in humans.

OBJECTIVE: CT is considered the gold standard imaging modality for measurement of visceral adipose tissue area. However, as CT imaging exposes subjects to ionising radiation, a comparable imaging technique without this exposure is desirable, such as MRI. Therefore, we compared the agreement of measures of visceral adipose tissue and subcutaneous adipose tissue area from single-slice images obtained at the umbilicus using a 3 T MRI scanner with single-slice images obtained via CT scan. METHODS: 64 images were obtained from 27 subjects who underwent MRI and CT scanning on the same day, after 10-12 hours of fasting. Visceral and subcutaneous adipose tissue depots were manually separated and quantified using a multimodality image-processing software program. RESULTS: We found good agreement between CT and MRI for the measurement of both visceral adipose tissue and subcutaneous adipose tissue. Bland-Altman difference analysis demonstrated a mean bias of -2.9% (as a portion of total abdominal area) for visceral adipose tissue and +0.4% for subcutaneous adipose tissue, as measured by MRI compared with CT. CONCLUSION: MRI is a safe, accurate and precise imaging modality for measuring both visceral and subcutaneous adipose tissue, making it a favourable alternative to CT for quantification of these adipose depots.

Neuroimaging features of neurodegeneration with brain iron accumulation.

NBIA characterizes a class of neurodegenerative diseases that feature a prominent extrapyramidal movement disorder, intellectual deterioration, and a characteristic deposition of iron in the basal ganglia. The diagnosis of NBIA is made on the basis of the combination of representative clinical features along with MR imaging evidence of iron accumulation. In many cases, confirmatory molecular genetic testing is now available as well. A number of new subtypes of NBIA have recently been described, with distinct neuroradiologic and clinical features. This article outlines the known subtypes of NBIA, delineates their clinical and radiographic features, and suggests an algorithm for evaluation.

Brain functional magnetic resonance imaging response to glucose and fructose infusions in humans.

AIMS: In animals, intracerebroventricular glucose and fructose have opposing effects on appetite and weight regulation. In humans, functional brain magnetic resonance imaging (fMRI) studies during glucose ingestion or infusion have demonstrated suppression of hypothalamic signalling, but no studies have compared the effects of glucose and fructose. We therefore sought to determine if the brain response differed to glucose vs. fructose in humans independently of the ingestive process. METHODS: Nine healthy, normal weight subjects underwent blood oxygenation level dependent (BOLD) fMRI measurements during either intravenous (IV) glucose (0.3 mg/kg), fructose (0.3 mg/kg) or saline, administered over 2 min in a randomized, double-blind, crossover study. Blood was sampled every 5 min during a baseline period and following infusion for 60 min in total for glucose, fructose, lactate and insulin levels. RESULTS: No significant brain BOLD signal changes were detected in response to IV saline. BOLD signal in the cortical control areas increased during glucose infusion (p = 0.002), corresponding with increased plasma glucose and insulin levels. In contrast, BOLD signal decreased in the cortical control areas during fructose infusion (p = 0.006), corresponding with increases of plasma fructose and lactate. Neither glucose nor fructose infusions significantly altered BOLD signal in the hypothalamus. CONCLUSION: In normal weight humans, cortical responses as assessed by BOLD fMRI to infused glucose are opposite to those of fructose. Differential brain responses to these sugars and their metabolites may provide insight into the neurologic basis for dysregulation of food intake during high dietary fructose intake.

Inheritance of inflorescence architecture in sorghum.

The grass inflorescence is the primary food source for humanity, and has been repeatedly shaped by human selection during the domestication of different cereal crops. Of all major cultivated cereals, sorghum [Sorghum bicolor (L.) Moench] shows the most striking variation in inflorescence architecture traits such as branch number and branch length, but the genetic basis of this variation is little understood. To study the inheritance of inflorescence architecture in sorghum, 119 recombinant inbred lines from an elite by exotic cross were grown in three environments and measured for 15 traits, including primary, secondary, and tertiary inflorescence branching. Eight characterized genes that are known to control inflorescence architecture in maize (Zea mays L.) and other grasses were mapped in sorghum. Two of these candidate genes, Dw3 and the sorghum ortholog of ramosa2, co-localized precisely with QTL of large effect for relevant traits. These results demonstrate the feasibility of using genomic and mutant resources from maize and rice (Oryza sativa L.) to investigate the inheritance of complex traits in related cereals.

Fertility restorer locus Rf1 [corrected] of sorghum (Sorghum bicolor L.) encodes a pentatricopeptide repeat protein not present in the colinear region of rice chromosome 12.

With an aim to clone the sorghum fertility restorer gene Rf1, a high-resolution genetic and physical map of the locus was constructed. The Rf1 locus was resolved to a 32-kb region spanning four open reading frames: a plasma membrane Ca(2+)-ATPase, a cyclin D-1, an unknown protein, and a pentatricopeptide repeat (PPR13) gene family member. An approximately 19-kb region spanning the cyclin D-1 and unknown protein genes was completely conserved between sterile and fertile plants as was the sequence spanning the coding region of the Ca(2+)-ATPase. In contrast, 19 sequence polymorphisms were located in an approximately 7-kb region spanning PPR13, and all markers cosegregated with the fertility restoration phenotype. PPR13 was predicted to encode a mitochondrial-targeted protein containing a single exon with 14 PPR repeats, and the protein is classified as an E-type PPR subfamily member. To permit sequence-based comparison of the sorghum and rice genomes in the Rf1 region, 0.53 Mb of sorghum chromosome 8 was sequenced and compared to the colinear region of rice chromosome 12. Genome comparison revealed a mosaic pattern of colinearity with an approximately 275-kb gene-poor region with little gene conservation and an adjacent, approximately 245-kb gene-rice region that is more highly conserved between rice and sorghum. Despite being located in a region of high gene conservation, sorghum PPR13 was not located in a colinear position on rice chromosome 12. The present results suggest that sorghum PPR13 represents a potential candidate for the sorghum Rf1 gene, and its presence in the sorghum genome indicates a single gene transposition event subsequent to the divergence of rice and sorghum ancestors.

A molecular cytogenetic map of sorghum chromosome 1. Fluorescence in situ hybridization analysis with mapped bacterial artificial chromosomes.

We used structural genomic resources for Sorghum bicolor (L.) Moench to target and develop multiple molecular cytogenetic probes that would provide extensive coverage for a specific chromosome of sorghum. Bacterial artificial chromosome (BAC) clones containing molecular markers mapped across sorghum linkage group A were labeled as probes for fluorescence in situ hybridization (FISH). Signals from single-, dual-, and multiprobe BAC-FISH to spreads of mitotic chromosomes and pachytene bivalents were associated with the largest sorghum chromosome, which bears the nucleolus organizing region (NOR). The order of individual BAC-FISH loci along the chromosome was fully concordant to that of marker loci along the linkage map. In addition, the order of several tightly linked molecular markers was clarified by FISH analysis. The FISH results indicate that markers from the linkage map positions 0.0-81.8 cM reside in the short arm of chromosome 1 whereas markers from 81.8-242.9 cM are located in the long arm of chromosome 1. The centromere and NOR were located in a large heterochromatic region that spans approximately 60% of chromosome 1. In contrast, this region represents only 0.7% of the total genetic map distance of this chromosome. Variation in recombination frequency among euchromatic chromosomal regions also was apparent. The integrated data underscore the value of cytological data, because minor errors and uncertainties in linkage maps can involve huge physical regions. The successful development of multiprobe FISH cocktails suggests that it is feasible to develop chromosome-specific "paints" from genomic resources rather than flow sorting or microdissection and that when applied to pachytene chromatin, such cocktails provide an especially powerful framework for mapping. Such a molecular cytogenetic infrastructure would be inherently cross-linked with other genomic tools and thereby establish a cytogenomics system with extensive utility in development and application of genomic resources, cloning, transgene localization, development of plant "chromonomics," germplasm introgression, and marker-assisted breeding. In combination with previously reported work, the results indicate that a sorghum cytogenomics system would be partially applicable to other gramineous genera.

Reanalysis of multislice (1)H MRSI in amyotrophic lateral sclerosis.

The goal of this work was to reexamine previously published (1) brain spectroscopy data of abnormal metabolite ratios in amyotrophic lateral sclerosis (ALS). Toward this goal, (1)H MR spectroscopic imaging data from 10 ALS and nine control subjects were reanalyzed using improved data analysis techniques, including automated curve fitting and tissue-volume correction. In the motor cortex of ALS, N-acetyl aspartate (NAA) was 23% (P = 0.004) lower than in controls, and in the posterior internal capsule of ALS choline compounds (Cho) were 20% (P = 0.02) higher. This demonstrates that the metabolite ratio changes in ALS were due to NAA loss in the motor cortex (as expected) and Cho increase in the posterior internal capsule (not expected). Magn Reson Med 45:513-516, 2001.

Resting brain metabolic activity in a 4 tesla magnetic field.

MRI is a major tool for mapping brain function; thus it is important to assess potential effects on brain neuronal activity attributable to the requisite static magnetic field. This study used positron emission tomography (PET) and (18)F-deoxyglucose ((18)FDG) to measure brain glucose metabolism (a measure of brain function) in 12 subjects while their heads were in a 4 T MRI field during the (18)FDG uptake period. The results were compared with those obtained when the subjects were in the earth's field (PET scanner), and when they were in a simulated MRI environment in the PET instrument that imitated the restricted visual field of the MRI experiment. Whole-brain metabolism, as well as metabolism in occipital cortex and posterior cingulate gyrus, was lower in the real (4 T) and simulated (0 T) MRI environments compared with the PET. This suggests that the metabolic differences are due mainly to the visual field differences characteristic of the MRI and PET instruments. We conclude that a static magnetic field of 4 T does not in itself affect this fairly sensitive measure of brain activity.

Determination of the MRI contrast agent concentration time course in vivo following bolus injection: effect of equilibrium transcytolemmal water exchange.

For bolus-tracking studies, it is commonly assumed that CR concentration bears a linear relationship with the measured (usually longitudinal) (1)H(2)O relaxation rate constant, R*(1) identical with(T(1) *)(-1). This requires that equilibrium transcytolemmal water exchange be in the fast exchange limit (FXL). However, though systems remain in fast exchange, the FXL will not usually obtain. Here, the consequences are considered: 1) the measurement of R(1) * itself can be affected, 2) the resultant non-linear [CR]-dependence causes significant error by assuming FXL, 3) the thermodynamic [CR] (based on the space in which CR is actually distributed) can be determined, 4) transcytolemmal water permeability may be estimated, and 5) the pharmacokinetic parameters can be factored. For a 30-sec, 0.17 mmol/kg dose of GdDTPA(2-), the FXL assumption underestimates the [CR] maximum in rat thigh muscle by a factor of almost two. Similar results are obtained for a rat brain GS-9L gliosarcoma tumor model.

4.0 T water proton T1 relaxation times in normal human brain and during acute ethanol intoxication.

BACKGROUND: It has been reported that acute ethanol intoxication decreases the brain water proton magnetic resonance T1 values, an effect that has been interpreted to indicate brain dehydration during this condition. Because water macromolecular interactions largely determine tissue water T1, another possible explanation for reduced brain water proton T1 values is that the interaction between water and brain macromolecules is altered by ethanol. METHODS: A 4.0 T magnetic resonance imaging (MRI) instrument was used to measure brain water proton T1 relaxation times before, during, and after ethanol intoxication (dose, 0.75 mg/kg) in healthy controls. RESULTS: The T1 relaxation times as assessed with MRI were highly reproducible. The mean, paired ethanol-induced differences in T1 were -0.004 +/- 0.007 sec (mean +/- standard deviation) for white matter and 0.010 +/- 0.015 sec for internal gray matter structures, neither of which was significant. CONCLUSIONS: This reasonably sensitive measurement does not support the view that tissue water content or water macromolecule interactions are significantly altered in the brain during acute alcohol intoxication in otherwise healthy subjects.

1H MRSI comparison of white matter and lesions in primary progressive and relapsing-remitting MS.

OBJECTIVE: To compare brain metabolite levels in patients with primary progressive (PP) and relapsing remitting (RR) MS and controls. HYPOTHESES: (1) creatine (Cr), a putative marker of gliosis, is elevated and N-acetylaspartate (NAA), a putative marker of axonal density and functional integrity, is reduced in PPMS lesions and normal appearing white matter (NAWM) compared to control white matter; (2) The pattern of metabolite change in PPMS is different than in RRMS. METHODS: MRI and proton magnetic resonance spectroscopic imaging (1H MRSI) were collected from 15 PPMS patients, 13 RRMS patients, and 20 controls. RESULTS: Cr was increased in PPMS NAWM compared to controls (P=0.035), and compared to RRMS NAWM (P=0.038). Cr was increased in focal MRI lesions from PPMS compared to lesions from RRMS (P=0.044) and compared to control white matter (P=0.041). NAA was similarly reduced in PPMS and RRMS NAWM compared to control. NAA was similarly reduced in PPMS and RRMS lesions, compared to control white matter. CONCLUSIONS: Creatine is higher in PPMS than RRMS NAWM and focal lesions. This observation is consistent with the notion that progressive disability in PPMS reflects increased gliosis and axonal loss whereas disability in RRMS reflects the cumulative effects of acute inflammatory lesions and axonal loss.

Activity of antifungal proteins against mold in sorghum caryopses in the field.

Sorghums were stressed with pathogenic fungi and sprinkling to determine relationships between changes in chitinase and sormatin in caryopses and grain mold resistance. Panicles of 10 cultivars differing in mold resistance and accumulation of antifungal proteins (AFPs) were inoculated at anthesis with Fusarium moniliforme and Curvularia lunata spores. Panicles were sampled at 30 and 50 days after anthesis, and caryopses were evaluated for chitinase and sormatin using western blots. Sprinkling panicles (to mimic rainfall) decreased sormatin and chitinase in most cultivars. Inoculation decreased AFPs in susceptible cultivars, but resistant cultivars maintained or increased AFPs in caryopses. Grain mold resistance corresponded to induction of AFP synthesis in response to sprinkling, fungal stress, and/or adverse field conditions. Sormatin and chitinase appear to be an active part of the defense mechanism of the caryopsis against grain mold.

Antifungal proteins and grain mold resistance in sorghum with nonpigmented testa.

Levels of four antifungal proteins (AFPs) were determined in mature caryopses (40-45 days after anthesis) of eight grain mold resistant (GMR) and eight susceptible (GMS) sorghum lines using the immunoblot technique. These 16 lines came from the same cross and were selected for high and low grain mold resistance. The 16 lines were grown in eight environments over three years. In the environments with grain mold incidence, levels of sormatin, chitinases, and ribosomal inactivating proteins (RIP) in the GMR group were higher than those in the GMS group. In a grain mold-free environment, the GMR group had higher RIP and lower beta-1,3-glucanase levels than the GMS group. Unlike the GMS group, chitinase, sormatin, and RIP levels in the GMR group were higher in the environments with grain mold than in the mold-free environment. AFPs correlated among themselves and with grain mold resistance. Grain mold infection pressure caused GMR lines to induce and/or retain more AFPs compared to GMS lines. The coexpression of these four AFPs may be a necessary prerequisite for resistance to grain mold in sorghums without a pigmented testa.