Dream: powerful differential expression analysis for repeated measures designs

A microbiota-directed complementary food intervention in 12-18-month-old Bangladeshi children improves linear growth

BACKGROUND

Globally, stunting affects ∼150 million children under five, while wasting affects nearly 50 million. Current interventions have had limited effectiveness in ameliorating long-term sequelae of undernutrition including stunting, cognitive deficits and immune dysfunction. Disrupted development of the gut microbiota has been linked to the pathogenesis of undernutrition, providing potentially new treatment approaches.

METHODS

124 Bangladeshi children with moderate acute malnutrition (MAM) enrolled (at 12-18 months) in a previously reported 3-month RCT of a microbiota-directed complementary food (MDCF-2) were followed for two years. Weight and length were monitored by anthropometry, the abundances of bacterial strains were assessed by quantifying metagenome-assembled genomes (MAGs) in serially collected fecal samples and levels of growth-associated proteins were measured in plasma.

FINDINGS

Children who had received MDCF-2 were significantly less stunted during follow-up than those who received a standard ready-to-use supplementary food (RUSF) [linear mixed-effects model, βtreatment group x study week (95% CI) = 0.002 (0.001, 0.003); P = 0.004]. They also had elevated fecal abundances of Agathobacter faecis, Blautia massiliensis, Lachnospira and Dialister, plus increased levels of a group of 37 plasma proteins (linear model; FDR-adjusted P < 0.1), including IGF-1, neurotrophin receptor NTRK2 and multiple proteins linked to musculoskeletal and CNS development, that persisted for 6-months post-intervention.

INTERPRETATION

MDCF-2 treatment of Bangladeshi children with MAM, which produced significant improvements in wasting during intervention, also reduced stunting during follow-up. These results suggest that the effectiveness of supplementary foods for undernutrition may be improved by including ingredients that sponsor healthy microbiota-host co-development.

FUNDING

This work was supported by the BMGF (Grants OPP1134649/INV-000247).

X-chromosome inactivation in human iPSCs provides insight into X-regulated gene expression in autosomes

BACKGROUND

Variation in X chromosome inactivation (XCI) in human-induced pluripotent stem cells (hiPSCs) can impact their ability to model biological sex biases. The gene-wise landscape of X chromosome gene dosage remains unresolved in female hiPSCs. To characterize patterns of de-repression and escape from inactivation, we performed a systematic survey of allele specific expression in 165 female hiPSC lines.

RESULTS

XCI erosion is non-random and primarily affects genes that escape XCI in human tissues. Individual genes and cell lines vary in the frequency and degree of de-repression. Bi-allelic expression increases gradually after modest decrease of XIST in cultures, whose loss is commonly used to mark lines with eroded XCI. We identify three clusters of female lines at different stages of XCI. Increased XCI erosion amplifies female-biased expression at hypomethylated sites and regions normally occupied by repressive histone marks, lowering male-biased differences in the X chromosome. In autosomes, erosion modifies sex differences in a dose-dependent way. Male-biased genes are enriched for hypermethylated regions, and de-repression of XIST-bound autosomal genes in female lines attenuates normal male-biased gene expression in eroded lines. XCI erosion can compensate for a dominant loss of function effect in several disease genes.

CONCLUSIONS

We present a comprehensive view of X chromosome gene dosage in hiPSCs and implicate a direct mechanism for XCI erosion in regulating autosomal gene expression in trans. The uncommon and variable reactivation of X chromosome genes in female hiPSCs can provide insight into X chromosome's role in regulating gene expression and sex differences in humans.

Systems biology dissection of PTSD and MDD across brain regions, cell types, and blood

The molecular pathology of stress-related disorders remains elusive. Our brain multiregion, multiomic study of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) included the central nucleus of the amygdala, hippocampal dentate gyrus, and medial prefrontal cortex (mPFC). Genes and exons within the mPFC carried most disease signals replicated across two independent cohorts. Pathways pointed to immune function, neuronal and synaptic regulation, and stress hormones. Multiomic factor and gene network analyses provided the underlying genomic structure. Single nucleus RNA sequencing in dorsolateral PFC revealed dysregulated (stress-related) signals in neuronal and non-neuronal cell types. Analyses of brain-blood intersections in >50,000 UK Biobank participants were conducted along with fine-mapping of the results of PTSD and MDD genome-wide association studies to distinguish risk from disease processes. Our data suggest shared and distinct molecular pathology in both disorders and propose potential therapeutic targets and biomarkers.

Uncoupling interferons and the interferon signature explains clinical and transcriptional subsets in SLE

Systemic lupus erythematosus (SLE) displays a hallmark interferon (IFN) signature. Yet, clinical trials targeting type I IFN (IFN-I) have shown variable efficacy, and blocking IFN-II failed to treat SLE. Here, we show that IFN type levels in SLE vary significantly across clinical and transcriptional endotypes. Whereas skin involvement correlated with IFN-I alone, systemic features like nephritis associated with co-elevation of IFN-I, IFN-II, and IFN-III, indicating additive IFN effects in severe SLE. Notably, while high IFN-II/-III levels without IFN-I had a limited effect on disease activity, IFN-II was linked to IFN-I-independent transcriptional profiles (e.g., OXPHOS and CD8+GZMH+ cells), and IFN-III enhanced IFN-induced gene expression when co-elevated with IFN-I. Moreover, dysregulated IFNs do not explain the IFN signature in 64% of patients or clinical manifestations including cytopenia, serositis, and anti-phospholipid syndrome, implying IFN-independent endotypes in SLE. This study sheds light on mechanisms underlying SLE heterogeneity and the variable response to IFN-targeted therapies in clinical trials.

Multiomic foundations of human prefrontal cortex tissue function

The prefrontal cortex (PFC) is a region of the brain that in humans is involved in the production of higher-order functions such as cognition, emotion, perception, and behavior. Neurotransmission in the PFC produces higher-order functions by integrating information from other areas of the brain. At the foundation of neurotransmission, and by extension at the foundation of higher-order brain functions, are an untold number of coordinated molecular processes involving the DNA sequence variants in the genome, RNA transcripts in the transcriptome, and proteins in the proteome. These "multiomic" foundations are poorly understood in humans, perhaps in part because most modern studies that characterize the molecular state of the human PFC use tissue obtained when neurotransmission and higher-order brain functions have ceased (i.e., the postmortem state). Here, analyses are presented on data generated for the Living Brain Project (LBP) to investigate whether PFC tissue from individuals with intact higher-order brain function has characteristic multiomic foundations. Two complementary strategies were employed towards this end. The first strategy was to identify in PFC samples obtained from living study participants a signature of RNA transcript expression associated with neurotransmission measured intracranially at the time of PFC sampling, in some cases while participants performed a task engaging higher-order brain functions. The second strategy was to perform multiomic comparisons between PFC samples obtained from individuals with intact higher-order brain function at the time of sampling (i.e., living study participants) and PFC samples obtained in the postmortem state. RNA transcript expression within multiple PFC cell types was associated with fluctuations of dopaminergic, serotonergic, and/or noradrenergic neurotransmission in the substantia nigra measured while participants played a computer game that engaged higher-order brain functions. A subset of these associations - termed the "transcriptional program associated with neurotransmission" (TPAWN) - were reproduced in analyses of brain RNA transcript expression and intracranial neurotransmission data obtained from a second LBP cohort and from a cohort in an independent study. RNA transcripts involved in TPAWN were found to be (1) enriched for RNA transcripts associated with measures of neurotransmission in rodent and cell models, (2) enriched for RNA transcripts encoded by evolutionarily constrained genes, (3) depleted of RNA transcripts regulated by common DNA sequence variants, and (4) enriched for RNA transcripts implicated in higher-order brain functions by human population genetic studies. In PFC excitatory neurons of living study participants, higher expression of the genes in TPAWN tracked with higher expression of RNA transcripts that in rodent PFC samples are markers of a class of excitatory neurons that connect the PFC to deep brain structures. TPAWN was further reproduced by RNA transcript expression patterns differentiating living PFC samples from postmortem PFC samples, and significant differences between living and postmortem PFC samples were additionally observed with respect to (1) the expression of most primary RNA transcripts, mature RNA transcripts, and proteins, (2) the splicing of most primary RNA transcripts into mature RNA transcripts, (3) the patterns of co-expression between RNA transcripts and proteins, and (4) the effects of some DNA sequence variants on RNA transcript and protein expression. Taken together, this report highlights that studies of brain tissue obtained in a safe and ethical manner from large cohorts of living individuals can help advance understanding of the multiomic foundations of brain function.

Integrating population and single-cell variations in vaccine responses identifies a naturally adjuvanted human immune setpoint

Multimodal single-cell profiling methods can capture immune cell variations unfolding over time at the molecular, cellular, and population levels. Transforming these data into biological insights remains challenging. Here, we introduce a framework to integrate variations at the human population and single-cell levels in vaccination responses. Comparing responses following AS03-adjuvanted versus unadjuvanted influenza vaccines with CITE-seq revealed AS03-specific early (day 1) response phenotypes, including a B cell signature of elevated germinal center competition. A correlated network of cell-type-specific transcriptional states defined the baseline immune status associated with high antibody responders to the unadjuvanted vaccine. Certain innate subsets in the network appeared "naturally adjuvanted," with transcriptional states resembling those induced uniquely by AS03-adjuvanted vaccination. Consistently, CD14+ monocytes from high responders at baseline had elevated phospho-signaling responses to lipopolysaccharide stimulation. Our findings link baseline immune setpoints to early vaccine responses, with positive implications for adjuvant development and immune response engineering.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

JAK-STAT signaling maintains homeostasis in T cells and macrophages

Immune cells need to sustain a state of constant alertness over a lifetime. Yet, little is known about the regulatory processes that control the fluent and fragile balance that is called homeostasis. Here we demonstrate that JAK-STAT signaling, beyond its role in immune responses, is a major regulator of immune cell homeostasis. We investigated JAK-STAT-mediated transcription and chromatin accessibility across 12 mouse models, including knockouts of all STAT transcription factors and of the TYK2 kinase. Baseline JAK-STAT signaling was detected in CD8+ T cells and macrophages of unperturbed mice-but abrogated in the knockouts and in unstimulated immune cells deprived of their normal tissue context. We observed diverse gene-regulatory programs, including effects of STAT2 and IRF9 that were independent of STAT1. In summary, our large-scale dataset and integrative analysis of JAK-STAT mutant and wild-type mice uncovered a crucial role of JAK-STAT signaling in unstimulated immune cells, where it contributes to a poised epigenetic and transcriptional state and helps prepare these cells for rapid response to immune stimuli.

Multi-omics Analysis Reveals Immune Features Associated with Immunotherapy Benefit in Patients with Squamous Cell Lung Cancer from Phase III Lung-MAP S1400I Trial

PURPOSE

Identifying molecular and immune features to guide immune checkpoint inhibitor (ICI)-based regimens remains an unmet clinical need.

EXPERIMENTAL DESIGN

Tissue and longitudinal blood specimens from phase III trial S1400I in patients with metastatic squamous non-small cell carcinoma (SqNSCLC) treated with nivolumab monotherapy (nivo) or nivolumab plus ipilimumab (nivo+ipi) were subjected to multi-omics analyses including multiplex immunofluorescence (mIF), nCounter PanCancer Immune Profiling Panel, whole-exome sequencing, and Olink.

RESULTS

Higher immune scores from immune gene expression profiling or immune cell infiltration by mIF were associated with response to ICIs and improved survival, except regulatory T cells, which were associated with worse overall survival (OS) for patients receiving nivo+ipi. Immune cell density and closer proximity of CD8+GZB+ T cells to malignant cells were associated with superior progression-free survival and OS. The cold immune landscape of NSCLC was associated with a higher level of chromosomal copy-number variation (CNV) burden. Patients with LRP1B-mutant tumors had a shorter survival than patients with LRP1B-wild-type tumors. Olink assays revealed soluble proteins such as LAMP3 increased in responders while IL6 and CXCL13 increased in nonresponders. Upregulation of serum CXCL13, MMP12, CSF-1, and IL8 were associated with worse survival before radiologic progression.

CONCLUSIONS

The frequency, distribution, and clustering of immune cells relative to malignant ones can impact ICI efficacy in patients with SqNSCLC. High CNV burden may contribute to the cold immune microenvironment. Soluble inflammation/immune-related proteins in the blood have the potential to monitor therapeutic benefit from ICI treatment in patients with SqNSCLC.

Benchmark of cellular deconvolution methods using a multi-assay reference dataset from postmortem human prefrontal cortex

Background

Cellular deconvolution of bulk RNA-sequencing (RNA-seq) data using single cell or nuclei RNA-seq (sc/snRNA-seq) reference data is an important strategy for estimating cell type composition in heterogeneous tissues, such as human brain. Computational methods for deconvolution have been developed and benchmarked against simulated data, pseudobulked sc/snRNA-seq data, or immunohistochemistry reference data. A major limitation in developing improved deconvolution algorithms has been the lack of integrated datasets with orthogonal measurements of gene expression and estimates of cell type proportions on the same tissue sample. Deconvolution algorithm performance has not yet been evaluated across different RNA extraction methods (cytosolic, nuclear, or whole cell RNA), different library preparation types (mRNA enrichment vs. ribosomal RNA depletion), or with matched single cell reference datasets.

Results

A rich multi-assay dataset was generated in postmortem human dorsolateral prefrontal cortex (DLPFC) from 22 tissue blocks. Assays included spatially-resolved transcriptomics, snRNA-seq, bulk RNA-seq (across six library/extraction RNA-seq combinations), and RNAScope/Immunofluorescence (RNAScope/IF) for six broad cell types. The Mean Ratio method, implemented in the DeconvoBuddies R package, was developed for selecting cell type marker genes. Six computational deconvolution algorithms were evaluated in DLPFC and predicted cell type proportions were compared to orthogonal RNAScope/IF measurements.

Conclusions

Bisque and hspe were the most accurate methods, were robust to differences in RNA library types and extractions. This multi-assay dataset showed that cell size differences, marker genes differentially quantified across RNA libraries, and cell composition variability in reference snRNA-seq impact the accuracy of current deconvolution methods.

Sex differences in brain cell-type specific chromatin accessibility in schizophrenia

Our understanding of the sex-specific role of the non-coding genome in serious mental illness remains largely incomplete. To address this gap, we explored sex differences in 1,393 chromatin accessibility profiles, derived from neuronal and non-neuronal nuclei of two distinct cortical regions from 234 cases with serious mental illness and 235 controls. We identified sex-specific enhancer-promoter interactions and showed that they regulate genes involved in X-chromosome inactivation (XCI). Examining chromosomal conformation allowed us to identify sex-specific cis- and trans-regulatory domains (CRDs and TRDs). Co-localization of sex-specific TRDs with schizophrenia common risk variants pinpointed male-specific regulatory regions controlling a number of metabolic pathways. Additionally, enhancers from female-specific TRDs were found to regulate two genes known to escape XCI, (XIST and JPX), underlying the importance of TRDs in deciphering sex differences in schizophrenia. Overall, these findings provide extensive characterization of sex differences in the brain epigenome and disease-associated regulomes.

Genetic regulatory effects in response to a high-cholesterol, high-fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWASs), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using non-human primate models.

Comprehensive proteomics of CSF, plasma, and urine identify DDC and other biomarkers of early Parkinson's disease

Parkinson's disease (PD) starts at the molecular and cellular level long before motor symptoms appear, yet there are no early-stage molecular biomarkers for diagnosis, prognosis prediction, or monitoring therapeutic response. This lack of biomarkers greatly impedes patient care and translational research-L-DOPA remains the standard of care more than 50 years after its introduction. Here, we performed a large-scale, multi-tissue, and multi-platform proteomics study to identify new biomarkers for early diagnosis and disease monitoring in PD. We analyzed 4877 cerebrospinal fluid, blood plasma, and urine samples from participants across seven cohorts using three orthogonal proteomics methods: Olink proximity extension assay, SomaScan aptamer precipitation assay, and liquid chromatography-mass spectrometry proteomics. We discovered that hundreds of proteins were upregulated in the CSF, blood, or urine of PD patients, prodromal PD patients with DAT deficit and REM sleep behavior disorder or anosmia, and non-manifesting genetic carriers of LRRK2 and GBA mutations. We nominate multiple novel hits across our analyses as promising markers of early PD, including DOPA decarboxylase (DDC), also known as L-aromatic acid decarboxylase (AADC), sulfatase-modifying factor 1 (SUMF1), dipeptidyl peptidase 2/7 (DPP7), glutamyl aminopeptidase (ENPEP), WAP four-disulfide core domain 2 (WFDC2), and others. DDC, which catalyzes the final step in dopamine synthesis, particularly stands out as a novel hit with a compelling mechanistic link to PD pathogenesis. DDC is consistently upregulated in the CSF and urine of treatment-naïve PD, prodromal PD, and GBA or LRRK2 carrier participants by all three proteomics methods. We show that CSF DDC levels correlate with clinical symptom severity in treatment-naïve PD patients and can be used to accurately diagnose PD and prodromal PD. This suggests that urine and CSF DDC could be a promising diagnostic and prognostic marker with utility in both clinical care and translational research.

BHLHE40/41 regulate microglia and peripheral macrophage responses associated with Alzheimer's disease and other disorders of lipid-rich tissues

Genetic and experimental evidence suggests that Alzheimer's disease (AD) risk alleles and genes may influence disease susceptibility by altering the transcriptional and cellular responses of macrophages, including microglia, to damage of lipid-rich tissues like the brain. Recently, sc/nRNA sequencing studies identified similar transcriptional activation states in subpopulations of macrophages in aging and degenerating brains and in other diseased lipid-rich tissues. We collectively refer to these subpopulations of microglia and peripheral macrophages as DLAMs. Using macrophage sc/nRNA-seq data from healthy and diseased human and mouse lipid-rich tissues, we reconstructed gene regulatory networks and identified 11 strong candidate transcriptional regulators of the DLAM response across species. Loss or reduction of two of these transcription factors, BHLHE40/41, in iPSC-derived microglia and human THP-1 macrophages as well as loss of Bhlhe40/41 in mouse microglia, resulted in increased expression of DLAM genes involved in cholesterol clearance and lysosomal processing, increased cholesterol efflux and storage, and increased lysosomal mass and degradative capacity. These findings provide targets for therapeutic modulation of macrophage/microglial function in AD and other disorders affecting lipid-rich tissues.

Convergence of the dysregulated regulome in schizophrenia with polygenic risk and evolutionarily constrained enhancers

Enhancers play an essential role in the etiology of schizophrenia; however, the dysregulation of enhancer activity and its impact on the regulome in schizophrenia remains understudied. To address this gap in our knowledge, we assessed enhancer and gene expression in 1,382 brain samples comprising cases with schizophrenia and unaffected controls. Dysregulation of enhancer expression was concordant with changes in gene expression, and was more closely associated with schizophrenia polygenic risk, suggesting that enhancer dysregulation is proximal to the genetic etiology of the disease. Modeling the shared variance of cis-coordinated genes and enhancers revealed a gene regulatory program that was highly associated with genetic vulnerability to schizophrenia. By integrating coordinated factors with evolutionary constraints, we found that enhancers acquired during human evolution are more likely to regulate genes that are implicated in neuropsychiatric disorders and, thus, hold potential as therapeutic targets. Our analysis provides a systematic view of regulome dysregulation in schizophrenia and highlights its convergence with schizophrenia polygenic risk and human-gained enhancers.

Baseline Expression of Immune Gene Modules in Blood is Associated With Primary Response to Anti-TNF Therapy in Crohn's Disease Patients

BACKGROUND AND AIMS

Anti-tumour necrosis factor [anti-TNF] therapy is widely used for the treatment of inflammatory bowel disease, yet many patients are primary non-responders, failing to respond to induction therapy. We aimed to identify blood gene expression differences between primary responders and primary non-responders to anti-TNF monoclonal antibodies [infliximab and adalimumab], and to predict response status from blood gene expression and clinical data.

METHODS

The Personalised Anti-TNF Therapy in Crohn's Disease [PANTS] study is a UK-wide prospective observational cohort study of anti-TNF therapy outcome in anti-TNF-naive Crohn's disease patients [ClinicalTrials.gov identifier: NCT03088449]. Blood gene expression in 324 unique patients was measured by RNA-sequencing at baseline [week 0], and at weeks 14, 30, and 54 after treatment initiation [total sample size = 814].

RESULTS

After adjusting for clinical covariates and estimated blood cell composition, baseline expression of major histocompatibility complex, antigen presentation, myeloid cell enriched receptor, and other innate immune gene modules was significantly higher in anti-TNF responders vs non-responders. Expression changes from baseline to week 14 were generally of consistent direction but greater magnitude [i.e. amplified] in responders, but interferon-related genes were upregulated uniquely in non-responders. Expression differences between responders and non-responders observed at week 14 were maintained at weeks 30 and 54. Prediction of response status from baseline clinical data, cell composition, and module expression was poor.

CONCLUSIONS

Baseline gene module expression was associated with primary response to anti-TNF therapy in PANTS patients. However, these baseline expression differences did not predict response with sufficient sensitivity for clinical use.

WITHDRAWN: The effect of ischemia on expression quantitative trait loci (eQTL) in human myocardium and insights into myocardial injury etiology

26 February, 2024. Research Square has withdrawn this preprint as it was submitted and made public without the full consent of all the authors and without the full consent of the principle investigator of the registered clinical trial. Therefore, this work should not be cited as a reference.

A host-microbial metabolite interaction gut-on-a-chip model of the adult human intestine demonstrates beneficial effects upon inulin treatment of gut microbiome

Background: The gut and its microbiome have a major impact on many aspects of health and are therefore also an attractive target for drug- or food-based therapies. Here, we report on the added value of combining a microbiome screening model, the i-screen, with fresh intestinal tissue explants in a microfluidic gut-on-a-chip model, the Intestinal Explant Barrier Chip (IEBC). Methods: Adult human gut microbiome (fecal pool of 6 healthy donors) was cultured anaerobically in the i-screen platform for 24 h, without and with exposure to 4 mg/mL inulin. The i-screen cell-free culture supernatant was subsequently applied to the luminal side of adult human colon tissue explants (n = 3 donors), fixed in the IEBC, for 24 h and effects were evaluated. Results: The supplementation of the media with inulin promoted the growth of Anaerostipes, Bifidobacterium, Blautia, and Collinsella in the in vitro i-screen, and triggered an elevated production of butyrate by the microbiota. Human colon tissue exposed to inulin-treated i-screen cell-free culture supernatant or control i-screen cell-free culture supernatant with added short-chain fatty acids (SCFAs) showed improved tissue barrier integrity measured by a 28.2%-34.2% reduction in FITC-dextran 4000 (FD4) leakage and 1.3 times lower transport of antipyrine. Furthermore, the release of pro-inflammatory cytokines IL-1β, IL-6, IL-8, and TNF-α was reduced under these circumstances. Gene expression profiles confirmed these findings, but showed more profound effects for inulin-treated supernatant compared to SCFA-supplemented supernatant. Conclusion: The combination of i-screen and IEBC facilitates the study of complex intestinal processes such as host-microbial metabolite interaction and gut health.

Alzheimer's disease transcriptional landscape in ex-vivo human microglia

Microglia are resident immune cells of the brain and are implicated in the etiology of Alzheimer's Disease (AD) and other diseases. Yet the cellular and molecular processes regulating their function throughout the course of the disease are poorly understood. Here, we present the transcriptional landscape of primary microglia from 189 human postmortem brains, including 58 healthy aging individuals and 131 with a range of disease phenotypes, including 63 patients representing the full spectrum of clinical and pathological severity of AD. We identified transcriptional changes associated with multiple AD phenotypes, capturing the severity of dementia and neuropathological lesions. Transcript-level analyses identified additional genes with heterogeneous isoform usage and AD phenotypes. We identified changes in gene-gene coordination in AD, dysregulation of co-expression modules, and disease subtypes with distinct gene expression. Taken together, these data further our understanding of the key role of microglia in AD biology and nominate candidates for therapeutic intervention.

Multi-tissue metabolic and transcriptomic responses to a short-term heat stress in swine

BACKGROUND

Heat stress (HS) is an increasing threat for pig production with a wide range of impacts. When submitted to high temperatures, pigs will use a variety of strategies to alleviate the effect of HS. While systemic adaptations are well known, tissue-specific changes remain poorly understood. In this study, thirty-two pigs were submitted to a 5-day HS at 32 °C.

RESULTS

Transcriptomic and metabolomic analyses were performed on several tissues. The results revealed differentially expressed genes and metabolites in different tissues. Specifically, 481, 1774, 71, 1572, 17, 164, and 169 genes were differentially expressed in muscle, adipose tissue, liver, blood, thyroid, pituitary, and adrenal glands, respectively. Regulatory glands (pituitary, thyroid, and adrenal) had a lower number of regulated genes, perhaps indicating an earlier sensitivity to HS. In addition, 7, 8, 2, and 8 metabolites were differentially produced in muscle, liver, plasma, and urine, respectively. The study also focused on the oxidative stress pathway in muscle and liver by performing a correlation analysis between genes and metabolites.

CONCLUSIONS

This study has identified various adaptation mechanisms in swine that enable them to cope with heat stress (HS). These mechanisms include a global decrease in energetic metabolism, as well as changes in metabolic precursors that are linked with protein and lipid catabolism and anabolism. Notably, the adaptation mechanisms differ significantly between regulatory (pituitary, thyroid and adrenal glands) and effector tissues (muscle, adipose tissue, liver and blood). Our findings provide new insights into the comprehension of HS adaptation mechanisms in swine.

Activity-Dependent Transcriptional Program in NGN2+ Neurons Enriched for Genetic Risk for Brain-Related Disorders

BACKGROUND

Converging evidence from large-scale genetic and postmortem studies highlights the role of aberrant neurotransmission and genetic regulation in brain-related disorders. However, identifying neuronal activity-regulated transcriptional programs in the human brain and understanding how changes contribute to disease remain challenging.

METHODS

To better understand how the activity-dependent regulome contributes to risk for brain-related disorders, we profiled the transcriptomic and epigenomic changes following neuronal depolarization in human induced pluripotent stem cell-derived glutamatergic neurons (NGN2) from 6 patients with schizophrenia and 5 control participants.

RESULTS

Multiomic data integration associated global patterns of chromatin accessibility with gene expression and identified enhancer-promoter interactions in glutamatergic neurons. Within 1 hour of potassium chloride-induced depolarization, independent of diagnosis, glutamatergic neurons displayed substantial activity-dependent changes in the expression of genes regulating synaptic function. Depolarization-induced changes in the regulome revealed significant heritability enrichment for schizophrenia and Parkinson's disease, adding to mounting evidence that sequence variation within activation-dependent regulatory elements contributes to the genetic risk for brain-related disorders. Gene coexpression network analysis elucidated interactions among activity-dependent and disease-associated genes and pointed to a key driver (NAV3) that interacted with multiple genes involved in axon guidance.

CONCLUSIONS

Overall, we demonstrated that deciphering the activity-dependent regulome in glutamatergic neurons reveals novel targets for advanced diagnosis and therapy.

Bioactive glycans in a microbiome-directed food for children with malnutrition

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

The documentation of Plasmodium vivax malaria across Africa especially in regions where Duffy negatives are dominant suggests possibly alternative erythrocyte invasion mechanisms. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa is unclear. In this study, we examined the expression of 4,404 gene transcripts belong to 12 functional groups and 43 erythrocyte binding gene candidates in Ethiopian isolates and compared them with the Cambodian and Brazilian P. vivax transcriptomes. Overall, there were 10-26% differences in the gene expression profile amongst geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian isolates. Members of the reticulocyte binding protein PvRBP2a and PvRBP3 expressed six-fold higher than Duffy binding protein PvDBP1 and 60-fold higher than PvEBP/DBP2 in the Ethiopian isolates. Other genes including PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression patterns were observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, gametocyte genes including PvAP2-G, PvGAP (female gametocytes), and Pvs47 (male gametocytes) were highly expressed across geographical samples.

Decoding glycomics with a suite of methods for differential expression analysis

Glycomics, the comprehensive profiling of all glycan structures in samples, is rapidly expanding to enable insights into physiology and disease mechanisms. However, glycan structure complexity and glycomics data interpretation present challenges, especially for differential expression analysis. Here, we present a framework for differential glycomics expression analysis. Our methodology encompasses specialized and domain-informed methods for data normalization and imputation, glycan motif extraction and quantification, differential expression analysis, motif enrichment analysis, time series analysis, and meta-analytic capabilities, synthesizing results across multiple studies. All methods are integrated into our open-source glycowork package, facilitating performant workflows and user-friendly access. We demonstrate these methods using dedicated simulations and glycomics datasets of N-, O-, lipid-linked, and free glycans. Differential expression tests here focus on human datasets and cancer vs. healthy tissue comparisons. Our rigorous approach allows for robust, reliable, and comprehensive differential expression analyses in glycomics, contributing to advancing glycomics research and its translation to clinical and diagnostic applications.

Different Genetic Signatures of Small-Cell Lung Cancer Characterize Anti-GABAB R and Anti-Hu Paraneoplastic Neurological Syndromes

OBJECTIVE

Small-cell lung cancer (SCLC) is the malignancy most frequently associated with paraneoplastic neurological syndromes (PNS) and can trigger different antibody responses against intracellular (Hu) or neuronal surface (GABAB R) antigens. Our aim was to clarify whether the genomic and transcriptomic features of SCLC are different in patients with anti-GABAB R or anti-Hu PNS compared with SCLC without PNS.

METHODS

A total of 76 SCLC tumor samples were collected: 34 anti-Hu, 14 anti-GABAB R, and 28 SCLC without PNS. The study consisted of 4 steps: (1) pathological confirmation; (2) next generation sequencing using a panel of 98 genes, including those encoding the autoantibodies targets ELAVL1-4, GABBR1-2, and KCTD16; (3) genome-wide copy number variation (CNV); and (4) whole-transcriptome RNA sequencing.

RESULTS

CNV analysis revealed that patients with anti-GABAB R PNS commonly have a gain in chromosome 5q, which contains KCTD16, whereas anti-Hu and control patients often harbor a loss. No significantly different number of mutations regarding any onconeural genes was observed. Conversely, the transcriptomic profile of SCLC was different, and the differentially expressed genes allowed effective clustering of the samples into 3 groups, reflecting the antibody-based classification, with an overexpression of KCTD16 specific to anti-GABAB R PNS. Pathway analysis revealed that tumors of patients with anti-GABAB R encephalitis were enriched in B-cell signatures, as opposed to those of patients with anti-Hu, in which T-cell- and interferon-γ-related signatures were overexpressed.

INTERPRETATION

SCLC genetic and transcriptomic features differentiate anti-GABAB R, anti-Hu, and non-PNS tumors. The role of KCTD16 appears to be pivotal in the tumor immune tolerance breakdown of anti-GABAB R PNS. ANN NEUROL 2023;94:1102-1115.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

RNA-sequencing of peripheral whole blood of individuals at ultra-high-risk for psychosis - A longitudinal perspective

BACKGROUND

The peripheral blood is an attractive source of prognostic biomarkers for psychosis conversion. There is limited research on the transcriptomic changes associated with psychosis conversion in the peripheral whole blood.

STUDY DESIGN

We performed RNA-sequencing of peripheral whole blood from 65 ultra-high-risk (UHR) participants and 70 healthy control participants recruited in the Longitudinal Youth-at-Risk Study (LYRIKS) cohort. 13 UHR participants converted in the study duration. Samples were collected at 3 timepoints, at 12-months interval across a 2-year period. We examined whether the genes differential with psychosis conversion contain schizophrenia risk loci. We then examined the functional ontologies and GWAS associations of the differential genes. We also identified the overlap between differentially expressed genes across different comparisons.

STUDY RESULTS

Genes containing schizophrenia risk loci were not differentially expressed in the peripheral whole blood in psychosis conversion. The differentially expressed genes in psychosis conversion are enriched for ontologies associated with cellular replication. The differentially expressed genes in psychosis conversion are associated with non-neurological GWAS phenotypes reported to be perturbed in schizophrenia and psychosis but not schizophrenia and psychosis phenotypes themselves. We found minimal overlap between the genes differential with psychosis conversion and the genes that are differential between pre-conversion and non-conversion samples.

CONCLUSION

The associations between psychosis conversion and peripheral blood-based biomarkers are likely to be indirect. Further studies to elucidate the mechanism behind potential indirect associations are needed.

The Effect of Dietary Advice Aimed at Increasing Protein Intake on Oral Health and Oral Microbiota in Older Adults: A Randomized Controlled Trial

Nutrition and oral health are closely related, especially in older adults in whom poor nutrition may lead to oral microbial perturbations, exacerbating poor oral health. In a 6-month randomized controlled trial, we evaluated the effects on oral microbiota and on oral health of dietary advice aimed at increasing protein intake to ≥1.2 g/kg adjusted body weight/day (g/kg aBW/d) in community-dwelling older adults with low habitual protein intake (<1.0 g/kg aBW/d). Food intake was measured via 24 h dietary recalls, oral health was measured via questionnaires, and oral microbial composition was assessed via the 16S rRNA sequencing of tongue swabs. Mean baseline protein intake was 0.8 g/kg aBW/day in both groups. In the high protein group (n = 47), participants increased their protein intake to mean 1.2 g/kg aBW/day at the 6-month follow-up. Protein intake in the control group (n = 43) remained at 0.9 g/kg a BW/day. The intervention did not affect self-reported oral health. While it caused moderate shifts in oral microbiota alpha- and beta-diversity measures, abundances of individual bacterial taxa were not affected. In conclusion, our intervention did not affect self-reported oral health within a period of 6 months, nor did it substantially affect the tongue microbiota composition.

Unravelling the sexual developmental biology of Cystoisospora suis, a model for comparative coccidian parasite studies

Introduction

The apicomplexan parasite Cystoisospora suis has global significance as an enteropathogen of suckling piglets. Its intricate life cycle entails a transition from an asexual phase to sexual development, ultimately leading to the formation of transmissible oocysts.

Methods

To advance our understanding of the parasite's cellular development, we complemented previous transcriptome studies by delving into the proteome profiles at five distinct time points of in vitro cultivation through LC/MS-MS analysis.

Results

A total of 1,324 proteins were identified in the in vitro developmental stages of C. suis, and 1,082 proteins were identified as significantly differentially expressed. Data are available via ProteomeXchange with identifier PXD045050. We performed BLAST, GO enrichment, and KEGG pathway analyses on the up- and downregulated proteins to elucidate correlated events in the C. suis life cycle. Our analyses revealed intriguing metabolic patterns in macromolecule metabolism, DNA- and RNA-related processes, proteins associated with sexual stages, and those involved in cell invasion, reflecting the adaptation of sexual stages to a nutrient-poor and potentially stressful extracellular environment, with a focus on enzymes involved in metabolism and energy production.

Discussion

These findings have important implications for understanding the developmental biology of C. suis as well as other, related coccidian parasites, such as Eimeria spp. and Toxoplasma gondii. They also support the role of C. suis as a new model for the comparative biology of coccidian tissue cyst stages.

The human microglia responsome: a resource to better understand microglia states in health and disease

Microglia, the immune cells of the brain, are increasingly implicated in neurodegenerative disorders through genetic studies. However, how genetic risk factors for these diseases are related to microglial gene expression, microglial function, and ultimately disease, is still largely unknown. Microglia change rapidly in response to alterations in their cellular environment, which is regulated through changes in transcriptional programs, which are as yet poorly understood. Here, we compared the effects of a set of inflammatory and restorative stimuli (lipopolysaccharide, interferon-gamma, resiquimod, tumor necrosis factor-alpha, adenosine triphosphate, dexamethasone, and interleukin-4) on human microglial cells from 67 different donors (N = 398 samples) at the gene and transcript level. We show that microglia from different anatomical brain regions show distinct responses to inflammatory stimuli. We observed a greater overlap between human stimulated microglia and human monocytes than with mouse microglia. We define specific microglial signatures across conditions which are highly relevant for a wide range of biological functions and complex human diseases. Finally, we used our stimulation signatures to interpret associations from Alzheimer's disease (AD) genetic studies and microglia by integrating our inflammatory gene expression profiles with common genetic variants to map cis -expression QTLs (eQTLs). Together, we provide the most comprehensive transcriptomic database of the human microglia responsome.

Highlights

RNA-sequencing of 398 human microglial samples exposed to six different triggers.Microglia from different anatomical regions show distinct stimulation responses.Responses in human microglia show a greater overlap with human monocytes than murine microglia.Mapping of response Quantitative Trait Loci identifies interactions between genotype and effect of stimulation on gene expression.Our atlas provides a reference map for interpreting microglia signatures in health and disease.

The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

The relationship between regulatory changes in cis and trans and the evolution of gene expression in humans and chimpanzees

BACKGROUND

Comparative gene expression studies in apes are fundamentally limited by the challenges associated with sampling across different tissues. Here, we used single-cell RNA sequencing of embryoid bodies to collect transcriptomic data from over 70 cell types in three humans and three chimpanzees.

RESULTS

We find hundreds of genes whose regulation is conserved across cell types, as well as genes whose regulation likely evolves under directional selection in one or a handful of cell types. Using embryoid bodies from a human-chimpanzee fused cell line, we also infer the proportion of inter-species regulatory differences due to changes in cis and trans elements between the species. Using the cis/trans inference and an analysis of transcription factor binding sites, we identify dozens of transcription factors whose inter-species differences in expression are affecting expression differences between humans and chimpanzees in hundreds of target genes.

CONCLUSIONS

Here, we present the most comprehensive dataset of comparative gene expression from humans and chimpanzees to date, including a catalog of regulatory mechanisms associated with inter-species differences.

Viruses in Laboratory Drosophila and Their Impact on Host Gene Expression

Drosophila melanogaster has one of the best characterized antiviral immune responses among invertebrates. However, relatively few easily transmitted natural virus isolates are available, and so many Drosophila experiments have been performed using artificial infection routes and artificial host-virus combinations. These may not reflect natural infections, especially for subtle phenotypes such as gene expression. Here, to explore the laboratory virus community and to better understand how natural virus infections induce changes in gene expression, we have analysed seven publicly available D. melanogaster transcriptomic sequencing datasets that were originally sequenced for projects unrelated to virus infection. We have found ten known viruses-including five that have not been experimentally isolated-but no previously unknown viruses. Our analysis of host gene expression revealed that numerous genes were differentially expressed in flies that were naturally infected with a virus. For example, flies infected with nora virus showed patterns of gene expression consistent with intestinal vacuolization and possible host repair via the upd3 JAK/STAT pathway. We also found marked sex differences in virus-induced differential gene expression. Our results show that natural virus infection in laboratory Drosophila does indeed induce detectable changes in gene expression, suggesting that this may form an important background condition for experimental studies in the laboratory.

Uncoupling interferons and the interferon signature explain clinical and transcriptional subsets in SLE

Interferons (IFN) are thought to be key players in systemic lupus erythematosus (SLE). The unique and interactive roles of the different IFN families in SLE pathogenesis, however, remain poorly understood. Using reporter cells engineered to precisely quantify IFN-I, IFN-II and IFN-III activity levels in serum/plasma, we found that while IFNs play essential role in SLE pathogenesis and disease activity, they are only significant in specific subsets of patients. Interestingly, whereas IFN-I is the main IFN that governs disease activity in SLE, clinical subsets are defined by the co-elevation of IFN-II and IFN-III. Thus, increased IFN-I alone was only associated with cutaneous lupus. In contrast, systemic features, such as nephritis, were linked to co-elevation of IFN-I plus IFN-II and IFN-III, implying a synergistic effect of IFNs in severe SLE. Intriguingly, while increased IFN-I levels were strongly associated with IFN-induced gene expression (93.5%), in up to 64% of cases, the IFN signature was not associated with IFN-I. Importantly, neither IFN-II nor IFN-III explained IFN-induced gene expression in patients with normal IFN-I levels, and not every feature in SLE was associated with elevated IFNs, suggesting IFN-independent subsets in SLE. Together, the data suggest that, unlike the IFN signature, direct quantification of bioactive IFNs can identify pathogenic and clinically relevant SLE subsets amenable for precise anti-IFN therapies. Since IFN-I is only elevated in a subset of SLE patients expressing the IFN signature, this study explains the heterogeneous response in clinical trials targeting IFN-I, where patients were selected based on IFN-induced gene expression rather than IFN-I levels.

Individual and Group-Based Effects of In Vitro Fiber Interventions on the Fecal Microbiota

The development of microbiome-targeted strategies is limited by individual differences in gut microbiome composition and metabolic responses to interventions. In vitro models that can replicate this variation allow us to conduct pre-clinical studies and assess efficacy. This study describes the exposure of 16 individual fecal microbiota samples to 5 different fibers using an in vitro system for the anaerobic cultivation of bacteria. The individual microbiota differed in composition and metabolite profiles (short-chain fatty acids and branched-chain fatty acids) after incubation with the fibers. Furthermore, microbiota composition after fiber incubation was significantly different between subjects with good intestinal health and subjects with Inflammatory Bowel Disease (IBD). α-diversity was differently affected by dietary fibers; for example, exposure to psyllium resulted in increased diversity in the healthy group and in decreased diversity in the IBD group. Instead, the functional metabolic profile did not differ between the two groups. Finally, the combination of all fibers, tested on the microbiota from IBD subjects, resulted in stronger overall effects on both microbiota composition and metabolite production compared to the single fibers. These results confirm that incubation with dietary fiber results in different compositional and functional effects on individual microbiota and that in vitro models represent successful tools for studying individual fiber effects.

A study of gene expression in the living human brain

A goal of medical research is to determine the molecular basis of human brain health and illness. One way to achieve this goal is through observational studies of gene expression in human brain tissue. Due to the unavailability of brain tissue from living people, most such studies are performed using tissue from postmortem brain donors. An assumption underlying this practice is that gene expression in the postmortem human brain is an accurate representation of gene expression in the living human brain. Here, this assumption - which, until now, had not been adequately tested - is tested by comparing human prefrontal cortex gene expression between 275 living samples and 243 postmortem samples. Expression levels differed significantly for nearly 80% of genes, and a systematic examination of alternative explanations for this observation determined that these differences are not a consequence of cell type composition, RNA quality, postmortem interval, age, medication, morbidity, symptom severity, tissue pathology, sample handling, batch effects, or computational methods utilized. Analyses integrating the data generated for this study with data from earlier landmark studies that used tissue from postmortem brain donors showed that postmortem brain gene expression signatures of neurological and mental illnesses, as well as of normal traits such as aging, may not be accurate representations of these gene expression signatures in the living brain. By using tissue from large cohorts living people, future observational studies of human brain biology have the potential to (1) determine the medical research questions that can be addressed using postmortem tissue as a proxy for living tissue and (2) expand the scope of medical research to include questions about the molecular basis of human brain health and illness that can only be addressed in living people (e.g., "What happens at the molecular level in the brain as a person experiences an emotion?").

Treatment-associated mRNA co-expression changes in monocytes of patients with posttraumatic stress disorder

PTSD is a prevalent mental disorder that results from exposure to extreme and stressful life events and comes at high costs for both the individual and society. Therapeutic treatment presents the best way to deal with PTSD-the mechanisms underlying change after treatment, however, remain poorly understood. While stress and immune associated gene expression changes have been associated with PTSD development, studies investigating treatment effects at the molecular level so far tended to focus on DNA methylation. Here we use gene-network analysis on whole-transcriptome RNA-Seq data isolated from CD14⁺ monocytes of female PTSD patients (N = 51) to study pre-treatment signatures of therapy response and therapy-related changes at the level of gene expression. Patients who exhibited significant symptom improvement after therapy showed higher baseline expression in two modules involved in inflammatory processes (including notable examples IL1R2 and FKBP5) and blood coagulation. After therapy, expression of an inflammatory module was increased, and expression of a wound healing module was decreased. This supports findings reporting an association between PTSD and dysregulations of the inflammatory and the hemostatic system and mark both as potentially treatment sensitive.

Chronic airway epithelial hypoxia exacerbates injury in muco-obstructive lung disease through mucus hyperconcentration

Unlike solid organs, human airway epithelia derive their oxygen from inspired air rather than the vasculature. Many pulmonary diseases are associated with intraluminal airway obstruction caused by aspirated foreign bodies, virus infection, tumors, or mucus plugs intrinsic to airway disease, including cystic fibrosis (CF). Consistent with requirements for luminal O2, airway epithelia surrounding mucus plugs in chronic obstructive pulmonary disease (COPD) lungs are hypoxic. Despite these observations, the effects of chronic hypoxia (CH) on airway epithelial host defense functions relevant to pulmonary disease have not been investigated. Molecular characterization of resected human lungs from individuals with a spectrum of muco-obstructive lung diseases (MOLDs) or COVID-19 identified molecular features of chronic hypoxia, including increased EGLN3 expression, in epithelia lining mucus-obstructed airways. In vitro experiments using cultured chronically hypoxic airway epithelia revealed conversion to a glycolytic metabolic state with maintenance of cellular architecture. Chronically hypoxic airway epithelia unexpectedly exhibited increased MUC5B mucin production and increased transepithelial Na+ and fluid absorption mediated by HIF1α/HIF2α-dependent up-regulation of β and γENaC (epithelial Na+ channel) subunit expression. The combination of increased Na+ absorption and MUC5B production generated hyperconcentrated mucus predicted to perpetuate obstruction. Single-cell and bulk RNA sequencing analyses of chronically hypoxic cultured airway epithelia revealed transcriptional changes involved in airway wall remodeling, destruction, and angiogenesis. These results were confirmed by RNA-in situ hybridization studies of lungs from individuals with MOLD. Our data suggest that chronic airway epithelial hypoxia may be central to the pathogenesis of persistent mucus accumulation in MOLDs and associated airway wall damage.

Development of MAPT S305 mutation models exhibiting elevated 4R tau expression, resulting in altered neuronal and astrocytic function

Due to the importance of 4R tau in the pathogenicity of primary tauopathies, it has been challenging to model these diseases in iPSC-derived neurons, which express very low levels of 4R tau. To address this problem we have developed a panel of isogenic iPSC lines carrying the MAPT splice-site mutations S305S, S305I or S305N, derived from four different donors. All three mutations significantly increased the proportion of 4R tau expression in iPSC-neurons and astrocytes, with up to 80% 4R transcripts in S305N neurons from as early as 4 weeks of differentiation. Transcriptomic and functional analyses of S305 mutant neurons revealed shared disruption in glutamate signaling and synaptic maturity, but divergent effects on mitochondrial bioenergetics. In iPSC-astrocytes, S305 mutations induced lysosomal disruption and inflammation and exacerbated internalization of exogenous tau that may be a precursor to the glial pathologies observed in many tauopathies. In conclusion, we present a novel panel of human iPSC lines that express unprecedented levels of 4R tau in neurons and astrocytes. These lines recapitulate previously characterized tauopathy-relevant phenotypes, but also highlight functional differences between the wild type 4R and mutant 4R proteins. We also highlight the functional importance of MAPT expression in astrocytes. These lines will be highly beneficial to tauopathy researchers enabling a more complete understanding of the pathogenic mechanisms underlying 4R tauopathies across different cell types.

Host transcriptomic plasticity and photosymbiotic fidelity underpin Pocillopora acclimatization across thermal regimes in the Pacific Ocean

Heat waves are causing declines in coral reefs globally. Coral thermal responses depend on multiple, interacting drivers, such as past thermal exposure, endosymbiont community composition, and host genotype. This makes the understanding of their relative roles in adaptive and/or plastic responses crucial for anticipating impacts of future warming. Here, we extracted DNA and RNA from 102 Pocillopora colonies collected from 32 sites on 11 islands across the Pacific Ocean to characterize host-photosymbiont fidelity and to investigate patterns of gene expression across a historical thermal gradient. We report high host-photosymbiont fidelity and show that coral and microalgal gene expression respond to different drivers. Differences in photosymbiotic association had only weak impacts on host gene expression, which was more strongly correlated with the historical thermal environment, whereas, photosymbiont gene expression was largely determined by microalgal lineage. Overall, our results reveal a three-tiered strategy of thermal acclimatization in Pocillopora underpinned by host-photosymbiont specificity, host transcriptomic plasticity, and differential photosymbiotic association under extreme warming.

Kimma: flexible linear mixed effects modeling with kinship covariance for RNA-seq data

MOTIVATION

The identification of differentially expressed genes (DEGs) from transcriptomic datasets is a major avenue of research across diverse disciplines. However, current bioinformatic tools do not support covariance matrices in DEG modeling. Here, we introduce kimma (Kinship In Mixed Model Analysis), an open-source R package for flexible linear mixed effects modeling including covariates, weights, random effects, covariance matrices, and fit metrics.

RESULTS

In simulated datasets, kimma detects DEGs with similar specificity, sensitivity, and computational time as limma unpaired and dream paired models. Unlike other software, kimma supports covariance matrices as well as fit metrics like Akaike information criterion (AIC). Utilizing genetic kinship covariance, kimma revealed that kinship impacts model fit and DEG detection in a related cohort. Thus, kimma equals or outcompetes current DEG pipelines in sensitivity, computational time, and model complexity.

AVAILABILITY AND IMPLEMENTATION

Kimma is freely available on GitHub https://github.com/BIGslu/kimma with an instructional vignette at https://bigslu.github.io/kimma_vignette/kimma_vignette.html.

Efficient differential expression analysis of large-scale single cell transcriptomics data using dreamlet

Advances in single-cell and -nucleus transcriptomics have enabled generation of increasingly large-scale datasets from hundreds of subjects and millions of cells. These studies promise to give unprecedented insight into the cell type specific biology of human disease. Yet performing differential expression analyses across subjects remains difficult due to challenges in statistical modeling of these complex studies and scaling analyses to large datasets. Our open-source R package dreamlet (DiseaseNeurogenomics.github.io/dreamlet) uses a pseudobulk approach based on precision-weighted linear mixed models to identify genes differentially expressed with traits across subjects for each cell cluster. Designed for data from large cohorts, dreamlet is substantially faster and uses less memory than existing workflows, while supporting complex statistical models and controlling the false positive rate. We demonstrate computational and statistical performance on published datasets, and a novel dataset of 1.4M single nuclei from postmortem brains of 150 Alzheimer's disease cases and 149 controls.

Comparative transcriptomics reveal differential gene expression in Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

Plasmodium vivax uses Duffy binding protein (PvDBP1) to bind to the Duffy Antigen-Chemokine Receptor (DARC) to invade human erythrocytes. Individuals who lack DARC expression (Duffy-negative) are thought to be resistance to P. vivax. In recent years, P. vivax malaria is becoming more prevalent in Africa with a portion of these cases detected in Duffy-negatives. Apart from DBP1, members of the reticulocyte binding protein (RBP) and tryptophan-rich antigen (TRAg) families may also play a role in erythrocyte invasion. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa and more specifically the expression level of several erythrocyte binding gene candidates as compared to DBP1 are largely unknown. This paper characterized the first P. vivax transcriptome in Africa and compared with those from the Southeast Asian and South American isolates. The expression of 4,404 gene transcripts belong to 12 functional groups including 43 specific erythrocyte binding gene candidates were examined. Overall, there were 10-26% differences in the gene expression profile amongst the geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian P. vivax. Erythrocyte binding genes including PvRBP2a and PvRBP3 expressed six-fold higher than PvDBP1and 60-fold higher than PvEBP/DBP2. Other genes including PvRBP1a, PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression was observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, the standard biomarker for detecting female gametocytes, PvAP2-G (PVP01_1440800), GAP (PVP01_1403000), and Pvs47 (PVP01_1208000) were highly expressed across geographical samples. These findings provide an important baseline for future comparisons of P. vivax transcriptomes from Duffy-negative infections and highlight potential biomarkers for improved gametocyte detection.

BHLHE40/41 regulate macrophage/microglia responses associated with Alzheimer's disease and other disorders of lipid-rich tissues

Background

Genetic and experimental evidence strongly implicates myeloid cells in the etiology of AD and suggests that AD-associated alleles and genes may modulate disease risk by altering the transcriptional and cellular responses of macrophages (like microglia) to damage of lipid-rich tissues (like the brain). Specifically, recent single-cell/nucleus RNA sequencing (sc/nRNA-seq) studies identified a transcriptionally distinct state of subsets of macrophages in aging or degenerating brains (usually referred to as disease-associated microglia or DAM) and in other diseased lipid-rich tissues (e.g., obese adipose tissue, fatty liver, and atherosclerotic plaques). We collectively refer to these subpopulations as lipid-associated macrophages or LAMs. Importantly, this particular activation state is characterized by increased expression of genes involved in the phagocytic clearance of lipid-rich cellular debris (efferocytosis), including several AD risk genes.

Methods

We used sc/nRNA-seq data from human and mouse microglia from healthy and diseased brains and macrophages from other lipid-rich tissues to reconstruct gene regulatory networks and identify transcriptional regulators whose regulons are enriched for LAM response genes (LAM TFs) across species. We then used gene knock-down/knock-out strategies to validate some of these LAM TFs in human THP-1 macrophages and iPSC-derived microglia in vitro, as well as mouse microglia in vivo.

Results

We nominate 11 strong candidate LAM TFs shared across human and mouse networks (BHLHE41, HIF1A, ID2, JUNB, MAF, MAFB, MEF2A, MEF2C, NACA, POU2F2 and SPI1). We also demonstrate a strong enrichment of AD risk alleles in the cistrome of BHLHE41 (and its close homolog BHLHE40), thus implicating its regulon in the modulation of disease susceptibility. Loss or reduction of BHLHE40/41 expression in human THP-1 macrophages and iPSC-derived microglia, as well as loss of Bhlhe40/41 in mouse microglia led to increased expression of LAM response genes, specifically those involved in cholesterol clearance and lysosomal processing, with a concomitant increase in cholesterol efflux and storage, as well as lysosomal mass and degradative capacity.

Conclusions

Taken together, this study nominates transcriptional regulators of the LAM response, experimentally validates BHLHE40/41 in human and mouse macrophages/microglia, and provides novel targets for therapeutic modulation of macrophage/microglia function in AD and other disorders of lipid-rich tissues.

Influenza vaccination reveals sex dimorphic imprints of prior mild COVID-19

Acute viral infections can have durable functional impacts on the immune system long after recovery, but how they affect homeostatic immune states and responses to future perturbations remain poorly understood1-4. Here we use systems immunology approaches, including longitudinal multimodal single-cell analysis (surface proteins, transcriptome and V(D)J sequences) to comparatively assess baseline immune statuses and responses to influenza vaccination in 33 healthy individuals after recovery from mild, non-hospitalized COVID-19 (mean, 151 days after diagnosis) and 40 age- and sex-matched control individuals who had never had COVID-19. At the baseline and independent of time after COVID-19, recoverees had elevated T cell activation signatures and lower expression of innate immune genes including Toll-like receptors in monocytes. Male individuals who had recovered from COVID-19 had coordinately higher innate, influenza-specific plasmablast, and antibody responses after vaccination compared with healthy male individuals and female individuals who had recovered from COVID-19, in part because male recoverees had monocytes with higher IL-15 responses early after vaccination coupled with elevated prevaccination frequencies of 'virtual memory'-like CD8+ T cells poised to produce more IFNγ after IL-15 stimulation. Moreover, the expression of the repressed innate immune genes in monocytes increased by day 1 to day 28 after vaccination in recoverees, therefore moving towards the prevaccination baseline of the healthy control individuals. By contrast, these genes decreased on day 1 and returned to the baseline by day 28 in the control individuals. Our study reveals sex-dimorphic effects of previous mild COVID-19 and suggests that viral infections in humans can establish new immunological set-points that affect future immune responses in an antigen-agnostic manner.

Surfaceome Profiling of Cell Lines and Patient-Derived Xenografts Confirm FGFR4, NCAM1, CD276, and Highlight AGRL2, JAM3, and L1CAM as Surface Targets for Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in children. The prognosis for patients with high-grade and metastatic disease is still very poor, and survivors are burdened with long-lasting side effects. Therefore, more effective and less toxic therapies are needed. Surface proteins are ideal targets for antibody-based therapies, like bispecific antibodies, antibody-drug conjugates, or chimeric antigen receptor (CAR) T-cells. Specific surface targets for RMS are scarce. Here, we performed a surfaceome profiling based on differential centrifugation enrichment of surface/membrane proteins and detection by LC-MS on six fusion-positive (FP) RMS cell lines, five fusion-negative (FN) RMS cell lines, and three RMS patient-derived xenografts (PDXs). A total of 699 proteins were detected in the three RMS groups. Ranking based on expression levels and comparison to expression in normal MRC-5 fibroblasts and myoblasts, followed by statistical analysis, highlighted known RMS targets such as FGFR4, NCAM1, and CD276/B7-H3, and revealed AGRL2, JAM3, MEGF10, GPC4, CADM2, as potential targets for immunotherapies of RMS. L1CAM expression was investigated in RMS tissues, and strong L1CAM expression was observed in more than 80% of alveolar RMS tumors, making it a practicable target for antibody-based therapies of alveolar RMS.

Species interactions, stability, and resilience of the gut microbiota - Helminth assemblage in horses

The nature and strength of interactions entertained among helminths and their host gut microbiota remain largely unexplored. Using 40 naturally infected Welsh ponies, we tracked the gut microbiota-cyathostomin temporal dynamics and stability before and following anthelmintic treatment and the associated host blood transcriptomic response. High shedders harbored 14 species of cyathostomins, dominated by Cylicocyclus nassatus. They exhibited a highly diverse and temporal dynamic gut microbiota, with butyrate-producing Clostridia likely driving the ecosystem steadiness and host tolerance toward cyathostomins infection. However, anthelmintic administration sharply bent the microbial community. It disrupted the ecosystem stability and the time-dependent network of interactions, affecting longer term microbial resilience. These observations highlight how anthelmintic treatments alter the triangular relationship of parasite, host, and gut microbiota and open new perspectives for adding nutritional intervention to current parasite management strategies.

Assessment of the Molecular Responses of an Ancient Angiosperm against Atypical Insect Oviposition: The Case of Hass Avocados and the Tephritid Fly Anastrepha ludens

Anastrepha spp. (Diptera: Tephritidae) infestations cause significant economic losses in commercial fruit production worldwide. However, some plants quickly counteract the insertion of eggs by females by generating neoplasia and hindering eclosion, as is the case for Persea americana Mill., cv. Hass (Hass avocados). We followed a combined transcriptomics/metabolomics approach to identify the molecular mechanisms triggered by Hass avocados to detect and react to the oviposition of the pestiferous Anastrepha ludens (Loew). We evaluated two conditions: fruit damaged using a sterile pin (pin) and fruit oviposited by A. ludens females (ovi). We evaluated both of the conditions in a time course experiment covering five sampling points: without treatment (day 0), 20 min after the treatment (day 1), and days 3, 6, and 9 after the treatment. We identified 288 differentially expressed genes related to the treatments. Oviposition (and possibly bacteria on the eggs' surface) induces a plant hypersensitive response (HR), triggering a chitin receptor, producing an oxidative burst, and synthesizing phytoalexins. We also observed a process of cell wall modification and polyphenols biosynthesis, which could lead to polymerization in the neoplastic tissue surrounding the eggs.