Geographical variations in maternal lifestyles during pregnancy associated with congenital heart defects among live births in Shaanxi province, Northwestern China

In this study, we aimed to explore regional differences in maternal lifestyle during pregnancy related to congenital heart defects (CHD) in Shaanxi province, Northwestern China. A large-scale epidemiologic survey of birth defects among infants born during 2010–2013, was conducted in Shaanxi province. Non-spatial and geographic weighted logistic regression models were used for analysis. The spatial model indicated that passive smoking frequency was positively associated with CHD for 43.3% of participants (P < 0.05), with the highest OR in North Shaanxi and the lowest in South Shaanxi. Approximately 49.2% of all mothers who ever drink tea were more likely to have an infant with CHD (P < 0.05), with the highest OR values observed in North and Central Shaanxi. Additionally, maternal alcohol intake frequency ≥ 1/week was significantly correlated with CHD among about 24.7% of participants (P < 0.05), with OR values ranging from 0.738 (Central Shaanxi) to 1.198 (North Shaanxi). The rates of unhealthy maternal lifestyles during pregnancy associated with CHD differed in various areas of the province. The role of geographical variations in these factors may provide some possible clues and basis for tailoring site-specific intervention strategies.

Scientific RepoRtS | (2020) 10:12958 | https://doi.org/10.1038/s41598-020-69788-0 www.nature.com/scientificreports/ for about half of the population in northwest China, is facing more serious BD than other regions of China 17 . Simultaneously, due to extremely unbalanced economic development and huge differences in socio-culture among Shaanxi regions, maternal lifestyle during pregnancy is inevitably disparate across different surveyed areas, and spatial heterogeneity may play an important role in the correlation between the incidence of CHD and maternal lifestyles during pregnancy. Therefore, we hypothesized that the relationship between maternal lifestyle during pregnancy and CHD varied significantly from region to region in Shaanxi. In most previous studies, however, conventional regression models were used to produce average parameters to explore the relationship between maternal lifestyle during pregnancy and CHD over the whole study area, but spatial heterogeneity was neglected [17][18][19][20] . In the present study, we used the Geographically Weighted Logistic Regression (GWLR) model to generate local coefficients to account for geographical variations in evaluating the correlation between CHD in infants and maternal lifestyles during pregnancy. GWLR allows for the examination of local changes and provides a map for the visualization of spatial heterogeneity 21 . This geographically weighted regression model was helpful in identifying distinct patterns that appeared in different villages and townships, so that we could clearly visualize differences in maternal lifestyles during pregnancy related to CHD among disparate regions. Thus, the current study aimed to (1) describe the geographic distribution of maternal lifestyles during pregnancy across different regions, (2) and explore the spatially varying relationship between maternal lifestyles during pregnancy and CHD adjusting for other possible confounding factors.

Results
Baseline characteristics of the participants. In this study, a total of 1605 villages were chosen within 272 townships of 30 counties in Shaanxi province, and 29,098 live infants were included in the final analysis. Among the infants included in the study, average age was 16.9 ± 11.3 months (range 0-46 months); 54.4% were male; 98.8% were singletons. Among the mothers included in the study, average age was 28.01 ± 4.86 years; average parity was 1.65 ± 0.82. The proportions of study participants from northern, central, and southern Shaanxi, were 25.7%, 53.9% and 20.3%, respectively. The baseline characteristics of the participants across different areas in Shaanxi province are presented in Table 1.
The overall prevalence of CHD among the study respondents was 0.760%, with a range from 0 to 4.4% among Shaanxi Province. Notably, the rate of CHD among subjects was higher in South Shaanxi (1.1%) than the corresponding rate in North Shaanxi (0.5%), and Central Shaanxi (0.8%) ( Supplementary Fig. S1). The rate of CHD among infants increased markedly in association with increases in fetal number or infant parity (P for trend < 0.001). A family history of CHD was associated with increased incidence of CHD among infants (P for trend < 0.001). The results obtained showed that the prevalence of CHD increased with maternal age (P for trend < 0.01). The rate of CHD was higher in the floating population than in the permanent population (P = 0.001). The results obtained showed that maternal lifestyle habits during pregnancy, such as second-hand smoking (P for trend < 0.01) and tea consumption (P for trend < 0.01), may increase risk for CHD.
The proportion of mothers with a history of passive smoking in pregnancy was highest in North Shaanxi (32.1%), followed by Central Shaanxi (20.8%), and South Shaanxi (20.4%) ( Supplementary Fig. S2). With regard to the consumption of alcohol, tea and coffee during pregnancy, South Shaanxi was the area with the highest rate (3.0%, 6.5% and 1.4%, respectively), compared to Central Shaanxi (1.0%, 1.6% and 0.8%, respectively), and North Shaanxi (0.5%, 0.9% and 0.4%, respectively) (Supplementary Fig. S3-S5). non-spatial logistic regression. Table 2 summarized the significant correlations between CHD and maternal lifestyle-related factors using a logistic regression model adjusting for other possible risk factors. Compared with single fetal births, twin and multi-fetal births were associated with a substantial increase in CHD incidence (OR 5.190,). CHD was more likely to occur in the children of mothers ≥ 30 years of age, compared with children of mothers 18-24 years of age (OR 2.353, 95% CI 1.338-4.140). The rate of CHD was higher in the floating population (OR 1.933, 95% CI 1.261-2.963) than in the permanent population. Participants who had a family history of CHD, were more likely to suffer from CHD (OR 11.250, 95%CI 4.675-27.075). Compared to those living in South Shaanxi, those living in Central Shaanxi and North Shaanxi had decreased risk for CHD (OR 0.440, 95%CI 0.251-0.773 and OR 0.591, 95%CI 0.382-0.916, respectively).
Compared to infants born to mothers without a history of passive smoking during pregnancy, infants born to women with a history of passive smoking during pregnancy had increased risk for CHD (OR 2.449, 95%CI 1.691-3.548). However, other maternal lifestyles during pregnancy, including the frequency of alcohol intake, and consumption of tea and coffee, were not significantly associated with CHD in surveyed areas.
Multivariate spatial logistic regression. Based on the results of GWLR, we determined that the optimal bandwidth size was 27,875,000 (persons). When the results of global regression model fitting were viewed in combination with the results from GWLR, we concluded that the geographically weighted regression had better fitness. The effects of spatial variation on local estimates for CHD risk factors, as determined using the GWLR model are presented in Table 3.
With regard to maternal lifestyles during pregnancy, the frequency of passive smoking was found to be positively associated with CHD for 43.3% of participants in South and North Shaanxi province (P < 0.05), with the highest OR observed in North Shaanxi and the lowest OR observed in South Shaanxi (Table 3 and Fig. 1). Approximately 49.2% of all mothers who ever drank tea during pregnancy were more likely to have an infant with CHD, compared to mothers who never drank tea during pregnancy (P < 0.05). For this trend, OR values were highest in North and Central Shaanxi (Table 3 and Fig. 2). Maternal alcohol intake more than once per week was significantly correlated with CHD occurrence among 24.7% of participants (P < 0.05), all of whom Table 1. Baseline characteristics of study population by different areas in Shaanxi province. Values are given as number and proportion of the study population; differences in socio-demographic characteristics between 2009 and 2013 were tested using χ 2 tests.

Discussion
Using the results from a large-scale cross-sectional survey conducted in Shaanxi province in northwestern China, we explored the geographic variation of some maternal lifestyle-related factors during pregnancy in relation to the risk of CHD after controlling for other baseline characteristics using spatial regression models. The results of a study conducted previously in China indicated that the rate of passive smoking among pregnant women ranged from 26.6% to 45.3%, which were 17.1% in Shanghai, 34.88% in Shenzhen and 33.96% in Zhengzhou respectively 22 . Similarly, our study showed that the proportion of mothers with a history of passive smoking in pregnancy was 23.6% (range 5.5-57.5%) in Shaanxi, which was highest in North Shaanxi (32.1%). These results showed a significant regional variation in Shaanxi. Simultaneously, the rates of the maternal consumption of alcohol, tea and coffee during pregnancy were higher in South Shaanxi in comparison to other   www.nature.com/scientificreports/ areas of Shaanxi. These findings were consistent with regional distribution of CHD incidence, which showed a decreasing trend from north to south in Shaanxi province. Global regression analysis revealed that mothers who had a history of passive smoking during pregnancy were more likely to give birth to a baby with CHD. Studies conducted in the UK and the USA reported similar findings 23,24 . The GWLR model showed that the OR for passive smoking during pregnancy was highest in North Shaanxi and lowest in Central Shaanxi. Our survey data indicated that passive smoking was higher in North Shaanxi than in Central Shaanxi or South Shaanxi. Most parts of Central Shaanxi are developed urban areas (e.g., Xi'an), where tobacco control measures have been implemented. In comparison with other areas of Shaanxi, Central Shaanxi has more public places that have become true smoke-free environments. Second-hand smoke may therefore be less of an issue for pregnant women in Central Shaanxi. In contrast, the area surveyed in North Shaanxi comprise small townships and villages, where efforts at tobacco control have not been sufficient to reduce the harmful effects of tobacco. One research study conducted in Shaanxi found that a high proportion of participants, especially less educated mothers in North Shaanxi and South Shaanxi, were exposed to environmental tobacco smoke for more than 15 min on at least 1 day per week 17 . These results remind us that it is extremely urgent to take adequate measures to control tobacco smoking, especially in North Shaanxi. Additional studies will be required to corroborate these findings and to elucidate possible reasons for this geographical variations.
Numerous previous studies have found that tea consumption during pregnancy could adversely affect fetal growth and development [25][26][27] . In our study, the non-spatial logistic regression models showed no significant association between tea consumption during pregnancy and CHD. However, the GWLR model revealed that tea consumption during pregnancy may increase the risk for CHD in 49.2% of all participants (mostly in Central and South Shaanxi province). When we focus on the map drawn by OR values, it is easy to see that OR values www.nature.com/scientificreports/ were highest in areas near the eastern part of Shaanxi, with a trend of gradual decline from east to west. On account of the high-altitude landform and warm and humid climate environment in South Shaanxi, the region is suitable for the growth of tea. Local residents therefore consume more tea than inhabitants of other areas. The results of our survey showed that the rate of tea consumption was highest in southeast Shaanxi (6.5%), followed by Central Shaanxi (1.6%), and then by North Shaanxi (0.9%). Household economic level may play a role in tea consumption. The socioeconomic stratification of residents in Shaanxi province has been described in previous studies 17,18 . As a consumer's income increases, he or she is more likely to drink tea 28 . However, more research is needed to unravel the geographic patterns modulating the harmful effects of tea in relation to risk for CHD. Although the non-global regression results showed no significant correlation between maternal history of alcohol intake and the likelihood of having a child with CHD, the spatial regression model revealed a positive association between maternal history of alcohol intake and the risk of CHD for 24.7% of mothers residing in North Shaanxi. On the map used to visualize these results, the areas with the highest OR values for CHD were mainly distributed in North Shaanxi, followed by South and Central Shaanxi. The relationship between maternal alcohol consumption and CHD has long been studied, with conflicting results. Studies in animal models have identified a significant association between maternal exposure to alcohol during pregnancy and an increased incidence of heart anomalies in offspring 29,30 . Research conducted by Polygenis et al. revealed no effect of maternal alcohol consumption on fetal malformations 31 . A review by Henderson et al. found that only one study reported a significant increase in risk for fetal malformations among mothers who consumed alcohol during pregnancy 32 . Geographical variations may explain, at least in part, the inconsistent results reported for the associations between maternal exposure to alcohol during pregnancy and CHD.
However, several limitations should be considered when interpreting our results. Because this research was based on a cross-sectional study of the population, the causal relationship between maternal lifestyle during pregnancy and CHD in offspring could not be demonstrated. Secondly, due to the limited number of areas sampled (only 30 counties and 175 townships), the variation trends in relevant factors were difficult to isolate. www.nature.com/scientificreports/ Thirdly, the reported information on maternal lifestyle during pregnancy such as that related to alcohol intake, passive smoking, and consumption of tea or coffee was self-reported by mothers. The duration of time that had elapsed between pregnancy and the maternal self-report may have contributed to recall bias. Therefore, we consciously adopted systematic approaches to maximize power and minimize bias in the study. For example, all questionnaires and procedures were pre-tested, and detailed interviewer guides were developed before the formal survey. During maternal interviews, the standardized questionnaire and detailed supporting material (e.g., calendars) were provided to participants to ensure that complete and accurate responses was obtained.
Reporting accuracy was ensured for mothers through highly structured interviews by skilled field staff from Xi'an Jiaotong University Health Science Center and experienced doctors from local maternal and children's hospitals who were blinded to the study hypotheses. Fourthly, the study subject to other potential confounders such as maternal illness, drug use, genetics, maternal obesity, and environmental factors. Despite these limitations, this study was the first exploration of the geographic varying association of maternal lifestyle factors with CHD in northwestern China, and filled in some of the gaps.
In conclusion, the prevalence of CHD among live births and the rates of the maternal consumption of alcohol, tea and coffee were higher in South Shaanxi, while the proportion of mothers with a history of passive smoking in pregnancy was highest in North Shaanxi. Maternal alcohol consumption and passive smoking during pregnancy may be important risk factors for CHD, especially in North Shaanxi. Tea consumption during pregnancy appeared to increase risk for CHD in North and Central Shaanxi.

Materials and methods
Study design and participants. From August to November 2013, a large-scale epidemiologic survey of BD among infants born during 2010-2013, was conducted in Shaanxi province in northwestern China. Because of the imbalance in population density between rural and urban areas, a stratified multi-stage sampling method was used to determine sampling units. While China's rural areas, counties, townships, and villages have a threelevel administrative structure, the administrative structure in urban areas consists of districts, streets, and communities. Twenty counties and ten districts were randomly selected for inclusion in the study. Six townships in each county were chosen, and 6 villages in each selected township were randomly selected. In urban areas, 3 streets in the chosen district and 6 communities on the chosen street were selected randomly. In the end, a completely random sampling method was adopted to select 30 live births and their mothers for each selected village and 60 live births and their mothers for each selected community. The response rate was 92.7%, and a total of 29,098 live babies were ultimately enrolled, representing 9.00% of the province's total infant population 18 .

Data collection.
Prior to the formal survey, all questionnaires and procedures were pre-tested, and detailed interviewer guides were developed. The final questionnaire had four parts: birth defects questionnaire, family questionnaire, reproductive history questionnaire, and family history questionnaire. In the survey, all the major data collected was reported by mothers with live births. After written informed consent was obtained, faceto-face interviews were performed to extract data on sociodemographic characteristics and maternal lifestyle during pregnancy 18 .
During the investigation, ten survey teams were established for field data collection in the areas surveyed, with each group consisting of 10-12 investigators, a supervisor from Xi'an Jiaotong University Health Science Center, and a pediatrician from the local maternal and child health hospital. Before the start of the survey, all investigators were trained via lectures and field exercises to standardize implementation of the questionnaire. The supervisors were required to report any questionnaires with errors and/or incomplete information. The doctors in each team helped to collect information on BD. The study was strongly supported by local hospitals, administrative health services, and the Shaanxi Provincial Ministry of Health. All data collection occurred at local village clinics and community health service centers. The study was reviewed and approved by the Ethics Committee of Xi'an Jiaotong University Health Science Center, and all methods involving humans were carried out in accordance with relevant guidelines (Declaration of Helsinki).
The diagnosis of CHD was a multistage process. As a first step, infants were screened for CHD by a doctor from the local maternal and child health hospital. Isolated patent foramen ovale and patent ductus arteriosus among neonates < 28 days of age were considered to be normal neonatal findings and excluded from the study. As soon as infants were found to have CHD, their mothers were interviewed. The attending physician reviewed the medical record and obtained information pertaining to diagnosis time, hospital, and type of CHD. For new cases of suspected CHD, a special neonatal echocardiogram and electrocardiogram were required for a definitive diagnosis at the First Affiliated Hospital of Xi'an Jiaotong University 18 .
In this study, the vector map of Shaanxi province was obtained through the "Cropping" function of ArcMap, from the Chinese vector map. Meanwhile, geographic information for each township/street in Shaanxi province was acquired from the Geographic Information Comprehensive Service Website (https ://www.tiand itu.com/). The information obtained included geographic coordinates, boundaries, and altitude. The basic geographic units in this study were townships and streets. Geographic location was defined as the geographic coordinates (i.e., latitude/longitude) for the township/street center where the participant lived.
Statistical analysis. The analysis procedure is shown in detail below. The geographic distribution of maternal lifestyle during pregnancy were shown on a map using ArcGIS 10.0 software (Environmental Systems Research Institute, Inc., Redlands, CA, US). Non-spatial logistic regression analysis was performed to explore the association between maternal lifestyle during pregnancy and CHD controlling for other confounders by means of IBM SPSS Statistics 25. Then, all explanatory variables were entered into the GWLR model. To confirm geographic discrepancy in the relationship between maternal lifestyle during pregnancy and CHD after adjusting for other factors, local average estimates (OR) and P values for each individual were displayed on a map using ArcGIS 10.0 software. The GWLR formula is as follows 33 : The equation assumes y i is any CHD for each individual i, x ij is a set of independent variables (j = 1, …, k) for individual i, (u i , v i ) is the x-y coordinates of individual i; β ij is the estimated effect of independent variable j for individual i.
The GWLR model was estimated with the iterative reweighted least squares method using GWR 4.0 software (https ://geoda cente r.asu.edu/softw are/downl oads/gwr_downl oads). For modeling, a GWLR equation was estimated for each township/street based on the observations for nearby townships/streets. A distance-based weighting scheme was used to allocate weights to each township/street. The kernel type and function for geographic weighting to estimate local coefficients for each township/street and bandwidth size was adaptive bisquare. The Corrected Akaike Information Criterion (AICc) was used for the golden section search method to select bandwidth in adaptive kernels. The variability of each coefficient in geographic space was tested by comparing different 34 . In the GWLR model, exponentiation of the explanatory variable coefficients was ultimately calculated to acquire the OR corresponding to the unit change in the variable 35 .