Elena Flowers
Alka M. Kanaya3,4
Li Zhang- 1University of California, San Francisco, Department of Physiological Nursing, San Francisco, CA, United States
- 2University of California, San Francisco, Institute for Human Genetics, San Francisco, CA, United States
- 3University of California, San Francisco, Department of Epidemiology and Biostatistics, San Francisco, CA, United States
- 4University of California, San Francisco, Department of Medicine, Division of General Internal Medicine, San Francisco, CA, United States
- 5New York University Bluestone Center for Clinical Research, New York, NY, United States
- 6New York University Department of Oral and Maxillofacial Surgery, New York, NY, United States
Prior studies focused on circulating microRNAs and the risk for complex diseases have shown inconsistent findings. The majority of studies focused on European and East Asian racial or ethnic groups, however, ancestry was not typically reported. We evaluated the risk for type 2 diabetes as an exemplar to show that race and ethnic group may contribute to inconsistent validation of previous findings of associations with microRNAs.
MicroRNAs (miRs) are emerging as “intermediary” biomarkers that capture the combined effect of gene-environment interactions. Circulating miRs found in serum and plasma are easily measured in blood and may have utility for determining risk for the development of complex diseases. More than 37 million people (11.3% of the population) in the United States have type 2 diabetes (T2D) and an additional 96 million (38% of the population) have prediabetes. (Centers for Disease Control and Prevention, 2022). However, there are marked differences in prevalence across racial and ethnic groups in the United States and globally. Given the complex etiology of T2D, both genetic (at times reflecting ancestry) risk factors as well as social determinants likely contribute to these observed differences. A number of prior studies identified associations between circulating miRs and risk for T2D. However, the findings from these studies are disparate. As with genetic association studies, one reason for these disparate findings includes demographic differences in the population sample (e.g., ancestral or genetic, racial or ethnicity). If expression of miRs is influenced by genetic ancestry as estimated by ancestry informative markers (AIMs), then defining a T2D miR “signature” may also reflect information about underlying genetic ancestry. In the absence of genetic information about a study sample, this signature could provide the opportunity to control for confounding based on differences in ancestry, discrete from differences in race or ethnicity, which would help parse out the separate impacts of genetic risk and social determinants (e.g., racism) on complex diseases when AIMs are not available. Given that miRs are responsive to environmental and behavioral factors, these biomarkers may also provide insights as to how ancestry impacts responses to social determinants across racial and ethnic groups. The purpose of this study was to test whether there were differences in associations between miRs and risk for T2D by racial or ethnic group that may reflect underlying confounding due to the impact of social determinants of health.
This study was a secondary analysis of data and biospecimens from participants in the Diabetes Prevention Program (DPP) trial that tested two approaches to prevention of T2D after 2-years. DPP sample characteristics, trial design, and methods have been described in detail previously. (The Diabetes Prevention Program, 1999; The Diabetes Prevention Program, 2000; The Diabetes Prevention Program, 2002). This analysis includes a randomly selected subset (n = 1,000) equally stratified by the intervention arm. We performed a literature review to identify miRs that were previously associated with prediabetes [100 mg/dl < fasting blood glucose (FBG) < 126 mg/dl)]. MicroRNAs are named in sequential order of entry into an annotation database, with the -3p or -5p suffix referring to which arm of the stem loop structure is represented (Ambros et al., 2003). From all miRs reported, we selected those that were statistically significant in at least two prior studies. Of those, miR-126-3p, miR-126-5p, miR-144, miR-192, and miR-23a were present in our existing dataset and were included in this study. The Multiplex Circulating MicroRNA Assay (Abcam, MA) was used to directly quantify miRs from plasma. Descriptive statistics were calculated to examine and evaluate the demographic and clinical characteristics of participants. (R, 2019) Expression of individual miRs was normalized using the set of miR probes (i.e., hsa-miR-92a-3p, hsa-mir-93-5p, hsa-miR-17-5p) identified by the geNorm algorithm for each experiment (Ambros et al., 2003). Spearman’s correlation coefficients were calculated to determine the associations between individual miRs and FBG (R, 2019). Subgroup tests were performed by sex, race, and ethnicity.
Sample characteristics are shown in Table 1. None of the miRs that we assessed were significantly associated with FBG in the study sample overall (Figure 1). However, when we stratified by sex, race, and ethnicity, we observed significant associations for miR-126-3p (ρ = 0.3, p < 0.05) in Hispanic men, miR-192 (ρ = 0.26, p < 0.05) in Black women, and miR-23a in males overall (ρ = −0.11, p < 0.05) and White males (r2 = 0.16, p < 0.05).
TABLE 1. Demographic and clinical characteristics.
FIGURE 1. Rows show race/ethnic group. Columns show microRNAs assessed. Color intensity shows the magnitude of the correlation. Statistically significant values are denoted with an asterisk.
The scientific community is reaching a consensus that race and ethnicity should be operationalized as a social rather than biological construct (Flanagin et al., 2021). This is a critical step forward for conducting research focused on the impact of racism and other social determinants on health outcomes, distinct from the impact of genomic characteristics on health outcomes. While both categories of risk factors (i.e., social, genetic) are important for understanding the etiology of complex diseases like T2D, and which are often correlated, they should not be conflated. Because miRs are encoded from the genome, variability in the genetic sequence that is observed among ancestral groups may also be relevant to the genomic regions encoding for miRs, and therefore differences in expression and function of individual miRs. Prior studies focused on miRs have not typically applied robust approaches to account for this variability using AIMs. A major challenge to designing studies that control for ancestry is the need for genetic data. In contrast, one of the most promising features of circulating miRs is that they can be easily measured from blood collected in clinical settings and do not require nuclear genetic material. As has been used using protein-encoding genes, expression levels of a subset of miRs that reflects ancestral genetic influences. This “fingerprint” could then be used to control for underlying genetic differences between individuals that may impact observed associations between miRs and risk for complex traits, including T2D. This would augment the established ease of measuring miRs with methods to account for underlying genetic differences confounding observed statistical associations. While we did not validate prior findings on associations between individual miRs and risk for T2D, we showed that some of these associations were significant when we stratified by racial or ethnic group and sex. Prior studies have not consistently reported these characteristics, and no studies reported ancestral characteristics. As is well established in genetic association studies, differences in underlying genetic characteristics is a likely contributor to inconsistent validation of miRs across studies.
The strongest observed association was between miR-126-3p and FBG in Hispanic men. MiR-126 was one of the first miRs identified in a study focused on T2D, (Zampetaki et al., 2010), and numerous studies showed that this miR is likely related to endothelial function and complications from T2D and not the etiology of T2D itself (Zampetaki et al., 2010). Prior studies on T2D-related complications showed inconsistent relationships with vascular complications in men who identify as Hispanic (Haw et al., 2021). Our prior study did not provide evidence that miR-126 was significantly associated with prediabetes versus normal glucose tolerance in independent studies of United States Latin and Mexican origin participants (Flowers et al., 2021). Failure to adequately control for ancestry may contribute to these disparate findings, and there may be disproportionate risk in this group because of shared ethnic-based exposures that could be confounded by genetic (ancestral) differences between distinct groups of men who identify as Hispanic.
Similarly, we saw a moderate correlation between miR-192 and FBG in Black women. A recent review showed evidence that miR-192 is related to many diseases, including two studies of T2D, by targeting the forkhead box 01 (FOX O 1) gene to repress insulin resistance (Ren et al., 2021). One was performed on the European continent, although participants’ racial and ethnic characteristics were not described (Jaeger et al., 2018). The other was performed in the United States, and 83% of the sample self-reported as White with the rest of the sample categorized as Other (Nunez Lopez et al., 2019). Thus, we cannot determine how ancestry, race, or ethnicity are related to expression levels of miR-192 and risk for T2D.
This study also showed that miR-23a was associated with FBG in men overall and also White men. A prior study found associations between miR-23 in adipose tissue and insulin resistance (Lozano-Bartolomé et al., 2018). The study took place in Spain and approximately half the participants were reported as male. Several prior studies identified miR-23a associated with T2D-related retinopathy (Zhao et al., 2016; Santovito et al., 2021) and β-cell destruction in type 1 diabetes (Zhao et al., 2016; Grieco et al., 2017). Most of these studies were conducted on the European and Asian continents. Recent studies have shown that T2D is a heterogeneous condition, with several subtypes that are characterized by distinct clusters of characteristics (e.g., insulin resistance vs low insulin secretion). (Ahlqvist et al., 2018; Zou et al., 2019; Anjana et al., 2020; Herder et al., 2021). These subtypes feature varying risk for T2D-related complications like retinopathy (Grieco et al., 2017) and evidence for genetic differences between these subgroups of patients (Ahlqvist et al., 2018). Future studies that characterize genetic ancestry may help to inform whether there are genetic differences in miRs associated with T2D subgroups by ancestral origin.
We limited our assessment to the subset of miRs identified in at least two prior studies of prediabetes; additional miR are certain to be relevant. Historically, studies have often used self-reported race and ethnicity as a proxy for genetic ancestry. A more accurate etiological model would differentiate race and ethnicity as a social construct, discrete from genetic ancestry as a biological construct, when assessing for potential confounding (Borrell et al., 2021). The observed associations between miRs and FBG for some subgroups by race and ethnicity and sex in this study cannot parse the underlying risk (i.e., social determinants, genetic, both) without AIMS. For population-based studies of T2D that do not have AIM data, the establishment of a miR “fingerprint” that reflects underlying genetic differences would prove instrumental to miR studies that can estimate the contributions of social vs biological impacts on risk for T2D. Cohorts that possess AIM and miR data can be used to develop and validate such a miR fingerprint.
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 human participants were reviewed and approved by University of California, San Francisco. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
EF was the Principal Investigator of the grant that supported this project. EF conceived of and executed the overall study and design, oversaw all molecular data collection, directed the analysis plan, and wrote and revised the manuscript. AK was the Principal Investigator of the PRYSMS trial and oversaw all clinical data collection, contributed to the design and interpretation of results for this study, and approved the final manuscript. LZ contributed to the analysis plan, contributed to the interpretation of the results, and approved the final manuscript. BA provided overall scientific guidance for this study, contributed to interpretation of the results, and approved the final manuscript.
Funding
This study was supported by the National Institute for Diabetes, Digestive and Kidney Disease grant number R01DK124228. Kanaya is supported by National Heart Lung, and Blood Institute of the National Institutes of Health grant number K24HL112827.
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.
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.
References
Ahlqvist, E., Storm, P., Käräjämäki, A., Martinell, M., Dorkhan, M., Carlsson, A., et al. (2018). Novel Subgroups of Adult-Onset Diabetes and Their Association with Outcomes: a Data-Driven Cluster Analysis of Six Variables. Lancet Diabetes Endocrinol. 6, 361–369. doi:10.1016/s2213-8587(18)30051-2
Ambros, V., Bartel, B., Bartel, D. P., Burge, C. B., Carrington, J. C., Chen, X., et al. (2003). A Uniform System for microRNA Annotation. RNA 9 (3), 277–279. doi:10.1261/rna.2183803
Anjana, R. M., Baskar, V., Nair, A. T. N., Jebarani, S., Siddiqui, M. K., Pradeepa, R., et al. (2020). Novel Subgroups of Type 2 Diabetes and Their Association with Microvascular Outcomes in an Asian Indian Population: a Data-Driven Cluster Analysis: the INSPIRED Study. BMJ Open Diabetes Res. Care 8, e001506. doi:10.1136/bmjdrc-2020-001506
Borrell, L. N., Elhawary, J. R., Fuentes-Afflick, E., Witonsky, J., Bhakta, N., Wu, A. H. B., et al. (2021). Race and Genetic Ancestry in Medicine - A Time for Reckoning with Racism. N. Engl. J. Med. 384, 474–480. doi:10.1056/nejmms2029562
Centers for Disease Control and Prevention (2022). National Diabetes Statistics Report Website. Available at: https://www.cdc.gov/diabetes/data/statistics-report/index.html (Accessed February 7, 2022).
Flanagin, A., Frey, T., and Christiansen, S. L.AMA Manual of Style Committee (2021). Updated Guidance on the Reporting of Race and Ethnicity in Medical and Science Journals. JAMA 326, 621–627. doi:10.1001/jama.2021.13304
Flowers, E., Ramírez-Mares, J. D., Velazquez-Villafaña, M., Rangel-Salazar, R., Sucher, A., Kanaya, A. M., et al. (2021). Circulating microRNAs Associated with Prediabetes and Geographic Location in Latinos. Int. J. Diabetes Dev. Ctries 41 (4), 570–578. doi:10.1007/s13410-020-00917-1
Grieco, F. A., Sebastiani, G., Juan-Mateu, J., Villate, O., Marroqui, L., Ladrière, L., et al. (2017). MicroRNAs miR-23a-3p, miR-23b-3p, and miR-149-5p Regulate the Expression of Proapoptotic BH3-Only Proteins DP5 and PUMA in Human Pancreatic β-Cells. Diabetes 66, 100–112. doi:10.2337/db16-0592
Haw, J. S., Shah, M., Turbow, S., Egeolu, M., and Umpierrez, G. (2021). Diabetes Complications in Racial and Ethnic Minority Populations in the USA. Curr. Diab Rep. 21, 2. doi:10.1007/s11892-020-01369-x
Herder, C., Maalmi, H., Strassburger, K., Zaharia, O. P., Ratter, J. M., Karusheva, Y., et al. (2021). Differences in Biomarkers of Inflammation between Novel Subgroups of Recent-Onset Diabetes. Diabetes 70 (5), 1198–1208. doi:10.2337/db20-1054
Jaeger, A., Zollinger, L., Saely, C. H., Muendlein, A., Evangelakos, I., Nasias, D., et al. (2018). Circulating microRNAs -192 and -194 Are Associated with the Presence and Incidence of Diabetes Mellitus. Sci. Rep. 8, 14274. doi:10.1038/s41598-018-32274-9
Lozano-Bartolomé, J., Llauradó, G., Portero-Otin, M., Altuna-Coy, A., Rojo-Martínez, G., Vendrell, J., et al. (2018). Altered Expression of miR-181a-5p and miR-23a-3p Is Associated with Obesity and TNFα-Induced Insulin Resistance. J. Clin. Endocrinol. Metab. 103, 1447–1458. doi:10.1210/jc.2017-01909
Nunez Lopez, Y. O., Retnakaran, R., Zinman, B., Pratley, R. E., and Seyhan, A. A. (2019). Predicting and Understanding the Response to Short-Term Intensive Insulin Therapy in People with Early Type 2 Diabetes. Mol. Metab. 20, 63–78. doi:10.1016/j.molmet.2018.11.003
Ren, F.-j., Yao, Y., Cai, X.-y., and Fang, G.-y. (2021). Emerging Role of MiR-192-5p in Human Diseases. Front. Pharmacol. 12, 614068. doi:10.3389/fphar.2021.614068
Santovito, D., Toto, L., De Nardis, V., Marcantonio, P., D’Aloisio, R., Mastropasqua, A., et al. (2021). Plasma microRNA Signature Associated with Retinopathy in Patients with Type 2 Diabetes. Sci. Rep. 11, 4136. doi:10.1038/s41598-021-83047-w
The Diabetes Prevention Program (DPP) (2002). The Diabetes Prevention Program (DPP): Description of Lifestyle Intervention. Diabetes Care 25, 2165–2171. doi:10.2337/diacare.25.12.2165
The Diabetes Prevention Program (2000). The Diabetes Prevention Program: Baseline Characteristics of the Randomized Cohort. The Diabetes Prevention Program Research Group. Diabetes Care 23, 1619–1629. doi:10.2337/diacare.23.11.1619
The Diabetes Prevention Program (1999). The Diabetes Prevention Program. Design and Methods for a Clinical Trial in the Prevention of Type 2 Diabetes. Diabetes Care 22, 623–634. doi:10.2337/diacare.22.4.623
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate Normalization of Real-Time Quantitative RT-PCR Data by Geometric Averaging of Multiple Internal Control Genes. Genome Biol. 3, RESEARCH0034. doi:10.1186/gb-2002-3-7-research0034
Zampetaki, A., Kiechl, S., Drozdov, I., Willeit, P., Mayr, U., Prokopi, M., et al. (2010). Plasma microRNA Profiling Reveals Loss of Endothelial miR-126 and Other microRNAs in Type 2 Diabetes. Circ. Res. 107, 810–817. doi:10.1161/circresaha.110.226357
Zhao, S., Li, T., Li, J., Lu, Q., Han, C., Wang, N., et al. (2016). miR-23b-3p Induces the Cellular Metabolic Memory of High Glucose in Diabetic Retinopathy through a SIRT1-dependent Signalling Pathway. Diabetologia 59, 644–654. doi:10.1007/s00125-015-3832-0
Keywords: microRNA, diabetes, fasting blood glucose, biomarker, prediabetes, race/ethnicity, ancestry, diabetes prevention program
Citation: Flowers E, Kanaya AM, Zhang L and Aouizerat BE (2022) The Role of Racial and Ethnic Factors in MicroRNA Expression and Risk for Type 2 Diabetes. Front. Genet. 13:853633. doi: 10.3389/fgene.2022.853633
Received: 12 January 2022; Accepted: 07 March 2022;
Published: 18 March 2022.
Edited by:
Yusuf Tutar, University of Health Sciences (Turkey), TurkeyReviewed by:
Mehmet Taşpinar, Aksaray University, TurkeyBurak Barut, Karadeniz Technical University, Turkey
Copyright © 2022 Flowers, Kanaya, Zhang and Aouizerat. 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: Elena Flowers, elena.flowers@ucsf.edu