Cell types differ in global coordination of splicing and proportion of highly expressed genes

Augmenting fibronectin levels in injured adult CNS promotes axon regeneration in vivo

In an attempt to repair injured central nervous system (CNS) nerves/tracts, immune cells are recruited into the injury site, but endogenous response in adult mammals is insufficient for promoting regeneration of severed axons. Here, we found that a portion of retinal ganglion cell (RGC) CNS projection neurons that survive after optic nerve crush (ONC) injury are enriched for and upregulate fibronectin (Fn)-interacting integrins Itga5 and ItgaV, and that Fn promotes long-term survival and long-distance axon regeneration of a portion of axotomized adult RGCs in culture. We then show that, Fn is developmentally downregulated in the axonal tracts of optic nerve and spinal cord, but injury-activated macrophages/microglia upregulate Fn while axon regeneration-promoting zymosan augments their recruitment (and thereby increases Fn levels) in the injured optic nerve. Finally, we found that Fn's RGD motif, established to interact with Itga5 and ItgaV, promotes long-term survival and long-distance axon regeneration of adult RGCs after ONC in vivo, with some axons reaching the optic chiasm when co-treated with Rpl7a gene therapy. Thus, experimentally augmenting Fn levels in the injured CNS is a promising approach for therapeutic neuroprotection and axon regeneration of at least a portion of neurons.

Nfe2l3 promotes neuroprotection and long-distance axon regeneration after injury in vivo

Nuclear factor erythroid 2 like (Nfe2l) gene family members 1-3 mediate cellular response to oxidative stress, including in the central nervous system (CNS). However, neuronal functions of Nfe2l3 are unknown. Here, we comparatively evaluated expression of Nfe2l1, Nfe2l2, and Nfe2l3 in singe cell RNA-seq (scRNA-seq)-profiled cortical and retinal ganglion cell (RGC) CNS projection neurons, investigated whether Nfe2l3 regulates neuroprotection and axon regeneration after CNS injury in vivo, and characterized a gene network associated with Nfe2l3 in neurons. We showed that, Nfe2l3 expression transiently peaks in developing immature cortical and RGC projection neurons, but is nearly abolished in adult neurons and is not upregulated after injury. Furthermore, within the retina, Nfe2l3 is enriched in RGCs, primarily neonatally, and not upregulated in injured RGCs, whereas Nfe2l1 and Nfe2l2 are expressed robustly in other retinal cell types as well and are upregulated in injured RGCs. We also found that, expressing Nfe2l3 in injured RGCs through localized intralocular viral vector delivery promotes neuroprotection and long-distance axon regeneration after optic nerve injury in vivo. Moreover, Nfe2l3 provided a similar extent of neuroprotection and axon regeneration as viral vector-targeting of Pten and Klf9, which are prominent regulators of neuroprotection and long-distance axon regeneration. Finally, we bioinformatically characterized a gene network associated with Nfe2l3 in neurons, which predicted the association of Nfe2l3 with established mechanisms of neuroprotection and axon regeneration. Thus, Nfe2l3 is a novel neuroprotection and axon regeneration-promoting factor with a therapeutic potential for treating CNS injury and disease.

The identification of small molecule inhibitors with anthelmintic activities that target conserved proteins among ruminant gastrointestinal nematodes

Gastrointestinal nematode (GIN) infections are a major concern for the ruminant industry worldwide and result in significant production losses. Naturally occurring polyparasitism and increasing drug resistance that potentiate disease outcomes are observed among the most prevalent GINs of veterinary importance. Within the five major taxonomic clades, clade Va represents a group of GINs that predominantly affect the abomasum or small intestine of ruminants. However, the development of effective broad-spectrum anthelmintics against ruminant clade Va GINs has been challenged by a lack of comprehensive druggable genome resources. Here, we first assembled draft genomes for three clade Va species (Cooperia oncophora, Trichostrongylus colubriformis, and Ostertagia ostertagi) and compared them with closely related ruminant GINs. Genome-wide phylogenetic reconstruction showed a relationship among ruminant GINs structured by taxonomic classification. Orthogroup (OG) inference and functional enrichment analyses identified 220 clade Va-specific and Va-conserved OGs, enriched for functions related to cell cycle and cellular senescence. Further transcriptomic analysis identified 61 taxonomically and functionally conserved clade Va OGs that may function as drug targets for new broad-spectrum anthelmintics. Chemogenomic screening identified 11 compounds targeting homologs of these OGs, thus having potential anthelmintic activity. In in vitro phenotypic assays, three kinase inhibitors (digitoxigenin, K-252a, and staurosporine) exhibited broad-spectrum anthelmintic activities against clade Va GINs by obstructing the motility of exsheathed L3 (xL3) or molting of xL3 to L4. These results demonstrate valuable applications of the new ruminant GIN genomes in gaining better insights into their life cycles and offer a contemporary approach to discovering the next generation of anthelmintics.IMPORTANCEGastrointestinal nematode (GIN) infections in ruminants are caused by parasites that inhibit normal function in the digestive tract of cattle, sheep, and goats, thereby causing morbidity and mortality. Coinfection and increasing drug resistance to current therapeutic agents will continue to worsen disease outcomes and impose significant production losses on domestic livestock producers worldwide. In combination with ongoing therapeutic efforts, advancing the discovery of new drugs with novel modes of action is critical for better controlling GIN infections. The significance of this study is in assembling and characterizing new GIN genomes of Cooperia oncophora, Ostertagia ostertagi, and Trichostrongylus colubriformis for facilitating a multi-omics approach to identify novel, biologically conserved drug targets for five major GINs of veterinary importance. With this information, we were then able to demonstrate the potential of commercially available compounds as new anthelmintics.

Upregulation of developmentally-downregulated miR-1247-5p promotes neuroprotection and axon regeneration in vivo

Numerous micro-RNAs (miRNAs) affect neurodevelopment and neuroprotection, but potential roles of many miRNAs in regulating these processes are still unknown. Here, we used the retinal ganglion cell (RGC) central nervous system (CNS) projection neuron and optic nerve crush (ONC) injury model, to optimize a mature miRNA arm-specific quantification method for characterizing the developmental regulation of miR-1247-5p in RGCs, investigated whether injury affects its expression, and tested whether upregulating miR-1247-5p-mimic in RGCs promotes neuroprotection and axon regeneration. We found that, miR-1247-5p is developmentally-downregulated in RGCs, and is further downregulated after ONC. Importantly, RGC-specific upregulation of miR-1247-5p promoted neuroprotection and axon regeneration after injury in vivo. To gain insight into the underlying mechanisms, we analyzed by bulk-mRNA-seq embryonic and adult RGCs, along with adult RGCs transduced by miR-1247-5p-expressing viral vector, and identified developmentally-regulated cilial and mitochondrial biological processes, which were reinstated to their embryonic levels in adult RGCs by upregulation of miR-1247-5p. Since axon growth is also a developmentally-regulated process, in which mitochondrial dynamics play important roles, it is possible that miR-1247-5p promoted neuroprotection and axon regeneration through regulating mitochondrial functions.

Experimental upregulation of developmentally downregulated ribosomal protein large subunits 7 and 7A promotes axon regeneration after injury in vivo

Ribosomal proteins are involved in neurodevelopment and central nervous system (CNS) disease and injury. However, the roles of specific ribosomal protein subunits in developmental axon growth, and their potential as therapeutic targets for treating CNS injuries, are still poorly understood. Here, we show that ribosomal protein large (Rpl) and small (Rps) subunit genes are substantially (56-fold) enriched amongst the genes, which are downregulated during maturation of retinal ganglion cell (RGC) CNS projection neurons. We also show that Rpl and Rps subunits are highly co-regulated in RGCs, and partially re-upregulated after optic nerve crush (ONC). Because developmental downregulation of ribosomal proteins coincides with developmental decline in neuronal intrinsic axon growth capacity, we hypothesized that Rpl/Rps incomplete re-upregulation after injury may be a part of the cellular response which attempts to reactivate intrinsic axon growth mechanisms. We found that experimentally upregulating Rpl7 and Rpl7A promoted axon regeneration after ONC in vivo. Finally, we characterized gene networks associated with Rpl/Rps, and showed that Rpl7 and Rpl7A belong to the cluster of genes, which are shared between translational and neurodevelopmental biological processes (based on gene-ontology) that are co-downregulated in maturing RGCs during the decline in intrinsic axon growth capacity.

Experimental gene expression of developmentally downregulated Crmp1, Crmp4, and Crmp5 promotes axon regeneration and retinal ganglion cell survival after optic nerve injury

Collapsin response mediator proteins (Crmps) play roles in neuronal development and axon growth. However, neuronal-specific roles of Crmp1, Crmp4, and Crmp5 in regeneration of injured central nervous system (CNS) axons in vivo are unclear. Here, we analyzed developmental and subtype-specific expression of Crmp genes in retinal ganglion cells (RGCs), tested whether overexpressing Crmp1, Crmp4, or Crmp5 in RGCs through localized intralocular AAV2 delivery promotes axon regeneration after optic nerve injury in vivo, and characterized developmental co-regulation of gene-concept networks associated with Crmps. We found that all Crmp genes are developmentally downregulated in RGCs during maturation. However, while Crmp1, Crmp2, and Crmp4 were expressed to a varying degree in most RGC subtypes, Crmp3 and Crmp5 were expressed only in a small subset of RGC subtypes. We then found that after optic nerve injury, Crmp1, Crmp4, and Crmp5 promote RGC axon regeneration to varying extents, with Crmp4 promoting the most axon regeneration and also localizing to axons. We also found that Crmp1 and Crmp4, but not Crmp5, promote RGC survival. Finally, we found that Crmp1, Crmp2, Crmp4, and Crmp5's ability to promote axon regeneration is associated with neurodevelopmental mechanisms, which control RGC's intrinsic axon growth capacity.

Excessive Sodium Intake Leads to Cardiovascular Disease by Promoting Sex-Specific Dysfunction of Murine Heart

Background

Globally, a high-salt diet (HSD) has become a threat to human health as it can lead to a high risk of cardiac damage. Although some studies investigating HSD have been carried out, the majority has been conducted in males, and there are few female-specific studies, thereby ignoring any effects of sex-specific damage on the heart. In this study, we determined how HSD induces different pathways of cardiovascular diseases through sex-specific effects on cardiac damage in mice.

Methods

An HSD murine model of male and female C57BL/6J mice was fed with sodium-rich chow (4% NaCl). After 8 weeks, cardiac tissues were collected, and the whole gene transcriptome of the hearts of male and female mice was characterized and analyzed using high-throughput RNA sequencing. Immunohistochemistry staining was used to further assess the harmful effects of HSD on protein expression of genes associated with immunity, fibrosis, and apoptosis in male and female mice.

Results

HSD drastically altered the cardiac transcriptome compared to that of the normal heart in both male and female mice and had a sex-specific effect on the cardiac composition in the transcriptome. HSD produced various differentially expressed genes and affected different KEGG pathways of the transcriptome in male and female mice. Furthermore, we found that HSD induced different pathways of cardiovascular disease in the male mice and female mice. The pathway of hypertrophic cardiomyopathy is significantly enriched in HSD-treated male mice, while the pathway of dilated cardiomyopathy is significantly enriched in HSD-treated female mice. Finally, metabolism, immunity, fibrosis, and apoptosis in the mouse heart showed sex-specific changes predicting cardiac damage.

Conclusion

Our results demonstrate that HSD adversely impacts cardiac structure and function by affecting the metabolism, immunity, fibrosis, and apoptosis in the murine heart and induces the mouse to suffer from sex-specific cardiovascular disease. This study provides a new perspective and basis for the differences in the pharmacology and interventional treatment of sex-specific cardiovascular diseases induced by HSD in men and women.

Lactation Associated Genes Revealed in Holstein Dairy Cows by Weighted Gene Co-Expression Network Analysis (WGCNA)

Simple Summary

Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In the past few years, studies on the gene expression changes of dairy cow mammary glands were only based on transcriptome comparisons between two lactation stages. Few studies focused on the relationships between gene expression of the dairy mammary gland and lactation stage or milk composition in a lactation cycle. In this study, we detected milk yield and composition in a lactation cycle. For the first time, we constructed a gene coexpression network using WGCNA on the basis of 18 gene expression profiles during six stages of a lactation cycle by transcriptome sequencing, generating 10 specific modules. Genes in each module were performed with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Module–trait relationship analysis showed a series of potential candidates related to milk yield and composition. The current study provides an important theoretical basis for the further molecular breeding of dairy cows.

Abstract

Weighted gene coexpression network analysis (WGCNA) is a novel approach that can quickly analyze the relationships between genes and traits. In this study, the milk yield, lactose, fat, and protein of Holstein dairy cows were detected in a lactation cycle. Meanwhile, a total of 18 gene expression profiles were detected using mammary glands from six lactation stages (day 7 to calving, −7 d; day 30 post-calving, 30 d; day 90 post-calving, 90 d; day 180 post-calving, 180 d; day 270 post-calving, 270 d; day 315 post-calving, 315 d). On the basis of the 18 profiles, WGCNA identified for the first time 10 significant modules that may be related to lactation stage, milk yield, and the main milk composition content. Genes in the 10 significant modules were examined with gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The results revealed that the galactose metabolism pathway was a potential candidate for milk yield and milk lactose synthesis. In −7 d, ion transportation was more frequent and cell proliferation related terms became active. In late lactation, the suppressor of cytokine signaling 3 (SOCS3) might play a role in apoptosis. The sphingolipid signaling pathway was a potential candidate for milk fat synthesis. Dairy cows at 315 d were in a period of cell proliferation. Another notable phenomenon was that nonlactating dairy cows had a more regular circadian rhythm after a cycle of lactation. The results provide an important theoretical basis for the further molecular breeding of dairy cows.

Dual SMAD inhibition and Wnt inhibition enable efficient and reproducible differentiations of induced pluripotent stem cells into retinal ganglion cells

Glaucoma is a group of progressive optic neuropathies that share common biological and clinical characteristics including irreversible changes to the optic nerve and visual field loss caused by the death of retinal ganglion cells (RGCs). The loss of RGCs manifests as characteristic cupping or optic nerve degeneration, resulting in visual field loss in patients with Glaucoma. Published studies on in vitro RGC differentiation from stem cells utilized classical RGC signaling pathways mimicking retinal development in vivo. Although many strategies allowed for the generation of RGCs, increased variability between experiments and lower yield hampered the cross comparison between individual lines and between experiments. To address this critical need, we developed a reproducible chemically defined in vitro methodology for generating retinal progenitor cell (RPC) populations from iPSCs, that are efficiently directed towards RGC lineage. Using this method, we reproducibly differentiated iPSCs into RGCs with greater than 80% purity, without any genetic modifications. We used small molecules and peptide modulators to inhibit BMP, TGF-β (SMAD), and canonical Wnt pathways that reduced variability between iPSC lines and yielded functional and mature iPSC-RGCs. Using CD90.2 antibody and Magnetic Activated Cell Sorter (MACS) technique, we successfully purified Thy-1 positive RGCs with nearly 95% purity.

Co-Expression Network Analysis of Spleen Transcriptome in Rock Bream (Oplegnathus fasciatus) Naturally Infected with Rock Bream Iridovirus (RBIV)

Rock bream iridovirus (RBIV) is a notorious agent that causes high mortality in aquaculture of rock bream (Oplegnathus fasciatus). Despite severity of this virus, no transcriptomic studies on RBIV-infected rock bream that can provide fundamental information on protective mechanism against the virus have been reported so far. This study aimed to investigate physiological mechanisms between host and RBIV through transcriptomic changes in the spleen based on RNA-seq. Depending on infection intensity and sampling time point, fish were divided into five groups: uninfected healthy fish at week 0 as control (0C), heavy infected fish at week 0 (0H), heavy mixed RBIV and bacterial infected fish at week 0 (0MH), uninfected healthy fish at week 3 (3C), and light infected fish at week 3 (3L). We explored clusters from 35,861 genes with Fragments Per Kilo-base of exon per Million mapped fragments (FPKM) values of 0.01 or more through signed co-expression network analysis using WGCNA package. Nine of 22 modules were highly correlated with viral infection (|gene significance (GS) vs. module membership (MM) |> 0.5, p-value < 0.05). Expression patterns in selected modules were divided into two: heavy infected (0H and 0MH) and control and light-infected groups (0C, 3C, and 3L). In functional analysis, genes in two positive modules (5448 unigenes) were enriched in cell cycle, DNA replication, transcription, and translation, and increased glycolysis activity. Seven negative modules (3517 unigenes) built in this study showed significant decreases in the expression of genes in lymphocyte-mediated immune system, antigen presentation, and platelet activation, whereas there was significant increased expression of endogenous apoptosis-related genes. These changes lead to RBIV proliferation and failure of host defense, and suggests the importance of blood cells such as thrombocytes and B cells in rock bream in RBIV infection. Interestingly, a hub gene, pre-mRNA processing factor 19 (PRPF19) showing high connectivity (kME), and expression of this gene using qRT-PCR was increased in rock bream blood cells shortly after RBIV was added. It might be a potential biomarker for diagnosis and vaccine studies in rock bream against RBIV. This transcriptome approach and our findings provide new insight into the understanding of global rock bream-RBIV interactions including immune and pathogenesis mechanisms.

Mycobacterial metabolic characteristics in a water meter biofilm revealed by metagenomics and metatranscriptomics

Mycobacteria represent one of the most persistent bacterial populations in drinking water distribution system (DWDS) biofilm communities; however, mycobacterial in situ metabolic profiles are largely unknown. In this study, the metabolic characteristics of mycobacteria in a household water meter biofilm were unveiled using a coupled metagenomic/metatranscriptomic approach. The water meter biofilm appeared to express nitrogenase genes (nifDKH) and a full complement of genes coding for several carbon-fixation pathways, especially the Calvin cycle, suggesting the CO₂ sequestration and dinitrogen fixation potential of the biofilm. These findings indicate that it may be difficult to prevent the formation of DWDS biofilms simply by controlling the availability of organic carbon or nitrogen. The composite genome of mycobacteria (CG-M) was reconstructed based on the obtained omics data. CG-M shared similar genome phylogeny and virulence-factor profiles with Mycobacterium avium complex, suggesting that population CG-M might represent a member of mycobacteria with pathogenicity. According to the gene expression patterns, population CG-M showed the metabolic potential to assimilate CO₂ via the Calvin cycle and/or anaplerotic reactions, and even to grow autotrophically with CO as the sole carbon and energy source. This suggests that organic carbon may not be a limiting factor for mycobacterial growth in DWDSs. Moreover, our results suggest that mycobacterial aromatic degradation is primarily achieved through the catechol meta-cleavage pathway, and biofilm mycobacteria could prefer phosphate as the phosphorus source.

Expression and activity of the calcitonin receptor family in a sample of primary human high-grade gliomas

Background

Glioblastoma (GBM) is the most common and aggressive type of primary brain cancer. With median survival of less than 15 months, identification and validation of new GBM therapeutic targets is of critical importance.

Results

In this study we tested expression and performed pharmacological characterization of the calcitonin receptor (CTR) as well as other members of the calcitonin family of receptors in high-grade glioma (HGG) cell lines derived from individual patient tumours, cultured in defined conditions.

Previous immunohistochemical data demonstrated CTR expression in GBM biopsies and we were able to confirm CALCR (gene encoding CTR) expression. However, as assessed by cAMP accumulation assay, only one of the studied cell lines expressed functional CTR, while the other cell lines have functional CGRP (CLR/RAMP1) receptors. The only CTR-expressing cell line (SB2b) showed modest coupling to the cAMP pathway and no activation of other known CTR signaling pathways, including ERK_1/2 and p38 MAP kinases, and Ca²⁺ mobilization, supportive of low cell surface receptor expression.

Exome sequencing data failed to account for the discrepancy between functional data and expression on the cell lines that do not respond to calcitonin(s) with no deleterious non-synonymous polymorphisms detected, suggesting that other factors may be at play, such as alternative splicing or rapid constitutive receptor internalisation.

Conclusions

This study shows that GPCR signaling can display significant variation depending on cellular system used, and effects seen in model recombinant cell lines or tumour cell lines are not always reproduced in a more physiologically relevant system and vice versa.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5369-y) contains supplementary material, which is available to authorized users.

Identification and characterization of SSR, SNP and InDel molecular markers from RNA-Seq data of guar (Cyamopsis tetragonoloba, L. Taub.) roots

Background

Guar [Cyamopsis tetragonoloba, L. Taub.] is an important industrial crop because of the commercial applications of the galactomannan gum contained in its seeds. Plant breeding programmes based on marker-assisted selection require a rich resource of molecular markers. As limited numbers of such markers are available for guar, molecular breeding programmes have not been undertaken for the genetic improvement of this important crop. Hence, the present work was done to enrich the molecular markers resource of guar by identifying high quality SSR, SNP and InDel markers from the RNA-Seq data of the roots of two guar varieties.

Results

We carried out RNA-Seq analysis of the roots of two guar varieties, namely, RGC-1066 and M-83. A total of 102,479 unigenes with an average length of 1016 bp were assembled from about 30 million high quality pair-end reads generated by an Illumina HiSeq 2500 platform. The assembled unigenes had 86.55% complete and 97.71% partially conserved eukaryotic genes (CEGs). The functional annotation of assembled unigenes using BLASTX against six databases showed that the guar unigenes were most similar to Glycine max. We could assign GO terms to 45,200 unigenes using the UniProt database. The screening of 102,479 unigenes with MISA and SAMtools version 1.4 softwares resulted in the identification of 25,040 high-confidence molecular markers which consisted of 18,792 SSRs, 5999 SNPs and 249 InDels. These markers tagged most of the genes involved in root development, stress tolerance and other general metabolic activities. Each of the 25,040 molecular markers was characterized, particularly with respect to its position in the unigene. For 71% of the molecular markers, we could determine the names, products and functions of the unigenes. About 80% of the markers, from a random sample of molecular markers, showed PCR amplification.

Conclusions

We have identified and characterized 25,040 high confidence SSR, SNP and InDel molecular markers in guar. It is expected that these markers will be useful in molecular breeding programmes and will also be helpful in studying molecular mechanisms of root development, stress tolerance and gum synthesis in guar.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5205-9) contains supplementary material, which is available to authorized users.

Superheat: An R package for creating beautiful and extendable heatmaps for visualizing complex data

The technological advancements of the modern era have enabled the collection of huge amounts of data in science and beyond. Extracting useful information from such massive datasets is an ongoing challenge as traditional data visualization tools typically do not scale well in high-dimensional settings. An existing visualization technique that is particularly well suited to visualizing large datasets is the heatmap. Although heatmaps are extremely popular in fields such as bioinformatics, they remain a severely underutilized visualization tool in modern data analysis. This paper introduces superheat, a new R package that provides an extremely flexible and customizable platform for visualizing complex datasets. Superheat produces attractive and extendable heatmaps to which the user can add a response variable as a scatterplot, model results as boxplots, correlation information as barplots, and more. The goal of this paper is two-fold: (1) to demonstrate the potential of the heatmap as a core visualization method for a range of data types, and (2) to highlight the customizability and ease of implementation of the superheat R package for creating beautiful and extendable heatmaps. The capabilities and fundamental applicability of the superheat package will be explored via three reproducible case studies, each based on publicly available data sources.

Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes

Retinal ganglion cells (RGCs) convey the major output of information collected from the eye to the brain. Thirty subtypes of RGCs have been identified to date. Here, we analyze 6225 RGCs (average of 5000 genes per cell) from right and left eyes by single-cell RNA-seq and classify them into 40 subtypes using clustering algorithms. We identify additional subtypes and markers, as well as transcription factors predicted to cooperate in specifying RGC subtypes. Zic1, a marker of the right eye-enriched subtype, is validated by immunostaining in situ. Runx1 and Fst, the markers of other subtypes, are validated in purified RGCs by fluorescent in situ hybridization (FISH) and immunostaining. We show the extent of gene expression variability needed for subtype segregation, and we show a hierarchy in diversification from a cell-type population to subtypes. Finally, we present a website for comparing the gene expression of RGC subtypes.

Replication Study: Transcriptional amplification in tumor cells with elevated c-Myc

As part of the Reproducibility Project: Cancer Biology, we published a Registered Report (Blum et al., 2015), that described how we intended to replicate selected experiments from the paper 'Transcriptional amplification in tumor cells with elevated c-Myc' (Lin et al., 2012). Here we report the results. We found overexpression of c-Myc increased total levels of RNA in P493-6 Burkitt's lymphoma cells; however, while the effect was in the same direction as the original study (Figure 3E; Lin et al., 2012), statistical significance and the size of the effect varied between the original study and the two different lots of serum tested in this replication. Digital gene expression analysis for a set of genes was also performed on P493-6 cells before and after c-Myc overexpression. Transcripts from genes that were active before c-Myc induction increased in expression following c-Myc overexpression, similar to the original study (Figure 3F; Lin et al., 2012). Transcripts from genes that were silent before c-Myc induction also increased in expression following c-Myc overexpression, while the original study concluded elevated c-Myc had no effect on silent genes (Figure 3F; Lin et al., 2012). Treating the data as paired, we found a statistically significant increase in gene expression for both active and silent genes upon c-Myc induction, with the change in gene expression greater for active genes compared to silent genes. Finally, we report meta-analyses for each result.

RNA-Seq Mouse Brain Regions Expression Data Analysis: Focus on ApoE Functional Network

ApoE expression status was proved to be a highly specific marker of energy metabolism rate in the brain. Along with its neighbor, Translocase of Outer Mitochondrial Membrane 40 kDa (TOMM40) which is involved in mitochondrial metabolism, the corresponding genomic region constitutes the neuroenergetic hotspot. Using RNA-Seq data from a murine model of chronic stress a significant positive expression coordination of seven neighboring genes in ApoE locus in five brain regions was observed. ApoE maintains one of the highest absolute expression values genome-wide, implying that ApoE can be the driver of the neighboring gene expression alteration observed under stressful loads. Notably, we revealed the highly statistically significant increase of ApoE expression in the hypothalamus of chronically aggressive (FDR < 0.007) and defeated (FDR < 0.001) mice compared to the control. Correlation analysis revealed a close association of ApoE and proopiomelanocortin (Pomc) gene expression profiles implying the putative neuroendocrine stress response background of ApoE expression elevation therein.

Prediction of Chromatin Accessibility in Gene-Regulatory Regions from Transcriptomics Data

The epigenetics landscape of cells plays a key role in the establishment of cell-type specific gene expression programs characteristic of different cellular phenotypes. Different experimental procedures have been developed to obtain insights into the accessible chromatin landscape including DNase-seq, FAIRE-seq and ATAC-seq. However, current downstream computational tools fail to reliably determine regulatory region accessibility from the analysis of these experimental data. In particular, currently available peak calling algorithms are very sensitive to their parameter settings and show highly heterogeneous results, which hampers a trustworthy identification of accessible chromatin regions. Here, we present a novel method that predicts accessible and, more importantly, inaccessible gene-regulatory chromatin regions solely relying on transcriptomics data, which complements and improves the results of currently available computational methods for chromatin accessibility assays. We trained a hierarchical classification tree model on publicly available transcriptomics and DNase-seq data and assessed the predictive power of the model in six gold standard datasets. Our method increases precision and recall compared to traditional peak calling algorithms, while its usage is not limited to the prediction of accessible and inaccessible gene-regulatory chromatin regions, but constitutes a helpful tool for optimizing the parameter settings of peak calling methods in a cell type specific manner.