- Research article
- Open Access
- Open Peer Review
Diabetes-specific genetic effects on obesity traits in American Indian populations: the Strong Heart Family Study
© Franceschini et al; licensee BioMed Central Ltd. 2008
- Received: 08 January 2008
- Accepted: 14 October 2008
- Published: 14 October 2008
Body fat mass distribution and deposition are determined by multiple environmental and genetic factors. Obesity is associated with insulin resistance, hyperinsulinemia, and type 2 diabetes. We previously identified evidence for genotype-by-diabetes interaction on obesity traits in Strong Heart Family Study (SHFS) participants. To localize these genetic effects, we conducted genome-wide linkage scans of obesity traits in individuals with and without type 2 diabetes, and in the combined sample while modeling interaction with diabetes using maximum likelihood methods (SOLAR 2.1.4).
SHFS recruited American Indians from Arizona, North and South Dakota, and Oklahoma. Anthropometric measures and diabetes status were obtained during a clinic visit. Marker allele frequencies were derived using maximum likelihood methods estimated from all individuals and multipoint identity by descent sharing was estimated using Loki. We used variance component linkage analysis to localize quantitative trait loci (QTLs) influencing obesity traits. We tested for evidence of additive and QTL-specific genotype-by-diabetes interactions using the regions identified in the diabetes-stratified analyses.
Among 245 diabetic and 704 non-diabetic American Indian individuals, we detected significant additive gene-by-diabetes interaction for weight and BMI (P < 0.02). In analysis accounting for QTL-specific interaction (P < 0.001), we detected a QTL for weight on chromosome 1 at 242 cM (LOD = 3.7). This chromosome region harbors the adiponectin receptor 1 gene, which has been previously associated with obesity.
These results suggest distinct genetic effects on body mass in individuals with diabetes compared to those without diabetes, and a possible role for one or more genes on chromosome 1 in the pathogenesis of obesity.
- Marker Allele Frequency
- American Indian Population
- Genetic Standard Deviation
- Obesity Trait
- ADIPOR1 Gene
Body fat mass distribution and deposition are determined by multiple environmental and genetic factors. Obesity is associated with insulin resistance, hyperinsulinemia, and incident type 2 diabetes [1, 2]. Moreover, insulin has anabolic effects on fat metabolism leading to fat deposition and obesity . Obesity and diabetes are highly prevalent in American Indian populations. Among individuals between 45 and 74 years who were recruited into the Strong Heart Study, approximately 35% were obese  and type 2 diabetes prevalence varied by recruiting center from 33% in South and North Dakota men to 72% in Arizona women . The estimated age-adjusted rates of diabetes were higher in recruiting centers with a higher prevalence of obesity .
Diabetes and obesity are associated with increased morbidity and mortality in the general population  and among American Indian subjects, and they account for increased incident hypertension , and fatal and nonfatal cardiovascular disease (CVD) events . The identification of genes for obesity in the environmental and genetic milieu of type 2 diabetes may bring us one step closer to the identification of the functional genes that influence susceptibility to obesity, a major public health problem among American Indian populations.
The genetic component of obesity-related traits among American Indians has been previously described [9–11]. In our previous work , we identified evidence for distinct genetic effects on obesity related traits in diabetic compared to non-diabetic Strong Heart Family Study (SHFS) participants. A significant additive genotype-by-diabetes interaction was observed for body mass index (BMI), waist-to-hip ratio (WHR) and percent body fat. The purpose of the present study is to localize these genetic effects on obesity-related traits while accounting for interaction with diabetes status.
The SHFS recruitment and genotyping methods have been described elsewhere [12, 13]. Briefly, American Indian subjects from 13 tribes were recruited from centers located in Arizona, Oklahoma and North and South Dakota. Clinical information and laboratory measures were obtained during a study exam. Anthropometric measures of body weight (kg) and height (m) were used to estimate BMI (kg/m2). WHR (in cm) was also recorded. Percentage body fat was measured using Quantum II bioelectric impedance (RJL Systems, Clinton Township, MI) . The percentage of body fat was calculated using a formula based on total body water . Type 2 diabetes was defined as a fasting blood glucose of 126 mg/dl or higher, history of diabetes or use of diabetic medications . The SHFS protocols were approved by the Indian Health Service Institutional Review Board (IRB), the IRBs from participating institutions and by the 13 tribes participating in these studies. The study was conducted in accordance with the principles of the Declaration of Helsinki. Written consent was obtained from all individuals of this study.
Participants, who are members of extended families, were genotyped for approximately 400 microsatellite markers spaced at average 10 cM intervals across the autosomal chromosomes. Pedigree relationships were verified using the Pedigree Relationship Statistical Tests program , which employs likelihood-based inference statistics for genome-wide identity-by-descent (IBD) allele sharing. Mendelian inconsistencies and spurious double recombinants were detected using the SimWalk2 program . The overall blanking rate for both types of errors was less than 1% of the total number of genotypes for Arizona, Dakota and Oklahoma. Marker allele frequencies were estimated using maximum likelihood methods estimated from all individuals  and multipoint IBD sharing was estimated using Loki software .
We conducted a genome scan of obesity traits while modeling interaction with diabetes status using maximum likelihood methods implemented in SOLAR 2.1.4. We studied 976 individuals from extended families enrolled between 1998 and 1999 (SHS phase 3). We examined the distribution of the obesity-related quantitative traits of body weight, BMI, WHR and percentage of body fat, and excluded observations greater than 2.5 standard deviations above or below the mean (< 5 observations excluded). The traits were adjusted for the effects of age, age2, sex, age-by-sex interaction and study center using linear regression models and an alpha = 0.05. The WHR trait was scaled for analysis by multiplying by 10.
We examined the evidence for additive genotype-by-diabetes interaction on obesity traits using maximum likelihood variance decomposition methods extended to allow for the genetic covariance between diabetic and non-diabetic subjects [19, 20]. We tested the hypothesis of interaction using the likelihood ratio test comparing a model including genotype-by-diabetes interaction to restricted models in which such interaction is excluded. We estimated the genetic correlation between diabetic and non-diabetic individuals (ρG) and the genetic standard deviations for diabetic and non-diabetic subjects. To test for differential additive genetic effects among diabetic and non-diabetic subjects, we compared a model in which ρG (diabetic, non-diabetic) =1 to a model without constraints. To test for differences in the magnitude of the genetic effects among diabetic and non-diabetic subjects, we compared a model in which the genetic standard deviations of the two groups are constrained to be equal to a model without these constraints. We considered evidence for genotype-by-environment interaction if the genetic correlation between diabetics and non-diabetics was different than 1 and/or the genetic variance for obesity traits in diabetics was significantly different than in non-diabetics.
We also tested for QTL-specific diabetes-by-obesity interaction. For these analyses, we compared a model in which the QTL standard deviations are constrained to be equal, to a model without constraints, using the likelihood ratio statistics. When comparing models with variances constrained to be equal, interpretation of significant differences were based on the assumption of an asymptotic chi-square distribution for the likelihood test statistic. For models that restricted the genetic correlation to one (the upper boundary of the parameter space, i.e., ρG = 1.0), the test statistic is asymptotically distributed as a 50:50 mixture of a chi-square and a point mass at 0 . The gene-by-environment analysis has an additional QTL variance component compared with the standard linkage model, therefore we corrected for the increase in degrees of freedom and all reported LOD scores are corrected LOD score .
For the obesity traits showing significant additive genotype-by-diabetes interaction at alpha = 0.05 (ρG ≠ 1), we performed genome scans accounting for the diabetes interaction. We adjusted the obtained LOD scores for the difference in the number of parameters between the polygenic and linkage versions of the interaction model, taking into account that each of the QTL standard deviations has a boundary at the null hypothesis.
Descriptive characteristics of American Indian subjects by diabetes status
Mean age (SD)
Mean Weight (kg) (SD)
Mean BMI (kg/m2) (SD)
Mean % body fat (SD)
Mean WHR (SD)
Number of relative-pairs
Additive genotype-by-diabetes interaction for obesity related quantitative traits in American Indian participants of the Strong Heart Family Study
Genetic standard deviation
Variance due to covariates*
Residual kurtosis after covariates
Genetic correlation (ρG) ± SE among diabetics and nondiabetics (P-value)
0.56 ± 0.16 (0.01)
0.45 ± 0.17 (0.0002)
0.67 ± 0.18 (0.08)
0.72 ± 0.15 (0.02)
We then evaluated the evidence for genotype-by-diabetes interaction at the QTL level, using the combined samples of diabetic and non-diabetic individuals and the maximum LOD score QTL identified on the diabetes-stratified analysis. The QTL- specific genotype-by-diabetes interaction was significant for weight on chromosome 1 at 242 cM (p for interaction < 0.001).
This study identified a diabetes-specific QTL for body weight on chromosome 1 among American Indians, which was less prominent for BMI and for the fat distribution trait of WHR. Diabetes and obesity are strongly associated in American Indian and other populations, and environmental factors such as diet and lifestyle play an important role in both conditions. Our findings suggest a role for diabetes-specific genes influencing body fat deposition in American Indians. Further studies are required to explore the generalizability of these findings to populations of different ethnic/racial and cultural backgrounds, in which the prevalence of obesity and diabetes vary from American Indians. Nevertheless, our findings suggest that the genetic susceptibility to obesity and type 2 diabetes may share common genetic pathways, by complex interactions of genes influencing obesity in the milieu of diabetes.
Several prior studies have showed linkage to the chromosome 1 region for diabetes-related traits and obesity , although genetic interactions have not been explored at this QTL. The 1q32 region was identified in linkage analysis of type 2 diabetes in African American participants of the Hypertension Genetic Epidemiology Network (HyperGEN) (LOD = 1.9) . Among non-diabetic individuals, linkage to this QTL was described in studies of fasting glucose among Pima Indians (LOD = 1.5)  and among participants of the Framingham Offspring Study (LOD = 1.8) , and in studies of fasting insulin among HyperGEN males (LOD = 2.1) . In addition, a West African study of type 2 diabetic sib-pairs described a LOD score of 1.6 for BMI and 1.5 for fat mass (estimated by bioelectric impedance), between the 1q25.1-q31.1 regions, which are within the 1-LOD score support region of our linkage peak . Participants in this West African study were about 5 BMI units heavier than the population from which they were enrolled .
Within the 1-LOD support region lays the ADIPOR1 gene. The adiponectin 1 receptor mediates adiponectin effects on fatty acid oxidation and glucose uptake through activation of AMP-activated protein kinase and peroxisome proliferator-activated receptor-α . ADIPOR1 deficient (-/-) mice show increased adiposity and decreased glucose tolerance, associated with decreased locomotor activity and energy expenditure . ADIPOR1 variants have been associated with type 2 diabetes, pre-diabetes states and measures of obesity although not in all studies [28–33]. For example, among Caucasian obese individuals (mean BMI 32 kg/m2) with impaired glucose tolerance, ADIPOR1 SNPs were associated with measures of body size including weight, waist, hip circumference and BMI . In addition, a gene-gene interaction of ADIPOR1 and ADIPONECTIN was observed in increased overall and abdominal adiposity and higher substrate oxidation . Therefore, ADIPOR1 is a potential candidate gene for pleiotropic effects on obesity and diabetes susceptibility in humans.
It is difficult to gauge the LOD score evidence reported herein in the face of the number of comparisons made in conducting gene-by-diabetes genome scans. In total, we conducted three genome scans (a minimum model, a maximal model and an interaction model). If we were to use a too conservative Bonferoni correction, we would need to observe a P = 0.05/3 * 380 markers: P = 0.00004 or a corresponding LOD score of 3.7. However, given the correlation and non-independence of the three genome scans of each obesity trait, a systematic correction would be extremely difficult.
Multiple interacting genes may underlie this QTL on chromosome 1, and therefore, densely genotyping this region is the next logical step to be taken to identify these genes influencing obesity in the context of diabetes. Our findings are limited to American Indian populations and should be validated in replication samples of other populations.
In summary, we identified a diabetes-specific QTL on chromosome 1 at 242 cM for weight in American Indians, in the region where the ADIPOR1 gene is located. This gene has been associated with multiple measures of body size and diabetes-related traits in humans and experimental models and, therefore, it is likely a good candidate for susceptibility to fat deposition in diabetic individuals. Importantly, our findings support that this region of the genome likely harbors susceptibility variants for obesity that are importantly different in type 2 diabetes individuals and those with normal glucose metabolic profile. Thus, future studies should more comprehensively evaluate the ADIPOR1 gene, especially with respect to the complex relationship between body weight and type 2 diabetes.
We thank the Strong Heart Family Study participants, the Indian Health Service hospitals and clinics staff, the directors of the SHS clinics and the many collaborators and staff of the Strong Heart Study. The opinions expressed in this paper are those of the author(s) and do not necessarily reflect the views of the Indian Health Service.
This research was funded by a cooperative agreement that includes grants U01 HL65520, U01 HL41642, U01 HL41652, U01 HL41654, and U01 HL65521 from the National Heart, Lung, and Blood Institute. Development of SOLAR and the methods implemented in it are supported by US Public Health Service grant MH059490 from the National Institutes of Health. Dr Franceschini is supported by the AHA #0675001N grant.
- Ford ES, Williamson DF, Liu S: Weight change and diabetes incidence: findings from a national cohort of US adults. Am J Epidemiol. 1997, 146 (3): 214-222.View ArticlePubMedGoogle Scholar
- Resnick HE, Valsania P, Halter JB, Lin X: Relation of weight gain and weight loss on subsequent diabetes risk in overweight adults. Journal of epidemiology and community health. 2000, 54 (8): 596-602. 10.1136/jech.54.8.596.View ArticlePubMedPubMed CentralGoogle Scholar
- Unger RH: Lipid overload and overflow: metabolic trauma and the metabolic syndrome. Trends in endocrinology and metabolism: TEM. 2003, 14 (9): 398-403. 10.1016/j.tem.2003.09.008.View ArticlePubMedGoogle Scholar
- Welty TK, Lee ET, Yeh J, Cowan LD, Go O, Fabsitz RR, Le NA, Oopik AJ, Robbins DC, Howard BV: Cardiovascular disease risk factors among American Indians. The Strong Heart Study. Am J Epidemiol. 1995, 142 (3): 269-287.PubMedGoogle Scholar
- Lee ET, Howard BV, Savage PJ, Cowan LD, Fabsitz RR, Oopik AJ, Yeh J, Go O, Robbins DC, Welty TK: Diabetes and impaired glucose tolerance in three American Indian populations aged 45–74 years. The Strong Heart Study. Diabetes Care. 1995, 18 (5): 599-610. 10.2337/diacare.18.5.599.View ArticlePubMedGoogle Scholar
- Must A, Spadano J, Coakley EH, Field AE, Colditz G, Dietz WH: The disease burden associated with overweight and obesity. Jama. 1999, 282 (16): 1523-1529. 10.1001/jama.282.16.1523.View ArticlePubMedGoogle Scholar
- Wang W, Lee ET, Fabsitz RR, Devereux R, Best L, Welty TK, Howard BV: A longitudinal study of hypertension risk factors and their relation to cardiovascular disease: the Strong Heart Study. Hypertension. 2006, 47 (3): 403-409. 10.1161/01.HYP.0000200710.29498.80.View ArticlePubMedGoogle Scholar
- Howard BV, Lee ET, Cowan LD, Devereux RB, Galloway JM, Go OT, Howard WJ, Rhoades ER, Robbins DC, Sievers ML, et al: Rising tide of cardiovascular disease in American Indians. The Strong Heart Study. Circulation. 1999, 99 (18): 2389-2395.View ArticlePubMedGoogle Scholar
- Almasy L, Goring HH, Diego V, Cole S, Laston S, Dyke B, Howard BV, Lee ET, Best LG, Devereux R, et al: A novel obesity locus on chromosome 4q: the Strong Heart Family Study. Obesity (Silver Spring). 2007, 15 (7): 1741-1748. 10.1038/oby.2007.207.View ArticleGoogle Scholar
- Diego VP, Goring HH, Cole SA, Almasy L, Dyer TD, Blangero J, Duggirala R, Laston S, Wenger C, Cantu T, et al: Fasting insulin and obesity-related phenotypes are linked to chromosome 2p: the strong heart family study. Diabetes. 2006, 55 (6): 1874-1878. 10.2337/db05-0668.View ArticlePubMedGoogle Scholar
- North KE, MacCluer JW, Williams JT, Welty TK, Best LG, Lee ET, Fabsitz RR, Howard BV: Evidence for distinct genetic effects on obesity and lipid-related CVD risk factors in diabetic compared to nondiabetic American Indians: the Strong Heart Family Study. Diabetes Metab Res Rev. 2003, 19 (2): 140-147. 10.1002/dmrr.355.View ArticlePubMedGoogle Scholar
- North KE, Howard BV, Welty TK, Best LG, Lee ET, Yeh JL, Fabsitz RR, Roman MJ, MacCluer JW: Genetic and environmental contributions to cardiovascular disease risk in American Indians: the strong heart family study. Am J Epidemiol. 2003, 157 (4): 303-314. 10.1093/aje/kwf208.View ArticlePubMedGoogle Scholar
- Lee ET, Welty TK, Fabsitz R, Cowan LD, Le NA, Oopik AJ, Cucchiara AJ, Savage PJ, Howard BV: The Strong Heart Study. A study of cardiovascular disease in American Indians: design and methods. Am J Epidemiol. 1990, 132 (6): 1141-1155.PubMedGoogle Scholar
- Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Diabetes Care. 1997, 20 (7): 1183-1197.Google Scholar
- Sun L, Wilder K, McPeek MS: Enhanced pedigree error detection. Hum Hered. 2002, 54 (2): 99-110. 10.1159/000067666.View ArticlePubMedGoogle Scholar
- Sobel E, Papp JC, Lange K: Detection and integration of genotyping errors in statistical genetics. Am J Hum Genet. 2002, 70 (2): 496-508. 10.1086/338920.View ArticlePubMedPubMed CentralGoogle Scholar
- Boehnke M: Allele frequency estimation from data on relatives. Am J Hum Genet. 1991, 48 (1): 22-25.PubMedPubMed CentralGoogle Scholar
- Heath SC, Snow GL, Thompson EA, Tseng C, Wijsman EM: MCMC segregation and linkage analysis. Genet Epidemiol. 1997, 14 (6): 1011-1016. 10.1002/(SICI)1098-2272(1997)14:6<1011::AID-GEPI75>3.0.CO;2-L.View ArticlePubMedGoogle Scholar
- Williams JT, Van Eerdewegh P, Almasy L, Blangero J: Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. I. Likelihood formulation and simulation results. Am J Hum Genet. 1999, 65 (4): 1134-1147. 10.1086/302570.View ArticlePubMedPubMed CentralGoogle Scholar
- Towne B, Siervogel RM, Blangero J: Effects of genotype-by-sex interaction on quantitative trait linkage analysis. Genet Epidemiol. 1997, 14 (6): 1053-1058. 10.1002/(SICI)1098-2272(1997)14:6<1053::AID-GEPI82>3.0.CO;2-G.View ArticlePubMedGoogle Scholar
- Blangero J, MacCluer JW, Kammerer CM, Mott GE, Dyer TD, McGill HC: Genetic analysis of apolipoprotein A-I in two dietary environments. Am J Hum Genet. 1990, 47 (3): 414-428.PubMedPubMed CentralGoogle Scholar
- Chen G, Adeyemo AA, Johnson T, Zhou J, Amoah A, Owusu S, Acheampong J, Agyenim-Boateng K, Eghan BA, Oli J, et al: A genome-wide scan for quantitative trait loci linked to obesity phenotypes among West Africans. International journal of obesity (2005). 2005, 29 (3): 255-259. 10.1038/sj.ijo.0802873.View ArticleGoogle Scholar
- Avery CL, Freedman BI, Heiss G, Kraja A, Rice T, Arnett D, Miller MB, Pankow JS, Lewis CE, Myers RH, et al: Linkage analysis of diabetes status among hypertensive families: the Hypertension Genetic Epidemiology Network study. Diabetes. 2004, 53 (12): 3307-3312. 10.2337/diabetes.53.12.3307.View ArticlePubMedGoogle Scholar
- Pratley RE, Thompson DB, Prochazka M, Baier L, Mott D, Ravussin E, Sakul H, Ehm MG, Burns DK, Foroud T, et al: An autosomal genomic scan for loci linked to prediabetic phenotypes in Pima Indians. J Clin Invest. 1998, 101 (8): 1757-1764. 10.1172/JCI1850.View ArticlePubMedPubMed CentralGoogle Scholar
- Meigs JB, Panhuysen CI, Myers RH, Wilson PW, Cupples LA: A genome-wide scan for loci linked to plasma levels of glucose and HbA(1c) in a community-based sample of Caucasian pedigrees: The Framingham Offspring Study. Diabetes. 2002, 51 (3): 833-840. 10.2337/diabetes.51.3.833.View ArticlePubMedGoogle Scholar
- North KE, Franceschini N, Borecki IB, Gu CC, Heiss G, Province MA, Arnett DK, Lewis CE, Miller MB, Myers RH, et al: Genotype-by-sex interaction on fasting insulin concentration: the HyperGEN study. Diabetes. 2007, 56 (1): 137-142. 10.2337/db06-0624.View ArticlePubMedGoogle Scholar
- Kadowaki T, Yamauchi T, Kubota N, Hara K, Ueki K, Tobe K: Adiponectin and adiponectin receptors in insulin resistance, diabetes, and the metabolic syndrome. J Clin Invest. 2006, 116 (7): 1784-1792. 10.1172/JCI29126.View ArticlePubMedPubMed CentralGoogle Scholar
- Damcott CM, Ott SH, Pollin TI, Reinhart LJ, Wang J, O'Connell JR, Mitchell BD, Shuldiner AR: Genetic variation in adiponectin receptor 1 and adiponectin receptor 2 is associated with type 2 diabetes in the Old Order Amish. Diabetes. 2005, 54 (7): 2245-2250. 10.2337/diabetes.54.7.2245.View ArticlePubMedGoogle Scholar
- Stefan N, Machicao F, Staiger H, Machann J, Schick F, Tschritter O, Spieth C, Weigert C, Fritsche A, Stumvoll M, et al: Polymorphisms in the gene encoding adiponectin receptor 1 are associated with insulin resistance and high liver fat. Diabetologia. 2005, 48 (11): 2282-2291. 10.1007/s00125-005-1948-3.View ArticlePubMedGoogle Scholar
- Siitonen N, Pulkkinen L, Mager U, Lindstrom J, Eriksson JG, Valle TT, Hamalainen H, Ilanne-Parikka P, Keinanen-Kiukaanniemi S, Tuomilehto J, et al: Association of sequence variations in the gene encoding adiponectin receptor 1 (ADIPOR1) with body size and insulin levels. The Finnish Diabetes Prevention Study. Diabetologia. 2006, 49 (8): 1795-1805. 10.1007/s00125-006-0291-7.View ArticlePubMedGoogle Scholar
- Hara K, Horikoshi M, Kitazato H, Yamauchi T, Ito C, Noda M, Ohashi J, Froguel P, Tokunaga K, Nagai R, et al: Absence of an association between the polymorphisms in the genes encoding adiponectin receptors and type 2 diabetes. Diabetologia. 2005, 48 (7): 1307-1314. 10.1007/s00125-005-1806-3.View ArticlePubMedGoogle Scholar
- Wang H, Zhang H, Jia Y, Zhang Z, Craig R, Wang X, Elbein SC: Adiponectin receptor 1 gene (ADIPOR1) as a candidate for type 2 diabetes and insulin resistance. Diabetes. 2004, 53 (8): 2132-2136. 10.2337/diabetes.53.8.2132.View ArticlePubMedGoogle Scholar
- Vaxillaire M, Dechaume A, Vasseur-Delannoy V, Lahmidi S, Vatin V, Lepretre F, Boutin P, Hercberg S, Charpentier G, Dina C, et al: Genetic analysis of ADIPOR1 and ADIPOR2 candidate polymorphisms for type 2 diabetes in the Caucasian population. Diabetes. 2006, 55 (3): 856-861. 10.2337/diabetes.55.03.06.db05-0665.View ArticlePubMedGoogle Scholar
- Loos RJ, Ruchat S, Rankinen T, Tremblay A, Perusse L, Bouchard C: Adiponectin and adiponectin receptor gene variants in relation to resting metabolic rate, respiratory quotient, and adiposity-related phenotypes in the Quebec Family Study. Am J Clin Nutr. 2007, 85 (1): 26-34.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/9/90/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.