《1. Introduction》

1. Introduction

Currently, bone-status assessment is mainly based on X-ray densitometry techniques such as quantitative computed tomography (QCT), dual-energy X-ray absorptiometry (DXA), and peripheral quantitative computed tomography (pQCT) [1,2]. The bone mineral density (BMD) provided by DXA is widely used in clinics to evaluate bone quality and predict the risk of bone fracture [3−5]. However, one disadvantage of these X-ray-based methods is the ionizing radiation that restricts their application in specialized populations such as pregnant women and newborn babies. Therefore, developing an alternative technique is of great significance. Quantitative ultrasound (QUS) has drawn great attention as a noninvasive tool for bone-status evaluation [6−11]. Besides the advantages in safety, non-ionizing radiation, portability, and low cost, the ultrasound can reflect information about the bone structure and elastic modulus. QUS for cancellous bone evaluation is mainly classified into two methods: through-transmission and backscattering. Through-transmission measurements use two transducers, one as the transmitter and one as the receiver. These types of measurements are typically performed on skeletal sites that are accessible by two transducers, such as the calcaneus. The broadband ultrasound attenuation (BUA) and speed of sound (SOS) are the two most common through-transmission parameters [12−14]. Numerous clinical QUS instruments have been developed based on the through-transmission method [7,15−17]. Over the last few decades, researchers have focused on the ultrasonic backscatter technique for cancellous bone evaluation [9−11,18]. Unlike through-transmission measurements, the backscatter technique uses a single transducer to both transmit and receive ultrasonic signals and can provide easier access to osteoporosis-sensitive bones, such as the hip and spine. Several studies have demonstrated that the backscatter signal is closely related to bone characteristics including number density, BMD, ultimate strength, size, shape, and elastic properties [19−28]. The backscatter signal is sensitive to the microstructure of cancellous bone, making it feasible to extract trabecular microstructural information from the backscatter signal [29,30]. Recently, some studies performed in vivo backscatter measurements and investigated the feasibility of using the backscatter signal to evaluate bone status in the clinic [11,28,31−34]. In a previous study, we investigated the influence of the selected backscatter signal on measured values of the backscatter parameters and proposed a signal-selection standard for the backscatter measurement of cancellous bone in vitro [35]. The in vivo signal-selection standard for cancellous bone evaluation is waiting to be analyzed.

The bone nutritional status in neonates and early childhood stages is significantly related to skeletal development later in life [32,36]. Lacking the ability to store calcium and phosphorus at birth, neonates have a high risk for metabolic bone disease [37−39]. Therefore, it is of great interest to monitor bone nutritional status in order to prevent metabolic bone diseases in neonates. Currently, neonatal bone evaluation is based on a blood test for specific biochemical markers. However, the biochemical markers (serum calcium, phosphorous, and alkaline phosphatase) are only weakly correlated with BMD [40,41]. Furthermore, blood tests are invasive and not usually preferred for neonates, especially preterm infants. X-ray-based densitometry techniques are not typically used in newborn infants because of the ionizing radiation and relatively low sensitivity [42]. Compared with the routine X-rays and biochemical markers, QUS is safer and can easily be performed on neonates. Previously, researchers have investigated the feasibility of QUS for neonatal bone assessment [32,39].

We developed a novel ultrasonic backscatter bone diagnostic (UBBD) instrument based on the ultrasonic backscatter method. The UBBD provides several advantages, including noninvasiveness, non-ionizing radiation, portability, and simplicity. Recently, the UBBD was used to perform backscatter measurements on adults and newborn babies. Jiang et al. demonstrated that the calcaneus backscatter signal measured with the UBBD was significantly correlated with the BMD in hip and spine ( R = 0.63−0.75, P<0.05) [31]. Zhang et al. measured the backscatter signal in neonates with the UBBD, and found a significant correlation between the backscatter coefficient (BSC) and gestational age and birth length ( R = 0.43−0.47, P<0.001) [32]. These studies demonstrate the utility of the UBBD in the backscatter measurement of cancellous bone in vivo.

In this study, we performed backscatter measurements on 467 neonates (268 males and 199 females) at the left calcaneus using the UBBD and investigated the feasibility of using the backscatter signal to characterize neonatal bone status. The influence of the selected backscatter signal of interest (SOI) on the neonatal bone assessment was analyzed. This study also optimized the SOI selection and attempted to propose an explicit standard for in vivo SOI selection in neonatal cancellous bone assessment.

《2. Material and methods》

2. Material and methods

《2.1. UBBD instrument》

2.1. UBBD instrument

The UBBD was developed based on a Linux-embedded system and has two modules: the ultrasonic signal transmit and receive (USTR) module, and the data-processing and display (DPD) module, as shown in Figure 1(a). The USTR module was based on the field-programmable gate array (FPGA). In the USTR module, a short pulse with a voltage of approximately±50 V was produced to excite the ultrasonic transducer. The backscatter signal received by the same transducer underwent pre-processing and amplification, and was passed through gain-compensation circuits, followed by high-speed analog to digital conversion. Next, the digital backscatter signal was concurrently transmitted to the DPD module. An ARM9 processer was used as the central processing unit in the DPD module. The backscatter signal was analyzed and the backscatter parameters were calculated in the DPD module. Results were displayed on a liquid crystal display (LCD) with a touch screen (9 in, 1 in= 2.54 cm). An external keyboard/mouse and mini-printer were also employed in this instrument.

An interactive interface was developed for the UBBD (Figure 1(b)). Transducers with different central frequencies were interchangeable. The system emitted a short (2 µs) single-cycle sinusoidal pulse for pulse-echo tissue insonification. The echoes were received by the same transducer, amplified and digitized with 12-bit accuracy at data rates as high as 40 000 000 samples per second. The backscatter signal was obtainable within 5 s. To reduce random noise, 128 waveforms were averaged in the time domain and then stored in the system for analysis.

《Fig. 1》

     

Fig.1 (a) The inner functional modules; (b) the surface appearance of the UBBD.

《2.2. Subjects and anthropometry》

2.2. Subjects and anthropometry

A total of 467 neonates (268 males and 199 females) hospitalized in the Neonatal Department in the Children’s Hospital of Fudan University, Shanghai, China, participated in this study. The following exclusion criteria were adopted: congenital malformation, inherited metabolic diseases, abnormalities of the digestive system, or bone diseases. Gestational age was recorded, along with the length, weight, and head circumference at the time of birth. Table 1 lists the anthropometric characteristics of the participated subjects.

《Tab.1》

Tab.1 Anthropometric characteristics of the subjects.

Note: Values are mean ± standard deviation (SD).

This study was approved and conducted by the ethics committee of the Children’s Hospital of Fudan University, and informed consent was obtained from all parents.

《2.3. Backscatter measurements》

2.3. Backscatter measurements

Backscatter measurements were performed with the UBBD. An unfocused transducer with a central frequency of 3.5 MHz (V546, Olympus-Panametrics Inc., Waltham, MA, USA) was used. The diameter of the transducer was 0.25 in and the −6 dB bandwidth was 3.07 MHz (1.62−4.69 MHz). As shown in Figure 2(a), the transducer was placed vertically on the inner part of the heel and the backscatter signal was measured at the left calcaneus for each subject. Ultrasonic gel (Echo Jelly; Aloka Medical Equipment Co., Shanghai, China) was used to couple the transducer to the soft tissue atop the calcaneus of the subject. The sampling frequency of the UBBD was set to 32.0 MHz. A typical backscattered signal acquired from a neonatal calcaneus at 3.5 MHz is shown in Figure 2(b). The signal reflected by a polished steel plate placed in pure water was obtained as the reference signal (shown in Figure 2(c)). After recording calcaneus backscatter signals in neonates, the UBBD system processed and analyzed the data.

《Fig. 2》

     

Fig.2 (a) Backscatter measurement performed at neonatal calcaneus; (b) the measured backscatter signal; (c) the reference signal.

《2.4. Signal analysis》

2.4. Signal analysis

The signals shown in Figure 2(b) consist of the excited signal (duration: 0−2 µs), the signals reflected by the soft tissue and the cortical bone (approximately 6−10 µs), and the backscatter signal from the trabecular bone (10−20 µs). Analysis of the backscattered signal of interest required excluding the signals reflected and backscattered by the soft tissue and cortical bone. As illustrated in Figure 2(b), a rectangular window was used to select the backscatter SOI. The time window was positioned with a time delay after 6 µs. The delayed amount was defined as T1 and the length of the time window was defined as T2. The backscatter parameters were calculated with the selected SOI.

The apparent backscatter coefficient (ABC( f)) is a function of the frequency and was defined as follows [23,24]:

where SSOI( f) and SR ( f) are the amplitude spectrum of the backscattered SOI and the reference signal, respectively.

The apparent integrated backscatter (AIB) is the averaged or integrated value of ABC in the effective frequency band [ fmin, fmax]:

Linear regression was performed on the ABC( f) in the frequency band of interest. The frequency slope of apparent backscatter (FSAB) and zero frequency intercept of apparent backscatter (FIAB) were determined as the slope and the zero frequency intercept of the line fitted to the ABC( f), respectively [20].

The spectral centroid shift (SCS) was defined as follows [28,31]:

where SSOI( f) is the amplitude spectrum of the backscattered SOI; [ fmin, fmax] corresponds to the −6 dB frequency band of the transducer; and f0 is the central frequency of the transducer. The effective frequency band was 1.62−4.69 MHz in this analysis.

Manipulation of the delay ( T1) and duration ( T2) of the time window allowed for the calculation of AIB, FSAB, FIAB, and SCS with different SOIs.

《2.5. Statistical analysis》

2.5. Statistical analysis

The mean and standard deviation (SD) values of the anthropometric parameters of the subjects were calculated (Table 1). The relationship between the four backscatter parameters under different SOI selections and each of the four anthropometric parameters was examined using scatter plots. Linear regression analysis was performed to determine the correlation of the backscatter parameters with the anthropometric parameters. Statistical significance was tested using one-way analysis of variance (ANOVA). A P<0.01 was considered to be statistically significant.

《3. Results》

3. Results

The backscatter SOI was varied through changes in T1 and T2. The correlation between the corresponding backscatter parameters (AIB, FIAB, FSAB, and SCS) and each of the neonatal anthropometric parameters was calculated. The relationship among the neonatal anthropometric parameters was also analyzed. The influence of the SOI selection on the backscatter parameters for neonatal bone evaluation was determined.

《3.1. The relationship among the four anthropometric parameters》

3.1. The relationship among the four anthropometric parameters

Table 2 shows the correlation among the anthropometric parameters. Significantly strong correlations ( R = 0.81−0.92, P<0.001) were observed between the gestational age, birth weight, head circumference, and birth length. The anthropometric parameters were used to characterize the neonatal nutritional status and reflected the bone nutritional status as well. Since the four anthropometric parameters were highly correlated with each other, the gestational age was used as the representative anthropometric for the following analysis.

《3.2. The influence of SOI selection on AIB for neonatal bone evaluation》

3.2. The influence of SOI selection on AIB for neonatal bone evaluation

Figure 3 shows the relationship between AIB and gestational age under two typical SOI selections. The AIB was found to be positively correlated with the gestational age ( R = 0.45, P<0.001) for T1 = 6.50 µs and T2 = 2.34 µs. A negative correlation ( R = −0.56, P<0.001) was observed between AIB and gestational age with T1 = 12.75 µs and T2 = 7.03 µs. The influence of SOI selection on the observed correlation coefficient (red for 1 and blue for −1) for AIB versus gestational age and AIB versus birth weight are shown in Figure 4(a) and Figure 4(b), respectively. In general, positive correlations (| R| up to 0.45, P<0.001) were observed when T1 was short (<8 µs), while negative correlations (| R| up to 0.56, P<0.001) were more common for T1 greater than 10 µs. The phenomena for AIB versus birth weight were quite similar to the correlations between AIB and gestational age, with AIB versus birth weight showing moderate correlation (| R| up to 0.55, P<0.001).

《Tab.2》

Tab.2 Correlations between the anthropometric parameters of the neonates.

     

Note: **: P < 0.001.

《Fig. 3》

Fig.3 The relationship between AIB and gestational age for two typical SOIs. (a) T1 = 6.50 µs, T2 = 2.34 µs; (b) T1 = 12.75 µs, T2 = 7.03 µs.

《Fig. 4》

     

Fig.4 The effect of SOI selection on the correlations between AIB and (a) gestational age, (b) birth weight.

《3.3. The influence of SOI selection on FIAB and FSAB for neonatal bone evaluation》

3.3. The influence of SOI selection on FIAB and FSAB for neonatal bone evaluation

The correlation coefficients for FIAB and gestational age plotted as a function of the T1 and T2 are shown in Figure 5(a) and Figure 5(b), respectively. In Figure 5(a), the correlation was positive when T1 was below 8 µs. The correlation coefficient decreased and eventually turned to negative with increasing T1. Negative correlations ( R = −0.50, P<0.001) were observed when T1 was approximately 13 µs. In Figure 5(b), the correlation coefficient tended to decrease as T2 increased such that the positive correlations weakened while the negative correlations strengthened with increasing T2.

《Fig. 5》

     

Fig.5 The correlation coefficient between FIAB and gestational age plotted as a function of (a) T1 and (b) T2.

《3.4. The influence of SOI selection on FSAB for neonatal bone evaluation》

3.4. The influence of SOI selection on FSAB for neonatal bone evaluation

Figure 6 shows the correlation coefficients between FSAB and gestational age. Generally speaking, significantly positive correlations were observed when T1 was longer than 10 µs. The correlation coefficient was 0.45 when T1 was 12.75 µs and T2 was 5.45 µs. However, when the T1 was short (<10 µs), the correlations were quite weak and insignificant. T2 had little effect on the correlations, and the correlations did not change when T2 was relatively long (>3 µs).

《Fig. 6》

     

Fig.6 The effect of SOI selection on the correl­ations between FSAB and gestational age.

《3.5. The influence of SOI selection on SCS for neonatal bone evaluation》

3.5. The influence of SOI selection on SCS for neonatal bone evaluation

Figure 7(a) shows the correlation coefficient for SCS and gestational age plotted as a function of T1 with respect to some T2 values. The correlations were quite weak and insignificant when T1 was short (<10 µs). With increasing T1, the correlations turned negative and increased. Moderate correlations ( R = 0.41, P<0.001) were observed when T1 was around 13 µs. Figure 7(b) plots the correlation coefficient for SCS and gestational age as a function of T2. T2 has a complex effect on the observed correlations, which means that as the analyzed SOI length varied, the backscatter spectrum also varied.

《Fig. 7》

     

Fig.7 The correlation coefficient between SCS and gestational age plotted as a function of (a) T1 and (b) T2.

《4. Discussion》

4. Discussion

It is well known that bone development is significantly correlated with gestational age, birth weight, and birth length [43−45]. These anthropometric characteristics reflect the skeletal development of neonates and can be used as an index of neonatal bone status. This study analyzed the relationship between the backscatter signal and gestational age, birth weight, birth length, and head circumference. The results showed that the backscatter signal could be significantly correlated with the anthropometric characteristics of the neonates.

Studies have indicated that the microstructure of cancellous bone reflects the bone nutritional status in neonates. Neonates with metabolic bone diseases have significant abnormalities in the microstructure of cancellous bone [45]. Therefore, measuring cancellous bone could help assess bone nutritional status in neonates. The calcaneus was chosen as the site of interest in this study as the bone tissue in the calcaneus is mainly cancellous bone and is easily accessible. The outer cortical shell and the overlying soft tissue in the calcaneus are relatively thin and flat, making it straightforward to place the transducer.

The apparent backscatter was not compensated for the frequency-dependent diffraction and attenuation effects. Therefore, in addition to the backscatter effect, the apparent backscatter parameters (AIB, FSAB, and FIAB) were sensitive to attenuation. The correlation between the apparent backscatter parameters and gestational age resulted from the combined effects of backscatter and attenuation. The correlation (negative or positive) depended primarily on which effect was dominant. When T1 was short (i.e., ultrasonic propagation length was small), the attenuation effect was lower than the backscatter effect and the backscatter effect dominated the observed correlation. For neonates with longer gestational age, the dominance of the backscatter effect was more pronounced than that of the attenuation effect, so the AIB and FIAB values were larger. Therefore, positive correlations were observed for AIB and FIAB with gestational age when T1 was short. In contrast, when T1 was sufficiently long, the attenuation effect (increasing exponentially with the ultrasonic propagation length) dominated the observed correlation. Neonates with longer gestational age had smaller AIB and FIAB values and larger FSAB. This observation is the reason for negative correlations for AIB and FIAB with gestational age and positive correlation for FSAB with gestational age when T1 was long. However, it was found that T1 should not be>16 µs because of heavy attenuation yielding a small SOI with a poor signal-to-noise ratio, which may influence the accuracy of the cancellous bone assessment. The T2 was the SOI length, corresponding to the ultrasonically interrogated trabecular volume. Because cancellous bone is anisotropic and inhomogeneous, the backscatter cross-section varied and the backscatter signal fluctuated as T2 changed. Therefore, the correlations for the backscatter parameters with gestational age fluctuated as T2 varied. Our recent work analyzed the influence of the SOI on the apparent backscatter measurement of cancellous bone in vitro [35]. Significantly positive correlations between AIB and bone volume fraction (bone volume/total volume, BV/TV) were observed for short T1 and negative correlations were observed when T1 was long for the frequency range 0.5−10 MHz. These in vitro results are consistent with the present in vivo findings.

Ultrasonic waves are attenuated as a function of frequency during propagation through biological tissues. In general, the frequency-dependent attenuation has an effect similar to that of a low-pass filter: The signals at higher frequencies are attenuated more than the signals at lower frequencies. As a result, the centroid of the backscattered spectrum is downshifted to a lower value compared with the original spectrum. When a Gaussian pulse propagates in a medium with linear frequency dependence of attenuation, the Gaussian spectrum retains its Gaussian form while the SCS is proportional to the product of the attenuation coefficient, the propagating distance, and the square of the bandwidth [46]. The propagating distance corresponded to T1 in this study. As T1 increased, the SCS was more pronounced. This may explain why significant correlations were observed for SCS and gestational age only when T1 was sufficiently long (approximately 12−16 µs). Note that the backscattered SCS is affected by all of the frequency-dependent factors: attenuation in soft tissue, cortical bone and trabecular bone, transmission coefficients of cortical bone, and multiple scattering effects.

The present study found that the apparent backscatter was significantly correlated with gestational age ( |R| up to 0.45, P<0.001) for short T1 (<8 µs), and negative correlations (| R| up to 0.56, P<0.001) for long T1 (>10 µs). These findings are consistent with our previous study, which found the ABC to be positively correlated with gestational age ( R = 0.19−0.47) for T1<2 µs and T2 = 8 µs [32]. Zhang et al. also found that gender had no influence on bone status at birth [32]. Skeletal development differences between male and female adults are more closely associated with hormone levels, exercise status, and other physiological and environmental factors [43,44]. Therefore, this study did not analyze data based on gender. We did not consider phase interference in the analysis of the backscatter parameters. Phase interference usually occurs in ultrasonic backscatter signals. The constructive or destructive interferences of ultrasonic waves might introduce random errors in the measured backscatter parameters. Phase interference might be one explanation of the weak correlations observed in this study. Due to heavy attenuations in cancellous bone, transducers with lower frequencies (1 MHz or 2.25 MHz) are preferred for cancellous bone assessment for the adult. Considering that the neonatal calcaneus is quite tiny and cancellous bone is not well-developed, transducers with relatively higher frequencies (3.5 MHz and 5 MHz) are usually recommended [32]. Therefore, we used a transducer with a central frequency of 3.5 MHz in this study. We only measured the backscatter signal at the calcaneus for each neonate and did not perform the spatial scanning and average. The highly statistical variance of the cancellous bone structure might have introduced some random errors. If spatial averaging was P performed, the observed correlations between backscatter parameters and anthropometric characteristics of neonates might be stronger.

Other authors have performed ultrasonic transmission measurements in neonates and found that QUS parameters were closely correlated with gestational age and birth weight [39,43−45]. Teitelbaum et al. obtained a positive correlation between SOS in cortical bone and birth weight ( R = 0.3; P<0.001) [39]. Rack et al. found that QUS values from the first week of life were significantly correlated with gestational age ( R = 0.5, P<0.001), and follow-up measurements correlated positively with weight during the week of measurement ( R = 0.4; P = 0.001) [47]. These studies demonstrate that QUS transmission parameters could be used to evaluate bone status in neonates. The present study performed backscatter measurements and demonstrated that the apparent backscatter could be more conveniently used to evaluate neonatal bone status with moderate correlations. This study also analyzed the SOI influence on apparent backscatter measurements of cancellous bone in neonates. Table 3 summarizes the optimization of the observed correlations with explicit SOI selections. The results can provide constructive suggestions for in vivo SOI selection in neonatal cancellous bone assessment.

《Tab.3》

Tab.3 Optimization of the observed correlations with explicit SOI selections.

This study introduced the UBBD, an ultrasonic backscatter instrument for cancellous bone evaluation. Other studies have also used the UBBD for the backscatter measurements. Jiang et al. measured the backscatter signal at the adult calcaneus with UBBD and found that backscatter signal was correlated with hip and spine BMD ( R = 0.63−0.75, P<0.05) [31]. Zhang et al. also performed backscatter measurement in neonatal calcaneus with UBBD and demonstrated that the BSC could be used to assess bone nutritional status in neonates [32]. These studies demonstrate the utility of the UBBD in backscatter measurement of cancellous bone in vivo.

《5. Summary》

5. Summary

The present study performed in vivo backscatter measurements of the calcaneus in neonates with the self-developed UBBD. Through changes in SOI selection, the backscatter parameters were calculated and the influence of SOI selection on the apparent backscatter measurement of neonatal cancellous bone was investigated. We found that the SOI selection significantly affected the correlations between the backscatter and the anthropometric characteristics of neonates. The AIB and FIAB were positively correlated with gestational age ( R2 up to 0.20, P<0.001) for short T1 (<8 µs), while negative correlations ( R2 up to 0.31, P<0.001) were commonly observed for T1>10 µs. Moderate positive correlations ( R2 up to 0.20, P<0.001) were observed for FSAB and SCS with gestational age when T1 was long (>10 µs). Based on the results, this study proposed an explicit standard for in vivo SOI selection and neonatal cancellous bone assessment.

《Acknowledgements》

Acknowledgements

This work was supported by the National Natural Science Foundation of China (11174060, 11327405, and 11504057), the Science and Technology Support Program of Shanghai (13441901900), the PhD Programs Foundation of the Ministry of Education of China (20130071110020), and the China Postdoctoral Science Foundation (2015M571490).

《Compliance with ethics guidelines》

Compliance with ethics guidelines

Chengcheng Liu, Rong Zhang, Ying Li, Feng Xu, Dean Ta, and Weiqi Wang declare that they have no conflict of interest or financial conflicts to disclose.