Table of Contents Author Guidelines Submit a Manuscript
Evidence-Based Complementary and Alternative Medicine
Volume 2011, Article ID 215653, 7 pages
Research Article

Resequencing DCDC5 in the Flanking Region of an LD-SNP Derived from a Kidney-Yang Deficiency Syndrome Family

Department of Fundamental Medicine, Chengdu University of Traditional Chinese Medicine, Chengdu 610075, China

Received 15 July 2010; Revised 17 February 2011; Accepted 24 February 2011

Copyright © 2011 Li Ping 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.


Objective. To explore the genetic traits of Kidney-yang deficiency syndrome (KDS). Design. Twelve KDS subjects and three spouses from a typical KDS family were recruited. Their genomic DNA samples were genotyped by Affymetrix 100K single-nucleotide polymorphism (SNP) arrays. The linkage disequilibrium (LD) SNPs were generated using GeneChip DNA analysis software (GDAS, Affymetrix). Genes located within 100 bp of the flanks of LD SNPs were mined via GeneView. 29 exons of the doublecortin domain containing 5 (DCDC5), a representative gene within the flank of an LD SNP, were resequenced. Results. Five LD SNPs display midrange linkage with KDS. Two genes with established functions, DCDC5 and Leucyl-tRNA synthetase, were mined in the flanks of LD SNPs. Resequencing of DCDC5 revealed a nonsynonymous variation, in which 3764T/A was replaced by C/G. Accordingly, the Ser1172 was substituted by Pro1172. The S1172P substitution effect was evaluated as “possibly damaging” by PolyPhen. Conclusion. We have identified a genomic variation of DCDC5 based on the LD SNPs derived from a KDS family. DCDC5 and other genes surrounding these SNPs display some relationships with key symptoms of KDS.

1. Introduction

The genetic features of kidney-yang deficiency syndrome (KDS) still remain obscure. However, there have been few investigations that have focused on genetic alterations that might correlate with this disease. Wang et al. [1] suggested that the cause and progression of KDS involves not only environmental factors but genetic elements as well. Pan et al. [2] analyzed the transcriptomic patterns of KDS subjects utilizing microarrays. The KDS transcriptome module was further analyzed by Yang et al. [3] and Wu et al. [4]. However, these studies did not identify any genetic alterations associated with this disease. Thus it was the purpose of this study to identify possible genetic alterations that could contribute to the development of KDS.

Re-sequencing is a promising method that is particularly valuable for the genetic studies of complex diseases. For example, Dapprich et al. [5] described a simple and automated method of targeting specific sequences of genomic DNA linked to SNPs of interest and associated with certain complex traits or drug responses. Smyth et al. [6] tested the shared and distinct non-HLA loci between two inflammatory disorders, type 1 diabetes, and celiac disease, indicating that a genetic susceptibility to both diseases shares common alleles. Wermter et al. [7] employed resequencing and relevant methods to analyze the mutations of melanin-concentrating hormone receptor 1 gene and reported that several alleles and/or haplotypes might be involved in juvenile-onset obesity. Indeed, nucleotide sequencing, originally described by Sanger and colleagues, is still the gold standard for the identification of genetic mutations and polymorphisms.

KDS and other TCM syndromes share many features with complex diseases. TCM, by its sheer nature, is based on the analysis of the simultaneous occurrence of pathological processes on a macroscopic level and thus provides a holistic approach complementary to the trend toward disaggregation of complex phenotypes as is practiced in conventional medicine [8]. For example, a typical complex disease, diabetes, has been traditionally named by TCM as thirst-and-emaciated syndrome, exactly manifesting the major symptoms of diabetes. Hence, we propose a hypothesis that TCM syndromes, just like hypertension, diabetes, and other complex diseases defined by Western medicine, possess intricate genetic traits involving DNA structure variations, such as SNP alleles, SNP linkage disequilibrium (LD), and gene mutations.

The purpose of this paper is to explore the DNA structure variations of KDS by a resequencing approach. Based on the importance of the resequencing method upon the genetic traits of complex diseases, and the shared features between complex diseases and TCM syndromes, we adopted this method to investigate genomic variations of KDS. We firstly employed a classic method universally applied in human genetics, pedigree approach, to report that KDS patients undergo overall attenuated functions in the mass-energy-information-carrying network [9]. Second, we utilized a genome-wide method, SNP array, to retrieve the LD SNPs within the KDS family [8]. In the present paper, we employed resequencing to identify the defined DNA mutations of genes located within 100 bp of the flanks of LD SNPs, and we report the results of LD SNPs and sequencing polymorphisms of doublecortin domain containing 5 (DCDC5), an important gene that may play an essential role in the pathology of KDS.

2. Materials and Methods

2.1. Collection of Genomic DNA Samples from a KDS Family

This work was approved by the Ethics Board, Chengdu University of TCM, and all participants signed the informed consent form before examination [8]. The process of collecting KDS pedigree was described previously [8, 9]. Briefly, we used a 40-item scoring table of KDS as a criteria to recruit KDS subjects in Chengdu, southwest China, and defined participants as KDS subjects (total score ≥12 points) or healthy individuals (total score ≤5 points). In order to improve the validity and reproducibility of the TCM diagnosis, every participant was diagnosed by five independent TCM physicians, and the examinations were repeated for three consecutive years (from 2003 through 2005), all on the first Saturday of December [8, 9].

Three mL of blood sample was collected from each subject. The genomic DNA was isolated by conventional phenol/chloroform extraction, washed with 95% ethanol, re-suspended in TE buffer, and stored at −80°C until further use.

2.2. Genotyping, LD Estimating, and Mining the Genes around LD SNPs

The detailed description of these procedures has been reported previously [8, 9]. Briefly, DNA samples were genotyped with the Affymetrix GeneChip Mapping 100K kit, at a commercial laboratory (National Engineering Center for Biochip at Shanghai, China). Genotype data of all the subjects were generated using a GeneChip DNA analysis Software (GDAS, Affymetrix). The pedigree information, allele frequencies, and map position of the SNPs were combined with the genotype data. To determine the LD SNPs from this KDS family, we performed a simulation study using Agilent Varia software [10, 11]. The following steps were automatically executed in the computer workstation of the laboratory. (a) likelihoods, ( ), are calculated for each possible value ranging from 0 to 0.5, in increments of 0.02, that is, 0, 0.02, 0.04,…, 0.46, 0.48, 0.5. (b) Choose which maximizes the likelihood ratio, ( )/ (1/2), of the occurrence of linkage to no linkage. (c) Calculate the lod adds (LOD) score for this .

Genes located within 100 bp of the flanks of LD SNPs were mined via an authoritative web tool, GeneView ( Their biological functions were then annotated by employing FatiGO (http://fatigo.bioinfo.ochoa.fib/es).

2.3. Re-Sequencing of Doublecortin Domain Containing 5 (DCDC5)

All of the 29 exons of DCDC5 were covered by 29 pairs of PCR primers (Table 1). These primers were designed based on the DCDC5 version issued on 26 June, 2007 (NM_198462.2). A website primer design platform, Primer3 (, was employed under the default primer selection conditions. All samples, including twelve KDS patients and three spouses from the KDS pedigree were resequenced and compared for variation discovery in DCDC5. A detailed laboratory protocol can be found on the SeattleSNPs website ( PCR products were sequenced using each pair of above PCR primers and BigDye Terminator Cycle Sequencing v2.0 kit (Applied Biosystems, Applera, USA). Sequencing reactions were electrophoresed on an ABI 377 automated sequencers. Finally, the substitution effect was evaluated by a free software, PolyPhen (

Table 1: Primers designed for the resequencing of DCDC5. A website primer design platform, Primer3, was employed under the default primer selection conditions to design these primers based on the genbank accession number NM_198462.2. All samples, including twelve KDS patients and three spouses from the KDS family were resequenced and electrophoresed on an ABI 377 automated sequencers.

3. Results

3.1. General Conditions of a Representative KDS Family

A typical KDS pedigree, including 17 KDS patients distributed in four generations who live in the area surrounding Chengdu, China, was selected as the subject for the present study. Their KDS symptoms were confirmed by repeated TCM diagnosis and clinic biochemistry assays [8]. Twelve available KDS subjects (Figure 1) were collected for genetic analysis. Simultaneously, three healthy spouses, whose life environment and ages were best matched with those of the KDS subjects, were recruited as non-KDS controls.

Figure 1: Pedigree tree of the collected KDS family. This family comprise four generations and 35 members that display an autosomal dominant model of inheritance. The KDS subjects were evaluated by a 40-item scoring table based on TCM criteria system of KDS promulgated by the Health Department of China (GB/T15657-1995). Dark color indicates KDS. The mark “?” indicates those were not evaluated due to death, nonaccessibility, or babies who were too young to be easily checked. The mark “*” indicates those genome SNP were checked by SNP chip. KDS: kidney-yang deficiency syndrome. TCM: traditional chinese medicine. SNP: single nucleotide polymorphism.
3.2. Two Genes with Established Functions Were Mined in the Flanks of LD SNPs

Five SNPs obtained from the KDS family were identified as LD SNPs. The sequences of alleles of LD SNPs identified from the KDS family were displayed in Table 2. It has been defined that indicates tight linkage; while , midrange linkage and , loose linkage [10]. Results displayed that before their values achieved 0.2, all LOD values of these SNPs were greater than 1.0. Hence, these LD SNPs demonstrated midrange linkage with KDS.

Table 2: Five SNPs of linkage disequilibrium derived from a KDS family. A conventional approach, lod adds Score method (LODs), was employed to calculate the levels of recombination condition. The LOD values (or Z value) were reckoned at different values ( ). The criteria standard of linkage are, : support linkage, : linkage. Accordingly, indicates tight linkage, , midrange linkage, , loose linkage.

We then examined the genes within 100 base pairs of flanks of these LD SNPs (Figure 2), and observed two genes with established functions. One gene, doublecortin domain containing 5 (DCDC5, NM_198462), is referred to rs514207. DCDC5 is located in chromosome 11p14.1-p13. The established functions of DCDC5 are sugar binding, interacting selectively with mono-, di-, or trisaccharide carbohydrate, and involving certain intracellular signaling cascades [1217]. Thus, DCDC5 might be one of the key factors that monitor saccharides and energy metabolism. Another gene, encoding Leucyl-tRNA synthetase (LARS), is located near rs1054020. It is also a pleiotropic gene that involved in cysteine-tRNA ligase activity, tRNA aminoacylation, ATP binding, protein biosynthesis, and other biological processes [18, 19].

Figure 2: Chromatograms of the LD SNPs. The upper part illustrates chromosome 11p14.1-p13 and flanks with respect to the gene DCDC5. The lower part is the result of linkage analysis corresponding to the assay of SNP arrays based on a KDS family.
3.3. SER/PRO Alteration Was Observed in Exon 27 of DCDC5

Results of resequencing of DCDC5 show several variations and/or mutations (Table 3). It is worth noting that in eight of the twelve KDS participants, the 3764T/A was replaced by C/G. As a nonsynonymous variation, the wild-type Ser1172 of DCDC5 was substituted by Pro1172 in the vicinity of C-terminal (Figure 3). This novel variation occurred within exon 27 of DCDC5. This mutation has not been previously reported and does not show any relationship with the established variants of DCDC5 [1217]. Results of the tolerable or deleterious probability for S1172P substitution were +1.040 for score 1, and −0.462 for score 2, which could be marked as “possibly damaging” by PolyPhen.

Table 3: Resequencing results of DCDC5 from a KDS pedigree.
Figure 3: Schematic representation of the Ser/Pro substitution occurred in DCDC5. Totally that there are 1258 residues of amino acids. The first and second DCX domains are located at 4–94 and 491–570, respectively. The Rincin-B-Lectin (RBL) domain is located at 238–376. DCDC5: doublecortin domain containing 5. DCX: doublecortin.

4. Discussion

A resequencing approach, based on the results of LD SNPs, is a promising strategy for exploring the genetic traits of complex diseases. Most functionally important SNPs are located within genes and other critical gene regulatory regions, such as the promoters, 5′ and 3′ untranslated regions, and predicted splice sites [2023]. Therefore, SNPs can alter gene  expression and phenotypes through a variety of mechanisms. Thus far, only a few SNPs have been identified as LD SNPs representing strong candidate biomarkers for certain diseases [8, 10, 22, 23]. However, given the extreme complexity of genomic structure, resequencing and relevant methods that can identify the specific variants in a sequence of interest compared with a known genomic sequence is a monumental task [22]. Hence, combined analysis of LD SNPs and resequencing can focus on the targets and identify the defined variations, especially in the exploration of genetic characteristics of complex diseases. Results of this present work also reveal that sequence analysis associated with LD SNP analysis is likely to be a useful tool to substantially increase our understanding of the genetic basis of KDS.

Results of resequencing reveal a novel variation occurring in the vicinity of the C-terminal of DCDC5. Nagase et al. [13] cloned DCDC5 from a fetal brain cDNA library. This gene is located in chromosome 11p14.1-p13 and is weakly expressed in male reproductive system, testies and eyes. Thus far, four transcripts, several amino acid variants and about 100 SNPs of DCDC5 have been documented [12]. In the present work, a novel genomic variant of DCDC5 was identified in family members suffering from KDS. Accordingly, within the NCBI reference sequence NP_940864.2, the Ser1172 was replaced by Pro1172. This variant has not been found in any established domains of DCDC5, nor is it located in any active enzyme site (e.g., enzyme site for serine/threonine protein kinase). This substitution could possibly damage the structure and/or function of DCDC5. Thus, functional consequences of this causative variant in DCDC5 associated with susceptibility to KDS remains to be confirmed. Notwithstanding a functional change in the protein, our findings do document a new variant associated with KDS, characterized by a novel SNP within the genomic region of DCDC5. To our knowledge, there have been no previous studies reporting genomic variations associated with TCM syndrome.

As a relatively new studied gene, DCDC5 is of importance in saccharide binding, energy metabolism and diverse signal transduction pathways. Recent reports [12, 13] revealed that DCDC5 located in chromosome 11, one of the most gene- and disease-enrich chromosomes in the human genome. Human DCDC5 contains two DCX domains separated by a ricin-like domain (RICIN_B_LECTIN) (Figure 2). DCX domains serve as protein-interaction platforms and are likely to mediate signal transduction pathways. Mutations in DCX genes result in abnormal neuronal migration, epilepsy, mental retardation, inherited blindness, and dyslectic reading disabilities [12]. However, the quest for how DCX domain-containing proteins exert their multiple functions is still unresolved [15]. In contrast, RICIN_B_LECTIN domains interact selectively with mono-, di-, or trisaccharide carbohydrates [16] and involve signaling cascades such as the regulation of cytoskeleton [17, 18]. Therefore, DCDC5 is an important and pleiotropic gene involving diverse functions such as saccharides binding, energy metabolism, cell cycles, and intracellular signaling cascades.

DCDC5 is likely to be involved in the physiological functions and pathological process of KDS. From the viewpoint of TCM, the major biological functions of DCDC5 could be included into the TCM concept of the kidney. TCM kidney, as a critical monitor of metabolic, immune, and reproductive systems, works closely with the neural-endocrine-immune network [8, 9, 24]. Consequently, many vital functions and activities of KDS patients might be globally attenuated or disrupted. The major symptoms of KDS, such as dislike of cold, persistent cold extremities, and favor of warmth are mostly due to the insufficient energy metabolism, thus resulting in an abnormal low level of energy output [2427]. Apparent overlap in major physiological functions and pathology between DCDC5 and KDS suggests certain relationships that still need to be determined. In other words, DCDC5 is an important but less studied gene that is worthy of consideration for further analysis by scientists investigating the pathology of KDS.

The main limitation of this work is that we only examined exons of DCDC5. Variations located within other critical sites of DCDC5, such as the introns, promoters, 5′ and 3′ untranslated regions, and predicted splice sites, have not been resequenced. Further, we only checked twelve KDS subjects from a typical KDS family. Therefore, analysis of other KDS families and individuals with KDS are necessary.

In conclusion, we report in this paper a genomic variation of DCDC5, which could possibly damage this protein, based on the analysis of LD SNPs from a KDS family. DCDC5 and other genes located in the flanks of these SNPs display relationships with key symptoms of KDS. However, further work is required to study their function in order to establish the detailed relationships of DCDC5 to KDS. The merits of linkage analysis combined by resequencing may warrant further exploration of the genetic features of TCM syndromes. Regardless, our current data provides new and relevant information on gene associations that may correlate with the development and progression of KDS.


The authors thank the National Scientific Fund Committee (NSFC) of China (no. 30371709 and no. 30582288) for financial support. Special gratitude is extended to Dr. Howard Young, National Cancer Institute-Frederick, Frederick, MD, USA for the language editing. The authors also wish to express their appreciation to all the patients who generously volunteered to participate in this study.


  1. J. G. Wang, L. Pan, B. Wu, and M. Wang, “Familial characteristics of kidney-yang deficiency and cold syndrome,” Journal of Toxicology and Environmental Health—Part A, vol. 69, no. 21, pp. 1939–1950, 2006. View at Publisher · View at Google Scholar · View at Scopus
  2. L. Pan, M. Wang, J. G. Wang, B. Wu, and K. M. Hui, “Clinical and molecular evaluation of warming and tonic herb treatment for sibling patients of a typical kidney-yang deficiency family,” American Journal of Chinese Medicine, vol. 34, no. 3, pp. 387–400, 2006. View at Publisher · View at Google Scholar · View at Scopus
  3. L. Yang, M. Wang, W. Wu, and L. Zhang, “Transcriptome analysis of cold syndrome using microarray,” American Journal of Chinese Medicine, vol. 35, no. 4, pp. 609–620, 2007. View at Publisher · View at Google Scholar · View at Scopus
  4. B. Wu, M. Wang, J. G. Wang, L. Pan, and K. M. Hui, “Establishment of a discriminant mathematical model for diagnosis of deficiency-cold syndrome using gene expression profiling,” Journal of Alternative and Complementary Medicine, vol. 12, no. 8, pp. 751–761, 2006. View at Publisher · View at Google Scholar · View at Scopus
  5. J. Dapprich, D. Ferriola, E. E. Magira, M. Kunkel, and D. Monos, “SNP-specific extraction of haplotype-resolved targeted genomic regions,” Nucleic Acids Research, vol. 36, no. 15, article e94, 2008. View at Publisher · View at Google Scholar · View at Scopus
  6. D. J. Smyth, V. Plagnol, N. M. Walker et al., “Shared and distinct genetic variants in type 1 diabetes and celiac disease,” The New England Journal of Medicine, vol. 359, no. 26, pp. 2767–2777, 2008. View at Publisher · View at Google Scholar · View at Scopus
  7. A. K. Wermter, K. Reichwald, T. Büch et al., “Mutation analysis of the MCHR1 gene in human obesity,” European Journal of Endocrinology, vol. 152, no. 6, pp. 851–862, 2005. View at Publisher · View at Google Scholar · View at Scopus
  8. W. J. Ding, Y. Z. Zeng, W. H. Li et al., “Identification of linkage disequilibrium SNPs from a Kidney-yang deficiency syndrome pedigree,” American Journal of Chinese Medicine, vol. 37, no. 3, pp. 427–438, 2009. View at Publisher · View at Google Scholar · View at Scopus
  9. W. J. Ding, S. L. Yan, Y. Z. Zeng et al., “Insufficient activity of MAPK pathway is a key monitor of Kidney-Yang Deficiency Syndrome,” Journal of Alternative and Complementary Medicine, vol. 15, no. 6, pp. 653–660, 2009. View at Google Scholar · View at Scopus
  10. B. D. Horne, J. F. Carlquist, L. A. Cannon-Albright et al., “High-resolution characterization of linkage disequilibrium structure and selection of tagging single nucleotide polymorphisms: application to the cholesteryl ester transfer protein gene,” Annals of Human Genetics, vol. 70, no. 4, pp. 524–534, 2006. View at Publisher · View at Google Scholar · View at Scopus
  11. A. Middeldorp, S. Jagmohan-Changur, Q. Helmer et al., “A procedure for the detection of linkage with high density SNP arrays in a large pedigree with colorectal cancer,” BMC Cancer, vol. 7, article 6, 2007. View at Publisher · View at Google Scholar · View at Scopus
  12. O. Reiner, F. M. Coquelle, B. Peter et al., “The evolving doublecortin (DCX) superfamily,” BMC Genomics, vol. 7, article 188, 2006. View at Publisher · View at Google Scholar · View at Scopus
  13. T. Nagase, R. Kikuno, K. I. Ishikawa, M. Hirosawa, and O. Ohara, “Prediction of the coding sequences of unidentified human genes. XVII. The complete sequences of 100 new cDNA clones from brain which code for large proteins in vitro,” DNA Research, vol. 7, no. 2, pp. 143–150, 2000. View at Google Scholar · View at Scopus
  14. C. J. Hoggart, J. C. Whittaker, M. De Iorio, and D. J. Balding, “Simultaneous analysis of all SNPs in genome-wide and re-sequencing association studies,” PLoS Genetics, vol. 4, no. 7, Article ID e1000130, 2008. View at Publisher · View at Google Scholar · View at Scopus
  15. J. M. Weimer and E. S. Anton, “Doubling up on microtubule stabilizers: synergistic functions of doublecortin-like kinase and doublecortin in the developing cerebral cortex,” Neuron, vol. 49, no. 1, pp. 3–4, 2006. View at Publisher · View at Google Scholar · View at Scopus
  16. Y. Liu, A. J. Chirino, Z. Misulovin et al., “Crystal structure of the cysteine-rich domain of mannose receptor complexed with a sulfated carbohydrate ligand,” Journal of Experimental Medicine, vol. 191, no. 7, pp. 1105–1115, 2000. View at Publisher · View at Google Scholar · View at Scopus
  17. Q. Liu, J. Zhou, S. P. Daiger et al., “Identification and subcellular localization of the RP1 protein in human and mouse photoreceptors,” Investigative Ophthalmology and Visual Science, vol. 43, no. 1, pp. 22–32, 2002. View at Google Scholar · View at Scopus
  18. E. Seiradake, W. Mao, V. Hernandez et al., “Crystal structures of the human and fungal cytosolic Leucyl-tRNA synthetase editing domains: a structural basis for the rational design of antifungal benzoxaboroles,” Journal of Molecular Biology, vol. 390, no. 2, pp. 196–207, 2009. View at Publisher · View at Google Scholar · View at Scopus
  19. S. W. Lue and S. O. Kelley, “A single residue in leucyl-tRNA synthetase affecting amino acid specificity and tRNA aminoacylation,” Biochemistry, vol. 46, no. 15, pp. 4466–4472, 2007. View at Publisher · View at Google Scholar · View at Scopus
  20. M. Ingman and U. Gyllensten, “SNP frequency estimation using massively parallel sequencing of pooled DNA,” European Journal of Human Genetics, vol. 17, no. 3, pp. 383–386, 2009. View at Publisher · View at Google Scholar · View at Scopus
  21. E. J. Peters, S. L. Slager, G. D. Jenkins et al., “Resequencing of serotonin-related genes and association of tagging SNPs to citalopram response,” Pharmacogenetics and Genomics, vol. 19, no. 1, pp. 1–10, 2009. View at Publisher · View at Google Scholar · View at Scopus
  22. P. S. Tarpey, R. Smith, E. Pleasance et al., “A systematic, large-scale resequencing screen of X-chromosome coding exons in mental retardation,” Nature Genetics, vol. 41, no. 5, pp. 535–543, 2009. View at Publisher · View at Google Scholar · View at Scopus
  23. C. Tocharoentanaphol, S. Promso, D. Zelenika et al., “Evaluation of resequencing on number of tag SNPs of 13 atherosclerosis- related genes in Thai population,” Journal of Human Genetics, vol. 53, no. 1, pp. 74–86, 2008. View at Publisher · View at Google Scholar · View at Scopus
  24. M. Q. Wang, B. Wu, R. Feng et al., “An investigation report on the cold syndrome family,” Journal of Integrated TCM and Western Medicine, vol. 12, no. 22, pp. 2385–2386, 2003. View at Google Scholar
  25. H. J. Kim, S. Y. Hwang, J. H. Kim et al., “Association between genetic polymorphism of multidrug resistance 1 gene and sasang constitutions,” Evidence-Based Complementary and Alternative Medicine, vol. 6, no. 1, pp. 73–80, 2009. View at Publisher · View at Google Scholar · View at Scopus
  26. I. Lampronti, M. T. H. Khan, M. Borgatti, N. Bianchi, and R. Gambari, “Inhibitory effects of Bangladeshi medicinal plant extracts on interactions between transcription factors and target DNA sequences,” Evidence-Based Complementary and Alternative Medicine, vol. 5, no. 3, pp. 303–312, 2008. View at Publisher · View at Google Scholar · View at Scopus
  27. K. Nagashima, T. Yoda, T. Yagishita, A. Taniguchi, T. Hosono, and K. Kanosue, “Thermal regulation and comfort during a mild-cold exposure in young Japanese women complaining of unusual coldness,” Journal of Applied Physiology, vol. 92, no. 3, pp. 1029–1035, 2002. View at Google Scholar · View at Scopus