An iteration normalization and test method for differential expression analysis of RNA-seq data

High quality genome annotation and expression visualisation of a mupirocin-producing bacterium

Pseudomonas strain NCIMB10586, in the P. fluorescens subgroup, produces the polyketide antibiotic mupirocin, and has potential as a host for industrial production of a range of valuable products. To underpin further studies on its genetics and physiology, we have used a combination of standard and atypical approaches to achieve a quality of the genome sequence and annotation, above current standards for automated pathways. Assembly of Illumina reads to a PacBio genome sequence created a retrospectively hybrid assembly, identifying and fixing 415 sequencing errors which would otherwise affect almost 5% of annotated coding regions. Our annotation pipeline combined automation based on related well-annotated genomes and stringent, partially manual, tests for functional features. The strain was close to P. synxantha and P. libaniensis and was found to be highly similar to a strain being developed as a weed-pest control agent in Canada. Since mupirocin is a secondary metabolite whose production is switched on late in exponential phase, we carried out RNAseq analysis over an 18 h growth period and have developed a method to normalise RNAseq samples as a group, rather than pair-wise. To review such data we have developed an easily interpreted way to present the expression profiles across a region, or the whole genome at a glance. At the 2-hour granularity of our time-course, the mupirocin cluster increases in expression as an essentially uniform bloc, although the mupirocin resistance gene stands out as being expressed at all the time points.

RALBP1 in Oxidative Stress and Mitochondrial Dysfunction in Alzheimer's Disease

The purpose of our study is to understand the role of the RALBP1 gene in oxidative stress (OS), mitochondrial dysfunction and cognition in Alzheimer's disease (AD) pathogenesis. The RALPB1 gene encodes the 76 kDa protein RLIP76 (Rlip). Rlip functions as a stress-responsive/protective transporter of glutathione conjugates (GS-E) and xenobiotic toxins. We hypothesized that Rlip may play an important role in maintaining cognitive function. The aim of this study is to determine whether Rlip deficiency in mice is associated with AD-like cognitive and mitochondrial dysfunction. Brain tissue obtained from cohorts of wildtype (WT) and Rlip+/- mice were analyzed for OS markers, expression of genes that regulate mitochondrial fission/fusion, and synaptic integrity. We also examined mitochondrial ultrastructure in brains obtained from these mice and further analyzed the impact of Rlip deficiency on gene networks of AD, aging, stress response, mitochondrial function, and CREB signaling. Our studies revealed a significant increase in the levels of OS markers and alterations in the expression of genes and proteins involved in mitochondrial biogenesis, dynamics and synapses in brain tissues from these mice. Furthermore, we compared the cognitive function of WT and Rlip+/- mice. Behavioral, basic motor and sensory function tests in Rlip+/- mice revealed cognitive decline, similar to AD. Gene network analysis indicated dysregulation of stress-activated gene expression, mitochondrial function and CREB signaling genes in the Rlip+/- mouse brain. Our results suggest that Rlip deficiency-associated increases in OS and mitochondrial dysfunction could contribute to the development or progression of OS-related AD processes.

The acute phase protein lactoferrin is a key feature of Alzheimer's disease and predictor of Aβ burden through induction of APP amyloidogenic processing

Amyloidogenic processing of the amyloid precursor protein (APP) forms the amyloid-β peptide (Aβ) component of pathognomonic extracellular plaques of AD. Additional early cortical changes in AD include neuroinflammation and elevated iron levels. Activation of the innate immune system in the brain is a neuroprotective response to infection; however, persistent neuroinflammation is linked to AD neuropathology by uncertain mechanisms. Non-parametric machine learning analysis on transcriptomic data from a large neuropathologically characterised patient cohort revealed the acute phase protein lactoferrin (Lf) as the key predictor of amyloid pathology. In vitro studies showed that an interaction between APP and the iron-bound form of Lf secreted from activated microglia diverted neuronal APP endocytosis from the canonical clathrin-dependent pathway to one requiring ADP ribosylation factor 6 trafficking. By rerouting APP recycling to the Rab11-positive compartment for amyloidogenic processing, Lf dramatically increased neuronal Aβ production. Lf emerges as a novel pharmacological target for AD that not only modulates APP processing but provides a link between Aβ production, neuroinflammation and iron dysregulation.

Higher plane of nutrition pre-weaning enhances Holstein calf mammary gland development through alterations in the parenchyma and fat pad transcriptome

Background

To reduce costs of rearing replacement heifers, researchers have focused on decreasing age at breeding and first calving. To increase returns upon initiation of lactation the focus has been on increasing mammary development prior to onset of first lactation. Enhanced plane of nutrition pre-weaning may benefit the entire replacement heifer operation by promoting mammary gland development and greater future production.

Methods

Twelve Holstein heifer calves (< 1 week old) were reared on 1 of 2 dietary treatments (n = 6/group) for 8 weeks: a control group fed a restricted milk replacer at 0.45 kg/d (R, 20% crude protein, 20% fat), or an accelerated group fed an enhanced milk replacer at 1.13 kg/d (EH, 28% crude protein, 25% fat). At weaning (8 weeks), calves were euthanized and sub-samples of mammary parenchyma (PAR) and mammary fat pad (MFP) were harvested upon removal from the body. Total RNA from both tissues was extracted and sequenced using the Illumina HiSeq2500 platform. The Dynamic Impact Approach (DIA) and Ingenuity Pathway Analysis (IPA) were used for pathway analysis and functions, gene networks, and cross-talk analyses of the two tissues.

Results

When comparing EH vs R 1561 genes (895 upregulated, 666 downregulated) and 970 genes (506 upregulated, 464 downregulated) were differentially expressed in PAR and MFP, respectively. DIA and IPA results highlight a greater proliferation and differentiation activity in both PAR and MFP, supported by an increased metabolic activity. When calves were fed EH, the PAR displayed transcriptional signs of greater overall organ development, with higher ductal growth and branching, together with a supportive blood vessel and nerve network. These activities were mediated by intracellular cascades, such as AKT, SHH, MAPK, and Wnt, probably activated by hormones, growth factors, and endogenous molecules. The analysis also revealed strong communication between MFP and PAR.

Conclusion

The transcriptomics and bioinformatics approach highlighted key mechanisms that mediate the mammary gland response to a higher plane of nutrition in the pre-weaning period.

Electronic supplementary material

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

A Hypothesis Testing Based Method for Normalization and Differential Expression Analysis of RNA-Seq Data

Next-generation sequencing technologies have made RNA sequencing (RNA-seq) a popular choice for measuring gene expression level. To reduce the noise of gene expression measures and compare them between several conditions or samples, normalization is an essential step to adjust for varying sample sequencing depths and other unwanted technical effects. In this paper, we develop a novel global scaling normalization method by employing the available knowledge of housekeeping genes. We formulate the problem from the hypothesis testing perspective and find an optimal scaling factor that minimizes the deviation between the empirical and the nominal type I error. Applying our approach to various simulation studies and real examples, we demonstrate that it is more accurate and robust than the state-of-the-art alternatives in detecting differentially expression genes.

A Human-Specific α7-Nicotinic Acetylcholine Receptor Gene in Human Leukocytes: Identification, Regulation and the Consequences of CHRFAM7A Expression

The human genome contains a variant form of the α7-nicotinic acetylcholine receptor (α7nAChR) gene that is uniquely human. This CHRFAM7A gene arose during human speciation and recent data suggests that its expression alters ligand tropism of the normally homopentameric human α7-AChR ligand-gated cell surface ion channel that is found on the surface of many different cell types. To understand its possible significance in regulating inflammation in humans, we investigated its expression in normal human leukocytes and leukocyte cell lines, compared CHRFAM7A expression to that of the CHRNA7 gene, mapped its promoter and characterized the effects of stable CHRFAM7A overexpression. We report here that CHRFAM7A is highly expressed in human leukocytes but that the levels of both CHRFAM7A and CHRNA7 mRNAs were independent and varied widely. To this end, mapping of the CHRFAM7A promoter in its 5'-untranslated region (UTR) identified a unique 1-kb sequence that independently regulates CHRFAM7A gene expression. Because overexpression of CHRFAM7A in THP1 cells altered the cell phenotype and modified the expression of genes associated with focal adhesion (for example, FAK, P13K, Akt, rho, GEF, Elk1, CycD), leukocyte transepithelial migration (Nox, ITG, MMPs, PKC) and cancer (kit, kitL, ras, cFos cyclinD1, Frizzled and GPCR), we conclude that CHRFAM7A is biologically active. Most surprisingly however, stable CHRFAM7A overexpression in THP1 cells upregulated CHRNA7, which, in turn, led to increased binding of the specific α7nAChR ligand, bungarotoxin, on the THP1 cell surface. Taken together, these data confirm the close association between CHRFAM7A and CHRNA7 expression, establish a biological consequence to CHRFAM7A expression in human leukocytes and support the possibility that this human-specific gene might contribute to, and/or gauge, a human-specific response to inflammation.