RnaSeqSampleSize: real data based sample size estimation for RNA sequencing

Genome-wide RNA-seq, DNA methylation and small RNA-seq analysis unraveled complex gene regulatory networks in psoriasis pathogenesis

Psoriasis is a complex inflammatory skin disease characterized by reversible albeit relapsing red scaly plaques in the skin of a patient. In addition to the genetic predisposition, involvement of epigenetic and non-coding RNAs have also been liked with the disease. Nevertheless, any comprehensive study involving transcriptomic, small-RNA and DNA methylation at the genomic level from same patients is lacking. To investigate the complex regulation of molecular pathways in psoriasis, we carried out multi-omics integrative analysis of RNA-sequencing, small RNA-sequencing and DNA methylation profiling from the psoriatic and adjacent normal skin tissues. Our multi-omics analysis identified the genes and biological processes regulated either independently or in combination by DNA methylation and microRNAs. We identified miRNAs that specifically regulated keratinocyte hyper-proliferation, and cell cycle progression and checkpoint signaling in psoriasis. On contrary, DNA methylation was found to be more predominant in regulating immune and inflammatory responses, another causative factor in psoriasis pathogenesis. Many characteristic pathways in psoriasis e.g., Th17 cell differentiation and JAK-STAT signaling, were found to be regulated by both miRNAs and DNA methylation. We carried out functional characterization of a downregulated miRNA hsa-let-7c-5p, predicted to target upregulated genes in psoriasis involved in cell cycle processes, Th17 cell differentiation and JAK-STAT signaling pathways. Overexpression of hsa-let-7c-5p in keratinocytes caused the downregulation of its target genes, resulting in reduced cell proliferation and migration rates, demonstrating potential of miRNAs in regulating psoriasis pathogenesis. In conclusion, our findings identified distinct and shared gene-networks regulated by DNA methylation and miRNAs of a complex disease with reversible phenotype.

Endometrial responsiveness to interferon-tau and its association with subsequent reproductive performance in dairy heifers

Our objectives were to evaluate the endometrial responsiveness of dairy heifers to an intrauterine infusion of recombinant bovine interferon-tau (rbIFN-τ) and to associate endometrial responses to rbIFN-τ with subsequent reproductive performance. In experiments 1 and 2, cyclic heifers were enrolled in a program for synchronization of the estrous cycle, and blood sampling and ultrasonography examinations were performed on d 0, 4, 7, 11, and 14 of the estrous cycle. In experiment 1, heifers were randomly assigned to receive an intrauterine infusion containing 2 µg of rbIFN-τ (rbIFN-τ = 19) or saline control (CTRL = 19) into the uterine horn ipsilateral to the corpus luteum (CL) on d 14 of the estrous cycle. Then, 6 hours after the infusion, the infused uterine horn was flushed for sampling of the uterine luminal fluid (ULF) for composition analysis, and the endometrium was biopsied for transcriptomics. In experiment 2, 100 heifers received an intrauterine infusion of rbIFN-τ, and the same procedures for uterine sample collection were performed as described in experiment 1. After the intrauterine test, heifers were enrolled in a breeding program and classified as highly fertile (HF; pregnant at first AI) or subfertile (SF; not pregnant at first AI). Statistical analyses were performed using regression models, which included the effects of treatment (experiment 1: CTRL vs. rbIFN-τ) or fertility group (experiment 2: HF vs. SF) and block of samples. Intrauterine infusion of rbIFN-τ increased the expression of classical interferon-stimulated genes in the endometrium (e.g., ISG15, MX1, OAS2, IRF9, and USP18), and an antiviral response was predicted to be the main downstream effect of the transcriptome changes. In addition, rbIFN-τ increased the abundance of cholesterol, glycerol, and the overall concentration of oxylipins in the ULF. Analysis of endometrial transcriptome between HF and SF heifers revealed important differences in the expression of genes associated with cell signaling, metabolism, attachment, and migration, with a large representation of genes encoding extracellular matrix proteins. In general, differentially expressed genes were expected to be downregulated by IFN-τ but seemed to fail to be downregulated in SF heifers, resulting in higher expression in SF compared with HF heifers. Subfertile heifers had lower concentrations of glycerol and an altered profile of oxylipins in the ULF, with a lower abundance of oxylipins derived from arachidonic acid and dihomo-γ-linolenic acid, and a greater abundance of oxylipins derived from linoleic acid. Measurements of ovarian function did not differ between groups and, therefore, did not influence the observed results in uterine biology. Overall, the endometrial responsiveness to IFN-τ is variable among individuals and associated with subsequent fertility of heifers, indicating that communication between conceptus and endometrium is critical for the uterine receptivity and survival of pregnancy.

Impact of florfenicol dosing regimen on the phenotypic and genotypic resistance of enteric bacteria in steers

The food animal sector's use of antimicrobials is heavily critiqued for its role in allowing resistance to develop against critically important antimicrobials in human health. The WHO recommends using lower tier antimicrobials such as florfenicol for disease treatment. The primary objective of this study was to assess the differences in resistance profiles of enteric microbes following administration of florfenicol to steers using both FDA-approved dosing regimens and two different detection methods. Our hypothesis was that we would identify an increased prevalence of resistance in the steers administered the repeated, lower dose of florfenicol; additionally, we hypothesized resistance profiles would be similar between both detection methods. Twelve steers were administered either two intramuscular (20 mg/kg q 48 h; n = 6) or a single subcutaneous dose (40 mg/kg, n = 6). Fecal samples were collected for 38 days, and E. coli and Enterococcus were isolated and tested for resistance. Fecal samples were submitted for metagenomic sequencing analysis. Metagenomics revealed genes conferring resistance to aminoglycosides as the most abundant drug class. Most multidrug resistance genes contained phenicols. The genotypic and phenotypic patterns of resistance were not similar between drug classes. Observed increases in resistant isolates and relative abundance of resistance genes peaked after drug administration and returned to baseline by the end of the sampling period. The use of a "lower tier" antimicrobial, such as florfenicol, may cause an increased amount of resistance to critically important antimicrobials for a brief period, but these changes largely resolve by the end of the drug withdrawal period.

Polygenic risk for schizophrenia converges on alternative polyadenylation as molecular mechanism underlying synaptic impairment

Schizophrenia (SCZ) is a genetically heterogenous psychiatric disorder of highly polygenic nature. Correlative evidence from genetic studies indicate that the aggregated effects of distinct genetic risk factor combinations found in each patient converge onto common molecular mechanisms. To prove this on a functional level, we employed a reductionistic cellular model system for polygenic risk by differentiating induced pluripotent stem cells (iPSCs) from 104 individuals with high polygenic risk load and controls into cortical glutamatergic neurons (iNs). Multi-omics profiling identified widespread differences in alternative polyadenylation (APA) in the 3' untranslated region of many synaptic transcripts between iNs from SCZ patients and healthy donors. On the cellular level, 3'APA was associated with a reduction in synaptic density of iNs. Importantly, differential APA was largely conserved between postmortem human prefrontal cortex from SCZ patients and healthy donors, and strongly enriched for transcripts related to synapse biology. 3'APA was highly correlated with SCZ polygenic risk and affected genes were significantly enriched for SCZ associated common genetic variation. Integrative functional genomic analysis identified the RNA binding protein and SCZ GWAS risk gene PTBP2 as a critical trans-acting factor mediating 3'APA of synaptic genes in SCZ subjects. Functional characterization of PTBP2 in iNs confirmed its key role in 3'APA of synaptic transcripts and regulation of synapse density. Jointly, our findings show that the aggregated effects of polygenic risk converge on 3'APA as one common molecular mechanism that underlies synaptic impairments in SCZ.

Transcriptomic and metabolomic differences between banana varieties which are resistant or susceptible to Fusarium wilt

Background

Fusarium wilt, caused by Fusarium oxysporum f. sp. cubense race 4 (Foc4), is the most lethal disease of bananas in Asia.

Methods

To better understand the defense response of banana to Fusarium wilt, the transcriptome and metabolome profiles of the roots from resistant and susceptible bananas inoculated with Foc4 were compared.

Results

After Foc4 inoculation, there were 172 and 1,856 differentially expressed genes (DEGs) in the Foc4-susceptible variety (G1) and Foc4-resistant variety (G9), respectively. In addition, a total of 800 DEGs were identified between G1 and G9, which were mainly involved in the oxidation-reduction process, cell wall organization, phenylpropanoid biosynthesis, and lipid and nitrogen metabolism, especially the DEGs of Macma4_08_g22610, Macma4_11_g19760, and Macma4_03_g06480, encoding non-classical arabinogalactan protein; GDSL-like lipase; and peroxidase. In our study, G9 showed a stronger and earlier response to Foc4 than G1. As the results of metabolomics, lipids, phenylpropanoids and polyketides, organic acids, and derivatives played an important function in response to Fusarium wilt. More importantly, Macma4_11_g19760 might be one of the key genes that gave G9 more resistance to Foc4 by a lowered expression and negative regulation of lipid metabolism. This study illustrated the difference between the transcriptomic and metabolomic profiles of resistant and susceptible bananas. These results improved the current understanding of host-pathogen interactions and will contribute to the breeding of resistant banana plants.

Clonal dynamics limits detection of selection in tumour xenograft CRISPR/Cas9 screens

Transplantable in vivo CRISPR/Cas9 knockout screens, in which cells are edited in vitro and inoculated into mice to form tumours, allow evaluation of gene function in a cancer model that incorporates the multicellular interactions of the tumour microenvironment. To improve our understanding of the key parameters for success with this method, we investigated the choice of cell line, mouse host, tumour harvesting timepoint and guide RNA (gRNA) library size. We found that high gRNA (80-95%) representation was maintained in a HCT116 subline transduced with the GeCKOv2 whole-genome gRNA library and transplanted into NSG mice when tumours were harvested at early (14 d) but not late time points (38-43 d). The decreased representation in older tumours was accompanied by large increases in variance in gRNA read counts, with notable expansion of a small number of random clones in each sample. The variable clonal dynamics resulted in a high level of 'noise' that limited the detection of gRNA-based selection. Using simulated datasets derived from our experimental data, we show that considerable reductions in count variance would be achieved with smaller library sizes. Based on our findings, we suggest a pathway to rationally design adequately powered in vivo CRISPR screens for successful evaluation of gene function.

Five families of diverse DNA viruses comprehensively restructure the nucleus

Many viruses have evolved ways to restructure their host cell's nucleus profoundly and unexpectedly upon infection. In particular, DNA viruses that need to commandeer their host's cellular synthetic functions to produce their progeny can induce the condensation and margination of host chromatin during productive infection, a phenomenon known as virus-induced reorganization of cellular chromatin (ROCC). These ROCC-inducing DNA viruses belong to 5 families (herpesviruses, baculoviruses, adenoviruses, parvoviruses, and geminiviruses) that infect a wide range of hosts and are important for human and ecosystem health, as well as for biotechnology. Although the study of virus-induced ROCC is in its infancy, investigations are already raising important questions, such as why only some DNA viruses that replicate their genomes in the nucleus elicit ROCC. Studying the shared and distinct properties of ROCC-inducing viruses will provide valuable insights into viral reorganization of host chromatin that could have implications for future therapies that target the viral life cycle.

Deep RNA-seq of male and female murine sensory neuron subtypes after nerve injury

ABSTRACT

Dorsal root ganglia (DRG) neurons have been well described for their role in driving both acute and chronic pain. Although nerve injury is known to cause transcriptional dysregulation, how this differs across neuronal subtypes and the impact of sex is unclear. Here, we study the deep transcriptional profiles of multiple murine DRG populations in early and late pain states while considering sex. We have exploited currently available transgenics to label numerous subpopulations for fluorescent-activated cell sorting and subsequent transcriptomic analysis. Using bulk tissue samples, we are able to circumvent the issues of low transcript coverage and drop-outs seen with single-cell data sets. This increases our power to detect novel and even subtle changes in gene expression within neuronal subtypes and discuss sexual dimorphism at the neuronal subtype level. We have curated this resource into an accessible database for other researchers ( https://livedataoxford.shinyapps.io/drg-directory/ ). We see both stereotyped and unique subtype signatures in injured states after nerve injury at both an early and late timepoint. Although all populations contribute to a general injury signature, subtype enrichment changes can also be seen. Within populations, there is not a strong intersection of sex and injury, but previously unknown sex differences in naïve states-particularly in Aβ-RA + Aδ-low threshold mechanoreceptors-still contribute to differences in injured neurons.

Single-cell atlas reveals correlates of high cognitive function, dementia, and resilience to Alzheimer's disease pathology

Alzheimer's disease (AD) is the most common cause of dementia worldwide, but the molecular and cellular mechanisms underlying cognitive impairment remain poorly understood. To address this, we generated a single-cell transcriptomic atlas of the aged human prefrontal cortex covering 2.3 million cells from postmortem human brain samples of 427 individuals with varying degrees of AD pathology and cognitive impairment. Our analyses identified AD-pathology-associated alterations shared between excitatory neuron subtypes, revealed a coordinated increase of the cohesin complex and DNA damage response factors in excitatory neurons and in oligodendrocytes, and uncovered genes and pathways associated with high cognitive function, dementia, and resilience to AD pathology. Furthermore, we identified selectively vulnerable somatostatin inhibitory neuron subtypes depleted in AD, discovered two distinct groups of inhibitory neurons that were more abundant in individuals with preserved high cognitive function late in life, and uncovered a link between inhibitory neurons and resilience to AD pathology.

Prospective Blood Transcriptomics Study in a Motor Vehicle Collision Cohort Identified a Protective Function of the SAMD15 Gene Against Chronic Pain

Traumatic brain injuries following motor vehicle collisions (MVCs) are ubiquitous. Surprisingly, there are no correlates between concussion impact force and long-term pain outcomes. To study the molecular underpinnings of chronic pain after MVC, we assembled a prospective cohort of 36 subjects that experienced MVC and suffered documented mild traumatic brain injuries. For each participant, a first blood sample was drawn within 72 hours of the collision, then a second one at the 6-month mark. Pain was also assessed at the second blood draw to determine if pain became chronic or resolved. Blood samples enabled transcriptomics analyses for immune cells. At the transcriptome-wide level, we found that Sterile Alpha Motif Domain Containing 15 (SAMD15) mRNA was significantly upregulated with time in subjects who resolved their pain whereas unregulated in those with persistent pain. Using several large publicly available datasets, such as the UK Biobank and the GTeX portal, we then linked elevated SAMD15 gene expression, elevated neutrophils cell counts, and decreased risk for chronic pain to increased dosage of the T allele at SNP rs4903580, situated within SAMD15's gene locus. The causality between the components of our model was established and supported by Mendelian randomization. Overall, our results support the role of SAMD15 as a potential gene effector for neutrophil-dependent chronic pain development. PERSPECTIVE: This article highlights the potential protective role of the SAMD15 gene against chronic pain following a mild traumatic brain injury. The expression of the gene is associated with a SNP rs4903580, which is itself associated with neutrophils counts as well as chronic pain in large genetic studies.

Molecular analysis of duodenal eosinophilia

BACKGROUND

Eosinophilic duodenitis (EoD), characterized by nonspecific gastrointestinal symptoms and increased numbers of duodenal eosinophils, may be in the eosinophilic gastrointestinal disease spectrum. However, diagnostic thresholds and pathogenic processes of duodenal tissue eosinophilia are inadequately characterized.

OBJECTIVE

We aimed to define an EoD transcriptome and pathologic pathways.

METHODS

RNA sequencing and histologic features of human duodenal biopsy samples were analyzed as a function of duodenal eosinophils levels. For analyses, we defined EoD as more than 52 peak eosinophils/hpf (n = 8), duodenal eosinophilia as 30 to 52 eosinophils/hpf (n = 11), and normal controls as fewer than 30 eosinophils/hpf (n = 8). Associations between gene expression and histologic features were analyzed with Spearman correlation.

RESULTS

We identified 382 differentially expressed genes (EoD transcriptome) between EoD and normal controls (>2-fold change [adjusted P < .05]). The EoD transcriptome distinguished EoD from controls (duodenal eosinophilia and normal controls). The duodenal eosinophil count was correlated with a distinct EoD transcriptome when 50 to 60 peak eosinophils/hpf were present. The EoD transcriptome was enriched in genes involved in IL-4/IL-13 signaling, mast cells, and myeloid progenitor cells. Among duodenal histologic features, lamina propria eosinophil sheets was the most associated with transcriptomic changes (r = 0.66; P < .01). EoD gene signatures were shared with eosinophilic esophagitis and eosinophilic gastritis but not with eosinophilic colitis or celiac disease.

CONCLUSION

We have identified an EoD transcriptomic signature that emerges at 50 to 60 peak eosinophils/hpf and established EoD as part of a spectrum of upper eosinophilic gastrointestinal disorder associated with type 2 immunity and distinct from eosinophilic colitis and celiac disease. These findings provide a basis for improving diagnosis and treatment.

Statistical Power Analysis for Designing Bulk, Single-Cell, and Spatial Transcriptomics Experiments: Review, Tutorial, and Perspectives

Gene expression profiling technologies have been used in various applications such as cancer biology. The development of gene expression profiling has expanded the scope of target discovery in transcriptomic studies, and each technology produces data with distinct characteristics. In order to guarantee biologically meaningful findings using transcriptomic experiments, it is important to consider various experimental factors in a systematic way through statistical power analysis. In this paper, we review and discuss the power analysis for three types of gene expression profiling technologies from a practical standpoint, including bulk RNA-seq, single-cell RNA-seq, and high-throughput spatial transcriptomics. Specifically, we describe the existing power analysis tools for each research objective for each of the bulk RNA-seq and scRNA-seq experiments, along with recommendations. On the other hand, since there are no power analysis tools for high-throughput spatial transcriptomics at this point, we instead investigate the factors that can influence power analysis.

Predicting severity in COVID-19 disease using sepsis blood gene expression signatures

Severely-afflicted COVID-19 patients can exhibit disease manifestations representative of sepsis, including acute respiratory distress syndrome and multiple organ failure. We hypothesized that diagnostic tools used in managing all-cause sepsis, such as clinical criteria, biomarkers, and gene expression signatures, should extend to COVID-19 patients. Here we analyzed the whole blood transcriptome of 124 early (1-5 days post-hospital admission) and late (6-20 days post-admission) sampled patients with confirmed COVID-19 infections from hospitals in Quebec, Canada. Mechanisms associated with COVID-19 severity were identified between severity groups (ranging from mild disease to the requirement for mechanical ventilation and mortality), and established sepsis signatures were assessed for dysregulation. Specifically, gene expression signatures representing pathophysiological events, namely cellular reprogramming, organ dysfunction, and mortality, were significantly enriched and predictive of severity and lethality in COVID-19 patients. Mechanistic endotypes reflective of distinct sepsis aetiologies and therapeutic opportunities were also identified in subsets of patients, enabling prediction of potentially-effective repurposed drugs. The expression of sepsis gene expression signatures in severely-afflicted COVID-19 patients indicates that these patients should be classified as having severe sepsis. Accordingly, in severe COVID-19 patients, these signatures should be strongly considered for the mechanistic characterization, diagnosis, and guidance of treatment using repurposed drugs.

Stem transcriptome screen for selection in wild and cultivated pitahaya (Selenicereus undatus): an epiphytic cactus with edible fruit

Dragon fruit, pitahaya or pitaya are common names for the species in the Hylocereus group of Selenicereus that produce edible fruit. These Neotropical epiphytic cacti are considered promising underutilized crops and are currently cultivated around the world. The most important species, S. undatus, has been managed in the Maya domain for centuries and is the focus of this article. Transcriptome profiles from stems of wild and cultivated plants of this species were compared. We hypothesized that differences in transcriptomic signatures could be associated with genes related to drought stress. De novo transcriptome assembly and the analysis of differentially expressed genes (DEGs) allowed us to identify a total of 9,203 DEGs in the Hunucmá cultivar relative of wild Mozomboa plants. Of these, 4,883 represent up-regulated genes and 4,320, down-regulated genes. Additionally, 6,568 DEGs were identified from a comparison between the Umán cultivar and wild plants, revealing 3,286 up-regulated and 3,282 down-regulated genes. Approximately half of the DEGs are shared by the two cultivated plants. Differences between the two cultivars that were collected in the same region could be the result of differences in management. Metabolism was the most representative functional category in both cultivars. The up-regulated genes of both cultivars formed a network related to the hormone-mediated signaling pathway that includes cellular responses to auxin stimulus and to hormone stimulus. These cellular reactions have been documented in several cultivated plants in which drought-tolerant cultivars modify auxin transport and ethylene signaling, resulting in a better redistribution of assimilates.

Maize Breeding for Low Nitrogen Inputs in Agriculture: Mechanisms Underlying the Tolerance to the Abiotic Stress

Nitrogen (N) is essential for sustaining life on Earth and plays a vital role in plant growth and thus agricultural production. The excessive use of N fertilizers not only harms the economy, but also the environment. In the context of the environmental impacts caused by agriculture, global maize improvement programs aim to develop cultivars with high N-use efficiency (NUE) to reduce the use of N fertilizers. Since N is highly mobile in plants, NUE is related to numerous little-known morphophysiological and molecular mechanisms. In this review paper we present an overview of the morpho-physiological adaptations of shoot and root, molecular mechanisms involved in plant response to low nitrogen environment, and the genetic effects involved in the control of key traits for NUE. Some studies show that the efficiency of cultivars growing under low N is related to deep root architecture, more lateral roots (LR), and sparser branching of LR, resulting in lower metabolic costs. The NUE cultivars also exhibit more efficient photosynthesis, which affects plant growth under suboptimal nitrogen conditions. In this sense, obtaining superior genotypes for NUE can be achieved with the exploitation of heterosis, as non-additive effects are more important in the expression of traits associated with NUE.

Proteome Remodeling of the Eye Lens at 50 Years Identified With Data-Independent Acquisition

The eye lens is responsible for focusing and transmitting light to the retina. The lens does this in the absence of organelles, yet maintains transparency for at least 5 decades before onset of age-related nuclear cataract (ARNC). It is hypothesized that oxidative stress contributes significantly to ARNC formation. It is in addition hypothesized that transparency is maintained by a microcirculation system that delivers antioxidants to the lens nucleus and exports small molecule waste. Common data-dependent acquisition methods are hindered by dynamic range of lens protein expression and provide limited context to age-related changes in the lens. In this study, we utilized data-independent acquisition mass spectrometry to analyze the urea-insoluble membrane protein fractions of 16 human lenses subdivided into three spatially distinct lens regions to characterize age-related changes, particularly concerning the lens microcirculation system and oxidative stress response. In this pilot cohort, we measured 4788 distinct protein groups, 46,681 peptides, and 7592 deamidated sequences, more than in any previous human lens data-dependent acquisition approach. Principally, we demonstrate that a significant proteome remodeling event occurs at approximately 50 years of age, resulting in metabolic preference for anaerobic glycolysis established with organelle degradation, decreased abundance of protein networks involved in calcium-dependent cell-cell contacts while retaining networks related to oxidative stress response. Furthermore, we identified multiple antioxidant transporter proteins not previously detected in the human lens and describe their spatiotemporal and age-related abundance changes. Finally, we demonstrate that aquaporin-5, among other proteins, is modified with age by post-translational modifications including deamidation and truncation. We suggest that the continued accumulation of each of these age-related outcomes in proteome remodeling contribute to decreased fiber cell permeability and result in ARNC formation.

Integrated Molecular Analysis Reveals 2 Distinct Subtypes of Pure Seminoma of the Testis

Objective

Testicular germ cell tumors (TGCT) are the most common solid malignancy in adolescent and young men, with a rising incidence over the past 20 years. Overall, TGCTs are second in terms of the average life years lost per person dying of cancer, and clinical therapeutics without adverse long-term side effects are lacking. Platinum-based regimens for TGCTs have heterogeneous outcomes even within the same histotype that frequently leads to under- and over-treatment. Understanding of molecular differences that lead to diverse outcomes of TGCT patients may improve current treatment approaches. Seminoma is the most common subtype of TGCTs, which can either be pure or present in combination with other histotypes.

Methods

Here we conducted a computational study of 64 pure seminoma samples from The Cancer Genome Atlas, applied consensus clustering approach to their transcriptomic data and revealed 2 clinically relevant seminoma subtypes: seminoma subtype 1 and 2.

Results

Our analysis identified significant differences in pluripotency stage, activity of double stranded DNA breaks repair mechanisms, rates of loss of heterozygosity, and expression of lncRNA responsible for cisplatin resistance between the subtypes. Seminoma subtype 1 is characterized by higher pluripotency state, while subtype 2 showed attributes of reprograming into non-seminomatous TGCT. The seminoma subtypes we identified may provide a molecular underpinning for variable responses to chemotherapy and radiation.

Conclusion

Translating our findings into clinical care may help improve risk stratification of seminoma, decrease overtreatment rates, and increase long-term quality of life for TGCT survivors.

TOFA-PREDICT study protocol: a stratification trial to determine key immunological factors predicting tofacitinib efficacy and drug-free remission in psoriatic arthritis (PsA)

INTRODUCTION

Psoriatic arthritis (PsA) is a chronic, inflammatory, musculoskeletal disease that affects up to 30% of patients with psoriasis. Current challenges in clinical care and research include personalised treatment, understanding the divergence of therapy response and unravelling the multifactorial pathophysiology of this complex disease. Moreover, there is an urgent clinical need to predict, assess and understand the cellular and molecular pathways underlying the response to disease-modifying antirheumatic drugs (DMARDs). The TOFA-PREDICT clinical trial addresses this need. Our primary objective is to determine key immunological factors predicting tofacitinib efficacy and drug-free remission in PsA.

METHODS AND ANALYSIS

In this investigator-initiated, phase III, multicentre, open-label, four-arm randomised controlled trial, we plan to integrate clinical, molecular and imaging parameters of 160 patients with PsA. DMARD-naïve patients are randomised to methotrexate or tofacitinib. Additionally, patients who are non-responsive to conventional synthetic (cs)DMARDs continue their current csDMARD and are randomised to etanercept or tofacitinib. This results in four arms each with 40 patients. Patients are followed for 1 year. Treatment response is defined as minimal disease activity at week 16. Clinical data, biosamples and images are collected at baseline, 4 weeks and 16 weeks; at treatment failure (treatment switch) and 52 weeks. For the first 80 patients, we will use a systems medicine approach to assess multiomics biomarkers and develop a prediction model for treatment response. Subsequently, data from the second 80 patients will be used for validation.

ETHICS AND DISSEMINATION

The study was approved by the Medical Research Ethics Committee in Utrecht, Netherlands, is registered in the European Clinical Trials Database and is carried out in accordance with the Declaration of Helsinki. The study's progress is monitored by Julius Clinical, a science-driven contract research organisation.

TRIAL REGISTRATION NUMBER

EudraCT: 2017-003900-28.

Inflammation and Interferon Signatures in Peripheral B-Lymphocytes and Sera of Individuals With Fibromyalgia

Fibromyalgia (FM) is an idiopathic chronic disease characterized by widespread musculoskeletal pain, hyperalgesia and allodynia, often accompanied by fatigue, cognitive dysfunction and other symptoms. Autoimmunity and neuroinflammatory mechanisms have been suggested to play important roles in the pathophysiology of FM supported by recently identified interferon signatures in affected individuals. However, the contribution of different components in the immune system, such as the B-lymphocytes, in the progression to FM are yet unknown. Furthermore, there is a great need for biomarkers that may improve diagnostics of FM. Herein, we investigated the gene expression profile in peripheral B-cells, as well as a panel of inflammatory serum proteins, in 30 FM patients and 23 healthy matched control individuals. RNA sequence analysis revealed 60 differentially expressed genes when comparing the two groups. The group of FM patients showed increased expression of twenty-five interferon-regulated genes, such as S100A8 and S100A9, VCAM, CD163, SERPINA1, ANXA1, and an increased interferon score. Furthermore, FM was associated with elevated levels of 19 inflammatory serum proteins, such as IL8, AXIN1, SIRT2 and STAMBP, that correlated with the FM severity score. Together, the results shows that FM is associated with an interferon signature in B-cells and increased levels of a set of inflammatory serum proteins. Our findings bring further support for immune activation in the pathogenesis of FM and highlight candidate biomarkers for diagnosis and intervention in the management of FM.

De novo genome assembly and in natura epigenomics reveal salinity-induced DNA methylation in the mangrove tree Bruguiera gymnorhiza

Mangroves are adapted to harsh environments, such as high ultraviolet (UV) light, low nutrition, and fluctuating salinity in coastal zones. However, little is known about the transcriptomic and epigenomic basis of the resilience of mangroves due to limited available genome resources. We performed a de novo genome assembly and in natura epigenome analyses of the mangrove Bruguiera gymnorhiza, one of the dominant mangrove species. We also performed the first genome-guided transcriptome assembly for mangrove species. The 309 Mb of the genome is predicted to encode 34 403 genes and has a repeat content of 48%. Depending on its growing environment, the natural B. gymnorhiza population showed drastic morphological changes associated with expression changes in thousands of genes. Moreover, high-salinity environments induced genome-wide DNA hypermethylation of transposable elements (TEs) in the B. gymnorhiza. DNA hypermethylation was concurrent with the transcriptional regulation of chromatin modifier genes, suggesting robust epigenome regulation of TEs in the B. gymnorhiza genome under high-salinity environments. The genome and epigenome data in this study provide novel insights into the epigenome regulation of mangroves and a better understanding of the adaptation of plants to fluctuating, harsh natural environments.

BNT162b2 vaccination enhances interferon-JAK-STAT-regulated antiviral programs in COVID-19 patients infected with the SARS-CoV-2 Beta variant

Background

SARS-CoV-2 infection activates interferon-controlled signaling pathways and elicits a wide spectrum of immune responses and clinical manifestations in human patients.

Methods

Here, we investigate the impact of prior vaccination on the innate immune response of hospitalized COVID-19 patients infected with the SARS-CoV-2 Beta variant through RNA sequencing of peripheral blood immune cells. Four patients had received the first dose of BNT162b2 about 11 days prior to the onset of COVID-19 symptoms and five patients were unvaccinated. Patients had received dexamethasone treatment. Immune transcriptomes were obtained at days 7-13, 20-32 and 42-60 after first symptomology.

Results

RNA-seq reveals an enhanced JAK-STAT-mediated immune transcriptome response at day 10 in vaccinated patients as compared to unvaccinated ones. This increase subsides by day 35. Expression of the gene encoding the antiviral protein oligoadenylate synthetase (OAS) 1, which is inversely correlated with disease severity, and other key antiviral proteins increases in the vaccinated group. We also investigate the immune transcriptome in naïve individuals receiving their first dose of BNT162b2 and identify a gene signature shared with the vaccinated COVID-19 patients.

Conclusions

Our study demonstrates that RNA-seq can be used to monitor molecular immune responses elicited by the BNT162b2 vaccine, both in naïve individuals and in COVID-19 patients, and it provides a biomarker-based approach to systems vaccinology.

Analysis of Salivary Microbiome and Its Association With Periodontitis in Patients With Obstructive Sleep Apnea

Objectives

This study aimed to analyze the periodontal conditions of patients with obstructive sleep apnea (OSA) in relation to the salivary microbiome.

Materials and Methods

In total, 54 male adults (27 with OSA, 27 controls) completed this cross-sectional study. All participants were monitored by overnight polysomnography (PSG) and underwent full-mouth periodontal examination. Saliva samples were then collected, and the microbial 16S ribosomal RNA gene was sequenced. The data were analyzed to determine the microbial distribution and the community structure of the two groups.

Results

Demonstrated by alpha and beta diversity, the OSA group had a lower microbial richness and a lower observed species than the controls. There was no significant difference in the microbial species diversity or evenness between the OSA and the non-OSA groups. The OSA group had fewer operational taxonomic units (OTUs), and the distribution of microbiome showed that several gram-positive bacteria had higher abundance in the OSA group. As for periodontal pathogens, the relative abundance of Prevotella was significantly increased in the OSA group. No significant difference was observed in the relative abundance of other pathogens at either the genus or species level.

Conclusions

The salivary microbial community structure was altered in patients with OSA in terms of species richness and trans-habitat diversity, along with an increase in Prevotella, a specific periodontal pathogen. These findings might explain the high prevalence of periodontitis in OSA patients.

Airway microbial communities, smoking and asthma in a general population sample

BACKGROUND

Normal airway microbial communities play a central role in respiratory health but are poorly characterized. Cigarette smoking is the dominant global environmental influence on lung function, and asthma has become the most prevalent chronic respiratory disease worldwide. Both conditions have major microbial components that are incompletely defined.

METHODS

We investigated airway bacterial communities in a general population sample of 529 Australian adults. Posterior oropharyngeal swabs were analyzed by sequencing of the 16S rRNA gene. The microbiota were characterized according to their prevalence, abundance and network memberships.

FINDINGS

The microbiota were similar across the general population, and were strongly organized into co-abundance networks. Smoking was associated with diversity loss, negative effects on abundant taxa, profound alterations to network structure and expansion of Streptococcus spp. By contrast, the asthmatic microbiota were selectively affected by an increase in Neisseria spp. and by reduced numbers of low abundance but prevalent organisms.

INTERPRETATION

Our study shows that the healthy airway microbiota in this population were contained within a highly structured ecosystem, suggesting balanced relationships between the microbiome and human host factors. The marked abnormalities in smokers may contribute to chronic obstructive pulmonary disease (COPD) and lung cancer. The narrow spectrum of abnormalities in asthmatics encourages investigation of damaging and protective effects of specific bacteria.

FUNDING

The study was funded by the Asmarley Trust and a Wellcome Joint Senior Investigator Award to WOCC and MFM (WT096964MA and WT097117MA). The Busselton Healthy Ageing Study is supported by the Government of Western Australia (Office of Science, Department of Health) the City of Busselton, and private donations.

A biallelic SNIP1 Amish founder variant causes a recognizable neurodevelopmental disorder

SNIP1 (Smad nuclear interacting protein 1) is a widely expressed transcriptional suppressor of the TGF-β signal-transduction pathway which plays a key role in human spliceosome function. Here, we describe extensive genetic studies and clinical findings of a complex inherited neurodevelopmental disorder in 35 individuals associated with a SNIP1 NM_024700.4:c.1097A>G, p.(Glu366Gly) variant, present at high frequency in the Amish community. The cardinal clinical features of the condition include hypotonia, global developmental delay, intellectual disability, seizures, and a characteristic craniofacial appearance. Our gene transcript studies in affected individuals define altered gene expression profiles of a number of molecules with well-defined neurodevelopmental and neuropathological roles, potentially explaining clinical outcomes. Together these data confirm this SNIP1 gene variant as a cause of an autosomal recessive complex neurodevelopmental disorder and provide important insight into the molecular roles of SNIP1, which likely explain the cardinal clinical outcomes in affected individuals, defining potential therapeutic avenues for future research.

Characterization and comparative analysis of transcriptional profiles of porcine colostrum and mature milk at different parities

BACKGROUND

Porcine milk is a complex fluid, containing a myriad of immunological, biochemical, and cellular components, made to satisfy the nutritional requirements of the neonate. Whole milk contains many different cell types, including mammary epithelial cells, neutrophils, macrophages, and lymphocytes, as well nanoparticles, such as milk exosomes. To-date, only a limited number of livestock transcriptomic studies have reported sequencing of milk. Moreover, those studies focused only on sequencing somatic cells as a proxy for the mammary gland with the goal of investigating differences in the lactation process. Recent studies have indicated that RNA originating from multiple cell types present in milk can withstand harsh environments, such as the digestive system, and transmit regulatory molecules from maternal to neonate. Transcriptomic profiling of porcine whole milk, which is reflective of the combined cell populations, could help elucidate these mechanisms. To this end, total RNA from colostrum and mature milk samples were sequenced from 65 sows at differing parities. A stringent bioinformatic pipeline was used to identify and characterize 70,841 transcripts.

RESULTS

The 70,841 identified transcripts included 42,733 previously annotated transcripts and 28,108 novel transcripts. Differential gene expression analysis was conducted using a generalized linear model coupled with the Lancaster method for P-value aggregation across transcripts. In total, 1667 differentially expressed genes (DEG) were identified for the milk type main effect, and 33 DEG were identified for the milk type x parity interaction. Several gene ontology (GO) terms related to immune response were significant for the milk type main effect, supporting the well-known fact that immunoglobulins and immune cells are transferred to the neonate via colostrum.

CONCLUSIONS

This is the first study to perform global transcriptome analysis from whole milk samples in sows from different parities. Our results provide important information and insight into synthesis of milk proteins and innate immunity and potential targets for future improvement of swine lactation and piglet development.

The effects of aging on molecular modulators of human embryo implantation

Advancing age has a negative impact on female fertility. As implantation rates decline during the normal maternal life course, age-related, embryonic factors are altered and our inability to monitor these factors in an unbiased genome-wide manner in vivo has severely limited our understanding of early human embryo development and implantation. Our high-throughput methodology uses trophectoderm samples representing the full spectrum of maternal reproductive ages with embryo implantation potential examined in relation to trophectoderm transcriptome dynamics and reproductive maternal age. Potential embryo-endometrial interactions were tested using trophectoderm sampled from young women, with the receptive uterine environment representing the most 'fertile' environment for successful embryo implantation. Potential roles for extracellular exosomes, embryonic metabolism and regulation of apoptosis were revealed. These biomarkers are consistent with embryo-endometrial crosstalk/developmental competency, serving as a mediator for successful implantation. Our data opens the door to developing a diagnostic test for predicting implantation success in women undergoing fertility treatment.

Comparative Transcriptome Analysis of Milk Somatic Cells During Lactation Between Two Intensively Reared Dairy Sheep Breeds

In dairy sheep industry, milk production dictates the value of a ewe. Milk production is directly related to the morphology and physiology of the mammary gland; both being designated targets of breeding strategies. Although within a flock breeding parameters are mutual, large differences in milk production among individual ewes are usually observed. In this work, we tested two of the most productive dairy sheep breeds reared intensively in Greece, one local the Chios breed and one foreign the Lacaune breed. We used transcriptome sequencing to reveal molecular mechanisms that render the mammary gland highly productive or not. While highly expressed genes (caseins and major whey protein genes) were common among breeds, differences were observed in differentially expressed genes. ENSOARG00000008077, as a member of ribosomal protein 14 family, together with LPCAT2, CCR3, GPSM2, ZNF131, and ASIP were among the genes significantly differentiating mammary gland's productivity in high yielding ewes. Gene ontology terms were mainly linked to the inherent transcriptional activity of the mammary gland (GO:0005524, GO:0030552, GO:0016740, GO:0004842), lipid transfer activity (GO:0005319) and innate immunity (GO:0002376, GO:0075528, GO:0002520). In addition, clusters of genes affecting zinc and iron trafficking into mitochondria were highlighted for high yielding ewes (GO:0071294, GO:0010043). Our analyses provide insights into the molecular pathways involved in lactation between ewes of different performances. Results revealed management issues that should be addressed by breeders in order to move toward increased milk yields through selection of the desired phenotypes. Our results will also contribute toward the selection of the most resilient and productive ewes, thus, will strengthen the existing breeding systems against a spectrum of environmental threats.

Hierarchicell: an R-package for estimating power for tests of differential expression with single-cell data

Background

Study design is a critical aspect of any experiment, and sample size calculations for statistical power that are consistent with that study design are central to robust and reproducible results. However, the existing power calculators for tests of differential expression in single-cell RNA-seq data focus on the total number of cells and not the number of independent experimental units, the true unit of interest for power. Thus, current methods grossly overestimate the power.

Results

Hierarchicell is the first single-cell power calculator to explicitly simulate and account for the hierarchical correlation structure (i.e., within sample correlation) that exists in single-cell RNA-seq data. Hierarchicell, an R-package available on GitHub, estimates the within sample correlation structure from real data to simulate hierarchical single-cell RNA-seq data and estimate power for tests of differential expression. This multi-stage approach models gene dropout rates, intra-individual dispersion, inter-individual variation, variable or fixed number of cells per individual, and the correlation among cells within an individual. Without modeling the within sample correlation structure and without properly accounting for the correlation in downstream analysis, we demonstrate that estimates of power are falsely inflated. Hierarchicell can be used to estimate power for binary and continuous phenotypes based on user-specified number of independent experimental units (e.g., individuals) and cells within the experimental unit.

Conclusions

Hierarchicell is a user-friendly R-package that provides accurate estimates of power for testing hypotheses of differential expression in single-cell RNA-seq data. This R-package represents an important addition to single-cell RNA analytic tools and will help researchers design experiments with appropriate and accurate power, increasing discovery and improving robustness and reproducibility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07635-w.

Transcriptome analyses of 7-day-old zebrafish larvae possessing a familial Alzheimer's disease-like mutation in psen1 indicate effects on oxidative phosphorylation, ECM and MCM functions, and iron homeostasis

BACKGROUND

Early-onset familial Alzheimer's disease (EOfAD) is promoted by dominant mutations, enabling the study of Alzheimer's disease (AD) pathogenic mechanisms through generation of EOfAD-like mutations in animal models. In a previous study, we generated an EOfAD-like mutation, psen1Q96_K97del, in zebrafish and performed transcriptome analysis comparing entire brains from 6-month-old wild type and heterozygous mutant fish. We identified predicted effects on mitochondrial function and endolysosomal acidification. Here we aimed to determine whether similar effects occur in 7 day post fertilization (dpf) zebrafish larvae that might be exploited in screening of chemical libraries to find ameliorative drugs.

RESULTS

We generated clutches of wild type and heterozygous psen1Q96_K97del 7 dpf larvae using a paired-mating strategy to reduce extraneous genetic variation before performing a comparative transcriptome analysis. We identified 228 differentially expressed genes and performed various bioinformatics analyses to predict cellular functions.

CONCLUSIONS

Our analyses predicted a significant effect on oxidative phosphorylation, consistent with our earlier observations of predicted effects on ATP synthesis in adult heterozygous psen1Q96_K97del brains. The dysregulation of minichromosome maintenance protein complex (MCM) genes strongly contributed to predicted effects on DNA replication and the cell cycle and may explain earlier observations of genome instability due to PSEN1 mutation. The upregulation of crystallin gene expression may be a response to defective activity of mutant Psen1 protein in endolysosomal acidification. Genes related to extracellular matrix (ECM) were downregulated, consistent with previous studies of EOfAD mutant iPSC neurons and postmortem late onset AD brains. Also, changes in expression of genes controlling iron ion transport were observed without identifiable changes in the prevalence of transcripts containing iron responsive elements (IREs) in their 3' untranslated regions (UTRs). These changes may, therefore, predispose to the apparent iron dyshomeostasis previously observed in 6-month-old heterozygous psen1Q96_K97del EOfAD-like mutant brains.

Impact of Sequencing Depth and Library Preparation on Toxicological Interpretation of RNA-Seq Data in a "Three-Sample" Scenario

While RNA-sequencing (RNA-seq) has emerged as a standard approach in toxicogenomics, its full potential in gaining underlying toxicological mechanisms is still not clear when only three biological replicates are used. This "three-sample" study design is common in toxicological research, particularly in animal studies during preclinical drug development. Sequencing depth (the total number of reads in an experiment) and library preparation are critical to the resolution and integrity of RNA-seq data and biological interpretation. We used aflatoxin b1 (AFB1), a model toxicant, to investigate the effect of sequencing depth and library preparation in RNA-seq on toxicological interpretation in the "three-sample" scenario. We also compared different gene profiling platforms (RNA-seq, TempO-seq, microarray, and qPCR) using identical liver samples. Well-established mechanisms of AFB1 toxicity served as ground truth for our comparative analyses. We found that a minimum of 20 million reads was sufficient to elicit key toxicity functions and pathways underlying AFB1-induced liver toxicity using three replicates and that identification of differentially expressed genes was positively associated with sequencing depth to a certain extent. Further, our results showed that RNA-seq revealed toxicological insights from pathway enrichment with overall higher statistical power and overlap ratio, compared with TempO-seq and microarray. Moreover, library preparation using the same methods was important to reproducing the toxicological interpretation.

Evidence for the placenta-brain axis: multi-omic kernel aggregation predicts intellectual and social impairment in children born extremely preterm

BACKGROUND

Children born extremely preterm are at heightened risk for intellectual and social impairment, including Autism Spectrum Disorder (ASD). There is increasing evidence for a key role of the placenta in prenatal developmental programming, suggesting that the placenta may, in part, contribute to origins of neurodevelopmental outcomes.

METHODS

We examined associations between placental transcriptomic and epigenomic profiles and assessed their ability to predict intellectual and social impairment at age 10 years in 379 children from the Extremely Low Gestational Age Newborn (ELGAN) cohort. Assessment of intellectual ability (IQ) and social function was completed with the Differential Ability Scales-II and Social Responsiveness Scale (SRS), respectively. Examining IQ and SRS allows for studying ASD risk beyond the diagnostic criteria, as IQ and SRS are continuous measures strongly correlated with ASD. Genome-wide mRNA, CpG methylation and miRNA were assayeds with the Illumina Hiseq 2500, HTG EdgeSeq miRNA Whole Transcriptome Assay, and Illumina EPIC/850 K array, respectively. We conducted genome-wide differential analyses of placental mRNA, miRNA, and CpG methylation data. These molecular features were then integrated for a predictive analysis of IQ and SRS outcomes using kernel aggregation regression. We lastly examined associations between ASD and the multi-omic-predicted component of IQ and SRS.

RESULTS

Genes with important roles in neurodevelopment and placental tissue organization were associated with intellectual and social impairment. Kernel aggregations of placental multi-omics strongly predicted intellectual and social function, explaining approximately 8% and 12% of variance in SRS and IQ scores via cross-validation, respectively. Predicted in-sample SRS and IQ showed significant positive and negative associations with ASD case-control status.

LIMITATIONS

The ELGAN cohort comprises children born pre-term, and generalization may be affected by unmeasured confounders associated with low gestational age. We conducted external validation of predictive models, though the sample size (N = 49) and the scope of the available out-sample placental dataset are limited. Further validation of the models is merited.

CONCLUSIONS

Aggregating information from biomarkers within and among molecular data types improves prediction of complex traits like social and intellectual ability in children born extremely preterm, suggesting that traits within the placenta-brain axis may be omnigenic.

Dedicated transcriptomics combined with power analysis lead to functional understanding of genes with weak phenotypic changes in knockout lines

Systematic knockout studies in mice have shown that a large fraction of the gene replacements show no lethal or other overt phenotypes. This has led to the development of more refined analysis schemes, including physiological, behavioral, developmental and cytological tests. However, transcriptomic analyses have not yet been systematically evaluated for non-lethal knockouts. We conducted a power analysis to determine the experimental conditions under which even small changes in transcript levels can be reliably traced. We have applied this to two gene disruption lines of genes for which no function was known so far. Dedicated phenotyping tests informed by the tissues and stages of highest expression of the two genes show small effects on the tested phenotypes. For the transcriptome analysis of these stages and tissues, we used a prior power analysis to determine the number of biological replicates and the sequencing depth. We find that under these conditions, the knockouts have a significant impact on the transcriptional networks, with thousands of genes showing small transcriptional changes. GO analysis suggests that A930004D18Rik is involved in developmental processes through contributing to protein complexes, and A830005F24Rik in extracellular matrix functions. Subsampling analysis of the data reveals that the increase in the number of biological replicates was more important that increasing the sequencing depth to arrive at these results. Hence, our proof-of-principle experiment suggests that transcriptomic analysis is indeed an option to study gene functions of genes with weak or no traceable phenotypic effects and it provides the boundary conditions under which this is possible.

Knockout mice benefit the understanding of gene functions in mammals. However, it has proven difficult for many genes to identify clear phenotypes, related due to lack of sufficient assays. As Lewis Wolpert put it in a famous quote “But did you take them to the opera?”, thus metaphorically alluding to the need to extend phenotyping efforts. This insight led to the establishment of phenotyping pipelines that are nowadays routinely used to characterize knock-out lines. However, transcriptomic approaches based on RNA-Seq have been much less explored for such deep-level studies. We conducted here both, a theoretical power analysis and practical RNA-Seq experiments on two knockout lines with small phenotypic effects to investigate the parameters including sample size, sequencing depth, fold change, and dispersion. Our dedicated RNA-Seq studies discovered thousands of genes with small transcriptional changes and enriched in specific functions in both knockout lines. We find that it is more important to increase the number of samples than to increase the sequencing depth. Our work shows that a deep RNA-Seq study on knockouts is powerful for understanding gene functions in cases of weak phenotypic effects, and provides a guideline for the experimental design of such studies.

Near normalization of peripheral blood markers in HIV-infected patients on long-term suppressive antiretroviral therapy: a case-control study

OBJECTIVE

To explore the differences in peripheral blood markers between HIV well controlled patients on long-term suppressive antiretroviral therapy (HIV-group) and age-matched healthy controls, to evaluate the benefits of virological suppression in those patients.

METHODS

We performed a case-control study in 22 individuals in the HIV-group and 14 in the healthy control-group. RNA-seq analysis was performed from peripheral blood mononuclear cells. Peripheral blood T-cell subsets were evaluated by flow cytometry and plasma biomarkers by immunoassays. All P values were corrected by the false discovery rate (q values).

RESULTS

Only the serine/arginine repetitive matrix 4 gene, which is involved in alternative RNA splicing events, was differentially expressed between HIV and healthy control groups (q value ≤0.05 and fold-change ≥2). However, 147 differentially expressed genes were found with a more relaxed threshold (P value ≤0.05 and fold-change ≥1.5), of which 67 genes with values of variable importance in projection at least one were selected for pathway analysis. We found that six ribosomal genes represented significant ribosome-related pathways, all of them downregulated in the HIV-group, which may be a strategy to facilitate viral production. T cells subset and plasma biomarkers did not show significant differences after false discovery rate correction (q value >0.05), but a noncorrected analysis showed higher values of regulatory CD4 T cells (CD4CD25CD127), MCP-1, and sVEGF-R1 in the HIV-group (P value ≤0.05).

CONCLUSION

T-cell subsets, plasma biomarkers, and gene expression were close to normalization in HIV-infected patients on long-term suppressive combination antiretroviral therapy compared with healthy controls. However, residual alterations remain, mainly at the gene expression, which still reveals the impact of HIV infection in these patients.

Wnt-regulated lncRNA discovery enhanced by in vivo identification and CRISPRi functional validation

BACKGROUND

Wnt signaling is an evolutionarily conserved developmental pathway that is frequently hyperactivated in cancer. While multiple protein-coding genes regulated by Wnt signaling are known, the functional lncRNAs regulated by Wnt signaling have not been systematically characterized.

METHODS

We comprehensively mapped Wnt-regulated lncRNAs from an orthotopic Wnt-addicted pancreatic cancer model and examined the response of lncRNAs to Wnt inhibition between in vivo and in vitro cancer models. We further annotated and characterized these Wnt-regulated lncRNAs using existing genomic classifications (using data from FANTOM5) in the context of Wnt signaling and inferred their role in cancer pathogenesis (using GWAS and expression data from the TCGA). To functionally validate Wnt-regulated lncRNAs, we performed CRISPRi screens to assess their role in cancer cell proliferation both in vivo and in vitro.

RESULTS

We identified 3633 lncRNAs, of which 1503 were regulated by Wnt signaling in an orthotopic Wnt-addicted pancreatic cancer model. These lncRNAs were much more sensitive to changes in Wnt signaling in xenografts than in cultured cells. Our analysis suggested that Wnt signaling inhibition could influence the co-expression relationship of Wnt-regulated lncRNAs and their eQTL-linked protein-coding genes. Wnt-regulated lncRNAs were also implicated in specific gene networks involved in distinct biological processes that contribute to the pathogenesis of cancers. Consistent with previous genome-wide lncRNA CRISPRi screens, around 1% (13/1503) of the Wnt-regulated lncRNAs were found to modify cancer cell growth in vitro. This included CCAT1 and LINC00263, previously reported to regulate cancer growth. Using an in vivo CRISPRi screen, we doubled the discovery rate, identifying twice as many Wnt-regulated lncRNAs (25/1503) that had a functional effect on cancer cell growth.

CONCLUSIONS

Our study demonstrates the value of studying lncRNA functions in vivo, provides a valuable resource of lncRNAs regulated by Wnt signaling, and establishes a framework for systematic discovery of functional lncRNAs.

Blood pressure changes PVAT function and transcriptome: use of the mid-thoracic aorta coarcted rat

Perivascular adipose tissue (PVAT) modifies the contractile function of the vessel it surrounds (outside-in signaling). Little work points to the vessel actively affecting its surrounding PVAT. We hypothesized that inside-out arterial signaling to PVAT would be evidenced by the response of PVAT to changes in tangential vascular wall stress. Rats coarcted in the mid-thoracic aorta created PVAT tissues that would exemplify pressure-dependent changes (above vs. below coarctation); a sham rat was used as a control. Radiotelemetry revealed a ∼20 mmHg systolic pressure gradient across the coarctation 4 wk after surgery. Four measures (histochemical, adipocyte progenitor proliferation and differentiation, isometric tone, and bulk mRNA sequencing) were used to compare PVAT above versus below the ligature in sham and coarcted rats. Neither aortic collagen deposition in PVAT nor arterial media/radius ratio above coarctation was increased versus below segments. However, differentiated adipocytes derived from PVAT above the coarctation accumulated substantially less triglycerides versus those below; their relative proliferation rate as adipogenic precursors was not different. Functionally, the ability of PVAT to assist stress relaxation of isolated aorta was reduced in rings above versus below the coarctation. Transcriptomic analyses revealed that the coarctation resulted in more differentially expressed genes (DEGs) between PVAT above versus below when compared with sham samples from the same locations. A majority of DEGs were in PVAT below the coarctation and were enriched in neuronal/synaptic terms. These findings provide initial evidence that signaling from the vascular wall, as stimulated by a pressure change, influences the function and transcriptional profile of its PVAT.NEW & NOTEWORTHY A mid-thoracic aorta coarcted rat was created to generate a stable pressure difference above versus below the coarctation ligature. This study determined that the PVAT around the thoracic aorta exposed to a higher pressure has a significantly reduced ability to assist stress relaxation versus that below the ligature and appears to retain the ability to be anticontractile. At the same time, the PVAT around the thoracic aorta exposed to higher pressure had a reduced adipogenic potential versus that below the ligature. Transcriptomics analyses indicated that PVAT below the coarctation showed the greatest number of DEGs with an increased profile of the synaptic neurotransmitter gene network.

The noncoding and coding transcriptional landscape of the peripheral immune response in patients with COVID-19

BACKGROUND

COVID-19 is currently a global pandemic, but the response of human immune system to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection remains unclear. Noncoding RNAs serve as immune regulators and thus may play a critical role in disease progression.

METHODS

We performed multi-transcriptome sequencing of both noncoding RNAs and mRNAs isolated from the red blood cell depleted whole blood of moderate and severe COVID-19 patients. The functions of noncoding RNAs were validated by analyses of the expression of downstream mRNAs. We further utilized the single-cell RNA-seq data of COVID-19 patients from Wilk et al. and Chua et al. to characterize noncoding RNA functions in different cell types.

RESULTS

We defined four types of microRNAs with different expression tendencies that could serve as biomarkers for COVID-19 progress. We also identified miR-146a-5p, miR-21-5p, miR-142-3p, and miR-15b-5p as potential contributors to the disease pathogenesis, possibly serving as biomarkers of severe COVID-19 and as candidate therapeutic targets. In addition, the transcriptome profiles consistently suggested hyperactivation of the immune response, loss of T-cell function, and immune dysregulation in severe patients.

CONCLUSIONS

Collectively, these findings provide a comprehensive view of the noncoding and coding transcriptional landscape of peripheral immune cells during COVID-19, furthering our understanding and offering novel insights into COVID-19 pathogenesis.

RNA-seq analysis of galaninergic neurons from ventrolateral preoptic nucleus identifies expression changes between sleep and wake

BACKGROUND

Previous studies show that galanin neurons in ventrolateral preoptic nucleus (VLPO-Gal) are essential for sleep regulation. Here, we explored the transcriptional regulation of the VLPO-Gal neurons in sleep by comparing their transcriptional responses between sleeping mice and those kept awake, sacrificed at the same diurnal time.

RESULTS

RNA-sequencing (RNA-seq) analysis was performed on eGFP(+) galanin neurons isolated using laser captured microdissection (LCM) from VLPO. Expression of Gal was assessed in our LCM eGFP(+) neurons via real time qPCR and showed marked enrichment when compared to LCM eGFP(-) cells and to bulk VLPO samples. Gene set enrichment analysis utilizing data from a recent single-cell RNA-seq study of the preoptic area demonstrated that our VLPO-Gal samples were highly enriched with galanin-expressing inhibitory neurons, but not galanin-expressing excitatory neurons. A total of 263 genes were differentially expressed between sleep and wake in VLPO-Gal neurons. When comparing differentially expressed genes in VLPO-Gal neurons to differentially expressed genes in a wake-active neuronal region (the medial prefrontal cortex), evidence indicates that both systemic and cell-specific mechanisms contribute to the transcriptional regulation in VLPO-Gal neurons. In both wake-active and sleep-active neurons, ER stress pathways are activated by wake and cold-inducible RNA-binding proteins are activated by sleep. In contrast, expression of DNA repair genes is increased in VLPO-Gal during wakefulness, but increased in wake-active cells during sleep.

CONCLUSION

Our study identified transcriptomic responses of the galanin neurons in the ventrolateral preoptic nucleus during sleep and sleep deprivation. Data indicate that VLPO contains mainly sleep-active inhibitory galaninergic neurons. The VLPO galanin neurons show responses to sleep and wake similar to wake-active regions, indicating these responses, such as ER stress and cold-inducible RNA-binding proteins, are systemic affecting all neuronal populations. Region-specific differences in sleep/wake responses were also identified, in particular DNA repair. Our study expands knowledge about the transcriptional response of a distinct group of neurons essential for sleep.

Reference genome and transcriptome informed by the sex chromosome complement of the sample increase ability to detect sex differences in gene expression from RNA-Seq data

BACKGROUND

Human X and Y chromosomes share an evolutionary origin and, as a consequence, sequence similarity. We investigated whether the sequence homology between the X and Y chromosomes affects the alignment of RNA-Seq reads and estimates of differential expression. We tested the effects of using reference genomes and reference transcriptomes informed by the sex chromosome complement of the sample's genome on the measurements of RNA-Seq abundance and sex differences in expression.

RESULTS

The default genome includes the entire human reference genome (GRCh38), including the entire sequence of the X and Y chromosomes. We created two sex chromosome complement informed reference genomes. One sex chromosome complement informed reference genome was used for samples that lacked a Y chromosome; for this reference genome version, we hard-masked the entire Y chromosome. For the other sex chromosome complement informed reference genome, to be used for samples with a Y chromosome, we hard-masked only the pseudoautosomal regions of the Y chromosome, because these regions are duplicated identically in the reference genome on the X chromosome. We analyzed the transcript abundance in the whole blood, brain cortex, breast, liver, and thyroid tissues from 20 genetic female (46, XX) and 20 genetic male (46, XY) samples. Each sample was aligned twice: once to the default reference genome and then independently aligned to a reference genome informed by the sex chromosome complement of the sample, repeated using two different read aligners, HISAT and STAR. We then quantified sex differences in gene expression using featureCounts to get the raw count estimates followed by Limma/Voom for normalization and differential expression. We additionally created sex chromosome complement informed transcriptome references for use in pseudo-alignment using Salmon. Transcript abundance was quantified twice for each sample: once to the default target transcripts and then independently to target transcripts informed by the sex chromosome complement of the sample.

CONCLUSIONS

We show that regardless of the choice of the read aligner, using an alignment protocol informed by the sex chromosome complement of the sample results in higher expression estimates on the pseudoautosomal regions of the X chromosome in both genetic male and genetic female samples, as well as an increased number of unique genes being called as differentially expressed between the sexes. We additionally show that using a pseudo-alignment approach informed on the sex chromosome complement of the sample eliminates Y-linked expression in female XX samples.

Circulating miR-206 and Wnt-signaling are associated with cardiovascular complications and a history of preeclampsia in women

Women with a history of preeclampsia (PE) have increased risk of cardiovascular disease (CVD) later in life. However, the molecular determinants underlying this risk remain unclear. We sought to understand how circulating miRNA levels are affected by prior PE, and related to biological pathways underpinning cardiovascular disease. RNA sequencing was used to profile plasma levels of 2578 miRNAs in a retrospective study of women with a history of PE or normotensive pregnancy, in two independent cohorts with either acute coronary syndrome (ACS) (n = 17–18/group) or no ACS (n = 20/group). Differential miRNA alterations were assessed in relation to a history of PE (within each cohort) or ACS (across cohorts), and compared with miRNAs previously reported to be altered during PE. A history of PE was associated with altered levels of 30 and 20 miRNAs in the ACS and non-ACS cohorts, respectively, whereas ACS exposure was associated with alterations in 259 miRNAs. MiR-206 was identified at the intersection of all comparisons relating to past/current PE and ACS exposure, and has previously been implicated in atherogenic activities related to hepatocytes, vascular smooth muscle cells and macrophages. Integration of all differentially altered miRNAs with their predicted and experimentally validated targets in silico revealed a number of highly targeted genes with potential atherogenic functions (including NFAT5, CCND2 and SMAD2), and one significantly enriched KEGG biological pathway (Wnt signaling) that was shared between all exposure groups. The present study provides novel insights into miRNAs, target genes and biological pathways that may underlie the long-term cardiovascular sequelae of PE.

Harmonization of quality metrics and power calculation in multi-omic studies

Multi-omic studies combine measurements at different molecular levels to build comprehensive models of cellular systems. The success of a multi-omic data analysis strategy depends largely on the adoption of adequate experimental designs, and on the quality of the measurements provided by the different omic platforms. However, the field lacks a comparative description of performance parameters across omic technologies and a formulation for experimental design in multi-omic data scenarios. Here, we propose a set of harmonized Figures of Merit (FoM) as quality descriptors applicable to different omic data types. Employing this information, we formulate the MultiPower method to estimate and assess the optimal sample size in a multi-omics experiment. MultiPower supports different experimental settings, data types and sample sizes, and includes graphical for experimental design decision-making. MultiPower is complemented with MultiML, an algorithm to estimate sample size for machine learning classification problems based on multi-omic data.

Skeletal muscle reprogramming by breast cancer regardless of treatment history or tumor molecular subtype

Increased susceptibility to fatigue is a negative predictor of survival commonly experienced by women with breast cancer (BC). Here, we sought to identify molecular changes induced in human skeletal muscle by BC regardless of treatment history or tumor molecular subtype using RNA-sequencing (RNA-seq) and proteomic analyses. Mitochondrial dysfunction was apparent across all molecular subtypes, with the greatest degree of transcriptomic changes occurring in women with HER2/neu-overexpressing tumors, though muscle from patients of all subtypes exhibited similar pathway-level dysregulation. Interestingly, we found no relationship between anticancer treatments and muscle gene expression, suggesting that fatigue is a product of BC per se rather than clinical history. In vitro and in vivo experimentation confirmed the ability of BC cells to alter mitochondrial function and ATP content in muscle. These data suggest that interventions supporting muscle in the presence of BC-induced mitochondrial dysfunction may alleviate fatigue and improve the lives of women with BC.

Seq-ing answers: Current data integration approaches to uncover mechanisms of transcriptional regulation

Advancements in the field of next generation sequencing lead to the generation of ever-more data, with the challenge often being how to combine and reconcile results from different OMICs studies such as genome, epigenome and transcriptome. Here we provide an overview of the standard processing pipelines for ChIP-seq and RNA-seq as well as common downstream analyses. We describe popular multi-omics data integration approaches used to identify target genes and co-factors, and we discuss how machine learning techniques may predict transcriptional regulators and gene expression.

Obesity does not alter endometrial gene expression in women with endometriosis

RESEARCH QUESTION

Does obesity affect endometrial gene expression in women with endometriosis, specifically women with stage I disease?

DESIGN

Differential gene expression analysis was conducted on endometrium from women with and without endometriosis (n = 169). Women were diagnosed after surgical visualization and staged according to the revised American Society for Reproductive Medicine (stage I-IV). Women were grouped by body mass index (BMI) (kg/m²) as underweight, normal, pre-obese or obese. After accounting for menstrual cycle stage, endometrial gene expression was analysed by BMI (continuous and grouped) in women with endometriosis, and in non-endometriosis controls.

RESULTS

No significant interaction effect was found between BMI and endometriosis status on endometrial gene expression. We have previously reported that obese women with endometriosis have a reduced incidence of stage I disease; however, stratifying our analysis into stage I endometriosis versus combined II, III and IV endometriosis failed to reveal any differentially expressed endometrial genes between normal, pre-obese and obese patients.

CONCLUSIONS

Despite obesity having deleterious effects on endometrial gene expression in other gynaecological pathologies, e.g. endometrial cancer and polycystic ovary syndrome, our results do not support an association between BMI and altered endometrial gene expression in women with or without endometriosis.

Characterising the mechanisms underlying genetic resistance to amoebic gill disease in Atlantic salmon using RNA sequencing

Background

Gill health is one of the main concerns for Atlantic salmon aquaculture, and Amoebic Gill Disease (AGD), attributable to infection by the amoeba Neoparamoeba perurans, is a frequent cause of morbidity. In the absence of preventive measures, increasing genetic resistance of salmon to AGD via selective breeding can reduce the incidence of the disease and mitigate gill damage. Understanding the mechanisms leading to AGD resistance and the underlying causative genomic features can aid in this effort, while also providing critical information for the development of other control strategies. AGD resistance is considered to be moderately heritable, and several putative QTL have been identified. The aim of the current study was to improve understanding of the mechanisms underlying AGD resistance, and to identify putative causative genomic factors underlying the QTL. To achieve this, RNA was extracted from the gill and head kidney of AGD resistant and susceptible animals following a challenge with N. perurans, and sequenced.

Results

Comparison between resistant and susceptible animals primarily highlighted differences mainly in the local immune response in the gill, involving red blood cell genes and genes related to immune function and cell adhesion. Differentially expressed immune genes pointed to a contrast in Th2 and Th17 responses, which is consistent with the increased heritability observed after successive challenges with the amoeba. Five QTL-region candidate genes showed differential expression, including a gene connected to interferon responses (GVINP1), a gene involved in systemic inflammation (MAP4K4), and a positive regulator of apoptosis (TRIM39). Analyses of allele-specific expression highlighted a gene in the QTL region on chromosome 17, cellular repressor of E1A-stimulated genes 1 (CREG1), showing allelic differential expression suggestive of a cis-acting regulatory variant.

Conclusions

In summary, this study provides new insights into the mechanisms of resistance to AGD in Atlantic salmon, and highlights candidate genes for further functional studies that can further elucidate the genomic mechanisms leading to resistance and contribute to enhancing salmon health via improved genomic selection.

Combined Omics Reveals That Disruption of the Selenocysteine Lyase Gene Affects Amino Acid Pathways in Mice

Selenium is a nonmetal trace element that is critical for several redox reactions and utilized to produce the amino acid selenocysteine (Sec), which can be incorporated into selenoproteins. Selenocysteine lyase (SCL) is an enzyme which decomposes Sec into selenide and alanine, releasing the selenide to be further utilized to synthesize new selenoproteins. Disruption of the selenocysteine lyase gene (Scly) in mice (Scly-/- or Scly KO) led to obesity with dyslipidemia, hyperinsulinemia, glucose intolerance and lipid accumulation in the hepatocytes. As the liver is a central regulator of glucose and lipid homeostasis, as well as selenium metabolism, we aimed to pinpoint hepatic molecular pathways affected by the Scly gene disruption. Using RNA sequencing and metabolomics, we identified differentially expressed genes and metabolites in the livers of Scly KO mice. Integrated omics revealed that biological pathways related to amino acid metabolism, particularly alanine and glycine metabolism, were affected in the liver by disruption of Scly in mice with selenium adequacy. We further confirmed that hepatic glycine levels are elevated in male, but not in female, Scly KO mice. In conclusion, our results reveal that Scly participates in the modulation of hepatic amino acid metabolic pathways.

RNA Sequencing Data: Hitchhiker's Guide to Expression Analysis

Gene expression is the fundamental level at which the results of various genetic and regulatory programs are observable. The measurement of transcriptome-wide gene expression has convincingly switched from microarrays to sequencing in a matter of years. RNA sequencing (RNA-seq) provides a quantitative and open system for profiling transcriptional outcomes on a large scale and therefore facilitates a large diversity of applications, including basic science studies, but also agricultural or clinical situations. In the past 10 years or so, much has been learned about the characteristics of the RNA-seq data sets, as well as the performance of the myriad of methods developed. In this review, we give an overview of the developments in RNA-seq data analysis, including experimental design, with an explicit focus on the quantification of gene expression and statistical approachesfor differential expression. We also highlight emerging data types, such as single-cell RNA-seq and gene expression profiling using long-read technologies.

MicroRNA Profile of HCV Spontaneous Clarified Individuals, Denotes Previous HCV Infection

Factors involved in the spontaneous cleareance of a hepatitis C (HCV) infection are related to both HCV and the interaction with the host immune system, but little is known about the consequences after a spontaneous resolution. The main HCV extrahepatic reservoir is the peripheral blood mononuclear cells (PBMCs), and their transcriptional profile provides us information of innate and adaptive immune responses against an HCV infection. MicroRNAs regulate the innate and adaptive immune responses, and they are actively involved in the HCV cycle. High Throughput sequencing was used to analyze the miRNA profiles from PBMCs of HCV chronic naïve patients (CHC), individuals that spontaneously clarified HCV (SC), and healthy controls (HC). We did not find any differentially expressed miRNAs between SC and CHC. However, both groups showed similar expression differences (21 miRNAs) with respect to HC. This miRNA signature correctly classifies HCV-exposed (CHC and SC) vs. HC, with the has-miR-21-3p showing the best performance. The potentially targeted molecular pathways by these 21 miRNAs mainly belong to fatty acids pathways, although hippo signaling, extracellular matrix (ECM) interaction, proteoglycans-related, and steroid biosynthesis pathways were also altered. These miRNAs target host genes involved in an HCV infection. Thus, an HCV infection promotes molecular alterations in PBMCs that can be detected after an HCV spontaneous resolution, and the 21-miRNA signature is able to identify HCV-exposed patients (either CHC or SC).

Identification of key genes and its chromosome regions linked to drought responses in leaves across different crops through meta-analysis of RNA-Seq data

BACKGROUND

Our study is the first to provide RNA-Seq data analysis related to transcriptomic responses towards drought across different crops. The aim was to identify and map which genes play a key role in drought response on leaves across different crops. Forty-two RNA-seq samples were analyzed from 9 published studies in 7 plant species (Arabidopsis thaliana, Solanum lycopersicum, Zea mays, Vitis vinifera, Malus X domestica, Solanum tuberosum, Triticum aestivum).

RESULTS

Twenty-seven (16 up-regulated and 11 down-regulated) drought-regulated genes were commonly present in at least 7 of 9 studies, while 351 (147 up-regulated and 204 down-regulated) were commonly drought-regulated in 6 of 9 studies. Across all kind of leaves, the drought repressed gene-ontologies were related to the cell wall and membrane re-structuring such as wax biosynthesis, cell wall organization, fatty acid biosynthesis. On the other hand, drought-up-regulated biological processes were related to responses to osmotic stress, abscisic acid, water deprivation, abscisic-activated signalling pathway, salt stress, hydrogen peroxide treatment. A common metabolic feature linked to drought response in leaves is the repression of terpenoid pathways. There was an induction of AL1 (alfin-like), UGKYAH (trihelix), WRKY20, homeobox genes and members of the SET domain family in 6 of 9 studies. Several genes involved in detoxifying and antioxidant reactions, signalling pathways and cell protection were commonly modulated by drought across the 7 species. The chromosome (Chr) mapping of these key abiotic stress genes highlighted that Chr 4 in Arabidopsis thaliana, Chr 1 in Zea mays, Chr 2 and Chr 5 in Triticum aestivum contained a higher presence of drought-related genes compared to the other remaining chromosomes. In seedling studies, it is worth notice the up-regulation of ERF4 and ESE3 (ethylene), HVA22 (abscisic acid), TIR1 (auxin) and some transcription factors (MYB3, MYB94, MYB1, WRKY53 and WRKY20). In mature leaves, ERF1 and Alfin-like 1 were induced by drought while other transcription factors (YABBY5, ARR2, TRFL2) and genes involved phospholipid biosynthesis were repressed.

CONCLUSIONS

The identified and mapped genes might be potential targets of molecular breeding activities to develop cultivars with enhanced drought resistance and tolerance across different crops.