Coevolution combined with molecular dynamics simulations provides structural and mechanistic insights into the interactions between the integrator complex subunits.

Finding the 3D structure of large, multi-subunit complexes is difficult, despite recent advances in cryo-EM technology, due to remaining challenges to expressing and purifying subunits. Computational approaches that predict protein-protein interactions, including Direct Coupling Analysis (DCA), represent an attractive alternative for dissecting interactions within protein complexes. However, they are readily applicable only to small proteins due to high computational complexity and a high number of false positives. To solve this problem, we proposed a modified DCA approach, a powerful tool to predict the most likely interfaces of protein complexes. Since our modified approach cannot provide structural and mechanistic details of interacting peptides, we combine it with Molecular Dynamics (MD) simulations. To illustrate this novel approach, we predict interacting domains and structural details of interactions of two Integrator complex subunits, INTS9 and INTS11. Our predictions of interacting residues of INTS9/INTS11 are highly consistent with crystallographic structure. We then expand our procedure to two complexes whose structures are not well-studied: 1) The heterodimer formed by the Cleavage and Polyadenylation Specificity Factor 100-kD (CPSF100) and 73-kD (CPSF73); 2) The heterotrimer formed by INTS4/INTS9/INTS11. Experimental data supports our predictions of interactions within these two complexes, demonstrating that combining DCA and MD simulations is a powerful approach to revealing structural insights of large protein complexes.

Metabolomic and Lipidomic Signature of Skeletal Muscle with Constitutively Active Mechanistic Target of Rapamycin Complex 1.

BACKGROUND: Regulation of mechanistic target of rapamycin complex 1 (mTORC1) plays an important role in aging and nutrition. For example, caloric restriction reduces mTORC1 signaling and extends lifespan, whereas nutrient abundance and obesity increase mTORC1 signaling and reduce lifespan. Skeletal muscle-specific knockout (KO) of DEP domain-containing 5 protein (DEPDC5) results in constitutively active mTORC1 signaling, muscle hypertrophy and an increase in mitochondrial respiratory capacity. The metabolic profile of skeletal muscle, in the setting of hyperactive mTORC1 signaling, is not well known. OBJECTIVES: To determine the metabolomic and lipidomic signature in skeletal muscle from female and male wild-type (WT) and DEPDC5 KO mice. METHODS: Tibialis anterior (TA) muscles from WT and transgenic (conditional skeletal muscle-specific DEPDC5 KO) were obtained from female and male adult mice. Polar metabolites and lipids were extracted using a Bligh-Dyer extraction from 5 samples per group and identified and quantified by LC-MS/MS. Resulting analyte peak areas were analyzed with t-test, analysis of variance, and Volcano plots for group comparisons (e.g., WT compared with KO) and multivariate statistical analysis for genotype and sex comparisons. RESULTS: A total of 162 polar metabolites (organic acids, amino acids, and amines and acyl carnitines) and 1141 lipid metabolites were detected in TA samples by LC-MS/MS. Few polar metabolites showed significant differences in KO muscles compared with WT within the same sex group. P-aminobenzoic acid, ß-alanine, and dopamine were significantly higher in KO male muscle whereas erythrose-4-phosphate and oxoglutaric acid were significantly reduced in KO females. The lipidomic profile of the KO groups revealed an increase of muscle phospholipids and reduced triacylglycerol and diacylglycerol compared with the WT groups. CONCLUSIONS: Sex differences were detected in polar metabolome and lipids were dependent on genotype. The metabolomic profile of mice with hyperactive skeletal muscle mTORC1 is consistent with an upregulation of mitochondrial function and amino acid utilization for protein synthesis.

Inferring Gene Regulatory Networks from RNA-seq Data Using Kernel Classification.

Gene expression profiling is one of the most recognized techniques for inferring gene regulators and their potential targets in gene regulatory networks (GRN). The purpose of this study is to build a regulatory network for the budding yeast Saccharomyces cerevisiae genome by incorporating the use of RNA-seq and microarray data represented by a wide range of experimental conditions. We introduce a pipeline for data analysis, data preparation, and training models. Several kernel classification models; including one-class, two-class, and rare event classification methods, are used to categorize genes. We test the impact of the normalization techniques on the overall performance of RNA-seq. Our findings provide new insights into the interactions between genes in the yeast regulatory network. The conclusions of our study have significant importance since they highlight the effectiveness of classification and its contribution towards enhancing the present comprehension of the yeast regulatory network. When assessed, our pipeline demonstrates strong performance across different statistical metrics, such as a 99% recall rate and a 98% AUC score.

Diverse Clinical Manifestations and Challenges of Mucormycosis: Insights From Serial Cases.

Mucormycosis is a severe and potentially life-threatening infection caused by a group of fungi classified as mucormycetes within the scientific order Mucorales. These infections are characterized by rapid and invasive fungal growth, presenting significant treatment challenges. Here we present 5 cases encountered from 2018 to 2022 at the University of Texas Medical Branch in Galveston, Texas, including a novel Apophysomyces species. These cases illustrate the diverse clinical manifestations of mucormycosis, including pulmonary, rhino-cerebral, gastrointestinal, and soft tissue involvement. Our investigation incorporates information provided by a multidisciplinary team of clinical collaborators, emphasizing the findings from radiology, histopathology, and microbiology. Given the escalating global incidence of mucormycosis, it is crucial for clinicians to become familiar with associated clinical findings, comorbidities, and risk factors to facilitate prompt recognition, appropriate diagnostic testing, and timely initiation of treatment.

Respiratory syncytial virus infection changes the piwi-interacting RNA content of airway epithelial cells.

Piwi-interacting RNAs (piRNAs) are small non-coding RNAs (sncRNAs) of about 26-32 nucleotides in length and represent the largest class of sncRNA molecules expressed in animal cells. piRNAs have been shown to play a crucial role to safeguard the genome, maintaining genome complexity and integrity, as they suppress the insertional mutations caused by transposable elements. However, there is growing evidence for the role of piRNAs in controlling gene expression in somatic cells as well. Little is known about changes in piRNA expression and possible function occurring in response to viral infections. In this study, we investigated the piRNA expression profile, using a human piRNA microarray, in human small airway epithelial (SAE) cells infected with respiratory syncytial virus (RSV), a leading cause of acute respiratory tract infections in children. We found a time-dependent increase in piRNAs differentially expressed in RSV-infected SAE cells. We validated the top piRNAs upregulated and downregulated at 24 h post-infection by RT-qPCR and identified potential targets. We then used Gene Ontology (GO) tool to predict the biological processes of the predicted targets of the most represented piRNAs in infected cells over the time course of RSV infection. We found that the most significant groups of targets of regulated piRNAs are related to cytoskeletal or Golgi organization and nucleic acid/nucleotide binding at 15 and 24 h p.i. To identify common patterns of time-dependent responses to infection, we clustered the significantly regulated expression profiles. Each of the clusters of temporal profiles have a distinct set of potential targets of the piRNAs in the cluster Understanding changes in piRNA expression in RSV-infected airway epithelial cells will increase our knowledge of the piRNA role in viral infection and might identify novel therapeutic targets for viral lung-mediated diseases.

Strength of association between body mass index and physical function scores in paediatric burn patients: A National Institute on Disability, Independent Living, and Rehabilitation Research Burn Model System study.

OBJECTIVE: Increased body weight has been associated with reduced muscle wasting in the early catabolic phase after a severe burn. Yet, overweight and obese non-burn children often exhibit impaired musculoskeletal function, which may lead to poor physical function (PF). We aimed to determine the association between body mass index (BMI) at discharge and self-reported PF and caregiver proxy-reported PF during recovery of burned children. MATERIALS AND METHODS: This is a retrospective multisite longitudinal study in paediatric burn patients ((8-17 y old at time of burn). PF outcome measures were self-reported mobility, proxy-reported mobility, and upper extremity PF evaluated using PROMIS measures at 6-, 12-, and 24-months after injury. Primary exposure variable was BMI-for-age at discharge. RESULTS: A total of 118 paediatric patients, aged 11.7 ± 3.3 y, with burns covering 37.6 ± 18.8% of their total body surface area (TBSA) and BMI-for-age of 23.1 ± 5.4 kg/m2 at discharge were analyzed. BMI at discharge was not significantly associated with self-reported mobility scores 6 months after burn (beta coefficient =-0.23, p = 0.31), had a positive effect on mobility at 12 months (beta = 0.46, p = 0.05), and no effect at 24 months after injury (beta=-0.10, p = 0.60), when adjusted for burn size. BMI did not have a significant effect on proxy-reported mobility or upper extremity PF. CONCLUSION: A greater BMI at discharge was associated with improved self-reported PF at 12 months after burn but not at 6 months or 24 months, which suggests a faster recovery of PF in paediatric patients of larger body weight. Our data suggests that a larger body weight does not compromise the recovery of PF after burn.

Why a Constant Number of Vertebrae? Digital Control of Segmental Identity during Vertebrate Development: The Somite Cycle Controls a Digital, Chromatin-Based Counter That Defines Segmental Identity and Body Plans in Vertebrate Animals.

It is not understood how the numbers and identities of vertebrae are controlled during mammalian development. The remarkable robustness and conservation of segmental numbers may suggest the digital nature of the underlying process. The study proposes a mechanism that allows cells to obtain and store the segmental information in digital form, and to produce a pattern of chromatin accessibility that in turn regulates Hox gene expression specific to the metameric segment. The model requires that a regulatory element be present such that the number of occurrences of the motif between two consecutive Hox genes equals the number of segments under the control of the anterior gene. This is true for the recently discovered hydroxyl radical cleavage 3bp-periodic (HRC3) motif, associated with histone modifications and developmental genes. The finding not only allows the correct prediction of the numbers of segments using only sequence information, but also resolves the 40-year-old enigma of the function of temporal and spatial collinearity of Hox genes. The logic of the mechanism is illustrated in the attached animated video. How different aspects of the proposed mechanism can be tested experimentally is also discussed.

Inferring Causation in Yeast Gene Association Networks With Kernel Logistic Regression.

Computational prediction of gene-gene associations is one of the productive directions in the study of bioinformatics. Many tools are developed to infer the relation between genes using different biological data sources. The association of a pair of genes deduced from the analysis of biological data becomes meaningful when it reflects the directionality and the type of reaction between genes. In this work, we follow another method to construct a causal gene co-expression network while identifying transcription factors in each pair of genes using microarray expression data. We adopt a machine learning technique based on a logistic regression model to tackle the sparsity of the network and to improve the quality of the prediction accuracy. The proposed system classifies each pair of genes into either connected or nonconnected class using the data of the correlation between these genes in the whole Saccharomyces cerevisiae genome. The accuracy of the classification model in predicting related genes was evaluated using several data sets for the yeast regulatory network. Our system achieves high performance in terms of several statistical measures.

Cigarette Smoke Condensate Exposure Changes RNA Content of Extracellular Vesicles Released from Small Airway Epithelial Cells.

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.

Maximum parsimony interpretation of chromatin capture experiments.

We present a new approach to characterizing the global geometric state of chromatin from HiC data. Chromatin conformation capture techniques (3C, and its variants: 4C, 5C, HiC, etc.) probe the spatial structure of the genome by identifying physical contacts between genomic loci within the nuclear space. In whole-genome conformation capture (HiC) experiments, the signal can be interpreted as spatial proximity between genomic loci and physical distances can be estimated from the data. However, observed spatial proximity signal does not directly translate into persistent contacts within the nuclear space. Attempts to infer a single conformation of the genome within the nuclear space lead to internal geometric inconsistencies, notoriously violating the triangle inequality. These inconsistencies have been attributed to the stochastic nature of chromatin conformation or to experimental artifacts. Here we demonstrate that it can be explained by a mixture of cells, each in one of only several conformational states, contained in the sample. We have developed and implemented a graph-theoretic approach that identifies the properties of such postulated subpopulations. We show that the geometrical conflicts in a standard yeast HiC dataset, can be explained by only a small number of homogeneous populations of cells (4 populations are sufficient to reconcile 95,000 most prominent impossible triangles, 8 populations can explain 375,000 top geometric conflicts). Finally, we analyze the functional annotations of genes differentially interacting between the populations, suggesting that each inferred subpopulation may be involved in a functionally different transcriptional program.

Coevolution of Residues Provides Evidence of a Functional Heterodimer of 5-HT2AR and 5-HT2CR Involving Both Intracellular and Extracellular Domains.

Serotonin is a neurotransmitter that plays a role in regulating activities such as sleep, appetite, mood and substance abuse disorders; serotonin receptors 5-HT2AR and 5-HT2CR are active within pathways associated with substance abuse. It has been suggested that 5-HT2AR and 5-HT2CR may form a dimer that affects behavioral processes. Here we study the coevolution of residues in 5-HT2AR and 5-HT2CR to identify potential interactions between residues in both proteins. Coevolution studies can detect protein interactions, and since the thus uncovered interactions are subject to evolutionary pressure, they are likely functional. We assessed the significance of the 5-HT2AR/5-HT2CR interactions using randomized phylogenetic trees and found the coevolution significant (p-value = 0.01). We also discuss how co-expression of the receptors suggests the predicted interaction is functional. Finally, we analyze how several single nucleotide polymorphisms for the 5-HT2AR and 5-HT2CR genes affect their interaction. Our findings are the first to characterize the binding interface of 5-HT2AR/5-HT2CR and indicate a correlation between this interface and location of SNPs in both proteins.

i-BLESS is an ultra-sensitive method for detection of DNA double-strand breaks.

Maintenance of genome stability is a key issue for cell fate that could be compromised by chromosome deletions and translocations caused by DNA double-strand breaks (DSBs). Thus development of precise and sensitive tools for DSBs labeling is of great importance for understanding mechanisms of DSB formation, their sensing and repair. Until now there has been no high resolution and specific DSB detection technique that would be applicable to any cells regardless of their size. Here, we present i-BLESS, a universal method for direct genome-wide DNA double-strand break labeling in cells immobilized in agarose beads. i-BLESS has three key advantages: it is the only unbiased method applicable to yeast, achieves a sensitivity of one break at a given position in 100,000 cells, and eliminates background noise while still allowing for fixation of samples. The method allows detection of ultra-rare breaks such as those forming spontaneously at G-quadruplexes.

The NFκB subunit RELA is a master transcriptional regulator of the committed epithelial-mesenchymal transition in airway epithelial cells.

The epithelial-mesenchymal transition (EMT) is a multistep dedifferentiation program important in tissue repair. Here, we examined the role of the transcriptional regulator NF-κB in EMT of primary human small airway epithelial cells (hSAECs). Surprisingly, transforming growth factor ß (TGFß) activated NF-κB/RELA proto-oncogene, NF-κB subunit (RELA) translocation within 1 day of stimulation, yet induction of its downstream gene regulatory network occurred only after 3 days. A time course of TGFß-induced EMT transition was analyzed by RNA-Seq in the absence or presence of inducible shRNA-mediated silencing of RELA. In WT cells, TGFß stimulation significantly affected the expression of 2,441 genes. Gene set enrichment analysis identified WNT, cadherin, and NF-κB signaling as the most prominent TGFß-inducible pathways. By comparison, RELA controlled expression of 3,138 overlapping genes mapping to WNT, cadherin, and chemokine signaling pathways. Conducting upstream regulator analysis, we found that RELA controls six clusters of upstream transcription factors, many of which overlapped with a transcription factor topology map of EMT developed earlier. RELA triggered expression of three key EMT pathways: 1) the WNT/ß-catenin morphogen pathway, 2) the JUN transcription factor, and 3) the Snail family transcriptional repressor 1 (SNAI1). RELA binding to target genes was confirmed by ChIP. Experiments independently validating WNT dependence on RELA were performed by silencing RELA via genome editing and indicated that TGFß-induced WNT5B expression and downstream activation of the WNT target AXIN2 are RELA-dependent. We conclude that RELA is a master transcriptional regulator of EMT upstream of WNT morphogen, JUN, SNAI1-ZEB1, and interleukin-6 autocrine loops.

Predicting proteome dynamics using gene expression data.

While protein concentrations are physiologically most relevant, measuring them globally is challenging. mRNA levels are easier to measure genome-wide and hence are typically used to infer the corresponding protein abundances. The steady-state condition (assumption that protein levels remain constant) has typically been used to calculate protein concentrations, as it is mathematically convenient, even though it is often not satisfied. Here, we propose a method to estimate genome-wide protein abundances without this assumption. Instead, we assume that the system returns to its baseline at the end of the experiment, which is true for cyclic phenomena (e.g. cell cycle) and many time-course experiments. Our approach only requires availability of gene expression and protein half-life data. As proof-of-concept, we predicted proteome dynamics associated with the budding yeast cell cycle, the results are available for browsing online at http://dynprot.cent.uw.edu.pl/ . The approach was validated experimentally by verifying that the predicted protein concentration changes were consistent with measurements for all proteins tested. Additionally, if proteomic data are available as well, we can also infer changes in protein half-lives in response to posttranslational regulation, as we did for Clb2, a post-translationally regulated protein. The predicted changes in Clb2 abundance are consistent with earlier observations.

Respiratory Syncytial Virus Infection Changes Cargo Composition of Exosome Released from Airway Epithelial Cells.

Exosomes are microvesicles known to carry biologically active molecules, including RNA, DNA and proteins. Viral infections can induce profound changes in exosome composition, and exosomes have been implicated in viral transmission and pathogenesis. No information is current available regarding exosome composition and function during infection with Respiratory Syncytial Virus (RSV), the most important cause of lower respiratory tract infections in children. In this study, we characterized exosomes released from RSV-infected lung carcinoma-derived A549 cells. RNA deep sequencing revealed that RSV exosomes contain a diverse range of RNA species like messenger and ribosomal RNA fragments, as well as small noncoding RNAs, in a proportion different from exosomes isolated from mock-infected cells. We observed that both RNA and protein signatures of RSV were present in exosomes, however, they were not able to establish productive infection in uninfected cells. Exosomes isolated from RSV-infected cells were able to activate innate immune response by inducing cytokine and chemokine release from human monocytes and airway epithelial cells. These data suggest that exosomes may play an important role in pathogenesis or protection against disease, therefore understating their role in RSV infection may open new avenues for target identification and development of novel therapeutics.

A Conserved Structural Signature of the Homeobox Coding DNA in HOX genes.

The homeobox encodes a DNA-binding domain found in transcription factors regulating key developmental processes. The most notable examples of homeobox containing genes are the Hox genes, arranged on chromosomes in the same order as their expression domains along the body axis. The mechanisms responsible for the synchronous regulation of Hox genes and the molecular function of their colinearity remain unknown. Here we report the discovery of a conserved structural signature of the 180-base pair DNA fragment comprising the homeobox. We demonstrate that the homeobox DNA has a characteristic 3-base-pair periodicity in the hydroxyl radical cleavage pattern. This periodic pattern is significant in most of the 39 mammalian Hox genes and in other homeobox-containing transcription factors. The signature is present in segmented bilaterian animals as evolutionarily distant as humans and flies. It remains conserved despite the fact that it would be disrupted by synonymous mutations, which raises the possibility of evolutionary selective pressure acting on the structure of the coding DNA. The homeobox coding DNA may therefore have a secondary function, possibly as a regulatory element. The existence of such element may have important consequences for understanding how these genes are regulated.

Comparison between Timelines of Transcriptional Regulation in Mammals, Birds, and Teleost Fish Somitogenesis.

Metameric segmentation of the vertebrate body is established during somitogenesis, when a cyclic spatial pattern of gene expression is created within the mesoderm of the developing embryo. The process involves transcriptional regulation of genes associated with the Wnt, Notch, and Fgf signaling pathways, each gene is expressed at a specific time during the somite cycle. Comparative genomics, including analysis of expression timelines may reveal the underlying regulatory modules and their causal relations, explaining the nature and origin of the segmentation mechanism. Using a deconvolution approach, we computationally reconstruct and compare the precise timelines of expression during somitogenesis in chicken and zebrafish. The result constitutes a resource that may be used for inferring possible causal relations between genes and subsequent pathways. While the sets of regulated genes and expression profiles vary between different species, notable similarities exist between the temporal organization of the pathways involved in the somite clock in chick and mouse, with certain aspects (as the phase of expression of Notch genes) conserved also in the zebrafish. The regulated genes have sequence motifs that are conserved in mouse and chicken but not zebrafish. Promoter sequence analysis suggests involvement of several transcription factors that may bind these regulatory elements, including E2F, EGR and PLAG, as well as a possible role of G-quadruplex DNA structure in regulation of the cyclic genes. Our research lays the groundwork for further studies that will probe the evolution of the regulatory mechanism of segmentation across all vertebrates.

G-Quadruplexes Involving Both Strands of Genomic DNA Are Highly Abundant and Colocalize with Functional Sites in the Human Genome.

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.

Systematic Determination of Human Cyclin Dependent Kinase (CDK)-9 Interactome Identifies Novel Functions in RNA Splicing Mediated by the DEAD Box (DDX)-5/17 RNA Helicases.

Inducible transcriptional elongation is a rapid, stereotypic mechanism for activating immediate early immune defense genes by the epithelium in response to viral pathogens. Here, the recruitment of a multifunctional complex containing the cyclin dependent kinase 9 (CDK9) triggers the process of transcriptional elongation activating resting RNA polymerase engaged with innate immune response (IIR) genes. To identify additional functional activity of the CDK9 complex, we conducted immunoprecipitation (IP) enrichment-stable isotope labeling LC-MS/MS of the CDK9 complex in unstimulated cells and from cells activated by a synthetic dsRNA, polyinosinic/polycytidylic acid [poly (I:C)]. 245 CDK9 interacting proteins were identified with high confidence in the basal state and 20 proteins in four functional classes were validated by IP-SRM-MS. These data identified that CDK9 interacts with DDX 5/17, a family of ATP-dependent RNA helicases, important in alternative RNA splicing of NFAT5, and mH2A1 mRNA two proteins controlling redox signaling. A direct comparison of the basal versus activated state was performed using stable isotope labeling and validated by IP-SRM-MS. Recruited into the CDK9 interactome in response to poly(I:C) stimulation are HSPB1, DNA dependent kinases, and cytoskeletal myosin proteins that exchange with 60S ribosomal structural proteins. An integrated human CDK9 interactome map was developed containing all known human CDK9- interacting proteins. These data were used to develop a probabilistic global map of CDK9-dependent target genes that predicted two functional states controlling distinct cellular functions, one important in immune and stress responses. The CDK9-DDX5/17 complex was shown to be functionally important by shRNA-mediated knockdown, where differential accumulation of alternatively spliced NFAT5 and mH2A1 transcripts and alterations in downstream redox signaling were seen. The requirement of CDK9 for DDX5 recruitment to NFAT5 and mH2A1 chromatin target was further demonstrated using chromatin immunoprecipitation (ChIP). These data indicate that CDK9 is a dynamic multifunctional enzyme complex mediating not only transcriptional elongation, but also alternative RNA splicing and potentially translational control.

Analysis of the TGFß-induced program in primary airway epithelial cells shows essential role of NF-κB/RelA signaling network in type II epithelial mesenchymal transition.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.