Imaging Lung Function in Mice Using SPECT/CT and Per-Voxel Analysis

Chronic lung disease is a major worldwide health concern but better tools are required to understand the underlying pathologies. Ventilation/perfusion (V/Q) single photon emission computed tomography (SPECT) with per-voxel analysis allows for non-invasive measurement of regional lung function. A clinically adapted V/Q methodology was used in healthy mice to investigate V/Q relationships. Twelve week-old mice were imaged to describe normal lung function while 36 week-old mice were imaged to determine how age affects V/Q. Mice were ventilated with Technegas™ and injected with 99mTc-macroaggregated albumin to trace ventilation and perfusion, respectively. For both processes, SPECT and CT images were acquired, co-registered, and quantitatively analyzed. On a per-voxel basis, ventilation and perfusion were moderately correlated (R = 0.58±0.03) in 12 week old animals and a mean log(V/Q) ratio of −0.07±0.01 and standard deviation of 0.36±0.02 were found, defining the extent of V/Q matching. In contrast, 36 week old animals had significantly increased levels of V/Q mismatching throughout the periphery of the lung. Measures of V/Q were consistent across healthy animals and differences were observed with age demonstrating the capability of this technique in quantifying lung function. Per-voxel analysis and the ability to non-invasively assess lung function will aid in the investigation of chronic lung disease models and drug efficacy studies.


Introduction
Chronic respiratory diseases such as asthma and chronic obstructive pulmonary disease (COPD) are burdensome on healthcare systems and are associated with substantial morbidity and mortality worldwide [1]. Methodologies for three-dimensional imaging of the lung are important as they allow for better characterization and understanding of these obstructive diseases [2]. While volumetric CT densitometry is now widely used, it generally contributes only structural information while other imaging techniques can provide functional information such as the contribution of alveolar ventilation (V) and pulmonary capillary perfusion (Q) to gas exchange within the lung [3]. There are non-imaging methods of V/Q assessment [4], and the forced-expiration pulmonary function tests (PFTs) currently used for diagnosis of obstructive diseases are relatively simple to perform, but these methods cannot evaluate lung function in a regional manner. In addition, PFTs are unable to differentiate pathologies and are insensitive to early disease [5]. While imaging methods are technically more difficult than PFTs they can provide both global and regional functional information. Further, these volumetric techniques are necessary to address the heterogeneous underlying pathologies that make up diseases such as COPD [2,6,7].
In the field of clinical nuclear medicine, static V/Q imaging is routinely performed with single photon emission computed tomography (SPECT) along with Technegas TM and 99m Techetium-labelled macroaggregated albumin (MAA) to provide the distributions of V and Q, respectively. Images acquired by this method depict the relative contributions of the core pulmonary processes of V and Q within any particular region in the lung. When V/Q images are co-registered to CT, regional analysis can be coupled to structural information. Clinical research examining the pathophysiology of COPD has been shown to benefit significantly from the use of radiologist-scored V/Q SPECT/CT [6] but investigation into the use of per-voxel volumetric analysis has only just begun in this area [8,9].
Pre-clinical models illustrating the utility of SPECT/CT as a tool in respiratory research are currently scarce [10]. Noninvasive assessment of small animal models offers an approach to address a natural link between discoveries at the molecular level and the application of clinically relevant diagnostics or therapeutics [10,11,12]. 3D imaging has now been used in various pulmonary disease models including asthma and emphysema, among others [10,13,14]. In regards to preclinical functional lung imaging, V/Q scanning has only seen service in rats through use of magnetic resonance imaging (MRI) and has so far lacked whole lung assessment [15,16]. Global lung measurements of V/Q have also been made in rats using the multiple inert gas elimination technique [17], but as this methodology does not employ any imaging modality it cannot offer regional analysis. Although rats offer better resolution of lung structure due to their size, the ability to test molecular hypotheses is greater in mice due to the relative abundance of knockout strains and commercially available tests. The mouse has yet to be investigated by any V/Q methodology and no extensive quantitative V/Q evaluation of lung function has been performed in a murine model.
Utilizing V/Q SPECT/CT to measure functional consequences of known pathologies in murine models could improve our understanding of the pathophysiology of respiratory disease and translate knowledge gained to the clinical setting, where similar imaging methodologies can be applied in most nuclear medicine departments. In this study we have successfully adapted a common clinical V/Q imaging protocol, employing Technegas TM and 99m Tc-MAA, to mice. In addition, we have applied a threedimensional quantitative per-voxel analysis to characterize the naïve lung environment. Further, the sensitivity of this lung function measurement was tested by global and regional assessment of age-related changes. This work provides a foundation for further investigation of lung function in respiratory disease models.

Animals
Female BALB/c mice (Charles River, QC, Canada) were acclimatised to specific pathogen-free housing conditions for a period of 2-3 weeks, with a 12:12 hour light:dark cycle, prior to imaging. Experimentation was performed on 17 mice (12-16 weeks old, henceforth referred to as 12 w.o.) while another group of mice were housed for a further 24 weeks before imaging (henceforth referred to as 36 w.o.). The experiments described in this study were approved by the Animal Research Ethics Board of McMaster University (Animal Utilization Protocol #080514) in accordance with the guidelines of the Canadian Council on Animal Care.

Imaging Protocol
Healthy, untreated mice were tracheally intubated through the oral cavity with a 20-gauge catheter following anaesthetisation with ketamine/xylazine (90 mg/kg, 6 mg/kg). Animals were then ventilated (0.02 L/min, 125 strokes/min) on a Rodent Ventilator (Model 683, Harvard Apparatus, Holliston, MA, USA) with Technegas TM (Cyclomedica, Lucas Heights, NSW, Australia) (0.04-0.12 MBq/mL) for 15 minutes, alternating between Technegas TM and room air in intervals of 15-30 seconds. Technegas TM is an ultrafine aerosol composed of carbon particles labelled with 99m Tc and suspended in argon, with a mean aerosol diameter between 30 and 60 nm [18]. Following ventilation mice were moved and strapped securely to a heated bed for animal welfare and consistency in body temperature between animals. SPECT scans were acquired on an X-SPECT system (Gamma Medica, Northridge, CA, USA) using dual sodium iodide crystals in combination with low energy pinhole collimators with 1 mm aperture and a radius of rotation of 3.5 cm. The ventilation SPECT scan consisted of thirty-two 50 second projections and was followed immediately by the collection of four rotations of 1024 Xray projections for CT, also acquired on the X-SPECT system with x-ray tube characteristics of 75kVp and 220 mA. Following ventilation SPECT and CT scans, mice were injected with 11-15 MBq of 99m Tc-MAA via the tail vein. Care was taken not to shift the position of the animal during the injection. Perfusion imaging entailed a 1024-projection CT scan and a SPECT scan of thirty-two 40 second projections. Supplemental gaseous anaesthesia was used (isoflurane (1%, 1 L/min)) to ensure sedation throughout the imaging procedure (approximately 1 hour 20 minutes).
Mice were recovered from the imaging procedure and monitored until radioactivity had decayed to background levels. All imaging work was completed at the McMaster Centre for Preclinical and Translational Imaging (MCPTI) at McMaster University (Hamilton, ON, Canada). A simplified protocol, detailing the sequence of imaging methodology, is shown in Figure 1. Reconstruction of 3D SPECT/CT Images SPECT images for both ventilation and perfusion scans were produced by iterative reconstruction of their respective projections using FLEX TM -SPECT software (Gamma Medica) into 80680680 arrays (0.472 mm isotropic voxels). The projections from the four ventilation CT rotations were first summed, then images were reconstructed using a Feldkamp cone beam backprojection algorithm in COBRA (Exxim Sofware, Pleasanton, CA, USA) into 51265126512 arrays (0.115 mm isotropic voxels). This high quality CT image was used in the processing and analysis of V/Q data to provide anatomical features for SPECT coregistration and to also allow for a density-based assessment of the lung environment. The perfusion CT projections were reconstructed under the same conditions. Calibration of each CT image for Hounsfield Unit (HU) scaling was performed using empty airspace within the field of view and a water-filled tube included in each scan.

Fusion and Co-registration of Ventilation and Perfusion SPECT/CT Data
Ventilation and perfusion CT images (512 3 matrix) were first subjected to a Gaussian filter (s = 0.8) using Matlab version 7.5 (Mathworks, Natick, MA, USA). A 'Lung' label was produced for the CT images as described previously for rats [13] using Amira 5.1 software (Visage Imaging, Andover, MA, USA), with the exception that the trachea and main bronchi (up to the point where they enter lung tissue) were removed from the ventilation lung label. Briefly, a 250HU threshold-limited growing algorithm was used to define the airspace. This selection was then subjected to a simple volume-growing algorithm and another threshold of 100HU applied to exclude tissues and structures not belonging to the lung. A hole-filling algorithm was also applied to include internal structures such as major blood vessels and, finally, the label was manually examined to ensure the process was consistent across animals. CT images and lung labels were then compressed to a matrix of 256 3 for computational purposes. Representative axial, coronal, and sagittal images for the CT, lung label, ventilation, and perfusion scans can be seen in Figure 2 while a representative lung label is shown in Figure 3A.
The perfusion CT and SPECT were fused, with in-house software developed in Matlab, by a process that maximised mutual information (MI) within the lung region as defined by the label field [19]. Powell's multidimensional direction set method was used to maximise MI using a one dimensional search algorithm based on golden section search and parabolic interpolation [20]: SPECT images underwent rigid body transformation until a change of less than 0.01 mm (translation) or 0.01 degrees (rotation) was observed along or around each axis to obtain fusion parameters. During the process, the 80 3 matrix of the SPECT image was interpolated for comparison to the CT image and was resampled to a 256 3 matrix when the stop criteria were met; this corresponds to a change in voxel size from 0.40 mm 3 to 0.23 mm 3 . When maximised MI was reached the result was visually inspected to confirm fusion quality. Parameters obtained from fusion of the perfusion SPECT and CT were then applied to the ventilation SPECT and CT as the spatial relationship between these two data sets remains constant.
The ventilation and perfusion data sets were then co-registered by repeating the MI maximisation process through rigid body transformation of their respective CT images, again within the lung segmentation. Specifically, the perfusion CT was adjusted to fit the ventilation CT. These co-registration parameters were then applied to the perfusion SPECT data allowing for comparison of ventilation and perfusion SPECT images. A simplified representation of image processing is shown in Figure 1.

Quantitative Per-Voxel Assessment
Co-registration of images allows for analysis of SPECT and CT data on a per-voxel basis. Within the lung label all voxels had a HU, V, Q, and log(V/Q) value; presented data represent distributions of all voxels within the lung label averaged across all animals, unless otherwise stated. V and Q values were converted to relative frequencies by dividing the activity value of each voxel by the total activity in the lung label and then multiplying by 100. The relative frequency of both V and Q were used to calculate V/Q ratios and a base 10 logarithm was applied to provide a log-normal distribution. Analysis of the log(V/Q) data involved calculation of the mean, standard deviation, skewness, excess kurtosis, and percentage of total lung volume (TLV) greater than 62 averaged standard deviations from the averaged mean log(V/Q), as a measure of V/Q mismatching. Any voxel where only V equalled zero was set to a log(V/Q) ratio of 2' and any voxel where only Q equalled zero was set to a log(V/Q) ratio of +'; these values were not included in the calculation of the log(V/Q) mean, standard deviation, skew or kurtosis. Any voxels where both V and Q equalled zero were noted and given a log(V/Q) value of zero.

Regionalization of Images
Two different methods of regionalization were performed on V, Q, and log(V/Q) data. Apex, middle, and base divisions were made by finding the axial slices closest to the 33 and 66 percentiles of volume. Inner and outer divisions were made by eroding the outer boundary of the lung label by 3 voxels, using a six-neighbour structuring element in a volumetric manner, and subtracting the new inner data from the original whole lung to produce the outer data. Standard deviations of log(V/Q) were calculated for all voxels in a region in the same manner as described above.

Respiratory-gated CT Image Processing & Quantitation
Inspiratory-and expiratory-gated CT images were produced by applying RespGate software (RespGate, Hamilton, ON, Canada) [21] to all ventilation CT projections. Respiratory-gated images were then reconstructed and calibrated as described for the ungated image used in the V/Q analysis. Air volume information from respiratory-gated CTs was extracted by a method previously described in rats [13]. Briefly, a segmentation was made encompassing the thoracic cavity using Amira 5.1 software and a histogram of the densities (25HU/bin) within this region of interest was produced. The volume of air within the thoracic segmentation was calculated from the histogram by multiplying the number of voxels in a HU density bin against an air coefficient for that HU bin and converting the result into millilitres. The air coefficient was defined by the centre of the HU bin and represented the fraction of air in a voxel of that value; e.g., the air coefficient for 2500 HU is 0.5. Summation of the airhistogram from the end expiratory image allowed for quantitation of functional residual capacity (FRC) and tidal volume (TV) was calculated as the difference between air volume in the inspiratory and expiratory lung states.

Data Analysis
Images were analyzed and all measures were output with the aid of Matlab. Boxplots represent the lower quartile, median and upper quartile while whiskers represent the minimum and maximum values within the data set. Data are expressed as the mean6SEM and, where applicable, statistical significance was determined by unpaired two-tailed t-test in Graphpad Prism (Graphpad Software Inc, La Jolla, CA, USA); p,0.05 was considered statistically significant. Two of the seventeen 12 w.o. animals were excluded for the following reasons: one V/Q data set was a statistical outlier for the ratio of total perfusion activity to total ventilation activity, and another animal received an improper intravenous delivery of MAA, as noted at the time of injection. One 36 w.o. animal was also excluded as an outlier for increased FRC volume.

Results
Many lung function parameters can be measured from V/Q SPECT/CT: Ventilation and perfusion distribution within the lung can be determined; the relationship between these processes can be studied by measuring their correlation on a per-voxel basis; the degree of V/Q mismatching is illustrated in the log(V/Q) distribution and can be quantified using thresholds, such as two standard deviations from the mean; from the CT, tidal volume and FRC can be quantified in addition to changes in lung density. Investigating these parameters provides a detailed evaluation from which lung dysfunction can be assessed.
Typically, the distribution of deposition in the ventilation images was somewhat heterogeneous and typified by areas of increased intensities, predominately in apical regions, while decreased intensities were observed in basal regions (Figure 2). The lung perfusion distribution demonstrated greater uniformity of counts with less apparent skew towards the upper regions of the lung. Similarly, log(V/Q) images described a relatively homogeneous distribution with only small regions of obvious V/Q mismatch. The average total number of counts from ventilation with Technegas TM , measured within the lung, was 1.9660.33610 5 while the average total number of counts from perfusion with 99m Tc-MAA was 3.6160.31610 6 (data not shown) creating a mean ratio of 24.263.7 to 1 between perfusion and ventilation total activity. This ratio was unaltered in the 36w.o. animals.
The degree of ventilation and perfusion matching in the lung was determined by plotting the distribution of log(V/Q) values ( Figure 4A). The mean and standard deviation of log(V/Q) values were narrowly distributed with averages of 20.07160.006 and 0.3660.02, respectively ( Figure 4A & B). In addition, the log(V/Q) distribution was not normal in nature as it was negatively skewed and positively kurtotic. Approximately 6.6% of TLV had a low log(V/Q) ratio mismatch while approximately 1.0% of TLV had a high log(V/Q) ratio mismatch ( Figure 4A). Of the low log(V/Q) mismatched volume 3.2360.75% of TLV was made up of voxels with a log(V/Q) value of 2' but no significant proportion of voxels with values of +' were found ( Figure 4C).
In the comparison of 12w.o. to 36w.o. animals, discrete histograms were generated for both ventilation and perfusion in order to assess any difference due to age in these physiological processes; 36w.o. animals demonstrated a shift towards lower relative count values for both ventilation and perfusion when compared to their younger counterparts ( Figure 5A). The leftward shift of these distributions was not apparent when these processes were translated to a per-voxel log(V/Q) ratio, where the difference between 12w.o. and 36w.o. animals is much more subtle ( Figure 5B). However, older animals had significant changes in standard deviation, skewness, and excess kurtosis of the log(V/Q) distribution ( Figure 5C). Furthermore, values for positive and negative ' were significantly increased along with values for high log(V/Q) ( Figure 5D); these changes were larger in magnitude for negative values than positive.
Log(V/Q) images of 12 and 36w.o. animals indicated that V/Q mismatching was distributed throughout the lung ( Figure 6A). Regionalization of lung volume into apex, middle, and base divisions (approximately 33.3% of volume per division) as well as inner and outer divisions (approximately 45 and 55% of volume per division, respectively) allowed for further analysis of the distribution of V/Q heterogeneity and mismatch; no change in regional volume proportion was observed with age. In 12w.o. animals, the standard deviation of log(V/Q) values was greater in the basal region of the lung compared to the apex or middle regions ( Figure 6B). In 36w.o. animals increases in standard deviation over 12.w.o. values were observed in all three regions. Inner and outer regionalization demonstrated that the increase in standard deviation values with age was greater in the periphery of the lung ( Figure 6B). This age-related difference in the outer region was associated with a significant decrease in the percentage of both ventilation and perfusion counts in that region (data not shown), indicating that the increase in V/Q heterogeneity is associated with alterations in the distribution of both processes. This regional distribution of increased log(V/Q) heterogeneity indicates that mismatching is present throughout the lung and is primarily associated with the periphery.
The CT density distribution of the lung was described by a bimodal curve with a primary peak at approximately 2325HU, describing the air volume, and a lesser peak at approximately 225HU, describing tissue volumes within or surrounding the lung ( Figure 3B). Post-acquisition respiratory gating produced CT images of end inspiratory and end expiratory lung states. These curves are unimodal in nature due to the calculation of air content based on HU. There were notable differences in density distributions between inspiratory and expiratory states including a broadening of the peak and a shift towards lower HU values during inspiration ( Figure 7A). These curves were further shifted towards lower densities in older mice. For 12w.o. animals, FRC was determined to be 0.1260.01 mL while TV was determined to be 0.04860.004 mL. FRC was significantly increased in the 36w.o. group while TV remained constant ( Figure 7B). When FRC, denoting lung volume, was plotted against the median value for ventilation and perfusion distributions a high degree of correlation (R 2 of 0.77, p,0.0001 and 0.75, p,0.0001 respectively) was observed for all 12w.o. animals indicating that the leftward shift in these distributions in the 36w.o. animals was due to the increased FRC values in this group. However, no correlation was seen when FRC was plotted against the standard deviation of the log(V/Q) curve, as a general measure of V/Q heterogeneity, indicating that FRC is not a driving factor in the log(V/Q) mismatch observed in older animals. In addition to the increase in FRC, a small increase in body mass was observed in 36w.o. mice (Table 1).

Discussion
The capability to study the core processes of a vital organ, namely ventilation and perfusion in the lung, at a scale relevant to these processes represents an imperative step towards understanding lung function and how it changes in respiratory disease. This study demonstrates that we have successfully translated a common clinical V/Q imaging technique to the preclinical setting to measure lung function in mice. Using V/Q SPECT/CT as a biomarker of lung function will support preclinical investigations of the pathogenesis of respiratory diseases and screening of novel drug therapies, helping to bridge the gap to the clinical understanding of disease and drug efficacy. In addition, we have developed a preclinical method of image co-registration and data analysis to objectively evaluate V/Q at a per-voxel level so that regional assessments of lung function can be performed.
A V/Q ratio of 0.80, or a log(V/Q) of 20.097, is considered optimal for a healthy human lung [22]. We found a mean log(V/ Q) of 20.07160.006, corresponding to a V/Q ratio of 0.85, in our naïve mice indicating that healthy mammalian lungs have similar respiratory physiology. These values represent the gas exchange capability for the entire lung, but ventilation and perfusion, at a small scale, are somewhat heterogeneous processes. While early work suggested that gravity was a large determinant of heterogeneity in V/Q ratios [23], methodologies with better resolution describe greater variability at any particular axial height within the lung for both ventilation and perfusion [24]. It is now thought that the majority of this intrinsic variability can be attributed to factors related to the basic architecture of the lung [25]. To define the natural degree of heterogeneity in this model V/Q mismatching was defined as the percentage of total lung volume (TLV) with log(V/Q) ratios 62 standard deviations from the mean, as 95% of values should fall within this range for a normal distribution. For this purpose the average mean and standard deviation of the total data set was used with the notion that severely mismatched voxels would fall beyond the scope of this average ''normally functioning'' region. Only a small percentage of volume (,1%) was found above +2SD while a higher proportion (,7%) was found below 22SD. This population included voxels where ventilation values, but not perfusion values, equalled zero and thus had a log(V/Q) value of 2'. We found that approximately 3% of TLV had such a log(V/ Q). Most of these zero ventilation values appear at the perimeter of the lung-label, especially around the caudal-most edges, and likely represent tissues not associated with alveolar regions or marginal errors in the co-registration process that were included in the lung label. Such voxels that are non-lung could contribute to the 6' categories; however, since the segmentation employed is a standardized process any sampling errors between animals should be consistent. Overall, the variability in V/Q distribution within the naïve animals studied was low, with a log(V/Q) standard deviation of 0.3660.02, and therefore represents a typical degree of ventilation and perfusion heterogeneity in this system.
In humans, a decline in V/Q matching is known to occur with increased age [26,27,28] so we elected to study the effect of aging on lung function to determine the sensitivity of our methodology.
To provide an experimental model to this effect a subset of mice were allowed to age for approximately 6 months, which is a considerable amount of time in the life of a mouse. Though the mice at 36 weeks of age had significantly more mass, and obesity is known to impact lung function [29], these mice were not considered obese and their weight would have little, if any, impact on V/Q matching (R = 0.1426, p = 0.79). However, we found that older animals had a significant increase in V/Q mismatching, as defined by the standard deviation of the log(V/Q) data, over their younger counterparts suggesting that this measure is sensitive to minor changes in lung function. The greater FRC observed in 36w.o. animals is likely associated with a reduction in lung elasticity, which is in turn due to a general loss of alveolar surface area [26]. This change in alveolar structure is potentially responsible for the increased V/Q heterogeneity demonstrated in older mice as similar mechanisms have been proposed clinically [30]. Regionalization of images further indicated that the increase in standard deviation of the log(V/Q) data is spread throughout the lung, but is also primarily associated with the peripheral volumes. The increase in standard deviation values in the outer region was associated with decreases in the percentages of both ventilation and perfusion associated with this volume. As greater than 30% of ventilation and perfusion counts were still present in this region at 36 weeks of age, and both ventilation and perfusion were decreased, it is unclear from this data what process is causing the increased V/Q mismatch observed. A decrease in perfusion with increasing distance toward the periphery has been noted previously by other methods [31] but the peripheral involvement of V/Q mismatch with age requires further research.
This study employed a per-voxel analysis in order to better depict the intrinsic V/Q heterogeneity in the lungs. Although the voxels used (0.23 mm 3 ) were larger than the alveolar structure, a great number of V/Q measurements were made within the lung label, forty-six thousand on average for a mouse lung. The strength of this approach is that measurements of the whole lung, such as mean and standard deviation, can be broken down into constituent parts. Small areas of mismatched V/Q will not necessarily stand out in global measurements but even these regions, in a voxel-by-voxel analysis, are likely made up of many voxels. The location and extent of mismatching could be quantified and a better understanding gained of how the lung functions on a regional level. We elected to use a segmentation defined from the ungated CT to ensure that all lung-associated voxels within the ungated SPECT images were included in the analysis. Suga et al. [8], in a clinical study, employed a threshold of 10% of maximum in their ventilation image data to define a region of interest for the lungs, but this method fails to address any regions of low ventilation activity and therefore could potentially exclude regions of lung, depending on the disease or model being used. Inclusion of all lung-associated volume is preferable and anatomical imaging, such as CT, is therefore not only desirable but also necessary.
Compared to clinical V/Q SPECT, the method used to deliver Technegas TM had to be modified for this preclinical work; while human patients are capable of inspiring an entire dose of activity in one to three breaths, mice require mechanical ventilation over the course of several minutes with Technegas TM at a relatively high level of radioactivity. Anaesthetic use is a further difference, though we intentionally kept the doses of these agents at a minimum to decrease any possible effect on respiration; respiratory gating confirmed that breathing patterns were similar between animals. Despite the differences between a clinical V/Q and the preclinical technique demonstrated here the results of this study represent a comparable methodology that achieves consistent results.
An important consideration for this methodology is the radiation dose to the subject animal. Work by Travis et al. suggests that the lethal dose for 50% (LD50) for BALB/c mice is above 5Gy and that the dose rate required is well beyond that used in our study. Furthermore, no pathological changes were observed in the lungs or kidneys at a dose-rate of 100 mGy/min [32]. Based on previous whole body dosimetry measurements of mice using the same X-SPECT system, average CT radiation absorption dose can be approximated at 270 mGy [14], delivered over the course of approximately 25 minutes. Further, whole body dose from 99m Tc exposure can be approximated at 51 mGy [33]. Thus, we estimate that the total radiation dose to a mouse from V/Q SPECT/CT is 321 mGy, well below the radiation dose required to have a significant negative impact on biological processes in BALB/c mice. This total dose could be reduced substantially if basic CT scans were acquired for co-registration purposes only instead of the multiple rotations required for high quality respiratory-gated images.
The characterisation of V/Q relationships in the lungs of healthy mice is a foundation upon which lung function can be studied under various physiological conditions and disease models. Respiratory diseases such as COPD would greatly benefit from V/ Q imaging as our understanding of its progressive pathology and its impact on respiratory function could be considerably augmented by this lung function imaging methodology. We believe that the presented work demonstrates that we can reliably measure biomarkers of lung function using a clinically adapted method for performing V/Q SPECT in mice. The methodologies employed for delivery of radiolabelled tracer molecules, spatial alignment of images, and analysis of the resulting data provide a platform for investigations into the functional consequences of lung disease within experimental mouse models.