Statistical imputation of classical human leukocyte antigen (HLA) alleles is becoming an indispensable tool for fine-mappings of disease association signals from case-control genome-wide association ...studies. However, most currently available HLA imputation tools are based on European reference populations and are not suitable for direct application to non-European populations. Among the HLA imputation tools, The HIBAG R package is a flexible HLA imputation tool that is equipped with a wide range of population-based classifiers; moreover, HIBAG R enables individual researchers to build custom classifiers. Here, two data sets, each comprising data from healthy Japanese individuals of difference sample sizes, were used to build custom classifiers. HLA imputation accuracy in five HLA classes (HLA-A, HLA-B, HLA-DRB1, HLA-DQB1 and HLA-DPB1) increased from the 82.5-98.8% obtained with the original HIBAG references to 95.2-99.5% with our custom classifiers. A call threshold (CT) of 0.4 is recommended for our Japanese classifiers; in contrast, HIBAG references recommend a CT of 0.5. Finally, our classifiers could be used to identify the risk haplotypes for Japanese narcolepsy with cataplexy, HLA-DRB1*15:01 and HLA-DQB1*06:02, with 100% and 99.7% accuracy, respectively; therefore, these classifiers can be used to supplement the current lack of HLA genotyping data in widely available genome-wide association study data sets.
Renal hypouricemia is an inherited disorder characterized by impaired tubular uric acid transport. Impairment of the function of URAT1, the main transporter for the reabsorption of uric acid at the ...apical membrane of the renal tubules, causes renal hypouricemia. The G774A mutation in the SLC22A12 gene encoding URAT1 predominates in Japanese renal hypouricemia. From data on linkage disequilibrium between the G774 locus and the 13 markers flanking it (12 single nucleotide polymorphisms and 1 dinucleotide insertion/deletion locus), we here estimate the age of this mutation at approximately 6820 years 95% confidence interval (CI) 1860–11,760 years; median = 2460 years. This indicates that the origin of the G774A mutation dates back from between the time when the Jomon people predominated in Japan and the time when the Yayoi people started to migrate to Japan from the Korean peninsula. These data are consistent with a recent finding that this G774A mutation was also predominant in Koreans with hypouricemia and indicate that the mutation originated on the Asian continent. Thus, this mutation found in Japanese patients was originally brought by immigrant(s) from the continent and thereafter expanded in the Japanese population either by founder effects or by genetic drift (or both).
Background: An imbalance in cytokine homoeostasis is thought to have a key role in the neuropsychiatric syndromes of systemic lupus erythematosus (NPSLE), and recently, a role for chemokines has been ...noted. Objective: To compare concentrations of monocyte chemotactic protein-1 (MCP-1)/CCL2 in cerebral spinal fluid (CSF) of patients with SLE, and with and without neuropsychiatric symptoms. Methods: CSF was obtained from 185 patients with SLE: 96 with NPSLE and 89 patients with SLE without neuropsychiatric symptoms (non-NPSLE patients). MCP-1/CCL2 concentrations were measured with an ELISA. Results: The average concentration of CSF MCP-1/CCL2 in patients with NPSLE was 1959 pg/ml, and in non-NPSLE patients 712 pg/ml. The average MCP-1/CCL2 concentration was significantly higher in the NPSLE group than in the non-NPSLE group (p<0.001). In one representative patient with NPSLE, MCP-1/CCL2 levels in the CSF decreased in parallel with a decline in neuropsychiatric symptoms. Conclusions: CSF MCP-1/CCL2 levels are higher in patients with NPSLE than in non-NPSLE patients. MCP-1/CCL2 may have an important role in the expression of NPSLE. These results indicate that CSF MCP-1/CCL2 reflects an inflammatory activity in the brain, suggesting that it might be used as a diagnostic tool and a monitor for therapeutic responses in patients with NPSLE.
Malignant hypertension is a serious form of arterial hypertension in which the physiopathological mechanisms include increased activity of the sympathetic nervous system, renin angiotensin system, ...and endothelium dysfunction. Family history of hypertension is an important predictive factor for hypertension and is associated with metabolic and hemodynamic abnormalities. Studies of these abnormalities in malignant hypertensive offspring have not yet been published. Therefore, we studied 42 offspring of malignant hypertensive parents (OMH group: age, 22+/-7 years; 23 male subjects; 27 white) and 35 offspring of normotensive parents (ONT group: age, 21+/-4 years; 23 male subjects; 25 white). All subjects had blood pressure <140/90 mm Hg. We evaluated body mass index; office blood pressure; 24-hour ambulatory and continuous beat-to-beat blood pressure monitoring (Finapres); biochemical analysis, including total cholesterol and fractions, triglycerides, glucose, and insulin; and hormonal analysis, including plasma renin activity, aldosterone, and catecholamines. The subjects were also submitted to cold pressure test and handgrip measurements. The body mass index was significantly higher in the OMH group (24+/-5 kg/m(2)) than in the ONT group (22+/-4 kg/m(2)). The OMH group showed significantly higher blood pressure and heart rate in office and Finapres measurements (P<0.05). In 24-hour ambulatory monitoring, the OMH group presented higher 24-hour blood pressure and heart rate, higher blood pressure during the night, and higher heart rate variability during the day compared with those of the ONT group. They also presented lower HDL cholesterol, higher levels of plasma insulin and norepinephrine, and higher insulin-to-glucose ratio (P<0.05) than the ONT group. There were no differences in the other biochemical parameters measured. In conclusion, OMH subjects show early hemodynamic, neurohumoral, and metabolic alterations that are typical of hypertensive metabolic syndrome.
An algorithm for estimating haplotypes associated with several quantitative phenotypes is proposed. The concept of a receiver operating characteristic (ROC) curve was introduced, and a linear ...combination of the quantitative phenotypic values was considered. This set of values was divided into two parts: values for subjects with and without a particular haplotype. The goodness of its partition was evaluated by the area under the ROC curve (AUC). The AUC value varied from 0 to 1; this value was close to 1 when the partition had high accuracy. Therefore, the strength of association between phenotypes and haplotypes was considered to be proportional to the AUC value. In our algorithm, the parameters representing a degree of association between the haplotypes and phenotypes were estimated so as to maximize the AUC value; further, the haplotype with the maximum AUC value was considered to be the best haplotype associated with the phenotypes. This algorithm was implemented by using R language. The effectiveness of our algorithm was evaluated by applying it to real genotype data of the Calpine-10 gene obtained from diabetics. The results showed that our algorithm was more reasonable and advantageous for use with several quantitative phenotypes than the generalized linear model or the neural network model.
Objective
To determine the most sensitive scoring method for assessment of rheumatoid arthritis (RA) disease activity using the American College of Rheumatology Core Data Set.
Methods
The subjects ...were 4,530 patients with RA (mean age 57.9 years, mean disease duration 12.7 years) who participated in a large observational cohort study of RA patients. The 68 joints assessed were classified into 15 joint areas, and each joint variable was categorized based on the presence or absence of swelling or pain in these areas. Multiple linear regression and analysis of variance were used to evaluate the significance of effects of these 15 joint areas on variables for assessment of RA disease activity such as patient's assessment of pain on a visual analog scale (VAS), patient's and physician's global assessment of disease activity on a VAS, HAQ (Health Assessment Questionnaire), and Japanese HAQ.
Results
Although the 3 most frequently affected joints were the wrist, metacarpophalangeal joints, and proximal interphalangeal joints, the 5 joints with the largest contributions to all of the variables assessed for disease activity were the shoulder, elbow, and knee joints, followed by the wrist and ankle joints. The combination of shoulder, elbow, and knee joints accounted for approximately 70% of the contribution to all the variables, while addition of the wrist and ankle joints increased this value to approximately 90%.
Conclusion
Scoring for assessment of RA disease activity would be more sensitive if separate joints such as the shoulder, elbow, knee, wrist, and ankle joints were weighted differently.
Individual disease risk estimated based on the data from single or multiple genetic loci is generally calculated using the genotypes of a subject, frequencies of alleles of interest, odds ratios and ...the average population risk. However, it is often difficult to estimate accurately the average population risk, and therefore it is often expressed as an interval. To better estimate the risk of a subject with given genotypes, we built R scripts using the R environment and constructed graphs to examine the change in the estimated risk as well as the relative risk according to the change of the average population risk. In most cases, the graph of the relative risk did not cross the line of y=1, thereby indicating that the order of the relative risk for given genotypes and the population average risk does not change when the average risk increases or decreases. In rare cases, however, the graph of the relative risk crossed the line of y=1, thereby indicating that the order of the relative risk for given genotypes and the population average risk does change owing to the change in the population risk. We propose that the relative risk should be estimated for not only the point average population risk but also for an interval of the average population risk. Moreover, when the graph crosses the line of y=1 within the interval, this information should be reported to the consumer.