greenPipes: an integrated data analysis pipeline for greenCUT&RUN and CUT&RUN genome-localization datasets.

MOTIVATION: To study gene regulation through transcription factors and chromatin modifiers, a variety of genome-wide techniques are used. Recently, CUT&RUN-based technologies have become popular, but a pipeline for the comprehensive analysis of CUT&RUN datasets is currently lacking. Here, we present the "greenPipes" package, which includes fine-tuned parameters specifically for bioinformatic analyses of greenCUT&RUN and CUT&RUN datasets. greenPipes provides additional functionalities for data analysis and data integration with other -omics technologies, which are either not available in other pipelines developed for CUT&RUN datasets or scattered in the literature as individual packages. AVAILABILITY AND IMPLEMENTATION: Source code and a manual of the greenPipes are freely available on GitHub website (https://github.com/snizam001/greenPipes). The test datasets, comprehensive annotation files, and other datasets are available at https://osf.io/ruhj9/. CONTACT: n.sheikh@dkfz-heidelberg.de or m.timmers@dkfz-heidelberg.de. SUPPLEMENTARY INFORMATION: The handbook of greenPipes is available online at Bioinformatics as Supplementary text.

Epigenetic-focused CRISPR/Cas9 screen identifies (absent, small, or homeotic)2-like protein (ASH2L) as a regulator of glioblastoma cell survival.

BACKGROUND: Glioblastoma is the most common and aggressive primary brain tumor with extremely poor prognosis, highlighting an urgent need for developing novel treatment options. Identifying epigenetic vulnerabilities of cancer cells can provide excellent therapeutic intervention points for various types of cancers. METHOD: In this study, we investigated epigenetic regulators of glioblastoma cell survival through CRISPR/Cas9 based genetic ablation screens using a customized sgRNA library EpiDoKOL, which targets critical functional domains of chromatin modifiers. RESULTS: Screens conducted in multiple cell lines revealed ASH2L, a histone lysine methyltransferase complex subunit, as a major regulator of glioblastoma cell viability. ASH2L depletion led to cell cycle arrest and apoptosis. RNA sequencing and greenCUT&RUN together identified a set of cell cycle regulatory genes, such as TRA2B, BARD1, KIF20B, ARID4A and SMARCC1 that were downregulated upon ASH2L depletion. Mass spectrometry analysis revealed the interaction partners of ASH2L in glioblastoma cell lines as SET1/MLL family members including SETD1A, SETD1B, MLL1 and MLL2. We further showed that glioblastoma cells had a differential dependency on expression of SET1/MLL family members for survival. The growth of ASH2L-depleted glioblastoma cells was markedly slower than controls in orthotopic in vivo models. TCGA analysis showed high ASH2L expression in glioblastoma compared to low grade gliomas and immunohistochemical analysis revealed significant ASH2L expression in glioblastoma tissues, attesting to its clinical relevance. Therefore, high throughput, robust and affordable screens with focused libraries, such as EpiDoKOL, holds great promise to enable rapid discovery of novel epigenetic regulators of cancer cell survival, such as ASH2L. CONCLUSION: Together, we suggest that targeting ASH2L could serve as a new therapeutic opportunity for glioblastoma. Video Abstract.

Integrating quantitative proteomics with accurate genome profiling of transcription factors by greenCUT&RUN.

Genome-wide localization of chromatin and transcription regulators can be detected by a variety of techniques. Here, we describe a novel method 'greenCUT&RUN' for genome-wide profiling of transcription regulators, which has a very high sensitivity, resolution, accuracy and reproducibility, whilst assuring specificity. Our strategy begins with tagging of the protein of interest with GFP and utilizes a GFP-specific nanobody fused to MNase to profile genome-wide binding events. By using a GFP-nanobody the greenCUT&RUN approach eliminates antibody dependency and variability. Robust genomic profiles were obtained with greenCUT&RUN, which are accurate and unbiased towards open chromatin. By integrating greenCUT&RUN with nanobody-based affinity purification mass spectrometry, 'piggy-back' DNA binding events can be identified on a genomic scale. The unique design of greenCUT&RUN grants target protein flexibility and yields high resolution footprints. In addition, greenCUT&RUN allows rapid profiling of mutants of chromatin and transcription proteins. In conclusion, greenCUT&RUN is a widely applicable and versatile genome-mapping technique.

A Complete Absence of Missense Mutation in Myosin Regulatory and Essential Light Chain Genes of South Indian Hypertrophic and Dilated Cardiomyopathies.

BACKGROUND: Myosin is a hexameric contractile protein composed of 2 heavy chains associated with 4 light chains of 2 distinct classes - 2 regulatory light chains (MYL2) and 2 essential light chains (MYL3). The myosin light chains stabilize the long alpha helical neck of the myosin head and regulate the myosin ATPase activities. OBJECTIVES: Mutations in MYL2 and MYL3 are reported to be associated with cardiomyopathies. However, there is no study available on these genes in Indian cardiomyopathies, and therefore we planned to study them. METHOD: For the first time we sequenced MYL2 and MYL3 genes in a total of 248 clinically well-characterized cardiomyopathies consisting of 101 hypertrophic and 147 dilated cases along with 207 healthy controls from south India. RESULTS: Our study revealed a total of 10 variations - 7 in MYL2 and 3 in MYL3, of which 3 are novel variations observed exclusively in cases. However, the 15 causative missense mutations previously reported are totally absent in our study, which showed that the sequences of MYL2 and MYL3 are highly conserved in Indian cases/controls. CONCLUSIONS: MYL2 and MYL3 mutations are rare and the least cause of cardiomyopathies in Indians.

c.620C>T mutation in GATA4 is associated with congenital heart disease in South India.

BACKGROUND: Congenital heart diseases (CHDs) usually refer to abnormalities in the structure and/or function of the heart that arise before birth. GATA4 plays an important role in embryonic heart development, hence the aim of this study was to find the association of GATA4 mutations with CHD among the south Indian CHD patients. METHOD: GATA4 gene was sequenced in 100 CHD patients (ASD, VSD, TOF and SV) and 200 controls. Functional significance of the observed GATA4 mutations was analyzed using PolyPhen, SIFT, PMut, Plink, Haploview, ESE finder 3.0 and CONSITE. RESULTS: We observed a total of 19 mutations, of which, one was in 5' UTR, 10 in intronic regions, 3 in coding regions and 5 in 3' UTR. Of the above mutations, one was associated with Atrial Septal Defect (ASD), two were found to be associated with Tetralogy of Fallot (TOF) and three (rs804280, rs4841587 and rs4841588) were strongly associated with Ventricular Septal Defect (VSD). Interestingly, one promoter mutation (-490 to 100 bp) i.e., 620 C>T (rs61277615, p-value = 0.008514), one splice junction mutation (G>A rs73203482; p-value = 9.6e-3, OR = 6.508) and one intronic mutation rs4841587 (p-value = 4.6e-3, OR = 4.758) were the most significant findings of this study. In silico analysis also proves that some of the mutations reported above are pathogenic. CONCLUSION: The present study found that GATA4 genetic variations are associated with ASD, TOF and VSD in South Indian patients. In silico analysis provides further evidence that some of the observed mutations are pathogenic.

TAS2R38 bitter taste perception in the Konkani Sarasvata Brahmin population.

BACKGROUND: The TAS2R38 gene carries markers for phenylthiocarbamide (PTC) sensitivity. Various studies have investigated the genotype-phenotype association pattern for bitter tasting ability and other factors in different populations. However, a paucity of such information for endogamous Indian populations is the reason behind this study. OBJECTIVE: To study the association of phenylthiocarbamide (PTC) sensitivity with TAS2R38 gene variations in Konkani Sarasvata Brahmin population. METHODS: We studied the association of the alleles rs714598, rs1726866, rs10246939 with PTC sensitivity and other factors in the Konkani Sarasvata Brahmin population. DNA was extracted from 114 individuals belonging to the Konkani Sarasvata Brahmin community. The TAS2R38 gene was sequenced to find the genotype distribution pattern. The association between genotype and phenotype was checked using the Chi-Square test and multifactorial logistical regression. RESULTS: We observed a 58.8% frequency of the AVI haplotype, which is the most prevalent in European populations. A higher number of non-taster haplotypes and diplotypes were observed in Konkani Sarasvata Brahmins, with the allele rs10246939 showing a significant association with PTC bitter taste sensitivity in both allelic (p = 8.6 × 10-4; Allele-G, OR = 3.57 [95% CI = 1.66-7.69]) and genotype-based (p = 6.9 × 10-4; genotype-AG, OR = 3.11 [95% CI = 0.73-13.20]; genotype-GG, OR = 40 [95% CI = 3.58-447.03]) tests. CONCLUSION: Our results are in line with earlier studies, which report an association between PTC sensitivity and the TAS2R38 gene in different populations. In the global context, Konkani Sarasvata Brahmins, who are mostly distributed along the southwestern coast of India, show a PTC sensitivity pattern slightly similar to that of West Eurasian populations. Our findings suggest ancestry specific selection in TAS2R38 gene variations for taste sensitivity at global level.

Maternal Transmission of the PAX7 Single Nucleotide Polymorphisms among Indian Cleft Trios.

Cleft lip and/or cleft palate (CL/P) is one of the most common congenital anomalies of the human face with a complex etiology involving multiple genetic and environmental factors. Several studies have shown the association of the paired box 7 ( PAX7 ) gene with CL/P in different populations worldwide. However, the current literature reveals no reported case-parent trio studies to evaluate the association between the PAX7 gene and the risk of nonsyndromic cleft lip and/or palate (NSCL/P) in the Indian population. Hence, the purpose of this study was to assess the PAX7 gene single nucleotide polymorphisms (SNPs) in the etiology of NSCL/P among the Indian cleft trios. Forty Indian case-parent trios of NSCL/P were included. The cases and their parents' genomic DNA were extracted. The SNPs rs9439714, rs1339062, rs6695765, rs742071, and rs618941of the PAX7 gene were genotyped using the Agena Bio MassARRAY analysis. The allelic transmission disequilibrium test was performed using PLINK software while pair-wise linkage disequilibrium by the Haploview program. The SNP rs9439714 showed evidence of association ( p -value = 0.02, odds ratio = 3) with NSCL/P. Considering the parent-of-origin effects, the SNPs rs9439714 and rs618941 showed an excess maternal transmission of allele C at rs9439714 ( p -value = 0.05) and G allele at rs618941 ( p -value = 0.04). The results of the present study suggested that the SNPs rs9439714 and rs618941 showed an excess maternal transmission of alleles suggestive of the possible role of the PAX7 gene involvement in the etiology of NSCL/P in the Indian population.

Alternative mRNA Splicing Controls the Functions of the Histone H3K27 Demethylase UTX/KDM6A.

The UTX/KDM6A histone H3K27 demethylase plays an important role in development and is frequently mutated in cancers such as urothelial cancer. Despite many studies on UTX proteins, variations in mRNA splicing have been overlooked. Using Nanopore sequencing, we present a comprehensive analysis of UTX/KDM6A splicing events in human cell lines and in tissue samples from bladder cancer cases and normal epithelia. We found that the central region of UTX mRNAs encoded by exons 12 to 17 undergoes extensive alternative splicing. Up to half of all stable mRNAs (8-48% in bladder tissues and 18-58% in cell lines) are represented by the UTX canonical isoform lacking exon 14 encoding a nuclear localization sequence, and hence exon 14-containing UTX isoforms exclusively localize to the nucleus, unlike the cytonuclear localization of the canonical isoform. Chromatin association was also higher for exon-14-containing isoforms compared to the canonical UTX. Using quantitative mass spectrometry, we found that all UTX isoforms integrated into the MLL3 and MLL4, PR-DUB and MiDAC complexes. Interestingly, one of the novel UTX isoforms, which lacks exons 14 and 16, fails to interact with PR-DUB and MiDAC complex members. In conclusion, UTX mRNAs undergo extensive alternative splicing, which controls the subcellular localization of UTX and its interactions with other chromatin regulatory complexes.

Multi-omics analyses of MEN1 missense mutations identify disruption of menin-MLL and menin-JunD interactions as critical requirements for molecular pathogenicity.

BACKGROUND: Loss-of-function mutations of the multiple endocrine neoplasia type 1 (MEN1) gene are causal to the MEN1 tumor syndrome, but they are also commonly found in sporadic pancreatic neuroendocrine tumors and other types of cancers. The MEN1 gene product, menin, is involved in transcriptional and chromatin regulation, most prominently as an integral component of KMT2A/MLL1 and KMT2B/MLL2 containing COMPASS-like histone H3K4 methyltransferase complexes. In a mutually exclusive fashion, menin also interacts with the JunD subunit of the AP-1 and ATF/CREB transcription factors. RESULTS: Here, we applied and in silico screening approach for 253 disease-related MEN1 missense mutations in order to select a set of nine menin mutations in surface-exposed residues. The protein interactomes of these mutants were assessed by quantitative mass spectrometry, which indicated that seven of the nine mutants disrupt interactions with both MLL1/MLL2 and JunD complexes. Interestingly, we identified three missense mutations, R52G, E255K and E359K, which predominantly reduce the MLL1 and MLL2 interactions when compared with JunD. This observation was supported by a pronounced loss of binding of the R52G, E255K and E359K mutant proteins at unique MLL1 genomic binding sites with less effect on unique JunD sites. CONCLUSIONS: Our results underline the effects of MEN1 gene mutations in both familial and sporadic tumors of endocrine origin on the interactions of menin with the MLL1 and MLL2 histone H3K4 methyltransferase complexes and with JunD-containing transcription factors. Menin binding pocket mutants R52G, E255K and E359K have differential effects on MLL1/MLL2 and JunD interactions, which translate into differential genomic binding patterns. Our findings encourage future studies addressing the pathophysiological relevance of the separate MLL1/MLL2- and JunD-dependent functions of menin mutants in MEN1 disease model systems.

A Box of Chemistry to Inhibit the MEN1 Tumor Suppressor Gene Promoting Leukemia.

Targeting protein-protein interactions (PPIs) with small-molecule inhibitors has become a hotbed of modern drug development. In this review, we describe a new class of PPI inhibitors that block menin from binding to MLL proteins. Menin is encoded by the MEN1 tumor suppressor, but acts as an essential cofactor for MLL/KMT2A-rearranged leukemias. The most promising menin-MLL inhibitors belong to the thienopyrimidine class and have recently entered phase I/II clinical trials for treating acute leukemias characterized by MLL/KMT2A translocations or NPM1 mutations. As single agents, thienopyrimidine compounds eradicate leukemia in a xenograft models of primary leukemic cells belonging to the MLL-rearranged or NPM1-mutant subtypes. These compounds are well tolerated with few or no side effects, which is remarkable given the tumor-suppressor function of menin. The menin-MLL inhibitors highlight how leukemia patients could benefit from a targeted epigenetic therapy with novel PPI inhibitors obtained by directed chemical evolution.

CYP2C9 Variations and Their Pharmacogenetic Implications Among Diverse South Asian Populations.

INTRODUCTION: Allelic frequency distribution of drug metabolizing enzyme genes among populations is important to identify risk groups for adverse drug reaction and to select representative populations for clinical trials. Although India emerged as an important hub for clinical trials, information about the pharmacogenetic diversity for this region is still lacking. Here, we investigated genetic diversity of cytochrome-P450-2C9 (CYP2C9) gene which metabolizes wide range of drugs and is highly expressed in the human liver. METHODS: In total, 1278 individuals from 36 diverse Indian populations, 210 individuals from in-house data-repository and 489 other South Asian samples from the 1000 Genomes Project were selected. Variants observed in CYP2C9 gene were subjected to various statistical analyses. RESULTS: High frequency of CYP2C9*3 (~13%) and CYP2C9*3/*3 (~1%) was observed among South Asians, compared to 21 populations living outside the Indian subcontinent. The allelic/genotypic frequency does not correlate with geographical location or linguistic affiliation, except populations speaking Tibeto-Burmans language, who have lower frequency of CYP2C9*3 and CYP2C9*3/*3. Since, South Asians practice strict endogamy, presence of unique mutation and high frequency of homozygous genotypes not surprising. CYP2C9*3 has been associated with therapeutic response.The effect of CYP2C9*3/*3 is more pronounced compared to heterozygous and wild type homozygous genotypes as evident in many in vitro studies. As South Asians have high frequency, it would be interesting to explore potential of CYP2C9*3 as a marker for personalized therapy. Our study revealed several rare functional variants, which form eight novel and rare haplotypes of CYP2C9 (CYP2C9*63-*70). Of which, CYP2C9*64, *65, *66, *68, *69 and *70 haplotypes are South Asian-specific. CONCLUSION: Overall, we find high genetic heterogeneity within South Asians and identified South Asian-specific putative functional CYP2C9 haplotypes. High frequency of CYP2C9*3 and CYP2C9*3/*3 was observed in South Asian populations. Taken together, current study greatly enriches the knowledge of naturally occurring CYP2C9 variants and its diversity in South Asia, which are relevant to further CYP2C9-related functional research and for personalized medicine.

Pan-Cancer Analysis Shows Enrichment of Macrophages, Overexpression of Checkpoint Molecules, Inhibitory Cytokines, and Immune Exhaustion Signatures in EMT-High Tumors.

Epithelial-mesenchymal transition (EMT) is a highly dynamic process that occurs under normal circumstances; however, EMT is also known to play a central role in tumor progression and metastasis. Furthermore, role of tumor immune microenvironment (TIME) in shaping anticancer immunity and inducing the EMT is also well recognized. Understanding the key features of EMT is critical for the development of effective therapeutic interventions. Given the central role of EMT in immune escape and cancer progression and treatment, we have carried out a pan-cancer TIME analysis of The Cancer Genome Atlas (TCGA) dataset in context to EMT. We have analyzed infiltration of various immune cells, expression of multiple checkpoint molecules and cytokines, and inflammatory and immune exhaustion gene signatures in 22 cancer types from TCGA dataset. A total of 16 cancer types showed a significantly increased (p < 0.001) infiltration of macrophages in EMT-high tumors (mesenchymal samples). Furthermore, out of the 17 checkpoint molecules we analyzed, 11 showed a significant overexpression (p < 0.001) in EMT-high samples of at least 10 cancer types. Analysis of cytokines showed significant enrichment of immunosuppressive cytokines-TGFB1 and IL10-in the EMT-high group of almost all cancer types. Analysis of various gene signatures showed enrichment of inflammation, exhausted CD8+ T cells, and activated stroma signatures in EMT-high tumors. In summary, our pan-cancer EMT analysis of TCGA dataset shows that the TIME of EMT-high tumors is highly immunosuppressive compared to the EMT-low (epithelial) tumors. The distinctive features of EMT-high tumors are as follows: (i) the enrichment of tumor-associated macrophages, (ii) overexpression of immune checkpoint molecules, (iii) upregulation of immune inhibitory cytokines TGFB1 and IL10, and (iv) enrichment of inflammatory and exhausted CD8+ T-cell signatures. Our study shows that TIMEs of different EMT groups differ significantly, and this would pave the way for future studies analyzing and targeting the TIME regulators for anticancer immunotherapy.

Novel FCN2 Variants and Haplotypes are Associated with Rheumatic Heart Disease.

Ficolins are pattern recognition molecules that are involved in innate immune defense. Ficonin-2 (FCN2) has a specific affinity for lipoteichoic acid present in the cell wall of Streptococcus pyogenes, an etiological agent for rheumatic heart disease (RHD). We have estimated FCN2 serum levels and analyzed the functional variants of FCN2 in 400 RHD patients, 617 healthy controls, and 581 individuals belonged to various ethnic populations, who are inhabited in various geographical regions of India. Our study revealed that the FCN2 -986A and +6359T alleles were the risk factors for RHD susceptibility (p = 0.0007 for -986G>A; p = 0.0004 for +6359C>T). The haplotype AGGT (p = 0.0024) was observed to be a risk factor for RHD susceptibility, and the haplotype GGAC (p = 0.002) was found to confer protection against RHD. The level of serum FCN2 was significantly higher in controls (p < 0.0001) and in controls with GGAC haplotypes (p < 0.0001). The frequency of the risk alleles -986A and +6359T was found to be more prevalent in Northern and North-Western (Indo-European) India. The protective GGAC haplotype was found more prevalent in Eastern (Tibeto-Burman) and South-Western (Dravidian) India. Alleles -986A and +6359T were in positive correlation with the prevalence of RHD (regression coefficient = 1.84 and 1.94, respectively), whereas GGAC haplotype was in negative correlation with prevalence of RHD (regression coefficient = -1.71). In conclusion, we found that low level of serum ficolin-2 is significantly associated with RHD. Further, FCN2 -986A and +6359T alleles and AGGT haplotype are associated with increased susceptibility to RHD, while GGAC haplotype is associated with moderate protection against RHD.

Dissection of TAF1 neuronal splicing and implications for neurodegeneration in X-linked dystonia-parkinsonism.

X-linked dystonia-parkinsonism (XDP) is a monogenic neurodegenerative disorder of the basal ganglia, which presents as a combination of hyperkinetic movements and parkinsonian features. The underlying genetic mechanism involves the insertion of a SINE-VNTR-Alu retrotransposon within the TAF1 gene. Interestingly, alterations of TAF1 have been involved in multiple neurological diseases. In XDP, the SINE-VNTR-Alu insertion in TAF1 has been proposed to result in alternative splicing defects, including the decreased incorporation of a neuron-specific microexon annotated as 34'. This mechanism has become controversial as recent studies failed to provide support. In order to resolve this conundrum, we examined the alternative splicing patterns of TAF1 mRNAs in XDP and control brains. The impact of the disease-associated SINE-VNTR-Alu on alternative splicing of microexon 34' was further investigated in cellular assays. Subsequently, microexon 34' incorporation was explored by RT-PCR and Nanopore long-read sequencing of TAF1 mRNAs from XDP and control brains tissues. Using cell-based splicing assays, we demonstrate that presence of the disease-associated SINE-VNTR-Alu does not affect the inclusion of microexon 34'. In addition, we show that (1) microexon 34'-containing TAF1 mRNAs are detected at similar levels in XDP as in controls and that (2) the architecture of TAF1 transcripts is remarkably similar between XDP and controls brains. These results indicate that microexon 34' incorporation into TAF1 mRNA is not affected in XDP brains. Our findings shift the current paradigm of XDP by discounting alternative splicing of TAF1 microexon 34' as the molecular basis for this disease.

Variations in macrophage migration inhibitory factor gene are not associated with visceral leishmaniasis in India.

BACKGROUND: The host genetic factors play important role in determining the outcome of visceral leishmaniasis (VL). Macrophage migration inhibitory factor (MIF) is an important host cytokine, which is a key regulator of innate immune system. Genetic variants in MIF gene have been found to be associated with several inflammatory and infectious diseases. Role of MIF is well documented in leishmaniasis diseases, including Indian visceral leishmaniasis, where elevated level of serum MIF has been associated with VL phenotypes. However, there was no genetic study to correlate MIF variants in VL, therefore, we aimed to study the possible association of three reported MIF gene variants -794 CATT, -173G > C and non-coding RNA gene LOC284889 in Indian VL phenotype. METHODS: Study subjects comprised of 214 VL patients along with ethnically and demographically matched 220 controls from VL endemic regions of Bihar state in India. RESULTS: We found no significant difference between cases and controls in allelic, genotypic and haplotype frequency of the markers analysed [-794 CATT repeats (χ2=0.86; p=0.35; OR=0.85; 95% CI=0.61-1.19); -173 G>C polymorphism (χ2=1.11; p=0.29; OR=0.83; 95% CI=0.59-1.16); and LOC284889 (χ2=0.78; p=0.37; OR=0.86; 95% CI=0.61-1.20)]. CONCLUSION: Since we did not find any significant differences between case and control groups, we conclude that sequencing of complete MIF gene and extensive study on innate and adaptive immunity genes may help in identifying genetic variations that are associated with VL susceptibility/resistance among Indians.

TG haplotype in the LRP8 is associated with myocardial infarction in south Indian population.

Myocardial infarction (MI) is a complex multifactorial cardiovascular disease. India experiences a much greater burden of MI, also suggesting an experimental increase of this burden in the future. The absolute reasons for MI are context dependent and differ with different geographical settings. Several reports indicate that SNPs that are associated with certain diseases in other populations may not be associated with Indian population. It is, therefore, important to validate the association of SNPs. Low density lipoprotein receptor related protein 8 (LRP8) gene plays central role in human lipoprotein metabolism as it facilitates the clearance of bad cholesterol LDL, VLDL from plasma and is reported to be associated with MI in the western population. However, this gene has not been studied in the South Indian population. We aim to test the role of the LRP8 gene variants correlating with the lipid profile in MI patients in South Indian population. We sequenced regions of SNPs rs10788952, rs7546246, rs2297660 and rs5174 of LRP8 in 100 MI patients and 100 age-matched controls. Our result revealed a total of 4 variations. None of the SNPs were significantly associated with MI (p>0.973). Interestingly, haplotype based association analysis showed TG and CG of rs10788952 and rs7546246 significantly associated with MI (p<0.01 and p<0.00005) and in particular, haplotype TG was positively correlated with the risk of MI, as this increased the LDL and total cholesterol level in MI patients in south Indians. Our results suggest that haplotype TG is a risk factor for MI in South Indian population.

c.*84G>A Mutation in CETP Is Associated with Coronary Artery Disease in South Indians.

BACKGROUND: Coronary artery disease (CAD) is one of the leading causes of mortality worldwide. It is a multi-factorial disease and several studies have demonstrated that the genetic factors play a major role in CAD. Although variations in cholesteryl ester transfer protein (CETP) gene are reported to be associated with CAD, this gene has not been studied in South Indian populations. Hence we evaluated the CETP gene variations in CAD patients of South Indian origin. METHODS: We sequenced all the exons, exon-intron boundaries and UTRs of CETP in 323 CAD patients along with 300 ethnically and age matched controls. Variations observed in CETP were subjected to various statistical analyses. RESULTS AND DISCUSSION: Our analysis revealed a total of 13 variations. Of these, one3'UTRvariant rs1801706 (c.*84G>A) was significantly associated with CAD (genotype association test: OR = 2.16, 95% CI: 1.50-3.10, p = 1.88x10-5 and allelic association test: OR = 1.92, 95% CI: 1.40-2.63, p = 2.57x10-5). Mutant allele "A" was observed to influence the higher concentration of mRNA (p = 7.09×10-3, R2 = 0.029 and ß = 0.2163). Since expression of CETP has been shown to be positively correlated with the risk of CAD, higher frequency of "A" allele (patients: 22.69% vs.controls: 13%) reveals that c.*84G>A is a risk factor for CAD in South Indians. CONCLUSIONS: This is the first report of the CETP gene among South Indians CAD patients. Our results suggest that rs1801706 (c.*84G>A) is a risk factor for CAD in South Indian population.

Dissecting the influence of Neolithic demic diffusion on Indian Y-chromosome pool through J2-M172 haplogroup.

The global distribution of J2-M172 sub-haplogroups has been associated with Neolithic demic diffusion. Two branches of J2-M172, J2a-M410 and J2b-M102 make a considerable part of Y chromosome gene pool of the Indian subcontinent. We investigated the Neolithic contribution of demic dispersal from West to Indian paternal lineages, which majorly consists of haplogroups of Late Pleistocene ancestry. To accomplish this, we have analysed 3023 Y-chromosomes from different ethnic populations, of which 355 belonged to J2-M172. Comparison of our data with worldwide data, including Y-STRs of 1157 individuals and haplogroup frequencies of 6966 individuals, suggested a complex scenario that cannot be explained by a single wave of agricultural expansion from Near East to South Asia. Contrary to the widely accepted elite dominance model, we found a substantial presence of J2a-M410 and J2b-M102 haplogroups in both caste and tribal populations of India. Unlike demic spread in Eurasia, our results advocate a unique, complex and ancient arrival of J2a-M410 and J2b-M102 haplogroups into Indian subcontinent.

IL10 Variant g.5311A Is Associated with Visceral Leishmaniasis in Indian Population.

BACKGROUND: Visceral leishmaniasis (VL) is a multifactorial disease, where the host genetics play a significant role in determining the disease outcome. The immunological role of anti-inflammatory cytokine, Interleukin 10 (IL10), has been well-documented in parasite infections and considered as a key regulatory cytokine for VL. Although VL patients in India display high level of IL10 in blood serum, no genetic study has been conducted to assess the VL susceptibility / resistance. Therefore, the aim of this study is to investigate the role of IL10 variations in Indian VL; and to estimate the distribution of disease associated allele in diverse Indian populations. METHODOLOGY: All the exons and exon-intron boundaries of IL10 were sequenced in 184 VL patients along with 172 ethnically matched controls from VL endemic region of India. RESULT AND DISCUSSION: Our analysis revealed four variations; rs1518111 (2195 A>G, intron), rs1554286 (2607 C>T, intron), rs3024496 (4976 T>C, 3' UTR) and rs3024498 (5311 A>G, 3' UTR). Of these, a variant g.5311A is significantly associated with VL (χ2=18.87; p =0.00001). In silico approaches have shown that a putative micro RNA binding site (miR-4321) is lost in rs3024498 mRNA. Further, analysis of the above four variations in 1138 individuals from 34 ethnic populations, representing different social and linguistic groups who are inhabited in different geographical regions of India, showed variable frequency. Interestingly, we have found, majority of the tribal populations have low frequency of VL ('A' of rs3024498); and high frequency of leprosy ('T' of rs1554286), and Behcet's ('A' of rs1518111) associated alleles, whereas these were vice versa in castes. Our findings suggest that majority of tribal populations of India carry the protected / less severe allele against VL, while risk / more severe allele for leprosy and Behcet's disease. This study has potential implications in counseling and management of VL and other infectious diseases.