A novel approach to design of cis-acting DNA structural element for regulation of gene expression in vivo

Taking advantage of the degeneracy of the genetic code we have developed a novel approach to introduce, within a gene, DNA sequences capable of adopting unusual structures and to investigate the role of such sequences in regulation of gene expression in vivo. We used a computer program that generates alternative codon sequences for the same amino-acid sequence to convert a stretch of nucleotides into an inverted-repeat sequence with the potential to adopt cruciform structure. This approach was used to replace a 51-base-pair EcoRI-HindIII segment in the N-terminal region of the beta-galactosidase gene in plasmid pUC19 with a 51-bp synthetic oligonucleotide sequence with the potential to adopt a cruciform structure with 18 bp in the stem region. In selecting the 51-bp sequence, care was taken to include those codons that are preferred in E. coli. E. coli DH5-alpha cells harbouring the plasmid containing the redesigned sequence showed drastic reduction in expression of the beta-galactosidase gene compared to cells harbouring the plasmid with the native sequence. This approach demonstrates the possibility of introducing DNA secondary-structure elements to alter regulation of gene expression in vivo.

Simulation and graph mining tools for improving gene mapping efficiency

Gene mapping is a systematic search for genes that affect observable characteristics of an organism. In this thesis we offer computational tools to improve the efficiency of (disease) gene-mapping efforts. In the first part of the thesis we propose an efficient simulation procedure for generating realistic genetical data from isolated populations. Simulated data is useful for evaluating hypothesised gene-mapping study designs and computational analysis tools. As an example of such evaluation, we demonstrate how a population-based study design can be a powerful alternative to traditional family-based designs in association-based gene-mapping projects. In the second part of the thesis we consider a prioritisation of a (typically large) set of putative disease-associated genes acquired from an initial gene-mapping analysis. Prioritisation is necessary to be able to focus on the most promising candidates. We show how to harness the current biomedical knowledge for the prioritisation task by integrating various publicly available biological databases into a weighted biological graph. We then demonstrate how to find and evaluate connections between entities, such as genes and diseases, from this unified schema by graph mining techniques. Finally, in the last part of the thesis, we define the concept of reliable subgraph and the corresponding subgraph extraction problem. Reliable subgraphs concisely describe strong and independent connections between two given vertices in a random graph, and hence they are especially useful for visualising such connections. We propose novel algorithms for extracting reliable subgraphs from large random graphs. The efficiency and scalability of the proposed graph mining methods are backed by extensive experiments on real data. While our application focus is in genetics, the concepts and algorithms can be applied to other domains as well. We demonstrate this generality by considering coauthor graphs in addition to biological graphs in the experiments.

The effect of progesterone replacement on gene expression in the corpus luteum during induced regression and late luteal phase in the bonnet monkey (Macaca radiata)

Background: In higher primates, although LH/CG play a critical role in the control of corpus luteum (CL) function, the direct effects of progesterone (P4) in the maintenance of CL structure and function are unclear. Several experiments were conducted in the bonnet monkey to examine direct effects of P4 on gene expression changes in the CL, during induced luteolysis and the late luteal phase of natural cycles. Methods: To identify differentially expressed genes encoding PR, PR binding factors, cofactors and PR downstream signaling target genes, the genome-wide analysis data generated in CL of monkeys after LH/P-4 depletion and LH replacement were mined and validated by real-time RT-PCR analysis. Initially, expression of these P4 related genes were determined in CL during different stages of luteal phase. The recently reported model system of induced luteolysis, yet capable of responsive to tropic support, afforded an ideal situation to examine direct effects of P4 on structure and function of CL. For this purpose, P4 was infused via ALZET pumps into monkeys 24 h after LH/P4 depletion to maintain mid luteal phase circulating P4 concentration (P4 replacement). In another experiment, exogenous P4 was supplemented during late luteal phase to mimic early pregnancy. Results: Based on the published microarray data, 45 genes were identified to be commonly regulated by LH and P4. From these 19 genes belonging to PR signaling were selected to determine their expression in LH/P-4 depletion and P4 replacement experiments. These 19 genes when analyzed revealed 8 genes to be directly responsive to P4, whereas the other genes to be regulated by both LH and P4. Progesterone supplementation for 24 h during the late luteal phase also showed changes in expression of 17 out of 19 genes examined. Conclusion: These results taken together suggest that P4 regulates, directly or indirectly, expression of a number of genes involved in the CL structure and function.

Post-transcriptional Repair of a Split Heat Shock Protein 90 Gene by mRNA trans-Splicing

Heat shock protein 90 participates in diverse biological processes ranging from protein folding, cell cycle, signal transduction and development to evolution in all eukaryotes. It is also critically involved in regulating growth of protozoa such as Dictyostelium discoideum, Leishmania donovani, Plasmodium falciparum, Trypanosoma cruzi, and Trypanosoma evansi. Selective inhibition of Hsp90 has also been explored as an intervention strategy against important human diseases such as cancer, malaria, or trypanosomiasis. Giardia lamblia, a simple protozoan parasite of humans and animals, is an important cause of diarrheal disease with significant morbidity and some mortality in tropical countries. Here we show that the G. lamblia cytosolic hsp90 ( glhsp90) is split in two similar sized fragments located 777 kb apart on the same scaffold. Intrigued by this unique arrangement, which appears to be specific for the Giardiinae, we have investigated the biosynthesis of GlHsp90. We used genome sequencing to confirm the split nature of the giardial hsp90. However, a specific antibody raised against the peptide detected a product with a mass of about 80 kDa, suggesting a post-transcriptional rescue of the genomic defect. We show evidence for the joining of the two independent Hsp90 transcripts in-trans to one long mature mRNA presumably by RNA splicing. The splicing junction carries hallmarks of classical cis-spliced introns, suggesting that the regular cis-splicing machinery may be sufficient for repair of the open reading frame. A complementary 26-nt sequence in the ``intron'' regions adjacent to the splice sites may assist in positioning the two pre-mRNAs for processing. This is the first example of post-transcriptional rescue of a split gene by trans-splicing.

Advantage of the ether linkage between the positive charge and the cholesteryl skeleton in cholesterol-based amphiphiles as vectors for gene delivery

Twelve novel cationic cholesterol derivatives with different linkage types between the cationic headgroup and the cholesteryl backbone have been developed. These have been tested for their efficacies as gene transfer agents as mixtures with dioleoyl phosphatidylethanolamine (DOPE). A pronounced improvement in transfection efficiency was observed when the cationic center was linked to the steroid backbone using an ether type bond. Among these, cholest-5-en-3b-oxyethane-N, N,N-trimethylammonium bromide (2a) and cholest-5-en-3b-oxyethane-N, N-dimethyl-N-2-hydroxyethylammonium bromide (3d) showed transfection efficiencies considerably greater than commercially available reagents such as Lipofectin or Lipofectamine. To achieve transfection, 3d did not require DOPE. Increasing hydration at the headgroup level for both ester- and ether-linked amphiphiles resulted in progressive loss of transfection efficiency. Transfection efficiency was also greatly reduced when a 'disorder'-inducing chain like an oleyl (cis-9-octadecenyl) segment was added to these cholesteryl amphiphiles. Importantly, the transfection ability of 2a with DOPE in the presence of serum was significantly greater than for a commercially available reagent, Lipofectamine. This suggests that these novel cholesterol-based amphiphiles might prove promising in applications involving liposome-mediated gene transfection. This investigation demonstrates the importance of structural features at the molecular level for the design of cholesterol-based gene delivery reagents that would aid the development of newer, more efficient formulations based on this class of molecules.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.

Efficient Dendrimer-DNA Complexation and Gene Delivery Vector Properties of Nitrogen-Core Poly(propyl ether imine) Dendrimer in Mammalian Cells

Dendrimers as vectors for gene delivery were established, primarily by utilizing few prominent dendrimer types so far. We report herein studies of DNA complexation efficacies and gene delivery vector properties of a nitrogen-core poly(propyl ether imine) (PETIM) dendrimer, constituted with 22 tertiary amine internal branches and 24 primary amines at the periphery. The interaction of the dendrimer with pEGFPDNA was evaluated through UV-vis, circular dichroism (CD) spectral studies, ethidium bromide fluorescence emission quenching, thermal melting, and gel retardation assays, from which most changes to DNA structure during complexation was found to occur at a weight ratio of dendrimer:DNA similar to 2:1. The zeta potential measurements further confirmed this stoichiometry at electroneutrality. The structure of a DNA oligomer upon dendrimer complexation was simulated through molecular modeling and the simulation showed that the dendrimer enfolded DNA oligomer along both major and minor grooves, without causing DNA deformation, in 1:1 and 2:1 dendrimer-to-DNA complexes. Atomic force microscopy (AFM) studies on dendrimer-pEGFP DNA complex showed an increase in the average z-height as a result of dendrimers decorating the DNA, without causing a distortion of the DNA structure. Cytotoxicity studies involving five different mammalian cell lines, using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide] (MTT) assay, reveal the dendrimer toxicity profile (IC50) values of similar to 400-1000 mu g mL(-1), depending on the cell line tested. Quantitative estimation, using luciferase assay, showed that the gene transfection was at least 100 times higher when compared to poly(ethylene imine) branched polymer, having similar number of cationic sites as the dendrimer. The present study establishes the physicochemical behavior of new nitrogen-core PETIM dendrimer-DNA complexes, their lower toxicities, and efficient gene delivery vector properties.

A 16-Gene Signature Distinguishes Anaplastic Astrocytoma from Glioblastoma

Anaplastic astrocytoma (AA; Grade III) and glioblastoma (GBM; Grade IV) are diffusely infiltrating tumors and are called malignant astrocytomas. The treatment regimen and prognosis are distinctly different between anaplastic astrocytoma and glioblastoma patients. Although histopathology based current grading system is well accepted and largely reproducible, intratumoral histologic variations often lead to difficulties in classification of malignant astrocytoma samples. In order to obtain a more robust molecular classifier, we analysed RT-qPCR expression data of 175 differentially regulated genes across astrocytoma using Prediction Analysis of Microarrays (PAM) and found the most discriminatory 16-gene expression signature for the classification of anaplastic astrocytoma and glioblastoma. The 16-gene signature obtained in the training set was validated in the test set with diagnostic accuracy of 89%. Additionally, validation of the 16-gene signature in multiple independent cohorts revealed that the signature predicted anaplastic astrocytoma and glioblastoma samples with accuracy rates of 99%, 88%, and 92% in TCGA, GSE1993 and GSE4422 datasets, respectively. The protein-protein interaction network and pathway analysis suggested that the 16-genes of the signature identified epithelial-mesenchymal transition (EMT) pathway as the most differentially regulated pathway in glioblastoma compared to anaplastic astrocytoma. In addition to identifying 16 gene classification signature, we also demonstrated that genes involved in epithelial-mesenchymal transition may play an important role in distinguishing glioblastoma from anaplastic astrocytoma.

Predicting gene ontology annotations of orphan GWAS genes using protein-protein interactions

Background: The number of genome-wide association studies (GWAS) has increased rapidly in the past couple of years, resulting in the identification of genes associated with different diseases. The next step in translating these findings into biomedically useful information is to find out the mechanism of the action of these genes. However, GWAS studies often implicate genes whose functions are currently unknown; for example, MYEOV, ANKLE1, TMEM45B and ORAOV1 are found to be associated with breast cancer, but their molecular function is unknown. Results: We carried out Bayesian inference of Gene Ontology (GO) term annotations of genes by employing the directed acyclic graph structure of GO and the network of protein-protein interactions (PPIs). The approach is designed based on the fact that two proteins that interact biophysically would be in physical proximity of each other, would possess complementary molecular function, and play role in related biological processes. Predicted GO terms were ranked according to their relative association scores and the approach was evaluated quantitatively by plotting the precision versus recall values and F-scores (the harmonic mean of precision and recall) versus varying thresholds. Precisions of similar to 58% and similar to 40% for localization and functions respectively of proteins were determined at a threshold of similar to 30 (top 30 GO terms in the ranked list). Comparison with function prediction based on semantic similarity among nodes in an ontology and incorporation of those similarities in a k nearest neighbor classifier confirmed that our results compared favorably. Conclusions: This approach was applied to predict the cellular component and molecular function GO terms of all human proteins that have interacting partners possessing at least one known GO annotation. The list of predictions is available at http://severus.dbmi.pitt.edu/engo/GOPRED.html. We present the algorithm, evaluations and the results of the computational predictions, especially for genes identified in GWAS studies to be associated with diseases, which are of translational interest.

Targeted Modifications in Adeno-Associated Virus Serotype 8 Capsid Improves Its Hepatic Gene Transfer Efficiency In Vivo

Recombinant adeno-associated virus vectors based on serotype 8 (AAV8) have shown significant promise for liver-directed gene therapy. However, to overcome the vector dose dependent immunotoxicity seen with AAV8 vectors, it is important to develop better AAV8 vectors that provide enhanced gene expression at significantly low vector doses. Since it is known that AAV vectors during intracellular trafficking are targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal machinery, we modified specific serine/threonine kinase or ubiquitination targets on the AAV8 capsid to augment its transduction efficiency. Point mutations at specific serine (S)/threonine (T)/lysine (K) residues were introduced in the AAV8 capsid at the positions equivalent to that of the effective AAV2 mutants, generated successfully earlier. Extensive structure analysis was carried out subsequently to evaluate the structural equivalence between the two serotypes. scAAV8 vectors with the wild-type (WT) and each one of the S/T -> Alanine (A) or K-Arginine (R) mutant capsids were evaluated for their liver transduction efficiency in C57BL/6 mice in vivo. Two of the AAV8-S -> A mutants (S279A and S671A), and a K137R mutant vector, demonstrated significantly higher enhanced green fluorescent protein (EGFP) transcript levels (similar to 9- to 46-fold) in the liver compared to animals that received WT-AAV8 vectors alone. The best performing AAV8 mutant (K137R) vector also had significantly reduced ubiquitination of the viral capsid, reduced activation of markers of innate immune response, and a concomitant two-fold reduction in the levels of neutralizing antibody formation in comparison to WT-AAV8 vectors. Vector bio-distribution studies revealed that the K137R mutant had a significantly higher and preferential transduction of the liver (106 vs. 7.7 vector copies/mouse diploid genome) when compared to WT-AAV8 vectors. To further study the utility of the K137R-AAV8 mutant in therapeutic gene transfer, we delivered human coagulation factor IX (h. FIX) under the control of liver-specific promoters (LP1 or hAAT) into C57BL/6 mice. The circulating levels of h. FIX: Ag were higher in all the K137R-AAV8 treated groups up to 8 weeks post-hepatic gene transfer. These studies demonstrate the feasibility of the use of this novel AAV8 vectors for potential gene therapy of hemophilia B.

A novel gene THSD7A is associated with obesity

Body mass index (BMI) is a non-invasive measurement of obesity. It is commonly used for assessing adiposity and obesity-related risk prediction. Genetic differences between ethnic groups are important factors, which contribute to the variation in phenotypic effects. India inhabited by the first out-of-Africa human population and the contemporary Indian populations are admixture of two ancestral populations; ancestral north Indians (ANI) and ancestral south Indians (ASI). Although ANI are related to Europeans, ASI are not related to any group outside Indian-subcontinent. Hence, we expect novel genetic loci associated with BMI. In association analysis, we found eight genic SNPs in extreme of distribution (P <= 3.75 x 10(-5)), of which WWOX has already been reported to be associated with obesity-related traits hence excluded from further study. Interestingly, we observed rs1526538, an intronic SNP of THSD7A; a novel gene significantly associated with obesity (P = 2.88 x 10(-5), 8.922 x 10(-6) and 2.504 x 10(-9) in discovery, replication and combined stages, respectively). THSD7A is neural N-glycoprotein, which promotes angiogenesis and it is well known that angiogenesis modulates obesity, adipose metabolism and insulin sensitivity, hence our result find a correlation. This information can be used for drug target, early diagnosis of obesity and treatment.

Regulation of Growth, Cell Shape, Cell Division, and Gene Expression by Second Messengers (p)ppGpp and Cyclic Di-GMP in Mycobacterium smegmatis

The alarmone (p)ppGpp regulates transcription, translation, replication, virulence, lipid synthesis, antibiotic sensitivity, biofilm formation, and other functions in bacteria. Signaling nucleotide cyclic di-GMP (c-di-GMP) regulates biofilm formation, motility, virulence, the cell cycle, and other functions. In Mycobacterium smegmatis, both (p) ppGpp and c-di-GMP are synthesized and degraded by bifunctional proteins Rel(Msm) and DcpA, encoded by rel(Msm) and dcpA genes, respectively. We have previously shown that the Delta rel(Msm) and Delta dcpA knockout strains are antibiotic resistant and defective in biofilm formation, show altered cell surface properties, and have reduced levels of glycopeptidolipids and polar lipids in their cell wall (K. R. Gupta, S. Kasetty, and D. Chatterji, Appl Environ Microbiol 81:2571-2578, 2015, http://dx.doi.org/10.1128/AEM.03999-14). In this work, we have explored the phenotypes that are affected by both (p) ppGpp and c-di-GMP in mycobacteria. We have shown that both (p) ppGpp and c-di-GMP are needed to maintain the proper growth rate under stress conditions such as carbon deprivation and cold shock. Scanning electron microscopy showed that low levels of these second messengers result in elongated cells, while high levels reduce the cell length and embed the cells in a biofilm-like matrix. Fluorescence microscopy revealed that the elongated Delta rel(Msm) and Delta dcpA cells are multinucleate, while transmission electron microscopy showed that the elongated cells are multiseptate. Gene expression analysis also showed that genes belonging to functional categories such as virulence, detoxification, lipid metabolism, and cell-wall-related processes were differentially expressed. Our results suggests that both (p) ppGpp and c-di-GMP affect some common phenotypes in M. smegmatis, thus raising a possibility of cross talk between these two second messengers in mycobacteria. IMPORTANCE Our work has expanded the horizon of (p) ppGpp and c-di-GMP signaling in Gram-positive bacteria. We have come across a novel observation that M. smegmatis needs (p) ppGpp and c-di-GMP for cold tolerance. We had previously shown that the Delta rel(Msm) and Delta dcpA strains are defective in biofilm formation. In this work, the overproduction of (p) ppGpp and c-di-GMP encased M. smegmatis in a biofilm-like matrix, which shows that both (p) ppGpp and c-di-GMP are needed for biofilm formation. The regulation of cell length and cell division by (p) ppGpp was known in mycobacteria, but our work shows that c-di-GMP also affects the cell size and cell division in mycobacteria. This is perhaps the first report of c-di-GMP regulating cell division in mycobacteria.

Caracterização de alterações epigenéticas no gene JARID1C e desequilíbrios genéticos como causas do retardo mental ligado ao x de etiologia idiopática

O retardo mental (RM) é caracterizado por um funcionamento intelectual significantemente abaixo da média (QI<70). A prevalência de RM varia entre estudos epidemiológicos, sendo estimada em 2-3% da população mundial, constituindo assim, um dos mais importantes problemas de saúde pública. Há um consenso geral de que o RM é mais comum no sexo masculino, um achado atribuído às numerosas mutações nos genes encontrados no cromossomo X, levando ao retardo mental ligado ao X (RMLX). Dentre os genes presentes no cromossomo X, o Jumonji AT-rich interactive domain IC (JARID1C) foi recentemente identificado como um potencial candidato etiológico do RM, quando mutado. O JARID1C codifica uma proteína que atua como uma desmetilase da lisina 4 da histona H3 (H3K4), imprescindível para a regulação epigenética. Tão recente como a identificação do gene JARID1C, é a descoberta de que mudanças no número de cópias de sequências de DNA, caracterizadas por microdeleções e microduplicações, poderiam ser consideradas como razões funcionalmente importantes de RMLX. Atualmente, cerca de 5-10% dos casos de RM em homens são reconhecidos por ocorrerem devido a estas variações do número de cópias no cromossomo X. Neste estudo, investigamos mutações no gene JARID1C, através do rastreamento dos éxons 9, 11, 12, 13, 15 e 16, em 121 homens de famílias com RM provavelmente ligado ao X. Paralelamente, realizamos a análise da variação do número de cópias em 16 genes localizados no cromossomo X através da técnica de MLPA no mesmo grupo de pacientes. Esta metodologia consiste em uma amplificação múltipla que detecta variações no número de cópias de até 50 sequências diferentes de DNA genômico, sendo capaz de distinguir sequências que diferem em apenas um nucleotídeo. O DNA genômico foi extraído a partir de sangue periférico e as amostras foram amplificadas pela técnica de PCR, seguida da análise por sequenciamento direto. Foram identificadas três variantes na sequência do gene JARID1C entre os pacientes analisados: a variante intrônica 2243+11 G>T, que esteve presente em 67% dos pacientes, a variante silenciosa c.1794C>G e a mutação inédita nonsense c.2172C>A, ambas presentes em 0,82% dos indivíduos investigados. A análise através do MLPA revelou uma duplicação em um dos pacientes envolvendo as sondas referentes ao gene GDI1 e ao gene HUWE1. Este trabalho expande o estudo de mutações no gene JARID1C para a população brasileira ereforça a importância da triagem de mutações neste gene em homens portadores de RM familiar de origem idiopática, assim como, é primeiro relato científico relativo à investigação de variações no número de cópias de genes localizados no cromossomo X em homens brasileiros com RM, através da técnica de MLPA.