Effects on skin dose from unwanted air gaps under bolus in an MR-guided linear accelerator (MR-linac) system

Bolus is commonly used in MV photon radiotherapy to increase superficial dose and improve dose uniformity for treating shallow lesions. However, irregular patient body contours can cause unwanted air gaps between a bolus and patient skin. The resulting dosimetric errors could be exacerbated in MR-Linac treatments, as secondary electrons generated by photons are affected by the magnetic field. This study aimed to quantify the dosimetric effect of unwanted gaps between bolus and skin surface in an MR-Linac. A parallel-plate ionization chamber and EBT3 films were utilized to evaluate the surface dose under bolus with various gantry angles, field sizes, and different air gaps. The results of surface dose measurements were then compared to Monaco 5.40 Treatment Planning System (TPS) calculations. The suitability of using a parallel-plate chamber in MR-Linac measurement was validated by benchmarking the percentage depth dose and output factors with the microDiamond detector and air-filled ionization chamber measurements in water. A non-symmetric response of the parallel-plate chamber to oblique beams in the magnetic field was characterized. Unwanted air gaps significantly reduced the skin dose. For a frontal beam, skin dose was halved when there was a 5 mm gap, a much larger difference than in a conventional linac. Skin dose manifested a non-symmetric pattern in terms of gantry angle and gap size. The TPS overestimated skin dose in general, but shared the same trend with measurement when there was no air gap, or the gap size was larger than 5 mm. However, the calculated and measured results had a large discrepancy when the bolus-skin gap was below 5 mm. When treating superficial lesions, unwanted air gaps under the bolus will compromise the dosimetric goals. Our results highlight the importance of avoiding air gaps between bolus and skin when treating superficial lesions using an MR-Linac system.

Longitudinal acquisition repeatability of MRI radiomics features: An ACR MRI phantom study on two MRI scanners using a 3D T1W TSE sequence

PURPOSE

The purpose of this study was to quantitatively assess the longitudinal acquisition repeatability of MRI radiomics features in a three-dimensional (3D) T1-weighted (T1W) TSE sequence via a well-controlled prospective phantom study.

METHODS

Thirty consecutive daily datasets of an ACR-MRI phantom were acquired on two 1.5T MRI simulators using a 3D T1W TSE sequence. Images were blindly segmented by two observers. Post-acquisition processing was minimized but an intensity discretization (fixed bin size of 25). One hundred and one radiomics features (shape n = 12; first order n = 16; texture n = 73) were extracted. Longitudinal repeatability of each feature was evaluated by Pearson correlation and coefficient of variance (CV_68% ). Interobserver feature value agreement was also quantified using intraclass correlation coefficient (ICC) and Bland-Altman analysis. A most repeatable radiomics feature set on both scanners was determined by feature coefficient of variance (CV_68% <5%), ICC (>0.75), and the ratio of the interobserver difference to the interobserver mean δ<5%.

RESULTS

No trend of radiomics feature value changed with time. Longitudinal feature repeatability CV_68% ranged 0.01-38.60% (mean/median: 12.5%/9.9%), and 0.01-40.47%, (8.49%/7.34%) on the scanners A and B. Shape features exhibited significantly better repeatability than first-order and texture features (all P < 0.01). Significant longitudinal repeatability difference was observed in texture features (P < 0.001) between the two scanners, but not in shape and first-order features (P > 0.30). First-order and texture features had smaller interobserver-dependent variation than acquisition-dependent variation. They also showed good interobserver agreement on both scanners (A:ICC = 0.80 ± 0.23; B:ICC = 0.80 ± 0.22), independent of acquisition repeatability. The repeatable radiomics features in common on both scanners, including 12 shape features, 0 first-order features, and 3 texture features, were determined as the most repeatable MRI radiomics feature set.

CONCLUSIONS

Radiomics features exhibited heterogeneous longitudinal repeatability, while the shape features were the most repeatable, in this phantom study with a 3D T1W TSE acquisition. The most repeatable radiomics feature set derived in this study should be helpful for the selection of reliable radiomics features in the future clinical use.

Reliability of radiomics features due to image reconstruction using a standardized T2 -weighted pulse sequence for MR-guided radiotherapy: An anthropomorphic phantom study

PURPOSE

To prospectively investigate the impact of image reconstruction on MRI radiomics features.

METHODS

An anthropomorphic phantom was scanned at 1.5 T using a standardized sequence for MR-guided radiotherapy under SENSE and compressed-SENSE reconstruction settings. A total of 93 first-order and texture radiomics features in 10 volumes of interest were assessed based on (1) accuracy measured by the percentage deviation from the reference, (2) robustness on reconstruction in all volumes of interest measured by the intraclass correlation coefficient, and (3) repeatability measured by the coefficient of variance over the repetitive acquisitions. Finally, reliable and unreliable radiomics features were comprehensively determined based on their accuracy, robustness, and repeatability.

RESULTS

Better accuracy and robustness of the radiomics features were achieved under SENSE than compressed-SENSE reconstruction. The feature accuracy under SENSE reconstruction was more affected by acceleration factor than direction, whereas under compressed-SENSE reconstruction, accuracy was substantially impacted by the increasing denoising levels. Feature repeatability was dependent more on feature types than on reconstruction. A total of 45 reliable features and 13 unreliable features were finally determined for SENSE, compared with 22 reliable and 26 unreliable features for compressed SENSE. First-order and gray-level co-occurrence matrix features were generally more reliable than other features.

CONCLUSION

Radiomics features could be substantially affected by MRI reconstruction, so precautions need to be taken regarding their reliability for clinical use. This study helps the guidance of the preselection of reliable radiomics features and the preclusion of unreliable features in MR-guided radiotherapy.

Deep Learning-Based Decision-Tree Classifier for COVID-19 Diagnosis From Chest X-ray Imaging

The global pandemic of coronavirus disease 2019 (COVID-19) has resulted in an increased demand for testing, diagnosis, and treatment. Reverse transcription polymerase chain reaction (RT-PCR) is the definitive test for the diagnosis of COVID-19; however, chest X-ray radiography (CXR) is a fast, effective, and affordable test that identifies the possible COVID-19-related pneumonia. This study investigates the feasibility of using a deep learning-based decision-tree classifier for detecting COVID-19 from CXR images. The proposed classifier comprises three binary decision trees, each trained by a deep learning model with convolution neural network based on the PyTorch frame. The first decision tree classifies the CXR images as normal or abnormal. The second tree identifies the abnormal images that contain signs of tuberculosis, whereas the third does the same for COVID-19. The accuracies of the first and second decision trees are 98 and 80%, respectively, whereas the average accuracy of the third decision tree is 95%. The proposed deep learning-based decision-tree classifier may be used in pre-screening patients to conduct triage and fast-track decision making before RT-PCR results are available.

A fast volumetric 4D-MRI with sub-second frame rate for abdominal motion monitoring and characterization in MRI-guided radiotherapy

Background

To propose a fast volumetric 4D-MRI based on 3D pulse sequence acquisition for abdominal motion monitoring and characterization in MRI-guided radiotherapy (MRgRT).

Methods

A 3D spoiled gradient echo sequence volumetric interpolated breath-hold examination (VIBE) [repetition time/echo time (TR/TE) =0.53/1.57 ms, flip-angle =5°, receiver bandwidth (RBW) =1,400 Hz/voxel] based 4D-MRI acquisition, accelerated by 4-fold controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA), named CAIPIRINHA-VIBE 4D-MRI, was implemented on a 1.5T MRI simulator (MR-sim) and applied for abdominal imaging of nine healthy volunteers under free breathing. One hundred and forty-four dynamics of the entire abdomen volume (56 slices), in total 8,064 (144×56) images with a voxel size of 2.7×2.7×4.0 mm³, were acquired in 89 s for 4D-MRI. This CAIPIRINHA-VIBE 4D-MRI was qualitatively compared with a 2D half-Fourier acquisition single-shot turbo spin-echo (2D-HASTE) based 4D-MRI. The motions of liver dome, kidney and spleen were analyzed using the CAIPIRINHA-VIBE 4D-MRI data. The kidney motion was quantitatively characterized in terms of motion range and the correlations between left and right kidneys.

Results

CAIPIRINHA-VIBE 4D-MRI was successfully conducted in all subjects. CAIPIRINHA-VIBE 4D-MRI exhibited much higher effective volumetric temporal resolution (0.615 vs. ~5 s/volume) and better reconstructed volume consistency than 2D-HASTE 4D-MRI. CAIPIRINHA-VIBE 4D-MRI was able to characterize the respiratory motion of abdominal organs simultaneously in three orthogonal directions, and could potentially be used for whole abdomen deformable motion tracking. Renal motion range was most pronounced in superior-inferior (SI) direction (L: 10.03±2.65 mm; R: 10.38±2.80 mm), significantly larger (P<0.001) than that in anterior-posterior (AP) and the least in left-right (LR) directions. Right kidney had significantly larger mobility (4.18±2.19 vs. 2.32±1.34 mm, P=0.045) than left kidney in AP, but not in LR and SI directions. The Pearson correlation coefficients r between left and right kidney motion were 0.5063 (P=0.164), 0.6624 (P=0.052) and 0.5752 (P=0.105) in LR, AP and SI correspondingly. The correlation of renal motion in SI and AP was found significant in right kidney (r=0.843, P=0.004) but not in left kidney (r=0.467, P=0.205).

Conclusions

A fast volumetric 4D-MRI was implemented for abdominal motion monitoring in MRgRT. A sub-second volumetric temporal resolution of 0.615 s, covering the entire abdomen, was demonstrated for respiratory motion monitoring and characterization. This technique holds potentials for MRgRT applications.

A pilot study of highly accelerated 3D MRI in the head and neck position verification for MR-guided radiotherapy

Background

To evaluate the performance of a highly accelerated 3D MRI on inter-fractional positional measurement for MR-guided radiotherapy (MRgRT) in the head and neck (HN).

Methods

Fourteen healthy volunteers received 159 scans on a 1.5 T MR-sim to simulate MRgRT fractions. MRI acquisition included a high-resolution (HQI-MRI, voxel-size =1.05×1.05×1.05 mm³, duration =5 min) and a highly-accelerated low-resolution (true-LQI-MRI, acceleration-factor =9, voxel-size =1.4×1.4×1.4 mm³, duration =86 s) T1w spin-echo sequence (TR/TE =420/7.2 ms). The first session HQI-MRI was used as the reference to mimic planning MRI. Other HQI-MRI was also retrospectively down-sampled in K-space and GRAPPA reconstructed to generate pseudo-LQI-MRI. Inter-sessional positional shift calculated from HQI-MRI, true-LQI-MRI and pseudo-LQI-MRI rigidly registering to the reference were analyzed and compared in the overall HN and the sub-regions of brain, nasopharynx, oropharynx and hypopharynx.

Results

The calculated SD of systematic errors (Σ) from HQI-MRI/pseudo-LQI-MRI/true-LQI-MRI images for overall HN were 1.11/1.14/1.08, 0.28/0.26/0.29, 0.43/0.44/0.60, and 0.77/0.79/0.74 mm for translation in LR, AP, SI and 3D, respectively; The corresponding RMS of random errors (σ) were 0.97/0.98/0.96, 0.28/0.27/0.26, 0.77/0.77/0.72, and 0.85/0.87/0.85 mm. For all sub-regions, brain showed the smallest Σ and σ in 3D. Other sub-regions showed direction-dependent error patterns, but the positioning results were consistent, independent of the datasets used for registration.

Conclusions

A highly-accelerated 3D-MRI could be used for MR-guided HN radiotherapy without compromising position verification accuracy.

A novel method for monitoring the constancy of beam path accuracy in CyberKnife

The aim of current work was to present a novel evaluation procedure implemented for checking the constancy of beam path accuracy of a CyberKnife system based on ArcCHECK. A tailor‐made Styrofoam with four implanted fiducial markers was adopted to enable the fiducial tracking during beam deliveries. A simple two‐field plan and an isocentric plan were created for determining the density override of ArcCHECK in MultiPlan and the constancy of beam path accuracy respectively. Correlation curves for all diodes involved in the study were obtained by analyzing the dose distributions calculated by MultiPlan after introducing position shifts in anteroposterior, superoinferior, and left–right directions. The ability of detecting systematic position error was also evaluated by changing the position of alignment center intentionally. The one standard deviation (SD) result for reproducibility test showed the RMS of 0.054 mm and the maximum of 0.263 mm, which was comparable to the machine self‐test result. The mean of absolute value of position errors in the constancy test was measured to 0.091 mm with a SD of 0.035 mm, while the root‐mean‐square was 0.127 mm with a SD of 0.034 mm. All introduced systematic position errors range from 0.3 to 2 mm were detected successfully. Efficient method for evaluating the constancy of beam path accuracy of CyberKnife has been developed and proven to be sensitive enough for detecting a systematic drift of robotic manipulator. Once the workflow is streamlined, our proposed method will be an effective and easy quality assurance procedure for medical physicists.

Performance evaluation of the CyberKnife system in real-time target tracking during beam delivery using a moving phantom coupled with two-dimensional detector array

The aim of the current study was to evaluate the tracking error of the Synchrony Respiratory Tracking system by conducting beam-by-beam analyses to determine the variation in the tracking beams measured during target motion. A moving phantom of in-house design coupled with a two-dimensional (2D) detector array was used to simulate respiratory motion in the superoinferior (SI) and anteroposterior (AP) direction. A styrofoam block with four implanted fiducial markers was placed on top of the detector to enable the fiducial-based respiratory tracking. Measurements were performed with the phantom under either stationary mode or sinusoidal motion of 6-s cycle and 15/20-mm amplitude at SI and AP direction. The measurement data were saved as movie files that were used to calculate the center shift of the beam with 100-ms sampling time. The tracking accuracy of the system was defined as the targeting error, which could be tracked with probability of > 95% (Ep95). The mean ± standard deviation of Ep95 was 0.28 ± 0.08 mm under stationary condition; 0.66 ± 0.23 mm (range: 0.28-1.22 mm) under sinusoidal respiratory motion. The maximum drift of the beam center for all beam paths was 2.7 mm. The tracking accuracy of CyberKnife Synchrony system was successfully evaluated using a moving phantom and 2D detector array; the maximum tracking error was < 1.5 mm for sinusoidal motion of amplitude ≤ 20 mm.

Assessment of positional reproducibility in the head and neck on a 1.5-T MR simulator for an offline MR-guided radiotherapy solution

Background

Recently, a shuttle-based offline magnetic resonance-guided radiotherapy (MRgRT) approach was proposed. This study aims to evaluate the positional reproducibility in the immobilized head and neck using a 1.5-T MR-simulator (MR-sim) on healthy volunteers.

Methods

A total of 159 scans of 14 healthy volunteers were conducted on a 1.5-T MR-sim with thermoplastic mask immobilization. MR images with isotropic 1.05³ mm³ voxel size were rigidly registered to the first scan based on fiducial, anatomical and gross positions. Mean and standard deviation of positional displacements in translation and rotation were assessed. Systematic error and random errors of positioning in the head and neck on the MR-sim were determined in the translation of, and in the rotation of roll, pitch and yaw.

Results

The systematic error (Σ) of translation in left-right (LR), anterior-posterior (AP) and superior-inferior (SI) direction was 0.57, 0.22 and 0.26 mm for fiducial displacement, 0.28, 0.10 and 0.52 mm for anatomical displacement, and 0.53, 0.22 and 0.49 mm for gross displacement, respectively. The random error (σ) in corresponding translation direction was 2.07, 0.54 and 1.32 mm for fiducial displacement, 1.34, 0.73 and 2.04 mm for anatomical displacement, and 2.24, 0.86 and 2.61 mm for gross displacement. The systematic error and random error of rotation were generally smaller than 1°.

Conclusions

Our results suggested that high gross positional reproducibility (<1 mm translational and <1° rotational systematic error) could be achieved on an MR-sim for the proposed offline MRgRT.

Development of a novel methodology for QA of respiratory-gated and VMAT beam delivery using Octavius 4D phantom

The objective of this study was to develop and evaluate a series of quality assurance (QA) techniques based on Octavius 4D phantom for testing of respiratory-gated treatment delivery, integrity of dose rate vs gantry speed in volumetric-modulated arc therapy (VMAT) commissioning, and multileaf collimator (MLC) positioning accuracy of a linear accelerator. An Octavius 4D phantom capable of rotating with the gantry and recording the detector signal with a sampling rate of 10 Hz was isocentrally set up and an inclinometer was also installed to measure the gantry angle simultaneously. A simple arc test was created and delivered with gating function activated to measure the timing accuracy of the gating window. A tailor-made dose rate vs gantry speed plan was also designed to test the accuracy of measured dose rate, gantry speed, and actual control points. All experiments were conducted while machine log files were collected for comparison. The variations of beam flatness, symmetry, and field size were analyzed as a function of gantry angle to evaluate the influence from the modulation of dose rate and gantry speed. MLC position accuracy was evaluated based on specific garden fence plans. The time of gating window was measured to be less than 10-millisecond deviation from the log data. Gantry backlash was observed and quantified to be 1.72° with an extra stabilization time of 1.16 seconds for a gating arc with gantry speed of 6°/s. In the dose rate vs gantry speed test, the mean deviation between measured gantry angle and log data was less than 0.2° after a time delay of 0.25 second was corrected. The measured dose rate agreed with the log data very well with a mean deviation of 0.05%, and even the transit of modulation was tracked successfully. There was a statistically significant difference on the variation of beam parameters between a VMAT plan and a simple arc plan. The induced MLC position errors were detected with an accuracy of 0.05 mm. The leaf position reproducibility was found to be better than 0.02 mm, whereas the routine MLC position accuracy was better than 0.1 mm. A time-resolved method using Octavius 4D phantom has been developed and proven to be convenient for respiratory gating QA, dose rate vs gantry speed test, and MLC QA. Gating time, dose rate, and gantry speed-induced leave position error could be directly measured with high accuracy after comparison with the machine log data. This study also highlights the capability of the phantom in quantifying the variation of flatness, symmetry, and field size during gantry rotation.

Correction to: Molecular characterization of Blastocystis sp. in captive wild animals in Qinling Mountains

There were errors in Fig. 1 of the originally published article. Correct fig. 1 is presented here.

Image quality assessment of a 1.5T dedicated magnetic resonance-simulator for radiotherapy with a flexible radio frequency coil setting using the standard American College of Radiology magnetic resonance imaging phantom test

BACKGROUND

A flexible RF coil setting has to be used on an MR-simulator (MR-sim) in the head and neck simulation scan for radiotherapy (RT) purpose, while the image quality might be compromised due to the sub-optimized flexible coil compared to the normal diagnostic radiological (DR) head coil. In this study, we assessed the image quality of an MR-sim by conducting the standard American College of Radiology (ACR) MRI phantom test on a 1.5T MR-sim under RT-setting and comparing it to DR-setting.

METHODS

A large ACR MRI phantom was carefully positioned, aligned and scanned 9 times for each under RT- and DR-setting on a 1.5T MR-sim, following the ACR scanning instruction. Images were analyzed following the ACR guidance. Measurement results under two coil settings were quantitatively compared. Inter-observer disagreements under RT-setting between two physicists were compared using Bland-Altman (BA) analysis and intra-class correlation coefficient (ICC).

RESULTS

The MR-sim with RT-setting obtained sufficiently good image quality to pass all ACR recommended criteria. No significant difference was found in phantom length accuracy, high-contrast spatial resolution, slice thickness accuracy, slice position accuracy, and percent-signal ghosting. RT-setting significantly under-performed in low-contrast object detectability, while better performed in image intensity uniformity. BA analysis showed that 95% limit of agreement and biases of phantom test measurement under RT-setting between two observers were very small. Excellent inter-observer agreement (ICC >0.75) was achieved in all measurements except for slice thickness accuracy (ICC =0.42, moderate agreement) under RT-setting.

CONCLUSIONS

Very good and highly reproducible image quality could be achieved on a 1.5T MR-sim with a flexible coil setting as revealed by the standard ACR MRI phantom test. The flexible RT-setting compromised in image signal-to-noise ratio (SNR) compared to the normal DR-setting, and resulted in reduced low-contrast object detectability.

Filmless methods for quality assurance of Tomotherapy using ArcCHECK

PURPOSE

Tomotherapy delivers an intensity-modulated radiation therapy (IMRT) treatment by the synchronization of gantry rotation, multileaf collimator (MLC), and couch movement. This dynamic nature makes the quality assurance (QA) important and challenging. The purpose of this study is to develop some methodologies using an ArcCHECK for accurate QA measurements of the gantry angle and speed, MLC synchronization and leaf open time, couch translation per gantry rotation, couch speed and uniformity, and constancy of longitudinal beam profile for a Tomotherapy unit.

METHODS

Four test plans recommended by AAPM Task Group 148 (TG148) and the manufacturer were chosen for this study. Helical and static star shot tests are used for checking the leaves opened at the expected gantry angles. Another helical test is to verify the couch traveled the expected distance per gantry rotation. The final test is for checking the couch speed constancy with a static gantry. ArcCHECK can record the detector signal every 50 ms as a movie file, and has a virtual inclinometer for gantry angle measurement. These features made the measurement of gantry angle and speed, MLC synchronization and leaf open time, and longitudinal beam profile possible. A shaping parameter was defined for facilitating the location of the beam center during the plan delivery, which was thereafter used to calculate the couch translation per gantry rotation and couch speed. The full width at half maximum (FWHM) was calculated for each measured longitudinal beam profile and then used to evaluate the couch speed uniformity. Furthermore, a mean longitudinal profile was obtained for constancy check of field width. The machine trajectory log data were also collected for comparison. Inhouse programs were developed in MATLAB to process both the ArcCHECK and machine log data.

RESULTS

The deviation of our measurement results from the log data for gantry angle was calculated to be less than 0.4°. The percentage differences between measured and planned leaf open time were found to be within 0.5% in all the tests. Our results showed mean values of MLC synchronization of 0.982, 0.983, and 0.995 at static gantry angle 0°, 45°, and 135°, respectively. The mean value of measured couch translation and couch speed by ArcCHECK had less than 0.1% deviation from the planned values. The variation in the value of FWHM suggested the couch speed uniformity was better than 1%. The mean of measured longitudinal profiles was suitable for constancy check of field width.

CONCLUSION

Precise and efficient methods for measuring the gantry angle and speed, leaf open time, couch translation per gantry rotation, couch speed and uniformity, and constancy of longitudinal beam profile of Tomotherapy using ArcCHECK have been developed and proven to be accurate compared with machine log data. Estimation of the Tomotherapy binary MLC leaf open time is proven to be precise enough to verify the leaf open time as small as 277.8 ms. Our method also makes the observation and quantification of the synchronization of leaves possible.

Technical note: Correlation between TQA data trends and TomoHD functional status

TomoTherapy Quality Assurance (TQA) is a software package developed to monitor certain aspects of machine performance. In this study, the TQA quantities or data trends most effective in monitoring energy drifts and magnetron stability were determined respectively. This retrospective study used data collected from three TomoHD units. The TQA modules investigated were Step‐Wedge Helical, Step‐Wedge Static, and Basic Dosimetry. First, the TQA quantities correlated with energy changes (|r|>0.85, where r is the Pearson's correlation coefficient) were found. The corresponding sensitivities to percentage depth dose (PDD) ratio changes were then calculated and compared. Second, the pulse‐by‐pulse dose stability was compared before and after each magnetron replacement using a nonparametric comparison test (Welch's t‐test), and the raw dose profiles were surveyed. In this study, exit detector flatness obtained in Basic Dosimetry was shown to be the most sensitive (r=0.945) to energy changes, followed by the energy differences in Step‐Wedge Static (r=0.942) and Step‐Wedge Helical (r=0.898). The three quantities could detect a PDD ratio change of 5.1×10−4,5.4×10−4, and 7.1×10−4, respectively. Pulse‐by‐Pulse Dose1 from Basic Dosimetry over a one‐week period before and after a magnetron replacement showed a significant difference (p<0.05) in only three of the nine instances. On the other hand, a raw output profile free from discontinuities, frequent dropped pulses and abnormal spikes was found to indicate that the magnetron would continue to function normally for a week 89% of the time.

PACS numbers: 87.56.bd, 87.56.Fc, 84.40.Fe

Accuracy in contouring of small and low contrast lesions: comparison between diagnostic quality computed tomography scanner and computed tomography simulation scanner-A phantom study

To evaluate the accuracy in detection of small and low-contrast regions using a high-definition diagnostic computed tomography (CT) scanner compared with a radiotherapy CT simulation scanner. A custom-made phantom with cylindrical holes of diameters ranging from 2-9 mm was filled with 9 different concentrations of contrast solution. The phantom was scanned using a 16-slice multidetector CT simulation scanner (LightSpeed RT16, General Electric Healthcare, Milwaukee, WI) and a 64-slice high-definition diagnostic CT scanner (Discovery CT750 HD, General Electric Healthcare). The low-contrast regions of interest (ROIs) were delineated automatically upon their full width at half maximum of the CT number profile in Hounsfield units on a treatment planning workstation. Two conformal indexes, CI(in), and CI(out), were calculated to represent the percentage errors of underestimation and overestimation in the automated contours compared with their actual sizes. Summarizing the conformal indexes of different sizes and contrast concentration, the means of CI(in) and CI(out) for the CT simulation scanner were 33.7% and 60.9%, respectively, and 10.5% and 41.5% were found for the diagnostic CT scanner. The mean differences between the 2 scanners' CI(in) and CI(out) were shown to be significant with p < 0.001. A descending trend of the index values was observed as the ROI size increases for both scanners, which indicates an improved accuracy when the ROI size increases, whereas no observable trend was found in the contouring accuracy with respect to the contrast levels in this study. Images acquired by the diagnostic CT scanner allow higher accuracy on size estimation compared with the CT simulation scanner in this study. We recommend using a diagnostic CT scanner to scan patients with small lesions (<1 cm in diameter) for radiotherapy treatment planning, especially for those pending for stereotactic radiosurgery in which accurate delineation of small-sized, low-contrast regions is important for dose calculation.

Characterization of Baylisascaris schroederi from Qinling subspecies of giant panda in China by the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA

In the present study, a total of 20 nematode isolates, (including 10 male and 10 female worms) representing Baylisascaris schroederi from 5 Qinling subspecies of giant pandas (Ailuropoda melanoleuca) in Shaanxi Province of China, were characterized and grouped genetically by the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA (rDNA). The rDNA fragment spanning 3' end of 18S rDNA, complete ITS-1 rDNA, and 5' end of 5.8S rDNA were amplified and sequenced. The sequence variability in ITS-1 rDNA was examined within B. schroederi and among parasites in order Ascaridata available in GenBank™, and their phylogenetic relationships were also reconstructed. The sequences of ITS-1 rDNA for all the B. schroederi isolates were 427 bp in length, with no genetic variation detected among these isolates. Phylogenetic analyses based on the ITS-1 rDNA sequences revealed that all the male and female B. schroederi isolates sequenced in the present study were posited into the clade of genus Baylisascaris, sistered to zoonotic nematodes in genus Ascaris, and the ITS-1 rDNA sequence could distinguish different species in order Ascaridata. These results showed that the ITS-1 rDNA provides a suitable molecular marker for the inter-species phylogenetic analysis and differential identification of nematodes in order Ascaridata.

Rarefied-gas heat transfer in micro- and nanoscale Couette flows

The physics of the heat conduction and viscous dissipation in rarefied gases is analyzed and discussed. A heat transfer model valid for arbitrary Knudsen numbers, defined as the ratio of the molecular mean free path to the characteristic length of channels, is derived by treating the heat transfer behavior in the slip and transition regimes as an intermediate function of continuum heat transfer model and free molecular heat transfer model. Comparison studies reveal that this model not only shows good agreement with the numerical results based on the direct simulation Monte Carlo method, but also has some unique features that can overcome the deficiencies existing in the previous models. Therefore, this model is capable to study the heat transport phenomena in very dilute gas Couette flows through micro/nanochannels more accurately.

Resistance to adefovir dipivoxil in lamivudine resistant chronic hepatitis B patients treated with adefovir dipivoxil

BACKGROUND

Adefovir dipivoxil (ADV) is a potent nucleotide analogue against both the wild-type and lamivudine (LMV) resistant hepatitis B virus (HBV). The cumulative incidence of ADV resistant mutations in the nucleoside/-tide treatment naive chronic hepatitis B patient (CHB) at weeks 48, 96, and 144 was 0, 0.8-3%, and approximately 5.9%, respectively.

AIMS

The aim of this study was to characterise the genotypic and phenotypic mutation profiles to ADV in 67 LMV resistant CHB patients who were treated with ADV.

METHODS

Serum HBV DNA was quantified by real time polymerase chain reaction. The ADV mutant was detected using matrix assisted laser desorption/ionisation time of flight mass spectrometry based genotyping assays, termed restriction fragment mass polymorphism (RFMP).

RESULTS

RFMP analysis revealed that a total of 11 amino acid substitutions developed in the rt domain of the HBV polymerase in nine patients. The cumulative incidence of genotypic ADV resistance at months 12 and 24 was 6.4% and 25.4%, respectively. The rtA181V, rtN236T, and rtA181T mutations were detected in five, four, and two of the 67 patients at treatment months 12-17, 3-19, and 7-20, respectively. Serial quantification of serum HBV DNA revealed that two patients with the rtA181V mutation, with or without the rtN236T mutation, and one patient with the rtA181T mutation displayed HBV DNA rebound.

CONCLUSION

Emergence of the ADV mutation in LMV resistant patients who are treated with ADV appeared to present earlier and more frequently than was reported in previous studies on nucleoside/-tide treatment naive patients.

MAPD--an objective way to select mAs for paediatric brain CT

CT is an advanced imaging modality, but the imaging parameters are normally selected subjectively. For standard head examinations, most of the parameters used are consistent amongst different centres, with the exception of large variations in the selection of the tube current-exposure time product (mAs). As a result, CT images may contain unacceptable levels of noise, or the patient may receive excessive radiation. In this study, the maximum anteroposterior diameter (MAPD) was shown to be a good criterion for mAs selection, and could be measured in a pilot view. 200 paediatric brain CT studies were randomly selected to determine the MAPD at the mid brain level. With knowledge of MAPD distribution, a phantom study was performed to determine the relationship between MAPD and the mAs required for consistent and acceptable image noise. It was found that the required mAs increased linearly with MAPD. Assuming the manufacturer's recommended value is "appropriate" for the average MAPD, the appropriate mAs value could be estimated. Using this method, appropriate mAs values were calculated retrospectively for a group of 240 randomly selected paediatric brain CT studies and compared with the actual mAs subjectively determined by the radiographer. Although their average values were similar, the difference between the calculated and actual values deviated markedly in some cases. When the actual mAs was smaller than the calculated value, higher image noise was observed. However, reduction of image noise was barely observed when the applied mAs was larger than the calculated value. Thus, this method is more objective and appropriate for determination of the mAs value for paediatric brain CT than the traditional subjective method.