(C) PLOS One [1]. This unaltered content originally appeared in journals.plosone.org. Licensed under Creative Commons Attribution (CC BY) license. url:https://journals.plos.org/plosone/s/licenses-and-copyright ------------ Risks of specific congenital anomalies in offspring of women with diabetes: A systematic review and meta-analysis of population-based studies including over 80 million births ['Tie-Ning Zhang', 'Department Of Clinical Epidemiology', 'Shengjing Hospital Of China Medical University', 'Shenyang', 'Clinical Research Center', 'Department Of Pediatrics', 'Xin-Mei Huang', 'Department Of Endocrinology', 'Shanghai Fifth People S Hospital', 'Fudan University'] Date: 2022-02 In this study, we observed an increased rate of CAs in offspring of women with diabetes and noted the differences for PGDM versus GDM. The RRs of overall CAs and CHDs in offspring of women with PGDM were higher than those in offspring of women with GDM. Screening for diabetes in pregnant women may enable better glycemic control, and may enable identification of offspring at risk for CAs. We searched PubMed and Embase from database inception to 15 October 2021 for population-based studies reporting on type-specific CAs in offspring born to women with PGDM (combined type 1 and 2) or GDM, with no limitation on language. Reviewers extracted data for relevant outcomes and performed random effects meta-analyses, subgroup analyses, and multivariable meta-regression. Risk of bias appraisal was performed using the Cochrane Risk of Bias Tool. This study was registered in PROSPERO (CRD42021229217). Primary outcomes were overall CAs and CHDs. Secondary outcomes were type-specific CAs. Overall, 59 population-based studies published from 1990 to 2021 with 80,437,056 participants met the inclusion criteria. Of the participants, 2,407,862 (3.0%) women had PGDM and 2,353,205 (2.9%) women had GDM. The meta-analyses showed increased risks of overall CAs/CHDs in offspring born to women with PGDM (for overall CAs, relative risk [RR] = 1.99, 95% CI 1.82 to 2.17, P < 0.001; for CHDs, RR = 3.46, 95% CI 2.77 to 4.32, P < 0.001) or GDM (for overall CAs, RR = 1.18, 95% CI 1.13 to 1.23, P < 0.001; for CHDs, RR = 1.50, 95% CI 1.38 to 1.64, P < 0.001). The results of the meta-regression analyses showed significant differences in RRs of CAs/CHDs in PGDM versus GDM (all P < 0.001). Of the 23 CA categories, excluding CHD-related categories, in offspring, maternal PGDM was associated with a significantly increased risk of CAs in 21 categories; the corresponding RRs ranged from 1.57 (for hypospadias, 95% CI 1.22 to 2.02) to 18.18 (for holoprosencephaly, 95% CI 4.03 to 82.06). Maternal GDM was associated with a small but significant increase in the risk of CAs in 9 categories; the corresponding RRs ranged from 1.14 (for limb reduction, 95% CI 1.06 to 1.23) to 5.70 (for heterotaxia, 95% CI 1.09 to 29.92). The main limitation of our analysis is that some high significant heterogeneity still persisted in both subgroup and sensitivity analyses. Pre-gestational diabetes mellitus (PGDM) has been known to be a risk factor for congenital heart defects (CHDs) for decades. However, the associations between maternal PGDM and gestational diabetes mellitus (GDM) and the risk of specific types of CHDs and congenital anomalies (CAs) in other systems remain under debate. We aimed to investigate type-specific CAs in offspring of women with diabetes and to examine the extent to which types of maternal diabetes are associated with increased risk of CAs in offspring. Previous meta-analyses have mainly focused on the associations between maternal diabetes and CHDs in offspring, and little is known about the influence of maternal diabetes on other specific types of CAs [ 13 , 14 ]. Additionally, new data from population-based studies of more than 36 million births have provided solid estimates of the risk of CHDs in offspring of women with diabetes [ 4 , 10 – 12 ]. This considerable amount of data could also be used to explore the association between maternal diabetes and other types of CAs. Currently, a quantitative summary of population-based studies on the associations between maternal diabetes (pre-gestational or gestational) and type-specific CAs in offspring is lacking. Comprehensive estimates of the risk of specific CAs associated with maternal diabetes are needed to counsel patients and for public health purposes. Moreover, it is essential that estimates are provided according to different types of maternal diabetes, given the diversity in etiology, treatment, and prognosis. Gestational diabetes mellitus (GDM), which is defined as any degree of glucose intolerance with onset or first recognition during pregnancy, is one of the most common complications of pregnancy and affects up to 9%–26% of the obstetric population [ 8 , 9 ]. Similar to PGDM, GDM also has a considerable impact on the health outcomes of the mother and infant during pregnancy, delivery, and beyond. Recently, an increasing number of studies have concentrated on evaluating the risks of specific types of CAs in offspring born to women with GDM [ 4 , 5 , 10 – 12 ]. The early period of organogenesis, which occurs during the third to eighth week of gestation, is an important time for organ development. However, hyperglycemia associated with GDM occurs after this critical early window for organogenesis. Therefore, the question as to whether there is an association between GDM and the risk of specific types of CAs in offspring remains. Currently, the global prevalence of diabetes is increasing among women of reproductive age [ 1 , 2 ]. A diabetic intrauterine environment can cause placental dysfunction and hormone alterations, which could lead to various congenital anomalies (CAs) in offspring of women with diabetes [ 1 ]. Notably, pre-gestational diabetes mellitus (PGDM, which includes type 1 and 2 diabetes) has been known to be a risk factor for congenital heart defects (CHDs) for decades [ 3 ]. However, there is controversy among current research regarding the association between maternal PGDM and the risk of specific types of CHDs and other CAs of the nervous, digestive, genitourinary, and musculoskeletal systems [ 4 – 7 ]. Further studies are thus needed for clarification of this risk. For studies that reported effect sizes separately, the results were pooled using a fixed effects model to obtain an overall estimate and then included in the pooled effect size in the meta-analysis. The effective count method proposed by Hamling et al. [ 20 ] was used to recalculate the effect sizes. If a selected study did not include an effect size, the unadjusted risk estimate and 95% CI were calculated from the raw data for simplicity using EpiCalc 2000 ( https://en.freedownloadmanager.org/Windows-PC/EpiCalc-2000-FREE.html ). Estimates were pooled using the DerSimonian and Laird random effects model to calculate summarized RRs and 95% CIs [ 21 ], in which I 2 values were calculated as indicators of heterogeneity. I 2 values of ≤25%, 26%–50%, 51%–74%, and ≥75% were considered to indicate no, low, moderate, and high heterogeneity between the included studies, respectively [ 22 ]. For the primary outcomes of the study, subgroup analyses were undertaken to explore causes of heterogeneity: by region (Europe, North/South America, or Asia-Pacific), year of enrollment (categorized using the median as the cutoff value: before 1997 or in or after 1997), number of participants (categorized using the median as the cutoff value: <282,260 or ≥282,260), and adjustment for confounders (i.e., maternal age, race/ethnicity, body mass index, education, smoking/alcohol consumption, parity, and pregnancy complications). Heterogeneity between subgroups was evaluated by meta-regression analysis if data were reported in more than 10 studies, following to the Cochrane guidelines [ 23 ]. Meta-regression analyses were also used to examine the extent to which the types of maternal diabetes (i.e., pre-gestational or gestational) are associated with increased risk of overall CAs/CHDs in offspring. Publication bias was examined by inspecting funnel plots for the outcomes and was further evaluated with Begg’s test [ 24 ] and Egger’s test [ 25 ] if sufficient studies existed (n ≥ 10) [ 17 ]. A sensitivity analysis was undertaken to explore the association of each study with the overall pooled estimate. Statistical analyses were conducted using Stata version 13.0 (StataCorp, College Station, Texas). A 2-tailed P value less than 0.05 was considered statistically significant. The risk of bias assessment was conducted by T-NZ and S-YG using the Risk of Bias in Non-randomized Studies–of Interventions (ROBINS-I) tool [ 19 ]. This tool comprises 7 domains: bias due to confounding, bias in the selection of participants, bias in the classification of interventions, bias due to deviations from intended interventions, bias due to missing data, bias in measurement of outcomes, and bias in the selection of the reported result. We rated the possible risk of bias in each of the 7 domains as low risk, moderate risk, serious risk, critical risk, or no information for each available outcome of each included study. Two authors (T-NZ and S-YG) independently performed data extraction according to the Cochrane Handbook guidelines [ 17 ]. Findings were reported according to PRISMA [ 15 ] and MOOSE guidance [ 16 ]. Any disagreement was settled by consensus among all authors. For studies that did not report any adjusted effect sizes, the crude risk estimate was used. If an included study reported several risk estimates, we extracted the fully adjusted effect sizes. For studies that reported the risk estimates of CAs stratified by isolated and multiple statuses, we used the effect sizes of the isolated CAs. Because odds ratios, prevalence rate ratios, and hazard ratios are excellent approximations of risk ratios in the case of rare outcomes [ 18 ], all risk estimates are referred to and reported as relative risks (RRs) for simplicity. If an included study lacked required data, we asked for help from the Library of Shengjing Hospital of China Medical University to get the missing information. Two independent authors (S-YG and T-NZ) reviewed the titles and abstracts to identify any relevant studies. The full texts of potentially eligible studies that appeared to meet the inclusion criteria were then obtained and independently evaluated by the 2 reviewers. Any disagreement was settled by consensus among all authors. If multiple studies were derived from the same dataset and reported the same associated outcome, the study with the most complete findings or the greatest number of participants was included for analysis. The literature review and study selection process referenced the PRISMA flowchart ( Fig 1 ). When information needed for inclusion in the analyses was missing, the Library of Shengjing Hospital of China Medical University helped us get full access of the relevant data. Population-based cross-sectional, case–control, and cohort studies that reported original data were eligible for inclusion if they (1) reported any CAs in offspring born to women with diabetes (i.e., pre-gestational [combined type 1 and 2] or gestational diabetes), (2) had a comparison group that included mothers without diabetes, and (3) provided sufficient data from which a risk estimate could be calculated if a risk estimate was not reported. All conference abstracts, guidelines, case reports, case series, commentaries, letters, and animal studies were excluded. The sensitivity analysis evaluated the effect of omitting 1 study at a time from each analysis. In the sensitivity analysis, we observed that the high I 2 value of 81.2% shown in the results for CHDs in offspring of mothers with GDM reduced to a moderate I 2 value of 55.0% when excluding the study by Billionnet et al. [ 47 ]. Although the increased risk association remained robust across scenarios, some moderate to significant heterogeneity still persisted and could not be reduced in sensitivity analyses. To explore the source of heterogeneity, we performed subgroup and meta-regression analyses in the predefined subgroups of study location, year of enrollment, study sample size, and adjustment for confounders (Tables 4 and 5 ). The findings of increased overall CA/CHD risk associated with maternal diabetes were consistently observed in most of the subgroup analyses. The results of the subgroup analyses suggested that differences in study sample size, population region, year of enrollment, and adjustment for confounders were major sources of heterogeneity. We observed that the high I 2 value of 81.2% observed in the results for CHDs in offspring of mothers with GDM was reduced to no (I 2 = 0%) or low (I 2 = 47.7%) heterogeneity after adjustment for race/ethnicity, body mass index, education, smoking/alcohol consumption, parity, and pregnancy complications ( Table 5 ). In addition, the results of meta-regression analyses showed statistically significant differences in the RRs of CAs/CHDs in PGDM versus GDM (all P meta-regression < 0.001) ( Fig 6 ). Maternal GDM was associated with an increased risk of CAs of the eye, ear, face, and neck (RR = 1.15, 95% CI 1.09 to 1.22, I 2 = 0%, P = 0.355) and musculoskeletal system (RR = 1.18, 95% CI 1.15 to 1.22, I 2 = 0%, P = 0.424) in offspring. In addition, maternal GDM also contributed to an increased risk of specific types of CAs in offspring, including hydrocephaly (RR = 1.34, 95% CI 1.16 to 1.54, I 2 = 0%, P = 0.960), holoprosencephaly (RR = 1.87, 95% CI 1.09 to 3.22, I 2 = 0%, P = 0.558), cleft lip with or without cleft palate (RR = 1.26, 95% CI 1.19 to 1.34, I 2 = 0%, P = 0.547), diaphragmatic hernia (RR = 1.21, 95% CI 1.08 to 1.37, I 2 = 0%, P = 0.779), omphalocele (RR = 1.21, 95% CI 1.05 to 1.40, I 2 = 0%, P = 0.743), and hypospadias (RR = 1.29, 95% CI 1.16 to 1.44, I 2 = 45.9%, P = 0.100) ( Table 3 ; Figs J1–J22 in S1 Fig ). We examined the associations between maternal PGDM and other type-specific CAs in offspring. Our results suggested that offspring of women with PGDM had an increased risk of CAs of the nervous system (RR = 2.54, 95% CI 1.73 to 3.73, I 2 = 94.8%, P < 0.001); eye, ear, face, and neck (RR = 3.14, 95% CI 2.90 to 3.39, I 2 = 0%, P = 0.444); digestive system (RR = 2.02, 95% CI 1.24 to 3.28, I 2 = 92.3%, P < 0.001); genitourinary system (RR = 1.73, 95% CI 1.35 to 2.21, I 2 = 89.2%, P < 0.001); and musculoskeletal system (RR = 1.98, 95% CI 1.45 to 2.72, I 2 = 94.4%, P < 0.001), as well as an increased risk of multiple CAs (RR = 3.06, 95% CI 2.36 to 3.96, I 2 = 39.6%, P = 0.158). The associations were statistically significant in 14 of 16 type-specific CA categories. The corresponding RRs ranged from 1.57 (for hypospadias, 95% CI 1.22 to 2.02, I 2 = 74.1%, P < 0.001) to 18.18 (for holoprosencephaly, 95% CI 4.03 to 82.06, I 2 = 66.3%, P = 0.085) ( Table 3 ; Figs H1–H25 in S1 Fig ). Eleven studies explored the relationship between GDM and CHDs in offspring [ 4 , 6 , 11 , 12 , 38 , 44 , 47 , 51 , 52 , 56 , 72 ]. Our results suggested that maternal GDM is associated with CHDs (RR = 1.50, 95% CI 1.38 to 1.64, I 2 = 81.2%, P < 0.001; Fig 5 ), with no evidence of publication bias (Begg’s P = 0.837, Egger’s P = 0.885; Fig F in S1 Fig ). Regarding specific types of CHDs, we found that offspring of women with GDM had an increased risk of heterotaxia (RR = 5.70, 95% CI 1.09 to 29.92, I 2 = 85.7%, P = 0.008), tetralogy of Fallot (RR = 1.41, 95% CI 1.20 to 1.66, I 2 = 0%, P = 0.600), left ventricular outflow tract defect (RR = 1.67, 95% CI 1.15 to 2.41, I 2 = 50.0%, P = 0.112), coarctation of aorta (RR = 1.50, 95% CI 1.23 to 1.83, I 2 = 35.4%, P = 0.213), right ventricular outflow tract defect (RR = 1.25, 95% CI 1.03 to 1.53, I 2 = 0%, P = 0.739), ventricular septal defect (RR = 1.31, 95% CI 1.24 to 1.38, I 2 = 0%, P = 0.960), and atrial septal defect (RR = 1.45, 95% CI 1.40 to 1.50, I 2 = 0%, P = 0.426) ( Table 2 ; Figs I1–I15 in S1 Fig ). A total of 18 studies reported on the association between maternal PGDM and CHDs in offspring [ 4 , 6 , 10 – 12 , 28 – 30 , 38 , 44 , 47 , 49 , 52 , 54 , 56 , 63 , 68 , 74 ]. Our results suggested that there is a statistically significant increase in risk of CHDs in offspring of women with PGDM (RR = 3.46, 95% CI 2.77 to 4.32, I 2 = 98.2%, P < 0.001; Fig 4 ), with no evidence of publication bias (Begg’s P = 0.60, Egger’s P = 0.85; Fig E in S1 Fig ). Similarly, maternal type 1 and type 2 diabetes were associated with increased risk of CHDs in offspring (type 1: RR = 3.75, 95% CI 1.86 to 7.57, I 2 = 99.1%, P < 0.001; Fig K2 in S1 Fig ; type 2: RR = 3.15, 95% CI 1.72 to 5.78, I 2 = 93.6%, P < 0.001; Fig J2 in S1 Fig ). Notably, we found that maternal PGDM was associated with increased risk of all specific types of CHDs available for examination in the present study. The RRs of specific types of CHDs ranged from 2.23 (for hypoplastic left heart, 95% CI 1.07 to 4.64, I 2 = 64.0%, P = 0.040) to 12.16 (for truncus arteriosus, 95% CI 7.52 to 19.68, I 2 = 0%, P = 0.866) ( Table 2 ; Figs G1–G18 in S1 Fig ). We assessed the risk of bias for 34 of 59 included studies using ROBINS-I. The assessments are summarized for primary outcomes in Figs A–D in S1 Fig . None of the included studies were rated with a low risk of bias in all domains. The main causes of serious or critical bias risk according to ROBINS-I were weaknesses in the confounding bias domain, selection of participant bias domain, and missing data bias domain. Discussion To the best of our knowledge, the present study is the first comprehensive systematic review and meta-analysis of population-based studies of over 80 million participants that shows an increased risk of type-specific CAs, especially CHDs, in offspring of women with pre-gestational or gestational diabetes. The study findings suggested that maternal PGDM was associated with a significant increase in the risk of CAs in offspring in 38 of 45 categories, while maternal GDM was associated with a small but significant increase in the risk of CAs in 16 of 45 categories. The corresponding RRs of overall CAs/CHDs in offspring of women with PGDM were higher than those in offspring of women with GDM, with no overlap in the 95% CIs. Although the exact pathophysiology of the relationship between maternal diabetes and CAs in offspring remains unclear, metabolic changes in women with diabetes could play a critical role in the development of CAs in their offspring. Sustained hyperglycemia is the main characteristic of diabetes; this activates multiple metabolic pathways that play a role in the formation of CAs [77]. Notably, a common mechanism behind diabetic complications is mitochondrial superoxide radical production [77]. The production of reactive oxygen species (ROS) is induced by hyperglycemia, which is the crucial process in diabetes mellitus pathogenesis, and oxidative stress (OS) is known to affect embryonic development [77–79]. Although OS does not cause a direct genotoxic effect, a previous study showed that OS affects the expression of some genes involved in the various stages of embryonic development, and each gene affected might have a specific sensitivity to hyperglycemic conditions and changes in the cellular redox state, thus mediating the formation of CAs in offspring [77]. Studies have suggested that the formation of CAs involves AMP-activated protein kinase (AMPK), an enzyme with kinase activity that is activated in response to an increase in adenosine monophosphate nucleoside levels [78,80]. AMPK could play a key role in the formation of CAs; it participates in the regulation of energy metabolism and, once activated, moves into the cell nucleus and phosphorylates multiple proteins, including hypoxia-inducible factor 1α, which could mediate the development of CAs [77]. However, further study is needed to discern the exact mechanisms involved in the activation of AMPK and the induction of CAs by ROS. An important finding from our meta-regression analyses is the statistically significant difference in the risk of overall CAs in offspring of women with PGDM versus offspring of women with GDM. That is, the risk of CAs in offspring was higher in women with PGDM than in those with GDM. Pregnancy begins at fertilization, and organogenesis begins during the third to eighth week post-conception and continues until birth. Therefore, the first trimester of pregnancy is the most critical period for organogenesis. In women with PGDM, there can be a lengthy period of sustained hyperglycemia before and during pregnancy, which can significantly impact organogenesis and contribute to CAs in offspring. This differs from GDM, which is usually diagnosed during the 24th to 28th week of gestation [8]. Therefore, in a woman with GDM, blood glucose levels could be normal or just slightly elevated during the first trimester, leading to minimal influence on organogenesis. This could partly explain why offspring of women with PGDM were at greater risk for CAs compared with offspring of women with GDM. However, women who develop GDM during pregnancy usually have evidence of metabolic dysfunction before pregnancy, such as pancreatic β-cell defects and increased insulin resistance [81,82], which may contribute to the development of hyperglycemia and thus increase the rate of malformations in infants, although further studies are needed to elucidate the potential mechanisms involved. Another key finding from our meta-regression analysis was the similar result observed regarding CHDs in offspring of women with PGDM and offspring of women with GDM. The heart is the first functional organ to develop and starts to beat and pump blood at around 22 days after fertilization [83]. The septum, including the interatrial septum, begins to form at 4 to 7 weeks of gestation. Hyperglycemia could have a more critical influence on heart development in the early stage of pregnancy than in the late stage of pregnancy. Therefore, screening for diabetes in pregnant women will enable better glycemic control, which might reduce the rate of malformations, especially during organogenesis. However, the exact mechanisms underlying the influence of diabetes on organogenesis in different stages remain unclear and require further study. In the present study, the results regarding specific types of CHDs in offspring with maternal PGDM were consistent with 2 previous meta-analyses [13,14]. A recent systematic review and meta-analysis conducted by Chen and colleagues involved a pooled analysis of 24 population-based studies and 18 hospital-based studies; the findings suggested that maternal GDM was significantly associated with the risk of most phenotypes of CHDs [13]. New data from population-based studies of more than 36 million births provided solid estimates of the associations between maternal GDM and specific types of CHDs in offspring [4,10–12]. However, these studies mainly focused on the association between different types of maternal diabetes and CHDs. Little is known about the association between maternal diabetes and other specific types of CAs in offspring or the extent to which types of maternal diabetes are associated with the increased risk of CAs. One recent meta-analysis by Parimi and colleagues explored the association between maternal diabetes and the risk in offspring of CAKUT, which refers to a range of structural and functional anomalies of the kidney, collecting system, bladder, and urethra [84]. Our findings were in line with the results from Parimi et al. [84] that showed that offspring of mothers with PGDM had an almost 2-fold increased risk of CAKUT; however, results regarding the association between maternal GDM and the risk of CAKUT were inconsistent. Our findings demonstrated associations between maternal diabetes and 23 CA categories (excluding CHD-related categories) in offspring and suggested that offspring of women with PGDM had an increased risk of 21 specific types of CAs, while increased risks of 9 specific types of CAs were observed in offspring of women with GDM. Strengths and limitations Our study has several strengths. The first strength is the large sample size of over 80 million births from population-based data, which provides robust evidence regarding the risk of CAs in offspring of women with diabetes and are widely generalizable. Second, our study examined the associations between maternal diabetes and various types of CAs across multiple categories of maternal diabetes. Unlike previous studies [13,14,84] that only assessed the risk of CHDs or CAKUT in maternal diabetes, the present study systematically and quantitatively summarized the associations between maternal diabetes and 45 type-specific CAs in offspring. Third, consistent results of the pooled RRs supported the robustness of the findings of our study. Finally, the current study examined the extent to which types of maternal diabetes (i.e., pre-gestational and gestational) are associated with increased risk of CAs in offspring. The relative consistency of associations observed appears to support the hypothesis that maternal diabetes, especially PGDM, increases the likelihood of type-specific CAs in offspring. However, several limitations should be noted. Although the increased risk association remained robust across various scenarios, some high levels of statistical heterogeneity generally persisted and could not be reduced in subgroup and sensitivity analyses. There were some causes of heterogeneity in the included studies. First, there is lack of consensus and uniformity in the screening standards and diagnostic criteria for GDM. Also, pre-pregnancy diabetes is sometimes unrecognized and discovered only during pregnancy as GDM, which could lead to overestimation of RRs associated with GDM. Second, the ascertainment of some CAs may vary substantially between studies. Some CAs are easy to ascertain (e.g., anencephaly), while some may not be recognized immediately after birth and may be discovered only later in infancy (e.g., milder atrial septal defects). This also contributes to the heterogeneity of the results. Third, most studies included live births only; the lack of information on stillbirths and terminations of pregnancy for fetal anomaly could introduce selection bias and lead to underestimation of the strength of the associations between maternal diabetes and risk of CAs in offspring. Fourth, there may be other unmeasured confounding factors in addition to those adjusted for in each study. In this regard, further study could be performed to reduce the aforementioned causes of heterogeneity in a more in-depth analysis. An additional limitation was that although we summarized and quantified the existing population-based data on overall CAs/CHD observed under maternal type 1 or type 2 diabetes, data on other type-specific CAs in offspring associated with maternal type 1 or type 2 diabetes are limited. Additional studies are needed to address this issue. Furthermore, information on treatment (e.g., insulin use) or how well-controlled blood glucose levels were in the study participants was not available in most of the studies included in the current study. Further work should strive to address this lack of information. Finally, we observed a negative association between maternal GDM and risk of gastroschisis. The reasons for why maternal GDM was inversely associated with the risk of gastroschisis are currently unknown; this finding warrants confirmation and further investigation in future studies. Residual confounding may contribute to the inverse association, but further confirmation is still needed. [END] [1] Url: https://journals.plos.org/plosmedicine/article?id=10.1371/journal.pmed.1003900 (C) Plos One. "Accelerating the publication of peer-reviewed science." Licensed under Creative Commons Attribution (CC BY 4.0) URL: https://creativecommons.org/licenses/by/4.0/ via Magical.Fish Gopher News Feeds: gopher://magical.fish/1/feeds/news/plosone/