Pseudogenes as Functionally Significant Elements of the Genome

A Genotype-Phenotype Model for Predicting Resistance Training Effects on Leg Press Performance

This study develops a comprehensive genotype-phenotype model for predicting the effects of resistance training on leg press performance. A cohort of physically inactive adults (N=193) underwent 12 weeks of resistance training, and measurements of maximum isokinetic leg press peak force, muscle mass, and thickness were taken before and after the intervention. Whole-genome genotyping was performed, and genome-wide association analysis identified 85 novel SNPs significantly associated with changes in leg press strength after training. A prediction model was constructed using stepwise linear regression, incorporating seven lead SNPs that explained 40.4% of the training effect variance. The polygenic score showed a significant positive correlation with changes in leg press strength. By integrating genomic markers and phenotypic indicators, the comprehensive prediction model explained 75.4% of the variance in the training effect. Additionally, five SNPs were found to potentially impact muscle contraction, metabolism, growth, and development through their association with REACTOME pathways. Individual responses to resistance training varied, with changes in leg press strength ranging from -55.83% to 151.20%. The study highlights the importance of genetic factors in predicting training outcomes and provides insights into the potential biological functions underlying resistance training effects. The comprehensive model offers valuable guidance for personalized fitness programs based on individual genetic profiles and phenotypic characteristics.

PTEN regulates expression of its pseudogene in glioblastoma cells in DNA methylation-dependent manner

Glioblastoma (GBM) is the most aggressive and frequent type of primary brain cancer in adult patients. One of the key molecular features associated with GBM pathogenesis is the dysfunction of PTEN oncosuppressor. In addition to PTEN gene, humans and several primates possess processed PTEN pseudogene (PTENP1) that gives rise to long non-coding RNA lncPTENP1-S. Regulation and functions of PTEN and PTENP1 are highly interconnected, however, the exact molecular mechanism of how these two genes affect each other remains unclear. Here, we analyzed the methylation level of the CpG islands (CpGIs) in the promoter regions of PTEN and PTENP1 in patient-derived GBM neurospheres. We found that increased PTEN methylation corelates with decreased PTEN mRNA level. Unexpectedly, we showed the opposite trend for PTENP1. Using targeted methylation and demethylation of PTENP1 CpGI, we demonstrated that DNA methylation increases lncPTENP1-S expression in the presence of wild type PTEN protein but decreases lncPTENP1-S expression if PTEN protein is absent. Further experiments revealed that PTEN protein binds to PTENP1 promoter region and inhibits lncPTENP1-S expression if its CpGI is demethylated. Interestingly, we did not detect any effect of lncPTENP1-S on the level of PTEN mRNA, indicating that in GBM cells PTENP1 is a downstream target of PTEN rather than its upstream regulator. Finally, we studied the functions of lncPTENP1-S and demonstrated that it plays a pro-oncogenic role in GBM cells by upregulating the expression of cancer stem cell markers and decreasing cell adhesion.

A human stomach cell type transcriptome atlas

BACKGROUND

The identification of cell type-specific genes and their modification under different conditions is central to our understanding of human health and disease. The stomach, a hollow organ in the upper gastrointestinal tract, provides an acidic environment that contributes to microbial defence and facilitates the activity of secreted digestive enzymes to process food and nutrients into chyme. In contrast to other sections of the gastrointestinal tract, detailed descriptions of cell type gene enrichment profiles in the stomach are absent from the major single-cell sequencing-based atlases.

RESULTS

Here, we use an integrative correlation analysis method to predict human stomach cell type transcriptome signatures using unfractionated stomach RNAseq data from 359 individuals. We profile parietal, chief, gastric mucous, gastric enteroendocrine, mitotic, endothelial, fibroblast, macrophage, neutrophil, T-cell, and plasma cells, identifying over 1600 cell type-enriched genes.

CONCLUSIONS

We uncover the cell type expression profile of several non-coding genes strongly associated with the progression of gastric cancer and, using a sex-based subset analysis, uncover a panel of male-only chief cell-enriched genes. This study provides a roadmap to further understand human stomach biology.

Construction and analysis of pseudogene-related ceRNA network in breast cancer

Breast cancer (BC) is one of the leading causes of cancer-related deaths in women. The present study explored the potential role of pseudogenes in BC via construction and analysis of a competing endogenous RNA (ceRNA) network through a three-step process. First, we screened differentially expressed genes in nine BC datasets. Then the gene-pseudogenes pairs (nine hub genes) were selected according to the functional enrichment and correlation analysis. Second, the candidate hub genes and interacting miRNAs were used to construct the ceRNA network. Further analysis of the ceRNA network revealed a crucial ceRNA module with two genes-pseudogene pairs and two miRNAs. The in-depth analysis identified the GBP1/hsa-miR-30d-5p/GBP1P1 axis as a potential tumorigenic axis in BC patients. In the third step, the GBP1/hsa-miR-30d-5p/GBP1P1 axis expression level was assessed in 40 tumor/normal BC patients and MCF-7 cell lines. The expression of GBP1 and GBP1P1 was significantly higher in the tumor compared to the normal tissue. However, the expression of hsa-miR-30d-5p was lower in tumor samples. Then, we introduced the GBP1P1 pseudogene into the MCF-7 cell line to evaluate its effect on GBP1 and hsa-miR-30d-5p expression. As expected, the GBP1 level increased while the hsa-miR-30d-5p level decreased in the GBP1P1-overexprsssing cell line. In addition, the oncogenic properties of MCF-7 (cell viability, clonogenicity, and migration) were improved after GBP1P1 overexpression. In conclusion, we report a ceRNA network that may provide new insight into the role of pseudogenes in BC development.

Post-transcriptional regulation of erythropoiesis

Erythropoiesis is a complex, precise, and lifelong process that is essential for maintaining normal body functions. Its strict regulation is necessary to prevent a variety of blood diseases. Normal erythropoiesis is precisely regulated by an intricate network that involves transcription levels, signal transduction, and various epigenetic modifications. In recent years, research on post-transcriptional levels in erythropoiesis has expanded significantly. The dynamic regulation of splicing transitions is responsible for changes in protein isoform expression that add new functions beneficial for erythropoiesis. RNA-binding proteins adapt the translation of transcripts to the protein requirements of the cell, yielding mRNA with dynamic translation efficiency. Noncoding RNAs, such as microRNAs and lncRNAs, are indispensable for changing the translational efficiency and/or stability of targeted mRNAs to maintain the normal expression of genes related to erythropoiesis. N6-methyladenosine-dependent regulation of mRNA translation plays an important role in maintaining the expression programs of erythroid-related genes and promoting erythroid lineage determination. This review aims to describe our current understanding of the role of post-transcriptional regulation in erythropoiesis and erythroid-associated diseases, and to shed light on the physiological and pathological implications of the post-transcriptional regulation machinery in erythropoiesis. These may help to further enrich our understanding of the regulatory network of erythropoiesis and provide new strategies for the diagnosis and treatment of erythroid-related diseases.

Progeria and Aging-Omics Based Comparative Analysis

Since ancient times aging has also been regarded as a disease, and humankind has always strived to extend the natural lifespan. Analyzing the genes involved in aging and disease allows for finding important indicators and biological markers for pathologies and possible therapeutic targets. An example of the use of omics technologies is the research regarding aging and the rare and fatal premature aging syndrome progeria (Hutchinson-Gilford progeria syndrome, HGPS). In our study, we focused on the in silico analysis of differentially expressed genes (DEGs) in progeria and aging, using a publicly available RNA-Seq dataset (GEO dataset GSE113957) and a variety of bioinformatics tools. Despite the GSE113957 RNA-Seq dataset being well-known and frequently analyzed, the RNA-Seq data shared by Fleischer et al. is far from exhausted and reusing and repurposing the data still reveals new insights. By analyzing the literature citing the use of the dataset and subsequently conducting a comparative analysis comparing the RNA-Seq data analyses of different subsets of the dataset (healthy children, nonagenarians and progeria patients), we identified several genes involved in both natural aging and progeria (KRT8, KRT18, ACKR4, CCL2, UCP2, ADAMTS15, ACTN4P1, WNT16, IGFBP2). Further analyzing these genes and the pathways involved indicated their possible roles in aging, suggesting the need for further in vitro and in vivo research. In this paper, we (1) compare "normal aging" (nonagenarians vs. healthy children) and progeria (HGPS patients vs. healthy children), (2) enlist genes possibly involved in both the natural aging process and progeria, including the first mention of IGFBP2 in progeria, (3) predict miRNAs and interactomes for WNT16 (hsa-mir-181a-5p), UCP2 (hsa-mir-26a-5p and hsa-mir-124-3p), and IGFBP2 (hsa-mir-124-3p, hsa-mir-126-3p, and hsa-mir-27b-3p), (4) demonstrate the compatibility of well-established R packages for RNA-Seq analysis for researchers interested but not yet familiar with this kind of analysis, and (5) present comparative proteomics analyses to show an association between our RNA-Seq data analyses and corresponding changes in protein expression.

Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa

BACKGROUND

The strawberry fleshy fruit is actually enlarged receptacle tissue, and the successful development of the embryo and endosperm is essential for receptacle fruit set. MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play indispensable regulatory roles in plant growth and development. However, miRNAs and phasiRNAs participating in the regulation of strawberry embryo and endosperm development have yet to be explored.

RESULTS

Here, we performed genome-wide identification of miRNA and phasiRNA-producing loci (PHAS) in strawberry seeds with a focus on those involved in the development of the early embryo and endosperm. We found that embryos and endosperm have different levels of small RNAs. After bioinformatics analysis, the results showed that a total of 404 miRNAs (352 known and 52 novel) and 156 PHAS genes (81 21-nt and 75 24-nt genes) could be found in strawberry seed-related tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. Based on refined putative annotation of PHAS loci, some auxin signal-related genes, such as CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that IAA biosynthesis is important for endosperm and embryo development during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified.

CONCLUSIONS

Taken together, these results expand our understanding of sRNAs in strawberry embryo and endosperm development and provide a genomic resource for early-stage fruit development.

The pseudogenes of eukaryotic translation elongation factors (EEFs): Role in cancer and other human diseases

The eukaryotic translation elongation factors (EEFs), i.e. EEF1A1, EEF1A2, EEF1B2, EEF1D, EEF1G, EEF1E1 and EEF2, are coding-genes that play a central role in the elongation step of translation but are often altered in cancer. Less investigated are their pseudogenes. Recently, it was demonstrated that pseudogenes have a key regulatory role in the cell, especially via non-coding RNAs, and that the aberrant expression of ncRNAs has an important role in cancer development and progression. The present review paper, for the first time, collects all that published about the EEFs pseudogenes to create a base for future investigations. For most of them, the studies are in their infancy, while for others the studies suggest their involvement in normal cell physiology but also in various human diseases. However, more investigations are needed to understand their functions in both normal and cancer cells and to define which can be useful biomarkers or therapeutic targets.

Single-cell sequencing reveals the new existence form of dermal papilla cells in the hair follicle regeneration of cashmere goats

The problem of human hair loss has caused widespread concern, however, such research is difficult because the periodicity is not obvious and the deeper levels knowledge of dermal papilla (DP) stem cells' differentiation are limited. Here, cashmere goats which have obvious periodicity of hair follicles were used, based on unbiased scRNA sequencing, we constructed DP cell lineage differentiation trajectory and revealed the key genes, signals and functions involved in cell fate decisions. And then we revealed the molecular landscape of hair follicle on regeneration. Revealed that DP cells differentiate into four intermediate cell states at different periodicity: Intermediate-cell-10 showed important functions in the growth and maintenance of cashmere; intermediate-cell-1 acting on apoptosis and cashmere shedding; intermediate-cell-0 initiated new follicular cycles, the migration of hair follicles and the occurrence of cashmere; and intermediate-cell-15 are suggested to be DP progenitor cells. In general, we provide new insights for hair regrowth. At the same time, it provides a new research ideas, directions and molecular landscape for the mechanism of dermal papilla cells.

Pseudogenes and their potential functions in hematopoiesis

Pseudogenes are DNA regions comprising defective copies of functional genes, the majority of which were generated by RNA- or DNA-level duplications. They exist across almost all forms of life and account for about one-quarter of the annotated genes in the human genome. Although these have been considered nonfunctional for decades, a growing number of pseudogenes have been found to be transcribed and to play crucial regulatory roles. Accumulating evidence indicates that they regulate gene expression through molecular interactions at the protein, RNA, and DNA levels. However, pseudogenes are often excluded in multiple genomewide analyses and functional screening, and their biological activities remain to be systematically disclosed. Here, we summarize the features of and progress of research on pseudogenes, in addition to discussing what is unknown about these genetic elements. Our previous findings, together with evidence of their poor conservation, prompted us to propose that pseudogenes may contribute to primate- or human-specific regulation, especially in hematopoiesis.

Pseudogene HSPB1P1 contributes to renal cell carcinoma proliferation and metastasis by targeting miR-296-5p to regulate HMGA1 expression

With the aid of next-generation sequencing technology, pseudogenes have been widely recognized as functional regulators in the development and progression of certain diseases, especially cancer. Our present study aimed to investigate the functions and molecular mechanisms of HSPB1-associated protein 1 pseudogene 1 (HSPB1P1) in renal cell carcinoma (RCC). HSPB1P1 expression at the mRNA levels was determined by quantitative real-time polymerase chain reaction, and its clinical significance was assessed. Cell viability was detected by Cell Counting Kit-8 assay. Cell migration and invasion were detected by transwell assays. The location of HSPB1P1 in RCC cells was detected by subcellular distribution analysis. The direct relationship between HSPB1P1 and miR-296-5p/HMGA1 axis was verified by dual-luciferase reporter assay and RNA immunoprecipitation assay. Our results identify the elevated expression of HSPB1P1 in RCC tissues and cell lines, which predicted advanced progression and poor prognosis in patients with RCC. Knockdown of HSPB1P1 suppressed cell proliferation, migration, and invasion, and reversed epithelial-mesenchymal transition process in RCC. HSPB1P1 was mostly enriched in the cytoplasm and functioned as a miRNA sponge for miR-296-5p and then regulated high-mobility group A1 expression. In conclusion, our study indicated that HSPB1P1 contributed to RCC progression by targeting the miR-296-5p/HMGA1 axis, and should be considered as a promising biomarker and therapeutic target for clinical applications.

Benchmarking and Testing Machine Learning Approaches with BARRA:CuRDa, a Curated RNA-Seq Database for Cancer Research

RNA-seq is gradually becoming the dominating technique employed to access the global gene expression in biological samples, allowing more flexible protocols and robust analysis. However, the nature of RNA-seq results imposes new data-handling challenges when it comes to computational analysis. With the increasing employment of machine learning (ML) techniques in biomedical sciences, databases that could provide curated data sets treated with state-of-the-art approaches already adapted to ML protocols, become essential for testing new algorithms. In this study, we present the Benchmarking of ARtificial intelligence Research: Curated RNA-seq Database (BARRA:CuRDa). BARRA:CuRDa was built exclusively for cancer research and is composed of 17 handpicked RNA-seq data sets for Homo sapiens that were gathered from the Gene Expression Omnibus, using rigorous filtering criteria. All data sets were individually submitted to sample quality analysis, removal of low-quality bases and artifacts from the experimental process, removal of ribosomal RNA, and estimation of transcript-level abundance. Moreover, all data sets were tested using standard approaches in the field, which allows them to be used as benchmark to new ML approaches. A feature selection analysis was also performed on each data set to investigate the biological accuracy of basic techniques. Results include genes already related to their specific tumoral tissue a large amount of long noncoding RNA and pseudogenes. BARRA:CuRDa is available at http://sbcb.inf.ufrgs.br/barracurda.

Pseudogene ANXA2P2 knockdown shows tumor-suppressive function by inhibition of the PI3K/PKB pathway in glioblastoma cells

The pseudogene annexin A2 pseudogene 2 (ANXA2P2) is highly expressed in glioblastoma (GBM). However, its role and mechanism involved in the progression of GBM remain poorly understood. ANXA2P2 messenger RNA expression was measured by quantitative reverse transcription-polymerase chain reaction. The protein levels were detected by Western blot. Cell viability was evaluated by the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide and lactate dehydrogenase (LDH) release assays. Cell invasive ability was investigated by the transwell assay and by epithelial-mesenchymal transition (EMT). Cell apoptosis was examined by flow cytometry. The results showed that ANXA2P2 expression was increased in GBM tissues and cells. Silencing of ANXA2P2 inhibited the activation of the phosphoinositide 3-kinase (PI3K)/protein kinase B (PKB) pathway in GBM cells. Knockdown of ANXA2P2 decreased cell viability, promoted LDH release, suppressed cell invasive ability, and EMT, and induced cell apoptosis in GBM cells. The addition of the PI3K/PKB activator 740Y-P abrogated the effects of ANXA2P2 knockdown on cell viability, LDH release, invasive ability, and apoptosis. In conclusion, knockdown of ANXA2P2 inhibited cell viability and invasion but promoted the apoptotic rate by suppressing the PI3K/PKB pathway in GBM cells. ANXA2P2 may represent a new target for the treatment of GBM.

Identification of Transcriptomic Differences between Lower Extremities Arterial Disease, Abdominal Aortic Aneurysm and Chronic Venous Disease in Peripheral Blood Mononuclear Cells Specimens

Several human tissues are investigated in studies of molecular biomarkers associated with diseases development. Special attention is focused on the blood and its components due to combining abundant information about systemic responses to pathological processes as well as high accessibility. In the current study, transcriptome profiles of peripheral blood mononuclear cells (PBMCs) were used to compare differentially expressed genes between patients with lower extremities arterial disease (LEAD), abdominal aortic aneurysm (AAA) and chronic venous disease (CVD). Gene expression patterns were generated using the Ion S5XL next-generation sequencing platform and were analyzed using DESeq2 and UVE-PLS methods implemented in R programming software. In direct pairwise analysis, 21, 58 and 10 differentially expressed genes were selected from the comparison of LEAD vs. AAA, LEAD vs. CVD and AAA vs. CVD patient groups, respectively. Relationships between expression of dysregulated genes and age, body mass index, creatinine levels, hypertension and medication were identified using Spearman rank correlation test and two-sided Mann-Whitney U test. The functional analysis, performed using DAVID website tool, provides potential implications of selected genes in pathological processes underlying diseases studied. Presented research provides new insight into differences of pathogenesis in LEAD, AAA and CVD, and selected genes could be considered as potential candidates for biomarkers useful in diagnosis and differentiation of studied diseases.

Genome-wide analysis of pseudogenes reveals HBBP1's human-specific essentiality in erythropoiesis and implication in β-thalassemia

The human genome harbors 14,000 duplicated or retroposed pseudogenes. Given their functionality as regulatory RNAs and low conservation, we hypothesized that pseudogenes could shape human-specific phenotypes. To test this, we performed co-expression analyses and found that pseudogene exhibited tissue-specific expression, especially in the bone marrow. By incorporating genetic data, we identified a bone-marrow-specific duplicated pseudogene, HBBP1 (η-globin), which has been implicated in β-thalassemia. Extensive functional assays demonstrated that HBBP1 is essential for erythropoiesis by binding the RNA-binding protein (RBP), HNRNPA1, to upregulate TAL1, a key regulator of erythropoiesis. The HBBP1/TAL1 interaction contributes to a milder symptom in β-thalassemia patients. Comparative studies further indicated that the HBBP1/TAL1 interaction is human-specific. Genome-wide analyses showed that duplicated pseudogenes are often bound by RBPs and less commonly bound by microRNAs compared with retropseudogenes. Taken together, we not only demonstrate that pseudogenes can drive human evolution but also provide insights on their functional landscapes.

FTH1 Pseudogenes in Cancer and Cell Metabolism

Ferritin, the principal intracellular iron-storage protein localized in the cytoplasm, nucleus, and mitochondria, plays a major role in iron metabolism. The encoding ferritin genes are members of a multigene family that includes some pseudogenes. Even though pseudogenes have been initially considered as relics of ancient genes or junk DNA devoid of function, their role in controlling gene expression in normal and transformed cells has recently been re-evaluated. Numerous studies have revealed that some pseudogenes compete with their parental gene for binding to the microRNAs (miRNAs), while others generate small interference RNAs (siRNAs) to decrease functional gene expression, and still others encode functional mutated proteins. Consequently, pseudogenes can be considered as actual master regulators of numerous biological processes. Here, we provide a detailed classification and description of the structural features of the ferritin pseudogenes known to date and review the recent evidence on their mutual interrelation within the complex regulatory network of the ferritin gene family.

HLA-J, a Non-Pseudogene as a New Prognostic Marker for Therapy Response and Survival in Breast Cancer

The human leukocyte antigen (HLA) genes are cell-surface proteins, essential for immune cell interaction. HLA-G is known for their high immunosuppressive effect and its potential as predictive marker in breast cancer. However, nothing is known about the HLA-J and its immunosuppressive, prognostic and predictive features, as it is assumed to be a "pseudogene" by in silico sequence interpretation. HLA-J, ESR1, ERBB2, KRT5 and KRT20 mRNA expression were analysed in 29 fresh frozen breast cancer biopsies and their corresponding resectates obtained from patients treated with neoadjuvant chemotherapy (NACT). mRNA was analysed with gene specific TaqMan-based Primer/Probe sets and normalized to Calmodulin 2. All breast cancer samples did express HLA-J and frequently increased HLA-J mRNA levels after NACT. HLA-J mRNA was significantly associated with overexpression of the ESR1 mRNA status (Spearman ρ 0,5679; p = 0.0090) and KRT5 mRNA (Spearman ρ 0,6121; p = 0.0041) in breast cancer core biopsies and dominated in luminal B subtype. Kaplan Meier analysis revealed that an increase of HLA-J mRNA expression after NACT had worse progression free survival (p = 0,0096), indicating a counterreaction of tumor tissues presumably to prevent elimination by enhanced immune infiltration induced by NACT. This counterreaction is associated with worse prognosis. To our knowledge this is the first study identifying HLA-J as a new predictive marker in breast cancer being involved in immune evasion mechanisms.

Functional characterization of long noncoding RNAs

PURPOSE OF REVIEW

Mounting evidence suggests that long noncoding RNAs (lncRNAs) are essential regulators of gene expression. Although few lncRNAs have been the subject of detailed molecular and functional characterization, it is believed that lncRNAs play an important role in tissue homeostasis and development. In fact, gene expression profiling studies reveal lncRNAs are developmentally regulated in a tissue-type and cell-type specific manner. Such findings have brought significant attention to their potential contribution to disease cause. The current review summarizes recent studies of lncRNAs in the heart.

RECENT FINDINGS

lncRNA discovery has largely been driven by the implementation of next generation sequencing technologies. To date, such technologies have contributed to the identification of tens of thousands of distinct lncRNAs in humans -- accounting for a large majority of all RNA sequences transcribed across the human genome. Although the functions of these lncRNAs remain largely unknown, gain-of-function and loss-of-function studies (in vivo and in vitro) have uncovered a number of mechanisms by which lncRNAs regulate gene expression and protein function. Such mechanisms have been stratified according to three major functional categories: RNA sponges (RNA-mediated sequestration of free miRNAs; e.g. H19, MEG3, and MALAT1); transcription-modulating lncRNAs (RNA influences regulatory factor recruitment by binding to histone modifiers or transcription factors; e.g. CAIF, MANTIS, and NEAT1); and translation-modulating lncRNAs (RNA modifies protein function via directly interacting with a protein itself or binding partners; e.g. Airn, CCRR, and ZFAS1).

SUMMARY

Recent studies strongly suggest that lncRNAs function via binding to macromolecules (e.g. genomic DNA, miRNAs, or proteins). Thus, lncRNAs constitute an additional mode by which cells regulate gene expression.

Protein-Coding Genes in Euarchontoglires with Pseudogene Homologs in Humans

An original bioinformatics technique is developed to identify the protein-coding genes in rodents, lagomorphs and nonhuman primates that are pseudogenized in humans. The method is based on per-gene verification of local synteny, similarity of exon-intronic structures and orthology in a set of genomes. It is applicable to any genome set, even with the number of genomes exceeding 100, and efficiently implemented using fast computer software. Only 50 evolutionary recent human pseudogenes were predicted. Their functional homologs in model species are often associated with the immune system or digestion and mainly express in the testes. According to current evidence, knockout of most of these genes leads to an abnormal phenotype. Some genes were pseudogenized or lost independently in human and nonhuman hominoids.

Brain cell somatic gene recombination and its phylogenetic foundations

A new form of somatic gene recombination (SGR) has been identified in the human brain that affects the Alzheimer's disease gene, amyloid precursor protein (APP). SGR occurs when a gene sequence is cut and recombined within a single cell's genomic DNA, generally independent of DNA replication and the cell cycle. The newly identified brain SGR produces genomic complementary DNAs (gencDNAs) lacking introns, which integrate into locations distinct from germline loci. This brief review will present an overview of likely related recombination mechanisms and genomic cDNA-like sequences that implicate evolutionary origins for brain SGR. Similarities and differences exist between brain SGR and VDJ recombination in the immune system, the first identified SGR form that now has a well-defined enzymatic machinery. Both require gene transcription, but brain SGR uses an RNA intermediate and reverse transcriptase (RT) activity, which are characteristics shared with endogenous retrotransposons. The identified gencDNAs have similarities to other cDNA-like sequences existing throughout phylogeny, including intron-less genes and inactive germline processed pseudogenes, with likely overlapping biosynthetic processes. gencDNAs arise somatically in an individual to produce multiple copies; can be functional; appear most frequently within postmitotic cells; have diverse sequences; change with age; and can change with disease state. Normally occurring brain SGR may represent a mechanism for gene optimization and long-term cellular memory, whereas its dysregulation could underlie multiple brain disorders and, potentially, other diseases like cancer. The involvement of RT activity implicates already Food and Drug Administration-approved RT inhibitors as possible near-term interventions for managing SGR-associated diseases and suggest next-generation therapeutics targeting SGR elements.

A mega-analysis of expression quantitative trait loci in retinal tissue

Significant association signals from genome-wide association studies (GWAS) point to genomic regions of interest. However, for most loci the causative genetic variant remains undefined. Determining expression quantitative trait loci (eQTL) in a disease relevant tissue is an excellent approach to zoom in on disease- or trait-associated association signals and hitherto on relevant disease mechanisms. To this end, we explored regulation of gene expression in healthy retina (n = 311) and generated the largest cis-eQTL data set available to date. Genotype- and RNA-Seq data underwent rigorous quality control protocols before FastQTL was applied to assess the influence of genetic markers on local (cis) gene expression. Our analysis identified 403,151 significant eQTL variants (eVariants) that regulate 3,007 genes (eGenes) (Q-Value < 0.05). A conditional analysis revealed 744 independent secondary eQTL signals for 598 of the 3,007 eGenes. Interestingly, 99,165 (24.71%) of all unique eVariants regulate the expression of more than one eGene. Filtering the dataset for eVariants regulating three or more eGenes revealed 96 potential regulatory clusters. Of these, 31 harbour 130 genes which are partially regulated by the same genetic signal. To correlate eQTL and association signals, GWAS data from twelve complex eye diseases or traits were included and resulted in identification of 80 eGenes with potential association. Remarkably, expression of 10 genes is regulated by eVariants associated with multiple eye diseases or traits. In conclusion, we generated a unique catalogue of gene expression regulation in healthy retinal tissue and applied this resource to identify potentially pleiotropic effects in highly prevalent human eye diseases. Our study provides an excellent basis to further explore mechanisms of various retinal disease etiologies.

The retina is a multilayered and highly specified neural tissue crucial for high-resolution visual perception and spatial orientation. Environmental and genetic insults to the retina result in many blinding diseases, such as age-related macular degeneration or glaucoma. Commonly, many of these diseases are age-related suggesting that minor changes are accumulating over a life-time, with little or no contribution of strong individual effects. Specifically, this is true for genetic factors known to underlie the etiology of complex diseases including the prevalent eye diseases. In our study, we searched for effects on gene expression due to genetic variation using 311 healthy post-mortem retinal tissue samples. We show that 3,007 of the 16,766 genes investigated are regulated in the retina by genetic variations. Of these, 80 genes are potentially associated to one or more of twelve complex eye diseases or retinal traits tested. Interestingly, 10 genes appear to be involved in the development of several eye traits suggesting that cellular mechanisms may act at a common point in the disease process. Consequently, our study provides the basis to further explore retinal disease pathways and is likely to highlight target molecules for future therapeutic applications.

The Oncogenic Roles of PTTG1 and PTTG2 Genes and Pseudogene PTTG3P in Head and Neck Squamous Cell Carcinomas

Background: Head and neck squamous cell carcinomas are a group of heterogeneous diseases that occur in the mouth, pharynx and larynx and are characterized by poor prognosis. A low overall survival rate leads to a need to develop biomarkers for early head and neck squamous cell carcinomas detection, accurate prognosis and appropriate selection of therapy. Therefore, in this paper, we investigate the biological role of the PTTG3P pseudogene and associated genes PTTG1 and PTTG2 and their potential use as biomarkers. Methods: Based on TCGA data and the UALCAN database, PTTG3P, PTTG1 and PTTG2 expression profiles and clinicopathological features with TP53 gene status as well as expression levels of correlated genes were analyzed in patients’ tissue samples. The selected genes were classified according to their biological function using the PANTHER tool. Gene Set Enrichment Analysis software was used for functional enrichment analysis. All statistical analyses were performed using GraphPad Prism 5. Results: In head and neck squamous cell carcinomas, significant up-regulation of the PTTG3P pseudogene, PTTG1 and PTTG2 genes’ expression between normal and cancer samples were observed. Moreover, the expression of PTTG3P, PTTG1 and PTTG2 depends on the type of mutation in TP53 gene, and they correlate with genes from p53 pathway. PTTG3P expression was significantly correlated with PTTG1 as well as PTTG2, as was PTTG1 expression with PTTG2. Significant differences between expression levels of PTTG3P, PTTG1 and PTTG2 in head and neck squamous cell carcinomas patients were also observed in clinicopathological contexts. The contexts taken into consideration included: T-stage for PTTG3P; grade for PTTG3, PTTG1 and PTTG2; perineural invasion and lymph node neck dissection for PTTG1 and HPV p16 status for PTTG3P, PTTG1 and PTTG2. A significantly longer disease-free survival for patients with low expressions of PTTG3P and PTTG2, as compared to high expression groups, was also observed. Gene Set Enrichment Analysis indicated that the PTTG3 high-expressing group of patients have the most deregulated genes connected with DNA repair, oxidative phosphorylation and peroxisome pathways. For PTTG1, altered genes are from DNA repair groups, Myc targets, E2F targets and oxidative phosphorylation pathways, while for PTTG2, changes in E2F targets, G2M checkpoints and oxidative phosphorylation pathways are indicated. Conclusions: PTTG3P and PTTG2 can be used as a prognostic biomarker in head and neck squamous cell carcinomas diagnostics. Moreover, patients with high expressions of PTTG3P, PTTG1 or PTTG2 have worse outcomes due to upregulation of oncogenic pathways and more aggressive phenotypes.

YRNAs and YRNA-Derived Fragments as New Players in Cancer Research and Their Potential Role in Diagnostics

YRNAs are a type of short, noncoding RNAs. A total of four different transcripts can be distinguished, which are YRNA1, YRNA3, YRNA4 and YRNA5. All YRNAs are relatively small, made up of about 100 nucleotides each. YRNAs are characterized by a stem-loop structure and each part of that structure carries a different function. YRNAs are transcribed in the nucleus by RNA polymerase III. Then, the YRNA molecule is bound to the polyuridine tail of the La protein responsible for both its nuclear retention and protection from degradation. They also bind to the Ro60 protein, making the molecule more stable. In turn, YRNA-derived small RNAs (YsRNAs) are a class of YRNAs produced in apoptotic cells as a result of YRNA degradation. This process is performed by caspase-3-dependent pathways that form two groups of YsRNAs, with lengths of either approximately 24 or 31 nucleotides. From all four YRNA transcripts, 75 well-described pseudogenes are generated as a result of the mutation. However, available data indicates the formation of up to 1000 pseudogenes. YRNAs and YRNA-derived small RNAs may play a role in carcinogenesis due to their altered expression in cancers and influence on cell proliferation and inflammation. Nevertheless, our knowledge is still limited, and more research is required. The main aim of this review is to describe the current state of knowledge about YRNAs, their function and contribution to carcinogenesis, as well as their potential role in cancer diagnostics. To confirm the promising potential of YRNAs and YRNA-derived fragments as biomarkers, their significant role in several tumor types was taken into consideration.

Comparative genomic analysis of the 3' UTR of human MDM2 identifies multiple transposable elements, an RLP24 pseudogene and a cluster of novel repeat sequences that arose during primate evolution

The MDM2 oncogene is a negative regulator of the p53 tumour suppressor. This relationship appears to have originated over a billion years ago. The human MDM2 gene encodes a variety of mRNAs with exceptionally long 3'UTRs (up to 5.7 kb); however, it was unclear whether MDM2 3'UTRs from other species are similarly long or conserved at the sequence level. Here, we report that all but one of the primate species most closely related to humans (greater and lesser apes) have similarly long 3'UTRs with high sequence similarity across their entire length. More distantly related species (Old world monkeys and new world monkeys) tend to have shorter MDM2 3'UTRs homologous to the corresponding position of the human MDM2 3'UTR while non-primate species exhibit little similarity at all. Remarkably, DNA sequences downstream of the shorter primate 3'UTRs are syntenic with distal regions in the human and other ape MDM2 3'UTRs. These homologous non-transcribed intergenic and transcribed 3'UTR-encoding regions are comprised of a variety of transposable elements, an RLP24 pseudogene and a cluster of novel repeat sequences suggestive of another unknown transposable element. Our analysis suggests that the primary difference between long and short MDM2 3'UTRs is a switch in polyA site usage to include conserved transposable elements that remain intergenic in more distantly related primates. It will be important to determine the relative contribution of these elements to post-transcriptional and translational regulation of MDM2 and hence p53-mediated tumour suppression.

Genomic Indexing by Somatic Gene Recombination of mRNA/ncRNA - Does It Play a Role in Genomic Mosaicism, Memory Formation, and Alzheimer's Disease?

Recent evidence indicates that genomic individuality of neurons, characterized by DNA-content variation, is a common if not universal phenomenon in the human brain that occurs naturally but can also show aberrancies that have been linked to the pathomechanism of Alzheimer’s disease and related neurodegenerative disorders. Etiologically, this genomic mosaic has been suggested to arise from defects of cell cycle regulation that may occur either during brain development or in the mature brain after terminal differentiation of neurons. Here, we aim to draw attention towards another mechanism that can give rise to genomic individuality of neurons, with far-reaching consequences. This mechanism has its origin in the transcriptome rather than in replication defects of the genome, i.e., somatic gene recombination of RNA. We continue to develop the concept that somatic gene recombination of RNA provides a physiological process that, through integration of intronless mRNA/ncRNA into the genome, allows a particular functional state at the level of the individual neuron to be indexed. By insertion of defined RNAs in a somatic recombination process, the presence of specific mRNA transcripts within a definite temporal context can be “frozen” and can serve as an index that can be recalled at any later point in time. This allows information related to a specific neuronal state of differentiation and/or activity relevant to a memory trace to be fixed. We suggest that this process is used throughout the lifetime of each neuron and might have both advantageous and deleterious consequences.

Genomic landscape of the signals of positive natural selection in populations of Northern Eurasia: A view from Northern Russia

Natural selection of beneficial genetic variants played a critical role in human adaptation to a wide range of environmental conditions. Northern Eurasia, despite its severe climate, is home to lots of ethnically diverse populations. The genetic variants associated with the survival of these populations have hardly been analyzed. We searched for the genomic signatures of positive selection in (1) the genome-wide microarray data of 432 people from eight different northern Russian populations and (2) the whole-genome sequences of 250 people from Northern Eurasia from a public repository through testing the extended haplotype homozigosity (EHH) and direct comparison of allele frequency, respectively. The 20 loci with the strongest selection signals were characterized in detail. Among the top EHH hits were the NRG3 and NBEA genes, which are involved in the development and functioning of the neural system, the PTPRM gene, which mediates cell-cell interactions and adhesion, and a region on chromosome 4 (chr4:28.7-28.9 Mb) that contained several loci affiliated with different classes of non-coding RNAs (RN7SL101P, MIR4275, MESTP3, and LINC02364). NBEA and the region on chromosome 4 were novel selection targets that were identified for the first time in Western Siberian populations. Cross-population comparisons of EHH profiles suggested a particular role for the chr4:28.7-28.9 Mb region in the local adaptation of Western Siberians. The strongest selection signal identified in Siberian sequenced genomes was formed by six SNPs on chromosome 11 (chr11:124.9-125.2 Mb). This region included well-known genes SLC37A2 and PKNOX2. SLC37A2 is most-highly expressed in the gut. Its expression is regulated by vitamin D, which is often deficient in northern regions. The PKNOX2 gene is a transcription factor of the homeobox family that is expressed in the brain and many other tissues. This gene is associated with alcohol addiction, which is widespread in many Northern Eurasian populations.

Unveiling ncRNA regulatory axes in atherosclerosis progression

Completion of the human genome sequencing project highlighted the richness of the cellular RNA world, and opened the door to the discovery of a plethora of short and long non-coding RNAs (the dark transcriptome) with regulatory or structural potential, which shifted the balance of pathological gene alterations from coding to non-coding RNAs. Thus, disease risk assessment currently has to also evaluate the expression of new RNAs such as small micro RNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), competing endogenous RNAs (ceRNAs), retrogressed elements, 3'UTRs of mRNAs, etc. We are interested in the pathogenic mechanisms of atherosclerosis (ATH) progression in patients suffering Chronic Kidney Disease, and in this review, we will focus in the role of the dark transcriptome (non-coding RNAs) in ATH progression. We will focus in miRNAs and in the formation of regulatory axes or networks with their mRNA targets and with the lncRNAs that function as miRNA sponges or competitive inhibitors of miRNA activity. In this sense, we will pay special attention to retrogressed genomic elements, such as processed pseudogenes and Alu repeated elements, that have been recently seen to also function as miRNA sponges, as well as to the use or miRNA derivatives in gene silencing, anti-ATH therapies. Along the review, we will discuss technical developments associated to research in lncRNAs, from sequencing technologies to databases, repositories and algorithms to predict miRNA targets, as well as new approaches to miRNA function, such as integrative or enrichment analysis and their potential to unveil RNA regulatory networks.

Tspan8-Tumor Extracellular Vesicle-Induced Endothelial Cell and Fibroblast Remodeling Relies on the Target Cell-Selective Response

Tumor cell-derived extracellular vesicles (TEX) expressing tetraspanin Tspan8-alpha4/beta1 support angiogenesis. Tspan8-alpha6/beta4 facilitates lung premetastatic niche establishment. TEX-promoted target reprogramming is still being disputed, we explored rat endothelial cell (EC) and lung fibroblast (Fb) mRNA and miRNA profile changes after coculture with TEX. TEX were derived from non-metastatic BSp73AS (AS) or metastatic BSp73ASML (ASML) rat tumor lines transfected with Tspan8 (AS-Tspan8) or Tspan8-shRNA (ASML-Tspan8kd). mRNA was analyzed by deep sequencing and miRNA by array analysis of EC and Fb before and after coculture with TEX. EC and Fb responded more vigorously to AS-Tspan8- than AS-TEX. Though EC and Fb responses differed, both cell lines predominantly responded to membrane receptor activation with upregulation and activation of signaling molecules and transcription factors. Minor TEX-initiated changes in the miRNA profile relied, at least partly, on long noncoding RNA (lncRNA) that also affected chromosome organization and mRNA processing. These analyses uncovered three important points. TEX activate target cell autonomous programs. Responses are initiated by TEX targeting units and are target cell-specific. The strong TEX-promoted lncRNA impact reflects lncRNA shuttling and location-dependent distinct activities. These informations urge for an in depth exploration on the mode of TEX-initiated target cell-specific remodeling including, as a major factor, lncRNA.

5-Hydroxymethylcytosine signatures in circulating cell-free DNA as diagnostic and predictive biomarkers for coronary artery disease

Background

The 5-hydroxymethylcytosine (5hmC) DNA modification is an epigenetic marker involved in a range of biological processes. Its function has been studied extensively in tumors, neurodegenerative diseases, and atherosclerosis. Studies have reported that 5hmC modification is closely related to the phenotype transformation of vascular smooth muscle cells and endothelial dysfunction. However, its role in coronary artery disease (CAD) has not been fully studied.

Results

To investigate whether 5hmC modification correlates with CAD pathogenesis and whether 5hmC can be used as a biomarker, we used a low-input whole-genome sequencing technology based on selective chemical capture (hmC-Seal) to firstly generate the 5hmC profiles in the circulating cell-free DNA(cfDNA) of CAD patients, including stable coronary artery disease (sCAD) patients and acute myocardial infarction (AMI) patients. We detected a significant difference of 5hmC enrichment in gene bodies from CAD patients compared with normal coronary artery (NCA) individuals. Our results showed that CAD patients can be well separated from NCA individuals by 5hmC markers. The prediction performance of the model established by differentially regulated 5hmc modified genes were superior to common clinical indicators for the diagnosis of CAD (AUC = 0.93) and sCAD (AUC = 0.93). Specially, we found that 5hmC markers in cfDNA showed prediction potential for AMI (AUC = 0.95), which was superior to that of cardiac troponin I, muscle/brain creatine kinase, and myoglobin.

Conclusions

Our results suggest that 5hmC markers derived from cfDNA can serve as effective epigenetic biomarkers for minimally noninvasive diagnosis and prediction of CAD.

Quality Control Strategy for CRISPR-Cas9-Based Gene Editing Complicated by a Pseudogene

CRISPR-Cas9 mediated gene editing in induced pluripotent stem cells became an efficient tool to investigate biological mechanisms underlying genetic-driven diseases while accounting for the respective genetic background. This technique relies on the targeting of a specific nucleotide sequence present in the gene of interest. Therefore, the gene editing of some genes can be complicated by non-coding pseudogenes presenting a high homology of sequence with their respective genes. Among them, GBA is raising special interest because of its implication as the most common genetic risk factor for Parkinson’s disease. In this study, we present an easy-to-use CRISPR-Cas9 gene editing strategy allowing for specific editing of point mutations in a gene without genetic alteration of its pseudogene exemplified by the correction or insertion of the common N370S mutation in GBA. A quality control strategy by combined fluorescence and PCR-based screening allows the early identification of correctly edited clones with unambiguous identification of the status of its pseudogene, GBAP1. Successful gene editing was confirmed by functional validation. Our work presents the first CRISPR-Cas9 based editing of a point mutation in GBA and paves the way for technically demanding gene engineering due to the presence of pseudogenes.

Screening for mouse genes lost in mammals with long lifespans

Background

Gerontogenes include those that modulate life expectancy in various species and may be the actual longevity genes. We believe that a long (relative to body weight) lifespan in individual rodent and primate species can be due, among other things, to the loss of particular genes that are present in short-lived species of the same orders. These genes can also explain the widely different rates of aging among diverse species as well as why similarly sized rodents or primates sometimes have anomalous life expectancies (e.g., naked mole-rats and humans). Here, we consider the gene loss in the context of the prediction of Williams’ theory that concerns the reallocation of physiological resources of an organism between active reproduction (r-strategy) and self-maintenance (K-strategy). We have identified such lost genes using an original computer-aided approach; the software considers the loss of a gene as disruptions in gene orthology, local gene synteny or both.

Results

A method and software identifying the genes that are absent from a predefined set of species but present in another predefined set of species are suggested. Examples of such pairs of sets include long-lived vs short-lived, homeothermic vs poikilothermic, amniotic vs anamniotic, aquatic vs terrestrial, and neotenic vs nonneotenic species, among others. Species are included in one of two sets according to the property of interest, such as longevity or homeothermy. The program is universal towards these pairs, i.e., towards the underlying property, although the sets should include species with quality genome assemblies. Here, the proposed method was applied to study the longevity of Euarchontoglires species. It largely predicted genes that are highly expressed in the testis, epididymis, uterus, mammary glands, and the vomeronasal and other reproduction-related organs. This agrees with Williams’ theory that hypothesizes a species transition from r-strategy to K-strategy. For instance, the method predicts the mouse gene Smpd5, which has an expression level 20 times greater in the testis than in organs unrelated to reproduction as experimentally demonstrated elsewhere. At the same time, its paralog Smpd3 is not predicted by the program and is widely expressed in many organs not specifically related to reproduction.

Conclusions

The method and program, which were applied here to screen for gene losses that can accompany increased lifespan, were also applied to study reduced regenerative capacity and development of the telencephalon, neoteny, etc. Some of these results have been carefully tested experimentally. Therefore, we assume that the method is widely applicable.

The Pancreatic Cancer-Initiating Cell Marker CD44v6 Affects Transcription, Translation, and Signaling: Consequences for Exosome Composition and Delivery

Pancreatic cancer-initiating cells (PaCIC) express CD44v6 and Tspan8. A knockdown (kd) of these markers hinders the metastatic capacity, which can be rescued, if the cells are exposed to CIC-exosomes (TEX). Additional evidence that CD44v6 regulates Tspan8 expression prompted us to explore the impact of these PaCIC markers on nonmetastatic PaCa and PaCIC-TEX. We performed proteome, miRNA, and mRNA deep sequencing analyses on wild-type, CD44v6kd, and Tspan8kd human PaCIC and TEX. Database comparative analyses were controlled by qRT-PCR, Western blot, flow cytometry, and confocal microscopy. Transcriptome analysis of CD44 versus CD44v6 coimmunoprecipitating proteins in cells and TEX revealed that Tspan8, several signal-transducing molecules including RTK, EMT-related transcription factors, and proteins engaged in mRNA processing selectively associate with CD44v6 and that the membrane-attached CD44 intracytoplasmic tail supports Tspan8 and NOTCH transcription. Deep sequencing uncovered a CD44v6 contribution to miRNA processing. Due to the association of CD44v6 with Tspan8 in internalization prone tetraspanin-enriched membrane domains (TEM) and the engagement of Tspan8 in exosome biogenesis, most CD44v6-dependent changes were transferred into TEX such that the input of CD44v6 to TEX activities becomes largely waved in both a CD44v6kd and a Tspan8kd. Few differences between CD44v6kd- and Tspan8kd-TEX rely on CD44v6 being also recovered in non-TEM derived TEX, highlighting distinct TEX delivery from individual cells that jointly account for TEX-promoted target modulation. This leads us to propose a model in which CD44v6 strongly supports tumor progression by cooperating with signaling molecules, altering transcription of key molecules, and through its association with the mRNA processing machinery. The association of CD44v6 with Tspan8, which plays a crucial role in vesicle biogenesis, promotes metastases by transferring CD44v6 activities into TEM and TEM-independently derived TEX. Further investigations of the lead position of CD44v6 in shifting metastasis-promoting activities into CIC-TEX may offer a means of targeting TEX-CD44v6 in therapeutic applications.

Genetic modifiers of rodent animal models: the role in cataractogenesis

Visual impairment leads to a decrease in quality of life. Cataract is the most commonly observed ocular disease in humans that causes vision disorders. The risk factors associated with cataract development include aging, infections, eye injuries, environmental causes, such as radiation and exposure to ultraviolet rays in sunlight, and genetic mutations. Additionally, several cataract patients display phenotypic heterogeneity, suggesting the role of genetic modifiers in the modulation of severity and onset time of cataractogenesis. However, the genetic modifiers associated with cataract have not been identified in humans yet. In contrast, the identification and mapping of genetic modifiers have been successfully carried out in mice and rats. In this review, we focus on the genetic modifiers of cataract in the rodent models.

Tubulin mRNA stability is sensitive to change in microtubule dynamics caused by multiple physiological and toxic cues

The localization, mass, and dynamics of microtubules are important in many processes. Cells may actively monitor the state of their microtubules and respond to perturbation, but how this occurs outside mitosis is poorly understood. We used gene-expression analysis in quiescent cells to analyze responses to subtle and strong perturbation of microtubules. Genes encoding α-, β, and γ-tubulins (TUBAs, TUBBs, and TUBGs), but not δ- or ε-tubulins (TUBDs or TUBEs), exhibited the strongest differential expression response to microtubule-stabilizing versus destabilizing drugs. Quantitative PCR of exon versus intron sequences confirmed that these changes were caused by regulation of tubulin mRNA stability and not transcription. Using tubulin mRNA stability as a signature to query the Gene Expression Omnibus (GEO) database, we find that tubulin genes respond to toxins known to damage microtubules. Importantly, we find many other experimental perturbations, including multiple signaling and metabolic inputs that trigger tubulin differential expression, suggesting their novel, to our knowledge, role in the regulation of the microtubule cytoskeleton. Mechanistic follow-up confirms that one important physiological signal, phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K) activity, indeed regulates tubulin mRNA stability via changes in microtubule dynamics. We propose that tubulin gene expression is regulated as part of many coordinated biological responses, with wide implications in physiology and toxicology. Furthermore, we present a new way to discover microtubule regulation using transcriptomics.