ORIGINAL RESEARCH article

Front. Psychol., 06 July 2026

Sec. Human Developmental Psychology

Volume 17 - 2026 | https://doi.org/10.3389/fpsyg.2026.1705627

Latent profile analysis and influence factors of neurodevelopment in children aged 0–6 years in China

  • 1. Capital Center for Children’s Health, Capital Medical University, Capital Institute of Pediatrics, Beijing, China

  • 2. Department of Maternal, Child and Adolescent Health, School of Public Health, Tianjin Medical University, Tianjin, China

  • 3. Fengtai Women & Children’s Health Hospital, Beijing, China

  • 4. Fangshan Maternal & Child Health Hospital, Beijing, China

  • 5. Nanmofang Community Health Service Center in Chaoyang District, Beijing, China

  • 6. Beijing Shijingshan Maternal and Child Health Hospital, Beijing, China

  • 7. Beijing Daxing Maternal and Child Care Hospital, Beijing, China

  • 8. Beijing Chaoyang District Maternal and Child Health Care Hospital, Beijing, China

  • 9. The First Affiliated Hospital of Tsinghua University, Beijing, China

  • 10. Haidian District Maternal and Child Health Care Hospital, Beijing, China

  • 11. Beijing Pinggu District Maternal and Child Health Hospital, Beijing, China

  • 12. Huairou District Maternal and Child Health Hospital in Beijing, Beijing, China

Abstract

Background:

Early childhood neuropsychological development exhibits substantial heterogeneity, yet the identification of distinct developmental profiles and their associated risk factors remains limited. This study aimed to delineate latent profiles of neuropsychological development in children aged 0–6 years and to identify factors influencing membership in these profiles.

Methods:

A total of 2,297 children (1,193 boys and 1,104 girls) were assessed using the Chinese Developmental Scale for Children – Second Edition (CDSC-II). Latent profile analysis (LPA) was conducted based on age-standardized percentile scores across five developmental domains: gross motor, fine motor, language, adaptive behavior, and social behavior. A one-step regression mixture model, adjusting for child age and sex, was employed to examine associations between profile membership and perinatal and demographic factors.

Results:

The LPA identified two distinct profiles: an Ordinary Development Group (28.6%, n = 656) and an Excellent Development Group (71.4%, n = 1,641), with good class separation (entropy = 0.818). In the regression mixture model, low birth weight (OR = 2.09, 95% CI: 1.19–3.67), preterm birth (OR = 5.23, 95% CI: 3.14–8.71), maternal pregnancy complications (OR = 1.43, 95% CI: 1.15–1.78), and older child age (OR = 1.02, 95% CI: 1.007–1.027) were significantly associated with increased odds of belonging to the Ordinary Development Group. Conversely, female sex (vs. male) was associated with lower odds of membership in the Ordinary Development Group (OR = 0.67, 95% CI: 0.54–0.83). Parenting style (parental vs. non-parental care) and the timing of complementary food introduction were not independently associated with profile membership after adjustment (both p > 0.05).

Conclusion:

These findings indicate that specific maternal environmental exposures during pregnancy and the postnatal period, alongside child demographic characteristics, are significantly associated with neurodevelopmental outcomes in early childhood. The results underscore the importance of early developmental screening and targeted intervention for high-risk children to optimize long-term neurodevelopmental trajectories.

1 Introduction

Prenatal and early life represent sensitive developmental windows during which environmental exposures can exert lasting effects on health trajectories (). The Developmental Origins of Health and Disease (DOHaD) paradigm emphasizes that early environmental exposures may permanently alter the structure and function of organs and tissues, with lifelong ramifications (; ; ; ). Characterizing these early, potentially modifiable factors offers opportunities for primordial prevention to promote public health.

In recent years, there has been growing interest in the impact of early-life environmental exposures on children’s neurodevelopmental outcomes (). For example, summarized that adverse early-life environments (from fetal period to childhood) are closely related to brain volume, microstructure, and connectivity, while positive early experiences may also buffer some adverse effects across generations. However, several gaps remain in the literature.

First, cultural differences are often overlooked in the assessment of children’s neuropsychological development (). Many neuropsychological assessments used in China are translated or adapted from foreign tools; there is an urgent need for locally developed measures that integrate cultural factors (). The Developmental Scale for Children Aged 0–6 Years (WS/T 580–2017), also known as the Chinese Developmental Scale for Children – Second Edition (CDSC-II) or the Children Neuropsychological and Behavior Scale (CNBS), was developed by the Capital Institute of Pediatrics. (). It is one of the most widely used clinical developmental diagnostic tests in China, with nationwide standardization and good reliability and validity (; ).

Second, previous research has paid insufficient attention to the heterogeneity of children’s neurodevelopment. According to classic developmental theory, development proceeds at varying rates across children and unevenly across different domains within the same child (). Therefore, when assessing neuropsychological development, it is important to consider potential heterogeneity across ability domains. Latent profile analysis (LPA) is a person-centered modeling approach () that classifies individuals into subgroups based on their response patterns on observed variables. LPA maximizes between-group variance and minimizes within-group variance, and statistical indices can be used to evaluate classification accuracy (; ). As such, LPA is well-suited to describe heterogeneity in child neurodevelopment (; ; ).

Third, prior studies have often ignored the measurement error associated with latent class assignments when examining associations between exposures and outcomes. Some studies have used traditional regression models based on most-likely class assignments derived from LPA (; ; ). However, this approach is suboptimal because latent class membership is estimated with uncertainty; treating it as observed introduces bias and typically underestimates the strength of associations (; ). Regression mixture models (RMMs) address this issue by estimating covariate effects and latent class membership simultaneously within a one-step framework, thereby properly accounting for classification uncertainty (; ; ).

Therefore, this study used the CDSC-II (CNBS), a locally developed neurodevelopmental assessment tool, to conduct latent profile analysis on the age-standardized percentile scores of five developmental domains (gross motor, fine motor, language, adaptive behavior, social behavior) in a sample of Chinese children aged 0–6 years. We then used a one-step regression mixture model to examine the associations between early-life environmental exposures and latent class membership.

We hypothesized that: (a) there are individual differences in children’s neuropsychological development, reflected in varying levels across domains; (b) child age and sex may be associated with class membership; and (c) specific early-life exposures (e.g., low birth weight, preterm birth, maternal pregnancy complications, feeding practices, parenting style) are associated with membership in a lower-functioning developmental profile.

2 Methods

2.1 Research object

The study’s dataset originates from cohort databases titled “Study on Developmental Behavior Monitoring and Appropriate Techniques for Intervention of Developmental Abnormalities in Children Aged 0–6 Years.” and “A study on self-assessment tools for psychological and behavioral development of 0–3 year old children”. Utilizing a multi-stage stratified random sampling technique, the study selected four districts in urban and four in suburban areas, comprising two urban and two suburban districts each. Within these districts, 1–2 medical and health service institutions were chosen as testing sites. The study analyzed test data collected from May 1, 2021, to August 30, 2023.

Inclusion Criteria: Children aged 0–6 years residing in the jurisdiction for more than half a year were included.

Exclusion Criteria: Children suffering from various acute and chronic diseases were excluded.

The study was conducted under The Code of Ethics of the World Medical Association (Declaration of Helsinki). We also obtained informed consent from the parents of participants before the experiment. The research protocol of this study was approved by the Medical Ethics Committee of Capital Institute of Paediatrics (KSSHERLLM 2017013; SHERLL2022051).

2.2 Neuropsychological development assessment tools

Children’s neuropsychological development was evaluated using the WS/T 580–2017 Developmental Behavior Assessment Scale for 0–6 Year-Old Children (CDSC-II) — also referred to as the Children Neuropsychological and Behavior Scale (CNBS) — and a standardized toolbox. This assessment covers five developmental domains: gross motor, fine motor, language, adaptive ability, and social behavior, along with the total development quotient. Developmental Age (DA) and Developmental Quotient (DQ) represent the child’s developmental level, with a DQ score below 70 indicating developmental delay, between 70 and 80 suggesting a developmental delay margin, and scores of 80 or above signifying normal development ().

2.3 Investigation of environmental exposures in early life

A custom questionnaire was employed for a face-to-face survey covering various aspects:

Demographic Data: Parental information including household registration, address, age, occupation, and monthly income per capita.

Maternal Health: Pregnancy history, diseases during pregnancy (e.g., anemia, hypertension), and maternal weight gain.

Birth Status: Delivery mode, gestational week, and newborn measurements (length, weight, head circumference).

Disease Status at Birth: Conditions like intrauterine asphyxia, birth distress, and hypoxic–ischemic encephalopathy.

Neonatal Period: Incidences of hyperbilirubinemia, neonatal pneumonia, sepsis, convulsions, and allergies.

Infancy Feeding: Primary caregivers, breastfeeding status, and introduction of complementary foods (binary yes/no indicator, assessing whether solid or semi-solid foods had been introduced at the time of assessment).

2.4 Introduction to latent profile analysis (LPA)

LPA is a modeling approach that identifies potential subpopulations within a population based on specific variables. The LPA assumes that people can, with varying degrees of probability, be typed into categories (subgroups) with different personal and/or environmental attributes of the configuration profile (). Its equation model can be expressed as:

Where μik and σik represent profile-specific (k) means and variances for variable i, and πk indicates profile density, or the proportion of N participants that belong to profile k. LPAs assume (a) samples drawn from a heterogeneous population produce data that are a mixture of K profile-specific distributions; (b) observed y indicator variables are distributed normally; and (c) the profile-specific mean vectors μk are the profile-specific (k) observed variable means. If both separate mean vectors (μk) and separate covariance matrices (∑k) were freely estimated for all of the K latent profiles contained in the data, the number of estimated parameters would quickly increase as the number of observed variables increased. Therefore, two LPA model constraints are commonly imposed to model y variation with a minimum number of estimated parameters. First, the local independence assumption states that conditional on correct latent profile extraction or enumeration, all y are uncorrelated within each k latent profile, and all k-specific off-diagonal covariance matrix elements are zero. Second, the homogeneity assumption states that profile-specific covariance matrix elements along the main diagonal are constrained to equality across all k (i.e., ∑k = ∑, Yk ~ N [μk, ∑]) (; ). Together, local independence and homogeneity assume that latent profile-specific (k) covariance matrices are diagonal and homogeneous, and that latent profiles differ only in their y-variable location (μk), not their y-variable relationship form (∑) ().

2.5 Statistical methods

Latent profile analysis indicators. The LPA was conducted using five continuous indicator variables corresponding to the five developmental domains of the CDSC-II: (1) Gross Motor, (2) Fine Motor, (3) Language, (4) Adaptive Behavior, and (5) Social Behavior. For each domain, we used age-standardized percentile scores as the LPA indicators. These percentile scores were calculated according to the CDSC-II manual as: Percentile = (child’s Developmental Age (DA) score / same-age-group normative median DA) × 100. The CDSC-II provides age-specific normative data based on a nationally representative sample of Chinese children (WS/T 580–2017). Using percentile scores rather than raw scores minimizes age-related confounding, as the percentiles represent a child’s relative standing within their age cohort.

Data preprocessing and missing data. The five domain percentile scores were examined for distributional properties prior to analysis. Skewness and kurtosis statistics were within acceptable ranges (all |skewness| < 2.0, |kurtosis| < 7.0), supporting the use of maximum likelihood estimation without additional transformations. Missing data were handled using Full Information Maximum Likelihood (FIML) under the assumption of missing at random (MAR). FIML uses all available information in the dataset to estimate model parameters without listwise deletion, preserving statistical power and producing unbiased estimates under MAR. The overall missing data rate was low (<5% for any single variable).

Model selection. We fitted models with 1 through 5 latent classes. Model fit was evaluated using the Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), adjusted BIC (aBIC), entropy, the Lo–Mendell–Rubin (LMR) test, and the bootstrap likelihood ratio test (BLRT). The smallest class proportion was also considered to ensure interpretability.

Regression mixture model. To examine the association between early-life environmental exposures and latent class membership, we used a one-step regression mixture model (RMM). In this approach, covariates predicting class membership are included directly in the LPA model, and all parameters (class profiles, class proportions, and covariate effects on class membership) are estimated simultaneously. This one-step approach properly accounts for classification uncertainty because it does not require assigning individuals to discrete classes; instead, each individual contributes to the likelihood through their posterior class membership probabilities ().

The following covariates were included as simultaneous predictors of class membership, selected a priori based on the Developmental Origins of Health and Disease (DOHaD) framework and prior literature: low birth weight (yes/no), preterm birth (yes/no), maternal history of pregnancy complications (yes/no), introduction of complementary foods (yes/no), parenting style (parental care vs. non-parental care), child sex (male/female), and child age in months (continuous). No stepwise variable selection was used; all covariates were entered together.

Software and estimation. All analyses were conducted using Mplus version 8.3 (). The LPA models were estimated using robust maximum likelihood (MLR) estimation, which provides standard errors and model fit statistics robust to non-normality and non-independence of observations. Missing data were handled using FIML as noted above. Model convergence was assessed by examining the iteration history and ensuring that all models achieved proper convergence (no negative residual variances, all parameter convergence criteria met). The key Mplus syntax for the final 2-class regression mixture model is provided in the Supplementary Appendix S1.

Quality control. Before the project was carried out, the test personnel were uniformly trained and all passed the assessment. Before carrying out the actual inspection, the Child Health Care Practitioners (CHCP) were taught in practice, and a consistency test was carried out on the CHCP. Only when the consistency test reached 0.9 could they obtain the qualification of the main examination. The data in this study were entered into the growth and development test system for children aged 0–6 years, and a new data set was established for analysis after logical error checking.

3 Results

3.1 Basic demographic characteristics of the subjects

A total of 2,297 children were included in the study, including 1,193 boys and 1,104 girls. The mean age of fathers was 32.29 ± 3.91 years, and that of mothers was 30.82 ± 2.92 years. Detailed information can be found in Table 1. The early environmental exposure of the subjects can be seen in Table 2.

Table 1

VariablesN%
Gender
 Male1,19351.93
 Female1,10448.07
Household registration
 Beijing1,56167.96
 Non-Beijing73632.04
Ethnicity
 Han2,24497.69
 Non-Han532.31
Father’s age
 21.0–29.9 years39217.07
 30.0–39.9 years177477.23
 ≥40 years1315.70
Mother’s age
 21.0–29.9 years85037.00
 30.0–39.9 years1,40361.80
 ≥40 years441.20
Father’s occupational type
 Government or Organizational Manager2199.53
 Professional/Technical Personnel40117.46
 Administrative or Relevant Personnel81935.66
 Social Production, Service, or Life Service Personnel1898.23
 Agriculture, Forestry, Livestock, Fishing, or Auxiliary Personnel1627.05
 Other Unclassified Occupations50722.07
Mother’s occupational type
 Government or Organizational Manager1978.58
 Professional/Technical Personnel36015.67
 Administrative or Relevant Personnel82135.74
 Social Production, Service, or Life Service Personnel1707.40
 Agriculture, Forestry, Livestock, Fishing, or Auxiliary Personnel160.70
 Other Unclassified Occupations73331.91
Father’s education level
 High School or Below1516.57
 Associate’s Degree or Bachelor’s Degree194084.46
 Graduate Degree2068.97
Mother’s education level
 High School or Below1305.66
 Associate’s Degree or Bachelor’s Degree197185.81
 Graduate Degree1968.53
Family per capita monthly income level
 < 3,000 yuan50.22
 3001–5,000 yuan96141.84
 >5,000 yuan1,33157.94

Sociodemographic characteristics of the sample (N = 2,297).

Table 2

VariablesN%
Maternal anemia during pregnancy
 Yes964.18
 No2,20195.82
Maternal hypertension during pregnancy
 Yes994.31
 No2,19895.69
Gestational diabetes
 Yes24810.8
 No204989.2
Maternal hypothyroidism
 Yes783.4
 No2,21996.6
Maternal abnormal pregnancy history *
 Abnormal Pregnancy1506.53
 Normal Pregnancy2,14793.47
Delivery method
 Natural Vaginal Delivery1,39060.51
 Cesarean Section86037.44
 Forceps or Vacuum Extraction472.05
Premature birth
 < 28 weeks70.3
 28-32 weeks331.44
 33-37 weeks30813.41
 ≥ 38 weeks194984.85
Birth weight
 < 1,500 g371.61
 1500-2499 g1054.57
 2500-4000 g202288.03
 > 4,000 g1335.79
Birth length
 < 46 cm873.79
 46-52 cm2,12792.6
 > 52 cm833.61
Head circumference at birth
 < 32 cm682.96
 32–34 cm173475.49
 > 34 cm49521.55
Singleton or multiple births
 Singleton2,25598.17
 Twins or Multiples421.83
Oxygen deprivation during delivery
 Yes200.87
 No2,27799.13
Severe neonatal infection
 Yes492.13
 No2,24897.87
Hyperbilirubinemia
 Yes903.92
 No2,20796.08
History of allergies
 Yes913.96
 No2,20696.04
Feeding method
 Breastfeeding1,22453.29
 Formula Feeding56724.68
 Mixed Feeding50622.03
Introduction of solid food
 Yes2,10791.73
 No1908.27
Primary caregiver
 Parents1,08647.28
 Grandparents1,19251.89
 Babysitters190.83

Early-life environmental exposures (N = 2,297).

*Abnormal pregnancy history includes previous miscarriage (first trimester), stillbirth (second trimester), or history of stillbirth during pregnancy or childbirth.

3.2 Latent profile analysis results

The study extracted 1–5 latent class models, and the fit indices are shown in Table 3. The AIC, BIC, and aBIC decreased monotonically with the increase in the number of classes. However, we selected the 2-class solution based on a balanced consideration of multiple criteria:

  • Statistical adequacy: The 2-class model demonstrated the highest entropy (0.805) among all fitted models, indicating the clearest class separation. While the 3-class model had a significant LMR test (p < 0.001), its entropy was slightly lower (0.797). The 4-class model showed a non-significant LMR test (p = 0.0245, which does not survive multiple comparison correction), and its entropy dropped further (0.755). The 5-class model had a non-significant LMR test (p = 0.299) and an unacceptably small smallest class size (4.18%, n ≈ 96), which would compromise statistical power and generalizability.

  • Parsimony and interpretability: The 2-class solution distinguished children with generally excellent development from those with ordinary development, which has clear clinical and public health relevance. The higher-order solutions (3-, 4-, and 5-class) did not reveal qualitatively distinct developmental profiles (i.e., domain-specific patterns); rather, they primarily subdivided the Ordinary Development Group by severity level, producing classes that differed quantitatively but not qualitatively in their domain patterns. Results for the 3- and 4-class solutions are provided in Supplementary Table S1 and Supplementary Figures S1–S3.

  • Replicability: Given that our sample, while large, was drawn from specific geographic regions, we prioritized a more parsimonious solution that would be more likely to replicate in independent samples.

Table 3

ModelAICBICaBICentropyLMRPB_LRTPMMCN
14,3924,4494,417
220772,1692,1180.805−2186.065<0.001−2186.065<0.00128.57%
31726185217820.797−1022.896<0.001−1022.896<0.00118.86%
41,2961,4571,3680.755−841.3060.0245−841.306<0.00117.42%
51,2221,4181,3100.758−620.3380.2989−620.338<0.0014.18%

Model fit indexes.

AIC, Akaike Information Criterion; BIC, Bayesian Information Criterion; aBIC, Adjusted Bayesian Information Criterion; BLRT, Bootstrap Likelihood Ratio Test; aLMR, Adjusted Lo–Mendell–Rubin Likelihood Ratio Test.

Based on the mean percentile scores of the five developmental domains across classes (Figure 1), the boundary between the two classes was clear. Because Class 1 had lower mean scores across all five domains compared to Class 2, we named Class 1 “Ordinary Development Group” and Class 2 “Excellent Development Group.” In the whole sample, 656 children (28.56%) were classified into the Ordinary Development Group, and 1,641 children (71.44%) were classified into the Excellent Development Group. In the Ordinary Development Group, there were 379 boys and 277 girls; in the Excellent Development Group, there were 814 boys and 827 girls. There was a statistically significant difference in the sex ratio between the two groups (χ2 = 12.533, p < 0.001), with a higher proportion of girls in the Excellent Development Group. The average age of the Excellent Development Group was 12.933 ± 10.065 months, and the average age of the Ordinary Development Group was 14.964 ± 12.115 months; the difference was statistically significant (t = 4.173, p < 0.001). Results for the 3-class and 4-class solutions, including model fit indices and profile plots, are provided in Supplementary Table S2–S3 and Supplementary Figures S1–S2.

Figure 1

3.3 Factors associated with different profiles of neuropsychological development

According to the one-step regression mixture model adjusting for child age and sex, five early-life factors were significantly associated with latent class membership. Low birth weight (OR = 2.09, 95% CI: 1.19–3.67), preterm birth (OR = 5.23, 95% CI: 3.14–8.71), maternal pregnancy complications (OR = 1.43, 95% CI: 1.15–1.78), and older child age (OR = 1.02, 95% CI: 1.007–1.027) were associated with higher odds of belonging to the Ordinary Development Group. Female sex (vs. male) was associated with lower odds of belonging to the Ordinary Development Group (OR = 0.67, 95% CI: 0.54–0.83). Parental care (parental vs. non-parental care) and introduction of complementary foods were not independently associated after adjustment (both p > 0.05). Detailed results are presented in Table 4.

Table 4

VariableBS.E.POR95% CI
Low birth weight (yes vs. no)0.7380.2860.0102.09[1.19, 3.67]
Preterm birth (yes vs. no)1.6540.260<0.0015.23[3.14, 8.71]
Maternal pregnancy complications (yes vs. no)0.3550.1120.0021.43[1.15, 1.78]
Introduction of complementary foods (yes vs. no)−0.1110.1180.3480.90[0.71, 1.13]
Parental care (vs. non-parental care)−0.4570.2620.0810.63[0.38, 1.06]
Child age (months, continuous)0.0170.0050.0011.02[1.007, 1.027]
Sex (female vs. male)−0.4070.111<0.0010.67[0.54, 0.83]

Early environmental factors associated with membership in the Ordinary Development Group (vs. Excellent Development Group) (N = 2,297).

*In this study, the mother’s history of abnormal pregnancy was a previous pregnancy, a miscarriage in the first trimester, a stillbirth in the second trimester, or a history of stillbirth during pregnancy or even childbirth. OR > 1 indicates higher odds of being in the Ordinary Development Group. All covariates were entered simultaneously. The model adjusted for child age and sex. B = log odds ratio; S.E. = standard error; CI = confidence interval.

4 Discussion

In recent years, both domestic and international scholars have increasingly adopted person-centered latent variable models to explore heterogeneity in child psychological development (; ). Notably, research utilizing the Bayley Scales of Infant and Toddler Development to analyze very preterm infants has identified four latent neurobehavioral development profiles, a methodology paralleling our study (). Additionally, Dong et al. have delineated high-level and low-level development groups in their exploration of committed compliance developmental trajectories in typically developing infants and toddlers, employing latent class analysis ().

Our investigation using the five age-standardized percentile scores of the CDSC-II revealed a two-class solution that primarily reflects a quantitative differentiation in overall developmental level rather than qualitatively distinct neurodevelopmental profiles. Specifically, children in the Excellent Development Group (Class 2) scored higher across all five domains (gross motor, fine motor, language, adaptive behavior, social behavior) in a relatively parallel fashion, while children in the Ordinary Development Group (Class 1) showed uniformly lower scores. This pattern is consistent with a general neurodevelopmental factor or g-factor commonly observed in developmental research (), where individual differences in development tend to be correlated across domains. Consequently, while latent profile analysis provides a useful empirical classification of children into higher- and lower-functioning groups, our findings do not support the existence of distinct, domain-specific developmental phenotypes in this sample. Future research with larger and more diverse samples may be needed to identify whether such phenotypes exist in the broader population of Chinese children.

This study indicates significant differences in gender and age within the two latent profiles, with girls more frequently classified into the Excellent Development Group, where members were also younger on average compared to the Ordinary Development Group. Historical research on gender’s impact on child psychological development suggests preschool-age girls possess advantages over boys in motor skills, language, and social development (; ; ; ). One proposed rationale involves differences in brain hemisphere development, affording girls inherent advantages in brain functional organization, notably in language and communication (). However, the argument that gender differences in brain structure or function directly translate to language task performance disparities remains contested, with limited evidence supporting gender differences in brain and language development (). This study’s observations of gender differences in neurodevelopment require further substantiation through additional methods, such as cognitive neurotests sensitive to age or neuroimaging studies.

Internationally, the consensus is that a child’s cognitive neurodevelopment level correlates strongly with age, reflected in age-specific developmental milestones (). Although we used age-standardized percentile scores to minimize the effect of age, the two latent classes still differed significantly in mean age, with the Ordinary Development Group being older on average (14.96 ± 12.12 months) than the Excellent Development Group (12.93 ± 10.07 months). This residual age difference may reflect the fact that older children in our sample had a broader range of developmental experiences and environmental exposures that could influence their relative standing on age-standardized percentiles. As a sensitivity analysis, we included child age (in months) as a continuous covariate in the regression mixture model; the pattern of associations remained substantively unchanged (see Supplementary Table S4), suggesting that the reported ORs are not primarily driven by age confounding. Nevertheless, residual age effects cannot be completely ruled out.

Using a one-step regression mixture model (RMM) that properly accounts for classification uncertainty, our study confirms associations (instead of causal effects) between several early-life environmental exposures and neurodevelopmental class membership. Children exposed to low birth weight (OR = 2.09, 95% CI: 1.19–3.67), preterm birth (OR = 5.23, 95% CI: 3.14–8.71), or maternal pregnancy complications (OR = 1.43, 95% CI: 1.15–1.78) were more likely to belong to the Ordinary Development Group. These findings align with the bulk of similar research (; ). Adverse maternal pregnancy histories may relate to a range of factors, including genetic inheritance, nutritional factors, and occupational hazards during pregnancy, potentially limiting offspring psychological development (; ; ; ; ). Further research is needed to elucidate the mechanisms through which maternal pregnancy history impacts offspring neurodevelopment.

In contrast, after adjustment for other covariates, introduction of complementary foods (OR = 0.90, 95% CI: 0.71–1.13, p = 0.348) and parental care (parental vs. non-parental caregiving; OR = 0.63, 95% CI: 0.38–1.06, p = 0.081) were not significantly associated with class membership. Although the direction of the point estimates is consistent with a potential protective effect, the lack of statistical significance suggests that these factors may not be independent contributors in this sample, or that their effects are smaller and require larger studies to be detected. Previous literature has emphasized the importance of early feeding practices and responsive caregiving (; ; ; ; ), but our null findings may reflect residual confounding by unmeasured socioeconomic or caregiving factors, or the relatively high baseline levels of complementary feeding (91.7% of children had been introduced) limiting variability.

Nevertheless, this study has several important limitations. First, the cross-sectional design precludes establishing temporal order between exposures and developmental outcomes. While many of the exposures (e.g., low birth weight, preterm birth, maternal pregnancy history) temporally precede the developmental assessment, other variables — notably parenting style and complementary food introduction — were measured concurrently with the outcome. Reverse causality is therefore a plausible alternative explanation: children with observed or perceived developmental delays may elicit different caregiving responses (e.g., more intensive parental care, altered feeding practices) or may be more likely to remain in parental rather than non-parental care. The direction of association cannot be determined from these data.

Second, several key exposures were assessed via caregiver self-report, which is subject to recall bias and social desirability bias. Maternal pregnancy history and feeding practices may be particularly susceptible to inaccurate recall, especially in a sample with children aged 0–6 years.

Third, the absence of detailed socioeconomic data (parental education, household income, maternal mental health) limits our ability to control for confounding by socioeconomic status, which is strongly associated with both the exposures and child developmental outcomes.

Fourth, although we used age-standardized percentile scores and additionally adjusted for age in sensitivity analyses, the significant age difference between the two latent classes suggests that residual age confounding cannot be entirely ruled out.

Fifth, as noted above, the two-class solution identified represents a general developmental level gradient rather than qualitatively distinct phenotypes. Whether this reflects a true absence of domain-specific profiles in this population or limited statistical power to detect such profiles deserves further investigation in larger samples.

5 Conclusion

This study used latent profile analysis to identify two quantitative developmental groups (Excellent and Ordinary) among Chinese children aged 0–6 years based on age-standardized CDSC-II percentile scores. The two-class solution primarily reflects a general developmental level gradient rather than qualitatively distinct domain-specific profiles. A one-step regression mixture model revealed that low birth weight, preterm birth, maternal pregnancy complications, older child age, and male sex are significantly associated with membership in the Ordinary Development Group (Class 1). Parenting style and introduction of complementary foods were not independently associated after adjustment. These findings highlight the importance of early-life environmental exposures for child neurodevelopment. However, due to the cross-sectional design, causal inferences cannot be drawn. Future longitudinal studies are needed to clarify the temporal direction of these associations and to explore whether more nuanced developmental phenotypes emerge in larger, more diverse samples.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Medical Ethics Committee of Capital Institute of Paediatrics (KSSHERLLM 2017013; SHERLL2022051). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

RL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Writing – original draft. RZ: Conceptualization, Formal analysis, Validation, Writing – original draft. XZ: Data curation, Investigation, Writing – original draft. HW: Data curation, Investigation, Writing – original draft. YuqG: Data curation, Investigation, Writing – original draft. CC: Data curation, Investigation, Writing – original draft. HY: Data curation, Investigation, Writing – original draft. LB: Data curation, Investigation, Writing – original draft. QL: Data curation, Investigation, Writing – original draft. JL: Data curation, Investigation, Writing – original draft. YuaG: Data curation, Investigation, Writing – original draft. XL: Data curation, Investigation, Writing – original draft. SZ: Data curation, Investigation, Writing – original draft. YZ: Data curation, Investigation, Writing – original draft. LG: Conceptualization, Formal analysis, Funding acquisition, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study received research grants from the China Disabled Persons Federation (Grant number: 2023CDPFHS-06) and the Capital Institute of Pediatrics Innovation Transformation Project (CX-2022-05). The sponsors had no role in the study design, data collection, analysis, interpretation, or manuscript preparation.

Acknowledgments

The authors thank the infants/children, parents and schools for their cooperation.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyg.2026.1705627/full#supplementary-material

References

Summary

Keywords

Chinese children, early-life environmental exposures, influencing factor, latent profile analysis, neuropsychological development

Citation

Li R, Zhang R, Zheng X, Wang H, Guo Y, Chen C, Yu H, Bai L, Li Q, Luo J, Guo Y, Li X, Zhang S, Zhang Y and Gao L (2026) Latent profile analysis and influence factors of neurodevelopment in children aged 0–6 years in China. Front. Psychol. 17:1705627. doi: 10.3389/fpsyg.2026.1705627

Received

15 September 2025

Revised

09 June 2026

Accepted

16 June 2026

Published

06 July 2026

Volume

17 - 2026

Edited by

Nataliia Kozhemiako, Brigham and Women's Hospital and Harvard Medical School, United States

Reviewed by

Adonay S. Nunes, Massachusetts General Hospital and Harvard Medical School, United States

Daigo Kato, Takatsuki General Hospital, Japan

Updates

Copyright

*Correspondence: Lei Gao,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics