Isolation of porcine circovirus 3 using primary porcine bone marrow-derived cells.

Porcine circovirus 3 (PCV3) was first reported in the United States in 2016; this virus is considered to be involved in diverse pathologies, such as multisystem inflammation, porcine dermatitis and nephropathy syndrome, and reproductive disorders. However, successful isolation of PCV3 using cultured cells has been rare. In this study, we aimed to isolate PCV3 using primary porcine bone marrow-derived cells. Mononuclear cells were isolated from the femur bones of clinically healthy pigs. These primary cells were cultured for 6-10 days post-seeding and infected with PCV3-containing tissue homogenates. The cells were cultured for up to 37 days, and the culture medium was changed every 3-4 days. The growth curve of PCV3 in porcine bone marrow cells revealed a decline in growth during the first 10 days post-infection, followed by an increase leading to > 1010 genomic copies/mL of the cell culture supernatant; moreover, the virus was capable of passaging. The indirect fluorescent antibody assay for PCV3 infection revealed the presence of PCV3 capsid protein in the cytoplasm and nuclei of infected cells. Bone marrow cells were passaged for more than 20 generations (over 5 months), and PCV3 persistently infected the cells. PCV3-infected bone marrow cells expressed mesenchymal markers. These results reflect that primary porcine bone marrow-derived mesenchymal cells are permissive to PCV3 and continuously replicate a high copy number of the PCV3 genome. These findings regarding the high replication rate of PCV3 in bone marrow-derived mesenchymal cells could enhance our understanding of PCV3 pathogenicity.

Expression of RNautophagy/DNautophagy-related genes is regulated under control of an innate immune receptor.

Double-stranded RNA (dsRNA) is a molecular pattern uniquely produced in cells infected with various viruses as a product or byproduct of replication. Cells detect such molecules, which indicate non-self invasion, and induce diverse immune responses to eliminate them. The degradation of virus-derived molecules can also play a role in the removal of pathogens and suppression of their replication. RNautophagy and DNautophagy are cellular degradative pathways in which RNA and DNA are directly imported into a hydrolytic organelle, the lysosome. Two lysosomal membrane proteins, SIDT2 and LAMP2C, mediate nucleic acid uptake via this pathway. Here, we showed that the expression of both SIDT2 and LAMP2C is selectively upregulated during the intracellular detection of poly(I:C), a synthetic analog of dsRNA that mimics viral infection. The upregulation of these two gene products upon poly(I:C) introduction was transient and synchronized. We also observed that the induction of SIDT2 and LAMP2C expression by poly(I:C) was dependent on MDA5, a cytoplasmic innate immune receptor that directly recognizes poly(I:C) and induces various antiviral responses. Finally, we showed that lysosomes can target viral RNA for degradation via RNautophagy and may suppress viral replication. Our results revealed a novel degradative pathway in cells as a downstream component of the innate immune response and provided evidence suggesting that the degradation of viral nucleic acids via RNautophagy/DNautophagy contributes to the suppression of viral replication.

Immunological features of a lung cancer patient achieving an objective response with anti-programmed death-1 blockade therapy.

The role of immune checkpoint inhibitors in metastatic lung cancer has been established in recent years and the pretherapeutic profiles of the tumor microenvironment in responders have been increasingly reported. The role of salvage surgery and the immune profiles of the posttherapeutic specimens in patients achieving an objective response have rarely been studied. We report a case of metastatic lung cancer treated by anti-programmed death-1 Ab followed by surgical resection. The immune status of the tumor was assessed, showing germinal center formation, memory B cell infiltration, and a high frequency of interferon gamma -secreting T cells.

An integrative network-based approach to identify novel disease genes and pathways: a case study in the context of inflammatory bowel disease.

BACKGROUND: There are different and complicated associations between genes and diseases. Finding the causal associations between genes and specific diseases is still challenging. In this work we present a method to predict novel associations of genes and pathways with inflammatory bowel disease (IBD) by integrating information of differential gene expression, protein-protein interaction and known disease genes related to IBD. RESULTS: We downloaded IBD gene expression data from NCBI's Gene Expression Omnibus, performed statistical analysis to determine differentially expressed genes, collected known IBD genes from DisGeNet database, which were used to construct a IBD related PPI network with HIPPIE database. We adapted our graph-based clustering algorithm DPClusO to cluster the disease PPI network. We evaluated the statistical significance of the identified clusters in the context of determining the richness of IBD genes using Fisher's exact test and predicted novel genes related to IBD. We showed 93.8% of our predictions are correct in the context of other databases and published literatures related to IBD. CONCLUSIONS: Finding disease-causing genes is necessary for developing drugs with synergistic effect targeting many genes simultaneously. Here we present an approach to identify novel disease genes and pathways and discuss our approach in the context of IBD. The approach can be generalized to find disease-associated genes for other diseases.

Purification and characterization of immunogenic recombinant virus-like particles of porcine circovirus type 2 expressed in silkworm pupae.

Porcine circovirus type 2 (PCV2) is a primary causative agent of postweaningmultisystemic wasting syndrome (PMWS), which has a significant economic impact on the swine industry. The capsid protein (Cap) encoded by ORF2 of the viral genome has been used effectively as a vaccine against PCV2 infection. The Cap protein can spontaneously assemble into virus-like particles (VLPs) that are safe and highly immunogenic for vaccine applications. Several expression systems, including bacteria, yeast and insect cells, have been utilized to produce PCV2 VLPs. However, in some cases, the recombinant Cap (rCap) proteins produced in bacteria and yeast do not assemble spontaneously. In this study, we expressed rCap protein using a silkworm-baculovirus expression vector system (silkworm-BEVS) for mass production of PCV2 VLPs and established a simple three-step protocol for its purification from pupae: extraction by detergent, ammonium sulfate precipitation and anion exchange column chromatography. Size-exclusion chromatography (SEC) analysis and transmission electron microscope (TEM) observation showed that purified rCap proteins formed VLPs with a similar morphology to that of the original virus. Furthermore, the VLPs produced in silkworms were capable of inducing neutralizing antibodies against PCV2 in mice. Our results demonstrated that the silkworm system is a powerful tool for the production of PCV2 VLPs and will be useful for the development of a reliable and cost-effective PCV2 vaccine.

Efficacy of genogroup 1 based porcine epidemic diarrhea live vaccine against genogroup 2 field strain in Japan.

BACKGROUND: Porcine epidemic diarrhea (PED) is a lethal infectious disease in suckling piglets with symptoms including watery diarrhea caused by PED virus (PEDV). Since the late 1990's, live vaccines based on genogroup 1 virus have been used in Japan, and a significant amount of the vaccine has been used even after new genogroups invaded in 2013. In this study, we evaluated the effect of a conventional PED live vaccine on a newly prevalent genogroup 2 field strain in experimental and field situations. METHODS: Two pregnant sows were administered twice the live vaccine before farrowing. A pregnant sow was served as a negative control. All newborn piglets were challenged with the genogroup 2 virus, and clinical signs were monitored for 7 days post challenge. PEDV-specific immune responses in serum and milk of the sows were assayed by virus neutralization assay. The efficacy of PED live vaccine in vaccinated or non-vaccinated farms was evaluated by comparing the mortality rate of suckling piglets after the onset of PED. RESULTS: The challenged piglets exhibited watery diarrhea with or without vaccination. However, the clinical score of piglets born from vaccinated sows significantly improved after the 4th day of the challenge. The survival rate of piglets in the vaccinated group at the end of the experimental period was 80%, whereas in the control group was 0%. Neutralizing antibody titers in serum and milk of control sow was negative throughout the experimental period, whereas high titers were observed in the vaccinated sows. The vaccinated farms significantly reduced the mortality rate of suckling piglets after the onset of PED, compared to farms not vaccinated. CONCLUSIONS: The conventional PED live vaccine induced the lactogenic immunity to vaccinated sows and showed partial protection against the genogroup 2 virus both under the experimental and field conditions.

Successful removal of porcine circovirus-1 from immunoglobulin G formulated in glycine solution using nanofiltration.

Porcine circovirus (PCV) is a potentially harmful virus that has been shown to contaminate biological products. The virus is resistant to many inactivation and/or removal procedures performed during manufacturing. Anion exchange chromatography has been shown to be useful for PCV type 1 (PCV1) removal; however, reduction of PCV1 using methods such as heat inactivation, low pH, nanofiltration, UV-C, and gamma irradiation has not been successful. Therefore, in this study, we evaluate various conditions for process solutions during nanofiltration using PCV1. The results indicated that PCV could be effectively removed from glycine solution at 0.1-0.3 M, pH 4.0 without IgG, using a nanofilter with a pore size of 19 nm (19-nm filter); log reduction values (LRVs) of ≥4.5 and ≥ 5.0, respectively, were obtained. In contrast, PCV1 was significantly removed (LRV: 2.2) in glycine solution at 0.3 M, pH 6.0 with 1.0% IgG using the 19-nm filter, but some virus genomes were detected in the filtrates. In summary, the use of a 19-nm filter in glycine solution with/without IgG is an appropriate condition for PCV removal.

Integrated pathway-based transcription regulation network mining and visualization based on gene expression profiles.

Conventionally, workflows examining transcription regulation networks from gene expression data involve distinct analytical steps. There is a need for pipelines that unify data mining and inference deduction into a singular framework to enhance interpretation and hypotheses generation. We propose a workflow that merges network construction with gene expression data mining focusing on regulation processes in the context of transcription factor driven gene regulation. The pipeline implements pathway-based modularization of expression profiles into functional units to improve biological interpretation. The integrated workflow was implemented as a web application software (TransReguloNet) with functions that enable pathway visualization and comparison of transcription factor activity between sample conditions defined in the experimental design. The pipeline merges differential expression, network construction, pathway-based abstraction, clustering and visualization. The framework was applied in analysis of actual expression datasets related to lung, breast and prostrate cancer.

Targeted Integration of RNA-Seq and Metabolite Data to Elucidate Curcuminoid Biosynthesis in Four Curcuma Species.

Curcuminoids, namely curcumin and its analogs, are secondary metabolites that act as the primary active constituents of turmeric (Curcuma longa). The contents of these curcuminoids vary among species in the genus Curcuma. For this reason, we compared two wild strains and two cultivars to understand the differences in the synthesis of curcuminoids. Because the fluxes of metabolic reactions depend on the amounts of their substrate and the activity of the catalysts, we analyzed the metabolite concentrations and gene expression of related enzymes. We developed a method based on RNA sequencing (RNA-Seq) analysis that focuses on a specific set of genes to detect expression differences between species in detail. We developed a 'selection-first' method for RNA-Seq analysis in which short reads are mapped to selected enzymes in the target biosynthetic pathways in order to reduce the effect of mapping errors. Using this method, we found that the difference in the contents of curcuminoids among the species, as measured by gas chromatography-mass spectrometry, could be explained by the changes in the expression of genes encoding diketide-CoA synthase, and curcumin synthase at the branching point of the curcuminoid biosynthesis pathway.

KNApSAcK Metabolite Activity Database for retrieving the relationships between metabolites and biological activities.

Databases (DBs) are required by various omics fields because the volume of molecular biology data is increasing rapidly. In this study, we provide instructions for users and describe the current status of our metabolite activity DB. To facilitate a comprehensive understanding of the interactions between the metabolites of organisms and the chemical-level contribution of metabolites to human health, we constructed a metabolite activity DB known as the KNApSAcK Metabolite Activity DB. It comprises 9,584 triplet relationships (metabolite-biological activity-target species), including 2,356 metabolites, 140 activity categories, 2,963 specific descriptions of biological activities and 778 target species. Approximately 46% of the activities described in the DB are related to chemical ecology, most of which are attributed to antimicrobial agents and plant growth regulators. The majority of the metabolites with antimicrobial activities are flavonoids and phenylpropanoids. The metabolites with plant growth regulatory effects include plant hormones. Over half of the DB contents are related to human health care and medicine. The five largest groups are toxins, anticancer agents, nervous system agents, cardiovascular agents and non-therapeutic agents, such as flavors and fragrances. The KNApSAcK Metabolite Activity DB is integrated within the KNApSAcK Family DBs to facilitate further systematized research in various omics fields, especially metabolomics, nutrigenomics and foodomics. The KNApSAcK Metabolite Activity DB could also be utilized for developing novel drugs and materials, as well as for identifying viable drug resources and other useful compounds.

Pilot study of incidence of gastroesophageal reflux after lung resection.

Patients undergoing lung resection may be at risk of gastroesophageal reflux (GER) and silent aspiration following surgery. Defining high-risk patients may lead to prevention strategies for silent aspiration and subsequent exacerbation of underlying pulmonary disease. A pilot study of 50 patients was performed to investigate postoperative gastroesophageal reflux disease (GERD) symptoms and the pepsin concentration in saliva. Patients answered a questionnaire concerning GERD symptoms before lung surgery and at the time of discharge. Saliva samples were obtained before surgery, on the third postoperative day and at discharge. Pepsin concentration was measured with Peptest. The pepsin concentration in saliva following resection was significantly elevated on postoperative day 3, but it returned to the baseline level at discharge. Patients undergoing resection of four or more lung subsegments had a continuously elevated pepsin concentration in saliva on postoperative day 3 [mean difference 65.63 ng/mL, 95% confidence interval (CI): 9.130-122.1] and at discharge (mean difference 76.22 ng/mL, 95% CI: 19.72-132.7). Patients with a >10% reduction of forced expiration volume in one second also had a continuous elevated pepsin concentration from the 3rd postoperative day. Lung resection resulted in elevated pepsin concentration in the saliva, which persisted in patients who received resections equivalent to or more than right middle lobectomy in volume. Resection of large volumes of lung may lead to anatomical changes and changes in breathing patterns and result in GER.

1,2,4-Triazol-3-ylidenes with an N-2,4-dinitrophenyl substituent as strongly π-accepting N-heterocyclic carbenes.

The synthesis and characterisation of a series of new Rh and Au complexes bearing 1,2,4-triazol-3-ylidenes with a N-2,4-dinitrophenyl (N-DNP) substituent are described. IR, NMR, single-crystal X-ray diffraction and computational analyses of the Rh complexes revealed that the N-heterocyclic carbenes (NHCs) behaved as strong π acceptors and weak σ donors. In particular, a natural bond orbital (NBO) analysis revealed that the contributions of the Rh→Ccarbene π backbonding interaction energies (ΔEbb ) to the bond dissociation energies (BDE) of the RhCcarbene bond for [RhCl(NHC)(cod)] (cod=1,5-cyclooctadiene) reached up to 63 %. The Au complex exhibited superior catalytic activity in the intermolecular hydroalkoxylation of cyclohexene with 2-methoxyethanol. The NBO analysis suggested that the high catalytic activity of the Au(I) complex resulted from the enhanced π acidity of the Au atom.

The involvement of annexin II in resistance to UVB-induced cell death and in the increased nucleotide excision repair capacity of UV-damaged DNA in human cells.

Annexin II, an HSP27-interacting protein, is involved in the protection of human cells against UVC. UVB is concerned with deleterious actions on human health. In this study, we attempted to confirm the anti-UVB effect of annexin II, and to elucidate the mechanisms underlying annexin II-involving UV resistance. The RSa cells were more sensitive to UVB lethality than the AP(r)-1 cells. Overproduction of annexin II in RSa cells resulted in increased resistance to UVB lethality, while annexin II siRNA-transfected AP(r)-1 cells were sensitized to UVB lethality. The excision capacity of the two major types (CPD and 6-4PP) of UVC- and UVB-damaged DNA in AP(r)-1 cells was greater than in RSa cells. The excision capacity of the RSa cells improved following upregulation of annexin II, while the capacity of the AP(r)-1 cells decreased after annexin II downregulation. Our results suggest that annexin II is involved in the UV resistance of human cells, via functioning in nucleotide excision repair.

[Susceptibilities of bacteria isolated from patients with lower respiratory infectious diseases to antibiotics (2006)].

From October 2006 to September 2007, we collected the specimen from 356 patients with lower respiratory tract infections in 14 institutions in Japan, and investigated the susceptibilities of isolated bacteria to various antibacterial agents and patients' characteristics. Of 414 strains that were isolated from specimen (mainly from sputum) and assumed to be bacteria causing in infection, 407 strains were examined. The isolated bacteria were: Staphylococcus aureus 64, Streptococcus pneumoniae 96, Haemophilus influenzae 87, Pseudomonas aeruginosa (non-mucoid) 52, P. aeruginosa (mucoid) 11, Klebsiella pneumoniae 20, and Moraxella catarrhalis 44. Of 64 S. aureus strains, those with 2 microg/ml or less of MIC of oxacillin (methicillin-susceptible S. aureus: MSSA) and those with 4 microg/ml or more of MIC of oxacillin (methicillin-resistant S. aureus: MRSA) were 27 (42.2%) and 37 (57.8%) strains, respectively. Against MSSA, imipenem had the most potent antibacterial activity and inhibited the growth of all strains at 0.063 microg/ml or less. Against MRSA, vancomycin and linezolid showed the most potent activity and inhibited the growth of all the strains at 1 microg/ml. Carbapenems showed the most potent activities against S. pneumoniae and in particular, panipenem inhibited the growth of all the strains at 0.063 microg/ml or less. Imipenem and faropenem also had a preferable activity and inhibited the growth of all the strains at 0.125 and 0.5 microg/ml, respectively. In contrast, there were high-resistant strains (MIC: over 128 microg/ml) for erythromycin (45.8%) and clindamycin (20.8%). Against H. influenzae, levofloxacin showed the most potent activity and its MIC90 was 0.063 microg/ml or less. Meropenem showed the most potent activity against P. aeruginosa (mucoid) and its MIC90 was 0.5 microg/ml. Against P. aeruginosa (non-mucoid), tobramycin had the most potent activity and its MIC90 was 2 microg/ml. Against K. pneumoniae, cefozopran was the most potent activity and inhibited the growth of all the strains at 0.063 microg/ml or less. Also, all the antibacterial agents except ampicillin generally showed a potent activity against M. catarrhalis and the MIC90 of them were 2 microg/ml or less. The approximately half the number (50.6%) of the patients with respiratory infection were aged 70 years or older. Bacterial pneumonia and chronic bronchitis accounted for 49.2% and 28.1% of all the respiratory infections, respectively. The bacteria frequently isolated from the patients with bacterial pneumonia were S. pneumoniae (29.2%), S. aureus (20.8%), and H. influenzae (12.9%). H. influenzae (25.0%) and P. aeruginosa (21.7%) also were frequently isolated from the patients with chronic bronchitis. Before the antibacterial agent administration, the bacteria frequently isolated from the patients were S. pneumoniae (27.5%) and H. influenzae (22.5%). The bacteria frequently isolated from the patients treated with macrolides was P. aeruginosa, and its isolation frequently was 39.4%.

Experimental inoculation of a tissue homogenate containing porcine circovirus type 3 obtained after two in vivo passages in NIBS miniature pigs.

Porcine circovirus type 3 (PCV3) is a novel porcine circovirus that has been detected in pigs showing various clinical and pathological conditions, as well as in many asymptomatic pigs. The pathogenesis of PCV3 infection in pigs remains unclear. To evaluate the in vivo growth and pathogenicity of PCV3, we performed two experiments on PCV3 infection in laboratory-grade miniature pigs with strictly controlled genetic backgrounds and microbiological status. A PCV3 passage experiment confirmed PCV3 genome detection in the sera and multiple organs via in vivo serial passage generations. PCV3 was successively passaged in miniature pigs by inoculating tissue homogenates from infected pigs supporting Koch's principles. In the PCV3 infection experiment, viremia was observed in all the inoculated pigs, and transient neurological signs were observed in one of the three pigs. Histopathologically, all three pigs in the PCV3 inoculation group exhibited lung disorders such as interstitial pneumonia and lymphoplasmacytic perivasculitis. In addition, one pig with neurological signs in the PCV3 inoculation group showed focal thrombosis in the meninges of the cerebellum. Vascular lesions in both the lungs and brain suggest that PCV3 may cause injury to vascular tissues. In situ hybridization (ISH)-RNA analysis demonstrated that the PCV3 genome was localized in the lymph nodes of pigs inoculated with PCV3. The PCV3 in vivo passage system in NIBS miniature pigs will help investigate the pathogenicity of PCV3.

Mutations in the spike gene of porcine epidemic diarrhea virus associated with growth adaptation in vitro and attenuation of virulence in vivo.

Previously, we have reported that a serial passage of 83P-5 strain of porcine epidemic diarrhea virus (PEDV) in Vero cells resulted in a growth adaptation of the virus in cultured cells at the 22nd passage. In this study, we further maintained the 83P-5 in Vero cells up to the 100th passage and analyzed changes in the spike (S), membrane (M), and nucleocapsid (N) gene sequences and pathogenicity of the virus at the 34th, 61st, and 100th passage levels. Sequence analyses revealed a strong selection for the S gene of 83P-5 in Vero cells, and virtually all mutations occurring at the 34th and 61st passages had been carried over to the 100th-passaged virus. In contrast, the viral M and N genes showed a strong conservation during the serial passage. Pigs experimentally infected with the 34th- or 61st-passaged virus, but not the 100th-passaged virus, exhibited diarrhea, indicating an attenuation of the 83P-5 at the 100th passage. Interestingly, S protein of the attenuated 100th-passaged 83P-5 showed a remarkable sequence similarity to that of previously reported DR-13 strain of attenuated PEDV that also had been established by serial passage in Vero cells. Further studies will be required to define whether the mutations in the S gene of 83P-5 that had been selected and accumulated during the serial passages are indeed the causalities of the growth adaptation in vitro and the attenuation of virulence in vivo.

Diffuse pleural thickening and thoracic contraction: An indistinguishable case from malignant pleural mesothelioma.

The differential diagnosis of reactive mesothelial hyperplasia and mesothelioma is difficult. We present a rare case of diffuse pleural thickening with thoracic contraction that was indistinguishable from mesothelioma. A 66-year-old woman with no history of asbestos exposure visited our hospital with a complaint of dyspnea. The clinical findings included circumferential pleural thickening on chest computed tomography image and a high concentration of hyaluronic acid in the pleural fluid. Pleural biopsies obtained by thoracoscopy under local anesthesia were pathologically consistent with mesothelioma, but the patient refused to take any kind of mesothelioma treatments. Four months later, she consented to a surgical pleural biopsy under general anesthesia to obtain larger tissue samples, which included typical proliferating polygonal cells positive for CAM5.2, calretinin, WT-1, D2-40, CK5/6, epithelial membrane antigen, and glucose transporter-1 and negative for carcinoembryonic antigen, BerEP4, and MOC31. The analysis was consistent with diagnosis of epithelioid mesothelioma. Fluorescence in situ hybridization, however, showed the presence of p16 gene, and the expression of BRCA1-associated protein-1 was detected by immunohistochemistry. Our final diagnosis was diffuse pleural thickening unrelated to asbestos exposure. Differential diagnosis of diffuse pleural thickening and malignant mesothelioma is thus difficult and routine immunohistochemical examinations are often insufficient for accurate diagnosis. Multiple diagnostic methods are required for correct diagnosis in a clinically marginal case.

Molecular characterization and expression of the equine M(1) and M(2)-pyruvate kinase gene.

To elucidate the molecular properties of the equine glycolytic enzymes equine M(1) (eM(1)) and M(2) (eM(2))-pyruvate kinase (PK), mRNAs were isolated from thoroughbred horse skeletal muscle and hair roots, respectively. The full-length eM(1) and eM(2)-PK cDNAs consist of 2,320 and 2,376 bp, respectively, containing a 1596 bp open reading frame. The cDNAs were mapped to equine chromosome 1, and the equine pyruvate kinase M (PKM) gene consists of twelve exons. Exon 9 of eM(1)-PK and exon 10 of eM(2)-PK were further investigated in five equine species. Out of 55 amino acids encoded by exon 9 in equines, Glu 418 and Lys 422, which are conserved in all PK isozymes among other vertebrates, were substituted by Gln 418 and His 422. Also, the transcriptional regulatory element(s), which have potential for involvement in alternative splicing between these exons, were completely conserved among the equines. In semi-quantitative RT-PCR analysis, strong expression of both eM(1) and eM(2)-PK mRNAs was found in skeletal muscle, heart, and brain of thoroughbred horses. In addition, the authors made the novel finding that eM(2)-PK derived from hair roots has a transcriptional start site different from that of other tissues and is more specific in its expression. These results suggest that eM(1) and eM(2)-PKs may have kinetic properties and transcriptional regulatory mechanisms different from those of other mammals.

Nanomanipulation and nanotechnology for future diagnostics.

Nanomanipulation is a technology to manipulate a small object sized in nanometer to submicron scale. Optical tweezers is one of nanomanipulation techniques, which can investigate pico-newton to femto-newton force exerted on microscopic objects. We have developed a cell palpation system by use of optical tweezers and performed palpation experiments on cells. With the cell palpation system, an operator manipulates a probe particle to touch a certain location of a cell and feels the strength of the cell by hand through a haptic device, which displays force calculated and generated by a computer. We expect this technique can be used in diagnostic purpose and utilized not only in research field but also in daily medicine.