Bioinformatyczna analiza struktury i oddziaływań z substratami metylotransferaz RNA odpowiedzialnych za oporność bakterii na antybiotyki

Wydział Biologii: Instytut Biologii Molekularnej i Biotechnologii

Cystic fibrosis gene repair: correction of ΔF508 using ZFN and CRISPR/Cas9 guide RNA gene editing tools

Cystic Fibrosis (CF) is an autosomal recessive monogenic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene with the ΔF508 mutation accounting for approximately 70% of all CF cases worldwide. This thesis investigates whether existing zinc finger nucleases designed in this lab and CRISPR/gRNAs designed in this thesis can mediate efficient homology-directed repair (HDR) with appropriate donor repair plasmids to correct CF-causing mutations in a CF cell line. Firstly, the most common mutation, ΔF508, was corrected using a pair of existing ZFNs, which cleave in intron 9, and the donor repair plasmid pITR-donor-XC, which contains the correct CTT sequence and two unique restriction sites. HDR was initially determined to be <1% but further analysis by next generation sequencing (NGS) revealed HDR occurred at a level of 2%. This relatively low level of repair was determined to be a consequence of distance from the cut site to the mutation and so rather than designing a new pair of ZFNs, the position of the existing intron 9 ZFNs was exploited and attempts made to correct >80% of CF-causing mutations. The ZFN cut site was used as the site for HDR of a mini-gene construct comprising exons 10-24 from CFTR cDNA (with appropriate splice acceptor and poly A sites) to allow production of full length corrected CFTR mRNA. Finally, the ability to cleave closer to the mutation and mediate repair of CFTR using the latest gene editing tool CRISPR/Cas9 was explored. Two CRISPR gRNAs were tested; CRISPR ex10 was shown to cleave at an efficiency of 15% and CRISPR in9 cleaved at 3%. Both CRISPR gRNAs mediated HDR with appropriate donor plasmids at a rate of ~1% as determined by NGS. This is the first evidence of CRISPR induced HDR in CF cell lines.

Identification of RNA editing in the human exome and development of DARNED DatabaseSF

RNA editing is a biological phenomena that alters nascent RNA transcripts by insertion, deletion and/or substitution of one or a few nucleotides. It is ubiquitous in all kingdoms of life and in viruses. The predominant editing event in organisms with a developed central nervous system is Adenosine to Inosine deamination. Inosine is recognized as Guanosine by the translational machinery and reverse-transcriptase. In primates, RNA editing occurs frequently in transcripts from repetitive regions of the genome. In humans, more than 500,000 editing instances have been identified, by applying computational pipelines on available ESTs and high-throughput sequencing data, and by using chemical methods. However, the functions of only a small number of cases have been studied thoroughly. RNA editing instances have been found to have roles in peptide variants synthesis by non-synonymous codon substitutions, transcript variants by alterations in splicing sites and gene silencing by miRNAs sequence modifications. We established the Database of RNA EDiting (DARNED) to accommo-date the reference genomic coordinates of substitution editing in human, mouse and fly transcripts from published literatures, with additional information on edited genomic coordinates collected from various databases e.g. UCSC, NCBI. DARNED contains mostly Adenosine to Inosine editing and allows searches based on genomic region, gene ID, and user provided sequence. The Database is accessible at http://darned.ucc.ie RNA editing instances in coding region are likely to result in recoding in protein synthesis. This encouraged me to focus my research on the occurrences of RNA editing specific CDS and non-Alu exonic regions. By applying various filters on discrepancies between available ESTs and their corresponding reference genomic sequences, putative RNA editing candidates were identified. High-throughput sequencing was used to validate these candidates. All predicted coordinates appeared to be either SNPs or unedited.

Premature ovarian aging in mice deficient for Gpr3

After becoming competent for resuming meiosis, fully developed mammalian oocytes are maintained arrested in prophase I until ovulation is triggered by the luteotropin surge. Meiotic pause has been shown to depend critically on maintenance of cAMP level in the oocyte and was recently attributed to the constitutive Gs (the heterotrimeric GTP-binding protein that activates adenylyl cyclase) signaling activity of the G protein-coupled receptor GPR3. Here we show that mice deficient for Gpr3 are unexpectedly fertile but display progressive reduction in litter size despite stable age-independent alteration of meiotic pause. Detailed analysis of the phenotype confirms premature resumption of meiosis, in vivo, in about one-third of antral follicles from Gpr3-/- females, independently of their age. In contrast, in aging mice, absence of GPR3 leads to severe reduction of fertility, which manifests by production of an increasing number of nondeveloping early embryos upon spontaneous ovulation and massive amounts of fragmented oocytes after superovulation. Severe worsening of the phenotype in older animals points to an additional role of GPR3 related to protection (or rescue) of oocytes from aging. Gpr3-defective mice may constitute a relevant model of premature ovarian failure due to early oocyte aging.

Traumatic brain injury exacerbates neurodegenerative pathology: improvement with an apolipoprotein E-based therapeutic.

Cognitive impairment is common following traumatic brain injury (TBI), and neuroinflammatory mechanisms may predispose to the development of neurodegenerative disease. Apolipoprotein E (apoE) polymorphisms modify neuroinflammatory responses, and influence both outcome from acute brain injury and the risk of developing neurodegenerative disease. We demonstrate that TBI accelerates neurodegenerative pathology in double-transgenic animals expressing the common human apoE alleles and mutated amyloid precursor protein, and that pathology is exacerbated in the presence of the apoE4 allele. The administration of an apoE-mimetic peptide markedly reduced the development of neurodegenerative pathology in mice homozygous for apoE3 as well as apoE3/E4 heterozygotes. These results demonstrate that TBI accelerates the cardinal neuropathological features of neurodegenerative disease, and establishes the potential for apoE mimetic therapies in reducing pathology associated with neurodegeneration.

Microbial inactivation of Pseudomonas putida and Pichia pastoris using gene silencing.

Antisense deoxyoligonucleotide (ASO) gene silencing was investigated as a potential disinfection tool for industrial and drinking water treatment application. ASOs bind with their reverse complementary mRNA transcripts thereby blocking protein translation. While ASO silencing has mainly been studied in medicine, it may be useful for modulating gene expression and inactivating microorganisms in environmental applications. In this proof of concept work, gene targets were sh ble (zeocin resistance) and todE (catechol-2,3-dioxygenase) in Pichia pastoris and npt (kanamycin resistance) in Pseudomonas putida. A maximum 0.5-fold decrease in P. pastoris cell numbers was obtained following a 120 min incubation with single-stranded DNA (ssDNA) concentrations ranging from 0.2 to 200 nM as compared to the no ssDNA control. In P. putida, a maximum 5.2-fold decrease was obtained after 90 min with 400 nM ssDNA. While the silencing efficiencies varied for the 25 targets tested, these results suggest that protein activity as well as microbial growth can be altered using ASO gene silencing-based tools. If successful, this technology has the potential to eliminate some of the environmental and health issues associated with the use of strong chemical biocides. However, prior to its dissemination, more research is needed to increase silencing efficiency and develop effective delivery methods.

permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.

MRP3: a molecular target for human glioblastoma multiforme immunotherapy.

BACKGROUND: Glioblastoma multiforme (GBM) is refractory to conventional therapies. To overcome the problem of heterogeneity, more brain tumor markers are required for prognosis and targeted therapy. We have identified and validated a promising molecular therapeutic target that is expressed by GBM: human multidrug-resistance protein 3 (MRP3). METHODS: We investigated MRP3 by genetic and immunohistochemical (IHC) analysis of human gliomas to determine the incidence, distribution, and localization of MRP3 antigens in GBM and their potential correlation with survival. To determine MRP3 mRNA transcript and protein expression levels, we performed quantitative RT-PCR, raising MRP3-specific antibodies, and IHC analysis with biopsies of newly diagnosed GBM patients. We used univariate and multivariate analyses to assess the correlation of RNA expression and IHC of MRP3 with patient survival, with and without adjustment for age, extent of resection, and KPS. RESULTS: Real-time PCR results from 67 GBM biopsies indicated that 59/67 (88%) samples highly expressed MRP3 mRNA transcripts, in contrast with minimal expression in normal brain samples. Rabbit polyvalent and murine monoclonal antibodies generated against an extracellular span of MRP3 protein demonstrated reactivity with defined MRP3-expressing cell lines and GBM patient biopsies by Western blotting and FACS analyses, the latter establishing cell surface MRP3 protein expression. IHC evaluation of 46 GBM biopsy samples with anti-MRP3 IgG revealed MRP3 in a primarily membranous and cytoplasmic pattern in 42 (91%) of the 46 samples. Relative RNA expression was a strong predictor of survival for newly diagnosed GBM patients. Hazard of death for GBM patients with high levels of MRP3 RNA expression was 2.71 (95% CI: 1.54-4.80) times that of patients with low/moderate levels (p = 0.002). CONCLUSIONS: Human GBMs overexpress MRP3 at both mRNA and protein levels, and elevated MRP3 mRNA levels in GBM biopsy samples correlated with a higher risk of death. These data suggest that the tumor-associated antigen MRP3 has potential use for prognosis and as a target for malignant glioma immunotherapy.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

An integrated approach to the prediction of chemotherapeutic response in patients with breast cancer.

BACKGROUND: A major challenge in oncology is the selection of the most effective chemotherapeutic agents for individual patients, while the administration of ineffective chemotherapy increases mortality and decreases quality of life in cancer patients. This emphasizes the need to evaluate every patient's probability of responding to each chemotherapeutic agent and limiting the agents used to those most likely to be effective. METHODS AND RESULTS: Using gene expression data on the NCI-60 and corresponding drug sensitivity, mRNA and microRNA profiles were developed representing sensitivity to individual chemotherapeutic agents. The mRNA signatures were tested in an independent cohort of 133 breast cancer patients treated with the TFAC (paclitaxel, 5-fluorouracil, adriamycin, and cyclophosphamide) chemotherapy regimen. To further dissect the biology of resistance, we applied signatures of oncogenic pathway activation and performed hierarchical clustering. We then used mRNA signatures of chemotherapy sensitivity to identify alternative therapeutics for patients resistant to TFAC. Profiles from mRNA and microRNA expression data represent distinct biologic mechanisms of resistance to common cytotoxic agents. The individual mRNA signatures were validated in an independent dataset of breast tumors (P = 0.002, NPV = 82%). When the accuracy of the signatures was analyzed based on molecular variables, the predictive ability was found to be greater in basal-like than non basal-like patients (P = 0.03 and P = 0.06). Samples from patients with co-activated Myc and E2F represented the cohort with the lowest percentage (8%) of responders. Using mRNA signatures of sensitivity to other cytotoxic agents, we predict that TFAC non-responders are more likely to be sensitive to docetaxel (P = 0.04), representing a viable alternative therapy. CONCLUSIONS: Our results suggest that the optimal strategy for chemotherapy sensitivity prediction integrates molecular variables such as ER and HER2 status with corresponding microRNA and mRNA expression profiles. Importantly, we also present evidence to support the concept that analysis of molecular variables can present a rational strategy to identifying alternative therapeutic opportunities.

Gene expression disruptions of organism versus organ in Drosophila species hybrids.

Hybrid dysfunctions, such as sterility, may result in part from disruptions in the regulation of gene expression. Studies of hybrids within the Drosophila simulans clade have reported genes expressed above or below the expression observed in their parent species, and such misexpression is associated with male sterility in multigenerational backcross hybrids. However, these studies often examined whole bodies rather than testes or had limited replication using less-sensitive but global techniques. Here, we use a new RNA isolation technique to re-examine hybrid gene expression disruptions in both testes and whole bodies from single Drosophila males by real-time quantitative RT-PCR. We find two early-spermatogenesis transcripts are underexpressed in hybrid whole-bodies but not in assays of testes alone, while two late-spermatogenesis transcripts seem to be underexpressed in both whole-bodies and testes alone. Although the number of transcripts surveyed is limited, these results provide some support for a previous hypothesis that the spermatogenesis pathway in these sterile hybrids may be disrupted sometime after the expression of the early meiotic arrest genes.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

A switch in the control of growth of the wing imaginal disks of Manduca sexta.

BACKGROUND: Insulin and ecdysone are the key extrinsic regulators of growth for the wing imaginal disks of insects. In vitro tissue culture studies have shown that these two growth regulators act synergistically: either factor alone stimulates only limited growth, but together they stimulate disks to grow at a rate identical to that observed in situ. It is generally thought that insulin signaling links growth to nutrition, and that starvation stops growth because it inhibits insulin secretion. At the end of larval life feeding stops but the disks continue to grow, so at that time disk growth has become uncoupled from nutrition. We sought to determine at exactly what point in development this uncoupling occurs. METHODOLOGY: Growth and cell proliferation in the wing imaginal disks and hemolymph carbohydrate concentrations were measured at various stages in the last larval instar under experimental conditions of starvation, ligation, rescue, and hormone treatment. PRINCIPAL FINDINGS: Here we show that in the last larval instar of M. sexta, the uncoupling of nutrition and growth occurs as the larva passes the critical weight. Before this time, starvation causes a decline in hemolymph glucose and trehalose and a cessation of wing imaginal disks growth, which can be rescued by injections of trehalose. After the critical weight the trehalose response to starvation disappears, and the expression of insulin becomes decoupled from nutrition. After the critical weight the wing disks loose their sensitivity to repression by juvenile hormone, and factors from the abdomen, but not the brain, are required to drive continued growth. CONCLUSIONS: During the last larval instar imaginal disk growth becomes decoupled from somatic growth at the time that the endocrine events of metamorphosis are initiated. These regulatory changes ensure that disk growth continues uninterrupted when the nutritive and endocrine signals undergo the drastic changes associated with metamorphosis.

The effects of sleep hypoxia on coagulant factors and hepatic inflammation in emphysematous rats.

OBJECTIVES: To develop a sleep hypoxia (SH) in emphysema (SHE) rat model and to explore whether SHE results in more severe hepatic inflammation than emphysema alone and whether the inflammation changes levels of coagulant/anticoagulant factors synthesized in the liver. METHODS: Seventy-five rats were put into 5 groups: SH control (SHCtrl), treated with sham smoke exposure (16 weeks) and SH exposure (12.5% O(2), 3 h/d, latter 8 weeks); emphysema control (ECtrl), smoke exposure and sham SH exposure (21% O(2)); short SHE (SHEShort), smoke exposure and short SH exposure (1.5 h/d); mild SHE (SHEMild), smoke exposure and mild SH exposure (15% O(2)); standard SHE (SHEStand), smoke exposure and SH exposure. Therefore, ECtrl, SHEShort, SHEMild and SHEStand group were among emphysematous groups. Arterial blood gas (ABG) data was obtained during preliminary tests. After exposure, hepatic inflammation (interleukin -6 [IL-6] mRNA and protein, tumor necrosis factor α [TNFα] mRNA and protein) and liver coagulant/anticoagulant factors (antithrombin [AT], fibrinogen [FIB] and Factor VIII [F VIII]) were evaluated. SPSS 11.5 software was used for statistical analysis. RESULTS: Characteristics of emphysema were obvious in emphysematous groups and ABGs reached SH criteria on hypoxia exposure. Hepatic inflammation parameters and coagulant factors are the lowest in SHCtrl and the highest in SHEStand while AT is the highest in SHCtrl and the lowest in SHEStand. Inflammatory cytokines of liver correlate well with coagulant factors positively and with AT negatively. CONCLUSIONS: When SH is combined with emphysema, hepatic inflammation and coagulability enhance each other synergistically and produce a more significant liver-derivative inflammatory and prothrombotic status.