Diabetes 53:3274-3285, 2004 © 2004 by the American Diabetes Association, Inc. Genetic Determinants of Energy Expenditure and Insulin Resistance in Diet-Induced Obesity in MiceFrom the Research Division, Joslin Diabetes Center and Department of Medicine, Harvard Medical School, Boston, Massachusetts
Diet-induced obesity is the primary determinant of the current epidemic of diabetes. We have explored the role of genetics in this phenomenon, using C57Bl/6 (B6), 129S6/SvEvTac (129), and intercross (B6 x 129)F2 mice on a low- or high-fat diet. Over an 18-week period, B6 and F2 mice gained more weight, had higher levels of insulin and leptin, and showed greater glucose intolerance than 129 mice, despite lower food intake. By contrast, metabolic rate and diet-induced thermogenesis were significantly higher in the 129 mice. Genome-wide scans identified several quantitative trait loci, including a quantitative trait locus that was linked with hyperinsulinemia/insulin resistance on chromosome 14 in a region similar to that seen in mice with genetically induced insulin resistance. Microarray analysis indicated significant changes in expression levels between B6 and 129 mice in the identified chromosomal area of Wnt5a and protein kinase C (PKC ). Thus, caloric efficiency, i.e., the "thrifty gene," is a dominant-acting genetic determinant of diet-induced obesity in mice and can be linked to a locus on chromosome 14, including genes linked to adipose development and insulin sensitivity.
There is a worldwide epidemic of obesity and type 2 diabetes (1). This is being driven, at least in part, by the westernization of diet with high-fat and higher caloric intake (2). The impact of these changes in dietary pattern, however, varies dramatically depending on the genetic background of the population exposed, with higher rates of diabetes and obesity in Hispanic, black, and Native-American groups. Defining the background genes that modify this response is central to understanding the pathogenesis of these conditions. Background genes play a major role in determining the phenotype of mice that are subjected to different diets (3,4), mice with naturally occurring mutations leading to obesity and diabetes (5,6), and mice that have been genetically altered by introduction of transgenes or knockout of endogenous genes (7,8). This complicates both the interpretation of genotype-phenotype relationships and the understanding of what is due to the primary genetic alteration and what is due to background genes or environmental modifiers (7,9). Several studies have suggested that C57Bl/6 (B6) mice are "obesity prone" as well as "susceptible to insulin resistance and glucose intolerance," whereas other strains, such as C3H/He, 129/Sv, and A/J mice, are resistant to obesity and diabetes (3,4,1013). We have previously shown that mice with a double-heterozygous deletion of the insulin receptor and insulin receptor substrate-1 (IRS-1) become markedly insulin resistant and hyperinsulinemic on the B6 background and that >90% of the mice develop diabetes by 6 months of age, whereas on the background of 129S6/SvEvTac (129), <2% of the double heterozygote mice become diabetic by 6 months, and the mice develop only a mild hyperinsulinemia (14). Using an intercross (B6 x 129)F2 generation, a genome-wide scan of the IR/IRS-1 double heterozygous mice revealed a locus with linkage to hyperinsulinemia on chromosome 14, a locus linked to hyperleptinemia on chromosome 7, and two loci with linkage to hyperglycemia on chromosomes 14 and 12, respectively (15). It is interesting that the donor strain in all four cases was B6. Recently, Colombo et al. (16) showed that mice with lipoatrophic diabetes have less insulin resistance in muscle and more insulin-resistant livers on a B6 background than on an FVB background, possibly as a result of differences in hepatic triglyceride clearance, although the exact loci contributing to this difference have not been mapped. In the present study, we explored in detail how high-fat dietinduced insulin resistance affects the B6 and 129 strains, as well as an F2 intercross between these strains, and asked whether the genes that control diet-induced insulin resistance are the same as those that control genetic insulin resistance and are dominant-acting genes and whether these genes could be identified using a combination of genetics and genomics.
F2 mice were obtained by breeding wild-type C57Bl/6 (B6) mice with 129 wild-type mice (strain 129S6/SvEvTac) from Taconic (Germantown, NY) to create an F1 outcross generation. The F1 generation was then intercrossed in a random manner to obtain a total of 100 F2 male mice. Fifty intercross F2 mice, 8 B6, and 8 129 male mice were maintained on a low-fat diet that contained 14% of calories from fat, 25% from protein, and 61% from carbohydrates (Taconic), and an equivalent number of mice were fed a high-fat diet containing 55% fat, 21% protein, and 24% carbohydrates (Harlan Teklad, Madison, WI). All mice were kept on the low-fat diet from the time of weaning until the end of puberty (at the age of 6 weeks), when the feeding study was initiated. The mice were maintained on the various diets for at least 18 weeks. Food intake was measured by placing the mice in individual cages and weighing the food every third day over 9 days. The mice were maintained on a 12-h light-dark cycle. All protocols for animal use were reviewed and approved by the Animal Care Committee of the Joslin Diabetes Center and were in accordance with the National Institutes of Health guidelines.
Phenotype analysis.
Indirect calorimetry and activity.
Genome-wide scan in 100 F2 intercross mice.
RNA isolation, cRNA preparation, array hybridization, and analysis.
Statistical analysis. Data analysis was performed using MapMaker/EXP 3.0 and MapMaker/QTL 1.1 (1820) as well as Map Manager version QTXb20 (21). MapMaker/QTL calculates the strength of the association between markers and traits as the log base 10 of the odds ratio (logarithm of odds [LOD] score). The permutation test (10,000 permutations) from Map Manager estimates an empirical genome-wide probability for a given likelihood ratio statistic (LRS). SPSS for Windows (version 11.5) was used for the other statistical analyses. P < 0.05 (two-tailed) was considered significant. Insulin levels in the F2 mice on both low- and high-fat diet deviated from a normal distribution and were log transformed for the quantitative trait locus (QTL) analysis.
Body weight gain. Parallel groups of 6-week-old B6 (n = 8), 129 (n = 8), and intercross (B6 x 129)F2 (n = 50) male mice were placed on a low-fat (14%) or high-fat (55%) diet (see RESEARCH DESIGN AND METHODS) for an 18-week period. Although starting body weights in all three groups were almost identical, during the 18 weeks of study, the F2 and B6 mice gained 18 g on the low-fat diet and 24 g on the high-fat diet, whereas the 129 mice gained only 6 g and 14.5 g on the two diets, respectively (Fig. 1). Indeed, the B6 and F2 mice gained more weight on a low-fat diet than the 129 mice gained even on the high-fat diet. The difference in weight gain between the low- and high-fat diets, however, was higher for the 129 strain (gained 35% more weight on the high-fat versus low-fat diet) than the B6 (gained 27% more weight on the high-fat versus low-fat diet). Thus, the B6 mouse gains more weight than the 129 independent of diet but is less sensitive to changes in dietary fat.
The variance is a measure of the spread of a distribution about its average value; thus, the larger the variance, the more heterogeneous is the group. At week 18 on a low-fat diet, the variance of the percentage weight gain was 307, 55, and 36% for the F2, B6, and 129 mice, respectively, whereas on the high-fat diet, the variance was 534, 118, and 132% for the F2, B6, and 129 mice, respectively. The larger variance in the F2 mice confirms the effects of the mixed genetic background. The larger variance on the high-fat versus low-fat diet also suggests a larger degree of heterogeneity among mice on factors such as individual dietary intake, activity, etc. The fact that the F2 mice show patterns of average weight gain almost identical to those of B6 suggests that B6 mice carry dominant genes in the control of body weight between these strains. In all cases, circulating leptin levels, which generally reflect adipose mass (22), paralleled the changes in body weight (Fig. 2A and B).
Effect of diet on glucose tolerance and insulin sensitivity. Both the F2 and the B6 mice showed an increase in insulin during the 18 weeks that was greater than that in the 129 mice, indicating higher levels of insulin resistance (Fig. 2C and D). This increase occurred faster when the mice were maintained on the high-fat diet. The variance of the insulin levels on low-fat diet at week 18 was 2,593, 1,251, and 0.4 ng/ml, and on high-fat diet it was 2,649, 1,242, and 1.3 ng/ml for the F2, B6, and 129 mice, respectively. Thus, the F2 mice show the effects of the mixed genetic background, and after 18 weeks, the spread of insulin was independent of the fat content of the diet. The random-fed plasma glucose levels on either diet were highest in the B6 mice ( 150 mg/dl), slightly lower in the F2 mice ( 140 mg/dl), and considerably lower in the 129 mice ( 105 mg/dl), and these did not change during the 18 weeks. Glucose tolerance tests performed after 18 weeks on the test diets revealed impaired glucose tolerance in B6 mice on both low- and high-fat diets. In contrast, the 129 mice showed impaired glucose tolerance only on the high-fat diet (Fig. 3A). The response to exogenous insulin during the insulin tolerance was also reduced in B6 mice on either diet, indicating greater insulin resistance in this strain (Fig. 3B).
Indirect calorimetry and activity. The increased weight gain of the B6 strain could potentially be due to a higher level of food intake or decreased energy expenditure. To decide which factors were most important, we subjected B6 and 129 mice to indirect calorimetry after they were maintained on the various diets for 6 weeks (age of 12 weeks) and caloric intake was assessed. Despite the higher level of weight gain, B6 mice ate less than the 129 on the low-fat diet and a similar number of calories on the high-fat diet (Fig. 3C), indicating the higher caloric efficiency (i.e., grams of weight gained per gram of food intake) of the B6 strain (Fig. 3D). Consistent with this, the metabolic rate, measured as O2 consumption, was 6% higher in the 129 mice than in the B6 mice during feeding on the low-fat diet, 5% higher during fasting on the high-fat diet (P = 0.01), and 14% higher (P < 0.001) during feeding on the high-fat diet (Fig. 4). In addition, metabolic rate showed significant increases in the 129 mice when changed from fasting to feeding (30% increase on low-fat diet and 10% increase on high-fat diet; Fig. 4A), whereas the B6 mice showed reduced thermic effect of food (Fig. 4). The higher metabolic rate of 129 mice was independent of activity. In fact, activity levels calculated by beam breaks over 24 h were 50% lower (P < 0.001) in the 129 than in the B6 mice under all conditions (Fig. 5A). The respiratory quotient, an indicator of the relative rates of substrate oxidation (23), was significantly lower in B6 than in 129 mice during feeding on both a low-fat (P = 0.027) and high-fat (P < 0.001) diet and during fasting on the low-fat diet (P < 0.001), indicative of a higher percentage of fat oxidation in the B6 mouse, whereas the respiratory quotient was similar in both strains during fasting on the high-fat diet (Fig. 5B). When both B6 and 129 mice were shifted from fasting to feeding, there was a change in substrate oxidation from fat to carbohydrate, as expected (Fig. 5B).
Identification of QTLs in intercross (B6 x 129)F2 mice. F2 mice (n = 100) were genotyped using a total of 114 polymorphic markers covering the 19 autosomal chromosomes and the X chromosome. Linkage maps for weight, weight gain, glucose, insulin, and leptin at the five time points on either the low- or high-fat diet were constructed by using MapMaker/EXP; QTLs were calculated using MapMaker/QTL and MapManager (18,21); and permutations test were performed using MapManager (24). The QTLs with the highest LOD scores are described below and in Table 1 and depicted in Fig. 6.
A linkage map of chromosome 14, based on genotyping data of 15 markers in F2 mice, revealed a QTL for insulin levels at week 18 in mice on a low-fat diet with a peak at marker D14Mit52 and a LOD score of 3.0, accounting for up to 25% of the variation of this trait (Fig. 6A, Table 1). The location of this QTL is virtually superimposable with the previously identified locus in the IR/IRS-1 double- heterozygous knockout mice (15). In the F2 generation, insulin levels of mice homozygous for the B6 allele at marker D14Mit52 were 8.9-fold higher than the insulin levels of mice homozygous for the 129 allele (P = 0.002), and mice heterozygous at the marker were intermediate between the other two genotypes, suggesting a gene- dosage effect (Fig. 6B). A QTL with linkage to weight gain at week 18 was also found on chromosome 14 around marker D14Mit192 with a LOD score of 3.0 (Fig. 6A, Table 1). Mice carrying two B6 alleles at marker D14Mit192 gained 5.2 g (or 22.5%) more weight at week 18 than mice carrying two 129 alleles (P = 0.007). A third QTL on chromosome 14 at marker D14Mit192 was found for glucose in mice that were maintained on a high-fat diet with a LOD score of 3.3 at week 12 (Table 1). For verifying the QTL analyses, tests with 10,000 permutations were carried out using MapManager (21,24). The permutation test estimates an empirical genome-wide probability for a given LRS (21). The permutation tests confirmed that the identified QTL with linkage to insulin at D14Mit52 as well as the QTLs at D14Mit192 were significant (LRS score at D14Mit52 for insulin = 11.3 with significance at 11.2, LRS score at D14Mit192 for weight gain = 13.0 with significance at 11.4, and LRS score at D14Mit192 for glucose = 14.8 with significance at 11.1). Our previous results have suggested a locus on chromosome 12 (D12Mit231) with linkage to hyperglycemia and another on chromosome 7 linked to leptin levels (15). In the present study, there was a QTL on chromosome 12 (D12Mit231) with linkage to leptin at week 2 in mice that were fed a low-fat diet (peak LOD score of 3.2; Fig. 6C). Surprisingly, mice homozygous for the 129 allele exhibited 1.9-fold higher leptin levels than B6 homozygous mice (P = 0.005; Table 1), whereas the previously identified locus on chromosome 7 linked with hyperleptinemia (15) was not statistically defined in these mice during diet-induced insulin resistance. Five other QTLs were identified in the genome-wide scan (Table 1). These included a hyperleptinemia locus on chromosome 10 (Fig. 6D), and on chromosome 3, we identified two separate loci linked with hyperinsulinemia and hyperleptinemia in mice that were fed a low-fat diet for 12 weeks (Table 1, Fig. 6E). Insulin and leptin levels, respectively, were significantly higher in mice homozygous for the B6 alleles at these markers (P = 0.002). On chromosome 11, a locus was significantly associated with insulin after 6 weeks on high-fat diet (Fig. 6F), and finally on chromosome 15 (Fig. 6G), we identified a locus associated with fasting plasma blood glucose.
Expression analysis of genes located in identified chromosomal areas. A similar enrichment of differentially expressed genes was observed in the region of marker D14Mit192, where there were seven significant changes at the P < 0.05 level of 36 measured defined genes. These included a 2.6-fold higher level of expression of the tyrosine kinase substrate Dok2 in fat of B6 over 129 mice, as well as bone morphogenetic protein 1 and esterase 10. Significant differences in fat of 129 over B6 included erythrocyte protein band 4.9 and protein phosphatase 3, and the latter gene was also upregulated in skeletal muscle, whereas the expression of SH3-domain protein 4 was higher in muscle of B6 than 129. The large number of changes on chromosome 14 between these two strains of mice suggests the action of some global cis-acting regulator of gene expression, such as altered patterns of DNA methylation, chromatic binding, or chromatic structure in these regions between these two strains of mice.
The development of diabetes and the associated metabolic syndrome in response to either genetic or acquired defects in both rodents and humans is dramatically affected by background genes. In mice, this was first demonstrated by Coleman et al. (5) for the ob/ob and db/db traits and has been subsequently observed in several studies, including analysis of the impact of heterozygosity for insulin receptor and IRS-1 deletion in mice (14,25), the phenotype of genetically induced lipodystrophy (16), and diet-induced obesity (26). In the insulin receptor/IRS-1 double-heterozygous knockout, we have previously shown that the background genes of B6 mice cause severe hyperinsulinemia and diabetes, whereas background genes of the 129 strain protect against diabetes (14). In the present study, we focused further on this problem by determining what accounts for the differences in response to diet-induced obesity in mice by defining the genetic background loci that contribute to the differences in phenotypes between the B6 and 129 mice that are subjected to environmentally induced insulin resistance and obesity and comparing those with the modifiers of genetically induced insulin resistance. We find that B6 mice are prone to becoming more obese on either high- or low-fat diets as they age compared with 129 mice. Although 129 mice are relatively lean, they are in fact highly susceptible to the content of fat in the diet, with a greater percentage increase in weight and greater deterioration in glucose tolerance than the B6 on a high- versus low-fat diet. Thus, mice of this strain are not obesity resistant, as previously described (4,13), but exhibit a complex interaction among genes, diet, age, and environment. In addition, we find that the higher body weight and tendency to glucose intolerance of the B6 mice are dominantly inherited when these two strains are crossed. The most obvious potential explanation for differences in body weight would be differences in food intake. Surprising, however, is that 129 mice ate significantly more calories than B6 on the low-fat diet and an equal number of calories on the high-fat diet but still remained leaner than their B6 counterparts. Thus, the B6 mice showed higher feeding efficiency, i.e., greater weight gain per calorie consumed, as compared with 129 mice. This would represent an example of the "thrifty gene" phenotype and represent an advantage during times of food shortage (27). In the absence of excessive food intake, low energy expenditure would be the major potential risk factor for body weight gain (28). Although B6 mice are actually more active than 129 mice, B6 mice have a lower metabolic rate than the 129 during feeding on both low- and high-fat diets and also during fasting on the high-fat diet. B6 mice also exhibit a smaller change in metabolism in response to a shift from a fasted to a fed state than 129 mice and a smaller change in response to the shift from a low- to a high-fat diet, i.e., reduced diet-induced thermogenesis. Whether the higher basal metabolic rate and higher diet-induced thermogenic response in 129 mice accounts entirely for the diminished weight is hard to predict. It is interesting that studies have shown that diet-induced thermogenesis is higher in lean humans compared with obese individuals (29). Furthermore, the higher metabolic rate in the 129 mouse in the basal state is independent of activity level, because these mice show approximately a 50% lower level of activity than the B6. Besides a low metabolic rate, a low-fat versus carbohydrate oxidation ratio, i.e., a high respiratory quotient, has been suggested to be a risk factor for body weight gain (23,30,31). Contrary to this dogma, B6 mice exhibit a slightly lower respiratory quotient than the 129 mice. An increase in body fat causes an increase in circulating leptin levels, which will normally decrease the energy intake and increase energy expenditure (32), yet in the B6 mice, the ratio of leptin to body weight and fat mass was much higher than that in the 129 mice. Thus, the B6 mice develop leptin resistance with age and obesity, mimicking humans with diet-induced obesity (15,33). It is interesting that leptin mRNA on high-fat diet is similar in B6 and 129 mice despite the differences in circulating leptin (manuscript in preparation), suggesting differences in either leptin translation efficiency or protein half-life between strains. In parallel with the higher body weight, the intraperitoneal glucose tolerance testing showed greater hyperglycemia and little diet-induced change in the B6 mice, whereas 129 mice were glucose tolerant on the low-fat diet but intolerant on the high-fat diet. Likewise, the B6 mice exhibited insulin resistance during the insulin tolerance test on both the low- and high-fat diets. In contrast, 129 mice were relatively insulin sensitive even on a high-fat diet. The development of glucose intolerance in 129 mice on a high-fat diet is in contrast to our previous studies that have shown that 129 mice do not develop glucose intolerance when subjected to genetically induced insulin resistance as a result of heterozygosity of IR and IRS-1 (14). The B6 mice, conversely, become insulin resistant and hyperinsulinemic on either diet and also are more prone to become diabetic when they carry the double heterozygous deletion for IR and IRS-1 (14). Characterization of the (B6 x 129)F2 intercross suggests that the differences between B6 and 129 mice are primarily determined by dominant-acting genes of the B6 mouse (24). In particular, we found that a locus on chromosome 14 between markers D14Mit55 and D14Mit52 previously identified as a modifier of genetically induced insulin resistance (15) also plays an important role in response to dietary insulin resistance.
Of the potential candidate genes in this area that are represented in the Affymetrix microarray, a disproportionate number of genes exhibited differential expression between B6 and 129 mice. These included several genes that are known to affect metabolism, such as PKC Wnt-5a is a member of a family of paracrine and autocrine factors that regulate cell growth and cell fate (38). Wnt signaling prevents preadipocytes from differentiating into adipocytes (38), and this could fit with the findings of the present study, showing that Wnt expression is upregulated in skeletal muscle, liver, and fat of the lean 129 strain as compared with the obese B6, in some cases by almost fivefold.
The expression level of transketolase was also significantly different between B6 and 129 mice; however, in this case, the expression was significantly higher in skeletal muscle of 129 versus B6, whereas it was significantly lower in liver of 129 versus B6, and the expression was similar in white fat. Transketolase is a ubiquitous enzyme involved in multiple metabolic pathways. Haploinsufficiency of transketolase results in reduced adipose tissue mass (39), whereas the deletion of PKC A separate QTL at marker D14Mit192 was identified with linkage to both weight gain and hyperglycemia and inherited by the B6 strain. On the basis of differential expression, Dok2 (also known as p56[Dok2], FRIP, or Dok-R) is an interesting candidate gene located in the area around marker D14Mit192, which is highly upregulated in epididymal fat of B6 mice in comparison with 129. Dok2 is a RasGAP-binding protein suggested to have an important role downstream of growth factor receptors as a potential negative regulator of signal transduction (42). Additional studies will be needed to determine whether this gene or other genes revealed in Table 2 play a role in the development of hyperinsulinemia or insulin resistance. Furthermore, the list of genes under the peaks given in the Table 2 may potentially be the quantitative trait genes; however, several unidentified candidate genes may still be discovered in the areas. A QTL with linkage to leptin at week 2 was detected around marker D12Mit231 and is inherited from the 129 genome. The locus overlaps the hyperglycemia locus that was identified in the insulin receptor/IRS-1 double-heterozygous mice that interacts significantly with a locus on chromosome 7 in the control of leptin (15). In the New Zealand Obese mouse, this same marker (D12Mit231) has been shown to affect adiposity, BMI, and leptin (43) and is a diabesity locus when interacting with several other loci (44). Loci linked with hyperleptinemia were also discovered on chromosomes 3 and 10 in the present model. Because of the low number of markers in the area of the QTL on chromosome 10 as well as chromosome 15, there is some uncertainty about the significance of those QTLs. Reed et al. (45) performed an analysis of an F2 intercross of C57BL/6ByJ and 129P3/J and identified QTLs for body weight, length, and adiposity on chromosome 2. We did not identify any significant QTLs on chromosome 2. This discrepancy may be due to small differences of the inbred mouse substrains being used in the two studies (Jackson versus Taconic, 129Ps/J versus 129S6/SvEvTac, and C57Bl/6 versus C57BL/6ByJ), differences in dietary fat intake (4.4 vs. 9% fat) that were used to provoke the insulin-resistance phenotype, and differences in the age at the time of study or time course of study or other unidentified variables.
Combining the mapping data with the microarray expression profile provides an opportunity to establish a potential link between genotype and phenotype and thereby improve the chances of finding new diabetogenes and pathways. If the link between QTL and gene expression can be proved directly, then this would suggest a cis-acting gene, i.e., one in which a DNA variation in the gene itself (most likely the promoter) controls the expression level. As more data on the complete sequence of these candidate genes in both strains of mice become available, this hypothesis can be tested. However, a LOD score can also be due to a trans-acting effect, i.e., a gene at the defined locus acts on other genes to affect their transcription. It has previously been suggested that it is easier to detect effects that are caused by DNA variation acting on the gene itself (46), and this is also suggested by our data as a result of the high number of genes with differential expression at the QTLs (Table 2). For example, around marker D14Mit192, In conclusion, B6 mice are prone to becoming obese, hyperleptinemic, insulin resistant, and glucose intolerant independent of the fat content in food, whereas 129 mice exhibit more modest weight gain and glucose intolerance only on a high-fat diet. This is due to higher caloric efficiency and reduced diet-induced thermogenesis in the B6 mouse, two manifestations of the "thrifty gene." This "thrifty" phenotype with weight gain, hyperinsulinemia, and hyperleptinemia is dominantly inherited when these strains are bred. QTL and gene-array expression data suggest that some of the major genes that control hyperinsulinemia/insulin resistance both as a result of genetic modification and obesity are located on murine chromosome 14 in the region of marker D14Mit52.
This work was supported by grants DK31036 and DK45935 to C.R.K. and by core support from the Diabetes Genome Anatomy Project (DK60837-02). We thank Scott Lannon, Laureen Mazzola, Katherine Gartrell, and Pei Lin for technical help and Drs. Andy Norris, Atul Butte, and Isaac Kohane for useful discussions.
K.A. is currently affiliated with Novo Nordisk, Bagsvaerd, Denmark. Address correspondence and reprint requests to C. Ronald Kahn, MD, Research Division, Joslin Diabetes Center, One Joslin Place, Boston, MA 02215. E-mail: c.ronald.kahn{at}joslin.harvard.edu Received for publication March 3, 2004 and accepted in revised form August 23, 2004
Abbreviations: IRS-1, insulin receptor substrate-1; LOD, logarithm of odds; LRS, likelihood ratio statistic; PKC, protein kinase C; QTL, quantitative trait locus
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||