Research Article | Open Access
Lian-ji Zhou, Da-wei Yang, Li-Na Ou, Xing-Rong Guo, Biao-liang Wu, "Circulating Expression Level of LncRNA Malat1 in Diabetic Kidney Disease Patients and Its Clinical Significance", Journal of Diabetes Research, vol. 2020, Article ID 4729019, 7 pages, 2020. https://doi.org/10.1155/2020/4729019
Circulating Expression Level of LncRNA Malat1 in Diabetic Kidney Disease Patients and Its Clinical Significance
Background. Long noncoding RNA MALAT1 is closely related to diabetes and kidney diseases and is expected to be a new target for the diagnosis and treatment of diabetic nephropathy. Objective. This study aimed to explore the circulating expression level and significance of lncRNA Malat1 in patients with type 2 diabetes mellitus (T2DM) and diabetic kidney disease (DKD). Methods. Quantitative real-time PCR (qPCR) was conducted to assess the expression of lncRNA Malat1 in 20 T2DM patients, 27 DKD patients, and 14 healthy controls, and then, the clinical significance was analyzed. Results. LncRNA MALAT1 expression in peripheral blood mononuclear cells (PBMC) was significantly upregulated in T2DM and DKD groups when compared to control. Pearson’s correlation analysis showed correlation of lncRNA MALAT1 levels with ACR, urine β2-microglobulin (β2-MG), urine α1-microglobulin (α1-MG), creatinine (Cr), and glycosylated hemoglobin (HbA1c), while negative with superoxide dismutase (SOD) (, ). Binary regression analysis showed that ACR, creatinine, α1-MG, and LncRNA Malat1 were the risk factors for diabetic nephropathy with OR value of 1.166, 1.031, 1.031, and 2.019 (). The area under ROC curve (AUC) of DKD identified by the above indicators was 0.914, 0.643, 0.807, and 0.797, respectively. The AUC of Joint prediction probability of DKD recognition was 0.914, and the sensitivity and specificity of DKD diagnosis were 1.0 and 0.806, respectively. (Take ≥0.251 as the diagnostic cutoff point). Conclusion. LncRNA Malat1 is highly expressed in DKD patients, and the combined detection of ACR, creatinine, α1-MG, and LncRNA Malat1 with diabetes mellitus may be the best way to diagnose diabetic nephropathy.
Diabetic kidney disease (DKD) has become the primary etiology for chronic kidney disease in China. The morbidity and mortality associated with DKD has been rapidly increased over the past 20 years [1, 2]. According to the latest epidemiological data, DKD is considered as the primary cause for the progression to dialytic end-stage kidney diseases in the urban population [3–5]. Long noncoding RNA (lncRNA) plays an important role in the modulation of various cellular responses, system development, and pathogenesis . Recent studies have confirmed the involvement of lncRNA in DKD processes such as renal fibrosis, extracellular matrix deposition, inflammatory response, oxidative stress, and apoptosis of renal tubular epithelial cells . It played a vital role in DKD pathogenesis and progression via a series of mechanisms. MALAT1, also known as noncoding nuclear-enriched abundant transcript 2 (NEAT2), was initially detected in the year 2003 by Ji et al.  in nonsmall cell lung cancer (NSCLC). It is regarded as one of the highly conservative lncRNAs in mammals. The latest studies have reported MALAT1 as one of the first lncRNAs detected and was significantly upregulated under high-glucose culture in retinal endothelial cells and in the retina of diabetic mice model [9, 10]. LncRNA MALAT1 modulates inflammatory response in diabetes-initiated microvascular complications, such as DKD and diabetic retinopathy [11–13]. Therefore, this study mainly aimed to explore the role of lncRNA MALAT1 in DKD and its potential pathogenesis, which could benefit from the identification of reliable biomarkers and novel therapeutic targets.
2. Materials and Methods
This study enrolled 47 diabetic patients who sought medical advice at the Affiliated Hospital, Youjiang Medical College for Nationalities between January 2017 and December 2017, which included 20 cases with type 2 diabetes mellitus (T2DM) (12 males and 8 females, with a mean age of years) and 27 cases were with DKD (16 males and 11 females, with a mean age of years). Meanwhile, 14 nondiabetic healthy volunteers were enrolled as control group who checked into our health examination center (8 males and 6 females, with a mean age of years). No statistical difference was detected in gender and age between DM, DKD, and control groups.
DKD patients should meet the enrollment criteria of DKD diagnosis standards of China T2DM prevention and treatment guidelines 2019 . T2DM patients (1) who meet the WHO diabetes mellitus diagnosis standards 1999; and (2) without other severe diseases that affect glycemia or albuminuria, such as pregnancy, infection, contagious disease, malignancy, and drug medication. Subjects in control group were enrolled by referring to the above 2 conditions, but missing diabetes mellitus diagnostic standards by glucose tolerance test.
Patients in DM and DKD groups were given insulin glargine (Sanofi-Aventis Deutschland Gmb H) combined with acarbose (Hangzhou Huadong Medicine Co., Ltd., China) to control blood glucose. In order to reduce the experimental error, the enrolled patients should be with only DKD; those with cardiovascular and cerebrovascular diseases were excluded.
This study was reviewed and approved by our hospital ethics committee and was conducted after gaining full informed consent from patients.
2.2. Anthropometrical and Biochemical Parameter Measurements
Study Subjects. General data were collected from all study subjects, which included peripheral blood sample, urine sample, body height, body weight, blood pressure, and heart rate. Also, body mass index (BMI) was calculated using the formula: (kg/m2). Mindray Automatic Biochemical Analyzer BS-2000M (Shenzhen, China) was used for measuring biochemical indexes of superoxide dismutase (SOD), creatinine (Cr), glycated hemoglobin (HbA1c), and blood lipids. Insulin was measured with magnetic particle chemiluminescence instrument (Autobio). Specific Proteins Analyzer BA400 (Bio Systems, Spain) was used for detecting ACR, urine β2-MG, and urine α1-MG. All these parameters were measured by our technicians who received specific training. SI units were applied accordingly.
Measurement of urine β2-MG, α1-MG, microalbumin (UMA) and creatinine, and ACR calculation.
Fasting midstream urine sample of 5 mL was collected and centrifuged at 3000 r/min for 10 min to collect the supernatant. This was followed by immunoturbidimetry analysis of UMA and urine Cr using Specific Proteins Analyzer BA400. Next, ACR was calculated using the formula: . ACR readout <3 mg/g was defined as normal while ≥3 mg/g abnormal; urine β2-MG readout <0.3 mg/L was normal while ≥0.3 mg/L abnormal; urine α1-MG readout <30 mg/L was normal while ≥30 mg/L abnormal. The normal range of blood creatinine was 50~106 μmol/L.
2.3. Real-Time PCR Analysis of LncRNA MALAT1 Expression
2.3.1. Total RNA Extraction and Reverse Transcription
5 mL venous blood was drawn from all patients and added to 9 mL of red cell lysis buffer, followed by mixing before proceeding to centrifugation at 2500 rpm for 5 minutes at 4°C in an ultralow temperature centrifuge. The supernatant was then discarded, and the extracted leukocytes were rinsed with 1 mL PBS buffer before preservation at -80°C in a freezer. The total RNA of leukocytes was extracted using Trizol kit (Pufei Bio, Shanghai, China), which was then proceeded to reverse transcription according to Promega M-MLV kit (RiboBio, Guangzhou, China) manual. The reverse transcription primers were also purchased from RiboBio. All the above procedures were strictly conducted with reference to the kit manual.
2.3.2. Quantitative Real-Time PCR (qRT-PCR)
PCR analysis was performed using SYBR Green (TIANGEN Bio, Beijing, China). PCR reaction system was as follows: 95°C for 15 min, 95°C for 10s, 60°C for 20s; 40 cycles. PCR reaction was evaluated with a melting curve. Each sample was performed in triplicate, and the average value was taken for analysis. LncRNA MALAT1 expression level was reflected by 2-ΔΔCt value, which was calculated using the formula: . The expression level of GADPH genes is the best reference to evaluate relative changes in gene activity in diabetic/high glucose exposed glomerular tissues .
Primer sequences were shown in Table 1.
2.4. Statistical Methods
Variance analysis was conducted using SPSS 25.0 software and expressed as . Chi-square test was used to analyze the relations of LncRNA Malat1 with each clinicopathological parameter. Pearson’s linear regression method was applied to analyze the correlation of LncRNA Malat1 with proteinuria, glomerular filtration rate, and other clinicopathological features. In order to explore the risk factors of DKD in patients with diabetes mellitus, the indexes with difference were included in the binary regression analysis, then drew the ROC curve of each factor in identifying DKD, and jointed multifactor ROC curve to calculate the area under the ROC curve (AUC), and analyzed the diagnostic cross-section points of each index and the sensitivity and specificity of DKD identification; was considered to be statistically significant.
3.1. Comparisons of General Clinical Data in Three Groups
All the three groups showed no statistical differences in either gender, age, disease duration, BMI, ACR, urine β2-MG, or triglycerides (), while statistical significance was observed in parameters such as urine α1-MG, CR, EGFR, SOD, FBG, and HbA1C (), shown in Table 2.
indicated statistical significance, , . Abbreviations: SOD: superoxide dismutase; BMI: body mass index; HbA1c: glycated hemoglobin; FBG: fasting blood-glucose; TG: triglyceride; LDL-C: low-density lipoprotein cholesterol; Cr: creatinine; eGFR: estimated glomerular filtration rate; ACR: urine albumin/creatinine; β2-MG: urine β2-microglobulin; α1-MG: urine α1-microglobulin (α1-MG).
3.2. Differential Expression of LncRNA MALAT1 in T2DM, DKD, and Healthy Control Groups
Compared with healthy control group, LncRNA MALAT1 expression was shown to be obviously upregulated in T2DM and DKD patients, with statistical significance (, , Figure1). Further analysis uncovered the marked rise of LncRNA MALAT1 expression in DKD patients over T2DM patients, showing statistical significance (, , Figure1).
3.3. LncRNA MALAT1 Expression in Injured Renal Tubules
The results of lncRNA MALAT1 expression in renal tubule urine β2-MG revealed obvious upregulation of serum lncRNA MALAT1 in DM patients with abnormal β2-MG when compared to those with normal urine β2-MG, showing statistically significant difference (, , Figure 2(a)).
The results of lncRNA MALAT1 expression in ACR indicated upregulation of serum lncRNA MALAT1 in DM patients with abnormal ACR over those with normal ACR, showing statistically significant difference (, , Figure 2(b)).
3.4. Correlation of lncRNA MALAT1 Expression in DM Patients with SOD, ACR, Urine β2-MG, Urine α1-MG, Cr, and HbA1c
Pearson’s linear correlation coefficient was used to analyze the correlation of lncRNA MALAT1 expression in DM patients with SOD, ACR, urine β2-MG, urine α1-MG, Cr, and HbA1c. The results revealed that lncRNA MALAT1 expression was positively correlated with ACR, urine β2-MG, urine α1-MG, Cr, and HbA1c (Figures 3(a)–3(e)), while a negative correlation was observed with SOD (Figure 3(f)).
3.5. Value of Joint Indicators and Joint Prediction Probability in Identifying DKD
Binary regression analysis showed that ACR, creatinine, α1-MG, and LncRNA Malat1 (2-ΔΔCt) were the risk factors of DKD, with OR value of 1.166, 1.031, 1.031, and 2.019 (). The area under ROC curve (AUC) of DKD identified by the above indicators were 0.914, 0.643, 0.807, and 0.797, respectively. The AUC of Joint prediction probability of DKD recognition was 0.914, the sensitivity and specificity of DKD diagnosis were 1.0 and 0.806, respectively. (Take ≥0.251 as the diagnostic cutoff point), shown in Table 3.
Draw the ROC curve of each index in identifying DKD, the area under the ROC curve (AUC) was calculated by multiparameter analysis and analyzed the diagnostic cross-section points of each index and the sensitivity and specificity of DKD identification (Figure 4).
DM is a chronic metabolic disease with its complications affecting the kidney, retina, and nerve. The number of global DM patients reaches to about 366 million by 2030 and 592 million by 2035, which is 1/10 of the earth’s population [16, 17]. DKD is a common DM complication that lasts long and is poorly controlled [18–20]. In China, about 20-40% of DM patients have DKD, and it has become the primary cause of CKD and end-stage kidney diseases [3, 21]. Despite clear understanding of the molecular mechanism of DKD, increasing evidences demonstrated the critical role of genetic regulation in its pathogenesis. More and more evidences reported the key role of lncRNA in diverse diseases . Previous studies have confirmed MALAT1 dysregulation in DKD progression, and it has a value of being a diagnostic target.
This study found significant upregulation of lncRNA MALAT1 in T2DM and DKD patients when compared to healthy controls, and this was in agreement with the study results of Zhang et al. . DKD deteriorates slowly and is regarded as the main early manifestation of kidney enlargement and rise in the glomerular filtration rate (GFR), without any clinical symptoms [24–27]. Along with disease progression, GFR increase is coupled with change in the filtration membrane charges and occurrence of urine microalbumin (mAlb). But GFR still appeared to be within normal range, with no relevant emergence of clinical symptoms. The sole abnormality involves urine mAlb. So far, mAlb is an important index for laboratory diagnosis of early DKD. In clinical practice, ACR is adopted to correct the fluctuations of random urine mAlb levels. In our study, the kidney function status of DKD and DM patients was evaluated by ACR, urine β2-MG, urine α1-MG, and Cr. The relative lncRNA MALAT1 expression was increased in DM patients with high ACR, showing a significant positive correlation. Similarly, the relative lncRNA MALAT1 expression was also increased in DM patients with high urine β2-MG, showing a significant positive correlation. In this study, the lncRNA MALAT1 expression showed a positive correlation with urine α1-MG, Cr, and HbA1c, implying lncRNA MALAT1 involvement in kidney impairment. The study conducted by Li et al. reported the role of MALAT1, miR-23c, and its target gene pyrophosphorylation-related protein ELAVL1 in renal tubular epithelial cells, and these played a potential role by affecting pyrophosphorylation cell death signaling pathway in DKD. LncRNA MALAT1 expression was upregulated in hyperglycemic conditions. Luciferase assay confirmed MALAT1 as the target gene of miR-23c. Therefore, downregulation of cellular miR-23c increased the expression of pyrophosphorylation-related protein ELAVL1 and NLRP3, followed by Caspase-1 activation, and subsequent secretion of proinflammatory chemokines IL-1β and IL-18, thereby promoting inflammatory response. In contrast, MALAT1 downregulation or miR-23c upregulation elicited suppression of ELAVL1 expression, which inhibited pyrophosphorylation signaling pathway, thereby modulating DKD initiation and development.
Our study also uncovered the negative correlation of lncRNA MALAT1 expression with SOD. Bulky production of activated oxygen and overwhelming oxidative stress response were considered as important pathological changes during DKD development . SOD, an antioxidative enzyme, catalyzes disproportionation of superoxide anion radical to hydrogen peroxide and oxygen, maintaining the balance of free radicals in the body . Previous reports demonstrated decreased SOD level in DM patients, which increased the peroxides in the glomerular microvascular system, which subsequently escalated the cytotoxic effect caused by peroxides, thereby resulting in DM renal cell damage . SOD level was negatively correlated with the severity of renal function impairment, indicating the importance of SOD during various stages of DKD development. SOD is regarded as a sensitive index for dynamic monitoring of DKD development.
Binary regression analysis showed that ACR, creatinine, α1-MG, and LncRNA Malat1 (2-ΔΔCt) were the risk factors of DKD, the sensitivity and specificity of DKD can be further improved to 100% and 80.6% by using the combined predictor variables, which is higher than that of single index. This study initially reported high expression of lncRNA MALAT1 in the circulating of DKD patients and showed correlation with ACR, urine β2-MG, urine α1-MG, and Cr. This study would further benefit from the theories of DKD pathogenesis and provide novel targets for DKD gene-targeted therapy. Also, our findings provide references to seek out biomarkers for diagnosis and outcome prediction of DKD and offer a novel strategy for managing clinically intractable DKD. Nonetheless, the influential factors of DKD are diverse. The specific pathogenesis, diagnosis, and outcome assessment should be further elucidated.
The inadequacies of this article are (1) insufficient sample size. (2) Long time span of sample collection, causing some RNA to degrade, Malat1 has decreased in the diabetic nephropathy group. (3) Since other nephropathy except diabetic nephropathies has been excluded in this study, this model is only applicable to the patients with diabetes mellitus history and who have been evaluated and considered as diabetic nephropathies.
The data for the current study are available from the corresponding author upon reasonable request.
All procedures performed in this study involving human participants were in accordance with the ethical standards of the institutional and national research committee and with the World Medical Association Declaration of Helsinki or comparable ethical standards.
Conflicts of Interest
The authors have no conflicts of interest to declare.
Lian-ji Zhou and Da-wei Yang contributed equally to this work.
This work was supported by the National Natural Science Foundation of China (grant number 81660797) and Guangxi Natural Science Foundation (2018GXNSFAA281124). We appreciate all of the patients for their agreement to participate in our study.
- H. J. L. Heerspink, H. H. Parving, D. L. Andress et al., “Atrasentan and renal events in patients with type 2 diabetes and chronic kidney disease (SONAR): a double-blind, randomised, placebo-controlled trial,” Lancet, vol. 393, no. 10184, pp. 1937–1947, 2019.
- S. Bell, E. H. Fletcher, I. Brady et al., “End-stage renal disease and survival in people with diabetes: a national database linkage study,” QJM, vol. 108, no. 2, pp. 127–134, 2015.
- L. Zhang, J. Long, W. Jiang et al., “Trends in chronic kidney disease in China,” The New England Journal of Medicine, vol. 375, no. 9, pp. 905-906, 2016.
- N. H. Cho, J. E. Shaw, S. Karuranga et al., “IDF diabetes atlas: global estimates of diabetes prevalence for 2017 and projections for 2045,” Diabetes Research and Clinical Practice, vol. 138, pp. 271–281, 2018.
- E. Ritz, I. Rychlík, F. Locatelli, and S. Halimi, “End-stage renal failure in type 2 diabetes: a medical catastrophe of worldwide dimensions,” American Journal of Kidney Diseases, vol. 34, no. 5, pp. 795–808, 1999.
- C. H. Li and Y. Chen, “Targeting long non-coding RNAs in cancers: progress and prospects,” International Journal of Biochemistry and Cell Biology, vol. 45, no. 8, pp. 1895–1910, 2013.
- H. Yi, R. Peng, L.-y. Zhang et al., “LincRNA-Gm4419 knockdown ameliorates NF-κB/NLRP3 inflammasome-mediated inflammation in diabetic nephropathy,” Cell Death & Disease, vol. 8, no. 2, p. e2583, 2017.
- P. Ji, S. Diederichs, W. Wang et al., “MALAT-1, a novel noncoding RNA, and thymosin _β_ 4 predict metastasis and survival in early-stage non-small cell lung cancer,” Oncogene, vol. 22, no. 39, pp. 8031–8041, 2003.
- T. Gutschner, M. Hammerle, and S. Diederichs, “MALAT1 -- a paradigm for long noncoding RNA function in cancer,” Journal of Molecular Medicine (Berlin, Germany), vol. 91, no. 7, pp. 791–801, 2013.
- R. Yoshimoto, A. Mayeda, M. Yoshida, and S. Nakagawa, “MALAT1 long non-coding RNA in cancer,” Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms, vol. 1859, no. 1, pp. 192–199, 2016.
- P. Puthanveetil, S. Chen, B. Feng, A. Gautam, and S. Chakrabarti, “Long non-coding RNA MALAT1 regulates hyperglycaemia induced inflammatory process in the endothelial cells,” Journal of Cellular and Molecular Medicine, vol. 19, no. 6, pp. 1418–1425, 2015.
- A. Thomas, A. Gordon, B. Feng, and S. Chakrabarti, “Alterations of Long Noncoding RNAs (lncRNA) Cause Pathogenetic Changes in Diabetic RetinopathyImage 4,” Canadian Journal of Diabetes, vol. 40, no. 5, p. S52, 2016.
- L. Zhou, D. Y. Xu, W. G. Sha, L. Shen, G. Y. Lu, and X. Yin, “Long non-coding MIAT mediates high glucose-induced renal tubular epithelial injury,” Biochemical and Biophysical Research Communications, vol. 468, no. 4, pp. 726–732, 2015.
- W. Jia, J. Weng, D. Zhu et al., “Standards of medical care for type 2 diabetes in China 2019,” Diabetes/Metabolism Research and Reviews, vol. 35, no. 6, article e3158, 2019.
- J. Biederman, J. Yee, and P. Cortes, “Validation of internal control genes for gene expression analysis in diabetic glomerulosclerosis,” Kidney International, vol. 66, no. 6, pp. 2308–2314, 2004.
- V. V. S. Uddandrao, P. Brahmanaidu, R. Ravindarnaik, P. Suresh, S. Vadivukkarasi, and G. Saravanan, “Restorative potentiality of S-allylcysteine against diabetic nephropathy through attenuation of oxidative stress and inflammation in streptozotocin-nicotinamide-induced diabetic rats,” European Journal of Nutrition, vol. 58, no. 6, 2019.
- F. Wagnew, S. Eshetie, G. D. Kibret et al., “Diabetic nephropathy and hypertension in diabetes patients of sub-Saharan countries: a systematic review and meta-analysis,” BMC Research Notes, vol. 11, no. 1, p. 565, 2018.
- V. Harjutsalo and P. H. Groop, “Epidemiology and risk factors for diabetic kidney disease,” Advances in Chronic Kidney Disease, vol. 21, no. 3, pp. 260–266, 2014.
- C. W. Park, “Diabetic kidney disease: from epidemiology to clinical perspectives,” Diabetes and Metabolism Journal, vol. 38, no. 4, pp. 252–260, 2014.
- L. Li, S. Jick, S. Breitenstein, and A. Michel, “Prevalence of diabetes and diabetic nephropathy in a large U.S. commercially insured pediatric population, 2002-2013,” Diabetes Care, vol. 39, no. 2, pp. 278–284, 2016.
- American Diabetes Association, “10. Microvascular complications and foot care,” Diabetes Care, vol. 40, Suppl 1, pp. S88–S98, 2016.
- E. Arvaniti, P. Moulos, A. Vakrakou et al., “Whole-transcriptome analysis of UUO mouse model of renal fibrosis reveals new molecular players in kidney diseases,” Scientific Reports, vol. 6, no. 1, p. 26235, 2016.
- Y. Zhang, H. Wu, F. Wang, M. Ye, H. Zhu, and S. Bu, “Long non-coding RNA MALAT1 expression in patients with gestational diabetes mellitus,” International Journal of Gynaecology and Obstetrics, vol. 140, no. 2, pp. 164–169, 2018.
- C. Qi, X. Mao, Z. Zhang, and H. Wu, “Classification and differential diagnosis of diabetic nephropathy,” Journal Diabetes Research, vol. 2017, article 8637138, pp. 1–7, 2017.
- N. Papadopoulou-Marketou, G. P. Chrousos, and C. Kanaka-Gantenbein, “Diabetic nephropathy in type 1 diabetes: a review of early natural history, pathogenesis, and diagnosis,” Diabetes/Metabolism Research and Reviews, vol. 33, no. 2, 2017.
- M. Bose, S. Almas, and S. Prabhakar, “Wnt signaling and podocyte dysfunction in diabetic nephropathy,” Journal of Investigative Medicine, vol. 65, no. 8, pp. 1093–1101, 2017.
- A. Flyvbjerg, “The role of the complement system in diabetic nephropathy,” Nature Reviews Nephrology, vol. 13, no. 5, pp. 311–318, 2017.
- X. Li, L. Zeng, C. Cao et al., “Long noncoding RNA MALAT1 regulates renal tubular epithelial pyroptosis by modulated miR-23c targeting of ELAVL1 in diabetic nephropathy,” Experimental Cell Research, vol. 350, no. 2, pp. 327–335, 2017.
- J. Sha, B. Sui, X. Su, Q. Meng, and C. Zhang, “Alteration of oxidative stress and inflammatory cytokines induces apoptosis in diabetic nephropathy,” Molecular Medicine Reports, vol. 16, no. 5, pp. 7715–7723, 2017.
- K. Asaba, A. Tojo, M. L. Onozato, A. Goto, and T. Fujita, “Double-edged action of SOD mimetic in diabetic nephropathy,” Journal of Cardiovascular Pharmacology, vol. 49, no. 1, pp. 13–19, 2007.
- H. Fujita, H. Fujishima, S. Chida et al., “Reduction of renal superoxide dismutase in progressive diabetic nephropathy,” Journal of the American Society of Nephrology, vol. 20, no. 6, pp. 1303–1313, 2009.
Copyright © 2020 Lian-ji Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.