Indranil Mukhopadhyay

Evaluation and comparison of gene clustering methods in microarray analysis

MOTIVATION Microarray technology has been widely applied in biological and clinical studies for simultaneous monitoring of gene expression in thousands of genes. Gene clustering analysis is found useful for discovering groups of correlated genes potentially co-regulated or associated to the disease or conditions under investigation. Many clustering methods including hierarchical clustering, K-means, PAM, SOM, mixture model-based clustering and tight clustering have been widely used in the literature. Yet no comprehensive comparative study has been performed to evaluate the effectiveness of these methods. RESULTS In this paper, six gene clustering methods are evaluated by simulated data from a hierarchical log-normal model with various degrees of perturbation as well as four real datasets. A weighted Rand index is proposed for measuring similarity of two clustering results with possible scattered genes (i.e. a set of noise genes not being clustered). Performance of the methods in the real data is assessed by a predictive accuracy analysis through verified gene annotations. Our results show that tight clustering and model-based clustering consistently outperform other clustering methods both in simulated and real data while hierarchical clustering and SOM perform among the worst. Our analysis provides deep insight to the complicated gene clustering problem of expression profile and serves as a practical guideline for routine microarray cluster analysis.

Association tests using kernel-based measures of multi-locus genotype similarity between individuals

In a genetic association study, it is often desirable to perform an overall test of whether any or all single‐nucleotide polymorphisms (SNPs) in a gene are associated with a phenotype. Several such tests exist, but most of them are powerful only under very specific assumptions about the genetic effects of the individual SNPs. In addition, some of the existing tests assume that the direction of the effect of each SNP is known, which is a highly unlikely scenario. Here, we propose a new kernel‐based association test of joint association of several SNPs. Our test is non‐parametric and robust, and does not make any assumption about the directions of individual SNP effects. It can be used to test multiple correlated SNPs within a gene and can also be used to test independent SNPs or genes in a biological pathway. Our test uses an analysis of variance paradigm to compare variation between cases and controls to the variation within the groups. The variation is measured using kernel functions for each marker, and then a composite statistic is constructed to combine the markers into a single test. We present simulation results comparing our statistic to the U‐statistic‐based method by Schaid et al. ([2005] Am. J. Hum. Genet. 76:780–793) and another statistic by Wessel and Schork ([2006] Am. J. Hum. Genet. 79:792–806). We consider a variety of different disease models and assumptions about how many SNPs within the gene are actually associated with disease. Our results indicate that our statistic has higher power than other statistics under most realistic conditions. Genet. Epidemiol. 34: 213–221, 2010.

Dopaminergic Mutations : Within-Family Association and Linkage in Multiplex Alcohol Dependence Families

Animal and human studies of addiction indicate that the D2 dopamine receptor (DRD2) plays a critical role in the mechanism of drug reward. D2 receptor density in the brains of alcoholics has been shown to be reduced relative to controls. Previous studies of DRD2 in association with alcohol dependence using variation in the TaqI A locus were highly controversial. Recently, a synonymous mutation, C957T, in the coding region of the human DRD2 gene has been identified which appears to have functional effects including alteration in receptor availability. In order to determine if susceptibility to alcohol dependence (AD) within multiplex alcohol dependence families would be altered by the C957T in the coding region of the D2 gene, within‐family association was studied in members of Caucasian multiplex alcohol dependence families. Members of control families with no personal alcohol or substance dependence history were included for case/control comparisons. Analyses performed to detect within‐family association showed evidence favoring an association for the C957T polymorphism (P = 0.038). Linkage analyses of polymorphisms in this region showed that only the C957T locus remained of interest (P = 0.015). Evidence for the C957T T allele having a role in AD susceptibility at the population level using a case/control comparison was statistically marginal (P = 0.062), but was consistent with the family data results. These results support a role for DRD2 as a susceptibility gene for alcohol dependence within multiplex families at high risk for developing alcohol dependence.

Transient receptor potential ankyrin 1 receptor activation in vitro and in vivo by pro-tussive agents: GRC 17536 as a promising anti-tussive therapeutic.

Cough is a protective reflex action that helps clear the respiratory tract which is continuously exposed to airborne environmental irritants. However, chronic cough presents itself as a disease in its own right and despite its global occurrence; the molecular mechanisms responsible for cough are not completely understood. Transient receptor potential ankyrin1 (TRPA1) is robustly expressed in the neuronal as well as non-neuronal cells of the respiratory tract and is a sensor of a wide range of environmental irritants. It is fast getting acceptance as a key biological sensor of a variety of pro-tussive agents often implicated in miscellaneous chronic cough conditions. In the present study, we demonstrate in vitro direct functional activation of TRPA1 receptor by citric acid which is routinely used to evoke cough in preclinical and clinical studies. We also show for the first time that a potent and selective TRPA1 antagonist GRC 17536 inhibits citric acid induced cellular Ca+2 influx in TRPA1 expressing cells and the citric acid induced cough response in guinea pigs. Hence our data provides a mechanistic link between TRPA1 receptor activation in vitro and cough response induced in vivo by citric acid. Furthermore, we also show evidence for TRPA1 activation in vitro by the TLR4, TLR7 and TLR8 ligands which are implicated in bacterial/respiratory virus pathogenesis often resulting in chronic cough. In conclusion, this study highlights the potential utility of TRPA1 antagonist such as GRC 17536 in the treatment of miscellaneous chronic cough conditions arising due to diverse causes but commonly driven via TRPA1.

Association of viral load with HPV16 positive cervical cancer pathogenesis: causal relevance in isolates harboring intact viral E2 gene.

We tested the hypothesis that cervical cancers (CaCx) harbor high HPV16 viral load compared to controls and this is influenced by E2 status and age of subjects. Viral load (natural log transformed values) per 100ng genomic DNA was estimated (152 cases and 87 controls) by Taqman assay. Median viral load was significantly higher (Mann-Whitney U test) among cases (17.21) compared to controls (9.86), irrespective of E2 status or upon considering E2 status as a covariate in logistic regression model (p<0.001). Viral load of E2 intact cases (17.80) was significantly higher (p<0.001) compared to those with disrupted E2 (9.78). At equivalent probability of being a case, viral load was higher among individuals (i) of lower age, irrespective of E2 status, and (ii) with intact E2 but of similar age as those with disrupted E2. Thus viral load in association with E2 status and/or age might be of causal relevance in CaCx pathogenesis.

A quest for miRNA bio-marker: a track back approach from gingivo buccal cancer to two different types of precancers.

Deregulation of miRNA expression may contribute to tumorigenesis and other patho-physiology associated with cancer. Using TLDA, expression of 762 miRNAs was checked in 18 pairs of gingivo buccal cancer-adjacent control tissues. Expression of significantly deregulated miRNAs was further validated in cancer and examined in two types of precancer (leukoplakia and lichen planus) tissues by primer-specific TaqMan assays. Biological implications of these miRNAs were assessed bioinformatically. Expression of hsa-miR-1293, hsa-miR-31, hsa-miR-31* and hsa-miR-7 were significantly up-regulated and those of hsa-miR-206, hsa-miR-204 and hsa-miR-133a were significantly down-regulated in all cancer samples. Expression of only hsa-miR-31 was significantly up-regulated in leukoplakia but none in lichen planus samples. Analysis of expression heterogeneity divided 18 cancer samples into clusters of 13 and 5 samples and revealed that expression of 30 miRNAs (including the above-mentioned 7 miRNAs), was significantly deregulated in the cluster of 13 samples. From database mining and pathway analysis it was observed that these miRNAs can significantly target many of the genes present in different cancer related pathways such as “proteoglycans in cancer”, PI3K-AKT etc. which play important roles in expression of different molecular features of cancer. Expression of hsa-miR-31 was significantly up-regulated in both cancer and leukoplakia tissues and, thus, may be one of the molecular markers of leukoplakia which may progress to gingivo-buccal cancer.

Some Novel Insights on HPV16 Related Cervical Cancer Pathogenesis Based on Analyses of LCR Methylation, Viral Load, E7 and E2/E4 Expressions

This study was undertaken to decipher the interdependent roles of (i) methylation within E2 binding site I and II (E2BS-I/II) and replication origin (nt 7862) in the long control region (LCR), (ii) expression of viral oncogene E7, (iii) expression of the transcript (E7-E1∧E4) that encodes E2 repressor protein and (iv) viral load, in human papillomavirus 16 (HPV16) related cervical cancer (CaCx) pathogenesis. The results revealed over-representation (p<0.001) of methylation at nucleotide 58 of E2BS-I among E2-intact CaCx cases compared to E2-disrupted cases. Bisulphite sequencing of LCR revealed overrepresentation of methylation at nucleotide 58 or other CpGs in E2BS-I/II, among E2-intact cases than E2-disrupted cases and lack of methylation at replication origin in case of both. The viral transcript (E7-E1∧E4) that produces the repressor E2 was analyzed by APOT (amplification of papillomavirus oncogenic transcript)-coupled-quantitative-RT-PCR (of E7 and E4 genes) to distinguish episomal (pure or concomitant with integrated) from purely integrated viral genomes based on the ratio, E7 CT/E4 CT. Relative quantification based on comparative CT (theshold cycle) method revealed 75.087 folds higher E7 mRNA expression in episomal cases over purely integrated cases. Viral load and E2 gene copy numbers were negatively correlated with E7 CT (p = 0.007) and E2 CT (p<0.0001), respectively, each normalized with ACTB CT, among episomal cases only. The k-means clustering analysis considering E7 CT from APOT-coupled-quantitative-RT-PCR assay, in conjunction with viral load, revealed immense heterogeneity among the HPV16 positive CaCx cases portraying integrated viral genomes. The findings provide novel insights into HPV16 related CaCx pathogenesis and highlight that CaCx cases that harbour episomal HPV16 genomes with intact E2 are likely to be distinct biologically, from the purely integrated viral genomes in terms of host genes and/or pathways involved in cervical carcinogenesis.

Association of IL1A and IL1B loci with primary open angle glaucoma

BackgroundRecent studies suggest that glaucoma is a neurodegenerative disease in which secondary degenerative losses occur after primary insult by raised Intraocular pressure (IOP) or by other associated factors. It has been reported that polymorphisms in the IL1A and IL1B genes are associated with Primary Open Angle Glaucoma (POAG). The purpose of our study was to investigate the role of these polymorphisms in eastern Indian POAG patients.MethodsThe study involved 315 unrelated POAG patients, consisting of 116 High Tension Glaucoma (HTG) patients with intra ocular pressure (IOP) > 21 mmHg and 199 non-HTG patients (presenting IOP < 20 mmHg), and 301 healthy controls from eastern India. Genotypes were determined by polymerase chain reaction and restriction digestion for three single nucleotide polymorphisms (SNPs): IL1A (-889C/T; rs1800587), IL1B (-511C/T; rs16944) and IL1B (3953C/T; rs1143634). Haplotype frequency was determined by Haploview 4.1 software. The association of individual SNPs and major haplotypes was evaluated using chi-square statistics. The p-value was corrected for multiple tests by Bonferroni method.ResultsNo significant difference was observed in the allele and genotype frequencies for IL1A and IL1B SNPs between total pool of POAG patients and controls. However, on segregating the patient pool to HTG and non-HTG groups, weak association was observed for IL1A polymorphism (-889C/T) where -889C allele was found to portray risk (OR = 1.380; 95% CI = 1.041-1.830; p = 0.025) for non-HTG patients. Similarly, 3953T allele of IL1B polymorphism (+3953C/T) was observed to confer risk to HTG group (OR = 1.561; 95% CI = 1.022-2.385; p = 0.039). On haplotype analysis it was observed that TTC was significantly underrepresented in non-HTG patients (OR = 0.538; 95% CI = 0.356- 0.815; p = 0.003) while TCT haplotype was overrepresented in HTG patients (OR = 1.784; 95% CI = 1.084- 2.937; p = 0.022) compared to control pool. However, after correction for multiple tests by Bonferroni method, an association of only TTC haplotype with non-HTG cases sustained (pcorrected = 0.015) and expected to confer protection.ConclusionThe study suggests that the genomic region containing the IL1 gene cluster influences the POAG pathogenesis mostly in non-HTG patients in eastern India. A similar study in additional and larger cohorts of patients in other population groups is necessary to further substantiate the observation.

Association between risk of oral precancer and genetic variations in microRNA and related processing genes

BackgroundMicroRNAs have been implicated in cancer but studies on their role in precancer, such as leukoplakia, are limited. Sequence variations at eight miRNA and four miRNA processing genes were studied in 452 healthy controls and 299 leukoplakia patients to estimate risk of disease.ResultsGenotyping by TaqMan assay followed by statistical analyses showed that variant genotypes at Gemin3 and mir-34b reduced risk of disease [OR = 0.5(0.3–0.9) and OR = 0.7(0.5–0.9) respectively] in overall patients as well as in smokers [OR = 0.58(0.3–1) and OR = 0.68(0.5–0.9) respectively]. Among chewers, only mir29a significantly increased risk of disease [OR = 1.8(1–3)]. Gene-environment interactions using MDR-pt program revealed that mir29a, mir34b, mir423 and Xpo5 modulated risk of disease (p < 0.002) which may be related to change in expression of these genes as observed by Real-Time PCR assays. But association between polymorphisms and gene expressions was not found in our sample set as well as in larger datasets from open access platforms like Genevar and 1000 Genome database.ConclusionVariations in microRNAs and their processing genes modulated risk of precancer but further in-depth study is needed to understand mechanism of disease process.

Increased Risk of Psoriasis due to combined effect of HLA-Cw6 and LCE3 risk alleles in Indian population

HLA-Cw6 is one of the most associated alleles in psoriasis. Recently, Late Cornified Envelop 3 (LCE3) genes were identified as a susceptibility factor for psoriasis. Some population showed epistatic interaction of LCE3 risk variants with HLA-Cw6, while some population failed to show any association. We determined the associations of a 32.2 kb deletion comprising LCE3C-3B genes and three SNPs (rs1886734, rs4112788; rs7516108) at the LCE3 gene cluster among the psoriasis patients in India. All three SNPs at the LCE3 gene cluster failed to show any association. In contrary, for patients with HLA-Cw6 allele, all three SNPs and the LCE3C-3B deletion showed significant associations. While, all five LCE3 genes were upregulated in psoriatic skin, only LCE3A showed significant overexpression with homozygous risk genotype compared to the non-risk genotype. LCE3B also showed significant overexpression in patients with HLA-Cw6 allele. Moreover, LCE3A showed significantly higher expression in patients bearing homozygous risk genotype in presence of HLA-Cw6 allele but not in those having non-risk genotype, demonstrating the combined effect of HLA-Cw6 allele and risk associated genotype near LCE3A gene. Integration of genetic and gene expression data thus allowed us to identify the actual disease variants at the LCE3 cluster among the psoriasis patients in India.