Literature DB >> 19893940

Doppler sonar analysis of swallowing sounds in normal pediatric individuals.

Cibele Fontoura Cagliari1, Ari Leon Jurkiewicz, Rosane Sampaio Santos, Jair Mendes Marques.   

Abstract

UNLABELLED: Among the methods for assessing swallowing sounds - videofluoroscopy modified barium study, fiberoptic swallowing endoscopy, neck auscultation through a microphone, accelerometer and, more recently, the Doppler sonar - we have chosen the latter. AIM: to analyze swallowing sounds by cervical auscultation using Doppler sonar, in a population between 2 and 15 years without oro-pharyngeal dysfunction. STUDY
DESIGN: cross-sectional historical cohort.
MATERIALS AND METHODS: we investigated 90 individuals in Curitiba (2006/2007). The population was separated by age into 3 groups: from 2 to 5 years, from 5 to 10 years of age and from 10 to 15 years of age. We obtained the average values for frequency, intensity and swallowing duration for saliva, liquid and pasty foods).
RESULTS: objective and measurable data were obtained. Significance related to gender was found in certain age groups and consistencies, under all the studied variables, except swallowing time.
CONCLUSION: neck auscultation using Doppler sonar is a sensitive method to detect swallowing sounds. There was swallowing interference associated with saliva and the other tested food types and with the biological development of the age range being studied. It is an easy to apply method, not expensive and non-invasive.

Entities:  

Mesh:

Year:  2009        PMID: 19893940      PMCID: PMC9442158     

Source DB:  PubMed          Journal:  Braz J Otorhinolaryngol        ISSN: 1808-8686


INTRODUCTION

An instrumental assessment together with the clinical exam of swallowing ultimately contributes to the definition of sequential approaches. Mechanisms which make this practice feasible require new methodologies to be created, tested and explored. Neck auscultation by Doppler sonar is an innovative technique, characterized as a non-invasive method, without radiation exposure, easy to apply and of low cost, enjoying much credibility gain in the clinical assessment of swallowing. All muscle, bone, cartilage and mucosa structures movement associated with food passage produces distinct and successive sounds. Neck auscultation is defined as a method used to hear swallowing sounds with an amplifier during the pharyngeal phase. This method involves the placement of a sensor in the individual's neck in order to hear and/or record the acoustic signals captured by a microphone in order to mainly assess pharyngeal swallowing competence and its interaction with breathing. Food transit sounds can be assessed, just like the mild larynx and pharynx movements which precede and succeed food passage. There is a lot of competitive noise in neck auscultation because the neck is a relatively small region, with constant and significant acoustic activity. Most of the noise is hydraulic in nature because of blood vessels, secretions, CSF and air passing through. Prior studies on neck auscultation included acoustic analyses, physiologic measures, processing and instrumental matters. One of the broadly investigated aspects associated with methodology is the place for pharyngeal swallowing sound detection, where there is a better signal-to-noise ratio. Using an accelerometer, we scanned 24 neck sites. Of these, three were referred to use during neck auscultation, and the area on the lateral trachea, immediately below the cricoid was selected as the best site for swallowing sound detection. The most prominent acoustic characteristic of the swallowing sound corresponding to food movement through the upper esophageal sphincter is found where the movements of the hyoid bone, the larynx and the epiglottis contributes to the swallowing acoustic signal. The respiratory action of the upper airways and the sudden changes in the respective muscles during the swallowing pharyngeal phase are also described as sound components. Changes in spectral characteristics are caused by an increase in food passing velocity, corresponding to the onset of pressurized flow to the stomach. Three swallowing sound components are considered: the first is a weak signal associated with larynx elevation and the food passing through the pharynx; the second, and stronger sound, is associated with the opening of the upper esophageal sphincter; and the third, a weak sound, is associated to the laryngeal downwards movement after swallowing. The swallowing acoustic landmark can be discussed in terms of acoustic signal frequency and duration, and sound wave amplitude. The analysis of such characteristics in the swallowing of liquid and solid consistencies must be carried out in relation to gender and age. The parameters which define swallowing and its disorders in adults can not be always used in children, because of differences in anatomical structure relations. The pediatric larynx is not simply a miniature of the adult larynx, because there are significant differences in size, its location in relation to the spine, in the make up of cartilage and soft tissue and adaptation to the environment, . Sound is a mechanical wave which moves through mediums, produced by a vibration source. Sound-frequencies higher than 25,000Hz, not hearable by humans, are known as ultrasounds. Ultrasonography is based on the phenomena of sound-tissue interaction, that is, based on the transmission of a sound wave through the medium one can see the mechanical properties of tissues. The transducer is the part of the equipment which is responsible for generating, transmitting and capturing ultrasound as it turns electric power into ultrasound power and vice-versa. One of the main peculiarities associated with ultrasound is the very possibility of doing a non-invasive study through the Doppler effect. The Doppler effect is defined as an alteration in the frequency sensation resulting from a situation in which the sound source is mobile, moving at a constant velocity and the receiver is stationary in some point of the trajectory. As the sound source approaches the receiver, the latter receives a larger number of waves per unit of time (higher frequency) and as it moves away it receives a lower number of waves (lower frequency). In continuous Doppler, the signal emitted can be represented by a continuous sinusoidal wave, which amplitude is associated to pressure variations in the propagation mean or to particle shift. The computer is capable of digitalizing the sounds and processing the noise produced by swallowing in wave format visual representations. The feasibility of using the continuous Doppler sonar as a tool to help assess swallowing sounds and to identify acoustic parameters was shown in 50 normal adults swallowing saliva, liquid and pasty food.

OBJECTIVE

Our goal with the present study was to analyze swallowing sounds in individuals from 2 to 15 years of age without oropharyngeal dysphagia obtained through neck auscultation with the Doppler sonar.

MATERIALS AND METHODS

We studied 90 individuals classified according to Chart 1. In order to classify them by age range, months were written as fractions of a year.
Chart 1

Individual distribution in relation to age and gender.

Age RangeFrom 2 to 5 yearsFrom 5 to 10 yearsFrom 10 to 15 years
GenderMFMFMF
n15 (16,66%)15 (16,66%)15 (16,66%)15 (16,66%)15 (16,66%)15 (16,66%)

M= male; F= female; n= number of individuals studied

Individual distribution in relation to age and gender. M= male; F= female; n= number of individuals studied Inclusion criteria were minimum age of 2 years and maximum age of 15 years and normal swallowing in all its phases, proven by a prior evaluation carried out by a speech therapist. Exclusion criteria were respiratory difficulty of any origin and structural or functional alterations which could impact swallowing, information obtained from parents and guardians and the assessment of the researcher. The continuous Doppler device was coupled to a Positivo notebook computer with an Intel Celeron M360 1.4GHz processor, 240MB of RAM memory, video card with integrated graphic accelerator Via UniChrome PRO IGP and 64MB of memory, Combo Drive (DVD player - CD recorder), integrated high definition ALC655 audio compatible with the AC'97, with Windows XP Professional operational system (Fig. 2).
Figure 2

Doppler sonar coupled to a notebook computer

Doppler sonar coupled to a notebook computer The acoustic signals were recorded and analyzed by the VoxMetria software, version 2.8h22, through the voice analysis file (Fig. 3).
Figure 3

Voxmetria software tracing

Voxmetria software tracing

Equipment and Software

We used a portable ultrasound detector (DF-4001 model from Martec) with a flat disc transducer of a single crystal (Figure 1), providing the interface for the Doppler. The ultrasound frequency by the Doppler effect was 2.5MHz, with a 10mW/cm2 output. The sound output power was 1W.
Figure 1

Sonar Doppler sonar

Sonar Doppler sonar

Standardizing Food Consistency

Liquid food consistency was defined as the swallowing of 2.5ml of water given in a disposable cup, and pasty food as the swallowing of 2.5ml of the Danoninho®[1] refrigerated yogurt-like product, given in a spoon.

Procedures

In order to capture swallowing sounds, the individual was instructed to remain seated and have his neck free for access. We made sure that everyone, specially the younger ones, were able to understand verbal commands which were given - under strict compliance with the protocol established. Participant collaboration was uniformly satisfactory. Subjects were instructed to swallow the entire volume offered in one single gulp. Swallowing took place immediately after the investigator's command, about three seconds after placing the transducer on the right neck lateral region. A minimum of three swallowing samples were collected from each one. Saliva swallowing sounds were the first captured, followed by those of liquid and then the pasty substance. Three minutes passed between the swallowing of each different substance. [1] Danoninho® nutritional make up: skim milk, strawberry mix (water, fruit-oligosaccharides, fructose, sugar, strawberry, calcium, zinc, modified starch, vitamins A and D, carmin cochonilha natural dye, citric acid acidulant, carrageen thickener and xantan gum, carboxymethylcellulose and guar gum stabilizers, potassium sorbate preservative and aromatizer), sugar, milk cream, calcium chloride, lacteous ferment, clotting agent and potassium sorbate preservative).

Swallowing sounds capture by the Doppler Sonar

The Doppler Sonar transducer was placed on the right lateral neck, on the tracheal lateral border, just below the cricoid cartilage. The Doppler device volume was set to One, for proper audio signal capture by the VoxMetria software and lower external noise interference. When the noise produced by the common carotid artery pulsation was captured, in some samples, the transducer was repositioned and a new sound capture ensued. The ultrasound energy beam emitted by the transducer was guided to form a 30° to 60° angle. We used contact Carbogel gel to reduce ultrasound scattering and to increase its eco and body transmission.

Acoustic analysis of the Swallowing Sound Signal

Five variables were measured for each saliva swallowing effort and for the liquid and pasty substances: Sound wave initial frequency (FI): frequency at acoustic signal onset, with the frequency window between 60 and 12000Hz; Sound wave peak frequency (FP): frequency at the point of highest acoustic signal shift, with the aforementioned frequency window (Fig. 4);
Figure 4

Sound wave peak and initial frequencies

Sound wave peak and initial frequencies Sound wave initial intensity (II): intensity at acoustic signal onset, with the window between 10 and 100dB; Sound wave peak intensity (IP): intensity at the highest point of the acoustic signal shift, with the aforementioned intensity window (Fig. 5);
Figure 5

Sound wave peak and initial intensities

Sound wave peak and initial intensities Swallowing time (T): time between the beginning and the end of the analyzed acoustic signal, measured by means of an audio signal, in seconds (Fig. 6).
Figure 6

Swallowing time

Swallowing time In order to plot the swallowing sounds from the individuals in terms of frequency, intensity and swallowing duration, we selected the sample considered the best audio and visual presentation. In order to specify saliva swallowing and that of liquid and pasty consistency foods, numbers and colors were standardized to the acronyms of the variables mentioned: number 1 (FI1, FP1, II1, IP1, T1) and the green color for saliva swallowing; number 2 (FI2, FP2, II2, IP2, T2) and the red color for liquid food; number 3 (FI3, FP3, II3, IP3, T3) and the blue color for pasty food. The statistical methodology was made up of descriptive (mean) and inferential (significance test) analysis approach. In order to analyze the significance of acoustic parameters between men and women at each age range for saliva, liquids and pasty foodstuff, the Student t test was used. The ANOVA test was used in order to check the significance within one age range for the three consistencies, and in the same consistency for the three age ranges based on the profile found at each one of the five variables. The significance level was equal to 0.05 where p<0.05 is considered significant. On the tables, significant values are stressed with an asterisk.

RESULT

Mean values and significance according to gender

Mean and significance values according to gender for saliva and liquid and pasty foodstuff for the different age ranges were checked in the initial and peak frequencies (Table 1, Table 2), initial and peak intensities (Table 3, Table 4) and swallowing times (Table 5).
Table 1

Mean values and results from the tests comparing initial frequencies according to gender.

CONSISTENCYAGE RANGES (YEARS)MALESFEMALESP
2 to 5721,91676,160,0394 *
Saliva5 to 10707,18705,280,9235
10 to 15707,42770,120,0003 *
2 to 5710,96704,570,7449
Liquid5 to 10746,85703,620,0223 *
10 to 15707,66751,600,0479 *
2 to 5715,18716,920,8629
Paste5 to 10731,41696,260,0428*
10 to 15700,77746,850,0392*
Table 2

Mean values and results from the peak frequency comparison tests according to gender.

CONSISTENCYAGE RANGES (YEARS)MALESFEMALESP
2 to 51078,631056,580,1472
Saliva5 to 101094,551081,480,0628
10 to 151092,411105,230,0434*
2 to 51086,951093,360,4402
Liquid5 to 101106,901100,010,0554
10 to 151105,001107,370,2124
2 to 51074,361073,410,9302
Paste5 to 101096,681075,550,0277 *
10 to 151091,941103,850,0636
Table 3

Mean values and results from the initial intensity comparison tests according to gender.

CONSISTENCYAGE RANGES (YEARS)MALESFEMALESP
2 to 562,0858,220,0394 *
Saliva5 to 1061,0960,940,9229
10 to 1561,1166,060,0003 *
2 to 561,4560,880,7444
Liquid5 to 1064,2260,810,0223*
10 to 1561,1364,600,0479*
2 to 561,5861,860,8631
Paste5 to 1063,0060,230,0428 *
10 to 1560,5864,220,0392 *
Table 4

Mean values and results from the peak intensity comparison tests according to gender.

CONSISTENCYAGE RANGES (YEARS)MALESFEMALESp
2 to 590,4188,580,1477
Saliva5 to 1091,6790,640,0631
10 to 1591,5092,520,0434 *
2 to 591,0791,580,4412
Liquid5 to 1092,6592,100,0559
10 to1592,5092,620,2149
2 to 590,0890,000,9300
Paste5 to 1091,8490,170,0277*
10 to 1591,4692,420,0635
Table 5

Mean values and results from the swallowing time comparison test according to gender.

CONSISTENCYAGE RANGES (YEARS)MALESFEMALESp
2 to 50,790,810,8223
Saliva5 to 100,790,760,7520
10 to 150,790,750,5396
2 to 50,870,960,3088
Liquid5 to 101,051,000,6274
10 to150,990,970,7250
2 to 50,810,830,7388
Paste5 to 100,850,990,1096
10 to 150,920,810,2605
Mean values and results from the tests comparing initial frequencies according to gender. Mean values and results from the peak frequency comparison tests according to gender. Mean values and results from the initial intensity comparison tests according to gender. Mean values and results from the peak intensity comparison tests according to gender. Mean values and results from the swallowing time comparison test according to gender.

Significance in the same age range for saliva, liquid and pasty consistencies

The type of consistency interference in each one of the age ranges in the five variables studied is depicted by the p values for males (Table 6) and females (Table 7).
Table 6

Significance by variable in the three consistencies in the same age range (males)

AGE RANGEFIFPIIIPT
2 to 5 years0,87410,56400,87430,56340,5807
5 to 10 years0,12230,0380*0,12250,0380*0,0233*
10 to 15 years0,93150,06890,93120,06790,1478

FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time

Table 7

Significance by variable in the three consistencies for the same age range (females)

AGE RANGEFIFPIIIPT
2 to 5 years0,09040,0085*0,09010,0085*0,1575
5 to 10 years0,85390,0139*0,85410,0139*0,0104*
10 to 15 years0,41270,62460,41250,62280,0014*

FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time

Significance by variable in the three consistencies in the same age range (males) FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time Significance by variable in the three consistencies for the same age range (females) FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time Significance of the same consistency in the three age ranges The p values for each one of the variables reflect the age interference in each one of the consistencies presented for males (Table 8) and females (Table 9).
Table 8

Significance for each consistency for the three ages studied (males).

CONSISTENCYFIFPIIIPT
Saliva0,72830,21510,72830,21510,9970
Liquid0,11710,0048*0,11730,0048*0,2140
Pasty0,33120,0276*0,33130,0276*0,4439

FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time

Tabela 9

Significancia para cada consistencia nas tres subfaixas etàrias (sexo feminino).

CONSISTENCYFIFPIIIPT
Saliva0,0000*0,0000*0,0000*0,0000*0,6318
Liquid0,0198*0,0033*0,0198*0,0033*0,6447
Pasty0,0321*0,0031*0,0321*0,0031*0,0827

FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time

Significance for each consistency for the three ages studied (males). FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time Significancia para cada consistencia nas tres subfaixas etàrias (sexo feminino). FI=initial frequency; FP=peak frequency; II=initial intensity; IP=peak intensity; T=swallowing time

DISCUSSION

In a literature review we did not find any paper mentioning the use of this methodology in the age range between two and 15 years. There was only one paper involving Doppler sonar as a tool to capture and analyze swallowing sounds in adults, using different parameters from ours. By accepting that the larynx and oropharyngeal tract structures vibrate and move as they work as a system of valves and pumps causing the swallowing sounds, the many characteristics of these sounds would initially depend on the food consistency. However, they would also depend on gender and age as determining factors of the “resonance box” volume and energy applied to those structures based on modulators of those sounds. The age range division, the use of different consistencies and the pairing among the individuals by gender aimed at trying to achieve a comparative analysis among the acoustic signal frequency and intensity variables and the sound wave time. If we assume that during swallowing at a higher frequency and intensity a shorter time corresponds to a more efficient performance of the swallowing process, one could expect to have linear series and a scaled performance parallel to age progression and neurological maturity. Nonetheless, this did not happen to the majority of the modification series obtained. Further studies, involving a larger number of participants may stress differences we did not find.

Gender-related significance

As we consider initial frequency and intensity mean values associated with swallowing saliva, we notice a gender-related significance in the 2-5 and 10-15 year age ranges. In liquids and paste food, the significance is similar in the 5-10 (higher mean values among males) and 10-15 years (higher mean values among females). The swallowing initial intensity and frequency mean values considering the three different food consistencies, in the age range of 10-15 years are always higher among females. Considering gender significance, this age range was the only one present in this variable among the three different food consistencies, which suggest better performance among girls (Table 1, Table 3). As far as gender is concerned, the frequency and peak intensity variables were significant in saliva consistency from 10-15 years, and pasty consistency among those between 5 and 10 years of age. The frequency and peak intensity mean values in all the three food consistencies studied are the lowest in the 2-5 years age range in comparison to the other ages, for both genders (Table 2, Table 4). When we consider that swallowing sounds involve laryngeal elevation we can notice that the lowest values found in these variables between 2-5 years in both genders would be caused by the lowest laryngeal dimension in that age range. Frequency and peak intensity mean values are always higher in the 5-10 age range in males and in the 10-15 age range in females, in all the consistencies studied. The high mean values in the 10-15 age range, in the peak and initial intensity and frequencies seem to point to a clear superiority in swallowing performance among females, with characteristics of an accelerated pace of development. Considering liquid, the peak intensity and frequency mean values were always higher when compared by age range in relation to the other consistencies, both in males and females (Table 2, Table 4). The high values seen among these two variables for liquid swallowing suggest the need for higher upper airway protection with a higher laryngeal elevation in relation to swallowing saliva and pasty food. There was no gender-related significance in the swallowing time variable in any of the age ranges or consistency type (Table 5). In relation to saliva swallowing mean times, we noticed they were equal in all age ranges for males, showing uniformity as far as age is concerned. Nonetheless, among females there is a lowering in age ranges from younger to older ages, expressing an enhancement as time passes. When we compared the swallowing time mean values among the three consistencies, considering the same sex and age range, saliva consistency swallowing requires a shorter time to be effective. When we analyze the mean values of this variable, we notice the inverse happening in relation to liquid consistency, where the highest values were found.

Significance in the same age range in saliva, liquid and pasty consistencies

Among males, we noticed significance in the frequency and peak intensity variables and swallowing time only among those between 5 and 10 years (Table 6). Among females, we had significant values in the peak intensity and frequency in the age range from 2-5 years and 5-10 years of age. In relation to swallowing time, age ranges 5-10 and 10-15 years showed significance (Table 7).

Significance for each food consistency in the three age ranges

The mean value variance within the same consistency, for the three age ranges, among males (Table 8) was seen for liquid and paste, only among frequency and peak intensity variables. Contrary to that, among females (Table 9), significance was marked for being present in all consistencies and in all the variables, except swallowing time.

FINAL REMARKS

The data obtained from normal children may serve as comparison basis for our testing of pediatric patients with swallowing abnormalities, in order to establish disease patterns. On the other hand, the alterations seen in the tests of these children with anatomical and physiological alterations can serve for a better interpretation of the results hereby mentioned. Swallowing sound interpretation and its correlation with physiological events which happen during this process do not seem yet clear. Being able to make this association will provide neck auscultation with greater credibility as instrumental diagnosis. The Doppler neck spectrography together with other swallowing assessment data such as videofluoroscopy and nasal-laryngeal fibroscopy can clarify the physiological causes of swallowing sounds.

CONCLUSION

The Doppler sonar, gathering visual and acoustic data to detect swallowing sound characteristics, provided objective and measurable data. The variance analysis of the mean values found in each variable - frequency, swallowing time and intensity - proved significant in relation to gender, age range and food consistencies, with indications of a direct interference relation of swallowing and the biological development of the age range studied.
  12 in total

1.  Sounds of swallowing following total laryngectomy.

Authors:  S L Hamlet; R L Patterson; S M Fleming; L A Jones
Journal:  Dysphagia       Date:  1992       Impact factor: 3.438

2.  Sound component duration of healthy human pharyngoesophageal swallowing: a gender comparison study.

Authors:  Sylvain Morinière; Patrice Beutter; Michèle Boiron
Journal:  Dysphagia       Date:  2006-07       Impact factor: 3.438

3.  Interpreting the sounds of swallowing: fluid flow through the cricopharyngeus.

Authors:  S L Hamlet; R J Nelson; R L Patterson
Journal:  Ann Otol Rhinol Laryngol       Date:  1990-09       Impact factor: 1.547

Review 4.  The physiologic cause of swallowing sounds: answers from heart sounds and vocal tract acoustics.

Authors:  J A Cichero; B E Murdoch
Journal:  Dysphagia       Date:  1998       Impact factor: 3.438

5.  Detection of swallowing sounds: methodology revisited.

Authors:  Julie A Y Cichero; Bruce E Murdoch
Journal:  Dysphagia       Date:  2002       Impact factor: 3.438

6.  Cervical auscultation of suckle feeding in newborn infants.

Authors:  F L Vice; J M Heinz; G Giuriati; M Hood; J F Bosma
Journal:  Dev Med Child Neurol       Date:  1990-09       Impact factor: 5.449

7.  Methodology for detecting swallowing sounds.

Authors:  K Takahashi; M E Groher; K Michi
Journal:  Dysphagia       Date:  1994       Impact factor: 3.438

8.  Stethoscope acoustics and cervical auscultation of swallowing.

Authors:  S Hamlet; D G Penney; J Formolo
Journal:  Dysphagia       Date:  1994       Impact factor: 3.438

9.  Collagen in the developing larynx. Preliminary study.

Authors:  S R Cohen; D T Cheung; M E Nimni; V Mahnovski; G Lian; N Perelman; A P Carranza
Journal:  Ann Otol Rhinol Laryngol       Date:  1992-04       Impact factor: 1.547

10.  Reliability and validity of cervical auscultation.

Authors:  Christiane Borr; Martina Hielscher-Fastabend; Andy Lücking
Journal:  Dysphagia       Date:  2007-04-25       Impact factor: 3.438

View more
  3 in total

1.  Characterization of swallowing sounds with the use of sonar Doppler in full-term and preterm newborns.

Authors:  Hellen Nataly Correia Lagos; Rosane Sampaio Santos; Edna Marcia da Silva Abdulmassih; Liliane Friedrich Gallinea; Mariangela Langone
Journal:  Int Arch Otorhinolaryngol       Date:  2013-10

2.  Development of an Ultrasonic Doppler Sensor-Based Swallowing Monitoring and Assessment System.

Authors:  Younggeun Choi; Minjae Kim; Baekhee Lee; Xiaopeng Yang; Jinwon Kim; Dohoon Kwon; Sang-Eok Lee; HyangHee Kim; Seok In Nam; Saewon Hong; Giltae Yang; Duk L Na; Heecheon You
Journal:  Sensors (Basel)       Date:  2020-08-13       Impact factor: 3.576

3.  Acoustic analysis of oropharyngeal swallowing using Sonar Doppler.

Authors:  Franciele Savaris Soria; Roberta Gonçalves da Silva; Ana Maria Furkim
Journal:  Braz J Otorhinolaryngol       Date:  2015-12-17
  3 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.