DAFS: a data-adaptive flag method for RNA-sequencing data to differentiate genes with low and high expression

The role of heterochronic gene expression and regulatory architecture in early developmental divergence

New developmental programs can evolve through adaptive changes to gene expression. The annelid Streblospio benedicti has a developmental dimorphism, which provides a unique intraspecific framework for understanding the earliest genetic changes that take place during developmental divergence. Using comparative RNAseq through ontogeny, we find that only a small proportion of genes are differentially expressed at any time, despite major differences in larval development and life history. These genes shift expression profiles across morphs by either turning off any expression in one morph or changing the timing or amount of gene expression. We directly connect the contributions of these mechanisms to differences in developmental processes. We examine F1 offspring - using reciprocal crosses - to determine maternal mRNA inheritance and the regulatory architecture of gene expression. These results highlight the importance of both novel gene expression and heterochronic shifts in developmental evolution, as well as the trans-acting regulatory factors in initiating divergence.

The BulkECexplorer compiles endothelial bulk transcriptomes to predict functional versus leaky transcription

Transcriptomic data can be mined to understand the molecular activity of cell types. Yet, functional genes may remain undetected in RNA sequencing (RNA-seq) experiments for technical reasons, such as insufficient read depth or gene dropout. Conversely, RNA-seq experiments may detect lowly expressed mRNAs thought to be biologically irrelevant products of leaky transcription. To represent a cell type's functional transcriptome more accurately, we propose compiling many bulk RNA-seq datasets into a compendium and applying established classification models to predict whether detected transcripts are likely products of active or leaky transcription. Here, we present the BulkECexplorer (bulk RNA-seq endothelial cell explorer) compendium of 240 bulk RNA-seq datasets from five vascular endothelial cell subtypes. This resource reports transcript counts for genes of interest and predicts whether detected transcripts are likely the products of active or leaky gene expression. Beyond its usefulness for vascular biology research, this resource provides a blueprint for developing analogous tools for other cell types.

Characterizing microglial gene expression in a model of secondary progressive multiple sclerosis

Multiple sclerosis (MS) is the most common inflammatory, demyelinating and neurodegenerative disease of the central nervous system in young adults. Chronic-relapsing experimental autoimmune encephalomyelitis (crEAE) in Biozzi ABH mice is an experimental model of MS. This crEAE model is characterized by an acute phase with severe neurological disability, followed by remission of disease, relapse of neurological disease and remission that eventually results in a chronic progressive phase that mimics the secondary progressive phase (SPEAE) of MS. In both MS and SPEAE, the role of microglia is poorly defined. We used a crEAE model to characterize microglia in the different phases of crEAE phases using morphometric and RNA sequencing analyses. At the initial, acute inflammation phase, microglia acquired a pro-inflammatory phenotype. At the remission phase, expression of standard immune activation genes was decreased while expression of genes associated with lipid metabolism and tissue remodeling were increased. Chronic phase microglia partially regain inflammatory gene sets and increase expression of genes associated with proliferation. Together, the data presented here indicate that microglia obtain different features at different stages of crEAE and a particularly mixed phenotype in the chronic stage. Understanding the properties of microglia that are present at the chronic phase of EAE will help to understand the role of microglia in secondary progressive MS, to better aid the development of therapies for this phase of the disease.

Single-chromosome dynamics reveals locus-dependent dynamics and chromosome territory orientation

Dynamic chromatin organization instantly influences DNA accessibility through modulating local macromolecular density and interactions, driving changes in transcription activities. Chromatin dynamics have been reported to be locally confined but contribute to coherent chromatin motion across the entire nucleus. However, the regulation of dynamics, nuclear orientation and compaction of subregions along a single chromosome are not well-understood. We used CRISPR-based real-time single-particle tracking and polymer models to characterize the dynamics of specific genomic loci and determine compaction levels of large human chromosomal domains. Our studies showed that chromosome compaction changed during interphase and that compactions of two arms on chromosome 19 were different. The dynamics of genomic loci were subdiffusive and dependent on chromosome regions and transcription states. Surprisingly, the correlation between locus-dependent nuclear localization and mobility was negligible. Strong tethering interactions detected at the pericentromeric region implies local condensation or associations with organelles within local nuclear microenvironments, such as chromatin-nuclear body association. Based on our findings, we propose a 'guided radial model' for the nuclear orientation of the long arm of chromosome 19.

Brain macrophages acquire distinct transcriptomes in multiple sclerosis lesions and normal appearing white matter

Multiple sclerosis (MS) is a disease of the central nervous system that is characterized by inflammation and focal areas of demyelination, ultimately resulting in axonal degradation and neuronal loss. Several lines of evidence point towards a role for microglia and other brain macrophages in disease initiation and progression, but exactly how lesion formation is triggered is currently unknown. Here, we characterized early changes in MS brain tissue through transcriptomic analysis of normal appearing white matter (NAWM). We found that NAWM was characterized by enriched expression of genes associated with inflammation and cellular stress derived from brain macrophages. Single cell RNA sequencing confirmed a stress response in brain macrophages in NAWM and identified specific microglia and macrophage subsets at different stages of demyelinating lesions. We identified both phagocytic/activated microglia and CAM clusters that were associated with various MS lesion types. These overall changes in microglia and macrophages associated with lesion development in MS brain tissue may provide therapeutic targets to limit lesion progression and demyelination.

Identification of distinct and age-dependent p16High microglia subtypes

Cells expressing high levels of the cyclin‐dependent kinase (CDK)4/6 inhibitor p16 (p16^High) accumulate in aging tissues and promote multiple age‐related pathologies, including neurodegeneration. Here, we show that the number of p16^High cells is significantly increased in the central nervous system (CNS) of 2‐year‐old mice. Bulk RNAseq indicated that genes expressed by p16^High cells were associated with inflammation and phagocytosis. Single‐cell RNAseq of brain cells indicated p16^High cells were primarily microglia, and their accumulation was confirmed in brains of aged humans. Interestingly, we identified two distinct subpopulations of p16^High microglia in the mouse brain, with one being age‐associated and one present in young animals. Both p16^High clusters significantly differed from previously described disease‐associated microglia and expressed only a partial senescence signature. Taken together, our study provides evidence for the existence of two p16‐expressing microglia populations, one accumulating with age and another already present in youth that could positively and negatively contribute to brain homeostasis, function, and disease.

Profiling Microglia From Alzheimer's Disease Donors and Non-demented Elderly in Acute Human Postmortem Cortical Tissue

Microglia are the tissue-resident macrophages of the central nervous system (CNS). Recent studies based on bulk and single-cell RNA sequencing in mice indicate high relevance of microglia with respect to risk genes and neuro-inflammation in Alzheimer's disease (AD). Here, we investigated microglia transcriptomes at bulk and single-cell levels in non-demented elderly and AD donors using acute human postmortem cortical brain samples. We identified seven human microglial subpopulations with heterogeneity in gene expression. Notably, gene expression profiles and subcluster composition of microglia did not differ between AD donors and non-demented elderly in bulk RNA sequencing nor in single-cell sequencing.

Custom selected reference genes outperform pre-defined reference genes in transcriptomic analysis

BACKGROUND

RNA sequencing allows the measuring of gene expression at a resolution unmet by expression arrays or RT-qPCR. It is however necessary to normalize sequencing data by library size, transcript size and composition, among other factors, before comparing expression levels. The use of internal control genes or spike-ins is advocated in the literature for scaling read counts, but the methods for choosing reference genes are mostly targeted at RT-qPCR studies and require a set of pre-selected candidate controls or pre-selected target genes.

RESULTS

Here, we report an R-based pipeline to select internal control genes based solely on read counts and gene sizes. This novel method first normalizes the read counts to Transcripts per Million (TPM) and then excludes weakly expressed genes using the DAFS script to calculate the cut-off. It then selects as references the genes with lowest TPM covariance. We used this method to pick custom reference genes for the differential expression analysis of three transcriptome sets from transgenic Arabidopsis plants expressing heterologous fungal effector proteins tagged with GFP (using GFP alone as the control). The custom reference genes showed lower covariance and fold change as well as a broader range of expression levels than commonly used reference genes. When analyzed with NormFinder, both typical and custom reference genes were considered suitable internal controls, but the custom selected genes were more stably expressed. geNorm produced a similar result in which most custom selected genes ranked higher (i.e. were more stably expressed) than commonly used reference genes.

CONCLUSIONS

The proposed method is innovative, rapid and simple. Since it does not depend on genome annotation, it can be used with any organism, and does not require pre-selected reference candidates or target genes that are not always available.

Urinary Sediment Transcriptomic and Longitudinal Data to Investigate Renal Function Decline in Type 1 Diabetes

Introduction: Using a discovery/validation approach we investigated associations between a panel of genes selected from a transcriptomic study and the estimated glomerular filtration rate (eGFR) decline across time in a cohort of type 1 diabetes (T1D) patients. Experimental: Urinary sediment transcriptomic was performed to select highly modulated genes in T1D patients with rapid eGFR decline (decliners) vs. patients with stable eGFR (non-decliners). The selected genes were validated in samples from a T1D cohort (n = 54, mean diabetes duration of 21 years, 61% women) followed longitudinally for a median of 12 years in a Diabetes Outpatient Clinic. Results: In the discovery phase, the transcriptomic study revealed 158 genes significantly different between decliners and non-decliners. Ten genes increasingly up or down-regulated according to renal function worsening were selected for validation by qRT-PCR; the genes CYP4F22, and PMP22 were confirmed as differentially expressed comparing decliners vs. non-decliners after adjustment for potential confounders. CYP4F22, LYPD3, PMP22, MAP1LC3C, HS3ST2, GPNMB, CDH6, and PKD2L1 significantly modified the slope of eGFR in T1D patients across time. Conclusions: Eight genes identified as differentially expressed in the urinary sediment of T1D patients presenting different eGFR decline rates significantly increased the accuracy of predicted renal function across time in the studied cohort. These genes may be a promising way of unveiling novel mechanisms associated with diabetic kidney disease progression.

The novel high-frequency variant of TRPV3 p.A628T in East Asians showing faster sensitization in response to chemical agonists

TRPV3, a member of the thermosensitive Ca²⁺-permeable TRPV channel subfamily expressed in skin and sensory nerves, is also activated by chemical agonists such as 2-aminoethyl diphenylborinate (2-APB). Repetitive stimuli induce sensitization of TRPV3 activation, characterized by the cumulative increase in current amplitude and linearization of current-voltage relation (I/V curve). Through genomic analysis of various populations, we found non-rare TRPV3 mutation (p.A628T) in East Asian people with an allele frequency of 0.249 while 0.007 in Caucasian. Slope conductance of unitary channel was not different between WT and p.A628T. Whole-cell patch clamp study of wildtype TRPV3 (WT) and p.A628T overexpressed in HEK293T cells showed similar sensitization by the repetitive increase in temperature from 23 to 37 °C, while slightly higher sensitization to 43 °C in p.A628T. In contrast, the repetitive application of 2-APB (10 μM) or carvacrol (100 μM) induced faster sensitization in p.A628T than WT. However, 1 μM farnesyl pyrophosphate, an intrinsic lipid metabolite agonist, induced similar level of slow activations in WT and p.A628T. In Fura-2 microspectrofluorimetry, the 2-APB pulses induced a faster increase of [Ca²⁺]_c in p.A628T than WT. In terms of ionic selectivity of channels, WT and p.A628T showed similar Ca²⁺ permeability (P_Ca/P_Na) calculated from the reversal potential of I/V curves. Taken together, p.A628T shows faster sensitization to chemical agonists that are reflected as higher [Ca²⁺]_c signaling. Based on the intriguing pharmacological sensitivity, the physiological implications of p.A628T in the East Asian population require further investigation.

Targeted Diet Modification Reduces Multiple Sclerosis-like Disease in Adult Marmoset Monkeys from an Outbred Colony

Experimental autoimmune encephalomyelitis (EAE) in common marmosets is a translationally relevant model of the chronic neurologic disease multiple sclerosis. Following the introduction of a new dietary supplement in our purpose-bred marmoset colony, the percentage of marmosets in which clinically evident EAE could be induced by sensitization against recombinant human myelin oligodendrocyte glycoprotein in IFA decreased from 100 to 65%. The reduced EAE susceptibility after the dietary change coincided with reduced Callitrichine herpesvirus 3 expression in the colony, an EBV-related γ1-herpesvirus associated with EAE. We then investigated, in a controlled study in marmoset twins, which disease-relevant parameters were affected by the dietary change. The selected twins had been raised on the new diet for at least 12 mo prior to the study. In twin siblings reverted to the original diet 8 wk prior to EAE induction, 100% disease prevalence (eight out of eight) was restored, whereas in siblings remaining on the new diet the EAE prevalence was 75% (six out of eight). Spinal cord demyelination, a classical hallmark of the disease, was significantly lower in new-diet monkeys than in monkeys reverted to the original diet. In new-diet monkeys, the proinflammatory T cell response to recombinant human myelin oligodendrocyte glycoprotein was significantly reduced, and RNA-sequencing revealed reduced apoptosis and enhanced myelination in the brain. Systematic typing of the marmoset gut microbiota using 16S rRNA sequencing demonstrated a unique, Bifidobacteria-dominated composition, which changed after disease induction. In conclusion, targeted dietary intervention exerts positive effects on EAE-related parameters in multiple compartments of the marmoset's gut-immune-CNS axis.

Evolution of Sex Chromosome Dosage Compensation in Animals: A Beautiful Theory, Undermined by Facts and Bedeviled by Details

Many animals with genetic sex determination harbor heteromorphic sex chromosomes, where the heterogametic sex has half the gene dose of the homogametic sex. This imbalance, if reflected in the abundance of transcripts or proteins, has the potential to deleteriously disrupt interactions between X-linked and autosomal loci in the heterogametic sex. Classical theory predicts that molecular mechanisms will evolve to provide dosage compensation that recovers expression levels comparable to ancestral expression prior to sex chromosome divergence. Such dosage compensating mechanisms may also, secondarily, result in balanced sex-linked gene expression between males and females. However, numerous recent studies addressing sex chromosome dosage compensation (SCDC) in a diversity of animals have yielded a surprising array of patterns concerning dosage compensation in the heterogametic sex, as well as dosage balance between sexes. These results substantially contradict longstanding theory, catalyzing both novel perspectives and new approaches in dosage compensation research. In this review, we summarize the theory, analytical approaches, and recent results concerning evolutionary patterns of SCDC in animals. We also discuss methodological challenges and discrepancies encountered in this research, which often underlie conflicting results. Finally, we discuss what outstanding questions and opportunities exist for future research on SCDC.

A Method for Cross-Species Visualization and Analysis of RNA-Sequence Data

In this methods article, I describe a computational workflow for cross-species visualization and comparison of mRNA-seq transcriptome profiling data. The workflow is based on gene set variation analysis (GSVA) and is illustrated using commands in the R programming language. I provide a complete step-by-step procedure for the workflow using mRNA-seq data sets from dog and human bladder cancer as an example.

Comparative evaluation of isoform-level gene expression estimation algorithms for RNA-seq and exon-array platforms

Given that the majority of multi-exon genes generate diverse functional products, it is important to evaluate expression at the isoform level. Previous studies have demonstrated strong gene-level correlations between RNA sequencing (RNA-seq) and microarray platforms, but have not studied their concordance at the isoform level. We performed transcript abundance estimation on raw RNA-seq and exon-array expression profiles available for common glioblastoma multiforme samples from The Cancer Genome Atlas using different analysis pipelines, and compared both the isoform- and gene-level expression estimates between programs and platforms. The results showed better concordance between RNA-seq/exon-array and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) platforms for fold change estimates than for raw abundance estimates, suggesting that fold change normalization against a control is an important step for integrating expression data across platforms. Based on RT-qPCR validations, eXpress and Multi-Mapping Bayesian Gene eXpression (MMBGX) programs achieved the best performance for RNA-seq and exon-array platforms, respectively, for deriving the isoform-level fold change values. While eXpress achieved the highest correlation with the RT-qPCR and exon-array (MMBGX) results overall, RSEM was more highly correlated with MMBGX for the subset of transcripts that are highly variable across the samples. eXpress appears to be most successful in discriminating lowly expressed transcripts, but IsoformEx and RSEM correlate more strongly with MMBGX for highly expressed transcripts. The results also reinforce how potentially important isoform-level expression changes can be masked by gene-level estimates, and demonstrate that exon arrays yield comparable results to RNA-seq for evaluating isoform-level expression changes.

N-acetylcysteine Counteracts Adipose Tissue Macrophage Infiltration and Insulin Resistance Elicited by Advanced Glycated Albumin in Healthy Rats

Background: Advanced glycation endproducts elicit inflammation. However, their role in adipocyte macrophage infiltration and in the development of insulin resistance, especially in the absence of the deleterious biochemical pathways that coexist in diabetes mellitus, remains unknown. We investigated the effect of chronic administration of advanced glycated albumin (AGE-albumin) in healthy rats, associated or not with N-acetylcysteine (NAC) treatment, on insulin sensitivity, adipose tissue transcriptome and macrophage infiltration and polarization.

Methods: Male Wistar rats were intraperitoneally injected with control (C) or AGE-albumin alone, or, together with NAC in the drinking water. Biochemical parameters, lipid peroxidation, gene expression and protein contents were, respectively, determined by enzymatic techniques, reactive thiobarbituric acid substances, RT-qPCR and immunohistochemistry or immunoblot. Carboxymethyllysine (CML) and pyrraline (PYR) were determined by LC/mass spectrometry (LC-MS/MS) and ELISA.

Results: CML and PYR were higher in AGE-albumin as compared to C. Food consumption, body weight, systolic blood pressure, plasma lipids, glucose, hepatic and renal function, adipose tissue relative weight and adipocyte number were similar among groups. In AGE-treated animals, insulin resistance, adipose macrophage infiltration and Col12a1 mRNA were increased with no changes in M1 and M2 phenotypes as compared to C-albumin-treated rats. Total GLUT4 content was reduced by AGE-albumin as compared to C-albumin. NAC improved insulin sensitivity, reduced urine TBARS, adipose macrophage number and Itgam and Mrc mRNA and increased Slc2a4 and Ppara. CD11b, CD206, Ager, Ddost, Cd36, Nfkb1, Il6, Tnf, Adipoq, Retn, Arg, and Il12 expressions were similar among groups.

Conclusions: AGE-albumin sensitizes adipose tissue to inflammation due to macrophage infiltration and reduces GLUT4, contributing to insulin resistance in healthy rats. NAC antagonizes AGE-albumin and prevents insulin resistance. Therefore, it may be a useful tool in the prevention of AGE action on insulin resistance and long-term complications of DM.

Transcriptomic analysis of purified human cortical microglia reveals age-associated changes

Microglia are essential for CNS homeostasis and innate neuroimmune function, and play important roles in neurodegeneration and brain aging. Here we present gene expression profiles of purified microglia isolated at autopsy from the parietal cortex of 39 human subjects with intact cognition. Overall, genes expressed by human microglia were similar to those in mouse, including established microglial genes CX3CR1, P2RY12 and ITGAM (CD11B). However, a number of immune genes, not identified as part of the mouse microglial signature, were abundantly expressed in human microglia, including TLR, Fcγ and SIGLEC receptors, as well as TAL1 and IFI16, regulators of proliferation and cell cycle. Age-associated changes in human microglia were enriched for genes involved in cell adhesion, axonal guidance, cell surface receptor expression and actin (dis)assembly. Limited overlap was observed in microglial genes regulated during aging between mice and humans, indicating that human and mouse microglia age differently.

The Non-Linear Child: Ontogeny, Isoniazid Concentration, and NAT2 Genotype Modulate Enzyme Reaction Kinetics and Metabolism

N-acetyltransferase 2 (NAT2) catalyzes the acetylation of isoniazid to N-acetylisoniazid. NAT2 polymorphism explains 88% of isoniazid clearance variability in adults. We examined the effects of clinical and genetic factors on Michaelis-Menten reaction kinetic constants of maximum velocity (V_max) and affinity (K_m) in children 0–10 years old. We measured the rates of isoniazid elimination and N-acetylisoniazid production in the blood of 30 children. Since maturation effects could be non-linear, we utilized a pharmacometric approach and the artificial intelligence method, multivariate adaptive regression splines (MARS), to identify factors predicting NAT2 V_max and K_m by examining clinical, genetic, and laboratory factors in toto. Isoniazid concentration predicted both V_max and K_m and superseded the contribution of NAT2 genotype. Age non-linearly modified the NAT2 genotype contribution until maturation at ≥ 5.3 years. Thus, enzyme efficiency was constrained by substrate concentration, genes, and age. Since MARS output is in the form of basis functions and equations, it allows multiscale systems modeling from the level of cellular chemical reactions to whole body physiological parameters, by automatic selection of significant predictors by the algorithm.

•

We identified the NAT2 K_m and V_max in children treated with isoniazid.

•

Artificial intelligence (AI) algorithms were used to find predictors of K_m and V_max.

•

Isoniazid concentration affected V_max and K_m, and superseded NAT2 genotype.

•

Age non-linearly modified NAT2 genotype contribution until maturation at ≥ 5.3 years.

•

AI output is in the form of equations that allow multiscale systems modeling.

The effects of maturation on drug metabolism have not been studied for the type phase II enzymes such as NAT2, which metabolizes the drug isoniazid. Genes have been found to control speed of isoniazid metabolism. Studies to characterize affinity and maximum velocity for isoniazid metabolism in people were last performed in two individuals' livers in the 1960s. We identified NAT2 affinity and maximum velocity in 30 tuberculosis children treated with isoniazid. Artificial intelligence methods found that metabolism was affected by the drug's concentration more than by genes, which were affected by age up to 5.3 years.

Evaluating the Stability of RNA-Seq Transcriptome Profiles and Drug-Induced Immune-Related Expression Changes in Whole Blood

Methods were developed to evaluate the stability of rat whole blood expression obtained from RNA sequencing (RNA-seq) and assess changes in whole blood transcriptome profiles in experiments replicated over time. Expression was measured in globin-depleted RNA extracted from the whole blood of Sprague-Dawley rats, given either saline (control) or neurotoxic doses of amphetamine (AMPH). The experiment was repeated four times (paired control and AMPH groups) over a 2-year span. The transcriptome of the control and AMPH-treated groups was evaluated on: 1) transcript levels for ribosomal protein subunits; 2) relative expression of immune-related genes; 3) stability of the control transcriptome over 2 years; and 4) stability of the effects of AMPH on immune-related genes over 2 years. All, except one, of the 70 genes that encode the 80s ribosome had levels that ranked in the top 5% of all mean expression levels. Deviations in sequencing performance led to significant changes in the ribosomal transcripts. The overall expression profile of immune-related genes and genes specific to monocytes, T-cells or B-cells were well represented and consistent within treatment groups. There were no differences between the levels of ribosomal transcripts in time-matched control and AMPH groups but significant differences in the expression of immune-related genes between control and AMPH groups. AMPH significantly increased expression of some genes related to monocytes but down-regulated those specific to T-cells. These changes were partially due to changes in the two types of leukocytes present in blood, which indicate an activation of the innate immune system by AMPH. Thus, the stability of RNA-seq whole blood transcriptome can be verified by assessing ribosomal protein subunits and immune-related gene expression. Such stability enables the pooling of samples from replicate experiments to carry out differential expression analysis with acceptable power.