New computational biology tools for the systematic analysis of the structure and expression of human genes

From the late 1980s, the automation of sequencing techniques and the computer spread gave rise to a flourishing number of new molecular structures and sequences and to proliferation of new databases in which to store them. Here are presented three computational approaches able to analyse the massive amount of publicly avalilable data in order to answer to important biological questions. The first strategy studies the incorrect assignment of the first AUG codon in a messenger RNA (mRNA), due to the incomplete determination of its 5' end sequence. An extension of the mRNA 5' coding region was identified in 477 in human loci, out of all human known mRNAs analysed, using an automated expressed sequence tag (EST)-based approach. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for GNB2L1, QARS and TDP2 and the consequences for the functional studies are discussed. The second approach analyses the codon bias, the phenomenon in which distinct synonymous codons are used with different frequencies, and, following integration with a gene expression profile, estimates the total number of codons present across all the expressed mRNAs (named here "codonome value") in a given biological condition. Systematic analyses across different pathological and normal human tissues and multiple species shows a surprisingly tight correlation between the codon bias and the codonome bias. The third approach is useful to studies the expression of human autism spectrum disorder (ASD) implicated genes. ASD implicated genes sharing microRNA response elements (MREs) for the same microRNA are co-expressed in brain samples from healthy and ASD affected individuals. The different expression of a recently identified long non coding RNA which have four MREs for the same microRNA could disrupt the equilibrium in this network, but further analyses and experiments are needed.

Functional expression of CCR1, CCR3, CCR4, and CXCR4 chemokine receptors on human platelets

Platelets are known to contain platelet factor 4 and beta-thromboglobulin, alpha-chemokines containing the CXC motif, but recent studies extended the range to the beta-family characterized by the CC motif, including RANTES and Gro-alpha. There is also evidence for expression of chemokine receptors CCR4 and CXCR4 in platelets. This study shows that platelets have functional CCR1, CCR3, CCR4, and CXCR4 chemokine receptors. Polymerase chain reaction detected chemokine receptor messenger RNA in platelet RNA. CCR1, CCR3, and especially CCR4 gave strong signals; CXCR1 and CXCR4 were weakly positive. Flow cytometry with specific antibodies showed the presence of a clear signal for CXCR4 and weak signals for CCR1 and CCR3, whereas CXCR1, CXCR2, CXCR3, and CCR5 were all negative. Immunoprecipitation and Western blotting with polyclonal antibodies to cytoplasmic peptides clearly showed the presence of CCR1 and CCR4 in platelets in amounts comparable to monocytes and CCR4 transfected cells, respectively. Chemokines specific for these receptors, including monocyte chemotactic protein 1, macrophage inflammatory peptide 1alpha, eotaxin, RANTES, TARC, macrophage-derived chemokine, and stromal cell-derived factor 1, activate platelets to give Ca(++) signals, aggregation, and release of granule contents. Platelet aggregation was dependent on release of adenosine diphosphate (ADP) and its interaction with platelet ADP receptors. Part, but not all, of the Ca(++) signal was due to ADP release feeding back to its receptors. Platelet activation also involved heparan or chondroitin sulfate associated with the platelet surface and was inhibited by cleavage of these glycosaminoglycans or by heparin or low molecular weight heparin. These platelet receptors may be involved in inflammatory or allergic responses or in platelet activation in human immunodeficiency virus infection.

Comparative quantitative assessment of GnRH- and LH-receptor mRNA expression in the urinary tract of sexually intact and spayed female dogs

Ovariectomy interrupts the regulatory loop in the hypothalamus-pituitary-gonad axis, leading to a several-fold increase in gonadotropin levels. This rise in hormonal secretion may play a causal role in ovariectomy-related urinary incontinence. The purpose of this study was to examine the effect of ovariectomy in bitches on the expression of GnRH- and LH-receptors in the lower urinary tract, and assess the relationship between receptor expression and plasma gonadotropin concentrations. Plasma gonadotropins were measured in 37 client-owned bitches. Biopsies were harvested from the mid-ventral bladder wall in all dogs, and from nine further locations within the lower urinary tract in 17 of the 37 animals. Messenger RNA of the LH and GnRH receptors was quantified using RT-PCR with the TaqMan Universal PCR Master Mix. Gonadotropins were measured with a canine-specific FSH-immunoradiometric assay and LH-radioimmunoassay. The hierarchical mixed ANOVA model using MINITAB, Mann-Whitney U-test, unpaired means comparison and linear regressions using StatView were applied for statistical analyses. Messenger RNA for both receptors was detected in all biopsy samples. Age was negatively correlated to mRNA expression of the LH and the GnRH receptors. A relationship between the mRNA values and the plasma gonadotropin concentrations was not established. Evaluation of results within each of the biopsy locations revealed greater LH-receptor expression in the proximal second quarter of the urethra in spayed bitches than in intact bitches (P=0.0481). Increased mRNA expression of LH receptors in this location could possibly play a role in the decrease in closing pressure of the urethra following ovariectomy.

Gingival fluid cytokine expression and subgingival bacterial counts during pregnancy and postpartum: a case series

OBJECTIVES The aim of this study was to assess gingival fluid (GCF) cytokine messenger RNA (mRNA) levels, subgingival bacteria, and clinical periodontal conditions during a normal pregnancy to postpartum. MATERIALS AND METHODS Subgingival bacterial samples were analyzed with the checkerboard DNA-DNA hybridization method. GCF samples were assessed with real-time PCR including five proinflammatory cytokines and secretory leukocyte protease inhibitor. RESULTS Nineteen pregnant women with a mean age of 32 years (S.D. ± 4 years, range 26-42) participated in the study. Full-mouth bleeding scores (BOP) decreased from an average of 41.2% (S.D. ± 18.6%) at the 12th week of pregnancy to 26.6% (S.D. ± 14.4%) at the 4-6 weeks postpartum (p < 0.001). Between week 12 and 4-6 weeks postpartum, the mean probing pocket depth changed from 2.4 mm (S.D. ± 0.4) to 2.3 mm (S.D. ± 0.3) (p = 0.34). Higher counts of Eubacterium saburreum, Parvimonas micra, Selenomonas noxia, and Staphylococcus aureus were found at week 12 of pregnancy than at the 4-6 weeks postpartum examinations (p < 0.001). During and after pregnancy, statistically significant correlations between BOP scores and bacterial counts were observed. BOP scores and GCF levels of selected cytokines were not related to each other and no differences in GCF levels of the cytokines were observed between samples from the 12th week of pregnancy to 4-6 weeks postpartum. Decreasing postpartum counts of Porphyromonas endodontalis and Pseudomonas aeruginosa were associated with decreasing levels of Il-8 and Il-1β. CONCLUSIONS BOP decreased after pregnancy without any active periodontal therapy. Associations between bacterial counts and cytokine levels varied greatly in pregnant women with gingivitis and a normal pregnancy outcome. Postpartum associations between GCF cytokines and bacterial counts were more consistent. CLINICAL RELEVANCE Combined assessments of gingival fluid cytokines and subgingival bacteria may provide important information on host response.

RNAi reduces expression and intracellular retention of mutant cartilage oligomeric matrix protein.

Mutations in cartilage oligomeric matrix protein (COMP), a large extracellular glycoprotein expressed in musculoskeletal tissues, cause two skeletal dysplasias, pseudoachondroplasia and multiple epiphyseal dysplasia. These mutations lead to massive intracellular retention of COMP, chondrocyte death and loss of growth plate chondrocytes that are necessary for linear growth. In contrast, COMP null mice have only minor growth plate abnormalities, normal growth and longevity. This suggests that reducing mutant and wild-type COMP expression in chondrocytes may prevent the toxic cellular phenotype causing the skeletal dysplasias. We tested this hypothesis using RNA interference to reduce steady state levels of COMP mRNA. A panel of shRNAs directed against COMP was tested. One shRNA (3B) reduced endogenous and recombinant COMP mRNA dramatically, regardless of expression levels. The activity of the shRNA against COMP mRNA was maintained for up to 10 weeks. We also demonstrate that this treatment reduced ER stress. Moreover, we show that reducing steady state levels of COMP mRNA alleviates intracellular retention of other extracellular matrix proteins associated with the pseudoachondroplasia cellular pathology. These findings are a proof of principle and the foundation for the development of a therapeutic intervention based on reduction of COMP expression.

Coordinated changes in mRNA turnover, translation, and RNA processing bodies in bronchial epithelial cells following inflammatory stimulation.

Bronchial epithelial cells play a pivotal role in airway inflammation, but little is known about posttranscriptional regulation of mediator gene expression during the inflammatory response in these cells. Here, we show that activation of human bronchial epithelial BEAS-2B cells by proinflammatory cytokines interleukin-4 (IL-4) and tumor necrosis factor alpha (TNF-alpha) leads to an increase in the mRNA stability of the key chemokines monocyte chemotactic protein 1 and IL-8, an elevation of the global translation rate, an increase in the levels of several proteins critical for translation, and a reduction of microRNA-mediated translational repression. Moreover, using the BEAS-2B cell system and a mouse model, we found that RNA processing bodies (P bodies), cytoplasmic domains linked to storage and/or degradation of translationally silenced mRNAs, are significantly reduced in activated bronchial epithelial cells, suggesting a physiological role for P bodies in airway inflammation. Our study reveals an orchestrated change among posttranscriptional mechanisms, which help sustain high levels of inflammatory mediator production in bronchial epithelium during the pathogenesis of inflammatory airway diseases.

Myocardial expression profiles of candidate molecules in patients with arrhythmogenic right ventricular cardiomyopathy/dysplasia compared to those with dilated cardiomyopathy and healthy controls.

BACKGROUND Arrhythmogenic right ventricular cardiomyopathy/dysplasia (ARVC/D) is mainly an autosomal dominant disease characterized by fibrofatty infiltration of the right ventricle, leading to ventricular arrhythmias. Mutations in desmosomal proteins can be identified in about half of the patients. The pathogenic mechanisms leading to disease expression remain unclear. OBJECTIVE The purpose of this study was to investigate myocardial expression profiles of candidate molecules involved in the pathogenesis of ARVC/D. METHODS Myocardial messenger RNA (mRNA) expression of 62 junctional molecules, 5 cardiac ion channel molecules, 8 structural molecules, 4 apoptotic molecules, and 6 adipogenic molecules was studied. The averaged expression of candidate mRNAs was compared between ARVC/D samples (n = 10), nonfamilial dilated cardiomyopathy (DCM) samples (n = 10), and healthy control samples (n = 8). Immunohistochemistry and quantitative protein expression analysis were performed. Genetic analysis using next generation sequencing was performed in all patients with ARVC/D. RESULTS Following mRNA levels were significantly increased in patients with ARVC/D compared to those with DCM and healthy controls: phospholamban (P ≤ .001 vs DCM; P ≤ .001 vs controls), healthy tumor protein 53 apoptosis effector (P = .001 vs DCM; P ≤ .001 vs controls), and carnitine palmitoyltransferase 1β (P ≤ .001 vs DCM; P = 0.008 vs controls). Plakophillin-2 (PKP-2) mRNA was downregulated in patients with ARVC/D with PKP-2 mutations compared with patients with ARVC/D without PKP-2 mutations (P = .04). Immunohistochemistry revealed significantly increased protein expression of phospholamban, tumor protein 53 apoptosis effector, and carnitine palmitoyltransferase 1β in patients with ARVC/D and decreased PKP-2 expression in patients with ARVC/D carrying a PKP-2 mutation. CONCLUSION Changes in the expression profiles of sarcolemmal calcium channel regulation, apoptosis, and adipogenesis suggest that these molecular pathways may play a critical role in the pathogenesis of ARVC/D, independent of the underlying genetic mutations.

GLUCAGON GENE EXPRESSION: ANALYSIS OF PREPROGLUCAGON

Glucagon is a 29 amino acid polypeptide hormone produced in the (alpha) cells of the pancreatic islets. The purpose of this research was to understand better the role of glucagon in the regulation of metabolic processes. As with other polypeptide hormones, the synthesis of glucagon is thought to involve a larger precursor, which is then enzymatically cleaved to the functional form. The specific research objectives were to obtain cloned copies of the messenger RNA (mRNA) for pancreatic glucagon, to determine their primary sequences, and from this coding information to deduce the amino acid sequence of the initial glucagon precursor. From this suggested preproglucagon sequence and prior information on possible proglucagon intermediate processing products, the overall objective of this research is to propose a possible pathway for the biosynthesis of pancreatic glucagon.^ Synthetic oligodeoxynucleotide probes of 14-nucleotides (14-mer) and 17-nucleotides (a 17-mer) complementary to codons specifying a unique sequence of mature glucagon were synthesized. The ('32)P-labeled-14-mer was hybridized with size-fractionated fetal bovine pancreatic poly(A('+))RNA bound to nitrocellulose. RNA fractions of (TURN)14S were found to hybridize specifically, resulting in an (TURN)10-fold enrichment for these sequences. These poly(A('+))RNAs were translated in a cell-free system and the products analyzed by gel electrophoresis. The translation products were found to be enriched for a protein of the putative size of mammalian preproglucagon ((TURN)21 kd). These enriched RNA fractions were used to construct a complementary DNA (cDNA) library is plasmid pBR322.^ Screening of duplicate colony filters with the ('32)P-labeled-17-mer and a ('32)P-labeled-17-mer-primed cDNA probe indicated 25 possible glucagon clones from 3100 colonies screened. Restriction mapping of 6 of these clones suggested that they represented a single mRNA species. Primary sequence analysis of one clone containing a 1200 base pair DNA insert revealed that it contained essentially a full-length copy of glucagon cDNA.^ Analaysis of the cDNA suggested that it encoded an initial translation product of 180 amino acids with an M(,r) = 21 kd. The first initiation codon (ATG, methionine) followed by the longest open reading frame of 540 nucleotides was preceded by a 5'-untranslated region of 90 nucleotides, and was followed by a longer 3'-untranslated region of 471 nucleotides, resulting in a total of 1101 nucleotides. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Subnuclear localization of the active variant surface glycoprotein gene expression site in Trypanosoma brucei

In Trypanosoma brucei, transcription by RNA polymerase II and 5′ capping of messenger RNA are uncoupled: a capped spliced leader is trans spliced to every RNA. This decoupling makes it possible to have protein-coding gene transcription driven by RNA polymerase I. Indeed, indirect evidence suggests that the genes for the major surface glycoproteins, variant surface glycoproteins (VSGs) in bloodstream-form trypanosomes, are transcribed by RNA polymerase I. In a single trypanosome, only one VSG expression site is maximally transcribed at any one time, and it has been speculated that transcription takes place at a unique site within the nucleus, perhaps in the nucleolus. We tested this by using fluorescence in situ hybridization. With probes that cover about 50 kb of the active 221 expression site, we detected nuclear transcripts of this site in a single fluorescent spot, which did not colocalize with the nucleolus. Analysis of marker gene-tagged active expression site DNA by fluorescent DNA in situ hybridization confirmed the absence of association with the nucleolus. Even an active expression site in which the promoter had been replaced by an rDNA promoter did not colocalize with the nulceolus. As expected, marker genes inserted in the rDNA array predominantly colocalize with the nucleolus, whereas the tubulin gene arrays do not. We conclude that transcription of the active VSG expression site does not take place in the nucleolus.

Inhibition of luciferase expression in transgenic Aedes aegypti mosquitoes by Sindbis virus expression of antisense luciferase RNA

A rapid and reproducible method of inhibiting the expression of specific genes in mosquitoes should further our understanding of gene function and may lead to the identification of mosquito genes that determine vector competence or are involved in pathogen transmission. We hypothesized that the virus expression system based on the mosquito-borne Alphavirus, Sindbis (Togaviridae), may efficiently transcribe effector RNAs that inhibit expression of a targeted mosquito gene. To test this hypothesis, germ-line-transformed Aedes aegypti that express luciferase (LUC) from the mosquito Apyrase promoter were intrathoracically inoculated with a double subgenomic Sindbis (dsSIN) virus TE/3′2J/anti-luc (Anti-luc) that transcribes RNA complementary to the 5′ end of the LUC mRNA. LUC activity was monitored in mosquitoes infected with either Anti-luc or control dsSIN viruses expressing unrelated antisense RNAs. Mosquitoes infected with Anti-luc virus exhibited 90% reduction in LUC compared with uninfected and control dsSIN-infected mosquitoes at 5 and 9 days postinoculation. We demonstrate that a gene expressed from the mosquito genome can be inhibited by using an antisense strategy. The dsSIN antisense RNA expression system is an important tool for studying gene function in vivo.

A low threshold level of expression of mutant-template telomerase RNA inhibits human tumor cell proliferation

The ribonucleoprotein telomerase synthesizes telomeric DNA by copying an intrinsic RNA template. In most cancer cells, telomerase is highly activated. Here we report a telomerase-based antitumor strategy: expression of mutant-template telomerase RNAs in human cancer cells. We expressed mutant-template human telomerase RNAs in prostate (LNCaP) and breast (MCF-7) cancer cell lines. Even a low threshold level of expression of telomerase RNA gene constructs containing various mutant templates, but not the control wild-type template, decreased cellular viability and increased apoptosis. This occurred despite the retention of normal levels of the endogenous wild-type telomerase RNA and endogenous wild-type telomerase activity and unaltered stable telomere lengths. In vivo tumor xenografts of a breast cancer cell line expressing a mutant-template telomerase RNA also had decreased growth rates. Therefore, mutant-template telomerase RNAs exert a strongly dominant-negative effect on cell proliferation and tumor growth. These results support the potential use of mutant-template telomerase RNA expression as an antineoplastic strategy.

Regulation of MAPK activation, AP-1 transcription factor expression and keratinocyte differentiation in wounded fetal skin

Fetal epithelium retains the ability to re-epithelialize a wound in organotypic culture in a manner not dependent on the presence of underlying dermal substrata. This capacity is lost late in the third trimester of gestation or after embryonic day 17 (E-17) in the rat such that embryonic day 19 (E-19) wounds do not re-epithelialize. Moreover, wounds created in E-17 fetuses in utero heal in a regenerative, scar-free fashion. To investigate the molecular events regulating re-epithelialization in fetal skin, the wound-induced expression profile and tissue localization of activator protein 1 (AP-1) transcription factors c-Fos and c-Jun was characterised in E-17 and E-19 skin using organotypic fetal cultures. The involvement of mitogen-activated protein kinase (MAPK) signaling in mediating wound-induced transcription factor expression and wound re-epithelialization was assessed, with the effect of wounding on the expression of keratinocyte differentiation markers determined. Our results show that expression of AP-1 transcription factors was induced immediately by wounding and localized predominantly to the epidermis in E-17 and E-19 skin. c-fos and c-jun induction was transient in E-17 skin with MAPK-dependent c-fos expression necessary for the re-epithelialization of an excisional wound in organotypic culture. In E-19 skin, AP-11 expression persisted beyond 12 h post-wounding, and marked upregulation of the keratinocyte differentiation markers keratin 10 and loricrin was observed. No such changes in the expression of keratin 10 or loricrin occurred in E-17 skin. These findings indicate that re-epithelialization in fetal skin is regulated by wound-induced AP-1 transcription factor expression via MAPK and the differentiation status of keratinocytes.

Oxygen-regulated gene expression in bovine blastocysts

Oxygen concentrations used during in vitro embryo culture can influence embryo development, cell numbers, and gene expression. Here we propose that the preimplantation bovine embryo possesses a molecular mechanism for the detection of, and response to, oxygen, mediated by a family of basic helix-loop-helix transcription factors, the hypoxia-inducible factors (HIFs). Day 5 compacting bovine embryos were cultured under different oxygen tensions (2%, 7%, 20%) and the effect on the expression of oxygen-regulated genes, development, and cell number allocation and HIFalpha protein localization were examined. Bovine in vitro-produced embryos responded to variations in oxygen concentration by altering gene expression. GLUT1 expression was higher following 2% oxygen culture compared with 7% and 20% cultured blastocysts. HIF mRNA expression (HIF1alpha, HIF2alpha) was unaltered by oxygen concentration. HIF2alpha protein was predominantly localized to the nucleus of blastocysts. In contrast, HIF1alpha protein was undetectable at any oxygen concentration or in the presence of the HIF protein stabilizer desferrioxamine (DFO), despite being detectable in cumulus cells following normal maturation conditions, acute anoxic culture, or in the presence of DFO. Oxygen concentration also significantly altered inner cell mass cell proportions at the blastocyst stage. These results suggest that oxygen can influence gene expression in the bovine embryo during postcompaction development and that these effects may be mediated by HIF2alpha.