Cystathionine beta-synthase variants : identification, characterization and modulation by thiol compounds

Tese de doutoramento, Farmácia (Bioquímica), Universidade de Lisboa, Faculdade de Farmácia, 2014

Mutagénèse semi-aléatoire au site actif de la DHFR humaine : création et caractérisation de variantes hautement résistantes au MTX.

La dihydrofolate réductase humaine (DHFRh) est une enzyme essentielle à la prolifération cellulaire. Elle réduit le dihydrofolate en tétrahydrofolate, un co-facteur impliqué dans la biosynthèse des purines et du thymidylate. La DHFRh est une cible de choix pour des agents de chimiothérapie comme le méthotrexate (MTX), inhibant spécifiquement l’enzyme ce qui mène à un arrêt de la prolifération et ultimement à la mort cellulaire. Le MTX est utilisé pour le traitement de plusieurs maladies prolifératives, incluant le cancer. La grande utilisation du MTX dans le milieu clinique a mené au développement de mécanismes de résistance, qui réduisent l’efficacité de traitement. La présente étude se penche sur l’un des mécanismes de résistance, soit des mutations dans la DHFRh qui réduisent son affinité pour le MTX, dans le but de mieux comprendre les éléments moléculaires requis pour la reconnaissance de l’inhibiteur au site actif de l’enzyme. En parallèle, nous visons à identifier des variantes plus résistantes au MTX pour leur utilisation en tant que marqueurs de sélection en culture cellulaire pour des systèmes particuliers, tel que la culture de cellules hématopoïétiques souches (CHS), qui offrent des possibilités intéressantes dans le domaine de la thérapie cellulaire. Pour étudier le rôle des différentes régions du site actif, et pour vérifier la présence d’une corrélation entre des mutations à ces régions et une augmentation de la résistance au MTX, une stratégie combinatoire a été dévelopée pour la création de plusieurs banques de variantes à des résidus du site actif à proximité du MTX lié. Les banques ont été sélectionnées in vivo dans un système bactérien en utilisant des milieux de croissance contenant des hautes concentrations de MTX. La banque DHFRh 31/34/35 généra un nombre considérable de variantes combinatoires de la DHFRh hautement résistantes au MTX. Les variantes les plus intéressantes ont été testées pour leur potentiel en tant que marqueur de sélection dans plusieurs lignées cellulaires, dont les cellules hématopoïétiques transduites. Une protection complète contre les effets cytotoxiques du MTX a été observée chez ces cellules suite à leur infection avec les variantes combinatoires. Pour mieux comprendre les causes moléculaires reliées à la résistance au MTX, des études de structure tridimensionnelle de variantes liées au MTX ont été entreprises. La résolution de la structure de la double variante F31R/Q35E lié au MTX a révélé que le phénotype de résistance était attribuable à d’importantes différences entre le site actif de la double variante et de l’enzyme native, possiblement dû à un phénomème dynamique. Une compréhension plus générale de la reconnaissance et la résistance aux antifolates a été réalisée en comparant des séquences et des structures de variantes de la DHFR résistants aux antifolates et provenant de différentes espèces. En somme, ces travaux apportent de nouveaux éléments pour la comprehension des intéractions importantes entre une enzyme et un ligand, pouvant aider au développement de nouveaux antifolates plus efficaces pour le traitement de diverses maladies. De plus, ces travaux ont généré de nouveaux gènes de résistance pouvant être utilisés en tant que marqueurs de sélection en biologie cellulaire.

Proportions of A1, A2, B and C β-casein protein variants in retail milk in the UK

The A1 variant protein of the β-casein family has been implicated in various disease states although much evidence is weak or contradictory. The primary objective was to measure, for the first time, the proportions of the key β-casein variant proteins in UK retail milk over the course of one year. In total, 55 samples of semi-skimmed milk were purchased from five supermarkets over the course of a year and the proportions of the A1, A2, B and C casein variant proteins were measured, using high resolution HPLC-MS. The results showed that β-casein in UK retail milk comprises approximately 0.58, 0.31, 0.07 and 0.03 A2, A1, B and C protein variants, respectively. The proportion of A2 is higher than some early studies would predict although the reasons for this and any implications for health are unclear

Protein variation in blood-dwelling schistosome worms generated by differential splicing of micro-exon gene transcripts

Schistosoma mansoni is a well-adapted blood-dwelling parasitic helminth, persisting for decades in its human host despite being continually exposed to potential immune attack. Here, we describe in detail micro-exon genes (MEG) in S. mansoni, some present in multiple copies, which represent a novel molecular system for creating protein variation through the alternate splicing of short (<= 36 bp) symmetric exons organized in tandem. Analysis of three closely related copies of one MEG family allowed us to trace several evolutionary events and propose a mechanism for micro-exon generation and diversification. Microarray experiments show that the majority of MEGs are up-regulated in life cycle stages associated with establishment in the mammalian host after skin penetration. Sequencing of RT-PCR products allowed the description of several alternate splice forms of micro-exon genes, highlighting the potential use of these transcripts to generate a complex pool of protein variants. We obtained direct evidence for the existence of such pools by proteomic analysis of secretions from migrating schistosomula and mature eggs. Whole-mount in situ hybridization and immunolocalization showed that MEG transcripts and proteins were restricted to glands or epithelia exposed to the external environment. The ability of schistosomes to produce a complex pool of variant proteins aligns them with the other major groups of blood parasites, but using a completely different mechanism. We believe that our data open a new chapter in the study of immune evasion by schistosomes, and their ability to generate variant proteins could represent a significant obstacle to vaccine development.

Dynamic Activation and Repression of the Plasmodium falciparum rif Gene Family and Their Relation to Chromatin Modification

The regulation of variant gene expression in Plasmodium falciparum is still only partially understood. Regulation of var genes, the most studied gene family involved in antigenic variation, is orchestrated by a dynamic pattern of inherited chromatin states. Although recent evidence pointed to epigenetic regulation of transcribed and repressed rif loci, little is known about specific on/off associated histone modifications of individual rif genes. To investigate the chromatin marks for transcribed and repressed rif loci, we cultivated parasites and evaluated the transcriptional status of chosen rif targets by qRT-PCR and performed ChIP assays using H3K9ac and H3K9me3 antibodies. We then monitored changes in the epigenetic patterns in parasites after several reinvasions and also evaluated the "poised'' mark in trophozoites and schizonts of the same erythrocytic cycle by ChIP using H3K4me2 specific antibodies. Our results show that H3K9 is acetylated in transcribed rif loci and trimethylated or even unmodified in repressed rif loci. These transcriptional and epigenetic states are inherited after several reinvasions. The poised modification H3K4me2 showed a tendency to be more present in loci in trophozoites that upon progression to schizonts strongly transcribe the respective locus. However, this effect was not consistently observed for all monitored loci. While our data show important similarities to var transcription-associated chromatin modifications, the observed swiftly occurring modifications at rif loci and the absence of H3K9 modification point to a different dynamic of recruitment of chromatin modifying enzymes.

A Z-DNA binding domain present in the human editing enzyme, double-stranded RNA adenosine deaminase

Editing of RNA changes the read-out of information from DNA by altering the nucleotide sequence of a transcript. One type of RNA editing found in all metazoans uses double-stranded RNA (dsRNA) as a substrate and results in the deamination of adenosine to give inosine, which is translated as guanosine. Editing thus allows variant proteins to be produced from a single pre-mRNA. A mechanism by which dsRNA substrates form is through pairing of intronic and exonic sequences before the removal of noncoding sequences by splicing. Here we report that the RNA editing enzyme, human dsRNA adenosine deaminase (DRADA1, or ADAR1) contains a domain (Zα) that binds specifically to the left-handed Z-DNA conformation with high affinity (KD = 4 nM). As formation of Z-DNA in vivo occurs 5′ to, or behind, a moving RNA polymerase during transcription, recognition of Z-DNA by DRADA1 provides a plausible mechanism by which DRADA1 can be targeted to a nascent RNA so that editing occurs before splicing. Analysis of sequences related to Zα has allowed identification of motifs common to this class of nucleic acid binding domain.

Use of a novel DNA melting profile assay for the identification of PCR-amplified genomic sequences encoding for variant-specific surface proteins from the clonal GS/M-83-H7 line of Giardia lamblia.

During infections, Giardia lamblia undergoes a continuous change of its major surface antigens, the variant-specific surface proteins (VSPs). Many studies on antigenic variation have been performed using G. lamblia clone GS/M-83-H7, which expresses surface antigen VSP H7. The present study was focused on the identification and characterization of vsp gene sequences within the genome of the clonal G. lamblia GS/M-83-H7 line. For this purpose, we applied a PCR which specifically amplified truncated sequences from the 3'-terminal region of the vsp genes. Upon cloning, most of the vsp gene amplification products were shown to be approximately identical in size and thus could not be distinguished from each other by conventional gel electrophoresis. In order to pre-estimate the sequence complexity within the large panel of vsp clones isolated, we elaborated a novel concept which facilitated our large-scale genetic screening approach: PCR products from cloned DNA molecules were generated and then subjected to a DNA melting profile assay based on the use of the LightCycler Instrument. This high-throughput assay system proved to be well suited to monitor sequence differences between the amplification products from closely related vsp genes and thus could be used for the primary, sequence-related discrimination of the corresponding clones. After testing 50 candidates, vsp clones could be divided into five groups, each characterized by an individual DNA melting profile of the corresponding amplification products. Sequence analysis of some of these 50 candidates confirmed data from the aforementioned assay in that clones were demonstrated to be identical within, but different between, the distinct groups. The nucleotide and deduced amino acid sequences of five representative vsp clones showed high similarities both among each other and also with the corresponding gene segment of the variant-specific surface antigen (VSP H7) expressed by the original GS/M-83-H7 variant type. Furthermore, three of the genomic vsp sequences turned out to be identical to vsp sequences that represented previously characterized transcription products from in vivo- or in vitro-switched GS/M-83-H7 trophozoites. In conclusion, the DNA melting profile assay seems to be a versatile tool for the PCR-based genotyping of moderately or highly diversified sequence orthologues.

G-protein β3 subunit gene (GNB3) variant in causation of essential hypertension

Essential hypertensives display enhanced signal transduction through pertussis toxin-sensitive G proteins. The T allele of a C825T variant in exon 10 of the G protein β3 subunit gene (GNB3) induces formation of a splice variant (Gβ3-s) with enhanced activity. The T allele of GNB3 was shown recently to be associated with hypertension in unselected German patients (frequency=0.31 versus 0.25 in control). To confirm and extend this finding in a different setting, we performed an association study in Australian white hypertensives. This involved an extensively examined cohort of 110 hypertensives, each of whom were the offspring of 2 hypertensive parents, and 189 normotensives whose parents were both normotensive beyond age 50 years. Genotyping was performed by polymerase chain reaction and digestion with BseDI, which either cut (C allele) or did not cut (T allele) the 268-bp polymerase chain reaction product. T allele frequency in the hypertensive group was 0.43 compared with 0.25 in the normotensive group (χ2=22; P=0.00002; odds ratio=2.3; 95% CI=1.7 to 3.3). The T allele tracked with higher pretreatment blood pressure: diastolic=105±7, 109±16, and 128±28 mm Hg (mean±SD) for CC, CT, and 7T, respectively (P=0.001 by 1-way ANOVA). Blood pressures were higher in female hypertensives with a T allele (P=0.006 for systolic and 0.0003 for diastolic by ANOVA) than they were in male hypertensives. In conclusion, the present study of a group with strong family history supports a role for a genetically determined, physiologically active splice variant of the G protein β3 subunit gene in the causation of essential hypertension.

A variant of the breast cancer type 2 susceptibility protein (BRC) repeat is essential for the RECQL5 Helicase to interact with RAD51 Recombinase for genome stabilization

The BRC repeat is a structural motif in the tumor suppressor BRCA2 (breast cancer type 2 susceptibility protein), which promotes homologous recombination (HR) by regulating RAD51 recombinase activity. To date, the BRC repeat has not been observed in other proteins, so that its role in HR is inferred only in the context of BRCA2. Here, we identified a BRC repeat variant, named BRCv, in the RECQL5 helicase, which possesses anti-recombinase activity in vitro and suppresses HR and promotes cellular resistance to camptothecin-induced replication stress in vivo. RECQL5-BRCv interacted with RAD51 through two conserved motifs similar to those in the BRCA2-BRC repeat. Mutations of either motif compromised functions of RECQL5, including association with RAD51, inhibition of RAD51-mediated D-loop formation, suppression of sister chromatid exchange, and resistance to camptothecin-induced replication stress. Potential BRCvs were also found in other HR regulatory proteins, including Srs2 and Sgs1, which possess anti-recombinase activities similar to that of RECQL5. A point mutation in the predicted Srs2-BRCv disrupted the ability of the protein to bind RAD51 and to inhibit D-loop formation. Thus, BRC is a common RAD51 interaction module that can be utilized by different proteins to either promote HR, as in the case of BRCA2, or to suppress HR, as in RECQL5.

Correlation between BRAF mutation and the clinicopathological parameters in papillary thyroid carcinoma with particular reference to follicular variant

Mutation of the BRAF gene is common in thyroid cancer. Follicular variant of papillary thyroid carcinoma is a variant of papillary thyroid carcinoma that has created continuous diagnostic controversies among pathologists. The aims of this study are to (1) investigate whether follicular variant of papillary thyroid carcinoma has a different pattern of BRAF mutation than conventional papillary thyroid carcinoma in a large cohort of patients with typical features of follicular variant of papillary thyroid carcinoma and (2) to study the relationship of clinicopathological features of papillary thyroid carcinomas with BRAF mutation. Tissue blocks from 76 patients with diagnostic features of papillary thyroid carcinomas (40 with conventional type and 36 with follicular variant) were included in the study. From these, DNA was extracted and BRAF V600E mutations were detected by polymerase chain reaction followed by restriction enzyme digestion and sequencing of exon 15. Analysis of the data indicated that BRAF V600E mutation is significantly more common in conventional papillary thyroid carcinoma (58% versus 31%, P = .022). Furthermore, the mutation was often noted in female patients (P = .017), in high-stage cancers (P = .034), and in tumors with mild lymphocytic thyroiditis (P = .006). We concluded that follicular variant of papillary thyroid carcinoma differs from conventional papillary thyroid carcinoma in the rate of BRAF mutation. The results of this study add further information indicating that mutations in BRAF play a role in thyroid cancer development and progression.

Follicular variant of papillary thyroid carcinoma : a diagnostic challenge for clinicians and pathologists

The follicular variant of papillary thyroid carcinoma (FVPTC) presents a type of papillary thyroid cancer that has created continuous diagnosis and treatment controversies among clinicians and pathologists. In this review, we describe the nomenclature, the clinical features, diagnostic problems and the molecular biology of FVPTC. It is important for clinicians to understand this entity as the diagnosis and management of this group of patient may be different from other patients with conventional PTC. The literature suggests that FVPTC behaves in a way similar, clinically, to conventional papillary thyroid carcinoma. However, there are some genotypic differences which may characterise this neoplasm. These parameters may account for the phenotypic variation described by some scientists in this type of cancer. Further understanding can only be achieved by defining strict pathological criteria, in-depth study of the molecular biology and long term follow-up of the optional patients with FVPTC.

Genome-wide association study identifies a variant in HDAC9 associated with large vessel ischemic stroke

Genetic factors have been implicated in stroke risk, but few replicated associations have been reported. We conducted a genome-wide association study (GWAS) for ischemic stroke and its subtypes in 3,548 affected individuals and 5,972 controls, all of European ancestry. Replication of potential signals was performed in 5,859 affected individuals and 6,281 controls. We replicated previous associations for cardioembolic stroke near PITX2 and ZFHX3 and for large vessel stroke at a 9p21 locus. We identified a new association for large vessel stroke within HDAC9 (encoding histone deacetylase 9) on chromosome 7p21.1 (including further replication in an additional 735 affected individuals and 28,583 controls) (rs11984041; combined P = 1.87 × 10 -11; odds ratio (OR) = 1.42, 95% confidence interval (CI) = 1.28-1.57). All four loci exhibited evidence for heterogeneity of effect across the stroke subtypes, with some and possibly all affecting risk for only one subtype. This suggests distinct genetic architectures for different stroke subtypes.

Localization of testis-variant histones in rat testis chromatin

Nucleosome core particles and oligonucleosomes were isolated by digesting rat testis nuclei with micrococcal nuclease to 20% acid-solubility, followed by fractionation of the digest on a Bio-Gel A-5m column. The core particles thus isolated were characterized on the basis of their DNA length of 151 +/- 5 base-pairs and sedimentation coefficient of 11.4S. Analysis of the acid-soluble proteins of the core particles indicated that histones TH2B and X2 are constituents of the core particles, in addition to the somatic histones H2A, H2B, H3 and H4. The acid-soluble proteins of the oligonucleosomes comprised all the histones, including both the somatic (H1, H2A, H2B, H3, H4 and X2) and the testis-specific ones (TH1 and TH2B). It was also observed that histones TH1 and H1 are absent from the core particles and were readily extracted from the chromatin by 0.6 M-NaCl, which indicated that both of them are bound to the linker DNA.

A genome-wide association study of caffeine-related sleep disturbance: Confirmation of a role for a common variant in the adenosine receptor

OBJECTIVES To identify common genetic variants that predispose to caffeine-induced insomnia and to test whether genes whose expression changes in the presence of caffeine are enriched for association with caffeine-induced insomnia. DESIGN A hypothesis-free, genome-wide association study. SETTING Community-based sample of Australian twins from the Australian Twin Registry. PARTICIPANTS After removal of individuals who said that they do not drink coffee, a total of 2,402 individuals from 1,470 families in the Australian Twin Registry provided both phenotype and genotype information. MEASUREMENTS AND RESULTS A dichotomized scale based on whether participants reported ever or never experiencing caffeine-induced insomnia. A factor score based on responses to a number of questions regarding normal sleep habits was included as a covariate in the analysis. More than 2 million common single nucleotide polymorphisms (SNPs) were tested for association with caffeine-induced insomnia. No SNPs reached the genome-wide significance threshold. In the analysis that did not include the insomnia factor score as a covariate, the most significant SNP identified was an intronic SNP in the PRIMA1 gene (P = 1.4 x 10(-)(6), odds ratio = 0.68 [0.53 - 0.89]). An intergenic SNP near the GBP4 gene on chromosome 1 was the most significant upon inclusion of the insomnia factor score into the model (P = 1.9 x 10(-)(6), odds ratio = 0.70 [0.62 - 0.78]). A previously identified association with a polymorphism in the ADORA2A gene was replicated. CONCLUSIONS Several genes have been identified in the study as potentially influencing caffeine-induced insomnia. They will require replication in another sample. The results may have implications for understanding the biologic mechanisms underlying insomnia.

A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure

The ratio of the lengths of an individual's second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.