Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa

Purpose: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. Methods: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. Results: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. Conclusions: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Transcription of Human Resistin Gene Involves an Interaction of Sp1 with Peroxisome Proliferator-Activating Receptor Gamma (PPAR gamma)

Background: Resistin is a cysteine rich protein, mainly expressed and secreted by circulating human mononuclear cells. While several factors responsible for transcription of mouse resistin gene have been identified, not much is known about the factors responsible for the differential expression of human resistin.Methodology/Principal Finding: We show that the minimal promoter of human resistin lies within similar to 80 bp sequence upstream of the transcriptional start site (-240) whereas binding sites for cRel, CCAAT enhancer binding protein alpha (C/EBP-alpha), activating transcription factor 2 (ATF-2) and activator protein 1 (AP-1) transcription factors, important for induced expression, are present within sequences up to -619. Specificity Protein 1(Sp1) binding site (-276 to -295) is also present and an interaction of Sp1 with peroxisome proliferator activating receptor gamma (PPAR gamma) is necessary for constitutive expression in U937 cells. Indeed co-immunoprecipitation assay demonstrated a direct physical interaction of Sp1 with PPAR gamma in whole cell extracts of U937 cells. Phorbol myristate acetate (PMA) upregulated the expression of resistin mRNA in U937 cells by increasing the recruitment of Sp1, ATF-2 and PPAR gamma on the resistin gene promoter. Furthermore, PMA stimulation of U937 cells resulted in the disruption of Sp1 and PPAR gamma interaction. Chromatin immunoprecipitation (ChIP) assay confirmed the recruitment of transcription factors phospho ATF-2, Sp1, Sp3, PPAR gamma, chromatin modifier histone deacetylase 1 (HDAC1) and the acetylated form of histone H3 but not cRel, C/EBP-alpha and phospho c-Jun during resistingene transcription.Conclusion: Our findings suggest a complex interplay of Sp1 and PPAR gamma along with other transcription factors that drives the expression of resistin in human monocytic U937 cells.

Stochastic simulations of the origins and implications of long-tailed distributions in gene expression

Gene expression noise results in protein number distributions ranging from long-tailed to Gaussian. We show how long-tailed distributions arise from a stochastic model of the constituent chemical reactions and suggest that, in conjunction with cooperative switches, they lead to more sensitive selection of a subpopulation of cells with high protein number than is possible with Gaussian distributions. Single-cell-tracking experiments are presented to validate some of the assumptions of the stochastic simulations. We also examine the effect of DNA looping on the shape of protein distributions. We further show that when switches are incorporated in the regulation of a gene via a feedback loop, the distributions can become bimodal. This might explain the bimodal distribution of certain morphogens during early embryogenesis.

Characterization of the human SLC22A18 gene promoter and its regulation by the transcription factor Sp1

SLC22A18, a poly-specific organic cation transporter, is paternally imprinted in humans and mice. It shows loss-of-heterozygosity in childhood and adult tumors, and gain-of-imprinting in hepatocarcinomas and breast cancers. Despite the importance of this gene, its transcriptional regulation has not been studied, and the promoter has not yet been characterized. We therefore set out to identify the potential cis-regulatory elements including the promoter of this gene. The luciferase reporter assay in human cells indicated that a region from -120 by to +78 by is required for the core promoter activity. No consensus TATA or CHAT boxes were found in this region, but two Sp1 binding sites were conserved in human, chimpanzee, mouse and rat. Mutational analysis of the two Sp1 sites suggested their requirement for the promoter activity. Chromatin-immunoprecipitation showed binding of Sp1 to the promoter region in vivo. Overexpression of Sp1 in Drosophila Sp1-null SL2 cells suggested that Sp1 is the transactivator of the promoter. The human core promoter was functional in mouse 3T3 and monkey COS7 cells. We found a CpG island which spanned the core promoter and exon 1. COBRA technique did not reveal promoter methylation in 10 normal oral tissues, 14 oral tumors, and two human cell lines HuH7 and A549. This study provides the first insight into the mechanism that controls expression of this imprinted tumor suppressor gene. A COBRA-based assay has been developed to look for promoter methylation in different cancers. The present data will help to understand the regulation of this gene and its role in tumorigenesis. (C) 2008 Elsevier B.V. All rights reserved.

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa

Purpose: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. Methods: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. Results: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease.

Characterization of the human SLC22A18 gene promoter and its regulation by the transcription factor Sp1

SLC22A18, a poly-specific organic cation transporter, is paternally imprinted in humans and mice. It shows loss-of-heterozygosity in childhood and adult tumors, and gain-of-imprinting in hepatocarcinomas and breast cancers. Despite the importance of this gene, its transcriptional regulation has not been studied, and the promoter has not yet been characterized. We therefore set out to identify the potential cis-regulatory elements including the promoter of this gene. The luciferase reporter assay in human cells indicated that a region from -120 by to +78 by is required for the core promoter activity. No consensus TATA or CHAT boxes were found in this region, but two Sp1 binding sites were conserved in human, chimpanzee, mouse and rat. Mutational analysis of the two Sp1 sites suggested their requirement for the promoter activity. Chromatin-immunoprecipitation showed binding of Sp1 to the promoter region in vivo. Overexpression of Sp1 in Drosophila Sp1-null SL2 cells suggested that Sp1 is the transactivator of the promoter. The human core promoter was functional in mouse 3T3 and monkey COS7 cells. We found a CpG island which spanned the core promoter and exon 1. COBRA technique did not reveal promoter methylation in 10 normal oral tissues, 14 oral tumors, and two human cell lines HuH7 and A549. This study provides the first insight into the mechanism that controls expression of this imprinted tumor suppressor gene. A COBRA-based assay has been developed to look for promoter methylation in different cancers. The present data will help to understand the regulation of this gene and its role in tumorigenesis. (C) 2008 Elsevier B.V. All rights reserved.

Array-based gene expression, CGH and tissue data defines a 12q24 gain in neuroblastic tumors with prognostic implication

Neuroblastoma has successfully served as a model system for the identification of neuroectoderm-derived oncogenes. However, in spite of various efforts, only a few clinically useful prognostic markers have been found. Here, we present a framework, which integrates DNA, RNA and tissue data to identify and prioritize genetic events that represent clinically relevant new therapeutic targets and prognostic biomarkers for neuroblastoma.

Part I: Enzyme replacement versus neurotrophic factor viral vector-mediated gene therapy of Parkinson’s disease, Part II: The effects of orally dosed tolcapone and entacapone on dopamine metabolism in the dorsal striatum and nucleus accumbens in rats

Part I: Parkinson’s disease is a slowly progressive neurodegenerative disorder in which particularly the dopaminergic neurons of the substantia nigra pars compacta degenerate and die. Current conventional treatment is based on restraining symptoms but it has no effect on the progression of the disease. Gene therapy research has focused on the possibility of restoring the lost brain function by at least two means: substitution of critical enzymes needed for the synthesis of dopamine and slowing down the progression of the disease by supporting the functions of the remaining nigral dopaminergic neurons by neurotrophic factors. The striatal levels of enzymes such as tyrosine hydroxylase, dopadecarboxylase and GTP-CH1 are decreased as the disease progresses. By replacing one or all of the enzymes, dopamine levels in the striatum may be restored to normal and behavioral impairments caused by the disease may be ameliorated especially in the later stages of the disease. The neurotrophic factors glial cell derived neurotrophic factor (GDNF) and neurturin have shown to protect and restore functions of dopaminergic cell somas and terminals as well as improve behavior in animal lesion models. This therapy may be best suited at the early stages of the disease when there are more dopaminergic neurons for neurotrophic factors to reach. Viral vector-mediated gene transfer provides a tool to deliver proteins with complex structures into specific brain locations and provides long-term protein over-expression. Part II: The aim of our study was to investigate the effects of two orally dosed COMT inhibitors entacapone (10 and 30 mg/kg) and tolcapone (10 and 30 mg/kg) with a subsequent administration of a peripheral dopadecarboxylase inhibitor carbidopa (30 mg/kg) and L- dopa (30 mg/kg) on dopamine and its metabolite levels in the dorsal striatum and nucleus accumbens of freely moving rats using dual-probe in vivo microdialysis. Earlier similarly designed studies have only been conducted in the dorsal striatum. We also confirmed the result of earlier ex vivo studies regarding the effects of intraperitoneally dosed tolcapone (30 mg/kg) and entacapone (30 mg/kg) on striatal and hepatic COMT activity. The results obtained from the dorsal striatum were generally in line with earlier studies, where tolcapone tended to increase dopamine and DOPAC levels and decrease HVA levels. Entacapone tended to keep striatal dopamine and HVA levels elevated longer than in controls and also tended to elevate the levels of DOPAC. Surprisingly in the nucleus accumbens, dopamine levels after either dose of entacapone or tolcapone were not elevated. Accumbal DOPAC levels, especially in the tolcapone 30 mg/kg group, were elevated nearly to the same extent as measured in the dorsal striatum. Entacapone 10 mg/kg elevated accumbal HVA levels more than the dose of 30 mg/kg and the effect was more pronounced in the nucleus accumbens than in the dorsal striatum. This suggests that entacapone 30 mg/kg has minor central effects. Also our ex vivo study results obtained from the dorsal striatum suggest that entacapone 30 mg/kg has minor and transient central effects, even though central HVA levels were not suppressed below those of the control group in either brain area in the microdialysis study. Both entacapone and tolcapone suppressed hepatic COMT activity more than striatal COMT activity. Tolcapone was more effective than entacapone in the dorsal striatum. The differences between dopamine and its metabolite levels in the dorsal striatum and nucleus accumbens may be due to different properties of the two brain areas.

Genetic variation in the LDL receptor-related protein 5 (LRP5) gene : Association with bone health and metabolic parameters

Bone mass accrual and maintenance are regulated by a complex interplay between genetic and environmental factors. Recent studies have revealed an important role for the low-density lipoprotein receptor-related protein 5 (LRP5) in this process. The aim of this thesis study was to identify novel variants in the LRP5 gene and to further elucidate the association of LRP5 and its variants with various bone health related clinical characteristics. The results of our studies show that loss-of-function mutations in LRP5 cause severe osteoporosis not only in homozygous subjects but also in the carriers of these mutations, who have significantly reduced bone mineral density (BMD) and increased susceptibility to fractures. In addition, we demonstrated for the first time that a common polymorphic LRP5 variant (p.A1330V) was associated with reduced peak bone mass, an important determinant of BMD and osteoporosis in later life. The results from these two studies are concordant with results seen in other studies on LRP5 mutations and in association studies linking genetic variation in LRP5 with BMD and osteoporosis. Several rare LRP5 variants were identified in children with recurrent fractures. Sequencing and multiplex ligation-dependent probe amplification (MLPA) analyses revealed no disease-causing mutations or whole-exon deletions. Our findings from clinical assessments and family-based genotype-phenotype studies suggested that the rare LRP5 variants identified are not the definite cause of fractures in these children. Clinical assessments of our study subjects with LPR5 mutations revealed an unexpectedly high prevalence of impaired glucose tolerance and dyslipidaemia. Moreover, in subsequent studies we discovered that common polymorphic LRP5 variants are associated with unfavorable metabolic characteristics. Changes in lipid profile were already apparent in pre-pubertal children. These results, together with the findings from other studies, suggest an important role for LRP5 also in glucose and lipid metabolism. Our results underscore the important role of LRP5 not only in bone mass accrual and maintenance of skeletal health but also in glucose and lipid metabolism. The role of LRP5 in bone metabolism has long been studied, but further studies with larger study cohorts are still needed to evaluate the specific role of LRP5 variants as metabolic risk factors.

Cytosines, but Not Purines, Determine Recombination Activating Gene (RAG)-induced Breaks on Heteroduplex DNA Structures

The sequence specificity of the recombination activating gene (RAG) complex during V(D)J recombination has been well studied. RAGs can also act as structure-specific nuclease; however, little is known about the mechanism of its action. Here, we show that in addition to DNA structure, sequence dictates the pattern and efficiency of RAG cleavage on altered DNA structures. Cytosine nucleotides are preferentially nicked by RAGs when present at single-stranded regions of heteroduplex DNA. Although unpaired thymine nucleotides are also nicked, the efficiency is many fold weaker. Induction of single- or double-strand breaks by RAGs depends on the position of cytosines and whether it is present on one or both of the strands. Interestingly, RAGs are unable to induce breaks when adenine or guanine nucleotides are present at single-strand regions. The nucleotide present immediately next to the bubble sequence could also affect RAG cleavage. Hence, we propose “C(d)C(S)C(S)” (d, double-stranded; s, single-stranded) as a consensus sequence for RAG-induced breaks at single-/double-strand DNA transitions. Such a consensus sequence motif is useful for explaining RAG cleavage on other types of DNA structures described in the literature. Therefore, the mechanism of RAG cleavage described here could explain facets of chromosomal rearrangements specific to lymphoid tissues leading to genomic instability.

Association of methylenetetrahydrofolate reductase gene polymorphisms & colorectal cancer in India

Background & objectives: Methylenetetrahydrofolate reductase (MTHFR) is a critical enzyme in folate metabolism and involved in DNA synthesis, DNA repair and DNA methylation. The two common functional polymorphisms of MTHFR, 677C -> T and 1298 A -> C have shown to impact several diseases including cancer. This case-control study was undertaken to analyse the association of the MTHFR gene polymorphisms 677 C -> T and 1298 A -> C and risk of colorectal cancer (CRC).Methods: One hundred patients with a confirmed histopathologic diagnosis of CRC and 86 age and gender matched controls with no history of cancer were taken for this study. DNA was isolated from peripheral blood samples and the genotypes were determined by PCR-RFLP. The risk association was estimated by compounding odds ratio (OR) with 95 per cent confidence interval (CI). Results: Genotype frequency of MTHFR 677 CC, CT and TT were 76.7, 22.1 and 1.16 per cent in controls, and 74,25 and 1.0 per cent among patients. The 'T' allele frequency was 12.21 and 13.5 per cent in controls and patients respectively. The genotype frequency of MTHFR 1298 AA, AC, and CC were 25.6, 58.1 and 16.3 per cent for controls and 22, 70 and 8 per cent for patents respectively. The 'C' allele frequency for 1298 A -> C was 43.0 and 45.3 per cent respectively for controls and patients. The OR for 677 CT was 1.18 (95% CI 0.59-2.32, P = 0.642), OR for 1298 AC was 1.68 (95% CI 0.92-3.08, P = 0.092) and OR for 1298 CC was 0.45(95% CI 0.18-1.12, P = 0.081). The OR for the combined heterozygous state (677 CT and 1298 AC) was 1.18(95% CI 0.52-2.64, P =0.697).Interpretation & conclusion: The frequency of the MTHFR 677 TT genotype is rare as compared to 1298 CC genotype in the population studied. There was no association between 677 C -> T and 1298 A -> C polymorphisms and risk of CRC either individually or in combination. The homozygous state for 1298 A -> C polymorphism appears to slightly lower risk of CRC. This needs to be confirmed with a larger sample size.

Drosophila indirect flight muscle specific Act88F actin mutants as a model system for studying congenital myopathies of the human ACTA1 skeletal muscle actin gene

Most human ACTA1 skeletal actin gene mutations cause dominant, congenital myopathies often with severely reduced muscle function and neonatal mortality. High sequence conservation of actin means many mutated ACTA1 residues are identical to those in the Drosophila Act88F, an indirect flight muscle specific sarcomeric actin. Four known Act88F mutations occur at the same actin residues mutated in ten ACTA1 nemaline mutations, A138D/P, R256H/L, G268C/D/R/S and R372C/S. These Act88F mutants were examined for similar muscle phenotypes. Mutant homozygotes show phenotypes ranging from a lack of myofibrils to almost normal sarcomeres at eclosion. Aberrant Z-disc-like structures and serial Z-disc arrays, ‘zebra bodies’, are observed in homozygotes and heterozygotes of all four Act88F mutants. These electron-dense structures show homologies to human nemaline bodies/rods, but are much smaller than those typically found in the human myopathy. We conclude that the Drosophila indirect flight muscles provide a good model system for studying ACTA1 mutations.

Methylenetetrahydrofolate reductase gene polymorphisms and risk of acute lymphoblastic leukemia in children

Methylenetetrahydrofolate reductase (MTHFR) is a critical enzyme in folate metabolism and is involved in DNA synthesis, DNA repair and DNA methylation. Genetic polymorphisms of this enzyme have been shown to impact several diseases, including cancer. Leukemias are malignancies arising from rapidly proliferating hematopoietic cells having great requirement of DNA synthesis. This case-control study was undertaken to analyze the association of the MTHFR gene polymorphisms 677 C"T and 1298 A"C and the risk of acute lymphoblastic leukemia in children. Materials and Methods: Eighty-six patients aged below 15 years with a confirmed diagnosis of acute lymphoblastic leukemia (ALL) and 99 matched controls were taken for this study. Analysis of the polymorphisms was done using the polymerase chain reaction -restriction fragment length polymorphism (PCR-RFLP) method. Results: Frequency of MTHFR 677 CC and CT were 85.9% and 14.1% in the controls, and 84.9% and 15.1% in the cases. The 'T' allele frequency was 7% and 7.5% in cases and controls respectively. The frequency of MTHFR 1298 AA, AC, and CC were 28.3%, 55.6% and 16.1% for controls and 23.3%, 59.3% and 17.4% for cases respectively. The 'C' allele frequency for 1298 A→C was 43.9% and 47% respectively for controls and cases. The odds ratio (OR) for C677T was 1.08 (95% CI 0.48- 2.45, p = 0.851) and OR for A1298C was 1.29(95% CI 0.65-2.29, p = 0.46) and OR for 1298 CC was 1.31 (95% CI 0.53-3.26, p =0.56). The OR for the combined heterozygous status (677 CT and 1298 AC) was 1.94 (95% CI 0.58 -6.52, p = 0.286). Conclusion: The prevalence of 'T' allele for 677 MTHFR polymorphism was low in the population studied. There was no association between MTHFR 677 C→T and 1298 A→C gene polymorphisms and risk of ALL, which may be due to the small sample size.

Regulation of cytochrome P-450b/e gene expression by a heme- and phenobarbitone-modulated transcription factor

The cloned DNA fragment of the cytochrome P-450b/e gene containing the upstream region from position -179 through part of the first exon is faithfully transcribed in freeze-thawed rat liver nuclei. Phenobarbitone treatment of the animal strikingly increases this transcription, and the increase is blocked by cycloheximide (protein synthesis inhibitor) or CoCl2 (heme biosynthetic inhibitor) treatment of animals. This picture correlates very well with the reported cytochrome P-450b/e mRNA levels in vivo and run-on transcription rates in vitro under these conditions. The upstream region (from position -179) was assessed for protein binding with nuclear extracts by nitrocellulose filter binding, gel retardation, DNase I treatment ("footprinting"), and Western blot analysis. Phenobarbitone treatment dramatically increases protein binding to the upstream region, an increase once again blocked by cycloheximide or CoCl2 treatments. Addition of heme in vitro to heme-deficient nuclei and nuclear extracts restores the induced levels of transcription and protein binding to the upstream fragment, respectively. Thus, drug-mediated synthesis and heme-modulated binding of a transcription factor(s) appear involved in the transcriptional activation of the cytochrome P-450b/e genes, and an 85-kDa protein may be a major factor in this regard.

Gene expression profile of epithelial cells and mesenchymal cells derived from limbal explant culture

Purpose: Limbal stem cell deficiency is a challenging clinical problem and the current treatment involves replenishing the depleted limbal stem cell (LSC) pool by either limbal tissue transplantation or use of cultivated limbal epithelial cells (LEC). Our experience of cultivating the LEC on denuded human amniotic membrane using a feeder cell free method, led to identification of mesenchymal cells of limbus (MC-L), which showed phenotypic resemblance to bone marrow derived mesenchymal stem cells (MSC-BM). To understand the transcriptional profile of these cells, microarray experiments were carried out.Methods: RNA was isolated from cultured LEC, MC-L and MSC-BM and microarray experiments were carried out by using Agilent chip (4x44 k). The microarray data was validated by using Realtime and semiquntitative reverse transcription polymerase chain reaction. Results: The microarray analysis revealed specific gene signature of LEC and MC-L, and also their complementary role related to cytokine and growth factor profile, thus supporting the nurturing roles of the MC-L. We have also observed similar and differential gene expression between MC-L and MSC-BM.Conclusions: This study represents the first extensive gene expression analysis of limbal explant culture derived epithelial and mesenchymal cells and as such reveals new insight into the biology, ontogeny, and in vivo function of these cells.