Comparative transcriptome analyses of Pseudomonas aeruginosa

One of the hallmarks of bacterial survival is their ability to adapt rapidly to changing environmental conditions. Niche adaptation is a response to the signals received that are relayed, often to regulators that modulate gene expression. In the post-genomic era, DNA microarrays are used to study the dynamics of gene expression on a global scale. Numerous studies have used Pseudomonas aeruginosa--a Gram-negative environmental and opportunistic human pathogenic bacterium--as the model organism in whole-genome transcriptome analysis. This paper reviews the transcriptome studies that have led to immense advances in our understanding of the biology of this intractable human pathogen. Comparative analysis of 23 P. aeruginosa transcriptome studies has led to the identification of a unique set of genes that are signal specific and a core set that is differentially regulated. The 303 genes in the core set are involved in bacterial homeostasis, making them attractive therapeutic targets.

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^

Avaliação funcional e estrutural de um novo peptídeo antimicrobiano do escorpião Tityus stigmurus

In Brazil, there is a high incidence of venomous animals. Among them, scorpions are highlighted by their medical importance, and for being their venom a source of several molecules with biological and pharmacological activity not yet fully understood, including several bioactive peptides. Antimicrobial peptides (AMPs) are components of the immune system in prokaryotes and eukaryotes, used in the first line of defense against microorganisms. In the present study, we characterized the first PAM previously identified through transcriptome of the venom gland of the scorpion Tityus stigmurus, named Stigmurin. The characteristics of Stigmurin were investigated by computational modeling and construction of dendrogram. In vitro tests investigated the antibacterial, antifungal, haemolytic and cytotoxic effects of crude venom and Stigmurin. In addition, the structural characteristics of Stigmurin were investigated by circular dochroism in water, 2, 2 , 2- trifluoethanol (TFE) and sodium dodecyl sulfate (SDS) and the models were refined by molecular dynamics simulations. The results showed that the selected sequence encodes a mature protein of 17 amino acid residues and the dendrogram reveals a case of convergent evolution. The crude venom showed no antimicrobial activity, however, Stigmurin exhibited a broad spectrum of antibacterial activity, with minimal inhibitory concentrations (MIC) ranging from 31.25 and 250 µg/mL for different strains, while the hemolytic activity at these concentrations was low. In cytotoxicity studies, the crude venom was unable to reduce cell viability in VERO E6 cells; in contrast, its activity in SiHa cells was significantly higher, corresponding to IC50 of 3.6 µg/mL. For Stigmurin the concentration sable to decrease cell viability of Vero E6 and SiHa cells in 50% were 275.67 µg/mL and 212.54 µg/mL, respectively. The dichroism spectra revealed the conformational flexibility, with predominating extended and β–sheet structures, as well as a remark able renaturation ability. The results suggest that Stigmurin could be considered as a potential antiinfective drug

Transcriptome analysis of the synganglion from the honey bee mite, Varroa destructor and RNAi knockdown of neural peptide targets

Date of Acceptance: 20/12/2015 This work was funded by BBSRC-LINK grant # BB/J01009X/1 and Vita Europe Ltd. We are grateful to the Scottish Beekeepers Association, especially Mr Phil McAnespie in supporting this work at its inception. We acknowledge partial funding from a Genesis Faraday SPARK Award, part of a Scottish Government SEEKIT project for the early part of this work. We are grateful to Prof David Evans for his advice on Varroa destructor viruses.

Single-Cell Transcriptome Analysis of Olfactory Sensory Neurons

Olfactory sensory neurons (OSNs), which detect a myriad of odorants, are known to express one allele of one olfactory receptor (OR) gene (Olfr) from the largest gene family in the mammalian genome. The OSNs expressing the same OR project their axons to the main olfactory bulb where they converge to form glomeruli. This “One neuron-one receptor rule” makes the olfactory epithelium (OE), which consists of a vast number of OSNs expressing unique ORs, one of the most heterogeneous cell populations. However, the mechanism of how the single OR allele is chosen remains unclear along with the question of whether one OSN only expresses a single OR gene, a hypothesis that has not been rigorously verified while we performed the experiments. Moreover, failure of axonal targeting to single glomerulus was observed in MeCP2 deficient OSNs where delayed development was proposed as an explanation for the phenotype. How Mecp2 mutation caused this aberrant targeting is not entirely understood.

In this dissertation, we explored the transcriptomes of single and mature OSNs by single-cell RNA-Seq to reveal their heterogeneity and further studied the OR gene expression from these isolated OSNs. The singularity of sequenced OSNs was ensured by the observation of monoallelic expression of X-linked genes from the hybrid samples from crosses between mice of different strains where strain-specific polymorphisms could be used to track the allelic origins of SNP-containing reads. The clustering of expression profiles from triplicates that originated from the same cell assured that the transcriptomic identities of OSNs were maintained through the experimental process. The average gene expression profiles of sequenced OSNs correlated well to the conventional transcriptome data of FACS-sorted Omp-positive cells, and the top-ranked expression of OR was conceded in the single-OSN transcriptomes. While exploring cellular diversity, in addition to OR genes, we revealed nearly 200 differentially expressed genes among the sequenced OSNs in this study. Among the 36 sequenced OSNs, eight cells (22.2%) showed multiple OR gene expression and the presences of additional ORs were not restricted to the neighbor loci that shared the transcriptional effect of the primary OR expression, suggesting that the “One neuron-one receptor rule” might not be strictly true at the transcription level. All of the inferable ORs, including additional co-expressed ORs, were shown to be monoallelic. Our sequencing of 21 Mecp2308 mutant OSNs, of which 62% expressed more than one OR genes, and the expression levels of the additional ORs were significantly higher than those in the wild-type, suggested that MeCP2 plays a role in the regulation of singular OR gene expression. Dual label in situ hybridization along with the sequence data revealed that dorsal and ventral ORs were co-expressed in the same Mecp2 mutant OSN, further implying that MeCP2 might be involved in regulation of OR territories in the OE. Our results suggested a new role of MeCP2 in OR gene choice and ratified that this multiple-OR expression caused by Mecp2 mutation did not accompany delayed OSN development that has been observed in the previous studies on the Mecp2 mutants.

Transcriptome and Functional Analysis of Carotid Body Glomus Cells

The carotid body (CB) is a major arterial chemoreceptor containing glomus cells that are activated by changes in arterial blood contents including oxygen. Despite significant advancement in the characterization of their physiological properties, our understanding on the underlying molecular machinery and signaling pathway in CB glomus cells is still limited.

To overcome these limitations, in chapter 1, I demonstrated the first transcriptome profile of CB glomus cells using single cell sequencing technology, which allowed us to uncover a set of abundantly expressed genes, including novel glomus cell-specific transcripts. These results revealed involvement of G protein-coupled receptor (GPCR) signaling pathway, various types of ion channels, as well as atypical mitochondrial subunits in CB function. I also identified ligands for the mostly highly expressed GPCR (Olfr78) in CB glomus cells and examined this receptor’s role in CB mediated hypoxic ventilatory response.

Current knowledge of CB suggest glomus cells rely on unusual mitochondria for their sensitivity to hypoxia. I previously identified the atypical mitochondrial subunit Ndufa4l2 as a highly over-represented gene in CB glomus cells. In chapter 2, to investigate the functional significance of Ndufa4l2 in CB function, I phenotyped both Ndufa4l2 knockout mice and mice with conditional Ndufa4l2 deletion in CB glomus cells. I found that Ndufa4l2 is essential to the establishment of regular breathing after birth. Ablating Ndufa4l2 in postnatal CB glomus cells resulted in defective CB sensitivity to hypoxia as well as CB mediated hypoxic ventilatory response. Together, our data showed that Ndufa4l2 is critical to respiratory control and the oxygen sensitivity of CB glomus cells.

Thyroid hormone and gland status of Greenland sledge dogs (Canis lupus familiaris) exposed to contaminated blubber

As a model of high trophic level carnivores, sledge dogs were fed from 2 to 18 months of age with minke whale blubber containing organohalogen compounds (OHC) corresponding to 128 µg PCB/day. Controls were fed uncontaminated porcine fat. Thyroid hormone levels were assessed in 7 exposed and 7 control sister bitches (sampled at age 6-18 months) and 4 exposed and 4 control pups, fed the same diet as their mothers (sampled age 3-12 months). Lower free and total T3 and T4 were seen in exposed vs. control bitches beyond 10 months of age, and total T3 was lower through 3-12 months of age in exposed pups. A negative correlation with thyroid gland weight was significant for SumDDT, as was a positive association with total T3 for dieldrin. This study therefore supports observational data that OHCs may adversely affect thyroid functions, and it suggests that OHC exposure duration of 10 months or more may be required for current OHC contamination levels to result in detectable adverse effects on thyroid hormone dynamics.

DNA in Amphibian and Reptile Venom Permits Access to Genomes Without Specimen Sacrifice

Amphibian defensive skin secretions and reptile venoms are rich sources of bioactive peptides with potential pharmacological/pharmaceutical applications. As amphibian and reptile populations are in rapid global decline, our research
group has been developing analytical methods that permit generation of robust molecular data from non-invasive skin secretion samples and venom samples. While previously we have demonstrated that parallel proteome and venom gland
transcriptome analyses can be performed on such samples, here we report the presence of DNA that facilitates the more widely-used applications of gene sequencing, such as molecular phylogenetics, in a non-invasive manner that circumvents specimen sacrifice. From this “surrogate” tissue, we acquired partial 12S and 16S rRNA gene sequences that are presented for illustration purposes. Thus from a single sample of amphibian skin secretion and reptile venom, robust and complementary proteome, transcriptome and genome data can be generated for applications in diverse scientific disciplines.

Comprehensive Transcriptome Analysis of Sex-Biased Expressed Genes Reveals Discrete Biological and Physiological Features of Male and Female Schistosoma japonicum

Schistosomiasis is a chronic and debilitating disease caused by blood flukes (digenetic trematodes) of the genus Schistosoma. Schistosomes are sexually dimorphic and exhibit dramatic morphological changes during a complex lifecycle which requires subtle gene regulatory mechanisms to fulfil these complex biological processes. In the current study, a 41,982 features custom DNA microarray, which represents the most comprehensive probe coverage for any schistosome transcriptome study, was designed based on public domain and local databases to explore differential gene expression in S. japonicum. We found that approximately 1/10 of the total annotated genes in the S. japonicum genome are differentially expressed between adult males and females. In general, genes associated with the cytoskeleton, and motor and neuronal activities were readily expressed in male adult worms, whereas genes involved in amino acid metabolism, nucleotide biosynthesis, gluconeogenesis, glycosylation, cell cycle processes, DNA synthesis and genome fidelity and stability were enriched in females. Further, miRNAs target sites within these gene sets were predicted, which provides a scenario whereby the miRNAs potentially regulate these sex-biased expressed genes. The study significantly expands the expressional and regulatory characteristics of gender-biased expressed genes in schistosomes with high accuracy. The data provide a better appreciation of the biological and physiological features of male and female schistosome parasites, which may lead to novel vaccine targets and the development of new therapeutic interventions.

Transcriptional profiling of the oesophageal gland region of male worms of Schistosoma mansoni

The intestinal tract of schistosomes opens at the mouth and leads into the foregut or oesophageal region that is lined with syncytium continuous with the apical cytoplasm of the tegument. The oesophagus is surrounded by a specialised gland, the oesophageal gland. This gland releases materials into the lumen of the oesophagus and the region is thought to initiate the lysis of erythrocytes and neutralisation of immune effectors of the host. The oesophageal region is present in the early invasive schistosomulum, a stage potentially targetable by anti-schistosome vaccines. We used a 44k oligonucleotide microarray to identify highly up-regulated genes in microdissected frozen sections of the oesophageal gland of male worms of S. mansoni. We show that 122 genes were up-regulated 2-fold or higher in the oesophageal gland compared with a whole male worm tissue control. The enriched genes included several associated with lipid metabolism and transmembrane transport as well as some micro-exon genes. Since the oesophageal gland is important in the initiation of digestion and the fact that it develops early after invasion of the mammalian host, further study of selected highly up-regulated functionally important genes in this tissue may reveal new anti-schistosome intervention targets for schistosomiasis control.

Sex steroid feedback regulation of pituitary gonadotropins during early secondary oocyte growth in coho salmon (Oncorhynchus kisutch) – a potential target of endocrine disruption.

Thesis (Ph.D.)--University of Washington, 2016-06

Étude et décontamination du transcriptome de novo du nématode doré Globodera rostochiensis

Le nématode doré, Globodera rostochiensis, est un nématode phytoparasite qui peut infecter des plantes agricoles telles la pomme de terre, la tomate et l’aubergine. En raison des pertes de rendement considérables associées à cet organisme, il est justifiable de quarantaine dans plusieurs pays, dont le Canada. Les kystes du nématode doré protègent les œufs qu’ils contiennent, leur permettant de survivre (en état de dormance) jusqu’à 20 ans dans le sol. L’éclosion des œufs n’aura lieu qu’en présence d’exsudats racinaires d’une plante hôte compatible à proximité. Malheureusement, très peu de connaissances sont disponibles sur les mécanismes moléculaires liés à cette étape-clé du cycle vital du nématode doré. Dans cet ouvrage, nous avons utilisé la technique RNA-seq pour séquencer tous les ARNm d’un échantillon de kystes du nématode doré afin d’assembler un transcriptome de novo (sans référence) et d’identifier des gènes jouant un rôle dans les mécanismes de survie et d’éclosion. Cette méthode nous a permis de constater que les processus d’éclosion et de parasitisme sont étroitement reliés. Plusieurs effecteurs impliqués dans le mouvement vers la plante hôte et la pénétration de la racine sont induits dès que le kyste est hydraté (avant même le déclenchement de l’éclosion). Avec l’aide du génome de référence du nématode doré, nous avons pu constater que la majorité des transcrits du transcriptome ne provenaient pas du nématode doré. En effet, les kystes échantillonnés au champ peuvent contenir des contaminants (bactéries, champignons, etc.) sur leur paroi et même à l’intérieur du kyste. Ces contaminants seront donc séquencés et assemblés avec le transcriptome de novo. Ces transcrits augmentent la taille du transcriptome et induisent des erreurs lors des analyses post-assemblages. Les méthodes de décontamination actuelles utilisent des alignements sur des bases de données d’organismes connus pour identifier ces séquences provenant de contaminants. Ces méthodes sont efficaces lorsque le ou les contaminants sont connus (possède un génome de référence) comme la contamination humaine. Par contre, lorsque le ou les contaminants sont inconnus, ces méthodes deviennent insuffisantes pour produire un transcriptome décontaminé de qualité. Nous avons donc conçu une méthode qui utilise un algorithme de regroupement hiérarchique des séquences. Cette méthode produit, de façon récursive, des sous-groupes de séquences homogènes en fonction des patrons fréquents présents dans les séquences. Une fois les groupes créés, ils sont étiquetés comme contaminants ou non en fonction des résultats d’alignements du sous-groupe. Les séquences ambiguës ayant aucun ou plusieurs alignements différents sont donc facilement classées en fonction de l’étiquette de leur groupe. Notre méthode a été efficace pour décontaminer le transcriptome du nématode doré ainsi que d’autres cas de contamination. Cette méthode fonctionne pour décontaminer un transcriptome, mais nous avons aussi démontré qu’elle a le potentiel de décontaminer de courtes séquences brutes. Décontaminer directement les séquences brutes serait la méthode de décontamination optimale, car elle minimiserait les erreurs d’assemblage.

Phylum-wide transcriptome analysis of oogenesis and early embryogenesis in selected nematode species

Oogenesis is a prerequisite for embryogenesis in Metazoa. During both biological processes important decisions must be made to form the embryo and hence ensure the next generation: (1) Maternal gene products (mRNAs, proteins and nutrients) must be supplied to the embryo. (2) Polarity must be established and axes must be specified. While incorporation of maternal gene products occurs during oogenesis, the time point of polarity establishment and axis specification varies among species, as it is accomplished either prior, during, or after fertilisation. But not only the time point when these events take place varies among species but also the underlying mechanisms by which they are triggered. For the nematode model Caenorhabditis elegans the underlying pathways and gene regulatory networks (GRNs) are well understood. It is known that there the sperm entry point initiates a primary polarity in the 1-celled egg and with it the establishment of the anteroposterior axis. However, studies of other nematodes demonstrated that polarity establishment can be independent of sperm entry (Goldstein et al., 1998; Lahl et al., 2006) and that cleavage patterns, symmetry formation and cell specification also differ from C. elegans. In contrast to the studied Chromadorea (more derived nematodes including C. elegans), embryos of some marine Enoplea (more basal representatives) even show no discernible early polarity and blastomeres can adopt variable cell fates (Voronov and Panchin 1998). The underlying pathways controlling the obviously variant embryonic processes in non-Caenorhabditis nematodes are essentially unknown. In this thesis I addressed this issue by performing a detailed unbiased comparative transcriptome analysis based on microarrays and RNA sequencing of selected developmental stages in a variety of nematodes from different phylogenetic branches with C. elegans as a reference system and a nematomorph as an outgroup representative. In addition, I made use of available genomic data to determine the presence or absence of genes for which no expression had been detected. In particular, I focussed on components of selected pathways or GRNs which are known to play essential roles during C. elegans development and/or other invertebrate or vertebrate model systems. Oogenesis must be regulated differently in non-Caenorhabditis nematodes, as crucial controlling components of Wnt and sex determination signaling are absent in these species. In this respect, I identified female-specific expression of potential polarity associated genes during gonad development and oogenesis in the Enoplean nematode Romanomermis culicivorax. I could show that known downstream components of the polarity complexes PAR-3/-6/PKC-3 and PAR-1/-2 are absent in non-Caenorhabditis species. Even PAR-2 as part of the polarity complex does not exist in these nematodes. Instead, transcriptomes of nematodes (including C. elegans), show expression of other polarity-associated complexes such as the Lgl (Lethal giant larvae) complex. This result could pose an alternative route for nematodes and nematomorphs to initiate polarity during early embryogenesis. I could show that crucial pathways of axis specification, such as Wnt and BMP are very different in C. elegans compared to other nematodes. In the former, Wnt signaling, for instance, is mediated by four paralogous beta-catenins, while other Chromadorea have fewer and Enoplea only one beta-catenin. The transcriptomes of R. culicivorax and the nematomorph show that regulators of BMP (e.g. Chordin), are specifically expressed during early embryogenesis only in Enoplea and the close outgroup of nematomorphs. In conclusion, my results demonstrate that the molecular machinery controlling oogenesis and embryogenesis in nematodes is unexpectedly variable and C. elegans cannot be taken as a general model for nematode development. Under this perspective, Enoplean nematodes show more similarities with outgroups than with C. elegans. It appears that certain pathway components were lost or gained during evolution and others adopted new functions. Based on my findings I can conjecture, which pathway components may be ancestral and which were newly acquired in the course of nematode evolution.