Literature DB >> 30697535

Executive Function in Relation to White Matter in Preterm and Full Term Children.

Irene M Loe1, Jenna N Adams2, Heidi M Feldman1.   

Abstract

Background: Executive function (EF) refers to cognitive abilities used to guide goal-directed behavior. Diffusion Tensor Imaging (DTI) provides quantitative characterization of white matter tracts in the brain. Children with preterm birth often have EF impairments and white matter injury. Aim: To examine the degree of association between EF scores and white matter fractional anisotropy (FA) as measured by DTI in children born preterm and term Study design: Cross-sectional study Subjects: Participants, 9-16 years of age, born preterm (n = 25; mean gestational age 28.6 weeks; mean birth weight 1,191 grams), and full term (n = 20) Outcome measures: White matter FA analyzed with Tract-Based Spatial Statistics, a technique that generates a skeleton representing the core of white matter tracts throughout the brain. Behavioral scores from EF tasks examining working memory, spatial memory capacity, and multiple skills from the Stockings of Cambridge.
Results: The groups performed comparably on all tasks. In both groups, unfavorable working memory strategy scores were associated with lower FA. Other measures of EF were not associated with whole skeleton FA in either group in either direction. Conclusions: Strategy score on a spatial working memory task was associated with FA in preterm and full term children, suggesting common underlying neurobiology in both groups. Associations were found in frontal-parietal connections and other major tracts. Lack of associations between other EF tasks and FA may be due to variation in how children accomplish these EF tasks. Future research is required to fully understand the neurobiology of EF in children born preterm.

Entities:  

Keywords:  diffusion tensor imaging; executive function; premature birth; preterm birth; white matter; working memory

Year:  2019        PMID: 30697535      PMCID: PMC6341022          DOI: 10.3389/fped.2018.00418

Source DB:  PubMed          Journal:  Front Pediatr        ISSN: 2296-2360            Impact factor:   3.418


Introduction

Executive Function (EF) refers to interrelated cognitive abilities (i.e., working memory, response inhibition, cognitive flexibility, and organization/planning) used to guide goal-directed behavior (1). EF skills are associated with academic achievement and social competence (2–4). Functional neuroimaging studies link EF skills to prefrontal, parietal, and temporal cortices (5–7). EF performance has also been associated with white matter connections among brain regions (8, 9). EF impairments are common in preterm children (10, 11). They are attributed, in part, to white matter injury (8, 9), which is common after preterm birth (12–15). Associations have been reported between EF skills in childhood and newborn white matter abnormalities on conventional magnetic resonance imaging (MRI) (16–19), which provides a qualitative impression of white matter properties. Diffusion tensor imaging (DTI) generates quantitative measures to characterize white matter (20). Fractional anisotropy (FA) measures the degree to which diffusion is anisotropic, reflecting directional coherence of axons within a white matter tract (21, 22). Reduced FA in preterm children compared to full term children has been found in multiple brain regions in the absence of cystic lesions (12, 14, 15) and has been used to infer subtle injury, though not all studies find differences in all tracts (23–25). Very few DTI studies of preterm children have focused on EF skills and white matter properties. One DTI study of very low birth weight adolescents found that low FA in the left cingulum and bilateral inferior fronto-occipital fasciculi (IFOF) was significantly correlated with low EF (9). Another DTI study found that diffusivity measures other than FA at term-equivalent age in the inferior occipital and cerebellar regions were associated with impaired EF at age 7 years (26). In contrast, a DTI study of 8-year-old very preterm children found no association between DTI measures and parent-rated EF (27). Techniques such as resting state functional MRI and structural connectivity may also be combined to better understand the relationship between white matter maturation and the development of functional connectivity among brain networks (28, 29). A study of typical development between ages 2 and 18 years found positive correlations between structural and functional connectivity and that the relationship strengthened with age (28). A study utilizing diffusion MRI with whole-brain mapping of connectivity in 11- to 31-month-old preterm infants found altered connection strength associated with age and preterm birth (30). Older children had stronger connections in tracts involving frontal lobe structures, whereas increased prematurity at birth was related to widespread reductions in connection strength in tracts involving all cortical lobes and several subcortical structures (30). This exploratory study used DTI in preterm children and adolescents to examine associations between EF skills and FA using Tract-Based Spatial Statistics (TBSS), a whole-brain, mean white matter FA skeleton approach. We focused on FA rather than other DTI metrics, such as axial diffusivity (AD) or radial diffusivity (RD), since both AD and RD are reflected in FA (i.e., FA = AD/RD2). In this sample, we previously reported group differences in FA (25), and significant associations were found between white matter FA in multiple regions and macroscopic injury (25), IQ (25), language (31), and parent-reported behavior symptoms (32). Therefore, we hypothesized that we would find significant associations between measures of EF skills and white matter FA. We suspected significant correlations in tracts that connect frontal and parietal regions, such as the superior longitudinal fasciculus and the arcuate fasciculus, given that other studies have found associations in these tracts or activations in frontal and parietal regions (33–36).

Methods

Participants

Participants were 9–16 year olds, born between 1992 and 1999, and participating in a larger multi-site study (37). The study was carried out in accordance with the recommendations of the Belmont Report, released by the National Commission for the Protection of Human Subjects in Biomedical and Behavioral Research. The protocol (#6985) was approved by the Stanford University IRB. All parents or guardians gave written informed consent and child and adolescent participants gave assent in accordance with the Declaration of Helsinki. Preterm subjects (n = 25) were born at ≤34 weeks gestation with birth weight <2,500 grams. Controls were born ≥37 weeks (n = 20). Participants included in the current study are those who had complete EF data and concurrent DTI data. Due to protocol changes in the larger study, not all subjects completed EF testing, resulting in a smaller number of controls. Exclusions included seizures; hydrocephalus; receptive vocabulary score <70; hearing loss; and non-English speaker. Controls were also excluded for language, learning, or psychiatric disorders. Three preterm subjects were excluded due to extremely enlarged ventricles (38). Preterm subjects were recruited by letters to families previously evaluated at the neonatal follow-up clinic. Controls were recruited through ads and group-matched to preterm children for age, sex, and race. There were no differences in maternal education, IQ scores, or handedness. See Table 1. Medical complications are in Supplemental Material.
Table 1

Participant demographics and behavior scores (N = 45).

Participant CharacteristicsPreterm (n = 25)Full Term (n = 20)
Number (n)%Number (n)%X2p
Boys124810500.8940.566
Maternal Education <High School3125250.4350.229
Non-white83210501.50.221
Left-handed or ambidextrous4164200.0810.775
Special education728006.630.010*
Mean (SD)RangeMean (SD)Rangetp
Age (years)12.7 (2.1)9.7–16.812.9 (2.1)9.4–16.90.2030.654
IQ109 (15.8)72–136113 (16.1)86–1420.7590.388
Gestational Age (weeks)28.6 (2.5)26–3439.3 (1.1)37–40300<0.001*
Birth Weight (grams)1191 (466)630–24103171 (403)2438–3771226<0.001*
EXECUTIVE FUNCTION OUTCOMES
Spatial Working Memory
Strategy Score31.6 (5.7)20–4033.8 (5.1)20–401.370.178
Forgetting Errors24.3 (17.1)0–5023.9 (13.9)4–49−0.080.936
Spatial Span
Span Length6.4 (1.3)4–96.7 (1.6)3–90.620.539
Stockings of Cambridge
Problems Solved7.9 (2.4)3–129.0 (1.9)6–121.710.094

indicates significance at p <0.05.

Participant demographics and behavior scores (N = 45). indicates significance at p <0.05.

Executive Function Measures

We utilized the Cambridge Neuropsychological Test Automated Battery (Cambridge Cognition, Ltd, Cambridge, UK), a computerized EF battery used in studies of preterm children (39). Spatial Working Memory (SWM), Spatial Span (SS), and Stockings of Cambridge (SOC) subtests captured different EF constructs. SWM measures the ability to retain and manipulate spatial information. Two outcomes include Strategy score [a low value indicates a systematic approach (40)] and Forgetting errors. SS measures spatial capacity and is a visual-spatial analog of the digit span task. Span length is the longest sequence successfully recalled. SOC measures response inhibition, planning, and organization skills. A high value on Problems solved indicates good planning/organization skills. Additional details are in Supplemental Material.

Statistical Analysis

Chi-square, independent t-tests, and Mann-Whitney U tests were used to evaluate between-group differences on demographic variables and outcome measures.

Imaging Protocols

Details for MRI and DTI protocols/processing were previously reported (32) and are in Supplemental Material. MRI data were acquired on a 3T Signa Excite (GE Medical Systems, Milwaukee, WI). For DTI, a diffusion-weighted, single-shot, spin-echo, echo-planar imaging sequence was used to acquire 60 slices, 2 mm thick, in 30 different diffusion directions (b = 900).

DTI Image Analysis and Post-processing

We used the voxel-wise “Tract Based Spatial Statistics” (TBSS) method from the Oxford Center for Functional MRI of the Brain (FMRIB) Diffusion Toolbox (41). TBSS identifies a core white matter “skeleton” that is anatomically equivalent across subjects. Advantages of TBSS include avoiding spatial smoothing (averaging voxels) and minimizing partial volume effects that can occur when more than one tract goes through a voxel, thereby leading to decreases in anatomic specificity (41). We analyzed FA, a ratio from 0 to 1. A tensor model was fit to each voxel (42), and FA images were calculated. FA images were registered to the most representative subject in the study (target) using the FMRIB Software Library's (FSL) non-linear registration tool (43). Each subject was aligned to the target; aligned FA images were averaged to create a mean FA map. An FA threshold ≥0.2 was designated to include major white matter pathways and to exclude peripheral tracts. The mean FA image was generated and thinned to create the mean FA-skeleton, a representation of the core of all major white matter tracts common to the group. For each subject, TBSS projected the highest FA on to the FA-skeleton, corresponding to the local tract core.

Statistical Analysis of DTI Data

Analysis of the FA values at the voxel level was conducted using a permutation-based tool for non-parametric statistical thresholding (“Randomize” function in FSL) to assess group related differences (44). This method generates non-parametric, two-sample, unpaired t-tests of FA in subjects compared with controls. We assessed the degree of association between EF scores and FA of the mean skeleton within the preterm group and the control group separately in both the positive and negative directions as well as group by task interactions, covarying by age. We set significance at p < 0.05 after correction for multiple comparisons using Threshold-Free Cluster Enhancement, a method that avoids using an arbitrary threshold for initial cluster-formation or identification of signal (45). Significant clusters of voxels on the FA-skeleton were inspected with reference to John Hopkins University DTI-based white matter atlases (46) to determine tract assignments.

Results

No significant behavioral differences were found between groups on raw scores on any EF measure (Table 1). On voxel-wise analysis, for both preterm and full term groups, SWM strategy scores were significantly negatively associated with the mean FA-skeleton; as strategy scores increased (indicating poor performance), FA decreased. There were no significant positive correlations. There were no significant associations in either direction for the other measures in either group. We also examined interaction effects across both groups; there were no significant group by task interactions. In both groups, tracts associated with strategy score included the superior longitudinal fasciculus (SLF), cingulum, inferior fronto-occipital fasciculus (IFOF), inferior longitudinal fasciculus (ILF), forceps major, and corticospinal tracts. The preterm group also had associations with the corpus callosum and forceps minor, whereas full term children had additional associations with the anterior thalamic radiations. Figure 1 shows the multiple regions of significant association on axial views of the FA-skeleton for the full term group (upper row) and the preterm group (lower row).
Figure 1

White matter regions on the fractional anisotropy (FA) skeleton (green) of statistically significant associations between FA and spatial working memory strategy score in full term (upper row, red) and preterm (bottom row, red) groups (p < 0.05, non-parametric permutation test, corrected for multiple comparisons) overlaid on axial slices of the mean FA image (grayscale). Tracts in which significant associations were found include: Both groups: a. cingulum; b. forceps major; c. inferior fronto-occipital fasciculus; d. inferior longitudinal fasciculus; e. superior longitudinal fasciculus; f. corticospinal tract. Preterm group only: g. corpus callosum/forceps minor. Full term group only: h. anterior thalamic radiation. L, left; R, right; A, anterior; P, posterior.

White matter regions on the fractional anisotropy (FA) skeleton (green) of statistically significant associations between FA and spatial working memory strategy score in full term (upper row, red) and preterm (bottom row, red) groups (p < 0.05, non-parametric permutation test, corrected for multiple comparisons) overlaid on axial slices of the mean FA image (grayscale). Tracts in which significant associations were found include: Both groups: a. cingulum; b. forceps major; c. inferior fronto-occipital fasciculus; d. inferior longitudinal fasciculus; e. superior longitudinal fasciculus; f. corticospinal tract. Preterm group only: g. corpus callosum/forceps minor. Full term group only: h. anterior thalamic radiation. L, left; R, right; A, anterior; P, posterior.

Discussion

In summary, we found that FA in multiple tracts was associated with a single measure of EF—strategy score. Lower FA was associated with poorer strategy in both preterm and full term groups. Similar tracts were identified for both groups, although unique tracts were also associated with each group. Our results are similar to another report of associations between lower FA and EF measures of problem solving and cognitive flexibility in a preterm sample within the cingulum and bilateral IFOF (9). The cingulum bundles course from the medial temporal lobe to the frontal lobes, and frontal regions have long been implicated in EF studies of typical adults (47, 48) and clinical populations using DTI (49–51). The IFOF courses from the occipital lobes through the external capsules to the inferior frontal cortex. A study of very preterm children at age 7 years found that diffusivity measures other than FA at term-equivalent age in similar inferior occipital regions were associated with poor planning ability (26). The associations with these tracts may reflect the importance of visual input to frontal regions during processing of the task. The association between FA and strategy score in both preterm and full term groups suggests a common neurobiological basis for memory strategy. In both groups, FA in the SLF was also associated with strategy score. These findings are similar to those in a study of typical children age 7–13 years that used the same CANTAB Spatial Working Memory task (34). Another study of typical children, ages 5–17 years, found that higher FA in the SLF was correlated with better set shifting/cognitive flexibility and inhibition (52). Associations of FA and strategy score were also found within the ILF, forceps major, and corticospinal tracts for both groups. The ILF courses from the occipital to temporal lobes. Adult DTI studies showed associations between the ILF and a memory composite score in healthy adults (48) and working memory in adults with diffuse axonal injury (49). The forceps major extends from the corpus callosum to connect homologous regions of the occipital lobe in each hemisphere. A tractography study of healthy aging in adults found associations between the splenium of the corpus callosum and memory and EF (53). The corticospinal tracts are associated with motor domains; associations may reflect motor planning and actions required to perform the task. In the preterm group, the corpus callosum and forceps minor were also significantly associated with SWM strategy score. The corpus callosum provides bi-hemispheric connections, and the forceps minor extend from the corpus callosum to connect homologous regions of the anterior frontal lobe in each hemisphere. A study of healthy 8–18 year olds found correlations between FA in frontal-parietal connections and neural activation on functional MRI in the superior frontal sulcus and inferior parietal lobe (54). In the full term group, the anterior thalamic radiations were identified. The ATR connects subcortical regions to the frontal lobe. A DTI study of patients with first-episode psychosis found associations of decreased FA in the ATR, among other regions, and impairments on a measure of cognitive flexibility (55). Subcortical-frontal connections are also important for EF in other DTI studies of children with Attention Deficit Hyperactivity Disorder (ADHD) (56, 57). We did not find associations between FA and SS or SOC in either the preterm or full term groups. These findings are in sharp contrast with other studies using the same group of participants and methods (25, 31, 32). These studies identified many positive associations between FA and other domains of function in PT children, including intelligence, language (i.e., linguistic processing speed), reading (i.e., syntactic comprehension and decoding), and behavior problems (i.e., inattention and internalizing symptoms) (25, 31, 32). The tracts associated with SWM strategy show some overlap with those associated with IQ, language, reading, and behavior symptoms, although the areas of correlation with EF are less widespread than those seen with language, reading and behavior outcomes. The reasons for the lack of findings with EF are unclear. Children may use a variety of different strategies for these tasks that differ across groups and/or represent higher individual variation within groups. These strategies may also vary with age and may have further increased the heterogeneity in this sample with wide age range. If these possibilities are true, the associations of EF skills and white matter of any specific tracts may be limited. Children may approach working memory strategy more consistently than they approach other EF skills. The lack of findings with other EF skills may also reflect methodological details. SWM strategy has a greater range than SS total span length or problems solved on the SOC. It is easier to detect correlations when the range of scores is high. Another possibility relates to limitations of TBSS, which focuses on the centers of white matter tracts and may not detect peripheral damage.

Limitations

The study utilized a small, heterogeneous convenience sample of preterm children who may not be representative of all preterm children. The EF tasks did not differentiate the groups behaviorally. The wide age range may have limited the ability to find brain-behavior relationships and group differences, as the groups may have different developmental trajectories for EF skills in middle childhood and early adolescence.

Conclusion

Use of a systematic strategy on a spatial working memory task was correlated with white matter FA in preterm and full term children, suggesting common underlying neurobiology in both groups. Better performance was associated with higher FA. Associations were found in frontal-parietal connections and other major tracts. Our study adds to the relatively small body of literature on EF and white matter in preterm children. Replication with a larger, more homogeneous sample may provide insight into the neurobiology underlying EF in children in general. Unexpectedly, no other associations between EF scores and FA were found in either group. Future work should focus on smaller age groups and adjust tasks and DTI methodology, including advanced imaging and analysis methods such as combined functional MRI and structural connectivity, to better understand the neurobiology of EF.

Author Contributions

IL conceptualized the study, analyzed and interpreted data, drafted the initial manuscript and critically reviewed, and revised the manuscript. JA contributed to conceptualization and interpretation of data, literature review, and critically reviewed and revised the manuscript. HF conceptualized the study, analyzed and interpreted data, and critically reviewed and revised the manuscript. All authors are responsible for the reported work and have approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
  55 in total

1.  Neurobehavioral evidence for working-memory deficits in school-aged children with histories of prematurity.

Authors:  M Luciana; L Lindeke; M Georgieff; M Mills; C A Nelson
Journal:  Dev Med Child Neurol       Date:  1999-08       Impact factor: 5.449

2.  The unity and diversity of executive functions and their contributions to complex "Frontal Lobe" tasks: a latent variable analysis.

Authors:  A Miyake; N P Friedman; M J Emerson; A H Witzki; A Howerter; T D Wager
Journal:  Cogn Psychol       Date:  2000-08       Impact factor: 3.468

3.  Increased brain activity in frontal and parietal cortex underlies the development of visuospatial working memory capacity during childhood.

Authors:  Torkel Klingberg; Hans Forssberg; Helena Westerberg
Journal:  J Cogn Neurosci       Date:  2002-01-01       Impact factor: 3.225

4.  Nonparametric permutation tests for functional neuroimaging: a primer with examples.

Authors:  Thomas E Nichols; Andrew P Holmes
Journal:  Hum Brain Mapp       Date:  2002-01       Impact factor: 5.038

5.  Combined analysis of DTI and fMRI data reveals a joint maturation of white and grey matter in a fronto-parietal network.

Authors:  Pernille J Olesen; Zoltan Nagy; Helena Westerberg; Torkel Klingberg
Journal:  Brain Res Cogn Brain Res       Date:  2003-12

6.  Object working memory deficits predicted by early brain injury and development in the preterm infant.

Authors:  Lianne J Woodward; Jamie O Edgin; Deanne Thompson; Terrie E Inder
Journal:  Brain       Date:  2005-09-08       Impact factor: 13.501

7.  Changes in white matter diffusion anisotropy in adolescents born prematurely.

Authors:  Torgil R Vangberg; Jon Skranes; Anders M Dale; Marit Martinussen; Ann-Mari Brubakk; Olav Haraldseth
Journal:  Neuroimage       Date:  2006-07-14       Impact factor: 6.556

Review 8.  Development of a superior frontal-intraparietal network for visuo-spatial working memory.

Authors:  Torkel Klingberg
Journal:  Neuropsychologia       Date:  2006-01-06       Impact factor: 3.139

9.  Brain areas involved in spatial working memory.

Authors:  Marieke van Asselen; Roy P C Kessels; Sebastiaan F W Neggers; L Jaap Kappelle; Catharina J M Frijns; Albert Postma
Journal:  Neuropsychologia       Date:  2005-11-21       Impact factor: 3.139

10.  Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data.

Authors:  Stephen M Smith; Mark Jenkinson; Heidi Johansen-Berg; Daniel Rueckert; Thomas E Nichols; Clare E Mackay; Kate E Watkins; Olga Ciccarelli; M Zaheer Cader; Paul M Matthews; Timothy E J Behrens
Journal:  Neuroimage       Date:  2006-04-19       Impact factor: 6.556

View more
  10 in total

1.  Neuroanatomical correlates of impulsive traits in children aged 9 to 10.

Authors:  Max M Owens; Courtland S Hyatt; Joshua C Gray; Joshua D Miller; Donald R Lynam; Sage Hahn; Nicholas Allgaier; Alexandra Potter; Hugh Garavan
Journal:  J Abnorm Psychol       Date:  2020-09-07

2.  White matter fiber-specific degeneration in older adults with metabolic syndrome.

Authors:  Christina Andica; Koji Kamagata; Wataru Uchida; Kaito Takabayashi; Keigo Shimoji; Hideyoshi Kaga; Yuki Someya; Yoshifumi Tamura; Ryuzo Kawamori; Hirotaka Watada; Masaaki Hori; Shigeki Aoki
Journal:  Mol Metab       Date:  2022-06-09       Impact factor: 8.568

3.  Associations Between White Matter Microstructures and Cognitive Functioning in 8-Year-Old Children: A Track-Weighted Imaging Study.

Authors:  Rajikha Raja; Xiaoxu Na; Alexandra Moore; Raymond Otoo; Charles M Glasier; Thomas M Badger; Xiawei Ou
Journal:  J Child Neurol       Date:  2022-03-07       Impact factor: 2.363

4.  Direct and Indirect Associations of Widespread Individual Differences in Brain White Matter Microstructure With Executive Functioning and General and Specific Dimensions of Psychopathology in Children.

Authors:  Carlos Cardenas-Iniguez; Tyler M Moore; Antonia N Kaczkurkin; Francisco A C Meyer; Theodore D Satterthwaite; Damien A Fair; Tonya White; Elisabet Blok; Brooks Applegate; Lauren M Thompson; Monica D Rosenberg; Donald Hedeker; Marc G Berman; Benjamin B Lahey
Journal:  Biol Psychiatry Cogn Neurosci Neuroimaging       Date:  2020-11-25

5.  White matter microstructure and cognitive outcomes in relation to neonatal inflammation in 6-year-old children born preterm.

Authors:  Sarah E Dubner; Cory K Dodson; Virginia A Marchman; Michal Ben-Shachar; Heidi M Feldman; Katherine E Travis
Journal:  Neuroimage Clin       Date:  2019-04-22       Impact factor: 4.881

6.  Altered brain metabolism contributes to executive function deficits in school-aged children born very preterm.

Authors:  Barbara Schnider; Ruth Tuura; Vera Disselhoff; Bea Latal; Flavia Maria Wehrle; Cornelia Franziska Hagmann
Journal:  Pediatr Res       Date:  2020-06-26       Impact factor: 3.756

Review 7.  Neurobehavioral Phenotype and Dysexecutive Syndrome of Preterm Children: Comorbidity or Trigger? An Update.

Authors:  Catherine Gire; Aurélie Garbi; Meriem Zahed; Any Beltran Anzola; Barthélémy Tosello; Valérie Datin-Dorrière
Journal:  Children (Basel)       Date:  2022-02-11

8.  Early Stopping in Experimentation With Real-Time Functional Magnetic Resonance Imaging Using a Modified Sequential Probability Ratio Test.

Authors:  Sarah J A Carr; Weicong Chen; Jeremy Fondran; Harry Friel; Javier Sanchez-Gonzalez; Jing Zhang; Curtis Tatsuoka
Journal:  Front Neurosci       Date:  2021-11-04       Impact factor: 4.677

9.  Diffusion tensor imaging in frontostriatal tracts is associated with executive functioning in very preterm children at 9 years of age.

Authors:  Hanna Kallankari; Virva Saunavaara; Riitta Parkkola; Leena Haataja; Mikko Hallman; Tuula Kaukola
Journal:  Pediatr Radiol       Date:  2020-09-01

Review 10.  Advanced Brain Imaging in Preterm Infants: A Narrative Review of Microstructural and Connectomic Disruption.

Authors:  Philippe Vo Van; Marianne Alison; Baptiste Morel; Jonathan Beck; Nathalie Bednarek; Lucie Hertz-Pannier; Gauthier Loron
Journal:  Children (Basel)       Date:  2022-03-04
  10 in total

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