In-House Validation of Multiplex PCR for Simultaneous Detection of Shiga Toxin-Producing Escherichia coli, Listeria monocytogenes and Salmonella spp. in Raw Meats.

The aim of the study was to perform in-house validation of the developed multiplex PCR (mPCR)-based alternative method to detect Shiga toxin-producing Escherichia coli (STEC), Listeria monocytogenes (L. monocytogenes) and Salmonella spp. in raw meats following the ISO 16140-2: 2016. A comparative study of the developed mPCR against the Bacteriological Analytical Manual (BAM) method was evaluated for inclusivity and exclusivity, sensitivity and the relative level of detection (RLOD). Inclusivity levels for each target bacterium were all 100%, while exclusivity for non-target bacteria was 100%. The sensitivity of the developed mPCR was calculated based on the analysis of 72 samples of raw meat. The sensitivity of the developed mPCR was 100%. The RLOD values of the developed mPCR for STEC, L. monocytogenes and Salmonella spp. were 0.756, 1.170 and 1.000, respectively. The developed mPCR showed potential as a tool for the fast, specific and sensitive detection of the three bacteria in the raw meat industry.

An Implementation Framework for Telemedicine to Address Noncommunicable Diseases in Thailand.

To maintain the continuity of noncommunicable disease (NCD) services and ascertain the health outcomes of patients with NCDs during the COVID-19 (coronavirus disease 2019) outbreak in Thailand, various telemedicine services have been developed. To achieve this determination, the implementation framework has been constructed based on recommendations from multidisciplinary experts (Thai NCD Collaboration Group). Within the framework, all key elements are illustrated with their priority and expected collaborations. Ultimately, active collaborations from multi-stakeholders are vitally important to ensure that telemedicine services for NCDs will finally become practical, successful, and sustainable.

Bacterial communities on facial skin of teenage and elderly Thai females.

The Human Microbiome Project was first established to understand the roles of human-associated microbes to human health and disease. This study presents preliminary findings of Thai female facial skin microbiome using three pooled samples from groups of skin microbiome profiles, namely (1) healthy and (2) acne-prone young adults (teenage.hea and teenage.acn) and (3) healthy elderly adults (elderly.hea) based on standard dermatological criteria. These samples were sequenced using 454-pyrosequencing targeting 16S rRNA (V3-V4 regions). Good's coverage index of greater than 92% shows sufficient sampling of our data for each group. Three unique OTUs for each microbiome profile (43, 258 and 59 for teenage.hea, teenage.acn and ederly.hea, respectively) were obtained with 134 shared OTUs among the three datasets. Based on Morisita-Horn similarity coefficient, age is the major factor that brings the community relationship factor closer. The comparison among the three datasets reveal majority of Gemmatimonadetes, Planctomycetes and Nitrospirae in the teenage.hea, whereas Firmicutes are more prevalent in teenage.acn and elderly.hea skin types. In addition, when comparing Thai facial microbial diversity with the 16S data from U.S. forehead female database, significant differences were found among orders of bacteria, pointing to possible differences in human ecto-flora.

Genetic analysis of Thai cattle reveals a Southeast Asian indicine ancestry.

Cattle commonly raised in Thailand have characteristics of Bos indicus (zebu). We do not know when or how cattle domestication in Thailand occurred, and so questions remain regarding their origins and relationships to other breeds. We obtained genome-wide SNP genotypic data of 28 bovine individuals sampled from four regions: North (Kho-Khaolampoon), Northeast (Kho-Isaan), Central (Kho-Lan) and South (Kho-Chon) Thailand. These regional varieties have distinctive traits suggestive of breed-like genetic variations. From these data, we confirmed that all four Thai varieties are Bos indicus and that they are distinct from other indicine breeds. Among these Thai cattle, a distinctive ancestry pattern is apparent, which is the purest within Kho-Chon individuals. This ancestral component is only present outside of Thailand among other indicine breeds in Southeast Asia. From this pattern, we conclude that a unique Bos indicus ancestor originated in Southeast Asia, and native Kho-Chon Thai cattle retain the signal of this ancestry with limited admixture of other bovine ancestors.

HLA-DRB1 and HLA-DQB1 Are Associated with Adult-Onset Immunodeficiency with Acquired Anti-Interferon-Gamma Autoantibodies.

Recently a newly identified clinical syndrome of disseminated non-tuberculous mycobacterial diseases (with or without other opportunistic infections in adult patients who were previously healthy, has been recognized in association with an acquired autoantibody to interferon-gamma. This syndrome is emerging as an important cause of morbidity and mortality, especially among people of Asian descent. Trigger for the production of this autoantibody remains unknown, but genetic factors are strongly suspected to be involved. We compared HLA genotyping between 32 patients with this clinical syndrome, and 38 controls. We found that this clinical syndrome was associated with very limited allele polymorphism, with HLA-DRB1 and DQB1 alleles, especially HLA-DRB1*15:01, DRB1*16:02, DQB1*05:01 and DQB1*05:02. Odds ratio of DRB1*15:01, DRB1*16:02, DQB1*05:01 and DQB1*05:02 were 7.03 (95% CI, 2.18-22.69, P<0.0001, 9.06 (95% CI, 2.79-29.46, P<0.0001), 6.68 (95% CI, 2.29-19.52, P = 0.0004), and 6.64 (95% CI, 2.30-19.20, P = 0.0004), respectively. Further investigation is warranted to provide better understanding on pathogenesis of this association.

Mycophenolic acid AUC in Thai kidney transplant recipients receiving low dose mycophenolate and its association with UGT2B7 polymorphisms.

BACKGROUND: Despite use of a lower mycophenolate dose in Thai kidney transplant patients, acceptable graft and patient outcomes can be achieved. We therefore examined the pharmacokinetics of mycophenolic acid (MPA) by area under the curve (AUC) and investigated genetic contribution in mycophenolate metabolism in this population. METHODS: Kidney transplant recipients with stable graft function who were receiving mycophenolate mofetil 1,000 mg/d in combination with either cyclosporine or tacrolimus, and prednisolone were studied. The MPA concentration was measured by fluorescence polarization immunoassay (FPIA), at predose and 1, 1.5, 2, 4, 6, 8, 10, and 12 hours after dosing. Genetic polymorphisms in UGT1A8, UGT1A9, and UGT2B7 were examined by denaturing high-performance liquid chromatography (DHPLC)-based single-base extension (SBE) analysis. RESULTS: A total 138 patients were included in study. The mean AUC was 39.49 mg-h/L (28.39-89.58 mg-h/L), which was in the therapeutic range. The correlation between the predose MPA concentration and AUC was poor. The mean AUC in the tacrolimus group was higher than that in the cyclosporine group. Polymorphisms in UGT2B7 showed significant association with AUC. CONCLUSION: Most of our patients with reduced mycophenolate dose had the AUC within the therapeutic range. Genetic polymorphisms in UGT2B7 may play a role in MPA metabolism in Thai kidney transplant patients.

Structural and functional diversity of free-living microorganisms in reef surface, Kra island, Thailand.

BACKGROUND: Coral reefs worldwide are being harmed through anthropogenic activities. Some coral reefs in Thailand remain well-preserved, including the shallow coral reefs along Kra island, Nakhon Si Thammarat province. Interestingly, the microbial community in this environment remains unknown. The present study identified biodiversity of prokaryotes and eukaryotes of 0.22-30 µm in sizes and their metabolic potentials in this coral reef surface in summer and winter seasons, using 16S and 18S rRNA genes pyrosequencing. RESULTS: The marine microbial profiles in summer and winter seasons comprised mainly of bacteria, in phylum, particular the Proteobacteria. Yet, different bacterial and eukaryotic structures existed between summer and winter seasons, supported by low Lennon and Yue & Clayton theta similarity indices (8.48-10.43% for 16S rRNA, 0.32-7.81% for 18S rRNA ). The topmost prokaryotic phylum for the summer was Proteobacteria (99.68%), while for the winter Proteobacteria (62.49%) and Bacteroidetes (35.88%) were the most prevalent. Uncultured bacteria in phyla Cyanobacteria, Planctomycetes, SAR406 and SBR1093 were absent in the summer. For eukaryotic profiles, species belonging to animals predominated in the summer, correlating with high animal activities in the summer, whereas dormancy and sporulation predominated in the winter. For the winter, eukaryotic plant species predominated and several diverse species were detected. Moreover, comparison of our prokaryotic databases in summer and winter of Kra reef surface against worldwide marine culture-independent prokaryotic databases indicated our databases to most resemblance those of coastal Sichang island, Chonburi province, Thailand, and the 3 tropical GOS sites close to Galapagos island (GS039, GS040 and GS045), in orderly. CONCLUSIONS: The study investigated and obtained culture-independent databases for marine prokaryotes and eukaryotes in summer and winter seasons of Kra reef surface. The data helped understand seasonal dynamics of microbial structures and metabolic potentials of this tropical ecosystem, supporting the knowledge of the world marine microbial biodiversity.

Microbial ecology of Thailand tsunami and non-tsunami affected terrestrials.

The effects of tsunamis on microbial ecologies have been ill-defined, especially in Phang Nga province, Thailand. This ecosystem was catastrophically impacted by the 2004 Indian Ocean tsunami as well as the 600 year-old tsunami in Phra Thong island, Phang Nga province. No study has been conducted to elucidate their effects on microbial ecology. This study represents the first to elucidate their effects on microbial ecology. We utilized metagenomics with 16S and 18S rDNA-barcoded pyrosequencing to obtain prokaryotic and eukaryotic profiles for this terrestrial site, tsunami affected (S1), as well as a parallel unaffected terrestrial site, non-tsunami affected (S2). S1 demonstrated unique microbial community patterns than S2. The dendrogram constructed using the prokaryotic profiles supported the unique S1 microbial communities. S1 contained more proportions of archaea and bacteria domains, specifically species belonging to Bacteroidetes became more frequent, in replacing of the other typical floras like Proteobacteria, Acidobacteria and Basidiomycota. Pathogenic microbes, including Acinetobacter haemolyticus, Flavobacterium spp. and Photobacterium spp., were also found frequently in S1. Furthermore, different metabolic potentials highlighted this microbial community change could impact the functional ecology of the site. Moreover, the habitat prediction based on percent of species indicators for marine, brackish, freshwater and terrestrial niches pointed the S1 to largely comprise marine habitat indicating-species.

Population structure of four Thai indigenous chicken breeds.

BACKGROUND: In recent years, Thai indigenous chickens have increasingly been bred as an alternative in Thailand poultry market. Due to their popularity, there is a clear need to improve the underlying quality and productivity of these chickens. Studying chicken genetic variation can improve the chicken meat quality as well as conserving rare chicken species. To begin with, a minimal set of molecular markers that can characterize the Thai indigenous chicken breeds is required. RESULTS: Using AFLP-PCR, 30 single nucleotide polymorphisms (SNPs) from Thai indigenous chickens were obtained by DNA sequencing. From these SNPs, we genotyped 465 chickens from 7 chicken breeds, comprising four Thai indigenous chicken breeds--Pradhuhangdum (PD), Luenghangkhao (LK), Dang (DA) and Chee (CH), one wild chicken--the red jungle fowls (RJF), and two commercial chicken breeds--the brown egg layer (BL) and commercial broiler (CB). The chicken genotypes reveal unique genetic structures of the four Thai indigenous chicken breeds. The average expected heterozygosities of PD=0.341, LK=0.357, DA=0.349 and CH=0.373, while the references RJF= 0.327, CB=0.324 and BL= 0.285. The F(ST) values among Thai indigenous chicken breeds vary from 0.051 to 0.096. The F(ST) values between the pairs of Thai indigenous chickens and RJF vary from 0.083 to 0.105 and the FST values between the Thai indigenous chickens and the two commercial chicken breeds vary from 0.116 to 0.221. A neighbour-joining tree of all individual chickens showed that the Thai indigenous chickens were clustered into four groups which were closely related to the wild RJF but far from the commercial breeds. Such commercial breeds were split into two closely groups. Using genetic admixture analysis, we observed that the Thai indigenous chicken breeds are likely to share common ancestors with the RJF, while both commercial chicken breeds share the same admixture pattern. CONCLUSION: These results indicated that the Thai indigenous chicken breeds may descend from the same ancestors. These indigenous chicken breeds were more closely related to red jungle fowls than those of the commercial breeds. These findings showed that the proposed SNP panel can effectively be used to characterize the four Thai indigenous chickens.

How bio-questionable are the different recombinant human erythropoietin copy products in Thailand?

PURPOSE: The high prevalence of pure red cell aplasia in Thailand has been associated with the sharp increase in number of recombinant human erythropoietin (rhEPO) copy products, based on a classical generic regulatory pathway, which have entered the market. This study aims to assess the quality of rhEPO copy products being used in Thailand. METHODS: Twelve rhEPO copy products were purchased from pharmacies in Thailand, shipped under controlled cold chain conditions to the Netherlands and characterized using (1) high performance size-exclusion chromatography, (2) asymmetrical flow field-flow fractionation, (3) sodium dodecyl sulfate polyacrylamide gel electrophoresis in combination with (4) Western blotting and additionally tested for (5) host cell protein impurities as well as (6) endotoxin contamination. RESULTS: Some of the tested rhEPO copy products showed high aggregate levels and contained a substantial amount of protein fragments. Also, one of rhEPO copy products had a high endotoxin level, exceeding the FDA limit. CONCLUSIONS: Our observations show that some of the tested copy products on the Thai market differ significantly from the originator rhEPO product, Epogen®. This comparison study supports a link between the quality attributes of copy rhEPO products and their immunogenicity.

iNJclust: Iterative Neighbor-Joining Tree Clustering Framework for Inferring Population Structure.

Understanding genetic differences among populations is one of the most important issues in population genetics. Genetic variations, e.g., single nucleotide polymorphisms, are used to characterize commonality and difference of individuals from various populations. This paper presents an efficient graph-based clustering framework which operates iteratively on the Neighbor-Joining (NJ) tree called the iNJclust algorithm. The framework uses well-known genetic measurements, namely the allele-sharing distance, the neighbor-joining tree, and the fixation index. The behavior of the fixation index is utilized in the algorithm's stopping criterion. The algorithm provides an estimated number of populations, individual assignments, and relationships between populations as outputs. The clustering result is reported in the form of a binary tree, whose terminal nodes represent the final inferred populations and the tree structure preserves the genetic relationships among them. The clustering performance and the robustness of the proposed algorithm are tested extensively using simulated and real data sets from bovine, sheep, and human populations. The result indicates that the number of populations within each data set is reasonably estimated, the individual assignment is robust, and the structure of the inferred population tree corresponds to the intrinsic relationships among populations within the data.

Insight into the peopling of Mainland Southeast Asia from Thai population genetic structure.

There is considerable ethno-linguistic and genetic variation among human populations in Asia, although tracing the origins of this diversity is complicated by migration events. Thailand is at the center of Mainland Southeast Asia (MSEA), a region within Asia that has not been extensively studied. Genetic substructure may exist in the Thai population, since waves of migration from southern China throughout its recent history may have contributed to substantial gene flow. Autosomal SNP data were collated for 438,503 markers from 992 Thai individuals. Using the available self-reported regional origin, four Thai subpopulations genetically distinct from each other and from other Asian populations were resolved by Neighbor-Joining analysis using a 41,569 marker subset. Using an independent Principal Components-based unsupervised clustering approach, four major MSEA subpopulations were resolved in which regional bias was apparent. A major ancestry component was common to these MSEA subpopulations and distinguishes them from other Asian subpopulations. On the other hand, these MSEA subpopulations were admixed with other ancestries, in particular one shared with Chinese. Subpopulation clustering using only Thai individuals and the complete marker set resolved four subpopulations, which are distributed differently across Thailand. A Sino-Thai subpopulation was concentrated in the Central region of Thailand, although this constituted a minority in an otherwise diverse region. Among the most highly differentiated markers which distinguish the Thai subpopulations, several map to regions known to affect phenotypic traits such as skin pigmentation and susceptibility to common diseases. The subpopulation patterns elucidated have important implications for evolutionary and medical genetics. The subpopulation structure within Thailand may reflect the contributions of different migrants throughout the history of MSEA. The information will also be important for genetic association studies to account for population-structure confounding effects.

Biomarker selection and classification of "-omics" data using a two-step bayes classification framework.

Identification of suitable biomarkers for accurate prediction of phenotypic outcomes is a goal for personalized medicine. However, current machine learning approaches are either too complex or perform poorly. Here, a novel two-step machine-learning framework is presented to address this need. First, a Naïve Bayes estimator is used to rank features from which the top-ranked will most likely contain the most informative features for prediction of the underlying biological classes. The top-ranked features are then used in a Hidden Naïve Bayes classifier to construct a classification prediction model from these filtered attributes. In order to obtain the minimum set of the most informative biomarkers, the bottom-ranked features are successively removed from the Naïve Bayes-filtered feature list one at a time, and the classification accuracy of the Hidden Naïve Bayes classifier is checked for each pruned feature set. The performance of the proposed two-step Bayes classification framework was tested on different types of -omics datasets including gene expression microarray, single nucleotide polymorphism microarray (SNParray), and surface-enhanced laser desorption/ionization time-of-flight (SELDI-TOF) proteomic data. The proposed two-step Bayes classification framework was equal to and, in some cases, outperformed other classification methods in terms of prediction accuracy, minimum number of classification markers, and computational time.

Distal effect of amino acid substitutions in CYP2C9 polymorphic variants causes differences in interatomic interactions against (S)-warfarin.

Cytochrome P450 2C9 (CYP2C9) is crucial in excretion of commonly prescribed drugs. However, changes in metabolic activity caused by CYP2C9 polymorphisms inevitably result in adverse drug effects. CYP2C9*2 and *3 are prevalent in Caucasian populations whereas CYP2C9*13 is remarkable in Asian populations. Single amino acid substitutions caused by these mutations are located outside catalytic cavity but affect kinetic activities of mutants compared to wild-type enzyme. To relate distal effects of these mutations and defective drug metabolisms, simulations of CYP2C9 binding to anti-coagulant (S)-warfarin were performed as a system model. Representative (S)-warfarin-bound forms of wild-type and mutants were sorted and assessed through knowledge-based scoring function. Interatomic interactions towards (S)-warfarin were predicted to be less favorable in mutant structures in correlation with larger distance between hydroxylation site of (S)-warfarin and reactive oxyferryl heme than wild-type structure. Using computational approach could delineate complication of CYP polymorphism in management of drug therapy.

Molecular dynamics simulations reveal structural instability of human trypsin inhibitor upon D50E and Y54H mutations.

Serine protease inhibitor Kazal type 1 (SPINK1) plays an important role in protecting the pancreas against premature trypsinogen activation that causes pancreatitis. Various mutations in the SPINK1 gene were shown to be associated with patients with pancreatitis. Recent transfection studies identified intracellular folding defects, probably caused by mutation induced misfolding of D50E and Y54H mutations, as a common mechanism that reduces SPINK1 secretion and as a possible novel mechanism of SPINK1 deficiency associated with chronic pancreatitis. Using molecular dynamics, we investigated the effects of D50E and Y54H mutations on SPINK1 dynamics and conformation at 300 K. We found that the structures of D50E and Y54H mutants were less stable than and were distorted from those of the wild type, as indicated by the RMSD plots, RMSF plots and DSSP series. Specifically, unwinding of the top of helices (the main secondary structures) and the distortion of the loops above the helices were observed. It may be possible that this distorted protein structure may be recognized as "non-native" by members of the chaperone family; it may be further retained and targeted for degradation, leading to SPINK1 secretion reduction and subsequently pancreatitis in patients as Király et al. (Gut 56:1433, 2007) proposed.

MetaSel: a metaphase selection tool using a Gaussian-based classification technique.

BACKGROUND: Identification of good metaphase spreads is an important step in chromosome analysis for identifying individuals with genetic disorders. The process of finding suitable metaphase chromosomes for accurate clinical analysis is, however, very time consuming since they are selected manually. The selection of suitable metaphase chromosome spreads thus represents a major bottleneck for conventional cytogenetic analysis. Although many algorithms have been developed for karyotyping, none have adequately addressed the critical bottleneck of selecting suitable chromosome spreads. In this paper, we present a software tool that uses a simple rule-based system to efficiently identify metaphase spreads suitable for karyotyping. RESULTS: The chromosome shapes can be classified by the software into four main classes. The first and the second classes refer to individual chromosomes with straight and skewed shapes, respectively. The third class is characterized as those chromosomes with overlapping bodies and the fourth class is for the non-chromosome objects. Good metaphase spreads should largely contain chromosomes of the first and the second classes, while the third class should be kept minimal. Several image parameters were examined and used for creating rule-based classification. The threshold value for each parameter is determined using a statistical model. We observed that the Gaussian model can represent the empirical probability density function of the parameters and, hence, the threshold value can be easily determined. The proposed rules can efficiently and accurately classify the individual chromosome with > 90% accuracy. CONCLUSIONS: The software tool, termed MetaSel, was developed. Using the Gaussian-based rules, the tool can be used to quickly rank hundreds of chromosome spread images so as to assist cytogeneticists to perform karyotyping effectively. Furthermore, MetaSel offers an intuitive, yet comprehensive, workflow to assist karyotyping, including tools for editing chromosome (split, merge and fix) and a karyotyping editor (moving, rotating, and pairing homologous chromosomes). The program can be freely downloaded from "http://www4a.biotec.or.th/GI/tools/metasel".

A novel maternally-derived insertional translocation resulting in partial trisomy 4q13.2-q22.1 with complex translocation t(8;20) in a family with intellectual disability.

We characterized the chromosomal aberration in family with intellectual disability, including two affected children and their affected mother. Initial standard karyotypes of the three individuals showed an apparently balanced translocation of chromosomes 8 and 20. Using molecular cytogenetic techniques, we observed complex structural chromosomal aberration comprising of reciprocal translocation between chromosomes 8 and 20 with pericentric inversion (8p11.12q22.3) and insertion of chromosome 4 segments into both der(8) and der(20). In particular, the insertion of chromosome 4 was complex. Two segments (4q13.2-q13.3 and 4q21.21-q22.1) were inserted into the der(8)t(8;20) breakpoint and one segment (4q13.3-q21.21) into the der(20)t(8;20) breakpoint. Both children inherited two normal chromosomes 4 from their parents and the der(8) and der(20) from the mother, resulting in partial trisomy of 4q13.2-q22.1. Interestingly, the mother, in addition to the same complex insertions and inversion, was founded to have a deletion of 4q13.2-q22.1 in one of her chromosomes 4, yielding no genetic imbalance but with potential disruption of intellectual dysfunction-related gene(s) at the breakpoints as the cause of her intellectual impairment. This family is the third case report of an insertional translocation mechanism causing partial trisomy 4q syndrome. Our study demonstrates that an insertion of an extra chromosomal segment, not primarily involving in translocation breakpoints, which results in partial trisomy, can be an unapparent cause of the abnormal phenotypes.

Small Ancestry Informative Marker panels for complete classification between the original four HapMap populations.

A protocol for the identification of Ancestry Informative Markers (AIMs) from genome-wide Single Nucleotide Polymorphism (SNP) data is proposed. The protocol consists of three main steps: identification of potential positive selection regions via F(ST) extremity measurement, SNP screening via two-stage attribute selection and classification model construction using a Naïve Bayes classifier. The two-stage attribute selection is composed of a newly developed round robin Symmetrical Uncertainty (SU) ranking technique and a wrapper embedded with a Naïve Bayes classifier. The protocol has been applied to the HapMap Phase II data. Two AIM panels, which consist of 10 and 16 SNPs that lead to complete classification between CEU, CHB, JPT and YRI populations, are identified. Moreover, the panels are at least four times smaller than those reported in previous studies. The results suggest that the protocol could be useful in a scenario involving a larger number of populations.

iLOCi: a SNP interaction prioritization technique for detecting epistasis in genome-wide association studies.

BACKGROUND: Genome-wide association studies (GWAS) do not provide a full account of the heritability of genetic diseases since gene-gene interactions, also known as epistasis are not considered in single locus GWAS. To address this problem, a considerable number of methods have been developed for identifying disease-associated gene-gene interactions. However, these methods typically fail to identify interacting markers explaining more of the disease heritability over single locus GWAS, since many of the interactions significant for disease are obscured by uninformative marker interactions e.g., linkage disequilibrium (LD). RESULTS: In this study, we present a novel SNP interaction prioritization algorithm, named iLOCi (Interacting Loci). This algorithm accounts for marker dependencies separately in case and control groups. Disease-associated interactions are then prioritized according to a novel ranking score calculated from the difference in marker dependencies for every possible pair between case and control groups. The analysis of a typical GWAS dataset can be completed in less than a day on a standard workstation with parallel processing capability. The proposed framework was validated using simulated data and applied to real GWAS datasets using the Wellcome Trust Case Control Consortium (WTCCC) data. The results from simulated data showed the ability of iLOCi to identify various types of gene-gene interactions, especially for high-order interaction. From the WTCCC data, we found that among the top ranked interacting SNP pairs, several mapped to genes previously known to be associated with disease, and interestingly, other previously unreported genes with biologically related roles. CONCLUSION: iLOCi is a powerful tool for uncovering true disease interacting markers and thus can provide a more complete understanding of the genetic basis underlying complex disease. The program is available for download at http://www4a.biotec.or.th/GI/tools/iloci.