A quantitative metric to identify critical elements within seafood supply networks

A theoretical basis is required for comparing key features and critical elements in wild fisheries and aquaculture supply chains under a changing climate. Here we develop a new quantitative metric that is analogous to indices used to analyse food-webs and identify key species. The Supply Chain Index (SCI) identifies critical elements as those elements with large throughput rates, as well as greater connectivity. The sum of the scores for a supply chain provides a single metric that roughly captures both the resilience and connectedness of a supply chain. Standardised scores can facilitate cross-comparisons both under current conditions as well as under a changing climate. Identification of key elements along the supply chain may assist in informing adaptation strategies to reduce anticipated future risks posed by climate change. The SCI also provides information on the relative stability of different supply chains based on whether there is a fairly even spread in the individual scores of the top few key elements, compared with a more critical dependence on a few key individual supply chain elements. We use as a case study the Australian southern rock lobster Jasus edwardsii fishery, which is challenged by a number of climate change drivers such as impacts on recruitment and growth due to changes in large-scale and local oceanographic features. The SCI identifies airports, processors and Chinese consumers as the key elements in the lobster supply chain that merit attention to enhance stability and potentially enable growth. We also apply the index to an additional four real-world Australian commercial fishery and two aquaculture industry supply chains to highlight the utility of a systematic method for describing supply chains. Overall, our simple methodological approach to empirically-based supply chain research provides an objective method for comparing the resilience of supply chains and highlighting components that may be critical.

Secondary structure element alignment kernel method for prediction of protein structural classes

In this paper, we aim at predicting protein structural classes for low-homology data sets based on predicted secondary structures. We propose a new and simple kernel method, named as SSEAKSVM, to predict protein structural classes. The secondary structures of all protein sequences are obtained by using the tool PSIPRED and then a linear kernel on the basis of secondary structure element alignment scores is constructed for training a support vector machine classifier without parameter adjusting. Our method SSEAKSVM was evaluated on two low-homology datasets 25PDB and 1189 with sequence homology being 25% and 40%, respectively. The jackknife test is used to test and compare our method with other existing methods. The overall accuracies on these two data sets are 86.3% and 84.5%, respectively, which are higher than those obtained by other existing methods. Especially, our method achieves higher accuracies (88.1% and 88.5%) for differentiating the α + β class and the α/β class compared to other methods. This suggests that our method is valuable to predict protein structural classes particularly for low-homology protein sequences. The source code of the method in this paper can be downloaded at http://math.xtu.edu.cn/myphp/math/research/source/SSEAK_source_code.rar.

Development and validation of the Emergency Department Assessment of Chest pain Score and 2h accelerated diagnostic protocol

Objective Risk scores and accelerated diagnostic protocols can identify chest pain patients with low risk of major adverse cardiac event who could be discharged early from the ED, saving time and costs. We aimed to derive and validate a chest pain score and accelerated diagnostic protocol (ADP) that could safely increase the proportion of patients suitable for early discharge. Methods Logistic regression identified statistical predictors for major adverse cardiac events in a derivation cohort. Statistical coefficients were converted to whole numbers to create a score. Clinician feedback was used to improve the clinical plausibility and the usability of the final score (Emergency Department Assessment of Chest pain Score [EDACS]). EDACS was combined with electrocardiogram results and troponin results at 0 and 2 h to develop an ADP (EDACS-ADP). The score and EDACS-ADP were validated and tested for reproducibility in separate cohorts of patients. Results In the derivation (n = 1974) and validation (n = 608) cohorts, the EDACS-ADP classified 42.2% (sensitivity 99.0%, specificity 49.9%) and 51.3% (sensitivity 100.0%, specificity 59.0%) as low risk of major adverse cardiac events, respectively. The intra-class correlation coefficient for categorisation of patients as low risk was 0.87. Conclusion The EDACS-ADP identified approximately half of the patients presenting to the ED with possible cardiac chest pain as having low risk of short-term major adverse cardiac events, with high sensitivity. This is a significant improvement on similar, previously reported protocols. The EDACS-ADP is reproducible and has the potential to make considerable cost reductions to health systems.

Patient reported adherence to hypertension treatment: A revalidation study

Background Adherence to hypertension management in patients with hypertension is known to influence their blood pressure control. It is important to measure patients’ adherence behaviours to assist with designing appropriate interventions to improve blood pressure control. Aims The purposes of this study were to use confirmatory factor analysis to revalidate the Therapeutic Adherence Scale for Hypertensive Patients (TASHP), and to calculate the cut-off score for classifying adherence behaviours into two groups: satisfactory and low adherence behaviours. Methods Systematic random sampling was used to recruit patients with hypertension in China. Demographic characteristics, the TASHP and blood pressure were collected. The psychometric tests of the TASHP included: construct validity, criteria-related validity, internal reliability, and split-half reliability. The area under the receiver operating characteristics curve and Youden index were used to identify the cut-off score of the TASHP for blood pressure control. Results This study involved 366 patients. Confirmatory factor analysis supported the four-component structure of the TASHP proposed in the original scale development study. The TASHP has a satisfactory internal reliability (Cronbach’s α > 0.7) and a satisfactory split-half reliability (Spearman–Brown coefficients > 0.7). The patients with overall scores of the TASHP ⩾ 109 points were considered to have satisfactory adherence behaviours. Conclusion The TASHP is a validated and reliable instrument to measure the adherence to hypertension management in Chinese patients with hypertension. The cut-off score of 109 points can be considered as an effective measure to classify the level of adherence into satisfactory and low adherence behaviours.

Postnatal depressive symptoms amongst women in Central Vietnam: A cross-sectional study investigating prevalence and associations with social, cultural and infant factors

Background This study investigated the prevalence and socio-cultural correlates of postnatal mood disturbance amongst women 18–45 years old in Central Vietnam. Son preference and traditional confinement practices were explored as well as factors such as poverty, parity, family and intimate partner relationships and infant health. Methods A cross-sectional study was conducted in twelve randomly selected Commune Health Centres from urban and rural districts of Thua Thien Hue Province, Vietnam. Mother-infant dyads one to six months postpartum were invited to participate. Questionnaires from 431 mothers (urban n = 216; rural n = 215) assessed demographic and family characteristics, traditional confinement practices, son preference, infant health and social capital. The Edinburgh Postnatal Depression Scale (EPDS) and WHO5 Wellbeing Index indicated depressive symptoms and emotional wellbeing. Data were analysed using general linear models. Results Using an EPDS cut-off of 12/13, 18.1 % (n = 78, 95 % CI 14.6 - 22.1) of women had depressive symptoms (20.4 % urban; 15.8 % rural). Contrary to predictions, infant gender and traditional confinement were unrelated to depressive symptoms. Poverty, food insecurity, being frightened of family members, and intimate partner violence increased both depressive symptoms and lowered wellbeing. The first model accounted for 30.2 % of the variance in EPDS score and found being frightened of one’s husband, husband’s unemployment, breastfeeding difficulties, infant diarrhoea, and cognitive social capital were associated with higher EPDS scores. The second model had accounted for 22 % of the variance in WHO5 score. Living in Hue city, low education, poor maternal competence and a negative family response to the baby lowered maternal wellbeing. Conclusions Traditional confinement practices and son preference were not linked to depressive symptoms among mothers, but were correlates of family relationships and wellbeing. Poverty, food insecurity, violence, infant ill health, and discordant intimate and family relationships were linked with depressive symptoms in Central Vietnam.

Bone health in children with inflammatory bowel disease : adjusting for bone age

OBJECTIVES: Clinical results of bone mineral density for children with inflammatory bowel disease are commonly reported using reference data for chronological age. It is known that these children, particularly those with Crohn disease, experience delayed growth and maturation. Therefore, it is more appropriate to compare clinical results with bone age rather than chronological age. MATERIALS AND METHODS: Areal bone mineral density (aBMD) was measured using dual energy x-ray absorptiometry, and bone age was assessed using the Tanner-Whitehouse 3 method from a standard hand/wrist radiograph. Results were available for 44 children ages 7.99 to 16.89 years. Areal bone mineral density measurements were converted to z scores using both chronological and bone ages for each subject. RESULTS: Areal bone mineral density z scores calculated using bone age, as opposed to chronological age, were significantly improved for both the total body and lumbar spine regions of interest. When subjects were grouped according to diagnosis, bone age generated z scores remained significantly improved for those with Crohn disease but not for those diagnosed with ulcerative colitis. Grouping of children with Crohn disease into younger and older ages produced significantly higher z scores using bone age compared with chronological for the older age group, but not the younger age group. CONCLUSIONS: Our findings, in accordance with those presented in the literature, suggest that aBMD results in children with Crohn disease should include the consideration of bone age, rather than merely chronological age. Bone size, although not as easily available, would also be an important consideration for interpreting results in paediatric populations. © 2009 by European Society for Pediatric Gastroenterology, Hepatology, and Nutrition and North American Society for Pediatric Gastroenterology, Hepatology, and Nutrition.

A new regulatory variant in the interleukin-6 receptor gene associates with asthma risk

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ∼20 ng ml -1 and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml -1. Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16 705 asthmatics and 30 809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.

Growth patterns in preterm infants born appropriate for gestational age

Aim: This study aimed to document the growth patterns of a contemporary cohort of preterm infants born appropriate for gestational age (AGA). It was hypothesised that preterm AGA (PT-AGA) infants would display poorer growth than full-term AGA (FT-AGA) infants. Methods: Sixty-four PT-AGA infants and 64 FT-AGA infants were assessed at 0, 4, 8 and 12 months of corrected age (CA). Measurements of weight and length were recorded at each of the specified ages. Centers for Disease Control and Prevention growth data were used to calculate Z-scores for weight and length based on CA. Results: The mean length and weight Z-scores of PT-AGA infants were found to be significantly less than those of FT-AGA infants at term, 4, 8 and 12 months of CA (P < 0.001). The mean weight Z-score of PT-AGA infants was found to be less than their mean length Z-score at each time point, though the differences were not significant. Conclusions: The results of this study suggest that PT-AGA infants are likely to display poorer growth than FT-AGA infants until at least 1 year of CA. Long-term growth monitoring in this population is recommended. © 2008 The Authors.

Identification of IL6R and chromosome 11q13.5 as risk loci for asthma

We aimed to identify novel genetic variants affecting asthma risk, since these might provide novel insights into molecular mechanisms underlying the disease. We did a genome-wide association study (GWAS) in 2669 physician-diagnosed asthmatics and 4528 controls from Australia. Seven loci were prioritised for replication after combining our results with those from the GABRIEL consortium (n=26 475), and these were tested in an additional 25 358 independent samples from four in-silico cohorts. Quantitative multi-marker scores of genetic load were constructed on the basis of results from the GABRIEL study and tested for association with asthma in our Australian GWAS dataset. Two loci were confirmed to associate with asthma risk in the replication cohorts and reached genome-wide significance in the combined analysis of all available studies (n=57 800): rs4129267 (OR 1·09, combined p= 2·4×10-8) in the interleukin-6 receptor (IL6R) gene and rs7130588 (OR 1·09, p=1·8×10-8) on chromosome 11q13.5 near the leucine-rich repeat containing 32 gene (LRRC32, also known as GARP). The 11q13.5 locus was significantly associated with atopic status among asthmatics (OR 1·33, p=7×10-4), suggesting that it is a risk factor for allergic but not non-allergic asthma. Multi-marker association results are consistent with a highly polygenic contribution to asthma risk, including loci with weak effects that might be shared with other immune-related diseases, such as NDFIP1, HLA-B, LPP, and BACH2. The IL6R association further supports the hypothesis that cytokine signalling dysregulation affects asthma risk, and raises the possibility that an IL6R antagonist (tocilizumab) may be effective to treat the disease, perhaps in a genotype-dependent manner. Results for the 11q13.5 locus suggest that it directly increases the risk of allergic sensitisation which, in turn, increases the risk of subsequent development of asthma. Larger or more functionally focused studies are needed to characterise the many loci with modest effects that remain to be identified for asthma. National Health and Medical Research Council of Australia. A full list of funding sources is provided in the webappendix. © 2011 Elsevier Ltd.

Genetic analyses in a sample of individuals with high or low BMD shows association with multiple Wnt pathway genes

INTRODUCTION Although the high heritability of BMD variation has long been established, few genes have been conclusively shown to affect the variation of BMD in the general population. Extreme truncate selection has been proposed as a more powerful alternative to unselected cohort designs in quantitative trait association studies. We sought to test these theoretical predictions in studies of the bone densitometry measures BMD, BMC, and femoral neck area, by investigating their association with members of the Wnt pathway, some of which have previously been shown to be associated with BMD in much larger cohorts, in a moderate-sized extreme truncate selected cohort (absolute value BMD Z-scores = 1.5-4.0; n = 344). MATERIALS AND METHODS Ninety-six tag-single nucleotide polymorphism (SNPs) lying in 13 Wnt signaling pathway genes were selected to tag common genetic variation (minor allele frequency [MAF] > 5% with an r(2) > 0.8) within 5 kb of all exons of 13 Wnt signaling pathway genes. The genes studied included LRP1, LRP5, LRP6, Wnt3a, Wnt7b, Wnt10b, SFRP1, SFRP2, DKK1, DKK2, FZD7, WISP3, and SOST. Three hundred forty-four cases with either high or low BMD were genotyped by Illumina Goldengate microarray SNP genotyping methods. Association was tested either by Cochrane-Armitage test for dichotomous variables or by linear regression for quantitative traits. RESULTS Strong association was shown with LRP5, polymorphisms of which have previously been shown to influence total hip BMD (minimum p = 0.0006). In addition, polymorphisms of the Wnt antagonist, SFRP1, were significantly associated with BMD and BMC (minimum p = 0.00042). Previously reported associations of LRP1, LRP6, and SOST with BMD were confirmed. Two other Wnt pathway genes, Wnt3a and DKK2, also showed nominal association with BMD. CONCLUSIONS This study shows that polymorphisms of multiple members of the Wnt pathway are associated with BMD variation. Furthermore, this study shows in a practical trial that study designs involving extreme truncate selection and moderate sample sizes can robustly identify genes of relevant effect sizes involved in BMD variation in the general population. This has implications for the design of future genome-wide studies of quantitative bone phenotypes relevant to osteoporosis.

Doctor’s knowledge and practices of traumatic brain injury management in Chinese prehospital settings

Objectives: The incidence and mortality of traumatic brain injury (TBI) has increased rapidly in the last decade in China. Appropriate ambulance service can reduce case-fatality rates of TBI significantly. This study aimed to explore the factors (age, gender, education level, clinical experience, professional title, organization, specialty before prehospital care, and training frequency) that could influence prehospital doctors’ knowledge level and practices in TBI management in China, Hubei Province. Methods: A cross-sectional questionnaire survey was conducted in two cities in Hubei Province. The self-administered questionnaire consisted of demographic information and questions about prehospital TBI management. Independent samples t-test and one-way ANOVA were used to analyze group differences in the average scores in terms of demographic character. General linear regression was used to explore associated factors in prehospital TBI management. Results: A total of 56 questionnaires were handed out and 52 (93%) were returned. Participants received the lowest scores in TBI treatment (0.64; SD=0.08) and the highest scores in TBI assessment (0.80; SD=0.14). According to the regression model, the education level was associated positively with the score of TBI identification (P=.019); participants who worked in the emergency department (ED; P=.011) or formerly practiced internal medicine (P=.009) tended to get lower scores in TBI assessment; participants’ scores in TBI treatment were associated positively with the training frequency (P=.011); and no statistically significant associated factor was found in the overall TBI management. Conclusion: This study described the current situation of prehospital TBI management. The prehospital doctors’ knowledge level and practices in TBI management were quantified and the influential factors hidden underneath were explored. The results indicated that an appropriate continuing medical education (CME) program enables improvement of the quality of ambulance service in China.

Influence of LRP5 polymorphisms on normal variation in BMD

Genetic studies based on cohorts with rare and extreme bone phenotypes have shown that the LRP5 gene is an important genetic modulator of BMD. Using family-based and case-control approaches, this study examines the role of the LRP5 gene in determining normal population variation of BMD and describes significant association and suggestive linkage between LRP5 gene polymorphisms and BMD in >900 individuals with a broad range of BMD. Introduction: Osteoporosis is a common, highly heritable condition determined by complex interactions of genetic and environmental etiologies. Genetic factors alone can account for 50-80% of the interindividual variation in BMD. Mutations in the LRP5 gene on chromosome 11q12-13 have been associated with rare syndromes characterized by extremely low or high BMD, but little is known about the contribution of this gene to the development of osteoporosis and determination of BMD in a normal population. Materials and Methods: To examine the entire spectrum of low to high BMD, 152 osteoporotic probands, their families (597 individuals), and 160 women with elevated BMD (T score > 2.5) were recruited. BMD at the lumbar spine, femoral neck, and hip were measured in each subject using DXA. Results: PAGE sequencing of the LRP5 gene revealed 10 single nucleotide polymorphisms (SNPs), 8 of which had allele frequencies of >5%, in exons 8, 9, 10, 15, and 18 and in introns 6, 7, and 21. Within families, a strong association was observed between an SNP at nucleotide C171346A in intron 21 and total hip BMD (p < 1 × 10-5 in men only, p = 0.0019 in both men and women). This association was also observed in comparisons of osteoporotic probands and unrelated elevated BMD in women (p = 0.03), along with associations with markers in exons 8 (C135242T, p = 0.007) and 9 (C141759T, p = 0.02). Haplotypes composed of two to three of the SNPs G121513A, C135242T, G138351A, and C141759T were strongly associated with BMD when comparing osteoporotic probands and high BMD cases (p < 0.003). An SNP at nucleotide C165215T in exon 18 was linked to BMD at the lumbar spine, femoral neck, and total hip (parametric LOD scores = 2.8, 2.5, and 2.2 and nonparametric LOD scores = 0.3, 1.1, and 2.2, respectively) but was not genetically associated with BMD variation. Conclusion: These results show that common LRP5 polymorphisms contribute to the determination of BMD in the general population.

The effect of transforming growth factor β1 gene polymorphisms in ankylosing spondylitis

Objectives. To determine whether genetic polymorphisms in or near the transforming growth factor β1 (TGFB1) locus were associated d with susceptibility to or severity of ankylosing spondylitis (AS). Methods. Five intragenic single-nucleotide polymorphisms (SNP) and three microsatellite markers flanking the TGFB1 locus were genotyped. Seven hundred and sixty-two individuals from 184 multiplex families were genotyped for the microsatellite markers and two of the promoter SNPs. One thousand and two individuals from 212 English and 170 Finnish families with AS were genotyped for all five intragenic SNPs. A structured questionnaire was used to assess the age of symptom onset, disease duration and disease severity scores, including the BASDAI (Bath Ankylosing Spondylitis Disease Activity Index) and BASFI (Bath Ankylosing Spondylitis Functional Index). Results. A weak association was noted between the rare TGFB1 + 1632 T allele and AS in the Finnish population (P = 0.04) and in the combined data set (P = 0.03). No association was noted between any other SNPs or SNP haplotype and AS, even among those families with positive non-parametric linkage scores. The TGFB1 +1632 polymorphism was also associated with a younger age of symptom onset (English population, allele 2 associated with age of onset greater by 4.2 yr, P = 0.05; combined data set, allele 2 associated with age of onset greater by 3.2 yr, P = 0.02). A haplotype of coding region SNPs (TGFB1 +869/ +915+1632 alleles 2/1/2) was associated with age of symptom onset in both the English parent-case trios and the combined data set (English data set, haplotype 2/1/2 associated with age of onset greater by 4.9 yr, P = 0.03; combined data set, haplotype 2/1/2 associated with greater age of onset by 4.2 yr, P = 0.006). Weak linkage with AS susceptibility was noted and the peak LOD score was 1.3 at distance 2 cM centromeric to the TGFB1 gene. No other linkage or association was found between quantitative traits and the markers. Conclusion. This study suggests that the polymorphisms within the TGFB1 gene play at most a small role in AS and that other genes encoded on chromosome 19 are involved in susceptibility to the disease.

Concordance of disease severity among family members with ankylosing spondylitis?

Objective. The heritability of disease activity and function in ankylosing spondylitis (AS) have been estimated at 0.51 and 0.63 (i.e., 51% and 63%), respectively. We examined the concordance of disease severity among family members in terms of disease activity, function, radiological change, prevalence of iritis, and juvenile onset. Methods. Disease activity and functional impairment due to AS were studied using the Bath AS Disease Activity Index (BASDAI) and Functional Index (BASFI) self-administered questionnaires; radiographic involvement was measured using the Bath AS Radiology Index (BASRI) scale. Familial correlation of BASDAI and BASFI was assessed in 406 families with 2 or more cases, using the program PAP. Parent-child and sibling-sibling concordance for iritis and juvenile AS were also studied in these families. Heritability of radiological disease severity based on the BASRI was assessed in 29 families containing 60 affected individuals using the program SOLAR. Results. Correlations between parent-child pairs for disease activity and function were 0.07 for both. Correlations between sibling pairs for disease activity and function were 0.27 and 0.36, respectively. The children of AS parents with iritis were more likely to develop iritis [27/71 (38%)] than children of non-iritis AS parents [13/70 (19%)] (p = 0.01). Parents with JAS were more likely to have children with JAS [17/30 (57%) compared to non-JAS parents 34/111 (30%)] (p = 0.002). The heritability of radiological disease severity based on the BASRI was 0.62. Conclusion. While correlation in severity between parent and child is poor, siblings do resemble each other in terms of severity, supporting the findings of segregation studies indicating significant genetic dominance in the heritable component of disease activity. Significant parent-child concordance for iritis and juvenile disease onset suggest that there are genetic risk factors for these traits independent of those determining the risk of AS itself. The finding of significant heritability of radiological change (BASRI) provides support using an objective measure for the observed heritability of the questionnaire-assessed disease severity scores, ASDAI and BASFI.

Suggestive linkage of the parathyroid receptor type 1 to osteoporosis

We have investigated the role of 23 candidate genes in the control of bone mineral density (BMD) by linkage studies in families of probands with osteoporosis (lumbar spine [LS] or femoral neck [FN] BMD T score < -2.5) and low BMD relative to an age- and gender-matched cohort (Z score < -2.0). One hundred and fifteen probands (35 male, 80 female) and 499 of their first- or second-degree relatives (223 males and 276 females) were recruited for the study. BMD was measured at the LS and FN using dual-energy X-ray absorptiometry and expressed as age- and gender-matched Z scores corrected for body mass index. The candidate genes studied were the androgen receptor, type I collagen A1 (COLIA1), COLIA2, COLIIA1, vitamin D receptor (VDR), colony-stimulating factor 1, calcium-sensing receptor, epidermal growth factor (EGF), estrogen receptor 1 (ESR1), fibrillin type 1, insulin-like growth factor 1, interleukin-1 alpha (IL-1α), interleukin-4 (IL-4), interleukin-6 (IL-6), interleukin-11 (IL-11), osteopontin, parathyroid hormone (PTH), PTH-related peptide, PTH receptor type 1 (PTHR1), transforming growth factor-beta 1, and tumor necrosis factors alpha and beta. Sixty-four microsatellites lying close to or within these genes were investigated for linkage with BMD. Using the program MapMaker/Sibs there was suggestive evidence of linkage between BMD and PTHR1 (maximum LOD score obtained [MLS] 2.7-3.5). Moderate evidence of linkage was also observed with EGF (MLS 1.8), COLIA1 (MLS 1.7), COLIIA1/VDR (MLS 1.7), ESR1 (MLS 1.4), IL-1α (MLS 1.4), IL-4 (MLS 1.2), and IL-6 (MLS 1.2). Variance components analysis using the program ACT, correcting for proband-wise ascertainment, also showed evidence of linkage (p ≤0.05) at markers close to or within the candidate genes IL- 1α, PTHR1, IL-6, and COLIIA1/VDR. Further studies will be required to confirm these findings, to refine the location of gene responsible for the observed linkage, and to screen the candidate genes targeted at these loci for mutations.