Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Noninvasive coronary vessel wall and plaque imaging with magnetic resonance imaging.

BACKGROUND: Conventional x-ray angiography frequently underestimates the true burden of atherosclerosis. Although intravascular ultrasound allows for imaging of coronary plaque, this invasive technique is inappropriate for screening or serial examinations. We therefore sought to develop a noninvasive free-breathing MR technique for coronary vessel wall imaging. We hypothesized that such an approach would allow for in vivo imaging of coronary atherosclerosis. METHODS AND RESULTS: Ten subjects, including 5 healthy adult volunteers (aged 35+/-17 years, range 19 to 56 years) and 5 patients (aged 60+/-4 years, range 56 to 66 years) with x-ray-confirmed coronary artery disease (CAD), were studied with a T2-weighted, dual-inversion, fast spin-echo MR sequence. Multiple adjacent 5-mm cross-sectional images of the proximal right coronary artery were obtained with an in-plane resolution of 0.5x1.0 mm. A right hemidiaphragmatic navigator was used to facilitate free-breathing MR acquisition. Coronary vessel wall images were readily acquired in all subjects. Both coronary vessel wall thickness (1.5+/-0.2 versus 1.0+/-0.2 mm) and wall area (21.2+/-3.1 versus 13.7+/-4.2 mm(2)) were greater in patients with CAD (both P:<0.02 versus healthy adults). CONCLUSIONS: In vivo free-breathing coronary vessel wall and plaque imaging with MR has been successfully implemented in humans. Coronary wall thickness and wall area were significantly greater in patients with angiographic CAD. The presented technique may have potential applications in patients with known or suspected atherosclerotic CAD or for serial evaluation after pharmacological intervention.

Three-dimensional patient-specific dosimetry in radioimmunotherapy with 90Y-ibritumomab-tiuxetan.

Aim of the present article was to perform three-dimensional (3D) single photon emission tomography-based dosimetry in radioimmunotherapy (RIT) with (90)Y-ibritumomab-tiuxetan. A custom MATLAB-based code was used to elaborate 3D images and to compare average 3D doses to lesions and to organs at risk (OARs) with those obtained with planar (2D) dosimetry. Our 3D dosimetry procedure was validated through preliminary phantom studies using a body phantom consisting of a lung insert and six spheres with various sizes. In phantom study, the accuracy of dose determination of our imaging protocol decreased when the object volume decreased below 5 mL, approximately. The poorest results were obtained for the 2.58 mL and 1.30 mL spheres where the dose error evaluated on corrected images with regard to the theoretical dose value was -12.97% and -18.69%, respectively. Our 3D dosimetry protocol was subsequently applied on four patients before RIT with (90)Y-ibritumomab-tiuxetan for a total of 5 lesions and 4 OARs (2 livers, 2 spleens). In patient study, without the implementation of volume recovery technique, tumor absorbed doses calculated with the voxel-based approach were systematically lower than those calculated with the planar protocol, with average underestimation of -39% (range from -13.1% to -62.7%). After volume recovery, dose differences reduce significantly, with average deviation of -14.2% (range from -38.7.4% to +3.4%, 1 overestimation, 4 underestimations). Organ dosimetry in one case overestimated, in the other underestimated the dose delivered to liver and spleen. However, both for 2D and 3D approach, absorbed doses to organs per unit administered activity are comparable with most recent literature findings.

Tertiary Himalayan structures and metamorphism in the Kulu Valley (Mandi-Khoksar transect of the Western Himalaya) - Shikar-Beh-nappe and crystalline nappe

The Crystalline Nappe of the High Himalayan Crystalline has been examined along the Kulu Valley and its vicinity (Mandi-Khoksar transect). This nappe was believed to have undergone deformation related only to its transport towards the SW essentially during the `'Main Central Thrust event''. New data has led to the conclusion that during the Himalayan orogeny, two distinctive phases, related to two opposite transport directions, characterize the evolution of this part of the chain, before the creation of the late NE-vergent backfolding. The first phase corresponds to an early NE-vergent folding and thrusting, creating the Tandi Syncline and the NE-oriented Shikar Beh Nappe stack, with a displacement amplitude of about 50 km. Two schistosities, together with a strong stretching lineation are developed at a deep tectonic level under amphibolite facies conditions (kyanite-staurolite-garnet-two mica schists). At a higher tectonic level and in the southern part of the section (Tandy Syncline and southern Kulu Valley between Kulu and Mandi) one or two schistosities are developed in the greenschist facies grade rocks (garnet-biotite and biotite schists). These structures and the associated Barrovian type metamorphism are all related to the NE-verging Shikar Beh Nappe. The creation of the NE-verging Shikar Beh Nappe may be explained by the reactivation of a SW dipping listric normal fault of the N Indian flexural passive margin, during the early stages of the Himalayan orogeny. In the second phase, the still hot metamorphic rocks of the Shikar Beh Nappe were folded and thrust towards the SW (mainly along the MBT and the MCT with a displacement in excess of 100 km) onto the cold, low-grade metamorphic rocks of the Larji-Kulu-Rampur Window or, near Mandi, on the non-metamorphic sandstones of the Ganges Molasse (Siwaliks). Sense of shear criteria and a strong NE-SW stretching-lineation indicate that the Crystalline Nappe has been overthrusted towards the SW. Thermometry on synkinematically crystallised garnet-biotite and garnet-hornblende pairs reveals the lower amphibolite facies temperature conditions related to the Crystalline Nappe formation. From the muscovite and biotite Rb-Sr cooling ages, the Shikar Beh Nappe emplacement occurred before 32 Ma and the southwestward thrusting of the Crystalline Nappe began before 21 Ma. Our model involving two opposite directions of thrusting goes against the conventional idea of only one main SW-oriented transport direction in the High Himalayan Crystalline Nappes.

Relationship between motion of coronary arteries and diaphragm during free breathing: lessons from real-time MR imaging.

OBJECTIVE: Diaphragmatic navigators are frequently used in free-breathing coronary MR angiography, either to gate or prospectively correct slice position or both. For such approaches, a constant relationship between coronary and diaphragmatic displacement throughout the respiratory cycle is assumed. The purpose of this study was to evaluate the relationship between diaphragmatic and coronary artery motion during free breathing. SUBJECTS AND METHODS: A real-time echoplanar MR imaging sequence was used in 12 healthy volunteers to obtain 30 successive images each (one per cardiac cycle) that included the left main coronary artery and the domes of both hemidiaphragms. The coronary artery and diaphragm positions (relative to isocenter) were determined and analyzed for effective diaphragmatic gating windows of 3, 5, and 7 mm (diaphragmatic excursions of 0-3, 0-5, and 0-7 mm from the end-expiratory position, respectively). RESULTS: Although the mean slope correlating the displacement of the right diaphragm and the left main coronary artery was approximately 0.6 for all diaphragmatic gating windows, we also found great variability among individual volunteers. Linear regression slopes varied from 0.17 to 0.93, and r2 values varied from .04 to .87. CONCLUSION: Wide individual variability exists in the relationship between coronary and diaphragmatic respiratory motion during free breathing. Accordingly, coronary MR angiographic approaches that use diaphragmatic navigator position for prospective slice correction may benefit from patient-specific correction factors. Alternatively, coronary MR angiography may benefit from a more direct assessment of the respiratory displacement of the heart and coronary arteries, using left ventricular navigators.

Renal arteries: navigator-gated balanced fast field-echo projection MR angiography with aortic spin labeling: initial experience.

A cardiac-triggered free-breathing three-dimensional balanced fast field-echo projection magnetic resonance (MR) angiographic sequence with a two-dimensional pencil-beam aortic labeling pulse was developed for the renal arteries. For data acquisition during free breathing in eight healthy adults and seven consecutive patients with renal artery disease, real-time navigator technology was implemented. This technique allows high-spatial-resolution and high-contrast renal MR angiography and visualization of renal artery stenosis without exogenous contrast agent or breath hold. Initial promising results warrant larger clinical studies.

A fast 3D approach for coronary MRA.

Two-dimensional (2D)-breath-hold coronary magnetic resonance angiography (MRA) has been shown to be a fast and reliable method to depict the proximal coronary arteries. Recent developments, however, allow for free-breathing navigator gated and navigator corrected three-dimensional (3D) coronary MRA. These 3D approaches have potential for improved signal-to-noise ratio (SNR) and allow for the acquisition of adjacent thin slices without the misregistration problems known from 2D approaches. Still, a major impediment of a 3D acquisition is the increased scan time. The purpose of this study was the implementation of a free-breathing navigator gated and corrected ultra-fast 3D coronary MRA technique, which allows for scan times of less than 5 minutes. Twelve healthy adult subjects were examined in the supine position using a navigator gated and corrected ECG triggered ultra-fast 3D interleaved gradient echo planar imaging sequence (TFE-EPI). A 3D slab, consisting of 20 slices with a reconstructed slice thickness of 1.5 mm, was acquired with free-breathing. The diastolic TFE-EPI acquisition block was preceded by a T2prep pre-pulse, a diaphragmatic navigator pulse, and a fat suppression pre-pulse. With a TR of 19 ms and an effective TE of 5.4 ms, the duration of the data acquisition window duration was 38 ms. The in-plane spatial resolution was 1.0-1.3 mm*1.5-1.9 mm. In all cases, the entire left main (LM) and extensive portions of the left anterior descending (LAD) and right coronary artery (RCA) could be visualized with an average scan time for the entire 3D-volume data set of 2:57 +/- 0:51 minutes. Average contiguous vessel length visualized was 53 +/- 11 mm (range: 42 to 75 mm) for the LAD and 84 +/- 14 mm (range: 62 to 112 mm) for the RCA. Contrast-to-noise between coronary blood and myocardium was 5.0 +/- 2.3 for the LM/LAD and 8.0 +/- 2.9 for the RCA, resulting in an excellent suppression of myocardium. We present a new approach for free-breathing 3D coronary MRA, which allows for scan times superior to corresponding 2D coronary MRA approaches, and which takes advantage of the enhanced SNR of 3D acquisitions and the post-processing benefits of thin adjacent slices. The robust image quality and the short average scanning time suggest that this approach may be useful for screening the major coronary arteries or identification of anomalous coronary arteries. J. Magn. Reson. Imaging 1999;10:821-825.

Endocytic and secretory traffic in Arabidopsis merge in the trans-Golgi network/early endosome, an independent and highly dynamic organelle.

Plants constantly adjust their repertoire of plasma membrane proteins that mediates transduction of environmental and developmental signals as well as transport of ions, nutrients, and hormones. The importance of regulated secretory and endocytic trafficking is becoming increasingly clear; however, our knowledge of the compartments and molecular machinery involved is still fragmentary. We used immunogold electron microscopy and confocal laser scanning microscopy to trace the route of cargo molecules, including the BRASSINOSTEROID INSENSITIVE1 receptor and the REQUIRES HIGH BORON1 boron exporter, throughout the plant endomembrane system. Our results provide evidence that both endocytic and secretory cargo pass through the trans-Golgi network/early endosome (TGN/EE) and demonstrate that cargo in late endosomes/multivesicular bodies is destined for vacuolar degradation. Moreover, using spinning disc microscopy, we show that TGN/EEs move independently and are only transiently associated with an individual Golgi stack.

Three-dimensional radiobiological dosimetry of kidneys for treatment planning in peptide receptor radionuclide therapy.

PURPOSE: Peptide receptor radionuclide therapy (PRRT) delivers high absorbed doses to kidneys and may lead to permanent nephropathy. Reliable dosimetry of kidneys is thus critical for safe and effective PRRT. The aim of this work was to assess the feasibility of planning PRRT based on 3D radiobiological dosimetry (3D-RD) in order to optimize both the amount of activity to administer and the fractionation scheme, while limiting the absorbed dose and the biological effective dose (BED) to the renal cortex. METHODS: Planar and SPECT data were available for a patient examined with (111)In-DTPA-octreotide at 0.5 (planar only), 4, 24, and 48 h post-injection. Absorbed dose and BED distributions were calculated for common therapeutic radionuclides, i.e., (111)In, (90)Y and (177)Lu, using the 3D-RD methodology. Dose-volume histograms were computed and mean absorbed doses to kidneys, renal cortices, and medullae were compared with results obtained using the MIRD schema (S-values) with the multiregion kidney dosimetry model. Two different treatment planning approaches based on (1) the fixed absorbed dose to the cortex and (2) the fixed BED to the cortex were then considered to optimize the activity to administer by varying the number of fractions. RESULTS: Mean absorbed doses calculated with 3D-RD were in good agreement with those obtained with S-value-based SPECT dosimetry for (90)Y and (177)Lu. Nevertheless, for (111)In, differences of 14% and 22% were found for the whole kidneys and the cortex, respectively. Moreover, the authors found that planar-based dosimetry systematically underestimates the absorbed dose in comparison with SPECT-based methods, up to 32%. Regarding the 3D-RD-based treatment planning using a fixed BED constraint to the renal cortex, the optimal number of fractions was found to be 3 or 4, depending on the radionuclide administered and the value of the fixed BED. Cumulative activities obtained using the proposed simulated treatment planning are compatible with real activities administered to patients in PRRT. CONCLUSIONS: The 3D-RD treatment planning approach based on the fixed BED was found to be the method of choice for clinical implementation in PRRT by providing realistic activity to administer and number of cycles. While dividing the activity in several cycles is important to reduce renal toxicity, the clinical outcome of fractionated PRRT should be investigated in the future.

Geology of the NW Indian Himalaya

This review paper deals with the geology of the NW Indian Himalaya situated in the states of Jammu and Kashmir, Himachal Pradesh and Garhwal. The models and mechanisms discussed, concerning the tectonic and metamorphic history of the Himalayan range, are based on a new compilation of a geological map and cross sections, as well as on paleomagnetic, stratigraphic, petrologic, structural, metamorphic, thermobarometric and radiometric data. The protolith of the Himalayan range, the North Indian flexural passive margin of the Neo-Tethys ocean, consists of a Lower Proterozoic basement, intruded by 1.8-1.9 Ga bimodal magmatites, overlain by a horizontally stratified sequence of Upper Proterozoic to Paleocene sediments, intruded by 470-500 Ma old Ordovician mainly peraluminous s-type granites, Carboniferous tholeiitic to alkaline basalts and intruded and overlain by Permian tholeiitic continental flood basalts. No elements of the Archaen crystalline basement of the South Indian shield have been identified in the Himalayan range. Deformation of the Himalayan accretionary wedge resulted from the continental collision of India and Asia beginning some 65-55 Ma ago, after the NE-directed underthrusting of the Neo-Tethys oceanic crust below Asia and the formation of the Andean-type 103-50 (-41) Ma old Ladakh batholith to the north of the Indus Suture. Cylindrical in geometry, the Himalayan range consists, from NE to SW, from older to younger tectonic elements, of the following zones: 1) The 25 km wide Ladakh batholith and the Asian mantle wedge form the backstop of the growing Himalayan accretionary wedge. 2) The Indus Suture zone is composed of obducted slices of the oceanic crust, island arcs, like the Dras arc, overlain by Late Cretaceous fore arc basin sediments and the mainly Paleocene to Early Eocene and Miocene epi-sutural intra-continental Indus molasse. 3) The Late Paleocene to Eocene North Himalayan nappe stack, up to 40 km thick prior to erosion, consists of Upper Proterozoic to Paleocene rocks, with the eclogitic and coesite bearing Tso Morari gneiss nappe at its base. It includes a branch of the Central Himalayan detachment, the 22-18 Ma old Zanskar Shear zone that is intruded and dated by the 22 Ma Gumburanjun leucogranite; it reactivates the frontal thrusts of the SW-verging North Himalayan nappes. 4) The late Eocene-Miocene SW-directed High Himalayan or ``Crystalline'' nappe comprises Upper Proterozoic to Mesozoic sediments and Ordovician granites, identical to those of the North Himalayan nappes. The Main Central thrust at its base was created in a zone of Eocene to Early Oligocene anatexis by ductile detachment of the subducted Indian crust, below the pre-existing 25-35 km thick NE-directed Shikar Beh and SW-directed North Himalayan nappe stacks. 5) The late Miocene Lesser Himalayan thrust with the Main Boundary Thrust at its base consists of early Proterozoic to Cambrian rocks intruded by 1.8-1.9 Ga bimodal magmatites. The Subhimalaya is a thrust wedge of Himalayan fore deep basin sediments, composed of the Early Eocene marine Subathu marls and sandstones as well as the up to 8'000 m-thick Miocene to recent Ganga molasse, a coarsening upwards sequence of shales, sandstones and conglomerates. The active frontal thrust is covered by the sediments of the Indus-Ganga plains.

Toward high-resolution myocardial tagging.

Magnetic resonance imaging with preceding tissue tagging is a robust method for assessing cardiac motion of the entire heartbeat cycle with a high degree of accuracy. One limitation of this technique, however, is the low resolution of the obtained displacement map of the labeled points within the myocardium. By a new tagging technique, which is based on the combination of two or more measurements of the same slice but with different grid positions, a highly improved resolution of cardiac motion data can be achieved. In combination with a multi-heart-phase echo-planar imaging sequence, such images with doubled grid frequency can be acquired in two short breath-hold periods.

Does the physical disector method provide an accurate estimation of sensory neuron number in rat dorsal root ganglia?

The physical disector is a method of choice for estimating unbiased neuron numbers; nevertheless, calibration is needed to evaluate each counting method. The validity of this method can be assessed by comparing the estimated cell number with the true number determined by a direct counting method in serial sections. We reconstructed a 1/5 of rat lumbar dorsal root ganglia taken from two experimental conditions. From each ganglion, images of 200 adjacent semi-thin sections were used to reconstruct a volumetric dataset (stack of voxels). On these stacks the number of sensory neurons was estimated and counted respectively by physical disector and direct counting methods. Also, using the coordinates of nuclei from the direct counting, we simulate, by a Matlab program, disector pairs separated by increasing distances in a ganglion model. The comparison between the results of these approaches clearly demonstrates that the physical disector method provides a valid and reliable estimate of the number of sensory neurons only when the distance between the consecutive disector pairs is 60 microm or smaller. In these conditions the size of error between the results of physical disector and direct counting does not exceed 6%. In contrast when the distance between two pairs is larger than 60 microm (70-200 microm) the size of error increases rapidly to 27%. We conclude that the physical dissector method provides a reliable estimate of the number of rat sensory neurons only when the separating distance between the consecutive dissector pairs is no larger than 60 microm.

Assembly of the Escherichia coli RuvABC resolvasome directs the orientation of holliday junction resolution.

Genetic recombination can lead to the formation of intermediates in which DNA molecules are linked by Holliday junctions. Movement of a junction along DNA, by a process known as branch migration, leads to heteroduplex formation, whereas resolution of a junction completes the recombination process. Holliday junctions can be resolved in either of two ways, yielding products in which there has, or has not, been an exchange of flanking markers. The ratio of these products is thought to be determined by the frequency with which the two isomeric forms (conformers) of the Holliday junction are cleaved. Recent studies with enzymes that process Holliday junctions in Escherichia coli, the RuvABC proteins, however, indicate that protein binding causes the junction to adopt an open square-planar configuration. Within such a structure, DNA isomerization can have little role in determining the orientation of resolution. To determine the role that junction-specific protein assembly has in determining resolution bias, a defined in vitro system was developed in which we were able to direct the assembly of the RuvABC resolvasome. We found that the bias toward resolution in one orientation or the other was determined simply by the way in which the Ruv proteins were positioned on the junction. Additionally, we provide evidence that supports current models on RuvABC action in which Holliday junction resolution occurs as the resolvasome promotes branch migration.

Free-breathing 3D coronary MRA: the impact of "isotropic" image resolution.

During conventional x-ray coronary angiography, multiple projections of the coronary arteries are acquired to define coronary anatomy precisely. Due to time constraints, coronary magnetic resonance angiography (MRA) usually provides only one or two views of the major coronary vessels. A coronary MRA approach that allowed for reconstruction of arbitrary isotropic orientations might therefore be desirable. The purpose of the study was to develop a three-dimensional (3D) coronary MRA technique with isotropic image resolution in a relatively short scanning time that allows for reconstruction of arbitrary views of the coronary arteries without constraints given by anisotropic voxel size. Eight healthy adult subjects were examined using a real-time navigator-gated and corrected free-breathing interleaved echoplanar (TFE-EPI) 3D-MRA sequence. Two 3D datasets were acquired for the left and right coronary systems in each subject, one with anisotropic (1.0 x 1.5 x 3.0 mm, 10 slices) and one with "near" isotropic (1.0 x 1.5 x 1.0 mm, 30 slices) image resolution. All other imaging parameters were maintained. In all cases, the entire left main (LM) and extensive portions of the left anterior descending (LAD) and the right coronary artery (RCA) were visualized. Objective assessment of coronary vessel sharpness was similar (41% +/- 5% vs. 42% +/- 5%; P = NS) between in-plane and through-plane views with "isotropic" voxel size but differed (32% +/- 7% vs. 23% +/- 4%; P < 0.001) with nonisotropic voxel size. In reconstructed views oriented in the through-plane direction, the vessel border was 86% more defined (P < 0.01) for isotropic compared with anisotropic images. A smaller (30%; P < 0.001) improvement was seen for in-plane reconstructions. Vessel diameter measurements were view independent (2.81 +/- 0.45 mm vs. 2.66 +/- 0.52 mm; P = NS) for isotropic, but differed (2.71 +/- 0.51 mm vs. 3.30 +/- 0.38 mm; P < 0.001) between anisotropic views. Average scanning time was 2:31 +/- 0:57 minutes for anisotropic and 7:11 +/- 3:02 minutes for isotropic image resolution (P < 0.001). We present a new approach for "near" isotropic 3D coronary artery imaging, which allows for reconstruction of arbitrary views of the coronary arteries. The good delineation of the coronary arteries in all views suggests that isotropic 3D coronary MRA might be a preferred technique for the assessment of coronary disease, although at the expense of prolonged scan times. Comparative studies with conventional x-ray angiography are needed to investigate the clinical utility of the isotropic strategy.

Direct three-dimensional myocardial strain tensor quantification and tracking using zHARP.

Images of myocardial strain can be used to diagnose heart disease, plan and monitor treatment, and to learn about cardiac structure and function. Three-dimensional (3D) strain is typically quantified using many magnetic resonance (MR) images obtained in two or three orthogonal planes. Problems with this approach include long scan times, image misregistration, and through-plane motion. This article presents a novel method for calculating cardiac 3D strain using a stack of two or more images acquired in only one orientation. The zHARP pulse sequence encodes in-plane motion using MR tagging and out-of-plane motion using phase encoding, and has been previously shown to be capable of computing 3D displacement within a single image plane. Here, data from two adjacent image planes are combined to yield a 3D strain tensor at each pixel; stacks of zHARP images can be used to derive stacked arrays of 3D strain tensors without imaging multiple orientations and without numerical interpolation. The performance and accuracy of the method is demonstrated in vitro on a phantom and in vivo in four healthy adult human subjects.