Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Pharmacol., 11 July 2025

Sec. Pharmacoepidemiology

Volume 16 - 2025 | https://doi.org/10.3389/fphar.2025.1530139

Factors associated with prescription of modern antidiabetics in newly diagnosed patients with type 2 diabetes. a real-world data study in a Spanish region

Irantzu Bengoa-Urrengoechea,&#x;Irantzu Bengoa-Urrengoechea1,2Isabel Aguilar-Palacio,,
&#x;Isabel Aguilar-Palacio2,3,4*María Jos Rabanaque,,&#x;María José Rabanaque2,3,4María Jesús Lallana,&#x;María Jesús Lallana2,5Adriana Gamba Cabezas&#x;Adriana Gamba Cabezas2Sara Malo,,&#x;Sara Malo2,3,4
  • 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
www.frontiersin.org

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).

VCisjk=logπsjk1πsjk=β0+Xβisjk+us+uj+uk+eisjk

The model had the following assumptions: first, the random effects us , uj and uk were normally distributed with mean 0 and variance σu2, second, the error component eisjk was also normally distributed with mean 0 and variance. σe2; Third, the random effects us, uj and uk and the error component eisjk were independent, and eisjk were all independent of each other. With i = 1, 2, n; s = 1, 2, three level of depopulation; j = 1, 2,3, 4 Quartiles deprivation index and k = 1,2,3,4 Quartiles family medicine pressure. Given the characteristics of this multilevel study, there was an intraclass correlation, which means that there were observations that are more similar to others in the same group than to those in other groups. Variance partition coefficients were calculated to see how much of the response variance belongs to each level. To evaluate statistical significance, a p-value smaller than 0.05 was used.

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.

Table 1
www.frontiersin.org

Table 1. Description of the total sample analyzed and comparison by sex.

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
www.frontiersin.org

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
www.frontiersin.org

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.

Figure 2
www.frontiersin.org

Figure 2. Random effects of basic healthcare area variables.

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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.

PubMed Abstract | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

PubMed Abstract | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

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.

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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-.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

PubMed Abstract | CrossRef Full Text | Google Scholar

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/.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, K. G., Hopkins, R., Shields, B. M., Thomas, N. J., McGovern, A. P., Dennis, J. M., et al. (2023). Recent UK type 2 diabetes treatment guidance represents a near whole population indication for SGLT2-inhibitor therapy. Cardiovasc Diabetol. 22 (1), 302–303. doi:10.1186/s12933-023-02032-x

PubMed Abstract | CrossRef Full Text | Google Scholar

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, Canada

Reviewed by:

Shaohua Qi, Houston Methodist Research Institute, United States
Petya 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

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