- 1Lozano Blesa University Clinical Hospital, Zaragoza, Spain
- 2Group in Health Services in Aragón (GRISSA), Fundación Instituto de Investigación Sanitaria de Aragón (IIS Aragón), Zaragoza, Spain
- 3Department of Preventive Medicine and Public Health, University of Zaragoza, Zaragoza, Spain
- 4Network for Research on Chronicity, Primary Care, and Health Promotion (RICAPPS), Instituto de Salud Carlos III (ISCIII), Madrid, Spain
- 5Aragon Health Service, Zaragoza, Spain
Aim: To describe the patterns of first prescription of antidiabetic drugs (AD) in patients with type 2 diabetes (T2D) and analyze the factors associated with the prescription of a modern one.
Methods: Observational longitudinal study conducted in the CArdiovascular Risk factors for HEalth Services research (CARhES) cohort. Individuals older than 15, resident in Aragón (Spain), diagnosed with T2D during 2018–2022 were selected and followed-up until 31st December 2022. Secondary use of data from the health system provided sociodemographic, clinical and pharmacological prescription information. We also considered additional variables by Basic Healthcare Area (BHA) of residence. AD were classified into “classical” and “modern” and their differences were described and compared. A multilevel methodology stratified by sex was developed, considering individual characteristics and characteristics of the BHA of residence, to analyze the factors associated to a modern AD.
Results: Our population-based cohort of 22,892 patients were mostly male, native, low-income and living in non-depopulated BHA. People who were younger, with heart failure, ischemic heart disease, chronic renal failure, obesity, with a previous major adverse cardiovascular event, higher socioeconomic level or lived in less deprived and more depopulated areas were more likely to get a modern AD prescription.
Conclusion: Our analyses showed that prescribing practices vary according to a range of sociodemographic, clinical and geographical characteristics. Knowledge of these factors is essential for implementing and improving equitable and person-centered approaches.
1 Introduction
The prevalence of diabetes is increasing at an unstoppable rate (Internacional Diabetes Federation, 2023). Type 2 diabetes (T2D) accounts for the vast majority (over 90%) of diabetes worldwide. Upon diagnosis of T2D, it is common practice to recommend hypoglycemic treatment in conjunction with lifestyle modification (Alemán et al., 2018). The timeliness of this treatment should be assessed by the physician based on the patient’s comorbidities and analytical results (mainly fasting plasma glucose (FPG), glycosylated hemoglobin (HbA1c) and glomerular filtration rate (GFR)) in order to select the most appropriate treatment for each patient.
The increasing availability of antidiabetic drugs (AD) for the management of T2D currently allows different options for achieving glycemic targets and choosing the therapy that best suits the individual characteristics of each patient (Dahlén et al., 2022; DeMarsilis et al., 2022). AD can be classified as “modern” or “classical” based on their mechanism of action. “Classical” correspond to those whose mechanism is based on β-cell dysfunction and insulin resistance (metformin, sulfonylureas, alpha-glucosidase inhibitors, thiazolidinediones, repaglinide, insulins) and “modern” those with an incretin effect and with a hypoglycaemic action at renal level (Dipeptidyl peptidase four inhibitors (DPP4i), Sodium-Glucose Transport Protein 2 (SGLT2) inhibitors (SGLT2i) and glucagon-like peptide-1 agonists (GLP1a)). After the commercialization of new drugs, especially SGLT2i and GLP1a, research analyzing their effectiveness and safety has demonstrated the increased benefits obtained versus the classical therapy, especially in individuals with certain comorbidities such as cardiovascular (CV) disease, chronic kidney disease or heart failure (Herrington et al., 2023; van der Aart-van der Beek et al., 2022; Palmer et al., 2021; Ussher and Drucker, 2023). These drugs have been shown to reduce CV mortality, slow the progression of renal impairment, and reduce the frequency of hospital readmissions for heart failure (McDonagh et al., 2022). Because of this, they have been gaining prominence as first-line drugs in certain pathologies and clinical conditions, which has been reflected in the evolution of treatment algorithms (Alemán et al., 2014; Mata et al., 2020; García Soidán, 2023).
It is necessary to analyze whether, following the introduction of these new drugs with new indications, they have become a reference in routine clinical practice. Furthermore, it is important to identify the factors associated with their prescription. In addition to the clinical characteristics of the patients and their analytical values, factors such as their socioeconomic level or geographical location must also be considered, in order to explore the existence of variations in clinical practice.
The objective of this study is to describe the pattern of AD prescription in newly diagnosed patients with T2D and to analyze the factors associated with the prescription of a modern AD.
2 Materials and methods
2.1 Study design and population
We conducted a retrospective, longitudinal, observational study based on the CARhES Cohort (CArdiovascular Risk factors for HEalth Services research), a dynamic cohort that collects information on patients with at least one CV risk factor (hypertension, T2D or dyslipidemia) in Aragón and was launched from 2017. The CARhES cohort was approved by the Research Ethics Committee of Aragón (CEICA, PI21/148) (Aguilar-Palacio et al., 2024). Aragon is located in the north-east of Spain and had around 1.3 million inhabitants in 2018. The Spanish health system is mainly tax financed and is based on universality, free healthcare, equity and fairness in funding (Bernal-Delgado et al., 2024).
For the present study, all the individuals of the CARhES cohort aged ≥16 years old with a new diagnosis of T2D between 1 January 2018 and 31 December 2022 were selected. The diabetes diagnosis was identified by a new diagnosis in Primary Care (PC), an emergency department or hospital, or with the first prescription of an AD. All patients with a first AD prescription prior to the diagnosis of T2D and those without treatment were excluded. The flow chart of the study population is shown in Figure 1.

Figure 1. Flow chart of the study population. n: number of patients; T2D: type 2 diabetes; AD: antidiabetic drug.
2.2 Variables of the study and data sources
Sociodemographic and clinical characteristics of all the individuals in the cohort were described. Regarding sociodemographic characteristics, we considered sex, age (categorized into 15–44, 45–64, 65–79 and 80 years or older), immigrant status (native/immigrant) and socioeconomic level. Socioeconomic level was calculated on the basis of pharmacy copayment levels and the type of economic activity, according to the type of user of the Aragón health service. From the combination of these two variables, six mutually exclusive categories were obtained (employed individuals earning less than 18,000 € per year, employed individuals earning 18,000 € per year or more, individuals receiving the unemployment allowance, individuals with a contributory pension of less than 18,000 € per year and individuals receiving free medicines (people with minimum integration income or who no longer receive the unemployment allowance), individuals with a contributory pension of 18,000 € per year or more, and others). The clinical information included was obtained from the adjusted morbidity groups (GMA) (Barrio-Cortes et al., 2020). This source of information considers all medical diagnoses available in primary healthcare, emergencies department and hospital discharge records. The morbidity burden is a summary index included in the GMA database. It is based on the clinical conditions present in the patient and the weight assigned to each condition according to the care and resources needed for its management. It was originally developed with data from the Spanish health system. From electronic medical records we collected major adverse CV events (MACE; myocardial infarction, stroke and CV death) and laboratory variables (FPG, HbA1c, GFR). Treatment data was obtained by collecting the ATC code(s) (A10) corresponding to the first prescription date as described in Supplementary Appendix S1.
Treatments were classified into “classical” and “modern”. “Classical” correspond to those whose mechanism is based on β-cell dysfunction and insulin resistance (ATC codes: A10BA2 (metformin), A10BB (sulfonylureas), A10BD (0–6) (combinations of oral blood AD), A10BF (alpha-glucosidase inhibitors), A10BG (thiazolidinediones), A10BX (repaglinide), A10A (insulins)) and “modern” those with an incretin effect and with a hypoglycemic action at renal level (ATC codes: A10BD (Palmer et al., 2021; Ussher and Drucker, 2023; McDonagh et al., 2022; Alemán et al., 2014; Mata et al., 2020; García Soidán, 2023; Aguilar-Palacio et al., 2024; Bernal-Delgado et al., 2024; Barrio-Cortes et al., 2020; Compés Dea et al., 2018; Spanish ministry, 2023; Bates et al., 2015; Moreno-Juste et al., 2020; Mata-Cases et al., 2016; Rojo-Martínez et al., 2020; Wang et al., 2019; Ali et al., 2022; Choe et al., 2018) (combinations of oral blood AD), A10BH (DPP4i), A10BJ (GLP1a), A10BK (SGLT2i)).
Sociodemographic and clinical variables were collected at the time of entry into the cohort with a T2D diagnosis. The recording of the variable obesity is not uniform, so that obesity at any time was taken to mean that they were obese. Analytical variables were collected at the closest time prior to the first prescription of AD (from 12 months before or, if this was not possible, up to 1 month after).
We also considered three additional variables related to the Basic Healthcare Area (BHA) of residence of each patient. The first variable was the BHA deprivation index categorized into four quartiles, from least (Q1) to highest (Q4) deprived. This deprivation index combines information of four indicators from the Population and Housing Census 2011 (last index available) (Compés Dea et al., 2018). Other variable obtained by BHA was the classification of the zone according to its depopulation level (assigned based on the criteria of the Spanish Ministry for the Ecological Transition and the Demographic Challenge (Spanish ministry, 2023): non-depopulated, mixed, depopulated). The last variable by BHA was family medicine pressure. This variable is the average number of visits per professional per working day and represents the pressure of care exerted by the general practitioner.
2.3 Statistical analyses
First, the study population was described. Means and standard deviations (SD) were used to describe continuous variables, and frequencies and percentages were used to describe categorical variables. The proportions of classical and modern prescriptions were obtained (%) and described by sex. To know which individual and area characteristics were associated with being firstly prescribed a modern vs. a classical treatment, bivariate analyses were conducted. Statistical differences were assessed using chi-square (categorical variables) and Mann–Whitney tests (continuous variable).
To analyze the factors associated to a modern AD prescription, a multilevel methodology stratified by sex was developed, considering individual characteristics and characteristics of the BHA of residence. We implemented a two-level model, with a cross-classification structure, where individuals could simultaneously belong to more than one group at a given hierarchical level. In this case, we analyzed the type of AD prescription, modern or classical, by their deprivation index (quartiles), depopulation level and family medicine pressure (quartiles), so all were considered random. Cross-random effects were used when each category of one factor co-exists with each category of the other factor (there is at least one category observation for both factors).
The model had the following assumptions: first, the random effects
Interactions between variables were systematically investigated and collinearity was demonstrated. Finally, the likelihood ratio test (LR test) was used to evaluate the final model. The significance of the fixed effects was also evaluated with the Wald Test. All analyses were performed using R statistical software (the R Foundation for Statistical Computing, Vienna, Austria) using a mixed-effects linear regression based on the lme4 package (Bates et al., 2015) in the statistical package R version 4.4.0.
3 Results
The study population was comprised of 22,892 patients with a new T2D diagnosis, 43.52% female, mean age of 64.54 ± 13.41 years (standard deviation) (Table1). The vast majority of the population (84.68%) were native. The 54.25% of the population lived in non-depopulated areas. In terms of socioeconomic level, more than a half of the population had income of less than €18,000 per year, with differences between men and women.
With regard to the presence of comorbidities, 73.89% of patients also presented hyperlipidemia, followed by 55.08% with hypertension. Women had a higher prevalence of heart failure, depression, chronic renal failure, dementia, hypertension and obesity than men. Male exhibited a higher prevalence of chronic obstructive pulmonary disease, ischemic heart disease, hyperlipidemia, and a history of previous MACE.
In consideration of the analytical values obtained prior to the initiation of the prescription, the majority of patients had elevated FPG (mean 154 ± 57.2) and HbA1c (mean 7.4%), in addition to a decreased GFR (mean 80.5 ± 18.5).
Initial treatment with monotherapy was observed in approximately 90% of cases. Metformin was the most commonly dispensed monotherapy drug, accounting for 56.88% of prescriptions, followed by combinations with metformin (12.58%); and SGLT2i (10.68%). Over one-third of patients (36.03%) started treatment with a modern AD, including DPP4i (5.92%), GLP1a (2.21%), or SGLT2i; without differences between women and men.
In Table 2, the differences between patients with a modern AD prescription and a classical one can be observed. The non-native patients and low-income people were less prescribed of modern AD. Patients with heart failure, chronic obstructive pulmonary disease, ischemic heart disease, chronic renal failure and dementia were prescribed more modern AD. Of all patients followed, 746 had had a MACE prior to the diagnosis of T2D (3.26%), of whom 415 were subsequently prescribed a modern AD. When analyzing laboratory data, in patients with poorer FPG, HbA1c and GFR control, modern drugs were prescribed more frequently. Of the prescriptions initiated in PC, the majority were for classical AD, whereas in specialized care, the majority of prescriptions were for a modern one (p = 0.00).

Table 2. Comparison of the characteristics of patients with a first prescription of a modern vs. classical antidiabetic drug. Bivariate analyses.
Different multilevel models were tested for the total population and stratified by sex in order to know which factors were associated with the prescription of a modern AD. The model with a higher explanatory capacity was the one that combined the deprivation index, depopulation level and family medicine pressure as the BHA variables. The results of this model can be found in Table 3. In the adjusted models, we observed that the odds of prescription of modern AD decreased with age, with the group aged from 65–79 years old being the lowest one of using modern AD (odds ratio (OR) 0.48; 95% confidence interval (95%CI) 0.41–0.58). No differences were observed by sex. There were differences by socioeconomic level. So, pensioners presented the lowest probability of modern AD, being statistically significant only in women. In both sexes, patients with heart failure, ischemic heart disease, chronic renal failure, obesity and a previous MACE had more probability of having a prescription of modern AD, in addition to patients with higher FPG, HbA1c and lower GFR.

Table 3. Factors associated with the prescription of a modern antidiabetic drug. Multilevel adjusted analyses for all the population and stratified by sex.
The final section of Table 3 presents the parameters associated with the model. σ2 represents the overall variance of the model, τ00 corresponds to the variances of the random effects, which are the grouping levels, ICC is the intra-class correlation index, indicating how similar the observations are within each class. In other words, it indicates that part of the total variability is due to the variability of the observations in the levels.
The influence of BHA variables can be observed in Figure 2. According to the depopulation level, non-depopulated areas presented the lowest risk of modern AD prescription. So, less deprived areas showed a higher probability of modern AD prescription, that decreases with the increase of deprivation until quartile 3, as is the case with family medicine pressure that increased the risk of modern AD prescription with the increasing pressure until quartile 3. There were not differences by sex as shown in Supplementary Appendix S2.
4 Discussion
In this population-based study conducted in patients with T2D newly diagnosed we observed differences at the individual level and by area of residence in the risk of receiving a modern AD prescription. People who were younger, had more comorbidities, higher socioeconomic level and lived in less deprived areas and depopulated areas were more likely to have a modern AD prescription.
Our study population of patients newly diagnosed with T2D were mostly male (table 1), with a mean age of 64.54, native-born and living in non-depopulated areas and with low incomes. The mean age at which diabetes was diagnosed in the subjects of our study was 64.54 years (SD = 13.41). This finding is consistent with those reported in other studies conducted in Aragón (Moreno-Juste et al., 2020) but higher than that reported in other studies both national (Mata-Cases et al., 2016; Rojo-Martínez et al., 2020) and international (Wang et al., 2019).
Our analyses revealed no significant discrepancies in the proportion of modern AD prescribed between male and female subjects. Nevertheless, the available evidence is inconclusive and exhibits inconsistencies across the studies examined (Ali et al., 2022; Choe et al., 2018) analyzing sex differences in prescribing pattern of AD. Moreover, as in our study, it has been previously described a lower prescription of SGLT2i and GLP1a in older, lower-income and immigrant populations (Ali et al., 2022; Lamprea-Montealegre et al., 2022; Montvida et al., 2018). It may be that factors such as difficulties in communication, adherence to treatment or the higher cost of modern drugs contribute to a lower frequency of prescription in these population groups (Ali et al., 2022; Lamprea-Montealegre et al., 2022; Montvida et al., 2018). Nevertheless, it is important to note that in Spain, all AD are classified as reduced-contribution treatments, whereby the patient is required to pay a fixed contribution of 10% of the retail price of the package. Furthermore, the maximum payable amount per package is set at 4.24€.
In our study, patients with comorbidities such as heart failure, ischemic heart disease, chronic renal failure, obesity and a previous MACE had a higher probability of having a prescription of modern AD. Novel AD SGLT2i and GLP1a, have shown impressive CV benefits, particularly in reducing heart failure hospitalization, CV death, and demonstrated beneficial effects on cardiac function (Nikolaidou et al., 2024). In adults with T2D, SGLT2i and GLP1a (but not DPP4i) reduce all-cause mortality and MACE compared with usual care. SGLT2i reduce chronic renal failure progression and heart failure hospitalization and GLP1a reduce stroke compared with usual care. Serious adverse events and severe hypoglycemia are less frequent with SGLT2i and GLP1a than with insulin or sulfonylureas (Qaseem et al., 2024).
In light of these considerations, clinical guidelines have been revised to designate these drugs as the first treatment option for this patient population (Davies et al., 2022). Until 2022, metformin was the primary treatment option, with treatment selection based on HbA1c levels (Bae et al., 2022). At present, comorbidities are given precedence, followed by consideration of the HbA1c level.
The majority of our patients started treatment in monotherapy and were most frequently prescribed metformin as similarly as described in a study in Catalonia (Ouchi et al., 2023). In this study also patients starting in combination therapy had higher HbA1c levels than those starting monotherapy. In clinical trials, initial combination therapy showed better glycemic control than monotherapy or a stepwise approach (Kim and Kim, 2024).
With the development of new classes of AD, the paradigm of therapeutic options for patients with glycemic and CV risk factors has shifted significantly. However, the way in which this has occurred in actual practice, especially in terms of the balance between older and newer ADs as initial and intensification therapeutic options, has undergone a much slower change (Montvida et al., 2018; Carney et al., 2022). According to a UK study (Young et al., 2023), the recent UK T2D treatment guideline represents a quasi-population indication for treatment with SGLT2i. In addition to these findings, an evaluation in the Netherlands showed that the prescription rate of SGLT2i was lower than would be indicated based on the CV risk of the patients in the absence of contraindications to their use (Hart et al., 2023). Both studies showed that these treatments were prescribed less than needed, despite being indicated for many more patients.
According to our results, the most deprived areas and non-depopulated areas presented less probability of prescription of a modern AD. Previous studies have shown that higher deprivation level is associated with higher prevalence and incidence of T2D (36,37), chronic diabetes complications, indicators of cardiometabolic risk (Laraia et al., 2012), mortality, and hospitalization due to diabetes complications (Choi et al., 2020). It is also related to worse disease control (Durfey et al., 2019) and higher treatment use (Frazer and Frazer, 2020; Di Filippo et al., 2022). It is suggested that more expensive drugs were favoured in areas with lower deprivation even for similar indications. Possible explanations include an unconscious prescriber bias towards newer and more expensive medications in affluent areas, or a more educated population likely to request newer medications (Frazer and Frazer, 2020). This phenomenon has been described for numerous chronic conditions (Di Filippo et al., 2022) and it can widen the gap in disease management and outcomes.
Regarding the multilevel model, it explained 8.5% of the variability, of which 7.3% corresponded to individual variables and 1.2% to BHA variables. It was similar when differentiating between women and men. This means that belonging to an area had very little effect on the prescription of a modern AD and that the variability was explained more by individual patient factors.
4.1 Strengths and limitations
The main strength of this study lies in the fact that we selected all the individuals diagnosed with T2D from 2018 to 2022, including data from administrative health data sources and electronic health records. A further key assumption is that the prescriptions we analyzed represent the majority burden of prescribing but it is important to consider that we could not collect private prescriptions. Nonetheless, they make up a minority of prescriptions, especially among chronic patients. According to the Health system review from European Observatory on Health Systems and Policies (Bernal-Delgado et al., 2024) it is described that the number of privately insured patients in 2022 in Spain only represented the 0.4% of the Spanish population aged 18 and over.
It is also necessary to take into account some limitations inherent to observational studies, such as the quality of the data, or the existence of incomplete cases. The analysis did not include diabetic patients who had been prescribed a modern drug prior to T2D diagnosis. Insidious onset of diabetes and difficulty in labelling patients with these diagnoses, may have misclassified some individuals, but this is expected to be a small number. It has also not been possible to differentiate between the different T2D subtypes. In addition, the patients included have different possible follow-up periods, ranging from 5 years to 0 days, depending on when they were diagnosed.
Modern AD have demonstrated their superiority in certain pathologies, such as heart failure, chronic renal failure or in patients with CV risk factors or obesity, to the extent that their prescription is even indicated as first-line therapy. Although these changes in the algorithms did not take place until 2023, their prescription was already more frequent with these pathologies in clinical practice. This study provides real-world evidence that the utilization pattern of AD in Patients with T2D in Aragon reflects the anticipation of clinical practice with regard to the updating of international guideline recommendations.
The findings largely reflect the prescribing patterns observed in previous research (Mata-Cases et al., 2016; Rojo-Martínez et al., 2020; Ouchi et al., 2023). The main innovation in this respect in our study is the focus on differences about deprivation index, depopulation level and family medicine pressure.
5 Conclusion
Sociodemographic, clinical and regional factors are associated with variability in the prescribing pattern. Our findings are helpful for understanding the use of these drugs and the predictors of their prescription, in order to better assess and understand the health outcomes of people with T2D.
The implementation of current clinical guidelines needs to be adapted to explicitly consider socio-economic and territorial factors. Although guidelines do not usually incorporate aspects such as deprivation, depopulation or care pressure, these determinants should be considered when applying them to ensure more equitable clinical practice. To this end, specific strategies such as prescription audits, targeted professional training and clinical decision support systems should be incorporated in areas with higher levels of inequality or family medicine pressure.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by the Research Ethics Committee of Aragón (CEICA, PI21/148). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
IB-U: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review and editing. IA-P: Conceptualization, Methodology, Supervision, Writing – original draft, Writing – review and editing. MR: Writing – review and editing. ML: Writing – review and editing. AG: Data curation, Formal Analysis, Methodology, Writing – review and editing. SM: Conceptualization, Methodology, Supervision, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. Supported by Instituto de Salud Carlos III (ISCIII), with group funds the Research Network on Chronicity, Primary Care and Health Promotion (RICAPPS, RD21/0016/0005) that is part of the Results-Oriented Cooperative Research Networks in Health (RICORS) (Carlos III Health Institute), funded by the European Union “NextGeneration EU/PRTR” and funded by Carlos III Health Institute (ISCIII), grant numbers PI22/01193 with the European Regional Development Fund (EDRF) “A way of making Europe” and co-funding by European Union. It is also supported by the Research Group in Health Services in Aragón (GRISSA) (B09-23R) that is part of the Department of Innovation, Research and University at the Government of Aragón (Spain).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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/fphar.2025.1530139/full#supplementary-material
References
Aguilar-Palacio, I., Rabanaque, M., Castel, S., Maldonado, L., González-García, J., Compés, L., et al. (2024). Cohort profile: the CArdiovascular risk factors for hEalth services research (CARhES) cohort study. Int. J. Epidemiol. 53 (2), dyae015. doi:10.1093/ije/dyae015
Alemán, J. J., Artola, S., Franch, J., Mata, M., Millaruelo, J. M. y, and Sangrós, J. ENDLR (2014). Recomendaciones para el tratamiento de la diabetes mellitus tipo 2: control glucémico. Algoritm Trat hiperglucemia diabetes tipo 2, 1–2.
Alemán, J. J., Menéndez, S. A., Lachica, L. Á., De La Puente, J. B., Barutell, L., Benito, B., et al. (2018). 5 Guía de diabetes tipo 2 para clínicos: recomendaciones de la redGDPS. Available online at: https://www.redgdps.org/gestor/upload/colecciones/GuiaDM2_web.pdf.
Ali, A. S., Khan, M., Ortega, D., Jaffery, S., Lamerato, L. E., and Budzynska, K. (2022). Prescription patterns of novel antihyperglycemic medications. J. Am. Board Fam. Med. 35 (2), 255–264. doi:10.3122/jabfm.2022.02.210360
Bae, J. P., Kadziola, Z. A., Liu, D., Chinthammit, C., Boye, K. S., and Mather, K. J. (2022). An early assessment of the real-world treatment patterns of type 2 diabetes: a comparison to the 2018 ADA/EASD consensus report recommendations. Diabetes Ther. 13 (8), 1499–1510. doi:10.1007/s13300-022-01289-x
Barrio-Cortes, J., del Cura-González, I., Martínez-Martín, M., López-Rodríguez, C., Jaime-Sisó, M. Á., and Suárez-Fernández, C. (2020). Atención Primaria comorbilidades de los pacientes crónicos según nivel. 52(2):86–95.
Bates, D., Mächler, M., Bolker, B., and Walker, S. (2015). Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67 (1 SE-Articles), 1–48. doi:10.18637/jss.v067.i01
Bernal-Delgado, E., Garcia-Armesto, S., Oliva, J., Sanchez Martinez, F. I., Repullo, J. R., Pena-Longobardo, L. M., et al. (2024). Spain: Health system review. Health Syst. transition 26, 1–179.
Carney, G., Kim, J. D., O’Sullivan, C., Thompson, W., Bassett, K., Levin, J., et al. (2022). Treatment pattern trends of medications for type 2 diabetes in British Columbia, Canada. BMJ Open Diabetes Res. Care 10 (6), e002995–12. doi:10.1136/bmjdrc-2022-002995
Choe, S. A., Kim, J. Y., Ro, Y. S., and Cho, S. I. (2018). “Women are less likely than men to achieve optimal glycemic control after 1 year of treatment: a multi-level analysis of a Korean primary care cohort,” Editor A. Palazón-Bru, 13. doi:10.1371/journal.pone.0196719
Choi, D. W., Lee, S. A., Lee, D. W., Joo, J. H., Han, K. T., Kim, S., et al. (2020). Effect of socioeconomic deprivation on outcomes of diabetes complications in patients with type 2 diabetes mellitus: a nationwide population-based cohort study of South Korea. BMJ open diabetes Res. care 8 (1), e000729. doi:10.1136/bmjdrc-2019-000729
Compés Dea, M. L., Olivan Bellido, E., Feja Solana, C., Aguilar Palacio, I., García-Carpintero Romero Del Hombrebueno, G., and Adiego Sancho, B. (2018). Construction of a deprivation index by basic healthcare area in Aragon using population and housing census 2011. Rev. Esp. Salud Publica 92, e201812087.
Dahlén, A. D., Dashi, G., Maslov, I., Attwood, M. M., Jonsson, J., Trukhan, V., et al. (2022). Trends in antidiabetic drug discovery: FDA approved drugs, new drugs in clinical trials and global sales. Front. Pharmacol. 12 (January), 1–16. doi:10.3389/fphar.2021.807548
Davies, M. J., Aroda, V. R., Collins, B. S., Gabbay, R. A., Green, J., Maruthur, N. M., et al. (2022). Management of hyperglycemia in type 2 diabetes, 2022. A consensus report by the american diabetes association (ADA) and the european association for the study of diabetes (EASD). Diabetes Care 45 (11), 2753–2786. doi:10.2337/dci22-0034
DeMarsilis, A., Reddy, N., Boutari, C., Filippaios, A., Sternthal, E., Katsiki, N., et al. (2022). Pharmacotherapy of type 2 diabetes: an update and future directions. Metabolism 137 (August), 155332. doi:10.1016/j.metabol.2022.155332
Di Filippo, A., Perna, S., Pierantozzi, A., Milozzi, F., Fortinguerra, F., Caranci, N., et al. (2022). Socio-economic inequalities in the use of drugs for the treatment of chronic diseases in Italy. Int. J. Equity Health 21 (1), 157. doi:10.1186/s12939-022-01772-8
Durfey, S. N. M., Kind, A. J. H., Buckingham, W. R., DuGoff, E. H., and Trivedi, A. N. (2019). Neighborhood disadvantage and chronic disease management. Health Serv. Res. 54 (S1), 206–216. doi:10.1111/1475-6773.13092
Frazer, J. S., and Frazer, G. R. (2020). GP prescribing in northern Ireland by deprivation index: retrospective analysis. Fam. Med. Community Heal 8 (3), e000376. Available online at: http://fmch.bmj.com/lookup/doi/10.1136/fmch-2020-000376.
García Soidán, F. J. (2023). Algoritmo de tratamiento de la diabetes mellitus tipo 2. RedGDPS 2023. Diabetes práctica. 14 (2), 43–52. doi:10.52102/diabet/pract.2023.2.art2
Guion, M., Mandereau-Bruno, L., Goria, S., Cosson, E., and Fosse-Edorh, S. (2024). Eleven-year trends in socioeconomic inequalities in the prevalence and incidence of pharmacologically treated type 2 diabetes in France, 2010–2020. Diabetes Metab. 50 (2), 101509. doi:10.1016/j.diabet.2024.101509
Hart, H. E., Kievits, O., Rutten, F. H., and Hollander, M. H. (2023). Evaluation of SGLT-2 inhibitor treatment in type 2 diabetes patients with very high cardiovascular risk. Prim. Care Diabetes 17 (2), 190–194. doi:10.1016/j.pcd.2023.02.001
Herrington, W. G., Staplin, N., Wanner, C., Green, J. B., Hauske, S. J., Emberson, J. R., et al. (2023). Empagliflozin in patients with chronic kidney disease. N. Engl. J. Med. 388 (2), 117–127. doi:10.1056/NEJMoa2204233
Internacional Diabetes Federation (2023). IDF diabetes atlas. IDF Diabetes Atlas 10th Ed. 537 million people Worldw. have diabetes 1, 31–88. Available online at: www.diabetesatlas.org.
Kim, J. Y., and Kim, N. H. (2024). Initial combination therapy in type 2 diabetes. Endocrinol. Metab. 39 (1), 23–32. Available online at: http://e-enm.org/journal/view.php?doi=10.3803/EnM.2023.1816.
Lamprea-Montealegre, J. A., Madden, E., Tummalapalli, S. L., Peralta, C., Neilands, T. B., Garcia, P. K., et al. (2022). Association of race and ethnicity with prescription of SGLT2 inhibitors and GLP1 receptor agonists among patients with type 2 diabetes in the veterans health administration system. JAMA 328 (9), 861–871. doi:10.1001/jama.2022.13885
Laraia, B. A., Karter, A. J., Warton, E. M., Schillinger, D., Moffet, H. H., and Adler, N. (2012). Place matters: neighborhood deprivation and cardiometabolic risk factors in the diabetes study of northern California (DISTANCE). Soc. Sci. Med. 74 (7), 1082–1090. doi:10.1016/j.socscimed.2011.11.036
Mata, C. M., Artola Menéndez, S., Díez Espino, J., Ezkurra Loiola, P., Franch Nadal, J., and García Soidán, F. J. (2020). Actualización de 2020 del algoritmo de tratamiento de la hiperglucemia en la diabetes mellitus tipo 2 de la redGDPS. Diabetes Práctica 11 (02), 41–76. Available online at: https://www.redgdps.org/algoritmo-de-tratamiento-de-la-.
Mata-Cases, M., Franch-Nadal, J., Real, J., and Mauricio, D. (2016). Glycaemic control and antidiabetic treatment trends in primary care centres in patients with type 2 diabetes mellitus during 2007-2013 in Catalonia: a population-based study. BMJ Open 6 (10), e012463–11. doi:10.1136/bmjopen-2016-012463
McDonagh, T. A., Metra, M., Adamo, M., Gardner, R. S., Baumbach, A., Böhm, M., et al. (2022). Guía ESC 2021 sobre el diagnóstico y tratamiento de la insuficiencia cardiaca aguda y crónica. Rev. Española Cardiol. 75 (6), 523.e1–523.e114. Available online at: https://linkinghub.elsevier.com/retrieve/pii/S0300893221005236.
Mezuk, B., Chaikiat, Å., Li, X., Sundquist, J., Kendler, K. S., and Sundquist, K. Depression, neighborhood deprivation and risk of type 2 diabetes. Available online at: https://europepmc.org/articles/PMC3773725.
Montvida, O., Shaw, J., Atherton, J. J., Stringer, F., and Paul, S. K. (2018). Long-term trends in antidiabetes drug usage in the U.S.: real-World evidence in patients newly diagnosed with type 2 diabetes. Diabetes Care 41 (1), 69–78. doi:10.2337/dc17-1414
Moreno-Juste, A., Poblador-Plou, B., Aza-Pascual-Salcedo, M., González-Rubio, F., Malo, S., Librero López, J., et al. (2020). Initial therapy, regimen change, and persistence in a Spanish cohort of newly treated type 2 diabetes patients: a retrospective, observational study using real-world data. Int. J. Environ. Res. Public Health 17 (10), 3742. doi:10.3390/ijerph17103742
Nikolaidou, A., Ventoulis, I., Karakoulidis, G., Anastasiou, V., Daios, S., Papadopoulos, S. F., et al. (2024). Hypoglycemic drugs in patients with diabetes mellitus and heart failure: a narrative review. Med. (B Aires) 60 (6), 912. doi:10.3390/medicina60060912
Ouchi, D., Giner-Soriano, M., Vilaplana-Carnerero, C., Monfa, R., Torres, F., and Morros, R. (2023). Longitudinal treatment patterns in patients recently diagnosed with type 2 diabetes mellitus in Catalonia. Diabetes Res. Clin. Pract. 202 (May), 110777. doi:10.1016/j.diabres.2023.110777
Palmer, S. C., Tendal, B., Mustafa, R. A., Vandvik, P. O., Li, S., Hao, Q., et al. (2021). Sodium-glucose cotransporter protein-2 (SGLT-2) inhibitors and glucagon-like peptide-1 (GLP-1) receptor agonists for type 2 diabetes: systematic review and network meta-analysis of randomised controlled trials. BMJ 372, m4573. Available online at: https://www.bmj.com/lookup/doi/10.1136/bmj.m4573.
Qaseem, A., Obley, A. J., Shamliyan, T., Hicks, L. A., Harrod, C. S., Crandall, C. J., et al. (2024). Newer pharmacologic treatments in adults with type 2 diabetes: a clinical guideline from the American college of physicians. Ann. Intern. Med. 177, 658–666. doi:10.7326/M23-2788
Rojo-Martínez, G., Valdés, S., Soriguer, F., Vendrell, J., Urrutia, I., Pérez, V., et al. (2020). Incidence of diabetes mellitus in Spain as results of the nation-wide cohortZGlAYmV0LmVzstudy. Sci. Rep. 10 (1), 1–9.
Spanish ministry (2023). Spanish ministry for the ecological transition and the demographic challenge. Madrid: Government of Spain. Demographic Challenge. Available online at: https://www.miteco.gob.es/es/reto-demografico/temas/.
Ussher, J. R., and Drucker, D. J. (2023). Glucagon-like peptide 1 receptor agonists: cardiovascular benefits and mechanisms of action. Nat. Rev. Cardiol. 20 (7), 463–474. doi:10.1038/s41569-023-00849-3
van der Aart-van der Beek, A. B., de Boer, R. A., and Heerspink, H. J. L. (2022). Kidney and heart failure outcomes associated with SGLT2 inhibitor use. Nat. Rev. Nephrol. 18 (5), 294–306. doi:10.1038/s41581-022-00535-6
Wang, C., Gao, Y., Zhu, L., Huang, M., Wu, Y., and Xuan, J. (2019). Treatment patterns in patients with newly diagnosed type 2 diabetes in China: a retrospective, longitudinal database study. Clin. Ther. 41 (8), 1440–1452. doi:10.1016/j.clinthera.2019.05.003
Keywords: antidiabetic drug, type 2 diabetes, prescribing pattern, multilevel analysis, real word data
Citation: Bengoa-Urrengoechea I, Aguilar-Palacio I, Rabanaque MJ, Lallana MJ, Gamba Cabezas A and Malo S (2025) Factors associated with prescription of modern antidiabetics in newly diagnosed patients with type 2 diabetes. a real-world data study in a Spanish region. Front. Pharmacol. 16:1530139. doi: 10.3389/fphar.2025.1530139
Received: 18 November 2024; Accepted: 26 June 2025;
Published: 11 July 2025.
Edited by:
Anick Bérard, Montreal University, CanadaReviewed by:
Shaohua Qi, Houston Methodist Research Institute, United StatesPetya Milushewa, Medical University of Sofia, Bulgaria
Slobodan Jankovic, University of Kragujevac, Serbia
Copyright © 2025 Bengoa-Urrengoechea, Aguilar-Palacio, Rabanaque, Lallana, Gamba Cabezas and Malo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Isabel Aguilar-Palacio, aWFndWlsYXJAdW5pemFyLmVz
†ORCID: Irantzu Bengoa-Urrengoechea, orcid.org/0000-0001-7283-7713; Isabel Aguilar-Palacio, orcid.org/0000-0001-7293-701X; María José Rabanaque, orcid.org/0000-0002-6671-5661; María Jesús Lallana, orcid.org/0000-0002-3803-5298; Adriana Gamba Cabezas, orcid.org/0009-0005-8707-1841; Sara Malo, orcid.org/0000-0002-7194-8275