Massive transcriptional start site analysis of human genes in hypoxia cells

An improved method for the highly specific detection of transcription start sites

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

Genome-Wide Mapping of 5' Isoforms with 5'-Seq

The transcriptome is far more complex than previously assumed. Transcripts from the same gene can differ in terms of transcription start site, transcription end site, or pattern of splicing, and growing evidence supports the functional importance of these distinct transcript isoforms. Easily identifying these isoforms experimentally via library construction and high-throughput sequencing is crucial. Current library construction methods for identifying transcription start sites (5' transcript isoforms) involve large number of steps and (expensive) reagents, utilization of cDNA intermediates for adapter ligation, and are less suitable for studying low-abundance isoforms. Here, I describe a quick protocol for the generation of sequencing libraries to define capped 5' isoforms (5'-Seq) of various abundances in yeast and suggest a 5' isoform data analysis pipeline. The protocol relies on the utilization of a dephosphorylation-decapping method (oligo-capping) to generate a sequencing library from mRNA fragments and is a simplification of previously published 5' isoform protocols in terms of the handling steps, time, and cost. This method is exemplified using Saccharomyces cerevisiae mRNA, but it can be applied to various cellular conditions to study the effects of 5' transcript isoforms on transcriptional and/or translational regulation. © 2023 Wiley Periodicals LLC. Basic Protocol: Construction of a DNA sequencing library from capped 5' isoforms Support Protocol: Sequencing data analysis.

Transcriptional Regulation of Liver-Type OATP1B3 (Lt-OATP1B3) and Cancer-Type OATP1B3 (Ct-OATP1B3) Studied in Hepatocyte-Derived and Colon Cancer-Derived Cell Lines

Due to alternative splicing, the SLCO1B3 gene encodes two protein variants; the hepatic uptake transporter liver-type OATP1B3 (Lt-OATP1B3) and the cancer-type OATP1B3 (Ct-OATP1B3) expressed in several cancerous tissues. There is limited information about the cell type-specific transcriptional regulation of both variants and about transcription factors regulating this differential expression. Therefore, we cloned DNA fragments from the promoter regions of the Lt-SLCO1B3 and the Ct-SLCO1B3 gene and investigated their luciferase activity in hepatocellular and colorectal cancer cell lines. Both promoters showed differences in their luciferase activity depending on the used cell lines. We identified the first 100 bp upstream of the transcriptional start site as the core promoter region of the Ct-SLCO1B3 gene. In silico predicted binding sites for the transcription factors ZKSCAN3, SOX9 and HNF1α localized within these fragments were further analyzed. The mutagenesis of the ZKSCAN3 binding site reduced the luciferase activity of the Ct-SLCO1B3 reporter gene construct in the colorectal cancer cell lines DLD1 and T84 to 29.9% and 14.3%, respectively. In contrast, using the liver-derived Hep3B cells, 71.6% residual activity could be measured. This indicates that the transcription factors ZKSCAN3 and SOX9 are important for the cell type-specific transcriptional regulation of the Ct-SLCO1B3 gene.

mRNA ageing shapes the Cap2 methylome in mammalian mRNA

The mRNA cap structure is a major site of dynamic mRNA methylation. mRNA caps exist in either the Cap1 or Cap2 form, depending on the presence of 2'-O-methylation on the first transcribed nucleotide or both the first and second transcribed nucleotides, respectively1,2. However, the identity of Cap2-containing mRNAs and the function of Cap2 are unclear. Here we describe CLAM-Cap-seq, a method for transcriptome-wide mapping and quantification of Cap2. We find that unlike other epitranscriptomic modifications, Cap2 can occur on all mRNAs. Cap2 is formed through a slow continuous conversion of mRNAs from Cap1 to Cap2 as mRNAs age in the cytosol. As a result, Cap2 is enriched on long-lived mRNAs. Large increases in the abundance of Cap1 leads to activation of RIG-I, especially in conditions in which expression of RIG-I is increased. The methylation of Cap1 to Cap2 markedly reduces the ability of RNAs to bind to and activate RIG-I. The slow methylation rate of Cap2 allows Cap2 to accumulate on host mRNAs, yet ensures that low levels of Cap2 occur on newly expressed viral RNAs. Overall, these results reveal an immunostimulatory role for Cap1, and that Cap2 functions to reduce activation of the innate immune response.

TSS-seq of Toxoplasma gondii sporozoites revealed a novel motif in stage-specific promoters

Toxoplasma gondii is one of the most common zoonotic protozoan parasites. It has three major infectious stages: rapidly multiplying tachyzoites (Tz), slowly replicating bradyzoites (Bz) and a resting/free-living stage, sporozoites (Sz). The regulatory mechanisms governing stage-specific gene expression are not fully understood. Few transcriptional start sites (TSS) are known for Sz. In this study, we obtained TSS of Sz using an oligo-capping method and RNA-seq analysis. We identified 1,043,503 TSS in the Sz transcriptome. These defined 38,973 TSS clusters, of which, 11,925 were expressed in Sz and 1535 TSS differentially expressed in Sz. Based on these data, we defined promoter regions and novel sporozoite stage-specific motifs using MEME. TGTANNTACA was distributed around -55 to -75 regions from each TSS. Interestingly, the same motif was reported in another apicomplexan, Plasmodium berghei, as a cis-element of female-specific gametocyte genes, implying the presence of common regulatory machinery. Further comparative analysis should better define the distribution and function of these elements in other members of this important parasitic phylum.

Global approaches for profiling transcription initiation

Transcription start site (TSS) selection influences transcript stability and translation as well as protein sequence. Alternative TSS usage is pervasive in organismal development, is a major contributor to transcript isoform diversity in humans, and is frequently observed in human diseases including cancer. In this review, we discuss the breadth of techniques that have been used to globally profile TSSs and the resulting insights into gene regulation, as well as future prospects in this area of inquiry.

Identification of VEGFA-centric temporal hypoxia-responsive dynamic cardiopulmonary network biomarkers

AIMS

Hypoxia, a pathophysiological condition, is profound in several cardiopulmonary diseases (CPD). Every individual's lethality to a hypoxia state differs in terms of hypoxia exposure time, dosage units and dependent on the individual's genetic makeup. Most of the proposed markers for CPD were generally aim to distinguish disease samples from normal samples. Although, as per the 2018 GOLD guidelines, clinically useful biomarkers for several cardio pulmonary disease patients in stable condition have yet to be identified. We attempt to address these key issues through the identification of Dynamic Network Biomarkers (DNB) to detect hypoxia induced early warning signals of CPD before the catastrophic deterioration.

MATERIALS AND METHODS

The human microvascular endothelial tissues microarray datasets (GSE11341) of lung and cardiac expose to hypoxia (1% O₂) for 3, 24 and 48 h were retrieved from the public repository. The time dependent differentially expressed genes were subjected to tissue specificity and promoter analysis to filtrate the noise levels in the networks and to dissect the tissue specific hypoxia induced genes. These filtered out genes were used to construct the dynamic segmentation networks. The hypoxia induced dynamic differentially expressed genes were validated in the lung and heart tissues of male rats. These rats were exposed to hypobaric hypoxia (simulated altitude of 25,000 or PO₂ - 282 mm of Hg) progressively for 3, 24 and 48 h.

KEY FINDINGS

To identify the temporal key genes regulated in hypoxia, we ranked the dominant genes based on their consolidated topological features from tissue specific networks, time dependent networks and dynamic networks. Overall topological ranking described VEGFA as a single node dynamic hub and strongly communicated with tissue specific genes to carry forward their tissue specific information. We named this type of VEGFAcentric dynamic networks as "V-DNBs". As a proof of principle, our methodology helped us to identify the V-DNBs specific for lung and cardiac tissues namely V-DNB^L and V-DNB^C respectively.

SIGNIFICANCE

Our experimental studies identified VEGFA, SLC2A3, ADM and ENO2 as the minimum and sufficient candidates of V-DNB^L. The dynamic expression patterns could be readily exploited to capture the pre disease state of hypoxia induced pulmonary vascular remodelling. Whereas in V-DNB^C the minimum and sufficient candidates are VEGFA, SCL2A3, ADM, NDRG1, ENO2 and BHLHE40. The time dependent single node expansion indicates V-DNB^C could also be the pre disease state pathological hallmark for hypoxia-associated cardiovascular remodelling. The network cross-talk and expression pattern between V-DNB^L and V-DNB^C are completely distinct. On the other hand, the great clinical advantage of V-DNBs for pre disease predictions, a set of samples during the healthy condition should suffice. Future clinical studies might further shed light on the predictive power of V-DNBs as prognostic and diagnostic biomarkers for CPD.

Characterization of cancer omics and drug perturbations in panels of lung cancer cells

To better understand the disruptions of transcriptional regulations and gene expression in lung cancers, we constructed a multi-omics catalogue of the responses of lung cancer cells to a series of chemical compounds. We generated and analyzed 3,240 RNA-seq and 3,393 ATAC-seq libraries obtained from 23 cell lines treated with 95 well-annotated compounds. To demonstrate the power of the created multi-omics resource, we attempted to identify drugs that could induce the designated changes alone or in combination. The basal multi-omics information was first integrated into co-expression modules. Among these modules, we identified a stress response module that may be a promising drug intervention target, as new combinations of compounds that could be used to regulate this module and the consequent phenotypic appearance of cancer cells have been identified. We believe that the multi-omics profiles generated in this study and the strategy used to stratify them will lead to more rational and efficient development of anticancer drugs.

FFPEcap-seq: a method for sequencing capped RNAs in formalin-fixed paraffin-embedded samples

The majority of clinical cancer specimens are preserved as formalin-fixed paraffin-embedded (FFPE) samples. For clinical molecular tests to have wide-reaching impact, they must be applicable to FFPE material. Accurate quantitative measurements of RNA derived from FFPE specimens is challenging because of low yields and high amounts of degradation. Here, we present FFPEcap-seq, a method specifically designed for sequencing capped 5′ ends of RNA derived from FFPE samples. FFPEcap-seq combines enzymatic enrichment of 5′ capped RNAs with template switching to create sequencing libraries. We find that FFPEcap-seq can faithfully capture mRNA expression levels in FFPE specimens while also detecting enhancer RNAs that arise from distal regulatory regions. FFPEcap-seq is a fast and straightforward method for making high-quality 5′ end RNA-seq libraries from FFPE-derived RNA.

CAGEfightR: analysis of 5'-end data using R/Bioconductor

BACKGROUND

5'-end sequencing assays, and Cap Analysis of Gene Expression (CAGE) in particular, have been instrumental in studying transcriptional regulation. 5'-end methods provide genome-wide maps of transcription start sites (TSSs) with base pair resolution. Because active enhancers often feature bidirectional TSSs, such data can also be used to predict enhancer candidates. The current availability of mature and comprehensive computational tools for the analysis of 5'-end data is limited, preventing efficient analysis of new and existing 5'-end data.

RESULTS

We present CAGEfightR, a framework for analysis of CAGE and other 5'-end data implemented as an R/Bioconductor-package. CAGEfightR can import data from BigWig files and allows for fast and memory efficient prediction and analysis of TSSs and enhancers. Downstream analyses include quantification, normalization, annotation with transcript and gene models, TSS shape statistics, linking TSSs to enhancers via co-expression, identification of enhancer clusters, and genome-browser style visualization. While built to analyze CAGE data, we demonstrate the utility of CAGEfightR in analyzing nascent RNA 5'-data (PRO-Cap). CAGEfightR is implemented using standard Bioconductor classes, making it easy to learn, use and combine with other Bioconductor packages, for example popular differential expression tools such as limma, DESeq2 and edgeR.

CONCLUSIONS

CAGEfightR provides a single, scalable and easy-to-use framework for comprehensive downstream analysis of 5'-end data. CAGEfightR is designed to be interoperable with other Bioconductor packages, thereby unlocking hundreds of mature transcriptomic analysis tools for 5'-end data. CAGEfightR is freely available via Bioconductor: bioconductor.org/packages/CAGEfightR .

DBTSS/DBKERO for integrated analysis of transcriptional regulation

DBTSS (Database of Transcriptional Start Sites)/DBKERO (Database of Kashiwa Encyclopedia for human genome mutations in Regulatory regions and their Omics contexts) is the database originally initiated with the information of transcriptional start sites and their upstream transcriptional regulatory regions. In recent years, we updated the database to assist users to elucidate biological relevance of the human genome variations or somatic mutations in cancers which may affect the transcriptional regulation. In this update, we facilitate interpretations of disease associated genomic variation, using the Japanese population as a model case. We enriched the genomic variation dataset consisting of the 13,368 individuals collected for various genome-wide association studies and the reference epigenome information in the surrounding regions using a total of 455 epigenome datasets (four tissue types from 67 healthy individuals) collected for the International Human Epigenome Consortium (IHEC). The data directly obtained from the clinical samples was associated with that obtained from various model systems, such as the drug perturbation datasets using cultured cancer cells. Furthermore, we incorporated the results obtained using the newly developed analytical methods, Nanopore/10x Genomics long-read sequencing of the human genome and single cell analyses. The database is made publicly accessible at the URL (http://dbtss.hgc.jp/).

Using RAMPAGE to Identify and Annotate Promoters in Insect Genomes

Application of Transcription Start Site (TSS) profiling technologies, coupled with large-scale next-generation sequencing (NGS) has yielded valuable insights into the location, structure, and activity of promoters across diverse metazoan model systems. In insects, TSS profiling has been used to characterize the promoter architecture of Drosophila melanogaster (Hoskins et al., Genome Res 21(2):182-192, 2011) and subsequently was employed to reveal widespread transposon-driven alternative promoter usage in the fruit fly (Batut et al., Genome Res 23:169-180, 2012).In this chapter we discuss the computational analysis of the experimental data derived from one of TSS profiling methods, RAMPAGE (RNA Annotation and Mapping of Promoters for Analysis of Gene Expression) that can be used for the precise, quantitative identification of promoters in insect genomes. We demonstrate this using the software tools GoRAMPAGE (Brendel and Raborn, GoRAMPAGE-A workflow for promoter detection by 5'-read mapping. https://github.com/BrendelGroup/GoRAMPAGE , 2016) and TSRchitect (Raborn and Brendel, TSRchitect: promoter identification from large-scale TSS profiling data. R Bioconductor package version 1.8.0 [Online]. Available: http://bioconductor.org/packages/release/bioc/html/TSRchitect.html , 2017), providing detailed instructions with the aim of taking the user from raw reads to processed results.

Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

Intussusceptive Vascular Remodeling Precedes Pathological Neovascularization

Supplemental Digital Content is available in the text.

Objective—

Pathological neovascularization is crucial for progression and morbidity of serious diseases such as cancer, diabetic retinopathy, and age-related macular degeneration. While mechanisms of ongoing pathological neovascularization have been extensively studied, the initiating pathological vascular remodeling (PVR) events, which precede neovascularization remains poorly understood. Here, we identify novel molecular and cellular mechanisms of preneovascular PVR, by using the adult choriocapillaris as a model.

Approach and Results—

Using hypoxia or forced overexpression of VEGF (vascular endothelial growth factor) in the subretinal space to induce PVR in zebrafish and rats respectively, and by analyzing choriocapillaris membranes adjacent to choroidal neovascular lesions from age-related macular degeneration patients, we show that the choriocapillaris undergo robust induction of vascular intussusception and permeability at preneovascular stages of PVR. This PVR response included endothelial cell proliferation, formation of endothelial luminal processes, extensive vesiculation and thickening of the endothelium, degradation of collagen fibers, and splitting of existing extravascular columns. RNA-sequencing established a role for endothelial tight junction disruption, cytoskeletal remodeling, vesicle- and cilium biogenesis in this process. Mechanistically, using genetic gain- and loss-of-function zebrafish models and analysis of primary human choriocapillaris endothelial cells, we determined that HIF (hypoxia-induced factor)-1α-VEGF-A-VEGFR2 signaling was important for hypoxia-induced PVR.

Conclusions—

Our findings reveal that PVR involving intussusception and splitting of extravascular columns, endothelial proliferation, vesiculation, fenestration, and thickening is induced before neovascularization, suggesting that identifying and targeting these processes may prevent development of advanced neovascular disease in the future.

Highly efficient 5' capping of mitochondrial RNA with NAD+ and NADH by yeast and human mitochondrial RNA polymerase

Bacterial and eukaryotic nuclear RNA polymerases (RNAPs) cap RNA with the oxidized and reduced forms of the metabolic effector nicotinamide adenine dinucleotide, NAD⁺ and NADH, using NAD⁺ and NADH as non-canonical initiating nucleotides for transcription initiation. Here, we show that mitochondrial RNAPs (mtRNAPs) cap RNA with NAD⁺ and NADH, and do so more efficiently than nuclear RNAPs. Direct quantitation of NAD⁺- and NADH-capped RNA demonstrates remarkably high levels of capping in vivo: up to ~60% NAD⁺ and NADH capping of yeast mitochondrial transcripts, and up to ~15% NAD⁺ capping of human mitochondrial transcripts. The capping efficiency is determined by promoter sequence at, and upstream of, the transcription start site and, in yeast and human cells, by intracellular NAD⁺ and NADH levels. Our findings indicate mtRNAPs serve as both sensors and actuators in coupling cellular metabolism to mitochondrial transcriptional outputs, sensing NAD⁺ and NADH levels and adjusting transcriptional outputs accordingly.

PHLDA1, another PHLDA family protein that inhibits Akt

The PHLDA family (pleckstrin homology‐like domain family) of genes consists of 3 members: PHLDA1, 2, and 3. Both PHLDA3 and PHLDA2 are phosphatidylinositol (PIP) binding proteins and function as repressors of Akt. They have tumor suppressive functions, mainly through Akt inhibition. Several reports suggest that PHLDA1 also has a tumor suppressive function; however, the precise molecular functions of PHLDA1 remain to be elucidated. Through a comprehensive screen for p53 target genes, we identified PHLDA1 as a novel p53 target, and we show that PHLDA1 has the ability to repress Akt in a manner similar to that of PHLDA3 and PHLDA2. PHLDA1 has a so‐called split PH domain in which the PH domain is divided into an N‐terminal (β sheets 1‐3) and a C‐terminal (β sheets 4‐7 and an α‐helix) portions. We show that the PH domain of PHLDA1 is responsible for its localization to the plasma membrane and binding to phosphatidylinositol. We also show that the function of the PH domain is essential for Akt repression. In addition, PHLDA1 expression analysis suggests that PHLDA1 has a tumor suppressive function in breast and ovarian cancers.

Waves of chromatin modifications in mouse dendritic cells in response to LPS stimulation

BACKGROUND

The importance of transcription factors (TFs) and epigenetic modifications in the control of gene expression is widely accepted. However, causal relationships between changes in TF binding, histone modifications, and gene expression during the response to extracellular stimuli are not well understood. Here, we analyze the ordering of these events on a genome-wide scale in dendritic cells in response to lipopolysaccharide (LPS) stimulation.

RESULTS

Using a ChIP-seq time series dataset, we find that the LPS-induced accumulation of different histone modifications follows clearly distinct patterns. Increases in H3K4me3 appear to coincide with transcriptional activation. In contrast, H3K9K14ac accumulates early after stimulation, and H3K36me3 at later time points. Integrative analysis with TF binding data reveals potential links between TF activation and dynamics in histone modifications. Especially, LPS-induced increases in H3K9K14ac and H3K4me3 are associated with binding by STAT1/2 and were severely impaired in Stat1-/- cells.

CONCLUSIONS

While the timing of short-term changes of some histone modifications coincides with changes in transcriptional activity, this is not the case for others. In the latter case, dynamics in modifications more likely reflect strict regulation by stimulus-induced TFs and their interactions with chromatin modifiers.

Caution needs to be taken when assigning transcription start sites to ends of protein-coding genes: a rebuttal

Naturally occurring stress-induced transcriptional readthrough is a recently discovered phenomenon, in which stress conditions lead to dramatic induction of long transcripts as a result of transcription termination failure. In 2015, we reported the induction of such downstream of gene (DoG) containing transcripts upon osmotic stress in human cells, while others observed similar transcripts in virus-infected and cancer cells. Using the rigorous methodology Cap-Seq, we demonstrated that DoGs result from transcriptional readthrough, not de novo initiation. More recently, we presented a genome-wide comparison of NIH3T3 mouse cells subjected to osmotic, heat, and oxidative stress and concluded that massive induction of transcriptional readthrough is a hallmark of the mammalian stress response. In their recent letter, Huang and Liu in contrast claim that DoG transcripts result from novel transcription initiation near the ends of genes. Their conclusions rest on analyses of a publicly available transcription start site (TSS-Seq) dataset from unstressed NIH3T3 cells. Here, we present evidence that this dataset identifies not only true transcription start sites, TSSs, but also 5′-ends of numerous snoRNAs, which are generally processed from introns in mammalian cells. We show that failure to recognize these erroneous assignments in the TSS-Seq dataset, as well as ignoring published Cap-Seq data on TSS mapping during osmotic stress, have led to misinterpretation by Huang and Liu. We conclude that, contrary to the claims made by Huang and Liu, TSS-Seq reads near gene ends cannot explain the existence of DoGs, nor their stress-mediated induction. Rather it is, as we originally demonstrated, transcriptional readthrough that leads to the formation of DoGs.

Comprehensive comparative analysis of 5'-end RNA-sequencing methods

Specialized RNA-seq methods are required to identify the 5' ends of transcripts, which are critical for studies of gene regulation, but these methods have not been systematically benchmarked. We directly compared six such methods, including the performance of five methods on a single human cellular RNA sample and a new spike-in RNA assay that helps circumvent challenges resulting from uncertainties in annotation and RNA processing. We found that the 'cap analysis of gene expression' (CAGE) method performed best for mRNA and that most of its unannotated peaks were supported by evidence from other genomic methods. We applied CAGE to eight brain-related samples and determined sample-specific transcription start site (TSS) usage, as well as a transcriptome-wide shift in TSS usage between fetal and adult brain.

Identification of Minimal p53 Promoter Region Regulated by MALAT1 in Human Lung Adenocarcinoma Cells

The MALAT1 long noncoding RNA is strongly linked to cancer progression. Here we report a MALAT1 function in repressing the promoter of p53 (TP53) tumor suppressor gene. p21 and FAS, well-known p53 targets, were upregulated by MALAT1 knockdown in A549 human lung adenocarcinoma cells. We found that these upregulations were mediated by transcriptional activation of p53 through MALAT1 depletion. In addition, we identified a minimal MALAT1-responsive region in the P1 promoter of p53 gene. Flow cytometry analysis revealed that MALAT1-depleted cells exhibited G1 cell cycle arrest. These results suggest that MALAT1 affects the expression of p53 target genes through repressing p53 promoter activity, leading to influence the cell cycle progression.

Identification of potential regulatory mutations using multi-omics analysis and haplotyping of lung adenocarcinoma cell lines

The functional relevancy of mutations occurring in the regulatory regions in cancers remains mostly elusive. Here, we identified and analyzed regulatory mutations having transcriptional consequences in lung adenocarcinoma-derived cell lines. We phased the mutations in the regulatory regions to the downstream heterozygous SNPs in the coding regions and examined whether the ChIP-Seq variant tags of the regulatory SNVs and the RNA-Seq variant tags of their target transcripts showed biased frequency between the mutant and reference alleles. We identified 137 potential regulatory mutations affecting the transcriptional regulation of 146 RefSeq transcripts with at least 84 SNVs that create and/or disrupt potential transcription factor binding sites. For example, in the regulatory region of NFATC1 gene, a novel and active binding site for the ETS transcription factor family was created. Further examination revealed that 31 of these disruptions were presented in clinical lung adenocarcinoma samples and were associated with prognosis of patients.

Transcription start sites at the end of protein-coding genes

Previous studies demonstrated that massive induction of transcriptional readthrough generates downstream of gene-containing transcripts (DoGs) in cells under stress condition. Here, we analyzed TSS-seq (transcription start site sequencing) data from the DBTSS database. We investigated TSS tags at the end of gene for all pan-stress and untreated-cell DoGs, in comparison with expression-matched non-DoGs. We observed significantly more TSS tags at the end of pan-stress and untreated-cell DoG genes than non-DoG genes, even though their TSS tags in the promoter is the same. Importantly, the median value of TSS tags at gene end normalized to gene promoter is significantly higher than the median expression ratio of short DoG to host gene and of long DoG to host gene. Our results indicate that downstream overlapping long non-coding RNAs derived from the TSS at the gene end may be an important source of DoGs.

Light Controls Protein Localization through Phytochrome-Mediated Alternative Promoter Selection

Alternative promoter usage is a proteome-expanding mechanism that allows multiple pre-mRNAs to be transcribed from a single gene. The impact of this mechanism on the proteome and whether it is positively exploited in normal organismal responses remain unclear. We found that the plant photoreceptor phytochrome induces genome-wide changes in alternative promoter selection in Arabidopsis thaliana. Through this mechanism, protein isoforms with different N termini are produced that display light-dependent differences in localization. For instance, shade-grown plants accumulate a cytoplasmic isoform of glycerate kinase (GLYK), an essential photorespiration enzyme that was previously thought to localize exclusively to the chloroplast. Cytoplasmic GLYK constitutes a photorespiratory bypass that alleviates fluctuating light-induced photoinhibition. Therefore, phytochrome controls alternative promoter selection to modulate protein localization in response to changing light conditions. This study suggests that alternative promoter usage represents another ubiquitous layer of gene expression regulation in eukaryotes that contributes to diversification of the proteome.

Relatively frequent switching of transcription start sites during cerebellar development

BACKGROUND

Alternative transcription start site (TSS) usage plays important roles in transcriptional control of mammalian gene expression. The growing interest in alternative TSSs and their role in genome diversification spawned many single-gene studies on differential usages of tissue-specific or temporal-specific alternative TSSs. However, exploration of the switching usage of alternative TSS usage on a genomic level, especially in the central nervous system, is largely lacking.

RESULTS

In this study, We have prepared a unique set of time-course data for the developing cerebellum, as part of the FANTOM5 consortium ( http://fantom.gsc.riken.jp/5/ ) that uses their innovative capturing of 5' ends of all transcripts followed by Helicos next generation sequencing. We analyzed the usage of all transcription start sites (TSSs) at each time point during cerebellar development that provided information on multiple RNA isoforms that emerged from the same gene. We developed a mathematical method that systematically compares the expression of different TSSs of a gene to identify temporal crossover and non-crossover switching events. We identified 48,489 novel TSS switching events in 5433 genes during cerebellar development. This includes 9767 crossover TSS switching events in 1511 genes, where the dominant TSS shifts over time.

CONCLUSIONS

We observed a relatively high prevalence of TSS switching in cerebellar development where the resulting temporally-specific gene transcripts and protein products can play important regulatory and functional roles.

Modeling Keratoconus Using Induced Pluripotent Stem Cells

PURPOSE

To model keratoconus (KC) using induced pluripotent stem cells (iPSC) generated from fibroblasts of both KC and normal human corneal stroma by a viral method.

METHODS

Both normal and KC corneal fibroblasts from four human donors were reprogramed directly by delivering reprogramming factors in a single virus using 2A "self-cleaving" peptides, using a single polycistronic lentiviral vector coexpressing four transcription factors (Oct 4, Sox2, Klf4, and Myc) to yield iPSC. These iPS cells were characterized by immunofluorescence detection using of stem cell markers (SSEA4, Oct4, and Sox2). The mRNA sequencing was performed and the datasets were analyzed using ingenuity pathways analysis (IPA) software.

RESULTS

The generated stem cell-like clones expressed the pluripotency markers, SSEA4, Oct4, Sox2, Tra-1-60, and also expressed pax6. Our transcriptome analysis showed 4300 genes, which had 2-fold change and 870 genes with a q-value of <0.05 in keratoconus iPSC compared to normal iPSC. One of the genes that showed difference in KC iPSC was FGFR2 (down-regulated by 2.4 fold), an upstream target of Pi3-Kinase pathway, was further validated in keratoconus corneal sections and also KC iPSC-derived keratocytes (down regulated by 2.0-fold). Both normal and KC-derived keratocytes expressed keratocan, signature marker for keratocytes. KC iPSC-derived keratocytes showed adverse growth and proliferation and was further confirmed by using Ly2924002, a PI3k inhibitor, which severely affected the growth and differentiation in normal iPSC.

CONCLUSIONS

Based on our result, we propose a model for KC in which inhibition FGFR2-Pi3-Kinase pathway affects the AKT phosphorylation, and thus affecting the keratocytes survival signals. This inhibition of the survival signals could be a potential mechanism for the KC-specific decreased cell survival and apoptosis of keratocytes.

Identification of Arabidopsis genic and non-genic promoters by paired-end sequencing of TSS tags

Information about transcription start sites (TSSs) provides baseline data for the analysis of promoter architecture. In this paper we used paired- and single-end deep sequencing to analyze Arabidopsis TSS tags from several libraries prepared from roots, shoots, flowers and etiolated seedlings. The clustering of approximately 33 million mapped TSS tags led to the identification of 324 461 promoters that covered 79.7% (21 672/27 206) of protein-coding genes in the Arabidopsis genome. In addition we identified intragenic, antisense and orphan promoters that were not associated with any gene models. Of these, intragenic promoters exhibited unique characteristics regarding dinucleotide sequences at TSSs and core promoter element composition, suggesting that these promoters use different mechanisms of transcriptional initiation. An analysis of base composition with regard to promoter position revealed a low GC content throughout the promoter region and several local strand biases that were evident for TATA-type promoters, but not for Coreless-type promoters. Most observed strand biases coincided with strand biases of single nucleotide polymorphism rate. Our analysis also revealed that transcription of a gene is supported by an average of 2.7 genic promoters, among which one specific promoter, designated as a top promoter, substantially determines the expression level of the gene.

The Interaction Between Human Papillomaviruses and the Stromal Microenvironment

Human papillomaviruses (HPVs) are small, double-stranded DNA viruses that replicate in stratified squamous epithelia and cause a variety of malignancies. Current efforts in HPV biology are focused on understanding the virus-host interactions that enable HPV to persist for years or decades in the tissue. The importance of interactions between tumor cells and the stromal microenvironment has become increasingly apparent in recent years, but how stromal interactions impact the normal, benign life cycle of HPVs, or progression of lesions to cancer is less understood. Furthermore, how productively replicating HPV impacts cells in the stromal environment is also unclear. Here we bring together some of the relevant literature on keratinocyte-stromal interactions and their impacts on HPV biology, focusing on stromal fibroblasts, immune cells, and endothelial cells. We discuss how HPV oncogenes in infected cells manipulate other cells in their environment, and, conversely, how neighboring cells may impact the efficiency or course of HPV infection.

TSS seq based core promoter architecture in blood feeding Tsetse fly (Glossina morsitans morsitans) vector of Trypanosomiasis

Background

Transcription initiation regulation is mediated by sequence-specific interactions between DNA-binding proteins (transcription factors) and cis-elements, where BRE, TATA, INR, DPE and MTE motifs constitute canonical core motifs for basal transcription initiation of genes. Accurate identification of transcription start site (TSS) and their corresponding promoter regions is critical for delineation of these motifs. To this end, the genome scale analysis of core promoter architecture in insects has been confined to Drosophila. The recently sequenced Tsetse fly genome provides a unique opportunity to analyze transcription initiation regulation machinery in blood-feeding insects.

Results

A computational method for identification of TSS in newly sequenced Tsetse fly genome was evaluated, using TSS seq tags sampled from two developmental stages namely; larvae and pupae. There were 3134 tag clusters among which 45.4 % (1424) of the tag clusters mapped to first coding exons or their proximal predicted 5′UTR regions and 1.0 % (31) tag clusters mapping to transposons, within a threshold of 100 tags per cluster. These 1393 non transposon-derived core promoters had propensity for AT nucleotides. The −1/+1 and 1/+1 positions in D. melanogaster, and G. m. morsitans had propensity for CA and AA dinucleotides respectively. The 1393 tag clusters comprised narrow promoters (5 %), broad with peak promoters (23 %) and broad without peak promoters (72 %). Two-way motif co-occurrence analysis showed that the MTE-DPE pair is over-represented in broad core promoters. The frequently occurring triplet motifs in all promoter classes are the INR-MTE-DPE, TATA-MTE-DPE and TATA-INR-DPE. Promoters without the TATA motif had higher frequency of the MTE and INR motifs than those observed in Drosophila, where the DPE motif occur more frequently in promoters without TATA motif. Gene ontology terms associated with developmental processes were overrepresented in the narrow and broad with peak promoters.

Conclusions

The study has identified different motif combinations associated with broad promoters in a blood-feeding insect. In the case of TATA-less core promoters, G.m. morsitans uses the MTE to compensate for the lack of a TATA motif. The increasing availability of TSS seq data allows for revision of existing gene annotation datasets with the potential of identifying new transcriptional units.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1921-6) contains supplementary material, which is available to authorized users.

The genome-wide identification of promoter regions in Toxoplasma gondii

Parasites change their transcriptional systems in different developmental stages and in response to environmental changes. To investigate the molecular mechanisms that underlie transcriptional regulation, it is essential to identify the exact positions of the transcriptional start sites (TSSs) and characterize the upstream promoter regions. However, it has been essentially impossible to obtain comprehensive information using conventional methods. Here, we introduce our TSS-seq method, which combines full-length technology, oligo-capping, and rapidly developing next-generation sequencing technology. TSS-seq has enabled identification of TSS positions and upstream promoter activities as digital TSS tag counts within a reasonable cost and time frame. In this chapter, we describe in detail the TSS-seq method for the identification and characterization of the promoters in Toxoplasma gondii.

Whole transcriptome analysis with sequencing: methods, challenges and potential solutions

Whole transcriptome analysis plays an essential role in deciphering genome structure and function, identifying genetic networks underlying cellular, physiological, biochemical and biological systems and establishing molecular biomarkers that respond to diseases, pathogens and environmental challenges. Here, we review transcriptome analysis methods and technologies that have been used to conduct whole transcriptome shotgun sequencing or whole transcriptome tag/target sequencing analyses. We focus on how adaptors/linkers are added to both 5' and 3' ends of mRNA molecules for cloning or PCR amplification before sequencing. Challenges and potential solutions are also discussed. In brief, next generation sequencing platforms have accelerated releases of the large amounts of gene expression data. It is now time for the genome research community to assemble whole transcriptomes of all species and collect signature targets for each gene/transcript, and thus use known genes/transcripts to determine known transcriptomes directly in the near future.

From discovery to function: the expanding roles of long noncoding RNAs in physiology and disease

Long noncoding RNAs (lncRNAs) are a relatively poorly understood class of RNAs with little or no coding capacity transcribed from a set of incompletely annotated genes. They have received considerable attention in the past few years and are emerging as potentially important players in biological regulation. Here we discuss the evolving understanding of this new class of molecular regulators that has emerged from ongoing research, which continues to expand our databases of annotated lncRNAs and provide new insights into their physical properties, molecular mechanisms of action, and biological functions. We outline the current strategies and approaches that have been employed to identify and characterize lncRNAs, which have been instrumental in revealing their multifaceted roles ranging from cis- to trans-regulation of gene expression and from epigenetic modulation in the nucleus to posttranscriptional control in the cytoplasm. In addition, we highlight the molecular and biological functions of some of the best characterized lncRNAs in physiology and disease, especially those relevant to endocrinology, reproduction, metabolism, immunology, neurobiology, muscle biology, and cancer. Finally, we discuss the tremendous diagnostic and therapeutic potential of lncRNAs in cancer and other diseases.

Genetic basis of the highly efficient yeast Kluyveromyces marxianus: complete genome sequence and transcriptome analyses

BACKGROUND

High-temperature fermentation technology with thermotolerant microbes has been expected to reduce the cost of bioconversion of cellulosic biomass to fuels or chemicals. Thermotolerant Kluyveromyces marxianus possesses intrinsic abilities to ferment and assimilate a wide variety of substrates including xylose and to efficiently produce proteins. These capabilities have been found to exceed those of the traditional ethanol producer Saccharomyces cerevisiae or lignocellulose-bioconvertible ethanologenic Scheffersomyces stipitis.

RESULTS

The complete genome sequence of K. marxianus DMKU 3-1042 as one of the most thermotolerant strains in the same species has been determined. A comparison of its genomic information with those of other yeasts and transcriptome analysis revealed that the yeast bears beneficial properties of temperature resistance, wide-range bioconversion ability, and production of recombinant proteins. The transcriptome analysis clarified distinctive metabolic pathways under three different growth conditions, static culture, high temperature, and xylose medium, in comparison to the control condition of glucose medium under a shaking condition at 30°C. Interestingly, the yeast appears to overcome the issue of reactive oxygen species, which tend to accumulate under all three conditions.

CONCLUSIONS

This study reveals many gene resources for the ability to assimilate various sugars in addition to species-specific genes in K. marxianus, and the molecular basis of its attractive traits for industrial applications including high-temperature fermentation. Especially, the thermotolerance trait may be achieved by an integrated mechanism consisting of various strategies. Gene resources and transcriptome data of the yeast are particularly useful for fundamental and applied researches for innovative applications.

Dynamics of enhancers in myeloid antigen presenting cells upon LPS stimulation

BACKGROUND

Recent studies have underscored the role of enhancers in defining cell type-specific transcriptomes. Cell type-specific enhancers are bound by combinations of shared and cell type-specific transcription factors (TFs). However, little is known about combinatorial binding of TFs to enhancers, dynamics of TF binding following stimulation, or the downstream effects on gene expression. Here, we address these questions in two types of myeloid antigen presenting cells (APCs), macrophages and dendritic cells (DCs), before and after stimulation with lipopolysaccharide (LPS), a potent stimulator of the innate immune response.

RESULTS

We classified enhancers according to the combination of TFs binding them. There were significant correlations between the sets of TFs bound to enhancers prior to stimulation and expression changes of nearby genes after stimulation. Importantly, a set of enhancers pre-bound by PU.1, C/EBPβ, ATF3, IRF4, and JunB was strongly associated with induced genes and binding by stimulus-activated regulators. Our classification suggests that transient loss of ATF3 binding to a subset of these enhancers is important for regulation of early-induced genes. Changes in TF-enhancer binding after stimulation were correlated with binding by additional activated TFs and with the presence of proximally located enhancers.

CONCLUSIONS

The results presented in this study reveal the complexity and dynamics of TF- enhancer binding before and after stimulation in myeloid APCs.

Methylation-specific PCR: four steps in primer design

Methylation-specific PCR (MSP) is still the method of choice for a single gene methylation study. The proper design of the primer pairs is a prerequisite for obtaining reliable PCR results. Despite numerous protocols describing the rules for MSP primer design, none of them provide a comprehensive approach to the problem. Our aim was to depict a workflow for the primer design that is concise and easy to follow. In order to achieve this goal, adequate tools for promoter sequence retrieval, MSP primer design and subsequent in silico analysis are presented and discussed. Furthermore, a few instructive examples regarding a good versus a poor primer design are provided. Finally, primer design is demonstrated according to the proposed workflow. This article aims to provide researchers, interested in a single gene methylation studies, with useful information regarding successful primer design.

Evaluation of sgRNA target sites for CRISPR-mediated repression of TP53

The CRISPR (clustered regularly interspaced short palindromic repeats) platform has been developed as a general method to direct proteins of interest to gene targets. While the native CRISPR system delivers a nuclease that cleaves and potentially mutates target genes, researchers have recently employed catalytically inactive CRISPR-associated 9 nuclease (dCas9) in order to target and repress genes without DNA cleavage or mutagenesis. With the intent of improving repression efficiency in mammalian cells, researchers have also fused dCas9 with a KRAB repressor domain. Here, we evaluated different genomic sgRNA targeting sites for repression of TP53. The sites spanned a 200-kb distance, which included the promoter, transcript sequence, and regions flanking the endogenous human TP53 gene. We showed that repression up to 86% can be achieved with dCas9 alone (i.e., without use of the KRAB domain) by targeting the complex to sites near the TP53 transcriptional start site. This work demonstrates that efficient transcriptional repression of endogenous human genes can be achieved by the targeted delivery of dCas9. Yet, the efficiency of repression strongly depends on the choice of the sgRNA target site.

DBTSS as an integrative platform for transcriptome, epigenome and genome sequence variation data

DBTSS (http://dbtss.hgc.jp/) was originally constructed as a collection of uniquely determined transcriptional start sites (TSSs) in humans and some other species in 2002. Since then, it has been regularly updated and in recent updates epigenetic information has also been incorporated because such information is useful for characterizing the biological relevance of these TSSs/downstream genes. In the newest release, Release 9, we further integrated public and original single nucleotide variation (SNV) data into our database. For our original data, we generated SNV data from genomic analyses of various cancer types, including 97 lung adenocarcinomas and 57 lung small cell carcinomas from Japanese patients as well as 26 cell lines of lung cancer origin. In addition, we obtained publically available SNV data from other cancer types and germline variations in total of 11,322 individuals. With these updates, users can examine the association between sequence variation pattern in clinical lung cancers with its corresponding TSS-seq, RNA-seq, ChIP-seq and BS-seq data. Consequently, DBTSS is no longer a mere storage site for TSS information but has evolved into an integrative platform of a variety of genome activity data.

The Babesia bovis gene and promoter model: an update from full-length EST analysis

BACKGROUND

Babesia bovis is an apicomplexan parasite that causes babesiosis in infected cattle. Genomes of pathogens contain promising information that can facilitate the development of methods for controlling infections. Although the genome of B. bovis is publically available, annotated gene models are not highly reliable prior to experimental validation. Therefore, we validated a preproposed gene model of B. bovis and extended the associated annotations on the basis of experimentally obtained full-length expressed sequence tags (ESTs).

RESULTS

From in vitro cultured merozoites, 12,286 clones harboring full-length cDNAs were sequenced from both ends using the Sanger method, and 6,787 full-length cDNAs were assembled. These were then clustered, and a nonredundant referential data set of 2,115 full-length cDNA sequences was constructed. The comparison of the preproposed gene model with our data set identified 310 identical genes, 342 almost identical genes, 1,054 genes with potential structural inconsistencies, and 409 novel genes. The median length of 5' untranslated regions (UTRs) was 152 nt. Subsequently, we identified 4,086 transcription start sites (TSSs) and 2,023 transcriptionally active regions (TARs) by examining 5' ESTs. We identified ATGGGG and CCCCAT sites as consensus motifs in TARs that were distributed around -50 bp from TSSs. In addition, we found ACACA, TGTGT, and TATAT sites, which were distributed periodically around TSSs in cycles of approximately 150 bp. Moreover, related periodical distributions were not observed in mammalian promoter regions.

CONCLUSIONS

The observations in this study indicate the utility of integrated bioinformatics and experimental data for improving genome annotations. In particular, full-length cDNAs with one-base resolution for TSSs enabled the identification of consensus motifs in promoter sequences and demonstrated clear distributions of identified motifs. These observations allowed the illustration of a model promoter composition, which supports the differences in transcriptional regulation frameworks between apicomplexan parasites and mammals.

Discovery of transcription start sites in the Chinese hamster genome by next-generation RNA sequencing

Chinese hamster ovary (CHO) cell lines are one of the major production tools for monoclonal antibodies, recombinant proteins, and therapeutics. Although many efforts have significantly improved the availability of sequence information for CHO cells in the last years, forthcoming draft genomes still lack the information depth known from the mouse or human genomes. Many genes annotated for CHO cells and the Chinese hamster reference genome still are in silico predictions, only insufficiently verified by biological experiments. The correct annotation of transcription start sites (TSSs) is of special interest for CHO cells, as these directly define the location of the eukaryotic core promoter. Our study aims to elucidate these largely unexplored regions, trying to shed light on promoter landscapes in the Chinese hamster genome. Based on a 5' enriched dual library RNA sequencing approach 6547 TSSs were identified, of which over 90% were assigned to known genes. These TSSs were used to perform extensive promoter studies using a novel, modular bioinformatics pipeline, incorporating analyses of important regulatory elements of the eukaryotic core promoter on per-gene level and on genomic scale.

CCL2 enhances pluripotency of human induced pluripotent stem cells by activating hypoxia related genes

Standard culture of human induced pluripotent stem cells (hiPSCs) requires basic Fibroblast Growth Factor (bFGF) to maintain the pluripotent state, whereas hiPSC more closely resemble epiblast stem cells than true naïve state ES which requires LIF to maintain pluripotency. Here we show that chemokine (C-C motif) ligand 2 (CCL2) enhances the expression of pluripotent marker genes through the phosphorylation of the signal transducer and activator of transcription 3 (STAT3) protein. Moreover, comparison of transcriptomes between hiPSCs cultured with CCL2 versus with bFGF, we found that CCL2 activates hypoxia related genes, suggesting that CCL2 enhanced pluripotency by inducing a hypoxic-like response.Further, we show that hiPSCs cultured with CCL2 can differentiate at a higher efficiency than culturing withjust bFGF and we show CCL2 can be used in feeder-free conditions [corrected]. Taken together, our finding indicates the novel functions of CCL2 in enhancing its pluripotency in hiPSCs.

A comparison of the rest complex binding patterns in embryonic stem cells and epiblast stem cells

We detected and characterized the binding sites of the representative Rest complex components Rest, Sin3A, and Lsd1. We compared their binding patterns in mouse embryonic stem (ES) cells and epiblast stem (EpiS) cells. We found few Rest sites unique to the EpiS cells. The ES-unique site features were distinct from those of the common sites, namely, the signal intensities were weaker, and the characteristic gene function categories differed. Our analyses showed that the Rest binding sites do not always overlap with the Sin3A and Lsd1 binding sites. The Sin3A binding pattern differed remarkably between the ES and EpiS cells and was accompanied by significant changes in acetylated-histone patterns in the surrounding regions. A series of transcriptome analyses in the same cell types unexpectedly showed that the putative target gene transcript levels were not dramatically different despite dynamic changes in the Rest complex binding patterns and chromatin statuses, which suggests that Rest is not the sole determinant of repression at its targets. Nevertheless, we identified putative Rest targets with explicitly enhanced transcription upon Rest knock-down in 143 and 60 common and ES-unique Rest target genes, respectively. Among such sites, several genes are involved in ES cell proliferation. In addition, we also found that long, intergenic non-coding RNAs were apparent Rest targets and shared similar features with the protein-coding target genes. Interestingly, such non-coding target genes showed less conservation through evolution than protein-coding targets. As a result of differences in the components and targets of the Rest complex, its functional roles may differ in ES and EpiS cells.

GJC2 promoter mutations causing Pelizaeus-Merzbacher-like disease

OBJECTIVE

Pelizaeus-Merzbacher-like disease is a rare hypomyelinating leukodystrophy caused by autosomal recessive mutations in GJC2, encoding a gap junction protein essential for production of a mature myelin sheath. A previously identified GJC2 mutation (c.-167A>G) in the promoter region is hypothesized to disrupt a putative SOX10 binding site; however, the lack of additional mutations in this region and contradictory functional data have limited the interpretation of this variant.

METHODS

We describe two independent Pelizaeus-Merzbacher-like disease families with a novel promoter region mutation and updated in vitro functional assays.

RESULTS

A novel GJC2 mutation (c.-170A>G) in the promoter region was identified in Pelizaeus-Merzbacher-like disease patients. In vitro functional assays using human GJC2 promoter constructs demonstrated that this mutation and the previously described c.-167A>G mutation similarly diminished the transcriptional activity driven by SOX10 and the binding affinity for SOX10.

INTERPRETATION

These findings support the role of GJC2 promoter mutations in Pelizaeus-Merzbacher-like disease. GJC2 promoter region mutation screening should be included in the evaluation of patients with unexplained hypomyelinating leukodystrophies.

Innate immunity interactome dynamics

Innate immune response involves protein–protein interactions, deoxyribonucleic acid (DNA)–protein interactions and signaling cascades. So far, thousands of protein–protein interactions have been curated as a static interaction map. However, protein–protein interactions involved in innate immune response are dynamic. We recorded the dynamics in the interactome during innate immune response by combining gene expression data of lipopolysaccharide (LPS)-stimulated dendritic cells with protein–protein interactions data. We identified the differences in interactome during innate immune response by constructing differential networks and identifying protein modules, which were up-/down-regulated at each stage during the innate immune response. For each protein complex, we identified enriched biological processes and pathways. In addition, we identified core interactions that are conserved throughout the innate immune response and their enriched gene ontology terms and pathways. We defined two novel measures to assess the differences between network maps at different time points. We found that the protein interaction network at 1 hour after LPS stimulation has the highest interactions protein ratio, which indicates a role for proteins with large number of interactions in innate immune response. A pairwise differential matrix allows for the global visualization of the differences between different networks. We investigated the toll-like receptor subnetwork and found that S100A8 is down-regulated in dendritic cells after LPS stimulation. Identified protein complexes have a crucial role not only in innate immunity, but also in circadian rhythms, pathways involved in cancer, and p53 pathways. The study confirmed previous work that reported a strong correlation between cancer and immunity.

Giardia lamblia transcriptome analysis using TSS-Seq and RNA-Seq

Giardia lamblia is a protozoan parasite that is found worldwide and has both medical and veterinary importance. We applied the transcription start sequence (TSS-seq) and RNA sequence (RNA-seq) techniques to study the transcriptome of the assemblage A WB strain trophozoite. We identified 8000 transcription regions (TR) with significant transcription. Of these regions, 1881 TRs were more than 500 nucleotides upstream of an annotated ORF. Combining both techniques helped us to identify 24 ORFs that should be re-annotated and 60 new ORFs. From the 8000 TRs, we were able to identify an AT-rich consensus that includes the transcription initiation site. It is possible that transcription that was previously thought to be bidirectional is actually unidirectional.

Genes in sport and doping

Genes control biological processes such as muscle production of energy, mitochondria biogenesis, bone formation, erythropoiesis, angiogenesis, vasodilation, neurogenesis, etc. DNA profiling for athletes reveals genetic variations that may be associated with endurance ability, muscle performance and power exercise, tendon susceptibility to injuries and psychological aptitude. Already, over 200 genes relating to physical performance have been identified by several research groups. Athletes’ genotyping is developing as a tool for the formulation of personalized training and nutritional programmes to optimize sport training as well as for the prediction of exercise-related injuries. On the other hand, development of molecular technology and gene therapy creates a risk of non-therapeutic use of cells, genes and genetic elements to improve athletic performance. Therefore, the World Anti-Doping Agency decided to include prohibition of gene doping within their World Anti-Doping Code in 2003. In this review article, we will provide a current overview of genes for use in athletes’ genotyping and gene doping possibilities, including their development and detection techniques.