The airway epithelial cell plays a central role in coordinating the pulmonary response to injury and inflammation. Here, transforming growth factor-β (TGFβ) activates gene expression programs to induce stem cell-like properties, inhibit expression of differentiated epithelial adhesion proteins and express mesenchymal contractile proteins. This process is known as epithelial mesenchymal transition (EMT); although much is known about the role of EMT in cellular metastasis in an oncogene-transformed cell, less is known about Type II EMT, that occurring in normal epithelial cells. In this study, we applied next generation sequencing (RNA-Seq) in primary human airway epithelial cells to understand the gene program controlling Type II EMT and how cytokine-induced inflammation modifies it. Generalized linear modeling was performed on a two-factor RNA-Seq experiment of 6 treatments of telomerase immortalized human small airway epithelial cells (3 replicates). Using a stringent cut-off, we identified 3,478 differentially expressed genes (DEGs) in response to EMT. Unbiased transcription factor enrichment analysis identified three clusters of EMT regulators, one including SMADs/TP63 and another NF-κB/RelA. Surprisingly, we also observed 527 of the EMT DEGs were also regulated by the TNF-NF-κB/RelA pathway. This Type II EMT program was compared to Type III EMT in TGFβ stimulated A549 alveolar lung cancer cells, revealing significant functional differences. Moreover, we observe that Type II EMT modifies the outcome of the TNF program, reducing IFN signaling and enhancing integrin signaling. We confirmed experimentally that TGFβ-induced the NF-κB/RelA pathway by observing a 2-fold change in NF-κB/RelA nuclear translocation. A small molecule IKK inhibitor blocked TGFβ-induced core transcription factor (SNAIL1, ZEB1 and Twist1) and mesenchymal gene (FN1 and VIM) expression. These data indicate that NF-κB/RelA controls a SMAD-independent gene network whose regulation is required for initiation of Type II EMT. Type II EMT dramatically affects the induction and kinetics of TNF-dependent gene networks.

The G-quadruplex is a non-canonical DNA structure biologically significant in DNA replication, transcription and telomere stability. To date, only G4s with all guanines originating from the same strand of DNA have been considered in the context of the human nuclear genome. Here, I discuss interstrand topological configurations of G-quadruplex DNA, consisting of guanines from both strands of genomic DNA; an algorithm is presented for predicting such structures. I have identified over 550,000 non-overlapping interstrand G-quadruplex forming sequences in the human genome—significantly more than intrastrand configurations. Functional analysis of interstrand G-quadruplex sites shows strong association with transcription initiation, the results are consistent with the XPB and XPD transcriptional helicases binding only to G-quadruplex DNA with interstrand topology. Interstrand quadruplexes are also enriched in origin of replication sites. Several topology classes of interstrand quadruplex-forming sequences are possible, and different topologies are enriched in different types of structural elements. The list of interstrand quadruplex forming sequences, and the computer program used for their prediction are available at the web address http://moment.utmb.edu/allquads.

Modulators (Ms) are proteins that modify the activity of transcription factors (TFs) and influence expression of their target genes (TGs). To discover modulators of NF-κB/RelA, we first identified 365 NF-κB/RelA-binding proteins using liquid chromatography-tandem mass spectrometry (LC-MS/MS). We used a probabilistic model to infer 8349 (M, NF-κB/RelA, TG) triplets and their modes of modulatory action from our combined LC-MS/MS and ChIP-Seq (ChIP followed by next generation sequencing) data, published RelA modulators and TGs, and a compendium of gene expression profiles. Hierarchical clustering of the derived modulatory network revealed functional subnetworks and suggested new pathways modulating RelA transcriptional activity. The modulators with the highest number of TGs and most non-random distribution of action modes (measured by Shannon entropy) are consistent with published reports. Our results provide a repertoire of testable hypotheses for experimental validation. One of the NF-κB/RelA modulators we identified is STAT1. The inferred (STAT1, NF-κB/RelA, TG) triplets were validated by LC-selected reaction monitoring-MS and the results of STAT1 deletion in human fibrosarcoma cells. Overall, we have identified 562 NF-κB/RelA modulators, which are potential drug targets, and clarified mechanisms of achieving NF-κB/RelA multiple functions through modulators. Our approach can be readily applied to other TFs.

The spatial organization of eukaryotic genomes is thought to play an important role in regulating gene expression. The recent advances in experimental methods including chromatin capture techniques, as well as the large amounts of accumulated gene expression data allow studying the relationship between spatial organization of the genome and co-expression of protein-coding genes. To analyse this genome-wide relationship at a single gene resolution, we combined the interchromosomal DNA contacts in the yeast genome measured by Duan et al. with a comprehensive collection of 1,496 gene expression datasets. We find significant enhancement of co-expression among genes with contact links. The co-expression is most prominent when two gene loci fall within 1,000 base pairs from the observed contact. We also demonstrate an enrichment of inter-chromosomal links between functionally related genes, which suggests that the non random nature of the genome organization serves to facilitate coordinated transcription in groups of genes.

Exposure to environmental tobacco smoke (ETS) is a known risk factor for the development of chronic lung diseases, cancer, and the exacerbation of viral infections. Extracellular vesicles (EVs) have been identified as novel mediators of cell–cell communication through the release of biological content. Few studies have investigated the composition/function of EVs derived from human airway epithelial cells (AECs) exposed to cigarette smoke condensate (CSC), as surrogates for ETS. Using novel high-throughput technologies, we identified a diverse range of small noncoding RNAs (sncRNAs), including microRNA (miRNAs), Piwi-interacting RNA (piRNAs), and transfer RNA (tRNAs) in EVs from control and CSC-treated SAE cells. CSC treatment resulted in significant changes in the EV content of miRNAs. A total of 289 miRNAs were identified, with five being significantly upregulated and three downregulated in CSC EVs. A total of 62 piRNAs were also detected in our EV preparations, with five significantly downregulated and two upregulated in CSC EVs. We used TargetScan and Gene Ontology (GO) analysis to predict the biological targets of hsa-miR-3913-5p, the most represented miRNA in CSC EVs. Understanding fingerprint molecules in EVs will increase our knowledge of the relationship between ETS exposure and lung disease, and might identify potential molecular targets for future treatments.

Background: In vertebrate development, the segmental pattern of the body axis is established as somites, masses of mesoderm distributed along the two sides of the neural tube, are formed sequentially in the anterior-posterior axis. This mechanism depends on waves of gene expression associated with the Notch, Fgf and Wnt pathways. The underlying transcriptional regulation has been studied by whole-transcriptome mRNA profiling; however, interpretation of the results is limited by poor resolution, noisy data, small sample size and by the absence of a wall clock to assign exact time for recorded points.Results: We present a method of Maximum Entropy deconvolution in both space and time and apply it to extract, from microarray timecourse data, the full spatiotemporal expression profiles of genes involved in mouse somitogenesis. For regulated genes, we have reconstructed the temporal profiles and determined the timing of expression peaks along the somite cycle to a single-minute resolution. Our results also indicate the presence of a new class of genes (including Raf1 and Hes7) with two peaks of activity in two distinct phases of the somite cycle. We demonstrate that the timeline of gene expression precisely reflects their functions in the biochemical pathways and the direction of causation in the regulatory networks.Conclusions: By applying a novel framework for data analysis, we have shown a striking correspondence between gene expression times and their interactions and regulations during somitogenesis. These results prove the key role of finely tuned transcriptional regulation in the process. The presented method can be readily applied to studying somite formation in other datasets and species, and to other spatiotemporal processes.