Reconstruction of the regulatory network for Bacillus subtilis and reconciliation with gene expression data

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb. 2016.00275

mRNA metabolism: nonsense mediated mRNA decay as a tool for gene therapy and the role of human DIS3L2 in transcript degradation

Tese de mestrado em Biologia Humana e Ambiente, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2015

Changes in reproductive roles are associated with changes in gene expression in fire ant queens.

Abstract In species with social hierarchies, the death of dominant individuals typically upheaves the social hierarchy and provides an opportunity for subordinate individuals to become reproductives. Such a phenomenon occurs in the monogyne form of the fire ant, Solenopsis invicta, where colonies typically contain a single wingless reproductive queen, thousands of workers and hundreds of winged nonreproductive virgin queens. Upon the death of the mother queen, many virgin queens shed their wings and initiate reproductive development instead of departing on a mating flight. Workers progressively execute almost all of them over the following weeks. To identify the molecular changes that occur in virgin queens as they perceive the loss of their mother queen and begin to compete for reproductive dominance, we collected virgin queens before the loss of their mother queen, 6 h after orphaning and 24 h after orphaning. Their RNA was extracted and hybridized against microarrays to examine the expression levels of approximately 10 000 genes. We identified 297 genes that were consistently differentially expressed after orphaning. These include genes that are putatively involved in the signalling and onset of reproductive development, as well as genes underlying major physiological changes in the young queens.

IFN stimulated gene expression in the liver is a better predictor of treatment response in chronic hepatitis c than the IL28b (IFN lambda 3) genotype

Background: Therapy of chronic hepatitis C (CHC) with pegIFNa/ribavirin achieves sustained virologic response (SVR) in ~55%. Pre-activation of the endogenous interferon system in the liver is associated non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. Methods: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. Results: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. Conclusions: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Bacterial community structure of a pesticide-contaminated site and assessment of changes induced in community structure during bioremediation.

The introduction of culture-independent molecular screening techniques, especially based on 16S rRNA gene sequences, has allowed microbiologists to examine a facet of microbial diversity not necessarily reflected by the results of culturing studies. The bacterial community structure was studied for a pesticide-contaminated site that was subsequently remediated using an efficient degradative strain Arthrobacter protophormiae RKJ100. The efficiency of the bioremediation process was assessed by monitoring the depletion of the pollutant, and the effect of addition of an exogenous strain on the existing soil community structure was determined using molecular techniques. The 16S rRNA gene pool amplified from the soil metagenome was cloned and restriction fragment length polymorphism studies revealed 46 different phylotypes on the basis of similar banding patterns. Sequencing of representative clones of each phylotype showed that the community structure of the pesticide-contaminated soil was mainly constituted by Proteobacteria and Actinomycetes. Terminal restriction fragment length polymorphism analysis showed only nonsignificant changes in community structure during the process of bioremediation. Immobilized cells of strain RKJ100 enhanced pollutant degradation but seemed to have no detectable effects on the existing bacterial community structure.

Gene duplication and neofunctionalization: POLR3G and POLR3GL.

RNA polymerase III (Pol III) occurs in two versions, one containing the POLR3G subunit and the other the closely related POLR3GL subunit. It is not clear whether these two Pol III forms have the same function, in particular whether they recognize the same target genes. We show that the POLR3G and POLR3GL genes arose from a DNA-based gene duplication, probably in a common ancestor of vertebrates. POLR3G- as well as POLR3GL-containing Pol III are present in cultured cell lines and in normal mouse liver, although the relative amounts of the two forms vary, with the POLR3G-containing Pol III relatively more abundant in dividing cells. Genome-wide chromatin immunoprecipitations followed by high-throughput sequencing (ChIP-seq) reveal that both forms of Pol III occupy the same target genes, in very constant proportions within one cell line, suggesting that the two forms of Pol III have a similar function with regard to specificity for target genes. In contrast, the POLR3G promoter-not the POLR3GL promoter-binds the transcription factor MYC, as do all other promoters of genes encoding Pol III subunits. Thus, the POLR3G/POLR3GL duplication did not lead to neo-functionalization of the gene product (at least with regard to target gene specificity) but rather to neo-functionalization of the transcription units, which acquired different mechanisms of regulation, thus likely affording greater regulation potential to the cell.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

CYP3A activity measured by the midazolam test is not related to 3435 C >T polymorphism in the multiple drug resistance transporter gene.

A recent study with 69 Japanese liver transplants treated with tacrolimus found that the MDR13435 C >T polymorphism, but not the MDR12677 G >T polymorphism, was associated with differences in the intestinal expression level of CYP3A4 mRNA. In the present study, over 6 h, we measured the kinetics of a 75 microg oral dose of midazolam, a CYP3A substrate, in 21 healthy subjects genotyped for the MDR13435 C >T and 2677 G >T polymorphism. No statistically significant differences were found in the calculated pharmacokinetic parameters between the three 3435 C >T genotypes (TT, CT and CC group, respectively: Cmax (mean +/- SD: 0.30 +/- 0.08 ng/ml, 0.31 +/- 0.09 ng/ml and 0.31 +/- 0.11 ng/ml; Apparent clearance: 122 +/- 29 l/h, 156 +/- 92 l/h and 111 +/- 35 l/h; t1/2: 1.9 +/- 1.1 h, 1.6 +/- 0.90 h and 1.7 +/- 0.7 h). In addition, the 30-min 1'OH midazolam to midazolam ratio, a marker of CYP3A activity, determined in 74 HIV-positive patients before the introduction of antiretroviral treatment, was not significantly different between the three 3435 C >T genotypes (mean ratio +/- SD: 3.65 +/- 2.24, 4.22 +/- 3.49 and 4.24 +/- 2.03, in the TT, CT and CC groups, respectively). Similarly, no association was found between the MDR12677 G >T polymorphism and CYP3A activity in the healthy subjects or in the HIV-positive patients. The existence of a strong association between the activity of CYP3A and MDR13435 C >T and 2677 G >T polymorphisms appears unlikely, at least in Caucasian populations and/or in the absence of specific environmental factors.

Burkholderia sartisoli sp. nov., isolated from a polycyclic aromatic hydrocarbon-contaminated soil.

A Gram-negative, rod-shaped, aerobic bacterium, designated strain RP007(T), was isolated from a polycyclic aromatic hydrocarbon-contaminated soil in New Zealand. Two additional strains were recovered from a compost heap in Belgium (LMG 18808) and from the rhizosphere of maize in the Netherlands (LMG 24204). The three strains had virtually identical 16S rRNA gene sequences and whole-cell protein profiles, and they were identified as members of the genus Burkholderia, with Burkholderia phenazinium as their closest relative. Strain RP007(T) had a DNA G+C content of 63.5 mol% and could be distinguished from B. phenazinium based on a range of biochemical characteristics. Strain RP007(T) showed levels of DNA-DNA relatedness towards the type strain of B. phenazinium and those of other recognized Burkholderia species of less than 30 %. The results of 16S rRNA gene sequence analysis, DNA-DNA hybridization experiments and physiological and biochemical tests allowed the differentiation of strain RP007(T) from all recognized species of the genus Burkholderia. Strains RP007(T), LMG 18808 and LMG 24204 are therefore considered to represent a single novel species of the genus Burkholderia, for which the name Burkholderia sartisoli sp. nov. is proposed. The type strain is RP007(T) (=LMG 24000(T) =CCUG 53604(T) =ICMP 13529(T)).

Characterization of pip5k3 fleck corneal dystrophy-linked gene in zebrafish.

PIKfyve is a kinase encoded by pip5k3 involved in phosphatidylinositols (PdtIns) pathways. These lipids building cell membranes have structural functions and are involved in complex intracellular regulations. Mutations in human PIP5K3 are associated with François-Neetens mouchetée fleck corneal dystrophy [Li, S., Tiab, L., Jiao, X., Munier, F.L., Zografos, L., Frueh, B.E., Sergeev, Y., Smith, J., Rubin, B., Meallet, M.A., Forster, R.K., Hejtmancik, J.F., Schorderet, D.F., 2005. Mutations in PIP5K3 are associated with François-Neetens mouchetee fleck corneal dystrophy. Am. J. Hum. Genet. 77, 54-63]. We cloned the zebrafish pip5k3 and report its molecular characterization and expression pattern in adult fish as well as during development. The zebrafish PIKfyve was 70% similar to the human homologue. The gene encompassed 42 exons and presented four alternatively spliced variants. It had a widespread expression in the adult organs and was localized in specific cell types in the eye as the cornea, lens, ganglion cell layer, inner nuclear layer and outer limiting membrane. Pip5k3 transcripts were detected in early cleavage stage embryos. Then it was uniformly expressed at 10 somites, 18 somites and 24 hpf. Its expression was then restricted to the head region at 48 hpf, 72 hpf and 5 dpf and partial expression was found in somites at 72 hpf and 5 dpf. In situ on eye sections at 3 dpf showed a staining mainly in lens, outer limiting membrane, inner nuclear layer and ganglion cell layer. A similar expression pattern was found in the eye at 5 dpf. A temporal regulation of the spliced variants was observed at 1, 3 and 5 dpf and they were also found in the adult eye.

Transcriptional regulation of CD4 gene expression by T cell factor-1/beta-catenin pathway.

By interacting with MHC class II molecules, CD4 facilitates lineage development as well as activation of Th cells. Expression of physiological levels of CD4 requires a proximal CD4 enhancer to stimulate basic CD4 promoter activity. T cell factor (TCF)-1/beta-catenin pathway has previously been shown to regulate thymocyte survival via up-regulating antiapoptotic molecule Bcl-xL. By both loss and gain of function studies, in this study we show additional function of TCF-1/beta-catenin pathway in the regulation of CD4 expression in vivo. Mice deficient in TCF-1 displayed significantly reduced protein and mRNA levels of CD4 in CD4+ CD8+ double-positive (DP) thymocytes. A transgene encoding Bcl-2 restored survival but not CD4 levels of TCF-1(-/-) DP cells. Thus, TCF-1-regulated survival and CD4 expression are two separate events. In contrast, CD4 levels were restored on DP TCF-1(-/-) cells by transgenic expression of a wild-type TCF-1, but not a truncated TCF-1 that lacks a domain required for interacting with beta-catenin. Furthermore, forced expression of a stabilized beta-catenin, a coactivator of TCF-1, resulted in up-regulation of CD4. TCF-1 or stabilized beta-catenin greatly stimulated activity of a CD4 reporter gene driven by a basic CD4 promoter and the CD4 enhancer. However, mutation of a potential TCF binding site located within the enhancer abrogated TCF-1 and beta-catenin-mediated activation of CD4 reporter. Finally, recruitment of TCF-1 to CD4 enhancer was detected in wild-type but not TCF-1 null mice by chromatin-immunoprecipitation analysis. Thus, our results demonstrated that TCF/beta-catenin pathway enhances CD4 expression in vivo by recruiting TCF-1 to stimulate CD4 enhancer activity.

The p53 gene expression and its developmental regulation in schistosomes

We have studied the gene expression, especially of the oncoproteins, and its regulation in schistosomes. Schistosomes have a complex life cycle with defined dimorphic lifestyle. The parasite are so far unique in biology in expressing oncogene products in their adult stage. In order to characterize the expression and developmental regulation, a lambda gt 11 cDNA library and lambda EMBL4 genomic DNA library of each growth stage of Schistosoma mansoni and S. japonicum was constructed, and was screened with various monoclonal antibodies against ongogene products. One positive plaque reacted to anti-p53 antibody (Ab-2, Oncogene Science, Inc.) was further analyzed. This fusion protein was about 120 KDa in molecular weights, and expressed as 1.4 Kb RNA in the adult stage. P53 gene is well-known as the negative regulator of the cell cicle, and the mutations in the gene are turning out to be the most common genetic alterations in human cancers. The comparison of the gene structure among species and stages were being conducted. Chromosome structures, C-band formation, and the results of in situ hybridization using the phage probe would be discussed.

Identification of two steroid-responsive promoters of different strength controlled by the same estrogen-responsive element in the 5'-end region of the Xenopus laevis vitellogenin gene A1.

A structural and functional analysis of the 5'-end region of the Xenopus laevis vitellogenin gene A1 revealed two transcription initiation sites located 1.8 kilobases apart. A RNA polymerase II binding assay indicates that both promoters form initiation complexes efficiently. In vitro, using a transcription assay derived from a HeLa whole-cell extract, the upstream promoter is more than 10-fold stronger than the downstream one. In contrast, both promoters have a similar strength in a HeLa nuclear extract. In vivo, that is in estrogen-stimulated hepatocytes, it is the downstream promoter homologous to the one used by the other members of the vitellogenin gene family, which is 50-fold stronger than the upstream promoter. Thus, if functional vitellogenin mRNA results from this latter activity, it would contribute less than 1% to the synthesis of vitellogenin by fully induced Xenopus hepatocytes expressing the four vitellogenin genes. In contrast, both gene A1 promoters are silent in uninduced hepatocytes. Transfection experiments using the Xenopus cell line B3.2 in which estrogen-responsiveness has been introduced reveal that the strong downstream promoter is controlled by an estrogen responsive element (ERE) located 330 bp upstream of it. The upstream promoter can also be controlled by the same ERE. Since the region comprising the upstream promoter is flanked by a 200 base pair long inverted repeat with stretches of homology to other regions of the X. laevis genome, we speculate that it might have been inserted upstream of the vitellogenin gene A1 by a recombination event and consequently brought under control of the ERE lying 1.5 kilobases downstream.

Manipulating keratinocyte stem cell proliferation and differentiation using RNA interference

ABSTRACT : The epidermis, the outermost compartment of the skin, is a stratified and squamous epithelium that constantly self-renews. Keratinocytes, which represent the main epidermal population, are responsible for its cohesion and barrier function. Epidermal renewal necessitates a fine equilibrium between keratinocyte proliferation and differentiation. The keratinocyte stem cell, located in the basal cell layer, is responsible for epidermal homeostasis and regeneration during the wound healing process. The transcription factor p63 structurally belongs to the p53 superfamily. It is expressed in the basal and supra-basal cell layers of stratified epithelia and is thought to be important for the renewal or the differentiation of keratinocyte stem cells (Yang et al., 1999; Mills et al., 1999). In order to better understand its function, we established an in vitro model of p63 deficient human keratinocyte stem cells using a shp63 mediated RNA interference. Knockdown of endogenous p63 induces downregulation of cell-adhesion genes as previously described (Carroll et al., 2006). Interestingly, the replating of attached p63-knockdown keratinocytes on a feeder layer results in a loss of attachment and proliferation. They are no longer clonogenic. However, if the same population are replated in a fibrin matrix, extended fibrinolysis is reported, a common process in wound healing, suggesting that p63 regulates the fibrinolytic pathway. This result was confirmed by Q-PCR and shows that the urokinase pathway, which mediates fibrinolysis, is upregulated. Altogether, these findings suggest a mechanism in which the fine tuning of p63 expression promotes attachment or release of the keratinocyte stem cell from the basement membrane by inducing genes of adhesion and/or of fibrinolysis. This mechanism may be important for epidermal self-renewal, differentiation as well as wound healing. Its misregulation may be partly responsible for the p63 knockout phenotype. The downregulation of p63 also induces a decrease in LEKTI expression. LEKTI (lymphoepithelial Kazal-type serine protease inhibitor) is a serine protease inhibitor encoded by the Spink5 gene. It is expressed and secreted in the uppermost differentiated layers of stratified epithelia and plays a role in the desquamation process. When this gene is disrupted, humans develop the Netherton syndrome (Chavanas et al., 2000b). It is a dermatosis characterized by hair dysplasias, ichtyosiform erythroderma and impairment in epidermal barrier function promoting inflammation similarly as in psoriasis with inflammatory infiltrate in excess. TNFα (tumor necrosis factor alpha) and EDA1 (ectodysplasin A1) are two transmembraneprecursors that belong to the TNF superfamily, which is involved in immune and inflammation regulation (Smahi et al., 2002). We suggest that the secreted serine protease inhibitor LEKTI plays a role in the regulation of TNFα and EDA1 precursor cleavage and absence of LEKTI induces excess of inflammation. To investigate this hypothesis, we induced downregulation of Spink5 expression in rat keratinocyte stem cells by using a shSpink5 mediated RNA interference approach. Interestingly, expression of TNFα and EDA1 is modified after knockdown of Spink5 by Q-PCR. Moreover, downregulation of Spink5 induces loss of cohesiveness between keratinocytes and colonies adopt a scattered phenotype. Altogether, these preliminary data suggest that downregulation of LEKTI may play a role in the inflammatory response in Netherton syndrome patients, by regulating TNFα expression.