Congenital structural anomalies, which range from heart defects and cleft palates to urinary tract malformations and limb abnormalities, are among the leading causes of infant morbidity, long-term disability, and pediatric surgical intervention worldwide. While genetics plays a central role, growing evidence points to ambient air pollution and broader environmental conditions as contributors to fetal developmental disruption. Most previous studies, however, have relied on global statistical models that implicitly assume the relationship between environmental exposure and health outcome is the same everywhere. The Yunnan study set out to test that assumption in one of China’s most geographically complex provinces.

The research team, led by Cheng-Hao Zhanghuang and colleagues, first painted a detailed epidemiological portrait of the inpatient cohort. Boys accounted for 67.68 percent of admissions, a male-to-female ratio of roughly 2.1 to 1, and cases were concentrated in early childhood, with toddlers aged one to three years forming the largest group at 28.89 percent. Digestive anomalies were the most common category, representing 29.63 percent of patients, followed by urogenital anomalies at 23.30 percent. Other structural anomalies, musculoskeletal anomalies, and circulatory anomalies made up the remainder. The most frequent individual diagnoses included congenital tongue anomalies, cryptorchidism, and polydactyly. Annual admissions rose steadily from 3,568 in 2014 to a peak of 6,237 in 2019, dipped during 2020, and climbed again to 6,017 by 2024.

To enable robust spatial modeling, the investigators filtered the cohort down to the most frequent conditions within each of five anomaly systems: circulatory, digestive, urogenital, musculoskeletal, and other structural anomalies. This yielded a spatial analysis dataset of 41,531 patients, a step designed to reduce statistical instability caused by counties with sparse case counts. Neurological anomalies were excluded because their numbers at the referral center were too small to support reliable spatial estimates. Patients with multiple anomalies were classified by their principal discharge diagnosis to keep categories mutually exclusive and reduce information bias.

The heart of the study lay in its environmental data assembly. The team compiled eleven county-level environmental and contextual variables averaged over 2014 to 2023, including carbon monoxide, sulfur dioxide, nitrogen dioxide, PM2.5, PM10, ozone, carbon dioxide, land surface temperature, elevation, population density, and the normalized difference vegetation index, a satellite-derived measure of green vegetation cover. Data came from sources such as the National Tibetan Plateau Data Center, NASA Earthdata, the LandScan population dataset, and the Emissions Database for Global Atmospheric Research. Variables with high multicollinearity were removed to ensure that each remaining predictor contributed independent information to the models.

Rather than relying solely on ordinary least squares regression, which produces a single average coefficient for the entire province, the researchers employed geographically weighted regression, or GWR. This technique allows regression coefficients to vary across space, estimating a separate local relationship for each county. Across all five anomaly systems, GWR consistently outperformed the global models, delivering higher coefficients of determination and lower corrected Akaike information criterion and cross-validation values. The authors interpret this as clear evidence of spatial non-stationarity: the strength and even the direction of environmental associations with hospital-based anomaly burden shift across the provincial landscape.

The specific patterns were striking. Carbon monoxide showed predominantly positive associations with referral-weighted institutional burden across anomaly systems, suggesting that counties with higher long-term CO levels tended to contribute more anomaly cases to the referral center. Sulfur dioxide, by contrast, exhibited pronounced spatial heterogeneity, with local coefficients flipping in both magnitude and direction depending on location. Vegetation coverage displayed a consistent negative association across all five systems, hinting that greener counties carried lower institutional anomaly burden, while population density showed positive but geographically variable relationships. The authors emphasize that these are contextual, area-level patterns rather than proof of individual-level causal effects.

Importantly, the researchers are careful about what their data can and cannot show. Because the study draws on a single referral center, the measured burden reflects healthcare-seeking behavior, referral pathways, transportation access, and institutional admission practices, not province-wide prevalence. Remote counties with poor road links or limited referral connections may be underrepresented even if their true anomaly burden is substantial. The lack of individual maternal residential histories also prevented trimester-specific prenatal exposure assessment, and genetic etiologies could not be reliably excluded. The authors explicitly frame their findings as descriptive and hypothesis-generating, requiring validation through population-based registries and multi-center studies before any policy conclusions are drawn.

Nevertheless, the methodological message is clear and potentially far-reaching. In regions marked by complex terrain, uneven economic development, and sharp urban-rural contrasts, one-size-fits-all environmental health models may obscure localized vulnerability. Spatially explicit approaches such as GWR can reveal where environmental associations are strongest, where they weaken, and where they reverse, offering surveillance programs a sharper tool for targeting resources. Proposed biological mechanisms linking prenatal air pollution exposure to congenital anomalies, including oxidative stress, placental dysfunction, and inflammatory disruption of embryonic signaling, remain speculative in this ecological context, but the mapped heterogeneity provides a concrete starting point for future mechanism-oriented investigation.

As congenital anomalies continue to impose a heavy surgical and developmental burden on pediatric health systems worldwide, the Yunnan study adds an important dimension to the evidence base: geography matters. The same pollutant may carry different weight in a densely populated basin than on a remote highland plateau, and greener landscapes may buffer contextual risk in ways that global models cannot capture. Whether these spatial patterns hold up in population-based data from other provinces and countries will determine whether geographically weighted thinking becomes a standard feature of environmental epidemiology for birth defects research.

Subject of Research: Spatial heterogeneity in area-level environmental associations with hospital-based congenital structural anomaly burden in Southwest China

Article Title: Spatial heterogeneity in area-level environmental context of hospital-based congenital structural anomaly burden in Southwest China: a retrospective study from a provincial pediatric referral center

Article References: Zhanghuang, C.-H., Ma, Y.-Y., Zheng, C.-L., Hu, X., Zhang, M.-X., Gao, Y.-P., Chen, J.-R., Yang, S.-W., Zhang, H., Dai, R.-T., Zhang, X.-C., Shen, J., Yan, B., & Wu, J. (2026). Spatial heterogeneity in area-level environmental context of hospital-based congenital structural anomaly burden in Southwest China: a retrospective study from a provincial pediatric referral center. World Journal of Pediatrics. https://doi.org/10.1007/s12519-026-01059-w

Image Credits: AI Generated

DOI: 10.1007/s12519-026-01059-w

Keywords: congenital anomalies, spatial epidemiology, geographically weighted regression, air pollution, Yunnan Province, pediatrics, birth defects, environmental health, vegetation index, population density, referral center, spatial heterogeneity