Identifying recent adaptations in large-scale genomic data.

Although several hundred regions of the human genome harbor signals of positive natural selection, few of the relevant adaptive traits and variants have been elucidated. Using full-genome sequence variation from the 1000 Genomes (1000G) Project and the composite of multiple signals (CMS) test, we investigated 412 candidate signals and leveraged functional annotation, protein structure modeling, epigenetics, and association studies to identify and extensively annotate candidate causal variants. The resulting catalog provides a tractable list for experimental follow-up; it includes 35 high-scoring nonsynonymous variants, 59 variants associated with expression levels of a nearby coding gene or lincRNA, and numerous variants associated with susceptibility to infectious disease and other phenotypes. We experimentally characterized one candidate nonsynonymous variant in Toll-like receptor 5 (TLR5) and show that it leads to altered NF-κB signaling in response to bacterial flagellin. PAPERFLICK:

Pyruvate kinase M knockdown-induced signaling via AMP-activated protein kinase promotes mitochondrial biogenesis, autophagy, and cancer cell survival.

Preferential expression of the low-activity (dimeric) M2 isoform of pyruvate kinase (PK) over its constitutively active splice variant M1 isoform is considered critical for aerobic glycolysis in cancer cells. However, our results reported here indicate co-expression of PKM1 and PKM2 and their possible physical interaction in cancer cells. We show that knockdown of either PKM1 or PKM2 differentially affects net PK activity, viability, and cellular ATP levels of the lung carcinoma cell lines H1299 and A549. The stable knockdown of PK isoforms in A549 cells significantly reduced the cellular ATP level, whereas in H1299 cells the level of ATP was unaltered. Interestingly, the PKM1/2 knockdown in H1299 cells activated AMP-activated protein kinase (AMPK) signaling and stimulated mitochondrial biogenesis and autophagy to maintain energy homeostasis. In contrast, knocking down either of the PKM isoforms in A549 cells lacking LKB1, a serine/threonine protein kinase upstream of AMPK, failed to activate AMPK and sustain energy homeostasis and resulted in apoptosis. Moreover, in a similar genetic background of silenced PKM1 or PKM2, the knocking down of AMPKα1/2 catalytic subunit in H1299 cells induced apoptosis. Our findings help explain why previous targeting of PKM2 in cancer cells to control tumor growth has not met with the expected success. We suggest that this lack of success is because of AMPK-mediated energy metabolism rewiring, protecting cancer cell viability. On the basis of our observations, we propose an alternative therapeutic strategy of silencing either of the PKM isoforms along with AMPK in tumors.

WITHDRAWN: Regulation of Pyruvate Kinase M Switch towards PKM1 by LKB1-AMPK Axis Tolerates Hypoglycemic Stress in Cancer Cells.

This article has been withdrawn by the authors. The PKM2 immunoblot in Fig 2E was reused as part of the Caspase-3 immunoblot in Fig 9C. The PKM2 immunoblot from 5 mM Glu, fractions 1-10 was reused as the PKM2 immunoblot from 1 mM Glu, fractions 1-10. The actin immunoblot from A549 cells from Fig 5A was reused as the actin blot from Fig 7C.

Apoptosis regulatory protein-protein interaction demonstrates hierarchical scale-free fractal network.

Dysregulation or inhibition of apoptosis favors cancer and many other diseases. Understanding of the network interaction of the genes involved in apoptotic pathway, therefore, is essential, to look for targets of therapeutic intervention. Here we used the network theory methods, using experimentally validated 25 apoptosis regulatory proteins and identified important genes for apoptosis regulation, which demonstrated a hierarchical scale-free fractal protein-protein interaction network. TP53, BRCA1, UBIQ and CASP3 were recognized as a four key regulators. BRCA1 and UBIQ were also individually found to control highly clustered modules and play an important role in the stability of the overall network. The connection among the BRCA1, UBIQ and TP53 proteins was found to be important for regulation, which controlled their own respective communities and the overall network topology. The feedback loop regulation motif was identified among NPM1, BRCA1 and TP53, and these crucial motif topologies were also reflected in high frequency. The propagation of the perturbed signal from hubs was found to be active upto some distance, after which propagation started decreasing and TP53 was the most efficient signal propagator. From the functional enrichment analysis, most of the apoptosis regulatory genes associated with cardiovascular diseases and highly expressed in brain tissues were identified. Apart from TP53, BRCA1 was observed to regulate apoptosis by influencing motif, propagation of signals and module regulation, reflecting their biological significance. In future, biochemical investigation of the observed hub-interacting partners could provide further understanding about their role in the pathophysiology of cancer.

ERK2-Pyruvate Kinase Axis Permits Phorbol 12-Myristate 13-Acetate-induced Megakaryocyte Differentiation in K562 Cells.

Metabolic changes that contribute to differentiation are not well understood. Overwhelming evidence shows the critical role of glycolytic enzyme pyruvate kinase (PK) in directing metabolism of proliferating cells. However, its role in metabolism of differentiating cells is unclear. Here we studied the role of PK in phorbol 12-myristate 13-acetate (PMA)-induced megakaryocytic differentiation in human leukemia K562 cells. We observed that PMA treatment decreased cancer-type anabolic metabolism but increased ATP production, along with up-regulated expression of two PK isoforms (PKM2 and PKR) in an ERK2-dependent manner. Interestingly, silencing of PK (PKM2 and PKR) inhibited PMA-induced megakaryocytic differentiation, as revealed by decreased expression of megakaryocytic differentiation marker CD61 and cell cycle behavior. Further, PMA-induced ATP production reduced greatly upon PK silencing, suggesting that PK is required for ATP synthesis. In addition to metabolic effects, PMA treatment also translocated PKM2, but not PKR, into nucleus. ERK1/2 knockdowns independently and together suggested the role of ERK2 in the up-regulation of both the isoforms of PK, proposing a role of ERK2-PK isoform axis in differentiation. Collectively, our findings unravel ERK2 guided PK-dependent metabolic changes during PMA induction, which are important in megakaryocytic differentiation.

Mapping of PARK2 and PACRG overlapping regulatory region reveals LD structure and functional variants in association with leprosy in unrelated indian population groups.

Leprosy is a chronic infectious disease caused by Mycobacterium Leprae, where the host genetic background plays an important role toward the disease pathogenesis. Various studies have identified a number of human genes in association with leprosy or its clinical forms. However, non-replication of results has hinted at the heterogeneity among associations between different population groups, which could be due to differently evolved LD structures and differential frequencies of SNPs within the studied regions of the genome. A need for systematic and saturated mapping of the associated regions with the disease is warranted to unravel the observed heterogeneity in different populations. Mapping of the PARK2 and PACRG gene regulatory region with 96 SNPs, with a resolution of 1 SNP per 1 Kb for PARK2 gene regulatory region in a North Indian population, showed an involvement of 11 SNPs in determining the susceptibility towards leprosy. The association was replicated in a geographically distinct and unrelated population from Orissa in eastern India. In vitro reporter assays revealed that the two significantly associated SNPs, located 63.8 kb upstream of PARK2 gene and represented in a single BIN of 8 SNPs, influenced the gene expression. A comparison of BINs between Indian and Vietnamese populations revealed differences in the BIN structures, explaining the heterogeneity and also the reason for non-replication of the associated genomic region in different populations.

Missense mutations in pyruvate kinase M2 promote cancer metabolism, oxidative endurance, anchorage independence, and tumor growth in a dominant negative manner.

The present study was designed to examine the functional relevance of two heterozygous mutations (H391Y and K422R), observed earlier by us in the Bloom syndrome condition. Cells stably expressing exogenous wild-type or mutant PKM2 (K422R or H391Y) or co-expressing both wild type and mutant (PKM2-K422R or PKM2-H391Y) were assessed for cancer metabolism and tumorigenic potential. Interestingly, cells co-expressing PKM2 and mutant (K422R or H391Y) showed significantly aggressive cancer metabolism as compared with cells expressing either wild-type or mutant PKM2 independently. A similar trend was observed for oxidative endurance, tumorigenic potential, cellular proliferation, and tumor growth. These observations signify the dominant negative nature of mutations. Remarkably, PKM2-H391Y co-expressed cells showed a maximal effect on all the studied parameters. Such a dominant negative impaired function of PKM2 in tumor development is not known; this study demonstrates for the first time the possible predisposition of Bloom syndrome patients with impaired PKM2 activity to cancer and the importance of studying genetic variations in PKM2 in the future to understand their relevance in cancer in general.

Moderate DNA damage promotes metabolic flux into PPP via PKM2 Y-105 phosphorylation: a feature that favours cancer cells.

Pyruvate kinase M2, an important metabolic enzyme, promotes aerobic glycolysis (Warburg effect) to facilitate cancer cell proliferation. Unravelling the status of this important glycolytic pathway enzyme under sub-lethal doses of etoposide, a commonly used anti-proliferative genotoxic drug to induce mild/moderate DNA damage in HeLa cells as a model system and discern its effect on: PKM2 expression, phosphorylation, dimer: tetramer ratio, activity and associated effects, was pertinent. Protein expression and phosphorylation of PKM2 from HeLa cells was estimated using Western blotting. Same protein lysate was also used to estimate total pyruvate kinase activity and the total dimer: tetramer content evaluated using glycerol gradient ultra-centrifugation. Intracellular PEP was estimated manually using standard curve; while NADPH was assessed by NADPH estimation kit. Unpaired t test and two-way-ANOVA was used for statistical analysis. A relative decrease in PKM2 expression and a subsequent dose and time dependent increase in Y105-phosphorylation were observed. A concomitant increase in PKM2 dimer content and Y105-phosphorylation responsible for reduced PKM2 activity promoted PEP accumulation and NADPH production, representing increased metabolic flux into PPP, a feature that favours cancer cells. It was apparent that the sub-lethal doses of etoposide induced inadequate damage to DNA in cancer cells in culture promoted pro-survival conditions due to Y105-phosphorylation of PKM2, its stable dimerization and inactivation, a unique association not known earlier, indicating what might happen in tumour revivals or recurrences.

Insulin enhances metabolic capacities of cancer cells by dual regulation of glycolytic enzyme pyruvate kinase M2.

BACKGROUND: Insulin is tightly associated with cancer progression; however, mechanistic insights into such observations are poorly understood. Recent studies show that metabolic transformation is critical to cancer cell proliferation. Here, we attempt to understand the role of insulin in promotion of cancer metabolism. To this end, the role of insulin in regulating glycolytic enzyme pyruvate kinase M2 (PKM2) was examined. RESULTS: We observed that insulin up-regulated PKM2 expression, through PI3K/mTOR mediated HIF1α induction, but significantly reduced PKM2 activity independent of this pathway. Drop in PKM2 activity was attributed to subunit dissociation leading to formation of low activity PKM2 oligomers, as assessed by density gradient centrifugation. However, tyrosine 105 phosphorylation of PKM2, known for inhibiting PKM2 activity, remained unaffected on insulin treatment. Interestingly, insulin-induced ROS was found responsible for PKM2 activity reduction. The observed changes in PKM2 status led to augmented cancer metabolism. Insulin-induced PKM2 up-regulation resulted in enhanced aerobic glycolysis as confirmed by PKM2 knockdown studies. Further, PKM2 activity reduction led to characteristic pooling of glycolytic intermediates and increased accumulation of NADPH; suggesting diversion of glucose flux towards macromolecular synthesis, necessary for cancer cell growth. CONCLUSION: The study identifies new PKM2-mediated effects of insulin on cancer metabolism, thus, advancing the understanding of insulin's role in cancer.

Invasive Breast Cancer: miR-24-2 Targets Genes Associated with Survival and Sensitizes MDA-MB-231 Cells to Berberine.

MicroRNA aberrations including that of miR-24-2 have been reported in various cancers. However, the target genes for miR-24-2 are yet to be identified and validated in invasive breast cancer and the triple-negative breast cancer (TNBC). Using in silico approaches and gene expression analyses, we identified and validated the target genes of miR-24-2 in invasive breast cancer, majority of which were TNBC. We studied the translational potential of these target genes using berberine in a TNBC cell line. Differentially expressed genes targeted by miR-24-2 were identified and analyzed for their survival effects using the The Cancer Genome Atlas-Breast Invasive Carcinoma (-BRCA) samples. Furthermore, we carried out protein-protein interaction, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, gene expression, and Kaplan-Meier survival analyses using common targets of miR-24-2 in invasive breast cancer/TNBC. We identified 11 biomarker candidate genes as crucial targets of miR-24-2. The survival of breast cancer patients was significantly associated with the low expressions of nine genes, including RACGAP1, KIAA1199, TIMM17A, LYRM7, IL1R1, SLC1A3, DTX4, L1CAM, and SAP30-like (SAP30L), and high expressions of two genes, SOD2 and HLA-DQB2. These in silico findings were validated by overexpressing miR-24-2 and assessing the expression pattern of these target genes in the TNBC MDA-MB-231 cells. miR-24-2 overexpression inhibited (by 20%; p < 0.001) cell proliferation and sensitized the anticancer effect of berberine. In all, this study reports on the novel target genes of miR-24-2 in invasive breast cancer/TNBC, and that miR-24-2 sensitizes MDA-MB-231 cells to berberine. These data lend evidence for the translational potentials of miR-24-2 for invasive breast cancer diagnostic and therapeutic innovation.

Association of variants in BAT1-LTA-TNF-BTNL2 genes within 6p21.3 region show graded risk to leprosy in unrelated cohorts of Indian population.

Host immune response against Mycobacterium leprae plays an important role in providing resistance to infection and disease progression. Genome-wide linkage and association studies suggest the possibility of multiple risk loci within HLA (6p21.3) region. Any systematic study of relevance within the histocompatibility complex of importance in host immune response would be pertinent because of non-replication of the known loci and unavailable information on some of the unexplored genes and regions. A systematic scan was performed of the selected region involving LTA-TNF-LTB genes within 6p21.3 with a resolution of 1SNP/127 bp; and the SNPs in flanking BAT1, NFKBIL and BTNL2-DRA genes on the basis of their tag status or their presence in promoter/exonic regions with MAF of >5%. Nine SNPs located in BAT1, LTA, TNF genes and BTNL2-DRA interval showed strong association with leprosy susceptibility in two independent sets of North Indian population which was replicated in a geographically distinct East Indian population. Conditional logistic regression showed at least one functional SNP remaining significant in each gene, suggesting an independent role of each of the disease associated SNPs. In vitro reporter assay revealed that two SNPs located at BAT1 promoter and 13 kb upstream to LTA gene affected the transcription factor binding site, hence the gene expression. We unravel the role of unexplored immunologically important genes, BAT1 and BTNL2, in addition to known LTA and TNF genes, and the haplotypes of the significantly associated SNPs therein, to understand susceptibility to the disease, leprosy and its differential severity.

Leprosy and the adaptation of human toll-like receptor 1.

Leprosy is an infectious disease caused by the obligate intracellular pathogen Mycobacterium leprae and remains endemic in many parts of the world. Despite several major studies on susceptibility to leprosy, few genomic loci have been replicated independently. We have conducted an association analysis of more than 1,500 individuals from different case-control and family studies, and observed consistent associations between genetic variants in both TLR1 and the HLA-DRB1/DQA1 regions with susceptibility to leprosy (TLR1 I602S, case-control P = 5.7 x 10(-8), OR = 0.31, 95% CI = 0.20-0.48, and HLA-DQA1 rs1071630, case-control P = 4.9 x 10(-14), OR = 0.43, 95% CI = 0.35-0.54). The effect sizes of these associations suggest that TLR1 and HLA-DRB1/DQA1 are major susceptibility genes in susceptibility to leprosy. Further population differentiation analysis shows that the TLR1 locus is extremely differentiated. The protective dysfunctional 602S allele is rare in Africa but expands to become the dominant allele among individuals of European descent. This supports the hypothesis that this locus may be under selection from mycobacteria or other pathogens that are recognized by TLR1 and its co-receptors. These observations provide insight into the long standing host-pathogen relationship between human and mycobacteria and highlight the key role of the TLR pathway in infectious diseases.

Genetic variations and interactions in anti-inflammatory cytokine pathway genes in the outcome of leprosy: a study conducted on a MassARRAY platform.

BACKGROUND: Mycobacterium leprae is the etiologic pathogen that causes leprosy. The outcome of disease is dependent on the host genetic background. METHODS: We investigated the association of 51 single-nucelotide polymorphisms (SNPs) in anti-inflammatory cytokines (IL-10, TGFB1, IL-6, IL-4, and IL-13) and receptors (IL-10RA, IL-10RB, TGFBR1, TGFBR2, IL-6R, IL-4R, IL-5RA, IL-5RB, and IL-13RA1) with susceptibility to leprosy in a case-control study from New Delhi in northern India. This was followed by replication testing of associated SNPs in a geographically distinct and unrelated population from Orissa in eastern India. The functional potential of SNPs was established with in vitro reporter assays. RESULTS: Significant associations (P < .05) were observed for 8 polymorphisms (rs1800871, rs1800872, and rs1554286 of IL-10; rs3171425 and rs7281762 of IL-10RB; rs2228048 and rs744751 of TGFBR2; and rs1800797 of IL-6) with leprosy. This association was replicated for 4 SNPs (rs1554286 of IL-10, rs7281762 of IL-10RB, rs2228048 of TGFBR2, and rs1800797 of IL-6). The interaction study revealed a significantly greater association with leprosy risk than was obtained for any SNP individually. CONCLUSIONS: This study provides an interesting insight on the cumulative polygenic host component that regulates leprosy pathogenesis.

HPV18 oncoproteins driven expression of PKM2 reprograms HeLa cell metabolism to maintain aerobic glycolysis and viability.

The molecular basis of human papillomavirus (HPV)-mediated cellular immortalization and malignant transformation has illustrated an indispensable role of viral E6/E7-oncoproteins. However, the impact of viral-oncoproteins on the metabolic phenotype of cancer cells remains ambiguous. We showed silencing of HPV18-encoded E6/E7-oncoprotein significantly reduced glucose consumption, lactate production, ATP level and viability. Silencing of HPV18-encoded E6/E7 in HeLa cells significantly down-regulated expression and activity of HK1, HK2, LDHA, and LDHB. Interestingly, there was an increased pyruvate kinase activity due to switch in expression from PKM2 isoform to PKM1. The switch in favor of alternatively spliced isoform PKM1, was regulated by viral-E6/E7-oncoprotein by inhibiting the c-Myc/hnRNP-axis. Further, the near absence of the PKM1 protein despite an adequate amount of PKM1 mRNA in HeLa cells was due to its proteasomal degradation. Our results suggests HPV18-encoded E6/E7 driven preferential expression of PKM2 is essential to support aerobic glycolysis and cell proliferation. Supplementary Information: The online version contains supplementary material available at 10.1007/s13337-022-00776-w.

Dominant negative mutations affect oligomerization of human pyruvate kinase M2 isozyme and promote cellular growth and polyploidy.

This study was designed to understand the mechanism and functional implication of the two heterozygous mutations (H391Y and K422R) of human pyruvate kinase M2 isozyme (PKM(2)) observed earlier in a Bloom syndrome background. The co-expression of homotetrameric wild type and mutant PKM(2) in the cellular milieu resulting in the interaction between the two at the monomer level was substantiated further by in vitro experiments. The cross-monomer interaction significantly altered the oligomeric state of PKM(2) by favoring dimerization and heterotetramerization. In silico study provided an added support in showing that hetero-oligomerization was energetically favorable. The hetero-oligomeric populations of PKM(2) showed altered activity and affinity, and their expression resulted in an increased growth rate of Escherichia coli as well as mammalian cells, along with an increased rate of polyploidy. These features are known to be essential to tumor progression. This study provides insight in understanding the modulated role of large oligomeric multifunctional proteins such as PKM(2) by affecting cellular behavior, which is an essential observation to understand tumor sustenance and progression and to design therapeutic intervention in future.

miR-24-2 controls H2AFX expression regardless of gene copy number alteration and induces apoptosis by targeting antiapoptotic gene BCL-2: a potential for therapeutic intervention.

INTRODUCTION: New levels of gene regulation with microRNA (miR) and gene copy number alterations (CNAs) have been identified as playing a role in various cancers. We have previously reported that sporadic breast cancer tissues exhibit significant alteration in H2AX gene copy number. However, how CNA affects gene expression and what is the role of miR, miR-24-2, known to regulate H2AX expression, in the background of the change in copy number, are not known. Further, many miRs, including miR-24-2, are implicated as playing a role in cell proliferation and apoptosis, but their specific target genes and the pathways contributing to them remain unexplored. METHODS: Changes in gene copy number and mRNA/miR expression were estimated using real-time polymerase chain reaction assays in two mammalian cell lines, MCF-7 and HeLa, and in a set of sporadic breast cancer tissues. In silico analysis was performed to find the putative target for miR-24-2. MCF-7 cells were transfected with precursor miR-24-2 oligonucleotides, and the gene expression levels of BRCA1, BRCA2, ATM, MDM2, TP53, CHEK2, CYT-C, BCL-2, H2AFX and P21 were examined using TaqMan gene expression assays. Apoptosis was measured by flow cytometric detection using annexin V dye. A luciferase assay was performed to confirm BCL-2 as a valid cellular target of miR-24-2. RESULTS: It was observed that H2AX gene expression was negatively correlated with miR-24-2 expression and not in accordance with the gene copy number status, both in cell lines and in sporadic breast tumor tissues. Further, the cells overexpressing miR-24-2 were observed to be hypersensitive to DNA damaging drugs, undergoing apoptotic cell death, suggesting the potentiating effect of mir-24-2-mediated apoptotic induction in human cancer cell lines treated with anticancer drugs. BCL-2 was identified as a novel cellular target of miR-24-2. CONCLUSIONS: mir-24-2 is capable of inducing apoptosis by modulating different apoptotic pathways and targeting BCL-2, an antiapoptotic gene. The study suggests that miR-24-2 is more effective in controlling H2AX gene expression, regardless of the change in gene copy number. Further, the study indicates that combination therapy with miR-24-2 along with an anticancer drug such as cisplatin could provide a new avenue in cancer therapy for patients with tumors otherwise resistant to drugs.

Functional implication of TRAIL -716 C/T promoter polymorphism on its in vitro and in vivo expression and the susceptibility to sporadic breast tumor.

Recently, TRAIL function has been elucidated beyond its known classical role of mediating cellular homeostasis and immune surveillance against transformed cells. Here, we show how CC genotype of -716 TRAIL promoter SNP rendered risk for sporadic breast cancer as compared to the CT and TT genotypes (P (recessive model) = 0.018, OR = 1.4, 95% CI = 1.1-1.9; P (allele model) = 0.010, OR = 1.3, 95% CI = 1.1-1.7). The in silico prediction of the introduction of core Sp1/Sp3-binding motif suggested the functional significance of the SNP variation. This functional implication was validated by luciferase assay in HeLa (P = 0.026), MCF-7 (P = 0.022), HepG2 (P = 0.024), and HT1080 (P = 0.030) cells and also by real-time expression studies on tumor tissues (P = 0.01), revealing the transcriptionally repressed status of -716 T when compared to -716 C allele. The SNP-SNP interactions reflected an enhanced protective effect of CT and TT genotypes with the protective genetic backgrounds of TP53-BRCA2 (P = 0.002, OR = 0.2, 95% CI = 0.1-0.6), IFNG (P = 0.0000002, OR = 0.3, 95% CI = 0.2-0.4), and common variant Casp8 (P = 0.0003, OR = 0.5, 95% CI = 0.3-0.7). Interestingly, a comparison with clinical parameters showed overrepresented CT and TT genotypes in progressing (P = 0.041) and ER/PR negative tumors (P = 0.024/0.006). This was explained by increased apoptotic index, calculated as a ratio of selected pro-apoptotic and anti-apoptotic gene expression profiles, in CC genotyped tumors, favoring either intrinsic (P = 0.008,0.018) or extrinsic (P = 0.025,0.217) pathway depending upon the ER/PR status. Our study reveals for the first time that a promoter SNP of TRAIL functionally modulates the gene and consequently its role in breast cancer pathogenesis, cautioning to consider the -716 TRAIL SNP status in patients undergoing TRAIL therapy.

Y chromosome diversity among the Iranian religious groups: a reservoir of genetic variation.

BACKGROUND: Iran is ethnically, linguistically and religiously diverse. However, little is known about the population genetics of Iranian religious communities. AIM: This study was performed in order to define the different paternal components of the Iranian gene pool. SUBJECTS AND METHODS: Fourteen Y chromosome bi-allelic markers were analysed in 130 male subjects from Assyrian, Armenian and Zoroastrian groups in comparison with 208 male subjects from three Iranian Muslim groups. RESULTS: Among the three Iranian Muslim groups, the Uromian people possessed a particularly close genetic relationship to the Armenian, whereas the Zoroastrian group was different from the Uromian, but had a close genetic relationship to the two other Muslim groups (Kermanian and Shirazian). The genetic results indicate a relationship between Armenian and Assyrian groups in Iran and a clear distinction of the former from the Zoroastrian group. However, Assyrians had elevated frequency (40%) of R*(xR1a) and low frequency (11%) of J. CONCLUSION: The results of this study may suggest that the Assyrian population either experienced Eurasian gene flow (possibly from Armenia) or that enforced relocations and expulsion of conquered people with different origin led to the integration of descendants with R haplogroup. This could also be due to genetic drift due to small population size and endogamy resulting from religious barriers.

Investigation of DNA damage response and apoptotic gene methylation pattern in sporadic breast tumors using high throughput quantitative DNA methylation analysis technology.

BACKGROUND: Sporadic breast cancer like many other cancers is proposed to be a manifestation of abnormal genetic and epigenetic changes. For the past decade our laboratory has identified genes involved in DNA damage response (DDR), apoptosis and immunosurveillance pathways to influence sporadic breast cancer risk in north Indian population. Further to enhance our knowledge at the epigenetic level, we performed DNA methylation study involving 17 gene promoter regions belonging to DNA damage response (DDR) and death receptor apoptotic pathway in 162 paired normal and cancerous breast tissues from 81 sporadic breast cancer patients, using a high throughput quantitative DNA methylation analysis technology. RESULTS: The study identified five genes with statistically significant difference between normal and tumor tissues. Hypermethylation of DR5 (P=0.001), DCR1 (P=0.00001), DCR2 (P=0.0000000005) and BRCA2 (P=0.007) and hypomethylation of DR4 (P=0.011) in sporadic breast tumor tissues suggested a weak/aberrant activation of the DDR/apoptotic pathway in breast tumorigenesis. Negative correlation was observed between methylation status and transcript expression levels for TRAIL, DR4, CASP8, ATM, CHEK2, BRCA1 and BRCA2 CpG sites. Categorization of the gene methylation with respect to the clinicopathological parameters showed an increase in aberrant methylation pattern in advanced tumors. These uncharacteristic methylation patterns corresponded with decreased death receptor apoptosis (P=0.047) and DNA damage repair potential (P=0.004) in advanced tumors. The observation of BRCA2 -26 G/A 5'UTR polymorphism concomitant with the presence of methylation in the promoter region was novel and emerged as a strong candidate for susceptibility to sporadic breast tumors. CONCLUSION: Our study indicates that methylation of DDR-apoptotic gene promoters in sporadic breast cancer is not a random phenomenon. Progressive epigenetic alterations in advancing tumors result in aberrant DDR-apoptotic pathway thereby promoting tumor development. We propose, since pathological epigenetic changes of the DDR-apoptotic genes are reversible modifications, these could further be targeted for therapeutic interventions.

The Indian origin of paternal haplogroup R1a1* substantiates the autochthonous origin of Brahmins and the caste system.

Many major rival models of the origin of the Hindu caste system co-exist despite extensive studies, each with associated genetic evidences. One of the major factors that has still kept the origin of the Indian caste system obscure is the unresolved question of the origin of Y-haplogroup R1a1*, at times associated with a male-mediated major genetic influx from Central Asia or Eurasia, which has contributed to the higher castes in India. Y-haplogroup R1a1* has a widespread distribution and high frequency across Eurasia, Central Asia and the Indian subcontinent, with scanty reports of its ancestral (R*, R1* and R1a*) and derived lineages (R1a1a, R1a1b and R1a1c). To resolve these issues, we screened 621 Y-chromosomes (of Brahmins occupying the upper-most caste position and schedule castes/tribals occupying the lower-most positions) with 55 Y-chromosomal binary markers and seven Y-microsatellite markers and compiled an extensive dataset of 2809 Y-chromosomes (681 Brahmins, and 2128 tribals and schedule castes) for conclusions. A peculiar observation of the highest frequency (up to 72.22%) of Y-haplogroup R1a1* in Brahmins hinted at its presence as a founder lineage for this caste group. Further, observation of R1a1* in different tribal population groups, existence of Y-haplogroup R1a* in ancestors and extended phylogenetic analyses of the pooled dataset of 530 Indians, 224 Pakistanis and 276 Central Asians and Eurasians bearing the R1a1* haplogroup supported the autochthonous origin of R1a1 lineage in India and a tribal link to Indian Brahmins. However, it is important to discover novel Y-chromosomal binary marker(s) for a higher resolution of R1a1* and confirm the present conclusions.